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

    帶導(dǎo)轉(zhuǎn)槽火箭筒發(fā)射過程內(nèi)流場及受力分析

    2022-12-30 12:09:26毛聰聰阮文俊孫雪明賈昊楠
    關(guān)鍵詞:發(fā)射筒火箭彈壁面

    毛聰聰,阮文俊,孫雪明,賈昊楠

    (1.南京理工大學(xué) 能源與動(dòng)力工程學(xué)院,江蘇 南京 210094;2.江西星火軍工工業(yè)有限公司,江西 南昌 331709;3.中國人民解放軍32181部隊(duì),河北 石家莊 050000)

    火箭彈發(fā)射過程中,發(fā)射筒會(huì)經(jīng)受高溫、高壓、高速燃?xì)馍淞鞯臎_擊。在對(duì)火箭發(fā)射筒進(jìn)行力學(xué)分析時(shí),也需將發(fā)射筒內(nèi)表面壓力數(shù)據(jù)作為邊界條件,為不同狀態(tài)下的發(fā)射筒有限元分析提供載荷數(shù)據(jù)[1-4]。尤其對(duì)帶導(dǎo)轉(zhuǎn)槽的火箭武器系統(tǒng)進(jìn)行發(fā)射過程分析時(shí),需要考慮由于導(dǎo)轉(zhuǎn)槽存在而導(dǎo)致的不對(duì)稱力帶來的影響,因此有必要對(duì)帶導(dǎo)轉(zhuǎn)槽的火箭發(fā)射筒三維內(nèi)流場及發(fā)射裝置內(nèi)的不對(duì)稱力進(jìn)行研究。

    國內(nèi)外許多科研人員對(duì)燃?xì)馍淞鞯臄?shù)值模擬做了大量研究。張磊等[5]通過大渦模擬研究了不同噴管口徑對(duì)燃?xì)馍淞髁鲌龅挠绊憽`嚧蝴惖萚6]建立了導(dǎo)彈發(fā)射過程二維流場計(jì)算模型,研究了不同時(shí)刻發(fā)射筒的內(nèi)流場和外流場特征。張俊等[7]通過研究某固體火箭發(fā)動(dòng)機(jī)燃?xì)饬鲌龅姆植迹玫搅肆鲌鲚S向與徑向的燃?xì)馍淞鲄?shù)分布規(guī)律。王華等[8]建立了多噴管流場計(jì)算模型,研究了燃?xì)馍淞髁鲌龅乃俣扰c溫度分布。Sinha等[9]對(duì)燃?xì)馍淞鳑_擊發(fā)射筒過程進(jìn)行數(shù)值模擬,將不同網(wǎng)格數(shù)量和湍流模型得到的計(jì)算結(jié)果進(jìn)行了對(duì)比分析。

    為了使數(shù)值模擬更加接近實(shí)際的發(fā)射過程,筆者采用動(dòng)網(wǎng)格技術(shù)對(duì)火箭彈發(fā)射時(shí)的運(yùn)動(dòng)過程進(jìn)行模擬[10-12]。通過建立火箭彈發(fā)射過程發(fā)射筒內(nèi)流場計(jì)算模型,使用SSTk-ω湍流模型對(duì)非穩(wěn)態(tài)流場進(jìn)行數(shù)值模擬,并對(duì)帶導(dǎo)轉(zhuǎn)槽的火箭發(fā)射筒內(nèi)流場和非對(duì)稱力進(jìn)行研究。

    1 理論基礎(chǔ)

    1.1 控制方程

    質(zhì)量守恒方程:

    (1)

    動(dòng)量方程:

    (2)

    能量方程:

    (3)

    1.2 湍流模型

    筆者采用的湍流模型為SSTk-ω湍流模型,與標(biāo)準(zhǔn)k-ω模型相比,該模型在湍流粘度定義中考慮了湍流剪切應(yīng)力的輸運(yùn)過程,是對(duì)標(biāo)準(zhǔn)k-ω模型的增強(qiáng)。表達(dá)式為

    (4)

    (5)

    式中:GK表示湍流的動(dòng)能;Gω為ω的方程;Γk和Γω分別代表k與ω的發(fā)散項(xiàng);Dω代表正交發(fā)散項(xiàng)。

    2 計(jì)算模型及網(wǎng)格劃分

    2.1 發(fā)射筒實(shí)體模型及計(jì)算模型

    某單管火箭武器發(fā)射筒的實(shí)體模型如圖1所示。將發(fā)射筒尾部的凸臺(tái)結(jié)構(gòu)簡化,并將導(dǎo)轉(zhuǎn)槽簡化為連續(xù)的螺旋槽后,根據(jù)實(shí)體模型建立的三維仿真計(jì)算物理模型如圖2所示。

    因?yàn)椴恍枰紤]火箭彈受到燃?xì)馍淞鞯臎_擊狀況,所以將火箭彈簡化為噴管。在發(fā)射筒力學(xué)分析時(shí)由于需要考慮導(dǎo)轉(zhuǎn)槽處的應(yīng)力集中,因此可以通過三維帶導(dǎo)轉(zhuǎn)槽的發(fā)射筒內(nèi)流場計(jì)算,得到螺旋導(dǎo)轉(zhuǎn)槽處不同時(shí)刻表面的壓力分布。

    2.2 初始及邊界條件

    流場初始化壓力大小為101.325 kPa,溫度為300 K。因涉密原因,文中噴管的入口壓力和火箭彈運(yùn)動(dòng)速度都經(jīng)過無量綱化處理。將噴管進(jìn)口設(shè)置為壓力入口,入口燃?xì)鉁囟葹? 500 K?;鸺龔棸l(fā)射時(shí)的入口壓力隨時(shí)間變化曲線由實(shí)驗(yàn)測定,選取實(shí)驗(yàn)曲線中壓力到達(dá)0.036 6的時(shí)刻作為計(jì)算的起始時(shí)刻,火箭彈運(yùn)動(dòng)出筒時(shí)刻作為計(jì)算的終止時(shí)刻,共歷時(shí)130 ms。從開始計(jì)算到火箭彈運(yùn)動(dòng)出筒的噴管入口壓力時(shí)間曲線如圖3所示。將發(fā)射筒壁面設(shè)定為絕熱無滑移壁面,發(fā)射筒尾端設(shè)置為壓力出口,壓力大小為101.325 kPa,溫度為300 K。根據(jù)噴管入口壓力隨時(shí)間變化曲線得到火箭彈運(yùn)動(dòng)速度隨時(shí)間變化的曲線,如圖4所示。

    2.3 計(jì)算條件設(shè)定

    流場計(jì)算選擇基于壓力的瞬態(tài)求解器,將燃?xì)猱?dāng)作理想氣體可壓縮流處理。不考慮火箭彈發(fā)射時(shí)的空氣阻力;不考慮燃?xì)饣瘜W(xué)反應(yīng)和固體顆粒的影響;不考慮燃?xì)馀c發(fā)射筒內(nèi)壁及周圍空氣的傳熱;不考慮火箭彈的旋轉(zhuǎn),僅考慮火箭彈沿水平方向的運(yùn)動(dòng)。流體粘性系數(shù)采用Sutherland公式計(jì)算。

    2.4 網(wǎng)格劃分

    在劃分網(wǎng)格時(shí),由于較難生成單塊的高質(zhì)量結(jié)構(gòu)化網(wǎng)格,因此采用分塊網(wǎng)格對(duì)計(jì)算域進(jìn)行劃分,4個(gè)網(wǎng)格劃分區(qū)域分別為噴管、尾部計(jì)算域、導(dǎo)轉(zhuǎn)槽和筒體。筆者采用Icem軟件劃分六面體結(jié)構(gòu)網(wǎng)格,對(duì)邊界區(qū)域添加了邊界層網(wǎng)格后,初始時(shí)刻的網(wǎng)格總數(shù)為60萬個(gè)。初始時(shí)刻網(wǎng)格劃分如圖5所示。

    2.5 動(dòng)網(wǎng)格設(shè)置

    為了處理火箭彈發(fā)射時(shí)的網(wǎng)格變化,在流場計(jì)算時(shí)通過動(dòng)網(wǎng)格對(duì)火箭彈運(yùn)動(dòng)進(jìn)行模擬。噴管網(wǎng)格水平移動(dòng),尾部計(jì)算域網(wǎng)格被拉長,筒體與導(dǎo)轉(zhuǎn)槽網(wǎng)格固定不動(dòng)。動(dòng)網(wǎng)格更新方法包括彈性光順法、動(dòng)態(tài)層法和局部重構(gòu)法,根據(jù)本文的網(wǎng)格劃分方式,采用動(dòng)態(tài)層法模擬火箭彈運(yùn)動(dòng)。該方法的基本思想是:根據(jù)與運(yùn)動(dòng)邊界相鄰網(wǎng)格的高度,判斷網(wǎng)格需要分割還是合并,當(dāng)網(wǎng)格高度滿足式(6)時(shí),網(wǎng)格將會(huì)被分割;而當(dāng)網(wǎng)格高度滿足式(7)時(shí),兩層網(wǎng)格將會(huì)合并。

    h≥(1+α)h0,

    (6)

    h≤βh0,

    (7)

    式中:α為網(wǎng)格層的分裂因子;β為網(wǎng)格層的合并因子;h0為理想單元網(wǎng)格高度。

    筆者取h0=3 mm,α=0.4,β=0.2。通過Profile文件將火箭彈的運(yùn)動(dòng)速度賦給噴管及尾部計(jì)算域,噴管不斷向前運(yùn)動(dòng),尾部計(jì)算域網(wǎng)格通過動(dòng)態(tài)層法進(jìn)行更新生成。初始時(shí)刻的運(yùn)動(dòng)域網(wǎng)格如圖6所示,火箭彈發(fā)射時(shí)運(yùn)動(dòng)域的網(wǎng)格變化如圖7所示。

    3 數(shù)值模擬結(jié)果及分析

    3.1 不同時(shí)刻流場縱截面壓力云圖

    計(jì)算時(shí)刻分別為50、90和130 ms時(shí)的流場縱截面壓力云圖,如圖8所示。

    因涉密原因,流場計(jì)算得到的壓力、速度及扭轉(zhuǎn)力都經(jīng)過了無量綱化處理。由于噴管內(nèi)壓力變化較大,因此為了更清楚顯示發(fā)射筒內(nèi)的壓力分布,將無量綱壓力在0.366以上的部分用相同顏色顯示。燃?xì)庠趪姽軆?nèi)由亞音速加速到超音速后,在噴管出口形成了超音速燃?xì)馍淞?。按膨脹程度劃分,燃?xì)馍淞鳛榍放蛎浬淞?。燃?xì)馍淞鲝膰姽艹隹诹鞒鰰r(shí)的壓力高于周圍環(huán)境壓力,所以一出噴管首先會(huì)膨脹加速。燃?xì)饨?jīng)過噴管出口兩側(cè)的膨脹扇區(qū)后,繼續(xù)膨脹加速同時(shí)壓力不斷降低。當(dāng)噴管出口的膨脹波系到達(dá)射流邊界后,向內(nèi)反射為一系列的壓縮波,壓縮波足夠強(qiáng)時(shí)會(huì)匯合形成斜激波。燃?xì)饨?jīng)過激波后,壓力增加、溫度升高同時(shí)速度降低。當(dāng)燃?xì)鈮毫υ黾拥礁哂诃h(huán)境壓力后,燃?xì)馀蛎浖铀偻瑫r(shí)壓力降低。隨著這樣膨脹壓縮過程的反復(fù)循環(huán),發(fā)射筒內(nèi)形成了一系列膨脹波與壓縮波結(jié)構(gòu)。

    3.2 中心軸線上的速度與壓力分布

    圖9和圖10分別為不同時(shí)刻流場中心軸線上的速度與壓力分布。從圖9和圖10中可以看出,燃?xì)馍淞鲝膰姽芰鞒龊?,噴管出口附近的膨脹波和壓縮波強(qiáng)度是最大的,因此燃?xì)庠趪姽芨浇乃俣扰c壓力波動(dòng)也最為明顯。在50 ms時(shí),燃?xì)庠谥行妮S線上經(jīng)歷了1~2個(gè)膨脹-壓縮波;在90 ms時(shí),燃?xì)庠谥行妮S線上經(jīng)歷了5~6個(gè)膨脹-壓縮波;在130 ms時(shí),燃?xì)庠谥行妮S線上經(jīng)歷了12~13個(gè)膨脹-壓縮波。因發(fā)射筒長度有限,所以當(dāng)燃?xì)獾竭_(dá)發(fā)射筒出口附近時(shí),燃?xì)獾乃俣扰c壓力仍有較大波動(dòng)。由于粘性效應(yīng)和周圍空氣的作用產(chǎn)生了湍流剪切層,湍流剪切層內(nèi)的粘性阻尼使得膨脹波與壓縮波結(jié)構(gòu)不斷瓦解,因此中心軸線上的燃?xì)馑俣扰c壓力呈現(xiàn)振蕩衰減的趨勢。

    3.3 不同時(shí)刻發(fā)射筒內(nèi)表面壓力分布

    分別取計(jì)算時(shí)刻為50、90和130 ms時(shí)的發(fā)射筒內(nèi)表面壓力云圖,如圖11所示。燃?xì)鈴膰姽芰鞒龊?,位于噴管出口后的發(fā)射筒區(qū)域直接受到燃?xì)馍淞鞯膹?qiáng)烈沖擊,因此其表面壓力明顯高于噴管前部未直接受到燃?xì)馍淞鳑_擊部分的壓力。

    當(dāng)燃?xì)饨?jīng)歷了在壁面附近的激波后壓力增加,其對(duì)發(fā)射筒壁面的壓力也隨之增加。50 ms時(shí),燃?xì)鈱?duì)發(fā)射筒壁面的最大壓力為0.078;90 ms時(shí),燃?xì)鈱?duì)發(fā)射筒壁面的最大壓力為0.065;130 ms時(shí),燃?xì)鈱?duì)發(fā)射筒壁面的最大壓力為0.049。當(dāng)燃?xì)鈮毫υ黾拥礁哂谥車諝獾膲毫?,燃?xì)馀蛎浖铀俨⑶覊毫Σ粩嘟档?,其?duì)壁面的壓力也逐漸降低。燃?xì)庠诎l(fā)射筒內(nèi)經(jīng)歷了一系列膨脹壓縮波,所以燃?xì)鈱?duì)發(fā)射筒內(nèi)表面的壓力在整體上呈現(xiàn)振蕩衰減的趨勢。

    3.4 不同時(shí)刻導(dǎo)轉(zhuǎn)槽內(nèi)的壓力及扭轉(zhuǎn)力分布

    不同時(shí)刻導(dǎo)轉(zhuǎn)槽表面壓力隨距離變化曲線如圖12所示。

    在發(fā)射筒右視圖中,位于導(dǎo)轉(zhuǎn)槽左側(cè)的壁面為本文的導(dǎo)轉(zhuǎn)槽左側(cè)面,位于導(dǎo)轉(zhuǎn)槽右側(cè)的壁面為本文的導(dǎo)轉(zhuǎn)槽右側(cè)面。50 ms時(shí),燃?xì)鈱?duì)導(dǎo)轉(zhuǎn)槽頂面、右側(cè)面和左側(cè)面的最大壓力分別為0.072,0.070和0.065;90 ms時(shí),燃?xì)鈱?duì)導(dǎo)轉(zhuǎn)槽頂面、右側(cè)面和左側(cè)面的最大壓力分別為0.065,0.060和0.055;130 ms時(shí),燃?xì)鈱?duì)導(dǎo)轉(zhuǎn)槽頂面、右側(cè)面和左側(cè)面的最大壓力分別為0.045,0.044和0.042。因50 ms時(shí)噴管入口壓力是3個(gè)時(shí)刻中最大的,故燃?xì)鈱?duì)導(dǎo)轉(zhuǎn)槽3個(gè)表面的最大壓力也是3個(gè)時(shí)刻中最大的。在不同時(shí)刻,導(dǎo)轉(zhuǎn)槽頂面受到的最大壓力值始終是3個(gè)表面中最大的,而導(dǎo)轉(zhuǎn)槽右側(cè)面受到的最大壓力略大于導(dǎo)轉(zhuǎn)槽左側(cè)面。

    從不同時(shí)刻導(dǎo)轉(zhuǎn)槽表面的壓力分布曲線中可以看出,導(dǎo)轉(zhuǎn)槽表面的壓力沿著中心軸線方向振蕩變化,但并未呈現(xiàn)振蕩衰減的趨勢。這是由于導(dǎo)轉(zhuǎn)槽表面會(huì)多次反射壓縮波,當(dāng)這些反射的壓縮波疊加后,燃?xì)庠俅谓?jīng)過這些壓縮波會(huì)增加對(duì)導(dǎo)轉(zhuǎn)槽表面的壓力,因此導(dǎo)轉(zhuǎn)槽表面的壓力并未沿著軸線方向振蕩衰減。

    導(dǎo)轉(zhuǎn)槽內(nèi)扭轉(zhuǎn)力隨時(shí)間變化的曲線如圖13所示。導(dǎo)轉(zhuǎn)槽兩個(gè)側(cè)面受到燃?xì)馍淞髦苯幼饔妹娣e隨著火箭彈向前運(yùn)動(dòng)而逐漸增加,故導(dǎo)轉(zhuǎn)槽兩個(gè)側(cè)面受到的燃?xì)庾饔昧σ苍诓粩嘣黾?。?dǎo)轉(zhuǎn)槽兩個(gè)側(cè)面上受到的燃?xì)馍淞髯饔昧Σ钪凳沟冒l(fā)射裝置內(nèi)產(chǎn)生了不對(duì)稱力即扭轉(zhuǎn)力,而右側(cè)面受到的燃?xì)馍淞髯饔昧κ冀K大于左側(cè)面,因此導(dǎo)轉(zhuǎn)槽內(nèi)扭轉(zhuǎn)力的方向始終垂直于右側(cè)面。扭轉(zhuǎn)力隨著時(shí)間呈現(xiàn)振蕩上升的趨勢,在火箭彈即將出筒時(shí)刻達(dá)到最大,最大的扭轉(zhuǎn)力值為769。

    4 結(jié)論

    筆者通過動(dòng)網(wǎng)格技術(shù)模擬了火箭彈從發(fā)射到出筒過程燃?xì)馍淞鲗?duì)發(fā)射筒的沖擊。通過流場計(jì)算得到不同時(shí)刻發(fā)射筒內(nèi)表面的壓力分布,將其加載到包含導(dǎo)轉(zhuǎn)槽結(jié)構(gòu)的三維發(fā)射筒有限元分析模型中進(jìn)行計(jì)算。可以得到發(fā)射筒在發(fā)射壓力作用下的應(yīng)力應(yīng)變情況,從而確定發(fā)射筒在發(fā)射壓力作用下的相對(duì)薄弱部位以及安全儲(chǔ)備。從本文的數(shù)值模擬中可以得到以下結(jié)論:

    1)超音速燃?xì)馍淞髋c環(huán)境大氣之間的壓力差,導(dǎo)致發(fā)射筒內(nèi)產(chǎn)生了一系列膨脹壓縮波。由于流體粘性效應(yīng)和周圍大氣的作用,膨脹波與壓縮波結(jié)構(gòu)在發(fā)射筒內(nèi)逐漸崩解。

    2)燃?xì)庠诎l(fā)射筒內(nèi)經(jīng)歷了一系列膨脹波和壓縮波,因此不同時(shí)刻中心軸線上的燃?xì)鈮毫?、燃?xì)馑俣群桶l(fā)射筒內(nèi)壁面上的壓力,在整體上呈現(xiàn)振蕩衰減的趨勢。

    3)導(dǎo)轉(zhuǎn)槽內(nèi)壓縮波的反射疊加,使得不同時(shí)刻導(dǎo)轉(zhuǎn)槽表面上的壓力并未沿著中心軸線振蕩衰減。導(dǎo)轉(zhuǎn)槽的存在導(dǎo)致火箭發(fā)射筒內(nèi)產(chǎn)生了扭轉(zhuǎn)力,并且隨著火箭彈向前運(yùn)動(dòng),扭轉(zhuǎn)力振蕩增加。

    猜你喜歡
    發(fā)射筒火箭彈壁面
    “鐵穹”反火箭彈系統(tǒng)
    二維有限長度柔性壁面上T-S波演化的數(shù)值研究
    某型增強(qiáng)型復(fù)合材料發(fā)射筒強(qiáng)度特性評(píng)估
    發(fā)射筒熱力耦合場下力學(xué)性能研究
    地空導(dǎo)彈發(fā)射筒焊接系統(tǒng)的校準(zhǔn)技術(shù)研究
    電子測試(2018年9期)2018-06-26 06:45:46
    降低某產(chǎn)品貯運(yùn)發(fā)射筒內(nèi)氣體濕度的工藝方法研究
    壁面溫度對(duì)微型內(nèi)燃機(jī)燃燒特性的影響
    顆粒—壁面碰撞建模與數(shù)據(jù)處理
    考慮裂縫壁面?zhèn)Φ膲毫丫a(chǎn)能計(jì)算模型
    色av中文字幕| 国产成年人精品一区二区| 国产精品一区二区精品视频观看| 一级毛片高清免费大全| 久久久久久久精品吃奶| 日韩欧美国产一区二区入口| 国产亚洲精品久久久久久毛片| 国产精品一区二区免费欧美| 性欧美人与动物交配| 国产精品香港三级国产av潘金莲| 亚洲午夜精品一区,二区,三区| 国产亚洲欧美精品永久| 国产午夜精品久久久久久| 午夜免费激情av| 国产精品免费一区二区三区在线| 精品国产国语对白av| 国产主播在线观看一区二区| 性欧美人与动物交配| 性色av乱码一区二区三区2| 久久中文字幕人妻熟女| 国产主播在线观看一区二区| 亚洲欧美激情在线| 欧美成人一区二区免费高清观看 | 女人高潮潮喷娇喘18禁视频| 12—13女人毛片做爰片一| 日韩中文字幕欧美一区二区| 久9热在线精品视频| 久久久水蜜桃国产精品网| 日韩欧美国产一区二区入口| 在线观看一区二区三区| 久99久视频精品免费| 高清毛片免费观看视频网站| 亚洲黑人精品在线| av网站免费在线观看视频| 少妇裸体淫交视频免费看高清 | 国产精品美女特级片免费视频播放器 | 免费无遮挡裸体视频| 高清在线国产一区| 一级作爱视频免费观看| 在线十欧美十亚洲十日本专区| or卡值多少钱| 欧美日本视频| 大码成人一级视频| 欧美另类亚洲清纯唯美| 法律面前人人平等表现在哪些方面| bbb黄色大片| cao死你这个sao货| 1024视频免费在线观看| 91在线观看av| 国产国语露脸激情在线看| 亚洲av电影不卡..在线观看| 人妻丰满熟妇av一区二区三区| 亚洲av电影不卡..在线观看| 热99re8久久精品国产| 91精品三级在线观看| 久久久久久人人人人人| 国产91精品成人一区二区三区| 色综合欧美亚洲国产小说| 亚洲成人精品中文字幕电影| 真人做人爱边吃奶动态| 无人区码免费观看不卡| 香蕉丝袜av| 宅男免费午夜| 国产成年人精品一区二区| 午夜a级毛片| 一级a爱片免费观看的视频| 亚洲精品一区av在线观看| 精品不卡国产一区二区三区| 妹子高潮喷水视频| 欧美国产精品va在线观看不卡| 久久精品亚洲精品国产色婷小说| 久久性视频一级片| 亚洲国产精品成人综合色| 亚洲自拍偷在线| 免费在线观看影片大全网站| 久久影院123| 老司机午夜福利在线观看视频| 99国产综合亚洲精品| 国产一区二区三区综合在线观看| 国产精品精品国产色婷婷| 两个人视频免费观看高清| 国产欧美日韩精品亚洲av| 日韩欧美一区二区三区在线观看| 涩涩av久久男人的天堂| 看黄色毛片网站| 亚洲视频免费观看视频| 日本 av在线| 国产91精品成人一区二区三区| 国产精品,欧美在线| 久久精品人人爽人人爽视色| 国产成人精品在线电影| 丝袜美腿诱惑在线| 啪啪无遮挡十八禁网站| 久久亚洲真实| 欧美日韩亚洲综合一区二区三区_| 露出奶头的视频| 日本a在线网址| 一级,二级,三级黄色视频| 欧美激情极品国产一区二区三区| 亚洲成人国产一区在线观看| 亚洲精品久久成人aⅴ小说| 搡老熟女国产l中国老女人| 宅男免费午夜| 美女高潮喷水抽搐中文字幕| cao死你这个sao货| 国产区一区二久久| 国产精品 欧美亚洲| 欧美大码av| 最近最新中文字幕大全电影3 | 成人国语在线视频| 亚洲欧美一区二区三区黑人| 成在线人永久免费视频| 性色av乱码一区二区三区2| 免费人成视频x8x8入口观看| 视频在线观看一区二区三区| 亚洲五月天丁香| 操出白浆在线播放| 高清在线国产一区| 熟妇人妻久久中文字幕3abv| 99久久久亚洲精品蜜臀av| 好男人电影高清在线观看| 午夜影院日韩av| videosex国产| 夜夜爽天天搞| 亚洲五月婷婷丁香| 18禁美女被吸乳视频| 日韩精品青青久久久久久| 日韩欧美一区视频在线观看| 成在线人永久免费视频| 18禁裸乳无遮挡免费网站照片 | 波多野结衣高清无吗| 日韩成人在线观看一区二区三区| 日本精品一区二区三区蜜桃| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 成人国语在线视频| 国产人伦9x9x在线观看| 精品卡一卡二卡四卡免费| 激情视频va一区二区三区| 精品少妇一区二区三区视频日本电影| 国产一级毛片七仙女欲春2 | 亚洲三区欧美一区| 九色亚洲精品在线播放| 国产亚洲精品久久久久久毛片| 午夜两性在线视频| 在线观看66精品国产| 国产一级毛片七仙女欲春2 | 久久久久久久久免费视频了| 性欧美人与动物交配| 性欧美人与动物交配| 精品久久蜜臀av无| 亚洲男人的天堂狠狠| 高清在线国产一区| 国产片内射在线| 女人精品久久久久毛片| 久久中文字幕人妻熟女| 国产一区二区激情短视频| avwww免费| 丰满的人妻完整版| 国产精品一区二区三区四区久久 | 成熟少妇高潮喷水视频| 亚洲精品国产一区二区精华液| 免费不卡黄色视频| 狂野欧美激情性xxxx| 久久影院123| 制服人妻中文乱码| 精品一区二区三区av网在线观看| 国产av一区在线观看免费| 欧美一级a爱片免费观看看 | 精品人妻在线不人妻| 波多野结衣av一区二区av| 村上凉子中文字幕在线| 老熟妇仑乱视频hdxx| 日韩 欧美 亚洲 中文字幕| 一a级毛片在线观看| 国产真人三级小视频在线观看| 视频区欧美日本亚洲| 两性夫妻黄色片| 亚洲天堂国产精品一区在线| 又黄又爽又免费观看的视频| 男女午夜视频在线观看| 亚洲欧美日韩无卡精品| 久久国产精品人妻蜜桃| svipshipincom国产片| 在线观看舔阴道视频| 久久久久久免费高清国产稀缺| 日本欧美视频一区| 亚洲人成77777在线视频| www.www免费av| 国产精品永久免费网站| 国产1区2区3区精品| www.999成人在线观看| 夜夜爽天天搞| 日韩欧美国产在线观看| 少妇 在线观看| 亚洲精品久久成人aⅴ小说| 国产精华一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 不卡av一区二区三区| 亚洲欧美日韩无卡精品| 国产一区二区在线av高清观看| 69av精品久久久久久| 成人三级做爰电影| 一级a爱视频在线免费观看| 国产私拍福利视频在线观看| 午夜精品在线福利| 免费在线观看日本一区| 大型黄色视频在线免费观看| 一个人免费在线观看的高清视频| 十分钟在线观看高清视频www| 午夜福利影视在线免费观看| 国产精品久久久av美女十八| 亚洲一码二码三码区别大吗| 久久伊人香网站| 91九色精品人成在线观看| 成人av一区二区三区在线看| √禁漫天堂资源中文www| 操出白浆在线播放| 精品久久蜜臀av无| 亚洲专区字幕在线| 美女免费视频网站| 国产高清有码在线观看视频 | 最近最新免费中文字幕在线| 亚洲欧美激情在线| 亚洲欧美激情在线| 99热只有精品国产| 午夜福利视频1000在线观看 | 久久国产精品影院| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 好男人电影高清在线观看| 成人国语在线视频| 日韩大尺度精品在线看网址 | 女人爽到高潮嗷嗷叫在线视频| 久久国产精品人妻蜜桃| 麻豆久久精品国产亚洲av| 免费无遮挡裸体视频| 国产午夜福利久久久久久| 精品国产一区二区三区四区第35| 宅男免费午夜| 黄频高清免费视频| 亚洲精品国产一区二区精华液| 亚洲精品一卡2卡三卡4卡5卡| 国产主播在线观看一区二区| 亚洲五月婷婷丁香| 午夜成年电影在线免费观看| 一区二区三区国产精品乱码| 一区在线观看完整版| 99riav亚洲国产免费| e午夜精品久久久久久久| 国产免费男女视频| 两个人免费观看高清视频| 亚洲五月天丁香| 国内久久婷婷六月综合欲色啪| 成年版毛片免费区| 黄片小视频在线播放| 国产欧美日韩一区二区精品| netflix在线观看网站| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 午夜老司机福利片| 中文字幕人妻熟女乱码| 国产不卡一卡二| av欧美777| 亚洲aⅴ乱码一区二区在线播放 | 夜夜看夜夜爽夜夜摸| 国产成人系列免费观看| 国产成年人精品一区二区| 久久天躁狠狠躁夜夜2o2o| 国产精品久久电影中文字幕| 欧美久久黑人一区二区| 亚洲av日韩精品久久久久久密| 欧美日韩乱码在线| 免费在线观看亚洲国产| 99精品在免费线老司机午夜| 此物有八面人人有两片| 91成年电影在线观看| 性少妇av在线| 50天的宝宝边吃奶边哭怎么回事| 色播在线永久视频| 亚洲欧美精品综合一区二区三区| 日本五十路高清| 18禁观看日本| 亚洲国产看品久久| 男女床上黄色一级片免费看| 国产熟女午夜一区二区三区| 欧美最黄视频在线播放免费| 午夜福利影视在线免费观看| 在线视频色国产色| 黄片小视频在线播放| 女人高潮潮喷娇喘18禁视频| 在线视频色国产色| 久久精品国产综合久久久| 一边摸一边做爽爽视频免费| 亚洲电影在线观看av| 十八禁网站免费在线| 国产极品粉嫩免费观看在线| 超碰成人久久| 深夜精品福利| 欧美色欧美亚洲另类二区 | 亚洲美女黄片视频| 久久久精品欧美日韩精品| www.精华液| 精品电影一区二区在线| 欧美午夜高清在线| 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 精品国产国语对白av| 国产精品亚洲美女久久久| 国产一区二区三区视频了| 婷婷六月久久综合丁香| 99国产综合亚洲精品| 亚洲成国产人片在线观看| 国产伦一二天堂av在线观看| 好看av亚洲va欧美ⅴa在| 久久九九热精品免费| 高潮久久久久久久久久久不卡| 亚洲专区字幕在线| 欧美老熟妇乱子伦牲交| 久久精品91无色码中文字幕| 国产激情欧美一区二区| 国产精品香港三级国产av潘金莲| 丝袜人妻中文字幕| 91九色精品人成在线观看| 人人妻人人爽人人添夜夜欢视频| 女警被强在线播放| 天堂影院成人在线观看| 制服人妻中文乱码| 欧美一区二区精品小视频在线| 老汉色av国产亚洲站长工具| 亚洲免费av在线视频| 亚洲av五月六月丁香网| 在线天堂中文资源库| 国产亚洲欧美98| 国产乱人伦免费视频| 18禁国产床啪视频网站| 麻豆国产av国片精品| 国产在线精品亚洲第一网站| 老汉色∧v一级毛片| 久久亚洲精品不卡| 一区二区三区精品91| 99精品在免费线老司机午夜| 丁香六月欧美| 神马国产精品三级电影在线观看 | 91老司机精品| 又紧又爽又黄一区二区| 男人舔女人下体高潮全视频| 国内毛片毛片毛片毛片毛片| 久久久国产成人精品二区| 亚洲国产日韩欧美精品在线观看 | 色哟哟哟哟哟哟| 日韩三级视频一区二区三区| 国产区一区二久久| 日日夜夜操网爽| 亚洲自拍偷在线| 俄罗斯特黄特色一大片| 亚洲自拍偷在线| 久久精品人人爽人人爽视色| aaaaa片日本免费| 曰老女人黄片| 精品国产美女av久久久久小说| 欧美色欧美亚洲另类二区 | or卡值多少钱| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 欧美成人一区二区免费高清观看 | 午夜激情av网站| 国产精品香港三级国产av潘金莲| 一边摸一边抽搐一进一小说| 亚洲情色 制服丝袜| 身体一侧抽搐| 免费av毛片视频| 一边摸一边抽搐一进一小说| 我的亚洲天堂| 最好的美女福利视频网| av超薄肉色丝袜交足视频| 中文字幕av电影在线播放| 美女午夜性视频免费| 欧美日韩一级在线毛片| 亚洲午夜精品一区,二区,三区| 亚洲成av片中文字幕在线观看| 免费不卡黄色视频| 亚洲熟女毛片儿| 一个人免费在线观看的高清视频| 日本三级黄在线观看| 18禁观看日本| 一个人观看的视频www高清免费观看 | 91麻豆精品激情在线观看国产| svipshipincom国产片| 国产主播在线观看一区二区| 伦理电影免费视频| 日韩欧美三级三区| 看免费av毛片| 嫩草影院精品99| 精品高清国产在线一区| 久久久久久大精品| 男男h啪啪无遮挡| 国产三级黄色录像| 精品高清国产在线一区| 亚洲九九香蕉| 99国产精品一区二区三区| 香蕉丝袜av| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| 精品不卡国产一区二区三区| 桃红色精品国产亚洲av| 亚洲精品国产区一区二| 免费在线观看黄色视频的| 老熟妇仑乱视频hdxx| 99精品久久久久人妻精品| 日韩免费av在线播放| 亚洲人成电影免费在线| 日本五十路高清| 很黄的视频免费| 国产精品久久久久久亚洲av鲁大| 久久久国产成人精品二区| 国产不卡一卡二| svipshipincom国产片| 亚洲激情在线av| 国产在线精品亚洲第一网站| 久久久久久大精品| 免费在线观看视频国产中文字幕亚洲| 满18在线观看网站| 精品少妇一区二区三区视频日本电影| 国产熟女午夜一区二区三区| 乱人伦中国视频| 极品人妻少妇av视频| 亚洲精品在线观看二区| 性欧美人与动物交配| 一区福利在线观看| 日韩欧美国产在线观看| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 久久久国产精品麻豆| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 熟妇人妻久久中文字幕3abv| 18禁美女被吸乳视频| 激情在线观看视频在线高清| 给我免费播放毛片高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 老鸭窝网址在线观看| 在线观看舔阴道视频| 亚洲精品美女久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产区一区二| 国产精品影院久久| 自线自在国产av| 俄罗斯特黄特色一大片| 免费观看精品视频网站| 国产精品香港三级国产av潘金莲| 日本撒尿小便嘘嘘汇集6| 在线观看www视频免费| 亚洲国产欧美一区二区综合| 国产精品乱码一区二三区的特点 | 桃色一区二区三区在线观看| 精品国产亚洲在线| 国产精品av久久久久免费| 日本三级黄在线观看| 欧美+亚洲+日韩+国产| 亚洲七黄色美女视频| 午夜成年电影在线免费观看| 999久久久国产精品视频| 99riav亚洲国产免费| 精品卡一卡二卡四卡免费| 在线免费观看的www视频| 亚洲男人的天堂狠狠| 国产男靠女视频免费网站| 国产亚洲av嫩草精品影院| 手机成人av网站| 成人亚洲精品一区在线观看| 久久人妻av系列| 97人妻精品一区二区三区麻豆 | 免费看美女性在线毛片视频| 久久青草综合色| 他把我摸到了高潮在线观看| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 人妻久久中文字幕网| 亚洲五月色婷婷综合| 日韩av在线大香蕉| 搞女人的毛片| 久久久水蜜桃国产精品网| 国产精品久久电影中文字幕| 亚洲欧美日韩无卡精品| 在线观看日韩欧美| 男女午夜视频在线观看| 亚洲欧美精品综合久久99| 午夜免费观看网址| 一区二区三区国产精品乱码| 国产成人精品无人区| 精品少妇一区二区三区视频日本电影| 亚洲人成电影免费在线| 美女大奶头视频| 久久久久国产一级毛片高清牌| 欧美黑人精品巨大| 看免费av毛片| 国产高清有码在线观看视频 | 精品日产1卡2卡| 国产精品,欧美在线| 国产av一区在线观看免费| 91麻豆精品激情在线观看国产| 国产精品一区二区免费欧美| 久久人人精品亚洲av| 露出奶头的视频| 波多野结衣巨乳人妻| 国产一区二区在线av高清观看| 亚洲国产欧美网| 99国产综合亚洲精品| 亚洲最大成人中文| 麻豆av在线久日| 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 人人妻人人澡欧美一区二区 | 欧美亚洲日本最大视频资源| av福利片在线| 嫩草影视91久久| 老司机午夜十八禁免费视频| 丁香欧美五月| √禁漫天堂资源中文www| av欧美777| 日韩欧美国产在线观看| 国产精品一区二区三区四区久久 | 亚洲精品国产色婷婷电影| 午夜免费观看网址| 久久国产精品人妻蜜桃| 国产av一区二区精品久久| 成人亚洲精品av一区二区| 国产男靠女视频免费网站| 国内精品久久久久久久电影| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 12—13女人毛片做爰片一| 欧美人与性动交α欧美精品济南到| 黑人巨大精品欧美一区二区mp4| а√天堂www在线а√下载| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 国产精品九九99| 久久青草综合色| 女生性感内裤真人,穿戴方法视频| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 欧美久久黑人一区二区| 国产成年人精品一区二区| 最新美女视频免费是黄的| 久久天堂一区二区三区四区| 久久午夜综合久久蜜桃| 久久久久久国产a免费观看| 欧美精品亚洲一区二区| 操出白浆在线播放| 国产蜜桃级精品一区二区三区| 国产成人免费无遮挡视频| 日本 欧美在线| 人人妻,人人澡人人爽秒播| 成人18禁高潮啪啪吃奶动态图| 国产精品自产拍在线观看55亚洲| 侵犯人妻中文字幕一二三四区| 日本撒尿小便嘘嘘汇集6| 日韩欧美三级三区| 亚洲最大成人中文| 免费在线观看黄色视频的| 男女做爰动态图高潮gif福利片 | 99国产精品免费福利视频| 亚洲国产精品合色在线| 三级毛片av免费| 国产欧美日韩一区二区三| 一区二区三区国产精品乱码| 99riav亚洲国产免费| 天天躁狠狠躁夜夜躁狠狠躁| 脱女人内裤的视频| av超薄肉色丝袜交足视频| 女人高潮潮喷娇喘18禁视频| 50天的宝宝边吃奶边哭怎么回事| а√天堂www在线а√下载| 悠悠久久av| 多毛熟女@视频| 午夜福利免费观看在线| 精品不卡国产一区二区三区| 久久精品国产99精品国产亚洲性色 | 嫩草影视91久久| 看免费av毛片| 亚洲片人在线观看| 日韩一卡2卡3卡4卡2021年| 国产乱人伦免费视频| а√天堂www在线а√下载| 色播在线永久视频| 欧美日韩亚洲综合一区二区三区_| 九色国产91popny在线| 亚洲,欧美精品.| 真人做人爱边吃奶动态| 亚洲第一青青草原| 精品无人区乱码1区二区| 国产99白浆流出| ponron亚洲| 国产精品永久免费网站| 亚洲精品美女久久av网站| 一二三四在线观看免费中文在| 亚洲人成伊人成综合网2020| 成人三级黄色视频| 一级毛片高清免费大全| 黑丝袜美女国产一区| 欧美亚洲日本最大视频资源| 亚洲国产欧美一区二区综合| 99久久久亚洲精品蜜臀av| www.精华液| 长腿黑丝高跟| 国产精品99久久99久久久不卡| 国产一区二区三区综合在线观看| 黄色成人免费大全| 亚洲aⅴ乱码一区二区在线播放 |