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

    Tungsten combustion in impact initiated W-Al composite based on W(Al) super-saturated solid solution

    2023-07-31 13:30:36KongxunZhaoXiaohongZhangXiaoranGuYuTangShunLiYicongYeLianZhuShuxinBai
    Defence Technology 2023年7期

    Kong-xun Zhao,Xiao-hong Zhang,Xiao-ran Gu,Yu Tang,Shun Li ,Yi-cong Ye,Li'an Zhu,Shu-xin Bai

    College of Aerospace Science and Engineering, National University of Defense Technology, Changsha, 410073, China

    Keywords:Tungsten combustion Reactive materials Super-saturated solid solution Shock-induced reactions W-Al composite

    ABSTRACT Element W can effectively improve the density of energetic structural materials.However,W is an inert element and does not combust in air.To change the reaction characteristics of W,60 at.% Al was introduced into W through mechanical alloying.XRD analysis shows that after 50 h of ball milling,the diffraction peak of Al completely disappears and W(Al60) super-saturated solid solution powder is obtained.Further observation by HAADF and HRTEM reveals that the W(Al60)super-saturated solid solution powder is a mixture of solid solution and amorphous phase.Based on the good thermal stability of W(Al60) alloy powder below 1000 °C,W(Al60)-Al composite was synthesized by hot pressing process.Impact initiation experiments suggest that the W(Al60)-Al composite has excellent reaction characteristics,and multiple types of tungsten oxides are detected in the reaction products,showing that the modified W is combustible in air.Due to the combustion of tungsten,the energy release rate of the W(Al60)-Al composite at speed of 1362 m/s reaches 2.71 kJ/g.

    1.Introduction

    Energetic structural materials (ESMs) are a new kind of energetic material with high strength.They have high stability under normal conditions,and can maintain structural integrity under detonation loading[1].However,when interacting with the target object,chemical reactions will occur between different components of the ESMs or between components and the air,resulting in combustion or even explosion reactions [2-5].Using ESMs to replace the inert parts of existing warheads,such as kinetic penetrator,shaped charge liners and reactive fragments,can not only impose traditional kinetic energy penetration on the target,but also bring additional chemical damage,such as high temperature and overpressure [6-8].Therefore,the development of ESMs has become an important technical way to enhance the damage effect of warheads.

    In order to obtain excellent energy release characteristics,most ESMs are composed of active metals,including Al,Zr,Ni,Mg,B,etc.[9,10].Al-based ESMs have high energy density and excellent energy release characteristics,so it has attracted much attention in the field of ESMs [11-13].Typical Al-based ESMs,such as Al-Ni[14,15],Al-Ta [16],Al-PTFE [17-19],Al-Fe2O3 [20,21],have been studied extensively.For example,drop hammer test shows that the Al/PTFE composite has a strong temperature destruction effect,the highest flame temperature of PTFE/Al/Si/CuO containing 20%PTFE/Si can reach 2589 K.Ren [18] fabricated Al/PTFE composites by mold pressing and sintering method and investigated the dynamic mechanical properties and ignition behavior of the materials under extreme conditions.The results show that Al/PTFE composite is an low threshold ESMs,and the reaction sensitivity increases with the increase of temperature.Xiong [22] studied the influence of the introduction of Cu and PTFE on the shock-induced chemical reaction properties of Al/Ni Composites.The results show that PTFE can improve the sensitivity of Al/Ni Composite under impact,while Cu can prolong the duration of the reaction.

    Al-based ESMs have good energy release characteristics,but its application is limited by low density,low strength and low permeability.Taking the Al/Ni composite as an example,its tensile strength is generally lower than 200 MPa,and its compressive strength is between 250 MPa and 500 MPa [23-25].For Al/PTFE ESMs,the quasi-static compressive and tensile strength of Al/PFTE is even less than 100 MPa,which makes it difficult to provide effective penetration capability [26-29].At present,the density and penetration ability of Al-based alloys are mainly improved by adding high density elements such as W[30,31].The experimental results of split Hopkinson bar indicate that the dynamic compressive strength of the W-PTFE-Al composites increases significantly with the content of W increased from 50 to 80 wt%.However,impact initiated experiments show that the reaction energy of WPTFE-Al composites decreases with the increase of W in oxygen and argon atmosphere[32].This is mainly because W has high thermal inertia and low adiabatic flame temperature in air,which impede its ignition and self-sustaining combustion [33,34].But in view of thermodynamics,the oxidation heat release per unit volume of W is 88.5 kJ/cm3,second only to boron and tantalum[35].

    Element Al is easy to burn and has a high combustion heat(84.0 kJ/cm3)[36].If a large number of Al atoms are dissolved into the W lattice and form a metastable W(Al) alloy powder,the high activity and high heat release of Al may induce the combustion of modified W.According to the Al-W binary phase diagram,the solid solubility limit of Al in W is about 15.9 at% at 1255°C,which may not be enough to realize the ignition of W.To achieve the dissolution of more Al atoms,mechanical alloying (MA) is a feasible technical means[37].Many scholars have pointed out that the solid solubility expansion of binary systems can be achieved through MA,such as W-Ag [38],Cu-Nb [39] Fe-Zr [40] and W-Al.Ouyang [41]first proposed that the solubility limit of Al in W could reach 50 at.%by MA.Rafiei[42]prepared and studied the formation mechanism of W-20 wt% Al,the result reveals that the formation of the fine layered microstructure play a great role in promoting the diffusion of Al atoms into W.Recent studies on W(Al) super-saturated solid solution powders mainly focus on the microstructure [43,44],oxidation resistance [45] and sinterability [46],and there is no relevant report on the combustion characteristics of W.

    To improve the reaction characteristics of W,this work attempts to introduce Al atoms into the W lattice and prepare W(Al) supersaturated solid solution powder by the MA process.The milling process,crystal structure and thermal stability of the W(Al) alloy powder are also studied.The reaction characteristics of the W(Al)-Al composite and the combustion characteristics of modified W(Al)alloy powder are tested through impact initiation experiments.

    2.Materials and methods

    In this study,W-60 at% Al (W(Al60)) super-saturated solid solution powder was prepared via mechanical alloying.Fig.1 shows the morphology of original materials.W particles (Tijo metal materials,3-4 μm,99.9%purity)and Al particles(Tijo metal materials,5-6 μm,99.5% purity) were blended to obtain the desired composition.This experiment used a YXQM-4L planetary ball mill(Changsha Miqi Instrument Equipment Co,Ltd) equipped with 500 mL stainless steel milling jar.The milling balls are stainless steel balls with diameters of 5 mm,8 mm and 10 mm,which are mixed in the proportion of 5:3:2.The experiments were conducted under an argon atmosphere at 250 rpm with a ball-to-powder mass ratio of 10:1.Alcohol was used as the process control agent(PCA)to prevent powder agglomeration.After each 10 min of milling,the planetary ball mill stopped for 10 min to avoid the volatilization of PCA induced by excessive temperature.The density of obtained alloy powder was measured by pycnometer method.The phase evolution during ball milling was detected by D/max2550 VB X-ray diffraction (XRD) with CuKα (λ=0.1542 nm) radiation in the 2θ range of 30-90°.The grain size and lattice strain of the alloy powder were calculated using the Scherrer equation and the Williamson-Hall method,respectively.The morphology evolution and elemental distribution of particles were characterized by scanning electron microscopy (SEM) equipped with energy dispersive spectroscopy(EDS).The crystal structures were analyzed by high-resolution transmission electron microscopy (HRTEM).In addition,STA449-F3 thermo-gravimetry differential scanning calorimeter (DSC) was used to study the thermal stability of alloy powder between 20 and 1000°C.

    Fig.1.Morphology of original powders: (a) W and (b) Al.

    Using 40 vol% pure Al powder as a binder,W(Al60)-Al composite was fabricated by hot pressing at 400°C for 2 h under 400 MPa.The quasi-static compression tests were conducted on a universal testing machine (WDW-100 B) under velocity control with a nominal strain rate of 10-3s-1.The impact initiation reaction of the W(Al60)-Al composite was verified by a ballistic gun test system.Fig.2 shows the schematic diagram of the impact initiation experimental device.14.5 mm ballistic gun,velocity measurement and 27 L test chamber were set in a straight line.The energetic projectiles were driven by ballistic gun.After penetrating the front steel plate and impacting the steel anvil,they fractured and began to react.By adjusting the mass of black powder,the initial velocity of the projectiles can be controlled.The reaction phenomenon in the chamber was recorded by a high speed camera,and the energy release characteristics were evaluated by pressure sensors.

    Fig.2.The schematic diagram of the experiment setup.

    3.Results and discussion

    3.1.Characterization of W(Al60) alloy powder

    The XRD patterns of W(Al60) mixed powders with different milling time are shown in Fig.3(a).After milling for 10 h,only three major diffraction peaks of Al remained,while the rest of Al peaks disappeared.By further increasing the milling time to 40 h,the diffraction peak of Al only left a faint (111) peak,and the Al (111)peak disappeared completely after 50 h of milling.Then,the diffraction peaks no longer changed with the extension of milling time,but the FWHM was slightly increased.In addition,the diffraction peak of W moved slightly to the higher angle during the MA process,which implies that Al atoms may dissolve into the W lattice and resulting in a decrease in the W lattice constant.Fig.3(b)illustrates the evolution trend of the grain size and lattice strain.The grain size of the W(Al60) mixed powder decreased monotonically with the extension of milling time,and the grain size decreased to 35 nm after 70 h ball milling.The reduction in the grain size will greatly increase the grain boundary area.The grain boundary is the bridge of atomic diffusion and the place where solute atoms gather first,the increase of the grain boundary area is beneficial to the expansion of solid solubility[38,47].Different from the grain size,the lattice strain first increased and then decreased with milling time.This means that amorphous phase may appear after ball milling for 40 h.Generally,amorphization will reduce interfaces and defects,which helps to release lattice strain [48].Further characterization of the crystal structure will be shown in HRTEM analysis.

    Fig.3.(a) XRD patterns of W(Al60) mixed powder with different milling times;(b) Evolution trend of grain size and lattice strain during ball milling.

    Fig.4(a)-Fig.4(c)provide the morphology evolution during the MA process.In the early stage,the powders undergo plastic deformation and become flakes (Fig.4(a)).Subsequently,the flake particles were cold welded and stacked together.Continuous plastic deformation led to severe work hardening.When the degree of work hardening exceeds a certain level,the particles begin to refine.Fig.4(b)and Fig.4(c)shows the morphology and EDS results of the powder milled for 50 h and 70 h,respectively.Considering the possible differences in the element ratio of small particle samples,we selected multiple small particles for EDS test to ensure the representativeness and repeatability of the date shown in Fig.4(b) and Fig.4(c).Although the diffraction peak of Al disappeared after ball milling for 50 h,the EDS results reveal that the element distribution was nonuniform.As shown in Fig.4(b),there was still a significant difference in the element ratio between region II and region III.However,when the milling time reached 70 h,the uniformity of the element distribution was obviously improved.In Fig.4(c),the W/Al atomic ratio in region I,II and III was close to 2:3,which is consistent with the design chemical ratio.

    Fig.4.SEM images of W(Al60)mixed powders with milling time of(a)10 h;(b)50 h and(c)70 h;(d)HAADF image and elemental maps of obtained W(Al60)alloy powder;(e)SAED;(f) HRTEM and (g) Corresponding IFFT image of obtained W(Al60) alloy powder.

    HRTEM observation was utilized to study the crystal structure of the alloy powder obtained by ball milling for 70 h.As shown in Fig.4(d),the HAADF image shows that the solvent and solute atoms have good dispersibility at the nanoscale,and the elemental maps of W and Al also demonstrated a good matching degree.Although the element distribution is uniform on the whole,it should be pointed out that the distribution of W still reflects certain local aggregation characteristics,indicating that the obtained W(Al60)alloy powder is not composed of a single-phase solid solution structure.Fig.4(e) gives the selected area electron diffraction(SAED)image of the prepared W(Al60)alloy powder.The diffraction rings correspond to different crystal plane groups of W,and the diffraction features of Al completely disappeared.In addition,the SAED image shows faint amorphous features.Through HRTEM image observation (Fig.4(f)),it can be confirmed that a small amount of amorphous phase existed around the crystal grains.Therefore,the prepared alloy powder with ball milling for 70 h was actually a mixture of solid solution and amorphous phase.The existence of amorphous phase also explains why the lattice strain in Fig.3(b) decreased after 40 h of ball milling.Under repeated impact and extrusion,the mixed powder experienced continuous deformation,and the lattice defect accumulated during this process.When the lattice defect exceeded the critical value,the lattice tended to collapse and transform into amorphous structure,and lattice strain in the crystalline region decreased.Besides,we also provide an inverse fast Fourier transform(IFFT)image to clarify the lattice defects in W(Al60)super-saturated solid solution powders.In the IFFT image,the lattice defects are reflected as the bending of the stripes.As shown in Fig.4(g),numerous lattice defects can be observed in the crystalline region.The lattice defects will accumulate with the extension of milling time;excessive lattice defects can lead to the formation of amorphous phase.Therefore,If the milling time is too long,the amorphous phase may gradually become the dominant phase of the alloy powder.

    An important prerequisite for the application of W(Al60) alloy powder is good thermal stability.Therefore,the mixed powders with different milling times were heated to 1000°C in an Ar atmosphere to test their stability at high temperature.As shown in Fig.5,for the original powder,there was an endothermic peak P1 and an exothermic peak P2,which corresponded to the melting of Al and the formation of Al4W,respectively.Extending the milling time to 30 h,the intensity of peak P1 was significantly weakened,while peak P2 disappeared completely.After 70 h of milling,there was no obvious endothermic or exothermic peak in the temperature range from room temperature to 1000°C,indicating that Al atoms were completely dissolved into the W lattice and could stay stable below 1000°C.The good thermal stability of W(Al60)supersaturated solid solution powder provides more options for the preparation of W(Al60)-based composites.

    Fig.5.DSC curves of W(Al60) mixed powder with milling time of 0 h,30 h and 70 h.

    3.2.Mechanical properties of the W(Al60)-Al composite

    The W(Al60)-Al composite based on W(Al60) super-saturated solid solution powder was fabricated by hot pressing.The density of W(Al60) alloy powder measured by pycnometer method is 9.30 g/cm3.According to the law of volume invariance,the theoretical maximum density (TMD) of the W(Al60)-Al composite is calculated to be 6.66 g/cm3.The actual density of the sample measured by Archimedes method was 6.44 g/cm3,which exceeded 96% of the TMD.XRD analysis shows that the W(Al60)-Al composite was a two-phase mixture of W(Al60) solid solution and Al matrix,and no intermetallic compounds were formed during the HP process (Fig.6(a)).The corresponding backscattered electron(BSE) image is given in Fig.6(b).It can be clearly seen that the W(Al60)-Al composite consists of two phases: the darker phase is the Al matrix,and the white phase is the W(Al) solid solution.Significantly,the fine W(Al60) particles exhibited obvious aggregation,and the split effect of W(Al60)particles on the Al matrix may affect the plasticity and strength of the W(Al60)-Al composite.

    Fig.6.(a)XRD pattern and(b)BSE image of the W(Al60)-Al composite;(c)True stress-strain curves and macroscopic cracks under quasi-static loading;(d)and(e)Microstructure of tensile fracture.

    Fig.6(c) shows the true stress-strain curves of the W(Al60)-Al composite under quasi-static loading.The quasi-static compressive and tensile strength of the composite are 320.4 MPa and 161.4 MPa,respectively.The macroscopic compression cracks propagate along the direction of 45°to the axial direction,indicating that the compression fracture mode is shear fracture.The cohesion strength of the weak interface between W(Al60) and Al matrix cannot effectively resist the shear load,which is the main reason for the lower compressive strength.According to Meyers's research on Al-Ta,Al-W,Al-Ni and other Al-based composite materials [12],the mechanical strength of composite materials with a “soft” matrix is generally low due to the lack of restriction on deformation and weak interfacial bonding strength.Similar conclusions were also obtained by Williamson and coworkers [49].The tensile properties of the materials are more sensitive to the interfacial bonding strength and original defects.As shown in Fig.6(d),numerous cracks occurred at the interface between the Al matrix and W(Al60).The separation of W(Al60) alloy powder from Al matrix is also the main inducer of tensile fracture.In addition,it was mentioned that some W(Al60) particles are gathered together by mechanical combination.The limited bonding strength makes these areas prone to cracks(Fig.6(e)).However,sufficient crushing is beneficial for the contact of the reaction components with air.The brittleness property may improve the sensitivity of the W(Al60)-Al composite under high speed impact [14].

    3.3.Impact initiated experiment results

    The impact initiated experiments were conducted through cylindrical projectiles against the test chamber.To study the response characteristics under different impact speeds,the tests were carried out at design speeds.According to the high-speed camera records,the impact induced initiation process of the W(Al60)-Al composite can be divided into four stages: perforate front plate,impact steel anvil,projectile fracture into debris clouds,chemical reaction and reaction product venting from the chamber.Fig.7(a)shows the reaction phenomenon of the W(Al60)-Al composite under an impact velocity of 1362 m/s.The violent flame vented from the chamber and the dazzling lights lasted for more than 200 ms.

    Loving Gordon doesn t mean in any way that we love you any less. There is plenty of room for two wonderful fathers in our lives. Gordon always encouraged contact with you, never spoke6 a word against you or undermined our feelings for you. We respect the fact that you never voiced negative feelings about Gordon.

    Fig.7.(a) Overpressure curves of W(Al60)-Al composite at different impact velocities;(b) Statistical diagram of fragment recovery rate and particle size distribution.

    Chemical reaction will produce an instantaneous blast wave and decline rapidly.The continuous reaction of remaining debris clouds continue to replenish the pressure and produce a global quasistatic pressure in the chamber.Generally,the overpressure value is used as a standard to measure the reaction characteristics of ESMs.Normally,the overpressure value of an inert projectile (20 steel) impacting the target plate at a speed of about 1200 m/s is about 0.03 MPa (Fig.7(a)).In contrast,the W(Al60)-Al composite exhibits a peak overpressure of 0.13 MPa at impact velocity of 988 m/s,indicating that it is a low-threshold ESMs.In addition,the overpressure values corresponding to impact velocities of 1181 m/s and 1362 m/s are 0.16 MPa and 0.28 MPa,respectively,showing a significant positive correlation with the impact velocity.This is mainly related to the fracture degree of the projectiles.In the impact initiated experiment,the projectile will break into debris clouds after impacting the steel anvil.Severe plastic deformation caused a rapid temperature rise.When small debris clouds with high temperature contact with air,deflagration will occur and release chemical energy.Therefore,the degree of fragmentation has a significant impact on the reaction characteristics.Fig.7(b)shows statistical data on the size of particles recovered after the reaction.The average size of the recovered particles was significantly smaller at an impact speed of 1362 m/s,and the particles below 150 μm accounted for more than 70%,indicating that the projectile was broken sufficiently.Sufficient fragmentation induced more violent chemical reactions.And the fine dust produced by combustion will vent from the test chamber with the gas.Therefore,the fragment recovery rate is negatively correlated with the impact velocity,and the fragment recovery rate at 1362 m/s is only 18.6%.

    3.4.Analysis of reaction products and mechanism

    To determine the reaction process of the W(Al60)-Al composite under high speed impact,the reaction products were analyzed by XRD and SEM.Fig.8 shows the XRD pattern of recovered fragments.Compared with the original XRD pattern of the W(Al60)-Al composite in Fig.6(a),the phase composition of the reaction product is significantly more complex,including not only W,Al2O3and Al4W,but also WO3,W3O and other atypical tungsten oxides WOx(e.g.WO2.9and W19O55).Among them,the appearance of Al2O3and Al4W phases can be expected,which is the typical reaction product of most W-Al composite materials.Surprisingly,the existence of a variety of tungsten oxides indicates that W is involved in the oxidation reaction.

    Fig.8.XRD patterns of reaction product.

    SEM observations provide direct evidence for this conclusion.Fig.9(a) shows the cross-sectional morphology of a larger recovered fragment.The dark region corresponds to the Al2O3phase formed by the Al matrix reaction.In addition,there are a large number of lath phases and dispersed spot phases.Fig.9(b) is the enlarged view of this area.EDS analysis suggests that the lath phase is Al4W,while the dispersed spots are W particles.The formation and distribution of Al4W indicate that Al has been precipitated during the impact reaction,and the dispersed fine W particles are evidence of Al precipitation.Because the size of this kind of particle is large and wrapped by the Al matrix,it cannot contact enough oxygen.Therefore,the precipitated Al reacted with W in situ to form Al4W.In addition,considering that W(Al60) super-saturated solid solution powder is stable below 1000°C,the precipitation of Al indicates that the impact initiation experiment may be accompanied by a high-temperature field.

    Fig.9(c) provides a representative morphology of the reaction product with a small size.Such regular spherical particles were formed by the solidification of molten liquid droplets,which confirms that a very high temperature was generated during the reaction process.According to the EDS results,the particle is mainly composed of three phases: the dark phase is Al2O3;the bright phase has a W/O element ratio of 2:3,so we think it is an insufficient combustion product of W;the light gray phase adjacent to Al2O3was also detected to be an oxide of tungsten,but the W/O element ratio of the bright phase is close to 1:3,which can be presumed to be the tungsten combustion product WO3(Fig.9(d)).Similar results can also be obtained in Fig.9(e)and Fig.9(f).Fig.9(e)shows an irregular large particle formed by the combustion and solidification of multiple small particles.A widely distributed WOXnetwork was found around the Al2O3phases.According to the EDS results,this phase may also be an insufficient combustion product of W.The above reaction product analyses prove that the W(Al60)super-saturated solid solution powder modified by MA has excellent reactivity.Driven by the impact kinetic energy and the combustion heat of the Al matrix,the W(Al60) particles successfully break the barrier of the high thermal inertia of W and deflagrate at the moment of contact air.In addition,the droplet-like reaction products indicate that the flame temperature is very high,which is conducive to the self-sustaining combustion of modified W(Al60)super-saturated solid solution powder.

    In traditional ESMs with W addition,the role of W is mainly reflected in the improvement of density and penetration ability,resulting in poor energy release characteristics of ESMs[50].For the W(Al60)-Al composite,the energy release rate is significantly improved due to the combustion of W.According to the study of Ames [51],the reaction products venting from the hole can be ignored since the time of reaction is short and the diameter of the hole on front plate is small.Thus,the chemical energy released in the chamber can be expressed as follows:

    where Δp is the peak value of quasi-static pressure,V is the chamber volume and is equal to 0.027 m3,and γ is the gas adiabatic index in the chamber,which is assumed to be 1.4.

    The energy release rate per unit mass iswheremis the mass of the projectile.Data on the energy release rate are tabulated in Table 1.The energy release rate corresponding to the impact velocity of 1362 m/s reaches 2.71 kJ/g,which is 65%of the energy released by TNT explosion (4.19 kJ/g).Besides,Fig.10 summarizes the reaction rate of some typical Al-base composites[52-55].It is worth mentioning that the density of the W(Al60)-Al composite is 6.44 g/cm3,which is higher than that of other Al-basedESMs shown in Fig.10.Generally,high density materials have better penetration characteristics.On this basis,the energy release rate of the W(Al60)-Al composite is still at a high level among all the Albased ESMs,even slightly better than that of the Al-Ni composites,which makes the W(Al60)-Al composite has greater competitive advantages in ESMs.

    Table 1 Energy release rate of W(Al60)-Al composite under different impact velocities.

    Fig.10.Energy release rate of typical Al-base composites.

    4.Conclusions

    This paper reports a W(Al60)-Al composite based on W(Al60)super-saturated solid solution powder.The density of the material has reached 6.44 g/cm3,which is higher than most Al-based ESMs.The impact initiation experiments show that the W(Al60)-Al composite is a low threshold ESMs and has excellent energy release characteristics.The energy release rate reaches 2.71 kJ/g at an impact velocity of 1362 m/s,which is even better than that of Al/Ni ESMs.In the reaction product analysis,a variety of tungsten oxides were detected,proving that the modified W based on supersaturated solid solution is combustible under impact load.The high activity and combustion heat of Al help W overcome the barrier of high thermal inertia,and effectively make up for the low adiabatic flame temperature of W.The droplet-like reaction products imply that the combustion process of W(Al60)alloy powder is accompanied by a high temperature field,which is of great significance for the self-sustaining combustion of W.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China,[Award number: 11972372] and [Award number: U20A20231].

    美女高潮到喷水免费观看| 国产亚洲av高清不卡| ponron亚洲| 1024香蕉在线观看| 高清毛片免费观看视频网站 | 精品福利观看| 韩国精品一区二区三区| 淫秽高清视频在线观看| 18禁黄网站禁片午夜丰满| 精品欧美一区二区三区在线| 无限看片的www在线观看| 免费在线观看影片大全网站| 欧美激情极品国产一区二区三区| 一级作爱视频免费观看| 亚洲色图综合在线观看| 18禁观看日本| 国产成年人精品一区二区 | 一进一出抽搐gif免费好疼 | 99久久99久久久精品蜜桃| 久久婷婷成人综合色麻豆| 青草久久国产| 女人精品久久久久毛片| 欧美人与性动交α欧美软件| 在线观看一区二区三区| 亚洲国产精品一区二区三区在线| 成人三级做爰电影| 国产免费男女视频| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲| 黄色片一级片一级黄色片| 国产在线精品亚洲第一网站| 精品国产超薄肉色丝袜足j| av在线播放免费不卡| 久久人妻福利社区极品人妻图片| 精品一区二区三区四区五区乱码| 久久久久久久久中文| 麻豆久久精品国产亚洲av | 一级毛片高清免费大全| 亚洲精品在线观看二区| 国内毛片毛片毛片毛片毛片| 国产精品亚洲av一区麻豆| 亚洲av片天天在线观看| 午夜影院日韩av| 精品久久久久久久毛片微露脸| 亚洲一卡2卡3卡4卡5卡精品中文| av国产精品久久久久影院| 亚洲 欧美 日韩 在线 免费| www.999成人在线观看| 超色免费av| 国产高清videossex| 欧美日韩精品网址| 亚洲视频免费观看视频| 91国产中文字幕| 成人影院久久| 在线观看www视频免费| 久久天堂一区二区三区四区| 欧美乱妇无乱码| 国产日韩一区二区三区精品不卡| 黄色 视频免费看| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 色播在线永久视频| 在线视频色国产色| 亚洲aⅴ乱码一区二区在线播放 | 色综合站精品国产| av片东京热男人的天堂| 欧美精品啪啪一区二区三区| 日本 av在线| cao死你这个sao货| av天堂久久9| 在线观看66精品国产| 婷婷六月久久综合丁香| 久久久久九九精品影院| 国产99久久九九免费精品| 亚洲精品国产区一区二| 久久精品人人爽人人爽视色| 成人国语在线视频| 欧美不卡视频在线免费观看 | 欧美大码av| 又紧又爽又黄一区二区| 亚洲精品美女久久av网站| 欧美成狂野欧美在线观看| 亚洲情色 制服丝袜| 午夜两性在线视频| 美女高潮到喷水免费观看| 在线观看免费午夜福利视频| 99久久99久久久精品蜜桃| 婷婷精品国产亚洲av在线| 精品久久久久久久毛片微露脸| 国产国语露脸激情在线看| 一级毛片精品| 亚洲少妇的诱惑av| 在线看a的网站| 在线av久久热| 亚洲自偷自拍图片 自拍| 欧美av亚洲av综合av国产av| 一个人观看的视频www高清免费观看 | 国产欧美日韩精品亚洲av| 午夜精品在线福利| 一级毛片高清免费大全| 丁香欧美五月| 这个男人来自地球电影免费观看| av天堂久久9| 久久久久久免费高清国产稀缺| 免费在线观看完整版高清| 一区二区三区精品91| 久久国产精品人妻蜜桃| 日日摸夜夜添夜夜添小说| 亚洲 国产 在线| 亚洲av电影在线进入| 极品人妻少妇av视频| 久久精品国产综合久久久| 女人高潮潮喷娇喘18禁视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成77777在线视频| 热99re8久久精品国产| 国产亚洲精品综合一区在线观看 | 欧美中文综合在线视频| 精品国产亚洲在线| 久久久国产精品麻豆| 国产国语露脸激情在线看| 丝袜美足系列| 一a级毛片在线观看| 国产精品 国内视频| 国产高清videossex| 日韩中文字幕欧美一区二区| 五月开心婷婷网| 久久这里只有精品19| 欧美乱码精品一区二区三区| 在线观看免费视频网站a站| 亚洲五月天丁香| 搡老乐熟女国产| 国产一区二区三区视频了| 欧美日本中文国产一区发布| 亚洲免费av在线视频| 超碰成人久久| 1024香蕉在线观看| 国产精品秋霞免费鲁丝片| 亚洲av片天天在线观看| 国产高清激情床上av| 日韩大码丰满熟妇| 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月| 在线观看午夜福利视频| 热re99久久国产66热| 国产男靠女视频免费网站| 免费女性裸体啪啪无遮挡网站| 精品午夜福利视频在线观看一区| 不卡一级毛片| 在线看a的网站| 久久久久精品国产欧美久久久| 在线观看www视频免费| 五月开心婷婷网| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 欧美日韩中文字幕国产精品一区二区三区 | cao死你这个sao货| 精品国内亚洲2022精品成人| 欧美av亚洲av综合av国产av| 青草久久国产| 两性午夜刺激爽爽歪歪视频在线观看 | 一区在线观看完整版| 久久国产亚洲av麻豆专区| 国产三级在线视频| 亚洲五月婷婷丁香| 国产成+人综合+亚洲专区| 国产精品久久久人人做人人爽| 91国产中文字幕| 黄色毛片三级朝国网站| 久久久国产欧美日韩av| 日韩高清综合在线| 韩国av一区二区三区四区| 久久国产精品男人的天堂亚洲| 午夜福利在线观看吧| 国产精品一区二区免费欧美| 涩涩av久久男人的天堂| 午夜免费激情av| 国产精品野战在线观看 | 国产精品 欧美亚洲| 亚洲人成网站在线播放欧美日韩| 一二三四社区在线视频社区8| 视频区图区小说| 亚洲午夜精品一区,二区,三区| 宅男免费午夜| 亚洲免费av在线视频| 老司机在亚洲福利影院| 国产免费现黄频在线看| 丰满人妻熟妇乱又伦精品不卡| 久久久久久免费高清国产稀缺| 视频区欧美日本亚洲| 91麻豆av在线| 久久婷婷成人综合色麻豆| 黑丝袜美女国产一区| 在线观看66精品国产| 国产高清激情床上av| 久久午夜亚洲精品久久| 韩国av一区二区三区四区| 一区福利在线观看| 热99re8久久精品国产| 老司机亚洲免费影院| 久久天堂一区二区三区四区| 国产97色在线日韩免费| 欧美日韩黄片免| 天天添夜夜摸| 久久影院123| 亚洲专区中文字幕在线| 亚洲精华国产精华精| 国产激情欧美一区二区| 日韩欧美一区视频在线观看| 国产欧美日韩综合在线一区二区| 在线观看免费视频日本深夜| 久久人妻福利社区极品人妻图片| 怎么达到女性高潮| 男女下面进入的视频免费午夜 | 88av欧美| 欧美成人免费av一区二区三区| 亚洲成人免费电影在线观看| 黄频高清免费视频| 一级毛片女人18水好多| 在线观看免费视频网站a站| 波多野结衣一区麻豆| 国产亚洲精品一区二区www| 日韩欧美一区视频在线观看| 国产亚洲欧美98| 麻豆成人av在线观看| 久久久国产精品麻豆| 精品久久久久久成人av| 神马国产精品三级电影在线观看 | 妹子高潮喷水视频| 大型av网站在线播放| 精品国内亚洲2022精品成人| 国产精品久久久av美女十八| 国产精品美女特级片免费视频播放器 | 女性被躁到高潮视频| 岛国视频午夜一区免费看| 日本vs欧美在线观看视频| 国产一区二区三区综合在线观看| 首页视频小说图片口味搜索| 亚洲精品成人av观看孕妇| 美女高潮到喷水免费观看| 欧美中文综合在线视频| 自线自在国产av| av超薄肉色丝袜交足视频| 国产成人免费无遮挡视频| 在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 麻豆成人av在线观看| 国产一区二区三区在线臀色熟女 | 国产黄a三级三级三级人| 精品欧美一区二区三区在线| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 亚洲色图综合在线观看| 国产一卡二卡三卡精品| 老司机福利观看| 午夜精品国产一区二区电影| 超色免费av| av电影中文网址| 欧美日本亚洲视频在线播放| 欧美日韩一级在线毛片| 在线观看免费视频日本深夜| 亚洲av第一区精品v没综合| 久久青草综合色| 男女下面进入的视频免费午夜 | 99热只有精品国产| 成熟少妇高潮喷水视频| 97碰自拍视频| 麻豆国产av国片精品| 一边摸一边抽搐一进一出视频| 国产三级黄色录像| 久久热在线av| 一区二区日韩欧美中文字幕| 在线av久久热| 嫩草影视91久久| a级毛片在线看网站| 看黄色毛片网站| 丰满的人妻完整版| 少妇的丰满在线观看| 麻豆久久精品国产亚洲av | 国产无遮挡羞羞视频在线观看| 欧美激情高清一区二区三区| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久 | 69av精品久久久久久| 美女福利国产在线| 国产精品一区二区三区四区久久 | 免费搜索国产男女视频| 国产免费男女视频| 黑丝袜美女国产一区| 黑人欧美特级aaaaaa片| 午夜a级毛片| aaaaa片日本免费| 亚洲熟女毛片儿| 一夜夜www| 欧美最黄视频在线播放免费 | 国产精品亚洲av一区麻豆| 丝袜在线中文字幕| 亚洲人成电影观看| 精品福利观看| 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 黄片播放在线免费| 久热这里只有精品99| 99久久精品国产亚洲精品| 久久精品91蜜桃| 一个人观看的视频www高清免费观看 | 久久久久精品国产欧美久久久| 国产野战对白在线观看| 日本wwww免费看| 变态另类成人亚洲欧美熟女 | 在线av久久热| 99国产极品粉嫩在线观看| 午夜精品在线福利| 国产乱人伦免费视频| 国产免费av片在线观看野外av| 国产熟女午夜一区二区三区| 欧美日韩精品网址| 亚洲精品美女久久av网站| 老司机午夜福利在线观看视频| 午夜精品国产一区二区电影| 欧美色视频一区免费| 亚洲精品一卡2卡三卡4卡5卡| 国产有黄有色有爽视频| 久久精品国产综合久久久| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 80岁老熟妇乱子伦牲交| 国产激情久久老熟女| 亚洲成人免费av在线播放| 99精品久久久久人妻精品| 中文亚洲av片在线观看爽| cao死你这个sao货| 一边摸一边抽搐一进一出视频| 亚洲九九香蕉| 欧美乱码精品一区二区三区| 黄色视频,在线免费观看| 国产成人av激情在线播放| 久热这里只有精品99| 可以免费在线观看a视频的电影网站| 日韩免费av在线播放| 日韩高清综合在线| 国产亚洲精品久久久久久毛片| 最近最新中文字幕大全免费视频| 国产免费男女视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩国产mv在线观看视频| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 在线观看免费视频网站a站| 亚洲第一欧美日韩一区二区三区| 日韩成人在线观看一区二区三区| 国产区一区二久久| 欧美日本亚洲视频在线播放| 国产亚洲欧美98| 男人舔女人的私密视频| 操出白浆在线播放| 三上悠亚av全集在线观看| 亚洲五月婷婷丁香| 日本黄色视频三级网站网址| 在线观看免费日韩欧美大片| 夜夜看夜夜爽夜夜摸 | 老司机深夜福利视频在线观看| 大陆偷拍与自拍| 国产激情欧美一区二区| av有码第一页| 99国产精品99久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 成人国语在线视频| 亚洲美女黄片视频| 国产区一区二久久| 天天躁狠狠躁夜夜躁狠狠躁| 在线av久久热| 成人三级黄色视频| 老司机午夜福利在线观看视频| 91成人精品电影| 无人区码免费观看不卡| 黑人巨大精品欧美一区二区蜜桃| 日本vs欧美在线观看视频| 亚洲 国产 在线| 在线观看免费午夜福利视频| 日韩精品青青久久久久久| 99香蕉大伊视频| 久久婷婷成人综合色麻豆| 青草久久国产| 女人被躁到高潮嗷嗷叫费观| 在线观看日韩欧美| 国产欧美日韩一区二区三区在线| 一级作爱视频免费观看| 纯流量卡能插随身wifi吗| 真人一进一出gif抽搐免费| 日韩中文字幕欧美一区二区| 后天国语完整版免费观看| 天堂俺去俺来也www色官网| 免费看十八禁软件| 亚洲精品一卡2卡三卡4卡5卡| 亚洲黑人精品在线| 国产精品免费视频内射| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区视频在线观看免费 | 中出人妻视频一区二区| 99re在线观看精品视频| 法律面前人人平等表现在哪些方面| 美女午夜性视频免费| 91麻豆av在线| 午夜a级毛片| 法律面前人人平等表现在哪些方面| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频| 男女做爰动态图高潮gif福利片 | 一级毛片高清免费大全| 国产精品爽爽va在线观看网站 | 97超级碰碰碰精品色视频在线观看| 成人国产一区最新在线观看| 精品久久久久久电影网| 亚洲少妇的诱惑av| 中文字幕色久视频| av欧美777| √禁漫天堂资源中文www| 中文字幕人妻熟女乱码| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 悠悠久久av| 在线观看日韩欧美| 99久久人妻综合| 亚洲人成伊人成综合网2020| 亚洲成人免费av在线播放| 亚洲成国产人片在线观看| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人国产一区在线观看| xxxhd国产人妻xxx| 亚洲av日韩精品久久久久久密| 丝袜美足系列| 51午夜福利影视在线观看| 岛国在线观看网站| 日韩高清综合在线| 麻豆久久精品国产亚洲av | 巨乳人妻的诱惑在线观看| 欧美成人免费av一区二区三区| 国产精品香港三级国产av潘金莲| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 中文字幕高清在线视频| 亚洲九九香蕉| 97人妻天天添夜夜摸| 又大又爽又粗| 久久婷婷成人综合色麻豆| 久99久视频精品免费| 久久欧美精品欧美久久欧美| 女人被狂操c到高潮| 久久国产精品男人的天堂亚洲| 校园春色视频在线观看| av视频免费观看在线观看| 最新美女视频免费是黄的| 日本黄色日本黄色录像| 9191精品国产免费久久| 亚洲av熟女| 变态另类成人亚洲欧美熟女 | 国产99久久九九免费精品| 久久精品人人爽人人爽视色| 在线十欧美十亚洲十日本专区| 日韩人妻精品一区2区三区| 丰满饥渴人妻一区二区三| 国产极品粉嫩免费观看在线| 国产野战对白在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲久久久国产精品| 91精品国产国语对白视频| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 欧美日韩黄片免| 精品国产一区二区久久| 99re在线观看精品视频| 国产成人av激情在线播放| 成人亚洲精品一区在线观看| 99国产精品一区二区三区| 久久精品91无色码中文字幕| 欧美日韩亚洲高清精品| 97超级碰碰碰精品色视频在线观看| 美女午夜性视频免费| 欧美大码av| 国产一区二区激情短视频| 国产熟女xx| 夜夜爽天天搞| 欧美亚洲日本最大视频资源| 夜夜夜夜夜久久久久| 亚洲精品久久成人aⅴ小说| 亚洲精品久久午夜乱码| 国内毛片毛片毛片毛片毛片| 一级a爱视频在线免费观看| 波多野结衣一区麻豆| av天堂在线播放| av网站在线播放免费| 80岁老熟妇乱子伦牲交| 久久青草综合色| 亚洲性夜色夜夜综合| 欧美日韩精品网址| 亚洲色图 男人天堂 中文字幕| 18禁观看日本| 一级片'在线观看视频| 日韩欧美三级三区| 男男h啪啪无遮挡| 午夜亚洲福利在线播放| 女同久久另类99精品国产91| 国产欧美日韩综合在线一区二区| 国产在线观看jvid| 成人18禁高潮啪啪吃奶动态图| 长腿黑丝高跟| 一本综合久久免费| 一边摸一边抽搐一进一小说| 精品电影一区二区在线| 国产成人欧美在线观看| 日本 av在线| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 日本黄色视频三级网站网址| 女人被狂操c到高潮| 欧美激情 高清一区二区三区| 91在线观看av| a在线观看视频网站| 日本五十路高清| 亚洲欧洲精品一区二区精品久久久| 久久精品人人爽人人爽视色| 一级a爱视频在线免费观看| 视频区图区小说| 久久亚洲精品不卡| 看片在线看免费视频| 久久久久久久精品吃奶| 久久久久久人人人人人| 精品久久久久久成人av| 琪琪午夜伦伦电影理论片6080| 美女国产高潮福利片在线看| 国产精品美女特级片免费视频播放器 | 午夜免费观看网址| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 亚洲少妇的诱惑av| 最新美女视频免费是黄的| 丝袜美腿诱惑在线| 日韩中文字幕欧美一区二区| 久久久国产成人精品二区 | 亚洲专区字幕在线| 青草久久国产| 国产又爽黄色视频| 亚洲 欧美一区二区三区| 69av精品久久久久久| 日本免费一区二区三区高清不卡 | 黑人巨大精品欧美一区二区mp4| 精品久久久久久电影网| 久久中文看片网| 一本大道久久a久久精品| 最近最新免费中文字幕在线| 亚洲av成人不卡在线观看播放网| 在线播放国产精品三级| 极品教师在线免费播放| 国产99白浆流出| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品美女久久久久99蜜臀| 欧美 亚洲 国产 日韩一| 两性夫妻黄色片| 欧美在线黄色| 99国产精品99久久久久| 这个男人来自地球电影免费观看| 91成人精品电影| 人人澡人人妻人| 国产亚洲精品久久久久5区| 欧美不卡视频在线免费观看 | 自线自在国产av| 国产精品 欧美亚洲| 在线永久观看黄色视频| 久久精品国产综合久久久| 亚洲九九香蕉| 麻豆久久精品国产亚洲av | 老司机靠b影院| 中文欧美无线码| 精品一区二区三卡| 精品高清国产在线一区| 女性生殖器流出的白浆| 一级黄色大片毛片| 制服诱惑二区| 日本 av在线| 欧美久久黑人一区二区| 国产伦人伦偷精品视频| 国产精品久久久久成人av| 99精品久久久久人妻精品| 国产蜜桃级精品一区二区三区| 丰满迷人的少妇在线观看| 18禁美女被吸乳视频| 老司机在亚洲福利影院| av片东京热男人的天堂| 亚洲欧美日韩无卡精品| 可以在线观看毛片的网站| 免费看十八禁软件| 日韩中文字幕欧美一区二区| 老司机福利观看| 色综合站精品国产| 好看av亚洲va欧美ⅴa在| 99在线视频只有这里精品首页| 自线自在国产av| 黑人欧美特级aaaaaa片| 交换朋友夫妻互换小说| 变态另类成人亚洲欧美熟女 | www.精华液| 80岁老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频| 欧美成人性av电影在线观看| 波多野结衣高清无吗| 女人精品久久久久毛片| 19禁男女啪啪无遮挡网站| 在线观看免费视频日本深夜| 亚洲人成电影免费在线| a级片在线免费高清观看视频| 一级毛片女人18水好多|