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的算法问题，求教！