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

    加筋板輪印載荷分布特性的試驗與數(shù)值分析

    2017-10-12 01:02:28程遠勝
    海洋工程 2017年4期
    關(guān)鍵詞:有限元變形

    劉 聰,程遠勝,張 攀,劉 均

    (華中科技大學 船舶與海洋工程學院,湖北 武漢 430074)

    加筋板輪印載荷分布特性的試驗與數(shù)值分析

    劉 聰,程遠勝,張 攀,劉 均

    (華中科技大學 船舶與海洋工程學院,湖北 武漢 430074)

    輪印載荷的分布對于直升機甲板加筋板結(jié)構(gòu)的設(shè)計至關(guān)重要。通過在橡膠塊上加壓的方式模擬輪印載荷,直接測量橡膠塊與加筋板結(jié)構(gòu)之間的接觸壓力,揭示輪印載荷分布特性,得到輪印載荷在加強筋和鋪板之間的分布具有較強的不均勻性。利用ABAQUS軟件進行輪印載荷分布的數(shù)值計算,結(jié)果表明,數(shù)值計算結(jié)果與試驗結(jié)果吻合較好,試驗工況下分布在加強筋上均布輪印載荷為接觸面上均布輪印載荷的1.279倍。研究了加筋板面板初始缺陷、加強筋高度以及加筋板面板厚度對輪印載荷分布的影響。數(shù)值結(jié)果表明:考慮面板初始缺陷、降低面板厚度和增加加強筋的高度使輪印載荷分配到加強筋上的比例增大,但在相對剛度變化相同的情況下,降低面板厚度和增加加強筋高度對輪印載荷分布的影響不同。

    輪印載荷分布;加筋板結(jié)構(gòu);接觸壓力;試驗;數(shù)值分析

    Abstract: The distribution characteristic of wheel load is critical to the design of stiffened helicopter deck structures. This paper presents an experiment to investigate the distribution characteristic of wheel load represented by contact pressure between the rubber wheel and stiffened plates. Experimental data reveal that the wheel load distribution is highly non-uniform. Simultaneously, the numerical software ABAQUS is employed to analyze the distribution of the wheel load. Results show that the numerical simulation agrees well with the experiment. The wheel load acting on stiffener is 1.279 times of the averaged load on the whole contacting region for the experimental case. Then, the effects of initial deflection of face plate, the height of stiffener and the thickness of face plate are studied using the validated numerical model. The wheel load acting on stiffener would increase with the decrease of the thickness of face plate and the increase of the height of stiffener. In addition, the initial deflection of face plate also results in the increase of wheel load on stiffener. Changing the same degree of the relative stiffness between plate and stiffener by decreasing the thickness of face plate or increasing the height of stiffener has different effect on the distribution characteristic of the wheel load.

    Keywords: wheel load distribution; stiffened plate; contact pressure; experiment; numerical simulation

    直升飛機作為現(xiàn)代海洋工程裝備和大型艦船上的設(shè)備,其在著艦和停機狀態(tài)下通過主輪組向甲板結(jié)構(gòu)傳遞載荷——輪印載荷。輪印載荷具有局部重載和位置不確定等特點,因此對結(jié)構(gòu)設(shè)計提出了較高的要求。

    目前,國內(nèi)外對于承受輪印載荷下的甲板結(jié)構(gòu)進行了較多的研究[1-8]。彭興寧等[9]進行了輪印載荷作用下的甲板結(jié)構(gòu)設(shè)計研究并給出了輪載作用下甲板結(jié)構(gòu)的設(shè)計圖譜;Romanoff等[10]采用理論分析和有限元相結(jié)合的方法研究了激光焊接I型(web-core)夾層板在輪印載荷作用下的變形和應(yīng)力;Jackson等[11]進行了輪印載荷作用下甲板結(jié)構(gòu)響應(yīng)的試驗。上述研究大都將載荷歸結(jié)為接觸面積內(nèi)的均布壓力,未考慮由于輪胎結(jié)構(gòu)非線性導致的輪印載荷在加強筋與鋪板之間的非均勻分布。

    輪印載荷下加筋板結(jié)構(gòu)設(shè)計一般是依據(jù)規(guī)范進行[12]。首先確定矩形輪印接觸面的大小,然后運用相應(yīng)公式進行面板和骨材的設(shè)計。但是,規(guī)范沒有給出當輪印塊跨兩根及以上骨材間距時的設(shè)計方法,且規(guī)范方法對于夾層結(jié)構(gòu)的設(shè)計并不完全適用。因此,需要深入探究輪印載荷的分布特性,為結(jié)構(gòu)設(shè)計提供精確的載荷依據(jù)。

    本文設(shè)計了一個加筋板結(jié)構(gòu),采用橡膠塊與加壓鋼塊組合的方式來模擬輪印載荷[13-14],運用試驗和數(shù)值計算相結(jié)合的方法,分析輪印載荷塊與結(jié)構(gòu)之間的接觸壓力分布,獲取輪印載荷在加筋板上的分布規(guī)律,并分析加筋板面板初始缺陷以及結(jié)構(gòu)參數(shù)對輪印載荷分布規(guī)律的影響。

    1 輪印載荷分布試驗

    1.1試驗?zāi)P?/p>

    圖1 加筋板試驗?zāi)P褪疽釬ig. 1 Sketch showing the stiffened plate

    以某直升機甲板為母型采用等效縮比的方式設(shè)計了一個加筋板結(jié)構(gòu)模型,模型的尺寸參數(shù):總長L=750 mm,總寬B=975 mm,板厚t=6 mm,寬度方向布置4根扁鋼形式的加強筋,間距δ=195 mm,扁鋼尺寸為35 mm×6 mm。模型四周設(shè)置由16b槽鋼(依據(jù)橫梁尺寸按慣性矩縮比)組成的邊框,槽鋼內(nèi)等間距設(shè)置肘板,試驗?zāi)P褪疽馊鐖D1所示。

    為了下文表述方便,現(xiàn)將坐標系定義如下:X軸為垂直于加強筋方向,橫向;Y軸為平行于加強筋方向,縱向;Z軸為加筋板的高度方向,垂向。且坐標原點位于加筋板面板中心處。

    將試驗?zāi)P头胖迷谝黄脚_上,四周邊框采用約束梁進行固定,運用液壓千斤頂施加壓力,試驗實物圖如圖2所示。

    圖2 試驗實物照片F(xiàn)ig. 2 Picture of tested model

    1.2測試系統(tǒng)

    試驗測試的物理量包括橡膠塊與結(jié)構(gòu)之間的接觸壓力、結(jié)構(gòu)特征點的應(yīng)力以及千斤頂施加的力。其中接觸壓力通過MFF系列薄膜壓力測試系統(tǒng)測量,結(jié)構(gòu)上的應(yīng)力通過DH-3816靜態(tài)應(yīng)變儀測量,施加的力通過力傳感器測量。

    MFF系列薄膜壓力測試系統(tǒng)由薄膜壓力傳感器、電荷放大器、數(shù)據(jù)采集器、數(shù)據(jù)處理軟件以及標定裝置組成,如圖3所示,其中壓力傳感器為美國Tekscan公司的FlexiForce A201型薄膜壓力傳感器,如圖3(c)所示,其參數(shù)如表1所示。使用過程中首先通過圖3(a)中的標定臺進行傳感器標定,然后將傳感器的感應(yīng)區(qū)定位在測量點上,通過導線將傳感器與電荷放大器相連接,然后再連接數(shù)據(jù)采集器,數(shù)據(jù)采集器通過USB線與計算機相連,通過專用的軟件采集并記錄測量到的壓力值。

    表1 A201型薄膜壓力傳感器參數(shù)Tab. 1 Parameters of A201 pressure sensor

    圖3 薄膜壓力測試系統(tǒng)Fig. 3 Test system of MFF

    1.3傳感器布置

    試驗所使用的傳感器包括A201型薄膜壓力傳感器、電阻式應(yīng)變片以及力傳感器。

    A201型薄膜壓力傳感器布置位置如圖4所示,傳感器的位置沿橫向分兩排布置,將加筋板相鄰骨材間距8等分,1~4號測點布置在縱向?qū)ΨQ軸上,6~8號測點布置在距離縱向?qū)ΨQ軸50 mm的平行線上,5號測點布置在2號和6號測點中間(距縱向?qū)ΨQ軸25 mm),沿橫向相鄰測點之間距離為24.4 mm。應(yīng)變片為中航工業(yè)電測儀器股份有限公司制造的型號為BE350-3BA(11)的電阻應(yīng)變片,其電阻為350.0±0.3 Ω。應(yīng)變片布置在加筋板背面用來測量各關(guān)注測點的局部應(yīng)力,其布置位置如圖5所示,圖中1~10號為雙向片,布置在板格中心和焊縫處;11~13號為單向片,布置在加強筋下端。力傳感器采用東華BLR-1型,量程為2 t。

    圖4 薄膜壓力傳感器布置示意Fig. 4 Sketch showing the arrangement of pressure sensors

    圖5 應(yīng)變片布置示意Fig. 5 Sketch showing the arrangement of strain gages

    1.4試驗工況以及結(jié)果

    根據(jù)輪胎接觸面積,以長、寬、高分別為A'=195 mm、B'=140 mm、H=27 mm的立方體橡膠塊模擬橡膠輪胎。試驗工況為:輪印塊X方向中心線與結(jié)構(gòu)縱對稱軸重合(y=0直線),輪印塊Y方向中心線與結(jié)構(gòu)加強筋方向重合(x=97.5 mm直線)。輪印塊長邊沿X方向,運用千斤頂施加力,力的大小為9 197 N。在該力作用下結(jié)構(gòu)上測量點的應(yīng)力值如表2所示,薄膜壓力傳感器各測點壓力如表3所示。

    表2 應(yīng)力的試驗結(jié)果與有限元結(jié)果對比Tab. 2 Comparison of stress results from simulation and experiment

    表3 接觸壓力的試驗結(jié)果與有限元結(jié)果對比Tab. 3 Comparison of contact pressure results from simulation and experiment

    2 有限元計算

    圖6 有限元模型Fig. 6 FEM model

    采用有限元軟件ABAQUS進行數(shù)值模擬。有限元模型如圖6所示。有限元仿真加筋板結(jié)構(gòu)參數(shù)與試驗?zāi)P蛶缀螀?shù)保持一致(不含邊框)。加筋板結(jié)構(gòu)采用S4R殼單元模擬,網(wǎng)格大小為5 mm,橡膠塊采用C3D8H實體單元模擬,網(wǎng)格大小為3 mm,上述網(wǎng)格劃分經(jīng)收斂性檢查為合適的網(wǎng)格劃分方案。加壓鋼塊約束為剛體,在中心參考點上加約束(UX、UY、ROTX、ROTY、ROTZ)和載荷(FZ=-9 197 N)。在鋼塊與橡膠塊上表面和橡膠塊下表面與加筋板面板之間建立接觸關(guān)系,法向定義為“hard”接觸,切向定義為庫倫摩擦模型[15]。加筋板材料為普通鋼,彈性模量為2.1×1011Pa,密度為7 850 kg/m3,泊松比為0.3,屈服極限為235 MPa。橡膠材料采用Mooney-Rivlin模型[16],彈性模量為2.798 MPa,泊松比為0.499。加筋板邊界條件如圖6所示,結(jié)構(gòu)四邊固支。

    2.1試驗工況計算結(jié)果

    根據(jù)試驗工況,計算所加載荷為9 197 N,表3給出了各測點應(yīng)力的有限元計算值。根據(jù)試驗值、有限元值以及兩者的誤差可以看出大部分測點誤差在10%以內(nèi),表明有限元結(jié)果與試驗結(jié)果吻合較好。F02號測點相對誤差較大主要是因為測點位于焊縫邊緣,應(yīng)力場比較復雜,運用殼單元無法模擬真實情況導致;F06號測點相對誤差大,一方面由于焊縫處應(yīng)力場的復雜性導致,另一方面由于載荷中心位于該測點處,加強筋的存在導致應(yīng)力變化梯度較大,應(yīng)力值受貼片位置的影響較大。

    提取接觸區(qū)域內(nèi)薄膜壓力傳感器測點所在線上的接觸壓力,即y=0 mm直線與y=50 mm直線上鋪板的接觸壓力,并與試驗結(jié)果進行對比,如圖7所示。表4給出了各測點接觸壓力的試驗結(jié)果和有限元結(jié)果,從表中可以看出大部分測點誤差在10%以內(nèi),表明有限元結(jié)果和試驗結(jié)果吻合的較好,4號測點和8號測點相對誤差較大是因為焊接加強筋導致加筋板面板在加強筋處產(chǎn)生幾何變形引起的。

    圖8給出了接觸區(qū)域三維接觸壓力云圖,接觸壓力在橡膠塊中點達到最大值,沿著X軸方向和Y軸方向逐漸衰減。繪出直線y=0 mm與直線x=97.5 mm上鋪板的接觸壓力曲線如圖9所示。接觸壓力的這種分布形態(tài)主要是因為,一方面,結(jié)構(gòu)在直線x=97.5 mm處設(shè)置有一根加強筋,加強筋的存在使得結(jié)構(gòu)在該處的變形較不含加強筋時小,從而該處的接觸壓力大;另一方面,橡膠塊與結(jié)構(gòu)之間切向接觸關(guān)系為庫倫摩擦模型,由于摩擦的存在而導致橡膠塊中部的接觸壓力大[17]。

    圖7 試驗結(jié)果與有限元結(jié)果對比Fig. 7 Comparison of the simulation results and the experiment results

    圖8 三維接觸壓力云圖Fig. 8 3D plot of contact pressure

    圖9 接觸壓力曲線Fig. 9 Contact pressure curves

    根據(jù)甲板縱骨、橫梁跨中最大彎矩計算公式[12],基于彎矩等效原則,計算垂直于加強筋方向的有效寬度Be,然后將有效寬度內(nèi)接觸壓力進行均勻化,得到分布在加強筋上的均布輪印載荷幅值ps,如圖10所示,進而可以求得ps與整個接觸面上的均布輪印載荷幅值pa的比值rp。上述各物理量按如下公式進行計算:

    式中:M0為跨中彎矩;a為構(gòu)件跨距;B'為輪印寬度(橡膠塊的寬度);KM0為跨中彎矩系數(shù);A為接觸面積;Ω為有效寬度范圍內(nèi)面積域。其中M0、KM0的計算式參見文獻[12]。

    圖10 均勻化過程示意Fig. 10 Sketch showing the homogenization process

    采用上述方法,計算有效寬度Be、平均接觸壓力ps和比值rp,結(jié)果如表4所示。由表4可知,考慮輪印載荷的非線性后,分布在加強筋上的均布輪印載荷較整個接觸面上的均布輪印載荷增加了1.279倍。

    表4 試驗工況均勻化結(jié)果Tab. 4 Homogenization of wheel load in experimental case

    2.2面板初始缺陷的影響

    2.1節(jié)有限元計算結(jié)果與1.4節(jié)試驗結(jié)果對比可見,4號測點和8號測點接觸壓力的相對誤差較大,而這兩個測點位于面板與加強筋的交接處,因此,認為誤差主要是由焊接加強筋引起的面板初始變形導致的。圖11給出了焊接完成后的試驗?zāi)P?,由圖11可見,加筋板的面板初始變形主要出現(xiàn)在加強筋與面板的交接焊縫處。由加筋板結(jié)構(gòu)焊接變形規(guī)律[18],認為加筋板面板初始變形形式如圖12所示,即采用半正弦初始變形形式。這里初始變形采用如下公式表達:

    其中,ω0為面板初始變形,B0為面板初始變形幅值,b為縱骨間距。根據(jù)實測值,面板初始變形幅值B0=0.8 mm。

    圖11 試驗?zāi)P虵ig. 11 Picture of experiment model

    圖12 加筋板面板初始變形形式Fig. 12 Initial deflection mode of the stiffened plate

    圖13為含面板缺陷的三維接觸壓力云圖。由圖13可見,由于面板初始變形的影響,導致加強筋處接觸壓力變大(較圖8)。圖14給出了試驗結(jié)果、不含缺陷有限元結(jié)果與含缺陷有限元結(jié)果的對比。由圖14可見,4號測點和8號測點的接觸壓力值變化較大,考慮初始變形后,4號測點計算值與試驗值的相對誤差由-19.63%變化為-6.40%,8號測點計算值與試驗值的相對誤差由-13.63%變化為8.47%,其余測點接觸壓力值變化不大。計算分布到加強筋上的均布輪印載荷幅值與接觸面內(nèi)的均布輪印載荷幅值的比值為1.365,與不考慮面板缺陷時的比值相差6.72%。

    圖13 含缺陷三維接觸壓力云圖Fig. 13 3D plot of contact pressure with initial deflection

    圖14 含缺陷有限元結(jié)果與試驗結(jié)果對比Fig. 14 Comparison of FEM results with initial deflection and experiment results

    提取結(jié)構(gòu)上部分關(guān)鍵測點的應(yīng)力值如表5所示,由表5可知,考慮面板的初始變形對結(jié)構(gòu)測點的應(yīng)力值影響不大。但考慮面板初始變形不僅影響加強筋處的接觸壓力值而且對分布到加強筋上的均布輪印載荷幅值也有較大的影響,因此在考慮輪印載荷分布問題時需要特別注意面板初始變形的影響。

    表5 含缺陷與不含缺陷應(yīng)力有限元結(jié)果Tab. 5 Stress results from FEM with & without initial deflection

    2.3加筋板結(jié)構(gòu)參數(shù)的影響

    本節(jié)通過數(shù)值仿真方法研究結(jié)構(gòu)參數(shù)對輪印載荷分布的影響,主要包括加強筋高度和面板板厚。計算工況選取為試驗工況,即輪印塊Y方向中心線與加強筋方向重合(x=97.5 mm直線),輪印塊X方向中心線與結(jié)構(gòu)縱對稱軸重合(y=0 mm直線),輪印塊長邊沿X方向。

    為了便于比較,將接觸壓力進行無量綱化[19],無量綱接觸壓力表示為p/pa。

    2.3.1 載荷大小

    在涉及結(jié)構(gòu)參數(shù)變化時,結(jié)構(gòu)所能承受的最大載荷將發(fā)生變化,而本文研究擬定各種工況下結(jié)構(gòu)最大應(yīng)力均處在210 MPa左右(與試驗應(yīng)力水平相同)。因此,需要先進行施加載荷大小的影響研究。以2.1節(jié)模型為研究對象,分別對加筋板結(jié)構(gòu)施加三種不同的載荷,即5 070 N、7 830 N和9 197 N,分別繪制y=0 mm直線和x=97.5 mm直線鋪板上的無量綱接觸壓力曲線如圖15所示。

    圖15 三種載荷下接觸壓力曲線Fig. 15 Comparison of contact pressure curves of the stiffened plate with three loads

    由圖15可見,同一工況不同載荷大小下輪印載荷的分布規(guī)律吻合。由文獻[12],該模型所能承受的輪印載荷為9 865 N,與試驗結(jié)果相當。所以可以認為結(jié)構(gòu)在彈性范圍內(nèi),在直升機甲板能承受的實際輪印載荷范圍內(nèi)施加載荷的大小對該無量綱接觸壓力分布無影響。

    2.3.2 加強筋高度

    根據(jù)文獻[12]中相對剛度系數(shù)的計算式:

    式中:γ1表示相對剛度系數(shù);a表示構(gòu)件跨距;t表示面板厚;b表示板格短邊;Is表示帶板縱骨或橫梁的慣性矩。

    通過改變加強筋高度h的方法改變Is來改變γ1,分別改變加強筋高度到60 mm和100 mm,使加強筋高度為35 mm、60 mm、100 mm模型的γ1分別為0.064 1、0.015 8和0.004 1。提取直線y=0 mm和直線x=97.5 mm上鋪板的無量綱接觸壓力曲線,如圖16所示。

    圖16 三種加強筋高度下接觸壓力曲線Fig. 16 Effect of different stiffener heights on contact pressure curves

    由圖16可見,接觸壓力的峰值均出現(xiàn)在輪印載荷塊的中心處。提取加筋板在y=0 mm直線上板格中心與加強筋處的節(jié)點位移,并計算兩者的相對位移(加強筋處位移減板格中心位移)。比較分析三種結(jié)構(gòu)的接觸壓力峰值以及輪印載荷的分布如表6所示。

    由表6可知,加強筋越高峰值無量綱接觸壓力和分布到加強筋上的均布輪印載荷越大,這主要是因為加強筋越高其慣性矩越大,使筋的變形越小,當加強筋高度由35 mm增加到60 mm時,加強筋處與板格中心處的相對位移由0.179 mm(位移差為正值主要是因為本工況輪印載荷中心線與加強筋重合,板格中心為載荷邊界位置,即考慮整個板架而言加強筋處為載荷中心,而加強筋的高度較小,對加強筋處板格支撐有限,從而導致加強筋處位移大于板格中心位移)變?yōu)?0.033 mm,而由圖16可以看出載荷更多的往加強筋處分布,使得該處無量綱峰值壓力增加了5.96%,載荷比值增加了5.47%。因此,增加加強筋的高度使輪印載荷更多的分布在加強筋上,但隨著高度的增加,載荷增加速度降低。

    表6 不同加強筋高度下輪印載荷分布Tab. 6 Effect of different stiffener heights on wheel load distribution

    2.3.3 面板厚度

    以試驗工況模型為參考,根據(jù)2.3.2節(jié)中的相對剛度計算方法,通過減小板厚的方法來降低結(jié)構(gòu)的相對剛度,分別減小板厚到3.5 mm和2.1 mm,使得兩個新模型的相對剛度分別為0.015 8和0.004 1。然后提取直線y=0 mm和直線x=97.5 mm上鋪板的無量綱接觸壓力曲線,如圖17所示。

    圖17 三種板厚下接觸壓力曲線Fig. 17 Effect of different plate thicknesses on contact pressure curves

    由圖17可見,隨著面板厚度的減小y=0 mm直線上鋪板接觸壓力曲線變得陡峭,x=97.5 mm直線上鋪板接觸壓力變大。取加筋板在y=0直線上板格中心與加強筋處的節(jié)點位移,并計算兩者的相對位移(加強筋處位移減板格中心位移)。比較分析三種結(jié)構(gòu)的接觸壓力峰值如表7所示。

    表7 不同面板厚度下輪印載荷分布Tab. 7 Effect of different plate thicknesses on wheel load distribution

    圖18 相對剛度與輪印載荷幅值比曲線Fig. 18 Relative stiffness verse amplitude ratio curves

    由表7可知,當面板厚度為6 mm時,相對位移為0.179 mm(同表6中H-35工況),而當面板厚度減小到3.5 mm時,板的剛度減小,板和加強筋的相對剛度減小,加強筋的支撐作用變強,從而導致相對位移變?yōu)樨撝?,且面板越薄,位移差越大。同時峰值無量綱接觸壓力和分布到加強筋上的輪印載荷越多。因此,減小面板厚度使輪印載荷更多的分布在加強筋上,且隨著面板厚度的減小,載荷增加速度增大。

    綜合表6和表7,以減小面板厚度和增加加強筋高度的形式來改變板和加強筋的相對剛度,得到在相對剛度變化相同的情況下輪印載荷分布到加強筋上的均布載荷幅值的比值,繪制比值隨相對剛度的變化曲線如圖18所示。由圖18可見,在相對剛度變化相同的情況下,改變面板厚度和改變加強筋高度對輪印載荷的分布影響不同,隨著板厚的降低輪印載荷分布到加強筋上的比例增加越多,而隨著加強筋高度的增加輪印載荷分布到加強筋上的比例增加越少。

    3 結(jié) 語

    采用試驗和有限元計算的方法,對輪印載荷在加筋板上的分布進行了研究。通過試驗方法直接測量了輪印接觸面積內(nèi)的接觸壓力以及結(jié)構(gòu)上的應(yīng)力。通過數(shù)值計算的方法分析了面板初始缺陷、加強筋高度以及加筋板面板厚度對輪印載荷分布的影響,得到以下結(jié)論:

    1) 有限元結(jié)果與試驗結(jié)果吻合較好,提出的數(shù)值仿真模型合理;

    2) 在試驗工況下,分布在加強筋上均布輪印載荷為接觸面上均布輪印載荷的1.279倍;

    3) 在接觸面內(nèi),面板的初始缺陷對接觸壓力的分布有很大影響,在進行輪印載荷分布試驗時需要特別注意初始缺陷的影響;

    4) 在直升機甲板能承受的實際輪印載荷范圍內(nèi),輪印載荷的大小對輪印載荷分布的影響可以忽略;

    5) 在相對剛度改變相同的情況下,改變面板厚度和改變加強筋高度對輪印載荷分布的影響不同。

    [1] GANGARAO H V, VALI A. Truck-tire steel grid deck contact pressure distributions [J]. Journal of Structural Engineering, 1990, 116(3): 791-808.

    [2] LUO R, EDLUND B. Ultimate strength of girders with trapezoidally corrugated webs under patch loading[J]. Thin-Walled Structures, 1996, 24(2): 135-156.

    [3] 楊坤榮.FPSO的直升機平臺設(shè)計及建造研究[D].廣州:華南理工大學,2015.(YANG Kunrong. Helideck design and construction for FPSO[D]. Guangzhou: South China University of Technology, 2015. (in Chinese))

    [4] 張玉龍, 程遠勝, 曾廣武, 等. 輪印載荷作用下波紋型夾層板格強度特性分析[J]. 中國造船, 2015, 56(1): 69-79.(ZHANG Yulong, CHENG Yuansheng, ZENG Guangwu, et al. Analysis on strength properties of corrugated-core sandwich plates under patch loading[J]. Ship Building of China, 2015, 56(1): 69-79.(in Chinese))

    [5] 蔣志巖, 古長江. 滾裝船車輛甲板強度分析[J]. 船舶, 2004(2): 20-23.(JIANG Zhiyan, GU Changjiang. Strength analysis for the vehicle decks of Ro-Ro ship[J]. Ship & Boat, 2004(2): 20-23.( in Chinese))

    [6] 鐘悅熾. 承受車輪載荷的滾裝通道設(shè)備結(jié)構(gòu)設(shè)計(下)[J]. 船舶, 1991(4): 36-44.(ZHONG Yuechi. Structural design of Ro-Ro ship channel devices under patch loading(lower)[J]. Ship & Boat, 1991(4): 36-44.( in Chinese))

    [7] 曾鳴. 登陸艦車輛甲板結(jié)構(gòu)設(shè)計和強度校核規(guī)范建立的研究[D]. 上海:上海交通大學, 2007.(ZENG Ming. Study of the structure design and the establishment of the strength criterion of the vehicle deck of the landing craft[D]. Shanghai: Shanghai Jiao Tong University, 2007. (in Chinese))

    [8] 曾鳴, 王福花, 王德禹. 氣墊船著陸墊與車輛甲板之間的非線性接觸分析[J]. 中國艦船研究, 2006, 1(3): 58-61.(ZENG Ming, WANG Fuhua, WANG Deyu. The non-linear analysis of the contact problem between the hovercraft cushion and the vehicle decks[J].Chinese Journal of Ship Research, 2006, 1(3):58-61.(in Chinese))

    [9] 彭興寧, 徐向東. 輪載作用下的甲板結(jié)構(gòu)設(shè)計[J]. 船舶工程, 2000(2): 26-29.(PENG Xingning, XU Xiangdong. Design of deck-platform by wheel load actions[J]. Ship Engineering, 2000(2): 26-29. (in Chinese))

    [10] ROMANOFF J, VARSTA P, REMES H. Laser-welded web-core sandwich plates under patch loading[J]. Marine Structures, 2007, 20(1): 25-48.

    [11] JACKSON R I. Design of deck structures under wheel loads[J]. Naval Architect, 1981(3): T119-T114.

    [12] STAINBACK J. Structural analysis of helicopter flight and hangar decks[R]. Naval Surface Warfare Center Carderock Division, Bethesda MD, 2001.

    [13] SALEEM M A, MIRMIRAN A, XIA J, et al. Experimental evaluation of aluminum bridge deck system[J]. Journal of Bridge Engineering, 2010, 17(1): 97-106.

    [14] BROWN D L, BERMAN J W. Fatigue and strength evaluation of two glass fiber-reinforced polymer bridge decks[J]. Journal of Bridge Engineering, 2010, 15(3): 290-301.

    [15] WANG H, AL-QADI I L, STANCIULESCU I. Effect of surface friction on tire-pavement contact stresses during vehicle maneuvering[J]. Journal of Engineering Mechanics, 2013, 140(4): 04014001.

    [16] 樓京俊, 朱石堅, 唐斯密. 在彈性體有限元的M—R模型中材料參數(shù)的逆向推算[J]. 噪聲與振動控制, 2012, 32(1): 29-33.(LOU Jingjun, ZHU Shijian, TANG Simi. Study on reverse deduction of two material parameters in mooney-revlin model of rubber isolator[J]. Noise and Vibration Control, 2012, 32(1): 29-33.(in Chinese))

    [17] 王吉忠, 莊繼德, 李日春. 輪胎胎面橡膠塊與剛性路面摩擦數(shù)值分析[J]. 農(nóng)業(yè)工程學報, 1998, 14(2): 104-108.(WANG Jizhong, ZHUANG Jide, LI Richun. Numerical analysis for friction contact of tread rubber-block with rigid road surface[J]. Transactions of the Chinese Society of Agricultural Engineering, 1998,14(2): 104-108.(in Chinese))

    [18] 陳哲超, 陳震. 加筋板結(jié)構(gòu)連續(xù)焊焊接變形規(guī)律[J]. 中國艦船研究, 2016, 11(3): 68-73.(CHEN Zhenchao, CHEN Zhen. Deformation rules of stiffened plates with continuous welding[J].Chinese Journal of Ship Research, 2016, 11(3):68-73.(in Chinese))

    [19] MAJUMDAR P K, LESKO J J, COUSINS T E, et al. Conformable tire patch loading for FRP composite bridge deck [J]. Journal of Composites for Construction, 2009, 13(6): 575-81.

    Experiment and numerical analyses on wheel load distribution on stiffened plate

    LIU Cong, CHENG Yuansheng, ZHANG Pan, LIU Jun

    (School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China)

    U663.6

    A

    10.16483/j.issn.1005-9865.2017.04.011

    1005-9865(2017)04-0084-10

    2016-07-31

    劉 聰(1991-),男,湖北京山人,碩士研究生,主要從事結(jié)構(gòu)設(shè)計與分析方面的研究。E-mail: hbjmlc@163.com

    程遠勝。E-mail: yscheng@hust.edu.cn

    猜你喜歡
    有限元變形
    談詩的變形
    中華詩詞(2020年1期)2020-09-21 09:24:52
    新型有機玻璃在站臺門的應(yīng)用及有限元分析
    基于有限元的深孔鏜削仿真及分析
    基于有限元模型對踝模擬扭傷機制的探討
    “我”的變形計
    變形巧算
    例談拼圖與整式變形
    會變形的餅
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    两性夫妻黄色片| 亚洲精品中文字幕在线视频| 日韩精品青青久久久久久| 黄色女人牲交| 免费搜索国产男女视频| e午夜精品久久久久久久| 久久久久久久久免费视频了| 9热在线视频观看99| 久久久久亚洲av毛片大全| 欧美午夜高清在线| 久久久久九九精品影院| 伦理电影免费视频| 夜夜看夜夜爽夜夜摸 | 国产精品成人在线| 淫妇啪啪啪对白视频| 一级a爱视频在线免费观看| 婷婷精品国产亚洲av在线| 亚洲av成人av| 男人的好看免费观看在线视频 | 黄色女人牲交| 男人的好看免费观看在线视频 | 亚洲av电影在线进入| 欧美色视频一区免费| av天堂久久9| 欧美日韩av久久| 丁香欧美五月| 亚洲国产精品一区二区三区在线| 亚洲精品国产一区二区精华液| 午夜福利影视在线免费观看| 久久国产乱子伦精品免费另类| 久久久精品欧美日韩精品| 中文字幕人妻熟女乱码| 亚洲国产欧美网| 成人av一区二区三区在线看| 国产男靠女视频免费网站| 在线观看www视频免费| 一级片'在线观看视频| 级片在线观看| 9色porny在线观看| 亚洲专区字幕在线| 亚洲成人久久性| 在线观看一区二区三区| 国内毛片毛片毛片毛片毛片| 别揉我奶头~嗯~啊~动态视频| av天堂久久9| 久久精品亚洲精品国产色婷小说| 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 乱人伦中国视频| 天堂动漫精品| 欧美激情高清一区二区三区| 久久精品成人免费网站| 69av精品久久久久久| 久久久久久久久久久久大奶| 99在线视频只有这里精品首页| 国产区一区二久久| 最近最新中文字幕大全电影3 | 男人舔女人下体高潮全视频| 亚洲熟妇熟女久久| 国产伦人伦偷精品视频| 这个男人来自地球电影免费观看| 首页视频小说图片口味搜索| 伦理电影免费视频| 在线观看免费高清a一片| 99国产精品一区二区三区| 欧美日韩黄片免| 国产亚洲精品第一综合不卡| 久久国产精品男人的天堂亚洲| 国产91精品成人一区二区三区| 国产免费现黄频在线看| 久久精品aⅴ一区二区三区四区| 激情在线观看视频在线高清| 亚洲专区中文字幕在线| www.999成人在线观看| 在线看a的网站| 如日韩欧美国产精品一区二区三区| 中出人妻视频一区二区| 国内久久婷婷六月综合欲色啪| 黄色怎么调成土黄色| 中文字幕av电影在线播放| 国产亚洲精品综合一区在线观看 | 亚洲欧美一区二区三区久久| 亚洲 欧美 日韩 在线 免费| 精品久久久久久,| 欧美+亚洲+日韩+国产| 午夜91福利影院| 亚洲人成网站在线播放欧美日韩| 美女国产高潮福利片在线看| 亚洲美女黄片视频| 在线免费观看的www视频| 51午夜福利影视在线观看| 成年人黄色毛片网站| 首页视频小说图片口味搜索| 日韩大尺度精品在线看网址 | 国产精品永久免费网站| 一区二区三区国产精品乱码| 亚洲精品成人av观看孕妇| 色在线成人网| 国产三级黄色录像| 免费搜索国产男女视频| 女人精品久久久久毛片| 欧美久久黑人一区二区| 欧美国产精品va在线观看不卡| 国产成人欧美| 久热这里只有精品99| 国产有黄有色有爽视频| 波多野结衣一区麻豆| 法律面前人人平等表现在哪些方面| 亚洲国产精品sss在线观看 | 亚洲色图综合在线观看| 亚洲成人免费av在线播放| 亚洲精品中文字幕一二三四区| 校园春色视频在线观看| 身体一侧抽搐| 一夜夜www| 免费久久久久久久精品成人欧美视频| 精品一区二区三区视频在线观看免费 | 日韩欧美一区二区三区在线观看| 香蕉国产在线看| 真人做人爱边吃奶动态| avwww免费| 亚洲人成电影观看| 看黄色毛片网站| 免费看a级黄色片| 久久人妻熟女aⅴ| 国产精品 欧美亚洲| 久久精品91蜜桃| 国产亚洲欧美在线一区二区| 一级片'在线观看视频| 女同久久另类99精品国产91| 99国产精品免费福利视频| 一区在线观看完整版| 男男h啪啪无遮挡| 天天躁狠狠躁夜夜躁狠狠躁| 大香蕉久久成人网| 午夜91福利影院| 久久狼人影院| 麻豆国产av国片精品| 丝袜美足系列| 久久狼人影院| 一进一出抽搐动态| 欧美日韩精品网址| 国产黄a三级三级三级人| 免费在线观看亚洲国产| 变态另类成人亚洲欧美熟女 | 成年人黄色毛片网站| 在线国产一区二区在线| 两个人免费观看高清视频| 久久九九热精品免费| 久久久久久久久中文| 欧美av亚洲av综合av国产av| 午夜精品在线福利| 好看av亚洲va欧美ⅴa在| 精品一区二区三卡| 麻豆国产av国片精品| 女人爽到高潮嗷嗷叫在线视频| 国产xxxxx性猛交| 视频在线观看一区二区三区| av国产精品久久久久影院| 男女午夜视频在线观看| 一级片免费观看大全| 久久人妻av系列| 亚洲 国产 在线| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久免费视频 | 欧美成狂野欧美在线观看| 男男h啪啪无遮挡| 国产不卡一卡二| 国产精品成人在线| 无限看片的www在线观看| 久久精品aⅴ一区二区三区四区| 男人操女人黄网站| 国产亚洲精品一区二区www| 国产蜜桃级精品一区二区三区| 香蕉国产在线看| 咕卡用的链子| 国产成年人精品一区二区 | 日韩 欧美 亚洲 中文字幕| 亚洲五月婷婷丁香| 男男h啪啪无遮挡| 午夜a级毛片| 一区在线观看完整版| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 精品久久久久久电影网| 精品一区二区三区视频在线观看免费 | 日韩有码中文字幕| 激情在线观看视频在线高清| 国产精品日韩av在线免费观看 | www日本在线高清视频| 18禁国产床啪视频网站| 色婷婷av一区二区三区视频| 久久精品91蜜桃| 亚洲九九香蕉| 人妻丰满熟妇av一区二区三区| 波多野结衣av一区二区av| 黄片大片在线免费观看| 不卡av一区二区三区| 在线免费观看的www视频| 桃色一区二区三区在线观看| 欧美精品亚洲一区二区| 久久人妻福利社区极品人妻图片| 亚洲精品国产色婷婷电影| 三上悠亚av全集在线观看| 久久99一区二区三区| 亚洲在线自拍视频| avwww免费| 99热国产这里只有精品6| 国产精品国产高清国产av| 天堂√8在线中文| 精品久久蜜臀av无| 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| 美国免费a级毛片| 亚洲精品一二三| 亚洲免费av在线视频| 亚洲三区欧美一区| 久久99一区二区三区| 日日夜夜操网爽| 国产主播在线观看一区二区| 性少妇av在线| 日韩中文字幕欧美一区二区| 久久精品国产亚洲av高清一级| 国产在线精品亚洲第一网站| 精品人妻在线不人妻| 99久久国产精品久久久| 成人特级黄色片久久久久久久| 免费高清在线观看日韩| 日韩大尺度精品在线看网址 | 国产精品久久视频播放| 老司机亚洲免费影院| 日韩三级视频一区二区三区| 天堂影院成人在线观看| 欧美在线黄色| 国产国语露脸激情在线看| 国产乱人伦免费视频| 国产一区二区激情短视频| 日韩欧美在线二视频| 黑人巨大精品欧美一区二区蜜桃| 国产色视频综合| 国产精品乱码一区二三区的特点 | 午夜老司机福利片| 国产无遮挡羞羞视频在线观看| 国产亚洲精品久久久久5区| 人人澡人人妻人| 国产在线观看jvid| 亚洲成a人片在线一区二区| 精品国产超薄肉色丝袜足j| 天天添夜夜摸| 免费一级毛片在线播放高清视频 | www.www免费av| 高清欧美精品videossex| 日本黄色日本黄色录像| 亚洲熟女毛片儿| 亚洲av第一区精品v没综合| 精品国内亚洲2022精品成人| 亚洲激情在线av| 桃色一区二区三区在线观看| 日韩欧美国产一区二区入口| 别揉我奶头~嗯~啊~动态视频| 国产精品98久久久久久宅男小说| 俄罗斯特黄特色一大片| 亚洲av片天天在线观看| 亚洲熟妇熟女久久| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站| 国产亚洲av高清不卡| 精品免费久久久久久久清纯| 看片在线看免费视频| 18禁国产床啪视频网站| 高清在线国产一区| 国产不卡一卡二| 桃色一区二区三区在线观看| 人人妻,人人澡人人爽秒播| 大香蕉久久成人网| 亚洲欧美激情在线| 大型av网站在线播放| 美女福利国产在线| 深夜精品福利| 中文字幕最新亚洲高清| 国产蜜桃级精品一区二区三区| 国产单亲对白刺激| 久久久国产欧美日韩av| 午夜老司机福利片| 久久人人97超碰香蕉20202| 超色免费av| 亚洲一区中文字幕在线| 日韩欧美免费精品| netflix在线观看网站| 亚洲熟妇中文字幕五十中出 | 亚洲中文av在线| 丰满迷人的少妇在线观看| 亚洲精品久久成人aⅴ小说| 首页视频小说图片口味搜索| 久久人人精品亚洲av| 操出白浆在线播放| 国产成人欧美| 波多野结衣一区麻豆| 国产精品亚洲一级av第二区| 久久精品成人免费网站| 一级黄色大片毛片| 麻豆av在线久日| 欧美 亚洲 国产 日韩一| 国产免费现黄频在线看| 久久国产亚洲av麻豆专区| 人人妻人人爽人人添夜夜欢视频| 最新美女视频免费是黄的| 国产高清激情床上av| 在线av久久热| 色在线成人网| 久久精品国产99精品国产亚洲性色 | 黑人巨大精品欧美一区二区mp4| а√天堂www在线а√下载| 韩国av一区二区三区四区| 欧美激情极品国产一区二区三区| 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 黑人猛操日本美女一级片| 丁香欧美五月| 老司机福利观看| 波多野结衣av一区二区av| 久久这里只有精品19| 国产欧美日韩一区二区三区在线| 日韩有码中文字幕| 好男人电影高清在线观看| 欧美乱妇无乱码| 午夜两性在线视频| 亚洲人成电影观看| 动漫黄色视频在线观看| 午夜精品在线福利| 悠悠久久av| 嫁个100分男人电影在线观看| 9热在线视频观看99| 亚洲欧美一区二区三区黑人| 欧美激情久久久久久爽电影 | 免费久久久久久久精品成人欧美视频| avwww免费| 亚洲欧美精品综合久久99| 亚洲精品一区av在线观看| 午夜老司机福利片| 精品人妻1区二区| 国产成人一区二区三区免费视频网站| 久久欧美精品欧美久久欧美| 老司机午夜十八禁免费视频| 久久欧美精品欧美久久欧美| 两人在一起打扑克的视频| 乱人伦中国视频| 国产成人一区二区三区免费视频网站| 婷婷精品国产亚洲av在线| 久久久精品国产亚洲av高清涩受| 婷婷精品国产亚洲av在线| 美女 人体艺术 gogo| 精品电影一区二区在线| x7x7x7水蜜桃| 12—13女人毛片做爰片一| 香蕉丝袜av| 国产成人欧美在线观看| 欧美人与性动交α欧美精品济南到| 91麻豆av在线| 国产一区二区三区视频了| 久久久久久久午夜电影 | 日韩免费av在线播放| 在线观看午夜福利视频| 国产伦一二天堂av在线观看| 精品国产美女av久久久久小说| 中文字幕人妻熟女乱码| 老司机在亚洲福利影院| 亚洲熟妇熟女久久| 熟女少妇亚洲综合色aaa.| 日本免费一区二区三区高清不卡 | 亚洲色图 男人天堂 中文字幕| 老司机深夜福利视频在线观看| 亚洲熟妇熟女久久| 亚洲国产中文字幕在线视频| 精品电影一区二区在线| 麻豆av在线久日| 亚洲一区二区三区不卡视频| 美女高潮喷水抽搐中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 12—13女人毛片做爰片一| 老鸭窝网址在线观看| 亚洲av成人一区二区三| 18禁国产床啪视频网站| 两人在一起打扑克的视频| 国产男靠女视频免费网站| 免费在线观看亚洲国产| 国产在线观看jvid| 免费av中文字幕在线| 久久精品国产清高在天天线| 人妻丰满熟妇av一区二区三区| 老汉色av国产亚洲站长工具| 一区二区三区国产精品乱码| 日韩欧美免费精品| 欧美日韩精品网址| 91国产中文字幕| 精品一区二区三区视频在线观看免费 | 中文字幕最新亚洲高清| 99国产精品免费福利视频| 亚洲 欧美一区二区三区| 国产亚洲欧美在线一区二区| 一本综合久久免费| 久久亚洲真实| 妹子高潮喷水视频| 国产欧美日韩综合在线一区二区| 午夜影院日韩av| 国产成人欧美在线观看| 久久精品国产亚洲av香蕉五月| 成人影院久久| 99国产精品免费福利视频| 可以在线观看毛片的网站| av有码第一页| 又黄又粗又硬又大视频| 99国产精品一区二区蜜桃av| 在线观看66精品国产| 国产色视频综合| 国产精品综合久久久久久久免费 | 午夜久久久在线观看| 国产欧美日韩一区二区三| 黄片播放在线免费| 久久香蕉国产精品| 国产成+人综合+亚洲专区| 国产成人精品久久二区二区91| 老熟妇乱子伦视频在线观看| 久久精品亚洲av国产电影网| 久久人人97超碰香蕉20202| e午夜精品久久久久久久| 巨乳人妻的诱惑在线观看| 午夜久久久在线观看| 国产一卡二卡三卡精品| 免费观看精品视频网站| 国产精品爽爽va在线观看网站 | 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯| 黑人巨大精品欧美一区二区mp4| 久久午夜亚洲精品久久| 国产深夜福利视频在线观看| 久久久久国产一级毛片高清牌| 国产精品香港三级国产av潘金莲| 啪啪无遮挡十八禁网站| 99国产精品一区二区三区| 久久久久久久久久久久大奶| 精品少妇一区二区三区视频日本电影| 久久久国产成人免费| 人人妻,人人澡人人爽秒播| 久久久久国产一级毛片高清牌| 久热这里只有精品99| 久久久久久久久免费视频了| 亚洲人成77777在线视频| 99国产精品一区二区三区| a级毛片在线看网站| 精品一区二区三区视频在线观看免费 | 亚洲色图综合在线观看| 欧美激情高清一区二区三区| 在线观看66精品国产| 后天国语完整版免费观看| 波多野结衣高清无吗| 亚洲第一欧美日韩一区二区三区| 国产精品久久久人人做人人爽| 黄片播放在线免费| 99精国产麻豆久久婷婷| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| 99久久人妻综合| 成熟少妇高潮喷水视频| 精品乱码久久久久久99久播| 一级毛片高清免费大全| 欧美成人性av电影在线观看| 日本黄色视频三级网站网址| 黄色片一级片一级黄色片| 国产亚洲精品第一综合不卡| 国产片内射在线| 欧美激情极品国产一区二区三区| 99国产精品一区二区三区| www.www免费av| 日韩免费av在线播放| www国产在线视频色| 丰满的人妻完整版| 真人做人爱边吃奶动态| 大型av网站在线播放| 动漫黄色视频在线观看| 夜夜夜夜夜久久久久| 80岁老熟妇乱子伦牲交| 一级毛片精品| 精品免费久久久久久久清纯| 天堂俺去俺来也www色官网| 免费在线观看黄色视频的| 在线观看一区二区三区激情| 久久人人爽av亚洲精品天堂| 亚洲av电影在线进入| 曰老女人黄片| 99精品久久久久人妻精品| 亚洲欧美日韩无卡精品| 久久人妻av系列| 亚洲人成77777在线视频| 亚洲色图av天堂| 后天国语完整版免费观看| 午夜影院日韩av| 精品久久久久久成人av| 啦啦啦 在线观看视频| 欧美 亚洲 国产 日韩一| av超薄肉色丝袜交足视频| 热re99久久精品国产66热6| 免费不卡黄色视频| 1024香蕉在线观看| 亚洲欧美一区二区三区久久| 久久久久国内视频| 日本五十路高清| 高清在线国产一区| 国产精品爽爽va在线观看网站 | 久久影院123| 精品久久久久久久毛片微露脸| 777久久人妻少妇嫩草av网站| 国产熟女午夜一区二区三区| 国产av在哪里看| 国产高清videossex| 美女高潮到喷水免费观看| 亚洲一区中文字幕在线| 中文字幕人妻熟女乱码| 在线观看一区二区三区| 中文字幕av电影在线播放| 国产成人精品无人区| 妹子高潮喷水视频| 精品久久蜜臀av无| 老汉色av国产亚洲站长工具| 欧美日韩瑟瑟在线播放| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 久久国产精品人妻蜜桃| 亚洲九九香蕉| 又紧又爽又黄一区二区| 国产精品偷伦视频观看了| 天堂俺去俺来也www色官网| 极品人妻少妇av视频| 久久亚洲精品不卡| 高清在线国产一区| 一级a爱片免费观看的视频| 大码成人一级视频| 婷婷六月久久综合丁香| 最近最新中文字幕大全电影3 | 亚洲欧美激情在线| 电影成人av| 夜夜看夜夜爽夜夜摸 | 久久人妻av系列| 99热只有精品国产| 最近最新中文字幕大全电影3 | 欧美久久黑人一区二区| 日本精品一区二区三区蜜桃| 黄色a级毛片大全视频| 欧美黄色片欧美黄色片| 91大片在线观看| 天堂动漫精品| 最近最新中文字幕大全电影3 | 一级片'在线观看视频| 久久精品aⅴ一区二区三区四区| 黄色丝袜av网址大全| 欧美老熟妇乱子伦牲交| 亚洲情色 制服丝袜| 国产一区二区三区综合在线观看| 国产精品野战在线观看 | 精品久久久精品久久久| 亚洲精品中文字幕一二三四区| 精品国产亚洲在线| 自拍欧美九色日韩亚洲蝌蚪91| 日韩三级视频一区二区三区| av免费在线观看网站| 日韩一卡2卡3卡4卡2021年| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 女生性感内裤真人,穿戴方法视频| 中文字幕最新亚洲高清| 久久精品亚洲熟妇少妇任你| 国产免费现黄频在线看| 91在线观看av| 国产有黄有色有爽视频| 777久久人妻少妇嫩草av网站| 18禁裸乳无遮挡免费网站照片 | 久久精品国产亚洲av香蕉五月| 午夜久久久在线观看| avwww免费| 日本vs欧美在线观看视频| 国产不卡一卡二| 级片在线观看| 亚洲熟妇熟女久久| 免费在线观看日本一区| 99精品欧美一区二区三区四区| 国产日韩一区二区三区精品不卡| 国产精品九九99| 免费在线观看视频国产中文字幕亚洲| 电影成人av| 免费在线观看影片大全网站| 日本wwww免费看| 在线视频色国产色| 级片在线观看| 亚洲久久久国产精品| netflix在线观看网站| 日韩av在线大香蕉| 中文亚洲av片在线观看爽| 亚洲五月婷婷丁香| 久久久久亚洲av毛片大全| 欧美一级毛片孕妇| 黄片小视频在线播放| 欧美最黄视频在线播放免费 | 又黄又粗又硬又大视频| 欧美日韩一级在线毛片| 国产精品久久久人人做人人爽| 中文字幕精品免费在线观看视频| 亚洲人成电影免费在线| 国产免费av片在线观看野外av| 高清欧美精品videossex| 久久久国产成人免费| 国产欧美日韩一区二区三| 在线十欧美十亚洲十日本专区| 一进一出抽搐动态| 人人妻人人爽人人添夜夜欢视频|