• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于靜力試驗(yàn)的魚(yú)雷殼體模態(tài)計(jì)算模型修正

    2018-06-28 11:16:52尹韶平張志民
    關(guān)鍵詞:艙段魚(yú)雷靜力

    王 升, 尹韶平, 郭 君, 王 中, 張志民

    ?

    基于靜力試驗(yàn)的魚(yú)雷殼體模態(tài)計(jì)算模型修正

    王 升1,2, 尹韶平1, 郭 君1, 王 中1,2, 張志民1

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

    為建立一種簡(jiǎn)單而又適用的魚(yú)雷殼體楔環(huán)連接結(jié)構(gòu)模態(tài)計(jì)算分析模型, 將楔環(huán)連接的魚(yú)雷殼體簡(jiǎn)化為2個(gè)殼體以及1個(gè)環(huán)形圈連接的有限元模型, 其中包含7對(duì)接觸面以及1個(gè)接觸剛度參數(shù)。通過(guò)魚(yú)雷殼體的靜力加載試驗(yàn)以及靜力仿真對(duì)比, 辨識(shí)出模型中的接觸剛度參數(shù)。對(duì)修正后的有限元模型進(jìn)行模態(tài)仿真計(jì)算, 仿真模態(tài)頻率與試驗(yàn)?zāi)B(tài)頻率誤差不大于2.65%, 結(jié)果表明了該模態(tài)計(jì)算方法的準(zhǔn)確性。可為魚(yú)雷殼體等含楔環(huán)連接結(jié)構(gòu)的有限元建模與模態(tài)分析提供借鑒。

    魚(yú)雷殼體; 楔環(huán)連接結(jié)構(gòu); 模態(tài)仿真; 模型修正

    0 引言

    魚(yú)雷殼體聲輻射對(duì)魚(yú)雷的綜合作戰(zhàn)效能影響重大, 對(duì)魚(yú)雷的攻擊隱蔽性、線(xiàn)導(dǎo)導(dǎo)引有效性、自導(dǎo)性能及發(fā)射平臺(tái)的安全性均有重要影響[1]。魚(yú)雷殼體的聲輻射研究中最主要的一環(huán)是全雷模態(tài)計(jì)算, 在魚(yú)雷復(fù)雜的力學(xué)環(huán)境影響下, 全雷模態(tài)計(jì)算不準(zhǔn)確會(huì)導(dǎo)致不可意料的共振現(xiàn)象發(fā)生, 共振往往會(huì)引起較大振動(dòng)能級(jí)的全雷振動(dòng)噪聲, 甚至造成結(jié)構(gòu)破壞或者設(shè)備元件的損壞[2]。

    在對(duì)魚(yú)雷做模態(tài)仿真計(jì)算時(shí), 受楔環(huán)連接結(jié)構(gòu)復(fù)雜的靜動(dòng)力學(xué)特性影響, 不能準(zhǔn)確地建立有限元分析模型, 會(huì)使得仿真計(jì)算結(jié)果誤差較大, 這種誤差會(huì)造成很多不可避免的損失與危險(xiǎn)。如何方便、準(zhǔn)確地對(duì)魚(yú)雷殼體連接結(jié)構(gòu)進(jìn)行簡(jiǎn)化處理, 如何正確地設(shè)置有限元模型的參數(shù), 將直接影響到魚(yú)雷模態(tài)仿真的準(zhǔn)確性。

    國(guó)內(nèi)學(xué)者針對(duì)不同的研究?jī)?nèi)容, 采取的楔環(huán)連接結(jié)構(gòu)模型的簡(jiǎn)化方式以及楔環(huán)連接結(jié)構(gòu)的仿真計(jì)算方法大不相同。黃鵬等[3]在楔環(huán)連接結(jié)構(gòu)2種有限元優(yōu)化設(shè)計(jì)方案研究中, 建立楔環(huán)連接結(jié)構(gòu)的參數(shù)化有限元模型, 以結(jié)構(gòu)最大應(yīng)力比為目標(biāo)函數(shù), 對(duì)楔環(huán)徑向和軸向尺寸進(jìn)行優(yōu)化仿真設(shè)計(jì)。王路等[4]在workbench響應(yīng)譜分析的魚(yú)雷振動(dòng)傳遞優(yōu)化中, 基于結(jié)構(gòu)參數(shù)化等效模型, 進(jìn)行優(yōu)化計(jì)算, 降低結(jié)構(gòu)振動(dòng)傳遞響應(yīng)。馬銳磊等[5]在魚(yú)雷楔環(huán)連接結(jié)構(gòu)等效剛度建模與模態(tài)分析中, 采用了等效彈性模量法對(duì)魚(yú)雷楔環(huán)連接結(jié)構(gòu)進(jìn)行仿真計(jì)算。嚴(yán)海等[6]在模態(tài)參數(shù)的水下航行器楔環(huán)結(jié)構(gòu)有限元模型修正中, 基于固有頻率和模態(tài)振型等模態(tài)參數(shù)修正有限元模型。趙榮國(guó)[7]在楔環(huán)連接結(jié)構(gòu)靜力接觸行為與動(dòng)力學(xué)特性研究中, 將楔環(huán)連接結(jié)構(gòu)簡(jiǎn)化處理后, 調(diào)整接觸面積以及接觸罰剛度系數(shù)修正有限元模型。

    文中通過(guò)楔環(huán)連接殼體的靜力加載試驗(yàn)以及靜力仿真, 對(duì)楔環(huán)連接結(jié)構(gòu)有限元模型進(jìn)行修正, 采用修正后的有限元模型進(jìn)行魚(yú)雷殼體的模態(tài)仿真計(jì)算, 并與修正前模型計(jì)算結(jié)果進(jìn)行比較, 結(jié)果表明了該計(jì)算方法的準(zhǔn)確性。

    1 模態(tài)理論

    多自由度一般粘性阻尼結(jié)構(gòu)系統(tǒng)動(dòng)力學(xué)微分方程[8]為

    設(shè)式(2)的特解為

    將式(3)代入式(2), 得

    2 魚(yú)雷殼體建模

    2.1 楔環(huán)連接結(jié)構(gòu)模型簡(jiǎn)化

    取研究對(duì)象為楔環(huán)連接的2段魚(yú)雷殼體, 材料為鋁合金。其中, 楔環(huán)連接部分由上殼體、2個(gè)楔環(huán)帶、下殼體、密封圈、填片和蓋板等組成。楔環(huán)結(jié)構(gòu)的接觸面對(duì)較多, 有諸多肋骨、槽、孔、倒角和圓角等特征, 不便于有限元網(wǎng)格的劃分, 且計(jì)算耗時(shí), 容易造成仿真計(jì)算不收斂。因此, 有必要對(duì)楔環(huán)結(jié)構(gòu)的有限元模型進(jìn)行簡(jiǎn)化處理, 簡(jiǎn)化了如小孔、槽、倒角及圓角等一些不影響整體結(jié)構(gòu)響應(yīng)的特征。由于2個(gè)楔環(huán)帶滿(mǎn)足自鎖條件, 可視為一個(gè)整體, 且不關(guān)心2個(gè)楔環(huán)帶之間的受力狀況, 所以將2個(gè)楔環(huán)帶看做1個(gè)環(huán)形圈, 方便建模, 簡(jiǎn)化后的模型物理參數(shù)如表1所示。

    表1 模型物理參數(shù)

    由表1可以看出, 簡(jiǎn)化后的模型與實(shí)際結(jié)構(gòu)主要參數(shù)相差不大, 簡(jiǎn)化模型可用于仿真分析, 確保仿真結(jié)果具有較高的可信度。

    2.2 基于Hypermesh的有限元模型網(wǎng)格劃分

    采用UG建立魚(yú)雷殼體幾何模型后, 在Hypermesh中進(jìn)行網(wǎng)格劃分[9], 大部分采用hex六面體單元, 局部采用penta6五面體單元。共計(jì)48 860個(gè)單元, 其中楔環(huán)有6 000個(gè)網(wǎng)格, 艙段殼體1有30 480個(gè)網(wǎng)格, 艙段殼體2有12 380個(gè)網(wǎng)格, 網(wǎng)格質(zhì)量檢查良好。有限元模型如圖1所示, 局部接觸部分網(wǎng)格如圖2所示。

    2.3 有限元模型設(shè)置

    將網(wǎng)格模型導(dǎo)入Workbench中, 建立有限元仿真模型[10]。

    接觸面和目標(biāo)面的選擇一般有如下原則:

    1) 根據(jù)接觸面材料屬性選擇, 材料較軟的應(yīng)當(dāng)選為目標(biāo)面, 材料硬的選為接觸面;

    2) 接觸面的網(wǎng)格一般比目標(biāo)面網(wǎng)格細(xì)密;

    3) 劃分了低階單元的表面應(yīng)為目標(biāo)面;

    4) 表面較大的應(yīng)為目標(biāo)面;

    5) 應(yīng)該選取平面或凹面為目標(biāo)面。

    依據(jù)以上原則, 共設(shè)置了7對(duì)接觸面, 如圖2所示, 分別為楔環(huán)與魚(yú)雷艙段殼體1之間的2對(duì)接觸面, 即為圖中編號(hào)①、②部分, 其中楔環(huán)為接觸體, 艙段殼體1為目標(biāo)體; 楔環(huán)與魚(yú)雷艙段殼體2之間的2對(duì)接觸面, 即圖中編號(hào)③、④部分, 其中楔環(huán)為接觸體, 艙段殼體2為目標(biāo)體; 魚(yú)雷艙段殼體1與魚(yú)雷艙段殼體2之間的3對(duì)接觸面, 即圖中編號(hào)⑤、⑥、⑦部分, 其中艙段殼體1為接觸體, 艙段殼體2為目標(biāo)體。接觸面與目標(biāo)面采用綁定方式, 非對(duì)稱(chēng)接觸行為, 采用pure penalty接觸算法, 通過(guò)高斯點(diǎn)(Gauss point)探測(cè)接觸狀態(tài), 罰剛度系數(shù)為待修正參數(shù), 其余保留默認(rèn)設(shè)置, 至此完成魚(yú)雷連接殼體有限元仿真模型的建立。

    3 魚(yú)雷艙段殼體靜力加載試驗(yàn)

    3.1 靜力加載試驗(yàn)

    通過(guò)靜力加載試驗(yàn), 測(cè)試楔環(huán)連接結(jié)構(gòu)在外載下測(cè)點(diǎn)的位移量, 用于修正有限元仿真模型。

    如圖3所示, 本次試驗(yàn)試件為楔環(huán)連接的2段魚(yú)雷艙段殼體, 試驗(yàn)工裝有殼體緊固工裝、轉(zhuǎn)接工裝、位移補(bǔ)償工裝以及加載工裝, 其中位移補(bǔ)償工裝固定在緊固工裝上, 用于補(bǔ)償試驗(yàn)中工裝位移變形而引起的試件位移變形量。試驗(yàn)機(jī)為電子萬(wàn)能試驗(yàn)機(jī)。

    試件上的測(cè)點(diǎn)布局如圖4所示, 測(cè)點(diǎn)A、B、C位于試件側(cè)面, 距離右端面距離依次為100 mm, 200 mm, 400 mm。

    載荷加載方式為1 mm/min的等速位移加載, 每個(gè)載荷步保載30 s, 記錄每個(gè)載荷下所有測(cè)點(diǎn)的位移。載荷步如下: 600 N→1200 N→1800 N→2400 N→3000 N→3600 N→4200 N→4800 N→5400 N→6000 N→6600N→7200 N→7800 N→8400 N。

    3.2 靜力加載試驗(yàn)結(jié)果

    將萬(wàn)能實(shí)驗(yàn)機(jī)的端力值和位移值清零, 百分表調(diào)零后開(kāi)始試驗(yàn), 試驗(yàn)現(xiàn)場(chǎng)如圖5所示。

    測(cè)點(diǎn)A、B和C的位移測(cè)量值如表2所示。

    表2 測(cè)點(diǎn)位移值

    4 魚(yú)雷艙段殼體有限元模型修正

    4.1 靜力試驗(yàn)仿真

    采用Hypermesh建立包含工裝夾具的魚(yú)雷艙段殼體有限元模型, 其中魚(yú)雷艙段殼體模型的網(wǎng)格劃分與2.2節(jié)相同, 殼體模型與緊固工裝、加載工裝共劃分有65 687個(gè)單元。將模型導(dǎo)入workbench中進(jìn)行仿真計(jì)算, 魚(yú)雷殼體的有限元模型設(shè)置同2.3節(jié), 殼體與緊固工裝、殼體與加載工裝的接觸部分均采用綁定方式, 并保留默認(rèn)設(shè)置。楔環(huán)連接的2段殼體以及緊固工裝材料為鋁合金材料, 加載工裝材料為結(jié)構(gòu)鋼, 靜力仿真有限元模型如圖6所示。

    約束緊固工裝下端面的全部自由度, 在加載工裝上施加8400 N的力載荷, 方向豎直向下。計(jì)算結(jié)構(gòu)的整體變形。整體變形如圖7所示, 圖中(a)圖為前視圖; (b)圖為俯視圖; (c)圖為右視圖。由圖可知, 試件的右端發(fā)生橢圓形形變, 造成試件的底端發(fā)生局部向上彎曲。

    4.2 試驗(yàn)結(jié)果與分析

    采用測(cè)點(diǎn)A、B、C的數(shù)據(jù)修正楔環(huán)接觸結(jié)構(gòu)接觸面的罰剛度系數(shù), 進(jìn)行仿真計(jì)算。經(jīng)過(guò)多次試算, 取罰剛度系數(shù)=0.000 08時(shí), 測(cè)點(diǎn)A、B、C的仿真結(jié)果與試驗(yàn)結(jié)果的比較如表3所示, 與試驗(yàn)數(shù)據(jù)擬合較好, 擬合曲線(xiàn)如圖8所示。

    圖8可知, 當(dāng)載荷力大于3600 N(認(rèn)為此時(shí)魚(yú)雷殼體連接結(jié)構(gòu)接觸狀態(tài)趨于穩(wěn)定)時(shí), 誤差均小于20%, 測(cè)點(diǎn)A、B、C的仿真數(shù)據(jù)與試驗(yàn)數(shù)據(jù)擬合得很好, 即罰剛度系數(shù)=0.000 08時(shí)的面接觸剛度與楔環(huán)連接剛度接近, 可用于有限元模型的模態(tài)仿真計(jì)算。

    表3 仿真結(jié)果與試驗(yàn)結(jié)果比較

    5 魚(yú)雷艙段殼體模態(tài)計(jì)算

    初步仿真計(jì)算時(shí), 罰剛度系數(shù)值采用默認(rèn)設(shè)置, 即=1。計(jì)算分析得前6階模態(tài)頻率, 與試驗(yàn)?zāi)B(tài)頻率的比較如表4所示。

    表4 模型修正前仿真計(jì)算結(jié)果

    由表4可知, 仿真結(jié)果與試驗(yàn)結(jié)果誤差較大, 須對(duì)有限元模型進(jìn)行修正, 采用靜力仿真得到的罰剛度系數(shù)FKN=0.000 08再次進(jìn)行模態(tài)計(jì)算, 得到的結(jié)果與試驗(yàn)?zāi)B(tài)頻率的比較如表5所示, 仿真模態(tài)振型如圖9所示。由表5可知, 仿真結(jié)果與試驗(yàn)結(jié)果的誤差均不大于2.65%, 遠(yuǎn)好于修正前模型計(jì)算結(jié)果。由圖9可知, 仿真振型與試驗(yàn)振型一致, 表明基于靜力試驗(yàn)修正后的有限元模型能夠更準(zhǔn)確地對(duì)魚(yú)雷殼體進(jìn)行模態(tài)仿真計(jì)算。

    表5 模型修正后仿真計(jì)算結(jié)果

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

    文中建立了一種只有7對(duì)接觸面的魚(yú)雷艙段殼體有限元模型, 基于魚(yú)雷殼體靜力加載試驗(yàn)結(jié)果, 修正楔環(huán)連接的有限元模型, 并用于模態(tài)仿真計(jì)算。修正后模態(tài)頻率結(jié)果與試驗(yàn)?zāi)B(tài)頻率誤差不大于2.65%, 準(zhǔn)確度較高。為魚(yú)雷的模態(tài)仿真計(jì)算提供了一種方便、適用、準(zhǔn)確的模型簡(jiǎn)化方式以及有限元模型修正方法。

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

    [2] 尹韶平. 魚(yú)雷減振降噪技術(shù)[M]. 北京: 國(guó)防工業(yè)出版社, 2016.

    [3] 黃鵬, 尹益輝, 莫軍. 楔環(huán)連接結(jié)構(gòu)兩種有限元優(yōu)化設(shè)計(jì)方案研究[J]. 機(jī)械工程學(xué)報(bào), 2006, 42(8): 205-209.Huang Peng, Yin Yi-hui, Mo Jun. Two Optimization Designs Research of Wedge-Ring Joint Structure Based on FEM[J]. Chinese Journal of Mechanical Engineering, 2006, 42(8): 205-209.

    [4] 王路, 尹韶平, 曹小娟, 等. 基于Workbench響應(yīng)普分析的魚(yú)雷振動(dòng)傳遞優(yōu)化[J]. 魚(yú)雷技術(shù), 2016, 24(1): 13-17.Wang Lu, Yin Shao-ping, Cao Xiao-juan, et al. Vibration Transfer Optimization of Torpedo Based on Workbench Response Spectrum[J]. Torpedo Technology, 2016, 24(1): 13-17.

    [5] 馬銳磊, 尹韶平, 曹小娟, 等. 魚(yú)雷楔環(huán)連接結(jié)構(gòu)等效剛度建模與模態(tài)分析[J]. 艦船科學(xué)技術(shù), 2014, 36(2): 143-147.Ma Rui-lei, Yin Shao-ping, Cao Xiao-juan, et al. Research on Finite Element Modeling and Modal Analysis of Torpedo’s Wedged-Ring Connection Structure[J]. Ship Science and Technology, 2014, 36(2): 143-147.

    [6] 嚴(yán)海, 馬銳磊, 梁躍, 等. 基于模態(tài)參數(shù)的水下航行器楔環(huán)結(jié)構(gòu)有限元模型修正[J]. 魚(yú)雷技術(shù), 2016, 24(2): 87-93.Yan Hai, Ma Rui-lei, Liang Yue, et al. Finite Element Model Correction for Wedged-Ring Connection Structure of Underwater Vehicle Based on Modal Parameters[J]. Torpedo Technology. 2016, 24(2): 87-93.

    [7] 趙榮國(guó). 楔環(huán)連接結(jié)構(gòu)靜力接觸行為與動(dòng)力學(xué)特性研究[D]. 綿陽(yáng): 中國(guó)工程物理研究院, 2004.

    [8] 曹樹(shù)謙, 張文德, 蕭龍翔. 振動(dòng)結(jié)構(gòu)模態(tài)分析[M]. 天津: 天津出版社, 2014.

    [9] 彭妍, 婁一青, 汪亞超, 等. 基于HyperMesh軟件的復(fù)雜結(jié)構(gòu)有限元建模[J]. 水電能源科學(xué), 2010, 28(8): 100-105.Peng Yan, Lou Yi-qing, Wang Ya-chao, et al. Finite Element Model of Complex Structure Based on Hypermesh Software[J]. Water Resources and Power. 2010, 28(8): 100-105.

    [10] Gu R J, Murty P, Zheng Q. Use of Penalty Variable in Finite Element Analysis of Contacting Objects[J]. Compu- ters & Structures, 2002, 80(31): 2449-2459.

    (責(zé)任編輯: 陳 曦)

    Correction of Modal Calculation Model for Torpedo Shell Based on Static Test

    WANG Sheng1,2, YIN Shao-ping1, GUO Jun1, WANG Zhong1,2, ZHANG Zhi-min1

    (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)

    The torpedo shell joined with wedge-ring connection is simplified as two shells and one annular ring. This simple structure is modeled to build a simple, applicative finite element analysis model, which contains seven contact surfaces and a parameter of contact stiffness coefficient. Static loading test and static simulation analysis are conducted to identify the contact stiffness of the model. The modal of the modified model is calculated via simulation, and the relative error of the modal frequency between simulation and modal test is not more than 2.65%, verifying the accuracy of the proposed modal simulation method. This study may provide a reference for finite element modeling and modal simulation of the structures with wedge-ring connection, such as torpedo shell.

    torpedo shell; structure with wedge-ring connection; modal simulation; model correction

    王升,尹韶平,郭君,等. 基于靜力試驗(yàn)的魚(yú)雷殼體模態(tài)計(jì)算模型修正[J].水下無(wú)人系統(tǒng)學(xué)報(bào), 2018, 26(3): 247-252.

    TJ630.3; TP391.9

    A

    2096-3920(2018)03-0247-07

    10.11993/j.issn.2096-3920.2018.03.010

    2017-08-29;

    2017-11-02.

    王 升(1992-), 男, 在讀碩士, 主要研究方向?yàn)轸~(yú)雷總體技術(shù).

    猜你喜歡
    艙段魚(yú)雷靜力
    空間站艙段運(yùn)輸專(zhuān)列
    魚(yú)雷造反
    基于有限元仿真電機(jī)軸的靜力及疲勞分析
    軍事島 魚(yú)雷人
    帶孔懸臂梁靜力結(jié)構(gòu)的有限元分析
    昆鋼科技(2021年6期)2021-03-09 06:10:24
    基于ABAQUS的叉車(chē)轉(zhuǎn)向橋靜力分析
    魚(yú)雷也瘋狂
    基于TwinCAT的艙段數(shù)字化柔性自動(dòng)對(duì)接平臺(tái)控制系統(tǒng)設(shè)計(jì)
    水下航行器電池艙段溫度場(chǎng)數(shù)值模擬
    靜力性拉伸對(duì)少兒短距離自由泳打腿急效研究
    久久99蜜桃精品久久| 欧美日韩av久久| 亚洲精品乱码久久久v下载方式| 中文字幕人妻丝袜制服| 国产探花极品一区二区| 国产极品天堂在线| 最黄视频免费看| a级毛片免费高清观看在线播放| 国产老妇伦熟女老妇高清| 18禁动态无遮挡网站| 国产伦理片在线播放av一区| 国产精品成人在线| 黄片播放在线免费| 妹子高潮喷水视频| 免费看光身美女| 亚洲精品自拍成人| 伦理电影免费视频| 三级国产精品片| 两个人免费观看高清视频| 99国产综合亚洲精品| 成人国产麻豆网| 色视频在线一区二区三区| 亚洲精品第二区| 一区二区日韩欧美中文字幕 | 亚洲人成网站在线观看播放| 男女啪啪激烈高潮av片| 91久久精品国产一区二区三区| 女性被躁到高潮视频| 精品人妻一区二区三区麻豆| 看十八女毛片水多多多| 亚洲伊人久久精品综合| 亚洲婷婷狠狠爱综合网| 久久久a久久爽久久v久久| 精品酒店卫生间| 国产 一区精品| 18禁在线播放成人免费| 飞空精品影院首页| 国产白丝娇喘喷水9色精品| 欧美人与性动交α欧美精品济南到 | 一个人免费看片子| 亚洲国产精品一区二区三区在线| 亚洲综合色网址| 国产精品久久久久成人av| 亚洲国产精品国产精品| 国产熟女欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| av电影中文网址| 久久久久久久久久久丰满| 免费大片黄手机在线观看| 欧美精品一区二区免费开放| 亚洲精品av麻豆狂野| 妹子高潮喷水视频| 国产免费现黄频在线看| 亚州av有码| 美女cb高潮喷水在线观看| 乱码一卡2卡4卡精品| 免费高清在线观看日韩| 国产精品人妻久久久影院| 我要看黄色一级片免费的| videos熟女内射| 久久久久久久久久人人人人人人| 色婷婷av一区二区三区视频| 在现免费观看毛片| 热99久久久久精品小说推荐| 免费观看a级毛片全部| 观看美女的网站| 亚洲精品中文字幕在线视频| 亚洲中文av在线| 日韩精品有码人妻一区| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| 夫妻性生交免费视频一级片| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说| 男人添女人高潮全过程视频| 亚洲精品久久久久久婷婷小说| 搡老乐熟女国产| 99热这里只有是精品在线观看| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 欧美日韩视频高清一区二区三区二| 久久精品久久久久久久性| 人妻一区二区av| 欧美日韩视频精品一区| 成人无遮挡网站| videos熟女内射| 制服丝袜香蕉在线| 在线观看国产h片| 国产亚洲一区二区精品| 夫妻性生交免费视频一级片| 精品久久久久久久久av| 婷婷色av中文字幕| 国产av精品麻豆| 97超视频在线观看视频| 午夜激情福利司机影院| 国产欧美亚洲国产| 亚洲伊人久久精品综合| 亚洲欧美一区二区三区国产| 十八禁网站网址无遮挡| 一区二区av电影网| 中文乱码字字幕精品一区二区三区| 男女无遮挡免费网站观看| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 18禁动态无遮挡网站| 亚洲成色77777| 中国三级夫妇交换| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 桃花免费在线播放| 午夜激情福利司机影院| 看十八女毛片水多多多| 国产精品国产三级国产av玫瑰| 一级毛片我不卡| 国产成人精品一,二区| 人成视频在线观看免费观看| 久久久久久久精品精品| 国产亚洲精品第一综合不卡 | 免费不卡的大黄色大毛片视频在线观看| 一本一本综合久久| 亚洲在久久综合| 少妇丰满av| 日韩欧美精品免费久久| 久久久精品免费免费高清| 亚洲精品,欧美精品| xxxhd国产人妻xxx| 精品人妻一区二区三区麻豆| 我的老师免费观看完整版| 欧美少妇被猛烈插入视频| 久久久精品区二区三区| 欧美日韩综合久久久久久| 中文精品一卡2卡3卡4更新| 最近的中文字幕免费完整| 91精品三级在线观看| 久久精品国产自在天天线| 能在线免费看毛片的网站| 日本欧美国产在线视频| 热99久久久久精品小说推荐| 一区在线观看完整版| 日韩三级伦理在线观看| 国产老妇伦熟女老妇高清| 亚洲第一av免费看| 日本免费在线观看一区| 成人漫画全彩无遮挡| 纯流量卡能插随身wifi吗| 人体艺术视频欧美日本| 一级二级三级毛片免费看| 三上悠亚av全集在线观看| 九九在线视频观看精品| 国产日韩一区二区三区精品不卡 | 大香蕉久久成人网| 插阴视频在线观看视频| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 日本欧美国产在线视频| 久久久久久久久久成人| 五月玫瑰六月丁香| 美女视频免费永久观看网站| 人体艺术视频欧美日本| 综合色丁香网| 日韩精品免费视频一区二区三区 | 美女国产视频在线观看| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人 | 精品久久久久久久久av| 老熟女久久久| 久久久久久久久久久久大奶| 人人妻人人添人人爽欧美一区卜| 精品久久久噜噜| 这个男人来自地球电影免费观看 | 夜夜看夜夜爽夜夜摸| 男女边吃奶边做爰视频| 欧美 日韩 精品 国产| 熟女电影av网| 精品少妇久久久久久888优播| 十八禁高潮呻吟视频| 中文天堂在线官网| 高清不卡的av网站| 久热这里只有精品99| 国产熟女午夜一区二区三区 | 汤姆久久久久久久影院中文字幕| 性色avwww在线观看| 欧美变态另类bdsm刘玥| 草草在线视频免费看| 人妻少妇偷人精品九色| 女性生殖器流出的白浆| 色哟哟·www| 91精品国产九色| 永久网站在线| 91精品伊人久久大香线蕉| 精品一品国产午夜福利视频| 亚洲熟女精品中文字幕| 水蜜桃什么品种好| 99精国产麻豆久久婷婷| av国产精品久久久久影院| 在线观看一区二区三区激情| 69精品国产乱码久久久| 国产精品熟女久久久久浪| 人人妻人人爽人人添夜夜欢视频| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 哪个播放器可以免费观看大片| 热99国产精品久久久久久7| 欧美激情国产日韩精品一区| 一级毛片电影观看| 精品久久久噜噜| 制服诱惑二区| 亚洲人与动物交配视频| 一级毛片黄色毛片免费观看视频| av在线播放精品| 一级爰片在线观看| 精品国产一区二区久久| 婷婷成人精品国产| 日韩一本色道免费dvd| 在线观看免费高清a一片| 人妻 亚洲 视频| 最新中文字幕久久久久| 青春草国产在线视频| 日本免费在线观看一区| 嫩草影院入口| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 国产免费又黄又爽又色| 亚洲国产av影院在线观看| 狂野欧美白嫩少妇大欣赏| 国产精品.久久久| 欧美另类一区| 精品久久久噜噜| 欧美bdsm另类| 国产亚洲一区二区精品| 日韩av在线免费看完整版不卡| 卡戴珊不雅视频在线播放| 久久久精品免费免费高清| 亚洲伊人久久精品综合| 国产69精品久久久久777片| 欧美最新免费一区二区三区| 校园人妻丝袜中文字幕| 免费人成在线观看视频色| 最新中文字幕久久久久| 99热全是精品| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 岛国毛片在线播放| 美女xxoo啪啪120秒动态图| 在线观看一区二区三区激情| 极品人妻少妇av视频| 99热这里只有是精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 久久久精品94久久精品| 一级毛片aaaaaa免费看小| 视频区图区小说| 少妇 在线观看| 99精国产麻豆久久婷婷| 男男h啪啪无遮挡| 在线播放无遮挡| 97超视频在线观看视频| 十八禁高潮呻吟视频| 亚洲色图 男人天堂 中文字幕 | 欧美 日韩 精品 国产| 精品人妻偷拍中文字幕| a级毛片黄视频| 老司机影院毛片| 亚洲国产精品一区三区| 国产又色又爽无遮挡免| 欧美国产精品一级二级三级| 国产色婷婷99| 日韩大片免费观看网站| 高清在线视频一区二区三区| 国产一区二区在线观看日韩| 亚洲精品成人av观看孕妇| 免费观看在线日韩| 国产精品久久久久久av不卡| 看免费成人av毛片| 视频中文字幕在线观看| 在线播放无遮挡| 一边亲一边摸免费视频| 免费少妇av软件| 国产伦精品一区二区三区视频9| 日本黄色片子视频| 国产黄色免费在线视频| 日韩 亚洲 欧美在线| 免费高清在线观看日韩| 亚洲成人手机| 大香蕉久久网| 久久99热这里只频精品6学生| 两个人免费观看高清视频| 亚州av有码| 黑人巨大精品欧美一区二区蜜桃 | 日本av手机在线免费观看| 久久这里有精品视频免费| 亚洲av.av天堂| 久热久热在线精品观看| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频| 少妇被粗大的猛进出69影院 | 色视频在线一区二区三区| 亚洲人成网站在线观看播放| 一级毛片电影观看| 欧美变态另类bdsm刘玥| 欧美少妇被猛烈插入视频| 免费少妇av软件| 亚洲av.av天堂| 狠狠精品人妻久久久久久综合| 久久99精品国语久久久| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠久久av| av不卡在线播放| 免费av不卡在线播放| 天美传媒精品一区二区| 在线精品无人区一区二区三| 成人漫画全彩无遮挡| 午夜视频国产福利| 日本猛色少妇xxxxx猛交久久| 精品亚洲成a人片在线观看| 免费播放大片免费观看视频在线观看| 女的被弄到高潮叫床怎么办| 69精品国产乱码久久久| 亚洲av不卡在线观看| 日本与韩国留学比较| 在线免费观看不下载黄p国产| 高清黄色对白视频在线免费看| 亚洲综合精品二区| 精品国产露脸久久av麻豆| 在线观看免费高清a一片| 十八禁网站网址无遮挡| 亚洲色图综合在线观看| 精品酒店卫生间| 久久青草综合色| 国产毛片在线视频| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 9色porny在线观看| 欧美+日韩+精品| 成人漫画全彩无遮挡| 18+在线观看网站| 免费av中文字幕在线| av在线app专区| 亚洲精品自拍成人| 自拍欧美九色日韩亚洲蝌蚪91| 午夜av观看不卡| 毛片一级片免费看久久久久| 18禁观看日本| 欧美日韩综合久久久久久| 高清av免费在线| 一边摸一边做爽爽视频免费| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄| 性色avwww在线观看| 久久狼人影院| 日韩av在线免费看完整版不卡| 天堂8中文在线网| 熟妇人妻不卡中文字幕| 最近中文字幕2019免费版| 国产精品三级大全| 亚洲成人一二三区av| 人人妻人人澡人人爽人人夜夜| 成人国语在线视频| 亚洲五月色婷婷综合| 熟女av电影| 日韩不卡一区二区三区视频在线| 国内精品宾馆在线| 亚洲精华国产精华液的使用体验| 视频区图区小说| 国产精品久久久久成人av| 亚洲中文av在线| 国产综合精华液| 两个人免费观看高清视频| 777米奇影视久久| 亚洲少妇的诱惑av| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 亚洲美女黄色视频免费看| 一级黄片播放器| 母亲3免费完整高清在线观看 | 人妻一区二区av| 男女高潮啪啪啪动态图| 久久久精品免费免费高清| 午夜激情久久久久久久| 人妻 亚洲 视频| a级毛片免费高清观看在线播放| 久久久欧美国产精品| 国产亚洲最大av| 另类亚洲欧美激情| 国产精品一区二区在线不卡| 麻豆成人av视频| 国产成人av激情在线播放 | 多毛熟女@视频| 两个人免费观看高清视频| 国产片内射在线| 国产精品国产av在线观看| 国产高清国产精品国产三级| 伊人久久精品亚洲午夜| 18禁在线无遮挡免费观看视频| 久久女婷五月综合色啪小说| 性色avwww在线观看| 免费观看在线日韩| 国产黄色视频一区二区在线观看| 国产精品成人在线| 国产爽快片一区二区三区| 26uuu在线亚洲综合色| 成年av动漫网址| 色婷婷av一区二区三区视频| 黄色欧美视频在线观看| 成人影院久久| 久久久久国产网址| 成人无遮挡网站| 免费观看在线日韩| 久久精品夜色国产| 哪个播放器可以免费观看大片| 亚洲精品国产av蜜桃| www.色视频.com| 欧美精品高潮呻吟av久久| 中文精品一卡2卡3卡4更新| 国产色婷婷99| 黑人高潮一二区| 在线观看www视频免费| 亚洲国产精品一区二区三区在线| 免费大片黄手机在线观看| 日韩精品有码人妻一区| 久久精品久久久久久久性| 国产免费一级a男人的天堂| 丝袜美足系列| 欧美人与善性xxx| 精品午夜福利在线看| 三上悠亚av全集在线观看| 亚洲精品aⅴ在线观看| 啦啦啦视频在线资源免费观看| 男人添女人高潮全过程视频| 精品久久久久久久久亚洲| 午夜激情av网站| 人人妻人人澡人人爽人人夜夜| 人人澡人人妻人| 在线观看一区二区三区激情| av福利片在线| 秋霞在线观看毛片| 我的老师免费观看完整版| 美女内射精品一级片tv| 伦精品一区二区三区| 高清av免费在线| 亚洲图色成人| 欧美日韩成人在线一区二区| 精品国产国语对白av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 熟妇人妻不卡中文字幕| 爱豆传媒免费全集在线观看| 国产成人午夜福利电影在线观看| 亚洲精品av麻豆狂野| 美女脱内裤让男人舔精品视频| 国产精品嫩草影院av在线观看| 亚洲一区二区三区欧美精品| 久热久热在线精品观看| 亚洲一级一片aⅴ在线观看| 一级毛片aaaaaa免费看小| 丰满迷人的少妇在线观看| 免费大片黄手机在线观看| 一本久久精品| 母亲3免费完整高清在线观看 | 成人亚洲欧美一区二区av| 美女福利国产在线| 制服丝袜香蕉在线| 丰满饥渴人妻一区二区三| 男人操女人黄网站| 国产深夜福利视频在线观看| 精品午夜福利在线看| 亚洲精华国产精华液的使用体验| 午夜福利视频精品| 高清黄色对白视频在线免费看| 麻豆乱淫一区二区| 国产伦精品一区二区三区视频9| 黄色欧美视频在线观看| 天堂8中文在线网| 夫妻性生交免费视频一级片| 丝袜脚勾引网站| 欧美精品一区二区免费开放| 精品人妻一区二区三区麻豆| √禁漫天堂资源中文www| 亚洲性久久影院| 国产免费一区二区三区四区乱码| 伊人亚洲综合成人网| 欧美亚洲日本最大视频资源| 18禁裸乳无遮挡动漫免费视频| 午夜免费男女啪啪视频观看| 男男h啪啪无遮挡| 人妻夜夜爽99麻豆av| 日韩一区二区三区影片| 免费播放大片免费观看视频在线观看| 美女中出高潮动态图| 亚洲国产精品一区三区| 99精国产麻豆久久婷婷| 免费黄色在线免费观看| 精品国产一区二区久久| 日韩欧美一区视频在线观看| 久久99热这里只频精品6学生| 在线天堂最新版资源| 狂野欧美白嫩少妇大欣赏| av福利片在线| 亚洲三级黄色毛片| 久久精品久久精品一区二区三区| 欧美激情极品国产一区二区三区 | 一级黄片播放器| 看非洲黑人一级黄片| 亚洲第一av免费看| 人人妻人人添人人爽欧美一区卜| 久久女婷五月综合色啪小说| 99九九在线精品视频| 成人亚洲欧美一区二区av| 人妻少妇偷人精品九色| 美女主播在线视频| 自线自在国产av| 日韩不卡一区二区三区视频在线| 亚洲欧美色中文字幕在线| 欧美性感艳星| 精品国产露脸久久av麻豆| 99热6这里只有精品| 日本wwww免费看| 两个人免费观看高清视频| 日日摸夜夜添夜夜爱| 亚洲欧美精品自产自拍| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 国产成人一区二区在线| 一本大道久久a久久精品| 亚洲精品,欧美精品| 精品亚洲成a人片在线观看| 中文字幕最新亚洲高清| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 亚洲美女视频黄频| 国产精品国产av在线观看| 美女内射精品一级片tv| 国产黄色免费在线视频| 亚洲少妇的诱惑av| 丰满迷人的少妇在线观看| 最近中文字幕2019免费版| 人妻系列 视频| 桃花免费在线播放| 亚洲性久久影院| 我要看黄色一级片免费的| 国产午夜精品久久久久久一区二区三区| 国产亚洲av片在线观看秒播厂| av天堂久久9| 人体艺术视频欧美日本| 日本黄色片子视频| 国产精品人妻久久久影院| 免费看光身美女| 男人爽女人下面视频在线观看| 日本色播在线视频| 久久精品熟女亚洲av麻豆精品| 啦啦啦视频在线资源免费观看| 美女大奶头黄色视频| 狠狠精品人妻久久久久久综合| 一级a做视频免费观看| 久久久久久久久久久免费av| 免费黄色在线免费观看| 女人精品久久久久毛片| 午夜福利,免费看| 亚洲综合精品二区| 亚洲天堂av无毛| 精品视频人人做人人爽| 国产精品久久久久久av不卡| freevideosex欧美| 91成人精品电影| 观看美女的网站| 我的老师免费观看完整版| 18禁观看日本| 日韩强制内射视频| 99国产综合亚洲精品| 热re99久久精品国产66热6| 日韩三级伦理在线观看| 人妻制服诱惑在线中文字幕| 欧美日韩av久久| 精品久久久久久电影网| 人妻一区二区av| 午夜老司机福利剧场| 91aial.com中文字幕在线观看| 欧美97在线视频| 久久韩国三级中文字幕| 免费高清在线观看视频在线观看| 丝瓜视频免费看黄片| 国产欧美另类精品又又久久亚洲欧美| 久久鲁丝午夜福利片| 九色亚洲精品在线播放| 久久久久久伊人网av| 亚洲,欧美,日韩| 国产精品熟女久久久久浪| 国产极品粉嫩免费观看在线 | 啦啦啦在线观看免费高清www| 国产精品一区二区三区四区免费观看| 欧美精品亚洲一区二区| 日韩制服骚丝袜av| 久久av网站| 亚洲精品日本国产第一区| 一区二区av电影网| 国产毛片在线视频| 曰老女人黄片| 久久久久久久久大av| 国产成人freesex在线| 精品人妻熟女av久视频| 久久久久久久久大av| 国产成人freesex在线| 狂野欧美白嫩少妇大欣赏| 成人国语在线视频| 黄色一级大片看看| 国产日韩欧美视频二区| 亚洲国产欧美在线一区| 十八禁网站网址无遮挡| 51国产日韩欧美| 日日撸夜夜添| 一个人看视频在线观看www免费| 午夜激情av网站| 亚洲一级一片aⅴ在线观看| 大香蕉久久网| 丰满饥渴人妻一区二区三|