镁合金壳体零件砂型铸造工艺分析及数值模拟.doc

镁合金壳体零件砂型铸造工艺分析及数值模拟 1、相关定义 1.1、VOF法的基本概念 为了描述充型合金液的自由表面,得出自由表面的形状和位置信息,对自 由表面就需要进行一定程度的处理。由于VOF法其具有运算速度快、占用空间 小等特点,在众多处理自由表面的方法中,应用范围最广且较适宜于处理自由 表面。 用VOF法处理自由表面时,需要通过对体积函数方程进行求解来获得得到 自由表面的状态信息。体积函数方程如下: 37 武汉理工大学硕士学位论文 . ,+塑+壁+ ! = 0’ (4-18) dx dy dz dt 用体积函数F(x, y , z , t)来描述整个流动场的自由表面,F(x, y, z, t)在点(X , y , z)处定义为: 流体的体积F(x+)、 单元内= (l-FD)dx (3) 如图4-5 (c),速度u的方向与自由表面垂直,且间隙在左侧,则此 时流体可能全部进入右边的网格单元,也可能只有部分进入: . uuxdt拳dF ={xdt (4.31”[) Fj3 X dx uxdtFDdx 对于u〈0的情况同理也可做上述的处理。但是用上述方法所确定的F值可 能出现比0小或者比1大的情况,因此,在完成流量计算后,必须把F值小于0 的网格单元重设为0,并对在计算中因这些调整所产生的流体体积积累的变化量 做记录。那么就需要对F值作一些调整,使自由表面单元的F值在0和1之间 变化。当在数值计算中,为提高程序的可读性即运行效率,设置一个判断因数 e,当单元的F值小于e时,即认为F为0,当F值比1-e大时,即认为网格 单元为满单元体。 ■ 对于三维情况做以下的处理,对体积函数方程在一个自由表面网格单元上 进行积分: 武汉理工大学硕士学位论文 F*。’p-Fp FF的压力;Ps为自由表面的压力;dc为两相 邻网格中心距离;d为自由表面网格中心与流体内部中心的距离。但是实际的自 由表面往往会出现相互交叉或相互重叠的现象,所以有可能同时出现多个差值 计算,为解决这个问题,我们一般用以下方法来计算自由表面单元中心的平均 压力?? Wn=l (4-41) 44 武汉理工大学硕士学位论文 式中:Pu心压力,k_表面单元中;F——p表面单元的流体体积分数;N ——相邻满单元的个数;P——d相邻满单元的压力;Ps——背压。 b.型壁 在有一个控制面与型壁相重合时,根据动量离散方程,此时压力项就不起 作用了,为简化计算就可以将相邻网格的压力值和边界网格的压力值设置成相 等,即- Pi+1’j’k =Pi’j’k (4-42) 4. 3. 2速度边界条件 a.型壁速度边界条件 如图4-7所示,为了求出阴影控制面上的速度,必须在此类边界的铸型一 侧设置假想网格使其能够参与迭代计算。在左边(i, j, k)为边界网格,(i+l, j, k)为添加的假想型壁网格,在假想网格上设置速度边界条件。一般存在两 种边界条件:无滑动边界以及自由滑动边界[39]。 mt* fs ‘1;. 1/ f V 图4-7 型壁速度边界条件 自由滑动型壁边界条件: Vwj43) .k =Vij.k (4-无滑动型壁边界条件: 45 武汉理工大学硕士学位论文 Vw 4) ,j,k =—Vi,j,k(4-4为了使计算结果更接近于实际边界条件的情况,对于上述的边界情况亦可 采取如下的修正方法进行修正: Uw,j,k=0 v=w,j,k9v,,j,k (4-45) Ww=,j,k 0Wi j k 式中e为修正系数,当9=-1时,那么此时就为无滑动边界条件;当6=1 时,为自由滑动边界条件;而-1〈0 , int mElt —emen; typedef struct ELEMENT { 53 武汉理工大学硕士学位论文 double vexO[3]; //网格坐标 char Gate; //标记浇口 t,不属于f char Mark; //标记铸件t,不属于f char Surface; //表面网格t,不是f double P; //网格单元的压力 double dp; //校正压力变化值 double Vof; //网格的新体积函数值 double Vofb; 〃网格当前的体积函数值 double dFx,dFy,dFz; //网格右,前,上控制面的流量改变率,流出为 负,流进为正 double u,v,w; 〃x,y,z方向上的速度,n+1时亥ij double ub,vb,wb; //x,y,z方向的上一时间步的速度,即n时刻 double density; 〃密度 int direction; //自由表面的位置 double D; //连续性方程值 } CElement; 使用方法如下: typedef UltraCArrayTEl_ Element 1