CN112883519A - Service life prediction method and device for cycloid wheel - Google Patents
Service life prediction method and device for cycloid wheel Download PDFInfo
- Publication number
- CN112883519A CN112883519A CN202110298067.8A CN202110298067A CN112883519A CN 112883519 A CN112883519 A CN 112883519A CN 202110298067 A CN202110298067 A CN 202110298067A CN 112883519 A CN112883519 A CN 112883519A
- Authority
- CN
- China
- Prior art keywords
- gear
- wheel
- cycloidal gear
- cycloidal
- tooth profile
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005299 abrasion Methods 0.000 claims abstract description 42
- 230000005540 biological transmission Effects 0.000 claims abstract description 42
- 230000004048 modification Effects 0.000 claims abstract description 20
- 238000012986 modification Methods 0.000 claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 230000001052 transient effect Effects 0.000 claims abstract description 8
- 230000015654 memory Effects 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 7
- 238000005461 lubrication Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 8
- 239000000463 material Substances 0.000 description 4
- 239000003638 chemical reducing agent Substances 0.000 description 3
- 238000009966 trimming Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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 relates to a service life prediction method and device for a cycloidal gear, and belongs to the technical field of service life prediction of cycloidal gears. The method comprises the steps of determining the abrasion loss of a cycloid wheel in a main working area at set time according to an Archard abrasion model and a Hertz contact theory through a modified tooth profile equation of the cycloid wheel, and taking the calculated abrasion loss as a first equidistant modification on a modified tooth profile to obtain a tooth profile equation after abrasion; and carrying out three-dimensional modeling and transient dynamics simulation on the worn tooth profile equation, substituting the corresponding rotation angle of the cycloidal gear and the pin wheel obtained by simulation into a calculation formula of the transmission error, calculating the transmission error of the cycloidal gear, and predicting the service life of the cycloidal gear according to the transmission error of the cycloidal gear. Through the process, the service life of the cycloid wheel can be accurately predicted according to the abrasion condition of the cycloid wheel.
Description
Technical Field
The invention relates to a service life prediction method and device for a cycloidal gear, and belongs to the technical field of service life prediction of cycloidal gears.
Background
The key component of the RV reducer is a cycloidal gear pin wheel planetary gear train, and the precision of a cycloidal gear in the cycloidal gear pin wheel planetary gear train directly determines the performance of the RV reducer.
Taking the geometric center of the cycloidal gear as an original point, selecting an axis which passes through the original point and is coincident with the symmetric axis of the cycloidal gear as a Y axis, and an axis which passes through the original point and is vertical to the Y axis as an X axis, wherein the tooth profile equation of the standard cycloidal gear is as follows:
in the formula: x represents the coordinates of the horizontal axis of the tooth profile of the cycloidal gear; y represents the coordinates of the vertical axis of the tooth profile of the cycloidal gear; rzRepresenting the radius of a pin wheel distribution circle;representing the rotation angle parameter of the cycloidal gear; a represents eccentricity; zaRepresenting the number of cycloidal gear teeth; zbIndicating the number of pin gear teeth; r iszRepresenting the pinwheel radius; k1Representing short-amplitude coefficients.
The standard cycloidal gear pin wheel transmission theoretically belongs to zero-clearance meshing, and the number of teeth participating in meshing is half of that of the pin wheel. Because the cycloid wheel is required to be shaped, the parts have processing errors and assembly errors, the cycloid wheel and the pin wheel can also generate elastic deformation after the gear train is loaded, transmission errors can exist in actual transmission, the transmission errors can be increased due to abrasion of the cycloid wheel, the accuracy of the RV reducer cannot meet the requirements due to overlarge transmission errors of the cycloid wheel, and the service life of the cycloid wheel is further influenced.
Disclosure of Invention
The invention aims to provide a service life prediction method and a service life prediction device for a cycloidal gear, so as to realize detection of the service life of the cycloidal gear.
The invention provides a service life prediction method of a cycloid wheel for solving the technical problems, which comprises the following steps:
1) determining the wear amount of the cycloidal gear in a main working area under set time according to an Archard wear model and a Hertz contact theory based on a modified tooth profile equation of the cycloidal gear, wherein the main working area of the cycloidal gear refers to a position between a tooth root and a tooth top of a gear tooth, and the position is close to the middle;
2) the calculated abrasion loss is regarded as that the modified tooth profile is modified once again in an equidistant manner, and a tooth profile equation after abrasion is determined;
3) carrying out three-dimensional modeling on the worn tooth profile equation, carrying out transient dynamics simulation to obtain the rotation angle of the cycloidal gear and the rotation angle of the pin wheel at each moment, and calculating a transmission error according to the rotation angle of the pin wheel and the rotation angle of the cycloidal gear at each moment and the rotation angle of the pin wheel and the rotation angle of the cycloidal gear when the pin wheel and the cycloidal gear are engaged at the engagement reference point;
4) and judging whether the transmission error meets the precision requirement, if so, increasing the set time to recalculate the abrasion loss, calculating the transmission error according to the simulation result of the tooth profile equation established by the new abrasion loss until the transmission error does not meet the precision requirement, and taking the set time for calculating the abrasion loss at the last time as the service life of the cycloid wheel.
The invention also provides a service life prediction device of the cycloidal gear, which comprises a processor and a memory, wherein the processor executes a computer program stored by the memory so as to realize the service life prediction method of the cycloidal gear.
The method comprises the steps of calculating the abrasion loss of the cycloidal gear in a main working area, and taking the calculated abrasion loss as a first equidistant modification on a modified tooth profile, so as to obtain a tooth profile equation after abrasion; and carrying out three-dimensional modeling and transient dynamics simulation on the worn tooth profile equation, determining a calculation formula of the transmission error, calculating the transmission error of the cycloidal gear according to the calculation formula, and predicting the service life of the cycloidal gear according to the transmission error of the cycloidal gear. Through the process, the service life of the cycloid wheel can be accurately predicted according to the abrasion condition of the cycloid wheel.
Further, the wear amount of the cycloid wheel in the step 1) is as follows:
hI=2a1λntεαIh
in the formula, the amount of wear hIExpressed in terms of wear depth at the point of engagement, in m, a1Contact half width, lambda is the lubrication state, n is the rotational speed, t is the run time, epsilonαIs degree of coincidence, IhIs the rate of wear.
Further, the worn tooth profile equation determined in step 2) is:
in the formulak″=azp/(rp+Δrp+ Δ h), wherein zcIndicating the number of teeth of the cycloid gear, a indicating the eccentricity, rpRepresenting the radius of the centre circle of the pin teeth, rrpRepresenting the radius of the teeth, Δ rrp、ΔrpAnd Δ θ is a profile modification parameter and Δ h is a wear amount function.
Further, the calculation formula of the transmission error is as follows:
in the formula10、φ20Respectively the pin wheel rotation angle and the cycloid wheel rotation angle phi when engaged at the engagement reference point1、φ2Respectively the angle of the pin wheel and the angle of the cycloidal gear at each moment, zpNumber of teeth of pin gear, zcThe number of teeth of the cycloid gear is shown.
Drawings
FIG. 1 is a flow chart of a method of predicting the life of a cycloidal gear of the present invention;
FIG. 2 is a schematic illustration of a gear tooth engagement in an embodiment of the present invention;
FIG. 3 is a schematic illustration of a gear tooth engagement contact pressure distribution in an embodiment of the present invention;
fig. 4 is a block diagram showing the structure of a service life predicting apparatus for a cycloid wheel according to the present invention.
Detailed Description
The following further describes embodiments of the present invention with reference to the drawings.
The method is based on the Hertz contact theory, the Archard wear theory and the differential geometry principle, the wear loss of the cycloidal gear in a main working area is calculated, a tooth profile equation of the cycloidal gear after the cycloidal gear works for a plurality of time periods under a specific working condition is obtained, transient dynamics simulation is carried out on a worn transmission pair through finite element software, a corresponding rotation angle relation of two wheels is obtained, the simulated rotation angles of the cycloidal gear and a pin wheel are brought into a calculation formula of a transmission error, the transmission error is calculated, and the service life is predicted according to the transmission error. The specific implementation process of the method is shown in fig. 1, and the specific implementation steps are as follows.
1. And determining the abrasion loss of the cycloidal gear in a main working area.
The primary working area of the gear teeth of the cycloidal gear is located near the middle between the tooth root and the tooth tip of the gear teeth, while the tooth root and the tooth tip are non-primary working areas. When two wheels of the cycloidal gear are engaged, wheel 1 is at contact point j, as shown in figure 2(1)'At j with wheel 2(2)'And (4) meshing. After dt times, the wheel 1 is at contact point j(1)”At contact point j with wheel 2(2)”Engagement, during dt times, the normal wear depth of wheel 1, wheel 2 at j point is dh respectively(1)And dh(2)The volume dropped by abrasion is dV(1)And dV(2)The friction distance is ds.
As can be seen from the archer wear formula, wheel 1 can be expressed as:
if the contact area at point j is dA and the pressure at the point of contact is p, the above equation can be rewritten as follows:
where K is the coefficient related to the wear of the material, and the hardness of the material is constant, the above equation can be simplified as follows:
in the formula IhWear rate, the depth of wear per unit friction distance; k is the wear coefficient; p is the contact pressure, and the integral of the wear amount h (normal wear depth at the meshing point) is expressed as:
h=∫kpds
for a cycloidal pin gear transmission pair, a pin gear is driven, a cycloidal gear is driven, and the sliding coefficient of the driving and driven wheels at the meshing point can be expressed as follows:
wherein v ist1、vt2Is the tangential speed of the driving and driven wheels at the point of engagement.
In an unmodified gear, the slip speed is usually solved by converting it into an expression relating to the gear ratio and the pressure angle, using the relationship between the tangential speed of the mesh point, the gear ratio and the pressure angle. However, for the cycloidal pin gear transmission pair, in actual use, the cycloidal gear needs to be modified to meet the working condition requirement, and modification generates tooth side clearance, so that the instantaneous transmission ratio is not constant. Therefore, in the invention, the influence generated by modification is considered, the method for accurately solving the meshing point is adopted, and the friction distance ds is directly used for solving the sliding coefficient.
The standard tooth profile equation of the cycloid wheel can be expressed as:
K1=azp/rp iH=zp/zc
in the formula, rpRepresenting the radius of the central circle of the needle teeth; r isrpRepresenting the needle tooth radius; a represents eccentricity; z is a radical ofpIndicating the number of pin teeth; k1Representing the short-amplitude coefficient; i.e. iHRepresents a gear ratio; z is a radical ofcRepresenting the number of cycloidal gear teeth; phi denotes the meshing phase angle.
For the trimming cycloid wheel, the trimming amount is directly superposed on a standard tooth profile equation, and the tooth profile equation after trimming is obtained as follows:
in the formulaWherein z isc、a、rp、rrpAs a basic parameter,. DELTA.rrp、ΔrpAnd Δ θ is a shape modification parameter. (these three shape modification parameters correspond to three shape modification methods of equidistance, displacement and corner respectively, and can be combined for use, also can be used singly, all have the value that is the combined shape modification, one has the value, and the other two take 0, is the single shape modification).
Expressing the tooth profile equations of the cycloidal gear and the pin wheel under a fixed coordinate system, establishing a TCA model, and enabling the position vector and the normal vector of the TCA model to be equal to each other, so that the accurate coordinates of each meshing point on the tooth profile of the cycloidal gear and the corresponding meshing point on the pin wheel can be obtained. Dispersing each meshing point on the cycloid wheel with equal arc length, and dividing any meshing point j(2)Two adjacent left and rightMesh point j-1(2)And j +1(2)The coordinates are used as the starting point and the end point of integration, the curve length is obtained by performing curve integration on the tooth profile of the cycloidal gear, and the curve length is used as the sliding distance ds of the cycloidal gear at the meshing point in dt time(2)Two corresponding meshing points j-1 on the pinwheel(1)And j +1(1)The length of minor arc enclosed on the gear tooth profile is used as the contact point j of the gear in dt times(1)Sliding distance ds of(1)Substituting the sliding distance to obtain the sliding coefficient lambda of the driving wheel and the driven wheel1And λ2。
In the transmission of the cycloid pin wheel, the size of the meshing force borne by each meshing point of the cycloid wheel is in direct proportion to each deformation quantity. Let the load torque be TcThe initial value of the maximum contact force of the pin teeth is as follows:
then, F is solved through successive iterationmaxThe distance l from the common normal line of the maximum stress point to the center of the cycloid wheelmax=azc。
The cycloidal pin gear transmission belongs to a typical Hertz contact, and the deformation of a maximum contact point under stress can be obtained by a Hertz formula as follows:
in the formulaρ represents a cycloid curvature radius; b represents the effective tooth width of the cycloid wheel; e1、E2Respectively representing the elastic modulus of materials of a cycloid wheel and a pinwheel; mu.s1、μ2Representing the poisson's ratio of the cycloidal and pinwheel materials, respectively.
Total deformation delta of each meshing point of cycloid wheel in normal directioniCan be expressed as
According to the above formula, the amount of deformation of the meshing point can be compared with the backlash Δ s (φ)i) Judging whether the needle teeth are meshed with the cycloid gear teeth or not according to the size relation of the gear teeth, and further determining the number of simultaneously meshed teeth. Assuming that the teeth numbered m to n are simultaneously engaged, the maximum contact force can be expressed as:
comparing the result of the above formula with the initially assumed maximum value, and if the values are equal, ending the solution; if not, F is requiredmaxIs given as Fmax0Repeating the iterative calculation until the two formulas have the same result, and outputting Fmax。
The normal engagement force applied to any engagement point can be expressed as:
the cycloidal pin gear transmission belongs to a typical Hertz contact, two tooth profiles are meshed and can be equivalently contacted with two cylinders at a meshing point, the pressure distribution is shown in figure 3, and the contact half width is as follows:
wherein W represents the tooth profile contact point normal pressure; r*Represents the equivalent radius; e*Denotes the equivalent elastic modulus, and has
W=Fi
In the formula, R1、R2The radii of curvature of the needle wheel and the cycloid wheel at the contact point are respectively shown, and v1 and v2 respectively show the poisson ratio of the cycloid wheel and the needle wheel.
The contact half width obtained by the above process is
The contact pressure of the meshing point is
In the formula, b represents an effective tooth width of the cycloid wheel.
The wear coefficient can be expressed as
The minimum oil film thickness of the cycloid pinwheel pair at the contact point is expressed as:
wherein λ ═ HminAnd/σ, σ represents the square root of the sum of the squares of the surface roughness of the two contact surfaces. The values of lambda are different, which represents that the lubrication states of the gears are different, and 3 ranges from small to large respectively represent 3 states of boundary lubrication, mixed lubrication and good lubrication. u. of1、u2Representing the scrolling speed of two wheels; eta is viscosity coefficient.
One point on the cycloid wheel participates in meshing for 1 time, and the generated friction distance is as follows:
S=2a1λ2
the total friction distance during operation is:
X=Sntεα
wherein t represents the operating time; n represents a rotation speed; epsilonαIndicating the degree of coincidence.
The amount of wear at the engagement point was:
hI=XIh
the amount of wear of the cycloid gear that can be obtained from the above process can be expressed as:
hI=2a1λntεαIh
in the formula, the wear amount is expressed in the unit of m in terms of the wear depth at the meshing point.
2. And establishing a tooth profile equation of the trimmed cycloid wheel after being worn based on the obtained wear loss.
The invention expresses the abrasion loss by the abrasion depth at the meshing point, so the abrasion loss is actually a depth which is overlapped inwards along the normal direction at the tooth profile point, the equidistant modification is also overlapped by a distance along the normal direction, the distance is a fixed value, the abrasion loss can be changed at different meshing points, the abrasion loss can be regarded as the equidistant modification, and the tooth profile equation after abrasion can be obtained by overlapping the equidistant modification in the same way according to the expression method of the equidistant modification in the tooth profile equation.
And (4) taking the abrasion loss as the modification tooth profile, and then performing equidistant modification once to obtain the tooth profile equation after abrasion. Change in the corresponding equation to be Δ rpBecomes Δ rp+ Δ h, the tooth profile equation after the traditional profile modification cycloid wheel is worn can be obtained as follows:
in the formulaWherein z isc、a、rp、rrpAs a basic parameter,. DELTA.rrp、ΔrpAnd delta theta respectively represent an equidistant modification parameter, a displacement modification parameter and a corner modification parameter, delta h is a function of the abrasion loss, and delta h and hIIt is the same, essentially the same, in order to be written in incremental form uniformly with other parameters.
3. Transient dynamics simulation is carried out on the worn transmission pair through finite element software, the corresponding corner relation of the two wheels is obtained, and the transmission error is calculated.
The transmission error formula is only related to a pin wheel corner, a cycloid wheel corner, a pin wheel tooth number, a cycloid wheel tooth number, a pin wheel corner and a cycloid wheel corner when meshing is carried out at a meshing reference point, however, the pin wheel corner and the cycloid wheel corner change along with the abrasion process, the tooth profile changes after abrasion, the corresponding relation between the pin wheel corner and the cycloid wheel corner also changes, transient dynamics simulation is carried out on the abraded tooth profile by utilizing finite element software to obtain the cycloid wheel corner and the pin wheel corner at any meshing position when the abraded cycloid wheel is transmitted, and the transmission error of any meshing position after abrasion can be obtained according to the transmission error formula.
Carrying out three-dimensional modeling on the established tooth profile equation after abrasion, carrying out transient dynamics simulation in finite element simulation software, and obtaining the rotating angle of the cycloid wheel and the pin wheel at each instant, so as to calculate the transmission error, wherein the calculation formula of the transmission error is as follows:
wherein phi10、φ20Respectively the pin wheel rotation angle and the cycloid wheel rotation angle phi when engaged at the engagement reference point1、φ2Respectively the angle of the pinwheel and the angle of the cycloidal gear, zpNumber of teeth of pin gear, zcThe number of teeth of the cycloid gear is shown.
4. And predicting the service life according to the transmission error.
The calculation of the abrasion loss is related to time, the longer the time is, the larger the abrasion loss is, the abrasion loss can change a tooth profile equation, when two wheels are meshed, the rotation angle of the cycloidal gear and the rotation angle of the pin wheel can both change, if the obtained transmission error meets the requirement, the tooth profile change caused by abrasion generated under the set time of the set condition is not enough to lose the precision, the time is increased, the abrasion loss is calculated again, the tooth profile equation after abrasion is obtained, the rotation angle is obtained through simulation again, the transmission error is calculated, the step is repeated until the obtained transmission error is larger than the set precision, the requirement is not met, the time used for calculating the abrasion loss at the last time is the time for ensuring the precision of the cycloidal gear, and the service life of the cycloidal gear is the service life of the cyclo.
Device embodiment
The apparatus proposed in this embodiment, as shown in fig. 4, includes a processor and a memory, where a computer program operable on the processor is stored in the memory, and the processor implements the method of the foregoing method embodiment when executing the computer program. That is, the method in the above method embodiment should be understood that the flow of the cycloidal wheel life prediction method may be implemented by computer program instructions. These computer program instructions may be provided to a processor such that execution of the instructions by the processor results in the implementation of the functions specified in the method flow described above.
The processor referred to in this embodiment refers to a processing device such as a microprocessor MCU or a programmable logic device FPGA; the memory referred to in this embodiment includes a physical device for storing information, and generally, information is digitized and then stored in a medium using an electric, magnetic, optical, or the like. For example: various memories for storing information by using an electric energy mode, such as RAM, ROM and the like; various memories for storing information by magnetic energy, such as hard disk, floppy disk, magnetic tape, magnetic core memory, bubble memory, and U disk; various types of memory, CD or DVD, that store information optically. Of course, there are other ways of memory, such as quantum memory, graphene memory, and so forth.
The apparatus comprising the memory, the processor and the computer program is realized by the processor executing corresponding program instructions in the computer, and the processor can be loaded with various operating systems, such as windows operating system, linux system, android, iOS system, and the like. As other embodiments, the device can also comprise a display, and the display is used for displaying the diagnosis result for the reference of workers.
Claims (5)
1. A service life prediction method of a cycloidal gear is characterized by comprising the following steps:
1) determining the wear amount of the cycloidal gear in a main working area under set time according to an Archard wear model and a Hertz contact theory based on a modified tooth profile equation of the cycloidal gear, wherein the main working area of the cycloidal gear refers to a position between a tooth root and a tooth top of a gear tooth, and the position is close to the middle;
2) the calculated abrasion loss is regarded as that the modified tooth profile is modified once again in an equidistant manner, and a tooth profile equation after abrasion is determined;
3) carrying out three-dimensional modeling on the worn tooth profile equation, carrying out transient dynamics simulation to obtain the rotation angle of the cycloidal gear and the rotation angle of the pin wheel at each moment, and calculating a transmission error according to the rotation angle of the pin wheel and the rotation angle of the cycloidal gear at each moment and the rotation angle of the pin wheel and the rotation angle of the cycloidal gear when the pin wheel and the cycloidal gear are engaged at the engagement reference point;
4) and judging whether the transmission error meets the precision requirement, if so, increasing the set time to recalculate the abrasion loss, calculating the transmission error according to the simulation result of the tooth profile equation established by the new abrasion loss until the transmission error does not meet the precision requirement, and taking the set time for calculating the abrasion loss at the last time as the service life of the cycloid wheel.
2. The method for predicting the life of a cycloid gear of claim 1 wherein the amount of wear of the cycloid gear in step 1) is:
hI=2a1λntεαIh
in the formula, the amount of wear hIExpressed in terms of wear depth at the point of engagement, in m, a1Contact half width, lambda is the lubrication state, n is the rotational speed, t is the run time, epsilonαIs degree of coincidence, IhIs the rate of wear.
3. The method for predicting the life of a cycloidal gear according to claim 1 or 2, wherein the worn tooth profile equation determined in step 2) is:
in the formulak″=azp/(rp+Δrp+ Δ h), wherein zcIndicating the number of teeth of the cycloid gear, a indicating the eccentricity, rpRepresenting the radius of the centre circle of the pin teeth, rrpRepresenting the radius of the teeth, Δ rrp、ΔrpAnd Δ θ is a profile modification parameter and Δ h is a wear amount function.
4. The method for predicting the life of a cycloidal gear according to claim 3, wherein the transmission error is calculated by the formula:
in the formula10、φ20Respectively the pin wheel rotation angle and the cycloid wheel rotation angle phi when engaged at the engagement reference point1、φ2Respectively the angle of the pin wheel and the angle of the cycloidal gear at each moment, zpNumber of teeth of pin gear, zcThe number of teeth of the cycloid gear is shown.
5. A cycloidal wheel life prediction device, characterized in that the device comprises a processor and a memory, the processor executing a computer program stored by the memory to implement the cycloidal wheel life prediction method according to any one of the preceding claims 1-4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110298067.8A CN112883519B (en) | 2021-03-19 | 2021-03-19 | Service life prediction method and device for cycloid wheel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110298067.8A CN112883519B (en) | 2021-03-19 | 2021-03-19 | Service life prediction method and device for cycloid wheel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112883519A true CN112883519A (en) | 2021-06-01 |
CN112883519B CN112883519B (en) | 2023-03-14 |
Family
ID=76041460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110298067.8A Active CN112883519B (en) | 2021-03-19 | 2021-03-19 | Service life prediction method and device for cycloid wheel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112883519B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013088121A (en) * | 2011-10-13 | 2013-05-13 | Ricoh Co Ltd | Apparatus with drive transmission mechanism |
CN105701348A (en) * | 2016-01-12 | 2016-06-22 | 芜湖赛宝信息产业技术研究院有限公司 | Method for predicting reliable life of harmonic reducer |
WO2018103160A1 (en) * | 2016-12-09 | 2018-06-14 | 深圳市荣德机器人科技有限公司 | Cycloidal-pin wheel speed reducer |
CN108730480A (en) * | 2018-07-27 | 2018-11-02 | 河南科技大学 | A kind of Cycloidal Wheel and its reverse active correction method of flank profil of RV retarders |
WO2018210350A1 (en) * | 2017-05-19 | 2018-11-22 | 北京工业大学 | Test apparatus for detecting rv reducer cycloidal wheel and needle bearing material abrasion |
CN111062103A (en) * | 2019-12-24 | 2020-04-24 | 安徽工业大学 | Metallurgical crane reducer gear service life assessment method with optimized parameters |
CN111059255A (en) * | 2019-12-27 | 2020-04-24 | 太原理工大学 | Calculation method for tooth surface abrasion of double-arc harmonic reducer |
CN112163325A (en) * | 2020-09-10 | 2021-01-01 | 浙江工业大学 | RV reducer service life prediction method based on digital twinning |
-
2021
- 2021-03-19 CN CN202110298067.8A patent/CN112883519B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013088121A (en) * | 2011-10-13 | 2013-05-13 | Ricoh Co Ltd | Apparatus with drive transmission mechanism |
CN105701348A (en) * | 2016-01-12 | 2016-06-22 | 芜湖赛宝信息产业技术研究院有限公司 | Method for predicting reliable life of harmonic reducer |
WO2018103160A1 (en) * | 2016-12-09 | 2018-06-14 | 深圳市荣德机器人科技有限公司 | Cycloidal-pin wheel speed reducer |
WO2018210350A1 (en) * | 2017-05-19 | 2018-11-22 | 北京工业大学 | Test apparatus for detecting rv reducer cycloidal wheel and needle bearing material abrasion |
CN108730480A (en) * | 2018-07-27 | 2018-11-02 | 河南科技大学 | A kind of Cycloidal Wheel and its reverse active correction method of flank profil of RV retarders |
CN111062103A (en) * | 2019-12-24 | 2020-04-24 | 安徽工业大学 | Metallurgical crane reducer gear service life assessment method with optimized parameters |
CN111059255A (en) * | 2019-12-27 | 2020-04-24 | 太原理工大学 | Calculation method for tooth surface abrasion of double-arc harmonic reducer |
CN112163325A (en) * | 2020-09-10 | 2021-01-01 | 浙江工业大学 | RV reducer service life prediction method based on digital twinning |
Non-Patent Citations (2)
Title |
---|
HAN LINSHAN等: "The Optimal Design of the 2K-V-type Cycloid-pin Reducer", 《2012 THIRD INTERNATIONAL CONFERENCE ON DIGITAL MANUFACTURING & AUTOMATION》 * |
张红等: "摆线针轮传动接触非线性疲劳寿命分析", 《国外电子测量技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112883519B (en) | 2023-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marafona et al. | Mesh stiffness models for cylindrical gears: A detailed review | |
US7599824B2 (en) | Gear cutting simulation method, gear cutting simulation program, and gear cutting simulation device | |
Jyothirmai et al. | A finite element approach to bending, contact and fatigue stress distribution in helical gear systems | |
CN102626809B (en) | Method for determining geometrical data of bevel pinion and bevel pinion drive with the bevel gear | |
Wang et al. | Simulating coupling behavior of spur gear meshing and fatigue crack propagation in tooth root | |
Chen et al. | Contact stress analysis of concave conical involute gear pairs with non-parallel axes | |
Guingand et al. | Analysis and optimization of the loaded meshing of face gears | |
Zschippang et al. | Face-gear drive: Assessment of load sharing, transmission characteristics and root stress based on a quasi-static analysis | |
Liu et al. | Investigation on contact and bending stress of face-hobbed and face-milled hypoid gear | |
CN112883519B (en) | Service life prediction method and device for cycloid wheel | |
Vishwakarma et al. | Finite element analysis of helical gear using three-dimensional cad model | |
CN109446710B (en) | Numerical calculation method for load transmission error of spiral bevel gear | |
Bair et al. | Mathematical model and characteristic analysis of elliptical gears manufactured by circular-arc shaper cutters | |
Li et al. | A Loaded Analysis Method for RV Cycloidal-pin Transmission Based on the Minimum Energy Principle. | |
Patel | Finite element analysis of stresses in involute spur & helical gear | |
Li et al. | Time-varying mesh stiffness calculation of spiral bevel gear with spalling defect | |
Xiang et al. | The meshing angular velocity and tangential contact force simulation for logarithmic spiral bevel gear based on Hertz elastic contact theory | |
CN113010978A (en) | Aviation straight gear shaping method based on dynamic simulation | |
Aziz et al. | Knowledge-based geometry generation for spur and helical gears | |
Hiltcher et al. | Numerical simulation and optimisation of worm gear cutting | |
Nachimowicz et al. | Modelling the meshing of cycloidal gears | |
CN113868801A (en) | Gear tooth contact analysis method based on full tooth surface search | |
Sell | Finite element modeling spur and helical gears in contact | |
Adnan et al. | Stress Analysis Validation for Gear Design | |
Shaker | Optimization of tooth-root profile for maximum load-carrying capacity: spur and bevel gears |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210601 Assignee: Ruihong Machinery Processing Factory in Xigong District, Luoyang City Assignor: HENAN University OF SCIENCE AND TECHNOLOGY Contract record no.: X2024980001395 Denomination of invention: A method and device for predicting the lifespan of cycloidal gears Granted publication date: 20230314 License type: Exclusive License Record date: 20240125 |
|
EE01 | Entry into force of recordation of patent licensing contract |