Thesis (Ph.D.) -- Virginia Commonwealth University, 2005. / Title from title-page of electronic thesis. Prepared for: Dept. of Mechanical Engineering. Bibliography: p. 112-123.

With the developments in computer technology, the numerical methods are used widely in geotechnical engineering. Finite element and finite difference are the most common methods used to simulate the behavior of soil and rock. Although the reliability of these methods are proven in several fields of application over the years, they are not equally satisfactory in every case and require sophisticated constitutive relations to model the discontinuous behavior of geomaterials since they assume the material is continuum or the location of discontinuum is predictable. The Discrete Element Method (DEM) has an intensive advantage to simulate discontinuity. This method is relatively new and still under development, yet it is estimated that it will replace of the continuum methods largely in geomechanics in the near feature. In this thesis, the theory and background of discrete element method are introduced, and its applicability in bearing capacity calculation of shallow foundations is investigated. The results obtained from discrete element simulation of bearing capacity are compared with finite element analysis and analytical methods. It is concluded that the DEM is a promising numerical analysis method but still have some shortcomings in geomechanical applications.

The effect of in-situ or redistributed stress on solute transport in fractured rocks is one of the major concerns for many subsurface engineering problems. However, it remains poorly understood due to the difficulties in experiments and numerical modeling. The main aim of this thesis is to systematically investigate the influences of stress on solute transport in fractured rocks, at scales of single fractures and fracture networks, respectively. For a single fracture embedded in a porous rock matrix, a closed-form solution was derived for modeling the coupled stress-flow-transport processes without considering damage on the fracture surfaces. Afterwards, a retardation coefficient model was developed to consider the influences of damage of the fracture surfaces during shear processes on the solute sorption. Integrated with particle mechanics models, a numerical procedure was proposed to investigate the effects of gouge generation and microcrack development in the damaged zones of fracture on the solute retardation in single fractures. The results show that fracture aperture changes have a significant influence on the solute concentration distribution and residence time. Under compression, the decreasing matrix porosity can slightly increase the solute concentration. The shear process can increase the solute retardation coefficient by offering more sorption surfaces in the fracture due to gouge generation, microcracking and gouge crushing. To study the stress effects on solute transport in fracture systems, a hybrid approach combing the discrete element method for stress-flow simulations and a particle tracking algorithm for solute transport was developed for two-dimensional irregular discrete fracture network models. Advection, hydrodynamic dispersion and matrix diffusion in single fractures were considered. The particle migration paths were tracked first by following the flowing fluid (advection), and then the hydrodynamic dispersion and matrix diffusion were considered using statistic methods. The numerical results show an important impact of stress on the solute transport, by changing the solute residence time, distribution and travel paths. The equivalent dispersion coefficient is scale dependent in an asymptotic or exponential form without stress applied or under isotropic compression conditions. Matrix diffusion plays a dominant role in solute transport when the hydraulic gradient is small. Outstanding issues and main scientific achievements are also discussed. / QC 20111011

Coatings utilizing the piezospectroscopic (PS) effect of alpha alumina could enable on the fly stress sensing for structural health monitoring applications. While the PS effect has been historically utilized in several applications, here by distributing the photo-luminescent material in nanoparticle form within a matrix, a stress sensing coating is created. Parallel to developing PS coatings for stress sensing, the multi-scale mechanics associated with the observed PS response of nanocomposites and their coatings has been applied to give material property measurements, providing an understanding of particle reinforced composite behavior. Understanding the nanoparticle-coating-substrate mechanics is essential to interpreting the spectral shifts for stress sensing of structures. In the past, methods to experimentally measure the mechanics of these embedded nano inclusions have been limited, and much of the design of these composites depend on computational modeling and bulk response from mechanical testing. The PS properties of Chromium doped alumina allow for embedded inclusion mechanics to be revisited with unique experimental setups that probe the particles state of stress under applied load to the composite. These experimental investigations of particle mechanics will be compared to the Eshelby theory and its derivative theories in addition to the nanocomposite coating mechanics. This work discovers that simple nanoparticle load transfer theories are adequate for predicting PS properties in an intermediate volume fraction range. With fundamentals of PS nanocomposites established, the approach was applied to selected experiments to prove its validity. In general it was observed that the elastic modulus values calculated from the PS response were similar to that observed from macroscale strain measurements such as a strain gage. When simple damage models were applied to monitor the elastic modulus, it was observed that the rate of decay for the elastic modulus was much higher for the PS measurements than for the strain gage. A novel experiment including high resolution PS maps with secondary strain maps from digital image correlation is reviewed on an open hole tension, composite coupon. The two complementary measurements allow for a unique PS response for every location around the hole with a spatial resolution of 400 microns. Progression of intermediate damage mechanisms was observed before digital image correlation indicated them. Using the PS nanocomposite model, elastic modulus values were calculated. Introducing an elastic degradation model with some plastic deformation allows for estimation of material properties during the progression of failure. This work is part of a continuing effort to understand the mechanics of a stress sensing PS coating. The mechanics were then applied to various experimental data that provided elastic property calculations with high resolution. The significance is in the experimental capture of stress transfer in particulate composites. These findings pave the way for the development of high resolution stress-sensing coatings.

The research presented in this work addresses open questions regarding (i) the fundamental understanding of powder compaction, and (ii) the complex mechanical response of particle-binder composites under large deformations. This work thus benefits a broad range of industries, from the pharmaceutical industry and its recent efforts on continuous manufacturing of solid tablets, to the defense and energy industries and the recurrent need to predict the performance of energetic materials. Powder compacts and particle-binder composites are essentially confined particulate systems with significant heterogeneity at the meso (particle) scale. While particle mechanics strategies for modeling evolution of mesoscale microstructure during powder compaction depend on the employed contact formulation to accurately predict macroscopic quantities like punch and die wall pressures, modeling of highly nonlinear, strain-path dependent macroscopic response without a distinctive yield surface, typical of particle-binder composites, requires proper constitutive modeling of these complex deformation mechanisms. Moreover, continued loading of particle-binder composites over their operational life may introduce significant undesirable changes to their microstructure and mechanical properties. These challenges are addressed with a combined effort on theoretical, modeling and experimental fronts, namely, (a) novel contact formulations for elasto-plastic particles under high levels of confinement, (b) a multi-scale experimental procedure for assessing changes in microstructure and mechanical behavior of particle-binder composites due to cyclic loading and time-recovery, and (c) a finite strain nonlinear elastic, endochronic plastic constitutive formulation for particle-binder composites.