This dissertation focuses mainly on three topics: (1) mixed-mode branching and segmentation of hydraulic fractures in brittle materials, (2) hydraulic fracture propagation in particulate materials, and (3) hydraulic fracturing in water flooding conditions. Mixed-mode loading is one of the primary causes of fracture branching and segmentation in brittle materials. We conducted the first laboratory experiments on the mixed mode I+III hydraulic fracturing. We found that a KIII/KI ratio as small as ~1% is sufficient for fracture front segmentation. In reality, such a small mode III component is always expected, for example, due to the small deviations of the fracture shape from planar. Thus, we concluded that fracture segmentation is likely to accompany growth of most, if not all, real hydraulic fractures. We also proposed a theoretical model that captures the main features of experimental observations and indicates the importance of the hydraulic effect of segmentation. Particulate materials often exhibit pronounced non-linear behavior and yielding even at relatively small loads. In order to adequately describe hydraulic fracturing in particulate materials with low or no cohesion, plasticity at the crack tip must be explicitly considered. We investigated the shear band mechanism of strain localization at the fracture front. This mechanism takes into account the fact that cohesionless material can not bear tension, and is in compression everywhere, including near the fracture front. To verify the shear band hypothesis, we conducted numerical simulations of the plastic deformation at the tip of a fracture in particulate material with strain softening. Our model describes the shear bands by properly placed and oriented dislocations. The model results are consistent with experimental observations. Water flooding, which in certain important cases, can result in processes resembly hydraulic fracturing by a low-viscosity fluid with extremely high leak-off. It is difficult to simulate this process in the laboratory. To investigate the fracture initiation mechanism in water flooding conditions, we conducted a numerical simulation of fluid injection into particulate material by using the discrete element code PFC2D. We also considered an analytical model of cavity initiation based on the fluidization mechanism. The estimates given by this model fit remarkably well with the numerical simulation results.

The soil-structure interface is fundamental to the performance of many geotechnical engineering systems; including penetration test devices, deep foundations, and retaining structures. In geotechnical engineering structures, the counterface may range from a polymer in the case of a geosynthetically reinforced earth retaining structure to steel for cone penetration testing or pile foundations. Interface strength is affected by many factors, among which surface roughness is the most dominant. To date common practice has been to characterize counterface surface roughness by a roughness parameter based on only its spatial properties and soil roughness separately by various incompatible means resulting in two roughness values unrelated to each other. The vast number of analyzing methods and developed parameters reveal the general confusion regarding this concept. Rather than analyzing the particulate and continuum media separately, it is compulsory to coalesce the analysis and quantify the relative nature of interface behavior. This can be accomplished by examining the particulate and continuum media through the same powerful tools. The motive of this study is to develop a unified approach to determining the index properties of particles and surfaces in a particle-surface interface. This is accomplished by examining several particle shape and surface roughness parameters in terms of their ability to uniquely describe and distinguish particulate medium and continuum roughness, respectively. In this study, surfaces are analyzed as derived particles by wrapping surface profiles and particles are evaluated as derived surfaces via unrolling particle outlines. In addition, particle shape parameters are modified to allow surface roughness analysis and surface roughness parameters are modified to characterize particle shape. A unified approach for particulate shape and continuum roughness would ultimately lead to a better understanding of micro-scale interaction mechanism and better quantification of macro-scale mobilized resistance for soil and engineering surface interaction.

Under loading conditions found in many geotechnical structures, it is common to observe failure in zones of high localized strain called shear bands. Existing models predict these localizations, but provide little insight into the micromechanics within the shear bands. This research captures the variation in microstructure inside and outside of shear bands that were formed in laboratory plane strain and two-dimensional discrete element method (DEM) biaxial compression experiments. Plane strain compression tests were conducted on dry specimens of Ottawa 20-30 sand to calibrate the device, assess global response repeatability, and develop a procedure to quantitatively define the onset of localization. A new methodology was employed to quantify and correct for the additional stresses imparted by the confining membrane in the vicinity of the shear band. Unsheared and sheared specimens of varying dilatancy were solidified using a two-stage resin impregnation procedure. DEM tests were performed using an innovative servo-controlled flexible lateral confinement algorithm to provide additional insights into laboratory results. The solidified specimens were sectioned and the resulting surfaces prepared for microstructure observation using bright field microscopy and morphological analysis. Local void ratio distributions and their statistical properties were determined and compared. Microstructural parameters for subregions in a grid pattern and along predefined inclined zones were also calculated. Virtual surfaces parallel to the shear band were identified and their roughnesses assessed. Similar calculations were performed on the DEM simulations at varying strain levels to characterize the evolution of microstructure with increasing strain. The various observations showed that the mean, standard deviation, and entropy of the local void ratio distributions all increased with increasing strain levels, particularly within regions of high local strains. These results indicate that disorder increases within a shear band and that the soil within the shear band does not adhere to the classical concept of critical state, but reaches a terminal void ratio that is largely a function of initial void ratio. Furthermore, there appears to be a transition zone between the far field and the fully formed shear block, as opposed to an abrupt delineation as traditionally inferred.