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

    Thermal stress damage mechanism in single-crystal germanium caused by 1080 nm laser irradiation

    2023-11-02 08:11:46YinChuanSha沙銀川ZeWenLi李澤文ZhiChaoJia賈志超BingHan韓冰andXiaoWuNi倪曉武
    Chinese Physics B 2023年10期
    關鍵詞:韓冰銀川

    Yin-Chuan Sha(沙銀川), Ze-Wen Li(李澤文)?, Zhi-Chao Jia(賈志超),Bing Han(韓冰), and Xiao-Wu Ni(倪曉武)

    1School of Science,Nanjing University of Science and Technology,Nanjing 210094,China

    2School of Materials Engineering,Jinling Institute of Technology,Nanjing 211169,China

    3School of Electronic Engineering and Optical Engineering,Nanjing University of Science and Technology,Nanjing 210094,China

    Keywords: thermal stress,single-crystal germanium,fracture,damage mechanism

    1.Introduction

    The study of laser interaction processes with semiconductor materials and the corresponding damage mechanisms has been a popular research area and is receiving increasing attention.As an important substrate and infrared window material,single-crystal germanium has important applications in photodetectors, infrared systems, solar cells and other fields.[1-3]With the development of laser processing technology,the use of lasers for punching, slotting and scribing germanium has also become a reality.[4-8]Compared with millisecond lasers,short-pulse lasers such as femtosecond and picosecond lasers have superior processing quality.However,the laser energy is easily shielded due to the plasma produced by the high power density, thus reducing the efficiency of the process.The use of millisecond lasers will solve these issues.However,singlecrystal germanium is subject to thermal stress damage during millisecond laser ablation due to its inherent brittleness,which reduces the yield rate of the process and the reliability of the subsequent product made of germanium.Therefore, study of the thermal stress damage process and mechanism of laser ablation of single-crystal germanium is of great importance for laser processing of germanium materials and products.

    Laser thermal stress damage to semiconductor silicon,such as slip and fracture, has been reported extensively and studied in various ways.Choiet al.used an optical deflection method to probe the ablation process of single-crystal silicon under continuous laser ablation at 1070 nm in real time and showed that the optical deflection method can effectively detect the onset of slip.[9]In addition,they used numerical simulations to obtain the moment of slip onset by comparing thermal stress and yield strength and verified their numerical results experimentally.[10]Liet al.investigated anisotropic thermal stresses in single-crystal silicon using experimental and numerical simulations.Their results showed that thermal slip occurs because the thermal stresses exceed the critical yield stress.[11]Experimental and numerical simulations of thermal stresses in millisecond laser-induced silicon by Xiet al.found that both melt damage and cracks(thermal stress damage)appeared on the sample surface when the laser power density reached 259 J·cm-2.[12]Jiaet al.conducted an experimental and numerical simulation study of laser ablation of singlecrystal silicon by using a high-speed charge-coupled device(CCD)to receive an amplified sample surface morphology in real time, demonstrating that thermal stress damage occurred both during and after the laser action,and that the second stage of thermal stress damage caused by the resolidification of the molten liquid was more severe.[13,14]

    The aforementioned studies mainly focused on the thermal stress damage to single-crystal silicon by a laser, while little research has been reported on the thermal stress damage to single-crystal germanium.In this paper, a high-speed CCD will be used to record the process of laser ablation of single-crystal germanium at a wavelength of 1080 nm in real time.The damage morphology of the ablated sample will be observed with an optical microscope.A three-dimensional numerical model based on Fourier’s heat conduction equation,Hooke’s law and the Alexander-Hasson(AH)equation will be developed to calculate the temperature, stress and dislocation fields during and after the laser action.The dislocation distribution on the sample surface after laser ablation at different power densities is quantitatively investigated and the dislocation density at specific locations is analyzed over time.

    2.Experimental setup

    The experimental setup for recording the stress damage due to 1080 nm laser ablation of single-crystal germanium in real time is shown in Fig.1.The experimental setup can be divided into an ablation light path section and a detection light path section.The ablation light path uses a continuous fiber laser with a wavelength of 1080 nm, which is focused on the sample surface by a convex lens with a focal length of 100 mm and a spot diameter of 2 mm, with a fixed ablation time of 10 ms.The laser beam has a Gaussian power distribution with M2 less than 1.30.The power density of the probe laser is negligible compared with that of the ablation laser.On the other side of the ablation path a high-speed CCD is used to receive the probe beam and an infrared filter is placed in front of the CCD to avoid interference from the irradiating laser.The highspeed CCD has a recording time of 1 s and records 500 frames per second.A DG535 time delay trigger is used to connect the ablation and detection light paths.In order to remove the oxygen from the sample surface and to minimize the effect of the black smoke from the oxidation of the sample on the realtime recording process,the laser irradiation is accompanied by a side blast of low-velocity nitrogen gas at~10 m·s-1in the tangential direction of the sample.

    Fig.1.Diagram of the experimental setup.

    The sample used in this experiment was a square singlecrystal germanium sheet with a side length of 10 mm and a thickness of 175 μm, without any doping.Five separate experiments were performed for each different laser energy density.The ablated samples were examined with an optical microscope to observe the damage morphology of the ablated samples.

    3.Experimental results

    Overall,the experimental results can be divided into three types based on the different damage morphologies.A typical comparison of the damage morphology of the sample surface after laser ablation at different energy densities is shown in Fig.2.

    Fig.2.Typical damage profile of a sample surface after laser ablation.

    In type 1 damage morphology the laser energy density was 0-12 J·cm-2, at which point no visible signs of damage occurred on the surface of the sample after laser ablation(Fig.2(a)).In type 2 damage morphology the laser energy density was 21-201 J·cm-2.In this case, the surface of the sample melted and a black molten pool area appeared on the surface of the sample after cooling at the end of laser ablation(Fig.2(b)).The high-speed CCD real-time recording process did not record the formation of cracks,fragments flying out or other phenomena.Type 3 damage morphology was seen for laser energy densities≥212 J·cm-2.In this case,high-speed CCD recording of the damage process shows that the laser ablation causes melting of the sample surface,and a few dozen to several hundred milliseconds after the end of ablation cracks suddenly appear in the sample surface,accompanied by fragments flying out and other phenomena.The damage to the sample surface after ablation is shown in Fig.2(c), with the black center being the melting pool area after cooling.Several irregular notches can be seen at the edge of the melting pool(red ovals in Fig.2(c)),which are the traces left by the flying fragments,also accompanied by the formation of cracks.

    Figure 3 is a high-speed CCD real-time recording of a set of damage morphologies with time, corresponding to a laser energy density of 212 J·cm-2,that is,type 3 damage morphology.The laser starts to emit light at 0 ms, Fig.3(a) corresponds to 0 ms, when the laser has not yet started to ablate the sample and only the detecting laser spot can be seen on the sample surface.Figure 3(b) corresponds to 12 ms, when the ablation laser has just finished its action.It is important to note that although the laser ablation time is 10 ms,the 10 ms image is overexposed due to the flames and flashes generated during the ablation process,so we have taken the 12 ms image as the surface profile of the sample at the end of the lasing.A black circular area can be seen inside the probe spot,which is the area of the melting pool left by laser ablation of the sample surface.At this point, cracks have not yet appeared and the faint black area at the edge of the melting pool represents some of the black smoke produced by the ablation process.Figure 3(c)corresponds to a time of 182 ms.The ablation area has not changed except for the black smoke floating above the surface of the sample.In Fig.3(d), taken at 184 ms, a black cracked area suddenly appears at the edge of the melting pool on the surface of the sample,whereas there was no such black cracked area in the previous frame.This indicates that the crack appeared very quickly and that their formation process took less than 2 ms.Figure 3(e)corresponds to 186 ms,when,in addition to the black cracked area that appeared before, a small black dot appears on the edge of the detection spot.Finally, Fig.3(f) represents the situation at 1000 ms when the surface of the sample has stabilized and there is no significant change in the melting pool or the crack,apart from some black smoke.

    The damage process described above shows that the melting pool appears during ablation of single-crystal germanium by a 1080 nm fiber laser, while the cracks appear during the cooling of the sample,a phenomenon similar to the stress damage caused by laser ablation of single-crystal silicon.[13]

    In order to further study the law of crack formation, the moments of crack emergence corresponding to different energy densities were analyzed and the results are shown in Fig.4.It can be seen from the figure that the moments when cracks appear vary greatly for the same energy density;there is no obvious pattern for the time of crack appearance for different energy densities(this conclusion is subject to further study for different ablation times and different sample sizes).This result shows that the moments at which cracks appear have a certain randomness,which is not related to the energy density of the laser but only to the nature of the material,such as defects and other factors.The defects themselves are random,which also leads to a randomness in the moment of crack appearance.It should be noted, however, that the randomness of the moment of crack occurrence does not mean that the appearance of cracks is completely random; from the results of the whole experiment,it is only when the laser energy density reaches or exceeds the damage threshold that cracks start to appear.Furthermore, as can be seen from Figs.2 and 3, the formation of cracks is accompanied by a process of fragments flying out.The flying out of the fragments should be part of a violent stress and energy release process,which is when the cracks first start to form, and then the cracks expand in the direction of the concentration of stress, so the location of the flying out of the fragments can be considered as the point of crack initiation.

    Fig.3.Real-time damage process on samples subjected to laser ablation at 212 J·cm-2.The laser ablation time was 10 ms.

    4.Numerical simulation

    In order to further investigate the physical process and damage mechanism of laser ablation of single-crystal germanium, a numerical simulation model as shown in Fig.5 was developed.The Gaussian incident laser irradiates vertically onto the (001)surface of the sample.The sample is a square single-crystal germanium sheet with a side length of 10 mm and a thickness of 175μm.TheX,YandZaxes correspond to the crystal directions[100],[010]and[001],respectively.

    Fig.5.Numerical simulation diagram of 1080 nm laser ablation of single-crystal germanium.

    The temperature fieldT(x,y,z,t) of the sample subjected to laser ablation can be derived by solving the three-dimensional Fourier heat conduction equation as follows:[15,16]

    whereρ,candkdenote the density,specific heat capacity and thermal conductivity of the material, respectively, andQrepresents the laser heat source.

    The absorption depth of single-crystal germanium for a 1080 nm laser at room temperature is very shallow and negligible compared with the spot diameter.Therefore,to simplify the calculations, this model simplifies the heat sourceQto a surface heat source,denoted as

    whereI0is the power density of the incident laser at the spot center,Rrepresents the reflectivity of the single-crystal germanium,andf(x,y)andg(t)denote the spatial and temporal distributions of the laser,respectively.For a Gaussian laser,I0andf(x,y)can be written as

    wherePdenotes the laser power at the sample surface andais the spot radius at the sample surface.

    When the temperature of the sample rises to the melting or boiling point,the material undergoes two phase change processes.In order to simulate the phase change processes, this model will use the equivalent specific heat capacity method to correct for the specific heat capacityc.The equivalent specific heat capacitycecan be expressed as[15]

    whereHmandHvrepresent the latent heat of melting and latent heat of vaporization of the material, respectively, andδmandδvare pulse smoothing functions that can be expressed as[15]

    whereTmandTvare the melting and boiling points of the material,respectively and ΔTm,vis the pulse width,set at 10 K.

    When a sample is thermally expanded, thermal stresses and thermal strains are developed within the material due to internal and external constraints.For single-crystal germanium, the stress{τ}T=[σx σy σz τyz τxz τxy] and strain{ε}T=[εxx εyy εzz γyz γxz γxy] relationship can be expressed by the generalized Hooke’s law[17-19]

    whereαTis the coefficient of thermal expansion,is the plastic strain tensor component, ΔTrepresents the temperature rise andE,ν,Gare Young’s modulus,Poisson’s ratio and the shear modulus, respectively.When the material temperature rises above the melting point, Young’s modulus is set to a value of 100 (~1011at room temperature) to describe the stress relief process caused by the phase change.[19]

    In addition to elastic deformation, plastic deformation also occurs during temperature rise of a material.[20-26]Microscopically, plastic deformation occurs when microdefects within the material proliferate, a phenomenon that is irreversible, and we can study this process by modelling the change in dislocation density.Therefore, the AH equation is introduced in this model to calculate the dislocation field during the temperature rise process.The dislocation multiplication process can be expressed as[20-25]

    whereNrepresents the dislocation density of the material and its initial value is set to 100 m-2,irepresents the current slip system,jrepresents the other slip systems,KandK?represent the material constants,Arepresents the hardening factor.νrepresents the slip velocity of the dislocation and is the effective stress,respectively expressed as

    whereQ0is the excitation energy of the dislocation,kBis the Boltzmann constant,ν0andmare material constants,τeffis the effective stress,τis the shear stress,bis the Burgers vector,andβandβ?are material constants which characterize the strength of the interaction.

    Single-crystal germanium is an anisotropic crystal with 12 slip systems.However,due to its symmetry and in order to simplify calculations,only the S1(111)[ˉ101]and S3(111)[01ˉ1]slip systems are considered in this model.The specific values of the relevant parameters and a more detailed modelling procedure can be found in Ref.[26].

    5.Numerical simulation results and discussion

    When the laser energy density is low, physical phenomena such as slipping and melting on the surface of semiconductor materials have been reported in a large body of literature from both experimental and numerical simulations.These results are broadly similar and the mechanisms are the same,and these will not be repeated in this paper.Here we first look at the overall trend in the dislocation field.Figure 6 shows a plot (from numerical simulations) of the laser energy density versus the maximum dislocation density in the sample at the end of the laser.As shown in this figure, the maximum dislocation density increases with increasing laser energy density.The experimental results show that when the laser energy density is≥212 J·cm-2,cracks will appear on the surface of the sample after the laser action.The increase in dislocation density will lead to more defects in the sample,which will reduce the fracture strength of the material and facilitate the formation of crack initiation points,thus providing the necessary conditions for fracture to occur.Although there is no order of magnitude difference in dislocation density between the simulated results at 212 J·cm-2and 201 J·cm-2, this increasing trend suggests that crack formation requires the accumulation of a certain amount of laser energy before it can occur.This process is similar to the generation of laser damage thresholds,so we propose dislocation fracture thresholds(Dth)to describe this cumulative process.Figure 6 corresponds to a dislocation fracture threshold of 2.4×1014m-2, which indicates that the maximum dislocation density in the sample needs to increase to 2.4×1014m-2in order for the sample to fracture under the parameters in this paper.It should be noted that the dislocation fracture threshold is dependent on the laser parameters and material parameters as well as temperature.Furthermore,there is no relationship between the moment of fracture and this threshold.These findings are in agreement with the experimental results.

    Fig.6.Variation of the maximum dislocation density in the sample with laser energy density at the end of the laser action.

    The fracture process of the material at the end of the laser action is analyzed below,using a laser power density of 245 J·cm-2as an example.

    Based on the results of previous studies it is clear that the slip of single-crystal germanium is dominated by the S1 slip system.[24]The dislocation distribution of the S1 slip system on the sample surface at the end of the 245 J·cm-2laser action is shown in Fig.7.It is clear from this figure that there are two symmetric regions about the center of the spot in the[ˉ120]direction and the[1ˉ20]direction,which have a much higher dislocation density than the other regions.Considering the symmetry,there are a total of eight such regions on the(001)plane.During cooling, the rapid proliferation of dislocation density in these regions leads to a rapid decrease in fracture strength and also provides fracture initiation points.Ignoring the effect of defects in the sample itself,fracture is most likely to occur first at one of these locations where the dislocation density is at its maximum,at which time the stress is released.The crack then extends outwards, and given that the fracture strength is lowest in the annular direction,the crack is most likely to extend along the annulus.Furthermore,although there are eight possible initiation points,the stresses at the other possible initiation points will not reach the fracture strength as the stresses are released after fracture,so that far fewer than eight initiation points are generally observed.This fracture process explains the fracture morphology of the samples in Figs.2 and 3 very well.

    Fig.7.Dislocation density distribution of the S1 slip system on the(001) face at the end of the laser action, with an action laser energy density of 245 J·cm-2.

    Figure 8 gives the relationship between stress and dislocation density as a function of time at the maximum dislocation density in the sample.As can be seen from the figure,the whole process can be divided into three stages.In the first stage, the laser irradiation stage, the stress increases rapidly at this point and so does the dislocation density,but the absolute value of dislocation density is of the order of 102 m-2.The dislocation proliferation at this point does not have much effect on the fracture strength of the sample,so fracture is unlikely to occur.In the second stage,around the end of the laser irradiation,the center of the laser spot is still transferring heat to the point,so the stress at the point is still increasing,and so is the dislocation density.As in the first stage,the dislocation density is not sufficient to cause fracture.Finally,in the third stage, the overall heat transfer outwards from the spot begins to decrease and the dislocation density begins to increase at an exponential rate.When the dislocation density increases to a certain level, a cracking point is formed, completing the entire fracture process as shown in Fig.3.From Fig.8, the fracture process should occur after 130 ms, while the probability of fracture decreases with time due to the decreasing stress values.Thus,the most likely time for fracture to occur is in the range of approximately 130-200 ms.Considering the influence of realistic sample defects and other physical processes, these conclusions can be considered to be in general agreement with the experimental results.

    Fig.8.Stress and dislocation density plotted against time, with points taken at the maximum dislocation density.

    6.Conclusion

    This paper investigates the effect of progressively increasing laser energy density at 1080 nm on the stress damage process in single-crystal germanium samples subjected to laser irradiation.The experimental results show that fracture does not occur until the laser energy density is greater than 212 J·cm-2;the fracture process captured by the high-speed CCD indicates that fracture occurs during the cooling process;the moment of fracture is not specifically related to the energy density.Numerical simulations show that: the sample dislocation fracture threshold is 2.4×1014m-2; the fracture process is influenced by both dislocation density and stress; there are eight possible fracture initiation points, but stress will be released after fracture occurs at one of the initiation points and the rest will have an impact on the crack path.The research in this paper will provide some theoretical references for laser processing of semiconductor materials and will deepen the understanding of the process of thermal stress damage generated in semiconductor materials under laser irradiation.

    猜你喜歡
    韓冰銀川
    數(shù)字銀川
    中國火炬(2023年8期)2023-08-16 14:44:40
    生命大于天
    房地產(chǎn)估價風險評估與防范
    基于正運動控制器的龍門機器人在鋼鐵行業(yè)的運用
    寧夏科協(xié)“健康養(yǎng)生,食品安全”大講堂在銀川順利舉辦
    寧夏科協(xié)“健康養(yǎng)生,食品安全”大講堂在銀川順利舉辦
    共享單車的“銀川模式”
    Preliminary Study about Narrative Art of the Micro Film
    銀川威力
    風能(2016年12期)2016-02-25 08:45:50
    從銀川跳到硅谷
    久久天堂一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 日韩中文字幕欧美一区二区| 99精品欧美一区二区三区四区| 亚洲av熟女| 99久久精品国产亚洲精品| 欧美最黄视频在线播放免费| 三级国产精品欧美在线观看 | 两人在一起打扑克的视频| 在线观看日韩欧美| 18美女黄网站色大片免费观看| 俄罗斯特黄特色一大片| 国产精品精品国产色婷婷| 制服丝袜大香蕉在线| 亚洲精品中文字幕一二三四区| 69av精品久久久久久| 久久精品91无色码中文字幕| 久久天躁狠狠躁夜夜2o2o| 欧美乱妇无乱码| 波多野结衣巨乳人妻| 欧美日韩国产亚洲二区| 国产精品久久久久久亚洲av鲁大| 成年女人永久免费观看视频| 成人av在线播放网站| 亚洲av成人精品一区久久| 黄色日韩在线| 亚洲,欧美精品.| 久久久久国内视频| 亚洲欧洲精品一区二区精品久久久| 色老头精品视频在线观看| 日韩精品青青久久久久久| 国产精品一及| 俺也久久电影网| 成人18禁在线播放| 精品国产亚洲在线| 一二三四社区在线视频社区8| 日韩高清综合在线| 日本一本二区三区精品| 免费在线观看影片大全网站| 国产成人福利小说| 黄色片一级片一级黄色片| 国产精品一区二区免费欧美| 国产极品精品免费视频能看的| 国产精品久久久久久精品电影| 欧美成人性av电影在线观看| 日本在线视频免费播放| 国产99白浆流出| 美女黄网站色视频| 午夜精品久久久久久毛片777| 日韩欧美国产在线观看| 亚洲精品一卡2卡三卡4卡5卡| 少妇丰满av| 99国产精品一区二区蜜桃av| av黄色大香蕉| 三级男女做爰猛烈吃奶摸视频| 桃色一区二区三区在线观看| 一a级毛片在线观看| av在线天堂中文字幕| 久久九九热精品免费| 一本精品99久久精品77| 身体一侧抽搐| 成年女人永久免费观看视频| 首页视频小说图片口味搜索| 亚洲va日本ⅴa欧美va伊人久久| 18禁黄网站禁片免费观看直播| 最近视频中文字幕2019在线8| 亚洲欧美激情综合另类| 国产成+人综合+亚洲专区| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区| 夜夜夜夜夜久久久久| 国产亚洲精品综合一区在线观看| 亚洲真实伦在线观看| 亚洲狠狠婷婷综合久久图片| xxxwww97欧美| 亚洲欧美日韩东京热| av中文乱码字幕在线| 国产一区二区三区在线臀色熟女| 成人亚洲精品av一区二区| 欧美绝顶高潮抽搐喷水| 在线视频色国产色| 久久久久国内视频| 脱女人内裤的视频| 亚洲精品一区av在线观看| 亚洲性夜色夜夜综合| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 中文字幕人成人乱码亚洲影| 禁无遮挡网站| 免费一级毛片在线播放高清视频| 女人高潮潮喷娇喘18禁视频| 日韩av在线大香蕉| 99久久成人亚洲精品观看| 在线观看免费视频日本深夜| 99在线人妻在线中文字幕| 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频| 精品国产亚洲在线| or卡值多少钱| 精品久久久久久成人av| 女同久久另类99精品国产91| 波多野结衣高清作品| 久9热在线精品视频| 99热这里只有是精品50| 日韩免费av在线播放| 日韩三级视频一区二区三区| aaaaa片日本免费| 国产精品99久久久久久久久| 97人妻精品一区二区三区麻豆| 51午夜福利影视在线观看| 色综合欧美亚洲国产小说| 欧美激情久久久久久爽电影| 99热只有精品国产| 真人一进一出gif抽搐免费| 国产探花在线观看一区二区| 久久久久久九九精品二区国产| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产99精品国产亚洲性色| 757午夜福利合集在线观看| 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 看免费av毛片| 亚洲成人精品中文字幕电影| 中国美女看黄片| 两个人看的免费小视频| av天堂在线播放| 一个人免费在线观看电影 | 久久久久久人人人人人| 亚洲精品粉嫩美女一区| 久久精品国产99精品国产亚洲性色| 欧美日韩瑟瑟在线播放| 91麻豆精品激情在线观看国产| 毛片女人毛片| 色综合婷婷激情| 99久久久亚洲精品蜜臀av| 国产亚洲欧美在线一区二区| 免费看日本二区| 熟妇人妻久久中文字幕3abv| 色哟哟哟哟哟哟| 女人高潮潮喷娇喘18禁视频| 久久精品夜夜夜夜夜久久蜜豆| 极品教师在线免费播放| 国产av一区在线观看免费| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 欧美成狂野欧美在线观看| 午夜成年电影在线免费观看| 久久精品91蜜桃| 亚洲成人中文字幕在线播放| 这个男人来自地球电影免费观看| 看片在线看免费视频| 久久人人精品亚洲av| 欧美+亚洲+日韩+国产| 九九热线精品视视频播放| 成年女人毛片免费观看观看9| 色老头精品视频在线观看| 色吧在线观看| 美女高潮的动态| 成人国产综合亚洲| 中文字幕人妻丝袜一区二区| 少妇的丰满在线观看| 久久国产精品影院| 舔av片在线| 欧美中文综合在线视频| 日韩欧美国产一区二区入口| 这个男人来自地球电影免费观看| 香蕉久久夜色| 国产精品一区二区三区四区久久| 久久精品国产99精品国产亚洲性色| 免费av不卡在线播放| 哪里可以看免费的av片| 精品日产1卡2卡| 午夜两性在线视频| 久久中文字幕人妻熟女| 亚洲欧美精品综合久久99| 亚洲中文字幕日韩| 99国产综合亚洲精品| 欧美成人免费av一区二区三区| 一个人看的www免费观看视频| 五月玫瑰六月丁香| 99久国产av精品| 欧美三级亚洲精品| 综合色av麻豆| 级片在线观看| 97人妻精品一区二区三区麻豆| 中亚洲国语对白在线视频| 免费观看精品视频网站| 男女下面进入的视频免费午夜| 嫁个100分男人电影在线观看| 一级毛片高清免费大全| www日本黄色视频网| 免费看十八禁软件| 看免费av毛片| 久久久国产欧美日韩av| 91久久精品国产一区二区成人 | 看免费av毛片| 搡老岳熟女国产| 男插女下体视频免费在线播放| 少妇的逼水好多| 国产精品久久久av美女十八| 成人性生交大片免费视频hd| 欧美日韩精品网址| 最近视频中文字幕2019在线8| www.熟女人妻精品国产| 亚洲自拍偷在线| 色播亚洲综合网| 亚洲欧美日韩高清在线视频| 国产伦精品一区二区三区视频9 | 亚洲国产中文字幕在线视频| 色哟哟哟哟哟哟| 日韩欧美免费精品| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月天丁香| 亚洲 欧美一区二区三区| 欧美av亚洲av综合av国产av| 最近在线观看免费完整版| 久久精品国产99精品国产亚洲性色| av天堂中文字幕网| 淫秽高清视频在线观看| 亚洲在线观看片| 美女免费视频网站| 成人av在线播放网站| 亚洲欧美日韩东京热| 99精品欧美一区二区三区四区| 国产视频内射| 国产午夜福利久久久久久| a级毛片a级免费在线| av黄色大香蕉| 国产不卡一卡二| 香蕉国产在线看| 韩国av一区二区三区四区| 午夜免费观看网址| 国产单亲对白刺激| 好男人电影高清在线观看| 亚洲精品中文字幕一二三四区| 国产成人av激情在线播放| 亚洲自拍偷在线| 免费观看的影片在线观看| netflix在线观看网站| 久久婷婷人人爽人人干人人爱| 国产真实乱freesex| 欧美日韩国产亚洲二区| 麻豆国产97在线/欧美| 欧美一级毛片孕妇| 欧美又色又爽又黄视频| 亚洲欧美日韩高清专用| 黄色成人免费大全| 九九热线精品视视频播放| 免费看a级黄色片| avwww免费| 精品一区二区三区视频在线 | 欧美+亚洲+日韩+国产| 欧美黑人巨大hd| 国产探花在线观看一区二区| 香蕉av资源在线| 精品一区二区三区四区五区乱码| av天堂中文字幕网| 国产午夜精品久久久久久| 精品熟女少妇八av免费久了| 欧美日韩亚洲国产一区二区在线观看| 中文资源天堂在线| 老汉色∧v一级毛片| 亚洲av成人av| 后天国语完整版免费观看| 天堂网av新在线| 国产极品精品免费视频能看的| 床上黄色一级片| 黄色片一级片一级黄色片| 久久久久久久久中文| 午夜免费成人在线视频| 搡老妇女老女人老熟妇| 欧美一级毛片孕妇| 欧美三级亚洲精品| 女生性感内裤真人,穿戴方法视频| 久久午夜亚洲精品久久| 波多野结衣巨乳人妻| 久久精品国产清高在天天线| 特大巨黑吊av在线直播| 男女午夜视频在线观看| 美女 人体艺术 gogo| 午夜两性在线视频| 国模一区二区三区四区视频 | 亚洲欧美一区二区三区黑人| 中文字幕高清在线视频| 亚洲国产欧美人成| 可以在线观看毛片的网站| 国产精品久久久人人做人人爽| www国产在线视频色| 亚洲 国产 在线| 给我免费播放毛片高清在线观看| 午夜两性在线视频| 亚洲精品久久国产高清桃花| 婷婷精品国产亚洲av| 91av网站免费观看| 香蕉久久夜色| 精品久久久久久久毛片微露脸| 亚洲欧美日韩高清在线视频| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 久久人人精品亚洲av| 熟女人妻精品中文字幕| 亚洲成人久久性| 嫁个100分男人电影在线观看| 亚洲午夜精品一区,二区,三区| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 黄片小视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 欧美在线黄色| 免费看日本二区| 午夜福利在线在线| 很黄的视频免费| ponron亚洲| 欧美乱妇无乱码| 女同久久另类99精品国产91| 免费看十八禁软件| 女同久久另类99精品国产91| 丝袜人妻中文字幕| 欧美日韩综合久久久久久 | 可以在线观看毛片的网站| 久久中文看片网| 亚洲男人的天堂狠狠| 亚洲美女黄片视频| 两人在一起打扑克的视频| 999久久久国产精品视频| 免费观看的影片在线观看| 黄色成人免费大全| 欧美xxxx黑人xx丫x性爽| 国产精品乱码一区二三区的特点| 在线观看美女被高潮喷水网站 | 亚洲欧美日韩高清专用| 三级国产精品欧美在线观看 | 99国产极品粉嫩在线观看| 在线观看日韩欧美| 18禁观看日本| 精品一区二区三区四区五区乱码| 婷婷精品国产亚洲av在线| 久99久视频精品免费| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲欧美98| 999久久久国产精品视频| 久9热在线精品视频| 国产极品精品免费视频能看的| 又大又爽又粗| 久久香蕉精品热| 老鸭窝网址在线观看| 国产精品一区二区三区四区免费观看 | 国产精品女同一区二区软件 | 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 国产高清视频在线观看网站| 怎么达到女性高潮| 国产一区二区在线av高清观看| 国产高清三级在线| 岛国视频午夜一区免费看| 国产综合懂色| 国产91精品成人一区二区三区| 在线播放国产精品三级| 1024香蕉在线观看| xxx96com| 99在线人妻在线中文字幕| 高潮久久久久久久久久久不卡| 免费电影在线观看免费观看| 最近最新中文字幕大全电影3| 亚洲激情在线av| 日本在线视频免费播放| 免费电影在线观看免费观看| 国产成人av激情在线播放| 亚洲精品色激情综合| 国产高清三级在线| 亚洲精品美女久久av网站| 中文字幕久久专区| 欧美日韩中文字幕国产精品一区二区三区| 色综合婷婷激情| 综合色av麻豆| 亚洲专区字幕在线| 91在线精品国自产拍蜜月 | 国产成人系列免费观看| 欧美日韩亚洲国产一区二区在线观看| 男女做爰动态图高潮gif福利片| 日本三级黄在线观看| 午夜a级毛片| 亚洲avbb在线观看| 久久婷婷人人爽人人干人人爱| 欧洲精品卡2卡3卡4卡5卡区| 国产91精品成人一区二区三区| 国产又黄又爽又无遮挡在线| 岛国视频午夜一区免费看| 亚洲在线观看片| 免费高清视频大片| 国产精品野战在线观看| 亚洲性夜色夜夜综合| 欧美日韩国产亚洲二区| 亚洲欧美日韩卡通动漫| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 中文亚洲av片在线观看爽| 俺也久久电影网| 精品一区二区三区视频在线观看免费| 精品熟女少妇八av免费久了| av天堂在线播放| 亚洲一区高清亚洲精品| 搡老岳熟女国产| 校园春色视频在线观看| 国产精品一及| 精品国产乱子伦一区二区三区| 波多野结衣高清作品| 国产欧美日韩精品一区二区| 视频区欧美日本亚洲| netflix在线观看网站| 偷拍熟女少妇极品色| 国产亚洲精品一区二区www| 搡老妇女老女人老熟妇| 一二三四在线观看免费中文在| 成人无遮挡网站| 国产视频一区二区在线看| 性色avwww在线观看| 亚洲 欧美 日韩 在线 免费| 久99久视频精品免费| 波多野结衣高清无吗| 狂野欧美白嫩少妇大欣赏| 日本五十路高清| 免费看日本二区| 最新中文字幕久久久久 | 亚洲国产精品999在线| 成熟少妇高潮喷水视频| av天堂在线播放| 香蕉国产在线看| 岛国在线观看网站| 久久久水蜜桃国产精品网| 欧美黄色淫秽网站| 午夜福利在线观看免费完整高清在 | 日韩欧美免费精品| 亚洲 欧美一区二区三区| 久久久成人免费电影| 日本黄色片子视频| 国产精品综合久久久久久久免费| 九色国产91popny在线| 国内揄拍国产精品人妻在线| 久久国产精品人妻蜜桃| 亚洲精品美女久久av网站| 人人妻人人看人人澡| 我要搜黄色片| 黄色成人免费大全| 欧美乱妇无乱码| 久久精品人妻少妇| 性欧美人与动物交配| 日本 av在线| 最近最新中文字幕大全电影3| 国产av不卡久久| netflix在线观看网站| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| 免费观看的影片在线观看| ponron亚洲| 美女高潮的动态| 精品国产三级普通话版| 国产成人av激情在线播放| 欧美性猛交黑人性爽| 欧美乱色亚洲激情| 在线观看午夜福利视频| 999久久久国产精品视频| 亚洲一区高清亚洲精品| 男女做爰动态图高潮gif福利片| 久久热在线av| 国产亚洲欧美在线一区二区| 亚洲最大成人中文| 男插女下体视频免费在线播放| 999久久久国产精品视频| 久久久国产成人精品二区| 国产高清激情床上av| 国产成人av激情在线播放| 淫秽高清视频在线观看| 成年版毛片免费区| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 黄色丝袜av网址大全| 亚洲精品美女久久av网站| 亚洲午夜理论影院| а√天堂www在线а√下载| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 91av网一区二区| 高清在线国产一区| 亚洲中文av在线| 欧美日韩国产亚洲二区| 欧美xxxx黑人xx丫x性爽| 1000部很黄的大片| 少妇的丰满在线观看| 在线永久观看黄色视频| 国产激情欧美一区二区| 搡老妇女老女人老熟妇| 亚洲成人免费电影在线观看| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 怎么达到女性高潮| 性色av乱码一区二区三区2| 一二三四在线观看免费中文在| 国产真实乱freesex| 2021天堂中文幕一二区在线观| 欧美另类亚洲清纯唯美| 亚洲av美国av| 欧美日韩国产亚洲二区| a在线观看视频网站| 日韩精品青青久久久久久| 欧美另类亚洲清纯唯美| 久久久水蜜桃国产精品网| 免费在线观看视频国产中文字幕亚洲| 这个男人来自地球电影免费观看| 国产三级中文精品| 日本三级黄在线观看| 久久香蕉国产精品| 成人永久免费在线观看视频| 18禁裸乳无遮挡免费网站照片| 97碰自拍视频| 成人18禁在线播放| 九九热线精品视视频播放| 婷婷亚洲欧美| 精品熟女少妇八av免费久了| 一a级毛片在线观看| 国产探花在线观看一区二区| 亚洲七黄色美女视频| 精品国产亚洲在线| 久久久国产成人精品二区| 丰满人妻一区二区三区视频av | 国产精品乱码一区二三区的特点| 毛片女人毛片| 亚洲成人精品中文字幕电影| 黑人欧美特级aaaaaa片| 免费大片18禁| 亚洲 欧美 日韩 在线 免费| 欧美黄色淫秽网站| 亚洲av美国av| 老司机深夜福利视频在线观看| 日韩高清综合在线| 午夜免费成人在线视频| 中文在线观看免费www的网站| 舔av片在线| 在线观看免费视频日本深夜| av片东京热男人的天堂| 黄片小视频在线播放| 亚洲无线观看免费| 桃红色精品国产亚洲av| 精品久久久久久,| 99热这里只有是精品50| 男女之事视频高清在线观看| 1000部很黄的大片| 美女扒开内裤让男人捅视频| 麻豆国产97在线/欧美| 国产精品国产高清国产av| 欧美色视频一区免费| 99re在线观看精品视频| 一级毛片精品| 国产精品 国内视频| 日韩精品青青久久久久久| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 一二三四在线观看免费中文在| 99热6这里只有精品| 久久精品国产亚洲av香蕉五月| 一级毛片精品| 麻豆av在线久日| 波多野结衣高清无吗| 中文在线观看免费www的网站| 偷拍熟女少妇极品色| 999久久久精品免费观看国产| 久久久国产成人精品二区| 波多野结衣高清作品| 色老头精品视频在线观看| 国产精品美女特级片免费视频播放器 | 狂野欧美白嫩少妇大欣赏| 国产精品久久久人人做人人爽| 91在线观看av| 男人的好看免费观看在线视频| 精品久久久久久,| 一级毛片女人18水好多| 99久久无色码亚洲精品果冻| 国产精品久久久久久亚洲av鲁大| 国产 一区 欧美 日韩| 久久中文字幕一级| 人妻夜夜爽99麻豆av| 欧美另类亚洲清纯唯美| 久久久久久大精品| 久久香蕉精品热| 久久草成人影院| 床上黄色一级片| 一二三四在线观看免费中文在| 桃色一区二区三区在线观看| 亚洲国产欧洲综合997久久,| 激情在线观看视频在线高清| 精品无人区乱码1区二区| 精品乱码久久久久久99久播| 超碰成人久久| 女警被强在线播放| 国产私拍福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲人成电影免费在线| 免费av毛片视频| 国产亚洲精品综合一区在线观看| 午夜福利视频1000在线观看| 国产一区二区在线观看日韩 | 又黄又爽又免费观看的视频| 热99re8久久精品国产| 日韩欧美三级三区| 成人国产综合亚洲| 老司机午夜十八禁免费视频| 中文字幕人妻丝袜一区二区| 白带黄色成豆腐渣| 两人在一起打扑克的视频| 长腿黑丝高跟| 亚洲专区国产一区二区| 国产单亲对白刺激| 欧美一区二区国产精品久久精品| 成人特级黄色片久久久久久久|