/PREP7
ET,1,SOLID45 !壳单元 类型1
ET,2,FLUID30 ! acoustic fluid element with ux & uy(0-Structure present at interface (unsymmetric element matrix))
ET,3,130 ! acoustic infinite line element
r,3,0.3,0,0
ET,4,FLUID30,,1,0 ! acoustic fluid element without ux & uy(1-No structure at the interface (symmetric element matrix))
! material properties
!材料1为铁
MP,EX,1,2.1E11
MP,DENS,1,7800
MP,NUXY,1,0.3
!材料2为空气介质
MP,DENS,2,1.21
MP,SONC,2,344
!建立模型
wprot,0,-90,0
SPHERE,,0.3,0,180 !SPHERE, RAD1, RAD2, THETA1, THETA2
wprot,0,90,0
BLOCK,0,0.1, -0.002, 0.002, -0.1, 0 !BLOCK, X1, X2, Y1, Y2, Z1, Z2
BLOCK,-0.1,0, -0.002, 0.002, -0.1, 0
VGLUE,2,3
VOVLAP,1,2 !Overlaps volumes
VOVLAP,3,4
!对编号为2的体划分网格
lesize,9,,,20,1
lesize,19,,,20,1
lesize,14,,,20,1
lesize,18,,,20,1
lesize,11,,,20,1
lesize,16,,,20,1
lesize,12,,,20,1
lesize,17,,,20,1
lesize,8,,,4,1
lesize,15,,,4,1
lesize,13,,,4,1
lesize,10,,,4,1
VSEL,s,,,2
VATT,1,,1
mshkey,1 !Specifies whether free meshing or mapped meshing should be used to mesh a model.1— Use mapped meshing.
VMESH,2
alls
!对编号为4的体划分网格
lesize,32,,,20,1
lesize,31,,,20,1
lesize,35,,,20,1
lesize,33,,,20,1
lesize,28,,,20,1
lesize,34,,,20,1
lesize,20,,,4,1
lesize,27,,,4,1
VSEL,s,,,4
VATT,1,,1
mshkey,1 !Specifies whether free meshing or mapped meshing should be used to mesh a model.1— Use mapped meshing.
VMESH,4
alls
!对编号为1的体划分网格(此步骤划分网格需要手动进行)
VSEL,s,,,1
VATT,2,,2 !指定的材料属性为空气介质(与结构耦合)
SMRT,1
VMESH,1
alls
!修改不与平板(结构)接触单元的材料类型类型为空气介质(不与结构耦合)
esel,s,type,,1
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodif,all,type,4
alls
!在边界生成吸收单元。(使用单元fluid130)
Asel,s,,,1,2
NSLA,S,1
type,3
real,3
mat,2
esurf !Generates elements overlaid on the free faces of existing selected elements.
alls
!建立流-固接触面
Asel,s,,,8
Asel,A,,,10
Asel,A,,,7
Asel,A,,,5
Asel,A,,,19
Asel,A,,,15
Asel,A,,,17
Asel,A,,,20
NSLA,S,1
esel,s,type,,2
sf,all,fsi,1
alls
/solu
antype,trans
trnopt,full,nodamp
da,6,all
da,18,all
time,0.02
nsubst,10
kbc,1
outpr,basic,1
outres,all,1
f,1,fy,10
alls
solve
/post26
nsol,2,4632,pres,pres
plvar,2
得出某节点声压随时间变化曲线,但是某时间点的声压级却全是灰色。
|