CN105956314A - Numerical test method capable of creating different earth-rock mixing ratios - Google Patents
Numerical test method capable of creating different earth-rock mixing ratios Download PDFInfo
- Publication number
- CN105956314A CN105956314A CN201610323741.2A CN201610323741A CN105956314A CN 105956314 A CN105956314 A CN 105956314A CN 201610323741 A CN201610323741 A CN 201610323741A CN 105956314 A CN105956314 A CN 105956314A
- Authority
- CN
- China
- Prior art keywords
- model
- unit
- node
- flac3d
- numerical
- 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.)
- Pending
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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a numerical test method capable of creating different earth-rock mixing ratios. The method comprises the following steps: an Flac3d model with N*N elements is created; after the earth-rock mixing ratios are determined, the initial model elements are divided into earth masses and rocks, and are written into an Flac3d calculation file; numerical tests are carried out through Flac3d software, and a test curve is simulated through statistics of axial loads and axial displacements. The method provided by the invention has the advantages that model creation is achieved; problems in real test parameter researches, such as high sampling difficulty, sample disturbance, and large discreteness of results, are solved; high modeling convenience and rapid convergence are achieve through the finite difference method. Therefore, the method can effectively serve engineering numerical simulation, and achieves high practicality.
Description
Technical field
The invention belongs to civil engineering numerical analysis field, be specifically related to a kind of number generating different soil stone mixing ratio
Value test method.
Background technology
Accumulation body belongs to a kind of discontinuous, uneven binary medium, and it is between the soil body and rock mass
Kind of special geobody, current people for its research also in exploring among.Accumulation body mechanical property not only takes
Certainly in the physico mechanical characteristic of soil with stone, more it is controlled by its inner structural features, such as soil stone mixing ratio, stone
The factors such as grating, stone shape, stone tendency all have considerable influence to accumulation body mechanical property.And laboratory is true
Determine that accumulation body physical and mechanical parameter exists sampling disturbance, each sample soil stone composite character is inconsistent, experimental error is big
Etc. shortcoming, causing result of the test discreteness relatively big, regularity is the strongest.Prior, laboratory test block is by yardstick
Limit, it is impossible to the well reflection soil stone composite character impact on accumulation body physico-mechanical properties.
Relatively rapider by numerical Method Research accumulation body parameter development in recent years, such method can well overcome
Sample difficulty, sample disturbance in actual experimental parameter study, result discreteness is big, test yardstick compared to accumulation
The problems such as internal portion structure dimension is less than normal.There is two large problems in numerical Method Research accumulation body parameter: (1) model
Set up problem, i.e. according to accumulation body every internal structure parameter (such as: soil stone mixing ratio, stone grating, stone
Shape, stone tendency etc.) generate the accumulation body test block meeting numerical experimentation requirement；(2) which kind of numerical computations is used
Means simulate real accumulation body parameter experiment.
Summary of the invention
The present invention is directed to the deficiencies in the prior art, it is provided that a kind of numerical experimentation side generating different soil stone mixing ratio
Method, modeling is convenient, fast convergence rate.
The present invention is by the following technical solutions:
A kind of Numerical Experimental Method generating different soil stone mixing ratio, comprises the following steps:
The first step: signal generating unit number is the Flac3d model of N*N: Definition Model unit number is by N*N side
Lattice becomes, and the size of model is 1*1*1/N (m), carries out the node that unit node coordinate is corresponding with each unit
Number and compose initial value, in the node coordinate write file Fracture.flac3d that assignment is completed；
Second step: initial model unit is divided into and is made up of the soil body and stone, and successively by node and correspondence thereof
Coordinate, unit and correspondence thereof node numbering write Flac3d calculation document in, generate different soil stone mixing
The numerical model of ratio；
3rd step: the Flac3d file that build is imported Flac3d software, is grouped, and to packet after
Model compose parameter, boundary condition is set, top, y direction and the bottom of model is retrained, then applies
Confined pressure, is loaded by displacement mode, adds up each and calculates xial feed and axial displacement, simulation test in step
Curve.
As preferably, model unit and node are numbered and comprise the following steps by the first step: element number be from
The bottom starts to increase the most successively, then the unit number of jth row i row is N*j+i, and each unit is corresponding
8 nodes, the coordinate of its correspondence is N* (j-1)+i, N* (j-1)+i+1；N* (j-1)+i+N, N* (j-1)+i+N+1,
N*N+N* (j-1)+i, N*N+N* (j-1)+i+1, N*N+N* (j-1)+i+N, N*N+N* (j-1)+i+N+1.
As preferably, the first step carries out assignment to unit node coordinate on model front view and comprises the following steps: set
Putting node x, the y of jth row i row, z coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively,
A [N*j+i] [2]=0.0, carries out assignment to unit node coordinate on the model back side, and its rule trend is jth row
Node x, the y of i row, z coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively,
A [N*j+i] [2]=1.0/ (N-1).
As preferably, second step is that random choose goes out sum unit on Flac3d model, and its attribute is fixed
Justice is the soil body, changes soil stone mixing ratio by the value changing sum, generates different numerical models.
As preferably, the 3rd step uses mohr-coulomb constitutive model to carry out numerical computations.
Beneficial effect: the present invention is determined by soil stone mixing ratio, quickly generates different soil stone mixing ratio model, and leads to
Cross Flac3d and carry out numerical experimentation simulation, be possible not only to solve model and set up problem, overcome actual experimental parameter
Research samples the problems such as difficulty, sample disturbance, result discreteness are big, employing finite difference calculus modeling side simultaneously
Just, fast convergence rate.
Accompanying drawing explanation
Fig. 1 is the design flow diagram of the present invention；
Fig. 2 is Flac3d hexahedro block grid；
Fig. 3 be soil-rock mixture rock-soil ratio be 50% illustraton of model；
Fig. 4 be soil-rock mixture rock-soil ratio be 60% illustraton of model；
Fig. 5 be soil-rock mixture rock-soil ratio be 90% illustraton of model；
Fig. 6 is different soil stone mixing ratio axial stress and axial strain graph of a relation.
Detailed description of the invention
The present invention is described in detail with detailed description of the invention below in conjunction with the accompanying drawings:
A kind of Numerical Experimental Method generating different soil stone mixing ratio, as it is shown in figure 1, comprise the following steps:
The first step: signal generating unit number is the Flac3d model of N*N, the size of model is 1*1*1/N (m).The most such as
Under:
Model unit and node are numbered by a, and element number rule is to start the most successively from the bottom
Increase, then the unit number of jth row i row is N*j+i, and corresponding 8 nodes of each unit, the coordinate of its correspondence is
N* (j-1)+i, N* (j-1)+i+1；N* (j-1)+i+N, N* (j-1)+i+N+1, N*N+N* (j-1)+i,
N*N+N* (j-1)+i+1, N*N+N* (j-1)+i+N, N*N+N* (j-1)+i+N+1；
B carries out assignment to unit node coordinate on model front view, and its rule trend is the node of jth row i row
X, y, z coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively, a [N*j+i] [2]=0.0,
Unit node coordinate on the model back side is carried out assignment, and its rule trend is the node x of jth row i row, y, z
Coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively, a [N*j+i] [2]=1.0/ (N-1).
Above-mentioned steps uses C language to program, as in figure 2 it is shown, Flac3d model node has the row of uniqueness
Column format, composes initial value to the node that unit node coordinate is corresponding with each unit, uses array
A [2*N*N] [3]={ [8]={ 0} represents 0} and node [N*N].
First unit node coordinate on model front view being carried out assignment, its programming is as follows:
Secondly carrying out model back-side unit node coordinate calculating assignment, its programming is as follows:
Finally, in the node coordinate write file Fracture.flac3d that assignment is completed, due to Flac3d
There is the data form of oneself, programmed as follows:
Second step: initial model unit is divided into and is made up of the soil body and stone, and it is written into Flac3d meter
Calculate in file, soil stone mixing ratio can be changed by the value changing sum, generate the numerical value of different soil stone mixing ratio
Model.Specific as follows:
First write node and the coordinate of correspondence thereof, then writing unit and the node numbering of correspondence thereof, the most laggard
Row packet, determines soil stone mixing ratio, uses random distribution to carry out being grouped (on the basis of the model that the first step generates
Random choose goes out sum unit, is the soil body by its attribute definition in follow-up process), for avoiding repeating choosing
Take unit, use recurrence thought again to be chosen by the unit of repetition.
Concrete programming is as follows:
3rd step: the Flac3d file built is imported Flac3d software, uses Flac3d software to carry out
Numerical experimentation, by statistics xial feed and axial shift simulation trial curve.Specific as follows:
First FLAC3d file is imported in model FLAC3d software, be grouped, and to the model after packet
Compose parameter.Choose and use mohr-coulomb constitutive model.
Then boundary condition is set, top, y direction and the bottom of model is retrained, then applies confined pressure,
Load by displacement mode.
Finally add up each and calculate xial feed and axial displacement in step, draw out stress～racking test curve.
Application example: said method is applied to soil-rock mixture model, and wherein the parameter of stone is: volume mould
Amount 25Gpa, modulus of shearing 12Gpa, density 2900kg/m3, cohesion is 4Gpa, angle of friction 45 °, anti-
Tensile strength 2Mpa, dilative angle 15 °；The parameter of the soil body is: bulk modulus 40Mpa, modulus of shearing 15Mpa,
Density 2200kg/m3, cohesion is 0.08Mpa, angle of friction 28 °, tensile strength 9Kpa, dilative angle 10
°。
Specifically comprise the following steps that
(1) different soil stone mixing ratio model is generated by C language programming as stated above.
(2) model generated is imported Flac3d software, use this model of mohr-coulomb, load
Speed is 5*10-5M/step loads.During calculating, statistics xial feed and axial displacement.Such as Fig. 3-5
Shown in, be respectively soil-rock mixture rock-soil ratio be the illustraton of model of 50%, 60% and 90%.
(3) result of step (2) is published picture display, as shown in Figure 6, display axial stress～axial strain
Relation.
Claims (5)
1. the Numerical Experimental Method generating different soil stone mixing ratio, it is characterised in that comprise the following steps:
The first step: signal generating unit number is the Flac3d model of N*N: Definition Model unit number is by N*N side
Lattice becomes, and the size of model is 1*1*1/N (m), carries out the node that unit node coordinate is corresponding with each unit
Number and compose initial value, in the node coordinate write file Fracture.flac3d that assignment is completed；
Second step: initial model unit is divided into and is made up of the soil body and stone, and successively by node and correspondence thereof
Coordinate, unit and correspondence thereof node numbering write Flac3d calculation document in, generate different soil stone mixing
The numerical model of ratio；
3rd step: the Flac3d file that build is imported Flac3d software, is grouped, and to packet after
Model compose parameter, boundary condition is set, top, y direction and the bottom of model is retrained, then applies
Confined pressure, is loaded by displacement mode, adds up each and calculates xial feed and axial displacement, simulation test in step
Curve.
The Numerical Experimental Method generating different soil stone mixing ratio the most according to claim 1, its feature exists
In, model unit and node are numbered and comprise the following steps by the first step: element number is from the beginning of the bottom
Increase the most successively, then the unit number of jth row i row is N*j+i, corresponding 8 nodes of each unit, its
Corresponding coordinate is N* (j-1)+i, N* (j-1)+i+1；N* (j-1)+i+N, N* (j-1)+i+N+1,
N*N+N* (j-1)+i, N*N+N* (j-1)+i+1, N*N+N* (j-1)+i+N, N*N+N* (j-1)+i+N+1.
The Numerical Experimental Method generating different soil stone mixing ratio the most according to claim 2, its feature exists
In, the first step carries out assignment to unit node coordinate on model front view and comprises the following steps: arrange jth row i
Node x, the y of row, z coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively,
A [N*j+i] [2]=0.0, carries out assignment to unit node coordinate on the model back side, and its rule trend is jth row
Node x, the y of i row, z coordinate is 1.0*i/ (N-1), a [N*j+i] [1]=1.0*j/ (N-1) respectively,
A [N*j+i] [2]=1.0/ (N-1).
The Numerical Experimental Method generating different soil stone mixing ratio the most according to claim 1, its feature exists
In, second step is that random choose goes out sum unit on Flac3d model, is the soil body by its attribute definition,
Change soil stone mixing ratio by the value changing sum, generate different numerical models.
The Numerical Experimental Method generating different soil stone mixing ratio the most according to claim 1, its feature exists
In, the 3rd step uses mohr-coulomb constitutive model to carry out numerical computations.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610323741.2A CN105956314A (en) | 2016-05-16 | 2016-05-16 | Numerical test method capable of creating different earth-rock mixing ratios |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610323741.2A CN105956314A (en) | 2016-05-16 | 2016-05-16 | Numerical test method capable of creating different earth-rock mixing ratios |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105956314A true CN105956314A (en) | 2016-09-21 |
Family
ID=56911663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610323741.2A Pending CN105956314A (en) | 2016-05-16 | 2016-05-16 | Numerical test method capable of creating different earth-rock mixing ratios |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105956314A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107084895A (en) * | 2017-04-27 | 2017-08-22 | 中国地质大学（武汉） | A kind of soil-rock mixture compaction test analogy method and system based on block stone form |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080319675A1 (en) * | 2007-06-22 | 2008-12-25 | Sayers Colin M | Method, system and apparatus for determining rock strength using sonic logging |
| CN102621009A (en) * | 2012-03-21 | 2012-08-01 | 武汉大学 | Test method for simulating long-term deformation of rockfill |
| CN104036538A (en) * | 2014-06-25 | 2014-09-10 | 清华大学 | Method and system for reconstruction and analysis of earth-rock aggregate three-dimensional microstructure |
-
2016
- 2016-05-16 CN CN201610323741.2A patent/CN105956314A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080319675A1 (en) * | 2007-06-22 | 2008-12-25 | Sayers Colin M | Method, system and apparatus for determining rock strength using sonic logging |
| CN102621009A (en) * | 2012-03-21 | 2012-08-01 | 武汉大学 | Test method for simulating long-term deformation of rockfill |
| CN104036538A (en) * | 2014-06-25 | 2014-09-10 | 清华大学 | Method and system for reconstruction and analysis of earth-rock aggregate three-dimensional microstructure |
Non-Patent Citations (2)
| Title |
|---|
| 廖秋林: "基于ANSYS平台复杂地质体FLAC3D模型的自动生成", 《岩石力学与工程学报》 * |
| 王琳: "土石混填地基隧道结构稳定性及补强措施研究", 《中国优秀硕士学位论文全文数据库-工程科技II辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107084895A (en) * | 2017-04-27 | 2017-08-22 | 中国地质大学（武汉） | A kind of soil-rock mixture compaction test analogy method and system based on block stone form |
| CN107084895B (en) * | 2017-04-27 | 2019-11-12 | 中国地质大学（武汉） | A kind of soil-rock mixture compaction test analogy method and system based on block stone form |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Farzam-Rad et al. | Isogeometric analysis of functionally graded plates using a new quasi-3D shear deformation theory based on physical neutral surface | |
| Liang et al. | A three-dimensional numerical investigation of the fracture of rock specimens containing a pre-existing surface flaw | |
| Shou et al. | 3D numerical simulation of initiation, propagation and coalescence of cracks using the extended non-ordinary state-based peridynamics | |
| CN113534291A (en) | Quantitative prediction method for different-scale fractures of low-permeability reservoir under constraint of rock mechanical layer | |
| CN110414183B (en) | Comprehensive pipe gallery anti-seismic analysis method with FRP (fiber reinforced plastic) ribs | |
| Pagani et al. | Free vibration analysis of composite plates by higher-order 1D dynamic stiffness elements and experiments | |
| Aslangil et al. | Numerical investigation of initial condition effects on Rayleigh-Taylor instability with acceleration reversals | |
| Wang et al. | Development and implementation of new triangular finite element based on MGE theory for bi-material analysis | |
| da Silva et al. | Analysis of the failure process by using the Lattice Discrete Element Method in the Abaqus environment | |
| Wang et al. | Improved maxwell model describing the whole creep process of salt rock and its programming | |
| CN106597542B (en) | A kind of prediction technique and device of reservoir characteristic parameter | |
| CN114547953A (en) | Fracturing construction parameter optimization method and system based on optimization design chart | |
| CN105443102A (en) | Horizontal well layer blast crack body modeling method | |
| Sha et al. | Modeling of Hydraulic Fracture of Concrete Gravity Dams by Stress‐Seepage‐Damage Coupling Model | |
| Wang et al. | Influence of material heterogeneity on the blast-induced crack initiation and propagation in brittle rock | |
| CN105956314A (en) | Numerical test method capable of creating different earth-rock mixing ratios | |
| Li et al. | A stepwise artificial boundary condition for wave propagation in elasto-plastic media | |
| Li et al. | Study on the influence of perforation parameters on hydraulic fracture initiation and propagation based on CDEM | |
| Fu et al. | Dual reciprocity boundary element based block model for discontinuous deformation analysis | |
| Liu et al. | Lognormal Distribution Function for Describing Seepage Damage Process of Single‐Cracked Rock | |
| Sołowski et al. | Material point method modelling of granular flow in inclined channels | |
| Lisjak et al. | 3D simulation of fluid-pressure-induced fracture nucleation and growth in rock samples | |
| Lu et al. | A new analysis of the complex two-dimensional multilayered anisotropic soil in time domain | |
| Yafeng et al. | Research on L-shaped concrete-filled steel tube core columns mechanical properties under eccentric | |
| Moresi et al. | An overview of numerical methods for Earth simulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160921 |
|
| RJ01 | Rejection of invention patent application after publication |