Papers by Author: Shui Lin Wang

Abstract: The purpose of this paper is to study the shear behavior of rock specimens containing joints with various distribution forms. Two sets of specimens are simulated by the rock failure process analysis code (RFPA2D). The friction-sliding failure pattern occurs with the lower undulation angle specimen, and the failure pattern turns to be tensile-shear failure mode gradually with the increase of undulation angle. The specimen possesses the highest peak shear load when the undulation angle is about 30º. And joint rock shear character also deteriorates with the increase of weak interlayer thickness. In the intermittent joint model, the unified connection ratio specimen’s peak shear load increases with rock bridge amount, and the multi-joint mode is beneficial to keep rock mass shear stiffness. This study comes to meaningful results to the expansion of joint rock strength evolution law with various joint distribution forms.

Abstract: the inverse analysis to material parameters is often translated into an optimization for an objective function, based on the correlation between the material parameters and the foregone information. But mostly because of the non-linear correlation, a good optimization algorithm with the capabilities to avoid being trapped by local optima is required during the process of optimization. So the present paper proposes a new global optimization algorithm, which couples the dynamic canonical descent algorithm and the improved Powell’s algorithm. The high efficiency of the new algorithm is shown on four known problems classically for testing optimization algorithms and finally, in the non-linear inverse analysis, the new algorithm is used for optimizing an objective function to get material parameters rightly.

Abstract: In geotechnical engineering, based on the theory of inverse analysis of displacement, the problem for identification of material parameters can be transformed into an optimization problem. Commonly, because of the non-linear relationship between the identified parameters and the displacement, the objective function bears the multimodal characteristic in the variable space. So to solve better the multimodal characteristic in the non-linear inverse analysis, a new global optimization algorithm, which integrates the dynamic descent algorithm and the modified BFGS (Brogden-Fletcher-Goldfrab-Shanno) algorithm, is proposed. Five typical multimodal functions in the variable space are tested to prove that the new proposed algorithm can quickly converge to the best point with few function evaluations. In the practical application, the new algorithm is employed to identify the Young’s modulus of four different materials. The results of the identification further show that the new proposed algorithm is a very highly efficient and robust one.

Abstract: This study is to present a numerical investigation on fragmentation and perforation of concrete slab by hard projectile using discrete particle approaches. Discrete Element Method (DEM) and Discontinuous Deformation Analysis (DDA), the two representative discrete particle approaches, are employed to simulate a normal perforation of concrete slab by a hard ogival-nose shaped projectile, and the phenomena of spalling, plugging and scabbing are reproduced.

Abstract: There are a large number of cracks, joints and layers with different scales and orientations in the rock masses. With the structural formation in it, the rock mass is unlike the isotropic and homogenous materials in physics and mechanics characteristics, and the failure of the rock mass is controlled by those discontinuities. This paper studies the stability of a jointed rock slope with bedding-planes. Universal Distinct Element Code (UDEC) is used to carry out the numerical analysis. The factor of safety (FOS) and failure mode are obtained by strength reduction technique before and after the steel bolts are installed in the slope. In addition, the forces along the bolts indicate that they reaches the maximum value in the potential slip surface of the slope.

Abstract: There are several different yield surfaces of Drucker-Prager yield criterion which corresponds to Mohr-Coulomb yield criterion in principal stress space. The different yield surfaces are determined by parameters in Drucker-Prager criterion. The influence of the different parameters on computational results is discussed in the paper, and the use of the equivalent Drucker-Prager criterion to Mohr-Coulomb criterion is suggested when elastoplastic analysis is performed in engineering problems.

Abstract: A numerical technique based on using manifold elements in finite element method, for modeling propagation of arbitrary cracks in solids, is described. When the region with crack(s)is subjected to external loading and the crack(s) starts to extend, the crack growth may intersect boundaries of nearby finite elements. Those intersected finite elements are replaced by manifold elements. The technique, by which the initial finite element mesh can be kept unchanged during the processes of crack propagation, is called manifold elements in finite element method. The crack growth is governed by the theories of linear elastic fracture mechanics. The stress intensity factors are computed by a contour integral technique and crack trajectory is determined by applying the maximum tangential stress criterion. Finally, test examples are given to verify the new method and the predicted trajectories are compared to experimentally obtained crack growth paths with good agreements

Abstract: A series of formulas about two-parameter parabolic Mohr strength criterion(2-PP Mohr criterion) are derived. Based on the results of uniaxial tension and uniaxial compression tests, the parameters involved in the criterion can be easily determined, then the criterion in terms of the major principal stress and the minor principal stress is derived, and the damage pattern is also discussed. At last, the formulas about the rupture angle and the friction angle are presented, and their relationship is also given. 2-PP Mohr criterion can describe not only shear but also tensile failure. In this criterion the ratio of the uniaxial compression strength and the uniaxial tension strength is not confined as in Griffith criterion. The formula about the rupture angle provides steady theoretical foundation for determining the direction of crack faces and damage patterns in the computation of macro crack propagation. In fact, Griffith criterion is only a special case of the two-parameter parabolic Mohr strength criterion proposed in this present paper.

Abstract: In finite element method, the order of complete polynomial of the interpolation function is related to the number of nodes in the element. This paper presents a four-node quadrilateral element with quadratic function on it. The presented displacement functions maintain C0-continuity. Meanwhile, the element stiffness matrix is derived from the displacement functions. Test problems show that high accuracy can be achieved by the use of the new displacement function on the element.