Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

Micro-fracture instabilities in granular solids

We study the aspect of unstable behavior (like strain localization bands) in elastic solids as a consequence of micro-fracturing. A two-scale approach of computational homogenization is considered. The macroscopic behavior is investigated by finite element computations on a unit cell. At the micro-level, we consider a granular structure with elastic grains. The inter-granular boundaries are modeled with cohesive laws, friction and unilateral contact. We show that decohesion between grains gives rise to macro-instabilities, indicated by the loss of ellipticity, typical for deformation localization bands. The relation between the microscopic softening on inter-granular boundaries and the onset of macro-instabilities is studied through numerical examples. The influence of the cohesive law and friction parameters is analyzed. For periodic distributions of granular structures, we prove the loss of periodicity by failure and the corresponding size dependence effect in the homogenized response. We present numerical examples of bifurcation of solutions for granular cell structures and of particular solutions specific to elementary volumes with periodic cell distribution. Size dependence appears in the unstable regime and is strongly influenced by cohesion and friction parameters.     

极化分析研究及其在波场分离中的应用

不同的波在空间的质点运动轨迹和极化方向是不同的，利用三分量或二分量记录，在特定的以ｔ０为中心，长为Ｎ△ｔ的滑动时窗中计算出描述波的质点运动轨迹和极化方向的特征参数。利用这些参数可以设计出合适的极化滤波函数，保留那些满足特定要求的质点运动，而衰减那些不满足设计要求的质点运动，利用某些算子可以把弹性波分解成Ｐ波和ＳＶ波，而且衰减面波；利用另外一些算子可以把Ｐ波和ＳＶ波衰减而使面波增强。这些算子甚至可以在嗓音和信号具有相同频率的情况下，成功地提高信嗓比。     

新型移动三维激光扫描技术在盾构管片椭圆度检测中的应用

铁运营阶段对隧道结构的变形监测保证了地铁运行的安全,而椭圆度检测是地铁隧道结构检测的重要工作。本文简单介绍了传统椭圆度检测的基本方法,分析了新型移动三维激光扫描检测系统基本原理及隧道椭圆度检测的方法和处理流程。通过工程案例实际应用以及对检测结果的综合分析,证明了移动三维激光扫描技术在盾构管片椭圆度检测中的优势。     

椭圆度-凹坑双缺陷海底管道局部屈曲特性研究

海底管道在制造、埋设以及使用过程中极易产生椭圆度-凹坑双缺陷,双缺陷影响管道局部屈曲,对含椭圆度-凹坑双缺陷海底管道的局部屈曲特性研究十分必要。现行规范中采用等效椭圆度对含椭圆度-凹坑海底管道进行评估,该方法无法准确评估不同缺陷形式的屈曲特性。采用形状系数对含椭圆度-凹坑双缺陷的海底管道进行评估,运用有限元软件ABAQUS进行数值模拟,并进行试验验证。在此基础上对含有不同凹坑深度、不同椭圆度的海底管道进行局部屈曲的数值模拟,计算不同形状椭圆度、含有不同凹坑深度海底管道的形状系数,对其进行敏感性分析。计算结果表明:形状系数对海底管道椭圆度、凹坑深度、径厚比敏感性较强;对凹坑宽度敏感性较弱。     

ELASTOPLASTIC RESPONSE OF PRESSURE SENSITIVE SOLIDS

Tensorially invariant constitutive relations are systematically derived for large strain elastoplastic response of geomaterials. The analysis centres on Mohr–Coulomb (MC) and Drucker–Prager (DP) models with arbitrary hardening and non-associated response. Both flow and deformation theories are constructed for each model with emphasis on linear incremental relations between the Eulerian strain rate tensor and the objective Jaumann stress rate tensor. Specifying the results for plane strain compression we find that deformation theory produces a much smaller tangent instantaneous shear modulus than flow theory. It follows that failure of ellipticity and onset of surface instabilities predicted by deformation theory for associated solids occur at much lower levels of strain than the corresponding flow theory results. On the other hand, flow theory predictions admit a considerable sensitivity to the level of non-associativity. In fact, at high levels of non-associativity flow theory predictions for loss of ellipticity can be at strains below those obtained from deformation theory. © 1997 John Wiley & Sons, Ltd.     

Structural softening,mesh dependence,and regularisation in non-associated plastic flow

A severe dependence of numerical simulations on the mesh density is usually attributed to the presence of strain softening in the constitutive relation. However, other material instabilities, like non-associated plastic flow, can also cause mesh sensitivity. Indeed, loss of ellipticity in quasi-static analyses is the fundamental cause of the observed mesh dependence. It has been known since long that non-associated plastic flow can cause loss of ellipticity, but the consequence for mesh sensitivity, and subsequently, for the difficulty of the equilibrium-finding iterative procedure to converge have remained largely unnoticed. We first demonstrate at the hand of a biaxial test structural softening and a marked mesh dependence for an ideally plastic material equipped with a non-associated flow rule. The phenomena are then analysed in depth using an infinitely long shear layer. Finally, it is shown that the mesh effect disappears when the standard continuum model is replaced by a Cosserat continuum, a well-known regularisation method for strain-softening constitutive relations.     

The Effect Of Zonal Tides On The DynamicalEllipticity Of The EarthAnd Its Influence On The Nutation

In this paper, the expressions of variations of the dynamical ellipticity and the principal moments of inertia due to the deformations produced by the zonal part of the tidal potential are obtained. Starting from these expressions, we have studied from equations related to Hamiltonian theory, their effects on the nutation and finally we have evaluated numerically such influences, with a level of truncation at 0.1 μas. Thus we have shown that some coefficients are quite large with respect to the usual accuracy of up-to-date observations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Core mantle boundary topography from short period PcP, PKP, and PKKP data

We use a total of 839,369 PcP, PKPab, PKPbc, PKPdf, PKKPab, and PKKPbc residual travel times from [Bull. Seism. Soc. Am. 88 (1998) 722] grouped in 29,837 summary rays to constrain lateral variation in the depth to the core-mantle boundary (CMB). We assumed a homogeneous outer core, and the data were corrected for mantle structure and inner-core anisotropy. Inversions of separate data sets yield amplitude variations of up to 5 km for PcP, PKPab, PKPbc, and PKKP and 13 km for PKPdf. This is larger than the CMB undulations inferred in geodetic studies and, moreover, the PcP results are not readily consistent with the inferences from PKP and PKKP. Although the source-receiver ambiguity for the core-refracted phases can explain some of it, this discrepancy suggest that the travel-time residuals cannot be explained by topography alone. The wavespeed perturbations in the tomographic model used for the mantle corrections might be too small to fully account for the trade off between volumetric heterogeneity and CMB topography. In a second experiment we therefore re-applied corrections for mantle structure outside a basal 290 km-thick layer and inverted all data jointly for both CMB topography and volumetric heterogeneity within this layer. The resultant CMB model can explain PcP, PKP, and PKKP residuals and has approximately 0.2 km excess core ellipticity, which is in good agreement with inferences from free core nutation observations. Joint inversion yields a peak-to-peak amplitude of CMB topography of about 3 km, and the inversion yields velocity variations of ±5% in the basal layer. The latter suggests a strong trade-off between topography and volumetric heterogeneity, but uncertainty analyses suggest that the variation in core radius can be resolved. The spherical averages of all inverted topographic models suggest that the data are best fit if the actual CMB radius is 1.5 km less than in the Earth reference model used (i.e. the average outer core radius would be 3478 km).