古典積層理論（CLT）
理论与物理
什么是CLT
老师，经典层合理论（CLT: Classical Lamination Theory）是什么？
CLT是描述纤维增强复合材料（FRP）层合板力学行为的基础理论。它是将基尔霍夫板理论扩展到各向异性层合板，通过各层的材料特性和纤维方向角来计算整个层合板的刚度。
和各向同性板理论有什么不同？
各向同性板中，一个参数——弯曲刚度 $D = Et^3/12(1-\nu^2)$ 就能决定一切。而在CLT中，需要一个称为ABD矩阵的6×6矩阵。它包含了面内刚度、弯曲刚度以及面内-弯曲耦合。
ABD矩阵
CLT的核心概念——ABD矩阵：
其中：
- $\{N\}$ — 面内力合力（$N_x, N_y, N_{xy}$）
- $\{M\}$ — 弯曲力矩合力（$M_x, M_y, M_{xy}$）
- $\{\varepsilon^0\}$ — 中面应变
- $\{\kappa\}$ — 曲率
$[A]$, $[B]$, $[D]$ 分别代表什么？
$[B]$ 很有意思。在各向同性板中 $[B] = 0$，对吧？
没错。在各向同性板或对称铺层（例如: $[0/90]_s$）中，$[B] = 0$。而在非对称铺层中，$[B] \neq 0$，会产生拉伸导致板翘曲这种奇特的现象。
ABD矩阵的计算
各矩阵的计算公式：
其中 $\bar{Q}_{ij}^{(k)}$ 是第 $k$ 层的转换刚度矩阵（考虑纤维方向角），$z_k$ 是层的位置。
根据各层的位置和纤维方向角，就能计算出整个层合板的ABD矩阵了。
是的。CLT是从“各层材料特性+纤维方向角+铺层顺序”推导出“层合板宏观刚度”的理论。在使用FEM壳单元处理复合材料时，内部使用的就是CLT。
铺层表示法
像“$[0/90/\pm 45]_s$”这样的表示法该怎么读？
按从下到上的顺序描述纤维方向角：
- $[0/90/\pm 45]_s$ = $[0/90/+45/-45/-45/+45/90/0]$
- $s$ 表示对称（symmetric）。以中面为基准镜像
- 对称铺层则 $[B] = 0$（无面内-弯曲耦合）
对称铺层是设计的基本原则吗？
基本上是的。非对称铺层在成型后会产生翘曲（固化翘曲），在制造上也有问题。飞机结构原则上都采用对称铺层。
总结
我来整理一下CLT的理论。
要点：
- ABD矩阵描述层合板的刚度 — $[A]$: 面内，$[B]$: 耦合，$[D]$: 弯曲
- 根据各层的 $\bar{Q}$ 和位置 $z_k$ 计算ABD矩阵 — 输入是层的材料特性和纤维方向角
- $[B] = 0$ 则无面内-弯曲耦合 — 通过对称铺层实现
- 对称铺层是设计的基本原则 — 防止翘曲
- FEM的复合材料壳单元内部使用CLT — CLT是FEM的基础理论
如果不理解CLT，就无法设置FEM的复合材料分析，也无法解读结果，对吧？
正是如此。CLT是复合材料设计的“字母表”。没有它，就无法谈论复合材料的FEM分析。
CLT的起源与Reissner和Mindlin的贡献
经典层合理论（CLT）由Whitney、Leissa等人在1960-70年代确立。CLT假设各铺层处于平面应力状态，基于基尔霍夫板理论（面外变形保持法线刚性）进行积分，求得整体刚度矩阵[ABD]。Reissner（1945年）的剪切变形理论和Mindlin的改进（1951年）则发展成了弥补CLT“无剪切刚度”局限的一阶剪切变形理论（FSDT）。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重，越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施力，加速度可忽略”的假设。但在冲击载荷或振动问题中，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——用相同的力拉铁棒和橡皮筋，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的性质就是杨氏模量 $E$，它决定了刚度。常见的误解是：“刚度高=强度高”。不对，刚度是“不易变形”，强度是“不易破坏”，这是两个不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力……这些都是外力。这里容易犯的错误是：搞错载荷方向。本想施加“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生这种情况。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃下去。实际上不会这样，所以设置适当的阻尼非常重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷、弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
FEM中的CLT实现
FEM的壳单元中是如何实现CLT的？
复合材料的壳单元中，根据各层的 $\bar{Q}$ 矩阵和铺层信息（纤维方向角、厚度、位置）构成ABD矩阵，并反映到单元的刚度矩阵中。厚度方向的积分点配置在各层内。
Nastran
```
PCOMP, 1, , , , SYM
, 1, 0.125, 0., YES,
, 1, 0.125, 90., YES,
, 1, 0.125, 45., YES,
, 1, 0.125, -45., YES
```
使用PCOMP卡片指定各层的材料ID、厚度、纤维方向角。SYM表示对称铺层。
Abaqus
```
*SHELL SECTION, COMPOSITE, ELSET=panel
0.125, 3, CFRP, 0.
0.125, 3, CFRP, 90.
0.125, 3, CFRP, 45.
0.125, 3, CFRP, -45.
```
每行定义一个铺层。依次为：厚度、积分点数、材料名、纤维方向角。
Ansys
在Workbench中，使用“ACP（Ansys Composite PrePost）”定义铺层。通过GUI操作可以直观地设置铺覆和纤维方向角。
为什么Nastran的PCOMP在航空航天领域被广泛使用？
PCOMP有自1980年代以来的历史，在航空航天认证（型号合格审定）中拥有大量的验证实绩。可以直接输出各层的应变、应力，并且与破坏判据（Tsai-Wu、Hashin等）的联动也很完善。
厚度方向的积分点
每层需要多少个积分点？
推荐每层厚度方向至少3个点（辛普森积分）。处理弹塑性时，需要5个点以上。
全部层共 $n$ 层 × 3点 = $3n$ 点。20层的层合板就是60点。厚度方向的积分点越多，计算成本越高，但对精度很重要。
材料坐标系
纤维方向角是相对于哪个坐标系定义的？
材料坐标系是为每一层单独定义的。通常以壳单元的面内方向（1轴方向）为基准指定纤维方向角。当单元方向改变时，纤维方向角的基准也会改变，因此确认材料方向至关重要。
要准确表现铺覆（曲面上纤维方向角的变化），需要为每个单元单独设置材料坐标系。ACP（Ansys）或Fibersim（Siemens）可以自动化这项工作。
总结
なった
詳しく
報告