在VL中数值计算的形式模拟阻抗管法测量样品的隔声量问题

年少才轻狂

我研究很长时间关于在VL中以数值计算的形式模拟阻抗管法测量样品隔声量的问题，但是还是没能解决问题，故特此再发一帖，寻求各位的指导和帮助，谢谢。

附：附件中有有关于这个阻抗管法的文献，我的仿真步骤就是参考文献进行。

模型：80x80x2mm铝板，采用自由边界条件，L1=L2=100mm，S1=S2=80mm。计算平面波垂直入射时样品材料的隔声量（传声损失）。铝板材料属性：密度：2730，弹性模量：7.76E10，泊松比：0.35192。

结合文献中提到的步骤，我用Virtual.lab仿真计算铝板的隔声量，但是仿真结果和文献的对不上，不知道问题出在哪里，现特把步骤列出，希望能帮我看出问题出在哪里。

步骤如下：

1）在ANSYS中建立被测样品的几何模型和驻波管流体几何模型（驻波前、后管的尺寸都是80x80x300mm），设置样品的材料属性，并进行划分网格（后面附有图形）；

2）将样品的材网格模型和驻波管前、后流体网格模型导入到Virtual.lab中，定义流体模型的材料属性（定义为空气）；

3）在流体模型中心轴线上分别定义四个场点：1（0,0,-181), 2（0,0,-101）,3（0,0,101）,4（0,0,181）.这四个场点分别代表四个传声器的位置；

4）设置驻波后管的尾端面为全吸声（Impedance（声阻抗），实部位416.5，虚部为0）；

5）定义耦合面（creat coupling surface）：分别定义样品模型与流体模型接触的两个面为耦合面，容差（tolerance)为5mm；

6）在驻波前管的端面定义平面波声源（0,0，-301），实部为1，虚部为0；

7）用Virtual.lab的FEM中的直接声振耦合求解驻波管中声场分布情况(0~2000Hz)；

8）分别读取4个场点的声压数据，导出到外部声压文件中。

9）在MATLAB中编辑公式1（P1 P2 P3 P4分别表示1 2 3 4个场点处的声压值），然后导入四场点声压值得到传递损失图。

年少才轻狂

之前忘记补充模型图，先截图如下：

superxjw

你中间结构网格给的什么约束？

lengxuef

模型传到百度网盘，然后把模型地址在这里发一下
另外，你说对不上，你把两者对比的曲线放上来。

年少才轻狂

lengxuef 发表于 2013-12-25 21:56
模型传到百度网盘，然后把模型地址在这里发一下
另外，你说对不上，你把两者对比的曲线放上来。 ...

我把网格模型，计算结果数据，以及我自己仿真过程的录像视频都放在百度网盘里面了，地址是http://pan.baidu.com/disk/home#dir，然后进入“在Virtual.lab中以数值计算的形式模拟阻抗管法测量样品的隔声量（传声损失）”，就可以看到所有数据。
附：这两种图片中，图片1是我用VL计算得到的隔声量曲线图，，图片2 是文献用Sysnoise计算得到的隔声量曲线图。我计算得到的和文献计算得到的完全不同，求指导？谢谢

年少才轻狂

年少才轻狂 发表于 2013-12-26 10:00
我把网格模型，计算结果数据，以及我自己仿真过程的录像视频都放在百度网盘里面了，地址是http://pan.bai ...

刚刚有点问题，我之前在图中指的是细实线，其实是粗实线代表80*80*2mm的隔声量。另外，粗虚线表示2mm无限大铝板的隔声量，细实线代表用经验公式计算得到的隔声量曲线。

年少才轻狂

superxjw 发表于 2013-12-25 21:26
你中间结构网格给的什么约束？

我对结构网格没有施加任何约束，也就是自由边界。因为我看到文献里是采用自由边界计算的铝板模态。

lengxuef

年少才轻狂 发表于 2013-12-26 10:00
我把网格模型，计算结果数据，以及我自己仿真过程的录像视频都放在百度网盘里面了，地址是http://pan.bai ...

网盘地址不对。。

年少才轻狂

lengxuef 发表于 2013-12-26 10:58
网盘地址不对。。

现在好了，谢谢。请进入http://pan.baidu.com/s/1jGr3FKi

superxjw

年少才轻狂 发表于 2013-12-26 10:13
我对结构网格没有施加任何约束，也就是自由边界。因为我看到文献里是采用自由边界计算的铝板模态。 ...

这种实验，试验件都是要夹持的，自由模态差距会很大的。。。

年少才轻狂

superxjw 发表于 2013-12-26 18:33
这种实验，试验件都是要夹持的，自由模态差距会很大的。。。

你好，我后来也算了采取约束的情况，但是计算结果显示，曲线走势一样，只是幅值的差别。