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

        ?

        潛器慣性類水動(dòng)力參數(shù)計(jì)算方法研究*

        2016-06-30 07:06:48于秋禮
        艦船電子工程 2016年6期
        關(guān)鍵詞:潛器

        高 峰 于秋禮

        (1.91439部隊(duì) 大連 116041)(2.中船重工第七一〇研究所 宜昌 443003)

        潛器慣性類水動(dòng)力參數(shù)計(jì)算方法研究*

        高峰1于秋禮2

        (1.91439部隊(duì)大連116041)(2.中船重工第七一〇研究所宜昌443003)

        摘要為精確獲取潛器水動(dòng)力參數(shù),對(duì)慣性類水動(dòng)力參數(shù)常用計(jì)算方法進(jìn)行了總結(jié),給出了方程解析法、半經(jīng)驗(yàn)公式法和水池試驗(yàn)法優(yōu)缺點(diǎn)。在對(duì)慣性水動(dòng)力參數(shù)計(jì)算機(jī)計(jì)算方法研究的基礎(chǔ)上,提出一種以Hess-Smith面元法為理論依據(jù),獲取潛器慣性類水動(dòng)力參數(shù)的方法。編寫慣性類水動(dòng)力參數(shù)數(shù)值計(jì)算程序,計(jì)算結(jié)果與圓球、橢球,等簡(jiǎn)單幾何體慣性類水動(dòng)力參數(shù)理論值對(duì)比,證明計(jì)算方法及所編程序的正確性,并探討了網(wǎng)格劃分形式和網(wǎng)格數(shù)量對(duì)計(jì)算精度的影響。最后以復(fù)雜構(gòu)型潛器為研究對(duì)象,計(jì)算其所受的慣性類水動(dòng)力參數(shù),與半經(jīng)驗(yàn)公式結(jié)果比較及定性分析,得知數(shù)值計(jì)算結(jié)果更具合理性。

        關(guān)鍵詞邊界元法; 水動(dòng)力參數(shù); 附加質(zhì)量; 潛器

        Class NumberTJ61+7

        1引言

        潛器近年來(lái)在軍事、勘探、救援等領(lǐng)域得到廣泛應(yīng)用,潛器研究受到更多關(guān)注,潛器水動(dòng)力特性是潛器研究重要課題之一。潛器的水動(dòng)力特性對(duì)其整體布局、航行控制、能源需求、水下作業(yè)甚至布放回收均具有重要影響;同時(shí)潛器水動(dòng)力參數(shù)獲取也是其六自由空間動(dòng)態(tài)仿真和運(yùn)動(dòng)控制的前提與基礎(chǔ);另外對(duì)潛器進(jìn)行減振降噪研究也需考慮潛器的水動(dòng)力特性。由此可見,潛器水動(dòng)力參數(shù)的精確對(duì)潛器的研發(fā)與設(shè)計(jì)具有重要意義。

        陳琪瑋給出了辨識(shí)航行體水動(dòng)力參數(shù)的智能計(jì)算方法[1],施小成給出利用遞推最小二乘法估計(jì)潛器水動(dòng)力參數(shù)數(shù)學(xué)模型[2],袁偉杰等提出了基于遺傳算法的自治水下機(jī)器人水動(dòng)力參數(shù)計(jì)算方法[3[4],李天森和S K Lee等分別給出基于水池試驗(yàn)的水動(dòng)力參數(shù)測(cè)量計(jì)算方法[5~6]。

        本文以復(fù)雜潛器慣性水動(dòng)力參數(shù)計(jì)算為研究對(duì)象,依據(jù)面元法理論,應(yīng)用計(jì)算機(jī)數(shù)值計(jì)算方法對(duì)潛器慣性水動(dòng)力參數(shù)獲取進(jìn)行了研究,可為潛器等相關(guān)裝備設(shè)計(jì)應(yīng)用提供參考。

        2常用慣性類水動(dòng)力計(jì)算方法

        潛器水動(dòng)力參數(shù)可分為兩種:一是和加速度有關(guān)的慣性類水動(dòng)力參數(shù),即附加質(zhì)量矩陣;二是和速度有關(guān)的粘性類水動(dòng)力參數(shù)[7~10]。粘性類水動(dòng)力參數(shù)可通過(guò)CFD軟件模擬潛器拖曳水池和懸臂水池試驗(yàn)求得,并已取得良好的數(shù)值計(jì)算精度[8~10]。

        潛器慣性類水動(dòng)力參數(shù)常用計(jì)算方法歸納總結(jié)主要包括以下幾種:方程解析、半經(jīng)驗(yàn)公式、水池試驗(yàn)測(cè)量、數(shù)值計(jì)算。

        2.1方程解析

        利用流體力學(xué)理論求解解析值,目前除了球體、橢球體、圓柱體等簡(jiǎn)單構(gòu)型幾何體的慣性類水動(dòng)力參數(shù)可直接由其繞流場(chǎng)的速度勢(shì)求得外,復(fù)雜構(gòu)型潛器的慣性類水動(dòng)力參數(shù)則很難求得。

        2.2半經(jīng)驗(yàn)公式

        半經(jīng)驗(yàn)公式估算是一種較為粗糙的計(jì)算方法,它提供的水動(dòng)力參數(shù)是在簡(jiǎn)單幾何形體理論計(jì)算公式的基礎(chǔ)上作一些修正得出的,或是根據(jù)在大量船模實(shí)驗(yàn)結(jié)果基礎(chǔ)上經(jīng)整理得出的半經(jīng)驗(yàn)的圖譜估算得到,在估算細(xì)長(zhǎng)體時(shí)效果很好,但復(fù)雜構(gòu)型潛器附體較多,外形各異,采用艦船、潛艇、魚雷等提供的半經(jīng)驗(yàn)公式,這不可避免地給潛器,尤其是復(fù)雜構(gòu)型潛器的水動(dòng)力參數(shù)估算帶來(lái)了很大的誤差。

        2.3水池試驗(yàn)測(cè)量

        水池試驗(yàn)測(cè)量慣性類水動(dòng)力參數(shù),常用的測(cè)量方式包括周期比較法、慣性法、電比擬法、平面運(yùn)動(dòng)機(jī)構(gòu)等。但就當(dāng)前來(lái)說(shuō),水池試驗(yàn)測(cè)量費(fèi)用高,周期長(zhǎng)、不利于潛器,尤其是低成本潛器的開發(fā)和研究。

        2.4數(shù)值求解

        隨著計(jì)算機(jī)計(jì)算能力的不斷提高,計(jì)算流體力學(xué)已成為流體力學(xué)研究中的一個(gè)重要分支和輔助解決工程問(wèn)題有效手段。計(jì)算機(jī)數(shù)值求解的方法又可分為兩種:一是運(yùn)用商用集成計(jì)算軟件求解,二是通過(guò)分析方程編程求解計(jì)算。

        2.4.1商用集成計(jì)算軟件

        用于求解慣性類水動(dòng)力參數(shù)的商業(yè)軟件有很多,其中比較常見的有以下幾種:AQWA、HYDROSTAR、MOSES、NEPTUNE、WADAM、WAMIT、WAVELOAD等,以上流體力學(xué)軟件包(Hydrodynamic Software Packages)多用以時(shí)域和頻域勢(shì)流理論為基礎(chǔ),計(jì)算艦船以及海上平臺(tái)波載作用下的水動(dòng)力問(wèn)題;慣性類水動(dòng)力參數(shù)多應(yīng)用頻域勢(shì)流理論求得[7]。

        大型計(jì)算流體動(dòng)力學(xué)軟件(CFD)fluent等,利用其動(dòng)網(wǎng)格技術(shù)使物體作振蕩運(yùn)動(dòng)求物體的受力,可用于潛器的慣性類水動(dòng)力參數(shù)的逐個(gè)計(jì)算[12~13]。但CFD建模相對(duì)復(fù)雜,劃分網(wǎng)格困難,且對(duì)網(wǎng)格質(zhì)量也有要求,否則使用動(dòng)網(wǎng)格時(shí)會(huì)錯(cuò)誤中止。而面元法網(wǎng)格只在潛器構(gòu)型表面劃分面元,建模相對(duì)簡(jiǎn)單,網(wǎng)格達(dá)到一定數(shù)量后,便可獲得良好的計(jì)算精度;計(jì)算簡(jiǎn)便,一次可完成潛器36個(gè)慣性類水動(dòng)力參數(shù)的全部計(jì)算,較CFD計(jì)算方法具有明顯的優(yōu)勢(shì)。

        此外氣動(dòng)力仿真軟件系統(tǒng)USAERO也可用于慣性類水動(dòng)力參數(shù)的獲得,并取得了良好的計(jì)算精度[14]。

        2.4.2編程數(shù)值計(jì)算求解

        編程數(shù)值計(jì)算求解慣性類水動(dòng)力參數(shù)的理論主要有兩種:1)切片理論(Strip Theory),基于波浪輻射與衍射沿潛器或艦船縱向方向的變化極其緩慢這一假定,對(duì)紐曼-開爾文方程(N-K方程)進(jìn)行一定合理的簡(jiǎn)化,應(yīng)用求解三維紐曼-開爾文方程(N-K方程)的近似方法計(jì)算。在假設(shè)潛器構(gòu)型為細(xì)長(zhǎng)體即潛器橫向尺寸遠(yuǎn)小于縱向尺寸的前提下,其計(jì)算結(jié)果可以滿足工程應(yīng)用。但對(duì)于橫截面沿縱向變化大的復(fù)雜構(gòu)型潛器其假設(shè)不成立,計(jì)算精度大為降低。2)邊界元法,將物體表面劃分單元,在每一個(gè)單元上布置強(qiáng)度待定的基本解(像源、匯、偶極子等),從而確定流場(chǎng)速度勢(shì),進(jìn)而得到慣性類水動(dòng)力參數(shù)。其中Hess-Smith面元法就是應(yīng)用邊界元法第三類方法:間接表達(dá)式進(jìn)行求解的。將積分方程完全用微分方程的單位奇異解表示,奇異解對(duì)應(yīng)的奇點(diǎn)以特定強(qiáng)度分布在感興趣的邊界上,奇點(diǎn)密度函數(shù)其本身并無(wú)具體的物理含義,但一旦從積分方程數(shù)值求出密度函數(shù),但只要作一些積分計(jì)算就可以得到物體內(nèi)任意一點(diǎn)處解參數(shù)的值。

        Hess-Smith面元法求解數(shù)值計(jì)算理論假定流體為無(wú)粘、無(wú)旋、不可壓縮的理想流體,這一假設(shè)與慣性類水動(dòng)力參數(shù)的定義一致。運(yùn)用Hess-Smith面元法計(jì)算復(fù)雜構(gòu)型潛器的慣性類水動(dòng)力參數(shù)在理論上可得到有力的支撐[15]。

        3Hess-Smith面元法慣性類水動(dòng)力計(jì)算方法

        采用Hess-Smith面元法編程,其理論依據(jù)以及數(shù)值模型詳細(xì)可參考文獻(xiàn)[15~16]。這一方法用多平面四邊形表面單元近似表示物面,每個(gè)單元上布置一強(qiáng)度未知的源,然后在物體表面的某些考察點(diǎn)上滿足法向速度為零的物面邊界條件,得到求單元源密度的線性代數(shù)方程組。求解方程組得到源密度分布,進(jìn)而可求流場(chǎng)內(nèi)任意點(diǎn)的速度、壓力等物理量。因此Hess-Smith面元法計(jì)算程序計(jì)算精度僅與數(shù)值模型建立的精準(zhǔn)有關(guān)。

        3.1Hess-Smith面元法慣性類水動(dòng)力程序計(jì)算流程

        慣性類水動(dòng)力數(shù)值計(jì)算程序以Hess-Smith面元法為計(jì)算核心,其計(jì)算流程如圖1所示:可大致分為前處理、求解計(jì)算以及后處理三部分。

        圖1 慣性類水動(dòng)力參數(shù)計(jì)算流程

        其中前處理包括幾何建模、網(wǎng)格劃分、數(shù)據(jù)輸入三部分:幾何建模主要是在CAD軟件中完成模型的建立;網(wǎng)格劃分主要是完成Hess-Smith計(jì)算程序所需的面元網(wǎng)格劃分;而數(shù)據(jù)輸入則完成面元網(wǎng)格的整理以及Hess-Smith計(jì)算程序所需格式數(shù)據(jù)的準(zhǔn)備工作。求解計(jì)算是整個(gè)程序的核心,綜合考慮計(jì)算機(jī)硬件和數(shù)值計(jì)算時(shí)間和計(jì)算精度等方面的因素,計(jì)算程序目前最大面元網(wǎng)格數(shù)為5500。后處理主要包括結(jié)果輸出和圖形顯示兩部分:結(jié)果輸出完成各面元網(wǎng)格的速度、壓力等信息以及潛器整體的慣性類水動(dòng)力參數(shù)及無(wú)量綱系數(shù)的輸出;圖形顯示完成各面元上速度以及壓力的顯示,進(jìn)而根據(jù)結(jié)果圖可在壓力或速度變化急劇的部位加密網(wǎng)格,對(duì)網(wǎng)格布局進(jìn)行優(yōu)化,進(jìn)一步提高計(jì)算精度。

        3.2Hess-Smith面元法慣性類水動(dòng)力程序計(jì)算精度

        為驗(yàn)證 Hess-Smith慣性類水動(dòng)力計(jì)算程序的有效性和計(jì)算精度,建立圓球和橢球網(wǎng)格模型,同時(shí)為比較不同網(wǎng)格劃分形式對(duì)計(jì)算精度的影響,圓球分別采用球?qū)ΨQ和柱對(duì)稱兩種網(wǎng)格劃分方式,其網(wǎng)格模型效果圖分別如圖2~圖4所示。計(jì)算誤差與網(wǎng)格數(shù)量之間的關(guān)系如圖5、圖6所示。

        圖2 圓球?qū)ΨQ網(wǎng)格模型

        圖3 圓球柱對(duì)稱網(wǎng)格模型

        圖4 橢球網(wǎng)格模型圖

        圖5 圓球球?qū)ΨQ模型與柱對(duì)稱模型精度對(duì)比

        圖6 橢球網(wǎng)格數(shù)目與計(jì)算精度關(guān)系圖

        從圖5可以看出,圓球柱對(duì)稱網(wǎng)格模型各坐標(biāo)軸上的網(wǎng)格劃分的差異性導(dǎo)致理論上各軸方向本應(yīng)相等的慣性水動(dòng)力參數(shù)出現(xiàn)不等的數(shù)值計(jì)算結(jié)果,而球?qū)ΨQ網(wǎng)格模型在網(wǎng)格劃分上考慮圓球三軸的對(duì)稱性,其各軸方向慣性水動(dòng)力參數(shù)結(jié)算結(jié)果相等;從整體趨勢(shì)上可以看出不管是球?qū)ΨQ網(wǎng)格模型還是柱對(duì)稱網(wǎng)格模型,其計(jì)算精度都隨著網(wǎng)格數(shù)目的增多而不斷提高;當(dāng)網(wǎng)格達(dá)到一定數(shù)量后,再行增加網(wǎng)格的數(shù)量對(duì)計(jì)算精度的提高將不會(huì)有太大影響。此外從圖6也可以看出,慣性類水動(dòng)力參數(shù)的計(jì)算精度隨著網(wǎng)格數(shù)量的增多而不斷提高。且其計(jì)算精度與理論值相比誤差不大于2%,計(jì)算精度良好,完全滿足工程應(yīng)用需要。

        綜上所述,可以得出以下兩條基本結(jié)論:1)網(wǎng)格模型的劃分要跟據(jù)實(shí)際模型的對(duì)稱性對(duì)稱劃分;2)在計(jì)算機(jī)硬件和計(jì)算時(shí)間允許的前提下盡可能增大模型網(wǎng)格的數(shù)量,盡量對(duì)模型進(jìn)行精細(xì)劃分,尤其是在壓力和速度變化急劇的面元部分,以提高數(shù)值計(jì)算精度。

        3.3潛器Hess-Smith慣性類水動(dòng)力計(jì)算結(jié)果以及分析

        以一潛器為具體模型實(shí)例,對(duì)其進(jìn)行慣性類水動(dòng)力參數(shù)求解計(jì)算,其網(wǎng)格模型如圖7所示,共計(jì)劃分5374個(gè)面元網(wǎng)格,經(jīng)計(jì)算其數(shù)值計(jì)算結(jié)果見表1。

        圖7 潛器網(wǎng)格模型

        圖8 潛器面元壓力分布

        表1 Hess-Smith面元法數(shù)值計(jì)算數(shù)據(jù)

        表2 半經(jīng)驗(yàn)公式計(jì)算數(shù)據(jù)

        由于該潛器目前尚未進(jìn)行水池水動(dòng)力測(cè)量試驗(yàn),其慣性類水動(dòng)力參數(shù)僅通過(guò)半經(jīng)驗(yàn)公式估算得到,但通過(guò)圓球、橢球計(jì)算結(jié)果精度以及上述定性的分析可知,Hess-Smith面元法數(shù)值計(jì)算結(jié)果數(shù)據(jù)應(yīng)更可取。

        4結(jié)語(yǔ)

        本文通過(guò)對(duì)慣性類水動(dòng)力參數(shù)計(jì)算方法優(yōu)缺點(diǎn)的對(duì)比與歸納,以Hess-Smith面元法為理論依據(jù)編寫數(shù)值計(jì)算程序,提供了一種精確計(jì)算潛器慣性類水動(dòng)力參數(shù)的手段;通過(guò)簡(jiǎn)單模型的數(shù)值計(jì)算結(jié)果與理論值的對(duì)比,驗(yàn)證了計(jì)算方法和程序的有效性和準(zhǔn)確性;通過(guò)對(duì)一潛器建模以及計(jì)算,得到了其慣性類水動(dòng)力參數(shù)即附加質(zhì)量矩陣,與半經(jīng)驗(yàn)公式結(jié)果對(duì)比分析,數(shù)值計(jì)算結(jié)果更可取。但是本文沒(méi)有考慮有限水深、自由液面、潛器附近物體等對(duì)潛器附加質(zhì)量的影響,后續(xù)可作進(jìn)一步分析與研究。

        參 考 文 獻(xiàn)

        [1] 陳瑋琪,顏開,王寶壽.辨識(shí)航行體水動(dòng)力參數(shù)的智能方法[J].船舶力學(xué),2011,15(4):359-363.

        [2] 施小成,趙俊成,周佳加.UUV水動(dòng)力參數(shù)估計(jì)及深度自抗擾控制方法[J].儀器儀表學(xué)報(bào),2014,46(11):1-6.

        [3] 袁偉杰,劉貴杰,朱紹峰.基于遺傳算法的自治水下機(jī)器人水動(dòng)和參數(shù)辨識(shí)方法[J].機(jī)械工程學(xué)報(bào),2010,35(12):96-100.

        [4] 劉建成,劉學(xué)敏,徐玉如.極大似然法在水下機(jī)器人系統(tǒng)辨識(shí)中的應(yīng)用[J].哈爾濱工程大學(xué)學(xué)報(bào),2001,22(5):1-4.

        [5] 李天森.魚雷操縱性[M].北京:國(guó)防工業(yè)出版社,1999:242-251.

        [6] S K Lee,S J Cheon,T H Joung,T S Jang, J H Lee. Evaluation of the Added Mass for a Spheroid-type UUV by VPMM Test[C]//Proceedings of the Nineteenth (2009) International Offshore and Polar Engineering Conference,2009:703-709.

        [7] 蔣新松,封錫盛.水下機(jī)器人[M].沈陽(yáng):遼寧科學(xué)技術(shù)出版社,2000:257-271.

        [8] 胡志強(qiáng),林揚(yáng),谷海濤.水下機(jī)器人粘性類水動(dòng)力數(shù)值計(jì)算方法研究[J].機(jī)器人,2007,29(2):145-150.

        [9] 谷海濤,林揚(yáng),胡志強(qiáng).帶槽道槳水下機(jī)器人阻力特性的數(shù)值分析[J].微計(jì)算機(jī)信息,2007,23(14):227-228.

        [10] 康濤,胡克,胡志強(qiáng)等.CFX與USAERO的水下機(jī)器人操縱性仿真計(jì)算研究[J].機(jī)器人,2005,27(6):535-538.

        [11] A P McCabe. An Appraisal of a Range of Fluid Modelling Software[D]. U.K.:Department of Engineering Lancaster University,2004:18-48.

        [12] 郭志彬,段文洋,趙彬彬.數(shù)值計(jì)算復(fù)雜構(gòu)型潛器附加質(zhì)量的方法[J]. 中國(guó)科技論文在線精品論文,2009, 2(5):428-435.

        [13] 馬燁,單雪雄.數(shù)值計(jì)算復(fù)雜外形物體附加質(zhì)量的新方法[J].計(jì)算機(jī)仿真,2007,24(5):75-78.

        [14] 程慶偉,張靖,李遺.氣動(dòng)力仿真軟件系統(tǒng)USAERO的計(jì)算精度研究[C]//飛行力學(xué)與飛行試驗(yàn)學(xué)術(shù)交流年會(huì)論文集.中國(guó)航空學(xué)會(huì).中國(guó)飛行試驗(yàn)研究院,2005:353-357.

        [15] J L Hess,A M O Smith. Calculation of Potential Flow about Arbitrary Bodies[J] . Progress in Aeronautical Science,1967,8(8):1-138.

        [16] 戴遺山.艦船在波浪中運(yùn)動(dòng)的頻域與時(shí)域勢(shì)流理論[M].北京:國(guó)防工業(yè)出版社,1998:15-54.

        Calculation Method of Inertial Hydrodynamic Parameters of Underwater Vehicle

        GAO Feng1YU Qiuli2

        (1.No.91439 Troops of PLA,Dalian116041)(2.No.710 R&D Institute, CSIC, Yichang443003)

        AbstractIn order to predict maneuverability and motion simulation of underwater vehicle with complex configuration,it is necessary to derive the hydrodynamic parameters in advance. So this paper aims to determine the inertial hydrodynamics for the governing equations of motion for underwater vehicle. By summarizing and comparing several computational techniques used to deriving inertial hydrodynamics, a code which is based on Hess-Smith panel elements method is programmed as a useful tool for inertial hydrodynamics numerical calculation. The numerical results of an ellipsoid and sphere calculated by the code were confirmed by the data of the theoretical value, which gives confidence that this approach will be useful for calculating inertial hydrodynamics. Essentially, this paper is looking for a general answer as to how many panels are enough over a prescribed geometry in order to accurately capture the flow effects based on different panel results. Thus inertial hydrodynamics calculated by the approach is presented and compared with the semi-empirical value, which shows that the approach can be used to calculate the inertial hydrodynamics of the underwater vehicle with complex configuration.

        Key Wordsboundary element method, hydrodynamic parameters, added mass, underwater vehicle

        *收稿日期:2015年12月6日,修回日期:2016年1月18日

        作者簡(jiǎn)介:高峰,男,工程師,研究方向:試驗(yàn)技術(shù)。

        中圖分類號(hào)TJ61+7

        DOI:10.3969/j.issn.1672-9730.2016.06.033

        猜你喜歡
        潛器
        基于改進(jìn)遞歸最小二乘估計(jì)的潛器軌跡預(yù)測(cè)
        潛器六自由度運(yùn)動(dòng)方程水動(dòng)力系數(shù)敏感性分析
        連續(xù)沖擊條件下水下潛器運(yùn)動(dòng)姿態(tài)估計(jì)方法
        潛器鉛壓載應(yīng)用技術(shù)研究
        內(nèi)孤立波作用下水下潛器的載荷特性數(shù)值分析
        水下維修載人潛器運(yùn)動(dòng)仿真
        載人潛器縱傾調(diào)節(jié)系統(tǒng)設(shè)計(jì)及FMEA分析
        潛器耐壓殼體圓度測(cè)量方法設(shè)計(jì)
        潛器在沖擊載荷下的運(yùn)動(dòng)和控制研究
        海洋工程(2014年1期)2014-10-11 06:19:28
        蛟龍?zhí)栞d人潛水器無(wú)動(dòng)力潛浮運(yùn)動(dòng)分析系統(tǒng)開發(fā)
        一个人在线观看免费视频www| a黄片在线视频免费播放| 国产无遮挡aaa片爽爽| a级特黄的片子| 精品无码AⅤ片| 国内大量揄拍人妻在线视频| 国产精品日本天堂| 国产女主播福利在线观看| 久久久久99精品成人片| 韩国精品一区二区三区无码视频| 久久国产免费观看精品| 日韩精品视频av在线观看| 成人自拍小视频在线看| 色狠狠色噜噜av天堂一区| 风流少妇又紧又爽又丰满| 亚洲精品国产精品av| 日本免费一二三区在线| 18禁无遮拦无码国产在线播放| 色综合天天网| 一级黄色一区二区三区视频| 黄片视频免费在线播放观看| 疯狂撞击丝袜人妻| 久久久久亚洲精品天堂| 日产一区一区三区区别| 国产猛烈高潮尖叫视频免费| 牛鞭伸入女人下身的真视频| 人妻中文字幕不卡精品| 在线观看国产激情视频| 97久久草草超级碰碰碰| 无码电影在线观看一区二区三区| 黄色中文字幕视频网站| 成人免费自拍视频在线观看| 无码ol丝袜高跟秘书在线观看| 欧美成人a在线网站| 国产三级在线观看高清| 精品人妻一区二区三区久久| 国产98在线 | 日韩| 99久久精品自在自看国产| 久久久婷婷综合亚洲av| 一本到亚洲av日韩av在线天堂| 日韩在线 | 中文|