共查询到20条相似文献，搜索用时 578 毫秒
1.
The aim of this paper is to obtain, from a numerical study using the discrete element method, quantitative information on the behaviour of a rockfill dam, particularly the effect of rock ageing on the movements of the dam. A numerical model is defined at three different scales: the block scale, the elementary representative volume scale and the dam scale. The model uses the discrete element code PCF2D, considering breakable clusters of 2D balls. Three rock-ageing laws, taking into account the time-evolving resistance for blocks, are proposed. The different parameters are determined from experimental data of laboratory tests performed on rock blocks. The relevance of the model seems acceptable, considering that the displacements obtained in the simulations are in the same order of magnitude as the displacements measured on the actual dam. 相似文献
2.
Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material. 相似文献
3.
颗粒破碎在岩土工程领域是很常见的现象,土工试验中无法显示颗粒破碎过程及其影响,本文采用离散单元软件PFC2D模拟了考虑颗粒破碎影响的粗粒土的直剪试验,给出了考虑颗粒破碎的粗粒土直剪试验的模拟方法,分析粗粒土的剪应力-剪切位移关系、剪胀和剪切强度等宏观力学行为,探讨基本颗粒间黏聚力、单颗粒孔隙率和粗粒土试样的孔隙率对剪切强度的影响。结果表明:颗粒破碎对剪切强度的破碎准则有影响,颗粒不破碎试样的剪切强度符合Mohr-Coulomb准则;颗粒破碎试样的强度包络线是幂函数关系。 相似文献
4.
The dam area of the SUOXI hydropower project shows high terrain undulation and complex geological conditions, containing 6 faults and 7 weak inter-beds. A geometric model developed to represent the geology and engineering structures should incorporate the geological realities and should allow suitable mesh generation to perform numerical stress analysis. This is an important precondition to perform rock mass stability analysis of a dam foundation based on a numerical stress analysis software such as FLAC3D. Using the modeling tools available in FLAC3D, it is difficult to construct a complex geological model even after performing a large amount of plotting and data analyses. The 3-D geological modeling technique suggested in this paper, named as Sealed Geological Modeling (SGM), is a powerful tool for constructing complex geological models for rock engineering projects that require numerical stress analysis. Applying this technique, first, the geological interfaces are constructed for the dam area of SUOXI hydropower project using various interpolation procedures including geostatistical techniques. Then a unitary wire frame is constructed and the interfaces are connected seamlessly. As the next step, a block tracing technique is used to build a geological model that consists of 130 seamlessly connected blocks. Finally, based on the Advancing Front Technique (AFT), each block is discretized into tetrahedrons and a mesh is generated including 57,661 nodes and 215,471 tetrahedrons which is suitable to perform numerical stress analysis using FLAC3D. 相似文献
5.
This paper describes the geomechanical back-analysis of the flexural failure of an overhanging limestone slab (volume = 68 m3) that collapsed on 26/01/1999 on the abandoned Cellina Valley state road SS 251 (western Friuli, NE Italy). The survey carried out on site examining the contact surfaces between the block and the slope has ascertained that the rock slab was connected to the stable rock mass by means of a single great rock bridge. The detailed reconstruction of the slab geometry and of the actual restraint conditions acting before the rockfall permitted an estimation of the stress state by using a 3D finite element code (Strand7). The finite element numerical model, referred to the restrained cross-section measured on site (Ares = 2.82 m2), gives a maximum value of σt = 5.19 MPa for the tensile stress, which is in good agreement with the estimated value for the characteristic strength mobilized at rupture by the intact rock material (T0 = 3.5–5.5 MPa). Considerably lower maximum tensile stress (σt = 2.25 MPa), smaller than one half, is obtained if the whole contact surface is assumed as resisting surface (Ares = 10.50 m2) and the rock bridge is neglected. 相似文献
6.
Three‐dimensional numerical analysis of concrete‐faced rockfill dam using dual‐mortar finite element method with mixed tangential contact constraints 
下载免费PDF全文
下载免费PDF全文 Concrete‐faced rockfill dam (CFRD) is a popular alternative to traditional dam types in the last two decades. The modelling of CFRD involves complex multi‐body contact and strong geometry and material nonlinearities. We present a numerical approach for the modelling of CFRDs in this paper. Based on the dual‐mortar finite element method, the presented approach considers different parts of rockfill and all concrete slabs as independent deformable continuum. The multi‐body contacts are modelled using Lagrange multipliers with a weak form segment‐to‐segment contact strategy. To alleviate instability induced by strong geometry nonlinearity in the slab–slab contact, we propose a mixed type of constraints for the tangential contact. A general transformation scheme is introduced to simplify the implementation of contact constraints. Three‐dimensional analysis of Tianshengqiao‐1 CFRD is performed. The nonlinear and time‐dependent deformation of the rockfill is considered. We study the influence of the rockfill deformation on the reliability of the concrete face. Three major concerns of the face, that is, the axial compression, the slab–slab separation and the face‐rockfill separation, are discussed in detail. The numerical results are compared with data from in‐situ observation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
7.
8.
《Geomechanics and Geoengineering》2013,8(3):175-181
Railway ballast particles undergo significant amount of breakage under repeated train load. Breakage of ballast particles, especially highly angular fresh ones, causes an increase in settlement, contributing to track degradation. The quantitative analysis of the influence of breakage on the stress-strain properties of ballast can be performed either experimentally or numerically. Numerical modeling has the advantage of simulating ballast breakage subject to various types of loading and different boundary conditions for a range of material properties. In this paper, ballast breakage under cyclic loading is simulated using a 2D discrete element method (DEM) utilizing the software PFC2D . A new subroutine is developed and incorporated in the PFC2D analysis to study ballast breakage and to quantify breakage in relation to particle size distribution. The influence of confining pressure on both breakage and permanent deformation is also studied and compared with laboratory observations. The findings of this paper provide an insight into the true ballast behavior under cyclic loading and are expected to assist railway practitioners in developing suitable design criteria for track stability. 相似文献
9.
A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
10.
数值流形方法是包含流形元、有限元及DDA在内的数值方法体系,建立流形元与DDA块体的接触方程,则可实现流形方法框架下的连续介质和散体系统共同作用模拟。针对填石路堤工程,编制了大型数值计算程序,采用块体随机生成、块体粒径控制及块体自然堆积的方法建立散体系统的DDA模型,对路堤的分层铺设、碾压及工后沉降变形等进行模拟分析。通过算例表明,在数值流形方法框架下,采用流形元与DDA共同作用的方法,可以很好地对同时存在连续变形和散体大变形的体系进行计算分析,其对该类问题的模拟更接近分析对象的实际情况,有助于从根本上揭示分析对象变形的细观机制和规律,并能考察更多因素对工程问题的影响。 相似文献
11.
This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
12.
Laboratory uniaxial tests on cylindrical specimens of a simulated rock are used to identify the influence of the diameter (D) and height (L) of the specimens on the damage parameters. The relations are expressed in terms of the characteristic parameter C = L2/D. 相似文献
13.
A combined finite-discrete element approach is used to simulate the complete 3D fracture process during conventional laboratory testing, including Brazilian indirect tension and uniaxial and biaxial compression. A typical granite rock type (based on the Lac du Bonnet granite) was simulated to investigate the fracture pattern and mechanical strength of brittle rock in the laboratory. Damage intensity parameters (D21 and D32) are introduced and utilized to characterize the induced damage in the models. These parameters provide an improved representation of the cumulative associated damage and facilitate a quantitative characterization of crack intensity during testing. The numerical simulations included both 3D and 2D models, and show that there is a good agreement between the strength response derived from simulations both in 3D and 2D and the considered rock material. A good correlation also exists between the fracture pattern in 3D and the equivalent 2D models. The influence of confinement on the biaxial strength and the associated damage in compression is investigated. While axial splitting is the dominant failure mode at low confinement, finite-discrete element simulations show that a shear failure mode tends to dominate as the confinement increases. The dependency of dilation upon the confining pressure is also demonstrated, the dilation angle decreasing with increased confinement. 相似文献
14.
H. Y. Liu S. Q. Kou P.-A. Lindqvist C. A. Tang 《Rock Mechanics and Rock Engineering》2007,40(2):107-144
Summary A series of numerical tests including both rock mechanics and fracture mechanics tests are conducted by the rock and tool
(R–T2D) interaction code coupled with a heterogeneous masterial model to obtain the physical–mechanical properties and fracture
toughness, as well as to simulate the crack initiation and propagation, and the fracture progressive process. The simulated
results not only predict relatively accurate physical–mechanical parameters and fracture toughness, but also visually reproduce
the fracture progressive process compared with the experimental and theoretical results. The detailed stress distribution
and redistribution, crack nucleation and initiation, stable and unstable crack propagation, interaction and coalescence, and
corresponding load–displacement curves can be proposed as benchmarks for experimental study and theoretical research on crack
propagation. It is concluded that the heterogeneous material model is reasonable and the R–T2D code is stable, repeatable and a valuable numerical tool for research on the rock fracture process. 相似文献
15.
Timothy P. Loomis Jibamitra Ganguly Stephen C. Elphick 《Contributions to Mineralogy and Petrology》1985,90(1):45-51
Data from experimentally-induced diffusion profiles at approximately 40 Kbar, 1,300–1,500° C in spessartine-almandine couples and a pyrope-almandine couple at 40 Kbar, 1,440° C, described in Part I, were used to derive tracer diffusion coefficients (D
*) of Fe, Mn and Mg in garnet. The experimental data were fitted by numerical simulations that model multicomponent, compositionally-dependent difussion, including the effects of nonideal thermodynamic mixing. The simulations use the formalism of irreversible thermodynamics and an eigenvector technique of solution. We were able to fit the asymmetrical spessartine-almandine profiles using constant D
* and either the Darken/Hartley-Crank or Manning-Lasaga models relating D
* and interdiffusion coefficients, and both models yielded D
Mg
*
consistent with the direct measurement of D
Mg
*
in by Cygan and Lasaga (1985) at lower temperatures (750–900° C). The results (equations 4.1–4.3 and Table 1) indicate that D
Fe
*
D
Mg
*
<D
Mn
*
and Q
FeQ
Mg>Q
Mn, where Q is the activation energy. In contrast, the asymmetry of pyrope-almandine profiles is too great to fit with either tracer model assuming constant D
* and indicates that D
Mg
*
is similar to its value in spessartine-almandine couples but D
Fe
*
is an order of magnitude less. The fit also suggests that D
Ca
*
< D
Fe
*
Mg
*
in pyrope-almandine couples. Synthesis of data from the two types of diffusion couples suggests that D
Mg
*
is insensitive to compositional changes, whereas D
Fe
*
is affected by Mn/Mg and Fe/Mg ratios and probably by other factors. These compositional effects on tracer coefficients are compatible with those documented by Morioka (1983) for cation diffusion in olivine. 相似文献
16.
Lorenzo Fedele Alberto Zanetti Vincenzo Morra Michele Lustrino Leone Melluso Riccardo Vannucci 《Contributions to Mineralogy and Petrology》2009,158(3):337-356
Trace element partition coefficients between clinopyroxenes and associated glassy matrix (Cpx/L
D) have been determined for 13 REE, HFSE4+,5+, U, Th, Sr, Pb, Sc and V from combined LA-ICP-MS/EMP analyses in selected trachytes and trachyphonolites from Campi Flegrei.
Composition of clinopyroxene and glass is pretty homogeneous in the trachyphonolites, pointing to an overall attainment of
the equilibrium conditions. In trachytes, conversely, phases show some compositional heterogeneity (due to the presence of
clinopyroxene xenocrysts) that requested a more careful petrographic and geochemical inspection of the samples to assess the
equilibrium clinopyroxene composition. In the trachyte clinopyroxenes, REE are compatible from Nd to Lu (Cpx/L
D up to 2.9), like Y, Ti, Sc and V. The Cpx/L
D for Eu is lower than those of the adjacent REE, highlighting Eu2+ contribution. High D values are also shown by U, Th, Pb, Zr, Hf, Nb and Ta relatively to basaltic and andesitic systems, whereas the D
Sr is roughly similar to that found for less evolved magmas. Trachyphonolites are characterized by an overall decrease of the
Cpx/L
D for highly-charged cations (with the exception of V), and by a slight increase of D
Sr. REE are still compatible from Nd to Lu (Cpx/L
D up to 2.1), like Ti, Y, Sc and V. This variation is also predicted for REE and Y by models based on the elastic strain theory,
being consistent with the slightly lower polymerization degree estimated for the trachyphonolites. However, the observed Cpx/L
D
(REE,Y) are matched by the modelled ones only considering very low T (≤825°C), which are believed unlikely. This mismatch cannot be attributed to effects induced by the water-rich composition
of the trachyte–trachyphonolite suite, since they would lower the observed Cpx/L
D
(REE,Y). Moreover, the anomalous inflections of measured Cpx/L
D for HREE suggests some crystal-chemical control, such as the entrance of these elements in a site distinct from M2. It is
concluded that the large Cpx/L
D determined for trachytes and trachyphonolites are likely induced by hitherto unconstrained changes of the Z3+ activities related to the composition of melt and/or solid. All these considerations strongly highlight the importance of
a direct characterization of trace element partitioning in natural samples from magmatic systems poorly characterized by experimental
studies.
Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. 相似文献
17.
Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density Dr=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
18.
Summary Discontinuities or joints in the rock mass have various shapes and sizes. Along with the joint orientation and spacing, the
joint persistence, or the relative size of the joint, is one of the most important factors in determining the block sizes
of jointed rock masses. Although the importance of joint persistence on the overall rock mass strength has long been identified,
the impact of persistence on rock strength is in most current rock mass classification systems underrepresented. If joints
are assumed to be persistent, as is the case in most designs, the sizes of the rock blocks tend to be underestimated. This
can lead to more removable blocks than actually exist in-situ. In addition, a poor understanding of the rock bridge strength
may lead to lower rock mass strengths, and consequently, to excessive expenditure on rock support.
In this study, we suggest and verify a method for the determination of the block sizes considering joint persistence. The
idea emerges from a quantitative approach to apply the GSI system for rock mass classification, in which the accurate block
size is required. There is a need to statistically analyze how the distribution of rock bridges according to the combination
of joint orientation, spacing, and persistence will affect the actual size of each individual block. For this purpose, we
generate various combinations of joints with different geometric conditions by the orthogonal arrays using the distinct element
analysis tools of UDEC and 3DEC. Equivalent block sizes (areas in 2D and volumes in 3D) and their distributions are obtained
from the numerical simulation. Correlation analysis is then performed to relate the block sizes predicted by the empirical
equation to those obtained from the numerical model simulation. The results support the concept of equivalent block size proposed
by Cai et al. (2004, Int. J. Rock Mech. Min. Sci., 41(1), 3–19). 相似文献
19.
N. P. Honkanadavar Nripendra Kumar Murari Ratnam 《Geotechnical and Geological Engineering》2014,32(4):1001-1015
Two types of modeled rockfill materials were collected from Renuka dam site, Himachal Pradesh, India and Salma dam site, Afghanistan. The rockfill material collected from Renuka dam site is rounded to sub-rounded in shape and the rockfill material collected from Salma dam site is angular to sub-angular in shape. The prototype gradation rockfill material consists maximum particle size larger than 1,000 mm. Therefore, for carrying out laboratory testing and modeling the bahaviour, the prototype rockfill material is scaled down to the maximum particle size (dmax) of 25, 50 and 80 mm for both projects material using parallel gradation technique. Triaxial compression and Index properties tests were conducted on both project rockfill materials and are presented. From the triaxial behaviour, it is observed that the stress–strain behaviour is non-linear, inelastic and stress dependent for both the materials. The material compresses during the initial shearing and shows dilation effect with further shearing. It is observed that the ?-value for alluvial rockfill material increases with increase in dmax and reverse trend is observed for blasted quarried rockfill material which shows the importance of the type of material. The stress–strain-volume change behaviour of both projects modeled rockfill material was predicted by using hierarchical single surface (HISS) model based on elasto plasticity and compared with the laboratory test results. From the comparison, it is observed that both results match closely. It is, therefore, suggested that the behaviour of both types of rockfill materials can be characterized successfully using HISS model. 相似文献
20.
Determination of groundwater solute transport parameters in finite element modelling using tracer injection and withdrawal testing data 总被引：1，自引：0，他引：1 下载免费PDF全文
下载免费PDF全文 Nguyen Van Hoang 《地下水科学与工程》2021,9(4):292-303
The groundwater tracer injection and withdrawal tests are often carried out for the determination of aquifer solute transport parameters. However, the parameter analyses encounter a great difficulty due to the radial flow nature and the variability of the temporal boundary conditions. An adaptive methodology for the determination of groundwater solute transport parameters using tracer injection and withdrawal test data had been developed and illustrated through an actual case. The methodology includes the treatment of the tracer boundary condition at the tracer injection well, the normalization of tracer concentration, the groundwater solute transport finite element modelling and the method of least squares to optimize the parameters. An application of this methodology was carried out in a field test in the South of Hanoi city. The tested aquifer is Pleistocene aquifer, which is a main aquifer and has been providing domestic water supply to the city since the French time. Effective porosity of 0.31, longitudinal dispersivity of 2.2 m, and hydrodynamic dispersion coefficients from D = 220 m2/d right outside the pumping well screen to D =15.8 m2/d right outside the tracer injection well screen have been obtained for the aquifer at the test site. The minimal sum of squares of the differences between the observed and model normalized tracer concentration is 0.00119, which is corresponding to the average absolute difference between observed and model normalized concentrations of 0.035 5 (while 1 is the worst and 0 is the best fit). 相似文献