简介ANSYS仿真分析PPT学习教案是一份面向工程技术人员与高校学生的有限元入门教学资料聚焦仿真前处理环节帮助读者理清静力/动力、线性/非线性分析的区别以及如何根据载荷特点选择分析类型。资源为1个PPTX文件共69页整个压缩包体积仅606KB内容却非常浓缩。教案系统讲解了模型精度控制方法离散误差、边界条件、单元形状并展开介绍结构降维、细节简化、形式变换、局部结构、对称性利用等实用处理手段同时覆盖各向同性、各向异性及复合材料的特点以及模型建立时单元选择、实常数设置、载荷施加限制等容易出错的要点有助于降低建模风险。已有101人学习下载适合希望在ANSYS中建立准确、高效分析模型的初学者参考对刚接触ANSYS的读者尤其友好。1. 先判断再建模静力/动力与线性/非线性怎么选很多刚接触 ANSYS 仿真分析的人拿到一个结构就急着导入几何、划分网格结果算出来的应力大得离谱还以为是网格不够密。实际排查下来大多数问题出在分析类型选错上把时间相关的振动载荷按静力算或者把已经进入塑性的结构按线弹性处理结果差一个数量级都很正常。ANSYS 分析前的准备阶段本质上是在回答四个问题载荷是否随时间快速变化结构刚度会不会因加载而显著改变关心的结果是变形还是应力以及需要多高的精度。这套判断逻辑并不复杂但漏掉任何一环后续所有操作都会建立在错误的地基上。这篇内容把分析类型选择、误差来源、结构简化、单元选型和网格控制这几条主线串起来适合刚入门的学生也适合做过几个例子但没系统梳理过建模决策的工程师。2. 误差控制离散误差、边界误差与计算误差的连锁效应2.1 模型误差的三个来源离散、边界与形状2.1.1 离散误差单元尺寸与阶次的博弈有限元方法用有限数量的单元去逼近连续体这种逼近本身就会带来离散误差。单元把原本光滑的位移场用多项式近似单元尺寸越大近似误差越大单元阶次越低对弯曲、应力梯度这类复杂响应的还原能力越差。典型的表现是细化网格后应力明显变化说明之前的结果没有收敛换用高阶单元后结果显著不同说明原单元阶次不足以描述真实变形模式。离散误差的控制手段很直接减小单元尺寸或提高单元阶次。但有代价计算量上升累计误差也可能增加所以存在一个甜点区间。ANSYS 里建立对比模型是常用的验证方式下面这组命令可以快速比较一阶和二阶平面单元在同一模型上的差异/PREP7 ET,1,PLANE182 ! 4节点一阶平面单元 ET,2,PLANE183 ! 8节点二阶平面单元 MP,EX,1,2.1E5 ! 弹性模量单位 MPa MP,PRXY,1,0.3 ! 泊松比 RECTNG,0,50,0,20 ! 50x20 的平面模型 ESIZE,2 ! 单元尺寸 2mm AMESH,ALL ! 划分平面网格这段命令先定义两种单元类型和材料属性再用矩形面建模并指定全局网格尺寸。PLANE182 是四节点单元每个节点 2 个平移自由度计算快但抗弯曲能力弱PLANE183 是八节点二次单元边界可以贴合曲线对弯曲主导的问题收敛速度快得多。实际操作中我会先用 PLANE182 跑一遍粗网格作为基线再用 PLANE183 加密一次两次结果的偏差如果超过 5%说明离散误差还没控制住。2.1.2 边界误差和单元形状误差边界误差来自约束和载荷施加位置与实际工况不符。常见问题包括把螺栓连接简化为全约束导致连接处刚度偏大把均布载荷简化成集中力导致局部应力奇异对称面取在最大应力位置导致峰值被约束硬化。这类误差和网格密度无关加密反而会让错误更明显因为约束处的应力集中会被放大。单元形状误差则是网格质量的问题。内角接近 0° 或 180° 的扁平单元、过度扭曲的单元会让刚度矩阵出现病态求解器要么报错要么给出失去物理意义的振荡结果。网格划分后应该检查长宽比和翘曲度指标ANSYS 在 Mesh Tool 里有单元质量检查选项也可以在 Preprocessor 里用SHPP命令设定形状检查阈值。2.2 计算误差舍入与截断在规模化后的代价计算误差和建模无关但会反过来限制你对模型规模的追求。舍入误差来自计算机浮点数的有限精度当模型有几十万甚至上百万自由度时求解过程中的加减运算是逐级累积的节点越多浮点误差叠加的次数越多。截断误差则来自求解算法的迭代终止条件ANSYS 默认的力收敛容差是 0.5%如果为了保证精度把容差设到 1e-6迭代次数会飙升而有效数字的收益在达到机器精度后会停止。这里有个反直觉的结论单元数和节点数不是越多越好。加密到一定程度后离散误差的下降速度会变慢而舍入误差的累积会持续存在两者叠加可能出现结果随网格加密反而轻微波动的现象。所以在工程实践中我一般用结果是否随网格变化收敛作为标准而不是一味追求最小网格尺寸。2.3 控制精度的四条路径与一组提取命令误差来源典型表现控制手段离散误差网格加密后结果持续大幅变化减小单元尺寸、提高单元阶次边界误差约束处应力异常集中或偏软按实际工况简化边界避免过约束单元形状误差局部应力出现锯齿状分布细化过渡区、重划畸形单元计算误差高网格密度下结果轻微漂移控制模型规模检查收敛容差设置误差控制路径总结下来就是提高单元阶次、增加单元数量、划分规则单元形状、建立与实际工况相符的边界条件。四条路径互相关联但优先级不同。边界条件错误时加密网格只会放大错误单元形状差时提高阶次收益也有限。求解完成后用下面这组命令快速提取最大等效应力/POST1 SET,LAST NSORT,S,EQV ! 按等效应力排序 *GET,SMAX,SORT,0,MAX ! 提取排序最大值 PRNSOL,S,COMP ! 打印节点应力分量NSORT把节点结果按等效应力降序排列*GET从排序结果里拿出最大值PRNSOL输出完整应力分量。判断误差是否可接受时我习惯把提取出的 SMAX 和网格尺寸画成曲线曲线进入平台的区域就是合理的网格密度区间。3. 结构简化的五条路径降维、简化、变换、局部与对称3.1 降维处理什么时候可以把实体变成杆梁或平面三维实体模型的计算成本远高于二维和线单元。降维处理的本质是判断结构的几何特征和受力模式允许在哪个方向上做简化。细长杆件只承受轴向力时可以用杆单元每个节点只有一个轴向自由度完全忽略弯曲梁构件承受弯矩剪力时用梁单元但要额外处理截面方位、节点偏移和自由度释放薄壁结构的面内受力可以用平面应力单元垂直于面的应力分量为零大坝、隧道这类长度方向不变的厚壁结构则适用平面应变假设。降维不是随便拍脑袋的。判断依据是尺寸比例长细比大于 10 的构件适合梁杆模型厚度远小于另外两个方向尺寸的适合壳或平面单元。ANSYS 里梁单元的截面方位由方向点或截面坐标控制节点偏移处理偏心连接自由度释放用于模拟铰接这些参数一个没设对结果就会偏离实际。3.2 细节简化与形式变换把几何复杂度转为分析效率细节简化的目标是把对整体分析影响不大的特征去掉典型对象包括远离关键部位的倒角、小孔、螺纹和凸台。这些特征的网格尺度往往只有主体结构的百分之一不简化的话局部网格会拖垮整个求解速度而且它们对整体刚度和远端应力的贡献可以忽略。判断标准有两个离关心区域的距离以及对整体刚度的贡献比例。形式变换解决的是另一类问题有些实体形状不规则直接划分网格会生成大量畸形单元不如把几何形式转换掉。教案里提到的做法是把加强筋转换为平面单元这样在保留主要传力路径的同时避免了三维网格在薄壁处的长细比失控。做这类变换时要注意等效刚度匹配筋的高度、宽度在转换成平面单元后要通过实常数或厚度参数还原。3.3 局部结构切割关心的集中受力区大型结构里如果关心的只是某个局部区域的应力没有必要对整个结构做精细网格。常见做法是把局部区域从整体结构中切割出来用切割面处的内力或位移结果作为边界条件施加到局部模型上。这样可以大幅减小模型规模同时把网格资源集中在真正需要高精度的地方。局部模型需要注意切割面到关心区的距离一般要大于一到两个特征尺寸否则切割面上的圣维南效应会污染目标区域的应力结果。3.4 对称性利用的边界条件与三个注意点对称性利用可以把模型缩小到原来的 1/2、1/4 甚至 1/8是降低计算成本最有效的手段之一。但对称分析有固定的前提几何对称、载荷对称、约束对称、材料对称四者缺一不可。满足前提后在对称面上施加对称约束边界条件即可。APDL 中的施加方式是DL,2, ,SYMMDL表示在线上施加自由度约束2是线编号SYMM是对称约束标记。这个命令会把线上节点的法向位移约束为零同时放开切向自由度和固定约束有本质区别不要混用。3.4.1 对称面上的载荷、刚度与最大应力教案里专门指出对称分析的三个注意点每一个都在实际项目里踩过坑。第一如果对称面上有作用的载荷对称分析时取载荷的 1/2否则等效载荷翻倍。第二若对称面上存在板或梁离散板和梁的单元所有节点均位于对称面上这时板的刚度应取整个单元刚度的 1/2而不是取 1/2 个单元的全部强度这两个做法虽然听着接近但结果完全不同。第三对称面不应该取在最大应力处最大应力点落在对称面上时对称约束会人为改变该点的变形模式应力结果失真。3.4.2 镜像命令在曲面上的禁区对称建模还容易犯一个错误用镜面生成对称部分。教案明确警告不能使用 ARSYSM、LSYMM 这类镜像技术来映射圆、圆柱、圆锥或球面原因是每个实常数的设置不能同时赋给多个基本原型段。曲面镜像后实常数对应的坐标系和节点方向会错乱壳厚、梁截面方向都可能在镜像后反转。对称部分应该通过复制加旋转或者直接重新建模实现而不是简单镜像。这一章的五种方法可以按下面的思路选用处理方法适用场景关键注意点降维处理长细比大的构件、薄壁结构梁截面方位、节点偏移细节简化远离关心区的小特征评估对整体刚度的影响形式变换不规则实体难划分网格等效刚度匹配局部结构大型结构局部高应力区切割面距离足够远对称性利用几何、载荷、约束全对称载荷取半、刚度减半、避开最大应力面4. 材料、单元与属性选型之前先看限制4.1 材料分类各向同性、各向异性与复合材料ANSYS 分析中的材料属性直接决定本构关系选错材料模型比单元选错更致命。各向同性材料在任意一点沿任何方向的力学、热学性能均相同所有金属材料都属此类ANSYS 中只需要定义弹性模量 EX、泊松比 PRXY 和密度 DENS 三个基本参数。各向异性材料的性能方向相关木材顺着纹理和横跨纹理的弹性模量可以差一个数量级合成纤维复合材料更是如此必须在材料定义中给出多组弹性常数。复合材料则是两种及以上材料的组合教案里将其分为单层和多层多层板在 ANSYS 中通常用壳单元的分层截面定义每一层有自己的厚度、材料和铺层角度。材料参数对仿真分析的影响在热力耦合场景里尤其明显如果只是做结构分析线弹性各向同性材料定义三行就够MP,EX,1,2.1E5 ! 弹性模量 210GPa MP,PRXY,1,0.3 ! 泊松比 MP,DENS,1,7.85E-9 ! 密度单位 t/mm^3MP是材料属性定义命令格式为MP,属性,材料编号,数值。注意单位体系的一致性如果在 mm 单位制下弹性模量用 MPa密度用 t/mm³力的单位就是 N。4.2 单元使用限制棱边载荷、面载荷、力矩与扭矩每种单元都有它的不能做的事教案里列出的限制非常实用实体单元不能施加棱边载荷因为实体单元的形函数只支持面力棱边载荷无法形成等效节点力轴对称单元不能施加面载荷只能施加在对应的线段上杆单元没有旋转自由度不能施加节点力矩和扭矩否则会得到零刚度响应。单元限制细节单元类型载荷限制建模提示实体单元不能施加棱边载荷载荷只能加在面或节点上轴对称单元不能施加面载荷载荷施加在对应线段杆单元不能施加节点力矩和扭矩只有轴向自由度无抗弯能力梁单元截面方位、节点偏移需显式控制自由度释放模拟铰接板单元不同厚度连接处节点必须重合网格过渡需要对齐节点板单元这里有个常见失误不同厚度的板连接时如果网格节点不重合计算结果会在连接处出现虚假的应力跳跃。解决办法是在厚度变化的边界处预先设置网格种子点或者用绑结接触把不匹配的网格连接起来。4.3 自由度协调壳单元 ROTZ 的虚拟刚度问题不同单元类型连接时自由度协调是一个隐蔽但致命的坑。杆单元只有轴向平动自由度梁单元有平动加转动自由度壳单元除了面内平动还有一个绕面法线的旋转自由度 ROTZ。问题在于壳单元的 ROTZ 与平面旋转刚度有关而这个刚度是虚拟的不是真实物理刚度所以梁单元和薄壳单元相连时对应的自由度并不协调可能出现局部零刚度求解器直接报奇异。教案里提到 SHELL43 和 SHELL63 单元在KEYOPT(3)2时会激活所有旋转自由度属于例外情况。实际建模里处理梁壳连接我的习惯是优先考虑几何上的合理过渡要么在连接区域使用实体单元做转接要么用 MPC 绑定约束把梁端节点和壳面节点关联起来这两种方式都比强行共用节点可靠。另外两个单元使用细节单元带中间节点时划分网格应让中间节点对齐否则单元边界的位移连续性被破坏相邻单元应当有相同的单元边节点数一边 3 节点一边 2 节点会导致变形不协调。关于非线性分析二次单元的积分点数量不一定比线性单元多但由于位移插值精度更高对接触和大变形问题的收敛性有实质帮助所以在塑性、接触、大变形场景下优先选用二次单元。4.4 单元属性在划分网格之前全部指定单元属性包括单元类型、实常数、材料编号和截面定义。很多新手习惯在生成网格后再回去补材料这样经常出现部分单元使用缺省属性导致结果异常。正确的流程是先定义好所有单元类型、材料、实常数划分网格前通过 Mesh Attributes 给每一类几何实体分配正确的属性。教案里强调了两点在实体模型上直接指定属性时将不考虑缺省属性这可以避免在网格划分操作中重新设置属性而且 ANSYS 的网格划分算法在一次对所有实体进行网格划分时效率更高清除实体模型上的网格不会删除指定的单元属性所以调整网格参数时不需要担心属性丢失。划分完网格后应该激活属性编号校核路径是 Utility Menu PlotCtrls Numbering把单元类型、材料编号、实常数编号都打开目视检查有没有属性错误的区域。壳单元和梁单元的截面定义示例ET,1,SHELL181 ! 壳单元 SECTYPE,1,SHELL ! 定义壳截面 SECDATA,2 ! 壳厚度 2mm ET,2,BEAM189 ! 梁单元 SECTYPE,2,BEAM,RECT ! 矩形梁截面 SECDATA,10,20 ! 截面宽 10mm、高 20mmSECTYPE声明截面类型和编号SECDATA输入截面几何参数。壳单元厚度方向的积分点个数直接影响弯曲应力计算精度积分点太少会低估壳表面的弹塑性响应太多则增加计算代价默认值对大多数薄壳分析是够用的。5. 网格划分的控制原则与无关性验证5.1 网格划分三步与三组选择原则网格划分在 ANSYS 里遵循固定三步定义单元属性、指定网格控制参数、生成网格。对应菜单路径分别是 Preprocessor Meshing Mesh Attributes、Preprocessor Meshing Size Control、以及最终的 Mesh 操作。在 Workbench 里这套流程对应为 Mesh 模块下的 Element Size、Local Sizing 和 Generate Mesh逻辑一致只是换成了参数化界面。网格粗细没有统一答案教案给出了三组判断原则。静力分析中只关心变形时可以画较少的网格因为位移是能量意义上的整体量对局部误差不敏感要计算应力或应变时则需要较多网格因为应力由位移求导得到对误差更敏感。固有属性分析中只算少数低阶模态可以选较少的网格高阶模态的振型复杂需要足够多的单元分辨局部变形。响应分析同理只关注某些部位的位移响应用粗网格即可要算应力响应则必须加密。分析场景粗网格可用必须加密静力分析只计算变形计算应力或应变模态分析只关心低阶模态计算高阶模态响应分析只关心位移响应计算应力响应网格无关性验证的操作流程是选定一组从粗到细的网格尺寸比如 10mm、5mm、2.5mm、1.25mm在每种网格下提取关心点的最大等效应力和最大位移比较相邻两组结果的变化率。步骤如下建立模型并确定目标结果如关键节点的等效应力峰值用最大网格尺寸求解一次提取结果逐次加密网格保持边界条件和载荷完全一致当相邻两次的应力变化小于 3% 时认为结果与网格无关。5.2 网格无关性验证细化到什么程度才合适加密网格到一定程度后结果会进入收敛平台继续细化只会增加求解时间甚至会因为累计误差出现轻微漂移。这时就该停止细化采用当前网格尺寸。每次求解完成后可以用下面的 APDL 命令提取关键结果/POST1 SET,LAST NSORT,S,EQV ! 按等效应力排序 *GET,SMAX,SORT,0,MAX ! 取最大等效应力 NSORT,U,SUM ! 按总位移排序 *GET,UMAX,SORT,0,MAX ! 取最大位移SET,LAST读取最后一个载荷步的结果NSORT对节点结果排序*GET从排序结果中提取极值。把两列数据记录成下面的格式对比网格尺寸 mm最大等效应力 MPa最大位移 mm应力变化率10248.50.532-5265.30.5476.8%2.5270.10.5511.8%1.25271.20.5520.4%到 2.5mm 时应力变化率已经降到 1.8%1.25mm 时降到 0.4%说明 2.5mm 的结果已经进入收敛平台。验证完成后实际生产计算采用 2.5mm 网格即可。如果关心的点恰好位于应力集中区还要顺带检查该点的等效应力是否落在不同单元尺寸下都保持稳定而不是只看全局最大值位置有没有漂移。网格无关性验证的目标不是追求绝对精确而是在计算成本和结果精度之间找到那个拐点。本文还有配套的精品资源点击获取 ANSYS仿真分析：从误差控制到网格无关性验证的建模要点
简介ANSYS仿真分析PPT学习教案是一份面向工程技术人员与高校学生的有限元入门教学资料聚焦仿真前处理环节帮助读者理清静力/动力、线性/非线性分析的区别以及如何根据载荷特点选择分析类型。资源为1个PPTX文件共69页整个压缩包体积仅606KB内容却非常浓缩。教案系统讲解了模型精度控制方法离散误差、边界条件、单元形状并展开介绍结构降维、细节简化、形式变换、局部结构、对称性利用等实用处理手段同时覆盖各向同性、各向异性及复合材料的特点以及模型建立时单元选择、实常数设置、载荷施加限制等容易出错的要点有助于降低建模风险。已有101人学习下载适合希望在ANSYS中建立准确、高效分析模型的初学者参考对刚接触ANSYS的读者尤其友好。1. 先判断再建模静力/动力与线性/非线性怎么选很多刚接触 ANSYS 仿真分析的人拿到一个结构就急着导入几何、划分网格结果算出来的应力大得离谱还以为是网格不够密。实际排查下来大多数问题出在分析类型选错上把时间相关的振动载荷按静力算或者把已经进入塑性的结构按线弹性处理结果差一个数量级都很正常。ANSYS 分析前的准备阶段本质上是在回答四个问题载荷是否随时间快速变化结构刚度会不会因加载而显著改变关心的结果是变形还是应力以及需要多高的精度。这套判断逻辑并不复杂但漏掉任何一环后续所有操作都会建立在错误的地基上。这篇内容把分析类型选择、误差来源、结构简化、单元选型和网格控制这几条主线串起来适合刚入门的学生也适合做过几个例子但没系统梳理过建模决策的工程师。2. 误差控制离散误差、边界误差与计算误差的连锁效应2.1 模型误差的三个来源离散、边界与形状2.1.1 离散误差单元尺寸与阶次的博弈有限元方法用有限数量的单元去逼近连续体这种逼近本身就会带来离散误差。单元把原本光滑的位移场用多项式近似单元尺寸越大近似误差越大单元阶次越低对弯曲、应力梯度这类复杂响应的还原能力越差。典型的表现是细化网格后应力明显变化说明之前的结果没有收敛换用高阶单元后结果显著不同说明原单元阶次不足以描述真实变形模式。离散误差的控制手段很直接减小单元尺寸或提高单元阶次。但有代价计算量上升累计误差也可能增加所以存在一个甜点区间。ANSYS 里建立对比模型是常用的验证方式下面这组命令可以快速比较一阶和二阶平面单元在同一模型上的差异/PREP7 ET,1,PLANE182 ! 4节点一阶平面单元 ET,2,PLANE183 ! 8节点二阶平面单元 MP,EX,1,2.1E5 ! 弹性模量单位 MPa MP,PRXY,1,0.3 ! 泊松比 RECTNG,0,50,0,20 ! 50x20 的平面模型 ESIZE,2 ! 单元尺寸 2mm AMESH,ALL ! 划分平面网格这段命令先定义两种单元类型和材料属性再用矩形面建模并指定全局网格尺寸。PLANE182 是四节点单元每个节点 2 个平移自由度计算快但抗弯曲能力弱PLANE183 是八节点二次单元边界可以贴合曲线对弯曲主导的问题收敛速度快得多。实际操作中我会先用 PLANE182 跑一遍粗网格作为基线再用 PLANE183 加密一次两次结果的偏差如果超过 5%说明离散误差还没控制住。2.1.2 边界误差和单元形状误差边界误差来自约束和载荷施加位置与实际工况不符。常见问题包括把螺栓连接简化为全约束导致连接处刚度偏大把均布载荷简化成集中力导致局部应力奇异对称面取在最大应力位置导致峰值被约束硬化。这类误差和网格密度无关加密反而会让错误更明显因为约束处的应力集中会被放大。单元形状误差则是网格质量的问题。内角接近 0° 或 180° 的扁平单元、过度扭曲的单元会让刚度矩阵出现病态求解器要么报错要么给出失去物理意义的振荡结果。网格划分后应该检查长宽比和翘曲度指标ANSYS 在 Mesh Tool 里有单元质量检查选项也可以在 Preprocessor 里用SHPP命令设定形状检查阈值。2.2 计算误差舍入与截断在规模化后的代价计算误差和建模无关但会反过来限制你对模型规模的追求。舍入误差来自计算机浮点数的有限精度当模型有几十万甚至上百万自由度时求解过程中的加减运算是逐级累积的节点越多浮点误差叠加的次数越多。截断误差则来自求解算法的迭代终止条件ANSYS 默认的力收敛容差是 0.5%如果为了保证精度把容差设到 1e-6迭代次数会飙升而有效数字的收益在达到机器精度后会停止。这里有个反直觉的结论单元数和节点数不是越多越好。加密到一定程度后离散误差的下降速度会变慢而舍入误差的累积会持续存在两者叠加可能出现结果随网格加密反而轻微波动的现象。所以在工程实践中我一般用结果是否随网格变化收敛作为标准而不是一味追求最小网格尺寸。2.3 控制精度的四条路径与一组提取命令误差来源典型表现控制手段离散误差网格加密后结果持续大幅变化减小单元尺寸、提高单元阶次边界误差约束处应力异常集中或偏软按实际工况简化边界避免过约束单元形状误差局部应力出现锯齿状分布细化过渡区、重划畸形单元计算误差高网格密度下结果轻微漂移控制模型规模检查收敛容差设置误差控制路径总结下来就是提高单元阶次、增加单元数量、划分规则单元形状、建立与实际工况相符的边界条件。四条路径互相关联但优先级不同。边界条件错误时加密网格只会放大错误单元形状差时提高阶次收益也有限。求解完成后用下面这组命令快速提取最大等效应力/POST1 SET,LAST NSORT,S,EQV ! 按等效应力排序 *GET,SMAX,SORT,0,MAX ! 提取排序最大值 PRNSOL,S,COMP ! 打印节点应力分量NSORT把节点结果按等效应力降序排列*GET从排序结果里拿出最大值PRNSOL输出完整应力分量。判断误差是否可接受时我习惯把提取出的 SMAX 和网格尺寸画成曲线曲线进入平台的区域就是合理的网格密度区间。3. 结构简化的五条路径降维、简化、变换、局部与对称3.1 降维处理什么时候可以把实体变成杆梁或平面三维实体模型的计算成本远高于二维和线单元。降维处理的本质是判断结构的几何特征和受力模式允许在哪个方向上做简化。细长杆件只承受轴向力时可以用杆单元每个节点只有一个轴向自由度完全忽略弯曲梁构件承受弯矩剪力时用梁单元但要额外处理截面方位、节点偏移和自由度释放薄壁结构的面内受力可以用平面应力单元垂直于面的应力分量为零大坝、隧道这类长度方向不变的厚壁结构则适用平面应变假设。降维不是随便拍脑袋的。判断依据是尺寸比例长细比大于 10 的构件适合梁杆模型厚度远小于另外两个方向尺寸的适合壳或平面单元。ANSYS 里梁单元的截面方位由方向点或截面坐标控制节点偏移处理偏心连接自由度释放用于模拟铰接这些参数一个没设对结果就会偏离实际。3.2 细节简化与形式变换把几何复杂度转为分析效率细节简化的目标是把对整体分析影响不大的特征去掉典型对象包括远离关键部位的倒角、小孔、螺纹和凸台。这些特征的网格尺度往往只有主体结构的百分之一不简化的话局部网格会拖垮整个求解速度而且它们对整体刚度和远端应力的贡献可以忽略。判断标准有两个离关心区域的距离以及对整体刚度的贡献比例。形式变换解决的是另一类问题有些实体形状不规则直接划分网格会生成大量畸形单元不如把几何形式转换掉。教案里提到的做法是把加强筋转换为平面单元这样在保留主要传力路径的同时避免了三维网格在薄壁处的长细比失控。做这类变换时要注意等效刚度匹配筋的高度、宽度在转换成平面单元后要通过实常数或厚度参数还原。3.3 局部结构切割关心的集中受力区大型结构里如果关心的只是某个局部区域的应力没有必要对整个结构做精细网格。常见做法是把局部区域从整体结构中切割出来用切割面处的内力或位移结果作为边界条件施加到局部模型上。这样可以大幅减小模型规模同时把网格资源集中在真正需要高精度的地方。局部模型需要注意切割面到关心区的距离一般要大于一到两个特征尺寸否则切割面上的圣维南效应会污染目标区域的应力结果。3.4 对称性利用的边界条件与三个注意点对称性利用可以把模型缩小到原来的 1/2、1/4 甚至 1/8是降低计算成本最有效的手段之一。但对称分析有固定的前提几何对称、载荷对称、约束对称、材料对称四者缺一不可。满足前提后在对称面上施加对称约束边界条件即可。APDL 中的施加方式是DL,2, ,SYMMDL表示在线上施加自由度约束2是线编号SYMM是对称约束标记。这个命令会把线上节点的法向位移约束为零同时放开切向自由度和固定约束有本质区别不要混用。3.4.1 对称面上的载荷、刚度与最大应力教案里专门指出对称分析的三个注意点每一个都在实际项目里踩过坑。第一如果对称面上有作用的载荷对称分析时取载荷的 1/2否则等效载荷翻倍。第二若对称面上存在板或梁离散板和梁的单元所有节点均位于对称面上这时板的刚度应取整个单元刚度的 1/2而不是取 1/2 个单元的全部强度这两个做法虽然听着接近但结果完全不同。第三对称面不应该取在最大应力处最大应力点落在对称面上时对称约束会人为改变该点的变形模式应力结果失真。3.4.2 镜像命令在曲面上的禁区对称建模还容易犯一个错误用镜面生成对称部分。教案明确警告不能使用 ARSYSM、LSYMM 这类镜像技术来映射圆、圆柱、圆锥或球面原因是每个实常数的设置不能同时赋给多个基本原型段。曲面镜像后实常数对应的坐标系和节点方向会错乱壳厚、梁截面方向都可能在镜像后反转。对称部分应该通过复制加旋转或者直接重新建模实现而不是简单镜像。这一章的五种方法可以按下面的思路选用处理方法适用场景关键注意点降维处理长细比大的构件、薄壁结构梁截面方位、节点偏移细节简化远离关心区的小特征评估对整体刚度的影响形式变换不规则实体难划分网格等效刚度匹配局部结构大型结构局部高应力区切割面距离足够远对称性利用几何、载荷、约束全对称载荷取半、刚度减半、避开最大应力面4. 材料、单元与属性选型之前先看限制4.1 材料分类各向同性、各向异性与复合材料ANSYS 分析中的材料属性直接决定本构关系选错材料模型比单元选错更致命。各向同性材料在任意一点沿任何方向的力学、热学性能均相同所有金属材料都属此类ANSYS 中只需要定义弹性模量 EX、泊松比 PRXY 和密度 DENS 三个基本参数。各向异性材料的性能方向相关木材顺着纹理和横跨纹理的弹性模量可以差一个数量级合成纤维复合材料更是如此必须在材料定义中给出多组弹性常数。复合材料则是两种及以上材料的组合教案里将其分为单层和多层多层板在 ANSYS 中通常用壳单元的分层截面定义每一层有自己的厚度、材料和铺层角度。材料参数对仿真分析的影响在热力耦合场景里尤其明显如果只是做结构分析线弹性各向同性材料定义三行就够MP,EX,1,2.1E5 ! 弹性模量 210GPa MP,PRXY,1,0.3 ! 泊松比 MP,DENS,1,7.85E-9 ! 密度单位 t/mm^3MP是材料属性定义命令格式为MP,属性,材料编号,数值。注意单位体系的一致性如果在 mm 单位制下弹性模量用 MPa密度用 t/mm³力的单位就是 N。4.2 单元使用限制棱边载荷、面载荷、力矩与扭矩每种单元都有它的不能做的事教案里列出的限制非常实用实体单元不能施加棱边载荷因为实体单元的形函数只支持面力棱边载荷无法形成等效节点力轴对称单元不能施加面载荷只能施加在对应的线段上杆单元没有旋转自由度不能施加节点力矩和扭矩否则会得到零刚度响应。单元限制细节单元类型载荷限制建模提示实体单元不能施加棱边载荷载荷只能加在面或节点上轴对称单元不能施加面载荷载荷施加在对应线段杆单元不能施加节点力矩和扭矩只有轴向自由度无抗弯能力梁单元截面方位、节点偏移需显式控制自由度释放模拟铰接板单元不同厚度连接处节点必须重合网格过渡需要对齐节点板单元这里有个常见失误不同厚度的板连接时如果网格节点不重合计算结果会在连接处出现虚假的应力跳跃。解决办法是在厚度变化的边界处预先设置网格种子点或者用绑结接触把不匹配的网格连接起来。4.3 自由度协调壳单元 ROTZ 的虚拟刚度问题不同单元类型连接时自由度协调是一个隐蔽但致命的坑。杆单元只有轴向平动自由度梁单元有平动加转动自由度壳单元除了面内平动还有一个绕面法线的旋转自由度 ROTZ。问题在于壳单元的 ROTZ 与平面旋转刚度有关而这个刚度是虚拟的不是真实物理刚度所以梁单元和薄壳单元相连时对应的自由度并不协调可能出现局部零刚度求解器直接报奇异。教案里提到 SHELL43 和 SHELL63 单元在KEYOPT(3)2时会激活所有旋转自由度属于例外情况。实际建模里处理梁壳连接我的习惯是优先考虑几何上的合理过渡要么在连接区域使用实体单元做转接要么用 MPC 绑定约束把梁端节点和壳面节点关联起来这两种方式都比强行共用节点可靠。另外两个单元使用细节单元带中间节点时划分网格应让中间节点对齐否则单元边界的位移连续性被破坏相邻单元应当有相同的单元边节点数一边 3 节点一边 2 节点会导致变形不协调。关于非线性分析二次单元的积分点数量不一定比线性单元多但由于位移插值精度更高对接触和大变形问题的收敛性有实质帮助所以在塑性、接触、大变形场景下优先选用二次单元。4.4 单元属性在划分网格之前全部指定单元属性包括单元类型、实常数、材料编号和截面定义。很多新手习惯在生成网格后再回去补材料这样经常出现部分单元使用缺省属性导致结果异常。正确的流程是先定义好所有单元类型、材料、实常数划分网格前通过 Mesh Attributes 给每一类几何实体分配正确的属性。教案里强调了两点在实体模型上直接指定属性时将不考虑缺省属性这可以避免在网格划分操作中重新设置属性而且 ANSYS 的网格划分算法在一次对所有实体进行网格划分时效率更高清除实体模型上的网格不会删除指定的单元属性所以调整网格参数时不需要担心属性丢失。划分完网格后应该激活属性编号校核路径是 Utility Menu PlotCtrls Numbering把单元类型、材料编号、实常数编号都打开目视检查有没有属性错误的区域。壳单元和梁单元的截面定义示例ET,1,SHELL181 ! 壳单元 SECTYPE,1,SHELL ! 定义壳截面 SECDATA,2 ! 壳厚度 2mm ET,2,BEAM189 ! 梁单元 SECTYPE,2,BEAM,RECT ! 矩形梁截面 SECDATA,10,20 ! 截面宽 10mm、高 20mmSECTYPE声明截面类型和编号SECDATA输入截面几何参数。壳单元厚度方向的积分点个数直接影响弯曲应力计算精度积分点太少会低估壳表面的弹塑性响应太多则增加计算代价默认值对大多数薄壳分析是够用的。5. 网格划分的控制原则与无关性验证5.1 网格划分三步与三组选择原则网格划分在 ANSYS 里遵循固定三步定义单元属性、指定网格控制参数、生成网格。对应菜单路径分别是 Preprocessor Meshing Mesh Attributes、Preprocessor Meshing Size Control、以及最终的 Mesh 操作。在 Workbench 里这套流程对应为 Mesh 模块下的 Element Size、Local Sizing 和 Generate Mesh逻辑一致只是换成了参数化界面。网格粗细没有统一答案教案给出了三组判断原则。静力分析中只关心变形时可以画较少的网格因为位移是能量意义上的整体量对局部误差不敏感要计算应力或应变时则需要较多网格因为应力由位移求导得到对误差更敏感。固有属性分析中只算少数低阶模态可以选较少的网格高阶模态的振型复杂需要足够多的单元分辨局部变形。响应分析同理只关注某些部位的位移响应用粗网格即可要算应力响应则必须加密。分析场景粗网格可用必须加密静力分析只计算变形计算应力或应变模态分析只关心低阶模态计算高阶模态响应分析只关心位移响应计算应力响应网格无关性验证的操作流程是选定一组从粗到细的网格尺寸比如 10mm、5mm、2.5mm、1.25mm在每种网格下提取关心点的最大等效应力和最大位移比较相邻两组结果的变化率。步骤如下建立模型并确定目标结果如关键节点的等效应力峰值用最大网格尺寸求解一次提取结果逐次加密网格保持边界条件和载荷完全一致当相邻两次的应力变化小于 3% 时认为结果与网格无关。5.2 网格无关性验证细化到什么程度才合适加密网格到一定程度后结果会进入收敛平台继续细化只会增加求解时间甚至会因为累计误差出现轻微漂移。这时就该停止细化采用当前网格尺寸。每次求解完成后可以用下面的 APDL 命令提取关键结果/POST1 SET,LAST NSORT,S,EQV ! 按等效应力排序 *GET,SMAX,SORT,0,MAX ! 取最大等效应力 NSORT,U,SUM ! 按总位移排序 *GET,UMAX,SORT,0,MAX ! 取最大位移SET,LAST读取最后一个载荷步的结果NSORT对节点结果排序*GET从排序结果中提取极值。把两列数据记录成下面的格式对比网格尺寸 mm最大等效应力 MPa最大位移 mm应力变化率10248.50.532-5265.30.5476.8%2.5270.10.5511.8%1.25271.20.5520.4%到 2.5mm 时应力变化率已经降到 1.8%1.25mm 时降到 0.4%说明 2.5mm 的结果已经进入收敛平台。验证完成后实际生产计算采用 2.5mm 网格即可。如果关心的点恰好位于应力集中区还要顺带检查该点的等效应力是否落在不同单元尺寸下都保持稳定而不是只看全局最大值位置有没有漂移。网格无关性验证的目标不是追求绝对精确而是在计算成本和结果精度之间找到那个拐点。本文还有配套的精品资源点击获取