一个声学仿真“问题”
小弟最近学习ＡＮＳＹＳ，想用ＡＮＳＹＳ做一个声学方面的仿真，可是求解完成后发现声压十分的有问题，附上命令流，斑竹和大大留步解答一下。里面的一些warning能不能解决？建模和加载有什么不妥？大家执行的时候可以把载荷步减小一点，网格化粗一些，因为很费时。各位大大和斑竹费心了。/PREP7
/TITLE,acoustic TL cal
ET,1,solid45 ! structural element
et,2,fluid30 !coupled interface
et,3,fluid130 !infinite element
r,3,1.9,0,0,0 !real set for infinite element
et,4,fluid30,,1,0 !general 3D acoustic element
SAVE
MP,EX,1,7e7 !material property
MP,DENS,1,2700
MP,NUXY,1,0.33
MP,SONC,1,3080
MP,DENS,2,1.21
MP,SONC,2,344
MP,MU,1,0
SAVE
K,1,0.5,0.5,0 !modal
K,2,0.5,-0.5,0
K,3,-0.5,-0.5,0
K,4,-0.5,0.5,0
LSTR, 4, 1
LSTR, 1, 2
LSTR, 2, 3
LSTR, 3, 4
SAVE
FLST,2,4,4
FITEM,2,1
FITEM,2,4
FITEM,2,3
FITEM,2,2
AL,P51X
FLST,2,1,5,ORDE,1
FITEM,2,1
VEXT,P51X, , ,0,0,0.0008,,,,
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , , ,-0.0004, , ,1
SAVE
CYL4,0,0,2, , , ,0.0008
FLST,3,1,6,ORDE,1
FITEM,3,2
VGEN, ,P51X, , , , ,-0.0004, , ,1
SAVE
VOVLAP,1,2
SAVE
sphere,0,1.5,0,360
SAVE
SPHERE,1.9,1.5,0,360,
SAVE
FLST,2,4,6,ORDE,2
FITEM,2,1
FITEM,2,-4
VOVLAP,P51X
SAVE
VSEL,S, , , 1 !meshing volume 1
VPLOT
TYPE, 1
MAT, 1
ESYS, 0
SECNUM,
SAVE
FLST,5,1,4,ORDE,1
FITEM,5,11
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
MSHAPE,0,3D
MSHKEY,1
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
SAVE
NSLA,S,1 !applying boundary to volume 1
NSEL,R,LOC,X,0.5
FLST,2,96,1,ORDE,10
FITEM,2,2
FITEM,2,17
FITEM,2,-31
FITEM,2,258
FITEM,2,273
FITEM,2,-287
FITEM,2,513
FITEM,2,-516
FITEM,2,577
FITEM,2,-636
D,P51X, , , , , ,ALL, , , , ,
SAVE
NSLA,S,1
NSEL,R,LOC,X,-0.5
FLST,2,96,1,ORDE,14
FITEM,2,1
FITEM,2,32
FITEM,2,47
FITEM,2,-60
FITEM,2,257
FITEM,2,288
FITEM,2,303
FITEM,2,-316
FITEM,2,517
FITEM,2,-520
FITEM,2,637
FITEM,2,-640
FITEM,2,697
FITEM,2,-752
D,P51X, , , , , ,ALL, , , , ,
SAVE
NSLA,S,1
NSEL,R,LOC,Y,0.5
FLST,2,96,1,ORDE,6
FITEM,2,1
FITEM,2,-16
FITEM,2,257
FITEM,2,-272
FITEM,2,513
FITEM,2,-576
D,P51X, , , , , ,ALL, , , , ,
SAVE
NSLA,S,1
NSEL,R,LOC,Y,-0.5
FLST,2,96,1,ORDE,10
FITEM,2,17
FITEM,2,32
FITEM,2,-46
FITEM,2,273
FITEM,2,288
FITEM,2,-302
FITEM,2,577
FITEM,2,-580
FITEM,2,637
FITEM,2,-696
D,P51X, , , , , ,ALL, , , , ,
SAVE
FLST,5,2,6,ORDE,2 !select zhangban
FITEM,5,9
FITEM,5,11
VSEL,S, , ,P51X
VPLOT
SAVE
TYPE, 1 !mesh zhangban
MAT, 1
ESYS, 0
SECNUM,
ESIZE,0.05,0,
MSHAPE,1,3D
MSHKEY,0
FLST,5,2,6,ORDE,2
FITEM,5,9
FITEM,5,11
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
SAVE
FLST,5,2,6,ORDE,2 !boundaries of zhang ban
FITEM,5,9
FITEM,5,11
VSEL,S, , ,P51X
VPLOT
FLST,2,8915,1,ORDE,8
FITEM,2,1
FITEM,2,-60
FITEM,2,257
FITEM,2,-316
FITEM,2,513
FITEM,2,-752
FITEM,2,1537
FITEM,2,-10091
D,P51X, , , , , ,ALL, , , , ,
SAVE
FLST,5,4,6,ORDE,2 !spacial acoustic et
FITEM,5,5
FITEM,5,-8
VSEL,S, , ,P51X
VPLOT
SAVE
TYPE, 2
MAT, 2
ESYS, 0
SECNUM,
ESIZE,0.05,0,
MSHAPE,1,3D
MSHKEY,0
FLST,5,4,6,ORDE,2
FITEM,5,5
FITEM,5,-8
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
SAVE
esel,s,type,,1 !modify acoustic et
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodify,all,type,4
alls
SAVE
ASEL,S, , , 23 !infinite boundary
nsla,s,1
type,3
real,3
mat,2
esurf
SAVE
ASEL,S, , , 24
nsla,s,1
type,3
real,3
mat,2
esurf
SAVE
ASEL,S, , , 26
nsla,s,1
type,3
real,3
mat,2
esurf
SAVE
ASEL,S, , , 27
nsla,s,1
type,3
real,3
mat,2
esurf
SAVE
VSEL,S, , , 6 !apply loads on surface
VPLOT
FLST,2,2,5,ORDE,2
FITEM,2,20
FITEM,2,-21
SFA,P51X,1,PRES,
SAVE
asel,s,,,1 !mark coulping surface
asel,a,,,2
asel,a,,,19
asel,a,,,22
asel,a,,,25
asel,a,,,28
NSLA,S,1
esel,s,type,,2
sf,all,fsi,1
alls
SAVE
ANTYPE,3 !loader step
HROPT,FULL
HROUT,ON
LUMPM,0
EQSLV,FRONT,1e-008,
PSTRES,0
FLST,2,1208,1,ORDE,2
FITEM,2,1537
FITEM,2,-2744
SAVE
HARFRQ,0,500,
NSUBST,5,
KBC,1
SAVE
GPLOT !choose all and save all
esel,all
nsel,all
SAVE
finish
/SOLU !solve
/STATUS,SOLU
SOLVE
页:
[1]