CN104333285B - Permagnetic synchronous motor standard is without sensing station Servocontrol device and method - Google Patents

Permagnetic synchronous motor standard is without sensing station Servocontrol device and method Download PDF

Info

Publication number
CN104333285B
CN104333285B CN201310602081.8A CN201310602081A CN104333285B CN 104333285 B CN104333285 B CN 104333285B CN 201310602081 A CN201310602081 A CN 201310602081A CN 104333285 B CN104333285 B CN 104333285B
Authority
CN
China
Prior art keywords
subtractor
axis current
given
actual
synchronous motor
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
CN201310602081.8A
Other languages
Chinese (zh)
Other versions
CN104333285A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201310602081.8A priority Critical patent/CN104333285B/en
Publication of CN104333285A publication Critical patent/CN104333285A/en
Application granted granted Critical
Publication of CN104333285B publication Critical patent/CN104333285B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

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

Abstract

The invention discloses a kind of permagnetic synchronous motor standard is without sensing station Servocontrol device, including permagnetic synchronous motor for being provided with HALL sensors etc.;HALL sensors are connected with quasi- sensorless strategy module, and quasi- sensorless strategy module is connected with each other with the 5th subtractor, location/velocity integrated controller, the 3rd subtractor, the 4th subtractor, pi regulator and IPark conversion modules respectively;5th subtractor is connected with location/velocity integrated controller;Location/velocity integrated controller and the 3rd subtractor are connected;Pi regulator is connected with IPark conversion modules with the 3rd subtractor, the 4th subtractor respectively;IPark conversion modules are connected with space vector modulation module, and space vector modulation module is connected with three-phase inverter, and three-phase inverter is connected with permagnetic synchronous motor.

Description

Permagnetic synchronous motor standard is without sensing station Servocontrol device and method
Technical field
The present invention relates to textile machine electromechanical integration technology, more particularly to a kind of suitable textile machine point-to-point motion Permagnetic synchronous motor standard is without sensing station Servocontrol device and method.
Background technology
With the continuous proposition of the requirements such as textile industry energy-conservation, high efficiency, high speed, high accuracy, permanent magnetic Ac servo system Unite have because of its own High Power Factor, high efficiency, high dynamic response, high reliability, exceed loading capability and height output torque The advantages of, being applied in textile machine has become a kind of trend, is one of core technology of following high volume applications, city Field prospect is boundless.But there are some special operation conditions in some textile machines:Will be in the crossing system time when running such as tricot machine The interior band yarn motion for driving sley bar to complete " stopping → motion → stopping ", it is desirable to have high response speed and positioning precision, drives The speed of dynamic system positioning time directly affects the raising of tricot machine speed of production;Sewing machine running electric motor starting stops Very frequently, belong to batch operation, its technical process requires that its main drive gear has in higher location control, hundreds of ms and realizes Startup from zero-speed to highest running speed or the shutdown from highest running speed to zero-speed, and stopping accuracy be less than ± 3 ° or Person is less;The concept of two axis servo driving is proposed currently for the control Chinese scholars of embroidery machine embroidery frame, with a DSP + CPLD controls the mode of two servomotors operation to drive embroidery frame to run, originally traditional position, speed, three ring of electric current Resource that control theory can show in programming realization is nervous so as to reducing servosystem.And embroidery machine embroidery frame is only to position Control is required, and with embroidery machine high efficiency, high-precision development, embroidery frame in order to mair motor coordinated operation, to position Response speed meeting more and more higher that control is required, etc..The characteristics of these textile machines have one jointly, it requires servo system The position response dynamic property of system is high, and tracking error is little, and with stronger capacity of resisting disturbance, and speed controlling is required It is not high, this point be more conform with servosystem make point-to-point (point position) move the characteristics of.And traditional location, speed, three ring control of electric current When the motion of point position is made, speed closed loop does not only play a role the positional servosystem (accompanying drawing 1) of structure processed, also as control Series connection link in structure reduces the dynamic response performance of alliance regulation.For this reason, it may be necessary to look for a kind of new method carry The position response performance of high system.With this simultaneously, textile machine (tricot machine as previously mentioned, embroidery machine etc.) usually may require that In working environments such as high temperature, low temperature, air pollutions and exist, physical sensors are deposited In the limitation for having instead resulted in system application.If can cancel physical sensors can undoubtedly expand systematic difference scope, together When can also improve the reliability and environmental suitability of system.And textile industry is a cost sensitive industry, is dropped as far as possible Low control cost is that all enterprises are pursued.For this purpose, the servosystem that some point position motion occasions are used is improving position sound The technology of sensorless strategy, such as patent documentation 1 (Patent No. 200910090386.9) and patent are needed on the basis of answering performance Document 2 (Patent No. 200510108039.6) proposes a kind of AC servo driver for being not required to current sensor, Ke Yishi The Current Sensorless speed closed loop vector controlled of existing non-salient pole permanent magnet synchronous motor, reduces the hardware cost of governing system, But patent documentation 1 and 2 remains position detection unit, total control cost is still higher, higher to cost requirement at some Weaving applications are simultaneously improper.
The content of the invention
The technical problem to be solved in the present invention is to provide the permagnetic synchronous motor that a kind of response speed puts position motion faster Standard is without sensing station servosystem and its control method.
In order to solve above-mentioned technical problem, the present invention provides a kind of permagnetic synchronous motor standard without sensing station SERVO CONTROL Device, including pi regulator, IPark conversion modules, space vector modulation module, three-phase inverter, permagnetic synchronous motor, The integration control of HALL15, quasi- sensorless strategy module, the 3rd subtractor, the 4th subtractor, the 5th subtractor and location/velocity Device processed;The HALL15 is arranged on permagnetic synchronous motor body, when permagnetic synchronous motor works, exports discrete by HALL15 HALL position signallingsThe quasi- sensorless strategy module receives discrete HALL position signallingsQ shaft voltages give uqU given with d shaft voltagesd, Jing after quasi- sensorless strategy module arithmetic, export actual rotor-position signal θm, actual d Axis current signal id, actual q axis current signal iqWith the loading moment T of permagnetic synchronous motorL;5th subtractor receives Given position signal θm *With rotor-position signal θm, and both are made with difference operation outgoing position error signal θm *m;Institute's rheme Put/speed integrated controller receives position error signal θm *m, actual q axis current signal iqWith loading moment TL, Jing fortune Calculate the given q axis current signals of output permagnetic synchronous motor3rd subtractor receives given q axis current signals With actual q axis current signal iq, and both are made with the error signal that difference operation exports q shaft currents4th subtraction Device receives given d axis current signalsWith actual d axis current signal id, and both are made with difference operation output d shaft currents Error signalThe pi regulator receives the error signal of d shaft currents respectivelyWith the error signal of q shaft currentsAnd Jing computings export the given u of q shaft voltages respectivelyqU given with d shaft voltagesd;The IPark conversion modules receive q axles electricity The given u of pressureqU given with d shaft voltagesd, the component of voltage u under static two phase coordinate system of Jing computings outputαWith component of voltage uβ;Institute State space vector modulation module and receive the component of voltage u under static two phase coordinate systemαWith component of voltage uβInput, Jing computings output Six tunnel control signals of three-phase inverter;The six tunnels control letter of the three-phase inverter receiving space Vector Modulation module output Number, drive permagnetic synchronous motor operation.
As to a kind of improvement of the standard without sensing station Servocontrol device of permagnetic synchronous motor:Location/velocity one The implementation method for changing controller is as follows;
The first step, obtain motor the equation of motion be:
Second step, takes the state variable of system:
Convolution (1), formula (2):
3rd step, orderF=-TL/ J, then the state equation of permanent magnetism synchronous electric machine position servo system be:
4th step, the sliding-mode surface for taking similar PD actuators are designed as:
S=c1x1+c2x2 (5)
5th step, chooses control law:
6th step, the control function for obtaining quadrature axis current is:
It is describedFor give q axis current signals,
The overall-in-one control schema of location/velocity is achieved that according to above formula;
The TeFor the output torque of permagnetic synchronous motor, J is motor inertia, and P is motor number of pole-pairs, ΨaFor permanent magnet magnetic Chain c1And c1For constant.
As a kind of further improvement of the standard without sensing station Servocontrol device of permagnetic synchronous motor:Standard is without sensor The phase current reconstruction module and full micr oprocessorism rotor position information estimation block that control module is connected by mutual signal is constituted;Institute State full micr oprocessorism rotor position information estimation block and receive actual q axis current signal iqWith discrete HALL position signallingsThe actual rotor-position signal θ of the full micr oprocessorism rotor position information estimation block outputm, loading moment TLAnd permanent magnetism Actual speed ω of synchronous motorm;The phase current reconstruction module receives actual speed ωm, the given u of q shaft voltagesqWith d shaft voltages Given ud;The actual d axis current signal i of the phase current reconstruction module outputd, actual q axis current signal iq
The implementation method of the phase current reconstruction module is as follows:
The La、Φa、RaThe respectively inductance of permagnetic synchronous motor, magnetic linkage and resistance;
The implementation method of the full micr oprocessorism rotor position information algorithm estimation block is as follows:
Equally by Electrical Motor volume, the kinematical equation of permagnetic synchronous motor is:
Jdωm/dt+Bmωm+TL=Tε (12)
m/ dt=ωm (13)
The BmFor viscous damping coefficient;
Because the sample frequency of controller is significantly larger than the transformation period of perturbing torque, disturbance load torque is used as one Individual state variable, it may be assumed that it is zero that it is a constant, i.e. disturbance load torque time differential;
If:
dTL/ dt=0 (14)
The dynamics state equation of motor is obtained by (12) and (13):
Wherein,
The input variable of formula is output torque Te, state variable is mechanical Angle Position (for the actual rotor-position letter of output Number θm), mechanical angular velocity (be actual speed ωm) and disturbance load torque (be loading moment TL), output variable is mechanical angle position Put (for the actual rotor-position signal θ of outputm);
Formula (5) can be written as formula (16):
By dynamics state equation (15), an omnidirectional vision can be set up, see formula (17):
Wherein,For estimative state variable,For feedback of status Gain battle array;
Obtained by (15)-(16):
Wherein,For observation error.Its characteristic equation is:
Det [sI- (A-KC)]=s3+(k1+Bm/J)s2+(k2+k1Bm/J)s-k3/ J=0 (19)
Appropriate K is selected, makes (A-KC) to have stable, appropriate eigenvalue,When, with x (t), u (t) andNothing Close;WithIt is unrelated;
According to specified expectation limit α, β, γ, then the desired character multinomial of observer is:
s3-(α+β+γ)s2+ (α β+β γ+γ α) s- α β γ=0 (20)
Obtained by (19) and (20):
Assume Bm=0, and α=β=γ, then formula (21) it is rewritable into:
Formula (17) it is rewritable into:
I.e.:
According to the desired characteristic of system, the position for selecting limit to be located, observer is constructed according to formula, you can observe load Torque TL, permagnetic synchronous motor actual speed ωmThe actual rotor-position signal θ with permagnetic synchronous motormValue.
A kind of permagnetic synchronous motor standard is without sensing station method of servo-controlling:
A, permagnetic synchronous motor are started working, after the motor shaft rotation of permagnetic synchronous motor, by HALL15 to standard without biography Sensor control module exports discrete HALL position signallings
B, quasi- sensorless strategy module receive the discrete HALL position signallings of HALL15And the q of pi regulator Shaft voltage gives uqU given with d shaft voltagesd, and by discrete HALL position signallingsQ shaft voltages give uqGive with d shaft voltages Determine udComputing is carried out, actual rotor-position signal θ after computing, is drawnm, actual d axis current signal id, actual q shaft currents Signal iqWith the loading moment T of permagnetic synchronous motorL;Quasi- sensorless strategy module exports actual rotor to the 5th subtractor Position signalling θm, actual q axis current signal i are exported to location/velocity integrated controllerqWith loading moment TL, subtract to the 3rd The actual q axis current signal i of musical instruments used in a Buddhist or Taoist mass outputq, actual d axis current signal i are exported to the 4th subtractord, become to IPark and change the mold The actual rotor-position signal θ of block outputm
C, by master system, to the position signalling θ that the input of the 5th subtractor is givenm *, then by standard without sensor control Molding block receives actual rotor-position signal θm;5th subtractor is again by given position signalling θm *With actual rotor position Confidence θmCarry out making difference operation, draw position error signal θm *m;5th subtractor is by position error signal θm *mOutput Give location/velocity integrated controller;
D, location/velocity integrated controller receive the position error signal θ of the 5th subtractor outputm *m, standard without biography The q axis current signal i of the reality of sensor control module outputqWith loading moment TL;Jing location/velocity integrated controller computings The given q axis current signals of permagnetic synchronous motor are drawn afterwardsLocation/velocity integrated controller is exported to the 3rd subtractor Given q axis current signals
The given q axis current signals of e, the 3rd subtractor receiving position/speed integrated controller outputAnd it is accurate The q axis current signal i of the reality of sensorless strategy module outputq;3rd subtractor is to given q axis current signalsAnd reality The q axis current signal i on borderqThe error signal of q shaft currents is drawn after carrying out computing3rd subtractor is exported to pi regulator The error signal of q shaft currents
F, the 4th subtractor receive the given d axis current signals that master system is providedAnd standard is without sensor control The q axis current signal i of the reality of molding block outputq;The d axis current signals that 4th subtractor will giveAnd the q of reality Axis current signal iqAfter carrying out making difference operation, the error signal of d shaft currents is drawn4th subtractor is exported to pi regulator The error signal of d shaft currents
G, pi regulator receive the error signal of the q shaft currents of the 3rd subtractor outputAnd the 4th subtractor it is defeated The error signal of the d shaft currents for going outAfter the computing of pi regulator, the given u of q shaft voltages is obtainedqIt is given with d shaft voltages ud;Pi regulator gives u to IPark conversion modules output q shaft voltagesqU given with d shaft voltagesd
H, IPark conversion module is respectively received the given u of q shaft voltages of pi regulator outputqU given with d shaft voltagesd, with And the rotor-position signal θ of the reality of quasi- sensorless strategy module outputm;Draw after the computing of IPark conversion modules Component of voltage u under static two phase coordinate systemαWith component of voltage uβ;IPark conversion modules export quiet to space vector modulation module The only component of voltage u under two phase coordinate systemsαWith component of voltage uβ
I, space vector modulation module receive the component of voltage under static two phase coordinate system of IPark conversion modules output uαWith component of voltage uβ;Six tunnel control signals of three-phase inverter are drawn Jing after space vector modulation module arithmetic;Space vector is adjusted Molding block exports six tunnel control signals of three-phase inverter to three-phase inverter 6;
J, three-phase inverter receive six tunnel control signals of the three-phase inverter of space vector modulation module output, and lead to The six tunnel control signals for crossing three-phase inverter drive the operation of permagnetic synchronous motor.
The invention has the beneficial effects as follows, permagnetic synchronous motor standard of the present invention compares biography without sensing station Servocontrol device System position, speed, the three-loop system control structure (accompanying drawing 1) of electric current, due to having lacked speed closed loop this cascaded structure, make new system The response speed to position command change unite faster;Eliminate the positions such as expensive photoelectric encoder, rotary transformer biography Sensor, and adopt particular algorithm to coordinate the HALL15 estimation rotor position informations of cheap low resolution, eliminate electric current biography Sensor, obtains phase current using the method that control theory is reconstructed, not only increases the environmental suitability and reliability of system, and Reduce the control cost of system.It is adapted to the textile machine applications moved in various environment, high efficiency point position.
Description of the drawings
Below in conjunction with the accompanying drawings the specific embodiment of the present invention is described in further detail.
Fig. 1 is traditional location, speed, electric current tricyclic structure servosystem theory diagram;
Fig. 2 is permagnetic synchronous motor theory diagram of the standard without sensing station Servocontrol device;
Theory diagrams of the Fig. 3 for location/velocity overall-in-one control schema module;
Fig. 4 is defined the theory diagram of sensorless strategy module.
Specific embodiment
Embodiment 1, Fig. 2 gives a kind of permagnetic synchronous motor standard without sensing station Servocontrol device;Adjust including PI Section device 3, IPark conversion modules 4, space vector modulation module 5, three-phase inverter 6, permagnetic synchronous motor 12, HALL15, quasi- nothing Sensor control block 14, the 3rd subtractor, the 4th subtractor, the 5th subtractor and location/velocity integrated controller 13; HALL is HALL sensors;
HALL15 is arranged on 12 body of permagnetic synchronous motor, when the motor shaft of permagnetic synchronous motor 12 rotates, is led to Cross HALL15 and export discrete HALL position signallingsQuasi- sensorless strategy module 14 receives discrete HALL position signallingsQ shaft voltages give uqU given with d shaft voltagesd;Discrete HALL position signallingsQ shaft voltages give uqGive with d shaft voltages Determine udJing after quasi- 14 computing of sensorless strategy module, actual rotor-position signal θ is exportedm, actual d axis current signal id、 Actual q axis current signal iqWith the loading moment T of permagnetic synchronous motor 12L;5th subtractor receives actual rotor-position letter Number θmWith given position signalling θm *, the 5th subtractor is to actual rotor-position signal θmWith given position signalling θm *Carry out Position error signal θ is drawn after making difference operationm *m;Location/velocity integrated controller 13 receives position error signal θm *m、 Actual q axis current signal iqWith loading moment TL;Jing after 13 computing of location/velocity integrated controller, permanent magnet synchronous electric is exported The given q axis current signals of machine 123rd subtractor receives given q axis current signalsBelieve with actual q shaft currents Number iq, the 3rd subtractor is to given q axis current signalsWith actual q axis current signal iqMake difference operation, and export q axles electricity The error signal of stream4th subtractor receives actual d axis current signal id, and received by master system given D axis current signalsBoth are made after difference operation by the 4th subtractor, export the error signal of d shaft currents Pi regulator 3 receives the error signal of q shaft currents respectivelyWith the error signal of d shaft currentsIt is defeated respectively Jing after computing Go out the given u of q shaft voltagesqU given with d shaft voltagesd;IPark conversion modules 4 receive the given u of q shaft voltagesqU given with d shaft voltagesd, The component of voltage u under static two phase coordinate system is exported Jing after computingαWith component of voltage uβ;Space vector modulation module 5 receives static Component of voltage u under two phase coordinate systemsαWith component of voltage uβInput, exports the six tunnels control letter of three-phase inverter 6 Jing after computing Number;Six tunnel control signals of the output of 6 receiving space Vector Modulation module of three-phase inverter 5, drive permagnetic synchronous motor 12 to run.
When real work, step is as follows:
1st, permagnetic synchronous motor 12 is started working, after the motor shaft rotation of permagnetic synchronous motor 12, by HALL15 to standard Sensorless strategy module 14 exports discrete HALL position signallings
2nd, quasi- sensorless strategy module 14 receives the discrete HALL position signallings of HALL15And pi regulator 3 The given u of q shaft voltagesqU given with d shaft voltagesd, and by discrete HALL position signallingsQ shaft voltages give uqWith d axles electricity The given u of pressuredComputing is carried out, actual rotor-position signal θ after computing, is drawnm, actual d axis current signal id, actual q axles Current signal iqWith the loading moment T of permagnetic synchronous motor 12L.Quasi- sensorless strategy module 14 is real to the output of the 5th subtractor The rotor-position signal θ on borderm, actual q axis current signal i are exported to location/velocity integrated controller 13qAnd loading moment TL, actual q axis current signal i are exported to the 3rd subtractorq, actual d axis current signal i are exported to the 4th subtractord, to IPark conversion modules 4 export actual rotor-position signal θm
3rd, by master system, to the position signalling θ that the input of the 5th subtractor is givenm *, then by standard without sensor control Molding block 14 receives actual rotor-position signal θm(described in step 1, quasi- sensorless strategy module 14 is to the 5th subtractor The actual rotor-position signal θ of outputm).5th subtractor is again by given position signalling θm *With actual rotor-position signal θmCarry out making difference operation, draw position error signal θm *m;5th subtractor is by position error signal θm *mExport to position/ Speed integrated controller 13.
4th, location/velocity integrated controller 13 receives the position error signal θ of the 5th subtractor outputm *m, quasi- nothing The q axis current signal i of the reality of the output of sensor control block 14qWith loading moment TL.Jing location/velocity integrated controllers The given q axis current signals of permagnetic synchronous motor 12 are drawn after 13 computingsLocation/velocity integrated controller 13 is to the 3rd The given q axis current signals of subtractor output
5th, the given q axis current signals of the 3rd subtractor receiving position/output of speed integrated controller 13And The q axis current signal i of the reality of the quasi- output of sensorless strategy module 14q.3rd subtractor is to given q axis current signals With actual q axis current signal iqThe error signal of q shaft currents is drawn after carrying out computing3rd subtractor is to pi regulator The error signal of 3 output q shaft currents
6th, the 4th subtractor receives the given d axis current signals that master system is providedAnd standard is without sensor control The q axis current signal i of the reality of the output of molding block 14q;The d axis current signals that 4th subtractor will giveAnd reality Q axis current signal iqAfter carrying out making difference operation, the error signal of d shaft currents is drawn4th subtractor is to pi regulator 3 The error signal of output d shaft currents
7th, pi regulator 3 receives the error signal of the q shaft currents of the 3rd subtractor outputAnd the 4th subtractor it is defeated The error signal of the d shaft currents for going outAfter the computing of pi regulator 3, the given u of q shaft voltages is obtainedqGive with d shaft voltages Determine ud;Pi regulator 3 gives u to the output q shaft voltages of IPark conversion modules 4qU given with d shaft voltagesd
8th, IPark conversion modules 4 are respectively received the given u of q shaft voltages of the output of pi regulator 3qU given with d shaft voltagesd, And the rotor-position signal θ of the reality of the quasi- output of sensorless strategy module 14m.After the computing of IPark conversion modules 4 Draw the component of voltage u under static two phase coordinate systemαWith component of voltage uβ;IPark conversion modules 4 are to space vector modulation module 5 Export the component of voltage u under static two phase coordinate systemαWith component of voltage uβ
9th, space vector modulation module 5 receives the voltage point under static two phase coordinate system of the output of IPark conversion modules 4 Amount uαWith component of voltage uβ.Six tunnel control signals of three-phase inverter 6 are drawn Jing after 5 computing of space vector modulation module.Swear in space Amount modulation module 5 exports six tunnel control signals of three-phase inverter 6 to three-phase inverter 6.
10th, three-phase inverter 6 receives six tunnel control signals of the three-phase inverter 6 of the output of space vector modulation module 5, And the operation of permagnetic synchronous motor 12 is driven by six tunnel control signals of three-phase inverter 6.
Above-described pi regulator 3, IPark conversion modules 4, space vector modulation module 5, three-phase inverter 6, permanent magnetism Synchronous motor 12, HALL15, the 3rd subtractor, the 4th subtractor, the 5th subtractor and permagnetic synchronous motor 12 are existing public affairs Know technology.
During the standard of above-described permagnetic synchronous motor is without sensing station Servocontrol device, by standard without sensor control Molding block 14 substitutes traditional based on biphase current sensor and photoelectric encoder 11 (as shown in Figure 1), rotary transformer equipotential The control program of sensor is put, position and the speed of traditional servo system are substituted using the monocyclic controller 13 of location/velocity integration Spend two ring control structures.
By the following description, quasi- sensorless strategy module 14 and location/velocity integrated controller are further described The operation method of 13 (the monocyclic controllers of location/velocity integration):
First, location/velocity integrated controller 13 of the invention, its implementation are as follows:
The first step, obtain motor (permagnetic synchronous motor 12) the equation of motion be:
Second step, the state for taking system (permagnetic synchronous motor of the present invention is accurate without sensing station Servocontrol device) become Amount:
Convolution (1), formula (2):
3rd step, orderF=-TL/ J, the then state equation of the positional servosystem of permagnetic synchronous motor 12 For:
4th step, the sliding-mode surface for taking similar PD actuators are designed as:
S=c1x1+c2x2 (5)
5th step, chooses control law:
6th step, the control function for obtaining quadrature axis current is:
Wherein,It is given for quadrature axis current,
The location/velocity overall-in-one control schema of the monocyclic controller of location/velocity integration 13 is achieved that according to above formula. Wherein, above-described TeFor motor (permagnetic synchronous motor 12) output torqueωmFor motor (permanent magnet synchronous electric Machine actual speed 12), TLFor loading moment, θmFor actual rotor-position signal, all by quasi- sensorless strategy module 14 Estimation is obtained;J is motor inertia, and P is motor number of pole-pairs, ΨaFor permanent magnet flux linkage, it is all the correlation ginseng of permagnetic synchronous motor 12 Number,For the position signalling for giving, c1And c1For constant.
Implement step as follows:
1st, by the 5th subtractor to actual rotor-position signal θmWith given position signalling θm *Position is drawn after computing Put error signal θm *m;By the 5th subtractor by position error signal θm *mGive location/velocity integrated controller 13;
2nd, the formula in accordance with the above of location/velocity integrated controller 13Calculate Go out given q axis current signals
2nd, quasi- sensorless strategy module 14 of the invention, its implementation are as follows:
Phase current reconstruction module 16 and full dimension that above-described quasi- sensorless strategy module 14 is connected by mutual signal Observer rotor position information estimation block 17 constitutes (as shown in Figure 4):
The implementation of phase current reconstruction module 16 is as follows:
Above-described ωmFor the actual speed of motor (permagnetic synchronous motor 12);La、Φa、RaRespectively motor (permanent magnetism Synchronous motor inductance 12), magnetic linkage and resistance;uqGiven, the u for q shaft voltagesdIt is given for d shaft voltages.
Implement step as follows:
1st, by q shaft voltages given uq, the given u of d shaft voltagesdExport with full micr oprocessorism rotor position information estimation block 17 Motor actual speed signal ωmIt is input in phase current reconstruction module 16;
2nd, phase current reconstruction module 16 calculates the quadrature axis current of 12 reality of permagnetic synchronous motor according to formula (8)-(11) Given idI given with actual direct-axis currentq
The computing formula of full micr oprocessorism rotor position information estimation block 17 is as follows:
Equally by Electrical Motor volume, the kinematical equation of motor (permagnetic synchronous motor 12) is:
Jdωm/dt+Bmωm+TL=Tε (12)
m/ dt=ωm (13)
Above-described ωmActual speed (i.e. rotor machinery angular velocity) for motor (permagnetic synchronous motor 12);θmFor reality The rotor-position signal (i.e. mechanical Angle Position) on border;TeFor the output torque of permagnetic synchronous motor 12;TL(bear for loading moment Carry perturbing torque);J is rotary inertia;BmViscous damping coefficient.
Because the sample frequency of controller is significantly larger than the transformation period of perturbing torque, TL(bear for loading moment Carry perturbing torque) as a state variable, it may be assumed that it is a constant, i.e. TLFor loading moment (i.e. disturbance load torque) Time differential is zero.
If:
dTL/ dt=0 (14)
The dynamics state equation of motor is obtained by (12) and (13):
Wherein,
Output torque T of the input variable of formula (15) for permagnetic synchronous motor 12e, state variable is actual rotor-position Signal θm(mechanical Angle Position), actual speed ω of permagnetic synchronous motor 12m(mechanical angular velocity) and TL(bear for loading moment Carry perturbing torque), output variable is actual rotor-position signal θm(mechanical Angle Position).Formula (5) can be written as formula (16);
By dynamics state equation (15), an omnidirectional vision can be set up, see formula (17):
Wherein,For estimative state variable,For feedback of status Gain battle array.
Obtained by (15)-(16):
Wherein,For observation error.Its characteristic equation is:
Det_sI- (A-KC)]=s3+(k1+Bm/J)s2+(k2+k1Bm/J)s-k3/ J=0 (19)
Appropriate K is selected, makes (A-KC) to have stable, appropriate eigenvalue,When, with x (t), u (t) andNothing Close;WithIt is unrelated.
According to specified expectation limit α, β, γ, then the desired character multinomial of observer is:
s3-(α+β+γ)s2+ (α β+β γ+γ α) s- α β γ=0 (20)
Obtained by (19) and (20):
Assume Bm=0, and α=β=γ, then formula (21) it is rewritable into:
Formula (17) it is rewritable into:
I.e.:
According to the desired characteristic of system, the position for selecting limit to be located, observer is constructed according to formula (24), you can observe Loading moment TL, actual rotor-position signal θmActual speed ω of (mechanical Angle Position) and permagnetic synchronous motor 12m(machinery Angular velocity) value.
The method realizes that step is as follows:
1st, HALL15 is by discrete HALL position signallingsFeed back to full micr oprocessorism rotor position information estimation block 17 In;
2nd, full micr oprocessorism rotor position information estimation block 17 estimates the reality of permagnetic synchronous motor 12 according to formula Rotor-position (actual rotor-position signal θm), actual speed ωmWith loading moment TL
3rd, q shaft voltages are given u by pi regulator 3qU given with d shaft voltagesdIt is input to (q axles in phase current reconstruction module 16 Voltage gives uqU given with d shaft voltagesdObtained by pi regulator 3, pi regulator 3 is existing known technology);
4th, phase current reconstruction module 16 estimates the d axis current signals of the reality of permagnetic synchronous motor 12 according to formula (24) idWith actual q axis current signal iq, and by actual q axis current signal iqFeed back to full micr oprocessorism rotor position information to estimate Calculate module 17.
Finally, in addition it is also necessary to it is noted that listed above is only a specific embodiment of the invention.Obviously, the present invention Above example is not limited to, there can also be many deformations.One of ordinary skill in the art can be straight from present disclosure The all deformations derived or associate are connect, protection scope of the present invention is considered as.

Claims (3)

1. permagnetic synchronous motor standard is characterized in that without sensing station Servocontrol device:Become including pi regulator (3), IPark Mold changing block (4), space vector modulation module (5), three-phase inverter (6), permagnetic synchronous motor (12), HALL (15), standard are without biography Sensor control module (14), the 3rd subtractor, the 4th subtractor, the 5th subtractor and location/velocity integrated controller (13); The HALL (15) is HALL sensors;
The HALL (15) is arranged on permagnetic synchronous motor (12) body, when permagnetic synchronous motor (12) works, by HALL (15) discrete HALL position signallings are exported;
The quasi- sensorless strategy module (14) receives the given u of discrete HALL position signallings, q shaft voltagesqGive with d shaft voltages Determine ud, Jing after the computing of quasi- sensorless strategy module (14), export actual rotor-position signal θm, actual d axis current signals id, actual q axis current signal iqWith the loading moment T of permagnetic synchronous motor (12)L
5th subtractor receives given position signal θm *With rotor-position signal θm, and both are made with difference operation outgoing position Error signal θm *m
The location/velocity integrated controller (13) receives position error signal θm *m, actual q axis current signal iqWith it is negative Carry torque TL, the given q axis current signals of Jing computings output permagnetic synchronous motor (12)
3rd subtractor receives given q axis current signalsWith actual q axis current signal iq, and fortune is differed to both Calculate the error signal of output q shaft currents
4th subtractor receives given d axis current signalsWith actual d axis current signal id, and it is poor to both Computing exports the error signal of d shaft currents
The pi regulator (3) receives the error signal of d shaft currents respectivelyWith the error signal of q shaft currentsAnd Jing Computing exports the given u of q shaft voltages respectivelyqU given with d shaft voltagesd
The IPark conversion modules (4) receive the given u of q shaft voltagesqU given with d shaft voltagesd, the static biphase seat of Jing computings output Component of voltage u under mark systemαWith component of voltage uβ
The space vector modulation module (5) receives the component of voltage u under static two phase coordinate systemαWith component of voltage uβInput, Jing Six tunnel control signals of computing output three-phase inverter (6);
The six tunnel control signals that the three-phase inverter (6) receiving space Vector Modulation module (5) exports, drive permanent magnet synchronous electric Machine (12) runs;
Phase current reconstruction module (16) and full micr oprocessorism rotor that quasi- sensorless strategy module (14) is connected by mutual signal Positional information estimation block (17) is constituted;The full micr oprocessorism rotor position information estimation block (17) receives actual q axles Current signal iqWith discrete HALL position signallings;Full micr oprocessorism rotor position information estimation block (17) output is actual Rotor-position signal θm, loading moment TLWith actual speed ω of permagnetic synchronous motor (12)m;The phase current reconstruction module (16) receive actual speed ωm, the given u of q shaft voltagesqU given with d shaft voltagesd;Phase current reconstruction module (16) output is real The d axis current signal i on borderd, actual q axis current signal iq
The implementation method of the phase current reconstruction module (16) is as follows:
i d = 1 L a s + R a ( u d + L a ω m i q ) - - - ( 8 )
i q = 1 L a s + R a ( u q - Φ a ω m - L a ω m i d ) - - - ( 9 )
i d ( k + 1 ) = ( 1 - α ) / R α z - α [ u d ( k ) + L a ω m ( k ) i q ( k ) ] - - - ( 10 )
i q ( k + 1 ) = ( 1 - α ) / R α z - α [ u q ( k ) - Φ a ω m ( k ) - L a ω m ( k ) i d ( k ) ] - - - ( 11 )
The La、Φa、RaThe respectively inductance of permagnetic synchronous motor (12), permanent magnet flux linkage and resistance;
The implementation method of the full micr oprocessorism rotor position information algorithm estimation block (17) is as follows:
Equally by Electrical Motor volume, the kinematical equation of permagnetic synchronous motor (12) is:
Jdωm/dt+Bmωm+TL=T is 12)
m/ dt=ωm (13)
The BmFor viscous damping coefficient;The TeFor the output torque of permagnetic synchronous motor (12), J is motor inertia;
Because the sample frequency of controller is significantly larger than the transformation period of perturbing torque, disturbance load torque is used as a shape State variable, it may be assumed that it is zero that it is a constant, i.e. disturbance load torque time differential;
If:
dTL/ dt=0 (14)
The dynamics state equation of motor is obtained by (12) and (13):
d x / d t = A x + B u y = C x - - - ( 15 )
Wherein,
A = 0 1 0 0 - B m / J - 1 / J 0 0 0 , B = 0 1 / J 0 , x = θ m ω m T L , C = 1 0 0 , u = T a , y = θ m ;
Formula (15) can be written as formula (16):
d d t θ m ω m T L = 0 1 0 0 - B m / J - 1 / J 0 0 0 θ m ω m T L + 0 1 / J 0 T a - - - ( 16 )
By dynamics state equation (15), an omnidirectional vision can be set up, see formula (17):
d x ^ / d t = A x ^ + B u + K ( y - y ^ ) y ^ = C x ^ - - - ( 17 )
Wherein,For estimative state variable, K=[k1 k2 k3]TFor state feedback gain matrix;
Obtained by (15)-(16):
d x ~ / d t = ( A - K C ) x ~ - - - ( 18 )
Its characteristic equation is:
Det [sI- (A-KC)]=s3+(k1+Bm/J)s2+(k2+k1 Bm/J)s-k3/ J=0 (19)
Appropriate K is selected, makes (A-KC) to have stable, appropriate eigenvalue,When, with x (t), u (t) andIt is unrelated;WithIt is unrelated;
According to specified expectation limit α, β, γ, then the desired character multinomial of observer is:
s3-(α+β+γ)s2+ (α β+β γ+γ α) s- α β γ=0 (20)
Obtained by (19) and (20):
k 1 = - ( α + β + γ ) - B m / J k 2 = ( α β + β γ + γ α ) + ( α + β + γ ) B m / J + ( B m / J ) 2 k 3 = α β γ J - - - ( 21 )
Assume Bm=0, and α=β=γ, then formula (21) it is rewritable into:
k 1 = - 3 α k 2 = 3 α 2 k 3 = α 3 J - - - ( 22 )
Formula (17) it is rewritable into:
d d t θ ^ m ω ^ m T ^ L = 0 1 0 0 - B m / J ^ - 1 / J ^ 0 0 0 θ ^ m ω ^ m T ^ L + 0 1 / J ^ 0 T a + k 1 k 2 k 3 ( θ m - θ ^ m ) - - - ( 23 )
I.e.:
d θ ^ m / d t = ω ^ + k 1 × ( θ m - θ ^ m ) d ω ^ m / d t = - B m / J ^ × ω ^ m - 1 / J ^ × T ^ L + 1 / J ^ × T a + k 2 × ( θ m - θ ^ m ) d T ^ L / d t = k 3 × ( θ m - θ ^ m ) - - - ( 24 )
According to the desired characteristic of system, the position for selecting limit to be located, observer is constructed according to formula (24), you can observe load Torque TL, permagnetic synchronous motor (12) actual speed ωmThe actual rotor-position signal θ with permagnetic synchronous motor (12)m's Value.
2. permagnetic synchronous motor standard, without sensing station Servocontrol device, is characterized in that according to claim 1:Position/ The implementation method of speed integrated controller (13) is as follows;
The first step, obtain motor the equation of motion be:
T a - T L = J dω m d t - - - ( 1 )
Second step, takes the state variable of system:
x 1 = θ m - θ m * x 2 = x · 1 = ω m - - - ( 2 )
Convolution (1), formula (2):
x · 1 = x 2 = ω m x · 2 = T a - T L J = 1 J [ 3 2 PΦ a i q - T L ] - - - ( 3 )
3rd step, orderF=-TL/ J, then the state equation of permanent magnetism synchronous electric machine position servo system be:
x · 1 x · 2 = 0 1 0 0 x 1 x 2 + 0 a i q + 0 F - - - ( 4 )
4th step, the sliding-mode surface for taking similar PD actuators are designed as:
S=c1x1+c2x2 (5)
5th step, chooses control law:
s · = - ϵ sgn s - k s - - - ( 6 )
6th step, the control function for obtaining quadrature axis current is:
i q * = - ϵ sgn s - k s - - - ( 7 )
It is describedFor give q axis current signals,
The overall-in-one control schema of location/velocity is achieved that according to above formula;
The P be motor number of pole-pairs, c1And c2For constant.
3. a kind of permagnetic synchronous motor standard for being carried out using claim 1 or 2 described devices is without sensing station SERVO CONTROL side Method, is characterized in that:
A, permagnetic synchronous motor (12) are started working, after the rotation of the motor shaft of permagnetic synchronous motor (12), by HALL (15) to Quasi- sensorless strategy module (14) exports discrete HALL position signallings;
B, quasi- sensorless strategy module (14) receive the discrete HALL position signallings and pi regulator (3) of HALL (15) The given u of q shaft voltagesqU given with d shaft voltagesd, and discrete HALL position signallings, q shaft voltages are given into uqGive with d shaft voltages Determine udComputing is carried out, actual rotor-position signal θ after computing, is drawnm, actual d axis current signal id, actual q shaft currents Signal iqWith the loading moment T of permagnetic synchronous motor (12)L;Quasi- sensorless strategy module (14) is real to the output of the 5th subtractor The rotor-position signal θ on borderm, actual q axis current signal i are exported to location/velocity integrated controller (13)qAnd carrying Square TL, actual q axis current signal i are exported to the 3rd subtractorq, actual d axis current signal i are exported to the 4th subtractord, Actual rotor-position signal θ is exported to IPark conversion modules (4)m
C, by master system, to the position signalling θ that the input of the 5th subtractor is givenm *, then by quasi- sensorless strategy mould Block (14) receives actual rotor-position signal θm;5th subtractor is again by given position signalling θm *With actual rotor position Confidence θmCarry out making difference operation, draw position error signal θm *m;5th subtractor is by position error signal θm *mOutput Give location/velocity integrated controller (13);
D, location/velocity integrated controller (13) receive the position error signal θ of the 5th subtractor outputm *m, standard without biography The q axis current signal i of the reality that sensor control module (14) is exportedqWith loading moment TL;Jing location/velocity integrated controllers (13) the given q axis current signals of permagnetic synchronous motor (12) are drawn after computingLocation/velocity integrated controller (13) To the q axis current signals that the output of the 3rd subtractor is given
The given q axis current signals that e, the 3rd subtractor receiving position/speed integrated controller (13) are exportedAnd it is accurate The q axis current signal i of the reality that sensorless strategy module (14) is exportedq;3rd subtractor is to given q axis current signals With actual q axis current signal iqThe error signal of q shaft currents is drawn after carrying out computing3rd subtractor is to pi regulator (3) export the error signal of q shaft currents
The given d axis current signals of f, the 4th subtractor system reception systemAnd quasi- sensorless strategy module (14) output Reality d axis current signal id;The d axis current signals that 4th subtractor will giveAnd the d axis current signal i of realityd After carrying out making difference operation, the error signal of d shaft currents is drawn4th subtractor exports d shaft currents to pi regulator (3) Error signal
G, pi regulator (3) receive the error signal of the q shaft currents of the 3rd subtractor outputAnd the 4th subtractor output D shaft currents error signalAfter the computing of pi regulator (3), the given u of q shaft voltages is obtainedqGive with d shaft voltages Determine ud;Pi regulator (3) gives u to IPark conversion modules (4) output q shaft voltagesqU given with d shaft voltagesd
H, IPark conversion module (4) is respectively received the given u of q shaft voltages that pi regulator (3) is exportedqU given with d shaft voltagesd, And the rotor-position signal θ of reality that quasi- sensorless strategy module (14) exportsm;Through the fortune of IPark conversion modules (4) The component of voltage u under static two phase coordinate system is drawn after calculationαWith component of voltage uβ;IPark conversion modules (4) are adjusted to space vector Molding block (5) exports the component of voltage u under static two phase coordinate systemαWith component of voltage uβ
I, space vector modulation module (5) receive the voltage point under static two phase coordinate system that IPark conversion modules (4) are exported Amount uαWith component of voltage uβ;Six tunnel control signals of three-phase inverter (6) are drawn Jing after space vector modulation module (5) computing;It is empty Between Vector Modulation module (5) to three-phase inverter 6 export three-phase inverter (6) six tunnel control signals;
J, three-phase inverter (6) receive the six tunnels control letter of the three-phase inverter (6) that space vector modulation module (5) is exported Number, and the operation of permagnetic synchronous motor (12) is driven by six tunnel control signals of three-phase inverter (6).
CN201310602081.8A 2013-11-21 2013-11-21 Permagnetic synchronous motor standard is without sensing station Servocontrol device and method Active CN104333285B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310602081.8A CN104333285B (en) 2013-11-21 2013-11-21 Permagnetic synchronous motor standard is without sensing station Servocontrol device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310602081.8A CN104333285B (en) 2013-11-21 2013-11-21 Permagnetic synchronous motor standard is without sensing station Servocontrol device and method

Publications (2)

Publication Number Publication Date
CN104333285A CN104333285A (en) 2015-02-04
CN104333285B true CN104333285B (en) 2017-04-05

Family

ID=52407956

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310602081.8A Active CN104333285B (en) 2013-11-21 2013-11-21 Permagnetic synchronous motor standard is without sensing station Servocontrol device and method

Country Status (1)

Country Link
CN (1) CN104333285B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102014212383A1 (en) * 2014-06-27 2015-12-31 Robert Bosch Gmbh Method for operating an electric machine, electric machine
CN105932919B (en) * 2016-06-07 2018-11-30 深圳市百盛传动有限公司 A kind of permanent magnet synchronous motor position control method based on motor model
CN106208839B (en) * 2016-07-25 2019-04-02 河南师范大学 A kind of permanent magnet synchronous motor voltage-type location measurement method used for wind power generation and device
KR101852754B1 (en) * 2016-10-27 2018-04-27 엘에스산전 주식회사 Sensoelss control system for permanent magnet synchronous machine
CN106292274B (en) * 2016-11-02 2019-05-10 湖南工业大学 A kind of estimation method cascading sliding mode observer and its adhesion coefficient and derivative
KR102676718B1 (en) * 2018-03-14 2024-06-19 현대자동차주식회사 Control method and control system of motor rotation speed
CN108880358A (en) * 2018-07-02 2018-11-23 哈尔滨理工大学 Method for controlling permanent magnet synchronous motor and device based on angular displacement without Time Delay Observer
CN108923712B (en) * 2018-08-28 2019-12-31 青岛大学 Single-loop control method, device and system for rotating speed of permanent magnet synchronous motor
CN110190793B (en) * 2019-05-31 2020-12-25 东南大学 Two-degree-of-freedom numerical control machine tool and control system and positioning method thereof
CN110323983B (en) * 2019-07-23 2021-05-11 广东工业大学 Current decoupling method, device, equipment and medium for permanent magnet synchronous motor
CN111416560B (en) * 2020-04-27 2024-04-05 廖铉泓 Permanent magnet synchronous motor control method and system without current sensor
CN115313950B (en) * 2022-08-19 2023-12-26 中山市瑞思科传动科技有限公司 Driving system of servo motor

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101867374A (en) * 2009-04-17 2010-10-20 索尼公司 Analog-digital converter, solid photographic device and camera arrangement
CN101931366A (en) * 2010-08-03 2010-12-29 浙江大学 Super capacitor-based energy-saving drive circuit of motor and control method
CN103401501A (en) * 2013-04-15 2013-11-20 湖南大学 Permanent magnet synchronous motor (PMSM) servo system control method based on fuzzy and active disturbance rejection control

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101867374A (en) * 2009-04-17 2010-10-20 索尼公司 Analog-digital converter, solid photographic device and camera arrangement
CN101931366A (en) * 2010-08-03 2010-12-29 浙江大学 Super capacitor-based energy-saving drive circuit of motor and control method
CN103401501A (en) * 2013-04-15 2013-11-20 湖南大学 Permanent magnet synchronous motor (PMSM) servo system control method based on fuzzy and active disturbance rejection control

Also Published As

Publication number Publication date
CN104333285A (en) 2015-02-04

Similar Documents

Publication Publication Date Title
CN104333285B (en) Permagnetic synchronous motor standard is without sensing station Servocontrol device and method
CN103997272B (en) The load disturbance compensation device of permagnetic synchronous motor and method
CN106208865B (en) More permanent magnet synchronous motor Virtual-shaft control methods based on Load Torque Observer
CN106887976A (en) Consider many permagnetic synchronous motor deviation coupling control methods of acceleration
CN105897097B (en) Permagnetic synchronous motor current predictive control method and device
CN108092567A (en) A kind of Speed control of permanent magnet synchronous motor system and method
CN107359837A (en) Torsion control system of synchronization generator with everlasting magnetic and method based on sliding mode observer and Active Disturbance Rejection Control
CN103051274B (en) Variable damping-based passive control method for two-degree-of-freedom permanent magnetic synchronous motor
CN103269199B (en) Electric car induction motor torque current setting device
CN106059419B (en) A kind of permanent magnet synchronous motor parallel connection vector control scheme
CN107370431A (en) A kind of industrial robot obscures Auto-disturbance-rejection Control with permagnetic synchronous motor
CN104885356B (en) Arrangement for controlling induction motor and inductance motor control method
CN110460280A (en) A kind of permasyn morot control method based on sliding formwork load torque observer
CN104953916A (en) Novel speed controller based on speed regulating system of permanent magnet synchronous motor
CN109217766A (en) The independent reversed decoupling control system of induction-type bearingless motor
CN106788045A (en) A kind of permagnetic synchronous motor model prediction PI changeable weight concurrency control methods
CN102843088A (en) Fuzzy sliding mode and neural network based composite control system of permanent synchronous motor
CN106992733A (en) Vehicle-mounted internal permanent magnet synchronous motor control method
CN106849809A (en) A kind of SVM DTC motor control algorithms based on magnetic linkage error vector method
CN105071730A (en) Bearing-free asynchronous motor stator orientated reverse closed-loop control system taking current dynamics into consideration
CN110401390A (en) The random command filtering fuzzy control method of permanent magnet synchronous motor based on observer
CN104852658A (en) Permanent magnet synchronous motor decoupling vector control device in two-phase stationary coordinate system and method thereof
CN106788053A (en) Electric motor car permagnetic synchronous motor systematic error compensation control method based on observer
CN106130429A (en) Bearing-free permanent magnet synchronous motor predictive controller and building method
CN109150045A (en) The independent Inverse Decoupling method of induction-type bearingless motor

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