求助：ansys仿真压电陶瓷振动在混凝土内的声场或应力分布问题

哪位大侠会用ANSYS做压电陶瓷振动时在混凝土内部辐射的声场分布啊？应力的分布也可以。我的压电陶瓷用的SOLID5，混凝土用的SOLID45，并把压电陶瓷和混凝土面接触的地方做了题粘接，现在已经建模和划分完网格了，但是不知道需要用哪种分析方法，怎么加载载荷和约束，希望大侠指点啊！代码如下：
fini
/clear
/prep7
k,1
k,2,,,50
k,3,,,50+2
l,1,2
l,2,3
Blc5,0,0,100,100
circle,1,6.25
csys,1
lsel,s,loc,z,0
lsel,r,loc,x,0,20
lgen,2,all,,,,45,,,,1
allsel
lplot
csys,0
L,7,9
L,6,8
L,4,10
L,5,11
allsel
lsel,u,line,,1
lsel,u,line,,2
asbl,all,all
CM,A1,area
Vdrag,A1,,,,,,1
csys,1
asel,s,loc,z,50
asel,r,loc,x,0,20
Vdrag,All,,,,,,2
csys,0
allsel
vglue,all
Vsel,s,loc,z,0,50
CM,V2,volume
Vsel,s,loc,z,50,52
CM,V1,volume
allsel
!赋属性
ET,1,solid5 !压电陶瓷
ET,2,solid45 !混凝土
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7500
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,1,,370
MPDATA,PERY,1,,370
MPDATA,PERZ,1,,635
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39e+011,7.78e+010,7.43e+010,0,0,0
TBDATA,,1.39e+011,7.43e+010,0,0,0,1.15e+011
TBDATA,,0,0,0,3.06e+010,0,0
TBDATA,,2.56e+010,0,2.56e+010,,,
TB,PIEZ,1,,,0
TBMODIF,1,1,0
TBMODIF,1,2,0
TBMODIF,1,3,-5.2
TBMODIF,2,1,0
TBMODIF,2,2,0
TBMODIF,2,3,-5.2
TBMODIF,3,1,0
TBMODIF,3,2,0
TBMODIF,3,3,15.1
TBMODIF,4,1,0
TBMODIF,4,2,0
TBMODIF,4,3,0
TBMODIF,5,1,0
TBMODIF,5,2,12.7
TBMODIF,5,3,0
TBMODIF,6,1,12.7
TBMODIF,6,2,0
TBMODIF,6,3,0
TBDE,PIEZ,1
TB,PIEZ,1,,,0
TBMODIF,1,1,0
TBMODIF,1,2,0
TBMODIF,1,3,-5.2
TBMODIF,2,1,0
TBMODIF,2,2,0
TBMODIF,2,3,-5.2
TBMODIF,3,1,0
TBMODIF,3,2,0
TBMODIF,3,3,15.1
TBMODIF,4,1,0
TBMODIF,4,2,0
TBMODIF,4,3,0
TBMODIF,5,1,0
TBMODIF,5,2,12.7
TBMODIF,5,3,0
TBMODIF,6,1,12.7
TBMODIF,6,2,0
TBMODIF,6,3,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,24000
MPDATA,PRXY,2,,0.2
TBDE,ANEL,1,,,
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39E+011,7.78E+010,7.43E+010,0,0,0
TBDATA,,1.39E+011,7.43E+010,0,0,0,1.15E+011
TBDATA,,0,0,0,3.06E+010,0,0
TBDATA,,2.56E+010,0,2.56E+010,,,
TB,CONC,2,1,9,
TBTEMP,0
TBDATA,,0.35,1,3.1125,-1,0,0
TBDATA,,0,0,0,,,
Vsel,s,Volume,,V1
VATT,1,,1,0
! the rubber matrix
Vsel,s,volume,,V2
VATT,2,,2,0
allsel
!网格划分
MSHAPE,0,3D
MSHKEY,1
csys,0
Lsel,all
Lsel,u,Line,,1
Lsel,u,Line,,2
Lsel,r,loc,z,1,49
Lesize,all,,,10
Lsel,all
Lsel,u,Line,,1
Lsel,u,Line,,2
Lsel,r,loc,z,50.1,51.9
Lesize,all,,,2
Lsel,s,Loc,z,0
Lsel,a,Loc,z,50
Lsel,a,Loc,z,52
Lesize,all,,,10
!two.划分网格
allsel
Vsweep,all