没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
为了研究水下工程结构辐射声场的识别问题,基于相位共轭方法采用离散阵元对水下有限长圆柱壳的辐射声场识别进行了数值仿真计算.采用不同的阵列形式分别在近场和远场对圆柱壳表面声压、法向速度、辐射声功率和声强进行了识别及重构误差分析,并讨论了激励频率变化对重构结果的影响.数值计算结果表明:只有在近场采用圆柱形阵列形式并基于声压梯度测量,使用偶极子源才能得到最好的声场重构结果;采用相位共轭方法可以较好地识别结构的表面声压、法向速度、辐射声功率和声强等;在远场使用相位共轭方法仅能够识别出圆柱壳表面声压最大值的位置;激励
资源推荐
资源详情
资源评论
书书书
第 卷第 期
哈 尔 滨 工 程 大 学 学 报
年 月
:
相位共轭法进行水下圆柱壳辐射声场的识别研究
刘松
,黎胜
,赵德有
( 大连理工大学 工业装备结构分析国家重点实验室,辽宁 大连 ; 大连理工大学 船舶工程学院,辽宁 大连 )
摘 要:为了研究水下工程结构辐射声场的识别问题,基于相位共轭方法采用离散阵元对水下有限长圆柱壳的辐射声场
识别进行了数值仿真计算 采用不同的阵列形式分别在近场和远场对圆柱壳表面声压、法向速度、辐射声功率和声强进
行了识别及重构误差分析,并讨论了激励频率变化对重构结果的影响 数值计算结果表明:只有在近场采用圆柱形阵列
形式并基于声压梯度测量,使用偶极子源才能得到最好的声场重构结果;采用相位共轭方法可以较好地识别结构的表面
声压、法向速度、辐射声功率和声强等;在远场使用相位共轭方法仅能够识别出圆柱壳表面声压最大值的位置;激励频率
发生变化时该方法仍能得到较精确的重构结果
相位共轭方法具有使用测点少,可自然使用局部重建方法,并且不存在
解不稳定问题等优点,具有很好的应用前景
关键词:相位共轭方法;水下圆柱壳;辐射声场;识别
中图分类号: 文献标识码: 文章编号:()
Study on the identification of radiation sound field from
underwater cylindrical shell by phase conjugation method
,
,
( , , ,;
, , ,)
Abstract:
, , ,
,
Keywords: ; ; ;
收稿日期:
基金项目:国家自然科学基金资助项目()
作者简介:刘松(),男,博士研究生,: ;
黎胜(),男,副教授,博士生导师
通信作者:刘松
水下弹性圆柱壳的振动声辐射问题一直是各国
学者广泛关注的研究课题,而对于水下弹性圆柱壳
辐射声场的重构问题研究的相对较少
目前对辐射
声场重构主要是近场声全息(
,)方法,很多学者
[]
围绕这个方法展
开了广泛的探讨
本文采用相位共轭(时间反转)方
法对水下圆柱壳的辐射声场进行识别 相位 共 轭
(时间反转)方法一直被广泛应用于光学和水声通
讯中,虽然国内外对该方法很关注,但在声学领域的
应用都只局限于对点声源的研究,对于将该方法用
于实际工程结构的声场重构却没有做深入的研究
时间反转方法实际上是利用了线性化波动方程中只
包含声压对时间的二阶导数这个特点 但时间反转
法用于声源识别和定位的能力受到声波衍射极限分
辨率的限制,理论上其分辨率最高为
波长大小,
因此并不能提供出足够的噪声源位置和噪声辐射特
性信息,一些学者对解决时间反转法中衍射极限的
方法
[]
,本文基于相位共轭方法采用离散阵元对
无限水深中的有限长圆柱壳的辐射声场识别和重构
进行了数值仿真计算,采用不同的阵列形式分别在
近场和远场对圆柱壳表面声压、法向速度、辐射声功
率和声强进行了识别及重构误差分析
结构声辐射和相位共轭理论
结构声辐射理论
结构在简谐力作用下考虑流体加载效应的有限
元形式运动方程为
[]
Zv = f
e
+ f
p
. ()
式中:Z (
ω
M
ω
C K)/ i
ω
为结构阻抗矩阵,
M、C 和 K 分别为结构质量矩阵、阻尼矩阵和刚度矩
阵,
ω
为激励圆频率;v 为结构速度向量;f
e
为外激
励力向量;f
p
为结构表面声压所引起的流体对结构
的作用力向量
结构振动声辐射的边界元法计算公式为
Ep = Dv
n
. ()
式中:E 和 D 阵为 N N 维的复系数影响矩阵,N 为
结构表面离散节点个数,p 为结构表面声压向量,v
n
为结构表面法向速度向量
由 v
n
和结构速度向量 v 之间的转换关系式
v
n
= G
v ()
以及 f
p
和 p 的关系式
f
p
= - GAp. ()
式中:G 为方向余弦转换矩阵,A
∫
N
NS,N 为形
状函数矩阵,可得求解结构 声耦合问题的方程为
(
Z + GAE
-
DG
)v = f
e
. ()
令
GAE
DG
C
A
(
ω
)
ω
M
A
(
ω
),则:M
A
和
C
A
(
ω
)分别为流体附加质量阵和流体阻尼阵
[]
. 将
其代入式(),即可得结构 声耦合自由振动方程,
求解其广义特征值问题,可得水下结构声固耦合的
固有频率.
按照式()求出结构速度向量 v 后,进而可求
出
v
n
和 p 等.
结构的辐射声功率和声强按照下式计算:
W =
∫
S
(pv
*
n
)S, ()
I =
(pv
*
n
). ()
式中:v
*
n
为结构法向振速的共轭,p 为结构表面声
压,( )表示取实部
相位共轭理论
线性化的齐次声学波动方程为
▽
p(r,t)-
c
a
p(r,t)
a
t
= . ()
式中:p(r,t)为声压,c 为流体介质中的声速. 时间反
转方法可实现声波的反向传播和自适应聚焦的原因
就在于线性化波动方程()中只包含声压对时间的
二阶导数,这样 p(r,t)和 p( r, t)都是波动方程的
解,其中 p(r,t)代表了由声源向外辐射正向传播的
声场,而 p( r, t)代表了反向传播(即向声源区传
播并汇聚于声源的声波)的声场. 时域中的时间反
转等价于频域中的相位共轭
[]
,即 p( r,t)和 p( r,
t)等价于 p( r,
ω
)和 p
*
( r,
ω
),此处 p
*
( r,
ω
)为
p(r,
ω
)的复共轭. 同样 p(r,
ω
)为正向传播的声场,
p
*
(r,
ω
)为反向传播的声场.
根据 积分方程,相位共轭阵
列测量到声压和声压梯度后,相位共轭声场为
[]
p
PC / P
(r,r
s
)=
∫
G(r,r')
a
p
*
(r',r
s
)
a
n
- p
*
(r',r
s
)
a
G(r,r')
a
[ ]
n
S'.
()
式中:
r
s
为源点位置,r 为场点位置,r'为阵列中阵元
的位置,
G(r,r')
kR
π
R
为自由场格林函数,R r
r' ,S 为闭合阵列表面积,n 为 S 的内法向单位矢
量,P
*
和
a
P
*
a
n
分别为相位共轭阵列测量得到的声压
和声压梯度
对包含 N 个阵元的离散的有限阵列,考虑阵元
面积,其相位共轭声场为
p
PC / P
(r,r
s
)=
∑
N
n =
G(r,r
n
)
a
p
*
(r
n
,r
s
)
a
n
- p
*
(r
n
,r
s
)
a
G(r,r
n
)
a
( )
n
s
[ ]
n
.
()
式中:r
n
为离散阵列中每个阵元的位置 离散的阵元
间距至少需要满足小于 波长以得到较高的重构
分辨率
还可采用以下 种相位共轭阵列进行识别和重
构 如基于测量声压使用单极子源发射,其相位共轭
声场为
p
PC / M
(r,r
s
)=
∑
N
n =
[(G(r,r
n
)p
*
(r
n
,r
s
))× S
n
].
()
如基于测量声压梯度使用偶极子源发射,其相
··
哈 尔 滨 工 程 大 学 学 报 第 卷
剩余6页未读,继续阅读
资源评论
weixin_38631331
- 粉丝: 5
- 资源: 907
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功