大别造山带地壳泊松比结构与超高压变质带—来自宽角反射与近垂直反射剖面的启示

综合宽角反射、近垂直反射的探测结果和有关地质资料，对大别造山带地壳结构和超高压变质带研究显示：大别造山带地壳具有层块结构特征。沿安义－庄墓剖面，上地壳有7个弹性块体，中地壳有5个，下地壳有4个。扬子与华北地块的主缝合带是超高压变质带，扬子地块与大别造山带的现今分界线是与郯庐断裂相交的太湖－马庙断裂，磨子潭－晓天断裂是大别造山带的北界，北淮阳构造带呈楔状向下插入，它与华北地块的分界是肥中断裂。郯庐断裂在中、上地壳近于直立，下地壳向西倾斜。超高压变质带的厚度为5－7km，产状向北插入到北大别块体之下，折返过程是构造就位，不是大别山的均匀抬升，折返的主运动面是水吼－五河高温剪切带。     

基于GM(1,1)优化模型的岩石边坡变形预报

岩石边坡系统是一典型灰色系统，其变形发展过程可以用灰色预测模型完成。尽管传统GM（1，1）模型预测有很多成功的实例，但是也存在一些预测偏差过大的情况，必须对其进行优化。逐步迭代法GM（1，1）模型不仅收敛速度快，而且与原始数据库序列的凹凸性保护一致。利用自编的计算程序对马步坎边坡预测测点G1沉降和开裂进行预测分析，结果表明逐步迭代法GM（1，1）优化模型计算精度较传统GM（1，1）模型和背景构造法GM（1，1）优化模型高，较好地反映了岩石边坡的变形趋势。     

Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.     

裂隙型单斜介质中弹性系数的计算及波的传播特性研究

根据Hudson等关于裂隙介质弹性系数计算的扰动理论及Bond变换矩阵原理，给出了各向同性介质中含多组垂直裂隙时等效弹性系数的计算方法。计算了含两组斜交的垂直裂隙形成的单斜各向异性介质中的等效弹性系数，并根据Christoffel方程推导、得出水平面内平面波传播的相速度和群速度随方位变化的特性。     

A non-linear stress-stiffness model for geomaterials at small to intermediate strains

A double exponential fitting model (DEFM) capable of expressing the non-linear stress-stiffness relationship of geomaterials has been proposed by Shibuya et al. (1997). The model comprises two material constants; the elastic stiffness at very small strains and the strength, together with other free parameters to determine the complete stress-stiffness relationship. In this paper, the capability of the original function used for DEFM in simulating the tangent stiffness-stress relationship of geomaterials is first discussed. Second, the methods for determining the free model parameters, as well as its conversion to obtain a stress-strain relationship are proposed. The applicability of DEFM to predicting non-linear stress-stiffness relationship is examined in detail in a total of forty-nine fitting cases of compression test data on sedimentary rock, artificial soft rock and soft clay. It is found that the DEFM is effective in expressing the non-linear stress-stiffness relationship of various kinds of geomaterials at small to intermediate strains, say less than 0.5%. The superiority of this model compared to other fitting models currently in use is also demonstrated in some of the fitting cases.     

Behaviour of a rupture of the 21 June 2000 earthquake in South Iceland as revealed in an asphalted car park

A major earthquake (M=6.6) occurred on 21 June 2000, in South Iceland. This paper presents an unusual example of left-lateral strike-slip displacement recorded in a newly asphalted car park surface through a mechanically consistent pattern of open fissures and pressure ridges resulting from simple shear and rotation. Measurement of these features allows accurate reconstruction of the local deformation. The behaviour of the asphalt layer resembles that of analogue physical models, especially in terms of rotations induced by shear deformation. It is finally shown that through a wide range of scales some basic patterns associating rotation and opposite senses of strike-slip exist in the South Iceland Seismic Zone.     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Static deformation of two welded monoclinic elastic half-spaces due to a long inclined strike-slip fault

Static deformation of two monoclinic elastic half-spaces in welded contact due to a long inclined strike-slip fault situated in one of the half-spaces is studied analytically and numerically. Closed-form algebraic expressions for the displacement at any point of the medium are obtained. The variation of the displacement at the interface with the horizontal distance from the fault is studied. The effect of anisotropy on the displacement field is examined. It is found that while the anisotropy of the source half-space has a significant effect on the displacement at the interface, the anisotropy of the other half-space has only a marginal effect.     

半无限空间中圆形孔洞周围SH波的散射

建立了求解在含有圆形孔洞的弹性半空间中SH波散射与圆形孔洞附近动应力集中问题的解析方法。利用SH波散射的对称性和多极坐标的方法，构造了一个可以预先满足半空间自由表面上应力自由边界条件的圆形孔洞对SH波散射的波函数。利用这一波函数，则可将该问题转化成对一个圆形孔洞散射的求解问题。该问题的解答最终又可归结为对一组无究代数方程组的求解问题，并可利用截断有限项的方法对其进行计算，最后给出了有关圆形也洞附近动应力集中问题的算例和数值结果，并讨论了波数与圆孔至自由边界距离变化对动应力集中的影响。     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。