3. 中面抽取与理想化：薄壁连接件的中面抽取方法，壳单元与实体单元的选用原则
做光机结构仿真，薄壁连接件是最让人头疼的零件之一。为什么？因为它的厚度远小于长宽尺寸，直接用实体单元去画网格，网格数量会爆炸，计算时间长得让人崩溃。我早年吃过这个亏——一个简单的铝合金支架，厚度才2mm，我用四面体实体网格一划，好家伙，30万个单元，算了一整夜还没收敛。
后来我学乖了。薄壁件，就该用壳单元。但壳单元的前提是——你得先把中面抽出来。
3.1 什么是中面？为什么要抽中面？
中面，说白了就是零件厚度方向上的那个“中间面”。对于等厚度的薄壁件，中面就是上下表面正中间的那个虚拟曲面。
用壳单元模拟薄壁件时，我们其实是在中面上建立网格，然后把厚度信息作为属性赋给这些网格。这样做的好处很明显：
- 网格数量大幅减少——从三维降到了二维，节点数能少一个数量级
- 计算效率提升——同样的模型，壳单元比实体单元快5~10倍
- 弯曲行为更准确——壳单元天生擅长模拟薄壁结构的弯曲变形
核心原则：当零件的厚度与最小特征尺寸之比小于1/10时，优先考虑使用壳单元。这个比值越小，壳单元的精度优势越明显。
3.2 中面抽取的三种主流方法
我常用的中面抽取方法有三种，具体用哪个，取决于你手里的软件和模型复杂度。
| 方法 | 适用场景 | 精度 | 操作难度 |
|---|---|---|---|
| 自动中面抽取 | 等厚度、形状规则的薄壁件 | 中等 | 低 |
| 偏置曲面法 | 不等厚度、曲率变化大的零件 | 高 | 中 |
| 手动构建法 | 复杂曲面、修补失败的区域 | 最高 | 高 |
3.2.1 自动中面抽取
大多数CAE软件（比如Ansys Workbench、Abaqus、HyperMesh）都自带自动中面抽取功能。操作很简单：选中零件，点一下“抽取中面”，软件会自动识别上下表面并生成中间曲面。
但这里有个坑——自动抽取对模型质量要求很高。如果零件上有小圆角、倒角、台阶面，或者模型是从其他格式转换过来的（比如STEP转IGES），自动抽取很容易失败。
我曾经遇到过：一个从SolidWorks导出的STEP文件，里面有个0.5mm的微小圆角。自动中面抽取直接报错，我排查了半天才发现是那个圆角在作怪。后来我把圆角抑制掉，中面就顺利抽出来了。所以，做中面抽取前，先清理模型上的小特征。
3.2.2 偏置曲面法
当零件厚度不均匀时，自动抽取就不太靠谱了。这时候我习惯用偏置曲面法：
- 先提取零件的一个外表面（通常是受力面或安装面）
- 测量该区域的局部厚度
- 将外表面沿法向偏置半个厚度距离
- 得到的就是中面
这个方法的好处是——你可以分段处理。比如一个支架，底板厚3mm，筋板厚2mm，你可以分别提取各自的外表面，然后按不同的厚度偏置。最后再把这几块中面拼接起来。
3.2.3 手动构建法
这是最后的杀手锏。当自动和偏置都搞不定时，我就手动来。具体做法是：
- 在关键截面位置画一条中心线
- 沿着中心线扫掠出一个曲面
- 用这个曲面去切割零件，得到中面
手动构建虽然慢，但精度最高。我一般只在连接件的关键区域（比如螺栓孔周围、焊缝附近）用手动法，其他区域用自动法，这样效率和质量都能兼顾。
3.3 壳单元与实体单元的选用原则
这个问题，我几乎每次培训都会被问到。其实没那么复杂，记住下面这张表就行：
| 判断条件 | 推荐单元类型 | 原因 |
|---|---|---|
| 厚度/特征尺寸 < 1/10 | 壳单元 | 弯曲主导，壳单元精度高、效率高 |
| 厚度/特征尺寸 > 1/5 | 实体单元 | 厚度方向应力不可忽略，实体更准确 |
| 关注局部应力集中（如螺栓孔） | 实体单元 | 壳单元无法准确模拟厚度方向的应力梯度 |
| 关注整体变形和刚度 | 壳单元 | 壳单元对整体弯曲刚度的计算非常精准 |
| 接触分析（面-面接触） | 壳单元 | 壳单元的接触算法更成熟，收敛性更好 |
| 需要输出厚度方向应力 | 实体单元 | 壳单元默认不输出厚度方向的应力分量 |
我的个人习惯：对于光机结构中的薄壁连接件，我通常先用壳单元做全局分析，找出高应力区域。然后把这些区域单独切出来，用实体单元做子模型分析。这样既保证了整体效率，又拿到了局部细节的准确结果。
3.4 中面抽取后的理想化处理
中面抽出来之后，别急着画网格。先做几步理想化处理，能省掉后面很多麻烦。
3.4.1 几何清理
中面上经常会有一些碎面、小孔、尖角。这些特征对计算结果影响很小，但会让网格质量变差。我的做法是：
- 删除直径小于3倍厚度的孔
- 合并相邻的碎面（角度差小于5°的可以直接合并）
- 修复中面上的裂缝和重叠区域
3.4.2 厚度赋值
中面本身没有厚度信息，你需要手动把厚度赋上去。对于等厚度零件，直接赋一个全局厚度值就行。对于变厚度零件，可以用“厚度场”的方式逐区域赋值。
; 示例：在HyperMesh中为变厚度中面赋值
; 区域1：底板，厚度3mm
*setproperty(pshell, 1, "thickness", 3.0)
; 区域2：筋板，厚度2mm
*setproperty(pshell, 2, "thickness", 2.0)
; 区域3：安装凸台，厚度4mm
*setproperty(pshell, 3, "thickness", 4.0)
3.4.3 偏移设置
壳单元的厚度偏移方式有三种：中面偏移、顶面偏移、底面偏移。默认是中面偏移，但有时候你需要调整。
举个例子：一个连接件通过螺栓固定在基座上，螺栓预紧力作用在顶面。如果你用中面偏移，力的作用点会偏离实际位置，导致局部应力计算不准。这时候应该用顶面偏移，让壳单元的顶面与实际受力面对齐。
注意：偏移设置会影响接触定义和约束施加。我建议在建模初期就确定好偏移方式，中途更改容易导致接触对错位。
3.5 知识体系总览
下面这张图把本章的核心逻辑串起来了。你对照着看，思路会更清晰。
3.6 避坑指南
最后，我把自己这些年踩过的坑总结一下，你遇到了可以直接拿来用：
- 中面与实体混合建模要谨慎——壳单元节点有6个自由度，实体单元节点只有3个。两者连接处需要做自由度耦合，否则算出来结果会飞掉。
- 壳单元的厚度方向要一致——同一个零件上，所有壳单元的厚度方向必须统一指向同一侧。方向不一致会导致接触定义混乱。
- 不要在中面上画太细的网格——壳单元的网格尺寸一般取厚度的3~5倍就够了。画太细不仅浪费计算资源，反而可能引入数值误差。
- 注意壳单元的剪切锁定——对于非常薄的壳（厚度<0.5mm），一阶壳单元容易出现剪切锁定，导致刚度偏大。这时候改用二阶壳单元或者减缩积分单元。
我曾经犯过一个低级错误：把一个2mm厚的支架用壳单元建模，但忘了赋厚度值。结果算出来的变形量是实际值的100倍。我查了三天才找到原因——壳单元默认厚度是1mm，而我实际需要的是2mm。从那以后，我每次建模都会在第一步检查材料属性和截面属性。
嗯，中面抽取和单元选用这块，核心就是这些。你只要记住：薄壁件用壳，厚壁件用实体，关键区域做子模型。这个思路能帮你解决90%以上的光机结构连接件建模问题。