声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5026|回复: 22

[Virtual.Lab] sysnoise中有关ATV的算法问题,求教!

[复制链接]
发表于 2007-10-31 10:34 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
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

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2007-10-31 18:55 | 显示全部楼层
Sum(ATV×Vn)
发表于 2007-10-31 18:57 | 显示全部楼层
ATV、节点振速以及辐射声压是对应上的,你的过程可能有问题。
 楼主| 发表于 2007-11-1 15:40 | 显示全部楼层
ATV输出文件的格式是Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag吗?ATV云图只有real和imag选项,云图和数据文件格式之间是如何对应的?谢谢
发表于 2007-11-1 17:06 | 显示全部楼层

回复 #4 dlchen 的帖子

Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag分别为各方向投影,ATV云图中应到是法线方向的投影。
 楼主| 发表于 2007-11-1 20:02 | 显示全部楼层
针对别的例子也进行了验证(计算过程应该没有问题),发现由输出的ATV数据和已知的BC边界振速,通过sum(ATVreal×Vreal)不等于sysnoise计算得到的声压实部Preal。

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

感谢w89986581!!
发表于 2007-11-2 09:18 | 显示全部楼层
声压实部等于ATVreal*Vreal-ATVimag*Vimag.
 楼主| 发表于 2007-11-2 09:48 | 显示全部楼层
非常感谢w89986581!声压实部对上了。声压虚部等于?
发表于 2007-11-2 10:21 | 显示全部楼层
ATVreal*Vimag+ATVimag*Vreal
另外,辐射声压要取模,模表示振动幅度,而不是实部,呵呵。
 楼主| 发表于 2007-11-2 10:52 | 显示全部楼层
该问题在w89986581的不懈帮助下已圆满解决,可以封口了,呵呵。再次感谢w89986581 !!!
发表于 2007-11-2 14:02 | 显示全部楼层
不用客气,希望有好的学习经验分享,呵呵。
发表于 2012-6-7 11:22 | 显示全部楼层
Axreal、Aximag、Ayreal、Ayimag、Azreal、Azimag
谢谢
发表于 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指点一二?
发表于 2013-5-6 11:15 | 显示全部楼层
发表于 2013-5-6 11:20 | 显示全部楼层
w89986581 发表于 2007-11-2 09:18
声压实部等于ATVreal*Vreal-ATVimag*Vimag.

请问教授,为什么按照你说的公式计算对不上呢,Preal=ATVreal*Vreal-ATVimag*Vimag;
Pimag=ATVreal*Vimag+ATVimag*Vreal; newshopping也提出了这个问题
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-17 17:54 , Processed in 0.074941 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表