Browsing by Author "Martin, James R. II"
Now showing 1 - 20 of 25
Results Per Page
Sort Options
- Applications of Roll-Along Electrical Resistivity Surveying in Conjunction with Other Geophysical Methods for Engineering and Environmental Site CharacterizationSayer, Suzanne (Virginia Tech, 1996-04-29)Roll-along electrical resistivity surveying was used with seismic refraction, magnetometer and gravity surveying in geophysical characterization of sites with a specific environmental or engineering problem. Three examples are presented where resistivity surveying provided vital constraints on acquisition and interpretation of other data in chaotic terrane. A commercially resistivity meter was used with prototype equipment designed, assembled, and tested at Virginia Tech. The equipment included a multiconductor cable consisting of interchangeable segments and a circuit allowing selection of numerous electrode configurations. The Sinking Creek Landfill, a 10-acre site, was used for disposal of municipal waste in the early 1970’s. Roll-along resistivity proved to be the most useful geophysical tool in ascertain its internal structure. Wenner configuration resistivity data, sensitive to both conductive leachate and ferrous metals, showed trenches within the landfill displayed in profile. Magnetic field measurements revealed anomalies over some trenches suggesting a method for discriminating between ferrous metal and leachate. Results of a resistivity survey can help planners of a refraction survey avoid low velocity “blind” layers. The Mid County Landfill borrow area is a 26 acre site situated within the Max Meadows Breccia and used for cover material for an adjacent landfill, The engineering problems were to measure the volume of rippable material, but travel time data were somewhat ambiguous. The refraction data interpreted using a) conventional 3-layer analysis b) horizontal 3-layer analysis of single shots, and c)continuous velocity gradient analysis of single end shots were compared with auger refusal depth. The single end horizontal analysis matched auger refusal depths most closely. Roll-along resistivity pseudo-sections made along the refraction lines proved to be effective for qualitatively imaging pinnacles and megaclasts. Excavation of fill material from a 75 acre river terrace in Pembroke exposed an antiform cut by high angle, near surface faults. Geophysical characterization was undertaken to determine the thickness of the alluvial deposit, and the relationship of the faults with structures in the underlying bedrock. Seismic refraction showed the terrace was as much as 134 feet thick. Resistivity pseudosections revealed a resistivity anomaly associated with the graben could be detected for a horizontal distance of several hundred feet. A gravity gradient paralleling the resistivity anomaly extends the feature more than 1000 ft from the exposed structure. Tenuous evidence of a bedrock escarpment beneath the near surface structure is found in a combination of seismic refraction, gravity, and electrical resistivity data. Roll-along resistivity has proved to be key to geophysical interpretation of these three areas. Images displayed on pseudosections reveal lateral inhomogeneity more clearly than could be discerned from seismic, gravity and magnetic data. Roll-along resistivity data can provide information for efficient siting of additional geophysical studies.
- Assessing the Seismic Hazard in Charleston, South Carolina: Comparisons Among Statistical ModelsStudent, Heather H. (Virginia Tech, 1997-01-27)Seismic hazard calculations for sites in eastern North America have traditionally assumed a Poisson process to describe the temporal behavior of earthquakes and have employed the Gutenberg-Richter relationship to define the frequency distribution of earthquake magnitude. For sites in areas where geological information indicates recurrent, large earthquakes, however, such data imply a rate for large events which often exceeds that predicted by the Gutenberg-Richter relationship. One way in which this discrepancy can be reconciled is to assume that the larger events occur as a time-dependent, or renewal, process and possess a "characteristic earthquake" magnitude distribution. The main purpose of this study is to make a quantitative comparison of seismic hazard estimates for Charleston of the influences of 1) the Poisson temporal model assuming the Gutenberg-Richter and characteristic earthquake magnitude recurrence relationships with 2) the renewal temporal model assuming the characteristic magnitude recurrence relationship. Other issues that are examined are the sensitivity of uncertainties of hazard model parameters such as maximum magnitude and seismic source delineation. Probabilistic seismic hazard calculations for the next 50 years were performed at Charleston for all potential seismic sources. The highest estimate of seismic hazard was obtained with the Poisson temporal model and characteristic earthquake recurrence relationship. The lowest hazard was obtained with the renewal temporal model and characteristic magnitude recurrence relationship. The results of this study are in good agreement with hazard estimates for Charleston in the most recent national seismic hazard maps.
- At-rest and compaction-induced lateral earth pressures of moist soilsIshihara, Katsuji (Virginia Tech, 1993-08-01)An instrumented oedometer was designed and constructed for the purpose of investigating at-rest and compaction-induced earth pressures in moist soils. The device has a split oedometer ring, and horizontal stresses are measured using load cells which support one half of the ring. Rapid cyclic loading was applied to compacted soil specimens, using a digital pressure regulator and a computer-based data acquisition system. The performance of the device was validated by performing tests on silicon rubber and Monterey sand.
- Critical height and surface deformation of column-supported embankmentsMcGuire, Michael Patrick (Virginia Tech, 2011-11-01)Column-supported embankments with or without basal geosynthetic reinforcement can be used in soft ground conditions to reduce settlement by transferring the embankment load to the columns through stress redistribution above and below the foundation subgrade level. Column-supported embankments are typically used to accelerate construction and/or protect adjacent facilities from additional settlement. The column elements consist of driven piles or formed-in-place columns that are installed in an array to support a bridging layer or load transfer platform. The bridging layer is constructed to enhance load transfer using several feet of compacted sand or sand and gravel that may include one or more layers of high-strength geotextile or geogrid reinforcement. Mobilization of the mechanisms of load transfer in a column-supported embankment requires some amount of differential settlement between the columns and the embankment as well as between the columns and the foundation soil. When the embankment height is low relative to the clear spacing between columns, there is the risk of poor ride quality due to the reflection of the differential foundation settlement at the surface of the embankment. The minimum embankment height where differential surface settlement does not occur for a particular width and spacing of column is the critical height. The conventional approach is to express critical height as a fixed ratio of the clear span between adjacent columns; however, there is no consensus on what ratio to use and whether a single ratio is applicable to all realistic column arrangements. The primary objective of this research is to improve the understanding of how column-supported embankments deform in response to differential foundation settlement. A bench-scale experimental apparatus was constructed and the equipment, materials, instrumentation, and test procedures are described. The apparatus was able to precisely measure the deformation occurring at the sample surface in response to differential settlement at the base of the sample. Critical heights were determined for five combinations of column diameter and spacing representing a wide range of possible column arrangements. In addition, tests were performed using four different column diameters in a single column configuration with ability to measure the load acting on the column and apply a surcharge pressure to the sample. In total, 183 bench-scale tests were performed over a range of sample heights, sample densities, and reinforcement stiffnesses. Three-dimensional numerical analyses were conducted to model the experiments. The critical heights calculated using the numerical model agreed with the experimental results. The results of the laboratory tests and numerical analyses indicate that critical height depends on the width and spacing of the columns and is not significantly influenced by the density of the embankment fill or the presence of reinforcement. A new method to estimate critical height was developed and validated against extensive case histories as well as experimental studies and numerical analyses performed by others.
- Development of an extended hyperbolic model for concrete-to-soil interfacesGómez, Jesús Emilio (Virginia Tech, 2000-07-19)Placement and compaction of the backfill behind an earth retaining wall may induce a vertical shear force at the soil-to-wall interface. This vertical shear force, or downdrag, is beneficial for the stability of the structure. A significant reduction in construction costs may result if the downdrag is accounted for during design. This potential reduction in costs is particularly interesting in the case of U.S. Army Corps of Engineers lock walls. A simplified procedure is available in the literature for estimating the downdrag force developed at the wall-backfill interface during backfilling of a retaining wall. However, finite element analyses of typical U.S. Army Corps of Engineers lock walls have shown that the magnitude of the downdrag force may decrease during operation of the lock with a rise in the water table in the backfill. They have also shown that pre- and post-construction stress paths followed by interface elements often involve simultaneous changes in shear and normal stresses and unloading-reloading. The hyperbolic formulation for interfaces (Clough and Duncan 1971) is accurate for modeling the interface response in the primary loading stage under constant normal stress. However, it has not been extended to model simultaneous changes in shear and normal stresses or unloading-reloading of the interface. The purpose of this research was to develop an interface model capable of giving accurate predictions of the interface response under field loading conditions, and to implement this model in a finite element program. In order to develop the necessary experimental data, a series of tests were performed on interfaces between concrete and two different types of sand. The tests included initial loading, staged shear, unloading-reloading, and shearing along complex stress paths. An extended hyperbolic model for interfaces was developed based on the results of the tests. The model is based on Clough and Duncan (1971) hyperbolic formulation, which has been extended to model the interface response to a variety of stress paths. Comparisons between model calculations and tests results showed that the model provides accurate estimates of the response of interfaces along complex stress paths. The extended hyperbolic model was implemented in the finite element program SOILSTRUCT-ALPHA, used by the U.S. Army Corps of Engineers for analyses of lock walls. A pilot-scale test was performed in the Instrumented Retaining Wall (IRW) at Virginia Tech that simulated construction and operation of a lock wall. SOILSTRUCT-ALPHA analyses of the IRW provided accurate estimates of the downdrag magnitude throughout inundation of the backfill. It is concluded that the extended hyperbolic model as implemented in SOILSTRUCT-ALPHA is adequate for routine analyses of lock walls.
- Development of Computational Tools for Characterization, Evaluation, and Modification of Strong Ground Motions within a Performance-Based Seismic Design FrameworkSyed, Riaz (Virginia Tech, 2003-12-09)One of the most difficult tasks towards designing earthquake resistant structures is the determination of critical earthquakes. Conceptually, these are the ground motions that would induce the critical response in the structures being designed. The quantification of this concept, however, is not easy. Unlike the linear response of a structure, which can often be obtained by using a single spectrally modified ground acceleration history, the nonlinear response is strongly dependent on the phasing of ground motion and the detailed shape of its spectrum. This necessitates the use of a suite (bin) of ground acceleration histories having phasing and spectral shapes appropriate for the characteristics of the earthquake source, wave propagation path, and site conditions that control the design spectrum. Further, these suites of records may have to be scaled to match the design spectrum over a period range of interest, rotated into strike-normal and strike-parallel directions for near-fault effects, and modified for local site conditions before they can be input into time-domain nonlinear analysis of structures. The generation of these acceleration histories is cumbersome and daunting. This is especially so due to the sheer magnitude of the data processing involved. The purpose of this thesis is the development and documentation of PC-based computational tools (hereinafter called EQTools) to provide a rapid and consistent means towards systematic assembly of representative strong ground motions and their characterization, evaluation, and modification within a performance-based seismic design framework. The application is graphics-intensive and every effort has been made to make it as user-friendly as possible. The application seeks to provide processed data which will help the user address the problem of determination of the critical earthquakes. The various computational tools developed in EQTools facilitate the identification of severity and damage potential of more than 700 components of recorded earthquake ground motions. The application also includes computational tools to estimate the ground motion parameters for different geographical and tectonic environments, and perform one-dimensional linear/nonlinear site response analysis as a means to predict ground surface motions at sites where soft soils overlay the bedrock. While EQTools may be used for professional practice or academic research, the fundamental purpose behind the development of the software is to make available a classroom/laboratory tool that provides a visual basis for learning the principles behind the selection of ground motion histories and their scaling/modification for input into time domain nonlinear (or linear) analysis of structures. EQTools, in association with NONLIN, a Microsoft Windows based application for the dynamic analysis of single- and multi-degree-of-freedom structural systems (Charney, 2003), may be used for learning the concepts of earthquake engineering, particularly as related to structural dynamics, damping, ductility, and energy dissipation.
- Disaggregated Seismic Hazard and the Elastic Input Energy Spectrum: An Approach to Design Earthquake SelectionChapman, Martin C. (Virginia Tech, 1998-06-25)The design earthquake selection problem is fundamentally probabilistic. Disaggregation of a probabilistic model of the seismic hazard offers a rational and objective approach that can identify the most likely earthquake scenario(s) contributing to hazard. An ensemble of time series can be selected on the basis of the modal earthquakes derived from the disaggregation. This gives a useful time-domain realization of the seismic hazard, to the extent that a single motion parameter captures the important time-domain characteristics. A possible limitation to this approach arises because most currently available motion prediction models for peak ground motion or oscillator response are essentially independent of duration, and modal events derived using the peak motions for the analysis may not represent the optimal characterization of the hazard. The elastic input energy spectrum is an alternative to the elastic response spectrum for these types of analyses. The input energy combines the elements of amplitude and duration into a single parameter description of the ground motion that can be readily incorporated into standard probabilistic seismic hazard analysis methodology. This use of the elastic input energy spectrum is examined. Regression analysis is performed using strong motion data from Western North America and consistent data processing procedures for both the absolute input energy equivalent velocity, (Vea), and the elastic pseudo-relative velocity response (PSV) in the frequency range 0.5 to 10 Hz. The results show that the two parameters can be successfully fit with identical functional forms. The dependence of Vea and PSV upon (NEHRP) site classification is virtually identical. The variance of Vea is uniformly less than that of PSV, indicating that Vea can be predicted with slightly less uncertainty as a function of magnitude, distance and site classification. The effects of site class are important at frequencies less than a few Hertz. The regression modeling does not resolve significant effects due to site class at frequencies greater than approximately 5 Hz. Disaggregation of general seismic hazard models using Vea indicates that the modal magnitudes for the higher frequency oscillators tend to be larger, and vary less with oscillator frequency, than those derived using PSV. Insofar as the elastic input energy may be a better parameter for quantifying the damage potential of ground motion, its use in probabilistic seismic hazard analysis could provide an improved means for selecting earthquake scenarios and establishing design earthquakes for many types of engineering analyses.
- Dynamic Analysis of Levee Infrastructure Failure Risk: A Framework for Enhanced Critical Infrastructure ManagementLam, Juan Carlos (Virginia Tech, 2012-05-29)Current models that assess infrastructure failure risk are "linear," and therefore, only consider the direct influence attributed to each factor that defines risk. These models do not consider the undeniable relationships that exist among these parameters. In reality, factors that define risk are interdependent and influence each other in a "non-linear" fashion through feedback effects. Current infrastructure failure risk assessment models are also static, and do not allow infrastructure managers and decision makers to evaluate the impacts over time, especially the long-term impact of risk mitigation actions. Factors that define infrastructure failure risk are in constant change. In a strategic manner, this research proposes a new risk-based infrastructure management framework and supporting system, Risk-Based Dynamic Infrastructure Management System (RiskDIMS), which moves from linear to non-linear risk assessment by applying systems engineering methods and analogs developed to address non-linear complex problems. The approach suggests dynamically integrating principal factors that define infrastructure failure risk using a unique platform that leverages Geospatial Information System services and extensions in an unprecedented manner. RiskDIMS is expected to produce results that are often counterintuitive and unexpected, but aligned to our complex reality, suggesting that the combination of geospatial and temporal analyses is required for sustainable risk-based decision making. To better illustrate the value added of temporal analysis in risk assessment, this study also develops and implements a non-linear dynamic model to simulate the behavior over time of infrastructure failure risk associated with an existing network of levees in New Orleans due to diverse infrastructure management investments. Although, the framework and RiskDIMS are discussed here in the context of levees, the concept applies to other critical infrastructure assets and systems. This research aims to become the foundation for future risk analysis system implementation.
- Effect of Variation of the Systemic Parameters on the Structural Response of Single Degree of Freedom Systems Subjected to Incremental Dynamic AnalysisDe, Samrat (Virginia Tech, 2003-09-01)This thesis presents the results of a study of the effect of variations of systemic parameters on the structural response of single degree of freedom systems subjected to Incremental Dynamic Analysis. The systemic parameters are mass, stiffness, damping, yield strength and geometric stiffness. Each of these parameters was varied one at a time while the other values were kept constant. For each variation of parameters a set of single-record IDA curves was obtained. Five to six ground motions were used for this study to generate the single-record IDA curves. These ground motions were scaled prior to their application on the structure. The scaling factor was based on the spectral acceleration at the fundamental frequency of the structure at 5% of critical damping. The scale factor is affected if the system parameters are changed. An important issue for this study was whether to persist with scaling corresponding to the median value from the range of the values of the parameter or to update the scaling according to the system. Based on some tests using both methods, the median scaling approach was found to be more suitable. The IDA curves for variation of parameters were then investigated to identify any trends that may help in qualitatively predicting the response of a system relative to another system. The response was measured by the peak displacement and the maximum base shear of the system. A clear trend was identified when the damping or the yield strength was varied. However, no definite trend was observed when the material stiffness or the geometric stiffness of the system was varied.
- The Effects Of Non-Plastic and Plastic Fines On The Liquefaction Of Sandy SoilsPolito, Carmine Paul (Virginia Tech, 1999-12-10)The presence of silt and clay particles has long been thought to affect the behavior of a sand under cyclic loading. Unfortunately, a review of studies published in the literature reveals that no clear conclusions can be drawn as to how altering fines content and plasticity actually affects the liquefaction resistance of a sand. In fact, the literature contains what appears to be contradictory evidence. There is a need to clarify the effects of fines content and plasticity on the liquefaction resistance of sandy soils, and to determine methods for accounting for these effects in engineering practice. In order to help answer these questions, a program of research in the form of a laboratory parametric study intended to clarify the effects which varying fines content and plasticity have upon the liquefaction resistance of sandy sands was undertaken. The program of research consisted of a large number of cyclic triaxial tests performed on two sands with varying quantities of plastic and non-plastic fines. The program of research also examined the applicability of plasticity based liquefaction criteria and the effects of fines content and plasticity on pore pressure generation. Lastly, a review of how the findings of this study may affect the manner in which simplified analyses are performed in engineering practice was made. The results of the study performed are used to clarify the effects of non-plastic fines content and resolve the majority of the inconsistencies in the literature. The effects of plastic fines content and fines plasticity are shown to be different than has been previously reported. The validity of plasticity based liquefaction criteria is established, the mechanism responsible for their validity is explained, and a new simplified criteria proposed. The effects of fines content and plasticity on pore pressure generation are discussed, and several recommendations are made for implementing the findings of this study into engineering practice.
- The Effects of Vibration on the Penetration Resistance and Pore Water Pressure in SandsBonita, John Anthony (Virginia Tech, 2000-07-28)The current approach for using cone penetration test data to estimate soil behavior during seismic loading involves the comparison of the seismic stresses imparted into a soil mass during an earthquake to the penetration resistance measured during an in-situ test. The approach involves an indirect empirical correlation of soil density and other soil related parameters to the behavior of the soil during the loading and does not involve a direct measurement of the dynamic behavior of the soil in-situ. The objective of this research was to develop an approach for evaluating the in-situ behavior of soil during dynamic loading directly through the use of a vibrating piezocone penetrometer. Cone penetration tests were performed in a large calibration chamber in saturated sand samples prepared at different densities and stress levels. A total of 118 tests were performed as part of the study. The piezocone penetrometer used in the investigation was subjected to a vibratory load during the penetration test. The vibratory units used in the investigations were mounted on top of a 1m section of drill rod that was attached at the lower end to the cone penetrometer. Pneumatic impact, rotary turbine, and counter rotating mass vibrators were used in the investigation. The vibration properties generated by the vibratory unit and imparted into the soil were measured during the penetration test by a series of load cells and accelerometers mounted below the vibrator and above the cone penetrometer, respectively. The tip resistance, sleeve friction and pore water pressure were also measured during the test by load cells and transducers in the cone itself. The vibration and cone data were compiled and compared to evaluate the effect of the vibration on the penetration resistance and pore water pressure in the soil mass. The results of the testing revealed that the influence of the vibration on the penetration resistance value decreased as the density and the mean effective stress in the soil increased, mainly because the pore water pressure was not significantly elevated throughout the entire zone of influence of the cone penetometer at the elevated stress and density conditions. An analysis of the soil response during the testing resulted in the generation of a family of curves that relates the soil response during the vibratory and static penetration to the vertical effective stress and density of the soil. The data used to generate the curves seem to agree with the proposed values estimated through the empirical relationship. An evaluation of the effects of the frequency of vibration was also performed as part of the study. The largest reduction in penetration resistance occurred when the input vibration approximated the natural frequency of the soil deposit, suggesting that resonance conditions existed between the input motion and the soil. An energy-based approach was developed to compare the energy imparted into the soil by the vibrator to the energy capacity of the soil. The input energy introduced into the soil mass prior to the reduction in penetration resistance agrees well with the energy capacity of the soil, especially in tests at the low effective stress level where a high excess pore water pressure was observed.
- Engineering Characteristics of Coal Combustion Residuals and a Reconstitution Technique for Triaxial SamplesLacour, Nicholas Alexander (Virginia Tech, 2012-06-19)Traditionally, coal combustion residuals (CCRs) were disposed of with little engineering consideration. Initially, common practice was to use a wet-scrubbing system to cut down on emissions of fly ash from the combustion facilities, where the ash materials were sluiced to the disposal facility and allowed to sediment out, forming deep deposits of meta-stable ash. As the life of the disposal facility progressed, new phases of the impoundment were constructed, often using the upstream method. One such facility experienced a massive slope stability failure on December 22, 2008 in Kingston, Tennessee, releasing millions of cubic yards of impounded ash material into the Watts Bar reservoir and damaging surrounding property. This failure led to the call for new federal regulations on CCR disposal areas and led coal burning facilities to seek out geotechnical consultants to review and help in the future design of their disposal facilities. CCRs are not a natural soil, nor a material that many geotechnical engineers deal with on a regular basis, so this thesis focuses on compiling engineering characteristics of CCRs determined by different researchers, while also reviewing current engineering practice when dealing with CCR disposal facilities. Since the majority of coal-burning facilities used the sluicing method to dispose of CCRs at one point, many times it is desirable to construct new "dry-disposal" phases above the retired ash impoundments; since in-situ sampling of CCRs is difficult and likely produces highly disturbed samples, a sample reconstitution technique is also presented for use in triaxial testing of surface impounded CCRs.
- EPOLLS: An Empirical Method for Prediciting Surface Displacements Due to Liquefaction-Induced Lateral Spreading in EarthquakesRauch, Alan F. (Virginia Tech, 1997-05-05)In historical, large-magnitude earthquakes, lateral spreading has been a very damaging type of ground failure. When a subsurface soil deposit liquefies, intact blocks of surficial soil can move downslope, or toward a vertical free face, even when the ground surface is nearly level. A lateral spread is defined as the mostly horizontal movement of gently sloping ground (less than 5% surface slope) due to elevated pore pressures or liquefaction in undelying, saturated soils. Here, lateral spreading is defined specifically to exclude liquefaction failures of steeper embankments and retaining walls, which can also produce lateral surface deformations. Lateral spreads commonly occur at waterfront sites underlain by saturated, recent sediments and are particularly threatening to buried utilities and transportation networks. While the occurrence of soil liquefaction and lateral spreading can be predicted at a given site, methods are needed to estimate the magnitude of the resulting deformations. In this research effort, an empirical model was developed for predicting horizontal and vertical surface displacements due to liquefaction-induced lateral spreading. The resulting model is called "EPOLLS" for Empirical Prediction Of Liquefaction-induced Lateral Spreading. Multiple linear regression analyses were used to develop model equations from a compiled database of historical lateral spreads. The complete EPOLLS model is comprised of four components: (1) Regional-EPOLLS for predicting horizontal displacements based on the seismic source and local severity of shaking, (2) Site-EPOLLS for improved predictions with the addition of data on the site topography, (3) Geotechnical-EPOLLS using additional data from soil borings at the site, and (4) Vertical-EPOLLS for predicting vertical displacements. The EPOLLS model is useful in phased liquefaction risk studies: starting with regional risk assessments and minimal site information, more precise predictions of displacements can be made with the addition of detailed site-specific data. In each component of the EPOLLS model, equations are given for predicting the average and standard deviation of displacements. Maximum displacements can be estimated using probabilities and the gamma distribution for horizontal displacements or the normal distribution for vertical displacements.
- An Examination of the Validity of Steady State Shear Strength Determination Using Isotropically Consolidated Undrained Triaxial TestsPorter, Jonathan R. (Virginia Tech, 1998-07-07)The assessment of the shear strength of soil deposits after the occurrence of large strains is an important issue for geotechnical engineers. One method for doing so, the steady state approach, is based on the assumption that the steady state undrained shear strength is a unique function of the in situ void ratio and effective stress. This method, which has been applied to liquefaction and flow failures, has been criticized because it may overestimate the in situ shear strength. The key to the steady state approach is accurate determination of the relationship between void ratio and effective stress at steady state. This is typically accomplished using conventional isotropically consolidated undrained (ICU) triaxial tests. The triaxial test was developed for measuring peak strengths, which typically occur at small strains, but steady state conditions typically occur at much larger strains. At large strain levels, the suitability of conventional triaxial testing procedures and error corrections is uncertain. The measured response at large strains may be inaccurate due to the influence of various testing errors. Furthermore, the true material response in the test specimen at large strains may not accurately represent in situ material behavior at large strains. This research effort consisted of an experimental and analytical study to examine the validity of steady state undrained shear strength determination using conventional ICU triaxial tests. The analytical study addressed triaxial testing errors and conventional corrections that are applied to test data and their influence on the measured steady state parameters. Finite element analyses were conducted to investigate the influence of variations in restraint at the end platens on stress distributions in the sample and measured stress-strain response. The finite element analyses incorporated axisymmetric interface elements to model the friction characteristics between the end platens and the specimen ends. The experimental study focused on several sands that are susceptible to liquefaction. An interface direct shear test program was conducted in order to evaluate various schemes for reducing end platen friction. ICU triaxial tests were conducted on each material using both conventional and lubricated end platens.
- An Experimental Study of the Dynamic Behavior of Slickensided SurfacesMeehan, Christopher Lee (Virginia Tech, 2006-01-25)When a clay soil is sheared, clay particles along the shear plane become aligned in the direction of shear, forming "slickensided" surfaces. Slickensided surfaces are often observed along the sliding plane in field landslides. Because the clay particles along a slickensided surface are already aligned in the direction of shear, the available shear resistance is significantly less than that of the surrounding soil. During an earthquake, ground shaking often causes landslide movement. For existing landslides or repaired landslides that contain slickensided rupture surfaces, it is reasonable to expect that the movement will occur along the existing slickensided surfaces, because they are weaker than the surrounding soil. The amount of movement that occurs is controlled by the dynamic resistance that can be mobilized along the slickensided surfaces. The objective of this study was to investigate, through laboratory strength tests and centrifuge model tests, the shearing resistance that can be mobilized on slickensided rupture surfaces in clay slopes during earthquakes. A method was developed for preparing slickensided rupture surfaces in the laboratory, and a series of ring shear tests, direct shear tests, and triaxial tests was conducted to study the static and cyclic shear resistance of slickensided surfaces. Two dynamic centrifuge tests were also performed to study the dynamic shear behavior of slickensided clay slopes. Newmark's method was used to back-calculate cyclic strengths from the centrifuge data. Test results show that the cyclic shear resistance that can be mobilized along slickensided surfaces is higher than the drained shear resistance that is applicable for static loading conditions. These results, coupled with a review of existing literature, provide justification for using cyclic strengths that are at least 20% larger than the drained residual shear strength for analyses of seismic stability of slickensided clay slopes. This represents a departure from the current state of practice, which is to use the drained residual shear strength as a "first-order approximation of the residual strength friction angle under undrained and rapid loading conditions" (Blake et al., 2002).
- Ground Improvement for Liquefaction Mitigation at Existing Highway BridgesCooke, Harry G. (Virginia Tech, 2000-07-13)The feasibility of using ground improvement at existing highway bridges to mitigate the risk of earthquake-induced liquefaction damage has been studied. The factors and phenomena governing the performance of the improved ground were identified and clarified. Potential analytical methods for predicting the treated ground performance were investigated and tested. Key factors affecting improved ground performance are the type, size, and location of the treated ground. The improved ground behavior is influenced by excess pore water pressure migration, ground motion amplification, inertial force phasing, dynamic component of liquefied soil pressure, presence of a supported structure, and lateral spreading forces. Simplified, uncoupled analytical methods were unable to predict the final performance of an improved ground zone and supported structure, but provided useful insights. Pseudostatic stability and deformation analyses can not successfully predict the final performance because of their inability to adequately account for the transient response. Equivalent-linear dynamic response analyses indicate that significant shear strains, pore water pressures and accelerations will develop in the improved ground when the treated-untreated soil system approaches resonance during shaking. Transient seepage analyses indicate that evaluating pore pressure migration into a three-dimensional improved zone using two-dimensional analyses can underestimate the pore pressures in the zone. More comprehensive, partially-coupled analyses performed using the finite difference computer program FLAC provided better predictions of treated ground performance. These two-dimensional, dynamic analyses based on effective stresses incorporated pore pressure generation, non-linear stress-strain behavior, strength reduction, and groundwater flow. Permanent movements of structures and improved soil zones were predicted within a factor of approximately two. Predictions of ground accelerations and pore water pressures were less accurate. Dynamic analyses were performed with FLAC for an example bridge pier and stub abutment on an approach embankment supported on shallow foundations and underlain by thick, liquefiable soils with and without improved ground zones. Ground improvement that restricted movements of the pier and stub abutment to tolerable levels included improved zones of limited size extending completely through the underlying liquefiable soils and formed through densification by compaction grouting or cementation by chemical grouting or jet grouting. A buttress fill at the abutment was unsuccessful.
- Liquefaction Susceptibility of Uncemented Calcareous Sands From Puerto Rico by Cyclic Triaxial TestingLaVielle, Todd Hunter (Virginia Tech, 2008-09-09)Laboratory tests were performed to investigate the liquefaction susceptibility of uncemented calcareous sands. A series of isotropically consolidated undrained monotonic and cyclic triaxial tests were performed using the Playa Santa sand from Porto Rico. Playa Santa sand is a poorly graded calcareous clean beach sand composed of angular particles with large intra-granular voids. A series of consolidated undrained triaxial tests were performed with the Playa Santa sand remolded to a variety of relative densities and consolidated under a range of confining pressures. In addition, cyclic triaxial tests were performed at a confining pressure of 100 kPa and three sets of relative densities (20%, 40% and 60%). Generation of excess pore pressure under different levels of cyclic loading was established. As a result, relationships were developed to relate the number of cycles required for triggering of liquefaction to cyclic stress ratio. It was seen that the Playa Santa sand was less susceptible liquefaction than quartzitic sands of the same relative density remolded and tested under similar conditions.
- A Numerical Investigation of the Seismic Response of the Aggregate Pier Foundation SystemGirsang, Christian Hariady (Virginia Tech, 2001-12-20)The response of an aggregate pier foundation system during seismic loading was investigated. The factors and phenomena governing the performance of the aggregate pier and the improved ground were identified and clarified. The key factors affecting the performance of the aggregate pier include soil density, stiffness modulus, and drainage capacity. The improved ground is influenced by soil stratification, soil properties, pore pressure dissipation, and earthquake time history. Comprehensive numerical modeling using FLAC were performed. The focus of the study in this research was divided into three parts: the studies of the ground acceleration, the excess pore water pressure ratio and the shear stress in soil matrix generated during seismic loading. Two earthquake time histories scaled to different peak acceleration were used in the numerical modeling: the 1989 Loma Prieta earthquake (pga = 0.45g) and the 1988 Saguenay earthquake (pga = 0.05g). The main results of the simulation showed the following effects of aggregate pier on liquefiable soil deposits: 1) The aggregate pier amplifies the peak horizontal acceleration on the ground surface (amax), 2) The aggregate pier reduces the liquefaction potential up to depth where it is installed, 3) Pore pressures are generally lower for soils reinforced with aggregate pier than unreinforced soils except for very strong earthquake, 4) The maximum shear stresses in soil are much smaller for reinforced soils than unreinforced soils. The excess pore water pressure ratio and the shear stress in the soil matrix calculated by FLAC were generally lower than those predicted by available procedures.
- Numerical modelling of transport of pollutant through soilsAhmad, Faheem (Virginia Tech, 1991-11-09)Prediction of subsurface migration of contaminant through soils involves analyses of unsaturated and saturated flow of water and advective dispersive transport of contaminant species. A finite element model is developed here for such an analysis. It is based on the transient nonlinear Richard's equation for the unsaturated flow and the mass transport equation using advective dispersive transport phenomenon. The model makes it possible to make advance predictions of the spread of the contaminant with respect to time and space, into the ground water system. The hydraulic properties of the unsaturated soils and the dispersion characteristics need to be obtained for such an analysis. The unsaturated flow parameters are obtained from a functional relationship between capillary pressure head and moisture saturation, and can be determined from laboratory tests on simple column samples of soils. A general expression is assumed to account for the effect of velocity dependence of the hydrodynamic dispersion coefficient in the mass transport problem. A computer program POLUT2D is developed based on the above assumptions. Pre and post processors for the computer program POLUT2D are also developed for interactive input of data and graphics displays of results. The computer program is first evaluated by comparing the results of a problem given in the literature with the results obtained by POLUT2D. The factors affecting the contaminant movement and distribution such as dispersivities, hydraulic conductivities and the effect of cutoff walls in controlling the spread of contaminant plume are studied. Also in this regard, the movement and spread of a contaminant at a landfill site in New Castle County, Delaware, is studied by comparing the simulated pattern of plume with the observed pattern.
- Passive Site Remediation for Mitigation of Liquefaction RiskGallagher, Patricia M. (Virginia Tech, 2000-10-27)Passive site remediation is a new concept proposed for non-disruptive mitigation of liquefaction risk at developed sites susceptible to liquefaction. It is based on the concept of slow injection of stabilizing materials at the edge of a site and delivery of the stabilizer to the target location using the natural groundwater flow. The purpose of this research was to establish the feasibility of passive site remediation through identification of stabilizing materials, a study of how to design or adapt groundwater flow patterns to deliver the stabilizers to the right place at the right time, and an evaluation of potential time requirements and costs. Stabilizer candidates need to have long, controllable gel times and low viscosities so they can flow into a liquefiable formation slowly over a long period of time. Colloidal silica is a potential stabilizer for passive site remediation because at low concentrations it has a low viscosity and a wide range of controllable gel times of up to about 100 days. Loose Monterey No. 0/30 sand samples (Dr = 22%) treated with colloidal silica grout were tested under cyclic triaxial loading to investigate the influence of colloidal silica grout on the deformation properties. Distinctly different deformation properties were observed between grouted and ungrouted samples. Untreated samples developed very little axial strain after only a few cycles and prior to the onset of liquefaction. Once liquefaction was triggered, large strains occurred rapidly and the samples collapsed within a few additional cycles. In contrast, grouted sand samples experienced very little strain during cyclic loading. What strain accumulated did so uniformly throughout loading and the samples remained intact after cyclic loading. In general, samples stabilized with 20 weight percent colloidal silica experienced very little (less than two percent) strain during cyclic loading. Sands stabilized with 10 weight percent colloidal silica tolerated cyclic loading well, but experienced slightly more (up to eight percent) strain. Treatment with colloidal silica grout significantly increased the deformation resistance of loose sand to cyclic loading. Groundwater and solute transport modeling were done using the codes MODFLOW, MODPATH, and MT3DMS. A "numerical experiment" was done to determine the ranges of hydraulic conductivity and hydraulic gradient where passive site remediation might be feasible. For a treatment are of 200 feet by 200 feet, a stabilizer travel time of 100 days, and a single line of low-head (less than three feet) injection wells, it was found that passive site remediation could be feasible in formations with hydraulic conductivity values of 0.05 cm/s or more and hydraulic gradients of 0.005 and above. Extraction wells will increase the speed of delivery and help control the down gradient extent of stabilizer movement. The results of solute transport modeling indicate that dispersion will play a large role in determining the concentration of stabilizer that will be required to deliver an adequate concentration at the down gradient edge. Consequently, thorough characterization of the hydraulic conductivity throughout the formation will be necessary for successful design and implementation of passive site remediation. The cost of passive site remediation is expected to be competitive with other methods of chemical grouting, i.e. in the range of $60 to $180 per cubic meter of treated soil, depending on the concentration of colloidal silica used.