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

    車載作用下公路橋梁耦合振動精細化建模及驗證分析

    2021-09-29 18:03殷新鋒晏萬里任厚乾劉揚
    湖南大學學報·自然科學版 2021年9期
    關(guān)鍵詞:響應面法

    殷新鋒 晏萬里 任厚乾 劉揚

    摘? ?要:現(xiàn)有車-橋耦合振動分析中車輛模型不能精確考慮車輛動力特性和柔性輪胎對車橋耦合振動響應的影響.為了進一步研究充氣輪胎胎壓對車-橋耦合振動的影響,基于LS-DYNA程序,采用線彈性橡膠材料模擬輪胎并定義輪胎內(nèi)氣壓,結(jié)合常用重載三軸汽車的結(jié)構(gòu)參數(shù),運用彈簧阻尼單元及梁、殼單元模擬車輛懸架系統(tǒng)的動力特性,建立可分析車輪氣壓的三維車輛模型;并基于實橋試驗結(jié)果及響應面法得到高精度有限元橋梁模型;通過顯式求解程序LS-DYNA內(nèi)置的接觸算法,將車輛子系統(tǒng)和橋梁子系統(tǒng)聯(lián)立耦合起來,形成顯式的車-橋耦合振動分析模型.計算結(jié)果與實測結(jié)果對比分析驗證了該方法的正確性,并分析了輪胎胎壓對橋梁振動的影響.

    關(guān)鍵詞:響應面法;有限元模型修正;連續(xù)剛構(gòu)橋;LS-DYNA;動力響應

    中圖分類號:U441.2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼:A

    Fine Modeling of Coupled Vibration of Highway Bridge under

    Vehicle Loading and Verification Analysis

    YIN Xinfeng1,YAN Wanli1,REN Houqian1,LIU Yang1,2

    (1. College of Civil Engineering,Changsha University of Science & Technology,Changsha 410114,China;

    2. College of Civil Engineering,Hunan University of Technology,Zhuzhou 412007,China)

    Abstract:In the existing coupling vibration analysis of vehicle bridge,the influence of vehicle dynamic characteristics and flexible tire on the vehicle bridge coupling vibration response can not be accurately considered. In order to further study the influence of pneumatic tire pressure on vehicle bridge coupling vibration,the linear elastic rubber material is used to simulate the tire and define the air pressure in the tire based on LS-DYNA program. Combined with the structural parameters of the commonly used heavy-duty three-axle vehicle,the dynamic characteristics of the vehicle suspension system are simulated by using the spring damping element,beam and shell element,and a three-dimensional vehicle model is established. Based on the test results of the real bridge and the response surface method,a high-precision finite element bridge model is obtained. The vehicle subsystem and the bridge subsystem are coupled together to form the explicit vehicle bridge coupling vibration analysis model through the contact algorithm built in the explicit solution program LS-DYNA. The correctness of the method is verified by comparing the calculation results with the measured results,and the influence of tire pressure on bridge vibration is analyzed.

    Key words:response surface method;finite element model updating;continuous rigid-frame bridge;LS-DYNA;dynamic response

    車輛在橋梁上行駛時產(chǎn)生的激勵導致車輛和橋梁產(chǎn)生相互振動,該振動稱為車橋耦合振動[1]. 國內(nèi)外學者對車橋耦合振動問題進行了大量研究,并取得非凡成果.夏禾[2]將車輛簡化為懸掛振動模型,研究了車-橋-墩相互作用的動力響應. Cai等[3-4]采用兩軸車輛模型,基于功率譜密度函數(shù)生成隨機橋面不平整度,分析了車輛荷載作用下橋梁振動響應. Huang等[5]采用三軸車輛有限元模型研究了簡支梁的車橋隨機振動. 韓萬水等[6]結(jié)合實測數(shù)據(jù)對橋梁模型進行修正以得到橋梁基準模型,并采用梁格法對車橋振動進行分析. 李奇等[7]考慮車體柔性的影響,分析了高速列車和簡支梁橋相互作用的影響. 鄧露等[8]運用動態(tài)稱重技術(shù)(BWIM)實時監(jiān)測車輛總重和軸重,進行中小型跨徑混凝土梁橋的振動研究.趙越等[9]基于等參映射與改進折半法在傳統(tǒng)車橋耦合分析的基礎(chǔ)上進一步提升其分析精度及計算效率,從而進行公路車橋耦合分析.但關(guān)于充氣輪胎胎壓對車-橋耦合振動影響的研究較少.主要原因為現(xiàn)有車-橋耦合振動分析中車輛模型多為簡化的質(zhì)量-彈簧-阻尼多自由度振動模型[10-13],該模型常將車輪簡化為點或者等效線面接觸,故不能精確考慮車輛動力特性和柔性輪胎對車橋耦合振動響應的影響.因此,建立精確車輛模型和橋梁模型是至關(guān)重要的,這樣才能提高數(shù)值模擬精度,為橋梁結(jié)構(gòu)的健康運營提供有利建議.

    本文以一座連續(xù)剛構(gòu)箱梁橋為工程背景.首先,基于LS-DYNA程序以車輛實際構(gòu)造及動力特性為基準,建立車輛精細三維有限元模型;然后,結(jié)合響應面法并依據(jù)實橋試驗結(jié)果對橋梁模型進行修正,以得高精度的橋梁有限元模型;最后,通過LS-DYNA程序?qū)蛄耗P秃蛙囕v模型進行耦合,求解車橋的振動響應,并通過改變車輪氣壓,分析該參數(shù)對橋梁振動響應的影響.

    1? ?車輛模型的建立及校驗

    1.1? ?車輛模型的建立

    在車橋耦合振動模型中,影響橋梁振動響應的主要因素為車輛模型的動力特性和荷載分布. 車輛模型需著重模擬懸架系統(tǒng)、車輪和軸重. 本文參照東風牌三軸載重自卸貨車,建立車輛有限元模型,其相關(guān)參數(shù)如表1所示.前中軸距為3.5 m,中后軸距為1.4 m,后輪距為1.8 m. 基于CAD軟件SolidWorks分別建立各部位的幾何模型并進行網(wǎng)格劃分,最后使用梁、殼、實體單元及離散單元等賦予不同部位網(wǎng)格不同的屬性,從而構(gòu)建車輛模型,如圖1所示.

    車輪由輪胎、輪盤和輪轂組成.輪胎采用線彈性橡膠材料,輪盤和輪轂采用線彈性鋼材材料,車輪模型中各部位連接均為剛性連接,邊界為剛性固態(tài)約束,有限元模型如圖2所示. 根據(jù)輪胎內(nèi)實際氣壓,使用LS-DYNA程序關(guān)鍵字*AIRBAG_SIMPLE_AIRBAG_MODEL定義由輪胎、輪轂組成的封閉體內(nèi)的氣壓,其氣體壓力值為0.6 MPa.

    由于車輛懸架系統(tǒng)結(jié)構(gòu)復雜,本文采用殼單元、剛體、彈簧阻尼單元和多點約束來模擬懸架系統(tǒng).前后懸架的彈簧和阻尼器的參數(shù)參考文獻[14],前懸架的彈簧剛度取800 N/mm,阻尼系數(shù)取20 Ns/mm;后懸架的彈簧剛度取1 200 N/mm,阻尼系數(shù)取25 Ns/mm. 前后懸架實體構(gòu)造及有限元模型如圖3所示.

    為了使車輪正常轉(zhuǎn)動,在輪盤與車軸之間設置旋轉(zhuǎn)鉸,再通過定義*INITIAL_VELOCTIY_GENERATION關(guān)鍵字設置車輪的轉(zhuǎn)動和平動速度,實現(xiàn)車輪滾動向前的狀態(tài),如圖4所示.

    1.2? ?車輛模型的校驗

    在車輛模型的3個軸上選擇6個節(jié)點,約束豎向位移. 然后對車輛施加重力荷載使車輛產(chǎn)生瞬時振動,再通過定義關(guān)鍵字*DAMPING_GLOBAL對車輛模型施加全局阻尼,使其快速達到穩(wěn)定狀態(tài),計算出約束反力,并與實測車輛軸重進行比較,最終得到車輛模型的軸載曲線如圖5所示.

    由表2知,車體總重偏差為0.87%,說明實測車輛軸重分布特性能體現(xiàn)于建立的車輛模型.

    通過對比車輛模型和實測車輛的自振頻率,驗證車輛模型動力特性的有效性. 如圖6所示,該車實測自振頻率為1.635 Hz,車輛有限元模型的自振頻率為1.647 Hz,顯然,兩者基頻相差很小.

    2? ?橋梁模型

    2.1? ?橋梁概況

    該橋主橋為三跨變截面預應力混凝土連續(xù)剛構(gòu)結(jié)構(gòu),跨徑布置為(65+120+65) m,其立面及斷面示意圖分別如圖7和圖8所示.

    2.2? ?實橋試驗

    利用有限元軟件ANSYS建立初始有限元模型如圖9所示. 主梁和墩身都用實體單元模擬,墩底采用固結(jié)連接. 混凝土密度取為2 500 kg/m3,主梁和墩的混凝土彈性模量分別為34.5 GPa和32.5 GPa.

    現(xiàn)場對橋梁進行靜力試驗,測試中加載車輛總數(shù)為6輛,每輛總重為35 t,車輛前軸重7 t,中軸和后軸各重14 t,車輛照片如圖10所示.

    在正式試驗前先進行預加載,以消除非彈性變形,確保試驗及設備處于良好工作狀態(tài). 正式試驗中將加載車分3級加載,每一級加載持荷2 ~ 3 min,待實測應變及撓度數(shù)據(jù)穩(wěn)定后進行數(shù)據(jù)的采集工作,再進行下一級的加載工作. 限于篇幅,僅就其中一個工況做詳要概述. 為了使加載截面承受最大正彎矩,在中跨跨中截面中心加載,相應的荷載效率系數(shù)為0.967.撓度測點布置如圖11所示.

    實橋試驗的撓度采用水準儀進行測量,現(xiàn)場測試照片如圖12所示.脈動試驗采用多通道數(shù)據(jù)采集分析系統(tǒng)NI公司的PXI系統(tǒng)進行試驗,采用8330B3型超低頻加速度傳感器進行數(shù)據(jù)采集,最低采樣頻率從0 Hz開始,采集主梁橫向和豎向振動數(shù)據(jù),再經(jīng)信號分析得到全橋的各階固有振動特性,現(xiàn)場采集照片如圖13所示.

    2.3? ?基于響應面法的有限元模型更新

    基于響應面法的有限元模型修正是用響應面函數(shù)來模擬實際結(jié)構(gòu)的響應函數(shù),將試驗設計與數(shù)理統(tǒng)計相結(jié)合,通過樣本選取、方差分析參數(shù)選取、響應面的擬合及采用優(yōu)化算法尋求響應面模型中的最優(yōu)解來進行有限元模型修正.

    據(jù)文獻[15-16]知,影響有限元模型與實測橋梁結(jié)構(gòu)差別的主要因素為結(jié)構(gòu)混凝土密度、主梁和墩混凝土彈性模量,因此取這三個因素作為變量. 為獲得響應與所選定的三參數(shù)之間的聯(lián)系,首先需要參數(shù)設計. 根據(jù)參數(shù)取值的變化規(guī)律,假設三參數(shù)的單位長度值為10%,則可得該三參數(shù)值的變化范圍見表3.

    以橋梁模態(tài)和靜力變形為目標來更新橋梁模型,選取橋梁結(jié)構(gòu)的第一階自振頻率(R1)、測點5撓度值(R2)和測點6撓度值(R3)作為目標函數(shù).這3個目標函數(shù)充分利用了現(xiàn)場的實橋試驗條件,且包含了橋梁靜、動力性能指標,可較為全面、準確地反映橋梁的力學性能.

    根據(jù)上述試驗設計,應用回歸分析技術(shù)對樣本數(shù)據(jù)進行響應面擬合,并用R2準則和R2adj準則進行響應面擬合精度的檢驗,其表達式見方程(1)(2).

    式中:SST = SSE + SSR表示模型的總方差;dT表示模型的總自由度.如果R2和R2adj的值都接近1且兩者差值很小,則表示響應面方程擬合得很好.

    采用二次多項式對上述樣本結(jié)果進行響應面擬合,并用R2和R2adj準則對響應面的精度進行檢驗,檢驗結(jié)果如表4所示.

    從表4知,R2和R2adj值均接近1,且兩者差值很小,說明模型擬合精度較高,能很好地表示各參數(shù)與目標函數(shù)的關(guān)系.

    響應面函數(shù)采取省略交差項影響的二次多項式,則響應面擬合方程如方程(3)~(5)所示:

    R1 = 0.821 9 - 0.000 3X1 + 0.012 7X2 +

    0.012 0X3 - 0.000 1X22 - 0.000 1X23? ? ? ?(3)

    R2 = 9.764 8 + 0.005 0X1 - 0.632 7X2 +

    0.067 9X3 + 0.008 7X22 - 0.001 5X23? ? ? ?(4)

    R3 = 16.557 2 + 0.008 6X1 - 1.050 3X2 +

    0.091 6X3 + 0.014 4X22 - 0.002 3X23? ? ? ?(5)

    該橋的第一階實測頻率為0.867 Hz;實橋跨中試驗加載的情形下,測點5、6的撓度值分別為5.687 mm和9.679 mm. 根據(jù)實測值和數(shù)值模擬值可以建立一個目標函數(shù):

    式中:R(i)為根據(jù)響應面法計算的響應;M(i)為實測值;coef(i)為權(quán)重系數(shù),此處取為[1.5 1 1].目標函數(shù)可以利用Matlab中的遺傳算法(Genetic Algorithm)工具箱優(yōu)化.更新后的上述三參數(shù)分別如表5所示.

    從表5可得,主梁和墩的混凝土彈性模量和密度都有所增加,說明依托工程的施工質(zhì)量較為可靠,同時參數(shù)的變動范圍仍在工程實踐的常見范圍之內(nèi),說明了修正結(jié)構(gòu)的可靠性.利用所更新的參數(shù)計算得到大橋的第一階振動頻率及上述工況的試驗加載測點撓度如表6所示.

    3? ?振動響應對比分析

    3.1? ?現(xiàn)場測點布置及荷載工況

    該連續(xù)剛構(gòu)橋主要對主箱梁在如圖14所示截面1、2處(即邊跨跨中、中跨跨中)的動位移與應變變化情況進行了分析,全橋共計布設動位移測點5個,應變測點18個.

    現(xiàn)場布置了3種荷載工況:工況1,一輛35 t重車置于橋縱向中軸線上;工況2,兩輛35 t重車并排對稱置于橋縱向中軸線兩側(cè),兩車橫向間距為1.3 m;工況3,四輛35 t重車以橋縱向中軸線對稱放置,車輛橫向間距為1.3 m,縱向間距為4 m. 所有工況,車輛均以40 km/h的速度從橋外駛?cè)氩⑼ㄟ^橋梁.限于篇幅,僅列出工況3的車輛布置圖,如圖15所示.

    3.2? ?結(jié)果對比分析

    3.2.1? ?動位移對比

    將上述各個工況進行對比分析,邊/中跨跨中的位移實測值與計算值對比如圖16和圖17所示.

    由圖16和圖17可知,計算值與實測值隨時間變化的波動規(guī)律一致,但存在同一時刻計算值與實測值吻合不是很好的情況,且計算值的曲線波峰稍大于實測值,兩者靜態(tài)豎向位移差值在1.71%~6.55%之間,振動位移(消除靜態(tài)位移值)幅值差值在3.71%~6.35%之間.其誤差主要原因可能是:1)橋梁模型雖經(jīng)過合理的參數(shù)更新,但與實測橋梁結(jié)構(gòu)特性還有一定差別;2)實橋試驗受車輛速度控制、行駛軌跡、儀器誤差等因素的影響,測量數(shù)據(jù)精度有較大的誤差.

    3.2.2? ?動應變對比

    將上述各個工況進行對比分析,邊、中跨跨中的動應變實測值與計算值對比如圖18和圖19所示.

    由圖18和圖19可知,計算值與實測值兩者的變化規(guī)律一致,且數(shù)值模擬結(jié)果與實測數(shù)據(jù)吻合得較好.但計算值的曲線波峰稍大于實測值,兩者動應變的差值在3.35%~6.94%之間,振動應變(消除靜態(tài)應變值)幅值差值在2.83%~5.96%.

    從上述對比情況可知,豎向振動位移與動應變的計算值與實測值隨時間的變化趨勢一致,且整體吻合較好,峰值接近,說明該法建立的橋梁模型和車輛模型是有效的,能應用于車橋耦合計算中.

    4? ?輪胎氣壓對橋梁振動的影響

    輪胎氣壓的大小與輪胎變形成幾何關(guān)系,因此,輪胎氣壓對跨中振動響應影響反映了輪胎變形對梁跨中振動的影響.以上述工況2為例,取3種分別為P1=0.4 MPa、P2=0.6 MPa、P3=0.8 MPa的不同氣壓,來分析車輪氣壓變化對中跨跨中動位移和動應變的影響.

    從圖20和圖21可得,當車輪氣壓由0.6 MPa減少到0.4 MPa時,車輪柔性變大,與橋面的接觸面增大,橋梁中跨跨中動位移峰值減小3.6%,動應變減小2.7%;當車輪氣壓由0.6 MPa增加到0.8 MPa時,車輪剛度變大,與橋面的接觸面減少,橋梁中跨跨中動位移峰值增大4.4%,動應變增大3.85%.

    5? ?結(jié)? ?論

    本文基于LS-DYNA程序建立了車-橋相互作用模型.設立多種荷載工況,通過計算值與實測值的對比,驗證了車輛模型的有效性及該法建立車橋耦合模型的可行性,并進一步分析了車輪氣壓變化對車橋振動的影響.結(jié)果表明:1)本文方法可較真實再現(xiàn)車-橋耦合振動情形,無需自編復雜的車橋耦合分析程序,只需定義LS-DYNA軟件自帶的面面接觸關(guān)鍵字,就可以實現(xiàn)車橋間的耦合作用;2)可方便獲得橋梁任意部位的應力、應變及其他響應數(shù)據(jù);3)可高精度分析車輪氣壓變化及車輛其他參數(shù)對耦合系統(tǒng)振動影響.

    參考文獻

    [1]? ? 殷新鋒. 汽車荷載作用下梁式橋與斜拉橋的動態(tài)響應分析[D].長沙:湖南大學,2010:6—20.

    YIN X F. Dynamic response analysis of girder bridge and cable-stayed bridge under vehicle load [D]. Changsha:Hunan University,2010:6—20. (In Chinese)

    [2]? ? 夏禾. 車輛與結(jié)構(gòu)動力相互作用[M]. 北京:科學出版社,2005:154—164.

    XIA H. Vehicle-structure interaction [M]. Beijing:Science Press,2005:154—164. (In Chinese)

    [3]? ? CHEN S R,CAI C S. Accident assessment of vehicles on long-span bridges inwindy environments[J]. Journal of Wind Engineering and Industrial Aerodynamics,2004,92(12):991—1024.

    [4]? ? CAI C S,SHI X M,ARAUJO M. Effect of approach span condition on vehicle-induced dynamic response [J]. ASCE Journal of Bride Engineering,2007,14(4):440—449.

    [5]? ? HUANG D Z,WANG T L. Impact analysis of cable-stayed bridges[J]. Computers and Structures,1992,43(5):897—908.

    [6]? ? 韓萬水,王濤,李永慶,等. 基于模型修正梁格法的車橋耦合振動分析系統(tǒng)[J]. 中國公路學報,2011,24(5):47—55.

    HAN W S,WANG T,LI Y Q,et al. Analysis system of vehicle-bridge coupling vibration with grillage method based on model updating[J]. China Journal of Highway and Transport,2011,24(5):47—55. (In Chinese)

    [7]? ? 李奇,吳定俊,邵長宇. 考慮車體柔性的車橋耦合系統(tǒng)建模與分析方法[J]. 振動工程學報,2011,24(1):41—47.

    LI Q,WU D J,SHAO C Y. Modeling and dynamic analysis method of vehicle-bridge coupling system considering car-body flexibility[J]. Journal of Vibration Engineering,2011,24(1):41—47. (In Chinese)

    [8]? ? 鄧露,李樹征,淡丹輝,等. 橋梁動態(tài)稱重技術(shù)在中小跨徑混凝土梁橋上的適用性研究[J]. 湖南大學學報(自然科學版),2020,47(3):89—96.

    DENG L,LI S Z,DAN D H,et al. Study on applicability of bridge weigh-in-motion technology in short-to medium-span concrete girder bridges[J].Journal of Hunan University (Natural Sciences),2020,47(3):89—96. (In Chinese)

    [9]? ? 趙越,黃平明,劉修平,等. 基于等參映射與改進折半法的公路車橋耦合分析系統(tǒng)[J]. 湖南大學學報(自然科學版),2020,47(11):29—37.

    ZHAO Y,HUANG P M,LIU X P,et al. Highway vehicle-bridge interaction analysis system based on isoparametric mapping and improved binary search [J]. Journal of Hunan University (Natural Sciences),2020,47 (11):29—37. (In Chinese)

    [10]? 李小珍,劉德軍,晉智斌. 大跨度鐵路懸索橋車-線-橋耦合振動分析[J]. 鋼結(jié)構(gòu),2010,25(12):6—12.

    LI X Z,LIU D J,JIN Z B. Analysis of train-track-bridge coupled vibration of a railway long-span suspension bridge[J]. Steel Construction,2010,25(12):6—12. (In Chinese)

    [11]? 朱志輝,王力東,龔威,等. 多種垂向輪軌關(guān)系的對比及改進的車-線-橋系統(tǒng)迭代模型的建立[J].中南大學學報(自然科學版),2017,48(6):1585—1593.

    ZHU Z H,WANG L D,GONG W,et al. Comparative analysis of several types of vertical wheel/rail relationship and construction of an improved iteration model for train-track-bridge system[J]. Journal of Central South University (Science and Technology),2017,48(6):1585—1593.(In Chinese)

    [12]? 譚國金,劉寒冰,程永春,等. 基于車-橋耦合振動的簡支梁橋沖擊效應[J]. 吉林大學學報(工學版),2011,41(1):62—67.

    TAN G J,LIU H B,CHENG Y C,et al. Analysis of impact of vehicle to simply supported beam bridge based on vehicle-bridge coupled vibration[J]. Journal of Jilin University (Engineering and Technology Edition),2011,41(1):62—67. (In Chinese)

    [13]? 李巖,吳志文,蔡明,等. 一種隨機車流與橋梁耦合振動的分析方法[J]. 哈爾濱工業(yè)大學學報,2018,50(3):46—52.

    LI Y,WU Z W,CAI M,et al. An analysis method for coupled vibration random traffic flow and bridge[J]. Journal of Harbin Institute of Technology,2018,50(3):46—52. (In Chinese)

    [14]? 陳水生,趙輝,桂水榮,等. 基于LS -DYNA公路橋車橋耦合的車輛模型研究[J]. 計算力學學報,2019,36(6):747—756.

    CHEN S S,ZHAO H,GUI S R,et al. Vehicle model research on vehicle bridge coupling of highway bridge based on LS -DYNA[J]. Chinese Journal of Computational Mechanics,2019,36(6):747—756. (In Chinese)

    [15]? DENG L. System identification of bridge and vehicle based on their coupled vibration[D]. Baton Rouge:Louisiana State University,2009:7—23.

    [16]? YIN X F,F(xiàn)ANG Z,CAI C S,et al. Non-stationary random vibration of bridges under vehicles with variable speed[J]. Engineering Structures,2010,32(8):2166—2174.

    猜你喜歡
    響應面法
    響應面法優(yōu)化超聲波輔助提取白頭翁總皂苷
    菌糠中產(chǎn)纖維素酶菌株的篩選、鑒定及增殖條件響應面法優(yōu)化
    玉米皮類黃酮超聲波提取工藝優(yōu)化及其抗氧化活性研究
    基于響應面法的聽小骨消聲器的優(yōu)化設計
    山藥多糖提取工藝優(yōu)化及其抗菌活性研究
    微波輔助提取白頭翁皂苷研究
    響應面法優(yōu)化超聲輔助提取蕎麥中蘆丁的工藝
    菊芋膳食纖維酸奶的工藝研究及營養(yǎng)分析
    船舶舷側(cè)構(gòu)件幾何尺寸優(yōu)化分析
    可溶性大豆多糖超聲波提取工藝及其抗氧化性研究
    国产熟女午夜一区二区三区| 亚洲成av人片免费观看| 美女国产高潮福利片在线看| 欧美另类亚洲清纯唯美| 两性夫妻黄色片| 欧美午夜高清在线| 在线观看舔阴道视频| 大码成人一级视频| 亚洲av电影在线进入| 男女午夜视频在线观看| e午夜精品久久久久久久| 国产精品精品国产色婷婷| 九色亚洲精品在线播放| www国产在线视频色| 中国美女看黄片| 亚洲精品美女久久久久99蜜臀| 日本欧美视频一区| 伊人久久大香线蕉亚洲五| 国产午夜福利久久久久久| 老司机福利观看| 丝袜人妻中文字幕| 日本a在线网址| 深夜精品福利| 婷婷六月久久综合丁香| 欧美乱码精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 久久久国产欧美日韩av| 久久精品91蜜桃| 亚洲电影在线观看av| 欧美绝顶高潮抽搐喷水| 久久人人爽av亚洲精品天堂| 人人妻人人澡人人看| 色综合站精品国产| 免费av毛片视频| 精品福利观看| 精品午夜福利视频在线观看一区| 桃色一区二区三区在线观看| 香蕉国产在线看| 免费人成视频x8x8入口观看| 88av欧美| 变态另类成人亚洲欧美熟女 | 一边摸一边做爽爽视频免费| 国内久久婷婷六月综合欲色啪| 国产精品自产拍在线观看55亚洲| 乱人伦中国视频| 精品久久蜜臀av无| 午夜福利,免费看| 精品久久久久久久人妻蜜臀av | 一个人观看的视频www高清免费观看 | 中文字幕高清在线视频| 亚洲人成伊人成综合网2020| 国产色视频综合| 国产亚洲av嫩草精品影院| 欧美色视频一区免费| 亚洲人成伊人成综合网2020| 男人的好看免费观看在线视频 | 搡老熟女国产l中国老女人| 国产亚洲精品av在线| 欧美中文日本在线观看视频| 18禁黄网站禁片午夜丰满| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 99精品在免费线老司机午夜| 黄色毛片三级朝国网站| 日日夜夜操网爽| 一级毛片高清免费大全| 啦啦啦观看免费观看视频高清 | 免费在线观看影片大全网站| 老汉色∧v一级毛片| 中文亚洲av片在线观看爽| 久久国产精品影院| av免费在线观看网站| 欧美精品啪啪一区二区三区| 中文字幕av电影在线播放| 亚洲精品国产色婷婷电影| 日日爽夜夜爽网站| 97人妻天天添夜夜摸| 大型黄色视频在线免费观看| 啦啦啦免费观看视频1| av视频免费观看在线观看| 国产高清激情床上av| 亚洲片人在线观看| 18禁国产床啪视频网站| 久久久久久国产a免费观看| 国产又爽黄色视频| 乱人伦中国视频| 啦啦啦免费观看视频1| 亚洲美女黄片视频| 成人国产一区最新在线观看| 一进一出抽搐gif免费好疼| 国产av在哪里看| 久久精品国产清高在天天线| 久久久久久大精品| 男女下面插进去视频免费观看| 极品教师在线免费播放| 国产精品久久视频播放| 夜夜夜夜夜久久久久| 我的亚洲天堂| 男男h啪啪无遮挡| 国产精品一区二区免费欧美| 韩国av一区二区三区四区| 黄网站色视频无遮挡免费观看| 日韩精品青青久久久久久| 亚洲av美国av| 国产精品av久久久久免费| 国产av精品麻豆| 亚洲精品美女久久久久99蜜臀| 久久久精品欧美日韩精品| 俄罗斯特黄特色一大片| 亚洲精品久久国产高清桃花| 非洲黑人性xxxx精品又粗又长| 操美女的视频在线观看| 制服人妻中文乱码| 多毛熟女@视频| 国产极品粉嫩免费观看在线| 91av网站免费观看| 久久精品国产亚洲av香蕉五月| 欧美一区二区精品小视频在线| 老司机在亚洲福利影院| 久久 成人 亚洲| 欧美久久黑人一区二区| 天堂影院成人在线观看| 国产精品一区二区免费欧美| 淫妇啪啪啪对白视频| 久久久久国产一级毛片高清牌| 久久这里只有精品19| 欧美日韩一级在线毛片| 一进一出抽搐动态| 国产精华一区二区三区| 久久亚洲真实| 性少妇av在线| 国产精品 欧美亚洲| 欧美大码av| 黄色a级毛片大全视频| 美女扒开内裤让男人捅视频| 国产高清videossex| 少妇的丰满在线观看| 色综合亚洲欧美另类图片| av在线播放免费不卡| 亚洲男人天堂网一区| 国产精品精品国产色婷婷| 久久九九热精品免费| 在线观看免费日韩欧美大片| 久99久视频精品免费| 国产亚洲av嫩草精品影院| 精品一区二区三区四区五区乱码| 两人在一起打扑克的视频| 国产精品美女特级片免费视频播放器 | 久久久精品国产亚洲av高清涩受| 黄色视频不卡| 乱人伦中国视频| 亚洲国产欧美日韩在线播放| 国产精品乱码一区二三区的特点 | 12—13女人毛片做爰片一| 天天一区二区日本电影三级 | 国产真人三级小视频在线观看| 麻豆成人av在线观看| 成人特级黄色片久久久久久久| 正在播放国产对白刺激| 97超级碰碰碰精品色视频在线观看| 色播亚洲综合网| 51午夜福利影视在线观看| 91麻豆精品激情在线观看国产| 美女大奶头视频| 亚洲少妇的诱惑av| 亚洲精品在线美女| 久久 成人 亚洲| 大码成人一级视频| 一区二区三区精品91| 狂野欧美激情性xxxx| 欧美乱码精品一区二区三区| 黄色片一级片一级黄色片| 90打野战视频偷拍视频| 国产精品久久久人人做人人爽| 不卡一级毛片| 午夜激情av网站| 精品福利观看| 国产免费男女视频| 国产高清videossex| 一个人免费在线观看的高清视频| 精品国产乱码久久久久久男人| 最近最新免费中文字幕在线| 国产欧美日韩综合在线一区二区| 欧美日韩亚洲国产一区二区在线观看| 男男h啪啪无遮挡| 久久久久国产精品人妻aⅴ院| 久久人人97超碰香蕉20202| 国产精品久久久人人做人人爽| 国产主播在线观看一区二区| 黑丝袜美女国产一区| 搡老妇女老女人老熟妇| av福利片在线| 成熟少妇高潮喷水视频| 看片在线看免费视频| 丝袜美足系列| 国产伦一二天堂av在线观看| 国产亚洲精品一区二区www| 欧美日韩瑟瑟在线播放| 精品一区二区三区视频在线观看免费| 国产真人三级小视频在线观看| 久久婷婷人人爽人人干人人爱 | 天天添夜夜摸| 性欧美人与动物交配| 黄色视频不卡| 国产欧美日韩精品亚洲av| 国产一区二区三区综合在线观看| 视频区欧美日本亚洲| 国产午夜精品久久久久久| 国产精品av久久久久免费| 高清毛片免费观看视频网站| 狂野欧美激情性xxxx| 亚洲精品国产色婷婷电影| 国产av一区在线观看免费| 夜夜躁狠狠躁天天躁| 欧美精品啪啪一区二区三区| 亚洲国产精品合色在线| 日本免费一区二区三区高清不卡 | 欧美成人一区二区免费高清观看 | 久久亚洲真实| 国产成人影院久久av| 欧美色欧美亚洲另类二区 | 久久精品亚洲熟妇少妇任你| 国产视频一区二区在线看| 国产成人系列免费观看| 自线自在国产av| 亚洲第一青青草原| 在线免费观看的www视频| av在线天堂中文字幕| 男人的好看免费观看在线视频 | 久久婷婷人人爽人人干人人爱 | 黄色毛片三级朝国网站| 国产精品98久久久久久宅男小说| 人人妻人人爽人人添夜夜欢视频| 日韩精品免费视频一区二区三区| 90打野战视频偷拍视频| 91国产中文字幕| 亚洲天堂国产精品一区在线| 女人爽到高潮嗷嗷叫在线视频| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 在线免费观看的www视频| 国产精品av久久久久免费| 中亚洲国语对白在线视频| 可以在线观看的亚洲视频| 久久草成人影院| 色播在线永久视频| 一进一出抽搐gif免费好疼| 国产精品免费一区二区三区在线| 亚洲精品中文字幕一二三四区| aaaaa片日本免费| 黄片播放在线免费| 侵犯人妻中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 搡老岳熟女国产| 女性生殖器流出的白浆| cao死你这个sao货| 精品欧美一区二区三区在线| 中文字幕高清在线视频| 又紧又爽又黄一区二区| 欧美一区二区精品小视频在线| 他把我摸到了高潮在线观看| 国产精品 国内视频| 色av中文字幕| 一本大道久久a久久精品| 日韩欧美国产一区二区入口| 88av欧美| 18美女黄网站色大片免费观看| 国产单亲对白刺激| 亚洲精品av麻豆狂野| 亚洲三区欧美一区| 男人舔女人的私密视频| 久久精品国产亚洲av香蕉五月| 午夜福利免费观看在线| 欧美日韩福利视频一区二区| 母亲3免费完整高清在线观看| 黑丝袜美女国产一区| 久久性视频一级片| 国产激情久久老熟女| 国产亚洲精品一区二区www| 亚洲成人免费电影在线观看| 给我免费播放毛片高清在线观看| 久久精品成人免费网站| av片东京热男人的天堂| 在线播放国产精品三级| 日韩国内少妇激情av| 免费在线观看影片大全网站| 国产主播在线观看一区二区| 国产成人精品在线电影| 午夜精品在线福利| www.精华液| 亚洲一区二区三区不卡视频| 18禁观看日本| 亚洲熟妇熟女久久| 一边摸一边抽搐一进一小说| 最新在线观看一区二区三区| 国产亚洲精品av在线| 午夜福利高清视频| 国产主播在线观看一区二区| 99国产综合亚洲精品| 后天国语完整版免费观看| 欧美中文综合在线视频| 中文字幕av电影在线播放| 一边摸一边抽搐一进一出视频| 色在线成人网| 亚洲九九香蕉| 香蕉国产在线看| 成人三级做爰电影| 一区二区三区国产精品乱码| 亚洲 欧美 日韩 在线 免费| 91成人精品电影| 免费高清视频大片| 亚洲精华国产精华精| av片东京热男人的天堂| 中文字幕久久专区| 亚洲av五月六月丁香网| 午夜视频精品福利| 久久精品aⅴ一区二区三区四区| 国产黄a三级三级三级人| 欧美日韩乱码在线| 十八禁人妻一区二区| 欧美一级毛片孕妇| 精品第一国产精品| 欧美乱码精品一区二区三区| 免费在线观看影片大全网站| 欧美老熟妇乱子伦牲交| 精品欧美国产一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 国产99久久九九免费精品| 欧美乱色亚洲激情| 国产精品98久久久久久宅男小说| 久久婷婷成人综合色麻豆| 99国产极品粉嫩在线观看| 久热这里只有精品99| 51午夜福利影视在线观看| 国产麻豆69| 国产亚洲精品av在线| 一区二区三区高清视频在线| 国产精品一区二区三区四区久久 | 一本大道久久a久久精品| 午夜两性在线视频| 亚洲视频免费观看视频| 亚洲国产欧美网| 在线国产一区二区在线| www.自偷自拍.com| 淫秽高清视频在线观看| 久久性视频一级片| 99国产精品一区二区蜜桃av| 日韩一卡2卡3卡4卡2021年| 成人欧美大片| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 国产亚洲精品久久久久久毛片| 在线永久观看黄色视频| 一个人观看的视频www高清免费观看 | 免费一级毛片在线播放高清视频 | 国产熟女午夜一区二区三区| 国产av一区在线观看免费| 国产精品99久久99久久久不卡| 99热只有精品国产| 久99久视频精品免费| 在线国产一区二区在线| 午夜影院日韩av| 国产亚洲精品久久久久5区| 黄色片一级片一级黄色片| 好看av亚洲va欧美ⅴa在| 亚洲国产看品久久| 一级,二级,三级黄色视频| 两个人看的免费小视频| 亚洲专区中文字幕在线| 长腿黑丝高跟| 国产麻豆成人av免费视频| 午夜免费鲁丝| 又紧又爽又黄一区二区| 啦啦啦免费观看视频1| 亚洲av第一区精品v没综合| bbb黄色大片| 人妻丰满熟妇av一区二区三区| 午夜免费观看网址| 亚洲九九香蕉| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美精品综合一区二区三区| 岛国视频午夜一区免费看| 人人妻人人澡人人看| 少妇熟女aⅴ在线视频| 日韩 欧美 亚洲 中文字幕| 老鸭窝网址在线观看| 国产在线精品亚洲第一网站| 久热爱精品视频在线9| 午夜免费成人在线视频| 国产精品自产拍在线观看55亚洲| 美女 人体艺术 gogo| 一区二区三区国产精品乱码| 亚洲全国av大片| 亚洲精品美女久久久久99蜜臀| 十八禁人妻一区二区| 亚洲av电影不卡..在线观看| 午夜福利影视在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 99国产精品免费福利视频| 在线观看免费日韩欧美大片| 精品第一国产精品| 天天躁夜夜躁狠狠躁躁| 亚洲第一青青草原| 一进一出抽搐动态| 女人爽到高潮嗷嗷叫在线视频| 日本五十路高清| 国产高清有码在线观看视频 | 三级毛片av免费| 中文字幕高清在线视频| av在线播放免费不卡| 91国产中文字幕| 亚洲国产日韩欧美精品在线观看 | 天堂动漫精品| 国产精品美女特级片免费视频播放器 | 久久久国产成人免费| 一级毛片高清免费大全| 他把我摸到了高潮在线观看| 国产熟女午夜一区二区三区| 欧美日本视频| 精品人妻在线不人妻| 99久久综合精品五月天人人| 嫩草影视91久久| 老鸭窝网址在线观看| 在线观看舔阴道视频| 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 韩国精品一区二区三区| 午夜免费激情av| 亚洲va日本ⅴa欧美va伊人久久| 成人亚洲精品一区在线观看| 人人澡人人妻人| 国产不卡一卡二| 欧美人与性动交α欧美精品济南到| 欧美日韩黄片免| 欧美亚洲日本最大视频资源| 久热这里只有精品99| 丝袜人妻中文字幕| 亚洲av片天天在线观看| 男人舔女人下体高潮全视频| 久久久国产成人免费| 色老头精品视频在线观看| 久久久久久久午夜电影| 91麻豆av在线| 精品电影一区二区在线| 中亚洲国语对白在线视频| 国产熟女xx| 宅男免费午夜| 午夜老司机福利片| 亚洲欧美精品综合久久99| 999精品在线视频| 久久天躁狠狠躁夜夜2o2o| 色综合婷婷激情| www.999成人在线观看| 法律面前人人平等表现在哪些方面| 亚洲少妇的诱惑av| 99精品久久久久人妻精品| 久久国产精品影院| 午夜亚洲福利在线播放| 成人精品一区二区免费| 亚洲成av片中文字幕在线观看| 久久伊人香网站| 午夜福利18| 精品无人区乱码1区二区| 黄色 视频免费看| 这个男人来自地球电影免费观看| 搞女人的毛片| 亚洲人成77777在线视频| 国产99久久九九免费精品| 久久天堂一区二区三区四区| 色av中文字幕| 久久国产乱子伦精品免费另类| 热99re8久久精品国产| 久久久久久免费高清国产稀缺| 亚洲黑人精品在线| 日韩国内少妇激情av| 俄罗斯特黄特色一大片| 国产1区2区3区精品| 88av欧美| 最新美女视频免费是黄的| 国产黄a三级三级三级人| 国产亚洲精品一区二区www| 亚洲黑人精品在线| 黄色 视频免费看| 操美女的视频在线观看| 国产av一区二区精品久久| 精品久久久久久久久久免费视频| 嫩草影院精品99| 手机成人av网站| 国产亚洲av嫩草精品影院| 黑人巨大精品欧美一区二区蜜桃| 校园春色视频在线观看| 天堂动漫精品| 久久精品亚洲精品国产色婷小说| 搡老妇女老女人老熟妇| 少妇 在线观看| 久99久视频精品免费| 亚洲国产精品sss在线观看| 久久伊人香网站| 国产一区二区三区综合在线观看| 午夜福利影视在线免费观看| 这个男人来自地球电影免费观看| 自线自在国产av| 亚洲av电影不卡..在线观看| 中文字幕人妻丝袜一区二区| 一本综合久久免费| 亚洲人成电影观看| 国产精品影院久久| 国产精品永久免费网站| 成人免费观看视频高清| 亚洲国产精品合色在线| 一区在线观看完整版| 国产亚洲欧美在线一区二区| 中文字幕精品免费在线观看视频| 丝袜人妻中文字幕| 乱人伦中国视频| 多毛熟女@视频| 久久精品影院6| 超碰成人久久| 午夜精品在线福利| 亚洲国产毛片av蜜桃av| 国产aⅴ精品一区二区三区波| 久久久久久国产a免费观看| 欧美 亚洲 国产 日韩一| 在线观看免费日韩欧美大片| 国产亚洲精品第一综合不卡| 久久午夜综合久久蜜桃| 极品教师在线免费播放| 国产亚洲精品久久久久久毛片| 精品无人区乱码1区二区| 两个人免费观看高清视频| 日日摸夜夜添夜夜添小说| 99国产精品一区二区蜜桃av| 99久久久亚洲精品蜜臀av| 久久草成人影院| 99精品欧美一区二区三区四区| 精品久久久精品久久久| 桃红色精品国产亚洲av| 久久久国产欧美日韩av| 人人澡人人妻人| 午夜亚洲福利在线播放| 国产一区二区三区视频了| 两人在一起打扑克的视频| 国产一卡二卡三卡精品| 精品卡一卡二卡四卡免费| 久久精品91无色码中文字幕| 国产精品久久久久久亚洲av鲁大| 午夜久久久在线观看| 国产一区二区激情短视频| 9热在线视频观看99| 国产极品粉嫩免费观看在线| 精品一品国产午夜福利视频| 波多野结衣巨乳人妻| 一级片免费观看大全| 欧美亚洲日本最大视频资源| 成熟少妇高潮喷水视频| 欧美激情极品国产一区二区三区| 免费无遮挡裸体视频| 人人妻人人澡人人看| 国产av一区二区精品久久| 午夜日韩欧美国产| 十八禁人妻一区二区| 久久人人精品亚洲av| 久久久久久久久免费视频了| 免费观看人在逋| 久久婷婷成人综合色麻豆| 成人三级做爰电影| 久久久久久国产a免费观看| 亚洲精品一卡2卡三卡4卡5卡| x7x7x7水蜜桃| 黄片小视频在线播放| 国产亚洲精品一区二区www| 成人手机av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲在线自拍视频| 一夜夜www| 一进一出抽搐gif免费好疼| 午夜成年电影在线免费观看| 成人亚洲精品一区在线观看| 国产视频一区二区在线看| av片东京热男人的天堂| 欧美黑人欧美精品刺激| 久久久久国产精品人妻aⅴ院| 精品欧美国产一区二区三| 搡老岳熟女国产| 黄色毛片三级朝国网站| 久久精品91无色码中文字幕| 最好的美女福利视频网| 久久久久久久久免费视频了| а√天堂www在线а√下载| 纯流量卡能插随身wifi吗| 欧美中文综合在线视频| 天天添夜夜摸| 不卡一级毛片| 国产激情久久老熟女| 长腿黑丝高跟| 一本综合久久免费| 色精品久久人妻99蜜桃| 国产成人欧美| 欧美黄色淫秽网站| 国产精品久久久人人做人人爽| 精品一区二区三区av网在线观看| 久久久久久久午夜电影| 亚洲最大成人中文| 桃红色精品国产亚洲av| 日韩欧美国产一区二区入口| 午夜福利成人在线免费观看| 岛国视频午夜一区免费看| 在线国产一区二区在线| 69av精品久久久久久| 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 日韩欧美一区视频在线观看| 免费在线观看日本一区| 国产亚洲精品久久久久5区|