隧道地震响应数值模拟研究

在进行隧道地震响应的数值模拟研究时, 横向计算范围和人工边界等对计算结果有很大的影响。本文以黄草坪隧道为研究对象, 应用有限差分程序FLAC3D对其进行地震响应的数值模拟研究, 将横向计算范围分别取为隧道洞径的5倍、6倍、7倍、8倍、9倍和10倍, 并分别采用FLAC3D中的截断边界、自由场地边界和粘性边界进行计算。研究结果表明, 当地震波为P波时, 横向计算范围取为洞径的7至8倍, 人工边界采用自由场地边界或粘性边界是比较合理的。      

郯庐断裂带(安徽段)及邻区的动力学分析与区域构造演化

依据区域构造层次划分,采用构造筛分法,层层深入,层层筛分,确定发生于各个不同时代地层/岩层内的断裂活动的同期及叠加的应力场特征。综合所有的同期应力场特征及辅以叠加的应力场特征来验证,从而确定了一个连续的、完整的断裂活动的应力场演化序列;结合区域构造变形特征分析,阐明郯庐断裂带(安徽段)的构造演化。应力场分析显示:晚三叠-早侏罗世应力场为北北西—南南东或近南北向挤压,属古特提斯构造域,断裂发生同造山走滑;早白垩世早期,应力场为北西—南东向挤压,断裂发生左行走滑运动,中国东部处于西环太平洋构造域;早白垩世晚期—古新世(始新世),区域发生北西—南东向伸展作用,断裂处于伸展断陷作用阶段;新生代,受区域上近东西向的挤压作用影响,断裂发生挤压逆冲兼右行走滑作用。     

Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.     

In situ stress field of eastern Australia

New in situ data based on hydraulic fracturing and overcoring have been compiled for eastern Australia, increasing from 23 to 110 the number of in situ stress analyses available for the area between and including the Bowen and Sydney Basins. The Bowen Basin displays a consistent north‐northeast maximum horizontal stress (σ_H) orientation over some 500 km. Stress orientations in the Sydney Basin are more variable than in the Bowen Basin, with areas of the Sydney Basin exhibiting north‐northeast, northeast, east‐west and bimodal σ_H orientations. Most new data indicate that the overburden stress (σ_V) is the minimum principal stress in both the Bowen and Sydney Basins. The Sydney Basin is relatively seismically active, whereas the Bowen Basin is relatively aseismic. Despite the fact that in situ stress measurements sample the stress field at shallower depth than the seismogenic zone, there is a correlation between the stress measurements and seismicity in the two areas. Mohr‐Coulomb analysis of the propensity for failure in the Sydney Basin suggests 41% of the new in situ stress data are indicative of failure, as opposed to 13% in the Bowen Basin. The multiple pre‐existing structural grains in the Sydney Basin further emphasise the difference between propensity for failure in the two areas. Previous modelling of intraplate stresses due to plate boundary forces has been less successful at predicting stress orientations in eastern than in western and central Australia. Nonetheless, stress orientation in the Bowen Basin is consistent with that predicted by modelling of stresses due to plate boundary forces. Variable stress orientations in the Sydney Basin suggest that more local sources of stress, such as those associated with the continental margin and with local structure, significantly influence stress orientation. The effect of local sources of stress may be relatively pronounced because stresses due to plate boundary forces result in low horizontal stress anisotropy in the Sydney Basin.     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Stress‐dependent hardening and failure surfaces of dry sand

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating the crush zone size under a cutting tool in coal

Summary As part of an effort to understand the mechanics of fine fragment formation in coal, which is important in studies of respirable dust due to mining, fracture toughness measurements and the strain energy density (SED) theory were applied to calculate the crush zone size under a cutting tool in coal. This zone is the major source of fine fragments in the 1 to 10 µm size range. The model used in these calculations is a boundary element program containing a failure criterion based on the SED theory. The boundary element program calculates linear elastic stresses at numerous points in the coal material ahead of a cutting bit. These stresses are then input to a subroutine called critical flaw length and orientation (CFLO) which uses the SED theory to determine the CFLO for a small crack at the boundary element stress computation point. The extent of crushing is based on earlier postulates about the role of inherent flaws in a fragmentation process. To form 1 to 10 µm fragments requires firstly a local stress strong enough to activate flaws with a characteristic length less than 1 to 10 µm and secondly, a flaw density sufficient to provide an average spacing between flaws also on the order of 1 to 10 µm. The locus of active 10 µm flaws represents the maximum possible extent of fine fragmentation in the 10 µm or less size range assuming that a sufficient inherent flaw density exists. The approach offers a first order approximation to the extent of crushing under a tool tip. The size and shape of the crush zone volume is affected by the attack angle and geometry of the tool.     

岩浆包裹体化学成分研究

岩浆包裹体化学成分研究难度较大,为了获得可信的数据,应当注意:1.非演化型岩浆包裹体的化学成分可以代表其初始成分。演化型岩浆包裹体应先均一、淬火后再行测定。2.均一演化型岩浆包裹体应严格遵守加热规则,否则过热作用会使包裹体壁部分熔化,造成淬火后所测包裹体成分与其真正的初始成分并不相当。3.实测资料证明,“边界层效应”对于岩浆包裹体化学成分影响微不足道。4.岩浆包裹体的化学成分只能代表其主矿物结晶时周围岩浆的成分,即仅相当于岩浆液相线上的一个点。5。把显微冷热台测温、激光喇曼探针和电子探针分析技术结合使用,对查明单个包裹体中挥发组分的性状和浓度具有很大的潜力和前途。     

边界单元法在电法勘探中的应用

本文用加权剩余法导出点源二维地电断面的边界单元法积分方程。通过付氏反变换,可获得三维空间的电位。若干模型上正演计算结果与解析解完全一致。对实测的ρ_s曲线进行了地形攻正及选择法反演拟合,解释结果与实际地质断面符合得很好。     

Active tectonics of the Himalaya

The Himalayan mountains are a product of the collision between India and Eurasia which began in the Eocene. In the early stage of continental collision the development of a suture zone between two colliding plates took place. The continued convergence is accommodated along the suture zone and in the back-arc region. Further convergence results in intracrustal megathrust within the leading edge of the advancing Indian plate. In the Himalaya this stage is characterized by the intense uplift of the High Himalaya, the development of the Tibetan Plateau and the breaking-up of the central and eastern Asian continent. Although numerous models for the evolution of the Himalaya have been proposed, the available geological and geophysical data are consistent with an underthrusting model in which the Indian continental lithosphere underthrusts beneath the Himalaya and southern Tibet. Reflection profiles across the entire Himalaya and Tibet are needed to prove the existence of such underthrusting. Geodetic surveys across the High Himalaya are needed to determine the present state of the MCT as well as the rate of uplift and shortening within the Himalaya. Paleoseismicity studies are necessary to resolve the temporal and spatial patterns of major earthquake faulting along the segmented Himalayan mountains.