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

    Insights into the regulation mechanism of ring-shaped magnetoelectric energy harvesters via mechanical and magnetic conditions?

    2021-10-28 07:15:52YangShi師陽(yáng)NiLi李妮andYongYang楊勇
    Chinese Physics B 2021年10期
    關(guān)鍵詞:楊勇

    Yang Shi(師陽(yáng)) Ni Li(李妮) and Yong Yang(楊勇)

    1School of Mechano-Electronic Engineering,Xidian University,Xi’an 710071,China

    2Research Center for Applied Mechanics,Xidian University,Xi’an 710071,China

    Keywords: magnetoelectric effect,boundary conditions,magneto-elastic coupling,external stimuli

    1. Introduction

    The emerging industrial revolution involves the internet of things(IoT),which is the key technology for enabling smart cities.[1]However, although this technology would promote innovation in society, as we would create a network between everything and realize remote distance control of electronic devices,some challenges exist in bringing the IoT to fruition.For example,this concept assumes the very wide implementation of distributed electronic devices to capture data,[2,3]each of which will require energy to run. Batteries in various forms are currently the primary power source for operating these devices. Despite many advantages,batteries are not the ultimate solution because they require periodic replacement or recharging. One of the promising solutions to these issues is harvesting ambient forms of energy,much of which would otherwise be wasted, to power the sensors. Among the various sources of energy,magnetic field is an attractive alternative source because they are available in most of the locations where the wireless sensor networks will be used.

    Recently, laminate magnetoelectric (ME) composites,formed by elastically combining piezoelectric and magnetostrictive materials, have been demonstrated as a potential choice for harvesting magnetic signals due to strong strain-mediated ME coupling.[4–6]In the ME effect, an electric field or voltage through the piezoelectric effect is resulted from the mechanical strain induced by an external magnetic field. The strength of the ME coupling in laminated composites depends mainly on the intrinsic property of components,[7]interfacial mechanical coupling,[8–10]and the boundary conditions.[11]Up to date,various investigations have been performed to achieve higher ME coupling coefficient by selecting magnetostrictive and piezoelectric components with excellent properties.[12,13]On the other hand,a great deal of researches have been devoted to further enhance ME coefficient via varying the configuration of ME composites.[14,15]Among various configurations, the diskring, ring-shaped and cylindrical composites show considerable advantages due to improved interfacial mechanical coupling. Unlike the conventional laminated plate or beam, the ring-shaped or cylindrical structures can transfer normal stress instead of shear stress via interface and thus improve the interfacial mechanical coupling.[16]In addition,some advanced fabrication methods, such as electro or electroless depositions,have been successfully used in preparing the ring-shaped composites.[17]This can eliminate the use of epoxy bonding layer or electrodes, and is thereby believed to be advantageous. Wanget al.[8,18]theoretically revealed that ME effect can be substantially enhanced by increasing the curvature of a piezoelectric/piezomagnetic cylinder, and an imperfect interface can annul small unwanted resonance for better device operation. Gao and Zhang[19]modeled the nonlinear ME coupling of a giant magnetostrictive/piezoelectric composite cylinder based on a standard square constitutive model. Most recently, Kumar and Arockiarajan[20]presented an analytical ME model in a cylindrical coordinate system incorporating the temperature effects on the individual composite constituents.These theoretical models have significantly contributed to the theory of ME effect in ring-shaped structures.

    Nearly, all prior works on ring-shaped ME composites focused on studying ME coupling with consideration of an open-circuit condition, thus they can only apply for sensing applications rather than harvesters. For the energy harvester using ring-shaped composites,however,the power conversion performance has not got enough attention. There are still considerable difficulties and challenges in the practical application of the ring-shaped magnetic energy harvester. Initially,it has been indicated that available ME effect can be observed only when the intrinsic frequency of structure closes to the frequency of ambient magnetic field.[21–23]However,the magnetic fields from the wireless sensor networks always have random distributions of frequency components, which makes it difficult to match the frequencies of the harvester and ambient magnetic field. Therefore, exploring new regulation mechanism of ME performance becomes one of important issue to be addressed. Mechanical boundary conditions may directly influence the strain distribution in the ring-shaped composites,and improve the performance of energy harvester. The insight into the ME coupling under different mechanical boundary conditions is currently lacking in any of the above mentioned references nor available in the state-of-the-art literature.Moreover,the material parameters(i.e.,elasticity modulus and piezomagnetic coefficient) of magnetostrictive materials are very sensitive to the external stimuli like bias magnetic and stress loadings.[24–28]This enables a new theoretical mechanism for tunable or active manipulation of ME performance of the energy harvester.[29]However, the nonlinear coupled effects of such external stimuli on ME coupling of the ringshaped composites have not been well investigated since most of the theoretical models were established using linear constitutive relations of constituent materials.

    The objective of this paper is to model ring-shaped ME energy harvesters and to reveal regulation mechanism via mechanical and bias magnetic conditions. The theoretical model for the energy harvesters based on ring-shaped ME composites is established and the basic equations are introduced in Section 2. Then analytical solutions for the ME coupling in ring-shaped composites and nonlinear field-dependent material parameters are derived in Section 3. In Section 4,some results are presented for tuning the ME performance via boundary conditions and external stimuli. Finally,some conclusions are made in Section 5.

    2. Problem description

    Figure 1 shows a magnetic energy harvester based on a ring-shaped ME composite. The ME composite is presumed to be assembled from an outer piezoelectric ring concentric to an inner magnetostrictive cylinder. Taking into account the symmetry of the cylindrical structure,a cylindrical coordinate system(r,?,z)is adopted, in which the coordinate directionr,?andzwill be designated as 1,2 and 3,respectively,so that the tensor coefficients can keep the same values in a cylindrical coordinate system as those in a Cartesian system. The piezoelectric material is polarized alongzdirection and magnetostrictive material is magnetized along the same direction.Herea,bandcare radiuses as denoted in Fig.1. The ambient ac magnetic field alongzdirection could induce a mechanical strain from the magnetostrictive phase that transfers to the piezoelectric phase via interfacial coupling. Then a change in polarization of the piezoelectric phase results from piezoelectric effect, which is expected to be harvested by an external resistanceRl. Note that in practical applications, the postprocessing electric circuit is very crucial for ME harvesters.Adequate power management circuits have to be applied in electromagnetic energy harvesting to produce appropriate supply performance. In the current work we intend to study clearly how mechanical and magnetic conditions affects ME coupling. Thus,only one-load electric circuit[30,31]is considered,which,however,may not hinder us from evaluating and comparing ME performance under different stress and magnetic loadings. As the height of the cylindrical structure is assumed to be much less than its radius in the current analysis,the only nonzero components of the stress tensor areσrrandσ??. Since for the harmonic motion, all the fields have the same time dependence exp(jωt). This factor will be dropped in the analysis. For macro-scale magneto-electro-elastic media,both flexo-effect and surface effect can be neglected.[32,33]In a cylindrical coordinate system(r,?,z),when such media were operated under the conditions that the body force, free charge density and current density are absent, the mechanical governing equations for the dynamics of the axisymmetric cylindrical ring can be expressed as[34]whereuris radial displacement, andρis the effective mass density of the ring-shaped composites.

    Fig.1. Schematic diagram of magnetic energy harvester based on ringshaped ME composites.

    The linear constitutive equations for the stress tensorpσi jand the electric displacementD3in the piezoelectric layer,and stress tensorsmσijin the magnetostrictive layer can be written in terms of displacement as follows:

    where the left superscriptspandmrepresent piezoelectric and magnetostrictive layers,respectively;ci jkl,eij,di j,anda33are elastic,piezoelectric,piezomagnetic,and dielectric constants,respectively;E3andH3are the components of electric field and magnetic field,respectively.

    3. ME coupling in ring-shaped composites

    3.1. General solutions for mechanical variables

    For solving mechanical variables,we substitute the above constitutive equations into Eq. (1), and rewrite the governing equation in the form of

    whereAiandBiare the integration constants determined from interfacial and boundary conditions.J1(·) andY1(·) are the first and second kind Bessel functions. The stresses and electric displacement in the composites are then obtained with combining the above solutions and constitutive relations.

    3.2. Mechanical interfacial and boundary conditions

    Since the thickness of the interfacial layer caused from adhesives may be very thin compared to the composites, the interface can be classified as a bonding layer-free interface.Thus, perfect interfacial conditions are considered in this study,given by

    Some conventional boundary conditions are listed in Table 1, which include stress-free and clamping (zero displacement) separately on the inner and outer radii of the magnetostrictive and piezoelectric layers,respectively. Note that one of the great advantages of the ring-shaped structure is that one can apply normal stress on the inner and outer boundaries,without the concern of interface failure caused by shear stress.In the current study,we propose new mechanical conditions by applying stressP0on the inner radii of magnetostrictive layer or the outer radii of piezoelectric layer. These mechanical conditions can be realized by assembling a concentric cylinder (annulus) with proper dimensions to the ME composites at inner(outer)boundary,since radical direction stress may be produced at boundaries via strong compression. As listed in Table 2, the proposed stress boundary conditions differ from those used in the overwhelming majority of existing analytical models.[19,35]The introduction of these conditions provides a utility for tuning the stress distribution in the composites,hence tuning the overall performance of the magnetic energy harvester. This highlights the first novelty of this study while accounting for possible practical operation issues.

    Table 1. List of stress-free and clamping boundary conditions and their mechanical interpretations.

    Table 2. List of applying stress boundary conditions and their mechanical interpretations.

    Substituting the expressions of the stress tensor and the displacements into the above mechanical interfacial and boundary conditions leads to eight equation sets with respect toAiandBi. These equation sets are then solved for unknown constants and the induced displacements in piezoelectric phase under different conditions can be determined.

    3.3. Analytical solutions for electrical variables

    In this section,the induced electrical variables are solved.It can be predicted from the above process that the displacement in piezoelectric phase has linear relations with the induced electric field, the applied magnetic field, and external stress applied at the boundaries. For convenience, the displacementpuris rewritten as

    The magnetic field harvester schematically displayed in Fig. 1 is operated under closed circuit with a resistanceRl.The induced ME voltage is denoted asVME.Combining Ohm’s law and the integral form of Gauss’s law, we can reach the following relation:

    Then using the expression of the electric displacement and considering the relation between the electric field and the voltage,we obtain the expression of the induced ME voltage under closed-circuit condition,given as

    Note that our solutions can also be used to predict the ME effect of ring-shaped composites used as magnetic sensors under open-circuit condition. In this case,the only change is setting the right part of Eq.(8)(i.e.,the induced current)to be zero. That is,Rl→∞yields equation that determines the open-circuit ME voltage as follows:

    Hence the ME voltage coefficientαMEand the output power densitypis obtained as

    whereIis the output current,and the operation Re is for taking the real part of a complex number.

    3.4. Nonlinear constitutive theory for magnetostrictive materials

    The solutions given in the above section is based on the linear constitutive relations of magnetostrictive materials.That is,material parameters of magnetostrictive materials,such as elastic tensor and piezomagnetic coefficient tensor,are considered as constants under fixed pre-stress and bias magnetic field. However, most magnetostrictive materials exhibit a complex nonlinear magneto-elastic coupling,[36–38]i.e., the material properties are significantly affected by the pre-stress and bias magnetic field. Considering the magneto-elastic coupled characteristic may result in some nonlinear ME behaviors,and set up new theoretical mechanism for tuning ME performance of the energy harvester. The nonlinear constitutive relations have been explored extensively and proposed on the basis of the Gibbs free energy function as follows:[39]

    wheresis the compliance coefficient tensor,Mis the magnetization tensors,Iis the second-rank unit tensor,λsandMsare the saturated magnetostriction and magnetization,kis relaxation factor,μ0is vacuum permeability, andf(x) is the Langevin function. The termε(σ,M) represents the magnetostriction related to magneto-elastic coupling. Readers can refer to Ref.[24]for the more detailed expressions ofs,ε(σ,M)and ˉσ.

    Let us consider now the ME energy harvester response to a combined magnetic fieldH=Hdc+h(f), whereHdcis the applied bias field andh(f)is ac magnetic field. Upon expanding the strains as a Taylor series in the vicinity of the applied biasHdc,we obtain

    Note here that the above linear-like form is fully magnetoelastic coupled since all the material parameters are strictly derived from the nonlinear magneto-elastic coupling constitutive relations. In order to unify the nonlinear and linear relations in form,the relation in Eq.(14)can be rearranged as

    Fig.2. (a)Elastic compliance coefficient and(b)piezomagnetic coefficient of Terfenol-D versus magnetic field under different pre-stresses.

    4. Results and discussions

    In this section, numerical calculations are performed for the ring-shaped composites and some significant predictions are presented based on our model. The piezoelectric material is chosen as BaTiO3, and the magnetostrictive materials for linear and nonlinear cases are chosen as CoFe2O4and Terfenol-D, respectively. The material parameters used for numerical calculations are listed in Tables 3 and 4.[44,45]The volume fraction of the piezoelectric phase is defined byv=(c?b)/(c?a). The inner and outer radii of the ring are fixed ata=9 mm andc=13 mm, respectively. The ac field isHac=1 Oe, the resistance isRl=50 k?, and the volume fraction isv=0.5. For studying the resonance ME coupling,a damping coefficientχis used to describe resonance loss.The damping coefficient is determined byχ=1/Q=?f/fr,whereQis the quality factor for the observed resonance,fris resonance frequency,and ?fis the bandwidth. Wanet al.[15]experimentally obtained the effective quality factor of~48 for the first bending resonance and~51 for the first longitudinal resonance, indicatingχ=0.02 may be believable in numerical calculation. Zhanget al.[46]and Zhaoet al.[47]also usedχ=0.02 to evaluate the resonance ME effects in multiferroic laminates plates and cylinders. In the current study,ω′=(1+0.02i)ωis taken into account in the numerical calculation.

    Table 3. Material parameters of BaTiO3 and CoFe2O4.

    Table 4. Material parameters of Terfenol-D.

    4.1. Model degeneration and validation

    It is noteworthy that the theoretical model in the above sections is established for ring-shaped structures,which,however, can be degenerated into the model for disk/ring structures. In the configurations of disk/ring structures, a piezoelectric disk is inserted into a magnetostrictive ring (denoted by P/M),or a magnetostrictive disk is inserted into a piezoelectric ring(denoted by M/P).Unlike the ring-shaped structures,the inner radiusaof the disk/ring structures is set to be zero,and the displacement at the center of disk should not be infinite,i.e.,one of the undetermined coefficients in Eq.(5)must be zero. Recently, the dynamic ME response of a Terfenol-D disk/PZT ring has been well reported experimentally by Li and Chen.[48]Here, we compare the predicted ME effects of Terfenol-D disk/PZT ring with experimental data in Ref.[48],and plot the open-circuit ME coefficient as a function of the driving frequency in Fig.3. Note that for the results presented in Fig.3,the inner diameter,the outer diameter,and the thickness are 10 mm,12 mm,and 0.95 mm,respectively. It may be observed that the proposed scheme is able to capture both the maximum ME coefficient and the resonance frequency. The two sets of results are in very close agreement,indicating that the theoretical model is sufficiently accurate for characterization of dynamic ME effect.

    Fig. 3. ME coefficient of Terfenol-D disk/PZT ring with varying the driving frequency.

    4.2. Regulation mechanism via mechanical boundary conditions

    The results presented in Fig. 4 report the ME voltage of ring-shaped energy harvester versus the driving frequency of magnetic field for both conventional and stress boundary conditions listed in Tables 1 and 2. The results for conventional boundary conditions are calculated by setting the applied stressP0in Eq. (9) to zero. Since the type of boundary conditions dictates the distribution of displacement field throughout the composite,the induced ME voltages under different boundary conditions change variously. For some of the conventional boundary conditions (i.e., F-F, F-C, and CF),two obvious resonance frequencies are obtained within the range of 0–1400 kHz,as shown in Fig.4(a). In addition,setting a stress-free condition makes it easy to induce a resonance response at low frequency, while setting a clamped condition is beneficial for tuning a larger resonance frequency. Thus,the condition F-F corresponds to the lowest resonance frequency,whereas the condition C-C to the largest one. Except for the type of boundary condition, the location at which a specific condition is applied may affect the transfer of displacement from one layer to the other, and therefore matters. It can be found that setting the outer boundary stress-free may result in notable ME coupling at a lower frequency. For example,the maximum ME voltages for F-F and C-F conditions are achieved at the fundamental frequency, while the F-C condition shows notable ME voltage at the second-order frequency.The expression in Eq.(9)indicates that ME voltage is significantly related to the displacements at the outer boundaryr=cand the interfacer=b. Thus,the effect of boundary condition on ME coupling is twofold. First, the type of boundary condition changes the radial displacement at the outer boundary.For example, the conditions F-C and C-C dictate a zero displacement at the outer boundary. Second,the solutions of the undetermined coefficients in Eq.(9),determined by the values and types of boundary conditions,will indirectly affect the displacement at the interface through the solutions of two Bessel differential equations. Based on these mechanisms, ME performance of the proposed ring-shaped energy harvester may be tuned by optimizing the prescribed value of the displacement on the boundaries.

    In the current theory,when the stress boundary conditions are applied,the only operation is to replace the zero stress by a nonzero one for a boundary. In other words, relaxing the conditions S-F or F-S from applied stress to stress-free results in the condition F-F. Similarly, the corresponding stress-free condition for S-C and C-S are F-C and C-F,respectively.Comparing Figs.4(a)and 4(b),we find that applying stress can improve ME voltage,without changing the resonance frequency of the corresponding conditions. The peak ME voltages and resonance frequencies under different boundary conditions are summarized in Fig.5. As can be seen, the peak ME voltages under conventional boundary conditions are significantly enhanced by applying stress at the free boundary. For example,the peak ME voltage under the condition F-F is increased from 0.63 to 1.80 V after applying a 6-kPa stress at the outer boundary. In essence, piezoelectric effect is improved by the additional mechanical strain when at least one of boundaries is subjected to applied stress,resulting in the enhanced output voltage.

    Fig.4. ME voltage of ring-shaped energy harvester versus the driving frequency of magnetic field for(a)conventional and(b)stress boundary conditions(P0=6 kPa).

    Figure 6 shows the ME voltage versus the driving frequency with varying the applied stress under different boundary conditions. For all the stress boundary conditions,the results with a zero stress reduce to those for the corresponding conventional conditions. It can be seen that as the stress increases,ME voltage is improved continually,demonstrating a monotonically increase. However, it should be noted that in actual operation the applied stress may be properly increased,always assuming that its valueP0should not exceed the allowable stress[σ]of constituent materials. The present model shows that the effect of stress boundary conditions on ME coupling is also twofold. First, once a stress is applied at inner or outer free boundary,it influences the solutions of displacement in Eq.(5)by changing the values of undetermined coefficientsAiandBi,In other words,stress boundary condition will change the strain distributions in the composites,and thus may improve both piezoelectric and piezomagnetic effects. Second,the applied stress increases the second term in the numerator of Eq. (9), and thereby additional voltage is directly obtained via the improved piezoelectric effect. Equation(9)also suggests that ME performance can be enhanced,only with the displacement related to stress boundary condition being of the same sign as that related to the applied magnetic field. This provides a theoretical guidance for a specific operation.

    Fig. 5. Maximum ME voltages and the corresponding response frequencies under different boundary conditions(P0=6 kPa).

    Figure 7 shows the output power density versus resistance under different boundary conditions. For any boundary condition, a maximum power density is obtained at an optimal resistanceRo. The optimal resistances are different for different conventional conditions. As is expected,peak power densities for conventional conditions F-F,F-C,and C-F are enhanced by applying stress at the original stress-free boundaries. Interestingly,but not surprisingly,the optimal resistance is unchanged after setting the applied stress. The peak power densities,and the corresponding optimal resistances and resonance frequencies are summarized in Fig. 8. The C-C boundary condition,being the one that sets the displacements at both boundaries to zero,may give rise to the stress built up in the composite,and thus exhibits notable ME performance. In addition,when the composite is constrained at least on one side, the optimal resistance appears to be reduced and the resonance frequency be increased. Consequently,ME harvester under the C-C boundary condition exhibits the maximum power density with the minimum optimal resistance and highest resonance frequency.In comparison, a remarkable fact is that the C-S boundary condition reduces the resonance frequency for the C-C boundary condition by about 40%with almost unchanging the peak power density at a larger optimal resistanceRo=15 k?.In all,there is still much room for improvement in ME performance since the applied stress in Fig.8 is only 6 kPa,which is far less than the allowable stress of constituent materials.

    Fig.6. ME voltage versus the driving frequency for different boundary conditions: (a)C-S,(b)S-C,(c)F-S,and(d)S-F.

    Fig.7. Output power density versus resistance for different boundary conditions: (a)C-S,(b)S-C,(c)F-S,and(d)S-F.

    Fig.8.Maximum power densities and the corresponding optimal resistances for different boundary conditions(P0=6 kPa).

    4.3. Regulation mechanism via multi-field stimuli

    In this section, we intend to study the nonlinear multifiled coupled mechanisms of ME energy harvester and tune ME performance via external stimuli. Experiments have demonstrated that a typical giant magnetostrictive material Terfenol-D exhibits strong magnetostrictive effect and complex nonlinear coupling characteristics when it is subjected to combined stress and magnetic loadings. Therefore, we select Terfenol-D as the magnetostrictive phase of the ringshaped composites and expect to obtain some novel phenomena. The material parameters used for numerical calculations are listed in Table 4. Unless otherwise mentioned,the applied pre-stress is fixed atσ0=?10 MPa. Figure 9 displays the ME voltage versus the applied bias magnetic field for different frequencies. Note that the results presented in Fig. 8(a)include ME voltages of Terfenol-D/BaTiO3using nonlinear model and CoFe2O4/BaTiO3using linear model. It can be seen that the ME voltage predicted from linear model is independent of the applied bias field at any frequency. However,when the nonlinear multi-field coupling of Terfenol-D is taken into account in the nonlinear model,an interesting phenomenon happens. That is, the ME voltage firstly increases and then decreases with bias magnetic field, indicating that an optimal bias field may be used to tune ME performance of the ring-shaped composites. This highlights the second novelty of this study while accounting for active control of ME energy harvesters. Generally, ME composites need to be operated at an external bias magnetic field for available ME response since piezomagnetic effect vanishes at zero field. This,however,inevitably increases the excess energy consumption.Most of the conventional ME composites used a pair of electromagnets to produce bias magnetic fields,because it is convenient to adjust the magnitude of bias field. However, the electromagnets may cause an increase in the complexity of the device and a compromise in the performance. In the current study we present a simple and effective way to design a self-biased ME energy harvester. Four pieces of plate-shaped NdFeB magnets are arranged in the 0°, 90°, 180°, and 270°directions in a ring-shaped aluminum mold cavity of dimensions the same as the magnetostrictive ring so as to create a bias magnetic field. With the NdFeB magnets in place, the Terfenol-D short fibers are aligned and then the epoxy is transferred into the mold. Therefore, a magnetically biased magnetostrictive composite ring is successfully fabricated. Based on the proposed nonlinear theoretical model,the optimal bias fields for different application scenarios can be accurately predicted in advance. Then,the precise static magnetic field may be produced by choosing NdFeB magnets with proper dimensions.It can be seen from Fig.9(a)that at comparatively lower frequencies only one peak appears in the ME voltage versus bias field curve of Terfenol-D/BaTiO3. As the frequency increases, the peak position moves slightly toward higher bias field, and the ME coupling effect is enhanced considerably.According to the interfacial conditions in Eq.(6),the termuHin Eq.(7)is proportional to(equivalent)piezomagnetic coefficient. Thus, we can deduce from Eq. (7) that ME voltage should be also proportional to the(equivalent)piezomagnetic coefficient. The solutions for nonlinear material parameters in section 3.4 show that the equivalent piezomagnetic coefficient depends strongly on the applied bias field. That is why ME voltage is markedly affected by the bias field. Essentially,the peak position of ME voltage corresponds to the optimal bias filedHoptat which the maximum value of the equivalent piezomagnetic coefficient of Terfenol-D occurs. Figure 9(a)also indicates that low frequency has a slight effect on theHopt. At frequencies from 400 to 430 kHz, the curve of ME coefficient varying with bias field exhibits a dual-peak phenomenon,namely,two peaks appear with the increase of bias field. A similar dual-peak phenomenon was recently experimentally observed in laminated ME composites by Zhaoet al.[49]and our group.[50]With increasing the frequency, the first peak moves toward lower bias field while the second one toward higher bias field. For example,the first peak appears at 216 Oe forf=400 kHz,and 76 Oe forf=430 kHz,respectively; the second peak appears at 560 Oe forf=400 kHz,and 896 Oe forf=430 kHz,respectively. At the same time,both the first and second peaks are reduced by increasing frequency. Since the dual-peak phenomenon appears at larger frequency and the optimal field corresponding the second peak is changed, we deduce that the origin and mechanism of this phenomenon may be related to both the frequency of ac field and the magnitude of bias field. To validate this, ME voltage versus the driving frequency is plotted in Fig.10,with the bias field being set at previously specified ones,i.e.,560,680,788,and 896 Oe, which are the values of the optimal fields at frequencies of 400, 410, 420, and 430 kHz, respectively. When the bias field is fixed, the ME voltage maximizes at a certain frequency due to the mechanical resonance of the ring-shaped composites under ac magnetic field. This indicates a one-toone correspondence between the resonance frequency and the bias field we call it the required field. It should be noted that the resonance frequency and the required field in Fig. 10 coincide exactly with those in Fig.9(b). Meanwhile,the second peak of ME voltage in Fig.9(b)at a fixed bias field are exactly the same as the peak ME voltage in Fig. 10 at the same bias field. Therefore,it has been proved that the occurrence of the second peak is caused by the resonance under a certain value of magnetic field. Equation (16) reveals that variation of the equivalent elastic coefficient with external bias field,known as?Eeffect,[51,52]may lead to the resonance frequency change with the bias field. That is why the peak ME voltage and the resonance frequency in Fig. 10 changes with the bias field.Based on this mechanism,a much broader range of resonance frequencies is expected to be achieved.

    Fig.9. ME voltage versus the applied bias magnetic field for different frequencies:(a)100,130,160,and 190 kHz,and(b)400,410,420,and 430 kHz(σ0=?10 MPa).

    Fig. 10. ME voltage versus the driving frequency at a bias magnetic field of 560,680,788,and 896 Oe.

    Figure 11 displays the ME voltage versus the applied bias magnetic field for different pre-stresses. As pre-stress increases,the first peak increases with the optimal field moving slightly toward right, while the second peak reduces with the required field moving toward right. Thereby, available ME coupling may be realized within a wide bias field. By way of example,the ME voltage withσ0=?10 MPa reaches 70%of the second peak when the bias field is in the range from 156 to 560 Oe. The variation of the optimal(required)fields and the peaks are plotted in Fig. 12. It can be seen that each keeps a linear relation with the external pre-stress,offering much convenience for tuning ME performance of the proposed energy harvester.

    Fig.11. ME voltage versus the applied bias magnetic field for different pre-stresses(f =400 kHz).

    Fig.12. The peak values and the corresponding optimal bias fields versus pre-stress:(a)the first peak,and(b)the second peak(f =400 kHz).

    Lastly,we study the effects of pre-stress on the resonance of the composites. Figure 13 shows the ME voltage versus the driving frequency for different pre-stresses. At a lower bias field(100 Oe),ME voltage is enhanced by increasing tensile stress, with the resonance frequency being reduced. On the contrary,the ME voltage is weakened by increasing compressive stress,with the resonance frequency being increased.However,at a high bias field(500 Oe),the ME voltage is improved via reducing tensile stress or increasing compressive stress, with the resonance frequency being reduced. This indicates that a critical bias field may occur, below and above which the effects of pre-stress on the peak ME voltage and the resonance frequency are completely opposite. This critical bias field is actually the same as the optimal field corresponding to extreme point of the equivalent coefficient of Terfenol-D.

    Fig. 13. ME voltage versus the driving frequency for different prestresses: (a)Hdc=100 Oe,and(b)Hdc=500 Oe.

    Figure 14 shows the contour plots of the peakME v oltage and resonance frequency as functions of pre-stress and bias field. The results presented further validate the different change laws of peak ME voltage and resonance frequency with the pre-stress at different bias magnetic fields. When the bias field is lower than 200 Oe,the peak ME voltage is reduced as pre-stress increases from?60 to 0 MPa. However, a nonmonotone variation is observed when the bias field exceeds 200 Oe. As a consequence, the maximum peak occurs in the region of (?30 MPa<σ0240 Oe). Figure 14(b)indicates that a broad range of resonance frequencies is obtained to be 385–425 kHz by coordinately control external pre-stress and magnetic stimuli.It is expected that the range of resonance frequencies can be further broadened by adjusting bias magnetic field. The critical bias filed for resonance frequency is about 275 Oe,below which the latter monotonically increases with pre-stress, or above which firstly reduces and then increases. Therefore, the predicted nonlinear behaviors under external stimuli build a new theoretical mechanism for tunable control of ME energy harvesters.

    Fig.14. Contour plots of(a)the peak ME voltage and(b)the resonance frequency as functions of pre-stress and bias field.

    5. Conclusions

    In summary, we have proposed a theoretical model for studying ME coupling of energy harvesters using ring-shaped composites, with introducing a new type of boundary conditions as well as considering the multi-field coupled material’s properties. This model can be applied to predict output ME voltage, resonance frequency and power density of applications including energy harvesting and sensing under external stimuli. By utilizing the theoretical model, some important conclusions can be drawn as follows:

    (1)Setting a stress-free condition at any boundary is beneficial for inducing mechanical resonance,and thereby the F-F boundary condition corresponds to the lowest resonance frequency and more harmonics within the range of frequencies(0–1400 kHz)investigated.

    (2)Clamped conditions on the inner and outer boundaries may give rise to the stress built up in the composite,and thus improve the maximum power density. Once the composite is constrained at least on one side,the optimal resistance reduces and the resonance frequency increases.

    (3)Stress boundary conditions strengthen strain distribution, and enhance ME coupling significantly, without changing the corresponding resonance frequency and optimal resistance. The applied stress in actual operation may be properly increased on the premise that its valuePshould not exceed the allowable stresses of constituent materials.

    (4) A theoretically revealed dual-peak phenomenon around resonance frequencies demonstrates the occurrence of two optimal bias fields corresponding to the strong ME coupling. The linear relation between pre-stress and the optimal bias field offering much convenience for tuning ME performance within a wide field range.

    (5) The peak ME voltage is increased from 1.4 to 3.2 V with applying proper pre-stress and bias field. A broad range of resonance frequencies can be obtained from 385 to 425 kHz by tuning pre-stress from 0 to 60 MPa and bias field from 100 to 400 Oe. It is expected that the mutual coordination of the applied bias magnetic field and pre-stress promotes commonly the improvement of ME performance.

    The present theoretical model provides more options for optimizing ME performance of ring-shaped composites via mechanical and magnetic conditions. Nonlinear coupling properties revealed in this paper may result in an emerging revolution for tunable, active, or even smart control of ME harvesters. In addition,the current modeling idea may be applied to solve nonlinear multi-field coupling problem of other smart materials and structures.

    猜你喜歡
    楊勇
    楊勇
    Activated dissociation of H2 on the Cu(001)surface:The role of quantum tunneling
    楊勇書(shū)法作品
    大型房企違約對(duì)金融部門(mén)的影響
    楊勇:我們是人民的勤務(wù)兵,不是官老爺
    Angular dependence of vertical force and torque when magnetic dipole moves vertically above flat high-temperature superconductor?
    Quantum nature of proton transferring across one-dimensional potential fields?
    凡人光陰
    牡丹(2021年5期)2021-03-19 13:28:03
    故鄉(xiāng)
    牡丹(2020年3期)2020-03-02 02:20:17
    家庭百事通(2017年9期)2017-09-11 07:28:52
    80岁老熟妇乱子伦牲交| 丰满人妻熟妇乱又伦精品不卡| 18禁黄网站禁片午夜丰满| 后天国语完整版免费观看| 麻豆av在线久日| 天堂中文最新版在线下载| 欧美黄色片欧美黄色片| 欧美老熟妇乱子伦牲交| 国产精品麻豆人妻色哟哟久久| av天堂在线播放| 啦啦啦啦在线视频资源| 国产精品一区二区精品视频观看| 99九九在线精品视频| 99九九在线精品视频| 夜夜骑夜夜射夜夜干| 一区二区三区四区激情视频| 亚洲欧美一区二区三区久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲自偷自拍图片 自拍| 夜夜骑夜夜射夜夜干| 在线亚洲精品国产二区图片欧美| 成年人免费黄色播放视频| 交换朋友夫妻互换小说| 99热网站在线观看| 男男h啪啪无遮挡| 在线观看www视频免费| 男人添女人高潮全过程视频| 天堂中文最新版在线下载| 国产麻豆69| 成人免费观看视频高清| 日韩伦理黄色片| 欧美精品av麻豆av| 国产精品免费视频内射| 亚洲精品一二三| cao死你这个sao货| 亚洲精品国产一区二区精华液| 99九九在线精品视频| 纯流量卡能插随身wifi吗| 欧美激情 高清一区二区三区| 欧美在线一区亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩国产mv在线观看视频| 欧美日韩亚洲高清精品| 国产精品九九99| 成人手机av| 亚洲欧洲精品一区二区精品久久久| 欧美 亚洲 国产 日韩一| 99久久人妻综合| 啦啦啦 在线观看视频| 嫁个100分男人电影在线观看 | 精品人妻1区二区| 欧美激情极品国产一区二区三区| 国产成人精品久久久久久| 国产亚洲精品久久久久5区| 19禁男女啪啪无遮挡网站| 国产成人一区二区在线| 午夜影院在线不卡| 97在线人人人人妻| 午夜福利视频在线观看免费| 国产高清视频在线播放一区 | 国产精品久久久久久精品电影小说| 国产激情久久老熟女| 亚洲国产欧美在线一区| 国产欧美日韩综合在线一区二区| 麻豆国产av国片精品| 国产黄色免费在线视频| 久久精品国产亚洲av高清一级| 男女床上黄色一级片免费看| 国产又色又爽无遮挡免| 久久毛片免费看一区二区三区| 十八禁网站网址无遮挡| 国产三级黄色录像| 午夜福利免费观看在线| 国产亚洲欧美在线一区二区| 成人亚洲精品一区在线观看| 精品国产一区二区三区久久久樱花| 国产精品av久久久久免费| 男女下面插进去视频免费观看| 18在线观看网站| 欧美日韩一级在线毛片| 国产男女内射视频| 丰满迷人的少妇在线观看| 一边摸一边做爽爽视频免费| 国产一区二区三区综合在线观看| 亚洲熟女精品中文字幕| 99热国产这里只有精品6| 伊人亚洲综合成人网| 亚洲自偷自拍图片 自拍| 亚洲av综合色区一区| 午夜影院在线不卡| 国产高清不卡午夜福利| 亚洲天堂av无毛| 18禁裸乳无遮挡动漫免费视频| av一本久久久久| 免费一级毛片在线播放高清视频 | 久久国产精品男人的天堂亚洲| 精品人妻熟女毛片av久久网站| 精品福利永久在线观看| 亚洲国产av影院在线观看| 一级毛片电影观看| 青青草视频在线视频观看| 人人澡人人妻人| 精品少妇内射三级| 美国免费a级毛片| 最近手机中文字幕大全| 亚洲中文字幕日韩| 一区二区日韩欧美中文字幕| 一级毛片黄色毛片免费观看视频| 久久精品久久久久久噜噜老黄| 国产高清videossex| 不卡av一区二区三区| 天天躁夜夜躁狠狠躁躁| 1024香蕉在线观看| 久久99一区二区三区| 国产成人av教育| 男女高潮啪啪啪动态图| 人妻人人澡人人爽人人| 国产在线一区二区三区精| 久热这里只有精品99| 国产男女超爽视频在线观看| 永久免费av网站大全| 免费看十八禁软件| 十八禁人妻一区二区| 欧美成狂野欧美在线观看| 女人精品久久久久毛片| 久久久久视频综合| 国产成人免费无遮挡视频| 亚洲欧美一区二区三区国产| 国产成人欧美| 午夜福利乱码中文字幕| 欧美日本中文国产一区发布| 人体艺术视频欧美日本| 男女免费视频国产| 久久综合国产亚洲精品| 亚洲综合色网址| 国产成人av激情在线播放| 欧美日韩视频高清一区二区三区二| 精品国产乱码久久久久久男人| 91国产中文字幕| 777米奇影视久久| 我的亚洲天堂| 午夜激情av网站| 男女午夜视频在线观看| 捣出白浆h1v1| 亚洲一区二区三区欧美精品| 国产欧美日韩综合在线一区二区| 久久久久精品国产欧美久久久 | 欧美老熟妇乱子伦牲交| 美国免费a级毛片| 久热爱精品视频在线9| 国产高清视频在线播放一区 | 欧美av亚洲av综合av国产av| 国产在线免费精品| 免费看十八禁软件| 又粗又硬又长又爽又黄的视频| 少妇猛男粗大的猛烈进出视频| 国产精品一二三区在线看| 日韩视频在线欧美| 亚洲人成77777在线视频| 亚洲美女黄色视频免费看| 夫妻午夜视频| 欧美 亚洲 国产 日韩一| 国产高清视频在线播放一区 | 黄片小视频在线播放| 久久精品国产亚洲av涩爱| 中国美女看黄片| 亚洲综合色网址| 亚洲欧美清纯卡通| 成人国产一区最新在线观看 | 精品少妇久久久久久888优播| 狂野欧美激情性bbbbbb| 天天影视国产精品| 国产精品一区二区在线观看99| 99香蕉大伊视频| 久久久久久久久免费视频了| 精品少妇一区二区三区视频日本电影| 国产真人三级小视频在线观看| 精品国产一区二区三区久久久樱花| 精品国产乱码久久久久久小说| 制服人妻中文乱码| 水蜜桃什么品种好| 亚洲国产看品久久| 国产熟女欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品av麻豆狂野| 蜜桃国产av成人99| 午夜福利乱码中文字幕| 性色av乱码一区二区三区2| 日本wwww免费看| 亚洲色图综合在线观看| 丝袜脚勾引网站| 免费少妇av软件| 国产成人精品在线电影| 亚洲欧美成人综合另类久久久| 脱女人内裤的视频| 一区二区三区激情视频| 国产精品久久久人人做人人爽| 久久人妻熟女aⅴ| 久久精品亚洲av国产电影网| 熟女少妇亚洲综合色aaa.| av在线播放精品| 一二三四在线观看免费中文在| 久热这里只有精品99| cao死你这个sao货| 2018国产大陆天天弄谢| 国产精品久久久久成人av| 韩国精品一区二区三区| 黄色视频不卡| 菩萨蛮人人尽说江南好唐韦庄| 精品国产国语对白av| 国产精品一二三区在线看| 国产主播在线观看一区二区 | 水蜜桃什么品种好| 国产不卡av网站在线观看| 免费高清在线观看视频在线观看| 一级毛片女人18水好多 | 黄色片一级片一级黄色片| 亚洲国产精品一区三区| 777久久人妻少妇嫩草av网站| 人人妻,人人澡人人爽秒播 | 精品免费久久久久久久清纯 | 免费av中文字幕在线| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 亚洲国产欧美日韩在线播放| 两个人免费观看高清视频| 久久女婷五月综合色啪小说| 如日韩欧美国产精品一区二区三区| 夫妻午夜视频| av电影中文网址| 黄网站色视频无遮挡免费观看| av欧美777| 国产女主播在线喷水免费视频网站| 99热国产这里只有精品6| 老司机亚洲免费影院| 熟女av电影| 亚洲色图 男人天堂 中文字幕| 极品人妻少妇av视频| 亚洲男人天堂网一区| 国产av一区二区精品久久| 18在线观看网站| 亚洲精品国产色婷婷电影| 热99国产精品久久久久久7| 亚洲专区中文字幕在线| 国产精品一国产av| 日本wwww免费看| 真人做人爱边吃奶动态| 国产亚洲一区二区精品| 国产野战对白在线观看| 一二三四社区在线视频社区8| 亚洲人成电影观看| 亚洲精品国产一区二区精华液| 国产片特级美女逼逼视频| 女人被躁到高潮嗷嗷叫费观| 十八禁网站网址无遮挡| 国产黄频视频在线观看| www.精华液| 手机成人av网站| 热99久久久久精品小说推荐| 美女大奶头黄色视频| 国产男女内射视频| 又黄又粗又硬又大视频| 99九九在线精品视频| 免费看十八禁软件| 99精品久久久久人妻精品| 亚洲国产欧美在线一区| 人妻人人澡人人爽人人| 涩涩av久久男人的天堂| 捣出白浆h1v1| 最新在线观看一区二区三区 | 精品高清国产在线一区| 99久久人妻综合| 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 午夜老司机福利片| 男女国产视频网站| 欧美乱码精品一区二区三区| 国产在线一区二区三区精| 欧美老熟妇乱子伦牲交| 母亲3免费完整高清在线观看| 国产成人欧美在线观看 | 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美一区二区综合| 精品一区在线观看国产| 亚洲五月色婷婷综合| 午夜免费男女啪啪视频观看| 欧美精品人与动牲交sv欧美| 男女国产视频网站| 美女视频免费永久观看网站| 欧美日本中文国产一区发布| 日韩av在线免费看完整版不卡| 国产欧美日韩综合在线一区二区| 精品少妇久久久久久888优播| 在线观看免费高清a一片| 每晚都被弄得嗷嗷叫到高潮| 午夜福利一区二区在线看| 亚洲欧美一区二区三区黑人| 久久久久国产一级毛片高清牌| 国产av国产精品国产| 欧美日韩亚洲高清精品| 免费在线观看完整版高清| 欧美日韩成人在线一区二区| 中国国产av一级| 亚洲欧美一区二区三区国产| 亚洲av美国av| 男人操女人黄网站| 制服人妻中文乱码| 精品少妇一区二区三区视频日本电影| 午夜免费观看性视频| 一级片免费观看大全| 国产亚洲午夜精品一区二区久久| 狠狠精品人妻久久久久久综合| 中文字幕人妻丝袜制服| www.精华液| 极品少妇高潮喷水抽搐| 别揉我奶头~嗯~啊~动态视频 | 少妇被粗大的猛进出69影院| 多毛熟女@视频| av在线老鸭窝| 国产97色在线日韩免费| 97在线人人人人妻| 少妇人妻久久综合中文| 国产欧美日韩一区二区三区在线| 成人黄色视频免费在线看| 午夜激情av网站| 男女之事视频高清在线观看 | 亚洲男人天堂网一区| 国产激情久久老熟女| 一个人免费看片子| 久久久久久亚洲精品国产蜜桃av| 精品人妻1区二区| 欧美日韩福利视频一区二区| 欧美日韩黄片免| 日韩av在线免费看完整版不卡| 五月天丁香电影| 久久人妻福利社区极品人妻图片 | 欧美xxⅹ黑人| 免费在线观看视频国产中文字幕亚洲 | 秋霞在线观看毛片| 最新的欧美精品一区二区| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 99国产综合亚洲精品| 欧美黄色淫秽网站| 天堂8中文在线网| 亚洲精品国产av蜜桃| kizo精华| 女警被强在线播放| av电影中文网址| 日韩欧美一区视频在线观看| 免费高清在线观看视频在线观看| 校园人妻丝袜中文字幕| 狠狠精品人妻久久久久久综合| 国产三级黄色录像| 欧美人与性动交α欧美精品济南到| 99国产精品99久久久久| 99久久人妻综合| 男女之事视频高清在线观看 | 水蜜桃什么品种好| 亚洲国产精品999| 久久天躁狠狠躁夜夜2o2o | 涩涩av久久男人的天堂| 免费看十八禁软件| 午夜日韩欧美国产| 午夜免费观看性视频| 久久久久国产一级毛片高清牌| 丁香六月天网| 各种免费的搞黄视频| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 飞空精品影院首页| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 中文欧美无线码| www.精华液| 国产xxxxx性猛交| 亚洲成国产人片在线观看| av天堂在线播放| 菩萨蛮人人尽说江南好唐韦庄| 丝瓜视频免费看黄片| 人人妻人人澡人人看| 999久久久国产精品视频| 在线观看免费午夜福利视频| 亚洲图色成人| 女人精品久久久久毛片| 波多野结衣av一区二区av| 桃花免费在线播放| 亚洲第一av免费看| 啦啦啦在线观看免费高清www| 午夜影院在线不卡| 亚洲国产中文字幕在线视频| 考比视频在线观看| 在线观看国产h片| 国产91精品成人一区二区三区 | 亚洲成国产人片在线观看| 久久久国产欧美日韩av| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频| 爱豆传媒免费全集在线观看| 国产精品香港三级国产av潘金莲 | 精品视频人人做人人爽| 亚洲天堂av无毛| 一区二区三区四区激情视频| 国产精品久久久久久精品古装| 国产有黄有色有爽视频| 亚洲情色 制服丝袜| 嫩草影视91久久| 久久99一区二区三区| 欧美久久黑人一区二区| 国产爽快片一区二区三区| 天堂中文最新版在线下载| 国产亚洲午夜精品一区二区久久| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产色婷婷电影| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| 国产1区2区3区精品| 搡老岳熟女国产| 日韩伦理黄色片| 亚洲自偷自拍图片 自拍| 国产在线视频一区二区| 亚洲国产中文字幕在线视频| 亚洲精品一区蜜桃| 两性夫妻黄色片| 午夜激情av网站| 免费观看av网站的网址| 色视频在线一区二区三区| 久久人人97超碰香蕉20202| 国产在线视频一区二区| 国产亚洲欧美精品永久| 久久 成人 亚洲| 妹子高潮喷水视频| 欧美激情高清一区二区三区| 国产精品三级大全| 亚洲少妇的诱惑av| 亚洲熟女精品中文字幕| 午夜激情av网站| 国产有黄有色有爽视频| 亚洲黑人精品在线| 午夜91福利影院| 成年美女黄网站色视频大全免费| 九色亚洲精品在线播放| 最新的欧美精品一区二区| 婷婷色av中文字幕| 日韩制服骚丝袜av| 男人添女人高潮全过程视频| 欧美少妇被猛烈插入视频| 欧美黄色片欧美黄色片| 在线观看免费高清a一片| 欧美日韩亚洲综合一区二区三区_| 超碰成人久久| 亚洲精品日韩在线中文字幕| 女人被躁到高潮嗷嗷叫费观| 国产97色在线日韩免费| 久久精品aⅴ一区二区三区四区| 青春草亚洲视频在线观看| 国产欧美亚洲国产| 91九色精品人成在线观看| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡 | 婷婷色综合大香蕉| 涩涩av久久男人的天堂| 国产精品 欧美亚洲| 国产伦理片在线播放av一区| 精品高清国产在线一区| 国产福利在线免费观看视频| 午夜两性在线视频| 国产极品粉嫩免费观看在线| 亚洲第一av免费看| 久久久久精品国产欧美久久久 | 99九九在线精品视频| 夜夜骑夜夜射夜夜干| 在线精品无人区一区二区三| 国产精品 欧美亚洲| 日本欧美视频一区| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| 性色av乱码一区二区三区2| 精品国产国语对白av| 亚洲成人国产一区在线观看 | 狠狠精品人妻久久久久久综合| 亚洲美女黄色视频免费看| 一级黄色大片毛片| 亚洲av欧美aⅴ国产| 国产精品一区二区精品视频观看| 99国产精品一区二区蜜桃av | 老司机在亚洲福利影院| 久久人妻福利社区极品人妻图片 | 女性生殖器流出的白浆| 黑人欧美特级aaaaaa片| 国产主播在线观看一区二区 | 精品人妻1区二区| 免费观看a级毛片全部| 秋霞在线观看毛片| 亚洲人成电影观看| 亚洲精品在线美女| 18在线观看网站| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 91九色精品人成在线观看| 岛国毛片在线播放| 99热网站在线观看| 无遮挡黄片免费观看| 亚洲免费av在线视频| 91成人精品电影| 操美女的视频在线观看| 婷婷色av中文字幕| 欧美成狂野欧美在线观看| 亚洲专区国产一区二区| 国产成人影院久久av| 亚洲精品第二区| av网站免费在线观看视频| 国产精品99久久99久久久不卡| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产有黄有色有爽视频| 好男人视频免费观看在线| 啦啦啦 在线观看视频| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 亚洲精品久久午夜乱码| 老汉色av国产亚洲站长工具| 啦啦啦视频在线资源免费观看| 另类亚洲欧美激情| 国产片内射在线| 久久av网站| 国产高清国产精品国产三级| a级片在线免费高清观看视频| 国产精品 欧美亚洲| 男人添女人高潮全过程视频| 国产片内射在线| 亚洲,欧美精品.| 久久综合国产亚洲精品| 国产一区二区三区综合在线观看| 一级毛片女人18水好多 | 久久久久久久精品精品| 亚洲精品美女久久av网站| 久久精品久久精品一区二区三区| 久久鲁丝午夜福利片| 咕卡用的链子| 男女之事视频高清在线观看 | 国产亚洲精品第一综合不卡| 日本vs欧美在线观看视频| 精品亚洲成国产av| 国产麻豆69| 老熟女久久久| 少妇人妻久久综合中文| 久久精品国产a三级三级三级| 国产麻豆69| 国产亚洲精品第一综合不卡| 777米奇影视久久| 国产又色又爽无遮挡免| 777米奇影视久久| avwww免费| 成人18禁高潮啪啪吃奶动态图| 久久精品aⅴ一区二区三区四区| 国产激情久久老熟女| 亚洲久久久国产精品| 一区二区三区四区激情视频| 岛国毛片在线播放| 国产成人a∨麻豆精品| 久久中文字幕一级| 无限看片的www在线观看| 亚洲国产欧美一区二区综合| 久久国产精品影院| 视频在线观看一区二区三区| 免费在线观看影片大全网站 | av国产久精品久网站免费入址| a级片在线免费高清观看视频| av又黄又爽大尺度在线免费看| 久久狼人影院| 日韩av免费高清视频| 精品亚洲乱码少妇综合久久| av福利片在线| 黑人巨大精品欧美一区二区蜜桃| 国产成人91sexporn| 国产日韩欧美在线精品| 久久午夜综合久久蜜桃| 视频区欧美日本亚洲| 在现免费观看毛片| 久久性视频一级片| 亚洲国产成人一精品久久久| 少妇人妻 视频| 巨乳人妻的诱惑在线观看| 国产人伦9x9x在线观看| 国产精品.久久久| 手机成人av网站| 亚洲熟女精品中文字幕| 91成人精品电影| 精品一区在线观看国产| 看十八女毛片水多多多| av不卡在线播放| 欧美日韩成人在线一区二区| 国产深夜福利视频在线观看| 日本vs欧美在线观看视频| 欧美变态另类bdsm刘玥| 国产精品欧美亚洲77777| 美女午夜性视频免费| 丰满人妻熟妇乱又伦精品不卡| 在线天堂中文资源库| 一本大道久久a久久精品| 电影成人av| 韩国高清视频一区二区三区| 少妇人妻 视频| 国产淫语在线视频| 日本av手机在线免费观看| 99精国产麻豆久久婷婷| 精品福利观看| 波多野结衣一区麻豆| 亚洲伊人色综图| 男女边吃奶边做爰视频| 日韩,欧美,国产一区二区三区| 亚洲 欧美一区二区三区|