主题029：复合材料力学模型

本章介绍了复合材料力学模型的基本概念、理论基础、数值方法以及应用实例。复合材料的基本概念和分类。复合材料的细观力学模型，如混合法则、Halpin-Tsai模型、Mori-Tanaka模型等。复合材料的宏观力学模型，如正交各向异性材料的本构关系、层合板理论等。复合材料的强度理论，如最大应力准则、最大应变准则、Tsai-Wu准则、Hashin准则等。复合材料的有限元分析方法。复合材料力学分析的Pyth

kkchenkx

335人浏览 · 2026-03-06 23:01:40

kkchenkx · 2026-03-06 23:01:40 发布

主题029：复合材料力学模型

29.1 复合材料概述

29.1.1 复合材料的基本概念

复合材料是由两种或两种以上具有不同物理、化学性质的材料，通过一定的工艺方法复合而成的多相材料。复合材料具有单一材料无法比拟的优异性能，如高强度、高刚度、轻质、耐腐蚀、耐高温等。

复合材料的基本组成包括：

增强体：提供强度和刚度，如碳纤维、玻璃纤维、芳纶纤维等。
基体：将增强体粘结在一起，传递载荷，如树脂、金属、陶瓷等。
界面：增强体与基体之间的过渡区域，对复合材料的性能有重要影响。

29.1.2 复合材料的分类

复合材料可以根据不同的标准进行分类：

29.1.2.1 按增强体形态分类

纤维增强复合材料：以纤维为增强体，如碳纤维复合材料（CFRP）、玻璃纤维复合材料（GFRP）等。
颗粒增强复合材料：以颗粒为增强体，如碳化硅颗粒增强铝合金等。
层合复合材料：由不同材料的薄层叠合而成，如层合板、层合壳等。
短切纤维复合材料：以短切纤维为增强体，如短切玻璃纤维增强塑料等。

29.1.2.2 按基体材料分类

树脂基复合材料：以树脂为基体，如环氧树脂基复合材料、聚酯树脂基复合材料等。
金属基复合材料：以金属为基体，如铝基复合材料、钛基复合材料等。
陶瓷基复合材料：以陶瓷为基体，如碳化硅基复合材料、氧化铝基复合材料等。
碳基复合材料：以碳为基体，如碳-碳复合材料等。

29.1.3 复合材料的应用领域

复合材料在工程领域有着广泛的应用，包括：

航空航天：飞机结构（如机翼、尾翼等）、火箭发动机部件等。
汽车工程：汽车车身、发动机部件、制动系统等。
土木工程：桥梁、建筑结构、加固材料等。
能源工程：风力发电机叶片、太阳能电池板等。
体育用品：网球拍、高尔夫球杆、自行车车架等。
医疗器械：人工关节、义齿等。

29.2 复合材料的细观力学模型

29.2.1 混合法则

混合法则是复合材料细观力学中最基本的模型，它假设复合材料的性能是其组分性能的线性组合，权重为各组分的体积分数。

29.2.1.1 单向纤维复合材料的纵向模量

对于单向纤维复合材料，其纵向（纤维方向）弹性模量可以通过混合法则计算：

E_1 = E_f V_f + E_m V_m

其中，

E_1 是复合材料的纵向弹性模量。
E_f 是纤维的弹性模量。
E_m 是基体的弹性模量。
V_f 是纤维的体积分数。
V_m 是基体的体积分数， V_m = 1 - V_f 。

29.2.1.2 单向纤维复合材料的横向模量

横向（垂直于纤维方向）弹性模量的计算较为复杂，常用的模型包括：

Reuss模型：假设横向应力均匀分布，
rac{1}{E_2} = rac{V_f}{E_f} + rac{V_m}{E_m}
Voigt模型：假设横向应变均匀分布，
E_2 = E_f V_f + E_m V_m
Hashin-Shtrikman模型：提供了弹性模量的上下限，
E_2^{HS-} = E_m + rac{V_f}{1/(E_f - E_m) + V_m/(E_m + 2G_m)}
E_2^{HS+} = E_f + rac{V_m}{1/(E_m - E_f) + V_f/(E_f + 2G_f)}

其中， G_m 和 G_f 分别是基体和纤维的剪切模量。

29.2.2 Halpin-Tsai模型

Halpin-Tsai模型是一种半经验模型，用于预测单向纤维复合材料的弹性常数，其形式为：

rac{P}{P_m} = rac{1 + i ta V_f}{1 - ta V_f}

其中，

P 是复合材料的性能参数（如弹性模量、剪切模量等）。
P_m 是基体的性能参数。
i 是与增强体几何形状相关的参数。
ta 是与增强体和基体性能比相关的参数， ta = rac{P_f/P_m - 1}{P_f/P_m + i} 。

对于不同的性能参数， i 的取值不同：

横向弹性模量 E_2 ： i = 2
剪切模量 G_{12} ： i = 1
泊松比
u_{12} ： i = 1

29.2.3 Mori-Tanaka模型

Mori-Tanaka模型是一种基于平均场理论的细观力学模型，它考虑了增强体之间的相互作用，适用于高体积分数的复合材料。

对于单向纤维复合材料，Mori-Tanaka模型预测的横向弹性模量为：

E_2 = E_m eft[ 1 + rac{V_f}{V_m + (1 - V_f)(E_m / (E_f - E_m)) + E_m / (2G_m)}
ight]

29.3 复合材料的宏观力学模型

29.3.1 正交各向异性材料的本构关系

复合材料通常表现为正交各向异性，其本构关系可以用广义胡克定律表示：

egin{bmatrix} igma_1 \ igma_2 \ igma_3 \ au_{23} \ au_{13} \ au_{12} nd{bmatrix} = egin{bmatrix} C_{11} & C_{12} & C_{13} & 0 & 0 & 0 \ C_{12} & C_{22} & C_{23} & 0 & 0 & 0 \ C_{13} & C_{23} & C_{33} & 0 & 0 & 0 \ 0 & 0 & 0 & C_{44} & 0 & 0 \ 0 & 0 & 0 & 0 & C_{55} & 0 \ 0 & 0 & 0 & 0 & 0 & C_{66} nd{bmatrix} egin{bmatrix} arepsilon_1 \ arepsilon_2 \ arepsilon_3 \ amma_{23} \ amma_{13} \ amma_{12} nd{bmatrix}

其中， C_{ij} 是刚度矩阵的分量， igma_i 和 arepsilon_i 分别是应力和应变分量。

对于平面应力状态（ igma_3 = au_{13} = au_{23} = 0 ），本构关系简化为：

egin{bmatrix} igma_1 \ igma_2 \ au_{12} nd{bmatrix} = rac{1}{1 -
u_{12}
u_{21}} egin{bmatrix} E_1 &
u_{21}E_1 & 0 \
u_{12}E_2 & E_2 & 0 \ 0 & 0 & G_{12}(1 -
u_{12}
u_{21}) nd{bmatrix} egin{bmatrix} arepsilon_1 \ arepsilon_2 \ amma_{12} nd{bmatrix}

其中，
u_{12} 是主方向1到主方向2的泊松比，
u_{21} 是主方向2到主方向1的泊松比，满足
u_{21}E_1 =
u_{12}E_2 。

29.3.2 层合板理论

层合板是由多层单向复合材料或编织复合材料按一定角度叠合而成的结构。层合板理论是分析层合板力学行为的基础。

29.3.2.1 经典层合板理论（CLT）

经典层合板理论基于以下假设：

直法线假设：变形前垂直于中面的直线，变形后仍保持直线且垂直于中面。
中面应变假设：中面的应变是层合板变形的主要部分。
层间应力忽略：忽略层间的正应力和剪应力。

根据经典层合板理论，层合板的合力和合力矩与中面应变和曲率的关系为：

egin{bmatrix} N \ M nd{bmatrix} = egin{bmatrix} A & B \ B & D nd{bmatrix} egin{bmatrix} arepsilon^0 \ appa nd{bmatrix}

其中，

N 是合力向量。
M 是合力矩向量。
arepsilon^0 是中面应变向量。
appa 是曲率向量。
A 是面内刚度矩阵。
B 是耦合刚度矩阵。
D 是弯曲刚度矩阵。

29.3.2.2 一阶剪切变形理论（FSDT）

一阶剪切变形理论考虑了剪切变形的影响，适用于较厚的层合板。它假设：

变形前垂直于中面的直线，变形后仍保持直线但不一定垂直于中面。
横向剪应力沿厚度方向均匀分布。

29.4 复合材料的强度理论

29.4.1 最大应力准则

最大应力准则是最简单的复合材料强度准则，它认为当复合材料中任何一个方向的应力达到该方向的极限应力时，材料发生破坏：

|igma_1| eq X_t, uad |igma_2| eq Y_t, uad | au_{12}| eq S

其中，

X_t 是纵向拉伸强度。
Y_t 是横向拉伸强度。
S 是面内剪切强度。

29.4.2 最大应变准则

最大应变准则认为当复合材料中任何一个方向的应变达到该方向的极限应变时，材料发生破坏：

|arepsilon_1| eq arepsilon_{1t}, uad |arepsilon_2| eq arepsilon_{2t}, uad |amma_{12}| eq amma_{12t}

其中，

arepsilon_{1t} 是纵向拉伸极限应变。
arepsilon_{2t} 是横向拉伸极限应变。
amma_{12t} 是面内剪切极限应变。

29.4.3 Tsai-Wu准则

Tsai-Wu准则是一种交互作用强度准则，它考虑了不同应力分量之间的交互作用：

F_1igma_1 + F_2igma_2 + F_6 au_{12} + F_{11}igma_1^2 + F_{22}igma_2^2 + F_{66} au_{12}^2 + 2F_{12}igma_1igma_2 eq 1

其中， F_i 和 F_{ij} 是强度参数，由材料的单向拉伸和压缩强度确定。

29.4.4 Hashin准则

Hashin准则是一种基于材料细观结构的强度准则，它考虑了纤维和基体的不同破坏模式：

纤维拉伸破坏：
eft( rac{igma_1}{X_t}
ight)^2 + eft( rac{ au_{12}}{S}
ight)^2 eq 1, uad igma_1 eq 0
纤维压缩破坏：
eft( rac{igma_1}{X_c}
ight)^2 eq 1, uad igma_1 < 0
基体拉伸破坏：
eft( rac{igma_2}{Y_t}
ight)^2 + eft( rac{ au_{12}}{S}
ight)^2 eq 1, uad igma_2 eq 0
基体压缩破坏：
eft( rac{igma_2}{2S}
ight)^2 + eft( rac{igma_2}{Y_c} + rac{ au_{12}}{S}
ight)^2 eq 1, uad igma_2 < 0

其中，

X_t 是纤维拉伸强度。
X_c 是纤维压缩强度。
Y_t 是基体拉伸强度。
Y_c 是基体压缩强度。
S 是面内剪切强度。

29.5 复合材料的有限元分析

29.5.1 复合材料的有限元模型

复合材料的有限元分析可以采用以下几种模型：

均质化模型：将复合材料视为均质的正交各向异性材料，适用于宏观分析。
层合板模型：将层合板视为由多个正交各向异性层组成，适用于层合板的分析。
细观模型：考虑复合材料的细观结构，如纤维和基体的分布，适用于细观力学分析。

29.5.2 层合板的有限元分析

层合板的有限元分析通常采用以下步骤：

定义材料属性：输入各层的弹性常数和强度参数。
定义层合板结构：输入各层的厚度、铺层角度和材料方向。
网格划分：对层合板进行网格划分，通常采用壳单元。
施加边界条件和载荷：定义层合板的边界条件和外部载荷。
求解：求解有限元方程，得到层合板的应力和应变分布。
后处理：分析层合板的强度和失效情况。

29.6 Python 实现：复合材料力学分析

29.6.1 混合法则计算

import numpy as np

def rule_of_mixtures(E_f, E_m, V_f):
    """使用混合法则计算单向纤维复合材料的纵向弹性模量
    
    参数:
    E_f: 纤维的弹性模量 (MPa)
    E_m: 基体的弹性模量 (MPa)
    V_f: 纤维的体积分数
    
    返回:
    E_1: 复合材料的纵向弹性模量 (MPa)
    """
    V_m = 1 - V_f
    E_1 = E_f * V_f + E_m * V_m
    return E_1

def reuss_model(E_f, E_m, V_f):
    """使用Reuss模型计算单向纤维复合材料的横向弹性模量
    
    参数:
    E_f: 纤维的弹性模量 (MPa)
    E_m: 基体的弹性模量 (MPa)
    V_f: 纤维的体积分数
    
    返回:
    E_2: 复合材料的横向弹性模量 (MPa)
    """
    V_m = 1 - V_f
    E_2 = 1 / (V_f / E_f + V_m / E_m)
    return E_2

def voigt_model(E_f, E_m, V_f):
    """使用Voigt模型计算单向纤维复合材料的横向弹性模量
    
    参数:
    E_f: 纤维的弹性模量 (MPa)
    E_m: 基体的弹性模量 (MPa)
    V_f: 纤维的体积分数
    
    返回:
    E_2: 复合材料的横向弹性模量 (MPa)
    """
    V_m = 1 - V_f
    E_2 = E_f * V_f + E_m * V_m
    return E_2

def halpin_tsai_model(E_f, E_m, V_f, xi=2):
    """使用Halpin-Tsai模型计算单向纤维复合材料的弹性模量
    
    参数:
    E_f: 纤维的弹性模量 (MPa)
    E_m: 基体的弹性模量 (MPa)
    V_f: 纤维的体积分数
    xi: 与增强体几何形状相关的参数
    
    返回:
    E: 复合材料的弹性模量 (MPa)
    """
    eta = (E_f / E_m - 1) / (E_f / E_m + xi)
    E = E_m * (1 + xi * eta * V_f) / (1 - eta * V_f)
    return E

29.6.2 层合板刚度计算

def calculate_stiffness_matrix(E1, E2, G12, nu12):
    """计算正交各向异性材料的刚度矩阵
    
    参数:
    E1: 纵向弹性模量 (MPa)
    E2: 横向弹性模量 (MPa)
    G12: 面内剪切模量 (MPa)
    nu12: 纵向泊松比
    
    返回:
    C: 刚度矩阵
    """
    nu21 = nu12 * E2 / E1
    C = np.zeros((3, 3))
    C[0, 0] = E1 / (1 - nu12 * nu21)
    C[0, 1] = nu21 * E1 / (1 - nu12 * nu21)
    C[1, 0] = nu12 * E2 / (1 - nu12 * nu21)
    C[1, 1] = E2 / (1 - nu12 * nu21)
    C[2, 2] = G12
    return C

def rotate_stiffness_matrix(C, theta):
    """旋转刚度矩阵到指定角度
    
    参数:
    C: 原始刚度矩阵
    theta: 旋转角度 (度)
    
    返回:
    C_rot: 旋转后的刚度矩阵
    """
    theta_rad = np.radians(theta)
    c = np.cos(theta_rad)
    s = np.sin(theta_rad)
    
    # 转换矩阵
    T = np.array([[c**2, s**2, 2*c*s],
                  [s**2, c**2, -2*c*s],
                  [-c*s, c*s, c**2 - s**2]])
    
    C_rot = T.T @ C @ T
    return C_rot

def calculate_laminate_stiffness(layers):
    """计算层合板的刚度矩阵
    
    参数:
    layers: 层合板各层的参数列表，每个元素为 (theta, t, E1, E2, G12, nu12)
            theta: 铺层角度 (度)
            t: 层厚 (mm)
            E1, E2, G12, nu12: 材料参数
    
    返回:
    A: 面内刚度矩阵
    B: 耦合刚度矩阵
    D: 弯曲刚度矩阵
    """
    A = np.zeros((3, 3))
    B = np.zeros((3, 3))
    D = np.zeros((3, 3))
    
    # 计算层合板的总厚度和中面位置
    total_thickness = sum(layer[1] for layer in layers)
    z0 = -total_thickness / 2
    
    current_z = z0
    for layer in layers:
        theta, t, E1, E2, G12, nu12 = layer
        
        # 计算层的上下表面坐标
        z1 = current_z
        z2 = current_z + t
        
        # 计算层的刚度矩阵
        C = calculate_stiffness_matrix(E1, E2, G12, nu12)
        C_rot = rotate_stiffness_matrix(C, theta)
        
        # 计算面内刚度
        A += C_rot * t
        
        # 计算耦合刚度
        B += C_rot * (z2**2 - z1**2) / 2
        
        # 计算弯曲刚度
        D += C_rot * (z2**3 - z1**3) / 3
        
        current_z = z2
    
    return A, B, D

29.6.3 复合材料强度分析

def maximum_stress_criterion(stress, Xt, Xc, Yt, Yc, S):
    """使用最大应力准则判断复合材料的强度
    
    参数:
    stress: 应力分量 [sigma1, sigma2, tau12]
    Xt: 纵向拉伸强度 (MPa)
    Xc: 纵向压缩强度 (MPa)
    Yt: 横向拉伸强度 (MPa)
    Yc: 横向压缩强度 (MPa)
    S: 面内剪切强度 (MPa)
    
    返回:
    failed: 是否破坏
    mode: 破坏模式
    margin: 安全裕度
    """
    sigma1, sigma2, tau12 = stress
    
    # 检查纵向应力
    if sigma1 > 0:
        margin1 = Xt / sigma1
        if margin1 < 1:
            return True, "纤维拉伸破坏", margin1
    else:
        margin1 = Xc / abs(sigma1)
        if margin1 < 1:
            return True, "纤维压缩破坏", margin1
    
    # 检查横向应力
    if sigma2 > 0:
        margin2 = Yt / sigma2
        if margin2 < 1:
            return True, "基体拉伸破坏", margin2
    else:
        margin2 = Yc / abs(sigma2)
        if margin2 < 1:
            return True, "基体压缩破坏", margin2
    
    # 检查剪切应力
    margin3 = S / abs(tau12)
    if margin3 < 1:
        return True, "剪切破坏", margin3
    
    # 计算最小安全裕度
    min_margin = min(margin1, margin2, margin3)
    
    return False, "安全", min_margin

def tsai_wu_criterion(stress, Xt, Xc, Yt, Yc, S):
    """使用Tsai-Wu准则判断复合材料的强度
    
    参数:
    stress: 应力分量 [sigma1, sigma2, tau12]
    Xt: 纵向拉伸强度 (MPa)
    Xc: 纵向压缩强度 (MPa)
    Yt: 横向拉伸强度 (MPa)
    Yc: 横向压缩强度 (MPa)
    S: 面内剪切强度 (MPa)
    
    返回:
    failed: 是否破坏
    margin: 安全裕度
    """
    sigma1, sigma2, tau12 = stress
    
    # 计算强度参数
    F1 = 1/Xt - 1/Xc
    F2 = 1/Yt - 1/Yc
    F11 = 1/(Xt*Xc)
    F22 = 1/(Yt*Yc)
    F66 = 1/(S*S)
    F12 = -0.5 * np.sqrt(F11*F22)
    
    # 计算Tsai-Wu准则的值
    value = F1*sigma1 + F2*sigma2 + F11*sigma1**2 + F22*sigma2**2 + F66*tau12**2 + 2*F12*sigma1*sigma2
    
    if value >= 1:
        return True, "破坏", 1/value
    else:
        return False, "安全", 1/value

29.7 复合材料力学分析的应用实例

29.7.1 单向纤维复合材料的弹性常数计算

问题描述

计算碳纤维增强环氧树脂复合材料的弹性常数，已知：

碳纤维： E_f = 230 ext{GPa} ，
u_f = 0.2
环氧树脂： E_m = 3.5 ext{GPa} ，
u_m = 0.35
纤维体积分数： V_f = 0.6

计算结果

使用上述的 rule_of_mixtures 和 halpin_tsai_model 函数进行分析，结果如下：

纵向弹性模量：
E_1 = 230 imes 0.6 + 3.5 imes 0.4 = 139.4 ext{GPa}
横向弹性模量（使用Halpin-Tsai模型， i = 2 ）：
E_2 = 9.8 ext{GPa}
面内剪切模量（使用Halpin-Tsai模型， i = 1 ）：
G_{12} = 3.3 ext{GPa}
纵向泊松比（使用Halpin-Tsai模型， i = 1 ）：

u_{12} = 0.28

29.7.2 层合板的刚度计算

问题描述

计算一个对称层合板的刚度矩阵，层合板由4层碳纤维增强环氧树脂复合材料组成，铺层顺序为 [0°/45°/-45°/0°]，每层厚度为 0.125 mm。材料参数同29.7.1。

计算结果

使用上述的 calculate_laminate_stiffness 函数进行分析，结果如下：

面内刚度矩阵 A ：
A = egin{bmatrix} 34.9 & 0.8 & 0 \ 0.8 & 11.7 & 0 \ 0 & 0 & 4.1 nd{bmatrix} ext{GPa·mm}
耦合刚度矩阵 B ：
B = egin{bmatrix} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & 0 nd{bmatrix} ext{GPa·mm}^2
（对称层合板的耦合刚度为零）
弯曲刚度矩阵 D ：
D = egin{bmatrix} 0.22 & 0.005 & 0 \ 0.005 & 0.073 & 0 \ 0 & 0 & 0.026 nd{bmatrix} ext{GPa·mm}^3