Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model          下载免费PDF全文

T. Saksala  D. Brancherie  A. Ibrahimbegovic 《国际地质力学数值与分析法杂志》2016,40(9):1339-1357

This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

On the formulation of a finite element method for the stiffened multi-layered airfoil/hydrofoil structure: Post buckling analysis for the wings of underwater gliders     

《Applied Ocean Research》2017

A finite element method is developed for the stiffened multi-layered airfoil/hydrofoil structure for the large deformation and finite strain problem. The kinematics of the airfoil/hydrofoil is set up. The Consistent Orthogonal Basis Function Space is applied for the airfoil/hydrofoil structure. Given the airfoil/hydrofoil configuration and boundary conditions, the basis function space can be uniquely determined, such that the diagonal mass matrix is obtained accurately and the basis functions are very identical with the mode shape functions of the structure. In order to satisfy the displacement compatibility condition between adjacent layers of the airfoil/hydrofoil, the traction degree-of-freedom is also induced.The post buckling analysis is presented for the wing (hydrofoil) structure of the underwater glider. The water pressure is applied on the outer surface and the critical buckling pressure is calculated. The post buckling equilibrium path is also given. The results are verified with ANSYS. The present study of the buckling analysis of the airfoil/hydrofoil under water pressure is helpful to the design of underwater glider.  相似文献   

Finite volume approach for finite strain consolidation          下载免费PDF全文

Rakesh Pratap Singh  Chandra Shekhar Prasad Ojha  Mahendra Singh 《国际地质力学数值与分析法杂志》2016,40(1):117-140

The paper presents the finite volume formulation and numerical solution of finite strain one‐dimensional consolidation equation. The equation used in this study utilises a nonlinear continuum representation of consolidation with varying compressibility and hydraulic conductivity and thus inherits the material and geometric nonlinearity. Time‐marching explicit scheme has been used to achieve transient solutions. The nonlinear terms have been evaluated with the known previous time step value of the independent variable, that is, void ratio. Three‐point quadratic interpolation function of Lagrangian family has been used to evaluate the face values at discrete control volumes. It has been shown that the numerical solution is stable and convergent for the general practical cases of consolidation. Performance of the numerical scheme has been evaluated by comparing the results with an analytical solution and with the piecewise piecewise‐linear finite difference numerical model. The approach seems to work well and offers excellent potential for simulating finite strain consolidation. Further, the parametric study has been performed on soft organic clays, and the influence of various parameters on the time ate consolidation characteristics of the soil is shown. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

西南极洲纳尔逊冰帽的表面运动          下载免费PDF全文

刘琛  秦大河  王晓军  钱嵩林 《极地研究》1991,3(4):11-17

本文通过对纳尔逊冰帽冰面速度和应变率的分析,得出冰帽边缘的最大夏季表面速度为2058m/a,最大年平均表面速度为15.5m/a,探求了冰面运动速度的空间分布和季节性变化,冰帽顶部的表面应变率为+0.0079/a和+0.0034/a。还讨论了冰面应变率从冰帽顶点到边缘的变化特点及其与冰面裂隙的关系。并将这些特征同我国大陆型冰帽进行了对比,从冰川动力学的角度阐明了纳尔逊冰帽的某些特征。  相似文献   

Interaction between different internal length scales in strain localization analysis of fully and partially saturated porous media—the 1‐D case     

B.A. Schrefler  H.W. Zhang  L. Sanavia 《国际地质力学数值与分析法杂志》2006,30(1):45-70

The paper analyses the interaction between two internal length scales during dynamic strain localization in multiphase porous materials. The first internal length is introduced in the mathematical model by the gradient‐dependent plasticity for the solid skeleton, while the second one is naturally contained in the multiphase model and is due to the seepage process of the water via Darcy's law, which induces a rate‐dependent behaviour of the solid skeleton. Numerical results of a one‐dimensional example of water saturated porous medium demonstrate the competing effect between these two length scales. The porous medium is here treated as a multiphase continuum, with the pores filled by water and air, the last one at constant atmospheric pressure. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

Inelastic constitutive properties and shear localization in Tennessee marble     

D. J. Holcomb  J. W. Rudnicki 《国际地质力学数值与分析法杂志》2001,25(2):109-129

The inelastic response of Tennessee marble is modelled by an elastic plastic constitutive relation that includes pressure dependence of yield, strain‐softening and inelastic volume strain (dilatancy). Data from 12 axisymmetric compression tests at confining pressures from 0 to 100 MPa are used to determine the dependence of the yield function and plastic potential, which are different, on the first and second stress invariants and the accumulated inelastic shear strain. Because the data requires that the strain at peak stress depends on the mean stress, the locus of peak stresses is neither a yield surface nor a failure envelope, as is often assumed. Based on the constitutive model and Rudnicki and Rice criterion, localization is not predicted to occur in axisymmetric compression although faulting is observed in the tests. The discrepancy is likely due to the overly stiff response of a smooth yield surface model to abrupt changes in the pattern of straining. The constitutive model determined from the axisymmetric compression data describes well the variation of the in‐plane stress observed in a plane strain experiment. The out‐of‐plane stress is not modelled well, apparently because the inelastic normal strain in this direction is overpredicted. In plane strain, localization is predicted to occur close to peak stress, in good agreement with the experiment. Observation of localization on the rising portion of the stress–strain curve in plane strain does not, however, indicate prepeak localization. Because of the rapid increase of mean stress in plane strain, the stress–strain curve can be rising while the shear stress versus shear strain curve at constant mean stress is falling (negative hardening modulus). Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

用变形分析方法研究鲜水河断裂带的地震形变     

刘本培  李建中  陈永奇  吴子安  孙桂芳  杜自钦 《大地测量与地球动力学》1995,(3)

提出了一种变形分析和用Ｒｏｂｕｓｔ原理进行位移场变换的方法，并用以计算了鲜水河断裂带的水平形变，给出了比较客观的位移场。按均匀应变模型计算了点和图形的均匀应变和分区应变的加权平均值，结果说明主应变轴方位与破裂线方位相关；在北西走向的主干断裂上，无论是破裂区或闭锁区，主压应变轴方位都是近东西向的。采用样条函数和截断多项式函数拟合跨断裂的点，用非线性动态平差模型求得了地震区的垂直形变模型，结果显示了道孚地震震前的加速形变、震时以震中为中心的下沉和震后的反向回升运动。  相似文献   

异向性岩芯冰的反复压缩-退火实验研究          下载免费PDF全文

黄茂桓  王文悌  李军  李纲 《极地研究》1990,2(4):21-27

在蠕变试验机上完成了三组(6块样品)单轴反复压缩-退火实验,试验机载荷精度为1％,实验温度为-2±0.2℃。实验样品取自南极洛多姆冰帽BHQ钻孔的冰岩芯。三组试样的受压方向分别与冰岩芯的轴向平行、成45°角和垂直。试验前试样的C轴组构为单极大型或大致为单极大型,晶粒大小不一。用0.8MPa的初始轴向应力对样品进行压缩,达到10％应变后卸载退火72小时。每组实验如此反复进行6次。结果表明,在较高的温度和大载荷下,样品经反复压缩-退火,其初始组构特征消失,代之以小环状组构和细的、等粒的晶体,某些情况下发展为多极大型组构。结构和组构分析表明,新组构的形成机制主要是再结晶作用。随着压缩-退火反复进行,六块样品在组构上的差异减少,它们的流变特性趋于一致,晶粒尺寸逐渐减少而趋于一个稳定值。  相似文献   

论基岩滑坡的临床聚能与启程弹冲     

胡广韬  毛彦龙  杨建宏 《地球科学与环境学报》1992,(4)

剧动式高速滑坡是一种大量存在于自然界中的地质灾害。本文以广家崖斜坡危岩体为例,对基岩滑坡的“临床聚能”和“启程弹冲加速效应”进行了理论分析和探讨。文中详细论证了斜坡岩体滑落破坏时所具有的“启程弹冲加速效应”,并提出了二维应力状态下的启程冲速度计算公式。 启程弹冲机理,主要是锁固段岩体内聚力作用和岩体聚集赋存弹性应变能释放的反应,弹冲速度与内聚力的大小有直接关系。岩体一旦剪断破坏,弹性应变能骤然释放,岩体表现出剧烈的弹性回冲,弹性应变能便瞬间转化为岩体滑动的弹冲动能。这便从理论上对高速滑坡提出了具有实际意义的新机理。  相似文献   

广西下、中三叠统砾质岩类特征及成因探讨     

张继淹 《南方国土资源》1997,(3)

文章分析论述了广西下、中三叠统砾质岩类所在层位及其特征，并将其划分为水道砾岩、崩塌—滑塌砾岩、沉积混杂岩３种成因类型，探讨了形成环境及其意义，指出它们均是三叠纪时期构造作用的沉积表现。  相似文献