三组结构面状态下不同产状结构面的力学效应研究

本文分析了三组结构面产状模型力学实验结果, 揭示了这种组合条件岩体中产生压剪区、张剪区及结构面、结构体的多种变形成分, 也揭示了压剪区和张剪区形成机理、及三组结构面组合下结构面产状的力学效应。     

工程岩体结构面间距空间分布的ARIMA模型及应用

结构面间距是反映结构面分布特征的重要指标之一,也是工程岩体质量级别划分的重要指标之一。与传统的结构面间距定义不同,本文认为只需按组测量结构面之间的间距同样可以分析结构面分布特征,并给出了结构面类型的4类划分法,依此建立结构面序列,利用时间序列的自回归滑动平均模型建立节理间距的AR IMA模型,提出了间距序列建模和模型检验方法,并以实例印证了AR IMA模型表示间距序列是合适的。     

Constitutive modeling of discontinuities by means of discrete and continuum approximations and damage models

In this paper the mathematical modeling of discontinuities using the discrete approximation and the continuum approximation with weak discontinuities is presented. First, the kinematics of discontinuities is discussed, then two constitutive models based on the continuum damage mechanics theory are developed. The first model is an isotropic damage model and is used in the discrete approximation. The second model is an anisotropic damage model and is used in the continuum approximation. These models are characterized for weighing the mode of failure in the failure criterion. An energy analysis is proposed to establish the equations that relate the parameters of both constitutive models; the fulfillment of the involved equations guarantee that both models are energetically equivalent. It is concluded that the proposed models are suitable to reproduce the constitutive behavior of discontinuities.     

三峡永久船闸高边坡关键块体的概率分析(英文)

The granitic rock mass existing in shiplock zone of the Three Gorges dam contains a number of major discontinuities,of them one hundred and thirty three major discontinuities have been mapped around the shiplock covering an area of 1740×600m.These major discontinuities were analyzed by using block theory.However,it is important to keep in mind that discontinuities are finite and some key blocks are actually infinite blocks.In this paper,a stochastic model is incorporated into the key block theory to study the slide probability of key blocks on the high slope of the Three-Gorges of the Yangtze River,and the size and density distribution of blocks are also analyzed.The results provide a theoretical basis for the reinforcement of the high slope on the Three-Gorges of the Yangtze River and are of great significance in engineering practice.     

基于数字钻孔影像的深部结构面类型识别及发育特征研究

深部岩体中不同类型的结构面对地下水封洞库的稳定性和水封条件产生不同的影响效应,开展深部岩体中结构面类型识别与发育特征研究具有非常重要的意义。结合黄岛地下水封洞库工程,基于数字钻孔影像技术,从曲线的形状、宽度、颜色、交错关系以及曲线周围岩石岩性等方面分析了洞库区深部发育的4类结构面的钻孔影像特征,并以ZK8中120m~-60m 段为例进行了深部结构面类型识别和发育特征研究。结果表明:(1)洞库区各类结构面具有不同的钻孔影像特征,可以借此开展深部结构面发育特征的分类研究; (2)洞库区构造节理和变质分异面发育数量沿纵深方向呈现逐渐减小的特征,且变质分异面仅发育在30m高程以上区域,岩脉侵入节理主要发育于局部较小范围,对洞室群的影响视其发育位置而异。因此应注重对深部结构面发育特征的研究,为工程建设提供翔实可靠的依据。     

岩石质量定量描述研究现状及趋向

结构面的存在破坏了岩石的完整性。定量描述结构面引起的岩石质量恶化程度是工程岩体勘测和设计的前提。本文评述了国内外岩石质量的定量描述方法，并指出今后的发展趋向。     

Modelling gas dynamics in 1D ducts with abrupt area change

Most gas dynamic computations in industrial ducts are done in one dimension with cross-section-averaged Euler equations. This poses a fundamental difficulty as soon as geometrical discontinuities are present. The momentum equation contains a non-conservative term involving a surface pressure integral, responsible for momentum loss. Definition of this integral is very difficult from a mathematical standpoint as the flow may contain other discontinuities (shocks, contact discontinuities). From a physical standpoint, geometrical discontinuities induce multidimensional vortices that modify the surface pressure integral. In the present paper, an improved 1D flow model is proposed. An extra energy (or entropy) equation is added to the Euler equations expressing the energy and turbulent pressure stored in the vortices generated by the abrupt area variation. The turbulent energy created by the flow–area change interaction is determined by a specific estimate of the surface pressure integral. Model’s predictions are compared with 2D-averaged results from numerical solution of the Euler equations. Comparison with shock tube experiments is also presented. The new 1D-averaged model improves the conventional cross-section-averaged Euler equations and is able to reproduce the main flow features.     

Modeling of hard disk drives for shock and vibration analysis – consideration of nonlinearities and discontinuities

This paper reviews recent advances (mostly after year 2000) in shock and vibration analysis of hard disk drives (HDD) considering the presence of nonlinearities and discontinuities. Components and dynamic phenomena in HDD where effects of mechanical nonlinearity and discontinuities are significant are discussed, e.g., head actuator suspension, dimple and slider, head–disk interface, fluid dynamic bearing, spinning disks, and load/unload dynamics. Ways to model these nonlinearities and discontinuities are reviewed in detail. Our research on modeling an entire HDD in operating mode subject to shock and vibration using a flexible multibody dynamics formulation is also presented. The numerical simulation of the shock response of a 1-in. form factor HDD is presented. A half-sinusoidal acceleration shock is applied at the base of the HDD. Response of the flying height for different shock amplitudes and duration times is simulated.     

Kinked crack analysis by a hybridized boundary element/boundary collocation method

In this research a two dimensional displacement discontinuity method (which is a kind of indirect boundary element method) using higher order elements (i.e. a source element with a cubic variation of displacement discontinuities having four sub-elements) is used to obtain the displacement discontinuities along each boundary element. In this paper, three kinds of the higher order boundary elements are used: the ordinary elements, the kink elements and the special crack tip elements.The boundary collocation technique is used for the calculation of the displacement discontinuities at the center of each sub-elements. Again a special boundary collocation technique is used to treat the kinked source elements occur in the crack analysis. Considering the two source elements (each having four sub-elements) joined at a corner (kink point). The collocation points in the cubic element model which are outside of the kink point are moved to the crack kink then the displacement discontinuities on the left and right sides of the kink are calculated. The displacement discontinuities of the kink point are obtained by averaging the corresponding values of its left and right sides. The special crack tip elements are also treated by the boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. Some simple example problems are solved numerically by the proposed method. The numerical results are compared with the corresponding results obtained by the previous methods cited in the literature. This comparison shows a very good agreement between the results and verify the accuracy and validity of the proposed method.     

Investigation of surfaces of discontinuity in perfectly conducting magnetizing media

Relationships on discontinuities in magnetizing perfectly conducting media in a magnetic field are investigated. The magnetic permeabilities before and after the discontinuity are assumed to be constant, but unequal, quantities. It is shown that shocks of two kinds, fast and slow, are possible in the formulation under consideration in the hydrodynamics of magnetizing media, as in magnetic hydrodynamics: It is shown that the entropy decreases on the rarefaction shocks diminishing the magnetic permeability, but can grow on the rarefaction shocks increasing the magnetic permeability, but such waves are not evolutionary. The relationships on discontinuities in the mechanics of a continuous medium are written down in general form in [1] with the electromagnetic field, polarization, and magnetization effects taken into account. Relationships on discontinuities in the ferrohydrodynamic and elec trohydrodynamic approximations were written down in [2] and [3–5], respectively, for the cases when the magnetic permeability and dielectric permittivity of the medium ahead of and behind the discontinuity are arbitrary functions of their arguments and are identical. A system of relationships on discontinuities propagated into a magnetizing perfectly conducting medium is investigated in this paper. The method proposed in [6] is used in the investigation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 104–110, January–February, 1976.We are grateful to A. A. Barmin for discussing the paper and for valuable remarks.     

Stability of surfaces of strong discontinuity in gases

The existence of solutions with surfaces of strong discontinuity is one of the principal features of the continua whose motions are described by systems of differential equations of hyperbolic type. Shock waves in gas dynamics, magnetohydrodynamics and in solids, detonation waves and combustion fronts, contact discontinuities, etc. are well-known examples of these surfaces. The discontinuities are usually investigated in accordance with the following scheme: 1) derivation of the boundary conditions on the discontinuity from the input system of differential equations in integral form; 2) verification of the fulfilment of the evolution conditions; 3) solution of the problem of the discontinuity structure and, when the occasion requires, obtaining supplementary boundary conditions; 4) investigation of the stability of the discontinuity. Only after obtaining positive results in all fours stages can we assert that the existence of the discontinuity is theoretically justified and that it can be used for constructing the solutions of particular boundary value problems. In the present paper attention will be concentrated on the problem of the stability of discontinuities, all the material, with the exception of the general results of Sec.1, being concerned with gas media and relating to discontinuities on whose surface the normal mass flow is nonzero. Having no way of exploring all the aspects of the problem of the stability of discontinuities in the same detail within the limited context of this paper, the authors hope to demonstrate the most general ideas and approaches which could subsequently be used to investigate the stability of discontinuities in various particular models of continua.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–22, March–April, 1996.     

结构面抗剪强度参数对岩质边坡稳定的影响

边坡工程中,结构面抗剪强度是非常重要的力学参数,合理地选择确定抗剪强度的方法,常常要分析c、φ值在边坡稳定中的作用。根据极限平衡原理,对10~30m高的平面滑动型边坡不同抗剪强度的控稳结构面c、φ值在边坡稳定中的作用进行了定量分析。研究表明c值对稳定系数的影响随坡高增大逐渐减弱,φ值对稳定系数的影响随坡高增大逐渐增强。得到了不同坡高下c、φ值边坡稳定作用等值曲线,该曲线很容易确定抗剪强度参数c、φ值在边坡稳定中的作用。     

Closed form solutions of Euler–Bernoulli beams with singularities

The problem of the integration of the static governing equations of the uniform Euler–Bernoulli beam with discontinuities is studied. In particular, two types of discontinuities have been considered: flexural stiffness and slope discontinuities. Both the above mentioned discontinuities have been modeled as singularities of the flexural stiffness by means of superimposition of suitable distributions (generalized functions) to a uniform one dimensional field. Closed form solutions of governing differential equation, requiring the knowledge of the boundary conditions only, are proposed, and no continuity conditions are enforced at intermediate cross-sections where discontinuities are shown. The continuity conditions are in fact embedded in the flexural stiffness model and are automatically accounted for by the proposed integration procedure. Finally, the proposed closed form solution for the cases of slope discontinuity is compared with the solution of a beam having an internal hinge with rotational spring reproducing the slope discontinuity.     

群速度控制格式及二维Riemann解

强激波和强接触间断的数值模拟一直是计算流体力学里一个富有挑战性的课题,它们是很多实际流动的基础。三阶迎风紧致格式是一种具有较高分辨率的高精度方法,但是在计算激波时仍有数值振荡产生。本文根据数值解的群速度特性,在三阶迎风紧致格式的基础上提出了一种群速度控制格式,使得能够正确模拟含有强激波和强接触间断的复杂流动。在此基础上构造了求解包含大压力比和密度比的二维界面问题的数值方法。计算结果表明,方法对激波和接触间断的分辨效果是令人满意的。     

Recent Development on the Conservation Property of Chimera

An estimate on the conservation error due to the non-conservative data interpolation scheme for overset grids is given in this paper. It is shown that the conservation error is a first-order term if second-order conservative schemes are employed for the Chimera grids and if discontinuities are located away from overlapped grid interfaces. Therefore in the limit of global grid refinement, valid numerical solutions should be obtained with a data interpolation scheme. In one demonstration case the conservation error in the original Chimera scheme was shown to affect flow even without discontinuities on coarse to medium grids. The conservative Chimera scheme was shown to give significantly better solutions than the original Chimera scheme on these grids with other factors being the same.     

Numerical investigation of the slope discontinuities in large deformation finite element formulations

In many multibody system applications, the system components are made of structural elements that can have different orientations, leading to slope discontinuities. In this paper, a numerical investigation of a new procedure that can be used to model structures with slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is presented. This procedure can be applied to model slope discontinuities in the case of commutative rotations of gradient deficient elements that are used for modeling thin beam and plate structures. An important special case to which the proposed procedure can be applied is the case of all planar gradient deficient ANCF finite elements. The use of the proposed method leads to a constant orthogonal element transformation that describes an arbitrary initial configuration. As a consequence, one obtains, in the case of large commutative rotations and large deformations, a constant mass matrix for structures which have complex geometry. The procedure used in this investigation to model slope discontinuities requires the use of the concept of the intermediate finite element coordinate system. For each finite element, a new set of gradient coordinates that define, at the discontinuity node, the element deformation with respect to the intermediate element coordinate system is introduced. These new gradient coordinates are assumed to be equal for the two finite elements at the point of intersection. That is, the change of the gradients of two elements at the intersection point from their respective intermediate initial reference configuration is assumed to be the same. This procedure leads to a set of linear algebraic equations that define the orthogonal transformation matrix for the finite element. Numerical examples are presented in order to demonstrate the use of the proposed procedure for modeling slope discontinuities.     

Transformation of hanging discontinuities into vortex systems in a stratified flow behind a cylinder

The mechanisms of formation of hanging discontinuities, vortex dipoles, and vortex arrays in the wave wake behind a cylinder moving at a constant velocity in a stratified fluid are investigated using various schlieren methods. The existence of discontinuities is attributable to the distortion of the internal-wave phase pattern in the shear flow and to the varying stratification and subsequent interaction of the waves with the appearing nonuniformities. Hanging discontinuities and vortex systems are low-velocity analogs of shock waves. An analysis of the internal-wave pattern indicates that the values of the normal velocity component differ on the upper and lower edges of the discontinuities. A regime diagram for flows of this kind is given.     

High-resolution algorithms of sharp interface treatment for compressible two-phase flows

In this paper, a kind of arbitrary high order derivatives (ADER) scheme based on the generalised Riemann problem is proposed to simulate multi-material flows by a coupling ghost fluid method. The states at cell interfaces are reconstructed by interpolating polynomials which are piece-wise smooth functions. The states are treated as the equivalent of the left and right states of the Riemann problem. The contact solvers are extrapolated in the vicinity of contact points to facilitate ghost fluids. The numerical method is applied to compressible flows with sharp discontinuities, such as the collision of two fluids of different physical states and gas–liquid two-phase flows. The numerical results demonstrate that unexpected physical oscillations through the contact discontinuities can be prevented effectively and the sharp interface can be captured efficiently.     

Non-homogeneous displacement jumps in strong embedded discontinuities

In this paper, strong discontinuities are embedded in finite elements to describe fracture in quasi-brittle materials. A new numerical formulation is introduced in which the displacement jumps do not need to be homogeneous within each finite element. Both the crack path and the displacement jumps are continuous across element boundaries. This formulation is compared with the discrete approach, in which interface elements are inserted to model the discontinuities, as well as with other embedded discontinuity approaches and with the partition of unity method. Numerical results have been obtained with relatively coarse meshes, which compare well with experimental results and with the results obtained from analyzes with interface elements.     

Three-dimensional interfacial fracture analysis of a one-dimensional hexagonal quasicrystal coating

In this paper, the three-dimensional (3D) interfacial fracture is analyzed in a one-dimensional (1D) hexagonal quasicrystal (QC) coating structure under mechanical loading. A planar interface crack with arbitrary shape is studied by a displacement discontinuity method. Fundamental solutions of interfacial concentrated displacement discontinuities are obtained by the Hankel transform technique, and the corresponding boundary integral-differential equations are constructed with the superposition principle. Green’s functions of constant interfacial displacement discontinuities within a rectangular element are derived, and a boundary element method is proposed for numerical simulation. The singularity of stresses near the crack front is investigated, and the stress intensity factors (SIFs) as well as energy release rates (ERRs) are determined. Finally, relevant influencing factors on the fracture behavior are discussed.