三维地形大地电磁场的边界元模拟方法

提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果.      

基于ILU-BICGSTAB算法的起伏地形MT二维正演模拟

针对起伏地形对异常体的探测效果的影响,利用有限单元法,基于MTLAB平台编写了大地电磁二维正演程序,为更好地模拟实际地形,采用矩形网格中进一步三角细化的剖分方式,并利用不完全LU分解的预条件BICGSTAB算法提高计算速度与精度,最后通过建立多种地电模型分析讨论了起伏地形对大地电磁场的影响,验证了该方法的实用性与准确性。     

三维复电阻率电磁场正演研究

将Cole—Cole模型表示的地下异常体各剖分单元复电阻率,引入到体积分方程求解均匀大地三维电磁场计算中,实现电偶源地面激发、接收复电阻率三维电磁场正演模拟。     

有限元计算机助三角剖分：介绍一种实用的三角剖分软件CAFEM

依据“鼠标人工给点，ＺＪＰ算法联网”的思路，设计了用于有限元法计算中剖分工作的机助三角单元剖分软件ＣＡＦＥＭ，克服了过去手工剖分费时、费力、出错率高的缺点。     

起伏地形条件下二维大地电磁各向异性正演

在大地电磁资料处理和解释中,大地电磁的各向异性正演一直是国内、外研究的前沿课题。这里首先从麦克斯韦方程组出发,得出各向异性介质二维大地电磁场的边值问题以及等价的变分问题,进而对计算区域进行完全非结构化三角形剖分,在单元内采用线性插值,将变分方程转化成线性方程组,求解出有限单元法数值解,获得各向异性介质下的电磁场值。讨论了电性主轴与笛卡尔坐标轴之间的夹角以及地形的起伏变化对视电阻率和相位值的影响,结果表明:起伏地形条件下,TE和TM两种极化模式的视电阻率和相位值都会受到影响,且TM模式受地形影响较大;由于介质的各向异性,视电阻率和相位断面图出现明显不同于各向同性的形态,导致横向分界面模糊。     

基于非结构化网格的大地电磁2D自适应有限元正演模拟

有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。     

基于Open GL条件随机模拟三维模型可视化

利用两阶段建模思路建立条件随机模拟模型。第一阶段在地质体表面利用Denaunay三角网剖分算法实现三角网网格化，第二阶段用Denaunay四面体剖分算实现地质体的四面体网格化，并对四面体单元进行三角面插值，进而得到三角网，最后用OpenGL实现地质体三维模型的可视化。     

潜水面变动的空间四面体单元网格自动生成方法

提出了一种由平面三角网格自动生成三维空间四面体单元网格的方法,其基本思路是将含水层平面分成三角形单元网格,将每个三角形沿垂直方向对应一个三棱柱,作为最初的三棱柱;依次将每个三棱柱按照潜水面位置和分层信息划分成不同数目的四面体,从而生成三维空间四面体单元网格。剖分结果表明,此方法综合考虑了潜水面的位置和各含水层顶底板的起伏情况,能够简洁地生成层面起伏的空间四面体网格,特别适用于刻画潜水面波动、含水层顶底板起伏和存在透镜体或岩性缺失等含水层结构复杂的情况下,自动地剖分并生成单元信息,为三维计算的实现奠定了基础。     

矿集区大地电磁噪声场源分析及噪声特点

结合长江中下游某矿集区大地电磁观测环境,分析大型矿山、通讯基站及其他主要电磁噪声源对视电阻率曲线的影响.依据时域波形将噪声分为脉冲噪声、周期噪声、三角波噪声、方波噪声、阶跃噪声、似充放电模式的三角噪声,便于在后期资料处理过程中针对不同的噪声类型对大地电磁数据进行降噪处理.     

无单元Galerkin法大地电磁三维正演模拟

无单元Galerkin法(EFGM)作为一种相对成熟的无网格方法,避免了网格剖分,其精度高,适用于复杂电导率分布和复杂边界形状的计算。本文将EFGM用于大地电磁三维正演,详述了三维EFGM形函数的构造过程,从大地电磁三维变分问题出发,利用Galerkin法结合高斯积分公式推导了相应的系统矩阵离散表达式,简述了边界条件的加载技术,研究了支持域尺寸对EFGM三维正演计算精度的影响,最后通过数值计算验证了EFGM三维算法的正确性。     

Risk evaluation of China’s natural disaster systems: an approach based on triangular fuzzy numbers and stochastic simulation

To deal effectively with the evaluation problem of natural disaster risk system affected by many uncertain factors, a multivariate connection number expression is presented. This expression is based on the index samples and evaluation grade criterions of natural disaster risk system and is capable of describing the hierarchy property and fuzziness of membership relationship between index samples and evaluation grade criterions. In this proposed method, the fuzzy evaluation grade criterion problem is resolved by combining triangular fuzzy numbers with multivariate connection number theory, and triangular fuzzy numbers are used to express the discrepancy degree coefficients of connection number and evaluation index weights. Accordingly, a connection number-based evaluation method for the natural disaster system of China (named CN-TFN for short) is established using triangular fuzzy numbers and stochastic simulation. The application results show that the spatial distribution of natural disaster risk grades of China has the trend of aggrandizement from west to east of China. The economically developed and densely populated coastal areas are very likely to have a high level of natural disaster risk grade or above; thus, these areas are the key regions of the natural disaster risk management of China. The results also show that the CN-TFN is able to reflect practical conditions of the evaluation problem of natural disaster system and to provide more reliability information as compared to the existing evaluation methods. This is as a result of its comprehensive usage of various information of subjective and objective uncertainties in the evaluation process of natural disaster risk system and its expression by confidence intervals. Due to the simplicity and generalization, the CN-TFM is applicable to comprehensive risk grade evaluation of various natural disaster systems.     

Numerical simulation of fracture flow with a mixed-hybrid FEM stochastic discrete fracture network model

We introduce a discrete fracture network model of stationary Darcy flow in fractured rocks. We approximate the fractures by a network of planar circle disks, which is generated on the basis of statistical data obtained from field measurements. We then discretize this network into a mesh consisting of triangular elements placed in three-dimensional space. We use geometrical approximations in fracture planes, which allow for a significant simplification of the final triangular meshes. We consider two-dimensional Darcy flow in each fracture. In order to accurately simulate the channeling effect, we assign to each triangle an aperture defining its hydraulic permeability. For the discretization we use the lowest order Raviart-Thomas mixed finite element method. This method gives quite an accurate velocity field, which is computed directly and which satisfies the mass balance on each triangular element. We demonstrate the use of this method on a model problem with a known analytical solution and describe the generation and triangulation of the fracture network and the computation of fracture flow for a particular real situation.     

基于混合常应力−光滑应变单元的极限分析方法

极限分析是岩土工程稳定性评价的重要方法之一。传统的有限元极限分析方法,采用低阶三角形单元时需要引入速度间断面并采用特殊网格布局,或者采用高阶三角形单元等措施来克服体积锁定问题和提高数值精度。在光滑有限元法（smoothed finite element method,简称SFEM）的基础上,提出了一种基于新型混合常应力−光滑应变单元的极限分析方法（mixed constant stress-smoothed strain element limit analysis,简称MCSE-LA方法）。在服从关联流动法则和Mohr-Coulomb屈服准则的基础上,MCSE-LA方法最终将数值极限分析转化为以应力和极限荷载乘子为基本未知量的二阶锥规划（second order cone programming,简称SOCP）问题。MCSE-LA方法具有形式简单、优化变量相对较少和无需显式的写出塑性内能耗散函数的优点,并且根据凸锥优化的对偶理论,可以从对偶问题中获得速度场和塑性乘子等信息。此外,还采用基于最大塑性剪应变率的网格自适应加密算法,该算法在塑性区细化网格,显著提高了新数值极限分析方法的计算效率和精度。最后通过边坡稳定分析的结果对比,验证了MCSE-LA方法的计算精度和效率均高于传统的有限元极限分析方法。     

Seismic response of 3D topographical irregularities under incoming elastic waves from point sources

We study some effects that produce 3D topographical irregularities under incoming elastic waves from point sources using the indirect boundary element method. This technique is based on the representation of elastic waves in terms of single-layer boundary sources. In this way reflected and diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. The field emited from the point source is computed with the moment tensor for a shear dislocation and analytical expressions of the elastodynamic 3D Green's functions. In this way we can construct a double couple with variable orientation. We compare the technique with that of Bouchon, who used the discrete wave number method for a shear dislocation in a halfspace with a triangular source function. We subsequently apply our method to simulate the seismic response of a mountain of ellipsoidal geometry. The source function that we have considered is a triangular pulse and we show results corresponding to the velocity registered over the surface of the irregularity and that of the halfspace. This motion is presented by means of snapshots showing the evolution of the wavefields that are present in the problem.     

基于统一标定的势接触力计算

传统的离散元类数值方法存在处理角-角接触时接触力方向不确定、三维情形时接触计算极其复杂等问题。Munjiza提出的基于势的接触力计算方法很好地解决了这一问题,但同时存在势的物理意义不明确,计算出的接触力与物理直观不符、相同嵌入量计算的接触力大小不一致等问题。在已有研究的基础上,提出了一种基于统一标定的势接触力计算方法,该方法重新定义了势函数,即三角形单元内一点的势与该点到3条边的最短距离成正比。采用这种势函数定义可以对嵌入量进行统一表征。在相同的嵌入量下,计算的接触力大小是一致的,同时保留基于势的接触力计算方法所有优点。新定义的势函数,解决了原有势函数所存在的重大缺陷,只要重叠区域不变,计算的总接触力是不变的,接触力计算具有局部特性,而与划分的三角形单元的整体形状和大小无关,即新的接触力计算方法,具有更好的网格鲁棒性。该方法使得基于势的接触力计算方法建立在坚实的物理基础之上,意义重大。     

A numerical solution of an irrigation flowfield

Two-dimensional seepage from an array of triangular channels into permeable soil underlain at a finite depth by a drain has been solved by using the Baiocchi method. The use of this approach allowed a simple formulation of the seepage problem through extension of the solution region to an a priori known region. The new problem was treated numerically by using the S.O.R. method combined with projection. This technique proved to be very accurate when compared with analytical results and efficient from a computer programming standpoint.     

露采矿山宕口坡面表面积计算方法探讨

介绍了露采矿山整治中山体坡面表面积的两种计算方法:剖面法和小三角锁法。剖面法是把宕口坡表面分割成无数个小梯形进行计算表面积的一种方法。小三角锁法是通过对宕口坡表面数字化测量,根据测量数据建立数字高程模型,并在表面建立小三角锁,对所有小三角的面积累加得到面积的一种方法。小三角锁法大大提高了计算的准确度和效率。     

Optimal lower bound of passive earth pressure using finite elements and non-linear programming

The present study pertains to the finding of the lower bound solution, formulating it as a non-linear programming problem using the generalized method developed by Lysmer with certain variations to incorporate the non-linear no-yield condition constraints directly in the analysis. The method considers the family of plane stress fields having the property that all stresses vary linearly within each triangular element of some mesh which covers the soil mass under study. For this type of stress field it is possible to express all equilibrium conditions as a set of linear constraints and the no-yield as a set of non-linear constraints. The boundary condition constraints may be of linear equality or inequality type. By expressing some of the design variables in terms of the remaining variables the linear equality constraints are implicity satisfied. Such a technique minimizes the complexity of the problem by eliminating the equality constraints and reduces the dimensionality of the problem, saving much, computational effort. The optimal lower bound is isolated by formulating it as a non-linear programming (NLP) problem subjected to both linear and non-linear inequality constraints. The sequential unconstrained minimization technique using the extended penalty function method as suggested by Kavlie has been used to isolate the optimal lower bound. The method has successfully been applied to the passive earth pressure and bearing capacity problem. Numerical results are obtained and compared with Lysmer's solution to show the effectiveness of the present approach.     

岩土参数随机场离散的三角形单元局部平均法

将不确定性岩土参数建模为随机场而非传统意义上的随机变量,基于随机场的局部平均理论,提出了用于二维随机场离散的三角形单元局部平均法。通过面积坐标变换和高斯数值积分,给出了三角形单元局部平均随机场协方差矩阵的解析计算方法和数值计算方法。采用算例再现了所提方法的分析过程和有效性,并与传统二维随机场四边形单元离散法进行了对比。结果表明:提出的二维随机场三角形单元离散法能与有限元三角形单元离散法完美结合,随机场单元与有限元单元的对应关系清晰,易于随机有限元程序的编制;对于随机场单元的均值,传统四边形单元离散法与所提方法的计算结果相同;对于随机场单元的方差,传统四边形单元离散法计算结果偏小,所提方法显得更加科学、合理。     

三角模糊数型边坡支护方案评价模型及其应用

由于边坡支护方案评价时掌握的信息不充分,资料不全,加之问题本身的复杂性、不确定性以及人类思维的模糊性,方案的评价值往往难以精确表达,用三角模糊数的形式表示更合理。针对属性值用三角模糊数表示的边坡支护方案优选问题,提出了一种基于集对分析联系数的三角模糊数型方案综合评价模型。首先,根据集对分析中的联系数理论,利用三角模糊数的中值及上下确界所限定的取值区间,将三角模糊数转化为联系数;然后,根据联系数的加法和乘法定义,建立了联系数决策模型;按照综合决策值大小即可得到最优方案。该模型既考虑了三角模糊数的中值,又兼顾了三角模糊数的上下确界,即联系数的差异度,更具有客观性、合理性。实例应用表明,该模型算理清晰,计算简便,结论与实际相符。