积层板 — CAE用语解说
积层板
老师，积层板在CFRP分析中经常出现。应该怎样处理呢？
积层板的理论基础
积层板的基本概念
积层板与简单的多层板叠放有什么区别？只是用胶黏剂粘合在一起而已吧。
根本区别在于各层（层片）的材料和纤维取向可以独立设计。例如，在飞机主翼蒙皮中，使用碳纤维强化塑料（CFRP），0度层负责拉伸刚性，±45度层负责剪切刚性，90度层负责横向刚性，各自担当其职。简单的板层叠放无法实现这种异向性的控制。
所谓异向性控制，就是说材料的特性因方向而大不相同？那用什么数学公式来表示呢？
正是这样。单层被视为正交异向性。用应力-应变关系的一般形式，通过刚性矩阵[C]表示为
"积分"是指沿板厚方向加起来的意思吗？具体是怎样的计算呢？
正是沿板厚方向的积分。结果是单位宽度的面内力{N}和弯曲矩{M}与板中面应变{ε⁰}和曲率{κ}相关联。这就是积层板的构成方程，表示为
积层板的数值计算手法
FEM中的建模手法
在FEM中分析积层板时，每一层都用3D单元按实际厚度模型化非常麻烦。通常怎样做呢？
你说得对。特别是在有数十层的情况下不实际。标准做法是给壳元素（例如Abaqus的S4R）分配"复合材料截面"或"积层截面"。此时单元表示中面，厚度方向的积分在单元公式内根据指定的层配置（打印顺序）进行。积分点通常在各层内设置多个（例如Simpson规则3点）。
各层有积分点是说某一层的下面可以计算应力吗？层间剥离（分层）怎样评价？
是的，通过各层积分点输出（Abaqus的SDV）可以得到下面、中间、上面的应力。对分层的直接评价，要么在层间插入黏着单元，要么使用基于应力的判定准则。后者的代表例是层间剪切应力（Interlaminar Shear Stress, ILSS）的评价。例如，航空航天业界常设层间剪切应力许容值为10MPa以下，超过这个值分层风险就会升高。
层间剪切应力可以从壳元素输出直接得到吗？还是要另外计算？
标准壳元素输出以面内应力为中心，层间剪切应力（τ_xz、τ_yz）通常不能直接输出。这些通常要通过后处理从平衡方程计算得到。具体是用
积层板的实务应用
分析工作流程与验证
开始积层板分析时，首先要定义什么？我不知道顺序。
有既定的工作流程。首先(1)定义单层（UD层）的所有9个材料特性。E1、E2、ν12、G12、G13、G23以及强度值（Xt、Xc、Yt、Yc、S12等）。然后(2)定义积层顺序（例：[0/45/90/-45]s）和各层厚度（例：各层0.125mm）。(3)定义单元坐标系相对于材料坐标系（纤维方向）的关系。这些任何一个搞错都会导致全部出错。
需要9个材料特性？我的教科书里只有E1、E2、ν12、G12这4个。
好问题。这4个是定义面内正交异向性的最少集合。但要正确评价弯曲和层间剪切，还需要面外剪切刚性G13和G23。这两个常常与面内剪切刚性G12的值不同。无谓地假设G12=G13=G23会导致层间剪切应力被低估，风险很大。例如某种CFRP，G12=5GPa，但G13=G23=4GPa。其他的是强度判定用。
分析结果怎样验证才能令人信服？应该从简单的拉伸试验模拟开始吗？
对，要从可验证的简单情况开始。具体步骤：1. 单向积层板（[0]8）拉伸，验证E1与材料输入值一致。2. [90]8拉伸确认E2。3. 对称积层板（如[0/90]2s）拉伸，与泊松比或面内剪切试验对比。这些都可以容易地与实验数据或古典积层板理论的手算结果对比。Ansys和Abaqus都提供这样的验证用教程。
实务中，还会进行积层顺序最优化吗？
当然会。特别是在轻量化被重视的航空航天和汽车领域是必要工程。使用工具（如Altair HyperStudy或Ansys OptiSlang）在约束条件（挠度、固有频率、屈曲荷载）下，以目标函数（重量或成本）最小为目标，把层角度和顺序作为变量进行最优化。但必定要加上制造性约束（连续同方向层不超过4层等）。
积层板软件比较
主要求解器的功能与特点
Ansys、Abaqus、MSC Nastran用于积层板分析时怎样选用呢？
按业界、工作流程、求解功能来分。Ansys（特别是ACP）的前后处理功能强大，对复杂曲面的积层定义和参数研究有优势，在汽车和风电业界受欢迎。Abaqus在非线性（材料、接触、大变形）分析方面压倒优势，在破坏和冲击模拟中不可或缺。MSC Nastran长期是航空航天业界事实上的标准，其求解器（SOL 105、400）和规范（NASA、ESA基准）的一致性是强项。
COMSOL怎样？说是多物理场软件，积层板分析可以吗？
可以，但应用领域不同。COMSOL的"复合材料模块"着眼于积层板的"多物理场行为"。例如热膨胀造成的翘曲（湿度问题）、积层板内压电层的驱动/传感功能、电磁波屏蔽效果分析等。纯粹的结构强度或屈曲的详细评价反而不是它的重点，而在处理这类复合现象时就能显其威力。
免费或低价软件（如CalculiX、Code_Aster）的积层板分析怎样？
核心求解功能（壳元素中的积层定义）是有的。CalculiX和Code_Aster也能用*SHELL SECTION或COQUE_3D定义层。但商用软件那样直观的前后处理GUI（层可视化、结果按层提取、破坏判定自动计算）几乎没有。所有输入数据（材料、积层顺序、积分点规则）都要用文本准确记述，学习成本和人工失误风险都高。适合研究和学习用。
积层板故障排除
常见错误与对策
积层板分析中施加弯曲时，变形形状有时会出异常。刚性太大还是太小，我难以判断。
首先怀疑"材料坐标系定义错误"。这是最常见的错误。本以为单元的1方向是纤维方向（0度），但由于单元方向的变化，实际材料坐标系可能旋转成意外角度，频繁发生。对策是分析前要把全部单元的材料方向可视化，用向量显示确认。Abaqus可用"Orientation"图、Ansys可用"Material Orientation"显示来检查。
还有，面外变形（如圆盘弯曲）的结果，与理论值或其他软件对比，变得极软。这是为什么？
这是"剪切锁定"或"面外剪切刚性低估"导致。首先确认使用的壳元素是否考虑厚度方向剪切变形（基于Reissner-Mindlin板理论）。其次确认面外剪切刚性G13、G23的输入值是否恰当。你没有用G12的值代替这些吧？再者，单元扭曲（宽纵比不好）时会发生剪切锁定，刚性被低估。改善网格正形或用低积分点单元（如Abaqus的S4R）有时能改善。
特征值屈曲分析中，即使计算到很高次的模式（如100次），临界荷载几乎不下降。这收敛了吗？
这是"局部屈曲模式"很多出现的典型现象。在薄壁无芯心的实心层板中，特别是90度层会个别波浪状屈曲，出现无数个模式。这些模式对实际结构行为没有贡献，只浪费计算资源。对策是明确分析目的。若想知道整体屈曲模式，就将模式形状可视化，只评价整体变形的低次模式（最初5～10个）。若要评价局部屈曲，就需要详细模型化那一层。
破坏判定（Tsai-Wu等）中，某些特定单元的破坏指数异常高。这是实际薄弱部分，还是计算噪声？
首先检查该单元是否在"应力奇异点"附近。例如，裂纹先端、陡峭的几何变化部、集中荷载作用点下方等。理论上这些点应力是无穷大，网格越细破坏指数越发散。实务中要参考"结构意义的区域"（离这些点稍远处）的值。积层板的自由端（边缘）层间应力也有奇异性，会出现同样问题。要评价边缘效应，需在边缘周围特别细化网格。
相关主题
价值
细节
错误