• <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男天堂| 在线观看美女被高潮喷水网站| 黑人高潮一二区| 亚洲一区二区三区欧美精品| 精品久久久噜噜| 狠狠精品人妻久久久久久综合| 麻豆成人午夜福利视频| a级毛片免费高清观看在线播放| 天堂中文最新版在线下载| 亚洲av成人精品一区久久| 香蕉精品网在线| 精品国产国语对白av| 精品国产国语对白av| 国产成人免费无遮挡视频| 亚洲情色 制服丝袜| 一级毛片电影观看| 啦啦啦在线观看免费高清www| 最近最新中文字幕免费大全7| 亚洲国产av新网站| av免费观看日本| 久久久午夜欧美精品| 久久国产精品大桥未久av | 亚洲av中文av极速乱| 久久免费观看电影| 亚洲av中文av极速乱| 成人毛片a级毛片在线播放| 免费观看av网站的网址| 国产有黄有色有爽视频| 精品国产露脸久久av麻豆| 国产成人精品福利久久| 亚洲久久久国产精品| 国产探花极品一区二区| 中文资源天堂在线| 80岁老熟妇乱子伦牲交| 亚洲精华国产精华液的使用体验| 国产伦精品一区二区三区四那| 国产精品麻豆人妻色哟哟久久| 久久久久久久国产电影| 久热久热在线精品观看| 久久人人爽人人片av| 日韩免费高清中文字幕av| 在线观看免费高清a一片| 国产成人a∨麻豆精品| 欧美国产精品一级二级三级 | av线在线观看网站| 晚上一个人看的免费电影| 女的被弄到高潮叫床怎么办| 婷婷色av中文字幕| 国产精品国产三级国产专区5o| 午夜激情福利司机影院| 免费观看的影片在线观看| 中国国产av一级| 91在线精品国自产拍蜜月| 亚洲精品国产成人久久av| 97在线视频观看| 国产亚洲午夜精品一区二区久久| 欧美精品国产亚洲| 午夜久久久在线观看| 精品久久国产蜜桃| 伦理电影免费视频| 久久久久人妻精品一区果冻| 搡老乐熟女国产| 视频中文字幕在线观看| 在线观看美女被高潮喷水网站| 少妇猛男粗大的猛烈进出视频| 国语对白做爰xxxⅹ性视频网站| 日韩一区二区视频免费看| 免费看日本二区| av一本久久久久| 男男h啪啪无遮挡| 亚洲久久久国产精品| 精品少妇久久久久久888优播| 少妇的逼水好多| kizo精华| 欧美激情国产日韩精品一区| 国产高清不卡午夜福利| 在线看a的网站| av在线观看视频网站免费| 日韩精品有码人妻一区| 色哟哟·www| av一本久久久久| 狂野欧美激情性xxxx在线观看| 啦啦啦视频在线资源免费观看| 毛片一级片免费看久久久久| 国产亚洲精品久久久com| 热99国产精品久久久久久7| 黄色视频在线播放观看不卡| 国产中年淑女户外野战色| 纯流量卡能插随身wifi吗| 色视频www国产| 一级黄片播放器| 妹子高潮喷水视频| 亚洲精品一区蜜桃| av国产精品久久久久影院| 9色porny在线观看| 日韩精品有码人妻一区| 精品久久久久久久久av| 国产成人免费观看mmmm| 多毛熟女@视频| 久久人人爽av亚洲精品天堂| 免费观看无遮挡的男女| 性色avwww在线观看| 亚洲四区av| 在线免费观看不下载黄p国产| 久久人人爽人人片av| 99热这里只有是精品在线观看| 三上悠亚av全集在线观看 | 成人亚洲精品一区在线观看| 免费av不卡在线播放| 久久人妻熟女aⅴ| 久久久久久久久久久久大奶| 午夜av观看不卡| 看非洲黑人一级黄片| 国产精品人妻久久久影院| 多毛熟女@视频| 色视频在线一区二区三区| 国产在视频线精品| a 毛片基地| 麻豆乱淫一区二区| 国产色爽女视频免费观看| kizo精华| 老司机影院毛片| 久久久亚洲精品成人影院| 日产精品乱码卡一卡2卡三| 国产成人91sexporn| 精品久久久久久久久av| 久久99热这里只频精品6学生| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区国产| 少妇熟女欧美另类| 麻豆精品久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 亚洲国产色片| 欧美精品亚洲一区二区| 亚洲精品中文字幕在线视频 | 欧美高清成人免费视频www| 色5月婷婷丁香| 久久久久久久久久久免费av| 美女大奶头黄色视频| 精品久久久久久电影网| 99九九线精品视频在线观看视频| 欧美日韩视频高清一区二区三区二| 99九九在线精品视频 | 国产成人精品久久久久久| 男女免费视频国产| 亚洲精品,欧美精品| 久久久久视频综合| 国产精品女同一区二区软件| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 国产日韩欧美视频二区| 亚洲不卡免费看| 80岁老熟妇乱子伦牲交| 欧美激情极品国产一区二区三区 | 搡女人真爽免费视频火全软件| 黑人高潮一二区| 99九九线精品视频在线观看视频| 嫩草影院入口| 久久久久精品性色| 久久久久久伊人网av| 久久精品国产自在天天线| 成人国产av品久久久| 国产熟女欧美一区二区| 寂寞人妻少妇视频99o| 国产一区二区三区av在线| 国产一区亚洲一区在线观看| 青春草亚洲视频在线观看| 国精品久久久久久国模美| av国产久精品久网站免费入址| 免费黄色在线免费观看| 丰满饥渴人妻一区二区三| 免费黄频网站在线观看国产| 人妻一区二区av| 久久久久久久久久人人人人人人| 久久免费观看电影| 亚洲美女黄色视频免费看| 国产欧美另类精品又又久久亚洲欧美| 丰满迷人的少妇在线观看| av卡一久久| 久久精品久久久久久久性| 国产色婷婷99| 黄色配什么色好看| 欧美精品一区二区大全| 亚洲国产av新网站| 久久久国产一区二区| 日本爱情动作片www.在线观看| 国产片特级美女逼逼视频| 国产精品伦人一区二区| 丝袜喷水一区| 亚洲中文av在线| 婷婷色麻豆天堂久久| 亚洲精品乱久久久久久| 一级毛片我不卡| 国产精品一二三区在线看| 久久精品国产鲁丝片午夜精品| 国产av一区二区精品久久| 久久鲁丝午夜福利片| av.在线天堂| 亚洲真实伦在线观看| 国产精品秋霞免费鲁丝片| 在线观看三级黄色| 国产av国产精品国产| 我要看黄色一级片免费的| 国产精品国产三级国产专区5o| 九九爱精品视频在线观看| 在线观看免费视频网站a站| 看免费成人av毛片| 国产成人精品久久久久久| 国产黄频视频在线观看| 精品午夜福利在线看| 美女xxoo啪啪120秒动态图| av卡一久久| 国产精品国产三级国产专区5o| 国产深夜福利视频在线观看| 成年女人在线观看亚洲视频| 精品人妻熟女毛片av久久网站| 男女免费视频国产| 欧美亚洲 丝袜 人妻 在线| av免费在线看不卡| 日本av手机在线免费观看| 亚洲精品日本国产第一区| 美女视频免费永久观看网站| 一个人免费看片子| 18+在线观看网站| 亚洲国产精品一区二区三区在线| 亚洲精品日韩av片在线观看| 日韩大片免费观看网站| 综合色丁香网| 下体分泌物呈黄色| 久久久久久久国产电影| 中文在线观看免费www的网站| 欧美变态另类bdsm刘玥| 五月开心婷婷网| 蜜臀久久99精品久久宅男| 一本—道久久a久久精品蜜桃钙片| 一级毛片久久久久久久久女| 啦啦啦视频在线资源免费观看| 麻豆乱淫一区二区| 亚洲精品一区蜜桃| 高清视频免费观看一区二区| 免费观看av网站的网址| 看非洲黑人一级黄片| 亚洲av电影在线观看一区二区三区| 日韩欧美精品免费久久| 亚洲人成网站在线播| 日韩在线高清观看一区二区三区| 永久网站在线| 精品国产一区二区久久| 久久ye,这里只有精品| 国产精品久久久久久精品电影小说| 中文字幕人妻熟人妻熟丝袜美| 亚洲久久久国产精品| 亚洲精品视频女| 大码成人一级视频| 中文字幕av电影在线播放| 99热全是精品| 乱码一卡2卡4卡精品| 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 国产精品国产三级国产av玫瑰| 国精品久久久久久国模美| 秋霞在线观看毛片| 人妻 亚洲 视频| 人妻夜夜爽99麻豆av| 国产男人的电影天堂91| 日韩免费高清中文字幕av| 麻豆成人午夜福利视频| 成人二区视频| 在线观看免费视频网站a站| 天天操日日干夜夜撸| 亚洲精品色激情综合| 国产成人a∨麻豆精品| 亚洲av在线观看美女高潮| 国产午夜精品一二区理论片| 99视频精品全部免费 在线| 亚洲一区二区三区欧美精品| √禁漫天堂资源中文www| 欧美xxxx性猛交bbbb| 久久久欧美国产精品| 男女边摸边吃奶| 91久久精品国产一区二区三区| 国产中年淑女户外野战色| 日韩一区二区视频免费看| 国产成人精品无人区| av又黄又爽大尺度在线免费看| 春色校园在线视频观看| 午夜影院在线不卡| 欧美日韩一区二区视频在线观看视频在线| 91成人精品电影| kizo精华| 国产午夜精品久久久久久一区二区三区| 一级毛片久久久久久久久女| 香蕉精品网在线| 国产白丝娇喘喷水9色精品| 欧美日韩精品成人综合77777| 啦啦啦视频在线资源免费观看| 国产一级毛片在线| 亚洲av二区三区四区| av免费在线看不卡| 亚洲四区av| 午夜免费观看性视频| 这个男人来自地球电影免费观看 | 亚洲第一av免费看| 国产成人精品久久久久久| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 日韩电影二区| 亚洲av欧美aⅴ国产| 一级毛片电影观看| 韩国av在线不卡| 中国美白少妇内射xxxbb| 只有这里有精品99| 免费观看性生交大片5| 内地一区二区视频在线| 肉色欧美久久久久久久蜜桃| 欧美精品国产亚洲| 老司机亚洲免费影院| 成年女人在线观看亚洲视频| 亚洲av欧美aⅴ国产| 国产亚洲最大av| 韩国高清视频一区二区三区| 久久鲁丝午夜福利片| 最近最新中文字幕免费大全7| 在线 av 中文字幕| 国产又色又爽无遮挡免| 亚洲美女黄色视频免费看| 久久 成人 亚洲| 老熟女久久久| 日本-黄色视频高清免费观看| 亚洲av福利一区| 乱系列少妇在线播放| 在线精品无人区一区二区三| 久久久亚洲精品成人影院| 免费观看a级毛片全部| 国产成人a∨麻豆精品| 久热久热在线精品观看| 亚洲经典国产精华液单| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 国产爽快片一区二区三区| 久久国产精品大桥未久av | 黄色视频在线播放观看不卡| 国产亚洲av片在线观看秒播厂| 欧美区成人在线视频| 丰满人妻一区二区三区视频av| 免费观看的影片在线观看| 国产色爽女视频免费观看| 国产女主播在线喷水免费视频网站| 少妇的逼好多水| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区综合在线观看 | 99热这里只有是精品在线观看| 内射极品少妇av片p| 免费播放大片免费观看视频在线观看| 久久精品久久精品一区二区三区| 国产精品久久久久久精品古装| 赤兔流量卡办理| 日本与韩国留学比较| 久久韩国三级中文字幕| 国产在线一区二区三区精| 在线天堂最新版资源| 国产成人免费观看mmmm| 中文乱码字字幕精品一区二区三区| 一级毛片aaaaaa免费看小| 一区二区三区精品91| 色视频在线一区二区三区| 丝袜在线中文字幕| 看免费成人av毛片| 亚洲av男天堂| 麻豆成人午夜福利视频| 久久久久久久亚洲中文字幕| 亚洲精品乱码久久久久久按摩| 综合色丁香网| 18禁动态无遮挡网站| 人体艺术视频欧美日本| 精品一区在线观看国产| 韩国高清视频一区二区三区| h视频一区二区三区| 久久99精品国语久久久| 一二三四中文在线观看免费高清| 乱人伦中国视频| 中文天堂在线官网| 久久综合国产亚洲精品| 亚洲精品中文字幕在线视频 | 在线观看av片永久免费下载| 亚洲欧洲精品一区二区精品久久久 | av免费观看日本| 王馨瑶露胸无遮挡在线观看| 中文天堂在线官网| 国产亚洲最大av| 成年av动漫网址| 全区人妻精品视频| 亚洲欧洲精品一区二区精品久久久 | freevideosex欧美| 亚洲怡红院男人天堂| 日韩亚洲欧美综合| 国产色婷婷99| 91久久精品电影网| 欧美精品国产亚洲| 色视频www国产| 伦理电影大哥的女人| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 亚洲无线观看免费| 国产一区二区在线观看日韩| 热re99久久国产66热| 极品少妇高潮喷水抽搐| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 欧美成人午夜免费资源| 交换朋友夫妻互换小说| 亚洲情色 制服丝袜| 亚洲精品乱码久久久v下载方式| 久久99一区二区三区| 精品国产一区二区久久| 国产成人免费无遮挡视频| 中文字幕久久专区| 日本午夜av视频| 国内精品宾馆在线| 日韩在线高清观看一区二区三区| 国内揄拍国产精品人妻在线| 国产 精品1| 日日摸夜夜添夜夜添av毛片| 日韩精品有码人妻一区| 国产老妇伦熟女老妇高清| 国产探花极品一区二区| 国产亚洲91精品色在线| 欧美丝袜亚洲另类| 亚洲av二区三区四区| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美日韩在线播放 | 啦啦啦在线观看免费高清www| 午夜91福利影院| 男人添女人高潮全过程视频| 一区二区三区精品91| 精品国产国语对白av| 女的被弄到高潮叫床怎么办| 中文资源天堂在线| 日本av免费视频播放| 中文精品一卡2卡3卡4更新| 久久久久国产精品人妻一区二区| 亚洲中文av在线| 亚洲自偷自拍三级| 亚洲欧美日韩另类电影网站| 国产在线免费精品| 日韩制服骚丝袜av| 亚洲熟女精品中文字幕| 一本久久精品| 99re6热这里在线精品视频| xxx大片免费视频| 极品教师在线视频| 最近中文字幕2019免费版| 噜噜噜噜噜久久久久久91| 又爽又黄a免费视频| 91aial.com中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 亚洲精华国产精华液的使用体验| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 熟女电影av网| 久久人人爽av亚洲精品天堂| 成人综合一区亚洲| 精品久久久精品久久久| 日韩欧美一区视频在线观看 | 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 97在线视频观看| 久久精品国产亚洲av天美| 亚洲av国产av综合av卡| 欧美日韩国产mv在线观看视频| 赤兔流量卡办理| 99久久精品热视频| 少妇人妻精品综合一区二区| av福利片在线观看| 成人漫画全彩无遮挡| 99热这里只有是精品50| 3wmmmm亚洲av在线观看| 国产精品99久久99久久久不卡 | 成年人午夜在线观看视频| 国产伦理片在线播放av一区| 国产欧美另类精品又又久久亚洲欧美| 免费看日本二区| 看十八女毛片水多多多| 狂野欧美激情性xxxx在线观看| 极品少妇高潮喷水抽搐| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产| 免费av不卡在线播放| 9色porny在线观看| 日日啪夜夜爽| 国产极品天堂在线| 六月丁香七月| 99久久精品一区二区三区| 日日撸夜夜添| 91精品国产九色| 只有这里有精品99| 99久国产av精品国产电影| 亚洲中文av在线| 久久人人爽av亚洲精品天堂| 在线观看一区二区三区激情| 3wmmmm亚洲av在线观看| 少妇 在线观看| 女性生殖器流出的白浆| 18禁在线播放成人免费| 精品一品国产午夜福利视频| av在线观看视频网站免费| 国语对白做爰xxxⅹ性视频网站| 日日摸夜夜添夜夜添av毛片| 久久久国产一区二区| 中文在线观看免费www的网站| 国内揄拍国产精品人妻在线| 狂野欧美激情性bbbbbb| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 精品酒店卫生间| 中文精品一卡2卡3卡4更新| 国产精品国产三级国产专区5o| 纵有疾风起免费观看全集完整版| 免费观看性生交大片5| kizo精华| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 草草在线视频免费看| 免费看日本二区| 亚洲精品,欧美精品| 精品人妻熟女av久视频| 插逼视频在线观看| 黄色怎么调成土黄色| 国产在线男女| 九色成人免费人妻av| 免费黄频网站在线观看国产| 一区二区av电影网| 欧美 日韩 精品 国产| 91午夜精品亚洲一区二区三区| av不卡在线播放| 亚洲成人av在线免费| 国产免费一区二区三区四区乱码| 国产探花极品一区二区| 欧美少妇被猛烈插入视频| 久久婷婷青草| 一级a做视频免费观看| 在线观看免费视频网站a站| 精品午夜福利在线看| 免费少妇av软件| 大片电影免费在线观看免费| 欧美97在线视频| 2018国产大陆天天弄谢| 日韩精品免费视频一区二区三区 | 97精品久久久久久久久久精品| 边亲边吃奶的免费视频| 插阴视频在线观看视频| 高清午夜精品一区二区三区| 国产在视频线精品| 全区人妻精品视频| 一本大道久久a久久精品| 你懂的网址亚洲精品在线观看| 人妻人人澡人人爽人人| 国产白丝娇喘喷水9色精品| 国模一区二区三区四区视频| 国产av精品麻豆| 夜夜骑夜夜射夜夜干| av天堂中文字幕网| 六月丁香七月| 人体艺术视频欧美日本| 春色校园在线视频观看| 肉色欧美久久久久久久蜜桃| 日本欧美国产在线视频| 国产高清不卡午夜福利| 亚洲国产最新在线播放| 肉色欧美久久久久久久蜜桃| 日本欧美国产在线视频| 国产女主播在线喷水免费视频网站| 女性被躁到高潮视频| 日韩欧美 国产精品| 亚洲欧美一区二区三区国产| 国产精品国产三级专区第一集| 国产视频内射| av免费观看日本| 九九爱精品视频在线观看| 日本vs欧美在线观看视频 | 国产伦精品一区二区三区视频9| 九九爱精品视频在线观看| 日本vs欧美在线观看视频 | 丰满人妻一区二区三区视频av| 久久久久久久亚洲中文字幕| 97超视频在线观看视频| 亚洲国产毛片av蜜桃av| 中文精品一卡2卡3卡4更新| 午夜视频国产福利| 久久综合国产亚洲精品| 亚洲久久久国产精品| 9色porny在线观看| 日韩,欧美,国产一区二区三区| 久久久久视频综合| 亚洲成人一二三区av| 精品人妻熟女av久视频| 日韩欧美 国产精品| 一级毛片我不卡| 99re6热这里在线精品视频| 日本爱情动作片www.在线观看| 大香蕉久久网| 有码 亚洲区| 久久婷婷青草| 欧美性感艳星| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 亚洲第一av免费看| 久久久国产一区二区| 免费看不卡的av| 国产精品国产三级国产av玫瑰| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 日韩视频在线欧美| 久久久久视频综合| 国产午夜精品一二区理论片|