与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 后处理结果