二维与三维问题的无界边界单元

本文通过对二维与三维问题的无界边界元的分析,指出二维问题中原有边界元上插值函数的缺点和不合理性,提出一种新型的无界边界元插值形式,对无界单元上的发散积分作了一些讨论。     

三维离散元法中的数据结构

论述了动态松弛三维离散元法的数据结构,提出了地下水、锚杆、边界力等数据的组织方法,为三维离散元程序的编写提供了依据。     

三度体重力矢量的有限单元法正演计算

以重力位在场源内部满足泊松方程为依据,以重力矢量满足第三类边界条件为切入点,推导了与三度体重力矢量满足的边值问题相对应的变分问题,进而利用有限单元法实现了对变分问题的求解.立方体模型试验结果表明:文中提出的新的系数矩阵存储方式较之传统方式能够更有效地节约存储空间,且为利用预条件共轭梯度技术更加快速地求解线性方程组提供了保障;重力矢量的计算精度与边界长度及单元网格的边长息息相关,其计算效率则主要取决于所要计算的节点总数和大型稀疏线性方程组求解算法的优劣;一般情况下,当单元的边长小于场源体边长的1/10、边界长度大于场源体长度的7.5倍时,能够获得理想的结果.     

二维扩散输移问题的一种新的有限体积算法

从空间离散格式上比较了有限元格式与有限体积格式的异同处,证明有限体积格式是有限元格式的一种特例,有限元格式也是一种守恒型格式。并根据有限元格式与有限体积格式的比较,提出了一种新的高精度、稳定性好的有限体积格式。通过模拟非定常的纯对流方程和二维对流扩散方程,证明了有限体积格式的优点。     

边界元法

前言边界积分方程法是一个比较新的数值应力分析方法,它具有一些优点,特别是在三维分析中这些优点显得更为突出。近几年来,边界元法在我国得到了初步的应用和发展,为了进一步推广,有的单位已举办了一些小型的训练班和座谈会。但是有关边界积分方程法在岩土力学中,特别是在岩土工程中的应用,所见报导不多。我们这个讲座除介绍基本理论和方法外,将着重介绍边界积分方程法用于工程的几个具体实例。     

边界元法(续)

二、积分方程的数值解积分方程(16)一般来说是不能得到解析解的,因此必须求助于数值方法,如象有限元法那样。但这里仅对区域表面进行离散化,而不必对整个物体这样做。     

分层伸展叠加变形二维离散元模拟

分层伸展叠加变形是裂陷盆地演化过程中普遍存在的一种复杂的构造地质现象,是裂陷盆地构造地质研究和油气成藏研究的热点。目前有关分层伸展叠加变形的研究仍处于几何学描述阶段,缺乏对其动力学演化过程的解析。本文采用构造地质研究领域中新兴的模拟方法——离散元,通过构建离散元模型模拟研究分层伸展叠加变形的断裂发育过程。结果表明:(1)上下两套伸展断裂系统之间是否存在韧性滑脱层(泥岩层或膏岩层)是分层叠加变形的主控物质因素;(2)中间韧性地层隔断了下部断层向上的继承性发育,从而形成上下两套断裂系统分层伸展,垂向上叠加的变形特征;(3)受韧性地层侧向牵引作用影响,上部断裂系统的横向展布明显扩大;(4)在相同岩性条件下,降低伸展速度会显著降低两侧断层发育的差异性,但不会改变垂向上分层伸展的变形特征。研究解析了分层伸展叠加变形的动力学演化过程,有助于进一步认识分层伸展叠加变形机制及其对油气运聚成藏的控制作用。     

三维结线动力无穷元及其特性研究

研究了文献[2]提出的结线动力无穷元的特性及其参数确定。由于这种单元将衰减函数定义在无穷向结线上,从而不仅可同时考虑多种波动形式,而且对于不同介质,不必通过改变节点坐标分布的形式实现介质间的过渡,而直接以结线上参数的改变来实现,因而可以更好地模拟无穷波动场,使用起来更方便．更容易推广。     

Three Dimensional Plate Kinematics in Romania

Since 1996 the Collaborative Research Center (CRC) 461``Strong Earthquakes: A Challenge for Geosciences andCivil Engineering' is funded by the DeutscheForschungsgemeinschaft (German Research Foundation) asa German contribution to the UN initiative``International Decade of Disaster Reduction' (IDNDR).A central project of this CRC is the subproject B1``Three Dimensional Plate Kinematics in Romania', whichis installed to detect borders of the geotectonicplates in Romania, to quantify their three-dimensionalmovement rates and to determine in detail thedeformation of each geotectonic unit in the VranceaRegion as a contribution to the research of thetectonic cause of the intermediate depth earthquakesin this region. These aims shall be accomplished byrepeated GPS-measurements in a network which has beeninstalled in 1997 and measured until now in 1997 and1998. This network is consisting of 26 stations,covering an area of 250 km × 380 km with the Vranceaarea as the centre. The geodetic frame is given by sixstations of the CEGRN (Central European GeodynamicRegional Network) of the CEI-initiative CERGOP, forwhich coordinates and velocities in a global tectonicscenario are available. The paper gives informationabout the tectonic background, the geodetic network,the GPS-measurements and the achieved accuracies andfirst results of deformation analyses. It is shownthat until now no significant deformations could bedetected. But the extension of the GPS network andremeasurements with time intervals of two years willallow statements about the recent crustal movementswithin the lifetime of the CRC 461.     

不同桩型情况下的桩基差异沉降三维有限元分析

近年来，以控制差异沉降为目标的桩基础变刚度调平方法在我国岩土工程界正受到越来越多的重视和应用。该文用三维有限元方法研究在同一基础下采用不同桩型情况时各桩沉降和刚度的关系，探讨通过改变桩型来调整刚度，最后控制差异沉降的方法，并通过一个工程实例阐述运用不同桩型进行设计的概念和方法，达到减小差异沉降。     

Boundary Element based Discontinuous Deformation Analysis

A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd.     

GPS三维监测试验网的两期成果分析

通过介绍GPS三维监测试验网的第二期观测成果,对两期观测成果进行了比较分析。结果表明,两期成果的尺度没有系统偏差,而GPS基线向量的坐标轴指向存在系统偏差。只要统一各期成果的坐标轴指向后GPS将有望以毫米级精度进行三维形变监测。     

SKUA基岩三维地质建模

以天津平原地区某个工区为例,介绍了三维地质模型建模的数据准备,SKUA建模的流程以及模型分析,共使用28条断层,6个地层的地质构造图,3条地质剖面构建了地质结构模型,表明SKUA建模是高效的和实用的。     

二维渗流有限元法在潮州水库浸没分析中的应用

通过对潮州市意东堤、北堤及东厢堤段的二维渗流计算,了解在潮州供水枢纽建成前后的潮州市城区地下水位变化情况,从而判断在枢纽建成后,在长期处于正常蓄水位的情况下,潮州市城区的浸没情况。同时复核在采取工程措施的情况下,对堤后浸没区的影响,确保在枢纽建成后,潮州市不产生浸没。     

Migration of Magnetotelluric Data in Two——Dimensional Model

Since the wave equation of magnetotelluric(MT)field is similar to the one of seis-mic,the migration techniques used in seismic can be applied to MT data.In thispaper we make use of the principle of reflector mapping(i.e.U/D imaging principle)to image MT data.That is,the MT wavefield observed on the surface of the earth canbe resolved into upgoing and downgoing waves,the waves are extrapolated downward bythe phase-shift method or the phase-shift plus interpolation(PSPI)method.Conductivi-ty interfaces of the medium could be found by using the time coincidence of the upgoingand downgoing waves.Theoretical calculations show that the migration technique of MTdata presented here is very effective.It can not only enhance the lateral resolution of MTdata,but also obtain the visual image of subsurface interfaces.As compared with theconventional 2-D inversion,this procedure is more simple in calculation and can be easi-ly put into practice on a personal computer and is able to obtain the MT depth section,which is     

苏北沿海三市三维地下水流数值模拟

20世纪80年代以来,苏北沿海三市（连云港、盐城和南通）地下水过量开采导致部分地区水位持续大幅下降,引发地面沉降、水质咸化等问题。本文基于研究区大量地质、水文地质资料,以及2005—2013年地下水开采和水位观测数据,经过模型识别、验证,建立了能刻画苏北沿海三市地下水流场演化的三维非均质各向异性地下水流数值模型,模拟了研究区2005—2013年开采量由约2.03亿m³逐渐减少到1.56亿m³条件下深部地下水系统流场（包括第Ⅱ、Ⅲ、Ⅳ承压含水层）的动态演化过程。结果表明,各观测孔的水位实测值与模拟值吻合良好,水位及水均衡模拟结果显示各含水层水位下降得到有效控制,但是该区域地下水仍处于超采状态。为进一步控制水位下降,防治咸水扩散、地面沉降等地质环境灾害,模型预测并比较了现状开采和限制开采2种不同开采方案下2014—2020年苏北沿海三市深层地下水流场的变化趋势。预报结果显示,2种开采方案下,各含水层水位下降速率均减小,尤其是限制开采方案下,第Ⅱ、Ⅲ和Ⅳ承压含水层的区域水位平均下降速率依次为0.15、0.16和0.15 m/a,分别比模拟期减小了46.43%、65.21%和48.28%,在开采强度降至最低的条件下储存量变化速率仍为负值,说明需要进一步加强地下水限制开采的力度。     

Multi-Region Boundary Element Analysis for Coupled Thermal-Fracturing Processes in Geomaterials

This paper describes a boundary element code development on coupled thermal–mechanical processes of rock fracture propagation. The code development was based on the fracture mechanics code FRACOD that has previously been developed by Shen and Stephansson (Int J Eng Fracture Mech 47:177–189, 1993) and FRACOM (A fracture propagation code—FRACOD, User’s manual. FRACOM Ltd. 2002) and simulates complex fracture propagation in rocks governed by both tensile and shear mechanisms. For the coupled thermal-fracturing analysis, an indirect boundary element method, namely the fictitious heat source method, was implemented in FRACOD to simulate the temperature change and thermal stresses in rocks. This indirect method is particularly suitable for the thermal-fracturing coupling in FRACOD where the displacement discontinuity method is used for mechanical simulation. The coupled code was also extended to simulate multiple region problems in which rock mass, concrete linings and insulation layers with different thermal and mechanical properties were present. Both verification and application cases were presented where a point heat source in a 2D infinite medium and a pilot LNG underground cavern were solved and studied using the coupled code. Good agreement was observed between the simulation results, analytical solutions and in situ measurements which validates an applicability of the developed coupled code.     

Three Dimensional Numerical Analysis of Underground Bifurcated Tunnel

As far as the bifurcated tunnel of underground engineering is concerned, it is usually used in the water conveyance system. Due to the complexity of underground rock masses and concrete lining, researches on mechanical characteristic and stability of the bifurcation tunnel have attracted more and more attention in the geotechnical field. In order to understand bifurcated tunnel in detail, three-dimensional (3D) numerical method is applied to solve the above key subjects by simulating a practical project. Furthermore, sub-model technology is applied to analyze the intersection position, corresponding deformation and stress results in the practical condition. Meanwhile, 3D excavation and support calculation under four conditions have been simulated based on 3D self-compiled code and Ansys software. In addition, the paper plays emphasis on the stress, displacement analysis considering different stress releasing ratio instead of rheological analysis, and the results corresponding with the fact indicate the feasibility of 3D elasto–visco–plastic code.