急！！！！！！！！请大家看看我的换能器水中与空气中分析（过程供大家借鉴）

yi1xiang

下面是我的命令流，建模中没有出现问题，求解设置都正常。
发送电压响应是错误的，我想请大家看看那里设置不合理。

yi1xiang

fini
/clea
/PREP7

con_1=0.02
con_2=0.01
con_3=0.00
con_4=0.00

r_0=0.800
r0=0.055
r1=0.005
r2=0.022
r3=0.005
r4=0.005
r5=0.037

h0=0.190
h1=0.008
h2=0.03
h3=0.003
h4=14*h3+0.02
h5=0.080
h6=0.053
h7=0.010
!r6=width
width_1=0.013

ET,1,13,7,,1
et,2,42,,,1

mp,DENS,1,7500
mp,perx,1,3.27e-9   
mp,perz,1,3.27e-9   
mp,pery,1,5.61e-9  
TBDE,ANEL,1,,,  
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39E+011,7.43E+010,7.78E+010,0,0,0
TBDATA,,1.15E+011,7.43E+010,0,0,0,1.39E+011
TBDATA,,0,0,0,2.56E+010,0,0
TBDATA,,2.56e10,0,3.06e10,,,
tb,piez,1   
tbdata,2,-5.2   
tbdata,5,15.1   
tbdata,8,-5.2   
tbdata,10,12.7  
tbdata,15,12.7
!硬铝  
UIMP,2,EX,,,7.15e10,
UIMP,2,DENS,,,2790,
UIMP,2,prxy,,,0.34,
!钢
UIMP,3,EX,,,21.6e10,
UIMP,3,DENS,,,7840,
UIMP,3,prxy,,,0.28,

k,108,r0,h2
k,6,r0,h1+h2
k,7,,h1+h2
k,109
k,10,r2,,

k,24,r2,-14*h3
kfill,10,24
k,25,r1,
k,39,r1,-14*h3
kfill,25,39

k,40,r3
k,41,r3,-h4+con_2
!k,42,r3+con_3,-h4+con_2
k,43,r4,-h4+con_2
k,44,r4,-h4
k,45,,-h4

!k,46,r4+con_4,-h4+h7
k,47,r4+con_4,-14*h3-h6
k,48,r5,-14*h3-h6
k,49,r5,-14*h3-h6+h5
k,50,r5-width_1,-14*h3-h6+h5
k,51,r5-width_1,-14*h3
k,52,r3+con_3,-14*h3

a,108,6,7,109,40,25,10
a,25,10,11,26
a,26,11,12,27
a,27,12,13,28
a,28,13,14,29
a,14,29,30,15
a,30,15,16,31
a,31,16,17,32
a,32,17,18,33
a,18,19,34,33
a,34,19,20,35
a,35,20,21,36
a,36,21,22,37
a,37,22,23,38
a,38,23,24,39
a,109,40,52,39,24,51,50,49,48,47,44,45
!a,43,44,47,48,49,50,51,24,39,52

esize,0.002
type,1
mat,1
amesh,2,15,1

type,2
mat,2
esize,0.0025
amesh,1
mat,3
esize,0.005
amesh,16
!定义电极
lsel,s,,,9,45,6
nsll,S,1
cp,1,volt,all
*get,n1,node,,num,min
d,n1,volt,1.0  
lsel,s,,,6,48,6
nsll,S,1
cp,2,volt,all   
d,all,volt,0.0
allsel,all
fini
/SOLU   
!*  
ANTYPE,2
!*

MODOPT,LANB,10   
EQSLV,SPAR  
MXPAND,10, , ,1  
!*  
MODOPT,LANB,10,10,50000, ,OFF

SOLVE   
FINISH  
/POST1  
SET,FIRST   
PLVECT,U, , , ,VECT,ELEM,ON,0

yi1xiang

fini
/clea
/PREP7

con_1=0.02
con_2=0.01
con_3=0.00
con_4=0.00

r_0=0.80
r0=0.055
r1=0.05
r2=0.022
r3=0.005
r4=0.005
r5=0.037

h0=0.190
h1=0.001
h2=0.04
h3=0.003
h4=14*h3+0.02
h5=0.075
h6=0.055
h7=0.010
!r6=width
width_1=0.0125
!定义单元类型，取轴对称
ET,1,13,7,,1
et,2,42,,,1
et,4,29,,,1             !包含结构自由度
et,5,29,,1,1            !不包含结构自由度
et,6,129,,,1

R,6,r_0,,,,
!压电陶瓷，极化方向y轴
mp,DENS,1,7500
mp,perx,1,3.27e-9   
mp,perz,1,3.27e-9   
mp,pery,1,5.61e-9  
TBDE,ANEL,1,,,  
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39E+011,7.43E+010,7.78E+010,0,0,0
TBDATA,,1.15E+011,7.43E+010,0,0,0,1.39E+011
TBDATA,,0,0,0,2.56E+010,0,0
TBDATA,,2.56e10,0,3.06e10,,,
!定义-[e]
tb,piez,1   
tbdata,2,5.2   
tbdata,5,-15.1   
tbdata,8,5.2   
tbdata,10,-12.7  
tbdata,15,-12.7
!硬铝  
UIMP,2,EX,,,7.15e10,
UIMP,2,DENS,,,2790,
UIMP,2,prxy,,,0.34,
!钢
UIMP,3,EX,,,21.6e10,
UIMP,3,DENS,,,7840,
UIMP,3,prxy,,,0.28,
!水
mp,dens,4,1030
mp,sonc,4,1500
mp,mu,4,0.0
mp,dens,5,1030
mp,sonc,5,1500
mp,mu,5,0.0
mp,dens,6,1030
mp,sonc,6,1500
mp,mu,6,1.0
!流体模型，以原点为中心，由1/4模型映射为1/2模型
CYL4,0.0,0, ,90,r_0,0  
BLC4,0.0,0.0,r0,h0/2
ASBA,1,2
aatt,5,,5
lesize,6,,,4,1
lesize,5,,,7,1
lesize,9,,,90,1
lesize,8,,,90,1
lesize,1,,,180,1
mshkey,0
mshape,0,2d
amesh,3
ARSYM,Y,3, , , ,0,0  
nummrg,all
!修改流体模型内边界
lrefine,4,7,1,3
lsel,s,,,4,7,1  
nsll,s,1  
esln,s,0
emodif,all,type,4
!修改流体模型外边界
csys,1  
nsel,s,loc,x,r_0
type,6  
real,6  
mat,6   
esurf
csys,0
   
!定义结构模型局部坐标系
local,11,0,0,h0/2-h1-h2
csys,11

k,108,r0,h2
k,109
k,10,r2,,

k,24,r2,-14*h3
kfill,10,24
k,25,r1,
k,39,r1,-14*h3
kfill,25,39

k,40,r3
k,41,r3,-h4+con_2
!k,42,r3+con_3,-h4+con_2
k,43,r4,-h4+con_2
k,44,r4,-h4
k,45,,-h4

!k,46,r4+con_4,-h4+h7
k,47,r4+con_4,-14*h3-h6
k,48,r5,-14*h3-h6
k,49,r5,-14*h3-h6+h5
k,50,r5-width_1,-14*h3-h6+h5
k,51,r5-width_1,-14*h3
k,52,r3+con_3,-14*h3

a,108,6,7,109,40,25,10
a,25,10,11,26
a,26,11,12,27
a,27,12,13,28
a,28,13,14,29
a,14,29,30,15
a,30,15,16,31
a,31,16,17,32
a,32,17,18,33
a,18,19,34,33
a,34,19,20,35
a,35,20,21,36
a,36,21,22,37
a,37,22,23,38
a,38,23,24,39
a,109,40,52,39,24,51,50,49,48,47,44,45
!a,43,44,47,48,49,50,51,24,39,52

esize,0.002
type,1
mat,1
amesh,4,17,1

type,2
mat,2
esize,0.0025
amesh,2

mat,3
esize,0.0025
amesh,18

!定义电极
lsel,s,,,17,53,6
nsll,S,1
cp,1,volt,all
*get,n1,node,,num,min
d,n1,volt,1.0             !正极
lsel,s,,,14,56,6
nsll,S,1
cp,2,volt,all
*get,n2,node,,num,min  
d,n2,volt,0.0             !负极
!于流体内边界除辐射面外，设置刚性边界条件
lsel,s, , ,5
nsll,s,1
d,all,ux,0
lsel,s, , ,7
nsll,s,1
d,all,ux,0
lsel,s, , ,4
nsll,s,1
d,all,uy,0
!于辐射面处流体模型，定义流固耦合
lsel,s,,,6
nsll,s,1
esln,s,0
esel,r,type,,4
sf,all,fsi,1
nsel,all
esel,all
!谐响应分析
/SOL
!*  
ANTYPE,3
!*  
!*  
HROPT,FULL  
HROUT,ON
LUMPM,0
!*  
EQSLV,FRONT,1e-008,
PSTRES,0
!*  
!*  
OUTPR,ALL,ALL,  
HARFRQ,5000,20000,  
NSUBST,100,  
KBC,1   
!*  
ALPHAD,0,   
BETAD,0,
DMPRAT,0.042, !定义阻尼4.2%  
!*  
ALLSEL,ALL  
/STATUS,SOLU

yi1xiang

水中的发送电压响应结果是锯齿状的，起伏很剧烈。声轴1.8米处的声压幅度、相位起伏也如此。远场条件已经满足了。

ymzhoull

我遇到过不少这种情况，有时候是129单元忘了给，有时候是svl的取值时没有取实部，不知道你是什么情况

shenda1304

楼主大虾，发射电压响应曲线应该怎样求啊，我知道那个求解公式但对于实现它还是一头雾水，你能贴出你做的代码让我们学习不，谢谢啦！

tzhdacy

需要做hormonic反应，然后对节点求声压解。记住流固耦合要用流固耦合单元，流体最外边边界用无反射边界条件。ok！