与Fluent 动网格+高铁制动盘制动过程仿真(一)相比,增加了一个模块,是用来划分固体域网格。
2.1 导入几何
与Fluent 动网格+高铁制动盘制动过程仿真(一)一致,因此不做过多阐述:





固体域区域需要注意,各部分命名如下图:



2.2 网格设置




3.1 General设置与网格导入
首先导入网格,由于是三部分网格,因此需要通过附加case的方式,将其余两部分网格导入,然后勾选稳态计算,具体设置如下图所示。


3.2 材料设置
此处对材料进行设置,采用air作为流体计算材料,具体设置如下图所示:



DEFINE_CG_MOTION(lc, dt, vel, omega, time, dtime){real t=CURRENT_TIME;omega[0]=0;omega[1]=0;omega[2]=-1*(100-t)/0.4;vel[0]=(100-t);vel[1]=0;vel[2]=0;}DEFINE_CG_MOTION(lc1, dt, vel, omega, time, dtime){real t=CURRENT_TIME;vel[0]=(100-t);vel[1]=0;vel[2]=0;}DEFINE_PROFILE(heatersource, thread, position){face_t f;real x[3];real eff = 0.85;real u = 0.3;real F = 20000;real s = 0.02;real w;real v = 100;real a = -1;w = 60 * 3.14 / 180.0;real tan_w = tan(w);begin_f_loop(f, thread){F_CENTROID(x, f, thread);real t = CURRENT_TIME;real temp_x0 = x[0] - (v * t + 0.5 * a * t * t);real r = sqrt(temp_x0 * temp_x0 + x[1] * x[1]);if (r > 0.15 && r < 0.3) {;}else{r=0;}real A = temp_x0 * tan_w - x[1];if (A >= 0.0 && x[1] >= 0.0) {F_PROFILE( f, thread, position) =(eff * F * u * (v + a * t) / 0.4) * (r / s);} else {F_PROFILE( f, thread,position) = 0.0;}}end_f_loop(f, thread);}
由于是多个域组成,因此需要通过interface将各个区域连接起来,具体设置如下图:

本案例最重要的便是动网格的设置,增加了固体域的动网格设置,具体设置如下图:


3.7 初始化设置
首先进行标准初始化设置,具体设置如下图:

3.8 计算设置
此处进行的计算设置如下:

4.1 后处理结果
