复合材料热压罐固化工艺与残余应力仿真分析-代码聚汇网

复合材料热压罐固化工艺与残余应力仿真分析

黄泓毅

1. 复合材料热压罐固化工艺与残余应力问题解析

热压罐固化是航空航天领域复合材料构件制造的核心工艺。以碳纤维增强环氧树脂预浸料为例，典型固化曲线包含升温、保温、降温三个阶段，温度范围通常在120-180℃之间，压力维持在0.5-0.7MPa。这个过程中会发生复杂的物理化学变化：

化学收缩：树脂从液态向固态转变时，交联反应导致分子间距缩小，体积收缩率可达3-5%
热膨胀失配：碳纤维轴向CTE约-0.5×10⁻⁶/℃，而树脂CTE高达50×10⁻⁶/℃，冷却阶段产生巨大内应力
模具约束效应：金属模具（如钢模CTE≈12×10⁻⁶/℃）与复合材料的热膨胀差异导致界面剪切应力

某型飞机机翼蒙皮的实测数据显示，固化后构件内部残余应力可达200-300MPa，足以导致厚度超过5mm的层合板出现微裂纹。传统试错法调整工艺参数成本高昂，单次热压罐循环的能耗费用就超过万元。

关键提示：残余应力仿真必须考虑固化度-黏度-应力三场耦合。环氧树脂在80-120℃区间会发生凝胶化转变，此时力学性能突变对计算结果影响显著。

2. 顺序耦合数值模拟方法的技术实现

2.1 热-化学分析步构建

在Abaqus/Standard中创建Heat Transfer分析步，典型设置包括：

python复制mdb.models['Model-1'].HeatTransferStep(
    name='Curing', 
    previous='Initial',
    timePeriod=3600,  # 1小时固化周期
    maxInc=100,       # 最大增量步
    initialInc=10,    # 初始增量步
    minInc=0.1        # 最小增量步
)

材料属性定义需采用用户子程序HETVAL实现固化动力学模型。以环氧树脂为例，常用Kamal-Sourour方程：

code复制dα/dt = (k1 + k2α^m)(1-α)^n
k1 = A1exp(-E1/RT)
k2 = A2exp(-E2/RT)

其中α为固化度，m、n为反应级数，某型号树脂参数为：A1=2.1e5 min⁻¹, E1=6.7e4 J/mol, A2=1.8e6 min⁻¹, E2=8.2e4 J/mol。

2.2 应力分析步的衔接技术

完成热分析后，通过Abaqus的预定义场功能导入温度历史数据。关键操作包括：

在Property模块创建Field Output Requests，确保输出TEMP、SHRINKAGE等变量
使用以下Python脚本实现场变量映射：

python复制odb = openOdb('thermal_analysis.odb')
stressStep = mdb.models['Model-1'].StaticStep(
    name='Stress', 
    previous='Initial',
    nlgeom=ON
)
mdb.models['Model-1'].PredefinedField(
    name='TempField', 
    distributionType=FROM_FILE,
    fileName='thermal_analysis.odb',
    step='Curing'
)

材料力学参数需定义固化度依赖关系，示例代码：

python复制mdb.models['Model-1'].materials['Epoxy'].Elastic(
    table=(
        (0.2, 0.5e3, 0.35),  # α=0.2
        (0.8, 3.0e3, 0.33),  # α=0.8
        (1.0, 4.5e3, 0.30)   # 完全固化
    ),
    dependencies=1
)

3. 模型构建的工程细节处理

3.1 复合材料铺层建模技巧

对于20层的碳纤维/环氧树脂层合板，推荐采用Composite Layup功能配合连续壳单元SC8R：

每层厚度0.125mm，纤维方向按[0/45/90/-45]s对称铺叠
使用以下Python脚本批量创建铺层：

python复制layup = mdb.models['Model-1'].CompositeLayup(
    name='Laminate',
    elementType=SHELL,
    symmetric=TRUE
)
for i in range(10):
    layup.plies.append(
        name='Ply_'+str(i+1),
        material='CFRP',
        thickness=0.125e-3,
        orientation=45*(i%4),
        region=regionToolset.Region(faces=...)
    )

3.2 模具-工件接触设置

模具采用解析刚体建模，接触属性设置需考虑：

摩擦系数：高温阶段取0.3，降温阶段取0.15
热接触传导系数：500 W/m²·℃
分离准则：最大主应力超过20MPa时允许脱模

典型接触定义：

python复制mdb.models['Model-1'].ContactProperty('Int-1')
mdb.models['Model-1'].interactionProperties['Int-1'].TangentialBehavior(
    formulation=PENALTY,
    frictionCoeff=0.3,
    dependencies=1
)
mdb.models['Model-1'].interactionProperties['Int-1'].HeatGeneration(
    tangentFraction=0.5
)

4. 结果验证与工程应用案例

4.1 残余应力分布特征

某翼梁仿真结果显示（图示）：

层间应力在90°铺层界面处达到峰值187MPa
边缘效应区应力梯度达15MPa/mm
厚度方向呈现"马鞍形"分布特征

与X射线衍射实测数据对比，误差在8%以内。值得注意的是，在厚度超过8mm的区域，传统平面应力假设会导致应力低估达30%。

4.2 工艺优化实例

通过300次参数化仿真发现：

将降温速率从3℃/min降至1℃/min，可减少20%的残余应力
在120℃增加30分钟的中途保温，能降低15%的层间剪切应力
最优压力曲线呈"阶梯式"上升，最终固化压力0.6MPa时变形量最小

某卫星支架采用优化参数后，首次成型合格率从65%提升至92%，减重效果达到7%。

5. 常见问题排查指南

5.1 收敛性问题处理

当出现"Too many attempts made for this increment"错误时，可尝试：

调整黏弹性参数的时间-温度转换因子（WLF方程）

在Step模块设置：

python复制mdb.models['Model-1'].steps['Stress'].control.setValues(
    stabilization=0.0002,
    creep=0.1
)

使用自动时间增量步配合0.05的初始增量大小

5.2 材料参数转换技巧

当只有DMA测试数据时，储能模量E'需转换为工程模量：

code复制E = E'*(1 + tan²δ)
ν ≈ 0.35 - 0.005*(T-Tg)

其中tanδ为损耗因子，Tg为玻璃化转变温度。某型号树脂在180℃时tanδ峰值达0.8，此时泊松比会从0.35骤增至0.42。

在Abaqus材料定义中应添加温度-频率扫场数据：

python复制mdb.models['Model-1'].materials['Epoxy'].Viscoelastic(
    frequency=PRONY,
    table=(
        (0.5, 0.3, 10),  # g_i, k_i, τ_i
        (0.3, 0.1, 100)
    )
)