力学仿真ansys

我做过一个悬臂梁弯曲的 ANSYS 仿真项目，客户想对比不同端部载荷下梁的应力分布和挠曲线，验证有限元结果和欧拉-伯努利梁理论是否一致。这个案例非常适合作为固体力学仿真的入门验证。

左侧图是悬臂梁在端部集中力 F 作用下的弯曲应力云图。梁左端固定，右端受向上（或向下）的集中力。颜色从蓝到红代表弯曲正应力 σ_x 从负到正，上下表面应力最大，中性轴处应力为零。图中应力范围从 −100 MPa 到 +100 MPa，说明最大弯曲应力出现在固定端上下表面。

右侧图是挠曲线 v(x)，也就是梁的竖向位移随位置变化。三条曲线分别对应 F=1 N、2 N、3 N。可以看到挠度和载荷成正比，F=3 N 时自由端位移最大，接近 −2.0 mm（ normalized 坐标）。图中还标注了欧拉-伯努利梁的解析解：v = Fx²(3L−x)/(6EI)，其中 L 是梁长，E 是弹性模量，I 是截面惯性矩。有限元结果和这条理论曲线重合得很好。

建模设置：梁长 6 m，截面设为矩形或圆形，材料弹性模量 200 GPa，泊松比 0.3。左端固定约束所有自由度，右端施加集中力。单元可以用 BEAM188，也可以用实体单元 SOLID185 再做对比。BEAM188 一维梁单元计算快，适合细长梁；SOLID185 能看局部应力集中，但网格要更细。

踩坑经验：第一个是单元选择。如果梁很细长（L/h > 10），BEAM188 足够；如果梁比较短粗，剪切变形不可忽略，要用 Timoshenko 梁单元或实体单元。第二个是载荷方向：集中力必须垂直于梁轴，如果施加成轴向力，结果会变成拉压而不是弯曲。第三个是网格密度：用实体单元时沿长度方向至少 20 个单元，沿高度方向至少 4 层，否则挠曲线会不光滑。

工程经验：悬臂梁是验证有限元设置正确性的标准模型。我每次换 ANSYS 版本或新电脑，都会先跑一遍悬臂梁，确认结果和解析解一致后再做复杂模型。另外，弯曲应力 σ = My/I 和挠度 v = FL³/(3EI) 这两个公式，客户现场问的时候可以直接用，非常直观。

更多 ANSYS 静力学、梁结构和材料力学仿真的实战案例，我整理在 [力学仿真ansys](https://www.keyanxueshu.com/category/fea/ansys/) 栏目里。[力学仿真ansys](https://www.keyanxueshu.com/)

悬臂梁仿真还有一个重要延伸：判断什么时候必须用 Timoshenko 梁理论。细长比 L/h>10 时，Euler-Bernoulli 理论足够；L/h 在 5-10 之间，剪切变形开始贡献 5%-15% 的挠度；L/h<5 时必须用 Timoshenko 或实体单元。我们这个案例梁长 6 m，如果截面高度 0.3 m，L/h=20，Euler-Bernoulli 完全适用。右侧三条挠曲线和 F 成正比，也验证了小变形线弹性假设。如果客户研究的是大挠度（δ>L/20），就要打开几何非线性，挠曲线会偏离三次抛物线，呈硬化特征。做报告时把这些适用范围讲清楚，能避免客户把结论用到错误场景。