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

    形狀記憶合金三維鼓包結(jié)構(gòu)熱變形的有限元模擬

    2014-12-09 09:06:06趙澎濤裘進(jìn)浩于慧臣
    機(jī)械工程材料 2014年3期
    關(guān)鍵詞:鼓包本構(gòu)馬氏體

    趙澎濤,裘進(jìn)浩,于慧臣

    (1.北京航空材料研究院航空材料檢測(cè)與評(píng)價(jià)北京市重點(diǎn)實(shí)驗(yàn)室,先進(jìn)高溫結(jié)構(gòu)材料重點(diǎn)實(shí)驗(yàn)室,北京100095;2.南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京210016)

    0 引 言

    形狀記憶合金(SMA)同時(shí)具備感知和驅(qū)動(dòng)功能,可以感知外界的溫度變化,并能夠?qū)崮苻D(zhuǎn)化為機(jī)械能,其以優(yōu)異的形狀記憶效應(yīng)、偽彈性[1-2]、生物相容性和阻尼特性,廣泛應(yīng)用于航空航天[3-4]、生物醫(yī)藥[5]、土木工程[6]等領(lǐng)域。但由于SMA 復(fù)雜的熱力學(xué)耦合關(guān)系,其本構(gòu)關(guān)系不能用通常的彈塑性模型描述,因此SMA 本構(gòu)行為的描述是目前研究的熱點(diǎn)和難點(diǎn)。為了實(shí)現(xiàn)SMA更廣泛的工程應(yīng)用,作者以SMA在三維自適應(yīng)發(fā)動(dòng)機(jī)進(jìn)氣道變體結(jié)構(gòu)中的應(yīng)用為研究背景,利用SMA所特有的形狀記憶效應(yīng),使用溫度控制的方式實(shí)現(xiàn)了鼓包結(jié)構(gòu)的變形。通過(guò)建立SMA三維本構(gòu)關(guān)系模型,采用回歸映射應(yīng)力更新算法進(jìn)行有限元實(shí)現(xiàn),并對(duì)SMA三維鼓包結(jié)構(gòu)在熱變形過(guò)程中的熱力學(xué)特性響應(yīng)進(jìn)行有限元分析,得到了一定溫度條件下該結(jié)構(gòu)的變形參數(shù),為SMA在發(fā)動(dòng)機(jī)進(jìn)氣道中的應(yīng)用提供支持。

    1 SMA的本構(gòu)關(guān)系模型

    采用Lagoudas關(guān)于SMA的三維唯象理論模型[7-8],即SMA在變形過(guò)程中產(chǎn)生的應(yīng)變主要由彈性應(yīng)變、熱膨脹應(yīng)變和相變引起的應(yīng)變等三部分構(gòu)成,假設(shè)SMA的Gibbs自由能(G)為

    式中:σ,εt,ξ,T 以及T0分別為真實(shí)應(yīng)力張量、相變應(yīng)變張量、馬氏體體積分?jǐn)?shù)、溫度以及參考溫度;S,α,ρ,c,s0和u0分別為等效柔度張量、等效熱膨脹系數(shù)、密度、等效比熱、參考溫度下的等效比熵和參考溫度下的等效內(nèi)能;f(ξ)為相變硬化函數(shù),表征了馬氏體相與母相或者馬氏體變體之間發(fā)生相互轉(zhuǎn)化而產(chǎn)生的應(yīng)變能。

    以上各等效材料特性參數(shù)可由材料中所含的馬氏體體積分?jǐn)?shù)以及純晶體相的特性參數(shù)計(jì)算得到,其具體的表達(dá)式如式(2)。式中物理量右上角處的A與M分別表示為奧氏體相和馬氏體相:

    結(jié)合熱力學(xué)相關(guān)理論,SMA的本構(gòu)關(guān)系及約束條件可以從Gibbs自由能推導(dǎo)得到,其表達(dá)式如式(3)和式(4):

    假設(shè)在馬氏體相變過(guò)程中忽略馬氏體變體的重定向過(guò)程,材料任何微觀狀態(tài)的改變都僅僅是馬氏體含量變化的結(jié)果。基于此假設(shè),SMA相變過(guò)程中的相變應(yīng)變?chǔ)舤和馬氏體體積分?jǐn)?shù)ξ的演化方程(類似于塑性力學(xué)理論中的流動(dòng)法則)可表示為

    式中:Λ為相變轉(zhuǎn)換張量,表征了相變應(yīng)變的方向。

    將相變演化方程式(5)代入式(4)的本構(gòu)關(guān)系約束條件中,可得到:

    式中:π為與馬氏體體積分?jǐn)?shù)ξ共軛的熱驅(qū)動(dòng)力。

    將Gibbs自由能的顯示表達(dá)式(1)代入上式,可得到:

    基于以上關(guān)于馬氏體正逆相變的假設(shè)以及邊界條件的表述,引入相變函數(shù)φ(相當(dāng)于塑性力學(xué)中的屈服函數(shù)),其表達(dá)式為

    式中:Y為相變臨界屈服值。

    基于不同相變過(guò)程的假設(shè),SMA唯象理論模型可以通過(guò)不同的相變硬化函數(shù)來(lái)表征。Tanaka指數(shù)形式理論模型[9]、Liang和 Rogers余弦模型[10]以及Lagoudas多項(xiàng)式模型[7]的硬化函數(shù)分別見(jiàn)式(9),(10),(11):

    2 UMAT材料子程序的編寫(xiě)及驗(yàn)證

    在Lagoudas的SMA三維本構(gòu)關(guān)系模型的基礎(chǔ)上,利用ABAQUS的二次開(kāi)發(fā)功能編寫(xiě)FORTRAN程序[11]實(shí)現(xiàn)SMA的有限元模擬。

    針對(duì)前面提到的SMA本構(gòu)關(guān)系模型,在UMAT中定義了13個(gè)變量(1為馬氏體體積分?jǐn)?shù)、2~7為相變應(yīng)變分量,8~13為個(gè)相變方向)和20個(gè)材料常數(shù)(1~2為彈性模量,3~4為熱膨脹系數(shù),5~8為相變溫度,9~10為應(yīng)力影響參數(shù),11~16為初始相變應(yīng)變,17為最大相變應(yīng)變,18為泊松比,19為容差,20為初始馬氏體體積分?jǐn)?shù))。

    為了驗(yàn)證編寫(xiě)的UMAT子程序的準(zhǔn)確性,在ABAQUS軟件中建立受單軸加載的單元模型,如圖1所示,一端固支約束,一端加載軸向力載荷。載荷加載歷史如圖2所示。模型中的SMA材料參數(shù)采用文獻(xiàn)[12]中的等原子鈦鎳合金參數(shù),即:奧氏體相彈性模量EA為70GPa,馬氏體相彈性模量EM為30GPa,泊松比νA和νM均為0.3,奧氏體相膨脹系數(shù)αA為22×10-6K-1,馬氏體相膨脹系數(shù)αM為22×10-6K-1,馬氏體相變開(kāi)始溫度Ms為291K(18℃),馬氏體相變結(jié)束溫度Mf為280K(7℃),奧氏體相變開(kāi)始溫度As為295K(22℃),奧氏體相變結(jié)束溫度Af為306K(33℃),最大相變應(yīng)變H為0.05。針對(duì)超彈性的驗(yàn)證,模型初始相為奧氏體,初始溫度為325K(52℃),單軸最大載荷為60kN。針對(duì)形狀記憶效應(yīng)的驗(yàn)證,模型初始相為馬氏體,初始溫度為292K(19℃),單軸最大載荷為30kN,加載結(jié)束后升溫至350K(77℃),高于奧氏體相變結(jié)束溫度。

    圖1 單軸加載模型示意Fig.1 Abridged general view of uniaxial loading model

    圖2 載荷加載歷史示意Fig.2 Abridged general view of applied load history

    從圖3,4可以看出,以上基于不同硬化函數(shù)的三種本構(gòu)關(guān)系模型都能夠較好地描述SMA的超彈性和形狀記憶效應(yīng)等特性,相變過(guò)程也遵循相應(yīng)的馬氏體奧氏體相變模型表征函數(shù),證明了該UMAT子程序的正確性,從而保證了在UMAT子程序基礎(chǔ)上進(jìn)行的鼓包熱變形模擬結(jié)果的準(zhǔn)確。

    圖3 基于不同硬化函數(shù)的SMA超彈性模型的模擬結(jié)果Fig.3 Simulated results of SMA pesudoelasticity model based different hardening functions

    圖4 基于不同硬化函數(shù)的SMA形狀記憶效應(yīng)的模擬結(jié)果Fig.4 Simulated results of SMA shape memory effect based different hardening fuctions

    3 有限元模擬結(jié)果與討論

    3.1 有限元模型的建立

    SMA鼓包熱變形模型是三維薄壁球體結(jié)構(gòu)的一部分,其尺寸參數(shù)包括:過(guò)中心點(diǎn)A的弧線在x-z平面上的投影長(zhǎng)(圖5中的BC線段,10cm)、鼓包厚度(1mm)、曲率半徑(42cm)。該鼓包模型在經(jīng)歷雙程記憶效應(yīng)訓(xùn)練之后,沿過(guò)中心點(diǎn)的徑長(zhǎng)方向具有高溫伸長(zhǎng)、低溫收縮的性質(zhì)。模型網(wǎng)格劃分采用了減縮積分的線性六面體單元(C3D8R),在厚度方向上劃分5個(gè)單元,以免產(chǎn)生沙漏效應(yīng)。材料參數(shù)與前面驗(yàn)證模型中采用的參數(shù)相同。

    圖5 三維鼓包熱變形模型的網(wǎng)格劃分Fig.5 Meshing of the three demensional bump thermal deformation model

    鼓包熱變形模型邊界條件為周邊固支約束,初始溫度為291K(18℃),初始狀態(tài)為馬氏體相。沿過(guò)中心點(diǎn)的弧長(zhǎng)方向具有初始的相變應(yīng)變,根據(jù)不同的訓(xùn)練效果分為-0.5%,-1%,-1.5%,-2%幾種情況。考慮鼓包沿厚度方向的初始相變應(yīng)變,初始馬氏體體積分?jǐn)?shù)為1。結(jié)構(gòu)的網(wǎng)格劃分如圖5所示。SMA鼓包表面加載溫度載荷,從初始溫度291K(18℃)加載至360K(87℃)并保持溫度,增量加載步數(shù)為100步,具體過(guò)程如圖6所示。

    圖6 溫度加載歷程示意Fig.6 Abridged general view of temperature loading history

    3.2 模擬結(jié)果與討論

    3.2.1 定初始相變的SMA鼓包相變特性

    從圖7中可以看出,鼓包模型在初始狀態(tài)時(shí)為純馬氏體相,體積分?jǐn)?shù)為1;隨著溫度的逐漸升高,在鼓包中部馬氏體相的體積分?jǐn)?shù)開(kāi)始下降,并逐步向鼓包的四周邊界擴(kuò)展;加載結(jié)束后鼓包的大部分區(qū)域基本完成了馬氏體逆相變,轉(zhuǎn)化為奧氏體相,只有在接近固支約束的邊界部分仍殘留有馬氏體相。

    從圖8中可以看出,圖5中的A、B兩點(diǎn)在297K(24℃)左右開(kāi)始向奧氏體相轉(zhuǎn)變,A點(diǎn)在313K(40℃)結(jié)束了馬氏體逆相變過(guò)程,馬氏體體積分?jǐn)?shù)為0,全部轉(zhuǎn)化為奧氏體相;而在B點(diǎn),由于邊界條件的約束,限制了鼓包的馬氏體逆相變過(guò)程,在加載結(jié)束后仍殘留有部分馬氏體相。

    3.2.2 定初始相變應(yīng)變的SMA鼓包熱變形力學(xué)特性

    由圖9可見(jiàn),A點(diǎn)的豎向撓度在初始階段增長(zhǎng)較慢,在23~40℃(296~313K)的溫度區(qū)間增長(zhǎng)較快,并且大部分的撓度變化發(fā)生在該區(qū)間內(nèi),而且該區(qū)間與馬氏體逆相變的溫度區(qū)間相近??芍?,鼓包豎向撓度的產(chǎn)生主要與SMA的相變有關(guān)。

    3.2.3 不同初始相變的SMA鼓包熱變形結(jié)果

    鼓包在加載過(guò)程中最大豎向撓度發(fā)生在加載結(jié)束溫度360K。由圖10可見(jiàn),最大豎向撓度與初始相變應(yīng)變絕對(duì)值之間基本呈線性正比關(guān)系。當(dāng)初始相變應(yīng)變絕對(duì)值為2%時(shí),鼓包最大的豎向撓度達(dá)到了5mm,相當(dāng)于鼓包模型BC線段長(zhǎng)度的5%。

    圖7 初始相變應(yīng)變?yōu)?1%時(shí),SMA鼓包模型在不同溫度下的馬氏體體積分?jǐn)?shù)云圖Fig.7 Martensite fraction of SMA bump model at initial phase-transformation strain of-1%and different temperatures

    圖8 SMA鼓包模型中不同位置處在不同溫度下的馬氏體體積分?jǐn)?shù)Fig.8 Volume fractions of martensite vs temperture for different locations in SMA bump model

    圖9 SMA鼓包中心A點(diǎn)處的豎向撓度隨溫度的變化Fig.9 Change of vertical deflection of the point A(in the middle of the SMA bump)with temperature:(a)vertical deflection of bump and(b)vertical deflection of point A

    圖10 SMA鼓包在加載過(guò)程中最大豎向撓度與初始相變應(yīng)變絕對(duì)值之間的關(guān)系Fig.10 The relation ofmaximum vertical deflection and initial phase-transformation strain during SMA bump

    從圖11中可以看出,在一定的初始相變應(yīng)變下,隨著溫度的升高,豎向高度曲線的形狀和最大高度也發(fā)生改變;豎向高度在313K(40℃)之后的變化較小,該溫度與鼓包中A點(diǎn)相變結(jié)束點(diǎn)的溫度基本吻合。

    圖11 SMA鼓包過(guò)A點(diǎn)弧線各點(diǎn)在不同初始相變下的豎向高度曲線Fig.11 Vertical height curves of the points which pass by the point A at different initial phase-transformation strains

    3 結(jié) 論

    (1)通過(guò)ABAQUS軟件材料子程序UMAT的接口,編寫(xiě)FORTRAN程序,將SMA的本構(gòu)關(guān)系模型導(dǎo)入到ABAQUS材料庫(kù)中,從而使有限元軟件具備了對(duì)SMA模型進(jìn)行數(shù)值模擬的能力;通過(guò)建立SMA單軸受力模型,證明了材料子程序UAMT能夠同時(shí)描述SMA的超彈性和形狀記憶效應(yīng)特性。

    (2)對(duì)于三維鼓包熱變形模型,馬氏體逆相變首先發(fā)生在鼓包的中部并向四周擴(kuò)展,除邊界部分未完成逆相變外,其余大部分完全轉(zhuǎn)化為奧氏體相;鼓包的撓度最高達(dá)到5mm,相當(dāng)于其弧線投影長(zhǎng)度的5%。

    致謝:

    感謝南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室的鼎力支持。

    [1]MIRZAEIFAR R,SHAKERI M,SADIGHI M.Nonlinear finite element formulation for analyzing shape memory alloy cylindrical panels[J].Smart Materials and Structures,2009,18(3):1-14.

    [2]LAGOUDAS D C.Shape memory alloys:modeling and engineering applications[M].New York:Springer,2008:1-43.

    [3]GAO X,TURNER T L,BURTON D,et al.Finite element analysis of adaptive-stiffening and shape-control SMA hybrid composites[J].Journal of Engineering Materials and Technology,2006,128(3):285-293.

    [4]REY N,TILLMAN G,MILLER R,et al.Shape memory alloy actuation for a variable area fan nozzle[C]//Proceedings of SPIE,Smart Structures and Materials:Industrial and Commercial Application of Smart Structures Technologies.[S.l.]:[s.n],2001:371-382.

    [5]徐祖耀,江伯鴻.形狀記憶材料[M].上海:上海交通大學(xué)出版社,2000:3-6.

    [6]WILLIAMS K,CHIU G,BERNHARD R.Adaptive-passive absorbers using shape-memory alloys[J].Journal of Sound and Vibration,2002,249(5):835-848.

    [7]BO Z,LAGOUDAS D C.Thermomechanical modeling of polycrystalline SMAs under cyclic loading-Part I:theoretical derivations[J].International Journal of Engineering Science,1999,37:1089-1140.

    [8]QIDWAI M A,LAGOUDAS D C.On thermomechanics and transformation surfaces of polycrystalline NiTi shape memory alloy material[J].International Journal of Plasticity,2000,16:1309-1343.

    [9]TANAKA K,NISHIMURA F,HAYASHI T,et al.Phenomenological analysis on subloops and cyclic behavior in shape memory alloys under mechanical and/or thermal loads[J].Mechanics of Materials,1995,19:281-92.

    [10]LIANG C,ROGERS C A.The multi-dimensional constitutive relations of shape memory alloys[J].Journal of Engineering Mathematics,1992,26:429-443.

    [11]莊茁,張帆,岑松,等.ABAQUS非線性有限元分析與實(shí)例[M].北京:科學(xué)出版社,2005.

    [12]LAGOUDAS D C,BO Z,QIDWAI M A.A unified thermodynamic constitutive model for SMA and finite element analysis of active metal matrix composite[J].Mechanics of Composite Materials & Structures,1996,3:153-79.

    猜你喜歡
    鼓包本構(gòu)馬氏體
    中低碳系列馬氏體不銹鋼開(kāi)發(fā)與生產(chǎn)
    山東冶金(2022年1期)2022-04-19 13:40:20
    激光制備預(yù)壓應(yīng)力超高強(qiáng)韌馬氏體層的組織與性能
    某滲濾液收集池底部防滲層鼓包的分析與治理
    一起鍋爐清洗后鼓包事故的原因分析及處理
    離心SC柱混凝土本構(gòu)模型比較研究
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    一種新型超固結(jié)土三維本構(gòu)模型
    Fe-C-Mn-Si-Cr的馬氏體開(kāi)始轉(zhuǎn)變點(diǎn)的熱力學(xué)計(jì)算
    上海金屬(2016年1期)2016-11-23 05:17:24
    換流站電阻冷卻面板鼓包分析
    湖北電力(2016年8期)2016-05-09 09:03:08
    關(guān)于Fe-1.4C合金馬氏體相變的研究
    搡老熟女国产l中国老女人| 又爽又黄a免费视频| 综合色av麻豆| 国产精品一区二区性色av| 国产精品一区二区免费欧美| 国产欧美日韩一区二区精品| av免费在线看不卡| 亚洲国产高清在线一区二区三| 日本三级黄在线观看| 日日摸夜夜添夜夜爱| 婷婷六月久久综合丁香| 国产精品亚洲一级av第二区| 久久国内精品自在自线图片| 欧美日韩在线观看h| 精品久久久久久久久久久久久| 毛片一级片免费看久久久久| 99久久精品热视频| av.在线天堂| 国产色爽女视频免费观看| 非洲黑人性xxxx精品又粗又长| 欧美日韩乱码在线| 久久精品久久久久久噜噜老黄 | 色噜噜av男人的天堂激情| 在线观看66精品国产| 国产精品一区二区免费欧美| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品成人综合色| 久久久成人免费电影| 在线a可以看的网站| 小蜜桃在线观看免费完整版高清| 99热全是精品| 综合色av麻豆| 草草在线视频免费看| 国产女主播在线喷水免费视频网站 | 久久人人精品亚洲av| 成人国产麻豆网| 熟女电影av网| 欧美性猛交╳xxx乱大交人| 免费观看精品视频网站| 在线看三级毛片| 一夜夜www| 丰满乱子伦码专区| 亚洲一级一片aⅴ在线观看| 2021天堂中文幕一二区在线观| 午夜福利在线观看免费完整高清在 | av中文乱码字幕在线| 身体一侧抽搐| 久久久久久久亚洲中文字幕| 久久久久久久久久久丰满| 两个人视频免费观看高清| 精品人妻视频免费看| 日韩精品青青久久久久久| 久久久久免费精品人妻一区二区| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 午夜久久久久精精品| 欧美色视频一区免费| 精品免费久久久久久久清纯| 女生性感内裤真人,穿戴方法视频| 99国产极品粉嫩在线观看| 国产白丝娇喘喷水9色精品| 99久久精品一区二区三区| 久久精品国产99精品国产亚洲性色| 插阴视频在线观看视频| 国产在线精品亚洲第一网站| 人妻夜夜爽99麻豆av| 免费不卡的大黄色大毛片视频在线观看 | 男女那种视频在线观看| 国产免费一级a男人的天堂| 精品乱码久久久久久99久播| 日本欧美国产在线视频| 黄色欧美视频在线观看| 永久网站在线| 欧美一区二区国产精品久久精品| 大又大粗又爽又黄少妇毛片口| 你懂的网址亚洲精品在线观看 | 国产午夜精品久久久久久一区二区三区 | 国模一区二区三区四区视频| 国产在视频线在精品| .国产精品久久| 草草在线视频免费看| 69av精品久久久久久| 免费看日本二区| 国产成人一区二区在线| 91av网一区二区| 91在线精品国自产拍蜜月| 成人鲁丝片一二三区免费| 91狼人影院| 午夜福利在线在线| 欧美zozozo另类| 精品日产1卡2卡| 久久精品国产亚洲av涩爱 | 亚洲精品粉嫩美女一区| 三级毛片av免费| 成人精品一区二区免费| 日韩欧美精品v在线| 俄罗斯特黄特色一大片| 男人狂女人下面高潮的视频| 免费av不卡在线播放| 99久久精品热视频| 国内精品宾馆在线| 淫妇啪啪啪对白视频| 欧美日韩一区二区视频在线观看视频在线 | 日韩一本色道免费dvd| 亚洲内射少妇av| 亚洲国产色片| 亚洲国产日韩欧美精品在线观看| eeuss影院久久| 最近中文字幕高清免费大全6| 久久草成人影院| 成人高潮视频无遮挡免费网站| 搡老岳熟女国产| 一个人免费在线观看电影| 国产精品一二三区在线看| 中文资源天堂在线| 国产高清激情床上av| 天天躁日日操中文字幕| 天堂网av新在线| 亚洲欧美精品自产自拍| 久久久久久国产a免费观看| 精品人妻熟女av久视频| 欧美人与善性xxx| 91av网一区二区| 亚洲aⅴ乱码一区二区在线播放| 亚洲av电影不卡..在线观看| 国产精品亚洲一级av第二区| 九色成人免费人妻av| 能在线免费观看的黄片| 两性午夜刺激爽爽歪歪视频在线观看| 久久这里只有精品中国| 国内精品宾馆在线| 午夜免费男女啪啪视频观看 | 亚洲精品乱码久久久v下载方式| 国产精品一区二区三区四区免费观看 | 久久久色成人| 国产伦一二天堂av在线观看| 精品一区二区免费观看| 久久久色成人| 特级一级黄色大片| 久久精品国产清高在天天线| 熟女人妻精品中文字幕| 天天一区二区日本电影三级| 99久久中文字幕三级久久日本| 亚洲精品国产av成人精品 | 淫秽高清视频在线观看| 精品99又大又爽又粗少妇毛片| 亚洲国产精品合色在线| 深夜精品福利| 成年免费大片在线观看| 欧美最新免费一区二区三区| 亚洲精品粉嫩美女一区| 男人狂女人下面高潮的视频| 免费观看在线日韩| 深夜精品福利| 美女内射精品一级片tv| 大又大粗又爽又黄少妇毛片口| 久久久久国产网址| 亚洲欧美成人综合另类久久久 | 国产高清有码在线观看视频| 欧美一级a爱片免费观看看| ponron亚洲| 中文字幕熟女人妻在线| 在线播放国产精品三级| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 联通29元200g的流量卡| 3wmmmm亚洲av在线观看| 人妻少妇偷人精品九色| 网址你懂的国产日韩在线| 黑人高潮一二区| a级毛色黄片| h日本视频在线播放| 亚洲熟妇熟女久久| 中文字幕免费在线视频6| av.在线天堂| 一级毛片电影观看 | 日韩av在线大香蕉| 午夜日韩欧美国产| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 熟女电影av网| 91在线精品国自产拍蜜月| 真实男女啪啪啪动态图| 国产精品无大码| 久久久久国内视频| 国产亚洲精品av在线| 3wmmmm亚洲av在线观看| 美女 人体艺术 gogo| 久久久久国产网址| 国产黄a三级三级三级人| 国产美女午夜福利| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 欧美日韩一区二区视频在线观看视频在线 | av专区在线播放| 日日啪夜夜撸| 亚洲精品国产av成人精品 | 97热精品久久久久久| 国产精品电影一区二区三区| 国产精品av视频在线免费观看| 在现免费观看毛片| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| 日本黄色片子视频| 国产精品一区二区性色av| 最近视频中文字幕2019在线8| 国产精品永久免费网站| 国产精品无大码| 免费观看的影片在线观看| 中国美白少妇内射xxxbb| 一进一出好大好爽视频| 一级av片app| 99视频精品全部免费 在线| 久久精品影院6| 一级av片app| 国内精品美女久久久久久| 综合色丁香网| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 麻豆国产av国片精品| 国产毛片a区久久久久| 波多野结衣高清作品| 成人鲁丝片一二三区免费| 99视频精品全部免费 在线| 听说在线观看完整版免费高清| 亚洲图色成人| 无遮挡黄片免费观看| 成人无遮挡网站| 人人妻人人看人人澡| 国产蜜桃级精品一区二区三区| 韩国av在线不卡| 女同久久另类99精品国产91| 淫秽高清视频在线观看| 久久久国产成人精品二区| 亚洲在线观看片| 日韩av在线大香蕉| 亚洲精品久久国产高清桃花| 99热这里只有精品一区| 国产av麻豆久久久久久久| 1000部很黄的大片| 久久久久久久久久黄片| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 国内精品久久久久精免费| 国产高清不卡午夜福利| 亚洲成人中文字幕在线播放| 欧美一区二区亚洲| 99热全是精品| 人人妻人人看人人澡| 亚洲精品亚洲一区二区| 亚洲精品粉嫩美女一区| 麻豆成人午夜福利视频| 亚洲av成人精品一区久久| 九九在线视频观看精品| 日韩欧美一区二区三区在线观看| 国产高清有码在线观看视频| a级毛色黄片| 在线天堂最新版资源| 少妇裸体淫交视频免费看高清| 国产又黄又爽又无遮挡在线| 日韩制服骚丝袜av| 免费在线观看成人毛片| 国产毛片a区久久久久| 丰满人妻一区二区三区视频av| 在线观看66精品国产| 联通29元200g的流量卡| 2021天堂中文幕一二区在线观| 免费观看人在逋| 国产 一区 欧美 日韩| 日本免费一区二区三区高清不卡| 国产精品精品国产色婷婷| 在线看三级毛片| 成人永久免费在线观看视频| 99热网站在线观看| 一进一出抽搐动态| 国内精品一区二区在线观看| 丰满人妻一区二区三区视频av| 最新在线观看一区二区三区| 久久久久国产精品人妻aⅴ院| 日本精品一区二区三区蜜桃| 国产精品一二三区在线看| 黄色配什么色好看| 国产精品人妻久久久影院| 久久久久性生活片| 啦啦啦观看免费观看视频高清| 久久久午夜欧美精品| 97在线视频观看| 国产人妻一区二区三区在| 91狼人影院| 日韩制服骚丝袜av| 久久午夜亚洲精品久久| 小蜜桃在线观看免费完整版高清| 美女大奶头视频| 在现免费观看毛片| 搞女人的毛片| 一进一出好大好爽视频| 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 精品久久久久久成人av| 91久久精品电影网| 六月丁香七月| 午夜视频国产福利| 色在线成人网| 亚洲av免费高清在线观看| 尾随美女入室| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看| 国产精品野战在线观看| 亚洲av.av天堂| 国产亚洲欧美98| 欧美三级亚洲精品| 观看美女的网站| 亚洲人成网站高清观看| 全区人妻精品视频| 超碰av人人做人人爽久久| 国产久久久一区二区三区| 久久热精品热| 色播亚洲综合网| 亚洲av中文字字幕乱码综合| 香蕉av资源在线| 一进一出抽搐动态| 免费av不卡在线播放| 亚洲精品国产成人久久av| av黄色大香蕉| 波多野结衣高清无吗| 18禁裸乳无遮挡免费网站照片| 亚洲欧美成人精品一区二区| 国产一区亚洲一区在线观看| 99久国产av精品| 国产高潮美女av| 午夜福利在线观看吧| 欧美日韩国产亚洲二区| 深夜精品福利| 国产精品久久久久久久久免| 一进一出抽搐gif免费好疼| 岛国在线免费视频观看| 免费一级毛片在线播放高清视频| 国产三级在线视频| 国产精华一区二区三区| 内射极品少妇av片p| 麻豆国产97在线/欧美| 免费看a级黄色片| 日韩av在线大香蕉| a级一级毛片免费在线观看| 男女做爰动态图高潮gif福利片| 99久久精品国产国产毛片| 亚洲人成网站在线观看播放| 日韩欧美 国产精品| 国产一区亚洲一区在线观看| 最近的中文字幕免费完整| 日本熟妇午夜| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩高清专用| 18+在线观看网站| 亚洲精品成人久久久久久| 亚洲最大成人手机在线| 亚洲丝袜综合中文字幕| 看非洲黑人一级黄片| 男人和女人高潮做爰伦理| 亚洲精华国产精华液的使用体验 | 身体一侧抽搐| av视频在线观看入口| 真人做人爱边吃奶动态| 精品国产三级普通话版| 一级黄片播放器| 99久久久亚洲精品蜜臀av| АⅤ资源中文在线天堂| 无遮挡黄片免费观看| 嫩草影院入口| 联通29元200g的流量卡| 男女啪啪激烈高潮av片| 激情 狠狠 欧美| 99久国产av精品| 国产精品一区二区性色av| 免费大片18禁| 国产精品永久免费网站| 日产精品乱码卡一卡2卡三| 六月丁香七月| 久久久欧美国产精品| 少妇丰满av| 99热6这里只有精品| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 一级av片app| 亚洲成人久久性| 免费观看人在逋| 中文字幕久久专区| 国产蜜桃级精品一区二区三区| 国产片特级美女逼逼视频| 国产白丝娇喘喷水9色精品| 卡戴珊不雅视频在线播放| 国语自产精品视频在线第100页| 一区二区三区免费毛片| 亚洲人成网站在线播放欧美日韩| 午夜精品在线福利| 亚洲精品亚洲一区二区| 成人二区视频| 九九久久精品国产亚洲av麻豆| 日韩欧美国产在线观看| 色噜噜av男人的天堂激情| 麻豆精品久久久久久蜜桃| 99热这里只有是精品50| 欧美极品一区二区三区四区| 精品久久久久久久人妻蜜臀av| 五月伊人婷婷丁香| 一本久久中文字幕| 成人美女网站在线观看视频| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 国产高清不卡午夜福利| 91在线精品国自产拍蜜月| 可以在线观看毛片的网站| 在线免费观看的www视频| 99热这里只有是精品在线观看| 欧美最新免费一区二区三区| 亚洲无线在线观看| 午夜久久久久精精品| 国产男靠女视频免费网站| 1024手机看黄色片| av福利片在线观看| 深爱激情五月婷婷| 成人美女网站在线观看视频| 日韩精品有码人妻一区| 熟女人妻精品中文字幕| 精品福利观看| 日韩成人av中文字幕在线观看 | 日韩欧美精品v在线| 在现免费观看毛片| .国产精品久久| 综合色丁香网| 亚洲人成网站在线观看播放| 草草在线视频免费看| 亚洲成a人片在线一区二区| 国产极品精品免费视频能看的| 国产精品女同一区二区软件| 淫妇啪啪啪对白视频| 99久久精品热视频| 秋霞在线观看毛片| 成年版毛片免费区| 色哟哟·www| 欧美一级a爱片免费观看看| www日本黄色视频网| 在线免费观看不下载黄p国产| 看十八女毛片水多多多| 美女xxoo啪啪120秒动态图| 免费观看精品视频网站| 午夜a级毛片| 啦啦啦啦在线视频资源| 中文字幕av成人在线电影| 最近手机中文字幕大全| h日本视频在线播放| 我要搜黄色片| 午夜精品一区二区三区免费看| 亚洲一级一片aⅴ在线观看| 精品人妻一区二区三区麻豆 | 国产精品一区二区三区四区久久| 午夜精品国产一区二区电影 | 亚洲成a人片在线一区二区| 亚洲av二区三区四区| 欧美一区二区国产精品久久精品| 寂寞人妻少妇视频99o| a级毛片a级免费在线| 中文字幕久久专区| 国产高清视频在线观看网站| 尤物成人国产欧美一区二区三区| av在线观看视频网站免费| 午夜激情福利司机影院| 国内少妇人妻偷人精品xxx网站| 91久久精品国产一区二区成人| 国产精品不卡视频一区二区| 免费不卡的大黄色大毛片视频在线观看 | 大型黄色视频在线免费观看| 偷拍熟女少妇极品色| 淫妇啪啪啪对白视频| 热99在线观看视频| 亚洲av不卡在线观看| 国产白丝娇喘喷水9色精品| 国产一区二区三区在线臀色熟女| 国产色婷婷99| 亚洲真实伦在线观看| 精品久久久久久久末码| 亚洲欧美日韩东京热| 在线观看美女被高潮喷水网站| 成人精品一区二区免费| 精品乱码久久久久久99久播| 国产人妻一区二区三区在| 亚洲av五月六月丁香网| 美女xxoo啪啪120秒动态图| 亚洲精品乱码久久久v下载方式| 日本五十路高清| 亚洲五月天丁香| 成人永久免费在线观看视频| 露出奶头的视频| 久久久久久久久中文| 日本 av在线| 成年免费大片在线观看| av中文乱码字幕在线| 舔av片在线| 看片在线看免费视频| eeuss影院久久| 亚洲成a人片在线一区二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲最大成人av| 18+在线观看网站| 天堂√8在线中文| 老司机午夜福利在线观看视频| 亚洲色图av天堂| av视频在线观看入口| 久久久久国产网址| 精品乱码久久久久久99久播| 国产v大片淫在线免费观看| 国产精品福利在线免费观看| 亚洲婷婷狠狠爱综合网| 国产精品乱码一区二三区的特点| 久久精品综合一区二区三区| 亚洲丝袜综合中文字幕| 一进一出好大好爽视频| 成年免费大片在线观看| 中文字幕av成人在线电影| 亚洲av熟女| 波多野结衣巨乳人妻| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久久久久| 欧美日韩精品成人综合77777| 少妇人妻精品综合一区二区 | av天堂中文字幕网| 亚洲五月天丁香| 国产高清不卡午夜福利| 日韩欧美精品免费久久| 亚洲一区二区三区色噜噜| 国国产精品蜜臀av免费| 日韩大尺度精品在线看网址| 成年女人永久免费观看视频| 免费一级毛片在线播放高清视频| 国产一区二区激情短视频| 天堂网av新在线| or卡值多少钱| 黄色欧美视频在线观看| 深夜精品福利| 午夜免费男女啪啪视频观看 | 99热精品在线国产| 国产一区二区亚洲精品在线观看| 天堂av国产一区二区熟女人妻| 亚洲av美国av| 成人综合一区亚洲| 长腿黑丝高跟| 成人精品一区二区免费| a级毛片免费高清观看在线播放| 日韩一区二区视频免费看| 俄罗斯特黄特色一大片| 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情在线99| 人人妻人人澡人人爽人人夜夜 | 久久午夜福利片| 国产老妇女一区| 老司机影院成人| 亚洲精品日韩在线中文字幕 | 久久久国产成人精品二区| av在线蜜桃| 最近在线观看免费完整版| 午夜福利高清视频| 久久国内精品自在自线图片| 精品一区二区三区视频在线观看免费| 午夜福利视频1000在线观看| 亚洲激情五月婷婷啪啪| 麻豆国产av国片精品| 成人av一区二区三区在线看| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区| 亚洲专区国产一区二区| 国产精品女同一区二区软件| 久久久欧美国产精品| 日韩成人av中文字幕在线观看 | 舔av片在线| 97超碰精品成人国产| 欧美性猛交╳xxx乱大交人| 少妇的逼水好多| 国产精品爽爽va在线观看网站| 日韩在线高清观看一区二区三区| 免费看美女性在线毛片视频| 黄色一级大片看看| 欧美zozozo另类| 国产精品久久久久久久电影| 乱系列少妇在线播放| 小说图片视频综合网站| 国产精品一区www在线观看| 成年女人看的毛片在线观看| 91在线精品国自产拍蜜月| 熟妇人妻久久中文字幕3abv| 中国国产av一级| 国产精品亚洲美女久久久| 91在线精品国自产拍蜜月| 亚洲av五月六月丁香网| 午夜精品在线福利| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲网站| 国产乱人视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩三级伦理在线观看| 中文字幕久久专区| 哪里可以看免费的av片| 国产精品电影一区二区三区| 亚洲性夜色夜夜综合| 伦理电影大哥的女人| 日韩亚洲欧美综合| 日本与韩国留学比较| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久久大av| 国产精品野战在线观看| 美女黄网站色视频| 免费一级毛片在线播放高清视频|