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

        ?

        車身零件預(yù)裝變形有限元虛擬評(píng)估分析

        2021-02-18 19:44:01程計(jì)棟鄧?yán)^濤石文
        計(jì)算機(jī)輔助工程 2021年4期
        關(guān)鍵詞:有限元

        程計(jì)棟 鄧?yán)^濤 石文

        摘要:? 在車身工藝同步工程階段,通過(guò)有限元仿真方法對(duì)車身制造過(guò)程中預(yù)裝零件的不同產(chǎn)品設(shè)計(jì)方案進(jìn)行虛擬評(píng)估分析,根據(jù)分析結(jié)果指導(dǎo)零件設(shè)計(jì)優(yōu)化。該仿真方法在生產(chǎn)實(shí)物階段得到驗(yàn)證,可提前識(shí)別制造過(guò)程質(zhì)量風(fēng)險(xiǎn)。與傳統(tǒng)的通過(guò)經(jīng)驗(yàn)評(píng)估判斷零件預(yù)裝變形的方法相比,在工藝同步工程中應(yīng)用有限元仿真技術(shù),為解決車身制造過(guò)程中的復(fù)雜工況變形問(wèn)題提供理論分析方法和依據(jù),可靠性高,并達(dá)到縮短制造周期、降低制造成本的目的。

        關(guān)鍵詞:? 車身制造; 零件預(yù)裝; 后端板; 同步工程; 有限元

        中圖分類號(hào):? U461;TB125文獻(xiàn)標(biāo)志碼:? B

        Car body part pre-install deformation finite element

        virtual assessment analysis

        CHENG Jidong, DENG Jitao, SHI Wen

        (Vehicle Manufacturing Department, SAIC-General Motors Co., Ltd., Shanghai 201201, China)

        Abstract: Finite element virtual assessment did for car body part pre-install deformation based on different product design proposals during process simultaneous engineering. The method that validated during car body pre-build can instruct product design and identify manufacturing quality risk in advance. Compared with traditional method that identify part preinstall deformation by experience, finite element simulation method can provide theoretical evidence to solve complicated manufacturing deformation issue in body process simultaneous engineering, and the reliability is high. What’s more, it is helpful to shorten manufacturing time and reduce cost.

        Key words: car body manufacturing; pre-install of parts; rear end panel; simultaneous engineering; finite element

        收稿日期:? 2021-11-09修回日期:? 2021-11-23

        作者簡(jiǎn)介: 程計(jì)棟(1981—),男,江蘇徐州人,高級(jí)工程師,研究方向?yàn)檐嚿砬捌跀?shù)字化規(guī)劃,(E-mail)Jidong_Cheng@saic-gm.com

        通信作者:鄧?yán)^濤(1991—),男,江西安義人,工程師,碩士,研究方向?yàn)檐嚿碇圃旃に?,(E-mail)Jitao_Deng@saic-gm.com0引言高速發(fā)展的汽車行業(yè)競(jìng)爭(zhēng)非常激烈,各大車企在不斷縮短產(chǎn)品研發(fā)周期、加快新車型導(dǎo)入的同時(shí),也在不斷通過(guò)技術(shù)創(chuàng)新尋求降低制造成本、提升制造質(zhì)量的突破。同步工程指從新產(chǎn)品開(kāi)發(fā)的最初階段開(kāi)始,就盡可能把設(shè)計(jì)與制造聯(lián)系起來(lái),要求整個(gè)產(chǎn)品開(kāi)發(fā)過(guò)程實(shí)施同步、一體化設(shè)計(jì),促使開(kāi)發(fā)者始終考慮從概念形成到用后處置的整個(gè)產(chǎn)品生命周期內(nèi)所有因素(包括工藝、質(zhì)量、成本、進(jìn)度和用戶要求)的一種系統(tǒng)方法[1]。同步工程技術(shù)在汽車開(kāi)發(fā)過(guò)程中已被廣泛應(yīng)用,如文獻(xiàn)[2-5]分別介紹同步工程在沖壓、車身、油漆、總裝工藝開(kāi)發(fā)中的應(yīng)用。據(jù)統(tǒng)計(jì),采用同步工程后,車企在開(kāi)發(fā)新產(chǎn)品時(shí)設(shè)計(jì)變更數(shù)量可減少近50%,開(kāi)發(fā)周期可縮短40%~60%,開(kāi)發(fā)費(fèi)用可降低30%~40%[6]。因此,同步工程技術(shù)的開(kāi)發(fā)和應(yīng)用,對(duì)于整車質(zhì)量風(fēng)險(xiǎn)識(shí)別,降低制造成本具有重要意義。

        1制造工況的有限元虛擬仿真有限元仿真技術(shù)已廣泛應(yīng)用在車身性能分析中,如對(duì)車身剛度、強(qiáng)度、碰撞等的分析[7-8],在沖壓、油漆工藝同步工程中也有一定的應(yīng)用,如鈑金零件的成形性分析[9]、油漆電泳膜厚度分析[10]等。車身制造是個(gè)復(fù)雜過(guò)程,包含上料、抓料、預(yù)裝、拼合、焊接、落位、運(yùn)輸、裝夾等各種工藝,工藝過(guò)程中的零件變形預(yù)測(cè)和控制是車身制造質(zhì)量保證的關(guān)鍵因素。車身性能的虛擬評(píng)估并不涉及制造工況,如車身模態(tài)、剛度和強(qiáng)度分析等,因此無(wú)法預(yù)測(cè)車身制造過(guò)程中零件變形造成的工藝失效等質(zhì)量風(fēng)險(xiǎn)。傳統(tǒng)的依靠經(jīng)驗(yàn)的評(píng)估方法可靠性不高,往往在產(chǎn)品實(shí)物驗(yàn)證階段才會(huì)發(fā)現(xiàn)問(wèn)題,造成制造成本的增加。因此,車身制造迫切需要在產(chǎn)品開(kāi)發(fā)早期,將虛擬仿真技術(shù)從產(chǎn)品性能虛擬評(píng)估拓展到制造過(guò)程評(píng)估,對(duì)復(fù)雜工況下零件變形、應(yīng)力狀態(tài)等進(jìn)行準(zhǔn)確的分析和辨識(shí),從而對(duì)產(chǎn)品設(shè)計(jì)進(jìn)行早期優(yōu)化和指導(dǎo),避免制造質(zhì)量問(wèn)題遺留到試生產(chǎn)或量產(chǎn)造車階段,以便提升整車制造質(zhì)量、縮短產(chǎn)品開(kāi)發(fā)周期、節(jié)約制造資源、降低開(kāi)發(fā)和制造成本。隨著計(jì)算機(jī)的發(fā)展,有限元仿真技術(shù)在工藝同步工程中發(fā)揮越來(lái)越重要的作用:陳丹丹等[11]從總體角度介紹有限元仿真技術(shù)在車身工藝同步工程中的各種應(yīng)用案例;趙偉鵬等[12]在工藝同步工程中采用有限元方法對(duì)汽車副儀表板進(jìn)行模流分析,指導(dǎo)產(chǎn)品拔模設(shè)計(jì),避免后期設(shè)計(jì)變更引起質(zhì)量、成本等問(wèn)題。目前,采用有限元方法對(duì)車身焊裝制造過(guò)程中鈑金的變形進(jìn)行分析預(yù)測(cè),并根據(jù)分析結(jié)果指導(dǎo)產(chǎn)品設(shè)計(jì)或工藝開(kāi)發(fā)的研究仍然較少。本文通過(guò)在車身工藝同步工程中應(yīng)用有限元虛擬仿真技術(shù),采用LS-DYNA有限元仿真軟件對(duì)車身制造過(guò)程中的零件預(yù)裝進(jìn)行動(dòng)態(tài)虛擬仿真分析,闡述車身零件預(yù)裝的具體分析方法和流程,為指導(dǎo)產(chǎn)品設(shè)計(jì)、降低制造成本和風(fēng)險(xiǎn)、提高制造質(zhì)量提供理論依據(jù)。

        2零件預(yù)裝車身制造是按照一定的裝配順序和尺寸定位要求,把單個(gè)沖壓零件通過(guò)電阻焊、弧焊、激光焊、螺栓、沖鉚、折邊、膠接等連接方式拼成一個(gè)完整的車體結(jié)構(gòu)的過(guò)程。零件預(yù)裝是車身制造過(guò)程常見(jiàn)的工藝形式,將車身零件或總成進(jìn)行預(yù)裝,然后在后道工藝進(jìn)行焊接,不僅提高工藝規(guī)劃的可操作性和機(jī)動(dòng)性,也適應(yīng)現(xiàn)階段車身制造高節(jié)拍、多車型混線生產(chǎn)的需求。常見(jiàn)的車身預(yù)裝零件,如側(cè)圍總成、車頂橫梁、后端板等見(jiàn)圖1。

        為適應(yīng)高節(jié)拍的生產(chǎn)需求,目前車身車間普遍采用高速滾床進(jìn)行零件或總成的機(jī)運(yùn)傳輸,零件在傳輸過(guò)程中會(huì)經(jīng)歷加速、減速甚至急停等工況,因此預(yù)裝零件在運(yùn)輸過(guò)程中存在變形、裝配失效甚至掉落的風(fēng)險(xiǎn)。尤其是后端板總成,其裝配方向與運(yùn)輸方向一致,在運(yùn)輸過(guò)程中承受較大的加速啟動(dòng)、減速急停過(guò)程,具有較高的運(yùn)輸風(fēng)險(xiǎn)。如圖2所示,在產(chǎn)品驗(yàn)證的試生產(chǎn)階段,某車型后端板預(yù)裝后在運(yùn)輸急停時(shí)掉落,不僅影響產(chǎn)品驗(yàn)證周期,而且會(huì)增加產(chǎn)品后期更改時(shí)的修模、物料、調(diào)試等制造成本。因此,在產(chǎn)品開(kāi)發(fā)早期階段,通過(guò)有限元虛擬仿真方法評(píng)估制造風(fēng)險(xiǎn)工況,并根據(jù)分析結(jié)果指導(dǎo)產(chǎn)品結(jié)構(gòu)優(yōu)化,對(duì)提升產(chǎn)品制造的可靠性、降低制造成本具有重要意義。

        后端板通過(guò)搭扣進(jìn)行預(yù)裝(見(jiàn)圖3),為評(píng)估零件預(yù)裝是否滿足車身制造要求,在車身前期工藝同步工程評(píng)估過(guò)程中,通過(guò)有限元仿真方法對(duì)零件預(yù)裝的有效性進(jìn)行評(píng)估,可以指導(dǎo)產(chǎn)品設(shè)計(jì),避免前期設(shè)計(jì)問(wèn)題遺留到制造階段,節(jié)約潛在的產(chǎn)品更改引起修模、調(diào)試等費(fèi)用,縮短制造周期。

        3預(yù)裝仿真流程零件預(yù)裝有限元仿真需要分析零件在傳輸過(guò)程中加速、急停等動(dòng)態(tài)工況下的運(yùn)動(dòng),因此可采用適合非線性動(dòng)態(tài)的有限仿真軟件LS-DYNA進(jìn)行分析求解,零件預(yù)裝有限元仿真流程見(jiàn)圖4。

        在車身工藝同步工程中,獲取產(chǎn)品性能分析的有限元網(wǎng)格模型并用于零件預(yù)裝仿真分析,既可以提高工作效率,又是產(chǎn)品與工藝協(xié)同的過(guò)程。獲取產(chǎn)品有限元模型后,由于零件存在設(shè)計(jì)和制造2種狀態(tài),需要在有限元前處理軟件HyperMesh中根據(jù)制造仿真工況進(jìn)行模型處理,如后端板搭扣在設(shè)計(jì)狀態(tài)并未扣上后端板,但在制造狀態(tài)下已扣上搭扣,因此需要根據(jù)制造狀態(tài)對(duì)模型進(jìn)行修改,將需要評(píng)估的后端板搭扣扣上。在完成制造工況模型處理后,需對(duì)模型進(jìn)行邊界條件的輸入,包括零件傳輸過(guò)程中底板運(yùn)輸孔的約束,設(shè)置運(yùn)輸過(guò)程的速度曲線或瞬時(shí)加速度等,以及模型的自接觸和車身搭扣與后端板的接觸等。在定義完接觸等相關(guān)的計(jì)算控制參數(shù)后,提交軟件進(jìn)行求解計(jì)算。提取求解及處理完畢的結(jié)果,可以查看零件運(yùn)輸過(guò)程的變形量及應(yīng)力、應(yīng)變狀態(tài),同時(shí)可結(jié)合變形動(dòng)畫等評(píng)估預(yù)裝零件在運(yùn)輸過(guò)程中是否存在掉落風(fēng)險(xiǎn),從而確認(rèn)產(chǎn)品設(shè)計(jì)是否滿足制造過(guò)程的要求。如果設(shè)計(jì)不滿足要求,可以根據(jù)仿真結(jié)果提出產(chǎn)品改進(jìn)意見(jiàn)、指導(dǎo)產(chǎn)品更改,并對(duì)更改后的產(chǎn)品進(jìn)行校核,直至滿足制造要求。對(duì)于簡(jiǎn)單的產(chǎn)品更改(如搭扣尺寸更改),也可通過(guò)網(wǎng)格的編輯直接對(duì)有限元模型進(jìn)行修改,從而快速完成新設(shè)計(jì)的校核,避免產(chǎn)品更改后再次進(jìn)行有限元建模,以便縮短評(píng)估周期、提高工作效率。

        4仿真模型及方案后端板預(yù)裝仿真采用CAE方法分析后端板搭扣設(shè)計(jì)是否滿足車身運(yùn)輸過(guò)程的要求,由于在該制造工況下零件的變形主要發(fā)生在后端板及搭扣,因此使用車身后部截?cái)嗄P瓦M(jìn)行仿真分析。后部截?cái)嗄P涂稍诒WC該工況計(jì)算精度的前提下,大大提高計(jì)算效率、節(jié)約計(jì)算資源。通過(guò)有限元軟件HyperMesh進(jìn)行前處理:采用平均大小為5 mm的網(wǎng)格尺寸,對(duì)關(guān)注的搭扣處網(wǎng)格細(xì)化至平均大小為2 mm,圓角處網(wǎng)格細(xì)化至平均大小為1 mm。模型總體單元數(shù)量為77 028個(gè),其中三角形單元2 776個(gè),占網(wǎng)格總數(shù)的3.6%。網(wǎng)格質(zhì)量采用如下標(biāo)準(zhǔn):長(zhǎng)寬比小于5,雅可比系數(shù)大于0.7,三角形單元角度范圍30°~120°,四邊形單元角度范圍45°~135°,最終建立的有限元仿真網(wǎng)格模型見(jiàn)圖5。

        前處理過(guò)程還包括零件厚度、焊點(diǎn)連接、材料模型、載荷設(shè)置等過(guò)程等。根據(jù)零件的厚度,采用section_shell關(guān)鍵字賦予零件板厚屬性;根據(jù)鈑金零件的材料特性,采用彈塑性材料模型賦予零件材料信息,包含彈性模量、泊松比、材料密度、屈服極限以及塑性階段材料的應(yīng)力-應(yīng)變曲線;采用實(shí)體單元建立焊點(diǎn)模型,并通過(guò)contact_spoweld關(guān)鍵字建立焊點(diǎn)與鈑金的接觸,通過(guò)load_body關(guān)鍵字模擬重力慣性載荷和機(jī)運(yùn)急停慣性載荷。同時(shí),后端板與車身之間設(shè)置面-面接觸,模擬預(yù)裝過(guò)程零件載荷的傳遞。

        后端板搭扣越多,后端板預(yù)裝越可靠,但搭扣的增加也增加產(chǎn)品的成本及生產(chǎn)操作的時(shí)間。為平衡產(chǎn)品制造、成本和預(yù)裝的可靠性,需要合理地選擇搭扣的數(shù)量、位置、尺寸等。根據(jù)前期產(chǎn)品和工藝的協(xié)同,最終確定后端板單側(cè)4個(gè)搭扣作為選擇方案(見(jiàn)圖5中的A、B、C和D),單側(cè)選擇2個(gè)搭扣作為最終方案。根據(jù)搭扣分布,上部和下部各選擇1個(gè)搭扣,最終選擇表1中的4種方案進(jìn)行后端板搭扣評(píng)估。在有限元模型設(shè)置過(guò)程中,選擇的2個(gè)搭扣應(yīng)扣上后端板,另2個(gè)搭扣應(yīng)處于放開(kāi)狀態(tài),圖5即為方案4所示模型,即搭扣B和C扣上后端板,搭扣A和D放開(kāi)。在有限元模型計(jì)算設(shè)置時(shí),扣緊的搭扣與后端板需要設(shè)置接觸及相關(guān)的接觸控制參數(shù)。

        5仿真結(jié)果及分析根據(jù)上述方案,對(duì)后端板預(yù)裝有限元模型分別進(jìn)行處理,并根據(jù)該產(chǎn)品制造工廠的運(yùn)輸條件進(jìn)行邊界條件的輸入,最后輸入LS-DYNA軟件進(jìn)行仿真計(jì)算,得出仿真結(jié)果。在HyperView中打開(kāi)d3plot文件,可提取出仿真結(jié)果,如零件的應(yīng)力云圖、變形云圖等。本文以方案1為例,對(duì)其仿真結(jié)果進(jìn)行分析。方案1的應(yīng)力云圖見(jiàn)圖6,可知搭扣A處零件的最大應(yīng)力已達(dá)到464.7 MPa,已超出材料的屈服極限340 MPa。方案1的變形云圖見(jiàn)圖7,可知后端的最大變形量達(dá)到67.56 mm,且搭扣已與后端板脫離。因此,根據(jù)零件的應(yīng)力及變形狀態(tài)可知,后端板在搭扣預(yù)裝運(yùn)輸時(shí)已掉落,不滿足制造要求。后端板搭扣各方案仿真結(jié)果見(jiàn)表2,可知:方案2和方案4滿足要要求,方案1和方案3不滿足要求。因此,為滿足制造要求必須選擇具有搭扣B的方案,搭扣C和D可二選一,不能選擇有搭扣A的方案。

        上部搭扣B對(duì)后端板的預(yù)裝效果更明顯,搭扣A雖然也在上部,但其處于斜面位置,與車身運(yùn)輸方向存在角度,故對(duì)于后端板的固定作用并不明顯,這在后續(xù)的生產(chǎn)實(shí)物驗(yàn)證中也得到驗(yàn)證,見(jiàn)圖8。

        綜合上述分析,通過(guò)對(duì)后端板預(yù)裝的有限元仿真,可指導(dǎo)搭扣位置和數(shù)量的設(shè)計(jì)。同時(shí),該方法可通過(guò)改變搭扣的尺寸,評(píng)估不同尺寸狀態(tài)下的搭扣是否滿足預(yù)裝要求,實(shí)現(xiàn)指導(dǎo)搭扣尺寸設(shè)計(jì)的目的。

        6結(jié)束語(yǔ)闡述車身零件預(yù)裝變形的有限元仿真方法和流程,在車身工藝同步工程階段應(yīng)用有限元仿真技術(shù),對(duì)后端板預(yù)裝搭扣進(jìn)行虛擬評(píng)估,指導(dǎo)搭扣數(shù)量和布置設(shè)計(jì),也可以指導(dǎo)搭扣尺寸的優(yōu)化。通過(guò)在車身工藝同步工程階段對(duì)車身制造工況進(jìn)行有限元虛擬評(píng)估,可以在產(chǎn)品開(kāi)發(fā)早期發(fā)現(xiàn)和解決車身制造問(wèn)題,提前識(shí)別制造質(zhì)量風(fēng)險(xiǎn),對(duì)于指導(dǎo)車身零件設(shè)計(jì),縮短車身制造周期,降低車身制造成本具有重要意義。參考文獻(xiàn):

        [1]黨紅霞. 焊裝SE分析在同步工程中的應(yīng)用[J]. 焊接技術(shù), 2011, 40(4): 24-28. DOI: 10.13846/j.cnki.cn12-1070/tg.2011.04.004.

        [2]羅仁平, 王燕, 趙金龍. 沖壓同步工程在汽車產(chǎn)品設(shè)計(jì)開(kāi)發(fā)中的應(yīng)用[J]. 汽車技術(shù), 2014(7): 5-7.

        [3]黃智. 汽車白車身制造工藝同步工程淺析[J]. 汽車工藝與材料, 2018(9): 49-54. DOI: 10.19710/J.cnki.1003-8817.20180139.

        [4]李欣, 鄔波, 易煒. 汽車涂裝同步工程概述[J]. 電鍍與涂飾, 2019, 38(15): 812-816. DOI: 10.19289/j.1004-227x.2019.15.010.

        [5]李成瑤. 同步工程之總裝工藝應(yīng)用[J]. 工業(yè)技術(shù), 2018(4): 75-76. DOI: 10.13612/j.cnki.cntp.2018.04.044.

        [6]李黎. 同步工程在長(zhǎng)安汽車開(kāi)發(fā)中的運(yùn)用[J]. 高科技與產(chǎn)業(yè)化, 2006(S2): 68-69.

        [7]郭毅, 劉子建, 秦歡. 基于碰撞性能的車身剛度鏈設(shè)計(jì)方法[J]. 汽車工程, 2018, 40(12):1426-1434. DOI: 10.19562/j.chinasae.qcgc.2018.012.008.

        [8]張立強(qiáng). 基于正碰CAE分析的車身縱梁結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 計(jì)算機(jī)輔助工程, 2018, 27(3): 69-72. DOI: 10.13340/j.cae.2018.03.015.

        [9]趙鴻鵠, 盧橋, 劉春雨, 等. 沖壓同步工程在某重型車項(xiàng)目的應(yīng)用[J]. 鍛壓裝備與制造技術(shù), 2019, 54(4): 77-82. DOI: 10.16316/j.issn.1672-0121.2019.04.018.

        [10]呂金濤, 陳夢(mèng)男, 李文碰, 等. 皮卡新車型涂裝同步工程(SE)分析方法研究[J]. 電鍍與涂飾, 2019, 38(12): 599-604. DOI: 10.19289/j.1004-227x.2019.12.007.

        [11]陳丹丹, 鄧?yán)^濤, 陳錫榮, 等. 有限元仿真技術(shù)在車身工藝同步工程中的應(yīng)用[J]. 制造業(yè)自動(dòng)化, 2017, 39(8): 69-73.

        [12]趙偉鵬, 朱凱慶, 趙宇馳. 同步工程在自主汽車開(kāi)發(fā)階段的應(yīng)用[J]. 汽車實(shí)用技術(shù), 2015(11): 144-146. DOI: 10.16638/j.cnki.1671-7988.2015.11.050.(編輯陳鋒杰)

        猜你喜歡
        有限元
        基于擴(kuò)展有限元的疲勞裂紋擴(kuò)展分析
        非線性感應(yīng)加熱問(wèn)題的全離散有限元方法
        TDDH型停車器制動(dòng)過(guò)程有限元分析
        新型有機(jī)玻璃在站臺(tái)門的應(yīng)用及有限元分析
        基于I-DEAS的履帶起重機(jī)主機(jī)有限元計(jì)算
        基于有限元模型對(duì)踝模擬扭傷機(jī)制的探討
        10MN快鍛液壓機(jī)有限元分析
        磨削淬硬殘余應(yīng)力的有限元分析
        基于SolidWorks的吸嘴支撐臂有限元分析
        箱形孔軋制的有限元模擬
        上海金屬(2013年4期)2013-12-20 07:57:18
        日韩一区二区三区久久精品| 91久久精品国产性色tv| 天码av无码一区二区三区四区 | 伊人精品在线观看| 连续高潮喷水无码| 精品久久亚洲一级α| 亚洲国产精品免费一区| 亚洲精品中文字幕导航| av大全亚洲一区二区三区| 精品国偷自产在线视频九色 | 一本色道久久88—综合亚洲精品 | 在线a免费观看| 日韩在线视精品在亚洲| 东京热无码人妻中文字幕| 国产在线欧美日韩精品一区二区| 日产乱码一区二区国产内射| 色噜噜色哟哟一区二区三区| 成人影院在线观看视频免费 | 久草视频福利| 国产亚洲精品性爱视频| 亚洲一区二区三区一区| 三级国产高清在线观看| 久久久国产精品123| 欧美熟妇另类久久久久久不卡 | 一区二区无码中出| 国产成人福利在线视频不卡| 亚洲av性色精品国产| 在线观看午夜视频一区二区| 欧美猛男军警gay自慰| 少妇的丰满3中文字幕| 精品国产91久久综合| 亚洲av网一区二区三区成人| 国产成人一区二区三区乱| 中文字幕日韩三级片| 日本黄页网站免费观看| 人妻丰满熟妇AV无码区HD| 国产精品99久久国产小草| 国产青青草自拍视频在线播放| 国产免费人成视频在线观看| 在线观看av网站永久| 粗大的内捧猛烈进出视频|