首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 203 毫秒
1.
 A thermoelastic model for calculating the high-pressure and high-temperature properties of isotropic solids is presented by extending the formalism by Thomsen and combining the resulting one with the Vinet model for static lattice and the Debye model for lattice vibration. Applying it to polycrystalline corundum, we have shown that the calculated values of entropy and heat capacity at constant pressure are in agreement with literature values to 2325 K at zero pressure and that the calculated values of thermal expansivity agree reasonably with experimental data to 1100 K at zero pressure. The model reproduces experimental data of sound velocities v p and v s of compressional and shear waves to 1825 K at zero pressure and those to 62 GPa at room temperature, and it reproduces also experimental shock-wave equation of state to 150 GPa. The velocity correlation (∂ln v s /∂ln v p ) S was found to have weak pressure and temperature dependences and the results under lower mantle conditions are compared with those of magnesian and calcium silicate perovskites and magnesiowüstite, and the PREM values of the Earth's lower mantle. Received: 12 February 2000 / Accepted: 15 July 2000  相似文献   

2.
In this paper, an attempt has been made to highlight the influence of different parameters such as number of cycles, confining pressure, void ratio, gradation, initial anisotropy and stress path on the dynamic properties of granular materials using Discrete Element Method (DEM). A series of strain controlled cyclic triaxial numerical simulations using three dimensional DEM have been carried out on an assembly of spheres. Dynamic properties such shear modulus (G) and damping ratio (D) were determined from the typical hysteresis loop obtained during cyclic triaxial test simulation. It has been observed from the test results that the numerical simulation using DEM has captured the variation of dynamic properties over a wide range of shear strain values for different parameters considered for the current investigation. Maximum shear modulus (G max) was found to be influenced by initial confining pressure, void ratio, gradation and initial anisotropy. Whereas, the damping ratio (D) was found to be influenced by number of cycles, initial confining pressure, gradation and stress path. Further it has been shown that the variation of shear modulus with shear strain can be divided into three distinct zones such as Isotropic Zone (IZ), Anisotropic Zone (AZ) and Stable Anisotropic Zone (SAZ). A drastic reduction of shear modulus with shear strain has been observed in the Anisotropic Zone (AZ). In addition, the results obtained using numerical simulations have been compared with the laboratory experimental values.  相似文献   

3.
Pressure is one of the most important parameters to be quantified in geological problems. However, in metamorphic systems the pressure is usually calculated with two different approaches. One pressure calculation is based on petrological phase equilibria and this pressure is often termed thermodynamic pressure. The other calculation is based on continuum mechanics, which provides a mean stress that is commonly used to estimate the thermodynamic pressure. Both thermodynamic pressure calculations can be justified by the accuracy and applicability of the results. Here, we consider systems with low‐differential stress (<1 kbar) and no irreversible volumetric deformation, and refer to them as conventional systems. We investigate the relationship between mean stress and thermodynamic pressure. We discuss the meaning of thermodynamic pressure and its calculation for irreversible processes such as viscous deformation and heat conduction, which exhibit entropy production. Moreover, it is demonstrated that the mean stress for incompressible viscous deformation is essentially equal to the mean stress for the corresponding viscous deformation with elastic compressibility, if the characteristic time of deformation is five times longer than the Maxwell viscoelastic relaxation time that is equal to the ratio of shear viscosity to bulk modulus. For typical lithospheric rocks, this Maxwell time is smaller than c. 10,000 years. Therefore, numerical simulations of long‐term (>10 kyr) geodynamic processes, employing incompressible deformation, provide mean stress values that are close to the mean‐stress value associated with elastic compressibility. Finally, we show that for conventional systems the mean stress is essentially equal to the thermodynamic pressure. However, mean stress and, hence, thermodynamic pressure can be significantly different from the lithostatic pressure.  相似文献   

4.
The role of the uppermost mantle strength in the pattern of lithosphere rifting is investigated using a thermo-mechanical finite-element code. In the lithosphere, the mantle/crust strength ratio (SM/SC) that decreases with increasing Moho temperature TM allows two strength regimes to be defined: mantle dominated (SM > SC) and crust dominated (SM < SC). The transition between the two regimes corresponds to the disappearance of a high strength uppermost mantle for TM > 700 °C. 2D numerical simulations for different values of SM/SC show how the uppermost mantle strength controls the style of continental rifting. A high strength mantle leads to strain localisation at lithosphere scale, with two main patterns of narrow rifting: “coupled crust–mantle” at the lowest TM values and “deep crustal décollement” for increasing TM values, typical of some continental rifts and non-volcanic passive margins. The absence of a high strength mantle leads to distributed deformations and wide rifting in the upper crust. These numerical results are compared and discussed in relation with series of classical rift examples.  相似文献   

5.
The fully temperature-dependent model of the effective pressure of the solid matrix and its related overpressure has been derived from the pressure balance equation, mass conservation, and Darcy’s law, and is directly useful in basin modeling. Application of the model in the Kuqa Depression of the Tarim Basin in western China proves that this overpressure model is highly accurate. The case of the present-day values of the calculated overpressure histories of Wells Kela2 and Yinan2 approach the field-measured data with mean absolute relative residuals of 3% and 5%, respectively. This indicates that the overpressure simulation is a practical alternative to using rock mechanics experiments for effective pressure measurement. Since calculation of the overpressure history uses the geohistory model and geothermal history model simulation outcomes, the relevant data used and the output of the two models of the Kela2 well are given as examples. The case studies show that the pore fluid density and viscosity used in the calculation of overpressures should be temperature-dependent, otherwise the calculation results would deviate far from the field-measured pressure data. They also show that the most sensitive parameter governing overpressure is permeability, and permeability can be calculated by using either the Kozeny–Carman formula or the porosity–power function. The Kozeny–Carman formula is better if accurate data for the specific surface area of the solid matrix (S a ) exists, otherwise, the porosity–power function is used. Furthermore, it is vital for calculating an accurate overpressure history that one can calibrate S a in the Kozeny–Carman formula, or index m in the porosity–power function by using field-measured pressure data as a constraint. In these specific case studies, the outcome found by using the Kozeny–Carman formula approaches the outcome found by using the porosity–power function with m=4, and both approach the field-measured pressure data.  相似文献   

6.
Masanori Kameyama   《Tectonophysics》2003,376(3-4):185-197
We investigated the similarity between thermal–viscous coupling (TVC) and frictional sliding, proposed by Kameyama and Kaneda [Pure Appl. Geophys. 159 (2002) 2011]. We consider a one-dimensional layer composed of viscous material, which is sandwiched and sheared by two thick elastic layers. The rate of viscous deformation depends on the temperature Tc in the viscous layer as well as shear stress τ. The temperature Tc changes owing to heating by viscous dissipation and conductive cooling. We carried out velocity-stepping tests for the steady-state deformation both numerically and analytically, and compared the temporal evolution of small perturbations with that of the spring-block model with rate- and state-dependent friction (RSF). We found that, as is the case of frictional slip stability, the manner of temporal evolution is classified into four regimes depending on whether it is stable or not and whether it is monotonous or oscillatory with time. By further interpreting TVC in terms of general RSF theory by Ruina [J. Geophys. Res. 88 (1983) 10359], we obtained the relations between the parameters appearing in the phenomenological RSF law and the nondimensional parameters which characterize the nature of TVC. A further improvement of this approach might be important for estimating the actual values of frictional constitutive parameters at the deeper portion of seismogenic faults of interplate or inland earthquakes where a ductile deformation is expected to be significant.  相似文献   

7.
成岩成矿深度的构造校正测算方法,是从测算压力中先消除掉构造附加静水压力之后再计算上覆岩石厚度,即成岩成矿深度的方法。该方法基于对地壳岩石处于固体应力状态的认识之上,采用弹性固体模型代替静止流体模型,对“上覆岩石重力在数值上等同于该点所承受的静水压力”这一通常的认识提出了不同见解,比沿袭至今单纯用压力/比重(或密度)方法得出深度更符合于实际情况。本文介绍了该方法的理论基础和野外地质研究方法——开展变形岩相形迹填图,在室内利用三维变形和古差应力测量,根据样品所处构造部位和性质,选择不同的参数换算成矿时差应力的众值。以胶东玲珑——焦家式金矿床为例,求得成矿深度仅3.5km或更浅,进而提出更深部位存在深部金矿富集带的预测意见。胶东几个大型金矿深部第二富集带已揭露的勘探资料证实这一认识比较符合实际情况。用这一方法测算出大别超高压带含柯石英榴辉岩形成深度仅32km多,而不是用压力/比重方法估算的100多公里,这为大别造山带的构造格局和演化历史的研究提出新的途径和方法。  相似文献   

8.
We have used incompressible Navier–Stokes in 2D finite element modelling to investigate rigid inclusion rotation under confined bulk simple shear flow. Confinement is defined as the ratio (S) between the channel width (H) and the inclusion's least axis (e2)(S=H/e2). The numerical results show that (i) inclusion rotation is strongly influenced by S and, when the confinement is effective, aspect ratio (R) and shape also play an important role. (ii) Back rotation is limited because inclusions reach a stable equilibrium orientation (se). (iii) There is also an unstable equilibrium orientation (ue), which defines an antithetic rotation field with se, and both se and ue depend on S, and inclusion R and shape.  相似文献   

9.
This paper presents a numerical investigation into mechanical behavior and strain localization in methane hydrate (MH) bearing sediments using the distinct element method (DEM). Based on the results of a series of laboratory tests on the bonded granules idealized by two glued aluminum rods and the available experimental data of methane hydrate samples, a pressure and temperature dependent bond contact model was proposed and implemented into a two-dimensional (2D) DEM code. This 2D DEM code was then used to numerically carry out a series of biaxial compression tests on the MH samples with different methane hydrate saturations, whose results were then compared with the experimental data obtained by Masui et al. [9]. In addition, stress, strain, void ratio and velocity fields, the distributions of bond breakage and averaged pure rotation rate (APR) as well as the evolution of strain localization were examined to investigate the relationships between micromechanical variables and macromechanical responses in the DEM MH samples. The numerical results show that: (1) the shear strength increases as methane hydrate saturation SMH increases, which is in good agreement with the experimental observation; (2) the strain localization in all the DEM MH samples develops with onset of inhomogeneity of void ratio, velocity, strain, APR, and distortion of stress fields and contact force chains; and (3) the methane hydrate saturation affects the type of strain localization, with one shear band developed in the case of 40.9% and 67.8% methane saturation samples, and two shear bands formed for 50.1% methane saturation sample.  相似文献   

10.
The characteristic relationships among relative permeability (K r ), saturation (S) and capillary pressure (P) of NAPLs are the important constitutive laws to simulate the NAPLs flow in the subsurface. In this study, a micro model was used to obtain the values of permeability, saturation for the multi-phase flow of five fluid-pairs. The perspective micro model allows one to clearly observe the multiphase flow and allow this study to precisely measure the fluid saturation by digital image analysis. The experimental results showed hysteresis phenomenon of relative permeability versus saturation and that was not interpreted by previous studies. By carefully examining the recorded images, this study found that the degree of the connectivity for the micro channel occupied by wetting phase fluid could influence the relative permeability. Therefore, for the same saturation, the relative permeability in the imbibition is higher than that in the drainage. The results of the K r –S experiments for five fluid-pairs also showed that the residual saturation of wetting phase fluid decreased with the wettability increasing but increased with the interfacial tension increasing. Those interpretations and experimental results are valuable references for groundwater remediation and oil reservoir development.  相似文献   

11.
To deeply understand the cracking mechanical behavior of brittle rock materials, numerical simulations of a rock specimen containing a single preexisting crack were carried out by the expanded distinct element method (EDEM). Based on the analysis of crack tips and a comparison between stress- and strain-based methods, the strain strength criterion was adopted in the numerical models to simulate the crack initiation and propagation processes under uniaxial and biaxial compression. The simulation results indicated that the crack inclination angle and confining pressure had a great influence on the tensile and shear properties, peak strength, and failure behaviors, which also showed a good agreement with the experimental results. If the specimen was under uniaxial compression, it was found that the initiation stress and peak strength first decreased and then increased with an increasing inclination angle α. Regardless of the size of α, tensile cracks initiated prior to shear cracks. If α was small (such as α ≤ 30°), the tensile cracks dominated the specimen failure, the wing cracks propagated towards the direction of uniaxial compression, and the propagation of shear cracks was inhibited by the high concentration of tensile stress. In contrast, if α was large (such as α ≥ 45°), mixed cracks dominated the specimen failure, and the external loading favored the further propagation of shear cracks. Analyzing the numerical results of the specimen with a 45° inclination angle under biaxial compression, it was revealed that lateral confinement had a significant influence on the initiation sequence and the mechanical properties of new cracks.  相似文献   

12.
The orientation of the straight internal foliation Si within large ( 5 mm) garnet porphyroblasts has been measured relative to the orientation of the external foliation Se around a single antiform of 0.5 m wavelength, which folds the dominant regional foliation. The internal foliation is not constant in orientation, but varies consistently both with position around the fold and with the porphyroblast ellipticity. The dip of Si (hinge dip taken as zero) is consistently less than the dip of Se; it increases with increasing dip of Se and with increasing ellipticity of the porphyroblasts. Si effectively defines a fold with an opening angle greater than that in the external foliation. The opening angle of this fold in Si decreases with increasing porphyroblast ellipticity. The observed variation in the orientation of Si can be explained qualitatively by a flattened flexural flow model for fold development, as could be expected for folding of a pre-existing, strongly anisotropic foliation. The measurements clearly demonstrate that rotation of porphyroblasts relative to geographical co-ordinates did occur during the development of this fold and that a model based on the classical theories of rotation of stiff inclusions in a weaker viscous matrix is most appropriate.  相似文献   

13.
Abstract Variation in the state of stress during heterogeneous deformation should be reflected in variation in the effective pressure of metamorphic reactions, whether this is mean stress or the normal stress acting across the reacting interface. The magnitude of this pressure variation will determine whether it is discernible in the preserved metamorphic mineral assemblages of heterogeneously deformed rocks. The magnitude of the mean stress difference across a non-slipping interface between two materials with viscosity ratio >c. 20:1 is effectively equal to the maximum shear stress for flow in the more viscous material. Progressive shortening of the interface results in a higher mean stress in the more competent material, whereas extension results in a lower mean stress. For high-P/low-T eclogite facies conditions, current experimental data indicate that clinopyroxene- and garnet-rich layers of eclogite should be very strong and that pressure differences of up to 800 MPa (8 kbar) between competent layer and weaker matrix may be possible. Such high values can be obtained in widely separated competent layers for values of bulk stress in the overall multilayer that are much lower (by a factor approaching the viscosity ratio). Extrusion of material between more rigid plates, which has been proposed as a regional mechanism of lateral ‘continental escape’for both the Alps and the Himalayas, should also be accompanied by a lateral gradient in effective pressure; otherwise extrusion could not occur. Maximum mean stresses with magnitudes that are many times the maximum shear stress required for plastic flow should develop for deformation zones that are long relative to their width (e.g. around 20 times for a width-to-thickness ratio of 10). Tectonic overpressure in progressively shortened competent layers, particularly in regions of extrusion between more rigid plates, might help explain the occurrence of isolated layers and pods of low-T eclogite (<550°C) with estimated peak pressures markedly in excess of those in the surrounding matrix. It cannot explain the occurrence of isolated high-T eclogites, because at temperatures >c. 550°C, the dramatic weakening of clinopyroxene in the power-law creep field precludes the development of significant overpressures in eclogite layers.  相似文献   

14.
为了研究处于自由面以下完全淹没状态的水平圆柱在强迫振荡运动时的水动力特性,采用基于黏性流理论建立的二维两相流数值波浪水槽模型,对不同液相黏性条件下强迫振荡水平圆柱的受力进行计算,并对压力、黏性切力和圆柱运动之间的相位关系特征进行对比和分析,进而结合流场分析解释黏性影响机理。结果表明:黏性切力和涡旋压力对流体作用力的贡献差别是导致不同流体黏性下流体作用力结果差异的主要原因;涡旋运动相对圆柱振荡运动的滞后性受流体黏性影响显著,导致不同流体黏性下压力之间有相位差;流体水质点相对于圆柱的滞后运动在大黏性流体中更为显著,这导致了其黏性切力的相位超前现象。  相似文献   

15.
The measurements of the parameters of split shear (S) waves from local deep-focus earthquakes recorded in 2005–2007 by a network of 12 seismic stations in Southern Sakhalin are presented. The results revealed the heterogeneous distribution of the anisotropic properties beneath Southern Sakhalin. The azimuths of the fast S-wave polarization beneath the stations in the central part of the peninsula are oriented along the NNW and NNE-NE directions normal to and along the Kuril Trench. Beneath the stations located along the western and eastern coasts, the azimuths of the fast S-wave polarization change their direction from NNW in the northern area to E-SE in the southern area. The highest anisotropy degree (up to 0.9–1.5%) is recorded beneath the central part of Southern Sakhalin. The maximum values of the discrepancy in the arrival time of the split S-waves are observed when the azimuth of the fast S-wave is oriented along the NNE beneath the active fault zones. The analysis of the variations of the S-wave lag time shows their weak depth dependence. The highest anisotropy is assumed in the upper layers of the medium (down to a depth of about 250 km). The results obtained for the dominating wave frequency of 1–5 Hz represent mainly the medium-scale anisotropy of the top of the studied region.  相似文献   

16.
为了研究处于自由面以下完全淹没状态的水平圆柱在强迫振荡运动时的水动力特性,采用基于黏性流理论建立的二维两相流数值波浪水槽模型,对不同液相黏性条件下强迫振荡水平圆柱的受力进行计算,并对压力、黏性切力和圆柱运动之间的相位关系特征进行对比和分析,进而结合流场分析解释黏性影响机理。结果表明:黏性切力和涡旋压力对流体作用力的贡献差别是导致不同流体黏性下流体作用力结果差异的主要原因;涡旋运动相对圆柱振荡运动的滞后性受流体黏性影响显著,导致不同流体黏性下压力之间有相位差;流体水质点相对于圆柱的滞后运动在大黏性流体中更为显著,这导致了其黏性切力的相位超前现象。  相似文献   

17.
Grüneisen’s parameters are central to studies of Earth’s interior because these link elastic data to thermodynamic properties through the equation of state and can be measured using either microscopic or macroscopic techniques. The original derivation requires that the mode Grüneisen parameter (γi) of the longitudinal acoustic (LA) mode equals the thermodynamic parameter (γth) for monatomic solids. The success of the Debye model indicates that γLA = γth is generally true. Available elasticity data for crystalline solids contain 30 reliable measurements, covering 10 structures, of the pressure derivatives of the bulk (KS) and the shear (G) moduli. For these phases, the measured values of γth and γLA agree well. Other solids in the database have disparate γLA values, suggesting large experimental uncertainties within which γLA = γth. This relationship allows inference of the pressure (P) derivative of the shear modulus (∂G/∂P = G′) from widely available measurements of γth, the isothermal bulk modulus (KT), ∂KT/∂P, and G. We predict G′ as 1.55 for stishovite, 1.6 to 2.15 for MgSiO3 ilmenite, 1.0 for γ-Mg1.2Fe0.8SiO4, and 0 for FeS (troilite). Similarly, G′ measured for MgSiO3 perovskite suggests that KS′ = 4, corroborating volume-pressure data. For many materials, pairs of G′ and KS′ = ∂KS/∂P from independent elasticity studies of a given phase define a nearly linear trend, suggesting systematic errors. Non-hydrostatic conditions and/or pressure calibration likely cause the observed variance in KS′ and G′. The best values for pressure derivatives can be ascertained because the trend defined by measured pairs of G′ with KS′ intersects the relationship of G′ to K′ defined by γLA = γth at a steep angle. Our results for isostructural series show linear correlations of KS′ with KS and of G′ with G. Values of KS′ are nearly 4 for high-pressure phases, which is consistent with the harmonic oscillator model, whereas G′ has a wide range of −1 to 3. Hence, inference of a detailed mineralogy inside the Earth is best constrained by comparing seismic determinations of shear moduli to laboratory measurements.  相似文献   

18.
An analysis of the vane test using an Arbitrary Lagrangian–Eulerian formulation within a finite element framework is presented. This is suitable for soft clays for which the test is commonly used to measure in situ undrained shear strength. Constitutive laws are expressed in terms of shear stress–shear strain rates, and that permits the study of time effects in a natural manner. An analysis of the shear stress distributions on the failure surface according to the material model is presented. The effect of the constitutive law on the shear band amplitude and on the position of the failure surface is shown. In general, the failure surface is found at 1–1·01 times the vane radius, which is consistent with some experimental results. The problem depends on two dimensionless parameters that represent inertial and viscous forces. For usual vane tests, viscous forces are predominant, and the measured shear strength depends mainly on the angular velocity applied. That can explain some of the comparisons reported when using different vane sizes. Finally, the range of the shear strain rate applied to the soil is shown to be fundamental when comparing experimental results from vane, triaxial and viscosimeter tests. Appart from that, an experimental relation between undrained shear strength and vane angular velocity has been reproduced by this simulation. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

19.
The Australian continent has an enigmatic present-day stress pattern with considerable regional variability in maximum horizontal stress (SHmax) orientations. Previous attempts to estimate the Australian SHmax orientation with geomechanical–numerical models indicate that plate boundary forces provide the major controls on the contemporary stress orientations. However, these models do not satisfactorily predict the observed stress orientation in major basins throughout eastern Australia, where the knowledge of the present-day crustal stresses is of vital importance for development and management of different types of geo-reservoirs. In addition, a new comprehensive stress-data compilation in Australia, which contains 2150 data records and is the key dataset for model calibration, provides motivation to construct a new geomechanical–numerical model for Australia. Herein, we present a 3D geomechanical–numerical model that predicts both the SHmax orientation and the relative stress magnitudes throughout the Australian continent. Our best-fit model, with mean absolute deviation of 15°, is in good agreement with observed SHmax orientations and the stress regime in most areas, and shows a much better fit in areas where the stress pattern was unable to be predicted by previous published attempts. Interestingly, the best-fit model requires a significant push from the western boundary of Australian continental model, which is possible supporting evidence for the east–west-oriented mantle drag postulated by state-of-the-art global convection models, or may be generated by the excess of gravitational potential energy from Tibetan Plateau, transferred through the Indo-Australian Plate. Hence, our modelling results provide a good first-order prediction of the stress field for areas where no stress information is currently available and can be used to derive initial and boundary conditions for local and reservoir-scale 3D geomechanical models across Australia.  相似文献   

20.
Plant roots can help to stabilise slopes. Existing analytical models to predict their mechanical contribution are however limited: they typically focus on the ultimate limit state, employ various empirical factors, and ignore much of the underlying root-soil interaction. A new model was developed based on large deflection Euler-Bernoulli elastic beam theory that can be used to study the mobilisation of root strength under various loading conditions (direct shear and pull-out). Both lateral and axial loading of the root by the soil were incorporated, based on existing methodologies for foundation piles (p-y and t-z curves). The model is able to take the key parameters into account (root biomechanical properties, root architectural properties, and soil properties) while remaining quick to solve using a numerical boundary value problem solver. The model was compared with experimental direct shear test data using various root analogues (rubber, plastic, and wood) in dry sand with various densities and effective stress levels and was able to accurately predict the measured shear force-displacement behaviour. Comparison with experimentally measured pull-out force-displacement curves using rubber and wooden root analogues with various architectures in dry and partially saturated sands was also satisfactory. In the future, this model can aid with addressing long-standing problems in the root-reinforcement community: quantifying the effect of (sequential) mobilisation of root strength in direct shear, the effect of the angle at which the root crosses a shear plane, the effect of root topology on root-reinforcement or the effect of root bending, and root shear shear forces on root-reinforcement.  相似文献   

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

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