CN113722891A - Harmonic reducer transmission system torsional rigidity fitting method - Google Patents

Harmonic reducer transmission system torsional rigidity fitting method Download PDF

Info

Publication number
CN113722891A
CN113722891A CN202110911816.XA CN202110911816A CN113722891A CN 113722891 A CN113722891 A CN 113722891A CN 202110911816 A CN202110911816 A CN 202110911816A CN 113722891 A CN113722891 A CN 113722891A
Authority
CN
China
Prior art keywords
torsional rigidity
model
harmonic gear
torque
torsional
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.)
Pending
Application number
CN202110911816.XA
Other languages
Chinese (zh)
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 University of Technology
Original Assignee
Beijing University of Technology
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 University of Technology filed Critical Beijing University of Technology
Priority to CN202110911816.XA priority Critical patent/CN113722891A/en
Publication of CN113722891A publication Critical patent/CN113722891A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Gears, Cams (AREA)

Abstract

The invention discloses a torsional rigidity fitting method of a harmonic gear transmission system, which aims at the problems that the expression of the traditional torsional rigidity model of harmonic gear transmission is complex and mutton is symmetrical, otherwise, the invention improves the torsional rigidity model as follows: the invention describes the nonlinearity of the torsional rigidity by using cubic Taylor series expansion, and describes the hysteresis characteristic of the torsional rigidity by using a hyperbolic function to obtain the relation between the torque T and the torsional angle. The torsional rigidity model established by the method ensures higher fitting precision, and parameters in the model are easy to identify.

Description

Harmonic reducer transmission system torsional rigidity fitting method
Technical Field
The invention belongs to the technical field of dynamic characteristics of harmonic gear transmission systems, and particularly relates to a torsional rigidity fitting model of a harmonic reducer transmission system
Background
The harmonic gear transmission has the advantages of large transmission ratio, small return difference, small volume, high transmission efficiency and the like, and is gradually applied to various manufacturing equipment and measuring equipment, such as industrial robots, machine tools, optical scanners, laser mirror positioning mechanisms and the like. Since harmonic drive is a key drive component of the above-described electromechanical system, its performance has a large impact on the manufacturing or measurement accuracy of the system. Harmonic gear drives, however, are driven due to nonlinear torsional stiffness. The transmission precision of the system is reduced, the transmission stability is deteriorated, and even the production can not meet the working requirement. A number of researchers have performed such as fitting using cubic odd functions, and also fitting using Maxwell hysteresis models. The former method has simple model, can not express the hysteresis characteristic of the torsional rigidity, and the latter model has the disadvantages of complexity, large calculation amount and low efficiency.
Disclosure of Invention
Based on the background technology, the invention uses cubic Taylor series expansion to describe the nonlinearity of the torsional rigidity, uses a hyperbolic function to describe the hysteresis characteristic of the torsional rigidity, and obtains the relation between the torque T and the torsional angle. And a constant is introduced into the model, so that the problem that the former model is completely symmetrical is solved, and the fitting precision is improved.
The modeling method comprises the following steps:
the model fitting process is mainly as follows:
firstly, expressing the nonlinear relation between torque and torsion angle by adopting Taylor series expansion, and in order to ensure the fitting precision and reduce the calculated amount, the invention adopts three times of Taylor expansion to express:
T=a+b·θ+c·θ2+d·θ3
wherein T is torque, a, b, c, d are polynomial coefficients, and theta is torsion angle
The hysteresis of torsional stiffness in harmonic gear drive systems is represented by hyperbolic functions. The following were used:
T=Atanh(B·θ)
where T is torque, A, B is coefficient, and θ is torsion angle
And superposing the two parts in the formula to obtain an integral torsional rigidity fitting model, wherein the model comprises the following steps:
T=a+b·θ+c·θ2+d·θ3+Atanh(B·θ)
from the torsional stiffness curves of harmonic gear transmissions, it can be seen that the graph consists of two closed curves, so the expression for the two curves can be expressed as:
T1=a1+b1·θ+c1·θ2+d1·θ3+A1 tanh(B1·θ)
T2=a2+b2·θ+c2·θ2+d2·θ3+A2 tanh(B2·θ)
wherein, T1、T2Respectively loaded and unloaded.
The corresponding relation of various values between the torque and the torsion angle of the harmonic gear transmission system can be obtained through experiments, parameters in the model can be identified by utilizing an nlnfit function in Matlab, and the functional relation between the torque and the torsion angle can be obtained.
Because of dT/d theta, the expression obtained by the above formula is subjected to derivation to obtain an expression of torsional rigidity of the harmonic gear.
Drawings
FIG. 1 is a graph of the fitting results of torsional stiffness according to the present invention
Detailed Description
Based on the background technology, the invention uses cubic Taylor series expansion to describe the nonlinearity of the torsional rigidity, uses a hyperbolic function to describe the hysteresis characteristic of the torsional rigidity, and obtains the relation between the torque T and the torsional angle. And a constant is introduced into the model, so that the problem that the former model is completely symmetrical is solved, and the fitting precision is improved.
The modeling method comprises the following steps:
the model fitting process is mainly as follows:
firstly, expressing the nonlinear relation between torque and torsion angle by adopting Taylor series expansion, and in order to ensure the fitting precision and reduce the calculated amount, the invention adopts three times of Taylor expansion to express:
T=a+b·θ+c·θ2+d·θ3
the hysteresis of torsional stiffness in harmonic gear drive systems is represented by hyperbolic functions. The following were used:
T=Atanh(B·θ)
and superposing the two parts in the formula to obtain an integral torsional rigidity fitting model, wherein the model comprises the following steps:
T=a+b·θ+c·θ2+d·θ3+Atanh(B·θ)
from the torsional stiffness curves of harmonic gear transmissions, it can be seen that the graph consists of two closed curves, so the expression for the two curves can be expressed as:
T=a1+b1·θ+c1·θ2+d1·θ3+A1 tanh(B1·θ)
T=a2+b2·θ+c2·θ2+d2·θ3+A2 tanh(B2·θ)
wherein, T1、T2Respectively loaded and unloaded.
The corresponding relation of various values between the torque and the torsion angle of the harmonic gear transmission system can be obtained through experiments, parameters in the model can be identified by utilizing an nlnfit function in Matlab, and the functional relation between the torque and the torsion angle can be obtained.
Because of dT/d theta, the expression obtained by the above formula is subjected to derivation to obtain an expression of torsional rigidity of the harmonic gear.

Claims (1)

1. A torsional rigidity fitting method of a harmonic gear transmission system is characterized by comprising the following steps: describing the nonlinearity of the torsional rigidity by using cubic Taylor series expansion, and describing the hysteresis characteristic of the torsional rigidity by using a hyperbolic function to obtain the relation between the torque T and the torsional angle; moreover, constants are introduced into the model, so that the problem that the former model is completely symmetrical is solved, and the fitting precision is improved;
the fitting procedure is as follows:
(1) firstly, expressing the nonlinear relation between the torque and the torsion angle by adopting Taylor series expansion, and in order to ensure the fitting precision and reduce the calculated amount, expressing by adopting three times of Taylor expansion:
T=a+b·θ+c·θ2+d·θ3
(2) a hyperbolic function is used for expressing a hysteresis phenomenon of torsional rigidity in a harmonic gear transmission system; the following were used:
T=A tanh(B·θ)
(3) and superposing the two parts in the formula to obtain an integral torsional rigidity fitting model, wherein the model comprises the following steps:
T=a+b·θ+c·θ2+d·θ3+A tanh(B·θ)
(4) the torsional stiffness curve of a harmonic gear drive is composed of two closed curves, so that the expression of the two curves can be expressed as:
T=a1+b1·θ+c1·θ2+d1·θ3+A1 tanh(B1·θ)
T=a2+b2·θ+c2·θ2+d2·θ3+A2 tanh(B2·θ)
(5) obtaining the corresponding relation of all values between the torque and the torsion angle of the harmonic gear transmission system through experiments, and identifying parameters in the model by utilizing an nlnfit function in Matlab to obtain the functional relation between the torque and the torsion angle;
(6) because of dT/d theta, the expression obtained by the above formula is subjected to derivation to obtain an expression of torsional rigidity of the harmonic gear.
CN202110911816.XA 2021-08-10 2021-08-10 Harmonic reducer transmission system torsional rigidity fitting method Pending CN113722891A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110911816.XA CN113722891A (en) 2021-08-10 2021-08-10 Harmonic reducer transmission system torsional rigidity fitting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110911816.XA CN113722891A (en) 2021-08-10 2021-08-10 Harmonic reducer transmission system torsional rigidity fitting method

Publications (1)

Publication Number Publication Date
CN113722891A true CN113722891A (en) 2021-11-30

Family

ID=78675308

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110911816.XA Pending CN113722891A (en) 2021-08-10 2021-08-10 Harmonic reducer transmission system torsional rigidity fitting method

Country Status (1)

Country Link
CN (1) CN113722891A (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110263492A (en) * 2019-07-15 2019-09-20 北京工业大学 A kind of harmonic speed reducer double circular arc tooth outline torsion stiffness calculation method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110263492A (en) * 2019-07-15 2019-09-20 北京工业大学 A kind of harmonic speed reducer double circular arc tooth outline torsion stiffness calculation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨聪彬等: "谐波传动柔软变形测量误差分析与补偿", 《光学精密工程》, vol. 29, no. 4, pages 793 - 801 *

Similar Documents

Publication Publication Date Title
CN104065322B (en) Method for controlling output feedback of motor position servo system
CN106845048B (en) Inner gearing gear shaft speed reducer nonlinear dynamics modeling method taking friction and tooth side clearance into account
CN104723340A (en) Impedance control method for flexibility joint mechanical arm based on connection and damping configuration
CN107944174B (en) Method for obtaining tooth direction load distribution coefficient of cylindrical gear
CN109766670B (en) Harmonic reducer reliability analysis method for industrial robot based on Copula function
Brauer Transmission error in anti-backlash conical involute gear transmissions: a global–local FE approach
CN103615501A (en) Small harmonic reducer and optimum design method thereof
Padmanabhan et al. Modeling of automotive gear rattle phenomenon: state of the art
Hu et al. Coupled translation-rotation vibration and dynamic analysis of face geared rotor system
CN110375054B (en) Asymmetric gear design method based on tooth profile inclination deviation
CN103092075A (en) Parameter determination method of machine driven system model with backlash characteristics
CN108107728B (en) Electro-hydraulic position servo system control method based on interference compensation
Liu et al. Nonlinear dynamic analysis of two external excitations for the gear system using an original computational algorithm
Gidado et al. Design, modeling and analysis of helical gear according bending strength using AGMA and ANSYS
CN112959321B (en) Robot flexible joint conversion error compensation method based on improved PI structure
CN113722891A (en) Harmonic reducer transmission system torsional rigidity fitting method
CN106227964A (en) Nonlinear systems with hysteresis parameter identification method based on extended state observer
CN109271655B (en) Material scale effect analysis method based on asymmetric finite element algorithm
CN116702549A (en) Method and system for calculating fatigue crack growth life of gear under consideration of random dynamic load
CN105784360A (en) Method for determining gear engagement dynamic stiffness based on engagement contact line length variation
Kawai et al. High back-drivable pseudo I-PD torque control using load-side torque observer with torsion torque sensor
CN102192303B (en) Method for calculating engagement noises of spiral bevel gear pair
Mote et al. Design and FEM analysis of helical gear
CN110263492B (en) Method for calculating torsional rigidity of double-arc tooth profile of harmonic reducer
Zhong et al. Low-delay admittance control of hydraulic series elastic actuator for safe human-robot collaboration

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