• <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)化分析
    可溶性大豆多糖超聲波提取工藝及其抗氧化性研究
    午夜激情久久久久久久| 亚洲精品aⅴ在线观看| 人妻一区二区av| av在线播放精品| 精品久久久精品久久久| 成人性生交大片免费视频hd| 亚洲精品影视一区二区三区av| 国产精品熟女久久久久浪| 水蜜桃什么品种好| 中文字幕亚洲精品专区| 青春草视频在线免费观看| 高清视频免费观看一区二区 | 亚洲自偷自拍三级| 久久久久久国产a免费观看| 亚洲国产欧美在线一区| 日韩av免费高清视频| 久久99热这里只有精品18| 九色成人免费人妻av| 午夜福利视频精品| 简卡轻食公司| 国产精品伦人一区二区| 午夜久久久久精精品| 日韩视频在线欧美| 一个人看视频在线观看www免费| 人人妻人人澡人人爽人人夜夜 | 人妻一区二区av| 欧美区成人在线视频| 日韩一本色道免费dvd| 日产精品乱码卡一卡2卡三| 男人舔奶头视频| 男人舔奶头视频| 国产一区二区在线观看日韩| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久久久黄片| 国产亚洲91精品色在线| 丰满乱子伦码专区| 搡老妇女老女人老熟妇| 性插视频无遮挡在线免费观看| 大香蕉97超碰在线| 成年人午夜在线观看视频 | 91精品一卡2卡3卡4卡| 成年版毛片免费区| 亚洲欧美精品自产自拍| 神马国产精品三级电影在线观看| 亚洲精品久久久久久婷婷小说| 亚洲精品久久久久久婷婷小说| 精品人妻熟女av久视频| 久久国产乱子免费精品| 麻豆国产97在线/欧美| 18禁在线播放成人免费| 婷婷色综合www| 久久人人爽人人片av| 午夜精品国产一区二区电影 | 久久久久网色| 日韩欧美精品免费久久| 日本黄色片子视频| 一级av片app| 午夜久久久久精精品| 日韩欧美一区视频在线观看 | 久久精品国产自在天天线| 26uuu在线亚洲综合色| 亚洲熟妇中文字幕五十中出| 我要看日韩黄色一级片| 日本av手机在线免费观看| 国产免费视频播放在线视频 | 亚洲精品视频女| 中文字幕av成人在线电影| 国产精品一及| 搡老乐熟女国产| 一区二区三区四区激情视频| 卡戴珊不雅视频在线播放| 亚洲精品久久午夜乱码| 亚洲av成人av| 少妇被粗大猛烈的视频| 欧美xxⅹ黑人| 精品少妇黑人巨大在线播放| 国产一区有黄有色的免费视频 | 精品久久久久久久人妻蜜臀av| videos熟女内射| 搞女人的毛片| 欧美不卡视频在线免费观看| 日韩人妻高清精品专区| 男女边摸边吃奶| 少妇丰满av| 亚洲国产欧美在线一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久久人人人人人人| 边亲边吃奶的免费视频| 久久精品国产亚洲网站| 汤姆久久久久久久影院中文字幕 | 99久久九九国产精品国产免费| 亚洲av日韩在线播放| 日本色播在线视频| 午夜激情久久久久久久| 偷拍熟女少妇极品色| 好男人视频免费观看在线| 国产片特级美女逼逼视频| 欧美丝袜亚洲另类| 丝袜喷水一区| 国产成人精品福利久久| 亚洲精品自拍成人| 一级毛片 在线播放| 久久精品综合一区二区三区| 嫩草影院入口| 春色校园在线视频观看| 国产精品美女特级片免费视频播放器| 精品久久久久久成人av| 老女人水多毛片| 天堂av国产一区二区熟女人妻| 亚洲性久久影院| 国产午夜福利久久久久久| 深夜a级毛片| 亚洲国产欧美在线一区| 久久97久久精品| 最近的中文字幕免费完整| 中文在线观看免费www的网站| 亚洲精华国产精华液的使用体验| 三级男女做爰猛烈吃奶摸视频| 日日啪夜夜撸| 日韩制服骚丝袜av| 18禁在线播放成人免费| 午夜精品在线福利| av免费在线看不卡| 两个人视频免费观看高清| 能在线免费观看的黄片| 国产精品麻豆人妻色哟哟久久 | 欧美一区二区亚洲| 亚洲最大成人手机在线| 精品少妇黑人巨大在线播放| 91狼人影院| 国产爱豆传媒在线观看| 亚洲18禁久久av| 午夜免费激情av| 乱码一卡2卡4卡精品| 夜夜看夜夜爽夜夜摸| 免费不卡的大黄色大毛片视频在线观看 | 街头女战士在线观看网站| 久久99热这里只有精品18| 日本av手机在线免费观看| 亚洲精品456在线播放app| 97超碰精品成人国产| 国产熟女欧美一区二区| 又大又黄又爽视频免费| 嫩草影院新地址| 国产69精品久久久久777片| 国产精品人妻久久久久久| 一区二区三区高清视频在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产午夜福利久久久久久| 国产成人免费观看mmmm| 一区二区三区免费毛片| 日本黄色片子视频| 真实男女啪啪啪动态图| 网址你懂的国产日韩在线| 国产高潮美女av| 18禁动态无遮挡网站| 久久久a久久爽久久v久久| 国产老妇女一区| 国产精品一区二区三区四区久久| 日韩av在线大香蕉| 人妻少妇偷人精品九色| 亚洲av成人精品一区久久| 国产熟女欧美一区二区| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频 | 男人和女人高潮做爰伦理| 中文字幕久久专区| 熟女电影av网| 亚洲综合色惰| 毛片女人毛片| 免费少妇av软件| 日本熟妇午夜| 人人妻人人澡欧美一区二区| 午夜激情久久久久久久| 国产黄色小视频在线观看| av国产免费在线观看| 我的女老师完整版在线观看| 久久综合国产亚洲精品| 精品一区二区三区视频在线| 亚洲美女搞黄在线观看| 亚洲人成网站在线观看播放| 人体艺术视频欧美日本| 亚洲自拍偷在线| 久久久精品94久久精品| 国产片特级美女逼逼视频| 亚洲av成人精品一区久久| 美女脱内裤让男人舔精品视频| 日韩精品有码人妻一区| 欧美成人a在线观看| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 少妇丰满av| 日日摸夜夜添夜夜爱| 久久久久久久大尺度免费视频| 大香蕉久久网| eeuss影院久久| 国精品久久久久久国模美| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 国产av码专区亚洲av| 久久99热这里只有精品18| 午夜福利视频精品| 高清毛片免费看| 免费看光身美女| 黄片无遮挡物在线观看| 中国国产av一级| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 天天躁夜夜躁狠狠久久av| 国产亚洲91精品色在线| 国产在视频线精品| 国产 一区 欧美 日韩| 天堂中文最新版在线下载 | ponron亚洲| 日韩成人av中文字幕在线观看| 国产精品久久久久久久久免| 日韩人妻高清精品专区| 一本久久精品| 欧美激情国产日韩精品一区| 国产成年人精品一区二区| 亚洲精品456在线播放app| 国内精品美女久久久久久| 免费观看性生交大片5| 中文字幕久久专区| 久久久久精品久久久久真实原创| 久久精品久久久久久久性| 精品熟女少妇av免费看| 美女内射精品一级片tv| 日日摸夜夜添夜夜添av毛片| 九色成人免费人妻av| 亚洲美女搞黄在线观看| 国产极品天堂在线| 精品一区二区三区视频在线| 少妇熟女aⅴ在线视频| 日韩三级伦理在线观看| 高清毛片免费看| 一区二区三区高清视频在线| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂 | 男人爽女人下面视频在线观看| 免费大片黄手机在线观看| 亚洲内射少妇av| 蜜桃久久精品国产亚洲av| 亚洲精品亚洲一区二区| 精品久久久精品久久久| 男人狂女人下面高潮的视频| 两个人的视频大全免费| 国产亚洲av嫩草精品影院| 成人高潮视频无遮挡免费网站| 国产精品人妻久久久影院| 欧美丝袜亚洲另类| 国产午夜福利久久久久久| av在线播放精品| 天堂av国产一区二区熟女人妻| 亚洲综合色惰| 久久久久久久大尺度免费视频| 日韩人妻高清精品专区| 五月天丁香电影| 国产白丝娇喘喷水9色精品| 亚洲最大成人av| 久久久久久久亚洲中文字幕| 美女国产视频在线观看| 最近2019中文字幕mv第一页| 成年女人看的毛片在线观看| 毛片一级片免费看久久久久| 中文资源天堂在线| 色综合站精品国产| av专区在线播放| 九九久久精品国产亚洲av麻豆| 久久精品久久久久久久性| 纵有疾风起免费观看全集完整版 | 91精品一卡2卡3卡4卡| av又黄又爽大尺度在线免费看| 又爽又黄无遮挡网站| 日本免费在线观看一区| 大片免费播放器 马上看| 麻豆成人av视频| 久久久国产一区二区| 又爽又黄无遮挡网站| 亚洲欧美一区二区三区黑人 | 看黄色毛片网站| 亚洲精品色激情综合| 99热这里只有精品一区| 97人妻精品一区二区三区麻豆| 97精品久久久久久久久久精品| 美女脱内裤让男人舔精品视频| 久久97久久精品| 国产精品不卡视频一区二区| 国产一区亚洲一区在线观看| 午夜激情欧美在线| 人人妻人人澡人人爽人人夜夜 | 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 国内精品宾馆在线| 在线天堂最新版资源| 日本爱情动作片www.在线观看| 日韩欧美一区视频在线观看 | 国产欧美日韩精品一区二区| 亚洲精品色激情综合| 亚洲精品国产av蜜桃| 国产精品福利在线免费观看| 一本一本综合久久| 伊人久久精品亚洲午夜| 久久久久久久久久黄片| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 精品酒店卫生间| 国产探花极品一区二区| 国产亚洲5aaaaa淫片| 视频中文字幕在线观看| 18禁动态无遮挡网站| 大陆偷拍与自拍| 天堂√8在线中文| 秋霞在线观看毛片| 在线播放无遮挡| 哪个播放器可以免费观看大片| 毛片一级片免费看久久久久| 啦啦啦韩国在线观看视频| 纵有疾风起免费观看全集完整版 | 国产精品久久视频播放| 秋霞伦理黄片| 亚洲精品成人久久久久久| 最近手机中文字幕大全| 欧美变态另类bdsm刘玥| 欧美日韩精品成人综合77777| videossex国产| 国产伦精品一区二区三区四那| 深夜a级毛片| 亚洲欧美清纯卡通| 少妇的逼水好多| 最新中文字幕久久久久| 亚洲精品国产av成人精品| 18+在线观看网站| 亚洲欧美成人综合另类久久久| 国内精品宾馆在线| 欧美高清成人免费视频www| 久久精品夜夜夜夜夜久久蜜豆| 国产片特级美女逼逼视频| 国产成人freesex在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美性感艳星| 内地一区二区视频在线| 91久久精品电影网| 两个人视频免费观看高清| 韩国av在线不卡| 日韩一区二区视频免费看| 日韩一本色道免费dvd| 五月玫瑰六月丁香| 国产黄色视频一区二区在线观看| 国产成人精品福利久久| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 国产一级毛片七仙女欲春2| 国产精品爽爽va在线观看网站| 午夜精品国产一区二区电影 | 国产又色又爽无遮挡免| 国产精品一及| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| 国产精品精品国产色婷婷| 欧美高清成人免费视频www| 在线播放无遮挡| 国产av不卡久久| 熟妇人妻不卡中文字幕| 国产三级在线视频| 男人和女人高潮做爰伦理| 国产精品.久久久| 亚洲欧洲日产国产| 国产精品一二三区在线看| 全区人妻精品视频| 国精品久久久久久国模美| 国产激情偷乱视频一区二区| 22中文网久久字幕| 亚洲国产精品成人综合色| 免费看不卡的av| 草草在线视频免费看| 免费观看性生交大片5| 老女人水多毛片| 久久久午夜欧美精品| 好男人视频免费观看在线| 国产亚洲精品av在线| www.色视频.com| 简卡轻食公司| 免费大片黄手机在线观看| 夫妻午夜视频| 在线播放无遮挡| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人久久爱视频| 亚洲无线观看免费| av免费在线看不卡| 婷婷色av中文字幕| 在线a可以看的网站| 噜噜噜噜噜久久久久久91| 中文字幕制服av| 亚洲av中文av极速乱| 午夜福利在线观看吧| 高清视频免费观看一区二区 | 免费看日本二区| 我的老师免费观看完整版| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 97人妻精品一区二区三区麻豆| 日韩亚洲欧美综合| 日本熟妇午夜| 精品久久久久久久久av| 国产久久久一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲精品影视一区二区三区av| 小蜜桃在线观看免费完整版高清| 精品一区二区免费观看| 美女高潮的动态| 我的老师免费观看完整版| av卡一久久| 亚洲av中文字字幕乱码综合| 国产成人免费观看mmmm| 欧美极品一区二区三区四区| 日韩精品青青久久久久久| 久久久色成人| 国产精品1区2区在线观看.| 午夜亚洲福利在线播放| 国产午夜精品一二区理论片| 在线免费观看的www视频| 国产v大片淫在线免费观看| 男女啪啪激烈高潮av片| 99久久精品国产国产毛片| 人妻一区二区av| 国产亚洲91精品色在线| 亚洲久久久久久中文字幕| 22中文网久久字幕| 亚洲国产最新在线播放| 欧美人与善性xxx| 亚洲国产欧美人成| 久久久a久久爽久久v久久| 日韩中字成人| 建设人人有责人人尽责人人享有的 | 亚洲av国产av综合av卡| 国语对白做爰xxxⅹ性视频网站| 九草在线视频观看| 18禁在线播放成人免费| 亚洲av中文字字幕乱码综合| 亚洲精品色激情综合| av线在线观看网站| 精华霜和精华液先用哪个| 婷婷色av中文字幕| 国产精品日韩av在线免费观看| 波多野结衣巨乳人妻| 视频中文字幕在线观看| 国产在线男女| 有码 亚洲区| 美女被艹到高潮喷水动态| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 国产 一区精品| 免费观看性生交大片5| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 少妇人妻精品综合一区二区| 欧美人与善性xxx| 久久综合国产亚洲精品| 国产v大片淫在线免费观看| 亚洲熟女精品中文字幕| 国产视频内射| 国产av在哪里看| a级毛色黄片| freevideosex欧美| 成人无遮挡网站| 久久久精品免费免费高清| 看非洲黑人一级黄片| 色视频www国产| 国产三级在线视频| 色视频www国产| 国产毛片a区久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 蜜臀久久99精品久久宅男| 免费观看的影片在线观看| 国产精品一区二区在线观看99 | 久久久久久久久久久丰满| 黄片wwwwww| 极品教师在线视频| 精品人妻视频免费看| 人妻系列 视频| 亚洲精品久久久久久婷婷小说| 国产黄色视频一区二区在线观看| 永久免费av网站大全| 中文乱码字字幕精品一区二区三区 | 久久久成人免费电影| 狠狠精品人妻久久久久久综合| 国产真实伦视频高清在线观看| 欧美丝袜亚洲另类| 国产极品天堂在线| 少妇熟女欧美另类| 神马国产精品三级电影在线观看| 最近的中文字幕免费完整| 男女边摸边吃奶| 午夜激情欧美在线| 亚洲精品乱码久久久久久按摩| 乱系列少妇在线播放| 91久久精品电影网| 能在线免费观看的黄片| 麻豆久久精品国产亚洲av| 亚洲最大成人手机在线| 夫妻性生交免费视频一级片| videossex国产| 欧美三级亚洲精品| 成人性生交大片免费视频hd| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 日本熟妇午夜| 日日摸夜夜添夜夜爱| 国内精品一区二区在线观看| 日本免费a在线| 亚洲无线观看免费| 成人漫画全彩无遮挡| 人妻一区二区av| 日本免费a在线| 国产一区二区三区综合在线观看 | 成年女人在线观看亚洲视频 | 在线观看美女被高潮喷水网站| 久久亚洲国产成人精品v| 偷拍熟女少妇极品色| 欧美激情国产日韩精品一区| 嫩草影院新地址| 国产精品一区二区三区四区免费观看| 久久97久久精品| 亚洲经典国产精华液单| 久久午夜福利片| 全区人妻精品视频| 男女边摸边吃奶| 中文字幕免费在线视频6| 国产精品一区二区三区四区久久| 十八禁网站网址无遮挡 | 亚洲熟女精品中文字幕| 日韩 亚洲 欧美在线| 日韩人妻高清精品专区| 人妻少妇偷人精品九色| ponron亚洲| 97精品久久久久久久久久精品| 男的添女的下面高潮视频| 亚洲人与动物交配视频| 久久热精品热| 国产乱来视频区| 亚洲成人中文字幕在线播放| 精品久久久久久成人av| 精品一区二区免费观看| 丝袜喷水一区| 国产亚洲91精品色在线| 日韩欧美一区视频在线观看 | 久久人人爽人人片av| 国产三级在线视频| 久热久热在线精品观看| 免费播放大片免费观看视频在线观看| 亚洲精品国产av蜜桃| 亚洲国产精品专区欧美| 99热这里只有是精品50| 午夜激情久久久久久久| 少妇的逼水好多| 日本免费a在线| 白带黄色成豆腐渣| 春色校园在线视频观看| 国产精品一区二区性色av| 久久久久精品性色| 亚洲精品自拍成人| 国产视频内射| 日韩不卡一区二区三区视频在线| 亚洲精品一区蜜桃| 日韩欧美三级三区| 夜夜爽夜夜爽视频| 淫秽高清视频在线观看| 99re6热这里在线精品视频| 最近中文字幕高清免费大全6| 网址你懂的国产日韩在线| videossex国产| 三级国产精品欧美在线观看| 插阴视频在线观看视频| 亚洲欧美精品自产自拍| 亚洲成人久久爱视频| 亚洲欧美精品专区久久| 久久久久久久大尺度免费视频| 亚洲av男天堂| 国产色婷婷99| 欧美精品国产亚洲| 一级毛片aaaaaa免费看小| 免费电影在线观看免费观看| 亚洲三级黄色毛片| av在线天堂中文字幕| 看十八女毛片水多多多| 2021少妇久久久久久久久久久| 91午夜精品亚洲一区二区三区| 男人狂女人下面高潮的视频| 99久久精品热视频| 青春草视频在线免费观看| 最近视频中文字幕2019在线8| 国产精品99久久久久久久久| 久久韩国三级中文字幕| 精品国产露脸久久av麻豆 | 国产午夜精品一二区理论片| 69人妻影院| 亚洲色图av天堂| 一本久久精品| 中文天堂在线官网| 干丝袜人妻中文字幕| 又大又黄又爽视频免费| 老师上课跳d突然被开到最大视频| 777米奇影视久久| 少妇的逼水好多| 女人十人毛片免费观看3o分钟| 全区人妻精品视频| 少妇熟女欧美另类| 在线观看一区二区三区| 亚洲在线观看片| 欧美性猛交╳xxx乱大交人| 欧美不卡视频在线免费观看| 51国产日韩欧美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 岛国毛片在线播放|