ANSYS弹塑性分析实战：从弹性到塑性收敛调试

简介：《2021年ANSYS弹性及塑性分析非常经典.docx》是一份面向ANSYS初、中级用户及机械、土木、航天等工程领域从业者的中文教程资料，聚焦弹性与塑性非线性分析的核心原理和软件操作方法，帮助读者理解材料屈服、塑性变形和ANSYS塑性选项的实际应用。整个压缩包仅含1个docx文档，大小108KB，内容以理论讲解、公式说明和操作流程为主，并以GUI与命令流两种方式提供了完整的塑性分析实例。文档从塑性基本概念、路径相关性与率相关性切入，对比工程应力应变与真实应力应变，随后系统介绍屈服准则、流动准则、强化准则以及BKIN、BISO、MKIN、MISO等常用塑性选项的参数设置，同时给出增强收敛性、查看结果等实用技巧，整体结构从基础概念到理论推导再到操作实例，便于读者边学边练。目前已有122人学习，适合初次接触ANSYS塑性分析的工程师快速建立知识框架，并按文档步骤完成带算例的实操练习。

1. 弹性及塑性分析：为什么一份经典案例能讲透ANSYS的取舍逻辑

一份被反复收藏的“2021年ANSYS弹性及塑性分析”案例材料，能成为很多结构工程师的入门经典，不是因为它包含了多复杂的装配体或非线性接触，而是它把两种力学性质完全不同的分析思路，放在同一个简单模型上完整跑了一遍。线弹性分析每个人都敢算，但算完心里没底；塑性分析理论上有数，但一调试就卡壳。这份材料真正的价值，是让你在同一个几何模型上同时看见弹性结果和塑性结果，再对着两条曲线思考结构到底该按哪种工况评估。

这篇实战笔记要解决的，就是这组问题：弹性分析怎么做才算完整，塑性分析又要打开哪些开关，同一结构两种算法差在哪，以及为什么你的弹塑性模型总是不收敛。适合刚接触有限元的应届生、做线性分析多年但没敢碰非线性的工程师，也适合需要给新人讲清楚弹塑性差异的仿真小组内部培训。

2. 先把线弹性分析跑通：前处理、求解与应力查看的最小闭环

线弹性分析是ANSYS里最便宜的一条路径，也是后续塑性分析的对标基线。它的前提是应力应变满足胡克定律，卸载后不留下任何残余变形，结构刚度恒定。绝大多数工程师拿到一个新结构，第一件事就是用线弹性算一遍，目的是快速定位高应力区，判断哪些位置需要更深入评估，再决定是否引入塑性分析。

2.1 材料参数与单元类型的选择依据

线弹性分析只需要两个材料常数：弹性模量和泊松比。常用钢材的弹性模量取2.06e5 MPa，泊松比0.3；铝合金一般取0.69e5 MPa和0.33。这里有一个容易忽略的前提：ANSYS本身不限定单位制，模型用mm、材料用MPa、力用N，得到的就是mm-MPa-N协调单位系；一旦几何用了mm、材料却填了Pa，结果会差出六个数量级。

单元类型的选择取决于结构形状。对于薄板受面内载荷这类问题，最常见的是PLANE183，这是一个8节点平面单元，支持平面应力、平面应变和轴对称三种行为。平面应力假设厚度方向的应力为零，适合薄板；平面应变假设厚度方向应变为零，适合厚截面或坝体类长结构。实体结构则常用SOLID185或SOLID186。三种单元的特性对比如下：

以一块带中心孔的平板为例：宽50mm、高50mm、厚5mm，孔半径5mm，左端约束，右端施加100MPa的均布拉应力。这是一个典型的平面应力问题，ANSYS经典界面前处理命令如下：

/PREP7 ET,1,PLANE183 ! 定义8节点平面单元 KEYOPT,1,3,2 ! 单元行为设为平面应力(2) MP,EX,1,2.06E5 ! 弹性模量，单位MPa MP,NUXY,1,0.3 ! 泊松比 RECTNG,0,50,0,50 ! 建立50x50mm矩形 CYL4,25,25,5 ! 建立半径5mm的孔，圆心(25,25) ASBA,1,2 ! 面1减去面2得到带孔板 ESIZE,1 ! 全局网格尺寸1mm AMESH,ALL ! 划分网格

KEYOPT,1,3,2这条命令决定了单元行为，数字3是单元关键项编号，末尾的2代表平面应力。如果把这里写错成平面应变，薄板拉板的应力结果会异常偏小。网格尺寸1mm在孔径10mm的模型上，孔边大约分布16个单元，足以捕捉基本应力集中。

2.2 网格密度对弹性结果的干扰与收敛性判断

线弹性分析中，网格密度的第一作用对象是应力峰值。在孔边、圆角根部、尖角位置，应力会随网格加密而单调上升，这就是常说的应力奇异现象。如果几何中存在理想尖角，无论网格细化到什么程度，应力都不会收敛到某个固定值。

判断网格是否收敛，我习惯的做法是：把全局网格尺寸从1mm改成0.5mm再算一遍，对比孔边应力峰值。如果两次结果差小于3%，基本可以接受；如果从310MPa跳到330MPa还在涨，说明网格未收敛，要继续细化或接受峰值偏大、改用第一层单元平均应力来评估。

在APDL里做局部细化比全局加密更实用，花费更少的单元数就能改善孔边应力分布：

LESIZE,6,,,32 ! 将孔边线段分成32份 AMESH,ALL ! 重新划分

这条命令把孔边线段细化为32个单元，而远离孔的区域仍保持1mm网格。注意线段编号6来自前处理的面生成顺序，不同版本可能不同，执行前最好先用LPLOT确认孔边的线号。

2.3 边界条件与载荷施加的稳定写法

边界条件施加方式直接影响求解稳定性。用线编号加载（DL、SFL）看上去简洁，但线号在不同建模顺序下会变。我一般更推荐用坐标选择节点来施加约束和载荷，这样命令可移植性更强：

FINISH /SOLU ALLSEL NSEL,S,LOC,X,0 ! 选择X=0处的所有节点 D,ALL,UX,0 ! 约束X方向位移 D,ALL,UY,0 ! 约束Y方向位移 ALLSEL NSEL,S,LOC,X,50 ! 选择X=50处节点 SF,ALL,PRES,-100 ! 施加100MPa拉应力，负号代表拉 SOLVE FINISH

这里NSEL,S,LOC,X,0按坐标选中左端面节点，ALLSEL在中间防止上次选择残留。SF命令给这些节点附加表面压力，正值表示压入表面，负值表示拉伸。如果不约束UY，平板存在刚体转动自由度，求解器会提示位移过大或直接算出荒谬结果；如果载荷方向填反，孔边应力云图将呈现负值异常。

求解完成后，进入后处理看Von Mises应力：

/POST1 PLNSOL,S,EQV ! 等效应力云图

云图红色区域集中在孔边上下两侧，这就是弹性分析给出的危险位置和应力峰值。记住这个值和危险区域，后面塑性分析就是拿它做对比的。

3. 从线性到非线性：塑性分析的核心物理量与屈服准则

结构局部应力超过屈服强度后，线弹性假设失效，应力和应变不再是唯一对应关系，卸载后也会留下永久变形。ANSYS处理塑性问题，需要引入非线性材料本构，并用增量迭代方式求解平衡方程，解到哪个载荷增量、材料处于哪个应力区间，都需要逐步追踪。

3.1 塑性分析必须打开的开关

塑性分析属于材料非线性，核心问题是材料本构描述。ANSYS中通过TB命令定义材料表，最常见的钢材弹塑性模型有两种：双线性各向同性强化模型（BISO）和双线性随动强化模型（BKIN）。两者都假设屈服前的弹性段和屈服后的强化段，只是强化行为的方向不同。BISO适合单调加载，屈服面等向扩大；BKIN考虑包辛格效应，适合往复加载，屈服面平移。

如果只做一次性强度分析，BISO够用；模拟反复加载或构件滞回行为，BKIN更贴近实际。还有一个开关经常被忽略：大变形选项。塑性是材料非线性，大变形是几何非线性，两者是独立开关。结构进入塑性后变形量通常增大，如果只定义材料本构却不开NLGEOM,ON，求解得到的位移场在量级上可能失真。

求解器也要从默认的线性求解切到牛顿迭代。APDL中用NROPT命令选择求解算法：

/SOLU NLGEOM,ON ! 打开大变形 NROPT,FULL,ON ! 全牛顿-拉弗森迭代

NROPT,FULL表示每一步更新刚度矩阵并做完全牛顿迭代，收敛更稳，代价是计算量大。对小模型这不是问题，但大模型需要考虑计算时间，可以改用修正牛顿法NROPT,MODI，更新频率降低、速度更快。

3.2 双线性材料模型参数的含义

仍以Q235钢为例，进入前处理器定义弹塑性材料：

/PREP7 MP,EX,1,2.06E5 MP,NUXY,1,0.3 TB,BISO,1,2 ! 定义双线性各向同性强化材料 TBDATA,1,235,2050 ! 屈服强度235MPa，切线模量2050MPa

TBDATA两个数值顺序不能颠倒：第一项是屈服应力，第二项是屈服后的切线模量。切线模量取自真实拉伸曲线屈服段之后的近似斜率，Q235取2050MPa已是一个保守估算，实际取值应参考材质证书或试验数据。如果TB,BISO下面跟3个数据，第三个是硬化参数，但对于双线性模型不需要。

用TBPLOT可以快速画出定义好的应力-应变曲线：

TBPLOT,BISO,1 ! 绘制材料本构曲线

如果曲线显示为两条直线且折点恰好位于235MPa附近，说明参数定义正确；如果载荷已加到100MPa但曲线仍是一条直线，多半是塑性表没有激活或材料表编号写错。

3.3 载荷步、子步与收敛准则的协同设置

塑性计算是否收敛，关键看载荷步和子步的分配。子步是载荷步内的增量段，ANSYS在每个子步内做平衡迭代。如果子步跨度过大，载荷增量一次越过屈服点，迭代法会在屈服点附近来回震荡；子步太小，计算时间成倍上涨。

推荐配置是打开自动时间步，并给出合理的子步上下限：

ALLSEL NSEL,S,LOC,X,0 D,ALL,UX,0 D,ALL,UY,0 ALLSEL NSEL,S,LOC,X,50 SF,ALL,PRES,-100 ! 施加目标载荷100MPa ALLSEL /SOLU AUTOTS,ON ! 自动时间步 NSUBST,20,100,5 ! 初始子步20，最大100，最小5 OUTRES,ALL,ALL ! 输出所有子步结果 SOLVE FINISH

NSUBST三个参数依次是初始子步数、最大子步数、最小子步数。自动时间步开启后，求解器会根据收敛情况动态调整步长：某子步迭代不收敛，自动劈小再试；连续收敛顺利，则自动增大子步加速。OUTRES,ALL,ALL是另一个容易被漏掉的设置，不写这条，后处理里很可能只看到最后一个子步的结果，塑性发展过程全丢。

4. 弹性与塑性结果对比：应力应变解读的差异与等效塑性应变

同一个带孔板，弹性分析和塑性分析会在后处理中呈现完全不同的画面。弹性结果的应力峰值往往很高，看起来吓人；塑性结果的应力峰值被"削平"，应变却明显增大。理解这组差异，才能对结构做出正确评估。

4.1 后处理中必须区分的三类应变

塑性分析后，总应变被分解为弹性应变和塑性应变两部分。等效塑性应变（Equivalent Plastic Strain）是塑性区累积损伤程度的关键指标，它不区分方向，只表征材料不可恢复变形的总量。

后处理命令的区别如下：

/POST1 PLNSOL,S,EQV ! Von Mises应力 PLNSOL,EPPL,EQV ! 等效塑性应变 PLNSOL,EPEL,EQV ! 等效弹性应变

一个典型结果是：孔边的等效应力显示在235MPa附近，不再是弹性分析算出的310MPa；等效塑性应变云图在孔边出现明显的红色高值区，向外迅速衰减。这代表孔边材料已屈服，原有线弹性解算出的高应力被塑性区重新分配，荷载由周围未屈服的材料分担。

产生这一现象后，如果依然用弹性分析结果评估强度，就会得到过于保守的结论，有些构件甚至因为弹性应力超标被误判为不合格。正确做法是对塑性区范围做检查：若塑性区局限在孔边很小的范围且等效塑性应变值低于材料断裂延伸率，结构整体仍然安全；若塑性区贯穿整个截面，则接近极限承载力状态。

4.2 用荷载-位移曲线判断塑性发展全过程

云图只能展示最终加载状态。想了解结构从弹性进入塑性的拐点，需要绘制荷载-位移曲线。先找到孔边或位移最大处的节点号，用POST26时间历程后处理提取该节点的位移，与施加的载荷值对应绘制。

操作可以这样做：

/POST26 NSOL,2,9,U,X ! 提取节点9的X方向位移

执行前先用NSORT或GUI拾取确认节点的实际编号。将节点位移作为横坐标、提取到的反力或应力作为纵坐标，曲线形态会呈现两段：前段是直线，斜率为结构的弹性刚度；后段发生弯折，斜率变缓，弯折点对应的荷载就是结构进入塑性状态的临界值。

这一弯折点通常低于线弹性分析中应力峰值达到屈服应力的那个载荷。因为局部屈服发生在靠近孔边的点，随后塑性区逐步扩展，最终使整个承载截面屈服。用这个曲线配合前面看过的等效塑性应变云图，能够直观判断结构的安全裕度和塑性储备，结论远比单看应力云图完整。

5. ANSYS弹性与塑性分析的避坑指南：翻车现场与排查思路

弹性分析容易在应力奇异上翻车，塑性分析容易在收敛和材料本构上翻车。下面整理几条高频踩坑记录，每一条都按现象、原因、解决的思路给出，方便你在自己的模型里对照排查。

5.1 应力越算越大：孔边应力峰值发散

现象：线弹性分析中，孔边应力峰值随网格加密持续上升，无论怎么做结果都稳定不下来。原因：几何中存在真实尖角或者载荷集中到单个节点，理论解在该点本身就是无穷大，有限元只是把无穷大逼近得越来越激进。解决：先判断这个位置在物理上是否真的存在理想尖角。如果只是建模简化，把尖角倒成小圆角再算；如果必须保留尖角，则不要读该点应力峰值，改读距尖角一个单元距离处的应力或采用第一层单元平均应力。

5.2 塑性分析一直不收敛：迭代发散排查顺序

现象：塑性分析能启动，但反复提示不平力超限或节点位移过大，自动时间步不断缩小仍无法收敛。原因：常见触发因素有四个，一是初始子步数太少导致载荷增量过大，二是有刚体位移没被约束，三是材料本构参数填错，四是切线模量填了负值。解决：先检查约束是否完整，用过约束查询确认模型没有自由度遗漏；再把AUTOTS关掉，用固定子步数50试算；如果仍然发散，检查TB参数，把切线模量改成0试试，看是否为材料定义问题。

5.3 塑性结果和弹性结果一模一样：载荷步设置不对

现象：材料定义了BISO，也开了大变形，算完发现应力云图和线弹性结果完全相同，塑性完全没有激活。原因：载荷只设了一个步，且子步数正好被求解器在第一步就全部加载完成，结构在第一个增量步内直接越过屈服点，迭代过程没有捕捉到塑性状态。解决：把载荷分成两到三个载荷步，第一步加载到80%目标载荷，第二步再加载到100%，配合NSUBST设置。这样材料在第一个载荷步内稳定进入屈服状态，后续迭代才会产生塑性应变。

5.4 材料曲线明显错误：参数量纲混用

现象：TBPLOT画出来的应力应变曲线在某个点剧烈跳变，或者屈服应力数值大得离谱，塑性阶段斜率陡升。原因：常见问题是弹性模量用了Pa而屈服强度用了MPa，两者相差三个数量级，导致本构曲线形状不合理。解决：养成检查量纲的习惯，弹性模量与屈服强度单位必须一致。还有一个快速核对方法：弹性段的应变值应该等于屈服强度除以弹性模量，比如235MPa/206000MPa约等于0.00114，如果曲线折点对应的应变值比这个数大出几十倍，材料参数一定填错了。

5.5 刚体位移提示但模型明明有约束

现象：求解过程提示"小的负主元"或位移过大，但检查边界约束时发现该约束的位置都有施加。原因：约束施加在了坐标选择后的节点子集上，后来某条ALLSEL命令可能把已施加的约束覆盖或选择范围外漏掉了关键节点。解决：求解前用D命令列表(DLIST)确认约束节点数量是否正确，用FSUM检查载荷合力是否平衡。这个习惯可以在大型装配体分析时省掉许多排查时间。

6. 收敛验证与参数微调：一个算例的完整调试过程

想快速验证ANSYS弹塑性配置是否正常，我最常用的是一个单杆拉伸算例，每次换新环境或新版本都先跑一遍它，几十秒就能确认材料定义、求解器设置、子步策略是否全部正确。

取杆长10mm，截面1mm×1mm，一端固定，另一端沿轴向拉伸2mm，材料为理想弹塑性，弹性模量200000MPa，屈服强度200MPa。理论结果是：应力先随应变线性上升，到200MPa后保持恒定，不再随加载位移增加。因为理想塑性模型的切线模量为0，应变0.2时总应变远大于屈服应变，应力必然停留在屈服值200MPa，杆件轴力为200N。

测试命令流如下：

/PREP7 ET,1,LINK180 ! 3D杆单元 MP,EX,1,2.0E5 MP,NUXY,1,0.3 TB,BISO,1,2 TBDATA,1,200,0 ! 理想弹塑性，切线模量为0 N,1,0,0,0 N,2,10,0,0 E,1,2 D,1,ALL,0 D,2,UY,0 D,2,UZ,0 D,2,UX,2 ! 拉伸位移2mm FINISH /SOLU NLGEOM,ON AUTOTS,ON NSUBST,50,200,20 SOLVE FINISH

求解完成后，后处理看杆的反力：

/POST1 ESOL,2,1,1,F,X ! 提取单元1在节点1的X方向反力 *GET,F1,NODE,1,RF,FX ! 保存反力值

如果提取到的反力约为200N，说明材料定义和求解路径正确；如果反力是0或者远小于200N，问题多半出在约束缺失或TB参数未激活。这个算例的另一个用途是摸索不同版本的推荐子步区间。把NSUBST从50分别改成10和200各跑一次，观察应力结果的数值抖动程度，就能建立新环境下的参数手感，之后再把它套用到含孔板这类稍大模型上，会少走很多弯路。

做完一个弹塑性项目后回头整理，我现在的习惯是至少保留三样东西：弹性基线结果、塑性最终云图、危险节点荷载位移曲线。工程评审时被问到为什么某个位置算出来应力没超标却担心失效，直接调出等效塑性应变云图和数据曲线，比空口解释“塑性重分布”有说服力得多。希望帮到你。

本文还有配套的精品资源，点击获取