2. OptiStruct复合材料建模基础
复合材料铺层优化能否得到可信结果，首先取决于有限元模型是否把“方向”和“厚度”表达正确。与各向同性金属不同，复合材料同一块层合板在不同方向上的刚度、强度差异可达数倍，因此 OptiStruct 中的单元类型、材料定义和坐标系设置三者必须相互匹配。本节围绕壳单元与实体单元的选择、各向异性材料卡片的填写、铺层坐标系与参考方向的建立展开，并给出可直接对照的建模示例。
一、壳单元与实体单元的选择依据
OptiStruct 中复合材料建模最常用的是壳单元（CQUAD4/CTRIA3）配合 PCOMP 属性，以及实体单元（CHEXA/CTETRA）配合多层实体属性。选择的核心判断依据是：厚度方向是否存在需要单独求解的应力分量，以及结构厚度与面内尺寸的比值。
- 壳单元 + PCOMP：适用于薄壁层合结构，如机翼蒙皮、车身覆盖件、风电叶片壳体。厚度方向应力近似为零，通过积分点存储每层材料方向与厚度，计算量小、铺层优化变量（各层角度、厚度）定义直观，是铺层优化的首选。
- 实体单元 + 多层实体：适用于厚截面、厚壁接头、层间应力显著的部位，如螺栓连接区、T 型接头、厚层合板边缘。厚度方向正应力与层间剪应力可被直接求解，但单元数量大，铺层优化变量数量随之增加。
- 混合建模：主体用壳单元，连接区局部用实体单元，通过 RBE3 或壳-实体耦合过渡，兼顾精度与效率。
一个实用的经验判据是：当结构最小面内尺寸与总厚度之比大于 20 时优先用壳单元；小于 10 且关注层间失效时用实体单元；介于两者之间需做壳/实体对比试算。示例中取一块 300 mm × 200 mm、总厚 2.4 mm 的层合板，比值远大于 20，因此采用壳单元建模。
二、各向异性材料定义
OptiStruct 用 MAT8 定义二维正交各向异性单层材料，关键字段包括：
MAT8 1 12000. 1000. .30 500. 500. 500.
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0
其中 E1、E2 为纤维方向与垂直纤维方向的弹性模量，NU12 为主泊松比，G12、G1Z、G2Z 为面内与面外剪切模量。示例中 E1=12000 MPa、E2=1000 MPa，说明该单层沿纤维方向刚度是横向的 12 倍，这正是必须正确设置参考方向的原因。若还需强度校核，应补充 XT、XC、YT、YC、S 等强度参数，并配合 PCOMP 中的 FT 失效准则字段。
PCOMP 属性卡用于定义铺层顺序，每行给出该层的厚度、材料编号和铺层角度：
PCOMP 10 2.4 0.0 0.0 0.0 0.0 0.0
+ 0.6 1 0.0 0.6 1 45.0
+ 0.6 1 -45.0 0.6 1 90.0
该示例为 [0/45/-45/90]s 对称层合板，单层厚 0.6 mm，共 8 层，总厚 2.4 mm。对称铺层可避免固化后的耦合翘曲，是工程中的常见做法。
三、铺层坐标系与参考方向设置
铺层角度是相对于参考方向定义的，因此必须先明确参考方向。OptiStruct 中常用两种方式：
- 单元坐标系：由 CQUAD4 的节点顺序决定，x 轴从节点 1 指向节点 2，z 轴为单元法向。角度 0° 表示纤维沿单元 x 轴。
- 材料坐标系：通过 MCID 字段引用 CORD2R/CORD2C 定义的局部坐标系，使参考方向与结构全局方向（如机身轴线、叶片展向）对齐，避免因网格划分方向不一致导致铺层角度混乱。
对于曲面结构，建议统一使用材料坐标系，并让坐标系 x 轴沿主受力方向。若网格单元 x 轴方向杂乱，即使铺层角度数值相同，实际纤维方向也会随单元变化，优化结果将失去物理意义。检查方法是：在 HyperMesh 中显示单元坐标系，观察箭头是否与预期纤维方向一致。
四、建模检查清单
完成上述设置后，建议按以下顺序检查模型，避免方向或厚度错误带入优化：
| 检查项 | 检查内容 | 判断依据 |
|---|---|---|
| 单元法向 | 所有壳单元法向一致朝外 | 法向翻转会导致铺层顺序颠倒 |
| 参考方向 | 单元坐标系或材料坐标系 x 轴与设计意图一致 | 显示箭头与纤维方向重合 |
| 铺层总厚 | PCOMP 各层厚度之和等于设计厚度 | 示例 8×0.6=2.4 mm |
| 材料参数 | E1、E2、NU12、G12 单位统一 | 建议统一用 MPa、mm |
| 对称性 | 对称铺层避免耦合 | [0/45/-45/90]s 满足 |
只有单元、材料、坐标系三者一致，后续的铺层角度和厚度优化变量才有明确的物理含义。实际建模中，建议先做一个单层板或简单层合板的试算，确认参考方向与预期一致后，再扩展到完整结构。