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

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

Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems

A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

重力坝应力-渗流相互作用的无单元耦合分析

应用两场交叉迭代解耦技术,建立了应力-渗流相互作用的无单元耦合分析计算模型。详细分析了东方红重力坝在不计应力-渗流耦合和考虑应力-渗流耦合两种情况下,正常蓄水时坝体的位移、应力和渗流位势分布的变化规律。分析表明,考虑耦合作用后,坝体的位移略小,但分布规律一致;坝体的水平正应力分布不同,垂直正应力变化较大,但分布规律一致;耦合作用条件下渗流对重力坝坝体应力的影响是不利的。计算成果符合一般规律。     

间断伽辽金法(DGM)求解弹性地基梁问题

间断伽辽金法使用节点位移一类未知数作为测试函数,削弱了内部单元边界上的一阶及n阶导数的连续性,大大降低了构造形函数的难度,特别适合控制方程为高阶微分方程问题的求解。基于间断伽辽金法的基本原理,推导了弹性地基梁四阶微分控制方程的积分“弱”形式,编制了计算程序,进行了数值计算和收敛性分析。计算结果表明：用间断伽辽金法求解弹性地基梁问题是十分有效率的。     

改进的无网格法及其在冻土区桩基温度场中的应用

采用带权正交基函数对传统的无网格法中的基函数进行了改进,避免了计算过程中可能出现的矩阵不可逆情况,而且编程容易实现且计算效率高。并将其应用到冻土区桩基包括了热传导和相变潜热的温度场中,同时考虑了混凝土水化热释放对桩周冻土的影响,将计算结果和有限元计算结果及现场实测数据进行了比较分析,从趋势和最大数值看都反映了实测曲线的趋势,又把不同深度处桩侧和不同桩径处温度随时间的变化规律计算结果和有限元计算结果做了对比,都验证了该方法的可行性和优越性。     

核废料贮存裂隙岩体中THM耦合过程的有限元分析

本文首先给出了核废料贮存裂隙岩体中THM耦合作用的数学模型,然后用加权残数法,建立了以结点位移、水压力和温度为求解量的有限元计算格式,提出了数值模拟渗流、附加应力和热载耦合作用的有限元分析方法,最后给出了一个BMT1问题的算例.     

Fully coupled hydraulic fracture simulation using the improved element partition method

The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation.     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

测定非饱和土土-水特征曲线的溢出水校正称量法

在利用压力板仪测定非饱和土的土-水特征曲线时,常用的取样称量法、溢出水称量法分别存在实测试样质量与实际试样质量不相等、实测排水量与实际排水量不相等等问题,从而造成试验结果失真。为此,溢出水校正称量法被引入用来消除取样称量法及溢出水称量法中存在的问题。通过取样称量法、溢出水称量法、溢出水校正称量法的对比研究发现,在同一基质吸力作用下,由取样称量法、溢出水称量法确定的试样饱和度均小于由溢出水校正称量法确定的试样饱和度,并对该现象进行了解释。随后,对取样称量法、溢出水称量法、溢出水校正称量法的适用范围进行了探讨。     

Global and local sensitivity analyses for coupled thermo–hydro–mechanical problems

This paper focuses on the sensitivity analysis for coupled thermo–hydro–mechanical problems employing both local and global sensitivity methods. A derivative‐based method is used in the local sensitivity approach, whereas the random balance designs method is used for the global sensitivity analysis. The main goal is to investigate the effect of uncertainties in the constitutive parameters on the results from nonlinear coupled thermo–hydro–mechanical analyses of unsaturated soil behavior whose modeling generally involves large sets of constitutive relations. Knowing the parameter sensitivity allows to qualitatively assess the validity of the results obtained by computational simulations of high‐risk situations, for example, emerging nuclear waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.     

全吸力范围内土-水特征曲线的快速测定方法

传统的宽（全）吸力范围内的土-水特征曲线（SWCC）测量方法动辄耗时数月。为了更好地使非饱和土力学理论指导工程实践,SWCC的快速测定显得尤为关键。为开发宽（全）吸力范围内SWCC的快速测定方法,联合压力板法、并行滤纸法、露点水势仪法进行了两种不同土体的宽吸力范围内SWCC的测定试验。试验结果表明：（1）并行滤纸法与串行滤纸法的试验结果能较好地吻合,故可采用并行滤纸法代替串行滤纸法以缩短测量时间;（2）当测量的试验数据满足SWCC的大致形状时（或满足工程设计需求时）,可适当减少试验数量,以提高测定效率,联合测定方法测量一条完整的SWCC仅需要10～12个点即可完成;（3）联合测定方法可将SWCC测量时间由数月缩短至7～10 d。所提出的联合测定方法实现了宽（全）吸力范围内SWCC的快速、精准测定,有望使得SWCC的测量成为常规的土工试验。     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

河网海湾水动力联网数学模型

在一维显式、二维隐式水动力联网数学模型的基础上,建立了一种一维、二维全隐河网海湾水动力联网数学模型。该一维河网模型采用Preissmann四点隐式格式,用节点水位控制法进行数值计算,二维海湾模型采用改进型双向隐式(DSI)法进行数值求解。在河口一维、二维连接处,水力因子通过接口断面法传递,无需重叠一个一维河段来传递水力因子,避免了口门处二维网格需取较小尺度的问题。模型在珠江三角洲河网及横门、洪奇沥口门海域做了检验,验证结果表明,建立的河网海湾水动力联网数学模型是可行的,可以用于河网、河口治理工程的数值研究。     

钻井船插桩对邻近桩影响的耦合欧拉-拉格朗日有限元方法研究

依据模拟钻井船在黏土层中插桩对邻近桩影响的离心模型试验结果,研究了通过耦合欧拉-拉格朗日（CEL）有限元计算并结合非线性地基梁有限元计算,分析钻井船插桩对邻近桩影响问题的可行性。CEL有限元方法将产生大变形的土体设为欧拉体,采用欧拉有限元方法计算该区域的变形响应,计算过程中,欧拉体的空间网格形状、大小位置保持不变,物质可在网格之间运动;其他土体设为拉格朗日体,采用拉格朗日有限元方法计算变形响应,计算过程中,物质的运动和网格的变形保持一致。运用罚函数方法实现欧拉体与拉格朗日体的耦合。通过CEL有限元计算,可以确定钻井船插桩导致的邻近桩桩身水平位移。进一步通过非线性地基梁有限元模型计算确定桩身弯矩。计算结果表明,利用CEL有限元方法并结合非线性地基梁有限元方法计算出的桩身位移和弯矩沿桩长的变化与离心模型试验结果基本一致。说明采用CEL有限元方法并结合非线性地基梁有限元方法分析黏土层中插桩对邻近桩的影响问题是可行的。CEL有限元模型中欧拉土体范围的设置对计算结果有明显影响。研究表明,若插桩深度小于0.75倍桩靴直径,可将欧拉土体范围设置成1.00～1.25倍桩靴直径;若插桩深度大于0.75倍桩靴直径,将欧拉土体范围取为插桩深度以下0.5倍桩靴直径是恰当的。     

A coupled time integration algorithm for discontinuous deformation analysis using the numerical manifold method

To improve the computational efficiency of the numerical manifold method for discontinuous deformation simulations, a spatial-domain coupled explicit-implicit time integration algorithm is proposed. A subdomain partition algorithm based on a super manifold element is developed for the numerical manifold method to simulate dynamic motions of blocky rock mass. In different subdomains, explicit or implicit time integration method is employed respectively based on its contact and motion status. These subdomains interact through assembling the corresponding explicit or implicit time integration-based matrices of different rock blocks. The computational efficiency of the discontinuity system under dynamic loading is improved by partially diagonalizing the global matrices. Two verification examples of a sliding block along an inclined plane under a horizontal acceleration excitation and a multiblock system acted on by dynamic forces are studied to examine the accuracy of the proposed numerical method, respectively. A highly fractured rock mass situated on an inclined slope subjected to seismic excitations is then studied to show the computational efficiency of the developed algorithm. The simulated results are in good agreement with those from the versions using purely implicit or explicit time integration algorithm for the numerical manifold method. The computational efficiency is shown to be higher using the proposed algorithm, which demonstrates its potential for application in dynamic analysis of highly fractured rock masses.     

A meshless method for axisymmetric problems in continuously nonhomogeneous saturated porous media

A meshless method based on the local Petrov–Galerkin approach is proposed to analyze 3-d axisymmetric problems in porous functionally graded materials. Constitutive equations for porous materials possess a coupling between mechanical displacements for solid and fluid phases. The work is based on the u–u formulation and the incognita fields of the coupled analysis in focus are the solid skeleton displacements and the fluid displacements. Independent spatial discretization is considered for each phase of the model, rendering a more flexible and efficient methodology. Both displacements are approximated by the moving least-squares (MLS) scheme. The paper presents in the first time a general meshless method for the numerical analysis of axisymmetric problems in continuously nonhomogeneous saturated porous media. Numerical results are given for boreholes in continuously nonhomogeneous porous medium with prescribed misfit and exponential variation of material parameters in the excavation zone.     

一种有限元转化为颗粒离散元的方法及其应用研究

为了充分发挥有限元与颗粒离散元各自的优势,提出了一种由有限元转化为颗粒流的方法。数值模型首先用较粗的有限元网格进行离散,并在单元上引入连续介质本构模型。力学计算开始后,实时跟踪各单元的应力状态。一旦某单元的应力满足Mohr-Coulomb准则或最大拉应力准则,删除该单元,同时创建具有一定数目、随机分布且微嵌入的颗粒簇。其后,该单元所在区域的非连续变形及失稳断裂由颗粒簇演化获得。各颗粒的质量、材料参数、速度、位移、接触力等信息根据插值从有限元单元中继承。为了实现有限元与颗粒流接触面的耦合计算,引入了点-棱（二维）及点-面（三维）接触模型,通过法向及切向弹簧实现接触力的计算。颗粒球与有限元板的碰撞分析、单轴压缩、岩石切割等案例展示了上述方法的精确性及合理性。