
简介《ANSYS有限元分析实例》以一块厚20mm、受平面内张力作用的带孔矩形板为典型算例完整演示ANSYS软件进行有限元分析的六步流程定义参数、创建几何模型、划分网格、加载数据、求解与结果分析。每个环节均给出具体的菜单路径、对话框选项和数值输入如工程名与标题设置、单位制选择、PLANE82单元定义、材料弹性模量与泊松比设定、矩形与圆形布尔运算、智能网格划分及后处理查看方式适合机械、土木、航空航天等相关专业初学者对照练习。压缩包内含1个PDF文档大小仅1.04MB内容紧凑、图文并茂可在PC端或移动端随时查阅。目前已有256人学习浏览对于希望快速掌握ANSYS经典界面基本操作流程的读者来说是一份实用的入门辅导材料。1. 带孔板拉伸的 ANSYS 有限元分析从实例 PDF 到能复现的完整套路这份 ANSYS 有限元分析实例 PDF演示的是机械专业课程设计里出现率最高的一个题目一个 200mm×100mm、厚 20mm 的带孔矩形板左边固定右边受 p20N/mm 的均布拉力求变形和应力分布。题目简单但把它完整跑通却要跨过 ANSYS 经典界面里的六道关卡定义参数、创建几何、划分网格、加载、求解、后处理。很多第一次用 ANSYS 的人都是在这份教程上栽跟头的——不是卡在菜单找不到就是算出来的变形方向跟物理直觉对不上。这篇笔记我把每一步的点击路径、参数含义和容易翻车的地方拆开讲新手可以照着复现熟手可以直接跳到第 5 章的避坑清单确认自己的操作习惯。2. 定义参数这一步单元类型、实常数与材料参数为什么必须连在一起设2.1 单元类型选 PLANE828 节点单元和 Plane strs w/thk 的含义打开 ANSYS 经典界面第一步不是画图而是先告诉软件“我要用什么单元”。教程里选的是 Structural 下的 Solid→8node 82也就是 PLANE82。这是 ANSYS 经典的 8 节点四边形平面单元带中间节点位移插值精度比 4 节点的 PLANE42 高一个量级。对于带孔板这种存在应力集中的问题孔边应力梯度大用 8 节点单元能在同样网格密度下拿到更平滑的应力结果。我第一次做这个例子时试过 PLANE42孔边最大应力的数值波动很明显换 PLANE82 后稳定了很多。选完单元类型后还有一个很多人直接点 OK 跳过的选项Element behavior 下拉框。这里必须选 Plane strs w/thk而不是默认的 Plane stress。两者的差别在于——默认平面应力假设板厚方向不受力按单位厚度计算而 Plane strs w/thk 会把板厚作为真实参数参与计算等效于带厚度的平面应力壳。这个题目里厚度 20mm 是明确给定的选了 w/thk 才能在后面对应实常数里输入 THK20。2.2 实常数 THK 与材料参数 EX/PRXY参数间的连带关系定义完单元类型紧接着要做两件事Real Constants实常数和 Material Props材料参数。这两步的顺序有讲究因为 PLANE82 在选了 Plane strs w/thk 之后实常数对话框里才会出现 THK 输入框。如果你在前面没选 w/thk这里的 THK 框可能显示不出来或者输入了也不参与计算。教程里实常数只填一个值THK20单位 mm。材料参数填两个值EX200000PRXY0.3。EX 是弹性模量200000 对应的是钢材常见的 200GPaPRXY 是泊松比0.3 同样是钢材的典型值。这里单位制的逻辑要捋清楚教程在第 1.2 步用/UNIT,SI声明了国际单位制但 ANSYS 内部并不会因为这条命令自动换算单位。你输入的 200000 和 20 到底代表 MPa 还是 Pa取决于你建模时的长度单位是 mm 还是 m。我们按 mm 建模、按 N 施加载荷那么应力单位就是 MPa弹性模量 200000 即 200000MPa完全自洽。参数项输入值含义常见错误THK实常数20板厚单位 mm忘记输入求解报“厚度未定义”EX弹性模量200000200GPa钢材输成 2e11量纲不匹配PRXY泊松比0.3钢材典型值输成 0.33 影响不大但没必要2.3 实际操作路径与命令流对照菜单路径是死的但如果你以后要批量做类似分析可以记下这套操作的 APDL 等价命令在输入窗口敲一遍比点菜单快得多FINISH /CLEAR /PREP7 ET,1,PLANE82 ! 定义 1 号单元为 PLANE82 KEYOPT,1,3,2 ! 设置第 3 个单元选项为 2即 Plane strs w/thk R,1,20 ! 实常数 1 号厚度 20mm MP,EX,1,200000 ! 材料 1 号弹性模量 200000MPa MP,PRXY,1,0.3 ! 材料 1 号泊松比 0.3这段命令里KEYOPT,1,3,2是最容易被忽略的一行。PLANE82 的第 3 个 KEYOPT 控制单元行为2 代表带厚度的平面应力如果不写这一行后面实常数的 THK 就形同虚设。R,1,20里的 1 是实常数编号20 是厚度值和单元类型编号、材料编号互相独立初学者经常把这三个编号搞混其实各管各的只要确保单元定义里引用的实常数编号存在即可。设置完这些前处理的“参数”环节才算真正闭环。3. 创建带孔几何模型矩形、圆孔与布尔减运算的点击顺序3.1 用 By 2Corners 和 Solid Circle 建立基础体素ANSYS 经典界面建模和现在主流的 CAD 软件思路不一样它是基于体素拼接的先画一个矩形再画一个圆然后用布尔减运算把圆从矩形里“抠”出去。几何尺寸按题目来——矩形直接用 By 2Corners 方式在对话框中输入 Width200、Height100起点默认放在全局坐标原点圆心位置教程截图没给具体值常见做法是放在矩形中心即 WP X100、WP Y50圆孔半径取 15mm 左右。如果你在做课程设计圆心和半径按自己题目的要求调整即可方法不变。这里注意一个细节By 2Corners 输入的 Width 和 Height 分别是沿 X 和 Y 方向的尺寸不是对角点坐标。我第一次用的时候以为是两点坐标填了 (0,0) 和 (200,100)出来的矩形位置就不对因为 ANSYS 的 By 2Corners 是用“左下角起点 宽高”来定位的起点默认从 WP工作平面原点开始。3.2 Subtract 布尔减选择顺序决定成败矩形和圆都画好后进入 Modeling→Operate→Booleans→Subtract→Areas。这一步是新手翻车重灾区。弹出的拾取对话框第一次让你选的是“被减的主体”也就是矩形选完矩形后弹出一个对话框让你确认点 OK然后再次弹出拾取框这次选的是“要减去的对象”也就是圆。ANSYS 的拾取逻辑是先选主体、后选客体两次选择之间对话框会切换。教程里特别提到要用 Next 按钮——当你在图形上点击圆时如果拾取窗口里同时显示多个面积需要点 Next 确认当前高亮的是圆而不是矩形外边。如果你反过来先选了圆再选矩形布尔运算会得到一个半月牙形状的残余几何关系直接错乱后面网格和求解全都白做。3.3 建模后的检查习惯布尔运算完成后图形窗口里应该是一个矩形中央带一个圆洞的平面。检查方法很简单Plot→Areas 重新显示面看洞口是否完整再 Plot→Lines 看边线数量矩形四边加上圆孔圆周应该被自动切分成若干段。这个圆孔圆周被分成几段直接关系到后面网格划分时“选圆边”的操作——ANSYS 默认把圆分成 4 条弧线所以后面第 4 章里要用两条圆边分别设置网格尺寸。如果发现布尔运算结果不对最干净的办法不是撤销而是用 Modeling→Delete→Area and Below 把错误的面删掉再重新画矩形和圆。反复撤销操作在经典界面里容易出现面编号混乱我一般直接删了重来比你去找后悔药快得多。4. 划分网格25mm 边尺寸与孔边 2 单元的粗细搭配策略4.1 矩形边 25mm、孔边 2 单元这套参数的来历网格划分是有限元分析里最“玄学”的环节但教程给的参数其实是有讲究的。矩形四边设置 Element edge length25意味着长边 200mm 被分成 8 段短边 100mm 被分成 4 段整个板的基础网格尺寸就是 25mm 见方。而圆孔的四条弧线不按尺寸控制改用 No.of element divisions2即每条弧线只分 2 个单元。为什么孔边要单独加密因为带孔板受拉时应力集中发生在孔边最大应力通常出现在孔的水平直径两端这里的应力梯度非常大。如果孔边网格和远处一样粗算出来的峰值应力会明显偏小而且随着网格加密数值一直波动。孔边加密到每条弧线 2 段整个圆周就是 8 个单元配合 PLANE82 的中间节点能在不大幅增加节点数的前提下捕捉到应力集中的峰值。这套“远处 25、孔边 2 段”的搭配是 ANSYS 教学例题里典型的粗细网格结合方案。4.2 Free 网格划分与孔边应力集中设置完线的尺寸后选择 Meshing→Mesh→Areas→Free拾取全部面点 Pick All。这里用的 Free 网格意味着 ANSYS 自动生成三角形和四边形混合网格对于带孔板这种几何不算复杂的模型完全够用。如果你用 Mapped映射网格虽然单元更规整但圆形边界会让映射控制变得麻烦需要额外设置边的对应关系教学案例没必要折腾。网格划分完成后马上做一件事Plot→Nodes 显示节点肉眼看孔边节点分布是否均匀。Free 网格在圆孔附近会自动加密这是 ANSYS 的智能尺寸传递在起作用。如果发现某个区域网格特别乱可以返工Meshing→Clear→Areas 先清除面网格调整边的 division 后再重新划分。4.3 网格质量检查的常见做法正规做法是划分后查看网格质量但教学案例中只要做到两点就基本达标第一孔边最小单元尺寸不能比厚度小太多否则单元长宽比畸形第二整个模型没有出现裂缝或重叠单元。可视化检查可以用 Utility Menu→PlotCtrls→Style→Size Shape打开单元形状显示再 Plot Elements 看单元拉伸情况。理论上网格数量越多精度越高本题目 Total Nodes 一般在一千级求解时间以秒计。如果你算出来的结果明显不合理先别怀疑算法回去看网格——尤其要看圆孔附近的网格有没有在布尔运算时被牵连产生畸形单元。这是 ANSYS 初学者最容易忽略的一点几何没问题、加载没问题但网格里埋了雷后处理数值怎么都对不上。5. 加载与求解避坑约束边选错、压力正负号和量纲不清是最大的三个坑5.1 施加位移约束与面压力加载阶段的操作路径本身不复杂但两个“选边”动作必须做对。先在 Preprocessor→Loads→Define Loads→Apply→Structural→Displacement→On Lines 中拾取矩形的左边线也就是固定端约束类型选 ALL DOF位移值填 0。这一步的物理含义是把左端“焊死”六个自由度全锁。然后施加压力载荷Apply→Structural→Pressure→On Lines这次拾取右边线Load PRES value 一项是关键。教程里填的是 -1但题设明明写的是 p20N/mm为什么是 -1这里需要说明两件事第一ANSYS 压力载荷的默认方向是“压向面内”正值为压缩、负值为拉伸右端受向右的拉力是拉应力所以一定要填负号第二数值上填 -1 还是 -20取决于你想让结果应力读多少。按题设 p20N/mm、板厚 20mm均布拉应力是 1MPa如果填 -20N/mm得到的应力是 1MPa 的 20 倍20MPa 级别教程填 -1 更多是为了演示流程、让云图数值好读真正做课程设计时建议直接填 -20让载荷的量纲和题设严格一致。5.2 压力正负号与量纲两个最容易翻车的点压力符号问题值得单独拎出来说。很多人第一次做这个题填了 1 而不是 -1结果求解出来模型的位移方向变成向左压缩和右端拉伸的物理图景完全相反。更隐蔽的问题是量纲如果你的 EX 输入的是 2e11Pa而几何尺寸用的还是 mm那么载荷单位也必须跟着换成 N/m 量级的等效值否则应力结果会差 6 个数量级。ANSYS 经典界面不会主动提醒你单位不一致这是它和 Workbench 最大的区别之一——Workbench 会做单位换算经典界面完全靠人工约定。教程第 1.2 步输入/UNIT,SI只是显示单位制不影响内部数值计算很多教程没强调这一点导致不少人以为输完这条命令单位就自动统一了。真实的单位制是靠“几何长度单位 材料参数单位 载荷单位”三者自洽来保证的。我通常的习惯是全部按 mm、N、MPa 体系来EX 输 200000 而不是 2e11压力输 -20N/mm这样后处理出来的位移单位是 mm应力单位是 MPa无需任何换算。5.3 求解前的检查清单点击 Solution→Solve→Current LS 之前建议花一分钟过一遍这个清单检查项正确状态错误后果左边线 ALL DOF0图上约束符号在左边缘模型刚体位移求解报错右边线压力为负值载荷箭头指向外变形方向反向实常数 THK20 已设单元选项为 w/thk厚度不参与计算材料 EX/PRXY 已设材料号正确刚度为零无法求解5.4 常见问题排查现象 1求解报错“Degree of freedom solution does not exist”或刚度矩阵奇异。原因约束没加或只加了部分自由度模型存在刚体位移。解决回到位移约束步骤确认左边线被完整选中ALL DOF 置 0。现象 2变形图显示整个板向左边飞出去。原因压力正负号填反了ANSYS 默认正压力指向面内你填的是正值所以模型被压缩。解决把 Load PRES value 改为负值。现象 3应力云图数值大得离谱。原因弹性模量单位输错了用的是 Pa 体系2e11但几何尺寸是 mm导致应力被放大。解决统一单位体系EX 填 200000载荷按 N/mm 输入。现象 4孔边应力看起来很均匀没有明显的应力集中。原因孔边网格太粗或者布尔减运算后圆孔边缘线没有单独布网格尺寸。解决重新给孔的四条弧线设置 No.of element divisions2并重新划分网格。现象 5选择边线时总是选不中想要的边。原因拾取对话框默认是单选模式图形窗口里有大量线重叠。解决把拾取模式从 Single 改为 Polygon或先放大视图再点选选多条线时按住鼠标左键拖出矩形框选。6. 后处理验证变形图、Von Mises 应力与孔边应力集中系数6.1 读入结果与显示变形求解完成后进入 General PostProc通用后处理器第一步是 Read Results→First Set把结果文件里的第一个载荷步读进来。如果不做这一步后处理菜单里很多选项是灰色的点不动。然后 Plot Results→Deformed Shape选择 Def undef edge这样变形后的形状和原始边界会叠加显示能直观看出位移方向是向右拉伸的左端不动、右端位移最大和物理预期一致。6.2 Von Mises 应力的正确打开方式应力云图用 Plot Results→Contour Plot→Nodal Solu 来做在列表里找到 Von Mises stress。这里有个细节理想情况下应选择 Nodal Solution 下的结果而不是 Element Solution。节点解是经过平均后的连续应力场云图平滑单元解在每个单元边界会出现不连续的台阶新手看到那种锯齿状云图会以为网格有问题其实是选错了结果类型。另外把 Interpolation Nodes 设置为 All applicable云图显示效果会更细腻这也是教程里特意提到的一个隐藏设置。6.3 和理论解对照孔边应力集中系数如果你按网格加密后的结果读取孔边最大 Von Mises 应力会发现数值明显大于远离孔区域的均匀应力。这正好验证了力学概念里的孔边应力集中。薄板圆孔受拉的应力集中系数理论值约 3但那是无限宽板本例题板宽 100mm、孔径约 30mm属于有限宽板基于 Howland 解的修正公式是 K_t 3.0 − 3.13(d/W) 3.66(d/W)² − 1.53(d/W)³代入 d/W0.3 后约为 2.2。网格粗细不同数值会在 1.82.3 之间波动能用这个范围判断你的网格是否足够。从那以后我每次做带孔板的 ANSYS 分析都会强制走一遍“理论系数估算 → 数值结果对照”的流程数值落入修正公式附近才敢说结果是可信的。希望这个习惯也能帮到你。6.4 顺手能用的两个小技巧后处理里还有两个高频操作值得记一下。一是用 Query Picker查询拾取器点击孔边节点直接读出该节点的应力值比云图上读色标准得多适用对象是带孔板、切口板这类找峰值应力的场景。二是把当前视角调整到孔边局部用 PlotCtrls→Pan Zoom Rotate 放大配合云图能清楚看到最大应力点是否出现在孔的水平直径两端——如果偏离了通常说明网格不对称或边界条件施加有问题。这个例子跑通后同样的流程可以直接迁移到 ANSYS Workbench 的静力分析模块甚至加一步模态分析来验证板的固有频率底层逻辑完全一致。本文还有配套的精品资源点击获取