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

    Advanced Frozen Permeability Technique and Applications in Developing High Performance Electrical Machines

    2016-11-17 07:33:45ZhuChu
    電工技術學報 2016年20期
    關鍵詞:謝菲爾德磁導率端電壓

    Zhu Z Q Chu W Q

    ?

    Advanced Frozen Permeability Technique and Applications in Developing High Performance Electrical Machines

    (University of Sheffield Sheffield S1 3JD U.K.)

    Magnetic saturation is one of the major challenges for all electrical machines and makes it much more difficult to design high performance electrical machines. Frozen permeability (FP) technique has recently been developed to predict the on-load components of electromagnetic qualities while considering the influence of magnetic saturation. In this paper, the FP technique and its applications in various electrical machines are comprehensively reviewed. It shows that the FP technique is able to accurately separate on-load magnetic field, flux linkage, and inductance components in all types of electrical machines. These on-load field components and parameters can be further utilized to investigate various on-load performance, such as on-load torque and torque ripple components, on-load back EMFs and voltages, flux weakening performance, radial force distribution, et al. The accurate on-load parameters can also be used to improve the machine models and controls. Hence, FP technique provides a holistic approach to the development of high performance electrical machines.

    Electrical machines, frozen permeability, on-load performance, permanent magnet

    0 Introduction

    Due to global energy security and environmental concerns, significant efforts have been made on exploiting renewable and sustainable resources of energy, such as wind power, solar power, and tidal power, etc. Electrical machines and drives represent one of the key enabling technologies for exploiting these energy resources. Meanwhile, electrical drives also exhibit the advantages of high efficiency, fast response and accuracy. These advantages coupled with recent advances in computational capabilities and power electronics have accelerated the pace of replacing conventional mechanical, hydraulic, or pneumatic propulsion and transmission systems by electrical systems for various applications[1-4].

    In electric drive systems, the system performance largely depends on the electrical machines employed. Progress in the materials, power electronics, and control technologies over the past few decades has enabled the developments of various electrical machines[1-12]. The major electrical machine techno- logies are summarized in Fig.1. Historically, electric drives were dominated by induction machines (IM) and electrically excited (EE) machines. However, more recently, new machine technologies such as permanent magnet (PM) machines have become more popular due to high torque density and efficiency, especially when high energy rare-earth permanent magnet materials are employed. Switched reluctance (SR) and synchronous reluctance machines are developed by utilizing reluctance torque. The torque densities of reluctance machines are similar to IM's and hence reluctance machines are usually considered as alternative options to IMs, albeit with a low power factor. However, electronic controllers are inevitably required to operate reluctance machines. Hybrid electrical machines are developed by uniquely integrating different machine technologies. Due to the variety ways of integration, a large number of hybrid electrical machines have been developed and are one of the most popular research topics under investi- gation[12].

    Despite of variety technologies, one of the major challenges for all electrical machines is the magnetic saturation. This is due to the nature of nonlinear soft magnetic materials. The magnetic saturation affects nearly all aspects of electrical machines, such as magnetic field distributions, flux linkages, voltages, average torque and torque ripple, overload capability, inductances, flux weakening capability, torque speed curve, iron losses and efficiency. More importantly, the magnetic saturation makes it much more difficult to predict the breakdown contributions and/or components of aforementioned performances. However, in order to develop high performance electrical machines, it will be greatly helpful to reveal the relationships between the design parameters and the breakdown contributions/components of aforemen- tioned performance under different loadings.

    (a) Classification

    (b) Induction (c) Electrically excited

    (d) Permanent magnet (e) Switched relectance

    (f) Synchronous reluctance (g) Hybrid

    Fig.1 Classification and cross-sections of different electrical machines

    On the other hand, with the progress of finite element software, frozen permeability (FP) technique has been developed and becomes available in the commercial software. The FP technique is able to calculate the electromagnetic qualities while con- sidering the influence of magnetic saturation and hence widely used to investigate the issues related with magnetic saturation. Thus, it can provide strong technical supports on developing high performance electrical machines.

    In this paper, the advanced FP technique and its applications in various electrical machines are comprehensively reviewed. The principle of FP technique is described in section 1. The sections 2~9 are devoted to the applications of FP technique on different topics of PM machines, such as on-load magnetic field, flux linkages and inductance estimation, separation of torque components, on-load cogging torque, on-load torque ripple, on-load back EMF, on-load voltage waveforms, flux weakening et al. The applications of FP technique to other technical issues and other types of electrical machines are highlighted in sections 10 and 11.

    1 Principle of Frozen Permeability

    The early publications on FP technique can be found in Ref.[13-16]. Currently, FP technique has been widely used in all types of electrical machines to investigate the issues related with magnetic saturation and cross-coupling as will be shown later. Taking PM machines as an example, the principle of FP technique is illustrated in Fig.2.

    Fig.2 Principle of frozen permeability technique

    When a PM machine is on load, the machine is excited by both the PM and the armature current. The resultant operating point can be presented by point 1 (allandall). It is well known that due to nonlinear-curve of soft magnetic material, althoughall=PM+i,all≠PM+i, wherePMrepresents the resultant flux density with the PM excitation (PM) only (point 2) andirepresents the resultant flux density with the armature excitation (i) only (point 3). Therefore, the on-load field components due to PM and armature cannot be accurately decomposed by using conventional calculations with the PM or armature excitation only.

    The FP technique is based on the fact that all excitations always share the same permeability distribution when electrical machines are on load. Hence, the on-load permeabilityallis obtained and stored by solving the on-load model (point 1). Then, two linear analyses represented by the points 4 and 5 are solved further based on the frozenallwith either the PM or armature excitation only to obtain(FP,PM) or(FPi), respectively. Since all three calculations are based on the same permeability (all),all=(FP,PM)+(FP,i). Therefore, the on-load PM and armature field components are decomposed. Further- more, sinceallvaries with the load conditions, the influence of magnetic saturation and cross-coupling has been inherently taken into account.

    Based on the principle of FP technique, its procedure can be summarized in Fig.3. Firstly, the nonlinear calculation with all excitations, referred as the whole model, is solved. Secondly, the permea- bility of the whole model, referred as the on-load permeability, in each element and at every rotor position is saved and frozen. Finally, linear calcula- tions are conducted based on the on-load permeability to obtain the on-load magnetic field components of the selected excitations.

    Fig.3 Procedure of frozen permeability technique

    Although the illustration in Fig.2 is based on a PM machine, the FP technique is applicable to all types of electrical machines. Also the selected excitations for the linear calculations can be tailored according to investigated topics. For example, to obtain the on-load PM or armature magnetic field component in PM machines, only the PM or armature excitation is applied in the linear calculation[13,16]. To study the cross-coupling inductances between the d- and q-axes in PM machines, only the d- or q-axis current is applied[13,16]. To investigate the mutual coupling between phases in SR machines, only one phase current is applied[17].

    The major challenge for the FP technique is the experimental validation, since the exact quantities that are calculated in the FP technique cannot be measured directly. Although it is tried in Ref.[18], various further work on experimental validation needs to be done. However, on the other hand, the FP technique has been widely accepted based on the fundamental magnetic circuit theory and used for on-load magnetic field separation and parameter predictions in various types of electrical machines.

    In following sections, the applications of the FP technique are comprehensively reviewed. Since the most of the applications are on different technical topics of PM machines, only one of following sections is focused on the applications on the other types of machines and will presented at the end. For the rest sections, the publications are grouped according to different research topics, such as on-load magnetic field separation, on-load flux linkages and inductance calculation, torque component separation, et al.

    2 On-Load Magnetic Field Separation

    On-load magnetic field separation is the basic purpose of the FP technique and the foundation for all further investigations. Nearly all investigations using FP technique involve the on-load magnetic field separation, although the on-load magnetic field results may not be presented.

    Fig.4 and Fig.5 show the on-load magnetic field components of an inset PM machines presented in Ref.[19] as an example. The peak phase currentais 4A, the current phase advance angle to the positive q-axis0°and the rotor positionis 0°(electrical). The main parameters of the inset PM machine are listed in Tab.1. It can be seen that the on-load magnetic field and permeability distributions in Fig.4 are neither symmetrical with the d-axis nor with the q-axis. This is due to the influence of armature field. Consequently, the on-load PM magnetic field with the PM excitation on the d-axis has the net q-axis flux linkage. The q-axis current also produces net d-axis flux linkage. Fig.5 shows the radial and tangential flux density components along the middle of the airgap in Fig.4 to illustrate the on-load magnetic flux density separation. It shows that for each sample point, the flux density in the whole model is the same as the mathematic summation of the two flux density components produced with either on-load PM or armature field only. This is also true for any other rotor position and load condition. Hence, the on-load magnetic field components can be successfully separated while the influence of magnetic saturation and cross-coupling has been taken into account.

    (a) Whole field (PM and current) (b) On-load permeability

    (c) On-load PM field (FP) (d) On-load armature field (FP)

    Fig.4 On-load field distributions of inset PM machine when Ia=4A, ?=0° and ??=0°[19]

    (a) Radial flux density components

    (b) Tangential flux density components

    Fig.5 Flux density components at middle of airgap of inset PM machine when Ia=4A, =0°and =0°[19]

    Tab.1 Main parameters of inset PM machine in Ref.[19]

    Other examples of separated on-load magnetic field distributions and flux densities in PM machines can be found in Ref.[20-34]. In Ref.[20], the on-load PM and armature magnetic field distributions and airgap flux densities in a large power linear PM machines are obtained to explain the reduction of trust force density. In Ref.[23], the on-load airgap PM and armature flux density components in a surface- mounted PM (SPM) machine with eccentricity are obtained by the FP finite element analyses to validate the analytical model based on improved conformal mapping method.

    Based on the on-load magnetic field separation, various on-load components of electromagnetic quantities such as flux densities and flux linkages can then be calculated for further investigations.

    3 On-Load Flux Linkages and Inductances Analyses

    For PM machines, the PM flux linkages and inductances are most important parameters in terms of the machine design, modelling and control. The inductances represent the armature flux linkages with the given armature currents. However, it is well known that the constant PM flux linkage and indu- ctance models may result in significant inaccuracy due to the influence of magnetic saturation. Since the FP technique is able to consider the influence of magnetic saturation, it is widely used for the calcu- lation of PM flux linkages and inductances.

    In Ref.[13,16,19, 24], it has been well illustrated that the average PM flux linkage and inductances per electric cycle can significantly change with load. Consequently, they are 2D functions of both d- and q-axis armature currents (or the current amplitudeaand the phase advance angle). Furthermore, there are cross-coupling between the d- and q-axes as well as the net average q-axis PM flux linkage. Hence, the full average flux linkage model for PM machines should be expressed as

    wheredandqare the total d- and q-axis flux linkages, respectively.d(PM),q(PM),d(i) andq(i) are the d- and q-axis flux linkages due to on-load PM and armature fields, respectively.dandqare the d- and q-axis currents, respectively.dd,qq, anddqare the d- and q-axis self- and mutual inductances, respectively.

    Please note thatq(PM) anddqare due to the cross-coupling between d- and q-axis. Furthermore, except the variation of average values with the current amplitudeaand phase advance angle, the flux linkages and inductances also change with the rotor position. Fig.6 shows the variations of on-load flux linkages of the inset PM machine listed in Tab.1 with the rotor position whena=4A and 0° as an example. It also confirms that not only the instant- aneousq(PM) but also its average are not zero. Theaveraged(i) is also not zero, although only q-axis current is applied.

    (a) d-axis flux linkages

    (b) q-axis flux linkages

    Fig.6 Variation of flux linkages with rotor position of inset PM machine when Ia=4A and ?=0°[19]

    The variation of average flux linkages and inductances withis further presented in Fig.7. The influence of magnetic saturation is represented by the variation ofd(PM),dd, andqq. When=90°, the d-axis magnetic saturation is the heaviest. Hence,d(PM) andddare at their minimum. When=0°,dis zero andqis maximum. The magnetic saturation is modest along the d-axis but the heaviest along the q-axis. Therefore,d(PM) andddare higher whileqqis the lowest. When=90°,dis at the negative maximum and the magnetic saturation along the d-axis is the lowest due to the maximum flux weakening.Therefore,d(PM) andddare at the maximum. The influence of cross-coupling is represented by the variation ofq(PM) anddq. The cross-coupling is zero and henceq(PM) anddqare zero when=±90°, since all excitations are applied to the d-axis. Althoughqis maximum when=0°, the cross-coupling is the heaviest when=30°. It is due to that the cross-coupling is aggravated by the magnetic saturation.

    (a) d-axis flux linkages

    (b) q-axis flux linkages

    (c) Inductances

    Fig.7 Variation of average flux linkages and inductances of inset PM machine with ? when Ia=4A[19]

    Some of the other examples of on-load flux linkage and inductance calculation using FP technique can be found in Ref.[15, 24, 31]. In Ref.[15], the d- and q-axis inductances are predicted based on the combination of FP and stored magnetic energy method to improve accuracy of torque estimation in a 36-slot/4-pole (36s/4p) IPM machine. In Ref.[24], the on-load PM flux linkages and inductances are calculated for the output torque estimation of a 30s/36p inset PM machine. It shows the torque variation with the current phase advance angle can be predicted more accurately when the cross-coupling term is taken into consideration. In Ref.[25], flux linkage components are obtained for the inductances and torque components comparison between two 5-phase 20-slot fault-tolerant PM machines. One has 18-pole spoke-type rotor and the other one has 20-pole V-shape rotor. In Ref.[26], the variations of average line inductance and its ripple with the load of a 24s/4p SPM machine driven by a 6-step 3-phase BLDC inverter are calculated. The FP technique is also employed to calculate the self- and mutual inductances between different phases in a tubular permanent magnet machine with transverse flux configuration for electromagnetic launch (EML) applications in Ref.[27,28], the d- and q-axis inductances of variable-flux flux-intensifying PM machine in Ref.[29], and the d-axis flux linkage components of stator slot PM machine with doubly salient structure in Ref.[30]. In Ref.[31], the on-load PM flux linkage of double-stator flux-switching PM machine using ferrite magnet is calculated by FP technique to explain the torque deterioration due to magnetic saturation under heavy load.

    4 Separation of Torque Components

    Based on the flux linkage components and inductances, it is then possible to discuss the torque component separation. Before the torque component separation, it is necessary to review the torque calculation theories.

    The Maxwell stress tensor and virtual work principle are two most widely used methods for torque calculation. It is well-known that these two methods are identical for torque calculation in normal FE simulations although attention should be paid to the mesh discretization[32,33]. This is also confirmed by Fig.8, which shows the instantaneous torque waveforms of the inset PM machine listed in Tab.1 whena=4.0A and=30°. The torque calculation using virtual work principle is given as

    Fig.8 Torque components of inset PM machine when Ia=4.0A and b =-30°[19]

    Since ?m/?m,d/d, andq/dhave no contribution to the average torque based on the fact that the magnetic field repeats every cycle, the average torque can also be calculated by the classical dq0 model in Equ.(5). Since Equ.(5) does not require derivative operation and the calculation of magnetic energy, it is much simpler and more widely used for torque calculation. However, it should be noted that the classical dq0 model in Equ.(5) is only able to predict the average torque but not the instantaneous torque.

    Combined with FP technique, there are two possible ways of torque component separation. One is based on the Maxwell stress tensor method[34,35]and the other one is based on the virtual work principle[13,33]. However, it is found in Ref.[33] that due to the influence of equivalent rotational magnetic saliency in the stator, a part of PM torque will be improperly attributed to the reluctance torque when the Maxwell stress tensor method is used together with the FP technique. This is eliminated by using the virtual work principle. Therefore, the on-load torque components separation should be based on the combination of virtual work principle and FP technique. Thus, the average torque com- ponents can be separated based on the flux linkage components obtained by the FP technique as

    wherePM,PM(d), andPM(q) are the PM torque and its components due tod(PM) andq(PM) respe- ctively.r,r(d-q), andr(dq) are the reluctance torque and its components due to the self-inductances and mutual inductance, respectively.

    As an example, the torque component variation with the current phase advanced angle of the inset PM machine whena=4.0A is shown in Fig.9.

    The average torque separation method suggested in Ref.[33] is widely accepted and used for various types of PM machines, such as 5-phase 20-slot spoke- type and V-shape fault-tolerant PM machines[25], 24s/4p SPM machine driven by a 6-step 3-phase BLDC inverter[36], 6s/4p asymmetric V-shape, surfaceinset, interior PM (IPM) machines[37-39], 12s/22p switched flux PM machine[40], 12-slot/10-pole SPM machine[41], modular linear flux reversal PM motors[42], and 12s/10p hybrid excited switched flux PM machines[43]. Although the torque component separation is also based on the flux linkage components, the torque components may be defined differently as shown in Ref.[16, 44, 45].

    (a) Main torque components

    (b) Torque components due to cross-coupling

    Fig.9 Torque components of inset PM machine when Ia=4.0A[19]

    Except the above mentioned publications, Ref.[46] suggests another possible way of torque component separation method based on the com- bination of the Maxwell stress tensor and the FP technique. Instead of two linear calculations with either PM or all armature currents, six calculations which are all possible combinations of PM,dandq, are carried out based on the on-load permeability. Then, six torque components are mathematically calculated from the six torque results of six linear analyses based on the Maxwell stress tensor. However, the improper torque attribution due to the influence of equivalent rotational magnetic saliency in the stator remains.

    5 On-load Cogging Torque Analyses

    For PM machines, cogging torque is another widely investigated topic. The cogging torque is produced by the interaction between the rotor PM and the stator magnetic reluctance variation. It is merely a pulsating torque and has no contribution to the average torque. For the open-circuit cogging torque, it can be calculated either by the Maxwell stress tensor and the virtual work principle as

    where0is the permeability of free air space.efis the effective axial length.is the radius of integration path.n(open) andt(open) are the normal and tangential airgap flux density components on open-circuit, respectively.

    For the on-load cogging torque, several calculation methods are suggested by analog with the open-circuit ones in Equ.(12) and Equ.(13). By analog with Equ.(12), the on-load cogging torque in Ref.[22, 34,35] is computed as

    wheremw(FP, PM) is the calculated torque based on the Maxwell stress tensor with on-load PM field only.n(FP, PM) andt(FP, PM) are the normal and tangential airgap flux density components due to the on-load PM field only.

    Being analog with Equ.(13), it is mentioned in Ref.[47] that the on-load cogging torque can be calculated as the derivative of the total magnetic energy with respect to the rotor position at constant current

    However, it is inappropriate since the total magnetic energy is also contributed by the armature field when the machine is on load.

    In Ref.[34], the on-load cogging torque is calculated in a slightly different way as

    wheremlcis the magnetic energy stored in the system less the magnetic energy stored in the coils and includes components as airgap, magnets, rotor and stator steel and any other motor regions, which are not parts of the coils.

    However, without using the FP method, tech- nically, it is impossible to separate the magnetic energy due to the armature field, since the on-load PM and armature fields exist in every region and are coupled each other. Hence, it is also inappropriate.

    In Ref.[48], an improved FP technique is suggested to calculate the on-load cogging torque based on the on-load PM field magnetic energy derivation using Equ.(17). This method is also used later in Ref.[40].

    wherem(FP, PM) is the total magnetic energy with on-load PM field only.

    Fig.10 compares the torque results calculated by Equ.(14) and Equ.(17) of the inset PM machine listedin Tab.1. It shows thatmw(FP,PM) has non-zero average torque and hence cannot be exactly the on-load cogging torque. This is due to the improper torque attribution caused by the equivalent rotating magnetic saliency in the stator[33,48]. By using Equ.(17), the improper torque attribution is inherently eliminated while the influence of electric load and magnetic saturation is still fully included. Hence, the on-load cogging torque can still be obtained properly and shown in Fig.10b.

    (a) By Equ.(14) based on Maxwell stress tensor

    (b) By Equ.(17) based on magnetic energy derivation

    Fig.10 Torque results of inset PM machine with on-load PM field only when Ia=4.0A[19]

    6 On-load Torque Ripple Analyses

    On-load torque ripple is undesirable since it could result in noise and vibration and is harmful to precise rotor position control. In PM machines, the on-load torque ripple is contributed by PM torque ripple, reluctance torque ripple and on-load cogging torque. However, the on-load torque ripple is often investigated as single resultant quantity. With the help of FP technique, it is possible to identify the contribution and variation of each on-load torque ripple component. This would greatly help the machine designs.

    In Ref.[22], the on-load torque ripple of a 12s/10p IPM machine with rotor shaping is further investigated based upon the calculation of the on-load cogging torque, the on-load back EMF and the inductance variation. In Ref.[41], all on-load torque ripple components are calculated using FP technique. The design trade-off between the open-circuit cogging torque and the on-load torque ripple is then investigated in a 12s/8p SPM machine having different slot openings. In Ref.[49], an analytical model is developed to predict the on-load transient torque of IPM machines by including not only the fundamental but also the harmonics of flux linkages. The flux linkage harmonics are segregated by FP technique to further explain the torque ripple variation. By way of example, on-load torque components as well as the variation of torque ripple components with the slot opening in a 12s/8p SPM machine are presented in Fig.11[41].

    (a) On-load torque components when slot opening is 1mm

    (b) Variation of torque ripple components with slot opening

    Fig.11 On-load torque components and variations of torque ripple components with slot opening in a 12s/8p SPM machine when Ia=10A and ??=0°[41]

    7 On-load Back EMF Analyses

    Back EMF, which is produced by PM, is one of the most important quantities regarding the terminal voltage and output torque, both their average value and ripples. Normally, the back EMFs are investigated under the open-circuit condition. Due to the influence of magnetic saturation, the on-load back EMF can be significantly different from the open-circuit EMF. With the help of the FP technique, the on-load back EMF can be obtained based on the on-load PM flux linkage variation with the rotor position.

    The relevant investigations on the on-load back EMFs can be found in Ref.[21,22] on a 12s/10p IPM machine with rotor shaping, Ref.[35] on a 12s/4p SPM machine, Ref.[50] on 24s/4p, 36s/6p, and 48s/8p multi-layer IPM machines as well as Ref.[51] on a 24s/4p SPM machine driven by a 6-step 3-phase BLDC inverter. All these investigations show that the on-load back EMFs are influenced by the electrical loading and can be significantly different from the open-circuit back EMFs. By way of example, the open-circuit and on-load back EMFs waveforms in Ref.[21] are presented in Fig.12.

    Fig.12 Open-circuit and on-load phase back EMF waveforms of 12s/10p IPM machine in Ref.[21]

    8 On-load Voltage Waveform Analyses

    When the load changes, not only the on-load back EMF but also the on-load voltage are affected. More importantly, the on-load voltage waveform may be greatly distorted from sinusoidal even when the open-circuit back EMFs are designed to be sinusoidal, as will be shown later. This on-load voltage distortion is harmful to various performance, especially the torque- speed curve. In order to investigate the mechanism of the on-load voltage distortion in various PM machines, the FP technique is employed in Ref.[52-57].

    The investigation in Ref.[53] is focused on a 12s/8p SPM machine with its parameters and cross-section given in Tab.2 and Fig.13. The influence of slot/pole number combination in SPM machine is carried out in Ref.[54]. The investigation in Ref.[55] is focused on the fractional slot IPM machine while the influence of slot/pole number combination in IPM machine is carried out in Ref.[56,57] is focused on the comparison of torque-speed curves between 42s/8p and 48s/8p IPM machines considering the on-load voltage distortions.

    Fig.14 and Fig.15 show the waveforms of the voltage and flux linkage components of the 12s/8p SPM machine in Ref.[52] as an example. The machine parameters and cross-section are given in Tab.2 and Fig.13. It can be seen from Fig.13 that although the open-circuit back EMF is sinusoidal, the on-load voltage waveform is greatly distorted from sinusoidal and exhibits sharp voltage spikes. Based on the FP technique, it can be seen from Fig.14 that the sharp voltage spikes are results of sudden changes of both the PM and armature flux linkages. The sudden flux linkage changes are further due to the variation local magnetic saturation in the tooth-tips as shown in Fig.16.

    Tab.2 Main parameters of 12s/8p SPM machine[52]

    Fig.13 Cross-section of 12s/8p SPM machine in Ref.[52]

    Fig.14 Open-circuit phase back EMF and on-load (Ia=10A and b =0°) terminal voltage waveforms of 12s/8p SPM machine[52]

    (a) Flux linkages

    (b) Voltages

    Fig.15 Flux linkage and voltage components of 12s/8p SPM machine when Ia=10A and ??=0°[52]

    9 Flux Weakening Performance Analyses

    Flux weakening control is an essential techniqueto achieve wide constant power speed range which is very often required for high performance electrical machines, such as electric vehicle applications. However, the constant parameter machine models could result in significant inaccuracy on the flux weakening performance prediction[15,58]. By using FP technique, the on-load flux linkages and inductances can be calculated. Then, it is possible to predict the flux weakening performance more accurately. The examples for flux weakening performance prediction based on FP technique calculated parameters can be found in Ref.[13,15,22,25,52-57,59]. However, it should be noted that the flux weakening performance prediction only based on the fundamental component of phase flux linkages (or average d- and q-axis flux linkages) could also result in large inaccuracy when the on-load voltage is distorted from sinusoidal[52-57]. It is due to that the voltage distortion will make the peak voltage much larger than the fundamental value. Consequently, the machine enters the flux weakening operation much earlier when the DC bus voltageDCis fixed.

    (a) Locations of detecting points

    (b) Variation of permeability in tooth-tips

    Fig.16 Relative permeability for tooth-tip bridges between each phase coils of 12s/8p SPM machine when Ia=10A and ??=0°[52]

    The influence of on-load voltage distortion on the flux weakening performance of the 12s/8p SPM machine listed in Tab.2, which is presented in Ref.[52], is shown in Fig.17 as an example. It can be seen that the actual base speed considering the voltage distortion is much lower than the one considering the fundamental components only. The maximum torque in the deep flux weakening area also drops much faster when considering the voltage distortion.

    (a) Torque and torque ripple against speed

    (b) d,q-axis currents against speed

    Fig.17 Variation of torque-speed characteristics for different calculation methods of 12s/8p SPM machine when Ia=10A and UDC=42V[52]

    10 Analyses on Other Issues

    Except the aforementioned parameters and per- formance, the FP technique has also been used to help investigating other issues in PM machines. For example, in Ref.[60,61], the radial fore distributions are analyzed based on the on-load airgap flux density components. In Ref.[62], the on-load PM leakage flux through the iron bridge of an IPM machine with segmented PMs is obtained to explain the enhance- ment of its flux weakening capability.

    11 Applications to Other Types of Electri- cal Machines

    With the foregoing sections focused on the applications in PM machines, this section is dedicated to the applications of the FP technique to the other types of electrical machines.

    The relevant applications to induction machines can be found in Ref.[63-65]. In Ref.[63], the FP technique is used to identify the influence of local magnetic saturation near the broken rotor bar on the stator current. In Ref.[64,65], the FP method is employed to calculate the equivalent circuit para- meters of induction machines, especially the leakage inductances.

    For electrically excited machines, the relevant publications can be found in Ref.[66-68]. In Ref.[66], the on-load magnetic field and flux linkage com- ponents are obtained to explain the terminal voltage increase of a PM-assisted salient pole electrically excited machine. In Ref.[67,68], variation of inductances with the load is calculated to improve the maximum torque per ampere control of electrically excited machine under heavy magnetic saturation operations.

    The investigations in Ref.[17] and Ref.[69] are on SR machines. In Ref.[17], the on-load phase flux linkage variations are obtained to explain the influ- ence of mutual coupling effect on torque production in SR machines having different winding topologies. In Ref.[69], the self- and mutual torques of conventional and mutually coupled SR machines are obtained based on the on-load self- and mutual flux linkages for performance comparison and selection of optimal current waveform.

    In Ref.[43], the FP technique is used on a hybrid excited switched flux PM machine to obtain flux linkages and inductances for further performance investigation including torque components and torque- speed curves. In Ref.[70], the torque separation is conducted by FP technique to assist the optimization of a hybrid rare-earth-free PM reluctance machine.

    12 Conclusion

    The FP technique and its applications in various electrical machines have been comprehensively reviewed in this paper. It shows that the FP technique is able to accurately separate on-load magnetic field, flux linkage and inductance components in all types of electrical machines. The on-load field components and parameters can be further utilized for com- prehensive on-load performance analyses, such as average torque and torque ripple, voltage, flux weakening and radial force distribution. The accurate parameters can also improve the machine models and controls. All these features are very useful on developing high performance electrical machines.

    Acknowledgments: The authors also thank Dr. Di Wu and Dr. Ziad Azar for their technical contributions.

    Reference

    [1] Chan C C. The state of the art of electric, hybrid, and fuel cell vehicles[J]. Proceedings of the IEEE, 2007, 95(4): 704-718.

    [2] Zhu Z Q, Howe D. Electrical machines and drives for electric, hybrid, and fuel cell vehicles[J]. Proceedings of the IEEE, 2007, 95(4): 746-765.

    [3] El-Refaie A M. Motors/generators for traction/ propulsion applications: a review[J]. IEEE Vehicular Technology Magazine, 2013, 8(1): 90-99.

    [4] Cao W P, Mecrow B C, Atkinson G J, et al. Overview of electric motor technologies used for more electric aircraft (MEA)[J]. IEEE Transactions on Industry Electronics, 2012, 59(9): 3523-3531.

    [5] Cheng M, Hua W, Zhang J, et al. Overview of stator- permanent magnet brushless machines[J]. IEEE Transactions on Industry Electronics, 2011, 58(11): 5087-5101.

    [6] Zhu Z Q. Switched flux permanent magnet machines- innovation continues[C]//International Conference on Electrical Machines and Systems, 2011: 1-10.

    [7] Eastham J F, Cox T, Proverbs J. Application of planar modular windings to linear induction motors by harmonic cancellation[J]. IET Electric Power App- lications, 2010, 4(3): 140-148.

    [8] Levi E, Bojoi R, Profumo F, et al. Multiphase induction motor drives—a technology status review[J]. IET Electric Power Applications, 2007, 1(4): 489- 516.

    [9] Bash M L, Pekarek S. Analysis and validation of a population-based design of a wound-rotor syn- chronous machine[J]. IEEE Transactions on Energy Conversion, 2012, 27(3): 603-614.

    [10] Vijayakumar K, Karthikeyan R, Paramasivam S, et al. Switched reluctance motor modeling, design, simulation, and analysis: a comprehensive review[J]. IEEE Transactions on Magnetics, 2008, 44(12): 4605-4617.

    [11] Bianchi N, Bolognani S, Bon D, et al. Torque harmonic compensation in a synchronous reluctance motor[J]. IEEE Transactions on Energy Conversion, 2008, 23(2): 466-473.

    [12] Amara Y, Vido L, Gabsi M, et al. Hybrid excitation synchronous machines: energy-efficient solution for vehicles propulsion[J]. IEEE Transactions on Vehicular Technology, 2009, 58(5): 2137-2149.

    [13] Bianchi N, Bolognani S. Magnetic models of saturated interior permanent magnet motors based on finite element analysis[C]//IEEE Industrial App- lication Conference, St. Louis, MO, USA, 1998: 27- 34.

    [14] Williamson S, Knight A M. Performance of skewed single-phase line-start permanent magnet motors[J]. IEEE Transactions on Industry Applications, 1999, 35(3): 577-582.

    [15] Kang G H, Hong J P, Kim G T, et al. Improved parameter modeling of interior permanent magnet synchronous motor based on finite element analysis[J]. IEEE Transactions on Magnetics, 2000, 36(4): 1867- 1870.

    [16] Hu J, Zou J, Liang W. Finite element calculation of the saturation dq-axes inductance for a direct drive PM synchronous motor considering cross-magneti- zation[C]//International Conference on Power Elec- tronics and Drive Systems, 2003: 677-681.

    [17] Walker J A, Dorrell D G, Cossar C. Effect of mutual coupling on torque production in switched reluctance motors[J]. Journal of Applied Physics, 2006, 99(10): 08R304(1-3).

    [18] Walker J A, Dorrell D G, Cossar C. Flux-linkage calculation in permanent-magnet motors using the frozen permeabilities method[J]. IEEE Transactions on Magnetics, 2005, 41(10): 3946-3948.

    [19] Chu W Q. Investigation of torque characteristics of permanent magnet and electrically excited machines[D]. UK: The University of Sheffield, 2013.

    [20] Kwak S Y, Kim J K, Jung H K. Characteristic analysis of multilayer-buried magnet synchronous motor using fixed permeability method[J]. IEEE Transactions on Energy Conversion, 2005, 20(3): 549-555.

    [21] Azar Z. Electromagnetic performance of fractional slot interior permanent magnet and synchronous reluctance brushless ac machines having non- overlapping concentrated windings[D]. UK: The University of Sheffield, 2012.

    [22] Azar Z, Zhu Z Q, Ombach G. Influence of electric loading and magnetic saturation on cogging torque, back-EMF and torque ripple of PM machines[J]. IEEE Transactions on Magnetics, 2012, 48(10): 2650-2658.

    [23] Alam F R, Abbaszadeh K. Magnetic field analysis in eccentric surface-mounted permanent-magnet motors using an improved conformal mapping method[J]. IEEE Transactions on Energy Conversion, 2016, 31(1): 333-344.

    [24] Schmidt E, Susic M. Parameter evaluation of permanent magnet synchronous machines with tooth coil windings using the frozen permeabilities method with the finite element analyses[C]//IEEE Canadian Conference on Electrical & Computer Engineering, Montreal, QC, 2012: 1-5.

    [25] Chen Q, Liu G, Zhao W, et al. Design and comparison of two fault-tolerant interior-permanent-magnet motors[J]. IEEE Transactions on Industrial Electronics, 2014, 61(12): 6615-6623.

    [26] Paula G T D, Monteiro J R B d A, Almeida T E P d, et al. Evaluation of surface mounted PM machine's parameters on load conditions using frozen permeabi- lity method. Part. I[C]//IEEE International Sym- posium on Industrial Electronics, Berlin, 2014: 806-811.

    [27] Wang Q, Zhao B, Xu Y, et al. Inductances and phase coupling analysis of tubular permanent magnet machines with transverse flux configuration[C]// International Symposium on Electromagnetic Launch Technology, 2014: 1-5.

    [28] Wang Q, Zhao B, Zhang J, et al. Inductances and phase coupling analysis of tubular permanent magnet machines with transverse flux configuration[J]. IEEE Transactions on Plasma Science, 2015, 43(5): 1232- 1235.

    [29] Sun A, Li J, Qu R. Inductance calculation in variable-flux flux-intensifying permanent magnet synchronous machines using improved frozen per- meability method[C]//IEEE Magnetics Conference, Beijing, 2015: 1-1.

    [30] Afinowi I, Zhu Z Q, Guan Y, et al. Electromagnetic performance of stator slot permanent magnet machines with/without stator tooth-tips and having single/double layer windings[J]. IEEE Transactions on Magnetics, 2016, 52(6): 8103410.

    [31] Kim D, Hwang H, Bae S, et al. Analysis and design of a double-stator flux-switching permanent magnet machine using ferrite magnet in hybrid electric vehicles[J]. IEEE Transactions on Magnetics, 2016, 52(7): 8106604.

    [32] Bianchi N, Alberti L. MMF harmonics effect on the embedded FE analytical computation of PM motors[J]. IEEE Transactions on Industry Applications, 2010, 46(2): 812-820.

    [33] Chu W Q, Zhu Z Q. Average torque separation in permanent magnet synchronous machines using frozen permeability[J]. IEEE Transactions on Mag- netics, 2013, 49(3): 1202-1210.

    [34] Ionel D M, Popescu M, McGilp M I, et al. Assessment of torque components in brushless permanent-magnet machines through numerical analysis of the elec- tromagnetic field[J]. IEEE Transactions on Industry Applications, 2005, 41(5): 1149-1158.

    [35] Abbaszadeh K, Alam F R. On-load field component separation in surface-mounted permanent- magnet motors using an improved conformal mapping method[J]. IEEE Transactions on Magnetics, 2016, 52(2): 5200112(1-12).

    [36] Paula G T D, Monteiro J R B d A, Almeida T E P d, et al. Evaluation of surface mounted PM machine's parameters on load conditions using frozen per- meability method. Part. II[C]//IEEE/IAS International Conference on Industry Applications, Berlin, 2014: 1-7.

    [37] Zhao W, Zhao F, Lipo T A, et al. Optimal design of a novel V-type interior permanent magnet motor with assisted barriers for the improvement of torque characteristics[J]. IEEE Transactions on Magnetics, 2014, 50(11): 8104504(1-4).

    [38] Zhao W, Chen D, Lipo T A, et al. Performance improvement of ferrite-assisted synchronous relu- ctance machines using asymmetrical rotor con- figurations[J]. IEEE Transactions on Magnetics, 2015, 51(11): 8108504(1-4).

    [39] Zhao W, Lipo T A, Kwon B I. Optimal design of a novel asymmetrical rotor structure to obtain torque and efficiency improvement in surface inset PM motors[J]. IEEE Transactions on Magnetics, 2015, 51(3): 8100704(1-4).

    [40] Fei W, Luk P C K, Miao D M, et al. Investigation of torque characteristics in a novel permanent magnet flux switching machine with an outer-rotor con- figuration[J]. IEEE Transactions on Magnetics, 2014, 50(4): 8100810(1-10).

    [41] Wu D, Zhu Z Q. Design tradeoff between cogging torque and torque ripple in fractional slot surface- mounted permanent magnet machines[J]. IEEE Transa- ctions on Magnetics, 2015, 51(11): 8108704(1-4).

    [42] Xu L, Liu G, Zhao W, et al. Analysis of new modular linear flux reversal permanent magnet motors[J]. IEEE Transactions on Magnetics, 2015, 51(11): 8109904(1-4).

    [43] Li G J, Zhu Z Q, Jewell G. Performance investigation of hybrid excited switched flux permanent magnet machines using frozen permeability method[J]. IET Electric Power Applications, 2015, 9(9): 586-594.

    [44] Tangudu J K, Jahns T M, El-Refaie A M, et al. Segregation of torque components in fractional-slot concentrated-winding interior PM machines using frozen permeability[C]//IEEE Energy Conversion Congress and Exposition, San Jose, CA, 2009: 3814- 3821.

    [45] Xia B, Fei W, Luk P. Analysis and design of V-spoke ferrite interior permanent magnet machine for traction applications[C]//International Conference on Power Electronics Systems and Applications, Hong Kong, 2015: 1-6.

    [46] Chen X, Wang J, Patel V I, et al. Reluctance torque evaluation for interior permanent magnet machines using frozen permeability[C]//IET International Conference on Power Electronics, Machines and Drives, Manchester, 2014: 1-6.

    [47] Popescu M, Ionel D M, Miller T J E, et al. Improved finite element computations of torque in brushless permanent magnet motors[J]. IEE Proceedings on Electrical Power Applications, 2005, 152(2): 271-276.

    [48] Chu W Q, Zhu Z Q. On-load cogging torque calculation in permanent magnet machines[J]. IEEE Transactions on Magnetics, 2013, 49(6): 2982-2989.

    [49] Pina A J, Pramod P, Islam R, et al. Extended model of interior permanent magnet synchronous motors to include harmonics in d- and q-axes flux linkages[C]// IEEE Energy Conversion Congress and Exposition, Montreal, QC, 2015: 1864-1871.

    [50] Seo J H, Kwak S Y, Jung S Y, et al. Investigation on EMF waveform in the interior permanent magnet synchronous machine considering load condition[C]// Biennial IEEE Conference on Electromagnetic Field Computation, Miami, FL, 2006: 321-321.

    [51] Paula G T D, Monteiro J R B d A, Almeida T E P d, et al. Evaluation of surface mounted PM machine's parameters on load conditions using frozen per- meability method. Part. II[C]//International Con- ference on Electrical Machines, Berlin, 2014: 156- 161.

    [52] Wu D. Effect of magnetic saturation in fractional slot PM machines with particular reference to terminal voltage distortion[D]. UK: The University of Sheffield, 2015.

    [53] Wu D, Zhu Z Q. On-load voltage distortion in fractional slot surface-mounted permanent magnet machines considering local magnetic saturation[J]. IEEE Transactions on Magnetics, 2015, 51(8): 8106410(1-10).

    [54] Wu D, Zhu Z Q. Influence of slot and pole number combinations on voltage distortion in surface- mounted permanent magnet machines with local magnetic saturation[J]. IEEE Transactions on Energy Conversion, 2015, 30(4): 1460-1471.

    [55] Zhu Z Q, Wu D. On-load voltage distortion in fractional-slot interior permanent magnet machines[J]. IEEE Transactions on Magnetics, 2015, 51(10): 1-9.

    [56] Zhu Z Q, Wu D, Ge X. Investigation of voltage distortion in fractional slot interior permanent magnet machines having different slot and pole number combinations[J]. IEEE Transactions on Energy Conversion, 2016, 31(3): 1192-1201.

    [57] Zhu Z Q, Wu D, Wu M C, et al. Influence of on-load voltage distortion on torque-speed characteristic of interior permanent magnet machines[C]//IEEE Energy Conversion Congress and Exposition, 2015: 760- 767.

    [58] Qi G, Chen J T, Zhu Z Q, et al. Influence of skew and cross-coupling on flux-weakening performance of permanent-magnet brushless AC machines[J]. IEEE Transactions on Magnetics, 2009, 45(5): 2110-2117.

    [59] Chai S H, Lee B H, Hong J P, et al. Design of IPMSM having high power density for position sensorless operation with high-frequency signal injection and the method of calculating inductance profile[C]// International Conference on Electrical Machines and Systems, Beijing, 2011: 1-5.

    [60] Dajaku G, Gerling D. Magnetic radial force density of the PM machine with 12-teeth/10-poles winding topology[C]//IEEE International Electric Machines and Drives Conference, Miami, FL, 2009: 1715-1720.

    [61] Xia X P, Zhu Z Q, Wu L J, et al. Comparison of radial vibration forces in 10-pole/12-slot fractional slot surface-mounted and interior PM brushless AC machines[C]//International Conference on Electrical Machines, Rome, 2010: 1-6.

    [62] Duan S, Zhou L, Wang J. Flux weakening mechanism of interior permanent magnet synchronous machines with segmented permanent magnets[J]. IEEE Transa- ctions on Apply Superconducting, 2014, 24(3): 1-5.

    [63] Sprooten J, Gyselinck J, Maun J C. Local and global effect of a broken bar in induction machines using fundamental electromagnetic laws and finite element simulations[C]//IEEE International Symposium on Diagnostics for Electric Machines, Power Electronics and Drives, Vienna , 2005: 1-6.

    [64] Ling Z, Zhou L, Guo S, et al. Equivalent circuit parameters calculation of induction motor by finite element analysis[J]. IEEE Transactions on Magnetics, 2014, 50(2): 833-836.

    [65] Jia L. Equivalent circuit parameters calculation of a wound rotor brushless doubly-fed machine based on finite element analysis[C]//IEEE Magnetics Con- ference, Beijing , 2015: 1.

    [66] Shirai K, Tokikuni Y, Shima K, et al. Causes of increase in the terminal voltage of a permanent- magnet-assisted salient-pole synchronous machine[C]// International Conference on Electrical Machines and Systems, Tokyo, 2009: 1-6.

    [67] Jeong I, Kim J, Kim Y, et al. Extended MTPA with cross coupling inductances for electrically excited synchronous motors[C]//IEEE Energy Conversion Congress and Exposition, Denver, CO, 2013: 867-873.

    [68] Jeong I, Gu B G, Kim J, et al. Inductance estimation of electrically excited synchronous motor via polynomial approximations by least square method[J]. IEEE Transactions on Industry Applications, 2015, 51(2): 1526-1537.

    [69] Li G J, Zhu Z Q, Ma X Y, et al. Comparative study of torque production in conventional and mutually coupled SRMs using frozen permeability[J]. IEEE Transactions on Magnetics, 2016, 52(6): 8103509.

    [70] Hwang H, Bae S, Lee C. Analysis and design of a hybrid rare-earth-free permanent magnet reluctance machine by frozen permeability method[J]. IEEE Transa- ctions on Magnetics, 2016, 52(7): 8106304(1-4).

    凍結磁導率先進技術及其在高性能電機研發(fā)中的應用

    諸自強 褚文強

    (謝菲爾德大學 謝菲爾德 S1 3JD 英國)

    如何考慮磁場飽和的影響是當前電機設計中一個亟待解決的重要問題。近年來提出的凍結磁導率技術越來越多地被應用到電機負載電磁性能和參數(shù)計算中。本文詳細闡述了凍結磁導率的原理及其在各種電機中的應用現(xiàn)狀。研究表明,凍結磁導率技術可以用來精確分離各種電機負載狀態(tài)下的電磁分量,如磁場、磁鏈、電感、轉矩、轉矩波動、反電勢和端電壓、弱磁性能以及徑向力等均可以用凍結磁導率技術預測,并用來輔助電機及其驅動控制系統(tǒng)的設計。因此,凍結磁導率技術為高性能電機的研發(fā)提供了一個全新的方法。

    電機 凍結磁導率 負載性能 永磁

    TM301.4

    June 20, 2016; Revised July 1, 2016.

    Note brief

    Professor Zhu Z Q, born in 1962, Fellow IEEE, Fellow IET, PhD, Professor at the University of Sheffield, UK. Major research interests include design, control, and applications of brushless permanent magnet machines and drives for applications ranging from automotive to renewable energy.

    E-mail: Z.Q.Zhu@Sheffield.ac.uk (Corresponding author)

    Dr Chu W Q, born in 1982, Principal Engineer, CSR Electric Drives Technology Research Centre, the University of Sheffield, UK. Major research interests include permanent magnet electrical machines and applications.

    E-mail: chu.wenqiang@gmail.com

    This work is partially supported by Guangdong Welling Motor Manufacturing Co. Ltd, Guangdong Innovative Research Team Program (2011N084), China, and the Royal Academy of Engineering/ Siemens Research Chair Program, UK.

    猜你喜歡
    謝菲爾德磁導率端電壓
    擊沉“海上堡壘”
    寬頻高磁導率R10k軟磁材料的開發(fā)
    山東冶金(2022年2期)2022-08-08 01:50:52
    基于FEMM的永磁電機動態(tài)凍結磁導率并行仿真及程序
    防爆電機(2021年3期)2021-07-21 08:12:50
    勵磁調節(jié)器PT回路故障后機端電壓變化分析
    人工心臟無位置傳感器無刷直流電動機非導通相端電壓分析
    微特電機(2016年3期)2016-11-29 09:43:45
    鈦管蒸發(fā)器在酸性溶液中的防腐防垢
    鋼板磁導率變化對船舶感應磁場的影響
    通信機房不間斷電源蓄電池的應用與維護
    江蘇通信(2014年4期)2014-04-14 06:18:42
    英國謝菲爾德大學實訓項目開關電源的設計與實施
    崇尚職業(yè)教育的大學典范——英國謝菲爾德哈勒姆大學
    亚洲精品日韩在线中文字幕| 人妻系列 视频| av线在线观看网站| 精品不卡国产一区二区三区| 久久这里只有精品中国| 丝袜美腿在线中文| 特大巨黑吊av在线直播| 欧美激情久久久久久爽电影| 91av网一区二区| 国国产精品蜜臀av免费| 欧美色视频一区免费| 久久久a久久爽久久v久久| 久久久久久久久久久免费av| 成人国产麻豆网| 乱人视频在线观看| 欧美3d第一页| 成人高潮视频无遮挡免费网站| av在线观看视频网站免费| 国产亚洲午夜精品一区二区久久 | 亚洲伊人久久精品综合 | 久久久欧美国产精品| av在线蜜桃| 爱豆传媒免费全集在线观看| 亚洲婷婷狠狠爱综合网| 久久人妻av系列| 亚洲成人av在线免费| 久久久久精品久久久久真实原创| 国产精品久久久久久av不卡| 又爽又黄a免费视频| 三级经典国产精品| 国产精品一及| 99在线视频只有这里精品首页| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 爱豆传媒免费全集在线观看| 午夜福利视频1000在线观看| 十八禁国产超污无遮挡网站| av在线蜜桃| 国产成人精品婷婷| 国产麻豆成人av免费视频| 久久精品熟女亚洲av麻豆精品 | 国产精品福利在线免费观看| 国产免费视频播放在线视频 | 黄色一级大片看看| 欧美日韩精品成人综合77777| 日本免费a在线| 久久99热6这里只有精品| 国产成人午夜福利电影在线观看| 少妇人妻一区二区三区视频| 亚洲av熟女| 老司机福利观看| 伊人久久精品亚洲午夜| 亚洲综合色惰| 国产真实乱freesex| 啦啦啦观看免费观看视频高清| 亚洲欧美精品自产自拍| 综合色av麻豆| 国产v大片淫在线免费观看| 国产探花极品一区二区| 五月伊人婷婷丁香| 亚洲五月天丁香| 久久精品国产亚洲av天美| 99久久九九国产精品国产免费| 国产一区二区在线观看日韩| 一本久久精品| 99久久精品一区二区三区| ponron亚洲| 色哟哟·www| 能在线免费观看的黄片| 熟妇人妻久久中文字幕3abv| 日韩精品有码人妻一区| a级毛片免费高清观看在线播放| 一级二级三级毛片免费看| 亚洲自拍偷在线| 国产精品三级大全| 18禁裸乳无遮挡免费网站照片| 国产麻豆成人av免费视频| 欧美97在线视频| 日本黄色视频三级网站网址| 精品久久久久久久久亚洲| 国产精品爽爽va在线观看网站| 男女边吃奶边做爰视频| 卡戴珊不雅视频在线播放| 大又大粗又爽又黄少妇毛片口| 亚洲成av人片在线播放无| 看十八女毛片水多多多| 亚洲精品乱码久久久久久按摩| 亚洲av.av天堂| 久久综合国产亚洲精品| 日本一二三区视频观看| 欧美成人精品欧美一级黄| 可以在线观看毛片的网站| 亚洲成色77777| 自拍偷自拍亚洲精品老妇| 精品午夜福利在线看| av视频在线观看入口| 最近2019中文字幕mv第一页| 在线播放国产精品三级| 国产成人一区二区在线| 国产免费一级a男人的天堂| 黄片wwwwww| 久久综合国产亚洲精品| 国产成人aa在线观看| 直男gayav资源| 久久久久久久久久成人| 男的添女的下面高潮视频| 日韩制服骚丝袜av| 综合色丁香网| 国产人妻一区二区三区在| 99久久人妻综合| 欧美+日韩+精品| 亚洲av日韩在线播放| 亚洲欧美成人综合另类久久久 | 麻豆av噜噜一区二区三区| 欧美97在线视频| 一级黄色大片毛片| 亚洲欧美日韩无卡精品| 麻豆国产97在线/欧美| 能在线免费看毛片的网站| 久久草成人影院| 亚洲人成网站在线观看播放| 日本-黄色视频高清免费观看| 久久精品国产自在天天线| 听说在线观看完整版免费高清| 午夜精品在线福利| 色尼玛亚洲综合影院| 国产免费男女视频| 日韩人妻高清精品专区| 日韩三级伦理在线观看| 午夜免费激情av| 中文字幕免费在线视频6| videos熟女内射| 搞女人的毛片| 久久人妻av系列| 日韩成人av中文字幕在线观看| 久久欧美精品欧美久久欧美| 亚洲欧美精品自产自拍| 国产精品野战在线观看| 国产免费又黄又爽又色| 99久久精品一区二区三区| 简卡轻食公司| 最近2019中文字幕mv第一页| 国产视频内射| 免费人成在线观看视频色| 亚洲在久久综合| av在线蜜桃| 国产三级在线视频| 99久久精品热视频| 国产一区二区在线av高清观看| 视频中文字幕在线观看| 国产精品久久视频播放| 国内精品宾馆在线| 国产69精品久久久久777片| 久久婷婷人人爽人人干人人爱| 九九久久精品国产亚洲av麻豆| 亚洲国产欧美在线一区| 亚洲18禁久久av| 嫩草影院精品99| 日日干狠狠操夜夜爽| 97热精品久久久久久| 精品人妻熟女av久视频| 极品教师在线视频| 中文字幕人妻熟人妻熟丝袜美| 波多野结衣巨乳人妻| 伊人久久精品亚洲午夜| 两个人的视频大全免费| 亚洲精品日韩av片在线观看| 99热网站在线观看| 九九热线精品视视频播放| 亚洲av熟女| 久久久久网色| 乱人视频在线观看| 22中文网久久字幕| 亚洲精品乱码久久久久久按摩| 亚洲国产精品sss在线观看| 久久久久久久久久黄片| 免费搜索国产男女视频| av女优亚洲男人天堂| 青春草国产在线视频| 国产91av在线免费观看| 天堂影院成人在线观看| 免费av观看视频| 一级黄色大片毛片| 亚洲精品自拍成人| 精华霜和精华液先用哪个| 国产欧美另类精品又又久久亚洲欧美| 色综合站精品国产| av在线播放精品| 一本久久精品| 国产单亲对白刺激| 日本一本二区三区精品| 亚洲精华国产精华液的使用体验| 国产精品国产高清国产av| 久久久久久九九精品二区国产| 91av网一区二区| 欧美精品国产亚洲| 欧美一区二区亚洲| 一区二区三区乱码不卡18| 中文字幕精品亚洲无线码一区| 亚洲人成网站在线播| 精品久久久久久久人妻蜜臀av| 小说图片视频综合网站| 蜜臀久久99精品久久宅男| 亚洲美女视频黄频| 我要看日韩黄色一级片| 午夜福利在线观看吧| 桃色一区二区三区在线观看| 成人美女网站在线观看视频| 国产精品国产高清国产av| 午夜福利在线在线| 美女黄网站色视频| 一个人免费在线观看电影| 精品一区二区三区视频在线| 国产精品久久电影中文字幕| 97超碰精品成人国产| 国产午夜福利久久久久久| 中文欧美无线码| 一边亲一边摸免费视频| 一区二区三区四区激情视频| 国产精品一区www在线观看| 亚洲国产成人一精品久久久| 啦啦啦韩国在线观看视频| 精品少妇黑人巨大在线播放 | 亚洲乱码一区二区免费版| 我的女老师完整版在线观看| 国产视频首页在线观看| 精品久久久久久久末码| 免费看av在线观看网站| 乱码一卡2卡4卡精品| 亚洲欧美成人精品一区二区| 国产黄片美女视频| 久久精品人妻少妇| 久久久成人免费电影| 亚洲国产欧洲综合997久久,| 在线观看美女被高潮喷水网站| 国产片特级美女逼逼视频| 精品99又大又爽又粗少妇毛片| 中文字幕久久专区| 黄色配什么色好看| 欧美一区二区精品小视频在线| 久久精品国产自在天天线| 国产黄a三级三级三级人| 精品免费久久久久久久清纯| 亚洲四区av| 成人二区视频| 美女被艹到高潮喷水动态| www.av在线官网国产| 久久热精品热| 少妇被粗大猛烈的视频| 欧美性猛交黑人性爽| 一边摸一边抽搐一进一小说| 一级黄片播放器| 十八禁国产超污无遮挡网站| 久久久久久久久久成人| 免费无遮挡裸体视频| 小说图片视频综合网站| 国产一区二区三区av在线| 国产欧美另类精品又又久久亚洲欧美| 日本av手机在线免费观看| 国产真实乱freesex| 亚洲av不卡在线观看| 青春草视频在线免费观看| 伊人久久精品亚洲午夜| 久久久亚洲精品成人影院| 久久这里只有精品中国| 嫩草影院精品99| 99久久精品国产国产毛片| 蜜桃亚洲精品一区二区三区| 国产淫语在线视频| av天堂中文字幕网| 精品国产三级普通话版| 听说在线观看完整版免费高清| 在现免费观看毛片| 成人亚洲精品av一区二区| 赤兔流量卡办理| 日本免费在线观看一区| 婷婷色综合大香蕉| 久久久久久国产a免费观看| 美女内射精品一级片tv| 久久精品国产亚洲网站| 级片在线观看| 不卡视频在线观看欧美| 亚洲av中文av极速乱| 色吧在线观看| 日本av手机在线免费观看| 简卡轻食公司| 国产av一区在线观看免费| 亚洲,欧美,日韩| 人妻夜夜爽99麻豆av| 亚洲av男天堂| 最近中文字幕高清免费大全6| 水蜜桃什么品种好| 亚洲欧美成人精品一区二区| 精品人妻偷拍中文字幕| 久久精品久久久久久噜噜老黄 | 国产 一区 欧美 日韩| 国产黄片美女视频| 哪个播放器可以免费观看大片| 久久综合国产亚洲精品| 久久久久性生活片| 欧美变态另类bdsm刘玥| 麻豆久久精品国产亚洲av| 日日撸夜夜添| 性插视频无遮挡在线免费观看| 麻豆av噜噜一区二区三区| 在线播放无遮挡| 久久精品人妻少妇| 国产亚洲91精品色在线| 国产片特级美女逼逼视频| 久久久a久久爽久久v久久| 看十八女毛片水多多多| 成人三级黄色视频| 我要搜黄色片| 干丝袜人妻中文字幕| 日韩成人伦理影院| 日韩av在线免费看完整版不卡| 亚洲国产欧洲综合997久久,| 久久久久九九精品影院| av线在线观看网站| 大香蕉久久网| eeuss影院久久| 一区二区三区四区激情视频| 国产精品久久久久久久久免| 中国美白少妇内射xxxbb| 午夜福利在线在线| 欧美性猛交黑人性爽| 欧美xxxx性猛交bbbb| 成人欧美大片| 欧美一区二区精品小视频在线| 一边亲一边摸免费视频| 国内精品宾馆在线| 精品国内亚洲2022精品成人| av.在线天堂| 亚洲精品乱久久久久久| 好男人在线观看高清免费视频| 丝袜喷水一区| 高清午夜精品一区二区三区| 国产精品,欧美在线| 一区二区三区免费毛片| 日韩欧美国产在线观看| 国语自产精品视频在线第100页| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 成人鲁丝片一二三区免费| 久久久久久久久久成人| 免费观看a级毛片全部| 两个人视频免费观看高清| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 婷婷六月久久综合丁香| 亚洲五月天丁香| 国产色爽女视频免费观看| 欧美日韩一区二区视频在线观看视频在线 | 成人毛片a级毛片在线播放| 日韩强制内射视频| 日韩一区二区三区影片| 亚洲成av人片在线播放无| 久久久午夜欧美精品| 综合色av麻豆| 久久精品国产亚洲网站| 欧美激情国产日韩精品一区| 免费观看a级毛片全部| 六月丁香七月| 真实男女啪啪啪动态图| 成人毛片60女人毛片免费| 成年女人永久免费观看视频| 乱系列少妇在线播放| 小蜜桃在线观看免费完整版高清| 亚洲av中文av极速乱| 综合色丁香网| 久久久久久久国产电影| 免费无遮挡裸体视频| 亚洲av.av天堂| 大香蕉97超碰在线| 欧美潮喷喷水| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 精品久久久久久成人av| 国语自产精品视频在线第100页| 精品一区二区免费观看| 一边摸一边抽搐一进一小说| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美精品自产自拍| 一边亲一边摸免费视频| 国产美女午夜福利| 国产精品人妻久久久久久| 91aial.com中文字幕在线观看| 少妇裸体淫交视频免费看高清| 亚洲国产精品久久男人天堂| 人人妻人人澡人人爽人人夜夜 | 女人久久www免费人成看片 | 久久久国产成人精品二区| 欧美人与善性xxx| 国产三级中文精品| 中国国产av一级| 日本黄色片子视频| 伊人久久精品亚洲午夜| 久久亚洲精品不卡| 最近手机中文字幕大全| 国产精品电影一区二区三区| 欧美日本亚洲视频在线播放| 搡老妇女老女人老熟妇| 99热这里只有精品一区| 夜夜看夜夜爽夜夜摸| 亚洲内射少妇av| 麻豆乱淫一区二区| 黄色一级大片看看| 亚洲av.av天堂| 国产真实乱freesex| 人人妻人人看人人澡| 国产午夜精品久久久久久一区二区三区| 国产黄色小视频在线观看| 99在线视频只有这里精品首页| 亚洲天堂国产精品一区在线| 三级经典国产精品| 神马国产精品三级电影在线观看| av在线天堂中文字幕| 久久精品熟女亚洲av麻豆精品 | 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 黄色欧美视频在线观看| 亚洲精品自拍成人| 日韩av在线大香蕉| 国产色爽女视频免费观看| 国产乱人视频| 一夜夜www| 久久人妻av系列| 亚洲av中文av极速乱| 69人妻影院| 高清在线视频一区二区三区 | 久久久成人免费电影| 日韩亚洲欧美综合| 午夜免费男女啪啪视频观看| av国产免费在线观看| 插阴视频在线观看视频| 成人欧美大片| 哪个播放器可以免费观看大片| 午夜老司机福利剧场| 午夜视频国产福利| 色噜噜av男人的天堂激情| 毛片女人毛片| 亚洲无线观看免费| 亚洲怡红院男人天堂| 久久久国产成人免费| 国产伦理片在线播放av一区| 男人舔女人下体高潮全视频| 日本欧美国产在线视频| 中文资源天堂在线| 99热网站在线观看| 欧美xxxx性猛交bbbb| 免费无遮挡裸体视频| 国产国拍精品亚洲av在线观看| 久久草成人影院| 欧美日韩一区二区视频在线观看视频在线 | 久久久久精品久久久久真实原创| 97人妻精品一区二区三区麻豆| 日本免费a在线| 99热这里只有是精品50| 村上凉子中文字幕在线| 性色avwww在线观看| 日日摸夜夜添夜夜添av毛片| 久久99热这里只有精品18| 亚洲人与动物交配视频| 亚洲中文字幕日韩| 国产高清有码在线观看视频| 亚洲在线自拍视频| 哪个播放器可以免费观看大片| 大香蕉97超碰在线| 欧美激情国产日韩精品一区| 日韩大片免费观看网站 | 大又大粗又爽又黄少妇毛片口| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看| 日日摸夜夜添夜夜添av毛片| 最近最新中文字幕免费大全7| 国产69精品久久久久777片| 你懂的网址亚洲精品在线观看 | a级一级毛片免费在线观看| 国产成人91sexporn| 国产精品一二三区在线看| 夜夜看夜夜爽夜夜摸| 18禁在线播放成人免费| 真实男女啪啪啪动态图| 男女下面进入的视频免费午夜| 国产人妻一区二区三区在| h日本视频在线播放| 舔av片在线| 国产成人a区在线观看| 精华霜和精华液先用哪个| 国产精品久久久久久精品电影小说 | 非洲黑人性xxxx精品又粗又长| 国产白丝娇喘喷水9色精品| 色哟哟·www| 高清在线视频一区二区三区 | 成年女人永久免费观看视频| 久久久久久久久久久丰满| 午夜福利在线在线| 国产淫片久久久久久久久| 国产亚洲一区二区精品| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧洲日产国产| 真实男女啪啪啪动态图| 国产精品av视频在线免费观看| 麻豆精品久久久久久蜜桃| 国产免费一级a男人的天堂| 男插女下体视频免费在线播放| www.av在线官网国产| 深爱激情五月婷婷| 赤兔流量卡办理| 国产精品一区二区三区四区久久| 亚洲国产精品久久男人天堂| 亚洲真实伦在线观看| 精品一区二区三区视频在线| 色网站视频免费| 欧美激情在线99| 中国美白少妇内射xxxbb| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 日韩中字成人| 国产精品国产三级国产专区5o | 午夜老司机福利剧场| 少妇熟女aⅴ在线视频| 91在线精品国自产拍蜜月| 97超碰精品成人国产| 床上黄色一级片| 欧美97在线视频| av免费在线看不卡| 国产精品久久电影中文字幕| 深夜a级毛片| 久久精品久久久久久久性| 日本免费一区二区三区高清不卡| 国产亚洲精品久久久com| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 日韩成人伦理影院| 亚洲国产欧美人成| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 国产乱人偷精品视频| 欧美bdsm另类| 又爽又黄无遮挡网站| 亚洲精品乱码久久久久久按摩| 国产片特级美女逼逼视频| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 99久久九九国产精品国产免费| 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 久久99热这里只频精品6学生 | 欧美日韩一区二区视频在线观看视频在线 | 卡戴珊不雅视频在线播放| 日韩,欧美,国产一区二区三区 | 久久久欧美国产精品| 一级毛片电影观看 | av.在线天堂| 亚洲精品国产av成人精品| 欧美最新免费一区二区三区| 色播亚洲综合网| 成年版毛片免费区| 少妇猛男粗大的猛烈进出视频 | 一夜夜www| 亚洲欧美日韩东京热| 日本免费a在线| 亚洲性久久影院| 免费大片18禁| 伦精品一区二区三区| 观看美女的网站| 不卡视频在线观看欧美| 国产黄色小视频在线观看| 麻豆成人av视频| 久久综合国产亚洲精品| 亚洲性久久影院| 国产高清有码在线观看视频| 九草在线视频观看| 国产成人a∨麻豆精品| 亚洲av二区三区四区| 久久午夜福利片| 国产一级毛片七仙女欲春2| 老司机影院成人| 1000部很黄的大片| 国产淫片久久久久久久久| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久久久免| av免费观看日本| 欧美极品一区二区三区四区| 97人妻精品一区二区三区麻豆| 一级av片app| 少妇丰满av| 别揉我奶头 嗯啊视频| 国产精品国产三级国产av玫瑰| av在线老鸭窝| 国产高清视频在线观看网站| 日韩视频在线欧美| 国产一区二区在线av高清观看| 国产精品野战在线观看| a级毛片免费高清观看在线播放| 国产成人午夜福利电影在线观看| 噜噜噜噜噜久久久久久91| 久久亚洲国产成人精品v| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 男女视频在线观看网站免费| 老司机影院毛片| 毛片女人毛片| 成人毛片a级毛片在线播放| 国产麻豆成人av免费视频| 日产精品乱码卡一卡2卡三| 国产视频内射| 哪个播放器可以免费观看大片| 成人亚洲精品av一区二区| 草草在线视频免费看| 国产综合懂色| 欧美精品国产亚洲| 麻豆一二三区av精品|