压电纵振动换能器,在水中，谐响应分时，如何建立水环境并施加约束？

aiyouguai

压电纵振动换能器，前盖板硬铝，压电陶瓷PZT-4，后盖板和预应力杆45钢。模型是面对称，而非轴对称，建立1/4模型。
问题：
1，在水中，谐响应分时，如何建立水环境并施加约束？
2，对称约束-----但运行时，只有一个轴的对称约束，另一个轴不知道是不显示还是无法施加
3，水环境中，无法提取压力值
4，水环境建模很差，望指点
程序如下：
/PREP7
!
ET,1,SOLID45,0
ET,2,SOLID5,3
ET,3,FLUID30
ET,4,FLUID130
R,4,0.2,0,0,0.072
ET,5,FLUID30,,1
!2(PZT-4)
MP,DENS,2,7.5E3
!弹性系数
TB,ANEL,2
TBDATA,1,13.2E10,7.1E10,7.3E10
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA,21,2.6E10
!介电常数
MP,PERX,2,804.6
MP,PERY,2,804.6
MP,PERZ,2,659.7
!压电常数
TB,PIEZ,2
TBDATA,3,-4.1
TBDATA,6,-4.1
TBDATA,9,14.1
TBDATA,14,10.5
TBDATA,16,10.5
!(4#FE)
MP,DENS,1,7.8E3
MP,EX,1,2E11
MP,PRXY,1,0.28
!硬铝
MP,DENS,3,2.7E3
MP,EX,3,7.28E11
MP,PRXY,3,0.33
!水
MP,DENS,4,1030
MP,SONC,4,1460
!模型建立
K,1,0,0
K,2,0,5E-3
K,3,0,10E-3,20E-3
K,4,0,10E-3,23E-3
K,5,0,0,23E-3
K,6,33.5E-3,0
L,1,6
A,1,2,3,4,5
VDRAG,1,,,,,,1
!
CYLIND,0.002,,0,0.023,,90
VSBV,1,2
CYLIND,0.01,0.0025,0.023,0.059,,90
CYLIND,0.002,,,0.059,,90
CONE,0.01,0.0475,0.059,0.070,,90
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST
WPROTA,0,0,-90
WPOFFS,,,-0.0335
VSBW,4
VDELE,5,,,1
WPROTA,0,90,
WPOFFS,,0.0335,-0.0335
VSBW,6
VDELE,4,,,1
WPROTA,,,90
WPROTA,-90,,
WPOFFS,0,-0.0335,
WPOFFS,,,0.07
BLOCK,,0.0335,,0.0335,,0.002
WPOFFS,,,0.002
SPHERE,0.2,,,90
VSBW,6
VDELE,7,,,1
WPOFFS,,,-0.072
VGLUE,ALL
!材料赋值
VSEL,S,,,10
VSEL,A,,,7
VATT,1,,1
VSEL,S,,,1
VATT,2,,2
VSEL,S,,,9
VATT,3,,1
VSEL,S,,,4
VATT,4,,3
VSEL,S,,,6
VATT,4,,5
ALLSEL
!网格化
VSEL,S,,,1
VSWEEP,ALL
MSHAPE,1,3D
MSHKEY,0
FLST,5,5,6,ORDE,5   
FITEM,5,4   
FITEM,5,6   
FITEM,5,-7  
FITEM,5,9  
FITEM,5,-10
VSEL, , , ,P51X
VMESH,ALL
ALLSEL
!定义无限边界
CSYS,2
ASEL,S,,,40
NSLA,S,1  
TYPE,4   
MAT,4
REAL,4   
esurf
ALLSEL
!对称约束-----但运行时，只有一个轴的对称约束，另一个轴不知道是不显示还是无法施加
NSEL,S,LOC,X,0
DSYM,SYMM,X      
NSEL,S,LOC,Y,0
DSYM,SYMM,Y
NSEL,ALL
ASEL,S,,,1
NSLA,S,1
CP,1,VOLT,ALL
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,0
ASEL,S,,,2
NSLA,S,1
CP,2,VOLT,ALL
*GET,N2,NODE,,NUM,MIN
D,N2,VOLT,20
NSEL,ALL
/SOLU
EQSLV,SPARCE
DMPR,0.01
ANTYPE,HARMIC
OUTRES,ALL,ALL
HARFRQ,7000,20000
NSUBST,26
KBC,1
SOLVE
/POST26
RFORCE,2,N1,AMPS
PI2=3.14159*2
CFACT,0,-1
PROD,3,2,1,,MHOS,,,,PI2
REALVAR,4,3,,,REAL3
IMAGIN,5,3,,,IMAG3
PLVA,4,5

aiyouguai

原来是没有定义结构和耦合接触面

hy1966

看了你的模型图，似乎水环境是1/8模型，最好把换能器壳体一起淹没在水中，做1/4模型来计算

aiyouguai

谢谢版主关注
建1/4水环境，就会有很多接触面，这样流固接触面该如何定义呀？
望赐教

hy1966

方法和前端面一样，可以模拟换能器各个方向的辐射

aiyouguai

无自由度水环境，ET,5,FLUID30,,1中，没有压力是怎么回事呀？

feitianlong

6# aiyouguai
在流固耦合接触面上需要施加流固耦合边界条件FSI

bouquet

楼主有没有施加FSI？结果如何？
我正在做类似的仿真，但是读不出声压，FSI设置没有问题，分析的结果也没有warning提示，就是不知道发生什么问题

aiyouguai

施加了FSI，但提取的前盖板位移，和30声压均为锯齿状，不知错在什么地方。

aiyouguai

hy1966老师，您好！加FSI流体-结构耦合如下，可是分析结果结构位移却是锯齿状，请问问题何在，能不能帮我稍稍修改稿下程序，谢谢！
CSYS,0
ASEL,S,,,32
NSLA,S,1
ESEL,S,TYPE,,3
SF,ALL,FSI,1
ALLSEL

aiyouguai

晕，刚才贴错了，是这一段
CSYS,0
NSEL,S,LOC,Z,0.069
ESEL,S,TYPE,,3
SF,ALL,FSI,1
ALLSEL

hy1966

为什么不用个简单模型试试呢

tzhdacy

需要修改成你需要的流固耦合单元
emodif,all
SF,all,fsi
CM,nfsi,NODE