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

    起落架有限元梁模型建立及其靜力學(xué)分析

    2016-10-09 08:32:39婁銳折世強黃海清
    航空工程進展 2016年3期
    關(guān)鍵詞:外筒起落架航向

    婁銳,折世強,黃海清

    (中航飛機起落架有限責(zé)任公司 工程技術(shù)研究中心,長沙 410200)

    ?

    起落架有限元梁模型建立及其靜力學(xué)分析

    婁銳,折世強,黃海清

    (中航飛機起落架有限責(zé)任公司 工程技術(shù)研究中心,長沙410200)

    為了獲得起落架載荷傳遞以及評估幾何非線性和靜不定對起落架載荷傳遞的影響,同時探索起落架結(jié)構(gòu)有限元梁模型的建模方法,應(yīng)用HyperMesh和Ansys聯(lián)合仿真的方法搭建某型具有局部靜不定特征的前起落架、主起落架有限元梁模型;使用此模型進行三種載荷工況的靜力學(xué)計算,獲得考慮靜不定和幾何非線性的起落架接頭和截面載荷,以及在單位載荷作用下起落架的航向、側(cè)向和扭轉(zhuǎn)剛度,并與傳統(tǒng)方法計算結(jié)果進行對比。結(jié)果表明:有限元方法和傳統(tǒng)方法獲得的接頭載荷有差異,部分差異超過了±10%,最大差異達到30%;計算前起落架外筒橫梁、斜梁截面內(nèi)力的傳統(tǒng)方法及其簡化具有合理性,對起落架支柱設(shè)計具有指導(dǎo)意義。

    起落架;有限元梁模型;靜不定;幾何非線性

    0 引 言

    獲得飛機地面載荷后,根據(jù)起落架結(jié)構(gòu)方案,準(zhǔn)確計算其載荷傳遞,包括接頭、運動副載荷以及零件截面內(nèi)力,對起落架設(shè)計和優(yōu)化具有重要意義。雖然起落架是單傳力形式,但某些位置存在局部靜不定并影響載荷傳遞,一般的處理方法是對結(jié)構(gòu)做出某種程度的假設(shè)和簡化,將靜不定問題轉(zhuǎn)化為靜定問題,從而可以通過平衡方程解算,所獲得的載荷是一個近似值。同時,起落架受載發(fā)生變形會產(chǎn)生附加彎矩,進一步影響載荷傳遞,行業(yè)規(guī)范對外載荷引起變形并對內(nèi)部載荷分布變化的影響有明確規(guī)定[1]。

    高澤迥[2]、劉銳琛[3]給出了計算靜定起落架結(jié)構(gòu)載荷傳遞的傳統(tǒng)方法,并指出有限元可以解決靜不定和變形附加載荷問題,但是并未給出有限元計算載荷傳遞的建模方法,以及傳統(tǒng)方法與有限元方法結(jié)果的差異;邵永起[4]研究了考慮起落架結(jié)構(gòu)變形的接頭載荷和系統(tǒng)內(nèi)力的計算方法,但此方法依賴于起落架支柱變形曲線的確定,并指出變形曲線應(yīng)由梁元素有限元法獲得;有限元方法針對起落架的靜力學(xué)應(yīng)用局限于零件和裝配件的三維實體有限元分析[5-6],研究目的是獲得應(yīng)力細(xì)節(jié),并不是載荷傳遞和變形模式。

    本文針對某型飛機前起落架、主起落架探索建立有限元梁模型,并進行幾何非線性靜力學(xué)分析和應(yīng)用,包括接頭載荷、零件截面內(nèi)力、初步屈曲分析、剛度隨行程的關(guān)系,從而獲得考慮局部靜不定和幾何非線性作用的起落架載荷傳遞規(guī)律,并針對傳統(tǒng)方法和有限元方法的計算結(jié)果進行對比分析,從而獲得兩種方法對載荷的影響。

    1 有限元方法

    有限元靜力學(xué)問題的平衡方程[7]為

    KU=Fa

    (1)

    式中:K為剛度矩陣;U為自由度列陣;Fa為外載荷列陣。

    K中元素包含所有靜定、靜不定結(jié)構(gòu)的剛度,通過求解式(1)可以獲得自由度列陣U,根據(jù)U計算接頭、運動副載荷和截面內(nèi)力。

    在Ansys軟件中,幾何非線性的求解使用N-R平衡迭代方法,幾何非線性的有限元控制方程[8-9]為

    (2)

    Ui+1=Ui+ΔUi

    (3)

    求解過程如下:通過式(2)獲得位移增量ΔUi,由式(3)得到新的節(jié)點位移Ui+1,式(2)右邊稱為殘余載荷列陣R,當(dāng)R所規(guī)定的范數(shù)滿足精度要求時,稱為載荷收斂,并獲得幾何非線性問題的收斂解U。

    2 起落架結(jié)構(gòu)及計算其結(jié)構(gòu)載荷傳遞的傳統(tǒng)方法

    某型飛機前起落架、主起落架結(jié)構(gòu)示意圖如圖1所示。前起落架主要組成部分如圖1(a)所示。支柱和阻力撐桿在航向所形成的三角形用于承受航向載荷,垂向載荷由支柱本身承受,外筒自身側(cè)向的三角形用于承受側(cè)向載荷。

    主起落架為上單翼翼根安裝的特型起落架,主要組成部分如圖1(b)所示。主梁用于連接航向桿、垂向桿和支柱,使三者成為一體,收放支架與機身連接,為了實現(xiàn)支柱的錐面收放,支柱上半部分連接叉形件,叉形件另一端與收放支架相連,為了增加收放支架的側(cè)向剛度,主梁和收放支架之間配置拉桿,下位鎖通過鎖臂和鎖殼之間的插銷實現(xiàn),航向桿和下位鎖承受航向載荷,垂向桿和主梁承受垂向載荷,下位鎖和主梁承受側(cè)向載荷。

    (a) 前起落架

    (b) 主起落架 圖1 起落架結(jié)構(gòu)示意圖Fig.1 Nose and main landing gear structure sketch

    前起落架外筒為三角結(jié)構(gòu),如圖2(a)所示。主接頭受到Fx、Fy和Fz載荷作用,載荷在外筒橫梁和斜梁傳遞具有靜不定特征,傳統(tǒng)方法需進行簡化和假設(shè),如圖2(b)所示。斜梁軸線方向較其他方向剛度大,可簡化為二力桿載荷P,橫梁與支柱軸線固定連接簡化為垂直于三角平面的旋轉(zhuǎn)副,二力桿載荷P由式(4)得到:

    P=(Fyd2+Fzd1)/d3

    (4)

    (a) 前起落架外筒接頭受力

    (b) 前起落架外筒橫梁、斜梁受載簡化 圖2 前起落架外筒傳力簡化示意圖Fig.2 Analysis simplification of load path of nose gear main fitting

    從圖1(b)可以看出:主起落架的主承力結(jié)構(gòu)與叉形件分開布置,叉形件與坐標(biāo)系的方向余弦全不為零,三個方向都能承受載荷,使得主起落架載荷傳遞具有靜不定特征,傳統(tǒng)方法需要進行簡化和假設(shè):將叉形件與收放支架斷開,同時將主梁與收放支架間的拉桿斷開。

    主起落架鎖銷受到的航向載荷Fx1和Fx2通過鎖殼傳遞至接頭也是一個局部靜不定問題,傳統(tǒng)方法將航向載荷在兩接頭處平分,同時假設(shè)側(cè)向載荷Fy1和Fy2在接頭處不產(chǎn)生航向載荷,如圖3所示。

    圖3 主起落架鎖殼受力與邊界Fig.3 Force analysis and boundary condition of main gear down lock shell

    3 梁模型建立

    梁模型坐標(biāo)系:原點任意,x軸逆航向為正,z軸沿支柱軸線向上為正,右手定則確定y軸。梁模型假設(shè):垂直載荷通過活塞桿完全傳遞給柱塞再傳遞給外筒。

    應(yīng)用HyperMesh前處理、Ansys求解的方法進行有限元梁模型的建立和求解[10]。對于具有細(xì)長形狀的零件,例如活塞桿、柱塞、航向桿應(yīng)用Beam188梁單元建模,梁單元的剛度由反映結(jié)構(gòu)截面特性的截面數(shù)據(jù)和材料屬性確定;對于具有復(fù)雜形狀的零件,例如叉形件、鎖殼與鎖臂則使用Shell281單元建模,殼單元的剛度由其厚度實常數(shù)和材料屬性確定;梁與梁之間的連接關(guān)系通過耦合節(jié)點自由度Couple DOF方式實現(xiàn),活塞桿與下支撐屬于此類連接方式;梁與殼的連接方式通過MPC裝配技術(shù)使用Conta175和Targe170單元并放開相關(guān)自由度實現(xiàn),鎖殼與鎖銷的連接屬于該類型,如圖4所示,與鎖銷配合的殼單元節(jié)點建立Conta175單元,鎖銷梁單元節(jié)點建立Targe170單元,機輪輪胎使用質(zhì)量單元Mass21模擬慣量。最后,定義接頭位移邊界條件,應(yīng)特別注意前起落架外筒與機身、前起落架上阻力撐桿與機身的側(cè)向約束屬于單邊約束,通過Conta178單元實現(xiàn)單邊約束定義。

    圖4 鎖殼與鎖銷的有限元連接模型Fig.4 FEA model of lock shell and pin connection

    根據(jù)前起落架結(jié)構(gòu)特點,前起落架使用全梁模型進行建模,如圖5所示。

    圖5 前起落架三維梁模型Fig.5 FEA beam model of nose gear

    根據(jù)主起落架結(jié)構(gòu)特點,主起落架采用梁殼混合模型,如圖6所示。

    圖6 主起落架梁殼模型Fig.6 FEA beam and shell mix model of main gear

    4 梁模型靜力學(xué)應(yīng)用

    獲得起落架接頭載荷、零件截面內(nèi)力、運動副載荷以及初步屈曲載荷是計算起落架強度和穩(wěn)定性的前提,獲得起落架輪軸中心航向、側(cè)向和扭轉(zhuǎn)剛度特性是起落架變形計算的前提。

    4.1接頭載荷

    起落架梁模型可以獲取考慮局部靜不定和幾何非線性條件下的接頭載荷。以主起落架為對象針對典型工況剎車MB64和轉(zhuǎn)彎MT64(具體載荷情況如表1所示)對接頭載荷進行比較,用于對比分析傳統(tǒng)方法與有限元方法接頭載荷計算結(jié)果的差異以及基于有限元方法的幾何非線性對接頭載荷的影響,計算結(jié)果如表2所示。

    傳統(tǒng)方法接頭載荷相對于有限元方法接頭載荷的相對變化量如圖7所示。由于MT64工況下F2、Fnx、Fmx、Fkx等于0,F(xiàn)ky載荷符號出現(xiàn)反向情況,不參與比較。

    基于有限元方法的不考慮幾何非線性的接頭載荷相對于考慮幾何非線性的接頭載荷相對變化量如圖8所示。由于MT64工況下載荷Fkx較小,不參與比較。

    圖7 不同工況條件下不同計算方法的 主起落架接頭載荷差異Fig.7 Main gear attachment load comparison between different models under different load cases

    圖8 不同工況條件下幾何非線性對主起落架 接頭載荷影響Fig.8 Geometry nonlinear influence on main gear attachment load under different load cases

    從圖7可以看出:部分接頭載荷差異超過了±10%,最大差異超過30%;大部分接頭載荷相對變化量為正,表明叉形件和拉桿承受了一部分地面載荷,起到減載作用;其余接頭載荷相對變化量為負(fù),表明幾何非線性變形引起的增載超過了靜不定的減載。

    從圖8可以看出:大部分接頭載荷相對變化量為負(fù),表明幾何非線性對大部分接頭載荷起到增加作用;接頭載荷差異介于±10%之間;MT64工況下載荷Fky差異達到-73%((6 442-23 903)/23 903×100%),并考慮到傳統(tǒng)方法與有限元方法Fky載荷符號出現(xiàn)反向的情況,表明載荷Fky對側(cè)向地面載荷較敏感。

    應(yīng)注意表2中MT64工況鎖殼x方向載荷Fnx和Fmx,因為轉(zhuǎn)彎工況輪軸Plx=Prx=0,所以傳統(tǒng)方法Fnx=0、Fmx=0,而有限元方法給出較大互為反向的x方向載荷,這是由于鎖殼受到的y方向載荷引起x方向變形引起的。

    表1 前起落架和主起落架典型工況

    表2 不同工況下、不同計算模型的主起落架接頭載荷

    4.2截面內(nèi)力

    起落架梁模型可以獲取考慮局部靜不定和幾何非線性條件下的截面內(nèi)力。針對前起落架典型工況起轉(zhuǎn)NS55,具體載荷情況如表1所示,對外筒橫梁、斜梁截面內(nèi)力進行比較,用于對比傳統(tǒng)方法與有限元方法處理局部靜不定的截面內(nèi)力計算結(jié)果差異,截面位置如圖9所示。根據(jù)式(4),可得傳統(tǒng)方法的二力桿載荷P=108 118 N,正號表示二力桿受拉。

    圖9 前起落架外筒不同截面位置與截面坐標(biāo)系Fig.9 Section locations and section coordinate systems of nose gear main fitting

    通過輸出單元截面內(nèi)力獲得梁模型結(jié)果并與傳統(tǒng)方法結(jié)果比較,如表3所示。

    表3 不同計算模型的前起落架外筒不同截面內(nèi)力

    從表3可以看出:傳統(tǒng)方法斜梁對于截面1只有軸向載荷Fx,而使用梁模型得出的截面1載荷除了具有較大的軸向載荷外,其余方向還承受了一定的載荷;對于截面1,傳統(tǒng)方法和梁模型獲得的截面內(nèi)力在截面設(shè)計的主要載荷Fx方面比較一致;對于截面2,傳統(tǒng)方法和梁模型獲得的截面內(nèi)力在截面設(shè)計的主要載荷Fx、Fy和Mz方面比較一致。因此,計算前起落架外筒橫梁、斜梁截面內(nèi)力的傳統(tǒng)方法及其簡化具有合理性,對前起落架支柱設(shè)計具有指導(dǎo)意義。

    4.3運動副載荷

    耦合兩節(jié)點自由度方式形成的運動副,以及MPC裝配技術(shù)形成的運動副可以通過“節(jié)點載荷”獲取運動副載荷,本文不展開討論。

    4.4特征值屈曲分析

    以主起落架為例,應(yīng)用主起落架有限元梁模型,進行全局特征屈曲分析[9],獲得載荷因子,為具有較大長細(xì)比的受壓零件設(shè)計提供初步屈曲分析,特征值屈曲分析要求有限元模型是線性的。

    在剎車工況MB64下,載荷因子f=2.146 49時的特征屈曲模式如圖10所示,為航向桿受壓失穩(wěn),則在此限制載荷工況下,航向桿受壓失穩(wěn)安全裕度MoS=2.146 49-1=1.146 49,由表2可以獲得航向桿受到的壓載荷為F2=-337 140 N,則航向桿的受壓失穩(wěn)臨界載荷為Fcr=f×F2=2.146 49×(-337 140)=-723 668 N,此載荷可用于航向桿受壓臨界設(shè)計條件。

    圖10 主起落架航向桿失穩(wěn)模式Fig.10 Longitudinal rod buckling mode of main gear

    4.5剛度特性

    以主起落架為例,分別建立行程0、0.2和0.4 m的梁模型,各自梁模型分別獲得單位載荷Fi作用下對應(yīng)的變形δi,其中i可以表示航向、側(cè)向和扭轉(zhuǎn)。根據(jù)公式Ki=Fi/δi分別計算航向、側(cè)向、扭轉(zhuǎn)剛度,假設(shè)剛度是行程的二次函數(shù),再通過二次擬合方法獲得全行程的剛度數(shù)據(jù),如圖11所示。

    圖11 主起落架輪軸中心剛度隨行程變化Fig.11 Main gear wheel axle center stiffness vs shock absorber travel

    5 結(jié) 論

    (1) 本文提供了一種新的載荷傳遞計算方法,所使用的有限元方法在獲得起落架載荷傳遞和整體剛度方面是一種有益的嘗試,此方法相對于傳統(tǒng)方法考慮了局部靜不定和幾何非線性條件對載荷傳遞的影響。

    (2) 本文比較了有限元方法和傳統(tǒng)方法獲得的接頭載荷和截面內(nèi)力,發(fā)現(xiàn)對于接頭載荷不同方法有較大差異,最大差異達到30%,這對起落架設(shè)計和確定零件截面是不利的,建議使用有限元方法獲得的載荷用于起落架設(shè)計。

    (3) 本文提出的有限元方法為仿真計算,計算結(jié)果應(yīng)通過靜力試驗來驗證。

    [1] 中國民用航空局. CCAR25-R4-2008運輸類飛機適航標(biāo)準(zhǔn)[S]. 北京: 中國民用航空局, 2008: 29.

    Civil Aviation Administration of China. CCAR25-R4-2008 Airworthiness standards: transport category airplanes[S]. Beijing: Civil Aviation Administration of China, 2008: 29.(in Chinese)

    [2] 高澤迥. 飛機設(shè)計手冊: 第14分冊——起飛著陸系統(tǒng)設(shè)計[M]. 北京: 航空工業(yè)出版社, 2002: 615-617.

    Gao Zejiong. Aircraft Design Manual: Part 14-Take off & landing system design[M]. Beijing: Aviation Industry Press, 2002: 615-617.(in Chinese)

    [3] 劉銳琛. 飛機起落架強度設(shè)計指南[M]. 成都: 四川科學(xué)技術(shù)出版社, 1989: 402-403.

    Liu Ruichen. Design guides for aircraft landing gear strength[M]. Chengdu: Sichuan Science and Technology Press, 1989: 402-403.(in Chinese)

    [4] 邵永起. 起落架彈性變形對結(jié)構(gòu)強度的影響[J]. 飛機設(shè)計, 2001(2): 38-43.

    Shao Yongqi. Landing gear elastic deformation impact on structure strength[J]. Aircraft Design, 2001(2): 38-43.(in Chinese)

    [5] 汪洪峰, 汪建利, 宋娓娓, 等. 飛機起落架主起外筒起飛降落過程中的應(yīng)力分布[J]. 機械設(shè)計, 2013, 30(6): 61-64.

    Wang Hongfeng,Wang Jianli,Song Weiwei,et al. Stress distribution of main outer cylinder of landing gear in process of aircraft tanking off and landing[J]. Journal of Machine Design, 2013, 30(6): 61-64.(in Chinese)

    [6] 黃銘敏, 左敦穩(wěn), 黎向鋒, 等. 四輪并式起落架焊接主起輪軸受載模型及分析[J]. 中國機械工程, 2012, 23(17): 2071-2074.

    Huang Mingmin, Zuo Dunwen, Li Xiangfeng, et al. Model establish and numerical analysis of welded wheel axle of four-wheel-juxtaposed landing gear[J]. China Mechanical Engineering, 2012, 23(17): 2071-2074.(in Chinese)

    [7] 曾攀. 有限元分析及應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2004: 96-98.

    Zeng Pan. Finite element analysis and applications[M]. Beijing: Tsinghua University Press, 2004: 96-98.(in Chinese)

    [8] 王新敏. 桿系結(jié)構(gòu)的幾何非線性分析綜述[J]. 石家莊鐵道學(xué)院學(xué)報, 1993, 6(3): 80-84.

    Wang Xinmin. On the geometric nonlinear analysis theory of frames[J]. Journal of Shijiazhuang Railway Institute, 1993, 6(3): 80-84.(in Chinese)

    [9] Ansys. Ansys help[M]. USA: Ansys company, 2013.

    [10] 賀李平, 龍凱, 肖介平, 等. ANSYS 13.0與Hyper Mesh 11.0聯(lián)合仿真有限元分析[M]. 北京: 機械工業(yè)出版社, 2012: 39-51.

    He Liping, Long Kai, Xiao Jieping, et al. ANSYS 13.0 and Hyper Mesh 11.0 cosimulation of finite element analysis[M]. Beijing: China Machine Press, 2012: 39-51.(in Chinese)

    (編輯:馬文靜)

    Establishing FEA Beam Model of Landing Gear and Its Analysis of Statics

    Lou Rui, She Shiqiang, Huang Haiqing

    (Engineering Technology Research Center, AVIC Landing Gear Advanced Manufacturing Corporation, Ltd., Changsha 410200, China)

    In order to obtain load path of landing gear and evaluate the load path impact of geometry nonlinear and static indeterminacy, meanwhile to explore method of building FEA beam model of landing gear, the FEA beam models of specific type nose and main landing gear with feature of local static indeterminacy are built by way of HyperMesh and Ansys co-simulation. Static calculation of three load cases is conducted using the FEA model. The joint load and section load of landing gear with feature of static indeterminacy and geometry nonlinear are obtained, and longitudinal, lateral and torsion stiffness of the gear are obtained under unit load, and compared with that by using the traditional method. Results show that there is difference between the joint load of the two methods. Some difference is exceeded by ±10% and the biggest difference achieved by 30%. The traditional method and simplification of calculating section load of level beam and inclined beam of nose strut is reasonable and is of guidance meaning in designing nose strut.

    landing gear; finite element beam model; static indeterminacy; geometry nonlinear

    2016-06-04;

    2016-07-23

    婁銳,lourui607@163.com

    1674-8190(2016)03-325-07

    V226

    A

    10.16615/j.cnki.1674-8190.2016.03.009

    婁銳(1985-),男,碩士,工程師。主要研究方向:起落架緩沖性能和載荷分析。

    折世強(1968-),男,研究員。主要研究方向:起落架緩沖性能和載荷分析。

    黃海清(1975-),男,高級工程師。主要研究方向:起落架強度計算與設(shè)計。

    猜你喜歡
    外筒起落架航向
    一種預(yù)熱型塑料吹膜機
    煤礦用掘進機截割臂伸縮外筒強度分析與優(yōu)化設(shè)計*
    知坐標(biāo),明航向
    尾管固井回接筒銑錐限位清銑及修復(fù)技術(shù)
    飛機秘密檔案
    一種多旋翼無人機起落架快速插接結(jié)構(gòu)
    液壓缸端蓋與外筒卡圈連接結(jié)構(gòu)探索
    考慮幾何限制的航向道模式設(shè)計
    基于干擾觀測器的船舶系統(tǒng)航向Backstepping 控制
    電子制作(2017年24期)2017-02-02 07:14:16
    使命:引領(lǐng)航向與保持穩(wěn)定
    法大研究生(2015年2期)2015-02-27 10:13:55
    亚洲成人精品中文字幕电影| 日韩大码丰满熟妇| 免费女性裸体啪啪无遮挡网站| 俄罗斯特黄特色一大片| 99riav亚洲国产免费| 老司机午夜福利在线观看视频| 国产精品一区二区免费欧美| 久热爱精品视频在线9| 国产野战对白在线观看| 香蕉丝袜av| 久久性视频一级片| 变态另类丝袜制服| 午夜福利成人在线免费观看| 99精品久久久久人妻精品| 国产色视频综合| av在线天堂中文字幕| 国产高清激情床上av| 精品一区二区三区视频在线观看免费| 又黄又爽又免费观看的视频| 中文字幕av电影在线播放| 久久 成人 亚洲| 又紧又爽又黄一区二区| 国产亚洲精品久久久久久毛片| 亚洲中文日韩欧美视频| 最近最新中文字幕大全电影3 | 91在线观看av| 黄色毛片三级朝国网站| 波多野结衣巨乳人妻| 搞女人的毛片| 亚洲少妇的诱惑av| 欧美午夜高清在线| 国产三级黄色录像| 久久人妻熟女aⅴ| 国产伦人伦偷精品视频| 久久久久国内视频| 69精品国产乱码久久久| 国产欧美日韩综合在线一区二区| 国产单亲对白刺激| 日本在线视频免费播放| 日日夜夜操网爽| 国产精品电影一区二区三区| 一本综合久久免费| 侵犯人妻中文字幕一二三四区| 在线观看66精品国产| 国产视频一区二区在线看| 麻豆成人av在线观看| 黄色成人免费大全| 欧美一区二区精品小视频在线| 少妇粗大呻吟视频| 男女午夜视频在线观看| 少妇被粗大的猛进出69影院| 在线天堂中文资源库| 国产亚洲欧美在线一区二区| 精品免费久久久久久久清纯| 在线观看午夜福利视频| 亚洲五月婷婷丁香| 中文字幕色久视频| 免费无遮挡裸体视频| 97人妻精品一区二区三区麻豆 | 亚洲欧美一区二区三区黑人| 国产av在哪里看| 欧美黑人欧美精品刺激| 久99久视频精品免费| 91老司机精品| 精品国产亚洲在线| 亚洲成av片中文字幕在线观看| 国产精品久久视频播放| 久久久久久久久免费视频了| 一进一出抽搐gif免费好疼| 午夜福利,免费看| 一级片免费观看大全| 精品日产1卡2卡| 熟妇人妻久久中文字幕3abv| 男女下面进入的视频免费午夜 | 在线观看午夜福利视频| 夜夜躁狠狠躁天天躁| 在线观看www视频免费| 99热只有精品国产| 99精品久久久久人妻精品| 久久青草综合色| 激情视频va一区二区三区| 国产免费av片在线观看野外av| 一夜夜www| 亚洲免费av在线视频| 精品国产国语对白av| 天天躁夜夜躁狠狠躁躁| 视频区欧美日本亚洲| 在线观看日韩欧美| 成人精品一区二区免费| 最好的美女福利视频网| 国产成人精品无人区| 狂野欧美激情性xxxx| 香蕉久久夜色| 久久人人精品亚洲av| 亚洲人成77777在线视频| 99精品久久久久人妻精品| 精品国内亚洲2022精品成人| 好男人电影高清在线观看| 国产黄a三级三级三级人| 一区在线观看完整版| 亚洲欧美日韩高清在线视频| 18禁观看日本| 国产精品av久久久久免费| 国产麻豆69| 亚洲色图综合在线观看| 国产伦一二天堂av在线观看| www国产在线视频色| 男人舔女人下体高潮全视频| 午夜影院日韩av| 人人妻人人澡欧美一区二区 | 香蕉国产在线看| 性少妇av在线| 正在播放国产对白刺激| 一级毛片精品| 午夜免费成人在线视频| 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 亚洲一区中文字幕在线| 99国产极品粉嫩在线观看| 国产又色又爽无遮挡免费看| 免费高清在线观看日韩| 后天国语完整版免费观看| 国产麻豆成人av免费视频| 99久久综合精品五月天人人| 波多野结衣巨乳人妻| 欧美不卡视频在线免费观看 | 国产单亲对白刺激| 中文字幕久久专区| 精品少妇一区二区三区视频日本电影| 国产又色又爽无遮挡免费看| 久久中文字幕一级| 久久久精品国产亚洲av高清涩受| 久久久久久久久免费视频了| 精品福利观看| 精品久久久久久,| 91精品三级在线观看| 51午夜福利影视在线观看| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| ponron亚洲| 国内毛片毛片毛片毛片毛片| 国内久久婷婷六月综合欲色啪| 午夜激情av网站| 亚洲成人国产一区在线观看| aaaaa片日本免费| 久久久久国产精品人妻aⅴ院| av视频免费观看在线观看| 成人手机av| 精品一区二区三区av网在线观看| 人人妻人人澡欧美一区二区 | 亚洲国产欧美日韩在线播放| 亚洲精品中文字幕一二三四区| 亚洲国产精品久久男人天堂| 天天躁夜夜躁狠狠躁躁| 亚洲五月天丁香| 黄色a级毛片大全视频| 成人国产一区最新在线观看| 69av精品久久久久久| 两个人看的免费小视频| 午夜视频精品福利| 日韩av在线大香蕉| 9色porny在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 一本久久中文字幕| 老司机深夜福利视频在线观看| 在线观看免费视频网站a站| 国产精华一区二区三区| 超碰成人久久| 老汉色av国产亚洲站长工具| 日本vs欧美在线观看视频| 亚洲一区中文字幕在线| 欧美激情极品国产一区二区三区| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| a级毛片在线看网站| 欧美一级a爱片免费观看看 | 少妇粗大呻吟视频| 电影成人av| 国产99久久九九免费精品| 国产男靠女视频免费网站| 人成视频在线观看免费观看| 国产1区2区3区精品| 日韩精品青青久久久久久| 国产精品日韩av在线免费观看 | 午夜影院日韩av| 久久人妻福利社区极品人妻图片| 丰满人妻熟妇乱又伦精品不卡| 纯流量卡能插随身wifi吗| 亚洲 欧美 日韩 在线 免费| 视频区欧美日本亚洲| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 久久久精品欧美日韩精品| √禁漫天堂资源中文www| 亚洲电影在线观看av| 天天添夜夜摸| 精品久久久久久久久久免费视频| 国产亚洲av高清不卡| avwww免费| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 麻豆久久精品国产亚洲av| 69av精品久久久久久| 欧美色欧美亚洲另类二区 | 午夜福利欧美成人| 日本 欧美在线| 中文字幕最新亚洲高清| 精品国产亚洲在线| 日韩欧美一区二区三区在线观看| 亚洲七黄色美女视频| 欧美色欧美亚洲另类二区 | 高清黄色对白视频在线免费看| 99精品久久久久人妻精品| 亚洲精品国产色婷婷电影| 久久精品国产综合久久久| 国产午夜福利久久久久久| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 99久久精品国产亚洲精品| 热re99久久国产66热| 妹子高潮喷水视频| 亚洲,欧美精品.| 久久久久九九精品影院| 国产精品久久电影中文字幕| 亚洲午夜理论影院| 99国产极品粉嫩在线观看| 99国产综合亚洲精品| 亚洲国产精品成人综合色| 午夜福利成人在线免费观看| 欧美激情久久久久久爽电影 | 国产激情久久老熟女| 国产精品乱码一区二三区的特点 | 可以免费在线观看a视频的电影网站| 亚洲欧美精品综合久久99| 妹子高潮喷水视频| 亚洲成人免费电影在线观看| 男女做爰动态图高潮gif福利片 | 国产成人一区二区三区免费视频网站| 一区二区三区国产精品乱码| av天堂在线播放| 俄罗斯特黄特色一大片| 少妇的丰满在线观看| 成人国产综合亚洲| 午夜福利欧美成人| 国产av在哪里看| 免费高清视频大片| 成人手机av| 午夜日韩欧美国产| 日韩国内少妇激情av| 一区福利在线观看| 黄片大片在线免费观看| 久久青草综合色| 亚洲精品一区av在线观看| 久热这里只有精品99| 日本 欧美在线| 精品久久久久久久久久免费视频| 亚洲欧美激情在线| 黑丝袜美女国产一区| 久久伊人香网站| 免费在线观看影片大全网站| 国产一区二区三区综合在线观看| 性少妇av在线| 国产精品亚洲美女久久久| 99re在线观看精品视频| 一级片免费观看大全| 亚洲欧美精品综合久久99| 免费无遮挡裸体视频| 亚洲激情在线av| 两性夫妻黄色片| 中文字幕av电影在线播放| 婷婷精品国产亚洲av在线| 99国产综合亚洲精品| 成人三级做爰电影| 国产精品一区二区精品视频观看| 国产一区在线观看成人免费| av超薄肉色丝袜交足视频| 最近最新中文字幕大全电影3 | 久久国产精品人妻蜜桃| 12—13女人毛片做爰片一| 久久久精品欧美日韩精品| 精品午夜福利视频在线观看一区| 欧美精品亚洲一区二区| 国产精品二区激情视频| 非洲黑人性xxxx精品又粗又长| 精品高清国产在线一区| 天天躁狠狠躁夜夜躁狠狠躁| 制服丝袜大香蕉在线| 美女免费视频网站| 久久久久久久久久久久大奶| 美女扒开内裤让男人捅视频| 欧美绝顶高潮抽搐喷水| 国产不卡一卡二| 日日干狠狠操夜夜爽| 一区在线观看完整版| 我的亚洲天堂| 国产成人欧美| 88av欧美| 女警被强在线播放| 一本大道久久a久久精品| 国产主播在线观看一区二区| 桃红色精品国产亚洲av| 黄色毛片三级朝国网站| 91成人精品电影| 亚洲av熟女| xxx96com| 亚洲全国av大片| 搡老妇女老女人老熟妇| www.熟女人妻精品国产| 国产午夜福利久久久久久| 亚洲无线在线观看| 91麻豆av在线| av网站免费在线观看视频| 黑人操中国人逼视频| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品久久男人天堂| 久久亚洲精品不卡| 久久草成人影院| 日本五十路高清| 可以在线观看的亚洲视频| 妹子高潮喷水视频| 亚洲成av人片免费观看| 亚洲中文字幕日韩| 一本久久中文字幕| 日本黄色视频三级网站网址| 黄片大片在线免费观看| 亚洲国产高清在线一区二区三 | 欧美激情高清一区二区三区| 香蕉久久夜色| av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av | 亚洲avbb在线观看| 国产一区在线观看成人免费| 国产亚洲欧美精品永久| 国产精品香港三级国产av潘金莲| 午夜福利成人在线免费观看| 男女下面进入的视频免费午夜 | 两个人看的免费小视频| 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 亚洲成av片中文字幕在线观看| 级片在线观看| 精品乱码久久久久久99久播| 成人18禁高潮啪啪吃奶动态图| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 亚洲片人在线观看| 欧美中文综合在线视频| 久久婷婷人人爽人人干人人爱 | 亚洲,欧美精品.| 亚洲精品国产色婷婷电影| 久久亚洲精品不卡| 在线播放国产精品三级| 宅男免费午夜| 黄色a级毛片大全视频| 亚洲第一青青草原| 亚洲在线自拍视频| 变态另类丝袜制服| 99国产精品99久久久久| 欧美激情 高清一区二区三区| 亚洲免费av在线视频| 99久久99久久久精品蜜桃| 热99re8久久精品国产| 中文字幕av电影在线播放| 精品人妻在线不人妻| 亚洲欧美激情综合另类| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女| √禁漫天堂资源中文www| 亚洲全国av大片| 黄色视频不卡| 欧美+亚洲+日韩+国产| 国产精品香港三级国产av潘金莲| 免费看十八禁软件| 久久人妻av系列| 亚洲精品在线观看二区| 一级,二级,三级黄色视频| 精品卡一卡二卡四卡免费| 欧美午夜高清在线| 日韩一卡2卡3卡4卡2021年| 9热在线视频观看99| 成人亚洲精品一区在线观看| 黄色丝袜av网址大全| 国产精品影院久久| 这个男人来自地球电影免费观看| 男男h啪啪无遮挡| 国产av一区二区精品久久| 非洲黑人性xxxx精品又粗又长| 少妇 在线观看| 亚洲狠狠婷婷综合久久图片| 91大片在线观看| 97人妻精品一区二区三区麻豆 | 亚洲av成人av| 欧美一级毛片孕妇| 日本 av在线| 国内久久婷婷六月综合欲色啪| 日本 av在线| 在线播放国产精品三级| 日本精品一区二区三区蜜桃| 国产1区2区3区精品| 亚洲无线在线观看| 搞女人的毛片| av超薄肉色丝袜交足视频| 侵犯人妻中文字幕一二三四区| 亚洲人成电影观看| 国产一区二区三区视频了| 久久香蕉激情| 波多野结衣高清无吗| 日韩一卡2卡3卡4卡2021年| 久久人人精品亚洲av| 可以免费在线观看a视频的电影网站| 啦啦啦免费观看视频1| 国产精品久久电影中文字幕| 欧美激情 高清一区二区三区| 淫妇啪啪啪对白视频| 99在线人妻在线中文字幕| 日日夜夜操网爽| 午夜福利免费观看在线| 成人永久免费在线观看视频| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久久久毛片| 国产色视频综合| 欧美一区二区精品小视频在线| 久久久精品国产亚洲av高清涩受| 在线视频色国产色| 欧美日韩亚洲国产一区二区在线观看| 亚洲视频免费观看视频| 一级毛片女人18水好多| 日韩视频一区二区在线观看| 欧美绝顶高潮抽搐喷水| 中国美女看黄片| 妹子高潮喷水视频| 午夜老司机福利片| 啦啦啦免费观看视频1| svipshipincom国产片| 日韩欧美三级三区| 成人欧美大片| 免费看美女性在线毛片视频| 9热在线视频观看99| 十八禁人妻一区二区| 麻豆久久精品国产亚洲av| 人人妻人人澡人人看| 91成人精品电影| 欧美另类亚洲清纯唯美| 人人妻人人爽人人添夜夜欢视频| 一区二区三区国产精品乱码| 亚洲av第一区精品v没综合| 国产精品一区二区精品视频观看| 精品久久久久久久毛片微露脸| 国产亚洲精品久久久久久毛片| 国产精品自产拍在线观看55亚洲| 淫妇啪啪啪对白视频| 国产男靠女视频免费网站| 日韩高清综合在线| 成年版毛片免费区| 亚洲人成77777在线视频| 亚洲情色 制服丝袜| 日韩精品免费视频一区二区三区| 国产极品粉嫩免费观看在线| 久热爱精品视频在线9| 久久精品国产亚洲av高清一级| 久久久久久久久中文| 免费高清视频大片| 亚洲男人的天堂狠狠| 日韩成人在线观看一区二区三区| 三级毛片av免费| 12—13女人毛片做爰片一| 一区二区三区精品91| 黄色成人免费大全| 国产亚洲精品久久久久5区| 国产亚洲精品av在线| 精品国产一区二区久久| 亚洲专区国产一区二区| cao死你这个sao货| 欧美国产精品va在线观看不卡| 久久中文字幕人妻熟女| 精品国产国语对白av| 国产亚洲欧美在线一区二区| 91在线观看av| 亚洲五月天丁香| 黄色 视频免费看| 日本撒尿小便嘘嘘汇集6| 午夜福利欧美成人| 高潮久久久久久久久久久不卡| 久久久久久久午夜电影| 一边摸一边抽搐一进一出视频| 99riav亚洲国产免费| 99久久久亚洲精品蜜臀av| 欧美黄色淫秽网站| 欧美午夜高清在线| 色综合婷婷激情| 欧美在线一区亚洲| 久久婷婷成人综合色麻豆| 国产精品秋霞免费鲁丝片| 日韩精品青青久久久久久| 很黄的视频免费| 神马国产精品三级电影在线观看 | 一区二区三区高清视频在线| 十八禁人妻一区二区| 女同久久另类99精品国产91| 热99re8久久精品国产| 在线观看一区二区三区| 免费在线观看日本一区| 午夜亚洲福利在线播放| 在线观看www视频免费| 精品国产一区二区久久| 两个人免费观看高清视频| 一级毛片女人18水好多| 亚洲欧美一区二区三区黑人| 一区二区三区激情视频| 操美女的视频在线观看| 精品国产美女av久久久久小说| 丁香六月欧美| 亚洲aⅴ乱码一区二区在线播放 | 女同久久另类99精品国产91| 精品国产一区二区久久| 久久精品国产亚洲av香蕉五月| 日日爽夜夜爽网站| 国产成人精品久久二区二区91| 亚洲中文字幕一区二区三区有码在线看 | 久久久久国产一级毛片高清牌| 18禁国产床啪视频网站| 欧美乱码精品一区二区三区| 久久久久九九精品影院| 免费看a级黄色片| 可以免费在线观看a视频的电影网站| 国产xxxxx性猛交| 成熟少妇高潮喷水视频| 成熟少妇高潮喷水视频| 少妇粗大呻吟视频| 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 欧美日韩乱码在线| 午夜福利影视在线免费观看| 黄频高清免费视频| 色婷婷久久久亚洲欧美| 久久人妻熟女aⅴ| 日韩国内少妇激情av| 欧美成人性av电影在线观看| 免费在线观看视频国产中文字幕亚洲| 搡老妇女老女人老熟妇| 又黄又爽又免费观看的视频| 亚洲黑人精品在线| 国产精品98久久久久久宅男小说| 麻豆av在线久日| 久久香蕉精品热| 两性夫妻黄色片| 麻豆久久精品国产亚洲av| 久99久视频精品免费| 亚洲专区中文字幕在线| 97超级碰碰碰精品色视频在线观看| 在线视频色国产色| 欧美日韩亚洲国产一区二区在线观看| 久久亚洲精品不卡| 最新美女视频免费是黄的| 欧美黄色片欧美黄色片| 国内精品久久久久精免费| 成人免费观看视频高清| 久久精品国产综合久久久| 亚洲av第一区精品v没综合| av免费在线观看网站| 午夜免费激情av| 老熟妇仑乱视频hdxx| 极品人妻少妇av视频| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 在线永久观看黄色视频| 黑人欧美特级aaaaaa片| 亚洲熟妇中文字幕五十中出| 日韩国内少妇激情av| 别揉我奶头~嗯~啊~动态视频| 午夜a级毛片| 午夜日韩欧美国产| 一本久久中文字幕| 丁香欧美五月| 亚洲精品av麻豆狂野| 国产精品久久视频播放| 最好的美女福利视频网| 午夜免费鲁丝| 久久久久九九精品影院| 欧美+亚洲+日韩+国产| 日韩欧美国产在线观看| 午夜亚洲福利在线播放| 国产精品亚洲一级av第二区| 欧美成狂野欧美在线观看| 亚洲av成人不卡在线观看播放网| 国产高清激情床上av| 女人精品久久久久毛片| 国产又色又爽无遮挡免费看| 国产成人av激情在线播放| 俄罗斯特黄特色一大片| 激情在线观看视频在线高清| 老鸭窝网址在线观看| 精品福利观看| 久久精品91无色码中文字幕| 夜夜看夜夜爽夜夜摸| 给我免费播放毛片高清在线观看| 国产高清视频在线播放一区| 国产麻豆69| 久久狼人影院| 国产精品一区二区在线不卡| 国产黄a三级三级三级人| 搡老熟女国产l中国老女人| 高清黄色对白视频在线免费看| 中文亚洲av片在线观看爽| 欧美日韩福利视频一区二区| 亚洲成人精品中文字幕电影| svipshipincom国产片| 免费高清视频大片| 搡老熟女国产l中国老女人| av网站免费在线观看视频| 国产99白浆流出| 亚洲欧美激情综合另类| 涩涩av久久男人的天堂| 精品日产1卡2卡| 亚洲成av片中文字幕在线观看|