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

    基于遺傳算法的分?jǐn)?shù)階鋰電池模型參數(shù)辨識(shí)

    2022-12-01 06:49:24吳俊楊維明姚靜黃秋安
    關(guān)鍵詞:階躍適應(yīng)度遺傳算法

    吳俊,楊維明,姚靜,黃秋安

    (1.湖北大學(xué)計(jì)算機(jī)與信息工程學(xué)院, 湖北 武漢 430062; 2.上海大學(xué)可持續(xù)能源研究院, 上海 200444)

    0 引言

    鋰離子電池由于其能量密度高,使用壽命長(zhǎng),記憶效應(yīng)低等特點(diǎn),通常被用于便攜式電子設(shè)備(如手機(jī)、藍(lán)牙耳機(jī))和電動(dòng)汽車(chē)等領(lǐng)域[1-2],并且在軍事和航空航天應(yīng)用中日益流行.鋰離子電池在其充放電過(guò)程中,鋰離子在正負(fù)電極上完成脫嵌和嵌入,在電解質(zhì)中完成轉(zhuǎn)移.由于傳統(tǒng)鋰電池使用液體電解質(zhì),而液體電解質(zhì)具有易燃、易爆、易揮發(fā)的特性,在受到外部撞擊或者其他極端情況時(shí),會(huì)發(fā)生電池短路,繼而起火燃燒,無(wú)法提供足夠的安全保障.最近的特斯拉汽車(chē)起火事件表明,必須對(duì)電池的運(yùn)行狀態(tài)進(jìn)行監(jiān)控和管理.

    單純使用SOCSOH等外特性參數(shù)描述電池,不足以揭示電池運(yùn)行狀態(tài).尤其不足以探究電池電氣特性的變化,當(dāng)電池內(nèi)阻變得很大時(shí),根據(jù)焦耳定律,容易造成電池局部過(guò)熱,進(jìn)而導(dǎo)致電池爆炸起火等重大安全事故.因此,要完成對(duì)電池運(yùn)行狀態(tài)的監(jiān)控管理,必須對(duì)電池的阻抗特性進(jìn)行描述.

    本文中主要工作是選用分?jǐn)?shù)階模型作為鋰電池的等效電路,主要完成三個(gè)工作:1)利用分?jǐn)?shù)階定義,給出該分?jǐn)?shù)階模型的任意序列輸入條件下的數(shù)值解;2)完成在該模型框架下的階躍響應(yīng);3)利用遺傳算法完成對(duì)該模型電路參數(shù)的辨識(shí)仿真,并分析其誤差.

    1 分?jǐn)?shù)階理論及改進(jìn)Randle模型的參數(shù)提取

    分?jǐn)?shù)階理論有三種定義,這里僅介紹一種常用的數(shù)值逼近方法,即Grüwald-Letnikov,簡(jiǎn)稱(chēng)GL定義.假設(shè)函數(shù)滿(mǎn)足處處連續(xù)且可導(dǎo)條件,根據(jù)整數(shù)階的Hospital法則,可擴(kuò)展得到其分?jǐn)?shù)階的導(dǎo)數(shù)定義為[3-5]:

    (1)

    說(shuō)明如下:

    1)「?表示向上取整,h表示步長(zhǎng),類(lèi)似整數(shù)階導(dǎo)數(shù)分母,越小越好,α表示分?jǐn)?shù)階階數(shù);

    2)t0表示時(shí)間起點(diǎn),即在該定義下,當(dāng)前函數(shù)值涉及從t0開(kāi)始到t之間所有時(shí)間點(diǎn)的函數(shù)值;

    GL定義下,其分?jǐn)?shù)階積分的拉氏變換見(jiàn)公式(2):

    (2)

    在經(jīng)典等效電路模型中,一階RC模型試圖模擬歐姆電阻和電荷轉(zhuǎn)移反應(yīng),對(duì)應(yīng)于高頻尾線和低頻直線,但沒(méi)有描述中頻范圍內(nèi)的擴(kuò)散行為.對(duì)于二階RC模型,可以模擬歐姆電阻和擴(kuò)散行為,但不能考慮電荷轉(zhuǎn)移過(guò)程.Yu Zhihao等[6]使用第三階或更高階的RC串聯(lián)模型,所有三個(gè)過(guò)程都可以產(chǎn)生.然而,這些模型有大量的參數(shù),不僅使與計(jì)算負(fù)擔(dān)相關(guān)的模型數(shù)學(xué)結(jié)構(gòu)復(fù)雜化,還增加系統(tǒng)校準(zhǔn)的工作量和過(guò)擬合的風(fēng)險(xiǎn).此外,正如Westerhoff等[7]所證明的,即使使用5個(gè)RC串聯(lián)網(wǎng)絡(luò),該模型也不能非常準(zhǔn)確地捕捉零相位.

    為了解決上述問(wèn)題,分?jǐn)?shù)階微積分(FOC)最近被用于鋰離子電池的建模,Wang J L等[8]用CPE單元替換一階RC模型并聯(lián)塊中的理想電容器,得到分?jǐn)?shù)階模型.這些分?jǐn)?shù)階等效電路模型可能包含一個(gè)或多個(gè)CPE元件.通過(guò)簡(jiǎn)單地將一階RC模型中的理想電容替換為分?jǐn)?shù)階單元,得到的分?jǐn)?shù)階模型FOM,如圖1所示.利用分?jǐn)?shù)階導(dǎo)數(shù),發(fā)現(xiàn)該模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的擬合有較好的精度,該分?jǐn)?shù)階模型的阻抗表達(dá)式為:

    Z(s)=R0+Zct(s)+Zd(s)

    (3)

    說(shuō)明如下:

    1)R0表示電解質(zhì)的歐姆電阻和連接電阻;

    2)Zct(s)表示電荷轉(zhuǎn)移電阻,通常由一個(gè)電阻Rp與CPE并聯(lián)組成,其中CPE表示電解質(zhì)和電極表面形成的非理想電容,也稱(chēng)為常相位單元[9-11];

    3)Zd(s)表示鋰離子在低頻條件下的擴(kuò)散電阻,被稱(chēng)為Warburg阻抗.

    常相位單元CPE的頻域表達(dá)式為:CPE(s)=1/Q1sα.其中,Q1是一個(gè)常數(shù),量綱是F·s1-α;α∈(0,1]表示非理想電容的分?jǐn)?shù)階階數(shù),當(dāng)α=1時(shí),CPE就退化成標(biāo)準(zhǔn)電容,這里τ=(RpQ1)-α,公式(3)中的電荷轉(zhuǎn)移電阻可表示為:

    (4)

    對(duì)于阻抗譜的實(shí)際測(cè)量,頻率范圍通常受到測(cè)試設(shè)備和物體特性的限制.因此有效的頻率范圍一般設(shè)置為:(100 mHz,100 kHz).本文中著重討論實(shí)際應(yīng)用情況下,半無(wú)限擴(kuò)散條件下的Warburg阻抗,即:

    (5)

    待辨識(shí)參數(shù)為:θ={R0,Rp,Q1,α,Q2}.

    2 改進(jìn)Randle模型下的分?jǐn)?shù)階響應(yīng)

    (6)

    此時(shí),除了上面提到的利用定義法求解分?jǐn)?shù)階響應(yīng)外,還可以借助文獻(xiàn)[12-13],利用反拉普拉斯變換,給出該階躍響應(yīng)的解析解,見(jiàn)公式(7)~(8),這里的α∈(0,1].文獻(xiàn)[12-13]中已經(jīng)證明,該解析解與定義法等價(jià).

    (7)

    (8)

    如圖2為改進(jìn)Randle模型下的階躍響應(yīng).階躍電流為1 A,橫軸表示時(shí)間,縱軸表示響應(yīng)電壓,參數(shù)θ={7.37,6.243 ,3.073 ,0.782 ,608.1},即參照下文中表1中的預(yù)設(shè)參數(shù)θ值取值.

    3 參數(shù)辨識(shí)及其結(jié)果分析

    辨識(shí)的基本步驟包括實(shí)驗(yàn)設(shè)計(jì)、數(shù)據(jù)采集、數(shù)據(jù)預(yù)處理、結(jié)構(gòu)辨識(shí)、參數(shù)估計(jì)、模型驗(yàn)證等.實(shí)驗(yàn)設(shè)計(jì)主要包括激勵(lì)信號(hào)的選擇,采樣周期的確定,激勵(lì)信號(hào)數(shù)據(jù)長(zhǎng)度,離線辨識(shí)還是在線辨識(shí)等,主要目的是使采集的輸入輸出數(shù)據(jù)盡可能多地包涵系統(tǒng)的內(nèi)在特性.因此,本仿真實(shí)驗(yàn)采用階躍電流和正弦電流作激勵(lì),希望盡可能與實(shí)際獲得的實(shí)驗(yàn)數(shù)據(jù)相仿或相近.

    參數(shù)辨識(shí)方法有很多,從信號(hào)的實(shí)時(shí)性看,可以分為兩大類(lèi),一是在線辨識(shí),二是離線辨識(shí).本文中采用離線辨識(shí)方法,利用遺傳算法,作參數(shù)優(yōu)化,得到最終的參數(shù)辨識(shí)結(jié)果.實(shí)際上除了一次完成算法外,其他無(wú)論哪種方法,其辨識(shí)流程大體均如圖3所示.

    如果用θ表示分?jǐn)?shù)階系統(tǒng)的參數(shù),用F(θ)表示待辨識(shí)模型的數(shù)學(xué)映射,當(dāng)輸入電流為i(t)時(shí),可以得到辨識(shí)模型計(jì)算值um和鋰電池模型的響應(yīng)電壓值u,計(jì)算目標(biāo)函數(shù)Jθ的值:

    Jθ=(u-um)2

    (9)

    更新并調(diào)整參數(shù)θ的值,使得目標(biāo)函數(shù)Jθ朝變小的方向靠近,當(dāng)輸出誤差函數(shù)小到某一指定精度ε,如當(dāng)Jθ<ε時(shí)即認(rèn)為辨識(shí)成功.可認(rèn)為在指定精度ε下,可以用辨識(shí)模型代替鋰電池系統(tǒng).

    3.1 遺傳算法原理遺傳算法[14-16]是美國(guó)教授John Hollond于上世紀(jì)60年代提出的一種受自然選擇過(guò)程啟發(fā)的元啟發(fā)式算法,屬于進(jìn)化智能算法的一種,通常依賴(lài)于生物啟發(fā)的算子,如選擇、交叉、變異等生成高質(zhì)量的優(yōu)化和搜索問(wèn)題的解決方案.遺傳算法本質(zhì)上是一種基于達(dá)爾文的“物競(jìng)天擇、適者生存”原則的仿生算法.目標(biāo)生物種群在捕食者和環(huán)境危害威脅等外在條件下,不能適應(yīng)的弱小生物個(gè)體被消滅,更能適應(yīng)生存環(huán)境的優(yōu)秀個(gè)體被保留.通過(guò)優(yōu)勝劣汰篩選后的生物成員往往比其他成員有更好的生存機(jī)會(huì).它們更有可能適應(yīng)進(jìn)化條件,它們產(chǎn)生的后代可能會(huì)遺傳并繼承他們的特征,學(xué)習(xí)和模仿它們的技能,從而產(chǎn)生更健康的后代.此外,基因突變?cè)谖锓N成員中隨機(jī)發(fā)生,其中一些突變可能會(huì)提高健康個(gè)體及其進(jìn)化后代長(zhǎng)期存活的機(jī)會(huì).最后存活下來(lái)得到的個(gè)體就是優(yōu)化問(wèn)題的最優(yōu)解.

    標(biāo)準(zhǔn)遺傳算法從隨機(jī)生成的可能解群體(個(gè)體)開(kāi)始.首先計(jì)算個(gè)體的適應(yīng)度,根據(jù)個(gè)體的適應(yīng)度值選擇部分個(gè)體作為父母,產(chǎn)生父群體.通過(guò)將交叉算子應(yīng)用于父群體,然后將變異算子應(yīng)用于它們的后代,產(chǎn)生了可能解的新群體(子群體).重復(fù)用新子代替換雙親(舊個(gè)體)進(jìn)行迭代,直到滿(mǎn)足停止標(biāo)準(zhǔn).

    本文中采用遺傳算法的流程圖如圖4所示.首先定義電池參數(shù)θ,隨機(jī)生成滿(mǎn)足要求的初始種群,計(jì)算初始種群的目標(biāo)函數(shù)值,通過(guò)遺傳算子(選擇、交叉、變異)得到新的候選解種群,重復(fù)計(jì)算新種群中目標(biāo)函數(shù)的值.當(dāng)某個(gè)體目標(biāo)函數(shù)值小于指定精度時(shí),停止迭代,該個(gè)體即為最優(yōu)個(gè)體輸出,否則繼續(xù)迭代.需要注意目標(biāo)函數(shù)值與適應(yīng)度一一對(duì)應(yīng),目標(biāo)函數(shù)需要作適應(yīng)度轉(zhuǎn)換,才能輸出該個(gè)體的適應(yīng)度.

    3.2 辨識(shí)結(jié)果分析本文中采用遺傳算法作為辨識(shí)算法,為了驗(yàn)證辨識(shí)算法的有效性,假定電池的既定參數(shù)如表頭所示,在激勵(lì)電流為階躍信號(hào)輸入的條件下,使用遺傳算法,取較好的4次結(jié)果展示在表1中.由于遺傳算法的收斂結(jié)果有一定波動(dòng)性,每次運(yùn)行結(jié)果都不盡相同,因此表中僅列出4組辨識(shí)結(jié)果.如表1中誤差均值MAE所示,每個(gè)點(diǎn)對(duì)應(yīng)的誤差均值精度達(dá)10-5,比原數(shù)據(jù)高出2個(gè)數(shù)量級(jí),結(jié)果正確有效.

    圖5是改進(jìn)Randle分?jǐn)?shù)階模型辨識(shí)結(jié)果(表1)的階躍響應(yīng)圖,橫軸是時(shí)間,縱軸是響應(yīng)電壓.其中圖5(a)是激勵(lì)時(shí)長(zhǎng)為0 ~ 100 s的階躍電流在4次參數(shù)辨識(shí)結(jié)果下的階躍響應(yīng).從圖5(a)的擬合效果看,4次辨識(shí)結(jié)果的擬合值基本與預(yù)設(shè)值響應(yīng)結(jié)果重合,擬合程度也“較好”.圖5(b)是對(duì)圖5(a)中從35 ~ 80 s時(shí)間片段(藍(lán)色虛線框部分)的局部放大.從圖5(b)可看出,除了“2#”外,其他3次辨識(shí)結(jié)果擬合曲線與預(yù)設(shè)值“完全重合”,擬合程度較高.

    表1 階躍響應(yīng)辨識(shí)結(jié)果

    圖6是“2#”個(gè)體遺傳進(jìn)化過(guò)程圖,即給出了“2#”遺傳算法優(yōu)化過(guò)程的實(shí)驗(yàn)結(jié)果,橫軸是遺傳算法執(zhí)行代數(shù),縱軸是適應(yīng)度值.由于響應(yīng)電壓的幅度值相對(duì)較小,在每個(gè)點(diǎn)上僅存在10-5的誤差,因此,需要進(jìn)行適應(yīng)度變換,這里使用公式ObiV=log(∑(u-um)2)作為適應(yīng)度函數(shù).由于本身電阻較小,擴(kuò)散系數(shù)很大,導(dǎo)致電壓在10-3量級(jí),因此在圖5(b)的微觀尺度上顯示出很大的差距.從實(shí)際的誤差均值上可看出,擬合后的響應(yīng)電壓與電壓真值之間差距很小.遺傳進(jìn)化過(guò)程圖(圖6)顯示,遺傳到180代時(shí),ObiV已經(jīng)小到-11個(gè)量級(jí),說(shuō)明均方誤差函數(shù)也已達(dá)到10-11,可以認(rèn)為對(duì)模型而言,已經(jīng)達(dá)到較好的精度,也進(jìn)一步驗(yàn)證圖5(b)中,看起來(lái)的差距很大的圖,實(shí)際并不大.

    4 結(jié)論

    本文中主要研究動(dòng)力鋰離子電池的改進(jìn)Randle模型,討論該模型在輸入信號(hào)為階躍電流情況下的解析解,并給出其階躍電壓響應(yīng)曲線,實(shí)現(xiàn)對(duì)該模型在階躍激勵(lì)下模型參數(shù)的提取.在此基礎(chǔ)上,利用遺傳算法對(duì)該系統(tǒng)模型的參數(shù)完成辨識(shí),參數(shù)辨識(shí)誤差均值接近10-5,擬合程度在98%以上.

    猜你喜歡
    階躍適應(yīng)度遺傳算法
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    基于階躍雙包層光纖的螺旋型光纖傳感器
    探討單位階躍信號(hào)的教學(xué)
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類(lèi)分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測(cè)
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    基于改進(jìn)的遺傳算法的模糊聚類(lèi)算法
    少數(shù)民族大學(xué)生文化適應(yīng)度調(diào)查
    一種階躍函數(shù)在矩形時(shí)間窗口頻域特性的分析方法
    亚洲精品第二区| 亚洲精品国产av蜜桃| 人妻 亚洲 视频| 性色av乱码一区二区三区2| 午夜av观看不卡| 久久久精品94久久精品| 免费日韩欧美在线观看| 精品视频人人做人人爽| 久久人妻福利社区极品人妻图片 | 中文字幕高清在线视频| 性色av一级| 日韩视频在线欧美| 国产熟女欧美一区二区| 国产成人av教育| 丰满迷人的少妇在线观看| 亚洲欧美清纯卡通| 亚洲国产毛片av蜜桃av| 成年美女黄网站色视频大全免费| 少妇粗大呻吟视频| 日韩精品免费视频一区二区三区| 永久免费av网站大全| 99热网站在线观看| av不卡在线播放| 伊人久久大香线蕉亚洲五| 免费看不卡的av| 高清黄色对白视频在线免费看| 人人妻人人澡人人看| 国产av国产精品国产| 免费看av在线观看网站| 久久av网站| 人成视频在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| 五月天丁香电影| videos熟女内射| 精品久久久久久久毛片微露脸 | 欧美黄色淫秽网站| 亚洲欧美一区二区三区久久| 国产一区有黄有色的免费视频| 大陆偷拍与自拍| 亚洲国产av新网站| 我要看黄色一级片免费的| 极品少妇高潮喷水抽搐| 美女脱内裤让男人舔精品视频| 国产成人av教育| 真人做人爱边吃奶动态| 91成人精品电影| 亚洲人成77777在线视频| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品国产区一区二| 欧美人与善性xxx| 91国产中文字幕| 久久国产精品人妻蜜桃| 欧美 日韩 精品 国产| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 亚洲中文日韩欧美视频| 久久精品国产a三级三级三级| 久久鲁丝午夜福利片| 免费观看a级毛片全部| 黄片播放在线免费| 老汉色∧v一级毛片| 91精品三级在线观看| 午夜av观看不卡| 久久青草综合色| 妹子高潮喷水视频| 在线天堂中文资源库| 麻豆国产av国片精品| 色视频在线一区二区三区| 国产精品久久久久久精品古装| 免费看不卡的av| 只有这里有精品99| 午夜福利视频精品| 久久精品国产a三级三级三级| 久久国产精品大桥未久av| 国产成人av教育| 99久久99久久久精品蜜桃| 国产深夜福利视频在线观看| 最近最新中文字幕大全免费视频 | 亚洲久久久国产精品| 男人操女人黄网站| 午夜日韩欧美国产| 精品亚洲成a人片在线观看| 国产熟女欧美一区二区| 激情五月婷婷亚洲| 999精品在线视频| 捣出白浆h1v1| 操美女的视频在线观看| 手机成人av网站| 两个人免费观看高清视频| 最新的欧美精品一区二区| 久久午夜综合久久蜜桃| 国产爽快片一区二区三区| 成人亚洲精品一区在线观看| av又黄又爽大尺度在线免费看| 五月天丁香电影| 久久久精品免费免费高清| 婷婷色麻豆天堂久久| 中国国产av一级| 久久久久久免费高清国产稀缺| 人妻一区二区av| 午夜av观看不卡| 91精品国产国语对白视频| 亚洲少妇的诱惑av| 天天躁夜夜躁狠狠久久av| 国产视频首页在线观看| 亚洲精品国产av蜜桃| 黑丝袜美女国产一区| 一本—道久久a久久精品蜜桃钙片| 超碰成人久久| 国产不卡av网站在线观看| 青草久久国产| 欧美日韩国产mv在线观看视频| 国产亚洲精品第一综合不卡| 国产成人精品久久二区二区91| 夫妻午夜视频| 纯流量卡能插随身wifi吗| 国产免费现黄频在线看| 大片电影免费在线观看免费| 黄色 视频免费看| 久久女婷五月综合色啪小说| 啦啦啦在线观看免费高清www| 亚洲欧美精品综合一区二区三区| 九草在线视频观看| 777米奇影视久久| 一级a爱视频在线免费观看| 亚洲,欧美精品.| 亚洲熟女精品中文字幕| 亚洲国产最新在线播放| 另类亚洲欧美激情| 菩萨蛮人人尽说江南好唐韦庄| 最近手机中文字幕大全| h视频一区二区三区| 免费高清在线观看视频在线观看| 亚洲精品av麻豆狂野| videos熟女内射| av网站免费在线观看视频| 国产高清国产精品国产三级| 丝袜在线中文字幕| 大话2 男鬼变身卡| 久久精品久久精品一区二区三区| 欧美精品高潮呻吟av久久| 日韩一本色道免费dvd| 久久人人爽av亚洲精品天堂| 久久久国产欧美日韩av| 91国产中文字幕| 精品少妇内射三级| 国产午夜精品一二区理论片| 国产麻豆69| 老司机亚洲免费影院| 国产精品国产av在线观看| 麻豆国产av国片精品| 丰满少妇做爰视频| 一区二区日韩欧美中文字幕| 少妇人妻 视频| 国产精品久久久久成人av| 久久亚洲精品不卡| 亚洲av成人精品一二三区| 国产精品欧美亚洲77777| 亚洲精品国产区一区二| 国产福利在线免费观看视频| 啦啦啦在线免费观看视频4| 亚洲av在线观看美女高潮| 亚洲七黄色美女视频| 日韩伦理黄色片| 精品少妇久久久久久888优播| 久久精品久久久久久久性| 男人添女人高潮全过程视频| 青草久久国产| 国产老妇伦熟女老妇高清| 亚洲第一青青草原| 免费在线观看黄色视频的| 我要看黄色一级片免费的| 欧美老熟妇乱子伦牲交| 爱豆传媒免费全集在线观看| 一级黄片播放器| 国产三级黄色录像| 亚洲精品国产一区二区精华液| 婷婷色av中文字幕| 精品福利永久在线观看| 久久 成人 亚洲| 国产日韩欧美在线精品| 日本一区二区免费在线视频| 国产精品免费视频内射| 亚洲九九香蕉| 国产精品一区二区免费欧美 | 午夜福利在线免费观看网站| 美女扒开内裤让男人捅视频| 久久精品国产a三级三级三级| 母亲3免费完整高清在线观看| xxxhd国产人妻xxx| 在线观看免费视频网站a站| 欧美人与性动交α欧美精品济南到| 亚洲欧美一区二区三区久久| 一二三四社区在线视频社区8| 成人国产av品久久久| 久久影院123| 国产在线一区二区三区精| av网站免费在线观看视频| 国产av国产精品国产| 精品一区二区三卡| 制服诱惑二区| 一二三四在线观看免费中文在| 女性生殖器流出的白浆| 欧美精品av麻豆av| 精品一区二区三卡| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看| 免费高清在线观看日韩| 精品一区二区三区av网在线观看 | 我要看黄色一级片免费的| 亚洲男人天堂网一区| 极品人妻少妇av视频| 国产精品免费视频内射| 51午夜福利影视在线观看| 久久人人爽av亚洲精品天堂| 国产主播在线观看一区二区 | 国产又爽黄色视频| 亚洲人成电影免费在线| 九色亚洲精品在线播放| 精品国产乱码久久久久久男人| 中文字幕最新亚洲高清| 国产淫语在线视频| 在线av久久热| 嫩草影视91久久| 久久久久久久久久久久大奶| 一本久久精品| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 午夜福利乱码中文字幕| 成年女人毛片免费观看观看9 | 亚洲一区二区三区欧美精品| 日韩免费高清中文字幕av| 人妻一区二区av| 久久久久国产一级毛片高清牌| 亚洲国产欧美日韩在线播放| 精品福利永久在线观看| 天天躁日日躁夜夜躁夜夜| 热99久久久久精品小说推荐| 亚洲综合色网址| 成人国产一区最新在线观看 | 久久久久久久国产电影| 视频在线观看一区二区三区| 中文乱码字字幕精品一区二区三区| 国产一区二区三区综合在线观看| 一边亲一边摸免费视频| 欧美+亚洲+日韩+国产| 国产精品 国内视频| 国产成人精品久久二区二区免费| 蜜桃在线观看..| 女人精品久久久久毛片| 丝瓜视频免费看黄片| 亚洲精品久久久久久婷婷小说| 婷婷色综合www| 一二三四社区在线视频社区8| 国产激情久久老熟女| 丝瓜视频免费看黄片| 嫁个100分男人电影在线观看 | 亚洲人成电影观看| 手机成人av网站| 最近中文字幕2019免费版| 丰满迷人的少妇在线观看| 久久国产精品男人的天堂亚洲| 亚洲人成电影免费在线| 亚洲精品中文字幕在线视频| 十八禁人妻一区二区| 午夜91福利影院| 精品少妇黑人巨大在线播放| 制服人妻中文乱码| 两个人免费观看高清视频| 在线av久久热| 十八禁高潮呻吟视频| 国产精品成人在线| 黄色视频在线播放观看不卡| 国产熟女欧美一区二区| 午夜91福利影院| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 久久亚洲精品不卡| 一边亲一边摸免费视频| 久久国产精品男人的天堂亚洲| 少妇人妻 视频| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 一边摸一边抽搐一进一出视频| 国产男女超爽视频在线观看| 热re99久久国产66热| 在线精品无人区一区二区三| 青春草亚洲视频在线观看| 国产又爽黄色视频| av又黄又爽大尺度在线免费看| www.自偷自拍.com| 日本一区二区免费在线视频| 中文字幕制服av| 欧美日韩成人在线一区二区| 国产黄色视频一区二区在线观看| 伊人亚洲综合成人网| 女人爽到高潮嗷嗷叫在线视频| 精品少妇黑人巨大在线播放| 国产精品一区二区免费欧美 | 汤姆久久久久久久影院中文字幕| 777久久人妻少妇嫩草av网站| 免费女性裸体啪啪无遮挡网站| 热re99久久国产66热| 国产欧美日韩精品亚洲av| 激情五月婷婷亚洲| 久久毛片免费看一区二区三区| 亚洲国产欧美在线一区| 亚洲第一av免费看| 18禁裸乳无遮挡动漫免费视频| 国产欧美日韩一区二区三区在线| 一区在线观看完整版| 亚洲专区国产一区二区| 国产成人精品无人区| 9191精品国产免费久久| 亚洲成色77777| 国产99久久九九免费精品| 又紧又爽又黄一区二区| cao死你这个sao货| 欧美日韩av久久| 中文字幕另类日韩欧美亚洲嫩草| 婷婷色综合www| 精品久久久久久久毛片微露脸 | 超碰成人久久| 精品免费久久久久久久清纯 | 18在线观看网站| a 毛片基地| 国产免费福利视频在线观看| 老司机深夜福利视频在线观看 | 亚洲精品乱久久久久久| 国产成人系列免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看人妻少妇| 在线av久久热| 搡老岳熟女国产| 首页视频小说图片口味搜索 | 一级毛片电影观看| 热99久久久久精品小说推荐| 久久午夜综合久久蜜桃| 免费在线观看完整版高清| av一本久久久久| 男女无遮挡免费网站观看| av又黄又爽大尺度在线免费看| 丁香六月天网| 久久久久国产一级毛片高清牌| 亚洲成人免费av在线播放| 免费不卡黄色视频| 亚洲成人免费av在线播放| 免费不卡黄色视频| av在线老鸭窝| 成年女人毛片免费观看观看9 | 午夜老司机福利片| 女人高潮潮喷娇喘18禁视频| 一本—道久久a久久精品蜜桃钙片| 十八禁高潮呻吟视频| 国产一区二区在线观看av| 亚洲精品国产av蜜桃| 美女主播在线视频| 狠狠精品人妻久久久久久综合| av福利片在线| 亚洲伊人色综图| 极品人妻少妇av视频| 久久精品国产a三级三级三级| 18禁黄网站禁片午夜丰满| 狠狠精品人妻久久久久久综合| 99香蕉大伊视频| 每晚都被弄得嗷嗷叫到高潮| 视频区图区小说| 国产无遮挡羞羞视频在线观看| 青草久久国产| 午夜视频精品福利| 精品一区二区三区四区五区乱码 | 麻豆国产av国片精品| 999精品在线视频| 亚洲天堂av无毛| 久久久精品区二区三区| 一区在线观看完整版| 欧美 日韩 精品 国产| 欧美人与性动交α欧美精品济南到| 欧美亚洲日本最大视频资源| 日本一区二区免费在线视频| 十八禁高潮呻吟视频| e午夜精品久久久久久久| 男女午夜视频在线观看| 大码成人一级视频| 国产av一区二区精品久久| 夫妻午夜视频| 欧美黄色片欧美黄色片| 国精品久久久久久国模美| 国产精品九九99| 日日摸夜夜添夜夜爱| 亚洲中文av在线| 亚洲成人手机| 肉色欧美久久久久久久蜜桃| 色婷婷久久久亚洲欧美| 男女下面插进去视频免费观看| 亚洲国产精品成人久久小说| 又大又爽又粗| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久电影网| 午夜福利视频在线观看免费| 中文精品一卡2卡3卡4更新| 亚洲国产看品久久| 丰满人妻熟妇乱又伦精品不卡| 久久久欧美国产精品| 搡老乐熟女国产| 亚洲男人天堂网一区| 一级毛片女人18水好多 | 亚洲精品日本国产第一区| 人成视频在线观看免费观看| 亚洲熟女毛片儿| 制服诱惑二区| av线在线观看网站| 爱豆传媒免费全集在线观看| 两个人免费观看高清视频| 在线观看国产h片| 欧美黑人欧美精品刺激| 99国产精品99久久久久| 99国产综合亚洲精品| 无限看片的www在线观看| 国产成人精品无人区| 一级a爱视频在线免费观看| 久久99精品国语久久久| 欧美成人精品欧美一级黄| 视频区图区小说| 丁香六月天网| 多毛熟女@视频| 亚洲,欧美精品.| 免费在线观看日本一区| 国产亚洲欧美在线一区二区| 欧美人与性动交α欧美软件| 国产精品 国内视频| 天堂俺去俺来也www色官网| 如日韩欧美国产精品一区二区三区| 晚上一个人看的免费电影| 亚洲欧美成人综合另类久久久| 精品亚洲成a人片在线观看| 婷婷丁香在线五月| a级片在线免费高清观看视频| 亚洲欧美日韩高清在线视频 | 亚洲成av片中文字幕在线观看| 日日爽夜夜爽网站| 性色av乱码一区二区三区2| 欧美日韩av久久| 嫁个100分男人电影在线观看 | 啦啦啦在线观看免费高清www| 50天的宝宝边吃奶边哭怎么回事| 日韩,欧美,国产一区二区三区| 精品少妇内射三级| 在线看a的网站| 久久人人爽人人片av| 香蕉国产在线看| 999久久久国产精品视频| 欧美黄色淫秽网站| 久久久欧美国产精品| 精品少妇久久久久久888优播| 国产成人精品在线电影| 免费日韩欧美在线观看| 一区二区三区精品91| 中文字幕制服av| 一级毛片 在线播放| 男女无遮挡免费网站观看| 青春草亚洲视频在线观看| 国产在线免费精品| 久久精品亚洲熟妇少妇任你| 只有这里有精品99| 国产在线一区二区三区精| 国产男人的电影天堂91| 亚洲,欧美,日韩| 王馨瑶露胸无遮挡在线观看| 大型av网站在线播放| 国产亚洲一区二区精品| 亚洲av成人精品一二三区| 亚洲第一av免费看| 黑人欧美特级aaaaaa片| 久久国产精品大桥未久av| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕| 人人妻,人人澡人人爽秒播 | 国产精品麻豆人妻色哟哟久久| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区综合在线观看| 我要看黄色一级片免费的| 久久精品成人免费网站| 日韩制服丝袜自拍偷拍| 亚洲精品国产一区二区精华液| 欧美人与善性xxx| 亚洲国产最新在线播放| 男女边摸边吃奶| 丝袜在线中文字幕| 精品久久蜜臀av无| 亚洲,欧美,日韩| 精品久久久精品久久久| 又大又爽又粗| av天堂在线播放| 日本猛色少妇xxxxx猛交久久| www.999成人在线观看| 美国免费a级毛片| 亚洲精品日韩在线中文字幕| 国产深夜福利视频在线观看| √禁漫天堂资源中文www| 捣出白浆h1v1| 老司机影院成人| 欧美精品一区二区免费开放| 午夜福利视频在线观看免费| 久久热在线av| 色网站视频免费| 少妇 在线观看| 亚洲欧美成人综合另类久久久| 国产淫语在线视频| 国产免费福利视频在线观看| 精品国产超薄肉色丝袜足j| av片东京热男人的天堂| 国产又色又爽无遮挡免| 大片免费播放器 马上看| 母亲3免费完整高清在线观看| 午夜视频精品福利| 日韩av免费高清视频| 十八禁网站网址无遮挡| 大片免费播放器 马上看| 老司机影院毛片| 欧美乱码精品一区二区三区| 日韩大码丰满熟妇| 99久久人妻综合| 国产成人免费观看mmmm| 国产日韩欧美在线精品| 欧美大码av| 99精国产麻豆久久婷婷| 亚洲熟女毛片儿| 一级毛片 在线播放| 亚洲第一青青草原| 中文字幕精品免费在线观看视频| 中文欧美无线码| 丝瓜视频免费看黄片| svipshipincom国产片| 在线天堂中文资源库| 成年人午夜在线观看视频| 欧美日韩视频精品一区| 欧美97在线视频| 亚洲 欧美一区二区三区| 51午夜福利影视在线观看| 美女视频免费永久观看网站| 国产精品人妻久久久影院| 精品少妇黑人巨大在线播放| 久久久久网色| 国产国语露脸激情在线看| 国产av一区二区精品久久| 亚洲国产欧美网| 久久天躁狠狠躁夜夜2o2o | √禁漫天堂资源中文www| 亚洲av日韩精品久久久久久密 | 亚洲欧美一区二区三区黑人| 飞空精品影院首页| 欧美中文综合在线视频| 性高湖久久久久久久久免费观看| 国产精品免费视频内射| 99久久精品国产亚洲精品| 国产伦理片在线播放av一区| 免费在线观看视频国产中文字幕亚洲 | av网站免费在线观看视频| 18禁观看日本| 亚洲欧洲精品一区二区精品久久久| 大陆偷拍与自拍| 久久精品久久久久久久性| av又黄又爽大尺度在线免费看| 日本黄色日本黄色录像| 一边摸一边抽搐一进一出视频| 日韩伦理黄色片| 亚洲人成网站在线观看播放| 女人爽到高潮嗷嗷叫在线视频| 久久性视频一级片| 日韩视频在线欧美| 美女福利国产在线| 中国美女看黄片| 国产亚洲精品第一综合不卡| 人人妻人人添人人爽欧美一区卜| cao死你这个sao货| 国产在线一区二区三区精| 欧美久久黑人一区二区| 在线看a的网站| 自线自在国产av| 久久影院123| 男女免费视频国产| 亚洲 国产 在线| 欧美精品av麻豆av| 成年美女黄网站色视频大全免费| 亚洲精品国产av蜜桃| 欧美精品av麻豆av| 精品免费久久久久久久清纯 | 亚洲精品国产一区二区精华液| 午夜福利免费观看在线| 欧美成人午夜精品| 可以免费在线观看a视频的电影网站| 亚洲精品美女久久久久99蜜臀 | 99香蕉大伊视频| 国产午夜精品一二区理论片| 又黄又粗又硬又大视频| 观看av在线不卡| 久久毛片免费看一区二区三区| 婷婷色综合www| 在线观看免费午夜福利视频| 亚洲视频免费观看视频| 亚洲欧美清纯卡通| 国产欧美日韩综合在线一区二区| 少妇被粗大的猛进出69影院| 一区二区三区乱码不卡18| 欧美人与性动交α欧美精品济南到| 人妻 亚洲 视频| 我要看黄色一级片免费的| 国产极品粉嫩免费观看在线| 国产精品国产av在线观看| 日本欧美国产在线视频| 亚洲av成人精品一二三区| 欧美黑人欧美精品刺激| 欧美精品亚洲一区二区| 七月丁香在线播放|