NAFEMS LE5：Z形截面悬臂梁扭转
基准算例的规格与理论背景
LE5验证的是什么
NAFEMS LE5是壳单元的验证问题：在Z形截面悬臂梁的自由端施加扭转力偶。有趣之处在于，载荷明明是扭转，合格判定却用根部附近的“轴向应力”来做。像Z形这样的开口截面，圣维南扭转刚度极小，只能靠翘曲（warping）受到约束这一机制——两侧翼缘被反向弯曲——来承担扭矩。其结果是翼缘中产生很大的轴向膜应力。LE5要考的，正是壳单元的膜性能能否正确再现“开口截面的约束扭转”这一结构力学要点。
问题设定
| 项目 | 内容 |
|---|---|
| 几何 | Z形截面悬臂梁：长10 m、腹板高2 m、翼缘宽1 m、板厚0.1 m（薄壁，用壳建模） |
| 材料 | \( E = 210 \) GPa、\( \nu = 0.3 \) |
| 载荷 | 自由端相当于扭矩1.2 MN·m：在上下翼缘的端边上各施加均匀分布的剪力0.6 MN且方向相反（构成力偶） |
| 约束 | 根部端面固定 |
| 目标值 | 点A（根部附近的规定点）中面（膜）轴向应力 \( \sigma_{xx} = -108 \) MPa |
为什么是壳的好题
明明是扭转问题，却用轴向应力判定，有点不可思议。这不是板弯曲的验证吗？
妙就妙在LE5的这个设计上。目标值取的是中面的膜应力，因此被验证的不是板的弯曲性能，而是作为膜的面内行为＋整个截面的力流。壳单元是“膜＋弯曲＋横向剪切”的合成体，单元列式的差别（减缩积分、面内转动自由度、抗闭锁措施）只有在复杂的力流中才会暴露出来。Z形截面的约束扭转会在翼缘与腹板之间来回传递膜力，形成一种恰到好处的复杂力流——单纯的拉伸或纯弯曲拉不开的单元实力差距，在这里会变成实实在在的数字。
实施步骤与读数要点
标准的实施步骤
- 建模——在中面上建立壳模型（板厚0.1 m以属性方式赋予）。确认翼缘与腹板的连接线是共享的
- 载荷——在上下翼缘自由端的边上，以沿边均匀的分布剪力（边载荷）各施加0.6 MN且方向相反。务必校核合力与合力矩
- 约束——根部端面全自由度固定
- 网格序列——从规格推荐的划分开始，再取2倍、4倍确认点A膜应力的收敛
- 读数——提取点A的中面（membrane/mid-surface）应力的轴向分量，报告其向−108 MPa的收敛过程
最常见的错误——表面、背面与中面弄混
壳的应力输出有上表面（top）、下表面（bottom）、中面（mid）三层。LE5的目标值是膜应力（mid）——若读取含板弯曲分量的top/bottom，就会系统性地偏离出一个弯曲量。确认工具默认输出的是哪一层（以及“表面”指的是哪一侧，即壳法线的朝向），是这道题的第一道关口。实务中“壳的应力是哪一层的数值”同样是报告书的必写事项——LE5也是把这条纪律刻进身体的题目。
载荷等效替换的规矩
把“扭矩1.2 MN·m”当作一个点上的扭矩施加是错误的，必须按规格作为翼缘端边上均匀分布剪力构成的力偶来给。集中扭矩会带来局部变形与奇异性，并影响到根部的应力。作为等效替换的校核，要在载荷汇总中确认：①所施加载荷的合力为零（纯力偶）；②合力矩为1.2 MN·m。这种“替换载荷后就校核合力与合力矩”的习惯，作为边界条件验证的载荷版，在实务全局都很有效。
实务中的活用方法
作为壳单元的验收试验
LE5最适合用于单元类型与列式的比较。在相同网格密度下把①1阶减缩积分壳、②1阶完全积分壳、③2阶壳并列起来，比较点A应力的收敛情况，就能把公司标准单元的实力（以及它在粗网格下的脾气）定量化。求解器更新、单元设置变更（面内转动刚度、沙漏控制参数）后的回归确认，也很适合用这道轻量的题目。把它与T4、LE11并列，把“传热、热应力、壳”这三件套作为验证台账的基本框架，是推荐的配置。
对结构设计的教训——开口截面不耐扭
LE5的物理与实际设计中的重要教训直接相连。开口截面（Z形、C形、L形）的圣维南扭转刚度只有闭口截面的几十分之一到几百分之一，扭矩要靠翘曲约束产生的轴向应力来承担，因而表现为“一扭就被弯曲应力破坏”。薄壁构件的抗扭设计中把截面封闭（做成箱形）为何如此有效，翘曲受约束的部位（根部、隔板处）为何会出现应力集中——解过LE5的经验，会把这些变成定量的直觉。基准算例不仅服务于验证，同时也是结构力学的活教材。
与梁理论（Vlasov）的对照
约束扭转可以用符拉索夫的薄壁梁理论（翘曲函数、双力矩）手算，把它与壳的解对照，就构成“壳FEM对经典理论”的相互验证。根部附近是圣维南原理不成立的区域（翘曲约束的影响占主导），因此与理论不可能完全一致，但用于确认应力的量级与分布趋势已经足够——作为“用另一套理论体系去校核FEM答案”的训练题，它同样有价值。
各工具的设置要点
设置对照
| 工具 | 边载荷 | 膜应力的输出 |
|---|---|---|
| Ansys Mechanical | 施加到边上的Force（按分布指定） | Shell中的Membrane Stress／指定中面 |
| Abaqus | *SHELL EDGE LOAD（分布剪力） | 指定截面点（SPOS/SNEG/中央），或由膜力SF换算 |
| Nastran系 | 施加到边节点的等效节点力（由均匀分布换算） | 由单元力（膜力）输出F1按 σ=F/t 换算最为可靠 |
| 支持壳单元的通用FEM | 若没有边载荷功能，则用等效节点力（端部节点取一半） | 确认“Z1/Z2/mid”之类的层指定 |
手工构造等效节点力时的注意事项
把边上的均匀分布载荷折算为节点力时，要遵循一致节点力的规则：1阶单元是内部节点等分、两端节点取一半（2阶单元则由中间节点分得更大权重）。给所有节点同一数值的“平均分摊”会在端部破坏分布形态，从而扭曲局部应力。这一点虽然大体还在“替换载荷后校核合力与合力矩”的保护范围内，但分布形态的失真会从这类校核中溜过去，因此必须对分配规则本身保持自觉。
壳验证的现状
壳列式的谱系与基准算例的角色
壳单元的发展史，就是剪切闭锁与膜闭锁对策（减缩积分＋沙漏控制、MITC系的混合插值、EAS法）的发展史本身，而包含LE5在内的标准基准算例群（受挤压圆筒、Scordelis-Lo屋盖、扭曲板等“壳的必做测试”）一直扮演着裁判的角色。评价新单元、新求解器时把这套经典题目跑一遍，是单元开发界的惯例；用户一侧也能用同一套题目掌握“自己所用单元的谱系与脾气”。
实体壳与免中面化的潮流
由于不愿花工时从CAD实体中抽取中面，用实体壳单元或厚度方向布置1个单元的实体直接求解薄壁件的做法正在增多。像LE5这样的基准算例，可以作为衡量“这种简化与单元选择能否保住膜行为精度”的尺子，因而也能用于建模方针的决策。用中面壳与实体壳分别求解LE5的内部对比研究，会成为方针转变时颇具说服力的依据。
接入自动验证
与T4、LE11一样，LE5也因轻量而适合CI与回归测试。单元库的更新、默认参数的变更都可能悄悄影响膜行为，因此把它作为“壳的基准器”并设置±1%容差的自动判定，会让环境更新时的安心感完全不同。
故障排查
“算不出−108 MPa”时的诊断表
| 症状 | 可能原因 | 对策 |
|---|---|---|
| 偏差达±几十% | 读取的是top/bottom应力（混入了弯曲分量） | 指定中面（膜）应力后重新提取 |
| 符号相反 | 点A取错（上翼缘与下翼缘弄反）、载荷方向不对 | 重新确认规格中点A的坐标与载荷方向 |
| 差了一倍或一半 | 力偶只施加了一侧／把0.6 MN误解为两侧的合计 | 在载荷汇总中校核合力矩=1.2 MN·m |
| 粗网格偏差很大且收敛缓慢 | 1阶减缩积分在粗网格下的脾气、翼缘宽度方向划分不足 | 宽度方向至少分4份，并与2阶单元比较 |
| 局部应力出现剧烈波动 | 集中扭矩、节点力分配不当、根部的奇异性 | 改用分布式边载荷。不要取根部角点的数值，而在点A的规定位置评价 |
| 工具之间残留几%的差异 | 面内转动刚度、剪切修正、应力恢复方法的差别 | 写明单元列式的设置再比较。用收敛值之间相互判定 |
验证台账“基本三件套”的完成
T4、LE11、LE5这三道题我都跑通了。接下来还应该补些什么呢？
帮助
更多
错误