CN106446439B - The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF - Google Patents
The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF Download PDFInfo
- Publication number
- CN106446439B CN106446439B CN201610886410.XA CN201610886410A CN106446439B CN 106446439 B CN106446439 B CN 106446439B CN 201610886410 A CN201610886410 A CN 201610886410A CN 106446439 B CN106446439 B CN 106446439B
- Authority
- CN
- China
- Prior art keywords
- model
- grout
- dimensional
- vof
- thixotropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
The invention discloses a kind of grout curtain method for numerical simulation that diffusion sharp side is solved based on cement grout thixotropy and VOF: establishing three-dimensional geological model and grouting pore model, couples three-dimensional geological model and grout hole model construction three dimensional unification model;Three dimensional unification model is imported into CFD software for calculation, divides and calculates grid;Cement grout thixotropic fluid model is established, VOF method captures slurries and spreads sharp side;It realizes and considers that slurries thixotropy and VOF calculate the dam foundation grout curtain Three-dimensional simulation on slurries diffusion sharp side under complex geological condition.The present invention realizes the three-dimensional numerical simulation of grouting for calculating under complex geological condition based on thixotropy model and VOF and determining slurries diffusion sharp side, it accurately learns under complex geological condition macro-effect after the filling of basement rock grout curtain, provides a kind of accurate numerical value means for prediction before the filling of Dam Heavy Curtain Grouting or filling post analysis.
Description
Technical field
The invention belongs to dam foundation grout curtain numerical simulation fields in hydraulic and hydroelectric engineering, more specifically, being to relate to
And a kind of grout curtain method for numerical simulation that diffusion sharp side is solved based on cement grout thixotropy and VOF.
Background technique
Grout curtain is the main means of dam foundation Anti-seeping technology, is risen emphatically to the safe operation for guaranteeing hydraulic structure
It acts on.Since seepage proof curtain is fully located at underground, there is concealment, complex geologic conditions, lead to slurries in the earth formation in addition
Motion process it is also extremely complex, bring difficulty for accurate description law of grouting fluid movement, continuity to curtain after construction and
Integrality grouting effect is difficult to be assessed.Therefore numerical simulation and analysis are carried out to grout curtain process, discloses geological condition
Relationship between size performance, determine grout curtain slurries sharp side overlap closure situation, to instruct Grouting Design and construction have
There are important theory significance and engineering application value.
Foreign countries are less about the research of numerical simulation of grouting[1], method is mostly using UDEC software simulation tunnel crack note
Slurry.UDEC is the two-dimensional discrete metaprogram of a processing discontinuous media, therefore numerical analysis is limited only to two dimension, while can not
Consider the rheological characteristic of the geologic(al) factor and cement grout near grout hole.
In country's grouting analog study, Lee's grain guiding principle etc.[2]Using the method for numerical simulation based on UDEC discrete unit program,
The flowing diffusion of diffusion profile and slurries to deep mine drilling grout curtain pressure is simulated;Yuan Jing is strong etc.[3]It is situated between based on granular media
The particle stream method of matter theory, with PFC2DSoftware has carried out mesomechanics simulation to soft stratum slip casting process, has studied not
With slip casting element to the affecting laws of grout spreading range and slip casting type;Liu Jian etc.[4]It is ground using comsol hydrodynamics module
Diffusion Law of the cement grout in hydrostatic and river channels lower plane crack is studied carefully;Han Lei etc.[5]With two-phase flow module in FLAC
Two-dimensional analog analysis has been carried out to the filling process of one 5m of certain earth and rockfill dam grouting hole section.The above method is two-dimension numerical calculation,
Geological conditions and the slurries rheological equationm of state are not considered.Wang Xiaoling etc.[6]Using mixed model to Bingham liquid slurries in goaf
Turbulent diffusion be simulated calculating.Deng Shaohui etc.[7]Regard cement grout as single homogeneous Bingham liquid to slurries big
Permeation grouting in dam foundation rock is simulated research.Though this seminar realizes grouting three-dimensional numerical value research, slurries are seen
Make Bingham liquid, the rheological equationm of state that slurries property is more met can be expressed not as good as thixotropy model, and be to be related to starching
The detailed calculation method on liquid diffusion sharp side.
In conclusion existing research majority is two-dimension numerical calculation both at home and abroad, and geological conditions and fluid properties are not considered;
Regard cement grout as Bingham liquid in extremely individual three-dimensional numerical value researchs, more accurately expresses cement slurry not as good as thixotropy model
Liquid rheological property, and be to be related to the calculation method on slurries diffusion sharp side.It is counted currently based on cement grout thixotropy model and VOF
The three-dimensional grouting numerically modeling for calculating the considerations of determining slurries diffusion sharp side complex geological condition has not been reported.
[1]Yong Xu,Weihong Peng.Research on Multiple Holes Grouting of
Fractured Rock Mass[J].Applied Mechanics and Materials,2013,547-551.
[2] Li Lianggang, Tang Ping, He Weishan wait the numerical simulation and application [J] mine locating work of the drilling grout curtain of depth mine
Journey, 2010,37 (12): 36-40.
[3] Yuan Jingqiang, Chen Weizhong, Tan Xianjun wait mesomechanics analog study [J] ground power of soft stratum slip casting
It learns, 2011,32 (S2): 653-659.
[4] Liu Jian, Liu Rentai, Zhang Xiao wait cement grout crack grouting Diffusion Law model test and numerical simulation [J]
Chinese Journal of Rock Mechanics and Engineering, 2012,31 (12): 2445-2452.
[5] research [J] rock-soil mechanics of Han Lei, Chen Jiansheng, Chen Liang grout curtain dilation angle and numerical simulation,
2012,33 (7): 2235-2240.
[6] Wang Xiaoling, Wang Qingsong, straight print, Ao Xuefei project of South-to-North water diversion goaf three-dimensional Bingham liquid turbulent flow mould
Quasi- [J] Journal of Hydraulic Engineering, 2013,11:1295-1302.
[7] the three of Deng Shaohui, Wang Xiaoling, Ao Xuefei, Ren Ping Yu, Li Ruijin dam foundation curtain Binghamian grouts grouting
Dimensional Simulation of Air Motion [J] Journal of Hydraulic Engineering, 2016,02:165-172+179.
Summary of the invention
Purpose of the invention is to overcome the shortcomings in the prior art, provide it is a kind of based on cement grout thixotropy and
VOF solves the grout curtain method for numerical simulation on diffusion sharp side, realizes under complex geological condition based on thixotropy model and VOF
The three-dimensional numerical simulation of grouting for determining slurries diffusion sharp side is calculated, accurately learns the filling of basement rock grout curtain under complex geological condition
Macro-effect afterwards provides a kind of accurate numerical value means for post analysis is predicted or filled before the filling of Dam Heavy Curtain Grouting.
The purpose of the present invention is what is be achieved through the following technical solutions.
The grout curtain method for numerical simulation that diffusion sharp side is solved based on cement grout thixotropy and VOF of the invention, packet
Include following steps:
(1) three-dimensional geological model and grouting pore model are established, three-dimensional geological model is coupled and grout hole model construction is three-dimensional
Unified model;
(2) the three dimensional unification model in step (1) is imported into CFD software for calculation, divides and calculates grid;
(3) cement grout thixotropic fluid model is established, VOF method captures slurries and spreads sharp side;
(4) it realizes and considers that slurries thixotropy and VOF calculate the dam foundation on slurries diffusion sharp side under complex geological condition
Grout curtain Three-dimensional simulation.
The foundation of three-dimensional geological model in the step (1): terrain modeling is carried out to earth's surface, utilizes the original money of exploration
These discrete topographic contours are formed the sealene triangle of interconnection with TIN model approach by material by certain regular partition
Network carries out dimensional topography reconstruct to earth's surface TIN model with NURBS technology, obtaining dimensional topography can using TIN model as intermediary
Depending on changing model.
The foundation of grouting pore model in the step (1): the data variable collection by establishing mesoscale grout hole, foundation
" point-line-face-body " modeling method realizes the building of grout hole threedimensional model.
The foundation of three dimensional unification model in the step (1): three-dimensional geological model and grouting pore model are subjected to boolean's meter
It calculates, grouting data is converted to three-dimensional space model on the basis of three-dimensional geological model, make three-dimensional geological model and grout hole
Model correctly couples.
Cement grout thixotropic fluid model in the step (3) are as follows:
In formula,For apparent viscosity, η0For the plastic viscosity of thixotropic fluid, τ0It is answered for the surrender of thixotropic fluid
Power, λ0It is the structural parameters of material,For shear strain rate, t is the time, and α is the dimensionless group for indicating material properties.
It is VOF method by solving each moment calculating grid that VOF method, which captures slurries diffusion front mask body, in the step (3)
Fluid volume determines the location and shape on slurries diffusion sharp side than function in unit, and the fluid volume compares function are as follows:
In formula, F is fluid volume than function, xiIndicate x, tri- different directions of y, z, uiIndicate x, tri- different directions of y, z
Speed.
Realize that consideration slurries thixotropy and VOF under complex geological condition calculate slurries and spread sharp side in the step (4)
Dam foundation grout curtain Three-dimensional simulation using continuity equation, couple VOF method biphase gas and liquid flow N-S equation, gas
Liquid two-phase mixtures density equation and mixing velocity equation, thixotropic fluid constitutive equation.
Compared with prior art, the beneficial effects brought by the technical solution of the present invention are as follows:
(1) present invention couples three-dimensional geological model and grouting pore model, provides one for the analysis of entire grouting works
A new visual angle, spatial distribution and geological conditions to grout hole form whole understanding, provide for grouting numerical analysis fine
Three-dimensional computations geometrical model;
(2) present invention is converted by the interface between three-dimensional geological model, grout hole model data and CFD model data,
It establishes CFD and calculates grid model, realize accurate expression of the true complex geological condition in CFD software, compensate for previous CFD and build
Drawing grid model by coordinate in mould makes the insufficient limitation of grid accuracy;
(3) present invention thixotropy model tormulation cement rheological equationm of state, it is contemplated that shearing rate and time in filling process
Two-phase factor can more reflect the time variation of cement grout than bingham model, capture curtain slurries using VOF method and spread sharp side,
Have the characteristics that amount of calculation is small, data storage capacity is small and is easily restrained to calculating;
(4) what the present invention established considers that slurries thixotropy and VOF calculate the curtain on slurries sharp side based on complex geological condition
Numerical simulation of grouting method realizes the porous point of sequence grouting fine analog of cement grout during curtain grouting in foundation of dam, simulation
Condition is more in line with the practical grouting condition in scene, can continuity to curtain after construction and integrality grouting effect effectively commented
Estimate.
Detailed description of the invention
Fig. 1 is the method flow diagram for coupling three-dimensional geological model and grout hole model construction three dimensional unification model;
Fig. 2 is to realize thixotropy slurries dam foundation grout curtain Three-dimensional simulation techniqueflow chart using CFD technology.
Specific embodiment
The invention will be further described with reference to the accompanying drawing.
The present invention proposes a kind of curtain filling that slurries thixotropy and VOF calculating slurries sharp side are considered based on complex geological condition
Starch Three-dimensional Numerical Simulation Method.Three-dimensional geological model and grout hole model construction three dimensional unification model are coupled first;Then pass through
Three dimensional unification model is imported CFD software, using polyhedron grid partitioning technology, to three dimensional unification model partition by data-interface
Fine grid blocks;The thixotropic fluid model that can really reflect its rheological property is finally established to cement grout, is calculated using VOF method
It captures slurries and spreads sharp side, realize and consider that slurries thixotropy and VOF calculate the big of slurries diffusion sharp side under complex geological condition
Dam foundation rock grout curtain Three-dimensional simulation.The curtain that diffusion sharp side is solved based on cement grout thixotropy and VOF of the invention
The techniqueflow of numerical simulation of grouting method is as depicted in figs. 1 and 2, in particular relates to include coupling three-dimensional geological model and filling
Starch pore model construction method, thixotropy modeling expression cement grout, VOF solve slurries diffusion sharp side method and based on
The dam foundation grout curtain numerical simulation of upper method, preferred forms are as follows:
One, three-dimensional geological model and grouting pore model are established, three-dimensional geological model is coupled and grout hole model construction is three-dimensional
Unified model, comprising:
(1) three-dimensional geological model is established, realization process can be summarized as follows: terrain modeling is carried out to earth's surface first, utilized
These discrete topographic contours are formed interconnection by certain regular partition with TIN model approach by the firsthand information such as contour
Irregular triangular network dimensional topography reconstruct is carried out to earth's surface TIN model with NURBS technology using TIN model as intermediary,
Obtain Three-dimensional Terrain Visualization model, Three-dimensional Terrain Visualization model be also used as all geologic body arithmetic operations space constraint and
Basis.1. establishing geology respectively using above-mentioned dimensional Modeling Technology according to the original geologic data of each geologic body itself
The upper and lower structural plane S1 and S2 of body;2. the spatial relation of structural plane S1 and S2 are utilized, respectively to receptor-earth's surface three-dimensional mould
Type carries out boolean operation, forms geologic body T1 and 3. finally carries out space topology analyzing to all unfavorable geologies.
(2) establish grouting pore model: mesoscale grouting works works, which refers to, to be distributed in rock mass space in practical grouting
Grouting carrier, according to " point-line-face-body " modeling method, realize and fill by establishing the data variable collection of mesoscale grout hole
Starch the building of hole threedimensional model.Comprehensive grout hole orifice position information (x, y, z), segment length's information L, orifice diameter D information, inclination angle
β information realizes the Visualization Modeling of grout hole.1. with the corresponding three-dimensional space point P1 of (x, y, z) Coordinate generation;2. P1 is
The center of circle, D/2 are radius, utilize rhinoceros Circle order (function ON_Circle::ON_Circle) formation curve C1;3. Z-direction
As axis, L is length, utilizes rhinoceros ExtrudeCurveStraight order (function
RhinoExtrudeCurveStraight the stretching setting-out for) realizing curve C1, to form curved surface S1;4. with rhinoceros Cap order
(function Cap) makes curved surface S1 form the constraint of space Z-direction, completes grout hole and turns from digital data to three-dimensional space model
Change, above step 1. -4. can be automatically performed with program.
(3) three-dimensional geological model and grout hole model construction three dimensional unification model are coupled: on the basis of three-dimensional geological model
On, realize grouting conversion of the data to three-dimensional space model, that is, three-dimensional geological model and grouting pore model are subjected to boolean
It calculates, grouting data is converted to three-dimensional space model on the basis of three-dimensional geological model, make three-dimensional geological model and grouting
Pore model correctly couples, and will provide a new visual angle for the analysis of entire grouting works.The three dimensional unification established using this method
Model provides fine computational geometry model basis for the dam foundation grout curtain CFD analysis of this research.
Two, three dimensional unification model is imported into CFD software for calculation, divides and calculates grid specifically:
(1) three dimensional unification model can be received compatible .stl format with CFD software to export, fluid (is calculated by CFD
Mechanics) software data-interface, by comprising truly quality shape and grouting hole data .stl formatted file imported into CFD software
In, realize truly expressed of the three dimensional unification model in CFD software.
(2) the three dimensional unification model based on importing establishes mould using polyhedron grid partitioning technology to different geologic provinces
Type: biggish size is selected in the overall situation, with zoning overlapping margins, the non-uniform grid of density degree is generated, so that net
The boundary and zoning boundary of lattice correspond;It is configured smaller size in grout hole regional area, makes mesh refinement, to the greatest extent
The geometry for maintaining grout hole may be retained, provide fine grid blocks model for the CFD calculating of grout curtain.
Three, cement grout thixotropic fluid model is established, VOF method captures slurries diffusion front mask body are as follows:
(1) cement grout thixotropic fluid model: fresh cement slurry shows as incompressible non-newtonian viscous fluid.In transient flow situation
Under, apparent viscosity is the function of shearing rate and time.In grouting works, reperfusion serosity becomes with grouting time property
Change, can be indicated with the thixotroping model in the case of transient flow.
The constitutive relation of thixotropic fluid are as follows:
In formula: τ is deviatoric stress tensor, Pa;For apparent viscosity, be λ andFunction, Pas;λ is the knot of material
Structure parameter,For shear strain rate;
To avoid numerical dissipation, thixotropic cement slurries apparent viscosity is improved are as follows:
In formula:For apparent viscosity, Pas;η0For the plastic viscosity of thixotropic fluid, Pas;τ0For thixotropy stream
The yield stress of body, Pa;λ0It is the structural parameters of material;For shear strain rate;T is the time;α is the nothing for indicating material properties
Dimensional parameters.
The above thixotropic fluid apparent viscosity (formula (2)) is programmed, exploitation, which is improved in CFD software, to be calculated,
It ensure that the truly expressed of fluid property.
(2) slurries diffusion sharp side is the interface of biphase gas and liquid flow, VOF model tracking gas-liquid interface, traditionally with one
A single game model, i.e. gas-liquid two-phase share a velocity field and density field.Basic ideas are as follows: by solve a set of equation of momentum and
Continuity equation simulates the movement for the fluid that two kinds do not blend mutually, tracks volume shared by every kind of fluid, the institute in each control volume
The volume function for having fluidised form is integer 1.VOF method is by fluid volume in each moment calculating grid cell of solution than function (stream
The ratio value function of the variation of body volume and grid cell volume itself) come determine slurries diffusion sharp side location and shape, institute
It states fluid volume and compares function are as follows:
In formula, F is fluid volume than function, xiIndicate x, tri- different directions of y, z, uiIndicate x, tri- different directions of y, z
Speed.
For the specified phase fluid in certain moment arbitrary mess unit, if F=0, show grid cell by other phases
Fluid is full of, that is, does not include specified phase fluid;If 0 < F < 1, show that there are the interfaces of multiphase substance in grid cell;If F
=1, show that grid cell is designated phase fluid and is completely filled with, is specified phase fluid unit.
Four, it realizes and considers that slurries thixotropy and VOF calculate the dam foundation on slurries diffusion sharp side under complex geological condition
Grout curtain Three-dimensional simulation specifically:
Fundamental equation includes continuity equation, biphase gas and liquid flow N-S equation, the gas-liquid two-phase hybrid density for coupling VOF method
Equation and mixing velocity equation, thixotropic fluid constitutive equation;Mixing velocity equation and hybrid density equation are to express seepage flow
The mixing velocity and hybrid density of gas-liquid two-phase in, mixing velocity equation and hybrid density equation pass through filtrational resistance source item
It is coupled in biphase gas and liquid flow N-S equation, the kinematic viscosity ν of the apparent viscosity η association N-S equation of thixotropic fluid is to express water
The thixotropic property of slurry liquid couples the N-S equation of VOF method to solve slurries diffusion sharp side.Wherein:
Continuity equation:
Biphase gas and liquid flow N-S equation:
In formula (4) and formula (5), φ is the porosity of porous media;U is the velocity vector of infinitesimal section, m/s, u and hole
The relationship of fluid true velocity u ' is u=φ u ';T is the time, and unit is s;F is unit mass force, and unit is N;P is that fluid is micro-
Pressure on first body, unit are Pa;ρ is fluid density, and unit is m3/s;ν is fluid kinematic viscosity, m2/ s, ν are apparent with fluid
The relationship of viscosities il is ν=η/ρ;S is filtrational resistance source item, is equal to inertia and loses item SiWith viscosity loss item SvThe sum of.For stream
Speed is lower to flow negligible inertia loss item, that is, thinks Si=0, viscosity loss item SvExpression formula are as follows:
In formula (6),For viscosity factor, unit is 1/m2, expression formula are as follows:
In formula (7), K is the infiltration coefficient of porous media, and unit is m/s.
After coupling VOF method, gas-liquid two-phase shares a set of equation, and ρ and ν of the formula (4) into formula (7) are no longer constants, but
The variable determined using the hybrid density equation and mixing velocity equation of gas-liquid two-phase, gas-liquid two-phase hybrid density equation and mixing
The expression formula of rate equation is respectively as follows:
Hybrid density equation ρ=Fgρg+(1-Fg)ρa (8)
Mixing velocity equation ν=Fgνg+(1-Fg)νa (9)
In formula (8) and formula (9): ρg、ρaRespectively the density of two kinds of slurries, gas fluids, unit are kg/m3;νg、νaRespectively
The coefficientof kinematic viscosity of two kinds of slurries, gas fluids, unit are m2/s;FgFor the volume fraction of slurries, zero dimension.
The constitutive relation of thixotropic fluid is introduced, i.e.,
In formula: τ is deviatoric stress tensor, Pa;For apparent viscosity, be λ andFunction, Pas;λ is the knot of material
Structure parameter,For shear strain rate;
To avoid numerical dissipation, thixotropic cement slurries apparent viscosity is improved are as follows:
In formula:For apparent viscosity, Pas;η0For the plastic viscosity of thixotropic fluid, Pas;τ0For thixotropy stream
The yield stress of body, Pa;λ0It is the structural parameters of material;For shear strain rate;T is the time;α is the nothing for indicating material properties
Dimensional parameters.
Design conditions include:
(1) condition for import: data are recorded according to the filling process of practical grouting works and calculate determining cement grout import
VELOCITY DISTRIBUTION, pressure distribution and corresponding volumetric concentration distribution, each variable, which takes, to be uniformly distributed, and inlet boundary condition is grout hole
The corresponding slurries VOF value for pouring into speed and import of the slurries on boundary;According to engineering design plan (EDP), grout hole is divided into I, II, III
Sequence hole, self-compiling program realize the inlet boundary condition temporally changed, and point sequence for carrying out I, II, III sequence hole is in the milk.
(2) exit condition: consider that dam foundation bottom is there may be lifting phenomenon in filling process, therefore take basement rock upper surface
For pressure export boundary.On the boundary of computational domain, pressure should meet Neumann condition, i.e. second kind boundary condition.
(3) wall boundary condition: two lateral boundaries of dam foundation upstream and downstream and bottom boundary are fricton-tight wall surface boundary, by solid
Law of the wall processing, uses non-slip condition: u=v=w=0 on wall surface, and wherein u, v, w are the speed in x, y, z direction.
(4) according to the geologic parameter of the geological stratifications such as packsand, siltstone, silty and argillaceous siltstoue, using not
Same porosity and resistance coefficient.
According to the continuity equation of above-mentioned determination, the biphase gas and liquid flow N-S equation for coupling VOF method, gas-liquid two-phase mixing speed
The design conditions for spending equation and hybrid density equation, thixotropy slurries constitutive equation and above-mentioned setting, pass through finite volume method
Carry out it is discrete, using SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm into
Row solves, and is expressed using thixotropy model cement grout fluid properties, finely captures slurries diffusion front using VOF method
Face, the porous point of sequence grouting of thixotropic cement slurries is fine during realizing curtain grouting in foundation of dam with CFD technology self-compiling program
Simulation analyzes curtain overlap joint closure situation and thixotropic cement slurries scattering and permeating process after the completion of grouting;It carries out simultaneously
The comparative study of grouting quantity and the analogue value of grouting time and field measurement value.
Although function and the course of work of the invention are described above in conjunction with attached drawing, the invention is not limited to
Above-mentioned concrete function and the course of work, the above mentioned embodiment is only schematical, rather than restrictive, ability
The those of ordinary skill in domain under the inspiration of the present invention, is not departing from present inventive concept and scope of the claimed protection situation
Under, many forms can also be made, all of these belong to the protection of the present invention.
Claims (6)
1. solving the grout curtain method for numerical simulation on diffusion sharp side based on cement grout thixotropy and VOF, which is characterized in that packet
Include following steps:
(1) three-dimensional geological model and grouting pore model are established, three-dimensional geological model and grout hole model construction three dimensional unification are coupled
Model;
(2) the three dimensional unification model in step (1) is imported into CFD software for calculation, divides and calculates grid;
(3) cement grout thixotropic fluid model is established, VOF method captures slurries and spreads sharp side;
Wherein, cement grout thixotropic fluid model are as follows:
In formula,For apparent viscosity, η0For the plastic viscosity of thixotropic fluid, τ0For the yield stress of thixotropic fluid, λ0It is
The structural parameters of material,For shear strain rate, t is the time, and α is the dimensionless group for indicating material properties;
(4) it realizes and considers that slurries thixotropy and VOF calculate the dam foundation curtain on slurries diffusion sharp side under complex geological condition
Be in the milk Three-dimensional simulation.
2. the grout curtain Numerical-Mode according to claim 1 for solving diffusion sharp side based on cement grout thixotropy and VOF
Quasi- method, which is characterized in that the foundation of three-dimensional geological model in the step (1): terrain modeling is carried out to earth's surface, utilizes exploration
Firsthand information, these discrete topographic contours are formed by certain regular partition by not advising for interconnection with TIN model approach
Then triangular net carries out dimensional topography reconstruct to earth's surface TIN model with NURBS technology, obtains three using TIN model as intermediary
Tie up terrain visualization model.
3. the grout curtain Numerical-Mode according to claim 1 for solving diffusion sharp side based on cement grout thixotropy and VOF
Quasi- method, which is characterized in that the foundation of grouting pore model in the step (1): the data by establishing mesoscale grout hole become
Quantity set realizes the building of grout hole threedimensional model according to " point-line-face-body " modeling method.
4. the grout curtain Numerical-Mode according to claim 1 for solving diffusion sharp side based on cement grout thixotropy and VOF
Quasi- method, which is characterized in that the foundation of three dimensional unification model in the step (1): by three-dimensional geological model and grouting pore model
Boolean calculation is carried out, grouting data are converted to three-dimensional space model on the basis of three-dimensional geological model, make three-dimensional geological mould
Type and grouting pore model correctly couple.
5. the grout curtain Numerical-Mode according to claim 1 for solving diffusion sharp side based on cement grout thixotropy and VOF
Quasi- method, which is characterized in that when VOF method capture slurries diffusion front mask body is that VOF method is each by solution in the step (3)
It carves and calculates the location and shape that fluid volume in grid cell determines slurries diffusion sharp side than function, the fluid volume compares letter
Number are as follows:
In formula, F is fluid volume than function, xiIndicate x, tri- different directions of y, z, uiIndicate x, the speed of tri- different directions of y, z
Degree.
6. the grout curtain Numerical-Mode according to claim 1 for solving diffusion sharp side based on cement grout thixotropy and VOF
Quasi- method, which is characterized in that realize that consideration slurries thixotropy and VOF under complex geological condition calculate slurry in the step (4)
Liquid spreads the dam foundation grout curtain Three-dimensional simulation on sharp side using continuity equation, the biphase gas and liquid flow of coupling VOF method
N-S equation, gas-liquid two-phase hybrid density equation and mixing velocity equation, thixotropic fluid constitutive equation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610886410.XA CN106446439B (en) | 2016-10-11 | 2016-10-11 | The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610886410.XA CN106446439B (en) | 2016-10-11 | 2016-10-11 | The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106446439A CN106446439A (en) | 2017-02-22 |
CN106446439B true CN106446439B (en) | 2019-05-24 |
Family
ID=58173864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610886410.XA Active CN106446439B (en) | 2016-10-11 | 2016-10-11 | The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106446439B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107401403B (en) * | 2017-09-06 | 2023-10-10 | 重庆科技学院 | Shale gas well multistage fracturing cement ring airtight complete visual evaluation device and method |
CN108629126B (en) * | 2018-05-09 | 2020-04-21 | 中国地质大学(北京) | Rock mechanics numerical modeling method considering macro-micro defect coupling |
CN109001078B (en) * | 2018-06-12 | 2021-02-09 | 北京城建道桥建设集团有限公司 | Method for establishing constitutive relation of II-type cement-based grouting material |
CN109859312B (en) * | 2019-02-01 | 2023-03-28 | 上海勘察设计研究院(集团)有限公司 | Refined three-dimensional geological model modeling method based on BIM technology |
CN110516366B (en) * | 2019-08-28 | 2023-04-07 | 北京工业大学 | Modeling method based on random microbead micromixer in ultra-high performance liquid chromatography analysis |
CN111735937A (en) * | 2020-07-27 | 2020-10-02 | 广东博智林机器人有限公司 | Method and device for detecting whether grouting is full in grouting equipment |
CN112069566B (en) * | 2020-08-04 | 2023-08-25 | 清华大学建筑设计研究院有限公司 | Method for finding patterns of cultural relic building protection shed |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034765A (en) * | 2012-12-14 | 2013-04-10 | 天津大学 | Goaf grouting dynamic overall process simulation method based on numerical simulation |
CN103984807A (en) * | 2014-04-29 | 2014-08-13 | 天津大学 | Numerical simulation method of dam grouting capable of coupling fine geological information and monitoring information |
CN104007045B (en) * | 2014-05-12 | 2016-03-23 | 河海大学 | A kind of slurry shield machine mud film forming method for numerical simulation |
-
2016
- 2016-10-11 CN CN201610886410.XA patent/CN106446439B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034765A (en) * | 2012-12-14 | 2013-04-10 | 天津大学 | Goaf grouting dynamic overall process simulation method based on numerical simulation |
CN103984807A (en) * | 2014-04-29 | 2014-08-13 | 天津大学 | Numerical simulation method of dam grouting capable of coupling fine geological information and monitoring information |
CN104007045B (en) * | 2014-05-12 | 2016-03-23 | 河海大学 | A kind of slurry shield machine mud film forming method for numerical simulation |
Non-Patent Citations (2)
Title |
---|
大坝基岩帷幕宾汉姆浆液灌浆的三维数值模拟;邓韶辉 等;《水利学报》;20160229;第165-172,179页 |
硅酸盐水泥--阳离子乳化沥青浆体的流变模型;彭建伟 等;《硅酸盐学报》;20140531;第42卷(第5期);第621-628页 |
Also Published As
Publication number | Publication date |
---|---|
CN106446439A (en) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106446439B (en) | The grout curtain method for numerical simulation on diffusion sharp side is solved based on cement grout thixotropy and VOF | |
CN102339326B (en) | Method for analyzing and simulating fluid flow of fracture-cavity oil reservoir | |
Blunt | Effects of heterogeneity and wetting on relative permeability using pore level modeling | |
CN107330191A (en) | The Numerical Analysis methods of Groundwater iron shield tunnel construction influence | |
Mariño et al. | Seepage and groundwater | |
CN104318015B (en) | Prediction earth and rockfill dam nodal region seepage field and the method for determining earth and rockfill dam dam body saturated surface | |
CN107060746A (en) | A kind of method of complex fracture oil deposit flow simulation | |
Acharya et al. | Porosity–permeability properties generated with a new 2-parameter 3D hydraulic pore-network model for consolidated and unconsolidated porous media | |
Chen et al. | Composite element method for the seepage analysis of rock masses containing fractures and drainage holes | |
Yao et al. | Fractured vuggy carbonate reservoir simulation | |
CN103034765A (en) | Goaf grouting dynamic overall process simulation method based on numerical simulation | |
CN103136436B (en) | Consider that the nucleic disperse and transport of adsorption of sediment analyzes method | |
Deng et al. | Simulation of grouting process in rock masses under a dam foundation characterized by a 3D fracture network | |
Wang et al. | Formation of a polymer thin wall using the level set method | |
Thorne | Processes of bank erosion in river channels | |
Xiang et al. | Seepage analysis of the fractured rock mass in the foundation of the main dam of the Xiaolangdi water control project | |
Li et al. | The theoretical and numerical analysis of water inrush through filling structures | |
Xu et al. | Numerical investigation of flow control technology for grouting and blocking of flowing water in karst conduits | |
CN113657048B (en) | Debris flow rainfall converging water source replenishment measuring and calculating method, debris flow dynamic process numerical simulation method and application | |
Xu et al. | A grouting simulation method for quick-setting slurry in karst conduit: The sequential flow and solidification method | |
Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
Zhang et al. | Multifield coupling numerical simulation of the seepage and stability of embankment dams on deep overburden layers | |
Yang et al. | Numerical investigation on the spewing mechanism of earth pressure balance shield in a high‐pressure water‐rich sand stratum | |
Ilyinsky et al. | Analytical solutions of seepage theory problems. Inverse method, variational theorems, optimization and estimates (a review) | |
Tao et al. | Grouting fractured coal permeability evolution based on industrial CT scanning |
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 |