亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)

        2024-07-05 05:04:30王子寧楊飛吳澤勛王德遠(yuǎn)齊建偉
        時(shí)代汽車 2024年10期
        關(guān)鍵詞:輕量化設(shè)計(jì)響應(yīng)面法

        王子寧 楊飛 吳澤勛 王德遠(yuǎn) 齊建偉

        摘 要:針對(duì)龍門下框架設(shè)計(jì)過(guò)于保守、材料利用率不高的問(wèn)題,提出了基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)方法。通過(guò)將中心復(fù)合試驗(yàn)設(shè)計(jì)與有限元數(shù)值模擬技術(shù)相結(jié)合,構(gòu)建了龍門下框架質(zhì)量、剛度、強(qiáng)度的二階多項(xiàng)式響應(yīng)面模型,并對(duì)其進(jìn)行了精度檢驗(yàn),驗(yàn)證了所建立的響應(yīng)面模型具有足夠的精度。在此基礎(chǔ)上,以各構(gòu)件厚度為設(shè)計(jì)變量,質(zhì)量為目標(biāo),剛度及強(qiáng)度為約束建立了龍門下框架輕量化設(shè)計(jì)數(shù)學(xué)模型,采用序列二次規(guī)劃算法進(jìn)行求解,并對(duì)優(yōu)化結(jié)果進(jìn)行仿真驗(yàn)證。結(jié)果表明,優(yōu)化后的龍門下框架滿足剛度和強(qiáng)度要求,重量降低了82.1kg,減重率達(dá)到了19.0%,輕量化效果顯著。

        關(guān)鍵詞:響應(yīng)面法 龍門下框架 輕量化設(shè)計(jì)

        1 引言

        近年來(lái),隨著能源和環(huán)境壓力的加大和國(guó)家新能源戰(zhàn)略的逐步實(shí)施,電動(dòng)汽車領(lǐng)域得到了快速發(fā)展,作為車輛電能供給載體的換電站大量涌現(xiàn)在各大城市中[1]。龍門舉升機(jī)構(gòu)作為換電站的重要組成部分,主要用于換電時(shí)完成車輛的舉升動(dòng)作,其主要分為龍門上框架和龍門下框架。目前,龍門下框架主要是基于類比方法進(jìn)行經(jīng)驗(yàn)式設(shè)計(jì),雖然能滿足性能要求,但所設(shè)計(jì)出的龍門下框架的材料利用率普遍較低,從而引起材料產(chǎn)生不必要的消耗,增加了材料成本。在此背景下,在滿足結(jié)構(gòu)可靠性的基礎(chǔ)上對(duì)龍門下框架進(jìn)行輕量化設(shè)計(jì)就顯得尤為必要。

        龍門下框架輕量化設(shè)計(jì)的基本思想為在保證結(jié)構(gòu)可靠性的前提下,盡可能地降低重量,所涉及的性能指標(biāo)包括質(zhì)量、剛度以及強(qiáng)度。通過(guò)理論推導(dǎo),可以獲得龍門下框架的質(zhì)量關(guān)于設(shè)計(jì)參數(shù)的表達(dá)式,卻很難得到其剛度及強(qiáng)度關(guān)于設(shè)計(jì)參數(shù)的表達(dá)式,而響應(yīng)面法恰恰可以有效解決這一難題。目前,響應(yīng)面法以精度高、成本低、計(jì)算方便、易于迭代等優(yōu)點(diǎn)被廣泛應(yīng)用在輕量化設(shè)計(jì)領(lǐng)域中。石磊等[2]通過(guò)拓?fù)鋬?yōu)化獲得了飛機(jī)起落架扭力臂重構(gòu)模型,在此基礎(chǔ)上基于響應(yīng)面法對(duì)其進(jìn)行了尺寸優(yōu)化設(shè)計(jì),實(shí)現(xiàn)了在保證剛度和強(qiáng)度的同時(shí),重量降低18.0%;王震虎等[3]利用響應(yīng)面法對(duì)白車身進(jìn)行了輕量化設(shè)計(jì),在扭轉(zhuǎn)剛度降低0.009%、彎曲剛度降低0.61%、一階扭轉(zhuǎn)模態(tài)降低0.49%、一階彎曲模態(tài)降低0.49%的情況下,實(shí)現(xiàn)減重1.73%;蔡德詠等[4]采用響應(yīng)面法對(duì)火箭炮發(fā)射箱箱體進(jìn)行了輕量化設(shè)計(jì),在滿足剛度、強(qiáng)度的基礎(chǔ)上,實(shí)現(xiàn)減重12%;朱兆一等[5]利用響應(yīng)面法對(duì)船舶開(kāi)孔板架結(jié)構(gòu)進(jìn)行了輕量化設(shè)計(jì),在抗屈曲能力提升31.7%的同時(shí),減重率達(dá)到了25.9%;李貝等[6]通過(guò)響應(yīng)面法對(duì)挖掘機(jī)動(dòng)臂進(jìn)行輕量化設(shè)計(jì),實(shí)現(xiàn)了在滿足剛度和強(qiáng)度的同時(shí),重量降低4.8%;伍建軍等[7]采用響應(yīng)面法對(duì)升降機(jī)滾筒線支架進(jìn)行了輕量化設(shè)計(jì),在滿足強(qiáng)度的前提下,實(shí)現(xiàn)減重42.5%;田碩[8]利用響應(yīng)面法對(duì)電動(dòng)旋耕機(jī)機(jī)架進(jìn)行了輕量化設(shè)計(jì),在滿足可靠性的前提下,實(shí)現(xiàn)減重10.4%。然而這些研究多集中在航空航天、汽車、軍工、船舶、工程機(jī)械、非標(biāo)設(shè)備、農(nóng)機(jī)等領(lǐng)域,將響應(yīng)面法應(yīng)用在換電站領(lǐng)域的研究卻鮮有報(bào)道。故本文提出了基于響應(yīng)面的龍門下框架輕量化設(shè)計(jì)方法,旨在降低龍門下框架重量,減少材料成本的同時(shí),為換電站的輕量化設(shè)計(jì)提供一種可借鑒的方法。

        本文針對(duì)龍門下框架的輕量化設(shè)計(jì)問(wèn)題,以各構(gòu)件厚度作為優(yōu)化變量,通過(guò)構(gòu)建龍門下框架各性能指標(biāo)關(guān)于厚度的二階多項(xiàng)式響應(yīng)面模型,來(lái)對(duì)其進(jìn)行基于尺寸優(yōu)化的輕量化設(shè)計(jì),并將輕量化結(jié)果與原方案進(jìn)行對(duì)比,以驗(yàn)證所提出的基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)方法的有效性。

        2 基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)流程

        基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)流程如圖1所示,首先根據(jù)實(shí)際服役情況對(duì)龍門下框架進(jìn)行有限元分析,以獲取龍門下框架的輕量化潛力信息,并根據(jù)工程經(jīng)驗(yàn)確定用于后續(xù)輕量化設(shè)計(jì)的設(shè)計(jì)變量;然后通過(guò)試驗(yàn)設(shè)計(jì)獲得樣本在設(shè)計(jì)區(qū)域內(nèi)的空間分布,基于這些樣本點(diǎn)更新有限元模型,并進(jìn)行相應(yīng)的有限元分析以獲取龍門下框架各樣本點(diǎn)的質(zhì)量、最大變形、最大等效應(yīng)力值,在此基礎(chǔ)上,構(gòu)建龍門下框架各性能響應(yīng)關(guān)于設(shè)計(jì)變量的二階多項(xiàng)式響應(yīng)面模型,并對(duì)其精度進(jìn)行檢驗(yàn);最后基于構(gòu)建的滿足精度要求的二階多項(xiàng)式響應(yīng)面模型,結(jié)合龍門下框架的輕量化需求及可靠性要求,建立基于響應(yīng)面法的龍門下框架的輕量化設(shè)計(jì)數(shù)學(xué)模型并進(jìn)行求解,之后對(duì)優(yōu)化結(jié)果進(jìn)行修正,并對(duì)修正值進(jìn)行仿真驗(yàn)證。

        3 龍門下框架靜力學(xué)分析

        龍門下框架主要由矩形管和鈑金件焊接而成,其組成部分主要包括上縱梁、下縱梁、橫梁、橫梁蓋板、加強(qiáng)筋、墊塊、固定板、連接板、護(hù)板、豎梁,其三維模型如圖2所示。

        將建立的三維模型導(dǎo)入ANSA中,設(shè)置材料為Q235鋼。根據(jù)工程經(jīng)驗(yàn),設(shè)置網(wǎng)格尺寸為10mm,網(wǎng)格盡量劃分為四邊形,在幾何特征復(fù)雜的地方使用三角形過(guò)渡,對(duì)于厚度小于等于12mm的構(gòu)件,采用殼單元進(jìn)行離散,對(duì)于厚度大于12mm的構(gòu)件,采用實(shí)體單元進(jìn)行離散,最終生成65990個(gè)單元和78498個(gè)節(jié)點(diǎn)。

        龍門下框架主要用于換電時(shí)完成車輛的舉升動(dòng)作,其上端通過(guò)4條鏈條與龍門上框架相連,下端通過(guò)四個(gè)墊塊來(lái)承載車輛重量,取龍門下框架舉升最大重量作為設(shè)計(jì)工況,最大設(shè)計(jì)載荷為5t,故在每個(gè)墊塊上分別施加12250N豎直向下的載荷,在固定板的8個(gè)螺栓孔處施加全約束。

        龍門下框架靜力學(xué)分析結(jié)果如圖3、圖4所示,由位移云圖可知,最大變形出現(xiàn)在橫梁底部護(hù)板中間位置,大小為4.6mm,小于設(shè)計(jì)要求的6mm。由應(yīng)力云圖可知,等效應(yīng)力最大值出現(xiàn)在豎梁與橫梁焊接處,大小為190.2MPa,小于材料屈服強(qiáng)度的235MPa,其他大部分應(yīng)力均小于50MPa。由此可知,龍門下框架雖然滿足剛度和強(qiáng)度要求,但設(shè)計(jì)過(guò)于保守,材料利用率并不高,具有較大的輕量化潛力。

        4 構(gòu)造龍門下框架響應(yīng)面模型

        4.1 響應(yīng)面法

        在工程實(shí)際中,我們非常希望能夠獲得所關(guān)注的性能指標(biāo)關(guān)于各設(shè)計(jì)參數(shù)間的解析表達(dá)式,從而可以高效開(kāi)展定性分析、定量分析以及優(yōu)化設(shè)計(jì)等工作。然而,性能指標(biāo)與設(shè)計(jì)參數(shù)的關(guān)系往往呈隱式非線性關(guān)系,無(wú)法通過(guò)理論推導(dǎo)得到,在此背景下,能夠快速、準(zhǔn)確地構(gòu)造性能指標(biāo)關(guān)于設(shè)計(jì)參數(shù)的表達(dá)式的響應(yīng)面法得到了快速發(fā)展。

        響應(yīng)面法是通過(guò)將試驗(yàn)設(shè)計(jì)與數(shù)理統(tǒng)計(jì)相結(jié)合,構(gòu)造出響應(yīng)關(guān)于設(shè)計(jì)變量的解析表達(dá)式,以此來(lái)代替響應(yīng)關(guān)于設(shè)計(jì)變量的實(shí)際關(guān)系的一種方法[9]。為了獲取龍門下框架結(jié)構(gòu)性能與設(shè)計(jì)變量之間的關(guān)系,本文采用二階多項(xiàng)式響應(yīng)面代理模型來(lái)構(gòu)造龍門下框架質(zhì)量、最大變形、最大等效應(yīng)力的響應(yīng)面模型。

        二階多項(xiàng)式響應(yīng)面模型可表示為:

        (1)

        式中,、、、為回歸系數(shù);為設(shè)計(jì)變量個(gè)數(shù);為誤差。

        4.2 設(shè)計(jì)變量的選取

        根據(jù)工程實(shí)際經(jīng)驗(yàn),選取上縱梁厚度、下縱梁厚度、橫梁厚度、橫梁蓋板厚度作為設(shè)計(jì)變量,其中,上縱梁、下縱梁、橫梁為矩形管,橫梁蓋板為鈑金件。根據(jù)GB/T 3094-2012、GB/T 6728-2017以及某司常用的矩形管、鈑金件厚度規(guī)格確定上縱梁厚度、下縱梁厚度、橫梁厚度、橫梁蓋板厚度的取值規(guī)格及范圍,其如表1所示。

        4.3 試驗(yàn)設(shè)計(jì)

        要想得到精確的優(yōu)化結(jié)果,則需要嚴(yán)格控制所構(gòu)建的響應(yīng)面模型的精度,而樣本點(diǎn)的數(shù)量及空間位置很大程度上決定了所建立的響應(yīng)面模型的精度,因此,在構(gòu)造響應(yīng)面模型的過(guò)程中,試驗(yàn)設(shè)計(jì)方法的選擇至關(guān)重要。由于中心復(fù)合設(shè)計(jì)是根據(jù)二次多項(xiàng)式的特點(diǎn)來(lái)構(gòu)造的,能以最少的試驗(yàn)次數(shù)提供關(guān)于設(shè)計(jì)變量和試驗(yàn)誤差的重要信息[10],且具有設(shè)計(jì)簡(jiǎn)單、預(yù)測(cè)性好、建模效率高等優(yōu)點(diǎn)[11],因此,本文采用中心復(fù)合試驗(yàn)設(shè)計(jì)方法在龍門下框架的設(shè)計(jì)空間中進(jìn)行樣本點(diǎn)的采集,基于這些樣本點(diǎn)更新有限元模型,并進(jìn)行相應(yīng)的有限元分析,從而獲得龍門下框架各樣本點(diǎn)的質(zhì)量、最大變形、最大等效應(yīng)力值,試驗(yàn)數(shù)據(jù)見(jiàn)表2。

        4.4 響應(yīng)面模型的建立

        為建立反映龍門下框架各性能指標(biāo)與設(shè)計(jì)變量關(guān)系的二階多項(xiàng)式響應(yīng)面模型,利用得到的25組數(shù)據(jù),采用最小二乘法對(duì)質(zhì)量、最大變形、最大等效應(yīng)力進(jìn)行擬合。對(duì)于本文的4個(gè)設(shè)計(jì)變量而言,響應(yīng)面的形式如果采用完全二次,則由式(1)可知,總項(xiàng)數(shù)為15。而項(xiàng)數(shù)過(guò)多會(huì)導(dǎo)致所構(gòu)建的響應(yīng)面模型較為復(fù)雜且預(yù)測(cè)精度不高,從而影響后續(xù)的優(yōu)化效率和求解精度[12]。因此,本文在構(gòu)建龍門下框架各響應(yīng)關(guān)于設(shè)計(jì)變量的二階響應(yīng)面模型的過(guò)程中,對(duì)各項(xiàng)的顯著性進(jìn)行檢驗(yàn),保留顯著項(xiàng),刪除不顯著項(xiàng),最終構(gòu)造的龍門下框架質(zhì)量、最大變形和最大等效應(yīng)力響應(yīng)面模型分別為:

        (2)

        (3)

        (4)

        4.5 精度檢驗(yàn)

        為了驗(yàn)證所構(gòu)建的龍門下框架各性能指標(biāo)響應(yīng)面模型是否滿足工程精度要求,需要對(duì)其進(jìn)行精度檢驗(yàn)。針對(duì)響應(yīng)面模型的精度檢驗(yàn)問(wèn)題,本文選擇調(diào)整的確定系數(shù)作為評(píng)價(jià)指標(biāo),一般來(lái)說(shuō),越接近于1,則響應(yīng)面模型精度越高,越接近于0,則精度越差??杀硎緸椋?/p>

        (5)

        式中:為仿真值;為擬合值;為仿真值均值;為樣本點(diǎn)的個(gè)數(shù);為響應(yīng)面模型的總項(xiàng)數(shù)。

        通過(guò)式(5)計(jì)算得到龍門下框架質(zhì)量、最大變形、最大等效應(yīng)力的調(diào)整的確定系數(shù)分別為0.9999、0.9987、0.9979,由此可知所構(gòu)建的龍門下框架各性能指標(biāo)響應(yīng)面模型均滿足精度要求。

        5 龍門下框架輕量化設(shè)計(jì)

        綜合考慮龍門下框架的輕量化需求及可靠性要求,基于構(gòu)建的滿足精度要求的質(zhì)量、最大變形以及最大等效應(yīng)力響應(yīng)面模型,以上縱梁厚度、下縱梁厚度、橫梁厚度、橫梁蓋板厚度為設(shè)計(jì)變量,以質(zhì)量極小化為目標(biāo),以最大變形小于等于6mm和最大等效應(yīng)力小于等于235MPa為約束條件,建立基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)數(shù)學(xué)模型如下:

        (6)

        采用尋優(yōu)收斂迅速且能簡(jiǎn)化計(jì)算復(fù)雜度的序列二次規(guī)劃算法[13]對(duì)龍門下框架輕量化數(shù)學(xué)模型進(jìn)行優(yōu)化求解,根據(jù)表1中設(shè)計(jì)變量的取值規(guī)格對(duì)優(yōu)化解進(jìn)行修正,基于修正結(jié)果更新有限元模型并重新進(jìn)行靜力學(xué)分析,得到的位移云圖及等效應(yīng)力云圖分別如圖5、圖6所示,輕量化設(shè)計(jì)前、后設(shè)計(jì)變量及性能指標(biāo)對(duì)比見(jiàn)表3。

        由表3可知,龍門下框架經(jīng)過(guò)輕量化設(shè)計(jì)后,從剛度來(lái)看:最大變形由優(yōu)化前的4.6mm增至5.8mm,雖有一定程度增大,但小于設(shè)計(jì)要求的6mm,滿足剛度要求;從強(qiáng)度來(lái)看:等效應(yīng)力最大值由優(yōu)化前的190.2MPa增至229.2MPa,雖有所增大,但小于材料屈服強(qiáng)度的235MPa,滿足強(qiáng)度要求;從輕量化效果來(lái)看:質(zhì)量由優(yōu)化前的432.1kg降為350.0kg,降低了82.1kg,減重率為19.0%,輕量化效果顯著。

        6 結(jié)論

        本文基于響應(yīng)面法對(duì)龍門下框架進(jìn)行了輕量化設(shè)計(jì)。構(gòu)建了龍門下框架質(zhì)量、最大變形以及最大等效應(yīng)力關(guān)于各構(gòu)件厚度的二階多項(xiàng)式響應(yīng)面模型,并驗(yàn)證了其具有足夠的精度,在此基礎(chǔ)上,以質(zhì)量為目標(biāo),以剛度及強(qiáng)度為約束,對(duì)各構(gòu)件的厚度進(jìn)行了優(yōu)化。結(jié)果表明,在滿足剛度及強(qiáng)度要求的前提下,質(zhì)量降低了82.1kg,減重率達(dá)到了19.0%,輕量化效果顯著,驗(yàn)證了所提出的基于響應(yīng)面法的龍門下框架輕量化設(shè)計(jì)方法的有效性,為換電站相關(guān)總成及子系統(tǒng)的輕量化設(shè)計(jì)提供了一種可借鑒的方法。

        參考文獻(xiàn):

        [1]馮淵,賀宏勝,李曉強(qiáng),等.純電動(dòng)汽車充換電站規(guī)模設(shè)計(jì)仿真研究[J].上海汽車,2012,259(03):7-9.

        [2]石磊,常嘯.基于拓?fù)渑c尺寸優(yōu)化的某型飛機(jī)起落架扭力臂輕量化設(shè)計(jì)[J].機(jī)電工程技術(shù),2021,50(02):70-75.

        [3]王震虎,周巧英,劉開(kāi)勇,等.基于響應(yīng)面模型的白車身多目標(biāo)輕量化設(shè)計(jì)[J].中國(guó)機(jī)械工程,2018,29(01):75-81.

        [4]蔡德詠,劉福軍,田柯文,等.火箭炮發(fā)射箱可靠性分析及輕量化設(shè)計(jì)[J].機(jī)械設(shè)計(jì),2017,34(05):81-85.

        [5]朱兆一,李曉文,李釗,等.基于響應(yīng)面法復(fù)合材料船舶開(kāi)孔板架優(yōu)化設(shè)計(jì)[J].艦船科學(xué)技術(shù),2022,44(09):48-52.

        [6]李貝,張倩.基于響應(yīng)面法的挖掘機(jī)動(dòng)臂優(yōu)化設(shè)計(jì)[J].工程機(jī)械,2021,52(05):79-85+11.

        [7]伍建軍,王子寧,任崇軒,等.基于響應(yīng)面法的升降機(jī)滾筒線支架輕量化設(shè)計(jì)[J].機(jī)械設(shè)計(jì)與制造,2020,354(08):263-266.

        [8]田碩. 新型整體式電動(dòng)旋耕機(jī)結(jié)構(gòu)設(shè)計(jì)及輕量化研究[D]. 成都: 電子科技大學(xué),2021.

        [9]Mostaani S, Singh D,F(xiàn)irouzbakhsh K,et al. Optimization of a passive vehicle suspension system for ride comfort enhancement with different speeds based on DOE method[J]. International Journal on Manufacturing & Material Science, 2011, 1(2): 5-10.

        [10]朱勝陽(yáng). 高速鐵路無(wú)砟軌道結(jié)構(gòu)傷損行為及其對(duì)動(dòng)態(tài)性能的影響[D]. 成都: 西南交通大學(xué), 2015.

        [11]賴宇陽(yáng). Isight參數(shù)優(yōu)化理論與實(shí)例詳解[M]. 北京: 北京航空航天大學(xué)出版社, 2012.

        [12]王子寧. 考慮穩(wěn)健性的滾筒線支架輕量化設(shè)計(jì)研究[D]. 贛州: 江西理工大學(xué),2019.

        [13]伍建軍,王子寧,金渶棋,等. 基于6穩(wěn)健性的運(yùn)輸工裝輕量化設(shè)計(jì)[J].機(jī)械設(shè)計(jì)與制造,2020,353(07):273-277.

        猜你喜歡
        輕量化設(shè)計(jì)響應(yīng)面法
        山藥多糖提取工藝優(yōu)化及其抗菌活性研究
        微波輔助提取白頭翁皂苷研究
        響應(yīng)面法優(yōu)化超聲輔助提取蕎麥中蘆丁的工藝
        菊芋膳食纖維酸奶的工藝研究及營(yíng)養(yǎng)分析
        軌道扣件彈性墊板結(jié)構(gòu)優(yōu)化設(shè)計(jì)
        船舶舷側(cè)構(gòu)件幾何尺寸優(yōu)化分析
        科技視界(2016年21期)2016-10-17 17:22:54
        乘用車車身結(jié)構(gòu)輕量化設(shè)計(jì)技術(shù)研究與實(shí)踐
        帽形截面薄壁焊接結(jié)構(gòu)軸向耐撞性的優(yōu)化設(shè)計(jì)方法研究
        飛機(jī)平尾普通肋的優(yōu)化設(shè)計(jì)
        汽車高強(qiáng)鋼大通風(fēng)孔車輪輪輻成形工藝優(yōu)化仿真
        成人免费a级毛片无码片2022| 97人妻蜜臀中文字幕| 日韩人妻一区二区中文字幕| 亚洲综合精品中文字幕| 国产农村妇女毛片精品久久| 亚洲区在线| 久久久精品人妻一区二区三区四| 亚洲gv白嫩小受在线观看| 国产精品麻豆A啊在线观看| 亚洲国产一区二区视频| 97一期涩涩97片久久久久久久 | 国产色在线 | 日韩| 无码av无码天堂资源网| 日本在线视频网站www色下载| 国产激情一区二区三区成人免费| 永久免费看黄网站性色| 日韩国产人妻一区二区三区| 性动态图av无码专区| 国产成人精品免费视频大全| 蜜桃夜夜爽天天爽三区麻豆av| 亚洲av午夜一区二区三| 久久99精品国产麻豆宅宅| 精品无码AV无码免费专区| 国产三级在线观看不卡| 华人在线视频精品在线| 一本色道久久88精品综合 | 成年性生交大片免费看| 99精品国产一区二区三区a片| 久久九九有精品国产尤物| 国产精品丝袜美腿诱惑| 少妇性l交大片7724com| 秋霞鲁丝片av无码| 国产精品美女久久久浪潮av| 日韩人妻系列在线观看| 男女上下猛烈啪啪免费看| 国产欧美日韩综合一区二区三区 | 亚洲第一幕一区二区三区在线观看| 超碰97资源站| 国产在线视欧美亚综合| 手机在线免费观看的av| 强奷乱码中文字幕|