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

    基于粗糙集理論的翼型參數(shù)優(yōu)選方法研究

    2016-03-24 08:08:46劉卓李自啟梁濱楊波魏斌斌高永衛(wèi)
    航空工程進(jìn)展 2016年1期
    關(guān)鍵詞:教練機(jī)粗糙集

    劉卓,李自啟,梁濱,楊波,魏斌斌,高永衛(wèi)

    (1.江西洪都航空工業(yè)集團(tuán)有限責(zé)任公司 650所,南昌 330024)

    (2.西北工業(yè)大學(xué) 航空學(xué)院,西安 710072)

    ?

    基于粗糙集理論的翼型參數(shù)優(yōu)選方法研究

    劉卓1,李自啟1,梁濱1,楊波2,魏斌斌2,高永衛(wèi)2

    (1.江西洪都航空工業(yè)集團(tuán)有限責(zé)任公司 650所,南昌330024)

    (2.西北工業(yè)大學(xué) 航空學(xué)院,西安710072)

    摘要:在設(shè)計機(jī)翼時,通常需要進(jìn)行多目標(biāo)決策,而教練機(jī)機(jī)翼翼型選取的先驗經(jīng)驗不足,需要進(jìn)一步研究,從而得出翼型的選擇規(guī)則。分析國內(nèi)外已有教練機(jī)機(jī)翼翼型使用情況,運(yùn)用粗糙集理論中基于重要度和條件信息熵結(jié)合的權(quán)重確定方法,提出在機(jī)翼方案設(shè)計階段翼型的選擇方法;對比分析三種分段方法在翼型性能指標(biāo)離散化時的優(yōu)缺點及其原因,并結(jié)合算例進(jìn)行適用性驗證。結(jié)果表明:在有限數(shù)據(jù)條件下,兩斷點方法是比較合適的屬性離散化方法。所提出的方法具有使用價值,對飛機(jī)設(shè)計方案階段翼型的選擇以及機(jī)翼的氣動設(shè)計具有參考意義。

    關(guān)鍵詞:翼型;教練機(jī);粗糙集

    0引言

    粗糙集(Rough Set)理論是波蘭數(shù)學(xué)家Z.Pawlak[1-2]在20世紀(jì)80年代初首先提出的一種可以分析模糊和不確定問題的數(shù)學(xué)理論,它的基本思想是通過案例庫分類歸納出概念和規(guī)則。經(jīng)過近三十年的發(fā)展,該理論已在機(jī)器學(xué)習(xí)、決策分析、模式識別、人工智能和數(shù)據(jù)挖掘等領(lǐng)域?qū)崿F(xiàn)了成功的應(yīng)用。粗糙集理論最大的優(yōu)點是在進(jìn)行數(shù)據(jù)處理的過程中不需要先驗知識,不受外界主觀因素的影響,僅利用現(xiàn)有數(shù)據(jù)樣本進(jìn)行分析處理,具有高度的客觀性,特別適用于進(jìn)行產(chǎn)品改型設(shè)計卻又缺乏系統(tǒng)理論支持的情況。粗糙集的上述特點相比模糊理論等其他需要依賴專家系統(tǒng)的決策方法具有較大優(yōu)勢。

    多目標(biāo)決策中,在進(jìn)行各目標(biāo)或?qū)傩?備選方案的特征、品質(zhì)或性能參數(shù))的評判和決策過程中通常會使用權(quán)重的概念,用權(quán)重反映各個因素在評判和決策過程中所占有的地位和所起的作用。因而,能否準(zhǔn)確地確定權(quán)重直接影響到最終的評判和決策結(jié)果。

    常用的權(quán)重確定方法有專家評分法、二元對比排序法和模糊統(tǒng)計法等,同時,一些新的研究方法也被運(yùn)用于確定權(quán)重,例如模糊判斷矩陣法、組合賦權(quán)法、誤差分析方法、灰色關(guān)聯(lián)度法和熵權(quán)法等。上述方法中,無論是依據(jù)決策者主觀判斷和經(jīng)驗的主觀賦權(quán)方法,還是基于數(shù)學(xué)理論和方法的客觀賦權(quán)方法,都需要依賴先驗知識。因此,粗糙集理論不需要先驗知識的優(yōu)點可以彌補(bǔ)傳統(tǒng)權(quán)重確定方法的不足,使權(quán)重確定更具客觀性。

    目前,已有不少基于粗糙集理論來求取屬性權(quán)重的研究方法,其中,基于正域重要度的權(quán)重確定方法和基于條件信息熵的權(quán)重確定方法應(yīng)用得較為廣泛[3-5]。然而,基于正域重要度的權(quán)重確定方法,由于其是冗余屬性,且決策表分辨能力不足,導(dǎo)致個別屬性權(quán)重為0;基于條件信息熵的權(quán)重確定方法,對冗余屬性的賦權(quán)過重,導(dǎo)致最終的權(quán)重值中冗余屬性大于非冗余屬性,偏離實質(zhì)。因此,結(jié)合上述兩種權(quán)重確定方法提出了將重要度和條件信息熵相結(jié)合的權(quán)重確定方法,提高了粗糙集方法的適用性[6-9]。

    機(jī)翼的空氣動力特性影響飛機(jī)的性能,翼型的選取和設(shè)計是機(jī)翼設(shè)計的基礎(chǔ)。翼型的外形參數(shù)主要包括:最大厚度和厚度分布、最大彎度和彎度分布、前后緣幾何參數(shù)等。性能參數(shù)主要包括:最大升力系數(shù)CLmax、失速迎角αs、最小阻力系數(shù)CDmin、最大升阻比(L/D)max、升力線斜率CLα和零升力矩系數(shù)CM0等。在設(shè)計飛機(jī)的機(jī)翼時,需要綜合考慮翼型的外形和性能參數(shù),即需要進(jìn)行多目標(biāo)決策。

    教練機(jī)與一般飛機(jī)的典型區(qū)別包括:①教練機(jī)的主要工況是機(jī)動飛行,而非巡航,然而由于缺乏先驗經(jīng)驗,教練機(jī)方案設(shè)計階段仍只能提出巡航設(shè)計要求;②飛機(jī)要有足夠的安全邊界,適合駕駛經(jīng)驗不足的操作者;③教練機(jī)要容易進(jìn)入尾旋,同時也要容易改出等。

    在公開的資料中,適合上述要求的機(jī)翼/翼型性能并無足夠的信息,即關(guān)于教練機(jī)機(jī)翼翼型選取的先驗經(jīng)驗不足,需要進(jìn)一步研究以得出選擇規(guī)則。

    本文運(yùn)用粗糙集理論中基于重要度和條件信息熵結(jié)合的權(quán)重確定方法,提出教練機(jī)機(jī)翼設(shè)計過程中,在機(jī)翼方案設(shè)計階段翼型的選擇方法,對比研究不同分段方法在翼型性能指標(biāo)離散化時的優(yōu)缺點及其原因,并采用風(fēng)洞實驗與CFD方法進(jìn)行算例驗證,以證明本文所用方法的適用性。

    1教練機(jī)翼型參數(shù)的確定

    1.1基于重要度和條件信息熵相結(jié)合的權(quán)重

    想要利用粗糙集方法,首先需建立決策信息表S=(U,A,V,f)。

    表1 決策信息表示例

    基于重要度和條件信息熵相結(jié)合的權(quán)重計算公式為[8]

    (1)

    式中:sig(c)為基于信息表示的條件屬性c的重要度;μ(c)為基于代數(shù)理論的條件屬性c的優(yōu)先級。

    將sig(c)和μ(c)結(jié)合起來定義權(quán)重,一方面考慮到各屬性權(quán)重不應(yīng)為0,另一方面考慮原為“非冗余”屬性指標(biāo)的權(quán)重應(yīng)高于原為“冗余”屬性指標(biāo)的權(quán)重,使得權(quán)重賦值更加合理。

    運(yùn)用基于重要度和條件信息熵結(jié)合的權(quán)重確定方法,得到在相應(yīng)決策變量下各個條件屬性的權(quán)重值,結(jié)果為一個權(quán)重向量:

    W=(w1,w2,…,wm)

    (2)

    1.2效用函數(shù)及其規(guī)則化

    初始決策信息表中所列條件屬性和決策屬性的值一般為具有一定意義的連續(xù)性數(shù)值,若要對其進(jìn)行粗糙集處理,需要將其賦予數(shù)學(xué)上的意義。第一步是對各個屬性規(guī)則化處理,即應(yīng)用因素效用函數(shù)方法將每列屬性值化為介于[0,1]的無量綱量,然后便可以使用統(tǒng)一的測度及加權(quán)來進(jìn)行定量評價與比較。

    若待評價的屬性為m個,且每個屬性對應(yīng)的樣本值有n個,則信息表可理解為一個m×n的屬性矩陣:

    C=(cij)m×n

    (3)

    (4)

    (5)

    則有以下三種數(shù)值規(guī)則化的情況:

    (1) 若屬性cij要求越大越好,則采用的效用函數(shù)為

    (6)

    (2) 若屬性cij要求越小越好,則采用的效用函數(shù)為

    (7)

    (3) 若屬性cij要求應(yīng)保持在[r1,r2]范圍內(nèi)時,則采用的效用函數(shù)為

    (i=1,2,…,m)

    (8)

    第三種情況存在特例,即當(dāng)屬性cij要求盡量接近某一定值r時,可將式(8)作特例處理,得到效用函數(shù):

    (i=1,2,…,m)

    (9)

    從屬性數(shù)值規(guī)則化的過程可以看出:決策對于各個屬性的偏好信息已經(jīng)體現(xiàn)在效用函數(shù)的規(guī)則化過程中,例如某些屬性需要越大越好,某些屬性需要趨近一個定值等。將偏好信息用于數(shù)值規(guī)則化進(jìn)而完成粗糙集的權(quán)重分析,便可達(dá)到?jīng)Q策所要求的效用,即因素效用函數(shù)的內(nèi)涵所在。完成各屬性的規(guī)則化后,數(shù)據(jù)仍然是連續(xù)性的數(shù)值,需要進(jìn)一步做離散化處理。

    在分析得到關(guān)于方案各個條件屬性(參數(shù))的權(quán)重后,將方案各條件屬性(參數(shù))按照效用函數(shù)規(guī)則化,得到效用值矩陣

    (10)

    各個備選方案的相對量化權(quán)重值的計算公式為

    (11)

    比較各個方案的最終綜合效用值,該值越大者即表明其所對應(yīng)的方案越好。

    1.3屬性的離散化

    要進(jìn)行數(shù)字化處理,則必須對連續(xù)變量進(jìn)行相應(yīng)的離散化處理。

    等距離散化方法是將每個屬性無量綱化后的值區(qū)間用一個或多個斷點等距分成多段,并將各段定義為不同等級,再用數(shù)字編碼組合表示,每個屬性的無量綱化值落在某分段內(nèi)則編碼為分段對應(yīng)的數(shù)字。

    在有限樣本數(shù)量的情況下,本文重點研究一斷點、兩斷點和三斷點等距離散化方法的優(yōu)缺點,并提出一種適用于翼型選型的最佳方法。

    一斷點、兩斷點和三斷點等距離散化方法對應(yīng)的離散等級分別為兩個等級(好、差)、三個等級(好、一般和差)和四個等級(好、較好、一般和差),分別用數(shù)字編碼組合表示為(2、1)、(3、2、1)以及(4、3、2、1)。為了得到教練機(jī)機(jī)翼翼型性能的選擇規(guī)則,搜集了20架相關(guān)教練機(jī)機(jī)翼翼根和翼梢的翼型作為初始數(shù)據(jù)樣本,忽略先驗經(jīng)驗,基于粗糙集理論進(jìn)行歸納分析。根據(jù)設(shè)計狀態(tài)為巡航狀態(tài)的要求,確定翼型優(yōu)選的決策變量為翼載和巡航馬赫數(shù)。列出20種教練機(jī)型號相應(yīng)的決策變量值,統(tǒng)計結(jié)果如表2所示。

    將關(guān)于翼型的條件屬性(參數(shù))與決策變量結(jié)合可得決策信息表,分別建立關(guān)于翼根翼型和翼梢翼型的決策信息表,兩者決策變量和對應(yīng)數(shù)值一致。

    運(yùn)用效用函數(shù)進(jìn)行條件屬性無量綱化處理,翼型各項性能參數(shù)的效用分別為:最大升力系數(shù)(CLmax)、失速迎角(αs)及最大升阻比((L/D)max)要求越大越好,最小阻力系數(shù)(CDmin)、升力線斜率(CLα)及零升力矩的絕對值要求越小越好。兩個決策變量中,翼載和巡航馬赫數(shù)均要求接近定值,分別為165.65kgf/m2和0.36(針對PC-7飛機(jī))。

    首先,按照上述要求進(jìn)行決策信息表各屬性(參數(shù))基于效用函數(shù)的無量綱化處理,然后,根據(jù)三種離散化方法分別建立離散決策信息表,其中兩個決策變量按照各0.5的權(quán)重線性加權(quán)為一個綜合評價指標(biāo)。所列翼型性能參閱文獻(xiàn)[10]。

    三種不同離散化方法下翼根翼型的離散決策信息表分別如表3~表5所示,限于篇幅,本文僅列出翼根翼型的相關(guān)數(shù)據(jù)。

    表5 教練機(jī)機(jī)翼翼根翼型一斷點離散決策信息表

    利用基于重要度和條件信息熵相結(jié)合的條件屬性權(quán)重確定方法,結(jié)合決策信息表(表3~表5),分別得到三種離散化方法下翼根翼型和翼梢翼型在決策變量約束下各條件屬性(氣動參數(shù))的權(quán)重,如表6~表7所示。

    表6 教練機(jī)機(jī)翼翼根翼型各氣動參數(shù)的權(quán)重值

    表7 教練機(jī)機(jī)翼翼梢翼型各氣動參數(shù)的權(quán)重值

    從表6~表7可以看出:三種不同離散化方法求得的條件屬性(氣動參數(shù))的權(quán)重存在差別,三斷點與兩斷點離散化方法的結(jié)果比較接近,一斷點離散化方法的結(jié)果偏差較大,這是因為各條件屬性的值域范圍較大,僅采用一個斷點劃分容易導(dǎo)致數(shù)據(jù)失真,離散化后的數(shù)據(jù)表不能真實反映各翼型對應(yīng)條件屬性的取值狀態(tài);兩斷點與三斷點離散化方法結(jié)果大致相同,在翼根翼型的分析中兩者最大升阻比和升力線斜率對應(yīng)的權(quán)重不同,其中決策變量對應(yīng)的狀態(tài)是巡航狀態(tài),最大升阻比所占權(quán)重應(yīng)當(dāng)較大,而三斷點離散化方法的結(jié)果顯示最大升阻比的權(quán)重較小,兩斷點離散化方法的結(jié)果更接近合理結(jié)果。因此,本文采用兩斷點等距離散化方法作為連續(xù)屬性離散化的方法。

    由兩斷點等距離散化方法得到的翼根翼型和翼梢翼型性能參數(shù)的權(quán)重可知:在決策變量的約束下,翼根翼型最小阻力系數(shù)、最大升阻比和零升力矩系數(shù)占比重比較大,故翼根翼型的選擇更看重上述三個參數(shù);翼梢翼型最小阻力系數(shù)和升力線斜率占比重比較大,故翼梢翼型選擇時更看重該兩個參數(shù)。

    2算例驗證

    為了進(jìn)一步驗證教練機(jī)機(jī)翼翼型優(yōu)選方法的可行性,根據(jù)PC-7飛機(jī)性能和上述翼型性能參數(shù)選擇規(guī)則,參照文獻(xiàn)[11-12]中的方法,設(shè)計兩組翼根、翼梢翼型并進(jìn)行機(jī)翼設(shè)計對比。限于篇幅,僅列出翼根翼型的設(shè)計結(jié)果,其氣動性能參數(shù)如表8所示。

    將表8按照效用函數(shù)方法做無量綱化處理,得到無量綱化后的效用值數(shù)據(jù)表,并計算各個備選翼型方案的相對量化權(quán)重值,如表9所示。

    表8 翼根備選翼型的氣動參數(shù)表

    表9 翼根備選翼型的氣動參數(shù)無量綱化效用值及相對量化權(quán)重

    從最終的計算結(jié)果可知:備選翼型中,翼根翼型以TR18-1的相對量化權(quán)重值高,翼梢翼型以TR12-1的相對量化權(quán)重值高。表明這一對翼型在翼根翼型和翼梢翼型決策中具有更多趨向于決策的優(yōu)勢。因此,新設(shè)計的機(jī)翼采用的翼根翼型為TR18-1,翼梢翼型為TR12-1。

    經(jīng)CFD和風(fēng)洞實驗驗證,新機(jī)翼比基準(zhǔn)機(jī)翼升阻比高2、最小阻力小約5%,表明根據(jù)上述準(zhǔn)則選擇的翼型取得了滿意的效果。

    3結(jié)論

    (1) 運(yùn)用粗糙集理論,基于重要度和條件信息熵結(jié)合的權(quán)重確定方法對翼型進(jìn)行優(yōu)選的方法是可行的。

    (2) 在教練機(jī)翼型的權(quán)重確定過程中,采用兩斷點法是合適的。

    (3) 在本文的研究范圍內(nèi),翼根翼型看重最小阻力系數(shù)、最大升阻比和零升力矩系數(shù),而翼梢翼型主要參考最小阻力系數(shù)和升力線斜率。

    參考文獻(xiàn)

    [1]PawlakZ.Roughsets[J].InternationalJournalofComputerandInformationSciences, 1982, 11(5): 341-356.

    [2]PawlakZ.Roughsettheoryanditsapplicationstodataanalysis[J].CyberneticsandSystems, 1998, 29(7): 661-688.

    [3] 譚宗鳳. 基于粗糙集的權(quán)重確定方法研究[D]. 桂林: 廣西師范大學(xué), 2012.

    TanZongfeng.Researchonmethodofattributeweightbasedonroughsetstheory[D].Guilin:GuangxiNormalUniversity, 2012.(inChinese)

    [4] 鮑新中, 張建斌, 劉澄. 基于粗糙集條件信息熵的權(quán)重確定方法[J]. 中國管理科學(xué), 2009, 17(3): 131-135.

    BaoXinzhong,ZhangJianbin,LiuCheng.Anewmethodofascertainingattributeweightbasedonroughsetsconditionalinformationentropy[J].ChineseJournalofManagementScience, 2009, 17(3): 131-135.(inChinese)

    [5] 徐澤水. 幾類多屬性決策方法研究[D]. 南京: 東南大學(xué), 2002.

    XuZeshui.Studyonmethodsformultipleattributedecisionmakingundersomesituations[D].Nanjing:SoutheastUniversity, 2002.(inChinese)

    [6] 曹秀英, 梁靜國. 基于粗集理論的屬性權(quán)重確定方法[J]. 中國管理科學(xué), 2002, 10(5): 98-100.

    CaoXiuying,LiangJingguo.Themethodofascertainingattributeweightbasedonroughsetstheory[J].ChineseJournalofManagementScience, 2002, 10(5): 98-100.(inChinese)

    [7] 陳鳳娟, 閆德勤. 知識粗糙性和條件信息熵的關(guān)系[J]. 計算機(jī)工程與應(yīng)用, 2009, 45(21): 160-162.

    ChenFengjuan,YanDeqin.Relationshipbetweenroughnessofknowledgeandconditionalinformationentropy[J].ComputerEngineeringandApplications, 2009, 45(21): 160-162.(inChinese)

    [8] 朱紅燦, 陳能華. 粗糙集條件信息熵權(quán)重確定方法的改進(jìn)[J]. 統(tǒng)計與決策, 2011(8): 154-156.

    ZhuHongcan,ChenNenghua.Theimprovementofroughsetconditionalentropyweightdeterminingmethod[J].Statisticsanddecision, 2011(8): 154-156.(inChinese)

    [9]ZiarkoW.Introductiontothespecialissueonroughsetsandknowledgediscovery[J].InternationalofComputationalIntelligence, 1995, 11(2): 223-226.

    [10] 飛機(jī)設(shè)計手冊總編委會. 飛機(jī)設(shè)計手冊: 第6冊——氣動設(shè)計[M]. 北京: 航空工業(yè)出版社, 2002.

    GeneralEditorialBoardofAircraft Design Manual.Aircraftdesignmanual:Vol6-Aerodynamicdesign[M].Beijing:AviationIndustryPress, 2002.(inChinese)

    [11]AbbottIH,AlbertE,VonDoenhoff.Theoryofwingsec-

    tions:including a summary of airfoil data[M]. New York: Dover Publications,1959.

    [12] Patankar S V, Spalding D B. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows[J]. International Journal of Heat and Mass Transfer, 1972, 15: 1787-1806.

    劉卓(1981-),男,高級工程師。主要研究方向:飛機(jī)總體設(shè)計。

    李自啟(1985-),男,碩士,工程師。主要研究方向:飛機(jī)總體設(shè)計。

    梁斌(1982-),男,工程師。主要研究方向:飛機(jī)總體設(shè)計。

    楊波(1989-),男,碩士研究生。主要研究方向:空氣動力學(xué)。

    魏斌斌(1990-),男,碩士研究生。主要研究方向:空氣動力學(xué)。

    高永衛(wèi)(1968-),男,博士,教授。主要研究方向:空氣動力學(xué)。

    (編輯:馬文靜)

    Method of Airfoil Selection in Wing Design Using Rough Set Theory

    Liu Zhuo1, Li Ziqi1, Liang Bin1, Yang Bo2, Wei Binbin2, Gao Yongwei2

    (1.650 Research Institute, Jiangxi Hongdu Aviation Industry GroupCorporation Limited, Nanchang 330024, China)(2.School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)

    Abstract:In a new design of an aircraft wing, the multiple objective decision making is needed. But in selecting training aircraft wings’ airfoils the prior experience is insufficient, so it is needed to further study and conclude the selection rules of the wing airfoil. The airfoils used in training aircraft wings are collected and a method of airfoil selecting for wing design by using rough set theory is provided. Weight reflecting importance of each factor is the key parameter in evaluating or decision-making. The method of weight determining combining the attribute importance with conditional information entropy, which is conformed by conducting a project of trainer wing design, is suitable for airfoil selection in wing design. The number of breaking point in discretization of parameter of airfoil is studied. It is advised that two breaking point method is preferable than one and three, when the number of airfoils in data base available is limited. The method used in this paper is believed to be useful for trainer aircraft wing’s design.

    Key words:airfoil; trainer aircraft; rough set theory

    作者簡介:

    中圖分類號:V271.6

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

    DOI:10.16615/j.cnki.1674-8190.2016.01.012

    文章編號:1674-8190(2016)01-087-07

    通信作者:魏斌斌,1274724740@qq.com

    收稿日期:2015-12-17;修回日期:2015-12-29

    猜你喜歡
    教練機(jī)粗糙集
    英國Aerospace Hawk 60系列教練機(jī)(紅箭飛行表演隊涂裝)
    第十頁 知識點 教-8教練機(jī)
    翱翔全球的K8教練機(jī)
    航空世界(2021年9期)2021-09-29 12:28:51
    基于Pawlak粗糙集模型的集合運(yùn)算關(guān)系
    基于二進(jìn)制鏈表的粗糙集屬性約簡
    優(yōu)勢直覺模糊粗糙集決策方法及其應(yīng)用
    多?;植诩再|(zhì)的幾個充分條件
    雙論域粗糙集在故障診斷中的應(yīng)用
    兩個域上的覆蓋變精度粗糙集模型
    《教練機(jī)》發(fā)刊詞
    教練機(jī)(2010年4期)2010-08-15 00:53:33
    国产精品熟女久久久久浪| 国产一区二区三区av在线| 精品久久久久久电影网| 一个人观看的视频www高清免费观看| 国产老妇伦熟女老妇高清| 国产69精品久久久久777片| 久久久a久久爽久久v久久| 看免费成人av毛片| 五月伊人婷婷丁香| 国产黄色小视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲在线自拍视频| 天堂中文最新版在线下载 | 国产爱豆传媒在线观看| 亚洲精品,欧美精品| 九九在线视频观看精品| 黄色一级大片看看| 亚洲最大成人中文| 看免费成人av毛片| 亚洲经典国产精华液单| 久久久久国产网址| 精品国产露脸久久av麻豆 | 午夜福利视频1000在线观看| 伊人久久国产一区二区| 女人久久www免费人成看片| 黄片无遮挡物在线观看| 高清av免费在线| 亚洲精品国产av成人精品| 18禁动态无遮挡网站| 赤兔流量卡办理| 国产一级毛片七仙女欲春2| 又大又黄又爽视频免费| 国产精品久久久久久av不卡| 久久精品国产亚洲av天美| 久久久久久久久久成人| 国产三级在线视频| 嫩草影院新地址| 久久人人爽人人片av| 国产淫语在线视频| 天天躁夜夜躁狠狠久久av| 国产欧美另类精品又又久久亚洲欧美| 男人和女人高潮做爰伦理| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美成人精品一区二区| 国产一区二区三区综合在线观看 | 大香蕉久久网| 免费大片黄手机在线观看| av天堂中文字幕网| 国内少妇人妻偷人精品xxx网站| av免费在线看不卡| 久久久精品欧美日韩精品| 亚洲自拍偷在线| 成人美女网站在线观看视频| 80岁老熟妇乱子伦牲交| 人妻制服诱惑在线中文字幕| 美女脱内裤让男人舔精品视频| 在线天堂最新版资源| 在线天堂最新版资源| 男女边摸边吃奶| 在线观看一区二区三区| 777米奇影视久久| 嫩草影院精品99| 免费观看无遮挡的男女| 国产亚洲av嫩草精品影院| 亚洲精品第二区| 狠狠精品人妻久久久久久综合| 十八禁国产超污无遮挡网站| 80岁老熟妇乱子伦牲交| 国产成人freesex在线| 91久久精品电影网| 中文字幕av成人在线电影| 国产精品无大码| 成人漫画全彩无遮挡| 人妻系列 视频| 天堂av国产一区二区熟女人妻| 成人毛片60女人毛片免费| 免费不卡的大黄色大毛片视频在线观看 | 日本wwww免费看| 哪个播放器可以免费观看大片| 一级av片app| 99久国产av精品国产电影| 国产精品久久久久久精品电影小说 | 国产精品美女特级片免费视频播放器| 最近2019中文字幕mv第一页| 国产91av在线免费观看| 久久精品国产亚洲av涩爱| 看免费成人av毛片| 26uuu在线亚洲综合色| 色综合亚洲欧美另类图片| 国产精品日韩av在线免费观看| 一级毛片 在线播放| 亚洲性久久影院| 亚洲国产日韩欧美精品在线观看| 伦理电影大哥的女人| 国产精品一区二区在线观看99 | 三级经典国产精品| 肉色欧美久久久久久久蜜桃 | 精华霜和精华液先用哪个| 国产精品一区二区三区四区久久| 18禁在线播放成人免费| 日韩欧美三级三区| 不卡视频在线观看欧美| 国产一区二区三区综合在线观看 | 久久精品熟女亚洲av麻豆精品 | 精华霜和精华液先用哪个| 亚洲精品乱久久久久久| 中文乱码字字幕精品一区二区三区 | 国产精品福利在线免费观看| 国产免费视频播放在线视频 | 精品午夜福利在线看| 99久国产av精品国产电影| 亚洲在线观看片| 我的老师免费观看完整版| 国产一区亚洲一区在线观看| 国产在视频线精品| 亚洲欧美清纯卡通| 蜜臀久久99精品久久宅男| 亚洲精品亚洲一区二区| 国产成年人精品一区二区| 女人十人毛片免费观看3o分钟| 国产精品国产三级国产专区5o| 乱码一卡2卡4卡精品| 国内精品一区二区在线观看| 亚洲精品亚洲一区二区| 欧美精品国产亚洲| 色播亚洲综合网| 午夜爱爱视频在线播放| 亚洲国产成人一精品久久久| 国产淫语在线视频| 夫妻性生交免费视频一级片| 日本三级黄在线观看| 熟妇人妻久久中文字幕3abv| 国产黄片美女视频| 免费看a级黄色片| 美女黄网站色视频| 成人漫画全彩无遮挡| 韩国av在线不卡| 亚洲三级黄色毛片| 中文资源天堂在线| 精品一区二区三区人妻视频| 免费看a级黄色片| 我要看日韩黄色一级片| 一本一本综合久久| 又爽又黄a免费视频| 久久鲁丝午夜福利片| 亚洲av中文av极速乱| 精品久久国产蜜桃| 亚洲av免费高清在线观看| 欧美 日韩 精品 国产| 国产午夜精品一二区理论片| 成人亚洲欧美一区二区av| 亚洲精品视频女| 国产精品精品国产色婷婷| 女人被狂操c到高潮| 免费无遮挡裸体视频| 久久99热这里只有精品18| 青青草视频在线视频观看| 搡老妇女老女人老熟妇| 国产乱人视频| 2022亚洲国产成人精品| 男人舔女人下体高潮全视频| 久久久久久久午夜电影| 大香蕉97超碰在线| 日本三级黄在线观看| 免费看光身美女| 亚洲乱码一区二区免费版| 欧美激情在线99| 精品99又大又爽又粗少妇毛片| 国产熟女欧美一区二区| 亚洲av免费在线观看| 精品久久久久久电影网| 免费av毛片视频| 国产高清国产精品国产三级 | av在线亚洲专区| av.在线天堂| av免费观看日本| 国产黄a三级三级三级人| 毛片一级片免费看久久久久| 一本久久精品| 激情五月婷婷亚洲| 久久精品国产鲁丝片午夜精品| 国产亚洲精品久久久com| 综合色丁香网| 我要看日韩黄色一级片| 麻豆成人午夜福利视频| 天堂影院成人在线观看| 看免费成人av毛片| 日本与韩国留学比较| 久久韩国三级中文字幕| 成人综合一区亚洲| 久久久久精品久久久久真实原创| 亚洲欧美一区二区三区国产| 大陆偷拍与自拍| 国产成人福利小说| 亚洲美女搞黄在线观看| 纵有疾风起免费观看全集完整版 | 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 免费播放大片免费观看视频在线观看| 三级毛片av免费| 少妇熟女aⅴ在线视频| 高清日韩中文字幕在线| 国产单亲对白刺激| 国产精品综合久久久久久久免费| 性插视频无遮挡在线免费观看| 精品人妻一区二区三区麻豆| 国产av在哪里看| 男女视频在线观看网站免费| 国产午夜精品一二区理论片| 免费电影在线观看免费观看| 久久草成人影院| 日本-黄色视频高清免费观看| 成人漫画全彩无遮挡| 尤物成人国产欧美一区二区三区| 亚洲成人精品中文字幕电影| 高清欧美精品videossex| 中文资源天堂在线| 97超视频在线观看视频| 中文字幕制服av| 1000部很黄的大片| 欧美日韩国产mv在线观看视频 | 国产综合懂色| 免费观看a级毛片全部| 婷婷色综合www| 亚洲无线观看免费| 亚洲内射少妇av| 精品99又大又爽又粗少妇毛片| 天堂中文最新版在线下载 | 高清av免费在线| 欧美日韩国产mv在线观看视频 | 国产视频首页在线观看| 99久久精品国产国产毛片| 国产v大片淫在线免费观看| 午夜免费观看性视频| av.在线天堂| 亚洲国产精品专区欧美| av国产久精品久网站免费入址| 菩萨蛮人人尽说江南好唐韦庄| 舔av片在线| 春色校园在线视频观看| av黄色大香蕉| 大话2 男鬼变身卡| 亚洲自偷自拍三级| 亚洲va在线va天堂va国产| 在线观看人妻少妇| 日韩欧美 国产精品| 国产精品日韩av在线免费观看| 亚洲国产色片| 成人欧美大片| 亚洲欧美精品专区久久| 欧美zozozo另类| av线在线观看网站| 我的老师免费观看完整版| 最近视频中文字幕2019在线8| 成年av动漫网址| av天堂中文字幕网| 蜜臀久久99精品久久宅男| 亚洲自偷自拍三级| 能在线免费观看的黄片| 国产在线男女| 三级国产精品欧美在线观看| 婷婷六月久久综合丁香| 午夜免费激情av| 国产色爽女视频免费观看| 天堂中文最新版在线下载 | 欧美潮喷喷水| 国产永久视频网站| 国产一区二区亚洲精品在线观看| 成人鲁丝片一二三区免费| 夫妻性生交免费视频一级片| 色综合亚洲欧美另类图片| 亚洲国产精品成人久久小说| 国产激情偷乱视频一区二区| 亚洲性久久影院| 亚洲成人精品中文字幕电影| 欧美成人a在线观看| 尤物成人国产欧美一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲av男天堂| 三级毛片av免费| 国产伦一二天堂av在线观看| 日韩一本色道免费dvd| 亚洲精品日韩av片在线观看| 91狼人影院| 亚洲欧美一区二区三区黑人 | 建设人人有责人人尽责人人享有的 | 成人鲁丝片一二三区免费| 高清视频免费观看一区二区 | 国产单亲对白刺激| av福利片在线观看| 黑人高潮一二区| 三级国产精品片| 亚洲av日韩在线播放| 嫩草影院入口| 国产伦理片在线播放av一区| 日韩强制内射视频| 成人午夜精彩视频在线观看| 亚洲国产欧美在线一区| 乱码一卡2卡4卡精品| 精品久久久久久久久亚洲| 国产精品美女特级片免费视频播放器| av在线蜜桃| 欧美日韩国产mv在线观看视频 | 成人欧美大片| 午夜激情久久久久久久| 久久久久久久亚洲中文字幕| 色综合亚洲欧美另类图片| 91精品国产九色| 国产黄色免费在线视频| 三级毛片av免费| 国产成人精品久久久久久| 亚洲婷婷狠狠爱综合网| 岛国毛片在线播放| 最近中文字幕高清免费大全6| 精品一区二区三卡| 日韩中字成人| 婷婷色综合www| 国产日韩欧美在线精品| 1000部很黄的大片| 久久精品综合一区二区三区| 在现免费观看毛片| 亚洲丝袜综合中文字幕| av播播在线观看一区| av卡一久久| 欧美潮喷喷水| 精品人妻视频免费看| 老司机影院成人| 欧美最新免费一区二区三区| 在线观看免费高清a一片| 国产成人一区二区在线| 听说在线观看完整版免费高清| 日韩一区二区视频免费看| 亚洲精品色激情综合| av线在线观看网站| 色综合亚洲欧美另类图片| 亚洲av电影不卡..在线观看| 日韩伦理黄色片| 超碰97精品在线观看| 日日啪夜夜撸| 99热这里只有是精品在线观看| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 老司机影院毛片| 卡戴珊不雅视频在线播放| 内射极品少妇av片p| 精品一区在线观看国产| 国产精品人妻久久久久久| 男女边吃奶边做爰视频| 国产成人精品一,二区| 高清日韩中文字幕在线| 免费黄频网站在线观看国产| 日韩制服骚丝袜av| 中文字幕人妻熟人妻熟丝袜美| 精品国产三级普通话版| 蜜桃久久精品国产亚洲av| av国产免费在线观看| 国产亚洲一区二区精品| 2021天堂中文幕一二区在线观| 久久人人爽人人片av| 亚洲av日韩在线播放| 日日摸夜夜添夜夜爱| 国产成人a区在线观看| 国产免费福利视频在线观看| 97在线视频观看| 亚洲一区高清亚洲精品| 国产亚洲精品av在线| 国产不卡一卡二| 熟女人妻精品中文字幕| 天天一区二区日本电影三级| 亚洲精品色激情综合| 蜜桃亚洲精品一区二区三区| 国产一级毛片七仙女欲春2| 能在线免费看毛片的网站| av天堂中文字幕网| 亚洲欧美清纯卡通| 精品一区二区三区人妻视频| 能在线免费观看的黄片| 婷婷色麻豆天堂久久| 午夜激情欧美在线| 日日啪夜夜撸| 亚洲精品中文字幕在线视频 | 亚洲精品影视一区二区三区av| 黄色配什么色好看| 日韩av在线大香蕉| 国产高清不卡午夜福利| 男插女下体视频免费在线播放| 色网站视频免费| 丰满少妇做爰视频| 欧美xxⅹ黑人| 国产精品一区www在线观看| 国产女主播在线喷水免费视频网站 | 不卡视频在线观看欧美| 91久久精品国产一区二区成人| 天天一区二区日本电影三级| 毛片一级片免费看久久久久| 国产精品.久久久| 国产成人a∨麻豆精品| 成人亚洲精品一区在线观看 | 最新中文字幕久久久久| 国产在线一区二区三区精| 欧美bdsm另类| 午夜激情欧美在线| 欧美不卡视频在线免费观看| 国内揄拍国产精品人妻在线| 亚洲欧美清纯卡通| 免费av观看视频| 26uuu在线亚洲综合色| 国产美女午夜福利| 一个人看的www免费观看视频| 国产视频首页在线观看| 观看免费一级毛片| eeuss影院久久| 国产麻豆成人av免费视频| 国产在视频线在精品| 婷婷色综合www| 日韩在线高清观看一区二区三区| 国产精品一及| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕 | 中文在线观看免费www的网站| 免费观看a级毛片全部| 午夜激情欧美在线| 插逼视频在线观看| 69av精品久久久久久| 中文字幕av在线有码专区| 久久久久久久久大av| 中文字幕制服av| 日日啪夜夜撸| 亚洲av免费高清在线观看| 大香蕉久久网| 69人妻影院| 中国美白少妇内射xxxbb| 国产成人精品婷婷| 久久久久久久大尺度免费视频| 国产精品综合久久久久久久免费| 18禁裸乳无遮挡免费网站照片| 国产黄a三级三级三级人| 精品欧美国产一区二区三| 久久久久久久久久成人| 久久久国产一区二区| 色尼玛亚洲综合影院| 中文天堂在线官网| 午夜老司机福利剧场| 成人高潮视频无遮挡免费网站| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 国产精品国产三级专区第一集| 日本免费a在线| 国产激情偷乱视频一区二区| 黄片wwwwww| 免费在线观看成人毛片| 亚洲无线观看免费| 亚洲精品成人久久久久久| 在现免费观看毛片| 国产日韩欧美在线精品| 国产黄频视频在线观看| 日韩亚洲欧美综合| 最近中文字幕2019免费版| 亚洲精品国产av蜜桃| 日韩,欧美,国产一区二区三区| 国产国拍精品亚洲av在线观看| 亚洲av电影不卡..在线观看| 汤姆久久久久久久影院中文字幕 | 亚洲精品乱久久久久久| 特大巨黑吊av在线直播| 寂寞人妻少妇视频99o| 少妇猛男粗大的猛烈进出视频 | 日日摸夜夜添夜夜添av毛片| 国产一区亚洲一区在线观看| 成年av动漫网址| 国产三级在线视频| 最近的中文字幕免费完整| av网站免费在线观看视频 | 久久久精品免费免费高清| 国产成人福利小说| 欧美97在线视频| 在线免费观看的www视频| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 一区二区三区乱码不卡18| 大陆偷拍与自拍| 一二三四中文在线观看免费高清| 亚洲性久久影院| 噜噜噜噜噜久久久久久91| 一级毛片aaaaaa免费看小| 在线观看一区二区三区| 一级黄片播放器| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 久久久久久久久久成人| 激情 狠狠 欧美| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 国产 一区 欧美 日韩| 免费电影在线观看免费观看| 亚洲精品一区蜜桃| 精品久久久久久电影网| 成人欧美大片| 国产高清有码在线观看视频| 午夜福利在线观看免费完整高清在| 三级毛片av免费| 久久精品国产自在天天线| 久久久久久久亚洲中文字幕| 99九九线精品视频在线观看视频| 亚洲经典国产精华液单| 国产精品一区www在线观看| 美女大奶头视频| 精品久久久久久久久久久久久| 九色成人免费人妻av| 熟妇人妻久久中文字幕3abv| 亚洲综合精品二区| 免费高清在线观看视频在线观看| 日本三级黄在线观看| 高清av免费在线| 久久精品国产亚洲av涩爱| 亚洲va在线va天堂va国产| 午夜日本视频在线| 亚洲图色成人| 久久久精品94久久精品| 女人久久www免费人成看片| 好男人在线观看高清免费视频| 街头女战士在线观看网站| 精品人妻一区二区三区麻豆| 国产亚洲最大av| 女人久久www免费人成看片| 亚洲精品日本国产第一区| 国产精品无大码| 国产综合懂色| 一级二级三级毛片免费看| 在线a可以看的网站| 最近的中文字幕免费完整| 最新中文字幕久久久久| 欧美潮喷喷水| 亚洲国产精品成人综合色| 国产女主播在线喷水免费视频网站 | 伊人久久精品亚洲午夜| 一本一本综合久久| 亚洲成人精品中文字幕电影| 国产高清三级在线| 欧美性感艳星| av.在线天堂| 久久99热6这里只有精品| 久久久久久久久大av| 免费av观看视频| 欧美成人一区二区免费高清观看| 国产乱人视频| 99热这里只有精品一区| 国产国拍精品亚洲av在线观看| 国产一区二区三区综合在线观看 | 青青草视频在线视频观看| 一级毛片黄色毛片免费观看视频| 麻豆av噜噜一区二区三区| 久久久久久伊人网av| 小蜜桃在线观看免费完整版高清| 亚洲婷婷狠狠爱综合网| 麻豆久久精品国产亚洲av| 乱系列少妇在线播放| 欧美三级亚洲精品| 成人毛片a级毛片在线播放| 亚洲最大成人手机在线| 日本免费在线观看一区| 亚洲四区av| 国产爱豆传媒在线观看| 国产精品熟女久久久久浪| 1000部很黄的大片| 成人av在线播放网站| 日韩av在线免费看完整版不卡| 久久人人爽人人片av| 欧美日韩在线观看h| 中文字幕亚洲精品专区| 国产精品久久久久久精品电影| 精品人妻一区二区三区麻豆| 精品午夜福利在线看| 美女被艹到高潮喷水动态| 国产亚洲最大av| 亚洲在久久综合| 在线观看免费高清a一片| 插逼视频在线观看| 中文乱码字字幕精品一区二区三区 | 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 精华霜和精华液先用哪个| 观看免费一级毛片| 国产高清国产精品国产三级 | 国语对白做爰xxxⅹ性视频网站| 两个人的视频大全免费| 欧美丝袜亚洲另类| 婷婷色综合大香蕉| 成人av在线播放网站| 亚洲精品视频女| 国语对白做爰xxxⅹ性视频网站| 国产精品综合久久久久久久免费| 亚洲精品中文字幕在线视频 | 夜夜爽夜夜爽视频| 一级毛片电影观看| 一个人看视频在线观看www免费| 欧美成人一区二区免费高清观看| 国产成人精品久久久久久| 99热6这里只有精品| 国产精品福利在线免费观看| 日本欧美国产在线视频| 你懂的网址亚洲精品在线观看| 成人无遮挡网站| 寂寞人妻少妇视频99o| 在线免费观看不下载黄p国产| 黄片无遮挡物在线观看| 欧美成人午夜免费资源| 国产一区亚洲一区在线观看| 80岁老熟妇乱子伦牲交| 搡老妇女老女人老熟妇| 精品不卡国产一区二区三区| 亚洲欧美中文字幕日韩二区| 91aial.com中文字幕在线观看|