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

    激光切割陶瓷加工過程的數(shù)值模擬研究*

    2019-12-09 06:44:56
    陶瓷 2019年11期
    關鍵詞:熱傳導溫度場徑向

    趙 婧

    (榆林市瀚霆化工技術開發(fā)有限公司 陜西 榆林 718100)

    前言

    激光加工指的是激光束作用于物體表面而引起的物體成形或改性的加工過程[1~2]。以激光作為加工能源,在硬脆性陶瓷加工方面的發(fā)展?jié)摿σ岩姸四遊3]:它可以實現(xiàn)無接觸式加工,減少了因接觸應力對陶瓷帶來的損傷;陶瓷對激光具有較高的吸收率(氧化物陶瓷對10.6 μm波長激光的最高吸收率可達80%以上),聚焦的高能激光束作用于陶瓷局部區(qū)域的能量可超過108CJ/cm2,瞬間就可使材料熔化蒸發(fā),實現(xiàn)高效率加工[4];由于聚焦光斑小,產(chǎn)生的熱影響區(qū)小,可以達到精密加工的要求。由于陶瓷是由共價鍵、離子鍵或2者混合化學鍵結合的物質,晶體間化學鍵方向性強,因而具有高硬度和高脆性的本征特性,相對于金屬材料,即使是高精密陶瓷,其顯微結構均勻度亦較差,嚴重降低了材料的抗熱震性,常溫下對剪切應力的變形阻力很大,極易形成裂紋、崩豁甚至于材料碎裂[5]。因此,高效無損傷激光切割陶瓷類高硬脆無機非金屬材料一直是一個亟待解決的問題[7~8]。

    研究激光與物質熱作用過程的理論方法主要有解析法和數(shù)值計算法2種。由于涉及到材料熔融和汽化等相變過程的問題非常復雜,理論解析本解比較困難,甚至于根本無法求解。因此一般采用數(shù)值方法來處理相關問題。隨著計算技術的飛速發(fā)展,人們逐漸提出了一些求解激光熔融問題控制方程的新方法,如將激光作用視為點熱源、作用對象視為半無限大介質、采用一維近似及假定一個簡單的速度場而形成的解析或者半解析法等。而其中的數(shù)值計算方法可已采用有限差分法(FDM)和有限單元法(FEM)等,其中有限單元法包含邊界元法、有限分析法。對激光與物質相互作用的數(shù)值計算就是從特定的物理模型出發(fā),用計算機進行數(shù)值計算或模擬,從而揭示激光與物質相互作用的某些性質和運動規(guī)律[9]。激光與物質的熱作用研究主要是計算激光輻照下物質溫度的變化,以及溫度變化引起的一些熱效應。通常情況下,由于解析或半解析方法必須作一些與實際情況相距甚遠的假設,計算結果可能與實際結果有較大偏差,但可作為進一步研究和純數(shù)值計算的基礎。主要研究內容如下:

    1)描述了Al2O3陶瓷的加熱過程,根據(jù)經(jīng)典熱傳導方程和熱應力方程,建立了不同模式激光輻照下中瞬態(tài)溫度場和熱應力場的物理模型,討論了激光輻射過程中陶瓷溫度場及產(chǎn)生的熱應力場分布。

    2)介紹了激光輻照Al2O3陶瓷的數(shù)學模型,利用軸對稱的熱傳導方程和伽遼金法得到了軸對稱的有限元的熱傳導方程,并介紹了空間軸對稱物體的溫度分布基本理論和有限元解法。

    3)計算得到了Al2O3陶瓷在不同模式的激光加熱時的溫度場以及在施加約束后相應的熱應力場,并進行了比較分析。

    1 激光輻照材料的溫度場理論

    激光照射靶材時,其能量被材料表層所吸收并轉變?yōu)闊?。該熱量通過熱傳導在靶材內擴散,從而形成溫度場,該溫度會導致靶材性質的變化。本節(jié)將主要對激光加熱過程中的熱傳導現(xiàn)象進行分析和計算,這對激光與材料相互作用的物理機理分析具有重要的意義。激光加載條件、物體形狀、初始和邊界條件以及物體的熱物理性質決定了物體溫度場的時空變化。首先給定激光加載條件,通常把被吸收的激光能量或光強作為面熱源(表層吸收)或體熱源(深層吸收)處理。有時邊界上還存在其他種類的熱源或冷卻條件,如表面輻射損失、氣動加熱或冷卻等。不同介質之間界面的熱接觸條件比較復雜,理想情況下界面兩邊的溫度和熱流量應當相等[17~18]。

    由于在物體某處溫度達到熔點時,溫度暫時不再上升,直到熱能達到能夠提供繼續(xù)融化的能量,才能繼續(xù)融化。這個過程是十分復雜的,在此處不予討論。只討論在熔融之前為發(fā)生相變時的升溫過程。在給定熱源、邊界和界面條件、初始溫度分布和物體的熱物理特性時,激光加熱各向同性物體的問題就可以歸結為不定常熱傳導方程的計算[19]。

    1.1 熱導方程與定解條件的理論公式

    三維熱傳導方程可寫為如下通用形式:

    式中:K——熱導率,W·m-1·K-1;

    ρ——材料密度,g/cm3;

    c——材料比熱容,kJ/kg;

    T——溫度,℃;

    t——時間變量,s;

    A(x, y, z, t)——每單位時間、單位體積傳遞熱給固體材料的加熱速率。

    由于材料的熱物理系數(shù)是溫度函數(shù),所以三維熱傳導方程是非線性的,其解非常復雜,很難得到解析解。然而事實上大部分材料的熱物理參數(shù)隨溫度變化并不明顯,故在一定條件下可假定其與溫度無關,在一定的溫度范圍內取其平均值進行計算,這樣方程(2.1)才可能得到解析解。若激光作用下材料是均勻各向同性的,則三維熱傳導方程可簡化為:

    式中:k=K/ρc——材料的熱擴散率,m2/s。

    由于求解熱傳導方程非常復雜,許多學者提出了一些熱模型,在求解熱傳導方程時通常的假定條件是[2,19]:

    1)被加熱材料是各向同性物質。

    2)材料的熱物理參數(shù)與溫度無關或取特定的平均值。

    3)忽略熱傳導中的輻射和對流,只考慮材料表面的熱傳導。

    1.2 激光輻照固體材料的理論模型

    假定激光束垂直入射于物體表面(即x=0),被加熱物體位于右半空間(x≥0),物體表面對激光的反射率為R,吸收系數(shù)為α,x=0處入射激光束的功率密度(光強)是I0(y,z:T),則物體內部的溫度場T可由熱傳導方程進行描述:

    t>0,x>0

    式中:熱物理常數(shù)ρ,c,k與位置及溫度有關; Q是其他體熱源項。上式右部是激光束深層吸收的體熱源,如果αx≈0或R≈1,用表面吸收率A代替(1-R),此體熱源項可改用邊界條件中的面熱源表示:

    光斑處激光束的強度通常可表示為空間分布I0(y,z)和無量綱時間波形B(t)的乘積。典型的B(t)波形有階梯波、矩形波、瞬時作用(δ函數(shù))、三角波和梯形波、指數(shù)衰減波形、Gauss波形和重復頻率波形等。其他面熱源以及常用的溫度、熱流和換熱邊界條件的提法與普通傳熱學問題相同。設t=0為激光開始作用的時刻,T0=(x,y,z:0)是給定的初始溫度分布。如果初始時物體與環(huán)境等溫,則熱傳導方程式計算的是物體的溫升△T=T-T0,考慮常溫即T0=300 ℃的情況。

    考慮被加熱物體厚度大于激光脈沖作用時間內熱量的傳播深度的情況,物體幾何形狀為圓柱形,表面光斑上的激光強度分布對于光束軸線是旋轉對稱的,物體的溫度場對于x軸也一定是旋轉對稱的。激光強度呈高斯分布,I=I0(t)exp(-r2/α2),a是光束的高斯半徑,表面吸收,表面反射率和吸收率分別為R和A,時間波形為連續(xù)函數(shù)B(t),得到溫度場分布為:

    2 激光輻照材料的應力場理論

    軸對稱物體是指其幾何形狀中心軸線,則物體內各處的應力外加載荷或溫度作用,以及約束情況都對稱于某一應變及位移分量也對稱于這一軸線。假設物體是連續(xù)的線性彈性且是各向均勻同性的,在外力或溫度作用下物體的位移和變形是微小的,物體在外加荷載或溫度作用之前處于自然狀態(tài),內部沒有應力存在。

    2.1 平衡微分方程

    在物體內點P外取一微元六面體PABC,六面體由兩個半徑差為dr的圓柱面、兩個夾角為dθ且包含z軸的垂直面和兩個相距為dz的水平面所圍成,作用在圓柱面的徑向正應力為σr。軸向正應力為σz。圓柱面上的剪應力τrz,環(huán)向正應力σθ。

    令體力的徑向分量為R,體力的軸向分量為Z。將微元體六個面的作用力都投影到半徑方向,則

    經(jīng)計算得空間軸對稱物體的平衡微分方程為:

    2.2 幾何方程

    在軸對稱的物體內每一點只能有2個方向的位移,即沿z軸方向的軸向位移及沿半徑:方向的徑向位移,由于對稱關系,物體內的任一點沒有沿圓周方向(環(huán)向)的位移。

    式中:{ε}——應變列向量;

    ur——徑向位移,mm;

    uz——軸向位移,mm。

    2.3 物理方程

    由廣義虎克定律可得:

    將應力分量表示為應變分量的函數(shù)為:

    式中:μ——泊松比;

    E——彈性模量,MPa。

    2.4 熱彈性方程[20]

    對于三維軸對稱模型,在忽略體積力和慣性力的前提下,材料體內的溫度場發(fā)生非均勻變化時產(chǎn)生的熱應變和熱應力可用如下的熱彈性方程來描述:

    式中:ur、uz——(r,z)點的位移在r,z方向上的分量;

    e——體應變,μm;

    β——熱膨脹系數(shù),1/k;

    E——楊氏模量,MPa;

    μ——泊松比。

    而應變、應力之間分別滿足如下的彈性力學關系式:

    e=εr+εz+εθ

    對應材料的簡支狀態(tài),熱彈性方程的自由邊界條件可取為:

    uz|r=a,z=d=0(a為半徑,d為高度)

    對應于靶材夾緊狀態(tài),熱彈性方程的固定邊界條件可取為:

    ur|r=a=uz|r=a=0

    2.5 控制方程、初始條件、邊界條件[21]

    2.5.1 控制方程

    假設陶瓷板為各向同性,其熱傳導方程為[22]:

    式中:c,ρ和k分別是陶瓷的質量熱容、密度和導熱系數(shù);P(x,y,z,t)為體熱源的生熱率。

    2.5.2 初始條件[23]

    t=0, T(t)=T0

    式中:T0——陶瓷的初始溫度,這里視作環(huán)境溫度300 K。

    2.5.3 邊界條件

    在激光加載面。熱傳導向陶瓷內部傳送的熱能,對流換熱向外界逸出的熱能以及表面向外界輻射的熱能,三者之和應該等于其吸收的激光能量,即

    αI(x,y,z,t)=

    式中:Ts——激光加載面的表面溫度,℃;

    h——對流傳熱系數(shù),W/(m2·K);

    ε——表面熱輻射系數(shù),W/(m2·℃);

    σ——波爾茲曼常數(shù),J/K;

    α——陶瓷對激光的吸收率,%;

    I(x,y,z,t)——光能密度,m·W/m2。

    在陶瓷板的前后左右和下表面,僅通過對流的方式與外界進行熱量交換,因此設置為對流的邊界條件[24],即

    式中:?T/?n——溫度沿表面外法線方向的偏導數(shù);

    Tn——玻璃表面的溫度,℃。

    激光輻射可以看作表面加熱源[25],因此P(x,y,z,t)=0。I(x,y,z,t)為面熱源的光能密度,也是激光的功率密度,當激光沿著陶瓷基板的正x軸方向運動時,激光在陶瓷板表面的功率密度可以用函數(shù)表示為:

    I(x,y,z,t)=I(x-vt,y)δ(z)

    式中:P0——激光功率,W;

    r——激光束半徑,mm。

    3 不同模式的激光輻照材料的溫度場分析

    3.1 高斯激光的溫度場分析

    建立了5 mm×2.5 mm的矩陣模型,為減少計算量,不影響模擬的情況下,劃分網(wǎng)格時,為了把更多的計算用于溫度變化明顯的區(qū)域,采用了不等劃分方法,在徑向分為100份,深度方向分為50份。如圖1所示。在計算完結果后可以用ANSYS軟件的3D還原功能將其還原,如圖2所示。

    圖1 柱坐標下軸截面的網(wǎng)格劃分

    圖2 柱坐標3D擴展后網(wǎng)格劃分效果圖

    3.1.1 熱傳導模型

    為了便于分析激光與陶瓷相互作用時的溫度場,做以下假設:

    1)激光束直接作用于材料表面(在材料內沒有新的熱源),這是因為激光的照射時間是短暫的。

    2)視陶瓷樣品表面對激光能量的吸收率為1,且不隨時間變化。事實上,吸收系數(shù)是根據(jù)溫度變化而變化的,但是對相同表面的的條件下,測量一個平均值是可行的。

    圖3 激光照射示意圖

    3)激光照射期間功率等參數(shù)不變。

    4)厚板的維度的簡化,選擇適當?shù)木S數(shù)進行簡化計算,長度,厚度,寬度當?shù)竭_一定尺度時,可以近似認為無窮大。

    5)略熱輻射與熱對流,認為空氣在照射期間是絕熱。

    6)忽略材料的密度隨溫度升高的變化。

    7)材料是各向同性的。

    8)激光束是穩(wěn)態(tài)的。

    在有了以上幾個假定后,就可以用一般熱傳到方程求解溫度場分布情況了。令激光垂直照射在玻璃板表面。

    3.1.2 計算結果溫度場分布圖及分析

    由圖4可以看出,在距離加熱點一定范圍內經(jīng)過激光輻照過的陶瓷,陶瓷表面(徑向)溫度傳播速度要比陶瓷內部(軸向)要快,超過該范圍,溫度趨于一致(即環(huán)境溫度)。中心點最高溫度達到1 867 K。

    圖4功率200 W光斑半徑0.5 mm時間步長5 ms高斯激光照射陶瓷界面溫度分布圖

    圖5徑向距離中心原點分別為0、0.5 mm、1.5 mm、2.0 mm隨時間溫度變化圖

    圖5說明了距離加熱中心越近的溫度越高且溫度上升的速度要比遠心點要快。

    在徑向上,距離加熱中心點越近的區(qū)域,溫度梯度越陡,且梯度隨距離的增大而減小,這是由于距離中心點的溫度上升速度要比遠心點要快,如圖6所示。在軸向上的情況基本差不多,但是相比徑向,軸向的梯度要比徑向要小,且梯度變化相對徑向要緩和,如圖7所示。

    圖6 徑向溫度梯度

    圖7 軸向溫度梯度

    在軸向上,溫度梯度要比徑向要順暢,在軸向上各點的溫度隨時間的增加溫度趨勢如半S形,如圖8所示。

    圖8 軸向不同深度隨時間變化溫度圖

    圖9 不同方向上的溫度變化

    在一定范圍內徑向溫度要比軸向溫度高,超出該范圍,溫度基本一致,即環(huán)境溫度300 K,如圖9所示。

    3.2 帽頂形激光的溫度場分析

    3.2.1 帽頂形激光的載人

    激光功率P=200 W 光斑半徑r=2.5 mm。

    3.2.2 計算結果、溫度分布圖及分析

    圖10 帽頂形激光示意圖

    圖11功率200 W光斑半徑0.5 mm時間5 ms帽頂形激光照射陶瓷界面溫度分布圖

    相比于圖4,帽頂形激光在軸向溫度上要比高斯形激光傳播要快,中心溫度達到2 105 K(假設此時陶瓷沒有熔化)。說明中心點的溫度要比高斯形激光產(chǎn)生的最高溫度要高得多。

    一定范圍內距離加熱中心越近的溫度越高且溫度上升的速度要比遠心點要快。超出該范圍,溫度基本趨于室溫恒定不變,故溫度基本無變化。

    圖12徑向r=0.5、1.0 mm、1.5 mm、2.0 mm、2.5 mm溫度隨時間變化圖

    圖13 徑向溫度梯度

    圖14 軸向溫度梯度

    圖15 軸向不同深度隨時間變化溫度變化情況

    圖16 不同方向上的溫度趨勢

    在徑向上,距離加熱中心點越近的區(qū)域,溫度梯度越陡,且梯度隨距離的增大而減小,也是由于距離中心點的溫度上升速度要比遠心點要快,如圖12所示。在軸向上的情況基本差不多,但是相比徑向,軸向的梯度要比徑向要小,且梯度變化相對徑向要緩和,如圖14所示。

    由圖15可以看出,在一定深度內,各部位的溫度變化近似為半S形,且在深度超過0.6 mm后,陶瓷溫度基本和環(huán)境溫度一致(300 K)。

    3.3 兩種模式的激光產(chǎn)生的溫度場的比較

    比較在相同的功率、時間步長、光束半徑相同的不同模式的激光對同一種模型加熱所產(chǎn)生的溫度。

    3.3.1 徑向上的溫度比較

    在一定的距離內,帽頂形激光的溫度傳播比高斯形激光傳播要快,超過該距離,帽頂形激光要比高斯形激光慢,直到溫度不變。不論是何種模式的激光,吸收激光的能量主要分布在光斑半徑以內。因此隨著時間的增長形成的徑向溫度梯度也越大。

    圖17 不同模式激光徑向溫度場比較

    圖18 不同模式激光在軸向上的溫度比較

    3.3.2 軸向上的溫度比較

    在軸向上,一定范圍內,帽頂形激光要比高斯激光溫度傳播要快。在加熱中心點(即光斑半徑范圍以內),帽頂形激光產(chǎn)生的溫度要比高斯形激光產(chǎn)生的溫度要高。

    4 激光輻照材料的應力場分析

    4.1 模型創(chuàng)建、載入激光

    本次模擬建立了20 mm×5 mm×1 mm的長方體模型,用的是一個三維耦合場固體元素SOLID5(該元素有八個節(jié)點,每個節(jié)點有六個自由度)。并對該模型在長度方向上分了100份,其模型及網(wǎng)格圖如圖19所示。這次設計模擬的是移動激光束對Al2O3陶瓷基板中心軸線的掃描過程其示意圖如圖20所示。

    圖19 應力模擬模型

    圖20 激光掃描陶瓷基本示意圖

    相關參數(shù):掃描速度v=2 m/s;時間步長T=10 ms;激光功率=200 W;激光束半徑r=0.5 mm。載入激光方程。

    4.2 計算應力結果及分析

    圖21 熱應力是Y方向產(chǎn)生的位移

    熱應力總是由于快速的升溫或降溫所產(chǎn)生的,當熱應力超過材料的臨界值的時候,材料就會產(chǎn)生斷裂。正如本次工作中,材料出現(xiàn)了斷裂,如圖21所示。

    在Y方向上產(chǎn)生的應力隨X變化而變化,如圖4.4所示,在一定距離內,陶瓷內部產(chǎn)生的是正應力(即應力方向背向陶瓷基板),超過該距離,內部便產(chǎn)生負應力(應力方向指向陶瓷基板)。過程比較復雜,本次工作不做討論。

    圖22 Y方向上的熱應力

    圖23 X方向的溫度

    圖24 等效應力顯示

    圖25 陶瓷基板表面的溫度分布

    溫度變化最快的部分如圖23所示,也就是產(chǎn)生斷裂的區(qū)域,該部位劇烈的熱對流產(chǎn)生了超過了陶瓷的臨界值,故而產(chǎn)生斷裂。

    在陶瓷基板上,產(chǎn)生了斷裂,最高溫度為1 026 K,要遠低于熔點2 050 K,之所以斷裂,由于陶瓷基板和環(huán)境之間的熱對流,在陶瓷基板的表面,之前的加熱熱量也會沿著陶瓷基板的正方方向傳播,而導致了最高溫度出現(xiàn)在中部,這樣會導致陶瓷基本中部很容易出現(xiàn)劇烈的熱對流現(xiàn)象,使陶瓷基板的部分區(qū)域發(fā)生了迅速的升溫和降溫,從而產(chǎn)生一個極大的應力足夠斷裂陶瓷基板。如圖24和圖25所示。

    5 結論與展望

    在不同模式激光加熱陶瓷的過程中,在相同激光功率、加熱時間步長、相同環(huán)境、模型以及材料的情況下:

    在陶瓷基板距加熱中心點的一定距離內,在基本表面上,帽頂形激光的溫度傳播比高斯形激光傳播要快,超過該距離,帽頂形激光要比高斯形激光慢,直到溫度都接近環(huán)境溫度。

    在陶瓷的軸向方向上,一定范圍內,帽頂形激光要比高斯激光溫度傳播要快,超出該范圍,溫度都趨于環(huán)境溫度。

    在加熱中心點,帽頂形激光產(chǎn)生的溫度要比高斯形激光瞬間產(chǎn)生的溫度要高。

    在激光掃描陶瓷基本的模擬實驗中,激光加熱可以使陶瓷基板沒有達到熔點而出現(xiàn)斷裂,斷裂一般出現(xiàn)在陶瓷基板的中后部,這是由于熱對流的不均勻導致了局部出現(xiàn)了劇烈的升溫或降溫而產(chǎn)生超過臨界值的熱應力導致斷裂出現(xiàn)。

    研究高斯形或者帽頂形激光與Al2O3陶瓷的相互作用過程對激光在工業(yè)、軍事等各個領域的應用都有非常重要的意義,且作用過程尚有許多問題和工作需要探索、分析和解決。在激光切割陶瓷方面,由于陶瓷本征的硬脆特性 ,使得該技術的發(fā)展也存在著許多亟待解決的問題,對陶瓷激光加工的作用機理更需要得到理論上的論證和支持,下面列出一些有待解決的問題和研究內容。

    考慮陶瓷表面熱對流、熱輻射,陶瓷材料的各向異性,熱物參數(shù)隨溫度的變化及材料對激光的吸收系數(shù)隨溫度的變化對陶瓷溫度場和熱應力場的影響。

    陶瓷對激光能量的吸收過程中是個復雜的過程,因此不同模式的激光對陶瓷作用時要綜合考慮界面吸收、雜質吸收和體吸收。

    需要建立更加符合實際的三維模型,根據(jù)陶瓷實際情況施加約束條件得到位移分布和應力場分布,并且可以研究應力場隨時間的變化,得到激光作用于陶瓷的應力場瞬態(tài)變化數(shù)據(jù)。

    研究帽頂形激光加載時可能造成的熔融,氣化反應和應力引起的損傷和斷裂情況。

    可以近一步研究移動激光束對陶瓷作用的溫度場和熱應力場變化,模擬激光焊接陶瓷過程。

    吸收激光的能量主要分布在光斑半徑以內,因此隨著加熱時間增加形成的徑向溫度梯度也越大。

    猜你喜歡
    熱傳導溫度場徑向
    一類三維逆時熱傳導問題的數(shù)值求解
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向對稱正解的存在性
    鋁合金加筋板焊接溫度場和殘余應力數(shù)值模擬
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    一類無窮下級整函數(shù)的Julia集的徑向分布
    基于紋影法的溫度場分布測量方法
    測控技術(2018年4期)2018-11-25 09:47:10
    MJS工法與凍結法結合加固區(qū)溫度場研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    熱傳導方程解的部分Schauder估計
    一類非線性反向熱傳導問題的Fourier正則化方法
    国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 成人18禁在线播放| 国产欧美日韩一区二区三| 国产av一区在线观看免费| 日韩欧美三级三区| 国产99白浆流出| netflix在线观看网站| 亚洲欧美精品综合久久99| 99精品久久久久人妻精品| 亚洲av熟女| 岛国视频午夜一区免费看| 一进一出抽搐gif免费好疼| 精品不卡国产一区二区三区| 亚洲av成人一区二区三| 脱女人内裤的视频| 国产视频内射| 男女视频在线观看网站免费 | av免费在线观看网站| 国产片内射在线| 国产成人欧美在线观看| 欧美日韩一级在线毛片| 亚洲aⅴ乱码一区二区在线播放 | 一级a爱视频在线免费观看| 可以在线观看毛片的网站| 一区二区三区精品91| 1024香蕉在线观看| 丝袜在线中文字幕| 欧美乱色亚洲激情| 91大片在线观看| 亚洲国产精品成人综合色| 丰满的人妻完整版| 久久亚洲精品不卡| av欧美777| 国产一级毛片七仙女欲春2 | 12—13女人毛片做爰片一| 亚洲欧美日韩高清在线视频| 欧美人与性动交α欧美精品济南到| 色综合站精品国产| 最好的美女福利视频网| 日本成人三级电影网站| 麻豆一二三区av精品| 午夜成年电影在线免费观看| 成年女人毛片免费观看观看9| 国产日本99.免费观看| 后天国语完整版免费观看| 久久精品91无色码中文字幕| 国产一区二区三区在线臀色熟女| 亚洲av成人av| 白带黄色成豆腐渣| 高潮久久久久久久久久久不卡| 中文字幕人妻熟女乱码| 日韩欧美一区视频在线观看| 国产精品1区2区在线观看.| 精品久久久久久久久久免费视频| 97超级碰碰碰精品色视频在线观看| 色综合站精品国产| 99久久久亚洲精品蜜臀av| 观看免费一级毛片| 人人妻人人澡人人看| 成在线人永久免费视频| 黄网站色视频无遮挡免费观看| 国内久久婷婷六月综合欲色啪| 制服丝袜大香蕉在线| 欧美在线黄色| 国产精品免费视频内射| 日韩欧美国产在线观看| 久久精品91无色码中文字幕| 国产精品,欧美在线| 99久久久亚洲精品蜜臀av| 美国免费a级毛片| bbb黄色大片| av有码第一页| 国产高清视频在线播放一区| 国产免费男女视频| 黄色丝袜av网址大全| 日韩高清综合在线| 欧美日韩福利视频一区二区| 在线永久观看黄色视频| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全电影3 | 精品无人区乱码1区二区| 久久久久久久午夜电影| 99re在线观看精品视频| 又紧又爽又黄一区二区| 观看免费一级毛片| 成年人黄色毛片网站| 免费看美女性在线毛片视频| 草草在线视频免费看| 欧美日韩精品网址| 男女那种视频在线观看| 最新在线观看一区二区三区| 一夜夜www| 国产精品香港三级国产av潘金莲| 免费无遮挡裸体视频| 久久久久精品国产欧美久久久| 免费在线观看黄色视频的| 国产精品av久久久久免费| 波多野结衣av一区二区av| 久久久久国产精品人妻aⅴ院| 国产精品一区二区免费欧美| 97超级碰碰碰精品色视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美在线黄色| 亚洲成人免费电影在线观看| 91在线观看av| 亚洲第一欧美日韩一区二区三区| 国产精品 国内视频| 中文字幕人妻丝袜一区二区| 在线视频色国产色| 国产精品一区二区精品视频观看| 国产亚洲欧美精品永久| 国产亚洲av嫩草精品影院| 成人手机av| 日本a在线网址| 久久久精品欧美日韩精品| 久久精品国产亚洲av高清一级| 黄频高清免费视频| 亚洲精品久久国产高清桃花| 日韩大码丰满熟妇| 一级作爱视频免费观看| 国产成人精品久久二区二区91| 一个人观看的视频www高清免费观看 | 精品久久蜜臀av无| a在线观看视频网站| 久久国产精品男人的天堂亚洲| 午夜福利高清视频| 美女国产高潮福利片在线看| 精品国产国语对白av| 久久久久精品国产欧美久久久| 黄色片一级片一级黄色片| 别揉我奶头~嗯~啊~动态视频| 在线观看午夜福利视频| 啪啪无遮挡十八禁网站| 亚洲精品中文字幕一二三四区| 免费观看人在逋| 久久久国产精品麻豆| 天堂√8在线中文| 亚洲精品粉嫩美女一区| 日韩欧美一区二区三区在线观看| 在线视频色国产色| 久久人人精品亚洲av| 国产高清视频在线播放一区| 天天一区二区日本电影三级| 人人妻人人看人人澡| 亚洲中文日韩欧美视频| 热re99久久国产66热| 啦啦啦韩国在线观看视频| 美国免费a级毛片| 欧美精品亚洲一区二区| 国产成人欧美| 亚洲一区中文字幕在线| 啦啦啦观看免费观看视频高清| 午夜久久久久精精品| 日本在线视频免费播放| 99热这里只有精品一区 | 亚洲成国产人片在线观看| 国产亚洲欧美精品永久| 精品国产乱子伦一区二区三区| 亚洲国产日韩欧美精品在线观看 | 一边摸一边抽搐一进一小说| www.999成人在线观看| 亚洲精品国产区一区二| 人人妻,人人澡人人爽秒播| 欧美日韩精品网址| 国产97色在线日韩免费| 国产真实乱freesex| 久久久久久久久久黄片| 一区二区日韩欧美中文字幕| 男女下面进入的视频免费午夜 | 久久久精品欧美日韩精品| 久久精品夜夜夜夜夜久久蜜豆 | 给我免费播放毛片高清在线观看| 嫩草影院精品99| 亚洲中文av在线| 久久国产乱子伦精品免费另类| 精品国产亚洲在线| 妹子高潮喷水视频| 亚洲国产精品合色在线| 中文字幕av电影在线播放| 啦啦啦观看免费观看视频高清| 精品久久久久久久久久久久久 | 亚洲av片天天在线观看| 一本精品99久久精品77| www国产在线视频色| www国产在线视频色| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品爽爽va在线观看网站 | 搡老妇女老女人老熟妇| 午夜免费观看网址| 国产主播在线观看一区二区| 97人妻精品一区二区三区麻豆 | 免费电影在线观看免费观看| 无遮挡黄片免费观看| 丝袜在线中文字幕| 亚洲精品国产区一区二| 久久香蕉精品热| 欧美日韩瑟瑟在线播放| 少妇粗大呻吟视频| 色在线成人网| 99国产极品粉嫩在线观看| 欧美三级亚洲精品| 成人永久免费在线观看视频| 91在线观看av| 人人妻人人看人人澡| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美精品综合久久99| 欧美成人免费av一区二区三区| 久久午夜亚洲精品久久| 在线观看免费日韩欧美大片| 男人的好看免费观看在线视频 | 亚洲人成77777在线视频| 欧美性猛交╳xxx乱大交人| 免费av毛片视频| 巨乳人妻的诱惑在线观看| 国产主播在线观看一区二区| 搡老岳熟女国产| 色在线成人网| 中文字幕人成人乱码亚洲影| 人人妻,人人澡人人爽秒播| 韩国精品一区二区三区| 欧美+亚洲+日韩+国产| 亚洲中文字幕一区二区三区有码在线看 | 无限看片的www在线观看| 午夜日韩欧美国产| 香蕉久久夜色| 久久久久免费精品人妻一区二区 | 激情在线观看视频在线高清| 天天一区二区日本电影三级| 精品人妻1区二区| av中文乱码字幕在线| 男人舔女人的私密视频| 一进一出好大好爽视频| 一个人免费在线观看的高清视频| 日韩大码丰满熟妇| 国产成人欧美在线观看| 久久国产精品人妻蜜桃| 日本熟妇午夜| 婷婷精品国产亚洲av在线| 国产免费av片在线观看野外av| 亚洲熟妇熟女久久| 19禁男女啪啪无遮挡网站| 国产精品影院久久| 免费在线观看影片大全网站| 亚洲全国av大片| 草草在线视频免费看| 国产精品久久久av美女十八| 一个人免费在线观看的高清视频| 在线观看一区二区三区| 久久久久久九九精品二区国产 | 成人三级做爰电影| 一本一本综合久久| 香蕉国产在线看| 久久亚洲真实| 1024手机看黄色片| 一卡2卡三卡四卡精品乱码亚洲| 欧美最黄视频在线播放免费| 欧美国产精品va在线观看不卡| 久久婷婷成人综合色麻豆| 一区二区日韩欧美中文字幕| 久久中文字幕人妻熟女| 高清毛片免费观看视频网站| 国产国语露脸激情在线看| 高清在线国产一区| 91在线观看av| 在线观看www视频免费| 久久久久久人人人人人| 国产一区在线观看成人免费| 一进一出好大好爽视频| videosex国产| 99riav亚洲国产免费| www.精华液| 搡老熟女国产l中国老女人| 精品电影一区二区在线| 岛国视频午夜一区免费看| 亚洲熟妇熟女久久| 一边摸一边抽搐一进一小说| 国产真实乱freesex| 精品日产1卡2卡| 亚洲精品av麻豆狂野| 日韩中文字幕欧美一区二区| www.自偷自拍.com| 久久婷婷人人爽人人干人人爱| 91av网站免费观看| 国内精品久久久久久久电影| 国产精品电影一区二区三区| 久久久久久免费高清国产稀缺| 久热这里只有精品99| 欧美日韩乱码在线| 99久久国产精品久久久| 免费在线观看日本一区| 亚洲一区中文字幕在线| 国产男靠女视频免费网站| 日韩欧美国产一区二区入口| 丰满的人妻完整版| 日韩三级视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 视频在线观看一区二区三区| 国产熟女午夜一区二区三区| 欧美国产日韩亚洲一区| 亚洲一区中文字幕在线| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 久久久久久国产a免费观看| 国产精品爽爽va在线观看网站 | 黄频高清免费视频| 日韩大码丰满熟妇| 日韩有码中文字幕| 精品无人区乱码1区二区| 日本在线视频免费播放| 黄色a级毛片大全视频| 9191精品国产免费久久| 777久久人妻少妇嫩草av网站| 一边摸一边做爽爽视频免费| 巨乳人妻的诱惑在线观看| 99热只有精品国产| 少妇裸体淫交视频免费看高清 | 欧美av亚洲av综合av国产av| 久久中文看片网| 亚洲av电影在线进入| 欧美黑人巨大hd| xxxwww97欧美| 在线播放国产精品三级| 成人国产综合亚洲| 每晚都被弄得嗷嗷叫到高潮| 国产又色又爽无遮挡免费看| 国产精品二区激情视频| 黑丝袜美女国产一区| 日本 欧美在线| 热99re8久久精品国产| 精品久久久久久久久久免费视频| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 久久国产乱子伦精品免费另类| 美女高潮喷水抽搐中文字幕| 日韩视频一区二区在线观看| 中国美女看黄片| 狂野欧美激情性xxxx| 欧美久久黑人一区二区| 精品少妇一区二区三区视频日本电影| 波多野结衣高清无吗| 国产精品九九99| 大型av网站在线播放| 精品无人区乱码1区二区| 国产片内射在线| 熟妇人妻久久中文字幕3abv| 久久精品91蜜桃| 精品国产一区二区三区四区第35| 琪琪午夜伦伦电影理论片6080| 99久久无色码亚洲精品果冻| 婷婷丁香在线五月| 看免费av毛片| 一区二区三区高清视频在线| 美女午夜性视频免费| 满18在线观看网站| 黄色片一级片一级黄色片| 窝窝影院91人妻| 亚洲欧美精品综合一区二区三区| 手机成人av网站| 黄网站色视频无遮挡免费观看| 18禁黄网站禁片免费观看直播| 国产精品二区激情视频| 啪啪无遮挡十八禁网站| 十分钟在线观看高清视频www| 淫妇啪啪啪对白视频| 国产亚洲av高清不卡| 日韩成人在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲五月婷婷丁香| 精品第一国产精品| 亚洲va日本ⅴa欧美va伊人久久| 伊人久久大香线蕉亚洲五| 国产私拍福利视频在线观看| 亚洲av熟女| 老司机午夜十八禁免费视频| 黄色毛片三级朝国网站| 成人特级黄色片久久久久久久| 国内揄拍国产精品人妻在线 | 国产色视频综合| 欧美另类亚洲清纯唯美| 91九色精品人成在线观看| 欧美一区二区精品小视频在线| 18禁国产床啪视频网站| 在线观看免费午夜福利视频| 亚洲 欧美一区二区三区| 神马国产精品三级电影在线观看 | 国产亚洲欧美在线一区二区| 国产97色在线日韩免费| 老司机在亚洲福利影院| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 老司机深夜福利视频在线观看| 国产又色又爽无遮挡免费看| 91国产中文字幕| 亚洲精品av麻豆狂野| 亚洲真实伦在线观看| 天天添夜夜摸| 国产精品 欧美亚洲| 夜夜看夜夜爽夜夜摸| 亚洲成人国产一区在线观看| 欧美激情高清一区二区三区| 午夜福利在线观看吧| 最近在线观看免费完整版| 久久久久精品国产欧美久久久| 亚洲五月色婷婷综合| 丰满的人妻完整版| 欧美日韩乱码在线| 精品欧美一区二区三区在线| 国产伦在线观看视频一区| 精华霜和精华液先用哪个| a级毛片在线看网站| 最近最新中文字幕大全电影3 | 一区福利在线观看| 久久国内精品自在自线图片| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 亚洲人成网站高清观看| 久久精品人妻少妇| 一边摸一边抽搐一进一小说| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 精品一区二区三区人妻视频| 此物有八面人人有两片| 综合色丁香网| 久久久a久久爽久久v久久| 老熟妇仑乱视频hdxx| 日本爱情动作片www.在线观看 | 日韩av不卡免费在线播放| 两个人视频免费观看高清| 99热这里只有是精品50| 婷婷精品国产亚洲av在线| 99久久无色码亚洲精品果冻| 国产av麻豆久久久久久久| 搡老岳熟女国产| 美女内射精品一级片tv| 天堂网av新在线| 国产真实伦视频高清在线观看| 久久九九热精品免费| 搞女人的毛片| avwww免费| 国产麻豆成人av免费视频| 看十八女毛片水多多多| 美女黄网站色视频| 日本黄色片子视频| 欧美三级亚洲精品| 超碰av人人做人人爽久久| 国产在线精品亚洲第一网站| 国产av在哪里看| 亚洲精品国产av成人精品 | 国产成人福利小说| 成年免费大片在线观看| h日本视频在线播放| 国产一区二区亚洲精品在线观看| 蜜桃亚洲精品一区二区三区| 午夜福利在线观看免费完整高清在 | 麻豆一二三区av精品| 国产不卡一卡二| 麻豆成人午夜福利视频| 久久久久国产精品人妻aⅴ院| 长腿黑丝高跟| 中文字幕免费在线视频6| 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| av在线亚洲专区| 国产精品国产三级国产av玫瑰| 国产精品一二三区在线看| 亚洲中文字幕日韩| 午夜亚洲福利在线播放| 午夜精品一区二区三区免费看| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 在线a可以看的网站| 亚洲av美国av| 欧美xxxx性猛交bbbb| 亚洲欧美日韩东京热| 亚洲综合色惰| 亚洲av电影不卡..在线观看| 国产老妇女一区| 亚洲国产欧美人成| 国产91av在线免费观看| 亚洲国产高清在线一区二区三| 狠狠狠狠99中文字幕| 99久久成人亚洲精品观看| 免费在线观看影片大全网站| 欧美性感艳星| 别揉我奶头 嗯啊视频| 亚洲精品成人久久久久久| 18禁裸乳无遮挡免费网站照片| 亚洲丝袜综合中文字幕| 男人舔女人下体高潮全视频| 国产精华一区二区三区| 伦理电影大哥的女人| 日本精品一区二区三区蜜桃| 秋霞在线观看毛片| 日本黄大片高清| 我要搜黄色片| 老女人水多毛片| 日日啪夜夜撸| 一区二区三区免费毛片| 亚洲国产精品合色在线| 精品午夜福利在线看| 欧美丝袜亚洲另类| 婷婷精品国产亚洲av| 久久人人爽人人爽人人片va| 欧美高清性xxxxhd video| 国产高潮美女av| 两个人视频免费观看高清| 国产精品野战在线观看| 美女免费视频网站| 一本一本综合久久| 国产免费男女视频| 亚洲自偷自拍三级| 亚洲av电影不卡..在线观看| 两个人的视频大全免费| av视频在线观看入口| 亚洲三级黄色毛片| 3wmmmm亚洲av在线观看| 熟女人妻精品中文字幕| 99久久九九国产精品国产免费| 色哟哟·www| 亚洲不卡免费看| 人人妻,人人澡人人爽秒播| 国产亚洲精品久久久com| 亚洲无线观看免费| 亚洲精品色激情综合| 久久精品国产99精品国产亚洲性色| 夜夜看夜夜爽夜夜摸| 色播亚洲综合网| 麻豆一二三区av精品| 亚洲美女视频黄频| 97超级碰碰碰精品色视频在线观看| 国产熟女欧美一区二区| 色播亚洲综合网| 人人妻人人看人人澡| 丰满人妻一区二区三区视频av| 校园春色视频在线观看| 精品人妻偷拍中文字幕| 久久久久国内视频| 日本熟妇午夜| 亚洲无线在线观看| 啦啦啦观看免费观看视频高清| 婷婷色综合大香蕉| 少妇人妻精品综合一区二区 | 成年女人毛片免费观看观看9| 我的老师免费观看完整版| 免费人成视频x8x8入口观看| 中文资源天堂在线| 波多野结衣高清作品| 国产久久久一区二区三区| 亚洲乱码一区二区免费版| 国产片特级美女逼逼视频| 熟女电影av网| 狠狠狠狠99中文字幕| 在线国产一区二区在线| 国产一区二区在线av高清观看| 日韩欧美国产在线观看| 日本精品一区二区三区蜜桃| 亚洲三级黄色毛片| 国产精品三级大全| 亚洲精华国产精华液的使用体验 | 18禁黄网站禁片免费观看直播| 亚洲人与动物交配视频| ponron亚洲| 男插女下体视频免费在线播放| 国产毛片a区久久久久| 午夜亚洲福利在线播放| 精品久久久久久久人妻蜜臀av| 亚洲av美国av| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩无卡精品| 欧美xxxx性猛交bbbb| 精品一区二区三区人妻视频| 久久久久国内视频| 国产黄片美女视频| av视频在线观看入口| 搡老熟女国产l中国老女人| 亚洲最大成人av| av专区在线播放| 国产高清视频在线观看网站| 日本一本二区三区精品| 菩萨蛮人人尽说江南好唐韦庄 | 麻豆一二三区av精品| 亚洲18禁久久av| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 中文字幕免费在线视频6| 日本色播在线视频| 欧美中文日本在线观看视频| 亚洲精品在线观看二区| 少妇人妻精品综合一区二区 | 国产av一区在线观看免费| 日本 av在线| 99热这里只有精品一区| 精品久久久久久久久亚洲| 成年女人毛片免费观看观看9| 午夜影院日韩av| 在线看三级毛片| 蜜臀久久99精品久久宅男| 一级黄片播放器| 日韩一本色道免费dvd| 久久精品夜色国产| 男女啪啪激烈高潮av片| 久久人人爽人人片av| 亚洲精品成人久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 99久久精品热视频| 99热这里只有是精品50| 国产一区二区在线av高清观看| 国产亚洲精品综合一区在线观看| 亚洲性夜色夜夜综合| 不卡视频在线观看欧美|