• <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片| 亚洲精品色激情综合| 日韩大片免费观看网站| 少妇人妻久久综合中文| 777米奇影视久久| 欧美精品国产亚洲| 在线观看美女被高潮喷水网站| 热99国产精品久久久久久7| 乱人伦中国视频| 国产在线视频一区二区| 国产午夜精品久久久久久一区二区三区| 国产精品免费大片| 亚洲成人手机| 纯流量卡能插随身wifi吗| 人妻一区二区av| 亚洲av中文av极速乱| 久久久午夜欧美精品| 欧美日韩精品成人综合77777| 一级二级三级毛片免费看| 久久精品夜色国产| av专区在线播放| 精品少妇黑人巨大在线播放| 免费高清在线观看视频在线观看| 免费av中文字幕在线| 最近手机中文字幕大全| 亚洲精品一区蜜桃| 免费av不卡在线播放| 少妇丰满av| 日本av手机在线免费观看| 亚洲精品第二区| 秋霞伦理黄片| 边亲边吃奶的免费视频| 亚洲国产欧美在线一区| 99久久精品国产国产毛片| av在线app专区| 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线 | 高清在线视频一区二区三区| 欧美最新免费一区二区三区| 一区二区av电影网| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 亚洲,欧美,日韩| 在线观看免费日韩欧美大片 | 亚洲中文av在线| 一级,二级,三级黄色视频| 又粗又硬又长又爽又黄的视频| 欧美 亚洲 国产 日韩一| 涩涩av久久男人的天堂| 高清午夜精品一区二区三区| 一级毛片 在线播放| 亚洲国产最新在线播放| av线在线观看网站| 亚洲经典国产精华液单| 亚洲国产av新网站| 美女福利国产在线| 九九爱精品视频在线观看| 大香蕉久久网| 又大又黄又爽视频免费| 午夜福利视频精品| 亚洲精品乱码久久久v下载方式| 精品酒店卫生间| 日韩伦理黄色片| 日日摸夜夜添夜夜爱| 纯流量卡能插随身wifi吗| 大香蕉久久网| 亚洲国产av影院在线观看| 美女内射精品一级片tv| 亚洲国产精品国产精品| 大陆偷拍与自拍| 99热网站在线观看| 夜夜骑夜夜射夜夜干| 亚洲国产毛片av蜜桃av| 永久网站在线| 2018国产大陆天天弄谢| 黄色视频在线播放观看不卡| 人妻 亚洲 视频| 婷婷色麻豆天堂久久| 国产视频首页在线观看| 建设人人有责人人尽责人人享有的| kizo精华| 国产精品成人在线| 黄色怎么调成土黄色| 国产精品蜜桃在线观看| 岛国毛片在线播放| 少妇人妻 视频| 久久午夜福利片| 国产乱人偷精品视频| 美女脱内裤让男人舔精品视频| 成人国产av品久久久| 99视频精品全部免费 在线| 桃花免费在线播放| 国产亚洲精品久久久com| 一级爰片在线观看| 男的添女的下面高潮视频| 男女高潮啪啪啪动态图| 性色av一级| 伊人久久精品亚洲午夜| 26uuu在线亚洲综合色| 精品卡一卡二卡四卡免费| av黄色大香蕉| 男女免费视频国产| 欧美另类一区| 性色av一级| www.av在线官网国产| 人人妻人人澡人人看| 精品国产国语对白av| 国产一区二区在线观看av| 大片免费播放器 马上看| 午夜福利,免费看| 国产高清不卡午夜福利| xxx大片免费视频| 国产成人aa在线观看| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 久久青草综合色| 亚洲精品一区蜜桃| 三级国产精品欧美在线观看| 91久久精品国产一区二区三区| 在线观看www视频免费| 日本黄色日本黄色录像| 久久久久久伊人网av| 国产亚洲精品久久久com| 午夜免费观看性视频| av在线观看视频网站免费| 久久青草综合色| 国产乱人偷精品视频| 久久99一区二区三区| 男女高潮啪啪啪动态图| 香蕉精品网在线| 夫妻性生交免费视频一级片| 五月玫瑰六月丁香| 一级黄片播放器| 免费看光身美女| 极品人妻少妇av视频| 免费观看的影片在线观看| 久久久精品免费免费高清| av在线观看视频网站免费| 黄色毛片三级朝国网站| tube8黄色片| 亚洲av欧美aⅴ国产| 久久午夜综合久久蜜桃| 熟妇人妻不卡中文字幕| 永久网站在线| 欧美性感艳星| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产欧美在线一区| 日韩一区二区三区影片| 日韩熟女老妇一区二区性免费视频| 亚洲综合色惰| 午夜激情av网站| 99国产精品免费福利视频| 少妇人妻久久综合中文| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 七月丁香在线播放| 欧美日韩av久久| 亚洲无线观看免费| 九九久久精品国产亚洲av麻豆| 国产精品嫩草影院av在线观看| av视频免费观看在线观看| 69精品国产乱码久久久| 我要看黄色一级片免费的| 国产午夜精品一二区理论片| 观看av在线不卡| 午夜老司机福利剧场| 日韩电影二区| 亚洲精品日韩av片在线观看| 一本久久精品| kizo精华| 精品一区二区免费观看| 国产精品久久久久久精品电影小说| 亚洲,欧美,日韩| 日日撸夜夜添| 狂野欧美白嫩少妇大欣赏| 看免费成人av毛片| 亚洲国产欧美在线一区| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| 国产国拍精品亚洲av在线观看| 免费看av在线观看网站| 国产有黄有色有爽视频| 久久久久人妻精品一区果冻| 乱人伦中国视频| 男的添女的下面高潮视频| 最近中文字幕高清免费大全6| 在线观看人妻少妇| 18在线观看网站| 毛片一级片免费看久久久久| 午夜视频国产福利| 人人妻人人澡人人看| 国产精品一区www在线观看| 免费大片黄手机在线观看| 免费不卡的大黄色大毛片视频在线观看| 免费黄色在线免费观看| 夫妻性生交免费视频一级片| 亚洲欧洲精品一区二区精品久久久 | 在线免费观看不下载黄p国产| 免费大片黄手机在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆精品久久久久久蜜桃| 国产男女超爽视频在线观看| 中文字幕最新亚洲高清| av黄色大香蕉| 在线观看免费日韩欧美大片 | 这个男人来自地球电影免费观看 | 午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 久久国产精品男人的天堂亚洲 | 91在线精品国自产拍蜜月| 夜夜爽夜夜爽视频| av免费观看日本| 亚洲精品国产av成人精品| 国产伦理片在线播放av一区| 成人国语在线视频| 久久99热6这里只有精品| 老司机影院毛片| 男女国产视频网站| 高清在线视频一区二区三区| 一区二区三区精品91| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 国产成人91sexporn| 久久精品国产亚洲网站| 日本免费在线观看一区| 亚洲情色 制服丝袜| av在线app专区| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 观看av在线不卡| a级毛片黄视频| 成年人免费黄色播放视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩亚洲高清精品| 国产精品 国内视频| 国产免费视频播放在线视频| 国产成人午夜福利电影在线观看| 久久亚洲国产成人精品v| 寂寞人妻少妇视频99o| 日本爱情动作片www.在线观看| 亚洲三级黄色毛片| av免费观看日本| a 毛片基地| 久久女婷五月综合色啪小说| 18禁在线播放成人免费| 成人综合一区亚洲| 精品国产乱码久久久久久小说| 国产成人aa在线观看| 老司机影院毛片| 久久人人爽av亚洲精品天堂| 久久狼人影院| 亚洲综合色惰| 亚洲国产欧美日韩在线播放| 伊人亚洲综合成人网| 99久久综合免费| 91精品一卡2卡3卡4卡| 国产精品国产三级国产av玫瑰| 蜜桃国产av成人99| 久久久久久久久久久久大奶| 一区二区三区免费毛片| a级片在线免费高清观看视频| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品古装| 亚洲国产精品一区二区三区在线| 国产精品久久久久久av不卡| 欧美激情 高清一区二区三区| 777米奇影视久久| 亚洲不卡免费看| av在线观看视频网站免费| 免费日韩欧美在线观看| 久久久国产精品麻豆| 中文字幕最新亚洲高清| 中文字幕久久专区| 日本vs欧美在线观看视频| 建设人人有责人人尽责人人享有的| 国产免费福利视频在线观看| 大片免费播放器 马上看| 少妇的逼水好多| 搡女人真爽免费视频火全软件| 春色校园在线视频观看| 亚洲精品久久成人aⅴ小说 | 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩视频高清一区二区三区二| 涩涩av久久男人的天堂| 亚洲av中文av极速乱| 天美传媒精品一区二区| 久久精品国产鲁丝片午夜精品| 中文天堂在线官网| 亚洲综合色网址| 狂野欧美白嫩少妇大欣赏| 国产成人一区二区在线| 一区二区av电影网| 女性被躁到高潮视频| 国产成人午夜福利电影在线观看| 一区在线观看完整版| videossex国产| 久久国产精品男人的天堂亚洲 | 欧美日韩精品成人综合77777| 午夜老司机福利剧场| 亚洲精品久久成人aⅴ小说 | 日日爽夜夜爽网站| 久久久a久久爽久久v久久| 22中文网久久字幕| 精品国产国语对白av| 人妻制服诱惑在线中文字幕| 久久午夜综合久久蜜桃| 国产av一区二区精品久久| 国产精品人妻久久久影院| 老熟女久久久| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 国产免费一级a男人的天堂| 欧美日韩国产mv在线观看视频| 国产探花极品一区二区| 久久午夜综合久久蜜桃| 国产午夜精品一二区理论片| 九色成人免费人妻av| 久久久a久久爽久久v久久| 亚洲怡红院男人天堂| 美女视频免费永久观看网站| 久久99蜜桃精品久久| 狂野欧美激情性bbbbbb| 欧美精品人与动牲交sv欧美| 女性被躁到高潮视频| 日本免费在线观看一区| 一本久久精品| 99国产综合亚洲精品| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美,日韩| 久久影院123| 国产又色又爽无遮挡免| 97在线视频观看| 男女无遮挡免费网站观看| 爱豆传媒免费全集在线观看| 91成人精品电影| www.色视频.com| 高清毛片免费看| 国产一区二区在线观看日韩| 18禁动态无遮挡网站| 男人操女人黄网站| 热re99久久精品国产66热6| 男女边摸边吃奶| 精品一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 好男人视频免费观看在线| 免费大片黄手机在线观看| 久久久久久久精品精品| 欧美+日韩+精品| 国产精品女同一区二区软件| 在线播放无遮挡| 91久久精品国产一区二区三区| 秋霞伦理黄片| a级毛片免费高清观看在线播放| 日韩一区二区视频免费看| 大片电影免费在线观看免费| 男人添女人高潮全过程视频| 又粗又硬又长又爽又黄的视频| 亚洲欧美色中文字幕在线| 成人国语在线视频| 久久久久久久亚洲中文字幕| 爱豆传媒免费全集在线观看| 天天影视国产精品| 国产男人的电影天堂91| 久久精品熟女亚洲av麻豆精品| 亚洲综合色网址| 精品一区二区三卡| 午夜激情久久久久久久| 国产淫语在线视频| 多毛熟女@视频| 国产av一区二区精品久久| 夫妻午夜视频| 少妇的逼水好多| 伦理电影免费视频| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜添av毛片| 色视频在线一区二区三区| 在线精品无人区一区二区三| 亚洲色图综合在线观看| 国产探花极品一区二区| 黄色毛片三级朝国网站| 黄色配什么色好看| 美女福利国产在线| 精品国产一区二区三区久久久樱花| 一区二区三区免费毛片| 日本免费在线观看一区| 精品一区二区三区视频在线| 男女边摸边吃奶| 一边摸一边做爽爽视频免费| 国产精品无大码| 国产精品久久久久久精品古装| 日本色播在线视频| 亚洲国产成人一精品久久久| 亚洲伊人久久精品综合| 老女人水多毛片| 一本—道久久a久久精品蜜桃钙片| 欧美精品高潮呻吟av久久| xxx大片免费视频| 免费人成在线观看视频色| 精品久久国产蜜桃| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图 | 多毛熟女@视频| 另类亚洲欧美激情| 人人澡人人妻人| 免费看光身美女| 少妇猛男粗大的猛烈进出视频| 久久女婷五月综合色啪小说| 久久国内精品自在自线图片| 久久久久久久久久久久大奶| 国产亚洲一区二区精品| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 美女国产视频在线观看| 国产免费又黄又爽又色| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩另类电影网站| 成人手机av| 国产精品蜜桃在线观看| 日本vs欧美在线观看视频| 少妇精品久久久久久久| 国产男女内射视频| 欧美老熟妇乱子伦牲交| 91在线精品国自产拍蜜月| 久久久国产精品麻豆| 人妻人人澡人人爽人人| 亚洲精品aⅴ在线观看| 少妇人妻精品综合一区二区| 一区在线观看完整版| 自线自在国产av| 中国三级夫妇交换| 国产免费又黄又爽又色| 成人影院久久| 免费大片黄手机在线观看| 一区二区三区免费毛片| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 91aial.com中文字幕在线观看| 中文字幕免费在线视频6| 国产欧美日韩综合在线一区二区| 91成人精品电影| 免费久久久久久久精品成人欧美视频 | 熟妇人妻不卡中文字幕| 久久久久视频综合| 精品一品国产午夜福利视频| 18禁在线播放成人免费| 亚洲成色77777| 一本一本综合久久| 亚洲精品美女久久av网站| 欧美国产精品一级二级三级| 人人妻人人爽人人添夜夜欢视频| 18禁在线无遮挡免费观看视频| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 免费观看a级毛片全部| 日日啪夜夜爽| 97在线视频观看| 午夜福利,免费看| 美女福利国产在线| 三级国产精品欧美在线观看| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 成人综合一区亚洲| 国产精品女同一区二区软件| 久久久久久久久久人人人人人人| 成年人午夜在线观看视频| 亚洲精品,欧美精品| 久久久久久人妻| 久久久久视频综合| 中国美白少妇内射xxxbb| 久久韩国三级中文字幕| 国产免费一级a男人的天堂| a级毛色黄片| 国产国语露脸激情在线看| 国产高清三级在线| 国产亚洲最大av| 久久久久国产网址| 国产精品熟女久久久久浪| 国产精品国产三级专区第一集| 纯流量卡能插随身wifi吗| 久久久久视频综合| 国产伦理片在线播放av一区| 国产免费一区二区三区四区乱码| 22中文网久久字幕| 丝袜在线中文字幕| 爱豆传媒免费全集在线观看| 日韩成人伦理影院| 卡戴珊不雅视频在线播放| 国产有黄有色有爽视频| 麻豆乱淫一区二区| 亚洲伊人久久精品综合| 少妇被粗大的猛进出69影院 | 五月天丁香电影| 国产精品成人在线| 91成人精品电影| 久久ye,这里只有精品| 国产精品.久久久| 成年人午夜在线观看视频| 制服诱惑二区| 国产一区有黄有色的免费视频| a级毛片免费高清观看在线播放| 黄色一级大片看看| 99久国产av精品国产电影| 日本与韩国留学比较| av黄色大香蕉| av天堂久久9| 波野结衣二区三区在线| 日韩 亚洲 欧美在线| 久久精品国产自在天天线| 夜夜爽夜夜爽视频| 啦啦啦中文免费视频观看日本| 永久网站在线| 丝袜在线中文字幕| 欧美精品国产亚洲| 欧美日本中文国产一区发布| 在线 av 中文字幕| 免费观看a级毛片全部| 亚洲第一av免费看| 久久人人爽av亚洲精品天堂| av.在线天堂| 热re99久久精品国产66热6| 国产成人91sexporn| 国产又色又爽无遮挡免| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| 老熟女久久久| 成人影院久久| 九九在线视频观看精品| 韩国av在线不卡| 亚洲综合色惰| 日日摸夜夜添夜夜添av毛片| 日韩免费高清中文字幕av| 久久99精品国语久久久| 国产色爽女视频免费观看| 美女xxoo啪啪120秒动态图| 高清午夜精品一区二区三区| 韩国高清视频一区二区三区| 亚洲精品色激情综合| 91精品三级在线观看| av不卡在线播放| 国产午夜精品久久久久久一区二区三区| 久久久a久久爽久久v久久| 在线观看美女被高潮喷水网站| 亚洲第一区二区三区不卡| 欧美亚洲日本最大视频资源| 久久人人爽人人爽人人片va| 国产成人精品福利久久| a 毛片基地| 亚洲美女视频黄频| 国产精品人妻久久久久久| 丰满少妇做爰视频| 日本与韩国留学比较| 国产一区有黄有色的免费视频| 校园人妻丝袜中文字幕| 纵有疾风起免费观看全集完整版| 熟妇人妻不卡中文字幕| 亚洲欧洲精品一区二区精品久久久 | 日韩制服骚丝袜av| 国产精品人妻久久久久久| 人妻人人澡人人爽人人| 国产一区二区三区av在线| 日韩视频在线欧美| 美女福利国产在线| 亚洲精品久久久久久婷婷小说| 乱码一卡2卡4卡精品| 免费大片18禁| 在线看a的网站| 国产爽快片一区二区三区| 激情五月婷婷亚洲| 99国产综合亚洲精品| 免费黄色在线免费观看| 成人国语在线视频| 日韩,欧美,国产一区二区三区| a 毛片基地| 少妇人妻久久综合中文| 伦精品一区二区三区| 久久久久网色| 国产伦精品一区二区三区视频9| 亚洲人成网站在线观看播放| 午夜福利影视在线免费观看| 日本色播在线视频| 80岁老熟妇乱子伦牲交| 国模一区二区三区四区视频| 中文欧美无线码| 精品视频人人做人人爽| 王馨瑶露胸无遮挡在线观看| 日本与韩国留学比较| 青春草视频在线免费观看| 热re99久久国产66热| 免费人妻精品一区二区三区视频| 欧美一级a爱片免费观看看| 在线观看免费高清a一片| 欧美亚洲 丝袜 人妻 在线| 亚洲第一av免费看| 男男h啪啪无遮挡| 一区二区三区精品91| 中文字幕久久专区| 好男人视频免费观看在线| 日本av免费视频播放| 精品一品国产午夜福利视频| 欧美日韩国产mv在线观看视频| 国产精品一国产av| 日韩欧美一区视频在线观看| 最后的刺客免费高清国语| 欧美人与性动交α欧美精品济南到 | 国产精品国产三级国产av玫瑰| 国产av一区二区精品久久| 免费高清在线观看日韩| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 看免费成人av毛片|