用积分方程法模拟各向异性地层中三维电性异常体的电磁响应

本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.     

基于面积-高程积分的叶尔羌河山区流域地貌发育特征及其控制因素研究

以叶尔羌河山区流域为研究对象, 利用ArcGIS中的水文分析工具获得1357个集水区的面积-高程积分(HI)值, 定量化研究帕米尔弧形构造带东缘叶尔羌河山区流域的地貌发育特征及其控制因素。结果表明: 1)叶尔羌河山区流域的1357个集水区面积与HI均值的相关性较差(R²=0.0936), 表明集水区面积对HI均值的影响较小; 2)HI值介于0.0717~0.6485, 均值为0.4056, 整体处于地貌发展的壮年期, 是构造、岩性和气候等因素共同作用的结果; 3)HI均值排序为断裂带区域(0.4148)>整个流域(0.4056)>非断裂带区域(0.4004), 表明构造活动是地貌发育的影响因素, 其影响主要表现在断裂带的规模和活动性; 4)非断裂带区域不同岩性的HI均值排序为: 变质岩(0.491353)>冰、雪及现代冰川(0.491351)>花岗岩(0.4426)>海相沉积岩(0.4098)>陆相沉积岩(0.3890)>海陆交互相沉积岩(0.3813)>第四系沉积物(0.3598), 表明岩性是地貌发育的重要影响因素, 且岩石的抗侵蚀能力差异是影响岩性与地貌的关系的重要因素; 5)气候主要通过降水和气温控制河流径流量来影响河流的侵蚀作用和搬运作用进而影响流域内地貌发育, 多年平均年降水对流域的地貌发育的影响较小, 多年平均年气温对流域的地貌发育有特别重要的影响, 是气候对地貌发育的主要控制因素。

     

流动环境中倾斜浮力射流的积分方法

给出了基于积分方法的一种倾斜浮力射流的数学模型,且流速、密度和浓度差剖面假定为高斯分布,并导出了用来封闭该数学模型的较为完全的掺入率公式.用计算精度较高的有限分析法计算这类定解问题,计算结果与试验资料基本吻合.     

瞬态弹性波散射的变换域积分方程法

针对瞬态弹性波散射的问题,从弹性动力学问题的积分表示定理出发,采用Laplace变换的方法,得到了变换域内均质体位移场的积分方程表示;在此基础上推导了适合瞬态弹性波对异质体散射求解的变换域位移场积分方程。     

Incremental elastoplastic FEM for simulating the deformation process of suction caissons subjected to cyclic loads in soft clays

This paper presents an incremental elastoplastic finite element method (FEM) to simulate the undrained deformation process of suction caisson foundations subjected to cyclic loads in soft clays. The method is developed by encoding the total-stress-based bounding surface model proposed by the authors in the ABAQUS software package. According to the model characteristics, elastoplastic stress states associated with the incremental strains of each iteration are determined using the sub-incremental explicit Euler algorithm, and the state parameters describing the cyclic accumulative rates of strains are updated by setting state variables during the calculations. The radial fallback method is also proposed to modify the stress states outside the bounding surface to the surface during determination of the elastoplastic stress states. The stress reversals of soil elements are judged by the angle between the incremental deviatoric stress and the exterior normal vector at the image stress point on the bounding surface to update the mapping centre and state variables during cyclic loading. To assess the general validity of the method, the reduced scale model tests and centrifuge tests of suction caissons subjected to cyclic loads are simulated using the method. Predictions are in relative good agreement with test results. Compared with the limit equilibrium and quasi-static methods, the method can not only determine the cyclic bearing capacity, but can also analyse the deformation process and the failure mechanisms of suction caisson under cyclic loads in soft clays.     

地震前电磁辐射信息特征研究

反映电磁辐射信息的地震前异常变化很多,但都有一个震前平静－加强－平静然后发震的特征,为了更好地体现这种特征,根据十几年来河北廊坊电磁台网的观测资料,通过给不同强度段的各台日信息量赋予不同的权值,并把当日各台权值累加来反映当日内多台信息量集中的强弱,然后把上述权累加和用一种不连续积分的方法反映异常的持续时间,通过对监测范围内４次Ｍｓ≥５．０地震震例分析,发现该方法反映震前异常特征明显,短临效果很好,     

A note on magnetic modelling with remanence

The general problem of magnetic modelling involves accounting for the effect of both remanent magnetization and the application of an external magnetic field. However, as far as the disturbing field of a magnetic body in a magnetic environment is concerned, there is an equivalence between the effects of these two causations that allows the remanence to be represented in terms of an equivalent primary magnetic H field. Moreover, due to the linearity of the magnetic field in terms of its causations, the general modelling problem involving an applied magnetic field in the presence of remanence can be simply and more efficiently dealt with in terms of an equivalent primary field acting in the absence of any remanent magnetization.     

Reaction-Diffusion Systems and Nonlinear Waves

The authors investigate the solution of a nonlinear reaction-diffusion equation connected with nonlinear waves. The equation discussed is more general than the one discussed recently by Manne et al. (2000). The results are presented in a compact and elegant form in terms of Mittag-Leffler functions and generalized Mittag-Leffler functions, which are suitable for numerical computation. The importance of the derived results lies in the fact that numerous results on fractional reaction, fractional diffusion, anomalous diffusion problems, and fractional telegraph equations scattered in the literature can be derived, as special cases, of the results investigated in this article.     

变分原理和非线性Lyapunov稳定性判据

本文通过构造斜压大气运动的原始方程模式和两层模式的Hamilton守恒量及积分不变的Casimir泛函,并使用变分法,得到了斜压大气运动的非线性Lyapunov稳定性判据,这个判据是Arnold-Blumen-Zeng判据的推广。     

The asymmetric elastic wavefields in a model comprising a liquid layer overlying an anisotropic solid seabed due to an arbitrary source within the solid

The asymmetric three-dimensional radiation pattern and resultant elastodynamic response of stress waves in a model comprising a compressible water column overlying a transversely isotropic seabed in which a time-harmonic source acts is theoretically investigated. The use of potential functions, the Hankel transform, and a Fourier series expansion are adopted to deal with the equations of motion for both media. Closed-form integral expressions are developed for the potentials and the stress/displacement components. The expressions and introduced procedure are sufficiently flexible to incorporate various types of source loads. To evaluate the field quantities, the residue method and a robust integration scheme are utilized to handle the poles and branch points within the integrand. Any possible number of dispersive propagation modes are taken into account in the integral evaluation. The deduced velocity dispersion curves depict the characteristics of the various modes. They also indicate the existing singular points (poles) for a specific dimensionless frequency and the surface wave type associated with each pole. Numerical results are presented for the hydrodynamic pressure and displacement in the liquid layer and stress and displacement components in the solid seabed due to distributed and concentrated source excitations. The formulation and the numerical scheme are valid for calculating the wavefield anywhere within the model including both far- and near-field effects. The sensitivity of the results to different parameters is also analyzed. Both analytical and numerical comparisons with existing solutions for simpler cases are made to confirm the validity of the results. The results are especially useful in seismic hazard assessment of submarine earthquakes, landslides, and tsunamis. They can also be extended to deal with the fluid-solid-structure interaction problems.