CN103780193A - Method for detecting initial position of rotor of power robot - Google Patents

Method for detecting initial position of rotor of power robot Download PDF

Info

Publication number
CN103780193A
CN103780193A CN201410063314.6A CN201410063314A CN103780193A CN 103780193 A CN103780193 A CN 103780193A CN 201410063314 A CN201410063314 A CN 201410063314A CN 103780193 A CN103780193 A CN 103780193A
Authority
CN
China
Prior art keywords
rotor
initial position
position angle
theta
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410063314.6A
Other languages
Chinese (zh)
Other versions
CN103780193B (en
Inventor
余海涛
孟高军
胡敏强
黄磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201410063314.6A priority Critical patent/CN103780193B/en
Publication of CN103780193A publication Critical patent/CN103780193A/en
Application granted granted Critical
Publication of CN103780193B publication Critical patent/CN103780193B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

The invention discloses a method for detecting the initial position of a rotor of a power robot. The method comprises the steps of firstly exerting three high-frequency low-voltage pulses on a motor by controlling a driving circuit, utilizing an average current oblique differential method for calculating wire inductance values Lab, Lbc and Lac, comparing the inductance values Lab, Lbc and Lac to determine the section of the position angle of the rotor, controlling an inverter to exert a high-frequency low-voltage pulse again according to the section of the position angle of the rotor, selecting the three optimal detection values in the four current peak value detecting processes for rotational coordinate transformation, and finally determining the initial position angle of the rotor. The method for detecting the initial position of the rotor of a permanent magnet synchronous motor can detect the initial position of the rotor of the permanent magnet synchronous motor quite accurately, and achieve smooth starting of the motor, the implementation process is simple, any motor parameters are not required to be depended on, the measuring result is accurate, and the method has the high engineering practical value.

Description

A kind of Power Robot method for detecting initial position of rotor
Technical field
The present invention relates to a kind of Power Robot method for detecting initial position of rotor, by detecting respectively after A, B, C phase inductance, calculate initial position of rotor by rotating coordinate transformation.
Background technology
At present, in the robot system of various structures, owing to adopting, the scheme efficiency of permagnetic synchronous motor (PMSM) is higher, and therefore this scheme has consequence.Particularly, in Power Robot and small scale robot, PMSM is because these advantages have obtained more application.And can carry out accurately estimating to be the precondition that magneto alternator High Performance Control Strategies (vector control or Direct Torque) and position-sensor-free operation realize to initial position of rotor, also be to be related to whether smoothly starting of robot, and can realize the key issue of maximum torque starting; Therefore, it is one of the focus of engineering technological research and difficulties that initial position of rotor detects always, especially in Power Robot, because what robot carried out is operating as basic is operation with high pressure, the circuit detecting is abnormally dangerous, inaccurate if initial position of rotor detects, and can cause the repercussion of Power Robot to rotate, result is likely because misoperation destroys whole power circuit and robot, the serious high voltage short circuit that even can cause.
Guarantee that motor does not rotate and the method for detection rotor initial position, current achievement in research is di/dt detection method and two kinds of high-frequency signal injections: di/dt detection method need to inject the series of voltage pulse that amplitude is identical, direction is different in motor, the size of detection more corresponding electric current is carried out estimated rotor initial position, this method is feasible, if but expected that accurate initial position of rotor need to apply the potential pulse of multiple different directions, comparatively complicated to the control of inverter; And employing High Frequency Injection, its basic principle is in motor, to inject specific high-frequency voltage signal, then detect current signal corresponding in motor to determine the initial position of rotor, but the method algorithm is comparatively complicated, implement comparatively difficulty, and need the extra hardware circuits such as low pass filter, increased cost.
In engineering, the method for comparatively conventional initial position of rotor is to force positioning mode, and this method need to be before robot operation, make by force robot rotate, this in Power Robot, be very unreliable be also impracticable.In addition some initial position detection methods, all need to rely on robot and drive parameter, and the accuracy that makes to operate reduces a lot.Therefore propose one and do not rely on robot driving parameter, and realize simply, initial position detection method is very necessary accurately.
Summary of the invention
Goal of the invention: in order to overcome the deficiencies in the prior art, the invention provides a kind of Power Robot method for detecting initial position of rotor, enough detect very accurately a kind of permanent-magnetic synchronous motor rotor initial position, realize the smooth starting of motor; Implementation process is simple, does not need to rely on any parameter of electric machine, and measurement result is accurate.
Technical scheme: for achieving the above object, the technical solution used in the present invention is:
A kind of Power Robot method for detecting initial position of rotor, first applies three high frequency low voltage pulses by controlling three-phase inverter (drive circuit) to magneto, and utilizes the oblique poor method of average current to calculate L ab, L bcand L ac, pass through L ab, L bcand L acrelatively determine between rotor position angle location; Then according between rotor position angle location, apply a high frequency low voltage pulse by controlling three-phase inverter to magneto again, in four line inductance detection processes, select three times optimum detected values to be rotated coordinate transform, finally determine initial position angle of rotor; Specifically comprise the steps:
(1) control mode of employing three-phase inverter conducting between two, be that each moment motor has two to be conducted, non-conduction mutually unsettled, apply three high frequency low voltage pulses to magneto, in the time that corresponding two are conducted, calculate corresponding two alternate inductance suc as formula (a) by line inductance computing module:
L ab ( θ r ) = 2 U d di ab / dt | on - di ab / dt | off L ac ( θ r ) = 2 U d di ac / dt | on - di ac / dt | off L bc ( θ r ) = 2 U d di bc / dt | on - di bc / dt | off - - - ( a )
Wherein: U dfor DC bus-bar voltage; L ab, L bcand L acbe respectively AB phase line inductance, BC phase line inductance and AC phase line inductance; Di/dt| onfor the rate of change of Current rise; Di/dt| offfor the rate of change of electric current decline; θ rfor initial position angle of rotor;
According to the saliency of magneto, detect current peak and present following Changing Pattern with the conversion of rotor position angle:
L bc ( θ r ) = L o + L g 2 cos 2 θ r L ac ( θ r ) = L o + L g 2 cos ( 2 θ r + 2 3 π ) L ab ( θ r ) = L o + L g 2 cos ( 2 θ r - 2 3 π )
Wherein: L ofor the mean value of every phase winding self-induction, be the existing parameter of motor, the mean value of each phase winding self-induction equates; L g2for the secondary harmonic amplitude of every phase winding self-induction, be the existing parameter of motor, the secondary harmonic amplitude of each phase winding self-induction equates;
The inductance value L calculating in conjunction with formula (a) ab, L bcand L ac, according to following criterion, tentatively judge initial position angle of rotor θ r:
If L bc>=L ab>L ac, θ rit between the lane place being positioned at, is 0 °~30 ° or 180 °~210 °;
If L ab>=L bc>L ac, θ rit between the lane place being positioned at, is 30 °~60 ° or 210 °~240 °;
If L ab>L ac>=L bc, θ rit between the lane place being positioned at, is 60 °~90 ° or 240 °~270 °;
If L ac>=L ab>L bc, θ rit between the lane place being positioned at, is 90 °~120 ° or 270 °~300 °;
If L ac>L bc>=L ab, θ rit between the lane place being positioned at, is 120 °~150 ° or 300 °~330 °;
If L bc>L ac>=L ab, θ rit between the lane place being positioned at, is 150 °~180 ° or 330 °~360 °;
(2) the initial position angle of rotor θ obtaining based on preliminary judgement r, according to θ rinterval difference, the magnetic direction same/opposite that can cause permanent magnet magnetic extreme direction and winding current to produce, thus make iron core more saturated/reduce saturated, magnetic permeability reduce/increase, line inductance value reduces/becomes greatly, thereby makes L in (a) formula ab, L bcand L acdiminish/become and be large, increase the impact of magnetic or degaussing for fear of electric machine iron core, we are according to θ rinterval again apply a high frequency low voltage pulse to magneto, carry out initial position angle of rotor θ rfinal judgement, rule as follows:
If (2a) initial position angle of rotor θ rinterval be 30 °~90 ° or 210 °~270 °, apply high frequency low voltage pulse to AB phase, C is mutually unsettled: when A phase winding is connected with forward dc busbar voltage, calculate now AB phase line inductance be according to formula (a) when B phase winding is connected with forward dc busbar voltage, calculate now AB phase line inductance according to formula (a) and be
Figure BDA0000469277360000033
Average
Figure BDA0000469277360000034
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 L bc L ac ΣL ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure BDA0000469277360000042
initial position angle of rotor θ rinterval be 30 °~90 °; If
Figure BDA0000469277360000043
initial position angle of rotor θ rinterval be 210 °~270 °;
If (2b) initial position angle of rotor θ rinterval be 0 °~30 °, 330 °~360 ° or 150 °~210 °, apply high frequency low voltage pulse to BC phase, A is mutually unsettled: when B phase winding is connected with forward dc busbar voltage, calculate now BC phase line inductance be according to formula (a)
Figure BDA0000469277360000044
when C phase winding is connected with forward dc busbar voltage, calculate now BC phase line inductance according to formula (a) and be
Figure BDA00004692773600000410
Average
Figure BDA0000469277360000045
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure BDA0000469277360000048
initial position angle of rotor θ rinterval be 150 °~210 °; If
Figure BDA0000469277360000049
initial position angle of rotor θ rinterval be 0 °~30 ° or 330 °~360 °;
If (2c) initial position angle of rotor θ rinterval be 90 °~150 ° or 270 °~330 °, apply high frequency low voltage pulse to AC phase, B is mutually unsettled: when A phase winding is connected with forward dc busbar voltage, calculate now AC phase line inductance be according to formula (a) when C phase winding is connected with forward dc busbar voltage, calculate now AC phase line inductance according to formula (a) and be
Figure BDA0000469277360000057
Average
Figure BDA0000469277360000052
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if initial position angle of rotor θ rinterval be 90 °~150 °; If initial position angle of rotor θ rinterval be 270 °~330 °;
Wherein, L αand L βfor being transformed into the inductance value under alpha-beta coordinate system.
Beneficial effect: Power Robot method for detecting initial position of rotor provided by the invention, can examine very accurately the method for detecting initial position of rotor of Power Robot drive motors, realize the smooth starting of robot; Implementation process is simple, does not need to rely on the parameter of any robot drive motors, and measurement result is accurate, has very strong engineering practical value.
Accompanying drawing explanation
Fig. 1 is implementing procedure figure of the present invention;
Fig. 2 arrives AB phase for applying high frequency low voltage pulse, when C is mutually unsettled, and the equivalent circuit diagram that A phase winding is connected with forward dc busbar voltage;
Fig. 3 arrives AB phase for applying high frequency low voltage pulse, when C is mutually unsettled, and the equivalent circuit diagram that B phase winding is connected with forward dc busbar voltage;
Fig. 4 is initial position angle of rotor while being 45 °, for applying high frequency low voltage pulse to AB phase, and when C is mutually unsettled, increasing magnetic chart when A phase winding is connected with forward dc busbar voltage;
Fig. 5 is initial position angle of rotor while being 45 °, for applying high frequency low voltage pulse to AB phase, and when C is mutually unsettled, increasing magnetic chart when B phase winding is connected with forward dc busbar voltage;
Fig. 6 is line inductance block plan.
Embodiment
Below in conjunction with accompanying drawing, the present invention is further described.
Be illustrated in figure 1 a kind of Power Robot method for detecting initial position of rotor, first apply three high frequency low voltage pulses by controlling three-phase inverter (drive circuit) to magneto, and utilize the oblique poor method of average current to calculate L ab, L bcand L ac, pass through L ab, L bcand L acrelatively determine between rotor position angle location; Then according between rotor position angle location, apply a high frequency low voltage pulse by controlling three-phase inverter to magneto again, in four line inductance detection processes, select three times optimum detected values to be rotated coordinate transform, finally determine initial position angle of rotor; Specifically comprise the steps:
(1) control mode of employing three-phase inverter conducting between two, be that each moment motor has two to be conducted, non-conduction mutually unsettled, apply three high frequency low voltage pulses to magneto, in the time that corresponding two are conducted, calculate corresponding two alternate inductance by line inductance computing module;
As shown in Figure 2, to apply high frequency low voltage pulse to AB phase (T 1and T 4conducting), C is mutually unsettled is example, calculates L ab, order:
L abr)=L aa+L bb+M ab+M ba
Wherein, L aafor A phase self-induction, L bbfor B phase self-induction, M abfor the mutual inductance of A phase and B phase, M bafor the mutual inductance of B phase and A phase;
Work as T 1and T 4conducting, when other power tubes turn-off, when A phase winding is connected with forward dc busbar voltage, B phase winding is while being connected with negative sense DC bus-bar voltage, voltage equation can be expressed as:
U d = 2 U T + 2 R i ab + L ab ( θ r ) di ab dt | on
Work as T 1and T 4turn-off, electric current is by diode VD 2and VD 3afterflow, voltage equation can be expressed as:
- U d = 2 U T + 2 R i ab + L ab ( θ r ) di ab dt | off
Wherein, U dfor DC bus-bar voltage, U ttube voltage drop during for power device conducting; R is the resistance of every phase winding, di/dt| onfor the rate of change of Current rise; Di/dt| offfor the rate of change of electric current decline;
Because power tube is very short from opening to time interval of shutoff, although the instantaneous value of line current is not identical, their mean value is almost identical, therefore, can suppose the U in above-mentioned two voltage equations tchange interval mean value with R pressure drop at power tube on off state and all equate respectively, can obtain so:
L ab ( θ r ) = 2 U d di ab / dt | on - di ab / dt | off
In like manner, can be in the hope of L bcand L ac, obtain formula (a):
L ab ( θ r ) = 2 U d di ab / dt | on - di ab / dt | off L ac ( θ r ) = 2 U d di ac / dt | on - di ac / dt | off L bc ( θ r ) = 2 U d di bc / dt | on - di bc / dt | off - - - ( a )
Wherein: L ab, L bcand L acbe respectively AB phase line inductance, BC phase line inductance and AC phase line inductance; θ rfor initial position angle of rotor;
According to the saliency of magneto, detect current peak and present following Changing Pattern with the conversion of rotor position angle:
L bc ( θ r ) = L o + L g 2 cos 2 θ r L ac ( θ r ) = L o + L g 2 cos ( 2 θ r + 2 3 π ) L ab ( θ r ) = L o + L g 2 cos ( 2 θ r - 2 3 π )
Wherein: L ofor the mean value of every phase winding self-induction, be the existing parameter of motor, the mean value of each phase winding self-induction equates; L g2for the secondary harmonic amplitude of every phase winding self-induction, be the existing parameter of motor, the secondary harmonic amplitude of each phase winding self-induction equates;
The inductance value L calculating in conjunction with formula (a) ab, L bcand L ac, according to following criterion, tentatively judge initial position angle of rotor θ r:
If L bc>=L ab>L ac, θ rit between the lane place being positioned at, is 0 °~30 ° or 180 °~210 °;
If L ab>=L bc>L ac, θ rit between the lane place being positioned at, is 30 °~60 ° or 210 °~240 °;
If L ab>L ac>=L bc, θ rit between the lane place being positioned at, is 60 °~90 ° or 240 °~270 °;
If L ac>=L ab>L bc, θ rit between the lane place being positioned at, is 90 °~120 ° or 270 °~300 °;
If L ac>L bc>=L ab, θ rit between the lane place being positioned at, is 120 °~150 ° or 300 °~330 °;
If L bc>L ac>=L ab, θ rit between the lane place being positioned at, is 150 °~180 ° or 330 °~360 °;
(2) the initial position angle of rotor θ obtaining based on preliminary judgement r, according to θ rinterval difference, the magnetic direction same/opposite that can cause permanent magnet magnetic extreme direction and winding current to produce, thus make iron core more saturated/reduce saturated, magnetic permeability reduce/increase, line inductance value reduces/becomes greatly, thereby makes L in (a) formula ab, L bcand L acdiminish/become and be large, increase the impact of magnetic or degaussing for fear of electric machine iron core, we are according to θ rinterval again apply a high frequency low voltage pulse to magneto, carry out initial position angle of rotor θ rfinal judgement, rule as follows:
If (2a) initial position angle of rotor θ rinterval be 30 °~90 ° or 210 °~270 °, apply high frequency low voltage pulse to AB phase, C is mutually unsettled: when A phase winding is connected (a, b are conducted) with forward dc busbar voltage, calculate now AB phase line inductance to be according to formula (a)
Figure BDA0000469277360000082
when B phase winding is connected (b, a are conducted) with forward dc busbar voltage, calculate now AB phase line inductance according to formula (a) to be
Figure BDA0000469277360000083
Because inductance is subject to increasing the impact of magnetic and degaussing, have influence on the detection of current peak, thereby can cause rotor position angle θ rerror, the therefore error of avoiding in order to try one's best, improve precision, average
Figure BDA0000469277360000084
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 L bc L ac ΣL ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure BDA0000469277360000093
initial position angle of rotor θ rinterval be 30 °~90 °; If
Figure BDA0000469277360000094
initial position angle of rotor θ rinterval be 210 °~270 °; Be analyzed as follows:
If θ rinterval be 30 °~90 °, in the time that a, b are conducted, the resultant magnetic field direction that permanent magnet magnetic extreme direction produces with winding current is identical, produce increase magnetic effect, thereby make iron core more saturated, magnetic permeability reduce, inductance value reduces, thereby makes the L calculating abdiminish; In the time that b, a are conducted, the resultant magnetic field opposite direction that permanent magnet magnetic extreme direction and winding current produce, produce demagnetizing effect, thereby it is saturated that iron core is moved back, and magnetic permeability increases, and inductance value increases, thereby makes the L calculating abbecome large; Therefore
Figure BDA0000469277360000095
If θ rinterval be 210 °~270 °, situation is contrary;
If (2b) initial position angle of rotor θ rinterval be 0 °~30 °, 330 °~360 ° or 150 °~210 °, apply high frequency low voltage pulse to BC phase, A is mutually unsettled: when B phase winding is connected (b, c are conducted) with forward dc busbar voltage, calculate now BC phase line inductance to be according to formula (a)
Figure BDA0000469277360000096
when C phase winding is connected (c, b are conducted) with forward dc busbar voltage, calculate now BC phase line inductance according to formula (a) to be
Figure BDA0000469277360000097
Because inductance is subject to increasing the impact of magnetic and degaussing, have influence on the detection of current peak, thereby can cause rotor position angle θ rerror, the therefore error of avoiding in order to try one's best, improve precision, average
Figure BDA0000469277360000098
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure BDA0000469277360000103
initial position angle of rotor θ rinterval be 150 °~210 °; If
Figure BDA0000469277360000104
initial position angle of rotor θ rinterval be 0 °~30 ° or 330 °~360 °; Be analyzed as follows:
If θ rinterval be 0 °~30 ° or 330 °~360 °, in the time that b, c are conducted, the resultant magnetic field opposite direction that permanent magnet magnetic extreme direction and winding current produce, produce demagnetizing effect, thereby it is saturated that iron core is moved back, magnetic permeability increases, inductance value increases, thereby makes the L calculating bcbecome large; In the time that c, b are conducted, the resultant magnetic field direction that permanent magnet magnetic extreme direction produces with winding current is identical, produces and increases magnetic effect, thereby make iron core more saturated, and magnetic permeability reduces, and inductance value reduces, thereby makes the L calculating bcdiminish; Therefore
Figure BDA0000469277360000105
If θ rinterval be 150 °~210 °, situation is contrary;
If (2c) initial position angle of rotor θ rinterval be 90 °~150 ° or 270 °~330 °, apply high frequency low voltage pulse to AC phase, B is mutually unsettled: when A phase winding is connected (a, c are conducted) with forward dc busbar voltage, calculate now AC phase line inductance to be according to formula (a)
Figure BDA0000469277360000106
when C phase winding is connected (c, a are conducted) with forward dc busbar voltage, calculate now AC phase line inductance according to formula (a) to be
Figure BDA0000469277360000107
Because inductance is subject to increasing the impact of magnetic and degaussing, have influence on the detection of current peak, thereby can cause rotor position angle θ rerror, the therefore error of avoiding in order to try one's best, improve precision, average
Figure BDA0000469277360000108
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure BDA0000469277360000112
initial position angle of rotor θ rinterval be 90 °~150 °; If
Figure BDA0000469277360000113
initial position angle of rotor θ rinterval be 270 °~330 °; Be analyzed as follows:
If θ rinterval be 90 °~150 °, in the time that a, c are conducted, the resultant magnetic field direction that permanent magnet magnetic extreme direction produces with winding current is identical, produce increase magnetic effect, thereby make iron core more saturated, magnetic permeability reduce, inductance value reduces, thereby makes the L calculating acdiminish; In the time that c, a are conducted, the resultant magnetic field opposite direction that permanent magnet magnetic extreme direction and winding current produce, produce demagnetizing effect, thereby it is saturated that iron core is moved back, and magnetic permeability increases, and inductance value increases, thereby makes the L calculating acbecome large; Therefore
Figure BDA0000469277360000114
If θ rinterval be 270 °~330 °, situation is contrary;
Wherein, L αand L βfor being transformed into the inductance value under alpha-beta coordinate system.
Be briefly described below in conjunction with an example.
A method for detecting initial position of rotor for Power Robot, specifically comprises the steps:
(1) apply three high frequency low voltage pulses by controlling three-phase inverter to motor, and utilize the oblique poor method of average current, calculate line inductance value L ab, L bcand L ac, by relatively determining between rotor-position location of three, according to Fig. 6, by comparing L ab, L bcand L actentatively judge θ rlane place between, suppose θ r=45 °, i.e. L ab>L bc>L acbe positioned at interval 2.
Because θ rwhile between the lane place being positioned at being 30 °~90 °, making alive pulse is conducted b, a, T 3, T 4conducting, calculates line inductance compare simultaneously
Figure BDA0000469277360000117
with
Figure BDA0000469277360000116
for the error of avoiding of trying one's best, in order to improve line inductance identification precision, average here according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 L bc L ac ΣL ab = L o cos 2 θ r L o sin 2 θ r
So, rotor position angle θ rcan be expressed as
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Now try to achieve θ r=45 ° or 225 °, compare simultaneously
Figure BDA0000469277360000123
with
Figure BDA0000469277360000124
can obtain
Figure BDA0000469277360000125
θ rbe positioned at 30 °~90 °, i.e. θ r=45 °.
The above is only the preferred embodiment of the present invention; be noted that for those skilled in the art; under the premise without departing from the principles of the invention, can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (1)

1. a Power Robot method for detecting initial position of rotor, is characterized in that: first apply three high frequency low voltage pulses by controlling three-phase inverter to magneto, and utilize the oblique poor method of average current to calculate L ab, L bcand L ac, pass through L ab, L bcand L acrelatively determine between rotor position angle location; Then according between rotor position angle location, apply a high frequency low voltage pulse by controlling three-phase inverter to magneto again, in four line inductance detection processes, select three times optimum detected values to be rotated coordinate transform, finally determine initial position angle of rotor; Specifically comprise the steps:
(1) control mode of employing three-phase inverter conducting between two, be that each moment motor has two to be conducted, non-conduction mutually unsettled, apply three high frequency low voltage pulses to magneto, in the time that corresponding two are conducted, calculate corresponding two alternate inductance suc as formula (a) by line inductance computing module:
L ab ( θ r ) = 2 U d di ab / dt | on - di ab / dt | off L ac ( θ r ) = 2 U d di ac / dt | on - di ac / dt | off L bc ( θ r ) = 2 U d di bc / dt | on - di bc / dt | off - - - ( a )
Wherein: U dfor DC bus-bar voltage; L ab, L bcand L acbe respectively AB phase line inductance, BC phase line inductance and AC phase line inductance; Di/dt| onfor the rate of change of Current rise; Di/dt| offfor the rate of change of electric current decline; θ rfor initial position angle of rotor;
According to the saliency of magneto, detect current peak and present following Changing Pattern with the conversion of rotor position angle:
L bc ( θ r ) = L o + L g 2 cos 2 θ r L ac ( θ r ) = L o + L g 2 cos ( 2 θ r + 2 3 π ) L ab ( θ r ) = L o + L g 2 cos ( 2 θ r - 2 3 π )
Wherein: L ofor the mean value of every phase winding self-induction, be the existing parameter of motor, the mean value of each phase winding self-induction equates; L g2for the secondary harmonic amplitude of every phase winding self-induction, be the existing parameter of motor, the secondary harmonic amplitude of each phase winding self-induction equates;
The inductance value L calculating in conjunction with formula (a) ab, L bcand L ac, according to following criterion, tentatively judge initial position angle of rotor θ r:
If L bc>=L ab>L ac, θ rit between the lane place being positioned at, is 0 °~30 ° or 180 °~210 °;
If L ab>=L bc>L ac, θ rit between the lane place being positioned at, is 30 °~60 ° or 210 °~240 °;
If L ab>L ac>=L bc, θ rit between the lane place being positioned at, is 60 °~90 ° or 240 °~270 °;
If L ac>=L ab>L bc, θ rit between the lane place being positioned at, is 90 °~120 ° or 270 °~300 °;
If L ac>L bc>=L ab, θ rit between the lane place being positioned at, is 120 °~150 ° or 300 °~330 °;
If L bc>L ac>=L ab, θ rit between the lane place being positioned at, is 150 °~180 ° or 330 °~360 °;
(2) the initial position angle of rotor θ obtaining based on preliminary judgement r, according to θ rinterval difference, the magnetic direction same/opposite that can cause permanent magnet magnetic extreme direction and winding current to produce, thus make iron core more saturated/reduce saturated, magnetic permeability reduce/increase, line inductance value reduces/becomes greatly, thereby makes L in (a) formula ab, L bcand L acdiminish/become and be large, according to θ rinterval again apply a high frequency low voltage pulse to magneto, carry out initial position angle of rotor θ rfinal judgement, rule as follows:
If (2a) initial position angle of rotor θ rinterval be 30 °~90 ° or 210 °~270 °, apply high frequency low voltage pulse to AB phase, C is mutually unsettled: when A phase winding is connected with forward dc busbar voltage, calculate now AB phase line inductance be according to formula (a)
Figure FDA0000469277350000022
when B phase winding is connected with forward dc busbar voltage, calculate now AB phase line inductance according to formula (a) and be
Figure FDA0000469277350000023
Average
Figure FDA0000469277350000024
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 L bc L ac ΣL ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure FDA0000469277350000032
initial position angle of rotor θ rinterval be 30 °~90 °; If
Figure FDA0000469277350000033
initial position angle of rotor θ rinterval be 210 °~270 °;
If (2b) initial position angle of rotor θ rinterval be 0 °~30 °, 330 °~360 ° or 150 °~210 °, apply high frequency low voltage pulse to BC phase, A is mutually unsettled: when B phase winding is connected with forward dc busbar voltage, calculate now BC phase line inductance be according to formula (a)
Figure FDA0000469277350000034
when C phase winding is connected with forward dc busbar voltage, calculate now BC phase line inductance according to formula (a) and be
Figure FDA00004692773500000310
Average
Figure FDA0000469277350000035
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if initial position angle of rotor θ rinterval be 150 °~210 °; If
Figure FDA0000469277350000039
initial position angle of rotor θ rinterval be 0 °~30 ° or 330 °~360 °;
If (2c) initial position angle of rotor θ rinterval be 90 °~150 ° or 270 °~330 °, apply high frequency low voltage pulse to AC phase, B is mutually unsettled: when A phase winding is connected with forward dc busbar voltage, calculate now AC phase line inductance be according to formula (a)
Figure FDA0000469277350000041
when C phase winding is connected with forward dc busbar voltage, calculate now AC phase line inductance according to formula (a) and be
Figure FDA0000469277350000047
Average
Figure FDA0000469277350000042
according to motor coordinate transform basic theories, the relation between rest frame A-B-C and alpha-beta coordinate system can be expressed as:
L α L β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 Σ L bc L ac L ab = L o cos 2 θ r L o sin 2 θ r
So, initial position angle of rotor θ rcan be expressed as:
θ r = 1 2 arctan L β L α or 1 2 arctan L β L α + π
Meanwhile, if
Figure FDA0000469277350000045
initial position angle of rotor θ rinterval be 90 °~150 °; If
Figure FDA0000469277350000046
initial position angle of rotor θ rinterval be 270 °~330 °;
Wherein, L αand L βfor being transformed into the inductance value under alpha-beta coordinate system.
CN201410063314.6A 2014-02-25 2014-02-25 A kind of Power Robot method for detecting initial position of rotor Expired - Fee Related CN103780193B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410063314.6A CN103780193B (en) 2014-02-25 2014-02-25 A kind of Power Robot method for detecting initial position of rotor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410063314.6A CN103780193B (en) 2014-02-25 2014-02-25 A kind of Power Robot method for detecting initial position of rotor

Publications (2)

Publication Number Publication Date
CN103780193A true CN103780193A (en) 2014-05-07
CN103780193B CN103780193B (en) 2016-04-06

Family

ID=50572113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410063314.6A Expired - Fee Related CN103780193B (en) 2014-02-25 2014-02-25 A kind of Power Robot method for detecting initial position of rotor

Country Status (1)

Country Link
CN (1) CN103780193B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017107105A1 (en) * 2015-12-23 2017-06-29 罗伯特·博世有限公司 Method and apparatus for online estimation of initial position of surface permanent magnet motor
CN109067283A (en) * 2018-07-27 2018-12-21 江苏大学 A kind of permanent-magnetic synchronous motor rotor initial position identification system and method
CN109462348A (en) * 2018-10-26 2019-03-12 湖南科技大学 The three-phase switch reluctance machine method for controlling position-less sensor and device extracted based on line inductive aspect point
CN112422012A (en) * 2020-09-28 2021-02-26 智新控制系统有限公司 Zero position identification method of zero-speed off-load permanent magnet synchronous motor
CN114268260A (en) * 2022-03-02 2022-04-01 浙江大学 Motor parameter identification method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050140324A1 (en) * 2003-11-27 2005-06-30 Siemens Aktiengesellschaft Method and device for identifying a faulty rotor position angle signal
CN101714844A (en) * 2009-11-10 2010-05-26 哈尔滨工业大学 Method for detecting initial position of magnetic pole of rotor of built-in permanent magnetic synchronous motor
CN102751937A (en) * 2012-07-20 2012-10-24 上海新时达电气股份有限公司 Method for determining initial position angle of rotor of permanent-magnetic synchronous motor

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050140324A1 (en) * 2003-11-27 2005-06-30 Siemens Aktiengesellschaft Method and device for identifying a faulty rotor position angle signal
CN101714844A (en) * 2009-11-10 2010-05-26 哈尔滨工业大学 Method for detecting initial position of magnetic pole of rotor of built-in permanent magnetic synchronous motor
CN102751937A (en) * 2012-07-20 2012-10-24 上海新时达电气股份有限公司 Method for determining initial position angle of rotor of permanent-magnetic synchronous motor

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何栋炜等: "内置式永磁同步电机转子初始位置估计方法", 《电机与控制学报》 *
王冉珺等: "永磁同步电机转子初始位置的检测方法", 《电机与控制学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017107105A1 (en) * 2015-12-23 2017-06-29 罗伯特·博世有限公司 Method and apparatus for online estimation of initial position of surface permanent magnet motor
US10601354B2 (en) 2015-12-23 2020-03-24 Robert Bosch Gmbh Method and apparatus for on-line estimation of initial position of surface permanent magnet electric machine
CN109067283A (en) * 2018-07-27 2018-12-21 江苏大学 A kind of permanent-magnetic synchronous motor rotor initial position identification system and method
CN109462348A (en) * 2018-10-26 2019-03-12 湖南科技大学 The three-phase switch reluctance machine method for controlling position-less sensor and device extracted based on line inductive aspect point
CN112422012A (en) * 2020-09-28 2021-02-26 智新控制系统有限公司 Zero position identification method of zero-speed off-load permanent magnet synchronous motor
CN112422012B (en) * 2020-09-28 2022-06-10 智新控制系统有限公司 Zero position identification method of zero-speed off-load permanent magnet synchronous motor
CN114268260A (en) * 2022-03-02 2022-04-01 浙江大学 Motor parameter identification method and device
CN114268260B (en) * 2022-03-02 2022-07-19 浙江大学 Motor parameter identification method and device

Also Published As

Publication number Publication date
CN103780193B (en) 2016-04-06

Similar Documents

Publication Publication Date Title
CN103780192B (en) Power Robot method for detecting initial position of rotor
US9577560B2 (en) Synchronous machine control device
CN103780193B (en) A kind of Power Robot method for detecting initial position of rotor
US8912739B2 (en) Synchronous machine control apparatus
EP3188359A1 (en) Device and control method for driving sensorless bldc motor
KR102286371B1 (en) Apparatus and Method for controlling temperature changes in a motor
US9236821B2 (en) Magnetic pole position estimating apparatus for electric motor, controlling apparatus for electric motor, and magnetic pole position estimating method for electric motor
EP2693628A1 (en) Method for identifying inductance parameters of synchronous electric machine and realization system thereof
US8395339B2 (en) Motor control device
CN106549618B (en) Electric automobile safety system based on magnetic field angle redundancy algorithm and control method thereof
JP5431465B2 (en) Method for determining inductance of permanent magnet synchronous machine
CN102332862B (en) The evaluation method of low salient pole embedded permanent magnet synchronous rotor initial position
CN103731083B (en) The method for detecting initial position of rotor of wind generator system
EP3876419B1 (en) Demagnetization sensing for permanent magnet synchronous motor drive and method
CN103795306A (en) Rotor position obtaining method and device based on permanent magnet synchronous transmission system
EP2704307A1 (en) Brushless motor control device and brushless motor control method
CN104967386A (en) Permanent magnet synchronous motor parameter identification method, device and control system
CN103236807A (en) Electro-magnetic doubly salient motor rotor position identification method based on three-phase six-state start
JP3687603B2 (en) PM motor magnetic pole position estimation method
JP6002643B2 (en) Control device and AC motor system using the same
US20180241333A1 (en) Electric motor control device
CN103997268A (en) Position-sensor-free control method for transformer substation patrol robot
US10862416B2 (en) Method for aligning a three-phase machine having soft start and three-phase machine
JP2004129410A (en) Magnetic pole position estimating apparatus for motor and and controller
US20230336098A1 (en) Motor control device and motor control method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160406

Termination date: 20200225