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

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

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.     

有限元法与伪谱法混合求解弹性波动方程

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

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

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

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

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.     

Analytic Solution of Interaction Between Surface Waves and Floating Elastic Plates with Elastic End-Restraints

Based on the dynamic theories of water waves and Mindlin plates,the analytic solution of interaction between surface waves and two-dimensional floating elastic plates with edge-restraint is constructed by use of the Wiener-Hopf technique.Firstly,without regard for elastic edge restraint,the wave-induced responses of elastic floating plate analyzed by the present method are in good agreement with the results from literature and experimental results.Therefore,it can be shown that the present method is valid.Secondly,three end-restraint cases(i.e.,the left-end elastic restraints,the both-end elastic restraints,and the right-end elastic restraints)are proposed to reduce the vibration of floating plates,in which the spring is used to connect the sea bottom and the floating plate's left(or right)edge.The relations between the spring stiffness and the parameters of wave-induced responses of floating plates are discussed.Moreover,the effective method to reduce the vibration of floating elastic plates can be obtained through comparison.     

对我国西南地区河谷深厚覆盖层成因机理的新认识

近年来,在我国水能资源开发过程中,发现各河流现代河床以下普遍堆积厚达数十米甚至上百米的松散堆积物。河谷深厚覆盖层的存在,不仅严重制约了工程坝址的选择,影响相关流域水电资源的开发利用,也给坝工设计带来巨大的困难。由于深厚覆盖层埋藏于现代河床之下,其形成年代一般先于一二级阶地,有悖于河流发育演化的常理,其成因一直令人费解。首次将河谷深切和深厚堆积事件与全球气候变化、海平面升降运动、地壳运动等有机地联系起来,并提出冰期、间冰期全球海平面大幅度升降,是导致河流深切成谷并形成深厚堆积的主要原因的新观点。在此基础上,引入层序地层学原理,从理论上较好地解释了全球气候变化导致海平面和河流侵蚀基准面大幅变化,并产生河谷深切和深厚堆积的原因和过程。最后,进一步将沿河大型古滑坡的孕育和发生与河谷深切事件相联系,提出沿河大型古滑坡是在河谷深切期因前缘临空较好而形成的新观点,从而对沿河古滑坡前缘剪出口高程往往低于现代河床数十米的原因给出了较合理的解释。     

基于弹性薄板的地表沉陷预计模型

在分析地表沉陷基本规律的基础上,依据弹性薄板理论,建立非充分采动条件下岩层和地表沉陷预计的一类新模型,并推导出地表任意点倾斜、曲率、水平移动以及水平变形的计算公式。该模型充分考虑到地质采矿因素（煤层倾角）及煤层上方各岩层的影响,克服传统预测方法的缺陷,特别是概率积分法关于拐点反对称要求。最后,应用实例证明该方法的应用效果。     

肋柱式锚杆支挡结构的徐变次内力分析

采用老化徐变理论和等效弹性模量法,对肋柱式锚杆支挡结构的徐变影响提出了具体计算方法,并通过实例计算了支挡结构的徐变次内力。计算结果表明,在逆作法条件下,随着开挖深度的加大,岩土侧向压力也增大,使后期施工的混凝土构件产生的徐变次内力也随之增大;考虑徐变后使结构设计更加安全可靠,尤其坡脚处后期浇筑的混凝土构件,更要关注徐变的影响。     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.