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

        ?

        基于多目標(biāo)遺傳算法的雷載計(jì)算機(jī)隔振系統(tǒng)優(yōu)化設(shè)計(jì)

        2017-05-06 02:16:28張溢文尹韶平王志杰吳培明崔鑫山
        關(guān)鍵詞:模態(tài)計(jì)算機(jī)優(yōu)化

        張溢文, 尹韶平, 王志杰, 郭 君, 吳培明, 崔鑫山

        ?

        基于多目標(biāo)遺傳算法的雷載計(jì)算機(jī)隔振系統(tǒng)優(yōu)化設(shè)計(jì)

        張溢文1,2, 尹韶平1, 王志杰1, 郭 君1, 吳培明1, 崔鑫山1,2

        (1. 中國(guó)船舶重工集團(tuán)公司第705研究所, 陜西西安, 710077; 2. 水下信息與控制重點(diǎn)實(shí)驗(yàn)室, 陜西西安, 710077)

        為解決復(fù)雜力學(xué)環(huán)境下雷載計(jì)算機(jī)的振動(dòng)隔離問(wèn)題, 提出了一種基于UG和ANSYS Workbench的聯(lián)合優(yōu)化設(shè)計(jì)方法。在UG中建立了隔振系統(tǒng)的參數(shù)化幾何模型, 并導(dǎo)入ANSYS Workbench建立有限元模型。計(jì)算了原始設(shè)計(jì)方案的模態(tài)參數(shù),通過(guò)與模態(tài)試驗(yàn)結(jié)果對(duì)比,驗(yàn)證了計(jì)算模型的正確性。隨后依據(jù)隔振設(shè)計(jì)理論確定了優(yōu)化目標(biāo)函數(shù),使用多目標(biāo)遺傳算法對(duì)隔振系統(tǒng)進(jìn)行了優(yōu)化設(shè)計(jì)。通過(guò)對(duì)比可知, 優(yōu)化方案較原始方案的最高階固有頻率降低了12.9%,頻率間隔降低了79.6%, 系統(tǒng)的有效隔振頻率增加了39.1 Hz,達(dá)到了良好的設(shè)計(jì)效果。

        雷載計(jì)算機(jī);優(yōu)化設(shè)計(jì);多目標(biāo)遺傳算法;隔振;模態(tài)理論

        0 引言

        輕型魚雷作為主要反潛武器, 可在多種武器平臺(tái)使用, 如可通過(guò)水面艦艇發(fā)射, 反潛飛機(jī)空投, 亦可通過(guò)助飛火箭進(jìn)行遠(yuǎn)程投送[1]。雷載計(jì)算機(jī)作為輕型魚雷的“大腦”, 負(fù)責(zé)全雷的信息處理及指令收發(fā)工作, 在雷載設(shè)備中占有重要的地位。輕型魚雷多平臺(tái)使用的特點(diǎn)決定了其雷載計(jì)算機(jī)在全生命周期中會(huì)經(jīng)歷復(fù)雜多樣的力學(xué)環(huán)境。因此, 依據(jù)隔振設(shè)計(jì)理論對(duì)雷載計(jì)算機(jī)隔振系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì), 增強(qiáng)雷載計(jì)算機(jī)的環(huán)境適應(yīng)性和可靠性具有重要的工程意義。

        針對(duì)電子設(shè)備振動(dòng)隔離問(wèn)題, 周亞?wèn)|等人[2]提出了一種慣導(dǎo)平臺(tái)隔振器的斜角布置方法, 并通過(guò)正則模態(tài)分析對(duì)比了改進(jìn)方案與傳統(tǒng)方案的模態(tài)特性。Kim Sang-Myeong等人[3]使用分散式速度反饋控制算法對(duì)多路主動(dòng)隔振系統(tǒng)進(jìn)行控制, 并通過(guò)試驗(yàn)驗(yàn)證了其隔振效能。A. Vulcan等人[4]為確保戰(zhàn)斗機(jī)雷達(dá)晶體振蕩器可以在高能量隨機(jī)振動(dòng)和高過(guò)載工況下正常工作, 設(shè)計(jì)了一種懸浮式新型晶體振蕩器??梢钥闯? 前人在振動(dòng)隔離問(wèn)題上的研究主要集中于提出新的隔振方法, 對(duì)已有方案的改進(jìn)優(yōu)化研究相對(duì)較少。文中在不改變?cè)挤桨赴惭b尺寸的前提下, 對(duì)輕型魚雷雷載計(jì)算機(jī)隔振系統(tǒng)進(jìn)行優(yōu)化。提出了一種基于3D設(shè)計(jì)工具和有限元分析工具的聯(lián)合優(yōu)化方法, 并使用多目標(biāo)遺傳算法對(duì)隔振系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)。

        1 雷載計(jì)算機(jī)隔振系統(tǒng)

        雷載計(jì)算機(jī)隔振系統(tǒng)如圖1所示, 4組8個(gè)橡膠隔振器以雷載計(jì)算機(jī)質(zhì)心為中心對(duì)稱分布, 將計(jì)算機(jī)懸浮支撐起來(lái)。由于文中針對(duì)雷載計(jì)算機(jī)整機(jī)進(jìn)行隔振設(shè)計(jì), 不考慮計(jì)算機(jī)內(nèi)部結(jié)構(gòu)的影響, 故將雷載計(jì)算機(jī)用質(zhì)量、質(zhì)心位置、裝配方式完全相同的質(zhì)量體代替, 以達(dá)到簡(jiǎn)化模型、減少計(jì)算量的目的。為便于施加邊界條件及說(shuō)明計(jì)算結(jié)果, 定義直角坐標(biāo)系如圖1所示。

        隔振系統(tǒng)的裝配方式如圖2所示, 橡膠隔振器為臺(tái)階圓柱結(jié)構(gòu), 臺(tái)階面與雷載計(jì)算機(jī)接觸, 內(nèi)圓柱面與大端面及套管接觸, 安裝螺釘穿過(guò)金屬套管將計(jì)算機(jī)與基座連接。

        2 建立參數(shù)化計(jì)算模型

        2.1 參數(shù)化幾何模型

        3D設(shè)計(jì)工具UG具有強(qiáng)大的參數(shù)化幾何建模能力, 并可通過(guò)WAVE幾何尺寸鏈接器建立部件間的尺寸鏈接關(guān)系, 實(shí)現(xiàn)裝配體的參數(shù)化建模[5]。隔振器是文中的優(yōu)化對(duì)象, 故在不改變安裝尺寸的前提下(圖3中的?7), 將隔振器的其余尺寸設(shè)為參數(shù)化尺寸。參數(shù)化尺寸的定義與取值范圍分別如圖3和表1所示。為了使優(yōu)化結(jié)果滿足工程設(shè)計(jì)需求, 將參數(shù)化尺寸定義為離散值, 在取值范圍內(nèi)以0.5 mm為步長(zhǎng)改變。

        表1 隔振器參數(shù)化尺寸取值范圍

        通過(guò)WAVE幾何尺寸鏈接器, 建立其余部件相關(guān)尺寸與參數(shù)化尺寸的函數(shù)關(guān)系, 即可得到隔振系統(tǒng)的參數(shù)化幾何模型。圖4為模型隨參數(shù)變化的情況。

        2.2 參數(shù)化模型有限元建模

        ANSYS Workbench可實(shí)現(xiàn)與UG的無(wú)縫對(duì)接, 通過(guò)Design Modeler模塊讀取UG幾何文件并激活其中的參數(shù)化尺寸, 即可建立ANSYS與UG聯(lián)合仿真的數(shù)據(jù)鏈接[6]。隨后對(duì)導(dǎo)入的幾何文件進(jìn)行有限元建模, 主要有邊界條件設(shè)置、材料屬性賦值和網(wǎng)格劃分3個(gè)步驟。

        邊界條件依據(jù)隔振系統(tǒng)的裝配關(guān)系定義(參見(jiàn)圖2), 金屬套管上、下端面設(shè)置為固定約束, 隔振器與其余部件的接觸關(guān)系設(shè)置為無(wú)摩擦接觸。

        雷載計(jì)算機(jī)隔振系統(tǒng)的材料參數(shù)如表2所示, 需要說(shuō)明的是, 在建模中雷載計(jì)算機(jī)等效模型按剛體處理, 故只需要輸入密度參數(shù)即可。

        使用Solid186六面體單元?jiǎng)澐志W(wǎng)格, 單元尺寸為1 mm, 劃分結(jié)果如圖5所示。由于雷載計(jì)算機(jī)等效模型為剛體, 僅在接觸面生成網(wǎng)格, 有效的減少了網(wǎng)格數(shù)量, 從而提高了優(yōu)化求解速度。

        表2 材料參數(shù)

        3 原始設(shè)計(jì)方案模態(tài)分析與試驗(yàn)驗(yàn)證

        隔振系統(tǒng)優(yōu)化設(shè)計(jì)與隔振效能評(píng)估均需要以獲取系統(tǒng)的固有頻率與振型等模態(tài)參數(shù)為基礎(chǔ)。本章將以原始設(shè)計(jì)方案為研究對(duì)象, 通過(guò)模態(tài)分析與模態(tài)試驗(yàn)對(duì)比的方式驗(yàn)證有限元模型, 確保隔振系統(tǒng)優(yōu)化設(shè)計(jì)與評(píng)估時(shí)使用模型的正確性。同時(shí), 分析結(jié)果也將用于與優(yōu)化方案的對(duì)比。

        3.1 模態(tài)分析

        原始設(shè)計(jì)方案的參數(shù)化尺寸取值如表3所示。輸入尺寸后刷新ANSYS項(xiàng)目, 即可更新有限元模型。模態(tài)分析使用ANSYS Workbench模態(tài)分析模塊, 選用迭代求解方式, 計(jì)算當(dāng)前方案隔振系統(tǒng)的前3階模態(tài)。

        表3 原始設(shè)計(jì)方案尺寸值

        3.2 模態(tài)試驗(yàn)

        模態(tài)試驗(yàn)采用錘擊法測(cè)量雷載計(jì)算機(jī)隔振系統(tǒng)在固支狀態(tài)下的模態(tài), 以帶螺紋孔的剛性基座為固支基礎(chǔ), 并使用一個(gè)轉(zhuǎn)接基座連接計(jì)算機(jī)隔振系統(tǒng)與剛性基座, 如圖6所示。4個(gè)3向加速度傳感器粘貼在計(jì)算機(jī)等效模型的上表面用于測(cè)量加速度數(shù)據(jù), 使用LMS Test. Lab測(cè)試系統(tǒng)采集與處理試驗(yàn)數(shù)據(jù)。

        3.3 模態(tài)分析與模態(tài)試驗(yàn)結(jié)果對(duì)比

        將模態(tài)分析與模態(tài)試驗(yàn)得出的隔振系統(tǒng)前3階模態(tài)參數(shù)進(jìn)行對(duì)比, 結(jié)果見(jiàn)表4。

        表4 固有頻率對(duì)比

        由圖7可知, 前3階模態(tài)振型分別為計(jì)算機(jī)沿、、方向的平動(dòng), 計(jì)算得出的振型與試驗(yàn)測(cè)出的振型存在一一對(duì)應(yīng)關(guān)系。由表4可知, 前3階模態(tài)的固有頻率誤差率較小, 最大誤差為6.21%。因此認(rèn)為, 雷載計(jì)算機(jī)隔振系統(tǒng)參數(shù)化計(jì)算模型建模方法正確, 可用于隔振系統(tǒng)優(yōu)化設(shè)計(jì)。

        4 隔振系統(tǒng)優(yōu)化設(shè)計(jì)目標(biāo)

        4.1 隔振設(shè)計(jì)基本理論

        單自由度隔振系統(tǒng)的動(dòng)力學(xué)模型如圖8所示。

        其中隔振體質(zhì)量為, 隔振體質(zhì)心位移為, 基座位移為, 隔振系統(tǒng)剛度為, 阻尼為。隔振系統(tǒng)動(dòng)力學(xué)方程如下

        (2)

        (4)

        對(duì)方程(1)進(jìn)行拉普拉斯變換得到系統(tǒng)的頻響函數(shù)

        進(jìn)行極坐標(biāo)轉(zhuǎn)換, 得到

        (6)

        稱為系統(tǒng)的幅頻特性, 將式(3)、式(4)代入式(5), 并由式(6)得到

        (8)

        將式(8)代入式(7)后便得到單自由度隔振系統(tǒng)的絕對(duì)傳遞率T()解析式

        T()的函數(shù)曲線如圖9所示。

        圖9 單自由度隔振系統(tǒng)絕對(duì)傳遞率曲線

        Fig. 9 Absolute transmission rate curves of vibration isolation system with single DOF

        4.2 隔振系統(tǒng)優(yōu)化設(shè)計(jì)目標(biāo)函數(shù)

        由隔振理論可知, 系統(tǒng)的隔振能力隨著頻率比的升高而提高。因此在外界激勵(lì)頻譜不變的情況下, 降低隔振系統(tǒng)的固有頻率即可提高隔振能力。對(duì)于多自由度隔振系統(tǒng)而言, 有效隔振頻率大于最高階固有頻率的倍, 故應(yīng)以降低最高階固有頻率作為優(yōu)化目標(biāo)。同時(shí), 若多自由度隔振系統(tǒng)滿足3向等剛度要求, 即隔振系統(tǒng)在,,3個(gè)方向的固有頻率相等或較為接近, 隔振系統(tǒng)在3個(gè)方向的振動(dòng)傳遞率曲線趨于重疊, 共振峰收窄, 亦可利于提高隔振系統(tǒng)的隔振能力。設(shè)隔振系統(tǒng)的前3階固有頻率為,,, 則隔振系統(tǒng)的優(yōu)化目標(biāo)如下:

        5 隔振系統(tǒng)優(yōu)化設(shè)計(jì)

        5.1 多目標(biāo)遺傳算法

        多目標(biāo)遺傳算法是基于達(dá)爾文進(jìn)化論和遺傳學(xué)機(jī)理, 通過(guò)模擬生物進(jìn)化過(guò)程搜索最優(yōu)解的仿生學(xué)算法。多目標(biāo)遺傳算法因其優(yōu)秀的全局尋優(yōu)能力和良好的魯棒性, 在優(yōu)化設(shè)計(jì)中被廣泛使用[8], 其求解流程如圖10所示。

        算法啟動(dòng)后將在設(shè)計(jì)變量取值范圍內(nèi)隨機(jī)生成個(gè)個(gè)體作為初始種群; 依據(jù)個(gè)體的適應(yīng)度從種群對(duì)種群中的個(gè)體進(jìn)行排序; 隨后借助自然遺傳學(xué)的交叉、變異算子, 產(chǎn)生下一代種群+1; 通過(guò)收斂準(zhǔn)則判斷是否終止運(yùn)算, 若終止, 則進(jìn)化中得到的具有最大適應(yīng)度的個(gè)體為最優(yōu)解, 若不終止, 則繼續(xù)進(jìn)化過(guò)程, 直至收斂[9]。

        5.2 求解參數(shù)設(shè)置及求解過(guò)程

        使用ANSYS Workbench Design Exploration模塊進(jìn)行隔振系統(tǒng)優(yōu)化設(shè)計(jì), 將2.1節(jié)中參數(shù)化尺寸取值范圍及增量步長(zhǎng)定義為約束條件, 將4.2節(jié)中的目標(biāo)函數(shù)設(shè)置為求解目標(biāo), 選擇多目標(biāo)遺傳算法求解。求解器參數(shù)設(shè)置[6, 10-12], 如表5所示。

        求解參數(shù)設(shè)置完成后啟動(dòng)求解器, 經(jīng)過(guò)7代進(jìn)化后種群滿足收斂準(zhǔn)則, 求解結(jié)束, 如圖11所示。目標(biāo)函數(shù)的變化情況如圖12所示。得出的最優(yōu)解尺寸參數(shù)如表6所示。

        5.3 優(yōu)化結(jié)果評(píng)估

        原始方案與優(yōu)化后的設(shè)計(jì)方案隔振性能對(duì)比如表7所示。

        表5 求解器參數(shù)

        表6 最優(yōu)解尺寸值

        振動(dòng)絕對(duì)傳遞率曲線對(duì)比如圖13所示。

        通過(guò)對(duì)比可知, 優(yōu)化方案較原始方案的最高階固有頻率降低了12.9%, 頻率間隔降低了79.6%, 系統(tǒng)的有效隔振頻率增加了39.1 Hz, 達(dá)到了增加有效隔振頻率的要求。優(yōu)化方案隔振系統(tǒng),,3個(gè)方向的振動(dòng)傳遞率曲線趨于重合, 達(dá)到了3向等剛度的設(shè)計(jì)要求。

        表7 隔振性能對(duì)比

        6 結(jié)束語(yǔ)

        針對(duì)雷載計(jì)算機(jī)隔振系統(tǒng)優(yōu)化問(wèn)題, 文中提出了一種基于3D設(shè)計(jì)軟件UG和有限元分析軟件ANSYS Workbench的聯(lián)合優(yōu)化設(shè)計(jì)方法。建立了聯(lián)合仿真參數(shù)化模型, 計(jì)算了原始設(shè)計(jì)方案的模態(tài)參數(shù), 通過(guò)與模態(tài)試驗(yàn)結(jié)果對(duì)比, 驗(yàn)證了計(jì)算模型的正確性。隨后依據(jù)隔振理論確定了目標(biāo)函數(shù), 使用多目標(biāo)遺傳算法對(duì)隔振系統(tǒng)進(jìn)行了優(yōu)化設(shè)計(jì), 達(dá)到了良好的設(shè)計(jì)效果。文中使用參數(shù)優(yōu)化方法實(shí)現(xiàn)了隔振器結(jié)構(gòu)優(yōu)化, 日后在進(jìn)行相關(guān)研究時(shí), 亦可使用無(wú)參數(shù)優(yōu)化方法, 開(kāi)展對(duì)結(jié)構(gòu)的拓?fù)湫螤顑?yōu)化研究。

        [1] 尹韶平, 劉瑞生. 魚雷總體技術(shù)[M]. 北京: 國(guó)防工業(yè)出版社, 2011.

        [2] 周亞?wèn)|, 雷宏杰, 董萼良, 等. 慣性導(dǎo)航平臺(tái)橡膠減振器斜角布置方法[J]. 振動(dòng)測(cè)試與診斷, 2014, 34(3): 447-451.Zhou Ya-dong, Lei Hong-jie, Dong E-liang, et al. Study on Placement Method for Rubber Dampers of Inertial Navigation Platform[J]. Journal of Vibration, Measurement & Diagnosis, 2014, 34(3): 447-451.

        [3] Kim S M, Elliott S J, Brennan M J.Decentralized Control for Multichannel Active Vibration Isolation[J]. IEEE Transactions on Control Systems Technology, 2001, 9(1): 93-100.

        [4] Vulcan A, Bloch M. A Low Noise Vibration Isolated Airborne Radar Synthesizer[C]//45th Annual Symposium on Frequency Control. [s.l.]: IEEE, 1991: 330-336.

        [5] 宋昌才. UG NX8.5標(biāo)準(zhǔn)教程[M]. 北京: 科學(xué)出版社, 2016.

        [6] ANSYS Inc. ANSYS Workbench Help Document[M]. Canonsburg: ANSYS Inc, 2014.

        [7] 丁文鏡. 減振理論[M]. 北京: 清華大學(xué)出版社, 2014.

        [8] Wang X D, Hirsch C, Kang S, et al. Multi-objective Optimization of Turbomachinery Using Improved NSGA-II and Approximation Model[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(9): 883-895.

        [9] 周明, 孫樹(shù)棟. 遺傳算法原理及應(yīng)用[M]. 北京: 國(guó)防工業(yè)出版社, 1999.

        [10] Yu L J, Liu S Y, Liu F M, et al. Energy Optimization of the Fin/Rudder Roll Stabilization System Based on the Multi-objective Genetic Algorithm[J]. Journal of Marine Science and Application, 2015(14): 202-207.

        [11] 李楠, 王明輝, 馬書根, 等. 基于多目標(biāo)遺傳算法的水陸兩棲可變形機(jī)器人結(jié)構(gòu)參數(shù)設(shè)計(jì)方法[J]. 機(jī)械工程學(xué)報(bào), 2012, 48(17): 10-20.Li Nan, Wang Ming-hui, Ma Shu-gen, et al. Mechanism- parameters Design Method of an Amphibious Transformable Robot Based on Multi-objective Genetic Algorithm[J]. Journal of Mechanical Engineering, 2012, 48(17): 10-20.

        [12] 李成冬, 孟春毅. 基于多目標(biāo)遺傳算法的汽車玻璃升降器結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 機(jī)械設(shè)計(jì)與制造, 2016(3): 253-256. Li Cheng-dong, Meng Chun-yi. Structure Research on Automotive Glass Riser with Multi-Objective Genetic Algorithm[J]. Machinery Design & Manufacture, 2016(3): 253-256.

        (責(zé)任編輯: 許 妍)

        Optimization Design of Vibration Isolation System for Torpedo Borne Computer Based on Multi-objective Genetic Algorithm

        ZHANG Yi-wen1,2, YIN Shao-ping1, WANG Zhi-jie1, GUO Jun1, WU Pei-ming1, CUI Xin-shan1,2

        (1. The 705 Research Institute, China Shipbuilding Industry Corporation, Xi′an 710077, China; 2. Science and Technology on Underwater Information and Control Laboratory, Xi′an 710077, China)

        To solve the vibration isolation problem of torpedo borne computer in complex mechanical environment, this paper proposes a joint optimization design method based on UG and ANSYS Workbench. A parametric geometry model is established by UG and imported into ANSYS Workbench for finite element modeling. Modal parameters of the original design scheme are calculated and compared with the modal test data to prove the correctness of the calculation model. Then the optimal objective function is determined based on vibration isolation theory, and optimization design of the vibration isolation system is carried out using multi-objective genetic algorithm. Compared with the original design, the optimization scheme reduces the maximum natural frequency by 12.9% and the frequency interval by 79.6%, and increases the effective vibration isolation frequency by 39.1 Hz.

        torpedo borne computer; optimization design; multi-objective genetic algorithm; vibration isolation; modal theory

        TJ630.3; TB115.2

        A

        2096-3920(2017)01-0057-07

        10.11993/j.issn.2096-3920.2017.01.005

        2016-11-19;

        2016-12-18.

        張溢文(1988-), 男, 在讀碩士, 主要研究方向?yàn)轸~雷總體設(shè)計(jì).

        [引用格式] 張溢文, 尹韶平, 王志杰, 等. 基于多目標(biāo)遺傳算法的雷載計(jì)算機(jī)隔振系統(tǒng)優(yōu)化設(shè)計(jì)[J]. 水下無(wú)人系統(tǒng)學(xué)報(bào), 2017, 25(1): 57-63.

        猜你喜歡
        模態(tài)計(jì)算機(jī)優(yōu)化
        超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
        計(jì)算機(jī)操作系統(tǒng)
        民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
        關(guān)于優(yōu)化消防安全告知承諾的一些思考
        一道優(yōu)化題的幾何解法
        基于計(jì)算機(jī)自然語(yǔ)言處理的機(jī)器翻譯技術(shù)應(yīng)用與簡(jiǎn)介
        科技傳播(2019年22期)2020-01-14 03:06:34
        信息系統(tǒng)審計(jì)中計(jì)算機(jī)審計(jì)的應(yīng)用
        國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
        基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
        Fresnel衍射的計(jì)算機(jī)模擬演示
        国产三级精品和三级男人| 99国产精品丝袜久久久久| 国产精品一区二区AV不卡| 国产自拍视频免费在线观看| 国产成人午夜高潮毛片| 男人扒开女人下面狂躁小视频 | 神马不卡一区二区三级| 亚洲国产精品久久性色av| 丰满少妇人妻无码| 东北寡妇特级毛片免费| 国产午夜精品久久久久99| 亚洲中字永久一区二区三区| 免费a级毛片18禁网站| 精品人妻人人做人人爽| 中文字幕人成人乱码亚洲| 少妇人妻av一区二区三区| 九九在线中文字幕无码| 久久久久亚洲av无码专区网站| 最新亚洲无码网站| 亚洲禁区一区二区三区天美| 精品人妻少妇嫩草av无码专区| 中文乱码人妻系列一区二区| 国产精品亚洲美女av网站| 日本在线观看不卡一区二区| 午夜精品久久久久久毛片| 亚洲中文字幕每日更新| 在线亚洲免费精品视频| 中文无码av一区二区三区| 国产精品免费大片| 最新手机国产在线小视频| 久久久精品亚洲人与狗| 日韩精品久久久久久久电影蜜臀| 国产啪精品视频网站丝袜| 一区二区免费国产a在亚洲| 丰满人妻一区二区三区视频| 人妻系列无码专区久久五月天| 久久久久久岛国免费网站| 国产的自拍av免费的在线观看| 天堂新版在线资源| 国产成人亚洲综合无码精品| 国产精品后入内射日本在线观看|