CN116882320B - Stability evaluation method, device and medium for offshore engineering service period foundation structure - Google Patents
Stability evaluation method, device and medium for offshore engineering service period foundation structure Download PDFInfo
- Publication number
- CN116882320B CN116882320B CN202311032722.0A CN202311032722A CN116882320B CN 116882320 B CN116882320 B CN 116882320B CN 202311032722 A CN202311032722 A CN 202311032722A CN 116882320 B CN116882320 B CN 116882320B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- fan
- soil body
- offshore engineering
- response
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 20
- 239000002689 soil Substances 0.000 claims abstract description 57
- 230000004044 response Effects 0.000 claims abstract description 39
- 238000005452 bending Methods 0.000 claims abstract description 38
- 238000012546 transfer Methods 0.000 claims abstract description 13
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 26
- 230000009471 action Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000009991 scouring Methods 0.000 claims description 10
- 230000001808 coupling effect Effects 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000007667 floating Methods 0.000 claims description 6
- 238000013097 stability assessment Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000009191 jumping Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 2
- 230000007774 longterm Effects 0.000 abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 15
- 239000013049 sediment Substances 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 13
- 239000012530 fluid Substances 0.000 description 10
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000011010 flushing procedure Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004660 morphological change Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 206010000496 acne Diseases 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses a stability evaluation method, a device and a medium for an offshore engineering service period foundation structure, belonging to the technical field of stability of offshore engineering, comprising the following steps: s1: calculating the response of the offshore engineering to external power; s2: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure; s3: calculating the self rigidity of the soil body; s4: and obtaining the stability evaluation result of the foundation structure in the offshore engineering service period according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body. The invention can effectively and accurately evaluate the conditions of bearing capacity, deformation, structural instability and the like of the fan under the long-term service condition, comprehensively judge the stability of the basic structure in the offshore engineering service period and improve the accuracy of the stability evaluation result.
Description
Technical Field
The invention relates to the technical field of ocean engineering stability, in particular to a stability evaluation method, a device and a medium for an offshore engineering service period foundation structure.
Background
Along with the construction of offshore engineering, the offshore engineering is increasingly built, such as cross-sea bridges, oil platforms, offshore wind farms and the like. Taking an offshore wind power plant as an example, under the long-term service of an offshore wind turbine, ocean power can cause the scouring of a wind turbine foundation and the weakening of soil mass to reduce the bearing capacity, and meanwhile, the lateral bending moment in the rotation process of the wind turbine blade can further aggravate the instability of the wind turbine, so that larger economic loss is caused. However, the process involves fan blade rotation, ocean impact load, soil erosion, six-degree-of-freedom deflection of the airframe structure, and is a complex system with intersecting multidisciplinary aspects such as soil mechanics, ocean dynamics, aerodynamics, structural mechanics, structural kinematics and the like.
At present, the stability research of offshore engineering is concentrated in a single field, such as calculating the lateral bending moment load generated by the rotation of a fan blade, calculating the wave force borne by the fan body, calculating the scouring depth of a fan foundation and the like, and most of related technologies are single generalization (two-dimensional linear generalization or three-dimensional nonlinearity but only single parameter analysis is considered), so that the conduction between all processes and the stability performance of the fan cannot be truly reflected, and no effective method for accurately predicting the bearing capacity, deformation and structural instability of the fan under the long-term service condition exists.
Therefore, how to provide a method, a device and a medium for evaluating the stability of an offshore engineering service period foundation structure is a problem to be solved by those skilled in the art.
Disclosure of Invention
In view of the above, the present invention provides a method, an apparatus and a medium for evaluating stability of an offshore engineering service phase infrastructure, which are used for solving the technical problems in the prior art.
In order to achieve the above object, the present invention provides the following technical solutions:
a stability evaluation method of an offshore engineering service period foundation structure comprises the following steps:
s1: calculating the response of the offshore engineering to external power, wherein the response to external power comprises the response to aerodynamic power, the response to wave power and the response to tide power;
s2: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
s3: calculating the self rigidity of the soil body;
s4: and obtaining the stability evaluation result of the foundation structure in the offshore engineering service period according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body.
Preferably, the S1 includes:
s1.1: calculating the rotating force of the fan under the action of aerodynamic force to obtain the lateral bending moment applied to the fan body;
s1.2: calculating the load force born by the fan body under the action of wave ocean current force;
s1.3: and calculating the scouring depth and range of the seabed of the fan foundation under the action of wave ocean current force.
Preferably, the S2 includes:
s2.1: and calculating the resultant force of wind power, wave current impact load force and pile foundation acting force caused by the bottom vortex after the seabed is flushed.
S2.2: and calculating the lateral comprehensive bending moment of the structure generated by the resultant force of the three components.
Preferably, the step S2 further includes: s2.3, calculating the anti-tilting torque generated by the soil body on the fan under the condition of acting force of the ocean soil on the buried end of the fan.
Preferably, the S3 includes:
s3.1: constructing a space fixed coordinate system by taking an origin o as a reference pointThe internal rigid body translational motion equation is obtainedTotal load applied to rigid body in coordinate system:
in the formula, xi= (xi 1, xi 2, xi 3) t Representing the fixed coordinate system of the o-point on the floating body in spaceDisplacement of (a);
ω＝(ω 1 ，ω 2 ，ω 3 ) t expressed in a coordinate systemLower rigid body angular velocity;
r g ＝(x g ，y g ，z g ) t representing a satellite coordinate systemA lower rigid body barycentric coordinate;
t is a transfer matrix between a satellite coordinate system and a space fixed coordinate system;
is a coordinate system->The total load of the lower rigid body;
s3.2: constructing a rotary motion equation in an oxyz of a satellite coordinate system by taking an origin o as a reference point, and obtaining the total moment of a rigid body relative to the point o under the oxyz coordinate system:
where Io represents the moment of inertia of the rigid body with respect to point o in the coordinate system oxyz; mo represents the total moment of the rigid body with respect to point o in the coordinate system oxyz.
Preferably, the S4 includes:
s4.1: if the self rigidity of the soil body is greater than or equal to the lateral comprehensive bending moment of the structure, jumping to the step S1, repeating the steps S1-S2, increasing the scouring depth of the seabed in the repeated process, and increasing the anti-tilting torque of the soil body on the fan;
s4.2: if the self rigidity of the soil body is smaller than the lateral comprehensive bending moment of the structure, judging that the fan overturns, wherein the lateral bending moment of the machine body is reduced in the overturning process, the fan is gradually stabilized, and if the self rigidity of the soil body is larger than the anti-tilting torque of the soil body to the fan, repeating the step S4.1;
and S4.3, under the coupling effect of ocean wind-wave-current load, obtaining the stability evaluation result of the offshore engineering service period foundation structure according to the process of the steps S4.1-S4.2.
Preferably, the step S4.3 includes: under the coupling effect of ocean wind-wave-current load, according to the process of the steps S4.1-S4.2, when the rigidity of the soil body is enough to support the fan overturning condition, the fan overturning state is stable at the moment, and the overturning angle and the overturning distance are obtained;
if the overturning angle and the overturning distance are both 0, judging that overturning does not occur;
if the overturning angle and/or the overturning distance are/is larger than the corresponding preset value, judging that the fan is completely separated from the soil body and is unstable and collapses.
A stability assessment apparatus for an offshore engineering service phase infrastructure, comprising:
the marine engineering system comprises a first calculation module, a second calculation module and a power generation module, wherein the first calculation module calculates the response of the marine engineering to external power, and the response to the external power comprises the response to aerodynamic power, the response to wave power and the response to tide power;
and an analysis module: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
a second calculation module: calculating the self rigidity of the soil body;
and a judging module: and obtaining the stability evaluation result of the foundation structure in the offshore engineering service period according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body.
A computer readable storage medium storing a computer program which when executed by a processor implements a method of stability assessment of an offshore engineering service infrastructure.
Compared with the prior art, the invention discloses a stability evaluation method, a device and a medium for a basic structure in the service period of offshore engineering, wherein a three-dimensional incompressible Navier-Stokes equation, an effective volume method and a three-dimensional structure unstructured mixed grid are adopted, firstly, the forms of a fan blade and a machine body can be finely generalized, secondly, the operation state of the fan blade can be accurately reflected by a superposition overlapping algorithm, thirdly, the ocean power process such as wave and tide is reflected by a fine capturing VOF algorithm, fourthly, pile foundation flushing depth in the ocean power process is reflected by a coupling sediment transport and soil liquefaction model, thirdly, fan overturning is calculated by a six-degree-of-freedom (6 DOF) algorithm, and thirdly, the action of a soil body exerted on the fan is reflected by a soil body bearing capacity model. The invention can effectively and accurately evaluate the conditions of bearing capacity, deformation, structural instability and the like of the fan under the long-term service condition, and comprehensively judge the stability of the basic structure in the offshore engineering service period.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of the method of the present invention;
FIG. 2 is a schematic view of the volume fraction of the present invention;
fig. 3 is a schematic representation of the superposition of a spatially-fixed coordinate system and a satellite coordinate system according to the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1, the embodiment of the invention discloses a stability evaluation method, a device and a medium of an offshore engineering service period foundation structure, which comprises the following steps:
s1: calculating the response of the offshore engineering to the external power, wherein the response to the external power comprises the response to the aerodynamic power, the response to the wave-current power and the response to the wave-current power;
s2: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
s3: calculating the self rigidity of the soil body;
s4: and obtaining the stability evaluation result of the foundation structure in the offshore engineering service period according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body.
In one embodiment, at S1: the response of the engineering to external power is calculated through the following steps, and the offshore wind turbine is taken as an example and comprises the following steps:
s1.1: calculating the rotating force of the fan under the action of aerodynamic force, and analyzing the lateral bending moment applied to the fan body;
specifically, the air force mainly acts on the fan blade to influence the rotation of the fan blade, and in general, the rotation of the fan blade aggravates the lateral bending moment of the fan, similar to a windmill.
S1.2: and calculating the load force born by the fan body under the action of wave ocean current force.
S1.3: and calculating the scouring depth and range of the seabed of the fan foundation under the action of wave ocean current force.
In one particular embodiment, for wave-fluid power solving, comprising: the basic control equation Navier-Stokes equation of fluid motion is deduced according to the law of conservation of mass and conservation of energy of fluid motion. For the standard Navier-Stokes equation, the continuous equation for incompressible fluid is:
▽·U＝0；
the momentum equation can be deduced from the law of conservation of momentum:
wherein U= (U, v, w) is a velocity vector in a Cartesian coordinate system, and v is a gradient operator;
p is an instantaneous pressure value, g is a gravitational acceleration, v is a kinematic viscosity, and is a relational expression about the kinematic viscosity μ and the density ρ:
specifically, the calculation formula of the gradient operator v is:
wherein x, y and z are 3 directions of a Cartesian coordinate system, and T is a transposition;
specifically, a Volume of fluid (VOF) function is used to capture and track the gas-liquid interface. The VOF algorithm is based on the fluid volume fraction as the lagrangian transport function of the interface, the basic idea is that each grid cell in the computational domain is defined as a fluid volume function α, and the value of α represents the ratio of the volume occupied by the liquid in the grid cell to the volume of the grid cell, and the volume fraction is schematically shown in fig. 2, and the value can be defined by the following formula:
in the classical VOF method, the volume fraction transport equation is defined as:
and correcting the transportation equation after correction as follows:
specifically, by the technical scheme, compared with the classical transport equation, the method has the following formulaThere are more manual compression terms, the last term. The modified VOF method can ensure that equation solution automatically meets the constraint at a water-gas interface, avoids interface blurring and has no influence on an external flow field, wherein U is as follows r Is the relative speed. The free liquid surface form in the interaction process of the wave tides and the wind power pile can be captured more accurately, and the acting force of power on the pile body is reflected accurately.
More specifically, the PIMPLE algorithm is adopted to realize control equation solving, the characteristics of the PISO algorithm (Pressure Implicit Split Operator) and the SIMPLE algorithm (Semi-Implicit Method for Pressure-Linked Equations) are fully combined, namely, the low-relaxation stable calculation is adopted in any time step to obtain the maximum time step (similar to the SIMPLE algorithm), the PISO algorithm is adopted to carry out multiple pressure correction in the iteration step to fully decompose the coupling problem of the pressure speed, and the calculation divergence problem caused by the overlarge change of the flow field in the PISO algorithm is effectively avoided.
In one particular embodiment, for a seabed scour, comprising:
the transportation of coastal sediment and the change of the section of the coastal beach are the results under the coupling action of a plurality of complex power factors, and the intrinsic action mechanism of the coastal sediment needs to be fully understood, and besides the characteristics of the coastal hydrodynamic force, the geological characteristics of the sandy beach, the coastal sediment transportation and the coastal engineering structures need to be deeply known. Early researches are mostly based on shallow water equations and Boussinesq equations, only two-dimensional generalized information averaged along the water depth can be given, water flow characteristics near a sand dam or a flushing pit cannot be fully described, and complex three-dimensional bed surface morphological changes near a structure and corresponding accurate flow field information are difficult to give.
Specifically, the tangential stress of the bed surface is a tangential component of the stress applied to the bed surface, and has two functions, namely, the quantification of the drag force of the bed surface on the unit area and the influence on the convection transportation of the shear flow of the bed surface. The body of water has the general characteristics of a newtonian fluid, namely its stress tensor sigma and strain tensor S ij With linear isotropic relationship between them, constitutive relationship σ=f ij S ij The only form of (2) is:
in the formula, the first term is a stress term caused by pressure, is a state term and has no direct relation with the deformation rate of the fluid, and the second term is a viscous stress caused by the deformation rate of the fluid motion, and is called a deflection stress tensor term. Considering the influence of wall stress tensors, the bed surface shear stress is pointed out as follows:
τ t ＝σ·e n ；
in particular, when the bed load moves near the bed surface, the velocity of the bed load is generally smaller than the local water flow velocity, so that the bed load moves to consume the average energy of the water flow. When the bed surface shear stress is larger than the sediment critical shear stress, the bed load moves, and when the bed load moves in the water body, the bed load is continuously exchanged with bed sand, and when the exchange of the bed load and the bed load is in a dynamic balance, the bed surface form is kept stable, and the bed surface form is as follows:
wherein ρ is the water density, ρ sed Density of sediment, g is gravity acceleration, d 50 Is the median particle diameter of the sediment, theta c Is the critical Hiltz number.
The critical hilz number is a criterion for starting silt particles, and the influence of gradient is considered to indicate:
when the water flow climbs up the slope (flow velocityUp a slope), θ c Will be greater than theta co While the water flow falls back along the slope (the flow rate is downward along the slope), θ c Will be less than theta co Wherein beta is the included angle between the sediment bed surface and the horizontal plane,angle of repose of sediment, θ co The calculation formula of the critical Hiltz number in the flat bottom seabed condition is as follows:
the bed load transport rates in different directions are:
wherein, C is an empirical coefficient, which is generally 1.5-2.3, and 1.5 is taken in the embodiment.
Wave force solving:
solving the lateral bending moment:
M＝F·l/2+M wind power ；
Wherein P is the point pressure, b is the width, H is the water body action range, F is the wave force, l is the length of the fan rod piece, and the comprehensive bending moment of the combined action of wind and wave current is obtained through the step.
M Wind power ＝1/2*P*D 3 *L*V；
Wherein M is Wind power The axial force moment is the lateral moment of the main shaft, P is the wind pressure born by the impeller, D is the diameter of the impeller, L is the length of the impeller of the fan, and V represents the rotating speed of the fan.
Specifically, the transportation of coastal sediment and the change of the section of the beach are the results under the coupling action of a plurality of complex power factors, and the intrinsic action mechanism of the coastal sediment needs to be fully understood, and the geological characteristics of the coastal sediment, the coastal sediment transportation and the coastal engineering structures need to be deeply known. In the prior art, based on a shallow water equation and a Boussinesq equation, only two-dimensional generalized information along the water depth average can be given, water flow characteristics near a sand dam or a flushing pit cannot be fully described, and complex three-dimensional bed surface morphological changes and corresponding accurate flow field information near a structure are difficult to give.
In one particular embodiment, S2: the parameter transfer quantitative analysis is carried out by utilizing the S1 result through the following steps of S2.1-2.3, and the method comprises the following steps:
s2.1: and calculating the resultant force of wind power, wave current impact load force and pile foundation acting force caused by the bottom vortex after the seabed is scoured.
S2.2: and S2.1, calculating the lateral comprehensive bending moment of the structure generated by the resultant force of the action.
S2.3: and calculating the anti-tilting torque generated by the soil body in the step S1.4 on the fan.
In one particular embodiment, S3: calculating the self stiffness of the soil body comprises: the six degrees of freedom calculation is as follows:
six degrees of freedom motion is performed in the process of overturning the fan, and two coordinate systems, namely a space fixed coordinate system, are used in the process of deriving a rigid motion equationAnd a satellite coordinate system oxyz. The satellite coordinate system origin o is located at the floating body gravity center position. When the float is in the initial state, the two coordinate systems coincide, see fig. 3. The incident wave is represented in a third spatially-fixed coordinate system ozz in which the OXY-plane is located on the still water surface and the Z-axis is vertically upward. Rigid body translational motion equation is in space fixed coordinate system +.>The inner expression, the revolving equation of motion is expressed in the random coordinate system oxyz, and the origin o is taken as a reference point. The equation of motion is shown in the following formula,
in the formula, xi= (xi 1, xi 2, xi 3) t Representing the fixed coordinate system of the o-point on the floating body in spaceDisplacement of (a);
representing the o-point on the floating body in a spatially fixed coordinate system +.>Acceleration in (a);
ω＝(ω 1 ，ω 2 ，ω 3 ) t representing the rigid body angular velocity in the coordinate system oxyz;
r g ＝(x g ，y g ，z g ) t representing the gravity center coordinates of the rigid body under the oxyz of the satellite coordinate system;
io represents the moment of inertia of the rigid body with respect to the point o in the coordinate system oxyz;
is a coordinate system->The total load of the lower rigid body;
M o representing the total moment of the rigid body relative to the point o in the coordinate system oxyz;
t is a transfer matrix between the satellite coordinate system and the spatially-fixed coordinate system.
In one particular embodiment, S4: the method for judging the rigidity of the soil body and the magnitude of the lateral bending moment born by the soil body by further carrying out parameter transfer analysis by utilizing the analysis result of the step S2 and the calculation result of the step S3.1-S3.3 comprises the following steps:
s4.1: if the self rigidity of the soil body is greater than or equal to the lateral comprehensive bending moment of the structure, jumping to the step S1, repeating the steps S1-S2, increasing the scouring depth of the seabed in the repeated process, and increasing the anti-tilting torque of the soil body on the fan;
s4.2: if the self rigidity of the soil body is smaller than the lateral comprehensive bending moment of the structure, the fan is judged to be overturned, the lateral bending moment of the machine body is reduced in the process of overturning, the fan is gradually stabilized, and if the self rigidity of the soil body is larger than the anti-tilting torque of the soil body to the fan, the step S4.1 is repeated;
and S4.3, under the coupling effect of ocean wind-wave-current load, obtaining the stability evaluation result of the offshore engineering service period foundation structure according to the process of the steps S4.1-S4.2.
Specifically, according to S4.1-S4.2, under the action of ocean wind-wave-current load coupling, the process can be carried out until the rigidity of soil is enough to support the fan overturning condition, and the fan overturning state is stable at the moment, and the overturning angle and the overturning distance can be given at the moment. If no overturning occurs, the overturning angle and the overturning distance are both 0, and the fan can be completely separated from the soil body in the process, so that the fan is unstable and collapses, and the overturning angle and the overturning distance are larger. For example: when the overturning angle and the overturning distance are close to 90 degrees and/or the overturning distance exceeds one time, the top of the fan is at a height from the seabed, the fan is completely separated from the soil body, and the fan is unstable and collapses.
A stability assessment apparatus for an offshore engineering service phase infrastructure, comprising:
the first calculation module is used for calculating the response of the offshore engineering to the external power, wherein the response to the external power comprises the response to the aerodynamic power, the response to the wave power and the response to the tide power;
and an analysis module: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
a second calculation module: calculating the self rigidity of the soil body;
and a judging module: and obtaining the stability evaluation result of the foundation structure in the offshore engineering service period according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body.
A computer readable storage medium storing a computer program which when executed by a processor implements a method of stability assessment of an offshore engineering service infrastructure.
Through the technical scheme, the three-dimensional incompressible Navier-Stokes equation, the effective volume method and the three-dimensional structure unstructured mixed grid are adopted, firstly, the fan blade and the machine body form can be finely generalized, secondly, the operation state of the fan blade is accurately reflected by a superposition overlapping algorithm, thirdly, the ocean power processes such as wave and tide are reflected by a fine capturing VOF algorithm, fourthly, pile foundation flushing depth in the ocean power process is reflected by a coupling sediment transportation and soil liquefaction model, thirdly, fan overturning is calculated by a six-degree-of-freedom (6 DOF) algorithm, and thirdly, the action of soil on the fan is reflected by a soil bearing capacity model.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (7)
1. The stability evaluation method of the offshore engineering service period foundation structure is characterized by comprising the following steps of:
s1: calculating the response of the offshore engineering to external power, wherein the response to external power comprises the response to aerodynamic power, the response to wave power and the response to tide power;
s2: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
s3: calculating the self stiffness of the soil body, comprising:
s3.1: constructing a space fixed coordinate system by taking an origin o as a reference pointThe internal rigid body translational motion equation is obtainedTotal load applied to rigid body in coordinate system:
in the formula, xi= (xi 1, xi 2, xi 3) t Representing the fixed coordinate system of the o-point on the floating body in spaceDisplacement of (a);
ω＝(ω 1 ，ω 2 ，ω 3 ) t expressed in a coordinate systemLower rigid body angular velocity;
r g ＝(x g ，y g ，z g ) t representing a satellite coordinate systemA lower rigid body barycentric coordinate;
t is a transfer matrix between a satellite coordinate system and a space fixed coordinate system;
is a coordinate system->The total load of the lower rigid body;
s3.2: constructing a rotary motion equation in an oxyz of a satellite coordinate system by taking an origin o as a reference point, and obtaining the total moment of a rigid body relative to the point o under the oxyz coordinate system:
where Io represents the moment of inertia of the rigid body with respect to point o in the coordinate system oxyz; mo represents the total moment of the rigid body with respect to point o in the coordinate system oxyz;
s4: according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body, a stability evaluation result of the foundation structure in the offshore engineering service period is obtained, and the method comprises the following steps:
s4.1: if the self rigidity of the soil body is greater than or equal to the lateral comprehensive bending moment of the structure, jumping to the step S1, repeating the steps S1-S2, increasing the scouring depth of the seabed in the repeated process, and increasing the anti-tilting torque of the soil body on the fan;
s4.2: if the self rigidity of the soil body is smaller than the lateral comprehensive bending moment of the structure, judging that the fan overturns, wherein the lateral bending moment of the machine body is reduced in the overturning process, the fan is gradually stabilized, and if the self rigidity of the soil body is larger than the anti-tilting torque of the soil body to the fan, repeating the step S4.1;
and S4.3, under the coupling effect of ocean wind-wave-current load, obtaining the stability evaluation result of the offshore engineering service period foundation structure according to the process of the steps S4.1-S4.2.
2. The method for evaluating the stability of an offshore engineering service foundation structure according to claim 1, wherein S1 comprises:
s1.1: calculating the rotating force of the fan under the action of aerodynamic force to obtain the lateral bending moment applied to the fan body;
s1.2: calculating the load force born by the fan body under the action of wave ocean current force;
s1.3: and calculating the scouring depth and range of the seabed of the fan foundation under the action of wave ocean current force.
3. The method for evaluating the stability of an offshore engineering service foundation structure according to claim 1, wherein S2 comprises:
s2.1: calculating the resultant force of wind power, wave current impact load force and pile foundation acting force caused by bottom vortex after scouring of the seabed;
s2.2: and calculating the lateral comprehensive bending moment of the structure generated by the resultant force of the three components.
4. The method for evaluating the stability of an offshore engineering service foundation structure according to claim 1, wherein the step S2 further comprises: s2.3, calculating the anti-tilting torque generated by the soil body on the fan under the condition of acting force of the ocean soil on the buried end of the fan.
5. The method for evaluating the stability of an offshore engineering service foundation structure according to claim 1, wherein S4.3 comprises: under the coupling effect of ocean wind-wave-current load, according to the process of the steps S4.1-S4.2, when the rigidity of the soil body is enough to support the fan overturning condition, the fan overturning state is stable at the moment, and the overturning angle and the overturning distance are obtained;
if the overturning angle and the overturning distance are both 0, judging that overturning does not occur;
if the overturning angle and/or the overturning distance are/is larger than the corresponding preset value, judging that the fan is completely separated from the soil body and is unstable and collapses.
6. A stability assessment apparatus for an offshore engineering service phase foundation structure, comprising:
the marine engineering system comprises a first calculation module, a second calculation module and a power generation module, wherein the first calculation module calculates the response of the marine engineering to external power, and the response to the external power comprises the response to aerodynamic power, the response to wave power and the response to tide power;
and an analysis module: carrying out parameter transfer quantitative analysis based on each response result to obtain a lateral comprehensive bending moment of the structure;
a second calculation module: calculating the self stiffness of the soil body, comprising:
s3.1: constructing a space fixed coordinate system by taking an origin o as a reference pointThe internal rigid body translational motion equation is obtainedTotal load applied to rigid body in coordinate system:
in the formula, xi= (xi 1, xi 2, xi 3) t Representing the fixed coordinate system of the o-point on the floating body in spaceDisplacement of (a);
ω＝(ω 1 ，ω 2 ，ω 3 ) t expressed in a coordinate systemLower rigid body angular velocity;
r g ＝(x g ，y g ，z g ) t representing a satellite coordinate systemA lower rigid body barycentric coordinate;
t is a transfer matrix between a satellite coordinate system and a space fixed coordinate system;
is a coordinate system->The total load of the lower rigid body;
s3.2: constructing a rotary motion equation in an oxyz of a satellite coordinate system by taking an origin o as a reference point, and obtaining the total moment of a rigid body relative to the point o under the oxyz coordinate system:
where Io represents the moment of inertia of the rigid body with respect to point o in the coordinate system oxyz; mo represents the total moment of the rigid body with respect to point o in the coordinate system oxyz;
and a judging module: according to the lateral comprehensive bending moment of the structure and the self rigidity of the soil body, a stability evaluation result of the foundation structure in the offshore engineering service period is obtained, and the method comprises the following steps:
s4.1: if the self rigidity of the soil body is greater than or equal to the lateral comprehensive bending moment of the structure, jumping to the step S1, repeating the steps S1-S2, increasing the scouring depth of the seabed in the repeated process, and increasing the anti-tilting torque of the soil body on the fan;
s4.2: if the self rigidity of the soil body is smaller than the lateral comprehensive bending moment of the structure, judging that the fan overturns, wherein the lateral bending moment of the machine body is reduced in the overturning process, the fan is gradually stabilized, and if the self rigidity of the soil body is larger than the anti-tilting torque of the soil body to the fan, repeating the step S4.1;
and S4.3, under the coupling effect of ocean wind-wave-current load, obtaining the stability evaluation result of the offshore engineering service period foundation structure according to the process of the steps S4.1-S4.2.
7. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements a method of evaluating the stability of an offshore engineering service infrastructure according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311032722.0A CN116882320B (en) | 2023-08-16 | 2023-08-16 | Stability evaluation method, device and medium for offshore engineering service period foundation structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311032722.0A CN116882320B (en) | 2023-08-16 | 2023-08-16 | Stability evaluation method, device and medium for offshore engineering service period foundation structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116882320A CN116882320A (en) | 2023-10-13 |
| CN116882320B true CN116882320B (en) | 2024-03-22 |
Family
ID=88268315
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311032722.0A Active CN116882320B (en) | 2023-08-16 | 2023-08-16 | Stability evaluation method, device and medium for offshore engineering service period foundation structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116882320B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115048708A (en) * | 2022-08-10 | 2022-09-13 | 中国长江三峡集团有限公司 | Vulnerability assessment method and device for single-pile foundation of offshore wind turbine |
| WO2023087601A1 (en) * | 2021-11-22 | 2023-05-25 | 中交天津港湾工程研究院有限公司 | Stability prediction method for deep water thin-walled steel cylinder |
-
2023
- 2023-08-16 CN CN202311032722.0A patent/CN116882320B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023087601A1 (en) * | 2021-11-22 | 2023-05-25 | 中交天津港湾工程研究院有限公司 | Stability prediction method for deep water thin-walled steel cylinder |
| CN115048708A (en) * | 2022-08-10 | 2022-09-13 | 中国长江三峡集团有限公司 | Vulnerability assessment method and device for single-pile foundation of offshore wind turbine |
Non-Patent Citations (5)
| Title |
|---|
| Stability prediction based on individual impedance measurement for distributed DC power systems;Hao Wang et al.;《IEEE Xplore》;20110707;全文 * |
| 导管架海洋平台整体结构海冰载荷环境下响应分析;邱少华;船海工程(第04期);全文 * |
| 海上风机复合单桩基础水平承载力数值分析;孔德森 等;科学技术与工程(第20期);全文 * |
| 海洋高桩基础水平单调及循环加载现场试验;朱斌 等;岩土工程学报(第06期);全文 * |
| 风浪流共同作用下海上风电基础与海床的动力响应分析;张小玲 等;防灾减灾工程学报(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116882320A (en) | 2023-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jung et al. | Effect of monopile foundation modeling on the structural response of a 5-MW offshore wind turbine tower | |
| Chen et al. | Static and dynamic loading behavior of a hybrid foundation for offshore wind turbines | |
| Coudurier et al. | Modelling of a tuned liquid multi-column damper. Application to floating wind turbine for improved robustness against wave incidence | |
| Wan et al. | A combined wind and wave energy-converter concept in survival mode: Numerical and experimental study in regular waves with a focus on water entry and exit | |
| Leble et al. | Demonstration of a coupled floating offshore wind turbine analysis with high-fidelity methods | |
| Yang et al. | Experimental and numerical studies on the low-frequency responses of a spar-type floating offshore wind turbine | |
| de Souza et al. | Design, structural modeling, control, and performance of 20 MW spar floating wind turbines | |
| Ding et al. | Numerical analysis of passive toroidal tuned liquid column dampers for the vibration control of monopile wind turbines using FVM and FEM | |
| Huang et al. | Hydrodynamic study and performance analysis of the OC4-DeepCwind platform by CFD method | |
| Chen et al. | Nonlinear wave loads on high-rise pile cap structures in the Donghai Bridge Wind Farm | |
| Utsunomiya et al. | Numerical modeling and analysis of a hybrid-spar floating wind turbine | |
| Yu et al. | Numerical modelling of wave run-up heights and loads on heaving buoy wave energy converter under the influence of regular waves | |
| Zhang et al. | An aero-hydro coupled method for investigating ship collision against a floating offshore wind turbine | |
| Connolly et al. | Fully coupled aero-hydro-structural simulation of new floating wind turbine concept | |
| Cheng et al. | Response of offshore wind turbine tripod suction bucket foundation to seismic and environmental loading | |
| Vijayan et al. | Shock transmission through universal joint of cutter suction dredger | |
| Zeng et al. | Nonlinear hydrodynamics of floating offshore wind turbines: A review | |
| Shi et al. | Aerodynamic performance of semi-submersible floating wind turbine under pitch motion | |
| Cao et al. | Experimental investigation on the dynamic response of an innovative semi-submersible floating wind turbine with aquaculture cages | |
| CN116882320B (en) | Stability evaluation method, device and medium for offshore engineering service period foundation structure | |
| Chen et al. | CCP-WSI blind test series 3: OpenFOAM simulation of focused wave interaction with a simplified wave energy converter | |
| Nevalainen | The effect of unsteady sea conditions on tidal stream turbine loads and durability | |
| Ko et al. | Numerical study for nonlinear hydrodynamic coefficients of an asymmetric wave energy converter | |
| Aliyar et al. | Experimental investigation of offshore crane load during installation of a wind turbine jacket substructure in regular waves | |
| Huang et al. | Wake interaction between two floating offshore wind turbines with blade deformation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |