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

    基于機(jī)器學(xué)習(xí)的高強(qiáng)鋼焊接等截面箱型柱整體穩(wěn)定性預(yù)測(cè)方法

    2024-01-03 11:49:31張營(yíng)營(yíng)徐浩陳培見(jiàn)馬俊周祎
    土木建筑與環(huán)境工程 2024年1期
    關(guān)鍵詞:預(yù)測(cè)模型機(jī)器學(xué)習(xí)

    張營(yíng)營(yíng) 徐浩 陳培見(jiàn) 馬俊 周祎

    摘要:目前,針對(duì)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性的研究多采用有限元建?;?qū)嶒?yàn)室試驗(yàn)方法,而基于機(jī)器學(xué)習(xí)的預(yù)測(cè)方法能夠顯著提升預(yù)測(cè)的準(zhǔn)確性和便捷性。為了準(zhǔn)確預(yù)測(cè)高強(qiáng)鋼焊接等截面箱型柱的整體穩(wěn)定性,提出使用纖維模型構(gòu)建數(shù)據(jù)庫(kù)并利用機(jī)器學(xué)習(xí)建立預(yù)測(cè)模型的方法。首先確定模型的輸入輸出參數(shù),并通過(guò)纖維模型方法建立數(shù)據(jù)庫(kù);接著,選用常見(jiàn)的3種不同類型的機(jī)器學(xué)習(xí)模型和現(xiàn)有規(guī)范中的經(jīng)驗(yàn)?zāi)P瓦M(jìn)行預(yù)測(cè),并依據(jù)評(píng)價(jià)指標(biāo)進(jìn)行性能對(duì)比;最后,根據(jù)可解釋算法分析機(jī)器學(xué)習(xí)模型的合理性。結(jié)果表明:大部分機(jī)器學(xué)習(xí)模型預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果吻合度略高于現(xiàn)有規(guī)范中的經(jīng)驗(yàn)?zāi)P停渲?,高斯過(guò)程回歸模型對(duì)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性的預(yù)測(cè)表現(xiàn)最優(yōu);機(jī)器學(xué)習(xí)預(yù)測(cè)模型中各類參數(shù)對(duì)構(gòu)件整體穩(wěn)定性的影響趨勢(shì)符合預(yù)期,驗(yàn)證了機(jī)器學(xué)習(xí)模型的合理性和可靠性;構(gòu)件的正則化長(zhǎng)細(xì)比對(duì)預(yù)測(cè)結(jié)果影響最大,而構(gòu)件初始缺陷的影響相對(duì)最小。

    關(guān)鍵詞:機(jī)器學(xué)習(xí);高強(qiáng)鋼;整體穩(wěn)定性;預(yù)測(cè)模型;纖維模型

    中圖分類號(hào):TU391? ? ?文獻(xiàn)標(biāo)志碼:A? ? ?文章編號(hào):2096-6717(2024)01-0182-12

    Machine learning method for overall stability of welded constant section box columns made of high strength steel

    ZHANG Yingying1a,1b, XU Hao1a, CHEN Peijian1a, MA Jun2, ZHOU Yi3

    (1a. School of Mechanics and Civil Engineering; 1b. Jiangsu Key Laboratory of Environmental Disaster and Structural Reliability of Civil Engineering, China University of Mining and Technology, Xuzhou 221116, Jiangsu, P. R. China; 2. South Branch of China Construction Eighth Engineering Bureau Co., Ltd, Shenzhen 518035, Guangdong, P. R. China; 3. College of Civil Engineering, Southwest Jiaotong University, Chendu 610031, P. R. China)

    Abstract: At present, finite element modeling or laboratory testing methods are generally used in the research of the overall stability of high-strength steel members. However, the prediction method based on machine learning (ML) has greatly improved the accuracy and convenience of component performance prediction. To accurately predict the overall stability of welded constant section box columns made of high strength steel,? ML method together with a database based on the fiber model is proposed in this paper. Firstly, the input and output parameters of the model are determined, and the database is provided. Then, three different ML models and empirical models in the existing specifications are selected for prediction, and the performance is compared according to the evaluation index. Finally, the rationality of ML models is analyzed according to interpretable algorithms. The results show that the prediction results of most ML models are in good agreement with the experimental results, which are slightly higher than the empirical models, and the Gaussian process regression model has the best prediction performance for the overall stability of high-strength steel members; the influential trend of various parameters on the overall stability of components meets the expectation, which verifies the rationality and reliability of the ML model; the regularized slenderness ratio has the greatest influence on the prediction results, while the initial defects have the least.

    Keywords: machine learning; high-strength steel; overall stability; prediction model; fiber model

    近年來(lái),隨著建筑產(chǎn)業(yè)的高速發(fā)展,人們開(kāi)始建造更多大跨度、超高層、大重量的建筑物。高強(qiáng)鋼材料具有更高的屈服強(qiáng)度和抗拉強(qiáng)度,能代替普通強(qiáng)度鋼作為大跨高層建筑的建筑材料來(lái)達(dá)到節(jié)約鋼材、減少造價(jià)、提高建筑效率等目的。隨著高強(qiáng)鋼產(chǎn)量需求的增加,高強(qiáng)鋼構(gòu)件的受力性能需要進(jìn)一步研究[1]。高強(qiáng)鋼軸心受壓構(gòu)件受力性能中的穩(wěn)定性一直是高強(qiáng)鋼應(yīng)用和研究的重點(diǎn)。近年來(lái),很多學(xué)者進(jìn)行了高強(qiáng)鋼構(gòu)件的整體穩(wěn)定性研究。薛加燁[2]、邱林波等[3]對(duì)460、550、690 MPa的箱型和H型截面構(gòu)件進(jìn)行軸壓試驗(yàn),并驗(yàn)證了高強(qiáng)鋼構(gòu)件整體穩(wěn)定承載力相對(duì)普通強(qiáng)度鋼的優(yōu)越性;趙金友等[4]對(duì)跨中無(wú)側(cè)向支撐的3個(gè)雙軸對(duì)稱和6個(gè)單軸對(duì)稱焊接工字形截面簡(jiǎn)支梁進(jìn)行了整體彎扭屈曲試驗(yàn)以研究Q460高強(qiáng)鋼焊接工字形截面簡(jiǎn)支梁的整體穩(wěn)定性能并對(duì)《高強(qiáng)鋼結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)》[5]中的簡(jiǎn)支梁整體穩(wěn)定系數(shù)計(jì)算公式進(jìn)行修正。Meng等[6]對(duì)結(jié)構(gòu)鋼圓形空心截面(CHS)梁柱的整體屈曲性能展開(kāi)試驗(yàn)與數(shù)值模擬研究,并對(duì)歐洲規(guī)范中擬議的設(shè)計(jì)方法進(jìn)行了可靠性評(píng)估。學(xué)者們?cè)趯?duì)高強(qiáng)鋼構(gòu)件的整體穩(wěn)定性進(jìn)行研究時(shí),大多采用有限元分析或現(xiàn)場(chǎng)試驗(yàn)的研究方法,無(wú)論從時(shí)間成本還是從人力成本上考慮,這兩種方法都在實(shí)際研究中都有不便之處。因此,有必要探尋一種高效且準(zhǔn)確的高強(qiáng)鋼構(gòu)件整體穩(wěn)定性能預(yù)測(cè)方法。

    作為人工智能領(lǐng)域分支之一的機(jī)器學(xué)習(xí)方法憑借其優(yōu)秀的預(yù)測(cè)精度、便捷的建模流程及可靠的理論支持,在土木領(lǐng)域迅速興起,并在建筑構(gòu)件的性能預(yù)測(cè)問(wèn)題上表現(xiàn)優(yōu)異。Limbachiya等[7]利用人工神經(jīng)網(wǎng)絡(luò)(ANN)預(yù)測(cè)蜂窩梁的腹板屈曲后剪切強(qiáng)度,驗(yàn)證了人工神經(jīng)網(wǎng)絡(luò)模型準(zhǔn)確性的同時(shí)說(shuō)明了機(jī)器學(xué)習(xí)預(yù)測(cè)模型與現(xiàn)有的設(shè)計(jì)規(guī)范模型相比性能更優(yōu);Wakjira等[8]對(duì)比不同的機(jī)器學(xué)習(xí)預(yù)測(cè)模型性能以提出織物增強(qiáng)水泥基加固梁的最佳預(yù)測(cè)模型,并對(duì)現(xiàn)有相關(guān)規(guī)范準(zhǔn)則模型進(jìn)行評(píng)價(jià);Sarothi等[9]首次編制了包含443個(gè)實(shí)驗(yàn)數(shù)據(jù)集組成數(shù)據(jù)庫(kù),結(jié)合機(jī)器學(xué)習(xí)預(yù)測(cè)結(jié)構(gòu)鋼雙剪螺栓連接的承載強(qiáng)度,并對(duì)其特征重要性進(jìn)行分析,驗(yàn)證了鋼材極限屈服強(qiáng)度比和螺栓排數(shù)對(duì)連接強(qiáng)度的影響顯著。機(jī)器學(xué)習(xí)能從海量數(shù)據(jù)中提取其中特征并尋求特征間的規(guī)律來(lái)達(dá)到數(shù)據(jù)預(yù)測(cè)的目的。當(dāng)遇到的任務(wù)非常復(fù)雜或者需要通過(guò)計(jì)算機(jī)自動(dòng)調(diào)整時(shí),就需要借助強(qiáng)大的機(jī)器學(xué)習(xí)技術(shù)來(lái)完成[10]。但用于預(yù)測(cè)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性的機(jī)器學(xué)習(xí)預(yù)測(cè)模型還尚待研究,且針對(duì)帶有初始缺陷的高強(qiáng)鋼構(gòu)件的整體穩(wěn)定性問(wèn)題尚缺乏可靠規(guī)范進(jìn)行指導(dǎo)。

    筆者研究開(kāi)發(fā)高強(qiáng)鋼焊接等截面箱型柱整體穩(wěn)定性的機(jī)器學(xué)習(xí)預(yù)測(cè)模型,對(duì)不同強(qiáng)度的高強(qiáng)鋼焊接等截面箱型柱,使用纖維模型算法創(chuàng)建數(shù)據(jù)庫(kù),參考現(xiàn)有鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范中對(duì)鋼結(jié)構(gòu)整體穩(wěn)定性的描述;將構(gòu)件總寬度、翼緣厚度、計(jì)算長(zhǎng)度、初始缺陷值、屈服強(qiáng)度、正則化長(zhǎng)細(xì)比、寬厚比、彈性模量這8個(gè)參數(shù)作為輸入?yún)?shù);以構(gòu)件極限承載力作為輸出參數(shù),并對(duì)比分析3種優(yōu)化后的常用的機(jī)器學(xué)習(xí)模型的預(yù)測(cè)結(jié)果。通過(guò)評(píng)價(jià)指標(biāo)評(píng)選出最優(yōu)預(yù)測(cè)模型并驗(yàn)證機(jī)器學(xué)習(xí)模型相較于現(xiàn)有規(guī)范中經(jīng)驗(yàn)?zāi)P偷膬?yōu)越性。最后,通過(guò)解釋性算法對(duì)機(jī)器學(xué)習(xí)預(yù)測(cè)模型做出解釋,在驗(yàn)證模型準(zhǔn)確性的基礎(chǔ)上,根據(jù)解釋內(nèi)容結(jié)合規(guī)范及已有研究驗(yàn)證模型的合理性。

    1 訓(xùn)練數(shù)據(jù)庫(kù)的建立

    1.1 纖維模型簡(jiǎn)介

    通常來(lái)說(shuō),對(duì)機(jī)器學(xué)習(xí)預(yù)測(cè)模型進(jìn)行訓(xùn)練時(shí)需要一個(gè)數(shù)據(jù)可靠且數(shù)據(jù)容量足夠大的數(shù)據(jù)庫(kù),以便機(jī)器學(xué)習(xí)模型進(jìn)行訓(xùn)練、驗(yàn)證和測(cè)試。常見(jiàn)的構(gòu)建數(shù)據(jù)庫(kù)的手段包括有限元模擬、現(xiàn)場(chǎng)試驗(yàn)以及搜集已發(fā)表論文中的成果數(shù)據(jù)。有限元的數(shù)值模擬方法雖然得到的數(shù)據(jù)足夠精確,但存在建模困難,耗時(shí)較長(zhǎng)等問(wèn)題?,F(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)獲取不僅需要耗費(fèi)大量人力和時(shí)間進(jìn)行試驗(yàn),且得到的數(shù)據(jù)可能會(huì)因?yàn)樵囼?yàn)環(huán)境等問(wèn)題產(chǎn)生噪聲數(shù)據(jù)。由于數(shù)據(jù)庫(kù)的數(shù)據(jù)需求量較大,搜集大量論文中的已有試驗(yàn)成果數(shù)據(jù)同樣費(fèi)時(shí)費(fèi)力。

    纖維模型的基本思想是將構(gòu)件長(zhǎng)度方向的截面劃分為網(wǎng)格狀的矩形纖維再通過(guò)數(shù)值積分的方法進(jìn)行計(jì)算。該方法需要分別建立截面內(nèi)力、外力和變形之間的關(guān)系。綜合考慮力的平衡與變形的協(xié)調(diào)來(lái)逐步得到荷載-位移曲線及其數(shù)值計(jì)算結(jié)果。該方法存在4個(gè)基本假定:不考慮局部屈曲和橫向扭轉(zhuǎn)屈曲、不考慮剪力影響、平截面假定、構(gòu)件屈曲后變形曲線為正弦曲線的半波[11]。高強(qiáng)鋼焊接箱型柱纖維模型的本構(gòu)模型在相關(guān)文獻(xiàn)中有所體現(xiàn)[12]。由于殘余應(yīng)力對(duì)軸向受壓結(jié)構(gòu)的屈曲行為有較大影響,所以纖維模型計(jì)算時(shí)將考慮班慧勇等[13]提出的焊接箱形截面統(tǒng)一殘余應(yīng)力模型。

    根據(jù)班慧勇等對(duì)于高強(qiáng)鋼焊接截面殘余應(yīng)力統(tǒng)一分布模型的研究可知,高強(qiáng)鋼焊接箱形截面的殘余應(yīng)力分布存在如下特點(diǎn):1)板件中心一定區(qū)域存在均勻分布的殘余壓應(yīng)力,而焊縫附件區(qū)域存在殘余拉應(yīng)力。2)焊縫附近區(qū)域的最大殘余拉應(yīng)力遠(yuǎn)小于鋼材屈服強(qiáng)度,靠近角部的殘余拉應(yīng)力數(shù)值相對(duì)偏小。3)板件中部區(qū)域的殘余壓應(yīng)力隨板件寬厚比增大而明顯減小。根據(jù)班慧勇等[13]的研究結(jié)果,鋼材焊接箱形截面的殘余應(yīng)力分布可以簡(jiǎn)化為如圖1所示的階梯狀分布形式在纖維模型中使用。

    考慮到纖維模型在運(yùn)算速度、建模難度方面的顯著優(yōu)勢(shì),很多學(xué)者都將纖維模型方法應(yīng)用于高強(qiáng)鋼構(gòu)件的研究中[14-15]。此外,纖維模型在分析高強(qiáng)鋼整體穩(wěn)定性能方面的適用性與優(yōu)越性均已經(jīng)得到了驗(yàn)證,有學(xué)者通過(guò)對(duì)比實(shí)驗(yàn)確定了纖維模型對(duì)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性預(yù)測(cè)結(jié)果與實(shí)際試驗(yàn)值的平均誤差僅為9.49%[16]。纖維模型方法不僅能夠快速構(gòu)建容量足夠的數(shù)據(jù)庫(kù),其中的數(shù)據(jù)質(zhì)量也能得到保證。因此,采用纖維模型方法作為數(shù)據(jù)庫(kù)的主要構(gòu)建方法,構(gòu)建了548組樣本數(shù)據(jù)供模型進(jìn)行訓(xùn)練。

    1.2 參數(shù)選取說(shuō)明

    考慮到鋼構(gòu)件的軸心受壓承載能力能直觀反映構(gòu)件的整體穩(wěn)定性,因此,選取構(gòu)件在荷載作用下的極限承載力作為輸出參數(shù)。

    根據(jù)現(xiàn)行《高強(qiáng)鋼結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)》(JGJ/T 483—2020)[5]可知,鋼構(gòu)件的極限承載力與構(gòu)件總寬、總高、翼緣厚度、腹板厚度、計(jì)算長(zhǎng)度、面積、慣性矩、回轉(zhuǎn)半徑、屈服強(qiáng)度、長(zhǎng)細(xì)比、正則化長(zhǎng)細(xì)比、寬厚比、彈性模量、初始缺陷等均有關(guān)聯(lián),因此,將上述參數(shù)作為輸入?yún)?shù)構(gòu)建纖維模型數(shù)據(jù)庫(kù)。在通過(guò)數(shù)據(jù)庫(kù)訓(xùn)練數(shù)據(jù)之前,采用Pearson相關(guān)系數(shù)法對(duì)輸入?yún)?shù)的相關(guān)性進(jìn)行分析,從而篩選掉相關(guān)性過(guò)大的參數(shù)(Pearson系數(shù)大于0.8)。上述涉及參數(shù)的相關(guān)性分析結(jié)果如圖2所示。

    方格顏色越深,數(shù)字越接近1,代表數(shù)據(jù)間的相關(guān)性越大。由于是對(duì)高強(qiáng)鋼焊接等截面箱型柱進(jìn)行研究,其對(duì)稱的特征也極大程度體現(xiàn)在部分幾何參數(shù)和材性參數(shù)中,所以相關(guān)性分析剔除掉了對(duì)稱數(shù)據(jù)中的某一方,例如,構(gòu)件總高和總寬相等,兩者只選構(gòu)件總寬數(shù)據(jù)作為輸入?yún)?shù)即可。依據(jù)相關(guān)性分析結(jié)果,最終選定構(gòu)件總寬、腹板厚度、計(jì)算長(zhǎng)度、初始缺陷值、屈服強(qiáng)度、正則化長(zhǎng)細(xì)比、寬厚比、彈性模量這8個(gè)參數(shù)作為輸入?yún)?shù)。

    此外,在構(gòu)建數(shù)據(jù)庫(kù)時(shí),還要注意保證訓(xùn)練樣本參數(shù)數(shù)據(jù)在一定范圍內(nèi)的離散性,只有訓(xùn)練參數(shù)足夠離散,才能使得機(jī)器學(xué)習(xí)模型的預(yù)測(cè)結(jié)果是考慮到所有的輸入?yún)?shù)得到的。數(shù)據(jù)集中特征參數(shù)的基本信息如表1所示。

    2 基于機(jī)器學(xué)習(xí)和規(guī)范準(zhǔn)則的模型

    2.1 機(jī)器學(xué)習(xí)模型

    機(jī)器學(xué)習(xí)(ML)模型是從一些不基于物理定律的數(shù)學(xué)方程發(fā)展而來(lái)的。機(jī)器學(xué)習(xí)模型訓(xùn)練、驗(yàn)證所依賴的數(shù)據(jù)庫(kù)均為上文提到的纖維模型數(shù)據(jù)庫(kù)。借助MATLAB軟件,使用隨機(jī)森林回歸(RF)、高斯過(guò)程回歸(GPR)和人工神經(jīng)網(wǎng)絡(luò)(ANN)這3種ML模型進(jìn)行訓(xùn)練,其中包括兩個(gè)單一模型和一個(gè)集成模型。在之前的研究中,大部分研究者只考慮了一種機(jī)器學(xué)習(xí)算法對(duì)于構(gòu)件性能預(yù)測(cè)的優(yōu)越性,主要是人工神經(jīng)網(wǎng)絡(luò)[16-18]。然而,所有的ML模型都有其優(yōu)點(diǎn)和缺點(diǎn),面對(duì)不同輸入?yún)?shù)及不同數(shù)據(jù)特征的構(gòu)件性能預(yù)測(cè)問(wèn)題,最適用的ML模型并不唯一。因此,有必要對(duì)眾多的ML模型進(jìn)行擇優(yōu)選取。選擇最合適的ML模型的方法就是基于評(píng)價(jià)指標(biāo)的試錯(cuò)過(guò)程。選取的3種機(jī)器學(xué)習(xí)方法在常規(guī)機(jī)器學(xué)習(xí)算法中屬較新研究成果,且在其他類型構(gòu)件性能的預(yù)測(cè)中也有良好表現(xiàn)。其他諸如SVM、線性回歸等簡(jiǎn)單機(jī)器學(xué)習(xí)方法,作者也進(jìn)行過(guò)預(yù)訓(xùn)練,但效果遠(yuǎn)不如所述3種方法,因此,僅對(duì)上述3種機(jī)器學(xué)習(xí)方法進(jìn)行分析研究。

    模型優(yōu)化方面,采用k-折疊交叉驗(yàn)證方法以避免訓(xùn)練過(guò)程中出現(xiàn)過(guò)度擬合的情況。k-折疊交叉驗(yàn)證就是將數(shù)據(jù)集隨機(jī)平均分為k份。輪流取其中1份作為測(cè)試集,其余數(shù)據(jù)作為訓(xùn)練集進(jìn)行試驗(yàn)。k值的選擇通常是5或10[19-20],當(dāng)數(shù)據(jù)庫(kù)內(nèi)樣本足夠多時(shí),k取5即可。超參數(shù)調(diào)整是模型優(yōu)化的一項(xiàng)重要工作,其目的是全面搜尋使得模型預(yù)測(cè)誤差最低的超參數(shù)或超參數(shù)組合。在一些優(yōu)秀的超參數(shù)算法提出之前,超參數(shù)只能通過(guò)人工調(diào)整、不斷試錯(cuò)的方法進(jìn)行調(diào)優(yōu),但此種方法不僅效率低下,得到的超參數(shù)組合也不一定能使得模型性能得到顯著提升。目前常見(jiàn)的實(shí)現(xiàn)自動(dòng)化超參數(shù)調(diào)優(yōu)的算法有網(wǎng)格搜索、隨機(jī)搜索和貝葉斯優(yōu)化。相較其他兩種算法,貝葉斯優(yōu)化算法的迭代次數(shù)較少,算法效率也相對(duì)較高,因此,使用貝葉斯優(yōu)化方法進(jìn)行超參數(shù)調(diào)優(yōu)。貝葉斯超參數(shù)調(diào)優(yōu)算法流程如圖3所示。圖4展示了模型的均方誤差隨著貝葉斯優(yōu)化迭代次數(shù)的增加逐步降低,同時(shí)模型訓(xùn)練性能逐步提高的趨勢(shì)。在第10次迭代之后,模型性能趨于穩(wěn)定。

    2.1.1 高斯過(guò)程回歸

    高斯過(guò)程回歸(GPR)模型預(yù)測(cè)步驟如圖5所示。高斯過(guò)程屬于一種特殊的二階矩隨機(jī)過(guò)程。通常認(rèn)為,自變量和因變量在有限維度下的概率分布都滿足一維或多維聯(lián)合高斯正態(tài)分布。多維高斯聯(lián)合分布由它的均值向量E(x)和協(xié)方差矩陣cov(x)所決定[21]。GP模型是服從聯(lián)合高斯分布的有限個(gè)隨機(jī)變量的集合,確定均值函數(shù)和協(xié)方差矩陣就能完整地確定一個(gè)高斯過(guò)程的性質(zhì)。協(xié)方差函數(shù)作為衡量數(shù)據(jù)點(diǎn)相互影響的標(biāo)準(zhǔn),是統(tǒng)計(jì)分析中的重要指標(biāo)之一,該函數(shù)表示兩個(gè)變量之間的協(xié)調(diào)變化率。常用協(xié)方差函數(shù)包括平方指數(shù)、Matern、各項(xiàng)同性指數(shù)等。該尺度的存在決定了具有一定差異的輸入特征是否被視為近似特征。此類超參數(shù)同樣被貝葉斯優(yōu)化器所考慮,并進(jìn)行優(yōu)化組合。

    2.1.2 人工神經(jīng)網(wǎng)絡(luò)

    人工神經(jīng)網(wǎng)絡(luò)(ANN)是一種模擬生物神經(jīng)系統(tǒng)的深度學(xué)習(xí)算法。算法的每個(gè)處理節(jié)點(diǎn)被稱為神經(jīng)元或節(jié)點(diǎn)。單個(gè)節(jié)點(diǎn)可以接收多個(gè)信息的輸入,但只返回一個(gè)輸出信息。從輸入到輸出的過(guò)程需要引入激活函數(shù)來(lái)提高神經(jīng)網(wǎng)絡(luò)的非線性表達(dá)能力。在神經(jīng)網(wǎng)絡(luò)中,每個(gè)神經(jīng)元都被賦予一個(gè)權(quán)重值,權(quán)重值的變化同樣影響著神經(jīng)網(wǎng)絡(luò)的性能,通過(guò)迭代訓(xùn)練樣本可以不斷調(diào)整神經(jīng)元的權(quán)重值,當(dāng)預(yù)測(cè)值與實(shí)際值之間的誤差達(dá)到最小值,則迭代停止。

    輸入層和輸出層之間的層稱為隱藏層,圖6即為單層隱藏層神經(jīng)網(wǎng)絡(luò)示意圖。除此之外,針對(duì)不同的數(shù)據(jù)特征,神經(jīng)網(wǎng)絡(luò)中隱藏層與神經(jīng)元的最優(yōu)設(shè)計(jì)數(shù)并不一樣。這些難以確定的超參數(shù)除了采用排列組合的方法試錯(cuò)外,還可以通過(guò)超參數(shù)調(diào)優(yōu)方法快速找到相對(duì)最優(yōu)解。

    2.1.3 隨機(jī)森林

    隨機(jī)森林是基于決策樹(shù)模型的一種集成算法,決策樹(shù)是一種非參數(shù)規(guī)則的算法,將一組數(shù)據(jù)點(diǎn)的特征空間進(jìn)行劃分,使得每一小塊空間區(qū)域都具有相似的響應(yīng)值。決策樹(shù)中的內(nèi)部節(jié)點(diǎn)可以針對(duì)輸入數(shù)據(jù)的某一特征進(jìn)行判斷,當(dāng)局部節(jié)點(diǎn)無(wú)法再根據(jù)特征進(jìn)行分類時(shí),則生成葉節(jié)點(diǎn)代表測(cè)試后的輸出。位于決策樹(shù)最高處的根節(jié)點(diǎn)則代表了數(shù)據(jù)中最重要的特征。簡(jiǎn)單決策樹(shù)模型示意圖如圖7所示。

    單一的決策樹(shù)模型通常會(huì)存在過(guò)度擬合的情況,且具有較高的方差。為了緩解這個(gè)問(wèn)題,以決策樹(shù)的分類回歸思想為基礎(chǔ),衍生出很多集成算法。集成學(xué)習(xí)方法將多個(gè)訓(xùn)練好的基本學(xué)習(xí)器通過(guò)一定的策略進(jìn)行結(jié)合,最終形成性能可靠的強(qiáng)學(xué)習(xí)器來(lái)進(jìn)行預(yù)測(cè),集成模型框架如圖8所示。

    隨機(jī)森林模型就是以決策樹(shù)為基本學(xué)習(xí)器的一個(gè)集成學(xué)習(xí)模型。如圖9所示,它包含多個(gè)由Bagging集成學(xué)習(xí)技術(shù)訓(xùn)練得到的決策樹(shù)。最終的預(yù)測(cè)結(jié)果由眾多具有較大差異的決策樹(shù)的輸出結(jié)果共同決定。隨機(jī)森林中的學(xué)習(xí)器個(gè)數(shù)與模型的學(xué)習(xí)效率和學(xué)習(xí)效果直接相關(guān),而學(xué)習(xí)器數(shù)量變化對(duì)兩者的影響呈現(xiàn)相反的表現(xiàn)趨勢(shì),此外,葉子節(jié)點(diǎn)上應(yīng)有的最少樣例數(shù)會(huì)決定模型是否更容易遭受噪聲數(shù)據(jù)的影響。諸如此類難以確定又對(duì)模型質(zhì)量有著重要影響的超參數(shù)將由貝葉斯優(yōu)化自動(dòng)確定,以使得模型性能達(dá)到其能力范圍內(nèi)的較高水準(zhǔn)。

    2.2 中國(guó)規(guī)范

    在中國(guó)《高強(qiáng)鋼結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)》(JGJ/T 483—2020)[5]中,軸心受壓構(gòu)件的極限承載力可由式(1)計(jì)算。

    式中:φ為穩(wěn)定系數(shù)。

    穩(wěn)定系數(shù)φ可根據(jù)截面類型、長(zhǎng)細(xì)比和鋼材屈服強(qiáng)度查詢規(guī)范中附錄表獲得。在最新的高強(qiáng)鋼設(shè)計(jì)規(guī)范中,箱型截面的高強(qiáng)鋼構(gòu)件在板厚小于40 mm時(shí),統(tǒng)一按b類截面設(shè)計(jì);當(dāng)板厚不小于40 mm時(shí),若板件寬厚比大于20,按b類截面設(shè)計(jì),否則按c類截面設(shè)計(jì)。截面分類方式的變化也是高強(qiáng)鋼設(shè)計(jì)規(guī)范與普通鋼設(shè)計(jì)規(guī)范中的主要區(qū)別。

    3 模型性能評(píng)價(jià)指標(biāo)

    為了研究開(kāi)發(fā)的ML模型的預(yù)測(cè)能力,將使用常見(jiàn)的不同性能指標(biāo)評(píng)估模型的性能,包括決定系數(shù)R2、均方誤差MSE、均方根誤差RMSE和平均絕對(duì)誤差MAE。

    1)決定系數(shù)

    4 模型驗(yàn)證與分析

    4.1 模型預(yù)測(cè)結(jié)果的對(duì)比

    選取4個(gè)評(píng)價(jià)指標(biāo),評(píng)估中國(guó)現(xiàn)有規(guī)范中的經(jīng)驗(yàn)?zāi)P秃?種常見(jiàn)機(jī)器學(xué)習(xí)模型對(duì)高強(qiáng)鋼焊接等截面箱型柱整體穩(wěn)定性的預(yù)測(cè)性能,并將相關(guān)數(shù)據(jù)進(jìn)行匯總。然而,模型在訓(xùn)練時(shí)的表現(xiàn)并不能代表其在模型測(cè)試時(shí)的表現(xiàn),需要加入模型在測(cè)試集上的表現(xiàn)來(lái)對(duì)模型進(jìn)行全面的評(píng)判。因此,采集了約70組中國(guó)學(xué)者對(duì)高強(qiáng)鋼焊接截面柱整體穩(wěn)定性進(jìn)行研究的真實(shí)試驗(yàn)數(shù)據(jù)[2,22-24],整理了其中25組相關(guān)數(shù)據(jù),見(jiàn)表2,并增加了模型對(duì)新的測(cè)試集數(shù)據(jù)的預(yù)測(cè)情況來(lái)進(jìn)一步驗(yàn)證模型性能與數(shù)據(jù)集質(zhì)量的優(yōu)劣,這些用于測(cè)試的數(shù)據(jù)需要保證未在模型的訓(xùn)練過(guò)程中出現(xiàn)。

    表3同時(shí)記錄了模型訓(xùn)練時(shí)的評(píng)價(jià)指標(biāo)情況以及機(jī)器學(xué)習(xí)模型或現(xiàn)有規(guī)范準(zhǔn)則模型在測(cè)試集數(shù)據(jù)上的表現(xiàn)?!陡邚?qiáng)鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范》(JGJ/T 483—2020)[5]所用模型的決定系數(shù)為0.83,RMSE值達(dá)到了1 176.6 kN,MAE也高達(dá)861.57 kN。模型預(yù)測(cè)精度在±861.6 kN,超出了測(cè)試集數(shù)據(jù)中平均極限承載力(3 149.4 kN)的27.3%。隨著鋼材免屈服強(qiáng)度的增大,設(shè)計(jì)規(guī)范模型逐漸表現(xiàn)出不適用性,且預(yù)測(cè)結(jié)果偏大。通過(guò)查閱相關(guān)文獻(xiàn),可以對(duì)此現(xiàn)象做出解釋,部分學(xué)者在對(duì)高強(qiáng)鋼焊接截面柱整體穩(wěn)定性能的研究中得出以下結(jié)論:采用中國(guó)規(guī)范設(shè)計(jì)時(shí),對(duì)于寬厚比小于20的Q690箱形柱,可選取a類柱子曲線[23];Q890和Q960鋼材焊接箱形軸壓構(gòu)件板厚小于40 mm時(shí)應(yīng)按a類曲線進(jìn)行設(shè)計(jì),區(qū)別于規(guī)范中建議使用b類曲線計(jì)算[24]。以上結(jié)論說(shuō)明中國(guó)現(xiàn)有規(guī)范對(duì)于高強(qiáng)鋼軸心軸壓構(gòu)件的截面分類并不完全準(zhǔn)確。且規(guī)范公式的提出考慮了設(shè)計(jì)師們進(jìn)行設(shè)計(jì)時(shí)的便捷性,通常是有過(guò)簡(jiǎn)化處理的。上述原因?qū)е卢F(xiàn)行規(guī)范在預(yù)測(cè)時(shí)的缺陷。

    反觀機(jī)器學(xué)習(xí)預(yù)測(cè)模型的表現(xiàn),機(jī)器學(xué)習(xí)預(yù)測(cè)模型中的高斯過(guò)程回歸模型整體表現(xiàn)明顯要優(yōu)于現(xiàn)有規(guī)范準(zhǔn)則模型。其在訓(xùn)練過(guò)程中的RMSE值為41.21 kN,MAE值為23.59 kN,決定系數(shù)為1.0。而在對(duì)測(cè)試集數(shù)據(jù)預(yù)測(cè)過(guò)程中,其RMSE值為265.13 kN,MAE值為163.17 kN,決定系數(shù)高達(dá)0.99,平均誤差率僅有5.6%。

    由表3可見(jiàn),機(jī)器學(xué)習(xí)模型對(duì)高強(qiáng)鋼焊接等截面箱型柱的整體穩(wěn)定性的預(yù)測(cè)效果相較現(xiàn)有規(guī)范有明顯的提升。但隨機(jī)森林模型在此類問(wèn)題的預(yù)測(cè)上表現(xiàn)很差,在訓(xùn)練過(guò)程中的決定系數(shù)為0.93,在測(cè)試集數(shù)據(jù)預(yù)測(cè)中決定系數(shù)僅有0.75。以決策樹(shù)為基本學(xué)習(xí)器的集成算法預(yù)測(cè)效果較差的原因可能是決策樹(shù)模型的決策邊界可能并不準(zhǔn)確且決策樹(shù)算法對(duì)個(gè)別數(shù)據(jù)比較敏感,進(jìn)而導(dǎo)致其集成模型的效果在回歸預(yù)測(cè)上的表現(xiàn)很差。

    對(duì)于機(jī)器學(xué)習(xí)預(yù)測(cè)模型在測(cè)試集上的表現(xiàn)相對(duì)模型訓(xùn)練時(shí)普遍變差的問(wèn)題,機(jī)器學(xué)習(xí)領(lǐng)域?qū)<艺J(rèn)為這是由于機(jī)器學(xué)習(xí)模型泛化能力不佳導(dǎo)致的,可通過(guò)調(diào)整超參數(shù)、使用K折交叉驗(yàn)證、歸一化數(shù)據(jù)等手段降低其影響,筆者研究中也有所涉及。其他的數(shù)據(jù)集本身的原因包括:1)纖維模型構(gòu)建的數(shù)據(jù)庫(kù)樣本與實(shí)際工程情況本就有著10%左右的誤差;2)在數(shù)據(jù)庫(kù)的構(gòu)建過(guò)程中,初始缺陷考慮得并不全面;3)實(shí)際試驗(yàn)數(shù)據(jù)的部分參數(shù)超出了數(shù)據(jù)庫(kù)樣本參數(shù)大小范圍約10%。能否對(duì)超出訓(xùn)練數(shù)據(jù)范圍的數(shù)據(jù)進(jìn)行準(zhǔn)確預(yù)測(cè)也是模型泛化能力強(qiáng)弱的體現(xiàn)。以上原因?qū)е铝藱C(jī)器學(xué)習(xí)模型在預(yù)測(cè)測(cè)試集數(shù)據(jù)時(shí)的性能下降。在后續(xù)研究中,將進(jìn)一步關(guān)注數(shù)據(jù)集和模型的質(zhì)量來(lái)提高機(jī)器模型在測(cè)試集上進(jìn)行預(yù)測(cè)時(shí)的表現(xiàn)。盡管如此,以高斯過(guò)程回歸模型為代表的眾多機(jī)器學(xué)習(xí)模型仍然展現(xiàn)出了優(yōu)于現(xiàn)有規(guī)范準(zhǔn)則模型及其他種類預(yù)測(cè)模型的性能。

    圖10為數(shù)據(jù)集在3個(gè)機(jī)器學(xué)習(xí)模型的中的真實(shí)值與預(yù)測(cè)值響應(yīng)圖,其中,藍(lán)色的點(diǎn)表示訓(xùn)練集中已知的輸出參數(shù),即軸心受壓極限承載力真實(shí)值。而黃色的點(diǎn)表示模型經(jīng)過(guò)訓(xùn)練后根據(jù)輸入?yún)?shù)得出的預(yù)測(cè)響應(yīng)值,即軸心受壓極限承載力預(yù)測(cè)值,當(dāng)兩點(diǎn)完全重合時(shí),表示機(jī)器學(xué)習(xí)模型在該點(diǎn)處的預(yù)測(cè)精度很高,重合的點(diǎn)越多,機(jī)器學(xué)習(xí)模型性能越好。另一種能直觀展現(xiàn)出機(jī)器學(xué)習(xí)模型預(yù)測(cè)效果的圖是圖11所展示的殘差圖,當(dāng)模型的性能足夠優(yōu)秀時(shí)其預(yù)測(cè)響應(yīng)與真實(shí)值相差無(wú)幾,在殘差圖中表現(xiàn)為所有預(yù)測(cè)點(diǎn)都接近于零線,從任意點(diǎn)到零線的垂直距離即為該點(diǎn)的殘差。

    綜合所有ML模型在模型性能可視化圖中的表現(xiàn)可以看出,高斯過(guò)程回歸模型和神經(jīng)網(wǎng)絡(luò)模型表現(xiàn)是這3種模型中的較優(yōu)模型。在預(yù)測(cè)響應(yīng)圖中,兩種模型的預(yù)測(cè)點(diǎn)與真實(shí)點(diǎn)重合度很高;在殘差圖中,兩種模型中各樣本點(diǎn)的殘差絕大多數(shù)都在100 kN以內(nèi),不超過(guò)訓(xùn)練樣本平均極限承載力(1 780.75 kN)的6%。說(shuō)明高斯過(guò)程回歸模型和神經(jīng)網(wǎng)絡(luò)模型在訓(xùn)練過(guò)程中的優(yōu)異性能。但結(jié)合表3可以看出,在預(yù)測(cè)測(cè)試集中新的實(shí)際工程數(shù)據(jù)時(shí),高斯過(guò)程回歸模型的穩(wěn)定性要略高于神經(jīng)網(wǎng)絡(luò)模型。因此,在高強(qiáng)鋼焊接箱形截面柱整體穩(wěn)定性的研究中,高斯過(guò)程回歸會(huì)是常見(jiàn)機(jī)器學(xué)習(xí)模型中的最優(yōu)模型,也是優(yōu)于現(xiàn)有規(guī)范經(jīng)驗(yàn)?zāi)P偷倪x擇。值得一提的是,高斯過(guò)程回歸模型預(yù)測(cè)25組真實(shí)實(shí)驗(yàn)數(shù)據(jù)時(shí),僅耗費(fèi)0.006 s,這樣的計(jì)算效率顯然要優(yōu)于更加精細(xì)的有限元模擬方法,盡管有限元建模能夠使模型的計(jì)算結(jié)果達(dá)到很高的精度,但其繁瑣的建模過(guò)程和較長(zhǎng)的計(jì)算時(shí)間使得效率更高的機(jī)器學(xué)習(xí)模型成為一種新的選擇。

    4.2 高斯過(guò)程回歸模型的可解釋性研究

    模型可解釋性是指對(duì)模型工作機(jī)制以及對(duì)模型預(yù)測(cè)結(jié)果的理解,機(jī)器學(xué)習(xí)模型的可解釋性越高,模型的可信任度就越高。諸如支持向量回歸等機(jī)器學(xué)習(xí)模型被稱為“黑盒”模型,從這些ML模型中探求輸入?yún)?shù)和輸出參數(shù)之的力學(xué)關(guān)系是很困難的。近年來(lái),陸續(xù)有學(xué)者研究并提出了幾種方法來(lái)提高M(jìn)L模型的可解釋性。部分依賴圖(PDP)和個(gè)體條件期望圖(ICE)就是用于對(duì)ML模型進(jìn)行解釋的常用可視化工具。

    部分依賴圖反映了在一個(gè)訓(xùn)練好的ML模型中輸入特征參數(shù)與輸出響應(yīng)參數(shù)之間的關(guān)系[25],通常包括線性關(guān)系、單調(diào)關(guān)系或者其他更復(fù)雜的關(guān)系。某個(gè)選定輸入特征參數(shù)的部分依賴性定義為忽略其他輸入變量的影響而獲得的平均預(yù)測(cè)值,或者說(shuō)部分相關(guān)性作為所選輸入特征參數(shù)的函數(shù),顯示了所選輸入特征參數(shù)對(duì)數(shù)據(jù)集的平均影響。各輸入特征的PDP圖如圖12所示。

    個(gè)體條件期望圖作為部分依賴圖的一種變體,代表了每個(gè)樣本的輸入特征參數(shù)和輸出響應(yīng)值之間的關(guān)系。部分相關(guān)性顯示了輸入和輸出之間的平均關(guān)系,但I(xiàn)CE圖將平均關(guān)系信息進(jìn)行了分解,并顯示了每次觀測(cè)時(shí)個(gè)體特征的依賴性[26]。

    選擇在上述研究中表現(xiàn)最好的高斯過(guò)程回歸模型,計(jì)算并繪制該模型的PDP圖,以確定8個(gè)輸入特征參數(shù)與高強(qiáng)鋼焊接等截面箱型柱極限承載力之間的關(guān)系。從圖12中的PDP圖可以看出,構(gòu)件的總高(總寬)、翼緣(腹板)厚度、材料屈服強(qiáng)度、材料彈性模量與高強(qiáng)鋼焊接等截面箱型柱的軸壓極限承載力在一定范圍內(nèi)呈正相關(guān)。構(gòu)件計(jì)算長(zhǎng)度、初始缺陷值、正則化長(zhǎng)細(xì)比、寬厚比與構(gòu)件軸壓極限承載力在一定范圍內(nèi)呈負(fù)相關(guān)。在對(duì)《高強(qiáng)鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范》(JGJ/T 483—2020)[5]中關(guān)于鋼結(jié)構(gòu)軸壓構(gòu)件整體穩(wěn)定性的相關(guān)規(guī)范公式進(jìn)行分析后發(fā)現(xiàn),機(jī)器學(xué)習(xí)模型做出的有關(guān)輸入?yún)?shù)對(duì)輸出參數(shù)影響趨勢(shì)的解釋均是合理的,這從側(cè)面體現(xiàn)了機(jī)器學(xué)習(xí)模型內(nèi)部運(yùn)作的合理性。但大部分趨勢(shì)都展現(xiàn)出線性或近線性,這與實(shí)際情況有所差距,這主要是由于訓(xùn)練樣本數(shù)據(jù)中的輸入?yún)?shù)范圍有限,如果將數(shù)據(jù)集參數(shù)范圍擴(kuò)大,將會(huì)得到更加豐富且貼合實(shí)際的部分依賴圖。在正則化長(zhǎng)細(xì)比部分依賴圖中可以看到更為明顯的趨勢(shì)特征,當(dāng)λn小于1時(shí),構(gòu)件的極限承載力以較大的梯度下降,當(dāng)λn到達(dá)1時(shí),構(gòu)件承載力下降速度大大減緩而進(jìn)入平臺(tái)期,當(dāng)λn大于1.5后,構(gòu)件承載力又恢復(fù)到λn數(shù)值較小時(shí)的下降梯度。

    4.3 神經(jīng)網(wǎng)絡(luò)模型的特征重要性研究

    使用Shapley Value方法制作特征參數(shù)重要性可視化圖像,一個(gè)特征的Shapley Value是該特征在所有的特征序列中的平均邊際貢獻(xiàn)。使用效果較好的神經(jīng)網(wǎng)絡(luò)模型計(jì)算8個(gè)輸入特征參數(shù)的特征重要性,其可視化的表現(xiàn)見(jiàn)圖13。然而,神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)結(jié)果會(huì)受其超參數(shù)的影響,因此,獲得的特征重要性排序并非絕對(duì),只作為建議來(lái)指導(dǎo)進(jìn)行工程實(shí)際。

    特征重要性分析表明,相對(duì)于其他輸入?yún)?shù),構(gòu)件的正則化長(zhǎng)細(xì)比、總高(總寬)、腹板厚度(翼緣厚度)、屈服強(qiáng)度是影響高強(qiáng)鋼焊接等截面箱型柱軸心受壓承載力的主要參數(shù)。由中國(guó)現(xiàn)有規(guī)范中對(duì)高強(qiáng)鋼構(gòu)件軸心受壓的整體穩(wěn)定性描述公式可知,高強(qiáng)鋼構(gòu)件整體穩(wěn)定性與構(gòu)件截面面積和材料屈服強(qiáng)度線性相關(guān),而特征重要性分析中得到的重要參數(shù)恰好就包含了屈服強(qiáng)度和面積(面積可由構(gòu)件總高與腹板厚度求得),這從雙向驗(yàn)證了規(guī)范與機(jī)器學(xué)習(xí)模型都具備較高的合理性。相比之下,構(gòu)件的初始缺陷和彈性模量對(duì)其極限承載力影響較小,這可能是由于數(shù)據(jù)庫(kù)中初始缺陷取值均僅為構(gòu)件計(jì)算長(zhǎng)度的千分之一以下,取值偏小導(dǎo)致初始缺陷對(duì)構(gòu)件極限承載力的影響不大。

    5 結(jié)論

    開(kāi)發(fā)了兩個(gè)單一模型和一個(gè)集成模型,用于預(yù)測(cè)高強(qiáng)鋼焊接等截面箱型柱的軸心受壓極限承載力以驗(yàn)證其整體穩(wěn)定性。將3種機(jī)器學(xué)習(xí)模型的預(yù)測(cè)性能和現(xiàn)有規(guī)范中的經(jīng)驗(yàn)?zāi)P瓦M(jìn)行比較,并對(duì)最優(yōu)機(jī)器學(xué)習(xí)模型進(jìn)行可解釋性研究。研究結(jié)果表明,機(jī)器學(xué)習(xí)技術(shù)在預(yù)測(cè)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性方面的可行性和高精度,并為將數(shù)據(jù)驅(qū)動(dòng)模型納入設(shè)計(jì)規(guī)范的修訂和完善提供了重要參考。得到以下主要結(jié)論:

    1)與其他ML模型和基于物理的方程相比,高斯過(guò)程回歸模型預(yù)測(cè)得到的結(jié)果最為準(zhǔn)確,該模型對(duì)測(cè)試集數(shù)據(jù)預(yù)測(cè)的RMSE值為265.13 kN,MAE值為163.17 kN,決定系數(shù)高達(dá)0.99。與25組已有試驗(yàn)結(jié)果相比,平均絕對(duì)誤差僅有5.6%。此外,隨著未來(lái)數(shù)據(jù)庫(kù)數(shù)據(jù)愈加豐富,參數(shù)數(shù)值包含范圍愈加廣泛,高斯過(guò)程回歸模型的性能可以進(jìn)一步提高。這項(xiàng)工作的結(jié)果證明了機(jī)器學(xué)習(xí)技術(shù)在預(yù)測(cè)高強(qiáng)鋼構(gòu)件整體穩(wěn)定性方面的可行性和高精度,并為將數(shù)據(jù)驅(qū)動(dòng)模型納入下一代國(guó)際設(shè)計(jì)規(guī)范做出了引導(dǎo)。

    2)根據(jù)模型的可解釋性研究發(fā)現(xiàn),作為最優(yōu)機(jī)器學(xué)習(xí)預(yù)測(cè)模型的高斯過(guò)程回歸模型,其中,輸入?yún)?shù)對(duì)預(yù)測(cè)結(jié)果的影響趨勢(shì)與現(xiàn)有規(guī)范準(zhǔn)則中的經(jīng)驗(yàn)公式一致,在證明了機(jī)器學(xué)習(xí)模型準(zhǔn)確性的基礎(chǔ)上又驗(yàn)證了模型的合理性。

    3)根據(jù)輸入特征參數(shù)重要性排序可知,正則化長(zhǎng)細(xì)比、材料屈服強(qiáng)度以及和構(gòu)件截面面積相關(guān)的幾何特征是對(duì)高強(qiáng)鋼軸心受壓承載力影響較大的參數(shù)。這與中國(guó)現(xiàn)有規(guī)范準(zhǔn)則中所述接近,而在機(jī)器學(xué)習(xí)模型的解釋中正則化長(zhǎng)細(xì)比作為對(duì)輸出結(jié)果影響最大的重要參數(shù)參與預(yù)測(cè)運(yùn)算。

    4)機(jī)器學(xué)習(xí)模型應(yīng)用于高強(qiáng)鋼構(gòu)件力學(xué)性能預(yù)測(cè)的工作中具有能夠快速甄別關(guān)鍵參數(shù)與力學(xué)性能之間的關(guān)系優(yōu)勢(shì),但現(xiàn)階段其計(jì)算精度相較于精細(xì)的有限元分析尚有差距,且計(jì)算結(jié)果的合理性依賴于模型的選取和大量數(shù)據(jù)的收集,因此,現(xiàn)階段成果僅作為高強(qiáng)鋼構(gòu)件力學(xué)性能分析的參考與借鑒。

    參考文獻(xiàn)

    [1]? 周祎, 仵振, 潘毅, 等. 高強(qiáng)鋼焊接箱形柱力學(xué)性能研究(Ⅱ): 二階非彈性分析[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2022, 43(7): 165-174.

    ZHOU Y, WU Z, PAN Y, et al. Research on mechanical properties of high strength steel welded box columns (Ⅱ): Second-order inelastic analysis [J]. Journal of Building Structures, 2022, 43(7): 165-174. (in Chinese)

    [2]? 薛加燁. 高強(qiáng)度鋼材受壓構(gòu)件整體穩(wěn)定性能試驗(yàn)研究 [D]. 南京: 東南大學(xué), 2014.

    XUE J Y. Experimental research on the overall buckling behavior of high strength steel members under compression [D]. Nanjing: Southeast University, 2014. (in Chinese)

    [3]? 邱林波, 侯兆新, 薛素鐸, 等. Q550高強(qiáng)度鋼焊接H形柱軸心受壓承載力試驗(yàn)研究[J]. 建筑鋼結(jié)構(gòu)進(jìn)展, 2015, 17(3): 7-12, 18.

    QIU L B, HOU Z X, XUE S D, et al. Experimental study on axial compression capacity of welding H-shaped Q550 high-strength steel columns [J]. Progress in Steel Building Structures, 2015, 17(3): 7-12, 18. (in Chinese)

    [4]? 趙金友, 李晶, 王鈞, 等. Q460高強(qiáng)鋼焊接工字形截面簡(jiǎn)支梁整體穩(wěn)定性能與設(shè)計(jì)方法研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2021, 42(11): 61-70.

    ZHAO J Y, LI J, WANG J, et al. Design method and overall buckling behavior of simply-supported beams with Q460 high strength steel welded I-section [J]. Journal of Building Structures, 2021, 42(11): 61-70. (in Chinese)

    [5]? 高強(qiáng)鋼結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn): JGJ/T 483—2020 [S]. 北京: 中國(guó)建筑工業(yè)出版社, 2020.

    Standard for design of high strength stell structure: JGJ/T 483—2020 [S]. Beijing: China Architecture & Building Press, 2020. (in Chinese)

    [6]? MENG X, GARDNER L. Stability and design of normal and high strength steel CHS beam-columns [J]. Engineering Structures, 2022, 251: 113361.

    [7]? LIMBACHIYA V, SHAMASS R. Application of artificial neural networks for web-post shear resistance of cellular steel beams [J]. Thin-Walled Structures, 2021, 161: 107414.

    [8]? WAKJIRA T G, IBRAHIM M, EBEAD U, et al. Explainable machine learning model and reliability analysis for flexural capacity prediction of RC beams strengthened in flexure with FRCM [J]. Engineering Structures, 2022, 255: 113903.

    [9]? SAROTHI S Z, AHMED K S, KHAN N I, et al. Predicting bearing capacity of double shear bolted connections using machine learning [J]. Engineering Structures, 2022, 251: 113497.

    [10]? 薛紅新. 基于機(jī)器學(xué)習(xí)方法的分類與預(yù)測(cè)問(wèn)題研究[D]. 太原: 中北大學(xué), 2019.

    XUE H X. Research on the problems of classification and prediction based on machine learning methods [D]. Taiyuan: North University of China, 2019. (in Chinese)

    [11]? 陳驥. 鋼結(jié)構(gòu)穩(wěn)定理論與設(shè)計(jì)[M]. 2版. 北京: 科學(xué)出版社, 2003.

    CHEN J. Stability of steel structures theory and design [M]. 2nd edition. Beijing: Science Press, 2003. (in Chinese)

    [12]? LIN Y J, LIU K Q, LI T J, et al. Numerical study and design method of high strength steel welded box-columns [J]. Advanced Steel Construction, 2020, 17(2): 127-135.

    [13]? 班慧勇, 施剛, 石永久. 960MPa高強(qiáng)鋼焊接箱形截面殘余應(yīng)力試驗(yàn)及統(tǒng)一分布模型研究[J]. 土木工程學(xué)報(bào), 2013, 46(11): 63-69.

    BAN H Y, SHI G, SHI Y J. Experimental study on residual stress in 960MPa high strength steel welded box sections and unified model [J]. China Civil Engineering Journal, 2013, 46(11): 63-69. (in Chinese)

    [14]? 潘毅, 仵振, 周祎, 等. 高強(qiáng)鋼焊接箱形柱受力性能研究(Ⅰ): 殘余應(yīng)力統(tǒng)一分布模型[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2022, 43(3): 138-147.

    PAN Y, WU Z, ZHOU Y, et al. Study on mechanical properties of high strength steel welded box columns (Ⅰ): Unified model of residual stress [J]. Journal of Building Structures, 2022, 43(3): 138-147. (in Chinese)

    [15]? 羅婧文. 考慮初始缺陷和殘余應(yīng)力影響的高強(qiáng)鋼壓桿穩(wěn)定性分析與評(píng)價(jià) [D]. 成都:西南交通大學(xué), 2020.

    LUO J W. Stability analysis and evaluation of high-strength steel compression column considering the effect of initial defects and residual stress [D]. Chengdu: Southwest Jiaotong University, 2020. (in Chinese)

    [16]? TAHIR Z U R, MANDAL P, ADIL M T, et al. Application of artificial neural network to predict buckling load of thin cylindrical shells under axial compression [J]. Engineering Structures, 2021, 248: 113221.

    [17]? FERREIRA F P V, SHAMASS R, LIMBACHIYA V, et al. Lateral-torsional buckling resistance prediction model for steel cellular beams generated by artificial neural networks (ANN) [J]. Thin-Walled Structures, 2022, 170: 108592.

    [18]? 胡鑫. 基于人工神經(jīng)網(wǎng)絡(luò)的HPC強(qiáng)度預(yù)測(cè)[D]. 長(zhǎng)沙: 湖南大學(xué), 2014.

    HU X. Prediction of high performance concrete strength based on artificial neural network [D]. Changsha: Hunan University, 2014. (in Chinese)

    [19]? KANG M C, YOO D Y, GUPTA R, et al. Machine learning-based prediction for compressive and flexural strengths of steel fiber-reinforced concrete [J]. Construction and Building Materials, 2021, 266: 121117.

    [20]? DEGTYAREV V V. Neural networks for predicting shear strength of CFS channels with slotted webs [J]. Journal of Constructional Steel Research, 2021, 177: 106443.

    [21]? 張寅. 用于回歸預(yù)測(cè)的高斯過(guò)程模型研究[D]. 天津: 河北工業(yè)大學(xué), 2014.

    ZHANG Y. Research on Gaussian process for regression and prediction [D]. Tianjin: Hebei University of Technology, 2014. (in Chinese)

    [22]? 李國(guó)強(qiáng), 王彥博, 陳素文, 等. Q460高強(qiáng)鋼焊接箱形柱軸心受壓極限承載力參數(shù)分析[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2011, 32(11): 149-155.

    LI G Q, WANG Y B, CHEN S W, et al. Parametric analysis of ultimate bearing capacity of Q460 high strength steel welded box columns under axial compression [J]. Journal of Building Structures, 2011, 32(11): 149-155. (in Chinese)

    [23]? 李國(guó)強(qiáng), 李天際, 王彥博. Q690高強(qiáng)鋼焊接箱形軸壓構(gòu)件整體穩(wěn)定研究及設(shè)計(jì)建議[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2017, 38(10): 1-9.

    LI G Q, LI T J, WANG Y B. Overall buckling behavior and design of Q690 high-strength steel welded box-columns subjected to axial compression [J]. Journal of Building Structures, 2017, 38(10): 1-9. (in Chinese)

    [24]? 班慧勇, 施剛, 石永久. 高強(qiáng)鋼焊接箱形軸壓構(gòu)件整體穩(wěn)定設(shè)計(jì)方法研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2014, 35(5): 57-64.

    BAN H Y, SHI G, SHI Y J. Research on design method for overall buckling behavior of welded box columns fabricated from high-strength steels [J]. Journal of Building Structures, 2014, 35(5): 57-64. (in Chinese)

    [25]? FRIEDMAN J H. Greedy function approximation: A gradient boosting machine [J]. The Annals of Statistics, 2001, 29(5): 1189-1232.

    [26]? GOLDSTEIN A, KAPELNER A, BLEICH J, et al. Peeking inside the black box: Visualizing statistical learning with plots of individual conditional expectation [J]. Journal of Computational and Graphical Statistics, 2015, 24(1): 44-65.

    (編輯? 胡玲)

    收稿日期:2022?04?30

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(52278229)

    作者簡(jiǎn)介:張營(yíng)營(yíng)(1985- ),男,博士,教授,主要從事大跨空間結(jié)構(gòu)研究,E-mail:zhangyingying85@163.com。

    通信作者:周祎(通信作者),男,博士,E-mail:suzhouzhouyi@home.swjtu.edu.cn。

    Received: 2022?04?30

    Foundation item: National Natural Science Foundation of China (No. 52278229)

    Author brief: ZHANG Yingying (1985- ), PhD, associate professor, main research interest: long span spatial structure, E-mail: zhangyingying85@163.com.

    corresponding author:ZHOU Yi (corresponding author), PhD, E-mail: suzhouzhouyi@home.swjtu.edu.cn.

    猜你喜歡
    預(yù)測(cè)模型機(jī)器學(xué)習(xí)
    基于矩陣?yán)碚撓碌母咝=處熑藛T流動(dòng)趨勢(shì)預(yù)測(cè)
    東方教育(2016年9期)2017-01-17 21:04:14
    基于支持向量回歸的臺(tái)灣旅游短期客流量預(yù)測(cè)模型研究
    基于詞典與機(jī)器學(xué)習(xí)的中文微博情感分析
    基于神經(jīng)網(wǎng)絡(luò)的北京市房?jī)r(jià)預(yù)測(cè)研究
    商情(2016年43期)2016-12-23 14:23:13
    中國(guó)石化J分公司油氣開(kāi)發(fā)投資分析與預(yù)測(cè)模型研究
    基于機(jī)器學(xué)習(xí)的圖像特征提取技術(shù)在圖像版權(quán)保護(hù)中的應(yīng)用
    基于網(wǎng)絡(luò)搜索數(shù)據(jù)的平遙旅游客流量預(yù)測(cè)分析
    前綴字母為特征在維吾爾語(yǔ)文本情感分類中的研究
    基于支持向量機(jī)的金融數(shù)據(jù)分析研究
    基于IOWHA法的物流需求組合改善與預(yù)測(cè)模型構(gòu)建
    免费电影在线观看免费观看| 床上黄色一级片| 狂野欧美激情性xxxx| or卡值多少钱| 国产成人福利小说| 老汉色av国产亚洲站长工具| 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩卡通动漫| 老司机午夜十八禁免费视频| 国产欧美日韩精品一区二区| 欧美日韩黄片免| 热99在线观看视频| 制服丝袜大香蕉在线| 精品一区二区三区视频在线 | 成人永久免费在线观看视频| 成人三级黄色视频| 亚洲七黄色美女视频| 在线观看免费午夜福利视频| 国产亚洲欧美98| 日本撒尿小便嘘嘘汇集6| 99久久精品国产亚洲精品| 午夜精品久久久久久毛片777| 非洲黑人性xxxx精品又粗又长| 热99re8久久精品国产| www.www免费av| 99riav亚洲国产免费| 国产成人影院久久av| 久99久视频精品免费| 国模一区二区三区四区视频 | 欧美成人性av电影在线观看| 97碰自拍视频| 亚洲美女黄片视频| 亚洲av成人不卡在线观看播放网| 国产精品电影一区二区三区| 精品日产1卡2卡| 中文在线观看免费www的网站| 91在线观看av| 97碰自拍视频| 91在线精品国自产拍蜜月 | 好看av亚洲va欧美ⅴa在| 91字幕亚洲| 国产av不卡久久| 99久久无色码亚洲精品果冻| 午夜福利视频1000在线观看| 91久久精品国产一区二区成人 | 日韩精品中文字幕看吧| 久久久久国产精品人妻aⅴ院| 精品国产乱子伦一区二区三区| 国产成人一区二区三区免费视频网站| 后天国语完整版免费观看| 欧美日韩一级在线毛片| 国产成人av教育| 成人性生交大片免费视频hd| 99在线视频只有这里精品首页| 俺也久久电影网| 久久九九热精品免费| 中国美女看黄片| 日日干狠狠操夜夜爽| 久久久久久久精品吃奶| 狠狠狠狠99中文字幕| av视频在线观看入口| 欧美xxxx黑人xx丫x性爽| 欧美成狂野欧美在线观看| 99久久综合精品五月天人人| av福利片在线观看| 国产精品一区二区三区四区久久| 狂野欧美白嫩少妇大欣赏| 亚洲av中文字字幕乱码综合| 91av网一区二区| 欧美日韩乱码在线| 成人18禁在线播放| 1000部很黄的大片| 成在线人永久免费视频| 精品99又大又爽又粗少妇毛片 | 久久精品91无色码中文字幕| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 国产高清激情床上av| 亚洲精华国产精华精| 欧美国产日韩亚洲一区| 国产一级毛片七仙女欲春2| 长腿黑丝高跟| 色视频www国产| 亚洲av成人一区二区三| 一级a爱片免费观看的视频| 成在线人永久免费视频| 欧美乱妇无乱码| 亚洲人成伊人成综合网2020| 97超级碰碰碰精品色视频在线观看| 成人av在线播放网站| 国产成人欧美在线观看| 国产高清videossex| 国产免费男女视频| 成年版毛片免费区| 91麻豆精品激情在线观看国产| 又大又爽又粗| 免费一级毛片在线播放高清视频| av天堂在线播放| a级毛片a级免费在线| 色视频www国产| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 久久九九热精品免费| 亚洲 欧美一区二区三区| 日韩免费av在线播放| 国产精品九九99| 麻豆成人午夜福利视频| 亚洲一区二区三区不卡视频| 他把我摸到了高潮在线观看| 真人一进一出gif抽搐免费| 啦啦啦韩国在线观看视频| 国内少妇人妻偷人精品xxx网站 | 日本成人三级电影网站| 老汉色av国产亚洲站长工具| 精品久久蜜臀av无| 国产视频一区二区在线看| 亚洲av免费在线观看| 女警被强在线播放| 国产毛片a区久久久久| 久久婷婷人人爽人人干人人爱| 国产成人影院久久av| 看免费av毛片| 美女午夜性视频免费| 亚洲国产欧美一区二区综合| 欧美一级a爱片免费观看看| 国产精品综合久久久久久久免费| 黄频高清免费视频| 天堂网av新在线| 老熟妇仑乱视频hdxx| 久久久精品欧美日韩精品| 亚洲欧美精品综合久久99| 搡老熟女国产l中国老女人| 人妻夜夜爽99麻豆av| 国产三级中文精品| 久久国产乱子伦精品免费另类| 日韩欧美三级三区| 精品无人区乱码1区二区| 一个人看视频在线观看www免费 | 亚洲精品久久国产高清桃花| 成人欧美大片| 国产真人三级小视频在线观看| cao死你这个sao货| 久久午夜综合久久蜜桃| 精品人妻1区二区| cao死你这个sao货| 国产1区2区3区精品| 亚洲美女黄片视频| 久久九九热精品免费| 午夜日韩欧美国产| 国内揄拍国产精品人妻在线| 久99久视频精品免费| 免费在线观看视频国产中文字幕亚洲| 亚洲精品色激情综合| 精品欧美国产一区二区三| 人人妻人人看人人澡| 精华霜和精华液先用哪个| 亚洲国产欧美一区二区综合| 中文字幕久久专区| 国产一区在线观看成人免费| 国产成人影院久久av| 亚洲中文日韩欧美视频| 在线永久观看黄色视频| 99热6这里只有精品| 757午夜福利合集在线观看| 精品乱码久久久久久99久播| www国产在线视频色| 熟女电影av网| 精品国产乱码久久久久久男人| 国产精品久久电影中文字幕| 亚洲电影在线观看av| 久久久久久久久久黄片| 男女视频在线观看网站免费| 丰满的人妻完整版| 国产精品精品国产色婷婷| 99久国产av精品| 黄色日韩在线| 国产精品自产拍在线观看55亚洲| 观看免费一级毛片| 美女午夜性视频免费| 欧美乱妇无乱码| 久久久久免费精品人妻一区二区| 日韩欧美三级三区| 色视频www国产| 女同久久另类99精品国产91| 国产欧美日韩精品亚洲av| 母亲3免费完整高清在线观看| 国产欧美日韩一区二区精品| 欧美一级毛片孕妇| www.熟女人妻精品国产| 亚洲在线观看片| 在线观看66精品国产| 日韩欧美精品v在线| 天天添夜夜摸| 国产成人av激情在线播放| 亚洲真实伦在线观看| 老鸭窝网址在线观看| 免费在线观看日本一区| 日韩成人在线观看一区二区三区| 两性夫妻黄色片| 在线观看舔阴道视频| 久久99热这里只有精品18| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 噜噜噜噜噜久久久久久91| 日韩欧美在线二视频| 国产高潮美女av| 两个人视频免费观看高清| 精品久久久久久久久久免费视频| 国产一区在线观看成人免费| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| xxx96com| av国产免费在线观看| 最近在线观看免费完整版| 欧美日韩国产亚洲二区| 天堂动漫精品| av黄色大香蕉| 美女免费视频网站| 国产成人一区二区三区免费视频网站| 国产极品精品免费视频能看的| 国内毛片毛片毛片毛片毛片| 国产激情久久老熟女| 午夜两性在线视频| 一二三四社区在线视频社区8| 免费大片18禁| 国产真实乱freesex| 少妇的逼水好多| 色综合亚洲欧美另类图片| 很黄的视频免费| 日本成人三级电影网站| 91在线精品国自产拍蜜月 | 成年女人看的毛片在线观看| 亚洲九九香蕉| 久久久久精品国产欧美久久久| 男人舔奶头视频| 99riav亚洲国产免费| 看免费av毛片| 日韩欧美一区二区三区在线观看| 男插女下体视频免费在线播放| 中文资源天堂在线| 国产黄a三级三级三级人| 国产v大片淫在线免费观看| 我的老师免费观看完整版| 亚洲黑人精品在线| 亚洲一区二区三区不卡视频| 淫秽高清视频在线观看| 中文字幕久久专区| 色综合亚洲欧美另类图片| 首页视频小说图片口味搜索| 91字幕亚洲| 一边摸一边抽搐一进一小说| 日韩精品中文字幕看吧| a级毛片a级免费在线| 日本一本二区三区精品| 美女免费视频网站| 首页视频小说图片口味搜索| 精品人妻1区二区| 成年人黄色毛片网站| 老熟妇乱子伦视频在线观看| 欧美日韩乱码在线| 亚洲在线观看片| 悠悠久久av| 99国产极品粉嫩在线观看| 女生性感内裤真人,穿戴方法视频| 桃红色精品国产亚洲av| 午夜精品一区二区三区免费看| 国产精品综合久久久久久久免费| 舔av片在线| 蜜桃久久精品国产亚洲av| 小说图片视频综合网站| 日本a在线网址| 中文字幕高清在线视频| 一二三四社区在线视频社区8| 亚洲精品乱码久久久v下载方式 | 我要搜黄色片| 中文字幕高清在线视频| 免费高清视频大片| www.熟女人妻精品国产| 亚洲美女视频黄频| 性欧美人与动物交配| 日韩欧美 国产精品| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av| 观看美女的网站| 一级毛片精品| 九色成人免费人妻av| svipshipincom国产片| 亚洲美女黄片视频| 最新在线观看一区二区三区| 午夜两性在线视频| 亚洲中文av在线| 操出白浆在线播放| www.熟女人妻精品国产| 999精品在线视频| 久久久久久大精品| 国产视频内射| 精品国产超薄肉色丝袜足j| 男女视频在线观看网站免费| 成年版毛片免费区| 99久久无色码亚洲精品果冻| 老熟妇乱子伦视频在线观看| 两个人视频免费观看高清| 精品久久久久久久人妻蜜臀av| 欧美av亚洲av综合av国产av| 99riav亚洲国产免费| 日本在线视频免费播放| 手机成人av网站| 久久这里只有精品中国| 成人性生交大片免费视频hd| 久久久久久人人人人人| 精品一区二区三区视频在线 | 他把我摸到了高潮在线观看| 99久久成人亚洲精品观看| 美女午夜性视频免费| 色视频www国产| 亚洲av成人不卡在线观看播放网| av女优亚洲男人天堂 | 老司机午夜十八禁免费视频| 神马国产精品三级电影在线观看| 亚洲人与动物交配视频| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 亚洲欧美精品综合一区二区三区| 天堂√8在线中文| 亚洲成a人片在线一区二区| 亚洲人成网站高清观看| 久久香蕉精品热| 19禁男女啪啪无遮挡网站| 成熟少妇高潮喷水视频| 在线观看一区二区三区| 欧美中文综合在线视频| 国产精品香港三级国产av潘金莲| 蜜桃久久精品国产亚洲av| 国产真人三级小视频在线观看| 国产精品亚洲av一区麻豆| 此物有八面人人有两片| 久久久久性生活片| 午夜影院日韩av| 老司机午夜十八禁免费视频| 日本与韩国留学比较| 老司机午夜十八禁免费视频| 嫩草影视91久久| 日本 欧美在线| 国产伦精品一区二区三区四那| 国产精品久久久久久久电影 | 中文字幕人妻丝袜一区二区| 美女高潮的动态| 亚洲中文字幕日韩| av在线天堂中文字幕| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 又大又爽又粗| 一区二区三区激情视频| 久久性视频一级片| 身体一侧抽搐| 国产精品美女特级片免费视频播放器 | 亚洲国产看品久久| 法律面前人人平等表现在哪些方面| 成人三级黄色视频| 亚洲男人的天堂狠狠| 精品久久久久久,| 桃红色精品国产亚洲av| 日韩大尺度精品在线看网址| 我的老师免费观看完整版| 免费大片18禁| 欧美中文综合在线视频| 不卡一级毛片| 国产精品精品国产色婷婷| 亚洲天堂国产精品一区在线| 久久精品人妻少妇| 日本 欧美在线| 亚洲欧洲精品一区二区精品久久久| 99re在线观看精品视频| 丰满人妻熟妇乱又伦精品不卡| 在线永久观看黄色视频| 亚洲av成人精品一区久久| 亚洲18禁久久av| 怎么达到女性高潮| 90打野战视频偷拍视频| 国产精品99久久久久久久久| 最近最新中文字幕大全免费视频| 亚洲av五月六月丁香网| 一本久久中文字幕| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| 国产精品电影一区二区三区| 欧美一区二区精品小视频在线| 九九热线精品视视频播放| 日本黄色片子视频| 午夜福利在线观看免费完整高清在 | 久久久久国内视频| 欧美日本视频| 在线视频色国产色| 久久久精品大字幕| 男人舔奶头视频| 国产精品久久久久久人妻精品电影| 久久国产精品人妻蜜桃| 亚洲成人免费电影在线观看| 99在线人妻在线中文字幕| 欧美极品一区二区三区四区| 黄色 视频免费看| 亚洲欧美日韩卡通动漫| 亚洲人成网站高清观看| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 日本免费一区二区三区高清不卡| 在线看三级毛片| 日本黄色视频三级网站网址| 成在线人永久免费视频| 国产av麻豆久久久久久久| 久久天堂一区二区三区四区| 国产在线精品亚洲第一网站| 日韩三级视频一区二区三区| 美女被艹到高潮喷水动态| 视频区欧美日本亚洲| 欧洲精品卡2卡3卡4卡5卡区| 欧美黑人欧美精品刺激| 在线国产一区二区在线| 黄频高清免费视频| 又紧又爽又黄一区二区| 国产真实乱freesex| 国产精品影院久久| 三级毛片av免费| 黑人巨大精品欧美一区二区mp4| 日本精品一区二区三区蜜桃| 波多野结衣高清作品| 在线观看免费午夜福利视频| 欧美绝顶高潮抽搐喷水| www.999成人在线观看| 国产av一区在线观看免费| 午夜福利欧美成人| aaaaa片日本免费| 夜夜看夜夜爽夜夜摸| 久久久精品欧美日韩精品| 亚洲精品乱码久久久v下载方式 | 欧美日韩综合久久久久久 | 69av精品久久久久久| 国产黄片美女视频| 美女cb高潮喷水在线观看 | 精品久久蜜臀av无| 亚洲美女视频黄频| 啦啦啦韩国在线观看视频| 中文字幕熟女人妻在线| 天天一区二区日本电影三级| 一边摸一边抽搐一进一小说| 日韩免费av在线播放| 欧美黄色淫秽网站| 草草在线视频免费看| 日本一本二区三区精品| 久久久精品欧美日韩精品| 岛国在线观看网站| 亚洲成av人片免费观看| 国产精品久久久久久亚洲av鲁大| 熟女少妇亚洲综合色aaa.| 特级一级黄色大片| 偷拍熟女少妇极品色| 白带黄色成豆腐渣| 国产欧美日韩一区二区精品| 色在线成人网| 一本精品99久久精品77| 成人特级av手机在线观看| 亚洲最大成人中文| 美女cb高潮喷水在线观看 | 欧美中文综合在线视频| 老汉色av国产亚洲站长工具| 国产伦一二天堂av在线观看| 久久久久国内视频| 嫩草影视91久久| 精品午夜福利视频在线观看一区| 99热只有精品国产| 色综合婷婷激情| cao死你这个sao货| 免费在线观看影片大全网站| 少妇的丰满在线观看| 色综合亚洲欧美另类图片| 日日夜夜操网爽| 俄罗斯特黄特色一大片| 黄色日韩在线| 国产美女午夜福利| 在线观看一区二区三区| 看免费av毛片| 亚洲七黄色美女视频| 亚洲熟女毛片儿| 三级毛片av免费| 国产不卡一卡二| 天堂网av新在线| 精品国产乱码久久久久久男人| av天堂中文字幕网| 国产精品美女特级片免费视频播放器 | 欧美乱妇无乱码| av中文乱码字幕在线| 日本在线视频免费播放| 亚洲av日韩精品久久久久久密| 成人高潮视频无遮挡免费网站| 又大又爽又粗| 国产av在哪里看| 国产成人啪精品午夜网站| 国产精品一区二区三区四区久久| 特级一级黄色大片| 亚洲成av人片免费观看| 亚洲国产精品sss在线观看| 观看免费一级毛片| 国产精品女同一区二区软件 | 国内精品久久久久久久电影| 精品国产三级普通话版| 亚洲av中文字字幕乱码综合| 夜夜爽天天搞| 国产爱豆传媒在线观看| 成人三级做爰电影| 变态另类丝袜制服| 亚洲午夜精品一区,二区,三区| 麻豆国产av国片精品| 亚洲自拍偷在线| 久久人妻av系列| 国产一区二区三区在线臀色熟女| h日本视频在线播放| 日韩有码中文字幕| 日韩 欧美 亚洲 中文字幕| 亚洲av成人精品一区久久| 91av网站免费观看| 成人特级黄色片久久久久久久| 国产视频内射| 精品乱码久久久久久99久播| 久久精品亚洲精品国产色婷小说| 俺也久久电影网| 最好的美女福利视频网| 欧美大码av| 亚洲av电影在线进入| h日本视频在线播放| 国产麻豆成人av免费视频| 在线观看一区二区三区| 看片在线看免费视频| 免费人成视频x8x8入口观看| 国产视频一区二区在线看| 毛片女人毛片| 国产精品综合久久久久久久免费| 国产一区二区三区视频了| 岛国在线免费视频观看| 搡老熟女国产l中国老女人| 香蕉丝袜av| 国产伦一二天堂av在线观看| 国产伦人伦偷精品视频| 国产伦一二天堂av在线观看| 最近最新中文字幕大全免费视频| 亚洲色图av天堂| 亚洲五月天丁香| 小说图片视频综合网站| 欧美成狂野欧美在线观看| 久久久久久九九精品二区国产| 美女 人体艺术 gogo| 久久这里只有精品中国| 黑人操中国人逼视频| 欧美成人一区二区免费高清观看 | 久久性视频一级片| 18禁黄网站禁片免费观看直播| 91老司机精品| 欧美zozozo另类| 久久久久久人人人人人| 成人18禁在线播放| 久久中文字幕一级| 麻豆久久精品国产亚洲av| 欧美一级毛片孕妇| 久久久久久国产a免费观看| av天堂中文字幕网| 色视频www国产| 热99在线观看视频| 国产在线精品亚洲第一网站| 亚洲精品国产精品久久久不卡| 成人亚洲精品av一区二区| 在线国产一区二区在线| 国产91精品成人一区二区三区| www.熟女人妻精品国产| 宅男免费午夜| 精品国产乱子伦一区二区三区| 成年女人永久免费观看视频| 欧美日韩一级在线毛片| 日本五十路高清| 很黄的视频免费| 国产一区二区在线观看日韩 | 国产高清videossex| 免费观看人在逋| 99久久久亚洲精品蜜臀av| 日韩欧美在线二视频| 国产精品精品国产色婷婷| 国内精品美女久久久久久| 国产精品久久电影中文字幕| 成在线人永久免费视频| 亚洲一区二区三区色噜噜| 国产精品久久电影中文字幕| 听说在线观看完整版免费高清| 亚洲中文字幕一区二区三区有码在线看 | 12—13女人毛片做爰片一| 美女扒开内裤让男人捅视频| 天天添夜夜摸| 人妻夜夜爽99麻豆av| 成人18禁在线播放| 午夜亚洲福利在线播放| 日韩人妻高清精品专区| xxx96com| 免费搜索国产男女视频| 日本与韩国留学比较| 窝窝影院91人妻| 午夜福利在线观看吧| 亚洲乱码一区二区免费版| 国产精品电影一区二区三区| 这个男人来自地球电影免费观看| 美女大奶头视频| 国产激情欧美一区二区| 中国美女看黄片| 亚洲成人中文字幕在线播放| 国产成人精品久久二区二区免费| a在线观看视频网站| 少妇丰满av| netflix在线观看网站|