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

    Regressive approach for predicting bearing capacity of bored piles from cone penetration test data

    2015-02-10 02:45:59IydAlkrooshMohmmdBhdoriHmidNikrzAlirezBhdori

    Iyd S.Alkroosh,Mohmmd Bhdori,Hmid Nikrz,Alirez Bhdori

    aDepartment of Civil Engineering,College of Engineering,Al-Qadissiya University,Al Diwaniyah,Iraq

    bSchoolof Soil Science and Engineering,University of Tehran,Tehran,Iran

    dSchoolof Environment,Science and Engineering,Southern Cross University,Lismore,New South Wales,Australia

    Full length article

    Regressive approach for predicting bearing capacity of bored piles from cone penetration test data

    Iyad S.Alkroosha,Mohammad Bahadorib,Hamid Nikrazc,Alireza Bahadorid,*

    aDepartment of Civil Engineering,College of Engineering,Al-Qadissiya University,Al Diwaniyah,Iraq

    bSchoolof Soil Science and Engineering,University of Tehran,Tehran,Iran

    cDepartment of Civil Engineering,Curtin University,Perth,Australia

    dSchoolof Environment,Science and Engineering,Southern Cross University,Lismore,New South Wales,Australia

    A R T I C L E I N F O

    Article history:

    Received 27 April2015

    Received in revised form

    12 June 2015

    Accepted 30 June 2015

    Available online 5 August 2015

    Bored piles Cone penetration test(CPT) Bearing capacity Least square support vector machine (LSSVM) Training Validation

    In this study,the least square support vector machine(LSSVM)algorithm was applied to predicting the bearing capacity of bored piles embedded in sand and mixed soils.Pile geometry and cone penetration test(CPT)results were used as input variables for prediction of pile bearing capacity.The data used were collected from the existing literature and consisted of 50 case records.The application of LSSVM was carried out by dividing the data into three sets:a training set for learning the problem and obtaining a relationship between input variables and pile bearing capacity,and testing and validation sets for evaluation ofthe predictive and generalization ability of the obtained relationship.The predictions ofpile bearing capacity by LSSVM were evaluated by comparing with experimental data and with those by traditionalCPT-based methods and the gene expression programming(GEP)model.It was found that the LSSVMperforms wellwith coef fi cientofdetermination,mean,and standard deviation equivalent to 0.99, 1.03,and 0.08,respectively,for the testing set,and 1,1.04,and 0.11,respectively,for the validation set.The low values of the calculated mean squared error and mean absolute error indicated that the LSSVMwas accurate in predicting the pile bearing capacity.The results of comparison also showed that the proposed algorithm predicted the pile bearing capacity more accurately than the traditional methods including the GEP model.

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

    1.Introduction

    Bearing capacity is one of the most important factors that govern the design of pile foundations.Therefore,it has been the subject of interest for many researchers throughout the history of the geotechnical engineering profession.As a result,numerous theoretical and experimental procedures have been proposed to predict the pile behavior and bearing capacity.However,accurate evaluation ofpile bearing capacity and certain interpretation ofpile load transfer mechanism are still far from being accomplished due to the complexity of the problem.

    The theoretical solutions,which employ the theory of bearing capacity to calculate the pile shaft and tip resistance,involve shortcomings resulting from considerable uncertainty over the factors that in fl uence the bearing capacity.Among those factors are the effect of installation method,stress history and soil compressibility.For bored piles embedded in layered soil,the problem is more complex due to sensitivity of the factors that affect the behavior of the pile and the dif fi culty in quantifying those factors. For instance,the friction angle between pile and the surrounding soil cannot be exactly determined because of the effect of installation procedure and the dif fi culty in fi nding the real soil properties.

    The experimental solutions that correlate the results of in-situ tests such as standard penetration test(SPT)or cone penetration test(CPT)with pile bearing capacity also involve setbacks.That may be attributed to that the SPT has substantially inherent variability and does not re fl ect soil compressibility(Abu-Kiefa,1998).Moreover,the SPT results are affected by many factors,such as operator, drilling,hammer ef fi ciency,and rate of blows.Hence,the accuracy of the proposed correlations between SPT data and pile bearing capacity is not assured.Although the correlation between pile capacity and CPT data can be a better alternative to the SPT correlation,comparative studies of the available CPT-based methods carried out by a number of researchers(e.g.Briaud,1988;Roberston et al.,1988;Eslami,1997;Abu-Farsakh and Titi,2004;Cai et al., 2008)have shown that the capacity predictions can be verydifferent for the same case depending on the method employed.It is also found that these methods cannot provide consistent and accurate prediction of pile bearing capacity.

    Considering the limitations of the proposed procedures for predicting pile bearing capacity and the limited success that they have achieved in terms of providing accurate prediction of pile bearing capacity,further research is required to overcome the complications associated with the problem.Arti fi cial intelligence techniques may be better alternatives,due to the capability ofbeing able to deal with complex and highly nonlinear functions,and employing the considerable capacity of computers to perform enormously iterated work.The modeling advantage of these techniques is their ability to capture the nonlinear and complex relationships between the targeted output and the factors affecting it, without having to assume a priori formula describing this relationship.A number of researchers(e.g.Teh et al.,1997;Abu-Kiefa, 1998;Das and Basudhar,2006;Ardalan et al.,2009;Shahin,201 Ornek et al.,2012;Tarawneh,2013)have successfully applie arti fi cial neural network(ANN),which is a form of arti fi cial inte ligence,to solving engineering problems.Genetic programmin (GP),which is another form of arti fi cialintelligence,has been use successfully in solving engineering problems(Rezania and Javad 2007;Alavi et al.,2011;Alkroosh and Nikraz,2011a,b;Gandom 2011;Gandomi and Alavi,2012).Recently,an emerging alg rithm,i.e.the least square support vector machine(LSSVM),which is a developed version of support vector machine(SVM),has been found successfulin solving engineering problems(Das et al.,2011a, b;Samui and Kothari,2011).This study investigates the feasibility of using the LSSVM to predict the bearing capacity of bored piles embedded in sand and mixed soils more accurately than the available methods.

    2.Support vector machine(SVM)

    The SVM is a method developed using the statistical learning concept(Suykens and Vandewalle,1999).It has been widely used across the world(Cortes and Vapnik,1995;Bazzani et al.,2001; Suykens et al.,2002;Amendolia et al.,2003;Baylar et al.,2009; übeyli,2010;Chen et al.,2011;Chamkalani et al.,2013;Ra fi ee-Taghanaki et al.,2013;Shokrollahi et al.,2013).

    If we have training samples with given data xi∈Rnand result data yi∈R with labels-1 and 1,respectively,the SVMestimates the function shown below(Suykens and Vandewalle,1999;Suykens et al.,2002):

    whereΦ(x)is the function that maps x,and w and b are the weight vector and bias variable.When the data are separable,we will have (Suykens and Vandewalle,1999;Suykens et al.,2002):

    Eq.(2)is nearly equal to(Das et al.,2011a,b;Chamkalani et al., 2013):

    The further development of linear SVM to non-independent case was also created by Cortes and Vapnik(1995).Simply,it is done by presenting extra variables into Eq.(3)(Suykens and Vandewalle,1999;Suykens et al.,2002):

    whereζkis the deviation factor.

    The optimalseparating hyperplane is predicted using the vector w that minimizes the functional conditions using the constraints (Eq.(4))(Suykens and Vandewalle,1999;Suykens et al.,2002; übeyli,2010):

    where p is the upper limit,and C is a coef fi cient.

    In the SVM,optimal separating hyperplane is calculated using the quadratic method(Cortes and Vapnik,1995):

    where a is the adjustable parameter,αiandβiare the Lagrange multipliers(Suykens and Vandewalle,1999;Suykens et al.,2002).

    In contrast to the SVM,the LSSVM is developed using minimization of cost equation(Suykens and Vandewalle,1999;Suykens et al.,2002):

    To derive the dual problem for the nonlinear classi fi cation problem of LSSVM,the Lagrange function is de fi ned as(Suykens and Vandewalle,1999;Suykens et al.,2002):

    The conditions for optimality can be obtained as

    By de fi ning ZT=[ΦT(x1)y1,ΦT(x2)y2,…,ΦT(xN)yN],Y=[y1,y2,…,yN],[1,1,…,1],ζ=[ζ1,ζ2,…,ζN],α=[α1,α2,…,αN],Eq.(10)is fi nally converted into the below form (Minoux,1986;Suykens and Vandewalle,1999):

    where INis an N×N identity matrix,andΩ∈RN×Nis the kernel matrix de fi ned by

    Most extensively used kernelfunctions are radialbasis function (RBF)(Eq.(13))and polynomial function(Eq.(14))(Gunn,1998; Müller et al.,2001):

    whereσis an independent variable.

    This work uses the RBF(Baylar et al.,2009;Xavier-de-Souza et al.,2009;Deng and Yeh,2010;Chamkalani et al.,2013;Ra fi ee-Taghanaki et al.,2013;Shokrollahi et al., 2013)kernel,which is suitable for the LSSVM mathematical modeling.

    3.Application of LSSVM for predicting pile bearing capacity

    3.1.Database used

    The application of LSSVM in this work is based on pile load test and CPT data,reported in Alsamman(1995),Eslami(1997)and Alkroosh and Nikraz(2011a,b).The case records were collected from a wide spectrum of geographic locations all over the world. Most of the piles embedded in layered soil of sand and clay.The piles were tested under slow maintained compression loads.According to Ng et al.(2004),the piles were classi fi ed into small diameter piles(diameter<600 mm)and large diameter piles (diameter>600 mm).

    Fig.1.Summary sheet for selected case record:(a)pile geometry and soil pro fi le;(b)cone tip resistance pro fi le;(c)load-movement plot(Alkroosh and Nikraz,2011a,b).

    3.2.Input and output variables

    A proper estimation of bearing capacity of pile foundation requires the identi fi cation of the factors that in fl uence the pile-soil interaction.An extensive study of relevant existing literature was carried out and concluded that pile geometry and soil properties are the signi fi cant factors that affect the pile bearing capacity.

    (1)Pile geometry

    Allgeotechnicalengineering sources con fi rm that pile diameter and length have signi fi cant in fl uence on bearing capacity of pile foundations.Therefore,these factors were selected to represent pile geometry for input of LSSVM model(Eslami,1997;Teh et al., 1997;Abu-Kiefa,1998).

    (2)Soil properties

    The soil properties that have principal in fl uence on pile bearing capacity are mainly interpreted as the angle of internal friction,φ, relative density,Dr,earth pressure coef fi cient,K0,and cohesion,c. Many researchers have suggested correlations between these factors and CPT data,and proposed solutions for prediction of pile bearing capacity.However,indirect correlations of CPT data with pile bearing capacity have been found to give inaccurate estimate of pile bearing capacity due to considerable uncertainty resulting from the intermediate steps and correlations.Consequently,direct correlations of CPT results with pile bearing capacity have been suggested and found to give better results.

    Thus,in this work,the CPT results were directly used as input variables based on the suggestions provided by the current CPT-based methods(e.g.Alsamman,1995;Eslami,1997).According to the methods,the cone point resistance measurements are averaged within pile tip in fl uence zone,which is a function of pile diameter, and used for calculating pile tip resistance.For smalldiameter piles, the in fl uence zone was taken based on the work done by Eslami (1997),while for large diameter piles the in fl uence zone was taken based on Alsamman(1995).The average cone point resistance along pile shaft was included as input,as most ofthe methods also suggest it to calculate pile shaft resistance.

    In fact,the use of sleeve friction in calculation of bearing capacity is not recommended by number of commonly used CPT-based methods,such as Laboratoire Central des Pont et Chausees (LCPC),due to the unreliability of sleeve friction measurements (Alkroosh and Nikraz,2011a,b).

    The mechanical cone point resistance was converted to equivalent electric values using Eq.(15)which was developed by Kulhawy and Mayne(1990):

    where q c is the cone point resistance;p a is the atmospheric pressure,which is 101 kPa.

    3.3.Pile bearing capacity,Qu

    The interpreted failure load,Qu,was considered as the pile bearing capacity and calculated according to Alsamman(1995).The failure load was de fi ned as the measured axial load that corresponds to a displacement equal to 5%of pile diameter plus the elastic compression of the pile.Fig.1 presents the de fi nition of the failure load for a case record selected from the database.

    3.4.LSSVM for predicting pile bearing capacity

    Table 1 Summary of database used in the application of LSSVM and the output results(Alkroosh and Nikraz,2011a,b).

    The data were mapped into range[-1,1]using below equation:where x and xnare the original and normalized desired variables, respectively;xminand xmaxare the extreme values of the variable x. This pre-processing procedure was applied to obtain the parameters of the LSSVM algorithm.Later,these values were changed to their original values.

    In the next step,the database was divided into three subsets including the training,testing and validation sets.The division of database into three subsets was normally performed randomly but considering the statisticalconsistency ofthe data ofthe subsets.For this purpose,70%,15%and 15%of the main dataset were randomly selected for building the LSSVM models(i.e.training,testing and validation sets,respectively).The allocated data for each subset are shown in Table 1.The effect of the percent allocation of the three subsets from the database on the accuracy of the fi nal model has been studied(Gharagheiziet al.,2011).As for the distribution of the data through the three subsets,we generally performed many distributions to avoid the local accumulations of the data in the feasible region of the problem.As a result,the acceptable distribution is the one with homogeneous accumulations of the data on the domain of the three subsets(Eslamimanesh et al.,2011).

    4.Results and discussion

    The ability of the LSSVMto predict the pile bearing capacity was evaluated primarily during the training phase by comparing the measured values of pile bearing capacity with the predicted values obtained by the LSSVM.When training phase was completed,it was expected that the proposed algorithm would correctly reproduce the targeted output values presuming the error is minimal.However,this was insuf fi cient to conclude that the trained LSSVM can predict the pile bearing capacity accurately as this technique,in general,has a high tendency towards over-fi tting.Over-fi tting refers to the large error in predictions when new data are presented to the trained LSSVM(i.e.ability of the LSSVM to memorize rather than generalize the form of the relationship between input and output data).Generally,over-fi tting is expected when data points in training sets are scanty(Das,2013).To overcome this setback,a testing dataset,not included in training data,was used to examine the predictive ability of LSSVM.If it was found to give accurate predictions,further evaluation was made to verify its generalization ability.This was done by examining the predictions of LSSVM through the use of validation dataset which was neither included among training nor testing datasets.Once the training testing and validation steps were successfully completed,the trained LSSVM can be recommended as a design tool for this type of application.

    The predictions of pile bearing capacity by the LSSVMas well as other CPT-based methods including Schmertmann(1978), Bustamante and Gianeselli(1982),Alsamman(1995),and Alkroosh and Nikraz(2011a,b)were evaluated.A brief description of the compared methods is provided in Table 2.The results ofcomparison are shown numerically in Table 3 and graphically in Figs.2 and 3. Table 3 also presents the results ofstatisticalanalysis used to evaluate the accuracy ofthe proposed method.Five statisticalindices including coef fi cient of determination,R,mean,μ,standard deviation,σ,mean absolute error,MAE,and mean squared error,MSE,were used in the analysis.The coef fi cient ofdetermination is calculated by

    where r is the coef fi cient of correlation;Q m i and Q p i are the measured and predicted pile bearing capacity of case i,respec-are the averages of measured and predicted pile bearing capacity,respectively;n is the number of observations.

    The optimalvalue of R is unity,which means that a perfect fi t is achieved between predicted and measured values.Table 3 shows that the calculated R=1 for the training and validation sets and R=0.99 for the testing set indicate that the LSSVMcan achieve an accurate correlation between predicted and measured pile bearing capacity.Table 3 also shows that the proposed method is more accurate than other methods.

    Table 2 Some CPT methods versus the LSSVM.

    Table 3 Results of evaluating the performance of LSSVM.

    Then the mean and standard deviation were respectively calculated by

    If the mean value is unity,the predicted values of pile bearing capacity are on average equivalent to the measured values.A mean value of less than unity indicates that the method tends to underpredict the pile bearing capacity.Conversely,more than unity is an indication of over-prediction.The optimum value of standard deviation is zero;the closer the standard deviation approaches to zero,the greater the accuracy is.As presented in Table 3,the calculated mean values of the proposed method are 0.98,1.03 and 1.04 for the training,testing and validation sets,respectively,which indicate that the LSSVM possesses a high capability in predicting pile bearing capacity.The results also indicate that,on average,the proposed method may tend to over-predict the pile bearing capacity in testing and validation sets,whereas under-predict it in training set.In comparison with other methods,the calculated mean values indicate that the LSSVM provides better estimate for pile bearing capacity.

    Fig.2.Performance of the LSSVM model in training,testing and validation sets.

    The predictive ability of the LSSVM is also evaluated by calculating the error.Table 3 shows that the calculated MSE and MAE are low when using the proposed method to predict pile bearing capacity.The table also shows that the proposed method has the lowest error in comparison with other methods.

    As shown in Fig.2,the strong ability of LSSVM to predict the pile bearing capacity is obvious.All the points in training,testing and validations sets are situated on or very close to the line of equality,suggesting that the proposed method is accurate in predicting the pile bearing capacity.Fig.3 shows that the LSSVM achieves very high correlation between predicted and measured values and it has the lowest scatter around the line of equality in comparison with other methods.From Fig.3,it is most evident that the LSSVM provides best prediction of total pile bearing capacity.

    In order to verify whether or not the results of the proposed method agree with those in the existing literature and experimental results,the sensitivity analysis was carried out.In this analysis,values of one input variable were allowed to change within the range of the training data whereas the values of other input variables were set constant.The data were then input in the developed LSSVM and the outputs were evaluated.The results,as presented in Fig.4,have shown that the proportion between each input variable and the predicted pile bearing capacity is incremental,and the variations of pile diameter and average cone point resistance within pile tip infl uence zone and along pile shaft have great in fl uence on pile bearing capacity.The variation of pile length slightly affects pile bearing capacity.

    5.Conclusions

    Fig.3.Comparisons ofthe performance of the LSSVM modelwith CPT-based methods.

    The results of this study indicate that the proposed LSSVMcan accurately predict the bearing capacity of bored piles from CPT data.The statistical analysis of the results reveals that this technique can obtain coef fi cient of determination and mean close to unity for all training,testing and validation sets.It has also achieved standard deviation less than 0.2 and low error values. The comparisons of the predictive ability of LSSVM with traditional CPT-based methods as well as the GEP model have shown that this technique predicts pile bearing capacity better.The relationship of pile bearing capacity with each of pile diameter, length,average of cone point resistance is incremental.The average cone point resistance within pile tip zone is the most in fl uential factor on pile bearing capacity.The output of this study suggests that LSSVM can be used as a design tool forpredicting the bearing capacity of bored piles installed in sand and mixed soils.

    Fig.4.Effects of variations of input variables on pile bearing capacity:(a)pile diameter;(b)pile length;(c)average cone point resistance along shaft;(d)average cone point resistance within tip zone.

    Con fl ict of interest

    The authors con fi rm that there are no known con fl icts of interest associated with this publication and there has been no signi fi cant fi nancialsupportforthis work thatcould have in fl uencedits outcome.

    Abu-Farsakh M,Titi H.Assessment of direct cone penetration test methods for predicting the ultimate capacity offriction driven piles.Journalof Geotechnical and Environmental Engineering 2004;130(9):935-44.

    Abu-Kiefa MA.General regression neural networks for driven piles in cohesionless soils.Journalof Geotechnicaland GeoenvironmentalEngineering 1998;124(12): 1177-85.

    Alavi AH,Aminian P,Gandomi AH,Esmaeili MA.Genetic based modeling of uplift capacity of suction caissons.Expert Systems with Applications 2011;38(10): 12608-18.

    Alkroosh I,Nikraz H.Correlation of pile axial capacity and CPT data using gene expression programming.Geotechnical and Geological Journal 2011a;29(5): 725-48.

    Alkroosh I,Nikraz H.Predicting axialcapacity ofdriven piles in cohesive soils using intelligent computing.Engineering Applications of Arti fi cial Intelligence 2011b;25(3):618-27.

    Alsamman O.The use of the CPT for calculating axial capacity of drilled shafts.PhD Thesis.Urbana,USA:University of Illinois at Urbana;1995.

    Amendolia SR,Cossu G,Ganadu ML,Golosio B,Masala GL,Mura GM.A comparative study of K-nearest neighbour,support vector machine and multi-layer perceptron for Thalassemia screening.Chemometrics and Intelligent Laboratory Systems 2003;69(1-2):13-20.

    Ardalan H,EslamiA,Nariman-Zadeh N.Piles shaftcapacity from CPT and CPTU data by polynomial neural networks and genetic algorithms.Computers and Geotechnics 2009;36(4):616-25.

    Baylar A,Hanbay D,Batan M.Application ofleastsquare support vector machines in the prediction of aeration performance of plunging overfall jets from weirs. Expert Systems with Applications 2009;36(4):8368-74.

    Bazzani A,Bevilacqua A,Bollini D,Brancaccio R,Campanini R,Lanconelli N, Riccardi A,Romani D.An SVM classi fi er to separate false signals from microcalci fi cations in digital mammograms.Physics in Medicine and Biology 2001;46(6):1651.

    Briaud J.Evaluation of cone penetration test methods using 98 pile load tests.In: Proceedings of the 1st International Symposium on Penetration TestingVol.2. Rotterdam,Netherlands:A.A.Balkema;1988.p.687-97.

    Bustamante M,Gianeselli L.Pile bearing capacity prediction by means of static penetrometer CPT.In:Proceedings of the 2nd European Symposium on Penetration Testing,ESOPT-II,Amsterdam,The Netherland;1982.p.493-500.

    Cai G,Liu S,Tong L,Du G.Assessment of direct CPT and CPTU methods for predicting the ultimate bearing capacity of single piles.Engineering Geology 2008;104(3-4):211-22.

    Chamkalani A,Amani M,Kiani MA,Chamkalani R.Assessment of asphaltene deposition due to titration technique.Fluid Phase Equilibria 2013;339:72-80.

    Chen TS,Chen J,Lin YC,Tsai YC,Kao YH,Wu KS.A novel knowledge protection technique based on support vector machine model for anti-classifi cation.In:Electricalengineering and control.Berlin,Germany:Springer;2011. p.517-24.

    Cortes C,Vapnik V.Support-vector networks.Machine Learning 1995;20(3):273-97.

    Das S,Samui P,Kim D,Sivakugan N,Biswal R.Lateral displacement of liquefaction induced ground using least square support vector machine.International Journal of Geotechnical Earthquake Engineering 2011b;2(2):29-39.

    Das S,Samui P,Kothari D.Site characterization model using machine learning.In: Machine tools:design,reliability and safety.New York,USA:Nova Science Publishers Inc.;2011a.p.175-85.

    Das SK,Basudhar PK.Undrained lateral load capacity of piles in clay using arti fi cial neural network.Computers and Geotechnics 2006;33(8):454-9.

    Das SK.Arti fi cial neural networks in geotechnical engineering:modeling and application issues.In:Yan X,Gandomi AH,Talatahati S,Alavi AH,editors. Metaheuristics in water;2013.p.231-70.

    Deng S,Yeh TH.Applying least squares support vector machines to the airframe wing-box structural design cost estimation.Expert Systems with Applications 2010;37(12):8417-23.

    Eslami A.Bearing capacity of piles from cone penetration data.PhD Thesis.Ottawa, Canada:University of Ottawa;1997.

    Eslamimanesh A,Gharagheizi F,Mohammadi AH,Richon D.Phase equilibrium modeling of structure H clathrate hydrates of methane+water“insoluble”hydrocarbon promoter using group contribution-support vector machine technique.Industrial and Engineering Chemistry Research 2011;50(20): 12807-14.

    Gandomi AH,Alavi AH.A new multi-gene genetic programming approach to nonlinear system modeling.Part II:geotechnical and earthquake engineering problems.Neural Computing and Applications 2012;21(1):189-201.

    Gandomi AH.Nonlinear genetic based models for prediction of fl ow number of asphalt mixtures.Journal of Materials in Civil Engineering 2011;23(3):248-63. Gharagheizi F,Eslamimanesh A,Farjood F,Mohammadi AH,Richon D.Solubility parameters of nonelectrolyte organic compounds:determination using quantitative structure-property relationship strategy.Industrial and Engineering Chemistry Research 2011;50(19):11382-95.

    Gunn SR.Support vector machines for classi fi cation and regression.Technical report.Southampton,UK:Faculty of Engineering,Science and Mathematics, Schoolof Electronics and Computer Science,University of Southampton;1998. Kulhawy F,Mayne P.Manual on estimating soil properties for foundation design. Palo Alto,CA,USA:Electric Power Research Institute;1990.

    Minoux M.Mathematical programming:theory and algorithms.New York,USA: Wiley;1986.

    Müller KR,Mika S,R?tsch G,Tsuda K,Sch?lkopf B.An introduction to kernel-based learning algorithms.IEEE Transactions on Neural Networks 2001;12(2):181-201.

    Ng C,Simons W,Menzies B.Soil-structure engineering of deep foundations,excavations and tunnels.London,UK:Thomas Telford Ltd.;2004.

    Ornek M,Laman M,Demir A,Yildiz A.Prediction of bearing capacity of circular footings on soft clay stabilized with granular soil.Soils and Foundations 2012;52(1):69-80.

    Ra fi ee-Taghanaki S,Arabloo M,ChamkalaniA,AmaniM,ZargariMH,Adelzadeh MR. Implementation of SVM framework to estimate PVT properties of reservoir oil. Fluid Phase Equilibria 2013;346:25-32.

    Rezania M,Javadi A.A new genetic programming model for predicting settlement of shallow foundations.Canadian Geotechnical Journal 2007;44(12):1462-73. Roberston P,Campanella R,Davies M,Sy A.Axialcapacity ofdriven piles in deltaic soils using CPT.In:Proceedings of the 1st International Symposium on Penetration TestingVol.2.Rotterdam,Netherlands:A.A.Balkema;1988.p.919-28. Samui P,Kothari DP.Utilization of a least square support vector machine(LSSVM) for slop stability analysis.Scientia Iranica 2011;18(1):53-8.

    Schmertmann JH.Guidelines for cone penetration test,performance and design. Washington D.C.,USA:US Department of Transportation;1978.

    Shahin M.Intelligent computing for modeling axial capacity of pile foundations. Canadian Geotechnical Journal 2010;47(2):230-43.

    Shokrollahi A,Arabloo M,Gharagheizi F,Mohammadi AH.Intelligent model for prediction of CO2-reservoir oil minimum miscibility pressure.Fuel 2013;112: 375-84.

    Suykens JA,Van Gestel T,De Brabanter J,De Moor B,Vandewalle J.Least squares support vector machines.Singapore:World Scienti fi c Publishing Company; 2002.

    Suykens JA,Vandewalle J.Least squares support vector machine classi fi ers.Neural Processing Letters 1999;9(3):293-300.

    Tarawneh B.Pipe pile setup:database and prediction modelusing arti fi cialneural network.Soils and Foundations 2013;53(4):607-15.

    Teh CI,Wong KS,Goh AT,Jaritngam S.Prediction of pile capacity using neural networks.Journal of Computing in Civil Engineering 1997;11(2):129-38.

    übeyli ED.Least squares support vector machine employing model-based methods coef fi cients for analysis of EEG signals.Expert Systems with Applications 2010;37(1):233-9.

    Xavier-de-Souza S,Suykens JA,Vandewalle J,Bolle D.Coupled simulated annealing. IEEE Transactions on Systems,Man,and Cybernetics.Part B:Cybernetics 2009;40(2):320-35.

    Dr.Iyad Alkrooshbegan his career as a civilengineer after graduation from University of Technology Baghdad in 1987.He received his MSc and PhD degrees in civilengineering from Curtin University,Western Australia in 2000 and 2011,respectively.Currently he has been a lecturer at Department of Civil Engineering,College of Engineering, Al-Qadisiya University,Iraq.Since 2011,Dr.Alkroosh has been working as a consultant for the Scienti fi c and Engineering Bureau,Al-Qadisiya University.His fi eld ofinterest includes but not limited to modeling using arti fi cial intelligence,load transfer behavior of piles,soil improvement, and structural design of shallow and deep foundations. Dr.Alkroosh is the author of 15 publications in international journals and conferences.

    Mohammad Bahadoriis a research student in the Department of Soil Science and Engineering at the University of Tehran,Iran.He received his BSc degree from Shahaid Chamran University,Ahwaz,Iran and his MSc degree from the University of Tehran,Iran.

    Dr.Hamid Nikrazis a professor at the Department of Civil Engineering,School of Engineering,Curtin University, Australia.Prior to his academic career he had many years’experience working in industry in the Middle East,Europe and Australia.He has provided numerous consultancy services in his fi eld to over 200 organizations nationally and internationally.He obtained his PhD at Curtin University of Technology in 1989.Dr.Nikraz is a Fellow Member of the Institution of Engineers,Australia.He has guided 75 PhD theses,30 Master students and published over 500 publications in various journals,books and conference proceedings.Prof.Nikraz has initiated research works in the area of Geopolymers,Concrete,Soil Stabilisation and Pavement Engineering.Prof.Nikraz is also the reviewer of many national and international journals.

    Alireza Bahadori,PhD,P.Eng,MIEAust.,is a research staff member and lecturer in the School of Environment,Science and Engineering at Southern Cross University,Lismore,NSW,Australia.He received his PhD from Curtin University,Perth,Western Australia.During the past twenty years,Dr.Bahadori has held various process and petroleum engineering positions and was involved in many large-scale projects at National Iranian Oil Co. (NIOC),Petroleum Development Oman(PDO),and Clough AMEC Pty Ltd.He is the author of several articles and books.His books have been published by multiple major publishers,including Elsevier,Springer,John Wiley& Sons,and Taylor and Francis Group.Dr.Bahadori is the recipientofthe highly competitive and prestigious Australian Government’s Endeavour International Postgraduate Research Award as part of his research in oiland gas area.He also received a Top-Up Award from the State Government of Western Australia through Western Australia Energy Research Alliance in 2009.Dr.Bahadori serves as a member of the editorial board and reviewer for a large number ofjournals.He is member ofthe Institution of Engineers Australia as a professional engineer and IChemE,UK.

    *Corresponding author.Tel.:+61 2 6626 9347.

    E-mail address:Alireza.bahadori@scu.edu.au(A.Bahadori).

    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.06.011

    国产精品国产三级专区第一集| 久久久久网色| 国产男女内射视频| 精品一区在线观看国产| 毛片一级片免费看久久久久| 色视频在线一区二区三区| 午夜久久久在线观看| 日本wwww免费看| 美女主播在线视频| 亚洲激情五月婷婷啪啪| 高清在线视频一区二区三区| 香蕉精品网在线| 国产精品一区二区在线不卡| 日日撸夜夜添| 2018国产大陆天天弄谢| 高清av免费在线| 黑人欧美特级aaaaaa片| 人妻人人澡人人爽人人| 在现免费观看毛片| 国产乱来视频区| 妹子高潮喷水视频| 韩国高清视频一区二区三区| 晚上一个人看的免费电影| 国产一级毛片在线| 国产精品一二三区在线看| 国产精品一二三区在线看| 久热久热在线精品观看| 日韩制服丝袜自拍偷拍| 乱码一卡2卡4卡精品| 18+在线观看网站| 久久精品国产亚洲av涩爱| 成年人免费黄色播放视频| 波野结衣二区三区在线| 欧美最新免费一区二区三区| www.av在线官网国产| 国产免费福利视频在线观看| 22中文网久久字幕| 街头女战士在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 国产成人91sexporn| 国产亚洲午夜精品一区二区久久| 欧美精品国产亚洲| www.熟女人妻精品国产 | 激情视频va一区二区三区| 在线免费观看不下载黄p国产| 黑人巨大精品欧美一区二区蜜桃 | 2018国产大陆天天弄谢| 久久精品久久精品一区二区三区| 高清不卡的av网站| 又黄又粗又硬又大视频| 黄色 视频免费看| 日本91视频免费播放| 久久精品aⅴ一区二区三区四区 | 丝袜在线中文字幕| 国产一级毛片在线| 久热久热在线精品观看| 亚洲图色成人| 久久久精品免费免费高清| 日韩电影二区| 久久狼人影院| 欧美激情极品国产一区二区三区 | av不卡在线播放| av女优亚洲男人天堂| 天堂俺去俺来也www色官网| 精品国产一区二区三区四区第35| 国产av精品麻豆| 国产精品国产三级专区第一集| 亚洲av成人精品一二三区| 99久久中文字幕三级久久日本| 中文字幕最新亚洲高清| 99九九在线精品视频| 国产亚洲一区二区精品| 巨乳人妻的诱惑在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品99久久99久久久不卡 | 黑人高潮一二区| 国产亚洲欧美精品永久| 久久久久久久精品精品| 在现免费观看毛片| 亚洲图色成人| 欧美激情极品国产一区二区三区 | 国产综合精华液| 丰满乱子伦码专区| 午夜激情久久久久久久| 久久久久久伊人网av| 亚洲综合色惰| 亚洲欧洲日产国产| 日本av手机在线免费观看| 久久97久久精品| 观看美女的网站| 久久久久人妻精品一区果冻| 91精品国产国语对白视频| www日本在线高清视频| 9热在线视频观看99| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 国产成人精品久久久久久| 极品人妻少妇av视频| 国产精品一区二区在线不卡| 99久国产av精品国产电影| 成人二区视频| 巨乳人妻的诱惑在线观看| 亚洲精华国产精华液的使用体验| 90打野战视频偷拍视频| 精品福利永久在线观看| 国产一区二区在线观看日韩| 男女边吃奶边做爰视频| 99久国产av精品国产电影| av线在线观看网站| 蜜桃在线观看..| 男女下面插进去视频免费观看 | 丝袜在线中文字幕| 成年动漫av网址| 欧美最新免费一区二区三区| 国产精品国产三级专区第一集| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜美足系列| 激情视频va一区二区三区| 久久精品aⅴ一区二区三区四区 | 亚洲国产精品999| 亚洲精品自拍成人| 亚洲av男天堂| 咕卡用的链子| 高清毛片免费看| 中国美白少妇内射xxxbb| 在线看a的网站| 2022亚洲国产成人精品| 女人精品久久久久毛片| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 精品亚洲成国产av| 欧美+日韩+精品| 国产永久视频网站| 日本av手机在线免费观看| 亚洲av在线观看美女高潮| 欧美日韩一区二区视频在线观看视频在线| 国产精品人妻久久久影院| 久久精品夜色国产| 久久女婷五月综合色啪小说| 国产亚洲午夜精品一区二区久久| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| av在线播放精品| videos熟女内射| 亚洲精品aⅴ在线观看| 欧美精品av麻豆av| av在线app专区| 男女下面插进去视频免费观看 | 国产精品一区二区在线观看99| 国产精品国产三级国产av玫瑰| 国产精品一区www在线观看| 免费在线观看完整版高清| 久久精品国产自在天天线| 午夜福利影视在线免费观看| 人人妻人人爽人人添夜夜欢视频| xxx大片免费视频| 少妇人妻精品综合一区二区| 人人妻人人澡人人爽人人夜夜| 亚洲国产欧美在线一区| 晚上一个人看的免费电影| 国产精品久久久久久久电影| 国产精品女同一区二区软件| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| 97在线视频观看| 美国免费a级毛片| 国产69精品久久久久777片| 日韩精品免费视频一区二区三区 | 99香蕉大伊视频| 大片电影免费在线观看免费| 啦啦啦中文免费视频观看日本| 国产乱来视频区| 这个男人来自地球电影免费观看 | 精品亚洲成a人片在线观看| av视频免费观看在线观看| 免费观看在线日韩| 如日韩欧美国产精品一区二区三区| 欧美日韩视频精品一区| 夫妻午夜视频| 春色校园在线视频观看| 午夜精品国产一区二区电影| av.在线天堂| 国产精品一区二区在线不卡| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| av不卡在线播放| 亚洲,欧美,日韩| 色婷婷av一区二区三区视频| 亚洲av综合色区一区| 在线天堂最新版资源| 男女国产视频网站| av在线老鸭窝| 最新中文字幕久久久久| 99热这里只有是精品在线观看| 久久国产精品大桥未久av| 午夜久久久在线观看| 欧美成人精品欧美一级黄| 中文字幕人妻丝袜制服| 91在线精品国自产拍蜜月| tube8黄色片| av在线观看视频网站免费| 美女内射精品一级片tv| 深夜精品福利| 一边亲一边摸免费视频| 男人舔女人的私密视频| 大香蕉久久成人网| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 欧美少妇被猛烈插入视频| 晚上一个人看的免费电影| 久久99一区二区三区| 综合色丁香网| 久久精品国产亚洲av天美| 国产免费现黄频在线看| 色婷婷av一区二区三区视频| 咕卡用的链子| 国产深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 性高湖久久久久久久久免费观看| 精品一区二区免费观看| 亚洲高清免费不卡视频| 男人舔女人的私密视频| 亚洲经典国产精华液单| 91午夜精品亚洲一区二区三区| 色哟哟·www| 久久久久久人妻| 美女大奶头黄色视频| 亚洲精品av麻豆狂野| 久久狼人影院| 一本色道久久久久久精品综合| 在线观看www视频免费| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 久久婷婷青草| 亚洲综合色惰| 男人舔女人的私密视频| 制服诱惑二区| 在线 av 中文字幕| 人体艺术视频欧美日本| 在线观看国产h片| 国产亚洲精品第一综合不卡 | 亚洲国产欧美日韩在线播放| 成人毛片a级毛片在线播放| 在线观看免费高清a一片| 一区二区日韩欧美中文字幕 | 热re99久久精品国产66热6| 亚洲第一区二区三区不卡| 综合色丁香网| 欧美老熟妇乱子伦牲交| 97在线视频观看| videos熟女内射| 国产 一区精品| 国产在视频线精品| 亚洲丝袜综合中文字幕| 黄色怎么调成土黄色| 有码 亚洲区| 亚洲欧美日韩卡通动漫| 国产亚洲精品久久久com| 丝袜在线中文字幕| 亚洲图色成人| 亚洲第一区二区三区不卡| 精品午夜福利在线看| 久久久久久人人人人人| 成人手机av| 在线 av 中文字幕| 一级毛片黄色毛片免费观看视频| 亚洲成人手机| av有码第一页| 亚洲 欧美一区二区三区| 国产av一区二区精品久久| 18禁裸乳无遮挡动漫免费视频| 国产成人欧美| 只有这里有精品99| 秋霞在线观看毛片| 一级,二级,三级黄色视频| 最近手机中文字幕大全| 黄色一级大片看看| 久久99热这里只频精品6学生| 久热这里只有精品99| 桃花免费在线播放| 欧美成人精品欧美一级黄| 婷婷色麻豆天堂久久| 色婷婷久久久亚洲欧美| av.在线天堂| av卡一久久| 亚洲,欧美精品.| 一边亲一边摸免费视频| 国产免费又黄又爽又色| 99久久中文字幕三级久久日本| 三上悠亚av全集在线观看| 蜜臀久久99精品久久宅男| 国产精品嫩草影院av在线观看| 国产白丝娇喘喷水9色精品| 亚洲情色 制服丝袜| 秋霞伦理黄片| 有码 亚洲区| 国产高清三级在线| 18禁裸乳无遮挡动漫免费视频| 狠狠精品人妻久久久久久综合| 国产精品国产三级专区第一集| 男人爽女人下面视频在线观看| 大陆偷拍与自拍| av在线播放精品| 久久av网站| 97在线视频观看| 亚洲第一av免费看| 国产欧美另类精品又又久久亚洲欧美| tube8黄色片| 极品人妻少妇av视频| 久久久久久久大尺度免费视频| freevideosex欧美| 亚洲精品成人av观看孕妇| 满18在线观看网站| av不卡在线播放| 国产高清不卡午夜福利| av卡一久久| 在线天堂中文资源库| 9色porny在线观看| 国产成人午夜福利电影在线观看| 欧美精品av麻豆av| 宅男免费午夜| 久久热在线av| a级片在线免费高清观看视频| 国产成人精品婷婷| 久久精品熟女亚洲av麻豆精品| 日韩大片免费观看网站| 久久久国产欧美日韩av| 中国三级夫妇交换| 国产色爽女视频免费观看| 免费久久久久久久精品成人欧美视频 | 久久鲁丝午夜福利片| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久人人人人人| 日韩伦理黄色片| 婷婷色av中文字幕| 91aial.com中文字幕在线观看| 一二三四在线观看免费中文在 | 中文字幕av电影在线播放| 在线精品无人区一区二区三| 在线观看美女被高潮喷水网站| a级片在线免费高清观看视频| www.熟女人妻精品国产 | 中文字幕av电影在线播放| 99国产精品免费福利视频| 亚洲国产精品一区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久精品94久久精品| av视频免费观看在线观看| 亚洲综合色网址| 国产一区有黄有色的免费视频| 99久久中文字幕三级久久日本| 国产在线一区二区三区精| 成人国产麻豆网| 女人被躁到高潮嗷嗷叫费观| 国产国语露脸激情在线看| 久久久久久人人人人人| 少妇人妻 视频| 18禁动态无遮挡网站| 哪个播放器可以免费观看大片| 亚洲av综合色区一区| 欧美精品一区二区免费开放| 如何舔出高潮| 蜜臀久久99精品久久宅男| 韩国av在线不卡| 婷婷色麻豆天堂久久| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 麻豆精品久久久久久蜜桃| 婷婷色麻豆天堂久久| 99热6这里只有精品| 亚洲精品视频女| 国产精品99久久99久久久不卡 | 免费观看在线日韩| 亚洲国产色片| freevideosex欧美| 婷婷色av中文字幕| 国产国语露脸激情在线看| 午夜精品国产一区二区电影| 久久这里只有精品19| 亚洲国产精品国产精品| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 国产男人的电影天堂91| 久久这里有精品视频免费| 欧美另类一区| 国产伦理片在线播放av一区| 啦啦啦中文免费视频观看日本| 免费大片黄手机在线观看| 高清不卡的av网站| 另类亚洲欧美激情| 中国三级夫妇交换| 久久精品国产综合久久久 | 欧美成人午夜免费资源| 熟女电影av网| 亚洲av电影在线观看一区二区三区| 内地一区二区视频在线| 日韩av不卡免费在线播放| 观看美女的网站| 青春草视频在线免费观看| 久久久久精品久久久久真实原创| 曰老女人黄片| 看十八女毛片水多多多| 欧美精品av麻豆av| 国产成人一区二区在线| 免费观看无遮挡的男女| 熟女电影av网| 欧美日韩av久久| 婷婷色麻豆天堂久久| 欧美丝袜亚洲另类| 国产国语露脸激情在线看| 亚洲av在线观看美女高潮| 一级毛片 在线播放| 国产成人精品一,二区| 一边亲一边摸免费视频| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 老司机影院毛片| 丰满少妇做爰视频| 天天操日日干夜夜撸| 99热6这里只有精品| 亚洲美女黄色视频免费看| 久久久久久人人人人人| 七月丁香在线播放| 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久 | 欧美精品高潮呻吟av久久| 日韩,欧美,国产一区二区三区| 日本wwww免费看| 在线免费观看不下载黄p国产| 国产毛片在线视频| 精品一区二区免费观看| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 又粗又硬又长又爽又黄的视频| 自线自在国产av| av片东京热男人的天堂| 久久精品熟女亚洲av麻豆精品| 91久久精品国产一区二区三区| 亚洲,欧美,日韩| 美女福利国产在线| 亚洲国产精品一区二区三区在线| 亚洲精品第二区| 国产 一区精品| 99久久人妻综合| 国产精品人妻久久久影院| 考比视频在线观看| 天美传媒精品一区二区| 精品国产一区二区三区久久久樱花| 高清欧美精品videossex| 母亲3免费完整高清在线观看 | 最近2019中文字幕mv第一页| 日本色播在线视频| 精品亚洲成国产av| av一本久久久久| 免费看光身美女| 99久久综合免费| 免费女性裸体啪啪无遮挡网站| 国产黄色视频一区二区在线观看| 中文欧美无线码| 久久人人97超碰香蕉20202| 丝袜美足系列| 国产av精品麻豆| 9热在线视频观看99| 免费人妻精品一区二区三区视频| 纯流量卡能插随身wifi吗| 亚洲综合色惰| 考比视频在线观看| 免费观看无遮挡的男女| 国产成人精品久久久久久| 一级黄片播放器| 少妇被粗大猛烈的视频| 又黄又爽又刺激的免费视频.| 男女高潮啪啪啪动态图| www.av在线官网国产| 夜夜骑夜夜射夜夜干| 一级,二级,三级黄色视频| 国产一级毛片在线| 各种免费的搞黄视频| 国产又色又爽无遮挡免| 在线观看免费高清a一片| 欧美xxⅹ黑人| 热re99久久精品国产66热6| 精品一区二区三区四区五区乱码 | 午夜福利,免费看| 毛片一级片免费看久久久久| 人成视频在线观看免费观看| 999精品在线视频| 日日撸夜夜添| 蜜桃在线观看..| 午夜福利视频精品| 久久99一区二区三区| 久久久精品区二区三区| 七月丁香在线播放| 蜜桃在线观看..| 男男h啪啪无遮挡| 国产成人精品久久久久久| 黄网站色视频无遮挡免费观看| 日韩不卡一区二区三区视频在线| 成人国产av品久久久| 少妇熟女欧美另类| 又大又黄又爽视频免费| 男女无遮挡免费网站观看| 观看美女的网站| 国产有黄有色有爽视频| 韩国av在线不卡| 精品一区二区三区四区五区乱码 | 免费不卡的大黄色大毛片视频在线观看| 午夜视频国产福利| 久久 成人 亚洲| 国产免费一级a男人的天堂| 超碰97精品在线观看| 成人国产av品久久久| 高清毛片免费看| 国产激情久久老熟女| 亚洲在久久综合| 中文字幕精品免费在线观看视频 | 成年av动漫网址| 这个男人来自地球电影免费观看 | 欧美日韩av久久| 五月伊人婷婷丁香| 丝袜喷水一区| 91久久精品国产一区二区三区| 伦理电影大哥的女人| 亚洲精品日本国产第一区| 日本av免费视频播放| av在线老鸭窝| 视频中文字幕在线观看| 狠狠婷婷综合久久久久久88av| 亚洲五月色婷婷综合| 精品久久久久久电影网| 精品人妻在线不人妻| 免费观看无遮挡的男女| 伊人亚洲综合成人网| 国产精品蜜桃在线观看| 国产一区二区在线观看日韩| 高清毛片免费看| 免费av中文字幕在线| 日韩熟女老妇一区二区性免费视频| 丝袜脚勾引网站| 久久精品久久久久久噜噜老黄| 在线观看免费视频网站a站| 黄色毛片三级朝国网站| 精品人妻偷拍中文字幕| 老熟女久久久| 久久久久人妻精品一区果冻| 国产熟女欧美一区二区| 啦啦啦啦在线视频资源| 亚洲伊人久久精品综合| 中国美白少妇内射xxxbb| 亚洲综合精品二区| 日本av免费视频播放| 免费少妇av软件| av在线播放精品| av网站免费在线观看视频| 久久久久久久大尺度免费视频| 日韩av免费高清视频| 女性生殖器流出的白浆| 在线免费观看不下载黄p国产| 亚洲国产精品专区欧美| 色5月婷婷丁香| 国内精品宾馆在线| 看非洲黑人一级黄片| 80岁老熟妇乱子伦牲交| 99九九在线精品视频| 亚洲第一区二区三区不卡| a 毛片基地| 91国产中文字幕| 下体分泌物呈黄色| 晚上一个人看的免费电影| 少妇高潮的动态图| 99久久综合免费| 最近中文字幕2019免费版| 99久国产av精品国产电影| 中文天堂在线官网| a级毛色黄片| 日韩中字成人| 丁香六月天网| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 欧美日韩视频高清一区二区三区二| 黄色配什么色好看| 丰满少妇做爰视频| 人成视频在线观看免费观看| 国产亚洲一区二区精品| 又大又黄又爽视频免费| 国产色爽女视频免费观看| 2018国产大陆天天弄谢| 日韩成人伦理影院| 久久国产亚洲av麻豆专区| 国产av码专区亚洲av| a 毛片基地| 91国产中文字幕| 亚洲精品国产av蜜桃| 草草在线视频免费看| 日韩视频在线欧美| 多毛熟女@视频| 成年人免费黄色播放视频| 一二三四中文在线观看免费高清| 90打野战视频偷拍视频| 亚洲精品日本国产第一区| 国产一区二区在线观看日韩| 另类亚洲欧美激情| 亚洲高清免费不卡视频| 人成视频在线观看免费观看| 人妻少妇偷人精品九色| av播播在线观看一区| 国产老妇伦熟女老妇高清| 日韩制服丝袜自拍偷拍| 美女国产视频在线观看| 久久婷婷青草| 成人18禁高潮啪啪吃奶动态图| 在线天堂中文资源库| 嫩草影院入口| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美中文字幕日韩二区|