Short-peak wave sea wave dynamic response analysis method and system based on finite volume method

Abstract

The invention provides a short-peak wave dynamic response analysis method and system based on a finite volume method, and the method comprises the steps: building a short-peak wave model based on the wind wave characteristic parameters under different sea conditions and the motion state parameters of a ship, and carrying out the calculation under ground coordinates to obtain a first wave surface equation; a first wave surface equation located in a ground coordinate system is converted into a ship body stable coordinate system through a coordinate and frequency conversion method, a second wave surface equation suitable for the surrounding environment of the ship is obtained, then a fluctuation pressure model caused by sea waves is established, the fluctuation pressure model is corrected through a Smith correction term method, and the ship body stable coordinate system is obtained. A multi-wave superposition fluctuation pressure model is obtained through a harmonic superposition analysis method; finally, the generalized disturbing force is constructed based on a finite volume method, the dynamic response of the ship in stormy waves can be predicted more accurately, modeling and simulation can be effectively conducted on the wave force borne by a complex ship body, and the calculation precision is effectively improved.

Description

Technical Field

[0001] The present invention relates to the technical field of ship engineering and ocean engineering, and in particular to a method and system for analyzing the dynamic response of short-peaked waves based on a finite volume method. Background Art

[0002] In modern naval warfare, air power has become an indispensable and important support for naval ships. With the continuous development and growth of my country's navy in recent years, its demand for ship-borne aviation power has also gradually increased. At the same time, the maturity of ship-borne helicopter technology and the continuous development of new vertical take-off and landing aircraft have made non-through-deck ships such as frigates and destroyers play an increasingly important role in the field of ship-borne aviation force delivery.

[0003] However, the ocean environment dominated by wind and waves and ocean weather changes, and the ship environment dominated by ship wake and deck disturbances, will directly affect the flight performance of ship-borne helicopters and other aircraft, and even directly endanger flight safety. Especially in high sea conditions, due to the violent shaking of the ship and the harsh flow field environment, ship-borne helicopters and other aircraft often find it difficult to perform take-off and landing missions normally, which poses a severe challenge to the ship-borne aviation combat capability.

[0004] At present, the existing technologies have certain limitations in dealing with these problems. For example, the existing interference force and interference torque calculation method uses the square ship assumption, that is, the hull is approximated as a rectangular parallelepiped. However, in actual engineering, the hull has a complex shape. Obviously, the above method will cause simulation distortion and inaccurate calculation results.

[0005] Therefore, there is an urgent need for a method that can effectively model and simulate the interference forces and moments acting on a complex hull in order to improve the safety and efficiency of shipborne aviation operations. Summary of the invention

[0006] The present invention aims at the problems that the calculation results of the disturbance force and moment on the complex hull are not accurate enough and the simulation will be distorted. The present invention provides a short-peak wave dynamic response analysis method based on the finite volume method. The short-peak wave model is established based on the characteristic parameters of the wind and waves and the motion state parameters of the ship, and the wave surface equation is obtained by using a specific calculation method and coordinate transformation method, and the multi-wave superposition fluctuation pressure model is obtained by a specific correction method. The generalized disturbance force (dynamic response) is constructed based on the finite volume method, which can more accurately predict the dynamic response of the ship in the wind and waves, and can effectively model and simulate the wave force on the complex hull, effectively improving the calculation accuracy. The present invention also relates to a short-peak wave dynamic response analysis system based on the finite volume method.

[0007] The technical solution of the present invention is as follows:

[0008] A method for analyzing the dynamic response of short-peaked waves based on a finite volume method, characterized in that it comprises the following steps:

[0009] Parameter acquisition and coordinate system establishment steps: obtain wind and wave characteristic parameters and ship motion state parameters under different sea conditions, and establish a ground coordinate system with any point on the sea level as the origin, the ship motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establish a hull stable coordinate system with the ship's center of mass as the origin, the ship's motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establish a hull coordinate system with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency;

[0010] The steps of constructing the short-peaked wave wind and wave model are as follows: in the ground coordinate system, a plane wave equation about the wave height is established based on the wave amplitude, wave number, main propagation direction of wind and wave, propagation direction of wind and wave harmonics, and the time of wind and wave shape transformation, and the short-peaked wave wind and wave model is obtained by recursion based on the plane wave equation;

[0011] The wave surface equation and wave height calculation steps are as follows: the directional spectrum of short-peaked waves is calculated according to the propagation direction of wind and wave harmonics, and a first expression about the energy of long-peaked waves is constructed according to the wave amplitude and seawater density under ground coordinates, an angular frequency interval is established according to the wind and wave angular frequency and the bandwidth set based on the wind and wave angular frequency, and then a second expression about the energy of long-peaked waves within the angular frequency interval is obtained, and a third expression about the energy spectrum of long-peaked waves is recursively obtained from the second expression, and then the long-peaked wave energy spectrum is calculated; a short-peaked wave wind and wave spectrum is calculated according to the directional spectrum of short-peaked waves and the energy spectrum of long-peaked waves, a fourth expression about the energy of short-peaked waves is established based on the short-peaked wave wind and wave spectrum, and the fourth expression of the energy of short-peaked waves is calculated by the fourth expression and the short-peaked wave wind and wave model. A wave surface equation is established; then a first relational expression between the encounter period and the encounter angular frequency is established, and the angle between the ship speed direction and the main propagation direction of the wind and waves is taken as the ship encounter angle, and a second relational expression about the encounter period is established based on the ship encounter angle, the natural period of the wind and waves, the natural wave speed of the wind and waves, and the ship speed, and a third relational expression about the encounter angular frequency is established based on the ship encounter angle, the angular frequency of the wind and waves, the natural wave speed of the wind and waves, and the ship speed; based on the first relational expression, the second relational expression and the third relational expression, and using the coordinate and frequency transformation method, the first wave surface equation located in the ground coordinate system is converted into the hull stable coordinate system, and a second wave surface equation applicable to the surrounding environment of the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation;

[0012] Steps for obtaining the multi-wave superposition fluctuation pressure model: in the hull stable coordinate system, the sea surface depth of a certain underwater point from the free liquid surface is calculated according to the liquid surface depth under still water conditions and the wave height encountered by the ship during navigation, and a fifth expression for the pressure at a certain underwater point is constructed according to the sea surface depth and the atmospheric pressure; a sixth expression for the pressure at a certain underwater point under short-peak wave and sea wave conditions is calculated according to the fifth expression and the short-peak wave wind and wave model; a fluctuation pressure model caused by sea waves is obtained according to the sixth expression, and the fluctuation pressure model is corrected by the Smith correction term method to obtain a corrected fluctuation pressure model; and the corrected fluctuation pressure model is expressed as a superposition of multiple cosine waves by the harmonic superposition analysis method to obtain a multi-wave superposition fluctuation pressure model;

[0013] The steps of constructing the generalized wave disturbance force are as follows: in the stable coordinate system of the hull, the ship is divided into multiple triangular face elements based on the finite volume method, and the vector pointing from the ship's center of mass to the centroid of a certain triangular face element is used as the radius vector of the triangular face element; the wave pressure received by the triangular face element is calculated according to the multi-wave superposition wave pressure model, and the wave force received by the triangular face element is calculated according to the area, normal vector and wave pressure of the triangular face element; the wave moment received by the triangular face element is calculated according to the radius vector of the triangular face element and the calculated wave force; based on the wave forces and wave moments received by all triangular face elements, the wave forces and wave moments of the entire ship are calculated respectively; and the coordinate and frequency transformation methods are used to convert the wave forces and wave moments of the entire ship in the stable coordinate system of the hull into the hull coordinate system, and the component arrays of the overall wave forces of the ship and the component arrays of the overall wave moments of the ship are obtained respectively, and then the generalized wave disturbance force is constructed to realize the dynamic response analysis of short-peak waves.

[0014] Preferably, in the wave surface equation and wave height calculation step, when the bandwidth in the third expression tends to zero, the integral expression of the long-peak wave energy is calculated according to the first expression and the third expression, and the relationship between the wave harmonic amplitude and the spectral density function is calculated according to the integral expression.

[0015] Preferably, in the parameter acquisition and coordinate system establishment steps, the wind and wave characteristic parameters also include wavelength, and the wave number is calculated based on the wavelength.

[0016] Preferably, the coordinate and frequency transformation method includes coordinate transformation based on rotation matrix and transformation of wave frequency to encounter frequency.

[0017] A short-peak wave dynamic response analysis system based on the finite volume method is characterized by comprising a parameter acquisition and coordinate system establishment module, a short-peak wave wind and wave model construction module, a wave surface equation and wave height calculation module, a multi-wave superposition wave pressure model acquisition module and a generalized wave interference force calculation module connected in sequence.

[0018] The parameter acquisition and coordinate system establishment module acquires wind and wave characteristic parameters and ship motion state parameters under different sea conditions, and establishes a ground coordinate system with any point on the sea level as the origin, the ship motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establishes a hull stable coordinate system with the ship's center of mass as the origin, the ship's motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establishes a hull coordinate system with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency;

[0019] The short-peaked wave wind and wave model construction module establishes a plane wave equation about the wave height in the ground coordinate system based on the wave amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, and time of wind and wave shape transformation, and recursively obtains the short-peaked wave wind and wave model according to the plane wave equation;

[0020] The wave surface equation and wave height calculation module calculates the directional spectrum of short-peaked waves according to the propagation direction of wind and wave harmonics, and constructs a first expression for the energy of long-peaked waves according to the wave amplitude and seawater density under ground coordinates, establishes an angular frequency interval according to the wind and wave angular frequency and the frequency band width set based on the wind and wave angular frequency, and then obtains a second expression for the energy of long-peaked waves within the angular frequency interval, and recursively obtains a third expression for the energy spectrum of long-peaked waves from the second expression, and then calculates the energy spectrum of long-peaked waves; calculates the short-peaked wave wind and wave spectrum according to the directional spectrum of short-peaked waves and the long-peaked wave energy spectrum, establishes a fourth expression for the energy of short-peaked waves based on the short-peaked wave wind and wave spectrum, and calculates the short-peaked wave wind and wave model by the fourth expression. The first wave surface equation; then establish a first relationship between the encounter period and the encounter angular frequency, and take the angle between the ship speed direction and the main propagation direction of the wind and waves as the ship encounter angle, and establish a second relationship about the encounter period based on the ship encounter angle, the natural period of wind and waves, the natural wave speed of wind and waves, and the ship speed, and establish a third relationship about the encounter angular frequency based on the ship encounter angle, the angular frequency of wind and waves, the natural wave speed of wind and waves, and the ship speed; based on the first relationship, the second relationship, and the third relationship, and using the coordinate and frequency transformation method, the first wave surface equation located in the ground coordinate system is converted into the hull stable coordinate system, and the second wave surface equation applicable to the surrounding environment of the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation;

[0021] The multi-wave superposition fluctuation pressure model acquisition module calculates the sea surface depth of a certain underwater point from the free liquid surface according to the liquid surface depth under calm water conditions and the wave height encountered by the ship during navigation in a stable coordinate system of the hull, and constructs a fifth expression for the pressure of a certain underwater point according to the sea surface depth and the atmospheric pressure; calculates a sixth expression for the pressure of a certain underwater point under short-peak wave and sea wave conditions according to the fifth expression and the short-peak wave wind and wave model; obtains the fluctuation pressure model caused by the sea waves according to the sixth expression, and corrects the fluctuation pressure model by using the Smith correction term method to obtain the corrected fluctuation pressure model; and expresses the corrected fluctuation pressure model as the superposition of multiple cosine waves by using the harmonic superposition analysis method to obtain the multi-wave superposition fluctuation pressure model;

[0022] The generalized wave disturbance force calculation module divides the ship into multiple triangular face elements based on the finite volume method in a stable hull coordinate system, and uses the vector from the ship's center of mass to the centroid of a certain triangular face element as the radius vector of the triangular face element, calculates the wave pressure on the triangular face element according to the multi-wave superposition wave pressure model, calculates the wave force on the triangular face element according to the area, normal vector and wave pressure of the triangular face element, and calculates the wave moment on the triangular face element according to the radius vector of the triangular face element and the calculated wave force; calculates the wave force and wave moment of the entire ship based on the wave force and wave moment on all triangular face elements; and converts the wave force and wave moment of the entire ship in the stable hull coordinate system into the hull coordinate system by using the coordinate and frequency transformation method, and obtains the component array of the entire ship's wave force and the component array of the entire ship's wave moment, thereby constructing a generalized wave disturbance force to realize the dynamic response analysis of short-peak waves.

[0023] Preferably, in the wave surface equation and wave height calculation module, when the bandwidth in the third expression tends to zero, the integral expression of the long-peak wave energy is calculated according to the first expression and the third expression, and the relationship between the wave harmonic amplitude and the spectral density function is calculated according to the integral expression.

[0024] Preferably, in the parameter acquisition and coordinate system establishment step module, the wind and wave characteristic parameters also include wavelength, and the wave number is calculated based on the wavelength.

[0025] Preferably, the coordinate and frequency transformation method includes coordinate transformation based on rotation matrix and transformation of wave frequency to encounter frequency.

[0026] The beneficial effects of the present invention are:

[0027] The present invention provides a short-peak wave dynamic response analysis method based on the finite volume method. Based on the characteristic parameters of wind and waves under different sea conditions and the motion state parameters of the ship, a short-peak wave wind and wave model including multiple characteristic parameters such as amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics is established, so that the behavior of wind and waves and their influence on the ship can be described more accurately and quickly, thereby improving the accuracy of the model. Based on the established ground coordinate system and hull stable coordinate system, the wind and wave model is converted from the ground coordinate system to the hull coordinate system through the coordinate and frequency transformation method, so that the actual motion state of the ship in the waves can be more accurately reflected. Then, a second wave surface equation suitable for the surrounding environment of the ship is calculated by adopting a specific calculation method and constructed expressions and relations, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation, so that the actual wave conditions can be reflected more accurately, especially for accurate wave height prediction under complex sea conditions, thereby reducing the risk of the ship encountering extreme sea conditions. Thus, the safety of the ship and its crew is protected; the Smith correction term method is then used to correct the wave pressure model obtained by establishing each expression to obtain a corrected wave pressure model to improve the calculation accuracy, which can effectively and accurately evaluate the force of the ship in the wind and waves; and the multi-wave superposition wave pressure model is obtained by the harmonic superposition analysis method. The multi-wave superposition wave pressure model takes into account the pressure change at a certain point underwater, the interaction between wind and waves and the nonlinear effect, effectively improves the accuracy and reliability of the wave pressure calculation, and makes the prediction of the dynamic response of the ship more accurate; finally, based on the finite volume method and using a specific calculation method, a generalized wave interference force is constructed, which can more accurately predict the dynamic response of the ship in the wind and waves, and can effectively model and simulate the wave force on the complex hull, effectively improve the calculation accuracy, and can more accurately evaluate the take-off and landing safety of the ship-borne helicopter under high sea conditions, which is of great significance to ensure the safe take-off and landing of the ship-borne helicopter under high sea conditions. The present invention can effectively improve the take-off and landing capability and safety of carrier-based aircraft under complex sea conditions, thereby enhancing the overall combat effectiveness of carrier-based aviation forces.

[0028] The present invention also relates to a short-peak wave ocean wave dynamic response analysis system based on the finite volume method. The system corresponds to the above-mentioned short-peak wave ocean wave dynamic response analysis method based on the finite volume method, and can be understood as a system that implements the above-mentioned short-peak wave ocean wave dynamic response analysis method based on the finite volume method, including a parameter acquisition and coordinate system establishment module, a long-peak wave ocean wave model construction module, a wave surface equation and wave height calculation module, a multi-wave superposition fluctuation pressure model acquisition module and a generalized wave interference force construction module that are connected in sequence. The modules cooperate with each other to establish a short-peak wave wind and wave model based on the characteristic parameters of wind and waves and the motion state parameters of the ship, and adopt a specific calculation method and coordinate conversion method to obtain the wave surface equation, and obtain the multi-wave superposition fluctuation pressure model through a specific correction method, and then construct a generalized wave interference force based on the finite volume method and a specific calculation method, which can more accurately predict the dynamic response of the ship in wind and waves, and can effectively model and simulate the wave force exerted on the complex hull, effectively improving the calculation accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 It is a flow chart of the ship dynamic stability modeling method based on multiple external force coupling of the present invention.

[0030] Figure 2 It is a schematic diagram of the ground coordinate system of the present invention.

[0031] Figure 3 It is a schematic diagram of the hull stability coordinate system of the present invention.

[0032] Figure 4 It is a schematic diagram of the hull coordinate system of the present invention.

[0033] Figure 5 It is a schematic diagram of the conversion relationship between the hull stable coordinate system and the hull coordinate system of the present invention.

[0034] Figure 6 It is a schematic diagram of the plane wave shape of the present invention.

[0035] Figure 7 It is a schematic diagram of the ship of the present invention sailing in wind and waves.

[0036] Figure 8 It is a schematic diagram of dividing a ship based on the finite volume method according to the present invention. DETAILED DESCRIPTION

[0037] The present invention will be described below in conjunction with the accompanying drawings.

[0038] The present invention relates to a method for analyzing the dynamic response of short-peaked waves based on a finite volume method. The flow chart of the method is as follows: Figure 1 As shown, the following steps are included in sequence:

[0039] Parameter acquisition and coordinate system establishment steps: obtain wind and wave characteristic parameters and ship motion state parameters under different sea conditions, and establish a ground coordinate system with any point on the sea level as the origin, the ship motion direction as the longitudinal axis, the direction of the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establish a hull stable coordinate system with the ship's center of mass as the origin, the ship's motion direction as the longitudinal axis, the direction of the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establish a hull coordinate system with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency.

[0040] Specifically, the ground coordinate system E-ξηζ is Figure 2 As shown in the figure, it is a fixed inertial coordinate system fixed to the earth, which is used to describe the motion state of the waves. Among them, the origin E of the coordinate system can be taken as any point on the sea level, the center of the earth is the vertical axis, that is, the Eζ axis points to the center of the earth, the direction of the ship's motion is the longitudinal axis, that is, the Eξ axis points to the direction of the ship's motion, the direction to the starboard of the ship and parallel to the sea level is the transverse axis, that is, the Eη axis points to the starboard of the ship and parallel to the sea level, the Eξ axis and the Eη axis are located in the sea level and are perpendicular to each other, and the distribution of the three coordinate axes conforms to the right-hand rule. The hull stable coordinate system Ox s y s z s like Figure 3 As shown in Figure 1, it is a coordinate system that moves with the ship and is used to describe the motion of the ship hull in a stable state. The origin of the coordinate system O is the center of mass of the ship, and the direction of motion of the ship is the longitudinal axis, that is, Ox s The axis points to the direction of the ship's motion, and the direction to the starboard side of the ship and parallel to the sea level is the transverse axis, that is, Oy s Pointing to the starboard side of the ship and parallel to the sea level, with the center of the earth as the vertical axis, i.e. Oz s The axis points to the center of the earth, Ox s y s The plane is parallel to the sea level, and the distribution of the three coordinate axes conforms to the right-hand rule. b y b z b ,like Figure 4 As shown in FIG. 1 , it is a motion coordinate system fixed to the hull, used to describe the motion state of the ship hull. The origin of the coordinate system O is the center of mass of the ship, and the direction of the bow of the ship is the longitudinal axis, that is, Ox b The axis coincides with the ship's roll axis and points to the bow, and the direction of the starboard side of the ship is the transverse axis, namely Oy bThe axis coincides with the pitch axis of the ship and points to the starboard side of the ship. The direction of the bottom of the ship is the vertical axis, i.e. Oz b The axis coincides with the bow axis of the ship and points to the bottom of the ship. The selection and definition of the above variables, that is, the symbols of the ship motion and force variables are shown in Table 1.

[0041] Table 1

[0042]

[0043] The three translational position quantities are swell x, swell y and heave (heave) z. The three translational position quantities are swell x, swell y and heave (heave) z. The hull coordinate system can be obtained by three simple rotations of the hull stable coordinate system, such as Figure 5 As shown, the three rotations are: around z s The yaw angle ψ is rotated by the axis, and O-x1y1z is obtained s ; Rotate the pitch angle θ around the y1 axis to get Ox b y1z1; around x b Axis rotation roll angle Get Ox b y b z b .

[0044] The steps of constructing the short-peak wave wind and wave model are as follows: based on the wave amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics and the time of wind and wave shape transformation, a plane wave equation about the wave height is established, and the short-peak wave wind and wave model is obtained by recursive calculation based on the plane wave equation.

[0045] When a ship is sailing under high sea conditions, it will be disturbed by various marine environments such as waves, sea breezes, and ocean currents. Among them, random waves are the most important disturbance factor. Therefore, before analyzing the motion state of the ship under high sea conditions, it is very necessary to understand the distribution of random waves under high sea conditions. Specifically, wave simulation has always been a very complex problem. In real situations, the wind and waves on the sea surface are extremely complex irregular random three-dimensional waves, so it is very difficult to accurately model the waves. However, in the process of studying the problem of ship motion, it is considered that the wind and waves generally propagate in one direction along the wind direction, that is, the wind and waves only propagate along their main propagation direction (that is, the ξ axis of the ground coordinate system). Therefore, it can be assumed that the crests and troughs of the waves are parallel to each other and perpendicular to the wave propagation direction. At this time, assuming that the main propagation direction of the wind and waves is the ξ axis, the wave height (wave surface height) ζ at a fixed point on the sea surface is expressed as a binary function of ξ and time t ζ = ζ(ξ,t). According to the Longuet-Higgins model, the long-peak wave model of irregular wind waves can be regarded as the superposition of countless waves with different amplitudes (amplitudes), wavelengths and initial phases in the vertical plane of the main propagation direction of the wind waves. According to the following fluid mechanics assumptions:

[0046] a) The fluid is incompressible;

[0047] b) The fluid has only potential velocity;

[0048] c) The wind and waves are slight.

[0049] The plane wave propagating along the ξ axis can be expressed as:

[0050] ζ=ζ a cos(kξ-ωt) (1)

[0051] Where ξ is the wave height (also known as the distance the wave surface deviates from the sea level, or the height of the wave surface relative to the mean sea level); ζ a is the amplitude; k is the wave number; λ is the wavelength, and the wave number is related to the wavelength: k = 2π / λ; ω is the angular frequency of wind waves. t is the time when the angular frequency of wind waves and the shape of wind waves change, describing the change of wind waves over time. When t = 0, its plane wave shape is as follows Figure 6 shown.

[0052] According to the plane wave formula, the long crest wave model (also called the Longuet-Higgins model of long crest waves) can be derived, as shown below:

[0053]

[0054] Among them, k i ,ω i , ai are the wave number, wind wave angular frequency and amplitude of the i-th harmonic respectively; ε i is the random initial phase of the i-th harmonic between 0 and 2π.

[0055] However, the above long-peak wave model still has problems such as being unable to accurately capture the propagation characteristics of waves and slow calculation speed because it only considers the main propagation direction of wind waves as the ξ axis. Therefore, the short-peak wave model of wind waves is introduced. This model considers that the propagation direction of different harmonics of wind waves will propagate along different propagation angles μ. It is defined that μ=0 when propagating along the ξ axis, and the angle μ is positive when rotating positively around the ζ axis according to the right-hand rule. Therefore, in the short-peak wave model of wind waves, the wave height ζ at a fixed point on the sea surface should be expressed as a three-variable function of ξ, η and time t ζ=ζ(ξ,η,t). The wave equation of the plane propagating at the propagation angle μ is re-described in the three-dimensional ground coordinate system E-ξηζ as follows:

[0056] ζ=ζ a cos[k(ξcosμ j -ηsinμ j )-ωt] (3)

[0057] Among them, a is the amplitude, k is the wave number, ζ is the wave height, μ represents the propagation angle, the propagation direction of different harmonics of wind and waves is along different propagation angles μ, the main propagation direction of wind and waves along the ξ axis, the main propagation direction of wind and waves along the η axis, ω is the angular frequency of wind and waves.

[0058] The short-peaked wave model (also called the Longuet-Higgins model of short-peaked wave) derived from the plane wave equation is as follows:

[0059]

[0060] The short-peaked wave model takes into account the main propagation direction of wind and waves along the ξ-axis and η-axis, and introduces the propagation direction of wind and wave harmonics as propagation along different propagation angles μ, so that the short-peaked wave model can more accurately capture the propagation characteristics of waves (including the propagation direction and energy distribution of wind and waves). And because the short-peaked wave model contains more degrees of freedom, it can more carefully characterize the irregularity of waves, thereby improving the accuracy of prediction, which is crucial for studying wave dynamics, marine engineering design and safety assessment. At the same time, the short-peaked wave model is also applicable to various complex marine environments, including wave conditions in different sea areas, seasons and weather conditions, and can better simulate the wind and wave behavior in the real world, providing a more reliable foundation for scientific research and engineering practice in related fields.

[0061] The wave surface equation and wave height calculation steps are as follows: the directional spectrum of short-peaked waves is calculated according to the propagation direction of wind and wave harmonics, and a first expression about the energy of long-peaked waves is constructed according to the wave amplitude and seawater density under ground coordinates, an angular frequency interval is established according to the wind and wave angular frequency and the bandwidth set based on the wind and wave angular frequency, and then a second expression about the energy of long-peaked waves within the angular frequency interval is obtained, and a third expression about the energy spectrum of long-peaked waves is recursively obtained from the second expression, and then the long-peaked wave energy spectrum is calculated; a short-peaked wave wind and wave spectrum is calculated according to the directional spectrum of short-peaked waves and the energy spectrum of long-peaked waves, a fourth expression about the energy of short-peaked waves is established based on the short-peaked wave wind and wave spectrum, and the fourth expression of the energy of short-peaked waves is calculated by the fourth expression and the short-peaked wave wind and wave model. A wave surface equation is established; then a first relationship between the encounter period and the encounter angular frequency is established, and the angle between the ship speed direction and the main propagation direction of the wind and waves is taken as the ship encounter angle, and a second relationship about the encounter period is established based on the ship encounter angle, the natural period of wind and waves, the natural wave speed of wind and waves and the ship speed, and a third relationship about the encounter angular frequency is established based on the ship encounter angle, the angular frequency of wind and waves, the natural wave speed of wind and waves and the ship speed; based on the first relationship, the second relationship and the third relationship, the coordinate and frequency transformation method is used to transform the first wave surface equation located in the ground coordinate system into the stable coordinate system of the hull, and the second wave surface equation applicable to the environment surrounding the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation.

[0062] Specifically, the directional spectrum of short-peaked waves is first calculated according to the propagation direction of wind wave harmonics, and the calculation is performed according to the following formula:

[0063]

[0064] The 12th ITTC suggested that n in M(μ) be temporarily set to 1, then:

[0065]

[0066] Then, according to the knowledge of fluid mechanics, long-crest waves have energy. Therefore, in the ground coordinates, the first expression of the energy E of long-crest waves is constructed according to the wave amplitude and seawater density, which is expressed as follows:

[0067]

[0068] Where i is the harmonic order, ζ ai is the amplitude of the ith harmonic, ρ is the seawater density, and g is the gravitational acceleration.

[0069] According to the angular frequency ω of the wind and waves and the bandwidth dω set based on the angular frequency of the wind and waves, an angular frequency interval (ω, ω+dω) is established, and then a second expression for the energy of long-peak waves in the angular frequency interval is obtained, which is expressed as follows:

[0070]

[0071] The third expression about the long-peak wave energy spectrum is derived from the second expression, that is, equation (8), which is defined as: ζ (ω) is shown in the following formula:

[0072]

[0073] According to formula 9, the long peak wave energy spectrum S can be calculated ζ (ω).

[0074] Preferably, when the bandwidth in the third expression tends to zero, that is, when dω→0, the integral expression of the long crest wave energy is calculated according to the first expression and the third expression, as shown below:

[0075]

[0076] According to the integral expression of long crest wave energy, the wave harmonic amplitude ζ can be obtained: ai With the spectral density function S ζ (ω) is related to:

[0077]

[0078] The short-peak wave wind wave spectrum is in the long-peak wave energy spectrum S ζ The directional spectrum of short-peaked waves M(μ) is added to the basis of (ω). Therefore, the short-peaked wave wind wave spectrum is calculated according to the directional spectrum of short-peaked waves and the energy spectrum of long-peaked waves, which is expressed as follows:

[0079] S ζ (ω,μ)=S ζ (ω)M(μ) (10)

[0080] The fourth expression for the energy of short-peaked waves is established based on the short-peaked wave spectrum, which is expressed as follows:

[0081]

[0082] Discretizing the above equation 11 and substituting it into the short crest wave wind wave model can obtain the first wave surface equation of the short crest wave wind wave, which is expressed as follows:

[0083]

[0084] In summary, the state of the waves at a certain moment can be obtained (the state includes wave height, energy distribution, waveform and phase distribution).

[0085] When a ship is sailing on the sea, due to its certain speed, the wave frequency encountered by the ship is not the same as the wave frequency at a fixed point on the sea. When the ship is sailing against the waves, the wave frequency it encounters is higher than the wave frequency at a fixed point on the sea; conversely, when the ship is sailing with the waves, the wave frequency it encounters is lower than the wave frequency at a fixed point on the sea. Therefore, when studying the problem of ships sailing in wind and waves, we first define the wave period encountered by the ship as the encounter period T e , the corresponding angular frequency is called the encounter angular frequency ω e , which is related to the encounter period, establish the first relationship between the encounter period and the encounter angular frequency, that is, In addition, if Figure 7 As shown, the encounter period of the ship is T e and the encounter angular frequency ω e The natural period of the sea wave T, the angular frequency of the wind and wave ω (also called the natural angular frequency ω), the ship speed V, and the ship encounter angle μ e (the angle between the ship speed V and the wind and wave propagation direction) and the inherent wave speed C of the wind and waves, and the relationship can be expressed as:

[0086]

[0087]

[0088] According to the above equations (13), (14) and the encounter period T e and the encounter angular frequency ω e The geometric relationship between the wave surface and the wave frequency is obtained by converting the first wave surface equation in the ground coordinate system, that is, equation (12), into the stable coordinate system of the hull through the coordinate and frequency transformation method (coordinate transformation based on the rotation matrix and the transformation from wave frequency to encounter frequency). The second wave surface equation suitable for the short-peaked waves around the ship is obtained, which is expressed as follows:

[0089]

[0090] According to the second wave surface equation of short-peaked wind waves, that is, according to formula (15), the wave height ζ encountered by the ship during navigation can be calculated: e , which indicates the shape of the wave surface encountered by the ship during navigation.

[0091] Steps for obtaining a multi-wave superposition fluctuation pressure model: in a stable coordinate system of the hull, calculate the sea surface depth of a certain underwater point from the free liquid surface according to the liquid surface depth under still water conditions and the wave height encountered by the ship during navigation, and construct a fifth expression for the pressure at a certain underwater point according to the sea surface depth and the atmospheric pressure; calculate a sixth expression for the pressure at a certain underwater point under short-peak wave and sea wave conditions according to the fifth expression and the short-peak wave wind and wave model; obtain a fluctuation pressure model caused by sea waves according to the sixth expression, and use the Smith correction term method to correct the fluctuation pressure model to obtain a corrected fluctuation pressure model; and express the corrected fluctuation pressure model as a superposition of multiple cosine waves through the harmonic superposition analysis method to obtain a multi-wave superposition fluctuation pressure model.

[0092] Specifically, this step can also be understood as a step for studying the distribution of the seawater pressure field under the sea surface. During the research process, it is generally believed that the interference force and torque of the waves on the hull are caused by the fluctuation of the pressure field distribution under the liquid surface. Therefore, before analyzing the wind and wave forces and torques on the ship, the distribution of the seawater pressure field under the sea surface is first studied. From the basic water pressure formula, it can be seen that the fifth expression, that is, the pressure P at a certain point underwater, can be expressed as:

[0093] P=P0+ρgh (16)

[0094] Where P0 is the atmospheric pressure; ρ is the density of seawater; g is the acceleration of gravity; and h is the depth of a certain point underwater (also called a research point) from the free surface of the liquid. Under wave conditions, the sea surface depth h can be expressed as the liquid surface depth z under still water conditions and the wave height ζ encountered by the ship during navigation. e The sum is expressed as follows:

[0095] h=z+ζ e (17)

[0096] According to the fifth expression and the short-peak wave wind wave model, the sixth expression of the pressure at a certain point underwater under the condition of short-peak wave sea waves is calculated. That is, the sixth expression of the pressure at a certain point underwater under the condition of long-peak wave sea waves is obtained by substituting equation (4) into equation (16), as shown in the following equation:

[0097] P=P0+ρgz+ρgζ a cos(kξ-ωt) (18)

[0098] Based on the sixth expression, the pressure fluctuation model caused by the waves is defined as ΔP, then:

[0099] ΔP=ρgζ a cos(kξ-ωt) (19)

[0100] According to the Smith effect, the amplitude of the wave surface fluctuation is affected by the depth, so the Smith correction term is used to correct the fluctuation pressure model to obtain the corrected fluctuation pressure model, as shown in the following formula:

[0101]

[0102]

[0103] ΔP=e -kz ρgζ a cos(kξ-ωt) (22)

[0104] Expanded in the hull stable coordinate system, we have:

[0105] ΔP=e -kz ρgζ a cos(kx s cosμ e +ky s sinμ e -ωe t ) (twenty three)

[0106] The modified fluctuating pressure model is expressed as the superposition of multiple cosine waves through the harmonic superposition analysis method, and the multi-wave superposition fluctuating pressure model is obtained, that is, formula (23) is expressed as the cosine wave superposition form, as shown in the following formula:

[0107]

[0108] The generalized wave disturbance force construction step can also be understood as the disturbance force and moment analysis step of the ship. By dividing the complex hull into surface meshes (triangular face elements), and by solving and superimposing the wave forces on each triangular face element, the wave forces acting on the complex hull can be effectively modeled and simulated. In the stable coordinate system of the hull, the ship is divided into multiple triangular face elements based on the finite volume method, and the vector pointing from the ship's center of mass to the centroid of a certain triangular face element is taken as the radius vector of the triangular face element. The wave pressure exerted on the triangular face element is calculated according to the multi-wave superposition wave pressure model, and the wave force exerted on the triangular face element is calculated according to the area, normal vector and wave pressure of the triangular face element. The wave moment exerted on the triangular face element is calculated according to the radius vector of the triangular face element and the calculated wave force. Based on the wave forces and wave moments exerted on all triangular face elements, the wave force and wave moment of the entire ship are calculated respectively. The coordinate and frequency transformation method is used to transform the wave force and wave moment of the entire ship in the stable coordinate system of the hull into the hull coordinate system, and the component arrays of the overall wave force of the ship and the component arrays of the overall wave moment of the ship are obtained respectively, thereby constructing the generalized wave interference force.

[0109] Specifically, Figure 8 As shown in the figure, in the stable coordinate system of the hull, the ship is divided into N triangular face elements based on the finite volume method, and the vector pointing from the ship's center of mass O to the centroid of the i-th triangular face element is taken as the radius vector r of the triangular face element. i , then the wave pressure and wave force can be calculated at the centroid of each face element through the multi-wave superposition wave pressure model. Among them, the wave pressure ΔP of the triangular face element is calculated according to the multi-wave superposition wave pressure model. i , and calculate the wave force on the i-th triangular surface element according to the area, normal vector and wave pressure of the triangular surface element, and calculate according to the following formula:

[0110] ΔF i =ΔP i n i S i (25)

[0111] In the above formula, the subscript i represents the i-th triangle face element; n i is the normal vector of the i-th triangle face; S i is the area of ​​the i-th face element; ΔP i is the fluctuating pressure ΔP on the i-th triangle element i , which is x s ,y s ,z,t constitutes the quaternion function ΔP(x s ,y s ,z,t), where x s ,y s ,z three parameters can be directly obtained from the ship's surface mesh file.

[0112] The vector pointing from the ship's mass center O to the centroid of the i-th triangle element is defined as the radius r of the triangle element. i , then according to the radius vector of the i-th triangle element and the calculated wave force ΔF i The wave moment on the i-th triangular surface element can be calculated according to the following formula:

[0113] ΔM i =r i ×ΔF i (26)

[0114] Then, based on the wave forces and wave moments on all triangular face elements, the wave force F and wave moment M of the entire ship can be calculated respectively. That is, the wave force (moment) of the entire ship can be expressed as the sum of the wave forces (moments) on each triangular face element as shown below:

[0115]

[0116] The coordinate and frequency transformation method is used to transform the wave force F and wave moment M of the whole ship in the stable coordinate system of the hull into the hull coordinate system, and the component arrays (F) of the whole wave force in the coordinate system of the ship body are obtained respectively. b , and the component array of the overall wave moment in the ship body coordinate system (M) b , and then through (F) b and (M) b The generalized wave disturbance force w is constructed and expressed as follows:

[0117]

[0118] Among them, (F) zb 9F) b The third term of is the force in the z-axis direction in the ship's coordinate system, and (M) xb, (M0 xb For (M) b The first and second items in the equation are the moments in the x- and y-axis directions in the ship's coordinate system. In the subsequent six-degree-of-freedom dynamics study, the generalized wave disturbance force w is recorded as follows:

[0119]

[0120] The present invention also relates to a short-peak wave ocean wave dynamic response analysis system based on the finite volume method. The system corresponds to the short-peak wave ocean wave dynamic response analysis method based on the finite volume method, and can be understood as a system for implementing the above method. The system includes a parameter acquisition and coordinate system establishment module, a short-peak wave wind and wave model construction module, a wave surface equation and wave height calculation module, a multi-wave superposition fluctuation pressure model acquisition module, and a generalized wave interference force calculation module, which are connected in sequence. Specifically,

[0121] The parameter acquisition and coordinate system establishment module acquires wind and wave characteristic parameters and ship motion state parameters under different sea conditions, and establishes a ground coordinate system with any point on the sea level as the origin, the ship motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establishes a hull stable coordinate system with the ship's center of mass as the origin, the ship's motion direction as the longitudinal axis, the direction of the starboard side of the ship parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establishes a hull coordinate system with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency;

[0122] The short-peaked wave wind and wave model construction module establishes a plane wave equation about the wave height in the ground coordinate system based on the wave amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, and time of wind and wave shape transformation, and recursively obtains the short-peaked wave wind and wave model according to the plane wave equation;

[0123] The wave surface equation and wave height calculation module calculates the directional spectrum of short-peaked waves according to the propagation direction of wind and wave harmonics, and constructs a first expression for the energy of long-peaked waves according to the wave amplitude and seawater density under ground coordinates, establishes an angular frequency interval according to the wind and wave angular frequency and the frequency band width set based on the wind and wave angular frequency, and then obtains a second expression for the energy of long-peaked waves within the angular frequency interval, and recursively obtains a third expression for the energy spectrum of long-peaked waves from the second expression, and then calculates the energy spectrum of long-peaked waves; calculates the short-peaked wave wind and wave spectrum according to the directional spectrum of short-peaked waves and the long-peaked wave energy spectrum, establishes a fourth expression for the energy of short-peaked waves based on the short-peaked wave wind and wave spectrum, and calculates the short-peaked wave wind and wave model by the fourth expression. The first wave surface equation; then establish a first relationship between the encounter period and the encounter angular frequency, and take the angle between the ship speed direction and the main propagation direction of the wind and waves as the ship encounter angle, and establish a second relationship about the encounter period based on the ship encounter angle, the natural period of wind and waves, the natural wave speed of wind and waves, and the ship speed, and establish a third relationship about the encounter angular frequency based on the ship encounter angle, the angular frequency of wind and waves, the natural wave speed of wind and waves, and the ship speed; based on the first relationship, the second relationship, and the third relationship, and using the coordinate and frequency transformation method, the first wave surface equation located in the ground coordinate system is converted into the hull stable coordinate system, and the second wave surface equation applicable to the surrounding environment of the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation;

[0124] The multi-wave superposition fluctuation pressure model acquisition module calculates the sea surface depth of a certain underwater point from the free liquid surface according to the liquid surface depth under calm water conditions and the wave height encountered by the ship during navigation in a stable coordinate system of the hull, and constructs a fifth expression for the pressure of a certain underwater point according to the sea surface depth and the atmospheric pressure; calculates a sixth expression for the pressure of a certain underwater point under short-peak wave and sea wave conditions according to the fifth expression and the short-peak wave wind and wave model; obtains the fluctuation pressure model caused by the sea waves according to the sixth expression, and corrects the fluctuation pressure model by using the Smith correction term method to obtain the corrected fluctuation pressure model; and expresses the corrected fluctuation pressure model as the superposition of multiple cosine waves by using the harmonic superposition analysis method to obtain the multi-wave superposition fluctuation pressure model;

[0125] The generalized wave disturbance force calculation module divides the ship into multiple triangular face elements based on the finite volume method in a stable hull coordinate system, and uses the vector from the ship's center of mass to the centroid of a certain triangular face element as the radius vector of the triangular face element, calculates the wave pressure on the triangular face element according to a multi-wave superposition wave pressure model, calculates the wave force on the triangular face element according to the area, normal vector and wave pressure of the triangular face element, and calculates the wave moment on the triangular face element according to the radius vector of the triangular face element and the calculated wave force; calculates the wave force and wave moment of the entire ship based on the wave force and wave moment on all triangular face elements; and converts the wave force and wave moment of the entire ship in the stable hull coordinate system into the hull coordinate system by using a coordinate and frequency transformation method, and obtains a component array of the entire ship's wave force and a component array of the entire ship's wave moment, thereby constructing a generalized wave disturbance force.

[0126] Preferably, in the wave surface equation and wave height calculation module, when the bandwidth in the third expression tends to zero, the integral expression of the long-peak wave energy is calculated according to the first expression and the third expression, and the relationship between the wave harmonic amplitude and the spectral density function is calculated according to the integral expression.

[0127] Preferably, in the parameter acquisition and coordinate system establishment step module, the wind and wave characteristic parameters also include wavelength, and the wave number is calculated based on the wavelength.

[0128] Preferably, the coordinate and frequency transformation method includes coordinate transformation based on rotation matrix and transformation of wave frequency to encounter frequency.

[0129] The present invention provides an objective and scientific short-peak wave dynamic response analysis method and system based on the finite volume method. A short-peak wave model is established based on the characteristic parameters of wind and waves and the motion state parameters of the ship, and a wave surface equation is obtained by using a specific calculation method and coordinate transformation method. A multi-wave superposition fluctuation pressure model is obtained by a specific correction method, and generalized interference force and torque (dynamic response) are constructed based on the finite volume method. The dynamic response of the ship in wind and waves can be predicted more accurately, and the wave force acting on the complex hull can be effectively modeled and simulated, effectively improving the calculation accuracy.

[0130] It should be noted that the above-described specific implementations can enable those skilled in the art to more fully understand the invention, but do not limit the invention in any way. Therefore, although this specification has described the invention in detail with reference to the drawings and embodiments, those skilled in the art should understand that the invention can still be modified or replaced by equivalents. In short, all technical solutions and improvements that do not deviate from the spirit and scope of the invention should be included in the protection scope of the patent for the invention.

Claims

1. A method for analyzing the dynamic response of short-peaked ocean waves based on the finite volume method, characterized in that: The following steps are involved: Parameter acquisition and coordinate system establishment steps: obtain wind and wave characteristic parameters and ship motion state parameters under different sea conditions, and establish a ground coordinate system with any point on the sea level as the origin, the ship motion direction as the longitudinal axis, the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establish a hull stable coordinate system with the ship's center of mass as the origin, the ship motion direction as the longitudinal axis, the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; A hull coordinate system is established with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency; The steps of constructing the short-peaked wave wind and wave model are as follows: in the ground coordinate system, a plane wave equation about the wave height is established based on the wave amplitude, wave number, main propagation direction of wind and wave, propagation direction of wind and wave harmonics, and the time of wind and wave shape transformation, and the short-peaked wave wind and wave model is obtained by recursion based on the plane wave equation; The wave surface equation and wave height calculation steps are as follows: the directional spectrum of short-peaked waves is calculated according to the propagation direction of wind and wave harmonics, and a first expression about the energy of long-peaked waves is constructed according to the wave amplitude and seawater density under ground coordinates, an angular frequency interval is established according to the wind and wave angular frequency and the bandwidth set based on the wind and wave angular frequency, and then a second expression about the energy of long-peaked waves within the angular frequency interval is obtained, and a third expression about the energy spectrum of long-peaked waves is recursively obtained from the second expression, and then the long-peaked wave energy spectrum is calculated; a short-peaked wave wind and wave spectrum is calculated according to the directional spectrum of short-peaked waves and the energy spectrum of long-peaked waves, a fourth expression about the energy of short-peaked waves is established based on the short-peaked wave wind and wave spectrum, and the fourth expression of the energy of short-peaked waves is calculated by the fourth expression and the short-peaked wave wind and wave model. A wave surface equation is established; then a first relational expression between the encounter period and the encounter angular frequency is established, and the angle between the ship speed direction and the main propagation direction of the wind and waves is taken as the ship encounter angle, and a second relational expression about the encounter period is established based on the ship encounter angle, the natural period of the wind and waves, the natural wave speed of the wind and waves, and the ship speed, and a third relational expression about the encounter angular frequency is established based on the ship encounter angle, the angular frequency of the wind and waves, the natural wave speed of the wind and waves, and the ship speed; based on the first relational expression, the second relational expression and the third relational expression, and using the coordinate and frequency transformation method, the first wave surface equation located in the ground coordinate system is converted into the hull stable coordinate system, and a second wave surface equation applicable to the surrounding environment of the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation; The steps of obtaining the multi-wave superposition fluctuation pressure model are as follows: in the hull stable coordinate system, the sea surface depth of a certain underwater point from the free liquid surface is calculated according to the liquid surface depth under calm water conditions and the wave height encountered by the ship during navigation, and the fifth expression about the pressure of a certain underwater point is constructed according to the sea surface depth and the atmospheric pressure; the sixth expression of the pressure of a certain underwater point under the short-peak wave and sea wave conditions is calculated according to the fifth expression and the short-peak wave wind and wave model; The fluctuating pressure model caused by the ocean waves is obtained according to the sixth expression, and the fluctuating pressure model is corrected by using the Smith correction term method to obtain the corrected fluctuating pressure model; and the corrected fluctuating pressure model is expressed as the superposition of multiple cosine waves by using the harmonic superposition analysis method to obtain the multi-wave superposition fluctuating pressure model; The steps of constructing the generalized wave disturbance force are as follows: in the stable coordinate system of the hull, the ship is divided into multiple triangular face elements based on the finite volume method, and the vector pointing from the ship's center of mass to the centroid of a certain triangular face element is used as the radius vector of the triangular face element; the wave pressure received by the triangular face element is calculated according to the multi-wave superposition wave pressure model, and the wave force received by the triangular face element is calculated according to the area, normal vector and wave pressure of the triangular face element; the wave moment received by the triangular face element is calculated according to the radius vector of the triangular face element and the calculated wave force; based on the wave forces and wave moments received by all triangular face elements, the wave forces and wave moments of the entire ship are calculated respectively; and the coordinate and frequency transformation methods are used to convert the wave forces and wave moments of the entire ship in the stable coordinate system of the hull into the hull coordinate system, and the component arrays of the overall wave forces of the ship and the component arrays of the overall wave moments of the ship are obtained respectively, and then the generalized wave disturbance force is constructed to realize the dynamic response analysis of short-peak waves.

2. The short-peak wave dynamic response analysis method based on the finite volume method according to claim 1 is characterized in that: In the wave surface equation and wave height calculation steps, when the bandwidth in the third expression tends to zero, the integral expression of the long-peak wave energy is calculated according to the first expression and the third expression, and the relationship between the wave harmonic amplitude and the spectral density function is calculated according to the integral expression.

3. The short crest wave dynamic response analysis method based on the finite volume method according to claim 1 is characterized in that: In the parameter acquisition and coordinate system establishment steps, the wind and wave characteristic parameters also include wavelength, and the wave number is calculated based on the wavelength.

4. The method for analyzing the dynamic response of short-peaked waves based on the finite volume method according to claim 1 is characterized in that: The coordinate and frequency transformation method includes coordinate transformation based on rotation matrix and transformation from wave frequency to encounter frequency.

5. A short-peak wave dynamic response analysis system based on the finite volume method, characterized in that: It includes a parameter acquisition and coordinate system establishment module, a short-peak wave wind and wave model construction module, a wave surface equation and wave height calculation module, a multi-wave superposition wave pressure model acquisition module and a generalized wave interference force calculation module. The parameter acquisition and coordinate system establishment module acquires the wind and wave characteristic parameters and the motion state parameters of the ship under different sea conditions, and establishes a ground coordinate system with any point on the sea level as the origin, the ship's motion direction as the longitudinal axis, the direction to the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; establishes a hull stable coordinate system with the ship's center of mass as the origin, the ship's motion direction as the longitudinal axis, the direction to the starboard side of the ship and parallel to the sea level as the transverse axis, and the center of the earth as the vertical axis; A hull coordinate system is established with the ship's center of mass as the origin, the direction of the ship's bow as the longitudinal axis, the direction of the starboard side of the ship as the transverse axis, and the direction of the bottom of the ship as the vertical axis; the wind and wave characteristic parameters include amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, wind and wave angular frequency, time of wind and wave shape transformation, wind and wave natural period, and wind and wave natural wave speed; the ship's motion state parameters include ship speed, encounter period, and encounter angular frequency; The short-peaked wave wind and wave model construction module establishes a plane wave equation about the wave height in the ground coordinate system based on the wave amplitude, wave number, main propagation direction of wind and waves, propagation direction of wind and wave harmonics, and time of wind and wave shape transformation, and recursively obtains the short-peaked wave wind and wave model according to the plane wave equation; The wave surface equation and wave height calculation module calculates the directional spectrum of short-peaked waves according to the propagation direction of wind and wave harmonics, and constructs a first expression for the energy of long-peaked waves according to the wave amplitude and seawater density under ground coordinates, establishes an angular frequency interval according to the wind and wave angular frequency and the frequency band width set based on the wind and wave angular frequency, and then obtains a second expression for the energy of long-peaked waves within the angular frequency interval, and recursively obtains a third expression for the energy spectrum of long-peaked waves from the second expression, and then calculates the energy spectrum of long-peaked waves; calculates the short-peaked wave wind and wave spectrum according to the directional spectrum of short-peaked waves and the long-peaked wave energy spectrum, establishes a fourth expression for the energy of short-peaked waves based on the short-peaked wave wind and wave spectrum, and calculates the short-peaked wave wind and wave model by the fourth expression. The first wave surface equation; then establish a first relationship between the encounter period and the encounter angular frequency, and take the angle between the ship speed direction and the main propagation direction of the wind and waves as the ship encounter angle, and establish a second relationship about the encounter period based on the ship encounter angle, the natural period of wind and waves, the natural wave speed of wind and waves, and the ship speed, and establish a third relationship about the encounter angular frequency based on the ship encounter angle, the angular frequency of wind and waves, the natural wave speed of wind and waves, and the ship speed; based on the first relationship, the second relationship, and the third relationship, and using the coordinate and frequency transformation method, the first wave surface equation located in the ground coordinate system is converted into the hull stable coordinate system, and the second wave surface equation applicable to the surrounding environment of the ship is obtained, and the wave height encountered by the ship during navigation is calculated according to the second wave surface equation; The multi-wave superposition fluctuation pressure model acquisition module calculates the sea surface depth of a certain underwater point from the free liquid surface according to the liquid surface depth under calm water conditions and the wave height encountered by the ship during navigation in the hull stable coordinate system, and constructs a fifth expression for the pressure of a certain underwater point according to the sea surface depth and the atmospheric pressure; calculates a sixth expression for the pressure of a certain underwater point under short-peak wave conditions according to the fifth expression and the short-peak wave wind and wave model; The fluctuating pressure model caused by the ocean waves is obtained according to the sixth expression, and the fluctuating pressure model is corrected by using the Smith correction term method to obtain the corrected fluctuating pressure model; and the corrected fluctuating pressure model is expressed as the superposition of multiple cosine waves by using the harmonic superposition analysis method to obtain the multi-wave superposition fluctuating pressure model; The generalized wave disturbance force calculation module divides the ship into multiple triangular face elements based on the finite volume method in a stable hull coordinate system, and uses the vector from the ship's center of mass to the centroid of a certain triangular face element as the radius vector of the triangular face element, calculates the wave pressure on the triangular face element according to the multi-wave superposition wave pressure model, calculates the wave force on the triangular face element according to the area, normal vector and wave pressure of the triangular face element, and calculates the wave moment on the triangular face element according to the radius vector of the triangular face element and the calculated wave force; calculates the wave force and wave moment of the entire ship based on the wave force and wave moment on all triangular face elements; and converts the wave force and wave moment of the entire ship in the stable hull coordinate system into the hull coordinate system by using the coordinate and frequency transformation method, and obtains the component array of the entire ship's wave force and the component array of the entire ship's wave moment, thereby constructing a generalized wave disturbance force to realize the dynamic response analysis of short-peak waves.

6. The short-peak wave dynamic response analysis system based on the finite volume method according to claim 5 is characterized in that: In the wave surface equation and wave height calculation module, when the bandwidth in the third expression tends to zero, the integral expression of the long-peak wave energy is calculated according to the first expression and the third expression, and the relationship between the wave harmonic amplitude and the spectral density function is calculated according to the integral expression.

7. The short-peak wave dynamic response analysis system based on the finite volume method according to claim 5 is characterized in that: In the parameter acquisition and coordinate system establishment step module, the wind and wave characteristic parameters also include wavelength, and the wave number is calculated based on the wavelength.

8. The short crest wave dynamic response analysis system based on finite volume method according to claim 5 is characterized in that: The coordinate and frequency transformation method includes coordinate transformation based on rotation matrix and transformation from wave frequency to encounter frequency.

Images