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

    提高軸孔連接結構承載能力的設計方法

    2017-08-27 05:36:11劉文章何景武嚴斌吳文偉
    船舶力學 2017年8期
    關鍵詞:接觸區(qū)軸孔耳片

    劉文章,何景武,嚴斌,吳文偉

    (1.中國船舶科學研究中心,江蘇無錫214082;2.北京航空航天大學航空科學與工程學院,北京100191)

    提高軸孔連接結構承載能力的設計方法

    劉文章1,何景武2,嚴斌1,吳文偉1

    (1.中國船舶科學研究中心,江蘇無錫214082;2.北京航空航天大學航空科學與工程學院,北京100191)

    文章針對軸孔連接結構提出了一種改善軸孔接觸特性來提高承載能力的方法。以工程中通用的螺栓連接結構為例,建立接觸分析模型,應用Persson接觸理論和優(yōu)化分析方法,通過優(yōu)化耳片孔接觸面的曲面形狀和曲面參數(shù)實現(xiàn)了降低軸孔接觸區(qū)最大接觸應力及改善其分布梯度的目的。針對優(yōu)化后的耳片孔接觸曲面設計結果,利用有限元素法建立了螺栓耳片連接結構模型,分析驗算了接觸應力。結果表明,優(yōu)化分析的解析法和有限元素法的計算結果基本一致。連接結構優(yōu)化設計后,接觸區(qū)最大應力降低50%以上,應力分布梯度明顯改善,提高了其承載能力。

    軸孔連接結構;Persson接觸理論;接觸應力;結構優(yōu)化;承載能力

    0 引言

    軸孔連接在船舶等工程結構設計中應用非常普遍,是機械結構中重要的傳力部件,比如螺栓與耳片。由于其本身的結構特點,承力時往往是外力通過軸直接施加在孔邊上,孔邊會因為擠壓而出現(xiàn)接觸應力分布集中的現(xiàn)象,導致其承載能力不高,而且在動態(tài)載荷作用下,容易造成局部點的接觸受損。因此孔邊接觸區(qū)域往往會成為影響整個連接結構承載能力的關鍵部位。

    但是,到目前為止,國內(nèi)外大部分的研究主要集中在接觸應力的求解方法上,對改善接觸應力分布的措施與方法研究不多。參考文獻[1]提出了一種利用能量守恒來計算接觸應力的方法。參考文獻[2-3]主要針對連接件接觸應力分布的求解領域。參考文獻[4-5]利用有限元方法來模擬計算接觸應力。參考文獻[6-7]提出了一種通過改變物體形狀來降低接觸應力的構想,但是沒有形成系統(tǒng)條理的改進方法。

    本文建立了一種系統(tǒng)地改善軸孔連接結構接觸區(qū)應力分布的方法。建立連接結構的平面接觸模型,確定接觸應力的理論求解方法,通過改變耳片孔接觸面的曲面形狀并優(yōu)化曲面參數(shù),使耳片孔邊接觸區(qū)的最大應力及應力分布梯度得到降低,如圖1所示,接觸區(qū)應力分布由初始圖線1變?yōu)槟繕藞D線2。計算軸孔連接結構接觸應力的有限元解,與理論解誤差對比,驗證優(yōu)化設計的準確性。通過上述過程來探究耳片孔接觸面幾何構型對接觸區(qū)應力分布的影響,達到提高承載能力的目的。

    圖1 優(yōu)化設計目標示意圖Fig.1 Schematic diagram of optimal target

    圖2 軸孔連接結構的接觸示意圖Fig.2 Schematic diagram of contact progress

    1 平面接觸模型的建立

    1.1 接觸過程分析

    通常在使用中,按照配合精度要求,螺栓軸和耳片孔是間隙配合,但是配合間隙很小,配合表面在無變形時已經(jīng)相當接近地貼合在一起。外力載荷P一般通過螺栓軸直接擠壓耳片而傳遞載荷,隨著外部載荷的逐步增加,軸孔接觸由線接觸(圖2虛線所示)迅速擴展到面接觸,使其與軸孔本身的有效幾何尺寸相當。這一邊界非線性的協(xié)調(diào)接觸過程,不適合彈性半空間的假設,不滿足Hertz接觸理論的適用條件[8-9],需要用Persson接觸理論求解。

    1.2 平面接觸模型

    圖3 平面接觸模型Fig.3 Model of plane contact

    圖4 實例幾何模型示意圖Fig.4 Geometric model

    在載荷作用下的接觸過程,沿軸孔連接結構的軸向基本不變,連接結構軸向尺寸相對較大且受力沿軸向變化較小,因此假定軸與孔的應力狀態(tài)為平面應變狀態(tài),建立如圖3的平面接觸模型,A為螺栓軸截面,B為耳片截面。以某結構物設計中所應用的螺栓耳片結構為例,如圖4所示。

    假定螺栓軸與耳片孔的配合精度為H8/f7[10],平面接觸模型中的幾何參數(shù)為:孔基本尺寸17,孔上偏差ES=0.027,軸下偏差ei=0.034,最大間隙△R=ES-ei=0.061 0,計算中設定螺栓軸半徑rA=16.939,耳片孔半徑rB=17.000。材料參數(shù)為:彈性模量EA=210 GPa,EB=195 GPa,泊松比νA=νB=0.3。外部載荷為:P=285.3 kN。

    2 接觸應力的理論求解

    應用改進的Persson接觸理論[11-12],求解軸孔連接結構的接觸應力分布規(guī)律。改進的Persson接觸理論,有以下兩個假設:(1)忽略接觸區(qū)徑向位移和切向位移間的耦合作用;(2)接觸區(qū)法向最大變形與徑向位移之間的關系滿足剛體接觸的幾何剛度條件。忽略接觸區(qū)域的摩擦作用[4],軸孔連接結構接觸應力的計算,符合Persson理論的兩個假設。

    軸A與孔B在點1處配合接觸,在A的軸心Oe處作用一外部集中力P,接觸區(qū)圓弧對應弧度為2ε,接觸區(qū)圓弧關于外力作用線對稱分布,如圖5所示。

    載荷P與接觸半弧度ε之間的關系[11-12]如下:

    圖5 Persson理論模型示意圖Fig.5 Schematic diagram of Persson model

    接觸壓力分布p(y)與變量角φ關系[11-12]如下:

    圖6 初始模型的接觸應力分布圖Fig.6 Contact stress distribution of initial model

    (1)式是關于ε的非線性方程,在外部載荷P已知的情況下,通過數(shù)值分析方法求得。代入(2)式即可計算得到接觸應力理論解。未優(yōu)化的初始模型接觸區(qū)應力p(y)隨角度φ的變化規(guī)律,如圖6所示。

    接觸弧度參數(shù)b=tan( ε/2)=0.163 0,接觸應力最大值pmax=33 292.0 kPa。

    3 耳片孔的優(yōu)化設計

    3.1 優(yōu)化設計模型

    接觸應力大小與兩個因素有關:(1)接觸點的曲率;(2)接觸點的數(shù)目。

    (1)接觸點曲率:如表1所示,rA不變,rB變小,接觸區(qū)域弧度變大,最大接觸應力變小,符合優(yōu)化方向。

    (2)接觸點數(shù)目:接觸點越多,承載區(qū)域會越大,接觸應力就會越分散,最大接觸應力也會越小,符合優(yōu)化方向。

    螺栓一般為標準件,選擇耳片孔為優(yōu)化對象。設計中選擇耳片孔邊接觸曲面為便于精加工的橢圓形狀,符合上述兩點優(yōu)化方向要求。建立優(yōu)化設計模型,通過優(yōu)化耳片孔的曲面參數(shù),使接觸區(qū)的接觸應力最大值pmax最小。

    優(yōu)化設計模型表述為:

    表1 接觸曲率變化對pmax影響表Tab.1 Effect of contact curvature on pmax

    如圖7所示,軸A圓的方程為:

    耳片孔B橢圓的方程為:

    式中:R為圓的半徑,c為圓心到X軸的距離,a為橢圓的短半軸,b1為橢圓的長半軸。

    軸A半徑尺寸R=16.939。根據(jù)配合精度,橢圓孔短軸a=17.00固定,選擇長軸b1為優(yōu)化設計變量,b1不同,對應接觸點的曲率就會變化。

    3.2 優(yōu)化設計方法

    圓與橢圓接觸,如圖7所示,最初接觸點為C和D。外部載荷P平均分配到兩個接觸點上。

    圖7 圓與橢圓接觸示意圖Fig.7 Schematic diagram of contact between circle and ellipse

    便于計算方便,設計變量由變量b1換算為接觸點和Y軸的角度θ。

    圓與橢圓在C點切線斜率相等,

    C點位于橢圓上,

    (7)式和方程(8)得到θ與b1關系式:

    接觸點處的橢圓曲率半徑為

    設計變量θ值變化,最初接觸點C和D都會變化,C和D處相應的曲率隨之變化,計算得到C和D處各自的接觸應力分布。

    θ值不同,C和D各自接觸區(qū)大小不同,分為C和D接觸區(qū)域是否重合兩種情況。

    (1)接觸區(qū)未重疊情況的接觸應力求解:

    初始接觸角θ大于等于各接觸點接觸弧度ε,即θ≥ε,如圖8所示。此時兩個接觸點的應力分布沒有重疊。

    兩個接觸區(qū)各自的應力分布分別為:

    θ=30°時的接觸應力分布圖,如圖9所示。

    圖8 未重疊接觸角度示意圖Fig.8 Schematic diagram of contact angle without overlap

    圖10 重疊接觸角度示意圖Fig.10 Schematic diagram of contact angle with overlap

    圖11 θ=22°時接觸應力分布圖Fig.11 θ=22°:Contact stress distribution

    忽略兩個接觸點接觸區(qū)的相互影響,計算兩個接觸區(qū)的應力分布。重疊部分,兩個區(qū)域接觸應力相加;未重疊部分,接觸應力保持不變。

    接觸區(qū)應力分布為:

    θ=22°時的接觸應力分布圖,如圖11。

    (3)優(yōu)化設計過程:

    隨著θ值變化,接觸區(qū)應力首先在兩側單獨存在,然后在中間部分重疊。最大應力由大變小,再變大,存在最小值。變化趨勢如圖12中箭頭方向所示。最大接觸應力最小情況發(fā)生在兩個接觸區(qū)域重疊的情況下。

    當x=0時的疊加應力p(0)與以C、D點為中心的單側最大應力相等時,最大接觸應力最小,滿足設計目標,即:p(0)=max p(y),計算可得最優(yōu)解θ值和最優(yōu)橢圓曲線參數(shù)b1。

    圖12 θ=22°-15°變化時接觸應力變化圖Fig.12 Change of contact stress when θ=22°-15°

    3.3 最優(yōu)結果

    如表2所示,最優(yōu)解θ=21.096 3°,最優(yōu)橢圓曲線參數(shù)b1=17.070 3,最大接觸應力值降低了58.7%,優(yōu)化效果非常顯著。

    表2 優(yōu)化設計結果Tab.2 The optimization design results

    4 有限元解和理論解誤差對比

    4.1 有限元建模[13-14]

    接觸問題是一種極度不連續(xù)的高度邊界非線性問題。兩個面接觸,必須施加相應的約束關系,并判斷接觸面積和接觸壓力。兩個面分開,則沒有接觸存在,就必須移走相應的約束。接觸問題的特點決定了其通常采用增量法求解,采用罰函數(shù)法將界面條件引入到求解方程。本文應用Abaqus軟件,采用Newton-Raphson法求解非線性方程。

    對初始未優(yōu)化模型(圓與圓接觸)和優(yōu)化后模型(圓與優(yōu)化設計得到最優(yōu)參數(shù)解的橢圓接觸)兩種情況進行建模。采用平面應變簡化,忽略接觸面摩擦[4],建立1/2對稱模型,對稱面上施加對稱邊界條件,即U1=UR2= UR3=0,其他周邊邊界為固支。

    統(tǒng)一采用CPE4R單元,結構性網(wǎng)格劃分。在接觸面附近對網(wǎng)格進行加密,以獲得良好的計算結果。選擇面面接觸,A為主面,B為從面,小滑移,無摩擦。加載方式為圓A圓心位置豎直外力P0=283.5 kN。圖13左圖未優(yōu)化模型中Circle A為圓A,單元類型為quad,單元數(shù)為526;Circle B為圓B,單元類型為quad,單元數(shù)為2 411。圖13右圖優(yōu)化模型中Circle A為圓A,單元類型為quad,單元數(shù)為848;Ellipse為橢圓B,單元類型為quad,單元數(shù)為3 308。

    圖13 有限元模型(左圖為未優(yōu)化,右圖為優(yōu)化)Fig.13 FEM models(left:not optimal,right:optimal)

    4.2 有限元計算結果

    計算分析得出兩種模型的接觸應力云圖,如圖14所示。計算得出兩種模型的接觸應力曲線圖,如圖15所示,p為接觸應力值,X為從圓最下端位置起沿圓弧向上的距離值。

    圖14 接觸應力云圖(左圖:未優(yōu)化,右圖:優(yōu)化)Fig.14 Contact stress contours(left:not optimal,right:optimal)

    圖15 接觸應力分布曲線圖Fig.15 Curve gragh of contact stress distribution

    圖14和圖15分析得出,圓與圓接觸區(qū)的最大應力發(fā)生在中間最下端位置,接觸圓弧較?。粓A與橢圓接觸區(qū)的最大應力也發(fā)生在中間最下端位置,接觸圓弧較大;圓與橢圓接觸情況下的接觸應力比圓與圓接觸情況下的接觸應力變化更加平緩,接觸區(qū)域面積更大,應力分布梯度更小,最大應力減小了54.8%。

    4.3 有限元解與理論解對比分析

    表3分別列出了初始模型和優(yōu)化模型的理論解析解和有限元解結果,有限元解和理論解基本吻合,數(shù)據(jù)誤差基本在5%左右,符合計算精度要求。兩方面都驗證了優(yōu)化分析與設計的準確性。

    表3 優(yōu)化前后最大接觸應力變化表/104kPaTab.3 Changes of maximum contact stress before and after optimization

    5 結論

    通過上述理論計算和數(shù)值分析的對比,可以得到如下幾點結論:

    (1)本文提出的優(yōu)化設計方法,適用于軸孔連接結構中軸的受力方向為單向受力且具備對孔進行精加工而間隙能夠精確控制的情況。對初始模型結構進行優(yōu)化設計后,軸孔接觸區(qū)的最大接觸應力降低50%以上,因此本文中的優(yōu)化設計方法可以用于提高軸孔連接結構承載能力的工程研究。

    (2)該方法從影響軸孔接觸區(qū)接觸應力大小的主要因素(接觸點數(shù)目及接觸點處的接觸曲率)出發(fā),通過改變孔接觸面的曲面形狀并優(yōu)化其曲面參數(shù),使得最初接觸點數(shù)目增多、接觸曲率得到優(yōu)化,從而使得外部載荷增加時軸孔的接觸順序從兩側逐步向中部擴展,增加了軸孔接觸面的分布范圍,降低了最大接觸應力。

    (3)與傳統(tǒng)的連接結構設計方法相比,本方法將接觸應力考慮到連接件的結構設計中,從連接件承載能力較低的根本原因出發(fā)來設計結構。利用此方法指導連接件設計,結構的利用率更高,對結構的設計工作有一定的推動作用。該設計思想也可以進一步推廣應用到其他有關連接結構中,具有廣泛的工程價值和應用前景。

    [1]楊東升,胡偉平,呂來清.銷釘耳片連接件接觸應力分析的能量解法[J].應用力學學報,2011,28(3):304-307. Yang Dongsheng,Hu Weiping,Lü Laiqing.The energy approach to contact stress analysis of pinned connections[J].Chinese Journal of Applied Mechanics,2011,28(3):304-307.(in Chinese)

    [2]Iyer K.Solutions for contact in pinned connections[J].International Journal of Solids and Structures,2001,38(50):9133-9148.

    [3]Yavari V,Rajabi I,Daneshvar F,et al.On the stress distribution around the hole in mechanically fastened joints[J].Mechanics Research Communications,2009,36(3):373-380.

    [4]顏東煌,劉雪鋒,田仲初,等.銷軸連接結構的接觸應力分析[J].工程力學,2008,25(1):229-234. Yan Donghuang,Liu Xuefeng,Tian Zhongchu,et al.Contact stress analysis of pinned connection[J].Engineering Mechanics,2008,25(1):229-234.(in Chinese)

    [5]楊生華.齒輪接觸有限元分析[J].計算力學學報,2003,20(2):189-193. Yang Shenghua.Finite element analysis of gear contact[J].Chinese Journal of Computational Mechanics,2003,20(2): 189-193.(in Chinese)

    [6]竇鵬,李友國,梁開明,汪長安.CVC熱軋機支承輥接觸應力有限元分析[J].清華大學學報,2005,45(12):1668-1671. Dou Peng,Li Youguo,Liang Kaiming,Wang Changan.Finite element analysis of contact stresses on the backup roll of CVC hot rolling mills[J].Journal of Tsinghua University,2005,45(12):1668-1671.(in Chinese)

    [7]魏大盛,王延榮.榫連結構接觸面幾何構形對接觸區(qū)應力分布的影響[J].航空動力學報,2010,25(2):407-411. Wei Dasheng,Wang Yanrong.Effects of profile of contact surfaces on the stress distribution for tenon jointing in blade disk assemblies[J].Journal of Aerospace Power,2010,25(2):407-411.(in Chinese)

    [8]Johnson K L.Contact mechanics[M].Cambridge:Cambridge University Press,1987.

    [9]Wriggers P,Zavarise G.Computational contact mechanics[M].John Wiley&Sons,Ltd,2002.

    [10]王寶忠.飛機設計手冊第10冊:結構設計[M].北京:航空工業(yè)出版社,1994:328-329. Wang Baozhong.Aircraft design manual vol10-structure design[M].Beijing:Aviation Industry Press,2000:410-411.(in Chinese)

    [11]Ciavarella M,Decuzzi P.The state of stress induced by the plane frictionless cylindrical contact.I.The case of elastic similarity[J].International Journal of Solids and Structures,2001,38(26):4507-4523.

    [12]Ciavarella M,Decuzzi P.The state of stress induced by the plane frictionless cylindrical contact.II.The general case(elastic dissimilarity)[J].International Journal of Solids and Structures,2001,38(26):4525-4533.

    [13]Sinclair G B,Cormier N G.Contact stresses in dovetail attachments:Physical modeling[J].Journal of Engineering for Gas Turbines and Power,2002,124(2):325-331.

    [14]Belytschko T,Liu W K,Moran B,et al.Nonlinear finite elements for continua and structures[M].John Wiley&Sons, 2013.

    Optimal design of improving the carrying capacity of shaft and hole connection structure

    LIU Wen-zhang1,HE Jing-Wu2,YAN Bin1,WU Wen-wei1
    (1.China Ship Scientific Research Center,Wuxi 214082,China;2.School of Aeronautic Science and Engineering,BeiHang University,Beijing 100191,China)

    In order to improve the contact characteristics of shaft and hole connection structure and improve its carrying capacity,one method was proposed.For a typical bolt and lug connection structure,the contact analysis model was established.By Persson contact theory and optimization method,the surface shape and its parameters of the contact surface of shaft and hole were optimized to lower the maximum contact stress and improve its distribution gradient.For the optimized design results of contact surface,using finite element analysis,the FEA model of bolt and lug connection structure was created to calculate the contact stress and verify the optimization results.The results show that the calculated results of analytical method and finite element method are basically the same.Through optimization,the maximum stress of the contact area is reduced by more than 50 percent and the stress distribution gradient is improved significantly.So the carrying capacity of shaft and hole connection structure is increased.

    shaft and hole connection structure;Persson contact theory;contact stress; structure optimization;carrying capacity

    O343.3

    A

    10.3969/j.issn.1007-7294.2017.08.009

    1007-7294(2017)08-1001-08

    2017-02-25

    劉文章(1987-),男,工程師,E-mail:articleliu@163.com;何景武(1961-),男,教授。

    猜你喜歡
    接觸區(qū)軸孔耳片
    非比例多軸載荷耳片的疲勞分析方法
    民用飛機耳片結構參數(shù)化設計方法研究
    基于三維模型耳片接頭幾何參數(shù)敏感性分析
    煉鋼轉(zhuǎn)爐托圈銷軸孔在線修復探討
    冶金設備(2021年2期)2021-07-21 08:44:28
    等高齒準雙曲面齒輪切齒控制方法的優(yōu)化試驗
    弧齒錐齒輪接觸斑點的試驗研究
    軸孔協(xié)調(diào)接觸建模與試驗研究
    含孔邊套筒的厚截面復合材料接頭耳片強度分析
    有關大跨距同軸孔系加工方法的研究
    接觸區(qū)中的跨文化接觸與交換
    av又黄又爽大尺度在线免费看| 日韩精品免费视频一区二区三区| 五月天丁香电影| 免费久久久久久久精品成人欧美视频| 国产黄色免费在线视频| 日韩电影二区| 91字幕亚洲| 国产精品国产三级国产专区5o| 超碰成人久久| 亚洲欧美日韩高清在线视频 | 亚洲欧美日韩高清在线视频 | 97精品久久久久久久久久精品| 99九九在线精品视频| 9色porny在线观看| 免费久久久久久久精品成人欧美视频| 亚洲欧美成人综合另类久久久| av网站免费在线观看视频| 久久鲁丝午夜福利片| 人人妻人人爽人人添夜夜欢视频| 国产亚洲av片在线观看秒播厂| 男女午夜视频在线观看| 日韩制服丝袜自拍偷拍| 亚洲国产精品国产精品| 久久精品国产a三级三级三级| 国产高清视频在线播放一区 | 少妇 在线观看| 亚洲午夜精品一区,二区,三区| 日本欧美国产在线视频| 成年人午夜在线观看视频| 婷婷色综合www| 99久久综合免费| 久久精品国产综合久久久| 午夜日韩欧美国产| 丝袜美足系列| 热re99久久国产66热| 最近手机中文字幕大全| 91九色精品人成在线观看| 亚洲av片天天在线观看| 亚洲成人手机| 日韩制服骚丝袜av| 美女主播在线视频| 国产精品av久久久久免费| 亚洲黑人精品在线| 欧美亚洲 丝袜 人妻 在线| 国产成人一区二区三区免费视频网站 | 午夜影院在线不卡| 男人添女人高潮全过程视频| 欧美av亚洲av综合av国产av| 久久精品aⅴ一区二区三区四区| 在线观看免费视频网站a站| 亚洲人成网站在线观看播放| 久热这里只有精品99| 久久精品熟女亚洲av麻豆精品| 国产精品香港三级国产av潘金莲 | 最新在线观看一区二区三区 | 国产在视频线精品| 老司机深夜福利视频在线观看 | 97人妻天天添夜夜摸| 国产精品国产三级国产专区5o| 超碰97精品在线观看| 欧美日韩黄片免| 麻豆乱淫一区二区| 亚洲精品国产区一区二| 亚洲 国产 在线| 嫩草影视91久久| 国产成人系列免费观看| 精品国产一区二区三区四区第35| 一区二区日韩欧美中文字幕| 人人妻人人爽人人添夜夜欢视频| 美女高潮到喷水免费观看| 久久精品国产a三级三级三级| 一区福利在线观看| 精品人妻熟女毛片av久久网站| 免费高清在线观看日韩| 精品免费久久久久久久清纯 | 国产亚洲午夜精品一区二区久久| 两个人看的免费小视频| 啦啦啦中文免费视频观看日本| 国产1区2区3区精品| 国产片特级美女逼逼视频| 建设人人有责人人尽责人人享有的| 国产精品一二三区在线看| 亚洲精品一区蜜桃| 午夜福利一区二区在线看| 精品人妻熟女毛片av久久网站| 18在线观看网站| 国产精品亚洲av一区麻豆| 一区二区三区四区激情视频| 亚洲少妇的诱惑av| 久久99一区二区三区| 2018国产大陆天天弄谢| 亚洲欧美精品综合一区二区三区| 在线精品无人区一区二区三| 极品人妻少妇av视频| 国语对白做爰xxxⅹ性视频网站| 丝袜喷水一区| 丝袜人妻中文字幕| 日日摸夜夜添夜夜爱| 好男人电影高清在线观看| 黑丝袜美女国产一区| 国产97色在线日韩免费| 亚洲精品第二区| 国产成人一区二区在线| 国产97色在线日韩免费| 99热全是精品| 午夜福利免费观看在线| 男女边吃奶边做爰视频| 中文字幕亚洲精品专区| 色网站视频免费| 首页视频小说图片口味搜索 | 免费人妻精品一区二区三区视频| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲| 电影成人av| 搡老乐熟女国产| 午夜两性在线视频| 日韩,欧美,国产一区二区三区| 日本一区二区免费在线视频| 高清欧美精品videossex| 国产xxxxx性猛交| 亚洲一区二区三区欧美精品| 午夜日韩欧美国产| 午夜免费男女啪啪视频观看| 日本av免费视频播放| 国产1区2区3区精品| 丁香六月天网| 又大又爽又粗| 精品高清国产在线一区| 午夜免费观看性视频| 国产精品国产三级专区第一集| 国产精品一区二区免费欧美 | 色94色欧美一区二区| 各种免费的搞黄视频| 9热在线视频观看99| 热99久久久久精品小说推荐| 国产麻豆69| 高清欧美精品videossex| 久久亚洲精品不卡| 日本av免费视频播放| 好男人视频免费观看在线| 天天影视国产精品| 欧美成狂野欧美在线观看| 亚洲成人手机| 精品人妻在线不人妻| 国产精品偷伦视频观看了| 国产1区2区3区精品| 国产精品麻豆人妻色哟哟久久| 天天躁夜夜躁狠狠久久av| 欧美久久黑人一区二区| 丁香六月天网| 国产精品一国产av| 我要看黄色一级片免费的| 午夜福利乱码中文字幕| 国产一区二区在线观看av| 色婷婷久久久亚洲欧美| 免费在线观看完整版高清| 国产成人av教育| 一本综合久久免费| 一本一本久久a久久精品综合妖精| 又大又黄又爽视频免费| 日韩 亚洲 欧美在线| 91精品三级在线观看| 精品一区二区三区四区五区乱码 | 久久女婷五月综合色啪小说| av天堂久久9| 三上悠亚av全集在线观看| 精品国产乱码久久久久久男人| 欧美少妇被猛烈插入视频| 亚洲精品中文字幕在线视频| 国产免费福利视频在线观看| 夫妻午夜视频| 男女午夜视频在线观看| 亚洲熟女精品中文字幕| 久久久久网色| 久久久国产精品麻豆| 午夜激情久久久久久久| 91字幕亚洲| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 精品人妻1区二区| 日韩一本色道免费dvd| 亚洲精品av麻豆狂野| 日本91视频免费播放| a级毛片在线看网站| 高清不卡的av网站| 考比视频在线观看| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 在线天堂中文资源库| 少妇的丰满在线观看| 激情视频va一区二区三区| 99精国产麻豆久久婷婷| 国产成人av激情在线播放| 国产野战对白在线观看| 日韩av在线免费看完整版不卡| 日日爽夜夜爽网站| 国产真人三级小视频在线观看| 久久精品亚洲av国产电影网| 五月天丁香电影| 久久精品国产亚洲av涩爱| 久久久久久亚洲精品国产蜜桃av| 悠悠久久av| 黑丝袜美女国产一区| 久久ye,这里只有精品| 自线自在国产av| 亚洲黑人精品在线| 一区二区三区乱码不卡18| 亚洲视频免费观看视频| 国产精品人妻久久久影院| 久久精品久久久久久久性| 97在线人人人人妻| 欧美成狂野欧美在线观看| 超碰成人久久| 国产成人精品在线电影| 欧美久久黑人一区二区| 少妇精品久久久久久久| 亚洲欧美成人综合另类久久久| 777久久人妻少妇嫩草av网站| 五月开心婷婷网| 亚洲国产中文字幕在线视频| 在线观看人妻少妇| 日本一区二区免费在线视频| 国产精品99久久99久久久不卡| 成人亚洲精品一区在线观看| 黄色a级毛片大全视频| 亚洲欧美色中文字幕在线| 久久精品亚洲av国产电影网| 九色亚洲精品在线播放| 老司机影院毛片| 亚洲精品国产色婷婷电影| 午夜免费成人在线视频| 国产熟女午夜一区二区三区| 99热全是精品| av网站免费在线观看视频| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| xxx大片免费视频| 老汉色∧v一级毛片| 欧美在线一区亚洲| 9191精品国产免费久久| 亚洲精品美女久久av网站| 一级片'在线观看视频| 亚洲精品国产av蜜桃| 国产精品久久久久久人妻精品电影 | 精品久久久久久久毛片微露脸 | 一级片'在线观看视频| 丝袜美足系列| 国产成人系列免费观看| 国产成人精品无人区| 多毛熟女@视频| 黄色视频在线播放观看不卡| 亚洲三区欧美一区| 一区福利在线观看| 国产精品久久久久久人妻精品电影 | 日韩av不卡免费在线播放| 欧美乱码精品一区二区三区| videosex国产| 亚洲,一卡二卡三卡| 国产精品久久久人人做人人爽| 菩萨蛮人人尽说江南好唐韦庄| 韩国高清视频一区二区三区| 在线看a的网站| 亚洲中文日韩欧美视频| 精品少妇黑人巨大在线播放| 久久国产精品人妻蜜桃| 久久毛片免费看一区二区三区| 一本色道久久久久久精品综合| 中文字幕制服av| 中国美女看黄片| 欧美性长视频在线观看| 一级毛片黄色毛片免费观看视频| av在线播放精品| 久久天堂一区二区三区四区| 成人国产av品久久久| 精品少妇内射三级| 岛国毛片在线播放| 午夜免费鲁丝| 高清视频免费观看一区二区| svipshipincom国产片| 精品亚洲成a人片在线观看| 麻豆国产av国片精品| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 一区二区三区四区激情视频| 国产在线视频一区二区| 黄片播放在线免费| 1024香蕉在线观看| 国产成人啪精品午夜网站| 国产野战对白在线观看| e午夜精品久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 一区在线观看完整版| 秋霞在线观看毛片| 又大又爽又粗| 一边摸一边做爽爽视频免费| 婷婷色麻豆天堂久久| 欧美性长视频在线观看| 精品免费久久久久久久清纯 | 天天躁夜夜躁狠狠久久av| 中文字幕精品免费在线观看视频| 欧美日韩亚洲综合一区二区三区_| 久久 成人 亚洲| 黄片小视频在线播放| 国产亚洲午夜精品一区二区久久| 在线观看免费视频网站a站| 一本—道久久a久久精品蜜桃钙片| 三上悠亚av全集在线观看| 青草久久国产| 建设人人有责人人尽责人人享有的| 人成视频在线观看免费观看| 9热在线视频观看99| 男女之事视频高清在线观看 | 十八禁网站网址无遮挡| 久久综合国产亚洲精品| 欧美激情 高清一区二区三区| 国产精品 欧美亚洲| 国产成人系列免费观看| 大码成人一级视频| 日韩视频在线欧美| xxx大片免费视频| 91成人精品电影| 18禁黄网站禁片午夜丰满| 亚洲精品在线美女| 午夜av观看不卡| 久久人人爽人人片av| 伊人亚洲综合成人网| 超碰成人久久| 国产成人精品久久二区二区91| 国产在视频线精品| 男女之事视频高清在线观看 | av有码第一页| 亚洲,一卡二卡三卡| 在线精品无人区一区二区三| 精品一区二区三区av网在线观看 | 妹子高潮喷水视频| 波多野结衣一区麻豆| 久久久精品区二区三区| 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 欧美成狂野欧美在线观看| 亚洲熟女毛片儿| 丝袜人妻中文字幕| 性色av一级| 欧美日韩成人在线一区二区| 免费看av在线观看网站| 日本午夜av视频| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看| 午夜两性在线视频| 制服诱惑二区| 亚洲精品国产av蜜桃| 最黄视频免费看| 久久久久久久大尺度免费视频| 亚洲美女黄色视频免费看| 国产精品国产三级国产专区5o| 视频在线观看一区二区三区| 国产亚洲一区二区精品| 欧美日韩视频精品一区| 久久久久久久国产电影| 极品人妻少妇av视频| 久久精品国产a三级三级三级| 亚洲av成人精品一二三区| 91字幕亚洲| 亚洲欧洲国产日韩| 极品少妇高潮喷水抽搐| 一本—道久久a久久精品蜜桃钙片| 丁香六月天网| 啦啦啦在线观看免费高清www| 国产免费又黄又爽又色| 夫妻午夜视频| 国产伦理片在线播放av一区| 日本av手机在线免费观看| 十八禁网站网址无遮挡| 啦啦啦视频在线资源免费观看| 国产主播在线观看一区二区 | 国产av一区二区精品久久| 亚洲精品国产一区二区精华液| 日韩精品免费视频一区二区三区| 久久久久精品国产欧美久久久 | 性高湖久久久久久久久免费观看| 精品高清国产在线一区| 久久中文字幕一级| 精品久久久精品久久久| 亚洲av日韩在线播放| a级毛片黄视频| 国产精品久久久久久人妻精品电影 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| 久久久精品94久久精品| 欧美黄色片欧美黄色片| 国产99久久九九免费精品| a级片在线免费高清观看视频| 青春草亚洲视频在线观看| 18禁国产床啪视频网站| av国产精品久久久久影院| 岛国毛片在线播放| 女人爽到高潮嗷嗷叫在线视频| 国产有黄有色有爽视频| 手机成人av网站| 国产精品 国内视频| 一边摸一边做爽爽视频免费| 免费观看a级毛片全部| 国产欧美日韩一区二区三区在线| 免费不卡黄色视频| 男人添女人高潮全过程视频| 久久ye,这里只有精品| 国产成人一区二区在线| 在线观看免费高清a一片| 免费高清在线观看日韩| www日本在线高清视频| 纵有疾风起免费观看全集完整版| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区免费开放| 伊人久久大香线蕉亚洲五| 亚洲伊人久久精品综合| 青春草视频在线免费观看| 欧美国产精品一级二级三级| 欧美黄色片欧美黄色片| 美女福利国产在线| 免费女性裸体啪啪无遮挡网站| 美女国产高潮福利片在线看| 亚洲五月婷婷丁香| 午夜精品国产一区二区电影| 七月丁香在线播放| 日韩精品免费视频一区二区三区| 亚洲av电影在线观看一区二区三区| 99re6热这里在线精品视频| 亚洲色图综合在线观看| 国产有黄有色有爽视频| 97人妻天天添夜夜摸| 久久99精品国语久久久| 中文字幕制服av| 免费不卡黄色视频| 国产主播在线观看一区二区 | 99精国产麻豆久久婷婷| 国产精品久久久久久精品电影小说| 中国美女看黄片| 999精品在线视频| 久久久精品国产亚洲av高清涩受| 久久久精品免费免费高清| 大片电影免费在线观看免费| 国产91精品成人一区二区三区 | 搡老乐熟女国产| 久久久久精品人妻al黑| 亚洲国产最新在线播放| 99国产精品一区二区蜜桃av | 亚洲伊人久久精品综合| 精品视频人人做人人爽| 91老司机精品| 亚洲成色77777| 欧美性长视频在线观看| 免费人妻精品一区二区三区视频| 国产欧美日韩一区二区三区在线| 中文字幕色久视频| 宅男免费午夜| 91精品伊人久久大香线蕉| 久久人妻熟女aⅴ| 九色亚洲精品在线播放| 中文字幕人妻熟女乱码| 久久久国产欧美日韩av| 国产人伦9x9x在线观看| videosex国产| 亚洲午夜精品一区,二区,三区| 免费人妻精品一区二区三区视频| 男人舔女人的私密视频| 国产一区亚洲一区在线观看| 国产亚洲精品第一综合不卡| av片东京热男人的天堂| 搡老岳熟女国产| 久久久久精品国产欧美久久久 | 国产又爽黄色视频| 亚洲情色 制服丝袜| 国产成人精品在线电影| 欧美日韩精品网址| 国产熟女午夜一区二区三区| 少妇精品久久久久久久| 国产日韩一区二区三区精品不卡| 国产免费现黄频在线看| 一边摸一边抽搐一进一出视频| 一区二区三区精品91| 日韩制服丝袜自拍偷拍| 午夜久久久在线观看| 深夜精品福利| 午夜福利影视在线免费观看| 18在线观看网站| 国产爽快片一区二区三区| 亚洲精品乱久久久久久| 国产伦理片在线播放av一区| 99九九在线精品视频| 电影成人av| 欧美日韩亚洲高清精品| 可以免费在线观看a视频的电影网站| 免费看不卡的av| 亚洲精品av麻豆狂野| 日本a在线网址| 建设人人有责人人尽责人人享有的| 国产精品秋霞免费鲁丝片| 男女无遮挡免费网站观看| 天堂8中文在线网| 免费高清在线观看日韩| 少妇 在线观看| 18禁黄网站禁片午夜丰满| 男人添女人高潮全过程视频| 免费在线观看日本一区| 51午夜福利影视在线观看| 一边摸一边做爽爽视频免费| 性少妇av在线| 91麻豆精品激情在线观看国产 | 丝袜喷水一区| 亚洲国产精品一区三区| 成在线人永久免费视频| 一边摸一边做爽爽视频免费| 性少妇av在线| 国产91精品成人一区二区三区 | 日日夜夜操网爽| 亚洲七黄色美女视频| 精品国产超薄肉色丝袜足j| 国产成人免费观看mmmm| 亚洲第一青青草原| 日韩免费高清中文字幕av| 十八禁人妻一区二区| 日本五十路高清| 欧美精品av麻豆av| 天天添夜夜摸| 又紧又爽又黄一区二区| 久久久精品国产亚洲av高清涩受| 考比视频在线观看| 亚洲av片天天在线观看| 最新的欧美精品一区二区| 亚洲国产精品国产精品| 久久午夜综合久久蜜桃| 欧美日韩亚洲综合一区二区三区_| 一边摸一边做爽爽视频免费| 国产97色在线日韩免费| 深夜精品福利| 免费日韩欧美在线观看| 99精品久久久久人妻精品| 男女无遮挡免费网站观看| 亚洲精品av麻豆狂野| 欧美在线一区亚洲| 天天躁夜夜躁狠狠躁躁| 国产精品三级大全| 视频区欧美日本亚洲| 一二三四在线观看免费中文在| 精品国产一区二区久久| 成人午夜精彩视频在线观看| 9191精品国产免费久久| 免费人妻精品一区二区三区视频| tube8黄色片| 精品少妇久久久久久888优播| 亚洲av日韩在线播放| 中文字幕制服av| 多毛熟女@视频| 亚洲av电影在线进入| 一级黄片播放器| 操美女的视频在线观看| 国产91精品成人一区二区三区 | 国产成人a∨麻豆精品| 日韩av在线免费看完整版不卡| 亚洲av在线观看美女高潮| 国产高清国产精品国产三级| 欧美激情高清一区二区三区| 老司机靠b影院| 国产激情久久老熟女| 亚洲成色77777| 日韩一区二区三区影片| 中文字幕色久视频| 国产在线免费精品| 中文字幕av电影在线播放| 母亲3免费完整高清在线观看| www.自偷自拍.com| 国产人伦9x9x在线观看| 免费观看av网站的网址| 老司机靠b影院| 男女边吃奶边做爰视频| 最近中文字幕2019免费版| 妹子高潮喷水视频| 一边亲一边摸免费视频| 成人免费观看视频高清| 国产视频一区二区在线看| 在线av久久热| 成年人免费黄色播放视频| 日韩 亚洲 欧美在线| 国产成人免费观看mmmm| 国产欧美日韩精品亚洲av| 丝袜脚勾引网站| 国产精品秋霞免费鲁丝片| 亚洲,欧美精品.| 老司机影院成人| 亚洲人成电影免费在线| 久久天堂一区二区三区四区| 超碰成人久久| 人人妻人人澡人人看| 热99久久久久精品小说推荐| 中文字幕精品免费在线观看视频| 亚洲一卡2卡3卡4卡5卡精品中文| 80岁老熟妇乱子伦牲交| 亚洲欧美一区二区三区久久| 波多野结衣一区麻豆| av视频免费观看在线观看| 午夜福利影视在线免费观看| 中文字幕人妻熟女乱码| kizo精华| 老司机亚洲免费影院| 久久久精品国产亚洲av高清涩受| 人妻一区二区av| a 毛片基地| 国产97色在线日韩免费| 成年人免费黄色播放视频| 亚洲五月色婷婷综合| 亚洲成国产人片在线观看| 在线天堂中文资源库| 国产视频一区二区在线看| 亚洲综合色网址| 久久久国产精品麻豆|