/batch,list
/verify,dv-1101c
/title, dv-1101c, Flotran Remeshing: Cylinder pass through a channel
/COM,****************************************************
/COM,* Cylinder pass through a channel (pure fluid flow)
/COM,* Verify partial domain re-meshing
/COM,****************************************************
r1 = 2.5
l1 = 5.0
l2= 20.0
l3=60.0
h = 10.0
rx0= l1+(l2-l1)/2.0
ry0 = h/2.0
nl1 = 5
nl3 = 10
nr=8
nc=8
/PREP7
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*SET,_FNCNAME,'DIS' !DIS以什么类型附给_FNCNAME 因为上面没有定义?
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),1.0
! /INPUT,force.func
*DIM,%_FNCNAME%,TABLE,6,3,1 !_FNCNAME 是7行4列的吗? 下面赋值确实是这样的!
!
! Begin of equation: a*{TIME}
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*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, -1, 0, 1, 17, 3, 1
*SET,%_FNCNAME%(0,2,1), 0.0, 99, 0, 1, -1, 0, 0
*SET,%_FNCNAME%(0,3,1), 0
! End of equation: a*{TIME}
!-->
RECTNG,0,l1,0,h,
RECTNG,l1,l2,0,h,
RECTNG,l2,l3,0,h,
aglue,all
CYl4,rx0,ry0,0,45,r1,135
CYl4,rx0,ry0,0,135,r1,225
CYl4,rx0,ry0,0,225,r1,315
CYl4,rx0,ry0,0,315,r1,405
asel,s,,,2
asel,a,,,3
asel,a,,,6
asel,a,,,7
allsel,below,area
aglue,all
allsel
FLST,3,4,5,ORDE,3
FITEM,3,2
FITEM,3,8
FITEM,3,-10
ASBA, 4,P51X
aplot
et,1,141
KEYOPT,1,4,1
type,1
mat,1
asel,s,,,1
esize,.5
mshape,1,2d
mshkey,0
amesh,all
et,2,141
KEYOPT,2,4,1
type,2
mat,1
asel,s,,,3
esize,.5
mshape,1,2d
mshkey,0
amesh,all
et,3,141
KEYOPT,3,4,1
type,3
mat,1
asel,s,,,5
esize,2.0
mshape,1,2d
mshkey,0
amesh,all
nsel,s,loc,x,0
D,ALL,VX,0.0,
D,ALL,VY, 0.0
D,ALL,UX,0.0,
D,ALL,UY, 0.0
nsel,s,loc,x,l3
D,ALL,pres,0.0
D,ALL,UX,0.0,
D,ALL,UY, 0.0
nsel,s,loc,y,0
D,ALL,vy,0.0
D,ALL,UX,0.0,
D,ALL,UY, 0.0
nsel,s,loc,y,h
D,ALL,vy,0.0
D,ALL,UX,0.0,
D,ALL,UY,0.0,
lsel,s,,,5
lsel,a,,,23
lsel,a,,,25
lsel,a,,,27
nsll,s,1
d,all,UX,%DIS% !这个%DIS% 是指 _FNCNAME里面的内容吗?
d,all,UY,0.0
d,all,VX,1.0
d,all,VY,0.0
allsel
cdwrite,db,fluid,cdb,
fini
! Flotran Setup
/solu
FLDATA30,QUAD,MOMD,2,
FLDATA30,QUAD,MOMS,2,
FLDATA30,QUAD,PRSD,2,
FLDATA30,QUAD,PRSS,2,
FLDATA30,QUAD,THRD,0,
FLDATA30,QUAD,THRS,0,
FLDATA30,QUAD,TRBD,0,
FLDATA30,QUAD,TRBS,2,
/solu
FLDATA4,TIME,NUMB,100000,
FLDATA4,TIME,TEND,44.0,
FLDA,SOLU,ALE,T ! ALE solution
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TRAN,1
!FLDATA1,SOLU,TURB,1
FLDATA2,TIME,GLOB,5
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,1.0
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,1.0
!FLDA,BULK,BETA,1.0e5
FLDA,TIME,STEP, 2.0
!
!!! Newmark method
!
FLDATA,OUTP,TAUW,T
FLDATA,TIME,METH,NEWM
FLDATA,TIME,DELT,0.5
!
!! Set ANSYS-STRUCTURE commands
SAVE
/COM
/COM Re-meshing Commands
/COM
FLDATA,REMESH,ELEM,ALL ! all defined element re-meshing
FLDATA,REMESH,ARMA,10.0 ! maximum aspect ratio
FLDATA,REMESH,VOCH,5.0 ! maximum volume change
FLDATA,REMESH,ARCH,5.0 ! maximum aspect ratio change
SOLVE
SAVE
finish
这是一个ANSYS里自带的例子,现在有些地方没看懂我用红色标出来了请高手帮忙指教一下吧!谢谢 |