Ice load prediction method and device in structure icebreaking process and electronic equipment

Abstract

The invention provides an ice load prediction method and device in the structure icebreaking process and electronic equipment, and belongs to the technical field of ice load prediction.The method comprises the steps that size parameters and material parameters of a target ice layer are obtained, and the material parameters comprise bending strength; calculating a target stress of a target point on the target ice layer according to the to-be-determined load, the coordinates of the current concentrated load acting point and the size parameters, wherein the target stress is represented by the to-be-determined load; determining the equivalent stress of the target ice layer according to the target stress and the material parameters, wherein the equivalent stress is represented by a to-be-determined load; and predicting the load to be determined according to the bending strength and the equivalent stress. According to the ice load prediction method in the structure icebreaking process provided by the invention, the technical problem that the limit load of any point of the ice layer cannot be accurately predicted in the ice load prediction method in the structure icebreaking process in the prior art is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of ice load prediction, and in particular to a method, device and electronic equipment for predicting ice load in a process of breaking ice on a structure. Background Art

[0002] When launching a structure in an icy ocean, the ice layer will break after reaching the limit load when the structure contacts it. This process will affect the trajectory of the structure. Therefore, the study of the interaction between structures and ice has received extensive attention, involving a variety of nonlinear problems such as structural failure, large deformation, and crack propagation. It is of great significance to establish an accurate prediction model for the interaction between structures and ice and to achieve rapid prediction of ice load characteristics during the icebreaking process. In related technologies, when the structure interacts with ice, it is impossible to accurately predict the limit load at any point on the ice layer.

[0003] It can be seen that the ice load prediction method of the structure breaking ice process in the related art has the technical problem of being unable to accurately predict the ultimate load at any point of the ice layer. Summary of the invention

[0004] The present invention provides a method, device and electronic equipment for predicting ice load during the ice-breaking process of a structure, so as to solve the technical problem that the ice load prediction method for the ice-breaking process of a structure in the related art cannot accurately predict the ultimate load at any point of the ice layer.

[0005] The present invention provides an ice load prediction method for a structure's icebreaking process, comprising the following steps: obtaining size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength; calculating a target stress of a target point on the target ice layer according to a load to be determined, coordinates of a current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are coordinates of a point on which the target ice layer is acted upon by a concentrated load of a structure, the load to be determined is a maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; determining an equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; and predicting the load to be determined according to the bending strength and the equivalent stress.

[0006] According to a method for predicting ice loads during icebreaking of a structure provided by the present invention, the target stress of a target point on the target ice layer is calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameter, including: determining the disturbance stress according to the load to be determined, the coordinates of the current concentrated load action point and the size parameter According to the disturbance Determine the target stress: ; in, is the normal stress of the target point along the X direction, is the normal stress of the target point along the Y direction, and together constitute the target stress; is the elastic modulus; z is the distance from the target point to the middle surface of the ice layer along the Z direction; is Poisson's ratio; is the disturbance; For the disturbance The second partial derivative with respect to x is For the disturbance The second partial derivative of y, x is the horizontal coordinate of the target point, and y is the vertical coordinate of the target point.

[0007] According to a method for predicting ice loads during ice breaking of a structure provided by the present invention, the disturbance is determined according to the load to be determined, the coordinates of the current concentrated load action point and the size parameter. , including: determining according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters : ; Based on the simply supported boundary condition, a double trigonometric series is used to represent : ; in, is the coefficient to be determined; a, b, h are the size parameters, a is the length of the target ice layer, b is the width of the target ice layer, and h is the thickness of the target ice layer; F is the load to be determined; is the horizontal coordinate of the current concentrated load action point coordinates, is the ordinate of the current concentrated load action point; m and n are the number of series expansion terms; x is the abscissa of the target point, and y is the ordinate of the target point.

[0008] According to a method for predicting ice load during ice breaking of a structure provided by the present invention, the Determining the target stress comprises: Substituting into the calculation formula of the target stress, we get The calculation formula is as follows: ; make , the target stress is: .

[0009] According to a method for predicting ice loads during icebreaking of a structure provided by the present invention, the equivalent stress of the target ice layer is determined according to the target stress and the material parameters, including: hydrostatic pressure , based on the small perturbation theory of elastic thin plates, , the hydrostatic pressure is determined according to the target stress: ; Based on the JH-2 constitutive model, the equivalent stress is determined according to the hydrostatic pressure and the material parameters; wherein, is the normal stress of the target point along the Z direction; P is the hydrostatic pressure.

[0010] According to a method for predicting ice loads during icebreaking of a structure provided by the present invention, the equivalent stress is determined based on the JH-2 constitutive model according to the hydrostatic pressure and the material parameters, including: normalizing the equivalent stress based on the JH-2 constitutive model for: ; When the target ice layer is not destroyed, is 0, according to the material parameters: ; The equivalent stress is determined based on the hydrostatic pressure: ; in, is the normalized equivalent stress of the target ice layer; is the normalized equivalent stress when the target ice layer is not destroyed; is the normalized equivalent stress when the target ice layer is completely destroyed; is the loss factor, when the target ice layer is not destroyed is 0, when the target ice layer is completely destroyed is 1; A, C, N are the material parameters; is the normalized hydrostatic pressure; is the normalized maximum tensile stress of the target ice layer; is the normalized strain rate; is the equivalent stress at the Hugoniot elastic limit of the target ice layer; is the hydrostatic pressure at the Hugoniot elastic limit; is the maximum tensile stress of the target ice layer; EPSO is the critical value of the quasi-static strain rate; the target ice layer is impacted by the structure, and the strain gradually increases until it reaches the deformation limit and is destroyed. The equivalent stress Reaching the maximum value is equal to the bending strength, Defined as .

[0011] The present invention also provides an ice load prediction device for a structure's icebreaking process, comprising the following modules: an acquisition module, used to acquire size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength; a calculation module, used to calculate the target stress of a target point on the target ice layer based on the small perturbation theory of an elastic thin plate according to the load to be determined, the coordinates of the current concentrated load point and the size parameters, wherein the coordinates of the current concentrated load point are the coordinates of the point where the target ice layer is subjected to the concentrated load of the structure, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; a determination module, used to determine the equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; a prediction module, used to predict the load to be determined based on the bending strength and the equivalent stress.

[0012] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, an ice load prediction method for an ice-breaking process of a structure as described above is implemented.

[0013] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method for predicting ice load in the ice-breaking process of a structure as described in any one of the above is implemented.

[0014] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the method for predicting ice load during an ice-breaking process of a structure as described in any one of the above is implemented.

[0015] The ice load prediction method, device and electronic equipment for the ice-breaking process of a structure provided by the present invention obtain the size parameters and material parameters of the target ice layer, wherein the material parameters include bending strength; calculate the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the action point of the concentrated load of the structure on the target ice layer, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; determine the equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; predict the load to be determined according to the bending strength and the equivalent stress; by calculating the relationship between the target stress of the target ice layer at the target point and the concentrated load of the structure on the target ice layer, the destruction process of ice being impacted by the structure can be more comprehensively reflected, and when the bending strength of the target ice layer is known, the ultimate load that the target ice layer can withstand at the target point can be solved. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced one by one below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 It is a flow chart of the ice load prediction method for the structure ice breaking process provided by the present invention.

[0018] Figure 2 It is one of the schematic diagrams of the ice load prediction method for the structure ice breaking process provided by the present invention.

[0019] Figure 3 This is the second schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention.

[0020] Figure 4 This is the third schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention.

[0021] Figure 5 This is the fourth schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention.

[0022] Figure 6 It is a structural schematic diagram of an ice load prediction device for a structure ice-breaking process provided by the present invention.

[0023] Figure 7 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0024] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0025] It should be noted that, in the description of the present invention, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "includes a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0026] When launching a structure in an icy ocean, the ice layer will break after reaching the limit load when the structure contacts the ice layer. This process will affect the trajectory of the structure. Therefore, the study of the interaction between the structure and ice has received extensive attention, involving a variety of nonlinear problems such as structural failure, large deformation, and crack propagation. It is of great significance to establish an accurate prediction model for the interaction between the structure and ice and to achieve rapid prediction of the ice load characteristics during the icebreaking process. In the related technology, when the structure interacts with the ice, it is impossible to accurately predict the limit load at any point of the ice layer. It can be seen that the ice load prediction method for the icebreaking process of the structure in the related technology has the technical problem of being unable to accurately predict the limit load at any point of the ice layer.

[0027] In order to at least solve some of the above problems, the following Figure 1-Figure 7 The present invention describes the ice load prediction method, device and electronic equipment for the structure ice breaking process.

[0028] The ice load prediction method for a structure during ice breaking provided in this embodiment can be applied to the scenario of predicting the ice load at any point of an ice layer during the process of ice breaking of a structure. The ice load prediction method for a structure during ice breaking in this embodiment can be executed by a server.

[0029] Figure 1 FIG. 1 is a flow chart of the ice load prediction method for the structure ice breaking process provided by the present invention, such as Figure 1 As shown, including but not limited to the following steps: Step 101, obtaining size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength.

[0030] Figure 2 is one of the schematic diagrams of the ice load prediction method for the structure ice breaking process provided by the present invention, such as Figure 2 As shown, Figure 2 Schematic diagram of the structure interacting with the target ice layer; Figure 3 FIG. 2 is a schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention. Figure 3 As shown, in this embodiment, the target ice layer is considered to be a thin plate with uniform thickness. The size parameters of the target ice layer include length, width, and thickness. Figure 2 In the coordinate system shown, the length, width and thickness of the target ice layer are a, b and h respectively, and the mid-plane of the ice layer refers to the horizontal plane at the thickness of h / 2 in the target ice layer.

[0031] The material parameters of the target ice layer are constants and are used to represent the physical properties of the target ice layer.

[0032] Step 102, calculating the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the point where the concentrated load of the structure acts on the target ice layer, and the load to be determined is the maximum concentrated load that the target point can withstand.

[0033] Optionally, based on the small perturbation theory of elastic thin plates, the target stress of the target point on the target ice layer is calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters.

[0034] The small perturbation theory of elastic thin plates, also known as the linear thin plate theory, is a theory in solid mechanics used to analyze the deformation and stress distribution of thin plates when subjected to loads perpendicular to the plate surface. The theory is based on the following assumptions: the small perturbation assumption, the perturbation (vertical displacement) of the plate is much smaller than the thickness of the plate, so the geometric nonlinear effect caused by the perturbation can be ignored; the straight normal assumption, before and after deformation, the straight line segment of the mid-plane (the middle plane when not deformed) of the plate remains a straight line, but may undergo translation and rotation; the plane stress state, along the thickness direction of the plate, the stress can be ignored, that is, the plate is considered to be subjected to only plane stress.

[0035] The load to be determined is the maximum concentrated load that the target point can bear, that is, if the load at the target point of the target ice layer exceeds the load to be determined, the target ice layer will break.

[0036] Based on the small perturbation theory of elastic thin plates, the target stress of the target point on the target ice layer can be calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters.

[0037] Step 103: determining the equivalent stress of the target ice layer according to the target stress and material parameters.

[0038] Optionally, based on the JH-2 constitutive model, the equivalent stress of the target ice layer is determined according to the target stress and material parameters.

[0039] The JH-2 constitutive model is an empirical model used to describe the dynamic mechanical behavior of materials under high pressure, high strain rate and high temperature conditions; based on the JH-2 constitutive model, the equivalent stress of the target ice layer can be determined according to the target stress and material parameters.

[0040] Step 104: predicting the load to be determined based on the bending strength and the equivalent stress.

[0041] After determining the equivalent stress, as the target ice layer is impacted by the structure, the strain gradually increases until it reaches the deformation limit and is damaged. At this time, the equivalent stress reaches its maximum value. The damage to the target ice layer caused by the impact is mainly bending damage, so the equivalent stress is equal to the bending strength at this time. Therefore, substituting the bending strength into the calculation formula of the equivalent stress can predict the value of the load to be determined.

[0042] Through the embodiments provided by the present application, the size parameters and material parameters of the target ice layer are obtained, and the material parameters include bending strength; the target stress of the target point on the target ice layer is calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the point where the target ice layer is acted upon by the concentrated load of the structure, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; the equivalent stress of the target ice layer is determined according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; the load to be determined is predicted according to the bending strength and the equivalent stress; the technical problem that the ice load prediction method of the structure breaking ice process in the related art cannot accurately predict the ultimate load of any point on the ice layer is solved, and the accuracy of ice load prediction in the structure breaking ice process is improved.

[0043] As an optional solution, based on the small perturbation theory of elastic thin plates, the target stress of the target point on the target ice layer is calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, including: Determine the disturbance according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters ; Based on the small perturbation theory of elastic thin plates, Determine the target stress: ; in, is the normal stress of the target point along the X direction, is the normal stress of the target point along the Y direction, and Together they form the target stress; is the elastic modulus; z is the distance from the target point to the middle surface of the ice layer along the Z direction; is Poisson's ratio; is the disturbance, For the disturbance The second partial derivative with respect to x is For the disturbance The second partial derivative of y, x is the horizontal coordinate of the target point, and y is the vertical coordinate of the target point.

[0044] Through this embodiment, the target stress can be accurately calculated, thereby improving the reliability and stability of the ice load prediction method during the icebreaking process of the structure.

[0045] As an optional solution, the disturbance is determined based on the load to be determined, the coordinates of the current concentrated load action point and the size parameters. ,include: Determined based on the load to be determined, the coordinates of the current concentrated load action point and the size parameters : ; Based on the simply supported boundary condition, a double trigonometric series is used to represent : ; in, is the coefficient to be determined; a, b, h are size parameters, a is the length of the target ice layer, b is the width of the target ice layer, h is the thickness of the target ice layer; F is the load to be determined; is the horizontal coordinate of the current concentrated load action point coordinates, is the ordinate of the current concentrated load action point; m and n are the number of series expansion terms; x is the abscissa of the target point, and y is the ordinate of the target point.

[0046] Through this embodiment, the disturbance can be accurately calculated , thereby improving the reliability and stability of the ice load prediction method during the structure's icebreaking process.

[0047] As an alternative, based on the small perturbation theory of elastic thin plates, Determine target stress, including: Will Substituting into the calculation formula of target stress, we get The calculation formula is as follows: ; make , then the target stress is: .

[0048] Through this embodiment, it is possible to accurately The target stress is determined, thereby improving the reliability and stability of the ice load prediction method during the icebreaking process of the structure.

[0049] As an optional solution, the equivalent stress of the target ice layer is determined based on the target stress and material parameters, including: Hydrostatic Pressure , based on the small perturbation theory of elastic thin plates, , determine the hydrostatic pressure based on the target stress: ; Based on the JH-2 constitutive model, the equivalent stress is determined according to the hydrostatic pressure and material parameters; in, is the normal stress of the target point along the Z direction; P is the hydrostatic pressure.

[0050] Through this embodiment, the equivalent stress of the target ice layer can be accurately determined according to the target stress and material parameters, thereby improving the reliability and stability of the ice load prediction method during the icebreaking process of the structure.

[0051] As an optional solution, the equivalent stress is determined based on the JH-2 constitutive model according to the hydrostatic pressure and material parameters, including: Based on the JH-2 constitutive model, normalized equivalent stress for: ; When the target ice layer is not broken, is 0, according to the material parameters: ; Determine the equivalent stress from the hydrostatic pressure: ; in, is the normalized equivalent stress of the target ice layer; is the normalized equivalent stress when the target ice layer is not destroyed; is the normalized equivalent stress when the target ice layer is completely destroyed; is the loss factor, when the target ice layer is not destroyed When the target ice layer is completely destroyed is 1; A, C, N are material parameters; is the normalized hydrostatic pressure; is the normalized maximum tensile stress of the target ice layer; is the normalized strain rate; is the equivalent stress at the Hugoniot elastic limit of the target ice layer; is the hydrostatic pressure at the Hugoniot elastic limit; is the maximum tensile stress of the target ice layer; EPSO is the critical value of the quasi-static strain rate; the target ice layer is impacted by the structure, and the strain gradually increases until it reaches the deformation limit and is destroyed. Reaching the maximum value Equal to the bending strength, Defined as .

[0052] Through this embodiment, the equivalent stress can be accurately determined according to the hydrostatic pressure and material parameters, thereby improving the reliability and stability of the ice load prediction method during the icebreaking process of the structure.

[0053] As an optional solution, based on the small perturbation theory of elastic thin plates, it can be concluded that: (1) in, is the elastic modulus; , , are the normal stresses of the target point along the x, y and z directions respectively; , is the positive strain along the x and y directions; is Poisson's ratio.

[0054] (2) in, , are the displacement components along the x and y directions respectively; is the disturbance; z is the distance from the target point to the middle surface of the ice layer along the z direction.

[0055] Substituting formula (2) into formula (1), we can obtain: (3) According to Kirchhoff thin plate theory, the elastic surface equation of ice can be expressed as: (4) in, is the load on the ice layer, D is the bending stiffness, , is the thickness of the ice layer; therefore, equation (4) can be written as: (5) The ice layer floats on the water surface, and the four sides can be considered as simply supported boundary conditions. According to formula (5), (6) Using double trigonometric series to express disturbance : (7) in, is the coefficient to be determined.

[0056] Substituting formula (7) into formula (4), we can obtain: (8) The solution is as follows: Using double trigonometric series representation : (9) Multiply both sides of equation (9) by , then From 0 to Integral; multiply by , then From 0 to integral, , For any positive integer, we can get: (10) Solving equation (10) we can obtain: (11) Combining equations (8) and (9), we get: (12) Both sides of formula (12) are trigonometric series. To maintain the equality, the coefficients must remain the same, that is: (13) Substituting formula (11) into formula (13), we can obtain: (14) When the structure impacts the target ice layer, the contact area is very small compared to the ice layer. The load on the ice layer can be approximated as a concentrated load. Assuming that the horizontal and vertical coordinates of the current concentrated load action point are , , the differential area at this position can be expressed as , then at the load acting position Can be written as , F is a concentrated load; therefore, It can be expressed as: (15) Substituting formula (15) into formula (7), we can obtain: (16) in, h is the thickness of the ice layer; (17) Substituting equations (16) and (17) into equation (3), we can obtain: (18) Formula (18) establishes the relationship between the stress generated at any point of the target ice layer and the concentrated load on the ice. However, since ice is a material with very complex mechanical properties and exhibits a high degree of nonlinearity during the destruction process, it is impossible to calculate when the ice is destroyed or obtain the maximum load that the ice can withstand simply by using Formula (18). Therefore, the relevant theories of the JH-2 constitutive model are introduced into the bending equation of the ice material.

[0057] The JH-2 constitutive model can accurately simulate the failure of materials under high strain rates. By setting strength parameters and damage parameters, the initial failure strength and failure process under different loading conditions can be determined. It is suitable for characterizing the destructive effects of ice after being impacted by structures.

[0058] The equivalent stress of the JH-2 constitutive model can be expressed as: (19) in, is the normalized equivalent stress; , are the normalized equivalent stress when the model is intact and completely destroyed, respectively. is the loss factor, which is 0 when the ice is not destroyed and 1 when the ice is completely destroyed.

[0059] (20) is the equivalent stress, is the equivalent stress at the Hugoniot elastic limit.

[0060] Under high strain rate loading, ice is a brittle material and will be destroyed when a small plastic deformation occurs. In order to simplify the calculation, the plastic deformation stage of ice is ignored and it is assumed that the destruction begins when the elastic deformation of ice reaches the elastic limit. , from formula (19), we can know that the stress generated by ice is: (twenty one) in, A , C , N is the material parameter of the target ice layer; is the normalized hydrostatic pressure; is the normalized maximum tensile stress; is the normalized strain rate.

[0061] (twenty two) (twenty three) EPSO is the critical value of quasi-static strain rate; is the equivalent plastic strain rate; is the hydrostatic pressure at the elastic limit.

[0062] (twenty four) According to the assumption of small perturbation of thin plates, , substituting formula (18) into formula (24) we can get: (25) Substituting equations (22), (23), (25) into equation (21), we can obtain: (26) The target ice layer is impacted by the structure, and the strain gradually increases until it reaches the deformation limit and is destroyed. When the ice reaches its maximum value, the damage caused by the impact is mainly bending damage, so at this time Equal to the flexural strength, defined as Combining equations (17), (20), and (26), we can obtain: (27) Formula (27) includes the size parameters and material parameters of the target ice layer, which can more comprehensively reflect the destruction process of the target ice layer under the impact of the structure. When the bending strength of the target ice layer is known, the load extreme value that the target point of the target ice layer can withstand can be solved. It should be noted that since Formula (27) contains a trigonometric series, different numbers of expansion terms need to be tried during the calculation to ensure the convergence of the final result.

[0063] Specifically, an example is given on how to calculate the ice load of the target ice layer based on the thin plate theory and the JH-2 constitutive model.

[0064] Assuming that the structure impacts the center of the ice layer along the negative direction of the z-axis, the length, width, and thickness of the ice layer are a, b, and h respectively. Determine the limit load that causes the ice layer to be destroyed.

[0065] (1) Ice material parameters are shown in Table 1: Table 1 Ice material parameters (2) According to formula (21) and (22), we can obtain It is 3438MPa.

[0066] (3) Since the structure is in vertical contact with the center of the ice layer along the negative direction of the z-axis, the coordinates of the contact point are .

[0067] (4) The ice layer usually breaks down from the contact point. The coordinates of the contact point are set as , , , Substitute the parameters obtained in steps (1), (2) and (3) into equation (27).

[0068] (5) Set the number of series expansion items m , n is 1, 10, 50, 100, 500, 1000. After programming calculation, the expressions obtained are: It can be seen that different series expansion terms, The coefficients of will be different; Figure 4 FIG. 3 is a schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention, such as Figure 4 As shown, when the number of expansion terms is 100, the coefficients are basically converged. Considering the accuracy of the calculation results and the calculation efficiency, the number of series expansion terms is finally taken as 500.

[0069] After calculation, the final result is: (28) Formula (28) can also be written as: (29) The relationship between the bending strength of the ice layer, the ultimate load that the ice plate can withstand, and the thickness of the ice layer is obtained; when the thickness and bending strength of the ice layer are known, the ultimate load that the ice layer can withstand can be obtained.

[0070] Figure 5 FIG. 4 is a schematic diagram of the ice load prediction method for the structure ice breaking process provided by the present invention, such as Figure 5 As shown in the test, the bending strength of the ice layer is The actual measured ice limit load is compared with the theoretical limit load obtained by formula (29); under the condition of ice layer thickness of 1m, the actual limit load of the ice layer is 18KN, the theoretical limit load is 19.9KN, and the difference between the theoretical value and the actual value is 10.6%. It can be seen that the results of the actual limit load and the theoretical limit load are very close, which proves that the theoretical calculation model proposed in this embodiment can better predict the ice load.

[0071] Figure 6 Schematic diagram of the structure of the ice load prediction device for the ice breaking process of the structure provided by the present invention. Figure 6 As shown, including but not limited to the following modules: An acquisition module 601 is used to acquire size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength; A calculation module 602 is used to calculate the target stress of the target point on the target ice layer based on the small perturbation theory of elastic thin plates according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the point where the concentrated load of the structure acts on the target ice layer, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; A determination module 603 is used to determine the equivalent stress of the target ice layer based on the JH-2 constitutive model according to the target stress and material parameters, where the equivalent stress is represented by a load to be determined; The prediction module 604 is used to predict the load to be determined according to the bending strength and the equivalent stress, wherein the material parameters include the bending strength.

[0072] Through the embodiments of the present application, the size parameters and material parameters of the target ice layer are obtained; based on the small perturbation theory of elastic thin plates, the target stress of the target point on the target ice layer is calculated according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the point where the target ice layer is acted upon by the concentrated load of the structure, and the load to be determined is the maximum concentrated load that the target point can withstand; based on the JH-2 constitutive model, the equivalent stress of the target ice layer is determined according to the target stress and material parameters; the load to be determined is predicted according to the bending strength and the equivalent stress, wherein the material parameters include the bending strength; the ice load prediction method for the structure breaking ice process in the related art solves the technical problem that the limit load of any point on the ice layer cannot be accurately predicted, thereby improving the accuracy of ice load prediction for the structure breaking ice process.

[0073] It should be noted that the ice load prediction device for the structure ice-breaking process provided by the present invention can execute the ice load prediction method for the structure ice-breaking process of any of the above-mentioned embodiments during specific operation, which will not be elaborated in this embodiment.

[0074] Figure 7 is a schematic diagram of the structure of the electronic device provided by the present invention, such as Figure 7As shown, the electronic device may include: a processor 710, a communication interface 720, a memory 730 and a communication bus 740, wherein the processor 710, the communication interface 720 and the memory 730 communicate with each other through the communication bus 740. The processor 710 may call the logic instructions in the memory 730 to execute the ice load prediction method for the structure breaking ice process, the method comprising: obtaining the size parameters and material parameters of the target ice layer, wherein the material parameters include bending strength; calculating the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters, wherein the coordinates of the current concentrated load action point are the coordinates of the point where the target ice layer is acted by the concentrated load of the structure, the load to be determined is the maximum concentrated load that the target point can bear, and the target stress is represented by the load to be determined; determining the equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; predicting the load to be determined according to the bending strength and the equivalent stress.

[0075] In addition, the logic instructions in the above-mentioned memory 730 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0076] On the other hand, the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the ice load prediction method for the ice breaking process of a structure provided in the above-mentioned embodiments, the method including: obtaining the size parameters and material parameters of the target ice layer, wherein the material parameters include bending strength; calculating the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load application point and the size parameters, wherein the coordinates of the current concentrated load application point are the coordinates of the point where the target ice layer is subjected to the concentrated load of the structure, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; determining the equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; predicting the load to be determined according to the bending strength and the equivalent stress.

[0077] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the ice load prediction method for the ice-breaking process of a structure provided in the above-mentioned embodiments, the method comprising: obtaining the size parameters and material parameters of the target ice layer, wherein the material parameters include bending strength; calculating the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load application point and the size parameters, wherein the coordinates of the current concentrated load application point are the coordinates of the point where the target ice layer is subjected to the concentrated load of the structure, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; determining the equivalent stress of the target ice layer according to the target stress and the material parameters, and the equivalent stress is represented by the load to be determined; and predicting the load to be determined based on the bending strength and the equivalent stress.

[0078] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, i.e., they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those of ordinary skill in the art may understand and implement it without creative effort.

[0079] Through the description of the above implementation modes, those skilled in the art can clearly understand that each implementation mode can be implemented by means of software plus a necessary general hardware platform, or of course by hardware. Based on such an understanding, the above technical solution can essentially or in other words be embodied in the form of a software product that contributes to the prior art. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, an optical disk, etc., and includes a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods of each embodiment or some parts of the embodiment.

[0080] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for predicting ice load during ice breaking of a structure, characterized in that: include: Acquiring size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength; Calculating a target stress of a target point on the target ice layer according to the load to be determined, the coordinates of a current concentrated load action point and the size parameter, wherein the coordinates of the current concentrated load action point are the coordinates of a point where the concentrated load of the structure acts on the target ice layer, the load to be determined is the maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; Determine the equivalent stress of the target ice layer according to the target stress and the material parameter, wherein the equivalent stress is represented by the load to be determined; The load to be determined is predicted according to the bending strength and the equivalent stress.

2. The ice load prediction method for a structure during ice breaking according to claim 1, characterized in that: The step of calculating the target stress of the target point on the target ice layer according to the load to be determined, the coordinates of the current concentrated load action point and the size parameter comprises: Determine the disturbance according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters ; According to the disturbance Determine the target stress: ； in, is the normal stress of the target point along the X direction, is the normal stress of the target point along the Y direction, and together constitute the target stress; is the elastic modulus; z is the distance from the target point to the middle surface of the ice layer along the Z direction; is Poisson’s ratio; is the disturbance; For the disturbance The second partial derivative with respect to x is For the disturbance The second partial derivative of y, x is the horizontal coordinate of the target point, y is the vertical coordinate of the target point.

3. The ice load prediction method for a structure during ice breaking according to claim 2, characterized in that: The disturbance is determined according to the load to be determined, the coordinates of the current concentrated load action point and the size parameter ,include: Determine according to the load to be determined, the coordinates of the current concentrated load action point and the size parameters : ； Based on simply supported boundary conditions, double trigonometric series are used to represent : ； in, is the coefficient to be determined; a, b, h are the size parameters, a is the length of the target ice layer, b is the width of the target ice layer, and h is the thickness of the target ice layer; F is the load to be determined; is the horizontal coordinate of the current concentrated load action point coordinates, is the ordinate of the current concentrated load action point; m and n are the number of series expansion items.

4. The ice load prediction method for a structure during ice breaking according to claim 3, characterized in that: The disturbance Determining the target stress includes: Will Substituting into the determination formula of the target stress, the target stress is: ； make , the target stress is: 。 5. The ice load prediction method for a structure during ice breaking according to claim 4, characterized in that: The determining the equivalent stress of the target ice layer according to the target stress and the material parameter includes: Hydrostatic Pressure , based on the small perturbation theory of elastic thin plates, , the hydrostatic pressure is determined according to the target stress: ； Based on the JH-2 constitutive model, determining the equivalent stress according to the hydrostatic pressure and the material parameters; in, is the normal stress of the target point along the Z direction; P is the hydrostatic pressure.

6. The ice load prediction method for a structure during ice breaking according to claim 5, characterized in that: The determining the equivalent stress based on the JH-2 constitutive model according to the hydrostatic pressure and the material parameters includes: Based on the JH-2 constitutive model, normalized equivalent stress for: ； When the target ice layer is not destroyed, is 0, according to the material parameters: ； The equivalent stress is determined based on the hydrostatic pressure: ； in, is the normalized equivalent stress of the target ice layer; is the normalized equivalent stress when the target ice layer is not destroyed; is the normalized equivalent stress when the target ice layer is completely destroyed; is the loss factor, when the target ice layer is not destroyed is 0, when the target ice layer is completely destroyed is 1; A, C, N are the material parameters; is the normalized hydrostatic pressure; is the normalized maximum tensile stress of the target ice layer; is the normalized strain rate; is the equivalent stress at the Hugoniot elastic limit of the target ice layer; is the hydrostatic pressure at the Hugoniot elastic limit; is the maximum tensile stress of the target ice layer; EPSO is the critical value of the quasi-static strain rate; the target ice layer is impacted by the structure, and the strain gradually increases until it reaches the deformation limit and is destroyed. The equivalent stress Reaching the maximum value is equal to the bending strength, Defined as .

7. An ice load prediction device for a structure during ice breaking, characterized in that: include: An acquisition module, used for acquiring size parameters and material parameters of a target ice layer, wherein the material parameters include bending strength; a calculation module, configured to calculate a target stress of a target point on the target ice layer according to a load to be determined, a coordinate of a current concentrated load action point, and the size parameter, wherein the coordinate of the current concentrated load action point is the coordinate of a point where the concentrated load of the structure acts on the target ice layer, the load to be determined is a maximum concentrated load that the target point can withstand, and the target stress is represented by the load to be determined; a determination module, configured to determine an equivalent stress of the target ice layer according to the target stress and the material parameter, wherein the equivalent stress is represented by the load to be determined; A prediction module is used to predict the load to be determined according to the bending strength and the equivalent stress.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method for predicting ice load during the ice-breaking process of a structure as claimed in any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for predicting ice load during the ice-breaking process of a structure as claimed in any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method for predicting ice load during the ice-breaking process of a structure as claimed in any one of claims 1 to 6 is implemented.

Images