第二章 顺序耦合场分析(9)

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

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库第二章 顺序耦合场分析(9)在线全文阅读。