CN113985835A - Method for estimating mechanical characteristic parameters of large-antenna radar servo system - Google Patents

Method for estimating mechanical characteristic parameters of large-antenna radar servo system Download PDF

Info

Publication number
CN113985835A
CN113985835A CN202111072947.XA CN202111072947A CN113985835A CN 113985835 A CN113985835 A CN 113985835A CN 202111072947 A CN202111072947 A CN 202111072947A CN 113985835 A CN113985835 A CN 113985835A
Authority
CN
China
Prior art keywords
friction
servo system
calculating
parameters
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.)
Granted
Application number
CN202111072947.XA
Other languages
Chinese (zh)
Other versions
CN113985835B (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.)
Beijing Institute of Radio Measurement
Original Assignee
Beijing Institute of Radio Measurement
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 Beijing Institute of Radio Measurement filed Critical Beijing Institute of Radio Measurement
Priority to CN202111072947.XA priority Critical patent/CN113985835B/en
Publication of CN113985835A publication Critical patent/CN113985835A/en
Application granted granted Critical
Publication of CN113985835B publication Critical patent/CN113985835B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0259Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterized by the response to fault detection
    • G05B23/0262Confirmation of fault detection, e.g. extra checks to confirm that a failure has indeed occurred
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/20Pc systems
    • G05B2219/24Pc safety
    • G05B2219/24065Real time diagnostics

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The embodiment of the invention discloses a method for estimating mechanical characteristic parameters of a large antenna radar servo system, which comprises the following steps: calculating friction torque by using the position steady-state error; and calculating an estimated value of the friction characteristic parameter and an estimated value of the rotational inertia. The method adopts an off-line estimation method, is simple and is suitable for a large-antenna radar servo system with high requirements on real-time performance and execution efficiency.

Description

Method for estimating mechanical characteristic parameters of large-antenna radar servo system
Technical Field
The invention relates to the field of radar servo systems. And more particularly, to a method for estimating mechanical characteristic parameters of a large antenna radar servo system.
Background
Mechanical characteristic parameters, including friction modeling parameters and antenna rotational inertia, are important parameters of a radar servo system. The friction can reduce the speed stability and the tracking precision of a servo system, and the rotational inertia of the antenna directly influences various performances of servo closed-loop control. At present, few researches are carried out on the estimation method of mechanical characteristic parameters of a large-antenna radar servo system, the problems that friction models are inaccurate, the parameters are nonlinear and difficult to estimate exist in the researches of other servo systems and the like are solved, the estimation of the rotational inertia generally depends on the design weight data of a load, and a certain error exists between the actual load rotational inertia and the rotational inertia. Most of the existing estimation methods are on-line estimation, and the used algorithm is complex and is not beneficial to engineering realization. From the technical and cost aspects, the method is not suitable for a large-antenna radar servo system with high requirements on real-time performance and execution efficiency.
Disclosure of Invention
In view of this, a first embodiment of the present invention provides a method for estimating mechanical characteristic parameters of a large antenna radar servo system, including:
calculating friction torque by using the position steady-state error;
and calculating an estimated value of the friction characteristic parameter and an estimated value of the rotational inertia.
In a specific embodiment, the method further comprises:
and acquiring preset parameters of a servo system control structure, and calculating a position steady-state error according to the preset parameters.
In one embodiment, the friction torque is:
F=-essKpKv
wherein ,essFor position steady state error, KpFor proportional control, KvIs the gain.
In one embodiment, the preset parameters of the servo system control structure include: position control parameters, speed loop parameters, and antenna loading, wherein,
setting the position control parameter to proportional control KpThe speed loop parameter is the gain KvAnd time constant T, equivalent to an inertia element
Figure BDA0003261070500000021
The antenna is loaded with
Figure BDA0003261070500000022
J is moment of inertia.
In one embodiment, the position steady state error is
Figure BDA0003261070500000023
wherein ,
Figure BDA0003261070500000024
in one embodiment, the calculating the friction characteristic parameter estimation value includes:
setting the rotation speed of the servo system to be omega respectivelyi(i is 1, …, n), corresponding n groups of friction torque can be obtained
Fi=-essiKpKv(i=1,…,n)
At n sets of friction moments Fi(i-1, …, n) to find the minimum value
Figure BDA0003261070500000025
Corresponding rotational speed
Figure BDA0003261070500000026
As the Stribeck velocity ωsIs recorded as
Figure BDA0003261070500000027
Is calculated as [ Fc,Fs,b]TLinear equation of friction torque as parameter
Fi=(1-μi)FciFs+bωi
Figure BDA0003261070500000028
wherein ,Fc、Fs
Figure BDA0003261070500000029
And b are four characteristic parameters of Coulomb friction moment, maximum static friction moment, Stribeck speed and viscous friction coefficient, omegaiIn order to control the speed of the servo system,
the regression matrix Φ is constructed as follows:
Figure BDA00032610705000000210
calculating [ F ]c,Fs,b]TIs recorded as
Figure BDA00032610705000000211
Figure BDA00032610705000000212
To obtain Fc、Fs、ωsAnd b estimated values of four friction characteristic parameters.
In one embodiment, the calculating the estimated value of the moment of inertia includes:
enabling the servo system to rotate according to a preset curve track, and setting the driving torque at T1 and T2 as Tt1 and Tt2If the preset acceleration is a, the estimated value of the moment of inertia is calculated
Figure BDA0003261070500000031
The invention has the following beneficial effects:
the method calculates the friction torque by utilizing the position steady-state error, calculates the estimated value of the friction characteristic parameter and the estimated value of the rotational inertia, and determines the mechanical characteristic parameter of the large-antenna radar servo system, thereby being convenient for the engineering application of the large-antenna radar servo system.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 shows a flow chart of a method for estimating mechanical characteristic parameters of a large antenna radar servo system according to an embodiment of the invention.
Fig. 2 shows a block diagram of a large antenna radar servo system according to an embodiment of the present invention.
Fig. 3 shows a schematic diagram of a preset curve trajectory according to an embodiment of the invention.
Detailed Description
In order to make the technical solutions and advantages of the present invention clearer, the following will describe embodiments of the present invention in further detail with reference to the accompanying drawings.
As shown in fig. 1, a method for estimating mechanical characteristic parameters of a large antenna radar servo system includes:
calculating friction torque by using the position steady-state error;
and calculating an estimated value of the friction characteristic parameter and an estimated value of the rotational inertia.
In a specific embodiment, the method further comprises:
and acquiring preset parameters of a large antenna radar servo system control structure, and calculating the position steady-state error according to the preset parameters.
The large antenna radar is a radar with an antenna diameter exceeding 10 meters, the servo system is shown in fig. 2, and the preset parameters of the servo system control structure include: position control parameters, speed loop parameters, and antenna loading, wherein,
setting the position control parameter to proportional control KpThe speed loop parameter is the gain KvAnd time constant T, which can be equivalent to an inertia element
Figure BDA0003261070500000041
The antenna is loaded with
Figure BDA0003261070500000042
J is moment of inertia.
When the servo system rotates at a constant speed, the friction torque F is a constant value, and a steady-state error e is caused by the constant valuessIs composed of
Figure BDA0003261070500000043
wherein ,
Figure BDA0003261070500000044
the friction torque is thus obtained as:
F=-essKpKv
wherein ,essFor position steady state error, KpFor proportional control, KvIs the gain.
Setting the rotation speed of the servo system to be omega respectivelyi(i is 1, …, n), corresponding n groups of friction torque can be obtained
Fi=-essiKpKv(i=1,…,n)
At n sets of friction moments Fi(i-1, …, n) to find the minimum value
Figure BDA0003261070500000045
Corresponding rotational speed
Figure BDA0003261070500000046
As the Stribeck velocity ωsIs recorded as
Figure BDA0003261070500000047
Is calculated as [ Fc,Fs,b]TLinear equation of friction torque as parameter
Fi=(1-μi)FciFs+bωi
Figure BDA0003261070500000048
wherein ,Fc、Fs
Figure BDA0003261070500000049
And b are four characteristic parameters of Coulomb friction moment, maximum static friction moment, Stribeck speed and viscous friction coefficient, omegaiIn order to control the speed of the servo system,
the regression matrix Φ is constructed as follows:
Figure BDA00032610705000000410
calculating [ F ]c,Fs,b]TIs recorded as
Figure BDA00032610705000000411
Figure BDA00032610705000000412
To obtain Fc、Fs、ωsAnd b estimated values of four friction characteristic parameters.
The calculating the estimated value of the moment of inertia comprises:
the servo system is made to rotate according to a preset curve track as shown in FIG. 3, and the driving torques at the time T1 and the time T2 are respectively set as Tt1 and Tt2If the preset acceleration is a, the estimated value of the moment of inertia is calculated
Figure BDA0003261070500000051
The estimated value of the moment of inertia J is obtained, mechanical characteristic parameters of the large-antenna radar servo system are completely clear, engineering application of the large-antenna radar servo system is facilitated, and the method adopts an off-line estimation method, is simple in calculation method, and is suitable for the large-antenna radar servo system with high requirements on instantaneity and execution efficiency.
It should be understood that the above-mentioned embodiments of the present invention are only examples for clearly illustrating the present invention, and are not intended to limit the embodiments of the present invention, and it will be obvious to those skilled in the art that other variations or modifications may be made on the basis of the above description, and all embodiments may not be exhaustive, and all obvious variations or modifications may be included within the scope of the present invention.

Claims (7)

1. A method for estimating mechanical characteristic parameters of a large antenna radar servo system is characterized by comprising the following steps:
calculating friction torque by using the position steady-state error;
and calculating an estimated value of the friction characteristic parameter and an estimated value of the rotational inertia.
2. The method of claim 1, further comprising:
and acquiring preset parameters of a servo system control structure, and calculating a position steady-state error according to the preset parameters.
3. The method of claim 1, wherein the friction torque is:
F=-essKpKv
wherein ,essFor position steady state error, KpFor proportional control, KvIs the gain.
4. The method of claim 2, wherein the preset parameters of the servo control structure comprise: position control parameters, speed loop parameters, and antenna loading, wherein,
setting the position control parameter to proportional control KpThe speed loop parameter is the gain KvAnd time constant T, equivalent to an inertia element
Figure FDA0003261070490000011
The antenna is loaded with
Figure FDA0003261070490000012
J is moment of inertia.
5. The method of claim 3, wherein the position steady state error is
Figure FDA0003261070490000013
wherein ,
Figure FDA0003261070490000014
6. the method of claim 2, wherein the calculating the friction characteristic parameter estimate comprises:
setting the rotation speed of the servo system to be omega respectivelyi(i is 1, …, n), corresponding n groups of friction torque can be obtained
Fi=-essiKpKv(i=1,…,n)
At n sets of friction moments Fi(i-1, …, n) to find the minimum value
Figure FDA0003261070490000015
Corresponding rotational speed
Figure FDA0003261070490000016
As the Stribeck velocity ωsIs recorded as
Figure FDA0003261070490000017
Is calculated as [ Fc,Fs,b]TLinear equation of friction torque as parameter
Fi=(1-μi)FciFs+bωi
Figure FDA0003261070490000018
wherein ,Fc、Fs
Figure FDA0003261070490000019
And b are four characteristic parameters of Coulomb friction moment, maximum static friction moment, Stribeck speed and viscous friction coefficient, omegaiIn order to control the speed of the servo system,
the regression matrix Φ is constructed as follows:
Figure FDA0003261070490000021
calculating [ F ]c,Fs,b]TIs recorded as
Figure FDA0003261070490000022
Figure FDA0003261070490000023
To obtain Fc、Fs、ωsAnd b estimated values of four friction characteristic parameters.
7. The method of claim 6, wherein the calculating the estimate of moment of inertia comprises:
enabling the servo system to rotate according to a preset curve track, and setting the driving torque at T1 and T2 as Tt1 and Tt2If the preset acceleration is a, the estimated value of the moment of inertia is calculated
Figure FDA0003261070490000024
CN202111072947.XA 2021-09-14 2021-09-14 Mechanical characteristic parameter estimation method for large antenna radar servo system Active CN113985835B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111072947.XA CN113985835B (en) 2021-09-14 2021-09-14 Mechanical characteristic parameter estimation method for large antenna radar servo system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111072947.XA CN113985835B (en) 2021-09-14 2021-09-14 Mechanical characteristic parameter estimation method for large antenna radar servo system

Publications (2)

Publication Number Publication Date
CN113985835A true CN113985835A (en) 2022-01-28
CN113985835B CN113985835B (en) 2023-09-29

Family

ID=79735781

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111072947.XA Active CN113985835B (en) 2021-09-14 2021-09-14 Mechanical characteristic parameter estimation method for large antenna radar servo system

Country Status (1)

Country Link
CN (1) CN113985835B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06225565A (en) * 1993-01-25 1994-08-12 Yaskawa Electric Corp Method of measuring load constant of motor drive system
CN202646656U (en) * 2012-05-07 2013-01-02 中国人民解放军63636部队 Gravity mismatching device for backlash elimination of precise measurement radar servo system
CN107102634A (en) * 2017-05-11 2017-08-29 北京理工大学 A kind of parameter Estimation and tracking and controlling method based on table servo system
CN110161862A (en) * 2019-06-05 2019-08-23 安徽大学 A kind of servo turntable friciton compensation control system and method based on switching system theory

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06225565A (en) * 1993-01-25 1994-08-12 Yaskawa Electric Corp Method of measuring load constant of motor drive system
CN202646656U (en) * 2012-05-07 2013-01-02 中国人民解放军63636部队 Gravity mismatching device for backlash elimination of precise measurement radar servo system
CN107102634A (en) * 2017-05-11 2017-08-29 北京理工大学 A kind of parameter Estimation and tracking and controlling method based on table servo system
CN110161862A (en) * 2019-06-05 2019-08-23 安徽大学 A kind of servo turntable friciton compensation control system and method based on switching system theory

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HUI YANG, ET AL: "Study on the friction torque test and identification algorithm for gimbal axis of an inertial stabilized platform", ORIGINAL ARTICLE, pages 1 - 10 *
向红标等: "基于Stribeck简化模型的伺服系统转动惯量的辨识", 第十七届中国科协年会——分3面向智能制造的先进测控技术学术会议论文集, pages 1 - 7 *
谭文斌等: "应用稳态误差分析辨识 LuGre 模型参数", 光学精密工程, pages 664 - 670 *

Also Published As

Publication number Publication date
CN113985835B (en) 2023-09-29

Similar Documents

Publication Publication Date Title
CN109927032B (en) Mechanical arm track tracking control method based on high-order sliding-mode observer
CA1098194A (en) Method and means for increasing the stiffness of limited frequency servo systems
CN101262199B (en) Motor controller and motor control system
CN111152225A (en) Uncertain mechanical arm fixed time trajectory tracking control method with input saturation
CN108803324A (en) Multi-joint industrial machinery arm contragradience finite time sliding-mode control
CN104111664B (en) In the method that speed ring overcomes motor dead band raising radar tracking precision
CN107561935A (en) Motor position servo system friciton compensation control method based on multilayer neural network
CN105171758B (en) Self-adaptive finite time convergence sliding-mode control method of robot
US11519386B2 (en) Individual pitch control for wind turbines
CN111788040A (en) Kinetic parameter identification method of robot, robot and storage device
CN110109352B (en) Fixed time self-adaptive attitude control method of three-degree-of-freedom four-rotor aircraft
CN104639003A (en) AC servo system rotation inertia identification method
CN111965976B (en) Robot joint sliding mode control method and system based on neural network observer
CN104699108A (en) Multi-rotor craft control allocation method
CN113128018B (en) Friction force calculation method and device, robot and readable storage medium
CN108469734A (en) Consider the motor servo system Auto-disturbance-rejection Control of state constraint
CN111572818A (en) Magnetic suspension control moment gyroscope frame rate servo system and control method
CN113985835A (en) Method for estimating mechanical characteristic parameters of large-antenna radar servo system
CN210380707U (en) Motor control system based on known parameters
CN108170030B (en) Power system control distribution method and related device
Magar et al. Direct adaptive torque control for maximizing the power captured by wind turbine in partial loading condition
CN114200842A (en) Flight control method of L1 self-adaptive dynamic inverse nonlinear fixed wing aircraft
CN109308008B (en) Active disturbance rejection control device with abnormality coping capability
CN113238572A (en) Preset-time quadrotor unmanned aerial vehicle attitude tracking method based on preset performance control
Ping et al. Improved LuGre-based friction modeling of the electric linear load simulator

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant