Method, device, equipment and medium for optimizing non-bonded structure of flexible riser for deep sea mining

Abstract

The application discloses a flexible riser non-bonding structure optimization method, device, equipment and medium for deep sea mining, relates to the technical field of deep sea engineering and composite material mechanics, and constructs a material parameter matrix and a geometric parameter matrix of a flexible riser of a non-metal non-bonding structure for deep sea mining at a preset depth; generates a flexibility matrix of each layer of the flexible riser; constructs a contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculates a radial Poisson's ratio; establishes each partial balance equation for representing each layer in the flexible riser, and calculates the overall radial stiffness of the flexible riser; and judges whether the overall radial stiffness is not less than a preset threshold value, if not, the non-bonding structure of the flexible riser is optimized by using an optimization parameter, the problems of accurate modeling of radial mechanical behavior, radial deformation prediction and structure optimization of a multilayer composite riser in an ultrahigh pressure environment are solved, the non-bonding characteristics are fully considered, and efficient and accurate calculation of the radial stiffness is realized.

Description

Technical Field

[0001] This invention relates to the fields of deep-sea engineering and composite material mechanics, and in particular to a method, apparatus, equipment and medium for optimizing the non-bonded structure of flexible risers in deep-sea mining. Background Technology

[0002] Flexible risers for 6000m deep-sea mining must withstand extreme internal transport pressures. Traditional metal risers are limited in application due to their excessive weight and susceptibility to corrosion. Traditional bonding layer stiffness calculation models cannot accommodate the interlaminar slip characteristics of unbonded interfaces, resulting in large deviations in radial stiffness calculations and failing to meet the pressure resistance design requirements of deep-sea risers. Furthermore, there is a lack of theoretical modeling for interlaminar void compression and nonlinear contact behavior under ultra-high transport pressures. Existing models do not consider the coupling relationship between pressure and contact stiffness, making it difficult to accurately predict the structural response under high pressure. In addition, the custom-made deep-sea radial experimental fixtures for unbonded structures are costly, and finite element simulations require the establishment of complex contact pairs, resulting in long calculation times for single-condition operations, which cannot support the rapid iterative design requirements of engineering projects.

[0003] As can be seen from the above, how to solve the problems of accurate modeling of radial mechanical behavior, prediction of radial deformation and structural optimization of multi-layer composite risers under ultra-high pressure, fully consider the unbonded characteristics, achieve efficient and accurate calculation of radial stiffness, and provide theoretical tools for the optimization of cross-sectional parameters and directional design of risers in deep-sea mining are problems to be solved in this field. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a method, apparatus, equipment, and medium for optimizing the unbonded structure of flexible risers in deep-sea mining. This method addresses the difficulties in accurately modeling radial mechanical behavior and predicting and optimizing the radial deformation of multi-layer composite risers under ultra-high pressure environments. It fully considers the unbonded characteristics, achieving efficient and accurate calculation of radial stiffness, and providing a theoretical tool for optimizing the cross-sectional parameters and directional design of risers in deep-sea mining. The specific solution is as follows:

[0005] In a first aspect, this application discloses an optimization method for the non-bonded structure of flexible risers in deep-sea mining, including:

[0006] Construct material parameter matrices and geometric parameter matrices for flexible risers of non-metallic, non-bonded structures for deep-sea mining at a preset depth;

[0007] The flexibility matrix of each layer of the flexible riser is generated based on the material parameter matrix and the geometric parameter matrix.

[0008] A contact balance equation is constructed to characterize the contact balance between two adjacent layers in a flexible riser, and the radial Poisson's ratio is calculated using the contact balance equation and the compliance matrix.

[0009] A plurality of sub-layer balance equations corresponding to layers of the flexible riser are established, and a plurality of radial rigidities of the layers are calculated by using the plurality of sub-layer balance equations, the geometric parameter matrix, the flexibility matrix and the radial Poisson's ratio, and the total radial rigidity of the flexible riser is obtained by superimposing the plurality of radial rigidities of the layers.

[0010] It is judged whether the total radial rigidity is not less than a preset threshold value, and if the total radial rigidity is not less than the preset threshold value, an optimization parameter is generated, and the flexible riser is optimized in a non-bonding structure by using the optimization parameter.

[0011] Optionally, the material parameter matrix and the geometric parameter matrix of the flexible riser of the non-metallic non-bonding structure for deep-sea mining at a preset depth are constructed, and the material parameter matrix and the geometric parameter matrix comprise:

[0012] A pressure-sensitive coefficient matrix is introduced for the flexible riser of the non-metallic non-bonding structure for deep-sea mining at a preset depth, and each layer of material parameters in the pressure-sensitive coefficient matrix is corrected under high pressure to construct a material parameter matrix.

[0013] An independent layer radius matrix is established, and the compression of the interlayer gap is combined with the independent layer radius matrix to construct a geometric parameter matrix.

[0014] The material parameter matrix is:

[0015] ;

[0016] ;

[0017] ;

[0018] wherein, is the elastic modulus, is the corrected elastic modulus, is the elastic modulus correction coefficient, is the Poisson's ratio, is the corrected Poisson's ratio, is the Poisson's ratio correction coefficient, is the elastic modulus, is the corrected elastic modulus, is the elastic modulus correction coefficient, and p is the internal pressure.

[0019] The geometric parameter matrix is:

[0020] ;

[0021] wherein, is the radius of the i-th layer under normal pressure, is the layer thickness under high pressure, is the initial gap thickness, and , is the void compressibility, and .

[0022] Optionally, the compliance matrix of each layer includes a first compliance matrix and a second compliance matrix.

[0023] The first compliance matrix is:

[0024] ;

[0025] ;

[0026] wherein, , , is the elastic modulus of the material in three directions, , , is the Poisson's ratio of the material in three directions, , , is the shear modulus of the material in three directions, , , , , , , , , , , , is the first compliance matrix coefficient, which is related to the modified material parameters;

[0027] The second compliance matrix is:

[0028] ;

[0029] ;

[0030] wherein a is the cosine value of the fiber laying angle, b is the sine value of the fiber laying angle, T is the coordinate conversion matrix in the off-axis coordinate system, is the transpose matrix of T, is the first compliance matrix.

[0031] Optionally, the method further includes constructing a contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating the radial Poisson's ratio by using the contact balance equation and the compliance matrix, including:

[0032] constructing a contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating the radial Poisson's ratio by using the contact balance equation, the first compliance matrix, and the second compliance matrix.

[0033] The contact balance equation is:

[0034] ;

[0035] wherein, Pi is the outer surface pressure of the i-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, Ki is the interlayer contact stiffness, and , hi is the initial void thickness, and ;

[0036] The radial Poisson's ratio calculation formula is:

[0037] ;

[0038] wherein, Pi is the outer surface pressure of the i-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, .

[0039] Optionally, the establishing of the respective layer balance equation for characterizing each layer in the flexible riser comprises:

[0040] establishing the respective layer balance equation for characterizing each layer in the flexible riser based on the second compliance matrix in the compliance matrix, the radial Poisson's ratio, and the geometric parameter matrix;

[0041] The layer balance equation is:

[0042] ;

[0043] wherein, Pi is the outer surface pressure of the i-th layer, , aij is an element of the second compliance matrix, Pi is the outer surface pressure of the i-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, Pi is the outer surface pressure of the i-th layer.

[0044] Optionally, the calculation formula of the radial stiffness of each layer is:

[0045] ;

[0046] wherein, Pi is the outer surface pressure of the i-th layer, Pi is the outer surface pressure of the i-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, Pi+1 is the outer surface pressure of the i+1-th layer, , is an element of the second flexibility matrix, is a radial Poisson's ratio;

[0047] The calculation formula for calculating the overall radial stiffness is:

[0048] ;

[0049] wherein, is a contact stiffness, and , is an interlayer contact stiffness, r is an average radius, and L is a pipeline length.

[0050] Optionally, the method further comprises: judging whether the overall radial stiffness is not less than a preset threshold value; and generating an optimization parameter if the overall radial stiffness is not less than the preset threshold value, the optimization parameter comprising:

[0051] generating the preset threshold value based on a hydrostatic pressure safety factor, a dynamic load coefficient, a material aging coefficient, and a non-bonding safety factor;

[0052] judging whether the overall radial stiffness is not less than the preset threshold value;

[0053] if the overall radial stiffness is not less than the preset threshold value, determining an interlayer gap thickness and an interlayer contact stiffness, and generating the optimization parameter by using the interlayer gap thickness and the interlayer contact stiffness;

[0054] if the overall radial stiffness is less than the preset threshold value, jumping to a flow of constructing a material parameter matrix and a geometric parameter matrix until the overall radial stiffness is not less than the preset threshold value.

[0055] In a second aspect, the present application discloses a flexible riser non-bonding structure optimization device for deep-sea mining, comprising:

[0056] a parameter matrix construction module, configured to construct a material parameter matrix and a geometric parameter matrix for a flexible riser of a non-metal non-bonding structure for deep-sea mining at a preset depth;

[0057] a flexibility matrix construction module, configured to generate a flexibility matrix of each layer of the flexible riser based on the material parameter matrix and the geometric parameter matrix;

[0058] a radial Poisson's ratio calculation module, configured to construct a contact balance equation for characterizing contact between two adjacent layers in the flexible riser, and calculate a radial Poisson's ratio by using the contact balance equation and the flexibility matrix;

[0059] The total radial stiffness calculation module is configured to establish a plurality of layer balance equations for characterizing the corresponding layers in the flexible riser, calculate radial stiffness of each layer by using the plurality of layer balance equations, the geometric parameter matrix, the compliance matrix, and the radial Poisson's ratio, and obtain the total radial stiffness of the flexible riser by superimposing the radial stiffness of each layer.

[0060] The unbonded structure optimization module is configured to determine whether the total radial stiffness is not less than a preset threshold, and if the total radial stiffness is not less than the preset threshold, generate an optimization parameter and perform unbonded structure optimization on the flexible riser by using the optimization parameter.

[0061] In a third aspect, the present application discloses an electronic device, comprising:

[0062] A memory is configured to save a computer program.

[0063] A processor is configured to execute the computer program to implement the deep-sea mining flexible riser unbonded structure optimization method.

[0064] In a fourth aspect, the present application discloses a computer storage medium configured to save a computer program, wherein the computer program is executed by a processor to implement the steps of the deep-sea mining flexible riser unbonded structure optimization method disclosed above.

[0065] It can be seen that the application provides a flexible riser unbonded structure optimization method for deep sea mining, which comprises the following steps: constructing a material parameter matrix and a geometric parameter matrix for a flexible riser of an unbonded structure for deep sea mining at a preset depth; generating a flexibility matrix of each layer of the flexible riser based on the material parameter matrix and the geometric parameter matrix; constructing a contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating a radial Poisson's ratio by using the contact balance equation and the flexibility matrix; establishing a corresponding layered balance equation for representing each layer in the flexible riser, and calculating a radial stiffness of each layer by using each layered balance equation, the geometric parameter matrix, the flexibility matrix and the radial Poisson's ratio; superimposing the radial stiffness of each layer to obtain a total radial stiffness of the flexible riser; and judging whether the total radial stiffness is not less than a preset threshold value, and if the total radial stiffness is not less than the preset threshold value, generating an optimization parameter and optimizing the unbonded structure of the flexible riser by using the optimization parameter. The application aims to solve the problem of accurate and efficient calculation of the radial stiffness of a 6000m deep sea mining flexible riser under an unbonded structure between layers, considers the gap that may occur between layers and the compression of the gap between layers under deep sea pressure, constructs a material parameter matrix and a geometric parameter matrix, generates a flexibility matrix of each layer of the flexible riser, constructs a contact balance equation for representing the contact between two adjacent layers in the flexible riser to calculate a radial Poisson's ratio, establishes a corresponding layered balance equation for representing each layer in the flexible riser to calculate a radial stiffness of each layer, and obtains a total radial stiffness of the flexible riser after superimposition. By constructing a contact mechanics model of an unbonded interface, a layered energy calculation framework and a deep sea pressure correction mechanism, efficient and accurate calculation of the radial stiffness is realized. If the total radial stiffness is not less than the preset threshold value, an optimization parameter is generated, and the unbonded structure of the flexible riser is optimized by using the optimization parameter. The application directly correlates complex geometric, material and other design parameters with the structural stiffness for optimization, so as to quickly provide an optimization parameter and optimize the unbonded structure of the flexible riser, shorten the traditional design cycle, improve the design efficiency of the flexible riser, and meet the needs of rapid iteration of engineering. BRIEF DESCRIPTION OF DRAWINGS

[0066] In order to more clearly illustrate the technical solutions in the embodiments of the application or the prior art, the following will briefly introduce the drawings needed in the embodiments or prior art description. Obviously, the drawings in the following description are only embodiments of the application, and for those skilled in the art, other drawings can be obtained without creative labor based on the provided drawings.

[0067] Figure 1 A flexible riser unbonded structure optimization method for deep sea mining disclosed by the application is shown in the flowchart;

[0068] Figure 2A flexible riser non-bonding structure optimization specific flow chart disclosed by the application;

[0069] Figure 3 A flexible riser non-bonding structure optimization system structure diagram for deep sea mining disclosed by the application;

[0070] Figure 4 A schematic diagram of conversion of a principal coordinate system of an orthotropic material and a cylindrical coordinate system disclosed by the application;

[0071] Figure 5 A flexible riser non-bonding structure optimization device structure schematic diagram for deep sea mining disclosed by the application;

[0072] Figure 6 An electronic device structure diagram provided by the application. DETAILED DESCRIPTION

[0073] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application. Obviously, the described embodiments are only part of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the application.

[0074] The 6000m deep sea mining flexible riser needs to withstand extreme internal conveying pressure, and the traditional metal riser is limited in application due to defects such as excessive self-weight and easy corrosion. The traditional bonding layer stiffness calculation model cannot adapt to the interlayer slip characteristics of the non-bonding interface, resulting in large calculation deviation of the radial stiffness, which cannot meet the design requirements of the deep sea riser pressure resistance; and the interlayer gap compression and contact nonlinear behavior under the super large conveying pressure lacks theoretical modeling, the existing model does not consider the coupling relationship of pressure-contact stiffness, and it is difficult to accurately predict the structural response under high pressure; in addition, the deep sea radial experiment tooling of the non-bonding structure is high in customization cost, and the finite element simulation needs to establish complex contact pairs, and the single working condition calculation is time-consuming, which cannot support the engineering rapid iterative design requirements. As can be seen from the above, how to solve the problems of accurate modeling of radial mechanical behavior, radial deformation prediction and structure optimization of multi-layer composite riser under super high pressure environment, fully considering the non-bonding characteristics, realizing efficient and accurate calculation of radial stiffness, and providing a theoretical tool for sectional parameter optimization and directional design of deep sea mining riser are problems to be solved in the field.

[0075] Referring to Figure 1 The embodiments of the application disclose a flexible riser non-bonding structure optimization method for deep sea mining, which can specifically include:

[0076] Step S11: constructing a material parameter matrix and a geometric parameter matrix for the flexible riser of the non-metal non-bonding structure for deep sea mining at a preset depth.

[0077] In the embodiment, a pressure sensitive coefficient matrix is introduced for the flexible riser of non-metallic non-bonded structure for deep sea mining at a preset depth, and the material parameter in each layer of the pressure sensitive coefficient matrix is high-pressure corrected to construct a material parameter matrix; an independent layer radius matrix is established, and the interlayer gap compression is combined with the independent layer radius matrix to construct a geometric parameter matrix.

[0078] The material parameter matrix is:

[0079]

[0080]

[0081]

[0082] wherein, Ei is the elastic modulus, Ei' is the corrected elastic modulus, Ei' is the corrected elastic modulus, μi is the Poisson's ratio, μi' is the corrected Poisson's ratio, μi' is the corrected Poisson's ratio, Ei is the elastic modulus, Ei' is the corrected elastic modulus, Ei' is the corrected elastic modulus, and p is the internal pressure;

[0083] The geometric parameter matrix is:

[0084]

[0085] wherein, ri is the radius of the i-th layer under normal pressure, hi is the layer thickness under high pressure, hi0 is the initial gap thickness, and hi0 is the initial gap thickness, and hi0 is the initial gap thickness, and .

[0086] The present application defines the material parameter matrix (elastic modulus and Poisson's ratio in different directions) of the non-metallic non-bonded flexible riser preset according to the high-pressure environment of deep sea, and defines the geometric parameter matrix (inner diameter and outer diameter of each layer) of each layer considering the possible gap between layers.

[0087] wherein, the material parameter matrix construction process is: introducing a pressure sensitive coefficient matrix , and high-pressure correcting the material parameters of each layer: wherein, p is the internal pressure, ​​​​The pressure sensitivity coefficient corresponding to different material properties. For the carbon fiber reinforced layer, the parameters such as elastic modulus E and Poisson's ratio v will change. Take the elastic modulus E as an example , and correct it to under high pressure . The typical value of the carbon fiber reinforced layer is . The pressure will change the ratio of the lateral strain to the axial strain of the material under stress, and the Poisson's ratio is corrected to , and the Poisson's ratio correction coefficient of the inner liner HDPE (High Density Polyethylene) is . More accurately reflect the influence of deep sea pressure on the mechanical properties of materials.

[0088] The geometric parameter matrix construction process is: considering that the interlayer gap will be compressed under deep sea pressure, the independent layer radius matrix is established . When calculating the radius of each layer, the interlayer gap compression is considered: . This formula comprehensively considers the influence of layer thickness and gap compression, accurately determines the actual radius of each layer under deep sea pressure, and provides the geometric parameter basis for accurately calculating the stiffness of each layer. In addition, the geometric parameter matrix also includes the laying angle of each layer of composite material .

[0089] Step S12: generating the flexibility matrix of each layer of the flexible riser based on the material parameter matrix and the geometric parameter matrix.

[0090] In this embodiment, the flexibility matrix of each layer includes a first flexibility matrix and a second flexibility matrix;

[0091] The first flexibility matrix is:

[0092] ;

[0093] ;

[0094] wherein, , , are the elastic moduli of the material in three directions, , , are the Poisson's ratios of the material in three directions, , , are the shear moduli of the material in three directions, , , , , , , , , 、 、 、 is the first compliance matrix coefficient, which is related to the corrected material parameters;

[0095] The second compliance matrix is:

[0096] ;

[0097] ;

[0098] wherein a is the cosine value of the fiber laying angle, b is the sine value of the fiber laying angle, T is the coordinate conversion matrix in the off-axis coordinate system, is the transpose matrix of T, is the first compliance matrix.

[0099] The present application defines the compliance matrix of each layer according to the material parameter matrix and the geometric parameter matrix, including the first compliance matrix and the second compliance matrix. The first compliance matrix is used to calculate the stress and strain in the on-axis coordinate system, while the second compliance matrix is used to calculate the stress and strain in the off-axis coordinate system. The first compliance matrix can be converted into the second compliance matrix through the coordinate conversion matrix.

[0100] The construction process of the first compliance matrix is:

[0101] According to the stress-strain relationship, it can be known that: The first compliance matrix is independently constructed for each layer Without considering the interlayer coupling, since there is no shear coupling effect between the non-bonded layers, the corresponding shear terms in the matrix are zero. This makes each layer can be considered independently when force analysis, which conforms to the characteristics of non-bonded structure, so only the radial-circumferential related terms are retained in the off-diagonal elements:

[0102] ;

[0103] The first compliance matrix coefficient is:

[0104] ;

[0105] Since the elastic properties of isotropic materials at different points in different directions are the same, only two parameters, elastic modulus and Poisson's ratio, are needed to express their material properties:

[0106] ;

[0107] The first compliance matrix coefficient of isotropic materials can be further simplified as:

[0108] ;

[0109] The construction process of the second flexibility matrix is as follows:

[0110] The fiber laying angle is known as , , The coordinate conversion matrix T in the off-axis coordinate system is:

[0111] ;

[0112] The coordinate conversion matrix T is used to convert the stress-strain relationship of the material in the off-axis coordinate system to the on-axis coordinate system for subsequent calculation, so the second flexibility matrix can be expressed as: .

[0113] ;

[0114] The determination of the elements is based on the geometric relationship of coordinate transformation and the operation of trigonometric functions. For example, the elements such as , in the matrix are obtained according to the square of the cosine and sine values of the angle, and the cross terms such as 2ab are obtained by product operation of trigonometric functions.

[0115] Step S13: Construct a contact balance equation for representing the contact balance between adjacent two layers in the flexible riser, and calculate the radial Poisson's ratio by using the contact balance equation and the flexibility matrix.

[0116] In this embodiment, a contact balance equation for representing the contact balance between adjacent two layers in the flexible riser is constructed, and the radial Poisson's ratio is calculated by using the contact balance equation, the first flexibility matrix, and the second flexibility matrix;

[0117] The contact balance equation is:

[0118] ;

[0119] wherein, is the outer surface pressure of the i-th layer, is the outer surface pressure of the i+1-th layer, is the interlayer contact stiffness, and , is the initial gap thickness, and ;

[0120] The radial Poisson's ratio is independently derived for each layer, and the radial Poisson's ratio calculation formula is:

[0121] ;

[0122] wherein, is the radial Poisson's ratio, Typical values for the composite layers are 0.3 for E and 0.3 for v, and .

[0123] The non-bonded flexible riser is separable between each layer, the most prominent feature is the radial interlayer contact pressure and the radial Poisson's ratio of two key parameters. According to the load constraints of the flexible riser, the interlayer contact pressure transfer calculation is carried out, and the radial Poisson's ratio of each layer is independently derived.

[0124] The present application assumes that the non-bonded interface only transmits radial pressure and does not transmit shear force, and establishes a contact balance equation. Under the action of deep sea pressure, each layer will deform radially, and since the layers are non-bonded, there will be a displacement difference between the layers. When there is relative displacement between the layers, the contact stiffness will work, and when the internal pressure p is known, the pressure distribution between the layers is adjusted according to this balance equation. In the non-bonded structure, due to the existence of interlayer slip, the transverse deformation will be larger than that of the traditional bonded structure, so a slip influence coefficient is introduced to correct the radial Poisson's ratio.

[0125] Step S14: Establishing each layer balance equation for characterizing each layer in the flexible riser, calculating the radial stiffness of each layer by using each layer balance equation, the geometric parameter matrix, the flexibility matrix, and the radial Poisson's ratio, superimposing the radial stiffness of each layer to obtain the total radial stiffness of the flexible riser.

[0126] In this embodiment, each layer balance equation for characterizing each layer in the flexible riser is established based on the second flexibility matrix in the flexibility matrix, the radial Poisson's ratio, and the geometric parameter matrix;

[0127] The layer balance equation is:

[0128] ;

[0129] Wherein, the total strain energy is , the element of the second flexibility matrix is the radial Poisson's ratio is the radius of the i-th layer under the deep sea pressure is the radius of the i+1-th layer under the deep sea pressure is the radial load is.

[0130] The radial stiffness of each layer is calculated by using each layer balance equation, the geometric parameter matrix, the flexibility matrix, and the radial Poisson's ratio, and the calculation formula of the radial stiffness of each layer is:

[0131] ;

[0132] Wherein, Ri is the radius of the ith layer under the deep-sea pressure, Ri is the radius of the ith layer under the deep-sea pressure, Ri is the radius of the ith layer under the deep-sea pressure, Ri is the radius of the ith layer under the deep-sea pressure, , Ri is the radius of the ith layer under the deep-sea pressure, Ri is the radius of the ith layer under the deep-sea pressure;

[0133] The calculation formula of the total radial stiffness is:

[0134] ;

[0135] wherein, is the contact stiffness, and , is the interlayer contact stiffness, r is the average radius, and L is the pipeline length.

[0136] The present application first establishes the equilibrium equation of each layer to solve the independent radial stiffness of each layer, and further deduces the stiffness expression of the non-bonded flexible riser in combination with the cross-section contact stiffness. The construction process of the layered equilibrium equation is as follows: the strain energy of each layer is calculated independently , without considering the interlayer cooperation, and the equilibrium equation is established according to the work done by the radial load . According to the energy balance principle, the work done by the radial load is equal to the increase of the strain energy of each layer, and the layered equilibrium equation is established through this relationship, wherein is the second flexibility matrix element of each layer, which is obtained through coordinate transformation. For a certain layer, such as a composite material layer, the inner and outer radii of the layer under the deep-sea pressure and are determined first, and then the integral formula is calculated according to the second flexibility matrix element of the layer. The integral here is the integral of the volume of the layer, and the total strain energy of the layer is obtained by calculating the integral of the strain energy density over the volume. Since it is a non-bonded structure, the strain energy of each layer is calculated independently, without considering the influence of interlayer cooperative deformation.

[0137] The calculation formula derivation process of the radial stiffness is as follows: the independent radial stiffness of each layer can be solved according to the independent equilibrium equation of each layer, and the total radial stiffness is deduced based on the equilibrium equation of each layer and the interface contact balance. First, the independent radial stiffness of each layer is calculated respectively, which can be deduced through the relationship between the strain energy and the displacement of each layer. Then, the stiffness of each layer is superimposed by considering the influence of the interlayer contact stiffness . This formula comprehensively considers the mechanical properties of each layer and the influence of interlayer contact, and can accurately calculate the radial stiffness of the non-bonded structure.

[0138] Step S15: judging whether the overall radial stiffness is not less than a preset threshold value, if the overall radial stiffness is not less than the preset threshold value, generating an optimization parameter, and using the optimization parameter to perform non-bonding structure optimization on the flexible riser.

[0139] In this embodiment, the preset threshold value is generated based on a hydrostatic pressure safety factor, a dynamic load coefficient, a material aging coefficient and a non-bonding safety factor; it is judged whether the overall radial stiffness is not less than the preset threshold value; if the overall radial stiffness is not less than the preset threshold value, the interlayer gap thickness and the interlayer contact stiffness are determined, and the optimization parameter is generated using the interlayer gap thickness and the interlayer contact stiffness; if the overall radial stiffness is less than the preset threshold value, the process of constructing the material parameter matrix and the geometric parameter matrix is jumped to, until the overall radial stiffness is not less than the preset threshold value.

[0140] The radial stiffness of the non-bonding flexible riser can be solved through S11-S14, it is judged whether the overall radial stiffness meets the design through the preset threshold value, if not, it is returned to S11 for continuous optimization, until the condition is met to end the cycle. In addition, the non-bonding parameters (gap thickness and contact stiffness) can also be selected and adjusted in the cycle.

[0141] Specifically, the preset threshold value is generated based on a hydrostatic pressure safety factor , a dynamic load coefficient , a material aging coefficient , and a non-bonding safety factor , considering the uncertainty of interface slip, the preset threshold value is taken as a safety criterion, and the following conditions are met:

[0142] ;

[0143] Wherein, design requirement stiffness. In deep sea environment, the riser not only bears the hydrostatic pressure, but also considers the influence of dynamic load (such as vibration caused by ocean current) and material aging on the stiffness. The non-bonding safety factor is specially introduced for the uncertainty of interlayer slip in non-bonding structure. Through the three safety criteria, it can be ensured that the calculated radial stiffness meets the safety requirements of deep sea mining riser under various working conditions. If the design requirement stiffness is not met, the parameters are adjusted (such as increasing the thickness of carbon fiber layer, optimizing the winding angle), and it is returned to S11 for iteration, until the optimization parameter meeting the 6000m deep sea environment is generated, which is finally used to realize the non-bonding structure optimization and design of the flexible riser.

[0144] In this embodiment, the optimization parameter includes the optimized interlayer gap thickness and the contact stiffness , and the typical scheme is as follows:

[0145] between the inner liner and the reinforcing layer ;

[0146] between the reinforcing layer and the outer cover .

[0147] interlayer gap thickness and contact stiffness The radial stiffness of the unbonded structure is significantly affected by these two parameters. By adjusting these two parameters, the contact state and mechanical properties between the layers can be changed. For example, reducing the gap thickness between the inner liner and the reinforcing layer can make the interlayer contact more closely, thereby increasing the radial stiffness; while increasing the contact stiffness can enhance the resistance to displacement difference between the layers, which will also affect the radial stiffness. By continuously trying different parameter values, the optimal parameters that meet the safety criteria are found and generated, realizing the optimization design of the unbonded structure.

[0148] The specific process of the flexible riser unbonded structure optimization in the present application is shown in Figure 2 , taking a 6000m deep-sea mining non-metallic unbonded flexible riser as an example: first, the initial design parameters are designed, the structure is: HDPE inner liner (2mm) + carbon fiber reinforcing layer (8 layers ± 55°) + PA outer cover (3mm); the initial gap between the layers is: inner liner-reinforcing layer , reinforcing-outer cover ; the target radial stiffness is: , then the unbonded stiffness calculation is carried out, the layered stiffness calculation: inner liner: =320 kN / m; reinforcing layer: =1200 kN / m; outer cover: =280 kN / m; contact stiffness calculation:

[0149] ;

[0150] System radial stiffness synthesis:

[0151] ;

[0152] Finally, the unbonded structure parameter optimization is carried out, the adjustment scheme is: increasing the reinforcing layer to 12 layers, the gap compression coefficient from 0.05 to 0.08; the inner liner-reinforcing layer gap is reduced to 0.1mm, and the contact stiffness is increased to 80MPa / mm. After optimization: ;

[0153] It should be noted that the safety factor needs to be reverified, and this is only an example of the method.

[0154] The present application also proposes an unbonded structure optimization system for a deep-sea mining flexible riser, and the structure is as followsFigure 3 As shown, it includes a preprocessing module for unbonded structural parameters, a layered stiffness calculation module, a contact stiffness synthesis module, an energy method calculation module, and an unbonded optimization module.

[0155] The non-bonding parameter preprocessing module includes:

[0156] Interlayer contact database: This database stores over 200 sets of unbonded interface parameters, covering different combinations of initial void thickness and contact stiffness. , Each set of parameters in the database contains detailed records of the corresponding simulation conditions, enabling the quick and accurate retrieval of suitable initial parameter values ​​based on specific deep-sea mining riser design requirements in practical applications, thus providing a reliable foundation for subsequent calculations.

[0157] The pressure-void model effectively handles the complex nonlinear relationship between pressure and void compressibility through its nonlinear mapping capability. During model training, a large amount of void thickness data under different pressure values ​​is collected as training samples. This data originates from high-precision numerical simulation results. By inputting this data into the model, the nonlinear model learns the inherent mapping law between pressure and void compressibility. This model can accurately predict the corresponding void compressibility based on the input deep-sea pressure value. In actual calculations, inputting the deep-sea pressure value into the model quickly yields the predicted void compressibility value, providing key parameters for subsequent calculations of the unbonded layer radius matrix and overall radial stiffness.

[0158] The core module for unbonded stiffness calculation includes a layered stiffness calculation module, a contact stiffness synthesis module, and an energy method calculation module. It also includes:

[0159] Layered Matrix Operation Engine: This engine boasts powerful parallel computing capabilities and is specifically optimized for stiffness matrix operations on unbonded layers. When processing multi-layered unbonded structures, traditional serial computation methods are inefficient because each layer has its own independent stiffness matrix. The layered matrix operation engine utilizes the multi-core processor architecture of modern computers, distributing the stiffness matrix calculation tasks of different layers across multiple processor cores for simultaneous processing. For example, for a non-metallic, unbonded flexible riser structure containing more than 100 layers, it can simultaneously calculate the first compliance matrix, second compliance matrix, etc., of each layer. By employing efficient matrix operation algorithms and parallel computing strategies, computation time is significantly reduced, enabling large-scale matrix operation tasks to be completed in a short time, providing strong support for the rapid calculation of radial stiffness in unbonded structures.

[0160] Contact iterative solver: Newton-Raphson method is adopted as the core algorithm for solving the interlayer pressure-displacement relationship. In the non-bonded structure, there is a complex nonlinear relationship between the interlayer pressure and displacement, which needs to be solved gradually through iteration. Newton-Raphson method has the advantages of fast convergence speed and high precision. In each iteration process, the gradient and Jacobian matrix of the function are calculated according to the current pressure and displacement values, and then the pressure and displacement values are updated by solving the linear equations to approach the true solution. Through repeated iterative calculation, the contact iterative solver can accurately solve the interlayer pressure-displacement relationship, thereby providing key intermediate results for accurate calculation of the radial stiffness of the non-bonded structure. In practical application, it can quickly and stably converge to the solution that meets the accuracy requirements, ensuring the accuracy and reliability of the entire calculation system.

[0161] Non-bonding optimization module, comprising:

[0162] Void-stiffness sensitivity analysis: automatically generated Curve, for example, when the interlayer void thickness of the inner liner and the reinforcing layer increases from 0.1 mm to 0.15 mm, the system automatically calculates 10 groups of equidistant sample points, solves the stiffness values of each working condition by layering energy method, and forms a continuous curve. Based on the inflection point of the curve, the critical void thickness is identified, and when it is greater than the threshold value, the radial stiffness decay rate significantly accelerates, which is used as the design threshold of the void thickness.

[0163] The non-metallic non-bonding flexible riser mainly includes an inner liner, a reinforcing layer, and an outer cover. The inner liner and the outer cover are made of non-metallic polymer materials such as PE (Polyethylene) and high-density polyethylene. The reinforcing layer is made of composite materials (aramid fiber, glass fiber, and carbon fiber, etc.). All functional layers are non-bonded, with typical orthotropic anisotropy characteristics. The conversion of the orthotropic anisotropy material principal coordinate system and the cylindrical coordinate system is shown in Figure 4 .

[0164] The innovation of the present application lies in breaking through the traditional bonding assumption, establishing a "contact stiffness-void compression" coupling model, and quantifying the interlayer slip effect through the radial stiffness formula, where the contact stiffness and the void compression coefficient constitute the core calculation logic of the patent protection. The layered cumulative method based on energy conservation is proposed, and the strain energy of each layer is calculated independently and coupled through the contact balance equation to solve the interlayer non-synchronous deformation problem of the non-bonded structure, which is specifically reflected in the independent integration and superposition mechanism of the contribution of each layer in the radial stiffness calculation formula. A three-dimensional correction system including elastic modulus correction coefficient, radial shrinkage coefficient, and Poisson's ratio correction coefficient is constructed to realize dynamic adjustment of material parameters under high pressure, such as the linear correction model of .

[0165] The differences between the present application and the prior art include the following points:

[0166] 1. Dimensional difference of mechanical model: The existing axial stiffness calculation patents mainly focus on the axial mechanical behavior, and construct the calculation system with axial load as the core. However, the present application focuses on the radial stiffness calculation, and considers the mechanical response of the non-bonding structure under the action of radial pressure under deep-sea high pressure. The two are different in mechanical dimension, and are different emphases in the study of the mechanical properties of non-bonding flexible risers. The present application solves the problem of accurate modeling of radial mechanical behavior by establishing a non-bonding multi-layer parameter matrix and deriving the radial coupling parameters of the non-bonding interface, fills the technical gap in the field of radial stiffness calculation, and perfects the research system of the mechanical properties of non-bonding flexible risers.

[0167] 2. Difference in interlayer relationship processing: The existing axial stiffness calculation related technologies assume that the interlayer is in a bonded state, and the deformation of each layer is consistent, and the axial load is borne together. However, the present application is aimed at a non-bonding flexible riser, and the interlayer realizes non-bonding contact through a lubricating layer or a gap, allowing the interlayer to separate and slip radially. The present application fully considers this non-bonding characteristic in the calculation process, accurately simulates the mechanical relationship between the layers by constructing an interlayer contact pressure transmission model and calculating the non-bonding radial Poisson's ratio, and is essentially different from the calculation model under the existing bonding assumption, and is more in line with the actual working state of the non-bonding flexible riser.

[0168] 3. Difference in environmental factor consideration: The existing axial stiffness calculation patents have relatively limited specific consideration of deep-sea environment, mainly focusing on mechanical calculation under normal working conditions. However, the present application is aimed at a 6000m deep-sea mining environment, and deeply studies the influence of 60MPa hydrostatic pressure on material parameters, interlayer gap and contact state. For example, by introducing a pressure-sensitive coefficient matrix to modify the material parameters in deep sea, considering the influence of pressure on the compression of the interlayer gap to define the cross-section geometry of the non-bonding layer, and fully and deeply integrating the deep-sea high-pressure environmental factors into every aspect of the radial stiffness calculation, the calculation results are more in line with the actual deep-sea working conditions.

[0169] In this embodiment, a material parameter matrix and a geometric parameter matrix are constructed for a flexible riser with a non-metallic, non-bonded structure for deep-sea mining at a preset depth. Based on the material parameter matrix and the geometric parameter matrix, a compliance matrix for each layer of the flexible riser is generated. A contact balance equation is constructed to characterize the contact balance equation between adjacent layers in the flexible riser, and the radial Poisson's ratio is calculated using the contact balance equation and the compliance matrix. Layer-specific balance equations are established to characterize each layer in the flexible riser, and the radial stiffness of each layer is calculated using the layer-specific balance equations, the geometric parameter matrix, the compliance matrix, and the radial Poisson's ratio. The radial stiffness of each layer is superimposed to obtain the overall radial stiffness of the flexible riser. It is determined whether the overall radial stiffness is not less than a preset threshold. If the overall radial stiffness is not less than the preset threshold, optimization parameters are generated, and the flexible riser is optimized for non-bonded structure using the optimization parameters. This application aims to solve the problem of accurately and efficiently calculating the radial stiffness of flexible risers in 6000m deep-sea mining under unbonded interlayer structures. Considering the potential gaps between layers and the compression of these gaps under deep-sea pressure, material parameter matrices and geometric parameter matrices are constructed to generate compliance matrices for each layer of the flexible riser. Contact equilibrium equations are constructed to characterize the contact equilibrium between adjacent layers in the flexible riser to calculate the radial Poisson's ratio. Layer-specific equilibrium equations are established to characterize the radial stiffness of each layer in the flexible riser, and these equations are superimposed to obtain the overall radial stiffness of the flexible riser. By constructing a contact mechanics model of the unbonded interface, a layered energy calculation framework, and a deep-sea pressure correction mechanism, efficient and accurate calculation of radial stiffness is achieved. If the overall radial stiffness is not less than a preset threshold, optimization parameters are generated. These optimization parameters are then used to optimize the unbonded structure of the flexible riser. This application directly correlates and optimizes complex geometric and material design parameters with structural stiffness, enabling rapid provision of optimization parameters and optimization of the unbonded structure of the flexible riser. This shortens the traditional design cycle, improves the design efficiency of flexible risers, and meets the needs of rapid engineering iteration.

[0170] See Figure 5 As shown in the figure, this invention discloses a flexible riser non-bonded structure optimization device for deep-sea mining, which may specifically include:

[0171] The parameter matrix construction module 11 is used to construct the material parameter matrix and geometric parameter matrix for the flexible riser of non-metallic, non-bonded structure for deep-sea mining at a preset depth.

[0172] The flexibility matrix construction module 12 is used to generate the flexibility matrix of each layer of the flexible riser based on the material parameter matrix and the geometric parameter matrix.

[0173] Radial Poisson's ratio calculation module 13 is used to construct a contact balance equation to characterize the contact balance between two adjacent layers in a flexible riser, and to calculate the radial Poisson's ratio using the contact balance equation and the compliance matrix.

[0174] a total radial stiffness calculation module 14 configured to establish a plurality of layering balance equations for characterizing each layer in the flexible riser, calculate radial stiffness of each layer by using the plurality of layering balance equations, the geometric parameter matrix, the compliance matrix, and the radial Poisson's ratio, and superimpose the radial stiffness of each layer to obtain a total radial stiffness of the flexible riser;

[0175] a non-bonded structure optimization module 15 configured to determine whether the total radial stiffness is not less than a preset threshold, generate an optimization parameter if the total radial stiffness is not less than the preset threshold, and perform non-bonded structure optimization on the flexible riser by using the optimization parameter.

[0176] In this embodiment, the material parameter matrix and the geometric parameter matrix of the flexible riser of the non-metallic non-bonded structure for deep-sea mining at a preset depth are constructed, the compliance matrix of each layer of the flexible riser is generated based on the material parameter matrix and the geometric parameter matrix, a contact balance equation for characterizing the contact between two adjacent layers in the flexible riser is constructed, the radial Poisson's ratio is calculated by using the contact balance equation and the compliance matrix, a plurality of layering balance equations for characterizing each layer in the flexible riser are established, the radial stiffness of each layer is calculated by using the plurality of layering balance equations, the geometric parameter matrix, the compliance matrix, and the radial Poisson's ratio, and the radial stiffness of each layer is superimposed to obtain the total radial stiffness of the flexible riser. It is determined whether the total radial stiffness is not less than a preset threshold. If the total radial stiffness is not less than the preset threshold, an optimization parameter is generated, and non-bonded structure optimization is performed on the flexible riser by using the optimization parameter. This application aims to solve the problem of accurate and efficient calculation of the radial stiffness of the flexible riser of the non-bonded structure for deep-sea mining at a depth of 6000 m. The gap between layers and the compression of the gap under deep-sea pressure are considered. The material parameter matrix and the geometric parameter matrix are constructed, the compliance matrix of each layer of the flexible riser is generated, the contact balance equation for characterizing the contact between two adjacent layers in the flexible riser is constructed to calculate the radial Poisson's ratio, the plurality of layering balance equations for characterizing each layer in the flexible riser are established to calculate the radial stiffness of each layer, and the total radial stiffness of the flexible riser is obtained after superimposition. By constructing the contact mechanics model of the non-bonded interface, the layering energy calculation framework, and the deep-sea pressure correction mechanism, efficient and accurate calculation of the radial stiffness is realized. If the total radial stiffness is not less than the preset threshold, an optimization parameter is generated, and non-bonded structure optimization is performed on the flexible riser by using the optimization parameter. This application directly correlates the complex geometric, material, and other design parameters with the structural stiffness for optimization, so as to quickly provide the optimization parameter and perform non-bonded structure optimization on the flexible riser, shorten the traditional design cycle, improve the design efficiency of the flexible riser, and meet the needs of rapid iteration in engineering.

[0177] In some embodiments, the parameter matrix construction module 11 can specifically include:

[0178] a high-pressure correction module for introducing a pressure-sensitive coefficient matrix for the flexible riser of the non-metallic non-bonded structure of the deep-sea mining at a preset depth, and correcting the material parameters in each layer of the pressure-sensitive coefficient matrix to construct a material parameter matrix;

[0179] a geometric parameter matrix construction module for establishing an independent layer radius matrix, and combining the interlayer gap compression with the independent layer radius matrix to construct a geometric parameter matrix;

[0180] The material parameter matrix is:

[0181] ;

[0182] ;

[0183] ;

[0184] wherein, is the elastic modulus, is the corrected elastic modulus, is the elastic modulus correction coefficient, is the Poisson's ratio, is the corrected Poisson's ratio, is the Poisson's ratio correction coefficient, is the elastic modulus, is the corrected elastic modulus, is the elastic modulus correction coefficient, and p is the internal pressure;

[0185] The geometric parameter matrix is:

[0186] ;

[0187] wherein, is the radius of the i-th layer under normal pressure, is the layer thickness under high pressure, is the initial gap thickness, and , is the gap compression coefficient, and .

[0188] In some embodiments, the layer flexibility matrix includes a first flexibility matrix and a second flexibility matrix;

[0189] The first flexibility matrix is:

[0190] ;

[0191] ;

[0192] wherein, , , are the elastic moduli of the material in three directions, , , are the Poisson's ratios of the material in three directions, , , are the shear moduli of the material in three directions, , , , , , , , , , , , is a first compliance matrix coefficient related to the modified material parameters;

[0193] The second compliance matrix is:

[0194] ;

[0195] ;

[0196] wherein a is a cosine value of the fiber laying angle, b is a sine value of the fiber laying angle, T is a coordinate conversion matrix in an off-axis coordinate system, is a transpose matrix of T, is the first compliance matrix.

[0197] In some specific embodiments, the radial Poisson's ratio calculation module 13 can specifically include:

[0198] a contact balance equation construction module, configured to construct a contact balance equation for characterizing contact balance between two adjacent layers in the flexible riser, and calculate the radial Poisson's ratio by using the contact balance equation, the first compliance matrix and the second compliance matrix.

[0199] The contact balance equation is:

[0200] ;

[0201] wherein, is an outer surface pressure of an i-th layer, is an outer surface pressure of an i+1-th layer, is an interlayer contact stiffness, and , is an initial gap thickness, and ;

[0202] The radial Poisson's ratio calculation formula is:

[0203] ;

[0204] wherein, is the radial Poisson's ratio, is a typical value of the composite material layer, and .

[0205] In some specific embodiments, the total radial stiffness calculation module 14 can specifically include:

[0206] A layered balance equation establishment module is configured to establish, based on the second flexibility matrix in the flexibility matrix, the radial Poisson's ratio, and the geometric parameter matrix, a corresponding layered balance equation for each layer in the flexible riser.

[0207] The layered balance equation is:

[0208] ;

[0209] wherein, is the total strain energy, , is an element of the second flexibility matrix, is the radial Poisson's ratio, is the radius of the i-th layer under the deep sea pressure, is the radius of the i+1-th layer under the deep sea pressure, is the radial load.

[0210] In some specific embodiments, the calculation formula of the radial stiffness of each layer is:

[0211] ;

[0212] wherein, is the radial stiffness of each layer, is the radial load, is the radius of the i-th layer under the deep sea pressure, is the radius of the i+1-th layer under the deep sea pressure, , is an element of the second flexibility matrix, is the radial Poisson's ratio;

[0213] The calculation formula for calculating the total radial stiffness is:

[0214] ;

[0215] wherein, is the contact stiffness, and , where r is the average radius, L is the length of the pipe.

[0216] In some embodiments, the non-bonded structure optimization module 15 can specifically include:

[0217] A preset threshold generation module is configured to generate a preset threshold based on a hydrostatic pressure safety factor, a dynamic load coefficient, a material aging coefficient, and a non-bonded safety factor.

[0218] A judgment module is configured to judge whether the overall radial stiffness is not less than the preset threshold.

[0219] An optimization parameter generation module is configured to, if the overall radial stiffness is not less than the preset threshold, determine an interlayer gap thickness and an interlayer contact stiffness, and generate an optimization parameter by using the interlayer gap thickness and the interlayer contact stiffness.

[0220] A flow jump module is configured to, if the overall radial stiffness is less than the preset threshold, jump to a flow of constructing a material parameter matrix and a geometric parameter matrix until the overall radial stiffness is not less than the preset threshold.

[0221] Figure 6 A structural schematic diagram of an electronic device provided by an embodiment of the present application is provided. The electronic device 20 can specifically include at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is configured to store a computer program, and the computer program is loaded and executed by the processor 21 to implement the related steps in the deep-sea mining flexible riser non-bonded structure optimization method performed by the electronic device disclosed in any of the preceding embodiments.

[0222] In the embodiment, the power supply 23 is configured to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol followed by the communication interface 24 can be any communication protocol applicable to the technical solution of the present application, which is not specifically limited here; the input / output interface 25 is configured to obtain external input data or output data to the outside world, and the specific interface type can be selected according to the specific application needs, which is not specifically limited here.

[0223] In addition, the memory 22 as a carrier for resource storage can be a read-only memory, a random access memory, a magnetic disk, or an optical disk, etc., and the resources stored thereon include an operating system 221, a computer program 222, and data 223, etc., and the storage mode can be temporary storage or permanent storage.

[0224] The operating system 221 is used to manage and control each hardware device on the electronic device 20 and the computer program 222, so as to realize the operation and processing of the processor 21 on the data 223 in the memory 22, and can be Windows, Unix, Linux, etc. The computer program 222 can further include a computer program capable of completing other specific work in addition to the computer program capable of completing the flexible riser unbonded structure optimization method for deep-sea mining performed by the electronic device 20 disclosed in any of the foregoing embodiments. The data 223 can include data received by the deep-sea mining flexible riser unbonded structure optimization device from an external device, data collected by the I / O interface 25, and the like.

[0225] The steps of a method or algorithm described in connection with the embodiments disclosed herein can be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module can reside in random access memory (RAM), flash memory, read-only memory (ROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0226] Further, the embodiments of the present application also disclose a computer readable storage medium, wherein the storage medium stores a computer program, and the computer program is loaded and executed by a processor to implement the steps of the flexible riser unbonded structure optimization method for deep-sea mining disclosed in any of the foregoing embodiments.

[0227] Finally, it should be noted that, in this document, the terms such as first and second are merely used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is any such actual relationship or order between these entities or operations. Moreover, the terms "include", "contain" 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 not only includes those elements, but also includes other elements not explicitly listed or inherent to such a process, method, article or device. Without more limitations, the element defined by the statement "including a" does not exclude the presence of another identical element in the process, method, article or device including the element.

[0228] The above describes in detail the flexible riser non-bonding structure optimization method, device, equipment and storage medium for deep sea mining provided by the present application, the principle and implementation mode of the present application are described by applying specific examples, and the above example is only used to help understand the method of the present application and its core idea; meanwhile, for those skilled in the art, according to the idea of the present application, the specific implementation mode and application range will be changed, and the above description should not be understood as the limitation of the present application.

Claims

1. A method of flexible riser unbonded structure optimization for deep sea mining, characterized in that, The method comprises the following steps: constructing a material parameter matrix and a geometric parameter matrix of a flexible riser of a non-metallic non-bonded structure for deep-sea mining at a preset depth; generating a layer flexibility matrix of the flexible riser based on the material parameter matrix and the geometric parameter matrix; constructing a contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating a radial Poisson's ratio by using the contact balance equation and the flexibility matrix; establishing a corresponding layer balance equation for representing each layer in the flexible riser, and calculating a radial stiffness of each layer by using the layer balance equation, the geometric parameter matrix, the flexibility matrix and the radial Poisson's ratio, and superimposing the radial stiffness of each layer to obtain a total radial stiffness of the flexible riser; judging whether the total radial stiffness is not less than a preset threshold value, and if the total radial stiffness is not less than the preset threshold value, generating an optimization parameter, and optimizing the non-bonded structure of the flexible riser by using the optimization parameter.

2. A method of flexible riser unbonded structure optimization for deep-sea mining according to claim 1, characterized in that, The method for constructing the material parameter matrix and the geometric parameter matrix of the flexible riser of the non-metallic non-bonded structure for deep-sea mining at the preset depth comprises the following steps: introducing a pressure sensitivity coefficient matrix to the flexible riser of the non-metallic non-bonded structure for deep-sea mining at the preset depth, and correcting each layer material parameter in the pressure sensitivity coefficient matrix to construct a material parameter matrix; establishing an independent layer radius matrix, and combining the interlayer gap compression with the independent layer radius matrix to construct a geometric parameter matrix; the material parameter matrix is: ； ； ； wherein, E is the elastic modulus, E is the corrected elastic modulus, E is the elastic modulus correction coefficient, v is the Poisson's ratio, v is the corrected Poisson's ratio, v is the Poisson's ratio correction coefficient, E is the elastic modulus, E is the corrected elastic modulus, E is the elastic modulus correction coefficient, p is the internal pressure; the geometric parameter matrix is: ； wherein Ri is the radius of the ith layer at normal pressure, hi is the layer thickness at high pressure, hi0is the initial void thickness, and , Kv is the void compressibility, and .

3. The method of flexible riser unbonded structure optimization for deep-sea mining of claim 1, wherein, the layer flexibility matrix comprises a first flexibility matrix and a second flexibility matrix; the first flexibility matrix is: ； ； wherein , , Ei are the elastic moduli of the material in the three directions, , , Pi are the Poisson's ratios of the material in the three directions, , , Gi are the shear moduli of the material in the three directions, , , , , , , , , , , , are the first matrix coefficients related to the modified material parameters; the second flexibility matrix is: ； ； Wherein, a is the cosine value of the fiber laying angle, b is the sine value of the fiber laying angle, T is the coordinate conversion matrix in the off-axis coordinate system, is the transpose matrix of T, is the first flexibility matrix.

4. A method of flexible riser unbonded structure optimization for deep sea mining according to claim 3, characterized in that, the method for constructing the contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating the radial Poisson's ratio by using the contact balance equation and the flexibility matrix comprises the following steps: constructing the contact balance equation for representing the contact between two adjacent layers in the flexible riser, and calculating the radial Poisson's ratio by using the contact balance equation, the first flexibility matrix and the second flexibility matrix; the contact balance equation is: ； wherein, is the pressure at the outer surface of the i-th layer, is the pressure at the outer surface of the i+1-th layer, is the interlayer contact stiffness, and , is the initial void thickness, and ; the formula for calculating the radial Poisson's ratio is: ； wherein, is the radial Poisson's ratio, is a typical value for composite layers, and .

5. The method of flexible riser unbonded structure optimization for deep-sea mining of claim 3, wherein, the method for establishing the layer balance equation for representing each layer in the flexible riser comprises the following steps: establishing the layer balance equation for representing each layer in the flexible riser based on the second flexibility matrix in the flexibility matrix, the radial Poisson's ratio and the geometric parameter matrix; the layer balance equation is: ； wherein, is the total strain energy, , is an element of the second flexibility matrix, is the radial Poisson's ratio, is the radius of the i-th layer under the deep-sea pressure, is the radius of the i+1-th layer under the deep-sea pressure, is the radial load.

6. The method of flexible riser unbonded structure optimization for deep sea mining of claim 3, wherein, the formula for calculating the radial stiffness of each layer is: ； wherein, is the radial stiffness of each layer, is the radial load, is the radius of the i-th layer under deep-sea pressure, is the radius of the i+1-th layer under deep-sea pressure, , is an element of the second flexibility matrix, is the radial Poisson's ratio; the formula for calculating the total radial stiffness is: ； wherein is the contact stiffness, and , is the interlayer contact stiffness, r is the average radius, and L is the pipe length.

7. A method of flexible riser unbonded structure optimization for deep sea mining according to any one of claims 1 to 6, characterized in that, the method for judging whether the total radial stiffness is not less than the preset threshold value, and generating the optimization parameter if the total radial stiffness is not less than the preset threshold value comprises the following steps: generating the preset threshold value based on a hydrostatic pressure safety factor, a dynamic load coefficient, a material aging coefficient and a non-bonding safety factor; judging whether the total radial stiffness is not less than the preset threshold value; if the total radial stiffness is not less than the preset threshold value, determining an interlayer gap thickness and an interlayer contact stiffness, and generating the optimization parameter by using the interlayer gap thickness and the interlayer contact stiffness; if the total radial stiffness is less than the preset threshold value, jumping to the process of constructing the material parameter matrix and the geometric parameter matrix until the total radial stiffness is not less than the preset threshold value.

8. A flexible riser unbonded structure optimisation apparatus for deep sea mining, characterised by, The method comprises the following steps: a parameter matrix construction module is configured to construct a material parameter matrix and a geometric parameter matrix for a flexible riser of a non-metallic non-bonded structure for deep-sea mining at a preset depth; a flexibility matrix construction module is configured to generate a flexibility matrix of each layer of the flexible riser based on the material parameter matrix and the geometric parameter matrix; a radial Poisson's ratio calculation module is configured to construct a contact balance equation for characterizing the contact between two adjacent layers in the flexible riser, and calculate a radial Poisson's ratio by using the contact balance equation and the flexibility matrix; a total radial stiffness calculation module is configured to establish a corresponding layer balance equation for characterizing each layer of the flexible riser, and calculate a radial stiffness of each layer by using each layer balance equation, the geometric parameter matrix, the flexibility matrix, and the radial Poisson's ratio, and obtain a total radial stiffness of the flexible riser by superimposing the radial stiffness of each layer; a non-bonded structure optimization module is configured to determine whether the total radial stiffness is not less than a preset threshold, and if the total radial stiffness is not less than the preset threshold, generate an optimization parameter and optimize the non-bonded structure of the flexible riser by using the optimization parameter.

9. An electronic device, comprising: The method comprises the following steps: a memory is configured to save a computer program; a processor is configured to execute the computer program to implement the method for optimizing a non-bonded structure of a flexible riser for deep-sea mining according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, a memory is configured to save a computer program; wherein the computer program is executed by a processor to implement the method for optimizing a non-bonded structure of a flexible riser for deep-sea mining according to any one of claims 1 to 7.

Images