传热与固化过程的耦合作用在复合材料成型工艺中具有重要影响。本文基于非结构化网格有限体积法提出一种解决结构形式复杂的复合材料成型过程中各向异性非稳态导热与固化反应动力学耦合问题的算法。各向异性复合材料的导热控制方程空间项采用非结构化网格有限体积法离散,时间项采用欧拉隐格式求解。固化动力学方程采用了欧拉显格式计算,并通过固化反应速率计算出能量方程的热源,同时根据温度分布及其变化更新固化反应速率,实现导热与固化反应的双向耦合。最后,通过文献中的几个经典算例验证了本文方法的正确性和可靠性,同时也表明本文方法可以处理复杂结构形式的复合材料成型过程的传热与化学反应的耦合模拟问题。
Thermal and curing reaction coupling has great influence on composites molding process.An algorithm of solving complicated structure composites curing process with coupled anisotropic transient thermal and curing reaction kinetics formulation based on unstructured gird finite volume method is presented.The spatial item of heat diffusion equation of anisotropic composites is discretized by unstructured grid finite volume method,and time item is discretized by Euler implicit scheme.Curing kinetics equation is calculated by Euler explicit scheme,and obtained the heat source of the energy equation through the curing reaction rate which subjected to temperature distribution and corresponding manufacturing process plans,thus realizing the two-way coupling of thermal and curing reaction.Finally,this method is verified correctness and reliability through several classical examples for which prior results and experimental data are available in the literature,at the same time,it shows that the method presented can deal with the numerical simulation of the coupled thermal and curing reaction kinetics of composites molding process with complicated structure.