rtm树脂流动模拟及其对模具设计的意义.docVIP

RTM树脂流动模拟及其对模具设计的意义 　　在复合材料加工中，模具的设计和制造对整个生产过程具有决定性影响。要设计和制造出合理的模具，仅仅依赖经验是不够的，因此国内外都开展了RTM工艺的数值模拟技术研究，利用数值模拟对模具设计方案检验和优化。近年来国外已推出了一些模拟软件，其中比较著名的软件有密西根州立大学的LMPC、福特公司开发的CFILL等[1]。在国内，孙玉敏[2]等利用RTM模拟仿真系统对风机叶片的RTM充模过程进行了三维模拟, 以优化成型工艺方案得到最佳的注射位置、合理的溢料口位置以及选择溢料口的依据，从而设计出合理的模具；邓京兰等[3-8]详细探讨了RTM成型工艺中树脂流动的充模过程, 利用多孔介质流体动力学基本理论和Darcy方程对RTM充模过程进行深入的数值分析,建立模拟RTM成型工艺中树脂等温填充过程的数学模型；王秀君[9]等针对三维异型零件，采用贴体坐标法与有限差分相结合的方法，研究RTM注射工艺的过程；Junying Yang[10]等提出了在纤维增强体束内采用Brinkman方程描述流动，在纤维束间采用Stokes方程进行描述，从而说明气泡的产生是由于微观和宏观树脂流动速度不一致导致两者的流动前锋不平行而形成的。上海交通大学塑性成形系[11]也根据流体动力学的基本理论和Darcy方程编写代码模拟非等温过程，取得一定进展。国外的相关研究更为广泛和深入，Akbar Shojaei[12-22]等利用Darcy方程、守恒方程和能量公式以及固化方程建立数学模型进行模拟，计算出树脂前锋位置以及压力场、温度场等，并与试验结果进行比较，吻合得较好。D.Antonelli[23]等根据固－液混合理论，建立多孔介质流体流动的温度耦合模型，从而能更精确地描述树脂在填充时的流动模式。通过这些数值模拟可以发现，纤维的渗透率、树脂的黏度和注射压力（或者注射流率）对制定合理的工艺和设计模具是非常重要的。本文运用非结构化四面体网格，结合控制体/有限元（CV/FEM）方法，根据流动控制方程，编写计算代码，对RTM工艺进行三维实体模拟，探讨其对模具设计的影响。 　　流动控制方程 　　为便于说明问题，采用三维等温树脂流动模型。 　　在RTM工艺中，流体的宏观流动过程可以用Darcy方程很好地描述。Darcy方程可写成以下形式[1]： 　　式中，v是流体的Darcy速度，是流体的压力梯度，[k]是增强体的渗透率张量，μ是树脂黏度。RTM工艺的流动过程为一准静态过程，质量守恒方程可写成如下形式： 　　 　　这样，方程（1）代入方程（2）得到流动控制方程： 　　其边界条件是： 　　（1）在注射口处，常压注射时P =P0，常流率注射时V=V0； 　　（2）在流动前沿处，P=0； 　　（3）在模具壁处，。 　　数值模拟 　　对于控制方程的离散，一般采用控制体/有限元（CV/FEM）的方法进行。这种方法最大的优点在于计算过程中无需网格的重新划分，大大节约了计算时间并降低了计算误差。在追踪流动前沿的过程中，用参数fi表示控制体是否被填充或者被部分填充。当fi=1时，表示控制体被完全填充；当 fi=0时，控制体完全没有被填充；当1fi0时，表示控制体正在被填充，此处为流体流动前沿。 　　通过填充系数，流动前锋位置就可以确定下来，控制体的填充的状态就可以确定下来。填充系数fi由下式确定： 　　式中，Qi(t)为流入该控制体的总的流率，△tmin为在所有即将被填充的控制体中花费最少的时间，fi0为前一次的填充系数，Vi为控制体的体积。 　　数值计算过程按以下步骤进行： 　　（1）模拟开始，计算单元参数并生成控制体模型与相关控制体参数； 　　（2）充模开始时，假定注射孔所在的控制体首先一次性充满，而在填充过程中，控制体被填充的程度是根据填充系数来确定的； 　　（3）根据控制方程以及边界条件，计算节点压力； 　　（4）根据Darcy方程确定流体运动速度，从而确定下一个时间步长，计算树脂流进每一个流动前锋控制体的体积，根据公式（4）更新相应的控制体填充系数； 　　（5）对新的被填充区域重复步骤（3）、（4），直至所有的控制体填充系数都为1.0，填充过程结束。 　　计算案例 　　根据流动控制方程和计算步骤，编写了计算程序。 　　在本案例中，使用的材料为石墨纤维/环氧树脂,模腔形状如图1所示，由于四面体网格对制品形状适应性很强，体现流体的三维实体流特点，故本例中划分了474个四面体单元。 　　案例1:各向同性的纤维渗透率张量kx=ky=kz=1.2×10-8m-2,树脂粘度μ=0.5Pa·s，恒压注射，注射压力P=0.1MPa。在流动模式中会发现流动不是按照制品中心对称的，需要注意。因此为便于观察树脂的流动模式，注射孔坐标为（44.57,22.67