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

    基于仿真的動態(tài)載荷下圓柱支筒屈曲特性

    2016-06-06 01:06:36張金鵬王雨時李作華張志彪
    探測與控制學(xué)報 2016年2期
    關(guān)鍵詞:仿真

    張金鵬,王雨時,李作華,聞 泉,張志彪

    (1.南京理工大學(xué)機(jī)械工程學(xué)院,江蘇 南京 210094)(2.遼寧華興機(jī)電有限公司,遼寧 錦州 121017)

    ?

    基于仿真的動態(tài)載荷下圓柱支筒屈曲特性

    張金鵬1,王雨時1,李作華2,聞泉1,張志彪1

    (1.南京理工大學(xué)機(jī)械工程學(xué)院,江蘇 南京210094)(2.遼寧華興機(jī)電有限公司,遼寧 錦州121017)

    摘要:針對目前對圓柱殼狀零件屈曲特性進(jìn)行仿真分析時忽略其上蓋、圓角等結(jié)構(gòu)的問題,為引信圓柱支筒保險機(jī)構(gòu)優(yōu)化設(shè)計和性能分析提供參考,利用ANSYS Workbench軟件仿真研究了美軍M739A1引信中的圓柱支筒在動態(tài)載荷作用下的屈曲特性,得到彈丸壽命期內(nèi)可信環(huán)境下該圓柱支筒形態(tài)、軸向變形量以及最大應(yīng)力隨時間變化規(guī)律。震動環(huán)境、跌落過載低于47 000 g、發(fā)射過載低于38 000 g以及落速小于80 m/s觸碰土壤目標(biāo)時支筒均不會屈曲。落速大于100 m/s觸碰土壤目標(biāo)時支筒屈曲,且隨著落速提高,屈曲臨界時間明顯減小,但臨界應(yīng)力增大不明顯。彈丸對地觸發(fā)落角越小,屈曲臨界應(yīng)力越小,屈曲臨界時間越長。支筒屈曲特性受彈丸旋轉(zhuǎn)環(huán)境影響很小,而圓角和上蓋影響不可忽略。

    關(guān)鍵詞:引信;仿真;動態(tài)特性;圓柱支筒;屈曲;保險機(jī)構(gòu)

    0引言

    圓柱殼是一種在工程應(yīng)用中十分常見的結(jié)構(gòu)形式。M557、M739系列等中大口徑火炮榴彈用觸發(fā)引信應(yīng)用圓柱支筒零件作為剛性保險件。此類圓柱支筒零件失效分析時不僅要考慮強(qiáng)度失效,還要考慮屈曲失穩(wěn)。外部刺激使結(jié)構(gòu)內(nèi)部的應(yīng)力達(dá)到某一臨界值時,若有一微小增量,則結(jié)構(gòu)的平衡位形將發(fā)生很大改變,這種情況稱為結(jié)構(gòu)失穩(wěn)或屈曲[1],該臨界值稱為臨界應(yīng)力,而開始發(fā)生屈曲的時刻稱為臨界時間。屈曲是一種十分復(fù)雜的過程,它由材料性質(zhì)、結(jié)構(gòu)尺寸、加載方式和結(jié)構(gòu)缺陷等多種因素決定[2]。文獻(xiàn)[3]通過理論分析近似得到圓柱殼體發(fā)生第一次屈曲的時間和軸向內(nèi)力。文獻(xiàn)[4]根據(jù)B-R準(zhǔn)則對圓柱殼進(jìn)行非線性動力學(xué)響應(yīng)分析,利用LS-NYNA3D軟件對圓柱殼結(jié)構(gòu)在沖擊載荷下的動態(tài)屈曲進(jìn)行仿真。文獻(xiàn)[5]和[6]參考圓柱殼狀零件的屈曲特性,利用實(shí)驗(yàn)和仿真軟件,得到了截錐殼的靜態(tài)和動態(tài)屈曲特性。目前對圓柱殼狀零件屈曲特性進(jìn)行仿真和理論研究時多忽略上蓋、圓角等結(jié)構(gòu),把模型簡化成上無蓋、下無底的純圓柱殼,即圓筒。尚未見到有對帶上蓋、圓角等結(jié)構(gòu)的圓柱筒的屈曲特性進(jìn)行有限元分析。本文針對此問題,利用ANSYS Workbench仿真軟件,以M739A1引信支筒為例,應(yīng)用有限元方法分析在軸向動態(tài)載荷下支筒屈曲特性。

    1M739A1引信支筒簡介

    美軍M739A1引信瞬發(fā)觸發(fā)機(jī)構(gòu)如圖1所示,其中的支筒結(jié)構(gòu)如圖2所示,支筒幾何參數(shù)如表1所列[7]。支筒由柱殼、圓角和帶有中心孔的上蓋三部分組成,由紫銅帶T2沖制。彈丸觸碰目標(biāo)后,擊針在目標(biāo)反力作用下,向下壓圓柱殼狀支筒,戳擊雷管,激發(fā)傳爆序列。關(guān)于M739A1支筒的動態(tài)屈曲特性,以往多以實(shí)驗(yàn)研究為主。文獻(xiàn)[8]利用分離式霍普金森壓桿裝置,通過改變槍彈射速得到不同沖擊速度下支筒的屈曲臨界載荷和臨界時間。文獻(xiàn)[9]也利用分離式霍普金森壓桿裝置研究了引信瞬發(fā)觸發(fā)機(jī)構(gòu)在軸向沖擊載荷作用下動態(tài)塑性屈曲特性。

    圖1 M739A1引信瞬發(fā)觸發(fā)機(jī)構(gòu)Fig.1 Contact firing device of M739A1 fuze

    D/mmt/mmd/mmh/mmr/mmT/mmΦ5.99-0.050.252±0.025Φ2.09+0.127.87-0.120.49+0.060.51±0.033

    圖2 M739A1引信支筒Fig.2 Support of M739A1 fuze

    2有限元仿真模型建立

    2.1有限元建模和網(wǎng)格劃分

    據(jù)M739A1引信壽命期內(nèi)可信環(huán)境,從以下方面考慮支筒動態(tài)特性仿真用載荷:

    1)震動環(huán)境:震動試驗(yàn)采用峰值230±34.5g、持續(xù)時間2±0.2 ms的1/2正弦波脈沖[10]。仿真時,分頭向上、頭向下兩種方位,均對仿真模型施加峰值264.5g、持續(xù)時間2.2 ms、脈沖形狀1/2正弦波的載荷。

    2)跌落環(huán)境:彈丸跌落時,引信內(nèi)部零件受到的沖擊過載波形、量值和作用時間會因?yàn)閺椡韬鸵诺男螤?、結(jié)構(gòu)和材料等不同而不同,在設(shè)計初期一般估算無包裝裸態(tài)彈丸頭部向上從1.5 m高度垂直跌落到鋼板時引信內(nèi)部零件所受的沖擊過載最大值為7 000~15 000g,作用時間約為150~250 μs。仿真時對仿真模型施加峰值15 000g,持續(xù)時間250 μs、脈沖形狀1/2正弦波的載荷。

    3)發(fā)射環(huán)境:美軍M739A1引信主要配用于105 mm口徑和155 mm口徑榴彈炮。兩者發(fā)射加速度隨時間變化的曲線如圖3所示,其中105 mm口徑榴彈炮數(shù)據(jù)由內(nèi)彈道學(xué)解算獲得,155 mm榴彈炮數(shù)據(jù)取自文獻(xiàn)[11]。從圖3可看出,105 mm口徑榴彈炮發(fā)射加速度明顯大于155 mm口徑榴彈炮的,故仿真時對仿真模型施加圖3中105 mm口徑榴彈炮發(fā)射加速度。

    4)觸碰目標(biāo)環(huán)境:美軍M739A1引信以300 m/s落速進(jìn)入沙石土壤能可靠發(fā)火[7],故確定配用M739A1引信的某大口徑炮彈以300 m/s落速,30°、60°和90°落角觸碰土壤目標(biāo)來考核圓柱殼支筒的動態(tài)屈曲特性。

    圖3 發(fā)射加速度曲線Fig.3 Acceleration curve in setback environment

    為提高仿真效率,把彈丸簡化成防雨帽、防雨桿架、外殼、擊針頭、擊針、支筒、雷管座、引信體、彈體和炸藥10種模塊,如圖4所示。支筒有限元模型尺寸參數(shù)采用表1中的平均值。對有限元模型作如下假設(shè):

    圖4 彈丸和引信頭部物理模型Fig.4 Physical model of the projectile and fuze headand mesh

    ①材料連續(xù)、均勻、各向同性,且內(nèi)部無殘余應(yīng)力,材料表面光滑,無曲率突變;

    ②支筒、雷管座、擊針、擊針頭均軸對稱,忽略形位誤差和表面粗糙度影響;

    ③載荷作用過程無熱效應(yīng);

    ④不考慮自重影響;

    ⑤載荷作用前各模塊應(yīng)力為零。

    采用四面體單元,應(yīng)用映射面網(wǎng)格劃分法劃分網(wǎng)格,共劃分60 905個網(wǎng)格、248 336個節(jié)點(diǎn)。有限元模型如圖5所示。

    圖5 有限元模型及網(wǎng)格劃分Fig.5 Finite element model

    2.2定義材料

    支筒、防雨帽、防雨桿架和彈體應(yīng)用Johnson/Cook材料模型,其余模塊應(yīng)用PLASTIC材料模型。材料仿真參數(shù)如表2和表3所列,表中系平均值,未計其散布。

    表2 支筒、防雨帽、防雨桿架和彈體材料參數(shù)[12-14]

    表3 土壤、外殼、擊針頭、擊針等材料參數(shù)[12-14]

    3仿真結(jié)果及分析

    3.1各工況下支筒動態(tài)特性

    仿真分別得到震動環(huán)境、跌落環(huán)境、發(fā)射環(huán)境以及觸碰目標(biāo)環(huán)境下支筒最大應(yīng)力單元的應(yīng)力值和支筒軸向變形量隨時間變化的動態(tài)特性,如圖6所示。

    圖6 四種情況下支筒動態(tài)特性Fig.6 Dynamic characteristics of support under four cases

    由圖6可看出,在震動環(huán)境、跌落環(huán)境以及發(fā)射(后坐)環(huán)境下,雖然支筒所受載荷大小、載荷波形和作用時間不同,但在載荷作用時間內(nèi)支筒最大應(yīng)力增大到最大值后逐步減小過程中,支筒軸向變形量也隨之同步變化。說明這三種情況下,支筒并未發(fā)生屈曲,處于彈性變形階段。彈丸以300 m/s落速、90°落角觸碰土壤目標(biāo)后,一開始防雨帽和防雨桿架受碰撞反力作用,支筒并未受力,此段時間大約13 μs。隨后支筒受防雨桿架作用,內(nèi)部最大應(yīng)力先是很快上升,后又逐步變慢。在50~64 μs時刻,最大應(yīng)力穩(wěn)定在770~785 MPa,平均778 MPa。隨后,支筒最大應(yīng)力緩慢減小。從支筒受力開始,支筒軸向變形量逐漸增加,直到80 μs左右,支筒完全壓垮,殘余高度約3.01 mm。在50 μs時支筒發(fā)生屈曲,故支筒屈曲臨界時間約為50-13=37 μs,支筒屈曲臨界應(yīng)力為778 MPa。

    彈丸60°落角時,支筒從15 μs開始受防雨桿架作用,57 μs后支筒應(yīng)力穩(wěn)定在771 MPa,說明支筒屈曲臨界應(yīng)力為771 MPa,屈曲臨界時間約為57-15=42 μs。彈丸30°落角時,支筒從20 μs開始受防雨桿架作用,73 μs后支筒應(yīng)力穩(wěn)定在715 MPa,說明支筒屈曲臨界應(yīng)力為715 MPa,屈曲臨界時間約為73-20=53 μs。彈丸落角越小,支筒屈曲臨界應(yīng)力也越小,屈曲臨界時間越長。

    3.2觸碰目標(biāo)環(huán)境下支筒變形情況

    仿真得彈丸以300 m/s落速、90°落角觸碰土壤目標(biāo)時引信支筒變形情況如圖6所示。

    從圖6(e)和圖7可知彈丸以90°落角觸碰土壤目標(biāo)后,前13 μs為防雨帽與防雨桿架作用時間。隨后支筒在軸向力作用下,內(nèi)部應(yīng)力和軸向變形量逐漸變大,在50 μs時,支筒壓縮量達(dá)1.522 mm(擊針達(dá)到戳擊雷管位置,支筒壓縮量為1.52 mm),可認(rèn)為此段過程占引信瞬發(fā)度170 μs中的50 μs。因仿真模型里未加入雷管,故擊針繼續(xù)壓縮支筒。在被壓縮過程中首先是上部軸向壓縮,徑向變粗,然后是支筒下部軸向壓縮,徑向變粗。同樣可得到在60°落角下,彈丸從觸碰目標(biāo)到支筒壓縮至刺發(fā)雷管位置的時間占瞬發(fā)度中的56 μs,而30°落角時則為96 μs。

    圖7 彈丸以90°落角觸碰目標(biāo)時引信頭部和支筒動態(tài)特性Fig.7 Dynamic characteristics of the head of fuze and support when impacting the soil target

    3.3仿真可行性說明

    霍普金森桿實(shí)驗(yàn)得該支筒屈曲臨界載荷平均為4 029 N,屈曲臨界時間約為47.1 μs[2,8-9]。據(jù)此推斷屈曲臨界應(yīng)力為810 MPa;仿真得落角90°時支筒屈曲臨界應(yīng)力為778 MPa,屈曲臨界時間為50 μs,兩者基本接近。另外,上述仿真結(jié)果與M739A1系列引信支筒在各種可信環(huán)境下的性能狀態(tài)和水平未見有不一致之處。

    4極限情況

    4.1跌落和發(fā)射環(huán)境

    跌落環(huán)境下,載荷作用時間變化不大,但幅值會隨彈丸和引信形狀、結(jié)構(gòu)和材料等不同有很大變化。對仿真模型施加持續(xù)時間250 μs、脈沖形狀1/2正弦波但峰值不同的載荷,載荷峰值達(dá)47 000g時,支筒發(fā)生屈曲,屈曲臨界載荷為452 MPa。

    發(fā)射環(huán)境下,假設(shè)載荷作用時間和脈沖形狀不變,增加載荷峰值至38 000g時,支筒明顯屈曲,屈曲臨界載荷為427 MPa。

    4.2彈丸非旋觸碰目標(biāo)

    支筒取最大實(shí)體狀態(tài)、最小實(shí)體狀態(tài)和平均尺寸狀態(tài),不考慮材料彈塑性力學(xué)特性散布,彈丸以不同落速、90°落角觸碰土壤目標(biāo)時,支筒動態(tài)特性仿真結(jié)果如表4所列。

    表4 不同落速下支筒屈曲特性

    從表4中可看出,彈丸落速越大,則支筒臨界應(yīng)力越大、臨界時間越小。但支筒在平均尺寸下,彈丸落速從200 m/s增加到300 m/s,臨界應(yīng)力則從736 MPa增加到779 MPa,只增加6%,臨界時間卻由84 μs減小到37 μs,減少56%。可見隨落速提升,支筒臨界應(yīng)力增加不大,但臨界時間減小明顯。M739A1引信支筒在彈丸落速大于100 m/s時,支筒會發(fā)生明顯屈曲,而在彈丸落速低于80 m/s時,支筒不會發(fā)生明顯屈曲。

    4.3彈丸旋轉(zhuǎn)觸碰目標(biāo)

    文獻(xiàn)[7]給出M739A1引信最大轉(zhuǎn)速為500 r/s。保守分析假設(shè)該轉(zhuǎn)速不衰減。在彈丸以300 m/s落速和500 r/s轉(zhuǎn)速觸碰土壤目標(biāo)時,支筒從10 μs開始軸向受力和變形,46 μs后,支筒應(yīng)力穩(wěn)定在788 MPa,故支筒屈曲臨界應(yīng)力為788 MPa,屈曲臨界時間為46-10=36 μs。這與前文非旋轉(zhuǎn)條件300 m/s落速觸碰土壤的屈曲臨界載荷778 MPa和屈曲臨界時間37 μs非常接近,說明彈丸觸碰目標(biāo)時,轉(zhuǎn)速對支筒的屈曲特性影響可忽略不計。

    4.4圓柱殼支筒是否帶上蓋和圓角的差別

    仿真得到純圓柱殼支筒以及帶上蓋和圓角支筒在上述極限發(fā)射環(huán)境下,最大應(yīng)力單元應(yīng)力分別是166 MPa和427 MPa,未帶上蓋的圓柱支筒沒有發(fā)生屈曲;而在以300 m/s落速、90°落角觸碰土壤目標(biāo)時,兩者屈曲臨界應(yīng)力分別是570 MPa和778 MPa,屈曲臨界時間分別為45 μs和37 μs??梢妶A角和上蓋對支筒的屈曲特性有較大影響,純圓柱殼支筒的承載能力更高。故若應(yīng)用純圓柱殼模型對M739A1引信支筒進(jìn)行簡化分析,會存在較大誤差。

    5結(jié)論

    1)震動環(huán)境下,支筒處于彈性變形區(qū)域,不會發(fā)生屈曲。

    2)跌落環(huán)境下,過載峰值高于47 000g時,平均尺寸狀態(tài)下支筒才會發(fā)生屈曲。

    3)發(fā)射環(huán)境下,過載峰值高于38 000g時,平均尺寸狀態(tài)下支筒才會發(fā)生屈曲。

    4)觸碰土壤目標(biāo)情況下,彈丸以300 m/s落速,90°落角觸碰土壤目標(biāo)時,支筒屈曲臨界應(yīng)力為778 MPa,屈曲臨界時間為37 μs。從彈丸著地到擊針壓潰支筒,刺發(fā)雷管,耗時約50 μs,即此段占引信瞬發(fā)度170 μs中的50 μs。

    5)彈丸以300 m/s落速觸碰土壤目標(biāo)時,落角越小,支筒屈曲臨界應(yīng)力越小,屈曲臨界時間越長。在未考慮材料彈塑性力學(xué)特性散布的前提下,落速小于80 m/s,支筒不會發(fā)生明顯屈曲,落速大于100 m/s,支筒會發(fā)生明顯屈曲,并且彈丸落速越大,引信支筒屈曲臨界時間越短,臨界應(yīng)力越大,但增大不明顯。

    6)彈丸旋轉(zhuǎn)對支筒的屈曲特性影響不大,故M739A1引信用于微旋彈時,觸發(fā)靈敏度和瞬發(fā)度不會有明顯的減退。

    7)不可忽略圓角和上蓋對支筒的屈曲特性的影響。

    8)M739引信頭部圓柱殼支筒保險機(jī)構(gòu)保險與“解除保險”設(shè)計閾值散布較小,裕度也較大,故其保險和觸發(fā)作用的可靠性均很高。

    參考文獻(xiàn):

    [1]陳鐵云,沈惠申.結(jié)構(gòu)的屈曲[M].上海:上??茖W(xué)技術(shù)文獻(xiàn)出版社,1993:21-64.

    [2]王國泰,譚惠民.引信支筒動態(tài)屈曲的實(shí)驗(yàn)研究[J].探測與控制學(xué)報(原《現(xiàn)代引信》),1988(1):29-33.

    [3]譚惠民,高世橋.圓柱殼受軸向沖擊作用時的塑性動態(tài)屈曲研究[J].北京:北京理工大學(xué)學(xué)報,1990(3):98-102.

    [4]李帥.沖擊載荷下圓柱殼非線性動力屈曲的數(shù)值研究[D].大連:大連理工大學(xué),2005.

    [5]劉茂全.引信后坐保險用截頂圓錐殼支筒力學(xué)特性[D].南京:南京理工大學(xué),2006.

    [6]劉茂全,王雨時.引信后坐保險用截頂圓錐殼形支筒屈曲特性研究[J].探測與控制學(xué)報,2007(s1):50-58.

    [7]MIL-HDBK-145C.MILITARY HANDBOOK ACTIVE FUZE CATALOG[M].USA:Department of Defense,2000

    [8]馬寶華,譚惠民,石巖,等.美國M739A1和M577A1引信反求工程研究報告[R].北京:北京理工大學(xué),1989.

    [9]王國泰,譚惠民,楊欣德.柱殼在軸向沖擊載荷作用下動態(tài)塑性屈曲的實(shí)驗(yàn)研究[J].應(yīng)用力學(xué)學(xué)報,1994(1):105-108.

    [10]GJB573A—引信環(huán)境與性能試驗(yàn)方法[S].北京:中國人民解放軍總裝備部,1998.

    [11]火炮技術(shù)研究開發(fā)中心.155mm自行加榴炮武器系統(tǒng)論文集[M].北京:航空工業(yè)出版社,2013:388-395.

    [12]《工程材料實(shí)用手冊》編寫委員會.工程材料實(shí)用手冊[M].北京:中國標(biāo)準(zhǔn)出版社,2002.

    [13]楊翔,王雨時,聞泉,等.單兵火箭彈引信前沖過載與觸發(fā)開關(guān)動態(tài)特性[J].探測與控制學(xué)報,2015(2):46-50.

    [14]陳勇軍,王雨時,聞泉.炮彈穿甲爆破彈彈底引信前沖過載系數(shù)仿真方法[J].探測與控制學(xué)報,2015(2):40-45.

    Buckling Characteristics Simulation of Cylindrical Support under Dynamic Load

    ZHANG Jinpeng1,WANG Yushi1,LI zuohua2,WEN Quan1,ZHANG Zhibiao1

    (1.School of Mechanical Engineering, Nanjing University of Science & Technology, Nanjing 210094, China;2.Liaoning Huaxing Electromechanical Co.Ltd,Jinzhou 121017,China)

    Abstract:Due to the problem that ignoring the structure of loam cake and rounded when doing the simulation of cylindrical shell buckling characteristics and in order to provides reference for optimizing the design and performance analysis of cylindrical support safe device, the ANSYS Workbench was used to do the study of simulation that the buckling characteristics of cylindrical support of M739 fuze under dynamic load. Getting the rule which cylindrical support’s form, axial deformation and maximum stress changed with time under the trusted environment of projectile lifetime. The support could not buckle under the vibration environment, less than 47 000 g fall overload, less than 38 000 g setback overload and impact the soil with the velocity less than 80 m/s. While the velocity was more than 100 m/s, the support buckle. The buckling critical time increased obviously with the velocity increasing, while the critical stress increases unobviously. The falling angle is smaller, the critical stress is less and critical time is longer. The influence of projectile spinning to support’s bucking characteristics could be ignored, while the fillet and loam cake cannot.

    Key words:fuze; simulation ; dynamic characteristics; cylindrical support, buckling; safety feature

    中圖分類號:TJ430

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1008-1194(2016)02-0037-06

    作者簡介:張金鵬(1990—),男,江蘇揚(yáng)州人,碩士研究生,研究方向:精密機(jī)械設(shè)計。E-mail:825163825@qq.com。

    基金項(xiàng)目:武器裝備預(yù)先研究項(xiàng)目資助(51305040101)

    *收稿日期:2015-09-17

    猜你喜歡
    仿真
    Proteus仿真軟件在單片機(jī)原理及應(yīng)用課程教學(xué)中的應(yīng)用
    工業(yè)機(jī)器人模擬仿真技術(shù)在職業(yè)教育中的應(yīng)用淺析
    一種幫助幼兒車內(nèi)脫險應(yīng)急裝置的仿真分析
    科技資訊(2016年18期)2016-11-15 20:09:22
    論虛擬仿真實(shí)訓(xùn)系統(tǒng)在口腔實(shí)驗(yàn)教學(xué)中的應(yīng)用
    科技資訊(2016年18期)2016-11-15 07:55:28
    基于機(jī)電設(shè)備電氣控制線路排故的仿真系統(tǒng)設(shè)計
    價值工程(2016年29期)2016-11-14 02:01:16
    航空電氣系統(tǒng)中故障電弧研究
    Buck開關(guān)變換器的基本參數(shù)設(shè)計及仿真分析
    試析PLC控制下的自動化立體倉庫仿真情況分析
    基于MADYMO的航空座椅約束系統(tǒng)優(yōu)化設(shè)計
    科技視界(2016年18期)2016-11-03 21:44:44
    中國體態(tài)假人模型與FAA Hybrid Ⅲ 型假人模型沖擊差異性分析
    科技視界(2016年18期)2016-11-03 20:31:49
    日日撸夜夜添| 久久午夜福利片| 国产精品一区二区在线不卡| 大片电影免费在线观看免费| 七月丁香在线播放| 亚洲综合色惰| 26uuu在线亚洲综合色| av国产免费在线观看| 婷婷色麻豆天堂久久| 日韩欧美 国产精品| 亚洲激情五月婷婷啪啪| 精品国产乱码久久久久久小说| 亚洲av成人精品一二三区| 国产成人91sexporn| 亚洲欧美成人综合另类久久久| 精品一区二区免费观看| 成人国产麻豆网| 99久久精品热视频| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合大香蕉| 精华霜和精华液先用哪个| 少妇丰满av| 国产一区二区三区av在线| 免费观看a级毛片全部| 国内揄拍国产精品人妻在线| 亚洲av男天堂| 一级片'在线观看视频| 麻豆国产97在线/欧美| 精品一区二区免费观看| 最近最新中文字幕免费大全7| 美女国产视频在线观看| 777米奇影视久久| 成人免费观看视频高清| 99视频精品全部免费 在线| 精品国产露脸久久av麻豆| 九九在线视频观看精品| 亚洲欧美一区二区三区黑人 | 美女内射精品一级片tv| 欧美成人一区二区免费高清观看| 国产免费一区二区三区四区乱码| 精品久久久久久久末码| 久久国内精品自在自线图片| 亚洲人与动物交配视频| 欧美丝袜亚洲另类| 18禁动态无遮挡网站| 中文字幕制服av| 亚洲人成网站在线播| 国产久久久一区二区三区| 国产乱人偷精品视频| 国产精品麻豆人妻色哟哟久久| 久久国产精品男人的天堂亚洲 | 欧美高清性xxxxhd video| 亚洲av在线观看美女高潮| 99热这里只有是精品50| 涩涩av久久男人的天堂| 中文在线观看免费www的网站| 成人一区二区视频在线观看| 一级毛片aaaaaa免费看小| 91精品一卡2卡3卡4卡| 久久av网站| 伦理电影大哥的女人| 国产熟女欧美一区二区| 男人舔奶头视频| 国产一区二区三区av在线| 国产伦精品一区二区三区视频9| 最近最新中文字幕免费大全7| 亚洲av在线观看美女高潮| 在线天堂最新版资源| 午夜福利在线观看免费完整高清在| 亚洲精品aⅴ在线观看| 只有这里有精品99| av免费观看日本| 一本—道久久a久久精品蜜桃钙片| 久久久国产一区二区| 一级二级三级毛片免费看| 精品久久久久久久末码| 乱系列少妇在线播放| 久热久热在线精品观看| 亚洲精品一区蜜桃| 成年女人在线观看亚洲视频| 亚洲性久久影院| 色网站视频免费| 国产黄色免费在线视频| 又大又黄又爽视频免费| 国产精品免费大片| 久久久久久久精品精品| 一本久久精品| 久久精品国产亚洲网站| 亚洲美女搞黄在线观看| 欧美一级a爱片免费观看看| 欧美日本视频| 亚洲av二区三区四区| 国产永久视频网站| 亚洲精品日韩av片在线观看| 高清日韩中文字幕在线| 国产精品精品国产色婷婷| 99久久综合免费| 免费在线观看成人毛片| 涩涩av久久男人的天堂| 亚洲精品一二三| 色综合色国产| 国产精品一区二区性色av| 国产真实伦视频高清在线观看| 久久久久久伊人网av| 亚洲精品国产成人久久av| 亚洲精品国产av成人精品| 久久99热6这里只有精品| 中文在线观看免费www的网站| 97超碰精品成人国产| 免费观看av网站的网址| 日产精品乱码卡一卡2卡三| 亚洲av成人精品一二三区| 色哟哟·www| 能在线免费看毛片的网站| 一级黄片播放器| 国产在线一区二区三区精| 亚洲精品一二三| 国产精品国产三级专区第一集| 婷婷色麻豆天堂久久| 一本色道久久久久久精品综合| 免费大片18禁| 亚洲欧美清纯卡通| 麻豆乱淫一区二区| 国产大屁股一区二区在线视频| 麻豆精品久久久久久蜜桃| 国产成人欧美| 男人添女人高潮全过程视频| 国产亚洲欧美在线一区二区| 一级,二级,三级黄色视频| 久久久久视频综合| 国精品久久久久久国模美| 欧美精品亚洲一区二区| 新久久久久国产一级毛片| 男女床上黄色一级片免费看| 国产91精品成人一区二区三区 | 国产激情久久老熟女| 高清不卡的av网站| 一级a爱视频在线免费观看| 乱人伦中国视频| 热99久久久久精品小说推荐| 成人亚洲精品一区在线观看| 丝袜人妻中文字幕| 色综合欧美亚洲国产小说| 高清不卡的av网站| 又紧又爽又黄一区二区| 97在线人人人人妻| 亚洲久久久国产精品| 亚洲黑人精品在线| 在线av久久热| 校园人妻丝袜中文字幕| 曰老女人黄片| 亚洲少妇的诱惑av| 欧美老熟妇乱子伦牲交| 日本a在线网址| 人人妻人人澡人人看| 精品一区在线观看国产| 日韩免费高清中文字幕av| 久久久久久免费高清国产稀缺| 日韩电影二区| 成在线人永久免费视频| 国产精品一国产av| 午夜久久久在线观看| 欧美乱码精品一区二区三区| 超碰97精品在线观看| 美女午夜性视频免费| 女性生殖器流出的白浆| 男的添女的下面高潮视频| 视频区图区小说| 久久影院123| 亚洲午夜精品一区,二区,三区| 国产深夜福利视频在线观看| 在线天堂中文资源库| 电影成人av| 一级片'在线观看视频| 亚洲国产欧美网| 亚洲三区欧美一区| 一区福利在线观看| 久久中文字幕一级| 男人操女人黄网站| 成人影院久久| 不卡av一区二区三区| 欧美日韩亚洲高清精品| 国产精品国产av在线观看| 大型av网站在线播放| 两人在一起打扑克的视频| av国产精品久久久久影院| 久久久久国产精品人妻一区二区| 男女午夜视频在线观看| 国产一区二区三区综合在线观看| 国产女主播在线喷水免费视频网站| 国产一区二区在线观看av| 黄色a级毛片大全视频| 国产日韩欧美视频二区| 青春草视频在线免费观看| 亚洲av电影在线进入| 欧美成狂野欧美在线观看| 欧美另类一区| 2021少妇久久久久久久久久久| 91成人精品电影| 久久精品人人爽人人爽视色| 久久性视频一级片| 精品一区在线观看国产| 丝袜脚勾引网站| tube8黄色片| 夫妻性生交免费视频一级片| 日本av免费视频播放| 久久久久精品国产欧美久久久 | 天堂8中文在线网| 狂野欧美激情性bbbbbb| 国产一区有黄有色的免费视频| 美女脱内裤让男人舔精品视频| 国产欧美日韩综合在线一区二区| 亚洲三区欧美一区| 精品一区在线观看国产| 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 久久鲁丝午夜福利片| 久久精品亚洲av国产电影网| 久久久国产一区二区| 午夜91福利影院| 亚洲国产成人一精品久久久| 一二三四社区在线视频社区8| 日韩av在线免费看完整版不卡| 中文精品一卡2卡3卡4更新| 欧美性长视频在线观看| 久久久久久久国产电影| 国产欧美亚洲国产| 成年人黄色毛片网站| 久久人人爽av亚洲精品天堂| 久久人人爽人人片av| 精品人妻在线不人妻| 纵有疾风起免费观看全集完整版| 成人国产av品久久久| 人妻 亚洲 视频| 大型av网站在线播放| 成人三级做爰电影| 日韩免费高清中文字幕av| 欧美日韩亚洲高清精品| 蜜桃国产av成人99| 午夜福利,免费看| 免费人妻精品一区二区三区视频| www.熟女人妻精品国产| 亚洲精品第二区| 午夜福利一区二区在线看| 亚洲国产欧美一区二区综合| 国产男女内射视频| 亚洲午夜精品一区,二区,三区| 久久狼人影院| av网站免费在线观看视频| 少妇的丰满在线观看| 欧美人与性动交α欧美软件| 国产免费福利视频在线观看| 女人高潮潮喷娇喘18禁视频| 捣出白浆h1v1| 后天国语完整版免费观看| 国产99久久九九免费精品| 国产精品一二三区在线看| 熟女av电影| 亚洲人成电影观看| 精品第一国产精品| 侵犯人妻中文字幕一二三四区| 午夜影院在线不卡| 欧美激情高清一区二区三区| 日本wwww免费看| 一本久久精品| 欧美在线黄色| 国产成人精品久久二区二区91| 婷婷丁香在线五月| 欧美精品啪啪一区二区三区 | 老司机深夜福利视频在线观看 | 久久精品熟女亚洲av麻豆精品| 国产免费福利视频在线观看| 亚洲精品国产一区二区精华液| 一边摸一边抽搐一进一出视频| 青春草视频在线免费观看| 久久久久久久久久久久大奶| 欧美另类一区| 激情视频va一区二区三区| 男女床上黄色一级片免费看| 午夜av观看不卡| 一区二区三区四区激情视频| 91字幕亚洲| 国产高清不卡午夜福利| 91麻豆av在线| 亚洲国产精品成人久久小说| 国产av一区二区精品久久| 亚洲精品日韩在线中文字幕| 久久久久精品国产欧美久久久 | 高清不卡的av网站| 中文字幕制服av| 啦啦啦 在线观看视频| 亚洲欧洲国产日韩| 午夜91福利影院| 国产免费又黄又爽又色| 欧美精品一区二区大全| 国产高清国产精品国产三级| 久久久精品区二区三区| 黄色视频不卡| 日韩制服骚丝袜av| 亚洲国产最新在线播放| 国产在线观看jvid| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美中文字幕日韩二区| 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 国产人伦9x9x在线观看| 蜜桃国产av成人99| 亚洲av日韩精品久久久久久密 | 精品卡一卡二卡四卡免费| av在线老鸭窝| 亚洲一码二码三码区别大吗| 男女国产视频网站| 午夜两性在线视频| 五月开心婷婷网| 一本—道久久a久久精品蜜桃钙片| 欧美国产精品va在线观看不卡| 亚洲精品久久成人aⅴ小说| 大香蕉久久成人网| 国产av国产精品国产| 成人国产一区最新在线观看 | 亚洲成人手机| 在线观看免费视频网站a站| 制服诱惑二区| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 亚洲av成人不卡在线观看播放网 | 人体艺术视频欧美日本| 午夜久久久在线观看| 一区二区三区激情视频| 黄色一级大片看看| 一级毛片 在线播放| 水蜜桃什么品种好| 性高湖久久久久久久久免费观看| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡| 亚洲国产成人一精品久久久| 大陆偷拍与自拍| 多毛熟女@视频| 国产高清视频在线播放一区 | 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 亚洲欧美激情在线| 久久精品熟女亚洲av麻豆精品| 搡老岳熟女国产| 热99久久久久精品小说推荐| 久热这里只有精品99| 欧美黑人精品巨大| 亚洲第一青青草原| 熟女av电影| 国产一区二区激情短视频 | 欧美中文综合在线视频| 丝袜脚勾引网站| 欧美黄色片欧美黄色片| 女警被强在线播放| 老司机影院成人| 黄片小视频在线播放| 国产欧美日韩综合在线一区二区| 天天影视国产精品| 精品国产超薄肉色丝袜足j| 天天添夜夜摸| 高清av免费在线| 亚洲成人免费电影在线观看 | 久久人人爽av亚洲精品天堂| 午夜精品国产一区二区电影| 狠狠婷婷综合久久久久久88av| 在线观看国产h片| 乱人伦中国视频| 亚洲精品自拍成人| 精品高清国产在线一区| 成年动漫av网址| 精品国产乱码久久久久久小说| 国产精品熟女久久久久浪| 亚洲成av片中文字幕在线观看| 日韩中文字幕视频在线看片| 亚洲三区欧美一区| 一级毛片黄色毛片免费观看视频| 国产精品一区二区在线不卡| 国产高清视频在线播放一区 | 一级,二级,三级黄色视频| 91精品国产国语对白视频| 亚洲精品久久成人aⅴ小说| 日本午夜av视频| 五月天丁香电影| 五月开心婷婷网| 国产激情久久老熟女| 午夜精品国产一区二区电影| 中文字幕人妻丝袜制服| 久久热在线av| 少妇精品久久久久久久| 成年人黄色毛片网站| 亚洲国产精品成人久久小说| 女性生殖器流出的白浆| 2018国产大陆天天弄谢| 久久久久久久国产电影| 精品福利永久在线观看| 亚洲中文日韩欧美视频| 亚洲视频免费观看视频| e午夜精品久久久久久久| 超碰97精品在线观看| 国产福利在线免费观看视频| 我的亚洲天堂| 亚洲国产看品久久| 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| 久久精品人人爽人人爽视色| av一本久久久久| 新久久久久国产一级毛片| 国产精品偷伦视频观看了| 亚洲自偷自拍图片 自拍| 久久精品国产亚洲av涩爱| 一级毛片电影观看| 男女午夜视频在线观看| 亚洲五月婷婷丁香| 最近手机中文字幕大全| 久久av网站| 一本色道久久久久久精品综合| 建设人人有责人人尽责人人享有的| 亚洲九九香蕉| 9色porny在线观看| 日韩视频在线欧美| 少妇人妻 视频| 亚洲成人免费电影在线观看 | 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 亚洲精品国产一区二区精华液| www.自偷自拍.com| 只有这里有精品99| 国产老妇伦熟女老妇高清| 人人妻人人澡人人爽人人夜夜| 国产片内射在线| 国产野战对白在线观看| 免费观看a级毛片全部| 人人妻人人添人人爽欧美一区卜| 免费看av在线观看网站| 国产成人精品久久久久久| 国产免费视频播放在线视频| 午夜福利乱码中文字幕| 性少妇av在线| 国产成人欧美| 男女边摸边吃奶| 丝袜美足系列| 国产精品二区激情视频| 午夜久久久在线观看| 久久久国产欧美日韩av| 在线观看免费午夜福利视频| 日韩 亚洲 欧美在线| 人体艺术视频欧美日本| 精品第一国产精品| 亚洲精品久久成人aⅴ小说| 国产成人免费观看mmmm| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 亚洲欧美一区二区三区久久| 尾随美女入室| 国产精品免费大片| 日韩制服丝袜自拍偷拍| 午夜av观看不卡| 亚洲三区欧美一区| 纵有疾风起免费观看全集完整版| 99久久99久久久精品蜜桃| 精品少妇黑人巨大在线播放| 香蕉国产在线看| 一本综合久久免费| 成人免费观看视频高清| 久久久久精品国产欧美久久久 | 好男人视频免费观看在线| 男女床上黄色一级片免费看| 欧美中文综合在线视频| 国产91精品成人一区二区三区 | 日韩视频在线欧美| 日韩制服骚丝袜av| 天天躁夜夜躁狠狠躁躁| 欧美人与性动交α欧美精品济南到| 777久久人妻少妇嫩草av网站| 国产熟女欧美一区二区| 国精品久久久久久国模美| 一二三四社区在线视频社区8| 激情五月婷婷亚洲| 男女无遮挡免费网站观看| 久久久国产一区二区| 搡老岳熟女国产| 在线av久久热| 国产野战对白在线观看| 免费不卡黄色视频| 丝瓜视频免费看黄片| 欧美日韩一级在线毛片| 自线自在国产av| 精品少妇黑人巨大在线播放| 欧美精品av麻豆av| av欧美777| 欧美日韩亚洲国产一区二区在线观看 | 色精品久久人妻99蜜桃| 久久久久网色| 免费观看av网站的网址| 欧美性长视频在线观看| 国产成人精品久久久久久| 欧美精品啪啪一区二区三区 | 免费高清在线观看日韩| 国产av国产精品国产| 大香蕉久久网| xxx大片免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 国产在线一区二区三区精| 18禁黄网站禁片午夜丰满| 久久久久久久大尺度免费视频| 国产熟女午夜一区二区三区| 久久精品久久久久久久性| 只有这里有精品99| 国产精品国产三级国产专区5o| 这个男人来自地球电影免费观看| a 毛片基地| 又大又爽又粗| 久久久久久免费高清国产稀缺| 亚洲一区二区三区欧美精品| 日韩 亚洲 欧美在线| 十八禁人妻一区二区| 亚洲av综合色区一区| 亚洲人成电影免费在线| 欧美+亚洲+日韩+国产| 高清视频免费观看一区二区| 日韩熟女老妇一区二区性免费视频| 亚洲久久久国产精品| 我的亚洲天堂| 在线天堂中文资源库| 丁香六月欧美| 国产成人精品在线电影| 国产成人免费无遮挡视频| 午夜福利在线免费观看网站| 天天躁日日躁夜夜躁夜夜| 51午夜福利影视在线观看| 亚洲国产精品成人久久小说| 老鸭窝网址在线观看| 97精品久久久久久久久久精品| 一级,二级,三级黄色视频| 一级黄片播放器| 国产91精品成人一区二区三区 | 国产免费现黄频在线看| 大片免费播放器 马上看| 99国产综合亚洲精品| 一本大道久久a久久精品| 满18在线观看网站| 欧美精品亚洲一区二区| 日韩一本色道免费dvd| 另类精品久久| 操出白浆在线播放| 狠狠精品人妻久久久久久综合| 国产麻豆69| 国产日韩欧美在线精品| 熟女少妇亚洲综合色aaa.| 另类亚洲欧美激情| 久久热在线av| 黄色一级大片看看| 亚洲成人手机| 叶爱在线成人免费视频播放| 我要看黄色一级片免费的| 精品福利永久在线观看| 一级黄片播放器| 观看av在线不卡| 新久久久久国产一级毛片| 国产精品人妻久久久影院| av欧美777| 精品第一国产精品| www.av在线官网国产| 亚洲国产欧美网| 欧美在线一区亚洲| 亚洲欧美激情在线| 麻豆国产av国片精品| av天堂在线播放| 91麻豆av在线| 欧美变态另类bdsm刘玥| 亚洲国产精品一区三区| 免费看不卡的av| 国产99久久九九免费精品| 日韩熟女老妇一区二区性免费视频| 高清欧美精品videossex| 欧美xxⅹ黑人| 又紧又爽又黄一区二区| 观看av在线不卡| 男女午夜视频在线观看| 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| 男女国产视频网站| 久久鲁丝午夜福利片| 99久久99久久久精品蜜桃| 51午夜福利影视在线观看| 国产黄色视频一区二区在线观看| 亚洲精品日韩在线中文字幕| 亚洲av综合色区一区| 性少妇av在线| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 黑丝袜美女国产一区| a级毛片黄视频| 只有这里有精品99| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 美女午夜性视频免费| 久久ye,这里只有精品| 日韩,欧美,国产一区二区三区| 男男h啪啪无遮挡| 亚洲欧美一区二区三区黑人| 国产视频一区二区在线看| 18禁黄网站禁片午夜丰满| 在线观看一区二区三区激情| 国产片内射在线| 女人久久www免费人成看片| 欧美精品啪啪一区二区三区 | 日本黄色日本黄色录像| 一本久久精品| 最新的欧美精品一区二区| 少妇粗大呻吟视频| 99久久精品国产亚洲精品|