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

    曼型干式煤氣柜動力特性分析及風(fēng)振響應(yīng)主要貢獻(xiàn)模態(tài)識別

    2016-01-15 02:58:38李正良,劉欣鵬,晏致濤
    振動與沖擊 2015年19期
    關(guān)鍵詞:動力特性

    曼型干式煤氣柜動力特性分析及風(fēng)振響應(yīng)主要貢獻(xiàn)模態(tài)識別

    李正良1,2, 劉欣鵬1, 晏致濤1,2, 焦紅偉1, 俞登科1

    (1.重慶大學(xué)土木工程學(xué)院,重慶400045;2.重慶大學(xué)山地城鎮(zhèn)建設(shè)與新技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,重慶400045)

    摘要:曼型干式煤氣柜不同于一般空間結(jié)構(gòu),其動力特性十分復(fù)雜,作為研究其頻域內(nèi)風(fēng)振響應(yīng)的基礎(chǔ),首先通過建立精細(xì)化有限元模型研究其在不同工況下的動力特性;其次,分別計算背景、共振響應(yīng)振型能量參與系數(shù),并按降序排列,選取結(jié)構(gòu)風(fēng)振響應(yīng)的主要模態(tài);最后,通過選取模態(tài)進(jìn)行頻域計算并與傳統(tǒng)頻域計算方法(完全二次振型疊加,CQC法)對比。結(jié)果表明:結(jié)構(gòu)自振頻率分布密集;動力特性受活塞位置及內(nèi)部氣壓影響較大;背景、共振響應(yīng)振型能量參與系數(shù)能準(zhǔn)確反映各階模態(tài)能量對響應(yīng)總能量的貢獻(xiàn)程度;通過振型能量參與系數(shù)法確定的主要模態(tài)進(jìn)行頻域分析具有計算效率高,計算精度可控的優(yōu)點(diǎn)。

    關(guān)鍵詞:動力特性;背景響應(yīng);共振響應(yīng);脈動風(fēng)荷載;振型能量參與系數(shù);頻域分析

    中圖分類號:TU973+.213;TU33+3

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.13465/j.cnki.jvs.2015.19.012

    Abstract:The dynamic characteristics of a MAN type dry gas storage tank are very complex, they are different from those of ordinary spatial structures. Here, as the basis of its wind-induced response analysis in frequency domain, its dynamic characteristics were analyzed under different working conditions by using the finite element method. Its modes were sorted with order-descending according to the calculation results of mode energy participation factors of background response and resonant response, respectively. Then, the modes with significant contributions to its wind-induced vibration were selected. Finally, the results calculated in frequency domain with a selected modes were compared with those calculated with CQC method. The results showed that the modal frequency distribution of the structure is dense; the location of piston and internal pressure have larger influences on the dynamic characteristics of the gas storage tank; the mode energy participation factors of background response and resonant one can reflect the contributions of different order modes’ energy to the total energy of the response; the analysis in frequency domain with the main modes determined using the mode energy participation factors has advantages of higher calculation efficiency and controllable computation accuracy.

    基金項(xiàng)目:國家自然科學(xué)基金(61004067)資助項(xiàng)目

    收稿日期:2014-08-01修改稿收到日期:2014-09-18

    Dynamic characteristics of a MAN type dry gas storage tank and recognition of main modes with contributions to its wind-induced vibration

    LIZheng-liang1,2,LIUXin-peng1,YANZhi-tao1,2,JIAOHong-wei1,YUDeng-ke1(1. School of Civil Engineering, Chongqing University, Chongqing 400045, China;2. Key Laboratory of New Technology for Construction of Cities in Mountain Area,Ministry of Education, Chongqing University, Chongqing 400045, China)

    Key words:dynamic characteristics; background response; resonant response; fluctuating wind loads; mode energy participation factor; analysis in frequency domain

    曼型干式煤氣柜(下文簡稱新型煤氣柜)是一種儲存可燃?xì)怏w的大型重要建筑物,廣泛用于冶金、石化和市政等行業(yè),其在節(jié)約能源和保護(hù)環(huán)境上發(fā)揮著重大作用。近年來,隨著我國鋼鐵行業(yè)的發(fā)展,國內(nèi)開始大量建造該類煤氣柜。

    新型煤氣柜作為空間結(jié)構(gòu)的代表,與傳統(tǒng)高層建筑結(jié)構(gòu)多側(cè)重抗震性能的研究不同[1],其風(fēng)振響應(yīng)已成為控制結(jié)構(gòu)設(shè)計的關(guān)鍵因素之一。因此對此類結(jié)構(gòu)風(fēng)振響應(yīng)分析成為了亟待解決的重要問題。從目前的研究狀況來看,在模態(tài)分解法的基礎(chǔ)上,利用隨機(jī)振動理論來求解結(jié)構(gòu)在脈動風(fēng)作用下的動力風(fēng)致響應(yīng)即頻域法,是分析空間結(jié)構(gòu)的首選方法。

    然而,與一般高層建筑結(jié)構(gòu)形式不同,新型煤氣柜具有高跨比較小;柜壁較薄,易產(chǎn)生局部振動;內(nèi)部活塞質(zhì)量巨大,對煤氣柜動力特性影響顯著等特點(diǎn)。國內(nèi)對于此類新型結(jié)構(gòu)形式的動力特性分析尚屬空白。因此作為頻域分析法的基礎(chǔ),對其動力特性的精細(xì)化分析尤為重要。

    另外,在頻域內(nèi)計算大跨度空間結(jié)構(gòu)的脈動風(fēng)振響應(yīng)一般采用兩種方法:基于荷載空間分布模式的Ritz向量疊加法[2]和完全二次振型疊加法(CQC)[3-4]。Ritz向量疊加法基本特點(diǎn)是在確定結(jié)構(gòu)主要模態(tài)時,考慮了外荷載的空間分布模式對激發(fā)結(jié)構(gòu)模態(tài)特性的影響,過濾掉與荷載分布向量正交的振型,從而在進(jìn)行頻域計算是大大減少了計算量,提高了計算效率,然而對于時空分布較復(fù)雜的脈動風(fēng)荷載來說,確定荷載分布模式較為困難,限制了Ritz向量疊加法在風(fēng)振響應(yīng)分析中的應(yīng)用。對于傳統(tǒng)的頻域計算方法,即完全二次振型疊加法,該方法自動包含了所有主要模態(tài)的耦合項(xiàng),計算精度較高。工程分析中,往往取前10~20階模態(tài)進(jìn)行計算,但這種做法,對于大型空間結(jié)構(gòu)來說,往往會忽略某些對風(fēng)振響應(yīng)貢獻(xiàn)較大的高階模態(tài)[5-10],因此如何合理考慮高階振型對結(jié)構(gòu)風(fēng)振響應(yīng)的影響是采用頻域分析法計算結(jié)構(gòu)風(fēng)振響應(yīng)的重要問題。對于曼型煤氣柜這種新型結(jié)構(gòu)形式,國內(nèi)對其動力特性及主要貢獻(xiàn)模態(tài)的研究尚屬空白。本文總結(jié)了振型疊加法計算結(jié)構(gòu)風(fēng)振響應(yīng)的公式[11],以及背景響應(yīng)與共振響應(yīng)的計算方法[12-14],并將背景—共振響應(yīng)振型能量參與系數(shù)法引入該類復(fù)雜空間結(jié)構(gòu)中,從能量的角度定量的選取風(fēng)振響應(yīng)中的主要貢獻(xiàn)模態(tài),解決了該類結(jié)構(gòu)風(fēng)振響應(yīng)計算精度的評價問題。

    1煤氣柜動力特性分析

    1.1煤氣柜簡介

    圖1 新型煤氣柜幾何模型圖 Fig.1 The geometry model of new type dry gas tank

    圖2 活塞幾何模型圖 Fig.2 The geometry model of new type dry gas tank piston

    曼型干式煤氣柜由外筒、柜頂、內(nèi)部活塞以及附屬結(jié)構(gòu)(通風(fēng)塔樓、外部電梯、外部樓梯等)組成,本文研究的煤氣柜的儲氣規(guī)模為30萬m3,其中該類煤氣柜的主要外形尺寸為:煤氣柜總高120m,柜體側(cè)壁高108m,外筒半徑66m,活塞總行程96m,柜頂與活塞的拱頂曲率半徑為88m,所有型鋼及鋼板為Q235鋼,見圖1~圖2。

    1.2無內(nèi)壓作用煤氣柜模態(tài)分析

    不考慮柜體內(nèi)外氣體壓力差,利用ANSYS有限元分析軟件分別對煤氣柜在三種工況(活塞高位、中位、低位)下的動力特性進(jìn)行分析。由于煤氣柜是軸對稱結(jié)構(gòu),因此在模態(tài)分析時均將振動形式相同的模態(tài)合并。

    1.2.1柜體在活塞高位狀態(tài)下的模態(tài)和頻率

    煤氣柜在活塞位于高位時(活塞底部距離地面97.2m),通過模態(tài)分析得到結(jié)構(gòu)的頻率和模態(tài)如表1和圖3(由于篇幅限制僅給出典型模態(tài),下同)。

    表1  新型煤氣柜結(jié)構(gòu)的自振頻率

    圖3 煤氣柜振型圖 Fig.3 Vibration chart of new type gas tank

    分析可知:當(dāng)活塞位于高位時,其自振頻率的分布比較密集,前40階模態(tài)頻率分布在1.3~7Hz;柜頂?shù)恼駝幽B(tài)僅出現(xiàn)在第14和25階,這一定程度說明柜頂?shù)膭偠纫黠@大于柜體;結(jié)構(gòu)前25階振型中,一、二階振型為x,y方向的懸臂梁式振動,第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動,其余各階均為殼體的局部振動。

    1.2.2柜體在活塞中位狀態(tài)下的模態(tài)和頻率

    煤氣柜在活塞位于中位時(活塞底部距離地面54m),通過模態(tài)分析得到結(jié)構(gòu)的頻率和振型如表2和圖4。

    表2 新型煤氣柜結(jié)構(gòu)的自振頻率

    圖4 煤氣柜振型圖 Fig.4 Vibration chart of new type gas tank

    分析可知:當(dāng)活塞位于中位時,煤氣柜模態(tài)頻率比高位模態(tài)頻率大,柜體整體剛度增大;與高位時模態(tài)一樣,自振頻率的分布比較密集,前40階振型的頻率分布在2~7Hz;結(jié)構(gòu)前25階振型中,一、二階振型為x,y方向的懸臂梁式振動,第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動,這與活塞高位時的振動形態(tài)相似,但其對應(yīng)的自振頻率與活塞高位差別較大。對于高階模態(tài),盡管兩種工況下(活塞高、中位)均為殼體的局部振動,但相應(yīng)階數(shù)的模態(tài)振動形式及其對應(yīng)自振頻率均變化較大,這是因?yàn)榛钊奈恢米兓瘜Y(jié)構(gòu)的動力特性影響很大。

    1.2.3柜體在活塞低位狀態(tài)下的模態(tài)和頻率

    煤氣柜在活塞低位狀態(tài)下(活塞底部距離地面3.6m),通過模態(tài)分析得到結(jié)構(gòu)的頻率和振型如表3和圖5。

    表3 新型煤氣柜結(jié)構(gòu)的自振頻率

    圖5 煤氣柜振型圖 Fig.5 Vibration chart of new type gas tank

    分析可知:結(jié)構(gòu)模態(tài)頻率比前兩種工況相應(yīng)模態(tài)頻率大;前40階振型的頻率分布在2.7~7Hz,自振頻率的分布相對高位和中位更密集;隨著活塞位置的繼續(xù)下降,煤氣柜前25階振型中,不同于高位和中位振型,低位的一階振型為環(huán)向四個波形的殼體振動,第二、三階振型為懸臂梁式振動,柜體的扭轉(zhuǎn)振動出現(xiàn)在第18階,其余為殼體的局部模態(tài)。

    1.3內(nèi)壓作用下煤氣柜結(jié)構(gòu)模態(tài)分析

    煤氣柜處于工作狀態(tài)時,內(nèi)部封閉有一定壓力的煤氣,30萬m3的煤氣柜設(shè)計內(nèi)壓為15kPa。內(nèi)壓的存在將會影響結(jié)構(gòu)的剛度,并且對結(jié)構(gòu)的動力特性產(chǎn)生影響。故本文主要針對活塞處于高位狀態(tài)時,內(nèi)壓作用與否對煤氣柜模態(tài)分析。

    表4 新型煤氣柜結(jié)構(gòu)的自振頻率

    表4為活塞在高位狀態(tài)下,考慮內(nèi)壓與未考慮內(nèi)壓煤氣柜結(jié)構(gòu)的頻率與振型特點(diǎn)描述(由于振型相似及篇幅原因內(nèi)壓作用下的振型截圖不再給出)。從表4可以看出,考慮內(nèi)壓對煤氣柜的平動振型和扭轉(zhuǎn)振型幾乎沒有影響,有內(nèi)壓和無內(nèi)壓的平動頻率和扭轉(zhuǎn)頻率均相同;而內(nèi)壓對柜體殼體局部振動影響較顯著,考慮內(nèi)壓的影響,殼體振動的頻率要比不考慮內(nèi)壓最大增加12%左右。由此可知,柜體內(nèi)部氣壓對其局部剛度影響顯著,而對整體剛度幾乎沒有影響。

    2背景-共振響應(yīng)振型能量參與系數(shù)

    2.1脈動風(fēng)振響應(yīng)振型疊加原理

    對于新型煤氣柜結(jié)構(gòu)在外施動力荷載作用下的動力方程為:

    (1)

    結(jié)構(gòu)振型矩陣Φ=[φ1,φ2,…,φn],且ΦTM Φ=I(其中I為單位矩陣),則ΦTK Φ=Ω2(其中Ω為結(jié)構(gòu)前n階圓頻率組成的對角矩陣,即diag[ω1,ω2,…,ωn]結(jié)構(gòu)在脈動風(fēng)作用下的位移響應(yīng)方差矩陣為[15]:

    LTΦH*T(ω)ΦTdω

    (2)

    式中,SPP(ω)測壓點(diǎn)脈動風(fēng)壓時程互譜矩陣;H(ω)為頻響函數(shù)矩陣,且H(ω)=[H1(ω),H2(ω),…,Hn(ω)],其中,Hj(ω)為第j階模態(tài)的頻響函數(shù)。

    2.2背景響應(yīng)及其能量參與系數(shù)的計算

    將脈動風(fēng)作用下的結(jié)構(gòu)振動看作靜力過程,即不考慮結(jié)構(gòu)的動力效應(yīng),可得到結(jié)構(gòu)總的背景響應(yīng)及各振型的背景響應(yīng)。

    (1)總背景響應(yīng)

    脈動風(fēng)作用下,總的背景響應(yīng)Xb可表示為:

    KXb=LPd(t)

    (3)

    Xb的方差可表示為:

    (4)

    式中,XP為靜力平衡方程KXb=L的解;b為背景響應(yīng)。

    (2)振型背景響應(yīng)

    設(shè)Xb=ΦQb,將其代入式(3),并對等式左右同時左乘ΦT,那么式(3)可改寫為:

    Ω2Qb=ΦTLPd(t)

    (5)

    得到:

    Qb=H0ΦTLPd(t)

    (6)

    則第j階模態(tài)背景響應(yīng)為:

    qb,j=Hj,0φTjLPd(t)

    (7)

    第j階模態(tài)對總背景響應(yīng)的貢獻(xiàn)為:

    Xb,j=φjHj,0φTjLPd(t)

    (8)

    振型位移Qb的方差為:

    (9)

    式中,第j階模態(tài)背景響應(yīng)的方差為:

    (10)

    (3)背景響應(yīng)振型能量參與系數(shù)

    脈動風(fēng)荷載在背景響應(yīng)上做的功為:

    (11)

    式中,∑diag[]表示矩陣的對角元素之和。

    對背景響應(yīng)功Wb取數(shù)學(xué)期望,得到

    (12)

    脈動風(fēng)荷載在第j階模態(tài)背景響應(yīng)上做的功為:

    (13)

    對在第j階模態(tài)背景響應(yīng)功Wb,j取數(shù)學(xué)期望,得到:

    (14)

    第j階振型背景響應(yīng)的能量參與系數(shù)是:

    (15)

    由式(15)可知,背景響應(yīng)的振型能量參與系數(shù)反映了振型位移對總背景響應(yīng)位移的平均貢獻(xiàn)率。背景響應(yīng)位移計算值的精度可通過振型累積參與系數(shù)以及振型參與系數(shù)來控制。

    2.4共振響應(yīng)及其能量參與系數(shù)的計算

    共振響應(yīng)是指脈動風(fēng)荷載在共振頻率處激起的結(jié)構(gòu)振動,對第j階振型的運(yùn)動方程兩側(cè)進(jìn)行傅里葉變換,得到第j階振型響應(yīng)的頻域解為:

    qj(ω)=Hj(ω)φTjLPd(ω)

    (16)

    第j階振型的共振響應(yīng)為:

    qr,j(ω)=Hj(ω)φTjLPd(ωj)

    (17)

    式中,下標(biāo)“r”表示共振響應(yīng);Pd(ωj)為脈動風(fēng)壓傅里葉譜值在第j階振型頻率ωj處的取值。

    第j階振型共振響應(yīng)的位移方差為:

    (18)

    式中,下標(biāo)“d”表示位移響應(yīng);ξj為第j階模態(tài)阻尼比;SPP(ωj)為脈動風(fēng)壓功率譜在第j階模態(tài)頻率ωj處的取值。

    第j階振型共振響應(yīng)的速度方差為:

    (19)

    式中,下標(biāo)“v”表示位移響應(yīng)。

    脈動風(fēng)荷載在第j階振型共振響應(yīng)上的功等于第j階振型的動能和彈性能之和:

    (20)

    對Wr,j取數(shù)學(xué)期望,得到:

    (21)

    共振響應(yīng)總能量,從理論上來說等于全部模態(tài)共振響應(yīng)的能量之和,但當(dāng)結(jié)構(gòu)自由度十分龐大時,共振總能量很難得到。因此,為了方便起見,將共振能量參與系數(shù)定義為各階共振響應(yīng)能量與背景響應(yīng)總能量之比,如下式:

    (22)

    3煤氣柜主要貢獻(xiàn)模態(tài)選取

    新型煤氣柜為軸對稱結(jié)構(gòu)體系,因此對其進(jìn)行風(fēng)振響應(yīng)頻域分析時僅考慮0°風(fēng)向角情況,本文運(yùn)用線性濾波AR模型法[16]對脈動風(fēng)速進(jìn)行模擬,模擬結(jié)果見圖6~圖7。

    圖6 脈動風(fēng)時程曲線 Fig.6 Turbulence wind velocity time-history curve

    圖7 模擬風(fēng)譜與目標(biāo)譜的比較 Fig.7 The target spectrum and simulation spectrum

    3.1煤氣柜背景響應(yīng)主要貢獻(xiàn)模態(tài)

    由于煤氣柜內(nèi)部活塞上下移動而造成結(jié)構(gòu)體系動力特性的時變性,因此本文首先給出活塞三個典型位置即高位、中位、低位前100階背景響應(yīng)能量參與系數(shù),見圖8。前文分析可知,煤氣柜內(nèi)部氣壓對結(jié)構(gòu)模態(tài)影響較大,因此以下工況分析均與實(shí)際情況相同即考慮內(nèi)部氣壓對結(jié)構(gòu)的影響。

    接下來分別對三種工況的背景響應(yīng)能量參與系數(shù)按照降序排列,并進(jìn)行累加計算。當(dāng)背景響應(yīng)振型能量參與系數(shù)<2%,并且累積系數(shù)使之>90%時,對背景響應(yīng)主要貢獻(xiàn)模態(tài)選取結(jié)束,篩選結(jié)果見圖9。

    圖8 背景響應(yīng)能量參與系數(shù) Fig.8 The mode energy participation factorof background response

    圖9 背景響應(yīng)主要貢獻(xiàn)模態(tài) Fig.9 The natural modes with significant contribution to background response

    3.2煤氣柜共振響應(yīng)主要貢獻(xiàn)模態(tài)

    分別給出三種工況下的共振能量參與系數(shù),見圖10。與背景響應(yīng)主要模態(tài)選取原則不同,對共振響應(yīng)主要貢獻(xiàn)模態(tài)選取時,當(dāng)共振能量參與系數(shù)<2%時即可結(jié)束篩選,這主要是因?yàn)樵谟嬎愎舱衲芰繀⑴c系數(shù)時,所用歸一化系數(shù)為背景響應(yīng)總能量值,這將造成共振累積能量很難達(dá)到90%以上,篩選結(jié)果見圖11。

    需要注意的是,本文風(fēng)荷載的采樣頻率為16Hz,所以在對共振主要模態(tài)的選取時,僅能考慮頻率低于16Hz的模態(tài),對于頻率大于16Hz的模態(tài),隨著模態(tài)階數(shù)的增加,其頻率越來越大,相對低頻的荷載,越來越難以激起高階模態(tài)的共振,因此高階模態(tài)共振響應(yīng)越來越小,可以忽略。

    圖10 共振響應(yīng)能量參與系數(shù) Fig.10 The mode energy participation factor of resonant response

    圖11 共振響應(yīng)主要貢獻(xiàn)模態(tài) Fig.11 The natural modes with significant contribution to resonant response

    對于共振響應(yīng)能量,由圖11分析可知,篩選出的主要模態(tài)對共振響應(yīng)總能量的貢獻(xiàn)都不大,這主要是兩方面造成的:結(jié)構(gòu)模態(tài)分布密集,脈動風(fēng)壓能夠激起多階模態(tài)的振動,這可以理解為多階模態(tài)均攤了共振響應(yīng)總能量;由于結(jié)構(gòu)局部模態(tài)剛度較大,各局部模態(tài)自振頻率與脈動風(fēng)荷載頻率相差較遠(yuǎn),使得脈動風(fēng)荷載無法激起某階模態(tài)的較大共振。

    為了驗(yàn)證上述模態(tài)選取方法的準(zhǔn)確性,本文首先采用完全二次振型疊加法(CQC法)按式(2)分別計算三種工況下煤氣柜在脈動風(fēng)作用下的總響應(yīng),并將前100階模態(tài)計算得到的位移均方根作為準(zhǔn)確值。接著按上述方法識別煤氣柜主要貢獻(xiàn)模態(tài),并通過其進(jìn)行頻域分析。最后,給出按上述兩種方法計算得到的結(jié)構(gòu)典型位置測點(diǎn)(高度100.8m)的脈動響應(yīng)均方根值,見圖12。

    圖12 CQC法與振型能量參與系數(shù)法的結(jié)果 Fig.12 The results of CQC method and mode energy participation factor method

    圖13 振型能量參與系數(shù)法的計算效率 Fig.13 The efficiency ofmode energy participation factor method

    結(jié)果表明,三種工況下,分別通過主要模態(tài)與前100階模態(tài)進(jìn)行頻域分析的結(jié)果吻合很好,證明主要貢獻(xiàn)模態(tài)選取合理,準(zhǔn)確。

    為了驗(yàn)證通過主要貢獻(xiàn)模態(tài)進(jìn)行結(jié)構(gòu)頻域分析的計算效率,本文給出了煤氣柜6個典型測點(diǎn)在三種工況下參與計算的主要貢獻(xiàn)模態(tài)數(shù)不同,其分析結(jié)果與標(biāo)準(zhǔn)值的比值,見圖13。

    結(jié)果表明,隨著主要貢獻(xiàn)模態(tài)數(shù)的不斷累加,分析結(jié)果迅速收斂,這比傳統(tǒng)的頻域分析方法效率更高,且計算精度可以通過調(diào)整篩選原則控制。對于這種自振頻率密集,參振模態(tài)多的結(jié)構(gòu)體系,振型能量參與系數(shù)能夠很好的識別主要模態(tài)。

    4結(jié)論

    (1)活塞位置對新型煤氣柜動力特性的影響顯著,隨著活塞位置降低,煤氣柜自振頻率逐漸變大。

    (2)煤氣柜在活塞高位和中位狀態(tài)下,一、二階振型為x,y方向的懸臂梁式振動,第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動;不同于高位和中位振型,低位的一階振型為環(huán)向四個波形的殼體振動,第二、三階振型為懸臂梁式振動,柜體的扭轉(zhuǎn)振動出現(xiàn)在第18階,其余各階局部振動。

    (3)內(nèi)壓對煤氣柜的平動振型和扭轉(zhuǎn)振型幾乎沒有影響,而內(nèi)壓對柜體殼體振動影響較顯著,考慮內(nèi)壓的影響,殼體振動的頻率要比不考慮內(nèi)壓最大增加12%左右。

    (4)脈動風(fēng)振響應(yīng)的計算精度與模態(tài)選取密切相關(guān),模態(tài)的選取取決于各階模態(tài)響應(yīng)對總響應(yīng)的貢獻(xiàn)大小。背景能量參與系數(shù)法能夠準(zhǔn)確識別出煤氣柜脈動風(fēng)振響應(yīng)的主要貢獻(xiàn)模態(tài)。

    (5)通過調(diào)整篩選原則,可以對煤氣柜頻域分析計算精度進(jìn)行控制,進(jìn)而能有效提高頻域分析效率。

    參考文獻(xiàn)

    [1]沈世釗.大跨空間結(jié)構(gòu)的發(fā)展—回顧與展望[J].土木工程學(xué)報,1998,31(3):5-14.

    SHEN Shi-zhao. Development Of Long-Span Structures—A Review And Prospect[J]. China Civil Engineering Journal,1998,31(3):5-14.

    [2]Wilson E L, YUAN Ming-wu. Dynamic analysis by direct superposition of Ritz vectors [J]. Earthquake Engineering and Structural Dynamics, 1982(10):813-821.

    [3]顧明,黃鵬,周毅,等.北京首都機(jī)場3號航站樓風(fēng)荷載和響應(yīng)研究[J].土木工程學(xué)報,2005,38(1): 40-44.

    GU Ming, HUANG Peng, ZHOU Yi,et al. A study on wind loads and responses of terminal 3 at beijing capital airport [J]. China Civil Engineering Journal, 2005,38(1):40-44.

    [4]王國硯,黃本才,林穎儒,等.基于CQC方法的大跨屋蓋結(jié)構(gòu)隨機(jī)風(fēng)振響應(yīng)計算[J].空間結(jié)構(gòu),2003,9(4):22-26.

    WANG Guo-yan, HUANG Ben-cai, LIN Ying-ru. Solution of wind induced random vibration of large span roof structures based on CQC method[J]. Spatial Structures,2003,9(4):22-26.

    [5]Nakayama M, Sasaky Y, Masuda K, et al. An efficient method for selection of vibration modes contributory to wind response on dome like roofs[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1998, 73(1):31-43.

    [6]Wu Y, Zhang J S, Chen B, et al. Identification of the dominant vibration modes of single-layer reticulated shells under wind action[J]. International Journal of Space Structures, 2007, 22(2):123- 132.

    [7]何艷麗,董石麟,龔景海. 空間網(wǎng)格結(jié)構(gòu)頻域風(fēng)振響應(yīng)分析模態(tài)補(bǔ)償法[J].工程力學(xué),2002,19(4):1-6.

    HE Yan-li, DONG Shi-lin, GONG Jing-hai. Wind-induced response of spatial structures with mode compensation in frequency domain[J].Engineering Mechanics,2002,19(4):1-6.

    [8]胡繼軍,網(wǎng)殼風(fēng)振及控制研究[D].上海:上海交通大學(xué),2000.

    [9]Mataki Y, Iwasa Y, Fukao Y, et al. Wind induced response of low-profile cable-reinforced air-supported structures [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1988, 29:253-262.

    [10]鄧華,董石麟,何艷麗,等.深圳游泳跳水館主館屋蓋結(jié)構(gòu)分析及風(fēng)振響應(yīng)計算[J].建筑結(jié)構(gòu)學(xué)報, 2004, 25(2):72-78.

    DENG Hua, DONG Shi-lin, HE Yan-li, et al. Behavior analysis and wind-induced vibration computationof mail hall roof structure of shenzhen swimming and diving hall[J]. Journal of Building Structures,2004,25(2):72-78.

    [11]田玉基,楊慶山,范重,等.國家體育場大跨度屋蓋結(jié)構(gòu)風(fēng)振系數(shù)研究 [J]. 建筑結(jié)構(gòu)學(xué)報, 2007, 28(2): 36-41.

    TIAN Yu-ji, YANG Qing-shan, FAN Zhong, et al. Gust factors for large-span roof of the national stadium[J]. Journal of Building Structures,2007,28(2):36-41.

    [12]Ashraf Ali M, Gould P L. On the resonant component of the response of single degree of freedom systems under wind loading [J]. Engineering Structures, 1985, 7(4): 280-292.

    [13]陳波, 武岳, 沈世釗. 背景響應(yīng)、共振響應(yīng)定義及其相關(guān)性分析方法 [J]. 振動工程學(xué)報, 2008, 21(2):140-145.

    CHEN Bo, WU Yue, SHEN Shi-zhao. Definitions and correlation analyses for background response and resonant response[J]. Journal of VibrationEngineering, 2008, 21(2):140-145.

    [14]田玉基,楊慶山. 大跨度屋蓋結(jié)構(gòu)脈動風(fēng)振響應(yīng)的參與振型[J]. 哈爾濱工業(yè)大學(xué)學(xué)報,2009(6):146-149.

    TIAN Yu-ji, YANG Qing-shan. Dominant modes of fluctuating wind-induced response for large-span roof [J]. Journal of Harbin Institute of Technology,2009(6):146-149.

    [15]方同.工程隨機(jī)振動[M].北京:國防工業(yè)出版社,1995.

    [16]舒新玲, 周岱.風(fēng)速時程AR模型及其快速實(shí)現(xiàn)[J].空間結(jié)構(gòu), 2003, 12(9): 27-32.

    SHU Xin-ling, ZHOU Dai. AR model of wind speed timeseries and its rapid implementation[J].Spatial Structures,2003,12(9):27-32.

    第一作者劉均男,副教授,1976年8月生

    猜你喜歡
    動力特性
    基于ANSYS的直立塔設(shè)備動力特性研究
    商情(2017年22期)2017-07-18 10:10:39
    在役鋼筋混凝土雙曲拱橋動力測試和有限元分析
    科技視界(2017年6期)2017-07-01 10:56:52
    大跨度軌道懸索橋健康監(jiān)測問題思考與對策
    魅力中國(2017年4期)2017-05-13 17:54:57
    AP1100汽機(jī)基礎(chǔ)動力特性研究
    單索面公軌兩用鋼桁梁獨(dú)塔斜拉橋施工動力特性分析
    斜拉—懸吊協(xié)作體系橋動力響應(yīng)分析
    干濕循環(huán)下水泥摻入對膨脹土動力特性影響的試驗(yàn)分析
    自升式海洋平臺動力特性分析
    高層建筑氣動彈性模型風(fēng)洞試驗(yàn)研究
    浮置板軌道結(jié)構(gòu)板下支承剛度研究
    亚洲成av人片免费观看| 国产激情欧美一区二区| 麻豆成人av在线观看| 日韩免费av在线播放| 亚洲片人在线观看| 久久亚洲真实| 亚洲男人的天堂狠狠| 在线国产一区二区在线| 日韩欧美国产一区二区入口| 亚洲精品色激情综合| 性色av乱码一区二区三区2| 在线观看免费午夜福利视频| 成人一区二区视频在线观看| 亚洲精品在线美女| 国产欧美日韩一区二区三| 国产成人精品久久二区二区91| 桃色一区二区三区在线观看| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清 | 日韩欧美国产在线观看| 国产黄片美女视频| 精品国产美女av久久久久小说| 成人午夜高清在线视频| 可以免费在线观看a视频的电影网站| 黑人操中国人逼视频| 国内少妇人妻偷人精品xxx网站 | 一级作爱视频免费观看| 日日摸夜夜添夜夜添小说| 精品少妇一区二区三区视频日本电影| 国产精品久久久久久精品电影| 国产高清有码在线观看视频 | 两个人的视频大全免费| 哪里可以看免费的av片| 真人一进一出gif抽搐免费| 久久午夜综合久久蜜桃| 久久精品影院6| 禁无遮挡网站| 看免费av毛片| 国产三级中文精品| 免费观看人在逋| 国产精品美女特级片免费视频播放器 | 亚洲片人在线观看| 久久久久九九精品影院| 国内精品久久久久精免费| 欧美三级亚洲精品| 免费高清视频大片| 国产麻豆成人av免费视频| 免费在线观看日本一区| 大型黄色视频在线免费观看| 叶爱在线成人免费视频播放| aaaaa片日本免费| 国产精品野战在线观看| 亚洲精品一区av在线观看| 亚洲狠狠婷婷综合久久图片| 国产区一区二久久| 亚洲一区二区三区不卡视频| 国产又色又爽无遮挡免费看| 日本 欧美在线| 亚洲国产中文字幕在线视频| 老熟妇乱子伦视频在线观看| 免费av毛片视频| cao死你这个sao货| 中文在线观看免费www的网站 | 麻豆国产av国片精品| 国产精品一区二区免费欧美| 日本黄大片高清| 在线观看一区二区三区| 亚洲国产欧美一区二区综合| 国产片内射在线| 国产精品99久久99久久久不卡| 久久精品国产亚洲av香蕉五月| 亚洲av熟女| 久久午夜综合久久蜜桃| 一进一出抽搐gif免费好疼| 免费高清视频大片| 桃色一区二区三区在线观看| 国产精品国产高清国产av| 国产成人系列免费观看| 国产精品 国内视频| 久久久国产成人精品二区| 一本大道久久a久久精品| 亚洲av成人不卡在线观看播放网| 男女之事视频高清在线观看| 91麻豆av在线| 中文亚洲av片在线观看爽| 久久中文看片网| 日韩高清综合在线| 两个人看的免费小视频| 免费在线观看影片大全网站| 久久精品综合一区二区三区| 夜夜夜夜夜久久久久| 日韩欧美在线乱码| 美女免费视频网站| 亚洲精华国产精华精| 天堂影院成人在线观看| 国产爱豆传媒在线观看 | 不卡av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 叶爱在线成人免费视频播放| 国产av一区在线观看免费| 国产精品99久久99久久久不卡| 国产av一区二区精品久久| 免费看十八禁软件| 国内久久婷婷六月综合欲色啪| 精品少妇一区二区三区视频日本电影| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清专用| 麻豆一二三区av精品| 神马国产精品三级电影在线观看 | 精品福利观看| 五月玫瑰六月丁香| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 亚洲欧美日韩无卡精品| 精品免费久久久久久久清纯| 久久国产精品人妻蜜桃| 18禁观看日本| 亚洲精品中文字幕一二三四区| 精品久久久久久成人av| 国产精品国产高清国产av| 久久精品aⅴ一区二区三区四区| 久久久久国内视频| 国内精品久久久久精免费| 欧美+亚洲+日韩+国产| 亚洲中文av在线| 亚洲欧美激情综合另类| 午夜激情av网站| 国产精品一及| 国产麻豆成人av免费视频| 久久精品国产综合久久久| 日本黄大片高清| 色综合亚洲欧美另类图片| 国产成人影院久久av| 亚洲人成网站高清观看| 人人妻,人人澡人人爽秒播| 午夜亚洲福利在线播放| www.熟女人妻精品国产| 12—13女人毛片做爰片一| 舔av片在线| 欧美成人免费av一区二区三区| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 最近在线观看免费完整版| 黄色毛片三级朝国网站| 一区福利在线观看| 日韩有码中文字幕| 国产精品一区二区精品视频观看| 欧美午夜高清在线| 一级作爱视频免费观看| 亚洲色图av天堂| 免费在线观看完整版高清| 国产麻豆成人av免费视频| 欧美性长视频在线观看| 男女床上黄色一级片免费看| 国产单亲对白刺激| 小说图片视频综合网站| 国产亚洲欧美在线一区二区| svipshipincom国产片| 不卡一级毛片| 成在线人永久免费视频| 亚洲国产精品久久男人天堂| 久久久久亚洲av毛片大全| 成人午夜高清在线视频| 国产精品电影一区二区三区| www.自偷自拍.com| 母亲3免费完整高清在线观看| 成人亚洲精品av一区二区| 午夜两性在线视频| 国内毛片毛片毛片毛片毛片| 中出人妻视频一区二区| 国产免费男女视频| 国产主播在线观看一区二区| 日韩大码丰满熟妇| 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| cao死你这个sao货| 欧美高清成人免费视频www| 国产三级在线视频| 一本一本综合久久| 免费无遮挡裸体视频| 中文字幕人妻丝袜一区二区| av片东京热男人的天堂| 亚洲国产高清在线一区二区三| 精品久久久久久久毛片微露脸| 午夜福利成人在线免费观看| 好男人电影高清在线观看| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 日本五十路高清| 夜夜爽天天搞| 色综合婷婷激情| 男女做爰动态图高潮gif福利片| 久久九九热精品免费| 精品高清国产在线一区| 日韩欧美三级三区| 久久久久久亚洲精品国产蜜桃av| 国产视频一区二区在线看| 女生性感内裤真人,穿戴方法视频| 老司机午夜十八禁免费视频| 色播亚洲综合网| 天天躁狠狠躁夜夜躁狠狠躁| 叶爱在线成人免费视频播放| 国产熟女午夜一区二区三区| 可以在线观看的亚洲视频| 国产精品一区二区免费欧美| 人妻夜夜爽99麻豆av| 亚洲无线在线观看| 亚洲av成人精品一区久久| av视频在线观看入口| aaaaa片日本免费| 中文字幕高清在线视频| 午夜影院日韩av| 夜夜爽天天搞| 国内精品久久久久精免费| 一进一出抽搐gif免费好疼| 午夜a级毛片| a级毛片在线看网站| 日本三级黄在线观看| 国产伦一二天堂av在线观看| 精品人妻1区二区| 亚洲性夜色夜夜综合| 夜夜夜夜夜久久久久| www日本在线高清视频| 巨乳人妻的诱惑在线观看| 亚洲一码二码三码区别大吗| 精品久久久久久成人av| 亚洲电影在线观看av| 久久久国产精品麻豆| 国产视频内射| 国产熟女xx| 91字幕亚洲| 亚洲中文av在线| 成熟少妇高潮喷水视频| 欧美 亚洲 国产 日韩一| 久久婷婷人人爽人人干人人爱| 色精品久久人妻99蜜桃| 欧美成人免费av一区二区三区| 黄色成人免费大全| 99热只有精品国产| 国产精品 国内视频| 久久草成人影院| 午夜日韩欧美国产| 亚洲精品在线观看二区| 99在线人妻在线中文字幕| 又大又爽又粗| 成人18禁在线播放| 757午夜福利合集在线观看| 久久久国产欧美日韩av| 男女之事视频高清在线观看| 婷婷丁香在线五月| 国产精品亚洲美女久久久| 亚洲国产精品999在线| 岛国在线免费视频观看| 在线观看日韩欧美| 亚洲中文日韩欧美视频| 日本免费一区二区三区高清不卡| 少妇被粗大的猛进出69影院| 99国产综合亚洲精品| 国产单亲对白刺激| 757午夜福利合集在线观看| 日本 av在线| 毛片女人毛片| 亚洲av片天天在线观看| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 一级毛片精品| 无遮挡黄片免费观看| 19禁男女啪啪无遮挡网站| 成人av一区二区三区在线看| 88av欧美| 欧美人与性动交α欧美精品济南到| 久久久久精品国产欧美久久久| 全区人妻精品视频| 国产av不卡久久| 麻豆成人午夜福利视频| 亚洲 欧美一区二区三区| 999久久久国产精品视频| 久久久久久大精品| 在线看三级毛片| 亚洲专区字幕在线| 精华霜和精华液先用哪个| 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清| 日韩欧美国产一区二区入口| 日本精品一区二区三区蜜桃| 亚洲人与动物交配视频| 91九色精品人成在线观看| 午夜福利欧美成人| 18美女黄网站色大片免费观看| 正在播放国产对白刺激| 一区福利在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日本五十路高清| 午夜精品久久久久久毛片777| 日本成人三级电影网站| 亚洲欧美一区二区三区黑人| 成人永久免费在线观看视频| av在线天堂中文字幕| 淫妇啪啪啪对白视频| 一区福利在线观看| 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区三区视频了| 亚洲精华国产精华精| 国产精品乱码一区二三区的特点| 香蕉丝袜av| 成人欧美大片| 19禁男女啪啪无遮挡网站| xxx96com| 久久中文看片网| 亚洲国产精品久久男人天堂| 一区福利在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产又色又爽无遮挡免费看| 校园春色视频在线观看| av欧美777| 久久精品国产99精品国产亚洲性色| 国产爱豆传媒在线观看 | 白带黄色成豆腐渣| 夜夜爽天天搞| 免费看a级黄色片| 国产成人aa在线观看| 999久久久精品免费观看国产| 久久精品91无色码中文字幕| 欧美乱色亚洲激情| 黄频高清免费视频| 欧美精品啪啪一区二区三区| 岛国在线免费视频观看| 激情在线观看视频在线高清| 人成视频在线观看免费观看| 久久久久国产精品人妻aⅴ院| 黄色视频,在线免费观看| 真人做人爱边吃奶动态| videosex国产| 男男h啪啪无遮挡| 日韩三级视频一区二区三区| 啦啦啦韩国在线观看视频| 日本精品一区二区三区蜜桃| 精品国产乱子伦一区二区三区| 亚洲成人国产一区在线观看| 国产精品一区二区免费欧美| 十八禁网站免费在线| 国产激情久久老熟女| 18禁黄网站禁片免费观看直播| aaaaa片日本免费| 精品电影一区二区在线| 一本一本综合久久| 一区二区三区高清视频在线| 亚洲一卡2卡3卡4卡5卡精品中文| 精品电影一区二区在线| 精品无人区乱码1区二区| 又粗又爽又猛毛片免费看| 好男人在线观看高清免费视频| 在线a可以看的网站| av在线天堂中文字幕| netflix在线观看网站| 一个人免费在线观看的高清视频| netflix在线观看网站| 91成年电影在线观看| 国产av麻豆久久久久久久| 亚洲精品美女久久久久99蜜臀| www.熟女人妻精品国产| 免费电影在线观看免费观看| xxx96com| 变态另类丝袜制服| 国产成人精品久久二区二区免费| 90打野战视频偷拍视频| 国产午夜精品论理片| 午夜福利免费观看在线| 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 真人做人爱边吃奶动态| 亚洲性夜色夜夜综合| 搡老熟女国产l中国老女人| 舔av片在线| 国产爱豆传媒在线观看 | 精品电影一区二区在线| 99热6这里只有精品| 亚洲精品在线观看二区| 欧美中文日本在线观看视频| 一a级毛片在线观看| 色av中文字幕| 国产高清激情床上av| 欧洲精品卡2卡3卡4卡5卡区| 18禁国产床啪视频网站| 国内精品久久久久精免费| 天堂√8在线中文| 欧美一级a爱片免费观看看 | 日韩精品中文字幕看吧| 岛国在线免费视频观看| 又紧又爽又黄一区二区| 久久久久久久精品吃奶| 精品少妇一区二区三区视频日本电影| 亚洲成av人片免费观看| 日本黄大片高清| 午夜精品久久久久久毛片777| ponron亚洲| 三级毛片av免费| 亚洲av成人不卡在线观看播放网| 女人被狂操c到高潮| 亚洲国产精品sss在线观看| 精品午夜福利视频在线观看一区| 色噜噜av男人的天堂激情| 美女扒开内裤让男人捅视频| 丰满的人妻完整版| www.熟女人妻精品国产| 国产高清视频在线观看网站| 国产1区2区3区精品| 久久香蕉国产精品| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| www日本在线高清视频| 夜夜爽天天搞| 亚洲精品在线美女| 三级毛片av免费| 免费看a级黄色片| 老司机午夜福利在线观看视频| 国产不卡一卡二| 在线播放国产精品三级| 波多野结衣高清作品| 久久 成人 亚洲| 久久久精品欧美日韩精品| 又黄又粗又硬又大视频| 男女视频在线观看网站免费 | 99热6这里只有精品| 色在线成人网| 国产精品 欧美亚洲| 宅男免费午夜| 日日爽夜夜爽网站| 91av网站免费观看| 欧美日韩瑟瑟在线播放| 嫁个100分男人电影在线观看| 黄色 视频免费看| bbb黄色大片| 女人高潮潮喷娇喘18禁视频| 精品久久久久久成人av| 深夜精品福利| 9191精品国产免费久久| 听说在线观看完整版免费高清| 国产伦一二天堂av在线观看| 亚洲五月天丁香| 亚洲成人免费电影在线观看| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 久久人妻av系列| 老汉色av国产亚洲站长工具| 搡老熟女国产l中国老女人| 大型av网站在线播放| 国产视频一区二区在线看| 99久久综合精品五月天人人| 亚洲国产看品久久| 观看免费一级毛片| 亚洲精品粉嫩美女一区| 国产亚洲精品一区二区www| 男女之事视频高清在线观看| 在线观看免费视频日本深夜| 女警被强在线播放| 少妇人妻一区二区三区视频| 国产久久久一区二区三区| 丰满人妻一区二区三区视频av | 国产av又大| 久久欧美精品欧美久久欧美| 熟女少妇亚洲综合色aaa.| 亚洲,欧美精品.| 丁香六月欧美| 51午夜福利影视在线观看| 黄色成人免费大全| 操出白浆在线播放| 99re在线观看精品视频| 亚洲在线自拍视频| 久久中文字幕一级| 成人永久免费在线观看视频| 免费看十八禁软件| 欧美+亚洲+日韩+国产| 欧美一级毛片孕妇| 在线免费观看的www视频| 国产激情偷乱视频一区二区| 少妇裸体淫交视频免费看高清 | 欧美日韩亚洲综合一区二区三区_| 757午夜福利合集在线观看| 久久久久九九精品影院| 亚洲精品在线观看二区| √禁漫天堂资源中文www| www.999成人在线观看| 亚洲中文日韩欧美视频| 99久久无色码亚洲精品果冻| 日韩欧美 国产精品| 日本撒尿小便嘘嘘汇集6| 一a级毛片在线观看| 亚洲在线自拍视频| 香蕉丝袜av| 黑人巨大精品欧美一区二区mp4| 最近在线观看免费完整版| 精品人妻1区二区| 国产精品 欧美亚洲| 九色成人免费人妻av| 亚洲一区二区三区不卡视频| 日韩高清综合在线| 午夜福利在线在线| 成人18禁在线播放| 亚洲av日韩精品久久久久久密| 色老头精品视频在线观看| 桃红色精品国产亚洲av| av天堂在线播放| 美女扒开内裤让男人捅视频| 国产黄a三级三级三级人| 一本大道久久a久久精品| 丝袜人妻中文字幕| 日本五十路高清| 在线播放国产精品三级| 91成年电影在线观看| 欧美激情久久久久久爽电影| 午夜精品在线福利| 国产野战对白在线观看| 中文在线观看免费www的网站 | 757午夜福利合集在线观看| 91国产中文字幕| av中文乱码字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品日韩av在线免费观看| 欧美日韩瑟瑟在线播放| 久久香蕉激情| 在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 亚洲欧美一区二区三区黑人| 亚洲 欧美 日韩 在线 免费| 日本撒尿小便嘘嘘汇集6| 他把我摸到了高潮在线观看| 9191精品国产免费久久| 日韩高清综合在线| 国产男靠女视频免费网站| 黄色丝袜av网址大全| 一个人免费在线观看的高清视频| 变态另类成人亚洲欧美熟女| 免费在线观看亚洲国产| 国产精品久久久人人做人人爽| 波多野结衣巨乳人妻| 国产精品 国内视频| 久久精品人妻少妇| 亚洲国产中文字幕在线视频| 亚洲av日韩精品久久久久久密| 国语自产精品视频在线第100页| 国产精品美女特级片免费视频播放器 | 亚洲在线自拍视频| 日韩欧美免费精品| 久久精品91蜜桃| 国产成人欧美在线观看| 国产精品一区二区三区四区久久| 又黄又粗又硬又大视频| 成年人黄色毛片网站| 亚洲成人中文字幕在线播放| 男女那种视频在线观看| 国产一区二区激情短视频| 国产精品一区二区三区四区免费观看 | 99久久久亚洲精品蜜臀av| 国产成人欧美在线观看| 亚洲熟妇中文字幕五十中出| 国产熟女xx| 久久九九热精品免费| 18禁黄网站禁片免费观看直播| 99热这里只有是精品50| 制服丝袜大香蕉在线| 午夜两性在线视频| 久久草成人影院| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 亚洲成人久久爱视频| 搡老岳熟女国产| 亚洲精品av麻豆狂野| 亚洲五月天丁香| 女人高潮潮喷娇喘18禁视频| 国产精品一区二区免费欧美| 一个人免费在线观看电影 | 老司机靠b影院| 很黄的视频免费| 精华霜和精华液先用哪个| 欧美日韩亚洲国产一区二区在线观看| 国产高清激情床上av| www.精华液| 久久精品国产清高在天天线| 精品国内亚洲2022精品成人| 亚洲无线在线观看| 日韩国内少妇激情av| 国产精品久久视频播放| 不卡一级毛片| 亚洲欧美日韩高清专用| 国产三级中文精品| 国产97色在线日韩免费| 午夜福利欧美成人| 国产又黄又爽又无遮挡在线| 国产精品,欧美在线| 亚洲欧美日韩高清在线视频| 国产一区二区激情短视频| 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 免费在线观看视频国产中文字幕亚洲| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 国产黄a三级三级三级人| 无限看片的www在线观看| 免费搜索国产男女视频| 亚洲专区国产一区二区| 好男人在线观看高清免费视频| 国产爱豆传媒在线观看 | 国产99白浆流出| x7x7x7水蜜桃| 美女午夜性视频免费| 在线播放国产精品三级| 天堂√8在线中文| 国产伦一二天堂av在线观看| cao死你这个sao货| 国产精品久久久久久久电影 | 色在线成人网|