System and method for fast prediction of maximum impact pressure of rectangular tank top cover with high liquid filling rate

Abstract

This invention discloses a rapid prediction system and method for the maximum impact pressure of the top cover of a rectangular tank with a high filling rate. The method includes the following steps: acquiring various parameters of the rectangular tank and the liquid; calculating the first-order natural angular frequency of the rectangular tank's swaying motion and the increment of the free surface wave height within a single half-excitation cycle; calculating the swaying period parameters corresponding to the first contact with the top; calculating the maximum swaying wave height; introducing an effective wet contact length coefficient to calculate the wet contact length; integrating the vertical velocity of the liquid along the wet contact liquid region to calculate the momentum of the wet contact liquid region; the impact force of the top cover is determined by the rate of change of the momentum of the wet contact liquid region with respect to time, and the dynamic impact pressure of the top cover is obtained by dividing by the wet contact area; and outputting the maximum dynamic impact pressure. This invention establishes the relationship between the momentum change of the wet contact liquid region and the impact pressure of the top cover, introduces swaying period parameters and an effective wet contact length coefficient, and rapidly estimates the maximum dynamic impact pressure of the top cover.

Description

Technical Field

[0001] This invention relates to the field of safety assessment technology for liquid storage structures, and in particular to a rapid prediction system and method for the maximum impact pressure of the top cover of a rectangular storage tank with a high filling rate. Background Technology

[0002] Under external dynamic excitation, the free liquid surface of a high-fill-rate rectangular tank is prone to significant sloshing. When the sloshing wave height exceeds the top cover clearance, the liquid will impact the top cover with a short-duration, high-amplitude force, forming dynamic impact pressure on the top cover. For passive containment cooling tanks, emergency water tanks, and other safety-related liquid storage structures in nuclear power engineering, the top cover impact pressure may affect the structural integrity of the top cover and connecting areas. Therefore, it is necessary to quickly determine the maximum dynamic impact pressure of the top cover during the engineering design and safety assessment stages.

[0003] While existing high-precision numerical simulation methods can describe the evolution of complex free liquid surfaces, they are complex to model, have high computational costs, and are sensitive to mesh or particle parameters, which is not conducive to the comparison and rapid verification of engineering schemes. Traditional simplified models are mostly focused on the dynamic water pressure of the sidewall or the response to small swaying, which is difficult to directly reflect the momentum change characteristics of the wet contact area impacting the top cover under high filling rate conditions.

[0004] Therefore, the technical problem to be solved by the present invention is to provide a method for rapid prediction of the maximum dynamic impact pressure of the top cover of a high-filling-rate rectangular tank with a simple calculation process, clear physical meaning of parameters, and calibration by combining experimental statistics, and to further use it for the impact resistance design verification of the top cover and engineering safety assessment. Summary of the Invention

[0005] Purpose of the invention: In order to overcome the shortcomings of the prior art, the present invention provides a rapid prediction system and method for the maximum impact pressure of the top cover of a rectangular tank with high filling rate. By establishing the relationship between the momentum change of the wet contact liquid area and the impact pressure of the top cover, and introducing the swaying period parameter and the effective wet contact length coefficient, the maximum dynamic impact pressure of the top cover can be rapidly estimated.

[0006] Technical Solution: To achieve the above objective, the present invention provides a rapid prediction method for the maximum impact pressure of the top cover of a high-fill-rate rectangular tank, comprising the following steps: S1: Obtaining various parameters of the rectangular tank and the liquid; S2: Calculating the first-order swaying natural angular frequency of the rectangular tank. And determine the external excitation angular frequency. With the first-order wobbling natural angular frequency The degree of proximity between them; S3: When the external excitation is in a frequency band that easily causes large sloshing of the liquid, calculate the increment of the free surface wave height within a single half-excitation cycle. ;like If the wave height is less than the preset threshold and the calculated wave height cannot reach the top cover, the output result is "no impact on the top cover" or "low impact risk". If the calculated wave height can reach the top cover, the swaying period parameter corresponding to the first impact is calculated. S4: Swaying period parameters based on the first touch-the-top. Calculate the maximum sway wave height S5: Introducing an effective wet contact length coefficient And calculate the wet contact length. S6: Vertical velocity of the liquid Integrating along the wet contact liquid region, combined with the maximum sloshing wave height and wet contact length The value of is used to calculate the momentum of the wet contact liquid region. S7: Characterize the impact process after the liquid comes into contact with the top cover as the change in momentum of the wet contact liquid region over time; Top cover impact force Momentum of the wet contact liquid region Regarding time The rate of change is determined, and the dynamic impact pressure of the top cover is determined. Depend on S8: Maximum dynamic impact pressure output, obtained by dividing by the wet contact area; .

[0007] Furthermore, in S1, the parameters that need to be obtained include the half-length of the storage tank. Storage tank width Liquid depth Freeboard height Liquid density and external excitation amplitude .

[0008] Furthermore, in S2, the first-order wobbling natural angular frequency The calculation formula is:

[0009] .

[0010] Furthermore, in S3, the free surface wave height increment within a single half-excitation cycle The calculation formula is:

[0011] ;

[0012] in, For gravitational acceleration, the external excitation angular frequency Take the first-order natural angular frequency of the oscillation The oscillation period parameters corresponding to the first touch of the top The calculation formula is:

[0013] .

[0014] Furthermore, in S4, considering the energy accumulation and impact peak hysteresis of the liquid under continuous excitation, a swaying period parameter is introduced for calculating the maximum dynamic impact pressure. , ,in For peak impact hysteresis correction term; maximum sway wave height The calculation formula is:

[0015] .

[0016] Furthermore, in S5, the effective wet contact length coefficient The relative wet contact length characterizes the effective impact between the liquid and the top cover, and is defined as follows: The wet contact length is calculated from this. .

[0017] Furthermore, in S6, the momentum of the wet contact liquid region... The calculation formula is as follows:

[0018] ;

[0019] Among them, the vertical velocity of the liquid The expression is:

[0020] ;

[0021] Momentum of the wet contact liquid region The calculation formula, when expanded, yields:

[0022] .

[0023] Furthermore, in the S7, the top cover impact force The calculation formula is:

[0024] ;

[0025] After unfolding, we get:

[0026] ;

[0027] Top cover dynamic impact pressure Depend on Divide by wet contact area The calculation formula is as follows:

[0028] .

[0029] Furthermore, in S8, the maximum dynamic impact pressure MDIP is calculated in the following parameterized form:

[0030] ;

[0031] For a rectangular tank with a high filling rate under conditions of rigid walls, horizontal harmonic excitation, and first-order swaying mode dominance, let , The expanded formula is obtained as follows:

[0032] .

[0033] Furthermore, the rapid prediction system for the maximum impact pressure of the top cover of a high-fill-rate rectangular tank includes a parameter input module for inputting tank geometric parameters, liquid parameters, and external excitation parameters; and a swaying parameter calculation module for calculating the first-order swaying natural angular frequency. Wave height increment First peak parameters and maximum wave height The parameter calibration module is used to determine the oscillation period parameters. and effective wet contact length coefficient The impact pressure calculation module is used to calculate the impact force of the top cover based on the rate of change of momentum of the wet contact liquid and convert it into MDIP; the safety verification module is used to compare MDIP with the top cover design limit, connection area bearing limit or specification limit; the result output module is used to output MDIP, parameter values, safety factor, risk level and design verification report.

[0034] Beneficial effects: The rapid prediction system and method for the maximum impact pressure of the top cover of a high-filling-rate rectangular tank of the present invention characterizes the impact process after the liquid comes into contact with the top cover as the change of momentum of the wet contact liquid area over time, and determines the impact force of the top cover by the momentum change rate, and then calculates the impact pressure of the top cover by combining the wet contact area. Thus, a parameterized expression of the maximum dynamic impact pressure of the top cover can be obtained, so that engineers can obtain the estimated result of the impact pressure of the top cover without conducting high-cost numerical simulations. Attached Figure Description

[0035] Appendix Figure 1 This is a two-dimensional schematic diagram of a rectangular storage tank. Detailed Implementation

[0036] The invention will now be further described with reference to the accompanying drawings.

[0037] As attached Figure 1 The method for rapidly predicting the maximum impact pressure of the top cover of a high-filling-rate rectangular tank includes the following steps.

[0038] S1: Obtain various parameters of the rectangular tank and the liquid. Specifically, the parameters to be obtained include the half-length of the tank. Storage tank width Liquid depth Freeboard height Liquid density and external excitation amplitude As attached Figure 1 In the diagram, label 1 indicates a rectangular tank structure, label 2 indicates the top cover of the rectangular tank, label 3 indicates the liquid sloshing inside the tank, and label 4 indicates the liquid level when the liquid is still inside the tank. Indicates the sloshing wave height of liquid in the tank.

[0039] S2: Calculate the first-order natural angular frequency of the rectangular tank's swaying motion. For horizontal excitation conditions dominated by the first-order sway mode, calculate or read the first-order sway natural angular frequency of the rectangular tank. And determine the external excitation angular frequency. With the first-order wobbling natural angular frequency The degree of closeness between them. Among them, the first-order wobbling natural angular frequency. The calculation formula is:

[0040] ;

[0041] in, This is the acceleration due to gravity.

[0042] S3: When the external excitation is in a frequency band that easily causes large sloshing of the liquid, calculate the free surface wave height increment within a single half-excitation cycle. .like If the wave height is less than the preset threshold and the calculated wave height cannot reach the top cover, the output result is "no top cover impact or low impact risk". If the calculated wave height can reach the top cover, the swaying period parameter corresponding to the first impact is calculated. Specifically, the increase in free surface wave height within a single half-excitation cycle. The calculation formula is:

[0043] ;

[0044] in, For gravitational acceleration, the external excitation angular frequency Take the first-order natural angular frequency of the oscillation The oscillation period parameters corresponding to the first touch of the top. The calculation formula is as follows:

[0045] .

[0046] S4: Based on the oscillation period parameters corresponding to the first touch-the-top Calculate the maximum sway wave height Because the energy accumulation and impact peak hysteresis of the liquid under continuous excitation need to be considered, a swaying period parameter for calculating the maximum dynamic impact pressure must be introduced. , ,in The impact peak hysteresis correction term represents the number of additional half-excitation cycles required for the liquid to continue developing to the state corresponding to the maximum dynamic impact pressure after the free surface of the liquid first touches the top. The preferred value is 2, but 3 can also be chosen based on safety margin or calibration error. Therefore, the maximum sway wave height... The calculation formula is:

[0047] .

[0048] This invention takes into account that the maximum dynamic impact pressure may not necessarily occur when the liquid first contacts the top cover, but may occur during several subsequent swaying peaks. Therefore, the initial contact period parameters are first determined based on the freeboard height and the free surface wave height increment within a single half-cycle. Then, through the impact peak lag correction term Determine the swaying period parameters for predicting maximum dynamic shock pressure. This describes the accumulation of liquid sloshing energy and the development of impact pressure under continuous external excitation.

[0049] The value of is not given arbitrarily; firstly This represents the sloshing period parameter corresponding to the theoretical first contact of the free surface with the top cover. At this point, the liquid has just touched the top cover, the wet contact length is small, and the liquid area involved in momentum exchange is limited. Therefore, the moment of first contact usually does not correspond to the maximum dynamic impact pressure. Under continuous harmonic excitation, the sloshing of the liquid still requires a development process from the first contact to the formation of the maximum impact pressure. As subsequent half-excitation cycles continue, the free surface rise, upward velocity, wet contact length, and momentum exchange in the wet contact liquid area will all further increase. Therefore, the maximum dynamic impact pressure often lags behind the moment of first contact. This lag can be corrected using an additional impact peak lag term. Characterization. Secondly, The value should not be too large. This is because actual liquid sloshing is not an ideal, undamped, linear growth process. Liquid viscous dissipation, local splashing, gas-liquid interaction, wall constraint, and energy loss after the impact of the top cover all weaken subsequent energy accumulation, preventing the impact pressure from increasing linearly. It increases indefinitely with the increase of [something]. Therefore, A finite correction value that reflects the continued accumulation of energy shortly after the initial contact with the surface should be selected. Based on extensive experimental results, the maximum dynamic impact pressure typically occurs during several peak wave periods following the initial contact with the surface. Error comparisons of different filling depths, excitation frequencies, and excitation amplitudes revealed that when… When the value is 2 or 3, the calculated value matches well with the maximum dynamic impact pressure measured in the experiment.

[0050] S5: Calculate wet contact length First, we introduce the effective wet contact length coefficient. , The relative wet contact length characterizes the effective impact force between the liquid and the top cover, and further characterizes the effective range of the liquid's actual action on the top cover, and is defined as follows: The wet contact length is calculated from this. ,in It is half the length of the storage tank. The recommended value range is 0.16 to 0.19, with an optimal value of 0.18.

[0051] When a large liquid sloshes and impacts the top cover, the actual pressure distribution exhibits significant localization; not all geometrically wetted contact areas contribute equally to the maximum dynamic impact pressure. The maximum impact pressure is primarily controlled by the high-speed upward-rushing liquid near the sidewalls or the impact leading edge. If the liquid size is too small, the equivalent contact area is too narrow, making it easy to miss liquid regions that still possess high vertical velocity and strong momentum exchange capacity; when When the effective contact area is too large, it extends excessively into the tank interior, incorporating liquid regions with lower vertical velocities, weaker top-contact synchronicity, and smaller impact contributions, thus weakening the ability to characterize local maximum impact pressure. Therefore, This can be understood as the reasonable range of the main high-momentum impact core area near the covered sidewall.

[0052] It is not a constant that remains unchanged for all rectangular storage tanks. The preferred value of 0.18 is obtained through back-calculation and error optimization based on multiple sets of experimental data under the conditions of a high-filling-ratio rectangular tank, rigid walls, horizontal harmonic excitation, and the dominance of the first-order swaying mode in this embodiment. This range is of reference significance for rectangular tanks with similar geometric proportions, filling ratios, freeboard heights, and excitation characteristics, but should not be construed as being applicable to all rectangular tanks. If there are significant changes in tank size proportions, liquid level, freeboard height, top cover type, liquid properties, or excitation conditions, recalibration should be performed through testing. .

[0053] S6: Increase the vertical velocity of the liquid Integrating along the wet contact liquid region, combined with the maximum sloshing wave height and wet contact length The value of is used to calculate the momentum of the wet contact liquid region. Specifically, the momentum of the wet contact liquid region. The calculation formula is:

[0054] ;

[0055] After unfolding it, we get:

[0056] .

[0057] Among them, the vertical velocity of the liquid The expression is:

[0058] .

[0059] In the formula, and This represents the two directions of a two-dimensional computational coordinate system. (See attached image.) Figure 1 As shown in the figure, a two-dimensional rectangular coordinate system O-XZ is established, where O is the midpoint of the inner surface of the rectangular tank bottom plate. The X-axis is arranged along the horizontal width of the tank, with the positive direction pointing to the right side wall; the Z-axis is arranged vertically, with the positive direction pointing from the bottom plate to the top cover. This coordinate system is used to describe the shape of the free surface of the liquid sloshing, the contact area between the liquid and the top cover, and the location of the impact, but does not represent the solid components in the tank. The velocity potential of the liquid is expressed as:

[0060] .

[0061] S7: The impact process after the liquid comes into contact with the top cover is characterized as the change of momentum of the wet contact liquid region over time, and then the impact force of the top cover is calculated. Dynamic impact pressure of the top cover Among them, the impact force of the top cover Momentum of the wet contact liquid region Regarding time The rate of change is determined, and the dynamic impact pressure of the top cover is determined. Then by Divide by the wet contact area to obtain. Specifically, the impact force of the top cover. The calculation formula is:

[0062] ;

[0063] After unfolding it, we get:

[0064] ;

[0065] Top cover dynamic impact pressure Depend on Divide by wet contact area The calculation formula is as follows:

[0066] .

[0067] Therefore, this invention does not directly use complex numerical flow field simulation methods, but characterizes the impact process after the liquid comes into contact with the top cover as the change process of momentum of the wet contact liquid area over time, and determines the impact force of the top cover by the momentum change rate, and then calculates the impact pressure of the top cover by combining the wet contact area.

[0068] S8: Output maximum dynamic impact pressure The maximum dynamic impact pressure MDIP is calculated using the following parametric form: ;in, For the oscillation period parameter, The effective wet contact length coefficient, The length of the storage tank is half that of the tank. For liquid depth, Freeboard height For the density of the liquid, To excite the amplitude, Let be the external excitation angular frequency. For a high-fill-rate rectangular tank with a rigid wall, horizontal harmonic excitation, and a first-order swaying mode dominating, let . , Then adjust the oscillation period parameter and effective wet contact length coefficient Substituting the value into the dynamic impact pressure of the top cover The calculation formula, when expanded, yields:

[0069] .

[0070] The present invention also provides a rapid prediction system for the maximum impact pressure of the top cover of a rectangular storage tank with a high filling rate, including a parameter input module, a swaying parameter calculation module, a parameter calibration module, an impact pressure calculation module, a safety verification module, and a result output module.

[0071] The parameter input module is used to input the tank's geometric parameters, liquid parameters, and external excitation parameters. The swaying parameter calculation module is used to calculate the first-order swaying natural angular frequency. Wave height increment First peak parameters and maximum wave height The parameter calibration module is used to determine the oscillation period parameters. and effective wet contact length coefficient The impact pressure calculation module calculates the impact force of the top cover based on the rate of change of momentum of the wet contact liquid, and converts it into MDIP. The safety verification module compares the MDIP with the top cover design limit, connection area load limit, or specification limit. The results output module outputs the MDIP, parameter values, safety factor, risk level, and design verification report.

[0072] In one embodiment, the tank is half the length 0.5m, tank width 0.03m, liquid depth 0.90m, freeboard height 0.1m, liquid density 1000 kg / m 3 Excitation amplitude The acceleration due to gravity is 0.050m. 9.8 m / s 2 .

[0073] First, the parameter input module acquires the aforementioned geometric parameters, liquid parameters, and excitation parameters, and calculates the first-order swaying natural angular frequency. .

[0074] .

[0075] External excitation angular frequency Take the first-order natural angular frequency of the oscillation ,Right now Subsequently, the swaying parameter calculation module calculates the wave height increment of the free surface during each half-excitation cycle. .

[0076] .

[0077] according to Obtain the oscillation period parameters corresponding to the first touch-the-top. .

[0078] Secondly, the parameter calibration module determines the parameters based on experimental statistical laws. At the same time, the effective wet contact length coefficient The value is set to 0.18.

[0079] Then, the maximum dynamic impact pressure MDIP of the top cover is calculated using the simplified analytical formula.

[0080] ;

[0081] After substituting, we get =48.27 kPa.

[0082] Finally, the safety verification module compares the calculated MDIP with the design allowable pressure of the top cover plate, stiffeners, or connection areas. If the MDIP is less than the design allowable value, the output indicates that the design requirements are met; if the MDIP is greater than the design allowable value, the output provides design suggestions to increase the top cover thickness, optimize the stiffener arrangement, improve connection strength, or adjust the liquid level operating range.

[0083] Existing technologies mainly include theoretical analytical methods, experimental methods, and numerical simulation methods. Theoretical analytical methods are usually based on potential flow theory, linear sloshing theory, or multimodal expansion methods, and are suitable for regular geometric containers and small-amplitude sloshing problems, and can predict the free surface morphology and sidewall hydrodynamic pressure relatively well. However, when liquid directly impacts the top cover under high filling rate conditions, the impact process involves free surface rise, expansion of wet contact area, local velocity abrupt change, and instantaneous momentum exchange. Traditional analytical methods cannot directly provide a stable, concise, and engineering-friendly expression for the maximum impact pressure of the top cover.

[0084] Experimental methods can directly obtain the impact pressure time history at the measuring point on the top cover, which is an important means of revealing the impact mechanism and verifying theoretical models. However, if design is based entirely on experimental data, a large number of repeated tests need to be carried out for different tank sizes, liquid depths, excitation frequencies, and excitation amplitudes. This is time-consuming, costly, and makes it difficult to conduct rapid parametric evaluation in the early stages of engineering design.

[0085] Numerical simulation methods include Computational Fluid Dynamics (CFD), Arbitrary Lagrangian-Eulerian (ALE) methods, and Smoothed Particle Hydrodynamics (SPH) methods. These methods can handle free surface deformation, local splashing, and strongly nonlinear impact processes, but they are computationally expensive and highly sensitive to mesh, interface reconstruction, particle resolution, boundary conditions, and fluid-structure interaction parameters. For rapid engineering evaluation of nuclear power plant tanks or large liquid storage structures, directly employing high-fidelity numerical simulations is often not economical.

[0086] This invention does not require the establishment of complex CFD or SPH numerical models. It can directly calculate MDIP based on tank geometric parameters and excitation parameters, making it suitable for scheme comparison and rapid verification in the early stages of design.

[0087] Furthermore, existing simplified calculation methods for roof impact pressure generally only provide theoretical expressions or require engineers to determine key parameters separately, lacking a parameter determination process that can simultaneously consider the impact stage, energy accumulation effect, and effective wet contact range. Therefore, in this invention, a swaying period parameter is introduced. and effective wet contact length coefficient , Characterizing the development stages and energy accumulation process of the impact, Characterizing the effective wet contact range facilitates understanding, calibration, and application by engineers, and also allows for the use of shaking table test or field monitoring data. and Calibration is performed to create an engineering model that is adapted to the specific tank conditions.

[0088] This invention achieves rapid, interpretable, and engineering-based assessment of the impact pressure on the top cover through a process of "tank parameter input—first-order sway parameter calculation—wet contact area determination—momentum change modeling—test calibration parameter correction—maximum dynamic impact pressure output—structural safety verification".

[0089] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A rapid prediction method for the maximum impact pressure of the top cover of a rectangular storage tank with a high filling rate, characterized in that: Includes the following steps: S1: Obtain various parameters of the rectangular storage tank and the liquid; S2: Calculate the first order sloshing natural angular frequency of the rectangular tank and determine the closeness between the external excitation angular frequency and the first order sloshing natural angular frequency ​ S3: When the external excitation is in the frequency band that easily causes large liquid sloshing, calculate the free surface wave height increment in a single half-excitation period ; if the calculated wave height is less than the preset threshold value and cannot reach the top cover, output a result that the top cover impact does not occur or the risk of impact is low; if the calculated wave height can reach the top cover, calculate the sloshing period parameter corresponding to the first top touch ; S4: calculate the maximum sway wave height based on the sway period parameter corresponding to the first touch top , calculate the maximum sway wave height ; S5: Introducing an effective wet contact length coefficient and calculating the wet contact length ; S6: liquid vertical velocity integrated along the wet contact liquid area, combined with the maximum sloshing wave height and the wet contact length to calculate the wet contact liquid area momentum ; S7: The impact process after the liquid comes into contact with the top cover is characterized as the change of momentum of the wet contact liquid region over time; top cover impact force. Momentum of the wet contact liquid region Regarding time The rate of change is determined, and the dynamic impact pressure of the top cover is determined. Depend on Divide by the wet contact area to obtain; S8: Output maximum dynamic impact pressure .

2. The method for rapid prediction of the maximum impact pressure of the top cover of a high-fill-rate rectangular storage tank according to claim 1, characterized in that: In S1, the parameters that need to be obtained include the half-length of the tank. Storage tank width Liquid depth Freeboard height Liquid density and external excitation amplitude .

3. The method for rapid prediction of the maximum impact pressure of the top cover of a high-filling-rate rectangular storage tank according to claim 2, characterized in that: In S2, the first-order wobbling natural angular frequency The calculation formula is: 。 4. The method for rapid prediction of the maximum impact pressure of the top cover of a high-filling-rate rectangular storage tank according to claim 3, characterized in that: In S3, the free surface wave height increment within a single half-excitation cycle The calculation formula is: ； in, For gravitational acceleration, the external excitation angular frequency Take the first-order natural angular frequency of the oscillation The oscillation period parameters corresponding to the first touch of the top The calculation formula is: 。 5. The method for rapid prediction of the maximum impact pressure of the top cover of a high-filling-rate rectangular storage tank according to claim 4, characterized in that: In S4, considering the energy accumulation and impact peak hysteresis of the liquid under continuous excitation, a swaying period parameter is introduced for calculating the maximum dynamic impact pressure. , ,in For peak impact hysteresis correction term; maximum sway wave height The calculation formula is: 。 6. The method for rapid prediction of the maximum impact pressure of the top cover of a high-filling-rate rectangular storage tank according to claim 5, characterized in that: In S5, the effective wet contact length factor The relative wet contact length characterizes the effective impact between the liquid and the top cover, and is defined as follows: The wet contact length is calculated from this. .

7. The method for rapid prediction of the maximum impact pressure of the top cover of a high-fill-rate rectangular storage tank according to claim 6, characterized in that: In S6, the momentum of the wet contact liquid region The calculation formula is as follows: ； Among them, the vertical velocity of the liquid The expression is: ； Momentum of the wet contact liquid region The calculation formula, when expanded, yields: 。 8. The method for rapid prediction of the maximum impact pressure of the top cover of a high-fill-rate rectangular storage tank according to claim 7, characterized in that: In the S7, the impact force of the top cover The calculation formula is: ； After unfolding, we get: ； Top cover dynamic impact pressure Depend on Divide by wet contact area The calculation formula is as follows: 。 9. The method for rapid prediction of the maximum impact pressure of the top cover of a high-fill-rate rectangular storage tank according to claim 8, characterized in that: In S8, the maximum dynamic impact pressure MDIP is calculated in the following parametric form: ； For a rectangular tank with a high filling rate under conditions of rigid walls, horizontal harmonic excitation, and first-order swaying mode dominance, let , The expanded formula is obtained as follows: 。 10. The rapid prediction system for the maximum impact pressure of the top cover of a high-fill-rate rectangular tank according to any one of claims 1 to 9, characterized in that: Includes a parameter input module for inputting tank geometry parameters, liquid parameters, and external excitation parameters; The sway parameter calculation module is used to calculate the first-order sway natural angular frequency. Wave height increment First peak parameters and maximum wave height ; The parameter calibration module is used to determine the oscillation period parameters. and effective wet contact length coefficient ; The impact pressure calculation module is used to calculate the top cover impact force based on the rate of change of momentum of the wet contact liquid and convert it into MDIP. The safety verification module is used to compare the MDIP with the top cover design limit, connection area load limit, or specification limit. The results output module is used to output MDIP, parameter values, safety factor, risk level, and design verification report.

Images