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.

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001 a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.

Weibull modulus for diverse strength due to sample-specificity

By sample specificity it is meant that specimens with the same nominal material parameters and tested under the same environmental conditions may exhibit different behavior with diversified strength. Such an effect has been widely observed in the testing of material failure and is usually attributed to the heterogeneity of material at the mesoscopic level. The degree with which mesoscopic heterogeneity affects macroscopic failure is still not clear. Recently, the problem has been examined by making use of statistical ensemble evolution of dynamical system and the mesoscopic stress re-distribution model (SRD). Sample specificity was observed for non-global mean stress field models, such as the duster mean field model, stress concentration at tip of microdamage, etc. Certain heterogeneity of microdamage could be sensitive to particular SRD leading to domino type of coalescence. Such an effect could start from the microdamage heterogeneity and then be magnified to other scale levels. This trans-scale sensitivity is the origin of sample specificity. The sample specificity leads to a failure probability Phi (N) with a transitional region 0 < (N) < 1, so that globally stable and catastrophic modes could co-exist. It is found that the scatter in strength can fit the Weibull distribution very well. Hence, the Weibull modulus is indicative of sample specificity. Numerical results obtained from the SRD for different non-global mean stress fields show that Weibull modulus increases with increasing sample span and influence region of microdamage.

Combined effect of surface tension, gravity and van der Waals force induced by a non-contact probe tip on the shape of liquid surface

Aiming at understanding how a liquid film on a substrate affects the atomic force microscopic image in experiments, we present an analytical representation of the shape of liquid surface under van der Waals interaction induced by a non-contact probe tip. The analytical expression shows good consistence with the corresponding numerical results. According to the expression, we find that the vertical scale of the liquid dome is mainly governed by a combination of van der Waals force, surface tension and probe tip radius, and is weekly related to gravity. However, its horizontal extension is determined by the capillary length.

Analysis of two-dimensional finite solids with microcracks

A general method is presented for solving the plane elasticity problem of finite plates with multiple microcracks. The method directly accounts for the interactions between different microcracks and the effect of outer boundary of a finite plate. Analysis is based on a superposition scheme and series expansions of the complex potentials. By using the traction-free conditions on each crack surface and resultant forces relations along outer boundary, a set of governing equations is formulated. The governing equations are solved numerically on the basis of a boundary collocation procedure. The effective Young's moduli for randomly oriented cracks and parallel cracks are evaluated for rectangular plates with microcracks. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method provides a direct and efficient approach to deal with finite solids containing multiple microcracks.

Relationships Between Initial Unloading Slope, Contact Depth, and Mechanical Properties for Conical Indentation in Linear Viscoelastic Solids

Using analytical and finite element modeling, we examine the relationships between initial unloading slope, contact depth, and mechanical properties for spherical indentation in viscoelastic solids with either displacement or load as the independent variable. We then investigate whether the Oliver-Pharr method for determining the contact depth and contact radius, originally proposed for indentation in elastic and elastic-plastic solids, is applicable to spherical indentation in viscoelastic solids. Finally, the analytical and numerical results are used to answer questions raised in recent literature about measuring viscoelastic properties from instrumented spherical indentation experiments.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

Cross-relaxation dynamics on anomalous saturation processes in low pressure supersonic cw HF chemical laser amplifier

It is assumed that both translational and rotational nonequilibrium cross-relaxations play a role simultaneoulsy in low pressure supersonic cw HF chemical laser amplifier. For two-type models of gas flow medium with laminar and turbulent flow diffusion mixing, the expressions of saturated gain spectrum are derived respectively, and the numerical calculations are performed as well. The numerical results show that turbulent flow diffusion mixing model is in the best agreement with the experimental result.

Reliable Validation Based on Optical Flow Visualization for CFD Simulations

A reliable validation based on the optical flow visualization for numerical simulations of complex flowfields is addressed in this paper. Several test cases, including two-dimensional, axisymmetric and three-dimensional flowfields, were presented to demonstrate the effectiveness of the validation and gain credibility of numerical solutions of complex flowfields. In the validation, images of these flowfields were constructed from numerical results based on the principle of the optical flow visualization, and compared directly with experimental interferograms. Because both experimental and numerical results are of identical physical representation, the agreement between them can be evaluated effectively by examining flow structures as well as checking discrepancies in density. The study shows that the reliable validation can be achieved by using the direct comparison between numerical and experiment results without any loss of accuracy in either of them.

Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Shear deformation analysis on laser-induced reverse bulging and plugging

A new kind of failure mode is observed in circular brass foils in which their peripheries are fixed and their surfaces are subjected to a long pulsed laser over a central region. The failure is classified into three stages; they are referred to as thermal bulging, localized shear deformation and perforation by plugging. A distinct feature of the failure mode is that bulging and plugging occurred in the direction opposite to the incident laser beam. To study the failure mode, we investigate the non-linear response of heated, non-homogeneous circular plates. Based on the large deflection equations of Berger [J. Appl. Mech. 22 (3), 465-472 (1965)], Ohnabe and Mizuguchi [Int. J. Non-Linear Mech. 28 (4), 365-372 (1993)] and the parabolic shear deformation theory of Bhimaraddi and Stevens [J. Appl. Mech. 51 (1), 195-198 (1984)], we have derived new coupled governing equations of shear deformation and deflection. The new equations are solved, for the plate with a clamped edge, by the Galerkin and iterative methods. The numerical results for the shear deformation distribution are in good agreement with the experimental observation.

An incremental micromechanical scheme for nonlinear particulate composites

A general incremental micromechanical scheme for the nonlinear behavior of particulate composites is presented in this paper. The advantage of this scheme is that it can reflect partly the effects of the third invariant of the stress on the overall mechanical behavior of nonlinear composites. The difficulty involved is the determination of the effective compliance tensors of the anisotropic multiphase composites. This is completed by making use of the generalized self-consistent Mori-Tanaka method which was recently developed by Dai et al. (Polymer Composites 19(1998) 506-513; Acta Mechanica Solida 18 (1998) 199-208). Comparison with existing theoretical and numerical results demonstrates that the present incremental scheme is quite satisfactory. Based on this incremental scheme, the overall mechanical behavior of a hard-particle reinforced metal matrix composite with progressive particle debonding damage is investigated.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

An Improved Ce/Se Scheme And Its Application To Detonation Propagation

An improved two-dimensional space-time conservation element and solution element ( CE/ SE) method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.