Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling discontinuous rock mass based on smoothed finite element method

This paper aims at developing a method for modeling rock mass with preexisting multiple discontinuities within the framework of the smoothed finite element method (SFEM). The discontinuity is simulated by an interface element with zero thickness, the stiffness matrix of which are derived explicitly based on the SFEM. An elastic damage constitutive relation with residual strength is introduced in order to describe the nonlinear mechanical behavior of the discontinuities. The computation codes of the present method were developed. The present method has been verified to be a sound approach for modeling discontinuous rock mass, inheriting the advantages of the SFEM.     

Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling

The initiation and propagation of directional hydraulic fracturing (DHF) was investigated based on true tri-axial experiment and finite element modeling. The influences of notch angle, notch length and injection rate on the DHF were investigated. The initiation and propagation of DHF was modeled by a 3D nonlinear finite element method. A comparison between experimental investigation and numerical modeling results indicates that there is a good correlation between unbalanced force (UF) and fracturing. UF can be used to predict the hydraulic fracture initiation and propagation.     

Coupled Model of Two-phase Debris Flow, Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper,we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth,the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress,and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition(Iverson and Ouyang,2015) is presented. In this formula,the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe's Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.     

构造物理模拟和PIV有限应变分析对构造裂缝预测的启示

油气储层中构造裂缝发育与有限应变状态关系密切,为了探索有限应变分析与构造裂缝预测的新技术方法,此次研究设计完成了一组单侧挤压收敛模型的物理模拟实验,并引入粒子图像测速(PIV,Particle Image Velocimetry)技术对实验过程进行了定量化分析。实验模型在垂向上为含粘性层的多层结构,实验结果形成了一个肉眼可见的箱状褶皱。通过PIV技术可以获取实验模型变形演化过程中各阶段的位移场数据,计算出各阶段的增量应变,实现从初始状态到褶皱形成之后整个变形过程的有限应变分析,探讨构造裂缝成因机制和分布规律,进行定量化裂缝预测。挤压变形过程初期,应变分布范围很广,有限应变较弱(约4%~8%),在挤压方向上的线应变表现为弱压应变,在垂向上的线应变表现为弱张应变,这种现象是褶皱和断层产生前平行层缩短和层增厚的纯剪变形结果,也是区域型张裂缝和剪裂缝形成的主要机制。褶皱和断层即将发育之时至发育之后,应变局限在断层发育的剪切带及附近区域,有限应变表现为较强(达20%)的剪切应变和剪切张应变,是断层面附近简单剪切变形作用的结果,也是局部型剪裂缝和张剪裂缝形成的主要机制。     

脐带缆在潜器下放过程中的运动建模与仿真

潜器在下放布置过程中需要考虑变长度的脐带缆对其运动的影响,在考虑流的情况下,对脐带缆的有限差分模型进行了改进,并采用最小二乘方法求解由时间和空间上的中心差分格式离散后的非线性方程组。为了验证模型改进的有效性,将潜器在均匀定常流中水平匀速直航时的数值解与该特殊情况下的解析解进行比较,两组解的吻合不仅证实了模型改进的有效性,而且表明上述数值计算方法是可靠且有效的。在潜器下放布置运动给定的情况下,脐带缆的运动仿真结果表明,潜器姿态角的调整会对拖曳点处的缆绳张力大小和变化趋势产生显著的影响,流在改变缆绳空间形状的同时引起了缆绳对潜器作用力和力矩的非线性和时变性,脐带缆总是阻止潜器自身动力对运动速度和姿态角的改变,且对纵荡运动的影响最大。     

The Effect of Concrete Deformation on Displacement of an Axially Loaded Drilled Shaft

This article presents the settlement of drilled shafts resulting from their structural deformations. Although drilled shafts are widely used as foundations for settlement-sensitive structures such as bridges and high-rise buildings, the structural deformations of drilled shafts are not typically taken into account in the design process. However, if unexpected structural deformations of drilled shafts cause additional settlement to the foundation, the serviceability of the superstructure can be jeopardized. Unfortunately, very few research efforts have been made to quantify the structural deformation of drilled shafts; this needs to be addressed to accurately predict the settlement of drilled shafts. In this study, we investigate the effect of structural deformation on displacement of axially loaded drilled shafts. Finite element analyses were performed to quantify the structural deformation of drilled shafts. The analysis results indicated that the structural deformation of drilled shafts could be quite significant for long drilled shafts. The main factors that affected the structural deformation of drilled shafts were found to be pile length, the material properties of drilled shafts, and the relative humidity of surrounding soil. An approximate equation is proposed to estimate the long-term deformation of drilled shafts.     

Multi-factor sensitivity analysis on the stability of submarine hydrate-bearing slope

There are many factors affecting the instability of the submarine hydrate-bearing slope (SHBS), and the interaction with hydrate is very complicated. In this paper, the mechanical mechanism of the static liquefaction and instability of submarine slope caused by the dissociation of natural gas hydrate (NGH) resulting in the rapid increase of pore pressure of gas hydrate-bearing sediments (GHBS) and the decrease of effective stress are analyzed based on the time series and type of SHBS. Then, taking the typical submarine slope in the northern South China Sea as an example, four important factors affecting the stability of SHBS are selected, such as the degree of hydrate dissociation, the depth of hydrate burial, the thickness of hydrate, and the depth of seawater. According to the principle of orthogonal method, 25 orthogonal test schemes with 4 factors and 5 levels are designed and the safety factors of submarine slope stability of each scheme are calculated by using the strength reduction finite element method. By means of the orthogonal design range analysis and the variance analysis, sensitivity of influential factors on stability of SHBS are obtained. The results show that the degree of hydrate dissociation is the most sensitive, followed by hydrate burial depth, the thickness of hydrate and the depth of seawater. Finally, the concept of gas hydrate critical burial depth is put forward according to the influence law of gas hydrate burial depth, and the numerical simulation for specific submarine slope is carried out, which indicates the existence of critical burial depth.     

页岩压裂影响因素三维模型数值分析

不同影响因素对页岩水力压裂效果有不同的影响。基于三维数值计算模型介绍了水力压裂的典型过程,针对水平最小主应力、页岩的弹性模量和渗透系数对压裂的影响进行了分析。结果表明:随着水平最小主应力的增加,裂缝高度也随之增加;随着页岩的弹性模量增加,裂缝的高度随之降低;随着渗透系数的增加,裂缝高度也随之增加;随着压裂液的黏度系数提高,裂缝的高度降低。