应力-应变曲线 — CAE术语解释
应力-应变曲线
在FEM中输入材料数据时,要求「应力-应变曲线」表,这就是拉伸试验的结果本身吗?
应力-应变曲线的理论基础
应力-应变曲线的基础
应力-应变曲线是拉伸试验结果的图表,对吧?但在CAE中设置材料模型时,我们使用的是这条曲线的「哪一部分」?是全部吗?
很好的问题。在最基本的线性弹性分析中,CAE只使用曲线的「初始斜率」,即弹性区的很小一部分。这个斜率就是杨氏模量E。例如,结构钢S45C根据JIS G 4051规定,E约为205 GPa。通过测量应力-应变曲线初始直线部分的斜率来确定这个值。在塑性分析中,屈服点之后的整条曲线就变得重要起来了。
屈服点就是曲线上明显弯折的地方,对吧?但我听说铝合金A5056这样的材料没有明显的弯折。CAE如何处理这种情况?
完全正确。许多有色金属和加工硬化钢没有明确的屈服点。在这种情况下,我们使用「0.2%耐力」这样的标准化定义。也就是说,将对应于塑性应变0.2%的应力视为屈服应力。在CAE软件的材料输入中,我们在Yield Stress字段输入这个0.2%耐力值,比如A5056约130 MPa,然后单独定义其后的硬化规则。
硬化规则?是用数学公式来表示屈服后曲线形状的东西吗?具体有哪些公式?
是的。最常见的是「幂次硬化规则」。它用以下公式近似塑性区的应力与塑性应变关系:
在拉伸试验中,试样会出现缩颈然后断裂,这个「缩颈」开始的点在应力-应变曲线上能看出来吗?CAE能预测断裂吗?
曲线上对应的是「最大应力点」。在那之后,公称应力下降,但由于缩颈部分的截面积急剧减小,「真应力」实际上继续上升。要在CAE中预测断裂,需要使用真应力-真应变曲线,还需要「破裂应变」等破坏准则。例如,Abaqus的Ductile Damage模型基于塑性应变和应力状态(应力三轴度)来计算损伤的发展。
应力-应变曲线的数值计算方法
FEM中的材料非线性处理
线性分析的支配方程是
表面形式相同,但实质完全不同。当材料变成非线性时,刚度矩阵
「逐步增加」是指把应力-应变曲线近似成阶梯形,对吧?这个步长如何确定?
正是那样。软件用多条直线段来近似曲线。步长(增量大小)通常是自动控制的。例如,Ansys Mechanical中,如果塑性应变的增量过大,就无法收敛,求解器会自动减小步长并重新尝试。初始步长通常设置为无量纲时间的0.01~0.1。
在CAE软件中输入应力-应变数据时,经常可以选择「公称」和「真」,应该输入哪一个?转换是自动进行的吗?
塑性分析的本质是计算不伴随体积变化的塑性应变,所以原则上应输入「真应力-真应变」数据。如果输入公称值,在大变形区会出现较大误差。许多软件(Abaqus、LS-DYNA等)有选项可以在输入公称值时自动转换为真值,但理解转换公式
数据点越多精度越高吗?10个点和100个点的结果会差很多吗?
在屈服点附近和曲率剧烈变化的区域需要细致的点。但是,在单调硬化区输入过多点数可能导致数值问题,例如相邻点间斜率反转。实际上,15到20个点通常足以捕捉曲线特征。软件会对这些点进行插值使用。
应力-应变曲线的实际应用
分析工作流程与验证
在实际项目中,如何获得材料数据?每次都做拉伸试验吗?
如果不是新材料,首先查找材料制造商的数据表。例如,ThyssenKrupp钢板或Kobe Steel铝合金,技术资料中会记载0.2%耐力、抗拉强度、伸长率的最低保证值。CAE需要整条曲线,所以要查阅公开数据库(如JAHM、MatWeb)或过去的内部试验数据。如果都没有,只能外包试验,成本约5-10万日元/样本。
获得的数据在输入CAE软件前,需要检查哪些项?
至少检查三点。首先,初始斜率(杨氏模量)是否对材料类型合理(钢约200 GPa)。其次,屈服后的曲线是否单调递增(没有应力下降的点)。第三,最大应力(抗拉强度)和断裂点的应变是否与数据表的公称值偏差不大。特别要注意是否有异常大的数据,如真应变超过1.0。
设置好材料模型后,有简便的方法来验证设置是否正确吗?
最有效的是「单轴拉伸分析」基准测试。在CAE软件中创建与试样相同的形状(例如JIS 13B试样),应用输入的材料模型进行拉伸分析。将输出的反力-位移曲线与从实验数据计算的公称应力-公称应变曲线比较。这样可以定量确认屈服点、加工硬化趋势、抗拉强度是否被重现。
复杂产品的分析中,各部件使用不同的材料。需要为所有材料准备非线性数据吗?
这是现实的难题。通常只对「主要能量吸收部件」或「大变形部位」的材料应用详细的非线性数据。例如,汽车碰撞分析中,车身骨架(高强度钢)必须有非线性数据,但刚性很高、几乎不变形的支架(铸铁)等往往用线性弹性(仅杨氏模量和泊松比)就足够了。平衡分析目的和成本来决定材料模型的详细程度,这是工程师的技巧所在。
应力-应变曲线的软件比较
主要软件的材料模型输入
Ansys Mechanical和Abaqus/CAE在应力-应变曲线的输入方法上有区别吗?
基本思想相同,但界面和术语不同。Ansys在Engineering Data源中选择「Multilinear Isotropic Hardening (MISO)」或「Bilinear Isotropic Hardening (BISO)」等硬化模型,在表格中输入应力-塑性应变的对应值。而Abaqus在材料属性中添加「Plastic」,直接定义包含「Yield Stress」和「Plastic Strain」列的表格。需要注意的是,Abaqus假设第一个点是屈服点(塑性应变0)。
我听说COMSOL Multiphysics有「充实的材料库」,那里的曲线可以直接用吗?
COMSOL的内置材料库中,线性弹性参数(E、ν)可靠,但非线性应力-应变曲线数据有限。而且,那只是某个特定牌号、热处理状态的「一个例子」。例如,「Structural Steel」的非线性数据可能来自SS400,也可能来自SM490,不清楚。实际设计应该按照自己的热处理条件输入数据,库数据只能作为学习和初期评估用。
LS-DYNA这样的显式求解器,输入方式不同吗?特别是考虑应变率效应时。
差异很大。LS-DYNA用数字选择材料模型(MAT_),例如弹塑性用「MAT_024 (Piecewise Linear Plasticity)」。在碰撞分析中,「应变率依存性」不能忽视。MAT_024可输入多条不同应变率(如0.001/s、1/s、100/s)的应力-应变曲线。钢材中,应变率增加10倍时,屈服应力约上升5-10%,这种效应可以通过输入表现。
自由CAE软件(CalculiX、Code_Aster等)中,高级材料模型的输入很困难吗?
功能不亚于商用软件,但用户界面门槛很高。CalculiX和Code_Aster需要直接编辑输入文件(.inp、.comm)。例如,在Code_Aster中定义多直线硬化,需在TRACTION关键字下用规定格式写出应变和应力列表。没有GUI,很难看到数据和发现错误。实际做法是在商用软件中确认设置后再移植。
应力-应变曲线的故障处理
常见收敛问题与对策
在包含非线性材料的分析中,荷载步骤中途停滞,显示「不收敛」错误。最先应该怀疑什么?
首先怀疑「材料模型输入错误」。特别是,检查应力-应变数据的起点(屈服点)的应力值是否过低,或者屈服后的曲线中应力是否出现下降(「软化」)。许多求解器在屈服后应力低于屈服应力时会变得物理上不稳定,导致无法收敛。把数据绘成图表,目视确认单调递增,这是第一步。
材料数据单调递增的情况下仍不收敛,下一步检查什么?
下一个是「网格」和「单元类型」。在大塑性变形区,单元形状极度扭曲会导致计算破裂。特别是一次单元(四面体、六面体的线性单元)容易发生体积锁定。塑性分析推荐使用二次单元,或者一次单元也可选用减积分单元(如Abaqus的C3D8R)。此外,应变集中区的网格细化也有效。
「体积锁定」具体是什么现象,对结果有什么影响?
当材料呈近似不可压缩性的塑性变形(泊松比ν=0.5)时,完全积分的一次单元会过度限制体积变化。结果单元刚度异常高,变形受阻。这就是体积锁定。后果是部件计算刚度过高,反力过大,或收敛性恶化。对策包括使用前述的减积分单元或Abaqus的「混合单元」(如C3D8H)。
用幂次硬化式定义材料数据和用多点表格输入,收敛性有差别吗?
有可能有差别。用硬化式定义时,应力对应变的导数(切线刚度)连续计算,牛顿法的收敛性较好。而多点表格在点间线性插值,切线刚度在点前后不连续变化,这会妨碍收敛,特别是迭代中跨越不连续点时。Ansys的「MISO」模型提供「平滑化」选项,可以平滑插值切线刚度,这有所帮助。
相关主题
价值
详细
错误