Liquefaction Centrifuge Modeling of Sands of Different Permeability

This paper presents the results of six centrifuge model tests of liquefaction and earthquake-induced lateral spreading of fine Nevada sand using an inclined laminar box. The centrifuge experiments simulate a gently sloping, 10 m thick stratum of saturated homogeneous sand of infinite lateral extent and relative densities ranging from 45 to 75%. Such idealized models approach some field situations and they provide significant general insight into the basic mechanisms and parameters influencing the lateral spreading phenomenon. The layer was subjected to lateral base shaking with prototype peak acceleration ranging from 0.20 to 0.41 g, a frequency of 2 Hz, and duration of approximately 22 cycles. The simulated field slope angle was 5°. The model deposits were all saturated with a viscous fluid 50 times more viscous than water, so that testing under the increased gravitational field (50 g) produced a deposit with the prototype permeability of the same fine-grained sand saturated with water in the field. Detailed discussions and comparisons of the six centrifuge tests are included. The observed effects of relative density Dr and input peak acceleration amax on the following measured parameters are summarized: thickness of liquefied soil H1, permanent lateral displacement DH, and ground surface settlement S. Comparisons and discussions are also presented on the effect of permeability for a Dr = 45% deposit. This is done by comparing the results reported herein using a viscous pore fluid, with other published centrifuge tests where a similar deposit using the same model soil, also tested at 50 g and shaken with the same input motion, was saturated with water, thus simulating a prototype sand having 50 times the permeability of the fine sand reported in this paper.     

Resonant Column Testing of Sands with Different Viscosity Pore Fluids

The effect of pore fluid viscosity on the stiffness, damping, and liquefaction characteristics of sands was investigated to assess the potential contributions to a centrifuge model seismic response for soils saturated with high-viscosity pore fluid. Resonant column tests with cyclic loading frequencies in the range of 20–45 Hz were performed on a variety of fluids and sand sizes. At a strain level less than 2 × 10?4, the damping increased with pore fluid viscosity and shear strain amplitude, and it decreased with sand particle size. However, at shear strain greater than about 2 × 10?4, the increased skeleton damping tended to mask any effect of additional damping due to fluid viscosity. The liquefaction tests on fine silica sand revealed that the increase in total energy dissipation was not more than 10% for 100 cS oil when compared with water at a driving frequency of 25 Hz. Based on the experimental results, a simple model is proposed to examine the dependency of viscous damping on pore fluid viscosity, loading frequency, particle size, and shear strain amplitude.     

Preferential Flow of a Nonaqueous Phase Liquid in Dry Sand

Geotechnical centrifuge testing is used to examine the preferential (fingered) flow of a nonaqueous phase liquid (NAPL) in a uniform dry sand. The results of nine experiments, containing a total of 87 observations of NAPL finger behavior, are analyzed. The observed finger tip velocities range from 0.01 to 0.3 cm/s, while the observed finger widths range from 0.3 to 3.6 cm. From the experimental data it is concluded that, asymptotically, the NAPL fingers are not fully saturated. For comparison, the behavior of water fingers is examined using the same experimental setup. In contrast to the NAPL fingers, and in agreement with other work reported in the literature, the water fingers are found to be fully saturated. In addition, it is confirmed that the water finger properties can be well predicted from known porous medium and fluid properties. A scaling analysis is presented that allows the NAPL finger properties to be inferred from models developed to describe water finger properties. The analysis predicts NAPL finger velocities to within 15% and NAPL finger widths to within 50% if both finger types are assumed saturated. By adjusting the analysis to account for the fact that the NAPL fingers are not fully saturated, NAPL finger widths can be predicted within to 10%, and NAPL finger velocities to within 30%.     

Behavior of Concrete in Water Subjected to Dynamic Triaxial Compression

To understand the behavior of concrete material in ambient water, a series of triaxial compressive tests of concrete cylindrical specimens (? 100×200?mm) was conducted on a large scale triaxial machine. The acting pattern of water, confining pressure, loading strain rate, and moisture content were chosen as test parameters. The water acting patterns on concrete were directly divided into mechanical loading and real water loading according to whether the specimens were directly exposed to water or not. The confining pressure ranged from 0–8 MPa and the strain rate included 10?5/s, 10?3/s, and 10?2/s. By testing dry and saturated specimens, the effect of moisture on concrete strength was also examined. The test results indicated that the compressive strengths of both dry and saturated concrete increase obviously with the confining pressure under mechanical confining pressure. However, the effect on the strengthened dry concrete specimens is more significant. The strength of dry concrete under real water loading decreased remarkably, even less than its uniaxial strength, whereas the compressive strength of the saturated concrete specimen under real water loading is close to its uniaxial compressive strength. The strength of concrete increases with strain rate, and this phenomenon becomes more apparent under water loading.     

Viscous Heating of Fluid Dampers under Small and Large Amplitude Motions: Experimental Studies and Parametric Modeling

This paper summarizes the results from a comprehensive experimental program in an effort to better understand the phenomenon of viscous heating of fluid dampers under small-stroke (wind loading) and large-stroke (earthquake loading) motions. Two dampers, one with 15-kip and one with 250-kip force output at peak design velocity, have been instrumented and tested under various amplitudes and frequencies. Temperature histories at different locations along the damper casing and within the silicon fluid that undergoes the shearing action have been recorded. Experimental data under small-stroke motions of the 250-kip damper showed that a single closed-form expression derived from first principles is capable of predicting the temperature rise at different locations of the damper with fidelity. The recorded data under long-stroke motions suggest a two-parameter law of cooling that allows the estimation of the internal temperature of the silicon oil once the external temperature on the damper casing is known. The presented cooling law is an extension of Newton’s law of cooling. The study concludes that for both dampers, the same values of the model parameters provide a good approximation of the cooling behavior. The study presents a valuable formula that can be used in practice to estimate the internal fluid temperature of the damper given the external shell temperature.     

Use of Exact Solutions of Wave Propagation Problems to Guide Implementation of Nonlinear Seismic Ground Response Analysis Procedures

One-dimensional nonlinear ground response analyses provide a more accurate characterization of the true nonlinear soil behavior than equivalent-linear procedures, but the application of nonlinear codes in practice has been limited, which results in part from poorly documented and unclear parameter selection and code usage protocols. In this article, exact (linear frequency-domain) solutions for body wave propagation through an elastic medium are used to establish guidelines for two issues that have long been a source of confusion for users of nonlinear codes. The first issue concerns the specification of input motion as “outcropping” (i.e., equivalent free-surface motions) versus “within” (i.e., motions occurring at depth within a site profile). When the input motion is recorded at the ground surface (e.g., at a rock site), the full outcropping (rock) motion should be used along with an elastic base having a stiffness appropriate for the underlying rock. The second issue concerns the specification of viscous damping (used in most nonlinear codes) or small-strain hysteretic damping (used by one code considered herein), either of which is needed for a stable solution at small strains. For a viscous damping formulation, critical issues include the target value of the viscous damping ratio and the frequencies for which the viscous damping produced by the model matches the target. For codes that allow the use of “full” Rayleigh damping (which has two target frequencies), the target damping ratio should be the small-strain material damping, and the target frequencies should be established through a process by which linear time domain and frequency domain solutions are matched. As a first approximation, the first-mode site frequency and five times that frequency can be used. For codes with different damping models, alternative recommendations are developed.     

Artificial Ground Freezing of Fully Saturated Soil: Viscoelastic Behavior

The transport and mechanical properties of saturated soil drastically change when temperatures drop below the freezing temperature of water. During artificial ground freezing, this change of properties is exploited in order to minimize deformations during construction work and for groundwater control. Whereas for the latter only the size of the frozen-soil body is relevant, which is obtained by solving the thermal problem, the design of the ground-freezing work for support purposes requires information about the mechanical behavior of frozen soil. In addition to the quantification of the improvement of mechanical properties during freezing, information about the dilation associated with the 9% volume increase of water during freezing is required in order to assess the risk of damage to surface infrastructure caused by frost heave. In this paper, a micromechanics-based model for the prediction of both the aforementioned phase-change dilation and the elastic and viscous properties of freezing saturated soil is presented. Hereby, the macroscopic material behavior is related to the behavior of the different constituents such as soil particles, water, and ice. Combined with the solution of the thermal problem, the proposed model provides the basis for predictions of the performance of support structures composed of frozen soil.     

Hydraulic Damping of Saturated Poroelastic Soils During Steady-State Vibration

A theoretical study of the steady-state response of a saturated poroelastic soil column during compressional and rotational harmonic vibrations is presented. Hydraulic damping due to Biot flow is evaluated for top-drained and double-drained boundary conditions and for compressional and rotational motions using the theory of a damped single-degree-of-freedom system. For compressional motions, the dynamic response of gravels and sands is highly influenced by the compressibility of the pore fluid. More hydraulic damping occurs as soil hydraulic conductivity increases and as the column boundary conditions change from top drained to double drained. On the other hand, hydraulic damping for rotational motions is significantly less than that for compressional motions and is dependent on a dimensionless hydraulic conductivity parameter Ks. For Ks within the range of 10?3–100, hydraulic damping may have an important contribution to total soil damping, especially at small strain levels.     

Dynamics of SDOF Systems with Nonlinear Viscous Damping

The effects of nonlinear viscous damping on the dynamic response of single-degree-of-freedom (SDOF) structural systems are analyzed. This kind of damping characterizes a special class of fluid viscous dampers recently utilized in the field of vibration control as base-isolation devices or viscoelastic elements included in steel braces of framed structures. The analytical relationship adopted to reproduce the mechanical behavior of the fluid viscous dampers is a fractional power-law of the velocity, the exponent of which ranges between 0.1 and 0.2. This function had been previously calibrated on the results of a special experimental survey carried out at the University of Florence. The dynamics of the classical linear-viscous SDOF oscillator is herein reformulated on the basis of the above-mentioned fractional viscous damping (FrVD) relationship. In particular, the transient and steady-state responses are examined in both free and forced vibration conditions. The magnification and transmissibility factors are analytically determined for different damping levels. Moreover, the relation between the viscous damping coefficient and the frequency ratio (i.e., the ratio of the dynamic load to the oscillator frequencies) is defined. The diagrams describing these functions provide direct correlations between the damping as well as the elastic properties of the system and the frequency content of the dynamic action.     

Mixture Theory for a Fluid-Saturated Isotropic Elastic Matrix

The theory for a fluid saturated linearly isotropic elastic matrix is still the basis for many geophysical applications, and commonly adopts Biot’s symmetric stress–strain laws for the matrix stress and fluid pressure. These involve a shear modulus and three elastic moduli governing the mixture and constituent compressions, in contrast to four compression moduli if Biot’s invalid potential energy argument is not applied. We now show that an energy argument applied to undrained loading also leads to three compression moduli, but distinct from those derived by Biot (Biot symmetry). However, there are two distinct solutions of this energy balance, corresponding to the Voigt and Reuss limits of the analogous theory of a linear two-phase elastic composite, whereas a unique undrained modulus not at either limit would be expected. It is proposed that an energy contribution is lost due to the idealised assumptions made for undrained loading, which therefore does not determine a further restriction, so that there are four independent compression moduli. The general and restricted combinations of the total pressure and fluid pressure (effective stress) governing the matrix compression are then presented, together with the alternative forms of the partial differential equations governing the deformation and flow.     

Summary of SAC Case Study Building Analyses

Summarized in this paper are the major findings from analytical studies of nine steel moment frame buildings conducted under Phase 1 of the SAC Steel Project. The buildings range in height from two to seventeen stories and most of them experienced damage to welded beam-column connections during the Northridge earthquake of 1994. Elastic response spectrum, inelastic static pushover, and elastic and inelastic time-history analyses were conducted using ground motion data representative of the Northridge earthquake to establish the loading∕deformation demands that the buildings experienced. The primary performance indices obtained from the analyses were demand-to-capacity ratios, interstory drift ratios, and inelastic hinge rotations. Maximum ratios of elastic member force demands to plastic strengths ranged between 1.0 and 2.0; maximum inelastic hinge rotations were 0.005–0.010 rad; and maximum interstory drift ratios were from 1 to 2%. These damage indices increased by 50%–150% under more severe ground motions recorded during the Northridge earthquake at the Sylmar site. Accuracy of the analyses is shown to be sensitive to a number of modeling parameters including finite joint size, joint panel behavior, composite beam action, strain hardening, second-order (P-Δ) effects, and three-dimensional response. Overall, there was only modest correlation between the frame performance indices and the observed connection damage, due largely to the fact that significant aspects of the connection fracture behavior are not captured in the frame analyses.     

Experimental Study of Particle Motion on a Smooth Bed under Shoaling Waves Using Particle Image Velocimetry

The motion of spherical particles (diameter 1.58 mm, specific gravity 2.5) on 2 and 3% plane slope was studied in a laboratory wave flume for shoaling wave conditions. The range of wave-height-to-water-depth ratio was 0.24相似文献   

Generalized Model for Geosynthetic-Reinforced Granular Fill-Soft Soil with Stone Columns

This paper pertains to the development of a mechanical model to predict the behavior of a geosynthetic-reinforced granular fill over soft soil improved with stone columns. The saturated soft soil has been idealized by Kelvin–Voight model to represent its consolidation behavior. The stone columns are idealized by stiffer springs. Pasternak shear layer and rough elastic membrane represent the granular fill and geosynthetic reinforcement layer, respectively. The nonlinear behavior of the granular fill and the soft soil is considered. Effect of consolidation of the soft soil due to inclusion of the stone columns has also been included in the model. Plane strain conditions are considered for the loading and reinforced foundation soil system. An iterative finite difference scheme is applied for obtaining the solution, and results are presented in nondimensional form. Comparison between the results from the present study and the analytical solution using theory of elasticity shows reasonable agreement. The advantage of using geosynthetic reinforcement is highlighted. Results indicate that inclusion of the geosynthetic layer effectively reduces the settlement. Nonlinearity in the behavior of the soft soil and the granular fill is reduced due to the use of geosynthetic reinforcement layer.     

SWAN-Mud: Engineering Model for Mud-Induced Wave Damping

This paper describes the implementation of a new dispersion relation and energy-dissipation equation obtained from a viscous two-layer model schematization in the state-of-the-art wave forecasting model SWAN to simulate wave damping in coastal areas by fluid mud deposits. This new dispersion relation is derived for a nonviscous, nonhydrostatic upper layer and a viscous, hydrostatic lower layer, covering most conditions encountered in nature. An algorithm is developed for a robust numerical solution of this new implicit dispersion relation through proper starting values in the iteration procedure. The implementation is tested against a series of analytical solutions and three schematic test cases. Next, four dispersion relations published in the literature are evaluated and compared with the new dispersion relation. The solution of the dispersion relations forms a multidimensional space. Comparison of the various model solutions through one-dimensional graphs can therefore become quite misleading, as shown in the discussion of a two-dimensional representation of the solution space, explaining for instance the variation in ambient conditions at which maximum wave damping is to be expected. The various models have been developed for a variety of conditions, such as shallow and deep water and shallow and thick mud layers; the various models agree well in their domain of applicability, but they show significant deviations when used outside their domain. Because the ambient and mud conditions may vary considerable in space and time at a particular site, the use of the new model is advocated because it covers most water depths and fluid mud thicknesses encountered in nature. The strength of the new SWAN-mud model lies in its large-scale applicability, assessing the two-dimensional evolution of wave fields in coastal areas. Therefore, the new implementation is evaluated with respect to the behavior of waves on a sloping seabed, representing real-world coasts. In all cases, the new SWAN-mud model behaves satisfactorily; a critical remaining issue, though, is the assessment of the relevant fluid mud parameters.     

Optimal Control of Earthquake Response Using Semiactive Isolation

The responses of two, low-rise, 2-degree-of-freedom base isolated structures with different isolation periods to a set of near-field earthquake ground motions are investigated under passive linear and nonlinear viscous damping, two pseudoskyhook semiactive control methods, and optimal semiactive control. The structures are isolated with a low damping elastic isolation system in parallel with a controllable damper. The optimal semiactive control strategy minimizes an integral norm of superstructure absolute accelerations subject to the constraint that the nonlinear equations of motion are satisfied and is determined through a numerical solution to the Euler–Lagrange equations. The optimal closed-loop performance is evaluated for a controllable damper and is compared to passive viscous damping and causal pseudoskyhook control rules. Results obtained from eight different earthquake records illustrate the type of ground motions and structures for which semiactive damping is most promising.     

Is P-Wave Velocity an Indicator of Saturation in Sand with Viscous Pore Fluid?

It is commonly assumed that within inundated sand the Skempton B value and P-wave velocity decrease with decrease in saturation. In centrifuge tests a common saturation procedure is to inundate the specimen with carbon dioxide while under a vacuum and then slowly introduce the viscous pore fluid. The B value and related saturation is difficult to measure in centrifuge models and P-wave velocity—saturation correlations have been used for this purpose. A laboratory emulation of centrifuge saturation procedures was made using a triaxial cell with top and bottom bender elements and a viscous methyl cellulose–water pore fluid. Contrary to expectations, the laboratory tests showed high P-wave velocities indicative of full saturation when B values were low. Numerical modeling of the laboratory tests indicated that if air bubbles within the pore fluid are numerous and closely spaced then there is a good correlation between saturation, B value, and P-wave velocity. However if the air bubbles are larger and only present in some of the pores then the P-wave velocity is not a good indicator of B value and average saturation. The laboratory tests also showed that placing the specimen under backpressure for several days increased saturation and related B values. It is suggested that this common laboratory procedure should be considered for saturating centrifuge test specimens.     

Combined Continua and Lumped Parameter Modeling for Nonlinear Response of Structural Frames to Impulsive Ground Shock

The response of a beam-column frame to impulsive ground shock, such as those induced by an underground explosion, has characteristics of both impact and natural earthquake responses. The critical effects may be governed by the dynamic response of individual elements as continuous mass systems, while to a certain extent the global vibration (as of lumped-mass systems) may also be involved. To incorporate both dynamic features, the present study proposes a combined continua and lumped parameter (CCLP) model, which consists of the basic beam-column element with distributed stiffness and mass, along with concentrated mass-springs at element ends to form the reduced dynamic system. To take into account of the shear deformation and rotary inertia which become important in the impulsive response, the governing equations are formulated based on the Timoshenko beam theory. The nonlinearities are described through three mechanisms, namely the distributed nonlinear flexural and diagonal shear behavior, and the direct sliding shear at the member ends. A generic restoring force model is adopted to describe the hysteretic behavior. Comparison with a scaled model test demonstrates that the CCLP model is capable of representing the primary dynamic features in a frame structure under impulsive ground shock. Extended parametric studies indicate that, with increase of the ground shock frequency, the failure tends to become shear dominant. For ground shocks of frequency at 20–30?Hz and above, the failure in a reinforced concrete column will require a peak ground velocity (PGV) on the order of 3?m/s, whereas failure in a beam would occur at PGV of about 1.5?m/s.     

Simulation of Upward Seepage Flow in a Single Column of Spheres Using Discrete-Element Method with Fluid-Particle Interaction

The interaction between soil particles and pore water makes the behavior of saturated and unsaturated soil complex. In this note, upward seepage flow through a granular material was idealized using a one-column particle model. The motion of the individual particles was numerically simulated using the discrete-element method taking the interaction with the fluid into account. The fluid behavior was simulated by the Navier-Stokes equation using the semi-implicit method for pressure-linked equation. This approach has already been applied in powder engineering applications. However, there are very few studies that have used this approach in geotechnical engineering. This note first describes the qualitative and quantitative validation of this method for hydraulic gradients below the critical one by comparing the results with an analytical solution. Then, the ability of the method to simulate the macroscale behavior due to the interaction between particles and pore water at hydraulic gradients exceeding the critical hydraulic gradient is discussed.     

Damping of Cables by a Transverse Force

A general asymptotic format is presented for the effect on the modal vibrations of a transverse damper close to the end of a cable. Complete locking of the damper leads to an increase of the natural frequencies, and it is demonstrated that the maximum attainable damping is a certain fraction of the relative frequency increase, depending on the type of damping device. The asymptotic format only includes a real and a complex nondimensional parameter, and it is demonstrated how these parameters can be determined from the frequency increase by locking and from an energy balance on the undamped natural vibration modes. It is shown how the asymptotic format can incorporate sag of the cable, and specific results are presented for viscous damping, the effect of stiffness and mass, fractional viscous damping, and a nonlinear viscous damper. The relation of the stiffness component to active and semiactive damping is discussed.     

Magnetic fluid separation of gold-containing products in the vibration field

New data on the process of the magnetic fluid (MF) separation, which is based on the ponderomotive effect of the magnetized separation medium—or ferromagnetic fluid (FMF) on the nonmagnetic bodies arranged in it—are obtained. The magnetization of the FMF in a nonuniform magnetic field increases the strength of the field of mass forces affecting the FMF and, as a consequence, the pressure gradient in the FMF. This phenomenon can be considered pseudoweighting of the FMF and, when controlling the magnetic field force, it can be used to separate nonmagnetic materials according to their specific weights. The behavior of the FMF in the vibration field is investigated theoretically, and the dependence of energy absorbed by it on the amplitude and frequency of vibrations is revealed. Under industrial conditions, a series of tests on separation of free gold from the products of washing the goldfields by the method of MF separation is performed. The results of these tests prove the prospects of including secondary Au-containing resources with difficult-to-recover gold in processing.