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

    基于精細(xì)積分法的結(jié)構(gòu)碰撞反應(yīng)譜研究

    2016-01-18 03:39:09尹俊紅,李青寧,張瑞杰
    振動與沖擊 2015年4期
    關(guān)鍵詞:阻尼比間隙

    第一作者 尹俊紅 女,博士生,1986年生

    通信作者 李青寧 男,教授,博導(dǎo),1952年生

    郵箱:lqn952@163.com

    基于精細(xì)積分法的結(jié)構(gòu)碰撞反應(yīng)譜研究

    尹俊紅,李青寧,張瑞杰,程麥理,閆磊,孫建鵬,韓春

    (西安建筑科技大學(xué) 土木工程學(xué)院,西安710055)

    摘要:碰撞是導(dǎo)致結(jié)構(gòu)地震破壞的重要原因之一,其與兩碰撞結(jié)構(gòu)的動力性能和地震作用密切相關(guān),應(yīng)給出與兩個結(jié)構(gòu)動力參數(shù)相關(guān)的反應(yīng)譜,以指導(dǎo)結(jié)構(gòu)考慮碰撞的抗震設(shè)計。本文給出相鄰結(jié)構(gòu)碰撞加速度反應(yīng)譜定義,推導(dǎo)了精細(xì)積分法碰撞計算公式,進(jìn)行了該反應(yīng)譜曲面的分析和探討。結(jié)果表明,精細(xì)積分法適用于結(jié)構(gòu)碰撞問題的計算,可得到高精度的解,無條件穩(wěn)定且可提高計算效率。碰撞會加大結(jié)構(gòu)的加速度響應(yīng);動力特性相同的兩相鄰結(jié)構(gòu)產(chǎn)生的響應(yīng)有差異;增大阻尼比和相鄰結(jié)構(gòu)阻尼比差別、以及設(shè)置足夠大的初始間隙,可有效降低結(jié)構(gòu)的碰撞響應(yīng)。

    關(guān)鍵詞:結(jié)構(gòu)碰撞;精細(xì)積分法;加速度反應(yīng)譜;阻尼比;間隙

    基金項目:國家自然科學(xué)基金資助項目(51078306);國家自然科學(xué)基金資助項目(51408453);高等學(xué)校博士學(xué)科點專項科研基金資助項目(20106120110004);高等學(xué)校博士學(xué)科點專項科研基金資助項目(20136120120022)

    收稿日期:2014-07-24修改稿收到日期:2014-09-30

    中圖分類號:U442.55文獻(xiàn)標(biāo)志碼: A

    Structural impact response spectra based on precise time-integration method

    YINJun-hong,LIQing-ning,ZHANGRui-jie,CHENGMai-li,YANLei,SUNJian-peng,HANChun(College of Civil Engineering, Xi’an University of Architecture and Technology, Xi’an 710055, China)

    Abstract:Collision is one of the main reasons leading to damage of structures under earthquake excitation,it is closely related to structures’ dynamic characteristics and earthquake action. Response spectra for two adjacent impacting structures related to their dynamic parameters should be deduced to guide their aseismic design considering impact. Here, acceleration response spectra for impact between two adjacent structures were defined, the calculation formulas of pounding based on precise time-integration method were derived,and the response spectra were analyzed and discussed. The study results indicated that the precise time-integration method with high accuracy, unconditionally stable and improved efficiency is suitable for the computation of structural impact problems; acceleration responses of structures increase due to pounding, the responses of two adjacent structures with the same dynamic characteristics are different; to increase damping ratios of the two structures and the difference between the two ratios, and to set up enough large initial gap between the two structures can effectively reduce their impact responses.

    Key words: structural pounding; precise time-integration method; acceleration response spectrum; damping ratio; gap

    碰撞會對結(jié)構(gòu)造成極大損害[1-2],碰撞問題被認(rèn)為是影響結(jié)構(gòu)地震反應(yīng)和抗震性能的一個重要因素。

    利用精細(xì)時程積分法無條件穩(wěn)定、精度高和受時間步長限制小的特點,將其引入求解結(jié)構(gòu)的碰撞問題,驗證精細(xì)積分法的計算精度和適用性。并基于精細(xì)積分法進(jìn)行了結(jié)構(gòu)的碰撞加速度反應(yīng)譜研究,反應(yīng)譜表明具有指定阻尼比和頻率的單自由度結(jié)構(gòu)對地面輸入的反應(yīng)程度[8],是根據(jù)地震對單自由度結(jié)構(gòu)的最大影響來表征地面運動特性。對于碰撞問題的反應(yīng)譜研究,Kawashima等[9]給出了相鄰結(jié)構(gòu)最大位移反應(yīng)譜,Ruangrassamee等[10]給出了碰撞影響下的相對位移反應(yīng)譜,王君杰等[11]給出橋梁結(jié)構(gòu)受撞擊的位移沖擊譜及統(tǒng)計規(guī)律,Jankowski R[12]給出結(jié)構(gòu)碰撞力反應(yīng)譜。相對位移譜可用于指導(dǎo)結(jié)構(gòu)防撞設(shè)計中初始間隙的設(shè)置;碰撞力反應(yīng)譜用于評估結(jié)構(gòu)因碰撞造成的潛在損害,并為防撞墻和緩沖裝置的設(shè)計直接提供碰撞力最大值。以上內(nèi)容對結(jié)構(gòu)的碰撞問題研究提供了有益的指導(dǎo),但目前還鮮見針對相鄰結(jié)構(gòu)碰撞加速度反應(yīng)譜的研究與探討,加速度譜可反映碰撞引起的結(jié)構(gòu)加速度響應(yīng)相對于地面輸入地震波的放大效果,并直接得到考慮碰撞的最大地震作用,而后即可根據(jù)反應(yīng)譜理論將其看做靜力作用對結(jié)構(gòu)進(jìn)行強度和剛度等驗算,以指導(dǎo)結(jié)構(gòu)進(jìn)行考慮碰撞的抗震設(shè)計?;诰?xì)積分法,考慮了結(jié)構(gòu)的阻尼比和初始間隙等參數(shù),進(jìn)行了結(jié)構(gòu)的碰撞加速度譜分析,為相鄰結(jié)構(gòu)的碰撞研究和抗震設(shè)計提供參考。

    1精細(xì)積分法在碰撞中的應(yīng)用

    1.1碰撞理論

    結(jié)構(gòu)碰撞理論分析模型見圖1,碰撞結(jié)構(gòu)體系的動力平衡方程為:

    (1)

    圖1 理論模型 Fig.1 Theoretical model

    在碰撞的理論以及模擬分析中,通常采用接觸單元模擬結(jié)構(gòu)的碰撞力。Jan-Hertz-damp模型(非線性-阻尼模型)[13]假定在相互接觸的階段,有可能發(fā)生塑性變形和局部損傷現(xiàn)象,且能量耗散僅在此過程發(fā)生;在分離階段,積累的彈性應(yīng)變能開始釋放,此階段沒有能量耗散。采用Jan-Hertz-damp模型模擬結(jié)構(gòu)的碰撞力,接觸過程中的碰撞力表示為:

    (2)

    (3)

    1.2精細(xì)積分法求解

    對式(1)采用精細(xì)積分法求解:

    (4)

    (5)

    (6)

    式(5)和式(6)代入式(1),有

    (7)

    整理為:

    (8)

    x=[U1(t)U2(t)p1p2]T

    (9)

    (10)

    (11)

    式(8)可化為:

    (12)

    式(12)即為結(jié)構(gòu)碰撞問題的狀態(tài)傳遞方程。一般解為[14]:

    (13)

    將荷載作用離散成步長Δt時間間隔,任意時刻表示為kΔt(k=0,1,2…),tk+1=tk+Δt,則

    (14)

    (15)

    (16)

    令T(Δt)=eHΔt

    (17)

    指數(shù)矩陣T(Δt)求解如下:

    (18)

    將指數(shù)矩陣展開為泰勒級數(shù):

    T(τ)=eHτ=I+Ta

    (19)

    (20)

    式中:l為截斷階數(shù)。

    由于Ta的元素值非常小,若它與直接矩陣I相加后進(jìn)行式(18)的運算,會由于計算機的舍入操作而喪失精度。為避免這一問題,式(18)計算過程如下:

    (Ta)t+1=2(Ta)i+(Ta)i·(Ta)i

    (i=1,2,…,N)

    (21)

    計算得到(Ta)N,則

    T(Δt)=l+(Ta)N

    (22)

    精細(xì)積分法的計算過程通過編程語言實現(xiàn),求解過程為:

    (2) 根據(jù)式(17)得離散時刻tk的狀態(tài)向量,其前兩項即為位移U1(tk)、U2(tk)。

    (4) 若相對位移不超過初始間隙,F(xiàn)c(tk)=0;若相對位移超過初始間隙,采用式(2)計算Fc(tk),并進(jìn)行下一步迭代。

    (6) 重復(fù)以上“(2)”~“(5)”步驟,直至地震動輸入完成。輸出計算結(jié)果。

    2精細(xì)積分法驗證

    為驗證精細(xì)積分法的精度和適用性,通過對兩小球自由彈性碰撞進(jìn)行模擬,圖1中間隙d=0,兩小球質(zhì)量均為 2 kg,自振圓頻率為10.25 rad/s,左側(cè)小球從初始位置向右側(cè)自由運動與右側(cè)小球發(fā)生彈性碰撞,初始時刻位U1(0)=0.1 m,U2(0)=0。假定碰撞過程中無能量損失。

    滿足動量守恒的結(jié)構(gòu)位移響應(yīng)理論解見圖2(a);采用Newmark-β法和Willson-θ法傳統(tǒng)時程分析方法,以及所提出的精細(xì)積分法求解的位移響應(yīng)見圖2(b)圖2(d),求解積分時間步長取0.000 1 s和0.000 2 s。由圖2可知,Newmark-β法和Willson-θ法采用0.000 1 s積分步長的計算結(jié)果與理論解有較大誤差,只有取更小的積分步長才可使其精度增加;且傳統(tǒng)時程積分方法受時間步長影響較大,隨著時間的推移,誤差在積分步長采用0.000 2 s時較0.000 1 s明顯增加,此時計算穩(wěn)定性較差。而精細(xì)積分法采用0.000 1 s積分步長的計算結(jié)果與理論解非常接近,表明與傳統(tǒng)積分法采用相同積分步長時,其計算精度較高;積分步長取0.000 2 s的計算結(jié)果仍然與理論解吻合,說明該方法在碰撞問題計算中也很穩(wěn)定,且總體計算誤差較其它兩種方法小很多,這也意味著選擇相對較大的積分步長不但保證精度和穩(wěn)定性,還可以使計算效率提高。

    采用精細(xì)積分法計算的碰撞力時程見圖3,可以看出,基于精細(xì)積分法求解的結(jié)構(gòu)位移時程與理論位移時程的峰值和相位一致,精細(xì)積分法求解的碰撞次數(shù)為6次,總體與理論相符,但碰撞位置有細(xì)微變化,并且隨著時間推移差異有些許的增加趨勢,這是因為理論解滿足動量守恒定律,碰撞是在瞬間完成的,而采用接觸模型的碰撞力總是有作用過程,這一過程隨著接觸碰撞剛度的增大而減小,其位移解也更接近理論解。

    圖2 質(zhì)點位移響應(yīng) Fig.2 Displacement response of two masses

    綜上,精細(xì)積分法求解碰撞問題是有效的、實用的,由式(20)和式(21)知,精細(xì)積分法理論上可以實現(xiàn)任何精度。但由于碰撞過程短暫,采用精細(xì)積分法理論上宜選擇較短的積分步長,這是由碰撞接觸問題的特點和機理決定,而非精細(xì)積分法本身算法機理決定,根據(jù)本算例驗證,在碰撞問題研究中選擇積分步長0.000 2 s較為合適的,尤其利用精細(xì)積分法的穩(wěn)定性特點,對于非撞擊時刻可選擇更大的積分步長,即變步長算法,可使計算效率大為提高。

    圖3 基于精細(xì)積分法的碰撞力時程 Fig.3 Pounding force time history used precise integration method

    3結(jié)構(gòu)碰撞反應(yīng)譜

    (k=1,2)

    (23)

    基于精細(xì)積分法進(jìn)行結(jié)構(gòu)碰撞加速度反應(yīng)譜曲線研究,圖1中左側(cè)為1#結(jié)構(gòu),右側(cè)為2#結(jié)構(gòu),碰撞參數(shù)選取如下:

    結(jié)構(gòu)質(zhì)量m1=m2=106kg,結(jié)構(gòu)阻尼比ξ1=ξ2=0.05,初始間隙d=0.05 m,相鄰結(jié)構(gòu)的自振周期T1和T2為0.05~3 s,對應(yīng)結(jié)構(gòu)剛度k1和k2為4.39×106~1.58×1010N/m,周期間隔為 0.05 s。輸入地震動采El-centro波(1940年,NS),加速度峰值PGA=3.069 4 m/s2。

    因碰撞結(jié)構(gòu)是相鄰兩個單自由度體系,需要采用三維圖形表示,為了更直觀的觀察對比,將結(jié)構(gòu)無碰撞反應(yīng)譜也表示為三維形式。

    3.1阻尼比影響

    圖4給出相鄰結(jié)構(gòu)阻尼比相等且一致變化的碰撞加速度反應(yīng)譜和無碰撞時的加速度反應(yīng)譜;圖5給出左側(cè)結(jié)構(gòu)阻尼比不變、 右側(cè)結(jié)構(gòu)阻尼比改變的碰撞加

    圖4 x結(jié)構(gòu)阻尼比ξ 1=ξ 2時的加速度反應(yīng)譜 Fig.4 Acceleration spectra for different values of identical damping ratios, ξ 1= ξ 2, of both structures

    速度反應(yīng)譜和無碰撞時的加速度反應(yīng)譜;以下左側(cè)圖為1#結(jié)構(gòu)加速度反應(yīng)譜,右側(cè)圖為2#結(jié)構(gòu)加速度反應(yīng)譜。表1給出不同阻尼比情況下的碰撞加速度譜峰值,可以看出,無碰撞情況下,加速度反應(yīng)譜響應(yīng)隨著周期的增加而減小,碰撞發(fā)生時,這一現(xiàn)象不在存在;碰撞加大了結(jié)構(gòu)的加速度響應(yīng),在結(jié)構(gòu)阻尼比都為0.02時,2#結(jié)構(gòu)的加速度響應(yīng)增加幅值為120%;最大加速度響應(yīng)對應(yīng)的周期點也發(fā)生改變。隨著阻尼比的增大,結(jié)構(gòu)加速度響應(yīng)增加幅度降低。

    圖5 1 #結(jié)構(gòu)阻尼比ξ 1=0.05,2 #結(jié)構(gòu)阻尼比ξ 2變化的加速度反應(yīng)譜 Fig.5 Acceleration spectra for different values of damping ratio of 2 # structure ξ 2and constant value of damping ratio of 1 # structure ξ 1=0.05

    由圖4及表1,阻尼比在0~0.02范圍時,2#結(jié)構(gòu)因碰撞引起的加速度響應(yīng)強于1#結(jié)構(gòu),阻尼比在0.02~0.2范圍,2#結(jié)構(gòu)的加速度響應(yīng)則減弱,低于1#結(jié)構(gòu)響應(yīng),總體差別在3%~26%。表明對于動力特性相同的兩相鄰結(jié)構(gòu),碰撞影響下產(chǎn)生的加速度響應(yīng)有差異。由圖5及表1,阻尼比差別越大,碰撞引起的加速度響應(yīng)隨之減小。

    由圖4~圖5,對于碰撞結(jié)構(gòu),結(jié)構(gòu)阻尼比的增大,使加速度響應(yīng)為零的區(qū)域增加,此區(qū)域相鄰結(jié)構(gòu)的自振周期相當(dāng);增大阻尼比,結(jié)構(gòu)不發(fā)生碰撞的區(qū)域增加,此區(qū)域結(jié)構(gòu)碰撞反應(yīng)譜與無碰撞反應(yīng)譜形狀一致。阻尼比在0~0.05范圍時,碰撞加速度反應(yīng)譜曲線多棱角和毛刺,在阻尼比為0.05~0.2范圍,譜曲線逐漸轉(zhuǎn)為圓滑,表明小阻尼結(jié)構(gòu)受碰撞影響產(chǎn)生的響應(yīng)較大阻尼結(jié)構(gòu)強烈,因阻尼增加使結(jié)構(gòu)的振動衰減加快,從而降低加速度響應(yīng),使譜曲線較為緩和。

    圖6 不同間隙的碰撞加速度反應(yīng)譜 Fig.6 Acceleration spectra for different values of gap size between structures

    表1 碰撞加速度反應(yīng)譜峰值(m/s 2)

    3.2間隙影響

    圖6給出初始間隙不同的相鄰結(jié)構(gòu)加速度反應(yīng)譜和無碰撞的加速度反應(yīng)譜,表2給出碰撞加速度反應(yīng)譜峰值??梢钥闯觯鲎布哟罅私Y(jié)構(gòu)的加速度響應(yīng),在間隙為0.25 m時,碰撞使2#結(jié)構(gòu)加速度響應(yīng)增加47%。

    隨著間隙的逐漸增加,結(jié)構(gòu)的加速度反應(yīng)譜峰值有增有減,在間隙分別為0和0.25 m時,1#和2#結(jié)構(gòu)碰撞加速度響應(yīng)總體降低幅度都為31%,表明設(shè)置足夠大的初始間隙可以減小結(jié)構(gòu)的碰撞響應(yīng)。間隙在00.15 m范圍時,2#結(jié)構(gòu)的加速度響應(yīng)強于1#結(jié)構(gòu),間隙在0.20.25 m范圍時,2#結(jié)構(gòu)響應(yīng)則弱于1#結(jié)構(gòu),總體1#和2#結(jié)構(gòu)因碰撞引起的響應(yīng)差別在07%,表明碰撞引起相鄰結(jié)構(gòu)的加速度響應(yīng)有差異,間隙改變對這種差異的影響并不顯著。

    表2 碰撞加速度反應(yīng)譜峰值(m/s 2)

    由圖6,間隙的增加,使結(jié)構(gòu)碰撞反應(yīng)譜曲線形狀與無碰撞時完全一致的區(qū)域逐漸增加,在間隙增大時,這一現(xiàn)象尤其明顯。在間隙在00.05m范圍,阻尼比相等,相鄰結(jié)構(gòu)周期較小時時碰撞響應(yīng)也相對較劇烈;間隙在0.10.2m范圍,轉(zhuǎn)為長周期時發(fā)生相對較大的加速度響應(yīng)。

    4結(jié)論

    基于精細(xì)積分法對結(jié)構(gòu)進(jìn)行了碰撞加速度反應(yīng)譜分析,精細(xì)積分法具有無條件穩(wěn)定與高精度解的特點,應(yīng)用于結(jié)構(gòu)的碰撞響應(yīng)分析中,提高計算效率和計算精度。主要結(jié)論如下:

    (1) 精細(xì)積分方法理論上可實現(xiàn)任意精度,基于精細(xì)積分法推導(dǎo)的碰撞理論公式適宜于結(jié)構(gòu)碰撞研究與應(yīng)用。由于碰撞發(fā)生時間很短,精細(xì)積分法宜采用變步長以更快速的求解結(jié)構(gòu)碰撞問題。但由積分步長選擇產(chǎn)生的精度誤差并非精細(xì)積分法本身造成的,而是由碰撞發(fā)生時間和機理決定。

    (2) 碰撞會加大結(jié)構(gòu)的加速度響應(yīng),響應(yīng)大小與地震作用和結(jié)構(gòu)動力特性有關(guān);碰撞對動力特性相同的兩相鄰結(jié)構(gòu)造成的響應(yīng)有差異。阻尼比差別引起的相鄰結(jié)構(gòu)的加速度差異較為顯著,間隙引起的響應(yīng)差異則相對小一些。

    (3) 相鄰結(jié)構(gòu)的阻尼比增大和阻尼比差別增加,在地震作用下因碰撞產(chǎn)生的加速度響應(yīng)會降低。

    (4) 若相鄰結(jié)構(gòu)設(shè)置足夠的初始間隙,可減小結(jié)構(gòu)碰撞的可能性和加速度響應(yīng)。

    參 考 文 獻(xiàn)

    [1] Kasai K, Maison B F. Building pounding damage during the 1989 Loma Prieta earthquake[J]. Engineering Structures,1997,19(3):195-207.

    [2] Rosenblueth E, Meli R. The 1985 earthquake: causes and effects in Mexico City[M].1986.

    [3] Newmark N M. A method of computation for structural dynamics[J]. Journal of Engineering Mechanics, 1959, 85(3): 249-260.

    [4] Wilson E L, et al. Nonlinear dynamic analysis of complex structures[J]. Earthquake Engineering and Structural Dynamics, 1973, 1(3): 241-252.

    [5] 鐘萬勰.結(jié)構(gòu)動力方程的精細(xì)時程積分[J].大連理工大學(xué)學(xué)報,1994, 34(4): 131-136.

    ZHONG Wan-xie. On precise time-intergration method for structural dynamics[J]. Journal of Dalian University of Technology, 1994, 34(4): 131-136.

    [6] 鐘萬勰. 暫態(tài)歷程的精細(xì)計算方法[J]. 計算結(jié)構(gòu)力學(xué)及其應(yīng)用, 1995, 12(1): 1-6.

    ZHONG Wan-xie. Precise time-intergration method for the transient process[J]. Computational structural Mechanics and Applications, 1995, 12(1): 1-6.

    [7] 鐘萬勰.計算結(jié)構(gòu)力學(xué)與最優(yōu)控制.[M].大連:大連理工大學(xué)出版社,1993.

    [8] 克拉夫 R,彭津 J.結(jié)構(gòu)動力學(xué)[M].王光遠(yuǎn),譯.北京:高等教育出版社,2006.

    [9] Kawashima K, Sato T. Relative displacement response spectrum and its application[C]. In: Eleventh World Conference on Earthquake Engineering,1996.

    [10] Ruangrassamee A, Kawashima K. Relative displacement response spectra with pounding effect[J]. Earthquake Engineering & Structural Dynamics, 2001, 30(10): 1511-1538.

    [11] 王君杰,喻志然.船撞設(shè)計沖擊譜研究[J].振動與沖擊,2014,33(14):11-14.

    WANG Jun-jie,YU Zhi-ran.Design impact spectrum for vessel-bridge collision analysis[J].Journal of Vibration and Shock,2014,33(14):11-14.

    [12] Jankowski R. Pounding force response spectrum under earthquake excitation[J].Engineering Structures,2006,28(8):1149-1161.

    [13] Jankowski R. Non-linear viscoelastic modelling of earthquake‐induced structural pounding[J]. Earthquake Engineering & Structural Dynamics, 2005,34(6):595-611.

    [14] 常春馨.現(xiàn)代控制理論基礎(chǔ)[M].北京:機械工業(yè)出版社,1988.

    [15] 胡聿賢.地震工程學(xué)[M].北京:地震出版社,2006.

    [16] Mesguita C.Impacto entre Estruturas Adjacentes devido a Accao Sismica, M. Sc. dissertation[M]. Instituto Superisor Tecnico, Lisbon, Portugal,1991.

    [17] Azevedo J, Bento R. Design criteria for buildings subjected to pounding[C]. Eleventh World Conference on Earthquake Engineering, Acapulco, Mexico, 1996.

    [18] van Mier J G M, Pruijssers A F, Reinhardt H W, et al. Load-time response of colliding concrete bodies[J]. Journal of Structural Engineering, 1991, 117(2): 354-374.

    猜你喜歡
    阻尼比間隙
    結(jié)構(gòu)阻尼比對超高層建筑順風(fēng)向風(fēng)振響應(yīng)的影響
    建材世界(2024年1期)2024-03-04 03:57:28
    間隙
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    飛行過載及安裝間隙對主安裝節(jié)推力測量的影響
    緊流形上的Schr?dinger算子的譜間隙估計
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    結(jié)構(gòu)構(gòu)件阻尼比對大跨度懸索橋地震響應(yīng)的影響
    淺談保護(hù)間隙的利弊與應(yīng)用
    廣西電力(2016年4期)2016-07-10 10:23:38
    臭氧發(fā)生器最佳放電間隙的確定及研究
    国产真实伦视频高清在线观看| 又爽又黄无遮挡网站| 亚洲婷婷狠狠爱综合网| 丝袜脚勾引网站| 51国产日韩欧美| 美女xxoo啪啪120秒动态图| 欧美国产精品一级二级三级 | 一区二区三区免费毛片| 国产中年淑女户外野战色| 啦啦啦啦在线视频资源| 欧美一级a爱片免费观看看| 联通29元200g的流量卡| av福利片在线观看| av国产免费在线观看| 18禁裸乳无遮挡动漫免费视频 | 国产综合懂色| 成人一区二区视频在线观看| 色5月婷婷丁香| 成人亚洲精品一区在线观看 | 99久久人妻综合| 亚洲三级黄色毛片| 亚洲经典国产精华液单| 国产精品国产三级专区第一集| 成人国产麻豆网| 五月玫瑰六月丁香| 亚洲aⅴ乱码一区二区在线播放| 久久久色成人| 高清毛片免费看| 内地一区二区视频在线| 亚洲国产精品999| 色网站视频免费| 亚洲av欧美aⅴ国产| 男女啪啪激烈高潮av片| 亚洲av男天堂| 国产成人免费无遮挡视频| 国产精品嫩草影院av在线观看| 欧美xxxx性猛交bbbb| 晚上一个人看的免费电影| 国产日韩欧美亚洲二区| 一个人看视频在线观看www免费| 亚洲国产高清在线一区二区三| 国产黄频视频在线观看| 日本猛色少妇xxxxx猛交久久| 国精品久久久久久国模美| 亚洲精品日本国产第一区| 蜜桃亚洲精品一区二区三区| 中文字幕久久专区| av天堂中文字幕网| 哪个播放器可以免费观看大片| av又黄又爽大尺度在线免费看| 中国国产av一级| videossex国产| 97热精品久久久久久| 毛片女人毛片| 高清欧美精品videossex| 国产成人a∨麻豆精品| 亚洲av免费高清在线观看| 日韩精品有码人妻一区| 日韩成人伦理影院| 国产精品麻豆人妻色哟哟久久| 亚洲不卡免费看| 99久久人妻综合| 午夜精品国产一区二区电影 | 97在线人人人人妻| 亚洲国产成人一精品久久久| 欧美高清性xxxxhd video| 麻豆精品久久久久久蜜桃| 男人和女人高潮做爰伦理| 一区二区三区乱码不卡18| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av天美| kizo精华| 永久网站在线| 亚洲美女视频黄频| 偷拍熟女少妇极品色| 禁无遮挡网站| 伊人久久国产一区二区| 免费看光身美女| 国产成人精品久久久久久| av国产免费在线观看| av一本久久久久| 亚洲va在线va天堂va国产| 国产视频内射| 精品久久久久久久人妻蜜臀av| 久久久久久久久久人人人人人人| 精品久久久久久久久亚洲| 欧美另类一区| 国产人妻一区二区三区在| videos熟女内射| 99re6热这里在线精品视频| 国内少妇人妻偷人精品xxx网站| 亚洲国产av新网站| 日本av手机在线免费观看| 国产精品嫩草影院av在线观看| 伊人久久国产一区二区| 一个人看的www免费观看视频| 日本熟妇午夜| 欧美人与善性xxx| 亚洲图色成人| 日韩一区二区三区影片| 成年人午夜在线观看视频| 日本一二三区视频观看| 国产精品一区www在线观看| 国产淫片久久久久久久久| 国产乱来视频区| 狂野欧美白嫩少妇大欣赏| 日韩欧美精品免费久久| 亚洲高清免费不卡视频| 18禁裸乳无遮挡免费网站照片| 嫩草影院入口| 国产av国产精品国产| 亚洲色图综合在线观看| 国产av国产精品国产| 久久精品久久精品一区二区三区| 国产高清有码在线观看视频| 视频中文字幕在线观看| 国产高清有码在线观看视频| 国产免费视频播放在线视频| 亚洲精品久久午夜乱码| 六月丁香七月| 日韩免费高清中文字幕av| 国产精品熟女久久久久浪| 水蜜桃什么品种好| 精品人妻熟女av久视频| 国产淫语在线视频| 成年女人看的毛片在线观看| 18禁裸乳无遮挡动漫免费视频 | 亚洲精品一区蜜桃| av福利片在线观看| 日日撸夜夜添| 婷婷色综合www| 国产91av在线免费观看| 少妇高潮的动态图| 日韩欧美精品v在线| 色哟哟·www| 高清欧美精品videossex| av线在线观看网站| av福利片在线观看| 久久精品久久久久久噜噜老黄| 搡老乐熟女国产| 欧美成人一区二区免费高清观看| 一本一本综合久久| 欧美丝袜亚洲另类| 亚洲不卡免费看| 亚洲国产精品国产精品| 熟妇人妻不卡中文字幕| 夜夜爽夜夜爽视频| 大又大粗又爽又黄少妇毛片口| 夜夜爽夜夜爽视频| 国产一区有黄有色的免费视频| 中文字幕久久专区| 亚洲最大成人中文| 国产熟女欧美一区二区| 蜜臀久久99精品久久宅男| 只有这里有精品99| 国产精品秋霞免费鲁丝片| 亚洲精品色激情综合| 国产黄a三级三级三级人| 久久ye,这里只有精品| 草草在线视频免费看| 亚洲国产欧美在线一区| 中国国产av一级| 欧美成人一区二区免费高清观看| 欧美xxxx黑人xx丫x性爽| 91精品伊人久久大香线蕉| 精品久久久精品久久久| 欧美xxxx黑人xx丫x性爽| av在线天堂中文字幕| 欧美一区二区亚洲| 欧美xxxx黑人xx丫x性爽| 中文字幕制服av| 久久精品国产亚洲av涩爱| 国产精品无大码| 日日撸夜夜添| 18禁裸乳无遮挡免费网站照片| 午夜精品国产一区二区电影 | 亚洲天堂国产精品一区在线| 久久精品国产亚洲av涩爱| 成年女人看的毛片在线观看| 成人国产av品久久久| 国产欧美日韩一区二区三区在线 | 天美传媒精品一区二区| 国产黄片视频在线免费观看| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久精品电影小说 | 黄色怎么调成土黄色| 国产成人免费无遮挡视频| 国内精品宾馆在线| 午夜激情福利司机影院| 性色av一级| 中文在线观看免费www的网站| 小蜜桃在线观看免费完整版高清| 一级爰片在线观看| 天堂网av新在线| 欧美+日韩+精品| 五月开心婷婷网| 亚洲第一区二区三区不卡| 白带黄色成豆腐渣| av在线蜜桃| 久久99蜜桃精品久久| 深夜a级毛片| 精品国产乱码久久久久久小说| 26uuu在线亚洲综合色| 日本午夜av视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美人成| 深爱激情五月婷婷| 又粗又硬又长又爽又黄的视频| 精品久久久久久电影网| 亚洲精品乱码久久久久久按摩| 十八禁网站网址无遮挡 | 亚洲丝袜综合中文字幕| 丝瓜视频免费看黄片| 18禁动态无遮挡网站| 直男gayav资源| 直男gayav资源| 午夜激情福利司机影院| 久久久精品欧美日韩精品| 国语对白做爰xxxⅹ性视频网站| 亚洲av免费高清在线观看| 国产精品秋霞免费鲁丝片| 国产免费福利视频在线观看| 精品视频人人做人人爽| 国产精品一区www在线观看| 国产精品久久久久久精品电影| 亚洲激情五月婷婷啪啪| kizo精华| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 好男人在线观看高清免费视频| 美女脱内裤让男人舔精品视频| 一级二级三级毛片免费看| 国产v大片淫在线免费观看| 深夜a级毛片| 日本欧美国产在线视频| 老司机影院毛片| 日韩不卡一区二区三区视频在线| 国产视频内射| 免费不卡的大黄色大毛片视频在线观看| 高清毛片免费看| 午夜福利视频1000在线观看| 国产精品一区二区三区四区免费观看| 99久久人妻综合| 午夜爱爱视频在线播放| 日韩欧美精品免费久久| 一个人观看的视频www高清免费观看| 日韩成人av中文字幕在线观看| 亚洲精品日本国产第一区| av免费在线看不卡| av免费在线看不卡| 久久鲁丝午夜福利片| 欧美最新免费一区二区三区| 观看免费一级毛片| 日本三级黄在线观看| 国产欧美日韩一区二区三区在线 | 久久精品国产亚洲网站| 性插视频无遮挡在线免费观看| 国产黄片视频在线免费观看| 中国国产av一级| 亚洲av一区综合| 观看美女的网站| 亚洲欧洲日产国产| 搡老乐熟女国产| 久久ye,这里只有精品| 久久久精品94久久精品| 中国国产av一级| 国产淫片久久久久久久久| 精品久久久噜噜| 欧美精品人与动牲交sv欧美| 亚洲最大成人av| 欧美3d第一页| 免费大片黄手机在线观看| 精品人妻视频免费看| 欧美性猛交╳xxx乱大交人| 国产高清有码在线观看视频| 欧美日韩亚洲高清精品| 精品少妇黑人巨大在线播放| 亚洲天堂av无毛| 美女内射精品一级片tv| 男女那种视频在线观看| 久久热精品热| 亚洲精品国产色婷婷电影| 啦啦啦啦在线视频资源| 国产成人精品福利久久| 久久久成人免费电影| 亚洲欧美清纯卡通| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 欧美极品一区二区三区四区| 美女脱内裤让男人舔精品视频| 欧美精品国产亚洲| 久久精品夜色国产| 伊人久久精品亚洲午夜| 一区二区av电影网| 日韩免费高清中文字幕av| 日韩电影二区| 久久精品国产自在天天线| 99久久精品国产国产毛片| 精品亚洲乱码少妇综合久久| 熟女av电影| 日产精品乱码卡一卡2卡三| 久久久久久国产a免费观看| 日韩亚洲欧美综合| 伊人久久国产一区二区| 女人久久www免费人成看片| 国产 一区精品| 色视频在线一区二区三区| 成人二区视频| 国产精品99久久久久久久久| 久久久亚洲精品成人影院| 大码成人一级视频| 久久国内精品自在自线图片| 中文字幕制服av| 国产伦在线观看视频一区| 亚洲一区二区三区欧美精品 | 内地一区二区视频在线| 一个人看视频在线观看www免费| av网站免费在线观看视频| 欧美精品国产亚洲| 亚洲最大成人手机在线| 99久久人妻综合| 精品酒店卫生间| 少妇人妻 视频| 一二三四中文在线观看免费高清| 如何舔出高潮| 日日摸夜夜添夜夜添av毛片| 日韩人妻高清精品专区| 精品久久久久久久久av| 男人爽女人下面视频在线观看| 日韩电影二区| 大话2 男鬼变身卡| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久久久午夜电影| 国产精品一区二区三区四区免费观看| 亚洲欧美精品自产自拍| av免费在线看不卡| 一区二区三区免费毛片| 老师上课跳d突然被开到最大视频| 欧美性猛交╳xxx乱大交人| 男人爽女人下面视频在线观看| 精品视频人人做人人爽| 精品一区二区三卡| av线在线观看网站| 国产av码专区亚洲av| 亚洲精品成人久久久久久| 日本av手机在线免费观看| 性色av一级| 特级一级黄色大片| 国产亚洲5aaaaa淫片| a级毛片免费高清观看在线播放| 色婷婷久久久亚洲欧美| 可以在线观看毛片的网站| 男的添女的下面高潮视频| 国产成年人精品一区二区| 午夜福利视频1000在线观看| 老女人水多毛片| 99久久九九国产精品国产免费| 97人妻精品一区二区三区麻豆| 超碰av人人做人人爽久久| 国产精品国产三级国产专区5o| 人妻一区二区av| 亚洲真实伦在线观看| 国产精品爽爽va在线观看网站| 在线免费十八禁| 99久久精品一区二区三区| 国产爱豆传媒在线观看| 亚洲国产欧美人成| 51国产日韩欧美| 看黄色毛片网站| 国产免费一级a男人的天堂| 精品午夜福利在线看| 国产精品熟女久久久久浪| 黄色欧美视频在线观看| 久久热精品热| 亚州av有码| 国产黄色视频一区二区在线观看| 国产精品国产av在线观看| av网站免费在线观看视频| 舔av片在线| 成年av动漫网址| 久久久久久久久久人人人人人人| 高清日韩中文字幕在线| 精品久久久精品久久久| 日韩一区二区视频免费看| av在线观看视频网站免费| 亚洲国产精品成人综合色| 国产色婷婷99| 中文字幕av成人在线电影| 美女xxoo啪啪120秒动态图| 色视频www国产| 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久 | 熟女av电影| 久久精品国产自在天天线| 插阴视频在线观看视频| 成年女人在线观看亚洲视频 | 国产美女午夜福利| 亚洲国产精品999| 搡女人真爽免费视频火全软件| 99久久精品国产国产毛片| 亚洲人成网站在线观看播放| 深夜a级毛片| 插阴视频在线观看视频| 午夜免费男女啪啪视频观看| 日韩中字成人| 久久久久九九精品影院| 久热这里只有精品99| 久久久午夜欧美精品| 免费观看无遮挡的男女| 成人美女网站在线观看视频| 国产精品国产三级专区第一集| 色播亚洲综合网| 亚洲av成人精品一区久久| 免费高清在线观看视频在线观看| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 亚洲av中文av极速乱| 91在线精品国自产拍蜜月| 熟妇人妻不卡中文字幕| 大话2 男鬼变身卡| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 午夜福利在线观看免费完整高清在| 综合色av麻豆| 中文字幕制服av| 人体艺术视频欧美日本| 国产精品爽爽va在线观看网站| 女人被狂操c到高潮| 伦理电影大哥的女人| 亚洲国产欧美在线一区| 麻豆乱淫一区二区| 国产一区二区三区综合在线观看 | 免费黄网站久久成人精品| 在线天堂最新版资源| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的 | 国产永久视频网站| 精品一区在线观看国产| 亚洲美女搞黄在线观看| 最近中文字幕2019免费版| 欧美国产精品一级二级三级 | 日韩欧美精品v在线| 18禁在线无遮挡免费观看视频| 精品久久久噜噜| 欧美bdsm另类| 久久久久精品性色| 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲成人一二三区av| 精品熟女少妇av免费看| 精品久久久久久久久av| 久久久久国产网址| 夜夜看夜夜爽夜夜摸| 中国国产av一级| 尤物成人国产欧美一区二区三区| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区成人| 成人欧美大片| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 大陆偷拍与自拍| 黄色怎么调成土黄色| 精品久久久噜噜| 五月伊人婷婷丁香| 水蜜桃什么品种好| 在线a可以看的网站| 国产精品女同一区二区软件| 乱码一卡2卡4卡精品| 国产精品成人在线| 在现免费观看毛片| 九草在线视频观看| 综合色av麻豆| 国产日韩欧美亚洲二区| 亚洲欧美日韩无卡精品| 国产精品嫩草影院av在线观看| 亚洲性久久影院| 国产成人精品婷婷| 国产爽快片一区二区三区| 久久精品人妻少妇| 国语对白做爰xxxⅹ性视频网站| 18+在线观看网站| 国产成年人精品一区二区| 色视频www国产| 中文天堂在线官网| 亚洲婷婷狠狠爱综合网| 内射极品少妇av片p| 国产精品久久久久久精品电影| 久久精品久久久久久久性| 亚洲欧美精品自产自拍| 国精品久久久久久国模美| 久久久精品免费免费高清| 久久精品国产自在天天线| 99热6这里只有精品| 成年女人看的毛片在线观看| 国产免费视频播放在线视频| 国产毛片在线视频| xxx大片免费视频| 高清在线视频一区二区三区| 成人黄色视频免费在线看| 在线免费十八禁| 岛国毛片在线播放| 成人漫画全彩无遮挡| 国产探花在线观看一区二区| 哪个播放器可以免费观看大片| 麻豆成人午夜福利视频| 成人免费观看视频高清| 国产精品精品国产色婷婷| 国产探花在线观看一区二区| 成人国产麻豆网| 亚洲成人一二三区av| 日韩av不卡免费在线播放| 能在线免费看毛片的网站| 国产探花在线观看一区二区| 日日撸夜夜添| 久久精品国产a三级三级三级| 亚洲成人中文字幕在线播放| 国产69精品久久久久777片| 神马国产精品三级电影在线观看| 免费观看性生交大片5| 亚洲最大成人中文| 综合色丁香网| 中文资源天堂在线| 精品久久国产蜜桃| 精品久久久噜噜| 一本色道久久久久久精品综合| 尤物成人国产欧美一区二区三区| 成人国产麻豆网| 亚洲精品国产av蜜桃| 日本午夜av视频| 热re99久久精品国产66热6| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说 | 国产乱来视频区| 久久久久久九九精品二区国产| 99久久精品国产国产毛片| 成人亚洲欧美一区二区av| 欧美成人a在线观看| 欧美激情久久久久久爽电影| 久久久精品欧美日韩精品| 精品久久久久久久人妻蜜臀av| 麻豆国产97在线/欧美| 亚洲丝袜综合中文字幕| 2018国产大陆天天弄谢| 午夜精品一区二区三区免费看| 国产精品国产三级国产av玫瑰| 欧美 日韩 精品 国产| 日本色播在线视频| 又爽又黄无遮挡网站| 91aial.com中文字幕在线观看| 久久精品国产亚洲网站| 精品人妻偷拍中文字幕| 在现免费观看毛片| 永久网站在线| 国产亚洲av嫩草精品影院| 国产熟女欧美一区二区| 国产免费一区二区三区四区乱码| 国产毛片a区久久久久| 亚洲欧美日韩卡通动漫| 亚洲经典国产精华液单| 日日摸夜夜添夜夜爱| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 国产成人精品福利久久| 欧美zozozo另类| 色综合色国产| 免费观看无遮挡的男女| 五月开心婷婷网| 少妇猛男粗大的猛烈进出视频 | 一级二级三级毛片免费看| 国产精品不卡视频一区二区| 男女国产视频网站| 精品视频人人做人人爽| 欧美精品国产亚洲| 免费少妇av软件| 亚洲内射少妇av| 亚洲在线观看片| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久精品久久久久真实原创| 香蕉精品网在线| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 哪个播放器可以免费观看大片| 一级黄片播放器| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说 | 内地一区二区视频在线| 国产精品av视频在线免费观看| 久久精品国产亚洲av天美| 国产精品秋霞免费鲁丝片| 极品少妇高潮喷水抽搐| 久久久久久伊人网av| 国产色爽女视频免费观看| 在线观看三级黄色| 欧美一级a爱片免费观看看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产男人的电影天堂91| 日韩电影二区| 久久久久久久午夜电影| 秋霞伦理黄片| 一级二级三级毛片免费看| 亚洲国产精品999| 最近中文字幕高清免费大全6| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 一本久久精品| 日本黄大片高清| 乱系列少妇在线播放| 久久99精品国语久久久| 如何舔出高潮| 久久99热这里只有精品18| 亚洲不卡免费看| 免费电影在线观看免费观看| 精品人妻一区二区三区麻豆|