Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Numerical Analyses of Discrete Gust Response for an Aircraft

Based on Navier-Stokes equations and structural and flight dynamic equations of motion, dynamic responses in vertical discrete gust flow perturbation are investigated for a supersonic transport model. A tightly coupled method was developed by subiterations between aerodynamic equations and dynamic equations of motion. First, under the assumption of rigid-body and single freedom of motion in the vertical plunging, the results of a direct-coupling method are compared with the results of quasi-steady model method. Then, gust responses for the one-minus-cosine gust profile arc analyzed with two freedoms of motion in plunging and pitching for the airplane configurations with and without the consideration of structural deformation.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

Component Assembling Model and Its Application to Quasi-Brittle Damage

Potential energy can be approximated by ‘‘pair-functional’’ potentials which is composed of pair potentials and embedding energy. Pair potentials are grouped according to discrete directions of atomic bonds such that each group is represented by an orientational component. Meanwhile, another kind of component, the volumetric one is derived from embedding energy. Damage and fracture are the changing and breaking of atomic bonds at the most fundamental level and have been reflected by the changing of these components’ properties. Therefore, material is treated as a component assembly, and its constitutive equations are formed by means of assembling these two kinds of components’ response functions. This material model is referred to as the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness and intrinsic induced anisotropy, etc.

A finite element method for cracked components of structures

In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element model of the global platform for the FE calculation of the structure in any environmental conditions. The stiffness matrix equation of the cracked component is derived by use of the finite variation principle and fracture mechanics. The equivalent parameters defining the element that simulates the cracked component are mathematically presented, and can be easily used for the FE calculation of large scale cracked structures together with any finite element program. The theories developed are validated by both lab tests and numerical calculations, and applied to the evaluation of crack effect on the strength of a fixed platform and a self-elevating drilling rig.

Characterization of a Laser-Discrete Quenched Steel Substrate/Chromium System by Dissolving Coatings

A laser-discrete quenched steel (LDQS) substrate/as-deposited chromium (top high-contraction (HC) and underlying low-contraction (LC) chromium) system was investigated by dissolving coatings in order to reveal the mechanism that the service life of the coated parts is largely improved using the hybrid technique of laser pre-quenching plus chromium post-depositing. It was found that the surface characteristics of the substrate, LC and HC chromium layer can be simultaneously revealed owing to the dissolution edge effect of chromium coatings. Moreover, the periodical gradient morphologies of the LDQS substrate are clearly shown: the surfaces of laser transformation-hardened regions are rather smooth; a lot of fine micro-holes exist in the transition zones; there are many micro-dimples in the original substrate. Furthermore, the novel method of dissolving coatings with sharp interfaces may be used to reveal the structural features of a substrate/coating system, explore the effect of the substrate on the initial microstructure and morphologies of coatings, and check the quality of the coated-parts.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

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.

Size Effects in Micro-bending Deformation of MEMS Devices Based on the Discrete Dislocation Theory

In the present research, the discrete dislocation theory is used to analyze the size effect phenomena for the MEMS devices undergoing micro-bending load. A consistent result with the experimental one in literature is obtained. In order to check the effectiveness to use the discrete dislocation theory in predicting the size effect, both the basic version theory and the updated one are adopted simultaneously. The normalized stress-strain relations of the material are obtained for different plate thickness or for different obstacle density. The prediction results are compared with experimental results.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Degradation failure features of chromium-plated gun barrels with a laser-discrete-quenched substrate

The effect of substrate laser-discrete quenching on the degradation failure of chromium-plated gun barrels was metallurgically investigated. The results show that substrate laser-discrete quenching changes the failure patterns of chromium coatings during firing, and some periodic through-thickness cracks in the fired chromium coatings are justly located at original substrate zones between two adjacent laser-quenched tracks. Moreover, chromium coatings and the laser-quenched zones on the substrate are simultaneously degraded in microstructure and property during firing. Furthermore, the periodic structure of the laser-discrete-quenched steel (LDQS) substrate near the breech remains after firing, and the hardness of the fired laser-quenched zones is still higher than that of original substrates. The specific failure features were utilized to illustrate the mechanism of the extended service life of chromium-plated gun barrels with the LDQS substrate. (c) 2007 Elsevier B.V All rights reserved.

Analysis of critical excavation depth for a jointed rock slope using a face-to-face discrete element method

The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D-s and the natural slope angle alpha, the joint inclination angle theta as well as the strength parameters of the joints c(r) ,phi(r) is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Microfluidic effects of transporting signaling components in cell coculture chips

A theoretical model has been developed to investigate the microfluidic transport of the signaling chemicals in the cell coculture chips. Using an epidermal growth factor (EGF)-like growth factor as the sample chemical, the effects of velocities and channel geometry were studied for the continuous-flow microchannel bioreactors. It is found that different perfusion velocities must be applied in the parallel channels to facilitate the communication, i.e., transport of the signaling component, between the coculture channels. Such communication occurs in a unidirectional way because the signaling chemicals can only flow from the high velocity area to the low velocity area. Moreover, the effect of the transport of the signaling component between the coculture channels on the growth of the monolayer cells and the multicellular tumor spheroid (MTS) in the continuous-flow coculture environment were simulated using 3D models. The numerical results demonstrated that the concentration gradients will induce the heterogeneous growth of the cells and the MTSs, which should be taken into account in designing the continuous-flow perfusion bioreactor for the cell coculture research.