999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

π型斷面超高斜拉橋渦振減振措施風(fēng)洞試驗研究

2018-04-24 08:07:17何旭輝王漢封劉夢婷
振動與沖擊 2018年7期
關(guān)鍵詞:風(fēng)速振動措施

李 歡, 何旭輝, 王漢封, 劉夢婷, 彭 思

(1. 中南大學(xué) 土木工程學(xué)院, 長沙 410075; 2. 中南大學(xué) 高速鐵路建造技術(shù)國家工程試驗室, 長沙 410075)

1955年第一座現(xiàn)代斜拉橋——斯曹姆松特橋建成至今,斜拉橋因其結(jié)構(gòu)簡單、受力明確、造型優(yōu)美、施工方便、跨越能力強而得到了長足的發(fā)展[1]。隨著斜拉橋跨度的進一步增加,為減輕結(jié)構(gòu)自重,降低造價,適應(yīng)復(fù)雜地區(qū)的施工條件,π型斷面主梁結(jié)構(gòu)形式在斜拉橋建造過程中得到了廣泛的應(yīng)用,特別是公路橋梁。π型斷面是滿足構(gòu)造要求最經(jīng)濟的斜拉橋主梁斷面形式之一,尤其是地形復(fù)雜施工困難的多山地區(qū)。但對于山區(qū)大跨度超高斜拉橋而言,其氣動穩(wěn)定性較差,很容易產(chǎn)生較大幅度的渦激振動。雖然渦激振動是一種限幅振動,不會導(dǎo)致橋梁發(fā)生毀滅性的破壞,但較大幅度的振動會引起結(jié)構(gòu)的疲勞損傷,縮短結(jié)構(gòu)壽命,降低橋上行人和車輛的舒適性[2-4],同時還可能引起其它災(zāi)難性的氣動不穩(wěn)定振動。π型斷面斜拉橋的渦激振動現(xiàn)象已經(jīng)引起人們的廣泛關(guān)注[5-7]。

為使π型梁在大跨度橋梁中安全使用,常常需要一定的減振措施來優(yōu)化其氣動特性,使其滿足抗風(fēng)要求。目前減振措施主要分為三類,即結(jié)構(gòu)措施、阻尼措施和氣動措施[8]。結(jié)構(gòu)措施復(fù)雜且實施難度大而較少采用;阻尼措施造價高、實際工程中維護困難也較少采用;氣動措施工作穩(wěn)定可靠、維護簡單和安裝方便已成為目前抑制渦振最常用的方法。Kubo等[9]針對簡化的π型斷面通過調(diào)整邊主梁的中心間距和橋面附屬設(shè)施的高度及邊主梁距主梁邊緣的距離提出了詳細的優(yōu)化方案,并通過識別主梁氣動阻尼和流場顯示解釋了該措施的減振機理,但無法為已完成初步設(shè)計的橋梁提供參考;Irwin[10]通過節(jié)段模型風(fēng)洞試驗和現(xiàn)場實測對比研究,采用豎向穩(wěn)定板成功抑制了塞文二橋的渦激振動;張志田等[11]采用穩(wěn)定板將某開口斷面斜拉橋的渦激振幅控制在規(guī)范允許的范圍內(nèi),但未對穩(wěn)定板的抑振規(guī)律進行系統(tǒng)研究;顏宇光等[12]研究了風(fēng)嘴、導(dǎo)流板、穩(wěn)定板、擾流板及其組合等多種氣動減振措施對開口斷面斜拉橋渦振控制效果,但無通用性的減振規(guī)律;楊光輝等[13]借助CFD(Compotational Fluid Dynamics)仿真技術(shù)分析了下中央穩(wěn)定板的減振機理,但數(shù)值模擬效果一般。由以上文獻研究結(jié)果可知,各種減振措施的減振效果參差不齊、通用性差,且減振機理尚不明確;不同橋梁使用同種減振措施抑振效果可能有較大差異;因此針對某一具體橋梁仍需進行詳細的研究,尤其是風(fēng)場環(huán)境復(fù)雜、結(jié)構(gòu)新穎的大跨度橋梁。

本文針對某擬建山區(qū)大跨度超高π型主梁斷面斜拉橋,通過節(jié)段模型試驗,分析比較了隔流板和下穩(wěn)定板等氣動措施對施工階段π型主梁斷面渦激振動性能的影響,并提出了改善渦激振動性能的抗風(fēng)措施。最后結(jié)合數(shù)值模擬的方法對渦振發(fā)生及減振措施的機理進行了初步探討。為同類型橋梁提供參考。

1 工程背景

某擬建山區(qū)大跨度超高斜拉橋為三塔雙索面疊合梁斜拉橋。橋梁跨徑布置為(249.5+550+550+249.5) m,全橋總長為1 599 m。三塔高度分別為320 m、328 m和298 m,中塔僅比法國米約大橋最高塔低15 m,屬世界第二,但該橋橋塔一次施工完成。橋面與河流平均水位垂直距離高達292.30 m。該橋主梁為π型梁,梁寬30.20 m,高3.58 m,寬高比為8.44;邊主梁高2.92 m,主梁懸挑長度2.25 m;除去斷面的懸挑長度,斷面的寬高比約為7.18;整橋布置及主梁斷面,如圖1所示。施工階段最大雙懸臂狀態(tài)塔高328 m,懸臂長度為270 m。

(a) 立面圖和平面圖(cm)

(b) 標準主梁斷面圖(mm)

2 節(jié)段模型渦激振動試驗

2.1 試驗設(shè)置

試驗在中南大學(xué)風(fēng)洞試驗室高速試驗段進行,該試驗段的幾何尺寸為長×寬×高=15.0 m×3.0 m×3.0 m,試驗風(fēng)速在0~94 m/s范圍內(nèi)連續(xù)可調(diào),湍流度小于0.3%,速度場不均性小于0.5%。試驗來流為均勻流。

根據(jù)實橋尺寸、風(fēng)洞斷面尺寸及規(guī)范要求,選取節(jié)段模型的縮尺比為1:40,模型采用加厚ABS板制作而成以保證其剛度要求,節(jié)段模型嚴格模擬了實橋的氣動外形,模型試驗參數(shù),如表1所示。節(jié)段模型彈性懸掛布置,如圖2所示。

表1 節(jié)段模型試驗參數(shù)

圖2 節(jié)段模型彈性懸掛狀態(tài)

試驗采用KEYENCE(LK-G150)激光位移計,其測量精度為0.04 mm。兩個激光位移計對稱布置在模型上下游,如圖2所示,通過同步測試獲得模型振動響應(yīng)信號。試驗采樣頻率為2 000 Hz,采樣時間為25 s。

2.2 試驗工況及減振措施

為研究主梁渦振性能及隔流板和下穩(wěn)定板兩種不同的氣動減振措施的減振效果,分別進行了不同風(fēng)攻角下原始斷面、增設(shè)隔流板和增設(shè)下穩(wěn)定板主梁斷面的渦振試驗。試驗工況,如表2所示。隔流板和下穩(wěn)定板兩種不同的氣動減振措施及安裝位置,如圖3所示。為保證試驗結(jié)果準確性,各減振措施均采用3 mm厚的ABS板由電腦雕刻而成,使其能與模型較好的契合。氣動措施幾何尺寸分別用邊主梁下翼緣寬度l=1.20 m和邊主梁高度h=2.92 m來無量綱化。

(a) 隔流板

(b) 下穩(wěn)定板

主梁斷面氣動措施尺寸風(fēng)攻角/(°)試驗阻尼比/%豎彎扭轉(zhuǎn)原始斷面-00.8600.841增設(shè)隔流板0.67l00.8730.8460.83l00.8730.8461.00l00.8730.846增設(shè)兩道下穩(wěn)定板0.48h-30.8620.8430.74h-30.8620.8431.10h-30.8620.8431.10h00.8620.8431.10h30.8620.843

3 結(jié)果及分析

3.1 無減振措施試驗結(jié)果分析

針對圖1所示,施工階段標準主梁斷面進行0°風(fēng)攻角的測試。由于該橋主梁為鋼-混疊合梁,依據(jù)規(guī)范[14]要求,試驗阻尼比的范圍為0.900%~1.100%。試驗時先采用較小的阻尼比,找到渦激共振鎖定風(fēng)速區(qū)間,然后調(diào)整豎彎阻尼比為0.860%,扭轉(zhuǎn)阻尼比為0.841%;由于橋位風(fēng)場環(huán)境復(fù)雜,橋塔高,最大懸臂長度大,因而采用了比規(guī)范稍小的阻尼比,使試驗結(jié)果偏安全。采用文獻[15]提供的計算公式將試驗結(jié)果換算到實橋,得到該橋的豎向渦振位移最大振幅修正系數(shù)為CV=1.61,扭轉(zhuǎn)渦振位移最大振幅修正系數(shù)為CT=1.85。并根據(jù)有限元模型動力計算結(jié)果由規(guī)范計算出實橋豎向渦激振動的允許振幅為[ha]=0.095 m,扭轉(zhuǎn)渦激振動的允許振幅為[θa]=0.219°。為了便于和規(guī)范限值比較,所有試驗結(jié)果均采用式(1)和式(2)來無量綱化,其中Chc為豎彎無量綱振幅,Cθc為扭轉(zhuǎn)無量綱振幅。

Chc=hc/[ha]=hc/0.095

(1)

Cθc=θc/[θa]=θc/0.219

(2)

圖4給出了換算到實橋的主梁振動位移隨風(fēng)速變化曲線。由圖4可知,施工階段0°風(fēng)攻角風(fēng)速16~23 m/s時出現(xiàn)了明顯的豎向渦激振動,最大振幅為規(guī)范限值的2.39倍;但無明顯的扭轉(zhuǎn)渦激振動。對渦振鎖定風(fēng)速下模型的位移時程曲線進行頻譜分析,如圖5所示。可見模型豎向渦振頻率與模型卓越頻率一致,均為6.280 Hz。

圖4 0°風(fēng)攻角無抗風(fēng)措施主梁渦振試驗結(jié)果

圖5 無抗風(fēng)措施主梁豎向渦激振動頻譜圖

3.2 減振措施試驗結(jié)果分析

根據(jù)表2所列減振措施在0°風(fēng)攻角下進行了三種寬度隔流板的渦振試驗,在-3°風(fēng)攻角下進行了三種長度兩道下穩(wěn)定板的渦振試驗。圖6為0°風(fēng)攻角增設(shè)不同寬度的隔流板主梁振動位移隨風(fēng)速變化曲線。由圖6可知,增設(shè)隔流板后主梁豎向渦振振幅降低了35%~50%,最大振幅仍超規(guī)范限值;隨著隔流板寬度的增大,其減振效果趨于穩(wěn)定,渦振起振風(fēng)速逐漸降低,但鎖定風(fēng)速區(qū)間大小幾乎保持不變,且對扭轉(zhuǎn)渦振幾乎沒有影響。

圖7為-3°風(fēng)攻角增設(shè)兩道不同長度的下穩(wěn)定板主梁振動位移隨風(fēng)速變化曲線。由圖7可知,增設(shè)兩道下穩(wěn)定板后,主梁豎向渦振振幅明顯減小。隨著下穩(wěn)定板長度的增加,其減振效果愈加顯著;當下穩(wěn)定板的長度增加到邊主梁高度的1.10倍時,豎向渦振振幅下降了80%左右,較好的抑制了主梁豎向渦激振動;然而隨著下穩(wěn)定板的增長,主梁扭轉(zhuǎn)渦振幅值逐漸增大,但振幅峰值仍滿足規(guī)范限值要求。

(a) 豎向振動位移幅值

(b) 扭轉(zhuǎn)振動位移幅值

(a) 豎向振動位移幅值

(b) 扭轉(zhuǎn)振動位移幅值

圖8為不同風(fēng)攻角增設(shè)兩道1.10h的下穩(wěn)定板主梁振動位移隨風(fēng)速變化曲線。由圖8可知,各攻角下主梁仍存在渦激共振現(xiàn)象,但振幅峰值均小于規(guī)范限值;雖然增設(shè)兩道1.10h的下穩(wěn)定板后主梁顫振臨界風(fēng)速由82.37 m/s下降到79.23 m/s,但仍滿足規(guī)范顫振檢驗風(fēng)速45.67 m/s的要求;因此施工階段可選取兩道1.10h的下穩(wěn)定板作為氣動減振措施。圖8表明渦激共振對風(fēng)攻角較為敏感,-3°攻角下穩(wěn)定板的減振效果最差,這為上文不同長度的下穩(wěn)定板選取-3°攻角進行風(fēng)洞試驗提供了依據(jù)。攻角的變化不改變渦振鎖定風(fēng)速區(qū)間的大小,但改變渦振起振風(fēng)速的大小,隨攻角的增大起振風(fēng)速有增大的趨勢。

(a) 豎向振動位移幅值

(b) 扭轉(zhuǎn)振動位移幅值

4 渦振機理初步探討

0°風(fēng)偏角下將π型主梁斷面簡化為二維斷面進行數(shù)值模擬已得到大家的認可[16-17],本次模擬采用類似方法。在FLUENT(16.0)CFD中使用SSTk-ω湍流模型,入口為速度邊界,出口為press-outlet邊界,橋梁表面為無滑移固體壁面邊界,計算域頂面、底面均為對稱界面。模型寬0.740 m、高0.091 m,計算域?qū)?0 m、高5 m,模型形心位于距入口3.5 m的計算域中央。計算域采用結(jié)構(gòu)化網(wǎng)格,在壁面處進行了加密處理,計算域及網(wǎng)格設(shè)置,如圖9所示。

(a) 計算域總體網(wǎng)格

(b) 近壁面處局部加密網(wǎng)格

渦振鎖定風(fēng)速下0°風(fēng)攻角時原始斷面周圍瞬時渦量如圖10所示。由圖10可知,與Shimada等[18]計算的寬高比為7的矩形斷面周圍渦量圖類似;且兩者的雷諾數(shù)分別為4.8×104和2.2×104,具有較高的可比性。該狀態(tài)下模型升力時程曲線,如圖11所示。數(shù)值計算升力系數(shù)平均值與風(fēng)洞試驗值0.09能較好地吻合。由以上對比分析結(jié)果可知,本次數(shù)值計算能夠較準確地模擬斷面周圍的擾流特征。

圖10 斷面周圍瞬時渦量圖

Fig.10 Instantaneous vorticity contours around the section

圖11 升力系數(shù)時程

各計算工況來流風(fēng)速均為渦振鎖定風(fēng)速,且采用來流風(fēng)速對流場進行無量綱化。圖12為原始斷面平均流場風(fēng)速矢量圖。由圖12可知,上表面氣流在模型迎風(fēng)邊緣分離,形成一個較小旋渦,并在模型上表面中部位置再附;下表面氣流在上游邊主梁下翼緣分離并部分附著在下游邊主梁的下翼緣上,在模型下部形成一個較大旋渦;下表面分離泡在上游邊主梁下翼緣形成并發(fā)展壯大,結(jié)合圖10可知,下表面旋渦結(jié)構(gòu)經(jīng)過下游邊主梁下翼緣進入尾渦區(qū)。

圖12 原始斷面風(fēng)速矢量圖

圖13為增設(shè)1.00l寬度的隔流板后主梁斷面平均流場風(fēng)速矢量圖。與圖12對比分析可知,增設(shè)隔流板前后斷面下表面旋渦結(jié)構(gòu)略為減小,前后邊主梁內(nèi)側(cè)均有小旋渦產(chǎn)生,下表面負壓值略有增大;上表面壓力場幾乎保持不變;雖然模型彈性懸掛時存在明顯的流固耦合作用,斷面繞流特性會有一些差異,但隔流板無法改變下表面旋渦結(jié)構(gòu)是其不能較好抑制渦振的主要原因。

圖13 增設(shè)1.00l隔流板斷面風(fēng)速矢量圖

圖14為不同風(fēng)攻角增設(shè)兩道1.10h長度的下穩(wěn)定板后主梁斷面平均流場風(fēng)速矢量圖。與圖12、圖13對比分析可知,0°風(fēng)攻角增設(shè)下穩(wěn)定板后主梁上表面旋渦結(jié)構(gòu)幾乎保持不變,下表面旋渦結(jié)構(gòu)在沿主梁下表面發(fā)展壯大的過程中被打碎為較小的旋渦,最后進入尾渦區(qū),且下表面負壓值明顯減小。雖然模型彈性懸掛時存在明顯的流固耦合作用,斷面繞流特性會有一些差異,但可以看出下穩(wěn)定板能較好的控制主梁下表面旋渦結(jié)構(gòu)的發(fā)展,有效抑制了渦振的幅值,但被打碎后的旋渦結(jié)構(gòu)不均勻的分布在主梁下表面的凹槽里,容易激發(fā)模型的扭轉(zhuǎn)渦激振動。-3°和3°攻角增設(shè)下穩(wěn)定板后主梁斷面周圍繞流形態(tài)與0°風(fēng)攻角類似,其抑振機理均為下穩(wěn)定板對主梁下表面旋渦結(jié)構(gòu)發(fā)展的抑制作用,但負攻角時下穩(wěn)定板對下表面旋渦發(fā)展的抑制作用有限,其減振效果略差。

(a) -3°攻角

(b) 0°攻角

(c) 3°攻角

5 結(jié) 論

對某擬建山區(qū)大跨度超高斜拉橋主梁節(jié)段模型進行了風(fēng)洞試驗,試驗結(jié)果表明,π型主梁斷面本身具有氣動不穩(wěn)定性,無抗風(fēng)措施時較易出現(xiàn)明顯的渦激振動現(xiàn)象。增設(shè)隔流板和下穩(wěn)定板能在不同程度上優(yōu)化主梁氣動特性,降低主梁渦激共振幅值。最后結(jié)合數(shù)值模擬,對渦振發(fā)生及減振措施的抑振機理進行了初步探討。結(jié)論如下:

(1) 增設(shè)一定寬度的隔流板雖然能在一定程度上抑制主梁渦激振動,但其減幅效果有限;設(shè)置隔流板后主梁渦振起振風(fēng)速降低,但鎖定風(fēng)速區(qū)間大小保持不變。

(2) 增設(shè)兩道一定長度的下穩(wěn)定板對豎彎渦振有明顯的抑制作用,一定長度范圍內(nèi)隨下穩(wěn)定板長度的增大減振效果愈加顯著,但較長的下穩(wěn)定板會激發(fā)扭轉(zhuǎn)渦振或放大扭轉(zhuǎn)渦振振幅;設(shè)置下穩(wěn)定板后主梁渦振起振風(fēng)速降低,鎖定風(fēng)速區(qū)間有減小的趨勢。

(3) 下穩(wěn)定板的減振效果對風(fēng)攻角較為敏感,正攻角時的減振效果優(yōu)于負攻角,且隨攻角的增大渦振起振風(fēng)速逐漸增大。

[1] 華文龍. 邊主梁斷面疊合梁斜拉橋渦振特性及減振措施研究[D]. 長沙:湖南大學(xué), 2013.

[2] FREIRE A M S, NEGRAO J H O, LOPES A V. Geometrical nonlinearities on the static analysis of highly flexible steel cable-stayed bridges[J]. Computers & Structures, 2006, 84 (31/32):2128-2140.

[3] 許福友, 丁威, 姜峰,等. 大跨度橋梁渦激振動研究進展與展望[J]. 振動與沖擊, 2010, 29(10):40-49.

XU Fuyou, DING Wei, JIANG Feng, et al. Development and prospect of study on vortex-induced vibration of long-span bridges[J]. Journal of Vibration & Shock, 2010, 29(10):40-49.

[4] 歐陽克儉, 陳政清, 韓艷,等. 橋面中央開口懸索橋渦激共振與制渦試驗研究[J]. 振動與沖擊, 2009, 28(7):199-202.

OUYANG Kejian, CHEN Zhengqing, HAN Yan, et al. Vortex-induced vibration of a suspension bridge with central-slotted box section and its control test study[J]. Journal of Vibration and Shock, 2009, 28(7):199-202.

[5] KUBO Y, SADASHIMA K, YAMAGUCHI E, et al. Improvement of aeroelastic instability of shallow π section[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2001, 98(14/15):1445-1457.

[6] 郭增偉, 趙林, 葛耀君,等. 基于橋梁斷面壓力分布統(tǒng)計特性的抑流板抑制渦振機理研究[J]. 振動與沖擊, 2012, 31(7):89-94.

GUO Zengwei, ZHAO Lin, GE Yaojun, et al. Mechanism analysis for vortex-induced vibration reduction of a flat streamlined steel box-shaped girder with airflow-suppressing board based on statistical property of surface pressure[J]. Journal of Vibration and Shock, 2012, 31(7):89-94.

[7] 錢國偉,曹豐產(chǎn),葛耀君. Ⅱ型疊合梁斜拉橋渦振性能及氣動控制措施研究[J]. 振動與沖擊, 2015, 34(2):176-181.

QIAN Guowei, CAO Fengchan, GE Yaojun. Vortex-induced vibration performance of a cable-stayed bridge with II shaped composite deck and its aerodynamic control measures[J]. Journal of Vibration and Shock, 2015, 34(2):176-181.

[8] 陳政清. 橋梁風(fēng)工程[M]. 北京: 人民交通出版社, 2005.

[9] KUBO Y, KIMURA K, SADASHIMA K, et al. Aerodynamic performance of improved shallow π shape bridge deck[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2002, 90(12/13/14/15):2113-2125.

[10] IRWIN P A. Bluff body aerodynamics in wind engineering[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96(6/7):701-712.

[11] 張志田,卿前志,肖瑋,等. 開口截面斜拉橋渦激共振風(fēng)洞試驗及減振措施研究[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版), 2011, 38(7):1-5.

ZHANG Zhitian, QING Qianzhi, XIAO Wei, et al. Vortex-induced vibration and control method for a cable-stayed bridge with open cross section[J]. Journal of Hunan University(Natural Sciences), 2011, 38(7):1-5.

[12] 顏宇光,楊詠昕,周銳. 開口斷面主梁斜拉橋的渦激共振控制試驗研究[J]. 中國科技論文, 2015, 10(7):760-764.

YAN Yuguang, YANG Yongxin, ZHOU Rui. Experimental study on vortex-induced vibration control measure for cable-stayed bridge with open sections[J]. China Sciencepaper, 2015, 10(7):760-764.

[13] 楊光輝,屈東洋,牛晉濤,等. π型截面渦激振動風(fēng)洞試驗及氣動抑制措施研究[J]. 石家莊鐵道大學(xué)學(xué)報(自然科學(xué)版), 2015, 38(1):34-39.

YANG Guanghui, QU Dongyang, NIU Jintao, et al. Researches on π-section vortex-induced vibration wind tunnel testing and aerodynamic suppression measures[J]. Jounal of Shijiazhuang Tiedao University (Natural Science), 2015, 38(1):34-39.

[14] 公路橋梁抗風(fēng)設(shè)計規(guī)范:JTG/T D60-01—2004[S]. 北京:人民交通出版社, 2004.

[15] 朱樂東. 橋梁渦激共振試驗節(jié)段模型質(zhì)量系統(tǒng)模擬與振幅修正方法[J]. 工程力學(xué), 2005, 22(5):204-208.

ZHU Ledong. Mass simulation and amplitude conversion of bridge section model test for vortex-induced oscillation[J]. Engineering Mechanics, 2005, 22(5):204-208.

[16] LEE S, LEE J S, KIM J D. Prediction of vortex-induced wind loading on long-span bridges[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1997(4):267-278.

[17] FUJIWARA A, KATAOKA H, ITO M. Numerical simulation of flow field around an oscillating bridge using finite difference method[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1993(1):567-575.

[18] SHIMADA K, ISHIHARA T. Application of a modifiedk-εmodel to the prediction aerodynamic characteristics of rectangular cross-section cylinders[J]. Journal of Fluids & Structures, 2002, 16(4):465-485.

猜你喜歡
風(fēng)速振動措施
振動的思考
放養(yǎng)雞疾病防治八措施
基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測
基于最優(yōu)TS評分和頻率匹配的江蘇近海風(fēng)速訂正
海洋通報(2020年5期)2021-01-14 09:26:54
振動與頻率
高中數(shù)學(xué)解題中構(gòu)造法的應(yīng)用措施
減少豬相互打斗的措施
中立型Emden-Fowler微分方程的振動性
夏季豬熱應(yīng)激及其防治措施
廣東飼料(2016年5期)2016-12-01 03:43:23
基于GARCH的短時風(fēng)速預(yù)測方法
主站蜘蛛池模板: 99精品国产自在现线观看| 亚洲自偷自拍另类小说| 高潮毛片无遮挡高清视频播放| 91免费精品国偷自产在线在线| 亚洲人成色在线观看| 亚洲手机在线| 91精品国产91久久久久久三级| 美女国内精品自产拍在线播放| 成人免费网站久久久| 综合人妻久久一区二区精品 | 日韩国产亚洲一区二区在线观看| 国产成人夜色91| 欧美一级99在线观看国产| 97在线公开视频| 香蕉久久永久视频| 在线高清亚洲精品二区| 在线亚洲小视频| 亚洲欧美国产五月天综合| 亚洲国产成人无码AV在线影院L| 伊人久久综在合线亚洲2019| 日本少妇又色又爽又高潮| 国产青榴视频| 黄色污网站在线观看| 亚洲中文字幕97久久精品少妇| 国产一级毛片在线| 亚洲一级无毛片无码在线免费视频| 色婷婷视频在线| 天天摸夜夜操| 色天天综合| 久草青青在线视频| 精品人妻系列无码专区久久| 国产一级妓女av网站| 99久久人妻精品免费二区| 五月六月伊人狠狠丁香网| 久久人搡人人玩人妻精品一| 久久婷婷色综合老司机| 永久在线精品免费视频观看| 久久亚洲美女精品国产精品| 亚洲熟女偷拍| 久久96热在精品国产高清| 久久永久免费人妻精品| 久久久波多野结衣av一区二区| 国内毛片视频| 欧美精品啪啪| 国产女人爽到高潮的免费视频| 亚洲午夜片| 久久亚洲国产最新网站| 99资源在线| 99久久国产综合精品女同| 亚洲天堂.com| 亚洲精品午夜无码电影网| 国产丰满大乳无码免费播放| 亚洲精品国产日韩无码AV永久免费网| 国产黄在线观看| 久久99精品久久久久久不卡| 91精品国产综合久久不国产大片| 亚洲性色永久网址| 91麻豆久久久| 久久99精品久久久大学生| 91精品视频在线播放| 欧美区日韩区| 色色中文字幕| 国产人在线成免费视频| 国产成人综合日韩精品无码不卡 | 久久精品丝袜| 亚洲成人在线网| 成人福利在线视频| 美女被操黄色视频网站| 亚洲男人在线天堂| 国产成人综合网| 伊人久热这里只有精品视频99| 女人av社区男人的天堂| 亚洲色图综合在线| 呦视频在线一区二区三区| 亚洲成AV人手机在线观看网站| 国产精品免费电影| 国产精品吹潮在线观看中文| 黄色片中文字幕| 狠狠v日韩v欧美v| 无码精品国产VA在线观看DVD| 国产网站免费观看| 亚洲欧美在线综合图区|