基于P1/P0四面体宏元的三维注塑充填速度场压力场有限元数值模拟  

3D numerical simulation with FEM of P1/P0 tetrahedral macro-element for velocity and pressure of plastic injection molding in filling stage

作　　者：刘永志[1] 赵振峰[1] 申长雨[1] 

机构地区：[1]郑州大学橡塑模具国家工程研究中心,河南郑州450002

出　　处：《化工学报》2011年第4期1124-1129,共6页CIESC Journal

基　　金：国家科技支撑计划项目(2006BAF04B11);长江学者和创新团队发展计划项目~~

摘　　要：计算效率和解的稳定性是影响三维注塑充填有限元数值模拟的关键因素。针对黏性不可压缩聚合物熔体三维充填过程的速度场和压力场,以提高计算机求解速度为出发点,分析了采用P1/P0四面体单元(速度线性,压力常数)得到的有限元方程的解不收敛的原因,提出一种采用P1/P0四面体宏元离散空间域的求解方案,从而降低了求解的自由度数量,提高了计算速度。通过模拟'圆管中定黏度流体的稳态流动'考察了该求解方案的模拟精度,通过模拟'圆管中定黏度流体的瞬态充填'比较了采用P1/P0四面体宏元和P2/P1四面体单元(二次速度,线性压力)的计算时间。最终将该方案应用到三维注塑充填的数值模拟中。Efficiency and stability are the key factors affecting the simulation results of three dimensional（3D）plastic injection molding with the finite element method（FEM）.In this paper tetrahedral element based on the Galerkin method was used for the 3D flow simulation of viscous incompressible plastic melt.As the scheme of P1/P0 tetrahedral element（linear velocity and constant pressure）produced unstable results,a new P1/P0 tetrahedral macro-element was used for the discretization of the spatial domain,which was stable and could reduce the CPU time by decreasing the degrees of freedom of the equation.As a numerical test,developed flow with constant viscosity in the cylindrical geometry was used to compare the accuracy of the results given by P1/P0 tetrahedral macro-element and P2/P1 tetrahedral element（quadratic velocity and linear pressure）.Besides,the higher efficiency of the P1/P0 tetrahedral macro-element was obviously witnessed by another test of simulating the filling process of the fluid with constant viscosity in cylindrical geometry.In the last part of this paper the scheme of P1/P0 tetrahedral macro-element was used for the simulation of plastic injection molding in filling stage.

关 键 词：注塑 充填 P1/P0四面体宏元 三维数值模拟 

分 类 号：TQ320.66[化学工程—合成树脂塑料工业]