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

    大型復雜結構安定性數值分析的新方法及其應用

    2016-06-23 08:39:29鄒宗園郭寶峰毛春燕
    中國機械工程 2016年3期
    關鍵詞:安定性彈塑性液壓缸

    鄒宗園 郭寶峰 金 淼 毛春燕

    1.先進鍛壓成形技術與科學教育部重點實驗室(燕山大學),秦皇島,0660042.太原重工股份有限公司,太原,030024

    大型復雜結構安定性數值分析的新方法及其應用

    鄒宗園1郭寶峰1金淼1毛春燕2

    1.先進鍛壓成形技術與科學教育部重點實驗室(燕山大學),秦皇島,0660042.太原重工股份有限公司,太原,030024

    摘要:大型復雜結構承載后普遍存在局部高應力區(qū),特別是某些巨型結構,由于受幾何尺寸和制造能力的限制,其局部高應力往往會接近甚至超過材料的屈服極限,因而成為彈性范圍內的設計難題。應用安定性分析設計準則,則可通過利用材料的塑性承載能力很好地解決這類問題。為此,在大型復雜結構的強度分析中引入安定性理論,應用Melan安定下限定理,通過將峰值載荷下的彈塑性應力場與彈性參考應力場相減,構造定理要求的自平衡應力場,并推導了相應的安定性判定條件。在此基礎上,結合彈塑性增量有限元技術,通過不同載荷水平下的加載-卸載計算,計算相應載荷水平下的殘余應力場。根據加載-卸載得到的殘余應力場與構造的自平衡場的關系,按照安定性判定條件確定安定極限載荷,從而建立了一種適合于大型復雜結構安定性數值分析的新方法,并應用該方法對某125 MN鍛造液壓機的主液壓缸缸體結構進行了安定性分析。

    關鍵詞:彈塑性;安定性;液壓缸;安定極限載荷;局部高應力

    0引言

    在工業(yè)制造領域,一些大型復雜承載結構,因其構件尺寸大、結構復雜,承載后應力分布不均勻,即使經過結構優(yōu)化設計也不可避免地在結構局部形狀和尺寸發(fā)生突變的位置產生應力集中,形成高應力區(qū)。這些局部高應力往往是結構中平均應力的2~5倍,常常接近甚至超過材料的屈服極限。例如,大型液壓機本體結構中的法蘭支撐式液壓缸,其法蘭根部圓角區(qū)域局部受力狀態(tài)復雜,且存在應力集中而出現局部高應力,所以其法蘭與缸筒的過渡部位一直是這種結構強度分析校核的關鍵控制部位,即使采用了經過優(yōu)化設計的過渡型線也無法避免局部高應力現象的發(fā)生。

    根據安定理論,對于彈性安定的結構,在初始承載循環(huán)后會產生局部塑性變形,但在結構安定之后便不再有新的塑性變形產生,結構表現為完全彈性響應。也就是說,即使在局部高應力區(qū)產生了塑性變形,結構也可能進入安定狀態(tài),并不會發(fā)生彈塑性強度破壞。安定理論合理地利用了材料的塑性承載能力,無需再將結構的最大應力控制在材料許用應力之內。

    顯然,若引入安定性理念,以結構的安定極限載荷作為結構的彈塑性承載能力界限,那么對于存在局部高應力區(qū)的大型復雜承載結構而言,就意味著能夠更充分地利用材料的彈塑性承載能力?,F階段,基于安定理論的彈塑性強度設計準則已經被寫入ASME Ⅷ-2[1]和EN-13445[2]等壓力容器建造規(guī)范。除了在壓力容器[3-4]的設計校核中得到應用外,也在核電工業(yè)構件[5-6]、航空航天構件[7]、鐵路道軌[8]、土工地基[9]、近海結構[10]及鋼結構框架[11]的分析設計中得到了應用。為促進安定理論的應用,已發(fā)展了許多安定性數值分析方法。這些方法包括直接循環(huán)法和基于安定定理的數值計算方法,如數學規(guī)劃法、彈性補償法等。然而直接循環(huán)法需進行不同載荷下的多次循環(huán)計算,繁瑣耗時[12];其他數值計算方法則需通過編制子程序,嵌入商業(yè)有限元軟件中實現分析[13-15],各有優(yōu)勢但又存在其應用的局限性,尚無法簡便有效地實現液壓機本體結構的安定性分析。

    本文針對存在局部高應力的大型復雜承載結構的強度分析問題,引入安定理論并根據Melan安定下限定理建立了一種應用彈塑性增量有限元技術求解結構安定極限載荷的數值分析新方法,并應用該方法對125 MN鍛造液壓機主缸進行了安定性分析。

    1安定性數值分析方法

    1.1安定性判定條件

    Melan安定定理描述為,對于理想彈塑性結構,如果存在一個時間無關的自平衡應力場,使得結構內的任何點在加載域內的所有載荷下,其與假設結構為理想彈性體時的彈性參考應力場疊加后得到的應力場滿足如下條件則結構安定:

    (1)

    根據上述定理,可以通過構造時間無關的自平衡應力場來求解滿足安定下限定理的最大外載荷,即安定極限載荷。

    液壓機本體結構件受到工作載荷的循環(huán)加載-卸載作用。液壓缸的循環(huán)工作載荷可表示為

    P(t)=P(t+T)

    (2)

    式中,T為載荷循環(huán)周期;P(t)∈[0,P]為載荷域。

    (3)

    由于載荷域P(t)∈[0,P]為一維加載域,且彈性參考應力是與載荷成正比的,所以在整個載荷域內的安定性可以通過載荷域兩個頂點處即谷值、峰值處進行判斷即可。因此若要滿足安定下限定理的安定條件(式(1)),則在循環(huán)載荷谷值、峰值處均需滿足式(1),即

    (4)

    (5)

    將式(3)代入式(5)可得

    (6)

    1.2數值方法的實現

    圖1 安定性數值分析方法的加載路徑

    對于圖1所示的峰值漸增的加載路徑,第一個載荷峰值一般可選為結構的設計載荷,載荷峰值的增幅ΔP則可根據計算精度要求進行確定。

    圖2 安定極限載荷的確定

    由于結構的安定性取決于其強度最薄弱部位的關鍵點,所以安定性數值分析只需監(jiān)測其關鍵點。那么,對于理想彈塑性材料結構,關鍵點的等效殘余應力σr首次等于材料的屈服極限σs時所對應的載荷Ps為安定極限載荷,如圖2所示。因此,通過從有限元分析結果中提取結構關鍵點在不同載荷峰值下加載-卸載結束時應力狀態(tài)的等效應力就可方便地確定結構的安定極限載荷。圖3所示為上述安定性數值分析流程。

    圖3 安定性有限元數值分析流程圖

    2液壓缸安定性分析

    2.1有限元模型

    圖4是某125 MN自由鍛造液壓機主液壓缸的缸體結構簡圖。該液壓缸為法蘭支撐結構,支撐法蘭位于缸口。缸體按照彈性范圍的強度設計準則設計完成。液壓缸的額定工作壓力為35 MPa。缸體材質為鍛鋼35,其屈服極限不小于235 MPa。缸體的有限元彈性強度分析結果表明,在缸筒與法蘭的圓角過渡區(qū)域存在局部高應力,最高應力值已接近材料的屈服極限,那么,根據彈性強度設計準則,其安全系數(材料屈服極限與最高應力的比值)略大于1,因此缸體的彈性強度安全裕度非常有限。但在液壓機投產以來3年多的運行實踐中,并未出現因液壓缸缸體強度不足而導致的異常情況。有鑒于此,應用上述安定性數值分析方法對該液壓缸缸體進行了彈塑性強度分析。

    (a)缸體簡圖

    (b)局部放大圖Ⅰ(c)局部放大圖Ⅱ圖4 液壓缸缸體結構簡圖

    由于螺栓孔及橫梁的彎曲變形對缸體高應力區(qū)影響很小,所以在建立缸體有限元模型時進行了簡化處理。但在上橫梁與法蘭及缸筒接觸的位置處,上梁的局部擠壓變形會對法蘭和缸筒相應位置的高應力區(qū)的受力狀態(tài)有較大影響,因此在建立有限元模型時需考慮上橫梁下面板局部變形的影響。上橫梁材質為ZG270-500,屈服極限取為270 MPa,其下面板局部實體模型的形狀和尺寸見圖4b,設置它與缸體的接觸計算條件,摩擦因數取為0.12。為簡單起見,假設缸體和上橫梁材料均為理想彈塑性材料,服從Mises屈服準則,有限元模型中缸體屈服極限取為235 MPa。

    有限元模型采用了平面四邊形四節(jié)點單元。對法蘭圓角、缸底圓角附近區(qū)域進行了局部網格細化處理。細化區(qū)形狀尺寸如圖4b和圖4c中虛線所示?;诰W格敏感性分析結果,細化和非細化區(qū)域的網格邊長分別取為1 mm和8 mm。缸體共163 785個單元,上橫梁共劃分了35 165個單元。

    在缸體內壁密封以上部分施加均勻分布的液體壓力,密封段總高240 mm,在密封高度范圍內按照線性分布規(guī)律從上到下施加遞減的液體壓力。在橫梁局部實體模型上施加相應的位移約束。缸體軸對稱模型的受力簡圖見圖5,缸體的軸對稱有限元分析模型如圖6所示。

    圖5 缸體的加載力學示意圖

    圖6 缸體的有限元模型

    應用本文提出的安定性分析新方法對該缸體進行分析計算。施加圖1所示的峰值逐漸增大的壓力載荷,第一個壓力載荷的峰值為35 MPa,之后壓力載荷的峰值逐漸增大,每次增幅為1 MPa。

    2.2計算結果與討論

    圖7是壓力為35 MPa時缸體的等效應力云圖。最高應力出現在法蘭圓角附近區(qū)域,等效應力高達230 MPa左右。其他區(qū)域的等效應力均小于150 MPa。

    圖7 壓力為35 MPa時的缸體等效應力云圖

    法蘭圓角區(qū)域為缸體的最高應力區(qū),所以為缸體安定性判定的關鍵位置。由載荷峰值逐漸增大的加載-卸載計算得到,缸體法蘭圓角區(qū)域的關鍵點的最大等效殘余應力與壓力載荷的關系如圖8所示。根據上文的安定載荷判定條件可知,該液壓缸缸體的安定極限壓力為68 MPa。

    圖8 關鍵點的等效殘余應力與壓力的關系曲線

    壓力載荷為68 MPa時缸體的等效應力分布如圖9所示。缸體法蘭圓角區(qū)域的局部應力達到屈服極限,產生局部塑性變形。為了進一步考察液壓缸多次循環(huán)承載后塑性變形的演化情況,按照從零到68 MPa再到零的載荷循環(huán)方式,對缸體進行了50次循環(huán)彈塑性增量有限元分析計算。每次加載結束時法蘭圓角區(qū)域的等效塑性應變率云圖見圖10。由于等效塑性應變率不為零的區(qū)域為塑性變形區(qū),所以由圖10可見,在初始循環(huán)階段法蘭圓角區(qū)域存在局部塑性變形區(qū),隨著循環(huán)次數的增加,塑性變形區(qū)逐漸減小,最后變成完全彈性響應區(qū),達到安定狀態(tài)。其中塑性應變最大點的總塑性應變隨載荷循環(huán)數的變化曲線如圖11所示。由圖11可見,經過初始幾個循環(huán)后塑性變形穩(wěn)定不變,不再有新的塑性變形產生,達到安定狀態(tài),穩(wěn)定后的總塑性應變?yōu)?.22%。

    圖9 壓力載荷為68 MPa時缸體的等效應力云圖

    由此可見,由于以安定極限載荷作為缸體的彈塑性強度界限,利用了缸體材料的塑性承載能力,所以缸體即便是在68 MPa的壓力下也能夠安全服役。

    由上述安定性分析結果可知,該液壓缸的安定極限壓力與額定工作壓力的比值約為1.94,說明該缸體的彈塑性強度安全裕度較為寬裕。

    3結論

    (1)提出了一種安定性數值分析新方法。該方法通過構造合適的自平衡應力場,根據Melan安定下限定理,給出了安定性判定條件,并結合彈塑性有限元技術實現了結構的安定性分析。

    (a)第1個循環(huán)

    (b)第5個循環(huán)

    (c)第8個循環(huán)

    (d)第18個循環(huán)圖10 不同載荷循環(huán)次數下法蘭圓角區(qū)域的等效塑性應變率云圖

    圖11 關鍵點的總等效塑性應變隨循環(huán)數的變化曲線

    (2)125 MN自由鍛液壓機主液壓缸缸體的安定性分析結果顯示,缸體的彈塑性安全裕度為1.94,可以在高于其額定工作壓力90%的情況下正常安全服役。

    (3)本文提出的數值方法可以簡便有效地實現大型復雜承載結構的安定性分析,可以更合理地解決存在局部高應力區(qū)的大型復雜承載結構的強度分析問題。

    參考文獻:

    [1]ASME. ASME-Ⅷ-2 Boiler and Pressure Vessel Code[S]. New York: the American Society of Mechanical Engineers, 2007.

    [2]European Commission. EN-13445 Unfired Pressure Vessels[S]. Brussels: European Committee for Standardization, 2002.

    [3]Abdalla H F. Shakedown Boundary Determination of a 90 Degrees Back-to-back Pipe Bend Subjected to Steady Internal Pressures and Cyclic In-plane Bending Moments[J]. International Journal of Pressure Vessels and Piping, 2014, 116: 1-9.

    [4]Simon J W, Chen Geng, Weichert D. Shakedown Analysis of Nozzles in the Knuckle Region of Torispherical Heads under Multiple Thermo-mechanical Loadings[J]. International Journal of Pressure Vessels and Piping, 2014, 116: 47-55.

    [5]Mishra A, Kumar R S, Chellapandi P. Progressive Deformation Behaviour of Thin Cylindrical Shell under Cyclic Temperature Variation Using Combined Hardening Chaboche Model[J]. Latin American Journal of Solids and Structures, 2014, 11(6): 980-992.

    [6]Abdalla H F. Elastic Shakedown Boundary Determination of a Cylindrical Vessel-nozzle Intersection Subjected to Steady Internal Pressures and Cyclic Out-of-plane Bending Moments[J]. Nuclear Engineering and Design, 2014, 267: 189-196.

    [7]Vermaak N, Valdevit L, Anthony E, et al. Implementations of Shakedown for Design of Actively Cooled Thermo Structural Panels[J]. Journal of Mechanics of Materials and Structures, 2011, 6: 9-10.

    [8]Brouzoulis J, Torstensson P T, Stock R, et al. Prediction of Wear and Plastic Flow in Rails—Test Rig Results, Model Calibration and Numerical Prediction[J]. Wear, 2011, 271: 92-99.

    [9]Tafreshi S N M, Khalaj O, Dawson A R. Repeated Loading of Soil Containing Granulated Rubber and Multiple Geocell Layers[J]. Geotextiles and Geomembranes, 2014, 42(1): 25-38.

    [10]Fadaee M J, Saffari H, Tabatabaei R. Shakedown Limit of Elastic-plastic Offshore Structures under Cyclic Wave Loading[J]. Ocean Engineering, 2008, 35(17/18): 1854-1861.

    [11]Atkociunas J, Venskus A. Optimal Shakedown Design of Frames under Stability Conditions According to Standards[J]. Computers & Structures, 2011, 89(3/4): 435-443.

    [13]Chen Shenshen, Liu Yinghua, Li Jun, et al. Performance of the MLPG Method for Static Shakedown Analysis for Bounded Kinematic Hardening Structures[J]. European Journal of Mechanics A/Solids, 2011, 30(2): 183-194.

    [14]Ure J, Chen Haofeng, Tipping D. Integrated Structural Analysis Tool Using the Linear Matching Method Part 1-Software Development[J]. International Journal of Pressure Vessels and Piping, 2014, 120:141-151.

    [15]Simon J W, Hoewer D, Weichert D. A Starting-point Strategy for Interior-point Algorithms for Shakedown Analysis of Engineering Structures[J]. Engineering Optimization, 2014, 46(5): 648-668.

    (編輯蘇衛(wèi)國)

    A New Shakedown Numerical Analysis Method for Huge Complex Strcuture and Its Applications

    Zou Zongyuan1Guo Baofeng1Jin Miao1Mao Chunyan2

    1.Key Laboratory of Advanced Forging & Stamping Technology and Science (Yanshan University),Ministry of Education of China,Qinhuangdao,Hebei,066004 2.Taiyuan Heavy Industry Co.,Ltd.,Taiyuan,030024

    Abstract:The local high stresses appeared commonly in a loaded huge complex structure. Some of them were usually close to or even exceed the material yield strength due to the restrictions of the geometric size and the practical manufacture. It was troublesome with the applications of the conventional elastic strength design rule. The elastic-plastic shakedown design rule might utilize the plastic bearing capacity of material, then to solve this kind of strength problem well. Thus the shakedown theorem was introduced, and based on Melan shakedown theorem, a self-equilibrium stress field was constructed by subtractting the elastic-plastic stress field and the elastic reference stress field of the peak load, and the corresponding shakedown condition was derived. Then with the applications of the elastic-plastic incremental finite element technique, the loading-unloading calculations for different load levels were conducted to obtain the residual stress fields. Taking consideration of the relation between the residual stress field and the constructed self-equilibrium stress field, the shakedown limit load was determined by the shakedown conditions. Then a shakedown numerical analysis method for huge complex structure was established. The shakedown analysis for the main hydraulic cylinder of a 125 MN forging hydraulic press was performed by this new method.

    Key words:elastic-plastic;shakedown;hydraulic cylinder body;shakedown limit load;local high stress

    收稿日期:2015-08-17

    基金項目:國家科技重大專項(2013ZX04003031);國家自然科學基金資助項目(51475408);河北省高校創(chuàng)新團隊領軍人才培育計劃資助項目 (LJRC012);河北省自然科學基金資助項目(E2015203223)

    中圖分類號:TG315.4

    DOI:10.3969/j.issn.1004-132X.2016.03.013

    作者簡介:鄒宗園,女,1986年生。燕山大學機械工程學院博士研究生。郭寶峰,男,1958年生。燕山大學機械工程學院教授、博士研究生導師。金淼(通信作者),男,1968年生。燕山大學機械工程學院教授、博士研究生導師。毛春燕,女,1978年生。太原重工股份有限公司技術中心工程師。

    猜你喜歡
    安定性彈塑性液壓缸
    冷軋雙切剪液壓缸螺栓斷裂分析
    冶金設備(2020年2期)2020-12-28 00:15:26
    矮塔斜拉橋彈塑性地震響應分析
    一種基于液壓缸負負載回油腔的調速控制系統(tǒng)
    重型機械(2020年2期)2020-07-24 08:16:20
    基于Android系統(tǒng)的汽油氧化安定性測定儀的設計
    電子制作(2018年11期)2018-08-04 03:25:50
    烴族組成對柴油儲存安定性影響研究
    彈塑性分析在超高層結構設計中的應用研究
    江西建材(2018年4期)2018-04-10 12:36:52
    液壓缸維修技術
    動載荷作用下冪硬化彈塑性彎曲裂紋塑性區(qū)
    關于水泥檢驗中安定性的時效問題
    一種新型旋轉直線組合式液壓缸設計
    国产黄片视频在线免费观看| 亚洲精品成人久久久久久| 亚洲中文字幕日韩| 波多野结衣高清无吗| 亚洲国产精品sss在线观看| 深夜a级毛片| 搡女人真爽免费视频火全软件| 插逼视频在线观看| 精华霜和精华液先用哪个| 身体一侧抽搐| 婷婷色综合大香蕉| 国产一区二区在线av高清观看| 搡老妇女老女人老熟妇| 最后的刺客免费高清国语| 色视频www国产| 99久久成人亚洲精品观看| 国产淫语在线视频| 国产精品一区二区三区四区免费观看| 久久久久精品久久久久真实原创| 国产亚洲最大av| 日本av手机在线免费观看| 97热精品久久久久久| 久久久久性生活片| 一级毛片aaaaaa免费看小| 亚洲成人av在线免费| 身体一侧抽搐| av免费在线看不卡| 久久久久久国产a免费观看| 国产精品嫩草影院av在线观看| 韩国高清视频一区二区三区| 内地一区二区视频在线| 国产av一区在线观看免费| av在线亚洲专区| 你懂的网址亚洲精品在线观看 | 嫩草影院新地址| 99久国产av精品国产电影| 少妇人妻精品综合一区二区| 久久久欧美国产精品| 欧美不卡视频在线免费观看| 婷婷色综合大香蕉| 免费看美女性在线毛片视频| 我要搜黄色片| 少妇人妻精品综合一区二区| 久久久久久久亚洲中文字幕| 久久亚洲国产成人精品v| 人人妻人人澡人人爽人人夜夜 | 欧美最新免费一区二区三区| 国产极品天堂在线| 美女xxoo啪啪120秒动态图| 插阴视频在线观看视频| 午夜福利成人在线免费观看| 精品国产一区二区三区久久久樱花 | 寂寞人妻少妇视频99o| 中文欧美无线码| 免费不卡的大黄色大毛片视频在线观看 | 欧美bdsm另类| АⅤ资源中文在线天堂| 少妇猛男粗大的猛烈进出视频 | 蜜臀久久99精品久久宅男| 国产精品av视频在线免费观看| 国产精品日韩av在线免费观看| av女优亚洲男人天堂| 免费黄网站久久成人精品| 亚洲人成网站在线观看播放| 日本-黄色视频高清免费观看| 国产精品,欧美在线| 国产黄片视频在线免费观看| 国产精品一二三区在线看| 黄色欧美视频在线观看| 日韩人妻高清精品专区| 少妇丰满av| 免费不卡的大黄色大毛片视频在线观看 | 纵有疾风起免费观看全集完整版 | 国产爱豆传媒在线观看| 爱豆传媒免费全集在线观看| 天天躁日日操中文字幕| 嘟嘟电影网在线观看| 夜夜爽夜夜爽视频| 国产av码专区亚洲av| 日韩一本色道免费dvd| 亚洲熟妇中文字幕五十中出| 嫩草影院入口| 日韩一区二区三区影片| 国产一区有黄有色的免费视频 | 男女那种视频在线观看| 中国美白少妇内射xxxbb| 国产午夜福利久久久久久| 五月玫瑰六月丁香| 成人无遮挡网站| videos熟女内射| 日本熟妇午夜| 欧美精品一区二区大全| 日本色播在线视频| 亚洲av男天堂| 欧美激情国产日韩精品一区| 国产又色又爽无遮挡免| 亚洲av中文字字幕乱码综合| 国产伦理片在线播放av一区| 国产日韩欧美在线精品| 久久这里有精品视频免费| 国产 一区 欧美 日韩| 晚上一个人看的免费电影| 亚洲国产精品成人久久小说| 色播亚洲综合网| 91久久精品国产一区二区三区| 99在线人妻在线中文字幕| www.色视频.com| 日本色播在线视频| 亚洲自拍偷在线| 国产午夜精品久久久久久一区二区三区| 亚洲成色77777| 国产成人一区二区在线| 麻豆成人av视频| 国产精品不卡视频一区二区| 国产精品不卡视频一区二区| 亚洲在线自拍视频| 久久草成人影院| 日本黄色片子视频| 国产一级毛片在线| 国产成人aa在线观看| 国产精品女同一区二区软件| 69人妻影院| 午夜亚洲福利在线播放| 中文精品一卡2卡3卡4更新| www日本黄色视频网| 欧美三级亚洲精品| 三级经典国产精品| www.av在线官网国产| 久久久a久久爽久久v久久| 欧美成人精品欧美一级黄| 极品教师在线视频| 中文字幕亚洲精品专区| 久久精品国产亚洲av天美| 久久人妻av系列| 日日摸夜夜添夜夜添av毛片| 国产成人a区在线观看| 中文字幕av在线有码专区| 白带黄色成豆腐渣| 精品熟女少妇av免费看| 三级国产精品欧美在线观看| 22中文网久久字幕| 亚洲乱码一区二区免费版| 国产黄片美女视频| 国产伦精品一区二区三区四那| 中文字幕久久专区| 国产又黄又爽又无遮挡在线| 国产中年淑女户外野战色| 国产黄色视频一区二区在线观看 | 久久久精品大字幕| 联通29元200g的流量卡| 午夜视频国产福利| 日韩成人伦理影院| 久久久欧美国产精品| 黄色一级大片看看| 欧美区成人在线视频| 亚洲色图av天堂| 在线天堂最新版资源| 天堂中文最新版在线下载 | 久久这里只有精品中国| 97超视频在线观看视频| 嫩草影院精品99| 69av精品久久久久久| 天天躁日日操中文字幕| 一级毛片久久久久久久久女| 特级一级黄色大片| 国产精品美女特级片免费视频播放器| 床上黄色一级片| 美女国产视频在线观看| 99久久精品热视频| 亚州av有码| 色综合站精品国产| 我要搜黄色片| 天天躁日日操中文字幕| 久久久久九九精品影院| 免费电影在线观看免费观看| 全区人妻精品视频| 欧美变态另类bdsm刘玥| 99久久无色码亚洲精品果冻| 国产精品日韩av在线免费观看| 日韩精品青青久久久久久| 欧美+日韩+精品| 丝袜喷水一区| 国产精品一区二区三区四区久久| 国产精品综合久久久久久久免费| 免费无遮挡裸体视频| 国产免费男女视频| 黄色欧美视频在线观看| 男人狂女人下面高潮的视频| 精品无人区乱码1区二区| 欧美激情在线99| 97热精品久久久久久| 99久国产av精品国产电影| 欧美zozozo另类| 亚洲国产精品成人综合色| 女人十人毛片免费观看3o分钟| 精品免费久久久久久久清纯| 舔av片在线| 视频中文字幕在线观看| 69av精品久久久久久| 日本一二三区视频观看| 99久国产av精品国产电影| 国产美女午夜福利| 日韩一区二区三区影片| 久99久视频精品免费| 精品人妻熟女av久视频| 欧美xxxx性猛交bbbb| 久久午夜福利片| 男女边吃奶边做爰视频| 激情 狠狠 欧美| 日韩三级伦理在线观看| 日韩三级伦理在线观看| 国产一区二区亚洲精品在线观看| av免费在线看不卡| 国产不卡一卡二| 亚洲精品成人久久久久久| www日本黄色视频网| 日韩人妻高清精品专区| 国产免费福利视频在线观看| 欧美日韩综合久久久久久| 三级经典国产精品| 久久精品91蜜桃| 欧美+日韩+精品| 亚洲av中文字字幕乱码综合| 亚洲av成人精品一区久久| videos熟女内射| 又爽又黄无遮挡网站| 中文字幕精品亚洲无线码一区| 久久久精品94久久精品| 日韩欧美三级三区| 国内揄拍国产精品人妻在线| 久久久亚洲精品成人影院| 亚洲成人av在线免费| 麻豆成人av视频| 99热这里只有是精品50| 中文字幕久久专区| 韩国高清视频一区二区三区| 麻豆成人午夜福利视频| 精品一区二区三区视频在线| 精品一区二区三区视频在线| 啦啦啦啦在线视频资源| 2021少妇久久久久久久久久久| 18禁在线播放成人免费| 日韩欧美三级三区| 亚洲一级一片aⅴ在线观看| 全区人妻精品视频| 亚洲国产高清在线一区二区三| 久久精品人妻少妇| 在线观看av片永久免费下载| 中文字幕av成人在线电影| 精品熟女少妇av免费看| 黄片无遮挡物在线观看| 九九久久精品国产亚洲av麻豆| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲真实伦在线观看| 2022亚洲国产成人精品| 久久欧美精品欧美久久欧美| 日韩欧美 国产精品| 欧美又色又爽又黄视频| 亚洲精品亚洲一区二区| 国产一级毛片七仙女欲春2| 日本免费在线观看一区| 国产精品一区二区性色av| 插逼视频在线观看| 亚洲av日韩在线播放| 国产精品一区二区在线观看99 | 搡女人真爽免费视频火全软件| 国产成年人精品一区二区| 嫩草影院精品99| 亚洲高清免费不卡视频| 水蜜桃什么品种好| 久久精品夜色国产| 日本一本二区三区精品| 小说图片视频综合网站| 国国产精品蜜臀av免费| 成人亚洲精品av一区二区| 亚洲成色77777| 日本wwww免费看| 好男人视频免费观看在线| 欧美97在线视频| 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 日韩中字成人| 水蜜桃什么品种好| 人人妻人人澡人人爽人人夜夜 | 久久国产乱子免费精品| 最近视频中文字幕2019在线8| 国产色爽女视频免费观看| 日韩av在线免费看完整版不卡| 少妇裸体淫交视频免费看高清| 亚洲av电影不卡..在线观看| 只有这里有精品99| 日韩av在线大香蕉| 可以在线观看毛片的网站| 亚洲精品亚洲一区二区| av国产久精品久网站免费入址| 久久这里有精品视频免费| 哪个播放器可以免费观看大片| 久久久久九九精品影院| 91狼人影院| 3wmmmm亚洲av在线观看| 成年免费大片在线观看| 久久精品91蜜桃| 听说在线观看完整版免费高清| 一卡2卡三卡四卡精品乱码亚洲| 婷婷色麻豆天堂久久 | 熟女电影av网| 最后的刺客免费高清国语| 国内精品一区二区在线观看| 爱豆传媒免费全集在线观看| 观看免费一级毛片| 国产av在哪里看| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 日本与韩国留学比较| 国产精品久久视频播放| 久久精品国产鲁丝片午夜精品| 亚洲自拍偷在线| 禁无遮挡网站| 亚洲,欧美,日韩| 综合色av麻豆| 久久久a久久爽久久v久久| 国产av在哪里看| 少妇猛男粗大的猛烈进出视频 | 狠狠狠狠99中文字幕| 女人被狂操c到高潮| 国产免费福利视频在线观看| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱| 国产精品一区二区三区四区久久| 国产亚洲av片在线观看秒播厂 | 国产精品国产三级国产av玫瑰| 色尼玛亚洲综合影院| 黄片wwwwww| 日本黄色片子视频| 深夜a级毛片| 久久99热6这里只有精品| 男的添女的下面高潮视频| 国产精品熟女久久久久浪| 亚洲成人av在线免费| 国产精品久久久久久av不卡| 国产免费又黄又爽又色| 午夜a级毛片| 国产成人aa在线观看| 晚上一个人看的免费电影| 国产亚洲5aaaaa淫片| 一区二区三区免费毛片| videos熟女内射| 一级黄色大片毛片| 亚洲成人中文字幕在线播放| 床上黄色一级片| 亚洲国产最新在线播放| 中文字幕制服av| 亚洲精华国产精华液的使用体验| 欧美三级亚洲精品| 国产淫语在线视频| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩精品一区二区| 亚洲成色77777| 久久精品久久久久久噜噜老黄 | 97热精品久久久久久| 禁无遮挡网站| 内地一区二区视频在线| 久久久精品大字幕| 精品国产一区二区三区久久久樱花 | 日日干狠狠操夜夜爽| 日韩视频在线欧美| 国产伦精品一区二区三区视频9| 精品无人区乱码1区二区| 欧美最新免费一区二区三区| 综合色av麻豆| 可以在线观看毛片的网站| 久久这里只有精品中国| 国产美女午夜福利| 久久精品人妻少妇| 男人和女人高潮做爰伦理| 国产精品国产高清国产av| 国产成人精品婷婷| 欧美97在线视频| 成人国产麻豆网| 久久99热6这里只有精品| 国产免费一级a男人的天堂| 一个人观看的视频www高清免费观看| 精品久久久久久久人妻蜜臀av| 精品国产三级普通话版| 色尼玛亚洲综合影院| 亚洲乱码一区二区免费版| 亚洲四区av| 1000部很黄的大片| 狠狠狠狠99中文字幕| 色尼玛亚洲综合影院| 国产精品.久久久| 亚洲人成网站高清观看| av在线亚洲专区| 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 乱码一卡2卡4卡精品| 99久久精品一区二区三区| 精品国产一区二区三区久久久樱花 | 久久99热6这里只有精品| 免费av观看视频| 只有这里有精品99| 高清av免费在线| 国语对白做爰xxxⅹ性视频网站| 国产精品一区www在线观看| 国产在视频线在精品| 久久精品91蜜桃| 中文精品一卡2卡3卡4更新| 国产精品熟女久久久久浪| 一级爰片在线观看| 国产极品精品免费视频能看的| 亚洲av成人av| 99在线视频只有这里精品首页| 91狼人影院| 久久久久精品久久久久真实原创| 插阴视频在线观看视频| 国产精品野战在线观看| 国产午夜精品久久久久久一区二区三区| 视频中文字幕在线观看| 久久精品国产鲁丝片午夜精品| 纵有疾风起免费观看全集完整版 | 人人妻人人看人人澡| av免费观看日本| 亚洲国产最新在线播放| 国产精品电影一区二区三区| 国产精品无大码| 18禁裸乳无遮挡免费网站照片| 国产欧美日韩精品一区二区| 97人妻精品一区二区三区麻豆| 国产高清国产精品国产三级 | 国产精品国产三级国产av玫瑰| 亚洲av免费高清在线观看| 男女那种视频在线观看| 精品熟女少妇av免费看| 色综合站精品国产| 99久国产av精品| 国产探花极品一区二区| 午夜日本视频在线| 我要看日韩黄色一级片| 热99re8久久精品国产| 亚洲,欧美,日韩| 日韩大片免费观看网站 | 日本一二三区视频观看| 99久久精品一区二区三区| 自拍偷自拍亚洲精品老妇| 国产高清有码在线观看视频| 尤物成人国产欧美一区二区三区| 日本黄大片高清| 亚洲国产欧洲综合997久久,| 观看美女的网站| 大香蕉97超碰在线| 国模一区二区三区四区视频| 91在线精品国自产拍蜜月| 亚洲av熟女| 99热6这里只有精品| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 99热这里只有精品一区| 国产高潮美女av| 亚洲欧美日韩卡通动漫| 亚洲三级黄色毛片| av国产免费在线观看| 亚洲最大成人手机在线| 国产伦一二天堂av在线观看| 美女国产视频在线观看| 久久久久性生活片| av福利片在线观看| 九九热线精品视视频播放| 色尼玛亚洲综合影院| 亚洲精品日韩在线中文字幕| videossex国产| 亚洲成av人片在线播放无| 国产精品久久久久久精品电影| 99久国产av精品国产电影| 久久精品国产亚洲av涩爱| 婷婷色综合大香蕉| 亚洲国产色片| 伊人久久精品亚洲午夜| 国产国拍精品亚洲av在线观看| 嫩草影院入口| 黄色欧美视频在线观看| 日本免费在线观看一区| 校园人妻丝袜中文字幕| 国产精品永久免费网站| 日本午夜av视频| 精品无人区乱码1区二区| 欧美最新免费一区二区三区| 性插视频无遮挡在线免费观看| 看十八女毛片水多多多| 99国产精品一区二区蜜桃av| 看免费成人av毛片| 亚洲自偷自拍三级| 国产在线一区二区三区精 | 观看美女的网站| 青春草视频在线免费观看| 欧美精品国产亚洲| 日本一本二区三区精品| 我的女老师完整版在线观看| 乱码一卡2卡4卡精品| 99久久精品热视频| 嘟嘟电影网在线观看| 国产亚洲5aaaaa淫片| 国产精品国产三级国产专区5o | 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 午夜免费激情av| 欧美最新免费一区二区三区| 91精品伊人久久大香线蕉| 白带黄色成豆腐渣| 日本爱情动作片www.在线观看| 搞女人的毛片| 欧美一区二区亚洲| 中文资源天堂在线| 青青草视频在线视频观看| 一个人看视频在线观看www免费| 日本免费a在线| 亚洲,欧美,日韩| 成人午夜高清在线视频| h日本视频在线播放| or卡值多少钱| 亚洲在久久综合| 免费观看的影片在线观看| 国产乱人视频| 最近2019中文字幕mv第一页| 午夜福利高清视频| 麻豆精品久久久久久蜜桃| 精品人妻视频免费看| 免费看日本二区| 亚洲国产精品合色在线| av免费在线看不卡| 国产精品麻豆人妻色哟哟久久 | 国产精品国产高清国产av| 91在线精品国自产拍蜜月| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站| 久久久久久久久中文| 国产精品爽爽va在线观看网站| 亚洲,欧美,日韩| 丝袜美腿在线中文| 国产高清不卡午夜福利| 一本久久精品| 亚洲成av人片在线播放无| av在线天堂中文字幕| 久久精品国产亚洲av涩爱| 少妇猛男粗大的猛烈进出视频 | 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 亚洲不卡免费看| a级毛片免费高清观看在线播放| 黄色日韩在线| 最新中文字幕久久久久| 成人综合一区亚洲| 女人十人毛片免费观看3o分钟| 午夜福利网站1000一区二区三区| 久久久久久久国产电影| 久久精品综合一区二区三区| 小说图片视频综合网站| 美女高潮的动态| 免费电影在线观看免费观看| 欧美激情国产日韩精品一区| 欧美色视频一区免费| 97人妻精品一区二区三区麻豆| 免费观看精品视频网站| 欧美激情在线99| 永久网站在线| 久久久久久久国产电影| 91久久精品国产一区二区三区| 成人亚洲欧美一区二区av| 桃色一区二区三区在线观看| 亚洲怡红院男人天堂| 成人美女网站在线观看视频| a级毛片免费高清观看在线播放| 国产精品乱码一区二三区的特点| 成人无遮挡网站| 久久精品久久久久久久性| 午夜福利在线观看免费完整高清在| 日韩高清综合在线| 美女脱内裤让男人舔精品视频| 日本欧美国产在线视频| 我的老师免费观看完整版| 日韩亚洲欧美综合| 国产极品天堂在线| 婷婷色av中文字幕| 精品久久久久久电影网 | 黄色日韩在线| 2021少妇久久久久久久久久久| 两个人的视频大全免费| 乱系列少妇在线播放| 婷婷色麻豆天堂久久 | 亚洲自拍偷在线| 别揉我奶头 嗯啊视频| 最近最新中文字幕免费大全7| 国产伦在线观看视频一区| 国产淫片久久久久久久久| 午夜a级毛片| 中文资源天堂在线| 久久久久网色| 久久亚洲精品不卡| 国产探花在线观看一区二区| 性插视频无遮挡在线免费观看| av在线观看视频网站免费| 成人综合一区亚洲| 免费搜索国产男女视频| 自拍偷自拍亚洲精品老妇| 大又大粗又爽又黄少妇毛片口| 六月丁香七月| 国产一区二区三区av在线| 久久这里有精品视频免费| 久久久久九九精品影院| 97人妻精品一区二区三区麻豆| 国产在线男女| 2021少妇久久久久久久久久久| 黄片无遮挡物在线观看| 99热这里只有是精品50| 男女视频在线观看网站免费| 国产av在哪里看|