高手近来看看这个结构声耦合的算例

anc

想模拟一个地板上4个角都固定的房间。当屋顶上有力作用时， 房间内部的声场。发现有2个问题：
1。力需要手动给，如果直接用f命令给无法计算。（）
2。手动给了力，计算结果是房间内部的声压没有值，但是表面有位移。

不知道什么原因，望高手指点。
前面建模和其他边界条件的命令如下：

/prep7
et,1,30             !Acoustic elements in contact with vibrating surface
et,2,SOLID45
et,3,30,,1          !Acoustic elements in interior (not in contact with walls)

MP,DENS,1,1.225     !Acoustic elements in contact With walls
MP,SONC,1,340
MP,MU,1,0.08

MP,DENS,2,1.225     !Acoustic elements in contact With floor
MP,SONC,2,340
MP,MU,2,0.04
   
MP,EX,3,2.1E10      ! Concrete part
MP,DENS,3,2.8E3
MP,PRXY,3,0.2
MP,MU,3,0.5E-2   

MP,EX,4,2.16E11     ! Steel part
MP,DENS,4,7.8E3
MP,PRXY,4,0.3   

MP,DENS,5,1.21      ! Acoustic elements in interior
MP,SONC,5,344
MP,MU,5,0.0

k,1,0.24,0.24,0.24
k,2,5.24,0.24,0.24
k,3,5.24,4.24,0.24
k,4,0.24,4.24,0.24
k,5,0.24,0.24,3.24
k,6,5.24,0.24,3.24
k,7,5.24,4.24,3.24
k,8,0.24,4.24,3.24
v,1,2,3,4,5,6,7,8    ! air volume

k,9,0.24,0.24,0
k,10,5.24,0.24,0
k,11,5.24,4.24,0
k,12,0.24,4.24,0
k,13,0.24,0.24,3.48
k,14,5.24,0.24,3.48
k,15,5.24,4.24,3.48
k,16,0.24,4.24,3.48
v,1,2,3,4,9,10,11,12  ! concrete part z=0~0.24
v,5,6,7,8,13,14,15,16 ! concrete part z=3.24~3.28

k,17,0,0.24,0.24
k,18,0,4.24,0.24
k,19,0,4.24,3.24
k,20,0,0.24,3.24
k,21,5.48,0.24,0.24
k,22,5.48,4.24,0.24
k,23,5.48,4.24,3.24
k,24,5.48,0.24,3.24
v,1,5,8,4,17,20,19,18  ! concrete part x=0~0.24
v,2,6,7,3,21,24,23,22  ! concrete part x=5.24~5.48

k,25,0.24,0,0.24
k,26,5.24,0,0.24
k,27,5.24,0,3.24
k,28,0.24,0,3.24
k,29,0.24,4.48,0.24
k,30,5.24,4.48,0.24
k,31,5.24,4.48,3.24
k,32,0.24,4.48,3.24
v,1,5,6,2,25,28,27,26 ! concrete part y=0~0.24
v,4,8,7,3,29,32,31,30 ! concrete part y=4.24~4.48

k,33,0,0,0
k,34,0.24,0,0
k,35,5.24,0,0
k,36,5.48,0,0
k,37,5.48,0.24,0
k,38,5.48,4.24,0
k,39,5.48,4.48,0
k,40,5.24,4.48,0
k,41,0.24,4.48,0
k,42,0,4.48,0
k,43,0,4.24,0
k,44,0,0.24,0

k,45,0,0,3.48
k,46,0.24,0,3.48
k,47,5.24,0,3.48
k,48,5.48,0,3.48
k,49,5.48,0.24,3.48
k,50,5.48,4.24,3.48
k,51,5.48,4.48,3.48
k,52,5.24,4.48,3.48
k,53,0.24,4.48,3.48
k,54,0,4.48,3.48
k,55,0,4.24,3.48
k,56,0,0.24,3.48

k,57,0,0,0.24
k,58,5.48,0,0.24
k,59,5.48,4.48,0.24
k,60,0,4.48,0.24
k,61,0,0,3.24
k,62,5.48,0,3.24
k,63,5.48,4.48,3.24
k,64,0,4.48,3.24

v,1,25,26,2,9,34,35,10  !  z=0~0.24 steel elements except corners
v,2,21,22,3,10,37,38,11 !  from x=0 , inverse clock
v,3,11,12,4,30,40,41,29
v,4,12,9,1,18,43,44,17
v,5,13,14,6,28,46,47,27 !  z=3.24~3.48 steel elements except corners
v,6,14,15,7,24,49,50,23 !  from x=0 , inverse clock
v,7,15,16,8,31,52,53,32
v,8,16,13,5,19,55,56,20

v,5,28,25,1,20,61,57,17  !four steel elements vertical to xy plane except corners
v,6,24,21,2,27,62,58,26
v,7,31,30,3,23,63,59,22
v,8,19,18,4,32,64,60,29

v,1,25,34,9,17,57,33,44  !four corners between z=0~0.24
v,2,26,35,10,21,58,36,37
v,3,22,38,11,30,59,39,40
v,4,29,41,12,18,60,42,43

v,5,28,46,13,20,61,45,56  !four corners between z=3.24~3.48
v,6,24,49,14,27,62,48,47
v,7,31,52,15,23,63,51,50
v,8,19,55,16,32,64,54,53

lesize,93,,,20
lesize,56,,,20
lesize,26,,,20
lesize,10,,,20
lesize,22,,,20
lesize,6,,,20
lesize,83,,,20
lesize,48,,,20
lesize,60,,,20
lesize,73,,,20
lesize,3,,,20
lesize,18,,,20
lesize,1,,,20
lesize,14,,,20
lesize,63,,,20
lesize,52,,,20

lesize,98,,,16
lesize,32,,,16
lesize,28,,,16
lesize,12,,,16
lesize,36,,,16
lesize,78,,,16
lesize,4,,,16
lesize,20,,,16
lesize,24,,,16
lesize,8,,,16
lesize,88,,,16
lesize,40,,,16
lesize,2,,,16
lesize,16,,,16
lesize,44,,,16
lesize,68,,,16

lesize,11,,,12
lesize,118,,,12
lesize,34,,,12
lesize,54,,,12
lesize,5,,,12
lesize,38,,,12
lesize,108,,,12
lesize,50,,,12
lesize,7,,,12
lesize,30,,,12
lesize,46,,,12
lesize,103,,,12
lesize,113,,,12
lesize,58,,,12
lesize,42,,,12
lesize,9,,,12

lesize,92,,,3
lesize,111,,,3
lesize,139,,,3
lesize,140,,,3
lesize,91,,,3
lesize,141,,,3
lesize,112,,,3
lesize,57,,,3
lesize,25,,,3
lesize,41,,,3
lesize,90,,,3
lesize,89,,,3

lesize,86,,,3
lesize,84,,,3
lesize,138,,,3
lesize,137,,,3
lesize,87,,,3
lesize,21,,,3
lesize,136,,,3
lesize,85,,,3
lesize,49,,,3
lesize,37,,,3
lesize,106,,,3
lesize,107,,,3

lesize,144,,,3
lesize,143,,,3
lesize,96,,,3
lesize,94,,,3
lesize,95,,,3
lesize,142,,,3
lesize,27,,,3
lesize,97,,,3
lesize,117,,,3
lesize,53,,,3
lesize,116,,,3
lesize,33,,,3

lesize,135,,,3
lesize,99,,,3
lesize,134,,,3
lesize,81,,,3
lesize,100,,,3
lesize,23,,,3
lesize,82,,,3
lesize,133,,,3
lesize,102,,,3
lesize,101,,,3
lesize,29,,,3
lesize,45,,,3

lesize,119,,,3
lesize,120,,,3
lesize,55,,,3
lesize,35,,,3
lesize,77,,,3
lesize,130,,,3
lesize,75,,,3
lesize,19,,,3
lesize,131,,,3
lesize,74,,,3
lesize,76,,,3
lesize,132,,,3

lesize,31,,,3
lesize,47,,,3
lesize,105,,,3
lesize,104,,,3
lesize,80,,,3
lesize,61,,,3
lesize,15,,,3
lesize,121,,,3
lesize,23,,,3
lesize,79,,,3
lesize,62,,,3
lesize,122,,,3

lesize,39,,,3
lesize,66,,,3
lesize,109,,,3
lesize,51,,,3
lesize,110,,,3
lesize,13,,,3
lesize,124,,,3
lesize,64,,,3
lesize,67,,,3
lesize,65,,,3
lesize,126,,,3
lesize,125,,,3

lesize,114,,,3
lesize,115,,,3
lesize,59,,,3
lesize,43,,,3
lesize,127,,,3
lesize,69,,,3
lesize,72,,,3
lesize,17,,,3
lesize,71,,,3
lesize,70,,,3
lesize,129,,,3
lesize,128,,,3

VMESH,ALL

!ELEMENT DEFINITION FOR STEEL BEAMS

NSEL,S,LOC,X,0,0.16
NSEL,R,LOC,Y,0,0.16
ESLN
NSEL,A,LOC,X,5.32,5.48
NSEL,R,LOC,Y,0,0.16
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,5.32,5.48
NSEL,R,LOC,Y,4.32,4.48
ESLN
NSEL,A,LOC,X,0,0.16
NSEL,R,LOC,Y,4.32,4.48
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,0,0.16
NSEL,R,LOC,Z,0,0.16
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,0,0.16
NSEL,R,LOC,Z,3.32,3.48
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,4.32,4.48
NSEL,R,LOC,Z,0,0.16
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,4.32,4.48
NSEL,R,LOC,Z,3.32,3.48
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,5.32,5.48
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,0,0.16
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,5.32,5.48
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,3.32,3.48
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0,0.16
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,0,0.16
ESLN
TYPE,2
MAT,4
EMODIF,ALL

NSEL,S,LOC,X,0,0.16
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,3.32,3.48
ESLN
TYPE,2
MAT,4
EMODIF,ALL

!ELEMENT DEFINITION FOR CONCRETE PLATE

NSEL,S,LOC,X,0.49,4.99
NSEL,R,LOC,Y,0.49,3.99
NSEL,R,LOC,Z,0,0.16
ESLN
TYPE,2
MAT,3
EMODIF,ALL

NSEL,S,LOC,X,0.49,4.99
NSEL,R,LOC,Y,0.49,3.99
NSEL,R,LOC,Z,3.32,3.48
ESLN
TYPE,2
MAT,3
EMODIF,ALL

NSEL,S,LOC,X,0,0.16
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,0.32,3.16
ESLN
TYPE,2
MAT,3
EMODIF,ALL

NSEL,S,LOC,X,5.32,5.48
NSEL,R,LOC,Y,0.32,4.16
NSEL,R,LOC,Z,0.32,3.16
ESLN
TYPE,2
MAT,3
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,0,0.16
NSEL,R,LOC,Z,0.32,3.16
ESLN
TYPE,2
MAT,3
EMODIF,ALL

NSEL,S,LOC,X,0.32,5.16
NSEL,R,LOC,Y,4.32,4.48
NSEL,R,LOC,Z,0.32,3.16
ESLN
TYPE,2
MAT,3
EMODIF,ALL

!ELEMENT DEFINITION FOR iNTERIOR ACOUSTIC
NSEL,S,LOC,X,0.74,4.74
NSEL,R,LOC,Y,0.74,3.74
NSEL,R,LOC,Z,0.74,2.74
ESLN
TYPE,3
MAT,5
EMODIF,ALL

!ELEMENT DEFINITION FOR ACOUSTIC ELEMENTS IN CONTRACT WITH FLOOR
NSEL,S,LOC,X,0.24,5.24
NSEL,R,LOC,Y,0.24,4.24
NSEL,R,LOC,Z,0.24,0.49
ESLN,S,1
TYPE,1
MAT,2
EMODIF,ALL

!FLUID-STRUCTURE INTERFACE DEFINATION
NSEL,S,LOC,X,0.24
NSEL,R,LOC,Y,0.24,4.24
NSEL,R,LOC,Z,0.24,3.24
ESLN
SF,ALL,FSI
NSEL,S,LOC,X,5.24
NSEL,R,LOC,Y,0.24,4.24
NSEL,R,LOC,Z,0.24,3.24
ESLN
SF,ALL,FSI
NSEL,S,LOC,X,0.24,5.24
NSEL,R,LOC,Y,0.24
NSEL,R,LOC,Z,0.24,3.24
ESLN
SF,ALL,FSI
NSEL,S,LOC,X,0.24,5.24
NSEL,R,LOC,Y,4.24
NSEL,R,LOC,Z,0.24,3.24
ESLN
SF,ALL,FSI
NSEL,S,LOC,X,0.24,5.24
NSEL,R,LOC,Y,0.24,4.24
NSEL,R,LOC,Z,0.24
ESLN
SF,ALL,FSI
NSEL,S,LOC,X,0.24,5.24
NSEL,R,LOC,Y,0.24,4.24
NSEL,R,LOC,Z,3.24
ESLN
SF,ALL,FSI

FINISH

/SOLU
ANTYPE,3      ! Harmonic analysis

!Boundary conditions
!Four corners of floor
NSEL,S,LOC,X,0,0.24
NSEL,R,LOC,Y,0,0.24
NSEL,R,LOC,Z,0,0.24
d,all,uz,0
d,all,uy,0
d,all,ux,0

NSEL,S,LOC,X,5.24,5.48
NSEL,R,LOC,Y,0,0.24
NSEL,R,LOC,Z,0,0.24
d,all,uz,0
d,all,uy,0
d,all,ux,0

NSEL,S,LOC,X,5.245.48
NSEL,R,LOC,Y,4.24,4.48
NSEL,R,LOC,Z,0,0.24
d,all,uz,0
d,all,uy,0
d,all,ux,0

NSEL,S,LOC,X,0,0.24
NSEL,R,LOC,Y,4.24,4.48
NSEL,R,LOC,Z,0,0.24
d,all,uz,0
d,all,uy,0
d,all,ux,0

alls

!Four applied forces
F,6025,FZ,-500

HARFRQ,0,120,                       ! Frequency range
NSUBST,10,       ! Number of frequency steps
KBC,1        ! Stepped loads

SOLVE
FINISH

andrew19810403

看了你的程序，主要原因是你的网格划的太少，出来的结果应该是红彤彤的一片吧，应该保证一个波长上至少有八个网格，这样的结果才会接近真实情况。

anc

我算低频 ，80Hz 对应波长4m左右， 我最大的网格是0.25m 一个波长远不止 8个网格。问题我倒是解决了，fsi 设置错了。

lunion

我的情况跟你差不多，是液固界面，液体里面没有声压，fsi应该怎么设置？

anc wrote:
我算低频 ，80Hz 对应波长4m左右， 我最大的网格是0.25m 一个波长远不止 8个网格。问题我倒是解决了，fsi 设置错了。

yi1xiang

我看了也是，FSI你弄得多了点

zhangbit

液体里面没有声压,错了吧,应该只有声压

suiyuanshihua

其实FSI弄多点也没有错,只是计算量会增加,因为增加了节点的自由度啊,;
他的错,fsi没有作用引起的声压没有从结构（固体）传递到空气单元（流体）。

2208001017

好复杂的东西，大家有没有GUI的资料介绍一个

tzhdacy

fsi的单元需要FLUID30，就是流固耦合的界面单元体即可

magiciansf

初学者，路过学习一下！祝你早日解决问题！

tzhdacy

建议一下，如果作者解决了问题，可以把方法和心得说一下，最好把结果图上传一下。我为人人，人人为我，谢谢！