发一个例子,一缸水收到冲击:
代码如下:
/prep7
!!!!! Define some dimensions and gravity
LX = 5.0 ! X length
LY = 2.5 ! Y length
DELX = 0.2 ! X spacing
DELY = 0.05 ! Y spacing
grav = 9.81 ! Gravity
et,1,141
!!!! meshing
rect,0,LX,0,LY
lsel,s,,,2,4,2
lesize,all,DELY
lsel,s,,,1,3,2
lesize,all,DELX
amesh,all
nsel,s,loc,y,0.0 ! Wall boundary conditions
nsel,a,loc,x,0.0
d,all,vx,0.0
d,all,vy,0.0
lsel,s,line,,2
nsll,s,1
d,all,vx
d,all,vy,
nsel,s,loc,x,0.0,5.0 ! Setup initial VFRC
nsel,r,loc,y,0.0,1.0
esln,s,1
ice,all,vfrc,1.0
nsel,s,loc,x,3,3.2 ! 冲击
nsel,r,loc,y,0.95,1.0
!d,all,vy,-90
d,all,pres,5e5
alls
/solu
!!!! FLOTRAN input
acel,0.0,grav,0
FLDA,NOMI,DENS,1.0e3
FLDA,NOMI,VISC,1.0d-3
/SOL
FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,0
FLDATA1,SOLU,TURB,0
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,1
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,0
FLDATA1,SOLU,RDSF,1
!*
/COM,,Transient Analysis,1
*SET,_z4,-1
*SET,_z14,1
FLDATA4,TIME,STEP,0.007,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,1.0e06
FLDATA4,TIME,TEND,0.14,
FLDATA4,TIME,GLOB,20,
FLDATA4,TIME,VX,0.01,
FLDATA4,TIME,VY,0.01,
FLDATA4,TIME,VZ,0.01,
FLDATA4,TIME,PRES,1e-006,
FLDATA4,TIME,TEMP,1e-006,
FLDATA4,TIME,ENKE,0.01,
FLDATA4,TIME,ENDS,0.01,
FLDATA4A,STEP,OVER,0
FLDATA4,TIME,OVER,0,
FLDATA4A,STEP,APPE,0
FLDATA4,TIME,APPE,0.01,
FLDATA4A,STEP,SUMF,0
FLDATA4,TIME,SUMF,1000000,
FLDATA4,TIME,BC,0
!*
FLDATA34,MIR,MOME,0.5,
FLDATA34,MIR,TEMP,0.5,
FLDATA34,MIR,TURB,0.5,
alls
solve
效果如下:
|