Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

华北克拉通破坏与岩石圈减薄

古太古代（约4.0 Ga）时地球上可能只有一个超级大陆, 它的岩石圈厚度高达400 km。在早元古代,这个超级大陆减薄、裂解成十几块,每块中心是太古宙岩石,边缘是元古宙岩石,且各块厚度不等（150~350km）。从元古宙之后这些被称之为稳定克拉通的大陆岩石圈就一直漂游在地幔软流圈之上。中国华北地块就是这些克拉通之一,与众不同的是它在中生代时遭受了第二次破坏,岩石圈厚度从古生代时的180~200 km 减少到现今的80~100 km。本文作者从流变学的视角出发,围绕华北克拉通破坏和岩石圈减薄这一核心问题,从     

PSO-RBFNN模型及其在岩土工程非线性时间序列预测中的应用

岩土工程受力变形演化是一个典型的非线性问题,其演化的高度非线性和复杂性,很难用简单的力学、数学模型描述,但可用粒子群优化径向基神经网络对岩土工程应力、位移非线性时间序列进行动态实时预测。网络径向基层的单元数通过均值聚类法确定后,所有其它参数：中心位置、形状参数、网络权值,均通过粒子群优化算法在全局空间优化确定。工程实例应用表明,该模型预测结果准确、精度高,有良好的应用前景。     

影响低渗透油藏低速非线性渗流的实验研究

以膨胀性的粘土为主要研究对象,通过粘土浆体的流变性实验和填砂管渗流实验,并使用一种新研制的微流量测量工具,研究了低渗岩心在较低流速下的渗流特性;利用微尺度流动理论,Einstein粘度定律等进行了理论分析,对引起低渗透油藏低速非线性渗流的影响因素进行了初步的探索.结论认为,粘土的膨胀性以及其水化后形成浆体的流变性是引起低速非线性渗流的主要因素之一.     

海洋石油管道单点提升分析

基于管道微单元体平衡建立了海管单点提升的非线性力学模型的控制微分方程组 ,使用变弧长的无量纲代换将动边界问题化为固定边界的两点边值问题 ,利用Maple环境下编制的两点边值问题的打靶法程序得到了该问题在各个提升阶段的数值解答和在单点提升过程中管道的极限弯矩约为 0 .71q1 3(EI) 2 3。     

A New Pursuit Mode of Nonlinear Wave Surface

The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bemoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. And the basis for establishing the numerical mode of nonlinear wave in time domain is to trace the position of wave free surface and to calculale the instantaneous surface height and surface potential function. This paper firstly utilizes the ‘0-1‘ combined BEM to separate the boundary by means of discretization of Green‘ s integral equation based on the Laplace equation, then separates the free surface of wave with FEM and derives the FEM equation of wave surface that satisfies the nonlinear boundary conditions. By jointly solving the above BEM and FEM equations, the wave potential and surface height could be obtained with iteration in time domain. Thus a new kind of nonlinear numerical mode is established for calculating wave transformation. The wave test in the numerical wave tank shows that the numerical simulation with this mode is of high accuracy.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

带乘性噪声的一类非线性系统的滤波算法

针对带乘性噪声的一类非线性系统,给出了1种带单重渐消因子的强跟踪状态滤波算法。该算法将非线性系统线性化后,采用了线性最小方差估计方法来进行状态估计,通过运用正交原理和引入渐消因子,使得滤波效果具有强跟踪的优良性能。该算法扩展了卡尔曼滤波在带乘性噪声非线性系统状态估计中的应用范围。仿真结果表明了该算法的有效性。