• <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嫩草精品影院| 脱女人内裤的视频| 全区人妻精品视频| 深夜精品福利| 禁无遮挡网站| 日韩欧美国产在线观看| 成人国产综合亚洲| avwww免费| 九九热线精品视视频播放| 欧美乱妇无乱码| 欧美另类亚洲清纯唯美| 久久久久国产精品人妻aⅴ院| 国产三级在线视频| 一夜夜www| 99热6这里只有精品| 怎么达到女性高潮| 又粗又爽又猛毛片免费看| 99热只有精品国产| 久9热在线精品视频| 伦理电影免费视频| 在线a可以看的网站| 色噜噜av男人的天堂激情| 国产伦人伦偷精品视频| 91麻豆av在线| 国产精品免费视频内射| 午夜免费观看网址| 亚洲成av人片免费观看| 亚洲 欧美一区二区三区| 全区人妻精品视频| 九色成人免费人妻av| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 成人av一区二区三区在线看| 麻豆国产97在线/欧美 | 精品久久久久久久久久免费视频| 国产亚洲av高清不卡| 日本 av在线| 欧美成人性av电影在线观看| 国产精品久久电影中文字幕| 久久香蕉精品热| 久久精品aⅴ一区二区三区四区| 狠狠狠狠99中文字幕| 国产区一区二久久| 婷婷亚洲欧美| 99久久精品热视频| 午夜免费激情av| 国产高清videossex| 国产乱人伦免费视频| 一a级毛片在线观看| 亚洲av熟女| 99久久综合精品五月天人人| 日本成人三级电影网站| 天堂动漫精品| 国内精品久久久久久久电影| 国产免费男女视频| 国产免费av片在线观看野外av| 亚洲aⅴ乱码一区二区在线播放 | 视频区欧美日本亚洲| 操出白浆在线播放| 91九色精品人成在线观看| av天堂在线播放| 露出奶头的视频| 美女 人体艺术 gogo| 波多野结衣巨乳人妻| 日本免费一区二区三区高清不卡| 叶爱在线成人免费视频播放| 性欧美人与动物交配| 欧美性长视频在线观看| 国产主播在线观看一区二区| 18禁裸乳无遮挡免费网站照片| 正在播放国产对白刺激| 老司机靠b影院| 亚洲18禁久久av| 亚洲精品美女久久av网站| 香蕉国产在线看| 夜夜躁狠狠躁天天躁| 999久久久精品免费观看国产| 99精品欧美一区二区三区四区| 久久性视频一级片| 亚洲国产欧美网| 亚洲精品在线观看二区| 我要搜黄色片| 欧美成人免费av一区二区三区| 日韩国内少妇激情av| 搡老岳熟女国产| 51午夜福利影视在线观看| 人妻久久中文字幕网| 欧美激情久久久久久爽电影| 国产精品,欧美在线| 成人手机av| 久久久精品大字幕| 天天躁夜夜躁狠狠躁躁| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 老司机深夜福利视频在线观看| 日韩欧美精品v在线| 18禁美女被吸乳视频| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 欧美黄色片欧美黄色片| 国产日本99.免费观看| 搡老熟女国产l中国老女人| 国产在线观看jvid| 欧美乱妇无乱码| 免费一级毛片在线播放高清视频| 巨乳人妻的诱惑在线观看| 一个人观看的视频www高清免费观看 | 亚洲成av人片免费观看| 一本精品99久久精品77| ponron亚洲| 国产精品精品国产色婷婷| 亚洲免费av在线视频| 岛国在线观看网站| 在线永久观看黄色视频| 亚洲人与动物交配视频| 日韩精品中文字幕看吧| 国产精品一区二区三区四区久久| e午夜精品久久久久久久| 岛国在线观看网站| 日韩三级视频一区二区三区| 亚洲欧美精品综合久久99| 欧美精品亚洲一区二区| 国产蜜桃级精品一区二区三区| 欧美乱妇无乱码| 成人精品一区二区免费| 啦啦啦免费观看视频1| 色精品久久人妻99蜜桃| 人人妻人人看人人澡| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品成人免费网站| 亚洲成av人片在线播放无| 欧美av亚洲av综合av国产av| 黄片大片在线免费观看| 黑人操中国人逼视频| 亚洲成人中文字幕在线播放| 久久天躁狠狠躁夜夜2o2o| 人妻丰满熟妇av一区二区三区| 欧美+亚洲+日韩+国产| 俄罗斯特黄特色一大片| 久久午夜综合久久蜜桃| 国产亚洲精品第一综合不卡| 青草久久国产| 亚洲欧美日韩东京热| 欧美黑人精品巨大| 国内精品久久久久久久电影| 久久香蕉精品热| 国产一区二区在线av高清观看| 成年免费大片在线观看| 午夜两性在线视频| 性色av乱码一区二区三区2| 成人18禁在线播放| 午夜免费成人在线视频| 丝袜人妻中文字幕| 欧美日本视频| 美女黄网站色视频| 成年免费大片在线观看| 麻豆成人av在线观看| 免费高清视频大片| 老司机福利观看| 精品少妇一区二区三区视频日本电影| 亚洲人成77777在线视频| 91麻豆av在线| 久久久久久人人人人人| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产精品成人综合色| 91大片在线观看| 日韩精品青青久久久久久| av视频在线观看入口| 一本精品99久久精品77| 嫩草影视91久久| 免费看美女性在线毛片视频| 欧美久久黑人一区二区| 日本五十路高清| cao死你这个sao货| 中文字幕av在线有码专区| 巨乳人妻的诱惑在线观看| 禁无遮挡网站| 国产午夜精品论理片| 草草在线视频免费看| 久久午夜亚洲精品久久| 两性午夜刺激爽爽歪歪视频在线观看 | 精品无人区乱码1区二区| 免费无遮挡裸体视频| 亚洲熟妇中文字幕五十中出| 欧美中文综合在线视频| 2021天堂中文幕一二区在线观| 久久亚洲真实| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区三| 18禁黄网站禁片午夜丰满| 久久久久精品国产欧美久久久| 国产亚洲精品av在线| 五月伊人婷婷丁香| 男人的好看免费观看在线视频 | 村上凉子中文字幕在线| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区精品| 精品高清国产在线一区| 国产精品影院久久| 在线观看美女被高潮喷水网站 | 最新在线观看一区二区三区| 国产探花在线观看一区二区| 精品电影一区二区在线| 免费看美女性在线毛片视频| 91麻豆av在线| 变态另类成人亚洲欧美熟女| 亚洲精品色激情综合| 99在线视频只有这里精品首页| 啦啦啦观看免费观看视频高清| 精品高清国产在线一区| 此物有八面人人有两片| 男女之事视频高清在线观看| 中出人妻视频一区二区| 少妇熟女aⅴ在线视频| 日韩精品免费视频一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 人人妻人人澡欧美一区二区| 90打野战视频偷拍视频| 亚洲人与动物交配视频| 91老司机精品| 999久久久精品免费观看国产| 国产精品久久久人人做人人爽| 在线观看日韩欧美| 欧美高清成人免费视频www| 亚洲成人免费电影在线观看| 亚洲av五月六月丁香网| 国产三级中文精品| 亚洲真实伦在线观看| 99在线人妻在线中文字幕| 给我免费播放毛片高清在线观看| 91大片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 他把我摸到了高潮在线观看| tocl精华| www日本在线高清视频| 波多野结衣高清作品| 欧美成人免费av一区二区三区| 日韩欧美 国产精品| 大型av网站在线播放| 中文字幕久久专区| 在线免费观看的www视频| 亚洲第一欧美日韩一区二区三区| 亚洲 国产 在线| 99riav亚洲国产免费| 香蕉久久夜色| 91成年电影在线观看| 国产片内射在线| av中文乱码字幕在线| 一本久久中文字幕| 国产精品国产高清国产av| 亚洲精品粉嫩美女一区| 精品一区二区三区av网在线观看| 他把我摸到了高潮在线观看| 国产精品久久电影中文字幕| 巨乳人妻的诱惑在线观看| 可以在线观看的亚洲视频| 国产乱人伦免费视频| av在线天堂中文字幕| 国产又黄又爽又无遮挡在线| 国产av一区二区精品久久| 午夜激情福利司机影院| 久久久久久九九精品二区国产 | 日本免费一区二区三区高清不卡| 两性夫妻黄色片| 手机成人av网站| 午夜福利免费观看在线| 高潮久久久久久久久久久不卡| 亚洲国产精品sss在线观看| 欧美在线一区亚洲| 99热这里只有是精品50| 日本免费一区二区三区高清不卡| 黄色视频,在线免费观看| 可以在线观看的亚洲视频| 免费在线观看日本一区| 91字幕亚洲| 淫秽高清视频在线观看| 日本一区二区免费在线视频| 国产精品电影一区二区三区| 老熟妇仑乱视频hdxx| 亚洲熟妇中文字幕五十中出| 一级作爱视频免费观看| 91在线观看av| 少妇人妻一区二区三区视频| cao死你这个sao货| 久久婷婷人人爽人人干人人爱| 久久九九热精品免费| 黑人欧美特级aaaaaa片| 一进一出抽搐gif免费好疼| 亚洲人成77777在线视频| 美女扒开内裤让男人捅视频| 不卡一级毛片| 久久香蕉国产精品| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 日韩欧美国产在线观看| 国产亚洲av嫩草精品影院| 一区二区三区国产精品乱码| 欧美午夜高清在线| 又爽又黄无遮挡网站| 12—13女人毛片做爰片一| 国产亚洲精品久久久久5区| 亚洲国产精品久久男人天堂| 日韩中文字幕欧美一区二区| 怎么达到女性高潮| 国产爱豆传媒在线观看 | 欧美中文综合在线视频| 黑人巨大精品欧美一区二区mp4| 老汉色av国产亚洲站长工具| 国产成人aa在线观看| 又黄又粗又硬又大视频| 国产精品香港三级国产av潘金莲| 人成视频在线观看免费观看| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉精品热| videosex国产| av片东京热男人的天堂| 又黄又粗又硬又大视频| 国产在线观看jvid| 亚洲精华国产精华精| 最近最新中文字幕大全电影3| 亚洲中文字幕一区二区三区有码在线看 | 又紧又爽又黄一区二区| 国产免费av片在线观看野外av| 免费看十八禁软件| 亚洲欧美日韩高清在线视频| 久久久久久人人人人人| 中文字幕熟女人妻在线| 精品欧美一区二区三区在线| 黑人巨大精品欧美一区二区mp4| 久久久久久久精品吃奶| 成人国产综合亚洲| 男女视频在线观看网站免费 | 免费人成视频x8x8入口观看| 精品国内亚洲2022精品成人| 免费在线观看视频国产中文字幕亚洲| 欧美不卡视频在线免费观看 | 亚洲天堂国产精品一区在线| 亚洲乱码一区二区免费版| www.www免费av| 真人做人爱边吃奶动态| 日韩精品免费视频一区二区三区| 高清毛片免费观看视频网站| 亚洲最大成人中文| 床上黄色一级片| 人妻久久中文字幕网| 亚洲专区中文字幕在线| 国产日本99.免费观看| 欧美日韩精品网址| 琪琪午夜伦伦电影理论片6080| 久久午夜综合久久蜜桃| 欧美在线黄色| 日韩有码中文字幕| 三级男女做爰猛烈吃奶摸视频| 嫩草影院精品99| 精品一区二区三区av网在线观看| 亚洲片人在线观看| 夜夜夜夜夜久久久久| 麻豆久久精品国产亚洲av| 母亲3免费完整高清在线观看| 在线国产一区二区在线| 91麻豆精品激情在线观看国产| 欧美黄色淫秽网站| 三级国产精品欧美在线观看 | 男人的好看免费观看在线视频 | 九九热线精品视视频播放| 黄色丝袜av网址大全| 最好的美女福利视频网| 在线观看午夜福利视频| 亚洲专区字幕在线| 国产精品久久久人人做人人爽| 在线观看免费午夜福利视频| www日本在线高清视频| 久久精品影院6| 非洲黑人性xxxx精品又粗又长| 老熟妇乱子伦视频在线观看| 在线永久观看黄色视频| 非洲黑人性xxxx精品又粗又长| av福利片在线观看| 亚洲熟妇中文字幕五十中出| 欧美黄色淫秽网站| 黑人巨大精品欧美一区二区mp4| 99国产综合亚洲精品| 亚洲欧洲精品一区二区精品久久久| 日韩精品青青久久久久久| 国产成人影院久久av| 91九色精品人成在线观看| 亚洲精品色激情综合| 亚洲av电影不卡..在线观看| 国产精品久久久久久亚洲av鲁大| 国产视频内射| 国产一区二区在线观看日韩 | 久久精品成人免费网站| 小说图片视频综合网站| 亚洲五月婷婷丁香| 日本免费一区二区三区高清不卡| 99久久精品国产亚洲精品| 2021天堂中文幕一二区在线观| 首页视频小说图片口味搜索| 男人舔女人下体高潮全视频| 久久午夜综合久久蜜桃| 19禁男女啪啪无遮挡网站| 一区二区三区国产精品乱码| 午夜亚洲福利在线播放| 777久久人妻少妇嫩草av网站| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽| 亚洲成av人片免费观看| 成人国产一区最新在线观看| 国产单亲对白刺激| 在线观看一区二区三区| 国产av一区在线观看免费| 亚洲,欧美精品.| 怎么达到女性高潮| 最近视频中文字幕2019在线8| 欧美乱色亚洲激情| 亚洲国产欧美一区二区综合| 亚洲中文字幕日韩| 国产麻豆成人av免费视频| 人人妻人人澡欧美一区二区| 欧美另类亚洲清纯唯美| 免费在线观看完整版高清| 国产探花在线观看一区二区| 天天添夜夜摸| 少妇熟女aⅴ在线视频| 一级片免费观看大全| 夜夜看夜夜爽夜夜摸| 欧美人与性动交α欧美精品济南到| 亚洲 国产 在线| 欧美丝袜亚洲另类 | а√天堂www在线а√下载| 精品福利观看| 精品欧美国产一区二区三| 久久久久久人人人人人| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 怎么达到女性高潮| 俄罗斯特黄特色一大片| 国产精品99久久99久久久不卡| 久久精品国产亚洲av高清一级| 久久精品国产综合久久久| 国产精品一区二区三区四区久久| 欧美高清成人免费视频www| 国产黄色小视频在线观看| 麻豆成人av在线观看| 1024视频免费在线观看| 免费在线观看影片大全网站| 最新在线观看一区二区三区| 91字幕亚洲| 欧美日韩一级在线毛片| 日本黄大片高清| 在线观看午夜福利视频| 国产伦一二天堂av在线观看| 女生性感内裤真人,穿戴方法视频| 久久久久精品国产欧美久久久| av在线播放免费不卡| 两人在一起打扑克的视频| 精品国产超薄肉色丝袜足j| 狠狠狠狠99中文字幕| 欧美乱码精品一区二区三区| 国产黄a三级三级三级人| 精品第一国产精品| 色综合亚洲欧美另类图片| 欧美成人免费av一区二区三区| 亚洲自拍偷在线| 久久久久性生活片| 好男人电影高清在线观看| 男人的好看免费观看在线视频 | 亚洲av第一区精品v没综合| 久久精品国产综合久久久| 视频区欧美日本亚洲| 国产精品电影一区二区三区| 国产成年人精品一区二区| 观看免费一级毛片| 久久天堂一区二区三区四区| 国产激情久久老熟女| 国内揄拍国产精品人妻在线| 美女扒开内裤让男人捅视频| 亚洲av美国av| 人人妻人人看人人澡| 岛国在线免费视频观看| 成人特级黄色片久久久久久久| 国产成人一区二区三区免费视频网站| 中文字幕人成人乱码亚洲影| 精品久久久久久久久久免费视频| 久久中文字幕人妻熟女| 国产精品98久久久久久宅男小说| 久99久视频精品免费| 国产在线观看jvid| 亚洲人成网站高清观看| 国内毛片毛片毛片毛片毛片| 国产亚洲精品久久久久久毛片| 亚洲人成电影免费在线| 五月伊人婷婷丁香| 精品国产乱子伦一区二区三区| 丝袜美腿诱惑在线| 亚洲人成网站在线播放欧美日韩| 亚洲人与动物交配视频| 国产主播在线观看一区二区| 日日摸夜夜添夜夜添小说| 国产高清有码在线观看视频 | a级毛片在线看网站| 每晚都被弄得嗷嗷叫到高潮| 欧美中文日本在线观看视频| 丰满人妻一区二区三区视频av | 久久久久久亚洲精品国产蜜桃av| 久久精品国产清高在天天线| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放 | 欧美成人免费av一区二区三区| 亚洲男人天堂网一区| 亚洲av片天天在线观看| 1024视频免费在线观看| 欧美又色又爽又黄视频| 级片在线观看| 久久久国产精品麻豆| 亚洲av五月六月丁香网| 日本一区二区免费在线视频| 国产精品久久久久久人妻精品电影| 亚洲精品中文字幕在线视频| 国产日本99.免费观看| 久久99热这里只有精品18| 99精品在免费线老司机午夜| cao死你这个sao货| 精品电影一区二区在线| 亚洲国产精品久久男人天堂| 视频区欧美日本亚洲| 国产精品一区二区精品视频观看| 99久久国产精品久久久| 丰满的人妻完整版| 中文字幕高清在线视频| 一级片免费观看大全| 又粗又爽又猛毛片免费看| 色av中文字幕| 国产av不卡久久| 日韩国内少妇激情av| tocl精华| 亚洲国产看品久久| 一级毛片精品| 亚洲精华国产精华精| 波多野结衣高清无吗| 久久久久精品国产欧美久久久| 亚洲av电影在线进入| 在线观看免费视频日本深夜| 精品熟女少妇八av免费久了| 草草在线视频免费看| 国产成人精品久久二区二区免费| 久久精品91无色码中文字幕| 午夜免费成人在线视频| 久久人妻av系列| 久久精品影院6| 亚洲中文日韩欧美视频| 国产精品av久久久久免费| netflix在线观看网站| 亚洲七黄色美女视频| 亚洲精品一卡2卡三卡4卡5卡| 久久精品成人免费网站| 在线观看舔阴道视频| 99热只有精品国产| 怎么达到女性高潮| 88av欧美| 久久久久精品国产欧美久久久| 亚洲国产精品sss在线观看| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 国产又黄又爽又无遮挡在线| 十八禁网站免费在线| 成人精品一区二区免费| 日韩欧美免费精品| 成年版毛片免费区| 国产成人影院久久av| 12—13女人毛片做爰片一| 免费在线观看成人毛片| 日韩成人在线观看一区二区三区| 国产激情久久老熟女| 丰满人妻一区二区三区视频av | 久久精品国产亚洲av香蕉五月| 两个人看的免费小视频| 天堂动漫精品| 亚洲午夜理论影院| 亚洲一码二码三码区别大吗| 午夜久久久久精精品| 精品电影一区二区在线| 热99re8久久精品国产| 波多野结衣高清无吗| 精品日产1卡2卡| 中文字幕人成人乱码亚洲影| 久久国产精品人妻蜜桃| 亚洲中文av在线| 757午夜福利合集在线观看| 一级片免费观看大全| 麻豆成人av在线观看| 久久香蕉激情| 亚洲自拍偷在线| 成人18禁在线播放| 亚洲熟妇中文字幕五十中出| 看黄色毛片网站| 久久人妻福利社区极品人妻图片| 又黄又爽又免费观看的视频| 日本 欧美在线| 人妻夜夜爽99麻豆av| 两人在一起打扑克的视频| 亚洲美女视频黄频| 久久久精品大字幕| 久久久精品欧美日韩精品| 曰老女人黄片| 久久午夜亚洲精品久久| 亚洲真实伦在线观看| 午夜成年电影在线免费观看| 美女高潮喷水抽搐中文字幕| 黄色毛片三级朝国网站| а√天堂www在线а√下载| 欧美成人午夜精品| 热99re8久久精品国产| 一a级毛片在线观看| 国内精品久久久久精免费|