CN114199438B - 一种基于船舶结构监测数据分析的载荷识别方法 - Google Patents
一种基于船舶结构监测数据分析的载荷识别方法 Download PDFInfo
- Publication number
- CN114199438B CN114199438B CN202111509674.0A CN202111509674A CN114199438B CN 114199438 B CN114199438 B CN 114199438B CN 202111509674 A CN202111509674 A CN 202111509674A CN 114199438 B CN114199438 B CN 114199438B
- Authority
- CN
- China
- Prior art keywords
- monitoring
- ship
- monitoring structure
- external
- stress
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0028—Force sensors associated with force applying means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明属于船体结构应力监测技术领域,具体涉及一种基于船舶结构监测数据分析的载荷识别方法。本发明通过结合船舶监测结构的结构分析结果,采用船舶监测结构的监测点优化布置方法,确定出船舶监测结构各监测点的具体位置;然后,通过结合由传感器获取的应力数据,采用对应的载荷识别方法,可依次判断给出船舶监测结构实际遭受的外部载荷作用情况,实现对船舶监测结构实际遭受的外部载荷形式的准确、快速识别,进而为实现船舶监测结构整体应力状态的快速获取提供了有力支撑。
Description
技术领域
本发明属于船体结构应力监测技术领域,具体涉及一种基于船舶结构监测数据分析的载荷识别方法。
背景技术
在船体结构应力监测中,通过采用在监测点位置布置传感器的方式,实现了对船舶结构进行健康监测的目的,但这仅能获得监测点位置的应力状态;而若想得到船舶监测结构更多位置的应力状态,就需要相应地增加传感器的布置数量。而布置的传感器数量过多,会带来许多问题,如:不利于经济性的要求;数据量成倍地增加,会给监测系统的数据处理带来很大困难等。因此,这就造成了目前的结构监测方法很难获得船舶监测结构整体应力状态的问题。事实上,由于船上人员作业、传感器的安装工艺和走线、监测点位置的结构复杂等因素的制约,会造成传感器的实际安装位置不在选定的监测点位置的情况;同时,对于在实际海洋环境中营运的船舶,我们也无法确定事先选定的监测点位置是否一直是结构的最危险位置,而这将会给船舶结构的安全性带来很大的隐患。事实上,若能通过由布置的有限数量传感器获得的监测数据,进而得到船舶监测结构的整体应力状态,就可以很好地解决上述问题。特别地,结构实际遭受外部载荷的准确、快速识别,是船舶监测结构整体应力状态快速获取的关键,应该给予足够的重视和研究。
发明内容
本发明的目的在于提供一种基于船舶结构监测数据分析的载荷识别方法。
一种基于船舶结构监测数据分析的载荷识别方法,包括以下步骤:
步骤1:根据船舶监测结构可能遭受m种形式外部载荷Fsj作用的情况,确定在船舶监测结构上布置监测点的数量m;j=1,2,…,m,m>2;
步骤2:对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;然后,分别对该离散模型施加上述各种形式外部载荷的单位载荷Fsj0,完成船舶监测结构的结构分析;结合结构分析结果,初步选定m个监测点位置ri;i=1,2,…,m;
步骤3:获取与每种形式外部载荷Fsj对应的各监测点位置的单位应力数据σ′ri(Fsj0),分别计算与各种形式外部载荷Fsj对应的监测点间的应力关系系数
步骤4:对于任意取定的p、q值,若存在则重新选择m个监测点位置ri,返回步骤3;若不存在/>则确定选定的m个监测点位置ri即为船舶监测结构的最终的监测点位置;
步骤5:获取布置在各监测点位置ri的传感器在遭受实际外部载荷的作用下的应力数据σ′ri(Fs),结合与每种形式外部载荷Fsj对应的监测点位置的单位应力数据σ′ri(Fsj0),计算在监测点ri位置与每种形式外部载荷Fsj对应的应力放大系数ksj-ri;
若存在ksj-r1=ksj-r2=…=ksj-rm,且对于任意n≠j,均不存在ksn-r1=ksn-r2=…=ksn-rm,则判定船舶监测结构此时仅遭受了外部载荷Fsj的作用;
若对于任意j,均不满足ksj-r1=ksj-r2=…=ksj-rm,则判定船舶监测结构此时遭受了多种形式外部载荷的组合作用;
步骤6:若船舶监测结构此时仅遭受了外部载荷Fsj的作用,则步骤5中计算得到的ksj-ri即为与该种形式外部载荷Fsj对应的应力放大系数ksj,则船舶监测结构实际遭受的外部载荷Fsj的数值大小为:Fsj=ksj·Fsj0;此时,由布置在监测点位置的传感器获得的应力数据σ′ri(Fs),即为与船舶监测结构实际遭受外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj);
步骤7:判断船舶监测结构同时遭受的外部载荷的数量x;
步骤7.1:初始化x=2;
步骤7.2:任取x种形式外部载荷判断船舶监测结构是否同时遭受了x种形式外部载荷/>的组合作用;
先由(1)、…、(x)联立求得与外部载荷对应的应力放大系数/>再分别由(x+1)、…、(m)验证所求的/>是否满足方程的等式关系;若方程(x+1)、…、(m)均满足等式关系,则判定船舶监测结构同时遭受了x种形式外部载荷/>的组合作用;
步骤7.3:若步骤7.2中α1、…、αx遍历其全部可能的取值后,仍无法判定船舶监测结构同时遭受了x种形式外部载荷的组合作用,且x<m-1,则令x=x+1,返回步骤7.2;若x=m-1,则判定船舶监测结构同时遭受了所有m种形式外部载荷的组合作用,令x=m,由步骤7.2中(1)、…、(m)联立求得与外部载荷/>对应的应力放大系数/>
步骤8:对于船舶监测结构同时遭受的x种形式外部载荷由步骤7.2中求得的/>即为与外部载荷/>对应的应力放大系数/>船舶监测结构实际遭受的x种形式外部载荷/>的数值大小为:
与船舶监测结构实际遭受的x种形式外部载荷对应的监测点ri位置的应力数据/>的计算方法为:
本发明的有益效果在于:
本发明通过结合船舶监测结构的结构分析结果,采用船舶监测结构的监测点优化布置方法,确定出船舶监测结构各监测点的具体位置;然后,通过结合由传感器获取的应力数据,采用对应的载荷识别方法,可依次判断给出船舶监测结构实际遭受的外部载荷作用情况,实现对船舶监测结构实际遭受的外部载荷形式的准确、快速识别。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明做进一步描述。
在现有的船体结构应力监测中,仅能获得结构在传感器布置位置的应力状态,而无法知晓船舶监测结构的整体应力状态。事实上,由于部分现实因素的制约,会造成传感器的实际安装位置不在选定的监测点位置的情况;同时,对于在实际海洋环境中营运的船舶,我们也无法确定事先选定的监测点位置是否一直是船舶监测结构的最危险位置。若出现监测点位置并非结构最危险位置的情况,就会给船舶的安全性带来很大的隐患。而若能通过由布置在监测点位置的传感器获得的应力数据,得到船舶监测结构的整体应力状态,就可以很好地解决上述问题。特别地,结构实际遭受外部载荷的准确、快速识别,是船舶监测结构整体应力状态快速获取的关键,应该给予足够的重视和研究。
基于此,本发明提出了一种基于船舶结构监测数据分析的载荷识别方法,具体流程见图1。该方法通过由布置在监测点位置的传感器获得的应力数据,结合具体的载荷识别方法,能够实现对船舶监测结构实际遭受外部载荷的准确、快速识别,进而为实现船舶监测结构整体应力状态的快速获取提供了有力支撑。
因此,针对船舶监测结构实际遭受外部载荷的识别研究,本发明采用了这样的思路:首先,在分析了结构可能遭受m(m>2)种形式外部载荷的作用后,确定在船舶监测结构上布置相应数量的监测点。然后,对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;分别对该离散模型施加上述各种形式外部载荷的单位载荷,完成船舶监测结构的结构分析;并由船舶监测结构的监测点优化布置方法,最终确定各监测点的位置。最后,结合由布置在监测点位置的传感器获得的应力数据,采用一种形式外部载荷单独作用对应的载荷识别方法、多种形式外部载荷组合作用对应的载荷识别方法、全部形式外部载荷组合作用对应的载荷识别方法等,具体判断船舶监测结构实际遭受的外部载荷作用情况。
首先,根据船舶监测结构可能遭受m(m>2)种形式外部载荷作用的情况,确定在船舶监测结构上布置相应数量的监测点。其次,对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;分别对该离散模型施加上述各种形式外部载荷的单位载荷,完成船舶监测结构的结构分析,并由船舶监测结构的监测点优化布置方法,最终确定各监测点的位置。然后,结合由布置在监测点位置的传感器获得的应力数据,由一种形式外部载荷单独作用对应的载荷识别方法,判断船舶监测结构是否实际仅遭受了一种形式外部载荷的单独作用。再后,若判定了船舶监测结构并非仅遭受一种形式外部载荷的单独作用,则由多种形式外部载荷组合作用对应的载荷识别方法,依次判断船舶监测结构是否同时遭受了两种形式、…、m-1种形式外部载荷的组合作用。最后,若直至最终判定了船舶监测结构也并非遭受了m-1种形式外部载荷的组合作用,则可确定船舶监测结构同时遭受全部形式外部载荷的组合作用。
本发明的原理如下:
1)船舶监测结构的监测点优化布置方法
由于船舶监测结构可能遭受m(m>2)种形式外部载荷Fsj(j=1,2,…,m)的作用,进而需要布置相应数量的监测点。首先,对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;分别对该离散模型施加上述各种形式外部载荷的单位载荷Fsj0,结合结构分析结果,初步选定m个监测点ri(i=1,2,…,m)位置,并获得与每种形式外部载荷对应的各监测点位置的单位应力数据σ′ri(Fsj0)。
然后,由公式(1)分别计算得到与上述各种形式外部载荷对应的监测点间的应力关系系数
最后,结合已获得的监测点间的应力关系系数对于任意取定的p、q值,判断是否存在这样的等式关系:/>若没有这样的等式关系存在,则确定初步选定的监测点位置为最终的监测点位置;若有这样的等式关系存在,则应考虑更换监测点的位置,并重复上述过程,直至获得满足要求的各监测点位置ri,以避免出现无法判断船舶监测结构实际遭受外部载荷的情况。
2)一种形式外部载荷单独作用对应的载荷识别方法
根据由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),结合与每种形式外部载荷对应的单位应力数据σ′ri(Fsj0),首先获得在监测点ri位置,与每种形式外部载荷对应的应力放大系数ksj-ri:
然后,结合由公式(2)获得的与各种形式外部载荷对应的应力放大系数ksj-ri,判断是否有如下关系的存在:
ks1-r1=ks1-r2=…=ks1-rm (3-1)
ks2-r1=ks2-r2=…=ks2-rm (3-2)
ks3-r1=ks3-r2=…=ks3-rm (3-3)
……
ksm-r1=ksm-r2=…=ksm-rm (3-m)
若由公式(2)获得的各应力放大系数ksj-ri满足(3-1),而不满足(3-2)、(3-3)、…、(3-m),则说明船舶监测结构此时仅遭受了外部载荷Fs1的作用,其余同理。在判断出船舶监测结构仅遭受了一种形式外部载荷Fsj的作用后,则由公式(2)计算得到的应力放大系数ksj-ri即为与该种形式外部载荷Fsj对应的应力放大系数ksj;相应地,由公式(4)可计算得到船舶监测结构实际遭受的外部载荷Fsj的数值大小。此时,由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),即为与船舶监测结构实际遭受的外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj)。
Fsj=ksj·Fsj0 (4)
若上述各应力放大系数ksj-ri均不满足(1)、(2)、(3)、…、(m)中的关系,则可判断船舶监测结构此时并非遭受了一种形式外部载荷的单独作用,而是遭受了多种形式外部载荷的组合作用。
3)多种形式外部载荷组合作用对应的载荷识别方法
首先,判断船舶监测结构是否同时遭受了两种形式外部载荷的组合作用。此时,根据由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),结合与每种形式外部载荷对应的监测点ri位置的单位应力数据,由方程(5)判断船舶监测结构是否实际遭受了两种形式外部载荷Fsα(j=α)、Fsβ(j=β)(β>α)的组合作用。
具体地,先由5-(a)、5-(b)联立求得与外部载荷Fsα、Fsβ对应的应力放大系数k′sα、k′sβ,再分别由5-(c)、…、5-(m)验证所求的k′sα、k′sβ是否满足方程的等式关系;若对应的方程均满足等式关系,则说明船舶监测结构同时遭受了两种形式外部载荷Fsα、Fsβ的组合作用,此时,上述所求得的k′sα、k′sβ即为与外部载荷Fsα、Fsβ对应的应力放大系数ksj(j=α,β)。相应地,由公式(4)可对应地计算得到船舶监测结构实际遭受的外部载荷Fsα、Fsβ的数值大小;由公式(6)可依次计算得到由传感器获得的、与外部载荷Fsj(j=α,β)对应的监测点ri位置的应力数据σ′ri(Fsj)。
σ′ri(Fsj)=ksj·σ′ri(Fsj0) (6)
若上述所求的k′sα、k′sβ使5-(c)、…、5-(m)中存在不满足等式关系的方程式,则考虑更新α、β的数值,直至α、β遍历其全部可能的取值;若此时,仍无法获得符合要求的外部载荷Fsα、Fsβ,则可判断船舶监测结构此时并非同时遭受了两种形式外部载荷的组合作用。
然后,若有m>3,则继续判断船舶监测结构是否同时遭受了三种形式外部载荷的组合作用。此时,根据由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),结合与每种形式外部载荷对应的监测点ri位置的单位应力数据,由方程(7)判断船舶监测结构是否实际遭受了三种形式外部载荷Fsα(j=α)、Fsβ(j=β)、Fsγ(j=γ)(γ>β>α)的组合作用。
具体地,先由7-(a)、7-(b)、7-(c)联立求得与外部载荷Fsα、Fsβ、Fsγ对应的应力放大系数k′sα、k′sβ、k′sγ,再分别由剩余方程验证所求的k′sα、k′sβ、k′sγ是否满足方程的等式关系;若对应的方程均满足等式关系,则说明船舶监测结构同时遭受了三种形式外部载荷Fsα、Fsβ、Fsγ的组合作用,此时,上述所求得的k′sα、k′sβ、k′sγ,即为与外部载荷Fsα、Fsβ、Fsγ对应的应力放大系数ksj(j=α,β,γ)。同样地,由公式(4)可对应地计算得到船舶监测结构实际遭受的外部载荷Fsα、Fsβ、Fsγ的数值大小;由公式(6)可依次计算得到由传感器获得的、与外部载荷Fsj(j=α,β,γ)对应的监测点ri位置的应力数据σ′ri(Fsj)。
若上述所求的k′sα、k′sβ、k′sγ使剩余方程中存在不满足等式关系的方程式,则考虑更新α、β、γ的数值,直至α、β、γ遍历其全部可能的取值;若此时仍无法获得符合要求的外部载荷Fsα、Fsβ、Fsγ,则可判断船舶监测结构此时并非遭受了三种形式外部载荷的组合作用。
顺序地,采用同样的方法,依次判断船舶监测结构是否同时遭受了其余多种形式外部载荷的组合作用,若最终判断船舶监测结构此时也并非遭受了m-1种形式外部载荷的组合作用,则说明船舶监测结构同时遭受了全部形式外部载荷的组合作用。
4)全部形式外部载荷组合作用对应的载荷识别方法
在判定了船舶监测结构同时遭受了全部形式外部载荷的组合作用后,根据由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),结合与每种形式外部载荷对应的监测点位置的单位应力数据,由方程(8)可求得与各种形式外部载荷Fsj对应的应力放大系数ksj:
在分别获得了与船舶监测结构实际遭受的各种形式外部载荷Fsj对应的应力放大系数ksj后,同样可由公式(4)对应地计算得到船舶监测结构实际遭受的外部载荷Fsj的数值大小;由公式(6)可依次计算得到由传感器获得的、与外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj)。
本发明的具体方法如下:
1)根据船舶监测结构可能遭受多种形式外部载荷Fsj(j=1,2,…,m,m>2)作用的情况,确定在船舶监测结构上布置m个监测点。
2)对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;然后,分别对该离散模型施加上述各种形式外部载荷的单位载荷Fsj0,完成船舶监测结构的结构分析;结合结构分析结果,初步选定m个监测点位置ri(i=1,2,…,m,m>2),并获得与每种形式外部载荷对应的各监测点位置的单位应力数据σ′ri(Fsj0)。
3)根据已获得的上述各单位应力数据σ′ri(Fsj0),由公式(1)分别计算得到与上述各种形式外部载荷对应的监测点间的应力关系系数对于任意取定的p、q值,判断是否存在有/>这样的等式关系。
4)结合3)中的判断结果,若存在有这样的等式关系,则需要重新确定监测点位置,并重复上述3)、4)的过程,直至获得满足要求的各监测点ri位置;若不存在有/>这样的等式关系,则可确定初步选定的监测点位置即为船舶监测结构的最终的监测点位置。
5)根据由布置在监测点ri位置的传感器获得的应力数据σ′ri(Fs),结合与每种形式外部载荷对应的监测点位置的单位应力数据σ′ri(Fsj0),由公式(2)计算得到在监测点ri位置,与每种形式外部载荷对应的应力放大系数ksj-ri;然后,结合一种形式外部载荷单独作用对应的载荷识别方法,判断船舶监测结构是否仅遭受了一种形式外部载荷的作用。
6)根据5)中的判断结果,若判定了船舶监测结构实际仅遭受了一种形式外部载荷的单独作用,则由公式(2)计算得到的应力放大系数ksj-ri,即为与该种形式外部载荷Fsj对应的应力放大系数ksj;进而,由公式(4)可计算得到船舶监测结构实际遭受的外部载荷Fsj的数值大小;此时,由布置在监测点位置的传感器获得的应力数据σ′ri(Fs),即为与船舶监测结构实际遭受外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj)。
7)在判定船舶监测结构实际遭受了多种形式外部载荷的组合作用后,首先,由方程(5)判断船舶监测结构是否同时遭受了两种形式外部载荷Fsα(j=α)、Fsβ(j=β)的组合作用;若确定了船舶监测结构实际遭受两种形式外部载荷的组合作用,则所求得的k′sα、k′sβ即为与外部载荷Fsα、Fsβ对应的应力放大系数ksj(j=α,β),进而,由公式(4)可对应地计算得到船舶监测结构实际遭受的外部载荷Fsα、Fsβ的数值大小;由公式(6)可计算得到由传感器获得的、与外部载荷Fsj(j=α,β)对应的监测点ri位置的应力数据σ′ri(Fsj)。
8)若α、β遍历其全部可能的取值后,仍无法获得符合要求的外部载荷Fsα、Fsβ;若有m>3,则由方程(6)判断船舶监测结构是否同时遭受了三种形式外部载荷Fsα(j=α)、Fsβ(j=β)、Fsγ(j=γ)的组合作用;若确定了船舶监测结构实际遭受三种形式外部载荷的组合作用,则所求得的k′sα、k′sβ、k′sγ,即为与外部载荷Fsα、Fsβ、Fsγ对应的应力放大系数ksj(j=α,β,γ);进而,由公式(4)可对应地计算得到船舶监测结构实际遭受的外部载荷Fsα、Fsβ、Fsγ的数值大小;由公式(6)依次计算得到由传感器获得的、与外部载荷Fsj(j=α,β,γ)对应的监测点ri位置的应力数据σ′ri(Fsj)。
9)若α、β、γ遍历其全部可能的取值后,仍无法获得符合要求的外部载荷Fsα、Fsβ、Fsγ,则顺序地,采用同样的方法,依次判断船舶监测结构是否同时遭受了其余多种形式的外部载荷组合作用。若最终判断船舶监测结构此时也并非遭受了m-1种形式外部载荷的组合作用,则说明船舶监测结构同时遭受了全部形式外部载荷的组合作用。
10)在判定船舶监测结构同时遭受了全部形式外部载荷的组合作用后,由方程(7)求得与各种形式外部载荷Fsj对应的应力放大系数ksj;然后,由公式(4)可对应地计算得到船舶监测结构实际遭受的外部载荷Fsj的数值大小;进而,由公式(6)可依次计算得到由传感器获得的、与外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj)。
本发明给出了一种基于船舶结构监测数据分析的载荷识别方法,根据结构可能遭受的外部载荷作用的情况,采用船舶监测结构的监测点优化布置方法,最终确定各监测点位置;进而,通过由布置在监测点位置的传感器获得的应力数据,结合相对应的载荷识别方法,识别给出船舶监测结构实际遭受的外部载荷的作用情况。其具体特点如下:
1)该方法能实现对船舶监测结构实际遭受的外部载荷形式的准确、快速识别。该方法具体通过船舶监测结构的监测点优化布置方法、一种形式外部载荷单独作用对应的载荷识别方法、多种形式外部载荷组合作用对应的载荷识别方法、全部形式外部载荷组合作用对应的载荷识别方法,实现了对船舶监测结构实际遭受的外部载荷形式的准确、快速识别。
2)该方法简便易行,实用性强。该方法主要通过结合船舶监测结构的结构分析结果,采用船舶监测结构的监测点优化布置方法,确定出船舶监测结构各监测点的具体位置;然后,通过结合由传感器获取的应力数据,采用对应的载荷识别方法,可依次判断给出船舶监测结构实际遭受的外部载荷作用情况。
3)该方法有助于船舶监测结构整体应力状态的获取。结构实际遭受的外部载荷形式的准确、快速识别,是船舶监测结构整体应力状态获取的关键。因而,在由本发明中的方法准确识别出了结构实际遭受的外部载荷形式后,结合相关船舶监测结构整体应力状态的快速获取方法,即可快速地获得船舶监测结构的整体应力状态。
进一步限定,本发明方法对于船舶监测结构可能遭受的外部载荷形式超过两种的情况具有很好的适用性。
进一步限定,对于初步选定的m个监测点,应分别对应于在各种形式外部载荷的单位载荷作用下,船舶监测结构的结构应力最大位置。特别地,对于有n(n<m)种形式外部载荷对应的结构应力最大位置出现重合的情况,应考虑适当更换其中n-1种形式外部载荷对应的监测点位置。
进一步限定,在需要更换初步选定的监测点位置时,应按照与相应形式外部载荷对应的结构应力由大到小的顺序进行,以保证选取的监测点位置具有较好的代表性。
进一步限定,若判断出船舶监测结构未遭受全部形式外部载荷的作用,则应保留由与船舶监测结构实际遭受外部载荷对应的监测点位置的应力数据,并依此进行后续研究,而删除由其余各监测点位置的传感器获得的应力数据,以解决出现的数据冗余问题。
进一步限定,若在监测点位置布置了三向传感器,则对于专利中所涉及的由布置在监测点位置的传感器获得的应力数据、与每种形式外部载荷对应的监测点位置的单位应力数据,均应为监测点位置的合成应力。
实施例1:
一种基于船舶结构监测数据分析的载荷识别方法,包括以下步骤:
步骤1:根据船舶监测结构可能遭受m种形式外部载荷Fsj作用的情况,确定在船舶监测结构上布置监测点的数量m;j=1,2,…,m,m>2;
步骤2:对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;然后,分别对该离散模型施加上述各种形式外部载荷的单位载荷Fsj0,完成船舶监测结构的结构分析;结合结构分析结果,初步选定m个监测点位置ri;i=1,2,…,m;
步骤3:获取与每种形式外部载荷Fsj对应的各监测点位置的单位应力数据σ′ri(Fsj0),分别计算与各种形式外部载荷Fsj对应的监测点间的应力关系系数
步骤4:对于任意取定的p、q值,若存在则重新选择m个监测点位置ri,返回步骤3;若不存在/>则确定选定的m个监测点位置ri即为船舶监测结构的最终的监测点位置;
步骤5:获取布置在各监测点位置ri的传感器在遭受实际外部载荷的作用下的应力数据σ′ri(Fs),结合与每种形式外部载荷Fsj对应的监测点位置的单位应力数据σ′ri(Fsj0),计算在监测点ri位置与每种形式外部载荷Fsj对应的应力放大系数ksj-ri;
若存在ksj-r1=ksj-r2=…=ksj-rm,且对于任意n≠j,均不存在ksn-r1=ksn-r2=…=ksn-rm,则判定船舶监测结构此时仅遭受了外部载荷Fsj的作用;
若对于任意j,均不满足ksj-r1=ksj-r2=…=ksj-rm,则判定船舶监测结构此时遭受了多种形式外部载荷的组合作用;
步骤6:若船舶监测结构此时仅遭受了外部载荷Fsj的作用,则步骤5中计算得到的ksj-ri即为与该种形式外部载荷Fsj对应的应力放大系数ksj,则船舶监测结构实际遭受的外部载荷Fsj的数值大小为:Fsj=ksj·Fsj0;此时,由布置在监测点位置的传感器获得的应力数据σ′ri(Fs),即为与船舶监测结构实际遭受外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj);
步骤7:判断船舶监测结构同时遭受的外部载荷的数量x;
步骤7.1:初始化x=2;
步骤7.2:任取x种形式外部载荷判断船舶监测结构是否同时遭受了x种形式外部载荷/>的组合作用;
先由(1)、…、(x)联立求得与外部载荷对应的应力放大系数/>再分别由(x+1)、…、(m)验证所求的/>是否满足方程的等式关系;若方程(x+1)、…、(m)均满足等式关系,则判定船舶监测结构同时遭受了x种形式外部载荷/>的组合作用;
步骤7.3:若步骤7.2中α1、…、αx遍历其全部可能的取值后,仍无法判定船舶监测结构同时遭受了x种形式外部载荷的组合作用,且x<m-1,则令x=x+1,返回步骤7.2;若x=m-1,则判定船舶监测结构同时遭受了所有m种形式外部载荷的组合作用,令x=m,由步骤7.2中(1)、…、(m)联立求得与外部载荷/>对应的应力放大系数/>
步骤8:对于船舶监测结构同时遭受的x种形式外部载荷由步骤7.2中求得的/>即为与外部载荷/>对应的应力放大系数/>船舶监测结构实际遭受的x种形式外部载荷/>的数值大小为:
与船舶监测结构实际遭受的x种形式外部载荷对应的监测点ri位置的应力数据/>的计算方法为:
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于船舶结构监测数据分析的载荷识别方法,其特征在于,包括以下步骤:
步骤1:根据船舶监测结构可能遭受m种形式外部载荷Fsj作用的情况,确定在船舶监测结构上布置监测点的数量m;j=1,2,…,m,m>2;
步骤2:对船舶监测结构进行离散处理,得到以传感器监测范围为单元尺寸的离散模型;然后,分别对该离散模型施加上述各种形式外部载荷的单位载荷Fsj0,完成船舶监测结构的结构分析;结合结构分析结果,初步选定m个监测点位置ri;i=1,2,…,m;
步骤3:获取与每种形式外部载荷Fsj对应的各监测点位置的单位应力数据σ′ri(Fsj0),分别计算与各种形式外部载荷Fsj对应的监测点间的应力关系系数
步骤4:对于任意取定的p、q值,若存在则重新选择m个监测点位置ri,返回步骤3;若不存在/>则确定选定的m个监测点位置ri即为船舶监测结构的最终的监测点位置;
步骤5:获取布置在各监测点位置ri的传感器在遭受实际外部载荷的作用下的应力数据σ′ri(Fs),结合与每种形式外部载荷Fsj对应的监测点位置的单位应力数据σ′ri(Fsj0),计算在监测点ri位置与每种形式外部载荷Fsj对应的应力放大系数ksj-ri;
若存在ksj-r1=ksj-r2=…=ksj-rm,且对于任意n≠j,均不存在ksn-r1=ksn-r2=…=ksn-rm,则判定船舶监测结构此时仅遭受了外部载荷Fsj的作用;
若对于任意j,均不满足ksj-r1=ksj-r2=…=ksj-rm,则判定船舶监测结构此时遭受了多种形式外部载荷的组合作用;
步骤6:若船舶监测结构此时仅遭受了外部载荷Fsj的作用,则步骤5中计算得到的ksj-ri即为与该种形式外部载荷Fsj对应的应力放大系数ksj,则船舶监测结构实际遭受的外部载荷Fsj的数值大小为:Fsj=ksj·Fsj0;此时,由布置在监测点位置的传感器获得的应力数据σ′ri(Fs),即为与船舶监测结构实际遭受外部载荷Fsj对应的监测点ri位置的应力数据σ′ri(Fsj);
步骤7:判断船舶监测结构同时遭受的外部载荷的数量x;
步骤7.1:初始化x=2;
步骤7.2:任取x种形式外部载荷判断船舶监测结构是否同时遭受了x种形式外部载荷/>的组合作用;
先由(1)、…、(x)联立求得与外部载荷对应的应力放大系数/>再分别由(x+1)、…、(m)验证所求的/>是否满足方程的等式关系;若方程(x+1)、…、(m)均满足等式关系,则判定船舶监测结构同时遭受了x种形式外部载荷/>的组合作用;
步骤7.3:若步骤7.2中α1、…、αx遍历其全部可能的取值后,仍无法判定船舶监测结构同时遭受了x种形式外部载荷的组合作用,且x<m-1,则令x=x+1,返回步骤7.2;若x=m-1,则判定船舶监测结构同时遭受了所有m种形式外部载荷的组合作用,令x=m,由步骤7.2中(1)、…、(m)联立求得与外部载荷/>对应的应力放大系数/>
步骤8:对于船舶监测结构同时遭受的x种形式外部载荷由步骤7.2中求得的/>即为与外部载荷/>对应的应力放大系数/>船舶监测结构实际遭受的x种形式外部载荷/>的数值大小为:
与船舶监测结构实际遭受的x种形式外部载荷对应的监测点ri位置的应力数据/>的计算方法为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111509674.0A CN114199438B (zh) | 2021-12-10 | 2021-12-10 | 一种基于船舶结构监测数据分析的载荷识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111509674.0A CN114199438B (zh) | 2021-12-10 | 2021-12-10 | 一种基于船舶结构监测数据分析的载荷识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114199438A CN114199438A (zh) | 2022-03-18 |
CN114199438B true CN114199438B (zh) | 2023-07-21 |
Family
ID=80652414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111509674.0A Active CN114199438B (zh) | 2021-12-10 | 2021-12-10 | 一种基于船舶结构监测数据分析的载荷识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114199438B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114878049A (zh) * | 2022-05-20 | 2022-08-09 | 中交天津港湾工程研究院有限公司 | 一种间接测量海上大型浮体浮运中支墩受力方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9310962D0 (en) * | 1993-05-27 | 1993-07-14 | Shipboard Informatics Limited | Hull monitoring apparatus and method |
KR20160143414A (ko) * | 2015-06-05 | 2016-12-14 | 현대중공업 주식회사 | 선박 모니터링 방법 및 장치 |
CN108240881A (zh) * | 2018-01-05 | 2018-07-03 | 哈尔滨工程大学 | 一种船体总纵强度应力监测方法 |
CN110646184A (zh) * | 2019-09-23 | 2020-01-03 | 中国船舶重工集团公司第七一九研究所 | 一种船舶动力管道疲劳应力的监测方法 |
CN112214831A (zh) * | 2020-09-02 | 2021-01-12 | 中国船舶重工集团公司第七研究院 | 一种破冰船船体结构的冰载荷分析装置 |
-
2021
- 2021-12-10 CN CN202111509674.0A patent/CN114199438B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9310962D0 (en) * | 1993-05-27 | 1993-07-14 | Shipboard Informatics Limited | Hull monitoring apparatus and method |
KR20160143414A (ko) * | 2015-06-05 | 2016-12-14 | 현대중공업 주식회사 | 선박 모니터링 방법 및 장치 |
CN108240881A (zh) * | 2018-01-05 | 2018-07-03 | 哈尔滨工程大学 | 一种船体总纵强度应力监测方法 |
CN110646184A (zh) * | 2019-09-23 | 2020-01-03 | 中国船舶重工集团公司第七一九研究所 | 一种船舶动力管道疲劳应力的监测方法 |
CN112214831A (zh) * | 2020-09-02 | 2021-01-12 | 中国船舶重工集团公司第七研究院 | 一种破冰船船体结构的冰载荷分析装置 |
Also Published As
Publication number | Publication date |
---|---|
CN114199438A (zh) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114199438B (zh) | 一种基于船舶结构监测数据分析的载荷识别方法 | |
US11842575B2 (en) | Method and system for vehicle analysis | |
CN109543143B (zh) | 非线性带偏差系统的多传感器融合估计方法 | |
JPH04312184A (ja) | 画像処理システムのエッジ検出方法及びその装置 | |
CN113722985A (zh) | 航空发动机健康状态评估和剩余寿命预测方法及系统 | |
CN112633567A (zh) | 一种等待时长的预测方法、设备及存储介质 | |
CN108985517A (zh) | 基于线性回归的短时交通流量预测方法 | |
CN114216599B (zh) | 一种基于传感器数据分析的船舶监测结构载荷识别方法 | |
CN114048642A (zh) | 一种航空发动机性能趋势分析方法 | |
CN116070985B (zh) | 危化品车辆装卸载过程识别方法 | |
CN114199436B (zh) | 一种基于传感器数据的船舶监测结构整体应力状态快速获取方法 | |
CN109567783B (zh) | 一种基于加权马尔可夫模型预测空间位置心磁信号的方法 | |
CN110598274A (zh) | 一种大型复杂结构损伤诊断的多主体与信息融合方法 | |
CN113722833B (zh) | 基于双通道长短时记忆网络的涡扇发动机剩余使用寿命预测方法 | |
CN115206455A (zh) | 基于深度神经网络的稀土元素组分含量预测方法及系统 | |
CN109946096B (zh) | 基于模型空间的高速列车风管故障诊断方法 | |
CN114383845A (zh) | 一种基于嵌入式零样本学习模型的轴承复合故障诊断方法 | |
CN113761795A (zh) | 一种飞机发动机故障检测方法和系统 | |
El-Badawy et al. | Correlation between different DNA period-3 signals: An analytical study for exons prediction | |
CN113902161A (zh) | 一种排水管网溢流预警智能分析的方法 | |
CN112487361A (zh) | 一种基于arima和小波变换的扇区短时交通流量预测方法 | |
CN111263326A (zh) | 基于多重融合卷积gru的车辆位置预测方法 | |
Konovalov et al. | Queueing systems with renovation vs. queues with red. supplementary material | |
CN118114118A (zh) | 一种基于cndt的典型武器装备故障诊断方法 | |
Han et al. | A study on the prediction of mechanical properties of titanium alloy based on adaptive fuzzy-neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |