Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application     

Daichao Sheng  David W. Smith  Scott W. Sloan  Antonio Gens 《国际地质力学数值与分析法杂志》2003,27(9):767-790

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

应用GPS观测青藏高原东北缘应力场变化   总被引：2，自引：0，他引：2  

张永志  胡斌  陈宏宾 《地球科学与环境学报》2003,25(4):80-84

采用各向同性弹性地球模型推导了地面位移场速率与地壳内任意点应力场变化的边界积分关系,同时利用青藏高原东北缘1999～2001年观测的GPS资料对观测区地壳深度为5 km和25 km的主应力和最大剪应力进行了计算分析.结果表明,青藏高原东北缘的主应力变化主要集中在祁连山断裂、海原断裂等,在1920～1954年间历史上发生过多次震级为7.0～8.5级强震的断裂附近,并具有主应力变化沿断层走向分布、最大剪应力沿断层走向交替变化等特征.  相似文献   

Dissolution and depletion of ferromagnesian minerals from Holocene tephra layers in an acid bog,New Zealand,and implications for tephra correlation     

A. P. W. Hodder  P. J. De Lange  D. J. Lowe 《第四纪科学杂志》1991,6(3):195-208

This study examines the depletion of ferromagnesian silicate minerals from a sequence of thin, distal, mainly rhyolitic tephra layers of Holocene age preserved in an acid peat bog (Kopouatai), North Island, New Zealand. The rate of such depletion has been fast, as indicated by the complete loss of biotite from one tephra layer (Kaharoa Tephra), in which it is normally dominant, in only ca. 770 yr. Chemical dissolution is advocated as the likely cause for the depletion, with amphiboles and other mineral grains commonly showing etch pits, microcaves, and other characteristic surface solution features. Theoretical thermodynamic and kinetic models show a marked increase in the rate of dissolution of all ferromagnesian minerals under conditions of low pH (< 4), but that where silica concentrations in solution are high the relative proportions of minerals remaining are unaffected. However, where concentrations of dissolved silica are low, as in most bog environments, the relative proportions of ferromagnesian minerals are affected as well as absolute amounts being decreased. Amphiboles are depleted relative to pyroxenes, consistent with kinetic studies. The results show that the identification and correlation of tephras on the basis of relative abundances of ferromagnesian minerals alone may be unreliable, and emphasise the need to use multiple criteria in such studies.  相似文献   

Microboudin structures of piemontite along the Sambagawa metamorphic belt, Japan: implications for lateral variation of differential stress during regional metamorphism     

T. Masuda  T. Shibutani  T. Ochiai  S. Akagi  H. Yamaguchi  Y. Kugimiya  N. Kimura  T. Miyake 《Journal of Metamorphic Geology》2004,22(3):199-205

Evidence is presented of a lateral variation in differential stress during metamorphism along a regional metamorphic belt on the basis of the proportion of microboudinaged piemontite grains (p) in a quartz matrix in metacherts. It is proposed that p is a practical indicator of relative differential stress. Analysis of 123 metacherts from the 800 km long Sambagawa metamorphic belt, Japan, reveals that p‐values range from < 0.01 to 0.7 in this region. Most samples from Wakayama in the mid‐belt area have p‐values of 0.4–0.6, whereas those from western Shikoku have p‐values of < 0.1. This difference cannot be explained by variations in metamorphic temperature, and is instead attributed to a regional, lateral variation in differential stress during metamorphism.  相似文献   

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion     

Chengqing Wu  Yong Lu  Hong Hao 《国际地质力学数值与分析法杂志》2004,28(1):93-109

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.  相似文献   

Viscoelastodynamics of A Stratified, Compressible Planet: Incremental Field Equations and Short- and Long-Time Asymptotes     

Detlef Wolf 《Geophysical Journal International》1991,104(2):401-417

  相似文献   

Shear and objective stress rates in hypoplasticity     

D. Kolymbas  I. Herle 《国际地质力学数值与分析法杂志》2003,27(9):733-744

  相似文献   

Application of stress release models to historical earthquakes from North China   总被引：5，自引：0，他引：5  

Xiao-Gu Zheng  David Vere-Jones 《Pure and Applied Geophysics》1991,135(4):559-576

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.  相似文献   

Focal mechanism inversion in the Giudicarie–Lessini seismotectonic region (Southern Alps, Italy): Insights on tectonic stress and strain     

Alfio Vigan  Gianni Bressan  Giorgio Ranalli  Silvana Martin 《Tectonophysics》2008,460(1-4):106-115

A set of 41 focal mechanisms (1989–2006) from P-wave first polarities is computed from relocated seismic events in the Giudicarie–Lessini region (Southern Alps). Estimated hypocentral depths vary from 3.1 to 20.8 km, for duration magnitudes (M_D) in the range 2.7–5.1. Stress and strain inversions are performed for two seismotectonic zones, namely G (Giudicarie) and L (Lessini). This subdivision is supported by geological evidence, seismicity distribution, and focal mechanism types. The available number of data (16 in G, 22 in L) does not make possible any further subdivisions. Seismotectonic zones G and L are undergoing different kinematic regimes: thrust with strike-slip component in G, and strike-slip in L. Principal stress and strain axes in each sub-region show similar orientations. The direction of maximum horizontal compressive stress is roughly perpendicular to the thrust fronts along the Giudicarie Belt in zone G, and compatible with right-lateral strike-slip reactivation of the faults belonging to the Schio-Vicenza system in zone L. On the whole, kinematic regimes and horizontal stress orientations show a good fit with other stress data from focal mechanisms and breakouts and with geodetic strain rate axes.  相似文献   

Diagenesis, palaeoclimate and tectono-sedimentary influences on clay mineralogy and stable isotopes from Upper Cretaceous marine successions of the Basque-Cantabrian Basin (N Spain)     

lvaro Jimnez Berrocoso  Maria Cruz Zuluaga  Javier Elorza 《Cretaceous Research》2008,29(3):386-404

Clay mineralogy and whole-rock stable isotopes (δ¹⁸O and δ¹³C) of Upper Cretaceous marly sediments on the Basque-Cantabrian Basin have been integrated to determine the main effects of diagenesis, palaeoclimate and tectono-sedimentary factors in sections belonging to deep- (Barrika) and platform-marine (Isla de Castro, Villamartín and Olazagutía) settings.The mean values for the clay assemblages and δ¹⁸O exhibit notable differences among the sections, partially explainable by the influence of diagenesis. The Barrika sediments, with more diagenetically advanced illite-smectite (I-S) mixed-layer (R1, 70% illite), authigenic chlorite, and low δ¹⁸O (−4.05‰ PDB), experienced higher diagenetic grade than Isla de Castro and Olazagutía, which have R0 I-S (20% illite) and heavier δ¹⁸O. Villamartín was also affected by higher diagenesis than Isla de Castro and Olazagutía, given the occurrence of R1 I-S (60% illite) and low δ¹⁸O (−4.11‰ PDB). However, the absence of other clays in Villamartín (e.g. authigenic chlorite) is indicative of less diagenetic grade than Barrika. These results show the useful integration of clay mineralogy and stable isotopes to detect different diagenetic grades in distinct marine successions of the same basin.Despite being influenced by diagenesis, the clay mineralogy partially preserves its inherited signature. This allows detection of major contents of I-S and mica, and minor kaolinite, interpreted as indicative of warm palaeoclimatic conditions. High kaolinite content in Villamartín and absence of kaolinite in Isla de Castro, though, are considered to be a product of neither diagenesis nor palaeoclimatic influences. Instead, tectono-sedimentary causes, related to unsuitable conditions for clay formation and transport from the local source areas, contributed to original clay differences. The inferred effects of diagenesis, palaeoclimate and tectono-sedimentary factors make this work important to show the potentially great variety of controls on the clay mineralogy of marine sections, which are often uncritically treated in studies concerning the Late Cretaceous.  相似文献