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

    低周反復(fù)試驗(yàn)軸力加載裝置的系統(tǒng)摩擦力

    2019-06-10 00:51楊紅郭超王志軍羅航
    土木建筑與環(huán)境工程 2019年2期
    關(guān)鍵詞:摩擦力

    楊紅 郭超 王志軍 羅航

    摘 要:在柱、剪力墻等的低周反復(fù)加載試驗(yàn)中,與施加軸壓力的豎向作動器相連的鋼鉸軸處或滑動導(dǎo)軌處存在的摩擦力會導(dǎo)致加載水平力的測量結(jié)果不同程度地較真實(shí)值偏高。通過試驗(yàn)和有限元計算,分析了該摩擦力的規(guī)律和量化取值。為了消除材料非線性并準(zhǔn)確考慮幾何非線性的影響,設(shè)計兩個型鋼柱進(jìn)行了材料彈性狀態(tài)下的低周反復(fù)試驗(yàn),并在OpenSees平臺上建立了能夠準(zhǔn)確計算二階效應(yīng)的有限元模型。根據(jù)試驗(yàn)測量結(jié)果、有限元分析結(jié)果對該試驗(yàn)的軸力加載裝置引起的摩擦力進(jìn)行計算,研究了摩擦系數(shù)的取值方法,并考察了摩擦力對柱試驗(yàn)結(jié)果的影響。研究結(jié)果表明:由軸力加載裝置引起的摩擦力主要受軸壓力值影響,柱頂位移的影響較小;隨著軸壓力增大,摩擦力增大、摩擦系數(shù)變小;最大摩擦力、最小摩擦系數(shù)均趨于較穩(wěn)定的取值;對于軸壓力為350~900 kN的鋼筋混凝土柱,由該裝置摩擦力引起的低周反復(fù)試驗(yàn)的柱頂水平力測量結(jié)果誤差約為5~13 kN。

    關(guān)鍵詞:低周反復(fù)試驗(yàn);豎向作動器;摩擦力;有限元模擬;二階效應(yīng)

    中圖分類號:TU375.3

    文獻(xiàn)標(biāo)志碼:A? 文章編號:2096-6717(2019)02-0122-11

    Abstract:In the cyclic loading test of column or shear wall with the axial compression load, frictional force exists in the roller of steel hinges or the sliding rail, which is connected to the vertical actuator for applying axial force. The measurement of lateral load would be higher than its actual value due to this friction. On the basis of experimental research and finite element simulation, the variation characteristics and quantization method of the friction were investigated. In order to eliminate the influence of material nonlinearity and calculate geometric nonlinearity accurately, two steel columns were designed and cyclic loading tests under elastic situation were carried out. Besides, a finite element model which can accurately calculate the second-order effects was established on the OpenSees platform. Based on experimental and numerical results, the frictional forces were calculated, the method of determining the friction coefficient was proposed, and the influence of friction on the experimental results was evaluated. It was revealed that the frictional force caused by the devices for applying axial force was mainly affected by the axial compression force, while the top displacement of columns was of weak influence. As the axial compression force increases, the frictional force increases and the friction coefficient decreases. However, the maximum frictional force or the minimum friction coefficient converges. In cyclic loading test of the reinforced concrete columns with axial force of 350~900 kN,the measurement error of column top horizontal force caused by friction is 5~13 kN approximately.

    Keywords:low cyclic loading test; vertical actuator; frictional force; finite element simulation; second-order effect

    采用截面纖維模型,基于力的非線性梁柱單元等細(xì)化數(shù)值模擬方法已經(jīng)能較好地對鋼筋混凝土柱在彈塑性階段的加卸載剛度、滯回特征等進(jìn)行有效預(yù)測。但陸新征等[1-3] 、戎賢等[4]、Moharrami 等[5]在研究中發(fā)現(xiàn),鋼筋混凝土(RC)柱的模擬承載力不同程度地低于試驗(yàn)結(jié)果,且隨著軸壓力增大,兩者的差異相應(yīng)加大。分析結(jié)果表明[6-7],這種系統(tǒng)性誤差主要與摩擦阻力有關(guān),在柱、剪力墻、梁柱組合體試件等的低周反復(fù)試驗(yàn)中,施加軸壓力的豎向作動器主要通過鋼鉸或滾軸導(dǎo)軌與反力架相連,試驗(yàn)過程中,鋼鉸和滾軸隨柱頂位移而滑動或滾動必然會引起相應(yīng)的摩阻作用,因此,該摩阻力將導(dǎo)致由水平作動器的力傳感器測量的水平力較柱承受的真實(shí)值偏大。

    《建筑抗震試驗(yàn)方法規(guī)程》[8] 并未充分考慮摩擦力影響,僅規(guī)定滑動導(dǎo)軌的摩擦系數(shù)不應(yīng)大于0.01。李林安等[9]通過機(jī)械試驗(yàn)分析了鋼鉸、軸承的滑動、滾動摩擦的性能,分別得出了與建筑抗震試驗(yàn)方法規(guī)程不同的鉸軸滑動摩擦系數(shù)解析式;《機(jī)械設(shè)計手冊》[10]的摩擦系數(shù)取值也與《建筑抗震試驗(yàn)方法規(guī)程》不同;Tanaka 等[11]研究了往復(fù)滑動中初始接觸面條件對鋼軸承摩擦系數(shù)的影響;Hanrahan等[12]研究了滾動摩擦系統(tǒng)摩擦系數(shù)取值;吳家杰等[13] 通過不同的豎向加載裝置分析了摩擦力對 RC 橋墩抗震性能影響。為了消除(減?。┠Σ亮σ鸬恼`差,Lavorato等[14]、Chang等[15]采用了不設(shè)置滑動支座的豎向力加載裝置進(jìn)行抗震性能試驗(yàn)研究。

    目前,學(xué)者們對低周反復(fù)試驗(yàn)中因軸力加載裝置引起的摩擦阻力如何計算缺乏研究。以型鋼柱為例,通過彈性階段的低周反復(fù)試驗(yàn)和相應(yīng)的有限元模擬,建立了一種確定軸力加載裝置引起摩擦力的計算方法;通過對RC柱的水平力測量值進(jìn)行摩擦修正,提高了試驗(yàn)結(jié)果的準(zhǔn)確性。

    1 摩擦力量化分析方法與試驗(yàn)設(shè)計

    1.1 摩擦力定量計算方法

    以圖1(a)所示懸臂柱為例,施加軸壓力的豎向作動器上端、下端均通過鋼鉸軸(圓柱形的鋼軸承穿過鋼鉸上擺、下擺的鉸孔,形成鉸接連接,以下簡稱鉸軸)與試驗(yàn)柱構(gòu)件、反力系統(tǒng)相連,隨著柱頂水平位移增大,鉸軸轉(zhuǎn)動將引起相應(yīng)的摩擦阻力。

    整個加載裝置主要在圖1(a)所示的4處位置產(chǎn)生摩擦阻力,分別是:豎向千斤頂與反力架連接的鉸軸(編號1);豎向千斤頂與構(gòu)件連接的鉸軸(編號2);柱頂與鋼鉸連接面之間可能存在滑動,形成相應(yīng)的滑動摩擦力(編號3);水平千斤頂與構(gòu)件連接的鉸軸(編號4)。筆者將上述由于整個加載裝置引起的摩擦阻力統(tǒng)稱為系統(tǒng)摩擦力(簡稱摩擦力)。

    在試驗(yàn)過程中,上述4處位置的摩擦力并不相同。即使對于編號1、編號2的單個鉸軸而言,其摩擦狀態(tài)也較為復(fù)雜。試驗(yàn)過程中,當(dāng)豎向作動器隨柱頂位移左右擺動時,鋼軸承與鋼鉸上、下擺產(chǎn)生相對轉(zhuǎn)動,同時鋼軸承也在上擺、下擺的鉸孔空隙處相對滑動(當(dāng)鋼軸承加工質(zhì)量偏低、鉸孔直徑偏大時,這種滑動更明顯),從而形成滾動摩擦和滑動摩擦混合作用的狀態(tài)。此外,鋼軸承和上下擺鉸孔的初始形狀并非理想圓形,隨著壓應(yīng)力加大,在軸承滾動時,軸承和鉸孔的接觸面形狀將發(fā)生微小變化,接觸面愈軟、形狀變化愈大,滾動摩擦力就愈大。綜上所述,整個加載裝置的系統(tǒng)摩擦力與鋼軸承、鋼鉸孔接觸面的光滑程度、制作工藝等有關(guān);此外,試驗(yàn)時啟動加載、采數(shù)停頓、峰值卸載等可能產(chǎn)生靜摩擦作用,其與水平位移行進(jìn)中的動摩擦交替作用,進(jìn)一步增加了加載裝置摩擦狀態(tài)的復(fù)雜性。

    在圖1(a)所示的鋼筋混凝土懸臂柱的低周反復(fù)試驗(yàn)中,直接測量的水平力位移曲線的水平力比真實(shí)值偏大,有兩方面的原因:一是豎向作動器隨柱頂水平移動逐漸傾斜,在水平方向形成的分力(見圖1(b));二是軸力加載裝置引起的摩擦力(見圖1(c))。

    采用鉸軸連接方式的試驗(yàn)中,隨柱頂水平位移增加,豎向作用器與豎直方向形成的夾角(θ)會逐漸加大,故豎向作用器的軸向力除存在豎向分量(即柱承受的軸壓力)外,還存在一個水平分量FN(見圖1(b)),F(xiàn)N將使實(shí)際測量的水平力FP比柱實(shí)際承受的水平力更小。這種試驗(yàn)誤差可以在試驗(yàn)完成后直接根據(jù)實(shí)測數(shù)據(jù)進(jìn)行修正,其修正方法是:首先計算各柱頂位移下的θ (根據(jù)柱頂水平位移、豎向作用器上下鉸軸的中心距實(shí)測結(jié)果進(jìn)行計算),然后計算軸壓力N的分力FN,將其與水平作動器的實(shí)測結(jié)果FP疊加,結(jié)果即為考慮摩擦力修正之前的RC柱頂承受的水平力。

    對于圖1(a)所示的軸力加載裝置(豎向作動器通過鉸軸與試件、反力系統(tǒng)連接),通過3個步驟定量計算該裝置引起的摩擦力。

    第1步,對型鋼柱進(jìn)行彈性階段的低周反復(fù)試驗(yàn),獲得水平力位移滯回曲線的測量結(jié)果。該實(shí)測水平力位移曲線包含了摩擦力Ff、軸力水平分量FN的影響。理論上看,型鋼柱彈性階段的水平力位移關(guān)系本是一條斜直線(若二階效應(yīng)的影響可近似忽略),但由于摩擦力等影響,實(shí)測滯回曲線會形成滯回圈。

    第2步,對于材料處于彈性狀態(tài)的型鋼柱,可通過有限元軟件精確地計算各水平位移下的柱頂水平力(FM)。FM的計算結(jié)果不受摩擦力Ff、軸力水平分量干擾FN,為無摩擦的理想狀態(tài)的結(jié)果。

    第3步,將試驗(yàn)結(jié)果、有限元計算結(jié)果相結(jié)合,可按式(1)計算各位移狀態(tài)下的摩擦力Ff。

    式中:Ff為施加軸力的豎向作動器等引起的系統(tǒng)摩擦力,F(xiàn)P為試驗(yàn)加載水平力(由與水平作動器相連的力傳感器直接測得),F(xiàn)N為豎向作動器軸向力的水平分量,F(xiàn)M為采用有限元模型計算的水平力,N為豎向作動器施加的軸向力(由與豎向作動器相連的力傳感器測得),θ為豎向作動器與豎直方向的夾角(見圖1(b))。分析時,式(1)中各水平力均以圖1(c)所示“加載正向”為正。

    在第1步中,采用型鋼柱試件進(jìn)行低周反復(fù)試驗(yàn)的目的是消除材料非線性的影響,但不利之處是鋼柱的二階效應(yīng)一般更為顯著。因此,在第2步計算時重點(diǎn)考慮了二階效應(yīng),建立能完善考慮二階效應(yīng)的有限元分析模型,并完成相應(yīng)計算。

    1.2 試驗(yàn)概況

    2個鋼柱試件編號為GZ-1、GZ-2,分別對應(yīng)文獻(xiàn)[16]中的RC柱SC-5和SC-6。應(yīng)強(qiáng)調(diào)的是,4個試件采用了同一套加載系統(tǒng)(包括反力架、豎向作動器兩端的連接鉸軸、水平和豎向作動器等),如圖2所示。

    合理設(shè)計型鋼試件,使其受力性能滿足如下要求:型鋼試件在整個實(shí)驗(yàn)過程中,始終在材料彈性階段且其彈性階段的最大頂點(diǎn)位移與相應(yīng)RC柱破壞時的最大水平位移相近,此外,型鋼試件目標(biāo)位移下的水平承載力也接近相應(yīng)RC柱的水平承載力。

    試件所用鋼材強(qiáng)度均為Q345B,試件尺寸見表1。

    按照《鋼結(jié)構(gòu)設(shè)計規(guī)范》[17]、《鋼結(jié)構(gòu)原理與設(shè)計》[18],完成了試件腹板與翼緣連接焊縫、柱腳和柱頂連接鉸以及其連接焊縫的設(shè)計、計算,結(jié)果如圖3所示。測得的鋼材基本力學(xué)性能見表2。

    采用擬靜力方法對試件進(jìn)行加載,加載裝置如圖2示。試驗(yàn)中,豎向軸力由150 t拉壓作動器施加,水平作用力由50 t拉壓作動器施加。水平荷載由水平作動器端部的50 t拉壓傳感器測量,試件柱頂位移由位移計測量。試驗(yàn)過程中通過采集的應(yīng)變數(shù)據(jù)隨時判斷試件是否處于彈性狀態(tài)。應(yīng)變片共設(shè)4排,各排應(yīng)變片距離柱底座上表面分別為50、150、250、350 mm。

    試件GZ-1分別在軸力等于288、576、928 kN下各進(jìn)行了3次彈性低周反復(fù)加載(共9次試驗(yàn)),試件GZ-2分別在軸力等于122 、244、 365 、576 、928、1 293 kN下各進(jìn)行了4次彈性低周反復(fù)加載(共24次試驗(yàn)),每次試驗(yàn)結(jié)束后、開始下次試驗(yàn)前,均將鉸軸與試件拆開并重新連接,使摩擦力歸零。各試驗(yàn)均采用位移控制,其加載制度見圖4,圖中Δ為型鋼柱柱頂位移,其各循環(huán)的位移幅值分別等于文獻(xiàn)[16]中RC柱屈服位移Δcy(RC柱SC-5和SC-6屈服位移Δcy均近似取為7 mm)的1、2、3倍(即7、14、21 mm)。

    1.3 試驗(yàn)荷載水平位移滯回曲線及修正

    正式試驗(yàn)開始之前,對各試件進(jìn)行了無軸力的驗(yàn)證試驗(yàn)。圖5為試件GZ-1脫開柱頂豎向作動器、未施加軸力時的水平力(FP)-位移(Δ)曲線測量結(jié)果,可以發(fā)現(xiàn),在無軸力狀態(tài)下其試驗(yàn)結(jié)果仍然存滯回圈。

    由前文分析可知,柱頂施加軸力的鉸軸脫開后,不會出現(xiàn)摩擦力Ff、軸力水平分量FN的干擾。經(jīng)仔細(xì)分析發(fā)現(xiàn),實(shí)測FP-Δ曲線形成滯回圈的原因是,試驗(yàn)過程中型鋼柱試件頂部鋼板相對于鋼鉸底板略有滑動(摩擦型高強(qiáng)螺栓直徑較螺栓孔更?。?。為了消除該因素影響,布置百分表測量該相對位移,并據(jù)其修正柱頂位移。如圖6所示,百分表通過磁性基座固定在鋼柱頂部,Δs為鋼柱頂部相對于鋼鉸的相對位移(受拉為正),Δt為位移計測得的柱頂位移值(受拉為正),則真實(shí)的柱頂水平位移Δ=Δt-Δs。

    圖5(b)表明,經(jīng)修正后,無軸力下FP-Δ曲線不再出現(xiàn)滯回圈,表明其結(jié)果是合理的。

    2 摩擦力計算結(jié)果與變化規(guī)律

    2.1 有限元模擬與二階效應(yīng)

    為了得到考慮二階效應(yīng)后試件的荷載位移曲線計算結(jié)果,在OpenSees平臺上,選取nonlinearBeamColumn單元建立型鋼柱的有限元模型。由于nonlinearBeamColumn單元不能直接考慮P-δ效應(yīng),OpenSees也難以精確計算P-Δ效應(yīng),故采用舒啟軍[19]的研究成果進(jìn)行計算。該方法根據(jù)Souza[20]的柔度法理論,通過在最終的剛度方程中直接按線彈性方式加入相應(yīng)項(xiàng),不考慮扭轉(zhuǎn)變形與其余變形間的耦合作用,建立能夠完全考慮二階效應(yīng)的柔度矩陣。

    考慮P-Δ效應(yīng)和P-δ效應(yīng)的截面力表達(dá)式為[19]

    張勇[21]和舒啟軍[19]通過算例和試驗(yàn)驗(yàn)證了采用上述方法建立的有限元模型對計算二階效應(yīng)的準(zhǔn)確性。計算時,鋼柱截面采用纖維模型,將試件的左右翼緣各劃分為5個纖維,腹板劃分為4個纖維;鋼材本構(gòu)采用Steel02模型,影響過渡曲線形狀的參數(shù)分別取20,0.925和0.15,應(yīng)變硬化率0.01。

    為了進(jìn)一步驗(yàn)證上述有限元模型的正確性,以2個試件的無軸力試驗(yàn)結(jié)果為依據(jù)對其進(jìn)行了校核。

    無軸力作用下FP-Δ曲線試驗(yàn)結(jié)果(已對柱頂連接螺栓滑動進(jìn)行了修正)和計算結(jié)果如圖7所示,表明無軸力時兩試件試驗(yàn)結(jié)果與模擬結(jié)果一致。

    2.2 摩擦力變化規(guī)律與計算結(jié)果

    2.2.1 荷載水平位移關(guān)系試驗(yàn)和有限元計算結(jié)果

    如前所述,試件GZ-1、GZ-2分別在不同軸壓力取值作用下進(jìn)行了數(shù)次獨(dú)立的低周反復(fù)加載試驗(yàn)。圖8和圖9給出了多次試驗(yàn)中的一組典型試驗(yàn)結(jié)果(均已對柱頂連接螺栓滑動、FN進(jìn)行了修正),以及與相應(yīng)有限元計算結(jié)果的對比。

    在圖8、圖9中,有限元計算時試件處于無摩擦力的理想狀態(tài),且二階效應(yīng)影響很小,故型鋼柱彈性階段的FM-Δ曲線計算結(jié)果近似為一條斜直線。對于試驗(yàn)結(jié)果,雖然試件材料處于彈性狀態(tài),但受Ff的影響,且摩擦力在加載、卸載過程中大小和方向均有所改變(具體分析詳后文),故試件的水平力測量結(jié)果相應(yīng)變化,并使柱頂水平力水平位移形成了彈性滯回圈。

    2.2.2 加載階段的摩擦力 圖8、圖9表明,在加載過程中(除加載初始階段外),試驗(yàn)滯回曲線基本與無摩擦力的模擬曲線的斜率相同(平行),即Ff在加載的大部分過程中基本保持不變,故近似選取±5、±7、±10、±14、±17、±21 mm等共6個位移點(diǎn)(其中,±7、±14、±21 mm為各循環(huán)的位移幅值)計算,并將其作為分析加載階段的代表值。因此,根據(jù)不同軸力下多次試驗(yàn)的結(jié)果、有限元模擬結(jié)果,按照式(1)即可計算摩擦力。

    以圖8 (c)中GZ-1試件件在928 kN下第1次低周反復(fù)試驗(yàn)結(jié)果為例,說明正向加載階段的摩擦力Ff計算方法。

    根據(jù)Ff的計算方法(式(1))可知,在圖8(c)中,紅色實(shí)線即為Ff。由圖4的加載制度可知,在每個峰值位移共循環(huán)兩次,則在5 mm位移處需計算并得到共6個Ff取值,將其平均值(記為F1f)作為5 mm位移處的摩擦力;同理,按照該方法可計算其它每個代表位移點(diǎn)處的摩擦力(5、7 mm處6個Ff的平均值,10、14 mm處4個Ff的平均值,17、21 mm處2個Ff的平均值),其具體計算過程見表3。

    根據(jù)表3中的摩擦力平均值(F1f)可繪制圖10中第1次正向加載的F1f曲線(第2次、第3次正向加載的摩擦力計算過程類似),負(fù)向加載的摩擦力采用與正向相同的計算過程,不再贅述。

    限于篇幅,圖10、圖11僅給出了試件GZ-1和GZ-2在部分軸力作用下,多次加載過程中的F1f的計算結(jié)果。其中,圖10正向加載、負(fù)向加載摩擦力不同的原因是,試件GZ-1制作時略有平面外偏斜,使得各軸力下正向、負(fù)向加載段摩擦力不同,且軸力為288 kN時兩者的差值最顯著。

    如前所述,每次試驗(yàn)開始前均將鉸軸與試件拆開并重新連接。圖10、圖11的試驗(yàn)?zāi)Σ亮1f存在一定離散性,說明即使是相同的軸力加載裝置(鋼鉸上下擺、鋼軸承、作動器等均相同)和柱試件,重新安裝后,其摩擦系數(shù)將發(fā)生改變。分析認(rèn)為,其主要原因是,制作的上下擺鉸孔和鋼軸承不是理想的圓形、上下擺和軸承的表面并不完全光滑、平整,很多局部位置隨機(jī)地存在微小的凹凸,因此,重新連接鉸軸后,鋼軸承與鋼鉸孔的接觸位置改變,會導(dǎo)致摩擦系數(shù)變化??傊?,摩擦系數(shù)除與鋼軸承、鋼鉸孔(或?qū)к墸┙佑|面的光滑程度、制作工藝等有關(guān)外,也與接觸位置的隨機(jī)性有關(guān)。

    為了更直觀表現(xiàn)摩擦力與軸力的關(guān)系,將試件GZ-1和GZ-2在每一軸力作用下的各次試驗(yàn)加載過程中的 F1f(見圖10和圖11)再次進(jìn)行平均,得到F2f的計算結(jié)果,如圖12所示。

    圖12表明:試件GZ-1和GZ-2在不同軸力下的加載段摩擦力平均值(F2f)較穩(wěn)定,且隨軸壓力值加大摩擦力也相應(yīng)增大;柱頂位移的變化對摩擦力平均值影響不明顯,這與羅小龍等[22]、吳二軍等[23]關(guān)于摩擦力的研究結(jié)果相符。

    為了消除試件微小傾斜所導(dǎo)致的試驗(yàn)測量的正向水平力、負(fù)向水平力存在一定差異(試件GZ-1較明顯)的不利影響,不考慮方向性(即Ff的正負(fù))可更客觀地分析Ff的大小與柱頂位移的關(guān)系。因此,取正向、負(fù)向柱頂位移絕對值相同的加載段Ff的絕對值(見圖12)進(jìn)行平均,得到加載段摩擦力平均值(記為F3f)與柱頂位移的關(guān)系(見圖13)。由圖13可見,相同軸力下不同柱頂位移點(diǎn)的摩擦力(F3f)較接近,摩擦力隨軸壓力的加大而增大。

    2.2.3 卸載階段的摩擦力

    為分析卸載階段摩擦力的變化規(guī)律,分別計算試件GZ-1、GZ-2 在軸力為928 kN下的卸載階段摩擦力? (僅以柱頂位移幅值為±7、±14、±21 mm的3個循環(huán)的卸載段為例),結(jié)果如圖14所示。

    首先,以位移幅值為±21 mm的正向卸載階段為例說明Ff的變化過程:在21 mm位移幅值時Ff為負(fù)值,在21~15 mm左右Ff會迅速變化為0,在15~10 mm繼續(xù)迅速增加,然后,隨著柱頂位移減小,趨于穩(wěn)定(等于加載階段的Ff)。

    上述分析表明,在21 mm位移幅值處,試件處于加載、卸載的臨界(轉(zhuǎn)換)狀態(tài),加載階段末Ff為負(fù)值,但開始卸載后,隨著柱頂位移減小,F(xiàn)f會經(jīng)歷一個“負(fù)向→絕對值下降→零→絕對值增加→正向”的變化過程,即摩擦力會改變方向,從而導(dǎo)致圖14中藍(lán)色曲線在21~10 mm經(jīng)歷了由負(fù)值變?yōu)檎档倪^程。

    卸載階段Ff改變方向的根本原因是,F(xiàn)f的方向總是與運(yùn)動方向相反。在加載過程中,F(xiàn)f始終與運(yùn)動方向相反,也就總是與作動器水平力FP方向相反,故卸載前的大小、指向基本不變。開始卸載的初期,F(xiàn)P方向不變(與加載時相同),但很快柱頂位移運(yùn)動方向變化,故Ff改變方向(數(shù)學(xué)符號變化),并隨后逐漸增大。

    圖14中只有幅值±21 mm的卸載段有較為明顯的穩(wěn)定段,幅值±7 mm、±14 mm的卸載段基本上只存在方向改變的過程。出現(xiàn)這種現(xiàn)象的原因是,在圖8、圖9的試驗(yàn)滯回曲線與模擬的FM-Δ線相交(Ff改變方向)的過程中,±7 mm、±14 mm循環(huán)測得的FP在卸載段的變化幅度一直有較大的增大或者減少(表現(xiàn)為卸載剛度變化),模擬計算的FM基本線性變化,即其剛度基本不變,因此,按式(1)進(jìn)行計算會導(dǎo)致Ff的正→負(fù)(或負(fù)→正)的變化過程較長,故幅值±7 mm、±14 mm的卸載段Ff的穩(wěn)定段需柱頂水平位移進(jìn)一步加大才能出現(xiàn)。

    基于以上分析,結(jié)合圖8、圖9可知,卸載過程是形成彈性滯回圈的過程,因卸載過程中試驗(yàn)滯回曲線與有限元計算的FM-Δ線相交,故按式(1)計算的卸載摩擦力正、負(fù)號會變化,即摩擦力方向變化。

    與加載階段的分析方法類似,選取0、±3 mm、±5 mm、±7 mm、±10 mm、±14 mm、±17 mm、±21 mm等8個位移點(diǎn)作為代表,將不同軸力作用下的多次試驗(yàn)的卸載段計算結(jié)果進(jìn)行平均,并將Ff的平均值記為F2′f(即為兩次平均的結(jié)果,與加載階段的摩擦力平均值F2f相對應(yīng)),所得結(jié)果如圖15、圖16所示(限于篇幅,僅給出幅值±21 mm的結(jié)果)。

    圖15、圖16表明,各軸力作用下不同加載幅值點(diǎn)的卸載段摩擦力平均值與圖14規(guī)律類似,即卸載段摩擦力仍隨軸力的增大而增大。

    以上分析結(jié)果表明,即使卸載段的柱頂位移相同,在不同位移幅值循環(huán)下其摩擦力會明顯不同(見圖14),難以對各循環(huán)的卸載摩擦力平均值進(jìn)行準(zhǔn)確定量。

    由于卸載摩擦力的穩(wěn)定段即等于加載摩擦力,且卸載階段摩擦力對修正試驗(yàn)結(jié)果較為次要,對其計算方法暫不深入分析。

    2.2.4 摩擦系數(shù) 如前所述,可取正向、負(fù)向加載段相同柱頂位移點(diǎn)(即5、7、10、14、17、21 mm)的摩擦力平均值(F3f,見圖13)作為計算摩擦系數(shù)μf的依據(jù),以消除試件可能存在的微小傾斜的影響。在圖17中,μf的“平均值”連線是根據(jù)圖13的摩擦力、相應(yīng)的軸壓力相除得到的計算結(jié)果。為了直觀地體現(xiàn)μf的離散性,在圖17中還給出了以每次加載過程中(不同軸力下,試件GZ-2、GZ-1分別進(jìn)行了3次、4次試驗(yàn))的摩擦力平均值F1f為依據(jù),除以相應(yīng)軸壓力所得到的μf計算結(jié)果(仍取正向、負(fù)向μf的平均值),具體詳見圖17中的各散點(diǎn)。

    圖17表明,不同柱頂位移點(diǎn)的μf接近,柱頂位移對軸力加載裝置引起的摩擦力影響較小。

    為分析軸壓力N取值對摩擦系數(shù)μf的影響,以同一軸壓力下的多次試驗(yàn)結(jié)果經(jīng)計算所得的摩擦力平均值F3f(見圖13)作為計算μf的依據(jù),建立N與加載段μf平均值的關(guān)系,結(jié)果如圖18所示,圖中各散點(diǎn)為不同軸力下多次試驗(yàn)分別計算所得的μf,其計算方法與圖17相同。

    圖18表明,軸力加載裝置引起的摩擦力主要與軸壓力N有關(guān)。隨著N增大,F(xiàn)f增大、μf減小,且最大Ff、最小μf均趨于較穩(wěn)定取值;相同N作用下Ff在一定范圍內(nèi)變化,離散性不大;對于試件GZ-2、GZ-1采用的試驗(yàn)裝置,μf可取約0.75%~2.45%。

    3 摩擦力對RC柱試驗(yàn)結(jié)果的影響

    根據(jù)試件GZ-2、GZ-1的分析結(jié)果確定μf之后,由于采用的試驗(yàn)加載裝置相同,可對文獻(xiàn)[16]中2個RC柱試件SC-5(軸壓比0.1)和SC-6(軸壓比0.3)的P-Δ曲線試驗(yàn)結(jié)果(FN的修正已完成)進(jìn)行摩擦力修正。

    取與試件SC-5、SC-6軸壓力對應(yīng)的μf平均值對上述兩RC柱的水平力測量結(jié)果FP進(jìn)行修正。RC柱SC-5的目標(biāo)軸力為365 kN,按圖18可取μf為1.45%,對應(yīng)的Ff為5.3 kN;柱SC-6目標(biāo)軸力為928 kN,其μf可取為1.325%(見圖18),對應(yīng)的Ff為12.3 kN。

    考慮摩擦力修正后的試驗(yàn)結(jié)果對比如圖19所示。

    圖19表明,試件SC-5、SC-6的軸壓力不同,軸力加載裝置引起的Ff分別約為最大水平力實(shí)測結(jié)果的3.06%、5.27%。

    軸力加載裝置引起的摩擦力的大小、方向在卸載階段均發(fā)生改變,規(guī)律較為復(fù)雜。因采用簡化方法,圖19中摩擦力修正后的水平力在卸載段(峰值點(diǎn)附近為摩擦力過渡段)存在誤差。

    4 結(jié)論

    基于型鋼柱試件彈性階段的低周反復(fù)試驗(yàn),結(jié)合精細(xì)有限元分析模型,建立了一種量化軸力加載裝置所引起摩擦阻力的計算方法,得到以下主要結(jié)論:

    1)在低周反復(fù)試驗(yàn)中,軸力加載裝置將引起摩擦阻力,該摩擦力主要與軸壓力大小有關(guān),試件頂部水平加載位移的影響較小。

    2)隨著軸壓力增大,軸力加載裝置引起的摩擦力增大、相應(yīng)的摩擦系數(shù)減小,最大摩擦力、最小摩擦系數(shù)趨于較穩(wěn)定的取值。該試驗(yàn)的軸力加載裝置的摩擦系數(shù)可取0.75%~2.45%。

    3)軸力加載裝置的摩擦系數(shù)還與鋼軸承、鋼鉸上下擺鉸孔接觸面的光滑程度、制作工藝等有關(guān),也與接觸位置的隨機(jī)性有關(guān)。

    4)對于軸壓力為350~900 kN的RC柱,該裝置的摩擦力引起的水平力測量誤差約為5~13 kN。

    參考文獻(xiàn):

    [1] 陸新征, 葉列平, 潘鵬, 等. 鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究及數(shù)值模擬競賽Ⅰ:框架試驗(yàn) [J]. 建筑結(jié)構(gòu), 2012, 42(11): 19-22.

    LU X Z,YE L P,PAN P, et al. Pseudo-static collapse experiments and numerical prediction competition of RC frame structure I: RC frame experiment [J]. Building Structure,2012, 42(11): 19-22. (in Chinese)

    [2] 陸新征, 葉列平, 潘鵬, 等. 鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究及數(shù)值模擬競賽Ⅱ:關(guān)鍵構(gòu)件試驗(yàn)[J]. 建筑結(jié)構(gòu), 2012, 42(11):23-26.

    LU X Z,YE L P,PAN P, et al. Pseudo-static collapse experiments and numerical prediction competition of RC frame structure II: Key elements experiment [J]. Building Structure,2012, 42(11): 23-26. (in Chinese)

    [3] 中國建筑學(xué)會建筑結(jié)構(gòu)防倒塌專業(yè)委員會. 鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究及數(shù)值模擬競賽Ⅲ:模擬結(jié)果分析 [J]. 建筑結(jié)構(gòu), 2012, 42(11): 27-30.

    Collapse Prevention Committee of Architectural Society of China. Pseudo-static collapse experiments and numerical prediction competition of RC frame structure Ⅲ: Simulation results analysis [J]. Building Structure, 2012, 42(11): 27-30. (in Chinese)

    [4] 戎賢,喬超男. 600 MPa級鋼筋混凝土十字形柱抗震性能試驗(yàn)研究與恢復(fù)力特性分析 [J]. 土木建筑與環(huán)境工程, 2018, 40(3): 139-146.

    RONG X, QIAO C N. Experimental analysis on seismic behavior and analysis of restoring force characteristics of 600 MPa RC cross-shaped columns [J]. Journal of Civil, Architectural & Environmental Engineering, 2018, 40(3): 139-146. (in Chinese)

    [5] MOHARRAMI M, KOUTROMANOS I. Finite element analysis of damage and failure of reinforced concrete members under earthquake loading [J]. Earthquake Engineering & Structural Dynamics, 2017, 46(15): 2811-2829.

    [6] 楊紅, 張洛, 張和平. 考慮縱筋屈曲及疲勞損傷的鋼筋混凝土柱抗震性能試驗(yàn)研究與非線性分析 [J]. 建筑結(jié)構(gòu)學(xué)報, 2013, 34(11): 130-140.

    YANG H,ZHANG L,ZHANG H P. Experiments and nonlinear analysis on seismic behavior of RC columns considering buckling and fatigue damage of reinforcing steel bar [J]. Journal of Building Structures , 2013, 34(11): 130-140. (in Chinese)

    [7] YANG H, XIE Q, ZHAO W,et al. Seismic behavior comparison of reinforced concrete interior beamcolumn joints based on different loading methods [J]. Engineering Structures,2018,166(7):31-45.

    [8] 中華人民共和國建設(shè)部. 建筑鋼結(jié)構(gòu)焊接技術(shù)規(guī)程:JGJ 81—2002 [S]. 北京:中國建筑工業(yè)出版社, 2002.

    Ministry of Construction of the PRC. Technical specification for welding of steel structure of building: JGJ81-2002 [S]. Beijing:China Architecture & Building Press, 2002. (in Chinese)

    [9] 李林安, 岳澄, 佟景偉, 等. 鉸軸滑板鋼支座力學(xué)性能的實(shí)驗(yàn)分析 [J]. 實(shí)驗(yàn)力學(xué), 2003, 18(1): 29-33.

    LI L A, YUE C, TONG J W, et al. An experimental analysis for the mechanical properties of steel bearing of hinge-axis sliding mode [J]. Journal Experimental Mechanics, 2003, 18(1): 29-33. (in Chinese)

    [10] 聞邦椿, 張義民, 鄂中凱,等. 機(jī)械設(shè)計手冊 [M]. 北京:機(jī)械工業(yè)出版社, 2015.

    WEN B C, ZHANG Y M, E Z K, et al. Machinerys handbook [M]. Beijing:Mechanical Industry Press, 2015.(in Chinese)

    [11] TANAKA Y, ENDO M, MORIYAMA S. Effect of initial contact surface condition on the friction and wear properties of bearing steel in cyclic reciprocating sliding contact [C] // Journal of Physics: Conference Series, 6th International Conference on Fracture Fatigue and Wear, Porto, Portugal,2017.

    [12] HANRAHAN B M, MISRA S, BEYAZ M I, et al. An adhesion-dominated rolling friction regime unique to micro-scale ball bearings [J]. Tribology Letters, 2014, 56(2): 215-221.

    [13] 吳家杰, 薛俊青, Briseghella B, 等. 豎向加載裝置對RC橋墩抗震性能影響試驗(yàn)研究[C] //第25屆全國結(jié)構(gòu)工程學(xué)術(shù)會議, 福建 廈門,2016.

    WU J J, XUE J Q, BRISEGHELLA B, et al .Experimental research on seismic behaviors of RC pier under different vertical load system [C] // The 25th National Academic Conference on Structural Engineering, Xiamen, Fujian,China, 2016. (in Chinese)

    [14] LAVORATO D, NUTI C. Pseudo-dynamic tests on reinforced concrete bridges repaired and retrofitted after seismic damage [J]. Engineering Structures, 2015,94(7): 96-112.

    [15] CHANG S Y, LI Y F, LOH C H. Experimental study of seismic behaviors of as-built and carbon fiber reinforced plastics repaired reinforced concrete bridge columns [J]. Journal of Bridge Engineering, 2004, 9(4): 391-402.

    [16] 張和平. 鋼筋混凝土柱抗震性能試驗(yàn)及優(yōu)化模擬分析 [D]. 重慶:重慶大學(xué), 2012.

    ZHANG H P. Experiment and optimized modeling on seismic performance of RC column [D]. Chongqing:Chongqing University, 2012. (in Chinese)

    [17] 中華人民共和國建設(shè)部. 鋼結(jié)構(gòu)設(shè)計規(guī)范 [M]. 北京:中國計劃出版社, 2003.

    Ministry of Construction of the PRC. Code for design of steel structures [M]. Beijing:China Planning Press, 2003. (in Chinese)

    [18] 趙風(fēng)華, 齊永勝. 鋼結(jié)構(gòu)原理與設(shè)計 [M]. 重慶:重慶大學(xué)出版社, 2010.

    ZHAO F H, QI Y S. Theory and design of steel structure [M]. Chongqing:Chongqing University Press, 2010. (in Chinese)

    [19] 舒啟軍. 細(xì)長柱動力二階效應(yīng)非線性分析及振動臺試驗(yàn)設(shè)計 [D]. 重慶:重慶大學(xué), 2011.

    SHU Q J. Dynamic nonlinear analysis of second-order effect and shaking table test design of slender column [D]. Chongqing:Chongqing University, 2011. (in Chinese)

    [20] SOUZA D M D. Force-based finite element for large displacement inelastic analysis of frames [D]. University of California, Berkeley, 2000.

    [21] 張勇. 水平地震下高橋墩單自由度體系二階效應(yīng)研究 [D]. 重慶:重慶大學(xué), 2010.

    ZHANG Y. Study on the second-order effect of SDOF system of high bridge pier under horizontal earthquake [D]. Chongqing:Chongqing University, 2010. (in Chinese)

    [22] 羅小龍, 周鋒, 童樂為, 等. 滾軸支座滾動摩擦性能的試驗(yàn)研究 [J]. 結(jié)構(gòu)工程師, 2011, 27(6): 95-99.

    LUO X L, ZHOU F, TONG L W, et al. Experimental study on rolling friction behavior of roller bearing [J]. Structural Engineers, 2011, 27(6): 95-99. (in Chinese)

    [23] 吳二軍, 郭彤, 李愛群. 滾軸滾動摩擦系數(shù)研究與建筑物遷移工程水平動力計算 [J]. 土木工程學(xué)報, 2006, 39(6): 57-60.

    WU E J, GUO T, LI A Q. A study on the rolling friction coefficient and the horizontal driving-force calculation in building monolithic moving engineering [J]. China Civil Engineering Journal, 2006, 39(6): 57-60. (in Chinese)

    (編輯 胡玲)

    猜你喜歡
    摩擦力
    計算摩擦力的兩種方法
    神奇的摩擦力
    有一種神奇的力量
    自不量力
    傳送帶上的摩擦力問題例析
    基于摩擦力突變問題的四種模型
    神奇的摩擦力
    “三步法”突破高中物理“摩擦力”難題
    摩擦力做功中需弄清的幾個問題
    摩擦力的三個誤區(qū)釋疑
    99久久人妻综合| 久久人妻熟女aⅴ| 乱人伦中国视频| 99久久人妻综合| 日韩制服丝袜自拍偷拍| 亚洲图色成人| 18禁观看日本| 99国产精品一区二区三区| 高清欧美精品videossex| 中文字幕最新亚洲高清| 中文字幕精品免费在线观看视频| 热re99久久精品国产66热6| 久久狼人影院| 老司机影院毛片| 十八禁人妻一区二区| 国产成人91sexporn| 一级毛片 在线播放| 欧美大码av| 精品少妇内射三级| 超碰97精品在线观看| 午夜影院在线不卡| 90打野战视频偷拍视频| 欧美精品一区二区免费开放| 亚洲精品国产一区二区精华液| svipshipincom国产片| 精品久久久精品久久久| 91九色精品人成在线观看| 麻豆乱淫一区二区| 秋霞在线观看毛片| 亚洲欧美清纯卡通| 美女国产高潮福利片在线看| 亚洲av综合色区一区| 嫩草影视91久久| 999久久久国产精品视频| 两个人免费观看高清视频| 国产精品二区激情视频| 国产有黄有色有爽视频| 中文字幕人妻丝袜一区二区| 免费女性裸体啪啪无遮挡网站| 纵有疾风起免费观看全集完整版| 男女之事视频高清在线观看 | 2018国产大陆天天弄谢| 久久女婷五月综合色啪小说| 少妇裸体淫交视频免费看高清 | 你懂的网址亚洲精品在线观看| 丝袜美腿诱惑在线| 国产伦人伦偷精品视频| 嫁个100分男人电影在线观看 | 热re99久久精品国产66热6| 亚洲精品国产区一区二| 国产一区二区激情短视频 | 色婷婷av一区二区三区视频| 午夜福利影视在线免费观看| 亚洲少妇的诱惑av| 国产有黄有色有爽视频| 欧美精品一区二区大全| 熟女少妇亚洲综合色aaa.| netflix在线观看网站| 中文字幕制服av| 久久久久国产精品人妻一区二区| 午夜福利视频在线观看免费| 国产福利在线免费观看视频| 国产亚洲午夜精品一区二区久久| 热99国产精品久久久久久7| 国产日韩欧美亚洲二区| 亚洲av电影在线进入| 国产成人av教育| 免费一级毛片在线播放高清视频 | 大码成人一级视频| 老司机在亚洲福利影院| 免费在线观看完整版高清| 视频在线观看一区二区三区| 国产一区二区 视频在线| 又紧又爽又黄一区二区| 欧美日本中文国产一区发布| 精品熟女少妇八av免费久了| 男人舔女人的私密视频| www.av在线官网国产| 亚洲九九香蕉| 国产成人91sexporn| www.熟女人妻精品国产| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀 | 欧美精品亚洲一区二区| 欧美黑人精品巨大| 只有这里有精品99| 成人三级做爰电影| 久久久久久久精品精品| 欧美日韩视频高清一区二区三区二| 婷婷丁香在线五月| 色婷婷久久久亚洲欧美| 中文字幕亚洲精品专区| 蜜桃国产av成人99| 亚洲第一青青草原| 午夜福利免费观看在线| 精品国产乱码久久久久久男人| 欧美在线黄色| a级片在线免费高清观看视频| 国产日韩欧美亚洲二区| 免费看不卡的av| 亚洲第一av免费看| 捣出白浆h1v1| 欧美精品一区二区大全| 亚洲欧美一区二区三区黑人| 波多野结衣一区麻豆| 国产亚洲午夜精品一区二区久久| 狂野欧美激情性xxxx| 黄色毛片三级朝国网站| 99久久人妻综合| 日韩电影二区| 久久中文字幕一级| 亚洲欧洲精品一区二区精品久久久| 99九九在线精品视频| 丝袜喷水一区| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区大全| 亚洲国产欧美在线一区| 亚洲成人手机| 国产免费又黄又爽又色| 视频区欧美日本亚洲| 少妇人妻久久综合中文| 老司机在亚洲福利影院| 亚洲男人天堂网一区| 久久九九热精品免费| 久久精品久久精品一区二区三区| 国产日韩欧美视频二区| 90打野战视频偷拍视频| 亚洲精品国产色婷婷电影| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 国产一区二区在线观看av| 91老司机精品| 97精品久久久久久久久久精品| 狠狠精品人妻久久久久久综合| 日本五十路高清| 一级毛片黄色毛片免费观看视频| 亚洲精品国产一区二区精华液| 欧美乱码精品一区二区三区| 国产主播在线观看一区二区 | 国产精品99久久99久久久不卡| 久久精品aⅴ一区二区三区四区| 久久久久久久久久久久大奶| 午夜免费鲁丝| 久久精品久久久久久久性| 欧美变态另类bdsm刘玥| 亚洲天堂av无毛| 午夜日韩欧美国产| 亚洲视频免费观看视频| 国产国语露脸激情在线看| 一本综合久久免费| 宅男免费午夜| 日本wwww免费看| 亚洲色图综合在线观看| 欧美在线黄色| 中文字幕人妻熟女乱码| 水蜜桃什么品种好| 精品少妇黑人巨大在线播放| e午夜精品久久久久久久| 丝袜脚勾引网站| 亚洲,欧美精品.| 国产欧美日韩综合在线一区二区| 午夜福利影视在线免费观看| 亚洲成av片中文字幕在线观看| 免费观看a级毛片全部| 丁香六月欧美| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久| 亚洲国产最新在线播放| 另类亚洲欧美激情| 国产成人91sexporn| 免费日韩欧美在线观看| 国产一区二区三区综合在线观看| 国产精品亚洲av一区麻豆| 一级a爱视频在线免费观看| 十八禁网站网址无遮挡| 一本—道久久a久久精品蜜桃钙片| 国产成人av激情在线播放| 亚洲av综合色区一区| 男女床上黄色一级片免费看| 久久精品国产亚洲av高清一级| 国产国语露脸激情在线看| 在线观看国产h片| 午夜激情av网站| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 午夜老司机福利片| 亚洲色图综合在线观看| 蜜桃国产av成人99| 国产精品二区激情视频| 欧美成狂野欧美在线观看| 色视频在线一区二区三区| 精品少妇一区二区三区视频日本电影| 久久久国产欧美日韩av| 午夜福利在线免费观看网站| 国产精品秋霞免费鲁丝片| 在线 av 中文字幕| 国产在线视频一区二区| 午夜久久久在线观看| 999精品在线视频| 男女之事视频高清在线观看 | 国产免费又黄又爽又色| 精品熟女少妇八av免费久了| 99热全是精品| 国产亚洲欧美在线一区二区| 国产片内射在线| 18禁观看日本| 国产色视频综合| 我要看黄色一级片免费的| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 男女边摸边吃奶| 看免费av毛片| 国产午夜精品一二区理论片| 最新在线观看一区二区三区 | 亚洲国产精品一区二区三区在线| 亚洲精品av麻豆狂野| 国产免费又黄又爽又色| 久久精品久久久久久噜噜老黄| 不卡av一区二区三区| 国产成人一区二区在线| 欧美另类一区| 涩涩av久久男人的天堂| 99精品久久久久人妻精品| 在线 av 中文字幕| 久久女婷五月综合色啪小说| 脱女人内裤的视频| 国产91精品成人一区二区三区 | 亚洲人成电影观看| 欧美精品啪啪一区二区三区 | 成人亚洲精品一区在线观看| 多毛熟女@视频| 欧美成人精品欧美一级黄| 国产精品av久久久久免费| 国产精品一区二区在线不卡| 欧美97在线视频| 久久精品亚洲av国产电影网| 狂野欧美激情性bbbbbb| 高清黄色对白视频在线免费看| 9色porny在线观看| av欧美777| 亚洲伊人久久精品综合| 黄色视频在线播放观看不卡| 免费观看人在逋| 免费在线观看日本一区| 国产真人三级小视频在线观看| 一区二区三区精品91| 又黄又粗又硬又大视频| 女人久久www免费人成看片| 日韩av免费高清视频| 纯流量卡能插随身wifi吗| 伊人亚洲综合成人网| 午夜91福利影院| 一本一本久久a久久精品综合妖精| 人人妻人人澡人人看| 国产亚洲av高清不卡| 国产精品av久久久久免费| 一级毛片女人18水好多 | 中国国产av一级| 午夜免费成人在线视频| 99精品久久久久人妻精品| 久久国产精品影院| 男女下面插进去视频免费观看| 国精品久久久久久国模美| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| 亚洲精品国产av成人精品| 一本一本久久a久久精品综合妖精| 一本色道久久久久久精品综合| 亚洲中文字幕日韩| 色精品久久人妻99蜜桃| 精品福利观看| 51午夜福利影视在线观看| 脱女人内裤的视频| www.999成人在线观看| 高清av免费在线| 亚洲专区中文字幕在线| 免费观看a级毛片全部| 曰老女人黄片| 国产有黄有色有爽视频| 青青草视频在线视频观看| 丰满饥渴人妻一区二区三| 欧美日韩亚洲综合一区二区三区_| 国产精品 欧美亚洲| 无遮挡黄片免费观看| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 91国产中文字幕| 一级毛片女人18水好多 | 国产人伦9x9x在线观看| 午夜福利一区二区在线看| 成人国产av品久久久| 黄片播放在线免费| 国产欧美日韩一区二区三 | 亚洲国产中文字幕在线视频| 国产免费又黄又爽又色| 国产成人免费观看mmmm| 中文字幕亚洲精品专区| 成年人黄色毛片网站| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 婷婷色综合www| 日本黄色日本黄色录像| av天堂久久9| 在线观看免费视频网站a站| 首页视频小说图片口味搜索 | xxx大片免费视频| 免费在线观看影片大全网站 | 超碰成人久久| 免费女性裸体啪啪无遮挡网站| 国精品久久久久久国模美| 国产欧美日韩一区二区三 | www.熟女人妻精品国产| 人妻 亚洲 视频| 成人国语在线视频| 亚洲专区中文字幕在线| 黄色片一级片一级黄色片| 久久久久久免费高清国产稀缺| 在线看a的网站| 一级毛片我不卡| 亚洲成av片中文字幕在线观看| 国产亚洲欧美精品永久| 免费少妇av软件| 免费在线观看黄色视频的| 婷婷色av中文字幕| 亚洲av国产av综合av卡| 国产黄色视频一区二区在线观看| 国产成人精品久久二区二区91| 好男人视频免费观看在线| 亚洲,欧美,日韩| 久久久精品94久久精品| 国产真人三级小视频在线观看| 精品一区二区三区四区五区乱码 | 人成视频在线观看免费观看| 女性被躁到高潮视频| 香蕉丝袜av| 亚洲国产欧美网| 性高湖久久久久久久久免费观看| av在线播放精品| a级片在线免费高清观看视频| 国产亚洲欧美精品永久| 一区福利在线观看| 国产男女内射视频| 国产真人三级小视频在线观看| 色婷婷av一区二区三区视频| 精品亚洲成国产av| 成年女人毛片免费观看观看9 | 91国产中文字幕| 欧美成人午夜精品| 国产色视频综合| 日韩中文字幕欧美一区二区 | 老司机靠b影院| 亚洲美女黄色视频免费看| 建设人人有责人人尽责人人享有的| 伊人亚洲综合成人网| a级毛片黄视频| 日本a在线网址| 2021少妇久久久久久久久久久| 亚洲成色77777| 国产精品一区二区免费欧美 | a级毛片在线看网站| 亚洲一区中文字幕在线| 亚洲国产精品一区二区三区在线| 高清欧美精品videossex| 亚洲人成77777在线视频| 中文字幕高清在线视频| 美女脱内裤让男人舔精品视频| 亚洲色图综合在线观看| 色94色欧美一区二区| 丝袜美足系列| 丝袜美腿诱惑在线| 老司机影院成人| 国产熟女欧美一区二区| 一本大道久久a久久精品| 少妇裸体淫交视频免费看高清 | 日韩 亚洲 欧美在线| 电影成人av| 久热爱精品视频在线9| 电影成人av| 国产精品免费视频内射| 国产精品偷伦视频观看了| 亚洲三区欧美一区| 亚洲熟女精品中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国产极品粉嫩免费观看在线| 久久久久久人人人人人| 黄色片一级片一级黄色片| 国产深夜福利视频在线观看| 性少妇av在线| 巨乳人妻的诱惑在线观看| 桃花免费在线播放| 在线观看一区二区三区激情| 高清av免费在线| 搡老岳熟女国产| 午夜激情久久久久久久| 欧美在线黄色| 国产片内射在线| a级片在线免费高清观看视频| 久久久久精品国产欧美久久久 | 国产在线免费精品| 天天躁日日躁夜夜躁夜夜| 欧美激情高清一区二区三区| 在线 av 中文字幕| 国产精品久久久久久人妻精品电影 | 午夜免费观看性视频| 91麻豆av在线| 免费看不卡的av| 久久ye,这里只有精品| 亚洲国产精品999| 黄片播放在线免费| 视频区欧美日本亚洲| 天天躁日日躁夜夜躁夜夜| 国产亚洲精品久久久久5区| av不卡在线播放| 天天操日日干夜夜撸| 9热在线视频观看99| 久久久久久久久久久久大奶| 精品国产超薄肉色丝袜足j| 波多野结衣一区麻豆| 天天躁夜夜躁狠狠躁躁| 操美女的视频在线观看| avwww免费| 搡老乐熟女国产| 国产又色又爽无遮挡免| av天堂久久9| 一二三四社区在线视频社区8| 天堂8中文在线网| 九草在线视频观看| 男男h啪啪无遮挡| 国产高清不卡午夜福利| 日本欧美国产在线视频| 婷婷色麻豆天堂久久| 一级黄色大片毛片| 亚洲熟女精品中文字幕| 久久久久久人人人人人| 日韩大片免费观看网站| 男人添女人高潮全过程视频| 欧美精品av麻豆av| 捣出白浆h1v1| 国产免费一区二区三区四区乱码| 最近中文字幕2019免费版| 人妻人人澡人人爽人人| 欧美在线一区亚洲| 另类精品久久| 国产成人91sexporn| 国产亚洲av片在线观看秒播厂| 一级毛片 在线播放| 国产一区有黄有色的免费视频| 99热网站在线观看| 久久精品aⅴ一区二区三区四区| 日韩中文字幕视频在线看片| 精品人妻1区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲伊人色综图| 成年女人毛片免费观看观看9 | 国产成人一区二区三区免费视频网站 | 九草在线视频观看| 飞空精品影院首页| 免费不卡黄色视频| av欧美777| 青青草视频在线视频观看| 99re6热这里在线精品视频| 中国国产av一级| 一二三四社区在线视频社区8| 啦啦啦视频在线资源免费观看| 免费高清在线观看日韩| 欧美日韩精品网址| 国产精品一国产av| 精品人妻1区二区| 国产男女内射视频| 日本黄色日本黄色录像| 高清欧美精品videossex| 人成视频在线观看免费观看| 亚洲精品中文字幕在线视频| 日韩制服骚丝袜av| 国产成人系列免费观看| 制服人妻中文乱码| 女性生殖器流出的白浆| a级毛片黄视频| 亚洲美女黄色视频免费看| 亚洲av在线观看美女高潮| 女人精品久久久久毛片| 午夜免费男女啪啪视频观看| 亚洲欧美日韩高清在线视频 | a级片在线免费高清观看视频| 亚洲精品久久午夜乱码| 成在线人永久免费视频| 久久国产精品影院| 国产成人一区二区三区免费视频网站 | 午夜福利在线免费观看网站| 黄色视频在线播放观看不卡| 99香蕉大伊视频| 国产精品久久久av美女十八| 国产日韩一区二区三区精品不卡| 久久鲁丝午夜福利片| 亚洲国产最新在线播放| 男的添女的下面高潮视频| 韩国高清视频一区二区三区| 亚洲欧美成人综合另类久久久| 少妇精品久久久久久久| 丰满迷人的少妇在线观看| 国产成人精品在线电影| 国产黄色视频一区二区在线观看| 一区在线观看完整版| 少妇精品久久久久久久| 黑人欧美特级aaaaaa片| 欧美日韩综合久久久久久| 热99国产精品久久久久久7| 国产精品偷伦视频观看了| 精品亚洲乱码少妇综合久久| 国产精品人妻久久久影院| 亚洲av美国av| 亚洲欧美一区二区三区久久| 又大又爽又粗| 新久久久久国产一级毛片| 欧美精品av麻豆av| 色网站视频免费| 下体分泌物呈黄色| 久久精品久久久久久噜噜老黄| 波多野结衣av一区二区av| 成人影院久久| 叶爱在线成人免费视频播放| 黄网站色视频无遮挡免费观看| 人人妻人人澡人人看| 一区福利在线观看| 真人做人爱边吃奶动态| 亚洲黑人精品在线| 91字幕亚洲| 国产精品 欧美亚洲| 成人免费观看视频高清| 一二三四在线观看免费中文在| 观看av在线不卡| 精品国产一区二区三区四区第35| 久久久久久久久免费视频了| 黄频高清免费视频| 亚洲三区欧美一区| 国产成人欧美在线观看 | 午夜福利视频精品| 精品免费久久久久久久清纯 | 18禁黄网站禁片午夜丰满| h视频一区二区三区| 久久国产亚洲av麻豆专区| 99国产精品一区二区蜜桃av | 国产精品国产三级专区第一集| 人妻人人澡人人爽人人| 国产色视频综合| 午夜影院在线不卡| 久久精品久久精品一区二区三区| 成年美女黄网站色视频大全免费| 天堂俺去俺来也www色官网| 操美女的视频在线观看| 人妻一区二区av| 色播在线永久视频| 只有这里有精品99| 涩涩av久久男人的天堂| 人人澡人人妻人| 国产成人精品久久久久久| 日韩一区二区三区影片| 免费不卡黄色视频| 欧美人与善性xxx| 欧美人与性动交α欧美软件| 日本色播在线视频| 成年动漫av网址| 日日夜夜操网爽| 亚洲国产欧美一区二区综合| 国产在线一区二区三区精| 黑人猛操日本美女一级片| 黑丝袜美女国产一区| 韩国高清视频一区二区三区| 一级毛片 在线播放| 亚洲精品一二三| 丰满少妇做爰视频| 国产精品国产三级专区第一集| 精品亚洲乱码少妇综合久久| 久久久久精品国产欧美久久久 | 久久久久久久久久久久大奶| 久久精品国产a三级三级三级| 欧美另类一区| 波野结衣二区三区在线| 国产精品一二三区在线看| e午夜精品久久久久久久| 亚洲成av片中文字幕在线观看| 亚洲自偷自拍图片 自拍| 婷婷色综合大香蕉| 成人午夜精彩视频在线观看| 热re99久久国产66热| 嫁个100分男人电影在线观看 | 大片电影免费在线观看免费| videosex国产| 国产亚洲av高清不卡| 可以免费在线观看a视频的电影网站| 18禁国产床啪视频网站| 欧美国产精品va在线观看不卡| 亚洲美女黄色视频免费看| 久久精品熟女亚洲av麻豆精品| 中文字幕av电影在线播放| 一级,二级,三级黄色视频| 在线观看人妻少妇| 欧美日韩视频高清一区二区三区二| 性色av乱码一区二区三区2| 每晚都被弄得嗷嗷叫到高潮| 丰满人妻熟妇乱又伦精品不卡| 亚洲av电影在线进入| 韩国高清视频一区二区三区| 各种免费的搞黄视频| 日韩,欧美,国产一区二区三区| 亚洲欧美中文字幕日韩二区| 国产成人av教育| 欧美黑人欧美精品刺激| 97人妻天天添夜夜摸| 99久久精品国产亚洲精品| 十八禁网站网址无遮挡| 少妇 在线观看| 大码成人一级视频| 亚洲人成77777在线视频|