CN103206515A - Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors - Google Patents

Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors Download PDF

Info

Publication number
CN103206515A
CN103206515A CN2013101369269A CN201310136926A CN103206515A CN 103206515 A CN103206515 A CN 103206515A CN 2013101369269 A CN2013101369269 A CN 2013101369269A CN 201310136926 A CN201310136926 A CN 201310136926A CN 103206515 A CN103206515 A CN 103206515A
Authority
CN
China
Prior art keywords
delta
error
gear
rightarrow
tooth
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
CN2013101369269A
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 CN2013101369269A priority Critical patent/CN103206515A/en
Publication of CN103206515A publication Critical patent/CN103206515A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention relates to a loaded tooth surface contact analysis method direct at epicycloid bevel gear errors, and belongs to the field of theory of non linear vibration. Gear pair elastic deformation, support system torsion, gear processing errors and machine tool setting errors are taken into consideration to derive a novel error loaded contact analysis method. When analyzing one loaded tooth surfaced contact, real situation can be closer. The loaded tooth surface contact analysis method is usually used for setting parameters and machine tool parameters via V-H-J (V is the longitudinal set parameter of a gear pair, and the upward direction of a small gear is active; H is a coordinate set parameter of the gear pair, and the direction facing towards the big end of the small gear is active; J is the adjusting parameter along a coordinate direction, and the direction, distant from a big gear, of the small gear is active) within a real contact area under the conditions of loading and delivering errors and simulated load distribution and processing, so that impact of the errors is compensated to improve delivering performance and correct size and position of the contact area.

Description

A kind of method that loads Tooth Contact Analysis at the overlikon spiral bevel gear error
Technical field
The present invention relates to a kind of method at overlikon spiral bevel gear error loading Tooth Contact Analysis, belong to gear nonlinear dynamics field.
Background technique
Crin Gen Beierge bevel gear (Klingelnberg Spiral Bevel Gear) is as one of two canine tooth systems of helical bevel gear, have characteristics such as stable drive, bearing capacity height, hard surface skiving technology, thereby being specially adapted to high-power and high pulling torque heavy load transmission field, is the core transmission component in the key areas such as heavy high-grade, digitally controlled machine tools, car transmissions, Aero-Space equipment.Steadily reliable in transmission process in order to guarantee gear, reduce vibration and noise, must carry out dynamic analysis to gear train assembly.Effectively setting up gear engagement mathematical model is the basis of dynamic analysis.
Domestic and international many scholars have carried out more extensive and deep research to gear teeth face contact analysis dynamics, but majority all is the Tooth Contact Analysis of considering only to comprise the lathe step-up error, and concentrate on straight toothed spur gear, and the loading flank of tooth Contact Analysis and Study of considering transmission error, shaft distortion, bearing deformation, tooth distortion, flank of tooth machining error and lathe step-up error at crin Gen Beierge bevel gear seldom.Truth can be pressed close to more when loading Tooth Contact Analysis to one when considering these sum of errors amount of deformation, and the size and location of transmitting performance and revising area of contact can be improved.
Summary of the invention
The purpose of this invention is to provide a kind of method that loads Tooth Contact Analysis at the overlikon spiral bevel gear error, this method can apply in the loading Tooth Contact Analysis of crin Gen Beierge bevel gear error, when loading Tooth Contact Analysis to one, can press close to truth more, this method can be improved the transmission performance of gear and revise the area of contact size simultaneously, also provides basis and foundation for overlikon spiral bevel gear vibration control, noise reduction and defects detection.
The present invention adopts following technological means to realize:
A kind of method at overlikon spiral bevel gear error loading Tooth Contact Analysis is characterized in that this method may further comprise the steps:
1) distortion and transmission error
1.1) transmission error under loading environment
Transmission error under the loading environment is made up of three parts: position deviation error delta E a, the transmission error Δ E that the gear resiliently deformable causes bThe transmission error Δ E that causes with the bending deflection of gear c, therefore, comprehensive transmission error Δ E is
ΔE=ΔE a+ΔE b+ΔE c (1)
1.2) load distribution and compatibility of deformation rule
Suppose to have the tooth of n to mesh, can obtain the load distribution relation
T = Σ i - 1 n T i = T 1 + T 2 + . . . + T n T i = 0.5 d i · F i · cos α i - - - ( 2 )
In the formula, T is resultant couple, T iBe the moment of the gear pair of i, d iBe the gear diameter of point of contact, F iBe the normal force of point of contact, α iLoading direction and tangential angle;
The right transmission deviation of each engaging tooth is identical, therefore, releases the gear transmission formula
Δ E 1 = Δ E 2 = . . . Δ E n Δ E i = Δ E ia = Δ E ib = Δ E ic - - - ( 3 )
In the formula, Δ E iIt is the transmission error of the gear pair of i;
1.3) axle distortion
The axle of gear reverses and is
Figure BDA00003071680300023
The bending deflection of gear is f Xr, f YrAnd f Zr, the axle relevant parameter of small gear is
Figure BDA00003071680300024
f Xl, f Yl, and f Zl
The supporting structure of overlikon spiral bevel gear can be reduced to an overhang or a support beam, uses Fi and Ti can calculate bending deflection and torsional deflection, and can calculate the displacement of point of contact
Figure BDA00003071680300025
f → ai = ( f xr , f yr , f zr ) - - - ( 4 )
1.4) bearing deformation
Consider supporting structure, Fi produces one
Figure BDA00003071680300027
Power because the effect of axial force, the axial displacement of small gear is f z, because effect and the bearing radial internal clearance of tangential force and radial force, be f along the displacement of the axial bearing centre of y y, because effect and the bearing axial internal clearance of axial force and radial force, be f along the displacement of the axial bearing centre of x x, therefore, the flank of tooth displacement that is produced by bearing deformation is
f → bi = ( f x , f y , f z ) f x = sign ( R zi ) ( δ z + u a ) f y = sign ( R yi ) ( δ y + u r / 2 ) cos ϵ f z = sign ( R xi ) ( δ x + u r / 2 ) sin ϵ - - - ( 5 )
Formula (5), sign (R) are sign functions; u a, u rBe respectively radial internal clearance and the axial internal clearance of bearing; ε is radially opposite force and the axially angle of opposite force; δ x, δ y, δ yBe axial deformation and the radial deformation that is caused by axial force and radial force effect;
1.5) the tooth distortion
1.5.1) be out of shape at the tooth at contact point of gear surface place
According to the Hertz formula, the tooth distortion at the contact point of gear surface place can expand into
δ c = λ 3 9 128 AF 2 ( 1 - μ 1 2 E 1 + 1 - μ 2 2 E 2 ) 2 - - - ( 6 )
In the formula, λ is a coefficient in the elastomechanics; E 1, E 2, μ 1, μ 2Be respectively Young's modulus and Poisson's ratio; A is the relative mean value of curvature; F is contact force;
1.5.2) tooth bending deflection
Calculate the bending deflection of tooth with the overhang method of Westinghouse, at the bending deflection value δ of contact point of gear surface wFor
δ w = FL 3 3 EI { 1 + 1.3 ( t / L ) + [ 0.25 + 0.75 ( 1 - μ ) ] · ( t / L ) 2 + 0.35 ( t / L ) 3 } - 1 - - - ( 7 )
In the formula, I is the quadrature moment of inertia in load cross section; T is the normal tooth thickness in load cross section; L is the tooth depth of load(ing) point; E is elasticity modulus of materials; μ is the material Poisson's ratio;
2) flank of tooth machining error
Consider tooth thickness error and tangential resultant error, the actual effect error is
δ jgr ′ = ( T → gr + k 1 · Δ F ir ′ · R r / r 1 ) · T → gr · u → r T → gr δ kgr ′ = ( T → gr + k 1 · ΔF ir ′ · R r / r 1 ) · T → gr · k → gr T → gr δ igr ′ = ( T → gr + k 1 · ΔE ir ′ · R r / r 1 ) · T → gr · i → gr T → gr - - - ( 8 )
In the formula,
Figure BDA00003071680300034
It is the tooth thickness error of gear;
Figure BDA00003071680300035
Be respectively axially, radially with tangential direction on unit an amount of; R rThe radially radius of tooth; Δ F ' IrIt is the tangential resultant error of gear; r 1It is the Pitch radius of gear; k 1Be conversion coefficient, coefficient k 1Be to use apportionment ratio, the relative compensating effect of tooth thickness error and tangential resultant error are determined, k 1Definition be associated with apportionment ratio and the total number of teeth in engagement of the actual tangential resultant error that meets normal distribution law;
3) lathe step-up error
The composition that lathe arranges parameter is: the error delta τ of cutter spacing setting angle τ; By cutterhead head center to carriage center M dThe lathe that produces of distance apart from error delta M dWith the error delta r in cutterhead head radius r;
3.1) cutter spacing sets angle τ and helixangle mBetween relation
τ=2sin -1(M d/ 2k), wherein k cutter range setting can be obtained setting the lathe of angular error generation apart from error delta M because of cutter spacing by this relation d
Δ M d = k cos τ 2 Δτ - - - ( 9 )
Cutter spacing sets angular error
Δτ = τ ′ = - R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) k M d Z p cos v - - - ( 10 )
In the formula, cos β m=M dZ p/ (2R), v=sin -1(0.5M dZ 0/ r) and R be pitch cone radius, R=R mAnd β=β m, lathe is apart from M d, helixangle m, cutterhead head radius r and pitch cone radius R; Cutterhead head number is Z 0, the crown gear number is Z p
3.2) cutterhead head radius r and helixangle mBetween relation
Be generally used for adjusting the error delta r of the blade pad of eccentric cutter and the combined error definition cutting end radius that slip is counted, this error can change the Δ β in the helix angle m, while R mAnd M dRemain unchanged,
Δτ = τ ′ = r R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) r 2 Z p cos v - r Z p R m sin ( β m - v ) cos v - R m Z 0 cos ( β m - v ) cos β m - - - ( 11 )
Following formula (11) has provided cutterhead head radius r and helixangle mBetween relation;
3.3) relation between lathe step-up error and the area of contact position
According to tooth surface equation and Tooth Contact Analysis equation, by the importing lathe parameter error relevant with them is set and sets up a new meshing condition equation
Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) - - - ( 12 )
In the formula, the three dimensional space distance vector of conjugation point of contact Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) , Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) Be to be respectively the left and right vector that revolves the conjugation point of contact in the processing engagement on the member flank of tooth; Δ r is the error of cutterhead head radius.
4) flank of tooth contact error analysis
The gear teeth face contact analysis is based on the rigidity of gear pair and support system, when considering the distortion of support system, the relative position of mesh tooth face will change in the space, this means the support system distortion and establishing error by the lathe that the lathe value of setting produces is equivalent in magnitude, according to formula (4), (5), (12) can obtain
Figure BDA00003071680300049
In the formula, α is the blade pressure angle; R rAnd R lIt is respectively the radius of gearwheel and pinion point of contact; δ r, δ lIt is respectively the generation angle of gearwheel and small gear; V be gear pair parameter vertically is set, and the upward direction of small gear is initiatively; H is that the coordinate axes of gear pair arranges parameter, and is initiatively towards the big extreme direction of small gear; J is the adjustment parameter along change in coordinate axis direction, and small gear is initiatively away from the direction of gearwheel.
Beneficial effect:
Characteristics of the present invention are only to have considered the lathe step-up error in the traditional loading Tooth Contact Analysis method of gear, and this method is at overlikon spiral bevel gear and consider the gear pair resiliently deformable, support system is reversed, Errors in Gear Processing and lathe step-up error have proposed a kind of new error and have loaded the contact analysis method.Can press close to truth more when being added in Tooth Contact Analysis by this method to one.The inventive method can apply to crin Gen Beierge bevel gear be used for to solve loading, transmission error, load distribute and machining simulation under real contact area by V-H-J parameter is set and lathe arranges parameter, the influence of compensating error improves the size and location of transmitting performance and revising area of contact.Summary of the invention comprises three parts.In first portion, traditional Tooth Contact Analysis formula is proposed; In the second portion, derive each sum of errors amount of deformation formula (transmission error, shaft distortion, bearing deformation, tooth distortion, flank of tooth machining error and lathe step-up error); In the third part, each sum of errors amount of deformation is introduced the loading Tooth Contact Analysis formula of deriving and making new advances in the Tooth Contact Analysis formula, proved the method that the present invention proposes by example at last.
By following description and accompanying drawings, the present invention can be more clear, and description of drawings is used for explaining the inventive method and embodiment.
Description of drawings
Fig. 1 summary flow chart of the present invention;
Fig. 2 coordinate system figure that rolls;
Fig. 3 gear tooth thickness error schematic diagram;
Fig. 4 lathe is apart from schematic diagram;
Fig. 5 cutter spacing sets angle τ and helixangle mBetween relation principle figure;
Fig. 6 cutterhead head radius r and helixangle mBetween relation principle figure;
The routing routing diagram of Fig. 7 error load Tooth Contact Analysis;
Fig. 8 is at the following transmission error curve of loading environment.
Embodiment
Shown in Fig. 1-8, a kind of method at overlikon spiral bevel gear error loading Tooth Contact Analysis of the present invention, its concrete steps are as follows:
1 Tooth Contact Analysis
The three dimensional space distance vector of conjugation point of contact
Figure BDA00003071680300061
Δ R d → = - H p l → + J p r → + V e l → - - - ( 1 )
In the formula, V be gear pair parameter vertically is set, and the upward direction of small gear is initiatively.H is that the coordinate axes of gear pair arranges parameter, and is initiatively towards the big extreme direction of small gear.J is the adjustment parameter along change in coordinate axis direction, and small gear is initiatively away from the direction of gearwheel, With
Figure BDA00003071680300064
It is the unit vector in the three dimensional space.
H = Δ R d → · ( e l → × p r → ) J = Δ R d → · ( e → l × p l → ) V = - Δ R d → · e l → - - - ( 2 )
Formula (1) and formula (2) have provided the Tooth Contact Analysis formula.
2 distortion and transmission errors
The distortion of system has comprised the resiliently deformable of gear, the distortion of gear shaft and the resiliently deformable of spring bearing.
2.1 the transmission error under loading environment
Transmission error under the loading environment is made up of three parts: position deviation error delta E a, the transmission error Δ E that the gear resiliently deformable causes bThe transmission error Δ E that causes with the bending deflection of gear cTherefore, comprehensive transmission error Δ E is
ΔE=ΔE a+ΔE b+ΔE c (3)
2.2 load distributes and the compatibility of deformation rule
Suppose to have the tooth of n to mesh, can obtain the load distribution relation
T = Σ i = 1 n T i = T 1 + T 2 + . . . + T n T i = 0.5 d i · F i · cos α i - - - ( 4 )
In the formula, T is resultant couple, T iBe the moment of the gear pair of i, d iBe the gear diameter of point of contact, F iBe the normal force of point of contact, α iLoading direction and tangential angle.
The right transmission deviation of each engaging tooth is identical.Therefore, release the gear transmission formula
Δ E 1 = Δ E 2 = . . . = Δ E n Δ E i = Δ E ia = Δ E ib = Δ E ic - - - ( 5 )
In the formula, Δ E iIt is the transmission error of the gear pair of i.
2.3 the distortion of axle
Fig. 2 has shown at spur roller gear system of coordinates O r(Ol) point of intersection of system axis: O BrBe the center of gear compound graduation circle,
Figure BDA00003071680300068
Be that point arrives point of contact p,
Figure BDA00003071680300069
Point is to the big end of gear, and the axle of gear reverses and is
Figure BDA000030716803000610
The bending deflection of gear is f Xr, f YrAnd f ZrThe reversing with the bending deflection parameter of axle of small gear is similar to gearwheel, and relevant parameter is
Figure BDA000030716803000611
f Xl, f Yl, and f Zl
The supporting structure of overlikon spiral bevel gear can be reduced to an overhang or a support beam.Use F iAnd T iCan calculate bending deflection and torsional deflection.Next step can calculate the displacement of point of contact
f → ai = ( f xr , f yr , f zr ) - - - ( 6 )
2.4 bearing deformation
Consider supporting structure, F iProduce one
Figure BDA00003071680300073
Power.Because the effect of axial force, the axial displacement of small gear is f zBecause effect and the bearing radial internal clearance of tangential force and radial force, be f along the displacement of the axial bearing centre of y yBecause effect and the bearing axial internal clearance of axial force and radial force, be f along the displacement of the axial bearing centre of x xTherefore, the flank of tooth displacement that is produced by bearing deformation is
f → bi = ( f x , f y , f z ) f x = sign ( R zi ) ( δ z + u a ) f y = sign ( R yi ) ( δ y + u r / 2 ) cos ϵ f z = sign ( R xi ) ( δ x + u r / 2 ) sin ϵ - - - ( 7 )
Formula (7), sign (R) are sign functions; u a, u rBe respectively radial internal clearance and the axial internal clearance of bearing; ε is radially opposite force and the axially angle of opposite force; δ x, δ y, δ yBe axial deformation and the radial deformation that is caused by axial force and radial force effect.
2.5 tooth distortion
1) tooth at the contact point of gear surface place is out of shape
According to the Hertz formula, the tooth distortion at the contact point of gear surface place can expand into
δ c = λ 3 9 128 AF 2 ( 1 - μ 1 2 E 1 + 1 - μ 2 2 E 2 ) 2 - - - ( 8 )
In the formula, λ is a coefficient in the elastomechanics; E 1, E 2, μ 1, μ 2Be respectively Young's modulus and Poisson's ratio; A is the relative mean value of curvature; F is contact force.
2) tooth bending deflection
Calculate the bending deflection of tooth with the overhang method of Westinghouse.Bending deflection value δ in contact point of gear surface wFor
δ w = FL 3 3 EI { 1 + 1.3 ( t / L ) + [ 0.25 + 0.75 ( 1 - μ ) ] · ( t / L ) 2 + 0.35 ( t / L ) 3 } - 1 - - - ( 9 )
In the formula, I is the quadrature moment of inertia in load cross section; T is the normal tooth thickness in load cross section; L is the tooth depth of load(ing) point; E is elasticity modulus of materials; μ is the material Poisson's ratio.
3 flank of tooth machining errors
In gear coordinate system (Fig. 3), can obtain the transverse tooth thickness in a method phase cross section by the vector conversion of contact points.Transverse tooth thickness
Figure BDA00003071680300077
For
S → r = R → br - R → br ′ - - - ( 10 )
In the formula,
Figure BDA00003071680300082
With
Figure BDA000030716803000812
Radial vector for gear.Represent that with vector symbol tooth thickness error is
T → gr 0 = T → gr S → r | S → r | - - - ( 11 )
In the formula, It is the tooth thickness error of gear.From the gear tooth thickness error of the rotational coordinates system of Fig. 6 .1, can obtain
T → gr = M ( - δ m ) j · T → gr 0 - - - ( 12 )
In the formula, δ mIt is working pressure angle; Tooth thickness error
Figure BDA00003071680300086
Axially, radially with tangential direction on influence the gap is set, also can be at V, H produces the variable of adjusted value on the J direction.Therefore, can be the error of gap variable as an actual effect
δ jgr ′ = T → gr · u → r δ kgr ′ = T → gr · k → gr δ igr ′ = T → gr · i → gr - - - ( 13 )
In the formula,
Figure BDA00003071680300088
Be respectively axially, radially with tangential direction on unit an amount of.Consider tooth thickness error and tangential resultant error, the actual effect error is
δ jgr ′ = ( T → gr + k 1 · Δ F ir ′ · R r / r 1 ) · T → gr · u → r T → gr δ kgr ′ = ( T → gr + k 1 · ΔF ir ′ · R r / r 1 ) · T → gr · k → gr T → gr δ igr ′ = ( T → gr + k 1 · ΔE ir ′ · R r / r 1 ) · T → gr · i → gr T → gr - - - ( 14 )
In the formula, Δ F ' IrIt is the tangential resultant error of gear; r 1It is the Pitch radius of gear; k 1It is conversion coefficient.Coefficient k 1Be to use apportionment ratio, the relative compensating effect of tooth thickness error and tangential resultant error are determined.k 1Definition be associated with apportionment ratio and the total number of teeth in engagement of actual tangential resultant error (meeting normal distribution law usually).
4 lathe step-up errors
The composition that lathe arranges parameter error is: the error delta τ of cutter spacing setting angle τ; By cutterhead head center to carriage center M dThe lathe that produces of distance apart from error delta M dWith the error delta r in cutterhead head radius r.
4.1 cutter M dAnd concern β between the helix angle m
Fig. 4 can obtain lathe apart from M d, helixangle m, cutterhead head radius r, and the relation between the pitch cone radius R.
M d 2 = R 2 + r 2 - 2 R · r sin ( β m - v ) - - - ( 15 )
sin v = M d Z 0 2 r = Z 0 cos β m Z p r - - - ( 16 )
In the formula, cos β m=M dZ p/ (2R), v=sin -1(0.5M dZ 0/ r) and R be pitch cone radius.Solved M dDifferential, R=R on a M mAnd β=β m, obtain
M dM d′=R mR′+rr′-R′rsin(β m-v)-R mr′sin(β m-v)-R mr(1-v′)cos(β m-v) (17)
v ′ = Z 0 [ r ( R ′ cos β m - R m sin β m ) - R m r ′ cos β m ] Z p r 2 cos v - - - ( 18 )
4.2 cutter spacing sets angle τ and helixangle mBetween relation
τ=2sin -1(M d/ 2k), wherein k is the cutter range setting, can be obtained setting the lathe of angular error generation apart from error delta M because of cutter spacing by this relation d
Δ M d = k cos τ 2 Δτ - - - ( 19 )
Lathe is Δ M apart from error dThe time, Fig. 5 has shown that the helix angle of tooth trace changes.But, R mWith r be constant, therefore
M dM d′=-R mr(1-v′)cos(β m-v) (20)
v ′ = - Z 0 R m sin β m Z p r cos v β ′ - - - ( 21 )
Three formula (19) (20) (21) according to the front can calculate cutter spacing and set angular error
Δτ = τ ′ = - R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) k M d Z p cos v - - - ( 22 )
In the formula, cutterhead head number is Z 0, the crown gear number is Z pAccording to top formula (22), the relation between central point and the helix angle sets angular error Δ τ definition by cutter spacing.
4.3 cutterhead head radius r and helixangle mBetween relation
Be generally used for adjusting the error delta r of the blade pad of eccentric cutter and the combined error definition cutting end radius that slip is counted.This error can change the Δ β in the helix angle m, as shown in Figure 6, while R mAnd M dRemain unchanged.
Got by formula (17) and formula (18)
v ′ = - Z 0 R m ( r sin β m + r ′ cos β m ) Z p r 2 cos v - - - ( 23 )
r′[r-R msin(β m-v)]=R mr(1-v′)cos(β m-v) (24)
Provide simultaneously
Δτ = τ ′ = r R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) r 2 Z p cos v - r Z p R m sin ( β m - v ) cos v - R m Z 0 cos ( β m - v ) cos β m - - - ( 25 )
Following formula (25) has provided cutterhead head radius Δ r and helixangle mBetween relation.
4.4 the relation between lathe step-up error and the area of contact position
According to tooth surface equation and Tooth Contact Analysis equation, this paper arranges the parameter error relevant with them by the lathe that imports and sets up a new meshing condition equation
Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) - - - ( 26 )
In the formula, the three dimensional space distance vector of conjugation point of contact Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) ,
Figure BDA00003071680300103
It is respectively the left and right vector that revolves the conjugation point of contact in the processing engagement on the member flank of tooth; Δ r is the error of cutterhead head radius.
5 flank of tooth contact error analysis (V-H-J Equivalent Calculation formula under the load effect)
The gear teeth face contact analysis is based on the rigidity of gear pair and support system.When considering the distortion of support system, the relative position of mesh tooth face will change in the space.This means the support system distortion and establishing error by the lathe that the lathe value of setting produces is equivalent in magnitude.According to formula (6), (7), (26) can obtain
Figure BDA00003071680300102
In the formula, α is the blade pressure angle; R rAnd R lIt is respectively the radius of gearwheel and pinion point of contact; δ r, δ lIt is respectively the generation angle of gearwheel and small gear.
The example of 6 error load Tooth Contact Analysis
Herein, we analyze a gear pair that processes with spiral AMK852 milling machine.Gear pair has method phase modulus 3.510,35 ° of helix angles, 20 ° of pressure angles, facewidth 20.00mm.Small gear and gearwheel are respectively left-handed and dextrorotation.Small gear has 21 teeth, and gearwheel has 26 teeth.We calculate a series of parameter according to the routing routing diagram of Fig. 7 in test.
Fig. 8 has provided the transmission error curve under loading environment, and line I, II and III have represented the transmission error of engagement beginning, transmission error in the engagement, the transmission error that engagement finishes respectively.The transmission error curve of the transmission error curve of Jia Zaiing and loading is not included on the comprehensive transmission error curve.From simulation result, it is more stable along with the increase of load and drive ratio to find to transmit performance.
Sum up by above instance analysis: this method is reasonably, effectively to the transmission characteristics that improves overlikon spiral bevel gear.The inventive method can apply to crin Gen Beierge bevel gear be used for to solve loading, transmission error, load distribute and machining simulation under real contact area by V-H-J parameter is set and lathe arranges parameter, the influence of compensating error improves the size and location of transmitting performance and revising area of contact.This method can be pressed close to truth more when loading Tooth Contact Analysis to one.

Claims (1)

1. one kind loads the method for Tooth Contact Analysis at the overlikon spiral bevel gear error, it is characterized in that this method may further comprise the steps:
1) distortion and transmission error
1.1) transmission error under loading environment
Transmission error under the loading environment is made up of three parts: position deviation error delta E a, the transmission error Δ E that the gear resiliently deformable causes bThe transmission error Δ E that causes with the bending deflection of gear c, therefore, comprehensive transmission error Δ E is
ΔE=ΔE a+ΔE b+ΔE c (1)
1.2) load distribution and compatibility of deformation rule
Suppose to have the tooth of n to mesh, can obtain the load distribution relation
T = Σ i - 1 n T i = T 1 + T 2 + . . . + T n T i = 0.5 d i · F i · cos α i - - - ( 2 )
In the formula, T is resultant couple, T iBe the moment of the gear pair of i, d iBe the gear diameter of point of contact, F iBe the normal force of point of contact, α iLoading direction and tangential angle;
The right transmission deviation of each engaging tooth is identical, therefore, releases the gear transmission formula
Δ E 1 = Δ E 2 = . . . Δ E n Δ E i = Δ E ia = Δ E ib = Δ E ic - - - ( 3 )
In the formula, Δ E iIt is the transmission error of the gear pair of i;
1.3) axle distortion
The axle of gear reverses and is
Figure FDA00003071680200013
The bending deflection of gear is f Xr, f YrAnd f Zr, the axle relevant parameter of small gear is
Figure FDA00003071680200014
f Xl, f Yl, and f Zl
The supporting structure of overlikon spiral bevel gear can be reduced to an overhang or a support beam, uses Fi and Ti can calculate bending deflection and torsional deflection, and can calculate the displacement of point of contact
Figure FDA00003071680200015
f → ai = ( f xr , f yr , f zr ) - - - ( 4 )
1.4) bearing deformation
Consider supporting structure, Fi produces one
Figure FDA00003071680200017
Power because the effect of axial force, the axial displacement of small gear is f z, because effect and the bearing radial internal clearance of tangential force and radial force, be f along the displacement of the axial bearing centre of y y, because effect and the bearing axial internal clearance of axial force and radial force, be f along the displacement of the axial bearing centre of x x, therefore, the flank of tooth displacement that is produced by bearing deformation is
f → bi = ( f x , f y , f z ) f x = sign ( R zi ) ( δ z + u a ) f y = sign ( R yi ) ( δ y + u r / 2 ) cos ϵ f z = sign ( R xi ) ( δ x + u r / 2 ) sin ϵ - - - ( 5 )
Formula (5), sign (R) are sign functions; u a, u rBe respectively radial internal clearance and the axial internal clearance of bearing; ε is radially opposite force and the axially angle of opposite force; δ x, δ y, δ yBe axial deformation and the radial deformation that is caused by axial force and radial force effect;
1.5) the tooth distortion
1.5.1) be out of shape at the tooth at contact point of gear surface place
According to the Hertz formula, the tooth distortion at the contact point of gear surface place can expand into
δ c = λ 3 9 128 AF 2 ( 1 - μ 1 2 E 1 + 1 - μ 2 2 E 2 ) 2 - - - ( 6 )
In the formula, λ is a coefficient in the elastomechanics; E 1, E 2, μ 1, μ 2Be respectively Young's modulus and Poisson's ratio; A is the relative mean value of curvature; F is contact force;
1.5.2) tooth bending deflection
Calculate the bending deflection of tooth with the overhang method of Westinghouse, at the bending deflection value δ of contact point of gear surface wFor
δ w = FL 3 3 EI { 1 + 1.3 ( t / L ) + [ 0.25 + 0.75 ( 1 - μ ) ] · ( t / L ) 2 + 0.35 ( t / L ) 3 } - 1 - - - ( 7 )
In the formula, I is the quadrature moment of inertia in load cross section; T is the normal tooth thickness in load cross section; L is the tooth depth of load(ing) point; E is elasticity modulus of materials; μ is the material Poisson's ratio;
2) flank of tooth machining error
Consider tooth thickness error and tangential resultant error, the actual effect error is
δ jgr ′ = ( T → gr + k 1 · Δ F ir ′ · R r / r 1 ) · T → gr · u → r T → gr δ kgr ′ = ( T → gr + k 1 · ΔF ir ′ · R r / r 1 ) · T → gr · k → gr T → gr δ igr ′ = ( T → gr + k 1 · ΔE ir ′ · R r / r 1 ) · T → gr · i → gr T → gr - - - ( 8 )
In the formula,
Figure FDA00003071680200025
It is the tooth thickness error of gear;
Figure FDA00003071680200026
Be respectively axially, radially with tangential direction on unit an amount of; R rThe radially radius of tooth; Δ F ' IrIt is the tangential resultant error of gear; r 1It is the Pitch radius of gear; k 1Be conversion coefficient, coefficient k 1Be to use apportionment ratio, the relative compensating effect of tooth thickness error and tangential resultant error are determined, k 1Definition be associated with apportionment ratio and the total number of teeth in engagement of the actual tangential resultant error that meets normal distribution law;
3) lathe step-up error
The composition that lathe arranges parameter is: the error delta τ of cutter spacing setting angle τ; By cutterhead head center to carriage center M dThe lathe that produces of distance apart from error delta M dWith the error delta r in cutterhead head radius r;
3.1) cutter spacing sets angle τ and helixangle mBetween relation
τ=2sin -1(M d/ 2k), wherein k is the cutter range setting, can be obtained setting the lathe of angular error generation apart from error delta M because of cutter spacing by this relation d
Δ M d = k cos τ 2 Δτ - - - ( 9 )
Cutter spacing sets angular error
Δτ = τ ′ = - R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) k M d Z p cos v - - - ( 10 )
In the formula, cos β m=M dZ p/ (2R), v=sin -1(0.5M dZ 0/ r) and R be pitch cone radius, R=R mAnd β=β m, lathe is apart from M d, helixangle m, cutterhead head radius r and pitch cone radius radius R; Cutterhead head number is Z 0, the crown gear number is Z p
3.2) cutterhead head radius r and helixangle mBetween relation
Be generally used for adjusting the error delta r of the blade pad of eccentric cutter and the combined error definition cutting end radius that slip is counted, this error can change the Δ β in the helix angle m, while R mAnd M dRemain unchanged,
Δτ = τ ′ = r R m ( r Z p cos v + R m Z 0 sin β m ) cos ( β m - v ) r 2 Z p cos v - r Z p R m sin ( β m - v ) cos v - R m Z 0 cos ( β m - v ) cos β m - - - ( 11 )
Following formula (11) has provided cutterhead head radius Δ r and helixangle mBetween relation;
3.3) relation between lathe step-up error and the area of contact position
According to tooth surface equation and Tooth Contact Analysis equation, by the importing lathe parameter error relevant with them is set and sets up a new meshing condition equation
Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) - - - ( 12 )
In the formula, the three dimensional space distance vector of conjugation point of contact Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) , Δ R d → ( Δr , Δτ ) = R br p → ( Δr , Δτ ) - R bl p → ( Δr , Δτ ) It is respectively the left and right vector that revolves the conjugation point of contact in the processing engagement on the member flank of tooth; Δ r is the error of cutterhead head radius.
4) flank of tooth contact error analysis
The gear teeth face contact analysis is based on the rigidity of gear pair and support system, when considering the distortion of support system, the relative position of mesh tooth face will change in the space, this means the support system distortion and establishing error by the lathe that the lathe value of setting produces is equivalent in magnitude, according to formula (4), (5), (12) can obtain
Figure FDA00003071680200041
In the formula, α is the blade pressure angle; R rAnd R lIt is respectively the radius of gearwheel and pinion point of contact; δ r, δ lIt is respectively the generation angle of gearwheel and small gear; V be gear pair parameter vertically is set, and the upward direction of small gear is initiatively; H is that the coordinate axes of gear pair arranges parameter, and is initiatively towards the big extreme direction of small gear; J is the adjustment parameter along change in coordinate axis direction, and small gear is initiatively away from the direction of gearwheel.
CN2013101369269A 2013-04-19 2013-04-19 Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors Pending CN103206515A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013101369269A CN103206515A (en) 2013-04-19 2013-04-19 Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013101369269A CN103206515A (en) 2013-04-19 2013-04-19 Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors

Publications (1)

Publication Number Publication Date
CN103206515A true CN103206515A (en) 2013-07-17

Family

ID=48753849

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013101369269A Pending CN103206515A (en) 2013-04-19 2013-04-19 Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors

Country Status (1)

Country Link
CN (1) CN103206515A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105893327A (en) * 2016-03-31 2016-08-24 重庆大学 Method for quickly computing elasticity deformation of deep groove ball bearing and angular contact ball bearing based on FFT (fast fourier transform)
CN106763642A (en) * 2017-02-22 2017-05-31 北京新能源汽车股份有限公司 Electric vehicle speed reducer noise-reduction method and electric vehicle speed reducer
CN107167105A (en) * 2017-06-23 2017-09-15 湖北工业大学 A kind of error compensating method of cycloid gear detection
CN107560583A (en) * 2017-08-25 2018-01-09 重庆市计量质量检测研究院 The axial line bearing calibration of cylindrical workpiece and its diameter measuring method for being segmented section
CN108628254A (en) * 2018-03-30 2018-10-09 重庆大学 Power causes error lower rolling tooth to process tooth surface parameters acquisition methods
CN109482983A (en) * 2018-11-09 2019-03-19 重庆理工大学 A kind of teeth grinding method of generating overlikon spiral bevel gear
CN109670276A (en) * 2019-02-19 2019-04-23 北京工业大学 Consider that the face gear driving error of processing and installation error measures simulating analysis
CN109931367A (en) * 2017-12-18 2019-06-25 浙江宏业农装科技股份有限公司 A kind of manufacturing method of the engaging structure of linear type bevel gear, speed reducer and the gear
CN110702048A (en) * 2019-10-29 2020-01-17 重庆齿轮箱有限责任公司 Contact spot detection method for fan-shaped thrust bearing

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1034616A (en) * 1988-01-29 1989-08-09 国家机械工业委员会成都工具研究所 Gear error engagement separating and measuring method and measuring machine
CN101025384A (en) * 2006-02-18 2007-08-29 山东理工大学 Tooth-face contact-area quantized testing method
CN101294868A (en) * 2008-06-27 2008-10-29 北京工业大学 Single tooth type gear global error measuring apparatus and method
CN101915667A (en) * 2010-07-23 2010-12-15 北京工业大学 Integrated error measuring technology and method thereof of gear pair

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1034616A (en) * 1988-01-29 1989-08-09 国家机械工业委员会成都工具研究所 Gear error engagement separating and measuring method and measuring machine
CN101025384A (en) * 2006-02-18 2007-08-29 山东理工大学 Tooth-face contact-area quantized testing method
CN101294868A (en) * 2008-06-27 2008-10-29 北京工业大学 Single tooth type gear global error measuring apparatus and method
CN101915667A (en) * 2010-07-23 2010-12-15 北京工业大学 Integrated error measuring technology and method thereof of gear pair

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘志峰等: "基于误差的摆线锥齿轮动态加载接触分析方法(上、下)", 《机床与液压》, no. 23, 31 December 2005 (2005-12-31) *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105893327A (en) * 2016-03-31 2016-08-24 重庆大学 Method for quickly computing elasticity deformation of deep groove ball bearing and angular contact ball bearing based on FFT (fast fourier transform)
CN105893327B (en) * 2016-03-31 2018-06-05 重庆大学 Deep groove ball bearing and angular contact ball bearing flexible deformation quick calculation method based on FFT
CN106763642A (en) * 2017-02-22 2017-05-31 北京新能源汽车股份有限公司 Electric vehicle speed reducer noise-reduction method and electric vehicle speed reducer
CN106763642B (en) * 2017-02-22 2020-09-04 北京新能源汽车股份有限公司 Noise reduction method for electric automobile speed reducer and electric automobile speed reducer
CN107167105A (en) * 2017-06-23 2017-09-15 湖北工业大学 A kind of error compensating method of cycloid gear detection
CN107560583B (en) * 2017-08-25 2019-04-26 重庆市计量质量检测研究院 The axial line bearing calibration of cylindrical workpiece and its diameter measuring method for being segmented section
CN107560583A (en) * 2017-08-25 2018-01-09 重庆市计量质量检测研究院 The axial line bearing calibration of cylindrical workpiece and its diameter measuring method for being segmented section
CN109931367A (en) * 2017-12-18 2019-06-25 浙江宏业农装科技股份有限公司 A kind of manufacturing method of the engaging structure of linear type bevel gear, speed reducer and the gear
CN108628254A (en) * 2018-03-30 2018-10-09 重庆大学 Power causes error lower rolling tooth to process tooth surface parameters acquisition methods
CN108628254B (en) * 2018-03-30 2020-03-10 重庆大学 Method for acquiring parameters of gear hobbing tooth surface under force-induced error
CN109482983A (en) * 2018-11-09 2019-03-19 重庆理工大学 A kind of teeth grinding method of generating overlikon spiral bevel gear
CN109482983B (en) * 2018-11-09 2020-01-31 重庆理工大学 method for grinding generating cycloidal-tooth bevel gear
CN109670276A (en) * 2019-02-19 2019-04-23 北京工业大学 Consider that the face gear driving error of processing and installation error measures simulating analysis
CN109670276B (en) * 2019-02-19 2023-05-02 北京工业大学 Face gear transmission error measurement simulation analysis method considering machining and installation errors
CN110702048A (en) * 2019-10-29 2020-01-17 重庆齿轮箱有限责任公司 Contact spot detection method for fan-shaped thrust bearing
CN110702048B (en) * 2019-10-29 2021-03-02 重庆齿轮箱有限责任公司 Contact spot detection method for fan-shaped thrust bearing

Similar Documents

Publication Publication Date Title
CN103206515A (en) Loaded tooth surface contact analysis method direct at epicycloid bevel gear errors
Deng et al. Parametric study of meshing characteristics with respect to different meshing rollers of the antibacklash double-roller enveloping worm gear
CN103105152B (en) A kind of Gear Processing Analysis of Surface Topography method based on simulation of gear machining model
CN101937482B (en) Design method for face gear slotting cutter tooth form
Deng et al. Modification design method for an enveloping hourglass worm gear with consideration of machining and misalignment errors
CN105522227A (en) Profile shaping method for machining cutter tooth main cutting edge according to cycloid bevel gear
Shih et al. A novel method for producing a conical skiving tool with error-free flank faces for internal gear manufacture
Simon Machine-tool settings to reduce the sensitivity of spiral bevel gears to tooth errors and misalignments
Kujawski et al. Analysis of generation capabilities of noncircular cogbelt pulleys on the example of a gear with an elliptical pitch line
CN105223814A (en) Involute gear formed grinding wheel computing method
Tsuji et al. Tooth contact analysis and manufacture on multitasking machine of large-sized straight bevel gears with equi-depth teeth
Zheng et al. Investigation on the wear of spur gears generated by modified cutter
Zeng et al. Modelling of straight bevel gear transmission and simulation of its meshing performance
Shih Mathematical model for face-hobbed straight bevel gears
Liang et al. Mathematical model and meshing analysis of internal meshing gear transmission with curve element
CN103268368A (en) Klingelnberg bevel gear contact regulating method
Uzun et al. Manufacturing the new type concave–convex profile involute gears modeled by CAD–CAM in CNC milling machines
CN107433454A (en) The ram flexibility compensating device and boring and milling machine of boring and milling machine
CN103207939A (en) Cycloidal tooth bevel gear dynamical modeling method based on variable transmission ratio
Tang et al. DYNAMIC CONTACT ANALYSIS AND TOOTH MODIFICATION DESIGN FOR EMU TRACTION GEAR.
Radzevich About hob idle distance in gear hobbing operation
Patel et al. Topology Optimization of Automotive Gear using Fea
Kawasaki et al. Accuracy measurement and evaluation of straight bevel gear manufactured by end mill using CNC milling machine
CN103793564A (en) System deformation calculation method of transmission gearbox
Khodaee et al. Finite element simulation as a tool to evaluate gear quality after gear rolling

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130717