Papers by Author: Shan Yong Wang

Abstract: A series of numerical model tests were performed to investigate the behaviour of the anisotropic rock surrounding circular excavations under high confining pressures. The aim was to provide information on the formation of fractures and failure around deep level rock tunnels under controlled conditions. Solid cubes containing a circular hole were confined to a vertical pressure with same as the confinement in the horizontal directions. In this modeling, the inhomogeneous rock is generated by using Weibull parameters which are related to the microstructural properties determined by crack size distribution and grain size. The fracture angle is assumed to be 45o. The observed failure zone around the excavation was simulated using both the maximum tensile strain criterion and Mohr-Coulomb criterion respectively (as the damage threshold). And RFPA (Realistic Failure Process Analysis) code was used as the calculating tool in this modelling, three opening modes are simulated and compared. Computational model predictions that include crack propagation and failure modes of rock show a good agreement with those of the observation in site. It is pointed out that the damage evolution of EDZ strongly depends on the inhomogeneous, the excavation mode, anisotropic property, and the various loading conditions. Concerning the existence of a weak plane, the amount of displacement at the side wall of the tunnel was quite large, since the shear deformation occurred in EDZ. The model is implemented in RFPA code and is able to represent the change in fracture patterns between the solid and jointed parts. This provides confidence for the application of the numerical model to the design of rock tunnels at great depth.

Abstract: An Elastic Plastic-Damage (EPD) model is developed to model the softening behaviour of the cement-soil admixture based on continuous damage mechanics. The softening behaviour is considered to be characteristic outcome of the material degradation due to damage in material. Material degradation is modelled by reducing progressively the stiffness and yield stress of the material when the damage variable has attained a critical index. The basic equations of the model are derived and presented. A Fortran program for this model has been developed and implemented into a finite element code ABAQUS. In order to evaluate the applicability of this model, several unconfined compression tests are simulated using ABAQUS with this model. The computed results are compared with measured data and good agreement is achieved.

Abstract: In this paper, a Material Failure Process Analysis code (MFPA2D) was employed to investigate the interaction of end effect zone of specimen with the wing crack propagation inside the brittle specimen containing pre-existing flaws under uniaxial compression comparing with the experimental results. The numerical results show that the shorter the distance between the pre-existing flaw and the specimen's end , the slower the crack propagation process and the shorter wing propagation length is , and vice versa. In addition, the end effect zone was also influenced by the wing crack propagation.

Abstract: Soil nailing is a widely used technique for stabilizing slopes and excavations. In all current design methods, the nail-soil interface shear strength, that is, the pull-out resistance of a soil nail is an important parameter which controls the design and safety assessment of the soil nailing system. The pressure grouting is a cost effective method for increasing the soil nail pull-out resistance and in turn improving the performance of the nailed structure. In this paper, a three dimensional (3-D) finite element (FE) model for pull-out tests is established and verified by comparing simulated results with measured data. This model is then used to simulate the effect of grouting pressure on the soil nail pull-out resistance.

Abstract: In numerical simulation of engineering problems, it is important to properly simulate the interface between two adjacent parts of the model. In finite element method, there are generally three methods for simulating interface problems: interface element method, surface based contact method and the method by using a thin layer of continuum elements. In this paper, simulation of interface problems is conducted using continuum elements and surface based contact methods. The results from each method are presented and compared with each other.

Abstract: One small-scale physical model test on the PVD (Prefabricated Vertical Drain) treated Hong Kong marine clay was simulated using finite element method (FEM) in this study. A User MATerial (UMAT) subroutine describing an Elastic Visco-Plastic (EVP) constitutive model was developed and incorporated into one commercial finite element code ABAQUS. A degressive permeability of the PVD strip was included to consider variations of its permeability during the consolidation process. The UMAT and the adopted reducing technique were demonstrated to be effective by good agreement between the observed consolidation settlement and excess pore water pressures and the simulated ones.

Abstract: In this paper, a coupled thermal-mechanical-damage model, Material Failure Process Analysis for Thermo code (abbreviated as MFPA-thermo), was applied to investigate the formation, extension and coalescence of cracks in FRCs, caused by the thermal mismatch of the matrix and the particles under uniform temperature variations. The effects of the thermal mismatch between the matrix and fibers on the stress distribution and crack development were also numerically studied. The influences of the material heterogeneity, the failure patterns of FRCs at varied temperatures are simulated and compared with the experimental results in the present paper. The results show that the mechanisms of thermal damage and fracture of the composite remarkedably depend on the difference between the coefficients of thermal expansion of the fibers and the matrix on a meso-scale. Meanwhile, the simulations indicate that the thermal cracking of the FRCs at uniform varied temperatures is an evolution process from diffused damage, nucleation, and finally linkage of cracks.

Abstract: Rock and concrete are typical heterogeneous material that the meso-scale heterogeneity may have a significant effect on their macro-scale mechanical responses. In this work, a digital image-based (DIB) technique is employed to characterize and quantify the heterogeneity of concrete, and the obtained data is directly imported into a numerical code named RFPA (Rock Failure Process Analysis) to study the effect of heterogeneity on the failure process of concrete. The upgraded RFPA is capable to simulate the progressive failure of brittle materials such as rock and concrete, representing both the growth of existing fractures and the formation of new fractures, obviating the need to identify crack tips and their interaction explicitly. The simulated results are in reasonable agreement with experimental measurements and phenomenological observations reported in previous studies.

Abstract: Many stiff clays forming part of natural slopes and earth dams exist in the fissured state. When these cracks are subjected to gravity induced normal and shear stresses they may propagate. The present discussion presents a numerical method to study the propagation direction of cracks under stress fields similar to those found in the field. Not only did the results on one crack propagation direction obtained from the numerical method and the analytical results agree well, but numerical results have been used to investigate the mechanisms of the whole process of two horizontal cracks initiation and propagation and coalescence in stiff soils.

Abstract: Based on cusp-type catastrophe theory, a sample rock-rock (hypocenter surrounding the rock) model for studying the pillar rockburst mechanism is presented in this paper. It is expounded theoretically that the stiffness ratio, K, of the roof and floor to the pillar plays an important role in the outbreak of instability. Using a newly developed numerical code, RFPA2D, the progressive failure process and associated microseismic behavior of the twin rock samples are simulated. The numerically simulated results also confirm that a soft roof and floor promotes an unstable failure or collapse of pillars. Additionally, the simulated results reproduced the deformation jump and the energy release that occur during a pillar rockburst. It is demonstrated that the proposed model properly simulates the pillar failure process.