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

過焊孔對H型鋼柱梁節(jié)點受力性能的影響

2021-12-02 11:25:56唐祖宏包恩和
科學技術與工程 2021年33期
關鍵詞:變形

曾 榕,唐祖宏,尹 霞,包恩和

(1.桂林理工大學土木與建筑工程學院,桂林 541004;2.廣西新能源與建筑節(jié)能重點實驗室,桂林 541004;3.貴州省巖土力學與工程安全重點實驗室,貴陽 550025)

鋼框架結構進行抗震設計時,需考慮“強柱弱梁”基本理念;根據《建筑抗震設計規(guī)范》(GB 50011—2010)[1],鋼框架在不同抗震等級(一級、二級、三級)時,其柱梁節(jié)點的強柱系數分別為1.15、1.10、1.05;根據《鋼結構建筑的抗震規(guī)定》(AISC—2016)[2],位于H型鋼柱強軸一側的節(jié)點,其強柱系數應大于1.1;文獻[3]指出為避免與“強柱弱梁”基本理念沖突,結構分析時不能把節(jié)點域作為框架的基本耗能結構。為了明確過焊孔對鋼框架柱梁節(jié)點的力學性能的影響,眾多學者進行了一系列研究,從不同節(jié)點域強度對H型鋼框架屈服強度、初始剛度、變形性能和極限承載力的影響[4-6],到構件連接方式對節(jié)點域的塑性和彈性剪切應變值的影響[7]。文獻[8]提出節(jié)點域剛度計算公式;鋼框架梁腹板較薄時,過焊孔底部應力集中松弛[9];研究過焊孔尺寸對鋼框架柱梁節(jié)點的影響[10];研究三種類型過焊孔在沖擊荷載下對鋼框架梁柱節(jié)點承載能力和變形性能的影響[11]。基于節(jié)點域影響的柱梁節(jié)點研究,研發(fā)了一種改進型開縫腹板梁的特殊抗彎工字形柱梁節(jié)點[12]。研究不等高鋼梁框架節(jié)點的地震損傷演化過程和力學性能[13-14]。通過有限元分析不同厚度節(jié)點域的變形能力,并提出了節(jié)點域的抗剪承載力計算公式[15];在柱中設計加強板來控制節(jié)點域的剪切變形[16]。通過有限元分析不同過焊孔對非加強型翼緣焊接-腹板焊接(welded unreinforced flange-welded web,WUF-W)節(jié)點循環(huán)特性的影響[17];通過加載試驗和有限元分析,研究了節(jié)點域內腹板和翼緣的厚度對連接件的抗剪性能的影響[18]。

綜上,基于過焊孔對H型鋼柱梁節(jié)點力學性能研究少的背景,且以往研究多以節(jié)點域高寬比,寬厚比為參量,現以節(jié)點域的強度為主要參量,設計5組材料為H型鋼、形狀為T形的柱梁節(jié)點試件,進行低周往復加載試驗,并建立模型進行有限元分析,來研究節(jié)點域強度變化和過焊孔對H型鋼柱梁節(jié)點力學性能的影響。

1 試件概況

1.1 試件設計

H型鋼柱梁節(jié)點試件的形狀為T形,梁和柱都采用H型鋼,其中梁為窄翼緣,柱為寬翼緣,梁尺寸(HN-400 mm×200 mm×8 mm×13 mm),柱尺寸(HW-300 mm×300 mm×10 mm×15 mm)。在框架節(jié)點處焊接不同厚度的鋼板來調整節(jié)點域強度,設計5個不同節(jié)點域強度的試件。試件柱梁節(jié)點的過焊孔采用雙弧線型,弧線半徑為R=35 mm及r=10 mm,過焊孔和柱梁翼緣焊接詳情如圖1所示。各試件均采用Q355鋼和E50焊條,柱翼緣與梁翼緣連接、節(jié)點域的水平加勁肋和加厚板與柱連接均采用對接焊縫形式,而梁腹板與柱翼緣連接采用角接焊縫形式。各試件的柱、梁、節(jié)點域的水平加勁板及加厚板等的參數如表1所示。

圖1 試驗過焊孔和柱梁翼緣焊接Fig.1 Weld access hole and welding of test beam-column joints

1.2 材料性質

梁和柱的翼緣板、腹板以及加勁肋等處各切取3組標準試驗片,進行常溫拉伸試驗后取強度和變形的平均值,如表2所示。

1.3 試件的柱、梁及節(jié)點域強度

表1 試件組成Table 1 Specimens composition

表2 材料力學性能試驗結果Table 2 Test results of material mechanical properties

(1)

δ1~δ9為布置的9個位移計,δ1測梁端水平位移,δ2測柱正面翼緣水平位移,δ3測柱背面翼緣水平位移,δ4、δ5測梁上下翼緣豎直位移,δ6、δ7測節(jié)點域對角線方向變形,δ8、δ9分別測柱左右兩端支座處的豎直位移,所測得位移分別記作δ1~δ9圖2 試件固定、加載及位移計布置Fig.2 Specimens fixation,loading and displacement meter arrangement

P為荷載;H為柱高;L為梁長;hb為節(jié)點域高度;hc為節(jié)點域寬度圖3 試件計算簡圖Fig.3 Calculation diagram of specimens

(2)

(3)

節(jié)點域與梁強度比為

(4)

強柱系數為

(5)

式中:lb、lc分別為梁端、柱端到節(jié)點域中心(柱梁軸線交點)的距離;hb、hc分別為節(jié)點域高度、寬度。

根據表1、表2及式(1)~式(5),計算試件各組成部分的強度、剪力、節(jié)點域與梁強度比和節(jié)點的強柱系數,試件的命名取決于節(jié)點域與梁強度比,用T-①~T-⑤表示,如表3所示。

表3 各試件命名和柱、梁及節(jié)點域強度Table 3 Specimen name and the strength of column,beam and panel zone

2 試驗加載方式與檢測

2.1 加載方式

為了便于加載,采用“橫柱豎梁”的方式(圖2)。柱高為3 m,兩端采用鉸接方式;梁的半跨長度為3.5 m,加載機制基于位移控制法,將低周往復荷載作用在梁端,用T形試件的位移角控制,且考慮抗震規(guī)范的結構層間彈塑性極限位移角為0.02 rad,為更符合強震下結構節(jié)點實際工況的柱梁節(jié)點變形及耗能等的判斷,控制位移角為±0.02、±0.04、±0.06 rad的三段,柱梁節(jié)點的位移角控制為結構層間彈塑性極限位移角的3倍,各段往復作用兩次,試件如果沒有達到最大荷載,就按±0.06 rad持續(xù)加載,當試件達到其最大荷載值的85%左右或加載不提升才停止加載。

注:Vp表示節(jié)點域剪力。

2.2 試件變形檢測及梁、柱和節(jié)點域變形幾何關系

試件梁、柱及節(jié)點域變形幾何關系如圖4所示。試驗中試件各組成部的轉角計算公式如下。

(1)試件位移角。

(6)

(2)梁轉角。

(7)

梁變形引起的梁端轉角θb如圖4(a)所示,公式為

(8)

bθ為梁轉角;θb為梁變形引起的梁端轉角;cθ為柱轉角;θc為柱變形引起的梁端轉角;γ為節(jié)點域轉角;θp為節(jié)點域剪切變形引起的梁端轉角圖4 試件位移角與梁、柱及節(jié)點域變形角的幾何關系Fig.4 Geometric relationship between displacement angles of specimens and deformation angles of beams,columns and panel zone

(3)柱轉角。

(9)

柱變形引起的梁端轉角θc如圖4(b)所示,公式為

(10)

(4)節(jié)點域轉角。

(11)

節(jié)點域剪切變形引起的梁端轉角θp如圖4(c)所示,公式為

(12)

3 有限元分析

為進一步把握節(jié)點域強度變化對試件各組成部受力-變形影響,采用ABAQUS平臺建立6組有限元分析模型。其中5組分析模型的各參量與試件相同,分別命名為分析模型①~模型⑤,簡稱AM-①~AM-⑤。本次試驗為考察各試件節(jié)點域強度與一般工程中假設的剛性節(jié)點域區(qū)別,建立1組剛性節(jié)點域分析模型,命名為分析模型⑥,簡稱AM-⑥。有限元分析模型考慮試件加載、固定情況及對稱性等因素,建立1/2模型,建模單元采用4點平面要素和2點梁要素,分析模型材料的本構關系基于母材試驗,服從Von Mises屈服準則,塑域中的本構關系為各向相同硬化;有限元分析模型(圖5)為不考慮反復加載的硬化作用對分析模型的過焊孔周圍受力的影響,單一分析過焊孔對梁端受力的影響,分析模型采用單向加載。

圖5 有限元分析模型Fig.5 Finite element analysis model

4 研究結果

4.1 試件荷載與變形關系及最終狀態(tài)

用梁端彎矩(Mb)作每個試件的基準荷載,把測量的結果,按式(6)~式(12)計算各個試件的位移角和柱、梁、節(jié)點域發(fā)生變形引起的梁端變形角。框架產生塑性變形的部位只有T-①~T-④4個試件的梁和節(jié)點域以及T-⑤試件的梁,5個試件的柱都是彈性狀態(tài)。

圖6(a)、圖6(d)和圖6(g)中虛線為梁全截面屈服彎矩計算值。試件T-① 的基準荷載-變形關系如圖6(a)~圖6(c)所示;試件T-③ 的基準荷載-變形關系如圖6(d)~圖6(f)所示;試件T-⑤ 的基準荷載-變形關系如圖6(g)~圖6(h)所示。

圖6(a)~圖6(c)中,把試件T-①的位移角控制在0.06 rad,加載至第5次正向循環(huán),梁端彎矩的最大值為478 kN·m;加載至相同位移角第13次正向循環(huán),當梁端彎矩達到最大值的85%左右,停止加載。T-①的主要變形由節(jié)點域變形引起,節(jié)點域變形占總累積塑性變形位移角(2.34 rad)的90%。

圖6(d)~圖6(f)中,把試件T-③的位移角控制在0.06 rad,加載至第1次正向循環(huán),梁翼緣發(fā)生局部屈曲,加載至相同位移角第1次負向循環(huán),梁端彎矩的最大值為607 kN·m;在相同位移角下試件繼續(xù)加載至第2次負向循環(huán),梁端局部屈曲繼續(xù)發(fā)展,當梁端彎矩達到最大值的85%左右,停止加載。T-③的主要變形由梁和節(jié)點域變形引起,梁和節(jié)點域變形分別占總累積塑性變形位移角(0.47 rad)的64%和36%。

圖6(g)~圖6(h)中,把試件T-⑤的位移角控制在0.04 rad,加載至第2次負向循環(huán)加載,梁翼緣發(fā)生局部屈曲變形,梁端彎矩的最大值為615 kN·m;繼續(xù)加載至0.06 rad第1次負向循環(huán),梁端局部屈曲繼續(xù)發(fā)展,當梁端彎矩達到最大值的85%左右,停止加載。T-⑤的塑性變形由梁變形引起,梁變形占總累積塑性變形位移角(0.3 rad)的100%,節(jié)點域保持彈性狀態(tài)。

圖6 試件的荷載-變形關系Fig.6 Load-deformation curve of specimens

增大節(jié)點域的強度,不利于試件的變形性能,但能提高試件的最大承載力。

試驗結束后,各試件裂縫主要分布在梁和引弧板之間的人工縫附近、梁翼緣外側中心厚度方向及過焊孔兩端。試件T-①的節(jié)點域產生剪切屈曲,如圖7(a)所示;試件T-③的梁端產生局部屈曲,如圖7(b)所示;試件T-⑤的梁端產生較嚴重的局部屈曲,如圖7(c)所示。

圖7 試件最終狀態(tài)Fig.7 The final state picture of specimen

4.2 試件的試驗與分析結果比較

基于試件T-①的正向滯回曲線繪制骨架曲線,與AM-①荷載-變形曲線對比,如圖8所示。

圖8 T-①與AM-①分析結果對比Fig.8 Results comparison between T-① and AM-①

分析結果和試驗結果的初始剛度吻合,而分析結果的屈服彎矩微高、第二剛度微小,塑域強度上升趨勢的分析與試驗結果較好吻合。可以利用有限元分析方法的優(yōu)勢追蹤試驗中無法把握的過焊孔的影響。

4.3 分析模型的梁受力分析

4.3.1 梁彎矩荷載分布

各分析模型在不同的梁端變形角比θb/pθb(θb為梁變形引起的梁端轉角;pθb為梁端屈服變形角)時梁腹板和梁翼緣彎矩荷載分布,如圖9和圖10所示。圖9和圖10的縱軸分別為梁腹板彎矩荷載(MW)和梁翼緣彎矩荷載(Mf),橫軸為梁端軸向0~1 m范圍;圖中虛線為梁腹板和梁翼緣全截面屈服彎矩計算值。圖9和圖10(a)~圖10(d)分別對應θb/pθb為0.2、1.4、3及5的狀態(tài)。

由圖9和圖10(a)、圖10(c)及圖10(d)可知,節(jié)點域強度的變化不影響梁腹板和梁翼緣的軸向彎矩分布曲線形狀;而圖9(b)和圖10(b)對應θb/pθb=1.4,節(jié)點域強度的變化影響梁彎矩荷載分布,隨著分析模型節(jié)點域強度的提升,梁腹板和翼緣承擔的彎矩荷載變大,尤其在柱梁節(jié)點過焊孔端部的彎矩荷載變化較大。

圖9 分析模型梁腹板彎矩荷載分布Fig.9 Bending moment distribution on the beam web of analysis model

圖10 分析模型梁翼緣彎矩荷載分布Fig.10 Bending moment distribution on the beam flange of analysis model

由圖9和圖10可知,柱梁節(jié)點的過焊孔會影響梁彎矩荷載分布曲線形狀,梁腹板軸向彎矩分布,因受過焊孔削弱截面的影響,梁端軸向0~80 mm范圍內的彎矩荷載相對小,其中,梁端軸向接近35 mm處的彎矩荷載最小;同時,因過焊孔端部易應力集中,梁端軸向0~10 mm及60~80 mm范圍內的彎矩荷載相對大。而梁翼緣軸向彎矩分布,梁端軸向0~80 mm范圍內,受節(jié)點過焊孔影響引起的梁腹板承載能力遞減、過焊孔周圍易應力集中及外載作用不變等原因,梁翼緣彎矩荷載相對大。

圖9和圖10(a)對應θb/pθb=0.2,AM-①~AM-⑥的變形在彈性范圍內,梁端軸向0~1 m范圍內作用的彎矩荷載遠小于梁翼緣及腹板的全截面屈服彎矩計算值。

圖9和圖10(b)對應θb/pθb=1.4,AM-①~AM-⑥的梁端發(fā)生塑變。節(jié)點域強度的變化影響梁彎矩荷載分布,隨著分析模型節(jié)點域強度的提升,梁腹板和翼緣承擔的彎矩荷載變大,尤其在柱梁節(jié)點過焊孔端部的彎矩荷載變化較大。梁腹板的最大彎矩荷載為115 kN·m左右,仍小于梁腹板全截面屈服彎矩計算值;AM-①~AM-⑥的梁翼緣承擔的彎矩荷載在梁端軸向0~80 mm范圍內,相對大,且大于梁翼緣全截面屈服彎矩計算值;梁端軸向80~350 mm范圍內,AM-④~AM-⑥處于塑性狀態(tài),其他分析模型保持彈性狀態(tài)。另外,強節(jié)點域分析模型AM-④、AM-⑤及AM-⑥的梁翼緣彎矩荷載分布基本一致。

圖9和圖10(c)對應θb/pθb=3.0,AM-②~AM-⑥的梁端發(fā)生較大的塑變。梁端軸向120~360 mm范圍內,梁腹板的彎矩荷載大于腹板全截面屈服彎矩計算值;AM-②~AM-④ 梁腹板的硬化大于AM-⑤和AM-⑥。梁端軸向0~360 mm范圍以外,隨著梁端軸向距離增大,荷載作用逐漸減小。而梁端軸向0~600 mm范圍內,AM-②~AM-⑥的梁翼緣彎矩荷載均大于梁翼緣全截面屈服彎矩計算值。梁端軸向0~600 mm范圍以外,隨著梁端軸向距離增大,梁翼緣彎矩荷載逐漸減小。另外,AM-①為節(jié)點域與梁強度比等于0.5的弱節(jié)點域模型,AM-①的變形集中于節(jié)點域,梁端變形角比θb/pθb小于3.0,不在此研究范圍內。

圖9和圖10(d)對應θb/pθb=5.0,AM-③~AM-⑥的梁彎矩荷載分布趨勢類似于圖9和圖10(c)對應θb/pθb=3.0時梁彎矩荷載分布趨勢。

4.3.2 梁端剪力荷載分析

圖11(a)中,AM-①為Rpb=0.5的弱節(jié)點域模型,θb/pθb<0.5時,AM-①處于彈性狀態(tài),梁端剪力荷載曲線和梁端腹板剪力荷載曲線基本保持直線形狀。當θb/pθb≥0.5時,AM-①的節(jié)點域進入塑域,在0.5<θb/pθb≤0.8范圍內,梁端腹板剪力荷載遞減;而0.8<θb/pθb≤1.0范圍內,梁端腹板剪力荷載有所提升。θb/pθb≤1.0時,AM-①的梁端剪力荷載逐漸增大,AM-⑥的梁端腹板剪力荷載保持直線形狀;當θb/pθb>1.0時,AM-①和AM-⑥的梁進入塑域,AM-①的梁端剪力荷載值和梁端腹板剪力荷載值及AM-⑥的梁端腹板剪力荷載值等基本保持不變。當θb/pθb值達到最大時,AM-①的梁端腹板剪力荷載值為AM-⑥的梁端腹板承擔剪力荷載值的1/3左右。另外,AM-②的梁端剪力荷載曲線和梁端腹板剪力曲線分布趨勢類似于AM-①。

圖11 分析模型的梁端剪力荷載-梁端變形角比Fig.11 Shear and deformation angle ratio curve of the beam-end in analysis model

AM-①、AM-③及AM-⑤的梁端翼緣剪力荷載值為梁端剪力荷載值與梁端腹板剪力荷載值之差;AM-⑥的節(jié)點域為剛性,梁端腹板剪力荷載不受節(jié)點域影響;因此,AM-⑥的梁端腹板剪力荷載值與AM-①、AM-③及AM-⑤的梁端腹板剪力荷載值之差為受節(jié)點域影響的梁端腹板剪力荷載值。

圖11(b)中,AM-③ 的Rpb=1.0,θb/pθb<1.0時,AM-③和AM-⑥處于彈性狀態(tài),梁端剪力荷載曲線和梁端腹板剪力荷載曲線基本保持直線形狀。當1.0 ≤θb/pθb≤ 2.0時,AM-③和AM-⑥的梁進入塑域,AM-③的梁端剪力荷載值呈遞增趨勢,而AM-③和AM-⑥的梁端腹板剪力荷載值呈遞減趨勢。當θb/pθb>2.0后,AM-③的梁端剪力荷載值繼續(xù)呈遞增趨勢,而AM-③和AM-⑥的梁端腹板剪力荷載值趨于平穩(wěn)。當θb/pθb值達到最大時,AM-③的梁端腹板剪力荷載值為AM-⑥的梁端腹板承擔剪力荷載值的2/3左右。另外,AM-④ 的梁端剪力荷載曲線和梁端腹板剪力曲線分布趨勢類似于AM-③。

圖11(c)中,AM-⑤的Rpb=1.4,θb/pθb<1.0時,AM-⑤和AM-⑥處于彈性狀態(tài),梁端剪力荷載曲線和梁端腹板剪力荷載曲線基本保持直線形狀。θb/pθb≥1.0時,AM-⑤的梁端剪力荷載值和梁端腹板剪力荷載值分布類似于AM-③;且AM-⑤的梁端腹板剪力荷載值大小與AM-⑥的梁端腹板剪力荷載值大小基本一致。由此可見,當Rpb=1.4時,AM-⑤可按剛性節(jié)點域考慮。

AM-①和AM-② 的屈服剪力取決于節(jié)點域強度,其屈服剪力值分別為85 kN和135 kN;AM-③的梁和節(jié)點域同時屈服,屈服剪力值為161 kN;AM-④和AM-⑤的屈服剪力取決于梁,其屈服剪力值也為161 kN。另外,節(jié)點域強度變化影響梁端腹板承擔的剪力荷載,隨著節(jié)點域強度的減小,梁端腹板承擔的剪力荷載也減小;例如,AM-①和AM-③ 的梁端腹板承擔剪力荷載值分別為AM-⑥的梁端腹板承擔剪力荷載值的1/3和2/3左右。

5 結論

基于節(jié)點域強度變化,且考慮過焊孔影響等因素,通過試驗和有限元分析,研究T形柱梁節(jié)點的各組成部的力學性能。結果如下。

(1)增大節(jié)點域的強度,不利于試件的變形性能,但能提高試件的最大承載力。

(2)模型處于彈性狀態(tài)或模型梁端變形角比較大(θb/pθb≥3)時,節(jié)點域強度變化對梁彎矩荷載分布無明顯的影響;而模型梁端變形角比等于1.4時,隨著節(jié)點域強度提升,梁彎矩荷載有增大趨勢,在弱節(jié)點域模型中此趨勢越明顯。

(3)柱梁節(jié)點的過焊孔會影響AM-①~AM-⑥模型的梁彎矩荷載分布。梁端軸向0~80 mm范圍內,受過焊孔削弱梁腹板截面的影響,梁腹板的彎矩荷載相對小,其中,梁端軸向接近35 mm處的彎矩荷載最小;而梁端軸向0~10 mm及80~120 mm范圍內的過焊孔端部,因應力易集中,梁腹板彎矩荷載偏大。在梁端軸向0~80 mm范圍內,受節(jié)點過焊孔影響引起的梁腹板承載能力遞減、過焊孔周圍易應力集中及外荷載作用不變等原因,梁翼緣彎矩荷載相對大;且梁腹板承擔的荷載越小,梁端軸向對應處的梁翼緣承擔的荷載越大。

(4)節(jié)點域強度變化影響AM-①~AM-⑥模型的梁端腹板承擔的剪力荷載,隨著節(jié)點域強度的減小,梁端腹板承擔的剪力荷載也減小。

猜你喜歡
變形
變形記
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
柯西不等式的變形及應用
“變形記”教你變形
不會變形的云
“我”的變形計
會變形的折紙
童話世界(2018年14期)2018-05-29 00:48:08
變形巧算
例談拼圖與整式變形
會變形的餅
主站蜘蛛池模板: 国产哺乳奶水91在线播放| m男亚洲一区中文字幕| 国产精品乱偷免费视频| 中文字幕av无码不卡免费| 国产精品自在自线免费观看| 国产乱子伦一区二区=| 精品国产Ⅴ无码大片在线观看81| 福利视频一区| 狂欢视频在线观看不卡| 国产精品无码一二三视频| 国产亚洲视频免费播放| 中文字幕首页系列人妻| 99精品在线看| A级毛片无码久久精品免费| 精品国产美女福到在线不卡f| 亚洲天堂免费在线视频| 国产成人麻豆精品| 免费人成黄页在线观看国产| 美女免费黄网站| 激情亚洲天堂| 91视频青青草| 美女啪啪无遮挡| 露脸真实国语乱在线观看| 色综合成人| 91最新精品视频发布页| 免费视频在线2021入口| 伊在人亞洲香蕉精品區| 又黄又湿又爽的视频| 毛片视频网| 国产拍揄自揄精品视频网站| 国产精品永久久久久| 欧美日韩一区二区在线免费观看| 欧美在线导航| 日本不卡在线| 97久久免费视频| 国产91高跟丝袜| 国产精品久久久久久久久久98 | 国产精品人人做人人爽人人添| 国产特一级毛片| 国产www网站| Aⅴ无码专区在线观看| 伊人福利视频| www欧美在线观看| 看你懂的巨臀中文字幕一区二区| 久久无码免费束人妻| 久久精品午夜视频| 国产一区免费在线观看| 91久久偷偷做嫩草影院免费看 | 国产成人一区二区| 国产综合欧美| 亚洲一级无毛片无码在线免费视频| 欧美一级特黄aaaaaa在线看片| 久久综合丝袜长腿丝袜| 亚洲av成人无码网站在线观看| 激情六月丁香婷婷四房播| 亚洲日本中文字幕乱码中文 | 亚洲欧美另类久久久精品播放的| 精品国产女同疯狂摩擦2| 国产高清精品在线91| 亚洲中文字幕无码mv| 特级aaaaaaaaa毛片免费视频| 久久美女精品国产精品亚洲| 精品国产中文一级毛片在线看 | 亚洲另类色| 超清无码熟妇人妻AV在线绿巨人| 热这里只有精品国产热门精品| 亚洲成aⅴ人在线观看| 性色生活片在线观看| 亚洲丝袜中文字幕| 日本国产在线| 中文国产成人久久精品小说| 在线高清亚洲精品二区| 色综合热无码热国产| 日韩无码视频网站| 欧美日本激情| 亚洲成人免费在线| 天堂在线亚洲| 免费av一区二区三区在线| 亚洲一区二区精品无码久久久| 2020极品精品国产| 99久久精品免费看国产免费软件| 无码网站免费观看|