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

        ?

        雷-靶碰撞結(jié)構(gòu)響應(yīng)仿真分析

        2019-03-15 01:49:58李靜肖陳彥勇
        水下無人系統(tǒng)學(xué)報 2019年1期
        關(guān)鍵詞:有限元結(jié)構(gòu)

        李靜肖, 陳彥勇, 張 濤, 許 達(dá)

        ?

        雷-靶碰撞結(jié)構(gòu)響應(yīng)仿真分析

        李靜肖, 陳彥勇, 張 濤, 許 達(dá)

        (中國船舶重工集團(tuán)公司 第705研究所昆明分部, 云南 昆明, 650118)

        針對雷-靶碰撞過程中靶體結(jié)構(gòu)及碰撞環(huán)境的特殊性與復(fù)雜性, 基于船舶碰撞內(nèi)部機(jī)理, 采用流固耦合及附加質(zhì)量法對魚雷撞擊目標(biāo)靶時的結(jié)構(gòu)響應(yīng)進(jìn)行了有限元仿真, 研究了碰撞過程中靶體的結(jié)構(gòu)損傷以及魚雷速度及加速度變化。仿真結(jié)果表明: 1) 撞擊過程具有很強(qiáng)的非線性特征, 撞擊角度的變化對雷體運(yùn)動及靶體變形有一定程度的影響; 2) 魚雷撞擊速度越大, 雷頭加速度響應(yīng)峰值越高, 靶體損傷程度越大。文中所做研究可為目標(biāo)靶的結(jié)構(gòu)設(shè)計及動態(tài)特性設(shè)計提供參考。

        雷-靶碰撞; 流固耦合; 附加質(zhì)量法; 有限元仿真; 結(jié)構(gòu)響應(yīng)

        0 引言

        隨著魚雷技術(shù)的不斷發(fā)展, 其先進(jìn)程度越來越高, 反潛訓(xùn)練對靶標(biāo)的使用要求也不斷提高。魚雷撞擊目標(biāo)靶時, 雷頭引信部位加速度需達(dá)到一定幅值及脈寬才能觸發(fā)引信動作, 因此, 靶體結(jié)構(gòu)剛度是否滿足觸發(fā)引信動作要求是靶體結(jié)構(gòu)設(shè)計的關(guān)鍵之一[1]。由于目標(biāo)靶結(jié)構(gòu)較為特殊, 同時碰撞環(huán)境中流體作用復(fù)雜, 目前國內(nèi)外尚未有對該類碰撞問題公開發(fā)表的研究。因此, 建立一套適用于該類雷-靶碰撞問題的研究方法, 并準(zhǔn)確預(yù)報碰撞過程中靶體的動力響應(yīng)對目標(biāo)靶結(jié)構(gòu)設(shè)計及以后該類問題的研究具有重要意義。對于雷-靶碰撞問題, 其碰撞機(jī)理與船舶碰撞類似, 因此可以參考船舶碰撞的分析方法, 同時考慮雷、靶自身結(jié)構(gòu)特點以及碰撞環(huán)境實現(xiàn)雷-靶碰撞結(jié)構(gòu)響應(yīng)分析研究。目前針對水面艦船以及潛器的碰撞問題, 主要有經(jīng)驗法[2]、解析法[3]、有限元數(shù)值仿真[4]和試驗方法[5]等4種方法, 文中利用ANSYS LS-DYNA非線性有限元軟件, 并結(jié)合流固耦合及附加質(zhì)量法開展雷-靶碰撞問題的研究。

        1 雷-靶碰撞仿真計算原理

        雷-靶碰撞過程中具有很強(qiáng)的非線性特性, 碰撞區(qū)域的結(jié)構(gòu)會迅速通過彈性階段進(jìn)入塑性流動狀態(tài), 并伴隨撕裂、屈曲等形式的破壞或失效[6]。因此可以采用顯式非線性有限元分析方法對該問題進(jìn)行分析。

        1.1 顯式時間積分算法

        考慮沙漏阻尼后, 雷-靶碰撞過程的動力學(xué)方程可表示為

        采用中心差分法, 可得到加速度、速度和位移的基本遞推格式

        在LS-DYNA中, 式(2)中的結(jié)構(gòu)質(zhì)量矩陣為集中質(zhì)量矩陣, 得到的動力學(xué)方程為解耦的, 因此無需進(jìn)行總體矩陣的計算以及平衡迭代, 但是該方法并不是無條件穩(wěn)定的, 為保證數(shù)值計算的穩(wěn)定性, 采用變步長的計算方法, 每一時刻的積分步長由當(dāng)前時刻的穩(wěn)定性條件進(jìn)行控制, 積分步長必須小于某一臨界值計算才能繼續(xù)進(jìn)行, 臨界時間步長由下式確定[7]

        1.2 接觸與摩擦計算

        在雷-靶碰撞的動態(tài)接觸問題中, 接觸與摩擦是2個關(guān)鍵問題, 即不同結(jié)構(gòu)界面碰撞時的接觸和相對滑動。碰撞過程中構(gòu)件的接觸通過接觸算法進(jìn)行處理, 文中采用LS-DYNA自帶的基于主-從面的動態(tài)接觸算法, 通過罰函數(shù)法計算接觸力, 接觸通過定義關(guān)鍵字實現(xiàn)自動面雙向接觸。仿真接觸中靶體設(shè)為主面, 魚雷設(shè)為從面。

        同時發(fā)生相對滑動的結(jié)構(gòu)件之間會產(chǎn)生摩擦力, 摩擦力可采用經(jīng)典庫倫摩擦進(jìn)行處理。摩擦力大小由下式確定

        1.3 水動作用力處理

        目標(biāo)靶采用濕式回轉(zhuǎn)結(jié)構(gòu), 碰撞過程中靶體內(nèi)外均有流體存在, 碰撞過程中不僅存在雷-靶結(jié)構(gòu)間的相互作用, 結(jié)構(gòu)與周圍流體也存在復(fù)雜的交互作用。結(jié)構(gòu)變形時會引起流體的運(yùn)動, 繼而導(dǎo)致流體作用于結(jié)構(gòu)的載荷大小及分布發(fā)生變化, 相反, 流體載荷的變化同樣會影響結(jié)構(gòu)的運(yùn)動。針對靶體結(jié)構(gòu)的特殊性以及流固耦合的復(fù)雜性, 采用流固耦合法和附加質(zhì)量法相結(jié)合的方式對水動作用力進(jìn)行求解。對靶體與其外表面周圍流體間的相互作用采用計算精度高的流固耦合法進(jìn)行求解; 對于靶體內(nèi)部濕式區(qū)域, 采用附加質(zhì)量法將流體介質(zhì)的動力影響以附加質(zhì)量的形式附加到靶體上, 即可保證碰撞區(qū)域的計算精度, 又可減少建模工作量, 提高計算效率。采用流固耦合法建模時, 結(jié)構(gòu)體采用Largrange單元模擬, 流體域采用Euler單元模擬, 并通過耦合控制方程將二者聯(lián)系, 使相互作用力通過耦合面發(fā)生傳遞與轉(zhuǎn)換, 所采用流固耦合控制方程[8]

        式中: 下標(biāo)和分別代表靶體及周圍流體; 為質(zhì)量矩陣, ;和分別為靶體位移向量和水壓; 為阻尼矩陣;為剛度矩陣, ;; 為壓力單元形函數(shù)矩陣;為結(jié)構(gòu)位移函數(shù); 為流體邊界法向量;是一個耦合矩陣, 表示流-固界面節(jié)點相互聯(lián)系的有效面積及界面上的法線節(jié)點力的矢量; 為結(jié)構(gòu)作用力; 為流體作用力。圖1為采用流固耦合法所建立的模型示意圖。

        附加質(zhì)量通過經(jīng)驗公式進(jìn)行計算, 根據(jù)相關(guān)文獻(xiàn)[9-11], 水下物體縱蕩和橫蕩運(yùn)動時的附加質(zhì)量分別為

        式中,為靶體質(zhì)量, 文中靶體縱蕩運(yùn)動時的附加質(zhì)量系數(shù)取0.04, 橫蕩運(yùn)動時的附加質(zhì)量系數(shù)取0.5。

        2 雷-靶碰撞數(shù)值仿真與分析

        2.1 計算模型

        2.1.1 目標(biāo)靶模型

        目標(biāo)靶為圓柱形回轉(zhuǎn)結(jié)構(gòu), 靶體板殼結(jié)構(gòu)采用殼單元進(jìn)行建模, 靶體內(nèi)部框架采用梁單元建模。靶體撞擊區(qū)域網(wǎng)格尺寸為60 mm, 非撞擊區(qū)域網(wǎng)格尺寸為150 mm。目標(biāo)靶體材料為雙線性各向同性硬化彈塑性材料模型。模型坐標(biāo)系以靶體長度方向為軸, 其中沿艏部方向為正; 靶體寬度方向為軸, 其中左舷為正方向; 型深方向為軸, 其中垂直于基線向上為正方向。靶體有限元模型見圖2。

        圖2 靶體有限元模型

        2.1.2 魚雷模型

        根據(jù)相關(guān)資料建立魚雷模型, 考慮到碰撞過程中只有雷頭及后續(xù)一部分結(jié)構(gòu)參與整個碰撞歷程, 因此在魚雷中后部建立實體, 通過調(diào)整密度使模型質(zhì)量與全雷實際質(zhì)量相同。魚雷艏端覆蓋有橡膠硫化層, 硫化層采用六面體實體單元模擬, 雷體金屬結(jié)構(gòu)采用六面體實體單元及六節(jié)點三棱柱單元模擬, 尾部采用六面體實體單元模擬, 魚雷有限元網(wǎng)格尺寸為30 mm。魚雷艏端硫化橡膠為超彈性材料模型, 中前段金屬為雙線性各向同性硬化彈塑性材料模型, 后段配重區(qū)域為剛體材料模型。魚雷有限元模型見圖3。

        圖3 魚雷有限元模型

        2.1.3 流場模型

        碰撞流體域采用圓柱區(qū)域, 直徑為目標(biāo)靶4倍直徑, 長30 m, 采用六面體網(wǎng)格劃分, 撞擊區(qū)域附近網(wǎng)格大小為60 mm, 非碰撞區(qū)域網(wǎng)格大小為100~150 mm。材料為NULL材料。魚雷、目標(biāo)靶及流場有限元網(wǎng)格總量約140萬。流場有限元模型見圖4, 碰撞系統(tǒng)有限元模型見圖5, 各區(qū)域材料參數(shù)見表1。

        圖4 流場有限元模型

        圖 5 碰撞有限元模型

        表1 仿真計算所用材料參數(shù)

        2.2 計算工況

        選取靶體2條環(huán)肋間隔最大處作為典型撞擊部位, 雷體分別以22 kn、30 kn和38 kn速度以30o、60o和90o角度撞擊靶體。

        2.3 計算結(jié)果分析

        通過分析碰撞過程中魚雷的速度時歷曲線以及加速度時歷曲線可得到碰撞過程中魚雷的運(yùn)動情況以及靶體各結(jié)構(gòu)件變形失效時序。

        2.3.1 靶體碰撞區(qū)域結(jié)構(gòu)損傷變形

        圖6反映了22 kn速度3種碰撞角度工況下靶體碰撞結(jié)構(gòu)損傷變形圖。從圖中可以看出, 結(jié)構(gòu)變形主要集中在撞擊過程中雷靶接觸區(qū)域, 距離碰撞區(qū)域較遠(yuǎn)的非碰撞區(qū)域并未出現(xiàn)形變, 類似剛體結(jié)構(gòu), 僅以質(zhì)量及慣性的形式對靶體運(yùn)動產(chǎn)生影響, 碰撞過程中靶體蒙皮變形程度隨碰撞角度的增加愈發(fā)明顯, 同時蒙皮損壞的形狀特征與雷頭的形狀有關(guān)。從90o碰撞變形圖可以發(fā)現(xiàn), 在此工況下, 靶體蒙皮區(qū)域仿真單元出現(xiàn)了失效, 表明碰撞區(qū)域蒙皮出現(xiàn)了撕裂及開口。

        圖6 22 kn速度下靶體蒙皮結(jié)構(gòu)損傷變形圖

        2.3.2 魚雷運(yùn)動

        圖7反映了碰撞過程中魚雷向速度分量的時歷曲線, 0時刻曲線與縱坐標(biāo)軸的交點表示碰撞開始瞬間魚雷的初速度分量。雷體材料與靶體材料相比剛度更大, 更不易變形, 僅在碰撞發(fā)生后短時間內(nèi)產(chǎn)生極小變形, 因此通過雷頭中后部監(jiān)測點速度變化來分析雷體運(yùn)動。從圖中可以看出在22 kn撞擊速度下, 雷頭接觸靶體蒙皮后蒙皮產(chǎn)生變形, 雷體速度急劇減小隨之衰減至0, 但由于靶體蒙皮及內(nèi)部框架的內(nèi)力釋放以及水動力的作用, 雷體從蒙皮表面彈開, 速度增加, 一定時間后, 由于慣性以及水動力的影響, 雷體速度趨于平穩(wěn), 在一定范圍內(nèi)震蕩。

        30 kn碰撞時, 30o和60o工況下魚雷速度變化與22 kn相同, 均是在碰撞前期迅速衰減后, 由于魚雷從蒙皮表面彈開, 速度反向增加。當(dāng)碰撞角度為90o時, 魚雷速度雖在碰撞前期迅速衰減, 但由于蒙皮強(qiáng)度達(dá)到極限, 發(fā)生失效破壞, 碰撞力卸載, 雷體穿透靶體向靶體內(nèi)部運(yùn)動, 在慣性及水動力作用下雷體速度震蕩并趨于穩(wěn)定。38 kn工況下, 30o碰撞時雷體同樣出現(xiàn)了反彈, 魚雷速度減小至0后反向增加。而60o及90o工況下雷體穿透蒙皮, 向靶體內(nèi)部運(yùn)動, 雷體速度在碰撞后期緩慢減小。

        圖7 不同工況下的魚雷速度時歷曲線

        對比不同撞擊角度結(jié)果發(fā)現(xiàn), 碰撞發(fā)生后魚雷速度迅速衰減, 隨著碰撞角度的增加, 雷體速度衰減會有一定程度的延后, 這是由于不同撞擊角度下, 雷頭觸及靶體后各結(jié)構(gòu)件變形程度不同, 從而導(dǎo)致魚雷速度衰減的程度不同。同時小角度撞擊下, 由于雷體與靶體蒙皮接觸面小, 雷體撞擊到蒙皮后產(chǎn)生反彈, 未穿透蒙皮, 隨撞擊角度的增加, 魚雷會穿透蒙皮繼續(xù)朝靶體內(nèi)部運(yùn)動。

        圖8反映了在22 kn速度90o工況時碰撞過程整個系統(tǒng)的能量變化。從圖中可以看出, 發(fā)生碰撞后, 系統(tǒng)動能轉(zhuǎn)化為內(nèi)能, 表現(xiàn)為動能曲線急劇減小, 內(nèi)能曲線增加, 當(dāng)碰撞結(jié)束后兩曲線趨于平穩(wěn), 極值點位置出現(xiàn)時間也與速度曲線相對應(yīng)。在碰撞過程中由于阻力、滑移及沙漏效應(yīng)等影響, 能量轉(zhuǎn)化過程出現(xiàn)損耗, 即為圖中其余能。

        圖8 能量變化時歷曲線

        圖9~圖 11反映了碰撞過程中雷頭監(jiān)測點部位的加速度時歷曲線, 由圖中曲線可知, 碰撞過程中加速度曲線呈現(xiàn)很強(qiáng)的非線性特征, 在不同的碰撞階段, 雷頭部位加速度出現(xiàn)不同程度的卸載。30°碰撞時, 碰撞前期雷頭加速度出現(xiàn)了明顯的大脈寬區(qū)域, 這說明在此工況下靶體蒙皮為明顯的塑性變形, 并且變形較慢。隨著碰撞速度的增大, 雷頭加速度在碰撞前期呈現(xiàn)出明顯的窄脈寬大幅值的卸載現(xiàn)象, 說明碰撞發(fā)生后蒙皮產(chǎn)生的塑性變形大、歷程短, 并伴有失效和破壞。碰撞結(jié)束后, 由于慣性和碰撞區(qū)域周圍流體的作用, 雷體加速度并未迅速衰減為零, 而是在一定幅值范圍內(nèi)震蕩, 并趨于平穩(wěn)。

        在碰撞發(fā)生開始后15 ms內(nèi), 雷頭監(jiān)測點部位加速度迅速增加, 出現(xiàn)較大幅值和脈寬區(qū)域, 該區(qū)域的幅值脈寬可包絡(luò)滿足魚雷觸發(fā)引信動作所要求的幅值及脈寬, 且具有較大余量。該結(jié)果說明在上述計算的工況范圍內(nèi), 目標(biāo)靶靶體剛度滿足魚雷觸發(fā)引信動作要求。同時, 隨著碰撞角度和速度的增加, 碰撞前期監(jiān)測點部位幅值及脈寬更大, 更能滿足魚雷觸發(fā)引信動作要求。

        圖9 22 kn工況下雷頭監(jiān)測點加速度時歷曲線

        圖10 30 kn工況下雷頭監(jiān)測點加速度時歷曲線

        圖11 38 kn工況下雷頭監(jiān)測點加速度時歷曲線

        3 結(jié)論

        文中針對魚雷與目標(biāo)靶碰撞問題, 基于流固耦合及附加質(zhì)量法, 利用LS-DYNA有限元分析軟件對不同撞擊速度及角度工況下的碰撞過程進(jìn)行了仿真, 得到了碰撞過程中靶體結(jié)構(gòu)損傷以及雷體運(yùn)動情況。

        1) 撞擊角度的變化對雷體運(yùn)動及靶體變形有一定程度的影響, 碰撞角度增大, 靶體蒙皮變形失效程度增加, 魚雷穿透蒙皮繼續(xù)向靶體內(nèi)部運(yùn)動的可能性增加。

        2) 魚雷撞擊速度越大, 雷頭加速度響應(yīng)峰值越高, 出現(xiàn)時間越早, 但其隨時間變化趨勢基本一致。且碰撞速度越大, 靶體損傷程度越大。

        [1] 肖文勇, 嚴(yán)謹(jǐn), 岳建軍, 等. 魚雷撞擊靶板動態(tài)仿真建模方法探討[J]. 艦船科學(xué)技術(shù), 2004, 26(3): 47-50.Xiao Wen-yong, Yan Jin, Yue Jian-jun, et al. Discussions on the Dynamic Simulation Modeling of a Torpedo Impacting the Plate[J]. Ship Science and Technology, 2004, 26(3): 47-50.

        [2] Minorsky. Analysis of Ship Collision Protection of Nuclear Powered Plant[J]. Journal of Ship Research, 1959, 3(2): 1-4.

        [3] Chang P Y, Seibold F, Thasanatorn C. A Rational Methodology for the Prediction of Structural Response due to Collisions of Ships[C]//Annual meeting of the Society of Naval Architects and Marine Engineers. New York, USA: Soc. Nav. Archit. Mar. Eng., Trans., 1980: 173-193.

        [4] 王自力, 顧永寧. 船舶碰撞動力學(xué)過程的數(shù)值仿真研究 [J]. 爆炸與沖擊, 2001, 21(1): 29-34.Wang Zi-li, Gu Yong-ning. Numerical Simulation of Ship Collisions[J]. Explosion and Shock Waves, 2001, 21(1): 29-34.

        [5] 朱新陽, 梅志遠(yuǎn), 吳梵. 潛艇典型結(jié)構(gòu)在撞擊載荷作用下動態(tài)響應(yīng)的試驗研究[J]. 船海工程, 2009, 38(4): 88-91. Zhu Xin-yang, Mei Zhi-yuan, Wu Fan. Research on Dynamic Response Test of Submarine Typical Structure Unit’s under Impact Load[J]. Ship & Ocean Engineering, 2009, 38(4): 88-91.

        [6] 楊桂通, 熊祝華. 塑性動力學(xué)[M]. 北京: 科學(xué)出版社, 1982.

        [7] 王自力. 船舶碰撞損傷機(jī)理與結(jié)構(gòu)耐撞性研究[D]. 上海: 上海交通大學(xué), 2000.

        [8] 汪玉, 周璞, 劉東岳. 考慮流固耦合作用的艦船抗沖擊仿真計算[J].振動與沖擊, 2005, 24(1): 74-80. Wang Yu, Zhou Pu, Liu Dong-yue. Numerical Simulation of Anti-shock Behavior of Ship Considering the Fluid-Structure Interaction in FEM[J]. Journal of Vibration and Shock, 2005, 24(1): 74-80.

        [9] Petersen M J. Dynamics of Ship Collisions[J]. Ocean Engineering, 1982, 9(4): 295-329.

        [10] 梅志遠(yuǎn). 基于MSC Dytran的潛艇結(jié)構(gòu)撞擊強(qiáng)度分析[J].計算機(jī)輔助工程, 2006, 15(z1): 71-74. Mei Zhi-yuan. Numerical Analysis Based on MSC Dytran Collision Strength of Submarine Structure[J]. Computer Aided Engineering, 2006, 15(z1): 71-74.

        [11] Chang P Y, Seibold F, Thasanatorn C. A Rational Methodology for the Prediction of Structural Response due to Collisions of Ships[J]. SNAME Transactions, 1980, 6(88): 173-193.

        Simulation on Structural Response of Torpedo-Target Collision

        LIJing-xiao, CHEN Yan-yong, ZHANG Tao, XU Da

        (Kunming Branch of the 705 Research Institute, China Shipbuilding Industry Corporation, Kunming 650118, China)

        In view of the particularity and complexity of target structure and collision environment in the collision process of a torpedo and a target, the internal mechanism of ship collision is utilized and the fluid-solid coupling and additional mass method are employed to carry out finite element simulation of the structural response of a torpedo when it hits a target. The structural damage of a target and the variation of torpedo velocity and acceleration during collision are studied. Simulation results show that: 1) The hitting process has strong nonlinear characteristics, and the hitting angle has a certain influence on the motion of a torpedo and the deformation of a target; 2) The higher the torpedo hitting velocity is, the higher the peak acceleration response of torpedo head and the damage degree of target become. This research may provide a reference for the structure design and dynamic characteristic design of targets.

        torpedo-target collision; fluid-solid coupling; additional mass method; finite element simulation; structural response

        李靜肖, 陳彥勇, 張濤, 等. 雷-靶碰撞結(jié)構(gòu)響應(yīng)仿真分析[J]. 水下無人系統(tǒng)學(xué)報, 2019, 27(1): 87-92.

        TJ630.3; TB122

        A

        2096-3920(2019)01-0087-06

        10.11993/j.issn.2096-3920.2019.01.015

        2018-09-20;

        2018-11-25.

        李靜肖(1994-), 男, 在讀碩士, 主要研究方向為魚雷總體技術(shù).

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

        猜你喜歡
        有限元結(jié)構(gòu)
        《形而上學(xué)》△卷的結(jié)構(gòu)和位置
        新型有機(jī)玻璃在站臺門的應(yīng)用及有限元分析
        基于有限元的深孔鏜削仿真及分析
        論結(jié)構(gòu)
        中華詩詞(2019年7期)2019-11-25 01:43:04
        基于有限元模型對踝模擬扭傷機(jī)制的探討
        新型平衡塊結(jié)構(gòu)的應(yīng)用
        模具制造(2019年3期)2019-06-06 02:10:54
        論《日出》的結(jié)構(gòu)
        創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長
        磨削淬硬殘余應(yīng)力的有限元分析
        基于SolidWorks的吸嘴支撐臂有限元分析
        天堂AV无码AV毛片毛| 亚洲色图第一页在线观看视频| 免费 无码 国产在线观看不卡| 久久久久国产精品片区无码| 日韩午夜三级在线视频| 国产精品一区二区三久久不卡| 天天狠天天添日日拍| 中文字幕亚洲乱码熟女一区二区| 久久精品re| 中文乱码字幕在线中文乱码| 人妖一区二区三区在线| 成人欧美日韩一区二区三区| 无码人妻丰满熟妇精品区| 一本无码人妻在中文字幕| 国产精品三级国产精品高| 国产免费人成视频在线 | 中文字幕无码日韩专区免费| 国内精品伊人久久久久av| 杨幂Av一区二区三区| 国产高清一区二区三区三州| 免费看又色又爽又黄的国产软件| 亚洲综合色区另类av| 青草热久精品视频在线观看| 中文字幕高清一区二区| 蜜桃av在线免费网站 | 亚洲国产中文在线二区三区免 | 国产香蕉尹人在线视频播放| 亚洲人妻av在线播放| 亚洲男女内射在线播放| 亚洲中文字幕无码久久| 国产欧美精品一区二区三区–老狼| 日本国主产一区二区三区在线观看| 国产日产韩国av在线| 无码人妻精品一区二区三| 久久久久久久99精品国产片| 亚洲av激情久久精品人| 婷婷开心五月亚洲综合| 老熟女富婆激情刺激对白| 国产真人无码作爱视频免费 | 亚洲尺码电影av久久| 天堂av在线一区二区|