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

    橋梁固定式防撞鋼套箱能力曲線的解析計(jì)算方法

    2020-12-06 10:42:45樊偉張澤文申?yáng)|杰孫洋
    關(guān)鍵詞:有限元分析

    樊偉 張澤文 申?yáng)|杰 孫洋

    摘? ?要:基于船-船碰撞解析研究,分析了鋼套箱豎向外板、橫肋、豎桁等構(gòu)件的變形機(jī)理與破壞模式,討論了鋼套箱在不同撞擊位置下的撞擊力差異. 結(jié)果表明:由于鋼套箱與船舶舷側(cè)結(jié)構(gòu)的尺寸差異,其橫豎桁排布密、間距小,使得撞擊點(diǎn)處相鄰構(gòu)件能夠很快參與受力,不同撞擊位置下的撞擊力差異并不顯著. 結(jié)合不同外形船艏撞擊下的鋼套箱受力情形,建立了典型的碰撞場(chǎng)景模型,提出了同時(shí)適用于帶球艏船舶、楔形艏船舶撞擊下的鋼套箱抗撞性能的解析計(jì)算方法. 采用提出的解析計(jì)算方法,計(jì)算了固定式鋼套箱在不同外形船舶撞擊下的撞擊力-撞深曲線,通過(guò)與精細(xì)有限元分析結(jié)果對(duì)比表明:解析計(jì)算方法的結(jié)果與精細(xì)有限元分析結(jié)果較好地吻合,證明了解析方法的有效性.

    關(guān)鍵詞:鋼套箱;解析計(jì)算方法;防撞裝置設(shè)計(jì);有限元分析

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

    Analytical Method for Estimating Anti-collision Capacity Curve

    of Fix-typed Steel Fenders of Bridge Structures

    FAN Wei1,2?,ZHANG Zewen1,SHEN Dongjie1,SUN Yang1

    (1. College of Civil Engineering,Hunan University,Changsha 410082,China;

    2. Key Laboratory for Wind and Bridge Engineering (Hunan University) of Hunan Province,Changsha 410082,China)

    Abstract:In this context, based on the analytical methodology of the ship-ship collision, the deformation mechanisms and failure modes of structural members in steel fenders (e.g., outer plate, transverse, and vertical ribs) were investigated. The characteristics of impact forces under different impact locations were carefully discussed. The analysis results indicate that the adjacent members at the impact location could work cooperatively at the same time because of the close arrangement and the small space inside the steel fender, and the differences of impact forces under various impact positions would be not significant. Based on the resistance characteristics of steel fenders under difference bow shapes, the typical collision scenarios were defined. Accordingly, the analytical method was developed to estimate the anti-collision performance of steel fenders, which is suitable for both ship vessels with a bulbous bow and those with a raked bow. The impact force versus crush depth curves of the fixed-type steel fenders were obtained by using the proposed analytical method for different ship impacts. The analytical results were compared with those obtained from the finite element analysis. The results show that the analytical results are in good agreement with the finite element results, indicating the rationality of the analytical method developed in this study.

    Key words:steel fender;analytical method;anti-collision device design;finite element analysis (FEA)

    隨著航運(yùn)業(yè)的快速發(fā)展,船橋矛盾日益突出,發(fā)生惡性船撞橋事故的概率也隨之升高. 據(jù)統(tǒng)計(jì),船舶撞擊是位列洪水之后航道橋梁倒塌的第二主因[1],船舶撞毀橋梁的問(wèn)題已成為橋梁工程界面臨的尖銳問(wèn)題之一[2]. 因此在設(shè)計(jì)通航水域橋梁時(shí),必須考慮船撞問(wèn)題[3]. 我國(guó)《公路橋涵設(shè)計(jì)通用規(guī)范》(JTG D60—2015)中規(guī)定:“可能遭受船舶或漂流物撞擊的橋墩,應(yīng)考慮船舶或漂流物的撞擊作用,并應(yīng)設(shè)置警示標(biāo)志和必要的防撞設(shè)施.”[4]為避免由船舶撞毀橋梁所帶來(lái)的經(jīng)濟(jì)損失、人員傷亡,保護(hù)有可能遭遇船舶碰撞的橋梁下部結(jié)構(gòu),給橋梁布置合理有效的防撞裝置顯得尤為重要. 樊偉[5]對(duì)國(guó)內(nèi)外橋梁工程實(shí)例調(diào)研統(tǒng)計(jì)得出:美國(guó)、加拿大等國(guó)家多采用獨(dú)立式防撞系統(tǒng)作為橋梁防撞的主要措施,然而與之不同的是,出于對(duì)經(jīng)濟(jì)性和適用性的考慮,我國(guó)越來(lái)越多的建設(shè)橋梁均采用鋼結(jié)構(gòu)護(hù)舷(鋼套箱)形式的防撞裝置,如蘇通大橋[6]、金塘大橋[7]、武漢天興洲大橋[8]、臺(tái)州椒江二橋[9]及杭州灣跨海大橋[10]等. 在船橋碰撞事故中,橋墩防撞裝置(鋼套箱)是抵御撞擊的重要結(jié)構(gòu),因此對(duì)鋼套箱的抗撞性能研究應(yīng)給予足夠重視.

    船撞橋問(wèn)題的研究起源于船船相撞的研究,其研究方法也大多借鑒船船相撞計(jì)算理論研究方法[11],大致分為:經(jīng)驗(yàn)靜力法、試驗(yàn)方法、精細(xì)有限元數(shù)值模擬方法和解析計(jì)算方法.

    由于橋墩結(jié)構(gòu)抗沖擊的足尺試驗(yàn)代價(jià)太大,難以展開(kāi),因此這方面的試驗(yàn)研究并不多見(jiàn)[12]. 目前,國(guó)內(nèi)外多采用數(shù)值模擬方法研究橋梁防撞擊問(wèn)題. 如許薛軍等[13]通過(guò)有限元軟件建立船舶撞擊橋墩的仿真模型,分析防撞套箱、橡膠護(hù)舷、橋墩之間的能量轉(zhuǎn)換及防撞套箱撞擊部位的受力與變形. Jiang等[14]利用有限元技術(shù)對(duì)一個(gè)浮式防撞鋼套箱進(jìn)行了全面深入的性能評(píng)估. Wang等[15]對(duì)一種新型柔性防撞鋼套箱展開(kāi)數(shù)值模擬分析,并通過(guò)對(duì)縮尺模型進(jìn)行一組落錘試驗(yàn)來(lái)驗(yàn)證數(shù)值模擬方法的正確性. 有限元數(shù)值模擬方法精度高、適用性廣,是目前鋼套箱性能研究的最常見(jiàn)方法.

    解析計(jì)算方法的優(yōu)點(diǎn)為計(jì)算快速、所需參數(shù)少且計(jì)算結(jié)果相對(duì)準(zhǔn)確. 采用精細(xì)有限元技術(shù)費(fèi)時(shí)費(fèi)力,因而在防撞裝置設(shè)計(jì)中應(yīng)用受到限制,尤其在設(shè)計(jì)前期方案過(guò)程中,更不適合用于開(kāi)展系統(tǒng)而詳細(xì)的影響因素分析以及設(shè)計(jì)優(yōu)化等[5]. 相比于精細(xì)有限元數(shù)值模擬方法,在船橋碰撞中解析計(jì)算方法目前應(yīng)用較少,尤其是針對(duì)橋梁防護(hù)裝置. 以往船船碰撞解析計(jì)算研究中所采用船舶排水量多為萬(wàn)噸以上[16-17],鋼套箱與此類船舶舷側(cè)結(jié)構(gòu)存在著尺寸差異,且實(shí)際工程中造型多為流線型,因此,通過(guò)在船側(cè)結(jié)構(gòu)碰撞機(jī)理的原有解析研究基礎(chǔ)上建立適合鋼套箱的簡(jiǎn)化分析方法,從而準(zhǔn)確快速地推導(dǎo)出船舶撞擊鋼套箱的撞擊力(能量)-撞深曲線,將在橋梁防撞裝置的初步設(shè)計(jì)階段,大大提高設(shè)計(jì)效率.

    近年來(lái)學(xué)者們大多針對(duì)帶球艏船舶的撞擊展開(kāi)分析研究,而楔形艏作為一種常見(jiàn)的艏部形狀,與其相關(guān)的碰撞研究未得到足夠重視. 且鋼套箱的變形損傷與撞擊船艏的形狀有著密切的關(guān)系,作為兩種最為常見(jiàn)的船艏構(gòu)造,研究?jī)深惒煌白矒粝聵蛄悍雷惭b置的抗撞性能具有重要意義.

    本文基于國(guó)內(nèi)外學(xué)者對(duì)船舶舷側(cè)結(jié)構(gòu)中板材在不同受力情形、不同破壞模式下的解析計(jì)算研究,提出了一套合理、完整的應(yīng)用于船舶撞擊橋梁鋼結(jié)構(gòu)護(hù)舷(鋼套箱)抗撞性分析的解析計(jì)算方法,并針對(duì)帶球艏船舶、楔形艏船舶兩類不同外形船艏結(jié)構(gòu)的撞擊開(kāi)展研究,擴(kuò)大了解析計(jì)算方法的適用范圍. 同時(shí),采用精細(xì)有限元技術(shù)對(duì)船舶撞擊鋼套箱進(jìn)行精細(xì)化建模計(jì)算. 通過(guò)比較有限元軟件計(jì)算出的撞擊力(能量)-撞深曲線與解析計(jì)算結(jié)果,在保證簡(jiǎn)化分析方法快速性的同時(shí),驗(yàn)證了該方法在橋梁防撞裝置初步設(shè)計(jì)階段的準(zhǔn)確性、適用性.

    1? ?鋼套箱構(gòu)件變形機(jī)理分析

    傳統(tǒng)的鋼套箱結(jié)構(gòu)多為板架構(gòu)件,通常由豎向外板、頂板、底板、縱橫艙壁、加勁肋等組成,當(dāng)船舶開(kāi)始接觸鋼套箱并向前行進(jìn)時(shí),隨著撞擊深度的增大,各個(gè)構(gòu)件參與碰撞并通過(guò)自身的變形損傷(撕裂、斷裂、壓彎和屈曲等)來(lái)充分吸收撞擊能量.

    在船舶撞擊鋼套箱解析計(jì)算中,分別對(duì)豎向外板的撕裂、縱橫艙壁的屈曲、加勁肋的彎曲斷裂等構(gòu)件破壞行為的變形機(jī)理進(jìn)行解析計(jì)算,最后將各個(gè)構(gòu)件的吸能性能整合在一起,得到整個(gè)鋼套箱在撞擊過(guò)程中的能量-撞深曲線.

    1.1? ?外板遭遇剛性球體擠壓阻力計(jì)算

    Wang等[18]曾根據(jù)球鼻艏外形運(yùn)用塑性力學(xué)中的上限定理提出了外板抵抗球艏船撞擊的簡(jiǎn)化力學(xué)模型,如圖1所示.

    邊界為2R的圓形外板受半徑為r的剛性球體撞擊下,在其未破裂前撞擊力和撞擊深度關(guān)系如下:

    式中:σ0 為材料的流動(dòng)應(yīng)力,通常認(rèn)為是材料的屈服強(qiáng)度和極限強(qiáng)度總和的一半;R為外板的變形區(qū)域半徑;r為球艏簡(jiǎn)化后得到的剛性球半徑;t為板的厚度;Δ為撞擊深度.

    當(dāng)r遠(yuǎn)小于R時(shí),可將剛性球體視為一點(diǎn),其荷載形式為點(diǎn)荷載,如圖2所示. 上述公式可簡(jiǎn)化為式(3)形式.

    通常情況下,球鼻艏的球鼻部分較尖銳,采用式(3)計(jì)算外板變形阻力更為方便快捷.

    1.2? ?外板遭遇楔形剛體擠壓阻力計(jì)算

    孫斌等[19]針對(duì)楔形艏撞擊船舶舷側(cè)外板進(jìn)行了研究,將外板損傷變形的能量耗散分為膜拉伸能量耗散和兩邊塑性鉸線的彎曲變形能量耗散,建立了外板在楔形船艏撞擊下的解析計(jì)算公式:

    式中:L1和L2為撞擊點(diǎn)距離外板變形區(qū)域左右邊界的距離;H為變形區(qū)的高度;各角度標(biāo)注如圖3所示;M0和N0分別為單位長(zhǎng)度板條梁的極限塑性彎矩和極限塑性力,大小分別為:

    1.3? ?板的面內(nèi)撕裂阻力解析計(jì)算

    Ohtsubo等[20]對(duì)船舶在接地過(guò)程中船底外板遇到楔形體撕裂問(wèn)題做了研究. 假設(shè)板材的斷裂是由韌性破壞引起的,認(rèn)為在板材撕裂過(guò)程中塑性變形和摩擦是耗能的主要因素,而斷裂是次要的,從而推導(dǎo)出板材撕裂問(wèn)題的一種上限解,解釋了板材撕裂過(guò)程中的破壞機(jī)理.

    外板在受到角度為2θ的楔形體撞擊下,其撕裂過(guò)程中的阻力分為瞬時(shí)撕裂力Fp和摩擦力Ff,其中Ff為Fp在摩擦位移方向上的投影. 最終獲得的板材面內(nèi)撕裂過(guò)程中力-位移關(guān)系為:

    式中:σ0為流動(dòng)應(yīng)力;t為板厚;l為撕裂長(zhǎng)度;θ等于楔形體角度的1/2;μ是摩擦因數(shù),其值為0.15~0.3.

    1.4? ?板的穿透模型阻力解析計(jì)算

    Wang等[21]經(jīng)過(guò)一系列的試驗(yàn)證實(shí)了在穿透過(guò)程中外板破裂后仍具有一定的承載力,且其破壞模式與板的面內(nèi)撕裂阻力解析機(jī)理極為相似,均涉及到了被撞擊物前端尖角處較大的膜拉伸以及較遠(yuǎn)端的顯著彎曲. 根據(jù)公式(7)得到板破裂后其承載力與撞擊深度的關(guān)系式如下:

    式中:l為每一道裂縫的長(zhǎng)度;n為裂縫的條數(shù);(n-2)π/n為相鄰裂縫之間形成的角度. 當(dāng)撞擊物前端過(guò)于尖銳時(shí)(例如銳利的球艏或楔形艏),相鄰兩道裂縫形成的角度依附于前端的夾角2θ,即(n-2)π/2n=θ. 一般尖銳的前端撞擊外板會(huì)形成的裂縫數(shù)目為n=2,此時(shí)(n-2)π/2n=0并不適用,將θ代入原公式較為合適[17].

    1.5? ?肋板折疊模型阻力解析計(jì)算

    圖4所示為肋板折疊模型簡(jiǎn)圖,Liu等[22]在原有研究的啟發(fā)下進(jìn)一步建立了更為合理的解析模型. 提出了新的肋板簡(jiǎn)化分析方法,給出了肋板在壓碎過(guò)程中初始及后續(xù)折疊更為精確的阻力計(jì)算公式,將一個(gè)完整的折疊變形區(qū)域設(shè)為3H0,當(dāng)b1=b2=b時(shí),第一道折疊過(guò)程中平均壓碎力大小為:

    1.6? ?十字板壓碎模型解析計(jì)算

    鋼套箱內(nèi)部有非常多的橫肋、豎桁相互交錯(cuò),從而形成十字板的結(jié)構(gòu). Zhang[11]對(duì)十字板的阻力進(jìn)行了解析計(jì)算研究. 其破壞模式類似于肋板壓碎破壞(如圖5所示),將第一道折疊區(qū)域設(shè)為2H1,其中H1 = 1.103(tb)0.5,在原有解析計(jì)算的基礎(chǔ)上推導(dǎo)出了X形、T形、L形多種不同截面形式的縱、橫梁組合結(jié)構(gòu)受壓下在第一道折疊過(guò)程中的平均變形阻力.

    b為十字板交點(diǎn)到邊界的距離,即縱(橫)肋長(zhǎng)度的一半. 當(dāng)十字板結(jié)構(gòu)在第N0道折疊中發(fā)生破裂后,假設(shè)阻力維持恒定,其后續(xù)的平均壓碎力為:

    通過(guò)給定鋼材的斷裂應(yīng)變,可以計(jì)算出構(gòu)件在第幾道折疊中破裂,即N0的值.

    2? ?鋼套箱抗撞性能的解析計(jì)算方法

    2.1? ?計(jì)算方法基本假定

    在解析計(jì)算方法中通常會(huì)適當(dāng)引入一些合理的假設(shè)[23],這些基本假設(shè)簡(jiǎn)化了結(jié)構(gòu)真實(shí)受力下的復(fù)雜性,便于快速有效地建立結(jié)構(gòu)理論模型. 結(jié)合橋梁防撞鋼套箱的構(gòu)造特點(diǎn),歸納出適合于鋼套箱抗撞性分析中的基本假定如下:

    1)保守地假定船艏為剛體,依靠鋼套箱的變形來(lái)吸收能量.

    2)假設(shè)船艏以恒定速度沿著撞擊方向持續(xù)前

    進(jìn),不考慮速度變化.

    3)不考慮鋼套箱內(nèi)部各構(gòu)件的相互作用,假定彼此之間的阻力是獨(dú)立的.

    4)各構(gòu)件的力學(xué)模型中,假定以沒(méi)有發(fā)生變形損傷的其他構(gòu)件作為局部受力破壞的邊界.

    5)計(jì)算撞深最大為鋼套箱內(nèi)、外板間距. 在船舶撞擊鋼套箱的過(guò)程中,內(nèi)板已與橋墩緊密貼合,沒(méi)有變形空間.

    2.2? ?不同船艏外形下的碰撞情景

    不同的船艏外形在與鋼套箱的碰撞過(guò)程中接觸到其內(nèi)部各個(gè)構(gòu)件時(shí)的撞擊深度不同,意味著各個(gè)構(gòu)件參與受力并耗散能量的撞深時(shí)刻會(huì)隨著船艏的外形而發(fā)生變化;另外不同外形船艏撞擊下鋼套箱的損傷面積也會(huì)受顯著影響,其參與耗能的構(gòu)件數(shù)目、各個(gè)構(gòu)件損傷程度都會(huì)有明顯差異. 為了便于研究,采用了文獻(xiàn)[11]中提出的假定球鼻艏與楔形艏模型對(duì)碰撞情景進(jìn)行說(shuō)明,兩種船艏的簡(jiǎn)化模型如圖6所示.

    假定的球鼻艏模型主要參數(shù)包括上甲板楔形夾角2θ、船艏高度Hdeck、首柱傾角φ以及球鼻處半橢球體3個(gè)半徑RL、RV、RH . 由于半橢球體是最先接觸且碰撞鋼套箱的部分,隨著船艏行進(jìn),擠壓并刺穿鋼套箱,因此與半橢球體尺寸有關(guān)的3個(gè)半徑是最重要的外形參數(shù),假定其值均與Hdeck成比例,關(guān)系如下:

    依據(jù)不同尺寸的球鼻部分,將球艏船大致分為大球艏與小球艏兩類,以此來(lái)討論球鼻外形的影響因素. 圖7所示為兩類球鼻艏撞擊鋼套箱的情形,大、小球鼻除了3個(gè)半徑RL、RV、RH有明顯差異外,小球鼻艏首柱傾角φ1也大于大球鼻艏首柱傾角φ2 . 小球鼻艏在碰撞行進(jìn)過(guò)程中,隨著位移的持續(xù)增大,球鼻以上楔形部分將接觸并擠壓鋼套箱上甲板. 這使得簡(jiǎn)化分析計(jì)算分成了兩個(gè)主要階段,即球鼻的初始接觸碰撞與楔形體參與的后續(xù)碰撞. 對(duì)于大球鼻艏來(lái)說(shuō),由于鋼套箱內(nèi)、外板間距有限,即使在最大撞深時(shí)刻也僅有球鼻部分參與碰撞,表明在碰撞全過(guò)程中,楔形部分始終未與鋼套箱頂板接觸,這也簡(jiǎn)化了解析計(jì)算流程. 因此,定義大、小球鼻艏的區(qū)別在于其上部楔形體是否會(huì)參與到擠壓鋼套箱的過(guò)程中.

    假定的楔形艏模型主要參數(shù)包括上甲板楔形夾角2θ、船艏高度Hdeck、首柱傾角φ. 與球鼻艏模型中球鼻部分的半橢球體水平半徑RH類似,對(duì)于楔形艏撞擊鋼套箱場(chǎng)景下,楔形夾角2θ是至關(guān)重要的,2θ? 值越小,船艏接觸到豎桁的位移時(shí)刻就會(huì)越大,如圖8所示. 這也將直接影響到船艏撞擊下鋼套箱的構(gòu)件參與數(shù)量及損傷面積.

    2.3? ?鋼套箱被撞位置不同的討論

    鋼套箱的組成包括頂板、底板、豎向外板、豎向內(nèi)板、豎桁與橫肋交錯(cuò)組成的縱橫艙壁及加勁肋. 船艏與鋼套箱初始接觸位置的不同,可能會(huì)導(dǎo)致構(gòu)件變形區(qū)域的改變、參與受力的位移順序以及提供阻力的構(gòu)件數(shù)量的不同,最終使得鋼套箱整體吸能效果有一些差異. 例如文獻(xiàn)[17]對(duì)海上浮式生產(chǎn)儲(chǔ)油船(FPSO)遭遇撞擊的研究中選取了3個(gè)典型被撞位置,得到了差異較大的撞擊力-位移曲線,因此對(duì)鋼套箱初始受壓位置的討論也很有必要.

    針對(duì)鋼套箱結(jié)構(gòu)的立面構(gòu)造,給出了如圖9所示的4種不同點(diǎn)荷載的作用位置:縱橫板交界處、橫肋中點(diǎn)、豎桁中點(diǎn)、豎桁和橫肋組成的艙壁內(nèi)部某點(diǎn). 初始碰撞位置的不同對(duì)鋼套箱豎向外板及內(nèi)部各構(gòu)件的受力順序及參與情況會(huì)帶來(lái)一定影響.

    一般來(lái)說(shuō),無(wú)論是帶球艏船舶亦或是楔形艏船舶,其尖銳的船艏結(jié)構(gòu),在初始接觸的時(shí)刻均可視為點(diǎn)荷載. 另外,相比于排水量大于10 000 DWT的船舶舷側(cè)結(jié)構(gòu),鋼套箱縱橫艙壁間距較小,這意味著套箱豎向外板由于受到未變形的橫肋、豎桁約束,變形區(qū)域較小,在較小的撞深下外板就會(huì)破裂并進(jìn)入撕裂模式,采用公式(3)計(jì)算外板破裂前變形阻力在一定計(jì)算精度范圍內(nèi)是合理的簡(jiǎn)化.

    基本假定中表明各構(gòu)件以沒(méi)有發(fā)生變形損傷的其他構(gòu)件作為局部受力破壞的邊界. 由于受到未變形的縱、橫肋約束,對(duì)于外板受力而言,不同的撞擊點(diǎn)位置會(huì)帶來(lái)不同的初始變形區(qū)域,如圖10中陰影區(qū)域所示4個(gè)大小不一的矩形面積. 在楔形船撞擊下的外板破裂前阻力計(jì)算中,邊界條件取矩形短邊的一半,即公式(4)中L1和L2的取值會(huì)隨著變形區(qū)域的不同發(fā)生變化. 對(duì)于破裂位移時(shí)刻,后續(xù)公式(17)和(18)中l(wèi)0值不同,也會(huì)使得外板破裂進(jìn)入撕裂模式的臨界撞深有所差別. 這意味著初始變形區(qū)域越大,外板阻力越大,破裂時(shí)的臨界撞深也越大,參與耗能更多.

    同樣,不同撞擊位置下鋼套箱內(nèi)部構(gòu)件參與受力的順序也會(huì)不同. 例如對(duì)于第一類位置,撞深零時(shí)刻套箱內(nèi)部十字板即開(kāi)始發(fā)生變形;對(duì)于第二、第三類,首先參與受力的則是單獨(dú)一道橫肋或豎桁;第四類撞擊位置發(fā)生在內(nèi)部一點(diǎn)時(shí),直到船艏接觸到矩形四個(gè)角點(diǎn)才會(huì)有交叉十字構(gòu)件參與受力,這意味著與前三種撞擊位置相比,內(nèi)部構(gòu)件參與耗能的位移時(shí)刻稍滯后,但在后續(xù)很短的撞深范圍內(nèi),會(huì)有多個(gè)T形、L形構(gòu)件同時(shí)或依次發(fā)生變形,參與耗能. 各構(gòu)件參與碰撞的情形在2.4節(jié)展開(kāi)討論.

    根據(jù)上述分析,通過(guò)后續(xù)精細(xì)化有限元技術(shù)模擬得到如圖11所示撞擊力-撞深曲線,其中橫坐標(biāo)撞深范圍為0 ~ 0.2 m. 從圖11可看出,位置1下的撞擊力因外板變形區(qū)域較大和開(kāi)始時(shí)內(nèi)部十字板構(gòu)件的參與,初始撞擊力斜率最大;位置4在撞擊初始時(shí)刻套箱內(nèi)部構(gòu)件未參與受力,初始撞擊力斜率最小. 另外,由于鋼套箱內(nèi)部構(gòu)件排布密、間距小的特點(diǎn),撞擊點(diǎn)處相鄰的構(gòu)件能很快參與受力,因此當(dāng)撞深達(dá)到0.2 m時(shí),不同位置下鋼套箱撞擊力趨于一致,與船舶舷側(cè)結(jié)構(gòu)相比,后續(xù)撞擊力差別相對(duì)較小.

    2.4? ?鋼套箱內(nèi)部構(gòu)件參與碰撞情況

    對(duì)于上述4類鋼套箱被撞位置,對(duì)其后續(xù)的構(gòu)件碰撞情形討論在思路上是一致的. 因此,選取橫肋、縱桁交界處為被撞位置,展開(kāi)鋼套箱內(nèi)部各構(gòu)件參與碰撞情況的分析.

    以圖12為例,IJ為鋼套箱頂板,圖中所示為5道豎桁、2道橫肋、1道頂板所組成的縱橫艙壁結(jié)構(gòu). 在解析計(jì)算分析中,力學(xué)模型邊界條件為距離撞擊處最近且未參與碰撞的縱、橫肋,即當(dāng)撞擊位置為點(diǎn)O時(shí),阻力計(jì)算過(guò)程中模型邊界即為四邊形ACEG.

    對(duì)于帶球鼻艏船舶來(lái)說(shuō),假設(shè)撞擊點(diǎn)為O點(diǎn),一旦船艏與鋼套箱接觸,HD與BF所組成的X形十字板開(kāi)始變形,如圖12(a)所示. 當(dāng)船艏的進(jìn)一步行進(jìn)導(dǎo)致HD與BF組成的X形十字板構(gòu)件壓碎后,隨著位移的持續(xù)增大,球鼻下一步將會(huì)接觸到點(diǎn)H和點(diǎn)D. 此時(shí),鋼套箱內(nèi)參與受力的構(gòu)件為由HN、HA、HG構(gòu)成的T形構(gòu)件和與之對(duì)稱的由DK、DC、DE構(gòu)成的T形縱、橫梁組合結(jié)構(gòu),如圖12 (b)所示.

    對(duì)于楔形艏船舶,假設(shè)楔形體最先接觸到B點(diǎn),那么參與受力的為由BA、BC、BO組成的T形構(gòu)件,如圖12 (c)所示. 隨著楔形船的進(jìn)一步前進(jìn)接觸至點(diǎn)O,參與碰撞的為由HO、OD、OF組成的T形構(gòu)件. 當(dāng)碰撞至點(diǎn)A時(shí),由于AB在之前的碰撞中參與受力變形,喪失承載能力及邊界約束,此時(shí)參與碰撞的為由AI、AH組成的L形構(gòu)件,如圖12 (d)所示. 值得一提的是,一般情況下兩道豎桁的間距與兩道橫肋之間的距離有所不同,因此在計(jì)算T形構(gòu)件阻力時(shí),帶球艏船舶與楔形艏船舶撞擊場(chǎng)景下公式(11)(12)(13)中參數(shù)c的選取有所不同,需加以區(qū)分.

    2.5? ?構(gòu)件撕裂臨界撞深的確立

    對(duì)于鋼套箱豎向外板,在破裂前后其阻力大為不同,破裂前主要以膜拉伸及彎曲能量耗散的方式抵抗撞擊,而破裂后的外板阻力大大降低,視為撕裂變形. 因此,找到外板撕裂的臨界撞深是合理、準(zhǔn)確預(yù)測(cè)鋼套箱抗撞性的重要一步. 一般情況下認(rèn)為,外板變形增大至材料的斷裂應(yīng)變 ε0時(shí),視為外板撕裂.圖2為點(diǎn)荷載作用下外板變形示意,其中點(diǎn)畫(huà)線總長(zhǎng)度L0為變形前外板水平方向長(zhǎng)度,實(shí)線總長(zhǎng)度L為撞深Δ時(shí)外板變形后的長(zhǎng)度,板的拉伸應(yīng)變?nèi)缦拢?/p>

    當(dāng)ε = ε0時(shí)外板破裂,由cos α與撞深Δ的幾何關(guān)系可推出此時(shí)的臨界撞深為:

    在單軸拉伸試驗(yàn)中,低碳鋼的斷裂應(yīng)變一般為0.2~0.35. 考慮到尺寸效應(yīng)和材料缺陷,文獻(xiàn)[16]提及全尺寸船側(cè)模型中的斷裂應(yīng)變?nèi)?.05~0.10. 設(shè)計(jì)者可考慮實(shí)際情況,根據(jù)設(shè)計(jì)標(biāo)準(zhǔn)或其他研究中的建議值自行選取,本文取ε0 = 0.10.

    對(duì)于肋板、十字板等受力后先折疊再壓碎的構(gòu)件,每道折疊深度為5H0或2H1,若假設(shè)一次折疊后板件進(jìn)入撕裂破壞,則在參與構(gòu)件較多的情形下,尤其是鋼套箱內(nèi)部結(jié)構(gòu)阻力計(jì)算會(huì)明顯不足. 對(duì)于十字板構(gòu)件,通過(guò)式(15)可確定十字板在第N0道折疊中構(gòu)件破碎,整個(gè)構(gòu)件折疊深度為D=(N0-1)×2H1 .

    2.6? ?加勁肋的簡(jiǎn)化計(jì)算方法

    由于加勁肋不是主要受力構(gòu)件,其分擔(dān)阻力的能力有限但也不可忽視,因此采用等效厚度法[24]簡(jiǎn)化計(jì)算加勁肋的阻力,將加勁肋的橫截面積均勻分布給外板,便可得到等效厚度如下:

    式中:t為外板自身厚度;A1為加勁肋橫截面積;d為相鄰兩道加勁肋的間距.

    2.7? ?解析計(jì)算方法流程

    綜合上述研究?jī)?nèi)容,給出鋼套箱抗撞性能研究的解析計(jì)算方法,匯總至圖13.

    3? ?計(jì)算模型介紹

    通過(guò)選取一般典型構(gòu)造,建立了簡(jiǎn)化船艏模型、簡(jiǎn)化鋼套箱模型,以便于確定解析計(jì)算方法的普遍適用性. 選取精細(xì)化的5 000 DWT球艏船模型和椒江二橋主橋鋼套箱模型,采用解析計(jì)算方法和數(shù)值模擬結(jié)果對(duì)比,驗(yàn)證適用于鋼套箱的解析計(jì)算方法在實(shí)際工程中的準(zhǔn)確性. 本節(jié)將介紹各船艏和鋼套箱尺寸設(shè)計(jì)、內(nèi)外部構(gòu)造及其相應(yīng)有限元模型.

    3.1? ?不同外形船艏模型

    由文獻(xiàn)[11]提出的簡(jiǎn)化船模設(shè)計(jì)尺寸,分別建立船艏模型Ⅰ和Ⅱ. 如圖14所示,假定球艏船的船艏高度Hdeck = 8 m,上甲板楔形夾角的一半θ = 17°,首柱傾角φ = 63°,球鼻橢球部分RL = 2.4 m、RV = 1 m、RH =0.4 m. 楔形船的船艏高度Hdeck = 6.2 m,上甲板楔形夾角的一半θ = 30°,首柱傾角φ = 74°.

    船艏Ⅰ和Ⅱ的有限元模型材料為Q235鋼,密度ρ = 7 850 kg/m3,彈性模量E = 2.06 × 105 MPa,泊松比ν = 0.3,假定船艏外殼為剛體,采用殼單元建立,厚度10 mm.

    4.4? ?算例四

    由算例一和算例二對(duì)比看出球鼻艏和楔形艏撞擊下的鋼套箱阻力差別,因此算例四計(jì)算了楔形艏撞擊下鋼套箱Ⅱ的阻力,結(jié)果如圖20所示. 解析計(jì)算出的力-位移曲線在發(fā)展趨勢(shì)上與數(shù)值模擬結(jié)果高度吻合. 由于楔形艏撞擊下鋼套箱內(nèi)部參與受力的構(gòu)件數(shù)量多,解析計(jì)算中忽略了構(gòu)件間的相互作用,因此略低估了鋼套箱的耗能性能,總體而言誤差較小,且出于偏安全的考慮,計(jì)算結(jié)果可以接受.

    5? ?結(jié)? ?論

    本文提出了適用于橋梁固定式鋼套箱能力分析的解析計(jì)算方法,給出了詳細(xì)的計(jì)算流程,得到了在給定撞深下力的變化趨勢(shì)和套箱耗能能力. 該方法簡(jiǎn)單、快速,能有效地分析得到鋼套箱的抗撞能力曲線,采用精細(xì)有限元分析對(duì)計(jì)算結(jié)果進(jìn)行了驗(yàn)證,得到主要結(jié)論如下:

    1)采用提出的方法計(jì)算得到的撞擊力-撞深曲

    線與數(shù)值模擬方法得到的結(jié)果較好地吻合,表明了該方法的有效性、適用性. 將其運(yùn)用于橋梁防撞鋼套箱初步設(shè)計(jì)中,將有效提高設(shè)計(jì)效率.

    2)對(duì)于帶球艏和楔形艏兩種類型艏部的船舶,

    其撞擊下的鋼套箱阻力存在顯著差別. 由于楔形艏船對(duì)鋼套箱的破壞程度更高,結(jié)果往往大于帶球艏船撞擊下的阻力.

    3)針對(duì)橋梁防撞鋼套箱整體尺寸小、內(nèi)部構(gòu)件密的構(gòu)造特點(diǎn),對(duì)不同撞擊位置下,初始撞深內(nèi)板件的變形特征做了較為詳細(xì)的分析與討論,分析假定的變形特征與數(shù)值模擬結(jié)果較好地吻合.

    4)解析計(jì)算了帶有一定弧度的鋼套箱碰撞阻

    力,并利用數(shù)值模擬方法對(duì)比結(jié)果,驗(yàn)證了鋼套箱解析計(jì)算方法在該類鋼套箱分析中的適用性.

    本文重點(diǎn)建立了解析計(jì)算方法的基本思路,僅進(jìn)行了有限情況的驗(yàn)證,后續(xù)應(yīng)結(jié)合不同類型船舶、不同套箱(如浮式套箱、組合結(jié)構(gòu)套箱)等進(jìn)一步驗(yàn)證方法的有效性,同時(shí),如何運(yùn)用其進(jìn)行設(shè)計(jì),有待進(jìn)一步深入研究.

    參考文獻(xiàn)

    [1]? ? LARSEN O D. Ship collision with bridges:the interaction between vessel traffic and bridge structures [M]. Geneva,Switzerland:IABSE,1993:1—6.

    [2]? ? FAN W,GUO W,SUN Y,et al. Experimental and numerical investigations of a novel steel-UHPFRC composite fender for bridge protection in vessel collisions [J]. Ocean Engineering,2018,165:1—21.

    [3]? ? 項(xiàng)海帆,范立礎(chǔ),王君杰. 船撞橋設(shè)計(jì)理論的現(xiàn)狀與需進(jìn)一步研究的問(wèn)題[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,30(4):386—392.

    XIANG H F,F(xiàn)AN L C,WANG J J. State of art of ship collision design for bridges and future research[J]. Journal of Tongji University (Natural Science),2002,30(4):386—392. (In Chinese)

    [4]? ? JTG D60—2015? 公路橋涵設(shè)計(jì)通用規(guī)范[S]. 北京:人民交通出版社,2015:66—67.

    JTG D60—2015 General code for design of highway bridge and culverts? [S]. Beijing:China Communication Press,2015:66—67. (In Chinese)

    [5]? ? 樊偉. 船撞下橋梁結(jié)構(gòu)動(dòng)力需求及樁承結(jié)構(gòu)防撞能力分析方法[D]. 上海:同濟(jì)大學(xué)土木工程學(xué)院,2012:18—24.

    FAN W. Dynamic demand of bridge structures and capacity of pile-supported protection structures under vessel impacts [D]. Shanghai: College of Civil Engineering,Tongji University,2012:18—24. (In Chinese)

    [6]? ? 游慶仲,董學(xué)武,吳壽昌. 蘇通大橋基礎(chǔ)工程的挑戰(zhàn)與創(chuàng)新[J].中國(guó)工程科學(xué),2007,9(6):22—26.

    YOU Q Z,DONG X W,WU S C. Challenge and innovation in Sutong bridge foundation project [J]. Engineering Science,2007,9(6):22—26. (In Chinese)

    [7]? ? 陳衛(wèi)國(guó),唐衡. 金塘大橋索塔墩1 600 t防撞鋼套箱安裝施工技術(shù)[J]. 公路,2009(1):133—138.

    CHEN W G,TANG H. Installation and construction technology of

    1 600 t anti-collision steel box for cable tower pier of Jintang bridge [J]. Highway,2009(1):133—138. (In Chinese)

    [8]? ? 陳煒,肖波,高宗余. 武漢天興洲公鐵兩用長(zhǎng)江大橋1號(hào)墩防撞設(shè)計(jì)與分析[J]. 世界橋梁,2006(2):78—80.

    CHEN W,XIAO B,GAO Z Y. Design and analysis of ship collision protection for pier no.1 of Wuhan Tianxingzhou Changjiang River road-cum-rail bridge [J]. World Bridges, 2006(2): 78—80. (In Chinese)

    [9]? ? FAN W,YUAN W,CHEN B. Steel fender limitations and improvements for bridge protection in ship collisions [J]. Journal of Bridge Engineering,2015,20(12):06015004.

    [10]? 袁濤,王小鋒. 杭州灣跨海大橋鋼套箱施工技術(shù)總結(jié)[J]. 公路,2006(9):76—83.

    YUAN T,WANG X F. Summary of construction technology of steel sleeve box of Hangzhou bay sea-crossing bridge [J]. Highway,2006(9):76—83. (In Chinese)

    [11]? ZHANG S. The mechanics of ship collisions [D]. Lyngby:Department of Naval Architecture and Offshore Engineering,Technical University of Denmark,1999:69—122.

    [12]? 王娟,錢(qián)江,周德源. 城市橋梁下部結(jié)構(gòu)抗重型車(chē)輛撞擊的數(shù)值仿真分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,43(7):88—95.

    WANG J,QIAN J,ZHOU D Y. Numerical simulation of urban bridge substructures impacted by heavy vehicles [J]. Journal of Hunan University (Natural Sciences), 2016, 43(7): 88—95. (In Chinese)

    [13]? 許薛軍,單成林. 夾層結(jié)構(gòu)曲面環(huán)形浮式橋墩防撞套箱碰撞分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,42(3):106—111.

    XU X J,SHAN C L. Impact analysis of the bridge pier anti-collision floating box sets made by sandwich structure with curved-shaped [J]. Journal of Hunan University (Natural Sciences),2015,42(3):106—111. (In Chinese)

    [14]? JIANG H,CHORZEPA M G. Case study:evaluation of a floating steel fender system for bridge pier protection against vessel collision [J]. Journal of Bridge Engineering,2016,21(11):05016008.

    [15]? WANG J J,SONG Y C,WANG W,et al. Calibrations of numerical models by experimental impact tests using scaled steel boxes [J]. Engineering Structures,2018,173:481—494.

    [16]? HARIS S,AMDAHL J. An analytical model to assess a ship side during a collision [J]. Ships and Offshore Structures,2012,7(4):431—448.

    [17]? GAO Z,HU Z,WANG G,et al. An analytical method of predicting the response of FPSO side structures to head-on collision [J].

    Ocean Engineering,2014,87:121—135.

    [18]? WANG G,OHTSUBO H,ARITA K. Large deflection of a rigid—plastic circular plate pressed by a sphere [J]. Journal of Applied Mechanics,1998,65(2):533—535.

    [19]? 孫斌,胡志強(qiáng),王晉. 楔形船艏撞擊舷側(cè)外板的結(jié)構(gòu)響應(yīng)分析[J]. 振動(dòng)與沖擊,2016,35(23):46—50.

    SUN B,HU Z Q,WANG J. Structural response analysis for a ship side plating impacted by raked bow [J]. Journal of Vibration and Shock,2016,35(23):46—50. (In Chinese)

    [20]? OHTSUBO H,WANG G. An upper-bound solution to the problem of plate tearing [J]. Journal of Marine Science and Technology,1995,1(1):46—51.

    [21]? WANG G,ARITA K,LIU D. Behavior of a double hull in a variety of stranding or collision scenarios [J]. Marine Structures,2000,13(3):147—187.

    [22]? LIU B,SOARES G C. Simplified analytical method for evaluating web girder crushing during ship collision and grounding [J]. Marine Structures,2015,42:71—94.

    [23]? HONG L. Simplified analysis and design of ships subjected to collision and grounding [D]. Trondheim,Norway:Norwegian University of Science and Technology,2008:17—19.

    [24]? PAIK J K,CHUNG J Y,CHUN M S. On quasi-static crushing of a stiffened square tube [J]. Journal of Ship Research,1996,40(3):258—267.

    [25]? FAN W,YUAN W C. Ship bow force-deformation curves for ship-impact demand of bridges considering effect of pile-cap depth [J]. Shock and Vibration,2014,2014(2):1—19.

    [26]? FAN W,YUAN W C. Numerical simulation and analytical modeling of pile-supported structures subjected to ship collisions including soil-structure interaction[J]. Ocean Engineering,2014,91:11—27.

    猜你喜歡
    有限元分析
    對(duì)于Pro/mechanica的數(shù)控車(chē)刀的有限元分析
    S型壓力傳感器形變的有限元分析
    橫向穩(wěn)定桿的側(cè)傾角剛度仿真與試驗(yàn)研究
    多自由度升降翻轉(zhuǎn)舞臺(tái)的運(yùn)動(dòng)分析
    演藝科技(2016年11期)2016-12-24 20:18:29
    高精度細(xì)深孔的鏜孔加工工藝分析
    基于LS—DYNA的某汽車(chē)乘客安全氣囊支架焊接強(qiáng)度的分析
    基于有限元分析談框剪結(jié)構(gòu)力學(xué)性能的影響因素
    自錨式懸索橋鋼箱梁頂推施工階段結(jié)構(gòu)分析
    隨機(jī)振動(dòng)載荷下發(fā)射裝置尾罩疲勞壽命分析
    航空兵器(2016年4期)2016-11-28 21:54:01
    有限元分析帶溝槽平封頭的應(yīng)力集中
    国产亚洲精品第一综合不卡 | 男女高潮啪啪啪动态图| 下体分泌物呈黄色| 18禁动态无遮挡网站| 欧美精品高潮呻吟av久久| 伦精品一区二区三区| 一级毛片黄色毛片免费观看视频| av片东京热男人的天堂| 精品卡一卡二卡四卡免费| 亚洲欧美日韩另类电影网站| av视频免费观看在线观看| 久热久热在线精品观看| 不卡视频在线观看欧美| 香蕉国产在线看| 免费av中文字幕在线| 亚洲色图综合在线观看| 一级黄片播放器| 国产xxxxx性猛交| 性高湖久久久久久久久免费观看| 制服诱惑二区| 久久ye,这里只有精品| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 亚洲精品美女久久久久99蜜臀 | 久久精品久久久久久久性| 精品视频人人做人人爽| 免费观看在线日韩| 久久久久久伊人网av| 亚洲高清免费不卡视频| 十八禁高潮呻吟视频| av有码第一页| 久久精品国产自在天天线| 美女xxoo啪啪120秒动态图| 在线 av 中文字幕| 亚洲精品久久成人aⅴ小说| 九九在线视频观看精品| 最新中文字幕久久久久| 午夜免费鲁丝| 久久久久久人人人人人| 久久av网站| 日韩av不卡免费在线播放| 一边亲一边摸免费视频| 97精品久久久久久久久久精品| 美女中出高潮动态图| 永久网站在线| 国产亚洲午夜精品一区二区久久| av女优亚洲男人天堂| 国产色婷婷99| 亚洲av欧美aⅴ国产| 国产成人精品婷婷| 国产精品99久久99久久久不卡 | 日日啪夜夜爽| 99国产综合亚洲精品| 99香蕉大伊视频| 国产成人av激情在线播放| 一边摸一边做爽爽视频免费| 成人综合一区亚洲| 天天躁夜夜躁狠狠久久av| 777米奇影视久久| videossex国产| 香蕉精品网在线| 午夜福利,免费看| 亚洲一区二区三区欧美精品| 亚洲欧洲精品一区二区精品久久久 | 乱码一卡2卡4卡精品| 久久精品人人爽人人爽视色| 99热国产这里只有精品6| 人人澡人人妻人| 亚洲人与动物交配视频| 在线观看免费视频网站a站| 大话2 男鬼变身卡| 2018国产大陆天天弄谢| 熟女人妻精品中文字幕| 80岁老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| 欧美性感艳星| 中文字幕制服av| 亚洲国产日韩一区二区| 极品人妻少妇av视频| 在线观看免费高清a一片| 三上悠亚av全集在线观看| 伦精品一区二区三区| 中文欧美无线码| 精品一区在线观看国产| 亚洲av综合色区一区| 国产毛片在线视频| 99热国产这里只有精品6| 日日爽夜夜爽网站| 一个人免费看片子| 日韩一本色道免费dvd| 高清av免费在线| 久久久久网色| 亚洲精品一二三| 边亲边吃奶的免费视频| 国产男人的电影天堂91| 九色亚洲精品在线播放| 18禁国产床啪视频网站| 国产欧美日韩综合在线一区二区| 亚洲高清免费不卡视频| 国产精品久久久av美女十八| 大香蕉久久成人网| 18禁在线无遮挡免费观看视频| 国产乱人偷精品视频| av.在线天堂| 国产一区亚洲一区在线观看| 欧美精品亚洲一区二区| 欧美激情极品国产一区二区三区 | 交换朋友夫妻互换小说| 免费看光身美女| 肉色欧美久久久久久久蜜桃| 日韩三级伦理在线观看| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 欧美日韩成人在线一区二区| 亚洲欧洲精品一区二区精品久久久 | 免费久久久久久久精品成人欧美视频 | 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产av成人精品| 精品亚洲乱码少妇综合久久| 久久久久久久久久人人人人人人| 国产成人午夜福利电影在线观看| 亚洲激情五月婷婷啪啪| 国产男人的电影天堂91| 亚洲四区av| 80岁老熟妇乱子伦牲交| 这个男人来自地球电影免费观看 | 日本-黄色视频高清免费观看| 视频区图区小说| 黄色怎么调成土黄色| 青春草国产在线视频| 国产不卡av网站在线观看| 天天操日日干夜夜撸| 伦理电影免费视频| 亚洲成人一二三区av| 亚洲,欧美精品.| 日本午夜av视频| 日本av手机在线免费观看| 欧美丝袜亚洲另类| 亚洲精品美女久久av网站| 精品国产一区二区久久| 好男人视频免费观看在线| 少妇熟女欧美另类| 黑丝袜美女国产一区| 草草在线视频免费看| 免费观看性生交大片5| 高清在线视频一区二区三区| 波野结衣二区三区在线| 久久久久视频综合| 亚洲精品第二区| 美女中出高潮动态图| 成人毛片a级毛片在线播放| 日韩av免费高清视频| 最后的刺客免费高清国语| 美女中出高潮动态图| 国产成人精品久久久久久| 亚洲精品色激情综合| 最近的中文字幕免费完整| 午夜视频国产福利| 亚洲精品乱码久久久久久按摩| 欧美另类一区| 亚洲成人一二三区av| 国产av精品麻豆| 女人精品久久久久毛片| 老司机亚洲免费影院| 免费少妇av软件| 高清黄色对白视频在线免费看| 在线观看三级黄色| 2021少妇久久久久久久久久久| 久久精品国产亚洲av天美| 波野结衣二区三区在线| 亚洲精品色激情综合| 国产精品一二三区在线看| 国产精品人妻久久久久久| freevideosex欧美| 久久精品久久精品一区二区三区| 久久人人97超碰香蕉20202| 九草在线视频观看| 精品一区二区免费观看| 亚洲成人一二三区av| 国产伦理片在线播放av一区| 免费人妻精品一区二区三区视频| 在线观看免费高清a一片| 亚洲国产成人一精品久久久| 日本免费在线观看一区| 国产精品女同一区二区软件| 在线亚洲精品国产二区图片欧美| 夫妻性生交免费视频一级片| xxx大片免费视频| 韩国av在线不卡| 精品人妻在线不人妻| 汤姆久久久久久久影院中文字幕| 啦啦啦视频在线资源免费观看| 纯流量卡能插随身wifi吗| 制服丝袜香蕉在线| 久久国产亚洲av麻豆专区| 九九爱精品视频在线观看| 男的添女的下面高潮视频| 丰满迷人的少妇在线观看| 亚洲精品一区蜜桃| 国产成人精品无人区| 男女国产视频网站| 婷婷色av中文字幕| 天堂中文最新版在线下载| 777米奇影视久久| 熟女人妻精品中文字幕| 夫妻性生交免费视频一级片| 夫妻午夜视频| 精品人妻一区二区三区麻豆| 在线天堂中文资源库| 亚洲国产av影院在线观看| 婷婷色综合www| 深夜精品福利| 我的女老师完整版在线观看| 水蜜桃什么品种好| 欧美日韩一区二区视频在线观看视频在线| 久久午夜综合久久蜜桃| 久久国产精品男人的天堂亚洲 | 午夜免费鲁丝| 欧美97在线视频| 亚洲激情五月婷婷啪啪| 中文字幕av电影在线播放| 久久狼人影院| av在线观看视频网站免费| 一区二区三区四区激情视频| 涩涩av久久男人的天堂| 男的添女的下面高潮视频| 亚洲国产精品一区三区| 久久人人97超碰香蕉20202| 日韩av在线免费看完整版不卡| 久久综合国产亚洲精品| 精品亚洲乱码少妇综合久久| 亚洲av国产av综合av卡| 在现免费观看毛片| 亚洲精品一二三| 日韩伦理黄色片| 国语对白做爰xxxⅹ性视频网站| tube8黄色片| 欧美日韩综合久久久久久| 爱豆传媒免费全集在线观看| 国产日韩欧美视频二区| 国产免费一级a男人的天堂| 色哟哟·www| 人妻系列 视频| 尾随美女入室| 女人精品久久久久毛片| 一边摸一边做爽爽视频免费| av.在线天堂| 欧美国产精品一级二级三级| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| www.色视频.com| 老司机影院成人| 女人久久www免费人成看片| 色婷婷av一区二区三区视频| 美国免费a级毛片| 在线观看免费视频网站a站| 男女啪啪激烈高潮av片| 国产av国产精品国产| 久久 成人 亚洲| 日日啪夜夜爽| 91国产中文字幕| 国产精品久久久久久精品电影小说| 亚洲人成网站在线观看播放| 国产一区二区在线观看日韩| 在现免费观看毛片| 日本黄色日本黄色录像| 欧美性感艳星| videos熟女内射| 亚洲第一区二区三区不卡| 午夜久久久在线观看| 性色av一级| 国产精品.久久久| 狂野欧美激情性xxxx在线观看| 久久人人97超碰香蕉20202| 永久网站在线| 啦啦啦中文免费视频观看日本| 熟妇人妻不卡中文字幕| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 免费人成在线观看视频色| av在线观看视频网站免费| 少妇猛男粗大的猛烈进出视频| 欧美精品高潮呻吟av久久| 1024视频免费在线观看| av女优亚洲男人天堂| 亚洲国产日韩一区二区| 国产男女内射视频| av线在线观看网站| 汤姆久久久久久久影院中文字幕| 大香蕉97超碰在线| 在线精品无人区一区二区三| 欧美精品国产亚洲| 午夜精品国产一区二区电影| 国产成人精品久久久久久| 久久午夜福利片| 国产又爽黄色视频| 免费观看a级毛片全部| 久久精品国产综合久久久 | 视频在线观看一区二区三区| 高清欧美精品videossex| 亚洲欧美成人精品一区二区| 久久久久久伊人网av| 丝袜脚勾引网站| 日韩成人伦理影院| 亚洲av男天堂| 热re99久久精品国产66热6| 亚洲经典国产精华液单| 男女高潮啪啪啪动态图| 亚洲av在线观看美女高潮| 免费观看a级毛片全部| 女性被躁到高潮视频| 国产成人一区二区在线| 自拍欧美九色日韩亚洲蝌蚪91| 少妇熟女欧美另类| 国产av码专区亚洲av| 高清不卡的av网站| 久久婷婷青草| 草草在线视频免费看| 免费黄色在线免费观看| 中文天堂在线官网| 啦啦啦中文免费视频观看日本| 国产成人精品一,二区| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 免费看光身美女| 桃花免费在线播放| 国产成人欧美| 国产永久视频网站| 美女内射精品一级片tv| 免费在线观看完整版高清| 国产男女内射视频| 日日啪夜夜爽| 国产 一区精品| 一本大道久久a久久精品| 新久久久久国产一级毛片| 草草在线视频免费看| 51国产日韩欧美| 国产高清不卡午夜福利| 人人妻人人添人人爽欧美一区卜| 久久狼人影院| 午夜激情久久久久久久| 韩国精品一区二区三区 | 国产精品久久久久久久久免| 精品人妻偷拍中文字幕| 97人妻天天添夜夜摸| 老司机影院毛片| 亚洲欧美日韩卡通动漫| 精品一品国产午夜福利视频| 美女国产视频在线观看| 水蜜桃什么品种好| 这个男人来自地球电影免费观看 | 亚洲内射少妇av| 久久久久久久大尺度免费视频| 成人黄色视频免费在线看| 国产 精品1| 久久久国产精品麻豆| 久久精品国产自在天天线| 最近2019中文字幕mv第一页| 在线看a的网站| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av涩爱| 亚洲欧洲精品一区二区精品久久久 | 免费观看在线日韩| 亚洲欧美清纯卡通| 丰满饥渴人妻一区二区三| 免费av不卡在线播放| 9191精品国产免费久久| 亚洲精品美女久久久久99蜜臀 | 婷婷色av中文字幕| 黑丝袜美女国产一区| 国产亚洲欧美精品永久| 久久精品国产a三级三级三级| 大香蕉久久成人网| 成人二区视频| 国产极品天堂在线| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看 | 天天躁夜夜躁狠狠久久av| 久久国内精品自在自线图片| 色94色欧美一区二区| 美女内射精品一级片tv| 精品酒店卫生间| 国产精品欧美亚洲77777| 日韩一区二区视频免费看| 色婷婷av一区二区三区视频| 丁香六月天网| 草草在线视频免费看| 免费看av在线观看网站| av在线观看视频网站免费| 亚洲欧美中文字幕日韩二区| av网站免费在线观看视频| 咕卡用的链子| 久久av网站| 精品久久国产蜜桃| 天天影视国产精品| 国产亚洲欧美精品永久| 精品国产一区二区三区四区第35| 色哟哟·www| 国产精品久久久av美女十八| 视频区图区小说| 伊人久久国产一区二区| 草草在线视频免费看| 男女午夜视频在线观看 | 成年人午夜在线观看视频| videosex国产| 91精品国产国语对白视频| 看十八女毛片水多多多| 国产av国产精品国产| 欧美变态另类bdsm刘玥| 久久ye,这里只有精品| 日韩精品免费视频一区二区三区 | 中文字幕av电影在线播放| 亚洲av日韩在线播放| av不卡在线播放| 亚洲经典国产精华液单| 美女福利国产在线| 夜夜爽夜夜爽视频| 欧美国产精品一级二级三级| 久久这里有精品视频免费| 久久久久视频综合| 欧美日韩亚洲高清精品| 国产精品成人在线| 韩国av在线不卡| 男女免费视频国产| 人妻系列 视频| 日韩一本色道免费dvd| 丝袜脚勾引网站| 18在线观看网站| 国产精品久久久久久精品电影小说| 黑丝袜美女国产一区| 美女福利国产在线| 亚洲国产精品一区二区三区在线| 最近手机中文字幕大全| 大陆偷拍与自拍| 欧美性感艳星| 亚洲成色77777| 国产免费视频播放在线视频| 不卡视频在线观看欧美| 久久影院123| 99re6热这里在线精品视频| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 日韩中文字幕视频在线看片| 97在线人人人人妻| 搡女人真爽免费视频火全软件| a级毛色黄片| 日韩伦理黄色片| 妹子高潮喷水视频| 久久99精品国语久久久| 亚洲国产欧美日韩在线播放| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品电影小说| 你懂的网址亚洲精品在线观看| 啦啦啦中文免费视频观看日本| 国产成人午夜福利电影在线观看| 中国国产av一级| 久久久久精品人妻al黑| 久久精品熟女亚洲av麻豆精品| 国产精品无大码| 晚上一个人看的免费电影| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站| 欧美bdsm另类| 黑人巨大精品欧美一区二区蜜桃 | 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到 | 国产男人的电影天堂91| 国产精品久久久av美女十八| 欧美国产精品一级二级三级| 午夜免费观看性视频| 国产激情久久老熟女| 极品少妇高潮喷水抽搐| 亚洲欧美成人综合另类久久久| 久热久热在线精品观看| 少妇被粗大猛烈的视频| 一个人免费看片子| 日韩三级伦理在线观看| 巨乳人妻的诱惑在线观看| 制服人妻中文乱码| 人妻一区二区av| 狂野欧美激情性xxxx在线观看| 熟妇人妻不卡中文字幕| 五月玫瑰六月丁香| 免费少妇av软件| 爱豆传媒免费全集在线观看| 午夜免费男女啪啪视频观看| 欧美国产精品一级二级三级| 18禁国产床啪视频网站| av国产久精品久网站免费入址| 日韩,欧美,国产一区二区三区| 黄色怎么调成土黄色| 亚洲国产日韩一区二区| 亚洲av免费高清在线观看| 日本免费在线观看一区| 春色校园在线视频观看| 最近中文字幕2019免费版| 咕卡用的链子| av.在线天堂| 日韩免费高清中文字幕av| 亚洲av中文av极速乱| 欧美成人午夜精品| 在线观看免费日韩欧美大片| 免费黄网站久久成人精品| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 久久这里只有精品19| 国产免费又黄又爽又色| 国产激情久久老熟女| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 日本黄色日本黄色录像| 美女内射精品一级片tv| 美女脱内裤让男人舔精品视频| 免费黄色在线免费观看| av在线观看视频网站免费| 妹子高潮喷水视频| 国产欧美另类精品又又久久亚洲欧美| 国产日韩欧美在线精品| 一边亲一边摸免费视频| 久久这里有精品视频免费| 免费观看无遮挡的男女| 国产白丝娇喘喷水9色精品| videossex国产| 日韩熟女老妇一区二区性免费视频| 久久青草综合色| 九草在线视频观看| 亚洲国产最新在线播放| freevideosex欧美| 欧美丝袜亚洲另类| av福利片在线| 午夜福利视频在线观看免费| 自拍欧美九色日韩亚洲蝌蚪91| 免费人成在线观看视频色| 亚洲av福利一区| 在线精品无人区一区二区三| 国精品久久久久久国模美| 成人免费观看视频高清| 97人妻天天添夜夜摸| 18在线观看网站| 亚洲av综合色区一区| av又黄又爽大尺度在线免费看| 国产色婷婷99| 欧美少妇被猛烈插入视频| av视频免费观看在线观看| 久久韩国三级中文字幕| 亚洲成人av在线免费| 成年人午夜在线观看视频| 看非洲黑人一级黄片| 久久精品国产a三级三级三级| 久久久久视频综合| 国产一区二区三区综合在线观看 | 亚洲在久久综合| 午夜免费鲁丝| 国产福利在线免费观看视频| 一边摸一边做爽爽视频免费| av免费观看日本| 在线观看人妻少妇| 国产亚洲午夜精品一区二区久久| 欧美老熟妇乱子伦牲交| 一本大道久久a久久精品| 美女脱内裤让男人舔精品视频| av女优亚洲男人天堂| 亚洲av欧美aⅴ国产| 青春草国产在线视频| 久久精品久久久久久噜噜老黄| 80岁老熟妇乱子伦牲交| xxxhd国产人妻xxx| 美女脱内裤让男人舔精品视频| 欧美3d第一页| 婷婷色综合大香蕉| 欧美精品亚洲一区二区| 亚洲综合色惰| 成年动漫av网址| 国产视频首页在线观看| 中文字幕人妻丝袜制服| 51国产日韩欧美| 视频在线观看一区二区三区| av黄色大香蕉| 2018国产大陆天天弄谢| 国产精品人妻久久久久久| 色网站视频免费| 亚洲av免费高清在线观看| 日韩伦理黄色片| 蜜桃在线观看..| 中文字幕精品免费在线观看视频 | 天天躁夜夜躁狠狠躁躁| 男女高潮啪啪啪动态图| 国产精品久久久久久久电影| 国产一级毛片在线| 国产激情久久老熟女| 亚洲人成77777在线视频| 啦啦啦啦在线视频资源| 日本猛色少妇xxxxx猛交久久| 久久国产精品大桥未久av| 激情视频va一区二区三区| 在线观看三级黄色| 亚洲经典国产精华液单| 啦啦啦啦在线视频资源| 国产男女内射视频| 欧美日韩综合久久久久久| 婷婷成人精品国产| 高清av免费在线| 99精国产麻豆久久婷婷| 国产精品秋霞免费鲁丝片| 国产精品人妻久久久久久| 国产精品熟女久久久久浪| 亚洲av成人精品一二三区| 亚洲性久久影院| 曰老女人黄片| 亚洲国产精品专区欧美| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 九色亚洲精品在线播放| 日韩大片免费观看网站| 少妇高潮的动态图| 日韩一区二区三区影片|