CN104021277B - Numerical analysis method for piping phenomenon - Google Patents
Numerical analysis method for piping phenomenon Download PDFInfo
- Publication number
- CN104021277B CN104021277B CN201410202868.XA CN201410202868A CN104021277B CN 104021277 B CN104021277 B CN 104021277B CN 201410202868 A CN201410202868 A CN 201410202868A CN 104021277 B CN104021277 B CN 104021277B
- Authority
- CN
- China
- Prior art keywords
- model
- alpha
- similarity criterion
- prototype
- piping
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a numerical analysis method for the piping phenomenon. The method comprises the steps of obtaining a seepage deformation similarity criterion on the basis of a fluid dynamics similarity criterion combined with alternate force acting item of medium on fluid; building a microscopic particle model on the basis of the seepage deformation similarity criterion combined with a particle flow method, wherein the model is consistent with a prototype in gravity horizontality; generating skeleton and filling particles under the conditions of zooming particle size distribution and maintaining the same porosity; conducting piping process simulation on the model, and outputting parameter information inside the model when the filling particles in the model reach the balance state or the skeleton particles are damaged; obtaining the actual numerical value by combining the seepage deformation similarity criterion with the parameter information inside the model, and feeding the actual numerical value back to the prototype so as to analyze the piping phenomenon. The numerical analysis method for the piping phenomenon conforms to the Darcy law and the non-Darcy law, the model can truly describe the prototype when two flow states exist at the same time, the problem that a conventional model or a centrifugal test cannot achieve same-ratio zooming of the particle diameter and the model is solved, and analysis is accurate.
Description
Technical field
The present invention relates to a kind of numerical analysis method of Piping phenomenon, belongs to technical field of civil engineering.
Background technology
China is to be subjected to one of flood damage country the most serious in the world, middle and lower reach of Yangtze River region of no relief big flood particularly frequency
It is numerous, serious.The destruction of dyke includes the various forms such as piping, Man Ding, collapse on the bank and overall collapse, and a large amount of big flood data show, dike
Base tube gushes harm maximum, and when dykes and dams occur failure by piping, seepage field has strong spatial character, at present to piping evolution
Seepage field spatial profiles in terms of research it is less, return its reason, on the one hand, because the current understanding to piping genesis mechanism is also
It is not deep enough；On the other hand, due to the complexity of engineering hydraulics, the analysis to seepage field is more difficult.
Currently, engineering circles and theoretical circles are to biasing toward the anti-permeability performance aspect research to soil the research of piping problem more, so
And with the development of piping research work, there is strong space-time characterisation gradually to have reached a kind of common recognition to seepage field.If retouching
The complicated change of material geometrical property and hydraulic characteristic(s) in flow event is stated, head, hydraulic gradient, particle displacement field and hole is obtained
The dynamic change of gap rate equivalent, then will provide strong support for the research of piping process.
Before the present invention, the research to piping is divided into theoretical research, laboratory test and numerical experimentation research, and permeates and become
Shape test in commonly assume that flow event meets Darcy's law, when occurring non-darcy flow in domain, using conventional model or from
There is certain deviation with prototype in the conclusion that heart model is obtained.In centrifugal test, such as tried using the grain diameter for reducing
Testing can face a problem：Prototype particle diameter range in the range of sand, and after particle diameter and model length reduce on year-on-year basis than chi,
Particle diameter may be already at clay scope, and when there is percolation phenomenon in clay and sand, the active force between soil and water faces
Larger difference, surface of clay there may be the strong and weak effect with reference to water, and in sand be not in then such case, namely adopt
Although similar with realizing in the method mechanism for reducing particle diameter on year-on-year basis, the physics of thin sight category, the impact of chemical action are very
Result and expected difference may be caused larger；Selected using reduced scale relation simultaneously, ensureing that it is similar same that model and prototype level are matched somebody with somebody
When can not eliminate prototype soil and model between differences of mechanical properties；The face particularly when from the fluid for meeting specified conditions
Face larger difficulty.Numerical analysis for piping adopts finite element modelling seepage field, continuity method that piping to occur mostly
The simulation in front relative steady seepage stage is feasible, after critical condition is reached, it is impossible to consider that particle is lost in the soil body geometry for causing
Characteristic and the interactional this complicated Soil water interactions of hydraulic characteristic(s), therefore cannot comprehensively explain piping genesis mechanism.Piping
Complicated Soil water interactions determine that the development of piping is a nonlinear dynamic process during generation, and public affairs are there is no at present
The most suitable research method recognized and theory.Existing scholar begins attempt to combine numeral using microimaging visualization tracking technique
Image recognition analysis means, to the whole process of sand piping record is tracked, and from thin aspect of seeing Soil water interactions are disclosed
Through the overall process of piping development, but the technology is limited only to laboratory test aspect, for actual piping seepage failure mistake
The analysis of journey does not still possess adaptability.
The content of the invention
The technical problem to be solved is to overcome the deficiencies in the prior art, there is provided one kind is accurate easily, energy
The numerical analysis method of the true a kind of Piping phenomenon for reflecting actual piping seepage failure process, based on existing hydrodynamics phase
It is dynamic based on complete porous media fluid like criterion, it is considered to the alternate power effect under darcy-non-darcy state between medium and fluid
Mechanical equation derives seepage failure similarity criterion, provides the similarity criterion that sand occurs to be followed during seepage deformation so that work as calculating
Occur darcy in region, in non-darcy flow, or domain while model can be carried out truly to prototype when there is two kinds of fluidised forms
Description.
The present invention specifically employs the following technical solutions solution above-mentioned technical problem：
A kind of numerical analysis method of Piping phenomenon, comprises the following steps：
Step (1), the alternate power effect item based on hydrodynamics similarity criterion binding medium to fluid, are expanded by introducing
The D-B-F equations of exhibition consider the darcy-non-Darcy effects of seepage flow, obtain the similarity criterion that sand occurs to be followed during seepage deformation；
Step (2), particle stream method is combined based on seepage deformation similarity criterion set up the thin sight granule die similar to prototype
Type；It is consistent with prototype that the thin sight granular model meets gravity horizontal, and according to identical with prototype in carefully sight granular model
Grading curve prototype soil sample particle is scaled by seepage deformation similarity criterion and keeps random under the conditions of same porosity
Generate skeleton and filler particles；
Step (3), to it is described it is thin sight granular model carry out piping process simulation, filler particles reach equilibrium-like in model
By model intrinsic parameter information output when state or skeleton particle are destroyed；
Step (4), the seepage deformation similarity criterion is combined acquisition actual numerical value with model intrinsic parameter information, and will be real
Border numeric feedback to prototype, to analyze Piping phenomenon.
Further, as a preferred technical solution of the present invention：Step (1) sand occurs to be abided by during seepage deformation
The similarity criterion followed is：
α in formulan、ασ、αγ、αt、αv、αl、αg、αρf、αμ、αk、αdRepresent respectively：Porosity, stress, severe, the time,
Speed, length, acceleration, fluid density, kinematic coefficient of viscosity, infiltration coefficient affinity constant, equivalent coefficient of permeability affinity constant,
Particle diameter zoom factor；A () formula characterization model must be consistent with prototype mesopore rate；(b) formula characterization model and prototype stress state
Unanimously, it is consistent comprising effective stress and total stress correspondence；(c)～(f) formulas characterize respectively hydrodynamics it is similar in it is humorous when it is accurate
Then, gravity similarity criterion, pressure similarity criterion, viscous force similarity criterion；G (), (h) characterize pore-fluid in laminar flow and turbulent flow
The alternate power exchange coefficient impact similar to fluid dynamic in porous media during state.
Further, as a preferred technical solution of the present invention：Thin granular model of seeing is by infiltration in the step (2)
Deformation similarity criterion scaling meets：Grain diameter and model scale on year-on-year basis with fluid density and moulded dimension scale on year-on-year basis, it is viscous
Coefficient ratio chi for model scale square.
Further, as a preferred technical solution of the present invention：Model intrinsic parameter delta data in the step (3)
Including model endoparticle speed, stress distribution, flow, flow velocity, pressure data.
Further, as a preferred technical solution of the present invention：The step (4) is by the similar standard of the seepage deformation
Then combine with model intrinsic parameter delta data and obtain the critical hydraulic gradient consistent with prototype, and reference passes through seepage deformation
The time scale ratio analyzer tube consistent with inertia time that similarity criterion is obtained comes through journey.
The present invention adopts above-mentioned technical proposal, can produce following technique effect：
The advantage and effect of the numerical analysis method of Piping phenomenon of the present invention be：
(1) Darcy's law is all based on when using conventional model and centrifugal model research Seepage problems, is occurred in a model
The related conclusions based on Darcy's law are still adopted to be inappropriate during non-darcy flow.Instant invention overcomes being tried using conventional model
Test or the limitation faced during sand seepage deformation problem is studied in centrifugal model test.
(2) based on existing porous media hydrodynamics similarity criterion, it is considered to which medium acts on item to the alternate power of fluid, gives
Go out the similarity criterion that sand occurs to be followed during seepage deformation, its core concept is to ensure that grain diameter is scaled on year-on-year basis with model,
And ensureing that the coefficient of viscosity is more consistent than chi with length than chi, the flow model in porous media set up as premise is similar to prototype.When calculating area
Occur darcy in domain, in non-darcy flow, or domain while model can truly be retouched to prototype when there is two kinds of fluidised forms
State, thus when in prototype fluidised form when darcy is converted to non-darcy, based on seepage deformation similarity criterion set up model for this
Process also can be described exactly.
(3) with particle as basic computational ele- ment, the scaling for simulating particle has original advantage to particle stream method, with reference to it
The advantage of research seepage deformation, overcoming conventional model test or centrifugal test can not realize that grain diameter is scaled on year-on-year basis with model
A difficult problem.
(4) research is carried out in analysis platform, and constraints is few, applied range.Avoid indoor style prepare,
The complicated processes such as instrument and equipment operation, improve test efficiency, have saved research cost.
Operation of the present invention is simple and clear, on the basis of particle flow theory, based on porous media hydrodynamics similarity criterion, examines
Consider effect item of the medium to fluid, give the similarity criterion that sand occurs to be followed during seepage deformation, set up model manipulation it is convenient,
Calculate simple, economical and efficient, analyze more accurate.
Description of the drawings
Fig. 1 be Piping phenomenon of the present invention numerical analysis method in set up and thin see granular model schematic diagram.
Fig. 2 sees granular model boundary condition schematic diagram for the present invention is thin.
Fig. 3 is the flow chart of the numerical analysis method of Piping phenomenon of the present invention.
Fig. 4 A are the comparing calculation of critical hydraulic gradient when carefully seeing aperture A generations piping in granular model and laboratory test
Result figure.
Fig. 4 B are the comparing calculation of critical hydraulic gradient when carefully seeing aperture B generations piping in granular model and laboratory test
Result figure.
Fig. 4 C are the comparing calculation of critical hydraulic gradient when carefully seeing aperture C generations piping in granular model and laboratory test
Result figure.
Fig. 5 A are the schematic diagram that sample I-1 of the present invention occurs seepage failure initial time when C apertures are opened.
Fig. 5 B are that sample I-1 of the present invention occurs the schematic diagram that seepage failure fine grained is lost in the moment when C apertures are opened.
Fig. 5 C are that sample I-1 of the present invention occurs the schematic diagram that seepage failure coarse granule is lost in the moment when C apertures are opened.
Fig. 5 D are the schematic diagram that sample I-1 of the present invention models when C apertures are opened destroy the moment.
Specific embodiment
Embodiments of the present invention are described with reference to Figure of description.
The present invention devises a kind of numerical analysis method of Piping phenomenon, the concrete steps of enforcement as shown in figure 3, including：
Step (1), the effect item based on hydrodynamics similarity criterion binding medium to fluid, derive and obtain sand generation
The similarity criterion followed during seepage deformation.It is specific as follows：
(1.1) based on existing porous media hydrodynamics similarity criterion, it is considered to alternate dynamic action of the medium to fluid
, introduce the D-B-F equations of extension；
(1.2) consider the darcy-non-Darcy effects of seepage flow, according to geometric similarity, kinematic similitude, dynamic similarity, derive sand
There is the similarity criterion followed during seepage deformation in soil.Concrete derivation is as follows：
In hydrodynamics, it is assumed that do not consider the change in temperature field, for incompressible fluid, solve by bulk solid group
Into porous media in flow field and have to comply with mass-conservation equation during pressure field：
In formula (1), n is the porosity of porous media,For fluid velocity,
Mean momentum conservation equation is：
Be different from N-S equations, the present invention considers effect item of the medium to fluid, refer to equation right-hand member last.
The Darcy-Brinkman-Forchheimer equations that formula (2) is referred to as extending by the present invention, left end Section 1 is represented
Local inertial force, Section 2 represents displacement inertia force, and two sums characterize inertia force, namely acceleration item；Equation right-hand member is every
The contribution of body force, pressure, viscous force and solid-liquid inter-phase forces to acceleration is characterized respectively.
N is the porosity of porous media in formula (2),For fluid velocity,For particle average movement velocity, τ is viscous
Power tensor, g is acceleration of gravity, and ρ f are fluid density, and p is pressure, and β is alternate power exchange coefficient, according to Ergun equations and
Wen and Yu (1966) proposes that formula determines, embodies as follows:
μ in formula (3)fTo characterize fluid viscous coefficient,For particle equivalent grain size.
Corresponding barometric gradient and apparent velocity can be expressed as：
Dp/dx characterizes barometric gradient in formula (4).
Equation Section 1 is Kozeny-Carman formula during laminar condition, and Section 2 represents barometric gradient and speed
Quadratic term it is proportional.When two-phase does not account for leading, between 1～2, this permeates and is the relation of barometric gradient and speed
Number ratios be gravity than chi x powers the reason for.As only considered velocity squared item when being completely in turbulent condition, ignore first
, formula (5) can be expressed as：
For convenience of illustrating, by formula (5)Equivalent coefficient of permeability is defined as, i is hydraulic gradient, and g is acceleration of gravity, is tied
Darcy's law expression formula v=ki during laminar condition is closed, k is infiltration coefficient.Make β1、β2The two of difference representative formula (3) right-hand member
, β1、β2Alternate force coefficient under respectively complete laminar flow and complete turbulent-flow conditions, then formula (2) can be expressed as：
Definition time, speed, length, acceleration, fluid density, kinematic coefficient of viscosity, alternate power exchange affinity constant
Respectively：αt、αv、αl、αg、αρf、αμ、αβ, affinity constant is substituted into into momentum conservation equation (2) and is obtained：
In formula (7), αdFor particle diameter zoom factor, the definition of τ is lateral fluid force.αnFor porosity affinity constant.
The similar basic demand of hydrodynamics is followed, the power ratio of the same name acted in respective point is calculated with unit volume should
When equal, convolution (6) has：
In formula (8)Alternate force coefficient affinity constant under respectively complete laminar flow and complete turbulent-flow conditions.
Formula (6) is substituted into into formula (8), and it is similar with reference to gravity, and porosity is similar, finally gives：
In formula in (9), αkFor infiltration coefficient affinity constant,For equivalent coefficient of permeability affinity constant.
(a) formula characterization model must be consistent with prototype mesopore rate in formula (9)；(b) formula characterization model and prototype stress shape
State is consistent, consistent comprising effective stress and total stress correspondence；(c)～(f) formulas characterize respectively hydrodynamics it is similar in it is humorous when
Criterion, gravity similarity criterion, pressure similarity criterion, viscous force similarity criterion；G (), (h) characterize pore-fluid in laminar flow and disorderly
The alternate power exchange coefficient impact similar to fluid dynamic in porous media during stream mode.
Step (2), particle stream method is combined based on seepage deformation similarity criterion set up the thin sight granule die similar to prototype
Type.It is specific as follows：
(2.1) indoor thin sight granular model is set up, as shown in figure 1, upstream boundary overlaps with the plane of x=0, in upstream edge
Boundary applies head, and other surfaces are set to impervious boundary, and A, B, C aperture (starting equal body of wall to cover) is reserved in upper surface, greatly
Little to follow model scale, to simulate the different head pieces that go out, orifice center position is set to zero-pressure force boundary, and simulation different orifice goes out stream
When delete the wall element of relevant position respectively, the model boundary condition schematic diagram is as shown in Fig. 2 model left side is head side
Boundary, aperture is zero-pressure force boundary, and other borders are set to the waterproof non-slip border of rigidity.
Design parameter is shown in Table 1；
(2.2) according to the seepage deformation similarity criterion of step 1 gained, with reference to grain composition in prototype, particle is carefully seen
Particle size range determines in model；
(2.3) and then, in a model according to prototype identical grading curve by prototype soil sample particle press similar standard
Then scale, scaling includes that grain diameter and model are scaled on year-on-year basis, fluid density and moulded dimension are scaled on year-on-year basis, and the coefficient of viscosity compares chi
For model scale square.To avoid contact with the generation of wash away phenomenon, model coboundary coefficient of friction is set to into higher value, it is raw
Coarse granule is first produced using PFC3D particles random generator into after boundary model, fine grained is produced afterwards, retaining hole in generating process
Gap rate is consistent with prototype, to keep same porosity under the conditions of random generate skeleton and filler particles.
(2.4) by carefully see granular model can perfection meet seepage deformation similarity criterion, model must be fulfilled for gravity horizontal with
Prototype is consistent, and by improving gravity horizontal model scale impact is offset.Guarantee that carefully seeing granular model reaches primary stress balance
State, simulates without original state under the conditions of seepage flow, and model reaches and make first after original steady-state model in the hydrostatic pressure.
Step (3), to it is thin see granular model and apply head boundary carry out seepage deformation, simulate piping process, record cast
Intrinsic parameter delta data.Specially：
(3.1) seepage deformation analysis is carried out to thin granular model applying head boundary of seeing, is separately turned on different apertures and (deletes
Except the overlying body of wall of relevant position), and the scope of freedom is set to, there is seepage failure in simulation different orifice, by monitoring hand
Section obtains particle volume, orifice flow and the unit internal porosity that relevant parameter such as aperture is flowed out in model, particle position, particle
Migration velocity, pressure, flow velocity etc., analysis model endoparticle speed, stress distribution and Penetration Signature, analysis fluid velocity, pressure
And the bonding force between particle and fluid, to obtain relevant parameter information in model in real time in simulation process.
(3.2) under every one-level hydraulic gradient, should when aperture flow velocity reaches stable state or can not keep stable state
Dialysis under level head terminates.As destroyed, apply next stage head, enable the model in hydrostatic pressure, again
Start analysis, to eliminate the impact that the lower partial particulate of upper level head effect has occurred and that migration, until simulated object can not be protected
Keep steady till fixed or orifice outflow stablize, otherwise continue slow applying next stage head.Up to the filler particles in model
Reach poised state or model skeleton particle is destroyed, by model intrinsic parameter data message output.
Step (4), by the seepage deformation similarity criterion combined with model intrinsic parameter delta data acquisition actual numerical value, root
Factually border numeric feedback is to prototype, to analyze Piping phenomenon.Specially：
The seepage deformation similarity criterion being derived by with step (1), to relevant parameter number in step (3) gained model
According to such as particle speed, stress distribution, flow, flow velocity data are analyzed, and based on the similarity criterion for deriving, obtain actual number
Value, according to actual numerical value Piping phenomenon is analyzed.Obtain actual numerical value and obtain correspondence parameter values as needed during analysis,
The preferred acquisition critical hydraulic gradient consistent with prototype in the present invention, and with reference to being obtained by seepage deformation similarity criterion and used
The time scale ratio of property time consistency, according to the change of piping evolution under same time scale ratio, is carried out to actual piping process
Characterize.
The seepage deformation similarity criterion is combined with model intrinsic parameter delta data, specifically with aperture discharge velocity in model
As a example by, according to gravity similarity criterion and viscous force similarity criterion in formula (9), it is consistent with flow velocity in prototype to be obtained model, such as
This, flow velocity and actual flow velocity is represented in model；Now as continued flow and flow velocity relation, according to：Q=VAt, binding model and original
The consistent conclusion of type flow velocity, is 1/N due to flowing through sectional area than chi2, time scale ratio is 1/N, can obtain model and prototype flow
It is 1/N than chi3.Thus the actual numerical value of prototype flow can be obtained, then is analyzed.
Thus, this method is based on porous media hydrodynamics similarity criterion, it is considered to which medium is acted on the alternate power of fluid
, the similarity criterion that sand occurs to be followed during seepage deformation is given, its core concept is to ensure that grain diameter with model on year-on-year basis
Scaling, the flow model in porous media set up as premise and prototype basic simlarity.When occurring darcy, non-darcy flow in zoning,
Or model can truly be described to prototype when in domain there is two kinds of fluidised forms simultaneously, thus when in prototype fluidised form by darcy
When converting to non-darcy, the model based on similarity criterion foundation also can be described exactly for the process, the number being consequently formed
Value analysis method can solve the problem that during existing piping analysis and cannot overcome while meet the technical problem of west, non-Darcy effects
Conventional model test or centrifugal test can not realize the difficult problem that grain diameter is scaled on year-on-year basis with model, can accurately apply to reality
Border piping analysis process.
In order to the numerical analysis method for verifying the Piping phenomenon of the present invention can simulate Piping phenomenon, to the infiltration in model
Deformation is analyzed, and spy is verified with following experimental data.
First, according to parameter shown in table 1, respectively 6 groups of soil samples are obtained by laboratory test and occurs to permeate broken in different orifice
The critical hydraulic gradient of bad when.Then according to the method for the present invention, sets up thin sight granular model analysis and obtains correspondence hydraulic gradient.
Apparent critical hydraulic gradient refers to Fig. 4 A, Fig. 4 B, Fig. 4 C with the variation relation of aperture A, B, C position.Can be with
Find out, both coincide substantially, meet critical hydraulic gradient with fine material content and fill the increase of packing and becoming for increasing
Gesture is actual with engineering to match.Fig. 5 A are the schematic diagram that sample I-1 occurs seepage failure initial time when C apertures are opened,
Starting stage, model endoparticle initial balance；Subsequently fine grained only a small amount of in sample is lost in, and is No. C as shown in Figure 5 B
There is the schematic diagram that fine grained in seepage failure is lost in the moment when aperture is opened；As fine grained is assembled to orifice area, thin
After grain is lost in the regular hour, individual coarse particles start to be taken by fine grained to wrap up in takes aperture out of, elapses over time, and coarse granule starts
Gradually it is lost in, is that the schematic diagram that seepage failure coarse granule is lost in the moment occurs when C apertures are opened as shown in Figure 5 C；Final apertures
There is larger pit in mouth region domain, is schematic diagram that model destroys the moment as shown in Figure 5 D, and it is actual with engineering to match.Together
When, whole numerical simulation is all performed by computer, it is to avoid the loaded down with trivial details mistake of sample, operating instrument is prepared in laboratory test
Journey, improves test efficiency, has saved research cost.
The model parameter of table 1
Thus experimental verification can enter Mobile state observation using the simulation failure by piping of particle stream method to each parameter in model,
The Meso process of failure by piping can be reflected.Prove that the present invention accurately, easily can be explained and analyzed to Piping phenomenon.Should
Understand that these embodiments are only illustrative of the invention and is not intended to limit the scope of the invention, after the present invention has been read, this
Modification of the art personnel to the various equivalent form of values of the present invention falls within the application claims limited range.
Claims (4)
1. a kind of numerical analysis method of Piping phenomenon, it is characterised in that comprise the following steps：
Step (1), the alternate power effect item based on hydrodynamics similarity criterion binding medium to fluid, by introducing extension
D-B-F equations consider the darcy-non-Darcy effects of seepage flow, obtain the similarity criterion that sand occurs to be followed during seepage deformation；
Step (2), particle stream method is combined based on seepage deformation similarity criterion set up the thin sight granular model similar to prototype；Institute
State it is thin see granular model to meet gravity horizontal consistent with prototype, and in carefully sight granular model according to prototype identical particle
Prototype soil sample particle is scaled and keeps random under the conditions of same porosity and generates bone by grading curve by seepage deformation similarity criterion
Frame and filler particles；
Step (3), to it is described it is thin sight granular model carry out piping process simulation, in model filler particles reach poised state or
By model intrinsic parameter information output when skeleton particle is destroyed；
Step (4), the seepage deformation similarity criterion is combined acquisition actual numerical value with model intrinsic parameter information, and by actual number
Value feeds back to prototype, to analyze Piping phenomenon；
Wherein, the similarity criterion that step (1) sand occurs to be followed during seepage deformation is：
α in formulan、ασ、αγ、αt、αv、αl、αg、αρf、αμ、αk、αdRepresent respectively：Porosity, stress, severe, the time, speed,
Length, acceleration, fluid density, kinematic coefficient of viscosity, infiltration coefficient affinity constant, equivalent coefficient of permeability affinity constant, particle diameter
Zoom factor；A () formula characterization model must be consistent with prototype mesopore rate；B () formula characterization model is consistent with prototype stress state,
It is consistent comprising effective stress and total stress correspondence；(c)～(f) formulas characterize respectively hydrodynamics it is similar in it is humorous when criterion, weight
Power similarity criterion, pressure similarity criterion, viscous force similarity criterion；G (), (h) characterize pore-fluid in laminar flow and turbulent condition
When the impact similar to fluid dynamic in porous media of alternate power exchange coefficient.
2. the numerical analysis method of Piping phenomenon according to claim 1, it is characterised in that：Thin sight in the step (2)
Grain model is met by seepage deformation similarity criterion scaling：Grain diameter scales on year-on-year basis same with moulded dimension with fluid density with model
Than scaling, the coefficient of viscosity than chi for model scale square.
3. the numerical analysis method of Piping phenomenon according to claim 2, it is characterised in that：In the step (3) in model
Parameters variation data include model endoparticle speed, stress distribution, flow, flow velocity, pressure data.
4. the numerical analysis method of Piping phenomenon according to claim 1, it is characterised in that：The step (4) is oozed described
Saturating deformation similarity criterion is combined with model intrinsic parameter delta data and obtains the critical hydraulic gradient consistent with prototype, and reference passes through
The time scale ratio analyzer tube consistent with inertia time that seepage deformation similarity criterion is obtained comes through journey.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410202868.XA CN104021277B (en) | 2014-05-14 | 2014-05-14 | Numerical analysis method for piping phenomenon |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410202868.XA CN104021277B (en) | 2014-05-14 | 2014-05-14 | Numerical analysis method for piping phenomenon |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104021277A CN104021277A (en) | 2014-09-03 |
| CN104021277B true CN104021277B (en) | 2017-04-19 |
Family
ID=51438029
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410202868.XA Expired - Fee Related CN104021277B (en) | 2014-05-14 | 2014-05-14 | Numerical analysis method for piping phenomenon |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104021277B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108088870A (en) * | 2017-12-08 | 2018-05-29 | 天津城建大学 | The method for building up of Nonlinear Heat Conduction model test similarity criterion |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104931401B (en) * | 2015-06-02 | 2017-10-20 | 中国科学院力学研究所 | Infiltration coefficient dynamic change model in a kind of sandy gravel soils piping erosion process |
| CN105225583A (en) * | 2015-11-17 | 2016-01-06 | 中国地质大学（武汉） | A kind of Seepage Experiment method in pore media and mould |
| TWI566106B (en) * | 2015-12-08 | 2017-01-11 | 黃良雄 | Simulation method of piping prevention |
| CN107808228B (en) * | 2017-09-27 | 2021-08-20 | 上海勘察设计研究院(集团)有限公司 | Groundwater extraction effectiveness evaluation method for low-permeability soil layer |
| CN108375535B (en) * | 2018-02-26 | 2020-06-23 | 河北建筑工程学院 | Method for establishing fluid-solid coupling test based on post-peak rock stress and permeability characteristics |
| CN108982324A (en) * | 2018-07-20 | 2018-12-11 | 河海大学 | It is a kind of to carry out under different compactness piping evaluation test device and test method in the soil body |
| CN109409367B (en) * | 2018-11-02 | 2021-09-21 | 四川大学 | Infrared image grading identification method based on rock temperature rise characteristics |
| CN111553098B (en) * | 2020-03-26 | 2023-06-02 | 山东大学 | Numerical simulation method, server and medium for seepage damage type gushing water process |
| CN111753414B (en) * | 2020-06-16 | 2023-03-24 | 河海大学 | Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law |
| CN113075106B (en) * | 2021-03-26 | 2022-11-11 | 桂林理工大学 | Method for determining hydraulic parameters and water inflow amount of gravel soil in erosion stage |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1320744A (en) * | 2001-05-14 | 2001-11-07 | 刘彦芝 | Rapid blocking method for hole or piping of flood control dike |
| CN201148608Y (en) * | 2008-01-21 | 2008-11-12 | 梁欣悦 | Piping-proof dam |
-
2014
- 2014-05-14 CN CN201410202868.XA patent/CN104021277B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1320744A (en) * | 2001-05-14 | 2001-11-07 | 刘彦芝 | Rapid blocking method for hole or piping of flood control dike |
| CN201148608Y (en) * | 2008-01-21 | 2008-11-12 | 梁欣悦 | Piping-proof dam |
Non-Patent Citations (5)
| Title |
|---|
| Forchheimer 型非达西渗流参数特征分析;张文娟，等.;《水电能源科学》;20140131;第32卷(第1期);第52-54页 * |
| Numerical simulation of non-Darcian flow through a porous medium;Duoxing Yang，et al.;《Particuology》;20091231(第3期);第193-198页 * |
| Numerical simulation of non-Darcian flows through spaces partially filled with a porous medium;V.A.F. Costa，et al.;《Computers and Structures》;20040520;第1535-1541页 * |
| Research on the Model Test System of Shaft Lining Interacting with Enclosure Rock under Deep and High Water Pressure;Zhao Xiao-min，et al.;《ICECE，2011 international conference on》;20111231;第3139-3143页 * |
| 渗透变形中相似性问题的颗粒流模拟;倪小东，等.;《土木建筑与环境工程》;20090630;第31卷(第3期);第55-60页 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108088870A (en) * | 2017-12-08 | 2018-05-29 | 天津城建大学 | The method for building up of Nonlinear Heat Conduction model test similarity criterion |
| CN108088870B (en) * | 2017-12-08 | 2020-06-19 | 天津城建大学 | Method for establishing non-linear heat conduction model test similarity criterion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104021277A (en) | 2014-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104021277B (en) | Numerical analysis method for piping phenomenon | |
| Ma et al. | Numerical simulation of tsunami waves generated by deformable submarine landslides | |
| González-Cao et al. | On the accuracy of DualSPHysics to assess violent collisions with coastal structures | |
| Greco et al. | A two-phase model for fast geomorphic shallow flows | |
| Hargreaves et al. | Validation of the volume of fluid method for free surface calculation: the broad-crested weir | |
| Zyserman et al. | Modelling morphological processes in the vicinity of shore-parallel breakwaters | |
| Asadollahi et al. | Numerical investigation of tsunami bore effects on structures, part I: Drag coefficients | |
| Guler et al. | Numerical assessment of tsunami attack on a rubble mound breakwater using OpenFOAM® | |
| Han et al. | Estimation of maximum local scour depths at multiple piles of sea/bay-crossing bridges | |
| Guanche et al. | Hybrid modeling of pore pressure damping in rubble mound breakwaters | |
| Meng et al. | Eulerian–Eulerian multiphase models for simulating collapse of submarine sediment column with rheological characteristics in air–water flow | |
| CN112485106B (en) | Physical model layered preparation and test method for controlling soil body state parameters | |
| Abareshi et al. | Equivalent pipe network model for a coarse porous media and its application to 2d analysis of flow through rockfill structures | |
| Karim et al. | Simulation of wave transformation in vertical permeable structure | |
| Cao et al. | Towards modelling wave-induced forces on an armour layer unit of rubble mound coastal revetments | |
| Duan et al. | Numerical investigation of the wave/current–induced responses of transient soil around a square mono-pile foundation | |
| Abdelrazek et al. | Simulation of the erosion and seepage failure around sheet pile using two-phase WC-SPH method | |
| Shigihara et al. | Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast | |
| Li et al. | Numerical simulation of random wave overtopping of rubble mound breakwater with armor units | |
| Issa et al. | Modelling a plunging breaking solitary wave with eddy-viscosity turbulent SPH models | |
| Yin et al. | 3d morphological impact modelling of offshore wind farms using les and hpc | |
| Goursand-Parente et al. | Numerical validation of an energetic dam-break induced flow on a bridge cross section using the SPH method | |
| Khoolosi et al. | Numerical simulation of impulsive water waves generated by subaerial and submerged landslides incidents in dam reservoirs | |
| Tao et al. | Numerical modeling of the mechanisms of piping erosion with coupled CFD and DEM method | |
| Tavouktsoglou et al. | CFD simulation of clear water scour at complex foundations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170419 Termination date: 20200514 |