Laplace domain calculation method for transient hydroelastic response of very large floating bodies

Abstract

The present application belongs to the field of marine engineering technology, and relates to a Laplace domain calculation method for the transient hydroelastic response of a super-large floating body, including the following steps: discretizing the super-large floating structure to obtain the overall stiffness matrix of the discretized super-large floating structure; calculating the discrete frequency response function of the floating hydrodynamic system based on the overall stiffness matrix of the discretized super-large floating structure; Laplace transforming the frequency response function of the floating hydrodynamic system to establish the pole-residue model matrix of the super-large floating structure system transfer function; obtaining the super-large floating time-varying hydroelastic response based on the pole-residue model of the super-large floating structure system transfer function. Through this application, the technical defects of the traditional frequency domain analysis method that can only solve the steady-state response and the time domain method numerical integration is very time-consuming are overcome, and a new Laplace domain method with both high precision and high efficiency is provided for the preliminary engineering analysis of super-large floating bodies.

Description

Technical Field

[0001] The present application relates to the field of marine engineering technology, and in particular to a Laplace domain calculation method for transient hydroelastic response of a super-large floating body. Background Art

[0002] In the late 1980s, super-large floating structures such as offshore floating airports, offshore military bases, offshore port terminals, and large floating bridges became a research hotspot in the marine engineering community. In recent decades, super-large floating structures have received great attention in marine utilization, and various concepts have been successfully put into practice. Compared with conventional marine floating structures, super-large floating bodies have the same order of magnitude of elastic deformation and rigid body displacement due to their huge size and relatively low structural stiffness. Therefore, the hydroelastic theory must be used to solve their displacement and force.

[0003] At present, the calculation of the hydroelastic response of very large floating bodies mainly adopts frequency domain or time domain methods. The traditional frequency domain method can only produce steady-state response and lacks the ability to capture transient behavior. In addition, the efficiency of the traditional time domain method is significantly lower, mainly because time domain methods, such as the direct solution method of the time domain Green's function and the Cummins method, can accurately predict the transient dynamic response of offshore floating structures through step-by-step integration, but the numerical method of the time domain Green's function has a huge amount of calculation and is very time-consuming to solve. The Cummins method is based on the frequency domain solution, which is more efficient than the direct solution method of the time domain Green's function, but it still needs to solve the radiation convolution term and the transient response of the floating body through step-by-step numerical integration, especially for long-term simulation with small time steps, and its calculation is still relatively time-consuming.

[0004] The current traditional frequency domain and traditional time domain methods cannot take into account both computational accuracy and efficiency, and it is difficult to effectively meet the engineering design and analysis needs of efficient and accurate prediction of time-varying hydroelastic responses of ultra-large floating bodies under multiple working conditions, for a long time, and with multiple steps. Therefore, it is urgent to develop an accurate and efficient hydroelastic analysis method for ultra-large floating structures to serve the engineering design and analysis of ultra-large floating structures, especially in the transient response analysis of floating structures with larger scale and more units. Summary of the invention

[0005] The embodiment of the present application provides a Laplace domain calculation method for the transient hydroelastic response of a very large floating body, so as to at least solve the problems of long calculation time and low efficiency in the calculation of the hydroelastic response of a very large floating body in the related art.

[0006] In a first aspect, an embodiment of the present application provides a Laplace domain calculation method for transient hydroelastic response of a very large floating body, comprising the following steps:

[0007] S1: Discretize the super-large floating structure to obtain the overall stiffness matrix of the discretized super-large floating structure ;

[0008] S2: Based on the discretized global stiffness matrix of the super-large floating structure, the discrete frequency response function of the floating hydrodynamic system is calculated. ;

[0009] S3: Perform Laplace transform on the frequency response function of the floating hydrodynamic system and establish the pole-residue model matrix of the transfer function of the super-large floating structure system ;

[0010] S4: Based on the pole-residue model of the transfer function of the super-large floating structure system, the transient hydroelastic response of the super-large floating body is obtained .

[0011] In some embodiments, step S1 comprises:

[0012] S11: uniformly discretize the super-large floating structure to obtain N rigid submodules;

[0013] S12: Create N-1 beam units and connect the N rigid submodules, wherein each beam unit is connected to two rigid submodules respectively;

[0014] S13: Create the stiffness matrix corresponding to each beam element based on the finite element method and the connected rigid submodules , the beam element stiffness matrix Combined to obtain the overall stiffness matrix of the discretized super-large floating structure .

[0015] In some embodiments, the beam element stiffness matrix is ​​established ,in,

[0016] ;

[0017] ;

[0018] ;

[0019] ;

[0020] in, is the length of the beam unit along the x-axis, For beam units along The length of the axis, For beam units along The length of the axis, is the cross-sectional area of ​​the beam unit in the width direction, for relative to Moment of inertia of the axis, for relative to Moment of inertia of the axis, is Young's modulus, is the shear modulus, , and is the correction factor.

[0021] In some embodiments, the correction factor , , The calculation method is: ; ; .

[0022] In some embodiments, step S2 comprises:

[0023] S21: Based on initial conditions , , set the transient hydroelastic response of the super-large floating body , establish the time domain control equation of the super-large floating body system:

[0024] ,

[0025] in, It's time. is a dummy time variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the external load in the time domain, is the kernel of the convolution term;

[0026] S22, based on the time domain control equation of the super-large floating body system, the frequency domain control equation of the motion analysis of the super-large floating body system is established:

[0027] ,

[0028] in, is the frequency value, is the frequency-dependent added mass matrix, is the frequency-dependent damping matrix, is the external frequency domain load, is an imaginary unit, is the frequency domain motion response;

[0029] S23, according to the potential flow theory, calculate , and ;

[0030] S24, Frequency domain control equations and motion analysis based on very large floating body systems , and , calculate the discrete frequency response function of the floating hydrodynamic system :

[0031] .

[0032] In some embodiments, the frequency-dependent added mass matrix is ​​calculated as:

[0033] .

[0034] In some embodiments, the frequency-dependent damping matrix is ​​calculated as:

[0035] .

[0036] In some of these embodiments, according to the frequency-dependent added mass matrix , get the mass matrix ;

[0037] According to the frequency-dependent damping matrix , get the kernel of the convolution term .

[0038] In some embodiments, step S3 comprises:

[0039] S31, discrete frequency response function based on floating hydrodynamic system , obtain the frequency response function of any component of the very large floating body ;

[0040] S32, yes Perform inverse Fourier transform to obtain the system impulse response function of each component of the super-large floating body Approximate expression:

[0041] ;

[0042] in, Indicates at time The impulse response function at coordinate At the input The output at is the frequency step size, It indicates frequency The number of variables, is the frequency The number of is an imaginary unit, is a natural constant, is a frequency variable;

[0043] S33, yes Perform Laplace transform to obtain the pole-residue model elements of the transfer function of the super-large floating structure system ,Depend on Combine to get the pole-residue model matrix :

[0044] ;

[0045] in, is a Laplace domain variable, yes The elements in the matrix represent the system transfer function in the Laplace domain at the coordinates At the input The output at are variables representing the number of poles and residues, is the number of poles and residues, yes The remainder in yes The extreme point in The system transfer function At the input Where The corresponding Group poles and residues.

[0046] In some embodiments, step S4 comprises:

[0047] S41, establish the governing equations of a very large floating body in the Laplace domain:

[0048] ;

[0049] in, is a Laplace domain variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the delay function in the Laplace domain, is the motion response in the Laplace domain, is the load outside the Laplace domain;

[0050] S42, in the Laplace domain, calculate the transient response at rest ;

[0051] in, is the pole-residue model matrix of the transfer function of the very large floating structure system;

[0052] S43, input external load signal , external load signal Performing Laplace transform, we get Laplace transform of :

[0053] ;

[0054] in, is the external load in the Laplace domain, indicating the coordinates in the Laplace domain The input at Indicated in Variables for the number of extreme points and residues, yes The remainder in yes The extreme point in Indicated in coordinates The external load at Group poles and residues, is the total number of extreme points and residues;

[0055] S44, Laplace transform based on external load signal and the pole-residue model of the transfer function of the very large floating structure system to calculate the hydroelastic response of the very large floating structure system in the Laplace domain;

[0056] ;

[0057] in, is the floating body in coordinates The Laplace domain motion response at are variables representing the number of poles and residues, yes The remainder in yes The extreme point in Indicated in coordinates At the input The motion response at Group poles and residues;

[0058] S45, based on , obtain a super large floating body Transient hydroelastic response ,Depend on Combined into .

[0059] Compared with the related technology, the embodiment of the present application provides a Laplace domain calculation method for the transient hydroelastic response of a very large floating body. By discretizing the very large floating body structure, a pole-residue model matrix of the transfer function of the very large floating body structure system is constructed to solve the response, and an analytical solution of the transient hydroelastic response of the very large floating body represented by a complex exponential form is obtained. This overcomes the technical defects of the traditional frequency domain analysis method that can only solve the steady-state response and the time domain method that the numerical integration is very time-consuming. This provides a new Laplace domain method with both high precision and high efficiency for the preliminary engineering analysis of very large floating bodies.

[0060] Details of one or more embodiments of the present application are set forth in the following drawings and description to make other features, objects, and advantages of the present application more readily apparent. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0062] Figure 1 A flow chart of a Laplace domain calculation method for a transient hydroelastic response of a very large floating body according to an embodiment of the present application;

[0063] Figure 2 This is a diagram of a numerical model of a super-large floating body according to an embodiment of the present application;

[0064] Figure 3 This is a discretization model diagram of a super-large floating body according to an embodiment of the present application;

[0065] Figure 4 A comparison diagram of the calculation results of the system impulse response function of an embodiment of the present application;

[0066] Figure 5 This is an impact load diagram of an embodiment of the present application;

[0067] Figure 6 This is a comparison chart of the transient heave response of the first submodule of the super-large floating body in an embodiment of the present application and the calculation results of the time domain method. DETAILED DESCRIPTION

[0068] In order to make the purpose, technical solutions and advantages of the present application clearer, the present application is described and illustrated below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not intended to limit the present application. Based on the embodiments provided in the present application, all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of the present application.

[0069] Obviously, the drawings described below are only some examples or embodiments of the present application. For ordinary technicians in this field, the present application can also be applied to other similar scenarios based on these drawings without creative work. In addition, it can also be understood that although the efforts made in this development process may be complicated and lengthy, for ordinary technicians in this field related to the content disclosed in this application, some changes in design, manufacturing or production based on the technical content disclosed in this application are just conventional technical means, and should not be understood as insufficient content disclosed in this application.

[0070] Reference to "embodiments" in this application means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those of ordinary skill in the art that the embodiments described in this application may be combined with other embodiments without conflict.

[0071] Unless otherwise defined, the technical terms or scientific terms involved in this application should be understood by people with ordinary skills in the technical field to which this application belongs. The words "one", "a", "a", "the" and the like involved in this application do not indicate a quantity limitation, and may indicate the singular or plural. The terms "include", "comprise", "have" and any of their variations involved in this application are intended to cover non-exclusive inclusions; for example, a process, method, system, product or device that includes a series of steps or modules (units) is not limited to the listed steps or units, but may also include steps or units that are not listed, or may also include other steps or units inherent to these processes, methods, products or devices. The words "connect", "connected", "coupled" and the like involved in this application are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The "multiple" involved in this application refers to two or more. "And / or" describes the association relationship of associated objects, indicating that there may be three relationships, for example, "A and / or B" can mean: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the objects before and after are in an "or" relationship. The terms "first", "second", "third", etc. involved in this application are only used to distinguish similar objects and do not represent a specific ordering of the objects.

[0072] This embodiment provides a Laplace domain calculation method for the transient hydroelastic response of a very large floating body. Figure 1 is a flow chart of a Laplace domain calculation method for a super-large floating body transient hydroelastic response according to an embodiment of the present application, such as Figure 1 As shown, the process includes the following steps:

[0073] S1: Discretize the super-large floating structure to obtain the overall stiffness matrix of the discretized super-large floating structure .

[0074] S2: Based on the discretized global stiffness matrix of the super-large floating structure, the discrete frequency response function of the floating hydrodynamic system is calculated. .

[0075] S3: Perform Laplace transform on the frequency response function of the floating hydrodynamic system and establish the pole-residue model matrix of the transfer function of the super-large floating structure system .

[0076] S4: Based on the pole-residue model of the transfer function of the super-large floating structure system, the transient hydroelastic response of the super-large floating body is obtained .

[0077] In some embodiments, step S1 comprises:

[0078] S11: Perform uniform discretization on the super-large floating structure to obtain N rigid submodules.

[0079] S12: Create N-1 beam units and connect the N rigid submodules, wherein each beam unit is connected to two rigid submodules.

[0080] S13: Create the stiffness matrix corresponding to each beam element based on the finite element method and the connected rigid submodules , the beam element stiffness matrix Combined to obtain the overall stiffness matrix of the discretized super-large floating structure .

[0081] N rigid submodules are connected using N-1 beam units to approximate the inherent structural deformation of the original flexible structure, while considering the coupling effect of the overall stiffness of the discretized multi-floating structure and the hydrodynamic coefficients of each submodule.

[0082] In some embodiments, the beam element stiffness matrix is ​​established ,in,

[0083] ;

[0084] ;

[0085] ;

[0086] ;

[0087] in, is the length of the beam unit along the x-axis, For beam units along The length of the axis, For beam units along The length of the axis, is the cross-sectional area of ​​the beam unit in the width direction, for relative to Moment of inertia of the axis, for relative to Moment of inertia of the axis, is Young's modulus, is the shear modulus, , and is the correction factor.

[0088] In the discretization of the global stiffness matrix of the very large floating structure The element stiffness matrix in To express, this formula represents the Nth beam unit stiffness matrix, a beam unit stiffness matrix The latitude is 12×12, discretizing the overall stiffness matrix of the super-large floating structure The latitude is .

[0089] In some embodiments, the correction factor , , The calculation method is: ; ; .

[0090] In some embodiments, step S2 comprises:

[0091] S21: Based on initial conditions , , set the transient hydroelastic response of the super-large floating body , establish the time domain control equation of the super-large floating body system:

[0092] ,

[0093] in, It's time. is a dummy time variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the external load in the time domain, is the kernel of the convolution term.

[0094] S22, based on the time domain control equation of the super-large floating body system, the frequency domain control equation of the motion analysis of the super-large floating body system is established:

[0095] ,

[0096] in, is the frequency value, is the frequency-dependent added mass matrix, is the frequency-dependent damping matrix, is the external frequency domain load, is an imaginary unit, is the frequency domain motion response.

[0097] S23, according to the potential flow theory, calculate , and .

[0098] S24, Frequency domain control equations and motion analysis based on very large floating body systems , and , calculate the discrete frequency response function of the floating hydrodynamic system :

[0099] .

[0100] In some embodiments, the frequency-dependent added mass matrix is ​​calculated as:

[0101] .

[0102] In some embodiments, the frequency-dependent damping matrix is ​​calculated as:

[0103] .

[0104] In some of these embodiments, according to the frequency-dependent added mass matrix , get the mass matrix .

[0105] According to the frequency-dependent damping matrix , get the kernel of the convolution term .

[0106] In some embodiments, step S3 comprises:

[0107] S31, discrete frequency response function based on floating hydrodynamic system , obtain the frequency response function of any component of the very large floating body .

[0108] S32, yes Perform inverse Fourier transform to obtain the system impulse response function of each component of the super-large floating body Approximate expression:

[0109] .

[0110] in, Indicates at time The impulse response function at coordinate At the input The output at is the frequency step size, It indicates frequency The number of variables, is the frequency The number of is an imaginary unit, is a natural constant, is a frequency variable.

[0111] S33, yes Perform Laplace transform to obtain the pole-residue model elements of the transfer function of the super-large floating structure system ,Depend on Combine to get the pole-residue model matrix :

[0112] .

[0113] in, is a Laplace domain variable, yes The elements in the matrix represent the system transfer function in the Laplace domain at the coordinates At the input The output at are variables representing the number of poles and residues, is the number of poles and residues, yes The remainder in yes The extreme point in The system transfer function At the input Where The corresponding Group poles and residues.

[0114] The value is usually 106.

[0115] In some embodiments, step S4 comprises:

[0116] S41, establish the governing equations of a very large floating body in the Laplace domain:

[0117] .

[0118] in, is a Laplace domain variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the delay function in the Laplace domain, is the motion response in the Laplace domain, is the load outside the Laplace domain.

[0119] S42, in the Laplace domain, calculate the transient response at rest .

[0120] in, It is the pole-residue model matrix of the transfer function of the very large floating structure system.

[0121] S43, input external load signal , for external load signal Performing Laplace transform, we get Laplace transform of :

[0122] ;

[0123] in, is the external load in the Laplace domain, indicating the coordinates in the Laplace domain The input at Indicated in Variables for the number of extreme points and residues, yes The remainder in yes The extreme point in Indicated in coordinates The external load at Group poles and residues, is the total number of extreme points and residues.

[0124] S44, Laplace transform based on external load signal and the pole-residue model of the transfer function of the very large floating structure system to calculate the hydroelastic response of the very large floating structure system in the Laplace domain;

[0125] .

[0126] in, is the floating body in coordinates The Laplace domain motion response at are variables representing the number of poles and residues, yes The remainder in yes The extreme point in Indicated in coordinates At the input The motion response at Group poles and residues;

[0127] S45, based on , obtain a super large floating body Transient hydroelastic response ,Depend on Combined into .

[0128] It should be noted that the steps shown in the above process or the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0129] In the embodiment, taking a certain super-large floating structure as an example, a continuous flexible box is selected, such as Figure 2 The flexible box is 300 m long, 60 m wide, 2 m high, has a draft of 0.5 m, and weighs 9.225 × 106 Kg. The xy plane of the coordinate system is located on the still water surface, and the z-axis is pointing upward. This embodiment calculates the instantaneous hydroelastic response of a super-large floating body under the excitation of a negative impact load along the z-axis at the center of gravity of the first submodule.

[0130] The super-large floating structure is discretized into 8 rigid submodules connected by 7 beam units, such as Figure 3 According to the finite element method and the correction factor, the beam unit stiffness matrix and the overall stiffness matrix of the discretized very large floating structure assembled from the unit stiffness matrix are obtained.

[0131] Based on the linear potential flow theory, the discrete frequency response function of the ultra-large floating system is calculated, where the frequency range is 0.1-2.5 rad / s and the frequency interval is = 0.001 rad / s, perform inverse Fourier transform and obtain the heave impulse response function ,like Figure 4 shown.

[0132] Consider external loads , calculate the poles and residues of the signal, the time domain signal obtained by the pole residue method is compared with the original load signal. Figure 5 The time domain signal obtained in the present application and the calculation method with different input and system poles provided in the present application can obtain an analytical solution for the heave response.

[0133] Figure 6 is the calculated heave response time history of the first submodule of the super-large floating body ,at the same time Figure 6 The response results calculated by the traditional time-domain numerical integration method are shown, and it can be seen that the two are in good agreement, proving that the present invention has high accuracy and effectiveness.

[0134] The technical features of the above-described embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above-described embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0135] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A Laplace domain calculation method for transient hydroelastic response of a very large floating body, characterized in that: The following steps are involved: S1: Discretize the super-large floating structure to obtain the overall stiffness matrix of the discretized super-large floating structure ; S2: Based on the discretized global stiffness matrix of the super-large floating structure, the discrete frequency response function of the floating hydrodynamic system is calculated. ; S3: Perform Laplace transform on the frequency response function of the floating hydrodynamic system and establish the pole-residue model matrix of the transfer function of the super-large floating structure system ; S4: Based on the pole-residue model of the transfer function of the super-large floating structure system, the transient hydroelastic response of the super-large floating body is obtained ; Wherein, the step S3 comprises: S31, discrete frequency response function based on floating hydrodynamic system , obtain the frequency response function of any component of the very large floating body ; S32, yes Perform inverse Fourier transform to obtain the system impulse response function of each component of the super-large floating body Approximate expression: ； in, Indicates at time The impulse response function at coordinate At the input The output at is the frequency step size, It indicates frequency The number of variables, is the frequency The number of is an imaginary unit, is a natural constant, is a frequency variable; S33, yes Perform Laplace transform to obtain the pole-residue model elements of the transfer function of the super-large floating structure system ,Depend on Combine to get the pole-residue model matrix : ； in, is a Laplace domain variable, yes The elements in the matrix represent the system transfer function in the Laplace domain at the coordinates At the input The output at are variables representing the number of poles and residues, is the number of poles and residues, yes The remainder in yes The extreme point in The system transfer function At the input Where The corresponding Group poles and residues; Wherein, the step S4 comprises: S41, establish the governing equations of a very large floating body in the Laplace domain: ； in, is a Laplace domain variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the delay function in the Laplace domain, is the motion response in the Laplace domain, is the load outside the Laplace domain; S42, in the Laplace domain, calculate the transient response at rest ; in, is the pole-residue model matrix of the transfer function of the very large floating structure system; S43, input external load signal , external load signal Performing Laplace transform, we get Laplace transform of : ； in, is the external load in the Laplace domain, indicating the coordinates in the Laplace domain The input at Indicated in Variables for the number of extreme points and residues, yes The remainder in yes The extreme point in Indicated in coordinates The external load at Group poles and residues, is the total number of extreme points and residues; S44, Laplace transform based on external load signal and the pole-residue model of the transfer function of the very large floating structure system to calculate the hydroelastic response of the very large floating structure system in the Laplace domain; ； in, is the floating body in coordinates The Laplace domain motion response at are variables representing the number of poles and residues, yes The remainder in yes The extreme point in Indicated in coordinates At the input The motion response at Group poles and residues; S45, based on , obtain a super large floating body Transient hydroelastic response ,Depend on Combined into .

2. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 1 is characterized in that: Step S1 includes: S11: uniformly discretize the super-large floating structure to obtain N rigid submodules; S12: Create N-1 beam units and connect the N rigid submodules, wherein each beam unit is connected to two rigid submodules respectively; S13: Create the stiffness matrix corresponding to each beam element based on the finite element method and the connected rigid submodules , the beam element stiffness matrix Combined to obtain the overall stiffness matrix of the discretized super-large floating structure .

3. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 2 is characterized in that: Establish beam element stiffness matrix ,in, ； ； ； ； in, is the length of the beam unit along the x-axis, For beam units along The length of the axis, For beam units along The length of the axis, is the cross-sectional area of ​​the beam unit in the width direction, for relative to Moment of inertia of the axis, for relative to Moment of inertia of the axis, is Young's modulus, is the shear modulus, , and is the correction factor.

4. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 3 is characterized in that: Correction Factor , , The calculation method is: ; ; .

5. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 1 is characterized in that: Step S2 includes: S21: Based on initial conditions , , set the transient hydroelastic response of the super-large floating body , establish the time domain control equation of the super-large floating body system: ， in, It's time. is a dummy time variable, is the mass matrix of the very large floating body system, is the added mass matrix, is the hydrostatic restoring coefficient matrix, is the external load in the time domain, is the kernel of the convolution term; S22, based on the time domain control equation of the super-large floating body system, the frequency domain control equation of the motion analysis of the super-large floating body system is established: ， in, is the frequency value, is the frequency-dependent added mass matrix, is the frequency-dependent damping matrix, is the external frequency domain load, is an imaginary unit, is the frequency domain motion response; S23, according to potential flow theory, calculate , and ; S24, Frequency domain control equations and motion analysis based on very large floating body systems , and , calculate the discrete frequency response function of the floating hydrodynamic system : 。 6. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 5 is characterized in that: The frequency-dependent added mass matrix is ​​calculated as: 。 7. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 5, characterized in that: The frequency-dependent damping matrix is ​​calculated as: 。 8. The Laplace domain calculation method for transient hydroelastic response of a very large floating body according to claim 5, characterized in that: According to the frequency-dependent additional mass matrix , get the mass matrix ; According to the frequency-dependent damping matrix , get the kernel of the convolution term .

Images