Search

A technique was developed for analyzing and developing mechanical properties for porous geomaterials subjected to the high pressures encountered in penetration and blast-type loadings. A finite element (FE) code was developed to verify laboratory test results or to predict unavailable laboratory test data for porous media loaded to multikilobar stress levels. This FE program eliminates a deficiency in the process of analyzing and developing mechanical properties for porous geomaterials by furnishing an advanced analysis tool to the engineer providing properties to material modelers or ground shock calculators. The FE code simulates quasi-static, axisymmetric, laboratory mechanical property tests, i.e., the laboratory tests are analyzed as boundary value problems. The code calculates strains, total and effective stresses, and pore fluid pressures for fully- and partially-saturated porous media. The time dependent flow of the pore fluid is also calculated. An elastic-plastic strain-hardening cap model calculates the time-independent skeletal responses of the porous solids. This enables the code to model nonlinear irreversible stress-strain behavior and shear-induced volume changes. Fluid and solid compressibilities were incorporated into the code, and partially-saturated materials were simulated with a "homogenized" compressible pore fluid. Solutions for several verification problems are given as proof that the program works correctly, and numerical simulations of limestone behavior under drained and undrained boundary conditions are also presented. Ph. D.

Remedial ground densification is commonly used to reduce the liquefaction susceptibility of loose, saturated sand deposits, wherein controlled liquefaction is typically induced as the first step in the densification process. Assuming that the extent of induced liquefaction is approximately equal to the extent of ground densification, the purpose of this research is to assess the feasibility of using earthquake liquefaction data in remedial ground densification design via energy-based concepts. The energy dissipated by frictional mechanisms during the relative movement of sand grains is hypothesized to be directly related to the ability of a soil to resist liquefaction (i.e., Capacity). This hypothesis is supported by energy-based pore pressure generation models, which functionally relate dissipated energy to residual excess pore pressures. Assuming a linearized hysteretic model, a "simplified" expression is derived for computing the energy dissipated in the soil during an earthquake (i.e., Demand). Using this expression, the cumulative energy dissipated per unit volume of soil and normalized by the initial mean effective confining stress (i.e., normalized energy demand: NED) is calculated for 126 earthquake case histories for which the occurrence or non-occurrence of liquefaction is known. By plotting the computed NED values as a function of their corresponding SPT penetration resistance, a correlation between the normalized energy capacity of the soil (NEC) and SPT penetration resistance is established by the boundary giving a reasonable separation of the liquefaction / no liquefaction data points. NEC is the cumulative energy dissipated per unit volume of soil up to initial liquefaction, normalized by the initial mean effective confining stress, and the NEC correlation with SPT penetration resistance is referred to as the Capacity curve. Because the motions induced during earthquake shaking and remedial ground densification significantly differ in amplitude, duration, and frequency content, the dependency of the derived Capacity curve on the nature of the loading needs to be established. Towards this end, the calibration parameters for energy-based pore pressure generation models are examined for their dependence on the amplitude of the applied loading. The premise being that if the relationship between dissipated energy and pore pressure generation is independent of the amplitude of loading, then the energy required to generate excess pore pressures equal to the initial effective confining stress should also be independent of the load amplitude. However, no conclusive statement could be made from results of this review. Next, first order numerical models are developed for computing the spatial distribution of the energy dissipated in the soil during treatment using the vibratory probe method, deep dynamic compaction, and explosive compaction. In conjunction with the earthquake-derived Capacity curves, the models are used to predict the spatial extent of induced liquefaction during soil treatment and compared with the predicted spatial extent of improvement using empirical expressions and guidelines. Although the proposed numerical models require further validation, the predicted extent of liquefaction and improvement are in very good agreement, thus giving credence to the feasibility of using the Capacity curve for remedial ground densification design. Although further work is required to develop energy-based remedial densification design procedures, the potential benefits of such procedures are as follows. By using the Capacity curve, the minimum dissipated energy required for successful treatment of the soil can be determined. Because there are physical limits on the magnitude of the energy that can be imparted by a given technique, such an approach may lead to improved feasibility assessments and initial designs of the densification programs. Ph. D.

Bridges without expansion joints are called "integral bridges." Eliminating joints from bridges crates concerns for the piles and the abutments of integral bridges because the abutments and the piles are subjected to temperature-induced cyclic lateral loads. As temperatures change daily and seasonally, the lengths of integral bridges increase and decrease, pushing the abutment against the approach fill and pulling it away. As a result the bridge superstructure, the abutment, the approach fill, the foundation piles and the foundation soil are all subjected to cyclic loading, and understanding their interactions is important for effective design and satisfactory performance of integral bridges. The ability of piles to accommodate lateral displacements is a significant factor in determining the maximum possible length of integral bridges. In order to build longer integral bridges, pile stresses should be kept low. This research project investigated the complex interactions that take place between the structural components of the integral bridge and the soil through experimental and analytical studies. A literature review was conducted to gain insight into the integral bridge/soil interactions, and to synthesize the information available about the cyclic loading damage to piles of integral bridges. The ability of the piles and the abutments to withstand cyclic loads was investigated by conducting large-scale cyclic load tests. Three pile types and three semi-integral abutments were tested in the laboratory. Experiments simulated 75 years of bridge life for each specimen by applying over 27,000 displacement cycles. Numerical analyses were conducted to investigate the interactions among the abutment, the approach fill, the foundation soil, and the piles. The original VDOT semi-integral abutment hinge experienced shear key failure as observed in two large-scale laboratory tests. The revised hinge detail did not exhibit any sign of damage. Both abutments tolerated 75-year worth of displacement cycles without any appreciable change in their behavior. Semi-integral abutments are recommended for longer integral bridges because they can reduce pile stresses. As the need to build longer integral bridges grows, the role of the semi-integral abutments is expected to become more important. The data from the experimental program indicates that steel H-piles are the best pile type for support of integral abutment bridges. Concrete piles are not recommended because under repeated lateral loads, tension cracks progressively worsen and significantly reduce vertical load carrying capacity of these piles. Pipe piles have high flexural stiffness, which results in an undesired condition for the shear stresses in the abutment. For this reason, stiff pipe piles are not recommended for support of integral bridges. Numerical analyses indicate that the interactions between the approach fill and the foundation soils create favorable conditions for stresses in piles supporting integral bridges. Because of these interactions, the foundation soil acts as if it were softer, resulting in reduction in pile stresses compared to a single pile in the same soil without the approach fill above it. Ph. D.

A soil-bentonite cutoff wall is a type of subsurface vertical barrier constructed by back-filling a trench with a mixture of soil, bentonite, and water. Although soil-bentonite cutoff walls are common, their mechanical behavior is not well understood. Current design procedures do not consider the final stress state of the consolidated soil-bentonite backfill or deformations in adjacent ground. The final stress state in the completed wall is important because it influences the hydraulic conductivity of the cutoff (Barrier 1995), the cutoff's susceptibility to hydraulic fracture, and the magnitude of deformations adjacent to the cutoff wall. Deformations adjacent to the cutoff wall can be significant and can cause damage to adjacent structures. The objectives of this research are to 1) add to the current body of knowledge of the properties of soil-bentonite mixtures, 2) evaluate constitutive models and select a model to represent soil-bentonite, 3) model a soil-bentonite cutoff wall using finite elements, and 4) investigate the influence of several factors on the deformations in adjacent ground. These objectives were met by first summarizing information from the literature on soil-bentonite properties and then performing a laboratory testing program on different soil-bentonite mixtures. Five constitutive models were evaluated to determine how well they match the data from the laboratory testing program. A model referred to as the RS model was chosen to best represent soil-bentonite, and provided a good match of the soil-bentonite behavior. The RS model, which is a special case of a more complicated existing model, is a non-associative Modified Cam Clay type model that has parameters to change the yield surface and plastic potential surface into ellipses of varying shapes. The RS model was implemented into the finite element program SAGE. A finite element model was developed using SAGE to simulate all stages of construction of a soil-bentonite cutoff wall including excavation of a trench under bentonite-water slurry, replacement of the bentonite-water slurry with soil-bentonite backfill, and consolidation of the soil-bentonite backfill. The model was calibrated with a well-documented case history, and predicted deformations in adjacent ground were close to measured deformations. Evaluation of the model indicates that there is good confidence in the prediction of deformations in adjacent ground, but there is lower confidence in the predicted stresses in the consolidated soil-bentonite and settlement of the backfill in the trench. A parametric study was then performed using the finite element model assuming sand sites of varying density and OCR. Deformations in adjacent ground were calculated for various soil conditions, soil-bentonite properties, and trench configurations. A correlation was found between maximum calculated settlement in adjacent ground and factor of safety against trench Ph. D.

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. Ph. D.

This dissertation describes enhancements made to the finite element program, SAGE, and research on the performance of deep excavations. SAGE was developed at Virginia Tech for analysis of soil-structure interaction problems (Morrison, 1995). The purpose of the work described in this text with SAGE was to increase the capabilities of the program for soil-structure analysis. The purpose of the research on deep excavations was to develop a deeper understanding of the behavior of excavation support systems. The significant changes made to SAGE during this study include implementation of Biot Consolidation, implementation of axisymmetric analysis, and creation of a steady state seepage module. These changes as well as several others are described. A new manual for the program is also included. A review of published studies of deep excavation performance and recent case histories is presented. Factors affecting the performance of excavation support systems are examined, and performance data from recent published case histories is compared to data from Goldberg et al.'s 1976 report to the Federal Highway Administration. The design, construction, and performance of the deep excavation for the Dam Number 2 Hydroelectric Project is described. Finite element analyses of the excavation that were performed with SAGE are presented and discussed. Ph. D.

This dissertation describes the development and verification of a general purpose computer program, CONSOL97, for analysis of one dimensional consolidation of multi-layered soil profiles. The program uses an elasto-visco-plastic model that can simulate both consolidation and creep in a single consistent analysis. The finite element program uses standard oedometer test data to model stress-strain-time relationships, and the effect of strain rate on preconsolidation stress observed in the laboratory and in the field. Validation of the computer program by simulating standard oedometer tests is described and the applicability of the program in predicting field behavior is examined. The oedometer test simulations indicate good agreement with stress-strain and strain-log time test results during loading. Unloading behavior produces excessive rebound. Well-instrumented field tests at Väsby, Sweden, Skå-Edeby, Sweden and Berthierville, Canada indicate that elasto-visco-plastic CONSOL97 analyses produce better predictions of field behavior than conventional elasto-plastic models. CONSOL97 results were in good agreement for the Väsby and Berthierville test fills but underestimated displacements and pore pressures near the center of the normally consolidated clay layer beneath the Skå-Edeby test fill. Ph. D.

A three-dimensional finite element model for the analysis of steady state seepage has been presented in this study. The theory of unsaturated flow has been used in the analysis of steady state seepage. The model applies the invariant mesh procedure in the finite element analysis. Galerkin's method is used in the formulation of the finite element equations. The pre and the post processor developed in the generation and viewing of the finite element mesh and the free surface has also been discussed in this study. The study presents the comparison of results obtained from the three-dimensional model with a previously validated two-dimensional model. Master of Science

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin has been verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTl 1M7/8HS preform and corn oil. Fluid/fiber interaction was measured through the use of the single fiber pull-out test. The surface tension of both the corn oil and Epon 815 was determined. The contact angle between these two fluids and glass and carbon fibers was also measured. These tests indicated that the glass fiber had a lower contact angle than the carbon fiber and therefore is wet out better than the carbon fiber by both fluids. This result is attributed to the sizing commonly used on the carbon fibers. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle. Master of Science

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. Ph. D.

Numerical methods and specifically the finite element method have improved significantly since their introduction in the 60's. These advances were mainly in: 1) introducing higher-order elements, 2) developing effective solution schemes, 3) developing sophisticated means of modeling the constitutive behavior of geotechnical materials, and 4) introducing iteration techniques to model material non-linearity. This thesis, on the other hand, deals with the topic of modeling the boundary conditions of the finite element problem. Typically, the boundary conditions will be approximated by specifying displacement constraints. such as restraining the bottom boundary of the finite element mesh against displacements in the horizontal and vertical directions (x- and y-directions). Where bedrock or dense residual soils underlie the soft foundation soil at a relatively shallow depth, this is a good assumption. However. when soft soil is encountered for large depths, the assumption of zero movement constraints for a mesh boundary at a shallower depth than the actual bedrock will result in a serious underestimation of stresses and displacements. By coupling boundary elements to the finite elements and using them to model the infinite extent of the foundation soil, a more realistic answer is obtained. Employing the coupled boundary element - finite element method, four cases were analyzed and the results compared to values of the pure finite element method. The results show that the coupled method indeed yielded higher stress- and displacement-values, indicating that the pure finite element method underestimates stresses and displacements when modeling very deep soils. Master of Science

Experimental and analytic investigations were performed to examine the influences of wall height, backfill behavior, and compaction on the magnitudes of backfill loads on rigid retaining walls. Measurements of lateral and vertical backfill loads were made during tests using the Virginia Tech instrumented retaining wall facility. The tests were performed with two soils, moist Yatesville silty sand and dry Light Castle sand. Two hand-operated compactors, a vibrating plate compactor and a rammer compactor, were used to compact the backfill. The backfill height was 6.5 feet in all of the tests. Analyses of backfill loads were made using a compaction- induced lateral earth pressure theory and a vertical shear force theory. The compaction-induced lateral earth pressure theory was revised from a previous theory. The revisions improved the accuracy with which the theory models the hysteretic stress behavior of the backfill during compaction. The theory was also extended to include the pore pressure response of moist backfill in a rational manner. A vertical shear force theory was also developed during this research. The theory is based on consideration of backfill compressibility and mobilization of interface shear strength at the contact between the backfill and the wall. The theory provides a useful basis for understanding how wall height, backfill compressibility, wall-backfill interface behavior, and compaction-induced lateral pressures affect the vertical shear forces on rigid walls. Studies were also made to investigate the cause of erratic pressure cell readings. An important cause of drift in pressure cell readings was found to be moisture changes in the concrete in which the pressure cells were mounted. It was found that this problem could be mitigated by applying a water-seal treatment to the face of the wall. Both the vibrating plate compactor and the rammer compactor were instrumented to measure dynamic forces and energy transfer during compaction. The force applied by the vibrating plate compactor was about one-quarter of the manufacturer’s rated force. The force applied by the rammer compactor was about twice the manufacturer’s rated force. The transferred energy measurements provided a basis for relating laboratory and field compaction procedures. Ph. D.

The purpose of this study was to develop a theory for clay deformation that models the behavior represented by Bjerrum's time lines (Bjerrum (1961). The theory models time-dependent compression (secondary compression) and the "aging" effect, or increase in preconsolidation pressure with time under sustained load. The theory appears to be consistent with many aspects of the consolidation behavior of real clays, including: (1) secondary compression, (2) aging, (3) different values of Cv determined by Casagrande's method and Taylor's method, and (4) faster settlement in the field than expected on the basis of conventional theory. It thus appears to have considerable potential for use in interpreting laboratory tests and predicting field settlements more accurately than has been possible up to the present time. In addition to development of the new theory, the study involved an extensive review of previous work and suggestions for possible future studies." 14. The results of this study lead to the conclusion that the OWL model did prepare the Shawsville children for first grade experiences and that their achievement would be at least as good as those children who participated in a traditional kindergarten program. Ph. D.

The results of this work are the closed-form expressions for the tangent and secant stiffness matrices of the 3-node Mindlin frame element. From the study of two solution methods and the behavior of the three elements, several conclusions have been made as follows: 1. Because the Mindlin frame element uses quadratic polynoms for the interpolation functions, at least two elements are needed for modelling a member of a plane frame structure. 2. The Mindlin frame element will give a more flexible structure due to the shear deformation effect, but the effect is not always significant. 3. The Mindlin frame element requires about twice degree-of-freedoms than the Bernoulli-Euler frame element requirement. 4. The formulation of tangent stiffness matrix can be simplified, in condition that the internal forces calculations are done using the appropriate nonlinear formulation. The calculation of the tangent stiffness matrix and the internal element forces can be done using different formulations. Master of Science

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. Master of Science

A deep-bed drying model simulating the drying of peanuts in a fixed bed is required for designing energy-efficient and automatically controlled dryers. A deep-bed drying model consists of a thin-layer drying model to calculate the moisture release from the material and a set of mass and energy balances. An experimental setup was constructed to determine drying rates of Virginia-type peanuts under 14 different drying air conditions. Selected empirical and semi-theoretical models available for modeling thin-layer drying rates were fitted to the collected data using nonlinear regression techniques. The modified Page's model and the two-term exponential model fitted the data better than other models considered. A deep-bed drying model PEATECH based on four coupled partial differential equations consisting of four variables, air temperature, peanut temperature, air humidity, and peanut moisture content was developed. Validation of the model was accomplished by using the data collected from 36 deep-bed drying experiments conducted using three laboratory dryers during 1987, 1988, and 1989. PEATECH predicted the variables within a peanut bed with an accuracy of less than ± 6%. The energy saving potential of exhaust-air recirculation was established by conducting simulated experiments using a modified version of PEATECH. Ph. D.

Cohesionless sands are known to be susceptible to failure by liquefaction when they are saturated and subjected to earthquake shaking. Considerable study has been directed towards this subject over the past 20 years in recognition of the possibility of large-scale property damage or loss of life due to this type of failure. Recent evidence has shown that small degrees of cementation in a sand significantly reduce the likelihood of liquefication. However, the work to date has been limited to studies with conventional testing devices and simple loading paths. These devices are suspected of inducing premature failure in cemented soils, and are not capable of modeling the effects of multiaxial loading. In this investigation, there were two major objectives. The first involved the development and fabrication of a new three~dimensional shear device with the capability of applying load to cemented sands with a minimum of stress concentration effects, and of using load paths which are more representative of the true effects of an earthquake than is possible in conventional equipment. The second concerned performance of a series of production tests to investigate the behavior of cemented sands under a range of earthquake loading paths. The production tests were largely performed using the new three—dimensional shear device. The test results showed that cemented soils have more resistance to earthquake loading than previously thought since stress concentration effects in conventional testing do induce premature failure through the effects of stress concentrations. On the other hand, it was found that either cemented or uncemented sands show less resistance to earthquake loadings if multiaxial stress conditions are applied to the sample as opposed to uniaxial loadings. This is important in explaining the fact that sites with seemingly similar conditions often show different behavior, since slightly different earthquake loading pattems can cause different responses. One factor explaining differences in response is found to be the mean normal stress, which is not the same for all loadings, and plays an important role in the pore pressures developed in the soil. Ph. D. incomplete_metadata

The finite element method is used in determining the load-capacity and behavior of a single concrete pile modeled after full scale instrumented test piles which were employed in the Tidewater area of Virginia. Specifically, this report describes the subsurface conditions and full scale pile tests that were performed during the design phase of the foundation for the Land Level Shipbuilding Facility in Newport News, Virginia. It also addresses the testing conducted by the author to determine the material properties required for finite element analyses, including: soil sampling and in-situ tests within the Yorktown formation; standard soil laboratory tests for index properties and strength-deformation characteristics of the soil along the pile perimeter and tip, plus interface tests to evaluate soil strength-deformation characteristics along the soil-pile interface. The data derived from the finite element method of analysis is then compared with the load-displacement and load-transfer data obtained from the full scale instrumented pile load tests. The comparative analysis shows favorable results leading to the conclusion that the finite element method shows promise for use in the design of deep foundations systems in the Tidewater area of Virginia. Master of Science

First-hand reports on the 1886 Charleston earthquake contain numerous accounts for the widespread occurrence of liquefaction related features in and near the meizoseismal zone. Recent geologic studies have found evidence for the repeated liquefaction of sandy soils in the Charleston area due to recurring large seismic events. In the course of this investigation 24 mini-cone penetration tests were performed at seven sites in or near the meizoseismal zone of the 1886 earthquake to determine the factors influencing ground failures due to liquefaction. These tests were supplemented with soil borings, sieve tests and a limited number of standard penetration tests to aid in characterization of the sandy soils. Additionally, soil boring records in Charleston were obtained which provided in-situ testing data in an area with documented historical damage. The range of sites at which testing was done, or information was available, represent locations experiencing various levels of liquefaction and distances from the zone of seismic energy release. Penetration data were used to evaluate resistance of the sandy soil to cyclic loading and fonned the basis for assessing the effects that the lateral extent and distribution of loose sand layers has on the surficial manifestation of liquefaction. With the absence of cementation and extensive soil development, soils as old as late Pleistocene age have been found to be very susceptible to liquefaction. At several sites these soils have undergone at least three episodes of liquefaction and presently exhibit low penetration resistances, indicating that the progressive densification of a liquefiable soil layer can be minor unless it is in very close proximity to a large venting feature. The size and density of occurrence of vents and sand blows has been found to be primarily dependent on the extent of both the liquefiable layer and any overlying resistant layers. Layered system relations utilized with field performance data, and historical and geologic evidence for the occurrence of liquefaction features to suggest that the near surface peak horizontal accelerations induced by the 1886 earthquake were approximately O.3g in the meizoseismal zone and O.2g in the city of Charleston. This is in contrast with previous estimates of seismic shaking all of which point toward values in the range of 0.5 to O.6g. The reason for the different acceleration estimates is not clear at this time, and will be further studied in future extension of this work. Master of Science