正交异向性 — CAE术语解说

分类: 术语集 | 2026-01-15

正交异向性的理论基础

正交异向性的物理意义

🧑‍🎓

"正交异向性"这个概念经常出现，但它与各向同性具体有什么不同呢？

🎓

根本区别在于材料特性是否随方向变化。各向同性材料在任何方向上的杨氏模量和泊松比都相同。而正交异向性材料在互相正交的3个主方向（例如纤维方向、面内垂直方向、厚度方向）上特性不同。以CFRP（碳纤维增强塑料）为例，纤维方向的杨氏模量约为150 GPa，而垂直于纤维的方向仅约10 GPa，相差10倍以上。

🧑‍🎓

明白了，三个方向特性不同。那用数学表达式来说，参数会比各向同性多吗？

🎓

完全正确。各向同性的线性弹性只需要杨氏模量E和泊松比ν两个参数。而正交异向性的完整线性弹性本构关系需要9个独立材料常数。包括3个方向的杨氏模量

$$ E_1, E_2, E_3 $$

、剪切弹性系数

$$ G_{23}, G_{31}, G_{12} $$

和泊松比

$$ \nu_{12}, \nu_{23}, \nu_{31} $$

。其中

$$ \nu_{ij} $$

表示j方向应力引起的i方向应变。

🧑‍🎓

那如果要通过实验测出这9个常数，工作量不是很大？实际设计中怎么处理？

🎓

这是个很好的问题。实际上，厚度方向的特性

$$ E_3, G_{23}, G_{31} $$

很难直接测量。因此在薄板结构（壳体单元）分析中，通常只关注面内特性（1、2方向）和面内剪切

$$ G_{12} $$

。厚度方向的剪切通常用"横剪切修正系数"简化处理。航空航天工业通常遵循ASTM D3039（拉伸试验）和ASTM D3518（面内剪切试验）等标准来获取数据。

正交异向性的数值计算方法

CAE软件中的输入与坐标系

🧑‍🎓

在CAE软件中定义正交异向性材料时，是不是直接有一栏输入这9个常数？

🎓

是的，主流CAE软件都有相应的界面。比如Abaqus，可以在"Engineering Constants"项中直接输入

$$ E_1, E_2, E_3, \nu_{12}, \nu_{13}, \nu_{23}, G_{12}, G_{13}, G_{23} $$

这9个参数。Ansys Mechanical也有"Orthotropic Elasticity"选项。但最关键的是，这些常数必须相对于"材料坐标系"定义。

🧑‍🎓

"材料坐标系"具体怎么设置？如果零件形状很复杂，纤维方向（1方向）在整体上很难统一定义啊。

🎓

这正是正交异向性分析的关键。每个软件有不同的做法，但都支持"局部坐标系分配"。Ansys中定义"Coordinate Systems"然后分配给材料；Abaqus用"Orientation"功能。例如在航空器复合材料面板中，面板面内方向为1、2方向，法向为3方向。对于曲面壳体，壳体的法向自动成为3方向。

🧑‍🎓

如果每个单元的坐标系都不同，求解器内部怎么处理全局刚度矩阵？是不是每次都要坐标变换？

🎓

完全正确。单元刚度矩阵

$$ [k_e] $$

首先在材料坐标系（局部坐标系）中计算，然后通过方向余弦矩阵

$$ [T] $$

从材料坐标系变换到全局坐标系：

$$ [k_e]_{global} = [T]^T [k_e]_{local} [T] $$

，最后装配到全局刚度矩阵。这样即使纤维方向沿零件形状弯曲，也能正确求解。

正交异向性的实际应用

分析工作流程与验证

🧑‍🎓

开始正交异向性材料分析前，首先要检查什么？

🎓

第一步要检查"材料常数的物理合理性"。输入的9个常数必须满足热力学允许条件，即任何加载下应变能都是正的。具体来说，杨氏模量和剪切弹性系数必须全为正

$$ E_i > 0, G_{ij} > 0 $$

，泊松比需满足

$$ |\nu_{ij}| < \sqrt{(E_i / E_j)} $$

。忽视这一点会导致分析发散或产生不现实的结果。

🧑‍🎓

材料常数没问题，但结果还是有异常，接下来该怀疑什么？

🎓

几乎肯定是"材料坐标系设置错误"。这是正交异向性分析最常见的问题。验证方法：对简单的拉伸或弯曲施加载荷，检查刚性较大方向（纤维方向）是否应力集中较高、变形较小。例如对CFRP板，沿纤维方向（1方向）拉伸的位移应该是沿垂直方向（2方向）拉伸的1/10以下。如果没有这种差异，说明坐标系有问题。

🧑‍🎓

对于像积层板这样每层纤维角度都不同的复合材料，怎么处理？

🎓

这种情况下，每个单层（lamina）都定义为正交异向性材料，然后指定层序（Ply Stacking Sequence）和各层的纤维角度（如0°、45°、-45°、90°）。Abaqus的"Composite Layup"、Ansys的"ACP (ANSYS Composite PrepPost)"或"Shell Layup"都支持这些功能。这些工具会自动按指定角度旋转各层的材料坐标系，并计算整个积层的等效刚度。

正交异向性的软件比较

主要软件的功能与特点

🧑‍🎓

Ansys、Abaqus、COMSOL在处理正交异向性材料上有区别吗？

🎓

基本的材料定义功能都有，但各有所长。Ansys在结构分析（Mechanical）的基础上，配备强大的复合材料专用前后处理模块"ACP"，特别在汽车和航空航天工业中广泛应用，甚至可进行详细的破坏分析。Abaqus的"Composite Layup"功能也很完善，在非线性行为（如损伤演化）的高级建模上有优势。

🧑‍🎓

COMSOL在这方面怎么样？

🎓

COMSOL Multiphysics的独特优势是可以轻松将正交异向性概念扩展到"结构"以外的物理场。例如热传导分析中，可定义纤维方向和垂直方向热传导率不同的"热正交异向性"。电导、扩散等现象也支持。此外，材料坐标系定义非常灵活，可用任意公式或函数指定方向，适合复杂织物结构建模。

🧑‍🎓

免费或低成本的CAE软件（如CalculiX、Code_Aster）能处理正交异向性吗？

🎓

可以的。CalculiX（免费求解器）在输入文件（.inp）中用*ELASTIC, TYPE=ORTHO关键字定义9个常数；Code_Aster（Salome-MECA环境）用ELAS命令的ORTH关键字。但它们通常不具备商用软件那样的图形化直观积层板定义和高级破坏准则，用户需要通过脚本或手工配置详细参数。

正交异向性的故障排除

常见错误及其原因

🧑‍🎓

分析运行时出现"负对角项"或"非正定矩阵"错误，求解器停止了。但我只输入了正数啊。

🎓

这是因为虽然各常数单独为正，但组合起来违反了热力学稳定条件（应变能可能为负）。问题通常出在泊松比。检查是否满足

$$ |\nu_{12}| < \sqrt{E_1/E_2} $$

。例如

$$ E_1=150GPa, E_2=10GPa $$

时，

$$ \sqrt{150/10} \approx 3.87 $$

，所以

$$ \nu_{12} $$

必须小于3.87。通常泊松比在0.5以下，但正交异向性中这个上限可能更大。

🧑‍🎓

应力集中区的剪切应力很高，超过了材料的剪切强度。这个结果能信吗？

🎓

这可能是"正交异向性材料的应力奇点"问题。异向性比（

$$ E_1/E_2 $$

）越大的材料，切口或孔边的应力梯度会比各向同性材料高得多，这在理论上是已知的。网格越细，剪切应力越高，甚至可能不收敛。在这种情况下，不能直接用线性弹性的应力结果进行强度评估，应该用平均化应力或从一开始就采用包含损伤准则的非线性分析。

🧑‍🎓

用壳体单元的积层板分析中，有人说面外剪切变形不能忽视。怎么处理？

🎓

这是"厚剪切变形"问题。通常的薄板理论（Kirchhoff-Love假设）忽略面外剪切，但复合材料的

$$ G_{13}, G_{23} $$

（面外剪切刚度）相对较小，不能忽视。有两个办法：第一，用考虑剪切变形的壳体单元（Abaqus的S4R、Ansys的SHELL181开启"Shear Correction"）。第二，若需更高精度，用连续体壳体单元或按厚度离散的实体单元，但计算成本会大幅增加。