【技术实现步骤摘要】
一种基于结构动刚度和声学边界元耦合的声场计算方法
[0001]本专利技术涉及声场预测
,尤其是一种基于结构动刚度和声学边界元耦合的声场计算方法。
技术介绍
[0002]大型船舶设备激励源众多、频谱特征各异,其诱发的结构振动和噪声既包含丰富的低频线谱成分,又有宽带连续谱成分,因此全频段振动噪声预报技术成为船舶结构声学设计各个阶段的现实需求。遗憾的是,现有成熟的数值技术,如有限元、边界元、统计能量分析等,因离散方式、分析变量等不同仅适用于特定频段。近年来发展起来的动刚度法,由于其位移函数采用解析解,且离散误差不依赖于频率,有望成为全频段振动噪声预报的强有力潜在手段。
[0003]边界元法是继有限元法之后的另一种高效实用的数值方法,特别适用于无限域和半无限域内流固耦合系统动力学预报问题。其将无限域上的三维微分方程转化为有限边界上的二维积分方程,通过在边界面上的离散实现维数下降,可以有效简化计算。与有限元法相比,边界元法具有降维、求解精度高、适用于无限域问题这三个优点。
[0004]结构有限元和声学边界元耦合方法已较为成熟,但鲜有人提出结构动刚度法和声学边界元耦合方法。这两种方法耦合难度较大,原因在于两种方法分别为解析法和数值方法,对网格离散要求不同,且动刚度以自由边位移作为自由度,边界元以网格节点位移作为自由度,无法做到像结构有限元
‑
声学边界元可以独立形成刚度或系数矩阵后再进行同步双向耦合。目前也没有公开文献报道结构动刚度和声学边界元直接双向耦合的数值模型。
技术实现思路
< ...
【技术保护点】
【技术特征摘要】
1.一种基于结构动刚度和声学边界元耦合的声场计算方法,其特征在于,所述方法包括:构建平板横向振动控制方程,并采用级数展开方式求解待定变量;根据所述待定变量推导平板横向振动动刚度方程;根据所述平板横向振动动刚度方程和光滑的结构流体交界面的Helmholtz积分方程,计算得到任意场点的声压值。2.根据权利要求1所述的基于结构动刚度和声学边界元耦合的声场计算方法,其特征在于,所述构建平板横向振动控制方程的方法包括:构建引入流体负载的平板横向振动控制方程,表达式为:构建引入流体负载的平板横向振动控制方程,表达式为:其中,
D=Eh3
/12(1
‑
ν2),ω是圆频率,E是杨氏模量,ν是泊松比,h是结构厚度;结构面密度m
s
=ρh,ρ是结构密度,ρ0是流体密度;w作为待定变量表示横向位移,所述横向位移是垂直于平板方向的位移;k
b
=(ρhω2/D)
1/4
,k0=ω/c,c为声波在流体中的传播速度。3.根据权利要求2所述的基于结构动刚度和声学边界元耦合的声场计算方法,其特征在于,所述采用级数展开方式求解待定变量的方法包括:采用改进Gorman叠加法,得到所述引入流体负载的平板横向振动控制方程在任意边界条件下的横向位移解为:w=w1+w2;;其中,a、b为平板的长和宽;M为cos函数展开项数;X
m
(x)、Y
m
(y)为待定卷积项,(x,y)为平板上任意点坐标。4.根据权利要求2所述的基于结构动刚度和声学边界元耦合的声场计算方法,其特征在于,根据所述待定变量推导平板横向振动动刚度方程的方法包括:根据所述横向位移确定平板四条自由边的位移列阵q
out
和外载荷列阵Q
out
;对所述位移列阵q
out
和外载荷列阵Q
out
分别在投影基上进行投影,得到所述位移列阵q
out
和外载荷列阵Q
out
分别关于待定系数列阵C
out
的关系方程为:其中,是位移列阵q
out
进行cos级数展开时的卷积系数;是外载荷列阵Q
out
进行
cos级数展开时的卷积系数;是的系数矩阵,是的系数矩阵,为已知量;联立这两个关系方程消去C
out
,得到引入流体负载的平板横向振动动刚度方程为:其中,K
out
是平板横向振动动刚度矩阵,且5.根据权利要求4所述的基于结构动刚度和声学边界元耦合的声场计算方法,其特征在于,根据所述平板横向振动动刚度方程和光滑的结构流体交界面的Helmholtz积分方程,计算得到任意场点的声压值的方法包括:将平板表面进行离散处理,划分得到多个网格单元;当外载荷已知时,根据所述引入流体负载的平板横向振动动刚度方程和网格单元的节点坐标,确定各个网格单元的平均振速对于光滑的结构流体交界面,相应的Helmholtz积分方程为:式中,S、E、V分别表示结构体的外表面、外域以及内域;Q为源点,P为场点;v
n
为结构体表面上节点处的外法向速度;p(P)为场点声压;D(Q,P)为Q点处的基本解,有:式中,r为源点到场点之间的距离,r=|P
‑
Q|;将平板表面声压数值离散后,常量声学边界元矩阵形式的方程为:其中,H、...
【专利技术属性】
技术研发人员:魏子天,殷学文,吴文伟,
申请(专利权)人:中国船舶科学研究中心,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。