三维FSM·DDM边界元法

对三维FSMDDM边界元数值模拟系统进行了基础研究,开发了单一介质的三维弹性问题的FSMDDM边界元数值模拟系统和多介质三维FSMDDM边界元数值模拟系统,通过验证和实际工程应用取得了满意的效果。     

Efficient solution of 3-D geomechanical problems by indirect BEM using iterative methods

The research herein primarily addresses to geomechanical problems of underground constructions in Mining and Civil Engineering. The problems are solved using the Indirect Boundary Element Method (IBEM). Although the geometry of the constructions themselves is usually very complicated, it will become much more complicated if we were to draw the existing joints. The computational problem therefore is how to deal with huge amount of equations and find out efficient methods of their formation and solution keeping in mind restraints of the computer memory and calculation time. Several approaches are used to enhance the performance of the Indirect Boundary Element Method. One of them deals with application of efficient equation solvers. It is shown that Krylov-type methods like CGS and GMRES with simple Jacoby preconditioning appear to be efficient and robust. In addition, adaptive integration on the boundary elements, together with diagonal dominance of equationsmake it possible to accelerate convergence of the iterative procedure. Some of the problems discussed allow a substantial reduction of matrix entries that leads to a very cheap iterative solution keeping reasonable accuracy of the results. © 1998 John Wiley & Sons, Ltd.     

Solving archaeological problems using techniques of soil magnetism

Archaeological investigations into the processes and factors controlling the magnetic enhancement of soils were initiated in the 1960s in conjunction with the application of magnetic prospection surveys. After a highly productive period of research, archaeological interest in soil magnetism experienced a downswing, which has just recently been reversed. Continuing research within the fields of rock magnetism and environmental magnetism has refined our understanding of the processes and products of magnetic enhancement and demonstrated the broad applicability of soil magnetic analyses. Using this information, together with recent archaeological studies, we update archaeologists on the variety of soil magnetic techniques and how they can be utilized to define, detail, and understand archaeological sites, features, and environments. A case study at the Cahokia Mounds Site in southwestern Illinois illustrates how soil magnetic techniques can be interfaced with other methods to understand archaeological landscapes and landscape change. © 1998 John Wiley & Sons, Inc.     

颗粒破碎的三维离散元模拟研究

采用PFC等离散元方法研究岩土材料的颗粒破碎已经成为热点。采用考虑局部应力集中的点荷载破碎准则,利用阿波罗填充和膨胀法保证破碎前、后颗粒之间的种群平衡,并引入尺寸因子来表征不同粒径的颗粒强度。在此基础上,开展了石英砂、钙质砂和萨克拉门托河砂3种不同破碎难易程度材料的数值试验,并与室内试验结果进行对比。结果表明:建立的三维颗粒破碎模型能够很好地描述破碎难易程度不同的颗粒材料的压缩特性;考虑应力集中效应的点荷载破碎准则比基于平均应力Mohr-Coulomb理论的颗粒破碎准则更能真实地反应颗粒材料的破碎现象。同时,所建立的模型能够揭示破碎对颗粒材料各向异性消散和级配曲线演化的影响规律。     

An artificial boundary approach for unbounded transient seepage problems

A new artificial boundary approach for transient seepage problems in unbounded domain is presented. The artificial boundary condition at the truncated boundary is derived from the analytical solutions for transient seepage problems in one dimension, including solutions, respectively, for flow in one‐dimensional infinite space and for radial flow in an infinite layer, and then it is tentatively applied for some two dimensional problems in addition to the one‐dimensional problems mentioned above. The boundary conditions derived relate the time‐dependent boundary flux with the time derivative of the hydraulic head at the truncated boundary, which makes the implementation much easier compared with the infinite element method. The accuracy and efficiency of the artificial boundary are validated by several numerical examples, which shows that the proposed boundary can give very good results for one‐dimensional transient seepage problems, as expected, whereas reasonable results can be also obtained for two‐dimensional problems, such as two‐dimensional axisymmetric flow and flow in an infinite plane. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite element simulations of 3D planar hydraulic fracture propagation using a coupled hydro-mechanical interface element

Two-dimensional hydraulic fracturing simulations using the cohesive zone model (CZM) can be readily found in the literature; however, to our knowledge, verified 3D cohesive zone modeling is not available. We present the development of a 3D fully coupled hydro-mechanical finite element method (FEM) model (with parallel computation framework) and its application to hydraulic fracturing. A special zero-thickness interface element based on the CZM is developed for modeling fracture propagation and fluid flow. A local traction-separation law with strain softening is used to capture tensile cracking. The model is verified by considering penny-shaped hydraulic fracture and plain strain Kristianovich‑Geertsma‑de Klerk hydraulic fracture (in 3D) in the viscosity- and toughness-dominated regimes. Good agreement between numerical results and analytical solutions has been achieved. The model is used to investigate the influence of rock and fluid properties on hydraulic fracturing. Lower stiffness tip cohesive elements tend to yield a larger elastic deformation around the fracture tips before the tensile strength is reached, generating a larger fracture length and lower fracture pressure compared with higher stiffness elements. It is found that the energy release rate has almost no influence on hydraulic fracturing in the viscosity-dominated regime because the energy spent in creating new fractures is too small when compared with the total input energy. For the toughness-dominated regime, the released energy during fracturing should be accurately captured; relatively large tensile strength should be used in order to match numerical results to the asymptotic analytical solutions. It requires smaller elements when compared with those used in the viscosity-dominated regime.     

Numerical simulation of drained triaxial test using 3D discrete element modeling

A discrete element modeling of granular material was carried out using a 3D spherical discrete model with a rolling resistance, in order to take into account the roughness of grains. The numerical model of Labenne sand was generated, and the desired porosity was obtained by a radius expansion method. Using numerical triaxial tests the micro-mechanical properties of the numerical material were calibrated in order to match the macroscopic response of the real material. Numerical simulations were carried out under the same conditions as the physical experiments (porosity, boundary conditions and loading). The pre-peak, peak and post-peak behavior of the numerical material was studied. The calibration procedure revealed that the peak stress of the sand sample does not only depend on local friction parameters but also on the rolling resistance. The larger the value of the applied rolling resistance, the higher the resulting stress peak. Furthermore, the deformational response depends strongly on local friction. The numerical results are quantitatively in agreement with the laboratory test results.     

Fully deformable discrete element analysis using a finite element approach

A method of analysis is presented for problems in which the deformation of the individual blocks play a significant role. The blocks are modeled as single quadrilateral elements and a constitutive model has been presented for computing the contact forces. In order to illustrate the influence of the deformation of the individual blocks two examples have been presented. The results of the analysis with deformable blocks and with rigid blocks are compared. These examples clearly demonstrate the importance of the deformation of blocks in the class of problems represented by these two examples.     

A 3D discrete FEM iterative algorithm for solving the water pipe cooling problems of massive concrete structures

Water pipe cooling has been widely used for the temperature control and crack prevention of massive concrete structures such as high dams. Because both under‐cooling and over‐cooling may reduce the efficiency of crack prevention, or even lead to great harm to structures, we need an accurate and robust numerical tool for the prediction of cooling effect. Here, a 3D discrete FEM Iterative Algorithm is introduced, which can simulate the concrete temperature gradient near the pipes, as well as the water temperature rising along the pipes. On the basis of the heat balance between water and concrete, the whole temperature field of the problem can be computed exactly within a few iteration steps. Providing the pipe meshing tool for building the FE model, this algorithm can take account of the water pipe distribution, the variation of water flow, water temperature, and other factors, while the traditional equivalent algorithm based on semi‐theoretical solutions can only solve problems with constant water flow and water temperature. The validation and convergence are proved by comparing the simulated results and analytical solutions of two standard second‐stage cooling problems. Then, a practical concrete block with different cooling schemes is analyzed and the influences of cooling factors are investigated. In the end, detailed guidance for pipe system optimization is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Site scale modeling of slow-moving landslides,a 3D viscoplastic finite element modeling approach

This paper presents an advanced 3D numerical methodology to reproduce the kinematics of slow active landslides, more precisely, to reproduce the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. For this purpose, finite element analyses are performed in a large area covering a long time-span (12 years), in order to exhibit different interacting slope movements. First, we perform a stability analysis using the shear strength reduction (SSR) technique with a Mohr-Coulomb failure criteria. It is done in order to compute factors of safety (FS) and to identify two different scenarios, the first one being stable (FS > 1) and the second one being unstable (FS < 1). In the studied test case, the Portalet landslide (Central Spanish Pyrenees), the first scenario corresponds to an initial stable configuration of the slope and the second one to an unstable excavated configuration. Second, taking the first scenario as an initial condition, a time-dependent analysis is performed using a coupled formulation to model solid skeleton and pore fluids interaction, and a simplified ground water model that takes into account daily rainfall intensity. In this case, a viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behavior and the delayed creep deformation due to the excavation. The fluidity parameter is calibrated to reproduce displacements measured by the monitoring systems. Our results demonstrate that 3D analyses are preferable to 2D ones for reproducing in a more realistic way the slide behavior. After calibration, the proposed model is able to simulate successfully short- and medium-term predictions during stages of primary and secondary creep.     

Elastoplastic analysis of jointed rocks using a coupled finite element and boundary element method

This paper presents a coupled, elastoplastic, finite element and boundary element method for the two-dimensional, non-linear analysis of anisotropic jointed rock. The non-linear and anisotropic behaviour of a jointed rock mass is simulated by representing the mass as an equivalent anisotropic, elastoplastic continuum, so that the influence of the jointing system is ‘smeared’ across the continuum, i.e. the individual joints are not modelled as discrete entities. Numerical examples have been solved to verify the capability, accuracy and efficiency of the present technique. The proposed technique has also been applied to the analysis of tunnel excavation problems in plane strain. The effects of anisotropy and non-linearity of the jointed rock mass during excavation have been investigated in some detail.     

Buckling analysis of 3D layered structures using a Cosserat continuum approach

This paper concerns the buckling analysis of three-dimensional (3D) layered continua using the Cosserat approach. The finite element method (FEM) is used to implement the adopted Cosserat formulation. As a result, the interfaces between the layers need not be explicitly modelled. Instead, the internal characteristic length, i.e., the layer thickness and also the interaction conditions between the layers are incorporated into the governing equations of the solution. This paper introduces a new 3D geometric stiffness matrix based on the principle of virtual work. The proposed geometric stiffness matrix is applied to the linear buckling analysis of a number of benchmark problems with various geometries, boundary conditions, and interaction conditions between the layers. In all cases, the FEM Cosserat solution exhibits a high level of consistency with the analytical solution.     

重磁2.5D/3D互相约束反演技术在三维地质填图中的应用研究

三维地质填图为我国启动的新一轮地质调查项目,为提高重磁资料在三维地质填图中的应用效果,笔者提出重磁资料2.5D/3D相互约束重磁反演技术方案:利用重磁资料2.5D剖面反演结果、3D物性反演结果作为彼此反演约束条件,并通过了理论模型试验。试验结果表明,该技术方案使反演结果中物性参数、空间位置更接近理论模型体。通过对本溪—临江地区思山岭铁矿磁异常及酸性岩体重力异常进行反演实践——估算磁异常铁矿资源量、研究酸性侵入岩深部展布形态的效果良好,可为大面积三维填图提供有效途径。     

非线性波浪变形计算的三维边界元方法

采用0-1混合型边界元剖分计算域边界面,以提高边界元方法模拟三维波场波浪变形问题的数值计算精度。采用Spongelayer阻尼消波和Sommerfeld放射条件相匹配的边界处理方式,以消除计算域内的反射波。借助时域内的波面位置追踪,实现了三阶Stokes波在数值波动水槽中的波浪变形计算,其结果吻合于理论值。     

A finite element approach to solve contact problems in geotechnical engineering

This paper presents a finite element approach to solve geotechnical problems with interfaces. The behaviours of interfaces obey the Mohr–Coulomb law. The FEM formulae are constructed by means of the principle of virtual displacement with contact boundary. To meet displacement compatibility conditions on contact boundary, independent degrees of freedom are taken as unknowns in FEM equations, instead of conventional nodal displacements. Examples on pressure distribution beneath a rigid strip footing, lateral earth pressure on retaining walls, behaviours of axially loaded bored piles, a shield‐driven metro tunnel, and interaction of a sliding slope with the tunnels going through it are solved with this method. The results show good agreement with analytical solutions or with in situ test results. Copyright © 2005 John Wiley & Sons, Ltd.     

A sequential partly iterative approach for multicomponent reactive transport with CORE2D

CORE^2D V4 is a finite element code for modeling partly or fully saturated water flow, heat transport, and multicomponent reactive solute transport under both local chemical equilibrium and kinetic conditions. It can handle coupled microbial processes and geochemical reactions such as acid–base, aqueous complexation, redox, mineral dissolution/precipitation, gas dissolution/exsolution, ion exchange, sorption via linear and nonlinear isotherms, and sorption via surface complexation. Hydraulic parameters may change due to mineral precipitation/dissolution reactions. Coupled transport and chemical equations are solved by using sequential iterative approaches. A sequential partly iterative approach (SPIA) is presented which improves the accuracy of the traditional sequential non-iterative approach (SNIA) and is more efficient than the general sequential iterative approach (SIA). While SNIA leads to a substantial saving of computing time, it introduces numerical errors which are especially large for cation exchange reactions. SPIA improves the efficiency of SIA because the iteration between transport and chemical equations is only performed in nodes with a large mass transfer between solid and liquid phases. The efficiency and accuracy of SPIA are compared to those of SIA and SNIA using synthetic examples and a case study of reactive transport through the Llobregat Delta aquitard in Spain. SPIA is found to be as accurate as SIA while requiring significantly less CPU time. In addition, SPIA is much more accurate than SNIA with only a minor increase in computing time. A further enhancement of the efficiency of SPIA is achieved by improving the efficiency of the Newton–Raphson method used for solving chemical equations. Such an improvement is obtained by working with increments of log concentrations and ignoring the terms of the Jacobian matrix containing derivatives of activity coefficients. A proof is given for the symmetry and non-singularity of the Jacobian matrix. Numerical analyses performed with synthetic examples confirm that these modifications improve the efficiency and convergence of the iterative algorithm. Changbing Yang is now at The University of Texas at Austin, USA.     

破碎岩体压实特性的三维离散元数值计算方法研究

研究破碎岩体的压实特性是矿井地下工程的基础工作之一,由于破碎岩体所处环境的隐蔽性与危险性,常采用实验室测试和数值计算的研究方法。提出了一种三维破碎岩体模型构建方法,即在3D Voronoi建立完整岩体数值模型的基础上,通过预定孔隙率,随机删除完整岩体中的块体反演破碎岩体结构,测定破碎岩体的压实特性。该方法可较真实地反映破碎岩体的块度特征、碎胀特性与压实特性,与现有研究方法有较高的吻合度,为矿山地下工程的安全控制提供了新的有效的研究方法。     

2D inversion of the magnetotelluric data from Mahallat geothermal field in Iran using finite element approach

The natural-field magnetotelluric (MT) method has proven very useful for mapping the geothermal fields as resistivity sections. The depth of investigation of the MT method is sufficiently large to penetrate deep into the upper crust. MT soundings along two transects across Mahallat geothermal field in Iran were carried out to determine the crustal structure in the region. The selected MT profiles in the region cross over the hydrothermally altered zones and different geological structures. Data were acquired along two profiles crossing the Mahallat hot springs with a total of 28 MT stations in a frequency range of 8,000 to 0.008 Hz. Spacing between stations was kept 500 m for a good resolution. We have used the code MT2DInvMATLAB for inversion using the method of finite elements for forward modeling. Apparent resistivity and phase data of transverse electric (TE), transverse magnetic (TM), and TE + TM modes along each profile were modeled. The geothermal fluid reservoir is resolved at 1,000 to 3,000 m depth and the geothermal resource is estimated to be located at 7,000 m or deeper.     

Modeling crack propagation path of anisotropic rocks using boundary element method

In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagation and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed‐mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Copyright © 2008 John Wiley & Sons, Ltd.