999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

某乘用車前懸下控制臂結(jié)構(gòu)性能評(píng)估與優(yōu)化設(shè)計(jì)

2024-05-13 00:00:00游永忠朱海燕吳國(guó)棟楊春輝
關(guān)鍵詞:模態(tài)優(yōu)化

摘要:【目的】為了評(píng)估某乘用車前懸下控制臂的結(jié)構(gòu)性能是否滿足設(shè)計(jì)要求?!痉椒ā渴紫?,基于有限元方法建立下控制臂網(wǎng)格模型,對(duì)其進(jìn)行自由模態(tài)仿真,獲取其前三階頻率。然后,建立前懸架多體動(dòng)力學(xué)模型,獲得下控制臂外聯(lián)點(diǎn)的力和扭矩,并采用慣性釋放方法對(duì)其進(jìn)行強(qiáng)度性能仿真,獲取各個(gè)工況的應(yīng)力水平。最后,對(duì)下控制臂的厚度進(jìn)行優(yōu)化,獲取最優(yōu)結(jié)構(gòu)?!窘Y(jié)果】下擺臂的前三階頻率均處于外部激勵(lì)頻率范圍之外,滿足振動(dòng)特性要求;各工況的最大應(yīng)力均低于材料屈服值,滿足強(qiáng)度屬性要求;優(yōu)化后的結(jié)構(gòu)性能符合設(shè)計(jì)標(biāo)準(zhǔn),并且其重量減輕了12.5%?!窘Y(jié)論】研究結(jié)果能夠?yàn)橥惍a(chǎn)品設(shè)計(jì)及優(yōu)化提供有效的借鑒和參考,具有重要的實(shí)際工程意義。

關(guān)鍵詞:下控制臂;模態(tài);頻率;強(qiáng)度;優(yōu)化

中圖分類號(hào):U463.33 文獻(xiàn)標(biāo)志碼:A

本文引用格式:游永忠,朱海燕,吳國(guó)棟,等. 某乘用車前懸下控制臂結(jié)構(gòu)性能評(píng)估與優(yōu)化設(shè)計(jì)[J]. 華東交通大學(xué)學(xué)報(bào),2024,41(6):106-111.

【研究意義】汽車懸架系統(tǒng)是車輛的重要組成部分,其保證了車輛的操作穩(wěn)定性和舒適性。下控制臂是前麥弗遜懸架的導(dǎo)向與傳力部件,通過(guò)球鉸和襯套與車輪和副車架相連,在車輛行駛過(guò)程中承受來(lái)自不同方向的力和扭矩,可能引起異響和開(kāi)裂,其結(jié)構(gòu)性能直接影響車輛的可靠性和穩(wěn)定性,與此同時(shí)結(jié)構(gòu)的輕量化有助于提升整車的燃油動(dòng)力性,因此下控制臂設(shè)計(jì)需要兼顧各項(xiàng)性能要求和輕量化要求。

【研究進(jìn)展】謝穎等[1]采用慣性釋放方法對(duì)某下控制臂進(jìn)行強(qiáng)度仿真,并對(duì)其進(jìn)行結(jié)構(gòu)優(yōu)化,最終滿足了強(qiáng)度要求。梁國(guó)棟等[2]為了解決下擺臂開(kāi)裂問(wèn)題,搭建了下擺臂有限元模型和前懸架動(dòng)力學(xué)模型,對(duì)其進(jìn)行強(qiáng)度仿真和臺(tái)架測(cè)試,并且進(jìn)行了疲勞壽命預(yù)測(cè),采用集成優(yōu)化平臺(tái)對(duì)其進(jìn)行優(yōu)化設(shè)計(jì),最終得到了最優(yōu)的結(jié)構(gòu)參數(shù),重量減小了6.5%。徐峰[3]基于建立的麥弗遜懸架模型,獲取了外聯(lián)點(diǎn)的載荷,得到了其強(qiáng)度應(yīng)力分布和模態(tài)頻率,對(duì)拓?fù)鋬?yōu)化結(jié)果進(jìn)行重構(gòu),最終實(shí)現(xiàn)了輕量化。Song 等[4]分別采用響應(yīng)面和克里金近似模型對(duì)某擺臂進(jìn)行多目標(biāo)輕量化分析,對(duì)比了兩種方法的實(shí)際效果。

【創(chuàng)新特色】下控制臂載荷的提取涉及理論計(jì)算與多體動(dòng)力學(xué)模型的搭建,通過(guò)求解計(jì)算和分解獲取下擺臂外連點(diǎn)的力和扭矩。【關(guān)鍵問(wèn)題】為了分析某乘用車前懸下控制臂的結(jié)構(gòu)性能,首先基于有限元方法對(duì)其進(jìn)行自由模態(tài)仿真,獲取其低階模態(tài)頻率和振型,再對(duì)其前懸架進(jìn)行動(dòng)力學(xué)仿真分析,得到其極限強(qiáng)度的力和扭矩,以此對(duì)其進(jìn)行強(qiáng)度仿真,再對(duì)其結(jié)構(gòu)進(jìn)行優(yōu)化,最后達(dá)到輕量化的目的。

1 下控制臂模態(tài)性能評(píng)估

1.1 模態(tài)分析原理

下控制臂模態(tài)性能評(píng)估是指其固有模態(tài)和振型進(jìn)行計(jì)算,下控制臂的運(yùn)動(dòng)微分方程為[5-7]

Mx?(t)+Cx?(t)+Kx(t)= f (t) (1)

式中:M,C和K為下控制臂的質(zhì)量矩陣、阻尼矩陣和剛度矩陣;x(t) 和f (t) 為下控制臂的位移向量和載荷向量。

通過(guò)傅里葉變換即得

(-ω ) 2M+ ωC + K x(ω)= f (ω) (2)

式中:ω 和x(ω) 分別為下控制臂的固有頻率及其振型。

1.2 建立有限元模型

某乘用車前懸下控制臂是鑄鋁件,重量為1.6 kg,其前點(diǎn)和后點(diǎn)與副車架相連,其外點(diǎn)與轉(zhuǎn)向節(jié)相連,將下控制臂的三維模型導(dǎo)入Hypermesh 軟件[8-10]中,首先對(duì)其進(jìn)行幾何清理和簡(jiǎn)化,刪除對(duì)結(jié)構(gòu)性能影響較小的部件和特征。由于該下控制臂外表面復(fù)雜,難以使用六面體單元建模,因此采用3 mm的CTRIA3單元對(duì)其表面進(jìn)行網(wǎng)格劃分,并且確保其最小單元大于1 mm,最大單元小于4 mm,長(zhǎng)寬比小于4,翹曲度小于15,扭曲度小于45,然后生成四面體網(wǎng)格單元。下控制臂材料為A6082,其彈性模量為7.3×104 MPa,泊松比為0.33,密度為2.7×103 kg/m3,屈服值為300 MPa,創(chuàng)建各項(xiàng)同性的材料及其屬性,并賦予下控制臂,以此建立下控制臂有限元網(wǎng)格模型,如圖1 所示。其中網(wǎng)格單元總共145 123 個(gè),節(jié)點(diǎn)總共34 318個(gè)。

1.3 模態(tài)分析結(jié)果

下控制臂的模態(tài)性能是其固有特性,有限元仿真時(shí)采用自由狀態(tài),不施加邊界與約束。采用Nastran軟件[11-13]導(dǎo)入下控制臂的網(wǎng)格模型,設(shè)定頻率計(jì)算范圍為1~300 Hz,基于Lanczos 方法提取其固有頻率及其振型。下控制臂的低階模態(tài)對(duì)其模態(tài)性能起決定性作用,因此只提取下控制臂的前三階模態(tài)進(jìn)行評(píng)估。

圖2 為下控制臂第一階振型。由圖2 可知,下控制臂第一階振型是頻率為115.3 Hz 的彎曲模態(tài),前后點(diǎn)振幅較大,最大振幅為60.1 mm。

圖3為下控制臂第二階模態(tài)振型。由圖3可知,下控制臂第二階振型是頻率為159.0 Hz的一階扭轉(zhuǎn)模態(tài),外點(diǎn)振幅較大,最大變形為77.6 mm。

圖4為下控制臂第三階模態(tài)振型。由圖4可知,下控制臂第三階振型是頻率為198.7 Hz的一階扭轉(zhuǎn)模態(tài),后點(diǎn)振動(dòng)比較劇烈,最大振幅為59.6 mm。

該乘用車發(fā)動(dòng)機(jī)怠速為800 r/min,通過(guò)換算即得其一階頻率為26.6 Hz,輪胎聲腔模態(tài)頻率通常為220 Hz左右,因此下控制臂的固有頻率有效避開(kāi)了發(fā)動(dòng)機(jī)激勵(lì)頻率與輪胎聲腔激勵(lì)頻率,不會(huì)引起耦合共振,符合振動(dòng)特性設(shè)計(jì)要求。

2 下控制臂強(qiáng)度性能評(píng)估

2.1 前懸架動(dòng)力學(xué)仿真分析

車輛在行駛過(guò)程中承受縱向載荷、橫向載荷和垂向載荷,因此其強(qiáng)度工況可以分為縱向制動(dòng)、橫向轉(zhuǎn)彎和垂向跳動(dòng)。根據(jù)整車參數(shù)計(jì)算各個(gè)工況下輪胎接地點(diǎn)的受力狀態(tài)。采用Adams軟件[14-16]建立穩(wěn)定桿、副車架和擺臂柔性體,并根據(jù)前懸架系統(tǒng)各個(gè)連接點(diǎn)的坐標(biāo)信息、整車參數(shù)、性能曲線等搭建前懸架多體動(dòng)力學(xué)剛?cè)狁詈夏P?,如圖5 所示??v向制動(dòng)(1.2g)時(shí)輪胎垂向載荷為12 132.8 N,輪胎縱向載荷為12 132.8 N。橫向轉(zhuǎn)彎(1.2g)時(shí)輪胎垂向載荷為為11 917.2 N,輪胎橫向載荷為11 917.2 N。垂向跳動(dòng)(3.5g)時(shí)輪胎垂向載荷為57 166.7 N。在輪胎接地點(diǎn)分別施加相應(yīng)的載荷,以此進(jìn)行動(dòng)力學(xué)仿真分析,最終得到下控制臂前點(diǎn)、后點(diǎn)和外點(diǎn)分別在縱向制動(dòng)、橫向轉(zhuǎn)彎和垂向跳動(dòng)極限工況下的力和扭矩。

2.2 強(qiáng)度分析結(jié)果

基于下控制臂有限元網(wǎng)格模型加載提取的力和扭矩,采用慣性釋放方法對(duì)其進(jìn)行強(qiáng)度性能仿真,以此獲取下控制臂在各個(gè)工況下的應(yīng)力水平和分布。

圖6 為下控制臂縱向制動(dòng)工況的應(yīng)力云圖。通過(guò)觀察圖6 可知,下控制臂的最大應(yīng)力為260.8 MPa,位于下控制臂外點(diǎn)附近,這是由于車輛在制動(dòng)時(shí)整車軸荷往前端轉(zhuǎn)移,導(dǎo)致其縱向受力比較大。

圖7 為下控制臂橫向轉(zhuǎn)彎工況的應(yīng)力云圖。通過(guò)圖7 可以看出,下控制臂的最大應(yīng)力為237.8 MPa,位于下控制臂后點(diǎn)區(qū)域,這是由于車輛在橫向轉(zhuǎn)彎時(shí)整車軸荷往左側(cè)偏移,導(dǎo)致其產(chǎn)生局部應(yīng)力集中。

圖8 為下控制臂垂向跳動(dòng)工況應(yīng)力云圖。通過(guò)觀察圖8 可知,下控制臂的最大應(yīng)力為221.2 MPa,位于下控制臂的內(nèi)側(cè)上端,這是由于車輛在垂向跳動(dòng)時(shí)作用力垂直向下,導(dǎo)致其內(nèi)側(cè)發(fā)生應(yīng)力集中現(xiàn)象。

綜上所述,下控制臂在3 種極限工況下的最大應(yīng)力均小于材料屈服值,滿足強(qiáng)度性能要求,符合工程設(shè)計(jì)標(biāo)準(zhǔn)。

3 優(yōu)化設(shè)計(jì)

3.1 優(yōu)化分析方法

通過(guò)下控制臂的模態(tài)性能和強(qiáng)度性能分析可知,其結(jié)構(gòu)特性還有一定的冗余,可以進(jìn)行適當(dāng)?shù)臏p重,以下控制臂的整體厚度值(初始值為8.0 mm)作為設(shè)計(jì)變量,將下擺臂的重量最輕作為目標(biāo)函數(shù),將其強(qiáng)度最大應(yīng)力小于300 MPa和第一階模態(tài)大于110 Hz 作為約束條件。第二代非劣排序遺傳算法探索性強(qiáng),計(jì)算效率高,能夠快速獲取最優(yōu)解。因此基于Isight 軟件并且調(diào)用第二代非劣排序遺傳算法[17-18]對(duì)下控制臂的厚度值進(jìn)行自動(dòng)迭代尋優(yōu),最終得到下控制臂的整體厚度最優(yōu)為7.2 mm。

3.2 優(yōu)化分析結(jié)果

圖9 為優(yōu)化之后下控制臂第一階振型。由圖9可知,優(yōu)化之后下控制臂第一階模態(tài)頻率為110.5Hz,符合動(dòng)態(tài)性能要求,其振型同樣表征為彎曲,后點(diǎn)振幅較大,最大變形為54.1 mm。與此同時(shí),優(yōu)化之后下控制臂第二、三階模態(tài)頻率分別為149.5 Hz和184.6 Hz。

圖10 為優(yōu)化之后下控制臂縱向制動(dòng)工況的應(yīng)力云圖。通過(guò)觀察圖10 可知,優(yōu)化之后,下控制臂的最大應(yīng)力達(dá)到了294.5 MPa,低于材料屈服值,符合強(qiáng)度特性要求。表1 為優(yōu)化前后參數(shù)對(duì)比,由表1 可知,下控制臂的重量降低至1.4 kg,達(dá)到了12.5%的輕量化,優(yōu)化效果比較明顯。優(yōu)化之后下控制臂的前三階模態(tài)頻率有所降低,強(qiáng)度應(yīng)力有所增大,但仍然能夠滿足設(shè)計(jì)要求。

4 結(jié)論

1)采用有限元方法并且基于Hypermesh 軟件建立下控制臂網(wǎng)格模型,采用Nastran 軟件對(duì)其進(jìn)行自由模態(tài)仿真,其前三階固有頻率分別為115.3,159.0 Hz和198.7 Hz,符合振動(dòng)特性要求。

2)基于前懸架多體動(dòng)力學(xué)模型提取下擺臂外聯(lián)點(diǎn)的力和扭矩,對(duì)其進(jìn)行極限強(qiáng)度性能仿真,得到其在縱向制動(dòng)、橫向轉(zhuǎn)彎和垂向跳動(dòng)工況時(shí)的最大應(yīng)力分別為260.8,237.8 MPa和221.2 MPa,均小于材料屈服值,滿足強(qiáng)度設(shè)計(jì)要求。

3)采用第二代非劣排序遺傳算法對(duì)下控制臂的厚度值進(jìn)行迭代優(yōu)化,最終得到下控制臂的整體厚度最優(yōu)為7.2 mm。并且下控制臂的重量減輕12.5%,在兼顧滿足結(jié)構(gòu)性能的前提下達(dá)到了輕量化的目的。

參考文獻(xiàn):

[1] 謝穎,廉俊盛,陳雷.汽車后懸架下控制臂的有限元分 析與優(yōu)化[J].內(nèi)燃機(jī)與配件,2020(3): 31-33.

XIE Y, LIAN J S, CHEN L. Finite element analysis and optimization of the lower control arm of automotive rear suspension[J]. Internal Combustion Engine amp; Parts, 2020(3): 31-33.

[2] 梁國(guó)棟,程樂(lè).某前懸架下擺臂疲勞失效分析與優(yōu)化[ J].機(jī)械設(shè)計(jì)與制造,2022(3): 292-295.

LIANG G D, CHENG L. Fatigue failure analysis and optimization of a front suspension lower arm[J]. Machinery Design amp; Manufacture, 2020(3): 292-295.

[3] 徐峰.麥弗遜懸架下擺臂拓?fù)鋬?yōu)化設(shè)計(jì)[J].內(nèi)燃機(jī)與配 件,2023(24): 59-61.

XU F. Topological optimization design of mcpherson suspension hem arm[J]. Internal Combustion Engine amp; Parts, 2023(24): 59-61.

[4] SONG X G, JUNG J H, SON H J, et al. Metamodel-based optimization of a control arm considering strength and durability performance[J]. Computers amp; Mathematics with Applications, 2010, 60(4): 976-980.

[5] 鄒喜紅,茍林林,傅雷,等.全域車懸架下控制臂有限元 分析及輕量化設(shè)計(jì)[J].公路交通科技,2023, 40(12): 227-235.

ZOU X H, GOU L L, FU L, et al. Finite element analysis and lightweight design for lower suspension control arm of all-terrain vehicle[J]. Journal of Highway and Transportation Research and Development, 2023, 40(12): 227-235.

[6] 胡軍,朱文豐,許凱樂(lè),等.基于模態(tài)分析的機(jī)床部件夾 層結(jié)構(gòu)阻尼參數(shù)的特性研究[J].機(jī)械強(qiáng)度,2020, 42 (3): 580-584.

HU J, ZHU W F. XU K L, et al. Research on characteristics of damping parameters of sandwich structure of machine tool components based on modal analysis[J]. Journal of Mechanical Strength, 2020, 42(3): 580-584.

[7] 黃海彬.基于有限元技術(shù)的儀表板綜合特性評(píng)估及結(jié) 構(gòu)輕量化[J].工程塑料應(yīng)用,2021, 49(8): 123-128.

HUANG H B. Comprehensive characteristic evaluation and structure lightweight of dashboard based on finite element technology[J]. Engineering Plastics Application, 2021,49(8): 123-128.

[8] 丁杰.地鐵車輛輔助變流器的噪聲仿真及測(cè)試[J].華東 交通大學(xué)學(xué)報(bào),2022, 39(1): 99-107.

DING J. Noise simulation and test of auxiliary converter for metro vehicle[J]. Journal of East China Jiaotong University, 2022, 39(1): 99-107.

[9] 丁杰,蔣莉,徐慶怡.基于HyperMesh的動(dòng)車組車體結(jié) 構(gòu)強(qiáng)度仿真[J].機(jī)電工程技術(shù),2023, 52(7): 268-271.

DING J, JIANG L, XU Q Y. Simulation of structural strength for emu body based on hypermesh[J]. Mechanical amp; Electrical Engineering Technology, 2023, 52(7): 268-271.

[10] 劉瓊,楊紀(jì)遠(yuǎn),張幕,等.基于Hypermesh的空氣彈簧成 型鼓模態(tài)分析[J].福建理工大學(xué)學(xué)報(bào),2023, 21(6): 544-550.

LIU Q, YANG J Y, ZHANG L, et al. Modal analysis of air spring forming drum based on Hypermesh[J]. Journal of FuJian University of Technology, 2023, 21(6): 544-550.

[11] 李翔翼.某輕卡蓄電池支架性能分析與輕量化設(shè)計(jì) [J].機(jī)械工程與自動(dòng)化,2023(6): 98-100.

LI X Y. Performance analysis and lightweight design of light truck battery bracket[J]. Mechanical Engineering amp; Automation, 2023(6): 98-100.

[12] 姜麗紅,吳慶捷.某前副車架性能分析與輕量化設(shè)計(jì) [J].機(jī)械強(qiáng)度,2020, 42(6): 1503-1508.

JIANG L H, WU Q J. Performance analysis and lightweight design of a front subframe[J]. Journal of Mechanical Strength, 2020, 42(6): 1503-1508.

[13] 陳佳,孫麗娟.基于NASTRAN的某商用輕卡后車門力 學(xué)性能仿真研究[J].汽車實(shí)用技術(shù),2021, 46(10): 8788.

CHEN J, SUN L J. Simulation study on the mechanical performance of a commercial light truck rear door[J]. Automobile AppliedTechnology, 2021, 46(10): 87-88.

[14] 劉克銘,袁鑫,戴羽森.某SUV前副車架強(qiáng)度與疲勞性 能仿真與試驗(yàn)研究[J].機(jī)械強(qiáng)度,2023, 45(4): 838-844.

LIU K M, YUAN X, DAI R S, et al. Simulation and experimental study on strength and fatigue performance of a suv front subframe[J]. Journal of Mechanical Strength, 2023,45(4): 838-844.

[15] 余志武,黎燕霞,徐磊,等.曲線導(dǎo)軌上高速磁浮列車運(yùn) 動(dòng)穩(wěn)定性[J].華東交通大學(xué)學(xué)報(bào),2023, 40(1): 25-33.

YU Z W, LI Y X, XU L, et al. Motion stability of highspeed maglev train on curved track[J]. Journal of East China Jiaotong University, 2023, 40(1): 25-33.

[16] 黃美婷.基于ADAMS的四輪獨(dú)立驅(qū)動(dòng)電動(dòng)汽車懸架 分析與優(yōu)化[J].河南工程學(xué)院學(xué)報(bào)(自然科學(xué)版), 2024, 36(1): 50-54.

HUANG M T. Suspension analysis and optimization of the four- wheel independent drive electric vehicle based on ADAMS[J]. Journal of Henan University of Engi-neering(Natural Science Edition), 2024, 36(1): 50-54.

[17] 石越卿,秦睿賢,陳秉智.列車防撞柱結(jié)構(gòu)耐撞性設(shè)計(jì) 與優(yōu)化[J].機(jī)械強(qiáng)度,2022, 44(5): 1134-1140.

SHI Y Q, QIN R X, CHEN B Z. Crash-worthiness design and optimization for collision post of train[J]. Journal of Mechanical Strength, 2022, 44(5): 1134-1140.

[18] 王康,曹永晟,賀啟才,等.基于多目標(biāo)優(yōu)化的白車身結(jié) 構(gòu)輕量化設(shè)計(jì)[J].時(shí)代汽車,2023(22): 151-154.

WANG K, CAO Y S, HE Q C, et al. Lightweight design of body-in-white structure based on multi-objective opti-mization[J]. Auto Time, 2023(22): 151-154.

通信作者:游永忠(1969—),男,講師,研究方向?yàn)榻煌ㄟ\(yùn)輸設(shè)備。E-mail:1493779073@qq.com。

(責(zé)任編輯:龔凱,李根)

基金項(xiàng)目:江西省教育廳科學(xué)技術(shù)研究項(xiàng)目(GJJ200647)

猜你喜歡
模態(tài)優(yōu)化
超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
車輛CAE分析中自由模態(tài)和約束模態(tài)的應(yīng)用與對(duì)比
國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
高速顫振模型設(shè)計(jì)中顫振主要模態(tài)的判斷
基于低碳物流的公路運(yùn)輸優(yōu)化
基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
主站蜘蛛池模板: 中文字幕在线播放不卡| 四虎影视国产精品| aaa国产一级毛片| 午夜毛片免费观看视频 | 色婷婷天天综合在线| 国产美女在线观看| 91无码视频在线观看| 国产精品自拍合集| 国产免费久久精品99re不卡| 欧美视频在线不卡| 一本色道久久88| 人与鲁专区| 四虎成人精品| 欧美激情第一欧美在线| 99视频在线免费| 无码高潮喷水在线观看| 国产地址二永久伊甸园| 国产欧美专区在线观看| 精品无码一区二区三区电影| 国产欧美日韩在线一区| 全裸无码专区| 日韩精品欧美国产在线| av在线人妻熟妇| 欧美成人A视频| 极品私人尤物在线精品首页 | 九色在线观看视频| 五月激情综合网| 欧美在线天堂| 91偷拍一区| 国产导航在线| 日韩天堂视频| 无码国内精品人妻少妇蜜桃视频| 久久中文字幕2021精品| 国产丝袜无码精品| 尤物午夜福利视频| 2019年国产精品自拍不卡| 精品91视频| 欧美亚洲激情| 亚洲一级毛片在线播放| 亚洲欧美日韩另类在线一| 91精品久久久久久无码人妻| 亚洲成人高清在线观看| 精品视频一区二区三区在线播| 99热国产这里只有精品9九 | 日本一区二区三区精品视频| 国产丝袜一区二区三区视频免下载| 1024国产在线| 在线视频一区二区三区不卡| 在线观看精品自拍视频| 亚洲欧美h| 国产哺乳奶水91在线播放| 欧美色丁香| 精品成人免费自拍视频| 久久精品人妻中文系列| 91香蕉视频下载网站| 91在线激情在线观看| 精品久久国产综合精麻豆| 亚洲精品老司机| 亚洲丝袜中文字幕| 午夜精品福利影院| 免费观看精品视频999| 国产欧美日韩视频怡春院| 在线观看免费国产| 久久福利片| 国产永久免费视频m3u8| 免费三A级毛片视频| 97精品久久久大香线焦| 国产在线自乱拍播放| 九九热精品在线视频| 亚洲AV免费一区二区三区| 欧美精品1区2区| 日韩在线中文| 色综合天天综合中文网| 欧美激情视频一区二区三区免费| 色一情一乱一伦一区二区三区小说| 最新国语自产精品视频在| 国产成人8x视频一区二区| 亚洲三级影院| 亚洲最大福利视频网| 国产激情在线视频| 亚洲欧美日韩另类| 国产激情在线视频|