材料的名义(Nominal)应力、应变是基于变形前的数据计算得到, 其中 为试件初始截面面积, 为试件初始长度。名义应力、应变也叫工程(Engineering)应力、应变。 CAE软件需要采用基于变形后的应力、应变,即真实的应力、应变。 其中 为试件当前截面面积, 为试件当前长度。 两种应力、应变的转化公式为: 下面来推导这两个公式。 一) 了解定积分的精确定义。 font.sans-serif'] = ['SimHei'] # 正常显示中文标签 plt.rcParams['axes.unicode_minus'] = False # 正常显示负号 # 名义应力应变 562.60897,571.00612,579.40327,587.800419,596.197567,604.594716,612.99186, 621.3890,629.78616,638.18331]) n = len(sigma_N) print(n) #真实应力应变 np.hstack( (v1, eps_T) ) sigma_T1 = np.hstack( (v1, sigma_T) ) fig, axs = plt.subplots(1, 1, figsize=(14,6)
应变是一个无量纲的量,通常以小数或百分数表示。二、Abaqus应力应变分析的流程基础概念2.1模型建立在Abaqus中进行应力应变分析,首先要建立准确的几何模型。 三、Abaqus应力应变分析中的关键概念3.1应力分量与应力张量在三维空间中,应力状态可以用应力张量来描述。应力张量包含9个分量,其中3个正应力分量和6个剪应力分量。 主应变是指材料在某一点处沿主应力方向的线应变;等效应变则用于描述材料在复杂应力状态下的总体变形程度,在塑性变形分析中具有重要意义。 四、Abaqus应力应变分析的应用场景4.1结构强度评估通过应力应变分析,可以评估结构在设计载荷下的应力分布情况,判断结构是否满足强度要求。 例如,通过模拟材料的拉伸、压缩、弯曲等试验,获取材料的应力应变曲线,验证材料模型的准确性,为材料的选用和性能改进提供依据。掌握Abaqus应力应变分析的核心概念,是开展高效、准确有限元分析的基础。
拉伸试验得到的应力应变,通常是指工程应力和工程应变,用于计算应力应变的横截面积和长度,是未变形的初始横截面积和初始长度(便于测量)。 与之对应的,还有真应力和真应变,用于计算应力应变的横截面积和长度,是变形后的横截面积和长度。 在应力低于比例极限的情况下,应力σ与应变ε成正比,即σ=Εε;式中E为常数,称为弹性模量或杨氏模量,是正应力与正应变的比值,弹性模量的单位与应力的单位相同。 剪切模量的定义与之类似,是切应力与切应变的比值。金属的应力应变曲线,通常分为四个阶段:弹性阶段、屈服阶段、应变硬化阶段和颈缩断裂阶段。 所以,对于有限元分析结果,仅仅根据应力是否超过抗拉强度来评价结构是否发生断裂,是不太准确的。脆性材料受到压缩时,会被压碎,但压缩强度极限要比拉伸强度极限大得多。
对于有限元分析(FEA)而言,应力应变计算是评估结构强度与安全性的核心工作。 孔边、拐角等应力集中部位,应力在极小范围骤变,即便加密网格,单元应力仍因 “尺寸效应” 失真,靠近集中点的应力远低于真实值(即 “应力锁死”)。 外推(如将内部积分点应力外推至表面)可得到真实表面应力,与许用应力对比完成校核。另一方面,规范要求的 “等效应力” 需外推整合。 单元应力应变外推不是 “优化步骤”,而是解决精度偏差、突破单元限制、满足设计标准的 “必要操作”,核心是通过数学方法修正误差,还原真实应力分布,让有限元结果从 “数值解” 变为 “工程可用解”。 实际应用中,工程师需结合结构类型(梁、板、实体)与分析目标(应力集中、疲劳寿命),选择线性、二次或子模型外推法,同时通过 “网格收敛性分析” 验证结果可靠性,为结构设计提供安全依据。
限于时间关系,这里只列出基本的应力提取程序,写入文本比较简单,随意百度即可实现,不再给出。 fieldOutputs['S'] #需要输出的场变量 sub_f_S1 = f_S1.getSubset(region=region) S_Value = sub_f_S1.values #读取所有单元的应力值
系统概述 该系统量测智能手表PCB板应力和应变,使用4组应变片,先将黏贴于PCB板上的应变片与PCLD-8813连接,透过此讯号调理板先将噪声滤除,再接到PCIE-1813应变量测卡,并依照应变片的规格设定 PCIE-1813的参数、取样频率、放大率与激励电压之后,即可开始采集应变片上的讯号应力应变。 应变又称“相对变形”,是物体由于外因(载荷、温度变化等)使它的几何形状和尺寸发生相对改变的物理量。 应力应变量测是利用应变片(Strain Gauge) 中金属导线阻值之变化来量测应变量。 ,为使应变量分析结果更为精确,数据采集卡的各通道之间必须能够同步取样,以避免不同通道间数据因取样时间延迟而使分析的结果产生偏差。 该系统通过能提供高精准应力量测的数据采集产品PCIE-1813,实现板高分辨率、高性价比量测系统,客户还可以依照不同的应变测试需求,搭配可弹性选择更多插槽数工业计算机,通过多卡同步功能快速建置更多通道的高精度应力应变量测系统
结构在线监测:工程师通过安装在飞机机翼的应变片,监测飞行器的结构变形等物理参数,进而获取飞行器的运行状态;2.热力耦合实验:对于高超声速飞行器,部分实验的迭代成本很高,通过物联网加数值仿真,给设计人员提供理论支撑 近来,结合实际项目对结构应力应变监测方案相关的技术进行汇总,具体如下所示: 图1 飞行器结构应力应变监测系统架构 试验方案 系统架构 图1展示了飞行器结构监测的系统架构,主要包含:1.传感器选型及安装 :购买市面上现有的应变片,后续采用胶水进行粘贴;2.数据采集系统:上图中采用了两种数据采集方案,其中,方案一(图1-f)内置有电源和4G通信模组,能够将传感器采集的数据直接上传到E-twin物联网平台; 实现方案 我们先展示方案二的具体实施方案,通过边缘计算将传感器采集的应变数据转化为有限元模型的边界条件。 f"发生其他错误: {e}") finally: # 断开连接 client.close() 上述代码能够远程获取传感器的实施结果,后续转化为力学工程师熟悉的数据形式,方便对问题进行分析定位
有限元计算结果主要包括位移和应力两个方面。在位移方面一般无需进行整理工作。但是由于应力在单元间不连续,一般要做一定的修正。比如,三角形常应变单元也是常应力单元。 这个CST单元建立的有限元模型得到的X方向的正应力(sigma_xx)分布如图所示: 为了由计算结果得到弹性体内某一结点的接近实际的应力值,通常可采用绕结点平均法或两单元平均法。 所谓绕结点平均法,就是把环绕某一结点的各单元常应力加以平均,用以表示该结点的应力。为了使由这样的平均得到的应力能够较好地表示结点处的实际应力,环绕该结点的各个单元的面积不应相差太大。 以下是绕结点平均法的代码: 计算结果如下: 得到结点应力后,用MATLAB画sigma_xx云图
(1)在Abaqus中定义塑性 利用Abaqus进行有限元模拟的过程中,我们有时候需要考虑材料的非线性,当定义材料塑性应力应变关系时,我们必须采用真实应力和真实应变[1]。 应力应变曲线 实验数据通常是工程应力和工程应变,因此,我们通常需要把工程应力(曲线A)转换成真实应力(曲线B)进行有限元模拟。 采用等效应变和等效应力的意义在于,等效应力将6个应力分量对变形的作用等效于一个单向拉伸力的作用,等效应变将6个应变分量,等效于一个单向拉伸力所产生的应变。 1、材料损伤起止判据 其中,等效塑性应变,失效初始应变 表1 Ti-6Al-4V的Johnson-Cook剪切失效参数 -0.09 0.25 -0.5 0.0014 3.87 2、材料失效判断 采用能量密度法作为完全失效的准则(材料断裂失效能是材料整个失效演化阶段应力应变的积分)。 [1]庄茁等. 基于ABAQUS的有限元分析和应用[M]. 清华大学出版社, 2009.
2、基于应变疲劳分析算法 稳定循环应力-应变迟滞曲线如下图,一般用Ramberg-Osgood方程表示, (1) 其中, 为弹性模量, 为循环硬化系数, 为循环应变硬化指数 图1 稳定的应力- 图5 ABAQUS材料设置 表1 材料属性 3.2 ABAQUS 边界条件 疲劳式样如图6,标注为英寸,需转为mm制。 图6 疲劳式样(单位为英寸) 图7 FEA边界条件 分别对不用应变幅进行有限元仿真,图8为应变幅1.72%的总应变分布,标距范围内均匀分布。 图8 应变幅1.72%的总应变分布 图9 中间位置单元的应力-应变响应 4、 FE-SAFE疲劳分析 FE-SAFE 中有两种方法来做弹塑性疲劳分析。 2) FEA中进行弹塑性分析,将应力应变导入到FE-SAFE中进行疲劳分析。 FE-SAFE疲劳分析分为以下几步,导入有限元模型、选择材料、设置分析组(Group)、设置载荷。
通过该模型,深入分析了在单向拉伸条件下,轮胎橡胶中裂纹的扩展方向。为了验证模拟结果的准确性,将所得计算结果与实际试验数据进行对比,发现两者高度一致,充分证明了本论文所建立的模型和分析方法的可靠性。 3 有限元计算结果及分析第一次加载使裂纹尖端单元应力小于橡胶材料单轴拉伸试验时对应的破坏应力;第二次加载使裂尖单元应力超过橡胶单向拉伸试验破坏时对应的应力大小。 本次试验单边裂纹试样预置裂纹长度为5mm,中心裂纹试样为6mm。裂纹应平滑,垂至于载荷轴线,前沿应锐利齐平,否则将影响实验结果。 通过对实验数据处理得出了单向拉伸的应力应变图7。图 7-9分别是对应中心裂纹,单边裂纹,中心斜裂纹的应力应变曲线。从图中可以看出单边裂纹试样断裂的应力应变最小,斜裂纹最大。 比较单轴拉伸应力应变图与带裂纹试件的应力应变图可以发现,虽然带有裂纹的式样尺寸比单轴拉伸试样的宽度大很多,但是从应力应变图上可以看出带有裂纹的拉断应力远小于没有裂纹的拉断应力。
不同的材料在不同载荷形式的作用下具有不同的疲劳破坏形式,因此,临界平面法综合参数的定义方式具有多种:(1)应力作用下的疲劳裂纹萌生;(2)应变作用下的疲劳裂纹萌生;(3)应力和应变的共同作用使裂纹萌生等 临界平面法根据单元积分点的应力和应变参量值决定临界平面的位置,具有一定的物理意义。 求解疲劳载荷一个循环周期内接触区域所有单元每个积分点的应力和应变; (2)采用Python提取ABAQUS的计算结果,把接触区域每个单元中每个积分点的应力应变输出到文件里边。 研究表明,应力和应变采用积分点比节点更加准确; (3)采用MATLAB,计算每个单元每个截面的应力应变值; (4)得到每个单元每个截面上的SWT参数,取最大值作为单元的SWT参数; (5)得到接触区域 图2 SWT参量估计微动疲劳裂纹萌生及寿命有限元模型示意图 压头以及试件的材料属性:在分析过程中,压头与试件都选用采用的弹塑性材料Ti-6Al-4V,弹性模量E1=116GPa,泊松比v1=0.342,
下面就其力学原理做大致的分析。 这里的应变是相同的,但由于材料的弹性模量突然发生变化,应力也会发生变化。 为了求这种复杂的应力分布,将复合梁转换为由单一材料制成梁。 如果将较"软"的材料2转变为较"硬"的材料1,分析过程同上。 一旦复合梁被转换成具有单一材料的梁,转换横截面上的正应力分布为连续线性,如图2g或图2h所示。 \sigma_{B^{'}} = \frac {2*10^6 *(170-36.38)}{9.358*10^{12}} = 28.6MPa C的正应力为 \sigma_{C} = \frac {2 *10^6 *36.38}{9.358*10^{12}} = 7.78MPa B点的正应力为 \sigma_{B} = n\sigma_{B^{'}} = 1.71MPa 实际正应力分布如图3d所示
特别是,大变形和大应变几何非线性经常与塑性联系在一起。如果你预测可能会发生大变形,你必须激活这些效应。NLGEOM,ON. 对于大应变分析,材料应力-应变属性应当使用应力+对数化的形变的方式输入。 不建议使用该选择于循环或高度不成比例载荷历史的小形变分析。建议用于大变形分析(large strain)。MISO选项可以包含多达20个不同温度曲线,每个曲线允许100个不同应力应变点。 你可以通过拟合材料拉伸应力-应变曲线来得到这四个常数。不同于MISO,不需要担心如何恰当选定应力-应变点来输入。但是该选项只是适用于如下图所示的拉伸曲线。该选项适合大应变分析。 4 非线性瞬态热应力分析中的重要命令 ①输出控制(结果输出到数据库),建立存储规格。 间接法计算热应力时,热分析的结果文件要作为结构分析的热载荷输入。 例1: NSUBST,6 OUTRES,ERASE 设置到默认值,对于静态和瞬态分析,默认的是输出每一载荷步的最后子步的所有结果;谐态分析是每一子步。
输出的单元应力应变曲线为什么和输入的不一样? 单元应力为什么比屈服强度还高? 单元应力超过定义的最大屈服应力后的发展趋势是怎样的?为什么会出现应力增大情况? 混凝土输入的是真实应力应变曲线还是名义应力应变曲线? 为了解决上面的这些问题,这里准备了一个非常简单的模型,一个平面的四边形单元,A点进行完全固定,B点约束Y方向位移,CD点施加-Y方向位移载荷,对单元形成单轴压缩分析;右侧的模型是对单元四个面额外施加P的围压 结论: 当模型选择平面应力单元时,得到的应力-应变曲线与输入的材料一致; 当模型选择平面应变单元时,应力-应变曲线各处均大于输入的材料性质,这是由于平面应变单元存在平面外法向约束,因此并不是纯粹的单轴压缩模型 ; 输入定义的末尾应变为0.37,后续的应力-应变曲线是水平延伸的,单轴条件下,不存在增大情况。
1)传感器系统 传感器系统的选择应与监测项目相匹配,一般包括: (1)结构静、动力反应监测传感器:用于监测桥梁在温度、不均匀沉降及荷载等作用下的应力、应变等反应,解决可靠性评估中的效应监测问题,同时也为自诊断系统提供模式样本及结构反演分析的原始数据 外场数据采集计算机系统对被测物理量量测结果进行预处理(如量测结果的修正换算,主应变计算等)并形成数据文件,通过通讯光缆传输至管理指挥中心。 3)结构分析和评估系统 主要包括结构分析和评估软件系统等。 1)桥梁应力监测测点布设 应力是运营桥梁局部响应最主要的表现量,对运营期活载作用下桥梁受力最不利截面、桥梁病害发展较快的部位等进行应力监测,了解桥梁关键截面应力情况,选取应变计作为应力监测仪器。 应变和应力分析:剔除徐变的应变和应力分析,应变和应力的分布范围统计分析,应变和应力的趋势变化分析和极值分析,应变和应力与温湿度之间的关系分析,应变和应力与车辆荷载之间的相关性分析,应变和应力与风荷载之间的相关性分析 模态分析:根据加速度传感器和应变传感器的测试结果,采用多种方法,进行模态分析,获得桥梁的频率、振型和阻尼比。
ABAQUSUMAT 初步 ABAQUS中虽然提供了很多的单元和材料的本构模型,但是在有限元模拟分析过程中实际情况是比较复杂的,因此,为了更好地跟实际情况相吻合,提高有限元模拟结果精度,我们需要对ABAQUS UMAT子程序的核心内容是定义材料在积分点处的雅克比矩阵,定义应力增量和应变增量的关系。 对于线性分析而言,我们可以直接采用一个物理方程中的矩阵D定义应力和应变之间的关系,但是在有限元计算中,我们遇到的很多情况都是非线性分析,因此,采用定义雅克比矩阵的方式,定义非线性过程中应力应变之间的关系 首先,弹性力学基本方程的矩阵形式为: 应力向量、应变向量的表述 例如:线弹性UMAT初步: SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD STRESS C - ------ - C PROPS(1) - E C PROPS(2) - NU CC C IF (NDI.NE.3) THEN C WRITE (6,
ANSYS Workbench有限元分析的基本流程即可。 1.5.1案例介绍 某如图1-24所示不锈钢钢板尺寸为320mmX50mmX20mm,其中一端为固定,另一端为自由状态,同时在一面上分布有均布载荷q=0.2MPa,请用ANSYS Workbench求解出应力与应变的分布云图 图1-26 创建分析项目 提示:本例是线性静态结构分析,创建项目时可直接创建项目B,而不创建项目A,几何体的导入可在项目B中的B3栏Geometry中导入创建。 本例的创建方法在对同一模型进行不同的分析时会经常用到。 对话框中选择文件路径,导入char01-01几何体文件,如图1-28所示,此时A2栏Geometry后的 变为 ,表示实体模型…… 文章来源:技术邻-大龙猫 全文链接:Workbench实例入门-悬臂梁的应力变形仿真分析
但是,两种材料性能相差巨大,比如钢筋的弹性模量比混凝土高出一个数量级,而且其应力-应变关系在受拉区和受压区对称分布,混凝土则不同。 考虑受拉刚化效应的方法总体上有三种: 根据粘结应力-滑移本构模型,建立粘结单元; 增大钢筋的刚度; 基于混凝土的平均应力和平均应变关系建立的,开裂后依然有一定的抗拉强度的模型; 如果考虑到裂缝间距 根据开裂后应力-应变关系形式的不同,有多种模型,如下图所示。 若分析素混凝土结构,(*tension stiffness)即素混凝土的单轴受拉应力-应变曲线,可直接按照规范附录C的建议取用;若分析钢筋混凝土结构,(*tension stiffness)为考虑了纵筋有利作用的混凝土受拉应力 ABAQUS帮助文档建议,对于配筋率相对较大的结构,当采用较细密的网格时,可以假设混凝土开裂后,应力从峰值点线性衰减为零,对应的应变为开裂应变的10倍。
固体力学中有三类变量:应力、应变和位移。 这三类变量通常有以下三种表示方法: 工程表示 正交张量表示 数学(矩阵)表示 在弹性范围内,这三种表示方法的等同的。 (1) 应力 一点的应力状态用6个独立的分量表示。 (直角坐标系) (2) 应变 一点的应变状态也用6个独立的分量表示。 (直角坐标) 笛卡尔坐标 剪应变的工程表示比张量表示差1/2 (3) 位移 一点的位移用3个独立的分量表示。 三维弹性理论问题的未知量有6个应力分量,6个应变分量以及3个位移分量。一共15个未知量。 实际上,应力、应变、位移都是弹性体内各点坐标的函数,即都是场变量。 数组是计算机语言的一种数据结构。在编程时,张量都要由数组来存储。比如,四阶张量通常由二维数组表示,二阶张量由一维数组表示。 应力张量 在程序中表示为 对于平面问题 在程序中表示为 应变张量 在程序中表示为 注意剪应变前面加系数2,意思是工程剪应变等于2倍的张量剪应变。更方便矩阵运算。