fini
/clear
/filname,jjj,1
/prep7
factor = 1000
elasfac = 1
vmax = 5
diam = 20.4/factor
radius = diam/2.
thick = 1.8/factor
length = 100./factor
mp,ex,2,3.0E5*elasfac
mp,dens,2,1150
mp,nuxy,2,0.3
et,1,141,,,1,1
et,2,42,,,1
!! Blood volume
rect,0,radius,0,length
lesize,1,,,10,-3
lesize,3,,,10,-3
lesize,4,,,60,2
lesize,2,,,60,1/2
type,1
amesh,1
!! Blood vessel wall
type,2
mat,2
real,2
rect,radius,radius+thick,0,length
lesize,6,,,50
lesize,8,,,50
lesize,5,,,2
lesize,7,,,2
amesh,2
alls
!!! The pressure pulse table
*DIM,pulse,TABLE,5,1,1,time
!*
*SET,PULSE(1,0,1) , 0
*SET,PULSE(2,0,1) , 0.003
*SET,PULSE(2,1,1) , 1666
*SET,PULSE(3,0,1) , .008
*SET,PULSE(3,1,1) , 1666
*SET,PULSE(4,0,1) , .010
*SET,PULSE(4,1,1) , 0.
*SET,PULSE(5,0,1) , 1.00
*SET,PULSE(5,1,1) , 0.
/number,1
/pnum,mat,1
fini
/solu
asel,s,loc,x,0,radius
allsel,below,area
nsla,s,1
cm,fluid,node
asel,s,loc,x,radius,radius+thick
allsel,below,area
nsla,s,1
cm,solid,node
nsel,s,loc,y,0
nsel,a,loc,y,length
cmsel,r,solid
d,all,uy,0
alls
cmsel,s,solid
nsel,r,loc,x,radius
sf,all,fsin,1
alls
cmsel,s,fluid
nsel,r,loc,x,0
d,all,ux,0
d,all,vx,0
alls
cmsel,s,fluid
nsel,r,loc,y,0
d,all,uy,0
d,all,pres,%pulse%
alls
cmsel,s,fluid
nsel,r,loc,y,length
d,all,uy,0
d,all,pres,0
alls
cmsel,s,fluid
nsel,r,loc,x,radius
d,all,vx,0
d,all,vy,0
sf,all,fsin,1
alls
FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,0
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,1
FLDATA1,SOLU,RDSF,1
FLDATA12,PROP,DENS,0
FLDATA13,VARY,DENS,0
FLDATA12,PROP,VISC,0
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,1050,
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,4e-3,
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS
FLDATA7,PROT,COND,CONSTANT
FLDATA8,NOMI,COND,-1,
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT
FLDATA8,NOMI,SPHT,-1,
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
FLDATA15,PRES,REFE,101350,
FLDATA16,BULK,BETA,36,
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,293,
FLDATA14,TEMP,TTOT,293,
FLDATA14,TEMP,BULK,293,
TOFFSET,0,
*SET,_z4,-999
*SET,_z14,1
FLDATA4,TIME,STEP,0.001,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,10,
FLDATA4,TIME,TEND,1.0e06
FLDATA4,TIME,GLOB,15,
FLDATA4,TIME,VX,0.01,
FLDATA4,TIME,VY,0.01,
FLDATA4,TIME,VZ,0.01,
FLDATA4,TIME,PRES,1e-004,
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,1.0e6
FLDATA4A,STEP,SUMF,10,
FLDATA4,TIME,SUMF,1.0e6
FLDATA4,TIME,BC,0
FLDATA4,TIME,TEND,1000000,
FSAN,1
FSOR,FLUID
FSTR,FLUID,TRAN
FSTR,SOLID,TRAN
FSIN,NONC
FSTI,0.025,0.001,
FSDT,0.001,
FSIT,20,
FSCO,ALL,0.01,
FSOU,0.1E-01
FSCO,ALL,0.01,
FSOU,1
FSRE,ALL,0.8,
solve |