CN116415463A - Analysis method of ship simulation status data - Google Patents

Analysis method of ship simulation status data Download PDF

Info

Publication number
CN116415463A
CN116415463A CN202310410457.9A CN202310410457A CN116415463A CN 116415463 A CN116415463 A CN 116415463A CN 202310410457 A CN202310410457 A CN 202310410457A CN 116415463 A CN116415463 A CN 116415463A
Authority
CN
China
Prior art keywords
ship
wave
roll
rolling
angle
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.)
Pending
Application number
CN202310410457.9A
Other languages
Chinese (zh)
Inventor
宋鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Shipbuilding Intelligent Technology Shanghai Co ltd
Original Assignee
China Shipbuilding Intelligent Technology Shanghai Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Shipbuilding Intelligent Technology Shanghai Co ltd filed Critical China Shipbuilding Intelligent Technology Shanghai Co ltd
Priority to CN202310410457.9A priority Critical patent/CN116415463A/en
Publication of CN116415463A publication Critical patent/CN116415463A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B71/00Designing vessels; Predicting their performance
    • B63B71/10Designing vessels; Predicting their performance using computer simulation, e.g. finite element method [FEM] or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及船舶仿真技术领域,尤其涉及船舶仿真状态数据分析方法,采用计算机建模软件搭建船舶有限元模型;获取船舶有限元模型的各节点坐标信息,建立船舶三维仿真图像;将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标;根据船舶排水量及浮心坐标获取船舶横摇回复力矩,从而计算出船舶的瞬时回复力;获取船舶海浪力,根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,以获取横摇概率及某一概率值下的横摇角数据将所获取的横摇角与事先设定的进水角进行对比,以判断航行状态;本发明通过仿真计算出船舶横摇角及横摇概率,以得到可靠的船舶运动下横摇角及横摇概率,预防降低船舶在航行中颠覆风险。

Figure 202310410457

The present invention relates to the field of ship simulation technology, in particular to a method for analyzing ship simulation state data, which uses computer modeling software to build a ship finite element model; obtains the coordinate information of each node of the ship finite element model, and establishes a three-dimensional simulation image of the ship; Convert the coordinates to the earth coordinate system to obtain the displacement of the ship and its buoyancy coordinates; obtain the rolling recovery moment of the ship according to the displacement of the ship and the coordinates of the buoyancy center, thereby calculating the instantaneous recovery force of the ship; Calculating the time history of rolling recovery moment and rolling motion time history to obtain the rolling probability and the rolling angle data under a certain probability value. Comparison is made to judge the sailing state; the present invention calculates the roll angle and roll probability of the ship through simulation, so as to obtain reliable roll angle and roll probability under ship motion, and prevent and reduce the risk of ship overturning during navigation.

Figure 202310410457

Description

船舶仿真状态数据分析方法Analysis method of ship simulation state data

技术领域technical field

本发明涉及船舶仿真技术领域,尤其涉及船舶仿真状态数据分析方法。The invention relates to the technical field of ship simulation, in particular to a method for analyzing state data of ship simulation.

背景技术Background technique

船舶作业作为目前重要的运输交通工具,其在海中的运输性能尤为重要,海浪环境会对船舶的横荡、纵荡、垂荡、横摇、纵摇和首摇六个自由度的运动产生干扰,在这六种运动中,只有垂荡、横摇、纵摇是纯粹的振荡运动,与另外三种运动区别在于,当扰动过大致使船舶偏离平衡位置时,震荡运动会受到复原力和力矩的作用。这种纯粹的振荡运动会产生较大的加速度,极易导致船载设备失灵及船员晕船现象,严重制约了船舶作业率,随着技术发展,对获取大量的船舶运动数据也变成现实,通过仿真船舶运动以解决或避免上述问题,在国内外的研究中,横摇角及其对应的横摇概率尤为重要,本申请针对该问题设计了一种船舶仿真状态数据分析方法,以预防降低船舶横摇运动中颠覆问题。Ship operation is an important means of transportation at present, and its transportation performance in the sea is particularly important. The sea wave environment will interfere with the six degrees of freedom of the ship's sway, surge, heave, roll, pitch and yaw motion. , among these six motions, only heave, roll, and pitch are purely oscillating motions. The difference from the other three motions is that when the disturbance is too large to make the ship deviate from the equilibrium position, the oscillating motion will be affected by the restoring force and moment. effect. This kind of pure oscillating motion will produce a large acceleration, which can easily lead to the failure of shipboard equipment and seasickness of the crew, which seriously restricts the ship's operating rate. With the development of technology, it has become a reality to obtain a large amount of ship motion data. Through simulation In order to solve or avoid the above problems, the roll angle and its corresponding roll probability are particularly important in the research at home and abroad. This application designs a ship simulation state data analysis method for this problem, so as to prevent the reduction of ship roll. Subvert the problem in the shake motion.

发明内容Contents of the invention

(一)发明目的(1) Purpose of the invention

为解决背景技术中存在的技术问题,本发明提出船舶仿真状态数据分析方法。In order to solve the technical problems existing in the background technology, the present invention proposes a ship simulation state data analysis method.

(二)技术方案(2) Technical solution

为解决上述问题,本发明提供了船舶仿真状态数据分析方法,包括:In order to solve the above problems, the invention provides a ship simulation state data analysis method, including:

采用计算机建模软件搭建船舶有限元模型;Using computer modeling software to build a ship finite element model;

获取船舶有限元模型的各节点坐标信息,建立船舶三维仿真图像;Obtain the coordinate information of each node of the ship finite element model, and establish a three-dimensional simulation image of the ship;

将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标;Converting the coordinates of each node to the earth coordinate system to obtain the displacement of the ship and the coordinates of its center of buoyancy;

根据船舶排水量及浮心坐标获取船舶横摇回复力矩,从而计算出船舶的瞬时回复力;According to the displacement of the ship and the coordinates of the buoyancy center, the restoring moment of the ship's roll is obtained, so as to calculate the instantaneous restoring force of the ship;

获取船舶海浪力,根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,以获取横摇概率及某一概率值下的横摇角数据;Obtain the ship's sea wave force, and calculate the roll restoration moment time history and the roll motion time history according to the instantaneous restoration force and the sea wave force, so as to obtain the roll probability and the roll angle data under a certain probability value;

将所获取的横摇角与事先设定的进水角进行对比,以判断航行状态。Compare the obtained roll angle with the preset water ingress angle to judge the sailing state.

作为本发明的一个实施方案,将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标,包括:As an embodiment of the present invention, the coordinates of each node are converted to a geodetic coordinate system to obtain the displacement of the ship and its buoyancy center coordinates, including:

沿海浪传递方向经转换至大地坐标系中的船舶节点划分至多个剖面中,所述剖面与海浪传递方向对应;The ship nodes converted to the earth coordinate system along the wave transmission direction are divided into multiple sections, and the sections correspond to the wave transmission direction;

获取船舶在海浪传递方向上的投影长度,根据所述投影长度计算出每个剖面面积及形心坐标信息。Obtain the projected length of the ship in the wave transmission direction, and calculate the area of each section and centroid coordinate information according to the projected length.

作为本发明的一个实施方案,将所述投影长度n等比进行分割,并设定每个被分割的投影长度为m,根据以下公式计算排水量V:As an embodiment of the present invention, the projection length n is divided into equal proportions, and each divided projection length is set as m, and the displacement V is calculated according to the following formula:

V=(S1+S2+……Sk)*m。V=(S1+S2+...Sk)*m.

作为本发明的一个实施方案,所述获取船舶海浪力,包括:As an embodiment of the present invention, the acquisition of ship sea wave force includes:

将船舶划分为左右对称的两个部分,记作左部分和右部分;Divide the ship into two symmetrical parts, denoted as the left part and the right part;

任取左部分和右部分设置为多个横剖面,并将所有横剖面按循序组合成平行四边形结构,对称化至未被选取的船舶部分,以获取整体船舶的所有平行四边形结构,并对所有平行四边形结构进行坐标转换。The left part and the right part are randomly selected as multiple cross-sections, and all cross-sections are sequentially combined into a parallelogram structure, and symmetrized to the unselected part of the ship to obtain all the parallelogram structures of the overall ship, and for all Parallelogram structure for coordinate transformation.

作为本发明的一个实施方案,对所述平行四边形结构进行三角形结构处理,将每个平行四边形结构拆分为两个三角形结构,对每个三角形结构进行进行计算,得出每个三角形结构对应的海浪力;As an embodiment of the present invention, the parallelogram structure is subjected to triangular structure processing, each parallelogram structure is split into two triangular structures, and each triangular structure is calculated to obtain the corresponding wave force;

对所有的三角形结构对应的海浪力积分整合,获取整体船舶的海浪力。Integrate and integrate the wave force corresponding to all the triangle structures to obtain the wave force of the whole ship.

作为本发明的一个实施方案,所述根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,包括:As an embodiment of the present invention, the calculation of the rolling recovery moment time history and the rolling motion time history according to the instantaneous restoring force and the sea wave force includes:

利用达朗贝尔原理计算出横摇回复力矩时程和横摇运动时程所对应的曲线;Calculate the curves corresponding to the time history of rolling restoring moment and rolling motion time history by using d'Alembert's principle;

跟所述横摇回复力矩时程和横摇运动时程曲线进行离散化分析,获取横摇概率。A discretization analysis is performed with the rolling recovery moment time history and the rolling motion time history curve to obtain the rolling probability.

作为本发明的一个实施方案,还包括:As an embodiment of the present invention, also includes:

构建随机海浪模型,对船舶的波倾角进行仿真;Build a random wave model to simulate the ship's wave inclination;

构建船舶横摇仿真模型,利用卷积定理及拉普朗反变换计算出船舶横摇角。A ship roll simulation model is constructed, and the ship roll angle is calculated by using the convolution theorem and Laplan inverse transform.

作为本发明的一个实施方案,所述构建随机海浪模型,对船舶的波倾角进行仿真,g为重力加速度,其中波高的表达公式如下:As an embodiment of the present invention, the random sea wave model is constructed to simulate the wave inclination angle of the ship, and g is the acceleration of gravity, wherein the expression formula of the wave height is as follows:

Figure BDA0004183044220000021
Figure BDA0004183044220000021

海浪谱函数P(k)表示为:The wave spectral function P(k) is expressed as:

Figure BDA0004183044220000022
其中J=0.78,/>
Figure BDA0004183044220000023
Figure BDA0004183044220000022
where J=0.78, />
Figure BDA0004183044220000023

波倾角能量谱密度函数Pa(k);Wave inclination energy spectral density function Pa(k);

Figure BDA0004183044220000031
则海浪的波高仿真为:
Figure BDA0004183044220000031
Then the wave height simulation of the ocean wave is:

Figure BDA0004183044220000032
设最大海浪有效波倾角为α=kUα
Figure BDA0004183044220000032
Assume that the maximum effective wave inclination of sea waves is α=kU α ,

那么波倾角的表达式为:Then the expression of wave inclination is:

Figure BDA0004183044220000033
Figure BDA0004183044220000033

作为本发明的一个实施方案,还包括显示端,用于将建立船舶三维仿真图像可视化数据显示。As an embodiment of the present invention, it also includes a display terminal for displaying the visualized data of the three-dimensional simulation image of the ship.

作为本发明的一个实施方案,所述根据船舶排水量及浮心坐标,还可以获取纵摇角和垂荡角;As an embodiment of the present invention, the pitch angle and heave angle can also be obtained according to the displacement of the ship and the coordinates of the center of buoyancy;

根据横摇角、纵摇角以及垂荡角计算相对误差,根据所述相对误差确定出误差峰值,根据所述误差峰值确定误差缩减最佳比例;Calculate the relative error according to the roll angle, pitch angle and heave angle, determine the error peak value according to the relative error, and determine the optimal ratio of error reduction according to the error peak value;

根据所述最佳缩减比例可相应调整缩减各输入参数数据。According to the optimal reduction ratio, each input parameter data can be adjusted and reduced accordingly.

本发明的上述技术方案具有如下有益的技术效果:The technical solution of the present invention has the following beneficial technical effects:

本发明通过仿真计算出船舶横摇角及横摇概率,通过海浪随机模型仿真出船舶的波倾角,以及产生船舶横摇时间历程和概率分布,以得到可靠的船舶运动下横摇角及横摇概率,预防降低船舶在航行中颠覆风险。The invention calculates the roll angle and roll probability of the ship through simulation, simulates the wave inclination angle of the ship through the random wave model, and generates the time history and probability distribution of the ship's roll, so as to obtain reliable roll angle and roll under ship motion. probability, prevent and reduce the risk of ship subversion during navigation.

附图说明Description of drawings

图1为本发明方法流程图;Fig. 1 is a flow chart of the method of the present invention;

图2为本发明方法中一实施例流程图。Fig. 2 is a flowchart of an embodiment of the method of the present invention.

具体实施方式Detailed ways

为使本发明的目的、技术方案和优点更加清楚明了,下面结合具体实施方式并参照附图,对本发明进一步详细说明。应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below in combination with specific embodiments and with reference to the accompanying drawings. It should be understood that these descriptions are exemplary only, and are not intended to limit the scope of the present invention. Also, in the following description, descriptions of well-known structures and techniques are omitted to avoid unnecessarily obscuring the concept of the present invention.

在航行的过程中,海浪以及船舶自身原因对船舶产生影响,会导致船舶出现复杂的摇摆运动,这在一定程度上会降低船舶的实用性和航海性能。而船舶的横摇运动时船舶坐标中六个自由度中对船舶影响较大的因素,因此船舶运动数学模型的建立对船舶横摇运动的一直和预报起到重要作用。During the sailing process, the waves and the ship's own reasons will affect the ship, which will cause the ship to appear complex rocking motions, which will reduce the practicability and navigating performance of the ship to a certain extent. The rolling motion of the ship is the factor that has a greater influence on the ship among the six degrees of freedom in the ship coordinates. Therefore, the establishment of the mathematical model of the ship's motion plays an important role in the constant and forecast of the ship's rolling motion.

如图1所示,本发明提出的船舶仿真状态数据分析方法,包括如下步骤:As shown in Figure 1, the ship simulation state data analysis method that the present invention proposes, comprises the following steps:

步骤S100,采用计算机建模软件搭建船舶有限元模型;步骤S200,获取船舶有限元模型的各节点坐标信息,建立船舶三维仿真图像;步骤S300,将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标;步骤S400,根据船舶排水量及浮心坐标获取船舶横摇回复力矩,从而计算出船舶的瞬时回复力;步骤S500,获取船舶海浪力,根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,以获取横摇概率及某一概率值下的横摇角数据;步骤S600,将所获取的横摇角与事先设定的进水角进行对比,以判断航行状态。作为示例而言,横摇角小于进水角,船舶正常航行;假设横摇角大于进水角则调整船舶航行状态,使得横摇角小于进水角。Step S100, using computer modeling software to build a ship finite element model; step S200, obtaining the coordinate information of each node of the ship finite element model, and establishing a three-dimensional simulation image of the ship; step S300, converting the coordinates of each node to the earth coordinate system, and obtaining The displacement of the ship and its center of buoyancy coordinates; step S400, obtain the rolling recovery moment of the ship according to the displacement of the ship and the coordinates of the center of buoyancy, thereby calculating the instantaneous recovery force of the ship; step S500, obtain the sea wave force of the ship, according to the instantaneous recovery force and The sea wave force calculates the time history of the rolling recovery moment and the rolling motion time history to obtain the rolling probability and the rolling angle data under a certain probability value; step S600, the obtained rolling angle is compared with the preset progress The water angle is compared to judge the navigation status. As an example, if the rolling angle is smaller than the flooding angle, the ship sails normally; if the rolling angle is greater than the flooding angle, then the sailing state of the ship is adjusted so that the rolling angle is smaller than the flooding angle.

或者是,假设此时横摇概率为3%,那么即可进一步降低横摇概率。Alternatively, assuming that the rolling probability is 3% at this time, the rolling probability can be further reduced.

在本实施例中,所述将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标,包括:沿海浪传递方向经转换至大地坐标系中的船舶节点划分至多个剖面中,所述剖面与海浪传递方向对应;获取船舶在海浪传递方向上的投影长度,根据所述投影长度计算出每个剖面面积及形心坐标信息。In this embodiment, the conversion of the coordinates of each node into the earth coordinate system to obtain the displacement of the ship and its buoyancy center coordinates includes: dividing the ship nodes into multiple sections along the direction of wave transfer into the earth coordinate system Among them, the section corresponds to the direction of wave transmission; the projected length of the ship in the direction of wave transmission is obtained, and the area and centroid coordinate information of each section is calculated according to the projected length.

进一步的,将所述投影长度n等比进行分割,并设定每个被分割的投影长度为m,根据以下公式计算排水量V:Further, the projection length n is divided into equal ratios, and each divided projection length is set to m, and the displacement V is calculated according to the following formula:

V=(S1+S2+……Sk)*m。V=(S1+S2+...Sk)*m.

参考图2,所述获取船舶海浪力,包括:将船舶划分为左右对称的两个部分,记作左部分和右部分;任取左部分和右部分设置为多个横剖面,并将所有横剖面按循序组合成平行四边形结构,对称化至未被选取的船舶部分,以获取整体船舶的所有平行四边形结构,并对所有平行四边形结构进行坐标转换。With reference to Fig. 2, the described acquisition of the ship's sea wave force includes: dividing the ship into two symmetrical parts, which are denoted as the left part and the right part; Sections are sequentially combined into parallelogram structures, symmetrized to unselected ship parts to obtain all parallelogram structures of the overall ship, and coordinate transformation is performed on all parallelogram structures.

可以理解的是,对所述平行四边形结构进行三角形结构处理,将每个平行四边形结构拆分为两个三角形结构,对每个三角形结构进行进行计算,得出每个三角形结构对应的海浪力;对所有的三角形结构对应的海浪力积分整合,获取整体船舶的海浪力。It can be understood that the triangular structure processing is performed on the parallelogram structure, and each parallelogram structure is split into two triangular structures, and each triangular structure is calculated to obtain the wave force corresponding to each triangular structure; Integrate and integrate the wave force corresponding to all the triangle structures to obtain the wave force of the whole ship.

所述根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,包括:利用达朗贝尔原理计算出横摇回复力矩时程和横摇运动时程所对应的曲线;跟所述横摇回复力矩时程和横摇运动时程曲线进行离散化分析,获取横摇概率。The calculation of the roll recovery moment time history and the roll motion time history according to the instantaneous recovery force and the sea wave force includes: calculating the roll recovery torque time history and the roll motion time history corresponding to the D’Alembert principle Curve; discretization analysis is carried out with the rolling recovery torque time history and rolling motion time history curve to obtain the rolling probability.

本发明还包括:The present invention also includes:

构建随机海浪模型,对船舶的波倾角进行仿真;构建船舶横摇仿真模型,利用卷积定理及拉普朗反变换计算出船舶横摇角。Build a random wave model to simulate the ship's wave inclination; build a ship roll simulation model, and use the convolution theorem and Laplan's inverse transform to calculate the ship's roll angle.

由于海浪的干扰是不规则的所以海浪产生的波长以及波高、周期都是随机的,因此在本实施例中假设随机波浪的构成由波幅、波长和随机相位均不一样的波元相互累加而成,故而可获取到随机波高表达方程如下。Because the interference of ocean waves is irregular, the wavelength, wave height, and cycle of ocean waves are random. Therefore, in this embodiment, it is assumed that random waves are composed of wave elements with different amplitudes, wavelengths, and random phases. , so the random wave height expression equation can be obtained as follows.

所述构建随机海浪模型,对船舶的波倾角进行仿真,g为重力加速度,其中波高的表达公式如下:The random sea wave model is constructed to simulate the wave inclination of the ship, g is the acceleration of gravity, and the expression formula of the wave height is as follows:

Figure BDA0004183044220000051
Figure BDA0004183044220000051

海浪谱函数P(k)表示为:The wave spectral function P(k) is expressed as:

Figure BDA0004183044220000052
其中J=0.78,/>
Figure BDA0004183044220000053
Figure BDA0004183044220000052
where J=0.78, />
Figure BDA0004183044220000053

波倾角能量谱密度函数Pa(k);Wave inclination energy spectral density function Pa(k);

Figure BDA0004183044220000054
则海浪的波高仿真为:
Figure BDA0004183044220000054
Then the wave height simulation of the ocean wave is:

Figure BDA0004183044220000055
设最大海浪有效波倾角为α=kUα
Figure BDA0004183044220000055
Assume that the maximum effective wave inclination of sea waves is α=kU α ,

那么波倾角的表达式为:Then the expression of wave inclination is:

Figure BDA0004183044220000056
Figure BDA0004183044220000056

但在实际航海中,海况个船舶入水深度以及船的自身体积形状都会对波倾角产生影响,因此,在实际的计算中需要对上述公式进行修正,其具体操作视需要而定,在此不做赘述。However, in actual navigation, the sea conditions, the depth of the ship's entry into the water and the shape of the ship's own volume will all affect the wave inclination angle. Therefore, the above formula needs to be corrected in the actual calculation. The specific operation depends on the need, and will not be done here. repeat.

进一步可以理解的是,本发明某些实施例中还包括显示端,用于将建立船舶三维仿真图像可视化数据显示。It can be further understood that some embodiments of the present invention also include a display terminal for displaying the visualized data of the three-dimensional simulation image of the ship.

本发明的上述步骤S400中,所述根据船舶排水量及浮心坐标,还可以获取纵摇角和垂荡角;根据横摇角、纵摇角以及垂荡角计算相对误差,根据所述相对误差确定出误差峰值,根据所述误差峰值确定误差缩减最佳比例;根据所述最佳缩减比例可相应调整缩减各输入参数数据。In the above step S400 of the present invention, the pitch angle and the heave angle can also be obtained according to the displacement of the ship and the coordinates of the buoyancy center; the relative error is calculated according to the roll angle, the pitch angle, and the heave angle, and according to the relative error The error peak value is determined, and the optimal ratio of error reduction is determined according to the error peak value; according to the optimal reduction ratio, each input parameter data can be adjusted and reduced accordingly.

具体的,建立垂荡角集合、横摇角集合以及纵摇角集合。Specifically, a heave angle set, a roll angle set, and a pitch angle set are established.

根据上述集合获取峰值误差,由确定的峰值误差进行参数缩减比例确认。The peak error is obtained according to the above collection, and the parameter reduction ratio is confirmed by the determined peak error.

以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的仅为本发明的优选例,并不用来限制本发明,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。The basic principles, main features and advantages of the present invention have been shown and described above. Those skilled in the art should understand that the present invention is not limited by the above-mentioned embodiments, and those described in the above-mentioned embodiments and description are only preferred examples of the present invention, and are not intended to limit the present invention, without departing from the spirit and scope of the present invention. Under the premise, the present invention will have various changes and improvements, and these changes and improvements all fall within the scope of the claimed invention.

Claims (10)

1.船舶仿真状态数据分析方法,其特征在于,包括:1. The ship simulation state data analysis method is characterized in that, comprising: 采用计算机建模软件搭建船舶有限元模型;Using computer modeling software to build a ship finite element model; 获取船舶有限元模型的各节点坐标信息,建立船舶三维仿真图像;Obtain the coordinate information of each node of the ship finite element model, and establish a three-dimensional simulation image of the ship; 将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标;Converting the coordinates of each node to the earth coordinate system to obtain the displacement of the ship and the coordinates of its center of buoyancy; 根据船舶排水量及浮心坐标获取船舶横摇回复力矩,从而计算出船舶的瞬时回复力;According to the displacement of the ship and the coordinates of the buoyancy center, the restoring moment of the ship's roll is obtained, so as to calculate the instantaneous restoring force of the ship; 获取船舶海浪力,根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,以获取横摇概率及某一概率值下的横摇角数据;Obtain the ship's sea wave force, and calculate the roll restoration moment time history and the roll motion time history according to the instantaneous restoration force and the sea wave force, so as to obtain the roll probability and the roll angle data under a certain probability value; 将所获取的横摇角与事先设定的进水角进行对比,以判断航行状态。Compare the obtained roll angle with the preset water ingress angle to judge the sailing state. 2.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,所述将所述各节点坐标转换至大地坐标系,获取船舶的排水量及其浮心坐标,包括:2. the ship simulation state data analysis method according to claim 1, is characterized in that, described described each node coordinate is converted to the earth coordinate system, obtains the displacement of ship and its center of buoyancy coordinate, comprising: 沿海浪传递方向经转换至大地坐标系中的船舶节点划分至多个剖面中,所述剖面与海浪传递方向对应;The ship nodes converted to the earth coordinate system along the wave transmission direction are divided into multiple sections, and the sections correspond to the wave transmission direction; 获取船舶在海浪传递方向上的投影长度,根据所述投影长度计算出每个剖面面积及形心坐标信息。Obtain the projected length of the ship in the wave transmission direction, and calculate the area of each section and centroid coordinate information according to the projected length. 3.根据权利要求2所述的船舶仿真状态数据分析方法,其特征在于,将所述投影长度n等比进行分割,并设定每个被分割的投影长度为m,根据以下公式计算排水量V:3. The ship simulation state data analysis method according to claim 2, characterized in that, the projection length n is divided into equal ratios, and each divided projection length is set as m, and the displacement V is calculated according to the following formula : V=(S1+S2+……Sk)*m。V=(S1+S2+...Sk)*m. 4.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,所述获取船舶海浪力,包括:4. the ship simulation state data analysis method according to claim 1, is characterized in that, described obtaining ship sea wave force, comprises: 将船舶划分为左右对称的两个部分,记作左部分和右部分;Divide the ship into two symmetrical parts, denoted as the left part and the right part; 任取左部分和右部分设置为多个横剖面,并将所有横剖面按循序组合成平行四边形结构,对称化至未被选取的船舶部分,以获取整体船舶的所有平行四边形结构,并对所有平行四边形结构进行坐标转换。The left part and the right part are randomly selected as multiple cross-sections, and all cross-sections are sequentially combined into a parallelogram structure, and symmetrized to the unselected part of the ship to obtain all the parallelogram structures of the overall ship, and for all Parallelogram structure for coordinate transformation. 5.根据权利要求4所述的船舶仿真状态数据分析方法,其特征在于,对所述平行四边形结构进行三角形结构处理,将每个平行四边形结构拆分为两个三角形结构,对每个三角形结构进行进行计算,得出每个三角形结构对应的海浪力;5. ship simulation state data analysis method according to claim 4, is characterized in that, described parallelogram structure is carried out triangular structure processing, each parallelogram structure is split into two triangular structures, for each triangular structure Perform calculations to obtain the wave force corresponding to each triangle structure; 对所有的三角形结构对应的海浪力积分整合,获取整体船舶的海浪力。Integrate and integrate the wave force corresponding to all the triangle structures to obtain the wave force of the whole ship. 6.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,所述根据所述瞬时回复力和海浪力计算出横摇回复力矩时程和横摇运动时程,包括:6. The ship simulation state data analysis method according to claim 1, characterized in that, calculating the rolling recovery moment time history and rolling motion time history according to the instantaneous restoring force and sea wave force, comprising: 利用达朗贝尔原理计算出横摇回复力矩时程和横摇运动时程所对应的曲线;Calculate the curves corresponding to the time history of rolling restoring moment and rolling motion time history by using d'Alembert's principle; 跟所述横摇回复力矩时程和横摇运动时程曲线进行离散化分析,获取横摇概率。A discretization analysis is performed with the rolling recovery moment time history and the rolling motion time history curve to obtain the rolling probability. 7.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,还包括:7. ship simulation state data analysis method according to claim 1, is characterized in that, also comprises: 构建随机海浪模型,对船舶的波倾角进行仿真;Build a random wave model to simulate the ship's wave inclination; 构建船舶横摇仿真模型,利用卷积定理及拉普朗反变换计算出船舶横摇角。A ship roll simulation model is constructed, and the ship roll angle is calculated by using the convolution theorem and Laplan inverse transform. 8.根据权利要求7所述的船舶仿真状态数据分析方法,其特征在于,所述构建随机海浪模型,对船舶的波倾角进行仿真,g为重力加速度,其中波高的表达公式如下:8. ship simulation state data analysis method according to claim 7, it is characterized in that, described building random sea wave model, the wave inclination angle of ship is simulated, and g is gravitational acceleration, and wherein the expression formula of wave height is as follows:
Figure FDA0004183044190000021
Figure FDA0004183044190000021
海浪谱函数P(k)表示为:The wave spectral function P(k) is expressed as:
Figure FDA0004183044190000022
其中/>
Figure FDA0004183044190000023
Figure FDA0004183044190000022
where />
Figure FDA0004183044190000023
波倾角能量谱密度函数Pa(k);Wave inclination energy spectral density function Pa(k);
Figure FDA0004183044190000024
则海浪的波高仿真为:
Figure FDA0004183044190000024
Then the wave height simulation of the ocean wave is:
Figure FDA0004183044190000025
设最大海浪有效波倾角为α=kUα
Figure FDA0004183044190000025
Assume that the maximum effective wave inclination of sea waves is α=kU α ,
那么波倾角的表达式为:Then the expression of wave inclination is:
Figure FDA0004183044190000026
Figure FDA0004183044190000026
9.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,还包括显示端,用于将建立船舶三维仿真图像可视化数据显示。9 . The method for analyzing ship simulation state data according to claim 1 , further comprising a display terminal for visually displaying the established three-dimensional ship simulation image data. 10 . 10.根据权利要求1所述的船舶仿真状态数据分析方法,其特征在于,所述根据船舶排水量及浮心坐标,还可以获取纵摇角和垂荡角;10. the ship simulation state data analysis method according to claim 1, is characterized in that, described according to ship displacement and buoyancy coordinates, can also obtain pitch angle and heave angle; 根据横摇角、纵摇角以及垂荡角计算相对误差,根据所述相对误差确定出误差峰值,根据所述误差峰值确定误差缩减最佳比例;Calculate the relative error according to the roll angle, pitch angle and heave angle, determine the error peak value according to the relative error, and determine the optimal ratio of error reduction according to the error peak value; 根据所述最佳缩减比例可相应调整缩减各输入参数数据。According to the optimal reduction ratio, each input parameter data can be adjusted and reduced accordingly.
CN202310410457.9A 2023-04-17 2023-04-17 Analysis method of ship simulation status data Pending CN116415463A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310410457.9A CN116415463A (en) 2023-04-17 2023-04-17 Analysis method of ship simulation status data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310410457.9A CN116415463A (en) 2023-04-17 2023-04-17 Analysis method of ship simulation status data

Publications (1)

Publication Number Publication Date
CN116415463A true CN116415463A (en) 2023-07-11

Family

ID=87054449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310410457.9A Pending CN116415463A (en) 2023-04-17 2023-04-17 Analysis method of ship simulation status data

Country Status (1)

Country Link
CN (1) CN116415463A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117390762A (en) * 2023-09-15 2024-01-12 中国人民解放军92942部队 Full-flow ship safety design analysis method
CN117556623A (en) * 2023-11-21 2024-02-13 中国舰船研究设计中心 Digital twinning-based large ship navigation kinematic model construction method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103387038A (en) * 2013-07-30 2013-11-13 大连理工大学 Analysis method for reducing rolling motion of ship
CN106772285A (en) * 2016-11-29 2017-05-31 公安部第三研究所 The preprocess method of boat-carrying X-band wave observation radar echo
CN108153976A (en) * 2017-12-25 2018-06-12 重庆华渝电气集团有限公司 The emulation mode and host computer of naval vessel rolling motion in a kind of wave
RU2720122C1 (en) * 2019-07-30 2020-04-24 Петр Геннадьевич Тенишев Method for automated control of ship movement at sea waves in order to reduce on-board and pitch motion and system implementing method (embodiments)

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103387038A (en) * 2013-07-30 2013-11-13 大连理工大学 Analysis method for reducing rolling motion of ship
CN106772285A (en) * 2016-11-29 2017-05-31 公安部第三研究所 The preprocess method of boat-carrying X-band wave observation radar echo
CN108153976A (en) * 2017-12-25 2018-06-12 重庆华渝电气集团有限公司 The emulation mode and host computer of naval vessel rolling motion in a kind of wave
RU2720122C1 (en) * 2019-07-30 2020-04-24 Петр Геннадьевич Тенишев Method for automated control of ship movement at sea waves in order to reduce on-board and pitch motion and system implementing method (embodiments)

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘鹏;籍艳;漆随平;王东明;宋文杰;孙佳;: "基于MATLAB的海浪及船舶横摇仿真模型研究", 山东科学, no. 06, 15 December 2012 (2012-12-15), pages 1 - 3 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117390762A (en) * 2023-09-15 2024-01-12 中国人民解放军92942部队 Full-flow ship safety design analysis method
CN117390762B (en) * 2023-09-15 2024-05-14 中国人民解放军92942部队 Full-flow ship safety design analysis method
CN117556623A (en) * 2023-11-21 2024-02-13 中国舰船研究设计中心 Digital twinning-based large ship navigation kinematic model construction method
CN117556623B (en) * 2023-11-21 2024-04-09 中国舰船研究设计中心 Digital twinning-based large ship navigation kinematic model construction method

Similar Documents

Publication Publication Date Title
CN116415463A (en) Analysis method of ship simulation status data
CN103387038B (en) Analytical Method for Reducing Rolling Motion of Ship
CN113010958B (en) Simulation system of self-propelled ship and operation method thereof
CN101246515A (en) Experimental Method of Plane Motion Mechanism of Digital Ship Model Based on CFD Software
CN103303433A (en) Virtual test system and test method for ship performance
CN113092065B (en) Analysis method for calculation of hydrodynamic characteristics of submerged net cage
CN111709151A (en) A method of creating a fast mathematical model of ships
CN111199095B (en) Method for constructing test environment for intelligent ship berthing and unberthing algorithm based on 3D simulation
WO2023216997A1 (en) Wind-induced drift prediction method and apparatus, and device and storage medium
CN110262526B (en) A physics-like numerical method for predicting the 6-DOF self-propulsion maneuvering of underwater robots
CN102521882A (en) Method for obtaining seabed terrain data based on discrete elevation and adaptive mixed weighting
CN114818511A (en) Semi-submersible floating type wind power plant dynamic layout optimization method, system and equipment
CN108820155A (en) The ship gyro based on Fuzzy Immune Controller, which subtracts, under a kind of low speed of a ship or plane shakes control method
CN115994496B (en) Numerical simulation method of high-resolution three-dimensional atmospheric CO2 concentration field in urban parks
CN110937082B (en) Ship overturning risk testing method based on random wind field and sea waves
CN111382537A (en) Analysis method of hydrodynamic characteristics of ocean buoy based on two-phase viscous flow theory
CN102662164B (en) Sea surface current information extraction method based on X-band radar image and particle swarm optimization
Liu et al. A motion tracking approach to position marine floating structures based on measured acceleration and angular velocity
CN118171395B (en) Digital twin system for unmanned sailing boat performance analysis
CN108153976B (en) Simulation method for rolling motion of ship in sea waves and upper computer
CN119167821A (en) A ship maneuvering motion simulation method with real-time human-computer interaction
CN117725670A (en) Method for optimizing overall dimensions of all parts of floating fan platform
CN105427371A (en) Method for keeping graphic object equal-pixel area display in three-dimensional perspective projection scene
CN116310221A (en) Ship swaying and visual wave motion matching method, terminal equipment and medium
CN108253965A (en) A kind of TLP platform stances orientation measuring system

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