CN104792447B - 一种大型船舶隔振装置动态耦合多载荷识别方法 - Google Patents
一种大型船舶隔振装置动态耦合多载荷识别方法 Download PDFInfo
- Publication number
- CN104792447B CN104792447B CN201410693028.8A CN201410693028A CN104792447B CN 104792447 B CN104792447 B CN 104792447B CN 201410693028 A CN201410693028 A CN 201410693028A CN 104792447 B CN104792447 B CN 104792447B
- Authority
- CN
- China
- Prior art keywords
- response
- point
- isolation mounting
- contact surface
- response point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Vibration Prevention Devices (AREA)
Abstract
本发明公开了一种大型船舶隔振装置动态耦合多载荷识别方法,属于船舶减振降噪技术领域,该方法包括如下步骤:一、将下层隔振器与基座连接的每个安装接触面均进行离散;二、在隔振装置未安装或拆除的情况下,在基座面板和船体结构上选取m个响应点,m>n;第k个安装接触面离散获得Nk个子面积;第i个子面积与第j个响应点之间的传递函数为Yij;所有安装接触面的分布传递函数矩阵为{H}m×n={H1 H2 ... Hn}T;其中Hk={Hk1 Hk2 ... Hkm},三、安装隔振装置,运转装置上设备,测试第j个响应点的振动速度自谱及赋相自谱,获得响应Vj,m个响应点的响应构成响应向量{V}m={V1 V2 ... Vm};四、则所有安装接触面上载荷{F}n为:
Description
技术领域
本发明属于船舶减振降噪技术领域,具体涉及一种大型船舶隔振装置动态耦合多载荷识别方法。
背景技术
机械振动噪声是船舶低速航行时的主要噪声源,机械设备进行弹性安装是保证达到目标噪声的重要手段。对于大型隔振装置,一般通过多个大型隔振器与基座相连,机械设备振动激励通过隔振装置传递至基座引起水下声辐射,隔振器与基座安装面则是振动传递的主要通道。隔振装置的动态载荷是船舶振动声辐射估算必不可少的输入参数;且可依据隔振设计要求对隔振装置动态载荷进行评估,这是工程中控制辐射噪声的重要手段。
目前船舶中机械设备动态载荷的频域识别方法包括两种:
1)基于隔振器阻抗特性的识别方法;
2)基于激励点与响应点之间传递函数求逆的识别方法。
原春晖(机械设备振动源特性测试方研究.武汉:华中科技大学博士论文,2006)对两种方法进行了试验验证。郭荣等(一种基于正则化和奇异值分解的结构载荷识别方法。发明专利,申请号201310073284.2)针对点接触载荷,采用Tikhonov正则化法对载荷进行识别。庞福振等(不平衡激励力间接测量方法。发明专利,申请号201210185489.5)针对点接触安装设备提供不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法。以上方法为将复杂问题简单化,通常认为设备安装机脚、隔振器及基座三者接触部位都是“点”而不是有一定面积的接触“面”,所考虑的力和响应的作用面积均被简化为一个点。这种点接触假设在接触面相对比较小的情况下成立,适应于辅机设备。
对于船舶大型浮筏隔振装置,其隔振器与基座接触面大(最大尺寸接近0.5米),安装面上载荷以及响应速度为空间分布形式,不再满足点接触假设,且多个安装面之间还存在动态耦合,急需新的动态载荷识别方法。
发明内容
有鉴于此,本发明提供了一种大型船舶隔振装置动态耦合多载荷识别方法,能够针对隔振器与基座接触面大的船舶大型浮筏隔振装置进行多载荷识别,并充分考虑到各安装接触面之间的动态耦合
一种大型船舶隔振装置动态耦合多载荷识别方法,大型船舶隔振装置中包括下层隔振器和上层隔振器,其中上层隔振器与设备相连,下层隔振器与基座相连,下层隔振器和基座之间具有n个安装接触面,其中最大尺寸与基座振动波长相近,该方法包括如下步骤:
第一步、对于每个安装接触面均进行离散,每个安装接触面离散后获得若干子面积;每个子面积测点满足点接触。
第二步、在不安装隔振装置的情况下,在基座面板和船体结构上选取m个响应点,m>n;对于第k个安装接触面,其离散获得Nk个子面积。
获取第i个子面积与第j个响应点之间的传递函数Yij,其中k=1,2,…,n;i=1,2,…,Nk;j=1,2,…,m。
则第k个安装接触面与第j个响应点之间的形成的传递函数为第k个安装接触面与m个响应点之间的形成的传递函数矩阵为Hk={Hk1 Hk2 ... Hkm};所有安装接触面形成的分布传递函数矩阵为{H}m×n={H1 H2 ... Hn}T。
第三步、在安装了隔振装置的情况下,运转船舶上设备,测试第j个响应点的振动速度自谱及赋相自谱,获得响应Vj,共m个响应点的响应构成响应向量{V}m={V 1V2 ...Vm}。
第四步、则隔振器和基座之间所有安装接触面上载荷{F}n为:为对{H}m×n采用奇异值分解求伪逆。
进一步地,子面积的尺寸△L为λ/5~λ/3,λ为最高分析频率时基座面板的弯曲波波长。
进一步地,响应点个数m比安装接触面个数n多4个以上。
进一步地,第二步中获取第i个子面积与第j个响应点之间的传递函数的方法为:在第i个子面积上选取一个测点A,采用力锤或激振机在测点A处进行激励,测试获得测点A与响应点之间传递函数,或者根据互易原理,在响应点处激励,测试测点A处的响应,获得测点A与响应点之间传递函数。
进一步地,在基座面板上选取响应点的位置包括加强腹板位置、交叉肋板以及无加强腹板位置;在船体结构上选取的响应点包括加强肋和非加强肋的位置。
有益效果:
1、本发明主要针对隔振器与基座接触面大的船舶大型浮筏隔振装置进行多载荷识别,通过接触面离散、测点选取以及建立分布式的传递函数矩阵的方式,由于各个接触面振动之间存在耦合,最终获得额传递函数矩阵满矩阵,该矩阵与响应向量相乘则可以进一步获得每个安装接触面上的总载荷;该方法利用了各个接触面之间的动态耦合的特性,以矩阵的方式解决了隔振装置对基座接触面积较大的情况下的动态载荷识别的问题,从而为进一步评估隔振装置隔振效果以及简化声学设计提供了便利。
2、本发明中在进行安装接触面的离散时,依据最高分析频率时基座面板的弯曲波波长给出了离散后的子面积尺寸,在该尺寸下能够保证每个子面积测点满足点接触。
附图说明
图1为本发明所针对的中大型船舶隔振装置图示,图中:设备1,筏架2,下层隔振器3,基座4和船体5。
图2为本本发明实施例2中的4个安装面的分布传递函数测试图示。图中:每个安装面6离散成2×2共4个子面积7,每个子面积7中心布置一个测点9(■),基座面板8上响应点10(●),船体5结构上响应点11(▲)。
图3为本发明实施例2中的分布传递函数识别方法图示。
图4为本发明实施例2中的设备运转下响应点布置图示。图中:基座面板8上响应点10(●),船体结构上响应点11(▲),相位参考点12(■)。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
实施例1、本实施例中提出的一种大型船舶隔振装置动态耦合多载荷识别方法,该方法应用的前提是:其所针对的大型船舶隔振装置包括下层隔振器和上层隔振器,其中上层隔振器与设备相连,如图1所示,其中包括设备1,筏架2,下层隔振器3,基座4和船体5;下层隔振器和基座之间具有n个安装接触面,其中最大尺寸与基座振动波长相近,形状一般为矩形或圆形,此时隔振器作用在基座上的动态载荷为分布载荷,且安装面内和周围响应速度也是空间分布。
该方法包括如下步骤:
第一步、在上述前提之下,隔振器与基座的接触部位为“面”,载荷作用面与响应点之间关系用分布传递函数描叙,不再是点对点的传递函数。借鉴“点”接触下载荷识别方法,提出基于载荷作用面与响应点之间分布传递函数矩阵求逆的识别方法。
由于接触面不再是点,需要将接触面进行离散才能获取分布传递函数,因此本步骤中对于每个安装接触面均进行离散,每个安装接触面离散后获得若干子面积;每个子面积测点满足点接触。
本实施例中每个离散后子面积的尺寸△L为λ/5~λ/3,λ为最高分析频率时基座面板的弯曲波波长,对于长为a,宽为b的接触面,离散后子面积数量分别为(a/ΔL)×(b/ΔL)。此时每个子面积测点是满足点接触的。
第二步、在不安装隔振装置的情况下,在基座面板和船体结构上选取m个响应点,m>n;对于第k个安装接触面,其离散获得Nk个子面积;
获取第i个子面积与第j个响应点之间的传递函数Yij,其中k=1,2,…,n;i=1,2,…,Nk;j=1,2,…,m。获取Yij的方法为:在第i个子面积上选取一个测点A,采用力锤或激振机在测点A处进行激励,测试获得测点A与响应点之间传递函数,或者根据互易原理,在响应点处激励,测试测点A处的响应,获得测点A与响应点之间传递函数。
船舶隔振装置下层隔振器相对于基座较软,此时隔振器作用在基座上的激励力可以近似为均布力激励。将安装接触面离散后,根据点传递函数合成分布传递函数,得到第k个安装接触面与第j个响应点之间的形成的传递函数为
第k个安装接触面与m个响应点之间的形成的传递函数矩阵为Hk={Hk1 Hk2 ...Hkm};所有安装接触面形成的分布传递函数矩阵为{H}m×n={H1 H2 ... Hn}T;
第三步、在安装了隔振装置的情况下,运转船舶上设备,测试第j个响应点的振动速度自谱及赋相自谱,获得响应Vj,共m个响应点的响应构成响应向量{V}m={V1 V2 ...Vm};
第四步、基座接触面上载荷{F}与测点响应{V}之间的关系为:{V}m={H}m×n{F}n。
其中载荷{F}为各个接触面上总载荷构成的向量。{H}m×n为分布传递函数矩阵,与点接触的载荷识别方法相比,此处的{H}m×n为分布传递函数矩阵而不是点传递函数矩阵,由于各个接触面振动之间存在耦合,导致{H}m×n为满矩阵。m为响应点的个数,n为载荷的个数(一个接触面只有一个总载荷),根据方程求解的要求,一般要求m>n。本实施例中m比n多4~6个以上。
对分布传递函数矩阵{H}m×n采用奇异值分解求伪逆,得到接触面上的载荷:
其中“+”表示伪逆。
则隔振器和基座之间所有安装接触面上载荷{F}n为:为对{H}m×n采用奇异值分解求伪逆。
实施例2、以下结合附图进一步说明本智力成果的实施例。
测点的布置如图2所示,图中:每个安装面6离散成2×2共4个子面积7,每个子面积7中心布置一个测点9(■),基座面板8上响应点10(●),船体5结构上响应点11(▲)。基座面板8上响应点10位置应靠近安装面6,基座面板上典型测点位置包括加强腹板位置、交叉肋板以及无加强腹板位置。响应点11位置应靠近基座在艇体上生根的附近,艇体上响应点包括加强肋和非加强肋的位置,且数量要比安装面6多。
一般情况下响应点个数m比载荷个数n多4~6个以上。
为了识别耦合多载荷,需对每个响应点与安装面之间的分布传递函数进行测试,共有m×n个分布传递函数。
测试步骤如下:
拆除隔振装置或在隔振装置安装之前对隔振器与基座连接的接触面6进行离散;
将接触面6离散后,在隔振器安装的基座面板8上和船体5上选择响应参考点10和11;
在没有隔振装置的情况下,如图3所示,采用力锤或激振机测试安装面6内子面积7的测点9与响应点10和11之间点传递函数。还可根据互易原理,在响应点10和11处激励,测试子面积7内点测点9的响应,依次测试所有子面积形成分布传递函数矩阵{H}m×n;
安装隔振装置后运转设备,如图4所示,测试各个响应点10和11振动速度自谱及赋相自谱,构成响应向量{V}m,测试中为排除振动干扰应停止其他无关设备;
利用所测速度数据和计算参数带入动态载荷识别公式,计算得到隔振器对基座上每个接触面的动态载荷。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种大型船舶隔振装置动态耦合多载荷识别方法,所述大型船舶隔振装置包括下层隔振器和上层隔振器,其中上层隔振器与机械设备相连,下层隔振器与基座相连,下层隔振器和基座之间具有n个安装接触面,其中最大尺寸大于或者等于基座振动波长的1/3,该方法的特征在于,包括如下步骤:
第一步、对于下层隔振器与基座连接的每个安装接触面均进行离散,每个安装接触面离散后获得若干子面积;每个子面积测点满足点接触;
第二步、在不安装隔振装置的情况下,在基座面板和船体结构上选取m个响应点,m>n;对于第k个安装接触面,其离散获得Nk个子面积;
获取第i个子面积与第j个响应点之间的传递函数Yij,其中k=1,2,…,n;i=1,2,…,Nk;j=1,2,…,m;
则第k个安装接触面与第j个响应点之间的形成的传递函数为第k个安装接触面与m个响应点之间的形成的传递函数矩阵为Hk={Hk1 Hk2 ... Hkm};所有安装接触面形成的分布传递函数矩阵为{H}m×n={H1 H2 ... Hn}T;
第三步、在安装了隔振装置的情况下,运转船舶上所述机械设备,测试第j个响应点的振动速度自谱及赋相自谱,获得响应Vj,共m个响应点的响应构成响应向量{V}m={V1 V2... Vm};
第四步、则下层隔振器和基座之间所有安装接触面上载荷{F}n为: 为对{H}m×n采用奇异值分解求伪逆。
2.如权利要求1所述的一种大型船舶隔振装置动态耦合多载荷识别方法,其特征在于,所述子面积的尺寸△L为λ5~λ3,λ为最高分析频率时基座面板的弯曲波波长。
3.如权利要求1所述的一种大型船舶隔振装置动态耦合多载荷识别方法,其特征在于,所述响应点个数m比安装接触面个数n多4个以上。
4.如权利要求1所述的一种大型船舶隔振装置动态耦合多载荷识别方法,其特征在于,所述第二步中获取第i个子面积与第j个响应点之间的传递函数的方法为:在第i个子面积上选取一个测点A,采用力锤或激振机在测点A处进行激励,测试获得测点A与响应点之间传递函数,或者根据互易原理,在响应点处激励,测试测点A处的响应,获得测点A与响应点之间传递函数。
5.如权利要求1所述的一种大型船舶隔振装置动态耦合多载荷识别方法,其特征在于,在基座面板上选取响应点的位置包括加强腹板位置、交叉肋板以及无加强腹板位置;在船体结构上选取的响应点包括加强肋和非加强肋的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693028.8A CN104792447B (zh) | 2014-11-26 | 2014-11-26 | 一种大型船舶隔振装置动态耦合多载荷识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693028.8A CN104792447B (zh) | 2014-11-26 | 2014-11-26 | 一种大型船舶隔振装置动态耦合多载荷识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104792447A CN104792447A (zh) | 2015-07-22 |
CN104792447B true CN104792447B (zh) | 2017-10-13 |
Family
ID=53557439
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410693028.8A Active CN104792447B (zh) | 2014-11-26 | 2014-11-26 | 一种大型船舶隔振装置动态耦合多载荷识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104792447B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108553A (zh) * | 2017-12-18 | 2018-06-01 | 中国舰船研究设计中心 | 一种舰船浮筏隔振系统摇摆稳定性计算方法 |
CN114199436B (zh) * | 2021-12-10 | 2023-07-21 | 哈尔滨工程大学 | 一种基于传感器数据的船舶监测结构整体应力状态快速获取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680270A (zh) * | 2012-06-07 | 2012-09-19 | 哈尔滨工程大学 | 不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法 |
CN102692331A (zh) * | 2012-06-07 | 2012-09-26 | 哈尔滨工程大学 | 垂向不平衡激励力作用下设备对船体结构激励载荷的间接测试方法 |
CN102779200A (zh) * | 2011-12-14 | 2012-11-14 | 西安交通大学 | 一种包含微观表面形貌的结合面接触性能分析方法 |
CN103196591A (zh) * | 2013-03-07 | 2013-07-10 | 同济大学 | 一种基于正则化和奇异值分解的结构载荷识别方法 |
CN203627643U (zh) * | 2013-12-26 | 2014-06-04 | 中国舰船研究设计中心 | 具有周期结构效应的机械隔振平台 |
CN103902834A (zh) * | 2014-04-14 | 2014-07-02 | 重庆大学 | 一种基于岭估计和l曲线法的结构损伤识别方法 |
-
2014
- 2014-11-26 CN CN201410693028.8A patent/CN104792447B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779200A (zh) * | 2011-12-14 | 2012-11-14 | 西安交通大学 | 一种包含微观表面形貌的结合面接触性能分析方法 |
CN102680270A (zh) * | 2012-06-07 | 2012-09-19 | 哈尔滨工程大学 | 不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法 |
CN102692331A (zh) * | 2012-06-07 | 2012-09-26 | 哈尔滨工程大学 | 垂向不平衡激励力作用下设备对船体结构激励载荷的间接测试方法 |
CN103196591A (zh) * | 2013-03-07 | 2013-07-10 | 同济大学 | 一种基于正则化和奇异值分解的结构载荷识别方法 |
CN203627643U (zh) * | 2013-12-26 | 2014-06-04 | 中国舰船研究设计中心 | 具有周期结构效应的机械隔振平台 |
CN103902834A (zh) * | 2014-04-14 | 2014-07-02 | 重庆大学 | 一种基于岭估计和l曲线法的结构损伤识别方法 |
Non-Patent Citations (2)
Title |
---|
"船舶机电设备基座面导纳及其应用";魏强等;《2005年船舶结构力学学术会议论文集》;20050930;第490-495页 * |
Identification of external structural loads from measured harmonic responses;Karlsson SES;《Jpurnal or Sound and Vibration》;19961231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104792447A (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN201408112Y (zh) | 一种风力发电机组风轮叶片加载试验装置 | |
Cruz et al. | Estimating the loads and energy yield of arrays of wave energy converters under realistic seas | |
CN102539101B (zh) | 力限控制振动试验系统及试验方法 | |
Day et al. | The 16-day wave in the Arctic and Antarctic mesosphere and lower thermosphere | |
CN104792447B (zh) | 一种大型船舶隔振装置动态耦合多载荷识别方法 | |
US9949025B2 (en) | Systems and methods for detecting transient acoustic signals | |
CN104792450B (zh) | 一种设备对安装基座激励力的间接测量方法 | |
Thys et al. | Real-time hybrid model testing of a semi-submersible 10MW floating wind turbine and advances in the test method | |
CN108007681A (zh) | 一种利用麦克风阵列进行机械故障检测的方法 | |
Oikonomou et al. | On the dynamics of an array of spar-buoy oscillating water column devices with inter-body mooring connections | |
Zhang et al. | Sensitivity analysis of an auto-correlation-function-based damage index and its application in structural damage detection | |
Zhao et al. | A diffusion approximation for ocean wave scatterings by randomly distributed ice floes | |
CN107782547A (zh) | 一种传动带的模态分析方法及装置 | |
CN106768767A (zh) | 一种基于频响函数的轴承座特性参数的测量系统及测量方法 | |
Si et al. | Random vibration of an elastic half-space subjected to a moving stochastic load | |
Nagel et al. | Ground motions induced by wind turbines | |
Maxit et al. | Numerical analysis of vibroacoustic beamforming gains for acoustic source detection inside a pipe conveying turbulent flow | |
Gong et al. | Experimental and numerical investigation on an optimization method of heaving buoy wave energy converter arrays based on a given target wave spectrum | |
Dunham et al. | Earthquake ruptures on rough faults | |
Zhu et al. | A modal test method using sound pressure transducers based on vibro-acoustic reciprocity | |
CN104330150B (zh) | 一种多面安装设备的多自由度振动激励力间接测试方法 | |
US8588030B2 (en) | Reducing noise in seismic surveying environment | |
Sutherland | Dynamic excitation of internal gravity waves in the equatorial oceans | |
Loh et al. | Experimental verification of building control using active bracing system | |
Meyer et al. | On the mooring methodology of heaving point absorber arrays |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |