mxlzhenzhu 发表于 2014-9-4 23:46

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

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

cquhdd 发表于 2014-9-6 11:43

本帖最后由 cquhdd 于 2014-9-6 11:59 编辑

按照詹老师的提示,我采用有限元法计算了一个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 发表于 2014-9-6 12:06

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

自己先顶一下~~

cquhdd 发表于 2014-9-7 10:51

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

lengxuef 发表于 2014-9-7 21:02

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

cquhdd 发表于 2014-9-7 22:32

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

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

cquhdd 发表于 2014-9-7 22:38

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

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

lengxuef 发表于 2014-9-8 00:09

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

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

cquhdd 发表于 2014-9-8 00:50

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

好的,明天尝试一下。

cquhdd 发表于 2014-9-8 02:09


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

cquhdd 发表于 2014-9-8 02:28

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

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

cquhdd 发表于 2014-9-8 02:28

本帖最后由 cquhdd 于 2014-9-8 03:03 编辑

lengxuef 发表于 2014-9-8 00:09
消声室把边界条件定义为无反射边界条件,貌似是有个数值416.5?记不清了,书上有的,也可以定义AML属性
你好,这里的无反射边界条件是指什么啊?
我查看帮助文件后,看到PML认为是一种无反射边界(a non-reflecting boundary condition in FEM),是指这个吗?
还有,461.5是根据空气的常数来确定的哈,如果不同的气体媒介,这个数值是稍有变化的哈?

cquhdd 发表于 2014-9-8 02:40

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

晕死,原来您就是阿伟老师啊!!!!{:{05}:}

cquhdd 发表于 2014-9-8 11:13

本帖最后由 cquhdd 于 2014-9-8 16:10 编辑

cquhdd 发表于 2014-9-8 02:40
晕死,原来您就是阿伟老师啊!!!!
按照您的建议,以及《Virtual.Lab Acoustics声学仿真计算高级应用实例(附光盘)》中封闭空间声场的计算流程,并将声阻抗率设置为416.5时,计算的不同平面的声场分布的结果如下。现在的有个疑问,在20m,18m平面的声场分布和自由声场很类似,但是在2m,6m,12m平面的声场分布就不像是自由声场的分布了。(个人理解,2m,6m,12m平面的声场分布也和20m,18m平面的声场分布类似,只是数值上有变化)。

cquhdd 发表于 2014-9-8 13:18

本帖最后由 cquhdd 于 2014-9-8 13:42 编辑

lengxuef 发表于 2014-9-8 00:09
消声室把边界条件定义为无反射边界条件,貌似是有个数值416.5?记不清了,书上有的,也可以定义AML属性
再请教一下,要模拟上述的封闭空间的消声室时,AML属性合适吗?PML/AML只适合于外声场分析吗?如何可以进行内声场的话,参数如何设置啊?number of layer 以及 那个厚度都表示啥意义啊?
页: 1 [2] 3
查看完整版本: vl中频域声压值问题