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

    基于元胞自動(dòng)機(jī)的鋁合金多坑腐蝕及疲勞壽命研究

    2021-10-18 11:21:06華磊劉雪峰常冬梅
    航空科學(xué)技術(shù) 2021年9期
    關(guān)鍵詞:元胞自動(dòng)機(jī)疲勞腐蝕

    華磊 劉雪峰 常冬梅

    摘要:將陰陽(yáng)極反應(yīng)分離的元胞自動(dòng)機(jī)法引入到鋁合金多坑腐蝕過(guò)程的模擬當(dāng)中。通過(guò)腐蝕模擬得到了較為符合實(shí)際特征的鋁合金材料多坑腐蝕形貌圖。通過(guò)計(jì)算不同腐蝕步數(shù)下的蝕坑形貌參數(shù),發(fā)現(xiàn)蝕坑深度近似于隨腐蝕時(shí)間線性增加,同時(shí)腐蝕失重也隨時(shí)間增加,但是增加的速度逐漸增加。將腐蝕形貌導(dǎo)入ANSYS軟件中進(jìn)行力學(xué)分析,給出了鋁合金材料腐蝕后的應(yīng)力分布,并且由此計(jì)算了腐蝕后的剩余疲勞壽命。發(fā)現(xiàn)隨著腐蝕時(shí)間的增加,剩余疲勞壽命大體上呈逐漸減少的趨勢(shì)。但是,由于蝕坑底端附近局部幾何形貌造成應(yīng)力集中的影響,這一趨勢(shì)并不一定總是單調(diào)的。

    關(guān)鍵詞:鋁合金;腐蝕;元胞自動(dòng)機(jī);疲勞

    中圖分類號(hào):TG171文獻(xiàn)標(biāo)識(shí)碼:ADOI:10.19452/j.issn1007-5453.2021.09.009

    基金項(xiàng)目:航空科學(xué)基金(201902067004);中央高?;究蒲袠I(yè)務(wù)費(fèi)項(xiàng)目(3122018C008)

    鋁合金材料具有重量(質(zhì)量)輕、比強(qiáng)度大等優(yōu)點(diǎn),因此在航空航天領(lǐng)域得到廣泛的應(yīng)用。但是鋁合金材料在服役過(guò)程中容易發(fā)生腐蝕,從而影響材料的可靠性。因此,對(duì)鋁合金腐蝕問(wèn)題的研究受到廣泛的關(guān)注。過(guò)去對(duì)于鋁合金腐蝕問(wèn)題的研究以試驗(yàn)方法為主,但是試驗(yàn)方法往往費(fèi)時(shí)費(fèi)力,成本較高,因此近些年在腐蝕數(shù)值模擬方面的研究越來(lái)越多。多數(shù)腐蝕數(shù)值模擬方法是針對(duì)一系列有關(guān)材料和環(huán)境參數(shù)方程的求解,因此為確定性方法[1]。確定性方法可以較準(zhǔn)確地模擬出材料的宏觀腐蝕進(jìn)程,但是,確定性方法沒(méi)有考慮到材料和環(huán)境在介觀或微觀尺度上的不均勻性,因此不適合模擬由此造成的材料較小尺度上的不規(guī)則腐蝕。而這種較小尺度的不規(guī)則腐蝕往往對(duì)材料的疲勞、斷裂等力學(xué)性質(zhì)有較大的影響。因此,在研究鋁合金的不規(guī)則腐蝕以及其造成的應(yīng)力集中等力學(xué)問(wèn)題時(shí),有必要引入非確定性腐蝕模擬方法。

    元胞自動(dòng)機(jī)方法是一種典型的非確定性算法,最早被用來(lái)模擬生物的生長(zhǎng)問(wèn)題,之后被應(yīng)用于植被生長(zhǎng)過(guò)程[2]、交通流[3]和火災(zāi)蔓延過(guò)程[4]等很多方面的模擬。近年來(lái),該方法作為一種有效的非確定性方法被引入到金屬材料的腐蝕模擬中,取得了較好的效果。2004年,Vautrin等[5]在元胞自動(dòng)機(jī)模型中允許陰極和陽(yáng)極反應(yīng)在不同位置同時(shí)發(fā)生反應(yīng),很好地模擬了真實(shí)的電化學(xué)腐蝕過(guò)程,并且由此驗(yàn)證了蝕坑內(nèi)的氫離子濃度指數(shù)(pH)不均衡現(xiàn)象。Aarao等[6]在2006年對(duì)這一模型的相關(guān)參數(shù)進(jìn)行了分析。之后出現(xiàn)了一系列基于這一模型的研究[7-9]。2016年,Pérez-Brokate等[10]將該模型引入到三維情況下,得到了更符合實(shí)際的模擬結(jié)果。2013年,Wang[11]和Han[12]將有限元與元胞自動(dòng)機(jī)腐蝕模型結(jié)合,研究了應(yīng)力作用下的腐蝕模擬問(wèn)題。隨后將研究引申到了雙蝕坑的問(wèn)題上。但是,到目前為止,基于元胞自動(dòng)機(jī)模擬得到的腐蝕形貌對(duì)材料力學(xué)性質(zhì)的研究仍然很少見(jiàn)到,而在這方面的研究有利于分析介觀到宏觀尺度腐蝕對(duì)材料力學(xué)性質(zhì)的影響。因此,本文將基于陰陽(yáng)極分離的元胞自動(dòng)機(jī)腐蝕模型模擬鋁合金多坑腐蝕形貌,并且在此基礎(chǔ)上研究鋁合金疲勞壽命隨腐蝕進(jìn)程的變化。

    1基于元胞自動(dòng)機(jī)的多坑腐蝕形貌

    1.1元胞自動(dòng)機(jī)模型

    在酸性堿性元胞的擴(kuò)散方面,假設(shè)A、B元胞在每個(gè)反應(yīng)步中會(huì)發(fā)生Ndiff次布朗運(yùn)動(dòng),當(dāng)其運(yùn)動(dòng)到E元胞的位置,則與之交換位置,而當(dāng)A元胞運(yùn)動(dòng)到B元胞的位置或相反時(shí),A、B元胞中和,兩個(gè)元胞均變?yōu)镋元胞。每個(gè)反應(yīng)步中,每個(gè)C元胞都有一定概率發(fā)生失效,該元胞和其下面的相鄰M元胞變?yōu)镋元胞。

    以上反應(yīng)的概率見(jiàn)表1。其中,NA-B為鄰居中A、B元胞個(gè)數(shù)差,NE為鄰居元胞中液態(tài)元胞(A、B、E)的個(gè)數(shù)。在每一步腐蝕模擬過(guò)程中,每次隨機(jī)選取一對(duì)可以發(fā)生反應(yīng)的R元胞和S元胞,按照式(6)~式(9)發(fā)生反應(yīng),反復(fù)選取,直到?jīng)]有可供反應(yīng)的元胞為止。然后進(jìn)行Ndiff步A、B元胞的擴(kuò)散模擬,每次隨機(jī)選取一個(gè)A或B元胞進(jìn)行擴(kuò)散,直到所有A、B元胞選取結(jié)束算作一個(gè)擴(kuò)散步。Ndiff步擴(kuò)散步結(jié)束后,算作一步腐蝕模擬結(jié)束,進(jìn)行下一步腐蝕模擬,直到結(jié)束所有模擬步。

    1.2模擬結(jié)果

    本文采用摩爾鄰居,元胞個(gè)數(shù)為1024×128個(gè),設(shè)PB為1/100000,PSSE為0.2,PRE為0,PPE為1/8,PPE為1/800,即為PPE的1/100。因此,在中性條件下,再鈍化產(chǎn)物的水解概率為酸性條件下的1/100。每個(gè)反應(yīng)步中的擴(kuò)散步數(shù)Ndiff為200。在MATLAB中編寫(xiě)程序,由此得到500、1000和1500反應(yīng)步后的腐蝕形貌如圖1(a)~圖1(c)所示,腐蝕形貌與實(shí)際腐蝕形貌特征相符。由圖可以明顯地看到,隨著計(jì)算步數(shù)的增加,蝕坑逐漸增大,而且有新的蝕坑逐漸產(chǎn)生。蝕坑邊界具有明顯的不規(guī)則性,接近于真實(shí)的蝕坑,這與確定性方法得到的結(jié)果是不同的。在第500步時(shí),橫坐標(biāo)0~100區(qū)域內(nèi)有三個(gè)小蝕坑。由圖1(b)可見(jiàn),到第1000步時(shí),這三個(gè)小蝕坑已經(jīng)合并成為一個(gè)大蝕坑,而且相比于其他蝕坑,這個(gè)合并的大蝕坑更不規(guī)則。在試驗(yàn)中經(jīng)常會(huì)看到類似的蝕坑合并現(xiàn)象。

    針對(duì)模型的隨機(jī)性,每個(gè)步數(shù)進(jìn)行5次計(jì)算取平均值,得到平均腐蝕失重率隨著步數(shù)的變化曲線如圖2所示。其中,腐蝕失重為A、B、E元胞個(gè)數(shù)除以總元胞個(gè)數(shù)。腐蝕失重隨時(shí)間逐漸增加,而且隨著時(shí)間的增加,腐蝕失重增加的速度有所增加。

    除了材料的腐蝕失重率以外,蝕坑的深度也是受關(guān)注的焦點(diǎn)。尤其是所有蝕坑中的最大深度,對(duì)材料的損傷情況有較大的影響。本文中,將R和P元胞作為蝕坑表面,由此將蝕坑表面的最低點(diǎn)與上表面的距離(元胞個(gè)數(shù))作為最大蝕坑深度。與前面相同,計(jì)算不同腐蝕步數(shù)情況下5次模擬的最大蝕坑平均深度如圖3所示??梢钥吹?,最大蝕坑平均深度近似隨著腐蝕步數(shù)線性增加。

    2疲勞壽命分析

    在本文中,利用有限元軟件ANSYS進(jìn)行鋁合金腐蝕后的疲勞壽命分析。首先需要將腐蝕形貌圖導(dǎo)入到ANSYS中。而ANSYS無(wú)法自動(dòng)識(shí)別腐蝕形貌,需要用其他軟件將蝕坑邊界信息提取出來(lái)再導(dǎo)入到ANSYS中建立模型。

    首先,識(shí)別蝕坑的邊界。如圖1可以看到,蝕坑的邊界實(shí)際上是由點(diǎn)描繪的圖形。這里采用軟件R2V提取材料的邊界,然后將結(jié)果導(dǎo)入AutoCAD,以圖1(b)中的1000步的腐蝕形貌圖為例,得到邊界圖如圖4所示。

    在AutoCAD中調(diào)整邊界尺寸,假設(shè)每個(gè)元胞的邊長(zhǎng)為2μm,則整個(gè)模型的尺寸為2048μm×256μm。將圖形保存后導(dǎo)入ANSYS進(jìn)行受力分析。選取的鋁合金材料楊氏模量為7.17×104MPa,泊松比0.3,左側(cè)固定,右側(cè)施加10N/μm的拉應(yīng)力,進(jìn)行計(jì)算得到x方向的應(yīng)力云圖如圖5所示。可以看到,在蝕坑底端附近出現(xiàn)了明顯的應(yīng)力集中現(xiàn)象,通常在較大的蝕坑底端的應(yīng)力集中更明顯。

    應(yīng)力集中會(huì)大大影響材料的疲勞壽命。以計(jì)算得到的應(yīng)力為基礎(chǔ),可以進(jìn)一步計(jì)算鋁合金材料腐蝕后的剩余疲勞壽命。基于ANSYS疲勞分析模塊,以模型的最大Mises應(yīng)力作為疲勞分析的最大應(yīng)力,以0作為最小應(yīng)力計(jì)算疲勞壽命,S—N曲線數(shù)值表見(jiàn)表2。其中參數(shù)由試驗(yàn)數(shù)據(jù)進(jìn)行擬合得到,由于本文中主要針對(duì)疲勞壽命的趨勢(shì)問(wèn)題進(jìn)行定性研究,因此具體的參數(shù)選擇僅取為典型值即可。

    在本例中,由此得到其對(duì)應(yīng)的疲勞壽命為2.43×105次。

    按照以上步驟,重新進(jìn)行一次腐蝕進(jìn)程模擬,分別將步數(shù)為500、600、800、1000、1200、1400、1500時(shí)的腐蝕形貌導(dǎo)入ANSYS中,得到材料疲勞壽命隨腐蝕步數(shù)的變化曲線如圖6所示。隨著腐蝕步數(shù)的增加,疲勞壽命大體上是逐漸減少的。但是,這一趨勢(shì)并不是線性的,甚至當(dāng)腐蝕步數(shù)為600步時(shí)的疲勞壽命還高于500步時(shí)。這是因?yàn)樵谶@種條件下的疲勞壽命變化主要與應(yīng)力集中有關(guān),而應(yīng)力集中除取決于蝕坑深度外,也與局部的腐蝕形貌有很大關(guān)系。在有些時(shí)候盡管蝕坑深度已經(jīng)較大,但是蝕坑底端的表面形貌比較光滑,因此局部應(yīng)力反而有所減小。例如,在本例中,腐蝕步數(shù)為500和600步時(shí),最大應(yīng)力分別為32.14MPa和30.03MPa,由于在500步時(shí)應(yīng)力集中造成的最大應(yīng)力更大,因此在500步時(shí)計(jì)算得到的疲勞壽命反而更少。但是由于蝕坑深度是單調(diào)增加的,因此疲勞壽命的整體趨勢(shì)隨腐蝕步數(shù)的增加而逐漸減少。

    3結(jié)論

    本文基于元胞自動(dòng)機(jī)法模擬了鋁合金材料的多坑電化學(xué)腐蝕過(guò)程,在算法中允許陰陽(yáng)極反應(yīng)發(fā)生在不同位置。通過(guò)模擬,成功得到了較符合實(shí)際特征的腐蝕形貌圖。在此基礎(chǔ)上,將腐蝕形貌進(jìn)行處理后導(dǎo)入ANSYS軟件中進(jìn)行力學(xué)分析,得到了鋁合金材料腐蝕后的應(yīng)力分布,并由此分析了鋁合金材料腐蝕后的疲勞壽命,得到了以下結(jié)論:

    (1)運(yùn)用陰陽(yáng)極反應(yīng)分離的元胞自動(dòng)機(jī)算法可以成功模擬鋁合金的腐蝕進(jìn)程,并且可以觀測(cè)到實(shí)際應(yīng)用中常見(jiàn)到的蝕坑合并現(xiàn)象。

    (2)在本算法模擬的條件下,鋁合金的腐蝕失重隨時(shí)間逐漸增加,而且增速逐漸加大,蝕坑的最大深度接近于隨時(shí)間線性增加。

    (3)鋁合金腐蝕后的剩余疲勞壽命大體上隨腐蝕時(shí)間減少,但是并不完全單調(diào)遞減。這是因?yàn)槠趬勖c應(yīng)力集中現(xiàn)象有關(guān),而應(yīng)力集中除取決于蝕坑深度以外,也受到蝕坑底部局部幾何形貌的影響,而這隨著腐蝕過(guò)程而劇烈變化。

    參考文獻(xiàn)

    [1]Walton J C,Cragnolino G,Kalandros S K. A numerical model of crevice corrosion for passive and active metals[J]. Corrosion Science,1996,38(1):1-18.

    [2]Dunn A G,Majer J D. Simulating weed propagation via hierarchical,patch-based cellular automata[C]//Computational Science-ICCS 2007:7th International Conference Proceedings,2007:762-769.

    [3]彭麟,譚惠麗,孔令江,等.開(kāi)放性邊界條件下雙車道元胞自動(dòng)機(jī)交通流模型耦合效應(yīng)研究[J].物理學(xué)報(bào),2003, 52(12): 3007-3013. Peng Lin, Tan Huili, Kong Lingjiang, et al. A study of coupling effect in cellular automata model of traffic flow for two-lane with open boundary conditions[J]. Acta Physica Sinica,2003, 52(12): 3007-3013.(in Chinese)

    [4]Encinas A H,Encinas L H,White S H,et al. Simulation of forest fire fronts using cellular automata[J]. Advances in Engineering Software,2007,38:372-378.

    [5]Vautrin-Ul C,Chausse A,Stafiej J,et al. Simulations of corrosion processes with spontaneous separation of cathodic and anodic reaction zones[J]. Polish Journal of Chemical Technology,2004,78(9):1795-1810.

    [6]Aarao F D A,Stafiej J,Badiali J P. Scaling theory in a model of corrosion and passivation[J]. The Journal of Physical Chemistry B,2006,110(35):17554-17562.

    [7]Vautrin-Ul C,Taleb A,Stafiej J,et al. Mesoscopic modelling of corrosion phenomena:coupling between electro chemical and mechanical processes,analysis of the deviation from the Faraday law[J]. Electrochimica Acta,2007,52(17):5368-5376.

    [8]Di Caprio D,Stafiej J. Cellular automata approach to corrosion and passivity phenomena[J]. Pure and Applied Chemistry,2010,85(1):62-67.

    [9]Bartosik?,Di Caprio D,Stafiej J. Cellular automata approach to corrosion and passivity phenomena[J]. Pure and Applied Chemistry,2013,85(1):247-256.

    [10]Pérez-Brokate C F,Di Caprio D,F(xiàn)éron D. Three dimensional discrete stochastic model of occluded corrosion cell[J]. Corrosion Science,2016,111:230-241.

    [11]Wang H,Han E. Simulation of metastable corrosion pit development under mechanical stress[J]. Corrosion Science,2013,90:128-134.

    [12]Wang H,Han E. Computational simulation of corrosion pit interactions under mechanochemical effects using a cellular automaton/finite element model[J]. Corrosion Science,2016,103:305-311.

    [13]Liu X F,Hua L,Chang D M. Simulation of pit interactions of multi-pit corrosion under an anticorrosive coating with a threedimensional cellular automata model[J]. Modelling and Simula tion in Materials Science and Engineering,2021,29(6):065018.

    Study on Multi-Pits Corrosion of Aluminum Alloy and Its Fatigue Life Based on Cellular Automata

    Hua Lei1,Liu Xuefeng1,Chang Dongmei2

    1. Civil Aviation University of China,Tianjin 300300,China

    2. Tianjin Key Laboratory of High Speed Cutting and Precision Machining,Tianjin University of Technology and Education,Tianjin 300222,China

    Abstract: In this paper, the cellular automata method of reaction separation between anode and cathode is introduced into the simulation of multi-pit corrosion process of aluminum alloy. Through corrosion simulation, the multipit corrosion morphology of aluminum alloy was obtained. It is found that the pit depth increases linearly with the corrosion time approximately, and the corrosion weight loss also increases with the corrosion time, but the increasing speed increases with time gradually. The corrosion morphology is imported into ANSYS software for mechanical analysis, the stress distribution of aluminum alloy after corrosion is obtained, and the residual fatigue life after corrosion is calculated. It is found that the residual fatigue life decreases with corrosion time. However, this trend is not always monotonous due to the effect of local geometry near the bottom of the pit on the stress concentration.

    Key Words: aluminum alloys; corrosion; cellular automata; fatigue

    猜你喜歡
    元胞自動(dòng)機(jī)疲勞腐蝕
    基于元胞自動(dòng)機(jī)的高速公路隧道群交通流建模與仿真
    淺析游泳訓(xùn)練的疲勞特點(diǎn)及恢復(fù)方法
    日本人開(kāi)始測(cè)量“疲勞”
    鎂合金的腐蝕特性及防護(hù)探討
    科技資訊(2016年19期)2016-11-15 08:38:47
    如何有效消除青少年中長(zhǎng)跑訓(xùn)練后的疲勞
    基于元胞自動(dòng)機(jī)的行人和車輛疏散機(jī)理研究
    論述天然氣埋地鋼管腐蝕與防護(hù)
    鍍鋅產(chǎn)品表面腐蝕現(xiàn)象研究
    基于元胞自動(dòng)機(jī)模擬滬金高速道路車流中的應(yīng)用
    考試周刊(2016年62期)2016-08-15 07:16:20
    日本人開(kāi)始測(cè)量“疲勞”
    男女国产视频网站| 麻豆av在线久日| 亚洲人成77777在线视频| 午夜免费观看性视频| 午夜福利一区二区在线看| 国产色婷婷99| 精品一品国产午夜福利视频| 免费在线观看完整版高清| 久久久欧美国产精品| 亚洲男人天堂网一区| 亚洲视频免费观看视频| 大香蕉久久网| 日韩一本色道免费dvd| 欧美精品人与动牲交sv欧美| 免费高清在线观看日韩| 亚洲人成77777在线视频| 在线观看免费日韩欧美大片| 在线天堂最新版资源| 亚洲一区二区三区欧美精品| a级毛片在线看网站| 日本wwww免费看| 街头女战士在线观看网站| 日韩制服骚丝袜av| 69精品国产乱码久久久| 午夜福利网站1000一区二区三区| 久久精品国产亚洲av涩爱| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片黄色毛片免费观看视频| 亚洲图色成人| 久久精品久久久久久久性| 国产在线一区二区三区精| 免费黄网站久久成人精品| 久久久久久久国产电影| 精品亚洲成a人片在线观看| 久久久久久人妻| 男人添女人高潮全过程视频| 看免费成人av毛片| 精品福利永久在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲第一区二区三区不卡| 国产黄色免费在线视频| 搡女人真爽免费视频火全软件| 日韩伦理黄色片| 亚洲,欧美,日韩| 蜜桃在线观看..| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻丝袜一区二区 | 99久久综合免费| 在现免费观看毛片| av在线老鸭窝| 边亲边吃奶的免费视频| 男女国产视频网站| 青春草国产在线视频| 国产精品国产av在线观看| 天天躁夜夜躁狠狠久久av| 欧美日韩亚洲国产一区二区在线观看 | 国产成人欧美| 久久鲁丝午夜福利片| 国产精品成人在线| 啦啦啦在线免费观看视频4| 老汉色∧v一级毛片| 国产伦理片在线播放av一区| 制服诱惑二区| 少妇猛男粗大的猛烈进出视频| 三上悠亚av全集在线观看| 久久鲁丝午夜福利片| 五月伊人婷婷丁香| 啦啦啦视频在线资源免费观看| 香蕉精品网在线| 在线天堂中文资源库| 观看av在线不卡| 亚洲精品久久午夜乱码| 成年av动漫网址| 国产精品.久久久| 国产欧美日韩综合在线一区二区| 咕卡用的链子| 天天影视国产精品| 欧美精品一区二区大全| 亚洲精品国产色婷婷电影| 一二三四在线观看免费中文在| 美女午夜性视频免费| 国产成人精品久久二区二区91 | 伦理电影免费视频| 亚洲成人一二三区av| 亚洲成av片中文字幕在线观看 | 超碰成人久久| 天天影视国产精品| 亚洲av中文av极速乱| 欧美另类一区| 一区二区三区激情视频| 久久久久精品人妻al黑| 日本色播在线视频| av在线播放精品| 午夜精品国产一区二区电影| 欧美日韩综合久久久久久| 一级毛片 在线播放| 亚洲欧美清纯卡通| 成年女人在线观看亚洲视频| 大香蕉久久成人网| 日日爽夜夜爽网站| 波多野结衣av一区二区av| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 亚洲五月色婷婷综合| 亚洲欧美成人精品一区二区| 26uuu在线亚洲综合色| 免费观看性生交大片5| 午夜精品国产一区二区电影| 天堂中文最新版在线下载| 亚洲av福利一区| 国产 精品1| 热99久久久久精品小说推荐| 在线观看免费日韩欧美大片| 色播在线永久视频| 日韩精品免费视频一区二区三区| 超色免费av| 色视频在线一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲经典国产精华液单| 欧美精品国产亚洲| 一边摸一边做爽爽视频免费| 久久精品国产亚洲av涩爱| 搡女人真爽免费视频火全软件| av电影中文网址| 天天影视国产精品| 香蕉国产在线看| 久久 成人 亚洲| 日日啪夜夜爽| 国产精品 国内视频| 精品第一国产精品| 亚洲精品一区蜜桃| 午夜福利网站1000一区二区三区| 成人影院久久| 亚洲国产成人一精品久久久| 夫妻性生交免费视频一级片| 精品午夜福利在线看| 侵犯人妻中文字幕一二三四区| 欧美激情 高清一区二区三区| 最黄视频免费看| 亚洲国产欧美网| 欧美av亚洲av综合av国产av | 久久久国产一区二区| 国产黄色视频一区二区在线观看| 久久鲁丝午夜福利片| 女人精品久久久久毛片| 免费黄色在线免费观看| 岛国毛片在线播放| 中文精品一卡2卡3卡4更新| 最近的中文字幕免费完整| 国产探花极品一区二区| 十分钟在线观看高清视频www| 久久久久精品性色| 欧美激情极品国产一区二区三区| 久久久久国产一级毛片高清牌| 亚洲欧美中文字幕日韩二区| 久久鲁丝午夜福利片| av免费在线看不卡| a级毛片黄视频| 久久久精品免费免费高清| 在线天堂最新版资源| 亚洲激情五月婷婷啪啪| 色网站视频免费| 最近中文字幕2019免费版| 亚洲av国产av综合av卡| 精品一区二区三区四区五区乱码 | 韩国精品一区二区三区| 免费少妇av软件| 亚洲国产日韩一区二区| 免费观看在线日韩| 美女主播在线视频| 久久久国产一区二区| av国产久精品久网站免费入址| 久久狼人影院| 人体艺术视频欧美日本| 在线观看国产h片| 久久久久久久大尺度免费视频| 国产av精品麻豆| 日本猛色少妇xxxxx猛交久久| 精品国产超薄肉色丝袜足j| 国产男女超爽视频在线观看| 国产精品女同一区二区软件| 国产精品麻豆人妻色哟哟久久| 国产精品 国内视频| 亚洲 欧美一区二区三区| 国产野战对白在线观看| av在线观看视频网站免费| 91精品伊人久久大香线蕉| 国产淫语在线视频| 日韩精品有码人妻一区| 97精品久久久久久久久久精品| av有码第一页| 国产精品欧美亚洲77777| 老汉色∧v一级毛片| 久久人妻熟女aⅴ| 久久午夜综合久久蜜桃| tube8黄色片| 街头女战士在线观看网站| 人妻人人澡人人爽人人| 赤兔流量卡办理| 欧美bdsm另类| 97精品久久久久久久久久精品| 亚洲中文av在线| 精品午夜福利在线看| 成人二区视频| 欧美日韩一区二区视频在线观看视频在线| 日本av手机在线免费观看| 老女人水多毛片| 这个男人来自地球电影免费观看 | 日韩制服丝袜自拍偷拍| 免费高清在线观看视频在线观看| 午夜福利影视在线免费观看| 91精品三级在线观看| av又黄又爽大尺度在线免费看| 国产亚洲一区二区精品| 永久网站在线| 久久韩国三级中文字幕| 亚洲国产精品999| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品亚洲一区二区| 国产av码专区亚洲av| 成人影院久久| 国产精品久久久久久久久免| 亚洲成色77777| 一级毛片电影观看| 国产亚洲一区二区精品| 性高湖久久久久久久久免费观看| 国产视频首页在线观看| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| 9热在线视频观看99| 国产欧美亚洲国产| 国产片内射在线| 国产日韩欧美亚洲二区| 日本色播在线视频| 人妻系列 视频| 最近中文字幕2019免费版| 两个人看的免费小视频| 国产av一区二区精品久久| 视频区图区小说| 国产黄色视频一区二区在线观看| 午夜福利在线观看免费完整高清在| 天堂8中文在线网| 久久久久久久亚洲中文字幕| 哪个播放器可以免费观看大片| 曰老女人黄片| 欧美最新免费一区二区三区| av国产精品久久久久影院| 男人添女人高潮全过程视频| 99久久中文字幕三级久久日本| 青春草国产在线视频| 制服丝袜香蕉在线| 午夜日韩欧美国产| 亚洲,欧美,日韩| 熟女少妇亚洲综合色aaa.| 丝瓜视频免费看黄片| 麻豆av在线久日| 婷婷色综合www| 91精品三级在线观看| 国产精品免费大片| 国产成人精品一,二区| 日韩伦理黄色片| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 不卡av一区二区三区| 黄色一级大片看看| 91久久精品国产一区二区三区| av在线app专区| 两性夫妻黄色片| 国产老妇伦熟女老妇高清| 免费大片黄手机在线观看| 亚洲精品国产色婷婷电影| 久久久久久久亚洲中文字幕| 国产97色在线日韩免费| 亚洲精品一区蜜桃| 99香蕉大伊视频| 秋霞在线观看毛片| 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 精品酒店卫生间| 久久影院123| 精品一区二区三卡| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线| 成人手机av| 男人爽女人下面视频在线观看| videosex国产| 亚洲色图综合在线观看| 久久久久精品人妻al黑| 日韩成人av中文字幕在线观看| 成人午夜精彩视频在线观看| 777米奇影视久久| 亚洲精品一区蜜桃| 国产麻豆69| 久热久热在线精品观看| 免费人妻精品一区二区三区视频| 久久久久精品性色| 亚洲精品aⅴ在线观看| 国语对白做爰xxxⅹ性视频网站| 不卡视频在线观看欧美| 国产精品女同一区二区软件| 尾随美女入室| 久久ye,这里只有精品| 亚洲第一av免费看| 中国三级夫妇交换| 欧美成人午夜免费资源| 人人妻人人澡人人看| 国产精品久久久久久精品古装| 美女大奶头黄色视频| 一区二区三区乱码不卡18| 日韩,欧美,国产一区二区三区| 伦理电影免费视频| 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 久久久久国产一级毛片高清牌| 国产乱人偷精品视频| 久久精品国产自在天天线| 国产乱来视频区| 肉色欧美久久久久久久蜜桃| 久久久精品94久久精品| 男人操女人黄网站| 欧美亚洲 丝袜 人妻 在线| 性色avwww在线观看| 成人午夜精彩视频在线观看| 26uuu在线亚洲综合色| 色视频在线一区二区三区| 韩国av在线不卡| 成年人免费黄色播放视频| 如日韩欧美国产精品一区二区三区| 夫妻午夜视频| 18禁动态无遮挡网站| 国产福利在线免费观看视频| 国产高清国产精品国产三级| 国产免费一区二区三区四区乱码| 久久狼人影院| 超色免费av| 国产亚洲精品第一综合不卡| 亚洲精品美女久久av网站| 精品一区二区免费观看| 亚洲欧美中文字幕日韩二区| 国产精品二区激情视频| 日本wwww免费看| 黄色视频在线播放观看不卡| 美女高潮到喷水免费观看| 又粗又硬又长又爽又黄的视频| 久热久热在线精品观看| 男女免费视频国产| 国产人伦9x9x在线观看 | 男女啪啪激烈高潮av片| 欧美日韩国产mv在线观看视频| 午夜免费男女啪啪视频观看| 久久久久精品久久久久真实原创| 欧美日韩精品网址| 久久久久国产网址| 一区二区三区四区激情视频| 欧美激情极品国产一区二区三区| 毛片一级片免费看久久久久| 亚洲av电影在线进入| 嫩草影院入口| 久久久久久久久久久免费av| 欧美日韩国产mv在线观看视频| 国产精品一区二区在线观看99| 啦啦啦在线免费观看视频4| 久久精品国产亚洲av涩爱| 综合色丁香网| 免费大片黄手机在线观看| av国产精品久久久久影院| 国产在线视频一区二区| 综合色丁香网| 美女高潮到喷水免费观看| 高清不卡的av网站| 国产免费现黄频在线看| 午夜福利影视在线免费观看| 久久99一区二区三区| 亚洲四区av| 大话2 男鬼变身卡| 最近2019中文字幕mv第一页| freevideosex欧美| 激情视频va一区二区三区| 成人手机av| 妹子高潮喷水视频| 婷婷色综合大香蕉| 妹子高潮喷水视频| 又大又黄又爽视频免费| 亚洲av在线观看美女高潮| 精品国产一区二区久久| 极品少妇高潮喷水抽搐| 亚洲美女视频黄频| 男女午夜视频在线观看| 色吧在线观看| 伊人久久国产一区二区| 精品国产乱码久久久久久小说| 国产精品国产三级国产专区5o| 国产黄色免费在线视频| 999精品在线视频| 成人黄色视频免费在线看| 丝袜在线中文字幕| 精品一区二区三区四区五区乱码 | 在线观看免费高清a一片| 亚洲中文av在线| 9191精品国产免费久久| 满18在线观看网站| 亚洲欧美精品综合一区二区三区 | av不卡在线播放| 香蕉丝袜av| 一本大道久久a久久精品| 寂寞人妻少妇视频99o| www日本在线高清视频| 亚洲一码二码三码区别大吗| 亚洲精品成人av观看孕妇| 亚洲精品久久久久久婷婷小说| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 国产在线一区二区三区精| 日韩大片免费观看网站| 人妻 亚洲 视频| 伦理电影免费视频| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 香蕉丝袜av| 亚洲一区中文字幕在线| 在线观看免费日韩欧美大片| 成年av动漫网址| av免费在线看不卡| 久热这里只有精品99| 精品第一国产精品| 成人亚洲精品一区在线观看| 日韩欧美精品免费久久| 久久99一区二区三区| 国产 精品1| 91国产中文字幕| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 波多野结衣av一区二区av| 亚洲内射少妇av| 熟女av电影| 国产av国产精品国产| 色网站视频免费| 亚洲精品日本国产第一区| a级毛片在线看网站| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 精品亚洲成国产av| 免费不卡的大黄色大毛片视频在线观看| 日日啪夜夜爽| 国产精品三级大全| 十八禁网站网址无遮挡| 精品国产国语对白av| 最近最新中文字幕大全免费视频 | 国产乱人偷精品视频| 亚洲美女视频黄频| 久久国产精品男人的天堂亚洲| 免费av中文字幕在线| 欧美黄色片欧美黄色片| 久久午夜福利片| 亚洲一码二码三码区别大吗| 亚洲精品在线美女| 国产一区二区在线观看av| 日本欧美视频一区| 久久午夜综合久久蜜桃| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 欧美激情高清一区二区三区 | 亚洲第一av免费看| 国产日韩一区二区三区精品不卡| 色94色欧美一区二区| 五月天丁香电影| 超色免费av| 伊人久久大香线蕉亚洲五| 日韩中文字幕视频在线看片| 18+在线观看网站| 久久99热这里只频精品6学生| 国产精品.久久久| 99久国产av精品国产电影| 中文字幕人妻丝袜制服| 亚洲三区欧美一区| 日本vs欧美在线观看视频| 少妇精品久久久久久久| 国产人伦9x9x在线观看 | 五月天丁香电影| 伊人亚洲综合成人网| 国产片特级美女逼逼视频| 亚洲精品久久成人aⅴ小说| 成年美女黄网站色视频大全免费| 综合色丁香网| 国产精品不卡视频一区二区| 夫妻午夜视频| 免费少妇av软件| 少妇的丰满在线观看| 伊人久久大香线蕉亚洲五| 久久精品国产自在天天线| 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 曰老女人黄片| 欧美xxⅹ黑人| 高清欧美精品videossex| videosex国产| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜爱| 不卡视频在线观看欧美| 日本av手机在线免费观看| 国产探花极品一区二区| 观看美女的网站| 不卡视频在线观看欧美| 综合色丁香网| 亚洲三级黄色毛片| 国产精品人妻久久久影院| 老女人水多毛片| 久久这里只有精品19| av线在线观看网站| 极品人妻少妇av视频| av在线app专区| 色吧在线观看| 久久这里有精品视频免费| 国精品久久久久久国模美| 免费在线观看视频国产中文字幕亚洲 | 精品午夜福利在线看| 国产精品二区激情视频| 超碰成人久久| 波多野结衣一区麻豆| 国产日韩欧美在线精品| 国产综合精华液| 亚洲经典国产精华液单| 亚洲欧美日韩另类电影网站| 色婷婷av一区二区三区视频| 久久毛片免费看一区二区三区| 日本爱情动作片www.在线观看| 永久网站在线| 亚洲国产欧美日韩在线播放| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩一区二区三区精品不卡| 99精国产麻豆久久婷婷| 秋霞在线观看毛片| 国产一区二区三区av在线| 久久精品国产a三级三级三级| 亚洲国产欧美在线一区| 精品少妇黑人巨大在线播放| 99久久精品国产国产毛片| 午夜日本视频在线| 99久久人妻综合| 黑人猛操日本美女一级片| 国产成人精品在线电影| 亚洲欧美精品综合一区二区三区 | 黄色 视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲,欧美,日韩| 免费日韩欧美在线观看| 哪个播放器可以免费观看大片| 欧美97在线视频| 日日摸夜夜添夜夜爱| a级毛片在线看网站| 乱人伦中国视频| 男女国产视频网站| 国产精品av久久久久免费| 久久人人爽人人片av| 国产极品天堂在线| tube8黄色片| 免费观看性生交大片5| 成人毛片a级毛片在线播放| 18禁动态无遮挡网站| 国产亚洲午夜精品一区二区久久| 久久鲁丝午夜福利片| 97人妻天天添夜夜摸| 91国产中文字幕| 国产又爽黄色视频| 国产成人免费无遮挡视频| 1024视频免费在线观看| 国产亚洲av片在线观看秒播厂| 欧美日韩视频高清一区二区三区二| 男女免费视频国产| 久久久久久人妻| 欧美精品av麻豆av| 免费观看av网站的网址| 成人二区视频| 午夜久久久在线观看| 国产白丝娇喘喷水9色精品| 狂野欧美激情性bbbbbb| 色吧在线观看| 久久鲁丝午夜福利片| 国产精品免费视频内射| 免费播放大片免费观看视频在线观看| 亚洲视频免费观看视频| 欧美 亚洲 国产 日韩一| 最黄视频免费看| 国产又色又爽无遮挡免| 超碰成人久久| 成人国产麻豆网| 尾随美女入室| 国产精品熟女久久久久浪| 日韩精品有码人妻一区| a级片在线免费高清观看视频| 国产在视频线精品| 男人操女人黄网站| 80岁老熟妇乱子伦牲交| 又黄又粗又硬又大视频| 免费在线观看黄色视频的| 日韩在线高清观看一区二区三区| 亚洲四区av| 搡女人真爽免费视频火全软件| 黄片小视频在线播放| 国产老妇伦熟女老妇高清| 精品久久久久久电影网| 精品少妇黑人巨大在线播放| 日韩 亚洲 欧美在线| 一本—道久久a久久精品蜜桃钙片| 亚洲国产日韩一区二区| 国产片特级美女逼逼视频| 国产一区二区在线观看av| 91精品三级在线观看| 狂野欧美激情性bbbbbb| 日韩三级伦理在线观看| 成人免费观看视频高清| 亚洲三区欧美一区| xxx大片免费视频|