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

伺服機械壓力機機身結構優化設計分析

2022-07-26 09:35:32王俊劉祥龐秋胡志力
精密成形工程 2022年7期
關鍵詞:有限元機械優化

王俊,劉祥,龐秋,胡志力

應用技術

伺服機械壓力機機身結構優化設計分析

王俊1a,1b,劉祥2,龐秋3,胡志力1a,1b

(1.武漢理工大學 a.現代汽車零部件技術湖北省重點實驗室;b.材料綠色精密成形技術與裝備湖北省工程中心,武漢 430070;2.東風(武漢)實業有限公司,武漢 430040;3.武漢科技大學 機械自動化學院,武漢 430081)

針對當前伺服機械壓力機機身質量與結構分配不合理、剛度不足的問題,以6 000 kN閉式伺服機械壓力機機身為研究對象,對該閉式伺服機械壓力機機身進行重新設計,以實現提高機身剛度和輕量化的目的。首先進行機身靜力學分析,確定拓撲優化空間,獲得壓力機機身優化分析邊界條件;然后采用變密度法對壓力機機身進行拓撲優化分析,選擇合適的密度閾值,獲得機身的拓撲優化結構。為了便于加工制造,減少制造成本,基于該優化結構并考慮可制造性重新設計機身結構。最后,通過有限元仿真分析和機身剛度測試試驗,對比優化前后機身的剛度。優化后的伺服機械壓力機機身質量減輕了10.9%,一階模態頻率提高了3.74%,機身剛度提高了約28%。通過對伺服機械壓力機機身結構進行優化設計,解決了高剛度機身設計的工程問題,為伺服機械壓力機生產制造提供了一定的理論和技術支撐。

伺服機械壓力機;機身;結構優化設計;OptiStruct

隨著機械制造業的快速發展,作為金屬成形設備的鍛壓設備也得到了越來越廣泛的應用,其中伺服機械壓力機是目前金屬成形加工領域中最先進的設備。伺服機械壓力機不僅機身整體結構復雜、使用的加工材料多,還具有高昂的制造成本,因此機身部件的整體質量和材料分布對整個伺服機械壓力機的工作性能有重要的影響[1]。伺服機械壓力機在工作時,產生的大部分力會由機身承載,因此機身會產生一定的變形,機身的變形不僅會影響零件的加工質量和生產效率,還會影響模具的使用壽命和整個壓力機的使用安全性[2]。對壓力機機身進行優化設計,不僅可以降低制造成本,實現輕量化,還可以提高機身的剛度和安全性能。為了保證使用性能,采用傳統方法設計的壓力機機身往往過于笨重,存在結構設計不合理和大量的材料浪費等問題[3]。

伺服機械壓力機機身整體性能與其結構設計和機身材料分布有著很大的關系[4]。傳統機身設計方法設計的機身存在結構不合理和質量過大等問題,因此需要在保證機身使用性能的前提下對機身進行結構優化設計,提高機身的整體性能[5]。為了對機身進行結構優化、提高機身的整體性能,許多高校學者對機身的優化設計進行了大量研究。劉運璽[6]通過尺寸優化的方法,以壓力機機身質量最小為目標函數,使機身質量減少了5.5%。錢抗抗[7]通過優化改變機身的板厚尺寸,以壓力機機身體積最小為目標函數,對機身進行尺寸優化,使機身的質量大大降低。上述學者對機身的優化設計主要還停留在尺寸設計方面,這雖然可以在一定程度上實現提高機身剛度和輕量化的目的,但機身整體的減重效果受到限制。

拓撲優化是通過給定模型在實際工況下的外部載荷和邊界條件,在一定的設計區域內尋找一種最優化的模型結構,以此來達到所設定的最大變形量等約束條件,使模型結構能夠達到一種或多種最優性能的設計方法[8-10]。王蘇號等[11]運用改進的優化算法,以位移、應力和頻率為約束條件對壓力機機身進行拓撲優化,實現了機身質量的減輕。饒柳生等[12]采用多目標拓撲優化,根據優化后的機身模型密度分布對機身進行優化設計,達到了機身動態和靜態特性提高的目的。

文中通過2種方案對壓力機機身進行拓撲優化,以機身有限元單元的密度為設計變量,通過采用不同的約束條件,以機身的質量最小和剛度最大為目標函數,得到了機身材料的最優分布,實現了提高機身剛度和輕量化的目的。

1 壓力機機身有限元建模及分析

1.1 原壓力機機身結構

文中的研究對象是某公司生產的6 000 kN伺服機械壓力機機身。伺服機械壓力機在工作時,整個機身承受來自傳動系統和工作臺的壓力,整個壓力機采用八面導向機構,具有良好的抗偏載性能和優異的導向性能。該伺服機械壓力機機身由底座、工作臺、立柱、橫梁等組成,主要由鋼板焊接而成,整個機身通過螺栓與地面連接。機身的三維模型是后續進行網格劃分的關鍵,網格劃分的好壞對后續的分析極為重要,良好的網格劃分可以提高運算的速度和精度,文中采用三維建模軟件SolidWorks,以機身實際結構尺寸為參考,建立機身的三維模型。為了提高網格劃分的質量、降低計算工作量,需要對機身進行一些合理的簡化,將底座、工作臺、立柱、橫梁視為一個整體,忽略一些機身倒角和小孔等微小細節以簡化形狀[13-14]。圖1為建立的原壓力機機身結構。

為保證伺服機械壓力機的工作精度,其機身應該滿足基本的工作要求,例如:具有良好的結構強度,機身工作時最大工作應力應在許用應力范圍內;具有符合使用要求的整體剛度,通常用額定載荷下的柔度大小來衡量,即柔度越小,整體剛度越好;具有較好的局部剛度,即機身受力部位的局部變形要在合理范圍內;具有良好的動態特性,機身的一階固有頻率應遠離伺服機械壓力機其他部件的頻率,以避免產生共振,一般要求低階固有頻率越高越好[15-16]。

圖1 原壓力機機身結構

1.2 壓力機機身有限元模型建立

首先通過三維建模軟件SolidWorks建立機身的三維模型,導出為通用的.IGS格式[17-18],然后將建立的簡化機身三維模型導入到HyperMesh中,建立機身結構優化分析的有限元模型,最后對建立的有限元分析模型進行分析優化。其中有限元建模主要包含以下幾個方面:(1)節點和單元,采用實體單元建模,共建立189 256個節點和826 882個四面體單元,將軸孔受力面的點耦合在一個節點,方便后面施加載荷;(2)載荷施加,在2個軸孔耦合點處分別施加集中載荷,大小都為3 000 kN,對工作臺面施加均布力,每個節點的力為2.3 kN;(3)材料屬性,材料采用Q235鋼,其屈服強度為235 MPa;(4)約束定義,實際工作時機身通過螺栓固定在地面上,因此要通過對機身模型底座設置全約束來限制機身底座各個方向的自由度。

建立有限元分析模型時還要注意網格的質量和數量,劃分網格的質量和數量對分析的結果非常重要[19-20]。在劃分網格時,要充分考慮機身的關鍵部位,應保證關鍵部位的網格劃分質量良好。圖2為建立的壓力機機身有限元分析模型。

圖2 壓力機機身有限元分析模型

1.3 壓力機機身有限元結果分析

通過求解器OptiStruct求解計算在HyperMesh中建立的有限元分析模型,并在HyperView中查看計算的結果,經過后處理得到機身的應變云圖和應力云圖,分別如圖3和圖4所示。提取壓力機機身有限元分析的最大應力max、體積和各階固有頻率。

計算結果顯示,機身最大變形量出現在工作臺處和上部橫梁處,為0.33 mm,機身的底座未出現明顯變形,最大應力值為119.5 MPa。通過分析結果可知,該機身能夠滿足使用要求,但整體質量過大,橫梁和立柱部分仍有較大的優化空間。因此,可以對機身進行結構優化,實現提高壓力機機身剛度和輕量化的目的。

圖3 壓力機機身應變分布

圖4 壓力機機身應力分布

2 壓力機機身結構拓撲優化

2.1 壓力機機身拓撲優化三維模型

根據現有壓力機機身的結構和尺寸大小,在三維建模軟件中對壓力機機身進行填充設計,建立壓力機機身拓撲優化的三維模型[21]。這一過程主要是對原壓力機機身底座區域和上橫梁區域進行填充,因為這些地方存在較大的設計潛力,通過填充設計建立機身拓撲優化的三維模型,圖5為建立的壓力機機身填充三維模型。

圖5 壓力機機身填充三維模型

2.2 壓力機機身拓撲優化有限元分析模型

對填充后的壓力機機身模型進行網格劃分,建立有限元分析模型,其約束和載荷的定義與1.2節壓力機機身分析模型中的保持一致,然后進行拓撲優化分析。

文中主要采用了2種拓撲優化方案對機身進行結構優化。第1種方案的優化目標為機身模型質量最小,約束條件為機身的最大變形量小于0.3 mm、一階頻率不小于29.12 Hz,設計變量為機身模型的單元密度。第2種方案的優化目標為機身模型剛度最大,約束條件為機身的最大變形量小于0.3 mm、一階頻率不小于29.12 Hz,設計變量為機身模型的單元密度。

2.3 壓力機機身拓撲優化結果分析

對壓力機機身有限元模型進行優化求解分析,通過HyperView查看優化結果,并將有限元網格單元的相對密度閾值設置為0.35,只顯示相對密度閾值大于0.35的單元。圖6為方案1的優化結果,圖7為方案2的優化結果。

圖6 方案1優化結果

圖7 方案2優化結果

對比優化后的結果發現,2種優化結果所去除材料的區域基本相似,機身頂端橫梁和立柱外側有部分單元消失,下端橫梁部分有大面積的單元消失。根據拓撲優化后的結果并考慮可制造性,對壓力機機身進行重新設計,圖8為重新設計的壓力機機身結構。

圖8 重新設計的壓力機機身結構

2.4 壓力機機身優化前后對比分析

為了驗證重新設計后的壓力機機身的結構剛度和強度,對重新設計的機身進行靜態性能和模態分析。圖9為優化后壓力機機身的應變分布圖,圖10為優化后壓力機機身的應力分布圖。

圖9 優化后壓力機機身應變分布

圖10 優化后壓力機機身應力分布

優化后壓力機機身的最大變形量為0.322 mm,壓力機機身的最大應力為126.1 MPa,一階模態頻率為30.21 Hz,符合使用要求,然后提取壓力機機身的性能參數。表1為原壓力機機身和優化后壓力機機身的性能參數對比,1—6分別為壓力機機身的前六階模態頻率。

通過對比表1的數據,可以得出:(1)優化后壓力機機身的質量減小了10.9%,實現了機身輕量化的目標;(2)優化后壓力機機身最大應力為126.10 MPa,仍小于材料的許用應力,滿足使用要求;(3)優化后壓力機機身的一階固有頻率提高了3.74%,機身的動態性能得到了改善。

表1 壓力機機身性能參數對比

Tab.1 Comparison of press frame performance parameters

3 壓力機機身優化試驗驗證分析

對2個壓力機機身模型進行同比例縮放,然后設置相同的3D打印參數打印出壓力機機身的實體塑料模型,圖11為優化后的壓力機機身實體塑料模型,圖12為原壓力機機身實體塑料模型。通過電子萬能試驗機對2個壓力機機身實體塑料模型進行壓縮試驗驗證,分別對2個模型施加從小到大的力,電子萬能試驗機的壓縮速度設置為1 mm/min,在力達到3 kN時,模型出現壓潰現象,停止壓縮,分別記錄2個壓力機機身實體塑料模型變形量和所受力的數據,并繪制成曲線進行比較。圖13為2個壓力機機身實體塑料模型所受壓力與變形量的曲線。

通過對比曲線可知,在2個壓力機機身實體塑料模型所受力相同的情況下,優化后的壓力機機身的變形量要略小于原壓力機機身的變形量,當力施加到3 kN時,原壓力機機身最大變形量為1.37 mm,優化后壓力機機身最大變形量為1.07 mm,機身剛度提高了28%。

圖11 優化后壓力機機身實體塑料模型

圖12 原壓力機機身實體塑料模型

圖13 壓力機機身所受力與變形量曲線

4 結論

1)基于不同的優化目標函數對伺服機械壓力機機身進行拓撲優化分析,并以分析結果為參考對伺服機械壓力機機身進行重新設計,在伺服機械壓力機機身性能滿足使用要求的前提下,機身質量減小了10.9%,達到了機身輕量化和節約成本的目的。

2)通過電子萬能試驗機對優化前后的壓力機機身進行壓縮承載試驗驗證,發現優化后的壓力機機身剛度提高了28%。

3)重新設計的壓力機機身材料總體呈現下重上輕的分布趨勢,機身的整體重心下移,有利于提高機身的動態性能。

[1] 陳聰. 高速壓力機機構優化設計及輕量化研究[D]. 大連: 大連理工大學, 2021: 9-13.

CHEN Cong. Optimization Design of Mechanism and Frame Structure of High-Speed Press[D]. Dalian: Dalian University of Technology, 2021: 9-13.

[2] 白玉冰, 侯佳欣, 張志國, 等. 基于有限元計算的冷鍛壓力機組合機身的結構優化[J]. 鍛壓技術, 2020, 45(6): 130-136.

BAI Yu-bing, HOU Jia-xin, ZHANG Zhi-guo, et al. Structural Optimization on Combination Body of Cold Forging Press Based on Finite Element Calculation[J]. Forging & Stamping Technology, 2020, 45(6): 130-136.

[3] 魏鳳凱, 黃慧, 單本軍, 等. 雙點機械壓力機機身有限元分析及優化[J]. 鍛壓裝備與制造技術, 2020, 55(3): 26-30.

WEI Feng-kai, HUANG Hui, SHAN Ben-jun, et al. Finite Element Analysis and Optimization of Frame for the Double-Point Mechanical Press[J]. China Metalforming Equipment & Manufacturing Technology, 2020, 55(3): 26-30.

[4] YI G D, WU J D, LI R S. Static and Dynamic Deformation Analysis and Optimization of a Fine Blanking Press Frame[C]// International Conference on Mechanical Design. Springer, Singapore, 2017.

[5] KORTA J, RANIOLO R, DANTI M, et al. Multi-Objective Optimization of a Car Body Structure [J]. SAE International Journal of Passenger Cars-Mechanical Systems, 2012, 5(3): 1143-1152.

[6] 劉運璽. HD-200高速壓力機工作機構的動力學仿真及其機身結構優化[D]. 濟南: 山東大學, 2013: 7-8.

LIU Yun-xi. Dynamics Simulation of the Operating Mechanism and Structure Optimization of the Frame of the HD-200High-Speed Press[D]. Jinan: Shandong University, 2013: 7-8.

[7] 錢抗抗. JG21-250A型壓力機機身的有限元分析及其優化設計[D]. 青島: 山東科技大學, 2008: 9-14.

QIAN Kang-kang. FEM Analysis and Optimum Design for Frame of JG21-250A Press[D]. Qingdao: Shandong University of Science and Technology, 2008: 9-14.

[8] 何健, 何猛, 夏涼, 等. 基于雙向漸進結構優化法的柔性機構設計[J]. 機械工程學報, 2021, 57(19): 39-47.

HE Jian, HE Meng, XIA Liang, et al. Design of Compliant Actuation Mechanisms by Evolutionary Structural Optimization Method[J]. Journal of Mechanical Engineering, 2021, 57(19): 39-47.

[9] XU X, WANG L, LIU G, et al. Optimal Design of Dual Derricks[J]. International Journal of Modelling Identification & Control, 2012, 17(1): 61-67.

[10] 黃志超, 程雯玉, 占金青, 等. 基于目標驅動優化技術的自沖鉚接機C型框架結構優化[J]. 鍛壓技術, 2013, 38(5): 104-109.

HUANG Zhi-chao, CHENG Wen-yu, ZHAN Jin-qing, et al. Structure Optimization of C-Frame Self-Piercing Riveting Machine Using Goal Driven Optimization[J]. Forging & Stamping Technology, 2013, 38(5): 104-109.

[11] 王蘇號, 賈方, 王興松. 基于Optistruct的伺服壓力機機身拓撲優化[J]. 鍛壓裝備與制造技術, 2007, 42(6): 34-36.

WANG Su-hao, JIA Fang, WANG Xing-song. Topology Optimization for Servo Press Body Based on Optistruct[J]. China Metal Forming Equipment & Manufacturing Technology, 2007, 42(6): 34-36.

[12] 饒柳生, 侯亮, 潘勇軍. 基于拓撲優化的機床立柱筋板改進[J]. 機械設計與研究, 2010, 26(1): 87-92.

RAO Liu-sheng, HOU Liang, PAN Yong-jun. Improvement of the Rib Plate in Machine Tools Column Based on Topology Optimization[J]. Machine Design & Research, 2010, 26(1): 87-92.

[13] 徐晶, 劉基岡, 張懿驊. 基于Hypermesh的液壓機機身結構優化設計[J]. 鍛壓技術, 2015, 40(11): 48-52.

XU Jing, LIU Ji-gang, ZHANG Yi-hua. Structure Optimization Design of Hydraulic Press Body Based on Hypermesh[J]. Forging & Stamping Technology, 2015, 40(11): 48-52.

[14] 朱從武, 張劍平, 王昊誠, 等. 開式壓力機機身優化分析[J]. 鍛壓裝備與制造技術, 2019, 54(2): 13-15.

ZHU Cong-wu, ZHANG Jian-ping, WANG Hao-cheng, et al. Optimized Analysis of Open Press Frame[J]. China Metalforming Equipment & Manufacturing Technology, 2019, 54(2): 13-15.

[15] ZHAO Xin-hao, LIU Yan-xiong, HUA Lin, et al. Finite Element Analysis and Topology Optimization of a 12000KN Fine Blanking Press Frame[J]. Structural and Multidisciplinary Optimization, 2016, 54(2): 375-389.

[16] 王中原, 霍光, 李宣, 等. 熱模鍛壓力機機身降重效率有限元分析[J]. 鍛壓裝備與制造技術, 2020, 55(3): 59-62.

WANG Zhong-yuan, HUO Guang, LI Xuan, et al. Finite Element Analysis on Weight Reduction Efficiency of Hot Die Forging Press Fuselage[J]. China Metalforming Equipment & Manufacturing Technology, 2020, 55(3): 59-62.

[17] YA B, AV V. Computer Aided Analysis of Screw Press Dynamics[J]. International Journal of Machine Tools and Manufacture, 1987, 27(2): 143-154.

[18] 韋源源, 劉子越, 李柱, 等. 基于ANSYS的大型壓力機機身結構分析與改進設計[J]. 機械工程與自動化, 2020(2): 45-46.

WEI Yuan-yuan, LIU Zi-yue, LI Zhu, et al. Structural Analysis and Improvement Design of Large Press Body Based on ANSYS[J]. Mechanical Engineering & Automation, 2020(2): 45-46.

[19] 魏鳳凱, 王四森, 李樂鑫, 等. SP-160開式壓力機機身的有限元分析及優化[J]. 鍛壓裝備與制造技術, 2016, 51(1): 17-20.

WEI Feng-kai, WANG Si-sen, LI Le-xin, et al. Finite Element Analysis and Optimization of Frame for SP-160 Open Press[J]. China Metal Forming Equipment & Manufacturing Technology, 2016, 51(1): 17-20.

[20] 姚菁琳, 張新國, 吳書森, 等. 基于有限元的大噸位壓力機立柱結構改進[J]. 鍛壓裝備與制造技術, 2013, 48(6): 54-57.

YAO Jing-lin, ZHANG Xin-guo, WU Shu-sen, et al. Improvement of Column Structure for Large-Tonnage Press Based on Finite Element Method[J]. China Metal Forming Equipment & Manufacturing Technology, 2013, 48(6): 54-57.

[21] 隋曉東, 戴挺, 高鋒. 壓力機機身振動特性的模態分析[J]. 噪聲與振動控制, 2010, 30(1): 10-14.

SUI Xiao-dong, DAI Ting, GAO Feng. Modal Analysis of Compressor Frame's Vibration[J]. Noise and Vibration Control, 2010, 30(1): 10-14.

Optimization Design and Analysis of Servo Mechanical Press Frame Structure

WANG Jun1a,1b, LIU Xiang2, PANG Qiu3, HU Zhi-li1a,1b

(1. a. Hubei Key Laboratory of Advanced Technology for Automotive Components; b. Hubei Engineering Research Center for Green & Precision Material Forming, Wuhan University of Technology, Wuhan 430070, China; 2. Dongfeng (Wuhan) Industrial Co., Ltd., Wuhan 430040, China; 3. School of Machinery and Automation, Wuhan University of Science and Technology, Wuhan 430081, China)

The work aims to take the 6 000 kN closed servo mechanical press frame as the research object, and redesign this closed servo mechanical press frame to solve the problems of unreasonable frame mass and structure distribution and insufficient stiffness in current servo mechanical press, thus improving the stiffness and lightening the weight of frame. Firstly, the static analysis of the frame was carried out and the space for topology optimization of the servo mechanical press frame was determined to obtain the boundary conditions for the optimization analysis of the press frame. Then, the variable density method was used to analyze the topology optimization of the press frame, and the appropriate density threshold was selected to obtain the topology optimization structure of the press frame. In order to facilitate manufacturing and reduce manufacturing costs, the frame structure was redesigned based on the optimized structure and manufacturability. Finally, the stiffness of the frame before and after the optimization was compared through the finite element simulation analysis and the frame stiffness test. From the results, the weight of the optimized servo mechanical press frame was reduced by 10.9%, the first-order modal frequency was increased by 3.74%, and the frame stiffness was increased by about 28%. The optimization design of servo mechanical press frame structure solves the engineering problem in the design of high stiffness frame and provides a certain theoretical and technical support for the production and manufacturing of servo mechanical press.

servo mechanical press;frame;structural optimization design; OptiStruct

10.3969/j.issn.1674-6457.2022.07.019

TG315.5

A

1674-6457(2022)07-0136-07

2022–01–15

國家重點研發計劃(2019YFB1704500);國家自然科學基金(52075400);湖北省重點研發計劃(2020BAB140)

王俊(1997—),男,碩士生,主要研究方向為輕量化技術。

龐秋(1979—),女,博士,副教授,主要研究方向為汽車輕量化設計制造技術。

責任編輯:蔣紅晨

猜你喜歡
有限元機械優化
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
調試機械臂
當代工人(2020年8期)2020-05-25 09:07:38
簡單機械
機械班長
按摩機械臂
磨削淬硬殘余應力的有限元分析
基于SolidWorks的吸嘴支撐臂有限元分析
主站蜘蛛池模板: 无码电影在线观看| 日本人妻丰满熟妇区| 91视频首页| 原味小视频在线www国产| 欧美三级视频在线播放| 日韩乱码免费一区二区三区| 不卡午夜视频| 日韩在线播放欧美字幕| 久久香蕉欧美精品| 日本手机在线视频| 91小视频在线| 久久国产成人精品国产成人亚洲| 午夜国产精品视频黄| 中文一级毛片| 四虎永久在线| 91久久偷偷做嫩草影院精品| 亚洲色图欧美在线| 国产青榴视频| 亚洲美女一区二区三区| a毛片在线免费观看| 久久www视频| 亚洲视频三级| 国产玖玖视频| 天天色综网| 无码丝袜人妻| 欧美一区国产| 亚洲欧美成人综合| 国产凹凸视频在线观看| 国产精品亚洲αv天堂无码| 国产欧美日韩专区发布| 亚洲一区二区无码视频| yjizz国产在线视频网| 欧美啪啪一区| 一区二区午夜| 国产真实乱子伦视频播放| 国产成人av一区二区三区| 伊人国产无码高清视频| 亚洲国产精品日韩av专区| 国产精品亚洲一区二区三区z| 亚洲一区波多野结衣二区三区| 国产sm重味一区二区三区| 女人18毛片水真多国产| 亚洲V日韩V无码一区二区| 国产精品美人久久久久久AV| 久久国产精品国产自线拍| 中文字幕佐山爱一区二区免费| 九九久久99精品| 国产免费久久精品99re不卡| 国产午夜一级毛片| 国产网站免费| 欧美精品v| 久久一本精品久久久ー99| 丁香婷婷久久| 日本一区二区不卡视频| 国产福利大秀91| 国产午夜不卡| 狠狠色狠狠综合久久| 天天色天天操综合网| 国产亚洲视频播放9000| 视频在线观看一区二区| 色综合天天综合中文网| 热热久久狠狠偷偷色男同| 毛片网站在线看| 亚洲婷婷丁香| 看国产毛片| 久久精品无码专区免费| 国产成人精品视频一区视频二区| 一本综合久久| 亚洲人成在线精品| 在线va视频| 女人18毛片久久| 日韩第一页在线| 国产精品私拍99pans大尺度| 蜜臀AV在线播放| 欧美成人综合在线| 一区二区三区国产| 色偷偷综合网| 国产精品亚洲一区二区三区z| 无码国产偷倩在线播放老年人| 91麻豆久久久| 欧美在线国产| 五月天久久综合|