asel,u,,,2,3 amesh,all eshape,0 asel,s,,,2,3 amesh,all
!-----------------
!!!!!Create element plot and write to the file gasket.grph asel,s,,,1,3 esla,s
/Title, Initial mesh for gasket and neighborhood eplot
/ZOOM,1,RECT,0.3,-0.6,0.4,-0.5 alls
!----------------- !!
!!!!!2. Create Physics Environment for the Fluid et,1,141
et,2,0! Gasket becomes the Null Element
vin=3.5e-1! Inlet water velocity (meters/second) !!CFD Solution Control flda,solu,flow,1 flda,solu,turb,1 flda,iter,exec,400 flda,outp,sumf,10
!!CFD Property Information flda,prot,dens,constant flda,prot,visc,constant
flda,nomi,dens,1000.! 1000 kg/m3 for density - water flda,nomi,visc,4.6E-4! 4.6E-4 kg-s/m (viscosity of water) flda,conv,pres,1.E-8! Tighten pressure equation convergence !!CFD Boundary Conditions (Applied to Solid Model) lsel,s,,,8,17,9 lsel,a,,,20
dl,all,,vx,0.,1! Centerline symmetry lsel,s,,,9 dl,all,,vx,0.,1
dl,all,,vy,vin,1! Inlet Condition lsel,s,,,2 lsel,a,,,18,19 lsel,a,,,21,22
dl,all,,vx,0.,1! Outer Wall dl,all,,vy,0.,1 lsel,s,,,1,3,2 lsel,a,,,6
dl,all,,vx,0.,1! Gasket dl,all,,vy,0.,1 lsel,s,,,15
dl,15,,pres,0.,1! Outlet pressure condition
!!! create named component of nodes at the bottom of gasket lsel,s,,,1 nsll,,1 cm,gasket,node
nlist ! List initial nodal positions of the bottom of the gasket /com, +++++++++ STARTING gasket coordinates -------- alls
/title,Fluid Analysis physics,write,fluid,fluid !!
!!!!!3.Create Physics Environment for the Structure !!
physics,clear
et,1,0! The Null element for the fluid region et,2,56! Gasket element - material 2 mp,ex,2,2.82E+6! Young\mp,nuxy,2,0.49967! Poisson\tb,mooney,2
tbdata,1,0.293E+6!Mooney-Rivlin Constants tbdata,2,0.177E+6!\lsel,s,,,2 nsll,,1 d,all,ux,0.
d,all,uy,0.! Fix the end of the gasket
alls
/title,structural analysis finish /solu
antype,static nlgeom,on cnvtol,f,,,,-1
physics,write,struc,struc physics,clear save !!
!!!!!4.Fluid-Structure Interaction Loop !!
loop=25! Maximum allowed number of loops
toler=0.005! Convergence tolerance for maximum displacement *dim,dismax,array,loop! Define array of maximum displacement values *dim,strcri,array,loop! Define array of convergence values *dim,index,array,loop
*do,i,1,loop!Execute fluid -> structure solutions /solu
physics,read,fluid! Read in fluid environment *if,i,ne,1,then
flda,iter,exec,100! Execute 100 global iterations for *endif!each new geometry
solve! FLOTRAN solution fini
! end of fluid portion
physics,read,struc! Read in structures environment
/assign,esave,struc,esav ! Files for restarting nonlinear structure /assign,emat,struc,emat
*if,i,gt,1,then! Structural restart loop
parsave,all! Save parameters for convergence check resume! Resume DB - to return original node positions parresume! Resume parameters needed for convergence check /prep7
antype,stat,rest fini *endif /solu solc,off
lsel,s,,,1,3,2! Select proper lines to apply fluid pressures lsel,a,,,6! to the entire gasket surface nsll,,1 esel,s,type,,2
ldread,pres,last,,,,,rfl ! Apply pressure surface load from Flotran alls
rescontrol,,none! Do not use multiframe restart for nonlinear solve