vl中频域声压值问题

mxlzhenzhu

cquhdd 发表于 2014-9-3 22:10
通过向 詹老师 求助，上述问题已经有点眉目了。后续会给出比较直观的对比。

希望解决问题后回来分享经验

cquhdd

按照詹老师的提示，我采用有限元法计算了一个2163立方米空混响室内的声场分布。具体参数设置如下：
1.   2163立方米的尺寸参数为：20.7m（高）×11m（长）×9.5（宽）。
2.   气体密度为1.251kg/m3. 氮气。声速344m/s。
3.   吸声边界条件处理。 设吸声系数为0.1， 根据法向吸声系数与声阻抗率的关系，计算得到声阻抗率为15816. 基于特征角的方法，在声学网格外表面建立了一组meshing group，并将其属性赋予为吸声板（absorption panel）。
4.  声源定义。声源采用单极子声源。设定了一个倍频程内的声源。倍频程内频率按照1HZ的步长进行等分，每个频率值上赋予的声功率为0.4694W。
5. 声学响应分析。
    1）计算参数设置。由于网格数比较大（网格尺寸为100mm），约为220万左右，因此在求解时采用了迭代法求解，参数设置为迭代步为30000， 收敛容差为10^(-9).
     2)  建立了两个平面场点。分别位于2m位置和6m位置。
     3） 采用线性扫频的方式计算。
     4） 结果处理。 在得到线性扫频的声学响应后，借助于后处理模块 vector processsing case， 将声学响应的结果处理为 倍频程内的结果。
6.  结果分析。
     下图分别是有限元法计算得到的2m平面和6m平面的中心频率为31.5Hz时的倍频程内声学响应值。
      可以看到在这种情况下，靠近房间中部位置的声压级差别都在5db以内，四周的声压级较大。（说明，声源布置在下图右下角的节点对应的房顶节点处），现在的问题是，靠近声源的这边（右边墙壁）区域出现了蓝色，即与左边的靠近墙壁部分的声压级分布不一致。这与解析法的结果不一致。解析法的结果完全是一种对称分布。解析法的结果如最后两张图。

     最后，这两种方法的结果为什么会不一致呢？有限元法的基本思路正确否？
     另，采用同样的思路，边界元法计算的结果与有限元的结果的声压分布的布局是一致的。

cquhdd

cquhdd 发表于 2014-9-6 11:43
按照詹老师的提示，我采用有限元法计算了一个2163立方米空混响室内的声场分布。具体参数设置如下：
1.   2 ...

自己先顶一下～～

cquhdd

居然没有人回复啊～{:{19}:}

lengxuef

自己没有做，不好说啊。只从描述上看，一点儿不成熟的思考，声源在右侧，左侧和右侧的声学响应会有些差别吧。。。

cquhdd

lengxuef 发表于 2014-9-7 21:02
自己没有做，不好说啊。只从描述上看，一点儿不成熟的思考，声源在右侧，左侧和右侧的声学响应会有些差别吧 ...

是这样子的，吸声系数较小的情况下解析法的左右差异不明显，较大的情况下，差异会比较直观。

cquhdd

lengxuef 发表于 2014-9-7 21:02
自己没有做，不好说啊。只从描述上看，一点儿不成熟的思考，声源在右侧，左侧和右侧的声学响应会有些差别吧 ...

再请教一个问题，vl中如何仿真消声室，按照《Virtual.Lab Acoustics声学仿真计算从入门到精通》封闭空间声场分布仿真分析的思路，把吸声系数设置为很大就可以了吗？

lengxuef

cquhdd 发表于 2014-9-7 22:38
再请教一个问题，vl中如何仿真消声室，按照《Virtual.Lab Acoustics声学仿真计算从入门到精通》封闭空间 ...

消声室把边界条件定义为无反射边界条件，貌似是有个数值416.5？记不清了，书上有的，也可以定义AML属性

cquhdd

lengxuef 发表于 2014-9-8 00:09
消声室把边界条件定义为无反射边界条件，貌似是有个数值416.5？记不清了，书上有的，也可以定义AML属性

好的，明天尝试一下。

cquhdd

按照上述的例子，将声阻抗率设置为0（实部和虚部都为0），计算的结果如下图，这该怎么解释啊？

cquhdd

cquhdd 发表于 2014-9-8 00:50
好的，明天尝试一下。

你好，这里的无反射边界条件是指什么啊？
PML是不是就是一种无反射的边界条件？

cquhdd

lengxuef 发表于 2014-9-8 00:09
消声室把边界条件定义为无反射边界条件，貌似是有个数值416.5？记不清了，书上有的，也可以定义AML属性

你好，这里的无反射边界条件是指什么啊？
我查看帮助文件后，看到PML认为是一种无反射边界（a non-reflecting boundary condition in FEM），是指这个吗？
还有，461.5是根据空气的常数来确定的哈，如果不同的气体媒介，这个数值是稍有变化的哈？

cquhdd

lengxuef 发表于 2014-7-29 23:40
点声源没法模拟混响室。

晕死，原来您就是阿伟老师啊!!!!

cquhdd

cquhdd 发表于 2014-9-8 02:40
晕死，原来您就是阿伟老师啊!!!!

按照您的建议，以及《Virtual.Lab Acoustics声学仿真计算高级应用实例(附光盘)》中封闭空间声场的计算流程，并将声阻抗率设置为416.5时，计算的不同平面的声场分布的结果如下。现在的有个疑问，在20m，18m平面的声场分布和自由声场很类似，但是在2m，6m，12m平面的声场分布就不像是自由声场的分布了。（个人理解，2m，6m，12m平面的声场分布也和20m，18m平面的声场分布类似，只是数值上有变化）。

cquhdd

lengxuef 发表于 2014-9-8 00:09
消声室把边界条件定义为无反射边界条件，貌似是有个数值416.5？记不清了，书上有的，也可以定义AML属性

再请教一下，要模拟上述的封闭空间的消声室时，AML属性合适吗?PML/AML只适合于外声场分析吗？如何可以进行内声场的话，参数如何设置啊？number of layer 以及 那个厚度都表示啥意义啊？