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

    樽海鞘算法優(yōu)化支持向量機的RC柱抗側移承載力預測

    2024-05-29 02:06:26歐陽謙駱歡
    地震研究 2024年3期
    關鍵詞:支持向量機特征選擇

    歐陽謙 駱歡

    摘要:現有鋼筋混凝土(RC)柱抗側移承載力預測模型缺乏泛化性能,延性柱抗彎承載力的預測模型不能用于非延性柱的抗剪承載力,反之亦然。機器學習(ML)方法能夠解決這一問題,但由于無法自動剔除冗余和不相關特征,使得ML模型復雜度高且容易過擬合。為此,提出一種樽海鞘算法優(yōu)化支持向量機(SSALS-SVM)方法,基于給定的數據集,SSALS-SVM能利用樽海鞘優(yōu)化算法(SSA)自動剔除冗余和不相關的特征,篩選最具代表性且各特征之間相關性弱的特征子集形成最優(yōu)特征組合,同時對控制模型非線性擬合能力的超參數進行優(yōu)化。優(yōu)化后的模型既能識別出影響延性和非延性RC柱抗側移承載力的設計變量,又能反映最優(yōu)特征組合與抗側移承載力間的非線性映射關系。為了驗證SSALS-SVM方法的泛化性能,基于248個RC柱抗側移承載力試驗數據,分別與現有的RC柱抗側移承載力預測模型進行對比,結果表明,SSALS-SVM比現有預測模型的泛化性能最高提升了83%。

    關鍵詞:鋼筋混凝土柱;抗側移承載力;支持向量機;樽海鞘優(yōu)化算法;特征選擇

    中圖分類號:TU973.2文獻標識碼:A文章編號:1000-0666(2024)03-0350-09

    doi:10.20015/j.cnki.ISSN1000-0666.2024.0051

    0引言

    地震是危害性極高的自然災害,強烈的地震會使建筑物的抗側移承載力下降、軸向承載能力喪失、非彈性形變機制發(fā)生改變,最終導致建筑物倒塌。在鋼筋混凝土(RC)框架結構中,RC柱主要承受豎向和橫向荷載,是框架結構主要的抗側力構件,其抗側移承載力包括延性柱的抗彎承載力和非延性柱的抗剪承載力(Noroozieh,Mansouri,2019)。RC柱發(fā)生破壞時,其失效模式分為三種:彎曲破壞、剪切破壞和彎剪破壞(于曉輝等,2022)。因此,準確預測RC柱的抗側移承載力有助于優(yōu)化RC框架的抗震設計。

    對于RC柱抗側移承載力預測,國內外許多學者開展了大量的研究工作(薛亦聰等,2020;張勤等,2014;甘丹等,2018;Pan,Li,2012)。在非延性柱抗剪承載力的預測方面,Priestley等(1994)提出了非延性柱抗剪承載力預測模型,將所得結果與剪切破壞柱試驗結果進行了比較,證明了該方法能準確預測非延性柱的抗剪承載力;鄧明科等(2018)基于桁架-拱模型對混凝土加固RC柱進行受力分析,并推導出了抗剪承載力計算公式;Sezen和Moehle(2004)基于非延性柱的試驗數據提出了一種新的柱抗剪承載力預測模型,并與現有的柱抗剪承載力預測模型進行對比,結果表明,模型的預測精度有所提高。數據驅動模型也被廣泛應用于預測RC柱抗剪承載力(Luo,Paal,2021;Aval et al,2017;Ketabdari et al,2020)。Kakavand等(2021)總結了現有基于力學計算公式中不同設計變量之間相互組合構成的輸入變量,基于497個非延性RC柱試件構成的數據集,按照考慮和不考慮位移延性的情況分別對這些輸入變量進行線性和非線性回歸分析,構建了用于預測非延性RC柱抗剪承載力的數據驅動模型。在延性柱抗彎承載力的預測方面,常用的方法是矩形應力塊法(Lee,Son,2000;Bae,Bayrak,2003),該方法需要估計兩個參數α1和β1來構建矩形應力塊。關于這些系數,《結構混凝土建筑規(guī)范要求》(ACI 318-22)與Ozbakkaloglu和Saatcioglu(2004)的研究給出了計算方法。雖然上述預測方法已被證明具有較好的預測性能,但是它們都存在著一定的局限性,即預測延性柱抗彎承載力的模型不能用于預測非延性柱的抗剪承載力,反之亦然。

    機器學習(ML)方法如支持向量機等已被證實能夠很好地解決上述問題,但其無法自動識別和剔除冗余與不相關特征,不能對數據集進行降維,這會增加模型的復雜度,容易導致最終形成的預測模型出現過擬合,進而影響模型的泛化性能。為此,許多研究者在如何將特征選擇與ML方法相結合的問題上開展了大量的研究工作(Guyon,Elisseeff,2003;Huang,Wang,2006)。Lin等(2008)將粒子群優(yōu)化算法與支持向量機相結合,最終形成的分類模型預測精度優(yōu)于網格搜索等方法。雖然上述研究成功地將特征選擇與ML方法相結合,但它們均聚焦在分類問題方面,而在回歸問題方面涉及較少。且RC柱抗側移承載力預測模型的構建過程屬于回歸問題,導致已有方法無法直接運用。

    基于以上研究結果,本文將樽海鞘優(yōu)化算法(SSA)(Mirjalili et al,2017)與ML方法中的支持向量機(LS-SVM)(Suykens et al,2002)相結合,提出了樽海鞘算法優(yōu)化支持向量機(SSALS-SVM)方法,收集了248個RC柱試件的試驗數據,建立了一組包含延性和非延性RC柱抗側移承載力的試驗數據集?;诖藬祿謩e與已有的計算模型進行了對比。

    1數學模型構建

    支持向量機(SVM)(Cortes,Vapnik,1995)是一種常用的監(jiān)督學習方法,其數學模型是基于1范數的損失函數,這使得其訓練過程需要求解復雜的二次規(guī)劃問題。而LS-SVM在繼承SVM優(yōu)點的情況下,將二次損失函數引入SVM,并將不等式約束條件替換為等式約束條件,極大地提高了計算效率。給定一組數據集{(xi,yi)}Ni,i=1,2,……,N,其中xi∈Rn為解釋變量,yi∈R為響應變量,LS-SVM通過將目標函數最小化來找出解釋變量與響應變量之間存在的非線性映射關系。其數學模型為:

    [HZ(]min[DD(X]w,b,ek[DD)]JP(w,ek)=[SX(]1[]2[SX)][WTHX]w[WTBX]T[WTHX]w[WTBX]+[SX(]1[]2[SX)]γ∑[DD(]N[]k=1[DD)][WTHX]e[WTBX]2k

    subject to y([WTHX]x[WTBX]k)=[WTHX]w[WTBX]Tφ([WTHX]x[WTBX]k)+b+[WTHX]e[WTBX]k[HZ)][JY](1)

    式中:[WTHX]e[WTBX]=[e1,e2,…,eN]T∈RN為誤差變量;[WTHX]w[WTBX]=[w1,w2,…,wh]∈Rh和b∈RN通過使目標函數最小化得到;γ為正則化參數;高維向量構建的矩陣[WTHX]φ(x)[WTBX]=[φ([WTHX]x[WTBX]1),φ([WTHX]x[WTBX]2),……,φ([WTHX]x[WTBX]N)]T∈RN×h,其中φ(·):Rn→Rh表示從n維到具有h維的高維度希伯特空間的映射函數;響應或輸出變量[WTHX]y[WTBX]=[y1,y2,……,yN]T∈RN。用拉格朗日乘子法將方程(1)的求解問題轉換為對偶問題求解,并引進拉格朗日乘子,以構建拉格朗日函數:

    L([WTHX]w[WTBX],b,ek,αk)=JP(w,[WTHX]e[WTBX]k)-∑[DD(]N[]k=1[DD)]αk[JB<2{][WTHX]w[WTBX]T[WTHX]φ(xk[WTBX])+b+eky-yk[JB>2}][JY](2)

    式中:[WTHX]α[WTBX]=[α1,α2,……,αN]T∈RN為拉格朗日乘子。利用Karush-Kuhn-Tucker(KKT)條件的最優(yōu)性求解式(2),即分別求拉格朗日函數對[WTHX]w[WTBX],b,ek,αk的偏導數,可得:

    (3)

    特征向量之間的內積在計算時較為復雜。為此,本文選用高斯核函數代替特征向量之間的內積,其計算模型如下:

    [WTHX]K[WTBX](xi,xj)=[WTHX]φ[WTBX]T([WTHX]x[WTBX]i)[WTHX]φ(x[WTBX]j)=exp[JB((]-[SX(][JB(=][WTHX]x[WTBX]i-[WTHX]x[WTBX]j[JB)=][]2σ2[SX)][JB))][JY](4)

    式中:σ2為核參數。LS-SVM模型的預測精度與兩個超參數(σ2和γ)的取值息息相關,因此,本文引入樽海鞘優(yōu)化算法(SSA)對其進行確定。此外,為了識別數據集中最能影響模型預測性能的解釋變量組合(即最優(yōu)特征組合),本文構建了一個控制解釋變量組合的參數f,并將其作為優(yōu)化參數耦合進SSA中,與前述兩個超參數同時進行優(yōu)化,進而形成樽海鞘算法優(yōu)化支持向量機(SSALS-SVM)的數學模型。具體過程如下所述:

    首先,給定f、σ2和γ 的取值范圍,隨機生成n個(f,σ2,γ)組合,其數學模型為:

    (fi,σ2i,γi)=rand*([WTHX]ub-lb)+lb[WTBX],i∈(1,2,…,n)[JY](5)

    式中:f為挑選的特征數目;[WTHX]ub[WTBX]為f、σ2和γ的上界;[WTHX]lb[WTBX]為f、σ2和γ的下界;d為參數的數目。根據隨機生成的f數值,隨機從數據集中選取f個解釋變量進行組合,作為最優(yōu)特征組合的候選者。此外,σ2和γ也需要同時進行優(yōu)化,因此需要建立f +2個決策變量,具體形式如圖1所示。

    如圖1所示,Pf為隨機挑選的解釋變量。將σ2、γ 以及Pf輸入目標函數中進行評估。對于目標函數,其建立過程如下所示:

    聯立式(3)中的各方程式,消除方程式中的[WTHX]w和e[WTBX]k,化簡可得:

    [JB([][HL(2]0[][WTHX]E[WTBX][WTHX]E[WTBX]T[][WTHX]K+I[WTBX]/γ[HL)][JB)]][JB([]b[WTHX]α[WTBX][JB)]]=[JB([]0[WTHX]y[WTBX][JB)]][JY](6)

    式中:[WTHX]E[WTBX]=[1,1,…,1]1×N;[WTHX]I[WTBX]為單位矩陣;[WTHX]K[WTBX]為核函數組成的核矩陣。

    通過求解矩陣方程(6),可得模型參數α和b,進而可建立對偶空間下的預測模型,其表達式為:

    為了判斷上述建立的f個解釋變量是否為最優(yōu)特征組合以及生成的σ2和γ值是否為最優(yōu)參數,需要使用一種評估函數對通過方程(7)得到的結果進行評估。為此,本文采用均方誤差(MSE)作為模型性能優(yōu)劣的評估標準。其數學模型為:

    MSE=[SX(]∑[DD(]n[]i=1[DD)](yi-Ay^G1i)2[]n[SX)][JY](8)

    因此,目標函數由方程(7)和方程(8)組成:方程(7)對結果進行預測,方程(8)對預測結果進行評估?;谝呀⒌哪繕撕瘮?,對n個(f,σ2,γ)組合進行評估,并對比n個MSE的大小,取其最小值作為本次評估的最佳適應度,對應的(f,σ2,γ)組合作為當前最優(yōu)位置F。

    將n個(f,σ2,γ)組合前一半設定為領導者,后一半設定為跟隨者。領導者的位置跟F相關,根據F,通過一次次位置更新,不斷向其靠近,位置更新公式為:

    (fi,σ2i,γi)=[JB({]F+ri[([WTHX]ub-lb[WTBX])×r2+[WTHX]lb[WTBX]] if r3≤0.5

    F-r1[([WTHX]ub-lb[WTBX])×r2+[WTHX]lb[WTBX]] if r3>0.5[JB)][JY](9)

    式中:r1=2e-([SX(]4t[]T[SX)])[HT4.5]2;t表示當前迭代次數;T表示最大迭代次數;r3和r2是一個0~1的隨機數,r3大于或小于0.5將決定領導者前進還是后退。跟隨者的位置則與前一個個體位置相關,其數學模型為:[KH*1D]

    (fi,σ2i,γi)=[SX(](fi,σ2i,γi)+(fi-1,σ2i-1,γi-1)[]2[SX)][JY](10)

    對于更新后的領導者和跟隨者,其位置可能會超出設定的范圍,因此需要將其重新限制到規(guī)定范圍內。此外,當領導者和跟隨者位置迭代更新后,需重新對其進行評估:將更新后的(f,σ2,γ)組合重新代入方程(4)、(6)和(7)中計算其預測值Ay^G1,并通過方程(8)計算其MSE。將更新后的每個(f,σ2,γ)組合的MSE與當前最佳適應度進行比較,若更新后的(f,σ2,γ)適應度優(yōu)于當前最佳適應度,則以MSE更優(yōu)的(f,σ2,γ)組合作為F。

    重復以上步驟,直到達到最大迭代次數或適應度達到最佳(MSE=0),輸出當前的F作為最優(yōu)(f,σ2,γ)組合,并輸出最優(yōu)特征組合(P1,…,Pf),進而可得到LS-SVM預測模型的最優(yōu)形式。因此,從理論分析上看,本文提出的SSALS-SVM方法既能夠優(yōu)化模型的超參數σ2和γ,又能自動篩選數據中最能影響模型預測性能的解釋變量x,形成最優(yōu)特征組合(P1,…,Pf)。

    給定一組數據集,圖2描述了本文提出的SSALS-SVM的構建過程,其大致分為以下幾個步驟

    ①隨機生成n個(f,σ2,γ)組合。

    ②將生成的n個(f,σ2,γ)組合輸入目標函數中進行評估,以確定F。該過程需要將控制解釋變量組合的參數f轉化為f個解釋變量構成的特征組合。

    ③更新領導者和跟隨者的位置,并重新對更新后的(f,σ2,γ)組合進行評估。

    ④重復迭代,直到達到最佳適應度或最大迭代次數,輸出最佳超參數組合和最優(yōu)特征組合,進而得到最優(yōu)LS-SVM模型。

    2RC柱抗側移承載力預測

    2.1RC柱抗側移承載力試驗數據集建立

    為了評估提出的SSALS-SVM方法對RC柱抗側移承載力預測的準確性以及泛化能力,本文搜集了248個彎曲破壞柱(延性柱)(194個)、剪切(18個)和彎曲-剪切(36個)破壞柱(非延性柱)的試驗數據集,具體信息見表1。表1包含數據集中具體的26個設計變量(輸入變量),以及各設計變量的平均值與標準差。

    2.2RC柱設計變量的最優(yōu)特征組合

    本節(jié)的主要目的是剔除無關變量,識別出既能影響延性柱抗彎承載力又能影響非延性柱抗剪承載力的設計變量,進而提高預測模型的泛化性能。為此,本文設定SSALS-SVM模型中f的取值范圍為[3,10],這是因為過多或者過少的特征容易導致模型出現過擬合或欠擬合。將σ2和γ的取值范圍設置為[10-6,10]和[10-6,1000 00],[KG)]這是因為σ2的取值過大時,會導致高斯核函數性能下降,而γ的取值越大,預測值與真實值的誤差越小。此外,將d設置為3,種群規(guī)模設置為30,迭代次數設置為100。隨機從26個設計變量中提取f個設計變量(f=3,4,…,10)組成數據集,對每一個f值,共執(zhí)行10次設計變量的隨機組合。為了確立設計變量的數量對RC柱抗側移承載力影響的不確定性,將248個試驗數據按7∶3劃分為訓練集和測試集,為了提高模型的預測精度,在訓練模型的過程中采用十折交叉驗證,即將訓練集十等分,每次選取九份作為開發(fā)子集,剩余一份作為驗證子集,并取十次訓練過程的平均MSE作為此次優(yōu)化的適應度。當達到最大迭代次數或最佳適應度時,輸出最優(yōu)特征組合及最優(yōu)σ2和γ?;谳敵龅淖顑?yōu)參數組合,重新建立SSALS-SVM模型,利用建立的SSALS-SVM模型預測測試集中柱的抗側移承載力,并計算R2、RMSE和MAE,其計算公式如下:

    R2=1-[SX(]∑npi=1(yi-Ay^G1i)2[]∑npi=1(yi-Ay-G1i)2[SX)][JY](11)

    RMSE=[KF(][SX(]∑[DD(]n[]i=1[DD)](yi-Ay^G1i)2[]n[SX)][KF)][JY](12)MAE=[SX(]1[]n[SX)]∑[DD(]n[]i=1[DD)][JB<2|]Ay^G1i-yi[JB>2|][JY](13)

    圖3為不同解釋變量數目f下模型預測結果的R2、RMSE、MAE箱型圖。該箱型圖是在f值一定時,基于10次不同的設計變量組合得到的。其中,箱型圖的四分位距離越長,說明上四分位數與下四分位數之間的差異越大,即不確定性越大,進而說明該f值下,無法建立最優(yōu)特征組合。

    如圖3所示,當f=3時,箱型圖的四分位距離最長,甚至出現異常值,這說明SSALS-SVM與下四分位數之間的差異越大,即不確定性越大,無法識別出柱設計變量的最優(yōu)特征組合。隨著f的取值變大,四分位距離逐漸縮小,即不確定性降低,模型的整體預測性能也在上升。當f=10時,箱型圖中R2的中位值和平均值達到最大,RMSE和MAE的中位值和平均值達到最小。此時,模型的整體預測性能達到最佳且不確定性最低。

    表2為不同挑選數目對應的最佳特征組合的R2、RMSE和MAE,其中,當挑選的特征數目為3時,模型性能較差,說明該數目下的特征組合不足以反映特征與響應變量之間的非線性關系。隨著挑選特征數目的不斷增加,當數目達到10時,模型性能達到最佳。因此,本文選取f=10時,箱型圖中性能最佳所對應的設計變量組合作為最優(yōu)特征組合。該組合包括軸向荷載、角部縱筋直徑、中部縱筋直徑、平行于加載位置處中部縱筋數量、箍筋肢數、箍筋直徑、箍筋數量、箍筋面積配筋率、跨高比和軸壓比。

    2.3RC柱抗側移承載力預測

    為了充分評估SSALS-SVM模型的泛化性能,基于上一節(jié)篩選的最優(yōu)特征組合,本文使用留一法交叉驗證對數據集中柱的抗側移承載力進行預測。留一法交叉驗證是將248組柱試驗數據劃分248份,每次取1份作為測試集,剩下247份作為訓練集,共執(zhí)行248次計算。因此,通過使用留一法交叉驗證,數據集中248個柱試件均作為測試集用于評估本文提出SSALS-SVM模型的泛化性能。使用留一法交叉驗證預測的結果可根據柱破壞類型(即彎曲破壞、剪切破壞和彎曲-剪切破壞)進行劃分,用于后續(xù)與傳統(tǒng)基于力學的計算公式和數據驅動模型得到的結果進行對比,具體過程見下節(jié)。

    3模型泛化性研究

    3.1SSALS-SVM模型與現有方法的對比分析

    基于上一節(jié)留一法預測的結果,本節(jié)主要用其與傳統(tǒng)基于力學的計算公式和數據驅動模型得到的預測結果進行對比分析。對應計算模型的具體數學公式見表3,其中“AMR”指不考慮位移延性時的數據驅動模型(Kakavand et al,2021),“AMRK”指考慮位移延性時的數據驅動模型(Kakavand et al,2021)。表4為不同預測模型得到的R2、RMSE、MAE值,其中粗體值表示結果的最優(yōu)值,“綜合”指的是將彎曲、剪切和彎曲-剪切破壞柱的預測結果合并(即延性柱的抗彎承載力和非延性柱的抗剪承載力),用于量化模型的泛化性能。從表4和圖4a中可知,ACI與OS模型預測性能較差,而相比于這兩種方法,數據驅動模型預測性能有明顯的提高,但SSALS-SVM模型的預測性能更加精確。從表4和圖4中可知,相比于數據驅動模型,SSALS-SVM模型的R2最高提升了近28%,RMSE最高降低了近69%,MAE最高降低了近78%。對于非延性柱,SM模型預測性能優(yōu)于PR模型,而數據驅動模型預測性能明顯高于這兩個模型,能較好地區(qū)分RC柱在地震荷載作用下的破壞模式(Luo,Paal,2018)。但是相比于數據驅動模型,SSALS-SVM模型在預測抗剪承載力時依舊表現出高預測精度,R2最高提升了近15%,RMSE最高降低了近53%,MAE最高降低了35%。

    3.2SSALS-SVM模型的泛化性能分析

    圖5為軸壓比、箍筋面積配筋率和跨高比在不同取值時對模型泛化性能的影響。圖中,預測值與試驗值比值越接近于1,則表明模型預測越準確。從圖中可以看出,SSALS-SVM模型預測值與試驗值比值絕大部分保持在1.0附近,這表明在本文所呈現的軸壓比、箍筋面積配筋率和跨高比范圍內,SSALS-SVM模型能夠準確地預測延性和非延性柱的抗側移承載力,具有較好的泛化性能。相比于其他模型,SSALS-SVM的泛化能力最高提升了近83%。[FL)]

    Priestley等(1994),簡稱PR[]抗剪承載力[]Vn=k[KF(]fc′[KF)]0.8Ag+[SX(]AstfytD'[]s[SX)]+[SX(]h-c[]2a[SX)]P

    其中,當μ≤2時,k=0.29;當μ≥4時,k=0.1;當2≤μ≤4時,k=0.48-0.095μ。

    Sezen和Moehle(2004),簡稱SM[]抗剪承載力[]Vn=[SX(]0.4k[KF(]fc′[KF)][]a/d[SX)][KF(]1+[SX(]P[]0.5[KF(]fc′[KF)]Ag[SX)][KF)]Ag+k[SX(]Astfytd[]s[SX)]

    其中,當μ≤2時,k=1;當2≤μ≤6時,k=1.15-0.075μ;當μ≥6時,k=0.7。

    AMR(2021)[]抗剪承載力[]Vmax=15000+0.44[KF(]fc′[KF)]bd[SX(]d[]a[SX)]+0.17[SX(]Astfytd[]s[SX)]+0.037Aslfyl+54.9[SX(]h-c[]2a[SX)]P[BHG5]

    AMRK(2021)[]抗剪承載力[]Vmax=16800+0.46k[KF(]fc′[KF)]bd[SX(]d[]a[SX)]+0.31k[SX(]Astfytd[]s[SX)]+0.032kAslfyl+57.2k[SX(]h-c[]2a[SX)]P

    其中,當μ≤2時,k=1;當2≤μ≤6時,k=1.1-0.05μ;當μ≥6時,k=0.8。

    Ozbakkaloglu和Saatcioglu(2004),簡稱OS[]抗彎承載力[](1)fc′<30 MPa時,α1=0.85,β1=0.85

    (2)fc′≥30 MPa時,α1=0.85-0.001(fc′-30)≥0.72,β1=0.85-0.002(fc′-30)≥0.67[BH]

    ACI318-22(2014),簡稱ACI[]抗彎承載力[](1)fc′<30 MPa時,α1=0.85,β1=0.85

    (2)fc′≥30 MPa時,α1=0.85,β1=0.85-0.008(fc′-30)≥0.65[BG)F]

    注:Ag為混凝土柱截面面積;Ast箍筋總截面面積;D′為兩側箍筋外邊緣之間的距離;c為中性軸深度;a為混凝土柱的剪跨長度;Asl為縱向鋼筋總截面面積。

    4結論

    本文提出了樽海鞘算法優(yōu)化支持向量機(SSALS-SVM)方法,用于剔除冗余和不相關特征。選取最具代表性、各特征之間相關性弱的特征子集,形成最優(yōu)特征子集,并解決了傳統(tǒng)的基于力學的計算公式和數據驅動模型在預測RC柱抗側移承載力時泛化性能不足的問題?;?48個RC柱試件抗側移承載力試驗數據,分別與傳統(tǒng)的基于力學的計算公式和數據驅動模型進行對比分析。主要結論如下:

    (1)相比于傳統(tǒng)機器學習方法,本文提出的SSALS-SVM方法能夠自動識別并剔除數據集中包含的冗余和不相關特征,篩選出既能影響延性柱抗彎承載力又能影響非延性柱抗剪承載力的設計變量,降低了模型的復雜度。

    (2)基于挑選的最優(yōu)特征組合,與傳統(tǒng)的基于力學的計算公式和數據驅動模型對比,本文提出的SSALS-SVM模型能夠準確預測RC柱抗側移承載力。相比于其他模型,其泛化能力最高提升了83%。

    (3)在本文所呈現的軸壓比、箍筋面積配筋率和跨高比范圍內,SSALS-SVM模型具有較好的泛化性能。

    參考文獻:

    鄧明科,張陽璽,胡紅波.2018.高延性混凝土加固鋼筋混凝土柱抗剪承載力計算[J].工程力學,35(3):159-166.Deng M K,Zhang Y X,Hu H B.2018.Experimental study and calculation of the shear capacity of RC columns strengthened with high ductile concrete[J].Engineering Mechanics,35(3):159-166.(in Chinese)

    甘丹,周緒紅,劉界鵬,等.2018.鋼管約束鋼筋混凝土柱受剪承載力計算[J].建筑結構學報,(9):96-103,111.Gan D,Zhou X H,Liu J P,et al.2018.Calculation for shear strength of reinforced-concrete columns constrained by steel tubes[J].Journal of Building Structures,(9):96-103,111.(in Chinese)

    薛亦聰,楊勇,于云龍.2020.型鋼混凝土短柱受剪承載力計算模型研究[J].建筑結構學報,41(S1):162-170.Xue Y C,Yang Y,Yu Y L.2020.Theoretical model on shear strength of steel reinforced concrete short columns[J].Journal of Building Structures,41(S1):162-170.(in Chinese)

    于曉輝,王猛,寧超列.2022.基于機器學習的鋼筋混凝土柱失效模式兩階段判別方法[J].建筑結構學報,43(8):220.Yu X H,Wang M,Ning C L.2022.A machine-learning-based two-step method for failure mode classification of reinforced concrete columns[J].Journal of Building Structures,43(8):220.(in Chinese)

    張勤,貢金鑫,馬穎.2014.單調和反復荷載作用下彎剪破壞鋼筋混凝土柱荷載-變形關系試驗研究及簡化模型[J].建筑結構學報,35(3):138.Zhang Q,Gong J X,Ma Y.2014.Study on lateral load-deformation relations of flexural-shear failure columns under monotonic and cyclic loading[J].Journal of Building Structures,35(3):138.(in Chinese)

    ACI 318-22,Building code requirements for structural concrete[S].

    Aval S B B,Ketabdari H,Gharebaghi S A.2017.Estimating shear strength of short rectangular reinforced concrete columns using nonlinear regression and gene expression programming[J].Structures,12:13-23.

    Bae S,Bayrak O.2003.Early cover spalling in high-strength concrete columns[J].Journal of Structural Engineering,129(3):314-323.

    Cortes C,Vapnik V.1995.Support vector machine[J].Machine Learning,20:273-297.

    Guyon I,Elisseeff A.2003.An introduction to variable and feature selection[J].Journal of Machine Learning Research,3(Mar):1157-1182.

    Huang C L,Wang C J.2006.A GA-based feature selection and parameters optimizationfor support vector machines[J].Expert Systems with Applications,31(2):231-240.

    Kakavand M R A,Sezen H,Taciroglu E.2021.Data-driven models for predicting the shear strength of rectangular and circular reinforced concrete columns[J].Journal of Structural Engineering,147(1):04020301.

    Ketabdari H,Karimi F,Rasouli M.2020.Shear strength prediction of short circular reinforced-concrete columns using soft computing methods[J].Advances in Structural Engineering,23(14):3048-3061.

    Lee J H,Son H S.2000.Failure and strength of high-strength concrete columns subjected to eccentric loads[J].ACI Structural Journal,97(1):75-85.

    Lin S W,Ying K C,Chen S C,et al.2008.Particle swarm optimization for parameter determination and feature selection of support vector machines[J].Expert Systems with Applications,35(4):1817-1824.

    Luo H,Paal S G.2018.Machine learning–based backbone curve model of reinforced concrete columns subjected to cyclic loading reversals[J].Journal of Computing in Civil Engineering,32(5):04018042.

    Luo H,Paal S G.2021.Metaheuristic least squares support vector machine-based lateral strength modelling of reinforced concrete columns subjected to earthquake loads[J].Structures,33:748-758.

    Mirjalili S,Gandomi A H,Mirjalili S Z,et al.2017.Salp swarm algorithm:A bio-inspired optimizer for engineering design problems[J].Advances in Engineering Software,114:163-191.

    Noroozieh E,Mansouri A.2019.Lateral strength and ductility of reinforced concrete columns strengthened with NSM FRP rebars and FRP jacket[J].International Journal of Advanced Structural Engineering,11(2):195-209.

    Ozbakkaloglu T,Saatcioglu M.2004.Rectangular stress block for high-strength concrete[J].ACI Structural Journal,101(4):475-483.

    Pan Z,Li B.2012.Truss-arch model for shear strength of shear-critical reinforced concrete columns[J].Journal of Structural Engineering,139(4):548-560.

    Priestley M J N,Verma R,Xiao Y.1994.Seismic shear strength of reinforced concrete columns[J].Journal of Structural Engineering,120(8):2310-2329.

    Sezen H,Moehle J P.2004.Shear strength model for lightly reinforced concrete columns[J].Journal of Structural Engineering,130(11):1692-1703.

    Suykens JVA,Van Gestel T,De Brabanter J.2002.Least squares support vector machines[M].Beijing:World Scientific Publishing House.

    Salp Swarm Algorithm-optimized Support Vector Machines For?Lateral Strength Prediction of RC Columns

    OUYANG Qian1,2,LUO Huan1,2

    (1.Hubei Geological Disaster Prevention and Control Engineering Technology Research Center,Yichang 443002,Hubei,China)

    (2.College of Civil Engineering &Architecture,China Three Gorges University,Yichang 443002,Hubei,China)

    Abstract?The existing methods of predicting the lateral strength of RC columns lack generalization performance,where methods of predicting the flexural strength of ductile columns cannot be used to predict the shear strength of non-ductile columns,and vice versa.While current machine learning methods can solve this problem,they cannot automatically remove a large number of redundant and irrelevant features from the dataset,which in turn increases the complexity of the ML model and leads to overfitting.To this end,this paper proposes a new method called salp swarm algorithm-optimized least squares support vector machines(SSALS-SVM),which can remedy the aforementioned problems.Based on a given data set,SSALS-SVM can adopt the salp swarm algorithm(SSA)to automatically eliminate redundant and irrelevant features and select the most representative feature subset with weak correlation among features to form an optimal feature combination,while the hyperparameters governing the nonlinear fitting ability of LS-SVM are also optimized.In this way,the optimized prediction model can not only identify the design variables that influence the lateral strength of ductile and non-ductile columns,but also reflect the nonlinear mapping relationship between the optimal feature combination and lateral strength of RC columns.The generalization performance of proposed SSALS-SVM for predicting the lateral strength of RC columns is verified by comparing with existing prediction models based on 248 experimental data of RC columns.Numerical results show that the generalization performance of proposed SSALS-SVM can be enhanced up to 83% higher than that of existing prediction models.

    Keywords:?reinforced concrete columns;lateral strength;support vector machines;salp swarm algorithm;feature selection

    *收稿日期:2023-09-20.

    基金項目:湖北省自然科學基金面上項目(2022CFB294);國家自然科學基金青年科學項目(52208485).

    第一作者簡介:歐陽謙(1999-),碩士研究生在讀,主要從事結構抗震與機器學習方面的交叉研究.E-mail:2411738467@qq.com.

    通信作者簡介:駱歡(1988-),博士,副教授,主要從事結構抗震與機器學習方面的交叉研究.E-mail:hluo@ctgu.edu.cn.[HT5"SS]歐陽謙,駱歡.2024.樽海鞘算法優(yōu)化支持向量機的RC柱抗側移承載力預測[J].地震研究,47(3):350-358,doi:10.20015/j.cnki.ISSN1000-0666.2024.0051.

    Ouyang Q,Luo H.2024.Salp swarm algorithm-optimized support vector machines for lateral strength prediction of RC columns[J].Journal of Seismological Research,47(3):350-358,doi:10.20015/j.cnki.ISSN1000-0666.2024.0051.

    猜你喜歡
    支持向量機特征選擇
    Kmeans 應用與特征選擇
    電子制作(2017年23期)2017-02-02 07:17:06
    基于改進支持向量機的船舶縱搖預報模型
    中國水運(2016年11期)2017-01-04 12:26:47
    基于SVM的煙草銷售量預測
    軟件導刊(2016年11期)2016-12-22 21:52:38
    動態(tài)場景中的視覺目標識別方法分析
    論提高裝備故障預測準確度的方法途徑
    價值工程(2016年32期)2016-12-20 20:36:43
    基于熵技術的公共事業(yè)費最優(yōu)組合預測
    價值工程(2016年29期)2016-11-14 00:13:35
    基于支持向量機的金融數據分析研究
    基于GA和ELM的電能質量擾動識別特征選擇方法
    聯合互信息水下目標特征選擇算法
    基于特征選擇聚類方法的稀疏TSK模糊系統(tǒng)
    亚洲精品av麻豆狂野| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区不卡视频| 午夜福利在线在线| 久久久久久大精品| 国产日本99.免费观看| 久久午夜综合久久蜜桃| 日本熟妇午夜| 亚洲第一av免费看| 两个人视频免费观看高清| 亚洲欧美激情综合另类| 黑人操中国人逼视频| 身体一侧抽搐| 久久久水蜜桃国产精品网| 18禁裸乳无遮挡免费网站照片 | 欧美 亚洲 国产 日韩一| 一进一出抽搐gif免费好疼| 国产精品,欧美在线| av超薄肉色丝袜交足视频| 成年免费大片在线观看| 欧美 亚洲 国产 日韩一| 免费看a级黄色片| 2021天堂中文幕一二区在线观 | 日日摸夜夜添夜夜添小说| 精品国产亚洲在线| 男女床上黄色一级片免费看| 欧美人与性动交α欧美精品济南到| 国产又黄又爽又无遮挡在线| 国产午夜福利久久久久久| 国产又色又爽无遮挡免费看| 亚洲国产欧美一区二区综合| 成熟少妇高潮喷水视频| 国产乱人伦免费视频| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 久久草成人影院| 香蕉久久夜色| 国产三级在线视频| 欧美日韩乱码在线| 亚洲av电影不卡..在线观看| 99国产精品99久久久久| 男人舔奶头视频| 哪里可以看免费的av片| 人人妻人人澡人人看| 丁香六月欧美| 97碰自拍视频| 少妇裸体淫交视频免费看高清 | 国产精品 欧美亚洲| 色婷婷久久久亚洲欧美| av欧美777| 最近最新中文字幕大全免费视频| √禁漫天堂资源中文www| 美女国产高潮福利片在线看| 精品国产美女av久久久久小说| 国产精品久久久久久人妻精品电影| 国产v大片淫在线免费观看| 国产成人精品久久二区二区免费| 99国产极品粉嫩在线观看| 国产精品日韩av在线免费观看| av福利片在线| 91九色精品人成在线观看| 日韩免费av在线播放| 成人免费观看视频高清| 久99久视频精品免费| 成人一区二区视频在线观看| 91字幕亚洲| 女生性感内裤真人,穿戴方法视频| 国产欧美日韩一区二区精品| 在线国产一区二区在线| 午夜福利18| 日韩欧美免费精品| 欧美中文综合在线视频| 日本黄色视频三级网站网址| 亚洲av日韩精品久久久久久密| 欧美日韩福利视频一区二区| 女人爽到高潮嗷嗷叫在线视频| 午夜免费成人在线视频| 亚洲七黄色美女视频| 精品欧美一区二区三区在线| 熟女电影av网| 99精品欧美一区二区三区四区| 国产色视频综合| 一级毛片女人18水好多| 久久欧美精品欧美久久欧美| 欧美成人性av电影在线观看| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 欧美日韩乱码在线| 亚洲第一电影网av| 国产高清视频在线播放一区| 听说在线观看完整版免费高清| 国产日本99.免费观看| 国产精品亚洲美女久久久| 亚洲成人精品中文字幕电影| 最近最新中文字幕大全免费视频| 正在播放国产对白刺激| 悠悠久久av| 亚洲av日韩精品久久久久久密| 在线看三级毛片| 免费在线观看视频国产中文字幕亚洲| 在线永久观看黄色视频| 中文字幕人妻熟女乱码| 亚洲国产日韩欧美精品在线观看 | 最新美女视频免费是黄的| 不卡一级毛片| 啦啦啦 在线观看视频| 日韩有码中文字幕| 无限看片的www在线观看| 成人三级黄色视频| 男女视频在线观看网站免费 | 黑人巨大精品欧美一区二区mp4| 久久天堂一区二区三区四区| 麻豆国产av国片精品| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 又大又爽又粗| 看黄色毛片网站| 午夜影院日韩av| 大型av网站在线播放| svipshipincom国产片| 国产精品亚洲一级av第二区| 女性生殖器流出的白浆| 特大巨黑吊av在线直播 | 在线观看日韩欧美| 欧美国产精品va在线观看不卡| 日韩精品免费视频一区二区三区| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| 女人高潮潮喷娇喘18禁视频| 高清毛片免费观看视频网站| 色综合婷婷激情| 搡老岳熟女国产| 久久精品成人免费网站| 国产熟女xx| 国产一区在线观看成人免费| 欧美成狂野欧美在线观看| 日日摸夜夜添夜夜添小说| 久久人人精品亚洲av| 夜夜爽天天搞| 国产在线精品亚洲第一网站| 国产三级在线视频| 视频在线观看一区二区三区| 精品高清国产在线一区| 一区福利在线观看| 亚洲人成电影免费在线| 亚洲国产精品合色在线| 男女床上黄色一级片免费看| 亚洲第一欧美日韩一区二区三区| 国产一级毛片七仙女欲春2 | 一区二区三区激情视频| 国产精品自产拍在线观看55亚洲| 国产精品永久免费网站| 午夜福利在线在线| 一本久久中文字幕| 热re99久久国产66热| 亚洲国产精品sss在线观看| 亚洲黑人精品在线| 日日爽夜夜爽网站| 久热这里只有精品99| 午夜久久久在线观看| 日本黄色视频三级网站网址| 亚洲av电影在线进入| 91麻豆av在线| 老司机靠b影院| 国产伦人伦偷精品视频| 在线播放国产精品三级| 悠悠久久av| 欧美日韩乱码在线| 免费观看精品视频网站| av欧美777| 麻豆成人午夜福利视频| 成人国产综合亚洲| 90打野战视频偷拍视频| 国产人伦9x9x在线观看| 午夜福利18| 久久久久国内视频| 可以在线观看毛片的网站| 国产三级在线视频| 久久久久九九精品影院| 亚洲一区二区三区不卡视频| 波多野结衣高清作品| 日韩一卡2卡3卡4卡2021年| 久久久久国产精品人妻aⅴ院| 国产一卡二卡三卡精品| 久久久国产成人精品二区| 黄网站色视频无遮挡免费观看| 亚洲精品国产精品久久久不卡| 亚洲成国产人片在线观看| 黄网站色视频无遮挡免费观看| 嫩草影院精品99| 午夜久久久在线观看| 老熟妇仑乱视频hdxx| 亚洲av熟女| 脱女人内裤的视频| 国产欧美日韩一区二区三| 国产精品美女特级片免费视频播放器 | 久久久国产成人免费| 国产精品亚洲美女久久久| 亚洲成人国产一区在线观看| www国产在线视频色| 久久香蕉精品热| 日韩视频一区二区在线观看| 精品久久久久久久久久免费视频| 国产三级在线视频| 给我免费播放毛片高清在线观看| 我的亚洲天堂| 99国产极品粉嫩在线观看| 母亲3免费完整高清在线观看| 成人国语在线视频| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| 久久国产精品人妻蜜桃| 久久人人精品亚洲av| 麻豆av在线久日| 黑人操中国人逼视频| 国产私拍福利视频在线观看| 免费看日本二区| 黄网站色视频无遮挡免费观看| 俄罗斯特黄特色一大片| 日日爽夜夜爽网站| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 熟妇人妻久久中文字幕3abv| 国产精品九九99| 亚洲avbb在线观看| 国产男靠女视频免费网站| 黄片大片在线免费观看| 久久久久国产精品人妻aⅴ院| 亚洲熟妇熟女久久| 在线播放国产精品三级| 中文字幕人成人乱码亚洲影| 无人区码免费观看不卡| 黄色成人免费大全| 国产黄色小视频在线观看| 人人澡人人妻人| 丁香欧美五月| 中文字幕精品亚洲无线码一区 | 久久久精品欧美日韩精品| 国产精品永久免费网站| 欧美在线黄色| 99久久精品国产亚洲精品| 亚洲人成网站在线播放欧美日韩| 亚洲黑人精品在线| 欧美性长视频在线观看| 一个人免费在线观看的高清视频| 欧美一级a爱片免费观看看 | 亚洲色图 男人天堂 中文字幕| 欧美性长视频在线观看| 国产精品精品国产色婷婷| 国产又黄又爽又无遮挡在线| 午夜福利一区二区在线看| 精品一区二区三区四区五区乱码| 久久婷婷成人综合色麻豆| 大型av网站在线播放| 亚洲国产欧美网| 免费电影在线观看免费观看| 美女扒开内裤让男人捅视频| 亚洲自拍偷在线| 久久性视频一级片| 亚洲人成网站在线播放欧美日韩| a级毛片在线看网站| 9191精品国产免费久久| 色综合站精品国产| 91大片在线观看| 美女午夜性视频免费| 亚洲成人久久爱视频| 国产亚洲精品久久久久久毛片| 久久人妻福利社区极品人妻图片| 国产激情偷乱视频一区二区| 变态另类丝袜制服| 精品国产乱子伦一区二区三区| 久久狼人影院| 特大巨黑吊av在线直播 | 色尼玛亚洲综合影院| 日本五十路高清| 青草久久国产| 国产精品99久久99久久久不卡| 中文字幕另类日韩欧美亚洲嫩草| 看黄色毛片网站| 国产精品一区二区三区四区久久 | 国产亚洲精品第一综合不卡| 一进一出抽搐动态| 999精品在线视频| 1024视频免费在线观看| 亚洲午夜精品一区,二区,三区| 91成年电影在线观看| 九色国产91popny在线| www日本黄色视频网| 可以在线观看的亚洲视频| www日本在线高清视频| 亚洲久久久国产精品| 少妇的丰满在线观看| 亚洲av成人不卡在线观看播放网| 一进一出抽搐动态| 国产精品综合久久久久久久免费| 久久九九热精品免费| 手机成人av网站| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| 亚洲久久久国产精品| 黄色丝袜av网址大全| 亚洲国产看品久久| 欧美一级毛片孕妇| 无遮挡黄片免费观看| 中文字幕精品亚洲无线码一区 | 亚洲av电影在线进入| 99精品在免费线老司机午夜| 国产蜜桃级精品一区二区三区| 欧美中文日本在线观看视频| 亚洲自拍偷在线| 午夜激情福利司机影院| 日本免费一区二区三区高清不卡| 国产熟女午夜一区二区三区| 亚洲一区二区三区色噜噜| 色精品久久人妻99蜜桃| 在线观看66精品国产| 亚洲aⅴ乱码一区二区在线播放 | 日韩 欧美 亚洲 中文字幕| 久久久久久亚洲精品国产蜜桃av| 最好的美女福利视频网| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 亚洲一卡2卡3卡4卡5卡精品中文| 每晚都被弄得嗷嗷叫到高潮| 男人操女人黄网站| 不卡一级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久人人人人人| 日本一区二区免费在线视频| 成人国产一区最新在线观看| 女警被强在线播放| 老汉色∧v一级毛片| 国产精华一区二区三区| 黄网站色视频无遮挡免费观看| 久久精品91无色码中文字幕| 国产精品免费视频内射| av超薄肉色丝袜交足视频| 日韩精品中文字幕看吧| 黑丝袜美女国产一区| 成人欧美大片| 国产精品一区二区三区四区久久 | av欧美777| 精品日产1卡2卡| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 日韩精品青青久久久久久| 久久久国产成人免费| 国产成人欧美| 日本一区二区免费在线视频| 丝袜在线中文字幕| 亚洲人成电影免费在线| 日韩精品青青久久久久久| 日本一区二区免费在线视频| 波多野结衣高清无吗| 亚洲人成电影免费在线| 波多野结衣高清无吗| 三级毛片av免费| 久久久久久九九精品二区国产 | 国产1区2区3区精品| 长腿黑丝高跟| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| 国产亚洲精品一区二区www| 淫秽高清视频在线观看| 91在线观看av| 18禁美女被吸乳视频| 中文资源天堂在线| 中文字幕另类日韩欧美亚洲嫩草| 欧美激情久久久久久爽电影| 中亚洲国语对白在线视频| 激情在线观看视频在线高清| 视频在线观看一区二区三区| 亚洲精品在线美女| 麻豆国产av国片精品| 亚洲成人久久性| 国产真实乱freesex| 久久久久久九九精品二区国产 | 亚洲国产欧美网| 天天添夜夜摸| 一本大道久久a久久精品| 丝袜美腿诱惑在线| 不卡av一区二区三区| 国产黄色小视频在线观看| 天天一区二区日本电影三级| 69av精品久久久久久| 色老头精品视频在线观看| 免费人成视频x8x8入口观看| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看 | 手机成人av网站| 国产精品 国内视频| 精品不卡国产一区二区三区| 精品欧美国产一区二区三| 国产高清视频在线播放一区| aaaaa片日本免费| 亚洲自拍偷在线| 国产v大片淫在线免费观看| 亚洲国产看品久久| 欧美大码av| 亚洲成av片中文字幕在线观看| 欧美最黄视频在线播放免费| 久久久国产精品麻豆| 中文字幕高清在线视频| 亚洲九九香蕉| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 欧美日韩精品网址| 成人手机av| 麻豆久久精品国产亚洲av| 最近在线观看免费完整版| 悠悠久久av| 久久久精品国产亚洲av高清涩受| 美女大奶头视频| 91麻豆av在线| 国产一卡二卡三卡精品| 成人欧美大片| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩精品网址| 桃红色精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 两人在一起打扑克的视频| 波多野结衣av一区二区av| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩无卡精品| 亚洲男人的天堂狠狠| 久久99热这里只有精品18| 日韩欧美国产在线观看| 国产一区二区在线av高清观看| 男男h啪啪无遮挡| 久久人妻福利社区极品人妻图片| 精品福利观看| 后天国语完整版免费观看| 亚洲 国产 在线| 在线看三级毛片| 国产真人三级小视频在线观看| 国产精品永久免费网站| 久久精品影院6| 国产亚洲精品综合一区在线观看 | 中出人妻视频一区二区| 香蕉av资源在线| 黄色a级毛片大全视频| videosex国产| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 国产精品久久久久久亚洲av鲁大| 成人三级做爰电影| 88av欧美| 韩国av一区二区三区四区| 国产成人精品无人区| 国产视频一区二区在线看| 欧美人与性动交α欧美精品济南到| 久久久久亚洲av毛片大全| 首页视频小说图片口味搜索| 国产精品香港三级国产av潘金莲| 精品一区二区三区av网在线观看| 天天添夜夜摸| 非洲黑人性xxxx精品又粗又长| 国产三级在线视频| 欧美国产日韩亚洲一区| 亚洲精品粉嫩美女一区| 99国产精品一区二区三区| 99久久无色码亚洲精品果冻| 在线观看www视频免费| 无人区码免费观看不卡| 男人舔女人下体高潮全视频| 黄色a级毛片大全视频| 一个人免费在线观看的高清视频| 国产精品国产高清国产av| 在线观看免费日韩欧美大片| 午夜福利在线在线| 国产精品 欧美亚洲| 天堂√8在线中文| 一本综合久久免费| 丝袜人妻中文字幕| 此物有八面人人有两片| 老司机午夜福利在线观看视频| 香蕉av资源在线| 国产1区2区3区精品| 久久人妻福利社区极品人妻图片| 在线十欧美十亚洲十日本专区| 久久天躁狠狠躁夜夜2o2o| 久久久久久免费高清国产稀缺| 精品无人区乱码1区二区| 在线观看一区二区三区| 国产在线精品亚洲第一网站| 久久中文字幕一级| 亚洲熟女毛片儿| 99热6这里只有精品| www日本黄色视频网| 男女那种视频在线观看| 丝袜美腿诱惑在线| 熟女电影av网| av片东京热男人的天堂| 日韩欧美一区视频在线观看| 中文字幕久久专区| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸| x7x7x7水蜜桃| 免费av毛片视频| 一级毛片高清免费大全| 最近最新中文字幕大全免费视频| 欧美亚洲日本最大视频资源| 成年女人毛片免费观看观看9| 色综合亚洲欧美另类图片| 久久国产亚洲av麻豆专区| 精品免费久久久久久久清纯| 国产99白浆流出| 在线av久久热| 1024香蕉在线观看| 午夜日韩欧美国产| 中国美女看黄片| 97碰自拍视频| 成人18禁在线播放| 成人免费观看视频高清| 777久久人妻少妇嫩草av网站| 1024手机看黄色片| 国产一区在线观看成人免费| 脱女人内裤的视频| cao死你这个sao货| 悠悠久久av| 精品日产1卡2卡| 老汉色av国产亚洲站长工具| 岛国在线观看网站| 制服诱惑二区| 美女午夜性视频免费| 啦啦啦观看免费观看视频高清| 韩国精品一区二区三区| 在线播放国产精品三级| 国产一区二区三区在线臀色熟女| 欧美精品亚洲一区二区| 精品无人区乱码1区二区| 国产精品久久久人人做人人爽| 成人精品一区二区免费| 中国美女看黄片| 99国产极品粉嫩在线观看| 伦理电影免费视频| 欧美性猛交╳xxx乱大交人| 午夜两性在线视频| 欧美中文日本在线观看视频| 一二三四社区在线视频社区8| 天天一区二区日本电影三级| 一级作爱视频免费观看| 久久精品夜夜夜夜夜久久蜜豆 | 婷婷亚洲欧美| 1024手机看黄色片| 国产精品香港三级国产av潘金莲| 久久 成人 亚洲| 欧美成人性av电影在线观看| 老熟妇仑乱视频hdxx| 午夜久久久在线观看| 久久久国产精品麻豆| 日本成人三级电影网站| 久久久精品国产亚洲av高清涩受| 久久天堂一区二区三区四区| 国产亚洲精品一区二区www| 青草久久国产| 国产一区在线观看成人免费| 国产1区2区3区精品| 中国美女看黄片| 国产精品久久久久久亚洲av鲁大| 日韩欧美国产在线观看| 极品教师在线免费播放| 国产99白浆流出| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久 | 叶爱在线成人免费视频播放| 亚洲三区欧美一区| 精品久久蜜臀av无| 国产三级黄色录像| 亚洲五月色婷婷综合| 亚洲一区高清亚洲精品| 长腿黑丝高跟| 手机成人av网站| 免费在线观看成人毛片| 欧美黑人欧美精品刺激| 激情在线观看视频在线高清| 热99re8久久精品国产| 草草在线视频免费看| 中文资源天堂在线| 男人操女人黄网站| 国产主播在线观看一区二区| 99精品在免费线老司机午夜| 国产精品久久视频播放| 亚洲成人精品中文字幕电影| 精品国内亚洲2022精品成人| 人成视频在线观看免费观看| av天堂在线播放| 国产又黄又爽又无遮挡在线| 狠狠狠狠99中文字幕| 欧美+亚洲+日韩+国产| 精品久久久久久久久久久久久 | 国产成人精品久久二区二区免费| 午夜免费激情av| 99re在线观看精品视频| 国产成人精品久久二区二区免费| 大香蕉久久成人网| 美女免费视频网站| 一区二区三区国产精品乱码| 亚洲五月婷婷丁香| 男人舔女人下体高潮全视频| 亚洲午夜理论影院| 男人操女人黄网站| 国产免费av片在线观看野外av| 无限看片的www在线观看| 久99久视频精品免费| 91麻豆精品激情在线观看国产| 无限看片的www在线观看| 操出白浆在线播放| 免费av毛片视频| 亚洲最大成人中文| 久久这里只有精品19| 中文字幕人妻丝袜一区二区| 国产精品98久久久久久宅男小说| 久久久久久亚洲精品国产蜜桃av| 亚洲精品粉嫩美女一区|