CN109033537A - The calculation method and system of rock-fill concrete casting process numerical simulation - Google Patents
The calculation method and system of rock-fill concrete casting process numerical simulation Download PDFInfo
- Publication number
- CN109033537A CN109033537A CN201810699310.5A CN201810699310A CN109033537A CN 109033537 A CN109033537 A CN 109033537A CN 201810699310 A CN201810699310 A CN 201810699310A CN 109033537 A CN109033537 A CN 109033537A
- Authority
- CN
- China
- Prior art keywords
- stone
- heap
- final
- analysis
- original state
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the present invention provides the calculation method and system of a kind of rock-fill concrete casting process numerical simulation, calculation method includes: the banking process that analysis heap stone original state threedimensional model is calculated according to finite element/discrete element coupling analytical method FEM/DEM and final accumulation form；After the banking process and final accumulation form for calculating analysis heap stone original state threedimensional model, it is based on smoothed particle method SPH method, calculates casting process and final occupied state of the analysis self-compacting concrete in the final accumulation form of heap stone.The calculation method and system of rock-fill concrete casting process numerical simulation provided in an embodiment of the present invention, it is capable of the casting process of accurate simulation self-compacting concrete, it can not only be visually observed that flow regime of the self-compacting concrete in enrockment hole, it can also be appreciated that the dense state that rock-fill concrete is final, a large amount of human cost, experimental cost and economic cost are saved, and avoids the problem that experimental method is difficult to online observation.
Description
Technical field
The present embodiments relate to rock-fill concrete pouring technique fields, pour more particularly, to a kind of rock-fill concrete
Build the calculation method and system of numerical simulation.
Background technique
Mass concrete occupies an important position in modern construction, especially hydraulic and hydroelectric engineering construction.China
Often it is only the mass concrete poured in hydraulic and hydroelectric engineering just at 10,000,000 cubes or more, in addition, harbour, airport are built
It builds and heavy-duty machine basis etc. is also tended to using mass concrete.(partial size is big by using big gauge block stone for rock-fill concrete
In 30cm) cement consumption is reduced, to effectively reduce the heat of hydration and reduce CO2Discharge, small with concrete structure shrinkage,
The features such as cracking resistance shear resistance improves, speed of application is fast, the heat of hydration is low, and temperature control is easy, construction quality is high, project cost is low, because
And it has broad application prospects, while also meeting the urgent need that China widelys popularize green low-carbon technology.
Rock-fill concrete in the construction process, the banking process of rockfill and last accumulation form, to heap after molding
The mechanical property of stone concrete has large effect, but due to the form pattern of rockfill in experiment cost height and Practical Project
Complexity, therefore be difficult to carry out more accurate experiment and Practical Project is analyzed, therefore carry out numerical computation method to enrockment
It is inevitable trend that concrete rockfill banking process, which carries out analytical calculation,.But up to the present, how accurately to simulate calculating
The pouring construction process of rock-fill concrete, there has been no any achievements and method for reference.
Summary of the invention
In view of the problems of the existing technology, the embodiment of the present invention provides a kind of rock-fill concrete casting process numerical simulation
Calculation method and system.
The embodiment of the present invention provides a kind of calculation method of rock-fill concrete casting process numerical simulation, comprising: according to having
Member/discrete element coupling analytical method FEM/DEM is limited, the banking process of analysis heap stone original state threedimensional model and final is calculated
Accumulation form；After calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation form, it is based on
Smoothed particle method SPH method calculates analysis self-compacting concrete pouring in the final accumulation form of heap stone
Process and final occupied state.
The embodiment of the present invention provides a kind of computing system of rock-fill concrete casting process numerical simulation, comprising: heap integral
Module is analysed, for calculating analysis heap stone original state three-dimensional mould according to finite element/discrete element coupling analytical method FEM/DEM
The banking process of type and final accumulation form；Analysis module is poured, for calculating the analysis heap stone original state three-dimensional
After the banking process of model and final accumulation form, it is based on smoothed particle method SPH method, it is mixed to calculate analysis self-compaction
Casting process and final occupied state of the solidifying soil in the final accumulation form of heap stone.
The embodiment of the present invention provides a kind of calculating equipment of rock-fill concrete casting process numerical simulation, comprising: at least one
A processor；And at least one processor being connect with the processor communication, in which: the memory is stored with can be by institute
The program instruction of processor execution is stated, the processor calls described program instruction to be able to carry out above-mentioned calculation method.
The embodiment of the present invention provides a kind of non-transient computer readable storage medium, the non-transient computer readable storage
Medium storing computer instruction, the computer instruction make the computer execute above-mentioned calculation method.
The calculation method and system of rock-fill concrete casting process numerical simulation provided in an embodiment of the present invention, pass through setting
Finite element and discrete element coupling analytical method carry out mechanical analysis, and it is progressive under stress can completely to depict heap stone
Banking process, calculated result are more accurate.Self-compacting concrete is analyzed in the final accumulation shape of heap stone by setting SPH method
Casting process and final occupied state in state, are capable of the casting process of accurate simulation self-compacting concrete, not only can be intuitive
Flow regime of the self-compacting concrete in enrockment hole is observed on ground, can also understand the final closely knit shape of rock-fill concrete
State.Present invention saves a large amount of human cost, experimental cost and economic costs, and avoid experimental method and be difficult to see online
The problem of survey.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the flow chart of the calculation method embodiment of rock-fill concrete casting process numerical simulation of the present invention；
Fig. 2 is the module map of the computing system embodiment of rock-fill concrete casting process numerical simulation of the present invention；
Fig. 3 is the frame signal of the calculating equipment of the rock-fill concrete casting process numerical simulation in the embodiment of the present invention
Figure.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Fig. 1 is the flow chart of the calculation method embodiment of rock-fill concrete casting process numerical simulation of the present invention, such as Fig. 1 institute
Show, comprising: it is three-dimensional to calculate analysis heap stone original state by S101, according to finite element/discrete element coupling analytical method FEM/DEM
The banking process of model and final accumulation form；S102, the accumulation for analyzing the heap stone original state threedimensional model is being calculated
After process and final accumulation form, it is based on smoothed particle method SPH method, calculates analysis self-compacting concrete in enrockment
Casting process and final occupied state in the final accumulation form of block.
It should be noted that due to the image of scale effect in physical model experiment, the logistics organizations of existing model dimension
It is unable to the intensity of correct response real structure, while prototype logistics organizations are limited by the limitation of experiment condition again.By that will simulate
Structural stress and the Finite Element of deformation and the discrete element method that bulk movement can be tracked combine referred to as finite element/
Discrete element coupling analytical method FEM/DEM, this method remain finite element and the respective advantage of discrete element, are able to solve more body fortune
Dynamic knowledge topic and mechanics problem is destroyed, during calculating the final accumulation form of heap stone, there are block stone body number is more, by
The characteristics of power is complicated, rockfill can generate deformation with stress variation, carries out power using finite element and discrete element coupling analytical method
Credit analysis, can completely depict progressive banking process of the heap stone under stress, calculated result is more accurate.
SPH (Smoothed Particle Hydrodynamics) method is the contracting of Smoothed Particle Hydrodynamics Method
It writes, is in a kind of non-mesh method gradually to grow up over nearly more than 20 years.The basic thought of this method is by continuous fluid
(or solid) is described with the Particle Group of interaction, carries various physical quantitys, including quality, speed etc. on each substance point,
By solving the kinetics equation of Particle Group and tracking the tracks of each particle, the mechanical behavior of whole system is acquired.This
Inventive embodiments are using SPH method casting process of the simulation self-compacting concrete in rockfill and and calculate final fluidised form.
The calculation method of rock-fill concrete casting process numerical simulation provided in an embodiment of the present invention, by the way that finite element is arranged
Mechanical analysis is carried out with discrete element coupling analytical method, progressive under stress of heap stone can be completely depicted and accumulate
Journey, calculated result are more accurate.Self-compacting concrete is analyzed in the final accumulation form of heap stone by setting SPH method
Casting process and final occupied state are capable of the casting process of accurate simulation self-compacting concrete, not only can intuitively be observed
To flow regime of the self-compacting concrete in enrockment hole, the final dense state of rock-fill concrete can also be understood.This hair
It is bright to save a large amount of human cost, experimental cost and economic cost, and avoid the difficulty that experimental method is difficult to online observation
Topic.
Based on the above embodiment, described according to finite element/discrete element coupling analytical method FEM/DEM, calculate analysis enrockment
The banking process of block original state threedimensional model and final accumulation form, before further include: it is three-dimensional to establish heap stone original state
Model, the heap stone include several block stone bodies.
Based on the above embodiment, described to establish heap stone original state threedimensional model, it specifically includes: soft using three-dimensional modeling
Part such as GID establishes heap stone original state threedimensional model according to the predetermined parameter of heap stone；Alternatively, to described block of stone
Body appearance profile carries out 3D and scans to form input data and be input in 3 d modeling software, establishes heap stone original state three-dimensional mould
Type.
Since the stone size, type, the quantity that need in each engineering or experiment are different, so according to heap
The actual parameter of stone establishes threedimensional model, has clearly a need for acquiring several actual parameters before foundation, such as the hardness, big of stone
Small, density etc..
The embodiment of the present invention can specifically use 3D scanner or profile measurement in the appearance profile of collection block stone body
Instrument is scanned, and obtains the formal parameter of block stone body, and the embodiment of the present invention is not specifically limited.
The present embodiment helps to obtain when subsequent progress finite element/discrete element coupling analytical method analysis banking process more quasi-
Really, true result.
Based on the above embodiment, the banking process for analyzing the heap stone original state threedimensional model and final of calculating
Accumulation form and it is described be based on smoothed particle method SPH method, calculate analysis self-compacting concrete in heap stone
Casting process and final occupied state in final accumulation form, between further include: establish self-compacting concrete model, it is described from
Air-free concrete model is Bingham rheological model.
Specifically, self-compacting concrete model are as follows:
Wherein, τ0For yield strength, μpFor plastic viscosity, τ is shear stress；γ is shear strain rate.
It should be noted that Bingham model (Bingham model): buffer, sliding block is in parallel, then connect with spring and
At reflection ground stickiness, elasticity, plastic flow feature and its process mechanical model.
Based on the above embodiment, described to be based on smoothed particle method SPH method, calculate analysis self-compacting concrete
Casting process and final occupied state in the final accumulation form of heap stone, specifically include: be based on SPH method, one by one when
Between step iteration update density, the position and speed of the fluid particles in the self-compacting concrete model；Based on the fluid
Density, the position and speed of particle analyze casting process of the self-compacting concrete in the final accumulation form of heap stone and most
Whole occupied state.
Function and equation needed for the present embodiment are made below and being further explained, the SPH in the embodiment of the present invention from
The form of dissipating are as follows:
Wherein, in which: N is the number that adjacent particles are total in support region；xjFor the position of particle j；mjAnd ρjRespectively particle
Quality, the density of j；Kernel function W (x-xj, h) and related with interparticle distance and smooth length h.
The present embodiment chooses kernel function are as follows:
Wherein, d is space dimensionality；H is smooth length.
Wherein, C is normaliztion constant.
The state equation that the present embodiment is chosen are as follows:
Wherein, B is as initial pressure, and for limiting the maximum change amount of density, the present invention is takenρ0It is particle
Initial density；γ is a constant coefficients；C is artificial velocity of sound, considers for numerical value computational stability and time step, usually takes
One value more much smaller than true velocity of sound generally takes 10 times of flow field maximum flow rate value, i.e.,H is that free surface is high
Degree.
The displacement relation formula of particle are as follows:
Wherein, r is the position vector of particle point, and u is the velocity vector of particle.
Boundary processing method in the embodiment of the present invention are as follows: Gu wall particle participates in the calculating of governing equation, but their positions
Fixed, displacement is set as 0 in the present invention.Lennard- generally can be used in interaction between boundary particle and fluid particles
Jones repulsive force method assumes when fluid particles distance Gu Bi is close enough, to prevent internal flow particle from passing through the side Gu Bi
It disperses behind boundary and causes to calculate collapse, the boundary particle fluid particles close to row apply a center repulsive force, repulsive force size
It is determined by mutual distance, projection speed, height and position etc..Free Surface is realized using the algorithm of the propositions such as Koshizuka
Dirichlet condition, interface particle pressure are assigned a value of 0 or certain external pressure force value.
Based on the above embodiment, the iteration of time step one by one updates the fluid grain in the self-compacting concrete model
Density, the position and speed of son, specifically include: time step one by one calculates the close of the fluid particles by continuity equation
Degree variation, is captured from by surface particle and carries out the density correction of Free Surface particle, obtain the self-compacting concrete model
In each moment density field；Time step one by one updates the rail of each fluid particles movement by fluid particles displacement equation
Mark obtains the position at any fluid particles each moment；Time step one by one determines any fluid particles by state equation
Pressure, the acceleration that calculation of pressure external force based on each fluid particles generates, the equation of momentum is passed through based on the acceleration based on
Calculate the velocity field at each moment in the self-compacting concrete model.
It should be noted that particle search is using tree search method or association chained list search technique.
The calculation method of rock-fill concrete casting process numerical simulation provided in an embodiment of the present invention, when by being arranged one by one
Between step iteration update density, the position and speed of the fluid particles in the self-compacting concrete model, can accurately calculate
The casting process of rock-fill concrete.
Based on the above embodiment, described according to finite element/discrete element coupling analytical method FEM/DEM, calculate analysis enrockment
The banking process of block original state threedimensional model and final accumulation form, specifically include: obtaining institute according to limited element analysis technique FEM
State each piece of stone body Continuum Mechanics behavior in heap stone original state threedimensional model；It is obtained according to DEM analysis method DEM
Non-continuum mechanics behavior in the heap stone original state threedimensional model between block stone body；By continuous Jie of all pieces of stone bodies
Matter mechanical behavior and non-continuum mechanics behavior substitute into the power balance equation of heap stone original state threedimensional model, obtain institute
State banking process and final accumulation form.
It should be noted that finite element analysis (FEA, Finite Element Analysis) is to utilize mathematical approach
Method simulates actual physical system (geometry and load working condition).Utilize simple and interaction element, so that it may use
The unknown quantity of limited quantity goes to approach the real system of unlimited unknown quantity.It is regarded as domain is solved by many referred to as finite elements
Small interconnection subdomain composition, assumes a suitable approximate solution to each unit, what then this domain of Derivation was total meets item
Part, to obtain the solution of problem.
Continuum Mechanics (Continuum mechanics) is a branch in physics, is processing including solid
The mechanics of so-called " continuous media " macroscopic property including body and fluid.For example, the conservation of mass, momentum and angular momentum theorem,
Conservation of energy etc..Continuum Mechanics thinks that the space that real fluid or solid are occupied can approximatively be regarded as continuously without sky
Unoccupied place is filled with " particle " (i.e. microcosmic upper sufficiently big, macroscopically sufficiently small micel).Macroscopic physical quantity possessed by particle
(such as quality, speed, pressure, temperature) meets all physical laws that should be followed, such as mass conservation law, newtonian motion
The transport properties such as law, law of conservation of energy, the law of thermodynamics and diffusion, viscosity and heat transfer.This hypothesis ignores substance
Specific microstructure (scope for belonging to condensed state physics to the Study on Microstructure of solid), and with one group of partial differential equation
To express macroscopic physical quantity (such as quality, several times, pressure etc.).
DEM analysis method is the method for numerical simulation for being specifically used to solve the problems, such as Discontinuous transmission.This method is heap stone
It is stereoscopic to allow block stone body to translate, rotation and deformation to be made of the joint plane between discrete block stone body and block stone body, and joint
Face can be compressed, separates or slide.Therefore, heap stone is counted as a kind of discontinuous discrete media.Big position may be present inside it
Shifting, rotation and sliding or even the separation of block, so as to more accurately simulate the Nonlinear Large Deformation feature in heap stone.
The general solution procedure of distinct element method are as follows: by solution room it is discrete be discrete element unit battle array, and according to practical problem
Adjacent Unit two is connected with reasonable connecting element；Relative displacement is basic variable between unit, by power and relative displacement
Relationship can be obtained Unit two between normal direction and tangential active force；Between the unit active force other units in all directions
And outer caused to unit effect of other physical fields makes every effort to resultant force and resultant moment, it can be in the hope of according to newton second law of motion
Obtain the acceleration of unit；Time integral is carried out to it, and then obtains the speed and displacement of unit.It is in office to obtain all units
The physical quantitys such as speed, acceleration, angular speed, displacement of the lines and the corner at meaning moment.
On the basis of the above embodiments, the present embodiment obtains institute according to limited element analysis technique FEM to further explain
State the detailed process of each piece of stone body Continuum Mechanics behavior in heap stone original state threedimensional model.
Specifically, the Continuum Mechanics behavior includes: interior force vector, grid caused by the column joints deformation of grid cell
The outer force vector and mass matrix of the node of unit.
Correspond ground, the continuous media power that each piece of stone body in threedimensional model is obtained according to limited element analysis technique
Scholarship and moral conduct is, specifically:
Grid dividing is carried out to heap stone original state threedimensional model, obtains the finite element model being made of grid cell.
The interior force vector is obtained according to grid cell volume, standard type function and cauchy stress tensor.
According to volume of grid cell when undeformed, density corresponding with volume of grid cell when undeformed and
The standard type function obtains the mass matrix.
According to the physical strength of block stone body, face power and cell list area, the outer force vector is obtained.
Further, interior force vector is calculated as follows:
Wherein, V indicates unit volume, and N is standard shape function, and T is cauchy stress tensor.
The outer force vector of node is calculated as follows:
Wherein, b indicates physical strength, and t indicates face power, and S is cell list area.
Mass matrix is calculated as follows:
Wherein, V0For volume of unit when undeformed, ρ0Expression and V0Corresponding density.
On the basis of the above embodiments, the present embodiment obtains institute according to DEM analysis method DEM to further explain
State the detailed process of the non-continuum mechanics behavior in heap stone original state threedimensional model between block stone body.
Specifically, non-continuum mechanics behavior includes: the contact force between block stone body.
Correspond ground, the Discontinuous transmission power obtained in threedimensional model between block stone body according to DEM analysis method
Scholarship and moral conduct is, specifically:
Searching for Contact Algorithm by non-binary trees, (contact judgment method then uses non-two that Munjiza and Andrews is proposed
Fork tree search Contact Algorithm, i.e., efficient NBS algorithm, this algorithm committed memory is small, for the feelings of Rigid Body Element adjoining dimensions
Condition may be implemented the Rapid contact between unit search) obtain threedimensional model in contact and objective body.
The contact and objective body in Contact Algorithm acquisition threedimensional model are searched for by non-binary trees；
According to formula:Obtain the normal direction contact between contact and objective body
Power；Wherein, βcAnd βtRespectively indicate contact and objective body, m and n are respectively indicated for the limited of discrete touch body and objective body
Number of unit,Indicate contact gesture, S indicates the boundary of mutual embedded part, and n indicates the exterior normal direction vector on the boundary.
According to formula:Obtain the tangential contact force between contact and objective body；Wherein ktIt is tangential
Rigidity, ηtFor tangential viscous damping coefficient, dtAnd vtRespectively indicate circumferentially opposite displacement and relative velocity.
If judgementThen according to Coulomb friction law, again according to formulaDescribed in calculating
Tangential contact force.
According to the vector sum of normal direction contact force and tangential contact force as the contact force.
On the basis of the above embodiments, the finite element discretization of the power balance equation of heap stone original state threedimensional model
Form are as follows:
X=X+u；
Wherein, M, C respectively indicate mass matrix and damping matrix, fintIndicate force vector in caused by deformation, fcExpression connects
Touch, fextIndicate all outer force vectors in addition to contact force,WithRespectively indicate the acceleration and speed of finite element node, X
For the initial position vector of finite element node, u is the motion vector of finite element node.
As a preferred embodiment, the calculating side of the rock-fill concrete casting process numerical simulation in the embodiment of the present invention
Method specifically includes:
Heap stone original state threedimensional model is established, the heap stone includes several block stone bodies.
According to finite element/discrete element coupling analytical method FEM/DEM, analysis heap stone original state threedimensional model is calculated
Banking process and final accumulation form.
The final accumulation form of heap stone original state threedimensional model is inputted in the form of boundary, and is established from close
Real concrete model, the self-compacting concrete model are Bingham rheological model.
Based on smoothed particle method SPH method, analysis self-compacting concrete is calculated in the final accumulation of heap stone
Casting process and final occupied state in form.
The embodiment of the present invention uses Paraview software as the poster processing soft, handles calculated result.It needs
It is bright, inherently a kind of mature software of numerical result processing of Paraview software, mainly analysis enrockment coagulation
The visualization processing and velocity field of native casting process, the analysis of stress field and the display for finally pouring state memorization compactness.
Heap stone original state threedimensional model is established, the heap stone includes several block stone bodies, specifically includes: using three
Modeling software is tieed up according to the predetermined parameter of heap stone, establishes heap stone original state threedimensional model；Alternatively, to described piece
Stone body appearance profile carries out 3D and scans to form input data and be input in 3 d modeling software, and it is three-dimensional to establish heap stone original state
Model.
According to finite element/discrete element coupling analytical method FEM/DEM, analysis heap stone original state threedimensional model is calculated
Banking process and final accumulation form, specifically include: it is three-dimensional to obtain the heap stone original state according to limited element analysis technique FEM
Each piece of stone body Continuum Mechanics behavior in model；The heap stone original state three is obtained according to DEM analysis method DEM
Non-continuum mechanics behavior in dimension module between block stone body；By the Continuum Mechanics behavior of all pieces of stone bodies and discontinuous Jie
Matter mechanical behavior substitutes into the power balance equation of heap stone original state threedimensional model, obtains the banking process and final accumulation
Form.
Based on smoothed particle method SPH method, analysis self-compacting concrete is calculated in the final accumulation of heap stone
Casting process and final occupied state in form, specifically include: SPH method are based on, one by one described in the update of time step iteration
Density, the position and speed of fluid particles in self-compacting concrete model；Density, position and speed based on the fluid particles
Degree analyzes casting process and final occupied state of the self-compacting concrete in the final accumulation form of heap stone.
Based on the above embodiment, Fig. 2 is the computing system embodiment of rock-fill concrete casting process numerical simulation of the present invention
Module map, as shown in Figure 2, comprising: accumulation analysis module 201, for according to finite element/discrete element coupling analytical method FEM/
DEM calculates the banking process for analyzing heap stone original state threedimensional model and final accumulation form；Analysis module 202 is poured, is used
In after calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation form, it is based on smooth particle
Fluid dynamics SPH method calculates casting process of the analysis self-compacting concrete in the final accumulation form of heap stone and most
Whole occupied state.
The short message dissemination system of the embodiment of the present invention can be used for executing rock-fill concrete casting process numerical value shown in FIG. 1
The technical solution of the calculation method embodiment of simulation, it is similar that the realization principle and technical effect are similar, and details are not described herein again.
Based on the above embodiment, Fig. 3 is the calculating of the rock-fill concrete casting process numerical simulation in the embodiment of the present invention
The block schematic illustration of equipment.Referring to FIG. 3, the embodiment of the present invention provides a kind of meter of rock-fill concrete casting process numerical simulation
Calculate equipment, comprising: processor (processor) 310,320, memory communication interface (Communications Interface)
(memory) 330 and bus 340, wherein processor 310, communication interface 320, memory 330 are completed mutually by bus 340
Between communication.Processor 310 can call the logical order in memory 330, to execute following method, comprising: according to limited
Member/discrete element coupling analytical method FEM/DEM calculates the banking process for analyzing heap stone original state threedimensional model and final heap
Product form；After calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation form, it is based on light
Sliding particle fluid dynamics SPH method calculates analysis self-compacting concrete pouring in the final accumulation form of heap stone
Journey and final occupied state.
The embodiment of the present invention discloses a kind of computer program product, and the computer program product is non-transient including being stored in
Computer program on computer readable storage medium, the computer program include program instruction, when described program instructs quilt
When computer executes, computer is able to carry out calculation method provided by above-mentioned each method embodiment, for example, according to limited
Member/discrete element coupling analytical method FEM/DEM calculates the banking process for analyzing heap stone original state threedimensional model and final heap
Product form；After calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation form, it is based on light
Sliding particle fluid dynamics SPH method calculates analysis self-compacting concrete pouring in the final accumulation form of heap stone
Journey and final occupied state.
Based on the above embodiment, the embodiment of the present invention provides a kind of non-transient computer readable storage medium, described non-temporary
State computer-readable recording medium storage computer instruction, it is real that the computer instruction makes the computer execute above-mentioned each method
Apply calculation method provided by example, for example, according to finite element/discrete element coupling analytical method FEM/DEM, calculate analysis heap
The banking process of stone original state threedimensional model and final accumulation form；The heap stone original state three-dimensional is analyzed calculating
After the banking process of model and final accumulation form, it is based on smoothed particle method SPH method, it is mixed to calculate analysis self-compaction
Casting process and final occupied state of the solidifying soil in the final accumulation form of heap stone.
Those of ordinary skill in the art will appreciate that: realize that above equipment embodiment or embodiment of the method are only schematic
, wherein the processor and the memory can be physically separate component may not be it is physically separated, i.e.,
It can be located in one place, or may be distributed over multiple network units.It can select according to the actual needs therein
Some or all of the modules achieves the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creative labor
In the case where dynamic, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as USB flash disk, mobile hard disk, ROM/RAM, magnetic disk, CD
Deng, including some instructions use is so that a computer equipment (can be personal computer, server or the network equipment etc.)
Execute method described in certain parts of each embodiment or embodiment.
The calculation method and system of rock-fill concrete casting process numerical simulation provided in an embodiment of the present invention, pass through setting
Finite element and discrete element coupling analytical method carry out mechanical analysis, and it is progressive under stress can completely to depict heap stone
Banking process, calculated result are more accurate.Self-compacting concrete is analyzed in the final accumulation shape of heap stone by setting SPH method
Casting process and final occupied state in state, are capable of the casting process of accurate simulation self-compacting concrete, not only can be intuitive
Flow regime of the self-compacting concrete in enrockment hole is observed on ground, can also understand the final closely knit shape of rock-fill concrete
State.Present invention saves a large amount of human cost, experimental cost and economic costs, and avoid experimental method and be difficult to see online
The problem of survey.By be arranged one by one time step iteration update the fluid particles in the self-compacting concrete model density,
Position and speed can accurately calculate the casting process of rock-fill concrete.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. a kind of calculation method of rock-fill concrete casting process numerical simulation characterized by comprising
According to finite element/discrete element coupling analytical method FEM/DEM, the accumulation of analysis heap stone original state threedimensional model is calculated
Process and final accumulation form；
After calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation form, it is based on smooth grain
Subflow body dynamics SPH method, calculate casting process of the analysis self-compacting concrete in the final accumulation form of heap stone and
Final occupied state.
2. calculation method according to claim 1, which is characterized in that described according to finite element/discrete element coupling analysis side
Method FEM/DEM calculates the banking process for analyzing heap stone original state threedimensional model and final accumulation form, before further include:
Heap stone original state threedimensional model is established, the heap stone includes several block stone bodies.
3. calculation method according to claim 1, which is characterized in that the heap stone original state three is analyzed in the calculating
The banking process of dimension module and final accumulation form and it is described be based on smoothed particle method SPH method, calculate analysis
Casting process and final occupied state of the self-compacting concrete in the final accumulation form of heap stone, between further include:
Self-compacting concrete model is established, the self-compacting concrete model is Bingham rheological model.
4. calculation method according to claim 3, which is characterized in that described to be based on the smoothed particle method side SPH
Method calculates casting process and final occupied state of the analysis self-compacting concrete in the final accumulation form of heap stone, specifically
Include:
Based on SPH method, time step iteration updates the density of the fluid particles in the self-compacting concrete model, position one by one
It sets and speed；
Density, position and speed based on the fluid particles, final accumulation form of the analysis self-compacting concrete in heap stone
In casting process and final occupied state.
5. calculation method according to claim 4, which is characterized in that the iteration of time step one by one updates described from close
Density, the position and speed of fluid particles in real concrete model, specifically include:
Time step one by one is calculated the variable density of the fluid particles by continuity equation, is captured from by surface particle simultaneously
The density correction for carrying out Free Surface particle, obtains the density field at each moment in the self-compacting concrete model；
Time step one by one is updated the track of each fluid particles movement by fluid particles displacement equation, obtains any fluid
The position at particle each moment；
Time step one by one determines the pressure of any fluid particles by state equation, the pressure gauge based on each fluid particles
The acceleration that external force generates is calculated, per a period of time in the self-compacting concrete model is calculated by the equation of momentum based on the acceleration
The velocity field at quarter.
6. calculation method according to claim 2, which is characterized in that it is described to establish heap stone original state threedimensional model,
It specifically includes:
Using 3 d modeling software according to the predetermined parameter of heap stone, heap stone original state threedimensional model is established；Or
Person,
3D is carried out to described piece of stone body appearance profile to scan to form input data and be input in 3 d modeling software, establishes heap stone
Original state threedimensional model.
7. calculation method according to claim 1, which is characterized in that described according to finite element/discrete element coupling analysis side
Method FEM/DEM calculates the banking process for analyzing heap stone original state threedimensional model and final accumulation form, specifically includes:
Each piece of stone body Continuum Mechanics in the heap stone original state threedimensional model is obtained according to limited element analysis technique FEM
Behavior；
The Discontinuous transmission power in the heap stone original state threedimensional model between block stone body is obtained according to DEM analysis method DEM
Scholarship and moral conduct is；
By the Continuum Mechanics behavior and non-continuum mechanics behavior substitution heap stone original state three-dimensional of all pieces of stone bodies
The power balance equation of model obtains the banking process and final accumulation form.
8. a kind of computing system of rock-fill concrete casting process numerical simulation characterized by comprising
Analysis module is accumulated, for it is initial to calculate analysis heap stone according to finite element/discrete element coupling analytical method FEM/DEM
The banking process of state threedimensional model and final accumulation form；
Analysis module is poured, for calculating the banking process for analyzing the heap stone original state threedimensional model and final accumulation
After form, it is based on smoothed particle method SPH method, calculates analysis self-compacting concrete in the final accumulation shape of heap stone
Casting process and final occupied state in state.
9. a kind of calculating equipment of rock-fill concrete casting process numerical simulation characterized by comprising
At least one processor；And
At least one processor being connect with the processor communication, in which:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to instruct energy
Enough calculation methods executed as described in claim 1 to 7 is any.
10. a kind of non-transient computer readable storage medium, which is characterized in that the non-transient computer readable storage medium is deposited
Computer instruction is stored up, the computer instruction makes the computer execute the calculation method as described in claim 1 to 7 is any.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810699310.5A CN109033537B (en) | 2018-06-29 | 2018-06-29 | Calculation method and system for numerical simulation in rock-fill concrete pouring process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810699310.5A CN109033537B (en) | 2018-06-29 | 2018-06-29 | Calculation method and system for numerical simulation in rock-fill concrete pouring process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109033537A true CN109033537A (en) | 2018-12-18 |
| CN109033537B CN109033537B (en) | 2020-12-01 |
Family
ID=65520981
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810699310.5A Expired - Fee Related CN109033537B (en) | 2018-06-29 | 2018-06-29 | Calculation method and system for numerical simulation in rock-fill concrete pouring process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109033537B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110750860A (en) * | 2019-09-11 | 2020-02-04 | 四川轻化工大学 | Soil slope landslide overall process analysis method |
| CN111191373A (en) * | 2020-01-03 | 2020-05-22 | 长沙有色冶金设计研究院有限公司 | High-concentration tailing filling slurry stirring process optimization method |
| CN111666699A (en) * | 2020-04-30 | 2020-09-15 | 山东大学 | Rock mass engineering cross-scale simulation calculation method based on REV full-area coverage |
| CN111914321A (en) * | 2020-06-09 | 2020-11-10 | 西安理工大学 | Method for establishing rock-fill concrete three-phase mesoscopic model |
| CN112052587A (en) * | 2020-09-02 | 2020-12-08 | 中国人民解放军陆军工程大学 | Three-dimensional microscopic discrete body model compacting method for sandy soil cushion |
| CN112364553A (en) * | 2021-01-13 | 2021-02-12 | 西南交通大学 | Method for evaluating coupling of finite element-discrete element of surface layer seepage corrosion of ballastless track foundation bed of high-speed rail |
| CN113158413A (en) * | 2021-02-05 | 2021-07-23 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | Diamond circular saw blade sawing system modeling and hard and brittle material sawing processing process numerical simulation method |
| CN113338218A (en) * | 2021-08-06 | 2021-09-03 | 西南交通大学 | Multi-scale multi-medium comprehensive inversion method for debris flow flexible protection |
| CN113515878A (en) * | 2021-07-07 | 2021-10-19 | 重庆交通大学 | Accumulation body discrete element three-dimensional modeling method based on shape and breakage of rock block |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN106960070A (en) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT |
-
2018
- 2018-06-29 CN CN201810699310.5A patent/CN109033537B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN106960070A (en) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT |
Non-Patent Citations (3)
| Title |
|---|
| 方坤河: "《碾压混凝土材料、结构与性能》", 28 February 2004, 武汉大学出版社 * |
| 李锡夔，万柯: "颗粒材料多尺度分析的连接尺度方法", 《力学学报》 * |
| 王永鹏: "光滑粒子流体动力学算法研究及其在流体力学中的数值模拟应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110750860A (en) * | 2019-09-11 | 2020-02-04 | 四川轻化工大学 | Soil slope landslide overall process analysis method |
| CN110750860B (en) * | 2019-09-11 | 2023-05-05 | 四川轻化工大学 | Whole-process analysis method for landslide |
| CN111191373A (en) * | 2020-01-03 | 2020-05-22 | 长沙有色冶金设计研究院有限公司 | High-concentration tailing filling slurry stirring process optimization method |
| CN111191373B (en) * | 2020-01-03 | 2023-04-07 | 长沙有色冶金设计研究院有限公司 | High-concentration tailing filling slurry stirring process optimization method |
| WO2021218070A1 (en) * | 2020-04-30 | 2021-11-04 | 山东大学 | Cross-scale simulation calculation method for rock mass engineering based on rev full area coverage |
| CN111666699A (en) * | 2020-04-30 | 2020-09-15 | 山东大学 | Rock mass engineering cross-scale simulation calculation method based on REV full-area coverage |
| CN111666699B (en) * | 2020-04-30 | 2023-06-02 | 山东大学 | Rock mass engineering trans-scale simulation calculation method based on REV full-area coverage |
| CN111914321A (en) * | 2020-06-09 | 2020-11-10 | 西安理工大学 | Method for establishing rock-fill concrete three-phase mesoscopic model |
| CN111914321B (en) * | 2020-06-09 | 2024-01-16 | 西安理工大学 | Method for establishing three-phase microscopic model of rock-fill concrete |
| CN112052587A (en) * | 2020-09-02 | 2020-12-08 | 中国人民解放军陆军工程大学 | Three-dimensional microscopic discrete body model compacting method for sandy soil cushion |
| CN112052587B (en) * | 2020-09-02 | 2023-12-01 | 中国人民解放军陆军工程大学 | Compaction method of three-dimensional microscopic discrete body model of sand cushion |
| CN112364553A (en) * | 2021-01-13 | 2021-02-12 | 西南交通大学 | Method for evaluating coupling of finite element-discrete element of surface layer seepage corrosion of ballastless track foundation bed of high-speed rail |
| CN113158413A (en) * | 2021-02-05 | 2021-07-23 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | Diamond circular saw blade sawing system modeling and hard and brittle material sawing processing process numerical simulation method |
| CN113515878A (en) * | 2021-07-07 | 2021-10-19 | 重庆交通大学 | Accumulation body discrete element three-dimensional modeling method based on shape and breakage of rock block |
| CN113515878B (en) * | 2021-07-07 | 2023-06-20 | 重庆交通大学 | Bulk discrete element three-dimensional modeling method based on block stone shape and breakage |
| CN113338218B (en) * | 2021-08-06 | 2021-10-26 | 西南交通大学 | Multi-scale multi-medium comprehensive inversion method for debris flow flexible protection |
| CN113338218A (en) * | 2021-08-06 | 2021-09-03 | 西南交通大学 | Multi-scale multi-medium comprehensive inversion method for debris flow flexible protection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109033537B (en) | 2020-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109033537A (en) | The calculation method and system of rock-fill concrete casting process numerical simulation | |
| Tsubota et al. | Simulation study on effects of hematocrit on blood flow properties using particle method | |
| Lv et al. | A novel coupled level set and volume of fluid method for sharp interface capturing on 3D tetrahedral grids | |
| CN108984829A (en) | The calculation method and system of rock-fill concrete rockfill banking process | |
| Ramalho et al. | The numerical simulation of crack propagation using radial point interpolation meshless methods | |
| Lii et al. | Ice melting simulation with water flow handling | |
| CN108717722A (en) | Fluid animation generation method and device based on deep learning and SPH frames | |
| CN105022928A (en) | Digitized real-time determination method for center-of-gravity position of fuel system of aircraft | |
| Biscarini et al. | Application of the lattice Boltzmann method for large‐scale hydraulic problems | |
| Tsubota | Short note on the bending models for a membrane in capsule mechanics: comparison between continuum and discrete models | |
| Darbani et al. | Meshless method for shallow water equations with free surface flow | |
| Kalateh et al. | Simulation of cavitating fluid–Structure interaction using SPH–FE method | |
| Andrade et al. | SPH fluids for viscous jet buckling | |
| Greaves | Simulation of interface and free surface flows in a viscous fluid using adapting quadtree grids | |
| KR20110072551A (en) | Method for analyzing shallow water flow using the two-dimensional river flow model with tensor-type eddy viscosity | |
| Ngo-Cong et al. | Incompressible smoothed particle hydrodynamics-moving IRBFN method for viscous flow problems | |
| Firoozjaee et al. | Element-free Galerkin method for numerical simulation of sediment transport equations on regular and irregular distribution of nodes | |
| Fujisawa et al. | Particle-based shallow water simulation with splashes and breaking waves | |
| Mirafzali et al. | Hydroelastic analysis of fully nonlinear water waves with floating elastic plate via multiple knot B-splines | |
| Ge et al. | Blending isogeometric and Lagrangian elements in three-dimensional analysis | |
| Zhao et al. | Comparison of depth-averaged concentration and bed load flux sediment transport models of dam-break flow | |
| dos Santos Brito et al. | Large viscoelastic fluid simulation on GPU | |
| Tsuruta | Improved particle method with high-resolution and computational stability for solid-liquid two-phase flows | |
| Schruff et al. | Numerical simulation of pore fluid flow and fine sediment infiltration into the riverbed | |
| Schuster | Large scale simulations of swirling and particle-laden flows using the Lattice-Boltzmann Method |
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 | ||
| 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: 20201201 Termination date: 20210629 |