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

    基于NSGA?Ⅱ的在役狀態(tài)TMD和結(jié)構(gòu)固有參數(shù)識別方法研究

    2025-04-10 00:00:00趙文韜溫青孫洪鑫王修勇華旭剛
    振動工程學(xué)報 2025年3期
    關(guān)鍵詞:結(jié)構(gòu)

    摘要: 調(diào)諧質(zhì)量阻尼器(tuned mass damper,TMD)的減振效率與結(jié)構(gòu)和TMD的固有參數(shù)相關(guān),準(zhǔn)確從結(jié)構(gòu)?TMD耦合系統(tǒng)響應(yīng)中識別結(jié)構(gòu)和TMD固有參數(shù)是對在役TMD減振性能評價的必要條件。提出了一種基于NSGA?Ⅱ(非支配排序遺傳算法,nondominated sorting genetic algorithm)的參數(shù)識別方法,從結(jié)構(gòu)?TMD耦合運動響應(yīng)中識別“裸結(jié)構(gòu)”和“裸TMD”的固有參數(shù),進(jìn)而實現(xiàn)對在役狀態(tài)TMD的減振性能評估。該方法構(gòu)建了結(jié)構(gòu)?TMD耦合運動方程,并將其減縮和轉(zhuǎn)化為結(jié)構(gòu)被控模態(tài)和TMD耦合的兩自由度系統(tǒng),借助系統(tǒng)狀態(tài)空間矩陣搭建兩個目標(biāo)函數(shù),通過遺傳算法尋找理論值與試驗值的最小誤差所對應(yīng)的最優(yōu)解,從而識別結(jié)構(gòu)和TMD的固有參數(shù)。開展了單自由度結(jié)構(gòu)?TMD耦合系統(tǒng)和多自由度結(jié)構(gòu)?TMD耦合系統(tǒng)參數(shù)識別數(shù)值仿真分析,結(jié)果表明:提出的方法可以從耦合系統(tǒng)動力響應(yīng)中準(zhǔn)確識別結(jié)構(gòu)和TMD的固有參數(shù)。

    關(guān)鍵詞: 參數(shù)識別; 結(jié)構(gòu)?TMD耦合系統(tǒng); NSGA?Ⅱ

    中圖分類號: TU311.3; TU352.1 " "文獻(xiàn)標(biāo)志碼: A " "文章編號: 1004-4523(2025)03-0507-10

    DOI:10.16385/j.cnki.issn.1004-4523.2025.03.007

    Identification method of in?service TMD and structural inherent parameters based on NSGA?Ⅱ

    ZHAO Wentao1, WEN Qing1, SUN Hongxin1, WANG Xiuyong1, HUA Xugang2

    (1.Hunan Provincial Key Laboratory of Structures for Wind Resistance and Vibration Control,

    Hunan University of Science and Technology, Xiangtan 411201, China;

    2.Key Laboratory for Wind and Bridge Engineering of Hunan Province, Hunan University, Changsha 410082, China)

    Abstract: The vibration reduction efficiency of a tuned mass damper (TMD) is closely related to the inherent parameters of the structure and TMD. Accurate identification of the structure and TMD inherent parameters from the response of the structure-TMD coupling system is necessary for the evaluation of the vibration reduction performance of in-service TMD. This paper offers a parameter identification approach based on NSGA-Ⅱ(Nondominated sorting genetic algorithm) that can identify the parameters of ‘bare structure’ and ‘bare TMD’ from the coupled structure-TMD response in order to solve the issue of state evaluation of TMDs in service. The structure-TMD coupling equation is constructed. It is reduced and transformed into a two-degree-of-freedom system of the structure-controlled mode and TMD coupling. Two objective functions are constructed by means of the system state space matrix. The genetic algorithm is used to find the optimal solution corresponding to the minimum error between the theoretical value and the experimental value, so as to identify the modal parameters of the structure and TMD. The numerical simulation analysis of parameter identification of single-degree-of-freedom structure-TMD coupling system and multi-degree-of-freedom structure-TMD coupling system is carried out. The results show that the proposed method can accurately identify the inherent parameters of structure and TMD from the dynamic system response of the coupling system.

    Keywords: parameter identification;structure-TMD coupling system;NSGA-Ⅱ

    調(diào)諧質(zhì)量阻尼器(tuned mass damper,TMD)對于結(jié)構(gòu)的人致振動和風(fēng)致振動都有良好的抑制效果[1?8],已被廣泛應(yīng)用于工程實際中。TMD主要由質(zhì)量塊、彈簧元件和阻尼元件三部分組成,其減振效率與質(zhì)量比、頻率、阻尼比等參數(shù)相關(guān)[9]。在實際工程中,TMD的參數(shù)會因為設(shè)計誤差、安裝誤差以及長期使用造成的疲勞損傷而偏離最優(yōu)設(shè)計參數(shù),影響其減振性能[10]。廣大學(xué)者對其減振性能評估問題展開了深入探討,現(xiàn)有的做法是通過對比鎖定前后結(jié)構(gòu)的動力響應(yīng)特性,進(jìn)而間接評估其減振性能。這種做法耗時耗力,且無法準(zhǔn)確評估其固有動力參數(shù),無法為現(xiàn)場精準(zhǔn)調(diào)諧提供理論指導(dǎo)。因此,探索從結(jié)構(gòu)?TMD耦合系統(tǒng)動力響應(yīng)中識別TMD和被控結(jié)構(gòu)固有參數(shù)的方法,具有重要的工程應(yīng)用價值,既可以避免頻繁鎖定TMD,也能為減振性能評估提供理論基礎(chǔ)。

    在役狀態(tài)TMD減振性能評估問題的難點在于從結(jié)構(gòu)?TMD耦合系統(tǒng)中評估“裸結(jié)構(gòu)”和“裸TMD”的固有動力參數(shù)特性。目前只有少數(shù)研究涉及此類參數(shù)估計方法,HWANG等[11]提出了一種基于H∞最優(yōu)模型降階和系統(tǒng)識別的參數(shù)識別方法,可以評估“裸結(jié)構(gòu)”的模態(tài)質(zhì)量;HAZRA等[12]提出了一種時域算法估計“裸結(jié)構(gòu)”的模態(tài)特性,并可以通過此算法對TMD進(jìn)行二次調(diào)諧;KANG等[13]采用系統(tǒng)辨識和逆模態(tài)風(fēng)荷載估計的方法探究TMD在臺風(fēng)作用下的動力特性,并對其減振性能進(jìn)行了評價;溫青等[14?15]提出了一種基于隨機子空間系統(tǒng)識別的模態(tài)參數(shù)識別方法,可以評估結(jié)構(gòu)和TMD的模態(tài)參數(shù),并通過單層框架試驗加以驗證。

    這些針對在役狀態(tài)TMD的性能評估研究存在一定的局限性,識別參數(shù)不全,無法實現(xiàn)TMD的現(xiàn)場調(diào)諧,或識別方法受主觀判斷影響較大。鑒于此,本文從信號處理和數(shù)據(jù)擬合的角度出發(fā),提出了一種基于NSGA?Ⅱ的結(jié)構(gòu)和TMD固有參數(shù)識別方法,首先,構(gòu)建多自由度結(jié)構(gòu)?TMD耦合運動方程,將其減縮和轉(zhuǎn)化為結(jié)構(gòu)被控模態(tài)和TMD耦合的兩自由度系統(tǒng)連續(xù)狀態(tài)空間模型;然后,構(gòu)建狀態(tài)空間方程與耦合運動響應(yīng)以及耦合頻率的函數(shù)關(guān)系,進(jìn)而建立兩個目標(biāo)函數(shù),并采用NSGA?Ⅱ算法尋找理論計算值與試驗值之間相差最小時所對應(yīng)的最優(yōu)解集,從而得出參數(shù)識別結(jié)果;最后采用數(shù)值仿真驗證該方法的可行性,并探究了噪聲以及多模態(tài)耦合對于該方法的影響程度。該方法理論簡單、識別誤差可以由目標(biāo)函數(shù)值直觀反映,且參數(shù)識別結(jié)果受人為判斷影響較小。

    1 參數(shù)識別方法

    1.1 結(jié)構(gòu)?TMD耦合系統(tǒng)狀態(tài)矩陣

    假定一個N自由度結(jié)構(gòu),采用單個TMD控制其第r階模態(tài)振動,為了保證良好的減振效果,將其放置在結(jié)構(gòu)第r階模態(tài)振型的最大位移處,由此組成具有N+1個自由度的結(jié)構(gòu)?TMD耦合系統(tǒng),其力學(xué)簡圖如圖1所示。為了方便計算,結(jié)構(gòu)阻尼采用Rayleigh阻尼,則上述系統(tǒng)的自由運動方程可以表示為:

    1.3 基于遺傳算法的參數(shù)優(yōu)化

    由圖2可知,本文提出的參數(shù)識別方法的關(guān)鍵工作之一是多參數(shù)目標(biāo)函數(shù)優(yōu)化。對于這種優(yōu)化參數(shù)較多且目標(biāo)函數(shù)較為復(fù)雜的模型,難以獲得目標(biāo)函數(shù)的導(dǎo)數(shù),不適應(yīng)于梯度優(yōu)化方法,傳統(tǒng)的最小二乘法受初值選取的影響較大,容易陷入局部最優(yōu),故本文選用具有全局搜索能力的遺傳算法進(jìn)行參數(shù)識別。

    遺傳算法是一種基于生物自然選擇與遺傳機理的概率搜索方法[16?17]。隨著遺傳算法的發(fā)展,遺傳算法也應(yīng)用于多目標(biāo)優(yōu)化領(lǐng)域,有眾多學(xué)者提出了基于Pareto最優(yōu)概念的多目標(biāo)遺傳算法,例如NSGA[18]、NSGA?Ⅱ[19?20]、NSGA?Ⅲ[21]等,其中,NSGA?Ⅱ通過引入快速非支配排序算法、精英策略、擁擠度比較算子,降低了算法的計算復(fù)雜度,保證了種群的多樣性,且具有較好的算法魯棒性,故本文采用NSGA?Ⅱ多目標(biāo)遺傳算法尋找最優(yōu)解,流程圖如圖3所示,其主要步驟如下:

    圖3 遺傳算法流程圖

    Fig.3 Flow chart of genetic algorithm

    (1)設(shè)定種群大小以及5個待識別參數(shù)的上下限范圍,隨機生成初始種群;其中ωr和ωd的參數(shù)上下限可以參考實測耦合模態(tài)頻率ω1和ω2,根據(jù)TMD減振理論,ωr和ωd在數(shù)值大小上總是處于ω1和ω2之間,且ω1和ω2之間的數(shù)值差異性會隨著結(jié)構(gòu)和TMD的阻尼比增大而減小。通常TMD的阻尼比較大,所以上述設(shè)定可以有效縮小ωr和ωd的搜索范圍,極大提高計算效率。結(jié)構(gòu)阻尼比ζr通常較小,可以設(shè)定在0~2%,TMD阻尼比ζd可以設(shè)定在2%~15%,質(zhì)量比μ可以設(shè)定在0~5%,上述參數(shù)設(shè)定范圍也可以根據(jù)尋優(yōu)結(jié)果逐步縮小搜索范圍。

    (2)計算每個種群個體所對應(yīng)的兩個目標(biāo)函數(shù)值Qopt1和Qopt2,依據(jù)目標(biāo)函數(shù)值對種群個體進(jìn)行快速非支配排序,并對每個非支配層的個體進(jìn)行擁擠度計算。

    (3)支配等級代表了種群個體的優(yōu)劣性,目標(biāo)函數(shù)值越小,支配等級越高,所對應(yīng)的次序號越??;而擁擠度值則反映了種群的分布性和多樣性,擁擠度值越高,說明種群個體的聚集密度越小。選擇初代種群中優(yōu)秀且聚集密度較小的個體執(zhí)行交叉和變異操作,產(chǎn)生新的一代種群,并與初代種群合并,再次進(jìn)行快速非支配排序和擁擠度計算,得到新一代種群。

    (4)設(shè)定最大進(jìn)化代數(shù)作為收斂準(zhǔn)則,滿足收斂條件則得到尋優(yōu)結(jié)果,否則繼續(xù)進(jìn)行循環(huán)。由于是多目標(biāo)尋優(yōu),而兩個目標(biāo)函數(shù)值難以同時達(dá)到最小,即沒有最優(yōu)解,只能得到兩個目標(biāo)函數(shù)相對均衡的Pareto解集,其中,Qopt1表明了理論響應(yīng)與實際響應(yīng)之間的差異性,差異性會隨著迭代次數(shù)增加而不斷被放大,而Qopt2僅表示理論耦合模態(tài)頻率與實際耦合模態(tài)頻率的差異性,因此,二者在數(shù)值上可能存在數(shù)量級的差異性。由于Qopt2數(shù)值極小,最優(yōu)個體之間以此為參照,意義較小。因此,本文選取Qopt1為最優(yōu)個體的評判標(biāo)準(zhǔn),且通過量化評判標(biāo)準(zhǔn)也可以減小人為判斷對尋優(yōu)結(jié)果的影響。

    2 數(shù)值仿真

    為了驗證本文提出的參數(shù)識別方法的可行性和有效性,分別以單自由度結(jié)構(gòu)?TMD耦合系統(tǒng)以及多自由度結(jié)構(gòu)?TMD耦合系統(tǒng)為研究對象,首先給定耦合系統(tǒng)的相關(guān)參數(shù),并計算其在自由振動情況下的動力響應(yīng),以此作為觀測數(shù)據(jù),然后采用上述參數(shù)識別方法從觀測數(shù)據(jù)中識別系統(tǒng)模態(tài)參數(shù),通過與設(shè)定值對比,驗證參數(shù)識別方法的可行性,同時探究了噪聲以及多模態(tài)耦合對于參數(shù)識別結(jié)果的影響。

    2.1 單自由度結(jié)構(gòu)?TMD耦合系統(tǒng)

    建立單自由度結(jié)構(gòu)?TMD耦合系統(tǒng),其力學(xué)簡圖如圖4所示。假定結(jié)構(gòu)的質(zhì)量m1、剛度k1和阻尼系數(shù)c1分別為50 kg、7895.70 N/m和6.28 N·s/m。由此可以得出結(jié)構(gòu)的頻率f1與結(jié)構(gòu)阻尼比ζ1分別為2 Hz和0.5%,按照Den Hartog提出的經(jīng)典TMD優(yōu)化理論確定TMD的設(shè)計參數(shù),假定質(zhì)量比為μ,則f0=f1/(1+μ),ζ_0=√(3μ/(8(1+μ))),然后根據(jù)質(zhì)量比以及TMD的頻率f0和阻尼比ζ0可以得到TMD的質(zhì)量m0、剛度k0以及阻尼系數(shù)c0。考慮到數(shù)值仿真的完備性,依據(jù)結(jié)構(gòu)與TMD是否調(diào)諧以及不同質(zhì)量比等情況共設(shè)定9組工況,具體參數(shù)如表1所示。在仿真分析時,設(shè)定采樣頻率為100 Hz,結(jié)構(gòu)和TMD的初始位移分別為0.01 m和0,初始速度均為0。

    以工況Ⅷ為例,依據(jù)設(shè)定的系統(tǒng)參數(shù),采用龍格?庫塔法計算結(jié)構(gòu)和TMD在自由振動下的動力響應(yīng)數(shù)值解,動力響應(yīng)及其功率譜圖如圖5所示。由功率譜圖中的峰值點對應(yīng)的橫坐標(biāo)可知,結(jié)構(gòu)?TMD耦合系統(tǒng)的兩階模態(tài)頻率實測值分別為1.8298和2.1098 Hz。以上述動力響應(yīng)數(shù)值解和模態(tài)頻率實測值作為觀測值,采用第1節(jié)參數(shù)識別方法識別結(jié)構(gòu)和TMD的模態(tài)參數(shù)。在仿真過程中,設(shè)定遺傳算法的種群大小為75,最優(yōu)個體系數(shù)為0.2,交叉率為0.8,變異率為0.05,最大進(jìn)化代數(shù)為300,模態(tài)參數(shù)(ω1、ω0、ζ1、ζ0、μ)的上下限分別為(1.8298×2π、18298×2π、0、0.03、0.005)和(2.1098×2π、2.1098×2π、0.02、0.15、0.05)。

    由于設(shè)定的種群大小為75,最優(yōu)個體系數(shù)為0.2,因此形成了由15個最優(yōu)個體組成的Pareto解集,如圖6所示。圖6(a)~(e)給出了工況Ⅷ中每個最優(yōu)個體的參數(shù)識別值及其誤差棒圖,可知:(1)種群中的15個最優(yōu)個體之間數(shù)值差異性較小,參數(shù)識別結(jié)果穩(wěn)定,既表明了種群中優(yōu)秀個體基因得以保留也表明了算法的穩(wěn)定性;(2)參數(shù)識別結(jié)果與真實值近乎重合,表明方法可以有效并準(zhǔn)確地識別結(jié)構(gòu)和TMD的模態(tài)參數(shù);(3)圖6(c)中結(jié)構(gòu)阻尼比的識別誤差棒較為明顯,由式(8)可知,結(jié)構(gòu)阻尼比ζr在系統(tǒng)矩陣A中只有1項,且識別精度與ωr、ωd、μ、ζd相關(guān),因此,這些參數(shù)的識別誤差都會傳遞給ζr,造成其識別誤差相對較大。圖6(f)中給出了每個最優(yōu)個體所對應(yīng)的Qopt1(目標(biāo)函數(shù)1)和Qopt2(目標(biāo)函數(shù)2)的值,由圖中的橫、縱坐標(biāo)軸可以直觀了解到參數(shù)識別的精度較高,且Qopt2在數(shù)值上遠(yuǎn)遠(yuǎn)小于Qopt1,所以選取Qopt1作為評判標(biāo)準(zhǔn)更為合理。

    在無噪聲情況下,以上9個工況的參數(shù)識別結(jié)果如表2所示??梢缘贸?,該方法在不同質(zhì)量比以及TMD頻率失諧或阻尼比失諧的情況下,都具有較高的參數(shù)識別精度。且TMD和結(jié)構(gòu)的頻率識別誤差小于1%,質(zhì)量比和TMD阻尼比識別誤差小于3%,結(jié)構(gòu)的阻尼比識別誤差小于6%,完全滿足TMD的現(xiàn)場調(diào)試需求,上述結(jié)果表明該方法理論正確、切實可行。

    在實際工程中,噪聲的影響不可避免,包括由外界環(huán)境激勵引起的激勵噪聲以及采集數(shù)據(jù)時儀器自身產(chǎn)生的測量噪聲。在仿真過程中,以零均值的高斯白噪聲模擬激勵噪聲R1和測量噪聲R2,將其疊加于結(jié)構(gòu)和TMD的響應(yīng)中,噪聲的強度由信噪比SNR確定,SNR=101g(Ps/Pn),其中Ps和Pn分別表示信號的能量和噪聲的能量,Ps可以由原信號與其平均值的差值平方和確定,Pn可以由原信號與帶噪聲的信號的差值平方和確定,信噪比越高,噪聲越小。為了探究噪聲對參數(shù)識別精度的影響,以工況Ⅷ為例,通過對比不同信噪比情況下的參數(shù)識別結(jié)果,分析測量噪聲和激勵噪聲對于參數(shù)識別精度的影響并評估該參數(shù)識別方法的抗干擾能力。表3給出了有噪聲情況下的參數(shù)識別結(jié)果,識別結(jié)果表明,隨著信噪比的減小,Qopt1無明顯變化趨勢,Qopt2有顯著的增大趨勢,主要因為噪聲的干擾增大了理論值與實際響應(yīng)的誤差,但即便在有較大噪聲干擾的情況下,該方法依然能準(zhǔn)確識別結(jié)構(gòu)和TMD的模態(tài)參數(shù),在信噪比大于30的范圍內(nèi),固有參數(shù)識別的誤差沒有因噪聲的影響而增大,表明提出的參數(shù)識別方法具有較好的抗噪能力。

    2.2 多自由度結(jié)構(gòu)?TMD耦合系統(tǒng)

    為了模擬該方法在實際工程結(jié)構(gòu)中應(yīng)用的可行性,以如圖7所示的3自由度結(jié)構(gòu)?TMD耦合系統(tǒng)為例開展了仿真分析。結(jié)構(gòu)的物理參數(shù)如下:m1=m2=m3=50 kg,ki=3000 N/m,ci=6 N·s/m,i=1,2,3,4。為了控制結(jié)構(gòu)的一階模態(tài)振動,將TMD安裝在結(jié)構(gòu)的1階振型的最大位移處,即質(zhì)量塊m2上,設(shè)定TMD的質(zhì)量比為2%,即m0為2 kg,TMD的頻率和阻尼比分別為0.92 Hz和8.57%,由此組成3自由度結(jié)構(gòu)?TMD耦合系統(tǒng),表4給出了原結(jié)構(gòu)以及附加TMD耦合系統(tǒng)的模態(tài)參數(shù),可知,TMD對于原結(jié)構(gòu)的第1階模態(tài)影響較大,且產(chǎn)生了較大的附加模態(tài)阻尼比,對原結(jié)構(gòu)的第2、3階模態(tài)影響較小。

    在理論推導(dǎo)過程中,式(5)的成立條件為振動必須只有TMD和被控模態(tài)參與,顯然在多自由度情況下不完全成立,多自由度結(jié)構(gòu)的動力響應(yīng)是由多模態(tài)疊加而成,因此會受到其他模態(tài)的干擾,為了探究該參數(shù)識別方法在多模態(tài)耦合情況下的適用性,設(shè)定4種工況:(1)TMD的初始位移x0=1 cm,質(zhì)量塊m1、m2、m3的初始位移x1、x2、x3皆為0;(2)TMD的初始位移x0=1 cm,質(zhì)量塊m1的初始位移x1=1 cm,其他質(zhì)量塊的初始位移為0;(3)TMD的初始位移x0=1 cm,質(zhì)量塊m2的初始位移x2=1 cm,其他質(zhì)量塊的初始位移為0;(4)TMD的初始位移x0=1 cm,質(zhì)量塊m3的初始位移x3=1 cm,其他質(zhì)量塊的初始位移為0。

    圖8給出了工況(2)~(4)結(jié)構(gòu)的位移響應(yīng)以及其功率譜圖,從圖中可以看出每個工況都至少有兩階模態(tài)參與其中。在不同的多模態(tài)參與下,各工況的識別結(jié)果如表5所示,從表中可以看出,多模態(tài)振動對參數(shù)識別結(jié)果影響較小,且TMD和結(jié)構(gòu)的頻率識別誤差小于1%,質(zhì)量比和TMD阻尼比識別誤差小于3%,TMD阻尼比識別誤差小于6%,結(jié)果表明在多自由度結(jié)構(gòu)多模態(tài)干擾的情況下,該方法依然適用。

    2.3 多自由度結(jié)構(gòu)?TMD有限元仿真

    為了探究上述參數(shù)識別方法在實際工程結(jié)構(gòu)中的適用性,以廈門某人行天橋為工程背景,建立其ANSYS有限元模型并進(jìn)行數(shù)值模擬分析。該橋主跨長90 m,橋面寬4.845 m,主跨采用張弦桁架結(jié)構(gòu)體系,上弦桿采用鋼桁架,下部設(shè)置雙柔索,中部與外部斜撐桿連接,形成整體應(yīng)力自平衡體系。主梁采用鋼桁架形式,由上弦桿、下弦桿、腹桿和中間斜桿組成,橋面鋪裝采用鋼格柵。圖9為橋梁有限元模型示意圖。模擬過程中,采用Beam188單元模擬鋼桁架桿件和斜撐,采用Link8單元模擬纜索,欄桿等附屬結(jié)構(gòu)的質(zhì)量通過Mass21單元施加到主梁上,橋面鋪裝采用Shell63單元模擬。由于直接采用橋梁的有限元模型計算時,自由度較多,工作量很大,為了減少計算自由度,引入了模態(tài)綜合技術(shù),首先通過模態(tài)分析求解橋梁的模態(tài)頻率和振型,然后利用振型的正交性將耦聯(lián)的方程解耦,得到互相獨立的模態(tài)方程,且由于橋梁的結(jié)構(gòu)振動響應(yīng)以低階為主,所以只需要考慮較少的振型就可以獲得較高的計算精度,極大地減小了計算量。本文提取橋梁的前20階振型,構(gòu)建橋梁的廣義坐標(biāo)方程,并采用3種不同參數(shù)的TMD分別安裝在橋梁前3階豎彎模態(tài)的振型最大位移處,進(jìn)而控制橋梁前3階豎彎模態(tài)振動,橋梁的前3階豎彎被控模態(tài)參數(shù)以及TMD的參數(shù)如表6所示。

    類似于在實橋試驗中激振TMD帶動橋梁結(jié)構(gòu)運動,在仿真過程中,可以分別給定TMD初始位移,通過數(shù)值仿真獲取結(jié)構(gòu)?TMD耦合動態(tài)響應(yīng),圖10給出了在不同模態(tài)階次下的結(jié)構(gòu)響應(yīng)及其功率譜圖。為了探究被控模態(tài)對于參數(shù)識別方法的影響,對不同模態(tài)階次下的結(jié)構(gòu)?TMD耦合系統(tǒng)分別做了參數(shù)識別,參數(shù)識別結(jié)果如表7所示,由表可知,對1階豎彎模態(tài)的參數(shù)識別精度明顯優(yōu)于對2、3階豎彎模態(tài)的識別精度,主要原因是由于耦合程度的差異性,從圖10可以看出,1階模態(tài)僅與第3階模態(tài)耦合,耦合程度較小,而2、3階模態(tài)的耦合程度較高,從而造成識別精度的差異性。但總體識別精度依然較高,對于TMD和結(jié)構(gòu)的頻率、質(zhì)量比最大參數(shù)識別誤差小于2%,對于TMD和結(jié)構(gòu)的阻尼比最大參數(shù)識別誤差小于5%。且從圖11可以看出,Pareto解集中個體差異性較小,有較好的算法穩(wěn)定性,由此反映出在多模態(tài)參與的情況下,該方法依然能準(zhǔn)確識別結(jié)構(gòu)和TMD的固有參數(shù)。

    3 結(jié) 論

    本文提出了一種基于NSGA?Ⅱ的在役狀態(tài)結(jié)構(gòu)和TMD固有參數(shù)識別方法,通過數(shù)值仿真分析了該方法的可行性,研究結(jié)果表明:

    (1)該方法可以有效識別結(jié)構(gòu)和TMD的固有參數(shù),固有頻率識別誤差小于1%,質(zhì)量比識別誤差小于3%,阻尼比的識別誤差不超過6%;

    (2)該方法具有良好的抗噪和抗干擾能力,在信噪比大于30的范圍內(nèi),參數(shù)識別精度不會因噪聲的存在而下降。

    論文通過數(shù)值仿真分析研究了一種基于NSGA?Ⅱ的在役狀態(tài)結(jié)構(gòu)和TMD固有參數(shù)識別方法的可行性,后續(xù)將開展相關(guān)的試驗和工程應(yīng)用研究。研究成果對在役狀態(tài)TMD的性能評價和減振效率提升具有重要意義。

    參考文獻(xiàn):

    [1] 單伽锃, 呂西林. 高層結(jié)構(gòu)調(diào)諧質(zhì)量阻尼器動力測試與評價研究進(jìn)展及工程實踐[J]. 建筑結(jié)構(gòu)學(xué)報, 2023, 44(3): 1-10.

    SHAN Jiazeng, LYU Xilin. Dynamic testing and evaluation of mass damper technology for high-rise buildings with state-of-the-art practice[J]. Journal of Building Structures, 2023, 44(3): 1-10.

    [2] 王浩, 鄒仲欽, 茅建校, 等. 基于壓重塊型TMD的大跨度斜拉橋減震控制[J]. 振動工程學(xué)報, 2018, 31(3): 391-398.

    WANG Hao, ZOU Zhongqin, MAO Jianxiao, et al. Seismic control of a long-span cable-stayed bridge using TMD-type counterweight[J]. Journal of Vibration Engineering, 2018, 31(3): 391-398.

    [3] 楊維青, 華旭剛, 溫青, 等. 杠桿式調(diào)諧質(zhì)量阻尼器對雙吊索尾流致振的減振優(yōu)化[J]. 振動工程學(xué)報, 2021, 34(4): 819-827.

    YANG Weiqing, HUA Xugang, WEN Qing, et al. Optimal design of lever-type tuned mass dampers for control of wakeinduced vibrations of coupled twin-cable hangers[J]. Journal of Vibration Engineering, 2021, 34(4): 819-827.

    [4] 汪志昊, 陳銀, 胡明祎, 等. 機器擾力作用下某廠房樓板豎向振動與TMD減振研究[J]. 振動工程學(xué)報, 2019, 32(6): 986-995.

    WANG Zhihao, CHEN Yin, HU Mingyi, et al. Vertical vibration and TMD mitigation of an industrial building floor subjected to machine excitation[J]. Journal of Vibration Engineering, 2019, 32(6): 986-995.

    [5] 華旭剛, 溫青, 陳政清, 等. 大跨度雙層曲線斜拉橋人致振動減振優(yōu)化與實測驗證[J]. 振動工程學(xué)報, 2016, 29(5): 822-830.

    HUA Xugang, WEN Qing, CHEN Zhengqing, et al. Design and experimental validation of structural vibration control of a curved twin-deck cable-stayed bridge subject to pedestrians[J]. Journal of Vibration Engineering, 2016, 29(5): 822-830.

    [6] 周晅毅, 林勇建, 顧明. 大跨屋蓋結(jié)構(gòu)MTMD風(fēng)振控制最優(yōu)性能研究[J]. 振動工程學(xué)報, 2015, 28(2): 277-284.

    ZHOU Xuanyi, LIN Yongjian, GU Ming. Optimum properties of MTMD for large-span roof subjected to wind loads[J]. Journal of Vibration Engineering, 2015, 28(2): 277-284.

    [7] LUO Y F, SUN H X, WANG X Y, et al. Wind induced vibration control and energy harvesting of electromagnetic resonant shunt tuned mass-damper-inerter for building structures[J]. Shock and Vibration, 2017, 2017(1): 4180134.

    [8] 練繼建, 趙悅, 練沖, 等. 基于電渦流-調(diào)諧質(zhì)量阻尼器的海上風(fēng)電筒型基礎(chǔ)結(jié)構(gòu)減振研究[J]. 振動與沖擊, 2019, 38(19): 20-25.

    LIAN Jijian, ZHAO Yue, LIAN Chong, et al. Vibration reduction of offshore wind turbine tube infrastructures based on EC-TMD[J]. Journal of Vibration and Shock, 2019, 38(19): 20-25.

    [9] ASAMI T, WAKASONO T, KAMEOKA K, et al. Optimum design of dynamic absorbers for a system subject to random excitation[J]. Transactions of the Japan Society of Mechanical Engineers Series C, 1990, 56(523): 619-627.

    [10] 王文熙, 華旭剛, 王修勇, 等. TMD系統(tǒng)在自身參數(shù)隨機偏離下的減振有效性和可靠性分析[J]. 振動與沖擊, 2016, 35(1): 228-234.

    WANG Wenxi, HUA Xugang, WANG Xiuyong, et al. Vibration reduction validity and reliability of a TMD system under random deviation of its own parameters[J]. Journal of Vibration and Shock, 2016, 35(1): 228-234.

    [11] HWANG J S, KIM H, KIM J. Estimation of the modal mass of a structure with a tuned-mass damper using H-infinity optimal model reduction[J]. Engineering Structures, 2006, 28(1): 34-42.

    [12] HAZRA B, SADHU A, LOURENCO R, et al. Re-tuning tuned mass dampers using ambient vibration measurements[J]. Smart Materials and Structures, 2010, 19(11): 115002.

    [13] KANG N, KIM H, SEONGWOO JO S C, et al. Performance evaluation of TMD under typhoon using system identification and inverse wind load estimation[J]. Computer-Aided Civil and Infrastructure Engineering, 2012, 27(6): 455-473.

    [14] 溫青, 華旭剛, 王修勇, 等. 基于耦合系統(tǒng)環(huán)境振動試驗的結(jié)構(gòu)和TMD參數(shù)識別[J]. 振動與沖擊, 2019, 38(10): 193-198.

    WEN Qing, HUA Xugang, WANG Xiuyong, et al. Parameter identification of structure and TMD based on ambient vibration tests of the coupled system[J]. Journal of Vibration and Shock, 2019, 38(10): 193-198.

    [15] 溫青, 華旭剛, 王修勇, 等. 利用耦合系統(tǒng)自由振動響應(yīng)識別結(jié)構(gòu)被控模態(tài)和TMD參數(shù)[J]. 振動工程學(xué)報, 2019, 32(4): 565-573.

    WEN Qing, HUA Xugang, WANG Xiuyong, et al. Parameter identification of controlled structural mode and TMD based on the free vibration tests of the coupled system[J]. Journal of Vibration Engineering, 2019, 32(4): 565-573.

    [16] 何勇, 王紅, 谷穗. 一種基于遺傳算法的VMD參數(shù)優(yōu)化軸承故障診斷新方法[J]. 振動與沖擊, 2021, 40(6): 184-189.

    HE Yong, WANG Hong, GU Sui. New fault diagnosis approach for bearings based on parameter optimized VMD and genetic algorithm[J]. Journal of Vibration and Shock, 2021, 40(6): 184-189.

    [17] 莊亮東, 楊悅, 吳楨灝. 遺傳算法在Y型偏心支撐組合框架抗震性能優(yōu)化中的應(yīng)用研究[J]. 工程力學(xué), 2023, 40(7): 185-195.

    ZHUANG Liangdong, YANG Yue, WU Zhenhao. The application of genetic algorithm in seismic performance optimization of Y-shape eccentrically braced composite frame[J]. Engineering Mechanics, 2023, 40(7): 185-195.

    [18] SRINIVAS N, DEB K. Muiltiobjective optimization using nondominated sorting in genetic algorithms[J]. Evolutionary Computation, 1994, 2(3): 221-248.

    [19] RAHIMBAKHSH H, KOHANSAL M E, TARKASHVAND A, et al. Multi-objective optimization of natural surveillance and privacy in early design stages utilizing NSGA?Ⅱ[J]. Automation in Construction, 2022, 143: 104547.

    [20] DEB K, AGRAWAL S, PRATAP A, et al. A fast elitist non-dominated sorting genetic algorithm for multi-objective optimization: NSGA?Ⅱ[C]// Parallel Problem Solving from Nature PPSN Ⅵ. Berlin, Heidelberg: Springer, 2000: 849-858.

    [21] HAN L, SHI X J, ZHAI Y Y. Test optimization selection method based on NSGA-3 and improved Bayesian network model[J]. Journal of Northwestern Polytechnical University, 2021, 39(2): 414-422.

    第一作者: 趙文韜(1994―),男,博士研究生。

    E-mail:zwt_hunankejidaxue@163.com

    通信作者: 溫 "青(1985―),男,博士,副教授。

    E-mail:cewenq@hnust.edu.cn

    猜你喜歡
    結(jié)構(gòu)
    DNA結(jié)構(gòu)的發(fā)現(xiàn)
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    循環(huán)結(jié)構(gòu)謹(jǐn)防“死循環(huán)”
    論《日出》的結(jié)構(gòu)
    縱向結(jié)構(gòu)
    縱向結(jié)構(gòu)
    我國社會結(jié)構(gòu)的重建
    人間(2015年21期)2015-03-11 15:23:21
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長
    亚洲精品日韩av片在线观看 | 最近最新中文字幕大全电影3| 性色av乱码一区二区三区2| 国产av在哪里看| 国产亚洲精品综合一区在线观看| 欧美中文日本在线观看视频| 欧美黑人巨大hd| 宅男免费午夜| 91字幕亚洲| 在线十欧美十亚洲十日本专区| 欧美另类亚洲清纯唯美| 熟女人妻精品中文字幕| 色精品久久人妻99蜜桃| АⅤ资源中文在线天堂| 精品久久久久久,| 国内精品美女久久久久久| 中文字幕av在线有码专区| 色综合欧美亚洲国产小说| 在线观看午夜福利视频| 国内揄拍国产精品人妻在线| 老汉色∧v一级毛片| 亚洲无线观看免费| 高清在线国产一区| 亚洲成a人片在线一区二区| 免费人成视频x8x8入口观看| 五月玫瑰六月丁香| 五月伊人婷婷丁香| 熟女电影av网| 看黄色毛片网站| 日本三级黄在线观看| 精品熟女少妇八av免费久了| 一边摸一边抽搐一进一小说| 欧美性猛交╳xxx乱大交人| 国产91精品成人一区二区三区| 久久精品人妻少妇| 国产精品国产高清国产av| 亚洲欧美日韩高清在线视频| 老司机福利观看| 别揉我奶头~嗯~啊~动态视频| 午夜亚洲福利在线播放| 午夜免费观看网址| 午夜福利18| 国产午夜福利久久久久久| 国产精品99久久久久久久久| 蜜桃久久精品国产亚洲av| 非洲黑人性xxxx精品又粗又长| 男女视频在线观看网站免费| 波多野结衣高清作品| 亚洲欧美日韩无卡精品| 一进一出抽搐gif免费好疼| 成人高潮视频无遮挡免费网站| 日日夜夜操网爽| 精品99又大又爽又粗少妇毛片 | 搡老熟女国产l中国老女人| 国产又黄又爽又无遮挡在线| 成年免费大片在线观看| 亚洲五月婷婷丁香| 超碰av人人做人人爽久久 | 在线看三级毛片| 午夜精品久久久久久毛片777| 叶爱在线成人免费视频播放| 亚洲18禁久久av| 亚洲18禁久久av| 久久亚洲精品不卡| 日本免费一区二区三区高清不卡| 不卡一级毛片| av在线蜜桃| 九九在线视频观看精品| 两性午夜刺激爽爽歪歪视频在线观看| 午夜精品一区二区三区免费看| 99久久成人亚洲精品观看| 嫩草影院入口| 极品教师在线免费播放| 无遮挡黄片免费观看| 午夜a级毛片| 在线观看一区二区三区| 草草在线视频免费看| 国产精品av视频在线免费观看| 国产三级在线视频| 99视频精品全部免费 在线| 久久久久免费精品人妻一区二区| 久久久久免费精品人妻一区二区| 又黄又粗又硬又大视频| 窝窝影院91人妻| 国产成人av激情在线播放| 色尼玛亚洲综合影院| 国产久久久一区二区三区| 色视频www国产| 一进一出抽搐动态| 午夜精品久久久久久毛片777| 日本一二三区视频观看| 十八禁人妻一区二区| 91在线精品国自产拍蜜月 | 日本三级黄在线观看| 99久久综合精品五月天人人| 在线看三级毛片| 日韩亚洲欧美综合| 久久九九热精品免费| 国内精品久久久久久久电影| 精品一区二区三区视频在线 | 亚洲色图av天堂| 国产精品一区二区三区四区免费观看 | 99在线人妻在线中文字幕| 亚洲无线在线观看| 成年女人永久免费观看视频| 五月伊人婷婷丁香| 国内揄拍国产精品人妻在线| 一进一出好大好爽视频| 亚洲片人在线观看| 最好的美女福利视频网| 成人国产一区最新在线观看| 欧美日韩综合久久久久久 | 欧美日韩综合久久久久久 | 成人三级黄色视频| 亚洲精品一区av在线观看| 精品国产超薄肉色丝袜足j| 一区二区三区激情视频| 国产一区二区三区视频了| 每晚都被弄得嗷嗷叫到高潮| 嫩草影视91久久| 五月伊人婷婷丁香| 久久久久久久久大av| 99久久精品热视频| 国产黄a三级三级三级人| 琪琪午夜伦伦电影理论片6080| 一个人观看的视频www高清免费观看| 久99久视频精品免费| 久久久久精品国产欧美久久久| a级一级毛片免费在线观看| 国产精品1区2区在线观看.| 伊人久久大香线蕉亚洲五| 美女黄网站色视频| 黄色视频,在线免费观看| 国产亚洲av嫩草精品影院| 麻豆成人午夜福利视频| 日本三级黄在线观看| 9191精品国产免费久久| 免费在线观看影片大全网站| 国内揄拍国产精品人妻在线| 舔av片在线| 国产高清有码在线观看视频| 免费搜索国产男女视频| 黄片大片在线免费观看| 久久久久九九精品影院| 亚洲五月天丁香| 亚洲五月天丁香| 国产精品 国内视频| 欧美成人a在线观看| 男女视频在线观看网站免费| 中亚洲国语对白在线视频| 老鸭窝网址在线观看| 少妇丰满av| 欧美日韩黄片免| 悠悠久久av| 在线观看av片永久免费下载| 日韩欧美国产一区二区入口| av中文乱码字幕在线| 亚洲av成人不卡在线观看播放网| 日韩免费av在线播放| 欧美+亚洲+日韩+国产| 久久久久亚洲av毛片大全| 久久精品人妻少妇| 欧美日本亚洲视频在线播放| 国产午夜精品久久久久久一区二区三区 | 给我免费播放毛片高清在线观看| 最近在线观看免费完整版| 最近最新免费中文字幕在线| 久久亚洲真实| 18美女黄网站色大片免费观看| 毛片女人毛片| 美女 人体艺术 gogo| 欧美黄色淫秽网站| 狂野欧美激情性xxxx| 欧美另类亚洲清纯唯美| 国产精品 欧美亚洲| 久久人妻av系列| 怎么达到女性高潮| 精品久久久久久久人妻蜜臀av| 99久久99久久久精品蜜桃| 国产精品亚洲美女久久久| 日本熟妇午夜| 国语自产精品视频在线第100页| 三级男女做爰猛烈吃奶摸视频| 国产精品亚洲一级av第二区| 亚洲中文字幕一区二区三区有码在线看| 国产极品精品免费视频能看的| ponron亚洲| 亚洲美女黄片视频| av国产免费在线观看| 国产精品久久久久久精品电影| 亚洲精品色激情综合| 波多野结衣高清作品| 90打野战视频偷拍视频| 51午夜福利影视在线观看| 亚洲五月天丁香| 国产亚洲精品久久久com| 国产成人aa在线观看| 国产精品亚洲av一区麻豆| 欧美国产日韩亚洲一区| 精品久久久久久久久久免费视频| 人人妻人人看人人澡| 国产免费男女视频| 99视频精品全部免费 在线| 老熟妇仑乱视频hdxx| 亚洲真实伦在线观看| 午夜福利欧美成人| 国产精品,欧美在线| 国产成人福利小说| 国产一区二区三区在线臀色熟女| 亚洲欧美一区二区三区黑人| 18美女黄网站色大片免费观看| 一进一出抽搐动态| 精品午夜福利视频在线观看一区| 亚洲国产中文字幕在线视频| 国产色爽女视频免费观看| 看黄色毛片网站| 欧美日韩精品网址| 国产高清视频在线播放一区| 日韩欧美精品v在线| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 99国产精品一区二区蜜桃av| 免费观看的影片在线观看| 午夜激情欧美在线| 国产成人欧美在线观看| xxx96com| 在线播放无遮挡| 国产精品98久久久久久宅男小说| 日韩欧美 国产精品| 无遮挡黄片免费观看| 久久中文看片网| 国产精品女同一区二区软件 | 日韩欧美在线乱码| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| av欧美777| 老汉色av国产亚洲站长工具| 国产高清有码在线观看视频| 在线观看日韩欧美| 91麻豆av在线| 啦啦啦韩国在线观看视频| 熟妇人妻久久中文字幕3abv| 国产极品精品免费视频能看的| 人人妻人人澡欧美一区二区| 欧美午夜高清在线| 51午夜福利影视在线观看| 国产精品嫩草影院av在线观看 | 国产伦精品一区二区三区四那| 黄色日韩在线| 欧美日本亚洲视频在线播放| 99热只有精品国产| 日日干狠狠操夜夜爽| 99在线视频只有这里精品首页| 久久精品国产亚洲av涩爱 | 制服丝袜大香蕉在线| 麻豆成人午夜福利视频| 色av中文字幕| 18+在线观看网站| 欧美日韩黄片免| 亚洲美女视频黄频| 岛国在线免费视频观看| 51国产日韩欧美| 国产在视频线在精品| 天天一区二区日本电影三级| 好男人电影高清在线观看| 国产精品免费一区二区三区在线| 日日干狠狠操夜夜爽| svipshipincom国产片| 精品不卡国产一区二区三区| 91av网一区二区| 19禁男女啪啪无遮挡网站| 一个人观看的视频www高清免费观看| 黄色片一级片一级黄色片| 国内久久婷婷六月综合欲色啪| 国产精品亚洲美女久久久| 日韩免费av在线播放| 深夜精品福利| 国产黄色小视频在线观看| 熟女人妻精品中文字幕| 精品人妻一区二区三区麻豆 | 熟女少妇亚洲综合色aaa.| 中文资源天堂在线| 国产乱人视频| 久久久久久久久大av| 成人午夜高清在线视频| 欧美黄色淫秽网站| 国产在线精品亚洲第一网站| 欧美不卡视频在线免费观看| 日韩欧美精品v在线| 丰满乱子伦码专区| 男人和女人高潮做爰伦理| 最近最新中文字幕大全电影3| 白带黄色成豆腐渣| 欧美日本亚洲视频在线播放| 久久精品影院6| 亚洲五月天丁香| 首页视频小说图片口味搜索| 啦啦啦韩国在线观看视频| 国产视频一区二区在线看| 亚洲精品日韩av片在线观看 | 性色avwww在线观看| 国产69精品久久久久777片| 国产色婷婷99| 精品不卡国产一区二区三区| www日本黄色视频网| 国产精品精品国产色婷婷| 在线观看免费午夜福利视频| 最近最新中文字幕大全电影3| 手机成人av网站| 久久国产乱子伦精品免费另类| 成年女人永久免费观看视频| 男插女下体视频免费在线播放| 人妻久久中文字幕网| 特级一级黄色大片| 亚洲精品一区av在线观看| 国产成人福利小说| 男插女下体视频免费在线播放| 亚洲第一电影网av| 黄色日韩在线| 最新中文字幕久久久久| 超碰av人人做人人爽久久 | 成人精品一区二区免费| 亚洲电影在线观看av| 草草在线视频免费看| 夜夜爽天天搞| 午夜福利成人在线免费观看| 亚洲无线观看免费| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 真人一进一出gif抽搐免费| 亚洲不卡免费看| 欧美日韩精品网址| www日本黄色视频网| 国内毛片毛片毛片毛片毛片| 一进一出抽搐动态| a级一级毛片免费在线观看| av在线蜜桃| 亚洲精品美女久久久久99蜜臀| 男插女下体视频免费在线播放| 尤物成人国产欧美一区二区三区| 岛国在线观看网站| 一级a爱片免费观看的视频| 乱人视频在线观看| 很黄的视频免费| 麻豆一二三区av精品| 日韩免费av在线播放| 19禁男女啪啪无遮挡网站| 久9热在线精品视频| 亚洲在线自拍视频| 国产老妇女一区| 久久婷婷人人爽人人干人人爱| 丝袜美腿在线中文| 精品99又大又爽又粗少妇毛片 | 午夜亚洲福利在线播放| 亚洲avbb在线观看| 白带黄色成豆腐渣| 草草在线视频免费看| 国产精品久久久久久人妻精品电影| 性欧美人与动物交配| 免费看光身美女| 美女免费视频网站| 欧美激情在线99| 亚洲av电影不卡..在线观看| 天天一区二区日本电影三级| 97人妻精品一区二区三区麻豆| 欧美在线一区亚洲| 久久伊人香网站| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 亚洲欧美日韩东京热| 不卡一级毛片| 久久亚洲真实| 国产激情欧美一区二区| 在线免费观看的www视频| svipshipincom国产片| 少妇人妻一区二区三区视频| xxxwww97欧美| 亚洲精品亚洲一区二区| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 国产黄色小视频在线观看| 免费av观看视频| 久久香蕉国产精品| 波多野结衣巨乳人妻| 老汉色av国产亚洲站长工具| 免费人成在线观看视频色| 国产精品久久久久久精品电影| 欧美中文综合在线视频| 国产精品久久久久久精品电影| 99久久成人亚洲精品观看| 免费无遮挡裸体视频| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 欧美zozozo另类| 国产亚洲精品av在线| 淫秽高清视频在线观看| 一本综合久久免费| 最近最新中文字幕大全免费视频| 精品福利观看| 18禁裸乳无遮挡免费网站照片| 欧美区成人在线视频| 国产成人影院久久av| 波野结衣二区三区在线 | 国产色爽女视频免费观看| 黄片小视频在线播放| 男女那种视频在线观看| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 十八禁网站免费在线| www.熟女人妻精品国产| 亚洲精品在线美女| 亚洲国产欧美人成| 久久精品综合一区二区三区| 国产主播在线观看一区二区| 亚洲精华国产精华精| 在线观看免费午夜福利视频| a级一级毛片免费在线观看| 日韩成人在线观看一区二区三区| 性色avwww在线观看| 国产探花极品一区二区| 亚洲内射少妇av| 蜜桃亚洲精品一区二区三区| 免费看光身美女| 美女黄网站色视频| 好男人在线观看高清免费视频| 小蜜桃在线观看免费完整版高清| 人人妻,人人澡人人爽秒播| 亚洲av第一区精品v没综合| 国产av一区在线观看免费| 久久久久久人人人人人| 国产黄色小视频在线观看| 国产高清videossex| av中文乱码字幕在线| 嫩草影院入口| 亚洲av成人精品一区久久| 久久久久亚洲av毛片大全| 久久精品91无色码中文字幕| 中文字幕久久专区| 国产真人三级小视频在线观看| 日本一二三区视频观看| 免费人成在线观看视频色| 啦啦啦免费观看视频1| 最近在线观看免费完整版| 网址你懂的国产日韩在线| 国产精品三级大全| 一级毛片女人18水好多| 亚洲美女视频黄频| 婷婷精品国产亚洲av在线| 国产男靠女视频免费网站| 欧美高清成人免费视频www| 欧美激情在线99| 国产91精品成人一区二区三区| 岛国在线观看网站| 757午夜福利合集在线观看| 国产伦在线观看视频一区| 好看av亚洲va欧美ⅴa在| 一级毛片女人18水好多| 99久久精品一区二区三区| 国产成人av激情在线播放| 最近最新中文字幕大全电影3| 一级黄片播放器| 美女被艹到高潮喷水动态| 午夜福利18| 免费观看的影片在线观看| www日本在线高清视频| 日本精品一区二区三区蜜桃| 天美传媒精品一区二区| 脱女人内裤的视频| 少妇丰满av| 高清毛片免费观看视频网站| 中文资源天堂在线| 国产真人三级小视频在线观看| 日韩欧美免费精品| 高潮久久久久久久久久久不卡| 欧美bdsm另类| 亚洲人与动物交配视频| 黄色成人免费大全| 香蕉av资源在线| 中出人妻视频一区二区| 好男人电影高清在线观看| 91在线精品国自产拍蜜月 | 国产野战对白在线观看| 天天躁日日操中文字幕| 国产成人a区在线观看| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 成年女人看的毛片在线观看| 免费搜索国产男女视频| 成人一区二区视频在线观看| 搡女人真爽免费视频火全软件 | 操出白浆在线播放| 亚洲av二区三区四区| 国产三级黄色录像| 人妻久久中文字幕网| 亚洲专区中文字幕在线| 日韩欧美精品免费久久 | 成人欧美大片| 五月伊人婷婷丁香| 国产精华一区二区三区| 亚洲av成人精品一区久久| 亚洲真实伦在线观看| 桃红色精品国产亚洲av| www日本在线高清视频| a级毛片a级免费在线| 九色国产91popny在线| www日本黄色视频网| 99热6这里只有精品| 国产99白浆流出| 久久天躁狠狠躁夜夜2o2o| 国产伦在线观看视频一区| 亚洲无线在线观看| 男女那种视频在线观看| 国产精品香港三级国产av潘金莲| 欧美色欧美亚洲另类二区| 国内精品美女久久久久久| 国产精品久久视频播放| 天天一区二区日本电影三级| 女生性感内裤真人,穿戴方法视频| 精品国内亚洲2022精品成人| 少妇裸体淫交视频免费看高清| 欧美最新免费一区二区三区 | e午夜精品久久久久久久| 国内揄拍国产精品人妻在线| 啦啦啦观看免费观看视频高清| 99国产综合亚洲精品| 精品福利观看| 在线观看日韩欧美| 2021天堂中文幕一二区在线观| 亚洲精品乱码久久久v下载方式 | 搞女人的毛片| 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 非洲黑人性xxxx精品又粗又长| avwww免费| 淫妇啪啪啪对白视频| 久久久精品欧美日韩精品| 国产精品国产高清国产av| 国产亚洲av嫩草精品影院| 欧美又色又爽又黄视频| 亚洲久久久久久中文字幕| 男女床上黄色一级片免费看| 在线播放无遮挡| 老鸭窝网址在线观看| 国产免费一级a男人的天堂| 国产69精品久久久久777片| 日本 av在线| 看免费av毛片| 搡老妇女老女人老熟妇| 久久精品91蜜桃| 欧美丝袜亚洲另类 | 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 十八禁网站免费在线| 国产午夜精品论理片| 欧美xxxx黑人xx丫x性爽| 日韩精品青青久久久久久| 欧美日韩亚洲国产一区二区在线观看| 高清日韩中文字幕在线| 人人妻人人看人人澡| 国产精品久久久久久人妻精品电影| 国产极品精品免费视频能看的| 最近在线观看免费完整版| 三级男女做爰猛烈吃奶摸视频| 观看美女的网站| 一级毛片女人18水好多| 国产蜜桃级精品一区二区三区| 美女 人体艺术 gogo| 国产精品1区2区在线观看.| 床上黄色一级片| 久久人妻av系列| 国产欧美日韩一区二区三| 精品99又大又爽又粗少妇毛片 | 国产在视频线在精品| 国产高清激情床上av| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 免费看日本二区| 少妇人妻一区二区三区视频| 国产亚洲欧美98| 91av网一区二区| 亚洲久久久久久中文字幕| 亚洲片人在线观看| 日韩欧美在线二视频| www日本黄色视频网| 天堂网av新在线| 国产精品久久久久久人妻精品电影| 91字幕亚洲| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久久毛片微露脸| 国产精品久久久久久久电影 | 亚洲,欧美精品.| 久久精品亚洲精品国产色婷小说| 1024手机看黄色片| 少妇丰满av| 亚洲成av人片在线播放无| 给我免费播放毛片高清在线观看| 91av网一区二区| 真人做人爱边吃奶动态| 久久久久九九精品影院| 亚洲国产精品sss在线观看| 人妻丰满熟妇av一区二区三区| 女人十人毛片免费观看3o分钟| 在线观看av片永久免费下载| 国产精品永久免费网站| 听说在线观看完整版免费高清| 日本黄色视频三级网站网址| 少妇的丰满在线观看| 国产一区二区亚洲精品在线观看| 亚洲精品一卡2卡三卡4卡5卡| aaaaa片日本免费| 美女高潮喷水抽搐中文字幕| 亚洲av五月六月丁香网| x7x7x7水蜜桃| 岛国视频午夜一区免费看| 老司机在亚洲福利影院|