Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Solution to the rotation of the elastic earth by method of rigid dynamics     

Yoshio Kubo 《Celestial Mechanics and Dynamical Astronomy》1990,50(2):165-187

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.  相似文献   

有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献   

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland     

Jacob C. Yde  N. Tvis Knudsen 《水文研究》2004,18(1):125-139

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

渤海新近系浅水三角洲沉积体系与大型油气田勘探   总被引：32，自引：8，他引：24  

朱伟林  李建平  周心怀  郭永华 《沉积学报》2008,26(4):575-582

新近纪渤海具有构造稳定、沉降缓慢、湖泊水域范围大、内部无分割、地形平缓、坡度小等特点,广泛发育浅水三角洲,从而形成渤海海域新近系良好的储盖组合和优越的油气成藏条件。浅水三角洲砂岩储层由于受河流作用的控制,砂体具有明显的方向性,存在明显的分叉现象,与曲流河相比,砂体连通性较好,砂体等厚图呈典型的朵叶状。湖相厚层泥岩普遍发育,与浅水三角洲砂体组成多套垂向优质储盖组合。浅水三角洲发育于渤海新近纪至今处于生烃的峰期主力富烃凹陷的有利区带,新构造运动造成新近纪以来凹陷主力源岩生烃峰期与新构造运动的同步耦合,为渤海大型浅水三角洲成因的优质储层油气成藏创造了优越条件,可以形成大型油气田。因此,重视并充分挖掘渤海浅水三角洲储层类型油气藏的勘探潜力将是追求浅层优质油藏、进一步增储上产的方向所在。  相似文献   

基于Poulos弹性理论的被动桩改进算法     

陈福全  黄伟达 《岩土力学》2008,29(4):905-910

传统的Poulos弹性理论仅适合于均质土中土体侧向位移时桩的性状分析,无法考虑土的层状特性。通过引入层状地基中作用一水平集中力的广义Mindlin解和地面作用有竖向荷载时的应力和位移通解,对Poulos方法进行了改进,使之扩展到多层土中,还用于研究堆载条件下的被动桩变形和受力响应。算例分析表明,改进弹性理论要比Poulos方法更为严密、合理,提高了计算精度,应用范围也更广。  相似文献   

深水沉积层序特点及构成要素   总被引：6，自引：0，他引：6  

蒋恕  王华  Paul  Weimer 《地球科学》2008,33(6):825-833

本文在回顾当前国际上深水沉积研究热点的基础上,结合在墨西哥湾深水研究的成果系统描述了深水沉积的定义、形成机理、深水沉积层序及深水沉积构成要素的特点.深水沉积主要是在重力流作用下深水环境的沉积,主要形成于相对水平面下降和早期上升的时期,主要分布在低位体系域中.深水层序以凝缩段为边界,块状搬运沉积最早形成并直接位于层序界面上,其上被河道-天然堤沉积所覆盖.典型深水沉积的要素主要由河道、天然堤及越岸沉积、板状砂、块状搬运沉积等构成,这些沉积要素时空上有序地分布.深水河道是物源的主要通道和沉积的重要场所,从上游至下游河道弯曲度增加,能量逐渐减弱.侧向迁移明显,垂向上由富砂的顺直河道演化为相对富泥的弯曲河道.天然堤及越岸沉积以泥质为主,天然堤沿河道呈楔状分布,其近端砂岩含量高,地层厚且倾角较陡;远端砂岩含量低,地层薄且平缓,侧向连续性好但垂向连续性差.板状砂主要为深水扇前缘非限制性沉积,可分为块型和层型.块型侧向连续性好,同时垂向连通性高.层型侧向连续性好,垂向连通性差.块状搬运沉积主要是低水位期坡上沉积物失稳形成的各类滑塌体及碎屑流,其对下伏地层侵蚀明显,分布广泛,变形构造常见,可作为油气良好的封盖层.  相似文献   

Reduction of hydroelastic responses of a very-long floating structure by a floating oscillating-water-column breakwater system     

D.C. Hong  S.Y. Hong  S.W. Hong 《Ocean Engineering》2006,33(5-6):610-634

The hydroelastic responses of a very-long floating structure (VLFS) placed behind a reverse T-shape freely floating breakwater with a built-in oscillating water column (OWC) chamber are analyzed in two dimensions. The Bernoulli–Euler beam equation is coupled with the equations of rigid and elastic motions of the breakwater and the VLFS. The interaction of waves between the floating rigid breakwater and the elastic VLFS is formulated in a consistent manner. It has been shown numerically that the structural deflections of the VLFS can be reduced significantly by a suitably designed reverse T-shape floating breakwater.  相似文献   

Interference fringes on GLORIA side-scan sonar images from the Bering Sea and their implications     

Q. J. Huggett  A. K. Cooper  M. L. Somers  A. R. Stubbs 《Marine Geophysical Researches》1992,14(1):47-63

GLORIA side-scan sonographs from the Bering Sea Basin show a complex pattern of interference fringes sub-parallel to the ship's track. Surveys along the same trackline made in 1986 and 1987 show nearly identical patterns. It is concluded from this that the interference patterns are caused by features in the shallow subsurface rather than in the water column. The fringes are interpreted as a thin-layer interference effect that occurs when some of the sound reaching the seafloor passes through it and is reflected off a subsurface layer. The backscattered sound interferes (constructively or desctructively) with the reflected sound. Constructive/destructive interference occurs when the difference in the length of the two soundpaths is a whole/half multiple of GLORIA's 25 cm wavelength. Thus as range from the ship increases, sound moves in and out of phase causing bands of greater and lesser intensity on the GLORIA sonograph. Fluctuations (or wiggles) of the fringes on the GLORIA sonographs relate to changes in layer thickness. In principle, a simple three dimensional image of the subsurface layer may be obtained using GLORIA and bathymetric data from adjacent (parallel) ship's tracks. These patterns have also been identified in images from two other systems; SeaMARC II (12 kHz) long-range sonar, and TOBI (30 kHz) deep-towed sonar. In these, and other cases world-wide, the fringes do not appear with the same persistence as those seen in the Bering Sea.  相似文献   

分层流体中内孤立波在潜浮式竖直薄板上透射和反射   总被引：2，自引：0，他引：2  

魏岗  尤云祥  缪国平  苏晓冰 《海洋工程》2007,25(1):1-8

采用边缘层理论研究了两层流体系统中内孤立波在潜浮式竖直薄板上的透射和反射问题,提出了非线性演化方程的“初值”条件,分析了内孤立波与薄板非线性相互作用的效应。研究表明:流体层的密度比以及薄板伸入上下层的深度对于反射和透射波结构具有显著的影响,薄板伸入下层越深、密度差越小,则薄板阻碍孤立波透射的效率越高;透射波通常演化为单峰孤立波和迅速衰减的尾波,反射波演化为缓慢衰减的尾波列;对于具有小密度差的跃层结构,内孤立波在潜浮式竖直薄板上的透射及其演化近乎是无障碍的。  相似文献