带障板的活塞声辐射

pengweicai

分析类型为 直接边界元


New Model 1 File model1.sdb Return
Option  BEM Direct Collocation Node Exterior Uncoupled Baffled NoAsymptotic Frequency  Return
Import Mesh Format Free File piston.fre Return
Boundary Velocity  Real .01 Imag 0
    Nodes All
    Return
Solve
    Frequency 1000
    Return






网格模型文件  piston.txt  后缀改为 .fre


2k.gif为2000HZ时场点声压值。

[ 本帖最后由 pengweicai 于 2008-7-21 18:37 编辑 ]

w89986581

命令流见附件test1.txt

榫眼

好东东！正需要，我对彭老大的感激之情，如滔滔之水，绵绵不绝啊！

minetree

昨天看了 w89986581的IBEM的例子，很受启发。在这里再深入讨论下。

仔细想了一下，无限大障板上的活塞声辐射（理想情况下是无厚度的膜的单侧辐射），从物理上来看，其空间的Green函数应该是源及其无限大板的镜像源的两个自由场的叠加。从这个出发，计算时，仅需要定义激励条件：活塞单元正面和反面的振动速度始终相等但反向就可以了。

该条件在Sysnoise中的IBEM Option下可以定义，具体方法是定义单元正面速度为v，负面速度为-v；这里单元负面速度的负号表明该速度与软件中定义的单元的正向（正面的外法向）相反，从而与单元正面的速度完全相反。

这样实现出来的总声场，当活塞的半径a/波长很小时，远场就相当于是两个靠的很近的同相点源的总辐射，即++模式；

同样原理，这时如果把单元负面速度值赋为v，即与正面速度同向等值，会怎么样呢？呵呵，想一想，其实就是真实的单一活塞在自由空间的辐射场，相当于IBEM中的Continuous速度激励条件。
这时当活塞的半径a/波长很小时，远场就相当于两个靠的很近的反相点源的总场，即+-模式，呵呵，就是偶级子了：）

为了方便对比，我把对比图发上来，左边的是无限大障板的活塞辐射（同相点源），右边是自由场单活塞辐射（偶级子），可以发现，右边中间垂直于轴线的面的声压为极小值0，总体声压为8字形分布，这都符合偶级子的特征。
而左边呢，P的分布在上下空间是完全对称的，也就是镜像的概念。这个和户外声传播中的无限大地面，是一回事。


相关的codes:

New Name 'model1_PVpNVn' Model 1 File model1_PVpNVn.sdb Return
Option BEM Indirect Frequency Return
Import Mesh Format Free File piston.fre Return
Set 1 Name  "Envelope" Envelope
    Elements all
    Return
Boundary Jump Pressure  Real 0 Imag 0
    Nodes Set 1
    Return
Boundary Velocity  Real 0.01 Imag 0 Positive
    Elements all
    Return
Boundary Velocity  Real -0.01 Imag 0 Negative
    Elements all
    Return
Save Return
Solve
    Frequency 1000
    Return
Point Plane  0 -4 4 To 0 -4 -4 Divide 50 To 0 4 4 Divide 50 Return
Save File .\IBEM_VBCsMdl1PVpNVn.sdb Return
PostProcess
    Points all
    Frequency 1000
    Near 2
    Far 5
    Quadrature 2 2 1
    Save Results Step 1
    Return

Save Return
close Return

New Name 'model2_PVpNVp' Model 1 File model2_PVpNVp.sdb Return
Option BEM Indirect Frequency Return
Import Mesh Format Free File piston.fre Return
Set 1 Name  "Envelope" Envelope
    Elements all
    Return
Boundary Jump Pressure  Real 0 Imag 0
    Nodes Set 1
    Return
Boundary Velocity  Real 0.01 Imag 0 Positive
    Elements all
    Return
Boundary Velocity  Real 0.01 Imag 0 Negative
    Elements all
    Return
Save Return
Solve
    Frequency 1000
    Return
Point Plane  0 -4 4 To 0 -4 -4 Divide 50 To 0 4 4 Divide 50 Return
Save File .\IBEM_VBCsMdl2PVpNVp.sdb Return
PostProcess
    Points all
    Frequency 1000
    Near 2
    Far 5
    Quadrature 2 2 1
    Save Results Step 1
    Return

Save Return
close Return

[ 本帖最后由 minetree 于 2008-7-22 15:37 编辑 ]

nmdb

peng主任您好，看到您的这个算例，有一个地方没有弄明白，就是带障板的活塞声辐射，按照边界元理论应该是间接边界元啊，rayleigh可以直接搞定。可是在计算的时候用直接边界元也可以。这为何啊？
直接边界元处理一般的连续的封闭曲面声辐射问题，可是对于此类不封闭的问题也搞的定，请赐教啊。

[ 本帖最后由 nmdb 于 2009-3-28 11:28 编辑 ]

pengweicai

直接边界元 只有对障板的处理可以不是封闭的。

rayleigh积分只是近似的算法， 有些误差。。

边界元计算精度高。。

hhhjjjhwc

不好意思问下你的文件里面怎么没有分析类型的选项啊？直接还是间接？还是有限，还是耦合？

w89986581

没有设置类型的话,sysnoise默认为间接边界元.

ryxc271

原帖由 pengweicai 于 2008-7-21 18:36 发表

                               
登录/注册后可看大图

分析类型为 直接边界元


New Model 1 File model1.sdb Return
Option  BEM Direct Collocation Node Exterior Uncoupled Baffled NoAsymptotic Frequency  Return
Import Mesh Format Free File piston.fre ...

eng



彭老大，我按你的题目要求做了，发现声压级的数量级小的问题，请指导，谢谢。

hhhjjjhwc

楼上的单位是不是弄成mm了？

ryxc271

hhhjjjhwc谢谢你的回复，我网格模型使用顶楼彭老大提供的，没有改变直接导入的，另外，后来我按比例放大网格模型，但是声压级的数量级还是没有改变，不晓得是是什么原因，谢谢高手来指点。

wenzhuhust

好好学习下，这些东西值得好好消化:lol :lol :lol

西瓜

下来看看啊

lhpvibration

非常及时！谢谢

cyh251

学习了,有灌水之嫌,勿怪