基于子结构导纳法的UUV动力舱段参数化建模研究
韩飞, 王敏庆    
西北工业大学 动力与能源学院, 陕西 西安 710072
摘要: 采用子结构导纳法建立了由燃料泵、主机、隔振圈以及舱段壳体组成的动力舱段耦合结构振动传递计算模型,计算了耦合结构的机械导纳,并通过试验验证了模型的有效性。以燃料泵受到的激励载荷作为输入参数,求解了舱段壳体目标点处的振动响应,以此开展了参数影响规律研究,结果表明:在300~1 500 Hz频率范围内,采用密度及杨氏模量较小的材料作为燃料泵选材能够有效降低舱段壳体的振动响应;摆盘箱结构参数对于舱段壳体振动响应的影响比气缸体结构参数更为显著,增大摆盘箱结构的厚度能够降低舱段壳体的振动响应。
关键词: 子结构     机械导纳     振动传递    

随着现代高科技的迅速发展,无人水下航行器(unmanned underwater vehicle,UUV)的隐蔽性和攻击力不断增强。UUV的声隐身设计已成为研发和建造过程中不可或缺的重要环节。国内众多学者采用有限元方法与统计能量方法围绕UUV结构振动与辐射噪声开展了大量的研究工作。由于UUV结构非常复杂,采用有限元方法求解存在计算量巨大、计算耗时过长等问题,为了提高分析频率上限,通常只能以牺牲舱段内部结构的建模精度为代价,例如将动力舱内部结构简化为集中质量或只保留安装隔板[1, 2]。采用统计能量方法能够有效解决结构宽带高频的振动噪声问题,但受模态密度等条件限制,不适用于低频结构振动问题的研究。此外,在统计能量分析模型中,研究者通常将舱段壳体的平均振动响应作为输入参数[3, 4, 5, 6],对结构振动及辐射噪声进行预报。在工程实际中,主机及各类辅机的振动是UUV结构机械噪声的主要振动源,其振动通过隔振圈传递给舱段壳体,引起舱段壳体的弯曲振动,进而形成结构水下辐射噪声。因此,如能对动力舱段内部结构进行有效建模,从结构振动源头入手,保留结构振动传递路径的完整性,建立一套动力舱耦合结构参数化理论计算模型,对于UUV结构的减振降噪设计具有重要参考价值。

子结构导纳法的计算精度受分析频带的影响较小,经过不断发展,在工程中得到了较为广泛的应用[7, 8, 9, 10]。本文分别建立燃料泵、主机以及隔振圈的参数化等效模型,并根据各子结构间的边界连续条件,建立上述子结构组成的动力舱段耦合结构振动传递参数化模型,获得耦合结构机械导纳,以燃料泵所受实际激励力作为输入参数,求解舱段壳的振动响应,开展参数影响规律研究,为UUV结构减振降噪设计提供参考。

1 子结构建模 1.1 辅 机

UUV动力舱段内部结构主要包括主机、隔振圈以及燃料泵等辅机结构,其中辅机结构通常较为复杂,难以使用解析法对其进行建模。本文以燃料泵为研究对象,采用有限元分析软件开展模态分析,发现在100~2 000 Hz频段内仅出现一阶模态,且该阶模态振型与梁的一阶弯曲振动模态振型非常接近。进一步结合燃料泵的结构特点,本文采用一端附有集中质量的环形截面空心梁模型对燃料泵进行等效建模。等效模型的参数获取流程如图 1所示。

图 1 燃料泵等效模型参数获取流程

通过ANSYS计算得到燃料泵输入导纳与等效模型计算所得结果对比如图 2所示。其中,本文各图所示机械导纳曲线均为机械导纳的幅值。

图 2 燃料泵输入导纳对比曲线

图 2所示对比曲线可以看出,等效模型计算结果与有限元仿真结果在模态频率附近频段较为吻合,随频率变化的趋势基本一致,说明等效模型能够较为准确地表征燃料泵实际结构的振动特性。除燃料泵外,在海水泵、滑油泵等辅机结构的等效建模研究中,上述方法亦具有一定适用性。

1.2 主 机

主机是动力舱内部振动传递路径中的核心结构。正常工作中,主机内部存在多种不平衡的往复力与回转惯性力,直接作用于气缸体和摆盘箱,是引起主机结构振动的主要原因。此外,燃料泵、海水泵等辅机直接安装于主机上,存在结构间的耦合作用力。由文献[11]给出的主机模态分析结果来看,主机的振动特性主要表现为气缸体端盖、摆盘箱箱体和摆盘箱端盖的局部弯曲振动,缸体部分的振动并不明显。与燃料泵结构相比,主机振动特性较为复杂,难以采用单一模型进行等效建模。

结合主机的外部激励情况以及模态分析结果,本文采用圆板模型对气缸体端盖以及摆盘箱端盖进行等效建模,采用集中质量模型对气缸体的缸体部分进行等效建模,采用圆柱壳模型对摆盘箱箱体进行等效建模。气缸体的缸体部分振动特性单一,根据缸体部分与气缸体端盖在振动传递中属并联关系的特点,将气缸体作为整体考虑。解耦后的各子结构受力情况如图 3所示。

图 3 主机各结构受力情况

Mrlb表示燃料泵与气缸体之间的耦合弯矩;Fq表示前隔振圈与气缸体之间环向分布的径向耦合力;FhMh表示后隔振圈与摆盘箱端盖之间环向分布的径向耦合力及耦合弯矩;Fr1Fz1M1分别表示气缸体与摆盘箱箱体在耦合边界1处环向分布的径向耦合力、轴向耦合力以及耦合弯矩;Fr2Fz2M2分别表示摆盘箱箱体与摆盘箱端盖在耦合边界2处环向分布的径向耦合力、轴向耦合力以及耦合弯矩。

假设结构面内振动与其弯曲振动相互独立,根据各子结构在耦合边界处力与位移的连续条件,联立得到耦合结构的振动方程

式中,A表示气缸体的3×3阶机械导纳矩阵,B11B12B21B22分别表示摆盘箱箱体在不同位置处的3×3阶机械导纳矩阵;C表示摆盘箱端盖的3×3阶机械导纳矩阵;AeBe1Be2Ce表示各子结构所受外部激励对应的3×3阶机械导纳矩阵。

令(1)式中FqMrlbFhMh任一项为1,其余项为0,求解(1)式即可得到对应单位激励力作用下各子结构间的耦合作用力,进而获取主机任意位置处的振动响应,最终得到主机的各类机械导纳参数。

1.3 隔振圈

主机工作时产生的绝大部分振动能量经前、后隔振圈向外传递。将“主机-隔振圈-舱段壳”等效为双振子系统,则该系统的运动方程为

式中,m为舱段壳质量,x1x2分别表示舱段壳和主机的径向位移;K为隔振圈的径向等效刚度,η为橡胶块材料的损耗因子。

对(2)式进行变形,得到隔振圈径向等效刚度的表达式为

式中,T为主机到舱段壳体的径向振动传递率,可通过四端参数法求解,参见文献[12]

通过上述方法将隔振圈结构等效为主机与舱段壳之间的连接刚度,在保证结构振动路径完整性及准确反映结构振动传递特性的前提下,能够有效减少子结构数目,简化耦合结构振动方程矩阵的维数。

2 耦合结构建模

由燃料泵1、主机2、隔振圈以及舱段壳3组成的动力舱耦合结构如图 4所示。

图 4 动力舱耦合模型示意图

其中a表示燃料泵与气缸体的连接位置,b、c分别表示主机与前、后隔振圈的环向连接边界,e表示燃料泵激励力的作用位置。

假设结构的面内振动与弯曲振动相互独立,根据各子结构在耦合边界处力与位移的连续条件,可得耦合结构的边界连续方程组

式中,K′a表示燃料泵与气缸体之间的扭转刚度;KbKc分别表示前、后隔振圈的径向刚度;K′c表示后隔振圈的扭转刚度;MaFbFcMc动分标表示连接边界处的耦合作用力或耦合作用弯矩;Y2cb表示激励力作用下子结构2由位置c到位置b处的位移机械导纳,如含上标“′”则表示角位移传递导纳,如含上标“~”则表示激励弯矩对应的机械导纳。

整理为矩阵形式

式中,K=diag{K′a,Kb,Kc,K′c}表示各子结构间的连接刚度矩阵,F=MaFbFcMcT表示各子结构间的耦合力矩阵,Ω=-FY′1ea000T,E为4×4阶单位矩阵,Y4×4表示子结构机械导纳矩阵

求解(6)式得到各子结构间的耦合作用力,进而计算得到耦合结构任意位置处的振动响应。

3 试验验证

本文以燃料泵激励位置到动力舱壳体滑油口位置处的传递导纳为考察对象,测试系统如图 5所示。

图 5 测试系统示意图

在计算机测试软件中设置白噪声激励信号,由多通道分析仪发出,经功率放大器输入至激振器并作用于辅机。通过阻抗头拾取激励位置处的力信号,由加速度传感器拾取在动力舱壳体在滑油口位置处的振动响应,各组信号传递至分析仪中通过力归一化处理获取耦合结构的传递导纳。试验测试结果与本文模型仿真结果对比如图 6所示。

图 6 耦合结构传递导纳对比曲线

在300~1 500 Hz频段内,根据本文模型计算所得耦合结构输入导纳以及传递导纳与测试结果随频率变化的趋势基本一致,能够反映出试验测试曲线中的大部分峰值,说明本文模型能够有效表征耦合结构的振动传递特性。理论计算结果与试验测试结果之间存在一定误差,这是由于实际结构本身非常复杂,在等效建模研究中难以避免。

4 参数仿真

本文所建动力舱段振动传递参数化计算模型能够以实际激励条件作为输入参数,更为直观地对实际结构的振动响应进行预报,快速实现结构参数改变对于耦合结构振动特性影响的分析。

以某工况下燃料泵所受激励力作为输入参数,其频谱特性如图 7所示,其中频率分辨率为4 Hz,频率范围为300~1 500 Hz,代入本文建立的动力舱耦合结构振动传递计算模型中,同样以滑油口作为目标点,仿真计算得到其振动响应如图 8所示。由于输入力信号的采样频率较高,仿真计算得到舱段壳的振动响应曲线包含的毛刺较多,为便于对比分析,在参数影响规律研究中,仿真结果以1/3倍频程的形式给出。

图 7 激励力频谱

图 8 动力舱壳体 振动响应

首先分析当燃料泵作为主要激励源时,其材料参数对于舱段壳体振动特性的影响。按照表 1中所示的3种材料作为燃料泵选材,分别计算舱段壳滑油口处的振动响应,结果对比如图 9所示。

表 1 材料参数列表
编号材料名称密度/(kg·m-3)杨氏模量/Pa
1铝合金2 6407.10×1010
2钛合金4 5001.07×1011
37 9002.16×1011

图 9 燃料泵材料参数对舱段壳体振动响应的影响

图 9a)可以看出,选用密度以及杨氏模量相对较小的铝合金材料能够有效降低由辅机振动引起的舱段壳振动响应。进一步考察单独改变材料密度及杨氏模量对于耦合结构振动特性的影响,仿真结果如图 9b)、图 9c)所示。舱段壳振动水平随材料密度及杨氏模量的减小均呈现逐渐降低的趋势,其中改变材料密度主要影响舱段壳在800 Hz以下低频段内的振动响应,随频率升高,材料密度的影响逐渐减小;改变材料杨氏模量能够影响舱段壳在整个分析频带内的振动水平。

分别改变主机气缸体端盖厚度、摆盘箱端盖厚度以及摆盘箱箱体厚度,考察其对于舱段壳振动特性的影响,仿真结果如图 10所示。从初步的仿真结果可以看出,改变气缸体端盖厚度对于目标点振动响应的影响较小;增大摆盘箱端盖厚度能够降低舱段壳的振动响应,其影响效果随频率的升高逐渐减小;增大摆盘箱箱体厚度在整个分析频段内均能能够显著降低舱段壳的振动响应。

图 10 主机结构厚度对舱段壳体振动响应的影响
5 结 论

本文采用子结构导纳法建立了由燃料泵、主机、隔振圈以及舱段壳体组成的UUV动力舱段耦合结构振动传递参数化模型,并根据所建模型对部分参数的影响规律进行了研究,主要结论如下:

1)本文所建模型具有较高精度,能够较为准确地反映出实际耦合结构的振动特性,为进一步建立包含更多构件的动力舱振动传递参数化计算模型奠定基础。

2)在300~1 500 Hz频率范围内,采用密度及杨氏模量较小的材料作为辅机选材能够有效降低舱段壳体的振动响应;摆盘箱结构参数对于舱段壳振动响应的影响比气缸体结构参数更为显著,增大摆盘箱结构厚度能够降低舱段壳体的振动响应。

参考文献
[1] 何惠江, 李楠. 基于APDL的鱼雷壳体结构参数化建模[J]. 鱼雷技术, 2010, 18(4): 246-248 He Huijiang, Li Nan. Parametric Modeling of Torpedo Shell Structure Based on APDL[J]. Torpedo Technology, 2010, 18(4): 246-248 (in Chinese)
Cited By in Cnki (4)
[2] 曹银萍, 石秀华. 基于ANSYS的鱼雷有限元建模与模态分析[J]. 弹箭与制导学报, 2009, 29(3): 289-292 Cao Yinping, Shi Xiuhua. Finite Element Modeling and Modal Analysis of Torpedo Based on ANSYS[J]. Journal of Projectiles, Rockets, Misseles and Guidance, 2009, 29(3): 289-292 (in Chinese)
Cited By in Cnki (7)
[3] 马锐磊, 尹韶平, 曹小娟, 等. 基于FEM的鱼雷低频振动环境预示方法研究[J]. 舰船电子工程, 2014, 34(5): 136-139 Ma Ruilei, Yin Shaoping, Cao Xiaojuan, et al. Torpedo's Vibration Environment Prediction Method of Low Frequency Based on FEM[J]. Ship Electronic Engineering, 2014, 34(5):136-139 (in Chinese)
Cited By in Cnki (1)
[4] 秦晓辉, 尹韶平. 统计能量分析在鱼雷结构振动及声辐射研究中的应用[J]. 鱼雷技术, 2006, 14(1): 24-27 Qin Xiaohui, Yin Shaoping. Application of Statistical Energy Analysis to Structure Vibration and Sound Radiation of Torpedo[J]. Torpedo Technology, 2006, 14(1): 24-27 (in Chinese)
Cited By in Cnki (6)
[5] 丁少春, 朱石坚, 楼京俊. 鱼雷结构振动与声辐射的统计能量分析[J]. 浙江大学学报:工学版, 2009, 43(7): 1222-1224 Ding Shaochun, Zhu Shijian, Lou Jingjun. Statistical Energy Analysis on Structure Vibration and Sound Radiation from Torpedo[J]. Journal of Zhejiang University Engineering Scicene, 2009, 43(7): 1222-1224 (in Chinese)
Cited By in Cnki (3)
[6] 刘凯, 朱石坚, 丁少春. 基于AutoSEA的鱼雷结构辐射噪声预报方法[J]. 鱼雷技术, 2010,18(2): 91-94 Liu Kai, Zhu Shijian, Ding Shaochun. A Prediction Method for Structural Noise Radiation of Torpedo Based on AutoSEA[J]. Torpedo Technology, 2010, 18(2): 91-94 (in Chinese)
Cited By in Cnki (9)
[7] Cuschieri J M. Structural Power-Flow Analysis Using a Mobility Approach of an L-Shaped Plate[J]. Journal of the Acoustical Society of America, 1990, 87(3): 1159-1165
Click to display the text
[8] Azimi S, Hamilton J F, Soedel W. The Receptance Method Applied to the Free Vibration of Continuous Rectangular Plates[J]. Journal of Sound and Vibration, 1984, 93(1): 9-29
Click to display the text
[9] Lee Y S, Choi M H. Free Vibrations of Circular Cylindrical Shells with an Interior Plate Using the Receptance Method[J]. Journal of Sound and Vibration, 2001, 248(3): 477-497
Click to display the text
[10] Lee Y S, Choi M H, Kim J H. Free Vibrations of Laminated Composite Cylindrical Shells with an Interior Rectangular Plate[J]. Journal of Sound and Vibration, 2003, 265(4): 795-817
Click to display the text
[11] 梁跃, 何长富, 彭博. 鱼雷热动力发动机机体振动模态分析[J]. 鱼雷技术, 2005, 13(4):18-20 Liang Yue, He Changfu, Peng Bo. Structural Vibration Modal Analysis of Torpedo Thermal Power Engine[J]. Torpedo Technology, 2005, 13(4):18-20 (in Chinese)
Cited By in Cnki (10)
[12] 高爱军. 基于四端参数分析法的鱼雷动力隔振技术研究[J]. 鱼雷技术, 2007, 15(1): 29-32 Gao Aijun. Research on Vibration Isolation fo Torpedo Power Based on Four-End Parameter Analysis Method[J]. Torpedo Technology, 2007, 15(4): 29-32 (in Chinese)
Cited By in Cnki (3)
Parametric Modeling of UUV Dynamic Cabin by Substructure Receptance Method
Han Fei, Wang Minqing     
School of Power and Energy, Northwestern Ploytechnical University, Xi'an 710072, China
Abstract: A parametric model of torpedo dynamic cabin, including auxiliary machinery, engine, ring isolator and cabin shell, was established based on substructure receptance method. The mechanical receptance of the coupling structure was calculated. The model was validated by mechanical receptance experiment. With the auxiliary machinery force as the input parameter, the vibration response of cabin shell was calculated. The influence of some major parameters was analyzed. Preliminary results show that choosing material of low Young modulus and low density as the material of auxiliary machinery can reduce the vibration response of cabin shell in the frequency band of 300~1 500 Hz; the parameters of swashplate box have a larger influence on the vibration of cabin shell than that of cylinder block and increasing the thickness of swashplate box can reduce the vibration response of cabin shell.
Key words: acceleration     autonomous underwater vehicles     calculations     computer software     dynamic response     errors     experiments     finite element method     mathematical models     matrix algebra     parameterization     stiffness     vibrations(mechanical)     mechanical receptance     substructure     vibration transmission    
西北工业大学主办。
0

文章信息

韩飞, 王敏庆
Han Fei, Wang Minqing
基于子结构导纳法的UUV动力舱段参数化建模研究
Parametric Modeling of UUV Dynamic Cabin by Substructure Receptance Method
西北工业大学学报, 2016, 34(1): 106-111
Journal of Northwestern Polytechnical University, 2016, 34(1): 106-111.

文章历史

收稿日期: 2015-10-08

相关文章

工作空间