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

    Optimization of the grain boundary diffusion process by doping gallium and zirconium in Nd-Fe-B sintered magnets

    2023-11-02 08:12:30ZhitengLi李之藤HaiboXu徐海波FengLiu劉峰RongshunLai賴榮舜RenjieWu武仁杰ZhibinLi李志彬YangyangZhang張洋洋andQiangMa馬強
    Chinese Physics B 2023年10期
    關鍵詞:馬強劉峰洋洋

    Zhiteng Li(李之藤), Haibo Xu(徐海波), Feng Liu(劉峰), Rongshun Lai(賴榮舜), Renjie Wu(武仁杰),Zhibin Li(李志彬), Yangyang Zhang(張洋洋), and Qiang Ma(馬強),?

    1School of Rare Earths,University of Science and Technology of China,Hefei 230026,China

    2Ganjiang Innovation Academy,Chinese Academy of Sciences,Ganzhou 341000,China

    3Baotou Yunsheng Strong Magnet Material Co.,LTD,Baotou 014000,China

    Keywords: Nd-Fe-B sintered magnet,ZrB2 phase,grain boundary diffusion,micromagnetic simulation

    1.Introduction

    The Nd-Fe-B alloy has been the most excellent magnet since it was invented by Sagawaet al.[1]In recent years,the rapid development of electric vehicles[2]has brought large demand for Nd-Fe-B magnets with high coercivity.Therefore,coercivity enhancement is a focus of Nd-Fe-B research.Substitution of heavy rare earth (HRE) elements, such as Tb/Dy, for Nd is a conventional strategy that is used to enhance the coercivity of Nd-Fe-B magnets, because of the higher anisotropy field (HA) of (Tb/Dy)2Fe14B than that of Nd2Fe14B.[3,4]However,the magnetic moments of the Tb/Dy and Fe atoms are antiparallel,[4,5]which causes a significant decrease in remanence(Br).

    In 2000,Parket al.[6]deposited Dy film on a sintered Nd-Fe-B magnet followed by heat treatment, which is the genesis of the grain boundary diffusion (GBD) process.During the GBD process, Tb/Dy atoms diffuse from the surface into the magnet through GB phases, and a core-shell structure is formed.[7,8]The coercivity mechanism of sintered Nd-Fe-B is generally considered as a nucleation mechanism,[9]and the defects and corners of the grains make it easy for the demagnetization domain to nucleate.[10]The higherHAof the shell can inhibit the nucleation at the grain surface and the thicker GB phase layers formed by GBD also contribute greatly to coercivity.[11,12]Compared with the addition of Tb/Dy,the decrease in remanence is much smaller.

    However, the coercivity increment by GBD tends to decrease when the magnet thickness increases.[13]Therefore,researchers have put considerable effort into increasing the depth of diffusion.At present,a lot of research is concerned with the development of diffusion sources, which can be divided into three types,[14]including HRE-based,[15-17]light rare-earthbased[18-21]and non rare earth based[22-24]diffusion sources.In addition to the diffusion sources, studies showed that the compositions in magnets have important effects on GBD.Baeet al.[25]reported that oxygen has a negative effect on GBD.A lower Al content in the magnet was found to be beneficial to GBD by Diaoet al.[26]And Zhonget al.[27]found that the Ti element in the Nd-Fe-B has a positive effect on GBD.However, research on the influence of the magnet on GBD is still insufficient.To gain a comprehensive understanding of GBD,it is necessary to determine how the magnet itself influences GBD.

    Ga and Zr are commonly used doping elements in Nd-Fe-B,which have been well studied.[28,29]Ga doping can improve the wettability of the Nd-rich phase to form a continuous GB layer,[9,30,31]and Zr doping can refine the grains by pinning in GB phases.[32]The main purpose of Ga and Zr doping is to optimize the GB phases and microstructure of magnets,and Xuet al.[33]reported that GB phases play an important role in GBD.Ga and Zr doping may have some effects on GBD,but they are still unknown.To optimize GBD,a comprehensive study of the factors that may affect GBD is necessary.Therefore,it is worth studying the influence of these two doping elements on GBD.

    In the field of permanent magnets, micromagnetic simulation is an effective way to study the relationship between microstructures and magnetic properties.[34-36]In this paper,micromagnetic simulation was used to further analyze the different structures of the magnet based on the experiment.

    2.Experimental and simulation methods

    2.1.Sample preparation and characterization

    Three kinds of samples with the compositions of

    were prepared via a powder metallurgical process, including strip-casting,hydrogen decrepitation and jet milling.The median particle sizes were controlled to about 3.7 μm for all samples.The fine powders were compressed and oriented under a 2.0 T magnetic field, followed by isostatic pressing at 195 MPa.Next,the green compacts were sintered(1030°C-1080°C)and annealed(900°C for 1st and 500°C-600°C for 2nd)at their respective optimal temperatures to achieve the optimum magnetic properties.These samples before GBD were marked as original,Ga-doped and Zr-doped,respectively.The magnets before GBD were cut into cylinders with a size ofφ10×6 mm.The samples were polished and then washed in distilled water ultrasonically.Next,the magnets were acidized in 3%nitric acid for 50 s.Then,the magnets were cleaned with distilled water and alcohol ultrasonically.Subsequently,all the magnets were covered by Tb film via a magnetron sputtering system.A commercial Tb target(99.9%,Hefei Kejing Materials Technology Co.,Ltd.) was used to fabricate the diffusion source.During the deposition, the argon pressure was maintained at 1 Pa,and the direct current power was set to 120 W.The film weights were controlled at 0.42 wt% of the magnet for all the samples.All three magnets were then treated at 900°C for 6 h in vacuum, and annealed at the optimal temperature(500°C)for 3 h in vacuum.Samples after GBD were marked as original GBD,Ga-doped GBD and Zr-doped GBD,respectively.

    The magnetic properties of the samples were measured by a B-H tracer(NIM-500C).The microstructures of the samples were observed using a backscattered electron scanning electronic microscope (BSE-SEM) (JEOL JSM 7200F).The element distributions were obtained via electron probe micro analysis (Shimadzu EPMA-1720H).The crystal structure of micro-areas in the triangular junction phase (TJP) was characterized using transmission electron microscopy (TEM, FEI Tecnai G2 F20).

    2.2.Micromagnetic simulation

    The parallel finite element micromagnetic package Magpar[37]was used for micromagnetic simulation to understand the effect of Tb-rich shell thickness on coercivity.In this work,a model with a polycrystalline microstructure based on Voronoi tessellation[38]was established by Neper[39]and Gmsh.[40]This model was a cube with a size of 300 nm×300 nm×300 nm containing 27 irregular grains, as shown in Fig.1.Due to the limitation of computing capability, the grain size in this model was smaller than the actual grain size.In addition, the model consisted of three parts, namely, the core of Nd2Fe14B,the shell of(Tb0.25Nd0.75)2Fe14B,and the GB phases.All the grains were surrounded by a Tb-rich shell and separated by GB phases.The width of the GB was set to about 5nm.The influence of the shell thickness on coercivity was studied by changing the shell thickness from 0 nm to 10 nm.The material parameters of Nd2Fe14B and Tb2Fe14B were taken from the literature.[4,41]The saturation magnetic polarizations (JS) were 1.61 T and 0.7 T for Nd2Fe14B and Tb2Fe14B.The anisotropy constants(K1)were 4.5 MJ/m3and 5.9 MJ/m3for Nd2Fe14B and Tb2Fe14B.And the exchange stiffness (A) was 12.5 pJ/m and 13.3 pJ/m for Nd2Fe14B and Tb2Fe14B.The magnetic properties for the GB phases were defined as amorphous soft magnetic phases[15]andJS,K1,andAwere 0.15 T,0 MJ/m3, and 5 pJ/m, respectively.The magnetic parameters for (Tb0.25Nd0.75)2Fe14B were obtained by linear interpolation.

    3.Results and discussion

    3.1.Magnetic properties

    Figure 2 shows the demagnetization curves and the comparison of the coercivity of the samples with different doping elements before and after GBD.The detailed magnetic properties are listed in Table 1.With Ga or Zr doping,the coercivity decreases slightly from 14.38 kOe to 13.41 kOe or 13.63 kOe,respectively.The decrease in magnetic properties may be due to the excessive addition of Ga or Zr.[42,43]The coercivities of the three GBD samples are all enhanced significantly after diffusion, but the coercivity increments of the three samples are quite different,as shown in Fig.2(b).The coercivity of the original magnet increases by only 23% through GBD,from 14.38 kOe to 17.62 kOe.The coercivity increment of the Ga-doped magnet is slightly higher, increasing by 32% from 13.41 kOe to 17.73 kOe.In contrast,the coercivity increment for the Zr-doped magnet is 7.97 kOe,increasing by 58%from 13.63 kOe to 21.60 kOe, which is quite huge for a magnet with a thickness of 6 mm.The results show that both Ga and Zr have a positive effect on GBD,but the effect of Zr is particularly significant.

    The performances of the GBD samples at high temperatures were also measured,as shown in Fig.3.The temperature coefficients of remanence(α)and coercivity(β)in the range ofTtoTiare calculated by the following equations:[44]

    The temperature coefficient of remanence of the Ga-doped and Zr-doped GBD magnets decreases slightly compared with the original GBD magnet, from-0.104 %/°C to-0.107 %/°C and-106%/°C,respectively.However, it can be found that the temperature coefficient of coercivity of the Zr-doped GBD magnet increases significantly to-0.530 %/°C, while it is only-0.622 %/°C and-0.630 %/°C for the original GBD magnet and Ga-doped GBD magnet.The results show that the demagnetization resistance of the Zr-doped GBD magnet at high temperature is greatly improved compared with the other two samples,which can be used in high-temperature applications.

    Fig.2.(a)Demagnetization curves and(b)the coercivity comparison of the original,Ga-doped and Zr-doped magnets before and after GBD at room temperature.

    Fig.3.The temperature dependence of remanence (a) and coercivity (b) of the original GBD magnet, Ga-doped GBD magnet and Zr-doped GBD magnet.

    3.2.Microstructure characterization and analysis

    The literature indicates that the coercivity increment by GBD is derived from the Tb-rich shell and the optimization of GB phases.[12]It implies that the Ga or Zr doping may have influences on the optimization of the Tb-rich shell structure or GB phase.Therefore,the microstructure evolution of magnets doped with different elements before and after GBD was analyzed to explore the reason for the differences in the coercivity increments.Figure 4 shows the BSE-SEM images of three samples before and after GBD.There are three different contrast areas, white, dark gray and light gray, which respectively represent the GB phases, (Pr, Nd)2Fe14B main phases,and (Tb, Pr, Nd)2Fe14B shells.It can be found in Fig.4(i)that there are some regions with different contrasts in the TJP(circled in red),which may be a new phase due to the addition of Zr.Figures 4(b), 4(f), and 4(j) display the morphologies of the corresponding three magnets after GBD at the depth of 20μm,respectively.Core-shell structures can be observed in all three magnets.For all three magnets after GBD, the shell thickness decreases with the increase in diffusion depth.It is also hard to observe the clearly visible shells in the original GBD and Ga-doped GBD magnets(Figs.4(c)and 4(g))at the depth of 100 μm.However, it is very interesting that there are still many thin shells in the Zr-doped GBD magnet at a depth of 100μm, and the shells were more uniform.Even at a depth of 200μm,there are still thin shells on the surface of some grains,as shown in Fig.4(l).In contrast,Tb-rich shells almost completely disappear and the GB phase is discontinuous in the original GBD and Ga-doped GBD magnets at the depth of 200μm.This means that there is more uniform distribution in the depth of Tb in the Zr-doped GBD magnet.It is worth noting that the shell thickness on the surface of the magnet decreases with the increase in the diffusion depth, as shown in Figs.4(b),4(f),and 4(j).The shells around the main phase in Figs.4(b)and 4(f)are much thicker with 2μm.But the shell thickness in Fig.4(j) is only about 0.7 μm.Therefore,the difference in the coercivity increment by GBD of the three magnets can mainly be attributed to the distribution of Tb atoms.

    Fig.4.BSE-SEM images of the original magnet(a),the Ga-doped magnet(e),the Zr-doped magnet(i)and original GBD magnet(b)-(d),the Ga-doped GBD magnet(f)-(h),and the Zr-doped GBD magnet(j)-(l)at depths of 20μm,100μm,and 200μm.

    To more accurately describe the Tb distribution of the three magnets after GBD,EPMA characterization was carried out, as shown in Fig.5.The results show that Tb atoms are enriched on the surface of the magnet.With the increase in depth,the proportion of the Tb-rich area and the concentration of Tb drop significantly.According to Fick’s law,the concentration gradient is the main driving force for diffusion.Therefore, the decrease in the Tb concentration gradient with the increase in depth is the main reason.When the depth exceeds 100μm,Tb elements are difficult to detect in the original GBD magnet.In contrast,Tb elements can still be observed in relatively low concentrations at the depth from 100μm to 200μm in the Ga-doped GBD magnet, but the green relatively lowconcentration Tb area basically disappears when the diffusion depth exceeds 200 μm, as shown in Fig.5(b).In particular,in Zr-doped GBD magnets, when the depth exceeds 200μm,the distribution of the green relatively low-concentration Tb area can still be detected,as shown in Fig.5(c).Moreover,the uniformity of the Tb distribution in the Zr-doped GBD magnet is significantly improved.Both the diffusion depth and distribution uniformity of the Tb elements are consistent with the increment of coercivity for GBD magnets.Accordingly, the Tb element distribution and diffusion depth shown by EPMA agree with the results of BSE-SEM.

    Fig.5.EPMA mappings of Tb in the original magnet(a),the Ga-doped magnet(b)and the Zr-doped magnet(c)after GBD at the depth of 0μm to 300μm.

    Figures 6(b),6(d),and 6(g)shows the high-magnification EPMA mappings of the Tb distribution of three GBD magnets at the depth of 100μm.The mapping shows that Tb is mainly distributed along the GB phase and exhibits a dispersion distribution.It is difficult to distinguish the core-shell structure in the original GBD magnet.In addition,the Ga-doped GBD magnet exhibits a discontinuous Tb-rich network.In contrast,the Tb-rich shells of the Zr-doped GBD magnet are complete and continuous.Similar results can be observed in the corresponding SEM in Figs.6(a), 6(c), and 6(f).Specifically, the concentration and distribution uniformity of Tb in the original GBD magnet,Ga-doped GBD magnet and Zr-doped GBD magnet at the depth of 100 μm are significantly improved in turn.This indicates that Zr doping can effectively optimize the distribution of Tb.The nucleation field of the reversed domain on grain surfaces is enhanced when it is covered by shell,and a complete shell can reduce the nucleation sites.[11]

    To understand how the doping elements Ga and Zr affect diffusion depth,the distributions of Ga and Zr were also analyzed via EPMA,as shown in Figs.6(e)and 6(h).It is found that Ga atoms are distributed in the TJP,as shown in Fig.6(e).The melting point of Ga is low, which can lower the melting point and improve the fluidity of GB phases.GB phases are the channels for GBD, and the decrease in the melting point of the GB can facilitate the migration of Tb atoms.Similar to the distribution of Ga,Zr also exists in the TJP in the Zr-doped GBD magnet,as shown in Fig.6(h).Figure 6(i)shows that the B elements are also enriched in the Zr-rich region, indicating the formation of the Zr-B compound, which corresponds to the phase with different contrast in Fig.4(i)mentioned above.

    Fig.6.SEM images[(a), (b), (c)]and EPMA mappings of Tb[(b), (d), (g)], Ga(e), Zr(h), and B(i)of the original GBD magnet, the Ga-doped GBD magnet and the Zr-doped GBD magnet at the depth of 100μm.

    The Zr-B compound in the TJP may be the main reason for the different distribution of Tb.To further analyze the role of Zr-B compounds in the TJP,TEM analysis was performed.Figure 7(a)shows the TJP that contains the Zr-B phase.The Zr-B phase has a strip morphology and is completely encapsulated inside the TJP.According to Fig.7(b),the selected area electron diffraction(SAED)of the Zr-rich area shows that this area in the TJP is consistent with the structural characteristics of hexagonal ZrB2,whose melting point is about 3000°C.Because the formation of ZrB2consumes boron,the consumption of boron leads to the reduction of Nd2Fe14B,[29]thus more GB phases are formed.The distribution of GB phases observed in Figs.4 and 6 also verifies this result.Therefore, the formation of ZrB2can provide a continuous GB structure for GBD.Moreover,we find that ZrB2is completely encased in the TJP.There are also many lattice dislocations between the TJP and ZrB2phase,as shown in Fig.7(c),resulting in a high concentration of defects at the interface,[27]which have a high Gibbs free energy.Therefore,Tb atoms need less diffusion activation energy at defect sites,which accelerates atomic diffusion.

    3.3.Micromagnetic simulation results

    Zr doping can significantly increase the diffusion depth,while the shell thickness near the surface also decreases, as shown in Fig.4.To find out whether the reduction in shell thickness would have a negative effect on coercivity, we carried out a micromagnetic simulation.Figure 8 shows that the coercivity is in accordance with the thickness of the shell.However, the effect of shell thickness on coercivity becomes small when the thickness exceeds 5 nm.In the range of 0 nm to 5 nm,the coercivity increases by 1975 Oe/nm,but by only 425 Oe/nm in the range of 5 nm to 10 nm.According to the simulation results, thin shells are enough to provide a large coercivity increase.The decrease in shell thickness within a certain range does not have a significant negative effect on coercivity.Therefore, the distribution of Tb in the Zr-doped GBD magnet is more advantageous to obtain higher coercivity, i.e., more Tb elements diffusing into the interior of the magnet to form thin shells.With the same amount of the Tb diffusion source(0.42%in this work),the thin shell also means that there can be more Tb surplus to provide the concentration gradient of Tb diffusion to the deeper part of the magnet.

    Fig.8.(a) Demagnetization curves and (b) the coercivity increment for models with shell thicknesses of 0 nm,2.5 nm,5 nm,7.5 nm,and 10 nm,calculated by Magpar.

    Fig.9.Magnetization reversal patterns on the outer surfaces(a)and crosssections of x=0.9(b)and z=0.93(c)of models with shell thicknesses of 0 nm,5 nm,and 10 nm.

    Moreover,the magnetization reversal patterns are shown in Fig.9.A clear difference can be observed between the magnets before and after GBD.In the model without Tb-rich shells, the magnetic moment begins to reverse at the GB and the edges of the grain surface (blue contrast).In the models with shells,the reversed domains move from the edges of the grain surface to the edge of the core not covered by the Tb-rich shells.The magnetization reversal patterns are also similar for the models with 5 nm shells and 10 nm shells,except that the nucleation site moves further into the grain in the model with 10 nm shells.Due to the influence of the demagnetizing field,the nucleation field at the grain surface is low.After GBD,the grain surface is covered by Tb-rich shell, which significantly improves the nucleation field at the grain surface.With a further increase in the shell thickness, the nucleation site moves further into the grain.However,from the surface of the magnet to the interior,the demagnetizing field gradually decreases and the distribution of this field inside the magnet is more uniform.[45]Thus,the nucleation field does not increase much when the nucleation site moves into the grain, which can explain why the coercivity increases little as the shell thickness increases further.Therefore, a 5 nm thin shell is enough to achieve high coercivity.Also,the thicker Tb-shells on the surfaces of the original GBD and Ga-doped GBD magnets result in a waste of Tb.

    4.Conclusion

    In summary, we found that Ga doping can decrease the melting point and facilitate the fluidity of GB phases,making it easier to form continuous GB phases.Therefore,Ga doping in magnets can slightly improve the coercivity increment by GBD, from 3.24 kOe to 4.32 kOe.The consumption of B in the Zr-doped magnet leads to an increase in the GB phases,which can promote the formation of continuous GB phases.More important is the increase in the defect concentration between ZrB2phases and GB phases in the TJP, causing high Gibbs free energy in the TJP, which makes the Tb atoms require less diffusion activation energy during migration.Therefore,Zr doping improves the diffusion depth and uniformity of Tb,which is the direct reason that allows Zr-doped magnets to achieve a 7.97 kOe increment by GBD technology.In the Zrdoped magnet, the increase in diffusion depth will lead to a decrease in shell thickness near the surface.In addition, the reduction in shell thickness within a certain range was found not to have a significant negative effect on coercivity by micromagnetic simulation.This study helps to understand the role of Ga and Zr in the GBD process,and may provide a new perspective for researchers to improve the effects of GBD.

    Acknowledgments

    Project supported by the National Natural Science Foundation of China(Grant No.52261037),self-deployed Projects of Ganjiang Innovation Academy, Chinese Academy of Sciences(Grant No.E055B002),the Project of Baotou City Science and Technology(Grant No.XM2022BT04),the Key Research Program of the Chinese Academy of Sciences (Grant No.ZDRW-CN-2021-3) and the Key Research Project of Jiangxi Province(Grant No.20203ABC28W006).

    The authors would like to thank Professor Zhenchen Zhong,Jiangxi University of Science and Technology,for providing guidance for this work.

    猜你喜歡
    馬強劉峰洋洋
    Effects of O2 addition on the plasma uniformity and reactivity of Ar DBD excited by ns pulsed and AC power supplies
    History of World Smile Day
    陌上
    洋洋兔 編繪
    一鳴驚人
    劉峰:從“發(fā)明狂人”到“創(chuàng)業(yè)達人”
    金色年華(2016年14期)2016-02-28 01:44:29
    竭澤而漁(下)
    相逢在海邊(花腔女高音獨唱)
    音樂生活(2014年11期)2014-04-25 12:43:32
    特殊的孩子
    故事林(2013年9期)2013-05-14 17:30:12
    這輩子
    藝術評鑒(2011年7期)2011-11-29 07:58:18
    久久精品91无色码中文字幕| 亚洲专区中文字幕在线| 欧美激情久久久久久爽电影 | 高清在线国产一区| 99香蕉大伊视频| 国产成人影院久久av| 91成年电影在线观看| 亚洲精品粉嫩美女一区| 久久亚洲精品不卡| 99久久人妻综合| 曰老女人黄片| 91九色精品人成在线观看| 亚洲 欧美一区二区三区| 日韩欧美在线二视频 | 精品国产美女av久久久久小说| 精品亚洲成a人片在线观看| 成人av一区二区三区在线看| 亚洲精品在线美女| 国产亚洲一区二区精品| 中文亚洲av片在线观看爽 | 午夜福利免费观看在线| 国产一区二区三区综合在线观看| 亚洲第一av免费看| 精品亚洲成国产av| 国产亚洲一区二区精品| 一边摸一边做爽爽视频免费| 国产亚洲精品一区二区www | 99re6热这里在线精品视频| 一区二区日韩欧美中文字幕| 国产真人三级小视频在线观看| 在线视频色国产色| videos熟女内射| 老熟女久久久| 十八禁高潮呻吟视频| www.自偷自拍.com| 天天躁夜夜躁狠狠躁躁| 动漫黄色视频在线观看| 人妻丰满熟妇av一区二区三区 | av不卡在线播放| 亚洲成国产人片在线观看| 亚洲熟女精品中文字幕| 国产亚洲欧美在线一区二区| 亚洲精品一二三| 在线永久观看黄色视频| 久久狼人影院| 亚洲精品久久成人aⅴ小说| 新久久久久国产一级毛片| ponron亚洲| 黄色 视频免费看| 婷婷成人精品国产| 国产99白浆流出| 亚洲精品中文字幕在线视频| 妹子高潮喷水视频| 美女视频免费永久观看网站| 亚洲欧美激情在线| 亚洲三区欧美一区| 又紧又爽又黄一区二区| 亚洲欧美精品综合一区二区三区| 黄色成人免费大全| 女同久久另类99精品国产91| 一级a爱视频在线免费观看| 老司机深夜福利视频在线观看| 韩国av一区二区三区四区| 女人爽到高潮嗷嗷叫在线视频| 男男h啪啪无遮挡| 精品国产乱子伦一区二区三区| 淫妇啪啪啪对白视频| 久久青草综合色| 久久青草综合色| 中文字幕高清在线视频| 9191精品国产免费久久| 欧美日韩国产mv在线观看视频| 免费一级毛片在线播放高清视频 | 日韩欧美三级三区| 亚洲久久久国产精品| 捣出白浆h1v1| 捣出白浆h1v1| 国产精品.久久久| 午夜福利在线免费观看网站| 亚洲人成77777在线视频| 日韩视频一区二区在线观看| 黄色视频,在线免费观看| 日韩熟女老妇一区二区性免费视频| 黄色视频,在线免费观看| 俄罗斯特黄特色一大片| 国产男女超爽视频在线观看| 99精国产麻豆久久婷婷| 欧美激情高清一区二区三区| 国产伦人伦偷精品视频| 高清黄色对白视频在线免费看| 一级片免费观看大全| 一级片免费观看大全| 这个男人来自地球电影免费观看| 久久久久国内视频| 欧美日韩福利视频一区二区| 亚洲国产看品久久| 午夜免费观看网址| 热99久久久久精品小说推荐| 黄片小视频在线播放| 在线天堂中文资源库| 岛国毛片在线播放| 一个人免费在线观看的高清视频| 狠狠婷婷综合久久久久久88av| 高清黄色对白视频在线免费看| 国产精品久久视频播放| 国产av又大| 无人区码免费观看不卡| aaaaa片日本免费| 最新在线观看一区二区三区| 久久草成人影院| 成人18禁在线播放| 精品久久久久久久久久免费视频 | 亚洲一区二区三区欧美精品| 中文欧美无线码| 欧美日韩福利视频一区二区| 精品午夜福利视频在线观看一区| 免费久久久久久久精品成人欧美视频| 国产精品一区二区在线不卡| 黄色a级毛片大全视频| 色综合婷婷激情| 免费日韩欧美在线观看| 午夜两性在线视频| 国产一区二区三区视频了| 国产深夜福利视频在线观看| 久久国产乱子伦精品免费另类| 日韩成人在线观看一区二区三区| 国产不卡一卡二| 日韩欧美一区视频在线观看| 国产免费av片在线观看野外av| 80岁老熟妇乱子伦牲交| 国产精品久久久av美女十八| 精品一品国产午夜福利视频| 国产精品国产av在线观看| 男女床上黄色一级片免费看| 黑丝袜美女国产一区| cao死你这个sao货| 国产欧美日韩一区二区精品| 精品国产一区二区久久| 精品午夜福利视频在线观看一区| 亚洲精品国产色婷婷电影| 黄色怎么调成土黄色| 91在线观看av| 国产精品1区2区在线观看. | 国产欧美日韩一区二区三| 天堂动漫精品| 午夜福利免费观看在线| 久久精品国产亚洲av高清一级| 69av精品久久久久久| 亚洲人成电影观看| 久久久国产成人精品二区 | 精品午夜福利视频在线观看一区| 欧美精品av麻豆av| 看黄色毛片网站| 亚洲色图av天堂| 老汉色∧v一级毛片| 午夜影院日韩av| 日韩欧美三级三区| 黄色视频,在线免费观看| 欧美精品人与动牲交sv欧美| 99久久国产精品久久久| 精品午夜福利视频在线观看一区| 久久热在线av| 王馨瑶露胸无遮挡在线观看| 女性被躁到高潮视频| 不卡一级毛片| 天天影视国产精品| 久久久久精品人妻al黑| 俄罗斯特黄特色一大片| 最新的欧美精品一区二区| 亚洲人成伊人成综合网2020| 国产精品欧美亚洲77777| 99国产精品99久久久久| 又黄又爽又免费观看的视频| 免费在线观看影片大全网站| 精品国产一区二区三区四区第35| 男女床上黄色一级片免费看| 99国产精品一区二区蜜桃av | 亚洲专区国产一区二区| 纯流量卡能插随身wifi吗| 波多野结衣av一区二区av| www.999成人在线观看| 欧美在线一区亚洲| 一进一出抽搐动态| 亚洲欧洲精品一区二区精品久久久| 欧美日韩国产mv在线观看视频| 不卡一级毛片| 精品国产一区二区三区四区第35| 国产一区二区激情短视频| 天天操日日干夜夜撸| 日本一区二区免费在线视频| 色尼玛亚洲综合影院| 精品久久久精品久久久| 精品国产亚洲在线| 精品亚洲成国产av| 99在线人妻在线中文字幕 | 亚洲午夜理论影院| 天堂俺去俺来也www色官网| 亚洲美女黄片视频| 久久久久精品人妻al黑| 色精品久久人妻99蜜桃| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品久久久久5区| 视频区图区小说| 欧美精品av麻豆av| 亚洲自偷自拍图片 自拍| 久久婷婷成人综合色麻豆| 国产xxxxx性猛交| av有码第一页| 桃红色精品国产亚洲av| 欧美一级毛片孕妇| xxxhd国产人妻xxx| 在线观看66精品国产| 在线视频色国产色| 国产欧美亚洲国产| 久久久久久久午夜电影 | av线在线观看网站| 1024香蕉在线观看| 中国美女看黄片| 亚洲成a人片在线一区二区| 脱女人内裤的视频| 国产亚洲av高清不卡| 久久国产精品男人的天堂亚洲| 美女高潮到喷水免费观看| 亚洲一区二区三区欧美精品| 久久精品国产99精品国产亚洲性色 | 51午夜福利影视在线观看| 69精品国产乱码久久久| 成人18禁在线播放| 欧美 亚洲 国产 日韩一| 悠悠久久av| 亚洲视频免费观看视频| cao死你这个sao货| 少妇的丰满在线观看| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区mp4| 人人妻,人人澡人人爽秒播| 亚洲男人天堂网一区| 老鸭窝网址在线观看| 精品久久久久久电影网| 精品国产乱子伦一区二区三区| www日本在线高清视频| 老熟女久久久| 国产精品二区激情视频| 亚洲一区二区三区欧美精品| 国产单亲对白刺激| 亚洲中文日韩欧美视频| 日本黄色视频三级网站网址 | 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 欧美日韩瑟瑟在线播放| 新久久久久国产一级毛片| 亚洲中文av在线| 国产色视频综合| 高清在线国产一区| 欧美不卡视频在线免费观看 | 99久久综合精品五月天人人| 美女福利国产在线| 国产一区二区三区视频了| 国产一卡二卡三卡精品| 一边摸一边抽搐一进一小说 | 在线视频色国产色| 黑丝袜美女国产一区| 国产熟女午夜一区二区三区| xxx96com| 黑人巨大精品欧美一区二区mp4| 国产蜜桃级精品一区二区三区 | 王馨瑶露胸无遮挡在线观看| 亚洲精品在线美女| 久久中文看片网| 日韩人妻精品一区2区三区| 精品久久蜜臀av无| 日韩免费高清中文字幕av| av线在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 视频在线观看一区二区三区| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 亚洲专区中文字幕在线| 国产欧美亚洲国产| 午夜精品久久久久久毛片777| 人人澡人人妻人| 国产亚洲一区二区精品| 热99久久久久精品小说推荐| 母亲3免费完整高清在线观看| 色综合欧美亚洲国产小说| 在线观看舔阴道视频| 精品电影一区二区在线| 亚洲精品国产精品久久久不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 久久青草综合色| 中文字幕制服av| 在线天堂中文资源库| 亚洲人成电影免费在线| 国产在线精品亚洲第一网站| 亚洲熟女毛片儿| 日韩欧美三级三区| 国产色视频综合| 午夜免费鲁丝| 亚洲精品中文字幕在线视频| 岛国毛片在线播放| 人人妻人人澡人人看| 啦啦啦视频在线资源免费观看| 精品电影一区二区在线| 捣出白浆h1v1| 制服诱惑二区| 淫妇啪啪啪对白视频| 亚洲精品国产区一区二| 亚洲av成人av| 亚洲人成77777在线视频| 久久精品国产99精品国产亚洲性色 | 亚洲情色 制服丝袜| 欧美成人午夜精品| 亚洲精品久久成人aⅴ小说| av网站免费在线观看视频| 法律面前人人平等表现在哪些方面| 曰老女人黄片| 人妻 亚洲 视频| 电影成人av| 制服人妻中文乱码| 一边摸一边做爽爽视频免费| 他把我摸到了高潮在线观看| 亚洲一区二区三区欧美精品| 成人影院久久| 精品福利观看| 丰满的人妻完整版| a在线观看视频网站| 亚洲精品久久午夜乱码| 精品国产超薄肉色丝袜足j| 一本一本久久a久久精品综合妖精| 国产av又大| 91精品国产国语对白视频| 老熟妇仑乱视频hdxx| 五月开心婷婷网| 精品久久久久久电影网| 超碰成人久久| 国产淫语在线视频| 一级片免费观看大全| 91大片在线观看| 成人av一区二区三区在线看| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 超碰成人久久| 国产高清激情床上av| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 18禁观看日本| 亚洲九九香蕉| 夜夜夜夜夜久久久久| 黄色丝袜av网址大全| 亚洲欧美精品综合一区二区三区| 涩涩av久久男人的天堂| 99精品欧美一区二区三区四区| 久久香蕉激情| 人人妻人人添人人爽欧美一区卜| 飞空精品影院首页| 国产成人精品在线电影| 在线av久久热| 9热在线视频观看99| 91麻豆精品激情在线观看国产 | 最近最新免费中文字幕在线| av一本久久久久| 999精品在线视频| 老司机影院毛片| 欧美一级毛片孕妇| 一级黄色大片毛片| 宅男免费午夜| 久久香蕉激情| 9191精品国产免费久久| 丁香欧美五月| 亚洲黑人精品在线| 久久热在线av| 黑人巨大精品欧美一区二区蜜桃| 精品国产乱码久久久久久男人| 男女高潮啪啪啪动态图| 老熟妇仑乱视频hdxx| 精品国产一区二区三区久久久樱花| 色综合欧美亚洲国产小说| 在线观看66精品国产| www.999成人在线观看| 很黄的视频免费| 免费在线观看视频国产中文字幕亚洲| 国产不卡一卡二| 国产欧美日韩一区二区精品| 午夜91福利影院| 在线观看免费视频网站a站| 波多野结衣一区麻豆| 看片在线看免费视频| videos熟女内射| 国产精品 国内视频| 黄色a级毛片大全视频| 婷婷精品国产亚洲av在线 | 好看av亚洲va欧美ⅴa在| 国产不卡一卡二| 亚洲色图 男人天堂 中文字幕| 色综合婷婷激情| 精品国产一区二区三区四区第35| 这个男人来自地球电影免费观看| 黄色怎么调成土黄色| av线在线观看网站| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 好男人电影高清在线观看| 少妇裸体淫交视频免费看高清 | 日本黄色视频三级网站网址 | 黄色a级毛片大全视频| 日本五十路高清| 大型av网站在线播放| 精品国内亚洲2022精品成人 | 国产精品一区二区在线观看99| 老司机福利观看| 久久久精品区二区三区| 18禁美女被吸乳视频| 我的亚洲天堂| 久久久国产欧美日韩av| 欧美精品人与动牲交sv欧美| 人人澡人人妻人| 免费少妇av软件| 91国产中文字幕| 美国免费a级毛片| 1024香蕉在线观看| 久久草成人影院| 成人特级黄色片久久久久久久| 久久国产精品大桥未久av| 国产精品免费大片| 一进一出抽搐gif免费好疼 | 18禁美女被吸乳视频| 久久久久久久久免费视频了| 男女午夜视频在线观看| a级毛片在线看网站| 成人精品一区二区免费| avwww免费| bbb黄色大片| 国产乱人伦免费视频| 国产午夜精品久久久久久| 欧美成人免费av一区二区三区 | 他把我摸到了高潮在线观看| 精品无人区乱码1区二区| 亚洲精品成人av观看孕妇| videosex国产| 成人18禁在线播放| 午夜免费成人在线视频| 成年人黄色毛片网站| 久久 成人 亚洲| 热99久久久久精品小说推荐| 三级毛片av免费| 一级a爱视频在线免费观看| 色94色欧美一区二区| 色尼玛亚洲综合影院| videos熟女内射| 国产成人影院久久av| 亚洲欧美一区二区三区久久| 80岁老熟妇乱子伦牲交| 很黄的视频免费| 午夜两性在线视频| 久久中文字幕一级| 人人妻人人添人人爽欧美一区卜| 在线视频色国产色| 中文字幕高清在线视频| 日韩欧美一区视频在线观看| 国产成人av教育| 久久香蕉精品热| 久久人妻熟女aⅴ| 人人澡人人妻人| 精品国产一区二区三区四区第35| 亚洲一卡2卡3卡4卡5卡精品中文| 手机成人av网站| 女性生殖器流出的白浆| www.精华液| 日韩制服丝袜自拍偷拍| 亚洲五月婷婷丁香| 超碰成人久久| 国产三级黄色录像| 午夜成年电影在线免费观看| 精品人妻1区二区| 99riav亚洲国产免费| 午夜免费成人在线视频| 在线观看一区二区三区激情| 99精品久久久久人妻精品| 亚洲精品一卡2卡三卡4卡5卡| 国产av精品麻豆| 国产在线精品亚洲第一网站| 亚洲国产精品sss在线观看 | 亚洲国产看品久久| 国产麻豆69| 1024视频免费在线观看| 脱女人内裤的视频| aaaaa片日本免费| 国产欧美亚洲国产| 精品人妻熟女毛片av久久网站| 99国产精品免费福利视频| 日本wwww免费看| 99久久精品国产亚洲精品| 精品一区二区三区av网在线观看| 久热爱精品视频在线9| www.999成人在线观看| 欧美丝袜亚洲另类 | 侵犯人妻中文字幕一二三四区| 亚洲精品粉嫩美女一区| 91麻豆av在线| 成人av一区二区三区在线看| 男女下面插进去视频免费观看| 亚洲国产欧美网| 在线天堂中文资源库| 久久久国产一区二区| 国产日韩欧美亚洲二区| 高清毛片免费观看视频网站 | 久久久久久久久免费视频了| 一区福利在线观看| 女人被狂操c到高潮| 久久婷婷成人综合色麻豆| 中文字幕另类日韩欧美亚洲嫩草| 国产成人av激情在线播放| 欧美午夜高清在线| 777米奇影视久久| 丁香欧美五月| 欧美丝袜亚洲另类 | 一个人免费在线观看的高清视频| 国产区一区二久久| 国产亚洲一区二区精品| av超薄肉色丝袜交足视频| 免费在线观看视频国产中文字幕亚洲| 国产精品一区二区在线观看99| 黄色怎么调成土黄色| 国产精品久久久人人做人人爽| 日本wwww免费看| 精品一区二区三区视频在线观看免费 | 成在线人永久免费视频| 丰满人妻熟妇乱又伦精品不卡| 91字幕亚洲| 91成年电影在线观看| 美女 人体艺术 gogo| 女性被躁到高潮视频| 91精品三级在线观看| aaaaa片日本免费| 国产亚洲精品久久久久久毛片 | 一级片'在线观看视频| 国产无遮挡羞羞视频在线观看| 在线免费观看的www视频| 女人精品久久久久毛片| 国产99白浆流出| 老鸭窝网址在线观看| 成人影院久久| 韩国精品一区二区三区| 免费人成视频x8x8入口观看| 夜夜躁狠狠躁天天躁| 18禁观看日本| 人成视频在线观看免费观看| 免费在线观看日本一区| 91精品三级在线观看| 精品高清国产在线一区| 老司机影院毛片| 91在线观看av| 久久久国产精品麻豆| 色播在线永久视频| 五月开心婷婷网| 国产成人精品久久二区二区免费| 天天添夜夜摸| 91九色精品人成在线观看| 欧美日韩亚洲综合一区二区三区_| 日日摸夜夜添夜夜添小说| 美女高潮到喷水免费观看| 久久 成人 亚洲| 国产男女超爽视频在线观看| 好看av亚洲va欧美ⅴa在| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩亚洲综合一区二区三区_| 中文字幕高清在线视频| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清在线视频| 国产一区二区激情短视频| 亚洲精品av麻豆狂野| 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清在线视频| 成年女人毛片免费观看观看9 | 成人18禁在线播放| 国产蜜桃级精品一区二区三区 | 少妇被粗大的猛进出69影院| 最新美女视频免费是黄的| a级片在线免费高清观看视频| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 久久青草综合色| 日韩熟女老妇一区二区性免费视频| 国产成+人综合+亚洲专区| 日韩欧美免费精品| 久久精品亚洲熟妇少妇任你| 一夜夜www| 亚洲精品一卡2卡三卡4卡5卡| 两性夫妻黄色片| av线在线观看网站| 久久久久国产一级毛片高清牌| 亚洲国产看品久久| 国产高清国产精品国产三级| 捣出白浆h1v1| 亚洲av成人一区二区三| 一夜夜www| 国产亚洲av高清不卡| 免费观看a级毛片全部| 亚洲 国产 在线| 国产成人欧美在线观看 | 两个人看的免费小视频| 国产欧美亚洲国产| 久久精品亚洲精品国产色婷小说| 狠狠狠狠99中文字幕| 亚洲欧美激情在线| 一级作爱视频免费观看| 亚洲精品在线观看二区| 美女扒开内裤让男人捅视频| 国产91精品成人一区二区三区| 久久影院123| 亚洲av日韩精品久久久久久密| 国产免费男女视频| 亚洲国产欧美一区二区综合| 不卡一级毛片| 女警被强在线播放|