dlchen 发表于 2007-10-31 10:34

sysnoise中有关ATV的算法问题,求教!

sysnoise用ATV计算场点声压的算法是ATV×Vn=p;
但具体算法又是如何实现的呢?如下是用ATV加权系数实部×边界振速实部求和,得到场点声压实部吗?结果好像对不上?请教了!

已知某一场点对边界元模型的ATV加权矩阵如下:

ATV                                                                           
30-OCT-2007 17:12:16
      90      10         9
ATV            FREQUENCY          FIELD POINT
         1      0.10000000E+01         1
         1         1      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.97607214E-16      0.23840003E-04
         2         2      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.39042886E-15      0.95360012E-04
         3         3      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.97607214E-16      0.23840003E-04
         4         4      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.39042886E-15      0.95360012E-04
         5         5      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.15617154E-14      0.38144005E-03
         6         6      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.39042886E-15      0.95360012E-04
         7         7      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.97607214E-16      0.23840003E-04
         8         8      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.39042886E-15      0.95360012E-04
         9         9      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.97607214E-16      0.23840003E-04

又已知边界节点振速为:
FREQUENCY
      1.00000000E+00
NODAL VELOCITY VALUES
         1         1      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.53072114E+00      0.41095023E-01
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         2         2      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.21228846E+01      0.16438009E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         3         3      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.53072114E+00      0.41095023E-01
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         4         4      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.21228846E+01      0.16438009E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         5         5      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.84915382E+01      0.65752037E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         6         6      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.21228846E+01      0.16438009E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         7         7      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.53072114E+00      0.41095023E-01
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         8         8      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00    -0.21228846E+01      0.16438009E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
         9         9      0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00   -0.53072114E+00      0.41095023E-01
                        0.00000000E+00      0.00000000E+00      0.00000000E+00
                        0.00000000E+00      0.00000000E+00      0.00000000E+00

在sysnoise中由ATV计算得该场点的声压为:
FREQUENCY
      0.10000000E+01
PRESSURE VALUES
         1         1      0.87488287E-02   -0.10244495E-01      0.13471885E-01
                         -0.49502577E+02      0.53557668E+02      0.95260613E-02

w89986581 发表于 2007-10-31 18:55

Sum(ATV×Vn)

w89986581 发表于 2007-10-31 18:57

ATV、节点振速以及辐射声压是对应上的,你的过程可能有问题。

dlchen 发表于 2007-11-1 15:40

ATV输出文件的格式是Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag吗?ATV云图只有real和imag选项,云图和数据文件格式之间是如何对应的?谢谢

w89986581 发表于 2007-11-1 17:06

回复 #4 dlchen 的帖子

Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag分别为各方向投影,ATV云图中应到是法线方向的投影。

dlchen 发表于 2007-11-1 20:02

针对别的例子也进行了验证(计算过程应该没有问题),发现由输出的ATV数据和已知的BC边界振速,通过sum(ATVreal×Vreal)不等于sysnoise计算得到的声压实部Preal。

由上述例子也可以看出,ATV系数都是正数,边界振速实部都是负数,它们乘积求和不应该等于正数(而计算声压实部却为正)。看来sysnoise中利用ATV求声压的算法应该不是乘积求和,或者ATV输出数据有误?不明白啊!(因为需要ATV的信息,因此希望能搞明白sysnoise中ATV的真实含义)。望赐教!

感谢w89986581!!

w89986581 发表于 2007-11-2 09:18

声压实部等于ATVreal*Vreal-ATVimag*Vimag.

dlchen 发表于 2007-11-2 09:48

非常感谢w89986581!声压实部对上了。声压虚部等于?

w89986581 发表于 2007-11-2 10:21

ATVreal*Vimag+ATVimag*Vreal
另外,辐射声压要取模,模表示振动幅度,而不是实部,呵呵。

dlchen 发表于 2007-11-2 10:52

该问题在w89986581的不懈帮助下已圆满解决,可以封口了,呵呵。再次感谢w89986581 !!!

w89986581 发表于 2007-11-2 14:02

不用客气,希望有好的学习经验分享,呵呵。

benbendog 发表于 2012-6-7 11:22

Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag
谢谢

newshopping 发表于 2012-6-14 17:51

以上的这个例子我按照:
Preal=ATVreal*Vreal-ATVimag*Vimag;
Pimag=ATVreal*Vimag+ATVimag*Vreal;
进行计算和
FREQUENCY
      0.10000000E+01
PRESSURE VALUES
         1         1      0.87488287E-02   -0.10244495E-01      0.13471885E-01
                         -0.49502577E+02      0.53557668E+02      0.95260613E-02
对不上
请求LZ指点一二?

song0404 发表于 2013-5-6 11:15

newshopping 发表于 2012-6-14 17:51 static/image/common/back.gif
以上的这个例子我按照:
Preal=ATVreal*Vreal-ATVimag*Vimag;
Pimag=ATVreal*Vimag+ATVimag*Vreal;


确实对不上啊,怎么回事儿

song0404 发表于 2013-5-6 11:20

w89986581 发表于 2007-11-2 09:18 static/image/common/back.gif
声压实部等于ATVreal*Vreal-ATVimag*Vimag.

请问教授,为什么按照你说的公式计算对不上呢,Preal=ATVreal*Vreal-ATVimag*Vimag;
Pimag=ATVreal*Vimag+ATVimag*Vreal; newshopping也提出了这个问题
页: [1] 2
查看完整版本: sysnoise中有关ATV的算法问题,求教!