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

        ?

        基于ADAMS的液壓支架動(dòng)力學(xué)特性仿真與優(yōu)化

        2017-09-11 16:30:13王亞軍
        科技視界 2017年12期
        關(guān)鍵詞:參數(shù)化液壓支架動(dòng)力學(xué)

        王亞軍

        【摘 要】液壓支架對(duì)于煤礦綜采是一種非常重要的支護(hù)設(shè)備。分別采用PRO/E和ADAMS建立液壓支架三維模型并完成動(dòng)力學(xué)分析。液壓支架模型的導(dǎo)入、約束和載荷的定義均基于接口軟件Mech/pro進(jìn)行設(shè)置,該模型的動(dòng)力學(xué)分析采用了參數(shù)化方法,在ADAMS中進(jìn)行仿真與模擬,計(jì)算結(jié)果表明:液壓支架的頂梁在運(yùn)動(dòng)過(guò)程中,橫向與垂向平均加速度的波動(dòng)形狀相似,但數(shù)值差異較大。在優(yōu)化研究方面,采用基于ADAMS的CAE優(yōu)化技術(shù),對(duì)液壓支架的四連桿機(jī)構(gòu)尺寸進(jìn)行優(yōu)化設(shè)計(jì),降低了頂梁端點(diǎn)橫向擺動(dòng)的角位移,從而使得液壓支架的支護(hù)效率得到明顯改善。

        【關(guān)鍵詞】動(dòng)力學(xué);液壓支架;參數(shù)化;優(yōu)化

        0 引言

        隨著煤礦綜采水平的提升,增強(qiáng)液壓支架的研發(fā)技術(shù)、提高液壓支架的設(shè)計(jì)水平、提升液壓支架的可靠性成為產(chǎn)品競(jìng)爭(zhēng)力的基本功能需求[1]。隨著機(jī)電技術(shù)的發(fā)展和進(jìn)步,液壓支架的動(dòng)態(tài)分析與優(yōu)化成為熱點(diǎn)研究問(wèn)題。良好的液壓支架設(shè)計(jì),可使得采煤工作面的移架和推溜速度得到大幅度提升,加強(qiáng)采煤工作面的頂板的支護(hù)強(qiáng)度,為煤礦工人提供更加優(yōu)越的采煤條件,并且能夠減小企業(yè)生產(chǎn)成本。對(duì)于液壓支架的結(jié)構(gòu)優(yōu)化研究,不但可以提升整機(jī)設(shè)備的支護(hù)效率,而且能夠在保證煤炭工人的工作環(huán)境和生產(chǎn)條件改善的同時(shí)降低提高生產(chǎn)效益。

        在煤礦生產(chǎn)中,井下工程技術(shù)人員控制著液壓支架的大部分功能,這部分人員熟練掌控著液壓支架的運(yùn)行和維護(hù),然而卻很少對(duì)液壓支架進(jìn)行深度的研究。一般地,在液壓支架的絕大部分生產(chǎn)廠商中,專(zhuān)業(yè)的從事液壓支架的研發(fā)工程技術(shù)人員很少把精力放在液壓支架的動(dòng)態(tài)特性的分析上,針對(duì)液壓支架的動(dòng)態(tài)特性進(jìn)行改善和優(yōu)化提高的內(nèi)容就更少,缺少人才和技術(shù)支撐,其設(shè)計(jì)研發(fā)內(nèi)容主要集中于結(jié)構(gòu)分析[2-3]。通常,技術(shù)人員提升研發(fā)效率的方法只是在原來(lái)的液壓支架的結(jié)構(gòu)基礎(chǔ)上通過(guò)簡(jiǎn)單的調(diào)整液壓系統(tǒng),將其性能提高和改善,沒(méi)有充分的理論研究和分析計(jì)算支撐,廠家在開(kāi)發(fā)新產(chǎn)品的過(guò)程中只能通過(guò)不斷制作新的實(shí)驗(yàn)樣機(jī),并且需要大量的重復(fù)的調(diào)試試驗(yàn),在樣機(jī)調(diào)試過(guò)程中任何微小問(wèn)題的出現(xiàn),都有極大的可能無(wú)法找到問(wèn)題出現(xiàn)的原因,從而再進(jìn)行反復(fù)的試驗(yàn),直到問(wèn)題解決才能更進(jìn)一步的設(shè)計(jì)和優(yōu)化,此種方法通過(guò)無(wú)數(shù)次不間斷的試驗(yàn),并在此基礎(chǔ)上做出無(wú)數(shù)次調(diào)試和改進(jìn),才能開(kāi)發(fā)和設(shè)計(jì)出新的產(chǎn)品。

        近幾年,CAD和CAE技術(shù)得到快速的發(fā)展和廣泛的應(yīng)用,對(duì)于液壓支架的設(shè)計(jì)與研究,也從經(jīng)驗(yàn)計(jì)算和校核這種最原始的方法過(guò)度到動(dòng)態(tài)數(shù)值模擬分析的方向上。為此,本文運(yùn)用數(shù)值模擬的方法對(duì)液壓支架的動(dòng)力學(xué)特性進(jìn)行分析和與優(yōu)化,提升設(shè)備的支護(hù)效率,增加此設(shè)備帶來(lái)的經(jīng)濟(jì)效益和社會(huì)效益。

        1 液壓支架系統(tǒng)的動(dòng)力學(xué)仿真

        1.1 ADAMS簡(jiǎn)介

        ADAMS是目前公認(rèn)程度最高的動(dòng)力學(xué)分析軟件之一,它是由美國(guó)MDI公司開(kāi)發(fā)出來(lái)的,通過(guò)多次不同版本的更新以及功能的逐步完善,現(xiàn)如今已經(jīng)成為廣泛應(yīng)用的虛擬樣機(jī)分析軟件[4]。ADAMS應(yīng)用到各個(gè)國(guó)家的各個(gè)不同的領(lǐng)域,占據(jù)了全球很大一部分的使用份額,無(wú)論從技術(shù)方面還是從經(jīng)濟(jì)方面而言,都取得了非常大的成功。ADAMS采用的是用戶交互式界面,集結(jié)了多種分析模塊,其主要功能有:(1)可以進(jìn)行系統(tǒng)的動(dòng)力學(xué)以及動(dòng)力學(xué)分析;(2)可以對(duì)液壓系統(tǒng)進(jìn)行模擬仿真;(3)可以對(duì)模型進(jìn)行振動(dòng)模擬的仿真。ADAMS還擁有十分強(qiáng)大的數(shù)據(jù)處理能力,可以根據(jù)用戶的不同需求,得出相應(yīng)的特性曲線,后處理模塊能夠保證分析結(jié)果的清晰度和準(zhǔn)確度。在動(dòng)力學(xué)分析中,多剛體系統(tǒng)動(dòng)力學(xué)理論中的核心為拉格朗日方程方法,可適用于多種復(fù)雜多變的剛性系統(tǒng)。利用ADAMS軟件,用戶可以十分方便地對(duì)虛擬系統(tǒng)進(jìn)行全面的靜力學(xué)、運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)分析,并且基于虛擬樣機(jī)的算法和技術(shù)進(jìn)行模型的二次分析、開(kāi)發(fā)和優(yōu)化。ADAMS具有開(kāi)放性特點(diǎn),具備多程序結(jié)構(gòu)和多種接口,使其成為一種多功能分析平臺(tái),分析結(jié)果可信程度高,特別適用于多自由度系統(tǒng)的動(dòng)態(tài)特性分析。

        1.2 參數(shù)化模型的建立

        在本文的研究中,液壓支架的參數(shù)化建模[5]主要基于三維軟件PRO/E完成。PRO/E是由美國(guó)的PTC公司研發(fā)出的一款不僅能進(jìn)行三維零件設(shè)計(jì),而且還能進(jìn)行三維零件裝配的造型軟件,能夠?qū)崿F(xiàn)裝配體的整體或者局部運(yùn)動(dòng)過(guò)程的數(shù)值模擬,是現(xiàn)階段最具有代表性三維CAD軟件之一。Pro/E在windows平臺(tái)下融合了參數(shù)化造型技術(shù)和變量化造型技術(shù)的三維實(shí)體造型系統(tǒng),其具有十分強(qiáng)大的的機(jī)械裝配設(shè)計(jì)和繪圖功能,能夠方便地與ANSYS、FLUENT、MATLAB等應(yīng)用分析軟件進(jìn)行數(shù)據(jù)間轉(zhuǎn)換和鏈接操作。Pro/E具有及其強(qiáng)大的建模功能,它不僅能完成復(fù)雜的造型設(shè)計(jì)和裝配設(shè)計(jì),將各類(lèi)零件、裝配件等實(shí)體進(jìn)行投影并最終生成二維工程圖,符合我國(guó)繪制標(biāo)準(zhǔn),投影精度很高,在很大的程度上簡(jiǎn)化了液壓支架及其零部件的設(shè)計(jì)過(guò)程。在Pro/E中,組件模塊提供了很多種的約束類(lèi)型和連接形式,約束類(lèi)型包括“匹配”、“對(duì)齊”、“插入”、“坐標(biāo)系”等,這些約束類(lèi)型可以完成相對(duì)靜止的參照設(shè)置;連接形式包括 “銷(xiāo)連接”、 “平面連接”等,這些鏈接類(lèi)型可以完成轉(zhuǎn)動(dòng)副或者平移副的定義。一般情況下,裝配體的約束類(lèi)型需要根據(jù)實(shí)際的自由度來(lái)確定連接形式。針對(duì)液壓支架的機(jī)構(gòu)特點(diǎn),參數(shù)化模型在建立的時(shí)候需要綜合運(yùn)用約束和連接,而且保證機(jī)構(gòu)處于完全約束狀態(tài),即驅(qū)動(dòng)構(gòu)件數(shù)量與自由度數(shù)量一致。

        本文以某型號(hào)的垂直梁式液壓支架為研究對(duì)象進(jìn)行建模分析,該類(lèi)型的液壓支架由插板、導(dǎo)桿、掩護(hù)梁、頂梁、底座、立柱等部件組成。為了將液壓支架的分析過(guò)程變得更加簡(jiǎn)便,本文將組成部件分為兩大類(lèi),一類(lèi)為組焊件,另一類(lèi)為組裝件。 其中,組焊件通過(guò)鋼板或圓管的焊接工藝而連接在一起,形成的連接為剛性連接,因此各部分之間相對(duì)位置是固定不變的,可直接按照裝配關(guān)系或連接關(guān)系進(jìn)行完全裝配;組裝件需要考慮運(yùn)動(dòng)自由度問(wèn)題,根據(jù)旋轉(zhuǎn)副和移動(dòng)副的特點(diǎn)進(jìn)行裝配,實(shí)現(xiàn)不同狀態(tài)下的液壓支架進(jìn)行動(dòng)態(tài)分析。endprint

        Mech/pro是MSC公司開(kāi)發(fā)的,為ADAMS實(shí)現(xiàn)數(shù)據(jù)無(wú)縫接入的連接軟件。文中液壓支架的參數(shù)化模型可以通過(guò)該軟件連接ADAMS與PRO/E,實(shí)現(xiàn)系統(tǒng)動(dòng)力學(xué)模型的參數(shù)化轉(zhuǎn)換,使虛擬樣機(jī)模型的動(dòng)力學(xué)和運(yùn)動(dòng)學(xué)計(jì)算更加高效。Mech/pro可以在 Pro/E中直接定義模型的剛體特性和相關(guān)的自由度約束,并將其傳遞到ADAMS中,實(shí)現(xiàn)全面的動(dòng)力學(xué)聯(lián)合仿真分析。

        在本文中,首先通過(guò)Pro/E三維軟件建立液壓支架的裝配圖模型,再通過(guò)Mech/pro接口軟件,完成虛擬樣機(jī)參數(shù)化轉(zhuǎn)化并導(dǎo)入到ADAMS中。在液壓支架的動(dòng)力學(xué)分析中,需要根據(jù)實(shí)際工作條件設(shè)置不同的約束類(lèi)型和連接形式,從而實(shí)現(xiàn)液壓支架的虛擬裝配和運(yùn)動(dòng)仿真。為了對(duì)液壓支架進(jìn)行更加精確地運(yùn)動(dòng)仿真,文中對(duì)四個(gè)立柱和各個(gè)油缸兩端的銷(xiāo)孔與其相關(guān)部件上耳座之間采用“銷(xiāo)連接”;活塞桿與缸體之間、插板與掩護(hù)梁之間、導(dǎo)桿與導(dǎo)向筒、前伸梁與頂梁之間采用“滑動(dòng)連接”。通過(guò)Mech/pro軟件,不僅仿真模型的相似度得到提升,工作效率也同樣得到提高。當(dāng)在接口軟件中施加約束時(shí),可以直接參選零件的各個(gè)特征,在導(dǎo)出的數(shù)據(jù)中包含有接口軟件中的所有數(shù)據(jù),真正地實(shí)現(xiàn)了無(wú)縫參數(shù)化轉(zhuǎn)換。為了操作方便,能在 Mech/pro 中添加約束的盡量在 Mech/pro 中添加。動(dòng)力學(xué)分析軟件ADAMS能夠直接識(shí)別的文件為cmd,通過(guò)該文件的導(dǎo)入,可得到參數(shù)化的液壓支架虛擬樣機(jī)模型,最終各個(gè)約束在ADAMS中如圖1所示。

        1.3 結(jié)果分析

        為了更好的模擬仿真液壓支架的運(yùn)動(dòng)學(xué)特性,需要完成邊界條件的定義。根據(jù)工況,文中定義頂梁驅(qū)動(dòng)電機(jī)的垂向速度為6mm/s,并設(shè)定支架從3000mm處開(kāi)始運(yùn)行。根據(jù)支護(hù)分析結(jié)果,對(duì)頂梁處和底座處添加載荷。在液壓支架虛擬樣機(jī)的載荷條件中,底座兩端施加的支撐力是由底座實(shí)際工況中的支撐力簡(jiǎn)化而來(lái)的,因此,為了更接近實(shí)際工況條件,文中將頂梁的壓力載荷按照結(jié)構(gòu)進(jìn)行分散處理。在ADAMS中,設(shè)置液壓支架升降運(yùn)動(dòng)為勻速運(yùn)動(dòng),并且建立支架動(dòng)載荷隨時(shí)間變化的關(guān)系,從而可以計(jì)算出液壓支架在不同高度下所對(duì)應(yīng)的動(dòng)力學(xué)特性,研究液壓支架的承載能力與支架高度變化之間的關(guān)系。在完成系統(tǒng)的設(shè)定之后,即可以進(jìn)行動(dòng)力學(xué)仿真分析,最終得到頂梁在空間(橫向與垂向)運(yùn)動(dòng)的曲線圖分別如圖2所示。

        根據(jù)圖2中的分析結(jié)果可知:液壓支架頂梁在橫向的運(yùn)動(dòng)較為穩(wěn)定,加速度變化較為平緩,液壓支架在移動(dòng)支護(hù)過(guò)程中,水平方向的慣性附加力比較小,慣性頻率較低;在豎直方向上,液壓支架的運(yùn)動(dòng)特性和橫向方向相似,但平均加速度在數(shù)值上為橫向方向的10倍以上,其慣性力相對(duì)于水平方向而言較大。為了適應(yīng)液壓支架的支護(hù)方向,需要進(jìn)一步對(duì)液壓支架的四連桿結(jié)構(gòu)進(jìn)行優(yōu)化分析和設(shè)計(jì)。

        2 液壓支架連桿系統(tǒng)優(yōu)化

        2.1 優(yōu)化原理

        本文對(duì)液壓支架連機(jī)構(gòu)的優(yōu)化方法主要基于ADAMS。在ADAMS中,采用優(yōu)化變量是應(yīng)用最多的參數(shù)化分析方法之一,它不但可以定義用戶本身的參數(shù)變量,同時(shí)還可以關(guān)聯(lián)仿真對(duì)象進(jìn)行結(jié)果分析[6]。在ADAMS/View中,為設(shè)計(jì)變量的選擇提供了一種可靠的設(shè)計(jì)方法,可以清晰的顯示各個(gè)參數(shù),這樣就能很方便地觀察、修改各個(gè)參數(shù)變量值。在給定初始參數(shù)后,通過(guò)參數(shù)化的選擇,可以將設(shè)計(jì)變量在一定的變化范圍內(nèi)連續(xù)取值,進(jìn)行相應(yīng)的分析,從而完成參變量的優(yōu)化。ADAMS主要提供了設(shè)計(jì)研究、試驗(yàn)設(shè)計(jì)和優(yōu)化分析等三種參變量的研究方法。

        通常,對(duì)于需要進(jìn)行優(yōu)化的變量在整個(gè)分析研究中是與其它變量相互關(guān)聯(lián)的,不同的變量之間存在著不同的函數(shù)關(guān)系。在優(yōu)化結(jié)果的計(jì)算中,必須使優(yōu)化變量保持對(duì)目標(biāo)函數(shù)的一致性,為此必須轉(zhuǎn)換參數(shù)變量,若不進(jìn)行參數(shù)變量的轉(zhuǎn)換,得出的優(yōu)化結(jié)果不僅使用戶無(wú)法明確的選擇參變量,甚至得到的優(yōu)化結(jié)果也不理想,無(wú)法有效地分析液壓支架的結(jié)構(gòu)性能。本文將液壓支架的四桿機(jī)構(gòu)作為優(yōu)化結(jié)構(gòu),優(yōu)化目標(biāo)多個(gè),因此計(jì)算結(jié)果也會(huì)出現(xiàn)不同的數(shù)值,每個(gè)優(yōu)化計(jì)算得出的數(shù)值都是獨(dú)立的,無(wú)法得出它們之間的關(guān)系,基于此種結(jié)果,必須通過(guò)某種方法使計(jì)算結(jié)果的數(shù)值唯一?,F(xiàn)今,在對(duì)模型進(jìn)行優(yōu)化時(shí),有三種主要方法可使用戶在計(jì)算中得到最終解:(1)生成法。這種方法是計(jì)算出無(wú)數(shù)的對(duì)分析優(yōu)化結(jié)果有用的數(shù)值,將這些數(shù)值構(gòu)成求解子集,然后通過(guò)篩選校正這些子集得到最優(yōu)解;(2)交互法。交互法的應(yīng)用對(duì)研究人員自身具備的經(jīng)驗(yàn)和學(xué)術(shù)水平要求較高,是研究人員直接對(duì)得到的多個(gè)解進(jìn)行有效地判斷和分析,從而得到最優(yōu)解;(3)轉(zhuǎn)化法。把多個(gè)優(yōu)化目標(biāo)變?yōu)橐粋€(gè)優(yōu)化目標(biāo)進(jìn)行優(yōu)化是轉(zhuǎn)換法的基本原理。在優(yōu)化計(jì)算最終解的過(guò)程中,通常都是運(yùn)用特定的優(yōu)化算法去計(jì)算結(jié)果,實(shí)現(xiàn)對(duì)目標(biāo)函數(shù)值的優(yōu)化,比如,多目標(biāo)進(jìn)化算法、多目標(biāo)粒子群算法和蟻群算法、模擬退火算法等。

        2.2 優(yōu)化模型的建立

        本文針對(duì)液壓支架四連桿結(jié)構(gòu),建立力學(xué)分析模型(結(jié)構(gòu)參數(shù)如圖3所示),橫向擺動(dòng)值最小的頂梁端點(diǎn)是遺傳算法參數(shù)化的主要優(yōu)化目標(biāo),即ω最小。

        把優(yōu)化后的參數(shù)在ADAMS中進(jìn)行仿真計(jì)算,可驗(yàn)證所得參數(shù)是否符合要求。文中所選用的遺傳方法即達(dá)到了結(jié)果的準(zhǔn)確又保證了計(jì)算的快速性。針對(duì)優(yōu)化目標(biāo),文中對(duì)四連桿中的l1、l2、l3、l4和 l5進(jìn)行優(yōu)化。

        2.3 結(jié)果分析

        在ADAMS/View中,設(shè)計(jì)變量既可以用于存儲(chǔ)數(shù)據(jù)又可以作為參數(shù)設(shè)計(jì)過(guò)程中的變量[7]。實(shí)現(xiàn)參數(shù)變量主要采用的方法是Create Design Variable,變量類(lèi)型及及取值的選擇是利用變量生成對(duì)話框來(lái)進(jìn)行的,本文中,設(shè)置設(shè)計(jì)變量分別為l1、l2、l3、l4以及l(fā)5長(zhǎng)度。

        在ADAMS/View中,通過(guò)多次連續(xù)迭代計(jì)算,得出三種優(yōu)化組合方案。通過(guò)反復(fù)的對(duì)比分析,文中采用效果更佳的優(yōu)化方案1,如表1所示,優(yōu)化后的頂梁端點(diǎn)軌跡橫向擺動(dòng)值明顯減小,這對(duì)于提升液壓支架的工作效率有著重要的意義。

        3 結(jié)語(yǔ)

        本文通過(guò)Pro/E三維軟件建立了液壓支架的三維模型,并通過(guò)Mech/pro導(dǎo)入到了ADAMS當(dāng)中,建立了剛性參數(shù)化的模型,實(shí)現(xiàn)動(dòng)態(tài)特性的數(shù)值模擬,結(jié)果表明,液壓支架頂梁在橫向和垂向的加速度在數(shù)值上有著非常大的差異,其慣性力相對(duì)于水平方向而言較大?;贏DAMS對(duì)液壓支架的主體四連桿結(jié)構(gòu)尺寸進(jìn)行了優(yōu)化,通過(guò)連續(xù)迭代計(jì)算,得出三種優(yōu)化組合方案,最終采用效果更佳的優(yōu)化方案提升了液壓支架的頂梁前端點(diǎn)的軌跡效率,使得橫向擺動(dòng)的角位移降至1.2°,得到良好的優(yōu)化效果。

        【參考文獻(xiàn)】

        [1]王國(guó)法.液壓支架技術(shù)體系研究與實(shí)踐[J].煤炭學(xué)報(bào),2010,35(11):1903-1908.

        [2]趙宏珠,宋秋爽.特大采高液壓支架發(fā)展與研究[J].采礦與安全工程學(xué)報(bào), 2007,24(3):265-269.

        [3]王國(guó)法.工作面支護(hù)與液壓支架技術(shù)理論體系[J].煤炭學(xué)報(bào),2014,39(8):1593-1601.

        [4]王立新,李震,王鎖吉.基于ADAMS的液壓支架運(yùn)動(dòng)學(xué)分析[J].煤礦機(jī)械, 2008,29(12):82-84.

        [5]蔡余龍,秦東晨,武紅霞,等.基于Pro/E和ADAMS的液壓支架虛擬樣機(jī)建模與仿真研究[J].煤礦機(jī)械,2009,30(11):36-38.

        [6]章強(qiáng),涂少伍,朱維柱,等.基于ADAMS的兩柱放頂煤液壓支架運(yùn)動(dòng)仿真分析[J].煤炭技術(shù),2015(11):301-302.

        [7]蔡文書(shū).基于SolidWorks和ADAMS的液壓支架動(dòng)力學(xué)仿真[J].中州煤炭, 2016(4):99-102.

        [責(zé)任編輯:朱麗娜]endprint

        猜你喜歡
        參數(shù)化液壓支架動(dòng)力學(xué)
        小天體環(huán)的軌道動(dòng)力學(xué)
        Pro/E的三維往復(fù)壓縮機(jī)參數(shù)化模型庫(kù)的建立
        淺析煤礦用液壓支架的輕量化設(shè)計(jì)
        一種懸架運(yùn)動(dòng)仿真快速建模方法研究
        針對(duì)優(yōu)化液壓支架立柱大修工作的對(duì)策探析
        綜采工作面電液控制支架在翟鎮(zhèn)煤礦的應(yīng)用
        基于Pro/ENGINEER的圓柱齒輪參數(shù)化虛擬設(shè)計(jì)
        綜采工作面設(shè)備井下倒安裝工藝的實(shí)踐應(yīng)用
        基于隨機(jī)-動(dòng)力學(xué)模型的非均勻推移質(zhì)擴(kuò)散
        基于MBD的飛機(jī)結(jié)構(gòu)件重量計(jì)算
        激情人妻中出中文字幕一区| 网站在线观看视频一区二区| 日本一区二区在线免费视频| 韩国三级中文字幕hd| 推油少妇久久99久久99久久| 国产高清天干天天视频| 国产精品人人爱一区二区白浆| 日本a爱视频二区三区| 国产福利永久在线视频无毒不卡| 老少交欧美另类| 亚洲va成无码人在线观看| 亚洲精品中文字幕不卡| 久久国产成人精品国产成人亚洲| 国产中文字幕乱人伦在线观看| 99精品国产高清一区二区麻豆| 麻豆av传媒蜜桃天美传媒| 国产自产拍精品视频免费看| 亚洲精品第四页中文字幕| 久久青青草原国产毛片| 女人下面毛多水多视频| 日韩亚洲制服丝袜中文字幕| 久久精品女人天堂av麻| 国产精品理论片在线观看| 女同性黄网aaaaa片| 亚洲av午夜成人片精品| 饥渴少妇一区二区三区| 人妻少妇猛烈井进入中文字幕 | av天堂手机在线免费| av中文字幕一区不卡| 久久国产精品-国产精品| 日产精品一区二区三区| 亚洲国产精品午夜一区| 麻豆精品国产av在线网址| 国产欧美日韩久久久久| 国产成人亚洲精品2020| 久久这里都是精品99| 久久性爱视频| 国产精品99久久免费| 国产高跟丝袜在线诱惑| 日本一区二区三区视频网站| 丰满人妻av无码一区二区三区|