水灌内流固耦合问题Ansys实现
该题实现了一水罐在减速状态下，罐壁的应力分布和流体（水）的压强与速度分布/PREP7K,1
K,2,,1,,
K,3,,-1,,
K,4,,,1,
K,5,,,-1,
LARC,2,5,1,1,
LARC,2,4,1,1,
LARC,4,3,1,1,
LARC,3,5,1,1,
FLST,2,4,4
FITEM,2,2
FITEM,2,1
FITEM,2,4
FITEM,2,3
AL,P51X
FLST,2,1,5,ORDE,1
FITEM,2,1
ARSCALE,P51X, , ,1,0.5,0.95, ,0,1
FLST,2,1,5,ORDE,1
FITEM,2,1
VEXT,P51X, , ,5.1,0,0,,,,
SPHERE,1, ,0,360,
wpro,,,90.000000
VSBW, 2
wpro,,,-90.000000
FLST,3,1,6,ORDE,1
FITEM,3,4
VGEN, ,P51X, , , , ,5.1, , ,1
VPLOT
VDELE, 4, , ,1
FLST,2,1,6,ORDE,1
FITEM,2,3
VLSCAL,P51X, , ,0.3,0.5,0.95, ,0,1
VPLOT
FLST,3,1,6,ORDE,1
FITEM,3,3
VSYMM,X,P51X, , , ,0,0
FLST,3,1,6,ORDE,1
FITEM,3,3
VGEN, ,P51X, , , 5.1, , , , ,1
SAVE
WPOFF,1.55,,,
WPRO,,,90.000000
ASBW,3
WPRO,,,-90.000000
WPOFF,2,,,
wpro,,,90.000000
ASBW,14
CSYS,0
WPAVE,0,0,0
CSYS,0
wpro,,,-90.000000
wpoff,,,0.4
ASBW,15
WPOFF,,,0.3
ASBW,16
wpoff,,,0.1
ASBW,15
CSYS,0
WPAVE,0,0,0
CSYS,0
wpoff,1.55,,
wpro,,,90.000000
ASBW,6
wpoff,,,2
ASBW,19
CSYS,0
WPAVE,0,0,0
CSYS,0
wpro,,,-90.000000
wpoff,,,-0.5
ASBW,20
CSYS,0
WPAVE,0,0,0
CSYS,0
wpoff,,,0.4
ASBW,4
wpoff,,,-0.8
ASBW,5
CSYS,0
WPAVE,0,0,0
CSYS,0
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
FINISH
/PREP7
et,1,142,
et,2,93
BTOL,0.10e-9
vadd,1,2
vadd,3,4
BTOL,0.10e-4
VPLOT
asel,s,,,1
asel,a,,,4
asel,a,,,7,8
asel,a,,,12,21
asel,a,,,25,27
asel,a,,,30
aatt,2,2,2 ! (material 2)
alls
vsel,s,,,1
vatt,1,1,1 ! (material 1)
alls
!
!划分网格(存在疑问)
eshape,1
asel,s,,,1
asel,a,,,4
asel,a,,,7,8
asel,a,,,12,21
asel,a,,,25,27
asel,a,,,30
amesh,all
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!设置流体环境
et,1,142
et,2,0 ! becomes the Null Element(流固耦合场计算处理方法)
!! CFD Solution Control
flda,solu,flow,1 !Solves the momentum and pressure equations
flda,solu,turb,1 !Activates the turbulence model
flda,iter,exec,400 !Sets iteration and output controls for steady state analyses
flda,outp,sumf,10 !set to be defout
!! 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
FLDATA18,METH,PRES,1 ! TDMA
FLDATA19,TDMA,PRES,200,
FLDA,RELX,PRES,0.5 !
FLDA,RELX,Vx,0.5
FLDATA24,TURB,MODL,5
!! CFD Boundary Conditions (Applied to Solid Model)
asel,s,,,7
asel,a,,,12
asel,a,,,25
asel,a,,,27
da,all,vx,0.,1
asel,s,,,16,21
asel,a,,,8
asel,a,,,14
asel,a,,,30
asel,a,,,26
da,all,vy,0.,1
da,all,vz,0.,1
alls
/title,Fluid Analysis
physics,write,fluid,fluid
!
!设置固体环境
physics,clear
et,1,0 ! The Null element for the fluid region
et,2,93 ! Gasket element - material 2
R,2,0.001, , , , , ,
mp,ex,2,3.2947E+9 ! Young's modulus for rubber
mp,nuxy,2,0.4270 ! Poisson's ratio for the rubber
asel,s,,,1
asel,a,,,4
nsla,,1 !Select all nodes associated with the selected areas
d,all,ux,0.
d,all,uy,0.
d,all,uz,0. ! Fix the bottom of the gasket
alls
/title,structural analysis
finish
/solu
antype,static
nlgeom,on
cnvtol,f,,,,-1 !Sets convergence(一致性) values for nonlinear analyses;here,no minimum is enforced
physics,write,struc,struc
physics,clear
SAVE
/solu
physics,read,fluid
flda,iter,exec,50
vsel,s,,,1
nslv,s,0
BF,all,FORC,3.68,,, !Defines a nodal body force load
alls
solve
!vsel,s,,,1
!nslv,s,0
!BF,all,FORC,3.68,,, !Defines a nodal body force load(另外两个加速度值:5.21,6.25)
!alls
!solve
vsel,s,,,1
nslv,s,0
BFDELE,all,forc
fini
physics,read,struc ! Read in structures environment
/assign,esave,struc,esav ! Files for restarting nonlinear structure
/assign,emat,struc,emat
/solu
asel,s,,,1
asel,a,,,4
asel,a,,,7,8
asel,a,,,12,21
asel,a,,,25,27
asel,a,,,30
nsla,,1
esel,s,type,,2
ldread,pres,last,,,2,,rfl ! Apply pressure surface load from Flotran
alls
rescontrol,,none ! Do not use multiframe restart for nonlinear
solve
FINISH 设置非线性迭代收敛残差， 那个CONVERGENCE是收敛的意思，不是一致性。晕
搞错了，不好意思
页:
[1]