3. 铺层定义与PCOMP/PCOMPG卡片
在 OptiStruct 的复合材料建模中,铺层信息最终都要落到壳单元的属性卡片上。最核心的两张卡片是 PCOMP 和 PCOMPG:前者是经典的分层壳属性卡,后者是带全局铺层编号(Global Ply)的增强版本。理解它们如何把“层数、单层厚度、铺层角度、材料 ID”这四类信息关联起来,是后续做自由尺寸优化(Free-size)和铺层顺序优化(Shuffle)的前提。
一、PCOMP 的字段结构与关联逻辑
PCOMP 采用“属性卡 + 若干层定义行”的组织方式。属性卡本身定义整体行为,例如 PID(属性编号)、Z0(参考面相对中面的偏移)、NSM(非结构质量)、SB(剪切修正系数)、FT(失效准则)以及 GE(阻尼)。真正描述铺层的是紧随其后的层定义行,每一行对应一个铺层,字段顺序为:
PCOMP PID Z0 NSM SB FT TREF GE
MID1 T1 THETA1 SOUT1
MID2 T2 THETA2 SOUT2
...
其中 MID 是材料 ID,指向 MAT8 正交各向异性材料卡;T 是该层厚度;THETA 是铺层角度;SOUT 是输出请求标志。四者的关联方式可以概括为:层数 = 层定义行的数量,每一行独立绑定一个材料 ID、一个厚度和一个角度。总厚度等于所有层厚度之和,中面位置由 Z0 与各层厚度共同决定。
这里有一个容易出错的约定:角度 THETA 是相对于壳单元局部坐标系(由 CORD2S 或单元自身方向定义)的,逆时针为正。如果单元局部坐标系与全局坐标系不一致,同一角度在不同单元上代表的实际纤维方向就不同。因此在建模阶段必须统一参照系,并在全文和图示中保持一致。
二、PCOMPG 的差异:全局铺层编号
PCOMPG 的字段与 PCOMP 基本一致,但多了一个关键字段 GP(Global Ply ID)。它的层定义行形如:
PCOMPG PID Z0 NSM SB FT TREF GE
GP1 MID1 T1 THETA1 SOUT1
GP2 MID2 T2 THETA2 SOUT2
...
GP 的作用是给每一层一个跨属性、跨单元的全局身份。这样在优化过程中,即使某层的厚度或角度发生变化,软件仍能追踪“这是第几号铺层”,从而支持铺层顺序优化、制造约束(如对称、平衡、连续层数限制)以及结果的后处理映射。简单说:PCOMP 关注“这一层是什么”,PCOMPG 额外关注“这一层在整个结构中是谁”。
三、两种定义方法的适用场景
选择 PCOMP 还是 PCOMPG,取决于优化目标和制造约束的复杂度。
- PCOMP 适用场景:铺层顺序和层数在设计阶段就已确定,只需要做尺寸优化(如每层厚度微调)或简单的材料替换;模型规模较小、不需要跨属性追踪某一层;或者使用传统 Nastran 风格输入、与旧模型兼容。
- PCOMPG 适用场景:需要做自由尺寸优化(Free-size),让软件在单元级别决定哪里需要材料;需要做铺层顺序优化(Shuffle),在满足制造约束的前提下重新排列铺层;需要施加对称、平衡、连续层数、丢层限制等制造约束;需要把优化结果按全局铺层编号映射回每个单元。
四、建模实操中的关键检查点
无论使用哪种卡片,建模时都应检查以下几点:
- 材料 ID 必须指向 MAT8,而不是 MAT1。MAT8 提供 E1、E2、NU12、G12、G1Z、G2Z 等正交各向异性参数,是复合材料铺层的材料基础。
- 厚度单位与模型单位一致。如果模型用 mm,单层厚度通常取 0.125 mm 或 0.2 mm;如果模型用 m,则要写成 0.000125。单位不一致会导致总厚度和刚度完全错误。
- 角度参照系要统一。THETA 相对于单元局部坐标系,逆时针为正。如果使用 CORD2S 定义局部坐标系,要确认其 XY 平面与壳单元法向一致。
- Z0 与中面位置。Z0 表示参考面相对中面的偏移,默认 0 表示参考面在中面。如果壳单元偏置,需要正确设置 Z0,否则弯曲刚度会出错。
- SOUT 输出标志。如果后续要看每层的应力、应变或失效结果,需要在对应层行设置 SOUT,否则结果文件里可能没有该层数据。
五、一个可参考的建模示例
假设一块 300 mm × 200 mm 的层合板,材料为 T300/ epoxy,单层厚度 0.125 mm,铺层顺序为 [0/45/-45/90]s,共 8 层。用 PCOMPG 定义时,可以写成:
PCOMPG 100 0.0 0.0 0.833 0.0 0.0 0.0
1 201 0.125 0.0 1
2 201 0.125 45.0 1
3 201 0.125 -45.0 1
4 201 0.125 90.0 1
5 201 0.125 90.0 1
6 201 0.125 -45.0 1
7 201 0.125 45.0 1
8 201 0.125 0.0 1
这里 GP 从 1 到 8 连续编号,MID=201 指向 MAT8,总厚度为 8 × 0.125 = 1.0 mm。对称铺层通过 GP 的镜像关系体现,后续做铺层顺序优化时,软件可以识别 1-8、2-7、3-6、4-5 的对称对,并在制造约束中保持这种关系。
总结来说,PCOMP 和 PCOMPG 的核心差异不在字段数量,而在于是否给每一层赋予全局身份。层数、厚度、角度、材料 ID 四者的关联方式在两张卡片上是一致的:每一层行独立绑定一组参数,层数由行数决定,总厚度由各层厚度累加。选择哪张卡片,取决于你是否需要跨属性追踪铺层、是否需要做自由尺寸和铺层顺序优化。理解这一点,后续的优化设置和结果解读就会顺畅很多。