• <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ā)生器最佳放電間隙的確定及研究
    法律面前人人平等表现在哪些方面| 特级一级黄色大片| 国产成人aa在线观看| 精品不卡国产一区二区三区| 日日干狠狠操夜夜爽| 特级一级黄色大片| 久久精品国产亚洲av高清一级| 女人爽到高潮嗷嗷叫在线视频| 国产精品爽爽va在线观看网站| 日日夜夜操网爽| 黄色视频不卡| 日本一本二区三区精品| 两个人免费观看高清视频| 色综合欧美亚洲国产小说| 在线观看一区二区三区| 午夜福利成人在线免费观看| 国产一区在线观看成人免费| 国产单亲对白刺激| 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| 亚洲精品一区av在线观看| 成年免费大片在线观看| 日本五十路高清| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 国产区一区二久久| 国内少妇人妻偷人精品xxx网站 | 人人妻,人人澡人人爽秒播| 老司机午夜十八禁免费视频| 国产三级在线视频| 欧美3d第一页| 男女下面进入的视频免费午夜| 午夜福利成人在线免费观看| 丰满的人妻完整版| 黄色a级毛片大全视频| 国产av不卡久久| 12—13女人毛片做爰片一| 又粗又爽又猛毛片免费看| 美女扒开内裤让男人捅视频| 男插女下体视频免费在线播放| 欧美性猛交黑人性爽| 久久这里只有精品中国| 一区二区三区国产精品乱码| 亚洲欧美一区二区三区黑人| 小说图片视频综合网站| 亚洲成人久久爱视频| 国产真人三级小视频在线观看| 国产aⅴ精品一区二区三区波| 19禁男女啪啪无遮挡网站| 香蕉久久夜色| 精品久久久久久久末码| 99久久精品热视频| 老汉色av国产亚洲站长工具| 两个人免费观看高清视频| 日韩免费av在线播放| 成年版毛片免费区| 国产99久久九九免费精品| 精品欧美一区二区三区在线| 国产麻豆成人av免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 麻豆一二三区av精品| 欧美日韩乱码在线| 熟女电影av网| 97超级碰碰碰精品色视频在线观看| avwww免费| 国产精品电影一区二区三区| 国产精品av久久久久免费| 黄色成人免费大全| 狂野欧美白嫩少妇大欣赏| 欧美在线黄色| 一级毛片女人18水好多| 日本免费a在线| 亚洲国产中文字幕在线视频| 99热只有精品国产| 亚洲国产欧美网| 午夜福利在线观看吧| 亚洲专区中文字幕在线| 欧美极品一区二区三区四区| 十八禁网站免费在线| 两性午夜刺激爽爽歪歪视频在线观看 | 观看免费一级毛片| 极品教师在线免费播放| 久热爱精品视频在线9| 黄色片一级片一级黄色片| avwww免费| 9191精品国产免费久久| 老鸭窝网址在线观看| 国产99白浆流出| 精品欧美一区二区三区在线| 老汉色av国产亚洲站长工具| 久久久久国产精品人妻aⅴ院| 1024手机看黄色片| 中文字幕人成人乱码亚洲影| 国产激情欧美一区二区| 男女视频在线观看网站免费 | 免费搜索国产男女视频| 午夜福利高清视频| 日韩欧美在线二视频| 久久久久久大精品| 亚洲人成网站在线播放欧美日韩| 亚洲免费av在线视频| 十八禁网站免费在线| 欧美日韩精品网址| 热99re8久久精品国产| 久久精品国产清高在天天线| 韩国av一区二区三区四区| 午夜福利视频1000在线观看| 欧美在线黄色| 精品高清国产在线一区| 桃红色精品国产亚洲av| www日本黄色视频网| 99国产精品99久久久久| 男插女下体视频免费在线播放| 97人妻精品一区二区三区麻豆| 免费观看精品视频网站| 男女视频在线观看网站免费 | 日本一二三区视频观看| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧美网| 久久香蕉激情| 国产精品一及| 欧美日韩精品网址| 国产成人欧美在线观看| videosex国产| 国产99白浆流出| 国产片内射在线| 久久精品国产亚洲av高清一级| 男人舔奶头视频| 免费一级毛片在线播放高清视频| 午夜免费激情av| 久久精品夜夜夜夜夜久久蜜豆 | 又紧又爽又黄一区二区| 色av中文字幕| 日韩欧美免费精品| 男人舔奶头视频| 成年版毛片免费区| 亚洲乱码一区二区免费版| 亚洲乱码一区二区免费版| 好男人电影高清在线观看| 亚洲一区二区三区色噜噜| 黑人巨大精品欧美一区二区mp4| 老司机深夜福利视频在线观看| 女人被狂操c到高潮| 久久中文字幕人妻熟女| 日韩精品青青久久久久久| 真人做人爱边吃奶动态| 18禁美女被吸乳视频| 波多野结衣高清无吗| 久久久久久久久中文| 99久久国产精品久久久| 日韩三级视频一区二区三区| 悠悠久久av| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区久久| 亚洲激情在线av| 亚洲美女视频黄频| av在线播放免费不卡| 18禁黄网站禁片免费观看直播| 在线视频色国产色| 国产av不卡久久| 亚洲av成人不卡在线观看播放网| 色综合婷婷激情| 亚洲欧美激情综合另类| 国产亚洲精品久久久久5区| 无遮挡黄片免费观看| 成人高潮视频无遮挡免费网站| 国产精品免费视频内射| 亚洲片人在线观看| а√天堂www在线а√下载| 美女 人体艺术 gogo| 天堂av国产一区二区熟女人妻 | 又爽又黄无遮挡网站| 黄频高清免费视频| 男人舔女人下体高潮全视频| 亚洲欧美精品综合一区二区三区| 日本在线视频免费播放| 国模一区二区三区四区视频 | 婷婷精品国产亚洲av在线| 亚洲av成人一区二区三| e午夜精品久久久久久久| 精品国产乱码久久久久久男人| 欧美日本亚洲视频在线播放| 搞女人的毛片| www.精华液| 国产精华一区二区三区| 19禁男女啪啪无遮挡网站| 99国产极品粉嫩在线观看| 特级一级黄色大片| 欧美国产日韩亚洲一区| 男人的好看免费观看在线视频 | 成熟少妇高潮喷水视频| 久久亚洲真实| 亚洲成人国产一区在线观看| 老熟妇仑乱视频hdxx| 精品高清国产在线一区| 精品久久蜜臀av无| 欧美成人一区二区免费高清观看 | 男人舔奶头视频| 99国产精品一区二区三区| 精品午夜福利视频在线观看一区| 老鸭窝网址在线观看| 波多野结衣高清作品| 国产精品影院久久| 天天一区二区日本电影三级| 国产av一区在线观看免费| 欧美人与性动交α欧美精品济南到| 国产亚洲欧美在线一区二区| 久久欧美精品欧美久久欧美| 两个人看的免费小视频| 国产精品 国内视频| 两个人的视频大全免费| 国产亚洲精品av在线| 熟女少妇亚洲综合色aaa.| 国产精华一区二区三区| a级毛片a级免费在线| 午夜福利18| 久久香蕉精品热| 看黄色毛片网站| 一卡2卡三卡四卡精品乱码亚洲| 给我免费播放毛片高清在线观看| 一个人免费在线观看的高清视频| 一a级毛片在线观看| 国产99白浆流出| 国产精品亚洲美女久久久| 级片在线观看| 国内精品久久久久久久电影| 成人午夜高清在线视频| 美女高潮喷水抽搐中文字幕| 午夜福利在线在线| 老熟妇乱子伦视频在线观看| 免费观看精品视频网站| 巨乳人妻的诱惑在线观看| 国产亚洲精品久久久久久毛片| 黄色毛片三级朝国网站| 成熟少妇高潮喷水视频| 国产真实乱freesex| 亚洲精品中文字幕一二三四区| 丰满人妻一区二区三区视频av | 精品人妻1区二区| 亚洲精品一卡2卡三卡4卡5卡| 天堂av国产一区二区熟女人妻 | 免费在线观看日本一区| 老司机午夜福利在线观看视频| 成年人黄色毛片网站| 亚洲精品中文字幕一二三四区| 亚洲熟妇中文字幕五十中出| 精华霜和精华液先用哪个| 两个人的视频大全免费| 手机成人av网站| 男男h啪啪无遮挡| 91大片在线观看| 人妻丰满熟妇av一区二区三区| ponron亚洲| 国产三级在线视频| 免费在线观看完整版高清| 黄色视频,在线免费观看| 日本 欧美在线| 女警被强在线播放| 又爽又黄无遮挡网站| 亚洲人成伊人成综合网2020| 嫁个100分男人电影在线观看| 性色av乱码一区二区三区2| 亚洲国产欧美人成| 国产精品免费一区二区三区在线| 少妇人妻一区二区三区视频| 18禁观看日本| 全区人妻精品视频| 国产精品乱码一区二三区的特点| 一级毛片高清免费大全| 亚洲国产精品久久男人天堂| 欧美丝袜亚洲另类 | 高潮久久久久久久久久久不卡| 国产野战对白在线观看| 精品免费久久久久久久清纯| 精品久久久久久久久久免费视频| 动漫黄色视频在线观看| 免费看美女性在线毛片视频| 午夜久久久久精精品| 97碰自拍视频| 亚洲av第一区精品v没综合| 亚洲人成网站高清观看| 国产午夜福利久久久久久| 亚洲性夜色夜夜综合| svipshipincom国产片| 妹子高潮喷水视频| 欧美性猛交╳xxx乱大交人| 色尼玛亚洲综合影院| 男人舔女人的私密视频| 成人午夜高清在线视频| 国产成人一区二区三区免费视频网站| 日韩大尺度精品在线看网址| 成人国产一区最新在线观看| 黄频高清免费视频| 最新美女视频免费是黄的| 99国产极品粉嫩在线观看| 国产成人精品无人区| 久久热在线av| 亚洲 欧美 日韩 在线 免费| 宅男免费午夜| 欧美色欧美亚洲另类二区| 少妇熟女aⅴ在线视频| 五月玫瑰六月丁香| 老司机福利观看| 成人av一区二区三区在线看| 久久亚洲真实| 亚洲国产精品久久男人天堂| 国内少妇人妻偷人精品xxx网站 | cao死你这个sao货| 亚洲国产中文字幕在线视频| 此物有八面人人有两片| 久久久久久久精品吃奶| 亚洲av中文字字幕乱码综合| 少妇人妻一区二区三区视频| 午夜精品久久久久久毛片777| 99久久99久久久精品蜜桃| 少妇被粗大的猛进出69影院| 精品一区二区三区四区五区乱码| 欧美激情久久久久久爽电影| 一二三四社区在线视频社区8| 国产97色在线日韩免费| 亚洲 欧美 日韩 在线 免费| 精品免费久久久久久久清纯| 日日夜夜操网爽| 嫩草影视91久久| 中国美女看黄片| 超碰成人久久| 琪琪午夜伦伦电影理论片6080| 蜜桃久久精品国产亚洲av| 成人国产一区最新在线观看| 老鸭窝网址在线观看| 亚洲av片天天在线观看| 精品国产美女av久久久久小说| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 欧美成人性av电影在线观看| 午夜精品一区二区三区免费看| 啦啦啦免费观看视频1| 亚洲国产欧美一区二区综合| 一区二区三区国产精品乱码| 在线免费观看的www视频| 老司机午夜福利在线观看视频| 不卡一级毛片| 老鸭窝网址在线观看| 看片在线看免费视频| 色综合亚洲欧美另类图片| 国产成人一区二区三区免费视频网站| 日本熟妇午夜| 99在线人妻在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 天天躁夜夜躁狠狠躁躁| 俺也久久电影网| 国产精品乱码一区二三区的特点| 成人午夜高清在线视频| 亚洲男人的天堂狠狠| 亚洲专区中文字幕在线| 久久精品91无色码中文字幕| 99久久99久久久精品蜜桃| 国产精品,欧美在线| 久久草成人影院| 成人国产一区最新在线观看| 亚洲国产欧美一区二区综合| 亚洲 欧美一区二区三区| 亚洲美女黄片视频| 嫁个100分男人电影在线观看| 成人高潮视频无遮挡免费网站| 久久久国产精品麻豆| 一本精品99久久精品77| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 成年版毛片免费区| 91av网站免费观看| 免费一级毛片在线播放高清视频| 色综合亚洲欧美另类图片| 可以免费在线观看a视频的电影网站| 熟女电影av网| 国产黄a三级三级三级人| 国产精品永久免费网站| 美女大奶头视频| 国产成人aa在线观看| 午夜两性在线视频| 亚洲人成77777在线视频| 精品久久久久久久久久久久久| 男人舔奶头视频| 精品久久久久久久末码| 好看av亚洲va欧美ⅴa在| 日韩三级视频一区二区三区| 天堂影院成人在线观看| 大型黄色视频在线免费观看| 1024视频免费在线观看| 亚洲国产欧美网| 亚洲乱码一区二区免费版| 男插女下体视频免费在线播放| 老司机深夜福利视频在线观看| 制服丝袜大香蕉在线| 老司机深夜福利视频在线观看| а√天堂www在线а√下载| 亚洲一区中文字幕在线| 亚洲精品久久国产高清桃花| 三级国产精品欧美在线观看 | 亚洲精品色激情综合| 国产精品久久电影中文字幕| 国产精品久久久久久亚洲av鲁大| 精品国产乱码久久久久久男人| or卡值多少钱| 日韩大尺度精品在线看网址| 亚洲,欧美精品.| 欧美国产日韩亚洲一区| 一区福利在线观看| 两人在一起打扑克的视频| 日本免费一区二区三区高清不卡| 亚洲五月婷婷丁香| 国产乱人伦免费视频| 亚洲一码二码三码区别大吗| 亚洲专区字幕在线| 国产精品久久久久久亚洲av鲁大| 嫁个100分男人电影在线观看| 欧美日韩瑟瑟在线播放| 亚洲av日韩精品久久久久久密| 精品久久久久久久久久久久久| 在线视频色国产色| aaaaa片日本免费| 久久久久久九九精品二区国产 | 女人爽到高潮嗷嗷叫在线视频| 久久久久久久午夜电影| 丰满人妻熟妇乱又伦精品不卡| 亚洲美女视频黄频| 亚洲乱码一区二区免费版| 真人一进一出gif抽搐免费| 久久久久国产精品人妻aⅴ院| 99精品欧美一区二区三区四区| 国产97色在线日韩免费| 黄色成人免费大全| 国产午夜精品久久久久久| 精品久久久久久久久久久久久| 欧美一区二区精品小视频在线| 少妇裸体淫交视频免费看高清 | 精品国产超薄肉色丝袜足j| 18美女黄网站色大片免费观看| 白带黄色成豆腐渣| 亚洲精品中文字幕一二三四区| 国产片内射在线| 国产精品av视频在线免费观看| 欧美极品一区二区三区四区| 国产精品影院久久| 亚洲人成伊人成综合网2020| 国产一区二区三区视频了| 中文字幕熟女人妻在线| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三| 高清在线国产一区| 久久99热这里只有精品18| 国产v大片淫在线免费观看| 成人国产一区最新在线观看| 十八禁人妻一区二区| 成人精品一区二区免费| 欧美在线黄色| 好男人在线观看高清免费视频| 久久精品亚洲精品国产色婷小说| 亚洲18禁久久av| 国产单亲对白刺激| 国产三级在线视频| 欧美色欧美亚洲另类二区| 女人高潮潮喷娇喘18禁视频| 成年免费大片在线观看| 男女午夜视频在线观看| 免费在线观看黄色视频的| 我的老师免费观看完整版| 国产真实乱freesex| 免费在线观看影片大全网站| 欧美色视频一区免费| 免费在线观看影片大全网站| 久久午夜亚洲精品久久| 久久精品91蜜桃| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 国产精品 国内视频| 一级黄色大片毛片| 国产精品久久视频播放| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放 | 视频区欧美日本亚洲| 日日夜夜操网爽| av有码第一页| 一本精品99久久精品77| 国产高清激情床上av| 亚洲国产中文字幕在线视频| 国产又色又爽无遮挡免费看| 精品无人区乱码1区二区| 999久久久国产精品视频| 久久香蕉激情| 免费在线观看影片大全网站| 给我免费播放毛片高清在线观看| 色老头精品视频在线观看| 国产黄色小视频在线观看| 久久 成人 亚洲| av有码第一页| 99精品在免费线老司机午夜| 亚洲国产欧洲综合997久久,| 一本综合久久免费| 午夜免费激情av| 午夜福利在线在线| 国产片内射在线| 久久久久久久午夜电影| 老司机午夜福利在线观看视频| 欧美午夜高清在线| 亚洲午夜理论影院| 中国美女看黄片| 男人舔奶头视频| 国产久久久一区二区三区| 亚洲精品粉嫩美女一区| 亚洲午夜精品一区,二区,三区| 1024手机看黄色片| 精华霜和精华液先用哪个| 99热这里只有是精品50| 俄罗斯特黄特色一大片| 午夜两性在线视频| 成年免费大片在线观看| 亚洲,欧美精品.| 欧美乱码精品一区二区三区| 又大又爽又粗| 国产亚洲精品第一综合不卡| xxx96com| 国产亚洲精品久久久久5区| 亚洲电影在线观看av| 国产97色在线日韩免费| 岛国视频午夜一区免费看| 好看av亚洲va欧美ⅴa在| 俺也久久电影网| 免费看a级黄色片| 国产又黄又爽又无遮挡在线| 国产精品精品国产色婷婷| 久久香蕉精品热| 国产区一区二久久| 日韩欧美三级三区| 搡老熟女国产l中国老女人| 麻豆av在线久日| 人人妻,人人澡人人爽秒播| 正在播放国产对白刺激| 极品教师在线免费播放| 在线观看免费视频日本深夜| 国产97色在线日韩免费| 男女床上黄色一级片免费看| 在线十欧美十亚洲十日本专区| 麻豆国产av国片精品| 亚洲国产看品久久| 窝窝影院91人妻| 一个人观看的视频www高清免费观看 | 国产在线精品亚洲第一网站| 少妇裸体淫交视频免费看高清 | 精品久久久久久久末码| 亚洲精品在线观看二区| 最近最新中文字幕大全电影3| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 性欧美人与动物交配| 日韩欧美三级三区| 天天一区二区日本电影三级| 欧美日韩福利视频一区二区| 午夜a级毛片| 国产成人精品久久二区二区免费| 丝袜人妻中文字幕| 人人妻人人澡欧美一区二区| 国产在线精品亚洲第一网站| 两个人视频免费观看高清| 免费在线观看影片大全网站| 蜜桃久久精品国产亚洲av| 91成年电影在线观看| 久久国产精品影院| 成人国语在线视频| АⅤ资源中文在线天堂| 午夜激情av网站| 日本撒尿小便嘘嘘汇集6| 岛国视频午夜一区免费看| 成人高潮视频无遮挡免费网站| 免费在线观看黄色视频的| 丰满的人妻完整版| 又黄又粗又硬又大视频| 国产熟女午夜一区二区三区| 男女做爰动态图高潮gif福利片| 亚洲成av人片在线播放无| 给我免费播放毛片高清在线观看| 99re在线观看精品视频| 看黄色毛片网站| 黄片大片在线免费观看| 五月伊人婷婷丁香| 国产激情偷乱视频一区二区| 久久久久国内视频| 美女免费视频网站| 国产激情偷乱视频一区二区| 18禁裸乳无遮挡免费网站照片| 在线视频色国产色| 免费一级毛片在线播放高清视频| 又大又爽又粗| 精华霜和精华液先用哪个| 波多野结衣巨乳人妻| 男人的好看免费观看在线视频 | 高潮久久久久久久久久久不卡| 人成视频在线观看免费观看| 99国产极品粉嫩在线观看| 黄色a级毛片大全视频| 久久亚洲真实| 村上凉子中文字幕在线| 欧美最黄视频在线播放免费| 免费看十八禁软件| 精华霜和精华液先用哪个| 中文字幕高清在线视频| 禁无遮挡网站| 国内揄拍国产精品人妻在线| 好看av亚洲va欧美ⅴa在| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久精品电影| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美不卡视频在线免费观看 | 1024香蕉在线观看| 欧美日韩瑟瑟在线播放|