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

    含初始撓度加筋板的極限承載能力分析

    2017-06-19 19:18:59劉春正牟金磊
    艦船科學(xué)技術(shù) 2017年5期
    關(guān)鍵詞:半波筋板加強(qiáng)筋

    劉春正,吳 梵,牟金磊

    (海軍工程大學(xué) 艦船工程系,湖北 武漢 430033)

    含初始撓度加筋板的極限承載能力分析

    劉春正,吳 梵,牟金磊

    (海軍工程大學(xué) 艦船工程系,湖北 武漢 430033)

    艦船長期服役,甲板結(jié)構(gòu)易產(chǎn)生初始撓度變形,這會(huì)對(duì)甲板承載能力帶來不利影響。加筋板作為船體甲板結(jié)構(gòu)的主要構(gòu)成單元,研究初始撓度變形對(duì)其極限承載力的影響具有重要意義。為了確定初始撓度變形對(duì)加筋板極限承載力的影響作用,根據(jù)實(shí)際情況假設(shè)初始撓度為雙三角級(jí)數(shù)形式,利用 Ansys 計(jì)算分析了整體初始撓度的幅值與半波數(shù)對(duì)極限載荷的影響和典型位置的應(yīng)力特性,并得到初始撓度對(duì)加筋板極限載荷的影響因子計(jì)算方法。計(jì)算結(jié)果分析表明,隨著初始撓度的幅值和半波數(shù)的增加,加筋板極限承載力逐漸減??;對(duì)于含有某種初始撓度的加筋板,其影響因子主要受加筋板的長寬比、厚度和加強(qiáng)筋間距等因素的影響。

    初始撓度;加筋板;極限承載力;應(yīng)力特性;影響因子

    0 引 言

    隨著我國造船業(yè)迅速發(fā)展,許多新式艦船紛紛建造下水,對(duì)長期服役的老舊艦船船體結(jié)構(gòu)技術(shù)狀態(tài)的鑒定工作也逐漸開啟,對(duì)艦船結(jié)構(gòu)強(qiáng)度的校核是技術(shù)狀態(tài)鑒定的重要組成部分。甲板作為艦船的重要結(jié)構(gòu)之一,對(duì)甲板結(jié)構(gòu)極限承載能力的分析是對(duì)艦船結(jié)構(gòu)強(qiáng)度評(píng)估的重要組成內(nèi)容。

    加筋板被廣泛應(yīng)用于許多工程結(jié)構(gòu)中,它大大提高了板的承載能力[1]。加筋板結(jié)構(gòu)為船體甲板結(jié)構(gòu)的基本構(gòu)件,其承載能力直接影響整體結(jié)構(gòu)的極限承載能力。艦船在中垂?fàn)顟B(tài)時(shí)甲板板架受總縱彎曲軸向壓力,甲板板架是距離中和軸最遠(yuǎn)的主船體受壓結(jié)構(gòu),所受的軸向壓力最大,只有保證甲板的承載能力,才能保證艦船整體結(jié)構(gòu)的安全性。加筋板作為船體甲板的主要結(jié)構(gòu)單元,研究其受壓狀態(tài)下的極限承載能力十分必要。一般地,受到軸向壓力作用的加筋板的可能破壞模式包括:總體屈曲后的總體破壞;加強(qiáng)筋之間的板格局部屈曲后引起的總體破壞;加強(qiáng)筋(無轉(zhuǎn)動(dòng))局部屈曲或屈服后引起的總體破壞;加強(qiáng)筋彎扭屈曲(側(cè)屈),整個(gè)加筋板隨之喪失承載能力;加筋板的整體屈服[2]。

    水面艦船在長期使用過程中,船體加筋板結(jié)構(gòu)由于各種載荷的作用,會(huì)產(chǎn)生撓度變形,對(duì)后續(xù)載荷而言,這種撓度變形在加載之前就已經(jīng)產(chǎn)生,故稱為初始撓度。初始撓度對(duì)加筋板的極限承載能力會(huì)產(chǎn)生不利影響。因此,在加筋板的極限承載能力分析中,考慮初始變形等因素的影響,具有重要的工程意義。

    國內(nèi)外考慮初始撓度影響下,對(duì)加筋板極限強(qiáng)度的研究比較多。Ueda 等[3]推導(dǎo)了受雙軸壓應(yīng)力、剪應(yīng)力,并考慮初始變形和焊接殘余應(yīng)力情況下,板有效寬度的計(jì)算公式;Paik 等[4]考慮焊接導(dǎo)致的初始變形,對(duì)加筋板進(jìn)行了非線性屈曲分析;萬育龍等[5]分析了考慮多種應(yīng)力組合模式及局部、整體初始缺陷組合模式下的加筋板的極限強(qiáng)度,歸納出了適用于 CSR-OT加筋板屈曲和極限強(qiáng)度計(jì)算的替代方法;羅剛等[6]研究了不同外力工況下加筋板極限強(qiáng)度的響應(yīng)特點(diǎn),過程中考慮了加筋板初始撓度的作用。這些研究雖然都考慮了初始撓度的影響作用,但都只是考慮了焊接產(chǎn)生的初始撓度,并且在有限元計(jì)算中,直接取一階失穩(wěn)模態(tài)作為加筋板的初始撓度模式,沒有考慮初始撓度的具體形態(tài)對(duì)加筋板承載能力的影響。本文從初始撓度的具體形態(tài)出發(fā),計(jì)算含初始撓度加筋板在加強(qiáng)筋方向單向受壓情況下的極限強(qiáng)度;考慮初始撓度幅值和半波數(shù)對(duì)加筋板極限承載力的影響,提出了計(jì)算公式修正系數(shù);并研究分析了含初始撓度加筋板典型位置剖面的應(yīng)力分布特性以及初始撓度對(duì)加筋板破環(huán)模式的影響。

    1 含初始撓度加筋板的穩(wěn)定性分析

    1.1 模型建立

    本文根據(jù)某艦船的甲板結(jié)構(gòu),選取甲板典型位置處的某一加筋板作為研究對(duì)象。該加筋板具有 6 條縱向加強(qiáng)筋,板尺寸 a × b 為 1 500 mm × 2 800 mm,加強(qiáng)筋間距 s = 400 mm,加強(qiáng)筋均為角鋼,角鋼尺寸為。

    研究中,利用有限元分析軟件 Ansys,根據(jù)加筋板的具體尺寸,建立合適的有限元模型,如圖 1 所示。采用 shell181 單元,單元長度在 20 ~ 30 mm 之間,計(jì)算中假定板的材料是理想彈塑性,彈性模量 E = 2.1 × 105MPa,泊松比 ν = 0.3。

    加筋板的邊界條件定義為四邊簡(jiǎn)支,具體定義方法可以參照相關(guān)文獻(xiàn)[7]:

    1)加筋板的四周邊界限制法向位移 Uz= 0,即為z 方向簡(jiǎn)支;

    2)長度方向兩端 AB,CD 邊(x = 0和 x = a 處)限制轉(zhuǎn)角位移 Rx= Rz= 0,即可沿 y 方向自由轉(zhuǎn)動(dòng);寬度方向兩端 AD,BC 邊(y = 0 和 y = b 處)限制轉(zhuǎn)角位移 Ry= Rz= 0,即可沿 x 方向自由轉(zhuǎn)動(dòng);

    3)AB,CD 邊線保持沿 x 方向的直邊約束,AD,BC 邊線保持沿 y 方向的直邊約束;

    4)AB、CD 邊線中點(diǎn)限制 y 方向位移,即 Uy= 0,AD、BC 邊線中點(diǎn)限制 x 方向位移,即 Ux= 0。

    加筋板的外載荷施加方式為在 AB,CD 兩端截面同時(shí)施加軸向壓力為 1 MPa 的集中力載荷,在有限元軟件中轉(zhuǎn)換施加到截面節(jié)點(diǎn)上。

    1.2 初始撓度設(shè)計(jì)

    初始撓度這種結(jié)構(gòu)上的缺陷可以分為整體和局部2 種形態(tài)。整體初始撓度是指板和加強(qiáng)筋同時(shí)產(chǎn)生撓度變形;局部初始撓度是指僅板產(chǎn)生的撓度變形。初始撓度形態(tài)主要包括兩方面因素,即幅值和半波數(shù)。在較長的跨度范圍內(nèi),初始撓度是由多個(gè)半波所組成的曲面,而在較短的范圍內(nèi)(2 個(gè)加強(qiáng)筋之間)一般只有一個(gè)半波。船體加筋板初始撓度形態(tài)可以看成由多種半波形式疊加而成,可用如下公式表示:

    也有用如下公式來表示焊接作用導(dǎo)致加筋板產(chǎn)生的撓度變形[8]:

    式中:w0為初始撓度,mm;wopl為加強(qiáng)筋間板格的最大初始撓度,mm;Boi為初始變形幅值,mm。

    加筋板結(jié)構(gòu)的一階屈曲模態(tài)是最容易產(chǎn)生的失穩(wěn)狀態(tài),對(duì)加筋板的穩(wěn)定性有重要影響。一般認(rèn)為,取一階模態(tài)作為結(jié)構(gòu)的初始缺陷比較合理。本文在考慮初始撓度時(shí),并不是考慮焊接影響產(chǎn)生的撓度變形,而是考慮長期外載荷作用下產(chǎn)生的撓度變形,這種撓度變形在加載之前就已經(jīng)存在,故稱為初始撓度。根據(jù)艦船甲板實(shí)際情況,假設(shè)初始撓度為:

    考慮半波數(shù)對(duì)彈性失穩(wěn)壓力的影響時(shí),限定 y 方向?yàn)橐粋€(gè)半波,只考慮 x 方向半波數(shù)的影響。

    本文擬采用式(3)的初始撓度形式模擬實(shí)際初始撓度情況,來研究整體初始撓度的幅值和半波數(shù)對(duì)加筋板極限承載力的影響。

    根據(jù)生產(chǎn)實(shí)踐,當(dāng)強(qiáng)力甲板厚度在 7~10 mm 之間時(shí),鋼板的絕對(duì)變形量達(dá)到 35 mm 時(shí),才會(huì)對(duì)鋼板進(jìn)行矯正修復(fù);從甲板變形測(cè)量數(shù)據(jù)的統(tǒng)計(jì)結(jié)果來看,甲板變形一般在 8 mm 左右。為了使研究更具有實(shí)用性和有效性,本文將初始撓度幅值設(shè)定在 0~30 mm 之間進(jìn)行研究。國內(nèi)關(guān)于甲板板架變形的半波數(shù)測(cè)量研究內(nèi)容較少,國外有針對(duì)船廠建造中船體變形情況測(cè)量的研究[9],初始撓度半波數(shù)一般在 3 個(gè)以內(nèi)。本文根據(jù)初始撓度實(shí)際情況,決定最多取 6 個(gè)半波作為研究對(duì)象。

    在有限元軟件 Ansys 中,可以通過對(duì)各個(gè)節(jié)點(diǎn)施加具體的法向位移來實(shí)現(xiàn)初始撓度的設(shè)計(jì)。圖 2 所示為具有 1 mm 幅值、一個(gè)半波初始撓度的加筋板。

    2 含初始撓度加筋板極限承載應(yīng)力特性分析

    當(dāng)計(jì)算模型不具有初始撓度時(shí),計(jì)算所得極限載荷 σult為 234.978 MPa。通過不斷改變初始撓度的幅值和半波數(shù),計(jì)算相應(yīng)的極限載荷,計(jì)算數(shù)據(jù)見表 1。

    2.1 初始撓度幅值對(duì)加筋板極限載荷的影響

    根據(jù)計(jì)算數(shù)據(jù)可以繪制出極限載荷與初始撓度幅值的關(guān)系圖,從而可以判定,在某一半波數(shù)情況下,極限載荷隨初始撓度幅值的變化規(guī)律。為了更加方便比較不同半波數(shù)的初始撓度情況下,極限載荷隨幅值的變化規(guī)律,將各個(gè)半波數(shù)下的曲線繪制到同一個(gè)坐標(biāo)系下,并考慮將初始撓度幅值和極限載荷轉(zhuǎn)化為無因次量,取初始撓度幅值 f 與加筋板厚度 t 的比值為自變量,取含初始撓度加筋板極限載荷 σ 與不含初始撓度加筋板極限載荷 σult比值作為因變量,如圖 3 所示。

    從圖 3 中可以看出,雖然初始撓度的半波數(shù)不同,但是加筋板的極限載荷與初始撓度幅值的變化規(guī)律相似。隨著初始撓度幅值的增大,極限載荷逐漸減小。當(dāng)初始撓度半波數(shù)為 1 和 2 時(shí),曲線變化近似于一次函數(shù),極限載荷隨初始撓度幅值的變化率基本保持不變;當(dāng)初始撓度半波數(shù)大于 2 時(shí),曲線變化近似于指數(shù)函數(shù),極限載荷的變化率隨初始撓度幅值的增大而減小。從整體來看,比值 σ/σult隨著初始撓度半波數(shù)的增大而減小。

    表1 含不同初始撓度加筋板的極限載荷(MPa)Tab. 1 P of the stiffened plate with different initial deflections(MPa)

    2.2 初始撓度半波數(shù)對(duì)加筋板極限載荷的影響

    由于計(jì)算數(shù)據(jù)較多,考慮初始撓度半波數(shù)對(duì)極限載荷的影響,只選取了含幅值為 5 mm,10 mm,15 mm,20 mm,25 mm 和 30 mm 初始撓度的數(shù)據(jù)進(jìn)行分析,如圖 4 所示。

    圖4 顯示了當(dāng)初始撓度幅值不變時(shí),增加初始撓度的半波數(shù),加筋板的極限載荷的變化情況。從圖中可看出,不同半波數(shù)的初始撓度情況下,極限載荷的變化規(guī)律相似,載荷都是隨著半波數(shù)的增加而減小。而且隨著初始撓度半波數(shù)的增大,極限載荷的變化率逐漸減小。

    2.3 典型剖面的應(yīng)力分布特點(diǎn)

    根據(jù)對(duì)中修船體甲板變形的勘驗(yàn)測(cè)量數(shù)據(jù)可知,船體加筋板出現(xiàn)的初始撓度幅值大部分在 10 mm 左右,具有一個(gè)半波形式,因此只討論含有幅值 10 mm、具有一個(gè)半波初始撓度的加筋板應(yīng)力分布特點(diǎn)。經(jīng)過有限元計(jì)算,當(dāng)含有一個(gè)半波、幅值 10 mm 初始撓度的加筋板達(dá)到其極限承載力時(shí),其結(jié)構(gòu)變形云圖如圖 5所示,結(jié)構(gòu)應(yīng)力分布云圖如圖 6 所示。

    從圖 5 和圖 6 可看出,加筋板 X 方向中橫剖面處 變形位移最大,且位移沿 Y 方向具有一定遞變規(guī)律。此剖面也是加筋板初始撓度最大幅值位置,應(yīng)力分布具有代表性。因此,取加筋板中橫剖面的節(jié)點(diǎn),觀察其應(yīng)力分布情況,如圖 7 所示。

    圖7 顯示了當(dāng)初始撓度為一個(gè)半波時(shí),不同撓度幅值下的加筋板中橫剖面的應(yīng)力分布情況。從圖 7 中可以看出,與不含初始撓度加筋板的應(yīng)力分布(圖形比例導(dǎo)致分布規(guī)律顯示為一條直線)相比,含初始撓度的加筋板中剖面的應(yīng)力水平成曲線變化,應(yīng)力分布成 U字型。中間部分所受合應(yīng)力較小,邊界附近所受的合應(yīng)力較大,應(yīng)力整體變化規(guī)律是由中間向兩邊逐漸增大。隨著初始撓度加筋板幅值增大,中橫剖面處應(yīng)力水平整體降低,且兩側(cè)邊界處應(yīng)力水平下降速率較大。

    選取中橫剖面中變形位移最大點(diǎn),即含初始撓度加筋板達(dá)到極限承載力時(shí)的位移最大點(diǎn),繪制該點(diǎn)的載荷位移曲線,以反映節(jié)點(diǎn)在加載過程中的位移變化,如圖 8 所示。從圖 8 可看出,隨著初始撓度的增大,極限承載點(diǎn)降低,到達(dá)極限承載點(diǎn)的位移逐漸變大,曲線的變化率逐漸降低,趨于平緩。這說明對(duì)于初始撓度幅值較大的加筋板,在加載過程中更容易產(chǎn)生結(jié)構(gòu)變形,更容易達(dá)到極限承載點(diǎn)。

    3 初始撓度影響因子確定

    由計(jì)算可知,含初始撓度加筋板的極限載荷與不 含初始撓度加筋板的極限載荷存在一定差距。后者與前者差值反映了初始撓度對(duì)加筋板極限載荷的影響,因此定義兩者差值與后者的比值的絕對(duì)值為影響因子λ(0 < λ < 1),即

    參考相關(guān)文獻(xiàn)對(duì)含有初始缺陷的加肋圓柱殼理論失穩(wěn)壓力修正系數(shù)研究的相關(guān)方法[10–11],可以首先初步假設(shè)影響因子與加筋板的長寬比、板厚、加強(qiáng)筋間距等幾何尺寸有關(guān),然后分析含初始撓度加筋板極限載荷與各主要尺寸參數(shù)的關(guān)系,計(jì)算模型是利用一節(jié)的有限元模型作為基礎(chǔ),進(jìn)行相應(yīng)的參數(shù)改變。下面將討論含有幅值 10 mm、一個(gè)半波初始撓度的加筋板的影響因子。

    3.1 加筋板長寬比與影響因子的關(guān)系

    根據(jù)相關(guān)對(duì)船體加筋板尺寸的統(tǒng)計(jì)數(shù)據(jù)[12],船體加筋板板長 a 與加強(qiáng)筋間距 s 的比值應(yīng)該滿足關(guān)系:2 ≤a/s ≤ 6。因此,控制加強(qiáng)筋間距不變,通過改變加筋板的長度 a 來改變加筋板的長寬比 a/b。

    已知加筋板計(jì)算模型加強(qiáng)筋間距 s = 400 mm,具有 6 根加強(qiáng)筋,加筋板寬度 b = 2 800 mm,所以加筋板長度取值范圍為 800 ≤ a ≤ 2 400,均勻取值:800 mm,1 000 mm,1 200 mm,1 400 mm,1 600 mm,1 800 mm,2 000 mm,2 200 mm,2 400 mm。經(jīng)過計(jì)算可以得到影響因子如表 2??梢詫⒈?2 數(shù)據(jù)放在直角坐標(biāo)系中顯示,如圖 9 所示。

    從圖 9 可看出,隨著加筋板長寬比 a/b 增大,影響 因子 λ 逐漸增大,且近似成二次函數(shù)規(guī)律變化。利用Matlab 可以擬合二次函數(shù)表達(dá)式:

    表2 λ 隨加筋板長寬比 a/b 的變化Tab. 2 λ varies with a/b

    擬合函數(shù)數(shù)據(jù)與原始數(shù)據(jù)最大誤差在 0.5% 以內(nèi),符合工程要求,擬合函數(shù)符合實(shí)際情況。由變化規(guī)律可知,當(dāng)加筋板長寬比 a/b 增大時(shí),加筋板初始撓度對(duì)其極限載荷影響增大。在滿足鋼板材料充分利用和結(jié)構(gòu)強(qiáng)度的情況下,應(yīng)盡量縮小加筋板的長寬比,以保持加筋板產(chǎn)生撓度變形時(shí),其對(duì)加筋板極限承載力的影響較小。利用式(5)可以計(jì)算確定其最小值,在最小值附近選取合適的長寬比數(shù)值,從而將初始撓度的影響降到最低。

    3.2 加筋板厚度與影響因子的關(guān)系

    船用加筋板的厚度一般在 8 mm 左右,對(duì)于特殊的大型船舶和軍船,加筋板的厚度可能會(huì)超過 10 mm。因此,加筋板的厚度在 5 ~ 14 mm 之間均勻取值進(jìn)行計(jì)算,計(jì)算結(jié)果見表 3。為了保持變量為無量綱因子,將加筋板厚度 t 與加筋板寬度 b 的比值作為因變量。將計(jì)算結(jié)果繪制成圖,如圖 10 所示。

    從圖 10 可看出,當(dāng)板厚較小時(shí),隨著加筋板 t/b增大,影響因子 λ 逐漸減小,且近似成一次函數(shù)規(guī)律變化;當(dāng)板厚增大到一定厚度,影響因子開始逐漸增大,但是增加幅度非常小,可忽略不計(jì)。所以,在結(jié)構(gòu)設(shè)計(jì)時(shí),增加板厚不僅可以保證結(jié)構(gòu)強(qiáng)度要求,而且可以減小初始撓度對(duì)加筋板極限承載力的影響。

    3.3 加強(qiáng)筋間距與影響因子的關(guān)系

    考慮加強(qiáng)筋間距對(duì)影響因子 λ 的影響,對(duì)加強(qiáng)筋間距的控制,主要是通過改變加強(qiáng)筋的數(shù)量來實(shí)現(xiàn),即保持加筋板的寬度 b 保持不變。因此,加筋板數(shù)量取 4 ~ 9 進(jìn)行計(jì)算,計(jì)算結(jié)果見表 4。為了保持變量為無量綱因子,將加強(qiáng)筋間距 s 與加筋板寬度 b 的比值作為因變量。將計(jì)算結(jié)果繪制成圖,如圖 11 所示。

    從圖 11 可看出,當(dāng)加強(qiáng)筋間距較小時(shí),影響因子隨其變化不明顯,近似為一條平行線;隨著加強(qiáng)筋間距增大,影響因子 λ 逐漸增大,且增大速率也逐漸增大。即當(dāng)加強(qiáng)筋數(shù)量較小時(shí),影響因子數(shù)值較高,加筋板的極限承載力易受初始撓度影響。對(duì)于密加筋板(加強(qiáng)筋數(shù)量 n ≥ 6),影響因子可以保持較小的數(shù)值,且此時(shí)加強(qiáng)筋數(shù)量對(duì)影響因子變化影響較小。所以在結(jié)構(gòu)設(shè)計(jì)時(shí),應(yīng)用密加筋板設(shè)計(jì)合適的加強(qiáng)筋數(shù)量,更有利于減小初始撓度對(duì)加筋板極限承載力的影響。

    表3 λ 隨 t/b 的變化規(guī)律Tab. 3 λ varies with t/b

    表4 λ 隨 s/b 的變化規(guī)律Tab. 4 λ varies with s/b

    4 結(jié) 語

    本文利用非線性數(shù)值計(jì)算方法,對(duì)含有初始撓度的加筋板在軸向壓力作用下的極限承載力進(jìn)行了非線性分析,并得到以下結(jié)論:

    1)初始撓度對(duì)加筋板的極限承載力影響較大,且隨著初始撓度的幅值和半波數(shù)的增大,加筋板的極限承載能力會(huì)逐漸降低,在加筋板長期使用過程中,應(yīng)盡量減小其撓度變形的幅值和半波數(shù);

    2)與不含有初始撓度的加筋板相比,當(dāng)達(dá)到結(jié)構(gòu)的極限承載力時(shí),含初始撓度加筋板中橫剖面內(nèi)的應(yīng)力分布成 U 字型,中間應(yīng)力水平最低,依次向邊界方向逐漸增大;

    3)含初始撓度加筋板影響因子與加筋板的長寬比、板厚和加強(qiáng)筋間距有關(guān);為了減小初始撓度對(duì)加筋板極限承載力的影響,可以在設(shè)計(jì)過程中,適當(dāng)減小加筋板的長寬比和加強(qiáng)筋間距,增大板厚。

    [1]張濤, 劉土光, 趙耀, 等. 初始缺陷加筋板的屈曲與后屈曲分析[J]. 船舶力學(xué), 2003, 7(1): 79–83. ZHANG Tao, LIU Tu-guang, ZHAO Yao, et al. Buckling and postbuckling of imperfect stiffened plates[J]. Journal of Ship Mechaniacs, 2003, 7(1): 79–83.

    [2]譚開忍, 李小平. 船體結(jié)構(gòu)極限強(qiáng)度研究進(jìn)展[J]. 船舶, 2006(5): 19–25. TAN Kai-ren, LI Xiao-ping. Research development of ultimate strength of hull structure[J]. Ship & Boat, 2006(5): 19–25.

    [3]UEDA Y, RASHED S M H, PAIK J K. Effective width of rectangular plates subjected to combined loads[J]. Journal of the Society of Naval Architects of Japan, 1986, 159: 269–281.

    [4]Nonlinear finite element method models for ultimate strength analysis of steel stiffened-plate structures under combined biaxial compression and lateral pressure actions —Part I: Plate elements.

    [5]萬育龍, 朱旭光. 加筋板屈曲和極限強(qiáng)度有限元計(jì)算方法研究[J]. 船海工程, 2013, 42(6): 17–21. WAN Yu-long, ZHU Xu-guang. Studies on the nonlinear finite element for buckling and ultimate strength of stiffened panels[J]. Ship & Ocean Engineering, 2013, 42(6): 17–21.

    [6]羅剛, 吳國民, 湯剛, 等. 側(cè)壓對(duì)加筋板極限強(qiáng)度的影響分析[J]. 武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版), 2014, 38(6): 1400–1403. LUO Gang, WU Guo-min, TANG Gang, et al. Analyse the influence of lateral pressure on the ultimate strength of the stiffen plate[J]. Journal of Wuhan University of Technology(Transportation Science & Engineering), 2014, 38(6): 1400–1403.

    [7]單成巍. 循環(huán)載荷作用下船體結(jié)構(gòu)的極限強(qiáng)度非線性有限元分析[D]. 武漢: 武漢理工大學(xué), 2013. SHAN Cheng-wei. Nonlinear finite element analysis for the ultimate strength of ship structures under cyclic loading[D]. Wuhan: Wuhan University of Technology, 2013.

    [8]楊淼. 基于HCSR的加筋板格非線性屈曲研究[D]. 大連: 大連理工大學(xué), 2011. YANG Miao. Analysis of non-linear buckling strength of stiffened plate based on harmonized common structural rules [D]. Dalian: Dalian University of Technology, 2011.

    [9]TEIXEIRA A P, IVANOV L D, SOARES C G . Assessment of characteristic values of the ultimate strength of corroded steel plates with initial imperfections[J]. Engineering Structures,2013, 56: 517–527.

    [10]DONNELL L H. Effect of imperfection on buckling of thin cylinders under external pressure[J]. Journal of Applied Mechanics. Trans. ASME, 1956, 78.

    [11]張二. 初始幾何缺陷對(duì)錐-環(huán)-柱結(jié)合殼力學(xué)性能影響研究[D]. 武漢: 海軍工程大學(xué), 2015. ZHANG Er. Research on the impact of initial imperfection to cone-toroid-cylinder combined shell’s mechanical property[D]. Wuhan: Naval University of Engineering, 2015.

    [12]張曉丹, 楊平. 加筋板在軸向壓力下的極限強(qiáng)度研究[J]. 武漢理工大學(xué)學(xué)報(bào), 2011, 35(2): 305–308. ZHANG Xiao-dan, YANG Ping. Ultimate strength of stiffened plate under axial compression[J]. Journal of Wuhan University of Technology, 2011, 35(2): 305–308.

    The ultimate bearing capacity analysis of the stiffened plate with initial deflection

    LIU Chun-zheng, WU Fan, MU Jin-lei
    (Naval University of Engineering, Department of Naval Architecture, Wuhan 430033, China)

    The initial deflection existed on the deck of the surface warship at long service and it brought bad effects on the bearing capacity. Study on the effect of initial deflection on the stiffened plate was meaningful since it was the main unit of the deck. For the research, the mode of the initial deflection was regarded as the double trigonometric series based on the practical truth and the effect of the amplitude and the number of half-waves on the critical stress entirely was analyzed by computing with Ansys. The stress state of the typical location was also discussed. The impact factor of the initial deflection was obtained. Results indicate that with the increasing of amplitude and the number of half-waves, the critical stress decreases. For the stiffened plate with one kind of initial deflection, the impact factor is up to the length-width ratio, the thickness of the stiffened plate and the space between stiffeners.

    initial deflection;stiffened plate;ultimate bearing capacity;stress state;impact factor

    U663.6

    A

    1672–7619(2017)05–0010–06

    10.3404/j.issn.1672–7619.2017.05.003

    2016–07–11;

    2016–08–22

    國家自然科學(xué)基金資助項(xiàng)目(51309231)

    劉春正(1991–),男,碩士研究生,研究方向?yàn)榇敖Y(jié)構(gòu)強(qiáng)度與振動(dòng)。

    猜你喜歡
    半波筋板加強(qiáng)筋
    薄壁扇形筋板擠壓成形開裂抑制及翻轉(zhuǎn)展寬策略
    Y 分支光學(xué)調(diào)制器的半波電壓特性研究
    汽車零件加強(qiáng)筋結(jié)構(gòu)設(shè)計(jì)
    汽車工程師(2019年7期)2019-08-12 07:25:04
    越野車側(cè)壁加強(qiáng)筋板工藝改進(jìn)
    模具制造(2019年3期)2019-06-06 02:10:52
    磁各向異性媒質(zhì)中半波天線的輻射功率
    半波帶法研究夫瑯禾費(fèi)單縫衍射的局限性分析*
    用半波帶法推導(dǎo)單縫夫瑯禾費(fèi)衍射暗紋條件的相關(guān)討論
    加強(qiáng)筋的截面尺寸對(duì)注塑制品翹曲變形的影響
    塑料制造(2016年5期)2016-06-15 20:27:39
    低速?zèng)_擊下復(fù)合材料加筋板的損傷阻抗性能
    基于板梁組合理論的正交加筋板聲振特性分析
    视频中文字幕在线观看| 国产亚洲5aaaaa淫片| 青春草国产在线视频| 日韩一区二区三区影片| 人妻少妇偷人精品九色| 亚洲精品456在线播放app| 精品人妻一区二区三区麻豆| 91精品一卡2卡3卡4卡| av福利片在线观看| 久久久亚洲精品成人影院| 亚洲精品亚洲一区二区| 久久久久精品久久久久真实原创| 亚洲国产欧美在线一区| 久久久久久久久大av| 99久久精品热视频| 欧美日韩在线观看h| 最黄视频免费看| 国产欧美另类精品又又久久亚洲欧美| 激情 狠狠 欧美| 国产成人a∨麻豆精品| 春色校园在线视频观看| 国产 精品1| 国产成人精品一,二区| 天堂俺去俺来也www色官网| 女的被弄到高潮叫床怎么办| 免费黄网站久久成人精品| 国产在线一区二区三区精| 亚洲第一av免费看| 亚洲av中文字字幕乱码综合| 欧美变态另类bdsm刘玥| 成人亚洲精品一区在线观看 | 亚洲av日韩在线播放| 欧美精品一区二区大全| av.在线天堂| 精品国产露脸久久av麻豆| 欧美日韩综合久久久久久| 在线观看av片永久免费下载| 少妇人妻 视频| a级毛片免费高清观看在线播放| av卡一久久| 久久精品夜色国产| 国产爽快片一区二区三区| 久久这里有精品视频免费| 黄色一级大片看看| 全区人妻精品视频| 在线观看免费日韩欧美大片 | 另类亚洲欧美激情| 日韩一区二区视频免费看| 国产伦在线观看视频一区| 你懂的网址亚洲精品在线观看| 日韩欧美 国产精品| 五月伊人婷婷丁香| 少妇的逼好多水| 久久久精品94久久精品| 国产视频首页在线观看| 中国国产av一级| 深爱激情五月婷婷| 黄色怎么调成土黄色| 嫩草影院新地址| 2018国产大陆天天弄谢| 亚洲,一卡二卡三卡| 热99国产精品久久久久久7| 国产精品蜜桃在线观看| 国产色爽女视频免费观看| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久久人人人人人人| 日本色播在线视频| 国产精品一区二区在线不卡| 女性被躁到高潮视频| 精品人妻一区二区三区麻豆| 色5月婷婷丁香| 欧美最新免费一区二区三区| 免费观看在线日韩| 欧美亚洲 丝袜 人妻 在线| 色婷婷av一区二区三区视频| 欧美少妇被猛烈插入视频| 日韩免费高清中文字幕av| 国产无遮挡羞羞视频在线观看| 最近最新中文字幕大全电影3| 亚洲久久久国产精品| 国产成人精品久久久久久| 一区二区av电影网| 少妇裸体淫交视频免费看高清| 18禁裸乳无遮挡免费网站照片| 国产无遮挡羞羞视频在线观看| 夜夜爽夜夜爽视频| 亚洲色图av天堂| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| 国产亚洲91精品色在线| 六月丁香七月| 最近手机中文字幕大全| 五月天丁香电影| 晚上一个人看的免费电影| 小蜜桃在线观看免费完整版高清| 欧美激情极品国产一区二区三区 | 一本—道久久a久久精品蜜桃钙片| 国产视频首页在线观看| 男人和女人高潮做爰伦理| 青春草国产在线视频| 成人一区二区视频在线观看| 黑丝袜美女国产一区| 国产黄频视频在线观看| 在线观看美女被高潮喷水网站| 亚洲欧美日韩卡通动漫| 乱系列少妇在线播放| 国产爽快片一区二区三区| 国产在线男女| 中文资源天堂在线| 国产一区二区三区综合在线观看 | 一级av片app| 一级av片app| 日日啪夜夜爽| 成人综合一区亚洲| 国产伦理片在线播放av一区| 波野结衣二区三区在线| 亚洲精品视频女| 亚洲欧美成人综合另类久久久| 啦啦啦在线观看免费高清www| 丰满乱子伦码专区| 黑人猛操日本美女一级片| 亚洲精品国产成人久久av| 青青草视频在线视频观看| 91久久精品国产一区二区成人| 麻豆乱淫一区二区| 亚洲婷婷狠狠爱综合网| 午夜福利影视在线免费观看| 少妇熟女欧美另类| 欧美 日韩 精品 国产| av黄色大香蕉| 在线观看一区二区三区激情| 97超碰精品成人国产| 高清视频免费观看一区二区| 中文天堂在线官网| 日产精品乱码卡一卡2卡三| 欧美bdsm另类| 美女xxoo啪啪120秒动态图| 国产亚洲一区二区精品| 国产日韩欧美亚洲二区| a级一级毛片免费在线观看| 三级经典国产精品| 精品人妻偷拍中文字幕| xxx大片免费视频| 一二三四中文在线观看免费高清| 成年免费大片在线观看| 丝袜喷水一区| 在线精品无人区一区二区三 | 黄色怎么调成土黄色| 亚洲国产精品成人久久小说| 国产乱来视频区| tube8黄色片| 在线精品无人区一区二区三 | 18禁裸乳无遮挡免费网站照片| 99热6这里只有精品| 日日啪夜夜爽| 另类亚洲欧美激情| av.在线天堂| 亚洲欧美精品专区久久| 国产欧美另类精品又又久久亚洲欧美| 天堂8中文在线网| 国产亚洲精品久久久com| 国产精品99久久99久久久不卡 | 日本免费在线观看一区| 亚洲精品,欧美精品| 国产亚洲91精品色在线| 久久韩国三级中文字幕| 99久久精品国产国产毛片| 亚洲国产成人一精品久久久| 日韩一区二区视频免费看| 中文资源天堂在线| 亚洲精品乱码久久久久久按摩| 亚州av有码| 91精品国产国语对白视频| 亚洲成人一二三区av| 久久久久久久亚洲中文字幕| 亚洲国产精品999| 全区人妻精品视频| 精华霜和精华液先用哪个| 在线观看美女被高潮喷水网站| 大香蕉97超碰在线| 国产乱来视频区| 在线观看一区二区三区激情| 久久久久久久久久成人| 国产熟女欧美一区二区| 欧美精品人与动牲交sv欧美| 精华霜和精华液先用哪个| 国国产精品蜜臀av免费| 18禁在线无遮挡免费观看视频| 男人狂女人下面高潮的视频| 99久久精品一区二区三区| 久久人妻熟女aⅴ| 亚洲美女视频黄频| 欧美性感艳星| 简卡轻食公司| 99热这里只有精品一区| 国产成人精品一,二区| 伦理电影大哥的女人| 中文字幕av成人在线电影| 日本vs欧美在线观看视频 | 在线观看国产h片| 毛片女人毛片| 伊人久久国产一区二区| 人人妻人人爽人人添夜夜欢视频 | 国产伦精品一区二区三区四那| 久久久久久九九精品二区国产| 黑人高潮一二区| 欧美日韩国产mv在线观看视频 | 亚洲激情五月婷婷啪啪| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片久久久久久久久女| 精品人妻熟女av久视频| 内地一区二区视频在线| 人妻一区二区av| 制服丝袜香蕉在线| 国产精品av视频在线免费观看| 欧美性感艳星| 一本一本综合久久| videos熟女内射| 亚洲欧美日韩东京热| 欧美亚洲 丝袜 人妻 在线| 国产男女内射视频| 精品久久久精品久久久| 婷婷色av中文字幕| 国产 一区 欧美 日韩| 欧美激情国产日韩精品一区| 中文字幕av成人在线电影| 一本色道久久久久久精品综合| 毛片一级片免费看久久久久| 99热这里只有是精品在线观看| 成年女人在线观看亚洲视频| 高清午夜精品一区二区三区| 国产精品国产三级专区第一集| 精品久久久久久久久av| 亚洲精品久久午夜乱码| 亚洲国产精品成人久久小说| 嫩草影院新地址| 亚洲精品中文字幕在线视频 | 亚洲欧美中文字幕日韩二区| 观看免费一级毛片| 亚洲av二区三区四区| 欧美日韩视频精品一区| av一本久久久久| 毛片女人毛片| 日韩视频在线欧美| 亚洲欧美一区二区三区黑人 | 热99国产精品久久久久久7| 精品一区二区三区视频在线| 日韩视频在线欧美| 黑丝袜美女国产一区| 久久精品国产自在天天线| 午夜老司机福利剧场| 精品亚洲乱码少妇综合久久| 亚洲国产日韩一区二区| 少妇丰满av| 97热精品久久久久久| 香蕉精品网在线| 麻豆成人av视频| 亚洲熟女精品中文字幕| 偷拍熟女少妇极品色| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 欧美日韩视频精品一区| 综合色丁香网| 午夜精品国产一区二区电影| 伦理电影大哥的女人| 亚洲精品,欧美精品| 舔av片在线| 晚上一个人看的免费电影| 内射极品少妇av片p| 久久精品人妻少妇| 我要看日韩黄色一级片| 欧美3d第一页| 欧美激情极品国产一区二区三区 | 日韩伦理黄色片| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 成人亚洲精品一区在线观看 | 国产成人免费无遮挡视频| 亚洲精品国产av蜜桃| 精品酒店卫生间| 搡女人真爽免费视频火全软件| av一本久久久久| 精品人妻熟女av久视频| 精品一区二区三卡| 女性生殖器流出的白浆| 免费播放大片免费观看视频在线观看| 少妇裸体淫交视频免费看高清| 国产无遮挡羞羞视频在线观看| 人妻制服诱惑在线中文字幕| 免费播放大片免费观看视频在线观看| 亚洲国产最新在线播放| 精品人妻一区二区三区麻豆| 在线精品无人区一区二区三 | 久久久久国产精品人妻一区二区| 尾随美女入室| 午夜日本视频在线| 亚洲欧洲日产国产| 国产精品99久久99久久久不卡 | 边亲边吃奶的免费视频| 久久人人爽人人爽人人片va| 日韩亚洲欧美综合| 精品亚洲成a人片在线观看 | 天美传媒精品一区二区| 极品教师在线视频| 三级经典国产精品| 日韩精品有码人妻一区| 精品久久久精品久久久| 涩涩av久久男人的天堂| 亚洲电影在线观看av| 亚洲国产欧美人成| 97精品久久久久久久久久精品| 18+在线观看网站| 这个男人来自地球电影免费观看 | 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕| 男女免费视频国产| 我的女老师完整版在线观看| 久久国内精品自在自线图片| 精品视频人人做人人爽| 身体一侧抽搐| 亚洲内射少妇av| 熟女人妻精品中文字幕| av视频免费观看在线观看| 日韩一本色道免费dvd| 日韩强制内射视频| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 亚洲怡红院男人天堂| 毛片女人毛片| 免费人成在线观看视频色| 久久久午夜欧美精品| 精品少妇久久久久久888优播| 久久人人爽人人片av| 人妻 亚洲 视频| 亚洲美女黄色视频免费看| 少妇裸体淫交视频免费看高清| 18禁动态无遮挡网站| av天堂中文字幕网| 欧美精品国产亚洲| 在线免费观看不下载黄p国产| 欧美+日韩+精品| 精品久久久久久久久亚洲| 婷婷色综合www| 久久久欧美国产精品| 少妇人妻一区二区三区视频| 亚洲成人手机| 免费看光身美女| 我要看日韩黄色一级片| 久久久欧美国产精品| 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 亚洲欧洲国产日韩| 欧美高清性xxxxhd video| 精品国产三级普通话版| 婷婷色av中文字幕| 久久99精品国语久久久| 免费在线观看成人毛片| 国产永久视频网站| 国产美女午夜福利| 亚洲av中文av极速乱| 在现免费观看毛片| 国产av码专区亚洲av| 亚洲aⅴ乱码一区二区在线播放| 91狼人影院| 欧美精品国产亚洲| 精品久久久噜噜| 99久久精品热视频| 精品亚洲成a人片在线观看 | 国产在线男女| 黑人猛操日本美女一级片| 联通29元200g的流量卡| 激情五月婷婷亚洲| 国产精品不卡视频一区二区| 在线观看国产h片| 国产av一区二区精品久久 | av黄色大香蕉| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产av蜜桃| 亚洲精品日韩av片在线观看| 亚洲av不卡在线观看| av播播在线观看一区| 秋霞在线观看毛片| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品在线观看| 中文字幕av成人在线电影| 国产片特级美女逼逼视频| 欧美区成人在线视频| 卡戴珊不雅视频在线播放| 国产在线视频一区二区| 日韩成人av中文字幕在线观看| 亚洲久久久国产精品| 久久青草综合色| 尾随美女入室| 超碰av人人做人人爽久久| 尾随美女入室| 久久ye,这里只有精品| 久久毛片免费看一区二区三区| 国产高清国产精品国产三级 | 日日摸夜夜添夜夜爱| 亚洲国产日韩一区二区| 欧美成人a在线观看| 最近中文字幕2019免费版| 18禁裸乳无遮挡免费网站照片| 精品国产露脸久久av麻豆| av国产久精品久网站免费入址| 看十八女毛片水多多多| 好男人视频免费观看在线| 国产黄片美女视频| 亚洲av免费高清在线观看| 永久网站在线| 日韩亚洲欧美综合| 久久久久久九九精品二区国产| 国产精品99久久99久久久不卡 | 熟女电影av网| 国产精品久久久久成人av| 亚洲国产高清在线一区二区三| 亚洲va在线va天堂va国产| 日本免费在线观看一区| 国产美女午夜福利| 亚洲欧美精品专区久久| 久久久久久久国产电影| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人精品一区久久| 精品久久久久久久久av| 赤兔流量卡办理| 一级毛片久久久久久久久女| 精品酒店卫生间| av国产久精品久网站免费入址| 精品国产露脸久久av麻豆| 日本欧美国产在线视频| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲网站| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 亚洲国产欧美在线一区| av视频免费观看在线观看| 免费大片18禁| 亚洲经典国产精华液单| 亚洲四区av| 午夜福利高清视频| 美女高潮的动态| 日韩中字成人| 精品亚洲成a人片在线观看 | 欧美激情国产日韩精品一区| 91精品国产国语对白视频| av天堂中文字幕网| 亚洲精品日本国产第一区| 这个男人来自地球电影免费观看 | 亚洲精品成人av观看孕妇| 欧美亚洲 丝袜 人妻 在线| 成人二区视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲第一av免费看| 少妇人妻一区二区三区视频| 亚洲国产成人一精品久久久| 国产精品一及| 日韩欧美一区视频在线观看 | 亚洲精品456在线播放app| 在线观看一区二区三区| 国产色爽女视频免费观看| 久久久精品94久久精品| 色5月婷婷丁香| 91精品伊人久久大香线蕉| 国产熟女欧美一区二区| 国产成人aa在线观看| 亚洲欧美日韩无卡精品| 亚洲国产精品国产精品| 一级片'在线观看视频| 另类亚洲欧美激情| 91在线精品国自产拍蜜月| 美女脱内裤让男人舔精品视频| 中文乱码字字幕精品一区二区三区| 日日啪夜夜撸| 如何舔出高潮| 一本—道久久a久久精品蜜桃钙片| 久久久精品免费免费高清| 免费大片18禁| 99视频精品全部免费 在线| 麻豆成人午夜福利视频| av在线观看视频网站免费| 在线亚洲精品国产二区图片欧美 | 看非洲黑人一级黄片| 国产在线一区二区三区精| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美精品专区久久| 亚洲国产精品国产精品| 精品久久久久久久久av| 免费大片黄手机在线观看| 插阴视频在线观看视频| 男女边摸边吃奶| 国产v大片淫在线免费观看| www.色视频.com| 一级毛片aaaaaa免费看小| 噜噜噜噜噜久久久久久91| 国产欧美亚洲国产| 中文乱码字字幕精品一区二区三区| 亚洲精品久久久久久婷婷小说| 熟妇人妻不卡中文字幕| 又粗又硬又长又爽又黄的视频| 人妻一区二区av| 美女内射精品一级片tv| 男的添女的下面高潮视频| 欧美 日韩 精品 国产| av免费在线看不卡| 免费看不卡的av| 日日摸夜夜添夜夜爱| 我的老师免费观看完整版| 成年免费大片在线观看| 日韩制服骚丝袜av| 午夜免费男女啪啪视频观看| 高清不卡的av网站| 1000部很黄的大片| 婷婷色av中文字幕| 五月玫瑰六月丁香| 国产精品99久久久久久久久| 一本久久精品| 国产精品成人在线| 老师上课跳d突然被开到最大视频| 亚洲精品国产色婷婷电影| 精品酒店卫生间| 久久久久精品性色| 日韩中字成人| 亚洲三级黄色毛片| 久久99热这里只频精品6学生| 国产熟女欧美一区二区| 成年免费大片在线观看| 七月丁香在线播放| 人妻一区二区av| 最后的刺客免费高清国语| 亚洲精品久久久久久婷婷小说| 欧美成人午夜免费资源| 色哟哟·www| 九草在线视频观看| 男人爽女人下面视频在线观看| 精华霜和精华液先用哪个| 99热这里只有是精品50| 精品少妇久久久久久888优播| 夜夜看夜夜爽夜夜摸| 美女高潮的动态| 色吧在线观看| 汤姆久久久久久久影院中文字幕| 精品亚洲乱码少妇综合久久| 我要看日韩黄色一级片| 高清午夜精品一区二区三区| 欧美丝袜亚洲另类| 七月丁香在线播放| 在线播放无遮挡| 午夜福利影视在线免费观看| 久久99精品国语久久久| av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 国产有黄有色有爽视频| 久久久久性生活片| 亚洲天堂av无毛| 九九爱精品视频在线观看| 国产真实伦视频高清在线观看| 日韩伦理黄色片| 中文字幕制服av| 久久久欧美国产精品| 色婷婷av一区二区三区视频| 高清毛片免费看| 日韩中文字幕视频在线看片 | 亚洲成人av在线免费| 亚洲av免费高清在线观看| 亚洲av欧美aⅴ国产| 超碰97精品在线观看| 国产成人91sexporn| 日韩欧美一区视频在线观看 | 99久久精品国产国产毛片| 波野结衣二区三区在线| 亚洲av免费高清在线观看| 国产精品成人在线| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频 | 日韩,欧美,国产一区二区三区| 国产av国产精品国产| 国产熟女欧美一区二区| 欧美xxxx性猛交bbbb| 亚洲美女黄色视频免费看| 纯流量卡能插随身wifi吗| 国产男人的电影天堂91| 国产精品99久久久久久久久| 天美传媒精品一区二区| 老司机影院毛片| 国产精品99久久99久久久不卡 | 欧美成人午夜免费资源| 高清av免费在线| 色5月婷婷丁香| 久久韩国三级中文字幕| 午夜福利在线观看免费完整高清在| 欧美xxxx性猛交bbbb| 久久精品久久久久久噜噜老黄| 国产日韩欧美在线精品| 国产精品女同一区二区软件| 九九在线视频观看精品| 麻豆乱淫一区二区| 国产成人91sexporn| 国产淫语在线视频| 久久ye,这里只有精品| 女性生殖器流出的白浆| 人妻少妇偷人精品九色| 久久99蜜桃精品久久| 欧美极品一区二区三区四区| 国产欧美另类精品又又久久亚洲欧美| 人妻 亚洲 视频| 免费观看无遮挡的男女| 在线看a的网站| 免费看av在线观看网站| 国产精品一及| 小蜜桃在线观看免费完整版高清| 舔av片在线| 精品酒店卫生间| 99久久精品国产国产毛片| 日日啪夜夜撸| 99热这里只有精品一区|