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

        ?

        空間微振動(dòng)模擬器柔性鉸鏈設(shè)計(jì)

        2014-12-19 01:46:46徐振邦吳清文陳立恒
        制造業(yè)自動(dòng)化 2014年18期
        關(guān)鍵詞:鉸鏈軸向柔性

        石 柱,徐振邦,吳清文,陳立恒

        SHI Zhu, XU Zhen-bang, WU Qing-wen, CHEN Li-heng

        (中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,長(zhǎng)春 130033)

        0 引言

        空間高分辨率光學(xué)遙感器有重要的商業(yè)價(jià)值和軍事價(jià)值,且迅速發(fā)展[1,2]??臻g設(shè)施的微振動(dòng)對(duì)其分辨率產(chǎn)生較大影響,對(duì)其振動(dòng)特性進(jìn)行研究從而實(shí)現(xiàn)有效隔振有重大意義。而空間微振動(dòng)有微振幅、多自由度、隨機(jī)、寬頻的運(yùn)動(dòng)特點(diǎn)。

        Stewart平臺(tái)作為一種并聯(lián)機(jī)構(gòu),由上平臺(tái)、下平臺(tái)和六根伸縮桿組成,下平臺(tái)固定,通過(guò)控制六根支腿的伸縮實(shí)現(xiàn)上平臺(tái)的六自由度運(yùn)動(dòng)[3],可以有效模擬空間微振動(dòng)。Stewart平臺(tái)上下平臺(tái)與伸縮桿關(guān)節(jié)處所用鉸鏈的特性對(duì)于整個(gè)機(jī)構(gòu)的動(dòng)態(tài)特性以及控制特性有著至關(guān)重要的影響。傳統(tǒng)剛性鉸鏈,會(huì)產(chǎn)生摩擦和間隙,而柔性鉸鏈具有無(wú)回退空程、無(wú)摩擦、高精度等優(yōu)點(diǎn),在實(shí)現(xiàn)精密定位和隔振平臺(tái)也有著眾多的應(yīng)用[4]。柔性鉸鏈可分為單向、雙向和萬(wàn)向??臻g微振動(dòng)模擬器要實(shí)現(xiàn)六自由度的運(yùn)動(dòng),根據(jù)并聯(lián)機(jī)器人機(jī)構(gòu)末端的運(yùn)動(dòng)特征是構(gòu)成該機(jī)構(gòu)的所有運(yùn)動(dòng)支鏈運(yùn)動(dòng)特征的交集理論[5],單腿必須實(shí)現(xiàn)六自由度的運(yùn)動(dòng)特征,因此選用萬(wàn)向轉(zhuǎn)動(dòng)柔性鉸鏈鏈。

        萬(wàn)向柔性鉸鏈在Stewart平臺(tái)上作為隔振器和精確指向機(jī)構(gòu)的應(yīng)用,都僅考慮到鉸鏈的轉(zhuǎn)向能力(旋轉(zhuǎn)剛度),而對(duì)軸向剛度的影響忽略[6]。對(duì)于微振動(dòng)模擬器在高頻振動(dòng)時(shí),不僅要求鉸鏈有較低的彎曲剛度和扭轉(zhuǎn)剛度,能有效實(shí)現(xiàn)關(guān)節(jié)的轉(zhuǎn)動(dòng),也必須有高的軸向剛度,提高關(guān)節(jié)耦合處剛度,增強(qiáng)系統(tǒng)基頻,使系統(tǒng)具有良好的動(dòng)態(tài)特性和控制特性。而對(duì)于萬(wàn)向鉸鏈,其軸向剛度、彎曲剛度及扭轉(zhuǎn)剛度特性變化趨勢(shì)一致(軸向剛度越大,彎曲剛度和扭轉(zhuǎn)剛度也就越大),這和鉸鏈設(shè)計(jì)要求相沖突。本文對(duì)直圓型萬(wàn)向柔性鉸鏈剛度計(jì)算進(jìn)行了深入研究,通過(guò)材料力學(xué)的公式計(jì)算分析鉸鏈結(jié)構(gòu)參數(shù)對(duì)各向剛度的影響,并輔助于有限元軟件仿真進(jìn)行驗(yàn)證,綜合考慮鉸鏈軸向剛度和旋轉(zhuǎn)剛度要求,選取合理結(jié)構(gòu)參數(shù)使得鉸鏈滿足指標(biāo)要求,有較高的軸向剛度和較低的旋轉(zhuǎn)剛度,對(duì)于萬(wàn)向柔性鉸鏈的剛度設(shè)計(jì)有一定的指導(dǎo)意義。

        1 設(shè)計(jì)思路與指標(biāo)

        依據(jù)仿真利用模型依據(jù)機(jī)構(gòu)的運(yùn)動(dòng)指標(biāo)要求提出柔性鉸鏈的各向剛度要求,然后建立模型,用公式推導(dǎo)出來(lái)剛度與鉸鏈結(jié)構(gòu)參數(shù)的關(guān)系公式,得到剛度特性與結(jié)構(gòu)參數(shù)的規(guī)律,利用有限元進(jìn)行仿真計(jì)算驗(yàn)證。根據(jù)運(yùn)動(dòng)指標(biāo)和材料強(qiáng)度參數(shù)確定柔性鉸鏈的結(jié)構(gòu)尺寸和剛度范圍,確定了最佳結(jié)構(gòu)參數(shù),然后檢驗(yàn)鉸鏈取此結(jié)構(gòu)是否滿足剛度指標(biāo)要求,之后對(duì)其強(qiáng)度進(jìn)行了校核,完成鉸鏈部分的結(jié)構(gòu)設(shè)計(jì)。

        微振動(dòng)模擬器結(jié)構(gòu)如圖1所示,下平臺(tái)固定,每個(gè)支腿都由一個(gè)音圈電機(jī)作為作動(dòng)器,支腿和上下平臺(tái)通過(guò)萬(wàn)向柔性鉸鏈連接。模擬器上平臺(tái)提供的振動(dòng)指標(biāo)要求為:六自由度、頻率30Hz~400Hz,加速度保持在2×10-3g。

        圖1 微振動(dòng)模擬器結(jié)構(gòu)圖

        當(dāng)模擬器高頻振動(dòng)時(shí),如果鉸鏈的軸向剛度低,使得支腿與平臺(tái)耦合處容易發(fā)生共振,影響控制精度,甚至導(dǎo)致控制失穩(wěn)。以極限高頻400Hz振動(dòng)為考慮目標(biāo),使關(guān)節(jié)耦合處不會(huì)發(fā)生共振,根據(jù)MATLAB理論模型換算出鉸鏈的軸向剛度應(yīng)高于3.4×104N/mm。

        對(duì)于旋轉(zhuǎn)剛度,較高則會(huì)產(chǎn)生力的干擾,影響動(dòng)態(tài)特性和控制特性。若模型使用理想鉸鏈,在平均工況(振動(dòng)頻率200Hz,電機(jī)力30N),鉸鏈的彎曲角度和扭轉(zhuǎn)角度分別為0.83’和0.037’。要求柔性鉸鏈能達(dá)到理想鉸鏈90%的轉(zhuǎn)向能力,基于此,根據(jù)Matlab理論模型換算出鉸鏈的彎曲剛度應(yīng)低于2.45×105N·mm/rad,扭轉(zhuǎn)剛度低于4.3×104N·mm/rad。

        2 柔性鉸鏈的剛度分析

        直圓形萬(wàn)向柔性鉸鏈,其結(jié)構(gòu)如圖2所示,在半徑R的圓柱體的中間部分切半徑r的槽,切槽部分厚度為t,三者關(guān)系為R=2r+t。由于設(shè)計(jì)的Stewart平臺(tái)的結(jié)構(gòu)尺寸限制,h取30mm,因此考慮R和t對(duì)鉸鏈剛度的影響即可。鉸鏈選取的材料為鈹銅合金,彈性模量E=128GPa,泊松比ν=0.35,強(qiáng)度極限σS=1035Mpa,安全系數(shù)取n=2,[σ]=517.5Mpa,[τ]=362.25Mpa。

        圖2 柔性鉸鏈結(jié)構(gòu)圖

        2.1 軸向剛度

        在鉸鏈上端施加垂直載荷P,其軸向剛度為K,軸向變形為l,l1為P引起的鉸鏈未切槽部分的軸向變形,l2為P引起的切槽部分(柔性部分)的軸向變形。根據(jù)材料力學(xué)公式:

        其中A表示未切槽部分橫截面積,ε表示正應(yīng)變,σ表示正應(yīng)力。

        利用MATLAB編寫程序?qū)ζ溥M(jìn)行計(jì)算[7]。綜合考慮R與t的影響,分別取R=5,6,7;t=1,1.5,2,2,5,3,3.5,4,4.5,5,5.5,6,8,10,單位為mm,作出R取不同值時(shí),柔性鉸鏈切槽部分和未切槽部分的軸向變形隨切槽部分厚度t的變化曲線如圖3所示。

        圖3 軸向變形隨切槽部分厚度t變化曲線

        通過(guò)我們可以得出:

        1)在h一定的情況下,柔性鉸鏈的軸向變形,幾乎不受半徑R的影響,由切槽部分的厚度t決定。

        2)在t大于5mm的時(shí)候,切槽部分的軸向變形和未切槽部分幾乎無(wú)差別,而在t小于5mm的時(shí)候,切槽部分的軸向變形大于未切槽部分。

        3)未切槽部分軸向變形幾乎不受t影響,而切槽部分軸向變形隨著t的變小呈指數(shù)型增長(zhǎng)。

        基于上述結(jié)論,取R=5mm。將切槽部分軸向變形和未切槽部分軸向變形相加得出其軸向變形l,根據(jù)公式(1)可以得出其軸向剛度K,作出其隨切槽部分厚度t變化的曲線。為了驗(yàn)證理論推導(dǎo)的正確性,同時(shí)通過(guò)有限元分析軟件Patran/Nastran對(duì)柔性鉸鏈進(jìn)行仿真。材料參數(shù)設(shè)置為E=128000MPa,ν=0.35,h=30mm,R=5mm,t分別1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,8,10(單位為mm),建立單元數(shù)為29820的有限元模型,在上表面施加total load 60N,對(duì)其進(jìn)行計(jì)算機(jī)仿真。然后提取軸向變形,并計(jì)算得到其軸向剛度,將其結(jié)果與理論結(jié)果相比較,得出結(jié)果如圖4所示。

        圖4 軸向剛度隨切槽厚度t變化曲線

        通過(guò)仿真結(jié)果與理論計(jì)算的結(jié)果的比較,兩條曲線基本一致,可以驗(yàn)證我們對(duì)于柔性鉸鏈軸向剛度受結(jié)構(gòu)參數(shù)變化的影響所作理論推導(dǎo)的正確性。分析曲線我們可以得出:軸向剛度隨著切槽厚度t的增加成近似線性提高。

        2.2 彎曲剛度

        假定在鉸鏈上端施加彎矩M,其彎曲剛度為Kθ,轉(zhuǎn)角為θ,C為倒數(shù)彎曲剛度(定義為M引起鉸鏈的轉(zhuǎn)角與M之比,C越大)。根據(jù)材料力學(xué)公式:

        上式中J(X)為鉸鏈各部分截面慣性矩,C1為未切槽部分倒數(shù)彎曲剛度,C2為切槽(柔性)部分倒數(shù)彎曲剛度,其值越大,表示這一部分的彎曲剛度越小。利用MATLAB編寫程序?qū)ζ溥M(jìn)行計(jì)算,分別取R=5,6,7;t=1,1.5,2,2,5,3,3.5,4,4.5,5,5.5,6,8,10(單位為mm),作出R取不同值時(shí),柔性鉸鏈切槽部分和未切槽部分的倒數(shù)彎曲剛度隨切槽部分厚度t的變化曲線,如圖5所示。

        圖5 倒數(shù)彎曲剛度隨切槽部分厚度t變化曲線

        通過(guò)上圖我們可以得出:

        在h一定的情況下,柔性鉸鏈的倒數(shù)彎曲剛度,幾乎不受半徑R的影響,由切槽部分的厚度t決定。

        在t大于4mm的時(shí)候,切槽部分的倒數(shù)彎曲剛度和未切槽部分都非常小,也就是鉸鏈在M的作用下幾乎沒(méi)有彎曲。而在t小于5mm的時(shí)候,切槽部分的倒數(shù)彎曲剛度相對(duì)較大。

        未切槽部分的倒數(shù)彎曲剛度幾乎為0,切槽部分倒數(shù)彎曲剛度隨著t的變小呈指數(shù)型增長(zhǎng)。

        基于上述結(jié)論,我們?nèi)=5mm,并將切槽部分軸向變形和未切槽部分倒數(shù)彎曲剛度相加得出C,由,計(jì)算出彎曲剛度的值,作出其隨t變化曲線。并利用Patran/Nastran有限元對(duì)其進(jìn)行仿真(t=2時(shí)有限元模型如圖7所示),給柔性鉸鏈一端施加彎矩M,提取出斷面的彎曲變形,再根據(jù)公式(4)求出其彎曲剛度,將仿真值與理論值相比較,結(jié)果如圖6所示。

        圖6 彎曲剛度隨切槽厚度t變化曲線

        通過(guò)仿真結(jié)果與理論計(jì)算的結(jié)果的比較,可以再一次驗(yàn)證對(duì)于柔性鉸鏈彎曲剛度因結(jié)構(gòu)參數(shù)變化影響理論分析的正確性。觀察彎曲剛度隨切槽厚度t變化曲線可以得知:彎曲剛度的大小隨著t的變大而幾乎成指數(shù)型增長(zhǎng)。

        2.3 扭轉(zhuǎn)剛度

        當(dāng)在鉸鏈一端施加扭矩T時(shí),鉸鏈的扭轉(zhuǎn)角為φ,根據(jù)材料力學(xué)公式:

        根據(jù)剪切模量G和彈性模量E的關(guān)系:G=E/2(1+ν),聯(lián)立公式(9)和公式(6),可得e=C(1+ν),最終得出扭轉(zhuǎn)剛度與彎曲剛度的關(guān)系式:

        上式中,IP為截面極慣性矩,G為剪切模量,Kφ為扭轉(zhuǎn)剛度,Kθ為彎曲剛度,ν為泊松比。由上公式我們可以看出扭轉(zhuǎn)剛度和彎曲剛度之間是一個(gè)線性關(guān)系。材料鈹銅合金ν=0.13,Kθ=Kθ/1.13。所以它隨結(jié)構(gòu)參數(shù)影響的變化趨勢(shì)與彎曲剛度一致。

        同樣取R=5,計(jì)算Kφ理論值,分析其隨t變化曲線,并同時(shí)利用Patran/Nastran有限元對(duì)其進(jìn)行仿真,并且將理論結(jié)果與有限元結(jié)果進(jìn)行比較,如圖7所示。

        圖7 扭轉(zhuǎn)剛度隨切槽厚度t變化曲線

        理論值與仿真值的比較基本吻合,進(jìn)一步說(shuō)明理論推導(dǎo)公式結(jié)果的正確性。扭轉(zhuǎn)剛度隨切槽厚度t的變化趨勢(shì)與彎曲剛度完全一致。

        在得出柔性鉸鏈的各向剛度隨切槽厚度t的變化規(guī)律之后,再結(jié)合我們對(duì)柔性鉸鏈的設(shè)計(jì)要求(較高的軸向剛度和較低的彎曲剛度與扭轉(zhuǎn)剛度),可以得到下列思路:切槽部分厚度t越大,軸向剛度越大,旋轉(zhuǎn)剛度也越大,這于結(jié)構(gòu)參數(shù)的選擇是一個(gè)矛盾。因此,要在軸向剛度和旋轉(zhuǎn)剛度的設(shè)計(jì)中選取一個(gè)平衡:即選取合適的參數(shù),使得旋轉(zhuǎn)剛度滿足要求的前提下,軸向剛度也盡量的大。

        3 參數(shù)確定

        在得出柔性鉸鏈的各向剛度隨切槽厚度t的變化規(guī)律之后,根據(jù)材料的許用應(yīng)力和鉸鏈的剛度指標(biāo)要求確定結(jié)構(gòu)的參數(shù)。

        Stewart平臺(tái)單腿電機(jī)驅(qū)動(dòng)力最大取P=60N,單腿的最大彎曲角度和扭轉(zhuǎn)角度為θ=±1°,φ=±1°,鉸鏈材料取彈性強(qiáng),強(qiáng)度高的鈹銅合金,E=128GPa,v=0.35,σs=1035Mpa,安全系數(shù)n取2,[σ]=517.5Mpa,[τ]=362.25Mpa。

        3.1 鉸鏈的正應(yīng)力

        柔性鉸鏈的正應(yīng)力由其軸向的電機(jī)驅(qū)動(dòng)力P和彎矩M產(chǎn)生??紤]運(yùn)動(dòng)和動(dòng)力指標(biāo)都取極限值的情況:P和M都取最大值,對(duì)其進(jìn)行正應(yīng)力強(qiáng)度校驗(yàn),根據(jù)材料力學(xué)的公式:

        分別將t取不同值時(shí)對(duì)應(yīng)的軸向剛度和彎曲剛度代入計(jì)算,由公式(4)計(jì)算出柔性鉸鏈的最大彎矩,代入公式(11),得到在運(yùn)動(dòng)指標(biāo)達(dá)到極限的時(shí)候,鉸鏈最大正應(yīng)力和切槽厚度t的變化曲線,如圖8所示。

        圖8 最大正應(yīng)力隨切槽厚度變化曲線

        3.2 鉸鏈的剪切應(yīng)力

        鉸鏈的剪切應(yīng)力由轉(zhuǎn)矩和扭矩產(chǎn)生,在運(yùn)動(dòng)指標(biāo)取極限值的情況下,將t取不同值時(shí),鉸鏈的彎曲剛度和扭轉(zhuǎn)剛度分代入下列公式計(jì)算,得到最大剪切應(yīng)力和切槽厚度t的變化曲線,如圖9所示。

        上式中τ 1、τ 2分別為彎矩和扭矩產(chǎn)生的剪切應(yīng)力。

        圖9 最大剪切應(yīng)力隨切槽厚度變化曲線

        由正應(yīng)力和剪切應(yīng)力的校驗(yàn)可以得知,要滿足柔性鉸鏈的強(qiáng)度要求,切槽厚度t,只能取1mm,1.5mm,2mm,8mm,10mm。結(jié)合第二節(jié)提出的指標(biāo):軸向剛度高于3.4×104N/mm彎曲剛度低于2.45×105N·mm/rad,扭轉(zhuǎn)剛度低于4.3×104N·mm/rad,t=2mm,為最終柔性鉸鏈參數(shù)選擇。其軸向剛度為8.50×104N/mm,彎曲剛度為3.72×104Nm/rad,扭轉(zhuǎn)剛度為2.75×104N·m/rad,均滿足指標(biāo)要求。

        因此,對(duì)于柔性鉸鏈的參數(shù)最終確定:h=30mm,R=5mm,t=2mm。

        4 有限元校核

        建立有限元模型后,給其同時(shí)施加極限動(dòng)力指標(biāo)和運(yùn)動(dòng)指標(biāo)情況下的軸向軸向載荷PMAX,MMAX,TMAX,然后查看應(yīng)力云圖如圖10所示。可以得知,最大應(yīng)力發(fā)生在鉸鏈切槽部分,其值約為300MPa,小于材料的許用應(yīng)力,即強(qiáng)度校核滿足要求,設(shè)計(jì)合理。

        圖10 強(qiáng)度校核結(jié)果

        5 結(jié)論

        本文基于Stewart平臺(tái)作空間微振動(dòng)模擬器的設(shè)計(jì)要求,選用萬(wàn)向柔性鉸鏈作為機(jī)構(gòu)關(guān)節(jié)。研究了萬(wàn)向柔性鉸鏈各向剛度特性與其結(jié)構(gòu)參數(shù)的關(guān)系,并結(jié)合整個(gè)Stewart平臺(tái)結(jié)構(gòu)運(yùn)動(dòng)學(xué)指標(biāo)、動(dòng)力學(xué)指標(biāo)和強(qiáng)度特性對(duì)柔性鉸鏈進(jìn)行結(jié)構(gòu)參數(shù)設(shè)計(jì),然后用有限元進(jìn)行校核,完成了用作空間微振動(dòng)模擬器的Stewart平臺(tái)的柔性鉸鏈設(shè)計(jì),且滿足指標(biāo)要求。

        萬(wàn)向柔性鉸鏈的剛度特性主要由切槽部分的厚度t決定,其軸向剛度和旋轉(zhuǎn)剛度都隨著切槽厚度t的增大而變大。對(duì)于Stewart平臺(tái)的關(guān)節(jié)處,要求鉸鏈的軸向剛度較高而旋轉(zhuǎn)剛度較低。因此對(duì)其關(guān)鍵參數(shù)t的確定,需要結(jié)合具體運(yùn)動(dòng)指標(biāo)、動(dòng)力指標(biāo)和材料強(qiáng)度特性進(jìn)行分析。一般在滿足材料強(qiáng)度特性的情況下,選取的切槽厚度t要使得旋轉(zhuǎn)剛度盡量小而軸向剛度盡量大。

        [1]ATAD E,et al..Atlast-9.2m:a large-aperture deployable space telescope[J].Proc.of SPIE,2010,7731,77312M1-10.

        [2]葉偉楠,董吉洪.大口徑主鏡輕量化結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計(jì)[J].中國(guó)光學(xué),2012,5(3):222-228.

        [3]Stewart D.A platform with six degrees of freedom[J].Proc.Inst.Mech.Eng,1965,180(15):371-378.

        [4]張定會(huì).柔性鉸鏈實(shí)現(xiàn)微位移的方法研究[J].工業(yè)儀表與自動(dòng)化裝置,1999,(5):11-13.

        [5]Byung JY,Heung YN,Goo B C etc.Design and Experiment of a3DOF Parallel Micro-Mechanism Utilizing Flexure Hinges[J].IEEE International conference on Robotics and Automation,2002.

        [6]高鵬,袁哲俊,姚英學(xué).基于柔性鉸鏈結(jié)構(gòu)的新型雙向微動(dòng)工作臺(tái)的研究[J].儀器儀表學(xué)報(bào),1998,19(2):192-195 .

        [7]劉偉,單軸柔性鉸鏈轉(zhuǎn)角剛度的計(jì)算機(jī)輔助計(jì)算仁[J].光學(xué)精密工程,2000,8(2):178-180.

        猜你喜歡
        鉸鏈軸向柔性
        一種柔性拋光打磨頭設(shè)計(jì)
        大型立式單級(jí)引黃離心泵軸向力平衡的研究
        灌注式半柔性路面研究進(jìn)展(1)——半柔性混合料組成設(shè)計(jì)
        石油瀝青(2021年5期)2021-12-02 03:21:18
        高校學(xué)生管理工作中柔性管理模式應(yīng)用探索
        基于虛擬鉸鏈打開(kāi)機(jī)構(gòu)的艙門提升機(jī)構(gòu)研究
        球鉸鏈防塵罩抱緊力優(yōu)化
        北京汽車(2019年6期)2019-02-12 05:19:50
        荒銑加工軸向切深識(shí)別方法
        汽車連接器帶鉸鏈護(hù)殼產(chǎn)品的塑料模具設(shè)計(jì)改進(jìn)
        汽車電器(2018年10期)2018-11-01 01:17:26
        微小型薄底零件的軸向車銑實(shí)驗(yàn)研究
        更換筆記本電腦鉸鏈
        CHIP新電腦(2014年5期)2014-05-14 20:54:03
        无码国产精品一区二区免费网曝| 一二区成人影院电影网| 日本理伦片午夜理伦片| 欧美成人猛片aaaaaaa| 久久不见久久见免费影院www| jjzz日本护士| 国产偷拍自拍在线观看| 国产久色在线拍揄自揄拍| 久久国产成人精品av| 日韩插啊免费视频在线观看| 国产哟交泬泬视频在线播放 | 中文字幕人妻一区二区二区| 午夜精品久久久久久久99老熟妇| 一区二区三区乱码在线 | 欧洲| 99久久精品免费看国产情侣 | 久久精品国产亚洲av高清蜜臀 | 国产做国产爱免费视频| 国产一区二区三区四区五区vm| 亚洲av综合色区在线观看| 99久久婷婷国产一区| 国产日产欧产精品精品| 国产精品第1页在线观看| 亚洲日韩成人无码不卡网站| 在线人妻va中文字幕| 亚洲熟妇少妇任你躁在线观看无码 | 亚洲av永久无码精品成人| 日本午夜精品一区二区三区| 国产综合无码一区二区辣椒| 国产欧美精品一区二区三区,| 国产精品黄色av网站| 中文字幕av伊人av无码av| 亚洲av无码久久寂寞少妇| 久久亚洲中文字幕精品一区四| 国产精品国产三级国产av18| 天天躁日日躁狠狠躁欧美老妇小说| 藏春阁福利视频| 日本一区二区三区在线播放| 久久精品国产自产对白一区| 国产精品妇女一二三区| 欧美成人小视频| 久久精品日韩免费视频|