CN105424276B - A kind of method and apparatus for obtaining motor rotary inertia - Google Patents
A kind of method and apparatus for obtaining motor rotary inertia Download PDFInfo
- Publication number
- CN105424276B CN105424276B CN201510946138.5A CN201510946138A CN105424276B CN 105424276 B CN105424276 B CN 105424276B CN 201510946138 A CN201510946138 A CN 201510946138A CN 105424276 B CN105424276 B CN 105424276B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- nth
- msup
- inertia
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 118
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 claims description 61
- 239000011159 matrix material Substances 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/10—Determining the moment of inertia
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Electric Motors In General (AREA)
- Centrifugal Separators (AREA)
Abstract
The embodiment of the present application provides a kind of method and apparatus for obtaining motor rotary inertia, for obtaining the rotary inertia of motor.Methods described includes:Obtain the n-th 1 actual speeds of the motor in the n-th 1 cycles;N is positive integer;Obtain the output of corresponding with motor adjustable model the n-th 1 estimate rotating speed;The n-th 1 presumption rotating speed rotating speed based on the n-th 1 wheel result of calculation outputs that is the adjustable model;According to the n-th 1 actual speed and the n-th 1 presumption rotating speed, the n-th wheel result of calculation is obtained;The n-th wheel result of calculation includes the n-th rotary inertia correlative;Judge whether the n-th 1 rotary inertia correlatives in the n-th rotary inertia correlative and the n-th 1 wheel result of calculation meet preset relation;When the n-th rotary inertia correlative and the n-th 1 rotary inertia correlative meet the preset relation, the actual rotation inertia that the n-th rotary inertia corresponding to the n-th rotary inertia correlative is the motor is determined.
Description
Technical Field
The invention relates to the technical field of electronics, in particular to a method and a device for obtaining the rotational inertia of a motor.
Background
The moment of inertia is an important parameter of the motor. At present, when the motor is controlled, the rotational inertia of the specification parameters provided by a motor manufacturer is generally directly adopted. However, the moment of inertia is not a constant value, and it is related to the actual operating conditions of the motor. For example, motors operate with and without fan blades, and the rotational inertia in both cases is very different. Therefore, the environment for testing the specification parameters is different from the actual application of the motor, which causes a larger or smaller error between the given moment of inertia and the actual moment of inertia in the specification parameters. And the error between the given moment of inertia and the actual moment of inertia can cause the technical problem of inaccurate control of the motor.
Disclosure of Invention
The embodiment of the application provides a method and a device for obtaining the rotational inertia of a motor, which are used for obtaining the rotational inertia of the motor.
In a first aspect, the present application provides a method for obtaining a rotational inertia of a motor, comprising:
acquiring the n-1 actual rotating speed of the motor in the n-1 period; n is a positive integer;
obtaining the n-1 th estimated rotating speed output by the adjustable model corresponding to the motor; the n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th round calculation result;
obtaining the nth round of calculation result according to the n-1 actual rotating speed and the n-1 estimated rotating speed; the nth round of calculation result comprises nth moment of inertia related quantity;
judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not;
and when the relevant quantity of the nth rotational inertia and the relevant quantity of the (n-1) th rotational inertia meet the preset relation, determining the nth rotational inertia corresponding to the relevant quantity of the nth rotational inertia as the actual rotational inertia of the motor.
Optionally, after determining whether the nth relevant rotational inertia amount and the nth-1 relevant rotational inertia amount in the n-1 th round of calculation result satisfy a preset relationship, the method further includes:
and when the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity do not meet the preset relation, inputting the nth calculation result into the adjustable model, so that the adjustable model outputs an nth estimated rotating speed different from the nth-1 estimated rotating speed in an n +1 th period based on the nth calculation result.
Optionally, obtaining an nth calculation result according to the n-1 th actual rotation speed and the n-1 st estimated rotation speed, where the nth calculation result includes:
obtaining the n-th round calculation result according to the following formula:
wherein,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the n-th load torque,Bm0is the initial viscous friction coefficient, J, of the adjustable model0Is the initial moment of inertia, T, of the adjustable modell0Is the initial load torque, Ω, of the adjustable modeliIth actual rotational speed, omega'iFor the i-th estimated rotation speed, Ωn-1Is the actual rotation speed of n-1, omega'n-1Estimating a rotation speed, T, for said n-1eiIs the ith torque, Ten-1Is the n-1 th torque, k1k2k3k4k5k6Not equal to 0, Δ T is the period.
Optionally, when the n-th round calculation result at least includes Jn'、Bmn' and Tln', the nth moment of inertia related quantity is JnIf yes, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the n-1 st round of calculation result meet a preset relation or not, wherein the judging step comprises the following steps of:
judging the n-1 th moment of inertia Jn-1' and Jn' is within a preset range; when J isn-1' and Jn' when the ratio is within the preset range, it indicates that the nth rotational inertia related quantity and the n-1 th rotational inertia related quantity satisfy the preset relationship; wherein, Jn-1' is the n-1 th moment of inertia related quantity.
Optionally, when the n-th calculation result at least includes an'、bn' and cn', the nth moment of inertia related quantity is bnIf yes, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the n-1 st round of calculation result meet a preset relation or not, wherein the judging step comprises the following steps of:
judgment bn-1' and bn' whether the absolute value of the difference is less than a threshold; when the absolute value is smaller than the threshold, the nth moment of inertia related quantity and the (n-1) th moment of inertia related quantity meet the preset relation; wherein, bn-1' is the n-1 th moment of inertia related quantity.
Optionally, the adjustable model outputs an nth estimated rotation speed different from the nth-1 estimated rotation speed in an n +1 th cycle based on the nth calculation result according to the following formula:
wherein omegan' is the n-th estimated rotation speed, TenIs the nth torque.
In another aspect, the present application provides an apparatus for obtaining a rotational inertia of a motor, including:
the actual rotating speed obtaining unit is used for obtaining the n-1 actual rotating speed of the motor in the n-1 period; n is a positive integer;
an estimated rotation speed obtaining unit for obtaining an n-1 th estimated rotation speed of the adjustable model output corresponding to the motor; the n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th round calculation result;
the calculation unit is used for obtaining the calculation result of the nth round according to the n-1 actual rotating speed and the n-1 estimated rotating speed; the nth round of calculation result comprises nth moment of inertia related quantity;
the judging unit is used for judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not;
and the determining unit is used for determining the nth moment of inertia corresponding to the nth moment of inertia relevant quantity as the actual moment of inertia of the motor when the nth moment of inertia relevant quantity and the nth-1 moment of inertia relevant quantity meet the preset relation.
Optionally, the apparatus further comprises:
and the input unit is used for inputting the calculation result of the nth round into the adjustable model after judging whether the correlation quantity of the nth moment of inertia and the correlation quantity of the nth-1 moment of inertia in the calculation result of the nth round-1 meet a preset relation or not, so that the adjustable model outputs an nth estimated rotating speed different from the nth-1 estimated rotating speed in an n +1 th period based on the calculation result of the nth round when the correlation quantity of the nth moment of inertia and the correlation quantity of the nth-1 moment of inertia do not meet the preset relation.
Optionally, the calculating unit is configured to obtain the nth calculation result according to the following formula:
wherein,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the n-th load torque,Bm0is the initial viscous friction coefficient, J, of the adjustable model0Is the initial moment of inertia, T, of the adjustable modell0Is the initial load torque, Ω, of the adjustable modeliIth actual rotational speed, omega'iFor the i-th estimated rotation speed, Ωn-1Is the actual rotation speed of n-1, omega'n-1Estimating a rotation speed, T, for said n-1eiIs the ith torque, Ten-1Is the n-1 th torque, k1k2k3k4k5k6Not equal to 0, Δ T is the period.
Optionally, when the n-th round calculation result at least includes Jn'、Bmn' and Tln', the nth moment of inertia related quantity is Jn' the judging unit is used for judging the n-1 th moment of inertia Jn-1' and Jn' is within a preset range; when J isn-1' and Jn' when the ratio is within the preset range, it indicates that the nth rotational inertia related quantity and the n-1 th rotational inertia related quantity satisfy the preset relationship; wherein, Jn-1' is the n-1 th moment of inertia related quantity.
Optionally, when the n-th calculation result at least includes an'、bn' and cn', the nth moment of inertia related quantity is bn' then, the judging unit is used for judging bn-1' and bn' whether the absolute value of the difference is less than a threshold; when the absolute value is smaller than the threshold, the correlation quantity of the nth moment of inertia and the nth moment of inertia are represented-1 the moment of inertia related quantity satisfies the preset relation; wherein, bn-1' is the n-1 th moment of inertia related quantity.
Optionally, the adjustable model outputs an nth estimated rotation speed different from the nth-1 estimated rotation speed in an n +1 th cycle based on the nth calculation result according to the following formula:
wherein omegan' is the n-th estimated rotation speed, TenIs the nth torque.
One or more technical solutions in the embodiments of the present application have at least one or more of the following technical effects:
in the technical scheme of the embodiment of the application, the n-1 actual rotating speed of the motor in the n-1 period and the n-1 estimated rotating speed output by the adjustable model corresponding to the motor are obtained firstly. The n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th calculation result; n is a positive integer. And then, obtaining an nth calculation result according to the (n-1) th actual rotating speed and the (n-1) th estimated rotating speed, wherein the nth calculation result comprises the nth moment of inertia related quantity. And then, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not, and when the preset relation is met, determining that the nth moment of inertia corresponding to the nth moment of inertia related quantity is the actual moment of inertia of the motor. Therefore, the actual rotational inertia of the motor is obtained through the technical scheme in the embodiment of the application, and the motor can be controlled more accurately according to the actual rotational inertia.
Drawings
FIG. 1 is a flow chart of a method for obtaining a rotational inertia of a motor according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of a model frame for obtaining the rotational inertia of a motor in the embodiment of the present application;
fig. 3 is a schematic structural diagram of a device for obtaining the rotational inertia of the motor in the embodiment of the present application.
Detailed Description
The embodiment of the application provides a method and a device for obtaining the rotational inertia of a motor, which are used for obtaining the rotational inertia of the motor.
The technical scheme provided by the application has the following general idea:
in the technical scheme of the embodiment of the application, the n-1 actual rotating speed of the motor in the n-1 period and the n-1 estimated rotating speed output by the adjustable model corresponding to the motor are obtained firstly. The n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th calculation result; n is a positive integer. And then, obtaining an nth calculation result according to the (n-1) th actual rotating speed and the (n-1) th estimated rotating speed, wherein the nth calculation result comprises the nth moment of inertia related quantity. And then, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not, and when the preset relation is met, determining that the nth moment of inertia corresponding to the nth moment of inertia related quantity is the actual moment of inertia of the motor. Therefore, the actual rotational inertia of the motor is obtained through the technical scheme in the embodiment of the application, and the motor can be controlled more accurately according to the actual rotational inertia.
The technical solutions of the present invention are described in detail below with reference to the drawings and specific embodiments, and it should be understood that the specific features in the embodiments and examples of the present invention are described in detail in the technical solutions of the present application, and are not limited to the technical solutions of the present application, and the technical features in the embodiments and examples of the present application may be combined with each other without conflict.
The term "and/or" herein is merely an association describing an associated object, meaning that three relationships may exist, e.g., a and/or B, may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
In a first aspect, the present application provides a method for obtaining a rotational inertia of a motor, please refer to fig. 1, including the following steps:
s101: and acquiring the n-1 actual rotating speed of the motor in the n-1 period.
S102: and obtaining the n-1 th estimated rotating speed output by the adjustable model corresponding to the motor.
S103: and obtaining the n-th round of calculation results according to the n-1 actual rotating speed and the n-1 estimated rotating speed.
S104: and judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation.
S105: and when the relevant quantity of the nth rotational inertia and the relevant quantity of the (n-1) th rotational inertia meet the preset relation, determining the nth rotational inertia corresponding to the relevant quantity of the nth rotational inertia as the actual rotational inertia of the motor.
In S101, the n-1 actual rotating speed of the motor in the n-1 period is obtained, specifically, the actual rotating speed of the motor in the n-1 period is obtained through the rotating speed sensor detection or the control software. In the embodiment of the present application, the period refers to a period calculated in each round, and the period length Δ T may be consistent with the period of the pulse modulation signal, for example, 1.4 × 10-4s to 2.5X 10-4Any length in s. In addition, n is a positive integer, for example, 3, 4, 10, or the like.
In S102, the n-1 th estimated rotation speed of the adjustable model output corresponding to the motor is obtained. In particular, during operation of the motor, a real model of the motor, i.e. the motor transfer function, is obtained. In the embodiment of the application, the transfer function of the adjustable model is consistent with the actual model form of the motor,and a function with adjustable parameters. Therefore, the adjustable model is corresponding to the motor. The estimated rotation speed is a rotation speed obtained from a function and an input of the adjustable model. For example, a realistic model of an electric machine isWherein J is the actual moment of inertia of the motor, Ω is the actual rotational speed of the output, TeFor the input torque command value, TlFor actual load torque, BmIs the actual viscous friction coefficient. The adjustable model is then in conformity with the actual model form of the motorWherein J 'is the moment of inertia of the adjustable model, omega' is the rotating speed output by the adjustable model, TeFor inputting the torque of the adjustable model, the same as the torque command value of the input motor, Tl' load Torque for Adjustable model, Bm' is the viscous friction coefficient of the adjustable model.
The n-1 th estimated rotation speed is the rotation speed output by the adjustable model based on the previous round of calculation, i.e., the n-1 th round of calculation. Specifically, the adjustable model adjusts parameters such as the moment of inertia of the adjustable model, the load torque of the adjustable model and/or the viscous friction coefficient of the adjustable model based on the n-1 th round calculation result, and then outputs the n-1 th estimated rotational speed based on the input.
In a specific implementation process, when performing the nth calculation to obtain the actual moment of inertia of the motor, S101 may be performed before S102, S102 may be performed before S101, or S101 and S102 may be performed simultaneously, which is not limited in this application.
After the n-1 actual rotating speed and the n-1 estimated rotating speed are obtained, in S103, the n-th round of calculation is carried out based on the n-1 actual rotating speed and the n-1 estimated rotating speed, and then the n-th round of calculation result is obtained. As shown in fig. 2, in the embodiment of the present application, the actual rotation speed Ω obtained from the motor and the estimated rotation speed Ω' output by the adjustable model are input into the calculation module for calculation in each calculation. And the calculation module processes according to a preset algorithm and further outputs the calculation result of the current round.
Specifically, the nth calculation result in the embodiment of the present application includes the nth rotational inertia related quantity. In the embodiment of the present application, the relevant quantity of the nth moment of inertia may be the nth moment of inertia itself, or may be a dependent quantity of the nth moment of inertia, for example, a reciprocal value of the nth moment of inertia or a value 3 times the nth moment of inertia, and the present application is not particularly limited.
After the nth moment of inertia related quantity is obtained, in S104, whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not is judged. In the embodiment of the present application, the preset relationship indicates that the parameters of the adjustable model are consistent with or close to the parameters of the actual model of the motor. Further, when the nth rotational inertia related quantity and the nth-1 rotational inertia related quantity satisfy the preset relationship, the nth rotational inertia corresponding to the nth rotational inertia related quantity is determined as the actual rotational inertia of the motor in S105.
Specifically, since the functional form of the adjustable model is consistent with the form of the motor actual model, and the motor actual model represents the actual state of the motor, when the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity satisfy the preset relationship, the parameter representing the adjustable model is adjusted to be consistent with or close to the parameter of the motor actual model, thereby representing that the adjustable model and the motor actual model at the moment are the same or close to the same. Therefore, the nth moment of inertia corresponding to the nth moment of inertia related quantity of the adjustable model is the same as or close to the actual moment of inertia of the motor, and therefore the nth moment of inertia can be used as the actual moment of inertia of the motor.
Therefore, as can be seen from the above description, when the nth relevant rotational inertia amount and the nth-1 th relevant rotational inertia amount satisfy the preset relationship, the adjustable model is equivalent to the actual model of the motor, and the nth rotational inertia amount of the adjustable model is equivalent to the actual rotational inertia amount of the motor, so that the actual rotational inertia amount of the motor is obtained by taking the nth rotational inertia amount as the actual rotational inertia amount.
In another embodiment of the present application, after S104, the method may further include:
and when the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity do not meet the preset relation, inputting the nth calculation result into the adjustable model, so that the adjustable model outputs an nth estimated rotating speed different from the nth-1 estimated rotating speed in an n +1 th period based on the nth calculation result.
Specifically, when it is determined in S104 that the nth relevant rotational inertia amount and the nth-1 relevant rotational inertia amount do not satisfy the preset relationship, it indicates that the parameter of the adjustable model is not adjusted to be consistent or nearly consistent with the motor at this time, and therefore, the adjustable model still needs to continuously adjust the parameter. Therefore, in the embodiment of the application, when the nth moment of inertia related quantity and the nth-1 th moment of inertia related quantity do not meet the preset relationship, the nth round of calculation result output by the calculation module is further input into the adjustable model, so that the adjustable model adjusts parameters, and the nth estimated rotating speed different from the nth-1 th estimated rotating speed is output in the (n + 1) th period.
In the embodiment of the application, if the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity do not meet the preset relation, inputting the nth round of calculation result into the adjustable model, and then executing S101 to S104 again until the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity meet the preset relation, and determining the actual moment of inertia of the motor.
Next, how to obtain the n-th round calculation result will be described in detail.
In the embodiment of the present application, the n-th round calculation result is obtained according to the following formula (1):
in particular,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the nth load torque.Bm0Initial viscous coefficient of friction for the adjustable model, J0Initial moment of inertia, T, for adjustable modelsl0Is the initial load torque of the adjustable model. B ism0、J0And Tl0And is arbitrarily set by those skilled in the art. Wherein, due to J0Will be denominator, so J is set0Any non-zero real number can be specifically set, and the application is not particularly limited. OmegaiThe ith actual rotating speed, namely the actual rotating speed of the motor in the ith period, omega'iFor the ith estimated rotation speed, the rotation speed output in the ith cycle may be modeled.
At the initial time, a ' (0), B ' (0), c ' (0) is used as the initial value, or B is used as the initial valuem0、J0And Tl0As an initial value, when n is 1, i.e., the 1 st round of calculation is performed, B is set to bem0'=Bm0、J0'=J0And Tl0'=Tl0Substitution intoOr a'0=a'(0),b'0B '(0) and c'0Substitution into c' (0)Obtaining the 0 th estimated rotation speed omega 'of the adjustable model'0. Detecting the 0 th actual rotating speed omega output by the motor0Further, Ω 'is added'0And Ω0Obtained by substituting the formula (1)So n in the embodiment of the present application starts from 1. To obtain b'1And c'1The process is similar.
When n is 2, i.e. the 2 nd calculation is made, a'1b'1And c'1Substitution intoObtaining a 1 st estimated rotation speed omega'1. And detects the 1 st actual rotation speed omega1Obtained by substituting the formula (1)
The subsequent processes are analogized, and are not repeated herein.
In the embodiment of the present application, k1And k2In particular to an' arbitrary value of convergence, k3And k4In particular to bn' arbitrary value of convergence, k5And k6In particular c isnAny value of convergence. At the same time, k1k2k3k4k5k6Not equal to 0. In the concrete implementation process, k1、k2、k3、k4、k5And k6May or may not be identical.
Ωn-1Is the actual rotation speed of n-1, omega'n-1The rotational speed is estimated for the n-1 th rotation. T iseiInputting the torque command values, T, of the adjustable model and the motor for the ith torque, i.e. the ith perioden-1The torque command values of the motor and the adjustable model are input for the (n-1) th torque, i.e., in the (n-1) th cycle. Δ T is the period.
In the embodiment of the application, the n-th round calculation result at least comprises Jn'、Bmn' and Tln', or at least includes an'、bn' and cn'. Therefore, when the nth round of calculation needs to be performed, Ω will ben-1And omega'n-1J is obtained by inputting the calculation module and then the calculation module executes the calculation based on the formula (1)n'、Bmn' and Tln', or an'、bn' and cn'。
In the embodiment of the present application, an electric machineThe actual model is as follows:and then the adjustable model isWherein J is the actual moment of inertia of the motor, Ω is the actual rotational speed of the output, TeAs input torque, BmFor actual viscous coefficient of friction, TlIs the actual load torque. J 'is the moment of inertia of the adjustable model, omega' is the rotation speed of the adjustable model output, TeAs input torque, Bm' viscous coefficient of friction, T, for Adjustable modelln' is the load torque of the adjustable model.
In a specific implementation process, there are various methods for obtaining the nth moment of inertia based on the (n-1) th actual rotational speed and the (n-1) th estimated rotational speed, for example, based on the lyapunov function, the bopft inequality, the least square method, and the like, and the present application is not limited specifically. In the embodiment of the present application, the derivation process of formula (1) will be obtained in detail by taking the boff function as an example.
First, letThe actual model of the motor can be deformed into:
the actual model of the motor can be deformed into:
wherein, in the formula (2)"·" denotes a differential operator. Converting equation (2) to a matrix form, yielding:
wherein a ═ a ], B ═ B ], and C ═ C.
Similarly, the adjustable model can be deformed into:
wherein,
Bm'is the viscous friction coefficient of the adjustable model to be determined, J' is the moment of inertia of the adjustable model to be determined, Tl' is the load torque of the adjustable model to be determined.
Then, equation (4) is subtracted from equation (3) to obtain
Wherein e- Ω'.
Then, add the linear compensator D to obtain the forward linear constant module
Wherein, I is the unit matrix, I (-w) is the input of the forward constant module, y is the output of the linear part and the input of the nonlinear part, and w is the output of the nonlinear part.
As can be seen from equations (5) and (6),
I(-w)=[A-A']Ω'+[B-B']Te+[C-C']. PublicFormula (7)
3 and 3 I 3 is 3 a 3 unit 3 matrix 3, 3 so 3 w 3 ═ 3 A 3- 3 A 3' 3]Ω'-[B-B']Te-[C-C']. Formula (8)
To make the system shown in fig. 2 a progressive hyperstable system requires a strict and positive transfer function of the equivalent forward block formed by equation (6). According to the positive theory of reality, the transfer function must satisfy the medium PA + ATP ═ Q, and PI ═ D, P and Q are arbitrary positive definite symmetric matrices.
Suppose P is selected as [1 ]]I.e. let P be the identity matrix, then we can conclude that D ═ 1],PA+ATP=-Q=[-2a]. WhileTherefore, it is not only easy to useQ is a positive definite matrix.
Further, the Bobo inequality is
γ0For a finite real number, T represents the time difference between the time of calculation and the initial time, and is specifically equal to n Δ T. From the above equation (6) and equation (8), equation (9) can be converted into:
as can be seen, it is desirable to make the final η (0, t) ≧ gamma0 2is established, needs eta1(0,t)≥-γ1 2,η2(0,t)≥-γ2 2,η3(0,t)≥-γ3 2,γ1、γ2、γ3Are finite real numbers.
Therefore, suppose
ThenIs composed of
Further, suppose thatThen when k is1When the pressure is higher than 0, the pressure is higher,
in addition, toThe two sides are derived to obtain the total,
here, take G2(t,τ)=k2(omega-omega ') omega', it is possible to make
Further, combining equation (11) and equation (13), it can be seen that,
η1(0,t)=η11(0,t)+η12(0,t)≥-γ11 2-γ12 2≥-γ1 2. Formula (14)
Therefore, the assumption of the above equation (10) holds, and thus
Converting d tau in equation (15) to dt to obtain the standard form of equation (15)
similarly, we can assume that b' deduces η2(0,t)≥-γ2 2let c' deduce η3(0,t)≥-γ3 2And then confirm that hypothetical b 'and c' hold, get
And
combining the formula (16), the formula (17) and the formula (18), the formula (17) can be obtained
Finally, discretizing the formula (19) to obtain
In the embodiment of the present application, the nth calculation result has two possibilities. First, at least Jn'、Bmn' and Tln' as the nth round of calculation. In a specific implementation process, the nth calculation result further may include an'、bn' and/or cn' etc., the present application is not particularly limited. Second, at least an'、bn' and cn' as the nth round of calculation. In a specific implementation process, the nth calculation result further includes Jn'、Bmn' and/or Tln'。
Because the nth calculation result has two possibilities, there are two specific implementation manners for determining whether the nth rotational inertia related quantity and the nth-1 rotational inertia related quantity satisfy the preset relationship in S104.
The first method comprises the following steps:
when the n-th calculation result at least includes Jn'、Bmn' and Tln', and the nth moment of inertia related quantity is JnWhen, S104 specifically includes:
judging the n-1 th moment of inertia Jn-1' and JnWhether the ratio of' is within a preset range.
When the relevant quantity of the nth moment of inertia is the nth moment of inertia Jn' by itself, the n-1 st moment of inertia J in the n-1 st round of calculation results is obtainedn-1'. In a first implementation, Jn-1' is the n-1 th moment of inertia related quantity. Then further obtain
In a first implementation manner, whether the nth rotational inertia related quantity and the nth-1 rotational inertia related quantity meet a preset relation is judged, specifically, whether the preset relation is judgedWhether it is within a preset range. In the embodiment of the present application, the predetermined range is a small range around 1, for example [0.995, 1.005 ]]Or [0.998, 1.002 ]]And the like. In a specific implementation process, the smaller the preset range is, the closer the nth moment of inertia, which is finally used as the actual moment of inertia, is to the actual moment of inertia. One of ordinary skill in the art to which this application pertains may rely on practiceThe present application is not particularly limited, as long as it is set.
In particular, whenWhen in a preset range, representsEqual to 1 or close to 1. WhileEqual to 1 or close to 1, indicating Ω'n-1Equal to or close to Ωn-1And further shows that the parameter of the adjustable model is equal to or close to the actual model of the motor at the moment, so that J can be adjustedn' is considered as the actual moment of inertia J of the machine. Therefore, whenAnd when the moment of inertia is within the preset range, the n-th moment of inertia related quantity and the n-1-th moment of inertia related quantity meet the preset relation.
Further, in the first implementation, since the nth moment of inertia related quantity is the nth moment of inertia itself, J is directly determined in S105n' is the actual moment of inertia.
And the second method comprises the following steps:
when the n-th calculation result at least includes an'、bn' and cn', and the nth moment of inertia related quantity is bnWhen, S104 specifically includes:
judgment bn-1' and bn' whether the absolute value of the difference is less than a threshold.
When the nth moment of inertia is related to quantity bn' then, the n-1 st moment of inertia b in the n-1 st round of calculation results is obtainedn-1'. In a second implementation, bn-1' is the n-1 th moment of inertia related quantity. Then further obtain | b'n-b'n-1L or | b'n-1-b'n|。
In a second implementation manner, it is determined whether the nth rotational inertia related quantity and the nth-1 rotational inertia related quantity meet a preset relationship, specifically, it is determined that | b'n-b'n-1L or | b'n-1-b'nIf | is less than the threshold. In the embodiment of the present application, the threshold is 0 or a number close to 0, such as 0, 0.05, or 0.1. In a specific implementation process, the smaller the threshold value is, the closer to 0 is, and finally, the closer to the actual moment of inertia is the nth moment of inertia which is the actual moment of inertia. Those of ordinary skill in the art to which the present application pertains may set the setting according to the practice, and the present application is not particularly limited.
Specifically, | b'n-b'n-1L or | b'n-1-b'nIf | is less than the threshold value, it represents | b'n-b'n-1L or | b'n-1-b'n| is equal to 0 or close to 0. And | b'n-b'n-1L or | b'n-1-b'nL is equal to 0 or close to 0, indicating Ω'n-1Equal to or close to Ωn-1And further shows that the parameter of the adjustable model is equal to or close to the actual model of the motor at the moment, so that J can be adjustedn' is considered as the actual moment of inertia J of the machine. Therefore, when | b'n-b'n-1L or | b'n-1-b'nAnd if the | is smaller than the threshold, the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity meet the preset relation.
Further, due toTherefore, in the second implementation manner, the nth moment of inertia corresponding to the nth moment of inertia related quantity determined in the step S105 is
Further, the adjustable model isSo that the n-1 th output is expressed by the formula (20)The rotation speed is estimated.
Wherein, Jn-1' calculated moment of inertia of n-1, T, for the n-1 th wheelln-1' is the n-1 th load torque, Bmn-1' is the n-1 th viscous friction coefficient, Ten-1Is the (n-1) th torque.
And the adjustable model outputs the nth estimated rotation speed according to the formula (21).
Wherein omegan' is the n-th estimated rotation speed, TenThe torque values of the motor and the adjustable model are input for the nth torque, i.e. the nth period.
Further, in a specific implementation, if the estimated rotation speed is obtained by a computer, the computer specifically performs calculation by using a formula after dispersion. Obtaining Ω, for example, using the formula discretized by formula (21)n', is specifically
Based on the same inventive concept as that of obtaining the rotational inertia of the motor in the foregoing embodiment, a second aspect of the present application further provides an apparatus for obtaining the rotational inertia of the motor, as shown in fig. 3, including:
an actual rotation speed obtaining unit 301, configured to obtain an n-1 th actual rotation speed of the motor in an n-1 th cycle; n is a positive integer;
an estimated rotation speed obtaining unit 302 for obtaining an n-1 th estimated rotation speed of the adjustable model output corresponding to the motor; the n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th calculation result;
a calculating unit 303, configured to obtain an nth calculation result according to the n-1 th actual rotation speed and the n-1 th estimated rotation speed; the nth round of calculation result comprises the nth moment of inertia related quantity;
a determining unit 304, configured to determine whether the nth rotational inertia related quantity and the (o-1) th rotational inertia related quantity in the n-1 th round of calculation result satisfy a preset relationship;
the determining unit 305 is configured to determine, when the nth relevant amount of rotational inertia and the nth-1 th relevant amount of rotational inertia satisfy a preset relationship, that the nth corresponding to the nth relevant amount of rotational inertia is an actual rotational inertia of the motor.
Further, the device in this application embodiment still includes:
and the input unit is used for inputting the nth round of calculation result into the adjustable model when the nth moment of inertia relevant quantity and the nth-1 moment of inertia relevant quantity do not meet the preset relation after judging whether the nth moment of inertia relevant quantity and the nth-1 moment of inertia relevant quantity in the nth round of calculation result meet the preset relation or not, so that the adjustable model outputs the nth estimated rotating speed different from the nth-1 estimated rotating speed in the (n + 1) th period based on the nth round of calculation result.
Specifically, the calculation unit 303 is configured to obtain the n-th calculation result according to the following formula:
in particular,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the n-th load torque,Bm0initial viscous coefficient of friction for adjustable model,J0Initial moment of inertia, T, for adjustable modelsl0Initial load torque, Ω, for the adjustable modeliIth actual rotational speed, omega'iFor the i-th estimated rotation speed, Ωn-1Is the actual rotation speed of n-1, omega'n-1Estimate the rotational speed, T, for the n-1 theiIs the ith torque, Ten-1Is the n-1 th torque, k1k2k3k4k5k6Not equal to 0, Δ T is the period.
When the n-th calculation result at least includes Jn'、Bmn' and Tln', n-th moment of inertia related quantity is Jn' time, the judgment unit 304 is used for judging the n-1 th moment of inertia Jn-1' and Jn' is within a preset range; when J isn-1' and JnWhen the ratio of' is in a preset range, the nth moment of inertia related quantity and the (n-1) th moment of inertia related quantity meet a preset relation; wherein, Jn-1' is the n-1 th moment of inertia related quantity.
When the n-th calculation result at least includes an'、bn' and cn', n-th moment of inertia related quantity is bn' time, the judging unit 304 is for judging bn-1' and bn' whether the absolute value of the difference is less than a threshold; when the absolute value is smaller than the threshold value, the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity meet a preset relation; wherein, bn-1' is the n-1 th moment of inertia related quantity.
Specifically, the adjustable model outputs an nth estimated rotation speed different from an nth-1 estimated rotation speed in an n +1 th period based on the nth calculation result according to the following formula:
wherein omegan' is the n-th estimated rotation speed, TenIs the nth torque.
Various modifications and specific examples of the method for obtaining the rotational inertia of the motor in the embodiments of fig. 1 to 2 are also applicable to the apparatus for obtaining the rotational inertia of the motor in the embodiment, and a person skilled in the art can clearly understand the implementation method of the apparatus for obtaining the rotational inertia of the motor in the embodiment from the foregoing detailed description of the method for obtaining the rotational inertia of the motor, so that the detailed description is omitted here for the sake of brevity of the description.
One or more technical solutions in the embodiments of the present application have at least one or more of the following technical effects:
in the technical scheme of the embodiment of the application, the n-1 actual rotating speed of the motor in the n-1 period and the n-1 estimated rotating speed output by the adjustable model corresponding to the motor are obtained firstly. The n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th calculation result; n is a positive integer. And then, obtaining an nth calculation result according to the (n-1) th actual rotating speed and the (n-1) th estimated rotating speed, wherein the nth calculation result comprises the nth moment of inertia related quantity. And then, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not, and when the preset relation is met, determining that the nth moment of inertia corresponding to the nth moment of inertia related quantity is the actual moment of inertia of the motor. Therefore, the actual rotational inertia of the motor is obtained through the technical scheme in the embodiment of the application, and the motor can be controlled more accurately according to the actual rotational inertia.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (12)
1. A method of deriving a rotational inertia of an electric machine, comprising:
acquiring the n-1 actual rotating speed of the motor in the n-1 period; n is a positive integer;
obtaining the n-1 th estimated rotating speed output by the adjustable model corresponding to the motor; the n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th round calculation result;
obtaining the nth round of calculation result according to the n-1 actual rotating speed and the n-1 estimated rotating speed; the nth round of calculation result comprises nth moment of inertia related quantity;
judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not;
and when the relevant quantity of the nth rotational inertia and the relevant quantity of the (n-1) th rotational inertia meet the preset relation, determining the nth rotational inertia corresponding to the relevant quantity of the nth rotational inertia as the actual rotational inertia of the motor.
2. The method according to claim 1, wherein after determining whether the nth moment of inertia related quantity and the (n-1) th moment of inertia related quantity in the n-1 th round of calculation result satisfy a preset relationship, the method further comprises:
and when the nth moment of inertia related quantity and the nth-1 moment of inertia related quantity do not meet the preset relation, inputting the nth calculation result into the adjustable model, so that the adjustable model outputs an nth estimated rotating speed different from the nth-1 estimated rotating speed in an n +1 th period based on the nth calculation result.
3. The method of claim 2, wherein obtaining the nth round of calculation results based on the nth-1 actual rotational speed and the nth-1 estimated rotational speed comprises:
obtaining the n-th round calculation result according to the following formula:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <msup> <mi>a</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>a</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mi>&Delta;</mi> <mi>T</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>b</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>b</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>i</mi> </mrow> </msub> <mi>&Delta;</mi> <mi>T</mi> <mo>+</mo> <msub> <mi>k</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>c</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>c</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>5</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>&Delta;</mi> <mi>T</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>6</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
wherein,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the n-th load torque,Bm0is the initial viscous friction coefficient, J, of the adjustable model0Is the initial moment of inertia, T, of the adjustable modell0Is the initial load torque, Ω, of the adjustable modeliIth actual rotational speed, omega'iFor the i-th estimated rotation speed, Ωn-1Is the actual rotation speed of n-1, omega'n-1Estimating a rotation speed, T, for said n-1eiIs the ith torque, Ten-1Is the n-1 th torque, k1k2k3k4k5k6Not equal to 0, and delta T is a period; k is a radical of1And k2To make an' arbitrary value of convergence, k3And k4To make bn' arbitrary value of convergence, k5And k6To make cnAn arbitrary value of convergence, and k1k2k3k4k5k6≠0。
4. The method of claim 3, wherein when the nth round of computation results includes at least Jn'、Bmn' and Tln', the nth moment of inertia related quantity is JnIf yes, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the n-1 st round of calculation result meet a preset relation or not, wherein the judging step comprises the following steps of:
judging the n-1 th moment of inertia Jn-1' and Jn' is within a preset range; when J isn-1' and Jn' when the ratio is within the preset range, it indicates that the nth rotational inertia related quantity and the n-1 th rotational inertia related quantity satisfy the preset relationship; wherein, Jn-1' is the n-1 th moment of inertia related quantity.
5. The method of claim 3, wherein when the nth calculation includes at least an'、bn' and cn', the nth moment of inertia related quantity is bnIf yes, judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the n-1 st round of calculation result meet a preset relation or not, wherein the judging step comprises the following steps of:
judgment bn-1' and bn' whether the absolute value of the difference is less than a threshold; when the absolute value is smaller than the threshold, the nth moment of inertia related quantity and the (n-1) th moment of inertia related quantity meet the preset relation; wherein, bn-1' is the n-1 th moment of inertia related quantity.
6. The method of claim 3, wherein the adjustable model outputs an nth estimated speed different from the nth-1 th estimated speed in an n +1 th cycle based on the nth round calculation result according to the following formula:
<mrow> <msup> <msub> <mi>J</mi> <mi>n</mi> </msub> <mo>&prime;</mo> </msup> <mfrac> <mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>&Delta;</mi> <mi>T</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>n</mi> </mrow> </msub> <mo>-</mo> <msup> <msub> <mi>T</mi> <mi>ln</mi> </msub> <mo>&prime;</mo> </msup> <mo>-</mo> <msup> <msub> <mi>B</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> </mrow>
wherein omegan' is the n-th estimated rotation speed, TenIs the nth torque.
7. An apparatus for obtaining rotational inertia of a motor, comprising:
the actual rotating speed obtaining unit is used for obtaining the n-1 actual rotating speed of the motor in the n-1 period; n is a positive integer;
an estimated rotation speed obtaining unit for obtaining an n-1 th estimated rotation speed of the adjustable model output corresponding to the motor; the n-1 th estimated rotating speed is the rotating speed output by the adjustable model based on the n-1 th round calculation result;
the calculation unit is used for obtaining the calculation result of the nth round according to the n-1 actual rotating speed and the n-1 estimated rotating speed; the nth round of calculation result comprises nth moment of inertia related quantity;
the judging unit is used for judging whether the nth moment of inertia related quantity and the nth-1 st moment of inertia related quantity in the calculation result of the nth-1 st round meet a preset relation or not;
and the determining unit is used for determining the nth moment of inertia corresponding to the nth moment of inertia relevant quantity as the actual moment of inertia of the motor when the nth moment of inertia relevant quantity and the nth-1 moment of inertia relevant quantity meet the preset relation.
8. The apparatus of claim 7, wherein the apparatus further comprises:
and the input unit is used for inputting the calculation result of the nth round into the adjustable model after judging whether the correlation quantity of the nth moment of inertia and the correlation quantity of the nth-1 moment of inertia in the calculation result of the nth round-1 meet a preset relation or not, so that the adjustable model outputs an nth estimated rotating speed different from the nth-1 estimated rotating speed in an n +1 th period based on the calculation result of the nth round when the correlation quantity of the nth moment of inertia and the correlation quantity of the nth-1 moment of inertia do not meet the preset relation.
9. The apparatus of claim 8, wherein the computing unit is configured to obtain the nth round of computation result according to the following formula:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <msup> <mi>a</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>a</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mi>&Delta;</mi> <mi>T</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>b</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>b</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>i</mi> </mrow> </msub> <mi>&Delta;</mi> <mi>T</mi> <mo>+</mo> <msub> <mi>k</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>c</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>=</mo> <msup> <mi>c</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>k</mi> <mn>5</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>&Delta;</mi> <mi>T</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>6</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&Omega;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
wherein,Bmn' is the n-th viscous friction coefficient, Jn' is the n-th moment of inertia, Tln' is the n-th load torque,Bm0is the initial viscous friction coefficient, J, of the adjustable model0Is the initial moment of inertia, T, of the adjustable modell0Is the initial load torque, Ω, of the adjustable modeliIth actual rotational speed, omega'iFor the i-th estimated rotation speed, Ωn-1Is the actual rotation speed of n-1, omega'n-1Estimating the rotation speed for the n-1 th,TeiIs the ith torque, Ten-1Is the n-1 th torque, k1k2k3k4k5k6Not equal to 0, and delta T is a period; k is a radical of1And k2To make an' arbitrary value of convergence, k3And k4To make bn' arbitrary value of convergence, k5And k6To make cnAn arbitrary value of convergence, and k1k2k3k4k5k6≠0。
10. The apparatus of claim 9, wherein when the nth round of computation results includes at least Jn'、Bmn' and Tln', the nth moment of inertia related quantity is Jn' the judging unit is used for judging the n-1 th moment of inertia Jn-1' and Jn' is within a preset range; when J isn-1' and Jn' when the ratio is within the preset range, it indicates that the nth rotational inertia related quantity and the n-1 th rotational inertia related quantity satisfy the preset relationship; wherein, Jn-1' is the n-1 th moment of inertia related quantity.
11. The apparatus of claim 9, wherein when the nth round of computation results includes at least an'、bn' and cn', the nth moment of inertia related quantity is bn' then, the judging unit is used for judging bn-1' and bn' whether the absolute value of the difference is less than a threshold; when the absolute value is smaller than the threshold, the nth moment of inertia related quantity and the (n-1) th moment of inertia related quantity meet the preset relation; wherein, bn-1' is the n-1 th moment of inertia related quantity.
12. The apparatus of claim 9, wherein the adjustable model outputs an nth estimated rotational speed different from the nth-1 th estimated rotational speed in an n +1 th cycle based on the nth round calculation result according to the following formula:
<mrow> <msup> <msub> <mi>J</mi> <mi>n</mi> </msub> <mo>&prime;</mo> </msup> <mfrac> <mrow> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> <mo>-</mo> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>&Delta;</mi> <mi>T</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>e</mi> <mi>n</mi> </mrow> </msub> <mo>-</mo> <msup> <msub> <mi>T</mi> <mi>ln</mi> </msub> <mo>&prime;</mo> </msup> <mo>-</mo> <msup> <msub> <mi>B</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>&prime;</mo> </msup> <msub> <msup> <mi>&Omega;</mi> <mo>&prime;</mo> </msup> <mi>n</mi> </msub> </mrow>
wherein omegan' is the n-th estimated rotation speed, TenIs the nth torque.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510946138.5A CN105424276B (en) | 2015-12-16 | 2015-12-16 | A kind of method and apparatus for obtaining motor rotary inertia |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510946138.5A CN105424276B (en) | 2015-12-16 | 2015-12-16 | A kind of method and apparatus for obtaining motor rotary inertia |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105424276A CN105424276A (en) | 2016-03-23 |
CN105424276B true CN105424276B (en) | 2018-03-23 |
Family
ID=55502632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510946138.5A Active CN105424276B (en) | 2015-12-16 | 2015-12-16 | A kind of method and apparatus for obtaining motor rotary inertia |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105424276B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424277B (en) * | 2015-12-16 | 2018-02-16 | 四川长虹电器股份有限公司 | A kind of method and apparatus for obtaining motor rotary inertia |
DE102017218398B3 (en) * | 2017-10-13 | 2019-03-07 | Robert Bosch Gmbh | Method for free-wheeling compensation of a drive unit, control unit and electric bicycle with this control unit and calibration method |
CN109245634B (en) * | 2018-10-30 | 2021-12-07 | 铜仁职业技术学院 | Control method for variable-rotation-inertia permanent magnet synchronous motor |
CN115655577B (en) * | 2022-12-13 | 2023-03-10 | 中海油田服务股份有限公司 | Method and device for measuring rotor rotational inertia and viscous friction coefficient of pulser |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3703824A (en) * | 1970-05-11 | 1972-11-28 | Dunlop Co Ltd | Apparatus for comparing moments of inertia of golf clubs |
CN104596702A (en) * | 2014-12-23 | 2015-05-06 | 北京首钢股份有限公司 | Measurement method for inherent rotational inertia |
CN104792461A (en) * | 2015-04-23 | 2015-07-22 | 西安交通大学 | Online measurement method for revolve-body high-precision moment of inertia |
CN104949798A (en) * | 2015-06-03 | 2015-09-30 | 东南大学 | Method of measuring rotational inertia of motor |
CN105424277A (en) * | 2015-12-16 | 2016-03-23 | 四川长虹电器股份有限公司 | Method and device for obtaining rotational inertia of motor |
-
2015
- 2015-12-16 CN CN201510946138.5A patent/CN105424276B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3703824A (en) * | 1970-05-11 | 1972-11-28 | Dunlop Co Ltd | Apparatus for comparing moments of inertia of golf clubs |
CN104596702A (en) * | 2014-12-23 | 2015-05-06 | 北京首钢股份有限公司 | Measurement method for inherent rotational inertia |
CN104792461A (en) * | 2015-04-23 | 2015-07-22 | 西安交通大学 | Online measurement method for revolve-body high-precision moment of inertia |
CN104949798A (en) * | 2015-06-03 | 2015-09-30 | 东南大学 | Method of measuring rotational inertia of motor |
CN105424277A (en) * | 2015-12-16 | 2016-03-23 | 四川长虹电器股份有限公司 | Method and device for obtaining rotational inertia of motor |
Also Published As
Publication number | Publication date |
---|---|
CN105424276A (en) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105424276B (en) | A kind of method and apparatus for obtaining motor rotary inertia | |
Zribi et al. | Controlling chaos in the permanent magnet synchronous motor | |
CN109249397B (en) | Six-degree-of-freedom robot dynamics parameter identification method and system | |
ITRM20080014A1 (en) | ROTATION AXIS POSITIONING CONTROL DEVICE. | |
EP3076261A1 (en) | Machinery control device and gain determination method for friction compensation | |
CN108227497B (en) | Control method of networked teleoperation system under condition of considering system performance limitation | |
CN107592048A (en) | The adaptive chaos control method of fractional order brushless DC motor system | |
CN105425587A (en) | Hysteresis nonlinear motor identification and control method | |
JP5723582B2 (en) | Method, program and system for configuring controller in sliding mode control system | |
CN106169897A (en) | A kind of motor speed real-time anti-saturation PID control method and device | |
JP2018072265A5 (en) | ||
CN104678763B (en) | Friction compensation and dynamic surface control method based on least squares support vector machine for electromechanical servo system | |
JP2006333594A (en) | Mechanical characteristic modeling apparatus and method, motor controller, and machine control system | |
CN105424277B (en) | A kind of method and apparatus for obtaining motor rotary inertia | |
CN105450124B (en) | A kind of method and apparatus for obtaining the parameter of electric machine | |
CN105553371B (en) | A kind of method and apparatus for obtaining the parameter of electric machine | |
CN105429544B (en) | A kind of method and apparatus obtaining the parameter of electric machine | |
JP4045860B2 (en) | Power transmission system test apparatus and control method thereof | |
JP4026310B2 (en) | Engine torque estimation method for engine bench system | |
JP6541407B2 (en) | Motor control device and program | |
JP6555143B2 (en) | Vehicle testing equipment | |
CN104038132B (en) | There is time-varying and measure the state observation method of the servomotor postponing output and noise | |
Kogut | Anti-lock braking system modelling and parameters identification | |
CN110208691B (en) | Motor load parameter detection method and device | |
Pichardo-Almarza et al. | Bond graph approach to build reduced order observers in linear time invariant systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |