首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
This paper presents a numerical method for continuum modelling of the dynamic bit–rock interaction process in percussive drilling. The method includes a constitutive model based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. The interaction between the drill bit and rock is modelled using contact mechanics by treating the bit as a rigid body. As the bit–rock interaction in percussive drilling is a transient event, the method is implemented in explicit dynamics FEM. The rock strength heterogeneity is characterized at the mesoscopic level statistically using the Weibull distribution. The bit–rock interaction is simulated under axisymmetric conditions using cylindrical and hemispherical buttons. The choice of the quite complex constitutive model accounting, e.g. for plastic compaction, viscoplastic shear and tensile failure along with induced damage and rate dependency is justified by numerical simulations. Moreover, the quasi‐static and dynamic cases are compared in plane strain simulations. Finally, some results clarifying the discrepancy of opinions found in the literature concerning the side (lateral) crack formation are obtained. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

2.
The aim of this study is to numerically model the fracture system at percussive drilling. Because of the complex behavior of rock materials, a continuum approach is employed relying upon a plasticity model with yield surface locus as a quadratic function of the mean pressure in the principal stress space coupled with an anisotropic damage model. In particular, Bohus granite rock is investigated, and the material parameters are defined based on previous experiments. This includes different tests such as direct tension and compression, three‐point bending, and quasi‐oedometric tests to investigate the material behavior at both tension and confined compression stress states. The equation of motion is discretized using a finite element approach, and the explicit time integration method is employed. Edge‐on impact tests are performed, and the results are used to validate the numerical model. The percussive drilling problem is then modeled in 3D, and the bit‐rock interaction is considered using contact mechanics. The fracture mechanism in the rock and the bit penetration‐ resisting force response are realistically captured by the numerical model. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
This paper develops a practical approach to simulating fracture propagation in rock and concrete based on an augmented virtual internal bond (VIB) method in which the cohesion of solid is modeled as material particles interconnected by a network of randomized virtual micro bonds described by the Xu–Needleman potential. The micro bond potential is used to derive macroscale constitutive relations via the Cauchy–Born rule. By incorporating different energy contributions due to stretch and shearing, as well as different energy levels under tension and compression of each micro bond, the derived macro constitutive laws are particularly useful for modeling quasi‐brittle materials such as rock and concrete which usually have different Poisson ratios and much higher compressive strength than tensile strength. The mesh‐size sensitivity associated with strain‐softening in the present constitutive model is addressed by adjusting material constants near the crack tip so that the biJ‐integral is kept equal to the intrinsic fracture energy of the material. Numerical examples demonstrate that the proposed VIB method is capable of simulating mixed mode fracture propagation in rock and concrete with results in consistency with relevant experimental observations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
Strain gradient implies an important characteristic in localized damage deformation, which can be observed in the softening state of brittle materials, and strain gradients constitute the basic behaviours of localization failure area of the materials. The most important point in strain gradient is its damaging function including an internal length scale, which can be used to express the scale effects of mechanical responses of brittle rock mass. By extending the strain gradient theory and introducing an intrinsic material length scale into the constitutive law, the authors develop an isotropic damage model as well as a micro‐crack‐based anisotropic damage model for rock‐like materials in this paper. The proposed models were used to simulate the damage localization under uniaxial tension and plain strain compression, respectively. The simulated results well illustrated the potential of these models in dealing with the well‐known mesh‐sensitivity problem in FEM. In the computation, elements with C1 continuity have been implemented to incorporate the proposed models for failure localization. When regular rectangle elements are encountered, the coupling between finite difference method (FDM) and conventional finite element method (FEM) is used to avoid large modification to the existing FEM code, and to obtain relatively higher efficiency and reasonably good accuracy. Application of the anisotropic model to the 3D‐non‐linear FEM analysis of Ertan arch dam has been conducted and the results of its numerical simulation coincide well with those from the failure behaviours obtained by Ertan geophysical model test. In this paper, new applications of gradient theories and models for a feasible approach to simulate localized damage in brittle materials are presented. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

6.
The strength parameter mi in the Hoek–Brown strength criterion is empirical and was developed by trial and error. To better understand the fundamental relationship between mi and the physical characteristics of intact rock, this paper presents a systematic study of mi by representing intact rock as a densely packed cemented particle material and simulating its mechanical behavior using particle flow modeling. Specifically, the three‐dimensional particle flow code (PFC3D) was used to conduct numerical true triaxial compression tests on intact rock and to investigate the effect of non‐spherical micro‐particle parameters on mi. To generate numerical intact rock specimens containing non‐spherical micro‐particles, a new genesis process was proposed, and a specific loop algorithm was used based on the efficiency of the process and the acceptability of generated specimens. Four main parameters—number, aspect ratio, size, and shape—of non‐spherical micro‐particles were studied, and the results indicated that they all have great effect on mi. The strength parameter mi increases when the number, aspect ratio, or size is larger or the shape becomes more irregular, mainly as a result of the higher level of interlocking between particles. This confirms the observations from engineering experience and laboratory experiments. To simulate the right strength parameter mi, it is important to use appropriate non‐spherical micro‐particles by controlling these four parameters. This is further demonstrated by the simulation of two widely studied rocks, Lac du Bonnet granite and Carrara marble. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

7.
Recognition of non‐linear constitutive rock/soil model from experimental results is often multi‐modal in the large parameter space. A genetic evolution algorithm is thus proposed for its recognition, including that of structure of the model and coefficients in the model. The structure of the model can be firstly determined according to mechanical mechanism if the mechanism is clearly understood or searched by using evolutionary algorithm. The coefficients to be determined are then searched in global optional space. With the new evolutionary algorithm, the non‐linear stress–strain–time constitutive law to describe strain softening behaviours of diatomaceous soil under consolidated and undrained state was recognized by learning stress–strain–time behaviour of an intact sample under consolidated pressure of σc=0.1 MPa and strain velocity ofa=0.175%/min. This model gave reasonable prediction for diatomaceous soils under varying consolidated pressures (0.1–3.5 MPa) and strain velocities (0.0044–1.75%/min). It indicates that the methodology proposed in this paper is robust enough and strongly attractive for recognition of non‐linear constitutive model of soil and rock materials. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

10.
《Resource Geology》2018,68(1):93-100
CO2 sequestration into saline aquifers is considered to be one of the most promising options for reducing industrial CO2 emissions to the atmosphere. However, there are still many uncertainties regarding the storage of CO2 in the subsurface because of a lack of knowledge about CO2–water–rock interaction within CO2 reservoirs and the potential risk of CO2 leakage. In this study, we construct a semi‐open type experimental system that can reproduce the interactions under conditions close to those of actual CO2 reservoirs. Using the system, we conduct CO2–water–rock interaction experiments for 8 months to monitor the long‐term reaction and the mobilization of harmful metal elements. Altered tuffaceous rock is used in the experiment because these tuffaceous rock formations (called “Green Tuff”) are a potential candidate for CO2 storage in Japan. The results show that the major‐element water composition will converge to the point where host rock dissolution and secondary mineral precipitation are balanced; then, the interaction will proceed under a certain groundwater composition. In addition, we found that groundwater contamination by some metal elements (Ni, Ba, and Mn) may reach unsafe levels for drinking water as a result of CO2‐water–rock interaction.  相似文献   

11.
A three dimensional constitutive model is formulated for deformation analysis of jointed rock masses containing up to three joint sets with arbitrary spatial configurations. A representative elementary volume (REV) that represents the deformational response of the rock mass is defined and the constitutive relationships are developed based on superposition of deformations of the REV components. By representing the constitutive relationships in a tensorial form, the model is able to implement deformation anisotropy of jointed rock masses. The Mohr-Coulomb failure criterion with tension cut-off is used for the intact rock and the joint sets. The model is implemented in FLAC3D and the deformations and strength values calculated by the model are compared with the results from a 3DEC model and analytical solutions. The model results are in good agreement with those obtained from 3DEC.  相似文献   

12.
This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

13.
  相似文献   

14.
The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K0‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

15.
  相似文献   

16.
Granulite facies cordierite–garnet–biotite gneisses from the southeastern Reynolds Range, central Australia, contain both orthopyroxene‐bearing and orthopyroxene‐free quartzofeldspathic leucosomes. Mineral reaction microstructures at the interface of gneiss and leucosome observed in outcrop and petrographically, reflect melt‐rock interaction during crystallization. Accessory monazite, susceptible to fluid alteration, dissolution and recrystallization at high temperature, is tested for its applicability to constrain the chemical and P–T–time evolution of melt‐rock reactions during crystallization upon cooling. Bulk rock geochemistry and phase equilibria modelling constrain peak pressure and temperature conditions to 6.5–7.5 kbar and ~850°C, and U–Pb geochronology constrains the timing of monazite crystallization to 1.55 Ga, coeval with the Chewings Orogeny. Modelling predicts the presence of up to 15 vol.% melt at peak metamorphic conditions. Upon cooling below 800°C, melt extraction and in situ crystallization of melt decrease the melt volume to less than 7%, at which time it becomes entrapped and melt pockets induce replacement reactions in the adjacent host rock. Replacement reactions of garnet, orthopyroxene and K‐feldspar liberate Y, REE, Eu and U in addition to Mg, Fe, Al, Si and K. We demonstrate that distinguishing between monazite varieties solely on the basis of U–Pb ages cannot solve the chronological order of events in this study, nor does it tie monazite to the evolution of melt or stability of rock‐forming minerals. Rather, we argue that analyses of various internal monazite textures, their composition and overprinting relations allow us to identify the chronology of events following the metamorphic peak. We infer that retrograde reactions involving garnet, orthopyroxene and K‐feldspar can be attributed to melt‐rock interaction subsequent to partial melting, which is reflected in the development of compositionally distinct monazite textural domains. Internal monazite textures and their composition are consistent with dissolution and precipitation reactions induced by a high‐T melt. Monazite rims enriched in Y, HREE, Eu and U indicate an increased availability of these elements, consistent with the breakdown of orthopyroxene, garnet and K‐feldspar observed petrographically. Our study indicates that compositional and textural analysis of monazite in relation to major rock‐forming minerals can be used to infer the post‐peak chemical evolution of partial melts during high‐ to ultrahigh‐temperature metamorphism.  相似文献   

17.
The Sanbagawa metamorphic belt of southwest Japan is one of the type localities of subduction‐related high‐P metamorphism. However, variable pressure–temperature (P–T) paths and metabasic assemblages have been reported for eclogite units in the region, leading to uncertainty about the subduction zone paleo‐thermal structure and associated tectonometamorphic conditions. To analyse this variation, phase equilibria modelling was applied to the three main high‐P metabasic rock types documented in the region – glaucophane eclogite, barroisite eclogite and garnet blueschist – with modelling performed over a range of P, T, bulk rock H2O and bulk rock ferric iron conditions using thermocalc . All samples are calculated to share a common steep prograde P–T path to similar peak conditions of ～16–20 kbar and 560–610 °C. The results establish that regional assemblage variation is systematic, with the alternation in peak amphibole phase due to peak conditions overlapping the glaucophane–barroisite solvus, and bulk composition effects stabilizing blueschist v. eclogite facies assemblages at similar P–T conditions. Furthermore, the results reveal that a steep prograde P–T path is common to all eclogite units in the Sanbagawa belt, indicating that metamorphic conditions were consistent along strike. All localities are compatible with predictions made by a ridge approach model, which attributes eclogite facies metamorphism and exhumation of the Sanbagawa belt to the approach of a spreading ridge.  相似文献   

18.
In many rock engineering applications such as foundations, slopes and tunnels, the intact rock properties are not actually determined by laboratory tests, due to the requirements of high quality core samples and sophisticated test equipments. Thus, predicting the rock properties by using empirical equations has been an attractive research topic relating to rock engineering practice for many years. Soft computing techniques are now being used as alternative statistical tools. In this study, artificial neural network models were developed to predict the rock properties of the intact rock, by using sound level produced during rock drilling. A database of 832 datasets, including drill bit diameter, drill bit speed, penetration rate of the drill bit and equivalent sound level (Leq) produced during drilling for input parameters, and uniaxial compressive strength (UCS), Schmidt rebound number (SRN), dry density (ρ), P-wave velocity (Vp), tensile strength (TS), modulus of elasticity (E) and percentage porosity (n) of intact rock for output, was established. The constructed models were checked using various prediction performance indices. Goodness of the fit measures revealed that recommended ANN model fitted the data as accurately as experimental results, indicating the usefulness of artificial neural networks in predicting rock properties.  相似文献   

19.
Orogenic gold mineralization at the Damang deposit, Ghana, is associated with hydrothermal alteration haloes around gold‐bearing quartz veins, produced by the infiltration of a H2O–CO2–K2O–H2S fluid following regional metamorphism. Alteration assemblages are controlled by the protoliths with sedimentary rocks developing a typical assemblage of muscovite, ankerite and pyrite, while intrusive dolerite bodies contain biotite, ankerite and pyrrhotite, accompanied by the destruction of hornblende. Mineral equilibria modelling was undertaken with the computer program thermocalc , in subsets of the model system MnO–Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–CO2–H2O–TiO2–Fe2O3, to constrain conditions of regional metamorphism and the subsequent gold mineralization event. Metapelites with well‐developed amphibolite facies assemblages reliably constrain peak regional metamorphism at ~595 °C and 5.5 kbar. Observed hydrothermal alteration assemblages associated with gold mineralization in a wide compositional range of lithologies are typically calculated to be stable within P–T–X(CO2) arrays that trend towards lower temperatures and pressures with increasing equilibrium fluid X(CO2). These independent P–T–X(CO2) arrays converge and the region of overlap at ~375–425 °C and 1–2 kbar is taken to represent the conditions of alteration approaching equilibrium with a common infiltrating fluid with an X(CO2) of ~0.7. Fluid‐rock interaction calculations with M–X(CO2) diagrams indicate that the observed alteration assemblages are consistent with the addition of a single fluid phase requiring minimum fluid/rock ratios on the order of 1.  相似文献   

20.
Niobium and Ta concentrations in ultrahigh‐pressure (UHP) eclogites and rutile from these eclogites and associated high pressure (HP) veins were used to study the behaviour of Nb–Ta during dehydration and fluid‐rock interaction. Samples were collected through a ～2 km profile at the Bixiling complex in the Dabie orogenic belt, Central‐Eastern China. All but one eclogite away from veins (EAVs) display nearly constant Nb/Ta ratios ranging from 16.1 to 19.2, with an average of 16.9 ± 0.8 (2 SE), similar to that of their gabbroic protolith from the Yangtze Block. Nb/Ta ratios of rutile from the EAVs range from 12.7 to 25.3 among different individual grains, with the average values close to those of the corresponding bulk rocks. These observations show that Nb and Ta were not significantly fractionated by prograde metamorphism up to eclogite facies when no significant fluid‐rock interaction occurs. In contrast, Nb/Ta ratios of rutile from eclogites close to veins (ECVs) are highly variable from 17.8 to 49.8, which are systematically higher (by up to 17) than those of rutile from the veins. These observations demonstrate that Nb and Ta were mobilized and fractionated during localized fluid flow and intensive fluid‐rock interaction. This is strongly supported by Nb/Ta zoning patterns in single rutile grains revealed by in situ LA‐ICP‐MS analysis. Ratios of Nb/Ta in the ECV‐hosted rutile decrease gradually from cores towards rims, whereas those in the EAV‐hosted rutile are nearly invariable. Furthermore, the vein rutile shows Nb/Ta zoning patterns that are complementary to those in rutile from their immediate hosts (ECVs), suggesting an internal origin for the vein‐forming fluids. The Nb/Ta ratios of such fluids evolved from low values at the early stage of subduction to higher values at later supercritical conditions with increased temperature and pressure. Quantitative modelling was conducted to constrain the compositional evolution of metamorphic fluids during dehydration and fluid‐rock interaction focusing on Nb–Ta distribution. The modelling results based on our proposed multistage fluid phase evolution path can essentially reproduce the natural observations reported in the present study.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号