CN113761775A - SPH simulation display method and system for overall process of disaster evolution of submarine landslide - Google Patents
SPH simulation display method and system for overall process of disaster evolution of submarine landslide Download PDFInfo
- Publication number
- CN113761775A CN113761775A CN202110863595.3A CN202110863595A CN113761775A CN 113761775 A CN113761775 A CN 113761775A CN 202110863595 A CN202110863595 A CN 202110863595A CN 113761775 A CN113761775 A CN 113761775A
- Authority
- CN
- China
- Prior art keywords
- particles
- slope
- particle
- landslide
- slope body
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 108
- 230000008569 process Effects 0.000 title claims abstract description 50
- 238000004088 simulation Methods 0.000 title claims abstract description 31
- 239000002245 particle Substances 0.000 claims abstract description 186
- 238000004364 calculation method Methods 0.000 claims abstract description 57
- 239000013535 sea water Substances 0.000 claims abstract description 43
- 230000009471 action Effects 0.000 claims abstract description 13
- 238000006073 displacement reaction Methods 0.000 claims abstract description 10
- 230000005484 gravity Effects 0.000 claims abstract description 8
- 230000000737 periodic effect Effects 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 20
- 239000000463 material Substances 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 10
- 239000007787 solid Substances 0.000 claims description 7
- 239000007790 solid phase Substances 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 2
- 230000005012 migration Effects 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 239000013618 particulate matter Substances 0.000 claims description 2
- 238000012876 topography Methods 0.000 claims description 2
- 230000005433 particle physics related processes and functions Effects 0.000 claims 1
- 239000012071 phase Substances 0.000 description 6
- 230000003993 interaction Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 229920000034 Plastomer Polymers 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000011343 solid material Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
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/25—Design optimisation, verification or simulation using particle-based methods
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention discloses an SPH simulation display method and system for the whole process of submarine landslide catastrophe evolution, which comprises the following steps: dispersing the calculation areas of the seabed slope body and the seawater into a set number of particles with different physical state attributes, and giving corresponding physical mechanical parameters to the physical state particles; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary; applying the infinitely-covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and the structural acting force applied to a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and simultaneously applying a displacement boundary condition; in each calculation time step, determining a catastrophe evolution stage of the seabed slope body, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step, and completing updating of particle physical mechanics information of solution domain; and displaying the physical and mechanical information of the particles at the set time step.
Description
Technical Field
The invention relates to the technical field of disaster simulation of a submarine landslide, in particular to an SPH simulation display method and system for the whole process of disaster evolution of the submarine landslide.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
Smooth Particle Hydrodynamics (smooth Particle dynamics) is a meshless method that describes the state of a system through a series of particles, which is an adaptive lagrange Particle method that simulates fluid flow. Compared with the traditional grid-based numerical method, the SPH method has self-adaptability, so that the formula structure of the SPH is not influenced by the randomness of particle distribution, and the problem of great deformation can be naturally treated. The SPH method takes particles as a computational framework for field variable approximation. In the calculation process, the SPH approximation does not need to use a grid defined in advance to provide connection information among particles, and even when no particles are subjected to thinning operation, the SPH approximation method still works well, so that the difficulty in the traditional FEM and FDM solving process is well avoided. Besides the non-grid property and the self-adaptive property, the SPH method has the characteristic that the Lagrange formula and the particle approximation method can be harmoniously combined. Compared with other non-grid methods in which non-grid nodes are only used as interpolation points, the particles in the SPH method simultaneously carry material properties and can move under the interaction of external force and internal force. The SPH particles have the function of both being an approximation point and a material composition, making the SPH particles more flexible. SPH has been shown to give stable numerical results by solving many problems with extreme deformation using arbitrarily distributed points. Although SPH is currently widely used in the fields of astrophysics, hydrodynamics, classical mechanics, and the like. However, the large-scale calculation research for ocean engineering, particularly offshore seabed slope engineering and the like is less.
Due to the complexity of the marine environment and the wireless coverage of seawater, when a landslide disaster occurs on the seabed, the landslide disaster cannot be detected by people unless tsunami is accompanied, and the landslide disaster is difficult to detect in many cases. The general submarine landslide has the advantages of long duration, long sliding distance, large carrying volume and strong destructiveness in the sliding process, and can transform submarine landforms and pose great threats to the safety of submarine engineering sites and underwater structures. The control of the catastrophe evolution process of the seabed landslide is beneficial to the stability evaluation of seabed artificial fields and the slope of the reclamation islands, and is also beneficial to the smoothness of seabed communication optical cables, petroleum engineering facilities and seabed channels.
The method is obviously different from land landslide in that in the process of occurrence of the submarine landslide, a landslide body interacts with seawater, the proportion of water and a slope body material is changed at any time, the physical and mechanical properties and the physical state of the landslide body are changed, and the catastrophe evolution mechanism of the whole submarine landslide process is difficult to effectively reveal only by adopting single slope body stability calculation or slip transport mechanism analysis. In the prior art, no effective measure is available for describing and displaying the evolution characteristics of a slope body in the catastrophe process of offshore submarine landslide.
Disclosure of Invention
In view of the above, the invention provides an SPH simulation display method and system for the whole process of the catastrophic evolution of the submarine landslide, which are used for respectively representing seawater and slope body materials (fluid phase and solid phase) in the submarine landslide process by SPH particles with different attributes, are easy to track the fluid surface and reconstruct the fluid surface, and can realize the visual display of the details and results of the mixed fluid of the landslide body and the seawater.
In order to achieve the above purpose, in some embodiments, the following technical solutions are adopted:
an SPH simulation display method for the whole process of submarine landslide catastrophe evolution comprises the following steps:
dispersing the calculation areas of the seabed slope body and the seawater into a set number of particles with different physical state attributes, and giving corresponding physical mechanical parameters to the physical state particles; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary;
applying the infinitely-covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and the structural acting force applied to a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and simultaneously applying a displacement boundary condition;
in each calculation time step, determining a catastrophe evolution stage of the seabed slope body, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step, and completing updating of particle physical mechanics information of solution domain;
and displaying the physical and mechanical information of the particles at the set time step.
In other embodiments, the following technical solutions are adopted:
an SPH simulation display system of the whole process of submarine landslide catastrophe evolution comprises:
the ionization module is used for dispersing the calculation area of the seabed slope body and the seawater into particles with set quantity and different physical state attributes, and endowing the particles with the corresponding physical and mechanical parameters to each physical state; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary;
the acting force applying module is used for applying the infinitely covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and structural acting force borne by a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and applying a displacement boundary condition;
the particle information updating module is used for determining a catastrophe evolution stage of the seabed slope body in each calculation time step, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step and completing particle physical and mechanical information updating of a solution domain;
and the catastrophe evolution overall process display module is used for displaying the particle physical and mechanical information of the set time step.
In other embodiments, the following technical solutions are adopted:
a terminal device comprising a processor and a computer-readable storage medium, the processor being configured to implement instructions; the computer readable storage medium is used for storing a plurality of instructions, and the instructions are suitable for being loaded by the processor and executing the SPH simulation display method of the overall process of the catastrophic evolution of the seabed landslide.
In other embodiments, the following technical solutions are adopted:
a computer readable storage medium, wherein a plurality of instructions are stored, and the instructions are suitable for being loaded by a processor of a terminal device and executing the SPH simulation display method of the whole process of the catastrophic evolution of the seabed landslide.
Compared with the prior art, the invention has the beneficial effects that:
(1) one or more embodiments of the invention introduce the core thought and the calculation advantages of smooth particle hydrodynamics into the field of ocean engineering, can capture the sliding accumulation characteristics of variable submarine landslide bodies and the complex sea wave action and surge forms, can realize effective simulation and visual display of the whole process of the offshore submarine landslide catastrophe evolution, and has higher calculation efficiency and precision.
(2) One or more embodiments of the invention can track the motion characteristics of the landslide body particles and the seawater particles in real time in the whole solving process, do not need equivalent calculation in a numerical method among different physical state particles, avoid complex coupling information interaction operation in the process of adopting different solving methods in different landslide stages, effectively ensure high-efficiency solving of a problem domain, and reproduce the gliding motion mode of the slope body, the debris flow carrying mechanism and the accumulation characteristics of the slope body materials in real time.
(3) One or more embodiments of the invention model the whole solution area, can simultaneously reproduce the sliding accumulation characteristics of a landslide body and the surge influence of seawater in the landslide process, can accurately track the wave height and action range of the sea wave, and can provide theoretical guidance for the safety protection of offshore marine engineering.
Advantages of additional aspects of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
FIG. 1 is a flow chart of an SPH simulation display method of the overall process of the catastrophic evolution of the seabed landslide in the embodiment of the invention;
FIG. 2 is a schematic diagram of a particle approximation method according to an embodiment of the present invention.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
The embodiments and features of the embodiments of the present invention may be combined with each other without conflict.
The method has obvious stage characteristics in the occurrence process of the submarine landslide disaster, such as motion stages of integral sliding, debris flow, turbidity current and the like. In the above process, the problem of ultra-large deformation is often accompanied, and the grid-based numerical method is no longer applicable. As a representative of the gridless method, the SPH method has been successfully applied to the simulation of large deformation of solid materials and complicated flow processes of mixed fluids. The sea water and slope body materials (fluid phase and solid phase) in the process of submarine landslide are respectively represented by SPH particles with different attributes, the fluid particles and the solid particles are interacted in a discrete form based on respective constitutive relations, and the coupling effect of the fluid particles and the solid particles only needs to ensure that the two conditions of no sliding at the boundary and consistent stress are met. Meanwhile, the SPH method can simulate fluid to obtain particle data results of calculation domains, is easier to track the fluid surface and reconstruct the fluid surface, is convenient for analyzing the mixed fluid details of the landslide body and the seawater and visualizing the results, and is an effective numerical simulation method for reproducing the motion process of the submarine landslide body.
Example one
Based on the description, in one or more embodiments, an SPH simulation display method for the whole process of the disaster evolution of the submarine landslide is disclosed, firstly, a calculation area is discretized into a certain number of SPH particles with different physical state attributes without grids, and the particles are subjected to assignment of physical and mechanical information, quality, stress, position, speed, period, amplitude and other physical and mechanical parameters; a non-slip boundary layer combining repulsive particles and virtual particles is arranged at the boundary, so that the influence of the boundary effect on a core calculation area is compensated; applying the gravity action and the wave action of the seawater on the seawater particles, converting the land slope acting force and the structural acting force on the slope body into a stress boundary condition of a slope body particle calculation area, and applying the stress boundary condition on the boundary layer particles, and simultaneously applying a displacement boundary condition; selecting a proper kernel function and a neighborhood solving radius to form a neighboring particle tight branch domain of all particles in a calculation domain; in each solving time step, performing particle traversal search in a tight branch region, judging the interaction relation between a slope material and seawater, obtaining the ratio of liquid phase and solid phase particle contents of research particles in a particle system, further determining a catastrophe evolution stage of the seabed slope, and simultaneously selecting a system constitutive equation set of a corresponding stage; solving a continuity equation, a momentum equation and an energy equation of the particles in the tight branch region by adopting a display integration method, and performing density summation, stress solution and energy change rate calculation of the particles in the region so as to complete the update of the information of the particles in the region in the solution within a single time step; and setting a time step interval for outputting the calculation result, and finishing post-processing display of the calculation result by adopting a post-processing program or software.
Referring to fig. 1, the SPH simulation display method for the overall process of the catastrophic evolution of the seafloor landslide of this embodiment specifically includes the following processes:
(1) the method comprises the steps of dispersing a submarine slope calculation area without grids into a certain number of particles with different physical attributes, dividing a submarine slope into units, creating slope model boundary subunits by combining slope surface topography and slope geometric boundaries, performing internal interpolation by means of slope profiles, uniformly inserting the model subunits, constructing unit groups of a whole slope model, selecting unit center points as coordinates of slope SPH particles, and simultaneously giving physical and mechanical parameters such as quality, density, stress, position, speed, motion period and the like to the SPH particles;
and similarly, creating seawater SPH particles.
(2) The layer number of the boundary layer particles can be increased, the boundary layer particles are laid in a staggered mode, and the repulsive particles are increased to prevent particles caused by excessive speed of real particles from escaping; and giving a certain speed to the virtual particles to realize the effect of no speed slip on the boundary, and creating a speed non-slip boundary particle layer by combining the boundary repulsive particles and the virtual particles.
In this embodiment, the speed non-slip boundary particle layer is set to weaken the influence of the boundary effect on the calculation region and ensure the accuracy of the calculation result of the calculation region. In the motion simulation after the instability of the seabed side slope, the bottom of the solid particle is often accelerated greatly due to excessive stress, and further the speed is excessive so as to escape from the boundary. By combining the reverse repulsive force of the repulsive particles and the initial velocity value of the virtual particles, the boundary slippage of the non-physical penetration boundary of the particles and the virtual particles can be prevented, the precision of the SPH approximation method at the boundary area is improved, and the accuracy of the result of the boundary simulation of the calculation area is effectively ensured. Boundary layer particles are not actually calculated regions and are only used for compensating and correcting the solution of the actual particles on the boundary of the calculation region.
(3) The gravity action of the seawater and the wave load are applied to the seawater particles, the land slope acting force and the structural acting force borne by the slope body are converted into the stress boundary condition of the slope body particle calculation area, and meanwhile, the constraint is converted into the displacement boundary condition.
The seabed side slope is usually in an infinite seawater coverage environment, the seawater depth also changes along with the extension of the side slope, the calculation cost and the calculation efficiency are limited, and the whole sea area cannot be subjected to discretization modeling, so that only a potential landslide area and a seawater layer with a certain height on a slope body are subjected to non-grid discretization particle modeling. Wave action, in addition to ocean bottom earthquakes, ocean bottom volcanic activity, often induces major factors in ocean bottom landslide. The SPH method has inherent advantages in capturing and tracking interfaces of different phase substances (different object states), so that loading of different waveform wave actions can be performed.
(4) And selecting a proper smooth kernel function and a particle neighborhood searching radius, performing particle pairing by adopting a linked list searching method, forming support domains and influence domains of all SPH particles, and determining the interaction relationship among the particles.
Integrating an arbitrary function and a smooth kernel function by applying an integral expression function method, replacing the kernel function delta (x-x ') of the delta function by the smooth function W (x-x', h), and expressing the integral expression of f (x) as follows:
f(x)≈∫Ωf(x′)W(x-x′，h)dx′；
w (x-x', h) is referred to in the SPH method as a smooth kernel, also known as a kernel. W (x-x', h) is defined in the tight branch, and in the smooth function, h is the smooth length defining the area of influence of the smooth function W. Typically, the support domain and the influence domain of a particle are the same.
The particle approximation of the function at particle i can be finally expressed as:
wherein, Wij＝W(xi-xjH), any function value at particle i can be approximated by applying a smooth function to the weighted average of the function values of all particles in its immediate branch.
(5) And in each calculation time step, traversing search in a particle tight branch region is carried out, the relative speed among the particles of the slope body and the content ratio of the particles of the sea water and the particles of the slope body in the tight branch region are obtained, the interaction relation between the slope body material and the sea water is judged, and the catastrophe evolution stage of the seabed slope body is further determined.
(5-1) the landslide process of the seabed is often accompanied with obvious stage characteristics, such as stages of overall slippage, debris flow, mud flow, turbidity flow and the like. Particulate matter concentration, volume concentration and landslide through seabed landslide body materialIndoor test of body material viscosity, calibrating seawater particle content ratio (liquid phase alpha) in particle system in different cataclysm evolution stages1) Slope particle content ratio (solid phase alpha)2) And the ratio beta of the contents of the two-phase particlesi；α1、α2After certain conversion, the concentration of the particulate matters is corresponded; beta is aiAfter certain conversion, the volume concentration is corresponding to the concentration; wherein the ratio of the contents of the two-phase particles is betaiThe specific calculation is as follows:
in this embodiment, β is the critical slip instability of the whole slopei＝β1Stage beta of slip mass debris flowi＝β2Slip mass mud flow phase betai＝β3Stage of turbulent flow of the landslide body betai＝β4。
(6) And selecting a constitutive equation of the particle system of the corresponding stage according to the catastrophe stage judgment in the particle tight branch region. In the stage of the overall movement and slippage of the landslide body, the relative movement speed between solid-phase particles of the landslide body is close to zero; aiming at offshore seabed sediment, field tests and indoor tests can determine that the slope material is likely to change from a solid state to a liquid state only when the water content of the slope reaches a certain value, relative displacement occurs between the slope materials, and the migration of a slip mass is in a debris flow mode. With the continuous sliding of the landslide body, the seawater interacts with the landslide body, and when the solid particle volume of the landslide body exceeds a certain value, the landslide body slides to present the mud flow characteristic of the non-Newtonian fluid. As the slope continues to slide, the water content of the slope material further increases and the slip of the slope will exhibit a turbulent flow pattern.
(6-1) for the overall slippage stage of the slope body, adopting Drucker-Prager yield criterion to describe the constitutive relation and the yield state of the slope body, wherein the yield condition can be expressed as follows:
wherein, I1And j2Being the first and second invariant of the stress tensor, αφAnd kcIs a constant and can be obtained by the cohesive force c and the internal friction angle phi of the slope material.
For plane strain problems, the constant αφAnd kcIs composed of
(6-2) aiming at the debris flow stage of the landslide body, the mixed high-density fluid of the slope body material and the seawater has obvious non-Newtonian fluid characteristics, a Bingham model (Bingham model) is adopted to represent the rheological characteristics of the debris flow,
ultimate shear stress τ0Also known as yield stress. The rheological parameter η is called plastic viscosity and is used to characterize the magnitude of viscosity of a material when flowing.
(6-3) aiming at the mud flow stage of the landslide body, the solid phase slope body particle content concentration is larger than that of the conventional sea current and turbidity current, the mud flow can be regarded as a non-Newtonian fluid with viscoplasticity characteristics, and a Herschel-Bulkly generalized Bingham plastomer model is adopted:
where k is the consistency coefficient, related to the fluid properties, in Pa · sn(ii) a n is a fluidity index and is dimensionless; tau is0In order to limit the dynamic shear stress,
in order to obtain a shear deformation speed,
the model integrates the characteristics of the Bingham model and the power law model.
(6-4) aiming at the turbid flow stage of the sliding mass, the sliding mass is gradually converted into a sand-carrying water body and has the characteristics of Newtonian fluid, the SPH water-sand coupling model is adopted to process the sliding mass at the stage, the essence of the sliding mass is consistent with that of the traditional SPH hydrodynamic model, and the momentum equation in the tensor form is as follows:
where u is the particle velocity, ρ is the particle density, g is the gravitational acceleration, and σ is the full stress tensor.
(7) In each time step, the particle density summation in the tight branch region is carried out by adopting a display integration method, the adjacent particles of the particles are searched, the particle density information at the moment is solved, and then the stress solution and the momentum change rate calculation are carried out by a system state equation; and then updating the particle physical and mechanical information of the solution domain, and entering the next calculation time step cycle after updating the particle system information.
In each time step, with reference to fig. 2, solving an SPH approximation of a continuity equation of the particle system by a kernel approximation and a particle approximation method to obtain a particle density change rate in a particle tight-branching domain; solving the stress change of the particle system according to the SPH approximate expression of the momentum equation; and obtaining the momentum change rate of the particle system by the SPH approximate expression of the energy equation. After the solving of the particle system equation is completed, particle information is updated, wherein the particle information mainly comprises the coordinates, stress, acceleration, momentum change rate and the like of the particles, and then the next time step is calculated.
(8) And setting the output time step interval of the particle system solving result, outputting the physical and mechanical information of the particle system after solving the time step at certain intervals, and carrying out batch processing on the output result through a post-processing program or software, so that the whole process display of the submarine landslide catastrophe evolution can be carried out in real time.
It should be noted that the process of performing the visualization is implemented by those skilled in the art according to the prior art, and is not described in detail herein.
Example two
In one or more embodiments, an SPH simulation display system for the overall process of catastrophic evolution of landslide is disclosed, comprising:
the ionization module is used for dispersing the calculation area of the seabed slope body and the seawater into particles with set quantity and different physical state attributes, and endowing the particles with the corresponding physical and mechanical parameters to each physical state; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary;
the acting force applying module is used for applying the infinitely covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and structural acting force borne by a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and applying a displacement boundary condition;
the particle information updating module is used for determining a catastrophe evolution stage of the seabed slope body in each calculation time step, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step and completing particle physical and mechanical information updating of a solution domain;
and the catastrophe evolution overall process display module is used for displaying the particle physical and mechanical information of the set time step.
It should be noted that specific implementation manners of the modules are already described in detail in the first embodiment, and are not described herein again.
EXAMPLE III
In one or more embodiments, a terminal device is disclosed, which includes a server, where the server includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to implement the SPH simulation display method of the overall process of the catastrophic evolution of the sea-bottom landslide in the first embodiment. For brevity, no further description is provided herein.
It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general purpose processors, digital signal processors DSP, application specific integrated circuits ASIC, off-the-shelf programmable gate arrays FPGA or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and so on. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include both read-only memory and random access memory, and may provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory. For example, the memory may also store device type information.
In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or instructions in the form of software.
Example four
In one or more embodiments, a computer-readable storage medium is disclosed, in which a plurality of instructions are stored, and the instructions are adapted to be loaded by a processor of a terminal device and execute the SPH simulation display method of the overall process of the catastrophic evolution of the sea-bottom landslide in the first embodiment.
Although the embodiments of the present invention have been described with reference to the accompanying drawings, it is not intended to limit the scope of the present invention, and it should be understood by those skilled in the art that various modifications and variations can be made without inventive efforts by those skilled in the art based on the technical solution of the present invention.
Claims (10)
1. An SPH simulation display method for the whole process of submarine landslide catastrophe evolution is characterized by comprising the following steps:
dispersing the calculation areas of the seabed slope body and the seawater into a set number of particles with different physical state attributes, and giving corresponding physical mechanical parameters to the physical state particles; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary;
applying the infinitely-covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and the structural acting force applied to a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and simultaneously applying a displacement boundary condition;
in each calculation time step, determining a catastrophe evolution stage of the seabed slope body, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step, and completing updating of particle physical mechanics information of solution domain;
and displaying the physical and mechanical information of the particles at the set time step.
2. The SPH simulation display method for the whole process of the catastrophic evolution of the submarine landslide according to claim 1, wherein the calculation area of the submarine slope body is dispersed into particles with set number and different physical-state attributes, and corresponding physical and mechanical parameters are given to the physical-state particles, and the method specifically comprises the following steps:
the method comprises the steps of dividing a seabed slope body into units, creating slope model boundary subunits by combining slope surface topography and slope geometric boundaries, performing internal interpolation by means of slope profiles, uniformly inserting the model subunits, constructing unit groups of the whole slope model, selecting unit center points as coordinates of slope SPH particles, and giving physical and mechanical parameters of quality, density, stress, position, speed and motion period to the SPH particles.
3. The SPH simulation display method for the whole process of the catastrophic evolution of the submarine landslide according to claim 1, wherein before determining the stage of the catastrophic evolution of the submarine slope, the method further comprises: and selecting a particle approximate kernel function and a smooth kernel radius, and performing particle pairing to form a tight branch domain and an influence domain of all SPH particles.
4. The SPH simulation display method of the whole process of the catastrophic evolution of the submarine landslide according to claim 3, wherein in each calculation time step, determining the catastrophic evolution stage of the submarine slope body specifically comprises:
and in each calculation time step, performing traversal search in a particle tight branch region, acquiring the relative speed among the particles of the slope body and the content ratio of the particles of the sea water and the particles of the slope body in the tight branch region, and determining the catastrophe evolution stage of the seabed slope body.
5. The SPH simulation display method for the whole process of the catastrophic evolution of the seabed landslide as claimed in claim 4, wherein the content ratio of seawater particles and slope particles in the particle system in different stages of the catastrophic evolution is calibrated through indoor test tests of the particulate matter concentration and the volume concentration of the seabed landslide body material and the viscosity of the landslide body material;
and obtaining a corresponding catastrophe evolution stage based on the obtained particle content ratio in the calculation time step.
6. The SPH simulation display method for the whole process of the catastrophic evolution of the submarine landslide according to claim 1, wherein the catastrophic evolution stage of the submarine slope body specifically comprises:
the stage of the integral movement and the sliding of the sliding mass, at the moment, the relative movement speed between solid phase particles of the sliding mass is close to zero;
when the water content of the slope body reaches a set value, the slope body material is changed from a solid state to a liquid state, relative displacement occurs among the slope body materials, and the migration of the sliding slope body is in a debris flow mode;
the seawater interacts with the landslide body along with the continuous sliding of the landslide body, and when the volume of solid particles of the landslide body exceeds a set value, the landslide body slides to present the mud flow characteristic of the non-Newtonian fluid;
and as the slope body continues to slide, the water content of the slope body material is further increased, and the sliding of the slope body presents a turbidity current flow mode.
7. The SPH simulation display method for the whole process of the catastrophic evolution of the submarine landslide as claimed in claim 1, wherein the calculation and the solution are performed within a single time step, and the updating of the particle physics mechanics information of the solution domain is completed, and the method specifically comprises the following steps:
in each time step, solving an SPH approximate expression of a continuity equation of the particle system through a kernel approximation method and a particle approximation method to obtain a particle density change rate in a tight branch domain; then, stress solution and momentum change rate calculation are carried out by a system state equation; and updating the particle coordinate, stress, acceleration and momentum change rate information of the solution domain.
8. The utility model provides a SPH simulation display system of seabed landslide catastrophe evolution overall process which characterized in that includes:
the ionization module is used for dispersing the calculation area of the seabed slope body and the seawater into particles with set quantity and different physical state attributes, and endowing the particles with the corresponding physical and mechanical parameters to each physical state; and a non-slip boundary layer combining repulsive particles and dummy particles is arranged outside the boundary;
the acting force applying module is used for applying the infinitely covered seawater gravity and periodic wave action on seawater particles, converting the land slope acting force and structural acting force borne by a slope body into a stress boundary condition of a slope body particle calculation area, applying the stress boundary condition on boundary particles consisting of repulsive particles and virtual particles, and applying a displacement boundary condition;
the particle information updating module is used for determining a catastrophe evolution stage of the seabed slope body in each calculation time step, selecting a system constitutive equation set corresponding to the catastrophe evolution stage, performing calculation solving in a single time step and completing particle physical and mechanical information updating of a solution domain;
and the catastrophe evolution overall process display module is used for displaying the particle physical and mechanical information of the set time step.
9. A terminal device comprising a processor and a computer-readable storage medium, the processor being configured to implement instructions; the computer readable storage medium is used for storing a plurality of instructions, wherein the instructions are suitable for being loaded by a processor and executing the SPH simulation display method of the whole process of the disaster evolution of the seabed landslide as set forth in any one of claims 1-7.
10. A computer readable storage medium having stored therein a plurality of instructions, wherein the instructions are adapted to be loaded by a processor of a terminal device and to execute the SPH simulation demonstration method of the whole process of the catastrophic evolution of the seafloor landslide of any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110863595.3A CN113761775B (en) | 2021-07-29 | 2021-07-29 | SPH simulation display method and system for whole disaster evolution process of submarine landslide |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110863595.3A CN113761775B (en) | 2021-07-29 | 2021-07-29 | SPH simulation display method and system for whole disaster evolution process of submarine landslide |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113761775A true CN113761775A (en) | 2021-12-07 |
| CN113761775B CN113761775B (en) | 2024-03-15 |
Family
ID=78788096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110863595.3A Active CN113761775B (en) | 2021-07-29 | 2021-07-29 | SPH simulation display method and system for whole disaster evolution process of submarine landslide |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113761775B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115618651A (en) * | 2022-11-22 | 2023-01-17 | 交通运输部天津水运工程科学研究所 | Physical simulation method for research on water intake entrainment effect of offshore power plant |
| CN115719012A (en) * | 2023-01-06 | 2023-02-28 | 山东科技大学 | Tailing pond ore drawing arrangement method based on unmanned aerial vehicle remote sensing and multiphase SPH algorithm |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819650A (en) * | 2012-08-16 | 2012-12-12 | 同济大学 | Computational simulation method of flow slide catastrophe of rock and soil material |
| CN104991999A (en) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | Dam bursting flood routing simulation method based on two-dimensional SPH |
| CN109284523A (en) * | 2018-07-19 | 2019-01-29 | 同济大学 | A kind of rock soil medium Progressive failure, class solid-liquid phase change behavior analogy method |
| CN112016224A (en) * | 2020-07-28 | 2020-12-01 | 西南大学 | SPH-based soil slope landslide analysis and determination method, system, terminal and medium |
-
2021
- 2021-07-29 CN CN202110863595.3A patent/CN113761775B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819650A (en) * | 2012-08-16 | 2012-12-12 | 同济大学 | Computational simulation method of flow slide catastrophe of rock and soil material |
| CN104991999A (en) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | Dam bursting flood routing simulation method based on two-dimensional SPH |
| CN109284523A (en) * | 2018-07-19 | 2019-01-29 | 同济大学 | A kind of rock soil medium Progressive failure, class solid-liquid phase change behavior analogy method |
| CN112016224A (en) * | 2020-07-28 | 2020-12-01 | 西南大学 | SPH-based soil slope landslide analysis and determination method, system, terminal and medium |
Non-Patent Citations (4)
| Title |
|---|
| G FOURTAKAS: "Modelling multi-phase liquid sediment scour and resuspension induced by raqid flows using smooth particle hydrodynamics accelerated with a Graphics Processing Unit", ADVANCES IN WATER RESOURCES * |
| 余和雨: "海底碎屑流运动特性研究的若干进展", 地质科技情报 * |
| 胡嫚: "基于弹塑性土体本构模型的滑坡运动过程SPH模拟", 岩土工程学报, vol. 38, no. 1 * |
| 颜东谊: "海底滑坡诱发海洋内波实验研究", 硕士电子期刊 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115618651A (en) * | 2022-11-22 | 2023-01-17 | 交通运输部天津水运工程科学研究所 | Physical simulation method for research on water intake entrainment effect of offshore power plant |
| CN115719012A (en) * | 2023-01-06 | 2023-02-28 | 山东科技大学 | Tailing pond ore drawing arrangement method based on unmanned aerial vehicle remote sensing and multiphase SPH algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113761775B (en) | 2024-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fuhrman et al. | Numerical simulation of wave-induced scour and backfilling processes beneath submarine pipelines | |
| Ren et al. | SPH-DEM modeling of the hydraulic stability of 2D blocks on a slope | |
| Paulsen et al. | An efficient domain decomposition strategy for wave loads on surface piercing circular cylinders | |
| Zakeri | Submarine debris flow impact on suspended (free-span) pipelines: Normal and longitudinal drag forces | |
| Dixen et al. | Numerical and experimental investigation of flow and scour around a half-buried sphere | |
| CN113761775B (en) | SPH simulation display method and system for whole disaster evolution process of submarine landslide | |
| Wu et al. | An implicit 2-D shallow water flow model on unstructured quadtree rectangular mesh | |
| Kong et al. | A high-resolution method for the depth-integrated solute transport equation based on an unstructured mesh | |
| Liu et al. | Two-dimensional scour simulations based on coupled model of shallow water equations and sediment transport on unstructured meshes | |
| Duan et al. | PORO-FSSI-FOAM: Seabed response around a mono-pile under natural loadings | |
| Finnegan et al. | Linear irregular wave generation in a numerical wave tank | |
| Zhang et al. | Coupling irregular particles and fluid: Complex dynamics of granular flows | |
| Leroyer et al. | Numerical strategies to speed up CFD computations with free surface—Application to the dynamic equilibrium of hulls | |
| CN113761774A (en) | Method and system for simulating landslide process of seabed | |
| Heß et al. | Debris flows with pore pressure and intergranular friction on rugged topography | |
| Bakhtyar et al. | Process-based model for nearshore hydrodynamics, sediment transport and morphological evolution in the surf and swash zones | |
| Yin et al. | Deep learning-based design model for suction caissons on clay | |
| Hejazi et al. | Numerical modeling of breaking solitary wave run up in surf zone using incompressible smoothed particle hydrodinamics (ISPH) | |
| Rattia et al. | Modelling local scour near structures with combined mesh movement and mesh optimisation | |
| Liu et al. | Numerical modeling of submarine turbidity currents over erodible beds using unstructured grids | |
| Toosi et al. | The influence of time scale in free surface flow simulation using Smoothed Particle Hydrodynamics (SPH) | |
| Ren et al. | A coupled metaball discrete element material point method for fluid–particle interactions with free surface flows and irregular shape particles | |
| Gushchin et al. | On one approach for mathematical modeling of the mixed zones dynamics in a stratified fluid | |
| Ginting et al. | Tapping the potential of shallow water model for wave simulations | |
| Tabasi et al. | A Copula-Based Simulation of Wave-Induced Pore Water Pressure Gradient and Local Acceleration Within Surf Zone for Natural and Laboratory Barred Beach Profiles |
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 |