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

    基于多保真度代理模型的混凝土泵車(chē)臂架銷(xiāo)軸輕量化設(shè)計(jì)

    2025-08-13 00:00:00李鵬李夢(mèng)聰肖立波王一棠宋學(xué)官陽(yáng)領(lǐng)
    中國(guó)機(jī)械工程 2025年4期
    關(guān)鍵詞:輕量化設(shè)計(jì)

    關(guān)鍵詞:混凝土泵車(chē);空心銷(xiāo)軸;輕量化設(shè)計(jì);多保真度代理模型

    中圖分類(lèi)號(hào):TH122

    DOI:10.3969/j.issn.1004-132X.2025.04.019 開(kāi)放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Lightweight Design of Concrete Pump Truck Boom Pins Based on Multi-fidelity Surrogate Model

    LI Peng1,2,3,4LI Mengcong1,4XIAO Libo2,3WANG Yitang1,4*SONG Xueguan1,4YANG Ling5 1.School of Mechanical Engineering,Dalian University of Technology,Dalian,Liaoning,116024 2.National Key Laboratory of Hoisting Machinery Key Technology,Changsha,410013 3.Zoomlion Heavy Industry Science amp; Technology Co.,Ltd.,Changsha,410013

    4.State Key Laboratory of High-performance Precision Manufacturing,Dalian,Liaoning,116024 5.Hunan College of Information,Changsha,410200

    Abstract: As a key connecting and supporting component in a concrete pump truck boom system, how to reduce the mass of the pin shafts was a noteworthy issue in boom design. A parameterized finite element model of the holow dumbbell pin shafts was constructed based on control parameters. Then,a more accurate feature mapping based MFS(FM-MFS) model was established through the mutual disturbance of high and low fidelity models and the reasonable allocation of high and low fidelity data.Genetic algorithm was used for optimization to obtain the optimal design scheme for the hollow dumbbell pins based on this surrogate model,achieving a weight reduction of 36% . Through theoretical,simulation,and experimental verification,it is shown that the hollow dumbbell pin shafts constructed based on the optimal parameters of the surrogate model may achieve lightweight and material consumption savings while still ensuring the physical and mechanical properties, which provides data support and reference for further applications in the future.

    Key words: concrete pump truck; hollow pin; lightweight design; multi-fidelity surrogate(MFS) model

    0 引言

    混凝土泵車(chē)作為基礎(chǔ)設(shè)施建設(shè)的重要工程機(jī)械,用于將混凝土輸送至施工現(xiàn)場(chǎng),在施工時(shí)能夠提高工程施工效率,降低人力和時(shí)間成本。隨著市場(chǎng)客戶的現(xiàn)實(shí)需求以及國(guó)家法規(guī)的日益嚴(yán)苛,超長(zhǎng)臂架以及重量合規(guī)已成為混凝土泵車(chē)發(fā)展的重要方向。超長(zhǎng)米段泵車(chē)在巨大的載荷作用下,要想實(shí)現(xiàn)重量合規(guī)是非常困難的,因此輕量化設(shè)計(jì)需求急迫且難度較大。

    臂架系統(tǒng)是混凝土泵車(chē)的關(guān)鍵組成部分,位于轉(zhuǎn)臺(tái)上方,擔(dān)負(fù)著輸送混凝土至指定位置的任務(wù)。臂架系統(tǒng)由臂架、連桿、銷(xiāo)軸、油缸、混凝土輸送管等零部件組成[1-2]。其中,銷(xiāo)軸用于連接臂架的各節(jié)臂,從而形成連續(xù)的結(jié)構(gòu),起到連接和傳遞力、彎矩、扭矩和運(yùn)動(dòng)的關(guān)鍵作用。常規(guī)混凝土泵車(chē)通常采用實(shí)心銷(xiāo)軸。在極限米段泵車(chē)臂架系統(tǒng)中,實(shí)心銷(xiāo)軸的總質(zhì)量可達(dá)1t以上。因此,在滿足銷(xiāo)軸強(qiáng)度、剛度、可靠性和使用壽命等要求下,實(shí)現(xiàn)其輕量化設(shè)計(jì)是十分有意義的。

    采用空心銷(xiāo)軸是實(shí)現(xiàn)結(jié)構(gòu)輕量化的一種方法??招匿N(xiāo)軸的設(shè)計(jì)允許將材料在軸心位置向邊緣移動(dòng),在相同的負(fù)荷條件下,通過(guò)合理的壁厚設(shè)計(jì),可以顯著節(jié)約材料,進(jìn)而減小零部件質(zhì)量,從而達(dá)到減小質(zhì)量的目標(biāo)。目前,雖已有學(xué)者對(duì)空心銷(xiāo)軸展開(kāi)設(shè)計(jì)[3-5],但研究較少且大多采用傳統(tǒng)結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法。雖然傳統(tǒng)方法可以滿足產(chǎn)品的功能和性能要求,但設(shè)計(jì)結(jié)構(gòu)未必是最合理和最輕量化的選項(xiàng),并且傳統(tǒng)設(shè)計(jì)方法存在著設(shè)計(jì)周期長(zhǎng)、改進(jìn)不明確等問(wèn)題[6-7]。代理模型技術(shù)在工程設(shè)計(jì)與優(yōu)化中扮演著重要的角色,它能顯著縮短優(yōu)化時(shí)間,并使復(fù)雜過(guò)程的優(yōu)化成為可能[8-9]。代理模型大致可分為單保真度代理模型和多保真度代理模型兩種,前者由單一高保真度樣本信息構(gòu)成,對(duì)維度高、非線性程度高的復(fù)雜問(wèn)題計(jì)算較為耗時(shí);后者可融合多種保真度的信息,能在保證精度的同時(shí)節(jié)省計(jì)算資源與時(shí)間成本。

    近年來(lái),多保真度代理模型引起了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注,并被應(yīng)用到結(jié)構(gòu)優(yōu)化設(shè)計(jì)中。LIU等[10]利用高斯過(guò)程和徑向基函數(shù)提出了一種基于多保真度代理模型的優(yōu)化框架,并將其用到了天線設(shè)計(jì)中。TAO等[11]提出了一種適用于馬赫數(shù)不確定條件下翼型和機(jī)翼的魯棒優(yōu)化的多保真度代理優(yōu)化框架。CHAKRABORTY等[12]針對(duì)魯棒設(shè)計(jì)優(yōu)化問(wèn)題,結(jié)合進(jìn)化計(jì)算與多保真度代理模型,提出了一種優(yōu)化框架并驗(yàn)證了方法的性能。

    本文利用多保真度代理模型技術(shù)對(duì)涉及多控制參數(shù)以及非線性接觸的泵車(chē)臂架銷(xiāo)軸進(jìn)行輕量化設(shè)計(jì)。經(jīng)理論計(jì)算、仿真和試驗(yàn)三重驗(yàn)證,在保有工件物理力學(xué)性能的同時(shí),成功地將銷(xiāo)軸質(zhì)量降低了 36% ,證明了多保真度代理模型相對(duì)于傳統(tǒng)代理模型的優(yōu)越性能,及其與中小型機(jī)械部件結(jié)構(gòu)優(yōu)化的適配性。

    1銷(xiāo)軸結(jié)構(gòu)參數(shù)化有限元模型

    空心啞鈴型銷(xiāo)軸有諸多變量,如圖1所示,僅憑重復(fù)嘗試難以確定全局最優(yōu)值。為獲得足夠擬合代理模型的數(shù)據(jù)集,需要建立幾何基體的參數(shù)化有限元模型。

    圖1空心啞鈴型銷(xiāo)軸幾何參數(shù)圖Fig.1 Geometric parameter diagramof hollow dumbbell shaped pinshaft

    空心啞鈴型銷(xiāo)軸各尺寸間的關(guān)系表述如下:

    式中: LT 為銷(xiāo)軸軸向總長(zhǎng); L4 為銷(xiāo)軸的直徑,因裝配需要,在此不予變動(dòng); L1 為銷(xiāo)軸與軸套接觸部位長(zhǎng)度; L2 為啞鈴段的長(zhǎng)度; L3 為啞鈴段的直徑; D 為內(nèi)部空心的直徑; R 為倒角尺寸; ρ 為比例因子。

    不難看出,如果將 L2,L3 和 ρ 作為控制參數(shù),則足夠引導(dǎo)參數(shù)模型的整體變化,因此選取這三個(gè)參數(shù)作為后續(xù)優(yōu)化設(shè)計(jì)的設(shè)計(jì)變量。

    銷(xiāo)軸和軸套的三維模型如圖2所示,據(jù)此建立銷(xiāo)軸及其配合軸套的參數(shù)化有限元模型。銷(xiāo)軸材料為 42CrMo ,屈服強(qiáng)度為 650MPa ,抗壓強(qiáng)度為 900MPa ,密度為 7850kg/m3 ,彈性模量為200MPa ,泊松比為0.3。對(duì)銷(xiāo)軸進(jìn)行自由網(wǎng)格劃分,對(duì)軸套進(jìn)行掃掠網(wǎng)格劃分,共生成兩種網(wǎng)格尺度的參數(shù)化有限元模型。這里,網(wǎng)格尺寸為 6mm 的模型稱(chēng)為高保真模型,如圖3所示;網(wǎng)格尺寸為10mm 的模型稱(chēng)為低保真模型,如圖4所示。

    圖2空心銷(xiāo)軸三維模型Fig.2 3Dmodel ofhollowpinshaft

    在銷(xiāo)軸與軸套間建立接觸對(duì),設(shè)定摩擦因數(shù)為0.2。在軸套與轉(zhuǎn)臺(tái)連接部位進(jìn)行固定約束,在軸套與油缸柄頭連接部位施加大小為 1425kN 的力載荷(兩側(cè)),對(duì)銷(xiāo)軸與轉(zhuǎn)臺(tái)配合端施加遠(yuǎn)端位移約束(只約束銷(xiāo)軸軸向位移,其他自由)??紤]銷(xiāo)軸的安裝方式,最終的邊界條件與載荷施加如圖5所示。

    圖3高保真有限元模型Fig.3High fidelity finite element model
    圖4低保真有限元模型Fig.4Low fidelityfinite elementmodel
    圖5 載荷與邊界條件Fig.5 Loadandboundaryconditions

    對(duì)于控制參數(shù) L2,L3 和 ρ ,在給定的區(qū)間內(nèi)任取一合理值,在同樣尺寸參數(shù)下對(duì)高、低保真模型進(jìn)行非線性接觸分析,對(duì)比觀察仿真結(jié)果的合理性。高、低保真模型的位移及應(yīng)力分布情況分別如圖6和圖7所示。

    經(jīng)測(cè)算,對(duì)于隨機(jī)給定設(shè)計(jì)參數(shù)的模型,運(yùn)行一次高保真仿真計(jì)算平均需要2640s左右,使用相同計(jì)算機(jī)運(yùn)行一次低保真仿真計(jì)算平均需要220s左右。另外,從圖6和圖7中可以看出,高保真模型的位移結(jié)果約為 0.98mm ,而低保真模型的位移結(jié)果為 1.18mm 。高保真模型的中間啞鈴段Mises應(yīng)力為 373.75MPa ,而低保真模型的中間啞鈴段Mises應(yīng)力為 442.87MPa 。這說(shuō)明,網(wǎng)格尺寸會(huì)影響非線性接觸計(jì)算,從而影響結(jié)構(gòu)傳力特性。選取網(wǎng)格尺寸 6~10mm 進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證,得到網(wǎng)格尺寸對(duì)仿真結(jié)果的影響如圖8所示,隨著網(wǎng)格的細(xì)化,仿真結(jié)果逐漸收斂和可信。綜上,有必要研究多保真代理模型,相較于完全采用高保真數(shù)據(jù)的單一保真度代理模型,多保真度代理模型可以在滿足一定計(jì)算精度的前提下,使得多參數(shù)和非線性計(jì)算成本更為低廉。

    圖6高保真模型的計(jì)算結(jié)果Fig.6 Calculationresultsofhighfidelitymodel
    圖7低保真模型的計(jì)算結(jié)果Fig.7 Calculationresults of low fidelity model
    圖8 網(wǎng)格無(wú)關(guān)性驗(yàn)證Fig.8 Verificationoffiniteelement meshindependence

    2 基于特征映射的多保真度代理模型

    代理模型的方法多樣,從建模數(shù)據(jù)可靠性角度可分為單一保真度代理模型和多保真度代理模型,其中單一保真度代理模型包括響應(yīng)面模型、克里金(KRG)模型、徑向基函數(shù)(RBF)模型、支持向量機(jī)回歸(SVR)模型等,多保真度代理模型包括MFS-KRG、MFS-RBF、基于特征映射的多保真度代理模型(feature mapping based multi-fidelitysurrogate,F(xiàn)M-MFS)等。多保真度代理模型由于可以用大量成本較低的低保真度數(shù)據(jù)混合高保真度數(shù)據(jù)進(jìn)行建模,因此在工程實(shí)際應(yīng)用中往往更具優(yōu)勢(shì)。主要體現(xiàn)在:在相同的數(shù)據(jù)成本下,多保真度代理模型精度更高;或者在給定的目標(biāo)精度下,多保真度代理模型成本更低。本文采用FM-MFS對(duì)空心啞鈴銷(xiāo)軸進(jìn)行優(yōu)化設(shè)計(jì),建立過(guò)程如圖9所示。

    圖9多保真度代理模型的建立過(guò)程Fig.9 The process of establishing multi-fidelity surrogatemodel

    首先通過(guò)低保真訓(xùn)練集構(gòu)建低保真模型 ,采用徑向基函數(shù)模型,具體計(jì)算公式如下:

    式中: x 為測(cè)試樣本; xli 為第 i 個(gè)低保真訓(xùn)練樣本; m 為低保真樣本數(shù): ?λi 為第 i 個(gè)基函數(shù)的權(quán)重系數(shù); 為歐氏距離; φ(?) 表示基函數(shù),選用二次函數(shù)。

    式(2)中權(quán)重系數(shù)通過(guò)插值條件求解,第 i 個(gè)訓(xùn)練樣本的低保真響應(yīng)為

    式中: 為 xli 對(duì)應(yīng)的低保真響應(yīng)。

    將式(2)代入式(4),權(quán)重系數(shù)可通過(guò)下式求解:

    式中: 為相關(guān)矩陣, 為低保真響應(yīng)構(gòu)成的向量。

    根據(jù)低保真模型得到的預(yù)測(cè)值,將高保真訓(xùn)練點(diǎn) 通過(guò)多項(xiàng)式特征進(jìn)行映射,得到過(guò)渡預(yù)測(cè):

    i=1,2,…,n

    式中: 為第 j 個(gè)單項(xiàng)式基; aj 為對(duì)應(yīng)的系數(shù);

    為高保真點(diǎn) xhi 通過(guò)低保真模型得到的預(yù)測(cè)值。

    式(6)可以改寫(xiě)為矩陣形式:

    其中, G 為 構(gòu)成的矩陣[13]。由于高保真樣本數(shù)量有限,因此式(6)中的項(xiàng)數(shù)有限,并且允許使用常數(shù)或低階多項(xiàng)式來(lái)近似過(guò)渡預(yù)測(cè)[14]

    下一步需要找到一個(gè)最優(yōu)向量 ,使過(guò)渡預(yù)測(cè)與高保真響應(yīng)之間的相關(guān)性最大化。使用皮爾遜相關(guān)系數(shù)計(jì)算兩者間的相關(guān)性,計(jì)算公式如下:

    式中: 為第 i 個(gè)高保真響應(yīng); 為高保真響應(yīng)的均值; 為過(guò)渡預(yù)測(cè)的均值, 為包含矩陣 每一列均值的行向量。

    為變量,以式(8)為目標(biāo)函數(shù),使過(guò)渡預(yù)測(cè)與高保真響應(yīng)的相關(guān)性最大化,即

    為便于求解,將 )的均值進(jìn)行規(guī)范化,即 μH=0 和 。規(guī)范化后的均值不會(huì)影響過(guò)渡預(yù)測(cè)與高保真響應(yīng)間的相關(guān)性。此時(shí),式(9)可改寫(xiě)為

    s.t.

    將式(10)改寫(xiě)為矩陣形式:

    過(guò)渡預(yù)測(cè) 采用由比例因子 ρ 和差異函數(shù)R(xH,xH)ω 組成的FM-MFS框架[15]進(jìn)行校正,其中 為基于高保真樣本構(gòu)建的相關(guān)矩陣, 為包含未知系數(shù)的向量。高保真響應(yīng)的計(jì)算公式為

    式中: 為優(yōu)化得到最優(yōu)解。

    為了得到未知參數(shù) ρ 和 ω ,將過(guò)渡預(yù)測(cè) )和關(guān)系矩陣 擴(kuò)展為增廣

    矩陣 ) R(xH,xH)] ,式(15)可變?yōu)?/p>

    FH=κβ

    其中, 為比例因子和系數(shù)向量組成的增廣向量, T]T。增廣矩陣 κ 為行滿秩,因?yàn)樵摼仃嚨闹葹?n (高保真樣本數(shù)),因此根據(jù)矩陣?yán)碚摵臀墨I(xiàn)[15-16],存在唯一的最小范數(shù)解 .

    然后可以得到比例因子 ρ 和系數(shù)向量 ,其中 的第一個(gè)元素是 ρ ,其余元素是 的分量。

    最終的多保真預(yù)測(cè)模型可表示為[17]

    式中: 為點(diǎn) xt 處的預(yù)測(cè); R(xt,xH) 為 xt 構(gòu)成的相關(guān)矩陣。

    3 銷(xiāo)軸參數(shù)化優(yōu)化及校核

    3.1 多保真度代理模型優(yōu)化

    以質(zhì)量最小為優(yōu)化目標(biāo),以結(jié)構(gòu)參數(shù)為設(shè)計(jì)變量,以結(jié)構(gòu)位移和應(yīng)力等為約束條件,對(duì)空心啞鈴銷(xiāo)軸結(jié)構(gòu)進(jìn)行輕量化優(yōu)化設(shè)計(jì)。

    空心啞鈴銷(xiāo)軸共有6個(gè)結(jié)構(gòu)參數(shù) (D,R,L1 ,L2,L3,L4) ,銷(xiāo)軸外徑 L4 因裝配原因不予變動(dòng),因此設(shè)計(jì)變量為 (D,R,L1,L2,L3) ,記為

    其中,三個(gè)控制變量的設(shè)計(jì)范圍為

    可推得各設(shè)計(jì)變量范圍為

    高、低保真度樣本的成本分配設(shè)計(jì)如下:模型的參數(shù)變動(dòng)會(huì)直接影響網(wǎng)格數(shù)量和求解時(shí)長(zhǎng),所以取隨機(jī)給定參數(shù)模型的計(jì)算時(shí)長(zhǎng)平均值,測(cè)得高、低保真樣本的成本比值(計(jì)算時(shí)長(zhǎng) )δ=12 。原單保真模型的高保真樣本點(diǎn)個(gè)數(shù)為 5ndv,ndv 為設(shè)計(jì)變量個(gè)數(shù),共5個(gè),它與計(jì)算成本的乘積即為初始計(jì)算成本。多保真模型的高、低保真樣本共享該成本,高保真樣本成本設(shè)占總量 80% ,因此兩類(lèi)樣本成本比值 θ=4 。容易推知下述關(guān)系:

    式中: nl,nh 為低保真度樣本和高保真樣本的數(shù)量。

    根據(jù)高低保真樣本的計(jì)算時(shí)長(zhǎng)比例以及式(22),可以計(jì)算并確保多保真代理模型和單保真代理模型的計(jì)算成本一致。

    以提高安全系數(shù)并減小質(zhì)量為目標(biāo),利用結(jié)構(gòu)參數(shù)與應(yīng)力等之間的代理模型,對(duì)銷(xiāo)軸進(jìn)行目標(biāo)優(yōu)化,尋找權(quán)衡的最佳決策方案。優(yōu)化問(wèn)題如下:

    式中: 為隨設(shè)計(jì)變量變化的銷(xiāo)軸質(zhì)量; Umax0 為銷(xiāo)軸最大允許變形量; σmax0 為銷(xiāo)軸許用彎曲應(yīng)力; τmax0 為銷(xiāo)軸許用剪切應(yīng)力; :Umax∝σmaxmax 分別為銷(xiāo)軸優(yōu)化期間的最大變形、彎曲應(yīng)力和剪切應(yīng)力; LbΩb 分別為設(shè)計(jì)變量上下區(qū)間。

    首先使用拉丁超立方采樣對(duì)高低保真樣本進(jìn)行試驗(yàn)設(shè)計(jì),獲取用于訓(xùn)練和測(cè)試的數(shù)據(jù);然后構(gòu)建設(shè)計(jì)參數(shù)和質(zhì)量、最大位移、最大應(yīng)力、最大剪應(yīng)力之間的FM-MFS代理模型。為了證實(shí)FM-MFS代理模型在面對(duì)多參數(shù)非線性復(fù)雜工程問(wèn)題時(shí)的卓越性能,使用單保真度RBF、KRG、SVR代理模型進(jìn)行對(duì)比,訓(xùn)練集和測(cè)試集與多保真度模型相同。采用均方根誤差 (RMSE)eRMS 和相對(duì)平均絕對(duì)誤差 (RMAE)eRMA 對(duì)比和評(píng)價(jià)各代理模型的精度,計(jì)算公式為

    式中: n 為測(cè)試點(diǎn)的數(shù)量; 分別為測(cè)試點(diǎn)的預(yù)測(cè)響應(yīng)和真實(shí)響應(yīng); 為真實(shí)響應(yīng)的平均值。

    誤差計(jì)算結(jié)果如圖10所示。RMSE反映代理模型的全局精度,RMAE則側(cè)重于評(píng)估模型局部精度,精度隨數(shù)值的下降而上升。由圖10可知,對(duì)于銷(xiāo)軸結(jié)構(gòu)的位移變形和von-Mises應(yīng)力預(yù)測(cè)模型,在耗費(fèi)相同成本(算時(shí))的情況下,多保真度代理模型精度顯著優(yōu)于RBF和SVR模型,對(duì)尤其適用于小樣本數(shù)據(jù)集的Kriging模型也有一定的精度優(yōu)勢(shì),對(duì)剪切應(yīng)力的預(yù)測(cè)精度則與單保真代理模型相當(dāng)。綜合來(lái)看,多保真代理模型彰顯了其面對(duì)高度非線性和多參數(shù)工程問(wèn)題時(shí)的卓越性能,其精度已經(jīng)完全可以支撐后續(xù)的優(yōu)化工作。

    采用高精度FM-MFS代理模型,結(jié)合NSGA-Ⅱ遺傳優(yōu)化算法,設(shè)置種群數(shù)量100,迭代次數(shù)500,交叉概率0.8,搜索優(yōu)化目標(biāo)的最優(yōu)方向,經(jīng)多次迭代后,最終結(jié)果在物理性能許可的范圍內(nèi),實(shí)現(xiàn)了泵車(chē)臂架的銷(xiāo)軸輕量化設(shè)計(jì)。為便于加工,優(yōu)化后的尺寸參數(shù)取整;暫不考慮制造和安裝誤差,因?yàn)榧庸ぶ圃煺`差范圍通常在 0.1mm 以內(nèi),安裝誤差引起的間隙通常在 0.6mm 以內(nèi),都不會(huì)明顯影響結(jié)構(gòu)應(yīng)力狀態(tài)。優(yōu)化后的銷(xiāo)軸結(jié)構(gòu)尺寸參數(shù)如表1所示。

    表1銷(xiāo)軸尺寸的最佳設(shè)計(jì)方案Tab.1 Optimaldesignschemeforpinshaftdimensionsmm

    3.2 優(yōu)化結(jié)果理論計(jì)算校核

    如前所述,銷(xiāo)軸材料為 42CrMo ,屈服強(qiáng)度為650MPa ,抗壓強(qiáng)度為 900MPa ,施加的載荷合力大小為 2.85MN 。受力方式為空心軸通軸雙剪切面,如圖11所示。

    銷(xiāo)軸的許用剪切應(yīng)力為

    [τ]=αaσbM

    式中: αa 為安全系數(shù),本文 αa=0.55;σb 為材料的極限強(qiáng)

    圖11 通軸兩剪切面銷(xiāo)軸結(jié)構(gòu)Fig.11 The double shearplane pin shaft structure

    度 ,σb=900MPa;γM 為材料安全系數(shù),疲勞強(qiáng)度下 γM= 1.1,靜態(tài)強(qiáng)度下 γM=1 。

    銷(xiāo)軸的許用彎曲應(yīng)力為

    [σ]=σs/1.25γM

    式中: σs 為材料的屈服強(qiáng)度,值為 650MPa 。

    圖11中AB截面位置的最大剪力為

    VABFF/2

    式中: F 為鉸點(diǎn)處所有銷(xiāo)軸受到某一部件作用力的合力,取值 2.85MN;γF 為載荷不均勻系數(shù),取值1.2。

    圖11中AB截面位置的最大彎矩為

    式中: t1 為外軸套長(zhǎng)度,取值 70mm;t2 為內(nèi)軸套長(zhǎng)度,取值 80mm;S12 為軸套間間隙,取值 1mm;c 為距離系數(shù),可通過(guò)彎矩試驗(yàn)進(jìn)行標(biāo)定。

    截面屬性系數(shù)為

    式中: Rw 為銷(xiāo)軸外徑的半徑,值為 L4/2;r 為空心內(nèi)徑的半徑,值為 D/2 。

    銷(xiāo)軸截面平均剪切應(yīng)力為

    式中: A 為截面面積, A=π(Rw2-r2) 。

    銷(xiāo)軸截面最大彎曲應(yīng)力為

    σ=MAB/W

    式中: W 截面抗彎系數(shù),按照?qǐng)A環(huán)計(jì)算。

    當(dāng)空心銷(xiāo)軸同時(shí)滿足 抗剪強(qiáng)度條件: τ/[τ]?1 抗彎強(qiáng)度條件 :σ/[σ]?1

    彎剪組合強(qiáng)度條件: (τ/[τ])2+(σ/[σ])2?1 時(shí),認(rèn)為該銷(xiāo)軸的設(shè)計(jì)是安全的。優(yōu)化后的空心銷(xiāo)軸抗剪強(qiáng)度、抗彎強(qiáng)度以及彎剪組合強(qiáng)度計(jì)算值如表2所示。根據(jù)理論計(jì)算空心銷(xiāo)軸的最大彎曲應(yīng)力為 371.27MPa ,小于許用彎曲應(yīng)力;最大剪切應(yīng)力為 257.83MPa ,小于許用剪切應(yīng)力。計(jì)算結(jié)果表明,該空心啞鈴銷(xiāo)軸符合理論強(qiáng)度校核要求。

    表2空心銷(xiāo)軸強(qiáng)度校核結(jié)果表Tab.2Strength verificationresults of hollow pin shaft

    3.3優(yōu)化結(jié)果仿真計(jì)算校核

    進(jìn)一步通過(guò)仿真計(jì)算進(jìn)行優(yōu)化結(jié)果校核。優(yōu)化后空心啞鈴銷(xiāo)軸的仿真計(jì)算結(jié)果如圖12所示。

    仿真結(jié)果與理論計(jì)算結(jié)果偏差較小。銷(xiāo)軸受壓時(shí)中間段最大Mises應(yīng)力為 361.21MPa ,小于許用彎曲應(yīng)力;兩端軸套處的最大剪切應(yīng)力為294.91MPa ,小于許用剪切應(yīng)力。仿真結(jié)果同樣表明,空心啞鈴銷(xiāo)軸在實(shí)現(xiàn)輕量化的同時(shí),依然能夠保證合格的物理力學(xué)性能。

    4空心啞鈴銷(xiāo)軸疲勞測(cè)試驗(yàn)證

    4.1 試驗(yàn)設(shè)備及加載過(guò)程

    試驗(yàn)加載設(shè)備為 103 t多功能結(jié)構(gòu)試驗(yàn)系統(tǒng),如圖13所示。加載系統(tǒng)最大垂向加載壓力10MN、最大垂向加載拉力 3MN ,標(biāo)準(zhǔn)試驗(yàn)空間1.4m×4.1m×3.0m 。該試驗(yàn)系統(tǒng)可采用荷載控制與位移控制兩種加載控制模式,可進(jìn)行靜態(tài)程控加載和動(dòng)態(tài)疲勞加載。測(cè)試樣件為原實(shí)心銷(xiāo)軸和優(yōu)化后的空心啞鈴銷(xiāo)軸,如圖14所示,原實(shí)心銷(xiāo)軸的質(zhì)量為 69.48kg ,優(yōu)化后的空心啞鈴銷(xiāo)軸質(zhì)量為 44.44kg ,減重 36% 。

    疲勞試驗(yàn)采用荷載控制模式進(jìn)行加載,每輪加載次數(shù)1000,加載策略為 50kN2850kN ,直到疲勞加載次數(shù)達(dá)11250。采用位移控制模式進(jìn)行卸載,卸載速度為 2mm/min ,卸載完成后拆下銷(xiāo)軸觀察表面磨損情況和疲勞開(kāi)裂情況。

    圖13加載設(shè)備Fig.13Fatigue test loading equipment
    圖14原實(shí)心銷(xiāo)軸和空心啞鈴銷(xiāo)軸Fig.14 Original solidpinshaft andhollow dumbbellpinshaft

    4.2 磨損情況和疲勞壽命

    疲勞加載11250次循環(huán)后,優(yōu)化后的空心啞鈴銷(xiāo)軸并未發(fā)現(xiàn)疲勞開(kāi)裂破壞現(xiàn)象,也未發(fā)生明顯的塑性變形和宏觀變形,但銷(xiāo)軸兩端存在輕微磨損現(xiàn)象,見(jiàn)圖15。試驗(yàn)結(jié)果表明,空心啞鈴銷(xiāo)軸通過(guò)了11250次疲勞試驗(yàn)加載,銷(xiāo)軸疲勞性能合格。

    圖15優(yōu)化銷(xiāo)軸加載11250次后的磨損情況Fig.15 Wearofhollowpinshaftafter11250 fatigue cycles

    將空心啞鈴銷(xiāo)軸與原實(shí)心銷(xiāo)軸的試驗(yàn)結(jié)果進(jìn)行對(duì)比發(fā)現(xiàn):優(yōu)化銷(xiāo)軸的表面損傷程度與原銷(xiāo)軸的表面損傷程度基本相同,但是兩類(lèi)試件的損傷面有所不同,優(yōu)化銷(xiāo)軸損傷面趨于三角形,而原實(shí)心銷(xiāo)軸損傷面趨于矩形,如圖16所示,產(chǎn)生該現(xiàn)象的原因在于加載過(guò)程中優(yōu)化銷(xiāo)軸的彎曲變形較原銷(xiāo)軸的彎曲變形略大。優(yōu)化銷(xiāo)軸和原銷(xiāo)軸在加載11250次后均未出現(xiàn)疲勞開(kāi)裂破壞,接觸面磨損程度基本相同,因此空心啞鈴優(yōu)化銷(xiāo)軸同樣具備足夠的工作性能和物理力學(xué)性能。

    圖16優(yōu)化銷(xiāo)軸與原實(shí)心銷(xiāo)軸的損傷對(duì)比Fig.16Comparison of damage between optimized pin shaftand original solid pinshaft

    5結(jié)論

    1)根據(jù)銷(xiāo)軸幾何尺寸關(guān)系,構(gòu)建了基于控制參數(shù)的參數(shù)化有限元模型??刂茀?shù)在給定的區(qū)間內(nèi)可取隨機(jī)值,進(jìn)行非線性接觸分析,位移和應(yīng)力趨勢(shì)分布合理,能夠支持后續(xù)的代理模型建立與優(yōu)化分析。

    2)基于多保真度代理模型理論,在相同的算時(shí)內(nèi),通過(guò)高、低保真度模型的相互擾動(dòng)和高、低保真數(shù)據(jù)的合理算時(shí)分配,獲得了相較于傳統(tǒng)RBF、Kriging等代理模型精度更高的FM-MFS代理模型。以該代理模型為基礎(chǔ),采用遺傳算法進(jìn)行優(yōu)化,最終獲得了空心啞鈴銷(xiāo)軸的最佳設(shè)計(jì)方案,實(shí)現(xiàn)了銷(xiāo)軸的輕量化,并通過(guò)了理論和仿真強(qiáng)度校核。

    3)搭建了銷(xiāo)軸疲勞性能測(cè)試臺(tái),經(jīng)歷11250次疲勞試驗(yàn)后,空心啞鈴優(yōu)化銷(xiāo)軸與原實(shí)心銷(xiāo)軸相比,表面磨損程度基本相同,且均未出現(xiàn)疲勞開(kāi)裂破壞現(xiàn)象,說(shuō)明基于代理模型最優(yōu)參數(shù)所構(gòu)建的空心啞鈴銷(xiāo)軸在實(shí)現(xiàn)輕量化和節(jié)省材料的同時(shí),其物理力學(xué)性能依然能夠得到保證,為后續(xù)進(jìn)一步應(yīng)用提供了數(shù)據(jù)支撐和參考。

    參考文獻(xiàn):

    [1]周振平,孫武和.混凝土機(jī)械構(gòu)造與使用維護(hù)[M]. 北京:化學(xué)工業(yè)出版社,2014. ZHOU Zhenping,SUN Wuhe. Construction and Maintenance of Concrete Machinery[M]. Beijing: Chemical Industry Press,2014.

    [2] 李鵬,王一棠,王佳茜,等.數(shù)據(jù)驅(qū)動(dòng)的混凝土泵車(chē)臂 架疲歷軾何諧編制LJ」.中國(guó)機(jī)械工桂,ZU24,35 (10):1881-1889. LI Peng,WANG Yitang,WANG Jiaxi,et al. Compilation of Data-driven Fatigue Load Spectrum for Concrete Pump Truck Boom[J]. China Mechanical Engineering,2024,35(10):1881-1889.

    [3]黃大為,田相玉,孫國(guó)香.空心銷(xiāo)軸的制備、計(jì)算及試 驗(yàn)[J].建設(shè)機(jī)械技術(shù)與管理,2016,29(5):60-65. HUANG Dawei, TIAN Xiangyu, SUN Guoxiang. Preparation,Calculation,and Testing of Hollow Pins[J]. Construction Machinery Technology and Management,2016,29(5):60-65.

    [4]朱奇,張?chǎng)?,黃宇華.長(zhǎng)臂架泵車(chē)空心銷(xiāo)軸裝配工裝 設(shè)計(jì)與應(yīng)用[J].建筑機(jī)械化,2021,42(9):66-67. ZHU Qi, ZHANG Xin,HUANG Yuhua. Design and Application of Assembly Fixture for Hollow Pin Shafts of Long Arm Pump Trucks[J]. Construction Mechanization,2021,42(9):66-67.

    [5]朱奇,何科英,蔣廉健.高端工程機(jī)械大長(zhǎng)徑比重型 空心銷(xiāo)軸拆拔工裝設(shè)計(jì)[J].建筑機(jī)械,2021(9):91- 93. ZHU Qi,HE Keying, JIANG Lianjian. Design of Dismantling and Drawing Fixtures for Heavy-duty Hollow Pin Shafts with Large Aspect Ratios in High-end Construction Machinery[J]. Construction Machinery,2021(9):91-93.

    [6]李政茂.商用車(chē)駕駛室骨架輕量化設(shè)計(jì)與研究[D]. 長(zhǎng)春:吉林大學(xué) 2010. LI Zhengmao. Lightweight Design and Research of Commercial Vehicle Driver's Cabin Skeleton[D]. Changchun:Jilin University,2010.

    [7]李龍.車(chē)身梁截面優(yōu)化方法的研究[D].長(zhǎng)春:吉林大 學(xué),2008. LI Long. Research on Optimization Method of Vehicle Beam Section[D]. Changchun:Jilin University, 2008.

    [8]何西旺,楊亮亮,冉仁杰,等.基于多評(píng)價(jià)標(biāo)準(zhǔn)的代理 模型綜合比較研究[J].機(jī)械工程學(xué)報(bào),2022,58 (16):403-419. HE Xiwang,YANG Liangliang,RAN Renjie,et al. A Comprehensive Comparative Study of Proxy Models Based on Multiple Evaluation Criteria[J]. Journal of Mechanical Engineering,2022,58(16): 403-419.

    [9]張偉,宋學(xué)官,石茂林,等.基于代理模型的機(jī)械式挖 掘機(jī)動(dòng)臂輕量化設(shè)計(jì)[J].機(jī)械設(shè)計(jì)與制造,2019 (4) :1-4. ZHANG Wei, SONG Xueguan, SHI Maolin,et al. Lightweight Design of Mechanical Excavator Arm Based on Agent ModelLJ」. Mechanical Design and Manufacturing,2019(4):1-4.

    [10]LIU B,KOZIEL S,ZHANG Q. A Multi-fidelity Surrogate-model-assisted Evolutionary Algorithm for Computationally Expensive Optimization Problems[J]. Journal of Computational Science,2016, 12:28-37.

    [11] TAO J,SUN G. Application of Deep Learning Based Multi-fidelity Surrogate Model to Robust Aerodynamic Design Optimization[J]. Aerospace Science and Technology,2019,92:722-737.

    [12] CHAKRABORTYS,CHATTERJEET, CHOWDHURY R, et al. A Surrogate Based Multi-fidelity Approach for Robust Design Optimization[J].Applied Mathematical Modelling,2017, 47:726-744.

    [13] PAPILA M,HAFTKA R T.Response Surface Approximations:Noise,Error Repair,and Modeling Errors[J].AIAA Journal,2000,38(12): 2336-2343.

    [14]ZHANG Y,KIM N H,PARK C,et al. Multifidelity Surrogate Based on Single Linear Regression [J].AIAA Journal,2018,56(12) :4944-4952.

    [15]SONG X,LYU L,SUN W,et al. A Radial Basis Function-based Multi-fidelity Surrogate Model:Exploring Correlation between High-fidelityand Low-fidelity Models[J]. Structural and Multidisciplinary Optimization,2019,60(3):965-981.

    [16]PETERSEN K B,PEDERSEN M S. The Matrix Cookbook[J]. Technical University of Denmark, 2008,7(15) :510.

    [17]WANG Y,LI K,LI Q,et al. Multi-fidelity Information Fusion with Hierarchical Surrogate Guided by Feature Mapping[J]. Knowledge-based Systems,2023,275:110693.

    (編輯袁興玲)

    作者簡(jiǎn)介:李鵬,男,1989年生,博士研究生。研究方向?yàn)榻Y(jié)構(gòu)振動(dòng)、結(jié)構(gòu)疲勞、結(jié)構(gòu)輕量化。E-mail;55131668@qq.com。王一橐*(通信作者),女,1993年生,博士研究生。研究方向?yàn)楣I(yè)大數(shù)據(jù)挖掘及數(shù)據(jù)驅(qū)動(dòng)的預(yù)測(cè)技術(shù)。E-mail:yitangwang @ mail.dlut.edu.cn。

    本文引用格式:

    李鵬,李夢(mèng)聰,肖立波,等.基于多保真度代理模型的混凝土泵車(chē)臂架銷(xiāo)軸輕量化設(shè)計(jì)[J].中國(guó)機(jī)械工程,2025,36(4):821-829.LIPeng,LIMengcong,XIAOLibo,etal.LightweightDesign ofConcrete Pump Truck BoomPins Based on Multi-fidelity Surro-gateModel[J].China Mechanical Engineering,2025,36(4):821-829.

    猜你喜歡
    輕量化設(shè)計(jì)
    FSC方程式賽車(chē)單體殼車(chē)架設(shè)計(jì)
    巴哈賽車(chē)車(chē)架仿真分析與優(yōu)化設(shè)計(jì)
    關(guān)于岸邊集裝箱起重機(jī)結(jié)構(gòu)輕量化設(shè)計(jì)及可靠性分析
    大型客機(jī)隔熱隔音棉的輕量化設(shè)計(jì)
    科技視界(2017年9期)2017-09-04 19:14:07
    伊朗設(shè)拉子地鐵車(chē)輛輕量化碳鋼車(chē)體設(shè)計(jì)
    科技視界(2017年2期)2017-04-18 09:07:27
    軌道扣件彈性墊板結(jié)構(gòu)優(yōu)化設(shè)計(jì)
    乘用車(chē)車(chē)身結(jié)構(gòu)輕量化設(shè)計(jì)技術(shù)研究與實(shí)踐
    FSAE賽車(chē)車(chē)架結(jié)構(gòu)拓?fù)鋬?yōu)化及輕量化設(shè)計(jì)研究
    帽形截面薄壁焊接結(jié)構(gòu)軸向耐撞性的優(yōu)化設(shè)計(jì)方法研究
    飛機(jī)平尾普通肋的優(yōu)化設(shè)計(jì)
    大型黄色视频在线免费观看| 性色avwww在线观看| 久久九九热精品免费| 国产单亲对白刺激| 亚洲在线观看片| 国产成人一区二区在线| 亚洲乱码一区二区免费版| 精品午夜福利在线看| 国产成人a∨麻豆精品| 国产日本99.免费观看| 大又大粗又爽又黄少妇毛片口| 变态另类成人亚洲欧美熟女| 草草在线视频免费看| 亚洲av第一区精品v没综合| 亚洲人与动物交配视频| 国产伦精品一区二区三区四那| 插逼视频在线观看| 国产成年人精品一区二区| 欧美xxxx黑人xx丫x性爽| 美女高潮的动态| 久久久久久久久中文| 亚洲成人久久性| 国产精品国产三级国产av玫瑰| 男女视频在线观看网站免费| 尤物成人国产欧美一区二区三区| 免费无遮挡裸体视频| 国产男靠女视频免费网站| 天堂动漫精品| 亚洲国产色片| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 日日撸夜夜添| 99riav亚洲国产免费| 高清毛片免费观看视频网站| 亚洲五月天丁香| 国产一区二区激情短视频| 日本精品一区二区三区蜜桃| 免费电影在线观看免费观看| 美女免费视频网站| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 欧美在线一区亚洲| 久久久色成人| 桃色一区二区三区在线观看| 亚洲国产精品国产精品| 可以在线观看毛片的网站| 精品少妇黑人巨大在线播放 | 亚洲七黄色美女视频| 禁无遮挡网站| 美女高潮的动态| 中文字幕免费在线视频6| 一个人看视频在线观看www免费| 国产精品福利在线免费观看| 99九九线精品视频在线观看视频| 国产精品一二三区在线看| 亚洲欧美成人精品一区二区| 人人妻人人澡欧美一区二区| 久久久久久久午夜电影| 欧美3d第一页| 两性午夜刺激爽爽歪歪视频在线观看| 久久这里只有精品中国| 欧美在线一区亚洲| 国产成人a∨麻豆精品| 日本免费a在线| 亚洲无线在线观看| 久久精品影院6| 国产黄色小视频在线观看| 男人的好看免费观看在线视频| 老司机福利观看| 蜜桃久久精品国产亚洲av| 国产精品野战在线观看| 亚洲av五月六月丁香网| 国内少妇人妻偷人精品xxx网站| 久久精品久久久久久噜噜老黄 | 夜夜看夜夜爽夜夜摸| 国产成人91sexporn| 一边摸一边抽搐一进一小说| 欧美最黄视频在线播放免费| 自拍偷自拍亚洲精品老妇| 噜噜噜噜噜久久久久久91| 免费在线观看成人毛片| 欧美绝顶高潮抽搐喷水| 久久久久久伊人网av| 亚洲五月天丁香| 中文在线观看免费www的网站| 三级男女做爰猛烈吃奶摸视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产蜜桃级精品一区二区三区| 亚洲无线在线观看| 在线免费观看不下载黄p国产| 久久人妻av系列| 国产精品1区2区在线观看.| 日韩欧美一区二区三区在线观看| 麻豆成人午夜福利视频| 97超碰精品成人国产| 无遮挡黄片免费观看| 一个人免费在线观看电影| 国产精品久久久久久久久免| 国产白丝娇喘喷水9色精品| 久久精品国产亚洲av香蕉五月| 18禁在线无遮挡免费观看视频 | 国产精品三级大全| 啦啦啦啦在线视频资源| 黄色配什么色好看| 精品午夜福利在线看| 可以在线观看的亚洲视频| 小蜜桃在线观看免费完整版高清| 成人三级黄色视频| 久久久久九九精品影院| 久久久久久久久中文| 亚洲美女黄片视频| 97热精品久久久久久| 国产高清三级在线| 亚洲精品一区av在线观看| 我要搜黄色片| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| 国产欧美日韩精品亚洲av| 国产成人aa在线观看| 国产精品人妻久久久久久| 91av网一区二区| 免费搜索国产男女视频| 精品国内亚洲2022精品成人| 精品人妻视频免费看| 一个人观看的视频www高清免费观看| 在线播放国产精品三级| 亚洲欧美精品自产自拍| 精品一区二区免费观看| 男女下面进入的视频免费午夜| 久久中文看片网| АⅤ资源中文在线天堂| 一级毛片久久久久久久久女| 国产久久久一区二区三区| 成人特级黄色片久久久久久久| 国产成人福利小说| 国产单亲对白刺激| 熟女电影av网| 国产一区二区亚洲精品在线观看| 欧美3d第一页| 国产老妇女一区| 免费高清视频大片| 久久精品国产99精品国产亚洲性色| 亚洲av五月六月丁香网| 久久久久久大精品| 久久久久久久久久久丰满| 国产一级毛片七仙女欲春2| 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 免费观看在线日韩| 日韩强制内射视频| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 狠狠狠狠99中文字幕| 老师上课跳d突然被开到最大视频| ponron亚洲| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 色哟哟哟哟哟哟| 成人特级av手机在线观看| 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 午夜免费激情av| 99精品在免费线老司机午夜| 久久国产乱子免费精品| 男人的好看免费观看在线视频| 看黄色毛片网站| 国产午夜精品久久久久久一区二区三区 | av女优亚洲男人天堂| 91在线观看av| 亚州av有码| 精品无人区乱码1区二区| 久久鲁丝午夜福利片| 国产精品人妻久久久影院| 少妇熟女欧美另类| 人人妻人人澡人人爽人人夜夜 | 亚洲人成网站高清观看| eeuss影院久久| 亚洲最大成人手机在线| 白带黄色成豆腐渣| 99九九线精品视频在线观看视频| 成人综合一区亚洲| 深夜a级毛片| 91在线精品国自产拍蜜月| 国产又黄又爽又无遮挡在线| 精华霜和精华液先用哪个| 性色avwww在线观看| 午夜日韩欧美国产| 成人国产麻豆网| 精品一区二区三区视频在线| 国产精品1区2区在线观看.| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆 | 99久国产av精品国产电影| 国产精品不卡视频一区二区| 国产欧美日韩精品亚洲av| 中出人妻视频一区二区| 在线观看av片永久免费下载| 亚洲va在线va天堂va国产| 国产色婷婷99| 久久精品国产亚洲av香蕉五月| 99精品在免费线老司机午夜| 国内揄拍国产精品人妻在线| 最新在线观看一区二区三区| 午夜福利在线观看免费完整高清在 | 91在线精品国自产拍蜜月| 精品久久久噜噜| 18禁黄网站禁片免费观看直播| 国产精品嫩草影院av在线观看| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 日本在线视频免费播放| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 少妇高潮的动态图| av在线天堂中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产av在哪里看| 日本一二三区视频观看| 日韩精品中文字幕看吧| 亚洲性夜色夜夜综合| 国产 一区 欧美 日韩| 亚洲精品日韩在线中文字幕 | 国产高清三级在线| 久久精品国产亚洲av香蕉五月| 一级黄片播放器| 欧美区成人在线视频| 免费大片18禁| 国产激情偷乱视频一区二区| 国产成年人精品一区二区| 如何舔出高潮| 国产精品福利在线免费观看| 午夜激情欧美在线| 亚洲成人久久性| 一个人观看的视频www高清免费观看| 亚洲人成网站在线观看播放| 欧美精品国产亚洲| 国产久久久一区二区三区| 亚洲精品在线观看二区| 亚洲国产色片| 亚洲欧美精品自产自拍| 99riav亚洲国产免费| 熟妇人妻久久中文字幕3abv| 国产一区二区三区在线臀色熟女| 日韩大尺度精品在线看网址| 97超视频在线观看视频| 国产91av在线免费观看| 国产男靠女视频免费网站| 给我免费播放毛片高清在线观看| 亚洲专区国产一区二区| av视频在线观看入口| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 亚洲美女搞黄在线观看 | 久久久久久久久久久丰满| 亚洲精品在线观看二区| 国产麻豆成人av免费视频| 搞女人的毛片| 舔av片在线| 日韩大尺度精品在线看网址| 97超视频在线观看视频| 精品久久久久久成人av| 久久国产乱子免费精品| 最新在线观看一区二区三区| 婷婷精品国产亚洲av在线| 美女cb高潮喷水在线观看| 国产精品乱码一区二三区的特点| 中国国产av一级| 亚洲欧美日韩高清在线视频| 国产精品一及| 欧美国产日韩亚洲一区| 男人舔女人下体高潮全视频| 欧美成人精品欧美一级黄| 日韩人妻高清精品专区| 一级毛片我不卡| 91狼人影院| 哪里可以看免费的av片| 亚洲国产精品久久男人天堂| 美女免费视频网站| 亚洲精品456在线播放app| 国产精品综合久久久久久久免费| 性插视频无遮挡在线免费观看| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 亚洲在线自拍视频| 国产精品,欧美在线| 欧美丝袜亚洲另类| 久久热精品热| 国产真实乱freesex| 在线观看66精品国产| 校园人妻丝袜中文字幕| 精品少妇黑人巨大在线播放 | 日本精品一区二区三区蜜桃| 真实男女啪啪啪动态图| 久久精品国产清高在天天线| 国产91av在线免费观看| avwww免费| 三级毛片av免费| 亚洲第一区二区三区不卡| 亚洲婷婷狠狠爱综合网| 三级经典国产精品| av卡一久久| 少妇熟女aⅴ在线视频| 欧美性感艳星| 在线播放无遮挡| 亚洲熟妇熟女久久| 久久久久精品国产欧美久久久| 日日啪夜夜撸| 久99久视频精品免费| 网址你懂的国产日韩在线| 听说在线观看完整版免费高清| 国内精品久久久久精免费| 国产精品久久久久久精品电影| 亚洲av电影不卡..在线观看| 精品久久久久久久久亚洲| 欧美丝袜亚洲另类| 国产午夜精品久久久久久一区二区三区 | 成人美女网站在线观看视频| 波多野结衣高清作品| 乱码一卡2卡4卡精品| 国产伦精品一区二区三区视频9| 欧美成人精品欧美一级黄| 国产午夜精品久久久久久一区二区三区 | 亚洲人成网站在线观看播放| 久久久久久大精品| 亚洲国产日韩欧美精品在线观看| 欧美性感艳星| 91狼人影院| 嫩草影院精品99| 午夜视频国产福利| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 欧美一区二区亚洲| 免费看光身美女| 日本免费a在线| 色5月婷婷丁香| a级毛片免费高清观看在线播放| 亚洲丝袜综合中文字幕| 有码 亚洲区| 99热这里只有是精品50| 国产高清视频在线观看网站| 99热这里只有是精品50| 99在线人妻在线中文字幕| 人妻久久中文字幕网| 国产精品无大码| 精品一区二区三区视频在线| 成人av一区二区三区在线看| 久久久久国内视频| 成人特级黄色片久久久久久久| 精品久久久久久久末码| 熟女电影av网| 禁无遮挡网站| 一区二区三区高清视频在线| 亚洲成人久久爱视频| 悠悠久久av| 国产精品嫩草影院av在线观看| 国产成人精品久久久久久| 春色校园在线视频观看| 国产成人精品久久久久久| 亚洲av免费高清在线观看| 夜夜看夜夜爽夜夜摸| 男插女下体视频免费在线播放| 亚洲经典国产精华液单| videossex国产| 国产精品1区2区在线观看.| 嫩草影院入口| av.在线天堂| 成人无遮挡网站| 国产免费男女视频| 亚洲国产欧美人成| 99热6这里只有精品| 日本免费a在线| 国产在线男女| 亚洲自偷自拍三级| 久久久色成人| 国产高清视频在线观看网站| 欧美色欧美亚洲另类二区| 全区人妻精品视频| 日本免费a在线| 亚洲人成网站在线观看播放| 一a级毛片在线观看| 亚洲欧美日韩高清专用| 欧美国产日韩亚洲一区| 亚洲精品在线观看二区| 久久国内精品自在自线图片| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 日韩一本色道免费dvd| 欧美xxxx性猛交bbbb| 国产三级中文精品| 麻豆av噜噜一区二区三区| 村上凉子中文字幕在线| 内地一区二区视频在线| 国产精品国产三级国产av玫瑰| 天堂网av新在线| 国产黄a三级三级三级人| 赤兔流量卡办理| 噜噜噜噜噜久久久久久91| 91狼人影院| 亚洲av第一区精品v没综合| 99久久九九国产精品国产免费| 午夜福利在线观看吧| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 国产乱人视频| 久久精品国产亚洲网站| 看非洲黑人一级黄片| 精品午夜福利视频在线观看一区| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 亚洲av中文字字幕乱码综合| 我要搜黄色片| .国产精品久久| 美女cb高潮喷水在线观看| 1000部很黄的大片| 99热这里只有是精品在线观看| 亚洲自拍偷在线| 欧美成人精品欧美一级黄| 国产高清激情床上av| 观看免费一级毛片| 有码 亚洲区| 日日摸夜夜添夜夜添av毛片| 欧美zozozo另类| 国产精品1区2区在线观看.| 春色校园在线视频观看| 91午夜精品亚洲一区二区三区| 最新在线观看一区二区三区| 男女做爰动态图高潮gif福利片| 国产aⅴ精品一区二区三区波| 免费高清视频大片| 91狼人影院| 少妇猛男粗大的猛烈进出视频 | 国产av在哪里看| 特大巨黑吊av在线直播| 国产私拍福利视频在线观看| 99精品在免费线老司机午夜| 亚洲性夜色夜夜综合| 色吧在线观看| 欧美最黄视频在线播放免费| 我要搜黄色片| 欧美日韩综合久久久久久| 久久久色成人| 午夜影院日韩av| 免费黄网站久久成人精品| 美女黄网站色视频| a级一级毛片免费在线观看| 免费大片18禁| 女同久久另类99精品国产91| 国产美女午夜福利| 国产精品无大码| 欧美国产日韩亚洲一区| 亚洲熟妇熟女久久| 欧美日韩国产亚洲二区| 久久午夜福利片| 三级国产精品欧美在线观看| 成人性生交大片免费视频hd| 精品人妻视频免费看| 九九爱精品视频在线观看| 国产成年人精品一区二区| 在线国产一区二区在线| 岛国在线免费视频观看| 高清午夜精品一区二区三区 | 中文字幕av在线有码专区| 男女之事视频高清在线观看| 国产精品亚洲一级av第二区| 欧美xxxx性猛交bbbb| 五月伊人婷婷丁香| 淫妇啪啪啪对白视频| 噜噜噜噜噜久久久久久91| 国产精品一区二区三区四区免费观看 | 国产精品不卡视频一区二区| 欧美成人a在线观看| 免费在线观看影片大全网站| 一级av片app| 精品一区二区免费观看| 国产欧美日韩精品一区二区| 日韩成人av中文字幕在线观看 | 狂野欧美白嫩少妇大欣赏| 亚洲成人精品中文字幕电影| 最近最新中文字幕大全电影3| 久久久久久久午夜电影| 亚洲av五月六月丁香网| 亚洲欧美中文字幕日韩二区| 97人妻精品一区二区三区麻豆| 精品人妻一区二区三区麻豆 | 日韩大尺度精品在线看网址| 综合色丁香网| 国产午夜福利久久久久久| 欧美色欧美亚洲另类二区| avwww免费| 啦啦啦韩国在线观看视频| 深爱激情五月婷婷| 女同久久另类99精品国产91| 秋霞在线观看毛片| 日韩欧美一区二区三区在线观看| 成人精品一区二区免费| 欧美成人a在线观看| 欧美一区二区精品小视频在线| 精品无人区乱码1区二区| 国产不卡一卡二| 日韩大尺度精品在线看网址| 日本爱情动作片www.在线观看 | 日韩 亚洲 欧美在线| 久久精品国产自在天天线| 国产成人a区在线观看| 国产精品精品国产色婷婷| 99久久九九国产精品国产免费| 好男人在线观看高清免费视频| 亚洲精品一区av在线观看| av黄色大香蕉| 噜噜噜噜噜久久久久久91| av福利片在线观看| 免费在线观看成人毛片| 国产私拍福利视频在线观看| 尾随美女入室| 美女cb高潮喷水在线观看| 国产精品福利在线免费观看| 男女那种视频在线观看| 欧美日韩乱码在线| 色哟哟哟哟哟哟| 欧美成人精品欧美一级黄| 国产欧美日韩精品亚洲av| 日韩,欧美,国产一区二区三区 | 亚洲aⅴ乱码一区二区在线播放| 国产免费一级a男人的天堂| 色视频www国产| 亚洲av.av天堂| 亚洲国产日韩欧美精品在线观看| 在线观看午夜福利视频| 日产精品乱码卡一卡2卡三| 亚洲天堂国产精品一区在线| 亚洲无线在线观看| 欧美日韩综合久久久久久| a级毛色黄片| 丰满人妻一区二区三区视频av| 中国美白少妇内射xxxbb| 成人无遮挡网站| 国产三级中文精品| 亚洲天堂国产精品一区在线| 超碰av人人做人人爽久久| 熟女电影av网| 男人狂女人下面高潮的视频| 久久精品影院6| 伦理电影大哥的女人| 人妻少妇偷人精品九色| 婷婷六月久久综合丁香| 久久亚洲国产成人精品v| 免费观看的影片在线观看| 美女大奶头视频| 成人一区二区视频在线观看| 欧美zozozo另类| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品| 中国国产av一级| 精品熟女少妇av免费看| 国产单亲对白刺激| 变态另类丝袜制服| 亚洲经典国产精华液单| 国产视频内射| 亚洲av免费在线观看| 内射极品少妇av片p| 亚洲国产精品成人久久小说 | 欧美成人a在线观看| 99国产极品粉嫩在线观看| av在线老鸭窝| 亚洲欧美日韩卡通动漫| 亚洲成人久久性| 别揉我奶头 嗯啊视频| 99九九线精品视频在线观看视频| 永久网站在线| 内地一区二区视频在线| 舔av片在线| 不卡一级毛片| 九九爱精品视频在线观看| 中国国产av一级| 婷婷精品国产亚洲av| 在现免费观看毛片| av中文乱码字幕在线| 国产精品国产三级国产av玫瑰| 亚洲国产精品久久男人天堂| 久久中文看片网| 国产激情偷乱视频一区二区| 欧美性感艳星| 黄色日韩在线| 亚洲国产精品成人久久小说 | 一个人看视频在线观看www免费| 69人妻影院| 波野结衣二区三区在线| 亚洲高清免费不卡视频| 我的老师免费观看完整版| 99热6这里只有精品| 五月伊人婷婷丁香| 精品一区二区三区av网在线观看| 亚洲精品国产成人久久av| 少妇高潮的动态图| 99热只有精品国产| 久久婷婷人人爽人人干人人爱| 精品午夜福利视频在线观看一区| 亚洲自偷自拍三级| 国产精品乱码一区二三区的特点| 三级毛片av免费| 日韩欧美三级三区| 国产精品乱码一区二三区的特点| 51国产日韩欧美| 99久久精品热视频| 日韩欧美在线乱码| 精品欧美国产一区二区三| 亚洲美女视频黄频| 精品少妇黑人巨大在线播放 | 淫秽高清视频在线观看| 国产成人精品久久久久久|