• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Lateral response of pile foundations in lique fi able soils

    2015-02-10 02:45:45AsskarJanalizadehAliZahmatkesh

    Asskar Janalizadeh,Ali Zahmatkesh

    Babol University of Technology,Babol,Iran

    Full length article

    Lateral response of pile foundations in lique fi able soils

    Asskar Janalizadeh,Ali Zahmatkesh*

    Babol University of Technology,Babol,Iran

    A R T I C L E I N F O

    Article history:

    Received 6 March 2015

    Received in revised form

    5 May 2015

    Accepted 6 May 2015

    Available online 11 June 2015

    Pile foundations Lateral spreading Liquefaction Pseudo-static method

    Liquefaction has been a main cause of damage to civilengineering structures in seismically active areas. The effects of damage ofliquefaction on deep foundations are very destructive.Seismic behavior ofpile foundations is widely discussed by many researchers for safer and more economic design purposes.This paper presents a pseudo-static method for analysis of piles in lique fi able soil under seismic loads.A freefi eld site response analysis using three-dimensional(3D)numerical modeling was performed to determine kinematic loads from lateral ground displacements and inertial loads from vibration of the superstructure.The effects of various parameters,such as soil layering,kinematic and inertial forces, boundary condition of pile head and ground slope,on pile response were studied.By comparing the numerical results with the centrifuge test results,it can be concluded that the use of the p-y curves with various degradation factors in lique fi able sand gives reasonable results.

    ?2015 Institute of Rock and Soil Mechanics,Chinese Academy of Sciences.Production and hosting by Elsevier B.V.All rights reserved.

    1.Introduction

    The liquefaction is one of the challenging issues in geotechnical engineering and it damages structures and facilities during earthquakes.This phenomenon was reported as the main cause of damage to pile foundations during the major earthquakes(Kramer, 1996).In many earthquakes around the world,extensive damage to piles of bridges and other structures due to liquefaction and lateral spreading has been observed(Boulanger et al.,2003).Failures were observed in both sloping and level grounds and were often accompanied with settlement and tilting of the superstructure (Adhikari and Bhattacharya,2008).The loss of soil strength and stiffness due to excess pore pressure in lique fi able soilmay develop large bending moments and shear forces in the piles.If the residual strength of the lique fi able soil is less than the static shear stresses caused by a sloping site or a free surface such as a river bank,signi fi cant lateral spreading or downslope displacements may occur. The moving soil can exert damaging pressures against the piles, leading to failure(Finn and Fujita,2002).The performance of structures above piles depends widely on the behavior of pile foundations under earthquake loading.During past earthquakes, because of inadequacy of the pile to sustain large shear forces and bending moments,the extensive damage in lique fi able soil has been caused due to both lateral ground movement and inertial loads transmitted to piles.Under earthquake loading,the performance of piles in lique fi ed ground is a complex problem due to the effects ofprogressive buildup of pore water pressures and decrease ofstiffness in the saturated soil(Liyanapathirana and Poulos,2005). These effects involve inertialinteraction between structure and pile foundation,signi fi cant changes in stiffness and strength ofsoils due to increase of pore water pressures,large lateral loads on piles, kinematic interaction between piles and soils,nonlinear response of soils to strong earthquake motions,kinematic loads from lateral ground displacements,and inertial loads from vibration of the superstructure(Bradley et al.,2009;Gao et al.,2011).

    Various approaches including shaking table and centrifuge tests and also various numerical methods have been developed for the dynamic response analysis of single pile and pile group.The soilpile-structure interaction has been investigated using the centrifuge test(e.g.Finn and Gohl,1987;Chang and Kutter,1989;Liu and Dobry,1995;Hushmand et al.,1998;Wilson,1998;Abdoun and Dobry,2002;Su and Li,2006)and shaking table test(e.g.Mizuno and Liba,1982;Yao et al.,2004;Tamura and Tokimatsu,2005; Han et al.,2007;Gao et al.,2011;Haeri et al.,2012).The obvious advantage of shaking table and centrifuge tests is the ability to obtain detailed measurements of response in a series of tests designed to physically evaluate the importance of varying earthquake characteristics(e.g.levelof shaking,frequency content),soil pro fi le characteristics,and/or pile-superstructure characteristics (Wilson,1998).However,some limitations exist in centrifuge tests, for example,sand grains in centrifuge tests correspond to bigger gravel particles in prototype(Towhata,2008).

    To simulate the piles in lique fi able soil layers,Finn and Fujita (2002),Klar et al.(2004),Oka et al.(2004),Uzuoka et al.(2007),Cheng and Jeremic(2009),Comodromos et al.(2009),and Rahmani and Pak(2012)used three-dimensional(3D)fi nite element method.The complexity and time-consuming nature of 3D nonlinear fi nite element method for dynamic analysis makes it useful only for very large practical projects or research and not feasible for engineering practice.However,it is possible to obtain reasonable solutions for nonlinear response of pile foundations with fewer computations by relaxing some of the boundary conditions in full 3D analysis(Finn and Fujita,2002).

    The simple approach for modeling and simulation of the piles in lique fi ed grounds is based on scaling of p-y springs,where p and y are the soil resistance per unit length of the pile and pile lateral displacement,respectively.Because of complexity and timeconsuming of two-dimensional(2D)and 3D numerical modeling, most of the designers and researchers prefer to use onedimensional(1D)Winkler method based on fi nite element or fi nite difference method for the seismic analysis ofpile foundations. In pseudo-static method,a static analysis is carried out to obtain the maximum response(de fl ection,shear force and bending moment) developed in the pile due to seismic loading.In Winkler models,p-y curves are used to de fi ne the behavior of the nonlinear spring at any depth.These p-y curves can be obtained from the results of model tests or fi eld(Liyanapathirana and Poulos,2005).The Winkler assumption is that the soil-pile interaction resistance at any depth is related to the pile shaft displacement at that depth only,independent of the interaction resistances above and below(Wilson, 1998).

    This pseudo-static method has been suggested early by Miura et al.(1989),Miura and O’Rourke(1991),Liu and Dobry(1995), JRA(1996),AIJ(1998)and recently by Liyanapathirana and Poulos (2005)and Elahi et al.(2010).This method for pile seismic analysis sometimes underestimates,and sometimes overestimates shears,moments and de fl ection of the piles.However,in many practical conditions,the results of pseudo-static method are reasonable(Tabesh,1997).

    In this paper,a pseudo-static method has been applied for estimation of the response of pile during dynamic loading.First,de finition of the geometry and the soil modeling parameters are presented.Next,the numerical model is verti fi ed by means of the centrifuge test.And then the effects ofvarious parameters,including soil layering,kinematic and inertial forces,boundary condition of pile head and ground slope,on the behaviors of piles are studied.

    2.Numerical analysis

    All simulations were conducted using the open-source computational platform OpenSees(McKenna and Fenves,2007).This platform allows for developing applications to simulate the performance ofstructuraland geotechnicalsystems subjected to static and seismic loadings.In this paper,the steps for calculation of pile response are summarized as follows:

    (1)A free-fi eld site response analysis was performed during the dynamic loading using 3D numerical modeling.From this analysis,time history of ground surface acceleration and the maximumground displacement along the length ofthe pile can be calculated.

    (2)The dynamic analysis was performed using the time history of ground surface acceleration calculated in Step 1 for pile length above ground and superstructure with a fi xed base.From this analysis,the maximum acceleration of superstructure can be calculated.

    (3)In 1D Winkler analysis,the maximum soildisplacement pro fi le calculated in Step 1 and the maximum acceleration of superstructure in Step 2 were applied to the pile as shown in Fig.1.

    First,the time history of the ground surface acceleration and the maximum ground displacement at each depth were obtained from the free-fi eld site response analysis.Taboada and Dobry(1993)and Gonzalez et al.(2002)showed that the pore pressure time histories recorded at the same elevation are identical,indicating the 1D behavior of the model.In free-fi eld analysis,the modelconsists ofa single column of 3D brick elements.The soil layers were modeled using cubic 8-noded elements with u-p formulation in which each node has four degrees of freedom:three for soil skeleton displacements and one for pore water pressure.To consider the effect of the laminar box in the numerical simulation,nodes at the same depths were constrained to have equal displacements in the horizontal and vertical directions.The pore water pressures were allowed to freely develop for all nodes except those at the surface and above the water table.The bottom boundary was assumed fi xed in all directions.

    The materialmodel plays a key role in the numericalsimulation of the dynamic behavior of lique fi able soils.The model in Dafalias and Manzari(2004),a critical state two-surface plasticity model, was used in this paper.This model requires fi fteen material parameters and two state parameters to describe the behavior ofsands and has been amply tested for simulating the behavior of granular soils subjected to monotonic and cyclic loadings(Jeremic et al., 2008;Taiebat et al.,2010;Rahmaniand Pak,2012).The key advantages of the model are that(1)it is relatively simple and(2)it has a unique calibration of input parameters.Thus,a single set of parameters independent of void ratio and effective consolidation stress levelwas used for the Dafalias and Manzari’s materialmodel. Table 1 presents the material parameters for Nevada sand.The additional parameters used for free-fi eld analysis are presented in Table 2.Itcan be noted thatatthe onsetofliquefaction,change ofsoil particles creates additionalpathways for water.This leads to a signi fi cant increase in permeability coef fi cient(Rahmani et al.,2012). In this study,the permeability coef fi cient value was increased 10 times the initialvalue(suggested by Rahmani et al.(2012)).

    For free-fi eld analysis,the simulations were carried out in two loading stages.At the fi rst stage,the soil skeleton and pore water weight were applied to soilelements.The values ofstress and strain in this stage were used as initialvalues for the next stage of loading. At the second stage,dynamic analysis was performed by application of an input motion to the model base.

    Fig.1.A beam on the nonlinear Winkler foundation(BNWF)model for pseudo-static analysis.

    Table 1 Materialparameters for Nevada sand(Rahmaniand Pak,2012).

    Table 2 Additional parameters for pseudo-static and free-fi eld analysis(Wilson,1998; Rahmani and Pak,2012).

    In the second step,after free-fi eld analysis,pile length above ground and superstructure were modeled.The pile was modeled as beam column elements with elastic section properties.The superstructure was modeled at the pile head.Generally,the superstructures above the pile foundations are multi-degree of freedom systems,but in the design of pile foundations,the superstructure was modeled as a single mass at the pile head to simplify the analysis.In this step,the base model was also fi xed.

    The modelconsidered for the third step(pseudo-static analysis) is shown in Fig.1.There are two versions ofthe pseudo-static BNWF method.These two methods are different in the way in which the lateral load on pile due to ground movement(kinematic load)is considered.The fi rst BNWF requires free-fi eld soilmovements as an input.The free-fi eld soil displacements are imposed on the free ends of the p-y springs due to lateraldilatation layers.In the second BNWF,the limit pressures over the depth of the lateral spreading soilwere applied and the p-y springs were removed in this interval. In case of limit pressure,interaction of the soil and pile was not modeled,because the analysis is simple and can be done by hand calculation.Inertia forces from superstructure are represented as static forces applied simultaneously with lateral spreading demands.When limit pressures are applied directly to the pile nodes, bending moments and cap displacements depend on acceleration records and are greatly overpredicted for smallto mediummotions. This can be explained by the fact that the lateral spreading displacements were not large enough to mobilize limit pressures and actualpressures are smaller than limit pressures.However,for large motions,the pile cap displacements were considerably underpredicted(Brandenberg et al.,2007).In this paper,the free-fi eld soil movement was used as an input.The cap mass,multiplied by the maximum acceleration of the superstructure obtained from Step 2 as a lateral force(F),was applied at the pile head.The material properties of p-y curves for non-lique fi ed sand were computed based on API(1987).These curves are de fi ned by the following equation:

    where Puis the ultimate bearing capacity at depth z,K is the initial modulus of subgrade reaction,and y is the lateral de fl ection.The initial tangent stiffness(Kin),based on Eq.(1),is obtained as Kin=Kz.The p-y curves were modeled as zero-length elements with PySimple1 materials.Under dynamic loading,the piles are in fl uenced by kinematic loads from lateral ground displacements and inertial loads from vibration of the superstructure.Fig.1 shows an idealized schematic of the BNWF model for kinematic(F)and inertial(Δs)loads.The loss of bearing capacity for piles in loose sandy soils(particularly vulnerable to liquefaction and lateral spreading during dynamic loading)also occurred.Therefore,the excessive forces imposed on the foundation due to ground displacement led to shearing ofthe piles and subsequent structural collapse of the superstructure.There are three methods for considering the in fl uence of liquefaction on p-y curves in sand.In the fi rst case,the lateralresistance oflique fi able sand is assumed to be zero.This method can lead to large design responses and high construction costs which may be very conservative(Rollins et al., 2005).Another approach is to treat lique fi able sand as undrained soft clay and use the p-y curves for soft clay.The undrained shear strength used in this case is obtained as a ratio of undrained shear strength to initialeffective overburden stress using fi eld data,and it is a function ofoverburden stress and relative density(Rollins et al., 2005;Varun,2010).The third method for the simulation of pile response in lique fi able soils is the use ofreduction factors,called pmultipliers.The p-y curves in lique fi able sand are multiplied by a factor usually between 0.01 and 0.3 to decrease the strength ofsand due to liquefaction(Rollins et al.,2005;Brandenberg et al.,2007; Varun,2010).In this paper,the third method(p-multipliers)was used.The free-fi eld soil displacement and lateralforce in head pile were imposed incrementally using a static load control integrator.

    3.Validation of the proposed method

    The performance and ability of the proposed approac h to simulate pile behavior in lique fi able soil have been demonstrated by comparison between the numerical simulations and centrifuge tests performed by Wilson(1998).In these tests,the soil pro fi le consisted of two horizontal layers of saturated uniformly graded Nevada sand(see Fig.2).On prototype scale,the lower layer was 11.4 m thick with relative density of80%(dense)and the upper was 9.1 mthick with relative density of35%(loose).The single pile was asteelpipe of 0.67 m in diameter,18.8 min length,and 19 mm in wall thickness.The pile tip was about 3.8 m above the container base. The superstructure mass(Ms)was 49.1 Mg.Properties of Nevada sand with Dr=35%and 80%are presented in Table 2.The Kobe acceleration record(Fig.3)was used as an input to shake model.

    Fig.2.Layout of the model for centrifuge test by Wilson(1998).

    It is important to specify stiffness and lateralresistance ofthe py curves in lique fi able soil as explained in Section 2.Three cases were considered to evaluate the effects ofvariations in stiffness and lateral resistance of the p-y curves on the testing results.These cases include:(1)use of the p-y curves without the in fl uence of liquefaction;(2)use of the p-y curves with a constant degradation factor in lique fi ed sand;and(3)use of the p-y curves with various degradation factors in lique fi able sand.In case(2),different degradation factors can be considered between 0.05 and 0.5 to reduce the strength of lique fi able soil.In this case,because of the smallstrength oflique fi able soil,especially that of surface ground,a degradation factor of 0.1 was considered.In case(3),variation in degradation factor with depth was taken from a smallvalue(top of lique fi able layer)to 1.0(bottom of lique fi able layer).An exponential decay function from bottom to top of lique fi ed layer is proposed as

    where R is a degradation factor,as a function of distance from the top of the layer(z);H is the lique fi able layer thickness;and R0is the degradation factor at the top of the layer.Both the ultimate resistance and initial stiffness of the p-y curves in the lique fi able layer were taken to be R%of their unreduced magnitudes.

    In free-fi eld analysis under the Kobe acceleration record scaled to 0.04g,0.12g and 0.22g,because of the ground level,lateral displacement(Δs)was less than about 5 cm and the effects of kinematic loads on seismic response are small.Therefore,in these tests,it is reasonable to only consider the inertial loads for calculation of the pile response.The superstructure displacement was calculated using acceleration of the superstructure obtained from the tests carried out by Wilson(1998).The comparison between the observed and simulated results of superstructure displacement is shown in Fig.4.In cases(2)and(3),by considering the centrifuge test results,the best predictions were obtained with degradation factor of 0.1.These results clearly illustrate that the performance and accuracy of BNWF mainly depend on the accuracy in selection of the correct curves.As seen in Fig.4,the de fl ections observed during the centrifuge test are much larger than those simulated without the in fl uence of liquefaction.This is due to the loss of bearing capacity for the piles in loose and medium sandy soils during dynamic loading.It can be noted that use of constant degradation factors at various depths gives unreasonable results. The pile head displacements have a good agreement with the available experimental data in case(3).

    Fig.3.Acceleration record of Kobe earthquake scaled to 0.22g used in the centrifuge test by Wilson(1998).

    Fig.4.Comparison of superstructure displacement in various cases with the centrifuge tests by Wilson(1998).

    Fig.5 compares the maximum bending moment recorded from the centrifuge tests with that obtained frompresent analysis.In this fi gure,variation of R is similar to case(3)and Fig.5 shows that the results obtained from the numerical analysis agree with the values recorded during the centrifuge test.It can be said that increasing the value of R from a smallvalue at the top of the layer to 1.0 at the bottom of the layer produces a reasonable response for the piles. Therefore,this method was used for subsequent analysis.

    4.Results and discussion

    Fig.5.Comparison of bending moment pro fi les with the centrifuge tests by Wilson (1998).

    In this section,the behavior of pile for various conditions is discussed.The soil pro fi le is the same as the one used in the centrifuge test performed by Wilson(1998).The pro fi le has two layers:the upper layer is lique fi able(relative density of 35%)whilethe lower layer(relative density of 80%)is not.Three different ground slopes of 1%,2%,and 4%were considered.The water table was supposed to be 1 m,2 m,3 m,and 4 m below the ground.This means that the thicknesses ofnon-lique fi able surface crust are 1 m, 2 m,3 m,and 4 m.The input motion for the model was a 20-cycles sinusoidalwave with a frequency of 2 Hz and the peak acceleration of 0.5g.It should be noted that the intensity,frequency content(e.g. predominant period)and the duration of strong shaking are important characteristics of an earthquake(Rathje et al.,1998). These characteristics affected the response of piles.The pile responses largely depend on the shaking amplitude.Increase in the shaking amplitude(because of more reduction of restraint on lique fi ed soil)resulted in a decrease in the restraint against bending under the lateral load,and the maximum bending moment in piles signi fi cantly increased(Gao et al.,2011).The frequency also had a signi fi cant effect on pile response.

    The free-fi eld analysis showed that displacement of level ground is signi fi cantly less than that of the sloping ground.Fig.6 compares the displacement of sloping ground when the thickness of the non-lique fi able surface layer is 1 m and 2 m with ground slope of 2%.This fi gure highlights the importance of non-lique fi able surface layer as a key parameter on ground displacement.When liquefaction occurs in sloping ground,because of displacements developing up to several meters,large lateral forces may act on the pile.This phenomenon is commonly called lateral spreading(Klar et al.,2004).In lateral spreading,the driving forces only exceed the resisting forces during those portions of the earthquake that impose net inertialforces in the downslope direction.Each cycle of net inertial forces in the downslope direction causes the driving forces to exceed the resisting forces along the slip surface,resulting in progressively and incrementally lateral movement(Day,2002). Based on the results of free-fi eld analysis,the displacement pro fi le can be matched with constant displacement across the upper soil layer,a linear variation across the lique fi able and non-lique fi able layers.

    Fig.6.Displacement pro fi le at two different thicknesses of non-lique fi able surface layer when the ground slope is 2%.

    Fig.7.Variations in bending moment along the pile in different sloping grounds(free head,without superstructure).

    The variations in bending moment along piles in different ground slopes for various conditions are presented in Figs.7-10. The results show that in sloping grounds,when a non-lique fi able soil layer overlies a lique fi able soil layer and piles are embedded in the non-lique fi able soil layer,the lateral spreading has more in fl uences on the damage of piles.An increase in bending moment occurred as the ground slope increased.After liquefaction,if the static shear stress caused by sloping ground is more than the shear strength of lique fi able soil,the non-lique fi able surface crust overlying a lique fi ed soil layer can slide with a considerable amount of displacement.In this condition(lateral spreading),the nonlique fi able surface layer was carried along with the underlying fully lique fi able soil and a large lateral force was imposed on the embedded piles(Ashour and Ardalan,2011).This force due to the lateral movement of the non-lique fi able layer has the potential to induce large bending moments in the piles leading to failure.

    Fig.8.Variations in bending moment along the pile in different sloping grounds(fi xed head,without superstructure).

    Fig.9.Variations in bending moment along the pile in different sloping grounds(free head,with superstructure).

    The boundary condition of the pile head has an important effect on the pile responses(moments,shear and de fl ections).In layered soil deposits,a lique fi able soil layer is overlain by a non-lique fi able layer;when the pile head is free,the maximum bending moment develops at a depth corresponding to the interface of lique fi able and non-lique fi able layers(see Figs.7 and 9).When the pile head is fi xed,the maximum bending moment develops at two locations: (1)at the pile head and(2)at the interface of the two layers(see Figs.8 and 10).

    Fig.10.Variations in bending moment along the pile in different sloping grounds (fi xed head,with superstructure).

    The dynamic effects during earthquake on deep foundations are critically important.These effects include the kinematic forces applied by the soil to the pile foundation and the inertial forces of the superstructure due to earthquake.The combination of cyclic horizontal kinematic loads due to ground displacements and inertialloads from the superstructure determines the criticalload for piles during the shaking phase(Cubrinovski et al.,2009).The kinematic loads depend on the magnitude of ground deformations and the stiffness of the soil during a given loading cycle.Due to the in fl uence of liquefaction on free-fi eld soil response and soil-pilestructure interaction,the magnitudes of inertial and kinematic loads change(Han et al.,2007).When the acceleration of ground surface or superstructure mass is large or lateral dynamic stiffness of pile group due to pile and/or soil stiffness is small,the inertial effects may become important(Elahi et al.,2010).Figs.7 and 8 illustrate that in the absence of a superstructure,the maximum bending moment near the pile head decreases signi fi cantly,and a major difference is observed at that location between piles with and without superstructure,but the values are approximately unchanged at large depths(Figs.9 and 10).In other words,the same kinematic forces were developed in the piles with and without superstructure.At greater depths,where the inertial effects from the superstructure are less signi fi cant,pile damage can occur due to the lateral loads arising from lateral spreading(the excessive ground movement).Both inertial and kinematic loads can cause damages at the pile head.

    The inertia effects of the superstructure before development of the pore water pressures and liquefaction are important and the kinematic effects can often be neglected.Ishihara(1997)stated that inertial forces are the cause of development of the maximum bending moment near the pile head.These forces are predominant before liquefaction and are con fi rmed by the results of Figs.7-10.If the shaking continues after liquefaction,the inertial forces are combined with kinematic forces on the pile foundation arising from large cyclic ground deformations.It can be said that the kinematic loading in the areas oflateralspreading with relatively strong nonlique fi able surface layers is important.Then,the pile failure near the bottom of the lique fi able layer is likely in fl uenced by kinematic loads from the lique fi able layer,while failure near the pile head is likely in fl uenced by inertial loads from the superstructure and kinematic loads from the non-lique fi able layer.It should be noted that the bending moments are underpredicted when the structural inertia forces are neglected.

    The pile response is sensitive to the thickness of the nonlique fi able surface layer(H)and thickness of the lique fi able layer (L).Fig.11 shows variations in pile head displacement when the ratio of thickness(H/L)for different sloping grounds is increased. The pile head displacement increases to a peak value and then decreases with subsequent increase in thickness ratio.As the thickness of the non-lique fi able surface layer increases,horizontal kinematic loads due to ground displacements also increase andresult in higher pile head displacement.However,when the nonlique fi able surface layer is thick,because of increasing effective stress on the lique fi able layer,ground displacements are decreased, resulting in smaller values of pile head displacement.In addition, when the non-lique fi able surface layer is thick,due to increase of static shear stress in ground with great slopes,ground displacement and pile head displacement also increase.

    Fig.11.Variations in pile head displacement against ratio of thickness(H/L).

    5.Conclusions

    This paper presents a method for analysis of piles in lique fi able soilunder seismic loads.Three steps for calculation ofpile response were:(1)free-fi eld response analysis using 3D numericalmodeling for calculation of ground surface acceleration and the maximum ground displacement along the length of the pile,(2)the dynamic analysis of pile length above ground and superstructure for calculation of the maximum acceleration of superstructure,and(3)1D Winkler analysis for calculation of pile response.All simulations in three steps were conducted using the open-source computational platform OpenSees.After veri fi cation ofthe numericalmodelusing a centrifuge test,analyses were carried out for various conditions.

    By comparing the numerical results with the centrifuge test,it can be concluded that using the p-y curves with various degradation factors in lique fi able sand produces reasonable results.In addition,the non-lique fi able surface layer especially in sloping ground plays a key role in ground displacement.When the pile head is free,the maximum bending moment develops at a depth corresponding to the interface of lique fi able and non-lique fi able layers.When the pile head is fi xed,there are two locations for developing the maximum bending moment:(1)at the pile head and(2)at the interface of the two layers.Moreover,at greater depths,where inertial effects from the superstructure are less signi fi cant,pile damage may occur due to lateralloads arising from lateral spreading.Both inertial and kinematic loads can cause damages at the pile head.

    Con fl ict of interest

    The author con fi rms that there are no known con fl icts of interest associated with this publication and there has been no signi fi cant fi nancial support for this work that could have in fl uenced its outcome.

    Abdoun T,Dobry R.Evaluation ofpile foundation response to lateralspreading.Soil Dynamics and Earthquake Engineering 2002;22(9-12):1051-8.

    Adhikari S,Bhattacharya S.Dynamic instability of pile-supported structures in lique fi able soils during earthquakes.Shock and Vibration 2008;15(6):665-85. American Petroleum Institute(API).Recommended practice for planning,designing and constructing fi xed offshore platforms(RP 2A-WSD).Washington,DC,USA: API;1987.

    Architectural Institutive of Japan(AIJ).Recommendation for design of building foundations.1998.

    Ashour M,Ardalan H.Piles in fully lique fi ed soils with lateral spread.Computers and Geotechnics 2011;38(6):821-33.

    Boulanger RW,Kutter BL,Brandenberg SJ,Singh P,Chang D.Pile foundations in lique fi ed and laterally spreading ground during earthquakes:centrifuge experiments and analyses.Davis,USA:Center for Geotechnical Modeling,University of California;2003.Report No.UCD/CGM-03/01.

    Bradley BA,Cubrinovski M,Dhakal RP,MacRae GA.Intensity measures for the seismic response of pile foundations.Soil Dynamics and Earthquake Engineering 2009;29(6):1046-58.

    Brandenberg J,Boulanger RW,Kutter BL,Chang D.Static pushover analyses of pile groups in lique fi ed and laterally spreading ground in centrifuge tests.Journalof Geotechnical and Geoenvironmental Engineering 2007;133(9):1055-66.

    Chang GS,Kutter BL.Centrifugal modeling of soil-pile-structure interaction.In: Engineering Geology and Geotechnical Engineering:Proceedings of the 25th Symposium,Reno,Nevada.Rotterdam,the Netherlands:A.A.Balkema;1989. p.327-36.

    Cheng Z,Jeremic B.Numerical modeling and simulation of pile in lique fi able soil. Soil Dynamics and Earthquake Engineering 2009;29(11-12):1405-16.

    Comodromos EM,Papadopoulou MC,Rentzepris IK.Pile foundation analysis and design using experimental data and 3-D numerical analysis.Computers and Geotechnics 2009;36(5):819-36.

    Cubrinovski M,Ishihara K,Poulos H.Pseudo-static analysis of piles subjected to lateral spreading.Bulletin of the New Zealand Society for Earthquake Engineering 2009;42(1):28-38.

    Dafalias YF,Manzari MT.Simple plasticity sand model accounting for fabric change effects.Journal of Engineering Mechanics 2004;130(6):622-34.

    Day RW.Geotechnicalearthquake engineering handbook.New York,USA:McGraw-Hill;2002.

    Elahi H,Moradi M,Poulos HG,Ghalandarzadeh A.Pseudostatic approach for seismic analysis of pile group.Computers and Geotechnics 2010;37(1-2):25-39.

    Finn WDL,Fujita N.Piles in lique fi able soils:seismic analysis and design issues.Soil Dynamics and Earthquake Engineering 2002;22(9-12):731-42.

    Finn WDL,Gohl WB.Centrifuge modelstudies of piles under simulated earthquake loading from dynamic response of pile foundations-experiment,analysis and observation.Geotechnical Special Publication American Society of Civil Engineering(ASCE)1987;11:21-38.

    Gao X,Ling XZ,Tang L,Xu P.soil-pile-bridge structure interaction in liquefying ground using shake table testing.Soil Dynamics and Earthquake Engineering 2011;31(7):1009-17.

    Gonzalez L,Abdoun T,Sharp MK.Modelling of seismically induced liquefaction under high con fi ning stress.International Journal of Physical Modelling in Geotechnics 2002;2(3):1-15.

    Haeri SM,Kavand A,Rahmani I,Torabi H.Response of a group of piles to liquefaction-induced lateral spreading by large scale shake table testing.Soil Dynamics and Earthquake Engineering 2012;38:25-45.

    Han JT,Kim SR,Hwang JI,Kim MM.Evaluation of the dynamic characteristics of soil-pile system in lique fi able ground by shaking table tests.In:The 4th International Conference on Earthquake Geotechnical Engineering.Thessaloniki, Greece;2007.Paper No.1340.

    Hushmand B,Scott RF,Crouse CB.Centrifuge liquefaction tests in a laminar box. Geotechnique 1998;38(2):253-62.

    Ishihara K.Terzaghi oration:geotechnical aspects of the 1995 Kobe earthquake.In: Proceedings ofthe 14th International Conference on Soil Mechanics and Foundation Engineering.Rotterdam,the Netherlands:A.A.Balkema;1997.p.2047-73. Japanese Road Association(JRA).Speci fi cation for highway bridges,part V:seismic design.Tokyo,Japan:JRA;1996.

    Jeremic B,Cheng Z,Taiebat M,Dafalias Y.Numerical simulation of fully saturated porous material.InternationalJournalfor Numericaland Analytical Methods in Geomechanics 2008;32(13):1636-60.

    Klar A,Baker R,Frydman S.Seismic soil-pile interaction in lique fi able soil.Soil Dynamics and Earthquake Engineering 2004;24(8):551-64.

    Kramer SL.Geotechnical earthquake engineering.Upper Saddle River,USA:Prentice-Hall Inc.;1996.

    Liu L,Dobry R.Effect ofliquefaction on lateralresponse ofpiles by centrifuge model tests.NCEER Bulletin 1995;9(1):7-11.

    Liyanapathirana DS,Poulos HG.Seismic lateral response of piles in liquefying soil. Journal of Geotechnical and Geoenvironmental Engineering 2005;131(12): 1466-79.

    McKenna F,Fenves GL.The OpenSees command language manual.Version 1.2. Berkeley,USA:Paci fi c Earthquake Engineering Research Center,University of California;2007.http://opensees.berkeley.edu.

    Miura F,O’Rourke TD.Lateral spreading effects on pile foundations.In:Proceedings of the 3rd US-Japan Workshop on Earthquake Resistant Design of Lifeline Facilities and Countermeasures for Soil Liquefaction.MCEER;1991.

    Miura F,Stewart HE,O’Rourke TD.Nonlinear analysis of piles subjected to liquefaction induced large ground deformation.In:Proceedings of the 2nd US-Japan Workshop on Liquefaction,Large Ground Deformations and Their Effect on Lifelines.MCEER;1989.

    Mizuno H,Liba M.Shaking table testing ofseismic building-pile-soilinteraction.In: Proceedings of the 5th Japan Earthquake Engineering Symposium,Tokyo, Japan;1982.p.1713-20.

    Oka F,Lu CW,Uzuoka R,Zhang F.Numerical study of structure-soil-group pile foundations using an effective stress based liquefaction analysis method.In: Proceedings of the 13th World Conference on Earthquake Engineering,Vancouver,Canada;2004.Paper No.3338(CD-ROM).

    RahmaniA,Fare OG,Pak A.Investigation ofthe in fl uence ofpermeability coef fi cient on the numerical modeling of the liquefaction phenomenon.Scientia Iranica 2012;19(2):179-87.

    Rahmani A,Pak A.Dynamic behavior of pile foundations under cyclic loading in lique fi able soils.Computers and Geotechnics 2012;40:114-26.

    Rathje EM,Abrahamson NA,Bray JD.Simpli fi ed frequency content estimates of earthquake ground motion.Journal of Geotechnical and Geoenvironmental Engineering 1998;124(2):150-9.

    Rollins K,Gerber T,Lane J,Ashford S.Lateral resistance of a full-scale pile group in lique fi ed sand.Journal of Geotechnical and Geoenvironmental Engineering 2005;131(1):115-25.

    Su D,Li XS.Effect ofshaking intensity on seismic response ofsingle-pile foundation in lique fi able soil.In:Ground modi fi cation and seismic mitigation.ASCE;2006. p.379-86.

    Tabesh A.Lateral seismic analysis of piles.PhD Thesis.Sydney,Australia:Department of Civil Engineering,University of Sydney;1997.

    Taboada VM,Dobry R.Experimental results of model No.1 at RPI.In: Arulanandan K,Scott RF,editors.Veri fi cation of numerical procedures for the analysis of soil liquefaction problems.Rotterdam,the Netherlands:A.A.Balkema;1993.p.3-18.

    Taiebat M,Jeremic B,Dafalias YF,Kaynia AM,Cheng Z.Propagation ofseismic waves through lique fi ed soils.Soil Dynamics and Earthquake Engineering 2010;30(4): 236-57.

    Tamura S,Tokimatsu K.Seismic earth pressure acting on embedded footing based on large-scale shaking table tests.In:Seismic performance and simulation of pile foundations in lique fi ed and laterally spreading ground.Reston,USA: American Society of Civil Engineers;2005.p.83-96.

    Towhata I.Geotechnical earthquake engineering.Berlin-Heidelberg,Germany: Springer-Verlag;2008.

    Uzuoka R,Sento N,Kazama M,Zhang F,Yashima A,Oka F.Three-dimensional numerical simulation of earthquake damage to group-piles in a lique fi ed ground. Soil Dynamics and Earthquake Engineering 2007;27(5):395-413.

    Varun V.A non-linear dynamic macroelement for soilstructure interaction analyses of piles in lique fi able soils.PhD Thesis.Atlanta,USA:School of Civil and Environmental Engineering,Georgia Institute of Technology;2010.

    Wilson DW.Soil-pile-superstructure interaction in liquefying sand and soft clay. PhD Thesis.Davis,USA:University of California;1998.

    Yao S,Kobayashi K,Yoshida N,Matsuo H.Interactive behavior of soil-pilesuperstructure system in transient state to liquefaction by means of large shake table tests.Soil Dynamics and Earthquake Engineering 2004;24(5):397-409.

    Ali Zahmatkeshwas born in 1984 in Ferdows,Khorasan, Iran.He received his M.Sc.degree in Geotechnical Engineering from Mazandaran University,Iran in 2010.Since 2012,he is a Ph.D.student at Babol University of Technology.His thesis topic is Analysis of Performance of Pile Foundations in Lique fi ed Soils.His research interests cover soil improvement,soil liquefaction and numerical modeling.

    *Corresponding author.Tel.:+98 9158342367.

    E-mail address:A.zahmatkesh@stu.nit.ac.ir(A.Zahmatkesh).

    Peer review under responsibility of Institute of Rock and Soil Mechanics,Chinese Academy of Sciences.

    1674-7755?2015 Institute of Rock and Soil Mechanics,Chinese Academy of Sciences.Production and hosting by Elsevier B.V.All rights reserved.

    http://dx.doi.org/10.1016/j.jrmge.2015.05.001

    精品久久久久久久久亚洲| 亚洲av中文av极速乱| 国产成人精品久久久久久| 久久热精品热| 亚洲自拍偷在线| 97人妻精品一区二区三区麻豆| 一级黄色大片毛片| 日韩欧美精品v在线| 国产精品日韩av在线免费观看| 久久天躁狠狠躁夜夜2o2o| 91精品国产九色| 日韩一区二区视频免费看| 国产视频一区二区在线看| 久久精品久久久久久噜噜老黄 | 99久国产av精品| 中出人妻视频一区二区| 亚洲一区高清亚洲精品| 亚洲一区高清亚洲精品| 99久久成人亚洲精品观看| 精品一区二区三区av网在线观看| 少妇的逼水好多| 特级一级黄色大片| 欧美日韩精品成人综合77777| 俄罗斯特黄特色一大片| 男人狂女人下面高潮的视频| 国产亚洲精品av在线| 亚洲一区二区三区色噜噜| 免费在线观看影片大全网站| 亚洲真实伦在线观看| 波多野结衣高清作品| 国产成人影院久久av| 久久精品国产亚洲av天美| 亚洲av五月六月丁香网| 精品乱码久久久久久99久播| av女优亚洲男人天堂| 精品一区二区三区av网在线观看| 日韩av不卡免费在线播放| 变态另类丝袜制服| 国产中年淑女户外野战色| 一级毛片aaaaaa免费看小| 亚洲欧美精品自产自拍| 国产一区二区激情短视频| 精品一区二区三区视频在线| 午夜视频国产福利| 午夜影院日韩av| 最近视频中文字幕2019在线8| 亚洲成人久久爱视频| 99热精品在线国产| 国产av在哪里看| 久久久久久大精品| 最近2019中文字幕mv第一页| 午夜日韩欧美国产| 美女大奶头视频| 少妇裸体淫交视频免费看高清| 日本黄色片子视频| 亚洲,欧美,日韩| 国产三级在线视频| 亚洲aⅴ乱码一区二区在线播放| 国产v大片淫在线免费观看| 欧美+日韩+精品| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 三级毛片av免费| 真人做人爱边吃奶动态| 亚洲自拍偷在线| 日本爱情动作片www.在线观看 | 国产精品久久久久久精品电影| 亚洲精品影视一区二区三区av| 国产精品久久电影中文字幕| 不卡一级毛片| 噜噜噜噜噜久久久久久91| 桃色一区二区三区在线观看| 免费看光身美女| 久久国产乱子免费精品| 99久久无色码亚洲精品果冻| 日本在线视频免费播放| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 久99久视频精品免费| 看片在线看免费视频| 天天躁夜夜躁狠狠久久av| 国国产精品蜜臀av免费| 国产日本99.免费观看| 亚洲国产日韩欧美精品在线观看| 久久这里只有精品中国| 午夜福利视频1000在线观看| 免费av观看视频| 国产男靠女视频免费网站| 欧美性感艳星| 高清毛片免费看| 国产老妇女一区| 91麻豆精品激情在线观看国产| 高清毛片免费观看视频网站| 成人欧美大片| 美女大奶头视频| 少妇猛男粗大的猛烈进出视频 | 插逼视频在线观看| 精品久久久久久久久av| 成人鲁丝片一二三区免费| 亚洲最大成人中文| 免费电影在线观看免费观看| 一进一出抽搐动态| 伊人久久精品亚洲午夜| 久久婷婷人人爽人人干人人爱| 欧美成人一区二区免费高清观看| 亚洲丝袜综合中文字幕| 无遮挡黄片免费观看| 搡老妇女老女人老熟妇| 久久久a久久爽久久v久久| 我的老师免费观看完整版| 有码 亚洲区| 国产日本99.免费观看| 日本熟妇午夜| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点| 国产精品免费一区二区三区在线| 欧美色欧美亚洲另类二区| 熟妇人妻久久中文字幕3abv| 久久精品国产清高在天天线| 少妇人妻精品综合一区二区 | 变态另类成人亚洲欧美熟女| 国产精品日韩av在线免费观看| 欧美丝袜亚洲另类| 国语自产精品视频在线第100页| 亚洲综合色惰| 国产麻豆成人av免费视频| 中文资源天堂在线| 久久精品国产清高在天天线| 在线看三级毛片| 欧美激情国产日韩精品一区| 久久久久久久久久黄片| 国产69精品久久久久777片| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 99久久九九国产精品国产免费| 91麻豆精品激情在线观看国产| 国产淫片久久久久久久久| 在线国产一区二区在线| 搡女人真爽免费视频火全软件 | 国产白丝娇喘喷水9色精品| 高清日韩中文字幕在线| a级毛色黄片| 亚洲三级黄色毛片| 成人av一区二区三区在线看| 淫秽高清视频在线观看| 日本在线视频免费播放| 欧美又色又爽又黄视频| 亚洲高清免费不卡视频| 国产亚洲精品综合一区在线观看| 久久久久久伊人网av| 午夜免费激情av| 亚洲丝袜综合中文字幕| 1000部很黄的大片| 91av网一区二区| 亚洲熟妇熟女久久| 最近中文字幕高清免费大全6| 99热6这里只有精品| 精品一区二区三区av网在线观看| 日本成人三级电影网站| 五月伊人婷婷丁香| 超碰av人人做人人爽久久| 亚洲国产色片| 成人美女网站在线观看视频| 成人特级黄色片久久久久久久| 99久久精品热视频| 久久6这里有精品| 夜夜夜夜夜久久久久| 在线国产一区二区在线| 午夜福利在线观看免费完整高清在 | 国产午夜福利久久久久久| 精品久久久久久久久av| 亚洲国产精品国产精品| 午夜福利高清视频| 日本a在线网址| 国产 一区 欧美 日韩| 国产爱豆传媒在线观看| 六月丁香七月| 欧美一级a爱片免费观看看| 国产精品福利在线免费观看| 亚洲精品国产成人久久av| 国产美女午夜福利| 国产 一区 欧美 日韩| 内射极品少妇av片p| 亚洲美女黄片视频| 日韩强制内射视频| 别揉我奶头~嗯~啊~动态视频| 久久久精品大字幕| 欧美日韩一区二区视频在线观看视频在线 | 久久久国产成人免费| 亚洲欧美日韩高清在线视频| 小说图片视频综合网站| 国产精品,欧美在线| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 99久久久亚洲精品蜜臀av| 久久久精品欧美日韩精品| 国产v大片淫在线免费观看| 精品福利观看| 国产成年人精品一区二区| 少妇熟女欧美另类| .国产精品久久| 九九热线精品视视频播放| 成熟少妇高潮喷水视频| 此物有八面人人有两片| 免费观看的影片在线观看| 91在线观看av| 国产私拍福利视频在线观看| 精品日产1卡2卡| 精品福利观看| 少妇裸体淫交视频免费看高清| 国产高清三级在线| 一a级毛片在线观看| 亚洲欧美中文字幕日韩二区| 黄色日韩在线| 观看美女的网站| av在线老鸭窝| 在线观看一区二区三区| 国产黄色视频一区二区在线观看 | 欧美日韩国产亚洲二区| 十八禁网站免费在线| 久久久久性生活片| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 成人美女网站在线观看视频| 热99在线观看视频| 最近的中文字幕免费完整| 精品福利观看| 搡老妇女老女人老熟妇| 级片在线观看| 身体一侧抽搐| eeuss影院久久| 免费大片18禁| 久久人人爽人人片av| 久久久欧美国产精品| 国产男靠女视频免费网站| 日本免费a在线| 又爽又黄无遮挡网站| 亚洲色图av天堂| 热99在线观看视频| 一本久久中文字幕| 国产在线精品亚洲第一网站| 日本黄大片高清| 麻豆国产97在线/欧美| 国产人妻一区二区三区在| 亚洲成av人片在线播放无| 国产 一区精品| 毛片一级片免费看久久久久| 亚洲第一电影网av| 日本黄色视频三级网站网址| 国产一区二区在线av高清观看| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠久久av| 欧美成人一区二区免费高清观看| 国产精品国产高清国产av| 久久久精品94久久精品| 日韩大尺度精品在线看网址| 99热只有精品国产| 久久久久久久久久久丰满| 免费看a级黄色片| 毛片女人毛片| .国产精品久久| 一级毛片久久久久久久久女| 黄色欧美视频在线观看| 色吧在线观看| 直男gayav资源| 村上凉子中文字幕在线| 国产日本99.免费观看| 精品一区二区三区人妻视频| 亚洲精品久久国产高清桃花| 亚洲av五月六月丁香网| 国产高清有码在线观看视频| 91狼人影院| 九色成人免费人妻av| 亚洲图色成人| 日本与韩国留学比较| 国内精品宾馆在线| 国产精品久久久久久久电影| 岛国在线免费视频观看| 久久久久久久久久成人| 久久久精品欧美日韩精品| 成人综合一区亚洲| 五月玫瑰六月丁香| 尤物成人国产欧美一区二区三区| 在线观看一区二区三区| 午夜福利成人在线免费观看| 无遮挡黄片免费观看| 能在线免费观看的黄片| 国产女主播在线喷水免费视频网站 | 久久人人爽人人片av| 成年女人看的毛片在线观看| 国产成人a∨麻豆精品| 欧美+日韩+精品| 久久久久久大精品| 小蜜桃在线观看免费完整版高清| 深夜精品福利| 最近最新中文字幕大全电影3| 村上凉子中文字幕在线| 欧美丝袜亚洲另类| 亚洲真实伦在线观看| 国产人妻一区二区三区在| 国产精品久久久久久久电影| 高清毛片免费观看视频网站| 熟女电影av网| 男插女下体视频免费在线播放| 伦理电影大哥的女人| 综合色丁香网| 亚洲,欧美,日韩| 九九热线精品视视频播放| 国产色婷婷99| 国产亚洲91精品色在线| 联通29元200g的流量卡| 天堂√8在线中文| 午夜久久久久精精品| eeuss影院久久| 欧美丝袜亚洲另类| 午夜福利18| 99热6这里只有精品| 国产精品综合久久久久久久免费| 国产欧美日韩一区二区精品| 寂寞人妻少妇视频99o| 午夜福利在线观看吧| 少妇熟女aⅴ在线视频| 美女xxoo啪啪120秒动态图| 久久久欧美国产精品| 精品久久国产蜜桃| 99久久精品一区二区三区| 国产高清视频在线观看网站| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 久久久久久久亚洲中文字幕| 亚洲av美国av| 亚洲人成网站高清观看| 久久久久久久久久黄片| 亚洲欧美日韩高清在线视频| 日韩三级伦理在线观看| 黄色欧美视频在线观看| 亚洲内射少妇av| 国产高清三级在线| 精品午夜福利在线看| 久久久a久久爽久久v久久| 亚洲av.av天堂| 在线播放国产精品三级| 91麻豆精品激情在线观看国产| 日本三级黄在线观看| avwww免费| 亚洲欧美精品综合久久99| 午夜福利18| 日本黄色片子视频| 日本色播在线视频| 欧美另类亚洲清纯唯美| 草草在线视频免费看| 中出人妻视频一区二区| 男女下面进入的视频免费午夜| 免费看日本二区| 人人妻人人澡欧美一区二区| 一级av片app| 亚洲真实伦在线观看| 韩国av在线不卡| 亚洲精品在线观看二区| 国产69精品久久久久777片| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 男女啪啪激烈高潮av片| 国产视频内射| 亚洲成人中文字幕在线播放| 校园春色视频在线观看| 久久这里只有精品中国| 老女人水多毛片| 丰满的人妻完整版| 国产在视频线在精品| 99久久精品热视频| 国产高清三级在线| 国产黄a三级三级三级人| 啦啦啦观看免费观看视频高清| 国产精品人妻久久久影院| 国产成人a区在线观看| 天堂动漫精品| 午夜视频国产福利| 欧美一区二区亚洲| 国产精品国产高清国产av| 免费人成视频x8x8入口观看| 精品福利观看| 亚洲欧美清纯卡通| 一区二区三区高清视频在线| 一进一出抽搐动态| 午夜福利18| 国产精品伦人一区二区| 国产高清激情床上av| 国产精品永久免费网站| 久久久久久久久久黄片| 99久久精品一区二区三区| 亚洲熟妇中文字幕五十中出| 蜜桃亚洲精品一区二区三区| 久久久久久九九精品二区国产| 中文字幕熟女人妻在线| 亚洲天堂国产精品一区在线| 美女免费视频网站| 久久精品夜色国产| 精品福利观看| 久久久久久大精品| 十八禁国产超污无遮挡网站| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| av福利片在线观看| 午夜亚洲福利在线播放| 波多野结衣巨乳人妻| 99热只有精品国产| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久亚洲av鲁大| 亚洲欧美清纯卡通| 看十八女毛片水多多多| 99九九线精品视频在线观看视频| 国产色婷婷99| 婷婷六月久久综合丁香| 美女高潮的动态| 午夜福利高清视频| 春色校园在线视频观看| 网址你懂的国产日韩在线| 一个人免费在线观看电影| 亚洲人成网站在线观看播放| 久久久久久久久久黄片| 国产伦精品一区二区三区四那| 久久久精品94久久精品| 国产免费男女视频| videossex国产| 国产精品久久久久久久电影| 国产黄色视频一区二区在线观看 | 国产高潮美女av| 日韩 亚洲 欧美在线| 国产三级在线视频| 久久精品国产亚洲av天美| 97超视频在线观看视频| 国产激情偷乱视频一区二区| 俄罗斯特黄特色一大片| 亚洲经典国产精华液单| 亚洲欧美日韩高清专用| 亚洲精品在线观看二区| 成年av动漫网址| 久久久久久久久大av| 淫妇啪啪啪对白视频| 成年女人永久免费观看视频| 欧美人与善性xxx| 男女视频在线观看网站免费| 精品熟女少妇av免费看| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 中文字幕人妻熟人妻熟丝袜美| 久久人人爽人人爽人人片va| 97超碰精品成人国产| av在线播放精品| 免费av不卡在线播放| 午夜精品在线福利| 自拍偷自拍亚洲精品老妇| 中文字幕免费在线视频6| 久久精品国产清高在天天线| 看片在线看免费视频| 综合色丁香网| 国产亚洲精品久久久久久毛片| 久久婷婷人人爽人人干人人爱| 午夜久久久久精精品| 久久人人爽人人爽人人片va| 久久久久久久久久成人| 久久久久国内视频| 亚洲激情五月婷婷啪啪| 精品日产1卡2卡| 国产精品免费一区二区三区在线| 亚洲中文日韩欧美视频| 少妇猛男粗大的猛烈进出视频 | 丝袜喷水一区| 午夜a级毛片| 成熟少妇高潮喷水视频| 男人舔女人下体高潮全视频| 最新在线观看一区二区三区| 此物有八面人人有两片| 国产成人aa在线观看| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 国产精品久久久久久av不卡| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 熟妇人妻久久中文字幕3abv| 偷拍熟女少妇极品色| 尾随美女入室| 听说在线观看完整版免费高清| 免费av观看视频| 国产一区二区三区在线臀色熟女| 午夜免费男女啪啪视频观看 | 国产精品电影一区二区三区| 国产国拍精品亚洲av在线观看| 又爽又黄a免费视频| 男女做爰动态图高潮gif福利片| 高清毛片免费看| 久久这里只有精品中国| 中文字幕免费在线视频6| av国产免费在线观看| 男人狂女人下面高潮的视频| 国产精品无大码| 久久久国产成人精品二区| 欧美日本亚洲视频在线播放| aaaaa片日本免费| 日本五十路高清| 久久人人爽人人片av| 婷婷精品国产亚洲av在线| 久久人人爽人人片av| av天堂在线播放| 少妇熟女欧美另类| av天堂在线播放| 国产亚洲欧美98| 深夜a级毛片| 亚洲精品456在线播放app| 日韩欧美精品v在线| 人妻久久中文字幕网| 亚洲精品成人久久久久久| 中文字幕免费在线视频6| 久久精品国产亚洲av香蕉五月| 日韩欧美国产在线观看| 亚州av有码| av专区在线播放| 亚洲国产欧美人成| 搞女人的毛片| 久久久色成人| 3wmmmm亚洲av在线观看| 国产精品电影一区二区三区| 欧美在线一区亚洲| 一区二区三区高清视频在线| 欧美另类亚洲清纯唯美| 男女边吃奶边做爰视频| 老熟妇乱子伦视频在线观看| 国产精品美女特级片免费视频播放器| 日本 av在线| 亚洲丝袜综合中文字幕| 99久久无色码亚洲精品果冻| 日韩一本色道免费dvd| 如何舔出高潮| 国产成人福利小说| 三级经典国产精品| 在线免费观看不下载黄p国产| 久久精品国产自在天天线| 大又大粗又爽又黄少妇毛片口| 精品免费久久久久久久清纯| 1000部很黄的大片| 精品99又大又爽又粗少妇毛片| 亚洲精品粉嫩美女一区| h日本视频在线播放| 精品人妻熟女av久视频| 草草在线视频免费看| 在线国产一区二区在线| 18禁裸乳无遮挡免费网站照片| av中文乱码字幕在线| 精品人妻一区二区三区麻豆 | 一级毛片aaaaaa免费看小| 亚洲aⅴ乱码一区二区在线播放| 国产爱豆传媒在线观看| 国产精品永久免费网站| 久久久精品欧美日韩精品| 全区人妻精品视频| 成年女人毛片免费观看观看9| 国产精品人妻久久久影院| 亚洲最大成人中文| 国产成人aa在线观看| 久久精品国产亚洲av香蕉五月| 亚洲精华国产精华液的使用体验 | 成人美女网站在线观看视频| 日日干狠狠操夜夜爽| 亚洲成人久久性| 欧美极品一区二区三区四区| 久久久精品94久久精品| 精品不卡国产一区二区三区| 成人鲁丝片一二三区免费| 波多野结衣高清作品| 亚洲av五月六月丁香网| 韩国av在线不卡| 亚州av有码| 看免费成人av毛片| 乱系列少妇在线播放| 久久久成人免费电影| 国产午夜精品论理片| 好男人在线观看高清免费视频| 少妇丰满av| 少妇的逼水好多| 亚洲欧美清纯卡通| eeuss影院久久| 亚洲人成网站在线观看播放| 91在线精品国自产拍蜜月| 插逼视频在线观看| 在线免费观看不下载黄p国产| 无遮挡黄片免费观看| 精品久久久久久成人av| 国产私拍福利视频在线观看| 99久久精品一区二区三区| 精品免费久久久久久久清纯| .国产精品久久| 麻豆国产av国片精品| 亚洲精品亚洲一区二区| 又爽又黄无遮挡网站| 国产伦精品一区二区三区视频9| 亚洲在线自拍视频| 我要看日韩黄色一级片| 婷婷亚洲欧美| 亚洲最大成人手机在线| 亚洲激情五月婷婷啪啪| 晚上一个人看的免费电影| 99久久中文字幕三级久久日本| 日韩一区二区视频免费看| 午夜亚洲福利在线播放| 99视频精品全部免费 在线| 一进一出抽搐动态| 精品久久久久久久久亚洲| 亚洲在线自拍视频| 美女xxoo啪啪120秒动态图| 变态另类成人亚洲欧美熟女| 国产精品一区www在线观看|