Nothing Special   »   [go: up one dir, main page]

CN1588793A - Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor - Google Patents

Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor Download PDF

Info

Publication number
CN1588793A
CN1588793A CNA2004100781411A CN200410078141A CN1588793A CN 1588793 A CN1588793 A CN 1588793A CN A2004100781411 A CNA2004100781411 A CN A2004100781411A CN 200410078141 A CN200410078141 A CN 200410078141A CN 1588793 A CN1588793 A CN 1588793A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
mfrac
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
CNA2004100781411A
Other languages
Chinese (zh)
Other versions
CN1283041C (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CNB2004100781411A priority Critical patent/CN1283041C/en
Publication of CN1588793A publication Critical patent/CN1588793A/en
Application granted granted Critical
Publication of CN1283041C publication Critical patent/CN1283041C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

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

Abstract

This invention relates to nonvelocity sensor permanent-magnet synchro-motor-air conditioning compressor system. It is characterized by that nonvelocity sensor vector control is used for permanent-magnet synchromotor, so, overcoming the shortcoming of the large rotating speed pulsation. At the same time, the torque command current composite control method is used to make further lowering rotating speed pulsation of the motor. It overcomes the problem of unsatability of the system during the procedure of regulating PI regulator parameter in traditional vector control system.

Description

Control method of permanent magnet synchronous motor-air conditioner compressor system without speed sensor
Technical Field
The invention belongs to the motor control technology, in particular to the technical field of application of a permanent magnet synchronous motor in a variable frequency air conditioner.
Background
With the coming of energy crisis worldwide, governments of various countries actively promote energy-saving and consumption-reducing technologies for the purpose of economic sustainable development. As one of the main appliances for home electricity, the conventional fixed frequency air conditioner is gradually coming out of the market due to its low operation efficiency. The new generation of inverter air conditioner has a series of advantages of obvious energy-saving effect, stable temperature regulation, low operation noise in the whole frequency range and the like, so that the inverter air conditioner is concerned by the market.
With the continuous improvement and perfection of the performance of the permanent magnet material and the gradual maturity of the research and development experience of the permanent magnet motor, the permanent magnet motor is more and more widely applied in the aspects of national defense, industrial and agricultural production, daily life and the like, and is developed in the aspects of high power, high functionalization and miniaturization. The permanent magnet motor has the characteristics of simple structure, small volume, light weight, low loss, high efficiency and the like, and is gradually valued by the industry.
Conventionally, the Permanent Magnet Motor can be divided into a Permanent Magnet Synchronous Motor (PMSM) and a brushless direct current Motor (BLDCM) according to the difference of the geometric shape of the rotor magnetic steel of the Permanent Magnet Motor. There are many similarities between the two, the biggest difference between them is: when the rotor rotates, the counter electromotive force generated in the stator has different waveforms, and the counter electromotive force of the magnetic synchronous motor is a sine wave, whereas the brushless dc motor is a trapezoidal wave. Therefore, the two motors differ in principle, model, and control method. A control system of a brushless dc motor is widely used in an air conditioner at present.
The vector control system of the permanent magnet synchronous motor can realize high-precision, high-dynamic performance and large-range speed and position control. Compared with a direct current motor, the motor has no mechanical commutator and brush; compared with asynchronous motor, it needs no reactive exciting current, so that it has high power factor, low stator resistance loss, measurable rotor parameter and good control performance. At present, a servo system composed of a permanent magnet synchronous motor is widely applied to the fields of flexible manufacturing systems, robots, office automation, numerical control machines and the like.
For compressor load in an air conditioner, the load torque of the motor periodically fluctuates due to changes in the pressure in the compressor cylinder during each rotation of the motor rotor. In a conventional control system using a brushless dc motor, the stator current of the brushless dc motor is a square wave, and each phase is turned on at 120 ° electrical angle and then turned off at 60 ° electrical angle. The brushless DC motor rotor position detector outputs a pulse every 60 DEG electrical angle because a switch changes state every 60 DEG electrical angle. In general, the rotor position angle of the brushless dc motor is unknown in each 60 ° electrical angle, so that the brushless dc motor correspondingly outputs six reference torque values to approximate the load torque in a six-step manner, and cannot accurately follow the change of the load torque. Fig. 1 is a schematic diagram of a method for controlling the load of a compressor and the torque of a brushless dc motor, and it can be seen that the set torque cannot track the load torque of the compressor well. This inevitably produces large rotational speed fluctuations, which affect the performance of the compressor.
The stator current of the permanent magnet synchronous motor is sine wave, and the position of a rotor needs to be continuously detected. The vector control system of the permanent magnet synchronous motor can realize high-precision, high-dynamic performance and large-range speed and position control, thereby achieving the effect of high-performance control in the variable frequency air conditioner. The vector control system of the permanent magnet synchronous motor can replace the original brushless direct current motor to quickly and accurately track load torque, reduce fluctuation of rotating speed and realize high-performance control. Meanwhile, because the compressor cannot be provided with a speed and position mechanical sensor, the method of using a speed sensor-free method to realize the control of the load of the high-performance compressor becomes a novel control mode of the air-conditioning compressor.
Vector control system of permanent magnet synchronous motor
In 1971, the vector control method of the alternating current motor proposed by Blaschke et al in Germany theoretically solves the problem of high-performance control of the torque of the alternating current motor. The basic idea of vector control stems from the rigorous simulation of dc motors. The direct current motor has good decoupling performance, and can achieve the purpose of controlling the torque of the motor by controlling the armature current and the exciting current of the direct current motor respectively. The ultimate goal of vector control is to improve the torque control performance of the motor, while the implementation still falls on controlling the stator current. The vector control divides the stator current into an excitation component and a torque component through the motor magnetic field orientation, and the excitation component and the torque component are respectively controlled, so that a good decoupling characteristic is obtained. The control method is firstly applied to the asynchronous motor and is quickly transplanted to the synchronous motor. In fact, vector control is easier to implement on permanent magnet synchronous machines. Since the machine does not induce slip frequency currents in the machine during vector control and control is less affected by parameters (mainly rotor parameters). Currently, the vector control technology is widely applied to permanent magnet synchronous motors.
Speed sensorless control
Vector control of PMSM typically controls stator current or voltage by detecting or estimating the position and magnitude of the motor rotor flux. Thus, the torque of the motor is only related to the magnetic flux and the current, and the control method is similar to the control method of the direct current motor, so that high control performance can be obtained. In the control of a conventional permanent magnet synchronous motor, in order to obtain the precise position and speed of the rotor, a mechanical sensor is generally required to be arranged on the shaft of the rotor to measure the speed and position of the motor. These mechanical sensors are often encoders, resolvers and tachogenerators. Mechanical sensors provide the rotor signals required by the motor, but also present problems to the governor system:
(1) the mechanical sensor increases the rotational inertia on the rotor shaft of the motor, increases the space size and the volume of the motor, increases a connecting wire and an interface circuit between the motor and a control system by using the mechanical sensor, enables the system to be easily interfered, and reduces the reliability.
(2) The system is not widely applicable to various occasions due to the limitation of the use conditions of the mechanical sensor, such as temperature, humidity, vibration and the like.
(3) The mechanical sensors and their auxiliary circuits increase the cost of the governor system, the price of some high-precision sensors can even be compared with the price of the motor itself.
In order to overcome the defects brought to the system by using the sensor, a plurality of scholars develop the research of the sensorless permanent magnet synchronous motor control system. The AC speed regulating system without mechanical sensor is to utilize the relevant electric signal in the motor winding to estimate the position and speed of the rotor by proper method to replace the mechanical sensor so as to realize the closed-loop control of the motor.
Control of PMSM-air conditioner compressor system
For a permanent magnet synchronous motor-air conditioner compressor load system without speed sensor vector control, speed fluctuation is caused due to the large-range periodical change of load torque. In the actual motor operation process of the system, parameters of the motor, including resistance, inductance, back electromotive force coefficient and the like, can change due to the change of the operation conditions (such as current, temperature, humidity and the like) of the motor, and meanwhile, due to errors in measurement, the parameters of the real motor cannot be accurately obtained. In addition, in order to reduce the cost, it is desirable to increase the length of the control cycle time as much as possible. The study on the motor parameter and the control period change shows that: for the compressor load in the air conditioner, in the process of each rotation of the motor rotor, due to the change of the pressure in the compressor cylinder, the load torque of the motor can fluctuate periodically, the identification of the rotating speed and the position and the operation performance of the system are easy to be subjected to parameter change, and similar problems also occur when the control period of the system is increased.
In the conventional vector control, the parameters of the PI regulator in the control block diagram are generally adjusted to make the response of the system meet the requirements. In the compressor load system, because of the periodic pulsation of the compressor load, the contradiction between the stability and the response speed in the PI parameter setting method is highlighted: when the PI gain is too small, the identified rotating speed is difficult to track the real rotating speed; the PI gain is too large, and the system is unstable and oscillates. Through the adjustment of the regulator parameters, the estimated speed can stably track the actual speed, however, the response speed of the system is reduced due to the reduction of the gain of the speed regulator. The output of the speed regulator, i.e., the torque command current, does not track the changes in load torque well, resulting in a significant lag of the motor electromagnetic torque with respect to the load torque. FIG. 2 shows the simulation result after the PI parameter of the speed regulator is set when the q-axis inductance fluctuates by 10%. Therefore, in each mechanical cycle of the rotor, the speed cannot be guaranteed to be a constant value, and large fluctuation occurs, so that the working performance of the compressor is influenced.
Disclosure of Invention
The invention aims to provide a control method of a permanent magnet synchronous motor-air conditioner compressor system without a speed sensor and with small fluctuation.
Compared with the brushless direct current motor control system which is widely used in the air conditioning system at present, the invention is characterized in that the speed sensorless vector control is applied to the permanent magnet synchronous motor-air conditioner compressor load system, thereby overcoming the defect that the output torque of the brushless direct current motor can not accurately follow the change of the load torque to cause larger rotating speed pulsation. Meanwhile, according to the characteristic of periodic pulsating load torque of the compressor, the invention provides an improved method for reducing the motor rotating speed pulsation in a permanent magnet synchronous motor-air conditioner compressor load system, namely a torque command current composite control method. The method can effectively overcome the contradiction between system stability and response speed in the parameter setting process of the PI regulator of the traditional vector control system, and improves the control performance of the permanent magnet synchronous motor in the variable frequency air conditioner.
Permanent magnet synchronous motor-air conditioner compressor system without speed sensor
The invention provides a permanent magnet synchronous motor-air conditioner compressor system, and a control schematic block diagram of the system is shown in figure 3. In the system, stator phase current of a permanent magnet synchronous motor is firstly input to a position estimation link through sampling and A/D conversion, the position and the speed of the motor are identified through a speed sensorless algorithm, and then the stator current is divided into an excitation component and a torque component by utilizing an existing vector control method and is respectively controlled. Wherein the exciting current idReference value idrefAnd (0), the inverter adopts a space vector PWM modulation method, and the position estimation link adopts a counter-electromotive-based embedded permanent magnet synchronous motor speed sensorless and position identification algorithm. The system has the innovation points that a permanent magnet synchronous motor vector control system is adopted to replace a brushless direct current motor system, and the rotating speed and the position of the motor are estimated by adopting a speed-sensorless method, so that the load torque is quickly and accurately tracked, the fluctuation of the rotating speed is reduced, and the high-performance control is realized.
Torque command current composite control method
A torque command current composite control method is provided for the parameter setting problem of a PI regulator when parameters change in a permanent magnet synchronous motor-air conditioner compressor load system. Compared with the traditional vector control strategy, the method can effectively reduce the rotation speed pulsation and improve the control performance of the permanent magnet synchronous motor in the air conditioning system.
Because there is a relatively strict correspondence between the torque and the position angle of the compressor load, feed forward compensation for the torque current can be added based on the estimated position. I adopted in the systemdA control method of 0, thus at steady state:
Tem=pn(iqψd-idψq)=pn(iqψ1+(Ld-lq)idiq)=pniqψr
in the formula, TemIs an electromagnetic torque, pnIs the number of pole pairs, id,iqFor stator winding d-q axis currents, #d,ψqIs d-q axis flux linkage, psirFor permanent magnet flux linkage of rotor, Ld,LqIs a d-q axis inductance.
Then, a position-torque table is formed based on the correspondence between the compressor torque and the motor rotor position. In each control period, the torque of the load is found according to the estimated position lookup table, and the required electromagnetic torque is obtained, so that the required torque current is directly calculated as follows:
<math> <mrow> <msup> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mi>em</mi> </msub> <mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> <msub> <mi>&psi;</mi> <mi>r</mi> </msub> </mrow> </mfrac> </mrow> </math>
FIG. 4 is a schematic block diagram of a torque command current compound control method. As shown in the figure, the method introduces the estimated angle into a torque current control link, namely a torque current i through a compressor load torque reference tableq_ref' and then calculating the torque current iq_ref' sum with the output of the speed PI regulator as the torque command current value (i)q_ref+iq_ref') the rest of the system is identical to the existing vector control system. FIG. 3 is a block diagram of a system control modified by the introduction of a torque command current compound control method.
The invention is characterized in that it comprises the following steps in sequence:
(1) system initialization
Inputting a user-provided compressor load torque-rotor position curve table to a Digital Signal Processor (DSP) in the system;
in the DSP, the following are set:
reference value omega of the rotational speed ref0, excitation current reference value id_ref=0;
Back emf estimation constant KeConstant of position estimation KθIs a set value;
the control period T is a set value;
setting proportional and integral constants K of speed regulatorp1,Ki1Proportional and integral constants K of torque current regulatorsp2,Ki2Proportional and integral constants K of excitation current regulatorp3,Ki3
Parameters of the permanent magnet synchronous machine provided by the user: number of pole pairs pnPermanent magnet linkage psi of rotorrDq axis inductance Ld,LqStator resistance R, back emf coefficient KE(ii) a At the same time, the initial value of the dq-axis voltage of the motor is set to zero, i.e., ud(0)=0,uq(0)=0;
(2) DSP detects stator three-phase current and calculates dq axle voltage of motor
DSP detects stator three-phase current i coming from the current transformer, the filter capacitor and the A/D converter on the stator side of the motor in sequencea,ib,ic
For u, n.gtoreq.1d,uqTaking the actual reference dq axis voltage calculated in the last digital control period T in the DSP, namely ud(n)=ud_ref(n-1),uq(n)=uq_ref(n-1);
(3) The DSP adopts a counter-electromotive force-based rotating speed and position identification method of an embedded permanent magnet synchronous motor non-speed sensor to identify the rotating speed and position of the motor in the nth periodAnd position thetaM(n);
Setting: a model motor with exactly the same parameters is built on the y δ axis, which is the estimated rotor axis and has an angular error Δ θ from the dq axis, Δ θ being θ - θMTheta is the angle between the d axis and the reference axis + A, thetaMIs the included angle between the gamma axis and the reference axis + A:
setting initial time, rotation speed <math> <mrow> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>,</mo> </mrow> </math> Position thetaM(0)=0,
Then:
<math> <mrow> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>K</mi> <mi>E</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mi>&theta;</mi> </msub> <mi>sgn</mi> <mo>{</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>}</mo> <msub> <mi>&Delta;i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math>
eM(n)=eM(n-1)-KeΔiδ(n)
Figure A20041007814100104
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>L</mi> <mi>d</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>R</mi> <mo>+</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <msub> <mi>pL</mi> <mi>&gamma;</mi> </msub> </mtd> <mtd> <mo>-</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&delta;</mi> </msub> <mo>-</mo> <mi>p</mi> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;</mi> </msub> <mo>-</mo> <msub> <mi>pL</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> <mtd> <mi>R</mi> <mo>-</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <msub> <mi>pL</mi> <mi>&delta;</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Δθ(n)=θ(n)-θM(n)
p is a differential operator;
when the estimation reaches the steady state, Δ θ ≈ 0, so that the formula Lγ,Lδ,LγδAre respectively approximately equal to LdL q0, said Lγ,Lδ,LγδAre respectively
Self-inductance of the gamma axis, self-inductance of the delta axis and mutual inductance between the gamma axis and the delta axis;
wherein,
n is the number of time periods after discretization,
θMis the estimated rotor position angle and is,
eMis the back-emf of the model motor,
e is the back-emf of the actual motor,
ω is the rotor speed, i.e. <math> <mrow> <mi>&omega;</mi> <mo>=</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mo>,</mo> </mrow> </math>
KeIs a back-emf estimation constant, KθThe position estimation constants are set values;
uγ,uδthe voltages on the gamma, delta axes for a real motor,
iγ(n-1),iδ(n-1) is the current response of the real motor on the gamma and delta axes,
i(n),i(n) is the current response of the model motor established by taking gamma and delta shafts as rotor shafts,
Δiγ(n),Δiδ(n) is the difference between the real current and the current of the model motor on the gamma and delta axes in each period;
(4) the DSP adopts a torque command current composite control method to calculate a dq axis voltage reference value, and the method sequentially comprises the following steps;
(4.1) comparing i in the current coordinate system in the abc coordinate systema,ib,icMultiplying by a coordinate transformation matrix to obtain a current value i under dq axisd、iq、i0
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
(4.2) reference value ω of rotation speed of DSPrefAnd (4) the estimated rotating speed obtained in the step (3) <math> <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math> The difference Δ ω (n) ═ ωref- ω (n) input speed regulator, finding q-axis referenceCurrent iq_ref(n), <math> <mrow> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>,</mo> </mrow> </math> Wherein, Kp1,Ki1Proportional and integral constants of the speed regulator;
(4.3) estimating the rotor position θ according to the step 3M(n) in each control cycle, looking up the rotor position-load torque curve table according to the estimated position, while the electromagnetic torque of the motor is equal to the load torque, thus obtaining the electromagnetic torque TemThereby calculating a required torque current iq_ref′:
<math> <mrow> <msup> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mi>em</mi> </msub> <mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> <msub> <mi>&psi;</mi> <mi>r</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow> </math>
(4.4) DSP to iq_ref′+iq_refWith actual q-axis current iqDifference value Δ i ofq(n)=iq_ref′(n)+iq_ref(n)-iq(n) inputting the torque current regulator to obtain a q-axis reference voltage uq_ref(n),
<math> <mrow> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> <msub> <mi>&Delta;i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>&Delta;i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> </mrow> </math>
Wherein, Kp2,Ki2Proportional and integral constants of the torque current regulator;
(4.5) DSP to id_ref0 and idDifference value Δ i ofd(n)=idref(n)-id(n) inputting the excitation current regulator to obtain d-axis reference voltage ud_ref(n), <math> <mrow> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>3</mn> </mrow> </msub> <msub> <mi>&Delta;i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>&Delta;i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>,</mo> </mrow> </math> Wherein, Kp3,Ki3Proportional and integral constants of the excitation current regulator;
(5) DSP according to uq_ref、uq_refThe voltage vector in the dq0 coordinate system is calculated,
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
(6) the DSP obtains switching signals of the inverter by using a space vector PWM method, obtains the on-off time of each switching tube of the inverter, and sends out corresponding PWM pulses to control the output voltage of the inverter, so as to realize the control of the permanent magnet synchronous motor: in each control period T, the on-off time T of each switching tube in three bridge arms of the invertera,Ta,T0
Figure A20041007814100125
T0=T-Ta-Tb
<math> <mrow> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>=</mo> <msqrt> <msup> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mn>2</mn> </msup> </msqrt> <mo>,</mo> </mrow> </math> <math> <mrow> <mi>&gamma;</mi> <mo>=</mo> <msup> <mi>tan</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mfrac> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> </mfrac> </mrow> </math>
The space of voltage vector distribution is divided into 6 regions, and the voltage u in each control periodα,uβThe resultant space vector magnitude under the α β axis is VsThe argument is gamma; and the action time of the switching vector is determined by V in each regionsT,VaTa,VbTbThe triangle formed is determined, where Va=Vb=EdcT is a control period, EdcIs a dc bus voltage; meanwhile, gamma also determines the area where the switching vector, i.e. the space vector, is located, i.e. the current operating state of the inverter.
By applying the torque command current composite control method, because the feedforward control of the torque command current is added, the torque current required by the torque changing in a large range can be calculated in advance according to the estimated position of the rotor, and meanwhile, the traditional PI controller can also adjust the feedback speed signal, so that the obtained electromagnetic torque can track the load torque quite well. FIG. 9 shows the system simulation result after the torque command current composite control method is adopted, wherein the q-axis inductance fluctuates by 10%, the control period is increased from 100us to 400 us. Therefore, the speed fluctuation of the motor in a steady state is greatly reduced (within +/-5%), and the system estimation algorithm is not influenced by parameters and can be continuously and stably carried out, so that the effectiveness of the method in improving the control performance of the permanent magnet synchronous motor in the variable frequency air conditioner is proved.
Drawings
FIG. 1 is a schematic diagram of a compressor load and brushless DC motor torque control method;
FIG. 2 is a simulation result after a PI parameter of a speed regulator is set when q-axis inductance fluctuates by 10%;
FIG. 3 is a diagram of a PMSM-air conditioner compressor control system;
FIG. 4. Torque command current compound control method;
FIG. 5 actual motor d-q axes and model motor γ - δ axes;
FIG. 6 software system control flow;
FIG. 7. software computation flow chart;
FIG. 8 is a schematic diagram of a space vector PWM method: a. inverter voltage vector and sector schematic; b, vector synthesis schematic diagram;
FIG. 9 shows that the q-axis inductance fluctuates by 10%, the control period is increased from 100us to 400us, and the system simulation result is obtained after the torque command current composite control method is adopted;
FIG. 10 is a block diagram of a hardware system architecture.
The symbols and variables in the figures are illustrated as follows:
a PMSM permanent magnet synchronous motor;
theta estimated rotor position angle;
the motor speed of ω estimation;
ωrefsetting the rotating speed of the motor;
idd-axis electricity of motor under dq0 coordinate systemA current, i.e. an excitation current;
iqthe q-axis current, namely the torque current, of the motor under a dq0 coordinate system;
iq_refa reference torque current;
iq_refa reference torque current obtained by a torque command current composite control method;
id_refan excitation current reference value;
ud_ref,uq_refd-axis reference voltage and q-axis reference voltage of the motor under a dq0 coordinate system respectively;
a conversion module from an abc-dq abc coordinate system to a dq0 coordinate system;
a conversion module from dq-alpha beta dq0 coordinate system to alpha beta 0 coordinate system;
SVPWM is a voltage space vector PWM modulation module, and PI is a proportional integrator module;
ud,uqstator winding d-q axis voltage;
id,iqstator winding d-q axis current;
ψd,ψqd-q axis flux linkage;
ψra rotor permanent magnet flux linkage;
r stator armature resistance;
Ld,Lqa d-q axis inductance;
Tem,TLelectromagnetic torque and load torque;
pnthe number of pole pairs;
moment of inertia of the J-rotor
p differential operator
ω rotor speed (also notedθ is the estimated rotor position angle);
eMis the back emf of the model motor;
e is the back emf of the actual motor.
Introduction of the coordinate system of the motors d, q, 0, α, β, 0:
the d, q, 0 coordinate system is a coordinate system which is placed on the rotor and rotates along with the rotor, the d axis is in the direction of the longitudinal axis of the rotor, the q axis leads the d axis by 90 degrees in electrical angle, and the 0 axis is an imaginary axis for ensuring the reversibility of transformation. The d, q, 0 coordinate system is a stationary coordinate system established on the stator. α, β, 0 can be seen as a d, q, 0 coordinate system with zero rotational speed.
The conversion relation of the voltage, the current and the flux linkage of the d, q, 0 coordinate system and the a, b, c coordinate system is as follows:
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>q</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>c</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
where θ is the rotor position angle.
The conversion relation of the voltage, the current and the flux linkage of the alpha, beta, 0 coordinate system and the a, b, c coordinate system is as follows:
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi></mi> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi></mi> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi></mi> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>c</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow> </math>
Detailed Description
FIG. 10 is a block diagram of the hardware system of the present invention. The experimental hardware of the invention adopts a motor control development system PE-Expert of Myway corporation of Japan, and the platform utilizes a DSP chip TMS320C32 of TI corporation and adopts C language programming. The hardware system mainly comprises a PC, a DSP board, an A/D, D/A converter, a PWM generator and a two-level voltage type inverter specially designed for vector control of the AC motor. The system of the invention detects the current of a PMSM stator loop and the direct current bus voltage of an inverter through sensors, carries out AD conversion by utilizing a development system PE-Expert, and carries out velocity-sensor-free vector control in a DSP thereof, and modules related to coordinate transformation, PI regulation and the like in figure 6. PWM pulses are formed by using a space vector modulation method to control an inverter, so that high-performance control of a permanent magnet synchronous motor in the variable-frequency air conditioner is realized.
The system control flow of the present invention is shown in fig. 6, and can be divided into the following steps:
1. firstly, initializing a software system, and setting a reference value omega of a rotating speedrefEstimation of rotational speedCalculation constant Ke,KθProportional and integral constants K of respective regulators of control period Tp,KiAnd inputting a user-provided compressor load torque-position curve table. Initializing motor parameters: number of pole pairs pnPermanent magnet linkage psi of rotorrD-q axis inductance Ld,LqStator resistance R.
2. Detecting three-phase currents of the stator and calculating the dq-axis voltage of the motor.
In a practical system, the stator current can be detected by a current sensor, and after being conditioned by filtering, the stator current is subjected to A/D conversion and then input into a Digital Signal Processor (DSP) for utilization. The initial value of the dq-axis voltage of the motor is zero, i.e. ud(0)=0,uq(0) 0; for u with n ≧ 1d,uqActual reference dq-axis voltage, i.e. u, calculated for the last digital control period in the DSPd(n)=ud_ref(n-1),uq(n)=uq_ref(n-1),n≥1。
3. And identifying the rotating speed and the position of the motor through a position estimation link.
The system adopts a counter-electromotive-force-based embedded permanent magnet synchronous motor speed sensorless rotation speed and position identification algorithm.
Since the rotor position angle θ is unknown, the precise d-q axis of the machine is not available, and all estimation algorithms can only be based on the estimated position angle θMI.e. assuming a model motor with exactly the same parameters is built on the estimated rotor axis (set as the y-delta axis), as shown in fig. 5, the initial moment <math> <mrow> <mi>&theta;</mi> <mo>=</mo> <mn>0</mn> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mo>=</mo> <mn>0</mn> <mo>.</mo> </mrow> </math>
Due to the existence of the rotor angle error delta theta, the equation of the actual motor on the gamma-delta axis is changed as follows:
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>R</mi> <mo>+</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <msub> <mi>pL</mi> <mi>&gamma;</mi> </msub> </mtd> <mtd> <mo>-</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&delta;</mi> </msub> <mo>-</mo> <mi>p</mi> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;</mi> </msub> <mo>-</mo> <msub> <mi>pL</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> <mtd> <mi>R</mi> <mo>-</mo> <msub> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mi>L</mi> </mrow> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <msub> <mi>pL</mi> <mi>&delta;</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Δθ(n)=θ(n)-θM(n)
Figure A20041007814100153
when the steady state is estimated to be reached, Δ θ is close to 0, so L in the formula (3)γ,Lδ,LγδAre respectively approximately equal to Ld,Lq,0。
In digital control, a discrete model of the motor is required to be used, a control period is recorded as T, and then the current response of the real motor on a gamma-delta axis in each control period is deduced as follows:
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>L</mi> <mi>d</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
meanwhile, the current response of the model motor established by taking the gamma-delta shaft as the rotor shaft is as follows,
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
subtracting the formulas (4) and (5) to obtain the difference between the currents of the real motor and the model motor in each period,
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
in each control period, a constant K is giveneAnd KθThe back emf is recurred and the rotor position angle is corrected as follows.
eM(n)=eM(n-1)-KeΔiδ(n)
<math> <mrow> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>K</mi> <mi>E</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mi>&theta;</mi> </msub> <mi>sgn</mi> <mo>{</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>}</mo> <msub> <mi>&Delta;i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math>
Wherein
Figure A20041007814100165
The motor speed can be estimated as follows:
<math> <mrow> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>T</mi> </mfrac> <mo>{</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>}</mo> <mo>=</mo> <mfrac> <mrow> <msub> <mi>e</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>K</mi> <mi>E</mi> </msub> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>K</mi> <mi>&theta;</mi> </msub> <mi>T</mi> </mfrac> <mi>sgn</mi> <mo>{</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>}</mo> <msub> <mi>&Delta;i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math>
4. and calculating a dq axis voltage reference value by a torque command current composite control method.
According to the estimated rotating speed and position obtained in the step 3, i in a current coordinate system under an abc coordinate systema,ib,icMultiplying by a coordinate transformation matrix to obtain a current value i under dq axisd、iq、i0
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
In the present system, the current and speed regulators both take the form of typical limiting plus PI (proportional integrator) regulators. The mathematical expression for the input and output of a conventional PI regulator is:
<math> <mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <mi>e</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mi>i</mi> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>e</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> </mrow> </math>
e (n) is the regulator input, y (n) is the regulator output, Kp,KiRespectively, proportional and integral constants. Reference value omega of the rotational speedrefDifference Δ ω (n) from the estimated rotation speed to ωref- ω (n) input speedA degree regulator for obtaining q-axis reference current i by proportional-integral operationq_ref(n), <math> <mrow> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>.</mo> </mrow> </math>
And (3) finding the magnitude of the torque of the load according to the estimated position lookup table in each control period by the rotor position estimated in the step (3), and obtaining the magnitude of the required electromagnetic torque, thereby directly calculating the magnitude of the required torque current as follows:
<math> <mrow> <msup> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mi>em</mi> </msub> <mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> <msub> <mi>&psi;</mi> <mi>r</mi> </msub> </mrow> </mfrac> </mrow> </math>
will iq_ref′+iq_refAnd iqDifference value Δ i ofq(n)=iq_ref′(n)+iq_ref(n)-iq(n) inputting the torque current regulator, outputting as q-axis reference voltage u by proportional-integral operationq_ref(n) that is <math> <mrow> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>.</mo> </mrow> </math> id_ref0 and idDifference value Δ i ofd(n)=id_ref(n)-id(n) inputting an excitation current regulator, and outputting as a d-axis reference voltage u through proportional-integral operationd_ref(n), <math> <mrow> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>3</mn> </mrow> </msub> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>&Delta;</mi> <msub> <mi>i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>.</mo> </mrow> </math>
5. The voltage vector in the α β axis is calculated. Conversion from dq0 coordinate system to α β 0 coordinate system according to uq_ref、uq_refThe voltage vector in the dq0 coordinate system is calculated.
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
6. And obtaining a switching signal of the inverter by using a space vector PWM method.
The switching states of the power devices in the three bridge arms of the inverter are respectively three switching states SA,SB,SCAnd (4) showing. Let Si(i ═ a, B, C) ═ 1 denotes inverters a, B, with the power devices on the C arm conducting, SiAnd (i ═ a, B, C) ═ 0 indicates that the power devices below the arms of the inverters a, B, C are on. According to SA,SB,SCThe inverter has 8 operating states in total. Corresponding to 8 operating states, there are 8 voltage space vectors V0-V7, wherein V1-V6 are non-zero vectors and have amplitude Edc(DC bus voltage), V0, V7 are zero vectors and the amplitude is 0. The spatial distribution of V0-V7 is shown in FIG. 8. a.
Dividing the voltage vector distribution space into 6 regions according to the voltage u in each control periodα,uβThe resultant space vector (amplitude of:. alpha. beta.) under the α. beta. axis is obtained <math> <mrow> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>=</mo> <msqrt> <msup> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mn>2</mn> </msup> </msqrt> <mo>,</mo> </mrow> </math> The argument is gamma, the angle of the lobe, <math> <mrow> <mo>(</mo> <mi>tan</mi> <mi>&gamma;</mi> <mo>=</mo> <mfrac> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> </mfrac> <mo>)</mo> </mrow> </math> ) The sector in which the switching vector is selected is suitable. Within each zone VsTwo non-zero vectors V of the region boundary can be utilizeda,VbAnd zero vector V0 or V7 by properly controlling their action time Ta,Tb,T0A reference voltage vector is obtained. As shown in fig. 8.b, V may be included in each regionsT,VaTa,VbTbThe action time T is obtained according to the sine theorem of the trianglea,Tb,T0
Figure A20041007814100182
T0=T-Ta-Tb
Wherein T is a control period, EdcIs the dc bus voltage.
Therefore, the on-off time of each switching tube is obtained, and the corresponding PWM pulse is sent out to control the output voltage of the inverter, so that the control of the permanent magnet synchronous motor is realized.
7. It is determined whether the control routine is suspended. If not, returning to the step 2 to start execution again; if yes, the data display waveform is recorded, and the control program is stopped.
The software calculation flow of the above system is shown in fig. 6.

Claims (1)

1. The control method of the speed sensorless permanent magnet synchronous motor-air conditioner compressor system is characterized by sequentially comprising the following steps of:
(1) system initialization
Inputting a user-provided compressor load torque-rotor position curve table to a Digital Signal Processor (DSP) in the system;
in the DSP, the following are set:
reference value omega of the rotational speedref0, excitation current reference value id_ref=0;
Back emf estimation constant KeConstant of position estimation KθIs a set value;
the control period T is a set value;
setting proportional and integral constants K of speed regulatorp1,Ki1Proportional and integral constants K of torque current regulatorsp2,Ki2Proportional and integral constants K of excitation current regulatorp3,Ki3
Parameters of the permanent magnet synchronous machine provided by the user: number of pole pairs pnPermanent magnet linkage psi of rotorr,dqShaft inductance Ld,LqStator resistance R, back emf coefficient KE(ii) a At the same time, the initial value of the dq-axis voltage of the motor is set to zero, i.e., ud(0)=0,uq(0)=0;
(2) DSP detects stator three-phase current and calculates dq axle voltage of motor
DSP detects stator three-phase current i coming from the current transformer, the filter capacitor and the A/D converter on the stator side of the motor in sequencea,ib,ic
For u with n ≧ 1d,uqTaking the actual reference dq axis voltage calculated in the last digital control period T in the DSP, namely ud(n)=ud_ref(n-1),uq(n)=uq_ref(n-1);
(3) The DSP adopts a counter-electromotive force-based rotating speed and position identification method of an embedded permanent magnet synchronous motor non-speed sensor to identify the rotating speed and position of the motor in the nth periodAnd position thetaM(n);
Setting: a model motor with exactly the same parameters is built on the y δ axis, which is the estimated rotor axis and has an angular error Δ θ from the dq axis, Δ θ being θ - θMTheta is the angle between the d axis and the reference axis + A, thetaMIs the included angle between the gamma axis and the reference axis + A;
setting initial time, rotation speed <math> <mrow> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>O</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>,</mo> </mrow> </math> Position thetaM(0)=0,
Then:
<math> <mrow> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>K</mi> <mi>E</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mi>&theta;</mi> </msub> <mi>sgn</mi> <mo>{</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>}</mo> <msub> <mi>&Delta;i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math>
eM(n)=eM(n-1)-KeΔiδ(n)
Figure A2004100781410002C4
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&Delta;i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mover> <mrow> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mo>&CenterDot;</mo> </mover> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>L</mi> <mi>d</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>M&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>L</mi> <mi>d</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mfrac> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> <mtd> <mn>1</mn> <mo>-</mo> <mfrac> <mi>R</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <mi>T</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <mrow> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>L</mi> <mi>q</mi> </msub> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mi>q</mi> </msub> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>d</mi> </msub> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfrac> <mi>T</mi> <msub> <mi>L</mi> <mi>q</mi> </msub> </mfrac> <msub> <mi>e</mi> <mi>M</mi> </msub> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>R</mi> <mo>+</mo> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <mi>p</mi> <msub> <mi>L</mi> <mi>&gamma;</mi> </msub> </mtd> <mtd> <mo>-</mo> <mover> <mrow> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mo>&CenterDot;</mo> </mover> <msub> <mi>L</mi> <mi>&delta;</mi> </msub> <mo>-</mo> <mi>p</mi> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msub> <mi>L</mi> <mi>&gamma;</mi> </msub> <mo>-</mo> <mi>p</mi> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> </mtd> <mtd> <mi>R</mi> <mo>-</mo> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msub> <mi>L</mi> <mi>&gamma;&delta;</mi> </msub> <mo>+</mo> <msub> <mi>pL</mi> <mi>&delta;</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&delta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>e</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&Delta;&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Δθ(n)=θ(n)-θM(n)
p is a differential operator;
when the estimation reaches the steady state, Δ θ ≈ 0, so that the formula Lγ,Lδ,LγδAre respectively approximately equal to Ld,Lq0, said Lγ,Lδ,LγδAre respectively
Self-inductance of the gamma axis, self-inductance of the delta axis and mutual inductance between the gamma axis and the delta axis;
wherein,
n is the number of time periods after discretization,
θMis the estimated rotor position angle and is,
eMis the back-emf of the model motor,
e is the back-emf of the actual motor,
ω is the rotor speed, i.e. <math> <mrow> <mi>&omega;</mi> <mo>=</mo> <mover> <msub> <mi>&theta;</mi> <mi>M</mi> </msub> <mo>&CenterDot;</mo> </mover> <mo>,</mo> </mrow> </math>
KeIs a back-emf estimation constant, KθThe position estimation constants are set values;
uγ,uδthe voltages on the gamma, delta axes for a real motor,
iγ(n-1),iδ(n-1) is the current response of the real motor on the gamma and delta axes,
i(n),i(n) is the current response of the model motor established by taking gamma and delta shafts as rotor shafts,
Δiγ(n),Δiδ(n) is the difference between the real current and the current of the model motor on the gamma and delta axes in each period;
(4) the DSP adopts a torque command current composite control method to calculate a dq axis voltage reference value, and the method sequentially comprises the following steps;
(4.1) comparing i in the current coordinate system in the abc coordinate systema,ib,icMultiplying by a coordinate transformation matrix to obtain a current value i under dq axisd、iq、i0
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&pi;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
(4.2) reference value ω of rotation speed of DSPrefAnd (4) the estimated rotating speed obtained in the step (3) <math> <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>&theta;</mi> <mo>&CenterDot;</mo> </mover> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </math> The difference Δ ω (n) ═ ωref- ω (n) input to the speed regulator, determining the q-axis reference current iq_ref(n), <math> <mrow> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mi>&Delta;&omega;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>,</mo> </mrow> </math> Wherein, Kp1,Ki1Proportional and integral constants of the speed regulator;
(4.3) estimating the rotor position θ according to the step 3M(n) in each control cycle, looking up the rotor position-load torque curve table according to the estimated position, while the electromagnetic torque of the motor is equal to the load torque, thus obtaining the electromagnetic torque TemThereby calculating a required torque current iq_ref′:
<math> <mrow> <msup> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mi>em</mi> </msub> <mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> <msub> <mi>&psi;</mi> <mi>r</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow> </math>
(4.4) DSP to iq_ref′+iq_refWith actual q-axis current iqDifference value Δ i ofq(n)=iq_ref′(n)+iq_ref(n)-iq(n) inputting the torque current regulator to obtain a q-axis reference voltage uq_ref(n),
<math> <mrow> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> <msub> <mi>&Delta;i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>&Delta;i</mi> <mi>q</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> </mrow> </math>
Wherein, Kp2,Ki2Proportional and integral constants of the torque current regulator;
(4.5) DSP to id_ref0 and idDifference value Δ i ofd(n)=idref(n)-id(n) inputting the excitation current regulator to obtain d-axis reference voltage ud_ref(n), <math> <mrow> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>K</mi> <mrow> <mi>p</mi> <mn>3</mn> </mrow> </msub> <msub> <mi>&Delta;i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>&Delta;i</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mi>T</mi> <mo>,</mo> </mrow> </math> Wherein, Kp3,Ki3Proportional and integral constants of the excitation current regulator;
(5) DSP according to uq_ref、uq_refThe voltage vector in the dq0 coordinate system is calculated,
<math> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mo>_</mo> <mi>ref</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
(6) the DSP obtains the switching signal of the inverter by using a space vector PWM method, obtains the on-off time of each switching tube of the inverter and sends out corresponding signalsThe PWM pulse controls the output voltage of the inverter to realize the control of the permanent magnet synchronous motor: in each control period T, the on-off time T of each switching tube in three bridge arms of the invertera,Ta,T0
Figure A2004100781410004C8
T0=T-Ta-Tb
<math> <mrow> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>=</mo> <msqrt> <msup> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> <mn>2</mn> </msup> </msqrt> <mo>,</mo> </mrow> </math> <math> <mrow> <mi>&gamma;</mi> <mo>=</mo> <msup> <mi>tan</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mfrac> <msub> <mi>u</mi> <mi>&alpha;</mi> </msub> <msub> <mi>u</mi> <mi>&beta;</mi> </msub> </mfrac> </mrow> </math>
The space of voltage vector distribution is divided into 6 regions, and the voltage u in each control periodα,uβThe resultant space vector magnitude under the α β axis is VsThe argument is gamma; and the action time of the switching vector is determined by V in each regionsT,VαTa,VbTbThe triangle formed is determined, where Vα=Vb=EdcT is a control period, EdcIs a dc bus voltage; meanwhile, gamma also determines the area where the switching vector, i.e. the space vector, is located, i.e. the current operating state of the inverter.
CNB2004100781411A 2004-09-17 2004-09-17 Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor Expired - Fee Related CN1283041C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2004100781411A CN1283041C (en) 2004-09-17 2004-09-17 Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2004100781411A CN1283041C (en) 2004-09-17 2004-09-17 Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor

Publications (2)

Publication Number Publication Date
CN1588793A true CN1588793A (en) 2005-03-02
CN1283041C CN1283041C (en) 2006-11-01

Family

ID=34604963

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004100781411A Expired - Fee Related CN1283041C (en) 2004-09-17 2004-09-17 Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor

Country Status (1)

Country Link
CN (1) CN1283041C (en)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100463355C (en) * 2007-03-07 2009-02-18 今创集团有限公司 Direct torque brushless DC servo control system and its working method
CN101025156B (en) * 2006-02-24 2010-04-21 海尔集团公司 DC frequency-changeable compressor driving apparatus and method
CN101247104B (en) * 2008-03-27 2010-06-02 无锡市飞象电子有限公司 Conjecture method of stator magnetic linkage oriented AC motor rotation velocity and rotor position
CN101794118A (en) * 2010-03-05 2010-08-04 华北电力大学 Excitation system parameter identification method based on system decoupling and sequence-optimization genetic algorithm
CN101443995B (en) * 2006-06-29 2011-04-27 三菱电机株式会社 Control device for AC rotary machinery
CN101299590B (en) * 2007-03-13 2011-06-08 通用汽车环球科技运作公司 Anti-windup control for a current regulator of a pulse width modulation inverter
CN102122916A (en) * 2011-04-18 2011-07-13 苏州秉立电动汽车科技有限公司 Compound control method based on vector control system of permanent magnet synchronous motor
CN102347719A (en) * 2010-07-28 2012-02-08 德昌电机(深圳)有限公司 Motor assembly and household appliance with same
CN102386838A (en) * 2011-11-08 2012-03-21 株洲南车时代电气股份有限公司 Motor control system and control method for alternating current transmission system
CN102577093A (en) * 2009-11-30 2012-07-11 爱信艾达株式会社 Control device of motor drive device
CN102594254A (en) * 2012-02-27 2012-07-18 美的集团有限公司 System for reducing noise of permanent magnet synchronous motor and control method
CN102748834A (en) * 2012-07-21 2012-10-24 中山大洋电机股份有限公司 HVAC (heating ventilation and air conditioning) control system of household central air conditioner
CN102843089A (en) * 2012-08-23 2012-12-26 四川长虹电器股份有限公司 Control method of permanent magnet synchronous motor (PMSM) for frequency conversion refrigerator
CN102868342A (en) * 2011-07-06 2013-01-09 尼得科Sr驱动有限公司 Control of electrical machine
CN103780185A (en) * 2012-10-19 2014-05-07 浙江盾安人工环境股份有限公司 Compressor torque control method and torque control device
CN103944480A (en) * 2014-04-21 2014-07-23 许继电气股份有限公司 Rotor magnetic flux linkage distinguishing method for permanent magnet motor
CN106130417A (en) * 2016-07-06 2016-11-16 苏州博联科技有限公司 Low-speed position evaluation method in a kind of permagnetic synchronous motor
CN106549610A (en) * 2017-01-13 2017-03-29 哈尔滨理工大学 It is a kind of to suppress torque pulsation of brushless DC motor control system and its method for suppressing torque ripple
CN107017813A (en) * 2017-05-17 2017-08-04 深圳市瀚路新能源汽车有限公司 Automobile drive motor control device and system
CN107947663A (en) * 2017-11-27 2018-04-20 奥克斯空调股份有限公司 A kind of low-frequency vibration suppressing method and device
CN109379016A (en) * 2018-11-10 2019-02-22 中山大洋电机股份有限公司 Based on vector control without position sensor permanent magnet synchronous motor constant moment of force control method
CN109412474A (en) * 2018-12-12 2019-03-01 上海航天控制技术研究所 A kind of control method of brushless direct current motor sine torque output
US10234165B2 (en) 2012-07-21 2019-03-19 Zhongshan Broad-Ocean Motor Co., Ltd. HVAC control system for household central air conditioning
CN110323986A (en) * 2019-06-04 2019-10-11 广东工业大学 A kind of permanent-magnet synchronous motor rotor position angle evaluation method
CN111102196A (en) * 2020-01-15 2020-05-05 珠海格力电器股份有限公司 Compressor torque compensation method and air conditioner
CN111277189A (en) * 2020-03-25 2020-06-12 海信(山东)空调有限公司 Compressor low-frequency vibration suppression method and system
CN111357188A (en) * 2017-11-16 2020-06-30 雷诺股份公司 Method for controlling a synchronous machine

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101252336B (en) * 2008-03-07 2011-12-28 清华大学 Permanent magnetism synchronous electric machine - compressor system high speed operation control method

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101025156B (en) * 2006-02-24 2010-04-21 海尔集团公司 DC frequency-changeable compressor driving apparatus and method
CN101443995B (en) * 2006-06-29 2011-04-27 三菱电机株式会社 Control device for AC rotary machinery
US8022660B2 (en) 2006-06-29 2011-09-20 Mitsubishi Electric Corporation Control apparatus for AC rotary machine
CN100463355C (en) * 2007-03-07 2009-02-18 今创集团有限公司 Direct torque brushless DC servo control system and its working method
CN101299590B (en) * 2007-03-13 2011-06-08 通用汽车环球科技运作公司 Anti-windup control for a current regulator of a pulse width modulation inverter
CN101247104B (en) * 2008-03-27 2010-06-02 无锡市飞象电子有限公司 Conjecture method of stator magnetic linkage oriented AC motor rotation velocity and rotor position
CN102577093A (en) * 2009-11-30 2012-07-11 爱信艾达株式会社 Control device of motor drive device
CN102577093B (en) * 2009-11-30 2014-08-20 爱信艾达株式会社 Control device of motor drive device
CN101794118A (en) * 2010-03-05 2010-08-04 华北电力大学 Excitation system parameter identification method based on system decoupling and sequence-optimization genetic algorithm
CN102347719A (en) * 2010-07-28 2012-02-08 德昌电机(深圳)有限公司 Motor assembly and household appliance with same
CN102122916A (en) * 2011-04-18 2011-07-13 苏州秉立电动汽车科技有限公司 Compound control method based on vector control system of permanent magnet synchronous motor
CN102122916B (en) * 2011-04-18 2013-02-13 苏州秉立电动汽车科技有限公司 Compound control method based on vector control system of permanent magnet synchronous motor
CN102868342B (en) * 2011-07-06 2016-12-21 尼得科Sr驱动有限公司 The control of motor
CN102868342A (en) * 2011-07-06 2013-01-09 尼得科Sr驱动有限公司 Control of electrical machine
CN102386838A (en) * 2011-11-08 2012-03-21 株洲南车时代电气股份有限公司 Motor control system and control method for alternating current transmission system
CN102386838B (en) * 2011-11-08 2014-11-19 株洲南车时代电气股份有限公司 Motor control system and control method for alternating current transmission system
CN102594254A (en) * 2012-02-27 2012-07-18 美的集团有限公司 System for reducing noise of permanent magnet synchronous motor and control method
US10234165B2 (en) 2012-07-21 2019-03-19 Zhongshan Broad-Ocean Motor Co., Ltd. HVAC control system for household central air conditioning
CN102748834B (en) * 2012-07-21 2014-08-13 中山大洋电机股份有限公司 HVAC (heating ventilation and air conditioning) control system of household central air conditioner
CN102748834A (en) * 2012-07-21 2012-10-24 中山大洋电机股份有限公司 HVAC (heating ventilation and air conditioning) control system of household central air conditioner
CN102843089B (en) * 2012-08-23 2014-11-05 四川长虹电器股份有限公司 Control method of permanent magnet synchronous motor (PMSM) for frequency conversion refrigerator
CN102843089A (en) * 2012-08-23 2012-12-26 四川长虹电器股份有限公司 Control method of permanent magnet synchronous motor (PMSM) for frequency conversion refrigerator
CN103780185A (en) * 2012-10-19 2014-05-07 浙江盾安人工环境股份有限公司 Compressor torque control method and torque control device
CN103944480A (en) * 2014-04-21 2014-07-23 许继电气股份有限公司 Rotor magnetic flux linkage distinguishing method for permanent magnet motor
CN106130417A (en) * 2016-07-06 2016-11-16 苏州博联科技有限公司 Low-speed position evaluation method in a kind of permagnetic synchronous motor
CN106549610A (en) * 2017-01-13 2017-03-29 哈尔滨理工大学 It is a kind of to suppress torque pulsation of brushless DC motor control system and its method for suppressing torque ripple
CN106549610B (en) * 2017-01-13 2018-12-18 哈尔滨理工大学 A kind of inhibition torque pulsation of brushless DC motor control system and method
CN107017813A (en) * 2017-05-17 2017-08-04 深圳市瀚路新能源汽车有限公司 Automobile drive motor control device and system
CN107017813B (en) * 2017-05-17 2024-04-09 合肥环信科技发展有限公司 Control device and system for driving motor for vehicle
CN111357188A (en) * 2017-11-16 2020-06-30 雷诺股份公司 Method for controlling a synchronous machine
CN107947663A (en) * 2017-11-27 2018-04-20 奥克斯空调股份有限公司 A kind of low-frequency vibration suppressing method and device
CN107947663B (en) * 2017-11-27 2020-01-24 奥克斯空调股份有限公司 Low-frequency vibration suppression method and device
CN109379016A (en) * 2018-11-10 2019-02-22 中山大洋电机股份有限公司 Based on vector control without position sensor permanent magnet synchronous motor constant moment of force control method
CN109412474A (en) * 2018-12-12 2019-03-01 上海航天控制技术研究所 A kind of control method of brushless direct current motor sine torque output
CN110323986A (en) * 2019-06-04 2019-10-11 广东工业大学 A kind of permanent-magnet synchronous motor rotor position angle evaluation method
CN110323986B (en) * 2019-06-04 2021-04-06 广东工业大学 Estimation method for rotor position angle of permanent magnet synchronous motor
CN111102196A (en) * 2020-01-15 2020-05-05 珠海格力电器股份有限公司 Compressor torque compensation method and air conditioner
CN111102196B (en) * 2020-01-15 2020-10-27 珠海格力电器股份有限公司 Compressor torque compensation method and air conditioner
CN111277189A (en) * 2020-03-25 2020-06-12 海信(山东)空调有限公司 Compressor low-frequency vibration suppression method and system
CN111277189B (en) * 2020-03-25 2022-02-25 海信(山东)空调有限公司 Compressor low-frequency vibration suppression method and system

Also Published As

Publication number Publication date
CN1283041C (en) 2006-11-01

Similar Documents

Publication Publication Date Title
CN1283041C (en) Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor
CN1278483C (en) Motor control device
CN1251394C (en) Motor control method and its apparatus
CN1463484A (en) Inverter control method and control device
CN1279688C (en) Speed controller for synchronous machine
CN1258257C (en) Electric motor driving device and freezing device using said motor driving device
CN1122356C (en) Synchronous motor driving method, compressor driving method device for these methods, and brushless DC motor driving device
CN1052834C (en) Controlling system for permanent-magnet synchronous electric motor
CN1082111C (en) Washing machine
CN1667942A (en) Rotor position presuming method and apparatus, motor control method, compressor and program
CN1212701C (en) Speed control device for AG electric motor
CN1042081C (en) system for controlling a brushless DC motor
CN1905351A (en) Motor controller, washing machine, air conditioner and electric oil pump
CN1211907C (en) Induction motor driver and its parameter estimation method
CN1516918A (en) Synchronuos reluctance motor control device
CN1748357A (en) Sensorless control system and method for a permanent magnet rotating machine
CN1171378C (en) Device for detecting rotor position of brushless d.c. motor
CN1976212A (en) Vector controller for a permanent magnet synchronous motor, inverter module, and permanent magnet synchronous motor constant display system
CN1437791A (en) Apparatus for controlling synchronous motor
CN1934781A (en) Motor controller
CN1956317A (en) Field weakening vector controller for permanent magnet synchronous motor and module
CN1835383A (en) Controlling device and its regulating method of synchronous motor
CN1819441A (en) Method and apparatus for calculating/controlling power generation torque
CN1533021A (en) Contrl device and control method for AC motor
CN1838518A (en) Power inverter system and method of correcting supply voltage of the same

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20061101

Termination date: 20091019