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

    基于協(xié)同進(jìn)化的航空高度單粒子翻轉(zhuǎn)故障生成方法研究

    2019-08-23 05:34王鵬蘆浩劉金枝薛茜男金志威
    現(xiàn)代電子技術(shù) 2019年16期
    關(guān)鍵詞:中子通量

    王鵬 蘆浩 劉金枝 薛茜男 金志威

    摘? 要: 大氣中子是航空高度下造成航空器機(jī)載電子系統(tǒng)翻轉(zhuǎn)故障的主要輻射源。面對(duì)航空單粒子效應(yīng)與航空高度影響關(guān)系,該文提出一種基于協(xié)同進(jìn)化的單粒子翻轉(zhuǎn)故障生成方法。生成的故障數(shù)據(jù)用于模擬實(shí)際中子通量隨高度變化的規(guī)律,為機(jī)載電子設(shè)備單粒子效應(yīng)加固和防護(hù)實(shí)驗(yàn)提供數(shù)據(jù)支持,同時(shí)也可以用于分析航空機(jī)組人員飛行期間所接收的中子輻射劑量。驗(yàn)證結(jié)果表明,該單粒子翻轉(zhuǎn)故障生成方法生成故障數(shù)據(jù)與真實(shí)高度變化下單粒子失效特征吻合,能夠滿足器件航空單粒子效應(yīng)加固測(cè)試的需求。

    關(guān)鍵詞: 協(xié)同進(jìn)化; 航空高度; 翻轉(zhuǎn)故障; 單粒子效應(yīng); 中子通量; 故障生成方法

    中圖分類號(hào): TN383+.3?34; V240.2? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼: A? ? ? ? ? ? ? ? ? ? ?文章編號(hào): 1004?373X(2019)16?0112?05

    0? 引? 言

    隨著航空器機(jī)載電子系統(tǒng)半導(dǎo)體集成度不斷提高,器件的特征尺寸越來(lái)越小,單粒子效應(yīng)會(huì)越來(lái)越嚴(yán)重,逐漸成為影響航空器可靠性和壽命的重要因素之一[1?3]。在北緯30° 30 000英尺高度飛行時(shí),航空器機(jī)載電子系統(tǒng)中應(yīng)用的特征尺寸為0.13 μm的100萬(wàn)門FPGA器件發(fā)生單粒子翻轉(zhuǎn)率為1.85×10-2次/天·器件,相當(dāng)于8×105 FIT,遠(yuǎn)遠(yuǎn)無(wú)法滿足航空器所需的可靠性與安全性要求[4]。單粒子效應(yīng)對(duì)航空器的威脅已經(jīng)引起了FAA(美國(guó)聯(lián)邦航空管理局)、EASA(歐洲航空安全局)和航空工業(yè)方的廣泛重視。

    目前對(duì)于單粒子效應(yīng)的研究手段主要有飛行搭載實(shí)驗(yàn)研究[5?6]、地面實(shí)驗(yàn)裝置模擬研究[7?8]和故障注入研究[9?11]。飛行搭載實(shí)驗(yàn)具有研究結(jié)果直觀、針對(duì)性強(qiáng)等優(yōu)點(diǎn),但其費(fèi)用高、周期長(zhǎng);地面實(shí)驗(yàn)裝置模擬具有輻射源選取靈活的優(yōu)點(diǎn),但是各種客觀因素的限制條件較多;采用故障注入的方法進(jìn)行單粒子效應(yīng)的測(cè)試研究是一種高效且低成本的方式,可以彌補(bǔ)上述方法的不足[12]。

    根據(jù)研究表明,在航空高度和地面層大氣中的單粒子效應(yīng)主要是由質(zhì)子和上層大氣分子互相碰撞產(chǎn)生的次級(jí)中子所引發(fā)的[13];在航空高度,次級(jí)宇宙射線輻射的中子成分貢獻(xiàn)了對(duì)機(jī)組人員約一半的劑量當(dāng)量[14]。如何真實(shí)模擬航空飛行高度下器件所遭受的中子單粒子效應(yīng),成為航空單粒子效應(yīng)防護(hù)實(shí)驗(yàn)首要解決的問(wèn)題。對(duì)于單粒子效應(yīng)的故障注入測(cè)試,大多數(shù)實(shí)驗(yàn)采用逐位翻轉(zhuǎn)的方式,雖然這種方式能夠測(cè)試所有的待測(cè)存儲(chǔ)位,但是并不能較好地模擬實(shí)際航空環(huán)境中的故障狀況,因而需要提出一種有效的試驗(yàn)方法模擬實(shí)際可能發(fā)生的故障。本文針對(duì)航空高度中子通量的數(shù)據(jù)變化特征,根據(jù)實(shí)際航線高度軌跡,進(jìn)行故障數(shù)據(jù)算法研究。生成的模擬故障數(shù)據(jù)符合實(shí)際的中子通量隨高度變化的規(guī)律,能夠滿足故障注入單粒子效應(yīng)研究的輸入需求,使輸入的數(shù)據(jù)更接近真實(shí)的動(dòng)態(tài)特征,為機(jī)載電子設(shè)備單粒子效應(yīng)加固和防護(hù)實(shí)驗(yàn)提供數(shù)據(jù)支持。同時(shí)模擬得出的數(shù)據(jù)能夠?yàn)楹娇諜C(jī)組人員所受的大氣層宇宙射線劑量的研究提供支撐。

    1? 航空高度單粒子效應(yīng)高度模型

    單粒子效應(yīng)是由輻射粒子直接或者間接引起的,航空高度下的輻射源主要是中子,其余有質(zhì)子、π介子、電子和光子,這些次級(jí)粒子是具有高能量的初級(jí)輻射粒子與大氣中的氮、氧原子碰撞而產(chǎn)生的。這些粒子具有較大的能量,如大氣中子的最高能量可達(dá)到1 000 MeV以上[14],而這些高能量的粒子可能會(huì)使航空器機(jī)載電子設(shè)備產(chǎn)生故障。機(jī)載電子設(shè)備單粒子效應(yīng)發(fā)生頻次正比于大氣中子通量(n/cm2·h)與單粒子效應(yīng)截面(cm2)的乘積。

    民用航空器短航線的飛行高度一般在20 000~32 000英尺之間,長(zhǎng)航線的飛行高度則是在26 000~42 000英尺之間,一些公務(wù)機(jī)的飛行高度可以達(dá)到50 000英尺。軍用航空器的飛行高度在民用航空器之上,第三代戰(zhàn)斗機(jī)的升限在55 000~65 000英尺,偵察機(jī)的升限會(huì)更高??紤]到軍用和民用航空器的飛行高度,需要對(duì)這個(gè)飛行高度范圍內(nèi)的中子通量進(jìn)行分析與計(jì)算,能夠在不同高度下隨機(jī)產(chǎn)生在容限范圍內(nèi)的中子通量值,進(jìn)而更加有效地測(cè)評(píng)航空高度下的單粒子翻轉(zhuǎn)率。中子通量與大氣深度密切相關(guān),在海平面上50 000~65 000英尺高度(大氣深度為55~125 g/cm2)達(dá)到極大值。當(dāng)大氣深度超過(guò)150 g/cm2時(shí),中子總通量隨深度變化呈指數(shù)衰減[14]。綜合考慮兩方面因素,選取80 000英尺作為研究高度。通過(guò)對(duì)文獻(xiàn)[15]中子通量數(shù)據(jù)的研究,進(jìn)行了算法處理,生成了高度與大氣中子通量的關(guān)系圖。圖1為0~80 000英尺范圍內(nèi)高度和大氣中子通量的關(guān)系[15]。

    本文依據(jù)之前的分析,在0~80 000英尺的航空高度范圍之間進(jìn)行模擬,以500英尺為節(jié)點(diǎn)依次取值,每個(gè)高度對(duì)應(yīng)生成1 000個(gè)中子通量值,最后生成的模擬數(shù)值可以用矩陣的形式表示。

    [J=J1J2?JN=μ0+σ1,1μ0+σ1,2…μ0+σ1,Sμ1+σ2,1μ1+σ2,2…μ1+σ2,S????μN(yùn)+σN,1μN(yùn)+σN,2…μN(yùn)+σN,S] (1)

    式中:[J1]表示某一高度;N=160表示高度個(gè)數(shù);S=1 000表示在這一高度下對(duì)應(yīng)生成的故障數(shù)據(jù)數(shù)量;μ表示每一高度下對(duì)應(yīng)的理論值;σi,j為某個(gè)高度下的第j個(gè)故障數(shù)據(jù)的誤差值。

    所需的故障注入數(shù)據(jù),在同一高度下均值是一致的,需要求解出誤差值對(duì)數(shù)據(jù)進(jìn)行處理,從而生成故障注入數(shù)組。

    對(duì)于故障注入數(shù)據(jù)進(jìn)行驗(yàn)證時(shí),需要盡可能保證生成的故障數(shù)據(jù)覆蓋到實(shí)際可能產(chǎn)生的各種數(shù)值,即保證故障數(shù)據(jù)的健壯性。對(duì)于每一個(gè)高度,均會(huì)生成1 000個(gè)故障數(shù)據(jù),這1 000個(gè)數(shù)據(jù)需服從一定的正態(tài)分布;同時(shí)每一組數(shù)據(jù)也應(yīng)該服從一定的分布,從而確保模擬產(chǎn)生的故障數(shù)據(jù)能夠覆蓋到實(shí)際可能收到的中子通量值。

    2? 航空高度單粒子翻轉(zhuǎn)故障生成算法流程

    要生成如上160×1 000的故障數(shù)據(jù)矩陣,如果采用常規(guī)的遺傳算法,常規(guī)的思路是將該矩陣轉(zhuǎn)換為一個(gè)一維矩陣,然后將其作為一個(gè)個(gè)體經(jīng)過(guò)一系列的選擇、交叉、變異等算子從而生成符合約束條件的數(shù)據(jù),其算法流程為:

    1) 初始化種群;

    2) 計(jì)算適應(yīng)度;

    3) 選擇、交叉、變異;

    4) 重組種群;

    5) 如果終止條件未達(dá)到,返回步驟2)。

    但是采用這種方式,即使采用實(shí)數(shù)編碼方式,每個(gè)個(gè)體的維度都將高達(dá)160 000,導(dǎo)致算法執(zhí)行過(guò)程中產(chǎn)生巨大的運(yùn)算量,解空間變大,算法的耗時(shí)將會(huì)很長(zhǎng),算法的效率低下。

    上述方法中每一個(gè)個(gè)體都是待解問(wèn)題的一個(gè)候選解,通過(guò)初始化生成大量的個(gè)體組成種群,并模仿生物學(xué)中的選擇、交叉、變異等過(guò)程來(lái)產(chǎn)生一個(gè)最優(yōu)解。

    針對(duì)航空高度單粒子效應(yīng)模擬的問(wèn)題,對(duì)應(yīng)的高度數(shù)據(jù)已經(jīng)確定,相應(yīng)的各個(gè)高度下的數(shù)據(jù)之間有著非常緊密的關(guān)系,每個(gè)高度下對(duì)應(yīng)生成的故障數(shù)據(jù)需要服從正態(tài)分布;同時(shí)為了保證每組數(shù)據(jù)總體分布的有效性,其也需要滿足相應(yīng)的分布。因此,本文基于協(xié)同進(jìn)化理論提出航空高度單粒子翻轉(zhuǎn)故障生成方法。

    首先生成一組航空高度數(shù)據(jù)矩陣,作為單粒子翻轉(zhuǎn)故障數(shù)據(jù)生成的基準(zhǔn),然后通過(guò)擬合生成的大氣中子通量與高度曲線模型,對(duì)應(yīng)提取出在1 h內(nèi)該航空高度下每平方厘米遭受到中子輻射的理論值。依據(jù)高度約束條件進(jìn)行故障數(shù)據(jù)種群初始化,在算法運(yùn)算過(guò)程中不使用變異算子,從而使對(duì)于高度的約束條件在整個(gè)運(yùn)算過(guò)程中始終成立。進(jìn)行適應(yīng)度的計(jì)算,對(duì)160個(gè)高度下的故障數(shù)據(jù)組進(jìn)行分類,對(duì)不符合約束條件的一類應(yīng)用交叉算子,之后在計(jì)算適應(yīng)度并進(jìn)行交叉算子,直到所有個(gè)體均符合約束條件,從而生成符合約束條件的單粒子翻轉(zhuǎn)故障數(shù)據(jù)。本文中的單粒子效應(yīng)故障生成均以單位面積計(jì)算。主要算法流程如下:

    1) 基于高度分布的初始化;

    2) 選擇;

    3) 離散重組;

    4) 如果有不符合約束條件的數(shù)據(jù)組,返回步驟2),否則進(jìn)入下一步;

    5) 故障數(shù)據(jù)隨機(jī)化。

    2.1? 種群初始化

    對(duì)于高度的編碼采用實(shí)數(shù)編碼,將高度數(shù)據(jù)[0,80 000]分成160個(gè)區(qū)間,通過(guò)區(qū)間的最大值來(lái)代表所處區(qū)間,如表1所示。在初始化種群時(shí),首先生成滿足高度需求的160×1的高度序列,利用中子通量和高度的關(guān)系,確定目標(biāo)高度下中子通量的理論值;之后生成翻轉(zhuǎn)故障值時(shí)均以這個(gè)高度所對(duì)應(yīng)的理論值為基準(zhǔn)值;最后進(jìn)行適應(yīng)度計(jì)算,從而更真實(shí)地模擬航空器在對(duì)應(yīng)高度下可能受到的中子通量大小。中子能量理論值如表2所示。

    2.2? 中子通量數(shù)值分布

    為了使生成的數(shù)據(jù)具有多樣性,且能夠?qū)赡艽嬖诘闹凶油恐颠M(jìn)行充分的驗(yàn)證,所以需要確保每個(gè)高度下生成的中子通量值滿足正態(tài)分布。將數(shù)據(jù)誤差分為5個(gè)區(qū)間,每個(gè)高度的誤差序列(σi,1,σi,2,…,σi,s)所屬的區(qū)間分布需服從一定的正態(tài)分布[N′1],N1~N(0,100)。由于Li∈[-100,100],所以需去掉N1中超出邊界的部分,所以各區(qū)間的分布數(shù)量:

    [F1(i)=Φ1(40i-100)-Φ1(40i-140)Φ1(100)-Φ1(-100)·S]? ?(2)

    式中,[Φ1(x)]為N1的概率分布函數(shù)。通過(guò)式(2)計(jì)算出各區(qū)間的分布數(shù),并為其增加一個(gè)±10的容差,各區(qū)間的分布如表3所示。

    2.3? 中子通量數(shù)據(jù)選擇

    選擇算子以適應(yīng)度為標(biāo)準(zhǔn),將N個(gè)高度分成兩類:符合高度約束條件的一類(Lf)和不符合高度約束條件(Lp)的一類。適應(yīng)度的計(jì)算以單個(gè)高度數(shù)據(jù)為單位。每個(gè)高度數(shù)據(jù)的適應(yīng)度的評(píng)價(jià)標(biāo)準(zhǔn)為中子通量值滿足正態(tài)分布。適應(yīng)度分類流程圖如圖2所示。首先判斷每個(gè)高度下的數(shù)據(jù)是否符合約束條件,以此為標(biāo)準(zhǔn)將160個(gè)高度分成兩類。之后按照誤差區(qū)間的不同,將1 000組數(shù)據(jù)依次分為5類,并統(tǒng)計(jì)各類的個(gè)數(shù)。判斷各區(qū)間內(nèi)的數(shù)量是否多于期望值,若多于期望值,將其歸為不符合約束的類別。

    2.4? 離散重組

    選擇算子將N個(gè)高度故障數(shù)據(jù)分成了兩類:Lf和Lp,重組算子只作用于Lp。其作用是生成新的Lp′,使Lp′中的故障數(shù)據(jù)個(gè)數(shù)逐漸減少,使盡可能多的故障數(shù)據(jù)滿足約束條件。重組算子就是將Lp′逐行隨機(jī)重新排序,即離散重組。采用多點(diǎn)交叉的方式既能夠較大程度地改變?cè)械腖p,又可以最大程度上滿足約束條件。

    3? 實(shí)驗(yàn)結(jié)果與分析

    基于本文協(xié)同進(jìn)化的航空高度單粒子翻轉(zhuǎn)故障生成算法,可以生成160×1 000故障數(shù)據(jù)矩陣。以其中30 000英尺高度為例,分析其生成的1 000個(gè)故障數(shù)據(jù)為是否滿足約束條件。通過(guò)對(duì)1 000個(gè)故障數(shù)據(jù)進(jìn)行分析處理,每個(gè)數(shù)據(jù)都代表該高度下每小時(shí)發(fā)生翻轉(zhuǎn)的次數(shù),可以得到其正態(tài)分布圖如圖3所示。

    從圖3可以看出生成的故障數(shù)據(jù)基本符合約束要求。在30 000英尺高度平均每小時(shí)中子通量為1 368.207 n/cm2·h,與標(biāo)準(zhǔn)數(shù)據(jù)對(duì)比的標(biāo)準(zhǔn)偏差為10.274 3,說(shuō)明了該算法生成的故障數(shù)據(jù)能夠很好地滿足需求,符合所需要的目標(biāo)要求。

    4? 實(shí)際航線的單粒子故障生成算法驗(yàn)證

    為了驗(yàn)證基于協(xié)同進(jìn)化的航空高度單粒子效應(yīng)故障數(shù)據(jù)生成的有效性。本文選擇了一條目前在飛的從洛杉磯飛往華盛頓的航線,整段航線時(shí)長(zhǎng)達(dá)4 h 25 min,包括滑行、起飛、爬升、巡航、下降、著陸等全程的飛行軌跡。飛行軌跡航線圖[16]見圖4。

    航線關(guān)鍵數(shù)據(jù)[17]見表4。飛機(jī)從07:01:00離開停機(jī)位,滑行時(shí)間為13 min,之后開始起飛、爬升,至07:43:50,飛機(jī)爬升至35 000英尺的巡航高度;在巡航高度飛機(jī)飛行3 h 23 min后開始下降,11:31:25飛機(jī)著陸完成航行。

    猜你喜歡
    中子通量
    用于乏燃料組件上下端部活化源項(xiàng)分析的截面庫(kù)研究
    控制棒運(yùn)動(dòng)對(duì)反應(yīng)堆功率探測(cè)器的影響機(jī)理研究
    基于氘氚中子源硼中子俘獲治療的中子慢化整形研究
    離散縱標(biāo)六角形節(jié)塊法及CMFD加速研究
    18個(gè)月?lián)Q料對(duì)CPR1000反應(yīng)堆壓力容器輻照監(jiān)督的影響
    裂變室輸出信號(hào)數(shù)字化處理的仿真研究
    修正快中子通量以提高碳氧測(cè)量精度的研究
    線性回歸方法在核數(shù)據(jù)處理中的應(yīng)用
    煤炭?jī)?nèi)中子通量與元素含量關(guān)系
    BNCT醫(yī)院中子照射器輻射場(chǎng)特性參數(shù)初步測(cè)量
    大片电影免费在线观看免费| av一本久久久久| 天天躁日日操中文字幕| 人妻系列 视频| 日韩视频在线欧美| 成人亚洲欧美一区二区av| 一区二区三区精品91| 日韩av免费高清视频| 日本一二三区视频观看| 久久久a久久爽久久v久久| 99久久九九国产精品国产免费| 久久久午夜欧美精品| 五月玫瑰六月丁香| 久久久久精品性色| 久热这里只有精品99| 日本黄大片高清| 菩萨蛮人人尽说江南好唐韦庄| 天美传媒精品一区二区| 干丝袜人妻中文字幕| 精品久久久久久久人妻蜜臀av| 国产成人精品久久久久久| 日本午夜av视频| freevideosex欧美| 身体一侧抽搐| 三级经典国产精品| 日韩国内少妇激情av| 日日摸夜夜添夜夜添av毛片| 国产探花极品一区二区| 亚洲成人av在线免费| 久久久亚洲精品成人影院| 欧美日韩在线观看h| 美女国产视频在线观看| 国产成人福利小说| 综合色丁香网| 插阴视频在线观看视频| 亚洲人与动物交配视频| 国产久久久一区二区三区| 国产探花在线观看一区二区| 91午夜精品亚洲一区二区三区| 又粗又硬又长又爽又黄的视频| 插阴视频在线观看视频| 成人综合一区亚洲| 精品人妻偷拍中文字幕| 18禁在线播放成人免费| 蜜桃亚洲精品一区二区三区| 网址你懂的国产日韩在线| 精华霜和精华液先用哪个| 超碰av人人做人人爽久久| 国产乱来视频区| 看黄色毛片网站| 国产黄色免费在线视频| 国产v大片淫在线免费观看| 国产人妻一区二区三区在| 国产一区有黄有色的免费视频| 在线天堂最新版资源| 精品一区二区免费观看| 国精品久久久久久国模美| 国产人妻一区二区三区在| 2021天堂中文幕一二区在线观| 在线免费观看不下载黄p国产| 亚洲高清免费不卡视频| 青春草亚洲视频在线观看| 一级毛片久久久久久久久女| 卡戴珊不雅视频在线播放| 日韩亚洲欧美综合| 国产精品嫩草影院av在线观看| 黄色视频在线播放观看不卡| 欧美日韩亚洲高清精品| 久久久久久久亚洲中文字幕| 91狼人影院| 亚洲国产高清在线一区二区三| 91aial.com中文字幕在线观看| 日韩制服骚丝袜av| 有码 亚洲区| 少妇猛男粗大的猛烈进出视频 | 成人欧美大片| 夫妻性生交免费视频一级片| 2021天堂中文幕一二区在线观| 综合色丁香网| 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 精品久久久噜噜| 国内精品美女久久久久久| 亚洲精品日韩在线中文字幕| 国产老妇女一区| 亚洲自拍偷在线| 久久99精品国语久久久| 亚洲精品久久久久久婷婷小说| 最新中文字幕久久久久| 国产精品无大码| 日韩三级伦理在线观看| 亚洲精品久久午夜乱码| 国内精品宾馆在线| 亚洲成人一二三区av| 精品一区二区三卡| 网址你懂的国产日韩在线| 神马国产精品三级电影在线观看| 草草在线视频免费看| 欧美国产精品一级二级三级 | .国产精品久久| 自拍偷自拍亚洲精品老妇| 又大又黄又爽视频免费| 欧美日韩国产mv在线观看视频 | 天堂中文最新版在线下载 | 中文字幕免费在线视频6| 人人妻人人澡人人爽人人夜夜| 亚洲国产av新网站| 男人爽女人下面视频在线观看| a级一级毛片免费在线观看| 免费高清在线观看视频在线观看| 偷拍熟女少妇极品色| 晚上一个人看的免费电影| 亚洲经典国产精华液单| 高清日韩中文字幕在线| 精品少妇久久久久久888优播| 97热精品久久久久久| 色婷婷久久久亚洲欧美| 高清在线视频一区二区三区| 亚洲av成人精品一区久久| 欧美+日韩+精品| 人体艺术视频欧美日本| 亚洲精品,欧美精品| 久久久久久伊人网av| 久久综合国产亚洲精品| 亚洲欧美中文字幕日韩二区| 亚洲天堂国产精品一区在线| 不卡视频在线观看欧美| 久久热精品热| 一本色道久久久久久精品综合| 一二三四中文在线观看免费高清| 少妇熟女欧美另类| 久久亚洲国产成人精品v| 一区二区三区精品91| 亚洲aⅴ乱码一区二区在线播放| 国产野战对白在线观看| 天天操日日干夜夜撸| 国产日韩一区二区三区精品不卡| 久久精品国产a三级三级三级| 精品免费久久久久久久清纯 | 秋霞在线观看毛片| 精品国产一区二区三区久久久樱花| 大码成人一级视频| 亚洲成人一二三区av| 亚洲人成77777在线视频| 爱豆传媒免费全集在线观看| 国产亚洲一区二区精品| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 亚洲欧美激情在线| 精品国产露脸久久av麻豆| 一边亲一边摸免费视频| 晚上一个人看的免费电影| 精品国产乱码久久久久久男人| 夫妻午夜视频| www.自偷自拍.com| av网站免费在线观看视频| 午夜免费鲁丝| 国产免费又黄又爽又色| 午夜av观看不卡| 9191精品国产免费久久| kizo精华| 国产精品亚洲av一区麻豆 | 操出白浆在线播放| 人妻人人澡人人爽人人| 久热这里只有精品99| 波多野结衣一区麻豆| av国产久精品久网站免费入址| 精品一区二区三卡| 69精品国产乱码久久久| 美女国产高潮福利片在线看| av有码第一页| 搡老乐熟女国产| 亚洲美女搞黄在线观看| 欧美在线一区亚洲| 中文字幕av电影在线播放| 国产一区有黄有色的免费视频| 国产极品天堂在线| 一区二区三区乱码不卡18| 一区二区日韩欧美中文字幕| 精品一区二区免费观看| 国产乱来视频区| 夫妻性生交免费视频一级片| 人人澡人人妻人| 日本91视频免费播放| 亚洲一区二区三区欧美精品| 咕卡用的链子| 另类亚洲欧美激情| 日日撸夜夜添| 亚洲精品在线美女| 中文天堂在线官网| av卡一久久| h视频一区二区三区| 哪个播放器可以免费观看大片| av电影中文网址| 人妻 亚洲 视频| 人人妻人人澡人人看| 下体分泌物呈黄色| 午夜福利免费观看在线| 欧美久久黑人一区二区| 亚洲第一av免费看| 日韩伦理黄色片| 曰老女人黄片| 精品国产露脸久久av麻豆| av天堂久久9| 国产乱人偷精品视频| 一级片免费观看大全| 卡戴珊不雅视频在线播放| 精品一品国产午夜福利视频| 最新在线观看一区二区三区 | av在线app专区| 亚洲精品日韩在线中文字幕| 99久国产av精品国产电影| 宅男免费午夜| 天天躁夜夜躁狠狠久久av| 成年美女黄网站色视频大全免费| 国产高清不卡午夜福利| 色吧在线观看| 男人爽女人下面视频在线观看| 一本一本久久a久久精品综合妖精| 看非洲黑人一级黄片| 午夜福利视频精品| 日韩制服丝袜自拍偷拍| 51午夜福利影视在线观看| 一级毛片电影观看| 激情视频va一区二区三区| 人人妻人人澡人人爽人人夜夜| 色吧在线观看| 亚洲精品成人av观看孕妇| 赤兔流量卡办理| 亚洲精品aⅴ在线观看| avwww免费| 免费高清在线观看视频在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲伊人色综图| 亚洲成色77777| 啦啦啦在线观看免费高清www| 一级a爱视频在线免费观看| 亚洲,欧美,日韩| 精品国产一区二区久久| 久久韩国三级中文字幕| 国产一区有黄有色的免费视频| 电影成人av| 欧美少妇被猛烈插入视频| 精品国产超薄肉色丝袜足j| 午夜福利在线免费观看网站| 久久久久国产一级毛片高清牌| 妹子高潮喷水视频| 国产精品av久久久久免费| 熟女少妇亚洲综合色aaa.| 悠悠久久av| 欧美日韩亚洲综合一区二区三区_| 日韩精品有码人妻一区| 午夜福利影视在线免费观看| 少妇精品久久久久久久| 最近中文字幕2019免费版| 操美女的视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| 午夜福利网站1000一区二区三区| 多毛熟女@视频| 天堂俺去俺来也www色官网| 精品人妻熟女毛片av久久网站| 国产极品粉嫩免费观看在线| 精品少妇内射三级| 国产精品久久久av美女十八| 欧美精品人与动牲交sv欧美| 精品一区在线观看国产| 韩国高清视频一区二区三区| 欧美激情极品国产一区二区三区| 日本猛色少妇xxxxx猛交久久| 伊人久久国产一区二区| 久久久久精品久久久久真实原创| 丝袜在线中文字幕| 两性夫妻黄色片| 国产 精品1| 国产高清不卡午夜福利| 精品少妇内射三级| 伦理电影大哥的女人| av天堂久久9| 久久狼人影院| 色综合欧美亚洲国产小说| 99久久人妻综合| 91国产中文字幕| 在线 av 中文字幕| 亚洲国产精品一区二区三区在线| 丝袜美足系列| 亚洲久久久国产精品| 欧美中文综合在线视频| 精品一区二区三卡| 成人三级做爰电影| 欧美在线黄色| 国产一区亚洲一区在线观看| 乱人伦中国视频| 中文欧美无线码| 国产精品免费视频内射| a级毛片在线看网站| 午夜福利影视在线免费观看| 美女脱内裤让男人舔精品视频| 看免费av毛片| 91精品三级在线观看| 免费av中文字幕在线| 一个人免费看片子| 女人精品久久久久毛片| 97在线人人人人妻| 日韩精品有码人妻一区| 熟妇人妻不卡中文字幕| 看非洲黑人一级黄片| 天堂俺去俺来也www色官网| 久久免费观看电影| av不卡在线播放| 色网站视频免费| 卡戴珊不雅视频在线播放| 99国产精品免费福利视频| 成人影院久久| kizo精华| 亚洲av电影在线观看一区二区三区| 久久狼人影院| 日本猛色少妇xxxxx猛交久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人欧美在线观看 | 美女中出高潮动态图| 中国三级夫妇交换| 黑人巨大精品欧美一区二区蜜桃| 亚洲av成人精品一二三区| 国产精品久久久久久精品古装| 男女边摸边吃奶| 中文乱码字字幕精品一区二区三区| 香蕉丝袜av| 99九九在线精品视频| 国产精品.久久久| 色视频在线一区二区三区| 熟女av电影| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| 国产片内射在线| 日韩视频在线欧美| 国产激情久久老熟女| 日本91视频免费播放| 操美女的视频在线观看| av在线老鸭窝| 精品久久蜜臀av无| 男人添女人高潮全过程视频| 午夜久久久在线观看| 国产色婷婷99| 一级片'在线观看视频| 精品亚洲成a人片在线观看| 午夜老司机福利片| 亚洲精品一区蜜桃| 日日撸夜夜添| 午夜影院在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人免费av在线播放| 七月丁香在线播放| 国产成人系列免费观看| 一边亲一边摸免费视频| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 久久精品熟女亚洲av麻豆精品| 自拍欧美九色日韩亚洲蝌蚪91| 国产一区有黄有色的免费视频| 少妇的丰满在线观看| av在线app专区| av片东京热男人的天堂| 日本91视频免费播放| 成人国产av品久久久| 国产精品成人在线| 青春草国产在线视频| 日韩av不卡免费在线播放| 亚洲国产精品成人久久小说| 久久天躁狠狠躁夜夜2o2o | 免费看不卡的av| 国产亚洲最大av| 午夜福利一区二区在线看| 男女高潮啪啪啪动态图| 国产成人一区二区在线| 日本爱情动作片www.在线观看| 国产在线视频一区二区| 国产又色又爽无遮挡免| 国产视频首页在线观看| 亚洲男人天堂网一区| 亚洲国产精品一区三区| 国产一区二区三区综合在线观看| 日日爽夜夜爽网站| 国产高清不卡午夜福利| 久久精品亚洲熟妇少妇任你| 一级片免费观看大全| 97精品久久久久久久久久精品| 这个男人来自地球电影免费观看 | 国产精品秋霞免费鲁丝片| 伊人久久大香线蕉亚洲五| 国产国语露脸激情在线看| 日韩制服丝袜自拍偷拍| 亚洲一级一片aⅴ在线观看| 国产深夜福利视频在线观看| 亚洲av男天堂| 久久国产精品男人的天堂亚洲| 精品午夜福利在线看| 男人添女人高潮全过程视频| 亚洲国产中文字幕在线视频| 日韩欧美一区视频在线观看| 青春草亚洲视频在线观看| 人人妻人人爽人人添夜夜欢视频| 精品少妇久久久久久888优播| 99re6热这里在线精品视频| xxxhd国产人妻xxx| 亚洲人成77777在线视频| 人妻人人澡人人爽人人| 交换朋友夫妻互换小说| 日韩中文字幕欧美一区二区 | 久久久久久久精品精品| 性高湖久久久久久久久免费观看| 咕卡用的链子| 纵有疾风起免费观看全集完整版| 欧美成人午夜精品| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品在线美女| bbb黄色大片| 日日啪夜夜爽| 久久久久视频综合| 国产免费又黄又爽又色| 国产亚洲最大av| 99re6热这里在线精品视频| 中文字幕av电影在线播放| 午夜福利在线免费观看网站| 欧美国产精品va在线观看不卡| 国产成人精品久久二区二区91 | 亚洲国产精品一区二区三区在线| 亚洲精品第二区| 如日韩欧美国产精品一区二区三区| 欧美 日韩 精品 国产| 久久这里只有精品19| 国产野战对白在线观看| 男女边吃奶边做爰视频| 精品少妇黑人巨大在线播放| 大片免费播放器 马上看| 肉色欧美久久久久久久蜜桃| 日韩 亚洲 欧美在线| 免费人妻精品一区二区三区视频| 久久久久久人妻| 18禁裸乳无遮挡动漫免费视频| 久久国产精品男人的天堂亚洲| 欧美少妇被猛烈插入视频| 人人妻人人爽人人添夜夜欢视频| 精品少妇内射三级| 久久99一区二区三区| 97精品久久久久久久久久精品| 中文欧美无线码| 男人舔女人的私密视频| 国产一区亚洲一区在线观看| 亚洲国产看品久久| 少妇的丰满在线观看| 秋霞伦理黄片| 校园人妻丝袜中文字幕| 人体艺术视频欧美日本| 亚洲欧美精品综合一区二区三区| 国产一区二区三区av在线| www.自偷自拍.com| 午夜激情久久久久久久| 久久国产精品男人的天堂亚洲| 伦理电影免费视频| av有码第一页| 久久久久网色| 国产精品无大码| 天天躁夜夜躁狠狠躁躁| 亚洲一级一片aⅴ在线观看| 久久精品aⅴ一区二区三区四区| 天堂8中文在线网| 成人免费观看视频高清| 日本欧美国产在线视频| 亚洲精品,欧美精品| 亚洲美女搞黄在线观看| 一边亲一边摸免费视频| 精品人妻熟女毛片av久久网站| 99九九在线精品视频| 人成视频在线观看免费观看| 超碰97精品在线观看| 一区二区av电影网| 搡老岳熟女国产| 国产片内射在线| 黄色视频不卡| 在线 av 中文字幕| 午夜av观看不卡| 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 看免费成人av毛片| 亚洲自偷自拍图片 自拍| 婷婷色麻豆天堂久久| 999久久久国产精品视频| 久久久久视频综合| 国产毛片在线视频| 一本久久精品| 最近最新中文字幕大全免费视频 | 国产精品免费视频内射| 欧美人与善性xxx| 人成视频在线观看免费观看| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 色播在线永久视频| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 亚洲第一区二区三区不卡| 色婷婷久久久亚洲欧美| 成年美女黄网站色视频大全免费| 女性被躁到高潮视频| 香蕉丝袜av| 欧美国产精品va在线观看不卡| 国产福利在线免费观看视频| 9热在线视频观看99| 国产免费福利视频在线观看| 老司机深夜福利视频在线观看 | 老司机在亚洲福利影院| 男女边摸边吃奶| 97精品久久久久久久久久精品| 色精品久久人妻99蜜桃| 九色亚洲精品在线播放| 美女大奶头黄色视频| 国产精品熟女久久久久浪| 亚洲欧美精品自产自拍| 日韩免费高清中文字幕av| 一级毛片黄色毛片免费观看视频| 久久精品国产综合久久久| 看十八女毛片水多多多| 亚洲国产最新在线播放| 一级毛片黄色毛片免费观看视频| 欧美 亚洲 国产 日韩一| 精品免费久久久久久久清纯 | 色视频在线一区二区三区| 老汉色av国产亚洲站长工具| 啦啦啦在线免费观看视频4| 看十八女毛片水多多多| 9191精品国产免费久久| 叶爱在线成人免费视频播放| 午夜91福利影院| 国产有黄有色有爽视频| 欧美成人午夜精品| 最新在线观看一区二区三区 | 人体艺术视频欧美日本| 建设人人有责人人尽责人人享有的| 少妇人妻久久综合中文| 亚洲,欧美精品.| 亚洲国产最新在线播放| 国产黄频视频在线观看| 丝袜喷水一区| 久久久久久久精品精品| 天天操日日干夜夜撸| 成年av动漫网址| 多毛熟女@视频| 午夜福利免费观看在线| 国产亚洲av高清不卡| 桃花免费在线播放| 男女床上黄色一级片免费看| 国产精品熟女久久久久浪| 亚洲一卡2卡3卡4卡5卡精品中文| 国产免费又黄又爽又色| 日日摸夜夜添夜夜爱| 伊人亚洲综合成人网| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 国产有黄有色有爽视频| 999精品在线视频| 大片电影免费在线观看免费| 成人免费观看视频高清| 一级黄片播放器| 十分钟在线观看高清视频www| 人人妻人人爽人人添夜夜欢视频| 国产精品嫩草影院av在线观看| 如日韩欧美国产精品一区二区三区| 国产一区有黄有色的免费视频| 亚洲精品久久午夜乱码| av网站免费在线观看视频| 亚洲专区中文字幕在线 | 欧美激情极品国产一区二区三区| 国产精品一区二区精品视频观看| 色婷婷av一区二区三区视频| 欧美黑人欧美精品刺激| av女优亚洲男人天堂| 国产精品秋霞免费鲁丝片| 亚洲成人国产一区在线观看 | 男人添女人高潮全过程视频| 91精品国产国语对白视频| 日韩一本色道免费dvd| 男女午夜视频在线观看| 久久久久久久精品精品| 国产精品 国内视频| 国产日韩欧美亚洲二区| 校园人妻丝袜中文字幕| 亚洲国产欧美一区二区综合| 精品人妻一区二区三区麻豆| 青青草视频在线视频观看| 亚洲精品aⅴ在线观看| 亚洲图色成人| 天堂8中文在线网| 交换朋友夫妻互换小说| 欧美成人精品欧美一级黄| 欧美国产精品一级二级三级| 叶爱在线成人免费视频播放| 飞空精品影院首页| 多毛熟女@视频| 精品国产乱码久久久久久男人| 欧美 日韩 精品 国产| 一级毛片电影观看| 欧美激情高清一区二区三区 | 国产一区亚洲一区在线观看| av有码第一页| 建设人人有责人人尽责人人享有的| 中文精品一卡2卡3卡4更新| 国产精品.久久久| 青春草国产在线视频| av网站在线播放免费| 国产高清不卡午夜福利| 丝袜在线中文字幕| 男女高潮啪啪啪动态图| 欧美xxⅹ黑人| 亚洲图色成人| 精品少妇久久久久久888优播|