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

    高聳煙囪結(jié)構(gòu)調(diào)諧質(zhì)量慣容阻尼器(TMDI)風(fēng)振控制方法及效果研究

    2022-11-05 10:27:40彭士濤洪寧寧
    工程力學(xué) 2022年11期
    關(guān)鍵詞:風(fēng)振最優(yōu)控制煙囪

    蘇 寧,彭士濤,洪寧寧

    (交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,天津 300456)

    高聳煙囪結(jié)構(gòu)對(duì)風(fēng)荷載作用極為敏感,尤其是橫風(fēng)向渦激共振會(huì)對(duì)結(jié)構(gòu)安全造成不利影響,其風(fēng)振控制是結(jié)構(gòu)抗風(fēng)優(yōu)化設(shè)計(jì)中的重點(diǎn)和難點(diǎn)。調(diào)諧質(zhì)量慣容阻尼器(tuned mass damper inerter,TMDI)是一種基于慣容器的新型的動(dòng)力吸振裝置。一種常見(jiàn)的齒輪齒條式慣容器示意圖如圖1 所示,慣容器提供的阻尼力fb與兩個(gè)連接端相對(duì)加速度u¨1-u¨2成正比。比例系數(shù)b稱為慣容系數(shù),其取值僅與齒輪和飛輪的機(jī)械參數(shù)有關(guān),通過(guò)齒輪和飛輪的合理選型,慣容器能夠自身重量成百上千倍的質(zhì)量慣性效應(yīng)[1],充分利用該效應(yīng)有助于實(shí)現(xiàn)動(dòng)力吸振器的輕量化設(shè)計(jì)[2-3],近年來(lái)得到廣泛關(guān)注。

    DEN HARTOG 基于不動(dòng)點(diǎn)理論[4],給出了傳統(tǒng)的調(diào)諧質(zhì)量阻尼器(TMD)的參數(shù)優(yōu)化及振動(dòng)控制性能經(jīng)驗(yàn)公式。后期學(xué)者們就主結(jié)構(gòu)阻尼比、激勵(lì)荷載形式等方面對(duì)上述理論進(jìn)行了完善[5-6],形成了完整系統(tǒng)的結(jié)構(gòu)振動(dòng)控制理論體系。另有學(xué)者們,在此基礎(chǔ)上結(jié)合等效線性化方法,提出了調(diào)諧液柱阻尼器(TLCD)[7-8]、非線性能量阱(NES)[9-11]等非線性吸振器的分析體系。

    TMDI 引入慣容器后,使得運(yùn)動(dòng)方程更為復(fù)雜,學(xué)者們對(duì)相應(yīng)的參數(shù)優(yōu)化及控制效率評(píng)估方法進(jìn)行了討論[12-13],但這些研究都是基于慣容器與地面相連進(jìn)行推導(dǎo)的,又稱為GTMDI(grounded TMDI)。而對(duì)于高層或高聳結(jié)構(gòu),慣容器接地難以在工程中實(shí)現(xiàn),一些學(xué)者將GTMDI 與基底隔震系統(tǒng)相結(jié)合[14],使其更具實(shí)際應(yīng)用價(jià)值。在應(yīng)用層面,TMDI 系統(tǒng)更容易安裝在高層結(jié)構(gòu)的層間[15-16],或在連體高層的振動(dòng)控制中發(fā)揮作用[17-18]。此外,還應(yīng)用于大跨度橋梁結(jié)構(gòu)的斷面間[19],以控制渦激共振的不利影響。

    研究表明,慣容器的連接位置不僅影響其在工程中的實(shí)施難度,而且對(duì)TMDI 的控制效果影響較大,且僅在一定范圍內(nèi),效果明顯優(yōu)于TMD[20],但文中尚未針對(duì)慣容器的位置的影響進(jìn)行參數(shù)化分析。后有學(xué)者建立了TMDI 控制廣義單自由度振動(dòng)的懸臂柱振動(dòng)響應(yīng)的運(yùn)動(dòng)方程,討論了截面變化對(duì)控制效果的影響[21],說(shuō)明TMDI 的控制效果與慣容器連接位置處的振型關(guān)系密切,但文中未給出相關(guān)的經(jīng)驗(yàn)公式供設(shè)計(jì)參考。

    本文將高聳煙囪簡(jiǎn)化為廣義單自由度結(jié)構(gòu),基于風(fēng)荷載頻譜的濾波表示,推導(dǎo)了TMDI 控制下結(jié)構(gòu)風(fēng)振響應(yīng)的解析解。在此基礎(chǔ)上,對(duì)TMDI最優(yōu)設(shè)計(jì)參數(shù)進(jìn)行了參數(shù)化分析,總結(jié)了相應(yīng)的經(jīng)驗(yàn)公式供設(shè)計(jì)參考。此外,通過(guò)對(duì)比TMDI 與TMD的風(fēng)振響應(yīng)控制效果,給出了慣容器起增強(qiáng)控制效果的判別準(zhǔn)則,以及TMDI 的等效TMD 質(zhì)量比計(jì)算公式,以指導(dǎo)動(dòng)力吸振器的輕量化設(shè)計(jì)。最后通過(guò)對(duì)某270 m 高混凝土煙囪風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行TMDI 風(fēng)振控制算例分析,驗(yàn)證了理論分析的有效性,并探討了TMDI 對(duì)高聳煙囪的風(fēng)振響應(yīng)控制效果。

    假設(shè)圓柱形煙囪內(nèi)、外半徑沿高度z的變化分別為Ri(z)和Ro(z),將高聳煙囪結(jié)構(gòu)分析模型假設(shè)為連續(xù)質(zhì)量的振動(dòng)體系,質(zhì)量、剛度沿高度分布m(z)、EI(z)分別為:

    1 理論推導(dǎo)

    1.1 運(yùn)動(dòng)方程

    式中:ρs、E分別為煙囪結(jié)構(gòu)的材料密度和彈性模量;me(z)為恒活荷載等效質(zhì)量代表值沿高度的分布。

    將無(wú)控的煙囪結(jié)構(gòu)(不含TMDI)進(jìn)一步簡(jiǎn)化為廣義單自由度振動(dòng)體系(圖2)。沿高度變化的振動(dòng)位移u(z,t)解耦表示為廣義位移響應(yīng)x(t)與歸一化振型函數(shù)Φ(z)的乘積:

    其中,Φ(H) = 1,則x(t)在數(shù)值上等同于頂點(diǎn)位移u(H,t)。結(jié)構(gòu)的廣義質(zhì)量M、廣義剛度K分別表示為:

    當(dāng)采用TMDI 進(jìn)行振動(dòng)控制時(shí),TMDI 的質(zhì)量為m,剛度為k,阻尼為c,慣容系數(shù)為b,慣容器連接位置(高度)χ。慣容器的連接位置(高度)χ與振動(dòng)控制效果具有很強(qiáng)關(guān)聯(lián)。一般來(lái)說(shuō),χ越小,慣容器的兩端距離越大,兩端相對(duì)加速度越大,振動(dòng)控制效果越好。同時(shí),χ的減小會(huì)增加慣容器的安裝難度。因此,χ的取值需要權(quán)衡控制效果與可實(shí)施性。在理論推導(dǎo)中,定義慣容器位置參數(shù)為其連接位置處的振型值,即φ= Φ(χ),這里,φ隨χ的單調(diào)增減。

    TMDI 的振動(dòng)頻率ωd和阻尼比ζd分別定義為:

    定義TMDI 調(diào)制比參數(shù),質(zhì)量比μ=m/M,慣質(zhì)比β =b/M,頻率比ν=ωd/ωs。根據(jù)虛功原理[22],可將體系運(yùn)動(dòng)方程組表示為:

    1.2 振動(dòng)響應(yīng)分析

    其中,Hx(ω)為煙囪主結(jié)構(gòu)位移響應(yīng)頻響函數(shù)的分子多項(xiàng)式,該多項(xiàng)式為偶次多項(xiàng)式,其偶次項(xiàng)系數(shù)為θj(j= 0, 1, 2)。

    TMDI 的響應(yīng)控制比J0定義為控制后的響應(yīng)標(biāo)準(zhǔn)差與無(wú)控響應(yīng)標(biāo)準(zhǔn)差之比,響應(yīng)控制率η 定義為控制后響應(yīng)X減小量占無(wú)控響應(yīng)X0的百分比,當(dāng)以脈動(dòng)響應(yīng)標(biāo)準(zhǔn)差(均方根)計(jì)算響應(yīng)控制率時(shí),η = 1 -J0。

    由式(14)、式(21)的頻域積分可得結(jié)構(gòu)在無(wú)控和控制狀態(tài)下的響應(yīng),當(dāng)假設(shè)風(fēng)荷載激勵(lì)假設(shè)為強(qiáng)度為δ 的白噪聲時(shí),即Sξ(ω) = δ,根據(jù)閉路積分定理[22]可得響應(yīng)的解析表達(dá)式如下:

    式中,N0和Nx為控制響應(yīng)的分子和分母表達(dá)式。則白噪聲激勵(lì)下響應(yīng)控制比為:

    將廣義風(fēng)荷載頻譜表示為如式(31)所示的白噪聲濾波的形式,則可求得式(14)、式(21)的解析表達(dá)式,能夠更為便捷高效地分析各參數(shù)對(duì)響應(yīng)及控制效果的影響[23]。

    1.3 風(fēng)荷載譜的濾波表示

    式中: Λa(?)、δa分別為順風(fēng)向風(fēng)荷載譜Sξa(ω)的濾波多項(xiàng)式和歸一化強(qiáng)度;ωa表示順風(fēng)向無(wú)量綱風(fēng)荷載譜峰值所對(duì)應(yīng)的圓頻率。圖3(a)給出了橫風(fēng)向無(wú)量綱廣義風(fēng)荷載譜風(fēng)洞試驗(yàn)數(shù)據(jù)和相關(guān)模型的對(duì)比,可見(jiàn),采用該濾波模型總體能夠較好地反映順風(fēng)向風(fēng)荷載譜在高頻段的衰減趨勢(shì)。

    對(duì)于橫風(fēng)向風(fēng)荷載,由于漩渦脫落引起的渦激振動(dòng),使得能量在斯托羅哈頻率附近集中,無(wú)量綱風(fēng)荷載譜曲線具有較為顯著的凸起,其濾波多項(xiàng)式可由二次函數(shù)來(lái)表示:

    根據(jù)閉路積分定理[22]可得,無(wú)控狀態(tài)下的順風(fēng)向及橫風(fēng)向風(fēng)振響應(yīng)為:

    值得說(shuō)明的是,采用式(35)計(jì)算橫風(fēng)向響應(yīng)時(shí),阻尼比ζT= ζs+ ζa為考慮氣動(dòng)阻尼ζa的總阻尼比。

    橫風(fēng)向風(fēng)振響應(yīng)σxa按式(40)~式(43)計(jì)算。式中,Dc0和Dc1為橫風(fēng)向響應(yīng)的分母和分子行列式,分別由式(41)、式(42)計(jì)算;κcj為橫風(fēng)向響應(yīng)濾波多項(xiàng)式Kc(?)=Λc(?)·Γ(?)的j次項(xiàng)系數(shù)(j=0, 1, 2, 3, 4, 5, 6),由式(41)確定。其中, γ′j為式(20)中的γj(j= 1, 2, 3),但結(jié)構(gòu)阻尼比ζs由總阻尼比ζT=ζs+ ζa替代。

    順風(fēng)向及橫風(fēng)向風(fēng)振響應(yīng)的控制比不僅隨阻尼器參數(shù)變化,還隨著風(fēng)荷載頻率與結(jié)構(gòu)頻率之比(ωma/ωs、ωmc/ωs)變化,一般地,當(dāng)風(fēng)荷載頻率與結(jié)構(gòu)頻率接近產(chǎn)生共振時(shí),響應(yīng)控制效果最好。定義順風(fēng)向、橫風(fēng)向的共振響應(yīng)控制比Ja、Jc來(lái)描述TMDI 對(duì)順風(fēng)向及橫風(fēng)向響應(yīng)的控制效果。

    圖4 給出了幾個(gè)典型工況下白噪聲激勵(lì)下的響應(yīng)控制比Jw(由式(30)計(jì)算得到)、順風(fēng)向及橫風(fēng)向共振控制比Ja、Jc(分別由式(44)、式(45)計(jì)算得到)隨TMDI 頻率比ν和阻尼比ζd變化的等值線圖。由圖4 可以看出,Jw、Ja、Jc在(ν,ζd)平面上存在一個(gè)極小值點(diǎn),即理論最優(yōu)參數(shù)。

    2 TMDI 參數(shù)優(yōu)化

    2.1 參數(shù)尋優(yōu)方法

    TMDI 的參數(shù)優(yōu)化指根據(jù)質(zhì)量、慣容參數(shù)(μ,β,φ),確 定 最 優(yōu) 的 頻 率 比 和 阻 尼 比(νopt,ζdopt)使得響應(yīng)控制效果最優(yōu),也就是,在特定的參數(shù)(μ,β,φ)下,尋找J0(或Jw、Ja、Jc)曲面的最小值點(diǎn)。第1 節(jié)給出了J0的解析計(jì)算方法,可大幅度提高計(jì)算效率和精度。由圖3 分析可知,給定參數(shù)(μ,β,φ)時(shí),一般存在唯一的理論最優(yōu)參數(shù)(νopt,ζdopt)可由下式得到:

    式中:J0也可以用Jw、Ja、Jc替代;Jopt為最優(yōu)參數(shù)(νopt,ζdopt)下的響應(yīng)控制比。

    理論上可以利用式(30)得到的白噪聲激勵(lì)下的Jw解析表達(dá)式推得最優(yōu)參數(shù)(νopt,ζdopt)隨(μ,β,φ)變化的表達(dá)式。但由于導(dǎo)函數(shù)表達(dá)形式過(guò)于復(fù)雜,無(wú)法得到解析解。因此,本文采用數(shù)值優(yōu)化方法—共軛梯度法,在不同的(μ,β,φ)條件下對(duì)(νopt,ζdopt)進(jìn)行尋優(yōu)。利用基于Den Hartog的不動(dòng)點(diǎn)理論得到的GTMDI 最優(yōu)解(式(47))[16]尋優(yōu)迭代的初值進(jìn)行計(jì)算。

    研究表明,結(jié)構(gòu)阻尼比ζs對(duì)最優(yōu)參數(shù)(νopt,ζdopt)結(jié)果影響較小[10],因此在參數(shù)分析時(shí),僅考慮ζs= 1%,分析參數(shù)(μ,β,φ)的變化范圍為:

    1)μ= 10^[-3∶0.1∶-1] (共41 種取值);

    2) β = [0, 0.01, 0.02, 0.05, 0.1: 0.1: 1.0] (共14 種取值);

    3)φ= [0: 0.1: 0.5, 0.60: 0.05: 0.80, 0.82: 0.02:1.00] (共21 種取值)。

    對(duì)白噪聲、順風(fēng)向、橫風(fēng)向3 種激勵(lì)模式下的控制比Jw、Ja、Jc進(jìn)行分析,共計(jì)36 162 個(gè)計(jì)算工況。

    圖5 給出了μ= 1%情況下,三種激勵(lì)模式最優(yōu)控制參數(shù)(νopt,ζdopt)和最優(yōu)控制比Jopt隨慣容器參數(shù)(β,φ)的變化??傮w看來(lái),隨著慣質(zhì)比β 增大、慣容器連接位置降低(φ減小)最優(yōu)頻率比νopt降低,最優(yōu)阻尼比ζdopt增大,這種變化趨勢(shì)隨著φ增大而減弱。慣容器連接位置降低(φ減小),最優(yōu)控制比Jopt減小,控制效果增強(qiáng)。當(dāng)φ=0 時(shí),GTMDI 最優(yōu)控制比隨慣質(zhì)比β 增大而減小。而當(dāng)φ= 1 時(shí),隨著慣質(zhì)比β 增大,Jopt呈現(xiàn)增大趨勢(shì)。雖然三種激勵(lì)模式下,雖然最優(yōu)控制參數(shù)(νopt,ζdopt)存在一些差異,但Jopt對(duì)激勵(lì)模式不敏感。

    2.2 最優(yōu)參數(shù)經(jīng)驗(yàn)公式

    根據(jù)上述參數(shù)分析結(jié)果,總結(jié)出TMDI 最優(yōu)控制參數(shù)(νopt,ζdopt)關(guān)于質(zhì)量、慣容參數(shù)(μ,β,φ)的經(jīng)驗(yàn)公式如下:

    圖5(a)、圖5(b)給出了最優(yōu)控制參數(shù)(νopt,ζdopt)經(jīng)驗(yàn)公式和分析結(jié)果的對(duì)比,可以發(fā)現(xiàn),經(jīng)驗(yàn)公式和分析結(jié)果趨勢(shì)一致,但局部存在一定差異。將采用經(jīng)驗(yàn)公式得到的最優(yōu)控制參數(shù)下的最優(yōu)控制比與分析結(jié)果的對(duì)比,如圖6 所示??梢园l(fā)現(xiàn),即便經(jīng)驗(yàn)公式估計(jì)最優(yōu)控制參數(shù)造成了一定誤差,但對(duì)最優(yōu)控制比結(jié)果影響不大。這可能是由于最優(yōu)控制參數(shù)(νopt, ζdopt)附近的控制比曲面J0(ν, ζd)形狀較為平坦,從一定程度上削弱了最優(yōu)控制比Jopt結(jié)果對(duì)最優(yōu)控制參數(shù)(νopt, ζdopt)誤差的敏感性。對(duì)Jopt的誤差進(jìn)一步統(tǒng)計(jì)分析發(fā)現(xiàn),97%以上的Jopt結(jié)果誤差不超過(guò)10%,說(shuō)明該經(jīng)驗(yàn)公式在的誤差是工程中可接受的,可以用于估算TMDI 的最優(yōu)參數(shù)。optdoptopt

    特別是,當(dāng)φ= 0,式(48)變化為GTMDI 優(yōu)化參數(shù)的取值,見(jiàn)式(49);進(jìn)一步,當(dāng)β = 0 時(shí),式(48)變化為T(mén)MD 最優(yōu)參數(shù)取值,見(jiàn)式(50)。

    3 最優(yōu)風(fēng)振控制效果分析

    3.1 與TMD 最優(yōu)控制效果的比較

    TMDI 研制的初衷之一是通過(guò)慣容器增強(qiáng)TMD 的振動(dòng)控制效果。通過(guò)上述研究發(fā)現(xiàn),慣容器的連接位置對(duì)TMDI 的振動(dòng)控制效果影響至關(guān)重要。由圖5(c)可知,慣容器接地時(shí)(φ= 0,GTMDI)的控制效果最優(yōu)。而對(duì)比圖4(a)和圖4(d)發(fā)現(xiàn),當(dāng)φ= 1 時(shí),其控制效果可能不如相同質(zhì)量比下的TMD。為比較相同質(zhì)量比μ下TMDI 與TMD 的控制效果,定義慣容器影響系數(shù)R為T(mén)MDI 與TMD的響應(yīng)控制比之比,即:

    當(dāng)R> 1 時(shí),TMDI 的控制效果不如TMD;當(dāng)R≤ 1 時(shí),慣容器起到增強(qiáng)振動(dòng)控制的效果,稱為慣容增強(qiáng)區(qū)。R值越小,增強(qiáng)效果越顯著。

    圖7 給出了2.1 節(jié)參數(shù)分析中μ= 1%時(shí)R隨慣容器參數(shù)(β,φ)的變化等值線圖??梢园l(fā)現(xiàn),慣容增強(qiáng)區(qū)內(nèi)(R≤ 1),隨著慣質(zhì)比β 增大、慣容器連接位置降低(φ減小),振動(dòng)控制增強(qiáng)效果越顯著。

    圖8 給出了μ= 0.1%、1%和10%時(shí)的慣容增強(qiáng)區(qū)邊界線??梢园l(fā)現(xiàn),隨著質(zhì)量比增大,慣容增強(qiáng)區(qū)域的邊界向著慣質(zhì)比β 增大、慣容器連接位置降低(φ減小)的方向移動(dòng)。為進(jìn)一步定量給出慣容增強(qiáng)區(qū)的條件,在2.1 節(jié)參數(shù)分析的基礎(chǔ)上對(duì)結(jié)構(gòu)阻尼比ζs= 0.1%、10%進(jìn)行了補(bǔ)充分析。將結(jié)果繪制在以R值為縱軸,β·(1 -φ)2/μ為橫軸的坐標(biāo)圖上,如圖9 所示。從圖中可以看出,一般地,當(dāng)β·(1 -φ)2/μ≥ 1時(shí),有R≤ 1。而R> 1 的情況大多出現(xiàn)在β·(1 -φ)2/μ< 1 時(shí)。因此,在工程上可以偏于保守地將慣容器起增強(qiáng)作用的判定條件表示為:

    為驗(yàn)證上述準(zhǔn)則的有效性,對(duì)其進(jìn)行統(tǒng)計(jì)假設(shè)檢驗(yàn)。分別繪制R> 1 和R< 1 時(shí),判別參數(shù)β·(1 -φ)2/μ的概率分布圖,如圖10 所示。由統(tǒng)計(jì)分析可得,判別式(52)的I 類(lèi)統(tǒng)計(jì)錯(cuò)誤(即在β·(1 -φ)2/μ<1 時(shí)R< 1)概率約為4.85%,但I(xiàn)I 類(lèi)統(tǒng)計(jì)錯(cuò)誤(即在β·(1 -φ)2/μ>1 時(shí)R> 1)的概率低于10-4,說(shuō)明判別式(52)的去偽率較高,具有較強(qiáng)的統(tǒng)計(jì)檢驗(yàn)效能。

    3.2 等效TMD 質(zhì)量比

    TMDI 的一個(gè)重要作用是通過(guò)慣容器提供足夠的質(zhì)量慣性效應(yīng)以實(shí)現(xiàn)動(dòng)力吸振器的輕量化設(shè)計(jì)。本節(jié)以最優(yōu)控制比為等效目標(biāo),重點(diǎn)給出慣容增強(qiáng)區(qū)內(nèi),TMDI 的等效TMD 質(zhì)量比,從而指導(dǎo)動(dòng)力吸振裝置的輕量化設(shè)計(jì)。

    其中:μe為T(mén)MDI 的等效質(zhì)量比,當(dāng)TMD 時(shí)(β =0),取為μ;GTMDI 時(shí)(φ= 0),取為μ+ β;在慣容增強(qiáng)區(qū)的TMDI 時(shí)(β(1 -φ)2/μ≥ 1),取為μ+β(1 -φ)2。

    則慣容器影響系數(shù)R進(jìn)一步表示為:

    將式(57)結(jié)果與3.1 節(jié)的慣容影響系數(shù)R分析結(jié)果進(jìn)行對(duì)比,如圖12 所示。從圖中可以發(fā)現(xiàn),經(jīng)驗(yàn)公式在慣容增強(qiáng)區(qū)(β(1 -φ)2/μ≥ 1)范圍內(nèi)誤差在10%以內(nèi),驗(yàn)證了等效質(zhì)量比的有效性。而當(dāng)β(1 -φ)2/μ< 1 時(shí),由于慣容器無(wú)增強(qiáng)控制效果,不推薦在工程中使用。

    4 算例分析

    本文以某H= 270 m 高混凝土煙囪為例,進(jìn)行TMDI 風(fēng)振控制算例分析。圖13 給出了該混凝土煙囪的截面尺寸(內(nèi)外半徑)隨高度的變化,煙囪的平均外直徑D= 25.34 m。采用梁?jiǎn)卧M煙囪該煙囪沿高度的剛度變化,建立多自由度有限元分析模型,經(jīng)過(guò)模態(tài)分析得到煙囪的前三階自振圓頻率和振型如圖13 所示,可得煙囪橫風(fēng)向渦激共振的臨界風(fēng)速為UCr= 49.4 m/s。

    為確定該煙囪風(fēng)荷載及風(fēng)振響應(yīng),對(duì)其進(jìn)行了剛性模型測(cè)壓風(fēng)洞試驗(yàn)及氣彈模型測(cè)振風(fēng)洞試驗(yàn)[29],風(fēng)洞試驗(yàn)?zāi)M美國(guó)規(guī)范ASCE 7-16[30]中的C 類(lèi)風(fēng)場(chǎng),在交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所大氣邊界層風(fēng)洞實(shí)驗(yàn)室中進(jìn)行。通過(guò)剛性模型測(cè)壓風(fēng)洞試驗(yàn)數(shù)據(jù)分析得到結(jié)構(gòu)順風(fēng)向及橫風(fēng)向的廣義風(fēng)荷載頻譜見(jiàn)圖3,通過(guò)氣彈模型風(fēng)洞試驗(yàn)得到煙囪風(fēng)振響應(yīng)時(shí)程,對(duì)橫風(fēng)向響應(yīng)進(jìn)行隨機(jī)減量法分析,得到自由振動(dòng)衰減曲線,進(jìn)而對(duì)自由衰減曲線的包絡(luò)線進(jìn)行參數(shù)擬合,得到煙囪結(jié)構(gòu)橫風(fēng)向氣動(dòng)阻尼及總阻尼比,如圖14 所示。

    對(duì)于風(fēng)振控制,以質(zhì)量比μ= 1%的TMD 為基準(zhǔn),TMDI 采用慣質(zhì)比β = 20%的慣容器,連接參數(shù)取φ= 0.9、0.75、0.5 和0.0,分別表示慣容增強(qiáng)區(qū)外、慣容增強(qiáng)區(qū)邊界附近、慣容增強(qiáng)區(qū)內(nèi)、GTMDI 4 種情況。采用本文給出的經(jīng)驗(yàn)公式計(jì)算最優(yōu)的頻率比和阻尼比(νopt, ζdopt)和等效質(zhì)量比μe,風(fēng)振控制計(jì)算工況阻尼器詳細(xì)參數(shù)見(jiàn)表1。

    表1 風(fēng)振控制計(jì)算工況阻尼器詳細(xì)參數(shù)Table 1 Detailed damper parameters of wind-induced vibration control cases

    風(fēng)振響應(yīng)計(jì)算分析分別采用多自由度結(jié)構(gòu)的動(dòng)力響應(yīng)時(shí)程分析法以及本文提出的濾波表示解析計(jì)算方法(式(36)、式(40))。

    表2 給出了典型風(fēng)速下(鎖定區(qū)前、臨界風(fēng)速附近、鎖定區(qū)外)的風(fēng)振響應(yīng)計(jì)算結(jié)果對(duì)比??梢钥闯鲲L(fēng)速較低時(shí),由于響應(yīng)本身較小,計(jì)算存在相對(duì)誤差角度,隨著響應(yīng)變大,計(jì)算誤差減小,在臨界風(fēng)速附近和鎖定區(qū)外本文方法和傳統(tǒng)時(shí)程分析相比誤差在±5%以內(nèi),說(shuō)明本文解析方法對(duì)單自由度振動(dòng)為主的高聳煙囪結(jié)構(gòu)具有較強(qiáng)的適用性。雖然表2 表明,傳統(tǒng)時(shí)域分析結(jié)果更接近氣彈試驗(yàn)結(jié)果,但采用本文的解析方法計(jì)算,無(wú)需復(fù)雜迭代計(jì)算,更為快捷,給出的解析公式更便于規(guī)范化表達(dá),同時(shí)也適用于項(xiàng)目初步設(shè)計(jì)階段的快速估算。

    表2 典型風(fēng)速下風(fēng)振響應(yīng)計(jì)算結(jié)果及誤差Table 2 Results and errors of wind-induced responses under typical wind speed cases

    圖16 給出了典型控制工況3、4 下煙囪在臨界風(fēng)速附近(U/UCr= 1.17)橫風(fēng)向共振響應(yīng)時(shí)程的分析結(jié)果,圖中對(duì)比了TMDI 和等效TMD 的響應(yīng)時(shí)程。從圖中可以發(fā)現(xiàn)二者吻合較好,進(jìn)一步驗(yàn)證了本文等效TMD 質(zhì)量比的有效性。根據(jù)ACI307-08 規(guī)范[31]綜合考慮順風(fēng)向極值響應(yīng)Xa、橫風(fēng)向組合響應(yīng)Xc。

    其中:xˉa為順風(fēng)向響應(yīng)的均值;ga、gc分別為順風(fēng)向和橫風(fēng)向響應(yīng)的峰值因子;X(U)為風(fēng)速U時(shí)的最不利響應(yīng)。由于橫風(fēng)向渦激共振受氣動(dòng)阻尼的影響,X(U)并不一定隨U增大而增大,因此下設(shè)計(jì)風(fēng)速Ud下的最不利響應(yīng)Xd取U≤Ud的包絡(luò)值。并利用最不利響應(yīng)Xd計(jì)算風(fēng)振響應(yīng)控制率η,分別見(jiàn)圖17、圖18。

    最不利響應(yīng)Xd表征了在設(shè)計(jì)風(fēng)速Ud下等效靜風(fēng)荷載的取值,控制率η 顯示了采用TMDI 后降低設(shè)計(jì)風(fēng)荷載的百分比。由圖18 可見(jiàn),當(dāng)設(shè)計(jì)風(fēng)速與臨界風(fēng)速相接近時(shí),響應(yīng)控制率達(dá)到最大。采用控制工況3 時(shí),設(shè)計(jì)風(fēng)荷載可降低越35%;控制工況4 可降低設(shè)計(jì)風(fēng)荷載約47%;GTMDI 可降低約59%,但在工程中難以實(shí)施。當(dāng)設(shè)計(jì)風(fēng)速超出渦激共振鎖定區(qū)(Ud> 2UCr)時(shí),采用控制工況3 和工況4 分別可降低設(shè)計(jì)風(fēng)荷載約26%和33%以上。

    5 結(jié)論

    本文從基本運(yùn)動(dòng)方程出發(fā),基于濾波表示法推導(dǎo)了廣義單自由度高聳煙囪結(jié)構(gòu)TMDI 風(fēng)振響應(yīng)控制效率解析表達(dá)式,通過(guò)參數(shù)分析,得到了TMDI 最優(yōu)參數(shù)、最優(yōu)控制比及等效TMD 質(zhì)量比的經(jīng)驗(yàn)公式。并以某270 m 高混凝土煙囪為例驗(yàn)證了理論推導(dǎo)及經(jīng)驗(yàn)公式的有效性,探討了TMDI對(duì)其風(fēng)振響應(yīng)的控制效果,主要結(jié)論如下:

    (1) 高聳結(jié)構(gòu)風(fēng)振響應(yīng)分析時(shí),可將順風(fēng)向和橫風(fēng)向的風(fēng)荷載譜表示為濾波形式,濾波多項(xiàng)式分別為一次和二次多項(xiàng)式,對(duì)于一階模態(tài)主導(dǎo)的結(jié)構(gòu),該方法得到的解析解與傳統(tǒng)時(shí)域分析結(jié)果較吻合。

    (2) 采用本文給出的TMDI 最優(yōu)控制參數(shù)(νopt,ζdopt)關(guān)于質(zhì)量、慣容參數(shù)(μ,β,φ)的經(jīng)驗(yàn)公式計(jì)算的最優(yōu)響應(yīng)控制比Jopt誤差在工程接受范圍內(nèi)。

    (3) 當(dāng)β·(1 -φ)2/μ≥ 1 時(shí),TMDI 的控制效果優(yōu)于相同質(zhì)量比μ下的TMD,說(shuō)明此時(shí)慣容器起到增強(qiáng)控制效果的作用,定義為慣容增強(qiáng)區(qū)。統(tǒng)計(jì)假設(shè)檢驗(yàn)表明,該準(zhǔn)則具有較強(qiáng)的統(tǒng)計(jì)檢驗(yàn)效能。

    (4) 在慣容增強(qiáng)區(qū)內(nèi)(β·(1 -φ)2/μ≥ 1),TMDI可等效為質(zhì)量比為μe的TMD,其中,等效質(zhì)量比μe= β·(1 -φ)2。

    (5) 高聳煙囪算例分析結(jié)果表明,TMDI 對(duì)橫風(fēng)向風(fēng)振響應(yīng)控制效果顯著。當(dāng)設(shè)計(jì)風(fēng)速與臨界風(fēng)速相接近時(shí),響應(yīng)控制率最大,可降低設(shè)計(jì)風(fēng)荷載45%以上;當(dāng)設(shè)計(jì)風(fēng)速超出渦激共振鎖定區(qū)(Ud> 2UCr)時(shí),TMDI 可降低設(shè)計(jì)風(fēng)荷載30%以上。

    猜你喜歡
    風(fēng)振最優(yōu)控制煙囪
    看不見(jiàn)的“煙囪”
    條件平均場(chǎng)隨機(jī)微分方程的最優(yōu)控制問(wèn)題
    帶跳躍平均場(chǎng)倒向隨機(jī)微分方程的線性二次最優(yōu)控制
    張力對(duì)接觸網(wǎng)系統(tǒng)風(fēng)振響應(yīng)影響的風(fēng)洞試驗(yàn)研究
    Timoshenko梁的邊界最優(yōu)控制
    高塔設(shè)備風(fēng)振失效原因分析及改善措施
    煙囪里的客房
    為什么工廠的煙囪都很高?
    采用最優(yōu)控制無(wú)功STATCOM 功率流的解決方案
    煙囪怪
    兒童繪本(2015年10期)2015-08-14 17:54:55
    欧美三级亚洲精品| www.999成人在线观看| 内射极品少妇av片p| 亚洲欧美激情综合另类| 女生性感内裤真人,穿戴方法视频| 国产精品一区二区免费欧美| 日日干狠狠操夜夜爽| 欧美在线黄色| 国产精品1区2区在线观看.| 99热6这里只有精品| 久久久久久国产a免费观看| 午夜福利在线观看免费完整高清在 | 91麻豆av在线| 女人十人毛片免费观看3o分钟| 啦啦啦免费观看视频1| 婷婷亚洲欧美| 久久久久精品国产欧美久久久| 偷拍熟女少妇极品色| 99久久久亚洲精品蜜臀av| 一区二区三区国产精品乱码| av片东京热男人的天堂| 亚洲欧美日韩卡通动漫| 91麻豆精品激情在线观看国产| 又紧又爽又黄一区二区| 国产综合懂色| 午夜精品久久久久久毛片777| 搡老妇女老女人老熟妇| 午夜福利18| 国产精品乱码一区二三区的特点| 欧美一级毛片孕妇| 欧美成人a在线观看| 亚洲激情在线av| 一进一出抽搐gif免费好疼| 欧美最黄视频在线播放免费| 非洲黑人性xxxx精品又粗又长| 久久6这里有精品| 男女床上黄色一级片免费看| 日韩精品青青久久久久久| 少妇高潮的动态图| xxx96com| 熟女少妇亚洲综合色aaa.| 嫁个100分男人电影在线观看| 国产精品三级大全| 亚洲av成人不卡在线观看播放网| 国产熟女xx| 久久九九热精品免费| 黄色视频,在线免费观看| 久久这里只有精品中国| 日韩人妻高清精品专区| 国产精品久久电影中文字幕| www.www免费av| 有码 亚洲区| 一级a爱片免费观看的视频| 国产一级毛片七仙女欲春2| 成人精品一区二区免费| 久久精品国产99精品国产亚洲性色| 在线观看午夜福利视频| 一个人免费在线观看的高清视频| avwww免费| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品乱码久久久v下载方式 | 亚洲成a人片在线一区二区| 美女高潮喷水抽搐中文字幕| 身体一侧抽搐| 看黄色毛片网站| 熟女少妇亚洲综合色aaa.| 成人av在线播放网站| 狠狠狠狠99中文字幕| 久久精品综合一区二区三区| eeuss影院久久| 国内精品美女久久久久久| 日日摸夜夜添夜夜添小说| 毛片女人毛片| 丰满乱子伦码专区| 亚洲自拍偷在线| 麻豆成人av在线观看| 久久久久久久久中文| 亚洲av不卡在线观看| 日本一二三区视频观看| 19禁男女啪啪无遮挡网站| 一a级毛片在线观看| 三级国产精品欧美在线观看| 国产探花极品一区二区| 久久精品国产99精品国产亚洲性色| www.999成人在线观看| 国产午夜精品久久久久久一区二区三区 | 在线播放国产精品三级| 色视频www国产| 国产成人影院久久av| 午夜亚洲福利在线播放| 高清日韩中文字幕在线| 亚洲自拍偷在线| 色哟哟哟哟哟哟| 免费一级毛片在线播放高清视频| 欧美中文综合在线视频| 亚洲人成伊人成综合网2020| 免费观看的影片在线观看| 国产伦精品一区二区三区视频9 | www日本在线高清视频| 欧洲精品卡2卡3卡4卡5卡区| 国产成人欧美在线观看| bbb黄色大片| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区国产精品久久精品| 在线观看一区二区三区| 精品一区二区三区视频在线观看免费| av专区在线播放| 少妇的丰满在线观看| 51午夜福利影视在线观看| 女生性感内裤真人,穿戴方法视频| 国产精品久久久人人做人人爽| 嫩草影视91久久| 国产乱人伦免费视频| 在线国产一区二区在线| 国产激情欧美一区二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 9191精品国产免费久久| 亚洲人成网站高清观看| 欧洲精品卡2卡3卡4卡5卡区| 国内精品久久久久久久电影| 亚洲最大成人手机在线| 欧美在线一区亚洲| 国产av不卡久久| www.色视频.com| 亚洲 国产 在线| 日本黄色视频三级网站网址| 两性午夜刺激爽爽歪歪视频在线观看| 脱女人内裤的视频| 日本在线视频免费播放| 老司机午夜十八禁免费视频| 亚洲精品久久国产高清桃花| 观看美女的网站| 亚洲精品国产精品久久久不卡| 精品久久久久久久毛片微露脸| 国产高清激情床上av| 国产成人福利小说| 超碰av人人做人人爽久久 | 一区福利在线观看| 热99在线观看视频| 国产蜜桃级精品一区二区三区| 五月伊人婷婷丁香| 日韩欧美三级三区| 国产伦精品一区二区三区视频9 | 免费电影在线观看免费观看| 成人三级黄色视频| 久久99热这里只有精品18| 国产三级黄色录像| 淫秽高清视频在线观看| av欧美777| 在线免费观看的www视频| 深夜精品福利| 婷婷六月久久综合丁香| 99精品欧美一区二区三区四区| 搡老岳熟女国产| 在线观看一区二区三区| 国产高清视频在线观看网站| 国产精品免费一区二区三区在线| 一个人看的www免费观看视频| 亚洲av二区三区四区| 757午夜福利合集在线观看| 日韩欧美 国产精品| 国产高清视频在线播放一区| 国产一区在线观看成人免费| 国产主播在线观看一区二区| 一a级毛片在线观看| 亚洲成人久久性| 亚洲av五月六月丁香网| 国产免费男女视频| 久久久久久久久中文| 啦啦啦韩国在线观看视频| 一本精品99久久精品77| 婷婷丁香在线五月| 国产av一区在线观看免费| 99热6这里只有精品| 十八禁人妻一区二区| 欧美zozozo另类| 久久99热这里只有精品18| 久久久久久人人人人人| 一区二区三区高清视频在线| 麻豆成人av在线观看| 日韩中文字幕欧美一区二区| 亚洲精品在线美女| 12—13女人毛片做爰片一| 男人舔奶头视频| 国产真实伦视频高清在线观看 | www日本在线高清视频| 一区二区三区激情视频| 亚洲成人精品中文字幕电影| 男女下面进入的视频免费午夜| 在线免费观看的www视频| 男女床上黄色一级片免费看| 97超视频在线观看视频| 亚洲av美国av| 亚洲真实伦在线观看| 精品人妻偷拍中文字幕| 亚洲中文字幕日韩| 国产精品99久久99久久久不卡| 757午夜福利合集在线观看| 两个人看的免费小视频| 免费人成在线观看视频色| 中文资源天堂在线| 亚洲国产精品合色在线| 久久久久久人人人人人| 亚洲精品日韩av片在线观看 | 国产午夜精品论理片| 中文资源天堂在线| 成年女人看的毛片在线观看| 身体一侧抽搐| 国产爱豆传媒在线观看| 欧美+亚洲+日韩+国产| 亚洲自拍偷在线| 国产高清激情床上av| 男女床上黄色一级片免费看| 在线观看美女被高潮喷水网站 | a在线观看视频网站| 啪啪无遮挡十八禁网站| av视频在线观看入口| 亚洲欧美一区二区三区黑人| 欧美性猛交╳xxx乱大交人| 琪琪午夜伦伦电影理论片6080| 欧美黑人欧美精品刺激| 午夜免费观看网址| 男女下面进入的视频免费午夜| 搡老岳熟女国产| ponron亚洲| 午夜福利高清视频| 成年版毛片免费区| 三级男女做爰猛烈吃奶摸视频| 婷婷精品国产亚洲av| 免费看光身美女| 国产色爽女视频免费观看| 亚洲不卡免费看| 亚洲在线自拍视频| 国产成人啪精品午夜网站| 国产单亲对白刺激| 黄色丝袜av网址大全| 欧美日韩亚洲国产一区二区在线观看| 女人被狂操c到高潮| 国产精品亚洲av一区麻豆| 女警被强在线播放| 欧美日韩一级在线毛片| 午夜激情福利司机影院| 可以在线观看的亚洲视频| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕精品亚洲无线码一区| 成人亚洲精品av一区二区| 观看免费一级毛片| 99热这里只有是精品50| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 欧美激情久久久久久爽电影| 日本黄大片高清| av国产免费在线观看| 99热只有精品国产| 综合色av麻豆| 久久精品综合一区二区三区| 亚洲国产精品999在线| 午夜福利欧美成人| 精华霜和精华液先用哪个| 亚洲精品成人久久久久久| 亚洲精品美女久久久久99蜜臀| 搡老妇女老女人老熟妇| 女人被狂操c到高潮| 哪里可以看免费的av片| 少妇裸体淫交视频免费看高清| 日韩欧美精品免费久久 | 欧美日韩综合久久久久久 | 特大巨黑吊av在线直播| 午夜精品一区二区三区免费看| 婷婷亚洲欧美| 亚洲精品在线美女| 精华霜和精华液先用哪个| 97人妻精品一区二区三区麻豆| 婷婷六月久久综合丁香| 久久精品影院6| 亚洲五月天丁香| 亚洲avbb在线观看| 一本精品99久久精品77| 久久久久久久亚洲中文字幕 | av视频在线观看入口| 亚洲国产色片| 亚洲成a人片在线一区二区| 欧美一级毛片孕妇| 国产又黄又爽又无遮挡在线| 久久久久精品国产欧美久久久| 最新美女视频免费是黄的| 男女那种视频在线观看| 法律面前人人平等表现在哪些方面| 两个人看的免费小视频| 99久久精品一区二区三区| 婷婷精品国产亚洲av| 欧美日韩精品网址| 一区二区三区高清视频在线| 欧美在线黄色| 黄色视频,在线免费观看| 老鸭窝网址在线观看| 久99久视频精品免费| 国产高清三级在线| 在线播放无遮挡| 手机成人av网站| 两个人的视频大全免费| 999久久久精品免费观看国产| 在线视频色国产色| 黄片小视频在线播放| 亚洲性夜色夜夜综合| 男插女下体视频免费在线播放| 久99久视频精品免费| 蜜桃久久精品国产亚洲av| 男女那种视频在线观看| 国产亚洲精品av在线| 在线观看舔阴道视频| 黄色日韩在线| 又粗又爽又猛毛片免费看| 亚洲成人精品中文字幕电影| 黄片小视频在线播放| 欧美性猛交╳xxx乱大交人| 国产v大片淫在线免费观看| 欧美黑人巨大hd| 中文字幕精品亚洲无线码一区| 日本一本二区三区精品| 午夜福利在线观看免费完整高清在 | 午夜两性在线视频| 久久99热这里只有精品18| 欧美黑人巨大hd| 在线播放国产精品三级| 国产激情偷乱视频一区二区| 日本免费一区二区三区高清不卡| 嫩草影院精品99| 国产美女午夜福利| 又黄又爽又免费观看的视频| 日本 av在线| 久久久久久久精品吃奶| 国产69精品久久久久777片| 91在线观看av| 婷婷精品国产亚洲av在线| 久久精品综合一区二区三区| 一进一出抽搐动态| 国产av麻豆久久久久久久| 中文字幕人成人乱码亚洲影| 国产精品美女特级片免费视频播放器| 757午夜福利合集在线观看| 国产一区二区在线av高清观看| 精品日产1卡2卡| 亚洲国产欧美网| 国产真人三级小视频在线观看| 欧美性猛交黑人性爽| 日韩欧美免费精品| 免费看日本二区| 搡老熟女国产l中国老女人| 在线观看av片永久免费下载| 久久久久久大精品| 亚洲人与动物交配视频| 国产探花在线观看一区二区| 亚洲专区中文字幕在线| 久久香蕉国产精品| 18禁黄网站禁片午夜丰满| 中文字幕高清在线视频| 男人舔奶头视频| 哪里可以看免费的av片| 操出白浆在线播放| 日韩欧美在线乱码| 亚洲 国产 在线| 亚洲欧美一区二区三区黑人| 99热这里只有精品一区| 少妇的逼水好多| 欧美不卡视频在线免费观看| 搡老岳熟女国产| 国模一区二区三区四区视频| 婷婷六月久久综合丁香| 黄色视频,在线免费观看| 无人区码免费观看不卡| 亚洲国产日韩欧美精品在线观看 | 成人特级av手机在线观看| 无人区码免费观看不卡| 精品久久久久久久人妻蜜臀av| 久久久久亚洲av毛片大全| 欧美性感艳星| www国产在线视频色| 欧美精品啪啪一区二区三区| 一级a爱片免费观看的视频| а√天堂www在线а√下载| 特大巨黑吊av在线直播| 婷婷亚洲欧美| 精品国产美女av久久久久小说| 一本一本综合久久| 精品国内亚洲2022精品成人| 日韩欧美 国产精品| 综合色av麻豆| 国产v大片淫在线免费观看| 激情在线观看视频在线高清| 国内揄拍国产精品人妻在线| a在线观看视频网站| 欧美一级a爱片免费观看看| 日本免费一区二区三区高清不卡| 中文字幕熟女人妻在线| 欧美日韩亚洲国产一区二区在线观看| 成年版毛片免费区| 亚洲激情在线av| 亚洲最大成人手机在线| 人妻丰满熟妇av一区二区三区| 国内少妇人妻偷人精品xxx网站| 综合色av麻豆| 欧美另类亚洲清纯唯美| 九九在线视频观看精品| 国产一区二区在线av高清观看| x7x7x7水蜜桃| 在线十欧美十亚洲十日本专区| 久久这里只有精品中国| 男人的好看免费观看在线视频| 精品久久久久久久久久免费视频| 一进一出抽搐gif免费好疼| 综合色av麻豆| 99久久精品国产亚洲精品| 熟妇人妻久久中文字幕3abv| 国产不卡一卡二| 中国美女看黄片| 免费在线观看日本一区| 亚洲成人久久性| 天天添夜夜摸| 美女高潮的动态| 岛国在线观看网站| 一进一出好大好爽视频| 麻豆国产av国片精品| 小蜜桃在线观看免费完整版高清| 午夜影院日韩av| 免费观看人在逋| 国产成人欧美在线观看| 手机成人av网站| 美女高潮喷水抽搐中文字幕| 狂野欧美白嫩少妇大欣赏| 国产伦一二天堂av在线观看| 热99在线观看视频| а√天堂www在线а√下载| 亚洲人成网站在线播| 熟女电影av网| 久久精品综合一区二区三区| 女人被狂操c到高潮| 国产乱人伦免费视频| 日韩欧美三级三区| bbb黄色大片| 欧美不卡视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 首页视频小说图片口味搜索| 欧美日韩精品网址| 深夜精品福利| 国产成人aa在线观看| 内射极品少妇av片p| 香蕉av资源在线| 18禁黄网站禁片午夜丰满| 日本熟妇午夜| 亚洲,欧美精品.| 免费在线观看成人毛片| 一个人看视频在线观看www免费 | 国产精品美女特级片免费视频播放器| 91九色精品人成在线观看| 婷婷六月久久综合丁香| 欧美色视频一区免费| 两性午夜刺激爽爽歪歪视频在线观看| 在线观看av片永久免费下载| 免费观看人在逋| 精品人妻一区二区三区麻豆 | 午夜视频国产福利| 色噜噜av男人的天堂激情| tocl精华| 亚洲va日本ⅴa欧美va伊人久久| 欧美性猛交黑人性爽| 黄色丝袜av网址大全| 村上凉子中文字幕在线| 亚洲av熟女| 国产色爽女视频免费观看| 欧美+亚洲+日韩+国产| 久久6这里有精品| 亚洲第一电影网av| 国内精品美女久久久久久| 少妇裸体淫交视频免费看高清| 中文字幕人妻熟人妻熟丝袜美 | 日韩av在线大香蕉| 亚洲男人的天堂狠狠| 亚洲成人免费电影在线观看| 国产伦精品一区二区三区视频9 | 首页视频小说图片口味搜索| 日本黄色视频三级网站网址| 亚洲无线在线观看| 欧美中文综合在线视频| 欧美+亚洲+日韩+国产| 中亚洲国语对白在线视频| 中文字幕高清在线视频| 久久精品国产亚洲av香蕉五月| 国产欧美日韩一区二区精品| 亚洲欧美日韩东京热| 长腿黑丝高跟| 91在线精品国自产拍蜜月 | 久久精品国产自在天天线| 精品欧美国产一区二区三| 婷婷亚洲欧美| 露出奶头的视频| 国产在视频线在精品| 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 久久久久久久久中文| 日韩欧美三级三区| 老司机福利观看| 1000部很黄的大片| 国产伦人伦偷精品视频| 精品一区二区三区av网在线观看| 草草在线视频免费看| 非洲黑人性xxxx精品又粗又长| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 免费观看精品视频网站| 国产伦精品一区二区三区视频9 | av专区在线播放| 久久99热这里只有精品18| 国产成人欧美在线观看| 两个人的视频大全免费| 天天躁日日操中文字幕| netflix在线观看网站| 波多野结衣高清作品| 亚洲男人的天堂狠狠| 人人妻人人澡欧美一区二区| 一本一本综合久久| 黄色成人免费大全| 日本a在线网址| 亚洲一区二区三区不卡视频| 国产成人福利小说| 制服人妻中文乱码| 91av网一区二区| 无人区码免费观看不卡| 国产极品精品免费视频能看的| 久久久久久久精品吃奶| 亚洲午夜理论影院| 日韩欧美三级三区| 日韩欧美在线二视频| 美女大奶头视频| 国内精品久久久久久久电影| av福利片在线观看| 欧美激情在线99| 欧美乱色亚洲激情| 亚洲熟妇熟女久久| 成人欧美大片| 天堂动漫精品| 尤物成人国产欧美一区二区三区| 日本 av在线| 91麻豆av在线| 2021天堂中文幕一二区在线观| 久久久久九九精品影院| 日韩精品中文字幕看吧| 叶爱在线成人免费视频播放| 看片在线看免费视频| 免费高清视频大片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜久久久久精精品| 国产一区在线观看成人免费| 两人在一起打扑克的视频| 久久精品人妻少妇| 亚洲欧美日韩高清专用| 全区人妻精品视频| 成人无遮挡网站| 白带黄色成豆腐渣| 欧美成人性av电影在线观看| 国产亚洲精品一区二区www| 国产精品一区二区三区四区久久| 国产高清videossex| 男女做爰动态图高潮gif福利片| 欧美日韩亚洲国产一区二区在线观看| 午夜福利18| 两人在一起打扑克的视频| 欧美又色又爽又黄视频| 熟女人妻精品中文字幕| av在线蜜桃| 精品欧美国产一区二区三| 婷婷亚洲欧美| 亚洲欧美日韩卡通动漫| 色老头精品视频在线观看| 叶爱在线成人免费视频播放| 欧洲精品卡2卡3卡4卡5卡区| 老熟妇仑乱视频hdxx| 国产在线精品亚洲第一网站| 一级a爱片免费观看的视频| 夜夜夜夜夜久久久久| 白带黄色成豆腐渣| 中国美女看黄片| 熟女电影av网| 欧美在线黄色| 色尼玛亚洲综合影院| 免费电影在线观看免费观看| 熟女人妻精品中文字幕| 久久久久久人人人人人| 久久精品综合一区二区三区| av视频在线观看入口| 97超级碰碰碰精品色视频在线观看| 国语自产精品视频在线第100页| 久久精品国产亚洲av香蕉五月| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 桃色一区二区三区在线观看| 亚洲国产高清在线一区二区三| 制服丝袜大香蕉在线| 亚洲欧美日韩东京热| 午夜免费观看网址| 色视频www国产| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 午夜两性在线视频| 少妇的逼好多水| 99热精品在线国产| 欧美一区二区国产精品久久精品| 久久精品国产99精品国产亚洲性色| 欧美日本亚洲视频在线播放| 久久性视频一级片| 国产高潮美女av| 亚洲精品影视一区二区三区av| 黄色日韩在线| 国产黄a三级三级三级人| a级一级毛片免费在线观看|