求助：分析圆管内声压分布 （附APDL）请高手指点...

siyola

用ANSYS做的圆管内声压分布的分析，具体描述为: 一个内径200，高1.5m的封闭圆管内充满水，底部由直径200的压电叠堆激励，求不同频率下的声压分布 频率范围为1kHz至20kHz。我现在的分析中没用将厚壁圆管建模，不知道这个思路对不对啊... 程序我直接粘贴在下边了 请高手帮忙指出问题...
/PREP7  
/TITLE,analysis of resonance   
ANTYPE,HARMIC
ET,1,FLUID29   !接触壁面和振动表面的声学单元
ET,2,FLUID29,,1  !内部声场单元（不接触壁面）
XDIV=10
YDIV=150! NUMBER OF DIVISIONS ALONG X,Y-AXIS
LEN=0.1
HGT=1.5
K,1
K,2,LEN
K,3,LEN,HGT
K,4,0,HGT
L,1,2,XDIV
L,2,3,YDIV
L,3,4,XDIV
L,4,1,YDIV
A,1,2,3,4
AMESH,ALL
MP,DENS,1,998
MP,SONC,1,1483
MP,MU,1,0.04 !管壁的边界导纳
MP,DENS,2,998!内部声场单元
MP,SONC,2,1483
MP,MU,2,0 ! 内部单元吸收系数为0
ESEL,ALL
NSEL,S,LOC,X,0,0.1
NSEL,A,LOC,Y,1.5,1.5
NSEL,A,LOC,X,0,0
NSEL,A,LOC,Y,0,1.5
NSEL,A,LOC,X,0.1
NSEL,A,LOC,Y,0,1.5
D,ALL,UX,,,,,UY !约束三个壁面位移为0
NSEL,S,LOC,X,0,0.1
NSEL,A,LOC,Y,0,1.5
ENODE
EINV
TYPE,2
MAT,2
EMODIF,ALL
ESEL,ALL !定义内部声场单元为2和材料2
!FSI 流固耦合的interface
NSEL,S,LOC,X,0,0.1
NSEL,A,LOC,Y,0,0
ENODE
SF,ALL,FSI
NROTAT,ALL
D,ALL,UY,1 !定义振动幅值
D,ALL,UX
NSEL,ALL
ESEL,ALL
FINISH
/SOLVE
HARFRQ,1000,1000
SOLVE
FINISH