复合材料气瓶固化的热场本质上可以认为包含两个阶段，复合材料气瓶表面和周围通过与空气的热对流换热，复合材料与内部塑料芯模以及金属接头的传热。在第一个阶段整体温度较低，热量从表面向复合材料层内部流入，此时固化速率很低，所以自发热量小于从外面传进来的热量，热流由外向内。当整体温度达到一定温度时进入第二个阶段，自发热使复合材料层温度高于环境，复合材料层开始向外换热，因此温度场内高外低。如果趋势正确，那么决定温差的就是：

• 复合材料的热导率

• 自发热热流

一、材料性能

本节是一个热传导问题，需要给出的材料性质有：复合材料的密度、比热、热导系数和热产生（即内热源，通过用户子程序HETVAL定义）。由于复合材料是正交各项异性材料，必须定义材料坐标系。表征反应程度的固化度场作为用户定义场（user defined field）也需要在材料属性中定义，需要用到用户子程序USDFLD。主要材料性能如下表 1所示。

表 1 固化-温度场有限元分析材料参数

二、复合材料气瓶模型建立

采用abaqus-wcm建立好的复合材料模型如图1所示，复合材料层的变角度以及变厚度如图2所示。

图 1 复合材料气瓶建模图

图2  复合材料气瓶封头缠绕变角度与变厚度图

对建立模型的复合材料添加材料属性，如下图所示，根据缠绕定律，对于每个网格赋予不同的材料属性，实现复合材料缠绕层变角度与变厚度的设置，如图2所示。

图3 复合材料缠绕层变角度/变厚度材料属性赋予

三、子程序详解

为了计算复合材料固化过程中的温度以及固化度，必须用到4个abaqus子程序。

用户子程序USDFLD定义表征化学反应程度的固化度场，FILD（1）即是每个增量步中每个积分点上固化度的值。通过STATEV（1）与定义热产生的用户子程序HETVAL传递数据。由于本项目使用的描述固化反应的固化动力学方程中固化度初始值不能为零，所以在子程序USDFLD中给固化度赋予一个初始值1×10-4。详细子程序内容如下：

用户子程序DISP，用于定义温度边界条件：理解为固化温度的来源即是加热源，曲线为固化得工艺温度曲线，如图4 所示。

具体子程序内容如下所示：

图 4 复合材料气瓶固化曲线

用户子程序FILM用于定义对流换热边界条件，H（1）给定表面对流换热系数，SINK给定与图4的工艺曲线一致的环境温度。用户子程序的详细代码如下。

用户子程序HETVAL，用于定义复合材料内部的反应热。STATEV（1）即为从用户子程序USDFLD中传递过来的固化度；STATEV（2）为固化率（固化度的变化率），由与固化度、温度、升温速率等变量相关的固化动力学方程确定，决定了反应热的放热速率；FLUX（1）表示反应热产生的单位时间、单位体积的热量。子程序如下所示：

将以上四个子程序采用Fortran语言编程成abaqus子程序进行调用。

点击完整阅读全文 ：

复合材料专场 | ABAQUS车载四型复合材料气瓶固化过程的数值模拟分析攻略

相关阅读推荐：

复合材料结构有限元分析实例