马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
~SATIN,12.5mm,SAT,,SOLIDS,0
/BATCH,LIST
/title, new analysize
/nopr
time =0.200
delta=0.001
apps=0.001
nts=0.2/delta
/prep7
*DEL,_FNCNAME
*DEL,_FNCMTID
*SET,_FNCNAME,'double'
! /INPUT,double.func
*DIM,%_FNCNAME%,TABLE,6,5,4
!
! Begin of equation: {TIME}
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 1, 0, 0
*SET,%_FNCNAME%(0,2,1), 0
*SET,%_FNCNAME%(0,3,1), 0
*SET,%_FNCNAME%(0,4,1), 0
*SET,%_FNCNAME%(0,5,1), 0
! End of equation: {TIME}
!
! Begin of equation: 0*{TIME}
*SET,%_FNCNAME%(0,0,2), 0.020, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 0, 0, 0, 1
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,3,2), 0, 99, 0, 1, -2, 0, 0
*SET,%_FNCNAME%(0,4,2), 0
*SET,%_FNCNAME%(0,5,2), 0
! End of equation: 0*{TIME}
!
! Begin of equation: 0.2*({TIME}-0.040)
*SET,%_FNCNAME%(0,0,3), 0.070, -999
*SET,%_FNCNAME%(2,0,3), 0.0
*SET,%_FNCNAME%(3,0,3), 0.0
*SET,%_FNCNAME%(4,0,3), 0.0
*SET,%_FNCNAME%(5,0,3), 0.0
*SET,%_FNCNAME%(6,0,3), 0.0
*SET,%_FNCNAME%(0,1,3), 1.0, -1, 0, 0.020, 0, 0, 1
*SET,%_FNCNAME%(0,2,3), 0.0, -2, 0, 1, 1, 2, -1
*SET,%_FNCNAME%(0,3,3), 0, -1, 0, 0.1, 0, 0, -2
*SET,%_FNCNAME%(0,4,3), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,5,3), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 0.2*({TIME}-0.040)
! Begin of equation: 0.2*(0.150-0.040) !!!!!!!!!算到稳态
*SET,%_FNCNAME%(0,0,4), 0.110, -999
*SET,%_FNCNAME%(2,0,4), 0.0
*SET,%_FNCNAME%(3,0,4), 0.0
*SET,%_FNCNAME%(4,0,4), 0.0
*SET,%_FNCNAME%(5,0,4), 0.0
*SET,%_FNCNAME%(6,0,4), 0.0
*SET,%_FNCNAME%(0,1,4), 1.0, -1, 0, 0.1, 0, 0, 0
*SET,%_FNCNAME%(0,2,4), 0.0, -2, 0, 0.050, 0, 0, -1
*SET,%_FNCNAME%(0,3,4), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,4), 0.0, 99, 0, 1, -3, 0, 0
*SET,%_FNCNAME%(0,5,4), 0
! End of equation: 50*0.040
!-->
!
et,1,142
VDELE, 1, , , 1
keyopt,1,4,1
vsel,s,,,2
type,1
mshape,1,3d
mshkey,0
esize,0.004
asel,s,,,49,64
aesize,all,0.002
vmesh,all
alls
!wall condition
asel,s,,,17
asel,a,,,19,40
asel,a,,,42,48
nsla,s,1
d,all,vx,0
d,all,vy,0
d,all,vz,0
d,all,ux,0
d,all,uy,0
d,all,uz,0
alls
!moving face
asel,s,,,49,64
nsla,s,1
d,all,ux,0.0
d,all,uy,0.0
!d,all,uz,0.0
d,all,uz,%double%
alls
! inlet condition
asel,s,,,41
nsla,s,1
d,all,ux,0.0
d,all,uy,0.0
d,all,uz,0.0
d,all,vx,1
!d,all,pres,1E+6
alls
!outlet
asel,s,,,18
nsla,s,1
d,all,ux,0.0
d,all,uy,0.0
d,all,uz,0.0
d,all,pres,1.0e+5
alls
! Flotran Setup
/solu
!kbc,0
!timint,on,fluid
NCNV,,1.0E25,0,0
flda,step,sumf,1
flda,time,numb,nts
flda,time,appe,apps
flda,time,tend,time
FLDA,SOLU,ALE,T ! ALE solution
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,TURB,1
FLDATA2,TIME,GLOB,5
FLDA,TIME,STEP,apps
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,1000
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,4.6
SOLVE
SAVE
FINISH |