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

    Induced dipole dominant giant electrorheological fluid

    2023-09-05 08:48:36RongShen沈容KunquanLu陸坤權ZhaohuiQiu邱昭暉andXiaominXiong熊小敏
    Chinese Physics B 2023年7期

    Rong Shen(沈容), Kunquan Lu(陸坤權),?, Zhaohui Qiu(邱昭暉), and Xiaomin Xiong(熊小敏)

    1Beijing National Laboratory for Condensed Matter Physics,Key Laboratory of Soft Matter and Biological Physics,Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

    2School of Physics,Sun Yat-sen University,Guangzhou 510275,China

    Keywords: electrorheological fluid,induced dipole,vacancies,high energy ball milling

    1.Introduction

    In 1949, Winslow[1]reported that the shear stress of a suspension composed of solid particles and insulating liquid would increase with the applied electric field, which was called the electrorheological (ER) effect.Because of the potential application prospect of this effect, great efforts have been made for the decades to obtain ER materials that can be used in practice.To date, however, ER fluids have still not been widely used in technology and engineering,due to no ER fluids with good enough comprehensive performance suitable for applications.

    In the research history of ER fluids,it has mainly experienced the developments of two generations of ER fluid materials.The first generation,termed as dielectric ER fluids,was prepared by mixing solid dielectric particles with insulating oil, of which the effect is based on the interactions of dielectric particles polarized in the electric field.[2–5]Thus the yield stressτymust be proportional to the square of electric fieldE,i.e.,τy∝E2.Despite extensive tests performed,the yield stress of the prepared ER fluids was only reach few kPa.In 1996,Maet al.made a theoretical calculation using the first principles,and concluded that the upper bound of the yield stress of this kind of ER fluids could only reach about 10 kPa.[6]The mechanism of traditional ER fluid determines that its yield stress is low,which is unable to meet the applications requirements.In 1998,it was found that the ER fluid prepared by surface modified Sr–Ti–O particles possessed a high yield stress,reaching about 30 kPa at an electric field of 3 kV/mm.[7,8]In 2003,an ER fluid consisted of Ba–Ti–O nano-particles with urea polar molecules coated prepared by Wenet al.had a yield stress of over 100 kPa, which was called giant ER(GER) fluid.[9]For GER materials the solid particles are generally prepared by wet chemical methods,and their surfaces are coated with polar molecules.[10–14]As an electric fieldEis applied,the polarized particles attract and get close to each other,and the local electric fieldElocbetween the particles can be 102–103times ofE.[15–19]Such a highEloccan cause the polar molecules on the particle surface in the region to orientate alongElocand interact with the polarized units on the surfaces of adjacent particles.These interactions are much larger than that between of polarized particles as in dielectric ER fluid,thus results in high shear stress up to several hundred kPa in GER.Because of the effect of ER originated from the orientation of inherent dipole moment the relationship between yield stress and electric field should be linear,that is,τy∝E,which is different from that ofτy∝E2as in conventional ER fluid.This GER is also known as polar molecule dominant ER fluids (PM-ER),[11,12]which can be considered as the second generation ER fluid with high yield stress.For the PM-ER fluids, however, the fatal flaw is that the coating or modification layer on the particle surface are easy to be worn and the polar molecules fail, resulting in the continuous decline of its shear stress with wearing time.Although the wear resistance can be slightly improved by dissolving polar molecules into the liquid,[20,21]it is still not suitable for application because polar molecules will lose their polarity due to the actions such as space charge.

    It can be seen from the above that the existing firstgeneration of the dielectric ER fluid and second-generation PM-ER fluid have not been widely used in practice so far due to their inherent deficiency.New ER fluids that can be applied in technology and engineering is eagerly anticipated.

    In the study, we found that the ER fluids prepared with two kinds of modified dielectric particles provide excellent ER effect: One is to introduce a large number of vacancies into the dielectric particles;The other is to embed conductive microclusters in the dielectric particles.Due to the action of the high local fieldElocat particle spacing, the vacancies or the conductive microclusters on particle surfaces can generate large induced dipole moment and cause strong interactions between particles,resulting in the high shear stress of the ER fluids.This new ER fluid can be referred to as an induced dipole-dominanted ER(ID-ER)fluid.

    Present paper reports the preparation method as well as the action principle for ID-ER fluids,which are different from that of the dielectric ER and PM-ER fluids.The particles in ID-ER fluids containing vacancies or conductive microclusters both are prepared by high-energy ball milling (HEBM)technique.The raw materials before HEBM are either pure dielectric particles or dielectric particles mixed with a small amount of conductive material.The obtained ID-ER fluids exhibit attractive comprehensive performance: high yield stress(>100 kPa), low current density (<20 μA/cm2), nice temperature stability and wear resistance,settlement resistance as well.Furthermore,the preparation method is simple,good repeatability, low cost.Therefore, ID-ER fluids should be the new type of ER fluids suitable for applications.[22]

    We have studied the TiO2particles containing oxygen vacancies treated by HEBM method as well as the properties of prepared ID-ER fluids with these particles, and will introduce the results in detail.The ID-ER fluids prepared with some other dielectric particles besides TiO2containing oxygen vacancies are briefly described as well.As for the IDER fluid consisted of the particles embedded with conductive microclusters,we have discussed the method and principle of preparing embedded carbon microclusters by heat treatment in previous article,[23]and the results for the particles processed by HEBM are supplemented briefly in this paper.

    This article is presented in the following order: firstly,the preparation method of particles used in ID-ER fluid; the characterization of morphology, structure, and properties of the particles modified by HNBM; the property test of ID-ER fluid;finally,the proposed physical mechanism for explaining the ID-ER effect.In addition, the deficiency shortcoming of ID-ER fluid is pointed out.

    2.Experimental details

    Market-purchased oxide dielectric particles,such as TiO2powder, were used as the raw material and heated for several hours in a furnace above 400?C to remove any possible adsorbents from the particles.The powder was then treated with HEBM and milled for more than 10 hours.On the other hand,by mixing small amount of graphite or metal powder to dielectric particles to perform HEBM,the particles with conductive microclusters embedded can be achieved.The longer ball milling time is beneficial to get sufficient oxygen vacancies for the prepared particles, or the particles with better conductive microcluster dispersion.

    Particle density was measured by a true density instrument, x-ray diffraction (XRD), scanning electron microscopy(SEM),high resolution transmission electron microscope(HRTEM),electron paramagnetic resonance(EPR),and x-ray absorption fine structure (XAFS) technique were used to study and compare the morphologies and microstructures of particles before and after HEBM.Thus the oxygen vacancies and conductive microclusters in the particles induced by HEBM were detected and analyzed.XAFS was performed at Beijing Synchrotron Radiation Facility(BSRF).

    The readied particles through HEBM were mixed with silicone oil or machine oil to prepare ID-ER fluid.Shear stress were measured by self-made flat plate rheometer or Anton Paar rheometer.The flat plate rheometer with the plates of rough surfaces can measure ER samples with shear stress up to more than 200 kPa.The yield stress was measured at a shear rate of 0.2 s?1.

    Response time and frequency dependence measurements for the ID-ER fluids at different shear rates were performed by using a home-made rotary rheometer connected to Trek 10/40A type frequency and waveform adjustable high voltage power supply.The change rate of the rising and falling edges of the square wave was 750 V/μs and the data acquisition rate was 1000 s?1or 4000 s?1.

    For testing the abrasion resistance of ER fluid,the sample was sealed in a rotating container composed of an inner roller with diameter of 40.5 mm and an outer cylinder with diameter of 42 mm,which rotated for a long period with a shear rate of 300 s?1without electric field applied.During the period the ID-ER samples were taken out several times to measure the changes of yield stress of them atEfield due to the abrasion.

    We mainly studied anatase type TiO2nano-particles treated by HEBM technique as well as the prepared ID-ER fluid.Other powders with high dielectric constant such as rutile-type TiO2, BaTiO3, CaTiO3, CaCu3Ti4O12as well as SiO2powders with lower dielectric constant were also processed by HEBM and the properties of ID-ER fluids prepared with these particles were investigated, respectively.In addition, several ID-ER fluids consisted of the particles with graphite or metal microclusters embedded were tested.

    3.Results and discussion

    3.1.Microstructure changes and oxygen vacancies in particles induced by HEBM

    HEBM or saying mechanical alloying (MA) technique was developed in the 1970s to prepare composite materials.In the milling process, a large number of defects and microstructures are generated through intense mechanical actions such as collision, extrusion, shear and friction, which make the material undergo non-equilibrium transformation.Many new materials which are difficult to be obtained by other conventional methods can be synthesized.[24–27]HEBM studies on TiO2nanopowders have shown that anatase TiO2can be transformed into rutile structure,[28,29]and a large number of oxygen vacancies can be created.[30,31]The simulation results show that the vacancies increase with the milling time and gradually reach saturation.[32]

    In our investigation,the initial particles were anatase-type TiO2with a diameter of about 20 nm–30 nm.The comparisons of observed results for TiO2powders before and after HEBM are given below.Figures 1 and 2 show STEM, HRTEM images,and x-ray diffraction spectra of TiO2samples before and after milling respectively.It can be seen that HEBM can cause the TiO2powder transformed from anatase structure to rutile phase dominated structure with a larger disorder and a small size about 10 nm–15 nm.The density measurements show that the initial density of anatase-TiO2particles is 3.85 g/cm3and then the density of rutile-TiO2after HEBM is changed to be about 3.55 g/cm3,which is decreased about 15%compared with 4.20 g/cm3for initial rutile-TiO2state.

    Figures 3(a) and 3(b) show that both O 1s and Ti 2p of XPS binding energy spectra for rutile-TiO2particles move towards high energy after HEBM process.The increase of O 1s binding energy is considered to be the contribution of oxygen vacancy, while the increase of that for Ti 2p is attributed to structural change and defect accumulation.[30]That literatures reported Ti3+appeared in TiO2particles after treatment of HEBM,[30,31]is not observed by us.This may be due to the small size of TiO2nanoparticls,and Ti3+is easy to be oxidized into Ti4+in the air.

    Although many reports have mentioned that HEBM can produce a large number of oxygen vacancies in particles,[28–32]it is still a challenge to observe oxygen vacancies directly and detect their concentration.Electron paramagnetic resonance (EPR) detection was performed on a sample to determine whether oxygen vacancies appeared in TiO2particles after HEBM.Figure 3(c)compares the change ofgfactor of EPR for TiO2particles before and after HEBM process.g=2.003 can be observed in milled TiO2particles,which is EPR characteristic structure of oxygen vacancy capturing an electron,whileg=2.000 for unmilled TiO2particles, which may be the contribution of free radicals or free electrons.[30,33]

    Fig.1.STEM(top)and HRTEM(bottom)images of TiO2 nanoparticles before(left)and after(right)HEBM.

    Fig.2.XRD patterns of anatase structure TiO2 powder before (black)and after(red)HEBM.The upper and lower are rutile and anatase TiO2 powder diffraction standards,respectively.

    The XAFS technique can study the atomic neighbor structure and is an effective method for studying oxygen vacancies quantitatively by analyzing the change of coordination number.We have applied the XAFS method to study the near neighbor structure of Ti atoms in TiO2particles before and after HEBM and then the oxygen vacancy content has been obtained by analyzing the change of coordination number of O atoms surround Ti atoms.For example, after 52 hours of HEBM, the average coordination number of Ti–O atoms decreases from 6 to 4.8, and then we can know that oxygen vacancies account for about 20% of the total O atoms.The investigations of XAFS for TiO2particles by HEBM under different conditions show that the oxygen vacancy accounts for about 10%–20%and even more of the total O atom number depending on the period of HEBM.Detailed results of the XAFS analysis will be reported in a separate paper.

    All these indicate that HEBM causes a large number of defects such as oxygen vacancies created in TiO2nanoparticles, of which the structure is changed from anatase to rutile.TiO2is a kind of non-stoichiometric compound expressed as TiO2?XwithXin a wider range.[33–37]According to our XAFS analysis, oxygen vacancies account for about 10%–20% of the total oxygen atoms.The density measurements show that the density of TiO2particles decreased by 15%after HEBM, which should be caused by a large number of grain boundary defects and oxygen vacancies in the particles.The existence of these oxygen vacancies in TiO2particles is the key factor to produce ID-ER effect.

    Fig.3.XPS spectra of O 1s(a)and Ti 2p(b)of rutile-TiO2 particles before and after HEBM.(c)The change of g factor of EPR for TiO2 particles before and after HEBM,the unit 1 Gs=10?4 T.

    3.2.Yield stress of ID-ER fluids

    The ER fluids consisted of the particles treated by HEBM to contain oxygen vacancies or conductive microclusters exhibit strong ER effect with yield stress up to several hundred kPa.In turn,the yield stress of the ER fluid made of unmilled TiO2particles is very low reaching only few kPa even for a sample with high volume fraction.In the following,the yield stress of two ID-GER fluids, vacancy type, and conducting microcluster type,varying with the electric field are given,respectively,while the latter is only described briefly.

    (i)Vacancy-type ID-ER fluid

    Figure 4(a) shows the relationship of yield stressτyand current densityρvarying with electric fieldEof ID-ER fluid,which was fabricated by mixing the TiO2particles containing oxygen vacancies to 32#machine oil with a volume fraction ofφ=43%.Because the yield stress was not very high,the slip effect of plate surface in the measurement can be negligible and then the measured values should truly reflect the relation betweenτyandE.In this caseτy–Ecurve can be well fitted byτy=a+bE2,showing the relationτy∝E2.In contrast,for an ID-ER fluid withφ=54%shown in the inset of Fig.4(a),the yield stress is high and some skid on the electrode plate was observed during measuring.The slippage results in the measured yield stress being lower than the actual value in the region of higher yield stress values deviating apparently fromτy∝E2relation.The principle of the relation ofτy versus Ewill be discussed later.The relationship between the yield stressτyand the volume fractionφof ID-ER fluids made by mixing the TiO2particles with 32#machine oil and 10#silicon oil are shown in Fig.4(b)respectively.

    Fig.4.(a)The dependence of τy and ρ on E of ID-ER fluid with φ=43%,which can be fitted by τy=a+bE2 with a=0.04 and b=0.96.Inset is the τy–E curve for φ=54%.(b)The relation of τy versus φ for the ID-ER fluids consisted of HEBM milled TiO2 nanoparticles and 10# silicone oil or 32# machine oil at 3 kV/mm and 5 kV/mm respectively.Inset is the photo of ID-ER sample with φ =54%.

    As seen in Fig.4(b)thatτyof ID-ER fluid can reach more than 200 kPa as the volume fraction is higher,and the sample still shows good fluidity.The leakage current density as low as less than 20 μA/cm2under a field of 5 kV/mm either by using silicon oil or machine oil.When preparing with silicone oil, because of its good temperature stability, the yield stress is almost unchanged in the range of 0?C–100?C,while it will change about 10% in machine oil.In addition, these ER fluids have good settlement resistance,and no obvious settlement has been observed in the sample placed for several months.

    The yield stresses of the ER fluids prepared with HEBM treated rutile-TiO2, BaTiO3, CaTiO3, CaCu3Ti4O12, LiB3O5,Al2O3,and SiO2particles are shown in Fig.5 respectively.All of them show good ER behavior except for SiO2particles.The low yield stress of ER fluid containing milled SiO2particles is due to the low dielectric constant of SiO2causing the inability for forming a high local electric field between particles.

    Fig.5.Relationship between the yield stress and electric field of ID-ER fluids prepared respectively by rutile-TiO2, BaTiO3, CaTiO3,CaCu3Ti4O12(CCTO),LiB3O5,Al2O3,SiO2 particles treated by HEBM.

    (ii)Conductive microcluster type ID-ER fluid

    In our previous report, TiO2particles inlaid with carbon microclusters were obtained by heating the Ti–O powder, which were originally made from wet chemistry, and the prepared ER fluid with these TiO2particles showed giant ER effect.[23]It is more convenient to prepare dielectric particles inlaid with conductive microclusters by HEBM technique,especially to disperse metal microclusters in the dielectric particles.For example, by mixing small amount graphite(0.67 wt%),Au(3 wt%),Ni(2 wt%)or stainless steel(6 wt%)powders to anatase-TiO2particles respectively to process by HEBM,the ER fluids composed with those particles perform good GER effects (i.e., conductive microcluster-type ID-ER fluids) as shown in Fig.6.The proportion of conductive microclusters embedded in the particles is low,no obvious effect on current density is able to be observed.

    In order to confirm the formation of conductor microclusters in TiO2particles by mixed small amount of conductors,the near neighbor structures of Ni and Au contained in TiO2has been analyzed by XAFS technique, respectively.The results show that Ni atoms are mostly oxidized with Ni–O coordination and only a small handful of Ni atoms form Ni–Ni coordination being dispersed in TiO2particles with the form of metal microclusters.Au atoms are not easy to be oxidized,and all of them form Au microclusters dispersed in TiO2particles keeping the particles being gray color.Some other easily oxidized metals,such as Al,Ti,Zn,and Cu,were mixed into TiO2to perform HEBM.All these metals were oxidized in HEBM process,or gradually oxidize after being placed in the air and then the color of TiO2particles became lighter.For example, XAFS analysis for TiO2particles mixed with 2-wt% Cu after HEBM treatment shows that there is only Cu–O coordination without Cu–Cu nearest neighbor existed.This indicates that the metal Cu atoms in TiO2particles after HEBM all are oxidized and there are no Cu metal microclusters.Detailed results of the above XAFS measurements and analyses will also be published separately.

    It should be pointed out that TiO2particles embedded into metal microclusters by HEBM also contain oxygen vacancies,and these oxygen vacancies will contribute to the ID-ER effect as well.However, the HEBM period required to inlaid metal microclusters to the particles is shorter.

    Fig.6.The relationship between yield stress and E of ID-ER fluid consisted of TiO2 particles containing graphite, Au, Ni, and stainless steel microclusters, respectively.The volume fractions of particles in ER fluids are about 50%.

    3.3.Dynamic response of shear stress to electric field E(t)for ID-ER fluid

    The response time of the shear stress to applied electric fieldE(t)for ER fluids is an important indicator of their performance.We have studied the response behaviors of ID-ER fluid in squareE(t)field and in sinusoidalE(t)field, respectively,as shown in Figs.7 and 8.The ID-ER sample was composed of TiO2particles with oxygen vacancies contained and machine oil(φ ≈48%).

    Fig.7.Response behavior of the shear stress to square wave E(t).(a) Response of shear stress to a square wave E(t) of Emax =4 kV/mm and frequency f =1 Hz(blue dash line),of which the pulse width is 400 ms and the sampling rate is 1000 s?1 at ˙γ =120 s?1.(b)Enlarged view of(a)for a single E(t)pulse.Panels(c)and(d)respectively show the normalized shear stresses at front edge and back edge of square wave E(t)for different shear rateswhere Emax =3 kV/mm (blue dotted line), f =1 Hz and sampling rate is 4000 s?1.The insets of panels (c) and (d) show the changes of tand t with ˙γ,respectively.

    The reasons for the relationship of response timetandwith ˙γcan be explained as follows.The electric polarization response is very fast and its contribution to relaxation time is negligible.Under the action of electric field, particle motion and configuration refactoring occur during shearing, resulting in the adjustment of force among all particles to cause the shear stress changing.Therefore, the correlation between shear rate and response time must be caused by the effect of shear motion on particle configuration.It was observed that the lamellar structures called shear bands can formed in ER fluids under electric field applied and shear action.[38,39]The higher the shear rate, the wider the shear bands.At the start edge ofEmax, a slightly higher shear rate is conducive to particle movement and promotes the formation of layered structures,which results in a shorter response time,as shown in Fig.7(c).However, the high shear rate has caused the existence of a wide shear band in the ER fluid asE(t)=Emax.WhenE(t)→0,it should take a long time for the structure to recover from the wide shear band to the uniform initial state and thusshould be slightly longer at higher shear rates,as illustrated in Fig.7(d).In the case of high shear rate, the configuration changes already are fast, and then the response time is not sensitive to the shear rate as showing in the insets of Figs.7(c)and 7(d).

    The response behaviors of ID-ER fluid in sinusoidalE(t)field fromf=1 Hz to 100 Hz were measured at ˙γ=50 s?1and selected plots forf=1 Hz, 5 Hz, and 50 Hz are shown in Figs.8(a)–8(c), respectively.It can be seen that with the increase of frequency, the distortion and delay of shear stress become serious,resulting in the amplitudes of shear stress decreasing with the increase of frequency.The maximum valuesτmaxand minimumτminof shear stress in sinusoidal fields from 1 Hz to 100 Hz are shown in Fig.8(d).For a comparison,the inset of Fig.8(d)shows the variation ofτmaxandτminwith the frequency for PM-ER fluid.[40]

    The above results indicate that the delay of shear stress in squareE(t)field leads to the shear stress in sine waveE(t)field decreasing with the increase of frequencyf.The amplitude of shear stress reduces by about 50% asf=50 Hz.However, compared with that in PM-ER liquid,[40,41]there is a significant improvement.

    In addition to the shear rate discussed above,the response of the shear stress of the ER to the electric field is also related to many other factors, such as volume fraction, electric field strength,oil viscosity,particle density and size,and the design of the ER device.

    Fig.8.Response of shear stress(red)with sinusoidal waves E(t)(blue)of Emax =3 kV/mm at different frequencies, shear rate ˙γ =50 s?1, sampling rate 4000 s?1.Panels(a)–(c)show the relationship between shear stress and sinusoidal E(t)for f =1 Hz,10 Hz,50 Hz,respectively.(d)The variation of τmax and τmin with f,inset shows the behavior of PM-ER fluids.

    3.4.Wear resistance

    The wear resistance of ER fluid is also an important indicator for their applications.Although PM-ER fluid can achieve high shear stress, the serious disadvantage is that the coated or modified layers on the particle surfaces are easy to wear and fall off, making its shear stress constantly reducing with service time.So that PM-ER fluid is also not suitable to be put into application.[20,21]In order to test the service life of ID-ER fluids, three different ER samples composed of TiO2particles treated by HEBM were used for wear resistance test,respectively.The samples were placed in sealed rotary cylinder, as illustrated in inset of Fig.9, to carry out friction test with a shear rate 300 s?1without electric field applied and were rubbed for 240 h,340 h,and 140 h respectively.During the period,ER samples were taken out several times for measuring the yield stress under applied electric field,and the results are shown by three cases(A),(B),and(C)in Fig.9 with red symbols.Where, case (A) is the yield stress for a conductive microcluster-type ID-ER fluid composed of TiO2particles containing 3-wt% stainless steel, cases (B) and (C) are for vacancy-type ID-ER fluids.The results show that the yield stresses of these ID-ER fluids are almost unchanged during the test and they have excellent wear resistance.By contrast, the yield stress of PM-ER fluid decreases rapidly with wear time as illustrated by cases (D) and (E) of Fig.9, represented the PM(S)-ER and PM(L)-ER fluids respectively.Here, PM(S)-ER was consisted of the solid particles with polar molecules coated while PM(L)-ER fluid was prepared by adding polar molecules to the liquid.

    Fig.9.Wear resistance of ER fluid samples: cases(A),(B),and(C)are the measured values of yield stress for three ID-ER samples during the period of wearing 240 h,340 h,and 140 h respectively(red color).The blue symbols are the wearing test results of PM-ER sample,where case(D)is PM(S)-ER fluid and case(E)is for PM(L)-ER fluid.

    Although the wear life of PM(L)-ER fluid is about five times longer than that of PM(S)-ER fluid,[21]it is still not up to the requirements of practical use.The failure of PM-ER fluid is caused either by the shedding of the polar molecules on the particles due to friction, or due to the loss of their polarities by attracting such as space charges during operation.

    It can be seen that ID-ER fluid has advantages of both high shear stress and wear resistance, which breaks through the difficulties that the ER fluid cannot be applied in practice for a long time.Therefore,the ID-ER fluid is expected to be a practical material for ER technology.

    3.5.Physical mechanism and related issues of ID-ER effect

    HEBM can significantly change the structure and microscopic morphology of TiO2nanoparticles, resulting in a large number of oxygen vacancies or embedded conductive microspheres, as shown in Fig.10(a).As for conductive microsphere-type ID-ER fluid, we have explained the mechanism in earlier publication,[23]and only that of vacancy-type ID-ER fluid will be discussed here.

    After passing the HEBM, the oxygen vacancies in TiO2particles can capture one or two electrons, or not, which can all be regarded as the local charged centers.Such charge centers,localized in the regions of vacancies,are neither like the tightly bound electrons in atoms, nor like the space charges spreading over a wide area in the medium.Oxygen vacancies can be polarized into induced dipole under an electric field,and their interactions are the source of vacancy type ID-ER effect.

    As an electric fieldEis applied,TiO2particles will be polarized and attract each other to get closer.Due to its high dielectric constant,the local electric fieldElocin the gap between two particles can be as high as 102–103times ofE.[15–19,39,42]Elocis nearly proportional toE,and its value is related to particle dielectric constant,particle size,particle spacing,and the position at the gap.Eloccan reach about 108V/m to 109V/m,whenEis 1×106V/m to 5×106V/m.

    Fig.10.(a)The internal structure and morphology of nanoparticles after HEBM.The left,the particles contain oxygen vacancies(gray);the right,the conductor microcluster (brown) embedded in the particle.(b) In the case of the particles containing vacancies.The right, an enlarged image shows that Eloc at the interval causes the induced dipole interaction from the surfaces of particles A and B.

    We can approximately calculate the induced dipole moment and the interaction force at the particle gap.The induced dipole momentμunder the action ofEloccan be expressed asμ=αEloc, whereαis the polarizability.Whileαcan be calculated by 4πεsε0r3,[16,39,43]whereεsis the dielectric constant of the solid medium,ε0is the dielectric constant of the vacuum, andris the radius of the oxygen vacancy.Asr=0.1 nm,[44]εs=100,Eloc=1×109V/m,thenμ ≈3.34 Debye with a dipole size of 0.2 nm.

    At the gaps between particles, the induced dipole on the surface of particle A can interact with that on the surface of particle B along the direction of the electric field.The force can be approximately expressed asfD?D=3μ2/2πεfε0d?D,[45]dD?Dis the center distance of the two induced dipoles,andεfis the dielectric constant of the oil.AsdD?D=0.4 nm andεf=2,thenfD?D≈1.32×10?10N with a relation offD?D∝μ2∝E2.The above estimation is only a rough approximation because the uneven distribution[15,19,42]ofElocis unconsidered.

    The interaction force between polarized particles in conventional ER fluid can be estimated roughly according to common method.[15,16]For TiO2particles with a radius of 10 nm and an interval ofδ= 0.2 nm inE= 5×106V/m field, the interaction force between two polarized particles is 4.9×10?12N.Obviously,fD?Dis more than an order of magnitude larger than the force induced by polarized particles.This can explain why the yield stress of ID-ER fluid is much higher than that of conventional ER fluid.

    The above analysis shows thatfD?Dis proportional toE2, that is the yield stressτyin ID-ER fluid should be proportional toE2, as shown in Fig.4(a).This behavior can be also observed in other ID-ER samples as shown in Figs.5 and 6.This relation ofτy∝E2of ID-ER fluid is different from that ofτy∝Ein PM-ER fluid, in which the polar molecules with intrinsic dipole moment adsorbed on the particle surface are oriented along theElocdirection and the interactions occur,and thus the forcef∝E.[9,11,12]On the other hand,the polar molecules adsorbed on the active particle surface in PM-ER fluid only exist on the particle surface and possess low adsorption energy,thus they are easily lost by friction.As shown in Fig.9,the yield stress of PM-ER fluid is reduced by about half after tens of hours of wear.Meanwhile,ID-ER fluid has good temperature stability,with the change of less than 10%in the temperature range of 25?C–170?C,[23]while the shear stress of PM-ER fluid will change by one times or even more in the 0?C–100?C range.In addition, the response time and frequency response of ID-ER fluid are improved compared with that of PM-ER fluid, maybe due to the response of induced polarization is faster than the reorientation of polar molecules.

    It is difficult to make all HEBM-prepared particles with identical function.For ER fluids with low volume fraction,the particles form chain structures when an electric field is applied.Each chain can contain about 5×105particles with a radius of 10 nm when the electrode plates are spaced 1 mm apart.The failure of individual particle can lead to entire chain ineffective.There may be not a few such failed chains between plates,resulting in too low shear stress for ID-ER fluid with low volume fraction.When the volume fraction reaches about 37% and the particles form network structure in ER suspension,[4]then the influence from these failed particles can be reduced.Therefore, only the ID-ER fluid with higher volume fraction possesses high shear stress,then the zero-field viscosity of the ID-ER fluid prepared by present process is bound to be higher.This shortcoming needs to be further improved through research.For example,the zero-field viscosity of an ID-ER fluid or the stress of it with lower volume fractions can be improved by changing the particle preparation method or changing the dispersion composition.The previous studies in PM-ER fluids have obtained successful results.[46,47]

    The vacancy-type ID-ER effect described in the present paper is due to the induced polarization of oxygen vacancies produced by HEBM technique.Can oxygen vacancies,or positive ion vacancies, produced by other methods produce similar effect? In principle,there may be similar results,but further confirmation is needed.Zhaoet al.[5]reported that doping rare earth or transition metal in oxide powder can enhance the ER effect,which may be caused by the polarization of oxygen vacancies or other defects introduced by doping.

    4.Conclusion

    The effect of traditional dielectric ER fluids originates from particle polarization interaction, which cannot meet the requirements of applications because of the low yield stress.Later on,the giant ER fluids or saying PM-ER fluids with high shear stress were developed,which were considered promising for application.However,because the particles were prepared by surface coating or modification, the polar molecules adsorbed on the surface of particles are vulnerable to wear and loss of efficacy.So that such PM-ER fluids are also unable to be widely used in practice.

    In our study,a new type of giant ER fluid called induced dipole-dominated ER(ID-ER)fluid has been developed.The preparation method,physical principle,and properties are introduced in this paper.Dielectric particles, such as TiO2, are processed by HEBM to generate a large number of oxygen vacancies or make conductor microclusters embedded in the particles.These oxygen vacancies or conductor microspheres can form induced dipoles in the electric field.Due to the high local electric fieldElocbetween particles, the dipole momentμ=αElocof the induced dipole can be large.The interactions of these induced dipoles between adjacent particles are very strong,resulting in the shear stress of the ID-ER fluid up to hundreds of kPa.The above mechanism results in the relationship between the shear stress and the electric field being asτy∝E2in ID-ER fluid.The experimental results show that in addition to high shear stress,ID-ER fluid also provides the advantages of low current density,fast response time,long wear resistance,good temperature stability,and anti-sedimentation.Moreover,the preparation method is simple and easy to repeat.

    In conclusion, the action principle and preparation method of ID-ER fluid are obviously different from either the dielectric-type ER fluid or the PM-ER fluid.This ID-ER fluid with excellent comprehensive performance should be a new generation of ER fluid suitable for practical applications.

    Acknowledgments

    Thanks to Wang Qiang, Xu Yanping, Liu Rui, Yu Richeng,Guo Qinlin,Wu Zhonghua,Zhong Jiajun,Duan Ran for their help in experiments,and to Zhou Luwei for his helpful discussion.

    Project supported by the National Key Research and Development Program of China (Grant No.2017YFA0403000)and the National Natural Science Foundation of China(Grant No.11874430).

    久久婷婷青草| 午夜福利,免费看| 性高湖久久久久久久久免费观看| 蜜臀久久99精品久久宅男| 亚洲av男天堂| 丰满乱子伦码专区| 亚洲欧美中文字幕日韩二区| 欧美精品高潮呻吟av久久| 免费看光身美女| 亚洲av中文av极速乱| 日本wwww免费看| 国产高清国产精品国产三级| a级毛片黄视频| videos熟女内射| www.av在线官网国产| 1024视频免费在线观看| 精品第一国产精品| 国产一级毛片在线| 大香蕉久久网| 久久精品久久久久久噜噜老黄| 人人澡人人妻人| 国产片特级美女逼逼视频| 乱码一卡2卡4卡精品| 亚洲成色77777| 26uuu在线亚洲综合色| 天美传媒精品一区二区| 一级毛片 在线播放| 国产精品无大码| 一区二区av电影网| 免费观看无遮挡的男女| 曰老女人黄片| 男人添女人高潮全过程视频| 亚洲av欧美aⅴ国产| 大码成人一级视频| av女优亚洲男人天堂| 又黄又爽又刺激的免费视频.| 国产欧美亚洲国产| 精品国产乱码久久久久久小说| 精品少妇久久久久久888优播| 国产精品99久久99久久久不卡 | 99视频精品全部免费 在线| 国产毛片在线视频| 日本黄大片高清| 免费观看a级毛片全部| 日韩精品有码人妻一区| 91久久精品国产一区二区三区| 9191精品国产免费久久| 美国免费a级毛片| 欧美3d第一页| 精品酒店卫生间| freevideosex欧美| 久久人人爽人人爽人人片va| 9191精品国产免费久久| 免费观看性生交大片5| www.熟女人妻精品国产 | 国产精品久久久久久久久免| 亚洲av综合色区一区| 国产免费又黄又爽又色| 老司机影院毛片| 如日韩欧美国产精品一区二区三区| 国产男人的电影天堂91| 婷婷色综合www| 纯流量卡能插随身wifi吗| 国产午夜精品一二区理论片| 人成视频在线观看免费观看| 2018国产大陆天天弄谢| 亚洲国产色片| 男女无遮挡免费网站观看| 国产男女超爽视频在线观看| 日本色播在线视频| 国产片内射在线| 我要看黄色一级片免费的| 欧美丝袜亚洲另类| 18禁国产床啪视频网站| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 久久国产亚洲av麻豆专区| 一边摸一边做爽爽视频免费| 亚洲精华国产精华液的使用体验| 一级片免费观看大全| 久久久久久伊人网av| 91在线精品国自产拍蜜月| 亚洲人与动物交配视频| 永久免费av网站大全| 桃花免费在线播放| 久久久久精品性色| 蜜桃在线观看..| 国产福利在线免费观看视频| 久久韩国三级中文字幕| 青春草国产在线视频| 国产xxxxx性猛交| 99视频精品全部免费 在线| 成人国产麻豆网| 校园人妻丝袜中文字幕| 亚洲av综合色区一区| 国产午夜精品一二区理论片| 日韩av在线免费看完整版不卡| 自拍欧美九色日韩亚洲蝌蚪91| 久久午夜福利片| 美女视频免费永久观看网站| 美女中出高潮动态图| 人体艺术视频欧美日本| 七月丁香在线播放| 18禁在线无遮挡免费观看视频| 国产成人欧美| 伦理电影大哥的女人| 婷婷色综合大香蕉| 91精品伊人久久大香线蕉| 热re99久久精品国产66热6| 日韩欧美一区视频在线观看| 午夜福利在线观看免费完整高清在| 少妇人妻久久综合中文| 午夜老司机福利剧场| 日韩制服丝袜自拍偷拍| 99精国产麻豆久久婷婷| www.色视频.com| www日本在线高清视频| 久久久国产欧美日韩av| 国产免费福利视频在线观看| 国产男女内射视频| 久久久久网色| 亚洲国产毛片av蜜桃av| 日本免费在线观看一区| xxx大片免费视频| 美国免费a级毛片| 久久精品国产鲁丝片午夜精品| 成年人免费黄色播放视频| 大香蕉久久成人网| 一级a做视频免费观看| 国产乱来视频区| 久久精品国产亚洲av涩爱| 视频在线观看一区二区三区| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 欧美人与善性xxx| 欧美日韩成人在线一区二区| 看免费av毛片| 亚洲国产成人一精品久久久| 久久国产精品大桥未久av| 欧美bdsm另类| 嫩草影院入口| 久久久久久人妻| www.av在线官网国产| 制服诱惑二区| 免费黄频网站在线观看国产| 精品少妇久久久久久888优播| 麻豆乱淫一区二区| 精品视频人人做人人爽| 成人毛片a级毛片在线播放| 亚洲精品色激情综合| 最近手机中文字幕大全| 尾随美女入室| 精品国产乱码久久久久久小说| 免费大片18禁| 热re99久久精品国产66热6| 赤兔流量卡办理| 精品福利永久在线观看| 国国产精品蜜臀av免费| 999精品在线视频| 91精品三级在线观看| 午夜av观看不卡| 亚洲成人手机| 久久精品熟女亚洲av麻豆精品| 满18在线观看网站| 国产亚洲最大av| 99热6这里只有精品| 一区二区日韩欧美中文字幕 | 最黄视频免费看| 免费大片18禁| 欧美97在线视频| 黄片无遮挡物在线观看| av卡一久久| 成年女人在线观看亚洲视频| 国产 一区精品| 丰满少妇做爰视频| 伦精品一区二区三区| 亚洲国产精品专区欧美| 另类精品久久| 国产女主播在线喷水免费视频网站| 黄片播放在线免费| 国产色婷婷99| 日韩不卡一区二区三区视频在线| 欧美日韩av久久| 免费看光身美女| 成人国产av品久久久| 亚洲av男天堂| 亚洲精品456在线播放app| 全区人妻精品视频| 看免费av毛片| 欧美变态另类bdsm刘玥| 26uuu在线亚洲综合色| www日本在线高清视频| 午夜福利,免费看| 亚洲av男天堂| 女的被弄到高潮叫床怎么办| 日韩人妻精品一区2区三区| 一区二区三区四区激情视频| 日韩视频在线欧美| 国产精品秋霞免费鲁丝片| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久这里只有精品19| 老司机亚洲免费影院| 午夜福利,免费看| 亚洲av.av天堂| 国产黄色视频一区二区在线观看| 99精国产麻豆久久婷婷| 啦啦啦中文免费视频观看日本| 亚洲国产毛片av蜜桃av| 日韩中字成人| av女优亚洲男人天堂| 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看| 国产日韩欧美视频二区| 永久网站在线| 青青草视频在线视频观看| 精品人妻一区二区三区麻豆| 交换朋友夫妻互换小说| 人妻人人澡人人爽人人| kizo精华| 精品一区二区三区四区五区乱码 | 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 99热6这里只有精品| 国产探花极品一区二区| 国产一区亚洲一区在线观看| 妹子高潮喷水视频| 在线天堂最新版资源| 高清av免费在线| 如何舔出高潮| 边亲边吃奶的免费视频| av电影中文网址| 成年动漫av网址| 最新中文字幕久久久久| 亚洲欧美一区二区三区国产| 中文字幕亚洲精品专区| 成人毛片60女人毛片免费| 视频中文字幕在线观看| 国产一区有黄有色的免费视频| 免费女性裸体啪啪无遮挡网站| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃 | 美女内射精品一级片tv| 熟女av电影| 免费观看a级毛片全部| av女优亚洲男人天堂| 免费在线观看完整版高清| 最黄视频免费看| 一二三四在线观看免费中文在 | 久久这里只有精品19| 丰满少妇做爰视频| 人人妻人人爽人人添夜夜欢视频| 亚洲成国产人片在线观看| 美女视频免费永久观看网站| 寂寞人妻少妇视频99o| 久久97久久精品| 午夜影院在线不卡| 国产精品一国产av| 人妻系列 视频| 国产在线免费精品| 美女xxoo啪啪120秒动态图| 精品99又大又爽又粗少妇毛片| 国产永久视频网站| 午夜福利视频精品| 亚洲欧美成人精品一区二区| 18禁裸乳无遮挡动漫免费视频| 妹子高潮喷水视频| 亚洲一码二码三码区别大吗| 美国免费a级毛片| 亚洲人成网站在线观看播放| 亚洲欧美成人综合另类久久久| 国内精品宾馆在线| 久久鲁丝午夜福利片| 国产激情久久老熟女| 久热久热在线精品观看| 最新的欧美精品一区二区| 高清不卡的av网站| 亚洲久久久国产精品| 桃花免费在线播放| 亚洲国产精品国产精品| 只有这里有精品99| 在现免费观看毛片| 伦理电影大哥的女人| 亚洲av电影在线观看一区二区三区| 国产爽快片一区二区三区| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 久久久国产欧美日韩av| 99热网站在线观看| 在现免费观看毛片| 天堂中文最新版在线下载| 伦理电影免费视频| 午夜视频国产福利| 亚洲av.av天堂| 在线亚洲精品国产二区图片欧美| 午夜福利网站1000一区二区三区| 高清欧美精品videossex| 久久av网站| 一边亲一边摸免费视频| 自线自在国产av| 超色免费av| 久久精品国产综合久久久 | av天堂久久9| 日韩视频在线欧美| 日韩精品有码人妻一区| 女人被躁到高潮嗷嗷叫费观| 1024视频免费在线观看| 欧美日韩av久久| 国产一区二区激情短视频 | 久久ye,这里只有精品| 美女视频免费永久观看网站| 国产精品一国产av| 一二三四在线观看免费中文在 | 两个人看的免费小视频| 巨乳人妻的诱惑在线观看| 大话2 男鬼变身卡| 国产免费一区二区三区四区乱码| 亚洲久久久国产精品| 丰满乱子伦码专区| 天天影视国产精品| 日日撸夜夜添| 欧美亚洲 丝袜 人妻 在线| 啦啦啦视频在线资源免费观看| 午夜福利乱码中文字幕| 国产激情久久老熟女| 老熟女久久久| 精品国产一区二区久久| 人妻系列 视频| 日韩人妻精品一区2区三区| 嫩草影院入口| 亚洲经典国产精华液单| 色视频在线一区二区三区| 久久毛片免费看一区二区三区| 国产成人91sexporn| 三上悠亚av全集在线观看| 成人毛片a级毛片在线播放| 一本久久精品| 女性生殖器流出的白浆| √禁漫天堂资源中文www| 丝瓜视频免费看黄片| 制服丝袜香蕉在线| 国产老妇伦熟女老妇高清| 亚洲,欧美,日韩| 色吧在线观看| 婷婷色综合www| 国产爽快片一区二区三区| 丝袜美足系列| 欧美成人精品欧美一级黄| 成年人免费黄色播放视频| 国产极品天堂在线| 又粗又硬又长又爽又黄的视频| 晚上一个人看的免费电影| 久久久久久久久久人人人人人人| 久久久久精品久久久久真实原创| 欧美xxxx性猛交bbbb| 亚洲欧美成人精品一区二区| 成人黄色视频免费在线看| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 哪个播放器可以免费观看大片| 黄色毛片三级朝国网站| 内地一区二区视频在线| 国产免费福利视频在线观看| 国产精品欧美亚洲77777| 黄网站色视频无遮挡免费观看| 九草在线视频观看| 成人国语在线视频| 久久久精品区二区三区| 久久久久久久久久人人人人人人| www.熟女人妻精品国产 | 亚洲精华国产精华液的使用体验| 26uuu在线亚洲综合色| 午夜激情av网站| 亚洲精品美女久久久久99蜜臀 | 亚洲av.av天堂| 亚洲欧美清纯卡通| 男人添女人高潮全过程视频| 最新的欧美精品一区二区| 成人影院久久| 精品亚洲成a人片在线观看| 日本黄大片高清| √禁漫天堂资源中文www| 国产成人欧美| 亚洲一码二码三码区别大吗| 免费大片18禁| 久久精品国产a三级三级三级| 九九在线视频观看精品| 国产成人精品婷婷| 欧美成人午夜精品| 欧美激情国产日韩精品一区| 草草在线视频免费看| 中文字幕精品免费在线观看视频 | 亚洲内射少妇av| 狂野欧美激情性bbbbbb| 99国产精品免费福利视频| av片东京热男人的天堂| 九色亚洲精品在线播放| 欧美日韩一区二区视频在线观看视频在线| 欧美另类一区| 如何舔出高潮| 国产av码专区亚洲av| 免费观看性生交大片5| 亚洲精品一区蜜桃| 欧美日韩成人在线一区二区| 夜夜骑夜夜射夜夜干| 巨乳人妻的诱惑在线观看| 免费高清在线观看视频在线观看| av卡一久久| 国产不卡av网站在线观看| a级毛片在线看网站| 9191精品国产免费久久| 青春草视频在线免费观看| 免费人妻精品一区二区三区视频| 成人漫画全彩无遮挡| 内地一区二区视频在线| 久久久亚洲精品成人影院| 综合色丁香网| 97在线人人人人妻| av网站免费在线观看视频| 免费观看性生交大片5| 精品人妻在线不人妻| 视频区图区小说| 亚洲国产色片| 国产在线视频一区二区| 草草在线视频免费看| 亚洲成人一二三区av| 99国产综合亚洲精品| 亚洲欧洲日产国产| 久久人人97超碰香蕉20202| 亚洲精品乱码久久久久久按摩| 亚洲高清免费不卡视频| 亚洲成色77777| 人人澡人人妻人| 大香蕉久久网| 午夜福利视频精品| 亚洲经典国产精华液单| 精品少妇内射三级| 亚洲欧美精品自产自拍| 美女国产高潮福利片在线看| 日本黄大片高清| 欧美xxxx性猛交bbbb| 22中文网久久字幕| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 狂野欧美激情性bbbbbb| 哪个播放器可以免费观看大片| 亚洲成人一二三区av| 搡女人真爽免费视频火全软件| 中文字幕制服av| 成人黄色视频免费在线看| av线在线观看网站| 日韩制服骚丝袜av| 夜夜爽夜夜爽视频| 大片电影免费在线观看免费| 午夜福利视频在线观看免费| 欧美国产精品一级二级三级| 成人午夜精彩视频在线观看| 老司机亚洲免费影院| 肉色欧美久久久久久久蜜桃| 亚洲精品日本国产第一区| 精品国产一区二区三区久久久樱花| 18禁在线无遮挡免费观看视频| 波野结衣二区三区在线| 老司机影院成人| 亚洲国产精品一区二区三区在线| 秋霞在线观看毛片| 街头女战士在线观看网站| 少妇熟女欧美另类| 日韩不卡一区二区三区视频在线| 高清在线视频一区二区三区| 日韩av免费高清视频| 热re99久久精品国产66热6| 亚洲av欧美aⅴ国产| 妹子高潮喷水视频| 亚洲熟女精品中文字幕| av在线app专区| 成人国语在线视频| 国产又色又爽无遮挡免| 人人澡人人妻人| 国产精品国产av在线观看| 亚洲,一卡二卡三卡| 国产一区二区激情短视频 | 中国三级夫妇交换| 日本黄色日本黄色录像| 国产一区二区三区av在线| 熟女人妻精品中文字幕| 制服丝袜香蕉在线| 美国免费a级毛片| 你懂的网址亚洲精品在线观看| 熟女人妻精品中文字幕| 亚洲精品久久久久久婷婷小说| 久久久久久人妻| 草草在线视频免费看| 老熟女久久久| 国产精品三级大全| 国产午夜精品一二区理论片| 男女国产视频网站| 在线精品无人区一区二区三| 看免费av毛片| 午夜福利在线观看免费完整高清在| 在线 av 中文字幕| 国产精品偷伦视频观看了| 黄色 视频免费看| 成人手机av| 亚洲av综合色区一区| 又大又黄又爽视频免费| 亚洲精品一区蜜桃| 免费av中文字幕在线| 亚洲欧美成人精品一区二区| 久久精品久久精品一区二区三区| av视频免费观看在线观看| 中文字幕最新亚洲高清| 免费在线观看完整版高清| 国产亚洲午夜精品一区二区久久| 九草在线视频观看| 国产精品久久久久久av不卡| 久久久久久久久久久久大奶| 建设人人有责人人尽责人人享有的| 国产麻豆69| 亚洲国产毛片av蜜桃av| 大码成人一级视频| 国产探花极品一区二区| 美女福利国产在线| 男女啪啪激烈高潮av片| 激情视频va一区二区三区| 桃花免费在线播放| 欧美97在线视频| 日韩成人av中文字幕在线观看| 成人亚洲精品一区在线观看| 精品国产国语对白av| 日韩人妻精品一区2区三区| 欧美成人午夜精品| 日日撸夜夜添| 成年女人在线观看亚洲视频| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 欧美 亚洲 国产 日韩一| 欧美精品一区二区免费开放| 男人舔女人的私密视频| 国产精品三级大全| 99re6热这里在线精品视频| 久久国产亚洲av麻豆专区| 成人亚洲欧美一区二区av| 国产精品久久久久久精品电影小说| 免费高清在线观看视频在线观看| 少妇人妻精品综合一区二区| 亚洲av成人精品一二三区| 国产精品久久久av美女十八| 校园人妻丝袜中文字幕| 大陆偷拍与自拍| 侵犯人妻中文字幕一二三四区| 成人手机av| 狠狠精品人妻久久久久久综合| 午夜久久久在线观看| 伦精品一区二区三区| 在现免费观看毛片| 日产精品乱码卡一卡2卡三| 日本91视频免费播放| 夫妻午夜视频| 久久久久久久久久人人人人人人| 美女国产高潮福利片在线看| 日韩精品有码人妻一区| 巨乳人妻的诱惑在线观看| 国产av国产精品国产| 欧美日韩综合久久久久久| 久久国内精品自在自线图片| 一级片免费观看大全| 建设人人有责人人尽责人人享有的| 亚洲一级一片aⅴ在线观看| 午夜视频国产福利| 精品久久久精品久久久| 天天躁夜夜躁狠狠久久av| 亚洲国产精品国产精品| 国产精品免费大片| 欧美最新免费一区二区三区| 亚洲精品日本国产第一区| 国产国语露脸激情在线看| 国产xxxxx性猛交| 亚洲精品av麻豆狂野| 午夜福利在线观看免费完整高清在| 丰满饥渴人妻一区二区三| 国产精品秋霞免费鲁丝片| 国产日韩欧美在线精品| 国产免费视频播放在线视频| 波野结衣二区三区在线| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 热99久久久久精品小说推荐| 国精品久久久久久国模美| 亚洲激情五月婷婷啪啪| 狠狠精品人妻久久久久久综合| 国产亚洲最大av| av有码第一页| 欧美成人午夜精品| 亚洲国产av新网站| 午夜福利,免费看| 午夜老司机福利剧场| 日本色播在线视频| 又黄又爽又刺激的免费视频.| 中文字幕av电影在线播放| 永久网站在线| 久久人人97超碰香蕉20202| 2018国产大陆天天弄谢| 老司机影院毛片| 97人妻天天添夜夜摸| www.色视频.com| 另类亚洲欧美激情| 大香蕉97超碰在线| 蜜臀久久99精品久久宅男| 日韩视频在线欧美| 热99久久久久精品小说推荐| 国产1区2区3区精品| 成年美女黄网站色视频大全免费| 熟妇人妻不卡中文字幕| 大码成人一级视频|