dyna做声学分析

superpoliceman

现在做一消声器的瞬态声场，老板要求用DYNA做，
有木有人做过？

superpoliceman

命令流如下
/units,si
/prep7
et,1,200
keyopt,1,1,6
keyopt,1,2,0
et,2,solid164
et,3,shell163
KEYOPT,3,1,9
R,1
RMODIF,1,1,5/6, ,0.01,0.01,0.01,0.01,   
!设定空气和刚体、膜材料
mp,dens,1,1.297
mp,ex,1,0
mp,nuxy,1,0
tb,eos,3,,,2,1
mp,dens,2,10e10
mp,ex,2,2e13
mp,nuxy,2,0
mp,dens,3,2500
mp,ex,3,3e10
mp,nuxy,3,0.01

r1=0.3
r2=0.84
m0=0.6
m1=1.2
!定义关键点和材料属性
k,1,
k,2,0,0,m0
k,3,,,m0+m1
k,4,,,m0+m1+m0
!
*do,i,1,3
lstr,i,i+1
*enddo
!
kwpave,2
pcirc,r1,,90,270,
local,11,1,0,0,m0,0,,0,1,1,
csys,11
k,101,r2,0
k,102,r2,90
k,103,r2,180
k,104,r2,270
a,101,102,103,104
kwpave,11
wpro,0,0,90    !绕Y轴逆时针旋转90
asbw,2
adele,3,,,1
asba,4,1,,,keep
vdrag,1,,,,,,1
vdrag,1,,,,,,2,3
vdrag,2,,,,,,2
vglue,all
csys,0
esize,0.1
type,1
mshkey,0       !自由分网
mshape,0,2d    !2维4边形
amesh,1
amesh,2
type,2
*do,i,1,3
vsweep,i
*enddo
vsweep,5
allsel,all
esel,s,type,,1
eplot
aclear,1
aclear,2
allsel
eplot


asel,s,,,7
aatt,3,1,3
amesh,all
allsel

csys,1
NSEL,S,LOC,X,R1
NSEL,U,LOC,Z,m0+0.001,m1+m0-0.001
esln
CM,ND1,elem
ALLSEL
NSEL,S,LOC,Z,m0
NSEL,R,LOC,X,r1-0.001,r2+0.001
esln
CM,ND2,elem
ALLSEL
NSEL,S,LOC,Z,m1+m0
NSEL,R,LOC,X,R1-0.001,R2+0.001
esln
CM,ND3,elem
ALLSEL
NSEL,S,LOC,X,R2
esln
CM,ND4,elem
CMGRP,NODALL,ND1,ND2,ND3,ND4
CMSEL,S,NODALL
mpchg,2,all
allsel
ldele,1,3,1,1
allsel
csys,0
NUMMRG,ALL

edpart,creat
partsel,'plot'

nsel,s,loc,x,-0.01,0.01
d,all,ux,0
allsel
finish
/solu
nsel,s,loc,z,0
cm,load,nodes,

*dim,time,array,3,1,1,
time(1,1,1)=0
time(2,1,1)=1.5
time(3,1,1)=3
*dim,compress,array,3,1,1,
compress(1,1,1)=0
compress(2,1,1)=1
compress(3,1,1)=0

edload,add,fz,0,load,time,compress,0,,
allsel,all
edgcale,1,1,1
edale,add,all,-1,0,0,0,0,0,1e20,1
edenergy,1,1,1,1
time,3
EDOPT,ADD,blank,BOTH
EDRST,100,  
EDHTIME,100,
EDDUMP,1,   
EDINT,5,0,  
SAVE
solve