结构辐射声功率和辐射效率问题

榫眼

                对四边简支矩形钢板new_ban.out进行结构辐射声功率和辐射效率计算
计算步骤：
      
          1.先在nastran中用模态频率响应法计算出结构模型的动力学响应（位移文件为new_ban_example_1.pch）
          2.在sysnoise中采用间接边界元法计算结构声辐射，把位移文件为new_ban_example_1.pch导人sysnoise中作为边界条件.
      
          3.在频率5-800hz间计算
Sysnoise计算代码(请大家看看计算过程有没有错？）
New Name 'bem_example' Model 1 File model1.sdb Return
Option
BEM Indirect Variational Uncoupled Unbaffled Frequency
Return
Import Mesh Format Patran File yroom\new_ban.out Return
Material Fluid

Sound Real 343
Rho Real 1.21

Return
Set 1 Name
"Envelope" Envelope

Elements All

Return
Boundary Jump Pressure
Real 0 Imag 0

Nodes Set 1

Return
Generate

Element All

From Displacements File
yroom\new_ban_example_1.pch Format Nastran

Frequency 5 To 800 LinStep 5

Norewind

Mesh File
yroom\new_ban.out Format Patran

Return
Solve

Frequency 5 To 800 LinStep 5

Return

榫眼

4.在sysnoise中不加自由边界条件和速度跳动量边界条件,其他如上保持不变
New Name 'bem_example' Model 1 File model1.sdb Return
Option
BEM Indirect Variational Uncoupled Unbaffled Frequency
Return
Import Mesh Format Patran File yroom\new_ban.out Return
Material Fluid

Sound Real 343
Rho Real 1.21

Return
Set 1 Name
"Envelope" Envelope

Elements All

Return
Generate

Element All

From Displacements File
yroom\new_ban_example_mode.pch Format Nastran

Frequency 5 To 800 LinStep 5

Norewind

Mesh File
yroom\new_ban.out Format Patran

Return
Solve

Frequency 5 To 800 LinStep 10

Return

榫眼

对嵌在无限挡板上的板状结构的辐射声功率可由原场Rayleigh积分求得：
            

榫眼

第一种方法加了边界条件，计算结果（输出辐射声功率实部）与理论值（引用某论文结果）相比，我计算出来的声功率，峰值在350HZ最大，而且输出声功率出现负值，不知是什么原因造成的？请大家指点一下！
第二种方法没加任何边界条件，输出辐射声功率实部与理论结果曲线有点类似，即在低频时辐射功率最大。
          请大家指教一下！谢谢

lyxh

你的计算中没有障板，当然和解析解不同。

lyxh

第二种方法等于在边界上加了载荷。实际情况中很难实现这种边界条件。

榫眼

那如何设立障板呢？请lyxh指点一下

榫眼

非常感谢lyxh的指点和解说！！！！！！！！！！！！

榫眼

1.在直接边界元分析模型中加了障板(baffled)和设置边界速度跳动量为0,以及设置了两场点a(0,0,0.2),b(0,0,-0.2) 分别置于障板的上下面,
结果与未加障板和未设置边界速度跳动量的结果一样,并没有任何改变,这是什么原因呢,而且测得两场点的声压值都
为0,


2.在直接边界元分析模型中加了障板(baffled)和设置边界声压跳动量为0,以及设置了两场点a(0,0,0.2),b(0,0,-0.2) 分别置于障板的上下面,
结果与未加障板和设置边界声压跳动量的结果一样,并没有任何改变,而且测得两场点的声压也都值都为0,
为什么设置的障板都没起作用呢?

pengweicai

应该是  间接边界元吧。

场点的计算 进行后处理没有。

模型的网格是不是在 障板的平面内啊？

榫眼

上面写错了，是间接边界元，模型的网格嵌在障板上，也进行过场点的计算，可是提取场点的声压值全都是0

榫眼

New Name 'bem_example' Model 1 File model1.sdb Return
Option  BEM Indirect Variational Uncoupled baffled Frequency  Return
Import Mesh Format Patran File yroom\new_ban.out Return
Material Fluid
    Sound Real 343  Rho Real 1.21
    Return
Set 1 Name  "Envelope" Envelope
    Elements All
    Return
Boundary Jump Pressure  Real 0 Imag 0
    Nodes Set 1
    Return
Point Line  0 0 0.2 To 0 0 -0.2 Divide 1 Return
Generate
    Element All
    From Displacements File  yroom\new_ban_example_mode.pch Format Nastran
    Frequency 5 To 800 LinStep 5
    Norewind
    Mesh File  yroom\new_ban.out Format Patran
    Return
Solve
    Frequency 5 To 800 LinStep 5
    Return   
PostProcess
    Points All
    Frequency 5 To 800 LinStep 5
    Near 2
    Far 5
    Quadrature 2 2 1
    Save Results Step 1
    Return
Store
Case Active
Results
    Frequency 5 To 800 LinStep 5
    Points All
  Return

Return

maxumoon

我也要做嵌在障板上的有限长圆柱壳的声辐射，理论解已经解出来了，想验证一下。LZ的问题解决了吗

w89986581

sysnoise是无法模拟无限长圆柱障板的,它只提供了无限大平面障板.你只能取有限长圆柱障板试试啦,呵呵.
提示:你的理论解应该有适用频率范围的.

yeyecoming

我想自己做个盒子，就挡住了往里面辐射的声波。这回就像个障板的作用了，然后只计算外声场。结果页没看出声辐射的结果和理论值有啥关系。一点都不一样