• <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在线а√下载| 精品一区二区三区视频在线观看免费| 日本在线视频免费播放| 琪琪午夜伦伦电影理论片6080| 欧美精品啪啪一区二区三区| 香蕉国产在线看| 久久精品人妻少妇| 精品一区二区三区av网在线观看| 1024视频免费在线观看| 国产在线观看jvid| 久久精品亚洲精品国产色婷小说| 久久精品国产99精品国产亚洲性色| 亚洲久久久国产精品| 日本成人三级电影网站| 91麻豆av在线| 久久亚洲真实| 成人三级黄色视频| 成人永久免费在线观看视频| 亚洲成国产人片在线观看| 成在线人永久免费视频| 日韩精品青青久久久久久| 成人永久免费在线观看视频| 老司机福利观看| 中国美女看黄片| 波多野结衣高清作品| 亚洲 欧美一区二区三区| 欧美激情极品国产一区二区三区| 精品国产超薄肉色丝袜足j| 日本三级黄在线观看| 男女之事视频高清在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲在线自拍视频| 日韩大码丰满熟妇| 国产一区在线观看成人免费| 少妇熟女aⅴ在线视频| 精品日产1卡2卡| 母亲3免费完整高清在线观看| 国产亚洲精品第一综合不卡| 好看av亚洲va欧美ⅴa在| 精品高清国产在线一区| 99热6这里只有精品| 欧美性猛交╳xxx乱大交人| 精华霜和精华液先用哪个| 免费在线观看视频国产中文字幕亚洲| 国产视频内射| 男女之事视频高清在线观看| 一区二区日韩欧美中文字幕| 日本黄色视频三级网站网址| 淫妇啪啪啪对白视频| 欧美日韩福利视频一区二区| 精品免费久久久久久久清纯| 日本 av在线| 色哟哟哟哟哟哟| 美国免费a级毛片| 成年女人毛片免费观看观看9| 日本a在线网址| 亚洲国产精品sss在线观看| 757午夜福利合集在线观看| 中文字幕最新亚洲高清| 一区二区三区精品91| 啦啦啦 在线观看视频| 自线自在国产av| 美国免费a级毛片| 成人18禁在线播放| 一本久久中文字幕| 欧美性猛交黑人性爽| 男女下面进入的视频免费午夜 | 成人永久免费在线观看视频| 日韩一卡2卡3卡4卡2021年| 欧美zozozo另类| 国产精品免费一区二区三区在线| 国产爱豆传媒在线观看 | 日韩欧美国产在线观看| 精品熟女少妇八av免费久了| 男女之事视频高清在线观看| 欧美黄色片欧美黄色片| 人人妻人人澡欧美一区二区| 久热爱精品视频在线9| 日韩大尺度精品在线看网址| 在线看三级毛片| а√天堂www在线а√下载| 久久亚洲精品不卡| 日韩欧美 国产精品| 曰老女人黄片| 亚洲精品av麻豆狂野| 成年女人毛片免费观看观看9| 欧美在线一区亚洲| 长腿黑丝高跟| 国产在线观看jvid| 中文亚洲av片在线观看爽| 91九色精品人成在线观看| 精品一区二区三区av网在线观看| av片东京热男人的天堂| 免费人成视频x8x8入口观看| 免费电影在线观看免费观看| 日韩精品青青久久久久久| 又黄又爽又免费观看的视频| 啦啦啦观看免费观看视频高清| 99热6这里只有精品| 一边摸一边抽搐一进一小说| 91麻豆精品激情在线观看国产| 久久久久久人人人人人| 中文资源天堂在线| √禁漫天堂资源中文www| 好男人在线观看高清免费视频 | 在线av久久热| 91麻豆av在线| av超薄肉色丝袜交足视频| 久久香蕉国产精品| 中文字幕人成人乱码亚洲影| 又紧又爽又黄一区二区| 精品卡一卡二卡四卡免费| 国内久久婷婷六月综合欲色啪| 一本久久中文字幕| 日本a在线网址| 日韩精品免费视频一区二区三区| 操出白浆在线播放| 欧美性长视频在线观看| 变态另类丝袜制服| 国产午夜福利久久久久久| 久久久久精品国产欧美久久久| 黄片大片在线免费观看| 免费看美女性在线毛片视频| 亚洲激情在线av| 欧美在线黄色| 亚洲成人久久性| 免费在线观看视频国产中文字幕亚洲| 夜夜躁狠狠躁天天躁| 国产成人欧美在线观看| 真人一进一出gif抽搐免费| 男女之事视频高清在线观看| 国产精品九九99| 午夜免费观看网址| 午夜福利在线在线| 999精品在线视频| 中文字幕精品免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三区四区五区乱码| a在线观看视频网站| cao死你这个sao货| 精品国产乱码久久久久久男人| 亚洲第一青青草原| 嫩草影视91久久| 国产一区二区三区视频了| 欧美激情久久久久久爽电影| 日本一本二区三区精品| 最新美女视频免费是黄的| 在线观看免费午夜福利视频| 国产精品日韩av在线免费观看| 最好的美女福利视频网| 又黄又粗又硬又大视频| 国产成人啪精品午夜网站| 很黄的视频免费| 搡老妇女老女人老熟妇| 一进一出好大好爽视频| 波多野结衣高清无吗| 精品久久久久久久人妻蜜臀av| xxxwww97欧美| 午夜福利视频1000在线观看| 国产亚洲精品久久久久5区| 97碰自拍视频| 91大片在线观看| bbb黄色大片| 一夜夜www| 亚洲熟妇中文字幕五十中出| 成人国产一区最新在线观看| 白带黄色成豆腐渣| 精品久久久久久久人妻蜜臀av| 色尼玛亚洲综合影院| 三级毛片av免费| 非洲黑人性xxxx精品又粗又长| 别揉我奶头~嗯~啊~动态视频| 一本久久中文字幕| www.www免费av| a在线观看视频网站| 夜夜爽天天搞| 老鸭窝网址在线观看| 国产伦人伦偷精品视频| 狠狠狠狠99中文字幕| 亚洲九九香蕉| 免费高清在线观看日韩| 老汉色av国产亚洲站长工具| 中文字幕人妻丝袜一区二区| 国产精品,欧美在线| 两个人视频免费观看高清| 很黄的视频免费| 高潮久久久久久久久久久不卡| 日本熟妇午夜| 亚洲美女黄片视频| 97人妻精品一区二区三区麻豆 | 999精品在线视频| 777久久人妻少妇嫩草av网站| 亚洲自偷自拍图片 自拍| 可以在线观看的亚洲视频| 两性夫妻黄色片| 大香蕉久久成人网| 色精品久久人妻99蜜桃| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久5区| 俄罗斯特黄特色一大片| 精品乱码久久久久久99久播| 草草在线视频免费看| 国产精品久久久久久精品电影 | 最近在线观看免费完整版| 久久午夜亚洲精品久久| 国产真实乱freesex| 久久午夜亚洲精品久久| av福利片在线| 国产麻豆成人av免费视频| e午夜精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看| 久久久久久久久中文| 熟妇人妻久久中文字幕3abv| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看视频国产中文字幕亚洲| 欧美 亚洲 国产 日韩一| 熟妇人妻久久中文字幕3abv| 女性被躁到高潮视频| 黑人巨大精品欧美一区二区mp4| 在线观看免费视频日本深夜| 久久久久久免费高清国产稀缺| 丝袜在线中文字幕| 亚洲国产欧美日韩在线播放| 午夜福利高清视频| 国产精品影院久久| 久久亚洲精品不卡| 哪里可以看免费的av片| 午夜福利免费观看在线| 少妇裸体淫交视频免费看高清 | 一区二区三区高清视频在线| av天堂在线播放| 午夜福利欧美成人| 日本精品一区二区三区蜜桃| 国产三级黄色录像| 99国产综合亚洲精品| 51午夜福利影视在线观看| 1024手机看黄色片| 老司机午夜十八禁免费视频| 亚洲av电影不卡..在线观看| 久久国产精品人妻蜜桃| 一级毛片精品| 亚洲午夜理论影院| 国产精品乱码一区二三区的特点| 又紧又爽又黄一区二区| 国产爱豆传媒在线观看 | 欧美日韩精品网址| 国产精品98久久久久久宅男小说| 国产亚洲精品久久久久久毛片| 亚洲 国产 在线| 精华霜和精华液先用哪个| 国产主播在线观看一区二区| 日本 av在线| 久久久久久久久久黄片| a级毛片a级免费在线| 久久久久久大精品| 男女下面进入的视频免费午夜 | 亚洲av成人av| 99热只有精品国产| 亚洲国产看品久久| www.自偷自拍.com| 亚洲国产毛片av蜜桃av| 美女高潮喷水抽搐中文字幕| 一夜夜www| 欧洲精品卡2卡3卡4卡5卡区| 女人爽到高潮嗷嗷叫在线视频| 亚洲专区国产一区二区| 久久久久亚洲av毛片大全| 国产真人三级小视频在线观看| 男女下面进入的视频免费午夜 | 老汉色av国产亚洲站长工具| 美女高潮到喷水免费观看| 黄色视频,在线免费观看| 国产成人精品久久二区二区91| 在线观看午夜福利视频| 国产av又大| x7x7x7水蜜桃| 欧美另类亚洲清纯唯美| 夜夜躁狠狠躁天天躁| 丁香六月欧美| 欧美日本视频| 国产激情偷乱视频一区二区| 99国产精品99久久久久| 天天一区二区日本电影三级| 欧美日韩精品网址| 不卡av一区二区三区| 婷婷丁香在线五月| 国产亚洲欧美精品永久| 国产日本99.免费观看| 久久精品91蜜桃| 亚洲一区中文字幕在线| 好看av亚洲va欧美ⅴa在| 女人爽到高潮嗷嗷叫在线视频| 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 中文字幕久久专区| 制服人妻中文乱码| 中文字幕最新亚洲高清| 亚洲人成网站高清观看| 亚洲 欧美一区二区三区| 国产成人精品久久二区二区免费| 老司机福利观看| 国产亚洲欧美精品永久| 免费在线观看日本一区| 国产成+人综合+亚洲专区| 12—13女人毛片做爰片一| 日本五十路高清| 亚洲真实伦在线观看| 亚洲精品在线美女| 国产人伦9x9x在线观看| 精品卡一卡二卡四卡免费| 成人欧美大片| 黄色毛片三级朝国网站| 欧美大码av| 国产精品久久电影中文字幕| 狂野欧美激情性xxxx| 国产激情久久老熟女| 啪啪无遮挡十八禁网站| 人成视频在线观看免费观看| 18禁裸乳无遮挡免费网站照片 | 午夜两性在线视频| 老汉色∧v一级毛片| 一进一出抽搐gif免费好疼| 久久久久久人人人人人| 精品一区二区三区视频在线观看免费| 不卡一级毛片| 亚洲成人久久性| 变态另类成人亚洲欧美熟女| av免费在线观看网站| 精品一区二区三区视频在线观看免费| 欧美又色又爽又黄视频| 成人av一区二区三区在线看| 麻豆久久精品国产亚洲av| 久久久久亚洲av毛片大全| 久久婷婷成人综合色麻豆| 国语自产精品视频在线第100页| 国产精品乱码一区二三区的特点| 国产亚洲av高清不卡| 亚洲av成人一区二区三| 色婷婷久久久亚洲欧美| 国产精品亚洲美女久久久| 又黄又粗又硬又大视频| 侵犯人妻中文字幕一二三四区| 成人av一区二区三区在线看| 欧美精品亚洲一区二区| www.自偷自拍.com| 亚洲国产看品久久| 午夜福利一区二区在线看| 国产一卡二卡三卡精品| 日本 av在线| 欧美日韩亚洲综合一区二区三区_| 国产又色又爽无遮挡免费看| 欧美日韩乱码在线| 美女国产高潮福利片在线看| 正在播放国产对白刺激| 久久久久精品国产欧美久久久| 欧美最黄视频在线播放免费| 亚洲国产日韩欧美精品在线观看 | 人人妻人人澡人人看| 午夜福利视频1000在线观看| 一本久久中文字幕| 别揉我奶头~嗯~啊~动态视频| 给我免费播放毛片高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一区福利在线观看| 久久狼人影院| 伊人久久大香线蕉亚洲五| 国产欧美日韩精品亚洲av| 亚洲av熟女| 国产成人影院久久av| a级毛片a级免费在线| 丰满的人妻完整版| 亚洲国产毛片av蜜桃av| 国产精品,欧美在线| 琪琪午夜伦伦电影理论片6080| 女性生殖器流出的白浆| 成人一区二区视频在线观看| 亚洲欧美一区二区三区黑人| 日韩高清综合在线| 可以在线观看的亚洲视频| 在线av久久热| 一本精品99久久精品77| 欧美一级a爱片免费观看看 | 国产激情欧美一区二区| 黄色成人免费大全| 日韩欧美 国产精品| 久久国产乱子伦精品免费另类| 日韩精品免费视频一区二区三区| 宅男免费午夜| 久久人妻福利社区极品人妻图片| 久久天躁狠狠躁夜夜2o2o| 成年人黄色毛片网站| 男女那种视频在线观看| 91九色精品人成在线观看| 成人永久免费在线观看视频| 熟女电影av网| 一级a爱片免费观看的视频| 久久久久久大精品| 久久久久国内视频| 国内久久婷婷六月综合欲色啪| 91在线观看av| 欧美激情久久久久久爽电影| 悠悠久久av| 国产又色又爽无遮挡免费看| 岛国视频午夜一区免费看| 黄频高清免费视频| 国产国语露脸激情在线看| 精品久久久久久久久久久久久 | 日韩av在线大香蕉| 在线观看日韩欧美| 天堂影院成人在线观看| 久久久久国产精品人妻aⅴ院| 老熟妇仑乱视频hdxx| 久久99热这里只有精品18| 午夜福利免费观看在线| 国产精华一区二区三区| 国产视频内射| 国产1区2区3区精品| 搡老熟女国产l中国老女人| 亚洲自偷自拍图片 自拍| 91字幕亚洲| 久久天堂一区二区三区四区| 久久久精品欧美日韩精品| 日韩大码丰满熟妇| aaaaa片日本免费| 中文字幕精品免费在线观看视频| 伊人久久大香线蕉亚洲五| 亚洲成av片中文字幕在线观看| 亚洲国产精品成人综合色| 午夜福利视频1000在线观看| 黄色a级毛片大全视频| 欧美成人午夜精品| av在线播放免费不卡| 叶爱在线成人免费视频播放| 成年免费大片在线观看| 久久精品亚洲精品国产色婷小说| 桃色一区二区三区在线观看| 在线av久久热| 正在播放国产对白刺激| 亚洲欧洲精品一区二区精品久久久| 一区二区三区国产精品乱码| 欧美不卡视频在线免费观看 | 亚洲精品粉嫩美女一区| 国产视频内射| 午夜日韩欧美国产| 美女免费视频网站| 天天添夜夜摸| 特大巨黑吊av在线直播 | 人妻久久中文字幕网| 男人操女人黄网站| 99精品久久久久人妻精品| 一区二区三区国产精品乱码| 黄色 视频免费看| av天堂在线播放| 美女大奶头视频| 手机成人av网站| 色老头精品视频在线观看| 欧美亚洲日本最大视频资源| 日韩 欧美 亚洲 中文字幕| 成在线人永久免费视频| 精品午夜福利视频在线观看一区| 手机成人av网站| 亚洲精品久久成人aⅴ小说| 一本久久中文字幕| 亚洲男人的天堂狠狠| 搞女人的毛片| www日本黄色视频网| 叶爱在线成人免费视频播放| 一进一出抽搐动态| 最新在线观看一区二区三区| 在线播放国产精品三级| 老司机在亚洲福利影院| 国产主播在线观看一区二区| 最近在线观看免费完整版| 国产v大片淫在线免费观看| 欧美乱码精品一区二区三区| 丁香六月欧美| 国产午夜福利久久久久久| 美女高潮喷水抽搐中文字幕| 国产精品av久久久久免费| 欧美日韩福利视频一区二区| 黄色丝袜av网址大全| av电影中文网址| 男女午夜视频在线观看| 女性生殖器流出的白浆| 亚洲成av人片免费观看| 国产在线精品亚洲第一网站| 757午夜福利合集在线观看| 久久久久亚洲av毛片大全| 国产精品一区二区三区四区久久 | 国产欧美日韩一区二区精品| 熟女少妇亚洲综合色aaa.| 一级毛片高清免费大全| 久久久久免费精品人妻一区二区 | 国产精品 欧美亚洲| 欧美精品亚洲一区二区| 91麻豆精品激情在线观看国产| 国产精品 国内视频| 精品久久久久久成人av| 欧美日韩精品网址| 亚洲激情在线av| aaaaa片日本免费| 国产免费男女视频| 国内揄拍国产精品人妻在线 | 国产真人三级小视频在线观看| 两个人视频免费观看高清| 精品久久久久久成人av| 一级黄色大片毛片| 免费av毛片视频| 别揉我奶头~嗯~啊~动态视频| 九色国产91popny在线| 国产精品综合久久久久久久免费| videosex国产| 免费看日本二区| 在线观看66精品国产| 中文字幕人成人乱码亚洲影| 性色av乱码一区二区三区2| av在线天堂中文字幕| 99热6这里只有精品| 精品久久久久久成人av| 黑人巨大精品欧美一区二区mp4| 欧美黑人欧美精品刺激| 亚洲成人久久爱视频| 色尼玛亚洲综合影院| 国产亚洲精品综合一区在线观看 | 亚洲精品国产区一区二| 男女下面进入的视频免费午夜 | 亚洲国产欧洲综合997久久, | 在线观看免费日韩欧美大片| 91老司机精品| 日韩免费av在线播放| 国产私拍福利视频在线观看| 两人在一起打扑克的视频| 国产精品国产高清国产av| 亚洲精品美女久久av网站| 亚洲免费av在线视频| 女同久久另类99精品国产91| 国产黄a三级三级三级人| 免费在线观看成人毛片| 99久久99久久久精品蜜桃| 搡老岳熟女国产| 亚洲精品久久国产高清桃花| 久99久视频精品免费| 国产黄片美女视频| 欧美黑人欧美精品刺激| 一区二区三区精品91| www.精华液| 欧美乱色亚洲激情| 亚洲精品美女久久久久99蜜臀| 欧美色欧美亚洲另类二区| 国产精品免费视频内射| 变态另类丝袜制服| 啪啪无遮挡十八禁网站| 国产精品久久电影中文字幕| 不卡一级毛片| 免费高清在线观看日韩| 欧美中文日本在线观看视频| 成人精品一区二区免费| 欧美在线黄色| 两性夫妻黄色片| 十分钟在线观看高清视频www| 欧美久久黑人一区二区| 欧美绝顶高潮抽搐喷水| 日本熟妇午夜| 国产精品99久久99久久久不卡| 一级毛片精品| xxx96com| 波多野结衣高清作品| 亚洲成人久久爱视频| 国产免费男女视频| 国产91精品成人一区二区三区| 麻豆成人午夜福利视频| 一本大道久久a久久精品| 欧美又色又爽又黄视频| 欧美黑人巨大hd| 91大片在线观看| 亚洲精品一区av在线观看| 最近最新免费中文字幕在线| 国内精品久久久久精免费| 午夜福利18| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 日韩欧美 国产精品| 国产精品久久久久久精品电影 | 国产又色又爽无遮挡免费看| 久久精品影院6| 久久久精品国产亚洲av高清涩受| 中文在线观看免费www的网站 | 中文字幕另类日韩欧美亚洲嫩草| 男人的好看免费观看在线视频 | 免费观看人在逋| svipshipincom国产片| 大型黄色视频在线免费观看| 亚洲一码二码三码区别大吗| 午夜久久久久精精品| 免费在线观看影片大全网站| 国产高清激情床上av| 亚洲五月天丁香| 亚洲成a人片在线一区二区| 白带黄色成豆腐渣| a级毛片a级免费在线| 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| 精品欧美国产一区二区三| 搞女人的毛片| 久久久久久国产a免费观看| 国产av一区在线观看免费| 女人高潮潮喷娇喘18禁视频| 日韩高清综合在线| 欧美 亚洲 国产 日韩一| 国产精品日韩av在线免费观看| 久久伊人香网站| 不卡av一区二区三区| 国产三级黄色录像|