U-transformation-finite element method in 3-dimensional analysis of orthotropic laminates

For an orthotropic laminate, an equivalent system with doubly cyclic periodicity is introduced. Then a 3-dimensional finite element model for the equivalent system is transformed into the unitary space, where the large finite element matrix equation is decoupled into some small matrix equations. Such a decoupling very efficiently reduces the computational effort. For an orthotropic laminate with four clamped edges, no exact elasticity solution is available, and the deflection values predicted by different methods have a considerable difference each other for a small length-to-thickness ratio. The present predictions are the largest because the present method is a full 3-dimensional finite element analysis without superfluous constraints. Illustrative numerical examples are presented to observe the distributions of stresses through the thickness of the laminates. (C) 2010 Elsevier Ltd. All rights reserved.

Fracture analysis for 3-dimensional bodies with a surface crack

This paper deals with fracture analyses in 3-dimensional bodies containing a surface crack. A general solution of stress-strain fields at crack tip is proposed. Based on the stress-strain fields obtained, a high-order 3-dimensional special element is established to calculate the stress intensity factors in a plate with a surface crack. The variation of stress intensity factors with geometric parameters is investigated.

三维离散元改进算法及其在三峡永久船闸高边坡工程中的应用

本文根据面－面接触模型，写出了三维离散元的物理方程及运动方程，对于接触刚度、阻尼系数、时步、计算控制及块体的滑动失稳判据等作了讨论。在NURBM－3D 长方体单元基础上改进了其划分单元的算法，依据岩体中三组真实的节理面划分平行六面体单元，实现了网格的自动剖分，从而可以将单元视为真实的岩块，更客观地体现了节理面对岩体稳定性的影响。运用改进后的三维离散元方法，编制了计算程序，对三峡工程永久船闸高连坡开挖进行了模拟。计算过程包括初始地应力场模拟、岩体开挖过程两个阶段。模拟计算首先给出了为水平、垂直节理面所切割的长方体单元情况下的应力场分析，考证了程序的可行性。将已有的三峡永久船闸实测节理简化为三组节理，完成了初始地应力场的模拟，给出了分四步开挖后的岩体位移场和应力场，计算出岩体开挖过程中产生滑动块体的位置，其结果可以定性地表明：沿节理面的块体滑动是边坡失稳的主要模式。编制的后处理程序提供了计算数据库与绘图软件Origin的接口，可以画出任意剖面或节理面上的应力等值线、位移矢量图及滑动块体分布图，由此可以直观判断岩体不稳定区域和高应力区范围。编制的计算程序在PC机上计算了10万块体单元，用于模拟大型工程，其计算结果与有限元及类似的连续介质模型相比更为合理，为用离散元法进一步研究可变形块体及岩体的蠕变提供了一种切实可行的工具。将离散元法应用于大型工程三维计算是可行的，但在本构模型建立、岩体破坏参数选择等方面均有等于更深入的研究，其结果也有待于更多工程的验证。

三峡永久船闸高边坡开挖三维离散元数值模拟:3-D Simulation of the Excavation of High Steep Slope of Three-Gorges Permanent Lock by Distinct Element Method

采用面一面接触的三维离散元刚性块体模型,从实测节理面中取出其中的三组,按照其倾向、倾角和节理间距将三峡永久船闸未开挖的区域划分为10~5个离散单元,通过施加力边界条件,给出了与实测初始地应力场接近的数值模拟结果；然后,分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势,个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Distinct Element Analysis on Propagation Characteristics of P-Wave in Rock Pillar with Finite length

以节理岩体等效刚度的概念为基础,讨论了离散元刚性块体模型中节理刚度的选取问题。采用面-面接触模型模拟了纵波在一维岩体中的传播,给出了纵波波形；研究了阻尼比、软弱夹层以及节理间是否可拉对波传播规律的影响。

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.

Study on cracked plates, shells and 3 dimensional bodies

This paper presents a summary of the authors' recent work in following areas: (1) The stress-strain fields at crack tip in Reissner's plate. (2) The calculations of the stress intensity factors in finite size plates. (3) The stress-strain fields at crack tip in Reissner's shell. (4) The calculations of the stress intensity factors and bulging coefficients in finite size spherical shells. (5) The stress-strain fields along crack tip in three dimensional body with surface crack. (6) The calculation of stress intensity factors in a plate with surface crack.

Short-surface waves radiated by a partially-immersed 3-dimensional oscillating body

The short-surface waves generated by a 3-D arbitrarily oscillating body floating onwater are discussed. In the far-field off the body, the phase and the amplitude functions ofthe radiated waves are determined by the ray method. An undetermined constant is includ-ed in the amplitude function. From the result of Ref. [1], the near-field boundary layersolution near the body waterline is obtained. The amplitude of this solution depends on thewhole wall shape of the body and the slope at the body waterline on the cross-sections per-pendicular to the waterline. By matching the far-field solution with the near-field bound-ary layer solution, the undetermined constant in the amplitude function of the far-fieldradiated waves is determined. For the special case of a half-submerged sphere which per-forms vertical oscillating motion, the result obtained in this paper is in agreement withthat of Ref. [ 2 ].

One-dimensional improvement and two-dimensional and 3-dimensional treatments for stable oscillation condition, mode structure and power output in high-speed-flow lasers

For high-speed-flow lasers, the one-dimensional and first-order approximate treatment in[1] under approximation of geometrical optics is improved still within the scope of approx-imation of geometrical optics. The strict accurate results are obtained, and what is more,two- and three-dimensional treatments are done. Thus for two- and three-dimensional cases, thestable oscillation condition, the formulae of power output and analytical expression of modesunder approximation of geometrical optics (in terms of gain function) are derived. Accord-ing to the present theory, one-and two-dimensional calculations for the typical case of Gerry'sexperiment are presented. All the results coincide well with the experiment and are better thanthe results obtained in [1].In addition, the applicable scope of Lee's stable oscillation condition given by [1] is ex-panded; the condition for the approximation of gcometrical optics to be applied to mode con-structure in optical cavity is obtained for the first time and the difference between thiscondition and that for free space is also pointed out in the present work.

SELF-CONSISTENT TREATMENT OF 3-DIMENSIONAL - 2-DIMENSIONAL AND 2-DIMENSIONAL - 2-DIMENSIONAL RESONANT-TUNNELING IN DOUBLE-BARRIER STRUCTURES

By using a transfer-matrix method on the basis of two-dimensional (2D) Bloch sums in accordance with a tight-binding scheme, a self-consistent calculation on the resonant tunneling in asymmetric double-barrier structures is presented, in which contributions to resonant tunneling from both three-dimensional (3D) electrons in the contacts and 2D electrons in the spacer or accumulation layers are considered simultaneously. The charge buildup effect on the current versus voltage (I-V) curves is evaluated systematically, showing quantitatively how it results in the I-V bistability and enhanced differences between I-V curves for positive and negative bias in an asymmetric double-barrier structure. Special attention is focused on the interaction between 3D-2D and 2D-2D resonant-tunneling processes, including the suppression of 2D-2D resonant tunneling by the charge buildup in the well accompanying the 3D-2D resonant tunneling. The effects of the emitter doping condition (doping concentration, spacer thickness) on the presence of two types of quasi-2D levels in the emitter accumulation layers, and on the formation of a potential bulge in the emitter region, are discussed in detail in relation to the tunneling process.

层状反倾向结构岩质边坡变形机理及稳定性研究

Toppling is a major failure model in anti-dip layered rock slopes. Because of the limited by testing means and experimental apparatus, present research on the deformation mechanism and stability analysis are mainly focus on the 2-Dimensional deformation, and the research really based on 3-Dimension is still limited. Therefore, based on the present research station, the article rely on the important hydroelectric project of typical anti-dip layered rock slopes -- The left bank slope of Long-tan hydropower-station in Guang Xi, China, and focused on the influencing factors, deformation mechanism and stability analysis of anti-dip layered rock slopes, three problems as follows are researched in this paper. (1) Deformation influencing factor analysis on ant-dip layered rock slopes Three influencing factors are included: geological factor, engineering factor and environmental factor. It is concluded that the toppling deformation of anti-dip layered rock slopes are more sensitive to geological and engineering factors, but less sensitive to environmental factor. In addition, the sensitivity of various factors to the rock toppling deformation is also arranged sequentially as follows: construction, gravitation, rainfall (underground water) and rock structure intensity, etc. (2) 3D deformation study on the anti-dip layered toppling rock slopes Used 3D Distinct Element Method (3DEC) analyzed the 3D deformation characteristic of anti-dip layered rock slops. It can be seen that the toppling characteristics are obvious when the inter-angle between slope direction and layer striking direction is under 20o, when the inter-angle is over 20o and equal or less than 40o，the toppling deformation characteristics decrease sharply with increase of inter-angle, when the inter-angle is over 40o , the slope deformation is not controlled by joints but influenced by other failure mode. Therefore, in order to quantify the toppling characteristics, a differential value of displacement vector angle between layered rock slope and block rock slope is proposed as a key index to distinguish failure model for anti-dip layered rock slopes, and it was applied to study the toppling of the rock slopes at Guangxi Long-tan hydropower station, China. The results indicated that the index was effective and instructive for analyzing the anti-dip layered rock slopes. (3) Stability analysis methods Because of the imperfection of some present slope analysis methods, based on slope failure mode and those three influencing factors, “slope stability entropy” method is defined in this paper, which makes good use of the sensitivity of relational matrix to influencing factors on slope stability and the qualification characteristics for information entropy to the irregularity of slope deformation. By this method, not only the randomness of geologic body on the base of dynamic analysis of slope failure mode is fully concerned, but also it makes the analysis time-saving and simple. Finally, the research findings were used to the engineering example successfully, and rational conclusion has been obtained.