Method and device for iteratively optimizing the crashworthiness of a marine structure under multiple load cases

Abstract

The application provides a marine structure anti-collision performance multi-working condition iterative optimization method and device, including: establishing a finite element analysis model containing soil layer information; applying operating period dead load to the finite element analysis model and performing static force solving to obtain an initial working condition of a ship collision load; obtaining multiple collision working conditions, calculating total ship collision energy and deformation control displacement; performing trial calculation on components of the finite element analysis model and setting the components as fiber beam units; calculating structure energy absorption of the multiple collision working conditions; calculating energy ratios of the multiple collision working conditions based on the structure energy absorption and the total ship collision energy, adjusting collision coefficients based on the energy ratios; calculating maximum displacements of the multiple collision working conditions, calculating displacement ratios of the multiple collision working conditions based on the maximum displacements and the deformation control displacement, adjusting diameters of rod members and wall thicknesses of the rod members based on the displacement ratios; and outputting final diameters of the rod members and wall thicknesses of the rod members. In the method, marine structure sizes can be automatically iteratively optimized, and the working efficiency of designers is improved.

Description

Technical Field

[0001] This invention relates to the field of marine structure simulation technology, and in particular to a method and apparatus for multi-condition iterative optimization of the collision resistance performance of marine structures. Background Technology

[0002] In recent years, with the rapid development of offshore wind power, the construction of marine engineering structures has also ushered in a new period of growth. As more and more marine engineering structures are built, the risk of ship collisions also increases. Therefore, avoiding accidents involving marine engineering structures caused by ship collisions can effectively reduce risks and prevent casualties and economic losses.

[0003] Existing technical solutions can use shell elements and hexahedral elements to simulate the impact force of a ship collision as static force, and nonlinear finite element analysis can effectively simulate the complex transient process from the start to the end of the collision. However, these solutions neglect material properties and plastic energy absorption, and cannot accurately simulate failure deformation under impact conditions. Furthermore, existing technical solutions can also use beam elements to construct the finite element analysis model, effectively reducing the computational load.

[0004] However, existing technical solutions still have the following drawbacks: First, there are shortcomings in terms of model realism. Existing technical solutions typically treat the lower boundary of the structure as fixed or simplified, ignoring the dynamic coupling and interaction between the pile foundation and the surrounding soil. For marine structures that rely on pile foundation bearing capacity, the constraint, damping, and energy absorption effects of the soil on the structure are crucial. Ignoring pile-soil coupling will lead to significant deviations in the prediction of collision force, structural response, and damage mode, failing to accurately reflect the actual working state of the structure.

[0005] Second, the accuracy of nonlinear material simulation needs to be improved. Using traditional beam elements or simplified constitutive models makes it difficult to accurately capture the entire process of gradual yielding, plastic hinge formation and development of component sections under impact loads, thus affecting the accurate assessment of the structure's ultimate bearing capacity and energy dissipation capacity.

[0006] Third, the design process is inefficient. Most existing technical solutions remain at the "simulation analysis-result evaluation" stage. If the results do not meet the requirements, designers need to manually adjust the structural dimensions and recalculate based on experience. This process is repetitive, time-consuming, and makes it difficult to quickly find the safe and economical optimal design solution under multiple working conditions (such as different tonnages, speeds, and impact angles), which seriously restricts design efficiency and quality.

[0007] In summary, there is an urgent need to develop a collision-resistant design method for marine structures that can accurately reflect the pile-soil coupling effect, precisely simulate nonlinear structural damage, and achieve automatic iterative optimization of dimensions under multiple working conditions, in order to fill the current technological gap and improve the safety and design intelligence level of marine structures. Summary of the Invention

[0008] In view of this, the purpose of this invention is to provide a multi-condition iterative optimization method and apparatus for the collision resistance performance of marine structures, so as to automatically iteratively optimize the dimensions of marine structures and improve the work efficiency of designers.

[0009] In a first aspect, embodiments of the present invention provide a multi-condition iterative optimization method for the collision resistance performance of marine structures. The method includes: establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports; wherein the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness; applying operational dead load to the finite element analysis model and performing static solution to obtain the initial conditions of ship collision load; wherein the operational dead load includes: gravity and equipment self-weight; acquiring multiple impact conditions and calculating the total ship collision energy and deformation control displacement of multiple impact conditions; wherein the parameters of the impact conditions include: number, initial impact load, impact coefficient, impact... Total duration and impact increment step; perform trial calculations on the components of the finite element analysis model; if the maximum stress of the target element of the component is greater than the preset threshold, set the target element as a fiber beam element; calculate the structural energy absorption under multiple impact conditions; where structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; calculate the energy ratio of multiple impact conditions based on structural energy absorption and total ship impact energy, and adjust the impact coefficient based on the energy ratio; calculate the maximum displacement under multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjust the member diameter and member wall thickness based on the displacement ratio; output the final member diameter and member wall thickness.

[0010] In an optional embodiment of this application, the above-mentioned step of adjusting the impact coefficient based on the energy ratio includes: if the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, adjusting the impact coefficient, and re-executing the steps of acquiring multiple impact conditions, calculating the total ship impact energy and deformation control displacement of the multiple impact conditions; wherein the second threshold is greater than the first threshold; if the energy ratio is greater than or equal to the first threshold and the energy ratio is less than or equal to the second threshold, performing the steps of calculating the maximum displacement of the multiple impact conditions, calculating the displacement ratio of the multiple impact conditions based on the maximum displacement and the deformation control displacement, and adjusting the rod diameter and rod wall thickness based on the displacement ratio.

[0011] In an optional embodiment of this application, the step of adjusting the impact coefficient if the energy ratio is less than a preset first threshold or greater than a preset second threshold includes: if the energy ratio is less than the preset first threshold or greater than the preset second threshold, adjusting the impact coefficient by using the reciprocal of the energy ratio as a correction coefficient.

[0012] In an optional embodiment of this application, the steps of adjusting the member diameter and member wall thickness based on the displacement ratio include: determining the maximum value from the displacement ratios of multiple impact conditions; if the maximum value of the displacement ratio is less than a preset third threshold or greater than a preset fourth threshold, adjusting the member diameter and member wall thickness, and re-executing the step of establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports; if the maximum value of the displacement ratio is greater than or equal to the third threshold and less than or equal to the fourth threshold, executing the step of outputting the final member diameter and member wall thickness.

[0013] In an optional embodiment of this application, the step of adjusting the rod diameter and rod wall thickness if the maximum value of the displacement ratio is less than a preset third threshold or greater than a preset fourth threshold includes: if the maximum value of the displacement ratio is less than a preset third threshold or greater than a preset fourth threshold, adjusting the rod diameter and rod wall thickness using the cube of the maximum value of the displacement ratio as a correction coefficient.

[0014] In an optional embodiment of this application, the steps of establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports include: classifying the soil based on marine structure drawings and geological survey reports, extracting key parameters of each soil layer and corresponding pile diameter, generating spring curves controlling the degrees of freedom of the soil layers, processing linear interpolation of variable diameter piles, and generating nonlinear curves of discrete force and displacement at each depth; fitting the nonlinear curves to continuous functions, fitting parameters for multiple degrees of freedom respectively, and recording the goodness of fit for quality control; creating a discrete element material model and defining the initial stiffness matrix of the discrete element; creating discrete elements at each soil layer depth, connecting pile nodes and fixed ground nodes, assigning nonlinear material values ​​to the discrete elements, performing nonlinear solutions, and outputting pile foundation displacement, internal force, and soil spring reaction force to obtain the finite element analysis model; wherein, while performing nonlinear solutions, the geometric nonlinearity and material nonlinearity of the pile foundation, the nonlinear force-displacement relationship of the soil spring, and the energy control criteria of ship collision loads are considered.

[0015] In an optional embodiment of this application, the step of obtaining multiple impact conditions includes: manually adding models of multiple impact conditions; and / or automatically adding models of multiple impact conditions.

[0016] In an optional embodiment of this application, the step of automatically adding multiple impact conditions to the model includes: determining the impact parameter design space of multiple impact conditions; and generating the condition matrix of multiple impact conditions using an experimental design approach.

[0017] In an optional embodiment of this application, the step of performing trial calculations on the components of the finite element analysis model and setting the target element as a fiber beam element if the maximum stress of the target element of the component is greater than a preset threshold includes: setting all elements of the components of the finite element analysis model as ordinary beam elements with linear elastic constitutive model; determining the input loads of multiple working conditions based on the initial impact load and impact coefficient; obtaining the maximum stress of each element through nonlinear finite element calculation; setting the target element as a fiber beam element if the maximum stress of the target element of the component is greater than a preset threshold; and setting all element materials of the components of the finite element analysis model as a bilinear constitutive model.

[0018] Secondly, embodiments of the present invention also provide a multi-condition iterative optimization device for the collision resistance performance of marine structures. The device includes: a finite element analysis model establishment module, used to establish a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports; wherein the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness; an initial state analysis module, used to apply operational dead load to the finite element analysis model and perform static solution to obtain the initial conditions of the ship collision load; wherein the operational dead load includes: gravity and equipment self-weight; and an impact condition and control deformation determination module, used to acquire multiple impact conditions and calculate the total ship collision energy and deformation control displacement of multiple impact conditions; wherein the parameters of the impact condition include: number, initial impact load, impact coefficient, total impact duration, and impact increment. The system includes the following modules: a fiber beam range determination module for performing trial calculations on components in the finite element analysis model; a structural energy absorption calculation module for calculating structural energy absorption under multiple impact conditions, including structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; an energy ratio calculation module for calculating the energy ratio of multiple impact conditions based on structural energy absorption and the total energy of the ship impact, and adjusting the impact coefficient based on the energy ratio; a displacement ratio calculation module for calculating the maximum displacement under multiple impact conditions, calculating the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjusting the member diameter and member wall thickness based on the displacement ratio; and a component parameter output module for outputting the final member diameter and member wall thickness.

[0019] The embodiments of the present invention bring the following beneficial effects: This invention provides a multi-condition iterative optimization method and apparatus for the collision resistance performance of marine structures. A finite element analysis model incorporating soil layer information is established based on marine structure drawings and geological survey reports. The parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness. An operational dead load is applied to the finite element analysis model, and static solutions are performed to obtain the initial conditions of the ship collision load. The operational dead load includes: gravity and equipment self-weight. Multiple impact conditions are acquired, and the total ship collision energy and deformation control displacement for each impact condition are calculated. The parameters of each impact condition include: number, initial impact load, impact coefficient, total impact duration, and... The impact increment step involves performing trial calculations on the components of the finite element analysis model. If the maximum stress of the target element of the component exceeds a preset threshold, the target element is set as a fiber beam element. The structural energy absorption under multiple impact conditions is calculated, including structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption. Based on the structural energy absorption and the total energy of the ship impact, the energy ratio of the multiple impact conditions is calculated, and the impact coefficient is adjusted based on the energy ratio. The maximum displacement under multiple impact conditions is calculated, and the displacement ratio of the multiple impact conditions is calculated based on the maximum displacement and deformation control displacement. The member diameter and member wall thickness are adjusted based on the displacement ratio. The final member diameter and member wall thickness are output. This method allows for automatic iterative optimization of marine engineering dimensions, improving the efficiency of designers.

[0020] Other features and advantages of this disclosure will be set forth in the following description, or some features and advantages may be inferred from the description or determined without doubt, or may be learned by practicing the techniques described above.

[0021] To make the above-mentioned objects, features and advantages of this disclosure more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0022] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0023] Figure 1 A flowchart of a multi-condition iterative optimization method for the collision resistance performance of marine structures provided in an embodiment of the present invention; Figure 2 A schematic diagram illustrating a multi-condition iterative optimization method for the collision resistance performance of marine structures provided in an embodiment of the present invention; Figure 3 A flowchart of another multi-condition iterative optimization method for the collision resistance performance of marine structures provided in an embodiment of the present invention; Figure 4 A schematic diagram of a marine engineering structure drawing provided for an embodiment of the present invention; Figure 5 A schematic diagram of the mechanical property curves of clay from 0.5m to 9.5m provided in an embodiment of the present invention; Figure 6 A schematic diagram of the mechanical property curves of sandy soil at a depth of 10.5m-14.5m provided for an embodiment of the present invention; Figure 7 This is a schematic diagram of a multi-condition iterative optimization device for the collision resistance performance of marine structures provided in an embodiment of the present invention; Figure 8 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0024] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0025] Currently, there is an urgent need to develop a collision-resistant design method for marine structures that can accurately reflect the pile-soil coupling effect, precisely simulate nonlinear structural damage, and achieve automatic iterative optimization of dimensions under multiple working conditions, in order to fill the current technological gap and improve the safety and intelligent design level of marine structures.

[0026] Based on this, the present invention provides a multi-condition iterative optimization method and apparatus for the impact resistance performance of marine structures. Specifically, it provides a multi-condition iterative optimization method for the impact resistance performance of marine structures based on a fiber beam-pile-soil coupling model, which can automatically iteratively optimize the dimensions of marine structures and improve the work efficiency of designers.

[0027] To facilitate understanding of this embodiment, a multi-condition iterative optimization method for the collision resistance performance of marine structures disclosed in this embodiment of the invention will first be described in detail.

[0028] Example 1: This invention provides a multi-condition iterative optimization method for the collision resistance performance of marine structures. (See also...) Figure 1 The flowchart shown illustrates a multi-condition iterative optimization method for the collision resistance performance of marine structures. This method includes the following steps: Step S102: Establish a finite element analysis model containing soil layer information based on marine structure drawings and geological survey report; wherein, the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness.

[0029] See Figure 2 This diagram illustrates a multi-condition iterative optimization method for the collision resistance performance of marine structures. This embodiment can establish a finite element analysis model. A finite element analysis model incorporating soil layer information is established based on complete marine structure drawings and geological survey reports. Parameters to be corrected include member diameter D and member wall thickness t. This embodiment can also group the finite element analysis model of the marine structure according to height to facilitate subsequent local parameter adjustments.

[0030] Step S104: Apply the operational dead load to the finite element analysis model and perform static solution to obtain the initial working condition of the ship collision load; wherein, the operational dead load includes: gravity and equipment self-weight.

[0031] like Figure 2 As shown, this embodiment can perform initial state analysis. To establish an analysis benchmark that conforms to actual working conditions, all operational dead loads, such as gravity and equipment self-weight, are applied to the established finite element analysis model, and static solutions are performed. The initial stress field and displacement field of the structure obtained from this solution are defined as the initial working condition of the ship collision load to ensure that the ship collision simulation can reflect the actual stress state of the structure before the collision.

[0032] Step S106: Obtain multiple impact conditions and calculate the total ship impact energy and deformation control displacement for multiple impact conditions; wherein the parameters of the impact conditions include: number, initial impact load, impact coefficient, total impact duration and impact increment step.

[0033] like Figure 2 As shown, this embodiment allows input of impact conditions and control of deformation. It is assumed that there are n impact conditions, with the corresponding impact locations numbered i = 1, 2, 3...n, and the initial impact load at each location is... Impact coefficient is The total impact duration is T, and the impact increment step is N; the ship weight corresponding to the i-th impact condition is... The speed at which the ship collided was Energy amplification factor The total energy of the ship collision at the corresponding impact point was calculated. The deformation control displacement for the i-th impact condition is: .

[0034] Step S108: Perform trial calculations on the components of the finite element analysis model. If the maximum stress of the target element of the component is greater than the preset threshold, set the target element as a fiber beam element.

[0035] like Figure 2 As shown, this embodiment can determine the range of fiber beams. Through trial calculations of the components, when the maximum stress of a target element is greater than a preset threshold, it can be determined that the maximum stress of the target element exceeds the yield strength by a certain proportion, and the target element can be set as a fiber beam element.

[0036] Step S110: Calculate the structural energy absorption under multiple impact conditions; wherein, structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption.

[0037] like Figure 2 As shown, this embodiment can calculate structural energy absorption. Structural energy absorption. ,in This indicates that structural deformation absorbs energy. This indicates that the component is concave and absorbs energy. This indicates that the ship deforms and absorbs energy, and the subscript i indicates the i-th impact condition.

[0038] This embodiment can be obtained through nonlinear finite element analysis. The calculation formula is: ,in This indicates that the structure resists the force generated by deformation. Represents the velocity of the node. This is determined by the impact over the total impact duration t. The energy absorbed by structural deformation can be obtained by integrating the integral.

[0039] This embodiment can be obtained through nonlinear finite element analysis. The calculation formula is: ,in The energy consumed by the component through the depression. This refers to the force that creates a dent (i.e., the impact force). The indentation depth is indicated by the relationship between the indentation force and the indentation depth of tubular structures such as CPEllinas (a type of tubular support member used in offshore facilities such as offshore platforms). ,in, Let be the material's yield stress, D be the pipe diameter, and t be the pipe wall thickness. This is determined by the impact dent depth. Internal The integral yields the energy absorption of the component's concavity.

[0040] This embodiment can be obtained through nonlinear finite element analysis. The calculation formula is: Where Rs is the load-deformation resistance function of the ship, characterizing the collision resistance generated by the ship during deformation. This represents the maximum deformation of the ship during a collision. The ship's load-deformation curve should be calculated by integration based on different operating conditions.

[0041] Step S112: Calculate the energy ratio of multiple impact scenarios based on structural energy absorption and total ship collision energy, and adjust the impact coefficient based on the energy ratio.

[0042] This embodiment can compare the total energy of ship collisions under various impact scenarios. and structural energy absorption Calculate the energy ratio .

[0043] In some embodiments, if the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, the impact coefficient is adjusted, and the steps of acquiring multiple impact conditions, calculating the total ship impact energy and deformation control displacement of the multiple impact conditions are re-executed; wherein, the second threshold is greater than the first threshold; if the energy ratio is greater than or equal to the first threshold and the energy ratio is less than or equal to the second threshold, the steps of calculating the maximum displacement of the multiple impact conditions, calculating the displacement ratio of the multiple impact conditions based on the maximum displacement and deformation control displacement, and adjusting the rod diameter and rod wall thickness based on the displacement ratio are executed.

[0044] like Figure 2 As shown, in this embodiment, the first threshold can be set to 1, and the second threshold can be set to 1.05. It can be determined whether... If yes, proceed to the next steps; otherwise, return to the steps for inputting the impact condition and controlling deformation, and modify the impact coefficient again. .

[0045] This embodiment can be achieved through energy comparison. Corrections are made to ensure that the energy absorbed by the components due to the impact and the total energy of the ship collision are within a reasonable range, providing a reasonable initial value for subsequent component corrections.

[0046] It should also be noted that when multiple impact conditions exist, if several impact conditions satisfy... or At that time, the other impact conditions meet the requirements. However, not all impact conditions require a step to return to the input impact condition and control deformation, and then make modifications. Only satisfy or The working condition returns to the input of the impact condition and the steps for controlling deformation and modification. .

[0047] In some embodiments, if the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, the reciprocal of the energy ratio is used as a correction coefficient to adjust the impact coefficient.

[0048] In this embodiment, it can be The impact coefficient is adjusted as a correction factor. .

[0049] Step S114: Calculate the maximum displacement of multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjust the rod diameter and rod wall thickness based on the displacement ratio.

[0050] This embodiment can calculate the displacement ratio: calculate the maximum displacement under various impact conditions. Through maximum displacement and deformation control displacement By comparison, the displacement ratio is obtained. .

[0051] In some embodiments, the maximum value can be determined from the displacement ratios of multiple impact conditions; if the maximum value of the displacement ratio is less than a preset third threshold or greater than a preset fourth threshold, the member diameter and member wall thickness are adjusted, and the step of establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports is re-executed; if the maximum value of the displacement ratio is greater than or equal to the third threshold and less than or equal to the fourth threshold, the step of outputting the final member diameter and member wall thickness is executed.

[0052] This embodiment can obtain maximum value ,like Figure 2 As shown, in this embodiment, the third threshold can be set to 0.95, and the fourth threshold can be set to 1. It can be determined whether... If yes, proceed to the next steps; otherwise, return to the steps of establishing the finite element analysis model and modify the member diameter D and member wall thickness t.

[0053] In some embodiments, if the maximum value of the displacement ratio is less than a preset third threshold or the maximum value of the displacement ratio is greater than a preset fourth threshold, the cube of the maximum value of the displacement ratio is used as a correction coefficient to adjust the rod diameter and rod wall thickness.

[0054] Considering that the deformation caused by the ship collision is a lateral displacement, the plastic section modulus of the circular tube cross-section... The calculation formula is: When unfolded, it becomes In the formula, D is the outer diameter of the circular tube, and t is the wall thickness, both of which are in the dimension of length. Therefore, the terms on the right side of the formula ( Since all values ​​are cubes of length, the dimension of the calculated result Z is also cubes of length. Therefore, the correction factors for the member diameter D and member wall thickness t can be taken as... .

[0055] Step S116: Output the final rod diameter and rod wall thickness.

[0056] like Figure 2As shown, this embodiment can output component parameters, including the rod diameter D and the rod wall thickness t. The rod diameter D is rounded to the nearest centimeter, and the rod wall thickness t is rounded to the nearest multiple of 2 mm.

[0057] Considering that in actual engineering projects, the diameter D of the member is usually measured in centimeters, and the wall thickness t is generally a multiple of 2 mm, to ensure engineering safety, this embodiment can round up the diameter D to the nearest centimeter (e.g., 109.4 cm is rounded up to 110 cm). Since the diameter already has a safety margin, the wall thickness calculated by t should be rounded up to the nearest multiple of 2 mm (e.g., 18.7 mm is rounded up to 18 mm, 19.1 mm is rounded up to 20 mm).

[0058] This invention provides a multi-condition iterative optimization method for the collision resistance performance of marine structures. A finite element analysis model incorporating soil layer information is established based on marine structure drawings and geological survey reports. The parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness. An operational dead load is applied to the finite element analysis model, and static solutions are performed to obtain the initial conditions of the ship collision load. The operational dead load includes: gravity and equipment self-weight. Multiple impact conditions are acquired, and the total ship collision energy and deformation control displacement for each impact condition are calculated. The parameters of each impact condition include: number, initial impact load, impact coefficient, total impact duration, and impact... The process involves incrementing the step size; performing trial calculations on the components of the finite element analysis model; if the maximum stress of the target element of the component exceeds a preset threshold, setting the target element as a fiber beam element; calculating the structural energy absorption under multiple impact conditions; where structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; calculating the energy ratio of multiple impact conditions based on structural energy absorption and the total energy of the ship impact, and adjusting the impact coefficient based on the energy ratio; calculating the maximum displacement under multiple impact conditions, calculating the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjusting the member diameter and member wall thickness based on the displacement ratio; and outputting the final member diameter and member wall thickness. This method allows for automatic iterative optimization of marine engineering dimensions, improving the work efficiency of designers.

[0059] Example 2: This invention provides another method for iterative optimization of the collision resistance performance of marine structures under multiple working conditions, see [link to relevant documentation]. Figure 3 The flowchart shown represents another method for iterative optimization of the collision resistance performance of marine structures under multiple working conditions. This method includes the following steps: Step S302: Establish a finite element analysis model containing soil layer information based on marine structure drawings and geological survey report; wherein, the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness.

[0060] In some embodiments, a finite element analysis model can be established through the following steps A1-A4, which can be referred to as a pile-soil coupling model: Step A1 can classify the soil based on marine engineering structure drawings and geological survey reports, extract key parameters of each soil layer and corresponding pile diameter, generate spring curves that control the degree of freedom of soil layers, process linear interpolation of variable diameter piles, and generate nonlinear curves of discrete force and displacement at each depth. Step A2: Fit the nonlinear curve to a continuous function, fit the parameters for multiple degrees of freedom, and record the goodness of fit for quality control. In this embodiment, the fit can be a Ramberg-Osgood (a classic theoretical model in solid mechanics used to describe the nonlinear stress-strain behavior of materials near the yield point) continuous function, with parameters fitted for each of the six degrees of freedom (DX, DY, DZ, RX, RY, RZ).

[0061] Step A3: Create a discrete element material model and define the initial stiffness matrix of the discrete element; Step A4: Create discrete elements at each soil depth, connect pile nodes to fixed ground nodes, assign nonlinear material values ​​to the discrete elements, perform nonlinear solutions, and output pile foundation displacement, internal forces, and soil spring reactions to obtain the finite element analysis model. During the nonlinear solution, the geometric nonlinearity and material nonlinearity of the pile foundation, the nonlinear force-displacement relationship of the soil spring, and the energy control criteria of the ship collision load are considered.

[0062] In this embodiment, the Newton-Raphson iteration can be used for nonlinear solution, while considering the geometric nonlinearity and material nonlinearity of the pile foundation, the nonlinear force-displacement relationship of the soil spring, and the energy control criterion of the ship collision load. Finally, the pile foundation displacement, internal force, and soil spring reaction force are output.

[0063] Step S304: Apply the operational dead load to the finite element analysis model and perform static solution to obtain the initial working condition of the ship collision load; wherein, the operational dead load includes: gravity and equipment self-weight.

[0064] Step S306: Manually add models for multiple impact scenarios; and / or automatically add models for multiple impact scenarios. Calculate the total ship impact energy and deformation control displacement for multiple impact scenarios; wherein the parameters of the impact scenarios include: number, initial impact load, impact coefficient, total impact duration, and impact increment step.

[0065] Assuming there are n impact conditions, the impact conditions in this embodiment can be applied manually and / or automatically.

[0066] In some embodiments, the impact parameter design space for multiple impact conditions can be determined; and the condition matrix for multiple impact conditions can be generated by experimental design.

[0067] If an automatic model addition method is used, this embodiment can determine the impact parameter design space: based on standard recommendations and project-specific risk analysis, the variation range and typical values ​​of key parameters such as ship tonnage (DWT), impact velocity (V), and impact angle (θ) are determined. This embodiment can also generate a test case matrix using experimental design methods: within the determined parameter design space, Latin hypercube sampling (LHS) is used to systematically generate a set (N) of representative impact case combinations. This method ensures that the generated impact cases uniformly cover the entire design space, avoiding the bias of manual selection.

[0068] Step S308: Perform trial calculations on the components of the finite element analysis model. If the maximum stress of the target element of the component is greater than the preset threshold, set the target element as a fiber beam element.

[0069] In some embodiments, the elements of the components in the finite element analysis model can be set as ordinary beam elements with linear elastic constitutive model; the input loads for multiple working conditions are determined based on the initial impact load and impact coefficient; the maximum stress of each element is obtained through nonlinear finite element calculation; if the maximum stress of the target element of the component is greater than a preset threshold, the target element is set as a fiber beam element; and the element material of the components in the finite element analysis model is set as a bi-segmented constitutive model.

[0070] In finite element analysis, not all elements enter the plastic state. Setting all elements as fiber beam elements would waste computational resources. This embodiment allows for trial calculations on components to determine the fiber beam range. For example, all elements can be set as ordinary beam elements using a linear elastic constitutive model. Initial impact loads can be set and / or generated. Impact coefficient The input loads for each impact condition are obtained. The maximum stress σ of each element is obtained through nonlinear finite element calculation. When σ exceeds the set value, the element is replaced with a fiber beam element. All element materials are then set to a bilinear constitutive model and incorporated into subsequent calculations.

[0071] Step S310: Calculate the structural energy absorption under multiple impact conditions; whereby structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption.

[0072] Step S312: Calculate the energy ratio of multiple impact scenarios based on structural energy absorption and total ship collision energy, and adjust the impact coefficient based on the energy ratio.

[0073] Step S314: Calculate the maximum displacement of multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjust the rod diameter and rod wall thickness based on the displacement ratio.

[0074] Step S316: Output the final rod diameter and rod wall thickness.

[0075] The method provided in the embodiments of the present invention mainly provides the following: 1. Through continuous exploration and experimentation, the introduction of nonlinear soil spring elements to accurately simulate pile-soil interaction effectively overcomes the errors caused by traditional fixed connections.

[0076] 2. Fiber beam elements are used instead of traditional beam or shell elements to simulate the main load-bearing components. Based on the material constitutive relationship, this element can automatically calculate the internal force distribution and deformation of the section under complex stress, thereby simulating the entire process of the component's development from elasticity and yielding to plasticity with high accuracy, while maintaining high computational efficiency.

[0077] In summary, the method provided by the embodiments of the present invention can construct a highly efficient closed-loop optimization system in terms of analysis process and efficiency by integrating parametric modeling, automated processes, and intelligent optimization algorithms: it automatically executes parallel calculations for multiple working conditions, completing the full analysis of multiple working conditions at once, significantly shortening the calculation cycle; it intelligently identifies the most unfavorable working condition, completely avoiding the tedious process and potential omissions of relying on manual judgment item by item in traditional methods, ensuring the objectivity and reliability of the evaluation benchmark; and it uses the most unfavorable working condition as the target working condition for optimization iteration, ensuring that the optimized design scheme can meet safety requirements in all possible collision scenarios, thereby achieving automatic optimization of component dimensions in a more economical and reasonable direction while ensuring the overall structural safety.

[0078] Example 3: This invention provides a specific example of a multi-condition iterative optimization method for the collision resistance performance of marine structures, which can be used to perform simulation analysis on a certain marine structure. The method mainly includes the following steps S1-S8.

[0079] Step S1, see Figure 4 The diagram shows a schematic of an offshore engineering structure. A finite element model was built based on the drawings. The upper structure of the model uses a single circular pipe, 10m long, with a pile diameter of 100cm and a wall thickness of 2.5cm, divided into 10 units, each 1m long. The lower structure also uses a single circular pipe, with a pile length of 15m, a pile diameter of 310cm, and a wall thickness of 9cm. The impact points are located at 003L, 005L, and 008L.

[0080] The superstructure and pile foundation are based on the beam element model constructed by Timoshenko beam and the von Mises linear follower hardening plastic constitutive model. The soil spring stiffness of each soil layer is simulated by load-displacement curves, thereby establishing the coupling relationship between the soil structure and the pile structure.

[0081] The relationship between the stiffness and degrees of freedom of each soil layer is as follows: taking the pile axis as the local coordinate X-axis, the Y and Z directions of the right-hand rule represent the horizontal direction of the pile, respectively. Therefore, the PY curve simulates the horizontal constraint of the pile (DY, DZ), the TZ curve simulates the axial constraint of the pile (DX), the QZ curve simulates the axial constraint of the pile bottom (DX), the MT (torsion) curve simulates the rotation constraint around the pile (RX), and the MT (BENDING) curve simulates the bending constraint around the pile (RY, RZ), etc.

[0082] Taking an ideal soil as an example, it consists of two layers: 10m thick clay and 5m thick sand. The main parameters of the clay are as follows: thickness is 10m, and buoyancy is 8.8kN / m³. 3 The comprehensive recommendations are as follows: undrained shear strength of 80 kPa, strain at half the ultimate strength of unconsolidated undrained compression test of clay of 0.07, standard value of ultimate skin friction of pile side of 75 kPa, pull-out coefficient of 0.54, and dimensionless support capacity coefficient of 20.

[0083] The main parameters of the sand are as follows: thickness is 5m, and buoyant unit weight is 9.93kN / m³. 3 The internal friction angle is 32°, and the initial response modulus of the foundation is 15400 kN / m. 3 The standard value of the ultimate skin friction of the pile side is 75 kPa, the standard value of the ultimate end resistance of the pile tip is 3900 kPa, the pull-out coefficient is 0.54, and the dimensionless support capacity coefficient is 20.

[0084] The spring, adjusted based on pile diameter and soil parameters, can be found in [reference needed]. Figure 5 A schematic diagram showing the mechanical property curves of clay from 0.5m to 9.5m is provided. Figure 6 The diagram shows a schematic of the mechanical property curves of sandy soil at depths of 10.5m-14.5m.

[0085] In this embodiment, QZ is the pile bottom constant set according to soil parameters, and the MT (BENDING) curve is a non-essential item, so neither is illustrated.

[0086] Step S2: Initialize the model by applying loads such as self-weight and buoyancy to the model.

[0087] Step S3: Apply three independent impact conditions to the model, as detailed in Table 1 below: Table 1

[0088] Step S4: Perform component trial calculations. In this case, the maximum stress of the member is 0.6 times the yield strength as the criterion for whether to use fiber beam elements. Through calculation, fiber beam elements are used for members numbered (000L~006L) in working condition 1; fiber beam elements are used for members numbered (000L~009L) in working condition 2; and fiber beam elements are used for members numbered (000L~004L) in working condition 3.

[0089] Step S5, calculate structural energy absorption, structural energy absorption under load condition 1. Working condition 2: Structural energy absorption Working condition 3: Structural energy absorption .

[0090] Step S6, calculate the structural energy ratio, for example: ; ; .

[0091] Among them, satisfying ,right Make a correction with a correction factor of 1.630, and return to step S3. Make corrections, after corrections Among them, those that satisfy the condition when The current load condition is retained without modification. (This condition is satisfied.) ,right Make a correction with a correction factor of 0.730, and return to step S3. Make corrections, after corrections .

[0092] Return to step S4 and recalculate the components for working conditions 1 and 3. Use the maximum stress of the member as 0.6 times the yield strength as the criterion for whether to use fiber beam elements. Through calculation, fiber beam elements are used for members numbered (000L~0006L) in working condition 1 and for members numbered (000L~0004L) in working condition 3.

[0093] Re-enter step S5 to calculate structural energy absorption; structural energy absorption under load case 1. MJ; Situation 3 structure energy absorption .

[0094] Re-enter step S6 to calculate the structural energy ratio: .

[0095] Among them, satisfying when The current load condition will be retained without modification.

[0096] Among them, satisfying when The current load condition will be retained without modification.

[0097] all All meet Proceed to step S7.

[0098] Step S7, proceed to calculate the displacement ratio: ; ; 6; 1.116.

[0099] Through calculation, it was found that The diameter D and wall thickness t of the member are corrected by a correction factor of . Repeat steps S1-S7 to obtain the results. The conditions are met. Enter S8.

[0100] Step S8: Output the diameter D and wall thickness t of each member. The diameter D is rounded to the nearest centimeter, and the wall thickness t is rounded to the nearest multiple of 2 mm.

[0101] This embodiment uses three random working conditions as general working conditions to verify the members at the impact points. The resulting member diameters are only the envelope results for three specific cases. Furthermore, this embodiment is only a single-pile structure; to maintain component consistency, the maximum value of the envelope result, i.e., D=110cm, t=28mm, can be used, as shown in Table 2. Subsequently, different ship parameters, different impact points, and other working conditions can be designed according to the complexity of the model, and the component dimensions can be optimized for different components.

[0102] Table 2

[0103] Example 4: Corresponding to the above method embodiments, this invention provides a multi-condition iterative optimization device for the collision resistance performance of marine structures. (See also...) Figure 7 The diagram shows a multi-condition iterative optimization device for the collision resistance performance of marine structures. This device includes: The finite element analysis model establishment module 71 is used to establish a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports; among which, the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness; The initial state analysis module 72 is used to apply the operational dead load to the finite element analysis model and perform static solution to obtain the initial working condition of the ship collision load; wherein, the operational dead load includes: gravity and equipment self-weight; The impact condition and control deformation determination module 73 is used to acquire multiple impact conditions and calculate the total ship impact energy and deformation control displacement of multiple impact conditions; the parameters of the impact condition include: number, initial impact load, impact coefficient, total impact duration and impact increment step; The fiber beam range determination module 74 is used to perform trial calculations on the components of the finite element analysis model. If the maximum stress of the target element of the component is greater than the preset threshold, the target element is set as a fiber beam element. The structural energy absorption calculation module 75 is used to calculate the structural energy absorption under multiple impact conditions; among which, structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption. The energy ratio calculation module 76 is used to calculate the energy ratio of multiple impact scenarios based on structural energy absorption and total ship collision energy, and to adjust the impact coefficient based on the energy ratio. The displacement ratio calculation module 77 is used to calculate the maximum displacement under multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjust the rod diameter and rod wall thickness based on the displacement ratio. The component parameter output module 78 is used to output the final member diameter and member wall thickness.

[0104] This invention provides a multi-condition iterative optimization device for the collision resistance performance of marine structures. A finite element analysis model incorporating soil layer information is established based on marine structure drawings and geological survey reports. The parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness. An operational dead load is applied to the finite element analysis model, and static solutions are performed to obtain the initial conditions of the ship collision load. The operational dead load includes: gravity and equipment self-weight. Multiple impact conditions are acquired, and the total ship collision energy and deformation control displacement for each impact condition are calculated. The parameters of each impact condition include: number, initial impact load, impact coefficient, total impact duration, and impact... The process involves incrementing the step size; performing trial calculations on the components of the finite element analysis model; if the maximum stress of the target element of the component exceeds a preset threshold, setting the target element as a fiber beam element; calculating the structural energy absorption under multiple impact conditions; where structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; calculating the energy ratio of multiple impact conditions based on structural energy absorption and the total energy of the ship impact, and adjusting the impact coefficient based on the energy ratio; calculating the maximum displacement under multiple impact conditions, calculating the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjusting the member diameter and member wall thickness based on the displacement ratio; and outputting the final member diameter and member wall thickness. This method allows for automatic iterative optimization of marine engineering dimensions, improving the work efficiency of designers.

[0105] The aforementioned energy ratio calculation module is used to adjust the impact coefficient and re-execute the steps of acquiring multiple impact conditions, calculating the total ship impact energy and deformation control displacement of multiple impact conditions based on whether the energy ratio is less than a preset first threshold or greater than a preset second threshold; wherein, the second threshold is greater than the first threshold; if the energy ratio is greater than or equal to the first threshold and less than or equal to the second threshold, the module executes the steps of calculating the maximum displacement of multiple impact conditions, calculating the displacement ratio of multiple impact conditions based on the maximum displacement and deformation control displacement, and adjusting the rod diameter and rod wall thickness based on the displacement ratio.

[0106] The aforementioned energy ratio calculation module is used to adjust the impact coefficient by using the reciprocal of the energy ratio as a correction coefficient if the energy ratio is less than a preset first threshold or greater than a preset second threshold.

[0107] The displacement ratio calculation module described above is used to determine the maximum value from the displacement ratios of multiple impact conditions. If the maximum value of the displacement ratio is less than a preset third threshold or greater than a preset fourth threshold, the member diameter and member wall thickness are adjusted, and the step of establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports is re-executed. If the maximum value of the displacement ratio is greater than or equal to the third threshold and less than or equal to the fourth threshold, the step of outputting the final member diameter and member wall thickness is executed.

[0108] The aforementioned displacement ratio calculation module is used to adjust the rod diameter and rod wall thickness by using the cube of the maximum displacement ratio as a correction coefficient if the maximum displacement ratio is less than a preset third threshold or greater than a preset fourth threshold.

[0109] The aforementioned finite element analysis model building module is used to classify the soil based on marine structural drawings and geological survey reports, extract key parameters of each soil layer and corresponding pile diameter, generate spring curves controlling the degrees of freedom of the soil layers, process linear interpolation of variable diameter piles, and generate nonlinear curves of discrete force and displacement at each depth. The nonlinear curves are fitted to continuous functions, parameters are fitted for multiple degrees of freedom, and the goodness of fit is recorded for quality control. A discrete element material model is created, and the initial stiffness matrix of the discrete element is defined. Discrete elements are created at each soil layer depth, connecting pile nodes and fixed ground nodes. Nonlinear material values ​​are assigned to the discrete elements, and nonlinear solutions are performed to output the pile foundation displacement, internal forces, and soil spring reactions, thus obtaining the finite element analysis model. During the nonlinear solution process, the geometric nonlinearity and material nonlinearity of the pile foundation, the nonlinear force-displacement relationship of the soil springs, and the energy control criteria for ship collision loads are considered.

[0110] The aforementioned impact condition and control deformation determination module is used to manually add models of multiple impact conditions; and / or automatically add models of multiple impact conditions.

[0111] The aforementioned impact condition and control deformation determination module is used to determine the impact parameter design space for multiple impact conditions; and to generate the condition matrix for multiple impact conditions using an experimental design approach.

[0112] The aforementioned fiber beam range determination module is used to set all the elements of the components in the finite element analysis model to ordinary beam elements with linear elastic constitutive model; determine the input loads for multiple working conditions based on the initial impact load and impact coefficient; obtain the maximum stress of each element through nonlinear finite element calculation; if the maximum stress of the target element of the component is greater than the preset threshold, set the target element to a fiber beam element; and set the element material of the components in the finite element analysis model to a bi-segmented constitutive model.

[0113] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the multi-condition iterative optimization device for the collision resistance performance of marine structures described above can be referred to the corresponding process in the aforementioned embodiments of the multi-condition iterative optimization method for the collision resistance performance of marine structures, and will not be repeated here.

[0114] Example 4: This invention also provides an electronic device for running the above-described multi-condition iterative optimization method for the collision resistance performance of marine structures; see also Figure 8 The diagram shows the structure of an electronic device, which includes a memory 100 and a processor 101. The memory 100 stores one or more computer instructions, which are executed by the processor 101 to implement the above-mentioned multi-condition iterative optimization method for the collision resistance performance of marine structures.

[0115] Furthermore, Figure 8 The electronic device shown also includes a bus 102 and a communication interface 103, with the processor 101, the communication interface 103 and the memory 100 connected via the bus 102.

[0116] The memory 100 may include high-speed random access memory (RAM) and may also include non-volatile memory, such as at least one disk storage device. Communication between this system network element and at least one other network element is achieved through at least one communication interface 103 (which can be wired or wireless), such as the Internet, wide area network, local area network, metropolitan area network, etc. The bus 102 may be an ISA bus, PCI bus, or EISA bus, etc. The bus can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 8 The symbol is represented by a single double-headed arrow, but this does not mean that there is only one bus or one type of bus.

[0117] Processor 101 may be an integrated circuit chip with signal processing capabilities. In implementation, each step of the above method can be completed by the integrated logic circuitry in the hardware of processor 101 or by instructions in software form. Processor 101 can be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), etc.; it can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. It can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of this invention. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the methods disclosed in the embodiments of this invention can be directly manifested as execution by a hardware decoding processor, or execution by a combination of hardware and software modules in the decoding processor. The software module can reside in a readily available storage medium in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, or registers. This storage medium is located in memory 100, and processor 101 reads information from memory 100 and, in conjunction with its hardware, completes the steps of the method described in the foregoing embodiments.

[0118] This invention also provides a computer-readable storage medium storing computer-executable instructions. When these computer-executable instructions are called and executed by a processor, they cause the processor to implement the above-described multi-condition iterative optimization method for the collision resistance performance of marine structures. For specific implementation details, please refer to the method embodiments, which will not be repeated here.

[0119] The computer program product of the multi-condition iterative optimization method and apparatus for the collision resistance performance of marine structures provided in this embodiment of the invention includes a computer-readable storage medium storing program code. The instructions included in the program code can be used to execute the methods in the preceding method embodiments. For specific implementation, please refer to the method embodiments, which will not be repeated here.

[0120] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the system and / or device described above can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0121] Furthermore, in the description of the embodiments of the present invention, unless otherwise explicitly specified and limited, the terms "installation," "connection," and "linking" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal connection of two components. Those skilled in the art can understand the specific meaning of the above terms in the present invention based on the specific circumstances.

[0122] If a function is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0123] In the description of this invention, it should be noted that the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," and "outer," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are used only for the convenience of describing the invention and for simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on the invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0124] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A multi-condition iterative optimization method for the collision resistance performance of marine structures, characterized in that, The method includes: A finite element analysis model containing soil layer information is established based on marine structure drawings and geological survey reports; wherein, the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness; The operational dead load is applied to the finite element analysis model and static solution is performed to obtain the initial working condition of the ship collision load; wherein, the operational dead load includes: gravity and equipment self-weight; Multiple impact scenarios are acquired, and the total ship impact energy and deformation control displacement for each impact scenario are calculated. The parameters of each impact scenario include: number, initial impact load, impact coefficient, total impact duration, and impact increment step. The component of the finite element analysis model is tested. If the maximum stress of the target element of the component is greater than a preset threshold, the target element is set as a fiber beam element. Calculate the structural energy absorption under multiple impact conditions; wherein, the structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; Based on the energy absorption of the structure and the total energy of the ship collision, the energy ratio of multiple impact scenarios is calculated, and the impact coefficient is adjusted based on the energy ratio. Calculate the maximum displacement of multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and the deformation control displacement, and adjust the rod diameter and rod wall thickness based on the displacement ratio; Output the final diameter and wall thickness of the rod.

2. The method according to claim 1, characterized in that, The step of adjusting the impact coefficient based on the energy ratio includes: If the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, the impact coefficient is adjusted, and the steps of acquiring multiple impact conditions and calculating the total ship impact energy and deformation control displacement of the multiple impact conditions are re-executed; wherein, the second threshold is greater than the first threshold; If the energy ratio is greater than or equal to the first threshold and the energy ratio is less than or equal to the second threshold, perform the steps of calculating the maximum displacement of multiple impact conditions, calculating the displacement ratio of multiple impact conditions based on the maximum displacement and the deformation control displacement, and adjusting the rod diameter and the rod wall thickness based on the displacement ratio.

3. The method according to claim 2, characterized in that, If the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, the step of adjusting the impact coefficient includes: If the energy ratio is less than a preset first threshold or the energy ratio is greater than a preset second threshold, the reciprocal of the energy ratio is used as a correction coefficient to adjust the impact coefficient.

4. The method according to claim 1, characterized in that, The steps of adjusting the diameter and wall thickness of the rod based on the displacement ratio include: The maximum value is determined from the displacement ratios of the multiple impact conditions; If the maximum value of the displacement ratio is less than the preset third threshold or the maximum value of the displacement ratio is greater than the preset fourth threshold, adjust the diameter of the rod and the wall thickness of the rod, and re-execute the step of establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports. If the maximum value of the displacement ratio is greater than or equal to the third threshold and the maximum value of the displacement ratio is less than or equal to the fourth threshold, then the step of outputting the final rod diameter and rod wall thickness is performed.

5. The method according to claim 4, characterized in that, If the maximum value of the displacement ratio is less than a preset third threshold or the maximum value of the displacement ratio is greater than a preset fourth threshold, the step of adjusting the diameter of the rod and the wall thickness of the rod includes: If the maximum value of the displacement ratio is less than a preset third threshold or the maximum value of the displacement ratio is greater than a preset fourth threshold, the cube of the maximum value of the displacement ratio is used as a correction coefficient to adjust the diameter of the rod and the wall thickness of the rod.

6. The method according to claim 1, characterized in that, The steps for establishing a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports include: Based on marine engineering structure drawings and geological survey reports, the soil is classified, key parameters of each soil layer and corresponding pile diameter are extracted, spring curves controlling the degree of freedom of soil layers are generated, linear interpolation of variable diameter piles is processed, and nonlinear curves of discrete force and displacement at each depth are generated. The nonlinear curve is fitted to a continuous function, and parameters are fitted for multiple degrees of freedom. The goodness of fit is recorded for quality control. Create a discrete element material model and define the initial stiffness matrix of the discrete element; Discrete elements are created at each soil depth, and pile nodes are connected to fixed ground nodes. Nonlinear material values ​​are assigned to the discrete elements, and nonlinear solutions are performed to output the pile foundation displacement, internal forces, and soil spring reactions, thus obtaining the finite element analysis model. The nonlinear solution considers the geometric nonlinearity and material nonlinearity of the pile foundation, the nonlinear force-displacement relationship of the soil spring, and the energy control criteria of the ship collision load.

7. The method according to claim 1, characterized in that, The steps to obtain multiple impact conditions include: Manually add models with multiple impact scenarios; And / or, automatically add multiple models of the aforementioned impact conditions.

8. The method according to claim 7, characterized in that, The steps for automatically adding multiple models of the aforementioned impact conditions include: Determine the impact parameter design space for multiple impact conditions; Multiple impact condition matrices were generated using experimental design.

9. The method according to claim 1, characterized in that, The step of performing trial calculations on the components of the finite element analysis model, and setting the target element as a fiber beam element if the maximum stress of the target element of the component is greater than a preset threshold, includes: The elements of the components in the finite element analysis model are all set as ordinary beam elements with linear elastic constitutive structure; The input loads for multiple operating conditions are determined based on the initial impact load and the impact coefficient. The maximum stress of each element is obtained through nonlinear finite element calculation; If the maximum stress of the target unit of the component is greater than a preset threshold, the target unit is set as a fiber beam unit; The element materials of the components in the finite element analysis model are all set as bilinear constitutive models.

10. A multi-condition iterative optimization device for the collision resistance performance of marine structures, characterized in that, The device includes: The finite element analysis model building module is used to build a finite element analysis model containing soil layer information based on marine structure drawings and geological survey reports; wherein, the parameters to be corrected in the finite element analysis model include: member diameter and member wall thickness; The initial state analysis module is used to apply the operational dead load to the finite element analysis model and perform static solution to obtain the initial working condition of the ship collision load; wherein, the operational dead load includes: gravity and equipment self-weight; The impact condition and control deformation determination module is used to acquire multiple impact conditions and calculate the total ship impact energy and deformation control displacement of the multiple impact conditions; wherein the parameters of the impact condition include: number, initial impact load, impact coefficient, total impact duration and impact increment step; The fiber beam range determination module is used to perform trial calculations on the components of the finite element analysis model. If the maximum stress of the target element of the component is greater than a preset threshold, the target element is set as a fiber beam element. The structural energy absorption calculation module is used to calculate the structural energy absorption of multiple impact conditions; wherein, the structural energy absorption includes: structural deformation energy absorption, component indentation energy absorption, and ship deformation energy absorption; An energy ratio calculation module is used to calculate the energy ratio of multiple impact scenarios based on the energy absorbed by the structure and the total energy of the ship collision, and to adjust the impact coefficient based on the energy ratio. The displacement ratio calculation module is used to calculate the maximum displacement of multiple impact conditions, calculate the displacement ratio of multiple impact conditions based on the maximum displacement and the deformation control displacement, and adjust the rod diameter and rod wall thickness based on the displacement ratio. The component parameter output module is used to output the final diameter and wall thickness of the rod.

Images