满库
/BATCH
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
!*
/prep7
!* define material proterties
mp,dens,1,2650. !mat 1 for dam
mp,ex,1,3.15e10
mp,prxy,1,.167
mp,dens,2,1000. !mat 2 for water
mp,sonc,2,1440
!* define element type
et,1,PLANE42,,,2
et,2,29
et,3,29,,1
!* define geometry
lwater=618.
k, 1,0.0,0.0,0.0
k, 2,70.2,0.,0.
k, 3,0.,66.5,0.
k, 4,21.9875,66.5,0.
k, 5,0.,103.,0.
k, 6,14.8,103.,0.
k, 7,-1.*lwater,103.,0.
k, 8,-1.*lwater,66.5,0.
k, 9,-1.*lwater,0.,0.
a, 1, 2, 4, 3
a, 3, 4, 6, 5
a, 8, 3, 5, 7
a, 9,1, 3, 8
asel,s,loc,x,0.,100.
aatt,1,,1
cm,adam,area
asel,s,loc,x,-1*lwater,0.
aatt,2,,3
cm,awater,area
allsel,all
!* mesh geometry
ESIZE,0,10
lsel,s,loc,y,0.1,66. !坝下部剖分分数
lesize,all,,,15
lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数
lsel,r,loc,y,-0.1,67.
lesize,all,,,40,0.5
lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数
lsel,r,loc,y,100.,104.
lesize,all,,,40,2.0
mshape,0,2D
mshkey,1
allsel,all
amesh,all
!更改与水体接触的单元类型
esel,s,type,,3
nsle,s
nsel,r,loc,x,-1.,1.
esln,r
emodif,all,type,2
allsel,all
finish
/solu
antype,modal
MODOPT,UNSYM,30
MXPAND,30, , ,0
/pnum,type,1
/number,1
nsel,s,loc,x,-1.,1. !库水坝体FSI面
esln,s
esel,r,mat,,2
sf,all,fsi
esel,s,mat,,1 !坝体约束
nsle,s
nsel,r,loc,y,-1.0,1.0
d,all,ux,0.
d,all,uy,0.
esel,s,mat,,2 !水体表面约束
nsle,s
nsel,r,loc,y,102.,104.
d,all,pres,0.
/pbc,all,,1
/psf,fsi,,2
allsel,all
gplot
save
solve |