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

    基于曲率模態(tài)信息熵和BP神經(jīng)網(wǎng)絡(luò)的簡支梁損傷識別方法

    2022-10-06 01:47:16項長生李凌云
    關(guān)鍵詞:信息熵廣義曲率

    項長生, 原 子, 周 宇, 李凌云

    (1. 蘭州理工大學(xué) a. 土木工程學(xué)院, b. 西部土木工程防災(zāi)減災(zāi)教育部工程研究中心, 蘭州 730050; 2. 安徽建筑大學(xué) 土木工程學(xué)院, 合肥 230601)

    結(jié)構(gòu)的損傷識別作為結(jié)構(gòu)健康監(jiān)測中的重要內(nèi)容,能夠為結(jié)構(gòu)剩余壽命評估提供重要參考依據(jù),有較高的研究價值.近年來,有關(guān)結(jié)構(gòu)損傷識別方法的種類繁多,成果豐碩[1-3].很多學(xué)者將結(jié)構(gòu)信息與不同算法相結(jié)合并應(yīng)用于實際工程的損傷診斷中,這印證了損傷識別的指標(biāo)及方法存在更高精度、更高效的趨勢[4-6].在這些方法中,將結(jié)構(gòu)動力指紋與神經(jīng)網(wǎng)絡(luò)算法相結(jié)合的成果較為豐碩[7-8].對于受彎構(gòu)件,曲率模態(tài)及其變形形式具有較高的精度及適用性,且其結(jié)合神經(jīng)網(wǎng)絡(luò)已被證明可以較好地識別損傷,但仍存在以下問題:曲率模態(tài)在靠近振型節(jié)點處精度不高;需要較密集測點,實際工程中難以實現(xiàn);利用曲率模態(tài)差等指標(biāo)進(jìn)行損傷識別需要結(jié)構(gòu)未損傷數(shù)據(jù),對部分既有橋梁難以實現(xiàn).

    信息熵作為控制論等領(lǐng)域的重要概念,在土木工程領(lǐng)域仍處于初步研究階段.Li等[9]將信息熵與結(jié)構(gòu)模態(tài)位移相結(jié)合,利用神經(jīng)網(wǎng)絡(luò)及信息融合技術(shù)提高了損傷識別的精度;程建曙[10]將信息熵用于混凝土裂縫擴(kuò)展機(jī)理的描述;叢培江等[11]闡明了“能量熵”的概念并解釋了其在土木工程應(yīng)用的原理.在此基礎(chǔ)上,Yang等[12-13]將結(jié)構(gòu)的模態(tài)看作波形并將基于模態(tài)位移的廣義局部信息熵應(yīng)用于梁式結(jié)構(gòu)的損傷識別中,結(jié)果證明該指標(biāo)優(yōu)于模態(tài)位移指標(biāo).

    本文針對曲率模態(tài)的不足,結(jié)合信息熵具有的較大研究潛力,對曲率模態(tài)進(jìn)行必要的改進(jìn)和擴(kuò)展,提出“模態(tài)信息熵”的概念使其可以彌補(bǔ)上述缺陷.通過有限元算例結(jié)合神經(jīng)網(wǎng)絡(luò)對損傷進(jìn)行定位、定量以驗證指標(biāo)的有效性,最后通過損傷鋼梁模態(tài)測試,結(jié)合神經(jīng)網(wǎng)絡(luò)進(jìn)行實際結(jié)構(gòu)損傷識別以驗證指標(biāo)的實用性.

    1 廣義局部曲率模態(tài)信息熵

    1.1 指標(biāo)引入及推導(dǎo)

    對于離散變量,設(shè)信號長度為n,信號i的信息量為xi,則信息熵可表示為

    H(X)=E[I(xi)]=E[-logapi]=

    (1)

    式中:I(xi)為xi所代表的自信息量;E[I(xi)]為所有自信息量的數(shù)學(xué)期望;pi為信號i出現(xiàn)的概率;下標(biāo)a通常取e.

    對于廣義局部信息熵,各個信號的和并不等于1,且利用測點i及其周圍測點的振動信號直接代表該測點所包含的信息量.基于以上考慮,并不能將該指標(biāo)稱作經(jīng)典意義上的信息熵,稱其為“廣義”;又因該指標(biāo)側(cè)重于表征測點所包含信息量占附近節(jié)點信息量的比重,因此稱其為“局部”.該指標(biāo)能夠進(jìn)行損傷識別的依據(jù)為:若該信號的長度中各個單位信號所包含信息量相等,不難證明pimax=1/n.該特性反映到結(jié)構(gòu)中則表現(xiàn)為:若結(jié)構(gòu)無損傷,結(jié)構(gòu)各個節(jié)點包含的信息量處處相等,指標(biāo)平滑;若結(jié)構(gòu)出現(xiàn)損傷,則損傷位置出現(xiàn)峰值,可定位損傷.

    根據(jù)“模態(tài)熵”的概念,將結(jié)構(gòu)的模態(tài)位移視作波形,結(jié)合曲率模態(tài)能夠定位、定量損傷的優(yōu)勢,本文構(gòu)造計算點i的廣義局部曲率模態(tài)信息熵的局部概率pi時,取i點和前后緊鄰的點共N個(N為正奇數(shù)且N≥3),并以每個點的曲率模態(tài)vi代表該點的信息量,則此時定義局部概率為

    (2)

    (3)

    式中:yi為i點的模態(tài)位移;h為測點間距.本文利用式(3)計算簡支梁橋的曲率模態(tài).i點的廣義局部曲率模態(tài)信息熵可表示為

    HGLvi=-pilnpi

    (4)

    1.2 概率窗口選取

    式(2)中,N為參與計算的點,宏觀上為滑動計算窗口.N的取值越大,則參與計算的點越多,對于曲率模態(tài),微小損傷造成的結(jié)構(gòu)特性變化有可能會埋沒在龐大的結(jié)構(gòu)信息里而無法表征,且滑動窗口可能包含不同測點的多個損傷信息,這使得窗口中所反映的損傷信息不完整,從而影響損傷識別的精度.同時,N的取值越大,則計算中結(jié)構(gòu)兩端因數(shù)據(jù)缺失,不滿足式(2)計算要求的點越多.綜上,N的取值不宜過大.本文取N=3,可得

    (5)

    聯(lián)合式(3)~(5)并結(jié)合信息熵定義,可得到結(jié)構(gòu)模態(tài)位移(各階振型分量)與廣義局部曲率模態(tài)信息熵的關(guān)系為

    (6)

    本文利用式(6)計算前兩階廣義局部曲率模態(tài)信息熵并以此指標(biāo)對結(jié)構(gòu)進(jìn)行損傷識別.簡支梁的二階豎向振型跨中存在振型節(jié)點,其存在使得損傷識別易發(fā)生誤判,因此有必要對該條件下?lián)p傷識別的指標(biāo)進(jìn)行改進(jìn).

    1.3 指標(biāo)改進(jìn)

    本文針對二階振型即跨中存在振型節(jié)點的振型進(jìn)行改進(jìn),探究利用其進(jìn)行損傷識別的可行性.限于篇幅,簡支梁的三階及更高階廣義局部曲率模態(tài)信息熵的計算方法及損傷識別可針對本文研究內(nèi)容進(jìn)行類比,此處不再贅述.

    利用這類振型數(shù)據(jù)進(jìn)行損傷識別時,根據(jù)式(4)、(5)則有

    (7)

    (8)

    式中:t為有限元單元劃分的個數(shù);p0為跨中振型節(jié)點處的局部概率;HGLv0為該處的廣義局部曲率模態(tài)信息熵.由式(8)可知,若曲率模態(tài)包含此類振型節(jié)點,則廣義局部曲率模態(tài)信息熵會有明顯的趨于0的峰值,容易使判別損傷位置時發(fā)生誤判.為了避免該情況發(fā)生,需對曲率模態(tài)進(jìn)行處理,處理的依據(jù)為:曲率模態(tài)無單位且為承彎結(jié)構(gòu)振動特性的特殊表現(xiàn)形式,其表達(dá)式為

    (9)

    則有

    (10)

    式中:y′i為第i點的結(jié)構(gòu)豎向振型分量;v′i為第i階變形前的結(jié)構(gòu)曲率模態(tài).處理后能夠有效地規(guī)避二階指標(biāo)產(chǎn)生趨于0的峰值,且移軸后結(jié)構(gòu)的振動趨勢未發(fā)生改變.

    2 簡支梁數(shù)值模擬算例

    2.1 MIDAS模型建立

    以一混凝土簡支梁橋為例,該橋的全長L=32 m,設(shè)計截面為箱形截面,如圖1a所示(單位:mm),單元屬性為梁單元,其中主梁結(jié)構(gòu)采用C50混凝土,彈性模量為3.45×107kN/m2,容重為25 kN/m3,泊松比為0.2.利用有限元軟件MIDAS/Civil建立簡支梁有限元模型,如圖1b所示,該簡支梁沿長度方向劃分為等長的32個單元.

    圖1 簡支梁橋幾何構(gòu)造及分析模型Fig.1 Geometric structure and analysis model of simply supported beam bridge

    為了驗證指標(biāo)的有效性,利用剛度折減法,通過彈性模量E的下降來模擬損傷,三個工況的損傷位置分別位于跨中、1/4跨和1/8跨,對應(yīng)的損傷單元號為4、8、16,損傷程度均為20%,如表1所示.

    表1 有限元模型損傷工況Tab.1 Damage conditions of finite element model

    2.2 指標(biāo)損傷定位結(jié)果

    簡支梁橋模型中有共計29個可計算測點.提取一階豎向振型并利用式(6)計算工況1-1至1-3結(jié)構(gòu)各個節(jié)點的前兩階廣義局部曲率模態(tài)信息熵,圖2為前兩階廣義局部曲率模態(tài)信息熵?fù)p傷識別結(jié)果.

    圖2 前兩階廣義局部曲率模態(tài)信息熵?fù)p傷識別結(jié)果Fig.2 Damage identification results of first two order generalized local curvature modal information entropy

    對于簡支梁,一階振型節(jié)點位于兩端支座處,二階振型節(jié)點位于簡支梁跨中.二階廣義局部曲率模態(tài)信息熵在該階振型節(jié)點(跨中)處損傷識別效果優(yōu)于一階指標(biāo),其峰值更明顯.但隨著損傷位置遠(yuǎn)離二階振型節(jié)點,靠近一階振型節(jié)點即支座處,一階廣義局部曲率模態(tài)信息熵峰值較二階更為明顯.說明不同階數(shù)的廣義局部曲率模態(tài)信息熵具有相同的性質(zhì),自階振型節(jié)點處對損傷較其他位置敏感.

    在跨中位置,前兩階指標(biāo)均可以準(zhǔn)確地定位損傷.在1/4跨及靠近支座的1/8跨處,二階指標(biāo)峰值較實際損傷位置有一個測點的偏移,但仍可以確定損傷的大致范圍.隨著單元網(wǎng)格劃分精度、單元數(shù)量的提升,可以減少該現(xiàn)象對損傷識別的影響;同時,二階廣義局部曲率模態(tài)信息熵較一階更為平滑.

    2.3 基于神經(jīng)網(wǎng)絡(luò)的損傷定量方法及驗證

    對工況1-1至1-3進(jìn)行損傷定量時采用三層BP神經(jīng)網(wǎng)絡(luò),根據(jù)損傷定位結(jié)果,取已知損傷位置附近3個節(jié)點的廣義局部曲率模態(tài)信息熵和曲率模態(tài)作為神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù),并構(gòu)造三種損傷位置下?lián)p傷程度為10%、15%、25%、30%的共計12個工況作為兩種指標(biāo)的訓(xùn)練樣本.根據(jù)三分法[14]試算神經(jīng)網(wǎng)絡(luò)的最優(yōu)隱含層節(jié)點數(shù)為5,輸出參數(shù)設(shè)置為損傷單元的殘余剛度[οi],神經(jīng)網(wǎng)絡(luò)的構(gòu)造均為3-5-1,各個網(wǎng)絡(luò)訓(xùn)練最大迭代步數(shù)均為50,訓(xùn)練精度均為10-7.得到的損傷識別結(jié)果如表2所示.

    由表2可知,在同樣的神經(jīng)網(wǎng)絡(luò)參數(shù)下,前兩階曲率模態(tài)損傷定量的精度隨著損傷位置靠近振型節(jié)點而降低.對于一階曲率模態(tài),在靠近支座的1/4跨和1/8跨處,損傷識別的相對誤差達(dá)到了4.88%和8.25%;廣義局部曲率模態(tài)信息熵?fù)p傷定量相對誤差僅為4.13%和1.38%;二階振型節(jié)點即簡支梁跨中處二階廣義局部曲率模態(tài)信息熵?fù)p傷程度識別結(jié)果誤差較小,相對誤差僅為1.1%,但若損傷位置在遠(yuǎn)離振型節(jié)點的1/8跨處,其損傷程度識別相對誤差達(dá)到了9.12%,高于曲率模態(tài).因此,兩種指標(biāo)在識別位置不同的損傷時有互補(bǔ)的特點.

    表2 神經(jīng)網(wǎng)絡(luò)損傷定量識別結(jié)果Tab.2 Quantitative damage identification results of BP neural network

    2.4 測點數(shù)量與指標(biāo)識別效果的關(guān)系

    利用曲率模態(tài)進(jìn)行損傷識別時往往需要數(shù)量較多的測點,以便能夠利用差分法進(jìn)行計算.為了對實際工程及模態(tài)測試的測點布置提供一定參考,以簡支梁模型為數(shù)值算例,將模型單元精細(xì)劃分為64和128個,以模擬測點數(shù)量的增加.以一階振型為研究對象,不同損傷程度下簡支梁損傷識別結(jié)果如圖3所示.

    圖3 不同測點數(shù)下1/8跨損傷識別結(jié)果Fig.3 Damage identification results of 1/8 span with different number of measuring points

    由圖3可以看出,不同損傷位置下,節(jié)點在數(shù)量為33時損傷識別結(jié)果清晰,損傷位置處表征損傷的指標(biāo)峰值較為明顯.隨著損傷程度的增加,其熵值減小并呈一定規(guī)律;但在測點數(shù)量為129時,損傷處有無法計算的熵值,且峰值處數(shù)值相差較大,損傷指標(biāo)較不規(guī)律;能夠?qū)p傷進(jìn)行定位,但無法進(jìn)行定量損傷.指標(biāo)在測點數(shù)量增多時呈現(xiàn)的不規(guī)律性對利用神經(jīng)網(wǎng)絡(luò)判定其損傷程度較為不利.本文中測點數(shù)為33時識別效果最佳.

    綜上所述,廣義局部曲率模態(tài)信息熵在測點數(shù)量的選擇上不同于曲率模態(tài),該指標(biāo)能夠較好避免曲率模態(tài)需要密集測點這一缺點,即無需密集測點即可獲得較好損傷識別結(jié)果.

    3 試驗驗證

    3.1 試驗?zāi)康?/h3>

    為了進(jìn)一步對本文所提出的廣義局部曲率模態(tài)信息熵、損傷定量識別方法、精度及應(yīng)用于實際工程中的識別結(jié)果進(jìn)行考量和驗證,基于上文對新?lián)p傷指標(biāo)的研究和總結(jié),本文將該指標(biāo)應(yīng)用于實際結(jié)構(gòu),通過對比試驗和數(shù)值模擬識別結(jié)果驗證該方法的可行性和實用性,以期為該指標(biāo)及方法在實際工程中的應(yīng)用提供一定的建議.

    3.2 試驗對象及測試設(shè)備

    以鋼梁為試驗對象,鋼梁橫截面尺寸及設(shè)計為:鋼梁梁體總長2 350 mm,計算跨徑1 750 mm,凈跨徑1 650 mm,截面尺寸100 mm×8 mm.為了與試驗室支墩緊密貼合,距離鋼梁兩端250 mm處設(shè)有長度為100 mm,截面尺寸為160 mm×8 mm的平面,螺栓固定位置位于鋼梁凈跨徑兩側(cè),預(yù)制螺孔以固定鋼梁,與下方支座用螺栓連接.兩支座中心間距為計算跨徑,鋼梁具體尺寸如圖4所示(單位:mm),命名為鋼梁-1.對鋼梁進(jìn)行材性試驗,測得該鋼材的平均屈服強(qiáng)度約為432 MPa;通過材性測定數(shù)據(jù),計算其彈性模量約為1.989×105N/mm2.

    圖4 鋼梁尺寸Fig.4 Steel beam dimension

    除去支墩、支座底板及螺栓若干,采集系統(tǒng)及加速度傳感器如圖5所示.

    圖5 加速度傳感器及采集分析系統(tǒng)Fig.5 Acceleration sensor and acquisition analysis system

    3.3 試驗工況設(shè)置

    鋼梁-1為預(yù)制無損梁,無損梁測試完畢后返廠加工并于30 cm處(以支座中心為原點)開口2.5 cm,定為損傷工況3-1,測試完畢后在87.5 cm處開口2.5 cm,定為損傷工況3-2.實際開口形式如表3所示.

    表 3 鋼梁損傷工況Tab.3 Damage conditions of steel beam

    3.4 測試方法及測點布置

    測點無需布置過于密集,但其數(shù)量應(yīng)保證能夠識別損傷,結(jié)合傳感器底座大小及鋼梁實際情況,測點間距設(shè)為5 cm,進(jìn)行數(shù)據(jù)分析時假定兩支座中心處振型分量與理論值相等,值為0.全梁共34個測點,其中測點1位于距固定鉸支座中心5 cm處.測點34位于距滾動支座中心5 cm處,其余各個測點均勻分布.

    試驗采用錘擊法進(jìn)行模態(tài)數(shù)據(jù)采集.試驗室共有6個高頻加速度傳感器,為了測得同一工況下結(jié)構(gòu)完整的振型,以5個測點及1個參考點為一組進(jìn)行測試,共需7組測試.通過力錘對參考點附近進(jìn)行激勵以觸發(fā)采樣,每組采樣采集三次,取平均值進(jìn)行分析.

    試算結(jié)果表明,結(jié)構(gòu)的一階振型峰值最大處在跨中部位,但該位置恰好為二階振型節(jié)點.結(jié)構(gòu)的二階振型在3/4處達(dá)到峰值且與一階振型方向相反.基于上述考慮,本文試驗中獲取結(jié)構(gòu)的一階振型時,將參考點設(shè)立在測點18處.獲取結(jié)構(gòu)二階振型時,則將參考點設(shè)立在靠近3/4處(131.25 mm)附近的測點26處.具體的測點分布形式如圖6所示(單位:mm).

    圖6 模態(tài)測試測點布置圖Fig.6 Layout of modal measuring points

    針對前兩階模態(tài),加速度響應(yīng)采集完畢后,利用隨機(jī)子空間分析法分析采集到的波形時域數(shù)據(jù),以參考點為基準(zhǔn)進(jìn)行振型歸一化和質(zhì)量歸一化處理,提取前兩階豎向振型分量.最后對數(shù)據(jù)進(jìn)行濾波預(yù)處理,確定損傷識別結(jié)果.試驗采集現(xiàn)場及采集過程如圖7所示.

    圖7 模態(tài)測試現(xiàn)場Fig.7 Modal test site

    3.5 試驗結(jié)構(gòu)有限元模型

    為了與實際試驗?zāi)P瓦M(jìn)行有效對比分析,確定試驗數(shù)據(jù)的準(zhǔn)確性,采用MIDAS/Civil和ANSYS兩種有限元軟件根據(jù)試驗對象的實際尺寸建立相應(yīng)有限元模型.

    MIDAS/Civil全模型沿梁方向共計235個長度為1 cm的梁單元,共2 350 mm.同時利用ANSYS進(jìn)行實體單元建模,單元類型為solid185.選用六面體單元對模型進(jìn)行掃掠劃分,網(wǎng)格大小選用4 mm.根據(jù)材性試驗結(jié)果對兩模型中單元賦予材料屬性,彈性模量為1.989×105N/mm2,泊松比為0.3,容重為7.698×10-5N/mm3.具體模型如圖8所示(以工況3-1為例).

    圖8 有限元損傷模型Fig.8 FEM damage model

    分別進(jìn)行模態(tài)分析,得到結(jié)構(gòu)的前兩節(jié)實測振型,并將ANSYS、MIDAS和實測振型進(jìn)行歸一化處理.將得到的振型代入式(3)~(10)計算結(jié)構(gòu)的前兩階曲率模態(tài)及廣義局部曲率模態(tài)信息熵.損傷識別結(jié)果如圖9所示.

    圖9 指標(biāo)實測、理論損傷位置Fig.9 Actual and theoretical damage index locations

    損傷定量神經(jīng)網(wǎng)絡(luò)依舊采用三層BP神經(jīng)網(wǎng)絡(luò),并構(gòu)造兩個損傷位置下開口深度為1.5、2.0、2.5、3.0 cm(含同時損傷)共計18個工況作為兩種指標(biāo)的訓(xùn)練樣本.各項參數(shù)的選取與上一節(jié)神經(jīng)網(wǎng)絡(luò)相同,輸出函數(shù)為該單元的殘余剛度[oi].確定隱含層最優(yōu)節(jié)點數(shù)為11,因此,神經(jīng)網(wǎng)絡(luò)的構(gòu)造為3-11-1.損傷定量結(jié)果如表4所示.

    表4 簡支鋼梁損傷定量識別結(jié)果Tab.4 Quantitative damage identification results of simply supported beam

    對比前兩階指標(biāo)可知,一、二階曲率模態(tài)指標(biāo)在各自階振型節(jié)點處的損傷定量識別精度均低于廣義局部曲率模態(tài)信息熵,靠近一階振型節(jié)點即30 cm處,定量誤差達(dá)到72.90%,靠近二階振型節(jié)點即87.5 cm處,曲率模態(tài)誤差達(dá)到34.80%,均高于本文提出指標(biāo).對于兩種損傷工況,指標(biāo)各自階振型節(jié)點處精度均低于遠(yuǎn)離振型節(jié)點處的損傷識別,這對前文結(jié)論進(jìn)行了驗證.

    4 結(jié) 論

    簡支梁數(shù)值算例中,本文提出的廣義局部曲率模態(tài)信息熵能夠較好進(jìn)行損傷定位.結(jié)合BP神經(jīng)網(wǎng)絡(luò)定量損傷時,對比曲率模態(tài),廣義局部曲率模態(tài)信息熵在判別位置靠近各階振型節(jié)點的損傷時,具有更高的精度.曲率模態(tài)在遠(yuǎn)離振型節(jié)點處識別精度較高,兩指標(biāo)互補(bǔ).本文提出的指標(biāo)具有無需密集測點即可得到效果良好損傷識別結(jié)果的優(yōu)勢.相比其他位置,指標(biāo)在識別振型節(jié)點附近損傷時,對測點數(shù)量的改變較為敏感.

    通過設(shè)置不同位置的參考點對鋼梁前兩階振型進(jìn)行了模態(tài)測試.限于試驗儀器精度及測試環(huán)境影響,本文所提指標(biāo)結(jié)合BP神經(jīng)網(wǎng)絡(luò)能夠進(jìn)行大致?lián)p傷定量,且在振型節(jié)點處指標(biāo)識別精度優(yōu)于曲率模態(tài),驗證了指標(biāo)在實際結(jié)構(gòu)中的實用性.

    在實際工程中,若損傷位置可測,可通過實際損傷位置,依據(jù)曲率模態(tài)和廣義局部曲率模態(tài)信息熵的互補(bǔ)特性選定采用的結(jié)構(gòu)模態(tài)采集階數(shù)和損傷指標(biāo).本文提出的指標(biāo)不需要結(jié)構(gòu)的初始模態(tài)即可進(jìn)行損傷識別,加之工程結(jié)構(gòu)中難以獲得精度較高的高階模態(tài)振型,使本文提出的方法更具有實際工程意義.

    猜你喜歡
    信息熵廣義曲率
    大曲率沉管安裝關(guān)鍵技術(shù)研究
    一類雙曲平均曲率流的對稱與整體解
    基于信息熵可信度的測試點選擇方法研究
    Rn中的廣義逆Bonnesen型不等式
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    從廣義心腎不交論治慢性心力衰竭
    基于信息熵的實驗教學(xué)量化研究
    電子測試(2017年12期)2017-12-18 06:35:48
    一種基于信息熵的雷達(dá)動態(tài)自適應(yīng)選擇跟蹤方法
    有限群的廣義交換度
    基于信息熵的IITFN多屬性決策方法
    伊人久久精品亚洲午夜| 国产av不卡久久| 男女国产视频网站| 日本av手机在线免费观看| av在线播放精品| 中文字幕av成人在线电影| 国产熟女欧美一区二区| 成人亚洲欧美一区二区av| 欧美激情久久久久久爽电影| 亚洲在久久综合| 18+在线观看网站| 亚洲av中文字字幕乱码综合| 国产爱豆传媒在线观看| 97超视频在线观看视频| 久久午夜福利片| 久久草成人影院| 欧美 日韩 精品 国产| 97超视频在线观看视频| 久久精品国产亚洲av涩爱| 色尼玛亚洲综合影院| 欧美高清性xxxxhd video| 免费黄色在线免费观看| 91久久精品国产一区二区成人| 久久久久久久久大av| 久久精品夜夜夜夜夜久久蜜豆| av在线观看视频网站免费| 国产一级毛片七仙女欲春2| 美女脱内裤让男人舔精品视频| 中国美白少妇内射xxxbb| 成人午夜精彩视频在线观看| 男女视频在线观看网站免费| 国产精品精品国产色婷婷| 精品一区在线观看国产| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 久久国产乱子免费精品| 久久6这里有精品| 国产麻豆成人av免费视频| 女的被弄到高潮叫床怎么办| 日本熟妇午夜| 日本与韩国留学比较| 亚洲国产欧美在线一区| videossex国产| 久久精品久久久久久久性| 夜夜爽夜夜爽视频| 欧美xxⅹ黑人| 国产成人福利小说| 精品一区二区三区视频在线| 精品人妻一区二区三区麻豆| 亚洲美女搞黄在线观看| 国产伦理片在线播放av一区| 日韩中字成人| 免费观看性生交大片5| 欧美丝袜亚洲另类| 久久久久免费精品人妻一区二区| 亚洲欧美精品专区久久| 成人漫画全彩无遮挡| 五月天丁香电影| 国产免费又黄又爽又色| 永久网站在线| 99久久中文字幕三级久久日本| 嫩草影院新地址| 好男人在线观看高清免费视频| 亚洲人成网站高清观看| 99视频精品全部免费 在线| 97人妻精品一区二区三区麻豆| 偷拍熟女少妇极品色| 人妻少妇偷人精品九色| 精品久久久久久久人妻蜜臀av| 午夜福利网站1000一区二区三区| 精品国产三级普通话版| 成人午夜高清在线视频| 国产在线一区二区三区精| 亚洲人成网站在线观看播放| 99热这里只有是精品50| 中文天堂在线官网| 国产精品一区二区三区四区免费观看| 国产精品久久久久久久电影| 成人毛片a级毛片在线播放| 国产免费福利视频在线观看| 国产亚洲5aaaaa淫片| 成人亚洲欧美一区二区av| 精品国产露脸久久av麻豆 | a级毛色黄片| 久久这里只有精品中国| 看非洲黑人一级黄片| 精品人妻偷拍中文字幕| 中文字幕免费在线视频6| 国产精品日韩av在线免费观看| 日日干狠狠操夜夜爽| 91久久精品国产一区二区三区| 亚洲精品久久午夜乱码| 午夜老司机福利剧场| 男人爽女人下面视频在线观看| 国产精品蜜桃在线观看| 日本一本二区三区精品| 日韩成人伦理影院| 日日摸夜夜添夜夜添av毛片| 十八禁网站网址无遮挡 | 欧美丝袜亚洲另类| 成人鲁丝片一二三区免费| 成人亚洲精品av一区二区| 亚洲av在线观看美女高潮| 真实男女啪啪啪动态图| 人人妻人人看人人澡| 内射极品少妇av片p| 欧美日韩一区二区视频在线观看视频在线 | 有码 亚洲区| 高清av免费在线| 色综合色国产| 久久精品夜夜夜夜夜久久蜜豆| 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| 亚洲国产高清在线一区二区三| 日韩制服骚丝袜av| 欧美性猛交╳xxx乱大交人| 国产精品一区www在线观看| 国产老妇伦熟女老妇高清| 自拍偷自拍亚洲精品老妇| 极品少妇高潮喷水抽搐| 久久精品熟女亚洲av麻豆精品 | 97超碰精品成人国产| 可以在线观看毛片的网站| 成人午夜精彩视频在线观看| 欧美人与善性xxx| 国产精品一区二区三区四区免费观看| 亚洲国产高清在线一区二区三| 一区二区三区高清视频在线| 亚洲精品日韩av片在线观看| 亚洲av男天堂| 国产黄a三级三级三级人| 91在线精品国自产拍蜜月| 欧美丝袜亚洲另类| 国产精品.久久久| 婷婷色综合大香蕉| 一本久久精品| 亚洲精品一二三| 亚洲精品久久久久久婷婷小说| 激情五月婷婷亚洲| 亚洲美女搞黄在线观看| 身体一侧抽搐| 久久久久久国产a免费观看| 成人av在线播放网站| av专区在线播放| 日产精品乱码卡一卡2卡三| 精品一区在线观看国产| 久热久热在线精品观看| 日本-黄色视频高清免费观看| av又黄又爽大尺度在线免费看| 一本一本综合久久| 18+在线观看网站| 亚洲欧美精品专区久久| 一夜夜www| 直男gayav资源| 九九在线视频观看精品| 可以在线观看毛片的网站| 2021天堂中文幕一二区在线观| 男女边吃奶边做爰视频| 久久99热6这里只有精品| 欧美人与善性xxx| av国产免费在线观看| 欧美日韩国产mv在线观看视频 | 国模一区二区三区四区视频| 亚洲图色成人| 伦理电影大哥的女人| 永久网站在线| 中文资源天堂在线| 九九在线视频观看精品| 26uuu在线亚洲综合色| 色哟哟·www| 日本免费a在线| 乱系列少妇在线播放| 久久久a久久爽久久v久久| 精品人妻偷拍中文字幕| 秋霞伦理黄片| 看十八女毛片水多多多| 日本色播在线视频| 中文精品一卡2卡3卡4更新| 在现免费观看毛片| av线在线观看网站| 中文资源天堂在线| 三级国产精品欧美在线观看| ponron亚洲| 天堂√8在线中文| 天天躁日日操中文字幕| 亚洲精品成人av观看孕妇| 观看免费一级毛片| 国产精品99久久久久久久久| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 欧美性猛交╳xxx乱大交人| 国产精品熟女久久久久浪| 如何舔出高潮| 中文字幕人妻熟人妻熟丝袜美| 看十八女毛片水多多多| 汤姆久久久久久久影院中文字幕 | 少妇高潮的动态图| 精品人妻视频免费看| 综合色av麻豆| 美女cb高潮喷水在线观看| 如何舔出高潮| 国产精品不卡视频一区二区| av天堂中文字幕网| 99久久精品热视频| 在线观看美女被高潮喷水网站| 亚洲精品成人av观看孕妇| 又黄又爽又刺激的免费视频.| 少妇高潮的动态图| 欧美最新免费一区二区三区| 精品熟女少妇av免费看| 亚洲欧美日韩无卡精品| 美女黄网站色视频| 两个人视频免费观看高清| 国产淫片久久久久久久久| 国产高清三级在线| 人体艺术视频欧美日本| 日本色播在线视频| 国产一区二区亚洲精品在线观看| 97热精品久久久久久| 中文字幕制服av| 国产毛片a区久久久久| 91aial.com中文字幕在线观看| 夫妻午夜视频| 精品欧美国产一区二区三| av在线亚洲专区| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说 | 日韩欧美三级三区| 成人无遮挡网站| av线在线观看网站| 国产精品嫩草影院av在线观看| 最后的刺客免费高清国语| 免费在线观看成人毛片| 亚洲人成网站在线观看播放| 国产免费福利视频在线观看| 久久久久久伊人网av| 69人妻影院| 少妇熟女欧美另类| 亚洲三级黄色毛片| 老司机影院毛片| 婷婷色综合大香蕉| 美女xxoo啪啪120秒动态图| 乱系列少妇在线播放| 成人美女网站在线观看视频| 麻豆久久精品国产亚洲av| 黑人高潮一二区| 欧美成人一区二区免费高清观看| 中文字幕制服av| 日韩国内少妇激情av| 国产国拍精品亚洲av在线观看| 国产精品国产三级国产专区5o| 99久久九九国产精品国产免费| 久久久精品94久久精品| av黄色大香蕉| 亚洲精品久久久久久婷婷小说| 欧美日韩国产mv在线观看视频 | 少妇人妻一区二区三区视频| 岛国毛片在线播放| 神马国产精品三级电影在线观看| 有码 亚洲区| 久久久久九九精品影院| 亚洲国产成人一精品久久久| 嫩草影院精品99| 国产亚洲av嫩草精品影院| 高清在线视频一区二区三区| 国产成人aa在线观看| 久久精品熟女亚洲av麻豆精品 | 国产午夜福利久久久久久| 日韩三级伦理在线观看| 日韩av在线大香蕉| 欧美3d第一页| 内射极品少妇av片p| 国产成年人精品一区二区| 秋霞在线观看毛片| 一级毛片黄色毛片免费观看视频| 青青草视频在线视频观看| videos熟女内射| 只有这里有精品99| 亚洲成人久久爱视频| 色综合色国产| 黑人高潮一二区| 成年av动漫网址| 男的添女的下面高潮视频| or卡值多少钱| 精品亚洲乱码少妇综合久久| 少妇熟女aⅴ在线视频| 国产成人精品久久久久久| 亚洲丝袜综合中文字幕| 免费看光身美女| 亚洲激情五月婷婷啪啪| 美女黄网站色视频| 能在线免费看毛片的网站| 亚洲人成网站高清观看| 亚洲熟女精品中文字幕| 97精品久久久久久久久久精品| 亚洲欧美成人精品一区二区| 五月伊人婷婷丁香| 丰满少妇做爰视频| 日韩电影二区| 久久久久久伊人网av| 最近2019中文字幕mv第一页| 亚洲成色77777| 嫩草影院精品99| 真实男女啪啪啪动态图| 99久久精品热视频| 国产精品国产三级国产专区5o| 1000部很黄的大片| 亚洲精品日本国产第一区| 久久亚洲国产成人精品v| 日韩av免费高清视频| 中文字幕亚洲精品专区| 亚洲精品日韩av片在线观看| 亚洲精品第二区| 一个人看视频在线观看www免费| 国产精品三级大全| av在线老鸭窝| 亚洲美女视频黄频| 亚洲一区高清亚洲精品| 久久热精品热| 国产精品一区www在线观看| 夫妻性生交免费视频一级片| 嫩草影院入口| 日韩国内少妇激情av| 亚洲av成人精品一二三区| 欧美日韩亚洲高清精品| 午夜福利在线在线| 国产亚洲一区二区精品| 成人毛片a级毛片在线播放| 久久6这里有精品| 99re6热这里在线精品视频| 中文精品一卡2卡3卡4更新| 免费看日本二区| 免费看a级黄色片| 男插女下体视频免费在线播放| 不卡视频在线观看欧美| 淫秽高清视频在线观看| 午夜精品一区二区三区免费看| 亚洲av男天堂| 国产久久久一区二区三区| 日产精品乱码卡一卡2卡三| 最近2019中文字幕mv第一页| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| 大香蕉久久网| 青青草视频在线视频观看| 亚洲欧美成人精品一区二区| 亚洲国产欧美人成| 国产在线男女| 777米奇影视久久| 99久久人妻综合| 最新中文字幕久久久久| 91精品国产九色| 人人妻人人澡欧美一区二区| 秋霞在线观看毛片| 精品人妻偷拍中文字幕| 国产成人aa在线观看| 全区人妻精品视频| 丝袜美腿在线中文| h日本视频在线播放| 一级毛片久久久久久久久女| 国产伦理片在线播放av一区| 久久久久久久午夜电影| 久久久久久久亚洲中文字幕| 一级毛片 在线播放| 午夜福利成人在线免费观看| 搡女人真爽免费视频火全软件| 日本黄色片子视频| 神马国产精品三级电影在线观看| 国产乱来视频区| 91久久精品电影网| 一级黄片播放器| 国产一级毛片七仙女欲春2| 青春草亚洲视频在线观看| 小蜜桃在线观看免费完整版高清| 美女主播在线视频| 嫩草影院新地址| 亚洲欧美成人综合另类久久久| 天天躁日日操中文字幕| 欧美成人午夜免费资源| 熟妇人妻久久中文字幕3abv| 日韩不卡一区二区三区视频在线| 春色校园在线视频观看| 日韩电影二区| 国产精品无大码| av在线天堂中文字幕| 久久久久久久久久久丰满| 国产高潮美女av| 18禁在线播放成人免费| 国产久久久一区二区三区| 纵有疾风起免费观看全集完整版 | 有码 亚洲区| 国产乱人偷精品视频| 亚洲乱码一区二区免费版| 久久久久精品性色| 亚洲欧美一区二区三区黑人 | 亚洲欧美日韩东京热| 听说在线观看完整版免费高清| 中文资源天堂在线| 18禁在线播放成人免费| 免费看不卡的av| 一级a做视频免费观看| 国产精品麻豆人妻色哟哟久久 | 高清在线视频一区二区三区| av又黄又爽大尺度在线免费看| videos熟女内射| 在线观看一区二区三区| 在线观看人妻少妇| 美女国产视频在线观看| 亚洲人成网站高清观看| 国产精品人妻久久久影院| 亚洲av成人精品一区久久| 国产精品一及| 精品99又大又爽又粗少妇毛片| 亚洲av国产av综合av卡| 欧美zozozo另类| 国产麻豆成人av免费视频| 中文字幕制服av| 国产有黄有色有爽视频| 国产亚洲午夜精品一区二区久久 | 国产成人午夜福利电影在线观看| 精华霜和精华液先用哪个| 国产人妻一区二区三区在| 色综合色国产| 2018国产大陆天天弄谢| 成年免费大片在线观看| 欧美区成人在线视频| 日韩视频在线欧美| av播播在线观看一区| 美女内射精品一级片tv| 中文字幕制服av| 国产精品国产三级专区第一集| 亚洲精品久久午夜乱码| 真实男女啪啪啪动态图| 亚洲乱码一区二区免费版| 亚洲国产成人一精品久久久| 日韩伦理黄色片| 乱码一卡2卡4卡精品| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站| 亚洲国产精品专区欧美| 大又大粗又爽又黄少妇毛片口| 激情五月婷婷亚洲| 免费大片黄手机在线观看| av在线播放精品| 欧美区成人在线视频| 丰满少妇做爰视频| 人体艺术视频欧美日本| 亚洲国产精品成人久久小说| 天天躁夜夜躁狠狠久久av| 干丝袜人妻中文字幕| 亚洲成人精品中文字幕电影| 午夜福利成人在线免费观看| 亚洲欧美成人综合另类久久久| 久久草成人影院| 亚洲成人精品中文字幕电影| 日日撸夜夜添| 久久午夜福利片| 三级国产精品片| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 日本免费在线观看一区| 国产片特级美女逼逼视频| 午夜福利高清视频| 大香蕉97超碰在线| 亚洲第一区二区三区不卡| 免费观看的影片在线观看| 久久鲁丝午夜福利片| 亚洲国产精品成人综合色| 国产单亲对白刺激| 一区二区三区高清视频在线| 18禁在线播放成人免费| 欧美zozozo另类| 男人爽女人下面视频在线观看| 内射极品少妇av片p| 777米奇影视久久| 欧美极品一区二区三区四区| 国产激情偷乱视频一区二区| 日韩中字成人| 2021少妇久久久久久久久久久| 十八禁网站网址无遮挡 | 国产伦精品一区二区三区四那| 26uuu在线亚洲综合色| av女优亚洲男人天堂| 欧美激情久久久久久爽电影| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 69人妻影院| 成人亚洲精品一区在线观看 | 18禁动态无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲一区高清亚洲精品| 2021天堂中文幕一二区在线观| 亚洲精品中文字幕在线视频 | 我的老师免费观看完整版| 简卡轻食公司| 亚洲欧美精品自产自拍| 亚洲精品第二区| 亚洲成人精品中文字幕电影| 国产人妻一区二区三区在| freevideosex欧美| 成人二区视频| 国产精品av视频在线免费观看| 亚洲精品一区蜜桃| 在线免费观看的www视频| 精品一区二区三卡| 激情 狠狠 欧美| 三级经典国产精品| 插阴视频在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 777米奇影视久久| 亚洲高清免费不卡视频| 亚洲不卡免费看| 国产黄片美女视频| 又粗又硬又长又爽又黄的视频| 中文字幕免费在线视频6| 看十八女毛片水多多多| 久久久久久伊人网av| 三级毛片av免费| 欧美最新免费一区二区三区| 国产精品爽爽va在线观看网站| av在线老鸭窝| 大陆偷拍与自拍| 午夜福利视频精品| 国产伦理片在线播放av一区| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 久久精品人妻少妇| 亚洲aⅴ乱码一区二区在线播放| 亚洲丝袜综合中文字幕| 日韩精品有码人妻一区| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 蜜桃亚洲精品一区二区三区| 国产精品一区二区三区四区久久| 国产精品人妻久久久久久| 欧美成人一区二区免费高清观看| 国内精品美女久久久久久| 成年人午夜在线观看视频 | 少妇的逼好多水| 久久亚洲国产成人精品v| 69人妻影院| 全区人妻精品视频| 免费观看的影片在线观看| 日韩人妻高清精品专区| 婷婷色av中文字幕| 黄片wwwwww| 天堂影院成人在线观看| 深爱激情五月婷婷| 免费av不卡在线播放| 69av精品久久久久久| 中文字幕av在线有码专区| eeuss影院久久| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 亚洲av中文av极速乱| 在线 av 中文字幕| 亚洲成人av在线免费| 国产精品久久视频播放| 午夜久久久久精精品| 日韩大片免费观看网站| 看黄色毛片网站| 久久99精品国语久久久| 91精品国产九色| 男人舔女人下体高潮全视频| 亚洲天堂国产精品一区在线| 美女主播在线视频| 亚洲av电影不卡..在线观看| 久久久久久久久大av| 欧美 日韩 精品 国产| 97人妻精品一区二区三区麻豆| 久久久色成人| 久久久久久久久久人人人人人人| 免费大片18禁| 久久这里有精品视频免费| 高清日韩中文字幕在线| 啦啦啦啦在线视频资源| 国产女主播在线喷水免费视频网站 | 床上黄色一级片| 99热这里只有精品一区| 在线观看一区二区三区| 禁无遮挡网站| 亚洲性久久影院| 少妇裸体淫交视频免费看高清| 国产成人freesex在线| 91久久精品电影网| 午夜精品在线福利| 美女脱内裤让男人舔精品视频| 亚洲熟妇中文字幕五十中出| 联通29元200g的流量卡| 高清av免费在线| 国产精品人妻久久久久久| 久久久精品94久久精品| 可以在线观看毛片的网站| eeuss影院久久| 99热这里只有精品一区| 亚洲精品亚洲一区二区| 少妇裸体淫交视频免费看高清| 免费看光身美女| 亚洲熟妇中文字幕五十中出| 国产一区有黄有色的免费视频 | 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 国产成人aa在线观看| 久久精品人妻少妇| 午夜激情福利司机影院| 特级一级黄色大片| 深爱激情五月婷婷| 久久久欧美国产精品| 国国产精品蜜臀av免费| 天堂网av新在线| 国产成年人精品一区二区| 黄色欧美视频在线观看| 国产伦精品一区二区三区视频9| 中文字幕久久专区|