biomecanique
Bioengineer
- Jun 21, 2004
- 2
Hello everybody,
I'm working in domain biomecanique, I have a problem for applying a displacement in following example. I don't understand the number declared in Table. Thanks for any answer.
/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'
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),1.0
! /INPUT,force.func
*DIM,%_FNCNAME%,TABLE,6,3,1
!
! 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%
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
I'm working in domain biomecanique, I have a problem for applying a displacement in following example. I don't understand the number declared in Table. Thanks for any answer.
/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'
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),1.0
! /INPUT,force.func
*DIM,%_FNCNAME%,TABLE,6,3,1
!
! 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%
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