Particleworks案例5-油箱晃荡

问题描述模拟泡罩的典型热成型情况，热成型泡罩、片材（蓝色）和模具（红色）如下。为了减少计算运行时间，并利用泡罩的对称特性，只有四分之一参与仿真计算。初始（1/4）薄膜具有以下尺寸：• 长度 = 15 mm• 宽度 = 5 mm• 初始厚度 = 0.35 mm 导入网格启动软件，选项设置为F.E.M. task、Time-dependent problem(s)、2D shell geometry。创建新mold，点击Adiabatic mold。计算域和边界选择Subdomain 1，点击Remove。设置接触条件，选择No contact along Subdomain 2和Contact。定义Mold运动。点击Mold motion type : fixed mold。点击菜单EVOL，设置mold的瞬态速度，选择velocity-x.同样操作，设置y-velocity，值为0。设置z-velocity，值为10 mm/s，时间关系为f(t) = Ramp function，给初值a为0.097，b、c、d的初值为1、0、0.103。创建sub-task，选择Shell model : Gen. Newtonian isothermal。选择Subdomain 2，点击Remove。设置Flow boundary conditions。选择Zero wall velocity (vn=vs=0) along Boundary 3边界，点击Plane of symmetry ( fs=0, vn=0 )，选择normal direction along X axis。选择Zero wall velocity (vn=vs=0) along Boundary 4边界，点击Plane of symmetry ( fs=0, vn=0 )，选择normal direction along Y axis。点击Inflation pressure imposed按钮，点击Constant，输入100000 Pa。点击菜单EVOL，打开Time dependence of inflation pressure面板，选择f(t) = Ramp function，a、b、c、d值分别为0.1、0、0.11、1。定义blister的接触。在Define contacts中点击Create a new contact problem。点击Select a contact wall，选择Mold 1 : Contact along Subdomain 2。点击Modify slipping coefficient，设置为1e10，点击Modify penalty coefficient，设置为1e10，点击Modify penetration accuracy，设置为0.05。定义mold的坐标。点击Specify mold side / cavity side，点击No。定义layers。点击创建Create a new layer，点击Inertia terms，选择Inertia will be taken into account。点击Modify coef. of permeability，设置5e-12。点击设置Initial thickness，点击Constant，设置0.35。求解设置设置Modify the min value of the time-step为0.0001s。Modify the max value of the time-step为0.01s。仿真结果显示厚度分布云图显示不同时刻下的厚度。来源：CFD饭圈