点源激励下的球壳水中的声辐射！

qiuchanglin

计算模型为以钢制球壳，半径1m，厚度1cm。中间充满空气，在球心处有一幅值为1的点源作用，球壳外面为水。
由于sysnoise不能处理双面耦合问题，基本思路为分步计算，如下：
1. 采用DBEM 设置内部流体为空气，然后加点源，算得球壳内表面的声压。
2. 采用DBEM COUPLE 设置外部流体为水，然后将第一步的声压作为压力边界条件加载计算。
不知上面的思路是否可行？？

    在计算第一步时，采用DBEM COUPLE 和DBEM UNCOUPLE 算得的壳体表面的声压一样。这和分析是一样的，因为内部是轻流体，因而是否耦合对于计算结果影响不大。还有另外一种方法就是第一步中把空气模型用有限元模型建出来，采用流体有限元和结构有限元耦合的方法计算得到内表面的声压。计算得到的壳体内表面处的声压和采用DBEM所得结果不一样，而且相差很大，请问这是什么原因啊？应该采用那种方法计算的结果更可信些？

w89986581

"还有另外一种方法就是第一步中把空气模型用有限元模型建出来，采用流体有限元和结构有限元耦合的方法计算得到内表面的声压。"这样做则忽略了水的耦合作用。

鉴于问题的特殊性，即中心对称，我建议设置阻抗边界条件，即脉动球的辐射阻抗，免去分析流固耦合。

Timoshenko

你的方法基本上是可以的，但是缺乏验证！

qiuchanglin

对于另外一种方法我的意思是用流体有限元（内部空气）和结构有限元耦合计算得到内表面的声压，然后再用直接边界元耦合的方法，将上面的声压作为边界条件导入然后计算球壳外的声场。 即在计算球壳内表面声压的时候采用了两种不同的方法：1.DBEM  2.流体有限元加结构有限元，结果发现这两种方法计算得到的内表面的声压完全不一样。相差很大。

w89986581

好久不见。

楼主：你的问题可以转化为一维问题，同时脉动球辐射阻抗是有解析解的，因此辐射声场也是有解析解的，呵呵。此例用来验证sysnoise分析的准确性再合适不过了。

个人意见，欢迎讨论。

qiuchanglin

确实如此，这个就是用来验证SYSNOISE的正确性，等我摸索出结果，贴上来。

w89986581

好的，欢迎！

qiuchanglin

这样的结果不知道算不算还行，前一段时间做的，由于比较忙没有贴上来。理论解见W u T W , DandapaniA. A boundary element solution
for sound t ransm ission th rough th in panels [ J ]. J　Sound V ib, 1994, 171: 145～ 157.　
由于理论解没有阻尼，在250ｈｚ附近趋向无穷大，ｓｙｓｎｏｉｓｅ中设置了模态阻尼为0.01。

qiuchanglin

理论公式

pengweicai

内场计算出的声压，有几千帕？？？

  还有 能对结果做些分析吗？  比如外场的声压 相比内场的要小不少啊。

qiuchanglin

确实有几千帕，可以参考附件结果也是这么大
内表面的声压明显大于外表面的声压。
球壳内表面的声压包括：直达声，散射声（反射和球壳振动向内部辐射），SYSNOISE计算内场时内部是空气，没有考虑耦合作用即认为结构内表面刚性，即向内部辐射的部分应该为零
球壳外表面的声压应该只有结构振动向外辐射的部分

[ 本帖最后由 qiuchanglin 于 2007-8-22 15:10 编辑 ]

realhappysong

没明白楼主到底是用了什么样的方法，在哪里设置的模态阻尼？
难道是采用流体有限元和结构有限元分析的内部声场，在耦合前计算了声模态与结构模态吗？

oydh2010

oydh2010

　　请勿讨论D版的下载与安装问题

[ 本帖最后由 w89986581 于 2008-4-12 20:02 编辑 ]