CN109724312B - Air conditioner compressor rotating speed control method - Google Patents

Air conditioner compressor rotating speed control method Download PDF

Info

Publication number
CN109724312B
CN109724312B CN201811529484.3A CN201811529484A CN109724312B CN 109724312 B CN109724312 B CN 109724312B CN 201811529484 A CN201811529484 A CN 201811529484A CN 109724312 B CN109724312 B CN 109724312B
Authority
CN
China
Prior art keywords
angular velocity
axis
harmonic
filtering
fluctuation
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.)
Active
Application number
CN201811529484.3A
Other languages
Chinese (zh)
Other versions
CN109724312A (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.)
Qingdao Haier Air Conditioner Gen Corp Ltd
Haier Smart Home Co Ltd
Chongqing Haier Air Conditioner Co Ltd
Original Assignee
Qingdao Haier Air Conditioner Gen Corp Ltd
Haier Smart Home Co Ltd
Chongqing Haier Air Conditioner Co Ltd
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 Qingdao Haier Air Conditioner Gen Corp Ltd, Haier Smart Home Co Ltd, Chongqing Haier Air Conditioner Co Ltd filed Critical Qingdao Haier Air Conditioner Gen Corp Ltd
Priority to CN201811529484.3A priority Critical patent/CN109724312B/en
Publication of CN109724312A publication Critical patent/CN109724312A/en
Application granted granted Critical
Publication of CN109724312B publication Critical patent/CN109724312B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

The invention discloses a method for controlling the rotating speed of an air conditioner compressor, which comprises the process of controlling the compressor according to real-time angular speed and moment; controlling the compressor according to the real angular velocity includes: filtering the shaft error to obtain an angular velocity compensation quantity; obtaining the compensated angular velocity output quantity according to the angular velocity compensation quantity; correcting the real-time angular speed according to the compensated angular speed output quantity and controlling a compressor; the control of the compressor according to the torque includes: calculating the difference between the target angular velocity fluctuation amount and the feedback angular velocity amount to obtain a first angular velocity difference value; filtering the first angular velocity difference to obtain a filtered angular velocity; inputting the filtering angular speed to a speed loop regulator to obtain an output torque; obtaining a moment compensation amount according to the first angular speed difference; and obtaining the compensated output torque according to the torque compensation amount and the output torque and controlling the compressor. By applying the invention, the effectiveness of inhibiting the fluctuation of the rotating speed of the compressor can be improved.

Description

Air conditioner compressor rotating speed control method
Technical Field
The invention belongs to the technical field of motor control, particularly relates to a compressor control technology, and more particularly relates to a method for controlling the rotating speed of an air conditioner compressor.
Background
When the compressor used by the air conditioner runs, the compressor is influenced by the working principle and the control technology of the air conditioner serving as a load, so that the load torque of the compressor is extremely unstable, large rotation speed fluctuation is easily caused, and the running of the compressor is not stable. The unstable operation of the compressor can cause the unstable operation of the whole air conditioning system, thereby causing various adverse effects. And unstable operation can also produce great operating noise, can not satisfy relevant noise standard requirement, influences air conditioner and uses the travelling comfort. This phenomenon is particularly serious in a single-rotor compressor.
Although the prior art also has a method for controlling the rotating speed of the compressor, the effect of inhibiting the rotating speed fluctuation is not ideal enough, and the problem of the rotating speed fluctuation of the compressor cannot be fundamentally solved.
Disclosure of Invention
The invention aims to provide a method for controlling the rotating speed of an air conditioner compressor, which improves the effectiveness of the fluctuation suppression of the rotating speed of the compressor.
In order to realize the purpose of the invention, the invention is realized by adopting the following technical scheme:
a method for controlling the rotating speed of an air conditioner compressor comprises a process of controlling the compressor according to a real-time angular speed and a process of controlling the compressor according to a moment;
the process of controlling the compressor according to the real-time angular velocity includes:
acquiring a shaft error Delta theta reflecting the deviation of the actual position and the estimated position of the compressor rotor;
filtering the shaft error delta theta to obtain a corrected shaft error delta theta 'and an angular speed compensation quantity P _ out corresponding to the corrected shaft error delta theta' after at least part of shaft error fluctuation is filtered;
compensating the angular velocity compensation quantity P _ out into an output angular velocity delta omega _ PLL of a phase-locked loop regulator in the phase-locked loop for controlling the compressor to obtain compensated angular velocity output quantity delta omega ', and obtaining delta omega' ═ P _ out + delta omega _ PLL;
correcting the real-time angular speed omega 1 for controlling the compressor according to the compensated angular speed output quantity delta omega', and controlling the compressor according to the corrected real-time angular speed omega 1;
the filtering processing is performed on the axis error Δ θ to obtain a corrected axis error Δ θ 'obtained by filtering out at least part of axis error fluctuation and an angular velocity compensation amount P _ out corresponding to the corrected axis error Δ θ', and specifically includes:
performing Fourier series expansion on the axis error delta theta to obtain the mechanical angle theta of the axis errormThe functional expression of (a);
the function expressions are respectively related to cos thetamnAnd-sin θmnAfter multiplication, extracting d-axis components and q-axis components of n-th harmonic waves of delta theta; thetamnMechanical angle for nth harmonic;
filtering d-axis components and q-axis components of partial harmonic waves by using an integrator to obtain a filtering result, and realizing filtering processing on the axis error delta theta;
respectively comparing the result of filtering out the d-axis component of partial harmonic and the result of filtering out the q-axis component of partial harmonic with cos (theta)mnshift-Pn) And-sin (theta)mnshift-Pn) Multiplying and performing inverse Fourier transform to obtain an angular velocity compensation quantity P _ out corresponding to the correction axis error delta theta' of the filtered partial harmonic component; thetashift-PnPhase compensation angle for nth harmonic;
the process of controlling the compressor according to the torque includes:
calculating the difference between the target angular velocity fluctuation amount and the feedback angular velocity amount to obtain a first angular velocity difference value; the feedback angular velocity amount is the sum of the direct current component P _ DC of the angular velocity compensation amount P _ out and the compensated angular velocity output amount Δ ω';
filtering the first angular velocity difference to obtain a filtered angular velocity at least part of which is filtered out of angular velocity fluctuation, and inputting the filtered angular velocity serving as an input quantity into a speed ring regulator in a speed ring for controlling a compressor to obtain an output torque of the speed ring regulator; meanwhile, executing moment compensation based on the first angular velocity difference to obtain moment compensation amount corresponding to part of angular velocity fluctuation in the first angular velocity difference;
compensating the torque compensation amount to the output torque of the speed ring regulator to obtain the compensated output torque;
and controlling the air-conditioning compressor according to the compensated output torque.
Compared with the prior art, the invention has the advantages and positive effects that: the method for controlling the rotating speed of the air-conditioning compressor provided by the invention comprises the steps of performing fluctuation filtering on a shaft error delta theta reflecting the deviation between the actual position and the estimated position of a compressor rotor, compensating an angular speed compensation quantity corresponding to a corrected shaft error after at least part of shaft error fluctuation is filtered into the output angular speed of a phase-locked loop regulator to obtain a compensated angular speed output quantity, correcting the real-time angular speed of the compressor according to the compensated angular speed output quantity, and enabling the variation quantity and the phase of a target rotating speed to be close to the variation quantity and the phase of the actual rotating speed when the compressor is controlled by the corrected real-time angular speed so as to enable the operation of the compressor to tend to be stable; moreover, because the fluctuation of the shaft error is a front end direct factor causing speed fluctuation, the periodical fluctuation of the shaft error is reduced by filtering the fluctuation of the shaft error at the front end, the speed fluctuation can be directly and quickly inhibited, and the effectiveness of speed control is improved. On the other hand, when extracting the harmonic component in the axis error Δ θ, the phase compensation angle is used to adjust the phase of the harmonic component, and the phase characteristics of the phase-locked loop are changed, so that the fluctuation suppression effect during the full-frequency-domain operation of the compressor can be improved, and the stability of the full-frequency-domain operation can be improved. In addition, the difference value between the feedback angular velocity amount and the target angular velocity fluctuation amount is subjected to filtering processing, and the filtered angular velocity with at least part of angular velocity fluctuation filtered out is input into the speed loop regulator as an input amount, so that the fluctuation of the output torque of the speed loop regulator can be reduced; meanwhile, a moment compensation amount is obtained based on the difference value between the output angular speed of the phase-locked loop regulator and the target angular speed fluctuation amount, the moment compensation amount is compensated into the output moment of the speed loop regulator, a compensated output moment is obtained, the compensated output moment reduces the difference moment between the motor moment and the load moment, and when the compressor is controlled according to the compensated output moment, the fluctuation of the rotating speed of the compressor can be obviously reduced based on the output moment of the speed loop regulator with reduced fluctuation and the compensated output moment with reduced difference moment, so that the compressor can run more stably; the compressor operates stably, and the effects of energy conservation and vibration reduction can be achieved.
Other features and advantages of the present invention will become more apparent from the following detailed description of the invention when taken in conjunction with the accompanying drawings.
Drawings
FIG. 1 is a partial flow chart of an embodiment of a method for controlling the rotational speed of an air conditioner compressor according to the present invention;
FIG. 2 is another partial flowchart of an embodiment of a method for controlling the rotational speed of an air conditioner compressor according to the present invention;
FIG. 3 is a control block diagram based on the method embodiment of FIGS. 1 and 2;
FIG. 4 is a logic block diagram of a specific example of the axis error fluctuation filtering algorithm of FIG. 3;
FIG. 5 is a logic block diagram of a specific example of the speed fluctuation extraction algorithm of FIG. 3;
FIG. 6 is a logic block diagram of a specific example of the torque compensation algorithm of FIG. 3.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings and examples.
Fig. 1 and 2 respectively show partial flowcharts of an embodiment of a method for controlling the rotational speed of an air conditioner compressor according to the present invention. Specifically, the rotational speed control method of the embodiment includes two processes: one is a process for controlling the compressor according to the real-time angular velocity, and the flow chart is shown in fig. 1; one is a process of controlling the compressor according to the torque, and the flowchart is shown in fig. 2. Specific implementations of these two processes are described below based on a control block diagram shown in fig. 1 and 2 in conjunction with fig. 3.
Referring to fig. 1, a partial flowchart of an embodiment of a method for controlling a rotational speed of an air conditioning compressor according to the present invention, specifically a flowchart for controlling a compressor according to a real-time angular velocity, is shown, and the embodiment adopts a process including the following steps to control a compressor according to a real-time angular velocity:
step 11: a shaft error Delta theta reflecting a deviation between an actual position and an estimated position of a compressor rotor is obtained.
In the control of the compressor, the phase of the compressor rotor can be locked to the target phase by a phase-locked loop (PLL) control technique, the control block of which is shown in fig. 3. In the prior art, a phase-locked loop regulator, typically a proportional-integral regulator, is included in the phase-locked loop of the compressor, see K of fig. 3P_PLLAnd KI_PLLand/S. Wherein, KP_PLL、KI_PLLIs the closed loop gain parameter of the phase locked loop. The axis error Δ θ is used as an input of the PLL regulator, and specifically, the axis error Δ θ is subtracted from a target angular fluctuation amount (0 shown in fig. 3), and the difference is input to the PLL regulator, and the output of the PLL regulator is an output angular velocity Δ ω _ PLL. Based on the output angular velocity Δ ω _ PLL of the phase-locked loop regulator, the phase-locked loop outputs a real-time angular velocity ω 1 for compressor control, and the rotor position is controlled using the real-time angular velocity ω 1.
The shaft error Δ θ, which reflects the deviation between the actual position and the estimated position of the compressor rotor, can be calculated by the following equation:
Figure GDA0003067163180000051
in the formula, the first step is that,
Figure GDA0003067163180000052
and
Figure GDA0003067163180000053
respectively a d-axis voltage set value and a q-axis voltage set value of the compressor, IdAnd IqReal-time d-axis current and real-time q-axis current, r, of the compressor, respectively*Is the resistance of the motor of the compressor,
Figure GDA0003067163180000054
is the q-axis inductance, omega, of the compressor1Is the real-time angular frequency of the compressor. Among the parameters, Id、IqAnd ω1The detection is carried out in real time by the detection means in the prior art, and other parameter values are known values.
Step 12: and filtering the axis error delta theta to obtain a corrected axis error delta theta' after at least part of axis error fluctuation is filtered.
Since the shaft error is used as an input to the phase locked loop, the real-time angular velocity of the compressor at the output of the phase locked loop is affected. If the shaft error fluctuation is large, the real-time angular speed output by the phase-locked loop is unstable, so that the rotor phase locking is unstable, and further, the compressor has faults of overcurrent, step loss and the like.
After the axial error Δ θ is obtained in step 11, filtering is performed on the axial error Δ θ to filter at least a part of fluctuation components, and a corrected axial error Δ θ' after at least a part of axial error fluctuation is filtered is obtained.
Step 13: an angular velocity compensation amount P _ out is obtained from the corrected axis error Deltatheta'.
The implementation of the above step 12 and step 13 is reflected in the control block diagram of fig. 3, and an axis error Δ θ fluctuation filtering algorithm is adopted to obtain the angular velocity compensation amount P _ out.
Specifically, the method for filtering the axis error Δ θ to obtain a corrected axis error Δ θ 'obtained by filtering out at least part of axis error fluctuation and an angular velocity compensation amount P _ out corresponding to the corrected axis error Δ θ' includes:
performing Fourier series expansion on the axis error delta theta to obtain the mechanical angle theta of the axis errormThe functional expression of (a);
respectively connecting the function expressions with cos thetamnAnd-sin θmnAfter multiplication, extracting d-axis components and q-axis components of n-th harmonic waves of delta theta; thetamnMechanical angle for nth harmonic;
filtering d-axis components and q-axis components of partial harmonic waves by using an integrator to obtain a filtering result, and realizing filtering processing on axis errors delta theta;
respectively comparing the result of filtering out d-axis component of partial harmonic and the result of filtering out q-axis component of partial harmonic with cos (theta)mnshift-Pn) And-sin (theta)mnshift-Pn) Multiplying and performing inverse Fourier transform to obtain an angular velocity compensation quantity P _ out corresponding to the correction axis error delta theta' of the filtered partial harmonic component; thetashift-PnThe phase compensation angle for the nth harmonic.
The more specific filtering process is described in detail later with reference to fig. 4.
Step 14: the angular velocity compensation amount P _ out is compensated to the output angular velocity delta omega _ PLL of the phase-locked loop regulator in the phase-locked loop for controlling the compressor, and the compensated angular velocity output amount delta omega' is obtained. Specifically, the compensated angular velocity output amount Δ ω' ═ P _ out + Δ ω _ PLL.
Step 15: and correcting the real-time angular speed omega 1 for controlling the compressor according to the compensated angular speed output quantity, and controlling the compressor according to the corrected real-time angular speed omega 1.
Specifically, the method of determining the real-time angular velocity corresponding to the target angular velocity fluctuation amount of 0 in the following velocity loop control is: referring to fig. 3, the compensated angular velocity output amount Δ ω' is added to the angular velocity command ω × in, and the real-time angular velocity ω 1 for controlling the compressor is output. The angular velocity command ω _ in is a given angular velocity value of the compressor control system, and the determination method of the value of the given angular velocity command ω _ in is implemented by using the prior art. The target angular velocity fluctuation quantity of the speed loop is 0, and the real-time angular velocity is determined based on the output angular velocity delta omega _ PLL of the phase-locked loop regulator and the given angular velocity command omega _ in, so that the compressor is controlled more accurately and stably.
Referring to a partial flow chart of an embodiment of the method for controlling the rotation speed of the air conditioner compressor according to the present invention shown in fig. 2, specifically, a flow chart for controlling the compressor according to the torque, the embodiment adopts a process including the following steps:
step 21: and calculating the difference between the target angular velocity fluctuation amount and the feedback angular velocity amount to obtain a first angular velocity difference value. The feedback angular velocity amount is the sum of the direct current component P _ DC of the angular velocity compensation amount P _ out and the compensated angular velocity output amount Δ ω'.
In compressor control, the rotational speed of the compressor rotor can be controlled to approach a set rotational speed by a speed loop (ASR) control technique. Referring to the block diagram of FIG. 3, the speed loop includes a speed loop regulator, typically a proportional integral regulator, see K of FIG. 3P_ASRAnd KI_ASR/S。
In this step, both the direct current component P _ DC of the angular velocity compensation amount P _ out and the compensated angular velocity output amount Δ ω' are used as inputs to the velocity loop. Specifically, a direct current component P _ DC of the angular velocity compensation amount P _ out is extracted, and the sum of the direct current component P _ DC and the compensated angular velocity output amount Δ ω 'is calculated to obtain a feedback angular velocity amount Δ ω 1, where Δ ω 1 is P _ DC + Δ ω'. The extraction of the DC component P _ DC of the angular velocity compensation amount P _ out can be realized by using a conventional technique for extracting a DC component, for example, by using a low-pass filter to extract the DC component P _ DC of the angular velocity compensation amount P _ out.
Then, a difference between the target angular velocity fluctuation amount and the feedback angular velocity amount Δ ω 1 is calculated, and the difference therebetween is determined as a first angular velocity difference Δ ω 2. Here, the target angular velocity fluctuation amount is a desired angular velocity fluctuation amount and is a known input amount. As a preferred embodiment, in this example, the target angular velocity fluctuation amount is 0.
Step 22: filtering the first angular velocity difference to obtain a filtered angular velocity at least part of which is filtered out of angular velocity fluctuation, and inputting the filtered angular velocity as an input quantity to a speed ring regulator in a speed ring for controlling the compressor to obtain an output torque of the speed ring regulator; and meanwhile, executing moment compensation based on the first angular velocity difference to obtain moment compensation amount corresponding to part of angular velocity fluctuation in the first angular velocity difference.
The first angular velocity difference Δ ω 2 is used as an input to the velocity loop regulator, affecting the output torque at the velocity loop output. If the first angular speed difference value fluctuates greatly, the fluctuation of the output torque is large, and further the fluctuation of the rotating speed of the compressor is large. After the first angular velocity difference is obtained in step 21, filtering is performed on the first angular velocity difference to filter out at least part of the angular velocity fluctuation component, so as to obtain a filtered angular velocity Δ ω _ K. The method for filtering the angular velocity can be implemented by adopting a filtering mode in the prior art, and more preferably, the filtering process is described in the following preferred embodiment. Then, the filtered angular velocity Δ ω _ K is input to the speed loop regulator as an input amount, and the output torque τ _ ASR of the speed loop regulator is obtained.
Meanwhile, a moment compensation algorithm is adopted, and moment compensation is performed based on the first angular velocity difference delta omega 2, so that moment compensation amount tau _ out corresponding to part of angular velocity fluctuation in the first angular velocity difference delta omega 2 is obtained. For the moment compensation algorithm, all possible solutions existing in the prior art may be adopted as long as it is ensured that the obtained moment compensation amount τ _ out corresponds to a part of angular velocity fluctuations in the first angular velocity difference Δ ω 2. The preferred moment compensation algorithm is described in the following preferred embodiments.
Step 23: and compensating the torque compensation amount to the output torque of the speed ring regulator to obtain the compensated output torque.
Specifically, the torque compensation amount τ _ out is added to the output torque τ _ ASR of the speed loop regulator to obtain the compensated output torque τM:τM=τ_out+τ_ASR。
Step 24: and controlling the air-conditioning compressor according to the compensated output torque. The specific control process refers to the prior art.
By adopting the method of the embodiment formed by the above-mentioned fig. 1 and fig. 2, the double loop control of the speed loop and the phase-locked loop of the compressor is realized. In addition, in the phase-locked loop control, the fluctuation filtering is carried out on the shaft error delta theta reflecting the deviation between the actual position and the estimated position of the compressor rotor, the angular speed compensation quantity corresponding to the corrected shaft error after at least part of the shaft error fluctuation is filtered is compensated to the output angular speed of the phase-locked loop regulator, the compensated angular speed output quantity is obtained, the real-time angular speed of the compressor is corrected according to the compensated angular speed output quantity, and when the compressor is controlled by the corrected real-time angular speed, the variation quantity and the phase of the target rotating speed can be close to the variation quantity and the phase of the actual rotating speed, so that the operation of the compressor tends to be stable. Moreover, because the fluctuation of the shaft error is a front end direct factor causing speed fluctuation, the periodical fluctuation of the shaft error is reduced by filtering the fluctuation of the shaft error at the front end, the speed fluctuation can be directly and quickly inhibited, and the effectiveness of speed control is improved. In the control of the speed loop, the difference value between the feedback angular velocity quantity and the target angular velocity fluctuation quantity is subjected to filtering processing, and the filtered angular velocity with at least part of angular velocity fluctuation filtered out is input into the speed loop regulator as an input quantity, so that the fluctuation of the output torque of the speed loop regulator can be reduced; obtaining a moment compensation amount based on a difference value between the output angular speed of the phase-locked loop regulator and the target angular speed fluctuation amount, compensating the moment compensation amount into the output moment of the speed loop regulator to obtain a compensated output moment, wherein the compensated output moment can reduce the difference moment between the motor moment and the load moment; then, when the compressor is controlled according to the compensated output torque, the fluctuation of the rotational speed of the compressor can be significantly reduced based on the output torque of the speed loop regulator with reduced fluctuation and the compensated output torque with reduced difference torque, so that the operation of the compressor tends to be smooth. The compressor operates stably, the technical effects of energy conservation and vibration reduction can be achieved, and the operation performance of the compressor is further improved.
In some other embodiments, the filtering process is performed on the axis error Δ θ to obtain a corrected axis error Δ θ' after at least part of the axis error fluctuation is filtered, and specifically includes: and performing filtering processing on the axis error delta theta, at least filtering d-axis components and q-axis components of first harmonic in the delta theta, realizing filtering of first harmonic components of the delta theta, and obtaining a corrected axis error delta theta' of which at least the first harmonic components are filtered. As a more preferable embodiment, the method for obtaining the corrected axis error Δ θ' after filtering at least part of the axis error fluctuation by filtering the axis error Δ θ further includes: and filtering d-axis components and q-axis components of second harmonic in the delta theta, realizing filtering of first harmonic components and second harmonic components of the delta theta, and obtaining corrected axis errors delta theta' for filtering the first harmonic components and the second harmonic components. Most of fluctuation components in the delta theta can be filtered out by filtering out the first harmonic component or the first harmonic component and the second harmonic component, the calculated amount is moderate, and the filtering speed is high.
Fig. 4 is a logic diagram showing a specific example of the axis error fluctuation filtering algorithm in fig. 3, specifically, a logic diagram showing a specific example of obtaining the angular velocity compensation amount P _ out corresponding to the corrected axis error Δ θ' after filtering the first harmonic component and the second harmonic component in the axis error Δ θ. According to the logic diagram shown in fig. 4, the specific process of the angular velocity compensation amount P _ out corresponding to the corrected axis error Δ θ' after filtering the first harmonic component and the second harmonic component in the axis error Δ θ is as follows:
firstly, Fourier series expansion is carried out on the axis error delta theta to obtain the axis error delta theta relative to the mechanical angle thetamIs used for the functional expression of (1). The method comprises the following specific steps:
Figure GDA0003067163180000091
wherein, Delta thetaDCIs the direct component of the axis error, θd_n=θpeak_ncosφn,θq_n=θpeak_nsinφn
Figure GDA0003067163180000092
△θpeak_nFor the n harmonic axis error fluctuation amplitude, thetam1Is the first harmonic mechanical angle, and the second harmonic mechanical angle thetam2Expressed as: thetam2=2θm1
And then, extracting a first harmonic component and a second harmonic component from the function expression, and filtering the first harmonic component and the second harmonic component by adopting an integrator to obtain a filtering result.
Specifically, the first harmonic component and the second harmonic component can be extracted from the functional expression by a low-pass filtering method or an integration method. With particular reference to FIG. 4, the functional expressions are each related to cos θm1And cos θm2After multiplication, a low-pass filter is used for filtering or an integrator is used for taking an integral average value in a period, and a d-axis component of a first harmonic and a d-axis component of a second harmonic of an axis error delta theta are extracted; respectively comparing the function expressions with-sin thetam1And-sin θm2After multiplication, the q-axis component of the first harmonic and the q-axis component of the second harmonic of the axis error delta theta are extracted by filtering through a low-pass filter or taking an integral average value in a period through an integrator. Then, the d-axis component and the q-axis component of the first harmonic and the d-axis component and the q-axis component of the second harmonic are respectively subtracted from 0, and the resultant is input to an integrator KI_PAnd performing integral filtering treatment in the/S, filtering d-axis components and q-axis components of the first harmonic and d-axis components and q-axis components of the second harmonic, obtaining filtering results of the first harmonic component and the second harmonic component, and realizing filtering treatment on the axis error delta theta. Also, the filtering result becomes an angular velocity.
Then, each filtering result is subjected to inverse Fourier transform, and an angular velocity compensation amount P _ out corresponding to the correction axis error Delta theta' of the first harmonic component and the second harmonic component which are filtered is obtained. Specifically, the filtering result of the d-axis component and the filtering result of the q-axis component are respectively equal to cos (θ)m1shift-P1) And-sin (theta)m1shift-P1) Multiplying the sum of the results after the Fourier inverse transformation to form an angular velocity compensation quantity P _ out1 corresponding to the correction axis error with the first harmonic component filtered; the filtering result of the d-axis component and the q-axis component are respectively compared with cos (theta)m2shift-P2) And-sin (theta)m2shift-P2) After multiplication and inverse Fourier transformForming an angular velocity compensation amount P _ out2 corresponding to the correction axis error with the second harmonic component filtered out; the sum of the two angular velocity compensation amounts forms an angular velocity compensation amount P _ out of P _ out1+ P _ out2 corresponding to the correction axis error Δ θ' from which the first harmonic component and the second harmonic component are filtered out. Wherein, thetashift-P1And thetashift-P2The phase compensation angle of the first harmonic and the phase compensation angle of the second harmonic are respectively. The angle numbers of the two phase compensation angles can be equal or unequal preset fixed values, and can also be variable angle values.
As a preferred embodiment, two phase compensation angles θshift-P1And thetashift-P2Equal and according to the closed-loop gain parameter K of the phase-locked loopP_PLL、KI_PLLAnd determining the angular speed command omega _ in of the phase-locked loop. Furthermore, it is necessary to satisfy: thetashift-Pn=(aKP_PLL+bKI_PLL+cKP_PLL/KI_PLL+ d ω in) pi. Wherein a, b, c and d are constant coefficients, and the constant coefficients are determined for a determined control system.
As a preferred embodiment, the control of harmonic filtering can also be achieved by adding an enable switch. Specifically, in the block diagram of fig. 4, Gain _1 and Gain _2 are enable switches for determining whether to turn on/off the filtering algorithm function. When the enable switch states of Gain _1 and Gain _2 are the functions of filtering the first harmonic and filtering the second harmonic, the angular velocity compensation amount P _ out corresponding to the correction axis error Δ θ' of filtering the first harmonic component and the second harmonic component is obtained as P _ out1+ P _ out 2. If the enable switch states of Gain _1 and Gain _2 are the functions of filtering the first harmonic and the second harmonic, the whole axis error filtering function is turned off, and the angular velocity compensation amount P _ out cannot be output. If one of the enable switches is in the state of turning on the filtering algorithm function, and the other enable switch is in the state of turning off the filtering algorithm function, the obtained angular velocity compensation quantity P _ out is only the angular velocity compensation quantity for filtering the first harmonic (the Gain _1 enable switch is in the state of turning on the filtering first harmonic function, and the Gain _2 enable switch is in the state of turning off the filtering second harmonic function), or is only the angular velocity compensation quantity for filtering the second harmonic (the Gain _1 enable switch is in the state of turning off the filtering first harmonic function, and the Gain _2 enable switch is in the state of turning on the filtering second harmonic function).
In the embodiment of filtering only the first harmonic component, the process of extracting the first harmonic component and filtering the first harmonic component in fig. 4 may be directly adopted. Of course, in the embodiment of filtering only the first harmonic component, the control of filtering the first harmonic component may also be implemented by adding an enable switch, and the specific implementation manner is also referred to fig. 4 and will not be repeated herein.
As a preferred embodiment, the filtering processing is performed on the first angular velocity difference Δ ω 2 to obtain a filtered angular velocity Δ ω _ K after at least part of the angular velocity fluctuation is filtered, and the method specifically includes: and extracting partial angular velocity fluctuation K _ out in the first angular velocity difference delta omega 2 by adopting a velocity fluctuation extraction algorithm, and calculating the difference between the first angular velocity difference delta omega 2 and the partial angular velocity fluctuation K _ out, wherein the difference is determined as the filtering angular velocity delta omega _ K.
In some other preferred embodiments, a speed fluctuation extraction algorithm is used to extract a partial angular velocity fluctuation in the first angular velocity difference, and a difference between the first angular velocity difference and the partial angular velocity fluctuation is calculated, where the difference is determined as a filtered angular velocity, and the method specifically includes: and adopting a speed fluctuation extraction algorithm to extract at least a first harmonic component in the first angular speed difference value as part of angular speed fluctuation, calculating the difference value between the first angular speed difference value and the first harmonic component, and determining the difference value as a filtering angular speed for filtering at least the first harmonic component. As a more preferable embodiment, the extracting a partial angular velocity fluctuation from the first angular velocity difference by using a velocity fluctuation extraction algorithm, and calculating a difference between the first angular velocity difference and the partial angular velocity fluctuation, where the difference is determined as a filtered angular velocity specifically includes: and extracting a first harmonic component and a second harmonic component in the first angular velocity difference value by adopting a velocity fluctuation extraction algorithm, taking the sum of the first harmonic component and the second harmonic component as part of angular velocity fluctuation, calculating the difference value between the first angular velocity difference value and the sum of the first harmonic component and the second harmonic component, and determining the difference value as the filtering angular velocity after the first harmonic component and the second harmonic component are filtered. Most of fluctuation components in the first angular velocity difference value can be filtered out by filtering out the first harmonic component in the first angular velocity difference value or filtering out the first harmonic component and the second harmonic component in the first angular velocity difference value, the calculated amount is moderate, and the filtering speed is high.
Fig. 5 is a logic block diagram showing a specific example of the speed fluctuation extraction algorithm in fig. 3, specifically, a logic block diagram showing a specific example of extracting the first harmonic component and the second harmonic component from the first angular velocity difference value to form a partial angular velocity fluctuation. Referring to fig. 5, this specific example obtains a partial angular velocity fluctuation containing a first harmonic component and a second harmonic component by the following method:
firstly, a Fourier series expansion is carried out on the first angular velocity difference delta omega 2 to obtain the first angular velocity difference delta omega 2 relative to the mechanical angle thetamIs used for the functional expression of (1). This process can be implemented using existing technology and is not described in detail here.
Then, the first harmonic component and the second harmonic component are extracted from the functional expression, respectively.
Specifically, as shown in FIG. 5, the functional expression is related to cos θm1After multiplication, pass through a low-pass filter
Figure GDA0003067163180000121
Filtering to obtain d-axis component omega of first harmonic before inverse transformationd1Then, carrying out inverse Fourier transform to obtain the d-axis component of the inverse transformed first harmonic; multiplying the functional expression by-sin θm1After multiplication, pass through a low-pass filter
Figure GDA0003067163180000122
Filtering to obtain q-axis component omega of the first harmonic before inverse transformationq1Then carrying out inverse Fourier transform to obtain the q-axis component of the inverse transformed first harmonic; then, the d-axis component and the q-axis component of the inverse-transformed first harmonic are added to obtain a first harmonic component K _ out1 in the first angular velocity difference value. Similarly, the functional expression is related to cos θm2After multiplication, pass through a low-pass filter
Figure GDA0003067163180000123
Filtering to obtain d-axis component omega of second harmonic before inverse transformationd2Then, carrying out inverse Fourier transform to obtain the d-axis component of the second harmonic after inverse transform; multiplying the functional expression by-sin θm2After multiplication, pass through a low-pass filter
Figure GDA0003067163180000124
Filtering to obtain q-axis component omega of second harmonic before inverse transformationq2Then carrying out inverse Fourier transform to obtain q-axis components of the second harmonic after inverse transform; then, the d-axis component and the q-axis component of the inverse-transformed second harmonic are added to obtain a second harmonic component K _ out2 in the first angular velocity difference value. Finally, the first harmonic component K _ out1 is added to the second harmonic component K _ out2, and the resulting sum forms part of the angular velocity fluctuation K _ out. Wherein, thetam1Mechanical angle of first harmonic, theta, in a functional expression developed as a Fourier seriesm2Mechanical angle of the second harmonic in a functional expression developed as a Fourier series, and θm2=2θm1,T_PD_filterIs the time constant of the low pass filter.
After obtaining a part of angular velocity fluctuation K _ out including a first harmonic component and a second harmonic component, calculating a difference between the first angular velocity difference Δ ω 2 and the part of angular velocity fluctuation K _ out as a filtered angular velocity Δ ω _ K, where the filtered angular velocity Δ ω _ K is the filtered angular velocity after the first harmonic component and the second harmonic component are filtered out.
As a preferred embodiment, the control of the harmonic extraction can also be achieved by adding an enable switch. Specifically, in the block diagram of fig. 5, Gain _1 and Gain _2 are enable switches for determining whether to turn on/off the extraction algorithm function. Under the condition that the enabling switch states of the Gain _1 and the Gain _2 are on, the functions of extracting the first harmonic and extracting the second harmonic are obtained, and partial angular velocity fluctuation formed by the first harmonic component and the second harmonic component is obtained: k _ out is K _ out1+ K _ out 2. If the enable switch states of Gain _1 and Gain _2 are the functions of extracting the first harmonic and extracting the second harmonic, the whole speed fluctuation extraction algorithm function is turned off, and part of the angular speed fluctuation is 0. If one of the enable switches is in the state of opening the extraction algorithm function, and the other enable switch is in the state of closing the extraction algorithm function, the obtained part of the angular speed fluctuation is only a first harmonic component in the first angular speed difference (the state of the Gain _1 enable switch is in the state of opening the extraction first harmonic function, and the state of the Gain _2 enable switch is in the state of closing the extraction second harmonic function) or only a second harmonic component in the first angular speed difference (the state of the Gain _1 enable switch is in the state of closing the extraction first harmonic function, and the state of the Gain _2 enable switch is in the state of opening the extraction second harmonic function).
In the embodiment of extracting only the first harmonic component, the process of extracting the first harmonic component in fig. 5 may be directly employed; of course, the control of the first harmonic extraction may also be implemented by adding an enable switch, and the specific implementation manner is also shown in fig. 5, which is not repeated herein.
Fig. 6 is a logic diagram showing a specific example of the torque compensation algorithm in fig. 3, specifically, a logic diagram showing a specific example of obtaining torque compensation amounts corresponding to the first harmonic component and the second harmonic component in the first angular velocity difference. Referring to fig. 6, this embodiment obtains the torque compensation amounts corresponding to the first harmonic component and the second harmonic component in the first angular velocity difference by using the following method:
firstly, a Fourier series expansion is carried out on the first angular velocity difference delta omega 2 to obtain the first angular velocity difference delta omega 2 relative to the mechanical angle thetamIs used for the functional expression of (1). This process can be implemented using existing technology and is not described in detail here.
Then, the d-axis correlation quantity and the q-axis correlation quantity of the first harmonic and the d-axis correlation quantity and the q-axis correlation quantity of the second harmonic are obtained from the functional expression. Specifically, the function expressions are respectively related to cos θm1And-sin θm1Multiplying to obtain d-axis correlation quantity and q-axis correlation quantity of the first harmonic in the first angular speed difference delta omega 2; respectively connecting the function expressions with cos thetam2And-sin θm2Multiplying to obtain the d-axis of the second harmonic in the first angular velocity difference delta omega 2The correlation quantity and the q-axis correlation quantity. Thetam1And thetam2The same as above.
Then, the d-axis related quantity and the q-axis related quantity of the first harmonic and the d-axis related quantity and the q-axis related quantity of the second harmonic are converted into d-axis torque and q-axis torque respectively.
Specifically to this example, as a preferred implementation, two steps of conversion to torque are employed: first, use the 1/T integratorIS is converted, TIRespectively converting d-axis correlation quantity and q-axis correlation quantity of the first harmonic and d-axis correlation quantity and q-axis correlation quantity of the second harmonic into d-axis initial torque delta tau 'of the first harmonic as a time constant of an integrator'd1And primary harmonic q-axis initial moment delta tau'q1D-axis initial moment delta tau 'of second harmonic'd2And q-axis initial moment delta tau 'of the second harmonic'q2. And then, respectively carrying out proportion adjustment on the d-axis initial moment and the q-axis initial moment, and determining the result after the proportion adjustment as the required d-axis moment and the q-axis moment. Specifically, the d-axis number f (ω)d1) D-axis initial moment delta tau 'for first harmonic'd1The proportion is adjusted to obtain the d-axis moment delta tau of the first harmonicd1. d number of axes f (ω)d1) From the d-axis component ω of the first harmonicd1And d-axis initial moment delta tau 'of the first harmonic'd1And (4) determining. Wherein the d-axis component ω of the first harmonicd1The d-axis correlation quantity of the first harmonic is determined, and specifically, the d-axis correlation quantity of the first harmonic is obtained after being filtered by a low-pass filter (see the description of fig. 5). According to q-axis number f (ω)q1) Q-axis initial moment delta tau 'to first harmonic'q1The proportion is adjusted to obtain the q-axis moment delta tau of the first harmonicq1. q axial number f (ω)q1) Q-axis component omega from the first harmonicq1And primary harmonic q-axis initial moment delta tau'q1And (4) determining. Wherein the q-axis component ω of the first harmonicq1Is determined according to the q-axis correlation quantity of the first harmonic, and specifically is obtained after the q-axis correlation quantity of the first harmonic is filtered by a low-pass filter (see the description of fig. 5). According to d-axis number f (ω)d2) D-axis initial moment delta tau 'for second harmonic'd2To make a ratioAdjusting to obtain d-axis torque delta tau of second harmonicd2. d number of axes f (ω)d2) From the d-axis component ω of the second harmonicd2And d-axis initial moment delta tau 'of second harmonic'd2And (4) determining. Wherein the d-axis component ω of the second harmonicd2The d-axis correlation quantity of the second harmonic is determined, and specifically obtained after the d-axis correlation quantity of the second harmonic is filtered by a low-pass filter (see the description of fig. 5). According to q-axis number f (ω)q2) Q-axis initial moment delta tau 'to second harmonic'q2The proportion is adjusted to obtain the q-axis moment delta tau of the second harmonicq2. q axial number f (ω)q2) From the q-component ω of the second harmonicq2And q-axis initial moment delta tau 'of the second harmonic'q2And (4) determining. Wherein the q-axis component ω of the second harmonicq2Is determined according to the q-axis correlation quantity of the second harmonic, and specifically is obtained after the q-axis correlation quantity of the second harmonic is filtered by a low-pass filter (see the description of fig. 5). In other embodiments, the d-axis related quantity and the q-axis related quantity can be directly converted into the corresponding d-axis moment and q-axis moment by only an integrator without proportional adjustment.
And finally, carrying out inverse Fourier transform on the moment to obtain a moment compensation quantity. Specifically, the d-axis moment and the q-axis moment of the first harmonic are respectively related to cos (theta)m1shift-K1) And-sin (theta)m1shift-K1) Summing the multiplication results after Fourier inverse transformation to form a moment compensation amount tau _ out1 corresponding to the first harmonic fluctuation in the first angular velocity difference delta omega 2; the d-axis moment and the q-axis moment of the second harmonic are respectively related to cos (theta)m2shift-K2) And-sin (theta)m2shift-K2) The results of the multiplication and the inverse fourier transform are summed to form a torque compensation amount τ _ out2 corresponding to the second harmonic fluctuation in the first angular velocity difference Δ ω 2. The sum of the two torque compensation amounts forms a torque compensation amount τ _ out corresponding to the first harmonic component and the second harmonic component τ _ out1+ τ _ out 2. Wherein, thetashift-K1And thetashift-K2A phase compensation angle of the first harmonic and a phase compensation angle of the second harmonic, respectively, the angle number of the two phase compensation angles being dependent on a given angular velocityAngular velocity phase determination in the command. The torque compensation amount is obtained through a phase compensation mode, and the output torque after compensation is obtained based on the torque compensation amount, so that the torque phase can deviate and deviate towards the load torque of the compressor, the difference torque of the motor torque and the load torque is further reduced, and the inhibition of the rotation speed fluctuation of the compressor is realized.
As a preferred embodiment, the control of the moment compensation can also be achieved by adding an enabling switch. Specifically, in the block diagram of fig. 6, Gain _1 and Gain _2 are enable switches for determining whether to turn on/off the torque compensation algorithm function. Under the condition that the enabling switch states of the Gain _1 and the Gain _2 are the first harmonic moment compensation and the second harmonic moment compensation, moment compensation quantities corresponding to the first harmonic component and the second harmonic component are obtained: τ _ out is τ _ out1+ τ _ out 2. If the enabling switch states of the Gain _1 and the Gain _2 are the conditions that the first harmonic moment compensation function and the second harmonic moment compensation function are closed, the whole moment compensation algorithm function is closed, and the moment compensation quantity is 0. If one of the enable switches is in the on-moment compensation algorithm function and the other enable switch is in the off-moment compensation algorithm function, the obtained moment compensation quantity is only the moment compensation quantity corresponding to the first harmonic component in the first angular speed difference (the Gain _1 enable switch is in the on-moment compensation function, and the Gain _2 enable switch is in the off-moment compensation function) or is only the moment compensation quantity corresponding to the second harmonic component in the first angular speed difference (the Gain _1 enable switch is in the off-moment compensation function, and the Gain _2 enable switch is in the on-moment compensation function).
In the embodiment of obtaining the torque compensation amount corresponding to the first harmonic component only, the process of obtaining the torque compensation amount corresponding to the first harmonic component in fig. 6 may be directly adopted; of course, the control of the first harmonic moment compensation can also be realized by adding an enable switch, and the specific implementation manner is also shown in fig. 6, which is not described in additional detail herein.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions.

Claims (10)

1. The method for controlling the rotating speed of the air conditioner compressor is characterized by comprising a process of controlling the compressor according to a real-time angular speed and a process of controlling the compressor according to a moment;
the process of controlling the compressor according to the real-time angular velocity includes:
acquiring a shaft error Delta theta reflecting the deviation of the actual position and the estimated position of the compressor rotor;
filtering the shaft error delta theta to obtain a corrected shaft error delta theta 'and an angular speed compensation quantity P _ out corresponding to the corrected shaft error delta theta' after at least part of shaft error fluctuation is filtered;
compensating the angular velocity compensation quantity P _ out into an output angular velocity delta omega _ PLL of a phase-locked loop regulator in the phase-locked loop for controlling the compressor to obtain compensated angular velocity output quantity delta omega ', and obtaining delta omega' ═ P _ out + delta omega _ PLL;
correcting the real-time angular speed omega 1 for controlling the compressor according to the compensated angular speed output quantity delta omega', and controlling the compressor according to the corrected real-time angular speed omega 1;
the filtering processing is performed on the axis error Δ θ to obtain a corrected axis error Δ θ 'obtained by filtering out at least part of axis error fluctuation and an angular velocity compensation amount P _ out corresponding to the corrected axis error Δ θ', and specifically includes:
performing Fourier series expansion on the axis error delta theta to obtain the axis error delta theta related to the mechanical angle thetamThe functional expression of (a);
the function expressions are respectively related to cos thetamnAnd-sin θmnAfter multiplication, extracting d-axis components and q-axis components of n-th harmonic waves of axis errors delta theta; thetamnMechanical angle for nth harmonic;
filtering d-axis components and q-axis components of partial harmonic waves by using an integrator to obtain a filtering result, and realizing filtering processing on the axis error delta theta;
respectively comparing the result of filtering out the d-axis component of partial harmonic and the result of filtering out the q-axis component of partial harmonic with cos (theta)mnshift-Pn) And-sin (theta)mnshift-Pn) Multiplying and performing inverse Fourier transform to obtain an angular velocity compensation quantity P _ out corresponding to the correction axis error delta theta' of the filtered partial harmonic component; thetashift-PnPhase compensation angle for nth harmonic;
the process of controlling the compressor according to the torque includes:
calculating the difference between the target angular velocity fluctuation amount and the feedback angular velocity amount to obtain a first angular velocity difference value; the feedback angular velocity amount is the sum of the direct current component P _ DC of the angular velocity compensation amount P _ out and the compensated angular velocity output amount Δ ω';
filtering the first angular velocity difference to obtain a filtered angular velocity at least part of which is filtered out of angular velocity fluctuation, and inputting the filtered angular velocity serving as an input quantity into a speed ring regulator in a speed ring for controlling a compressor to obtain an output torque of the speed ring regulator; meanwhile, executing moment compensation based on the first angular velocity difference to obtain moment compensation amount corresponding to part of angular velocity fluctuation in the first angular velocity difference;
compensating the torque compensation amount to the output torque of the speed ring regulator to obtain the compensated output torque;
and controlling the compressor according to the compensated output torque.
2. The method according to claim 1, wherein the filtering the axis error Δ θ to obtain a corrected axis error Δ θ' after filtering at least part of the axis error fluctuation, specifically comprises:
and filtering the axis error delta theta to at least filter a d-axis component and a q-axis component of a first harmonic in the axis error delta theta, so as to filter a first harmonic component of the axis error delta theta and obtain a corrected axis error delta theta' in which the first harmonic component is at least filtered.
3. The method according to claim 2, wherein said filtering said axis error Δ θ to obtain a corrected axis error Δ θ' after filtering out at least a portion of the axis error fluctuation, further comprises: and filtering d-axis components and q-axis components of second harmonics in the axis error delta theta, realizing filtering of first harmonic components and second harmonic components of the axis error delta theta, and obtaining corrected axis error delta theta' of filtering the first harmonic components and the second harmonic components.
4. The method of claim 1, wherein the nth harmonic phase is compensated by an angle θshift-PnAccording to the closed loop gain parameter K of the phase-locked loopP_PLL、KI_PLLAnd the angular speed command ω _ in of the phase-locked loop is determined and satisfies:
θshift-Pn=(aKP_PLL+bKI_PLL+cKP_PLL/KI_PLL+ d ω _ in) × pi, a, b, c, d are constant coefficients.
5. The method according to any one of claims 1 to 4, wherein the filtering the first angular velocity difference value to obtain a filtered angular velocity after at least part of angular velocity fluctuations are filtered out includes:
and extracting part of angular velocity fluctuation in the first angular velocity difference by adopting a velocity fluctuation extraction algorithm, calculating the difference between the first angular velocity difference and the part of angular velocity fluctuation, and determining the difference as the filtering angular velocity.
6. The method according to claim 5, wherein the extracting, by using a velocity fluctuation extraction algorithm, a partial angular velocity fluctuation in the first angular velocity difference, and calculating a difference between the first angular velocity difference and the partial angular velocity fluctuation, the difference being determined as the filtered angular velocity, specifically comprises:
and adopting a speed fluctuation extraction algorithm to extract at least a first harmonic component in the first angular speed difference value, taking the first harmonic component as the partial angular speed fluctuation, calculating the difference value between the first angular speed difference value and the first harmonic component, and determining the difference value as a filtering angular speed for filtering at least the first harmonic component.
7. The method according to claim 6, wherein the extracting the first harmonic component from the first angular velocity difference value by using a velocity fluctuation extraction algorithm specifically comprises:
performing Fourier series expansion on the first angular velocity difference to obtain a first angular velocity difference value related to a mechanical angle thetamThe functional expression of (a);
extracting a d-axis component and a q-axis component of the first harmonic from the function expression respectively;
and adding the d-axis component and the q-axis component of the first harmonic to obtain a first harmonic component in the first angular velocity difference.
8. The method of claim 6, wherein the extracting a portion of angular velocity fluctuations in the first angular velocity difference using a velocity fluctuation extraction algorithm further comprises: extracting a second harmonic component in the first angular velocity difference value by adopting a velocity fluctuation extraction algorithm, and taking the sum of the first harmonic component and the second harmonic component as the partial angular velocity fluctuation;
the calculating a difference between the first angular velocity difference and the partial angular velocity fluctuation, which is determined as the filter angular velocity, further includes: and calculating the difference value between the first angular velocity difference value and the sum of the first harmonic component and the second harmonic component, wherein the difference value is determined as the filtering angular velocity after the first harmonic component and the second harmonic component are filtered.
9. The method according to any one of claims 1 to 4, wherein the performing torque compensation based on the first angular velocity difference to obtain a torque compensation amount corresponding to a part of angular velocity fluctuation in the first angular velocity difference specifically includes:
performing Fourier series expansion on the first angular velocity difference to obtain a first angular velocity difference value related to a mechanical angle thetamThe functional expression of (a);
the function expressions are respectively related to cos thetamnAnd-sin θmnMultiplying to obtain d-axis correlation quantity and q-axis correlation quantity of the n-th harmonic of the first angular velocity difference; thetamnMechanical angle for nth harmonic;
converting the d-axis correlation quantity and the q-axis correlation quantity of the n-th harmonic into d-axis moment and q-axis moment of the n-th harmonic respectively;
respectively connecting the d-axis moment and the q-axis moment of the n-th harmonic with cos (theta)mnshift-Kn) And-sin (theta)mnshift-Kn) Multiplying and performing inverse Fourier transform to obtain moment compensation quantity of the n-th harmonic, and determining the moment compensation quantity as the moment compensation quantity corresponding to part of angular velocity fluctuation in the first angular velocity difference; thetashift-KnA phase compensation angle for the nth harmonic, the phase compensation angle being determined from the angular velocity phase in a given angular velocity command ω in;
the n-th harmonic is a first harmonic and a second harmonic; or, the nth harmonic is a first harmonic.
10. The method according to claim 9, wherein the converting the d-axis correlation quantity and the q-axis correlation quantity of the n-th harmonic into a d-axis moment and a q-axis moment of the n-th harmonic respectively comprises:
respectively converting the d-axis correlation quantity and the q-axis correlation quantity of the n-th harmonic into a d-axis initial moment and a q-axis initial moment of the n-th harmonic by adopting an integrator;
and respectively carrying out proportion adjustment on the d-axis initial moment and the q-axis initial moment of the n-th harmonic, and determining the result after the proportion adjustment as the d-axis moment and the q-axis moment of the n-th harmonic.
CN201811529484.3A 2018-12-13 2018-12-13 Air conditioner compressor rotating speed control method Active CN109724312B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811529484.3A CN109724312B (en) 2018-12-13 2018-12-13 Air conditioner compressor rotating speed control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811529484.3A CN109724312B (en) 2018-12-13 2018-12-13 Air conditioner compressor rotating speed control method

Publications (2)

Publication Number Publication Date
CN109724312A CN109724312A (en) 2019-05-07
CN109724312B true CN109724312B (en) 2021-09-21

Family

ID=66296078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811529484.3A Active CN109724312B (en) 2018-12-13 2018-12-13 Air conditioner compressor rotating speed control method

Country Status (1)

Country Link
CN (1) CN109724312B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112879297B (en) * 2021-01-28 2023-07-18 青岛海尔空调器有限总公司 Control method of air conditioner compressor and air conditioner

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004166408A (en) * 2002-11-13 2004-06-10 Yoichi Hayashi Permanent magnet synchronous motor control method
JP2007166690A (en) * 2005-12-09 2007-06-28 Hitachi Appliances Inc Motor control device
CN101330270A (en) * 2007-06-22 2008-12-24 三洋电机株式会社 Motor control device and compressor
JP2012005199A (en) * 2010-06-15 2012-01-05 Toshiba Corp Motor controller, compressor and heat pump device
CN103967794A (en) * 2013-02-05 2014-08-06 广东美的制冷设备有限公司 Vibration compensation method for single-rotor compressor and controller
CN104038127A (en) * 2013-03-07 2014-09-10 日立空调·家用电器株式会社 Motor control device
CN105515484A (en) * 2016-01-14 2016-04-20 广东美芝制冷设备有限公司 Rotary vibration inhibition method and device of compressor and compressor control system
CN106788071A (en) * 2017-01-06 2017-05-31 南京航空航天大学 A kind of method for improving permanent-magnet synchronous motor rotor position estimated accuracy

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004166408A (en) * 2002-11-13 2004-06-10 Yoichi Hayashi Permanent magnet synchronous motor control method
JP2007166690A (en) * 2005-12-09 2007-06-28 Hitachi Appliances Inc Motor control device
CN101330270A (en) * 2007-06-22 2008-12-24 三洋电机株式会社 Motor control device and compressor
JP2012005199A (en) * 2010-06-15 2012-01-05 Toshiba Corp Motor controller, compressor and heat pump device
CN103967794A (en) * 2013-02-05 2014-08-06 广东美的制冷设备有限公司 Vibration compensation method for single-rotor compressor and controller
CN104038127A (en) * 2013-03-07 2014-09-10 日立空调·家用电器株式会社 Motor control device
CN105515484A (en) * 2016-01-14 2016-04-20 广东美芝制冷设备有限公司 Rotary vibration inhibition method and device of compressor and compressor control system
CN106788071A (en) * 2017-01-06 2017-05-31 南京航空航天大学 A kind of method for improving permanent-magnet synchronous motor rotor position estimated accuracy

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于傅里叶变换的空调压缩机转速波动抑制方法;张国柱;《电器》;20131015(第2013(S1)期);491-496 *
电动汽车空调压缩机永磁电机无传感器控制;冯慧;《上海交通大学硕士学位论文》;20150123;17-18 *

Also Published As

Publication number Publication date
CN109724312A (en) 2019-05-07

Similar Documents

Publication Publication Date Title
CN109724330B (en) Method for suppressing fluctuation of rotation speed of compressor of air conditioner
CN109724297B (en) Method for controlling fluctuation of rotating speed of compressor
CN109698647B (en) Method for inhibiting fluctuation of rotating speed of air conditioner compressor
CN109510553B (en) Method for controlling fluctuation of rotating speed of air conditioner compressor
CN109586643B (en) Method for controlling speed fluctuation of single-rotor compressor
CN109724302B (en) Air conditioner compressor rotating speed control method
CN109751232B (en) Method for inhibiting fluctuation of rotating speed of air conditioner compressor
CN109458336B (en) Method for controlling rotation speed of single-rotor compressor
CN109724312B (en) Air conditioner compressor rotating speed control method
CN109742996B (en) Method for controlling fluctuation of rotating speed of air conditioner compressor
CN109724325B (en) Method for controlling rotation speed of air conditioner compressor
CN109713963B (en) Method for suppressing fluctuation of rotation speed of air conditioner compressor
CN109510554B (en) Method for suppressing fluctuation of rotation speed of air conditioner compressor
CN109724329B (en) Method for suppressing fluctuation of rotation speed of compressor of air conditioner
CN109724308B (en) Compressor rotating speed control method
CN109724327B (en) Method for controlling fluctuation of rotating speed of compressor of air conditioner
CN109724300B (en) Method for controlling rotational speed of compressor
CN109458339B (en) Method for controlling rotating speed of single-rotor compressor
CN109724310B (en) Method for controlling rotating speed of air conditioner compressor
CN109724311B (en) Method for controlling rotating speed of air conditioner compressor
CN109737063B (en) Method for controlling fluctuation of rotating speed of air conditioner compressor
CN109639208B (en) Method for controlling fluctuation of rotating speed of air conditioner compressor
CN109713964B (en) Method for controlling fluctuation of rotation speed of air conditioner compressor
CN109724318B (en) Method for controlling fluctuation of rotating speed of air conditioner compressor
CN109724319B (en) Air conditioner compressor rotating speed control method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20210823

Address after: No.1 Gangcheng South Road, Jiangbei District, Chongqing, 400026

Applicant after: CHONGQING HAIER AIR-CONDITIONER Co.,Ltd.

Applicant after: QINGDAO HAIER AIR CONDITIONER GENERAL Corp.,Ltd.

Applicant after: Haier Zhijia Co.,Ltd.

Address before: 266101 Haier Industrial Park, 1 Haier Road, Laoshan District, Shandong, Qingdao

Applicant before: QINGDAO HAIER AIR CONDITIONER GENERAL Corp.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant