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

    雪荷載作用下幾字型鋼日光溫室極限承載力分析

    2022-02-04 12:17:44姜迎春徐占洋白義奎王鐵良
    農業(yè)工程學報 2022年19期
    關鍵詞:承載力

    王 聰,姜迎春,徐占洋,張 峰,白義奎,王鐵良

    ·農業(yè)生物環(huán)境與能源工程·

    雪荷載作用下幾字型鋼日光溫室極限承載力分析

    王 聰1,姜迎春2,徐占洋1,張 峰1,白義奎1,王鐵良1※

    (1. 沈陽農業(yè)大學水利學院,沈陽 110866;2. 沈陽農業(yè)大學工程學院,沈陽 110866)

    雪災是導致日光溫室倒塌的主要原因之一。為探明幾字型鋼日光溫室在雪荷載作用下的失穩(wěn)機理,該研究采用有限單元法,以8 m跨日光溫室為研究對象,模擬其在雪荷載(均勻分布雪荷載和非均勻分布雪荷載)作用下的失穩(wěn)破壞過程,計算其極限承載力,并探究縱向系桿、初始幾何缺陷、截面參數(shù)對極限承載力的影響。結果表明:對于凈截面面積、上翼緣寬度、腹板高度和壁厚均相同的幾字型鋼和空心矩形鋼管,幾字型鋼日光溫室極限承載力稍高于空心矩形鋼管日光溫室極限承載力;相較于均勻分布雪荷載,日光溫室拱架對非均勻分布雪荷載更為敏感,非均勻分布雪荷載作用下的極限承載力約是均勻分布雪荷載作用下的28%,在日光溫室結構設計中,應重點考慮非均勻分布雪荷載工況;在非均勻分布雪荷載作用下,屋脊和后屋面支座處為危險截面,最先進入全截面屈服狀態(tài);縱向系桿的設置可有效抑制結構平面外變形,進而提高結構極限承載力,有縱向系桿約束條件下的結構極限承載力約是無縱向系桿約束條件下的1.25倍;該日光溫室拱架對初始幾何缺陷敏感度較低,當最大初始幾何缺陷幅值從5 mm增加到20 mm時,極限承載力降低約2%;在幾字型鋼截面選取時,在滿足規(guī)范要求寬厚比前提下,建議上翼緣寬度與翻邊寬度之比控制在4.17左右,腹板高度與翻邊寬度之比不大于9.25,下翼緣寬度與翻邊寬度之比不大于1.7,上翼緣寬度與下翼緣寬度之比控制在3.33左右,腹板高度與下翼緣寬度之比控制在4.67左右。該研究結果可為開口冷彎薄壁型鋼日光溫室拱架抗雪設計提供參考。

    溫室;承載力;荷載;有限單元法;失穩(wěn)破壞;冷彎薄壁型鋼;開口截面;截面優(yōu)化

    0 引 言

    日光溫室是一種具有中國特色、擁有自主知識產權的設施類型,具有結構簡單、節(jié)約能源、經(jīng)濟效益相對較高等優(yōu)點[1]。作為一種農業(yè)建筑,其不僅要給植物提供適宜的生長環(huán)境,還需具有抵抗風雪等外部荷載的能力[2]。雪災是北方地區(qū)主要氣象災害之一,2021年11月,遼寧出現(xiàn)特大暴雪天氣,造成大量日光溫室倒塌。從調查結果來看,日光溫室結構失效模式主要為失穩(wěn)破壞[3]。因此,了解日光溫室結構在雪荷載作用下的失穩(wěn)機理,對于提高日光溫室結構抗雪災能力具有重要意義。

    從日光溫室骨架形式來看,目前主要流行2種結構形式:桁架結構和單管結構。桁架結構將構件的內力轉化為桿件的軸力,提高了截面的利用效率,在承載能力方面具有顯著優(yōu)勢[4]。目前,國內學者針對桁架結構日光溫室的穩(wěn)定性能[5-7],結構優(yōu)化[8-10],節(jié)點設計[11]和動力響應[12-13]等方面進行了較多研究,這些研究結果為桁架結構日光溫室應用與推廣提供了理論依據(jù)。然而,桁架結構存在用鋼量大、標準化程度低、抗腐蝕性能差、焊接工作量大且質量不易保證等缺點。

    近年來,隨著日光溫室結構向輕簡化、裝配化的方向發(fā)展,單管結構日光溫室逐步開始發(fā)展并越來越多的應用在實際工程中[14-15]。相較于桁架結構,單管結構具有用鋼量少、標準化程度高、抗腐蝕性能好、安裝方便等優(yōu)點。目前,針對單管結構日光溫室受力性能的研究主要有:齊飛等[16]采用“全荷載組合”對平橢圓管日光溫室拱架進行結構分析;丁敏等[17-18]研究了覆蓋材料蒙皮效應對方管日光溫室整體穩(wěn)定性能的影響;閆冬梅等[19]運用3D3S軟件研究了拉桿布置位置對橢圓管日光溫室拱架內力分布的影響,并建議在不設置屋脊拉桿情況下,盡量把屋脊做成圓弧形;針對單管骨架截面剛度和強度較低等問題,齊飛等[20-21]基于臨時加固策略研究了前屋面豎向支撐對溫室骨架安全性能的影響,并給出了支撐數(shù)量和支撐位置的建議值。Wang等[22-24]運用ANSYS軟件對裝配落地式單管日光溫室拱架進行靜力特性與動力時程分析。白義奎等[25]對平橢圓單管日光溫室進行全尺寸加載試驗。以上研究成果為單管結構日光溫室防災減災和應用推廣提供了理論依據(jù)。

    值得指出的是,以上研究僅限于閉口截面的方管、矩形管或平橢圓管。相較于閉口截面鋼管,冷彎薄壁開口截面型鋼(如幾字型鋼)具有加工成本較低、制作工藝簡單、材料利用率較高等特點。因此,其在實際工程中得到廣泛應用[26-28]。然而,由于開口截面屬性,幾字型鋼構件抗扭能力弱,截面形心與剪心不重合,在豎向雪荷載作用下易發(fā)生失穩(wěn)破壞[29]。因此,探究幾字型鋼日光溫室拱架在雪荷載作用下的失穩(wěn)機理,計算其極限承載力,對于提高其抗雪災能力具有重要意義。

    為此,本研究以幾字型鋼日光溫室拱架為研究對象,采用ANSYS有限元分析軟件,模擬其在雪荷載作用下的失穩(wěn)破壞過程,分析其破壞模式,計算其極限承載力,并探究縱向系桿、初始幾何缺陷、截面參數(shù)對極限承載力的影響,從而尋求最優(yōu)截面尺寸參數(shù),為幾字型鋼日光溫室結構抗雪設計提供參考。

    1 溫室結構參數(shù)

    本文選取某跨度為8 m,脊高為3.8 m的日光溫室為研究對象,具體結構參數(shù)如圖1所示。前屋面水平投影6.5 m,后屋面水平投影1.5 m,后墻高2.3 m,后屋面角45°。為保證日光溫室整體穩(wěn)定性,在前屋面和屋脊處布置4道縱向系桿。溫室前后屋面拱架選用(××××)50 mm×70 mm×15 mm×12 mm×2 mm冷彎薄壁幾字型鋼,截面尺寸如圖2所示。

    圖1 日光溫室剖面示意圖

    注:a為下翼緣寬度,mm;b為上翼緣寬度,mm;c為翻邊寬度,mm;h為腹板高度,mm;t為壁厚,mm。

    2 有限元模型

    2.1 材料與單元類型

    日光溫室前后拱架均采用冷彎薄壁幾字型鋼,鋼材選用Q235,楊氏模量206 GPa,泊松比0.3,鋼材密度7 850 kg/m3,鋼材按理想彈塑性材料考慮。為考慮幾字型鋼翼緣和腹板可能發(fā)生局部屈曲問題,本研究采用Shell181殼單元模擬日光溫室拱架。同時,采用殼單元可精確模擬縱向系桿對日光溫室拱架側向約束位置。

    2.2 荷載與邊界條件

    日光溫室拱架自重由軟件自動計算,根據(jù)《農業(yè)溫室結構荷載規(guī)范》GB/T51183—2016[30],前屋面覆蓋15 mm厚發(fā)泡聚乙烯保溫被,荷載取0.007 kN/m2;后屋面采用100 mm厚金屬夾芯板,荷載取0.30 kN/m2;溫室內固定設備荷載取0.07 kN/m2;由于暫未確定吊掛方式,作物荷載按單位面積荷載計算,取0.15 kN/m2[30];作用于溫室拱架上翼緣雪荷載水平投影值按式(1)計算。

    式中k為作用于溫室拱架上翼緣雪荷載,kN/m2;r為屋面積雪分布系數(shù);t為加熱影響系數(shù),對于不加熱日光溫室,取值為1.0[30];0為基本雪壓,kN/m2。

    考慮到風致積雪漂移現(xiàn)象,本研究考慮2種雪荷載分布模式:均勻分布雪荷載和非均勻分布雪荷載。2種雪荷載分布模式下的積雪分布系數(shù)如圖3所示。對于均勻分布雪荷載:前屋面積雪分布系數(shù)r由前屋面水平投影(1=6.5 m)和脊高(=3.8 m)確定(r=2×1/(8×)=0.43);后屋面積雪分布系數(shù)r,b由后屋面角(45°)確定(r,b=0.8×(60°?)/30°=0.40)。對于非均勻分布雪荷載:日光溫室前屋面積雪分布系數(shù)呈三角形分布,且最大值在前屋面1/3處,覆蓋保溫被時,r,m最大值可取2.0[30];后屋面積雪分布系數(shù)為0.75r,b[30],計算得0.30。

    注:L1為前屋面水平投影長度,m;H為脊高,m;α為后屋面角度,(°);均勻分布雪荷載工況下,μr為前屋面積雪分布系數(shù),μr,b為后屋面積雪分布系數(shù);非均勻分布雪荷載工況下,μr,m為按照覆蓋保溫被工況下前屋面積雪分布系數(shù)。

    根據(jù)《農業(yè)溫室結構荷載規(guī)范》GB/T51183-2016[30],日光溫室按承載能力極限狀態(tài)設計時,荷載基本組合的效應設計值d按式(2)計算。

    式中d為荷載基本組合的效應設計值,kN/m2;G為永久荷載分項系數(shù),取值為1.0[30];Gk為永久荷載標準值計算的荷載效應值,kN/m2;Q1為主導可變荷載的分項系數(shù),對于雪荷載取值為1.2[30];Q1k為主導可變荷載標準值計算的荷載效應值,kN/m2;Qi為第個非主導可變荷載的分項系數(shù),對于屋面活荷載、作物荷載、移動設備荷載均取1.2[30];ci為第個非主導可變荷載的組合值系數(shù),取值為1.0[30];Qik為第個非主導可變荷載標準值計算的荷載效應值,kN/m2。

    假定本研究日光溫室拱架與基礎和后墻均為固接。對于幾字型鋼日光溫室,考慮縱向系桿對拱架側向位移的限制作用,在有限元模型中,約束縱向系桿與拱架下翼緣連接節(jié)點處面外自由度。

    2.3 模型分析方法

    本研究采用一致缺陷模態(tài)法考慮日光溫室拱架初始幾何缺陷。在ANSYS分析過程中,首先對結構進行特征值屈曲分析,將得到的一階屈曲模態(tài)按照一定比例疊加到理想模型中作為初始幾何缺陷。在彈塑性非線性分析過程中,打開大變形開關以考慮幾何非線性,為完整追蹤荷載-位移曲線極限點后的平衡路徑,即荷載-位移曲線下降段,采用弧長法計算日光溫室拱架整個破壞過程。

    3 結果與分析

    3.1 失穩(wěn)機理

    為了分析開口截面與閉口截面日光溫室拱架在雪荷載作用下失穩(wěn)機理,本研究選取凈截面面積、上翼緣寬度、腹板高度和壁厚均相等的幾字型鋼與空心矩形鋼管作對比,2種截面尺寸如圖4所示。經(jīng)計算,2種截面凈面積均為464 mm2;對于幾字型鋼,其繞軸的慣性矩為3.12×105mm4,繞軸的慣性矩為2.85×105mm4;對于空心矩形鋼管,其繞軸的慣性矩為3.27×105mm4,繞軸的慣性矩為1.94×105mm4。

    圖4 截面尺寸示意圖

    圖5展示了在雪荷載作用下,2種截面日光溫室拱架屋脊節(jié)點荷載-位移曲線。由圖5可知,在均勻和非均勻分布2種雪荷載工況下,拱架荷載-位移曲線變化規(guī)律相似。當荷載較小時,拱架主要發(fā)生繞軸的整體彎曲變形,由于空心矩形鋼管繞軸的慣性矩稍大于幾字型鋼繞軸的慣性矩,因此,在相同荷載作用下,空心矩形鋼管拱架位移低于幾字型鋼拱架。隨著荷載進一步增加,拱架截面出現(xiàn)局部屈曲,并伴隨發(fā)生繞軸的變形,由于幾字型鋼繞軸的慣性矩大于空心矩形鋼管繞軸的慣性矩,致使幾字型鋼拱架極限承載力稍高于空心矩形鋼管拱架極限承載力。在均勻分布雪荷載作用下,空心矩形鋼管和幾字型鋼日光溫室拱架極限承載力分別為4.13和4.29 kN/m2,2種截面日光溫室拱架達到極限承載力時所發(fā)生的位移分別為26.90和39.34 mm;在非均勻分布雪荷載作用下,空心矩形鋼管和幾字型鋼日光溫室拱架極限承載力分別為1.17和1.21 kN/m2,2種截面日光溫室拱架達到極限承載力時所發(fā)生的位移分別為28.17和40.70 mm;與均勻分布雪荷載作用下的極限承載力相比,此2種截面日光溫室拱架在非均勻分布雪荷載作用下的極限承載力僅為均勻分布雪荷載作用下的28%。因此,此類日光溫室拱架對非均勻分布雪荷載更為敏感,在溫室結構設計中,應重點關注非均勻分布雪荷載情況。在接下來的分析過程中,主要討論幾字型鋼日光溫室拱架在非均勻分布雪荷載作用下的失穩(wěn)機理。

    注:A代表結構開始屈服,B代表結構達到極限承載力,C代表結構破壞。

    為進一步了解幾字型鋼日光溫室拱架在非均勻分布雪荷載作用下的破壞過程,圖6展示了該溫室拱架在不同雪荷載水平下的等效應力發(fā)展過程。由圖6a可知,當雪荷載為0.92 kN/m2時,溫室屋脊和后屋面支座處開始屈服,而其他區(qū)域均仍處于彈性狀態(tài);隨著荷載的不斷增加,溫室屋脊和后支座處塑性區(qū)開始由翼緣向腹板擴展,同時,前屋面2/5跨度處(距離屋脊水平投影2.65 m)也開始進入屈服狀態(tài),當荷載達1.21 kN/m2時(圖6b),屋脊和后屋面支座處全截面進入屈服狀態(tài),結構達到極限承載力;當達到極限承載力之后,雖然降低雪荷載可以減小屋脊處等效應力(圖6c),然而,由于前屋面位移持續(xù)增加,致使前屋面2/5跨度處截面應力持續(xù)增加,當前屋面2/5跨度處全截面進入屈服狀態(tài)時,結構最終發(fā)生破壞。

    圖6 非均勻分布雪荷載作用下等效應力

    3.2 縱向系桿對極限承載力的影響

    縱向系桿不僅可以傳遞荷載,還可以為日光溫室拱架提供平面外側向支撐,從而提高日光溫室拱架平面外穩(wěn)定性。本研究日光溫室在前屋面和屋脊處布置4道縱向系桿(圖1)。在有限元分析過程中,通過約束縱向系桿與拱架連接節(jié)點處側向位移來考慮縱向系桿對拱架支撐作用。圖7為非均勻分布雪荷載作用下,幾字型鋼日光溫室屋脊節(jié)點處荷載-位移曲線。由圖7可知,對于幾字型鋼日光溫室拱架,縱向系桿可有效抑制拱架平面外變形,從而提高拱架極限承載力。有縱向系桿約束條件下的極限承載力約是無縱向系桿約束條件下極限承載力的1.25倍。因此,有效布置縱向系桿對提高幾字型鋼日光溫室拱架極限承載力具有重要意義。

    圖7 有無系桿約束下荷載-位移曲線

    3.3 初始幾何缺陷對極限承載力的影響

    日光溫室拱架在制作、運輸和組裝過程中不可避免產生初始幾何缺陷,由于初始幾何缺陷的存在,加劇了荷載的二次效應,進而降低結構極限承載力。為了探究初始幾何缺陷對拱架極限承載力的影響,本研究采用一致缺陷模態(tài)法獲得拱架初始幾何缺陷,即假設拱架初始幾何缺陷分布與其最低階屈曲模態(tài)一致,其初始幾何缺陷如圖8a所示。圖8b為不同初始幾何缺陷幅值下,幾字型鋼日光溫室拱架在非均勻分布雪荷載作用下的荷載-位移曲線。由圖8b可知,該日光溫室拱架極限承載力隨初始幾何缺陷幅值的增加而降低。當初始幾何缺陷幅值從5 mm增加到20 mm時,極限承載力降低了約2%,說明該日光溫室拱架對初始幾何缺陷敏感度較低。

    圖8 初始幾何缺陷對極限承載力的影響

    3.4 截面參數(shù)對極限承載力的影響

    當前,幾字型鋼截面尺寸的選取往往依據(jù)建造者經(jīng)驗,且缺乏相關規(guī)范參考。本節(jié)中,以截面尺寸=50 mm,=70 mm,=15 mm,=12 mm,=2 mm的幾字型鋼為原型,在保證凈截面面積,即用鋼量不變的條件下,通過改變任意2個截面參數(shù),來探究不同截面參數(shù)對幾字型鋼日光溫室拱架極限承載力的影響,進而為幾字型鋼日光溫室拱架截面選取提供參考。

    3.4.1 上翼緣寬度和翻邊寬度

    首先,在保證凈截面面積、腹板高度、下翼緣寬度和壁厚不變的情況下,探究上翼緣寬度和翻邊寬度對極限承載力的影響。經(jīng)計算,當腹板高度=70 mm,下翼緣寬度=15 mm,壁厚=2 mm,為保證凈截面面積不變,上翼緣寬度()與翻邊寬度()需滿足274 mm,當上翼緣寬度從30 mm增加到70 mm時,翻邊寬度則從22 mm降到2 mm,相應上翼緣寬度與翻邊寬度之比則從1.36增加到35。圖9為非均勻分布雪荷載作用下,上翼緣寬度與翻邊寬度之比對極限承載力的影響。由圖9可知,隨著上翼緣寬度與翻邊寬度之比從1.36增加到35,極限承載力呈先增加后降低的變化規(guī)律。當上翼緣寬度與翻邊寬度之比小于4.17時,極限承載力隨上翼緣寬度與翻邊寬度之比的增加而增加;當上翼緣寬度與翻邊寬度之比大于4.17時,極限承載力隨上翼緣寬度與翻邊寬度之比的增加而降低。這是由于當上翼緣寬度增加到一定程度時,上翼緣開始出現(xiàn)局部屈曲,此外,隨著翻邊寬度的減小,翻邊對下翼緣的約束作用減弱,致使極限承載力下降。當上翼緣寬度與翻邊寬度之比從4.17增加到35時,極限承載力下降約13%。因此,建議幾字型鋼上翼緣寬度與翻邊寬度之比宜控制在4.17左右。

    注:此分析為非均勻分布雪荷載工況。下同。

    3.4.2 腹板高度和翻邊寬度

    在保證凈截面面積、上翼緣寬度、下翼緣寬度和壁厚不變的情況下,探究腹板高度和翻邊寬度對極限承載力的影響。經(jīng)計算,當上翼緣寬度=50 mm,下翼緣寬度=15 mm,壁厚=2 mm,為保證凈截面面積不變,腹板高度()與翻邊寬度()需滿足82 mm,當腹板高度從60 mm增加到80 mm時,翻邊寬度則從22 mm降到2 mm,相應腹板高度與翻邊寬度之比則從2.72增加到40。圖10為非均勻分布雪荷載作用下,腹板高度與翻邊寬度之比對極限承載力的影響。由圖10可知,當腹板高度與翻邊寬度之比小于9.25時,極限承載力隨腹板高度與翻邊寬度之比的增加而增加;當腹板高度與翻邊寬度之比從2.72增加到9.25時,極限承載力增加約30%。而當腹板高度與翻邊寬度之比大于9.25時,極限承載力幾乎保持不變,這是由于當腹板增加到一定高度時,腹板開始出現(xiàn)屈曲,結構極限承載力由腹板局部屈曲荷載控制,因此,進一步增加腹板高度并不能提高結構極限承載力。建議幾字型鋼腹板高度與翻邊寬度之比不大于9.25。

    3.4.3 下翼緣寬度和翻邊寬度

    在保證凈截面面積、上翼緣寬度、腹板高度和壁厚不變的情況下,探究下翼緣寬度和翻邊寬度對極限承載力的影響。經(jīng)計算,當上翼緣寬度=50 mm,腹板高度=70 mm,壁厚=2 mm,為保證凈截面面積不變,下翼緣寬度()與翻邊寬度()需滿足27 mm,當下翼緣寬度從5 mm增加到25 mm時,翻邊寬度則從22 mm降到2 mm,相應下翼緣寬度與翻邊寬度之比則從0.23增加到12.5。圖11為非均勻分布雪荷載作用下,下翼緣寬度與翻邊寬度之比對極限承載力的影響。由圖11可知,當下翼緣寬度與翻邊寬度之比小于1.7時,極限承載力隨下翼緣寬度與翻邊寬度之比增加而增加;當下翼緣寬度與翻邊寬度之比從0.23增加到1.7時,極限承載力提高約41%。因此,在不增加凈截面面積、翼緣寬度和腹板高度的情況下,僅通過改變下翼緣寬度與翻邊寬度之比,即可有效提升拱架極限承載力。然而,當下翼緣寬度與翻邊寬度之比大于1.7時,極限承載力幾乎不再隨著下翼緣寬度與翻邊寬度之比增加而增加,原因為:隨著翻邊寬度的減小,翻邊對下翼緣的約束減弱,同時,隨著下翼緣寬度的增加,下翼緣出現(xiàn)局部屈曲現(xiàn)象,此時拱架極限承載力由下翼緣局部屈曲荷載控制。因此,建議幾字型鋼下翼緣寬度與翻邊寬度之比不大于1.7。

    圖10 腹板高度與翻邊寬度之比對極限承載力的影響

    圖11 下翼緣寬度與翻邊寬度之比對極限承載力的影響

    3.4.4 上翼緣寬度和下翼緣寬度

    在保證凈截面面積、腹板高度、翻邊寬度和壁厚不變的情況下,探究上翼緣寬度和下翼緣寬度對極限承載力的影響。經(jīng)計算,當腹板高度=70 mm,翻邊寬度=12 mm,壁厚=2 mm,為保證凈截面面積不變,上翼緣寬度()與下翼緣寬度()需滿足280 mm,當上翼緣寬度從30 mm增加到70 mm時,下翼緣寬度則從25 mm降到5 mm,相應上翼緣寬度與下翼緣寬度之比從1.2增加到14。圖12為非均勻分布雪荷載作用下,上翼緣寬度與下翼緣寬度之比對極限承載力的影響。由圖12可知,隨著上翼緣寬度與下翼緣寬度之比的增加,極限承載力先上升后下降。當上翼緣寬度與下翼緣寬度之比小于3.33時,極限承載力隨上翼緣寬度與下翼緣寬度之比的增加而增加;而當上翼緣寬度與下翼緣寬度之比大于3.33時,極限承載力隨上翼緣寬度與下翼緣寬度之比的增加而急劇下降,這是由于隨著上翼緣寬度的增加,腹板對上翼緣的約束作用減弱,同時,上翼緣寬度的增加致使其發(fā)生局部屈曲現(xiàn)象;此外,隨著下翼緣寬度的減小,下翼緣對腹板的約束作用減弱,致使腹板也出現(xiàn)局部屈曲現(xiàn)象,從而造成極限承載力急劇下降。當上翼緣寬度與下翼緣寬度之比從3.33增加到14時,極限承載力下降約34%。因此,建議幾字型鋼上翼緣寬度與下翼緣寬度之比宜控制在3.33左右。

    圖12 上翼緣寬度與下翼緣寬度之比對極限承載力的影響

    3.4.5 腹板高度和下翼緣寬度

    在保證凈截面面積、上翼緣寬度、翻邊寬度和壁厚不變的情況下,探究腹板高度和下翼緣寬度對極限承載力的影響。經(jīng)計算,當上翼緣寬度=50 mm,翻邊寬度=12 mm,壁厚=2 mm,為保證凈截面面積不變,腹板高度()與下翼緣寬度()需滿足85 mm,當腹板高度從60 mm增加到80 mm時,下翼緣寬度則從25 mm降到5 mm,相應腹板高度與下翼緣寬度之比從2.4增加到16。圖13為非均勻分布雪荷載作用下,腹板高度與下翼緣寬度之比對極限承載力的影響。由圖13可知,隨著腹板高度與下翼緣寬度之比的不斷增加,極限承載力先上升后下降。當腹板高度與下翼緣寬度之比小于4.67時,極限承載力隨腹板高度與下翼緣寬度之比增加而增加;然而,當腹板高度與下翼緣寬度之比大于4.67時,極限承載力隨腹板高度與下翼緣寬度之比增加而下降;當腹板高度與下翼緣寬度之比從4.67增加到16時,極限承載力下降約21%,這是由于隨著下翼緣寬度的減小,下翼緣對腹板的約束作用減弱,同時,由于腹板高度的增加,腹板出現(xiàn)局部屈曲現(xiàn)象,致使極限承載力降低。當腹板高度與下翼緣寬度之比為4.67時,極限承載力達到最大值。因此,建議幾字型鋼腹板高度與下翼緣寬度之比宜控制在4.67左右。

    圖13 腹板高度與下翼緣寬度之比對極限承載力的影響

    4 結論與建議

    本文選取某跨度為8 m,脊高為3.8 m的幾字型鋼日光溫室為研究對象,采用ANSYS有限元軟件,對其在雪荷載作用下的失穩(wěn)破壞過程進行模擬,并探究了縱向系桿、初始幾何缺陷和截面參數(shù)對其極限承載力的影響,得到結論如下:

    1)在凈截面面積、上翼緣寬度、腹板高度和壁厚相同的條件下,幾字型鋼日光溫室拱架極限承載力稍高于空心矩形鋼管日光溫室拱架極限承載力。

    2)相較于均勻分布雪荷載,該日光溫室拱架對非均勻分布雪荷載更為敏感。幾字型鋼日光溫室拱架在非均勻分布雪荷載作用下的極限承載力約是均勻分布雪荷載作用下的28%。因此,在日光溫室結構設計中,應重點關注非均勻分布雪荷載的情況。

    3)縱向系桿的設置可有效抑制幾字型鋼日光溫室拱架平面外變形,從而提高結構極限承載力。在非均勻分布雪荷載作用下,有縱向系桿約束的極限承載力約是無縱向系桿約束條件下的1.25倍。該幾字型鋼日光溫室拱架對初始幾何缺陷敏感度較低。

    4)建議幾字型鋼上翼緣寬度與翻邊寬度之比控制在4.17左右,腹板高度與翻邊寬度之比不大于9.25,下翼緣寬度與翻邊寬度之比不大于1.7,上翼緣寬度與下翼緣寬度之比控制在3.33左右,腹板高度與下翼緣寬度之比控制在4.67左右。

    [1] 魏曉明,周長吉,曹楠,等. 中國日光溫室結構及性能的演變[J]. 江蘇農業(yè)學報,2012,28(4):855-860.

    Wei Xiaoming, Zhou Changji, Cao Nan, et al. Evolution of structure and performance of Chinese solar greenhouse[J]. Jiangsu Journal of Agricultural Sciences, 2012, 28(4): 855-860. (in Chinese with English abstract)

    [2] Wang C, Nan B, Wang T L, et al. Wind pressure acting on greenhouses: A review[J]. International Journal of Agricultural and Biological Engineering, 2021, 14(2): 1-8.

    [3] 張琪靜,劉建斌,周朝輝,等. 保護地設施果樹雪災防范與減災技術[J]. 北方果樹,2022(1):30-33.

    [4] 周長吉. 周博士考察拾零(八十二)桁架結構日光溫室骨架及其構造[J]. 農業(yè)工程技術,2018,38(19):54-61.

    [5] 王斌,金寶宏,宋建夏. 豎向荷載下日光溫室鋼骨架承載力有限元分析[J]. 寧夏大學學報(自然科學版),2009,30(4):336-338.

    Wang Bin, Jin Baohong, Song Jianxia. Finite element analysis of load bearing capacity of steel skeletons for solar greenhouses under vertical loads[J]. Journal of Ningxia University (Natural Science Edition), 2009, 30(4): 336-338. (in Chinese with English abstract)

    [6] 唐中祺,劉麗霞,王瑞東,等. 兩種日光溫室鋼骨架結構安全性能分析[J]. 甘肅農業(yè)大學學報,2016,51(6):53-57.

    Tang Zhongqi, Liu Lixia, Wang Ruidong, et al. Safety performance of two steel frames of solar greenhouse[J]. Journal of Gansu Agricultural University, 2016, 51(6): 53-57. (in Chinese with English abstract)

    [7] 王洪義,祖歌,楊鳳軍,等. 高緯度地區(qū)多功能日光溫室設計[J]. 農業(yè)工程學報,2020,36(6):170-178.

    Wang Hongyi, Zu Ge, Yang Fengjun, et al. Design of multi-functional solar greenhouses in high latitude areas[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(6): 170-178. (in Chinese with English abstract)

    [8] 關瑞,吳堅勇,劉中華,等. 日光溫室鋼架結構分析及優(yōu)化[J]. 山西農業(yè)大學學報(自然科學版),2019,39(6):100-106.

    Guan Rui, Wu Jianyong, Liu Zhonghua, et al. Analysis and optimization of greenhouse steel frame structure[J]. Journal of Shanxi Agricultural University (Natural Science Edition), 2019, 39(6): 100-106. (in Chinese with English abstract)

    [9] 郄麗娟,趙獻輝,許廣義,等. 三連棟日光溫室拱式桁架結構靜力分析與結構優(yōu)化[J]. 中國農學通報,2014,30(35):242-247.

    Qie Lijuan, Zhao Xianhui, Xu Guangyi, et al. Static analysis of truss type steel structure for three-span sunlight greenhouse[J]. Chinese Agricultural Science Bulletin, 2014, 30(35): 242-247. (in Chinese with English abstract)

    [10] 袁余,趙驚鵬,張程,等. 日光溫室結構拓撲優(yōu)化設計[J]. 農機化研究,2013,35(10):55-59.

    Yuan Yu, Zhao Jingpeng, Zhang Cheng, et al. The topological optimization design of the structure of the solar greehouse[J]. Journal of Agricultural Mechanization Research, 2013, 35(10): 55-59. (in Chinese with English abstract)

    [11] 周鳳,何斌. 日光溫室新型裝配式骨架節(jié)點設計與有限元分析[J]. 北方園藝,2021(12):50-56.

    Zhou Feng, He Bin. Design and finite element analysis of new fabricated skeleton joints in solar greenhouse[J]. Northern Horticulture, 2021(12): 50-56. (in Chinese with English abstract)

    [12] 姜迎春,白義奎,王永剛,等. 考慮脈動風速的平面剛架日光溫室結構動力響應規(guī)律[J]. 中國農業(yè)大學學報,2019,24(4):136-147.

    Jiang Yingchun, Bai Yikui, Wang Yonggang, et al. Dynamic response analyses of plane frame solar greenhouse considering fluctuating wind speed[J]. Journal of China Agricultural University, 2019, 24(4): 136-147. (in Chinese with English abstract)

    [13] Wang C, Jiang Y C, Li X Y, et al. Wind-induced vibration response analysis of Chinese solar greenhouses[J]. Computers and Electronics in Agriculture, 2021, 181: 105954.

    [14] 周長吉. 中國日光溫室結構的改良與創(chuàng)新(四):溫室骨架形式的改良與創(chuàng)新[J]. 中國蔬菜,2020(7):1-6.

    [15] 周長吉. 周博士考察拾零(一百零一)單管骨架日光溫室結構及其構造[J]. 農業(yè)工程技術,2020,40(4):49-54.

    [16] 齊飛,閆冬梅,胡林,等. 日光溫室荷載組合方法及應用[J]. 農業(yè)工程學報,2022,38(7):264-271.

    Qi Fei, Yan Dongmei, Hu Lin, et al. Method and application of load combination in solar greenhouses[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(7): 264-271. (in Chinese with English abstract)

    [17] 丁敏,施旭棟,李密密,等. 薄膜承載力及其對日光溫室結構穩(wěn)定性能的影響[J]. 農業(yè)工程學報,2013,29(12):194-202.

    Ding Min, Shi Xudong, Li Mimi, et al. Load-bearing capacity of films and its effect on structure stability of Chinese solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(12): 194-202. (in Chinese with English abstract)

    [18] 丁敏,李密密,施旭棟,等. 考慮覆蓋材料蒙皮效應的溫室結構穩(wěn)定承載力計算[J]. 農業(yè)工程學報,2016,32(增刊1):224-232.

    Ding Min, Li Mimi, Shi Xudong, et al. Stable bearing capacity calculation of greenhouse structures considering skin effect of covering material[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(Supp.1): 224-232. (in Chinese with English abstract)

    [19] 閆冬梅,胡林,周長吉,等. 橢圓管單管拱架日光溫室結構性能分析[J]. 農業(yè)工程學報,2022,38(5):217-224.

    Yan Dongmei, Hu Lin, Zhou Changji, et al. Analysis of the performance of the oval tube arch of single-tube solar greenhouses[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(5): 217-224. (in Chinese with English abstract)

    [20] 齊飛,閆冬梅,魏曉明. 日光溫室前屋面支撐位置對實腹式骨架安全性的影響[J]. 農業(yè)工程學報,2020,36(16):174-181.

    Qi Fei, Yan Dongmei, Wei Xiaoming. Influences of south roof support position change on the skeleton structure safety in solid belly solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(16): 174-181. (in Chinese with English abstract)

    [21] 楊升華,劉雪迎,蔣秀根. 基于臨時加固策略的日光溫室極端風雪災害對策[J]. 江蘇大學學報(自然科學版),2022,43(1):45-53.

    Yang Shenghua, Liu Xueying, Jiang Xiugen. Countermeasures against extreme wind and snow disasters to solar greenhouses based on temporary reinforcement[J]. Journal of Jiangsu University (Natural Science Edition), 2022, 43(1): 45-53. (in Chinese with English abstract)

    [22] Wang C, Jiang Y C, Bai Y K, et al. Numerical study on static properties and failure mechanisms of landing assembled Chinese solar greenhouses[J]. Computers and Electronics in Agriculture, 2021, 188: 106347.

    [23] Wang C, Jiang Y C, Wang T L, et al. Analysis of wind-induced responses of landing assembled Chinese solar greenhouses[J]. Biosystems Engineering, 2022, 220: 214-232.

    [24] Li X Y, Wang C, Jiang Y C, et al. Dynamic response analysis of a whole steel frame solar greenhouse under wind loads[J]. Scientific Reports, 2022, 12(1): 5200

    [25] 白義奎,王鴻,王鐵良,等. 新型鋼骨架結構日光溫室設計與測試[J]. 沈陽農業(yè)大學學報,2013,44(5):542-547.

    Bai Yikui, Wang Hong, Wang Tieliang, et al. Design and test of new steel frame solar greenhouse[J]. Journal of Shenyang Agricultural University, 2013, 44(5): 542-547. (in Chinese with English abstract)

    [26] 張京社,柴文臣. 高寒區(qū)現(xiàn)代日光溫室結構設計與建造[J]. 山西農業(yè)科學,2012,40(8):870-873.

    Zhang Jingshe, Chai Wenchen. Modern solar greenhouse structure design and construction in high latitude and cold climate areas[J]. Journal of Shanxi Agricultural Sciences, 2012, 40(8): 870-873. (in Chinese with English abstract)

    [27] 周長吉. 周博士考察拾零(九十八)以EPS空腔模塊為圍護墻體的鋼骨架裝配式日光溫室[J]. 農業(yè)工程技術,2019,39(31):26-34.

    [28] 周長吉. 周博士考察拾零(一百二十七)一種自防水無拼縫橡塑保溫被圍護的主動儲放熱組裝結構日光溫室——記北京華美沃龍農業(yè)科技有限公司的創(chuàng)新[J]. 農業(yè)工程技術,2022,42(10):36-45.

    [29] 趙云,丁敏,韓盛柏,等. 考慮壓-彎-扭耦合作用的開口截面桿件非線性靜力模型[J]. 中國農業(yè)大學學報,2021,26(3):115-123.

    Zhao Yun, Ding Min, Han Shengbai, et al. Nonlinear static model for the bar with open-section considering pressure-bending-torsion coupling[J]. Journal of China Agricultural University, 2021, 26(3): 115-123. (in Chinese with English abstract)

    [30] 中華人民共和國住房和城鄉(xiāng)建設部,中華人民共和國國家質量監(jiān)督檢驗檢疫總局. 農業(yè)溫室結構荷載規(guī)范:GB/T 51183-2016[S]. 北京:中國計劃出版社,2016.

    Ultimate bearing capacity of the solar greenhouse with hat-shaped steel under snow loads

    Wang Cong1, Jiang Yingchun2, Xu Zhanyang1, Zhang Feng1, Bai Yikui1, Wang Tieliang1※

    (1.110866; 2.110866)

    Hat-shaped steel members are widely used in solar greenhouses, due to their low cost, fast construction, and high material efficiency. This study aims to determine the ultimate bearing capacity of a solar greenhouse with the hat-shaped steel under snow loads. The typical solar greenhouse with an 8 m span and 3.8 m ridge height was selected as the research object. The finite element method (FEM) under ANSYS software was used to analyze the instability mechanism and failure modal of the structure under snow loads (uniform and non-uniform snow loads). An investigation was made to clarify the effects of the longitudinal tie bars, initial geometric imperfections, and sectional dimensions on the ultimate bearing capacity of the structure under non-uniform snow loads. Both the material and geometrical nonlinearity were considered in the finite element model. A bilinear kinematic hardening model was adopted for the steel with a yield strength of 235 MPa, Young’s modulus of 206 GPa, and Poisson’s ratio of 0.3. The geometrical nonlinearity was activated using the ‘NLGEOM’ option. To consider the local buckling, the greenhouse skeletons were then modeled with the Shell181 element suitable for the large strains and rotations. Fixed hinge supports were used for both ends of the skeleton. An arc-length method was utilized to trace the nonlinear load-displacement curve, in order to calculate the ultimate bearing capacity of the structure under snow loads. The ultimate bearing capacity of the solar greenhouse with the hat-shaped steel was slightly higher than that of the hollow rectangular section under the same conditions of net section area, upper flange width, web depth, and wall thickness. The solar greenhouse was more sensitive to the non-uniform snow loads, compared with the uniform ones. The ultimate bearing capacity of the hat-shaped steel solar greenhouse under non-uniform snow loads was about 28% of that under uniform snow loads. Therefore, some suggestions were presented for the non-uniform snow loads in the design stage of the solar greenhouse structure. The roof ridge and north roof end were dangerous sections under non-uniform snow loads, which firstly entered the full section yield state. The longitudinal tie bars were expected to effectively improve the ultimate bearing capacity of the greenhouse structure. The ultimate bearing capacity of the structure with the longitudinal tie bars was about 1.25 times that without tie bars. The ultimate bearing capacity was only reduced by 2%, when the initial geometric imperfections amplitude increased from 5 to 20 mm. It infers that the solar greenhouse was not sensitive to the initial geometric imperfection. The cross-section size of a hat-shaped steel was recommended that the ratio of the upper flange width to the lip width, the upper flange width to the lower flange width, and the web depth to the lower flange width were about 4.17, 3.33, and 4.67, respectively, while, the ratio of the web depth to the lip width, and the lower flange width to lip width were less than 9.25 and 1.7, respectively. These findings can provide a strong reference for the solar greenhouse with the open cold-formed thin-walled steel under snow loads.

    greenhouse; bearing capacity; load; finite element method; buckling failure; cold-formed thin-walled steel; open section; sectional optimization

    10.11975/j.issn.1002-6819.2022.19.019

    S26; S625.1; TU261

    A

    1002-6819(2022)-19-0172-08

    王聰,姜迎春,徐占洋,等. 雪荷載作用下幾字型鋼日光溫室極限承載力分析[J]. 農業(yè)工程學報,2022,38(19):172-179.doi:10.11975/j.issn.1002-6819.2022.19.019 http://www.tcsae.org

    Wang Cong, Jiang Yingchun, Xu Zhanyang, et al. Ultimate bearing capacity of the solar greenhouse with hat-shaped steel under snow loads[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(19): 172-179. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.19.019 http://www.tcsae.org

    2022-08-30

    2022-09-30

    遼寧省重點研發(fā)項目(2021JH2/10200022);遼寧省科學研究經(jīng)費項目(LJKZ0689);國家建設高水平大學公派研究生項目(202108210256);沈陽農業(yè)大學研究生創(chuàng)新培育項目(2022YCXB02)

    王聰,博士生,研究方向為農業(yè)設施與環(huán)境工程。Email:congwang@stu.syau.edu.cn

    王鐵良,教授,博士生導師,研究方向為農業(yè)設施與環(huán)境工程。Email:tieliangwang@syau.edu.cn

    猜你喜歡
    承載力
    再生混凝土抗剪鍵接縫受剪性能及承載力計算
    結構工程師(2022年2期)2022-07-15 02:22:56
    基于單樁豎向承載力計算分析研究
    工程與建設(2019年2期)2019-09-02 01:33:52
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    南通沿海開發(fā)水資源承載力的實證研究
    PVA-ECC抗剪加固帶懸臂RC梁承載力計算研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    基于SAP2000的光伏固定支架結構承載力分析
    太陽能(2015年6期)2015-02-28 17:09:30
    FRP筋混凝土板抗彎承載力理論研究
    K形偏心支撐鋼框架的彈性抗側剛度與極限承載力
    潛艇極限承載力計算與分析
    国产午夜精品论理片| 免费一级毛片在线播放高清视频| 怎么达到女性高潮| 免费在线观看日本一区| 亚洲av五月六月丁香网| 亚洲美女视频黄频| 国产高清有码在线观看视频| 人妻久久中文字幕网| 夜夜躁狠狠躁天天躁| 亚洲av二区三区四区| 亚洲精品日韩av片在线观看 | 欧美成人性av电影在线观看| 中出人妻视频一区二区| 国产精品乱码一区二三区的特点| eeuss影院久久| 国产精品女同一区二区软件 | 男插女下体视频免费在线播放| 熟女少妇亚洲综合色aaa.| 国产亚洲精品一区二区www| 精品人妻偷拍中文字幕| 欧美一区二区国产精品久久精品| 免费在线观看亚洲国产| 无限看片的www在线观看| 久久久久亚洲av毛片大全| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 欧美国产日韩亚洲一区| 变态另类成人亚洲欧美熟女| 国产91精品成人一区二区三区| 色精品久久人妻99蜜桃| 九色国产91popny在线| 欧美最新免费一区二区三区 | 波多野结衣高清作品| 18禁裸乳无遮挡免费网站照片| 日本三级黄在线观看| 少妇高潮的动态图| 日本 av在线| 成人午夜高清在线视频| 日本三级黄在线观看| 国产乱人伦免费视频| 午夜福利欧美成人| 国模一区二区三区四区视频| 精品久久久久久久久久免费视频| 亚洲成人免费电影在线观看| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩卡通动漫| 成熟少妇高潮喷水视频| 亚洲专区国产一区二区| 岛国视频午夜一区免费看| 日韩欧美精品免费久久 | 一级黄色大片毛片| 久久精品国产自在天天线| 天堂√8在线中文| 丰满乱子伦码专区| 国产精品自产拍在线观看55亚洲| 欧美成人a在线观看| 免费在线观看亚洲国产| 少妇的丰满在线观看| 又爽又黄无遮挡网站| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩高清在线视频| av专区在线播放| 国产精品99久久久久久久久| 岛国在线观看网站| 久久久久久人人人人人| 99在线人妻在线中文字幕| 久久久国产成人免费| 国产成人av激情在线播放| 亚洲精品在线美女| 国产av不卡久久| 熟女人妻精品中文字幕| 熟妇人妻久久中文字幕3abv| 丁香欧美五月| 午夜福利在线观看吧| www.熟女人妻精品国产| 不卡一级毛片| 国产精品自产拍在线观看55亚洲| 又爽又黄无遮挡网站| 欧美精品啪啪一区二区三区| 丰满人妻一区二区三区视频av | 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品sss在线观看| 一夜夜www| 成年女人毛片免费观看观看9| 一个人看视频在线观看www免费 | 国产色婷婷99| 亚洲精品一区av在线观看| 99国产综合亚洲精品| 免费看日本二区| 亚洲一区二区三区不卡视频| 草草在线视频免费看| 欧美xxxx黑人xx丫x性爽| 欧美大码av| 精品一区二区三区人妻视频| 听说在线观看完整版免费高清| 欧美黄色淫秽网站| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 麻豆国产av国片精品| 99热精品在线国产| 最好的美女福利视频网| 久久6这里有精品| 在线播放国产精品三级| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产午夜福利久久久久久| 亚洲精品美女久久久久99蜜臀| 亚洲国产色片| 99国产精品一区二区三区| eeuss影院久久| 色精品久久人妻99蜜桃| 日韩大尺度精品在线看网址| 手机成人av网站| 一区二区三区高清视频在线| 18禁黄网站禁片免费观看直播| 热99在线观看视频| 日本精品一区二区三区蜜桃| 九九热线精品视视频播放| 91久久精品国产一区二区成人 | 亚洲欧美日韩高清在线视频| 欧美乱色亚洲激情| 午夜福利高清视频| 国产精品,欧美在线| 精品免费久久久久久久清纯| 精品欧美国产一区二区三| 亚洲精品456在线播放app | 日本撒尿小便嘘嘘汇集6| 99久国产av精品| aaaaa片日本免费| 99久久成人亚洲精品观看| 波多野结衣巨乳人妻| av天堂中文字幕网| 村上凉子中文字幕在线| 欧美午夜高清在线| 性色av乱码一区二区三区2| 丝袜美腿在线中文| 熟女人妻精品中文字幕| xxx96com| 岛国在线观看网站| 亚洲乱码一区二区免费版| 亚洲国产欧洲综合997久久,| 精品免费久久久久久久清纯| 18禁黄网站禁片午夜丰满| av黄色大香蕉| 亚洲欧美日韩东京热| 日韩精品青青久久久久久| 天堂动漫精品| 99riav亚洲国产免费| 亚洲不卡免费看| 亚洲最大成人中文| 国产黄a三级三级三级人| 欧美色欧美亚洲另类二区| 国产精品乱码一区二三区的特点| 美女高潮喷水抽搐中文字幕| 99久久综合精品五月天人人| 亚洲国产高清在线一区二区三| 亚洲精品国产精品久久久不卡| 桃色一区二区三区在线观看| 国产一区二区三区视频了| a级毛片a级免费在线| 蜜桃亚洲精品一区二区三区| av在线蜜桃| 天天添夜夜摸| 亚洲在线自拍视频| 精华霜和精华液先用哪个| 久久久久国内视频| 一区二区三区国产精品乱码| 一a级毛片在线观看| 狠狠狠狠99中文字幕| 欧美中文综合在线视频| 精品国产亚洲在线| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| 国产精品久久久久久亚洲av鲁大| 一本久久中文字幕| 国内精品一区二区在线观看| 麻豆国产97在线/欧美| 国产精品久久久久久久久免 | 中文亚洲av片在线观看爽| 国产一区二区亚洲精品在线观看| 色尼玛亚洲综合影院| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 国产高清videossex| 日韩精品青青久久久久久| 搡老妇女老女人老熟妇| 日韩有码中文字幕| 一本精品99久久精品77| 无限看片的www在线观看| 国内少妇人妻偷人精品xxx网站| 在线免费观看的www视频| 久久精品91无色码中文字幕| 国内精品久久久久久久电影| 亚洲av电影不卡..在线观看| 91在线观看av| 国产爱豆传媒在线观看| 偷拍熟女少妇极品色| 尤物成人国产欧美一区二区三区| 亚洲专区中文字幕在线| 久久亚洲真实| 成人精品一区二区免费| 国产精品,欧美在线| 97超级碰碰碰精品色视频在线观看| 欧美最新免费一区二区三区 | 久久久久久人人人人人| 亚洲精品久久国产高清桃花| 在线国产一区二区在线| 岛国视频午夜一区免费看| 美女高潮的动态| 午夜两性在线视频| 69av精品久久久久久| 怎么达到女性高潮| 色综合亚洲欧美另类图片| 老司机午夜十八禁免费视频| 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 少妇的逼水好多| 丰满的人妻完整版| 又爽又黄无遮挡网站| 999久久久精品免费观看国产| 黄色视频,在线免费观看| 天堂影院成人在线观看| 国产一区二区在线观看日韩 | 欧美区成人在线视频| 最近视频中文字幕2019在线8| 国产成年人精品一区二区| 国产淫片久久久久久久久 | 亚洲精品在线美女| 久久天躁狠狠躁夜夜2o2o| 国产伦精品一区二区三区四那| 国产精品av视频在线免费观看| 热99re8久久精品国产| 亚洲精品一卡2卡三卡4卡5卡| 欧美一区二区国产精品久久精品| 欧美一级a爱片免费观看看| 国产真人三级小视频在线观看| 精品熟女少妇八av免费久了| 在线天堂最新版资源| 亚洲性夜色夜夜综合| 国产综合懂色| 久久婷婷人人爽人人干人人爱| 午夜激情福利司机影院| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片| 欧美乱色亚洲激情| 噜噜噜噜噜久久久久久91| 两个人看的免费小视频| 午夜a级毛片| 亚洲成av人片免费观看| 国产成人av激情在线播放| 黄色成人免费大全| 精品人妻1区二区| 国产精品女同一区二区软件 | 久久中文看片网| 久久久久九九精品影院| 中文字幕人妻熟人妻熟丝袜美 | 欧美xxxx黑人xx丫x性爽| 国产在视频线在精品| 婷婷精品国产亚洲av| 99久久精品热视频| 亚洲五月婷婷丁香| 欧美最黄视频在线播放免费| 中文字幕av成人在线电影| 色精品久久人妻99蜜桃| 国产成+人综合+亚洲专区| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 免费av毛片视频| 久久草成人影院| 久久久久久人人人人人| 1024手机看黄色片| 全区人妻精品视频| 久久天躁狠狠躁夜夜2o2o| 欧美绝顶高潮抽搐喷水| 亚洲性夜色夜夜综合| 人妻丰满熟妇av一区二区三区| 一a级毛片在线观看| 日本撒尿小便嘘嘘汇集6| 国产免费一级a男人的天堂| 欧美黑人巨大hd| 国产精品亚洲美女久久久| 99在线人妻在线中文字幕| 精品人妻一区二区三区麻豆 | 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 无人区码免费观看不卡| 一卡2卡三卡四卡精品乱码亚洲| 免费av观看视频| 国产精品av视频在线免费观看| 久久午夜亚洲精品久久| 欧美高清成人免费视频www| 丁香欧美五月| 十八禁人妻一区二区| 国产激情欧美一区二区| 久久人人精品亚洲av| 欧美在线一区亚洲| 精品一区二区三区av网在线观看| 好男人电影高清在线观看| 免费观看人在逋| xxx96com| 亚洲成人久久性| 麻豆国产97在线/欧美| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 欧美+亚洲+日韩+国产| 观看美女的网站| 国产精品 欧美亚洲| 午夜福利在线观看免费完整高清在 | or卡值多少钱| 精品99又大又爽又粗少妇毛片 | 中出人妻视频一区二区| 亚洲人成伊人成综合网2020| 国产伦一二天堂av在线观看| 我要搜黄色片| 好看av亚洲va欧美ⅴa在| 欧美一区二区国产精品久久精品| 成人欧美大片| 欧美一区二区国产精品久久精品| 成人精品一区二区免费| 亚洲内射少妇av| 最近最新中文字幕大全电影3| 亚洲片人在线观看| 日本在线视频免费播放| 久久精品国产亚洲av香蕉五月| 狂野欧美激情性xxxx| 最好的美女福利视频网| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 亚洲人成伊人成综合网2020| 女人高潮潮喷娇喘18禁视频| 美女高潮的动态| 国产男靠女视频免费网站| 少妇人妻精品综合一区二区 | eeuss影院久久| www国产在线视频色| 此物有八面人人有两片| aaaaa片日本免费| 此物有八面人人有两片| e午夜精品久久久久久久| 欧美+亚洲+日韩+国产| 亚洲人成电影免费在线| 12—13女人毛片做爰片一| tocl精华| 国产高潮美女av| 精品一区二区三区视频在线 | 性欧美人与动物交配| 女人高潮潮喷娇喘18禁视频| 久久久国产成人精品二区| 亚洲精品一区av在线观看| 很黄的视频免费| 男人舔奶头视频| 国产成人av教育| 亚洲专区中文字幕在线| 免费在线观看成人毛片| 欧美日本亚洲视频在线播放| 国产精品久久久久久久久免 | 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3| 亚洲激情在线av| 在线播放国产精品三级| 久久久久久久久久黄片| 免费在线观看日本一区| 欧美一区二区国产精品久久精品| 99久久精品一区二区三区| 国产国拍精品亚洲av在线观看 | 国产私拍福利视频在线观看| 国产成人a区在线观看| 99热只有精品国产| 国产色爽女视频免费观看| 国产91精品成人一区二区三区| 一本精品99久久精品77| 十八禁网站免费在线| 啦啦啦观看免费观看视频高清| 床上黄色一级片| 九色国产91popny在线| 97超视频在线观看视频| 99精品在免费线老司机午夜| 琪琪午夜伦伦电影理论片6080| 波多野结衣巨乳人妻| 国产一区二区激情短视频| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 日日夜夜操网爽| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品综合一区在线观看| 免费看美女性在线毛片视频| 韩国av一区二区三区四区| 天天添夜夜摸| 欧美日韩瑟瑟在线播放| 欧美xxxx黑人xx丫x性爽| 在线观看免费视频日本深夜| 国产精华一区二区三区| 久久久久久九九精品二区国产| 成人性生交大片免费视频hd| 此物有八面人人有两片| 亚洲精品456在线播放app | av中文乱码字幕在线| 日本撒尿小便嘘嘘汇集6| 久久人人精品亚洲av| 少妇人妻精品综合一区二区 | 舔av片在线| 国产精品 欧美亚洲| 日韩免费av在线播放| 国产精品 国内视频| 成人精品一区二区免费| 黄色女人牲交| 老司机午夜十八禁免费视频| 在线天堂最新版资源| 少妇裸体淫交视频免费看高清| 色尼玛亚洲综合影院| 国产精品久久久久久亚洲av鲁大| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 88av欧美| 法律面前人人平等表现在哪些方面| 人人妻,人人澡人人爽秒播| 成人午夜高清在线视频| 一区二区三区国产精品乱码| 欧美性猛交╳xxx乱大交人| 亚洲成av人片在线播放无| 欧美日韩乱码在线| 国产乱人伦免费视频| 一本久久中文字幕| 国产乱人视频| 最新在线观看一区二区三区| 一区福利在线观看| 制服人妻中文乱码| 欧美一区二区国产精品久久精品| 老司机深夜福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 在线观看66精品国产| 精品福利观看| 男人的好看免费观看在线视频| 国产极品精品免费视频能看的| 日韩有码中文字幕| 美女大奶头视频| 日本三级黄在线观看| 日韩欧美三级三区| 看黄色毛片网站| 精品免费久久久久久久清纯| 小说图片视频综合网站| 操出白浆在线播放| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 88av欧美| 欧美日本视频| 嫩草影院入口| 91久久精品电影网| 国产97色在线日韩免费| 国产精品久久久人人做人人爽| 国产成人av激情在线播放| 色综合站精品国产| 亚洲欧美日韩东京热| 在线观看午夜福利视频| 国产一级毛片七仙女欲春2| 国产美女午夜福利| 亚洲精品成人久久久久久| 精品人妻一区二区三区麻豆 | 人妻久久中文字幕网| 香蕉av资源在线| 国产精品电影一区二区三区| 国产精品自产拍在线观看55亚洲| 成年版毛片免费区| 精品一区二区三区视频在线 | svipshipincom国产片| 中文字幕人妻丝袜一区二区| 99久久99久久久精品蜜桃| 18禁美女被吸乳视频| 精品不卡国产一区二区三区| bbb黄色大片| av国产免费在线观看| 国产综合懂色| 中文字幕熟女人妻在线| 婷婷精品国产亚洲av| 99在线视频只有这里精品首页| 欧美日本亚洲视频在线播放| www.熟女人妻精品国产| 欧美日韩中文字幕国产精品一区二区三区| 国内精品一区二区在线观看| 一本一本综合久久| 国产真实伦视频高清在线观看 | 久久精品国产99精品国产亚洲性色| 国产精品久久久人人做人人爽| 99精品在免费线老司机午夜| 51午夜福利影视在线观看| av天堂中文字幕网| 热99re8久久精品国产| 欧美bdsm另类| 免费观看精品视频网站| 性色av乱码一区二区三区2| 国产亚洲精品久久久久久毛片| 无人区码免费观看不卡| 久久精品国产自在天天线| 国产 一区 欧美 日韩| 色综合站精品国产| 免费无遮挡裸体视频| 国产欧美日韩一区二区精品| 91av网一区二区| 日本一本二区三区精品| 亚洲性夜色夜夜综合| av福利片在线观看| 在线免费观看不下载黄p国产 | av天堂在线播放| 日本在线视频免费播放| 成人特级av手机在线观看| 国产毛片a区久久久久| 精品欧美国产一区二区三| 亚洲成人精品中文字幕电影| 中文字幕人妻熟人妻熟丝袜美 | 国产精品一区二区三区四区免费观看 | 日本一本二区三区精品| 欧美一区二区亚洲| 国产精品亚洲一级av第二区| 很黄的视频免费| 国产免费av片在线观看野外av| a级一级毛片免费在线观看| 91在线精品国自产拍蜜月 | 国产成人av教育| eeuss影院久久| 一级毛片女人18水好多| 黄色视频,在线免费观看| 国产精品亚洲一级av第二区| 在线看三级毛片| 国产伦精品一区二区三区视频9 | 国产主播在线观看一区二区| 亚洲av中文字字幕乱码综合| 欧美成人a在线观看| 日本免费一区二区三区高清不卡| 亚洲avbb在线观看| 国产国拍精品亚洲av在线观看 | 久久久久久国产a免费观看| 手机成人av网站| 亚洲成人久久爱视频| 日韩 欧美 亚洲 中文字幕| 欧美最黄视频在线播放免费| 91久久精品国产一区二区成人 | 久久性视频一级片| 成人三级黄色视频| 夜夜躁狠狠躁天天躁| 精品免费久久久久久久清纯| 两个人的视频大全免费| 亚洲avbb在线观看| 午夜福利免费观看在线| av天堂中文字幕网| 两人在一起打扑克的视频| 草草在线视频免费看| 色哟哟哟哟哟哟| x7x7x7水蜜桃| 一进一出好大好爽视频| 香蕉久久夜色| 久久九九热精品免费| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩高清在线视频| 色综合亚洲欧美另类图片| 午夜免费男女啪啪视频观看 | 国产精品亚洲一级av第二区| 最新在线观看一区二区三区| www.999成人在线观看| 又黄又爽又免费观看的视频| 国产精品嫩草影院av在线观看 | 岛国在线观看网站| 国产美女午夜福利| 色综合亚洲欧美另类图片| 久久久久免费精品人妻一区二区| 亚洲成av人片免费观看| 中文字幕人成人乱码亚洲影| 看片在线看免费视频| 男女之事视频高清在线观看| 12—13女人毛片做爰片一| 亚洲不卡免费看| 嫩草影院精品99| 熟女电影av网| 中出人妻视频一区二区| 亚洲国产欧美人成| 黄色成人免费大全| 欧美丝袜亚洲另类 | av国产免费在线观看| 丁香六月欧美| 色精品久久人妻99蜜桃| 久久久久久大精品| 成人一区二区视频在线观看| 成人特级av手机在线观看| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 一个人观看的视频www高清免费观看| 一级黄色大片毛片| 桃红色精品国产亚洲av| 中文字幕av成人在线电影| 国产单亲对白刺激| 久久香蕉国产精品| 午夜两性在线视频| 婷婷丁香在线五月| 淫秽高清视频在线观看| 亚洲成人免费电影在线观看| 欧美成人免费av一区二区三区| 九九热线精品视视频播放| 动漫黄色视频在线观看| 欧美另类亚洲清纯唯美| 国产精品自产拍在线观看55亚洲| 国产99白浆流出| 亚洲不卡免费看| 我要搜黄色片| 亚洲专区中文字幕在线| 久久人妻av系列| 色播亚洲综合网| 村上凉子中文字幕在线| 丁香六月欧美| 国产精华一区二区三区| 一区福利在线观看| 欧美日韩国产亚洲二区| 日韩有码中文字幕| 久久久久久久久久黄片| 国产黄片美女视频| 亚洲一区高清亚洲精品|