1、设备管理与维修2023 2（上）图 2基座及激振机布置0引言在研究水下结构的振动和声辐射时，不可避免会涉及流固耦合问题，对流固耦合规律的认识程度是实现对水下结构振动和声辐射预报的关键。文献 1、文献 2 分别阐述了声学和流体力学的基本原理，文献 3 研究了有限元的声辐射，文献 4、文献 5 则提出附加阻尼质量法，采用 FORTRAN 代码计算外域流体附加质量和附加阻尼矩阵，采用结构有限元分析程序对结构和内域流体作有限元分析，采用 DMAP 代码将附加质量和附加阻尼同结构质量矩阵和阻尼矩阵相叠加，实现了流固耦合计算。文献 6 依据不同域的特点，结构和流体内域采用有限元方法、流体外域采用边界元法，
2、研发了计算水下大型结构流固耦合振动和声辐射的计算程序。一些学者认为水中结构流固耦合振动及其声辐射是一个非线性问题，然而大多数研究这一问题的文章均采用线性方程进行求解。本文着重对小振幅振动时的流固耦合系统进行研究，对激振力和辐射声压的关系进行数值和试验验证，深化对水下结构耦合振动规律的认识。1基本理论结构在流体中振动时，本质上就是结构振动和流体的耦合。忽略介质对声波的粘性吸收，假定介质为均匀的理想流体，略去介质的热传导作用，考虑介质中传播的是小振幅波。在定解条件：固体壁面上的边界条件、有界条件、辐射条件下，由这些条件可以得到以速度势 控制的方程如下5。2+k2=0limrr limrrr+ik(
3、)=0nxrS0=rj|（1）其中，为拉普拉斯算子，k=c0，c0为声速，n 为物面内法线方向单位矢量，rj为 Kronecker，S0为结构湿表面；第二式为有界条件，第三式为辐射条件，第四式为固体壁面上的边界条件。式（1）是一个关于速度势 的线性定解问题，所以其描述的结构与流体耦合振动所组成的是一个线性系统。2试验验证及数值仿真2.1试验模型介绍制造了一个加肋圆柱壳，在圆柱壳内部设计了基座和 9 条等间距分布肋骨，在基座上方安装激振机（图 1图 2）。其中，圆柱壳长 1000 mm、直径 800 mm，肋骨尺寸 5 mm25 mm，肋骨水下结构激振力与辐射声压关系的验证王永富（中国人民解放军
4、海军 91970 部队，辽宁大连116041）摘要：为研究水下结构耦合振动规律，验证水下结构振动时激振力和辐射声压之间的关系，制作一个加肋圆柱壳，在指定激振频率下通过改变激振机激振力测量相应的辐射声压值，最终获得激振力和辐射声压关系曲线。对加肋圆柱壳模型进行数值仿真，采用结构有限元耦合流体边界元方法并利用有限元软件 PATRAN/NASTAN 计算水下结构振动与声辐射，在指定激振频率下计算不同激振力对应的辐射声压。试验和数值仿真计算结果均表明，在指定激振频率下激振力和辐射声压呈线性关系，小振幅振动时流固耦合振动所组成系统的是一个线性系统。关键词：流固耦合；线性系统；试验研究；有限元耦合边界元中
5、图分类号：TB532文献标识码：BDOI：10.16621/ki.issn1001-0599.2023.02.05图 1加肋圆柱壳模型13设备管理与维修2023 2（上）图 5模型有限元网格图 4水听器布置图 8155 Hz 激振频率下激振力与辐射声压关系图 3试验模型水下布置图 6激振力示意间距 100 mm；基座面板 375 mm200 mm，端盖、壳体厚度分别为 10 mm、3 mm；所用材料密度为 7.8510-9N/mm3，杨氏模量为 206 000 N/mm2，阻尼、泊松比分别为 0.06 和 0.3。2.2模型水下布置对于水下结构声辐射的测量，比较理想的环境是消声水池或开阔的湖泊
6、。本次试验是在某水库中进行，背景噪声较低，是比较好的声场环境。用压载铁对圆柱壳进行压载，并用软绳连接至双体船使圆柱壳在处于悬浮状态。模型水下布置见图 3：圆柱壳轴线距离水面 5 m（Z），水听器距离模型中心6.8 m（r）。在试验船的首部和左舷放置水听器（图 4）。2.3数值仿真（1）对圆柱壳模型进行数值仿真，采用结构有限元耦合流体边界元方法并利用有限元软件 PATRAN/NASTAN计算水下结构振动与声辐射6，在指定激振频率下计算不同激振力幅值时的辐射声压。（2）有限元网格见图 5，在基座面板施加垂向激振力 F（图 6），数值计算采用的圆柱壳模型材料属性见上文。（3）流体密度 为 1020
7、kg/m3，流体中声速 c 为 1450 m/s。（4）模型放置的位置和声场计算点与试验状态完全相同（图3、图 4）。2.4数值仿真与试验结果对比试验时，在指定激振频率下通过改变激振力的幅值得到对应状态下的辐射声压。图 7 和图 8 分别是激振频率为 128 Hz 和155 Hz 时激振力与辐射声压响应曲线。由图可知，试验测得1#和 2#水听器的激振力辐射声压关系曲线基本符合线性关系。数值计算结果同样表明激振力与辐射声压符合线性关系。对比数值计算和试验测得的辐射声压值，换算为辐射声压级后最大不超过 6 dB，导致这种误差的因素主要有圆柱壳建造焊接过程中的工艺对阻尼因子的影响、水中背景噪声的影响
8、等。上述结果表明，激振力和辐射声压是线性关系，从而验证了理论分析的结果。3结论本文通过水下加肋圆柱壳振动与噪声试验，并采用有限元耦合流体边界元法进行了数值仿真及声辐射计算，结果表明，在指定激振频率下激振力和辐射声压呈很好的线性关系，验证了小振幅振动时壳体结构与流体耦合振动所组成的是一个线性系统这一理论。参考文献1杜功焕，朱哲民，龚秀芬.声学基础（第 2 版）M.南京：南京大学出版社，2001.2张兆顺，崔桂香.流体力学 M.北京：清华大学出版社 19983L.Lagnet，Guyader J L.Sound radiation by finite circular cylindricalrin
9、g-stiffed shells J.Sound and Vibration，1990（2）：173-191.4 ZHOU Q D，ZHANG W K,Joseph P F.A new method for deter-mining acoustic added mass and damping coefficients of fluid-structure interactionC.Proceedings of the eighth internationalsymposium on practical design of ships and other floating structures，Shanghai:2001：1185-1195.5 Q.Zhou,P.F.Joseph.“A numerical method for the calculation ofdynamic response and acoustic radiation from an underwater struc-ture”J.Sound and Vibration，2005：853-873.6纪刚，张纬康，周其斗.静水压力作用的水下结构振动及声辐射 J.中国造船，2006（3）：37-44编辑吴建卿图 7128 Hz 激振频率下激振力与辐射声压关系14