第一次接触有限元仿真的研究生十有八九会遇到这种场景导师丢过来一句“用有限元算一下这个结构的应力”你打开软件导入模型画网格点求解几分钟后屏幕上出现一张五彩斑斓的云图。那一刻感觉良好直到导师追问一句“你用的哪种单元约束怎么加的网格收敛了没有材料实验数据是哪来的”——你才发现自己什么都不确定。这不是个例。在带过的师弟师妹里会点软件按钮的人很多能在“算完之后”把每一步选择的依据讲清楚的人很少。有限元仿真的难点从来不在操作而在建模判断模型该简化成平面问题还是实体问题材料该用各向同性还是正交各向异性边界条件怎么加才不会人为改变受力路径网格加密到什么程度算收敛结果和理论估算差多少倍以内可以接受。这些问题软件不会替你想。作为研二师兄这篇文章想给的是一套能被复用的思路不是某个软件的按钮攻略。核心判断是有限元的本质是“带误差的计算”你做每一步都在引入误差只有把误差控住结果才敢用。全文会按这个顺序展开先画清有限元仿真的整体思维链路再给出软件选型和基础学习建议然后用一个碳纤维复合材料板拉伸仿真实例把“理论估算→材料常数→单元网格→边界条件→求解→验证”完整走一遍最后整理新手最常见的报错和工程实践建议。如果你正准备入门或者已经卡在“算出来的结果不敢用”的阶段这篇文章值得收藏。1. 这篇文章真正要解决的问题先说清楚这不是一篇软件操作手册。ANSYS、Abaqus、COMSOL 的官方教程和帮助文档比任何博客都详细照着点一遍谁都会。真正挡在初学者面前的是下面几类问题你可以先对号入座已经会导出云图但导师一追问“算得对不对”就答不上来知道有限元三个字怎么写但不知道建模时哪些东西能简化、哪些不能面对碳纤维复合材料这类各向异性材料不知道工程常数到底该怎么填、方向怎么对网格越画越密算得越来越慢仍然不知道收敛点在哪里模型报错或者结果明显不合理时没有一套系统的排查思路。这篇文章的判断标准很简单读完之后你能独立说清楚“我这个模型做了哪些简化、引入了哪些误差、结果在多大范围内可信”比你会用多少个软件都重要。仿真做的不是“算得漂亮”而是“算得可信”。后面所有内容都围绕“可信”两个字展开。2. 有限元的底层逻辑节点、单元、刚度矩阵初学者最容易犯的错是绕过原理直接学操作。结果就是换一个软件就不会用了。有限元方法再复杂底层也就是把一个连续体切成有限个小块在每个小块上做近似再组装成一个大方程[ [K]{u}{F} ]这三个符号对应了你建模时做的三件事([K]) 是刚度矩阵——由材料本构、单元类型、几何尺寸共同决定。你填的材料参数错了它一定错({F}) 是载荷向量——你施加的力、压力、温度、位移都会到这里({u}) 是节点位移——这是软件最后求出来的东西应力应变是在位移基础上推导的。很多人以为“我点一下 Solve软件就在算有限元”。其实你花大量时间做的建模本质上都是在组装这个方程。有一个判断值得早一点建立如果材料参数错了或者约束边界条件错了任何求解器都救不了你。软件能做的是解好你给的方程不是判断你给的方程是否符合物理。下表是几个典型错误对结果的污染程度错误类型对结果的影响软件会不会报错材料弹性模量填错数量级位移和应力全面偏离一般不会各向异性材料方向没对齐应力分布完全错误一般不会约束不足导致刚体位移求解失败或位移离谱可能报奇异错误载荷加错位置局部应力路径错误不会单元类型选错平面应力/平面应变薄板/厚板结果完全不同不会这五类错误里至少三类不会触发任何警告。这也是为什么“会操作”和“会仿真”是两回事。2.1 再往前一步单元、节点与自由度单元是有限元模型的基本组成块节点是单元的角点或边上的关键点。每个节点有几个“自由度”三维结构问题通常是三个平动自由度梁/壳单元还会带转动自由度。求解器最后求的是每个节点的自由度值单元内部的值靠形函数插值得到。新手要建立的一个基础直觉是单元越小插值误差越小但计算量越大。这个矛盾贯穿整个仿真过程后面第 5 节的网格收敛性检查就是专门处理它的。3. 软件选型与基础学习路径先判断工具边界经常收到类似问题“研二了学 ANSYS 还是 Abaqus还是学 COMSOLMATLAB 里的 Simscape 能算有限元吗”从实践角度看学习路径的正确顺序是“先定问题再选工具”而不是反过来。结构静力学、强度校核、模态分析ANSYS Mechanical 和 Abaqus 都很成熟接触、材料非线性、大变形问题Abaqus 的求解器经验更丰富多物理场耦合比如热-结构、流-固耦合COMSOL 的建模方式更顺手课程作业和简单验证开源的 CalculiX、FreeCAD FEM 也够用还能省去许可证问题至于 MATLAB 里的 Simscape它属于系统级建模工具把电机、液压、机械传动等部件搭成微分方程模型适合做“子系统控制逻辑”的仿真不适合回答“这个区域应力是多少”这类连续体力学的细节问题。换句话说Simscape 和你说的有限元仿真解决的不是同一层问题。研究一个板的应力分布用结构有限元软件研究一个机械臂的关节驱动和控制才考虑 Simscape。3.1 一条可执行的学习路线不要一上来就啃有限元理论大部头也不要一上来就做课题组最难的项目。比较稳妥的路线是材料力学、弹性力学基础概念先过关应力、应变、弹性模量、泊松比、平面应力与平面应变的区别找 3 个有解析解的经典问题练手悬臂梁弯曲、开孔板拉伸、复合材料层合板拉伸分别用软件算一遍和解析解对比每个问题都做一次网格收敛性检查体会“结果随网格加密而变化”的过程再进入课题组项目从别人做过的模型开始改而不是从零建模尽早学脚本化操作APDL、Python、Abaqus Python API让模型可复现、可追溯。这条路线里最关键的是第 2 步和第 3 步。很多新手直接跳到第 4 步结果模型一错就不知道错在哪因为没有“正确基准”可以对照。4. 完整实例碳纤维复合材料板拉伸仿真下面用一个非常典型的问题把完整流程走一遍。问题描述一块碳纤维增强树脂基复合材料层合板长 100 mm、宽 25 mm、总厚度 1 mm。一端固定另一端施加 10000 N 的拉伸力。要求计算板内的应力和伸长量。这个题目看起来简单但足够覆盖建模决策的所有关键环节。我们逐步处理。4.1 第一步先算理论解给仿真设一个“锚点”很多人的习惯是拿到问题直接建模。更推荐的做法是先花五分钟做理论估算哪怕只算一个量纲正确的近似值。这个值就是后面判断仿真结果是否可信的“锚点”。以单向 0° 铺层纤维沿拉伸方向为例最简单的情况可以手算截面积 (A 25 \times 1 25 \text{ mm}^2)名义应力 (\sigma F/A 10000 / 25 400 \text{ MPa})若材料纵向模量 (E_1 135 \text{ GPa})伸长量 (\Delta L FL/(EA) 10000 \times 100 / (135000 \times 25) \approx 0.296 \text{ mm})如果铺层换成 ([0/90]_s) 这样的正交铺层等效模量就不是 135 GPa 了需要用层合板理论估算。这里给一个用 Python 做预估算的最小脚本它能在你打开仿真软件之前先在命令行里得到一个可信的预期值。# 文件clt_quick_check.py # 作用仿真前用经典层合板理论快速估算等效模量 # 运行环境Python 3.8numpy import numpy as np # 单向碳纤维复合材料工程常数示例值真实数据以材料试验报告为准 E1 135e9 # 纤维方向拉伸模量, Pa E2 8.8e9 # 横向拉伸模量, Pa G12 4.47e9 # 面内剪切模量, Pa nu12 0.30 # 主泊松比 # 柔量矩阵分量 S11 1.0 / E1 S22 1.0 / E2 S12 -nu12 / E1 S66 1.0 / G12 # 二维刚度矩阵 Q denom S11 * S22 - S12**2 Q11 S22 / denom Q22 S11 / denom Q12 -S12 / denom Q66 1.0 / S66 # [0/90]s 对称层合板0° 与 90° 层等厚 A11_h 0.5 * (Q11 Q22) A12_h Q12 A22_h A11_h # 面内等效拉伸模量 Ex_lam A11_h - A12_h**2 / A22_h print(f单层板纵向模量 E1 {E1/1e9:8.2f} GPa) print(f单层板横向模量 E2 {E2/1e9:8.2f} GPa) print(f[0/90]s 面内等效模量 Ex {Ex_lam/1e9:8.2f} GPa) # 用等效模量估算 100mm×25mm×1mm 板在 10000N 下的变形 F 10000.0 # N L 100e-3 # m b 25e-3 # m t 1e-3 # m A b * t sigma F / A delta F * L / (Ex_lam * A) print(f宏观应力 sigma {sigma/1e6:8.2f} MPa) print(f宏观伸长量 delta_L {delta*1e3:8.4f} mm)运行后预期输出大致为单层板纵向模量 E1 135.00 GPa 单层板横向模量 E2 8.80 GPa [0/90]s 面内等效模量 Ex 72.32 GPa 宏观应力 sigma 400.00 MPa 宏观伸长量 delta_L 0.5531 mm这里真正要注意的问题有两个。第一复合材料层合板的“等效模量”不等于单层板的纵向模量它取决于铺层角度和顺序第二宏观应力是 400 MPa不代表每一层内部应力都是 400 MPa0° 层承担得多90° 层承担得少要看到层内应力必须用分层单元建模。这就是“理论锚点”的价值——它告诉你整体量级也告诉你模型的局限在哪。4.2 第二步材料工程常数怎么填碳纤维复合材料单层板属于正交各向异性材料三维情况下需要 9 个独立工程常数常数含义典型量级示例(E_1)纤维方向拉伸模量135 GPa(E_2)、(E_3)横向拉伸模量8~10 GPa(G_{12})、(G_{13})面内剪切模量4~5 GPa(G_{23})横观剪切模量3~4 GPa(\nu_{12})、(\nu_{13})主泊松比0.25~0.35(\nu_{23})横向泊松比0.35~0.45新手最容易踩的坑不是不知道填这些数而是不知道这些数必须来自材料试验报告或可靠的文献数据并且要特别注意坐标系对齐。碳纤维复合材料只有在纤维方向上的模量才是 (E_1)如果你的模型里纤维方向是 45°那材料坐标系就必须旋转 45°。在二维平面应力问题里只需要关注面内的 4 个常数(E_1)、(E_2)、(\nu_{12})、(G_{12})。但一旦升级为三维实体模型就必须补齐全部 9 个。这也是为什么建议新手从二维问题开始——因为变量少容易验证。4.3 第三步ANSYS APDL 完整示例下面用 ANSYS 经典环境APDL写一个最小批处理示例。为了方便和理论值对照这里直接用单向 (0^\circ) 铺层即材料纵向模量 (E_1135) GPa 沿 X 方向。单位制采用 mm、N、MPa三者必须配套。! 文件composite_plate_2d.inpAPDL 批处理示例 ! 单位制mm、N、MPa必须全程一致 ! 题目100mm×25mm 单向板左端固定右端受 10000N 拉伸 /PREP7 ! 四节点平面应力单元KEYOPT(3)2 表示平面应力 ET,1,PLANE182 KEYOPT,1,3,2 ! 正交各向异性材料示例工程常数MPa MP,EX,1,135000 MP,EY,1,8800 MP,PRXY,1,0.30 MP,GXY,1,4470 ! 几何100 mm × 25 mm 矩形板 RECTNG,0,100,0,25 ! 网格单元尺寸 2 mm ESIZE,2 AMESH,ALL ! 左端固定X0 处节点施加全部约束 NSEL,S,LOC,X,0 D,ALL,ALL,0 ALLSEL ! 右端施加节点力总拉力 10000 N方向 X NSEL,S,LOC,X,100 *GET,NCNT,NODE,,COUNT F,ALL,FX,10000/NCNT ALLSEL ! 求解 /SOLU SOLVE FINISH ! 后处理 /POST1 PLNSOL,U,X ! X 方向位移云图 PLNSOL,S,X ! X 方向正应力云图 PRRSOL ! 支反力汇总这段代码里有几个值得注意的建模决策平面应力单元板厚 1 mm相对面内尺寸很小属于薄板问题适合平面应力假设。如果你把它当成厚块就得用实体单元材料常数也得补全到 9 个载荷施加方式把 10000 N 均分到右端每个节点上这是初学者最容易引起局部应力集中的操作。按圣维南原理端部附近的应力会受加载方式影响但只要看远离两端的中部区域结果就接近理论均匀应力 400 MPa支反力检查PRRSOL会输出固定端支反力。理想情况下支反力在 X 方向的合力应当等于 -10000 N这是最基本的力平衡验证。如果要在 ANSYS Workbench 里完成同样的问题思路完全一致在 Engineering Data 里填正交各向异性材料几何模型用一块 100×25 mm 的平面板一端固定另一端施加 10000 N 的力然后在结果里读取应力和位移。APDL 的好处是过程可脚本化、可追溯更适合课题组做参数化研究。4.4 第四步求解 有限元仿真入门：建模判断与误差控制全指南
第一次接触有限元仿真的研究生十有八九会遇到这种场景导师丢过来一句“用有限元算一下这个结构的应力”你打开软件导入模型画网格点求解几分钟后屏幕上出现一张五彩斑斓的云图。那一刻感觉良好直到导师追问一句“你用的哪种单元约束怎么加的网格收敛了没有材料实验数据是哪来的”——你才发现自己什么都不确定。这不是个例。在带过的师弟师妹里会点软件按钮的人很多能在“算完之后”把每一步选择的依据讲清楚的人很少。有限元仿真的难点从来不在操作而在建模判断模型该简化成平面问题还是实体问题材料该用各向同性还是正交各向异性边界条件怎么加才不会人为改变受力路径网格加密到什么程度算收敛结果和理论估算差多少倍以内可以接受。这些问题软件不会替你想。作为研二师兄这篇文章想给的是一套能被复用的思路不是某个软件的按钮攻略。核心判断是有限元的本质是“带误差的计算”你做每一步都在引入误差只有把误差控住结果才敢用。全文会按这个顺序展开先画清有限元仿真的整体思维链路再给出软件选型和基础学习建议然后用一个碳纤维复合材料板拉伸仿真实例把“理论估算→材料常数→单元网格→边界条件→求解→验证”完整走一遍最后整理新手最常见的报错和工程实践建议。如果你正准备入门或者已经卡在“算出来的结果不敢用”的阶段这篇文章值得收藏。1. 这篇文章真正要解决的问题先说清楚这不是一篇软件操作手册。ANSYS、Abaqus、COMSOL 的官方教程和帮助文档比任何博客都详细照着点一遍谁都会。真正挡在初学者面前的是下面几类问题你可以先对号入座已经会导出云图但导师一追问“算得对不对”就答不上来知道有限元三个字怎么写但不知道建模时哪些东西能简化、哪些不能面对碳纤维复合材料这类各向异性材料不知道工程常数到底该怎么填、方向怎么对网格越画越密算得越来越慢仍然不知道收敛点在哪里模型报错或者结果明显不合理时没有一套系统的排查思路。这篇文章的判断标准很简单读完之后你能独立说清楚“我这个模型做了哪些简化、引入了哪些误差、结果在多大范围内可信”比你会用多少个软件都重要。仿真做的不是“算得漂亮”而是“算得可信”。后面所有内容都围绕“可信”两个字展开。2. 有限元的底层逻辑节点、单元、刚度矩阵初学者最容易犯的错是绕过原理直接学操作。结果就是换一个软件就不会用了。有限元方法再复杂底层也就是把一个连续体切成有限个小块在每个小块上做近似再组装成一个大方程[ [K]{u}{F} ]这三个符号对应了你建模时做的三件事([K]) 是刚度矩阵——由材料本构、单元类型、几何尺寸共同决定。你填的材料参数错了它一定错({F}) 是载荷向量——你施加的力、压力、温度、位移都会到这里({u}) 是节点位移——这是软件最后求出来的东西应力应变是在位移基础上推导的。很多人以为“我点一下 Solve软件就在算有限元”。其实你花大量时间做的建模本质上都是在组装这个方程。有一个判断值得早一点建立如果材料参数错了或者约束边界条件错了任何求解器都救不了你。软件能做的是解好你给的方程不是判断你给的方程是否符合物理。下表是几个典型错误对结果的污染程度错误类型对结果的影响软件会不会报错材料弹性模量填错数量级位移和应力全面偏离一般不会各向异性材料方向没对齐应力分布完全错误一般不会约束不足导致刚体位移求解失败或位移离谱可能报奇异错误载荷加错位置局部应力路径错误不会单元类型选错平面应力/平面应变薄板/厚板结果完全不同不会这五类错误里至少三类不会触发任何警告。这也是为什么“会操作”和“会仿真”是两回事。2.1 再往前一步单元、节点与自由度单元是有限元模型的基本组成块节点是单元的角点或边上的关键点。每个节点有几个“自由度”三维结构问题通常是三个平动自由度梁/壳单元还会带转动自由度。求解器最后求的是每个节点的自由度值单元内部的值靠形函数插值得到。新手要建立的一个基础直觉是单元越小插值误差越小但计算量越大。这个矛盾贯穿整个仿真过程后面第 5 节的网格收敛性检查就是专门处理它的。3. 软件选型与基础学习路径先判断工具边界经常收到类似问题“研二了学 ANSYS 还是 Abaqus还是学 COMSOLMATLAB 里的 Simscape 能算有限元吗”从实践角度看学习路径的正确顺序是“先定问题再选工具”而不是反过来。结构静力学、强度校核、模态分析ANSYS Mechanical 和 Abaqus 都很成熟接触、材料非线性、大变形问题Abaqus 的求解器经验更丰富多物理场耦合比如热-结构、流-固耦合COMSOL 的建模方式更顺手课程作业和简单验证开源的 CalculiX、FreeCAD FEM 也够用还能省去许可证问题至于 MATLAB 里的 Simscape它属于系统级建模工具把电机、液压、机械传动等部件搭成微分方程模型适合做“子系统控制逻辑”的仿真不适合回答“这个区域应力是多少”这类连续体力学的细节问题。换句话说Simscape 和你说的有限元仿真解决的不是同一层问题。研究一个板的应力分布用结构有限元软件研究一个机械臂的关节驱动和控制才考虑 Simscape。3.1 一条可执行的学习路线不要一上来就啃有限元理论大部头也不要一上来就做课题组最难的项目。比较稳妥的路线是材料力学、弹性力学基础概念先过关应力、应变、弹性模量、泊松比、平面应力与平面应变的区别找 3 个有解析解的经典问题练手悬臂梁弯曲、开孔板拉伸、复合材料层合板拉伸分别用软件算一遍和解析解对比每个问题都做一次网格收敛性检查体会“结果随网格加密而变化”的过程再进入课题组项目从别人做过的模型开始改而不是从零建模尽早学脚本化操作APDL、Python、Abaqus Python API让模型可复现、可追溯。这条路线里最关键的是第 2 步和第 3 步。很多新手直接跳到第 4 步结果模型一错就不知道错在哪因为没有“正确基准”可以对照。4. 完整实例碳纤维复合材料板拉伸仿真下面用一个非常典型的问题把完整流程走一遍。问题描述一块碳纤维增强树脂基复合材料层合板长 100 mm、宽 25 mm、总厚度 1 mm。一端固定另一端施加 10000 N 的拉伸力。要求计算板内的应力和伸长量。这个题目看起来简单但足够覆盖建模决策的所有关键环节。我们逐步处理。4.1 第一步先算理论解给仿真设一个“锚点”很多人的习惯是拿到问题直接建模。更推荐的做法是先花五分钟做理论估算哪怕只算一个量纲正确的近似值。这个值就是后面判断仿真结果是否可信的“锚点”。以单向 0° 铺层纤维沿拉伸方向为例最简单的情况可以手算截面积 (A 25 \times 1 25 \text{ mm}^2)名义应力 (\sigma F/A 10000 / 25 400 \text{ MPa})若材料纵向模量 (E_1 135 \text{ GPa})伸长量 (\Delta L FL/(EA) 10000 \times 100 / (135000 \times 25) \approx 0.296 \text{ mm})如果铺层换成 ([0/90]_s) 这样的正交铺层等效模量就不是 135 GPa 了需要用层合板理论估算。这里给一个用 Python 做预估算的最小脚本它能在你打开仿真软件之前先在命令行里得到一个可信的预期值。# 文件clt_quick_check.py # 作用仿真前用经典层合板理论快速估算等效模量 # 运行环境Python 3.8numpy import numpy as np # 单向碳纤维复合材料工程常数示例值真实数据以材料试验报告为准 E1 135e9 # 纤维方向拉伸模量, Pa E2 8.8e9 # 横向拉伸模量, Pa G12 4.47e9 # 面内剪切模量, Pa nu12 0.30 # 主泊松比 # 柔量矩阵分量 S11 1.0 / E1 S22 1.0 / E2 S12 -nu12 / E1 S66 1.0 / G12 # 二维刚度矩阵 Q denom S11 * S22 - S12**2 Q11 S22 / denom Q22 S11 / denom Q12 -S12 / denom Q66 1.0 / S66 # [0/90]s 对称层合板0° 与 90° 层等厚 A11_h 0.5 * (Q11 Q22) A12_h Q12 A22_h A11_h # 面内等效拉伸模量 Ex_lam A11_h - A12_h**2 / A22_h print(f单层板纵向模量 E1 {E1/1e9:8.2f} GPa) print(f单层板横向模量 E2 {E2/1e9:8.2f} GPa) print(f[0/90]s 面内等效模量 Ex {Ex_lam/1e9:8.2f} GPa) # 用等效模量估算 100mm×25mm×1mm 板在 10000N 下的变形 F 10000.0 # N L 100e-3 # m b 25e-3 # m t 1e-3 # m A b * t sigma F / A delta F * L / (Ex_lam * A) print(f宏观应力 sigma {sigma/1e6:8.2f} MPa) print(f宏观伸长量 delta_L {delta*1e3:8.4f} mm)运行后预期输出大致为单层板纵向模量 E1 135.00 GPa 单层板横向模量 E2 8.80 GPa [0/90]s 面内等效模量 Ex 72.32 GPa 宏观应力 sigma 400.00 MPa 宏观伸长量 delta_L 0.5531 mm这里真正要注意的问题有两个。第一复合材料层合板的“等效模量”不等于单层板的纵向模量它取决于铺层角度和顺序第二宏观应力是 400 MPa不代表每一层内部应力都是 400 MPa0° 层承担得多90° 层承担得少要看到层内应力必须用分层单元建模。这就是“理论锚点”的价值——它告诉你整体量级也告诉你模型的局限在哪。4.2 第二步材料工程常数怎么填碳纤维复合材料单层板属于正交各向异性材料三维情况下需要 9 个独立工程常数常数含义典型量级示例(E_1)纤维方向拉伸模量135 GPa(E_2)、(E_3)横向拉伸模量8~10 GPa(G_{12})、(G_{13})面内剪切模量4~5 GPa(G_{23})横观剪切模量3~4 GPa(\nu_{12})、(\nu_{13})主泊松比0.25~0.35(\nu_{23})横向泊松比0.35~0.45新手最容易踩的坑不是不知道填这些数而是不知道这些数必须来自材料试验报告或可靠的文献数据并且要特别注意坐标系对齐。碳纤维复合材料只有在纤维方向上的模量才是 (E_1)如果你的模型里纤维方向是 45°那材料坐标系就必须旋转 45°。在二维平面应力问题里只需要关注面内的 4 个常数(E_1)、(E_2)、(\nu_{12})、(G_{12})。但一旦升级为三维实体模型就必须补齐全部 9 个。这也是为什么建议新手从二维问题开始——因为变量少容易验证。4.3 第三步ANSYS APDL 完整示例下面用 ANSYS 经典环境APDL写一个最小批处理示例。为了方便和理论值对照这里直接用单向 (0^\circ) 铺层即材料纵向模量 (E_1135) GPa 沿 X 方向。单位制采用 mm、N、MPa三者必须配套。! 文件composite_plate_2d.inpAPDL 批处理示例 ! 单位制mm、N、MPa必须全程一致 ! 题目100mm×25mm 单向板左端固定右端受 10000N 拉伸 /PREP7 ! 四节点平面应力单元KEYOPT(3)2 表示平面应力 ET,1,PLANE182 KEYOPT,1,3,2 ! 正交各向异性材料示例工程常数MPa MP,EX,1,135000 MP,EY,1,8800 MP,PRXY,1,0.30 MP,GXY,1,4470 ! 几何100 mm × 25 mm 矩形板 RECTNG,0,100,0,25 ! 网格单元尺寸 2 mm ESIZE,2 AMESH,ALL ! 左端固定X0 处节点施加全部约束 NSEL,S,LOC,X,0 D,ALL,ALL,0 ALLSEL ! 右端施加节点力总拉力 10000 N方向 X NSEL,S,LOC,X,100 *GET,NCNT,NODE,,COUNT F,ALL,FX,10000/NCNT ALLSEL ! 求解 /SOLU SOLVE FINISH ! 后处理 /POST1 PLNSOL,U,X ! X 方向位移云图 PLNSOL,S,X ! X 方向正应力云图 PRRSOL ! 支反力汇总这段代码里有几个值得注意的建模决策平面应力单元板厚 1 mm相对面内尺寸很小属于薄板问题适合平面应力假设。如果你把它当成厚块就得用实体单元材料常数也得补全到 9 个载荷施加方式把 10000 N 均分到右端每个节点上这是初学者最容易引起局部应力集中的操作。按圣维南原理端部附近的应力会受加载方式影响但只要看远离两端的中部区域结果就接近理论均匀应力 400 MPa支反力检查PRRSOL会输出固定端支反力。理想情况下支反力在 X 方向的合力应当等于 -10000 N这是最基本的力平衡验证。如果要在 ANSYS Workbench 里完成同样的问题思路完全一致在 Engineering Data 里填正交各向异性材料几何模型用一块 100×25 mm 的平面板一端固定另一端施加 10000 N 的力然后在结果里读取应力和位移。APDL 的好处是过程可脚本化、可追溯更适合课题组做参数化研究。4.4 第四步求解