Three‐dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method

Failure in geotechnical engineering is often related to tension‐induced cracking in geomaterials. In this paper, a coupled meshless method and FEM is developed to analyze the problem of three‐dimensional cracking. The radial point interpolation method (RPIM) is used to model cracks in the smeared crack framework with an isotropic damage model. The identification of the meshless region is based on the stress state computed by FEM, and the adaptive coupling of RPIM and FEM is achieved by a direct algorithm. Mesh‐bias dependency, which poses difficulties in FEM‐based cracking simulations, is circumvented by a crack tracking algorithm. The performance of our scheme is demonstrated by two numerical examples, that is, the four‐point bending test on concrete beam and the surface cracks caused by tunnel excavation. Copyright © 2014 John Wiley & Sons, Ltd.     

复杂边界及实际地形上溃坝洪水流动过程模拟

建立了基于无结构三角网格下采用有限体积法求解的二维水动力学模型,用于模拟溃坝洪水在复杂边界及实际地形上的流动过程。该模型采用Roe格式的近似Riemann解计算界面水流通量,结合空间方向的TVD-MUSCL格式及时间方向的预测-校正格式,可使模型在时空方向具有二阶计算精度。模型中引入最小水深概念,提出了有效的干湿界面处理方法。模拟了理想条件下溃坝水流过程,研究不同最小水深取值对干河床上洪水演进的影响,并用两组简单溃坝水流的水槽试验资料对模型进行验证。采用该模型模拟了实际溃坝洪水的流动过程,所得计算结果与实测资料及已有模型计算结果较为符合。     

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.     

复杂区域自适应三角形网格全自动生成方法

为了使用自适应有限单元法求解岩土工程问题，本文在行被法的基础上，提出了一种适用于多介质复杂区域的三角形网格全自动生成方法。文中给出三个工程算例．     

通用三维GIS场数据模型研究与实践

随着三维地理信息系统(geographic information system,GIS)应用的发展,空中电磁信号场、空气与水体污染场、地下地质属性场等连续、非匀质的三维空间属性场数据进行建模、可视化与分析计算对GIS基础软件提出了新的挑战。已被广泛应用的TIN(triangulated irregular network)与栅格数据属于表面模型,通常用来表达二维场数据,其中由TIN组成的封闭表面可用于表达三维体,但两者都无法描述非匀质、连续变化的三维属性场。探讨了在通用GIS基础软件中基于TIN和Grid进行升维表达,实现不规则四面体网格(tetrahedralized irregular mesh,TIM)和体元栅格(voxel grid)两大类三维空间场数据模型,分析了两类数据模型的特点,探讨了三维场数据在可视化和分析计算方面的关键技术,实现了基于TIM与体元栅格表达和应用三维GIS场数据,完善了GIS三维数据模型体系,实现了空/天、地表、地下全空间的一体化表达、模拟与分析计算。     

A fluidized landslide on a natural slope by artificial rainfall

An experiment to induce a fluidized landslide by artificial rainfall was conducted on a natural slope at Mt. Kaba-san in the village of Yamato, Ibaraki Prefecture, Japan. The experimental slope was 30 m long, 5 m wide, and the average slope gradient was 33°. A landslide initiated 24,627.5 s (410 m/27.5 s) after the start of sprinkling at a rainfall intensity of 78 mm/h. The landslide mass was 14 m long and 1.2 m deep (at maximum). It first slid, then fluidized, and changed into a debris flow. The travel distance was up to 50 m in 17s. The apparent friction angle of the fluidized landslide was 16.7°. Formation of the sliding surface was detected by soil-strain probes. Motion of the surface of the failed landslide mass was determined by stereo photogrammetry.     

Centrifuge modeling of load-deformation behavior of rocking shallow foundations

Shallow foundations supporting building structures might be loaded well into their nonlinear range during intense earthquake loading. The nonlinearity of the soil may act as an energy dissipation mechanism, potentially reducing shaking demands exerted on the building. This nonlinearity, however, may result in permanent deformations that also cause damage to the building. Five series of tests on a large centrifuge, including 40 models of shear wall footings, were performed to study the nonlinear load-deformation characteristics during cyclic and earthquake loading. Footing dimensions, depth of embedment, wall weight, initial static vertical factor of safety, soil density, and soil type (dry sand and saturated clay) were systematically varied. The moment capacity was not observed to degrade with cycling, but due to the deformed shape of the footing–soil interface and uplift associated with large rotations, stiffness degradation was observed. Permanent deformations beneath the footing continue to accumulate with the number of cycles of loading, though the rate of accumulation of settlement decreases as the footing embeds itself.     

Load testing of anchors for wire mesh and cable net rockfall slope protection systems

This study presents results of field tests conducted on anchors used to support wire mesh and cable net rockfall protection systems. The load transfer and failure characteristics of these anchors are different from those used in most civil applications in that loads are often applied transversely to the top of tendon rather than axially. The study included vertical as well as horizontal series of tests conducted on some anchors widely used in wire mesh and cable net rockfall protect systems. It was found that the deformation characteristics of these anchors under vertical loading are nonlinear. They are approximated by a hyperbolic formulation and used to calculate the ultimate capacity. Top-downward progressive cracking of the grout was observed during loading and influences the deformation characteristics of these anchors under horizontal loading. The anchors deflected excessively before they could attain their ultimate capacity in the horizontal direction. Based on the field tests, it appears that the deformation under horizontal loading in the systems can be limited by using an enlarged grout zone at the top.     

A review of landslide hazards in Japan and assessment of their susceptibility using an analytical hierarchic process (AHP) method

In spite of its small size, Japan suffers many landslide disasters due to intense rainfall and earthquakes. This article describes the distribution and topography of these landslides, and a new method of evaluating the susceptibility, the analytical hierarchic process (AHP). The method assigns scores to each factor of micro-topography of landslide-prone areas identified in aerial photographs, and assesses the susceptibility of landslide from the total score. In addition, a method of simulating sliding mass runout is briefly presented for the designating sediment-related disaster warning areas.     

基于动力可靠性分析的滑移隔震体系的优化设计

本文探讨了滑移隔震体系在零均值高斯白噪声地震作用下的优化设计问题,首先利用编译等效线性化方法并按首次超超损坏泊松过程模型的双壁问题得到质点的滑移可靠性函数,然后采用拉格朗日乘子法进行了该滑移隔震体系基于动力可靠性分析的优化设计。