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

    平板型鋁合金格柵結(jié)構(gòu)支座節(jié)點的承載性能

    2016-06-17 06:42:10陳偉剛白光波董石麟朱忠義
    關(guān)鍵詞:有限元分析

    陳偉剛,鄧 華, 白光波, 董石麟, 朱忠義

    (1. 浙江大學(xué) 空間結(jié)構(gòu)研究中心,浙江 杭州 310058;2. 浙江省空間結(jié)構(gòu)重點實驗室,浙江 杭州 310058.3. 北京市建筑設(shè)計研究院有限公司,北京 100045)

    ?

    平板型鋁合金格柵結(jié)構(gòu)支座節(jié)點的承載性能

    陳偉剛1,2,鄧華1,2, 白光波1,2, 董石麟1,2, 朱忠義3

    (1. 浙江大學(xué) 空間結(jié)構(gòu)研究中心,浙江 杭州 310058;2. 浙江省空間結(jié)構(gòu)重點實驗室,浙江 杭州 310058.3. 北京市建筑設(shè)計研究院有限公司,北京 100045)

    摘要:通過有限元分析與足尺試驗考察平板型鋁合金格柵結(jié)構(gòu)支座節(jié)點的承載性能.加載試驗揭示了支座節(jié)點的位移、應(yīng)變發(fā)展特點及破壞模式.有限元分析給出了支座節(jié)點中的桿件、環(huán)槽鉚釘及蓋板的應(yīng)力分布狀況.研究表明:支座節(jié)點的破壞為位于工字鋁受拉翼緣的最外排鉚釘孔處發(fā)生斷裂;與受拉翼緣連接的下蓋板存在傳遞路徑明顯的高應(yīng)力區(qū),但由于分布區(qū)域較窄并未導(dǎo)致蓋板破壞,而上蓋板的應(yīng)力要低很多;受拉翼緣斷裂截面處的環(huán)槽鉚釘處于最不利的拉彎受力狀態(tài),但并未發(fā)生破壞;當(dāng)桿件斷裂時,承受的極限彎矩約為純彎狀態(tài)下強(qiáng)度設(shè)計值的0.9倍,剪彎比可達(dá)0.3;節(jié)點的有限元結(jié)果與試驗結(jié)果吻合較好,反映了數(shù)值模擬方法的有效性.

    關(guān)鍵詞:鋁合金結(jié)構(gòu);平板格柵;支座節(jié)點;斷裂破壞;有限元分析

    鋁合金單層網(wǎng)格結(jié)構(gòu)具有重量輕、外形美觀、耐腐蝕性好、易于維護(hù)等優(yōu)點[1],廣泛應(yīng)用于工業(yè)和民用建筑.該類結(jié)構(gòu)的桿件通常采用工字型截面擠壓型材(簡稱“工字鋁”).由于焊接會顯著降低鋁合金母材的強(qiáng)度,工字鋁桿件間通常采用較為可靠的機(jī)械連接形式[2],其中最常用的是板式節(jié)點.板式節(jié)點[3]是由美國Temcor公司研發(fā)的一種節(jié)點體系,具有構(gòu)造簡單,氣動張拉環(huán)槽鉚釘?shù)氖┕に俣瓤斓葍?yōu)點.但由于該節(jié)點僅為各桿件的翼緣板相連而腹板不連,抗剪承載能力較弱,故主要應(yīng)用于球面網(wǎng)殼等以薄膜內(nèi)力為主的單層網(wǎng)格結(jié)構(gòu)[4-5].

    近年來為滿足建筑設(shè)計的要求,國內(nèi)一些重要工程開始將鋁合金單層網(wǎng)殼結(jié)構(gòu)應(yīng)用于平板型的屋蓋格柵結(jié)構(gòu)中(如重慶國際博覽中心、常州西太湖花博會主場館等),并仍采用板式節(jié)點.由于平板型格柵結(jié)構(gòu)以受彎為主,且鋁合金材料的彈性模量較低,因此相比于以薄膜內(nèi)力為主的球面網(wǎng)殼結(jié)構(gòu),平板型格柵結(jié)構(gòu)的跨越能力和抗變形能力均較弱,一般需要設(shè)置相對較密的豎向支承點.出于構(gòu)造簡單、連接便利的考慮,這些支承點處支座節(jié)點的做法主要是遵循機(jī)械連接的原則通過板式節(jié)點改造而成.因支座節(jié)點對結(jié)構(gòu)安全至關(guān)重要,故在進(jìn)行平板型鋁合金格柵結(jié)構(gòu)設(shè)計時,有必要對這些以承受彎矩和剪力為主的支座節(jié)點的受力性能和承載能力進(jìn)行研究.

    本文結(jié)合實際工程,對平板型鋁合金格柵結(jié)構(gòu)的支座節(jié)點進(jìn)行試驗研究和有限元分析.介紹平板型鋁合金格柵結(jié)構(gòu)支座節(jié)點的構(gòu)造及其受力特點,通過試驗考察加載過程中支座節(jié)點的變形、應(yīng)力發(fā)展情況和破壞模式等.采用有限元方法對試件的加載過程進(jìn)行數(shù)值模擬,利用數(shù)值計算結(jié)果對模型試驗的精度和有效性進(jìn)行考察,并對該類型支座節(jié)點的受力特點、薄弱區(qū)域進(jìn)行深入分析,提出此類平板型鋁合金格柵結(jié)構(gòu)支座節(jié)點的分析設(shè)計建議.

    1支座節(jié)點構(gòu)造及受力特點

    1.1節(jié)點構(gòu)造

    某實際工程中平板型鋁合金格柵結(jié)構(gòu)是由2個正交方向桿件和一個斜向桿件構(gòu)成的三向網(wǎng)格體系,如圖1所示.格柵結(jié)構(gòu)由下部撐桿支承,支承點處支座節(jié)點的具體做法是在標(biāo)準(zhǔn)板式節(jié)點的上、下蓋板中心開圓孔,然后穿入一根帶螺紋的實心銷軸,通過擰緊上、下各2個螺母將銷軸與上下蓋板連接起來,在與斜向桿件垂直方向設(shè)置2個加勁撐,桿件、加勁撐與蓋板之間仍采用環(huán)槽鉚釘連接.根據(jù)實際需要,將銷軸下端機(jī)加工成單向鉸耳板或雙向鉸耳板與下部撐桿連接.如圖2所示.

    圖1 平板型鋁合金格柵結(jié)構(gòu)Fig.1 Flat aluminum lattice structure

    圖2 支座節(jié)點構(gòu)造Fig.2 Formation of bearing joint

    1.2受力特點

    與以承受薄膜內(nèi)力為主的單層球面網(wǎng)殼不同,平板型鋁合金格柵結(jié)構(gòu)的節(jié)點在跨中以承受正彎矩為主,在支座處同時承受負(fù)彎矩和較大剪力(支座反力),如圖3所示.

    圖3 豎向荷載下平板型格柵結(jié)構(gòu)的彎矩分布示意圖Fig.3 Illustrative moment distribution of flat aluminum lattice structure under vertical loads

    針對該平板型鋁合金格柵結(jié)構(gòu)實際工程,由設(shè)計單位根據(jù)整體結(jié)構(gòu)分析結(jié)果,提供了代表性支座節(jié)點的6根相連桿件的最不利內(nèi)力設(shè)計值,作為確定節(jié)點試驗加載量的參考值,如表1所示.表中,桿件剪力F與彎矩M間的相對大小通過剪彎比α=FH/M來描述,其中桿件截面高度H=300 mm.從式中可知,α值越大,說明剪力也越大,反之則剪力越小.

    表1 桿件最不利內(nèi)力設(shè)計值

    Tab.1Design values of the most unfavorable internal forces in members

    桿件編號M/(kN·m)F/kNαL151.427.20.16L239.325.90.20L346.724.40.16L456.431.90.17L518.12210.35L63817.40.14

    從表1可以看出,支座節(jié)點處的4根正交桿件(L1~L4)承受較大的彎矩,2根斜向桿件(L5、L6)的彎矩相對較小.節(jié)點域6根桿件的剪力設(shè)計值相差不大,剪彎比最大為0.35,最小為0.14.

    2試驗設(shè)計

    2.1試件和加載裝置

    選取該實際工程中的足尺支座節(jié)點作為試件進(jìn)行靜力加載試驗,如圖4所示.試件中的桿件、加勁撐和圓形蓋板均采用國產(chǎn)6061-T6鋁材,桿件截面規(guī)格和蓋板尺寸如表2所示.銷軸和螺母為Q345-B鋼材,銷軸直徑為100 mm.環(huán)槽鉚釘材質(zhì)為304HC不銹鋼,直徑為9.66 mm,孔徑為9.96 mm.桿件尺寸及編號如圖5所示,F1~F6為桿件L1~L6各加載點對應(yīng)的載荷值,桿件L1~L4的加載點位置相同:桿件L5~L6的加載點位置相同,試件基本參數(shù)均列于表2.

    為滿足該類節(jié)點試件的加載要求,專門設(shè)計了節(jié)點試驗加載架,如圖6所示.該加載架采用Q345鋼,主要由立柱、頂部連梁(包括交叉連梁和邊連梁)、桿端輔助短柱及中心輔助短柱4部分組成.在立柱上距地面700 mm處設(shè)牛腿,牛腿頂面標(biāo)高與桿端輔助短柱頂面標(biāo)高一致.中心輔助短柱為空腔型,腔內(nèi)放置液壓千斤頂.節(jié)點銷軸下端加工成半球狀,與腔內(nèi)千斤頂上設(shè)置的凹槽形墊塊相配合,以保證千斤頂頂緊時節(jié)點可轉(zhuǎn)動.4根立柱柱腳及2根桿端輔助短柱底部通過錨栓固定在試驗室的地槽上.正交桿件L1~L4外端擱置于4根立柱的牛腿上,斜向桿件L5、L6外端放置于2個輔助短柱上.在距牛腿頂面280 mm高處設(shè)置側(cè)向限位角鋼,防止L1~L4在加載過程發(fā)生端截面的扭轉(zhuǎn).利用連梁下方設(shè)置的6個300 kN液壓千斤頂實現(xiàn)加載,千斤頂采用吊架固定在連梁上,各加載點處均放置剛性墊塊,以保證該處桿件能均勻受力.

    圖4 支座節(jié)點試驗試件Fig.4 Experimental specimen

    節(jié)點構(gòu)件規(guī)格/mm數(shù)量材料L1~L4H300×150×8×1046061-T6L5、L6H300×150×6×1026061-T6蓋板Φ450×1026061-T6銷軸M1001Q345-B環(huán)槽鉚釘M9.66160304HC

    圖5 桿件編號及加載點位置Fig.5 Member numbering and loading points

    圖6 加載架和試驗裝置Fig.6 Loading frame and test equipments

    2.2內(nèi)力模擬

    根據(jù)支座節(jié)點的連接構(gòu)造形式及設(shè)計單位提供的支座節(jié)點內(nèi)力值可以看出,該連接近似剛接.同時,為確定試驗加載方案,假定節(jié)點為剛性連接,觀察外力作用下與節(jié)點相連各桿件桿端的內(nèi)力分配情況.如圖7所示為試件的簡化計算模型.

    6根桿件在板式支座節(jié)點中心處定義為剛接連接,桿件外端節(jié)點施加鉸接約束,而中心節(jié)點僅施加一個豎向位移約束.在加載點位置對正交桿件L1~L4施加單位荷載,斜向桿件L5~L6施加0.5倍單位荷載.如表3所示為不同約束條件下的桿件內(nèi)力,表中A、B分別指考慮中心節(jié)點約束和撤去中心節(jié)點約束.從表3中可以看出,當(dāng)施加中心節(jié)點約束時,節(jié)點域各桿件剪力在外荷載作用下的增長速度遠(yuǎn)大于其彎矩值;L1~L4桿件截面的剪彎比α為1.01,L5、L6桿端截面的剪彎比為0.63,和實際工程內(nèi)力相差較大.

    圖7 試件簡化計算模型Fig.7 Simplified analytical model of specimen

    桿件編號AM/(kN·m)F/kNαBM/(kN·m)F/kNαL1~L40.220.741.010.650.080.04L5~L60.100.210.630.500.160.10

    為使桿件截面承受較大的彎矩,撤去中心節(jié)點豎向約束,從表3中可以看出,此時桿件截面彎矩明顯提高,但桿件截面剪力和剪彎比卻又相對較小.經(jīng)綜合分析,最終決定采用兩階段加載法來模擬節(jié)點同時產(chǎn)生較大彎矩和剪力的情況.第1階段加載時先撤去中心節(jié)點的豎向約束,這樣可使桿件在節(jié)點域產(chǎn)生較大的桿端彎矩.當(dāng)彎矩達(dá)到一定量值時,施加中心節(jié)點豎向約束并進(jìn)入第2階段加載.該階段桿件截面彎矩增長有限,而桿件截面剪力增長迅速.

    2.3加載方案

    由我國現(xiàn)行《鋁合金結(jié)構(gòu)設(shè)計規(guī)范》[6]可知6061-T6鋁材的設(shè)計強(qiáng)度為210 MPa,由此可計算試件中工字鋁桿件在純彎狀態(tài)下的強(qiáng)度設(shè)計值為M0=108.78 kN·m.觀察表1,發(fā)現(xiàn)實際工程支座節(jié)點中最大的桿件截面彎矩為56.40 kN·m,約為0.5M0.

    為滿足本工程設(shè)計的要求并有效考察該類節(jié)點的受力性能和極限承載能力,綜合分析后確定試驗中第一加載階段正交桿件上加載點的最大加載值F1~F4取150 kN,斜向桿件上加載點加載值F5~F6取75 kN.由表3可知,此時對應(yīng)的桿件截面彎矩分別為116和97 kN·m.為便于表述,定義此組加載值為F0.

    試驗時采用6個加載點分級同步加載.加載初期每級荷載增量為F0/4.當(dāng)加載值超過3F0/4后,每級荷載增量減少為F0/8.加載至F0(即第1階段加載結(jié)束)后,頂升中心千斤頂開始第2階段加載.此間進(jìn)一步將每級加載值減小至F0/12,直至試件發(fā)生破壞.

    2.4測點布置

    在下蓋板的銷軸邊緣布置2個位移計D1、D2測量試件的豎向位移,如圖8所示.同時,對桿件、蓋板進(jìn)行應(yīng)變測量.應(yīng)變測點位置如圖9所示,分別為:

    圖8 位移測點Fig.8 Displacement gauging points

    1) 蓋板邊緣附近的工字鋁上、下翼緣各布置正應(yīng)變測點2個(所在截面稱為測點截面),主要測量該截面上的測點正應(yīng)力,并計算出桿件截面彎矩.

    2) 桿件測點截面中性軸處的腹板2面對稱布置應(yīng)變花,主要用于測量剪應(yīng)變并計算該截面的剪力.

    3)上、下蓋板表面對應(yīng)各桿件分別布置應(yīng)變測點6個.

    翼緣和腹板上測點編號分別為Sx-i和Wx-i,其中x為桿件編號,i為應(yīng)變片編號,如圖9(b).蓋板測點編號為PU(D)-i,其中U和D分別為上、下蓋板,見圖9(d).

    圖9 應(yīng)變測點Fig.9 Strain gauging points

    3試驗結(jié)果和分析

    3.1試驗現(xiàn)象

    由《鋁合金結(jié)構(gòu)設(shè)計規(guī)范》[6]可知,鋁合金材料的彈性模量為70×103MPa,因此易知T6061-T6型材設(shè)計強(qiáng)度對應(yīng)的應(yīng)變值為3.0×10-3.鄧華等[1]進(jìn)一步對6061-T6鋁材進(jìn)行拉伸試驗,確定該材料σ0.2對應(yīng)的屈服應(yīng)變參考值為3.5×10-3.對該試件進(jìn)行第1階段加載的過程中,鋁合金桿件和節(jié)點域并未出現(xiàn)明顯的變形和破壞現(xiàn)象,且各測點應(yīng)變值均小于屈服應(yīng)變,表明試件總體上處于彈性受力狀態(tài).

    頂升中心輔助短柱空腔內(nèi)的千斤頂直至與銷軸下端半球面貼緊,并開始第2階段加載.當(dāng)加載至約1.2F0時,突然聽到“嘭”的響聲,桿件L1在節(jié)點域內(nèi)的下翼緣沿前排鉚釘孔截面發(fā)生斷裂如圖10(a)所示,而后沿腹板發(fā)生45°左右的斜向斷裂如圖10(b)所示.相應(yīng)地F1加載點的千斤頂壓力迅速下降,且無法繼續(xù)加載,試驗終止.

    圖10 桿件L1斷裂破壞Fig.10 Failure with fracture of member L1

    3.2荷載-位移曲線

    圖11  測點D2荷載-位移曲線Fig.11 Load-deformation curves of gauging point D2

    如圖11所示為位移測點D2的荷載-位移曲線,圖中,Δ為位移;Fn為相對荷載值,F(xiàn)n=當(dāng)前荷載/F0.由于第2階段加載時中心千斤頂已頂升,故該階段測點D2的位移基本不再發(fā)生改變. 從圖11中可以看出,第1階段的荷載與位移呈明顯的線性關(guān)系,也說明試件總體處于彈性受力狀態(tài).

    3.3荷載-應(yīng)變曲線

    1) 蓋板測點

    圖12 下蓋板測點荷載-應(yīng)變曲線Fig.12 Load-strain curves of lower cover plate

    以下蓋板表面應(yīng)變測點的荷載-應(yīng)變ε曲線為對象進(jìn)行分析,如圖12所示.從圖12中可以看出,各測點在外荷載作用下的變化趨勢基本一致.加載前期,各測點均在彈性范圍內(nèi),且應(yīng)變曲線保持線性關(guān)系.然而,加載到約0.7F0時,多個測點開始表現(xiàn)出非線性特征,其中PD-3、PD-4尤為明顯.第1階段加載結(jié)束時,這2個測點的應(yīng)變值實際上已超過屈服應(yīng)變值,并隨荷載的增加進(jìn)一步快速增大.其余測點的應(yīng)變值在整個受力過程中均小于鋁合金的屈服應(yīng)變.相比之下,上蓋板各測點在加載過程中處于受壓狀態(tài),但應(yīng)變最大值僅為1.3×10-3,相比下蓋板要小很多.

    2. 桿件測點

    根據(jù)試件的對稱性,僅對L1、L4、L6這3根桿件翼緣測點應(yīng)變進(jìn)行分析.如圖13所示為3根桿件翼緣測點的荷載-應(yīng)變曲線.從圖13中可以看出,同根桿件上對稱測點的荷載-應(yīng)變曲線基本接近,且上翼緣測點受壓,下翼緣測點受拉.僅測點SL4-1和SL4-2可能受測量誤差的影響存在一定的差異.在第1階段加載過程中,各測點應(yīng)變曲線基本保持線性關(guān)系;進(jìn)入第2加載階段后,各測點的應(yīng)變值雖有所增加,但除測點SL1-3、SL4-4外均小于屈服應(yīng)變值(圖13中虛線對應(yīng)位置).

    圖13 桿件翼緣測點荷載-應(yīng)變曲線Fig.13 Load-strain curves of gauging points on member’s flanges

    3.4桿件測點截面內(nèi)力

    根據(jù)桿件同一測點截面上、下翼緣的4個測點正應(yīng)變值,由材料力學(xué)[7]公式可換算成該截面的彎矩.同理,根據(jù)中性軸處的測點剪應(yīng)變值,也可計算出該截面的剪力.具體計算公式如下:

    (1)

    (2)

    式中:σi為同一測點截面上、下翼緣測點的正應(yīng)力絕對值,i=1,2,3,4;τj為測點截面中性軸處剪應(yīng)力,j=1,2;Iz為工字鋁截面慣性矩;Sz為對工字鋁截面中性軸的靜矩;y為測點距截面中性軸距離;tw為工字鋁截面腹板厚度.

    在各級荷載作用下,桿件L1、L4、L6測點截面彎矩和剪力的變化分別如圖14和15所示.從圖14中可以看出,第1加載階段桿件截面彎矩增長速度較快.進(jìn)行第2階段加載后,彎矩較大桿件L1、L4的測點截面彎矩增速明顯放緩,增量非常有限.桿件L1的荷載-彎矩曲線在最后階段快速減小,這主要是由于該桿發(fā)生斷裂使得截面內(nèi)力迅速釋放.由圖15可知,桿件測點截面的荷載-剪力關(guān)系與彎矩變化正好相反.第1加載階段,剪力只是緩慢增加.當(dāng)中心節(jié)點豎向位移約束后,測點截面剪力隨荷載的增加而快速增加.如表4所示為在試件發(fā)生破壞前一級荷載作用下(該級加載下桿件截面彎矩達(dá)到最大),各桿測點截面的彎矩、剪力及由試驗值計算得到的剪彎比α,表中β1、β2分別為彎矩和剪力試驗值與表1內(nèi)力設(shè)計值的比值.δ1、δ2為誤差,分別為彎矩和剪力的試驗值減去有限元值除于試驗值.從表4中可以看出,發(fā)生斷裂的桿件L1能承受的最大彎矩約為表1中設(shè)計值的1.86倍,純彎狀態(tài)下強(qiáng)度設(shè)計值的0.9倍.桿件承受的最大剪力約為設(shè)計值的3.56倍,剪彎比為0.30.

    4有限元分析模型

    4.1試件材料參數(shù)及本構(gòu)關(guān)系

    支座節(jié)點的數(shù)值模擬涉及到鋁合金型材和蓋板、不銹鋼鉚釘3種材料.在一般情況下,不銹鋼和鋁合金材料的應(yīng)力-應(yīng)變(σ-ε)呈非線性關(guān)系,且2種材料的應(yīng)力-應(yīng)變關(guān)系均可采用Ramberg-Osgood模型[8]來描述,即

    (3)

    圖14 測點截面的荷載-彎矩曲線Fig. 14 Load-bending moment curves of measuring sections

    圖15 測點截面的荷載-剪力曲線Fig. 15 Load-shearing force curves of measuring sections

    式中:n為反映材料應(yīng)力-應(yīng)變關(guān)系非線性程度的系數(shù),稱為應(yīng)變硬化系數(shù);σ0.2為殘余應(yīng)變?yōu)?.2%時對應(yīng)的應(yīng)力;E為初始彈性模量.

    對于鋁合金材料,一般根據(jù)Steinhardt的建議[9]取n=σ0.2/10,其中σ0.2單位為MPa.對于不銹鋼材料,通常采用修正的Ramberg-Osgood模型[10]:

    (4)

    式中:n=ln20/ln(σ0.2/σ0.01);m=1+3.5σ0.2/σu,為應(yīng)力大于σ0.2時的材料應(yīng)變硬化系數(shù);σu為極限應(yīng)力;εu=1-σ0.2/σu為對應(yīng)σu的擬合應(yīng)變; E0.2=E/(1+0.002nE/σ0.2),為殘余應(yīng)變?yōu)?.2%時對應(yīng)的彈性模量,σ0.01為殘余應(yīng)變?yōu)?.01%時對應(yīng)的應(yīng)力.鄧華等[1]根據(jù)材性測試結(jié)果,分別計算出6061-T6鋁合金和304HC不銹鋼2種材料的主要物理和力學(xué)參數(shù),如表5所示.表中ν為泊松比.

    表52種材料的力學(xué)參數(shù)和應(yīng)變硬化系數(shù)

    Tab.5Mechanical parameters and strain hardening coefficients of two materials

    牌號σ0.01/MPaσ0.2/MPaσu/MPaE/GPanmν6061-T6—246.6283.469.824.7—0.3304HC290.1460.7724.21896.483.230.3

    4.2模型建立

    根據(jù)文獻(xiàn)[11]中所確立的鋁合金板件環(huán)槽鉚釘連接有限元模型建立方法,在ABAQUS中采用八節(jié)點六面體非協(xié)調(diào)模式單元(C3D8I)模擬節(jié)點試件的工字鋁桿件、蓋板、加勁撐及環(huán)槽鉚釘如圖16所示.與鋁合金板件相比,中心鋼銷軸具有較大的剛度.為便于分析,建立有限元模型時不再對鋼銷軸按實體元建模,而是將上、下蓋板上銷軸孔邊周圈節(jié)點的三向線位移一并耦合,即模擬上、下蓋板在剛性銷軸約束下的協(xié)同變形情況.6根桿件的加載點位置與試驗相同,采用和試驗相同的方式進(jìn)行加載.

    圖16 試件有限元模型Fig.16 Finite element model of specimen

    在支座節(jié)點有限元模型中,通過建立8個接觸對來考慮環(huán)槽鉚釘與鋁合金板件(蓋板、桿件翼緣)間、蓋板與桿件翼緣間的接觸摩擦以及鉚釘桿和孔壁之間的擠壓作用,接觸對定義及其示意分別如表6和圖17所示.根據(jù)接觸界面介質(zhì)不同將各接觸對的摩擦系數(shù)分為2組,即鋁合金材料之間的摩擦系數(shù)μ1和不銹鋼與鋁合金材料間的摩擦系數(shù)μ2. 在有限元模型中分別取2組摩擦系數(shù)為μ1=0.15、μ2=0.25.在ABAQUS中利用BoltLoad命令將預(yù)緊力分為多個分析步逐步施加.根據(jù)鄧華等[1]的實測結(jié)果,取環(huán)槽鉚釘?shù)念A(yù)緊力F=18 kN.

    表6  接觸對定義

    圖17 板式節(jié)點接觸面Fig.17 Contact interfaces of lap connection

    5有限元結(jié)果的對比與分析

    5.1桿件測點截面內(nèi)力

    為更直觀對比有限元計算結(jié)果與試驗結(jié)果的異同,將有限元分析得到的各桿件截面的荷載-彎矩曲線、荷載-剪力曲線一并列于圖14和圖15.從這2個對比圖可以看出,有限元計算結(jié)果與試驗結(jié)果總體吻合較好.對具體數(shù)值進(jìn)行對比(詳見表4)發(fā)現(xiàn),除桿件L4截面剪力存在18.44%的較大誤差外,其余桿件截面內(nèi)力的誤差均在10%左右.但應(yīng)該指出,以上誤差的產(chǎn)生除試驗中試件加工誤差、構(gòu)件安裝誤差、加載偏心、測量誤差等因素外,實際上由于有限元分析中也存在接觸、摩擦、塑性變形等較為復(fù)雜的非線性因素,在迭代和收斂過程中也存在一定的數(shù)值計算誤差.

    5.2變形及應(yīng)變比較

    1)試件變形

    將通過有限元分析得到的測點D2的荷載-位移曲線一并列于圖11,發(fā)現(xiàn)其與試驗曲線吻合很好.但有限元計算得到的荷載-位移曲線在加載初期有一段斜率較小的波動(交叉)段,這主要是由于蓋板與桿件上、下翼緣接觸面產(chǎn)生相對滑移所引起的.

    2)測點應(yīng)變

    重點關(guān)注受力較大的下蓋板和工字翼緣,將有限元分析得到的部分測點的荷載-應(yīng)變曲線與試驗結(jié)果如圖18所示.從圖18中可以看出,試驗結(jié)果和有限元計算結(jié)果均能較好地吻合,部分測點的應(yīng)變數(shù)值雖在一定程度上受測量誤差的影響,但差異總體上較小.

    5.3破壞形態(tài)分析

    重點關(guān)注桿件L1在桿端下翼緣與蓋板連接區(qū)域的Mises應(yīng)力情況,如圖19所示.從圖19中可以看出,其最大Mises應(yīng)力為310.47 MPa,約為鋁合金材料屈服強(qiáng)度246 MPa的1.3倍,且主要集中在前排鉚釘孔附近的位置,并向后排逐漸遞減,這與加載過程中該位置率先發(fā)生破壞(圖10)是一致的.

    5.4環(huán)槽鉚釘受力

    如圖20所示為連接桿件L1下翼緣與蓋板的環(huán)槽鉚釘Mises應(yīng)力云圖.從圖20中可以看出,同圖19中該位置的鉚釘孔受力情況相對應(yīng),位于前排的環(huán)槽鉚釘受力最大,并向后排逐漸遞減,其最大Mises應(yīng)力為496.28 MPa,略大于304HC不銹鋼材料屈服強(qiáng)度460 MPa.此外,從圖20中還可以看出,受節(jié)點試件承受彎矩和剪力的影響,環(huán)槽鉚釘處于拉彎受力狀態(tài).

    5.5蓋板應(yīng)力

    2塊蓋板在加載完成后的Mises應(yīng)力云圖如圖21所示.從圖21(a)中可以看出,下蓋板局部區(qū)域達(dá)到屈服應(yīng)力,主要集中在斜桿L5、L6相對應(yīng)的銷軸孔邊及L1~L4鉚釘孔連線截面處(圖中黑色箭頭線),且這些應(yīng)力較大區(qū)域也形成明顯的應(yīng)力傳遞路徑,將內(nèi)力經(jīng)銷軸傳至下部支撐.除以上區(qū)域外,下蓋板其他區(qū)域應(yīng)力均小于設(shè)計強(qiáng)度210 MPa,處于彈性范圍內(nèi).上蓋板圖21(b)除斜桿對應(yīng)的銷軸孔區(qū)域及鉚釘孔存在一定的應(yīng)力集中外,整體的應(yīng)力水平與下蓋板相比要低較多,均在彈性范圍內(nèi).

    圖18 典型測點的荷載-應(yīng)變曲線對比Fig. 18 Load-strain curves of some typical measure points

    圖19 L1鉚釘連接處Mises應(yīng)力云圖Fig.19 Mises stress distribution in area of lockbolts connection of member L1

    圖20 環(huán)槽鉚釘Mises應(yīng)力云圖Fig.20 Mises stress distribution of lockbolts

    圖21 蓋板Mises 應(yīng)力云圖Fig.21 Mises stresses distribution of cover plates

    6結(jié)論

    (1) 在試驗技術(shù)方面,兩階段加載方法可以有效模擬以承受彎矩和剪力為主的板式支座節(jié)點的受力狀態(tài).

    (2) 在彎矩和剪力共同作用下,板式支座節(jié)點的破壞形式為節(jié)點域內(nèi)工字鋁桿件受拉翼緣沿鉚釘孔截面的脆性斷裂破壞,斷裂前沒有明顯征兆.

    (3) 上下蓋板的受力存在較大的差異.與工字鋁受拉翼緣連接的下蓋板應(yīng)力較高,應(yīng)力傳遞路徑明顯,但分布區(qū)域較窄,并未導(dǎo)致蓋板發(fā)生破壞.上蓋板整體的應(yīng)力水平與下蓋板相比要低很多

    (4) 受拉翼緣斷裂截面處的環(huán)槽鉚釘受力最為不利,呈拉彎受力狀態(tài),但不銹鋼的強(qiáng)度較高,也未出現(xiàn)鉚釘破壞.

    (5) 斷裂桿件L1承受的最大彎矩約為工程彎矩設(shè)計值的1.86倍,截面純彎狀態(tài)下強(qiáng)度設(shè)計值的0.9倍.桿件能夠承擔(dān)的最大剪力約為工程剪力設(shè)計值的3.56倍.在極限狀態(tài)下,桿件的剪彎比可達(dá)0.30.

    (6)板式支座節(jié)點的有限元分析結(jié)果與試驗結(jié)果吻合較好,反映了本文的有限元模擬方法有效性.

    參考文獻(xiàn)(References):

    [1] 鄧華,陳偉剛,白光波,等. 鋁合金板件環(huán)槽鉚釘搭接連接受剪性能試驗研究[J]. 建筑結(jié)構(gòu)學(xué)報.2016(01): 143-149.

    DENG Hua, CHEN Wei-gang, BAI Guang-bo, et al. Experimental study on shearing behavior of lockbolted lap connection for aluminum alloy plates [J].Journal of Building Structures. 2016(01): 143-149.

    [2] 楊聯(lián)萍,韋申,張其林. 鋁合金空間網(wǎng)格結(jié)構(gòu)研究現(xiàn)狀及關(guān)鍵問題[J]. 建筑結(jié)構(gòu)學(xué)報. 2013, 34(2): 1-19, 60

    YANG Lian-ping, WEI Shen, ZHANG Qi-ling. Aluminum reticulated spatial structures: state of the art and key issues [J]. Journal of Building Structures, 2013, 34(2): 1-19, 60.

    [3] HWANG K J. Advanced investigations of grid spatial structures considering various connection systems [D]. Stuttgart: University of Stuttgart, 2010: 11-45.

    [4] 沈祖炎,郭小農(nóng),李元齊. 鋁合金結(jié)構(gòu)研究現(xiàn)狀簡述[J]. 建筑結(jié)構(gòu)學(xué)報, 2007, 28(6): 100-109.

    SHEN Zu-yan, GUO Xiao-nong, LI Yuan-qi. State of the arts of research on aluminum alloy structures [J]. Journal of Building Structures, 2007, 28(6): 100-109.

    [5] GUO Xiao-nong, XIONG Zhe, LUO Yong-feng, et al. Experimental investigation on the semi-rigid behaviour of aluminium alloy gusset joints [J]. Thin-Walled Structures. 2015, 87 (2): 30-40.

    [6] GB50429—2007. 鋁合金結(jié)構(gòu)設(shè)計規(guī)范[S]. 北京: 中國計劃出版社, 2007:14-17.

    GB50429—2007. Code for design of aluminum structures[S]. Beijing: China Planning Press, 2007: 14-17.

    [7] 孫訓(xùn)方,方孝淑,關(guān)來泰. 材料力學(xué)[M]. 北京: 高等教育出版社, 2001:90-124.

    [8] RAMBERG W, OSGOOD W R. Description of stress-strain curves by three parameters [R]. Washington DC: National advisory committee for aeronautics, 1943.

    [9] MAZZOLANI F M. Aluminum alloy structures [M]. 2th ed. Boca Raton : The Chemical Rubber Company Press, 1994:35-37.

    [10] RASMUSSEN K J R. Full-range stress-strain curves for stainless steel alloys [J]. Journal of Constructional Steel Research. 2003, 59(1): 47-61.

    [11] 陳偉剛. 平板型鋁合金格柵結(jié)構(gòu)板式節(jié)點的受力性能研究[D]. 杭州: 浙江大學(xué), 2015: 34-56.

    CHEN Wei-gang. Investigation on mechanical behavior of gusset joints of flat aluminum alloy lattice structures [D]. Hangzhou: Zhejiang University, 2015: 34-56.

    Load-bearing behavior of bearing joint of flat aluminum alloy lattice structures

    CHEN Wei-gang1,2, DENG Hua1,2, BAI Guang-bo1,2, DONG Shi-lin1,2, ZHU Zhong-yi3

    (1.SpaceStructuresResearchCenter,ZhejiangUniversity,Hangzhou310058,China;2.ZhejiangProvincialKeyLaboratoryofSpaceStructures,ZhejiangUniversity,Hangzhou310058,China;3.BeijingInstituteofArchitecturalDesign(Group)CompanyLimited,Beijing100045,China)

    Abstract:The load-bearing behaviors of bearing joint of flat aluminum alloy lattice structures were investigated by means of finite element analysis and full-scale experiment. The loading test revealed the characteristics of displacement and strain developments of the bearing joint as well as its damage modes. By adopting the finite element analysis, the stress distributions respectively in members, lockbolts and cover plates of the bearing joint were presented. Results show that the outermost row of bolt holes on the tensile flange of aluminum alloy I-beams ruptures and causes the damage of the bearing joint. An obvious high-stress zone exists on the lower cover plate connected with the tensile flange though it is too narrow to induce the fracture of cover plates, while stresses on the upper cover plate are much smaller. The lockbolts in the fracture section of the tensile flange lie in the most unfavorable tensile-bending state, but no damage happens. When the fracture appears in the I-beam, the corresponding ultimate moment is approximately 0.9 times the design strength of aluminum I-beams under pure bending, and the moment-shear ratio reaches 0.3. The results of numerical analysis well coincide with those of the loading test, indicating the validity of the suggested numerical simulation method.

    Key words:aluminum alloy structure; flat lattice structure; bearing joint; fracture failure; finite element analysis

    收稿日期:2015-05-20.浙江大學(xué)學(xué)報(工學(xué)版)網(wǎng)址: www.journals.zju.edu.cn/eng

    作者簡介:陳偉剛(1986-),男,博士生,從事空間結(jié)構(gòu)等研究.ORCID: 0000-0001-5025-5421. E-mail: wgchen@zju.edu.cn通信聯(lián)系人:鄧華(1971-),男,教授,博導(dǎo).ORCID: 0000-0002-0792-0518. E-mail: denghua@zju.edu.cn

    DOI:10.3785/j.issn.1008-973X.2016.05.004

    中圖分類號:TU 395

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1008-973X(2016)05-0831-10

    猜你喜歡
    有限元分析
    對于Pro/mechanica的數(shù)控車刀的有限元分析
    S型壓力傳感器形變的有限元分析
    橫向穩(wěn)定桿的側(cè)傾角剛度仿真與試驗研究
    多自由度升降翻轉(zhuǎn)舞臺的運(yùn)動分析
    演藝科技(2016年11期)2016-12-24 20:18:29
    高精度細(xì)深孔的鏜孔加工工藝分析
    基于LS—DYNA的某汽車乘客安全氣囊支架焊接強(qiáng)度的分析
    基于有限元分析談框剪結(jié)構(gòu)力學(xué)性能的影響因素
    自錨式懸索橋鋼箱梁頂推施工階段結(jié)構(gòu)分析
    隨機(jī)振動載荷下發(fā)射裝置尾罩疲勞壽命分析
    航空兵器(2016年4期)2016-11-28 21:54:01
    有限元分析帶溝槽平封頭的應(yīng)力集中
    国产精品永久免费网站| 色在线成人网| 国产精品免费一区二区三区在线| 91久久精品国产一区二区成人| 在线看三级毛片| 一本精品99久久精品77| 国内精品一区二区在线观看| 五月伊人婷婷丁香| 一夜夜www| 在线看三级毛片| 看片在线看免费视频| 精品一区二区三区视频在线观看免费| 日韩欧美三级三区| 午夜免费激情av| 人妻久久中文字幕网| 国产亚洲精品av在线| 国产爱豆传媒在线观看| a级毛片a级免费在线| 春色校园在线视频观看| 精品久久久久久久久久免费视频| 免费av不卡在线播放| 日韩中字成人| 国产黄色小视频在线观看| 国内精品久久久久久久电影| 亚洲成av人片在线播放无| 久久久久久九九精品二区国产| 久久久久久久久久黄片| 国产淫片久久久久久久久| 久久精品国产鲁丝片午夜精品 | 亚洲熟妇熟女久久| 久久久色成人| 国产av麻豆久久久久久久| 国产老妇女一区| 国产高清激情床上av| 欧美激情在线99| 久久久国产成人精品二区| 国产亚洲精品综合一区在线观看| 亚洲成人精品中文字幕电影| 天美传媒精品一区二区| 亚洲欧美激情综合另类| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 亚洲精品色激情综合| www.www免费av| 直男gayav资源| 又黄又爽又免费观看的视频| 一进一出好大好爽视频| 亚洲av免费高清在线观看| 亚洲人成网站高清观看| 一进一出抽搐gif免费好疼| 免费不卡的大黄色大毛片视频在线观看 | 国产成人福利小说| 久久久久性生活片| 天天一区二区日本电影三级| 免费av毛片视频| 男女视频在线观看网站免费| 国产色爽女视频免费观看| 久久香蕉精品热| 国产免费av片在线观看野外av| 国产精品久久久久久av不卡| 尤物成人国产欧美一区二区三区| 不卡一级毛片| a级毛片a级免费在线| 中文字幕精品亚洲无线码一区| 一进一出抽搐gif免费好疼| avwww免费| 一区二区三区四区激情视频 | 久久久色成人| 18+在线观看网站| 亚洲成人久久性| 久久人人爽人人爽人人片va| 99热6这里只有精品| 国产在线男女| 美女免费视频网站| 啪啪无遮挡十八禁网站| 成人美女网站在线观看视频| 国产精品久久久久久亚洲av鲁大| 免费看美女性在线毛片视频| 中文字幕免费在线视频6| 老熟妇乱子伦视频在线观看| 久久婷婷人人爽人人干人人爱| 最近最新免费中文字幕在线| 亚洲欧美日韩高清在线视频| 99久久久亚洲精品蜜臀av| 哪里可以看免费的av片| 97超级碰碰碰精品色视频在线观看| 成人欧美大片| 久久国产乱子免费精品| 岛国在线免费视频观看| 亚洲天堂国产精品一区在线| 亚洲精品亚洲一区二区| 欧美激情在线99| 亚洲第一电影网av| 欧美一区二区精品小视频在线| 久久精品综合一区二区三区| 中出人妻视频一区二区| 日韩中文字幕欧美一区二区| 成人特级av手机在线观看| 欧美色欧美亚洲另类二区| 韩国av一区二区三区四区| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 欧美国产日韩亚洲一区| 一级黄色大片毛片| 久久久久久国产a免费观看| 国产中年淑女户外野战色| 国产精品人妻久久久久久| 人人妻,人人澡人人爽秒播| 18禁在线播放成人免费| 免费黄网站久久成人精品| 波野结衣二区三区在线| 国产一区二区在线av高清观看| 三级毛片av免费| 亚洲国产日韩欧美精品在线观看| 久久亚洲精品不卡| 给我免费播放毛片高清在线观看| 亚洲成av人片在线播放无| 成人三级黄色视频| 亚洲欧美清纯卡通| 尾随美女入室| 在线观看一区二区三区| 国产淫片久久久久久久久| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 天天躁日日操中文字幕| 99久久精品一区二区三区| 精品人妻1区二区| 搡老妇女老女人老熟妇| 十八禁国产超污无遮挡网站| 色在线成人网| 午夜老司机福利剧场| 国产成人a区在线观看| 一进一出抽搐gif免费好疼| 亚洲狠狠婷婷综合久久图片| 91狼人影院| 五月伊人婷婷丁香| 99在线视频只有这里精品首页| 赤兔流量卡办理| 国产精品无大码| 久久亚洲真实| 啦啦啦韩国在线观看视频| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式| 亚洲中文字幕日韩| 噜噜噜噜噜久久久久久91| 亚洲精品456在线播放app | 日韩 亚洲 欧美在线| 亚洲性久久影院| 婷婷六月久久综合丁香| 真人一进一出gif抽搐免费| 国产在线男女| 黄色女人牲交| 三级毛片av免费| 人妻制服诱惑在线中文字幕| 国产精品一及| 最新在线观看一区二区三区| av天堂中文字幕网| 亚洲最大成人av| 亚洲欧美精品综合久久99| 三级男女做爰猛烈吃奶摸视频| 国产亚洲91精品色在线| 内地一区二区视频在线| 欧美+亚洲+日韩+国产| 嫩草影院新地址| 99久久成人亚洲精品观看| 成人毛片a级毛片在线播放| av在线天堂中文字幕| 精品久久久久久成人av| 精品一区二区三区视频在线观看免费| 91在线精品国自产拍蜜月| 午夜久久久久精精品| 日韩一本色道免费dvd| 国产黄色小视频在线观看| 级片在线观看| 亚洲国产精品sss在线观看| 欧美成人a在线观看| 国内精品久久久久精免费| 最好的美女福利视频网| 99在线视频只有这里精品首页| 欧美+日韩+精品| 久久久久久久久大av| 美女 人体艺术 gogo| 日韩,欧美,国产一区二区三区 | 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 观看美女的网站| 中文在线观看免费www的网站| 国产乱人视频| 日韩欧美精品免费久久| 韩国av一区二区三区四区| 国产免费av片在线观看野外av| 国产成人aa在线观看| 欧美日本视频| 国产亚洲91精品色在线| 欧美国产日韩亚洲一区| 亚洲欧美日韩无卡精品| 好男人在线观看高清免费视频| 男人舔奶头视频| 一进一出抽搐gif免费好疼| av中文乱码字幕在线| 亚洲av.av天堂| 久久精品人妻少妇| 色噜噜av男人的天堂激情| 我要看日韩黄色一级片| 最新在线观看一区二区三区| 中国美白少妇内射xxxbb| 国产午夜福利久久久久久| 91麻豆精品激情在线观看国产| 在线观看午夜福利视频| 亚洲欧美精品综合久久99| 看片在线看免费视频| 久久精品国产亚洲av涩爱 | 亚洲午夜理论影院| 日本成人三级电影网站| 婷婷精品国产亚洲av在线| 国产一区二区在线观看日韩| 亚洲最大成人av| 成人二区视频| 久久久久久九九精品二区国产| 黄色视频,在线免费观看| 日本在线视频免费播放| 免费看a级黄色片| 午夜老司机福利剧场| 少妇的逼水好多| 精品人妻视频免费看| 日本一二三区视频观看| 白带黄色成豆腐渣| 欧美性猛交黑人性爽| 别揉我奶头 嗯啊视频| 午夜福利在线在线| 久久久国产成人免费| 亚洲成av人片在线播放无| 久久久久久久久大av| 午夜激情欧美在线| 999久久久精品免费观看国产| 欧美性感艳星| 美女被艹到高潮喷水动态| 久久精品久久久久久噜噜老黄 | 最近中文字幕高清免费大全6 | 久久99热这里只有精品18| 欧美xxxx性猛交bbbb| 黄色日韩在线| 校园人妻丝袜中文字幕| 久久久久久久久中文| 久久久成人免费电影| 国内精品久久久久久久电影| 两个人的视频大全免费| 色精品久久人妻99蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 亚洲第一电影网av| 精华霜和精华液先用哪个| 亚洲国产色片| 欧美三级亚洲精品| 午夜福利在线观看吧| 在线播放无遮挡| 国产精品,欧美在线| 免费看av在线观看网站| 99国产极品粉嫩在线观看| 国产亚洲精品久久久com| 国产毛片a区久久久久| 色哟哟·www| 国内久久婷婷六月综合欲色啪| 波多野结衣高清无吗| 午夜福利视频1000在线观看| 欧美日韩瑟瑟在线播放| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 亚洲人与动物交配视频| 我的老师免费观看完整版| 欧美在线一区亚洲| 老司机深夜福利视频在线观看| 久久精品影院6| 久久精品国产鲁丝片午夜精品 | 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 美女 人体艺术 gogo| 熟女电影av网| 亚洲精品影视一区二区三区av| 午夜日韩欧美国产| 亚洲内射少妇av| 久久精品国产亚洲av天美| 午夜影院日韩av| 午夜精品久久久久久毛片777| 免费av不卡在线播放| 日韩欧美免费精品| 国产精品人妻久久久久久| 五月玫瑰六月丁香| 国产aⅴ精品一区二区三区波| 成人三级黄色视频| 99热这里只有是精品在线观看| av在线蜜桃| 中亚洲国语对白在线视频| av在线亚洲专区| 亚洲精品粉嫩美女一区| 干丝袜人妻中文字幕| 国产成人一区二区在线| 亚洲成a人片在线一区二区| 色尼玛亚洲综合影院| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 亚洲天堂国产精品一区在线| 香蕉av资源在线| 欧美极品一区二区三区四区| 成人三级黄色视频| 蜜桃亚洲精品一区二区三区| 日本色播在线视频| 欧美丝袜亚洲另类 | 午夜亚洲福利在线播放| 亚洲中文日韩欧美视频| 成年女人永久免费观看视频| 成人一区二区视频在线观看| 久久人人爽人人爽人人片va| 精品一区二区三区av网在线观看| 免费不卡的大黄色大毛片视频在线观看 | 欧美日韩乱码在线| 国产精品99久久久久久久久| 午夜亚洲福利在线播放| 伊人久久精品亚洲午夜| 一个人免费在线观看电影| 亚洲中文日韩欧美视频| 国产精品爽爽va在线观看网站| 啦啦啦观看免费观看视频高清| 国产免费av片在线观看野外av| 特大巨黑吊av在线直播| 搡老岳熟女国产| 欧美3d第一页| 亚洲成人精品中文字幕电影| 成人鲁丝片一二三区免费| 亚洲 国产 在线| 成人av在线播放网站| 色哟哟哟哟哟哟| 99riav亚洲国产免费| 黄色日韩在线| 国产成人影院久久av| 变态另类成人亚洲欧美熟女| 午夜福利在线观看吧| 亚洲色图av天堂| 午夜亚洲福利在线播放| 成人高潮视频无遮挡免费网站| 极品教师在线免费播放| 久久久久免费精品人妻一区二区| 亚洲性夜色夜夜综合| 免费一级毛片在线播放高清视频| 亚洲无线在线观看| 日日摸夜夜添夜夜添av毛片 | 国产黄片美女视频| 精品久久久久久久末码| 欧美日韩综合久久久久久 | 一个人看视频在线观看www免费| 亚洲精品国产成人久久av| 成人欧美大片| 琪琪午夜伦伦电影理论片6080| 啦啦啦韩国在线观看视频| 国产成年人精品一区二区| 国产 一区 欧美 日韩| 国产免费av片在线观看野外av| 亚洲精品乱码久久久v下载方式| 成人三级黄色视频| 中文字幕免费在线视频6| 美女免费视频网站| 亚洲avbb在线观看| 黄片wwwwww| 一级av片app| 国产一区二区三区av在线 | 熟女人妻精品中文字幕| 国产精品久久久久久亚洲av鲁大| 国产中年淑女户外野战色| 成人av一区二区三区在线看| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩无卡精品| 国产单亲对白刺激| eeuss影院久久| 91精品国产九色| 婷婷丁香在线五月| 国产真实伦视频高清在线观看 | 国产精品av视频在线免费观看| 成人亚洲精品av一区二区| 男人舔女人下体高潮全视频| 直男gayav资源| 热99在线观看视频| 级片在线观看| 五月玫瑰六月丁香| 亚洲av中文字字幕乱码综合| 91久久精品国产一区二区三区| 国产av麻豆久久久久久久| 国产一区二区在线av高清观看| av专区在线播放| 亚洲国产欧美人成| 欧美+日韩+精品| 免费观看人在逋| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 国产视频一区二区在线看| 毛片女人毛片| 久久精品人妻少妇| 免费看日本二区| 18禁黄网站禁片免费观看直播| 色av中文字幕| 久久久久久久久大av| 国产在线精品亚洲第一网站| 变态另类丝袜制服| 日日摸夜夜添夜夜添av毛片 | 成人特级黄色片久久久久久久| 亚洲精品久久国产高清桃花| 国产免费一级a男人的天堂| 亚洲成人久久爱视频| 欧美+亚洲+日韩+国产| 桃色一区二区三区在线观看| www日本黄色视频网| 欧美另类亚洲清纯唯美| 久久久久久伊人网av| 中文在线观看免费www的网站| 99热只有精品国产| 国产精品野战在线观看| 美女xxoo啪啪120秒动态图| 国产伦一二天堂av在线观看| 日韩av在线大香蕉| 久久香蕉精品热| 1024手机看黄色片| 午夜福利高清视频| 久久99热6这里只有精品| 婷婷色综合大香蕉| 国产av一区在线观看免费| www.色视频.com| 淫妇啪啪啪对白视频| 国产成人福利小说| 搡老熟女国产l中国老女人| 搡老妇女老女人老熟妇| 精品不卡国产一区二区三区| 亚洲国产色片| 久久人人精品亚洲av| 老司机午夜福利在线观看视频| 最新中文字幕久久久久| 男女那种视频在线观看| 精品午夜福利在线看| 可以在线观看毛片的网站| 国产精品爽爽va在线观看网站| 久久精品国产亚洲av涩爱 | 一夜夜www| 看免费成人av毛片| 国产精华一区二区三区| 亚洲 国产 在线| 女人被狂操c到高潮| 久久精品久久久久久噜噜老黄 | 国产一区二区亚洲精品在线观看| 熟妇人妻久久中文字幕3abv| 一边摸一边抽搐一进一小说| 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 午夜福利欧美成人| 又黄又爽又刺激的免费视频.| 久久久久久国产a免费观看| 男女啪啪激烈高潮av片| 桃红色精品国产亚洲av| 老熟妇乱子伦视频在线观看| 久久精品综合一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 最好的美女福利视频网| 麻豆成人午夜福利视频| 熟女电影av网| 丝袜美腿在线中文| 内射极品少妇av片p| 极品教师在线免费播放| 久久久久久久久中文| 精品久久久久久久久av| 老师上课跳d突然被开到最大视频| 热99在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| av在线观看视频网站免费| 熟女人妻精品中文字幕| 免费在线观看影片大全网站| 桃红色精品国产亚洲av| 小说图片视频综合网站| videossex国产| 亚洲精品一卡2卡三卡4卡5卡| 91久久精品国产一区二区成人| 内射极品少妇av片p| 国产黄色小视频在线观看| 亚洲无线观看免费| 嫩草影院精品99| 国产精品人妻久久久久久| 亚洲av日韩精品久久久久久密| 久久久久久久久中文| 欧美zozozo另类| 欧美日韩黄片免| 久久精品国产亚洲网站| 国内揄拍国产精品人妻在线| 99热精品在线国产| 欧美成人性av电影在线观看| 精品一区二区三区视频在线观看免费| 国产亚洲欧美98| 丰满的人妻完整版| 人人妻人人看人人澡| 亚洲七黄色美女视频| 亚洲中文字幕日韩| 国产亚洲欧美98| 99热网站在线观看| 天堂av国产一区二区熟女人妻| 免费观看人在逋| 国产熟女欧美一区二区| 99热6这里只有精品| 国产一区二区在线观看日韩| 免费在线观看日本一区| 校园人妻丝袜中文字幕| 不卡一级毛片| 国产单亲对白刺激| 亚洲18禁久久av| 伦精品一区二区三区| 99riav亚洲国产免费| 人人妻,人人澡人人爽秒播| 亚洲图色成人| 变态另类丝袜制服| 伊人久久精品亚洲午夜| 日韩亚洲欧美综合| 亚洲欧美日韩无卡精品| 亚洲在线观看片| 一级毛片久久久久久久久女| 国产精品亚洲一级av第二区| 韩国av在线不卡| 黄色女人牲交| 国产在线男女| 久久午夜福利片| 大型黄色视频在线免费观看| 亚洲图色成人| 国产一区二区三区在线臀色熟女| 亚洲专区国产一区二区| 国产中年淑女户外野战色| 黄色女人牲交| 亚洲成av人片在线播放无| 国模一区二区三区四区视频| 国产av在哪里看| 又爽又黄a免费视频| 综合色av麻豆| 一级黄色大片毛片| 麻豆精品久久久久久蜜桃| 国产精品av视频在线免费观看| 精品久久久久久久人妻蜜臀av| 嫩草影院入口| 国产精品综合久久久久久久免费| 最新在线观看一区二区三区| www日本黄色视频网| 精品久久久久久久末码| 国内精品久久久久精免费| 国产成人福利小说| 少妇的逼水好多| 一级毛片久久久久久久久女| 国产毛片a区久久久久| 69av精品久久久久久| 国产高清不卡午夜福利| 日本欧美国产在线视频| 日本免费一区二区三区高清不卡| 成人鲁丝片一二三区免费| 久久精品久久久久久噜噜老黄 | 国产精品自产拍在线观看55亚洲| 午夜日韩欧美国产| 午夜亚洲福利在线播放| 亚洲av五月六月丁香网| 亚洲精品影视一区二区三区av| 夜夜看夜夜爽夜夜摸| 午夜精品久久久久久毛片777| 日日啪夜夜撸| 特大巨黑吊av在线直播| 99热网站在线观看| 精品国内亚洲2022精品成人| 久久久久久大精品| 一个人看视频在线观看www免费| 久久亚洲真实| 亚洲中文字幕日韩| 国产人妻一区二区三区在| 国产三级中文精品| 免费看av在线观看网站| 久久这里只有精品中国| 午夜精品一区二区三区免费看| 久久精品综合一区二区三区| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久av不卡| 亚洲第一区二区三区不卡| 一进一出抽搐动态| 国产午夜精品久久久久久一区二区三区 | 91久久精品国产一区二区成人| 日韩av不卡免费在线播放| 中文字幕人妻熟人妻熟丝袜美| 超碰av人人做人人爽久久| 国产精品一区二区在线观看99| 99热这里只有是精品在线观看| xxx大片免费视频| 少妇高潮的动态图| 日韩一区二区三区影片| 一本色道久久久久久精品综合| 免费大片黄手机在线观看| 在线免费观看不下载黄p国产| 少妇人妻久久综合中文| 99国产精品免费福利视频| 国产人妻一区二区三区在| 99热网站在线观看| 国产一区二区三区综合在线观看 | 日韩精品有码人妻一区| 国产又色又爽无遮挡免| 亚洲色图综合在线观看| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区国产| 这个男人来自地球电影免费观看 | 亚洲欧美精品自产自拍| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 中文字幕制服av| 免费不卡的大黄色大毛片视频在线观看|