• <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è)教育的大學典范——英國謝菲爾德哈勒姆大學
    老司机靠b影院| 久久国产亚洲av麻豆专区| svipshipincom国产片| 欧美中文日本在线观看视频| 亚洲人成伊人成综合网2020| 日韩有码中文字幕| 最近最新免费中文字幕在线| 精品久久久久久久末码| 亚洲国产精品成人综合色| 亚洲欧美一区二区三区黑人| 嫩草影院精品99| 欧美性猛交╳xxx乱大交人| 日本免费a在线| av在线天堂中文字幕| 黄色视频,在线免费观看| 久久中文字幕一级| 亚洲狠狠婷婷综合久久图片| or卡值多少钱| 成在线人永久免费视频| 亚洲av电影在线进入| videosex国产| 久久久国产精品麻豆| 国产亚洲欧美在线一区二区| 长腿黑丝高跟| 亚洲中文av在线| 久久香蕉国产精品| 在线天堂中文资源库| 亚洲av成人av| 亚洲精品国产精品久久久不卡| 亚洲精品国产精品久久久不卡| 色综合站精品国产| av有码第一页| 黑丝袜美女国产一区| av在线天堂中文字幕| 久久中文字幕人妻熟女| 两个人视频免费观看高清| 国产高清有码在线观看视频 | 久久国产精品人妻蜜桃| 黑人欧美特级aaaaaa片| 一级毛片女人18水好多| 日日爽夜夜爽网站| 亚洲av成人一区二区三| 久久久久久久久中文| 麻豆一二三区av精品| 亚洲中文日韩欧美视频| 欧美 亚洲 国产 日韩一| 成人国产综合亚洲| 成人特级黄色片久久久久久久| 欧美亚洲日本最大视频资源| 身体一侧抽搐| 女人高潮潮喷娇喘18禁视频| 亚洲国产毛片av蜜桃av| 日本 av在线| 99热只有精品国产| 亚洲精品美女久久久久99蜜臀| 亚洲中文av在线| 极品教师在线免费播放| 成人午夜高清在线视频 | 欧美性长视频在线观看| 手机成人av网站| 国产一级毛片七仙女欲春2 | 中国美女看黄片| 国产欧美日韩一区二区精品| 很黄的视频免费| 老司机福利观看| 男女那种视频在线观看| 脱女人内裤的视频| 国产亚洲精品av在线| 久久久久国产精品人妻aⅴ院| 亚洲专区国产一区二区| 老熟妇仑乱视频hdxx| 美女大奶头视频| 一区二区日韩欧美中文字幕| 亚洲精品美女久久久久99蜜臀| 听说在线观看完整版免费高清| 日韩欧美免费精品| a级毛片在线看网站| 一区二区三区精品91| 午夜精品在线福利| 免费在线观看亚洲国产| 色综合欧美亚洲国产小说| 特大巨黑吊av在线直播 | 国产不卡一卡二| 免费av毛片视频| 美女高潮到喷水免费观看| 叶爱在线成人免费视频播放| 夜夜爽天天搞| 欧美色视频一区免费| 97超级碰碰碰精品色视频在线观看| av片东京热男人的天堂| 色av中文字幕| 国产在线精品亚洲第一网站| 亚洲精品av麻豆狂野| 超碰成人久久| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一卡2卡三卡4卡5卡| 夜夜看夜夜爽夜夜摸| 国产爱豆传媒在线观看 | 久久国产精品男人的天堂亚洲| 免费av毛片视频| 久久久久久国产a免费观看| 亚洲av熟女| 亚洲成人久久爱视频| 后天国语完整版免费观看| 1024手机看黄色片| 亚洲欧美激情综合另类| 午夜福利成人在线免费观看| 啦啦啦观看免费观看视频高清| 搞女人的毛片| 亚洲国产高清在线一区二区三 | 女生性感内裤真人,穿戴方法视频| 俄罗斯特黄特色一大片| 国产成人影院久久av| 人人妻人人澡人人看| 国产久久久一区二区三区| 啦啦啦观看免费观看视频高清| 欧美绝顶高潮抽搐喷水| 亚洲色图av天堂| 少妇粗大呻吟视频| 亚洲国产看品久久| 国产精品免费一区二区三区在线| 国产亚洲欧美在线一区二区| 国产成人精品久久二区二区91| 99在线视频只有这里精品首页| xxxwww97欧美| 夜夜躁狠狠躁天天躁| 欧美绝顶高潮抽搐喷水| 亚洲精品在线美女| 亚洲精品一区av在线观看| 久久99热这里只有精品18| 极品教师在线免费播放| 在线观看午夜福利视频| 欧美日韩一级在线毛片| 国产单亲对白刺激| 久久午夜综合久久蜜桃| 久久 成人 亚洲| 岛国在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 久久久久亚洲av毛片大全| 日本 av在线| 成人免费观看视频高清| 国产精品久久久久久人妻精品电影| 好男人在线观看高清免费视频 | 在线视频色国产色| 国产精品一区二区三区四区久久 | 亚洲欧美精品综合久久99| 成人一区二区视频在线观看| 无限看片的www在线观看| 中文亚洲av片在线观看爽| av中文乱码字幕在线| 少妇被粗大的猛进出69影院| 国产免费av片在线观看野外av| 99在线人妻在线中文字幕| 搡老妇女老女人老熟妇| av超薄肉色丝袜交足视频| 久久久国产成人免费| 国产区一区二久久| 亚洲男人的天堂狠狠| 人妻丰满熟妇av一区二区三区| 精品国产国语对白av| 变态另类成人亚洲欧美熟女| 悠悠久久av| 亚洲全国av大片| 久久久国产欧美日韩av| 午夜福利一区二区在线看| 亚洲精品粉嫩美女一区| 精品国产乱码久久久久久男人| a级毛片在线看网站| 在线观看免费午夜福利视频| 国产精品亚洲av一区麻豆| 最新美女视频免费是黄的| 一区二区三区国产精品乱码| 在线观看66精品国产| 嫩草影视91久久| 亚洲精品久久成人aⅴ小说| 亚洲精品粉嫩美女一区| 午夜亚洲福利在线播放| 视频在线观看一区二区三区| x7x7x7水蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 欧美最黄视频在线播放免费| 91字幕亚洲| 国产亚洲av高清不卡| 欧美成人午夜精品| 99国产极品粉嫩在线观看| 欧美在线黄色| 亚洲精品美女久久av网站| 叶爱在线成人免费视频播放| 在线观看免费日韩欧美大片| 国产成人欧美在线观看| 免费一级毛片在线播放高清视频| 久热爱精品视频在线9| 国产亚洲精品第一综合不卡| 欧美黄色片欧美黄色片| 国产亚洲精品久久久久5区| 日韩有码中文字幕| 午夜免费鲁丝| 精品国产超薄肉色丝袜足j| 国产不卡一卡二| 国产亚洲精品久久久久5区| 亚洲成人久久爱视频| 制服人妻中文乱码| 亚洲一区中文字幕在线| 国产男靠女视频免费网站| 国产亚洲欧美在线一区二区| 久久午夜亚洲精品久久| 国产片内射在线| 国产精品国产高清国产av| 自线自在国产av| 日本精品一区二区三区蜜桃| 色播亚洲综合网| 亚洲一区高清亚洲精品| 麻豆一二三区av精品| 精品一区二区三区视频在线观看免费| 99热只有精品国产| 日日干狠狠操夜夜爽| 可以在线观看的亚洲视频| 久久国产精品人妻蜜桃| av免费在线观看网站| 日韩成人在线观看一区二区三区| 国产一区二区三区在线臀色熟女| av福利片在线| 黄色a级毛片大全视频| 日本黄色视频三级网站网址| 久热这里只有精品99| 欧美日韩亚洲综合一区二区三区_| 久99久视频精品免费| 人妻丰满熟妇av一区二区三区| 日本成人三级电影网站| 啦啦啦免费观看视频1| 99国产综合亚洲精品| 99热6这里只有精品| 一级毛片精品| 老司机靠b影院| 深夜精品福利| 18禁裸乳无遮挡免费网站照片 | 中文字幕精品免费在线观看视频| 亚洲av成人不卡在线观看播放网| 深夜精品福利| 久久九九热精品免费| 欧美久久黑人一区二区| 婷婷精品国产亚洲av| 国产成人系列免费观看| 日韩视频一区二区在线观看| 久久久国产成人免费| 中文字幕久久专区| 亚洲精品中文字幕在线视频| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 色综合欧美亚洲国产小说| 超碰成人久久| 搡老岳熟女国产| 性色av乱码一区二区三区2| 精品熟女少妇八av免费久了| 怎么达到女性高潮| 精品少妇一区二区三区视频日本电影| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 欧美丝袜亚洲另类 | 精品国产超薄肉色丝袜足j| 日韩欧美一区视频在线观看| 日韩av在线大香蕉| 中出人妻视频一区二区| 在线十欧美十亚洲十日本专区| 搡老岳熟女国产| 男女视频在线观看网站免费 | 怎么达到女性高潮| 身体一侧抽搐| 久久香蕉国产精品| 精品无人区乱码1区二区| 91在线观看av| 男女下面进入的视频免费午夜 | 国产片内射在线| 激情在线观看视频在线高清| 美女大奶头视频| 一个人观看的视频www高清免费观看 | 最新美女视频免费是黄的| 精品国产超薄肉色丝袜足j| 久久精品影院6| 亚洲人成伊人成综合网2020| 中文资源天堂在线| 美女免费视频网站| cao死你这个sao货| 国产精品 欧美亚洲| 中文亚洲av片在线观看爽| 亚洲精品美女久久久久99蜜臀| 可以在线观看的亚洲视频| 黑人操中国人逼视频| 亚洲欧美日韩无卡精品| 亚洲一码二码三码区别大吗| www.精华液| 中文亚洲av片在线观看爽| 十八禁网站免费在线| 人人妻人人澡人人看| 韩国av一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 嫩草影院精品99| 精品久久蜜臀av无| 成人国产综合亚洲| 亚洲,欧美精品.| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲午夜理论影院| 成人午夜高清在线视频 | 亚洲精品国产区一区二| av欧美777| 国产成人精品久久二区二区91| 亚洲国产中文字幕在线视频| 侵犯人妻中文字幕一二三四区| 在线免费观看的www视频| 99riav亚洲国产免费| 一级毛片精品| 91成人精品电影| av超薄肉色丝袜交足视频| 啦啦啦观看免费观看视频高清| 亚洲全国av大片| 亚洲欧美一区二区三区黑人| 日韩欧美在线二视频| 亚洲电影在线观看av| 精品久久久久久久久久免费视频| 少妇 在线观看| 亚洲男人天堂网一区| 亚洲成人久久爱视频| 搡老熟女国产l中国老女人| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 国内精品久久久久精免费| 男人舔奶头视频| 亚洲av第一区精品v没综合| 午夜成年电影在线免费观看| 久久国产精品影院| 亚洲熟妇中文字幕五十中出| 在线观看免费午夜福利视频| av电影中文网址| 日韩欧美一区视频在线观看| 男人的好看免费观看在线视频 | 男女那种视频在线观看| 亚洲国产欧美日韩在线播放| 天堂√8在线中文| 国产成人一区二区三区免费视频网站| 欧美 亚洲 国产 日韩一| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品电影一区二区三区| 99国产极品粉嫩在线观看| 免费女性裸体啪啪无遮挡网站| 黄色丝袜av网址大全| 日韩欧美 国产精品| videosex国产| 可以免费在线观看a视频的电影网站| 白带黄色成豆腐渣| 丰满的人妻完整版| 久久香蕉激情| 人人妻人人看人人澡| 一区福利在线观看| 亚洲成av人片免费观看| 中文字幕av电影在线播放| 最好的美女福利视频网| 1024视频免费在线观看| 中文字幕av电影在线播放| 国产真人三级小视频在线观看| 丝袜人妻中文字幕| 精品高清国产在线一区| 精华霜和精华液先用哪个| 一级片免费观看大全| 50天的宝宝边吃奶边哭怎么回事| 国产精品影院久久| 日本 欧美在线| 高潮久久久久久久久久久不卡| 亚洲精品在线观看二区| 欧美zozozo另类| 啦啦啦韩国在线观看视频| 国产不卡一卡二| 午夜福利视频1000在线观看| 精品国产国语对白av| 久热这里只有精品99| 中文字幕人妻熟女乱码| 国产亚洲精品综合一区在线观看 | 在线视频色国产色| www.www免费av| 国产成人精品久久二区二区免费| 国内揄拍国产精品人妻在线 | 一级毛片精品| 中文字幕人成人乱码亚洲影| 国产熟女午夜一区二区三区| 国产免费男女视频| 黄片小视频在线播放| 精品一区二区三区四区五区乱码| 色老头精品视频在线观看| 亚洲国产欧洲综合997久久, | 国产精品久久久久久亚洲av鲁大| 一个人观看的视频www高清免费观看 | 丁香六月欧美| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 波多野结衣巨乳人妻| 国产蜜桃级精品一区二区三区| 在线观看免费午夜福利视频| 美女国产高潮福利片在线看| 国产午夜精品久久久久久| 99久久综合精品五月天人人| 精品国产超薄肉色丝袜足j| 18禁美女被吸乳视频| 日本一本二区三区精品| 国产精品久久久av美女十八| 制服诱惑二区| 又黄又粗又硬又大视频| 99精品欧美一区二区三区四区| 两人在一起打扑克的视频| 亚洲国产高清在线一区二区三 | 久久久国产成人精品二区| 亚洲精品在线美女| 日日干狠狠操夜夜爽| 一级毛片精品| 在线永久观看黄色视频| 国产亚洲欧美98| 色哟哟哟哟哟哟| 欧美黑人巨大hd| 精品国产国语对白av| 国产精品,欧美在线| 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 久久中文字幕人妻熟女| 国产在线精品亚洲第一网站| 久久国产乱子伦精品免费另类| 99在线视频只有这里精品首页| 色婷婷久久久亚洲欧美| 日韩三级视频一区二区三区| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 制服诱惑二区| 白带黄色成豆腐渣| 极品教师在线免费播放| 国产成人av激情在线播放| 亚洲国产看品久久| 欧美成狂野欧美在线观看| 亚洲午夜精品一区,二区,三区| 哪里可以看免费的av片| 国产成人啪精品午夜网站| 午夜免费成人在线视频| 男人操女人黄网站| 欧美精品啪啪一区二区三区| 国产一区二区三区视频了| 无遮挡黄片免费观看| 麻豆国产av国片精品| 黄色成人免费大全| 成人免费观看视频高清| 性欧美人与动物交配| 国产伦在线观看视频一区| 999精品在线视频| 国产91精品成人一区二区三区| 美女高潮到喷水免费观看| 免费人成视频x8x8入口观看| 国产极品粉嫩免费观看在线| 最近最新中文字幕大全电影3 | 午夜福利18| 亚洲精品久久国产高清桃花| 久久久久免费精品人妻一区二区 | 99在线视频只有这里精品首页| 777久久人妻少妇嫩草av网站| 久久人妻福利社区极品人妻图片| 精品一区二区三区av网在线观看| 巨乳人妻的诱惑在线观看| 99久久国产精品久久久| 一区二区三区激情视频| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一级毛片七仙女欲春2 | 免费无遮挡裸体视频| 91av网站免费观看| 国产精品野战在线观看| 欧美激情 高清一区二区三区| 国产精品 欧美亚洲| 99re在线观看精品视频| 国产视频内射| 欧美一区二区精品小视频在线| 日韩大码丰满熟妇| 欧美不卡视频在线免费观看 | 亚洲av熟女| 亚洲av成人一区二区三| 国产av不卡久久| 91老司机精品| 国产精品久久久人人做人人爽| 中出人妻视频一区二区| 亚洲精品一区av在线观看| 一进一出抽搐动态| 欧美三级亚洲精品| 亚洲精品久久国产高清桃花| 手机成人av网站| 男女之事视频高清在线观看| 久久精品国产综合久久久| 亚洲中文av在线| 国产99久久九九免费精品| 欧美激情 高清一区二区三区| 亚洲国产精品久久男人天堂| 日韩成人在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 身体一侧抽搐| 午夜福利视频1000在线观看| 久久香蕉精品热| 久久天躁狠狠躁夜夜2o2o| 成人国产综合亚洲| 日本a在线网址| 无限看片的www在线观看| 精品久久久久久久毛片微露脸| 成人欧美大片| 变态另类丝袜制服| 脱女人内裤的视频| 波多野结衣高清作品| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦免费观看视频1| 亚洲一区高清亚洲精品| 神马国产精品三级电影在线观看 | 午夜激情福利司机影院| 日本免费一区二区三区高清不卡| 制服诱惑二区| 999久久久精品免费观看国产| 99riav亚洲国产免费| 国产精品98久久久久久宅男小说| 大型av网站在线播放| 欧美日韩精品网址| 满18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 极品教师在线免费播放| 国产亚洲精品av在线| 嫩草影院精品99| 亚洲国产欧美网| 欧美性猛交黑人性爽| 久久久久久久久中文| 亚洲狠狠婷婷综合久久图片| 色综合婷婷激情| 国产精品精品国产色婷婷| 亚洲国产欧美网| 黄色视频,在线免费观看| 国产精品野战在线观看| 波多野结衣高清作品| 人妻丰满熟妇av一区二区三区| 亚洲精品美女久久久久99蜜臀| 日韩大尺度精品在线看网址| 亚洲国产欧洲综合997久久, | 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 久久精品国产亚洲av香蕉五月| 亚洲国产精品久久男人天堂| 日韩欧美一区视频在线观看| 国产成人啪精品午夜网站| 久久精品国产亚洲av香蕉五月| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品电影 | 国产精品香港三级国产av潘金莲| 久久久久国产一级毛片高清牌| 人人妻,人人澡人人爽秒播| 99精品久久久久人妻精品| 好男人电影高清在线观看| x7x7x7水蜜桃| 2021天堂中文幕一二区在线观 | 国产精品影院久久| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 亚洲成国产人片在线观看| 亚洲自拍偷在线| 婷婷丁香在线五月| 日韩欧美一区视频在线观看| 精品久久久久久成人av| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 国产91精品成人一区二区三区| 两个人看的免费小视频| 在线永久观看黄色视频| 欧美大码av| 一二三四社区在线视频社区8| 亚洲精品av麻豆狂野| 国产精品精品国产色婷婷| 国产精品亚洲美女久久久| av视频在线观看入口| 婷婷亚洲欧美| 日韩大码丰满熟妇| 妹子高潮喷水视频| 亚洲人成电影免费在线| 91字幕亚洲| 国产区一区二久久| 美女午夜性视频免费| 熟女少妇亚洲综合色aaa.| 一区福利在线观看| 日韩欧美国产在线观看| 一二三四社区在线视频社区8| 日日摸夜夜添夜夜添小说| 淫秽高清视频在线观看| 国产成+人综合+亚洲专区| 日韩av在线大香蕉| 黑人巨大精品欧美一区二区mp4| 一本一本综合久久| www.自偷自拍.com| 亚洲免费av在线视频| 制服人妻中文乱码| 亚洲av美国av| 高潮久久久久久久久久久不卡| 久久久水蜜桃国产精品网| 国产片内射在线| 亚洲欧美日韩无卡精品| 狂野欧美激情性xxxx| 日本 av在线| 久久久国产成人精品二区| 久久精品人妻少妇| 最新在线观看一区二区三区| 国产又色又爽无遮挡免费看| 级片在线观看| 天堂√8在线中文| 热re99久久国产66热| 日本免费a在线| 成年版毛片免费区| 大香蕉久久成人网| 日韩免费av在线播放| 欧美黄色淫秽网站| 婷婷六月久久综合丁香| 国产成人av激情在线播放|