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

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观测青藏高原东北缘应力场变化

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

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

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

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.     

FITS、BMP和SCR图象格式及相互转换

本文详细介绍图象的ＦＩＴＳ格式、ＢＭＰ格式和ＳＣＲ格式。ＦＩＴＳ格式已经是天文界的通用格式，几乎所有的天文软件包都支持这一格式，而ＢＭＰ格式在ＰＣ计算机上有广泛的运用，有大量ＰＣ软件支持ＢＭＰ图象的显示，处理和打印。云南天文台的ＳＣＲ格式是ＰＣ机采集ＣＣＤ图象所使用的格式。实现这三种格式的相互转换，就可以自由的将ＣＣＤ图象在ＰＣ机和工作站上进行各种处理。转换的包括图象头和图象数据的转换，其关键在于交换数据的高低字节。     

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

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.     

Application of stress release models to historical earthquakes from North China

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

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.