ANSYS弹塑性本构:切线模量填多少?以及6个让非线性算不收敛的坑
割线斜率取屈服强度的 8~10 倍——这一条经验值,够你少走半年弯路。
非线性算不收敛,是结构仿真最劝退的一关。但弹塑性这一类,八成的不收敛不是模型错了,而是三个开关没打开、一个参数填错了。
这篇整理自 ANSYS 结构力学实战零基础到就业课程《弹塑性本构》一讲:材料曲线怎么输、切线模量取多少、2D 与带梁 3D 两个完整案例,以及最后那节最值钱的非线性收敛调参清单。
一、非线性有三大类,本节讲第二类
ANSYS 结构非线性的三大类
① 几何非线性讲杆单元时已经涉及过(大变形 / 刚性硬化)
② 材料本构非线性← 本节内容:弹塑性
③ 接触非线性
二、弹塑性到底是什么
弹塑性是材料力学与固体力学中的核心概念,指的是:物体在外力作用下同时产生可恢复的弹性变形和不可恢复的塑性变形的性质。
当材料的受力超过弹性极限或屈服强度时,应力和应变呈非线性关系,卸载后会出现残余变形。
单轴拉伸曲线的三段行为
应力 < 比例极限 → 材料呈线性
应力 < 屈服点 → 材料呈弹性,载荷卸载后所有应变都可以恢复
应力 > 屈服点 → 材料呈塑性,载荷卸载后应变不能恢复
两个工程上的实用认知:
比例极限和屈服点非常接近,所以有些软件直接假设两者相等——你在不同软件里看到的"屈服强度"含义可能不完全一致。
结构的塑性响应是基于单轴实验获得的。通过单轴应力应变实验可以得到材料的比例极限、屈服极限(弹性极限)、应力强化。输入数据是塑性应变 + 真实应力,适用于大应变和循环载荷(这一点和橡胶材料的输入方式是一样的)。
图1 单轴拉伸的弹塑性应力-应变曲线
⚠ 一个版本坑:曲线的下降段输不进去
在目前版本的 ANSYS 里,输入应力应变曲线时,它的下降段是没办法输入的。
而在19 版本之前是可以输入的。
→ 如果你的材料曲线有明显软化段(下降段),要么用旧版本,要么想别的等效办法。这一点在做失效、断裂类分析时尤其要注意。
三、常用本构与那个关键的经验值
常用的弹塑性本构有:双线性本构、随动本构(随动强化)、多线性强化本构,另外还有各向同性本构。
而真正让人卡住的,往往是切线模量(割线斜率)该填多少。课程里给的经验值非常明确:
割线斜率(切线模量)经验值
普通结构钢 800 ~ 8000 MPa
微小塑性应变(< 0.5%)非常接近时 可设到 100 MPa
★ 通用规则:割线斜率 = 屈服强度的 8 ~ 10 倍
通常默认就按 8~10 倍来设
这条规则在后面案例里得到了印证:课程给的弹塑性参数是屈服强度 250 MPa、切线模量 2000 MPa——正好是 8 倍。
图2 常用弹塑性本构:双线性、随动强化、多线性强化
四、案例 1:2D 平面应力拉伸(含循环加载)
建模步骤
1 插入静力分析。弹性模量和泊松比是必须输入的(这一点和橡胶材料不同);单位改成 MPa
2 材料数据(网上下载的 Excel 表)→ 用曲线拟合得到本构
3 结构设为 2D 模型,单位选毫米;必须在 XY 平面,并且关于 Y 轴对称(关于 X 轴对称不做要求)
4 由草绘生成表面(100 × 20 的概念面)
5 确认是平面应力,厚度 5 mm,材料改成自定义材料
6 网格:选择线 → 分区给 25 个;尺寸调整选这两条线 → 分区改成 5 mm / 5 个
7 载荷给 10 步;打开弱弹簧;打开大变形
8 位移来自指定线、方向 X;数据表编好后统一改成 0 再 Ctrl+V 粘贴
9 一个面全约束;另一个面插入位移:X 自由、Y 不变
10 分析设置:10 步改成 100 步;时间步用自动(程序控制)
图3 坑:更新后发现网格没有映射(变成了扫掠形式),需要手动插入映射网格
★★★ 坑:累积塑性应变看不了
插入变形、应变、塑性应变、等效塑性应变之后,会发现累积塑性应变看不了。
解决办法:到「输出(Output)」里,把「非线性」打开。然后重新求解就能看到了。
这是弹塑性分析里最常见的"结果项灰着点不动"的原因。
图4 在「输出」里打开「非线性」,累积塑性应变才可用
结果查看要点
结果摘要里的可用求解项 往下拉能看到 EL = 弹性应变、EPPL = 塑性应变;可单独看 X 方向
应力-应变图表 Ctrl 同时选中应变和应力 → 在结果里插入图表,即可得到数据表
时间范围 时间结束是 1800 秒,若要全部显示,坐标至少给到 3600
趋势 等效应变随时间减小,累积塑性一直在增加
★ 反直觉现象 改成 400 步反而出现奇怪结果;有时候步数少、让自适应去调,反而能算出来
加硬化 vs 不加硬化:课程里做了对比——把结果复制过来(不来自上游),再插入一个弹塑性本构双线性并给参数,重新求解。加了硬化参数之后,塑性区的分布和应力水平都明显不同。硬化参数不是可选项,它直接决定结果。
五、案例 2:带梁的常规模型(对标 Abaqus)
★★★ 第一个坑:导入模型时线单元默认不确认
双击打开导入的模型,你会发现中间是没有任何东西的,而且这时候想去改类型也改不了,只能删掉重来。
原因:导入模型时,梁单元/线单元一般情况下默认是不确认的。
→ 导入带梁的模型时,一定要在这里把线单元勾选上。
建模步骤
1 单位选毫米;创建概念横截面:圆管,直径 2.5、壁厚 1.25 mm
2 形成共享节点
3 材料两套:3D 打印材料(E = 2.69E10、泊松比 0.32)与结构钢
4 插入弹塑性本构:屈服 250、切线模量 2000(= 8 倍)
5 本构材料从上游复制(copy)过来;类型选多线性
6 接触:默认没有接触 → 插入手动连接,先选上表面,隐藏体后框选线,摩擦系数 0.3;打开大变形,调整接触、必要时翻转;下表面同样操作
7 网格:默认方法、自适应,选 4
8远程位移:选面,Z 方向走 5 mm;另外两个方向选自由(因为会有来回晃动);下表面插入远程位移并全部固定
9 分析设置 50 步;热弹簧(弱弹簧)打开;大变形打开
10 计算慢时:把两个零件改成刚体、网格重新画,然后重新提交计算
图5 导入带梁模型:不勾选线单元,打开后中间空空如也且无法改类型
★ 第二个坑:新版本对输入值做了限制
课程里尝试把某些参数改成 0.22、−100、100 时,输入框会标黄(不允许输入)。
之前的版本没有这个问题,是可以输入的;现在不行了。遇到标黄就说明超出了当前版本允许的取值范围,要换一种表达方式。
结果对标
反力11kN,与论文中Abaqus的结果基本一致
梁的应力需要单独打开显示。
图6 结果对比:与论文中用 Abaqus 算的结果基本一致
六、★ 非线性收敛调参清单(最值钱的一节)
课程最后的总结很实在:对于 ANSYS 来说,弹塑性 + 几何非线性需要调的参数比较多。下面这几项是不收敛时该逐个检查的:
调什么 → 怎么调
非线性控制 → 一般选择非对称(Unsymmetric)
大变形 → 一般打开
线搜索(Line Search) → 也可以打开,改到 0.1
位移收敛 → 可以设置,例如改成 2%、或改成 100 N——这些都会提高收敛性
弱弹簧 / 热弹簧 → 模型约束不足时打开
时间增量 → 越小越容易收敛
收敛过程中,在「求解信息」里盯这三项
· 最大残差 比较大时说明还远没收敛
· 线性搜索 是否被激活、迭代了几次
· 时间增量 是否被不断二分(时间增量越小越容易收敛)
图7 非线性控制:一般选择「非对称」
图8 求解信息:盯住最大残差、线性搜索与时间增量
✅ 弹塑性分析速查
1. 三段行为:<比例极限 线性|<屈服点 弹性(可全恢复)|>屈服点 塑性(有残余变形);比例极限与屈服点很近,有些软件视为相等
2. 版本坑:当前版本应力应变曲线的下降段输不进去(19 版之前可以)
3. 常用本构:双线性、随动强化、多线性强化
4. 切线模量:结构钢 800~8000 MPa;微小塑性应变(<0.5%)可到 100 MPa;通用规则 = 屈服强度的 8~10 倍(250 → 2000)
5. 材料输入:弹性模量与泊松比必须输入(与橡胶不同);数据用曲线拟合;输入的是塑性应变 + 真实应力
6. 2D 模型:必须 XY 平面、关于 Y 轴对称;平面应力 + 厚度 5mm
7. 网格:分区后要确认是映射网格,不是扫掠
8. 累积塑性应变看不了 → 到「输出」里把非线性打开
9. 结果项:EL = 弹性应变,EPPL = 塑性应变
10. 步数:不是越多越好,400 步反而可能算不出,步数少 + 自适应常常更稳
11. 带梁模型导入:线单元默认不确认,必须手动勾选,否则打开是空的且改不了类型
12. 提速:把不参与变形的零件改成刚体并重画网格
13. 收敛调参:非对称 + 大变形 + 线搜索 0.1 + 位移收敛 2%/100N + 弱弹簧;时间增量越小越易收敛。
【配套资料领取】
扫码添加小助理,私信回复关键词【ANSYS弹塑性本构】。领取本篇配套资抖包:①核心知识点速查;②操作步骤checklist;③参数公式与排错表










评论0