2. OptiStruct复合材料建模基础
在OptiStruct中建立复合材料模型，核心工作可以归结为三件事：把单层材料的方向性刚度写清楚（MAT8），把多层单层的堆叠顺序与角度写清楚（PCOMP/PCOMPG），再把角度和厚度变成可被优化器驱动的设计变量。这三步分别对应材料卡、属性卡和优化设计变量，缺一不可。本节围绕这三张卡片展开，说明它们各自的字段含义、相互引用关系，以及参数化的具体做法。
一、MAT8：正交各向异性单层材料
复合材料单层在面内是正交各向异性的，即沿纤维方向（1方向）与垂直纤维方向（2方向）的刚度差异很大。OptiStruct用MAT8卡片描述这种材料，它属于二维正交各向异性材料，主要字段包括：
- E1、E2：1方向与2方向的弹性模量，单位与模型单位制一致。以碳纤维环氧单层为例，E1常在120000～140000 MPa量级，E2常在8000～10000 MPa量级，两者相差一个数量级，这正是正交各向异性的直观体现。
- NU12：主泊松比，即1方向受拉时2方向收缩与1方向伸长的比值，典型值约0.25～0.35。
- G12、G1Z、G2Z：面内剪切模量与两个横向剪切模量。G12典型值约4000～6000 MPa。
- RHO：密度，用于质量与模态分析。
- X、Y、S等强度字段：用于失效判据（如Tsai-Wu、Hill），若只做刚度优化可暂不填写，但一旦涉及强度约束就必须补齐。
需要特别强调的是方向约定：MAT8的1方向是材料坐标系下的纤维方向，而单元在网格中的方向由单元坐标系决定。因此单层的实际纤维方向 = 材料坐标系1方向 + PCOMP中的铺层角度。这个叠加关系是后续角度参数化的基础，必须在建模之初就统一约定，避免同一模型中出现两套方向定义。
二、PCOMP与PCOMPG：铺层属性
PCOMP是经典铺层属性卡，它把若干单层按顺序堆叠成一个层合板，并赋予壳单元。其关键字段包括：
- MIDi：第i层的材料ID，指向MAT8。
- Ti：第i层的厚度。
- THETAi：第i层的铺层角度，单位为度，相对于材料坐标系1方向逆时针为正。
- SOUT：可选输出请求，用于提取每层的应力应变。
- NSM：非结构质量，用于模拟胶层、涂层等。
PCOMP的层序约定是：从壳单元法向正方向看，第一层位于最外侧（通常是外蒙皮），最后一层位于最内侧。这个顺序直接影响弯曲刚度与中性面位置，建模时不能随意颠倒。
PCOMPG是PCOMP的推广形式，它允许为每一层指定全局层编号（GPi），从而在多个属性之间共享同一套层定义。当结构由多个区域拼接、但需要统一追踪某一层的连续性时，PCOMPG比PCOMP更方便。例如机翼蒙皮与翼梁缘条共固化，用PCOMPG可以让第3层在蒙皮和缘条上保持同一编号，便于后处理与优化约束的施加。
三、铺层角度与厚度的参数化方法
参数化的目标是让优化器能够调整角度和厚度，而不是把数值写死。OptiStruct提供两条主要路径：
1. 厚度参数化（尺寸优化）：把PCOMP中每层的Ti设为DESVAR设计变量，通过DVPREL1建立设计变量与属性字段的线性关系。例如定义设计变量T1对应第1层厚度，初值0.125 mm，上下限0.08～0.20 mm，再用DVPREL1把T1关联到PCOMP第1层的T字段。优化器在满足应力、位移、屈曲等约束的前提下，自动分配各层厚度。
2. 角度参数化（自由尺寸/铺层顺序优化）：角度不能像厚度那样连续插值，因为0°与90°之间线性过渡会经过45°，而45°的刚度特性与两端完全不同。OptiStruct的做法是把角度离散为候选集，例如{0, ±45, 90}，通过拓扑优化或铺层顺序优化在候选集中选择。具体实现时，可用DESVAR定义角度变量，配合DVPREL1关联THETA字段，并设置离散选项；也可使用PCOMPG配合铺层顺序优化，让优化器直接决定每层的角度与顺序。
下面这张图展示了从MAT8到PCOMP再到设计变量的完整数据流关系，以及角度与厚度两条参数化路径的分叉点。
四、建模实操要点与示例参数
下面给出一个可直接参考的建模片段。假设单层材料为碳纤维环氧，单层固化厚度0.125 mm，铺层为[0/45/-45/90]s，共8层，总厚度1.0 mm。MAT8与PCOMP的填写如下：
MAT8 1 135000. 9000. 0.30 5000. 5000. 5000.
0.0 1600. 50. 200. 80.
PCOMP 10 1.0 0.0 1 0.125 0.0 1
0.125 45.0 1 0.125 -45.0 1 0.125
90.0 1 0.125 90.0 1 0.125 -45.0
1 0.125 45.0 1 0.125 0.0 1
0.125 0.0 1 0.125其中MAT8第1行给出E1、E2、NU12、G12、G1Z、G2Z，第2行给出密度与强度参数；PCOMP第1行给出总厚度、SOUT、层数，后续每三个字段为一组（MID、Ti、THETAi）。注意PCOMP的总厚度字段应等于各层厚度之和，OptiStruct会以此校验层定义的一致性。
参数化时，若把第1层厚度设为设计变量T1，初值0.125，下限0.08，上限0.20，则DVPREL1可写为：
DVPREL1 100 PCOMP 10 T1 0.0 1.0
T1 0.125这表示设计变量T1以系数1.0线性影响PCOMP第10号属性第1层的厚度字段。角度参数化则把T字段换成THETA字段，并配合DESVAR的DDVAL离散选项，把候选角度限定为{0, 45, -45, 90}。
五、结果观察与判断依据
模型建好后，可通过以下方式确认参数化是否正确生效：查看优化迭代历史中设计变量的变化曲线，确认厚度变量在上下限内收敛、角度变量落在候选集内；提取每层应力时，检查SOUT是否已开启，并核对第1层与最后一层的应力分布是否符合铺层顺序的预期。若厚度变量始终贴在上限或下限，说明约束过紧或初值不合理，需要调整约束或放宽边界。
本节的核心逻辑可以概括为：MAT8定义单层刚度方向，PCOMP/PCOMPG定义堆叠顺序与角度，DESVAR+DVPREL1把厚度和角度变成优化器可调的量。三者通过MID和PID形成引用链，任何一环的方向约定或层序出错，都会在后续优化中放大为不可信的结果。因此，建模阶段的方向统一与层序核对，比优化算法本身更值得投入时间。