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

    鋼箱計(jì)算失效應(yīng)變的沖擊試驗(yàn)

    2014-08-08 13:02:03王君杰李軍孟德巍
    關(guān)鍵詞:鋼箱計(jì)算精度

    王君杰+李軍+孟德巍

    建筑科學(xué)與工程學(xué)報(bào)2014年文章編號(hào):16732049(2014)01005006

    收稿日期:20131107

    基金項(xiàng)目:國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(“九七三”計(jì)劃)項(xiàng)目(2013CB036305);國家自然科學(xué)基金項(xiàng)目(51278373);

    交通運(yùn)輸部西部交通建設(shè)科技項(xiàng)目(2007 318 822 34)

    作者簡介:王君杰(1962),男,遼寧本溪人,教授,博士研究生導(dǎo)師,工學(xué)博士,

    摘要:為研究船橋碰撞有限元分析中鋼板的網(wǎng)格尺寸與鋼材計(jì)算失效應(yīng)變之間的關(guān)系,進(jìn)行了3個(gè)鋼箱模型的落錘沖擊試驗(yàn)。采用LSDYNA軟件對(duì)試驗(yàn)?zāi)P瓦M(jìn)行了有限元建模和碰撞計(jì)算,并與試驗(yàn)結(jié)果進(jìn)行了對(duì)比。定義了一個(gè)相關(guān)系數(shù)來反映試驗(yàn)結(jié)果與計(jì)算結(jié)果之間的相關(guān)性,并據(jù)此定義了與網(wǎng)格尺寸相關(guān)的計(jì)算失效應(yīng)變合理取值區(qū)間。研究結(jié)果表明:為得到合理精度的計(jì)算結(jié)果,鋼板的計(jì)算失效應(yīng)變的取值應(yīng)隨鋼板網(wǎng)格尺寸變化,使用大的網(wǎng)格尺寸時(shí)應(yīng)采用小的失效應(yīng)變,使用小的網(wǎng)格尺寸時(shí)應(yīng)采用大的失效應(yīng)變;將計(jì)算失效應(yīng)變合理取值區(qū)間與自適應(yīng)網(wǎng)格剖分技術(shù)結(jié)合,可以在保證計(jì)算精度的同時(shí),提高計(jì)算效率。

    關(guān)鍵詞:鋼箱;網(wǎng)格尺寸;失效應(yīng)變;計(jì)算精度;自適應(yīng)網(wǎng)格剖分;沖擊試驗(yàn)

    中圖分類號(hào):U661.72文獻(xiàn)標(biāo)志碼:A

    Impact Test on Computational Failure Strain of Steel BoxesWANG Junjie, LI Jun, MENG Dewei

    (State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China)Abstract: In order to investigate the relations between the meshing size of steel plates and computational failure strain of steel in shipbridge collision finite element analysis, the drop hammer impact tests of three steel boxes were carried out. The finite element model was built and impact computation was conducted for test models using LSDYNA software, and the computational results were compared with the test results. A correlation coefficient was defined to obtain a reasonable failure strain interval related to the meshing size of the steel boxes. The research results show that in order to get reasonable computational accurcy, the values of computational failure strainshould vary with the meshing size of steel plates. The larger failure strain should be used for the smaller meshing size, and the smaller failure strain should be used for the larger meshing size. The computational accuracy and computational efficiency can be obtained at the same time through combining the reasonable computational failure strain interval and the adaptive meshing technology.

    Key words: steel box; meshing size; failure strain; computational accuracy; adaptive meshing; impact test

    0引言

    對(duì)于橋梁船撞設(shè)計(jì),能夠把握船撞橋的碰撞力特征和破壞形態(tài)具有重要的實(shí)際意義。隨著計(jì)算機(jī)技術(shù)的快速發(fā)展,特別是有限元技術(shù)的日益進(jìn)步和成熟,數(shù)值模擬分析在橋梁船撞設(shè)計(jì)和防撞設(shè)施設(shè)計(jì)上逐漸得到了廣泛的應(yīng)用。與此同時(shí),數(shù)值仿真計(jì)算結(jié)果的可靠性和如何使用有限元技術(shù)對(duì)此類碰撞問題進(jìn)行有效的模擬日益成為關(guān)注的焦點(diǎn)。建立有限元模型是有限元分析過程的關(guān)鍵,而網(wǎng)格劃分是建立有限元模型的中心工作,模型的合理性很大程度上可以通過所劃分的網(wǎng)格形式體現(xiàn)出來。

    對(duì)于固體碰撞問題,研究發(fā)現(xiàn)網(wǎng)格尺寸與材料失效應(yīng)變的取值存在一定的相關(guān)性[15]。Lehmann等[1]、Kitamura等[2]、Paik等[3]分別基于一系列的碰撞試驗(yàn)或鋼板拉伸試驗(yàn)總結(jié)得出:必要的數(shù)值失效應(yīng)變(在有限元模型中調(diào)整該值以匹配試驗(yàn)數(shù)據(jù))是網(wǎng)格尺寸的函數(shù),總的趨勢是大尺寸的網(wǎng)格需要定義較小的數(shù)值失效應(yīng)變。Pedersen[6],Kitamura等[2]和高震等[7]將這種現(xiàn)象解釋為:較大的單元降低了高應(yīng)力點(diǎn)的應(yīng)力,以致撕裂不能及時(shí)發(fā)生,提高了結(jié)構(gòu)抗力;大的網(wǎng)格尺寸使用較小的失效應(yīng)變是考慮了裂縫、侵蝕和沖擊荷載等的影響。然而,由于碰撞的實(shí)際行為狀態(tài)非常重要,因而一個(gè)算例或一個(gè)試驗(yàn)的結(jié)果很難直接應(yīng)用于其他碰撞情況。船橋碰撞中的不確定因素很多,而這些因素又會(huì)對(duì)碰撞結(jié)果產(chǎn)生極大的影響。換言之,當(dāng)分析的問題稍有不同時(shí),就可能需要對(duì)網(wǎng)格尺寸與失效應(yīng)變的關(guān)系進(jìn)行調(diào)整。

    本文中筆者基于3次碰撞試驗(yàn)研究了網(wǎng)格尺寸和失效應(yīng)變?nèi)≈抵g的關(guān)系,對(duì)碰撞破壞問題中的網(wǎng)格劃分方法進(jìn)行了進(jìn)一步的探討。

    1模型設(shè)計(jì)與測試

    本文中以船舶正向撞擊鋼結(jié)構(gòu)防撞設(shè)施為背景進(jìn)行了3次試驗(yàn),其初始能量(碰撞前系統(tǒng)的能量)比約為1∶2∶3。由于試驗(yàn)中沖頭的初始動(dòng)能基本全部由被撞鋼箱吸收,因此3次試驗(yàn)中被撞鋼箱吸收的能量比例也約為1∶2∶3(表1),并由此來實(shí)現(xiàn)被撞鋼箱結(jié)構(gòu)不同程度的破壞,試驗(yàn)安裝見圖1。試驗(yàn)中的主動(dòng)撞擊結(jié)構(gòu)由上導(dǎo)向板、落錘、下導(dǎo)向板及撞擊鐵塊組成。試驗(yàn)時(shí),由電動(dòng)葫蘆通過掛彈鉤將主動(dòng)撞擊結(jié)構(gòu)提升至預(yù)定的跌落高度,然后掛彈鉤通電釋放主動(dòng)撞擊結(jié)構(gòu),主動(dòng)撞擊結(jié)構(gòu)自由跌落并沖擊鋼結(jié)構(gòu)防撞設(shè)施縮尺模型,從而使被撞鋼箱出現(xiàn)不同程度的破壞。3次試驗(yàn)的結(jié)構(gòu)布置見圖2。

    試驗(yàn)中的被撞鋼箱具體尺寸參見文獻(xiàn)[8]。由于被撞鋼箱鋼板較薄(4.62 mm),焊接過程中產(chǎn)生的殘余應(yīng)力將對(duì)結(jié)構(gòu)的性能產(chǎn)生很大影響,故在焊接結(jié)束后對(duì)被撞鋼箱進(jìn)行了鋼板回火處理、焊縫質(zhì)

    表1試驗(yàn)工況

    Tab.1Testing Cases試驗(yàn)

    編號(hào)撞擊質(zhì)量/kg跌落高度/m沖擊速度/

    endprint

    (m·s-1)初始動(dòng)能/

    (kN·m)19 7681.85.94172 307.529 7684.18.97392 478.2314 0284.59.39618 634.8圖1試驗(yàn)安裝

    Fig.1Test Installation圖2試驗(yàn)的結(jié)構(gòu)布置

    Fig.2Structural Arrangements of Tests量檢查和鋼板的靜態(tài)拉伸試驗(yàn)等。

    試驗(yàn)中,對(duì)沖擊過程中主動(dòng)撞擊錘的加速度和被撞鋼箱上選取點(diǎn)的應(yīng)變響應(yīng)進(jìn)行測量,對(duì)沖擊結(jié)束后被撞鋼箱上一些選取的點(diǎn)進(jìn)行人工位移測量,使用高速攝影設(shè)備對(duì)沖擊過程中被撞鋼箱的變形情況進(jìn)行輔助性的記錄。測量系統(tǒng)見圖3。

    圖3測量系統(tǒng)

    Fig.3Measurement System2數(shù)值計(jì)算方法

    在大量分析和比較的基礎(chǔ)上,本文中使用如下的計(jì)算參數(shù):選擇在汽車碰撞分析中廣泛使用的多段線性塑性模型來代表鋼材在沖擊作用下的力學(xué)屬性;采用CowperSymonds公式[9][(σYd/σY=1+(ε/C)q,其中,σYd為材料動(dòng)態(tài)失效應(yīng)變,σY為材料靜態(tài)失效應(yīng)變,ε為應(yīng)變率,C,q均為材料常數(shù)]來考慮應(yīng)變率的影響,并結(jié)合粘塑性公式來減少考慮應(yīng)變率時(shí)的響應(yīng)噪聲,本試驗(yàn)中被撞鋼箱為Q235軟鋼,C=40.4,q=1/5;采用材料的實(shí)際應(yīng)力應(yīng)變關(guān)系;鋼材的失效模式采用最大有效塑性應(yīng)變失效模式,失效模式表述為εpeff≥εfailure,即當(dāng)有限元模型中單元的應(yīng)變超過設(shè)定值時(shí),單元失效,失效后的單元從模型中刪除。3數(shù)值計(jì)算結(jié)果符合度檢驗(yàn)方法

    由于沖頭相對(duì)于被撞鋼箱來說剛度極大,因此在試驗(yàn)中,沖頭可以近似地視為剛體,故在得到?jīng)_頭的加速度時(shí)程后,可以根據(jù)牛頓第二定律得到碰撞過程中的碰撞力時(shí)程。

    在進(jìn)行船舶碰撞的有限元仿真分析時(shí),主要關(guān)注以下2個(gè)方面的結(jié)果:①船舶與橋梁或防撞結(jié)構(gòu)之間的撞擊力特征,包括碰撞力輪廓和碰撞力峰值;②船舶或橋梁的破壞模式、破壞情形。

    為了討論不同網(wǎng)格尺寸對(duì)應(yīng)的最佳失效應(yīng)變?nèi)≈?,本文中定義描述試驗(yàn)加速度結(jié)果與數(shù)值仿真結(jié)果之間差別的2個(gè)指標(biāo)或準(zhǔn)則:①計(jì)算值與測量值之間的相關(guān)系數(shù),采用Pearson相關(guān)系數(shù)r;②峰值加速度計(jì)算值與測量值之間的相對(duì)誤差e。r,e的計(jì)算方法分別為

    r=XY-XYN/

    (X2-(X)2N)(Y2-(Y)2N)(1)

    e=amax-a′maxamax(2)

    式中:X,Y為2個(gè)數(shù)值序列;amax為最大加速度測量值; a′max為最大加速度計(jì)算值。

    3次試驗(yàn)測量得到的加速度試驗(yàn)結(jié)果見圖4,其中,g為重力加速度。試驗(yàn)結(jié)果表明,開始段均具有明顯的周期約為0.003 8 s的波動(dòng)段,且該波動(dòng)段在第1次試驗(yàn)中響應(yīng)最大,在第3次試驗(yàn)中響應(yīng)最小,

    圖4加速度試驗(yàn)結(jié)果及其修正

    Fig.4Acceleration Test Results and Amendment初步判定該響應(yīng)與試驗(yàn)場地地基有關(guān)[10]。本文中采用HilbertHuang[1112]原理修正了第1次和第2次試驗(yàn)的波動(dòng)段加速度試驗(yàn)結(jié)果。4網(wǎng)格剖分與失效應(yīng)變的關(guān)系

    考慮被撞鋼箱上的殼單元邊長為2 cm的情形,將單元的失效應(yīng)變?chǔ)谭謩e取為0.12,0.15,0.20,0.25,0.30,研究單元失效應(yīng)變對(duì)加速度數(shù)值模擬結(jié)果的影響,結(jié)果見圖5。

    圖5失效應(yīng)變對(duì)加速度計(jì)算結(jié)果的影響

    Fig.5Effects of Failure Strain on Acceleration

    Computational Results由圖5可以看出,失效應(yīng)變?nèi)≈挡煌?,得到的?jì)算結(jié)果差別很大,這說明失效應(yīng)變的合理取值對(duì)船橋碰撞數(shù)值模擬計(jì)算結(jié)果的合理性和可靠性是十分重要的。

    本文中對(duì)被撞鋼箱劃分了7種網(wǎng)格尺寸,網(wǎng)格單元邊長l分別為0.5,1,1.4,2,2.5,3,5 cm,并采用多種失效應(yīng)變(最小失效應(yīng)變0.05,最大失效應(yīng)變0.6)來研究失效應(yīng)變和網(wǎng)格尺寸的關(guān)系。相關(guān)系數(shù)與加速度峰值相對(duì)誤差的計(jì)算結(jié)果見圖6,7。

    圖6相關(guān)系數(shù)計(jì)算結(jié)果

    Fig.6Computational Results of Correlation Coefficient從圖6,7可以看出:

    (1)失效應(yīng)變?nèi)≈禐?.05時(shí)(很小時(shí)),不同網(wǎng)格尺寸下的計(jì)算結(jié)果均與試驗(yàn)結(jié)果差別很大,總體上來說,加速度計(jì)算結(jié)果與試驗(yàn)結(jié)果之間的相關(guān)系數(shù)較小,且可能出現(xiàn)非常大的峰值加速度,這說明失效應(yīng)變?nèi)≈颠^小不能合理預(yù)測鋼箱發(fā)生塑性損傷時(shí)的沖擊反應(yīng)。

    (2)隨著失效應(yīng)變的增大,大尺寸網(wǎng)格可以率先計(jì)算得到較好的結(jié)果,隨著失效應(yīng)變繼續(xù)增加,小尺寸的網(wǎng)格依次得到精度較好的計(jì)算結(jié)果。

    (3)存在一個(gè)失效應(yīng)變的取值區(qū)間[εLf,εUf],εLf為區(qū)間的下界,εUf為區(qū)間的上界,當(dāng)失效應(yīng)變?cè)诖藚^(qū)間取值時(shí)可獲得較高精度的計(jì)算結(jié)果。這個(gè)取值區(qū)間的范圍隨網(wǎng)格尺寸變化:網(wǎng)格尺寸大,失效應(yīng)變?nèi)≈祬^(qū)間寬度?。环粗?,網(wǎng)格尺寸小,則失效應(yīng)變區(qū)間寬度大。

    (4)當(dāng)網(wǎng)格尺寸足夠小時(shí),對(duì)于本試驗(yàn),網(wǎng)格尺寸不宜大于1 cm,失效應(yīng)變?nèi)?.35~0.50是合理的,但是同時(shí)也可以看出,合理的失效應(yīng)變?nèi)≈蹬c鋼箱遭受的沖擊破壞程度可能有關(guān)。

    圖8為最優(yōu)失效應(yīng)變范圍。從圖8可以看出,對(duì)于同一網(wǎng)格的不同失效應(yīng)變情形,相關(guān)系數(shù)總是先增大再減小。據(jù)此定義[εLf,εUf]為合理取值區(qū)間。合理取值區(qū)間的下界εLf和上界εUf根據(jù)相關(guān)系數(shù)r>0.9的條件確定。圖8中綜合了3次試驗(yàn)的加速度結(jié)果,得到了適用于3次試驗(yàn)的失效應(yīng)變合理取值區(qū)間的下界和上界。圖8中顯示合理區(qū)間的下界和上界均隨著網(wǎng)格尺寸的增大而逐步變小,上界變化相對(duì)較快,區(qū)間范圍隨網(wǎng)格尺寸的增大而變小。

    從圖8還可以看出,如果采用大網(wǎng)格模型進(jìn)行分析,則必須謹(jǐn)慎地定義失效應(yīng)變;如果采用小網(wǎng)格模型進(jìn)行分析,則可以適當(dāng)?shù)貙⑹?yīng)變?nèi)〉么笠稽c(diǎn),以計(jì)算結(jié)果逐步穩(wěn)定時(shí)為宜。

    有限元仿真模擬的另一項(xiàng)重要內(nèi)容是對(duì)結(jié)構(gòu)的變形模式或破壞形式進(jìn)行捕捉。圖9中給出了第3次試驗(yàn)中被撞鋼箱背面鋼板出現(xiàn)的褶皺和0.5 cm以及5 cm網(wǎng)格計(jì)算對(duì)該褶皺捕捉情況的對(duì)比。由此可見,0.5 cm網(wǎng)格對(duì)該褶皺的描述非常好,且結(jié)構(gòu)變形光滑平順,而5 cm網(wǎng)格則完全沒有表現(xiàn)出該處出現(xiàn)的褶皺。這表明當(dāng)需要研究結(jié)構(gòu)的變形模式或破壞模式時(shí),應(yīng)該劃分尺寸較小的網(wǎng)格。5結(jié)語

    (1)網(wǎng)格尺寸與失效應(yīng)變是相關(guān)的,即單元的失效應(yīng)變依賴于網(wǎng)格的尺寸,使用大尺寸網(wǎng)格時(shí)應(yīng)采用較小的失效應(yīng)變,且失效應(yīng)變的定義較為敏感;使用小尺寸網(wǎng)格時(shí)應(yīng)采用較大的失效應(yīng)變,且失效應(yīng)變的定義不敏感,因此,根據(jù)相關(guān)系數(shù)定義了合理失效應(yīng)變區(qū)間。網(wǎng)格大則合理失效應(yīng)變?nèi)≈祬^(qū)間小,網(wǎng)格小則合理失效應(yīng)變?nèi)≈祬^(qū)間大。

    (2)為保證計(jì)算結(jié)果的精度,建議計(jì)算時(shí)可逐步細(xì)化模型中網(wǎng)格的尺寸并對(duì)失效應(yīng)變的取值進(jìn)行較大幅度的變化,當(dāng)失效應(yīng)變的取值在較大范圍內(nèi)變圖7加速度峰值相對(duì)誤差計(jì)算結(jié)果

    Fig.7Computational Results of Relative Error of Peak Acceleration圖8最優(yōu)失效應(yīng)變范圍

    Fig.8Optimum Failure Strain Interval化而對(duì)計(jì)算結(jié)果的影響較小時(shí),可以認(rèn)為已經(jīng)獲得了較好的計(jì)算模型。

    (3)采用自適應(yīng)網(wǎng)格剖分并結(jié)合最優(yōu)失效應(yīng)變區(qū)間概念進(jìn)行計(jì)算分析,可以在獲得同等計(jì)算精度的同時(shí)節(jié)省大量的建模時(shí)間和計(jì)算分析時(shí)間。

    endprint

    (4)需要注意的是,本文中的結(jié)論是在縮尺模型試驗(yàn)中獲得的,其適用性還有待在足尺試驗(yàn)或?qū)嶋H船橋碰撞事故中得到進(jìn)一步驗(yàn)證。圖9試驗(yàn)與計(jì)算褶皺的對(duì)比

    Fig.9Comparisons of Folds Between Test and Computational Results參考文獻(xiàn):

    References:[1]LEHMANN E,PESCHMANN J.Energy Absorption by the Steel Structure of Ships in the Event of Collisions[J].Marine Structures,2002,15(4/5):429441.

    [2]KITAMURA O.FEM Approach to the Simulation of Collision and Grounding Damage[J].Marine Structures,2002,15(4/5):403428.

    [3]PAIK J K,AMDAHL J,BARLTROP N.Collision and Grounding[C]//MANSOURE A E,ERTEKIN R C.Proceedings of the 15th International Ship and Offshore Structures Congress.San Diego:ISSC,2003:71107.

    [4]SERVIS D,SAMUELIDES M,LOUKA T,et al.The Implementation of Finite Element Codes for the Simulation of Shipship Collision[J].Journal of Ship Research,2002,46(4):239247.

    [5]NAAR H,KUJALA P,SIMONSEN B C,et al.Comparison of the Crashworthiness of Various Bottom and Side Structures[J].Marine Structures,2002,15(4/5):443460.

    [6]PEDERSEN P T.Ship Impacts:Bow Collisions[J].International Journal of Impact Engineerin,1993,13(2):163187.

    [7]高震,顧永寧,胡志強(qiáng).結(jié)構(gòu)沖擊試驗(yàn)的校準(zhǔn)計(jì)算[J].船舶力學(xué),2005,9(2):7782.

    GAO Zhen,GU Yongning,HU Zhiqiang.Benchmark Study of Structural Impact Test[J].Journal of Ship Mechanincs,2005,9(2):7782.

    [8]李軍.沖擊數(shù)值模擬可靠性的試驗(yàn)檢驗(yàn)[D].上海:同濟(jì)大學(xué),2009.

    LI Jun.Experimental Examination of the Trustworthiness of Impact Numerical Simulation[D].Shanghai:Tongji University,2009.

    [9]COWPER G R,Symonds P S.Strain Hardening and Strain Rate Effects in the Impact Loading of Cantilever Beams[R].Providence:Brown University,1958.

    [10]華南理工大學(xué),東南大學(xué),浙江大學(xué),等.地基與基礎(chǔ)[M].2版.北京:中國建筑工業(yè)出版社,1991.

    South China University of Technology, Southest University,Zhejiang University,et al.Soils and Foundations[M].2nd ed.Beijing:China Architecture & Building Press,1991.

    [11]HUANG N E,WU M C,LONG S R,et al.A Confidence Limit for the Empirical Mode Decomposition and Hilbert Spectral Analysis[J].Proceedings of the Royal Society A,2003,459(2037):23172345.

    [12]HUANG N E,SHEN Z,LONG S R,et al.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Nonstationary Time Series Analysis[J].Proceedings of the Royal Society A,1998,454(1971):903995.(上接第15頁)

    associated with CSA S13607 will be greater than that of GB 50018—2002.

    (4) The difference of the nominal axial strength between the two standards is primarily influenced by the flange widthtothickness ratio. For typical Csection wall studs investigated herein, the difference on the nominal axial strength is primarily influenced by the flange widthtothickness ratio wf/t. If the flange widthtothickness ratio wf/t is not less than 17.8, the difference on the nominal axial strength is dominated by the difference of flange effective width, and the nominal axial strength evaluated by GB 50018—2002 is less than that of CSA S13607, with the maximum difference being 31.9%. However, when wf/t is approximately less than 17.8, then the difference on the nominal axial strength is primarily governed by the difference of the web effective width and the nominal axial strength evaluated by GB 50018—2002 is slightly greater than that of CSA S13607, with the maximum difference being 8.9%.References:[1]CSA S13607,North American Specification for the Design of Coldformed Steel Structural Members[S].

    [2]GB 50018—2002,Technical Code of Coldformed Thinwall Steel Structures[S].

    [3]JGJ 227—2011,Technical Specification for Lowrise Coldformed Thinwall Steel Buildings[S].

    [4]CHEN J.Stability of Steel Structures:Theory and Design[M].Beijing:Science Press,2008.

    [5]YU W W,LABOUBE R A.Coldformed Steel Design[M].New York:John Wiley & Sons,2010.

    [6]ZHOU X H,WANG S J.Stability Theory and Its Applications of Thinwalled Members[M].Beijing:Science Press,2009.

    [7]XU L.Advanced Structural Steel Design[R].Waterloo:University of Waterloo,2012.

    [8]Canadian Institute of Steel Construction.Handbook of Steel Construction[M].10th ed.Markham:Canadian Institute of Steel Construction,2010.

    [9]Canadian Sheet Steel Building Institute (CSSBI).Lightweight Steel Frame Metric Section Properties[R].Cambridge:Canadian Sheet Steel Building Institute,2011.

    endprint

    (4)需要注意的是,本文中的結(jié)論是在縮尺模型試驗(yàn)中獲得的,其適用性還有待在足尺試驗(yàn)或?qū)嶋H船橋碰撞事故中得到進(jìn)一步驗(yàn)證。圖9試驗(yàn)與計(jì)算褶皺的對(duì)比

    Fig.9Comparisons of Folds Between Test and Computational Results參考文獻(xiàn):

    References:[1]LEHMANN E,PESCHMANN J.Energy Absorption by the Steel Structure of Ships in the Event of Collisions[J].Marine Structures,2002,15(4/5):429441.

    [2]KITAMURA O.FEM Approach to the Simulation of Collision and Grounding Damage[J].Marine Structures,2002,15(4/5):403428.

    [3]PAIK J K,AMDAHL J,BARLTROP N.Collision and Grounding[C]//MANSOURE A E,ERTEKIN R C.Proceedings of the 15th International Ship and Offshore Structures Congress.San Diego:ISSC,2003:71107.

    [4]SERVIS D,SAMUELIDES M,LOUKA T,et al.The Implementation of Finite Element Codes for the Simulation of Shipship Collision[J].Journal of Ship Research,2002,46(4):239247.

    [5]NAAR H,KUJALA P,SIMONSEN B C,et al.Comparison of the Crashworthiness of Various Bottom and Side Structures[J].Marine Structures,2002,15(4/5):443460.

    [6]PEDERSEN P T.Ship Impacts:Bow Collisions[J].International Journal of Impact Engineerin,1993,13(2):163187.

    [7]高震,顧永寧,胡志強(qiáng).結(jié)構(gòu)沖擊試驗(yàn)的校準(zhǔn)計(jì)算[J].船舶力學(xué),2005,9(2):7782.

    GAO Zhen,GU Yongning,HU Zhiqiang.Benchmark Study of Structural Impact Test[J].Journal of Ship Mechanincs,2005,9(2):7782.

    [8]李軍.沖擊數(shù)值模擬可靠性的試驗(yàn)檢驗(yàn)[D].上海:同濟(jì)大學(xué),2009.

    LI Jun.Experimental Examination of the Trustworthiness of Impact Numerical Simulation[D].Shanghai:Tongji University,2009.

    [9]COWPER G R,Symonds P S.Strain Hardening and Strain Rate Effects in the Impact Loading of Cantilever Beams[R].Providence:Brown University,1958.

    [10]華南理工大學(xué),東南大學(xué),浙江大學(xué),等.地基與基礎(chǔ)[M].2版.北京:中國建筑工業(yè)出版社,1991.

    South China University of Technology, Southest University,Zhejiang University,et al.Soils and Foundations[M].2nd ed.Beijing:China Architecture & Building Press,1991.

    [11]HUANG N E,WU M C,LONG S R,et al.A Confidence Limit for the Empirical Mode Decomposition and Hilbert Spectral Analysis[J].Proceedings of the Royal Society A,2003,459(2037):23172345.

    [12]HUANG N E,SHEN Z,LONG S R,et al.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Nonstationary Time Series Analysis[J].Proceedings of the Royal Society A,1998,454(1971):903995.(上接第15頁)

    associated with CSA S13607 will be greater than that of GB 50018—2002.

    (4) The difference of the nominal axial strength between the two standards is primarily influenced by the flange widthtothickness ratio. For typical Csection wall studs investigated herein, the difference on the nominal axial strength is primarily influenced by the flange widthtothickness ratio wf/t. If the flange widthtothickness ratio wf/t is not less than 17.8, the difference on the nominal axial strength is dominated by the difference of flange effective width, and the nominal axial strength evaluated by GB 50018—2002 is less than that of CSA S13607, with the maximum difference being 31.9%. However, when wf/t is approximately less than 17.8, then the difference on the nominal axial strength is primarily governed by the difference of the web effective width and the nominal axial strength evaluated by GB 50018—2002 is slightly greater than that of CSA S13607, with the maximum difference being 8.9%.References:[1]CSA S13607,North American Specification for the Design of Coldformed Steel Structural Members[S].

    [2]GB 50018—2002,Technical Code of Coldformed Thinwall Steel Structures[S].

    [3]JGJ 227—2011,Technical Specification for Lowrise Coldformed Thinwall Steel Buildings[S].

    [4]CHEN J.Stability of Steel Structures:Theory and Design[M].Beijing:Science Press,2008.

    [5]YU W W,LABOUBE R A.Coldformed Steel Design[M].New York:John Wiley & Sons,2010.

    [6]ZHOU X H,WANG S J.Stability Theory and Its Applications of Thinwalled Members[M].Beijing:Science Press,2009.

    [7]XU L.Advanced Structural Steel Design[R].Waterloo:University of Waterloo,2012.

    [8]Canadian Institute of Steel Construction.Handbook of Steel Construction[M].10th ed.Markham:Canadian Institute of Steel Construction,2010.

    [9]Canadian Sheet Steel Building Institute (CSSBI).Lightweight Steel Frame Metric Section Properties[R].Cambridge:Canadian Sheet Steel Building Institute,2011.

    endprint

    (4)需要注意的是,本文中的結(jié)論是在縮尺模型試驗(yàn)中獲得的,其適用性還有待在足尺試驗(yàn)或?qū)嶋H船橋碰撞事故中得到進(jìn)一步驗(yàn)證。圖9試驗(yàn)與計(jì)算褶皺的對(duì)比

    Fig.9Comparisons of Folds Between Test and Computational Results參考文獻(xiàn):

    References:[1]LEHMANN E,PESCHMANN J.Energy Absorption by the Steel Structure of Ships in the Event of Collisions[J].Marine Structures,2002,15(4/5):429441.

    [2]KITAMURA O.FEM Approach to the Simulation of Collision and Grounding Damage[J].Marine Structures,2002,15(4/5):403428.

    [3]PAIK J K,AMDAHL J,BARLTROP N.Collision and Grounding[C]//MANSOURE A E,ERTEKIN R C.Proceedings of the 15th International Ship and Offshore Structures Congress.San Diego:ISSC,2003:71107.

    [4]SERVIS D,SAMUELIDES M,LOUKA T,et al.The Implementation of Finite Element Codes for the Simulation of Shipship Collision[J].Journal of Ship Research,2002,46(4):239247.

    [5]NAAR H,KUJALA P,SIMONSEN B C,et al.Comparison of the Crashworthiness of Various Bottom and Side Structures[J].Marine Structures,2002,15(4/5):443460.

    [6]PEDERSEN P T.Ship Impacts:Bow Collisions[J].International Journal of Impact Engineerin,1993,13(2):163187.

    [7]高震,顧永寧,胡志強(qiáng).結(jié)構(gòu)沖擊試驗(yàn)的校準(zhǔn)計(jì)算[J].船舶力學(xué),2005,9(2):7782.

    GAO Zhen,GU Yongning,HU Zhiqiang.Benchmark Study of Structural Impact Test[J].Journal of Ship Mechanincs,2005,9(2):7782.

    [8]李軍.沖擊數(shù)值模擬可靠性的試驗(yàn)檢驗(yàn)[D].上海:同濟(jì)大學(xué),2009.

    LI Jun.Experimental Examination of the Trustworthiness of Impact Numerical Simulation[D].Shanghai:Tongji University,2009.

    [9]COWPER G R,Symonds P S.Strain Hardening and Strain Rate Effects in the Impact Loading of Cantilever Beams[R].Providence:Brown University,1958.

    [10]華南理工大學(xué),東南大學(xué),浙江大學(xué),等.地基與基礎(chǔ)[M].2版.北京:中國建筑工業(yè)出版社,1991.

    South China University of Technology, Southest University,Zhejiang University,et al.Soils and Foundations[M].2nd ed.Beijing:China Architecture & Building Press,1991.

    [11]HUANG N E,WU M C,LONG S R,et al.A Confidence Limit for the Empirical Mode Decomposition and Hilbert Spectral Analysis[J].Proceedings of the Royal Society A,2003,459(2037):23172345.

    [12]HUANG N E,SHEN Z,LONG S R,et al.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Nonstationary Time Series Analysis[J].Proceedings of the Royal Society A,1998,454(1971):903995.(上接第15頁)

    associated with CSA S13607 will be greater than that of GB 50018—2002.

    (4) The difference of the nominal axial strength between the two standards is primarily influenced by the flange widthtothickness ratio. For typical Csection wall studs investigated herein, the difference on the nominal axial strength is primarily influenced by the flange widthtothickness ratio wf/t. If the flange widthtothickness ratio wf/t is not less than 17.8, the difference on the nominal axial strength is dominated by the difference of flange effective width, and the nominal axial strength evaluated by GB 50018—2002 is less than that of CSA S13607, with the maximum difference being 31.9%. However, when wf/t is approximately less than 17.8, then the difference on the nominal axial strength is primarily governed by the difference of the web effective width and the nominal axial strength evaluated by GB 50018—2002 is slightly greater than that of CSA S13607, with the maximum difference being 8.9%.References:[1]CSA S13607,North American Specification for the Design of Coldformed Steel Structural Members[S].

    [2]GB 50018—2002,Technical Code of Coldformed Thinwall Steel Structures[S].

    [3]JGJ 227—2011,Technical Specification for Lowrise Coldformed Thinwall Steel Buildings[S].

    [4]CHEN J.Stability of Steel Structures:Theory and Design[M].Beijing:Science Press,2008.

    [5]YU W W,LABOUBE R A.Coldformed Steel Design[M].New York:John Wiley & Sons,2010.

    [6]ZHOU X H,WANG S J.Stability Theory and Its Applications of Thinwalled Members[M].Beijing:Science Press,2009.

    [7]XU L.Advanced Structural Steel Design[R].Waterloo:University of Waterloo,2012.

    [8]Canadian Institute of Steel Construction.Handbook of Steel Construction[M].10th ed.Markham:Canadian Institute of Steel Construction,2010.

    [9]Canadian Sheet Steel Building Institute (CSSBI).Lightweight Steel Frame Metric Section Properties[R].Cambridge:Canadian Sheet Steel Building Institute,2011.

    endprint

    猜你喜歡
    鋼箱計(jì)算精度
    基于多矩形餅切函數(shù)的彈藥毀傷效能評(píng)估方法
    爆炸與沖擊(2024年3期)2024-11-01 00:00:00
    多跨連續(xù)鋼箱拱橋拱腳關(guān)鍵位置疲勞性能分析
    大跨度上承式鋼箱桁肋拱橋設(shè)計(jì)
    中承式鋼箱景觀拱橋總體設(shè)計(jì)與計(jì)算
    鋼箱內(nèi)部爆炸破壞的SPH 數(shù)值模擬*
    爆炸與沖擊(2019年5期)2019-06-21 02:54:42
    鋼箱加固鋼筋混凝土梁橋的有限元分析
    基于SHIPFLOW的散貨船阻力預(yù)報(bào)分析
    廣東造船(2018年4期)2018-10-09 06:08:18
    Abaqus多尺度有限元模型在橋梁工程中的應(yīng)用
    世界家苑(2018年1期)2018-04-27 11:42:06
    五跨連續(xù)中承式鋼箱拱橋設(shè)計(jì)
    上海公路(2018年4期)2018-03-21 05:57:28
    連續(xù)變寬箱形橋梁計(jì)算精度探究
    床上黄色一级片| 三级毛片av免费| 亚洲国产精品国产精品| av在线亚洲专区| 国产精品综合久久久久久久免费| 色综合亚洲欧美另类图片| 在线观看美女被高潮喷水网站| av女优亚洲男人天堂| 欧美极品一区二区三区四区| 国产欧美日韩精品一区二区| a级毛色黄片| 99热网站在线观看| 在线观看一区二区三区| 九九久久精品国产亚洲av麻豆| 九色成人免费人妻av| 久久综合国产亚洲精品| 国内精品宾馆在线| 久久久久性生活片| 国产久久久一区二区三区| 99久久成人亚洲精品观看| 日本一本二区三区精品| 国产午夜精品论理片| av专区在线播放| 中文字幕av在线有码专区| 伦精品一区二区三区| 免费av观看视频| 国产精品一区二区在线观看99 | 亚洲国产色片| 成人美女网站在线观看视频| ponron亚洲| 亚洲熟妇中文字幕五十中出| 久久久久九九精品影院| 久久久精品大字幕| 大型黄色视频在线免费观看| 噜噜噜噜噜久久久久久91| 最近2019中文字幕mv第一页| 久久久久久久久久成人| 精品久久久久久久久亚洲| 亚洲av一区综合| 夜夜夜夜夜久久久久| 亚洲欧美成人精品一区二区| 久久久午夜欧美精品| 久久久久久久午夜电影| 亚洲人与动物交配视频| 亚洲av中文字字幕乱码综合| 只有这里有精品99| 天美传媒精品一区二区| 人妻久久中文字幕网| 特级一级黄色大片| avwww免费| 少妇裸体淫交视频免费看高清| 日本撒尿小便嘘嘘汇集6| 亚洲精华国产精华液的使用体验 | 91精品国产九色| 国产成人一区二区在线| 亚洲成人中文字幕在线播放| 嫩草影院新地址| 成人二区视频| 亚洲色图av天堂| 婷婷色综合大香蕉| 国产精品久久久久久av不卡| 国产又黄又爽又无遮挡在线| 三级国产精品欧美在线观看| 成人特级av手机在线观看| 精品人妻熟女av久视频| 日韩在线高清观看一区二区三区| 欧美日韩精品成人综合77777| 国产精品嫩草影院av在线观看| 国产在线精品亚洲第一网站| 久久久久久九九精品二区国产| 国产成年人精品一区二区| 午夜精品在线福利| 欧美高清成人免费视频www| 亚洲精品久久久久久婷婷小说 | av卡一久久| 99久久人妻综合| 婷婷精品国产亚洲av| 少妇被粗大猛烈的视频| 99久国产av精品| 26uuu在线亚洲综合色| 亚洲国产欧美在线一区| 日本av手机在线免费观看| a级毛色黄片| 久久这里有精品视频免费| 日本成人三级电影网站| 午夜a级毛片| 中文精品一卡2卡3卡4更新| 日韩欧美精品免费久久| 蜜桃亚洲精品一区二区三区| 赤兔流量卡办理| 国产亚洲精品av在线| 精品一区二区三区人妻视频| 国产一区二区三区av在线 | 最近的中文字幕免费完整| 内地一区二区视频在线| 麻豆成人av视频| 日本黄大片高清| 免费电影在线观看免费观看| а√天堂www在线а√下载| 国产美女午夜福利| 男女那种视频在线观看| 国产毛片a区久久久久| 国产极品精品免费视频能看的| 男人舔女人下体高潮全视频| 我要搜黄色片| 最近2019中文字幕mv第一页| 又爽又黄a免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美+亚洲+日韩+国产| АⅤ资源中文在线天堂| 在现免费观看毛片| 欧美另类亚洲清纯唯美| 国产亚洲精品av在线| 国产极品精品免费视频能看的| 少妇被粗大猛烈的视频| 精品免费久久久久久久清纯| 久久精品夜色国产| 成人综合一区亚洲| av福利片在线观看| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 美女高潮的动态| 99精品在免费线老司机午夜| 亚洲国产高清在线一区二区三| 亚洲成人精品中文字幕电影| 少妇的逼好多水| 91久久精品国产一区二区成人| 人妻系列 视频| 大又大粗又爽又黄少妇毛片口| 少妇高潮的动态图| 日产精品乱码卡一卡2卡三| 色哟哟·www| 国产精品一区www在线观看| 午夜久久久久精精品| 免费看光身美女| 中文欧美无线码| 国产精品日韩av在线免费观看| 欧美激情在线99| 亚洲va在线va天堂va国产| 午夜激情福利司机影院| 亚洲国产精品国产精品| 日韩视频在线欧美| 亚洲性久久影院| av在线蜜桃| 深夜精品福利| 能在线免费看毛片的网站| 人妻系列 视频| 亚洲欧美中文字幕日韩二区| 久久久久网色| 国产精品一区www在线观看| 真实男女啪啪啪动态图| 国产伦理片在线播放av一区 | eeuss影院久久| 一进一出抽搐gif免费好疼| 欧美bdsm另类| 两个人视频免费观看高清| 国产亚洲欧美98| 成人午夜精彩视频在线观看| 国产老妇女一区| 国产 一区精品| 九色成人免费人妻av| 天天躁日日操中文字幕| 麻豆精品久久久久久蜜桃| 国产大屁股一区二区在线视频| 欧美激情国产日韩精品一区| 12—13女人毛片做爰片一| 久久精品国产自在天天线| 国产精品爽爽va在线观看网站| 国产真实乱freesex| 国内久久婷婷六月综合欲色啪| 国产免费男女视频| 日韩一本色道免费dvd| 久久精品国产亚洲av香蕉五月| av在线天堂中文字幕| 国产麻豆成人av免费视频| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| 少妇人妻精品综合一区二区 | 国产黄a三级三级三级人| 欧美丝袜亚洲另类| av天堂中文字幕网| 久久综合国产亚洲精品| 日本一本二区三区精品| 欧美bdsm另类| 欧美色欧美亚洲另类二区| 人人妻人人澡人人爽人人夜夜 | 日本黄大片高清| 青青草视频在线视频观看| 国产综合懂色| 久久精品国产自在天天线| 精品欧美国产一区二区三| 亚洲国产日韩欧美精品在线观看| 97超视频在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区三区四区免费观看| 久久久久久久午夜电影| 成年免费大片在线观看| 亚洲国产日韩欧美精品在线观看| 热99在线观看视频| 久久这里只有精品中国| 日韩制服骚丝袜av| 久久精品影院6| 免费人成视频x8x8入口观看| 国产黄色视频一区二区在线观看 | 国产成人aa在线观看| 天堂av国产一区二区熟女人妻| 亚洲无线观看免费| 日韩大尺度精品在线看网址| ponron亚洲| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 最后的刺客免费高清国语| 99热这里只有精品一区| 亚洲在久久综合| 久久精品91蜜桃| 九九久久精品国产亚洲av麻豆| 日韩av在线大香蕉| 欧美成人免费av一区二区三区| 波多野结衣巨乳人妻| 亚洲欧美日韩无卡精品| 亚洲最大成人手机在线| 久久亚洲国产成人精品v| 日韩欧美一区二区三区在线观看| 99久久无色码亚洲精品果冻| 国产高清三级在线| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av涩爱 | 成人特级av手机在线观看| 狠狠狠狠99中文字幕| 一级黄片播放器| a级一级毛片免费在线观看| 天堂中文最新版在线下载 | 特级一级黄色大片| 亚洲成人精品中文字幕电影| 又粗又硬又长又爽又黄的视频 | 性欧美人与动物交配| 中文资源天堂在线| 2021天堂中文幕一二区在线观| 一级二级三级毛片免费看| 一本久久中文字幕| 最近最新中文字幕大全电影3| 日韩成人伦理影院| 看非洲黑人一级黄片| 久久久久久伊人网av| 九九热线精品视视频播放| 美女xxoo啪啪120秒动态图| 深爱激情五月婷婷| 插阴视频在线观看视频| 欧美3d第一页| 日韩大尺度精品在线看网址| 亚洲欧美精品自产自拍| 91aial.com中文字幕在线观看| 欧美性猛交黑人性爽| 美女内射精品一级片tv| 久久99热6这里只有精品| 欧美激情国产日韩精品一区| 亚洲乱码一区二区免费版| 久久久久久久久大av| 欧美成人一区二区免费高清观看| 亚洲美女搞黄在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 成人国产麻豆网| 在线观看美女被高潮喷水网站| 夜夜夜夜夜久久久久| 日韩一区二区视频免费看| 少妇被粗大猛烈的视频| 干丝袜人妻中文字幕| 国产精品综合久久久久久久免费| 国内精品宾馆在线| 亚洲经典国产精华液单| 91久久精品国产一区二区成人| 亚洲成av人片在线播放无| 国产精品免费一区二区三区在线| 亚洲久久久久久中文字幕| av在线观看视频网站免费| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 99精品在免费线老司机午夜| 国产亚洲91精品色在线| 老师上课跳d突然被开到最大视频| 日本色播在线视频| 啦啦啦啦在线视频资源| 免费看a级黄色片| 国产老妇伦熟女老妇高清| 看十八女毛片水多多多| 美女xxoo啪啪120秒动态图| 国产探花在线观看一区二区| 深夜精品福利| 国产精品麻豆人妻色哟哟久久 | 亚洲av中文字字幕乱码综合| 97超视频在线观看视频| av黄色大香蕉| 黄片wwwwww| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 一区二区三区四区激情视频 | 色播亚洲综合网| 麻豆成人av视频| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 国产成年人精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 久久人人精品亚洲av| 三级毛片av免费| 1000部很黄的大片| 亚洲欧美日韩高清专用| 国产成人精品久久久久久| 欧美在线一区亚洲| 国产毛片a区久久久久| 我要搜黄色片| 久久精品国产亚洲av涩爱 | 大香蕉久久网| 欧美zozozo另类| 欧美人与善性xxx| 激情 狠狠 欧美| 国产精品野战在线观看| 国产成人福利小说| 国产欧美日韩精品一区二区| 男女边吃奶边做爰视频| 免费一级毛片在线播放高清视频| 国国产精品蜜臀av免费| 精品国内亚洲2022精品成人| АⅤ资源中文在线天堂| 成人综合一区亚洲| 床上黄色一级片| 精品久久久噜噜| 免费电影在线观看免费观看| 亚洲国产精品国产精品| 亚洲av男天堂| 国产精品麻豆人妻色哟哟久久 | 97人妻精品一区二区三区麻豆| 欧美色欧美亚洲另类二区| 久久精品人妻少妇| 国产中年淑女户外野战色| 12—13女人毛片做爰片一| 久久国内精品自在自线图片| 国内精品美女久久久久久| 一级毛片我不卡| 网址你懂的国产日韩在线| 国产黄片视频在线免费观看| 国产成人freesex在线| 免费观看人在逋| 尾随美女入室| 久久综合国产亚洲精品| 国产午夜精品论理片| 丰满乱子伦码专区| 久久这里有精品视频免费| 99在线人妻在线中文字幕| 久久久久久久亚洲中文字幕| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| АⅤ资源中文在线天堂| a级一级毛片免费在线观看| 少妇的逼好多水| 亚洲在线观看片| av免费在线看不卡| 国产 一区 欧美 日韩| 久久这里只有精品中国| 国产大屁股一区二区在线视频| 美女国产视频在线观看| 国产精品麻豆人妻色哟哟久久 | av在线观看视频网站免费| 秋霞在线观看毛片| 日韩欧美精品v在线| 女人被狂操c到高潮| 最好的美女福利视频网| 国产中年淑女户外野战色| 男女啪啪激烈高潮av片| 色视频www国产| 亚洲自拍偷在线| 国产单亲对白刺激| 男人舔女人下体高潮全视频| 高清在线视频一区二区三区 | 国产成人福利小说| 人人妻人人看人人澡| 一级毛片电影观看 | 夫妻性生交免费视频一级片| 日本一二三区视频观看| 国产一级毛片七仙女欲春2| 国内精品宾馆在线| 久久热精品热| 久久精品国产亚洲网站| 免费观看a级毛片全部| 精品日产1卡2卡| www日本黄色视频网| 久久草成人影院| www日本黄色视频网| 欧美另类亚洲清纯唯美| 菩萨蛮人人尽说江南好唐韦庄 | 人人妻人人看人人澡| 成人特级av手机在线观看| 国产成人午夜福利电影在线观看| 超碰av人人做人人爽久久| 国产午夜精品论理片| 日本av手机在线免费观看| 黄色配什么色好看| 国产精品久久视频播放| 久久精品91蜜桃| 99在线视频只有这里精品首页| 1000部很黄的大片| 久久精品国产亚洲网站| 日韩中字成人| 国产大屁股一区二区在线视频| 美女黄网站色视频| 成人欧美大片| or卡值多少钱| 久久中文看片网| 日韩视频在线欧美| 久久久精品欧美日韩精品| 亚洲无线观看免费| 卡戴珊不雅视频在线播放| 婷婷精品国产亚洲av| 高清毛片免费观看视频网站| 一进一出抽搐gif免费好疼| 亚洲成人久久爱视频| 99国产极品粉嫩在线观看| 亚洲成a人片在线一区二区| 99国产极品粉嫩在线观看| 欧美丝袜亚洲另类| 亚洲成人精品中文字幕电影| 欧美另类亚洲清纯唯美| 菩萨蛮人人尽说江南好唐韦庄 | 中文字幕久久专区| 亚洲内射少妇av| 99久久精品一区二区三区| 国产成人福利小说| 在线播放无遮挡| 精品免费久久久久久久清纯| 国产高清不卡午夜福利| 六月丁香七月| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费观看的www视频| 韩国av在线不卡| 国产精品久久久久久久电影| 99久久中文字幕三级久久日本| 麻豆乱淫一区二区| 欧美日韩一区二区视频在线观看视频在线 | 午夜视频国产福利| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 国产欧美日韩精品一区二区| av免费在线看不卡| 精品日产1卡2卡| 国产爱豆传媒在线观看| 能在线免费看毛片的网站| 精品少妇黑人巨大在线播放 | 欧美潮喷喷水| 久久久久久伊人网av| 啦啦啦观看免费观看视频高清| 国产精品三级大全| 免费看光身美女| 97热精品久久久久久| 亚洲欧美精品自产自拍| 99热这里只有是精品在线观看| 看非洲黑人一级黄片| 尤物成人国产欧美一区二区三区| 久久人人爽人人爽人人片va| av在线亚洲专区| 女人十人毛片免费观看3o分钟| 夫妻性生交免费视频一级片| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 美女内射精品一级片tv| 看免费成人av毛片| 久久综合国产亚洲精品| 国产伦理片在线播放av一区 | 美女 人体艺术 gogo| 一进一出抽搐gif免费好疼| 精品无人区乱码1区二区| eeuss影院久久| av在线老鸭窝| 国产视频首页在线观看| 久久久久久久久久久免费av| 麻豆国产av国片精品| 亚洲精品日韩在线中文字幕 | 免费不卡的大黄色大毛片视频在线观看 | 久久九九热精品免费| 久久热精品热| 成人二区视频| 免费不卡的大黄色大毛片视频在线观看 | 熟女人妻精品中文字幕| 久久99精品国语久久久| 亚洲久久久久久中文字幕| 国产精品一区二区三区四区久久| 日日摸夜夜添夜夜添av毛片| 69av精品久久久久久| 欧美三级亚洲精品| 日本撒尿小便嘘嘘汇集6| 26uuu在线亚洲综合色| 天堂网av新在线| 欧美性猛交╳xxx乱大交人| 成年版毛片免费区| 久99久视频精品免费| 国产av在哪里看| 色尼玛亚洲综合影院| 国产成年人精品一区二区| 在线播放无遮挡| 一区二区三区免费毛片| 久久精品人妻少妇| 人妻久久中文字幕网| 免费不卡的大黄色大毛片视频在线观看 | 欧美成人a在线观看| 99久久九九国产精品国产免费| 日本av手机在线免费观看| 国语自产精品视频在线第100页| av女优亚洲男人天堂| 91av网一区二区| av在线老鸭窝| 听说在线观看完整版免费高清| 国产精品一区二区在线观看99 | 国产男人的电影天堂91| 久久久欧美国产精品| 嫩草影院精品99| 国产精品1区2区在线观看.| 国产精品一区www在线观看| 久久欧美精品欧美久久欧美| 国产av一区在线观看免费| 床上黄色一级片| 国产成人freesex在线| 国内精品宾馆在线| 日本色播在线视频| 国产成人a区在线观看| 乱人视频在线观看| 久久久精品大字幕| 在线观看一区二区三区| 成人一区二区视频在线观看| 91在线精品国自产拍蜜月| 在线观看免费视频日本深夜| 99热6这里只有精品| 精华霜和精华液先用哪个| 99久久成人亚洲精品观看| 狂野欧美激情性xxxx在线观看| 国语自产精品视频在线第100页| 麻豆成人午夜福利视频| 在线观看一区二区三区| 麻豆一二三区av精品| 69av精品久久久久久| 男的添女的下面高潮视频| 国产精品久久久久久久久免| 国产av不卡久久| 精品少妇黑人巨大在线播放 | 在线a可以看的网站| 国产亚洲av嫩草精品影院| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全电影3| 听说在线观看完整版免费高清| 国产亚洲精品久久久久久毛片| 久久99热6这里只有精品| 日韩欧美三级三区| 国产精品不卡视频一区二区| 亚洲欧美清纯卡通| 欧美日本亚洲视频在线播放| 91精品一卡2卡3卡4卡| 能在线免费看毛片的网站| 2021天堂中文幕一二区在线观| 一本精品99久久精品77| 免费观看人在逋| 人妻夜夜爽99麻豆av| 中出人妻视频一区二区| 有码 亚洲区| 免费av不卡在线播放| 国产麻豆成人av免费视频| 久久婷婷人人爽人人干人人爱| 亚洲欧美成人精品一区二区| 免费人成视频x8x8入口观看| 国产爱豆传媒在线观看| 国内久久婷婷六月综合欲色啪| 久久久久久大精品| 人体艺术视频欧美日本| 干丝袜人妻中文字幕| 亚洲18禁久久av| 岛国毛片在线播放| 国产中年淑女户外野战色| 亚洲综合色惰| 精品少妇黑人巨大在线播放 | 亚洲人成网站在线观看播放| 中国美女看黄片| 免费看a级黄色片| 日本一本二区三区精品| 欧美zozozo另类| 亚洲av中文av极速乱| av黄色大香蕉| 嫩草影院新地址| 国内精品美女久久久久久| 99热这里只有精品一区| 听说在线观看完整版免费高清| 永久网站在线| 欧美成人精品欧美一级黄| 丰满乱子伦码专区| 亚洲电影在线观看av| 亚洲国产色片| 少妇人妻一区二区三区视频| 天堂影院成人在线观看| 乱人视频在线观看| 身体一侧抽搐| 99热这里只有是精品在线观看| 97人妻精品一区二区三区麻豆| 久久人妻av系列| 69人妻影院| 国产精品久久久久久亚洲av鲁大| 99热这里只有是精品50| 91在线精品国自产拍蜜月| 亚洲精品自拍成人| 国产色爽女视频免费观看| 岛国在线免费视频观看| 久久精品人妻少妇| 国产亚洲av片在线观看秒播厂 | 亚洲自拍偷在线| 国产探花在线观看一区二区| 欧美成人一区二区免费高清观看| 久久人人精品亚洲av| 日韩一区二区视频免费看| 久久午夜福利片| 欧美又色又爽又黄视频|