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

  相似文献   

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions     

Kostantinos Papadakis  Constantine Goudas 《Astrophysics and Space Science》2006,305(2):99-124

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.  相似文献   

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

多接收器电感耦合等离子体质谱法测定镉标准溶液中镉的同位素组成   总被引：1，自引：1，他引：0  

高博  梁细荣  刘颖  胡光黔  曾文  涂湘林 《岩矿测试》2008,27(5):321-324

建立了多接收器电感耦合等离子体质谱仪对Cd标准溶液同位素组成进行分析的实验方法。仪器的质量分馏校正采用Standard-Sample Bracketing法。实验结果用δ114/110Cd来表达。在此研究基础上,以SPEX Cd标准溶液为参考标准,对国外4种Cd标准溶液进行了测定。结果表明,实验测定的精度在0.07‰~0.13‰(δ114/110Cd),与目前文献报道的结果具有相似的精度。以最新SPEX Cd标准样品(δ114/110Cd=0)为基准,计算的δ114/110CdJMC、δ114/110CdSPEX-1、δ114/110CdBAM1012和δ114/110CdM櫣nster的值分别为0.55‰、0.56‰、-0.65‰和5.14‰,说明不同批次SPEX标准溶液的Cd同位素组成是明显不同的,最新的SPEX Cd标准溶液与SPEX-1 Cd的δ114/110Cd值存在着0.56‰的差别。将以SPEX Cd为参考标准的δ114/110CdBAM1012和δ114/110CdM櫣nster转化为以SPEX-1或JMC为参考标准后,得到的结果与文献报道的结果在误差范围内一致。  相似文献   

A new method of estimating the ratio between in situ rock stresses and tectonics based on empirical and probabilistic analyses   总被引：1，自引：0，他引：1  

L.I. Gonzlez de Vallejo  T. Hijazo 《Engineering Geology》2008,101(3-4):185-194

This paper describes a new procedure for assessing the ratio between in situ stresses in rock masses by means of K (K = σ_H / σ_v, being σ_H and σ_v principal stress) and tectonics for purposes of engineering geology and rock mechanics. The method combines the use of the logic decision tree and the empirical relationship between the Tectonic Stress Index, TSI, and a series of K in situ values obtained from an extensive database. The decision tree considers geological and geophysical factors affecting stress magnitudes both on the regional and local scale. The TSI index is defined by geological and geomechanical parameters. The method proposed provides an assessment of the magnitude of horizontal stresses of tectonic origin. Results for several regions of Europe are presented and the possible applications of the procedure are discussed.  相似文献   

引水工程秦巴段隧洞地应力模拟及工程地质问题     

李东林  吴树仁  韩金良  刘景儒 《地质与勘探》2008,44(5):81-86

利用Ansys有限元软件对三峡引水工程秦巴段线路不同深度、不同截面形态隧洞围岩的应力重分布情况进行模拟计算,得到圆形隧洞、城门形隧洞和马蹄形隧洞围岩的应力数值和等值线图.利用库仑一纳维尔强度准则,对花岗岩区及灰岩区隧洞围岩的剪切破坏进行分析;利用三种岩爆应力判别指标,对岩爆进行初步预测;进而对不同深度、不同截面形态的隧洞围岩的稳定性进行分析;最后对引水工程隧洞设计中截面形态的选择给出了初步建议.  相似文献   

基于Poulos弹性理论的被动桩改进算法     

陈福全  黄伟达 《岩土力学》2008,29(4):905-910

传统的Poulos弹性理论仅适合于均质土中土体侧向位移时桩的性状分析,无法考虑土的层状特性。通过引入层状地基中作用一水平集中力的广义Mindlin解和地面作用有竖向荷载时的应力和位移通解,对Poulos方法进行了改进,使之扩展到多层土中,还用于研究堆载条件下的被动桩变形和受力响应。算例分析表明,改进弹性理论要比Poulos方法更为严密、合理,提高了计算精度,应用范围也更广。  相似文献