两端开口的圆管声场是否也会存在非稳定频率的问题？

gazing

1。模型是一个两端开口的圆管，内有一系列分布点声源。管壁完全刚性。

2。分布点源沿一垂直于轴向的圆周分布，模拟一个旋转点声源。

3。用间接边界元方法，在两端自由边上设置自由边界条件。

4。计算结果，会在2阶谐频（BPFs）上比文献仿真结果大2倍左右，1阶和3阶都符合的比较好。

5。已经同样用间接边界元方法，在远场设一个很小的板来近似自由空间的情况，所得结果是符合理论解的，说明声源的设置没有问题。

对管道网格的划分满足小于1/6波长，尝试过减小网格和提高积分精度到3，3，2，结果变化都不大。

想到间接边界元存在非稳定频率的问题，但是按照书上说的，是封闭声腔的问题，在两端开口的圆管中都否也会存在？如果不存在，那可能会是什么方面的问题？请高手指点。

gazing

补充一下，点源是指点力源，相当于偶极子，在sysnoise5.5中我是用两个距离很近相位相反的单极子来模拟的。

pengweicai

对你的声源不是太懂。 最好把模型贴上来。

gazing

附件：1。两端开口连通自由空间的圆管。设置自由边边界条件。用间接边界元这样直接做开口圆管，应该没有什么问题吧？对Sysnoise的反复无常有点敬畏了。。。。
           2。管道内中部有声源，为旋转点力源。用一系列分布偶极子声源模拟。
           3。管道中部的分布声源。每一个点力源（偶极子）在sysnoise中用相位相反幅值相同的两单极子模拟。具体的声源设置，一时说不清楚，总之是相当于旋转的偶极子源，用read从外部导入命令流。

除去声源设置，整个模型的命令流很简单：

New Name 'test' Model 1 File test.sdb Return
Environment Section SETUP ANSYSREVISION '5.7' Return
Import Mesh Format Ansys File roundduct.cdb Return
Check Mesh Overdetermination Return
Coarse Skin Fix Return
Material Fluid
    Name 'air'
    Sound Real 3.4000e+002  Rho Real 1.2250e+000
    Return
Set Name  "Envelope" Envelope Return
Boundary Jump Pressure  Real 0 Imag 0
    Nodes Set 1
    Return
Read File  diprot1029.cmd
Table 1 —————————————————————————————以下是导入的具体的声源设置，省略
Type Frequency Name st1
0.000000 Magnitude 1.554247 Phase 0.000000
9.010000 Magnitude 3.108417 Phase -1.406250
18.020000 Magnitude 3.108183 Phase -2.812500
27.030000 Magnitude 3.107793 Phase -4.218750
Return
Source Name s1 Spherical
Amplitude Table 1
Position 0.499849 0.012271 0.000100
Return

Table 2
Type Frequency Name st2
0.000000 Magnitude 1.554247 Phase 0.000000
9.010000 Magnitude 3.108417 Phase -2.812500
18.020000 Magnitude 3.108183 Phase -5.625000
27.030000 Magnitude 3.107793 Phase -8.437500
Return
Source Name s2 Spherical
Amplitude Table 2
Position 0.499398 0.024534 0.000100
Return

Table 3
Type Frequency Name st3
0.000000 Magnitude 1.554247 Phase 0.000000
9.010000 Magnitude 3.108417 Phase -4.218750
18.020000 Magnitude 3.108183 Phase -8.437500
27.030000 Magnitude 3.107793 Phase -12.656250
Return
Source Name s3 Spherical
Amplitude Table 3
Position 0.498645 0.036782 0.000100
Return

。。。。。。。。。。。