CN113446960A - Tooth surface point cloud theoretical distribution modeling method and measuring method - Google Patents

Tooth surface point cloud theoretical distribution modeling method and measuring method Download PDF

Info

Publication number
CN113446960A
CN113446960A CN202110719148.0A CN202110719148A CN113446960A CN 113446960 A CN113446960 A CN 113446960A CN 202110719148 A CN202110719148 A CN 202110719148A CN 113446960 A CN113446960 A CN 113446960A
Authority
CN
China
Prior art keywords
point cloud
tooth surface
tooth
gear
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110719148.0A
Other languages
Chinese (zh)
Other versions
CN113446960B (en
Inventor
张祥雷
洪天野
周宏明
吴雨韩
陈体兴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wenzhou University
Original Assignee
Wenzhou University
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 Wenzhou University filed Critical Wenzhou University
Priority to CN202110719148.0A priority Critical patent/CN113446960B/en
Publication of CN113446960A publication Critical patent/CN113446960A/en
Application granted granted Critical
Publication of CN113446960B publication Critical patent/CN113446960B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2416Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures of gears

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention discloses a tooth surface point cloud theoretical distribution modeling method and a measuring method, wherein an established model is a method for measuring gear precision based on line structure light. The modeling method comprises the steps of firstly establishing a gear coordinate system and a sensor coordinate system respectively according to the position relation of a gear and a linear structured light sensor during actual measurement, deducing a light path equation of linear structured light by using a coordinate transformation method, and then solving the light path equation in series with a standard tooth surface equation to obtain an accurate tooth surface point cloud theoretical distribution position. The modeling method can calculate the theoretical distribution point cloud of the tooth surface before actual measurement, is convenient for optimizing the measurement position and angle of the line structured light sensor in advance, and can also be used for comparing the measured actual point cloud with the theoretical point cloud point by point to realize rapid denoising and rapid calculation of gear approximation error.

Description

Tooth surface point cloud theoretical distribution modeling method and measuring method
Technical Field
The invention belongs to the technical field of gear detection, and particularly relates to a tooth surface point cloud theoretical distribution modeling method and a measuring method.
Background
The current trend in gear measurement is fast, fully digital. In order to quickly acquire three-dimensional shape information of a gear tooth surface, more and more students begin to acquire a gear tooth surface point cloud by using various optical measurement methods to calculate gear errors, and linear structured light measurement is one of the methods. However, the following problems often exist in the process of actually measuring and calculating gear errors:
1. the measurement method is empirical. The measurement effect of the tooth surface point cloud cannot be accurately predicted before measurement, and the method (the position where a measuring head is placed and the measurement angle) used in the measurement of the linear structured light mainly depends on experience and test;
2. the method of calculating the error using the actual measured point cloud is complex. The current error calculation is mainly performed by calculating the difference between each point and the theoretical position of each point. However, because the analytic formulas of all parts of the gear tooth profile are different, the part to which the point belongs needs to be judged before point-by-point cloud calculation, so that the operations of judging an involute starting point, fitting the involute tooth profile and the like need to be carried out first, and the requirement on operators is high;
3. the method for preprocessing the inclinometer gear point cloud is complex. When the straight-tooth cylindrical gear is measured by linear structured light parallel to the axis of the gear, the obtained data has strong relevance, and because the intersection points of the light path and the tooth surface are all positioned on a straight line vertical to the axis, the measured point cloud is easy to fit the intersected straight line, and noise points are removed. However, when the helical gear is measured, the intersection point is located on the curve, the intersection curves of different corners are different, and the data correlation is not strong, so that the algorithm required when the noise point of the helical gear is removed is complex.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a tooth surface point cloud theoretical distribution modeling method and a measuring method
The technical scheme adopted by the invention is as follows: a tooth surface point cloud theoretical distribution modeling method for measuring gear precision based on line structured light comprises the following steps:
w1, modeling an optical path equation of linear structured light:
establishing a gear coordinate system Og-XgYgOptical path coordinate system OL-XLYL(ii) a Wherein, OLThe point is the initial position of the measuring head, and the light ray is from OLStarting to XLThe negative axis direction; the distance between the measuring head and the tooth center is set as D; xgThe included angle between the shaft and the connecting line of the measuring head and the gear center is sigma1And XLAngle of axis being sigma2(ii) a Then O isLPoint-on-gear coordinate system Og-xgygThe coordinates in (d) can be expressed as (Dcos σ @1,Dsinσ1) The initial slope of the optical path is k ═ tan σ2
The equation for the initial position of the optical path can be expressed as:
y-D sinσ1=-tanσ2(x-Dcosσ1) (1);
O′L-X′LY′Lis O after a rotation time tL-XLYLThe position is defined as the position with the rotation speed of omega and the rotation angle of W-omega; then O'LPoint-on-gear coordinate system Og-xgygThe coordinate in (D cos (W + σ) is1),D sin(W+σ1) The slope of the optical path is k' ═ -tan (σ)2-W);
The equation of the optical path after the rotation time t can be expressed as:
y-D sin(W+σ1)=-tan(σ2-W)[x-Dcos(W+σ1)] (2);
the equation for the optical path plane can be expressed as:
Figure BDA0003136281280000021
wherein B is the tooth width, tbIs the tooth width factor;
w2. three-dimensional mathematical model of standard involute gear:
a tooth surface equation can be deduced according to a generating mode of a spiral surface and an end surface tooth profile equation, theta represents the angle of the rotation of a bus around a rotating shaft, the right-hand rotation is positive, the left-hand rotation is negative, delta represents the initial included angle between the Yw axis of a fixed coordinate system Ow and the tooth thickness symmetrical line of the measured gear,
the equation for the involute helicoid is:
Figure BDA0003136281280000031
the flank equation for the root transition curve is:
Figure BDA0003136281280000032
the tooth surface equation of the addendum circle is:
Figure BDA0003136281280000033
the tooth surface equation of the root circle is:
Figure BDA0003136281280000034
ωAcan be taken from the value of alpha' in formula (5)
Figure BDA0003136281280000035
Calculating to obtain;
in the above four equations, B is the tooth width, tbV is more than or equal to 0 and less than or equal to 1.0; theta denotes an angle through which the bus bar rotates around the rotation axis,
Figure BDA0003136281280000036
theta is positive when the gear rotates rightwards, and theta is negative when the gear rotates leftwards; delta represents an angle for rotating the Yw axis of the fixed coordinate system Ow to the tooth thickness symmetrical line of the measured gear, and clockwise is positive and anticlockwise is negative;
when any tooth surface is worked out, the delta values in the above four groups of equations can be used
Figure BDA0003136281280000037
Substituting to obtain; the meanings and the value ranges of other variables in the formula are the same as those in the corresponding end face tooth profile equation, and only the end face parameters are required to replace the normal face parameters;in the formula, the method for taking +/-and +/-m' is also the same as the equation of the tooth profile of the end face;
w3. calculating the tooth surface distribution point cloud:
the intersection line of the light path and the tooth surface is a complex curve, so that the z-axis coordinate of each sampling point is firstly determined when the point cloud of the intersection point is calculated, and then the formula (3) and the formulas (4), (5), (6) and (7) are sequentially connected and solved according to the intersection sequence of the light path and each part of the tooth surface to obtain the point cloud of the tooth surface;
the number of equally spaced sampling points on the intersecting curve of the light path and the tooth surface is ntbThen the sampling interval is
Figure BDA0003136281280000041
The z-axis coordinate of the sampling point in each sampling is
Figure BDA0003136281280000042
Let the sampling frequency be f, and the tooth surface point cloud number measured after the rotation time t be Nt=ωtfntb
Because each tooth shape of the standard gear is the same, only the rotation angle needs to be calculated in the process of calculating the point cloud of the tooth surface
Figure BDA0003136281280000043
Namely, it is
Figure BDA0003136281280000044
And (3) carrying out point counting, rotating to obtain a tooth surface point cloud of another z-1 tooth, thus obtaining a complete tooth surface point cloud, wherein N is 2 pi fntbAnd (5) rotating to obtain the tooth surface point cloud of the other z-1 tooth, thus obtaining the complete tooth surface point cloud.
A gear precision detection method adopting the tooth surface point cloud theoretical distribution modeling method for measuring the gear precision based on the line structured light comprises the following steps:
the method comprises the following steps: establishing a tooth surface point cloud theoretical distribution model by adopting the tooth surface point cloud theoretical distribution modeling method based on the line structured light measurement gear precision;
step two: the tooth surface is measured through line structured light, and the specific process is as follows: the plane where the light is located is parallel to the axis of the gear, and the light path covers all the tooth sides; in the measuring process, the measuring head is fixed, the gear rotates for a circle at a constant speed, and the measuring head collects points on an intersection line segment of a light path and a tooth surface according to the frequency f to obtain a tooth surface point cloud;
step two: removing obvious noise points, which comprises the following specific steps: calculating the Euclidean distance between each actually measured point cloud and a theoretical point cloud in a tooth surface point cloud theoretical distribution model, wherein the difference value exceeding a threshold value is a noise point;
step three: the parameters which are completely the same when the method is used for actual measurement are input into a tooth surface point cloud theoretical distribution model for calculation, and the point cloud which is completely the same as the distribution of the actually measured point cloud can be obtained, so that the difference value between the actually measured point cloud and the theoretically distributed point cloud can be completely obtained point by point after the point cloud is registered, and the difference value can be approximately used as tooth profile deviation because each point and the corresponding point are completely positioned at the same z-axis height.
The invention has the following beneficial effects:
1. the tooth surface theoretical distribution point cloud of the specified measurement parameters can be directly calculated before actual measurement, and the measurement scheme (position and angle) can be conveniently optimized. For example, aiming at the shielding phenomenon of the peripheral teeth to the measuring teeth, parameters can be adjusted in the model, measuring parameters capable of measuring more tooth surfaces are selected, and then actual measurement is carried out.
2. The method is simple and convenient to remove obvious noise points in the process of calculating the precision of the linear structure light detection gear.
3. The error calculation after the actual point cloud measurement is simplified.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is within the scope of the present invention for those skilled in the art to obtain other drawings based on the drawings without inventive exercise.
FIG. 1 is a schematic diagram of a three-dimensional point cloud measurement system;
FIG. 2 is a schematic view of a structured light detection gear;
FIG. 3 is a cross-section of the measuring device at the lower end face of the gear;
FIG. 4 is a three-dimensional mathematical model of the tooth surface of the gear under test;
FIG. 5 is a cross point of a light path and a tooth surface when t is more than or equal to 0.6s and less than or equal to 0.9 s;
FIG. 6 is a complete tooth surface point cloud;
fig. 7 is a flow chart of a method.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings.
The invention provides a tooth surface point cloud theoretical distribution modeling method for measuring gear precision based on line structured light, which comprises the following steps:
w1, modeling an optical path equation of linear structured light:
the method for linear structured light measurement of the tooth surface comprises the following steps: the plane of the light is parallel to the gear axis, and the light path covers all the tooth sides. In the actual measurement process, the measuring head is fixed, the gear rotates for a circle at a constant speed, and the measuring head collects points on the intersection line segment of the light path and the tooth surface according to a certain frequency f to obtain the tooth surface point cloud. The process is subjected to mathematical modeling, in order to facilitate calculation, the measuring method is equivalent to the fact that the gear is still in the calculating process, and the measuring head rotates around the axis of the gear at a constant speed.
As shown in fig. 1, a conventional optical gear measuring table is taken as an example. The measuring table comprises a displacement platform, a rotary platform, a displacement control PC, a line structure light sensor and a measuring head parameter control PC. The linear structure light sensor is connected with the displacement platform, and the PC can be used for controlling X, Y, Z to move in three axial directions. The sensor is horizontally arranged, so that the plane of the line-structured light is parallel to the axis of the gear. The gear is horizontally placed on the rotary platform, and the C shaft can be controlled to rotate by the PC. The linear structure measures the positional relationship of the tooth surfaces as shown in fig. 2, the plane of the light is parallel to the gear axis, and the light path covers all the tooth sides.
FIG. 3 is a cross section of the measuring device at the lower end face of the gear, establishing a gear coordinate system Og-XgYgOptical path coordinate system OL-XLYL(ii) a Wherein, OLThe point is the initial position of the measuring head, and the light ray is from OLStarting to XLThe negative axis direction; the distance between the measuring head and the tooth center is set as D; xgThe included angle between the shaft and the connecting line of the measuring head and the gear center is sigma1And XLAngle of axis being sigma2(ii) a Then O isLPoint-on-gear coordinate system Og-xgygThe coordinates in (d) can be expressed as (Dcos σ @1,Dsinσ1) The initial slope of the optical path is k ═ tan σ2
The equation for the initial position of the optical path can be expressed as:
y-Dsinσ1=-tanσ2(x-Dcosσ1) (1);
O′L-X′LY′Lis O after a rotation time tL-XLYLThe position is defined as the position with the rotation speed of omega and the rotation angle of W-omega; then O'LPoint-on-gear coordinate system Og-xgygThe coordinate in (D cos (W + σ) is1),D sin(W+σ1) The slope of the optical path is k' ═ -tan (σ)2-W);
The equation of the optical path after the rotation time t can be expressed as:
y-D sin(W+σ1)=-tan(σ2-W)[x-D cos(W+σ1)] (2);
the equation for the optical path plane can be expressed as:
Figure BDA0003136281280000061
wherein B is the tooth width, tbIs the tooth width factor;
w2. three-dimensional mathematical model of standard involute gear:
the tooth surface equation of the involute gear can be derived from the tooth surface profile equation. The end tooth profile consists of the following parts: the AB section is a transition curve formed by cutting a tool fillet into a base circle, and the BC section is an involute generated by conjugate of generating motion of the tool and a workpiece; the other parts are tooth top circular arcs and tooth root circular arcs.
The tooth profile equation of the straight gear is deduced, and the tooth profile equation is directly written without being described in detail.
The coordinate of any point k on the tooth profile involute is as follows:
Figure BDA0003136281280000071
pressure angle alpha with independent variable of k pointkThe value range is (alpha)B,αC) The values are:
Figure BDA0003136281280000072
the method for taking the plus or minus sign in the formula is to take an upper end sign when the right side tooth profile of the y axis is calculated and calculate a lower end sign of the left side tooth profile.
The transition curve tooth profile equation is as follows:
Figure BDA0003136281280000073
the value range of alpha' is
Figure BDA0003136281280000074
When alpha' is taken, the transition curve is connected with the involute and is tangent to the point B; get
Figure BDA0003136281280000075
The transition curve is tangent to the dedendum circle at the point A. The method for taking the plus or minus sign in the formula is to take an upper end sign when the right side tooth profile of the y axis is calculated and calculate a lower end sign of the left side tooth profile.
The tooth form of the end surface of the helical gear is similar to that of a straight gear, but each parameter is no longer a standard value. The section perpendicular to the tangent line of the spiral line on the indexing cylinder of the bevel gear is called as a normal surface, and because the gear is cut along the spiral line direction, the tooth profile parameters of the normal surface are consistent with those of the cutter and are standard values. However, since the end face parameters are generally used for geometric calculation of the helical gear model, the end face parameters and the normal face parameters are converted as shown in table 1.
TABLE 1 calculation formula of bevel gear end face parameters
Bevel gear end face parameters (symbol) Formula (II)
Modulus of elasticity mt mt=mn/cosβ
Reference circle pressure angle αt αt=arctan(tanαn/cosβ)
Radius of reference circle rt rt=mtz=(mn/cosβ)z
Radius of base circle rbt rbt=rtcosαt
Coefficient of tooth crest height hat hat=hancosβ
Coefficient of tooth tip clearance ct ct=cncosβ
Radius of addendum circle rat rat=rt+hatmt
Root circle radius rft rat=rt-(hat+ct)mt
In the table, subscripts "n", "t" denote a normal plane and an end plane, respectively, and β is a pitch angle on the pitch circle.
Calculating the cutter parameters of the end face of the bevel gear as follows:
Figure BDA0003136281280000081
a tooth surface equation can be deduced by a generating mode of a spiral surface and an end surface tooth profile equation, as shown in figure 4, theta represents the angle of a bus rotating around a rotating shaft, right-hand rotation is positive, left-hand rotation is negative, delta represents an initial included angle between an Yw axis of a fixed coordinate system Ow and a tooth thickness symmetrical line of a measured gear,
the equation for the involute helicoid is:
Figure BDA0003136281280000082
the flank equation for the root transition curve is:
Figure BDA0003136281280000091
the tooth surface equation of the addendum circle is:
Figure BDA0003136281280000092
the tooth surface equation of the root circle is:
Figure BDA0003136281280000093
ωAcan be taken from the value of alpha' in formula (5)
Figure BDA0003136281280000094
Calculating to obtain;
in the above four equations, B is the tooth width, tbThe increment of the tooth width (v is more than or equal to 0 and less than or equal to 1.0); theta denotes an angle through which the bus bar rotates around the rotation axis,
Figure BDA0003136281280000095
theta is positive when the gear rotates rightwards, and theta is negative when the gear rotates leftwards; delta represents an angle for rotating the Yw axis of the fixed coordinate system Ow to the tooth thickness symmetrical line of the measured gear, and clockwise is positive and anticlockwise is negative;
when any tooth surface is worked out, the delta values in the above four groups of equations can be used
Figure BDA0003136281280000096
Substituting to obtain; the meanings and the value ranges of other variables in the formula are the same as those in the corresponding end face tooth profile equation, and only the end face parameters are required to replace the normal face parameters; in the formula, the method for taking +/-and +/-m' is also the same as the equation of the tooth profile of the end face;
w3. calculating the tooth surface distribution point cloud:
the intersection line of the light path and the tooth surface is a complex curve, so that the z-axis coordinate of each sampling point is firstly determined when the point cloud of the intersection point is calculated, and then the formula (3) and the formulas (4), (5), (6) and (7) are sequentially connected and solved according to the intersection sequence of the light path and each part of the tooth surface to obtain the point cloud of the tooth surface;
the number of equally spaced sampling points on the intersecting curve of the light path and the tooth surface is ntbThen the sampling interval is
Figure BDA0003136281280000101
The z-axis coordinate of the sampling point in each sampling is
Figure BDA0003136281280000102
Let the sampling frequency be f, and the tooth surface point cloud number measured after the rotation time t be Nt=ωtfntb
Because each tooth shape of the standard gear is the same, only the rotation angle needs to be calculated in the process of calculating the point cloud of the tooth surface
Figure BDA0003136281280000103
Namely, it is
Figure BDA0003136281280000104
And (3) carrying out point counting, rotating to obtain a tooth surface point cloud of another z-1 tooth, thus obtaining a complete tooth surface point cloud, wherein N is 2 pi fntbAnd (5) rotating to obtain the tooth surface point cloud of the other z-1 tooth, thus obtaining the complete tooth surface point cloud.
Example 1:
TABLE 2 helical gear parameters
Figure BDA0003136281280000105
TABLE 3 gauge head parameters
Parameter(s) Numerical value
Included angle sigma between initial position of measuring head and tooth center1/(°) 10
Initial angle of light path and xgAngle of axis sigma2/(°) 20
Distance d/mm between measuring head and tooth center 50
Gear speed omega (anticlockwise)/(rad/s) π/2
Sampling frequency f/Hz 2000
Number of sampling points n at a time 10
According to the input parameters, the complete tooth surface measurement needs 4 seconds in the actual measurement, and the minimum measurement time is
Figure BDA0003136281280000111
And second. Fig. 5 is a calculation result when t is set to be 0.6s or more and 0.9s or less, a red line is a light ray of a sampling point, 10 points are taken at equal intervals for each frequency on an intersecting curved surface of a light path and a tooth surface, that is, intersection points of 10 light rays and the tooth surface are calculated, blue points are intersection points of the light rays and the tooth surface, the total number of the intersection points is 6000, 4000 (that is, the number of intersection points of one tooth) of any continuous intersection points are taken, a tooth surface point cloud of the other z-1 tooth is calculated according to an included angle between two teeth, and 80000 complete tooth surface point clouds are finally obtained, as shown in fig. 6.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by relevant hardware instructed by a program, and the program may be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc.
The above disclosure is only for the purpose of illustrating the preferred embodiments of the present invention, and it is therefore to be understood that the invention is not limited by the scope of the appended claims.

Claims (2)

1. A tooth surface point cloud theoretical distribution modeling method for measuring gear precision based on line structured light is characterized by comprising the following steps:
w1, modeling an optical path equation of linear structured light:
establishing a gear coordinate system Og-XgYgOptical path coordinate system OL-XLYL(ii) a Wherein, OLThe point is the initial position of the measuring head, and the light ray is from OLStarting to XLThe negative axis direction; the distance between the measuring head and the tooth center is set as D; xgThe included angle between the shaft and the connecting line of the measuring head and the gear center is sigma1And XLAngle of axis being sigma2(ii) a Then O isLPoint-on-gear coordinate system Og-xgygThe coordinates in (d) can be expressed as (Dcos σ @1,Dsinσ1) The initial slope of the optical path is k ═ tan σ2
The equation for the initial position of the optical path can be expressed as:
y-Dsinσ1=-tanσ2(x-Dcosσ1) (1);
OL′-XL′YL' is O after a rotation time tL-XLYLThe position is defined as the position with the rotation speed of omega and the rotation angle of W-omega; then O isLPoint-on-gear coordinate system Og-xgygThe coordinate in (D) is (Dcos (W + σ)1),Dsin(W+σ1) The slope of the optical path is k' ═ -tan (σ)2-W);
The equation of the optical path after the rotation time t can be expressed as:
y-Dsin(W+σ1)=-tan(σ2-W)[x-Dcos(W+σ1)] (2);
the equation for the optical path plane can be expressed as:
Figure FDA0003136281270000011
wherein B is the tooth width, tbIs the tooth width factor;
w2. three-dimensional mathematical model of standard involute gear:
a tooth surface equation can be deduced according to a generating mode of a spiral surface and an end surface tooth profile equation, theta represents the angle of the rotation of a bus around a rotating shaft, the right-hand rotation is positive, the left-hand rotation is negative, delta represents the initial included angle between the Yw axis of a fixed coordinate system Ow and the tooth thickness symmetrical line of the measured gear,
the equation for the involute helicoid is:
Figure FDA0003136281270000021
the flank equation for the root transition curve is:
Figure FDA0003136281270000022
the tooth surface equation of the addendum circle is:
Figure FDA0003136281270000023
the tooth surface equation of the root circle is:
Figure FDA0003136281270000024
ωAcan be taken from the value of alpha' in formula (5)
Figure FDA0003136281270000025
Calculating to obtain;
in the above four equations, B is the tooth width, tbV is more than or equal to 0 and less than or equal to 1.0; theta denotes an angle through which the bus bar rotates around the rotation axis,
Figure FDA0003136281270000026
theta is positive when the gear rotates rightwards, and theta is negative when the gear rotates leftwards; delta represents an angle for rotating the Yw axis of the fixed coordinate system Ow to the tooth thickness symmetrical line of the measured gear, and clockwise is positive and anticlockwise is negative;
when any tooth surface is worked out, the delta values in the above four groups of equations can be used
Figure FDA0003136281270000027
Substituting to obtain; the meanings and the value ranges of other variables in the formula are the same as those in the corresponding end face tooth profile equation, and only the end face parameters are required to replace the normal face parameters; in the formula, "+/-)"
Figure FDA0003136281270000035
The method is also the same as the equation of the tooth profile of the end face;
w3. calculating the tooth surface distribution point cloud:
the intersection line of the light path and the tooth surface is a complex curve, so that the z-axis coordinate of each sampling point is firstly determined when the point cloud of the intersection point is calculated, and then the formula (3) and the formulas (4), (5), (6) and (7) are sequentially connected and solved according to the intersection sequence of the light path and each part of the tooth surface to obtain the point cloud of the tooth surface;
the number of equally spaced sampling points on the intersecting curve of the light path and the tooth surface is ntbThen the sampling interval is
Figure FDA0003136281270000031
The z-axis coordinate of the sampling point in each sampling is
Figure FDA0003136281270000032
Let the sampling frequency be f, and the rotation time be measuredThe number of tooth surface point clouds is Nt=ωtfntb
Because each tooth shape of the standard gear is the same, only the rotation angle needs to be calculated in the process of calculating the point cloud of the tooth surface
Figure FDA0003136281270000033
Namely, it is
Figure FDA0003136281270000034
And (3) carrying out point counting, rotating to obtain a tooth surface point cloud of another z-1 tooth, thus obtaining a complete tooth surface point cloud, wherein N is 2 pi fntbAnd (5) rotating to obtain the tooth surface point cloud of the other z-1 tooth, thus obtaining the complete tooth surface point cloud.
2. A gear accuracy detection method using the tooth surface point cloud theoretical distribution modeling method for measuring gear accuracy based on line structured light according to claim 1, characterized by comprising the steps of:
the method comprises the following steps: the tooth surface point cloud theoretical distribution modeling method for measuring the gear precision based on the line structured light is adopted to establish a tooth surface point cloud theoretical distribution model according to claim 1;
step two: the tooth surface is measured through line structured light, and the specific process is as follows: the plane where the light is located is parallel to the axis of the gear, and the light path covers all the tooth sides; in the measuring process, the measuring head is fixed, the gear rotates for a circle at a constant speed, and the measuring head collects points on an intersection line segment of a light path and a tooth surface according to the frequency f to obtain a tooth surface point cloud;
step three: removing obvious noise points, which comprises the following specific steps: calculating the Euclidean distance between each actually measured point cloud and a theoretical point cloud in a tooth surface point cloud theoretical distribution model, wherein the difference value exceeding a threshold value is a noise point;
step four: the parameters which are completely the same when the method is used for actual measurement are input into a tooth surface point cloud theoretical distribution model for calculation, and the point cloud which is completely the same as the distribution of the actually measured point cloud can be obtained, so that the difference value between the actually measured point cloud and the theoretically distributed point cloud can be completely obtained point by point after the point cloud is registered, and the difference value can be approximately used as tooth profile deviation because each point and the corresponding point are completely positioned at the same z-axis height.
CN202110719148.0A 2021-06-28 2021-06-28 Tooth surface point cloud theoretical distribution modeling method and measuring method Active CN113446960B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110719148.0A CN113446960B (en) 2021-06-28 2021-06-28 Tooth surface point cloud theoretical distribution modeling method and measuring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110719148.0A CN113446960B (en) 2021-06-28 2021-06-28 Tooth surface point cloud theoretical distribution modeling method and measuring method

Publications (2)

Publication Number Publication Date
CN113446960A true CN113446960A (en) 2021-09-28
CN113446960B CN113446960B (en) 2023-01-03

Family

ID=77813257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110719148.0A Active CN113446960B (en) 2021-06-28 2021-06-28 Tooth surface point cloud theoretical distribution modeling method and measuring method

Country Status (1)

Country Link
CN (1) CN113446960B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117029705A (en) * 2023-06-27 2023-11-10 苏州瑞威盛科技有限公司 Gear bar span measuring system and method based on non-contact 3D vision

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289534A (en) * 2011-05-26 2011-12-21 无锡永凯达齿轮有限公司 Method for modeling involute helical gear accurately
CN102567596A (en) * 2012-03-05 2012-07-11 西北工业大学 Parametric modeling method of face gear teeth
WO2012141235A1 (en) * 2011-04-13 2012-10-18 株式会社トプコン Three-dimensional point group position data processing device, three-dimensional point group position data processing system, three-dimensional point group position data processing method and program
JP2013221845A (en) * 2012-04-16 2013-10-28 Mitsubishi Heavy Ind Ltd Device and method for gear inspection
DE102015000975A1 (en) * 2015-01-23 2016-07-28 Liebherr-Verzahntechnik Gmbh Method for producing a workpiece with corrected tooth geometry and / or modified surface structure
CN107514979A (en) * 2017-07-20 2017-12-26 重庆大学 A kind of gear Integrated Measurement System and its method based on structure light
CN107588742A (en) * 2017-10-25 2018-01-16 北京工业大学 A kind of cylindrical gear profile bias measurement method based on line-structured light
CN107860313A (en) * 2017-10-25 2018-03-30 北京工业大学 It is a kind of based on the beveled gear teeth of line-structured light to bias measurement method
CN108050946A (en) * 2017-11-30 2018-05-18 北京工业大学 A kind of gear tooth thickness measuring method based on line-structured light
CN109978993A (en) * 2019-03-21 2019-07-05 中南大学 Involute helical face gear three-dimensional modeling method
CN110940283A (en) * 2019-12-11 2020-03-31 北京工业大学 High-precision measurement method for tooth pitch deviation and tooth profile deviation of tiny gear based on white light interferometer
CN112580160A (en) * 2020-12-03 2021-03-30 南京工业大学 Calibration method for on-machine measurement system of forming gear grinding machine
CN112935420A (en) * 2021-03-23 2021-06-11 重庆大学 Involute gear shaving cutter and three-dimensional modeling method and gear shaving processing method thereof

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012141235A1 (en) * 2011-04-13 2012-10-18 株式会社トプコン Three-dimensional point group position data processing device, three-dimensional point group position data processing system, three-dimensional point group position data processing method and program
CN102289534A (en) * 2011-05-26 2011-12-21 无锡永凯达齿轮有限公司 Method for modeling involute helical gear accurately
CN102567596A (en) * 2012-03-05 2012-07-11 西北工业大学 Parametric modeling method of face gear teeth
JP2013221845A (en) * 2012-04-16 2013-10-28 Mitsubishi Heavy Ind Ltd Device and method for gear inspection
DE102015000975A1 (en) * 2015-01-23 2016-07-28 Liebherr-Verzahntechnik Gmbh Method for producing a workpiece with corrected tooth geometry and / or modified surface structure
CN107514979A (en) * 2017-07-20 2017-12-26 重庆大学 A kind of gear Integrated Measurement System and its method based on structure light
CN107588742A (en) * 2017-10-25 2018-01-16 北京工业大学 A kind of cylindrical gear profile bias measurement method based on line-structured light
CN107860313A (en) * 2017-10-25 2018-03-30 北京工业大学 It is a kind of based on the beveled gear teeth of line-structured light to bias measurement method
CN108050946A (en) * 2017-11-30 2018-05-18 北京工业大学 A kind of gear tooth thickness measuring method based on line-structured light
CN109978993A (en) * 2019-03-21 2019-07-05 中南大学 Involute helical face gear three-dimensional modeling method
CN110940283A (en) * 2019-12-11 2020-03-31 北京工业大学 High-precision measurement method for tooth pitch deviation and tooth profile deviation of tiny gear based on white light interferometer
CN112580160A (en) * 2020-12-03 2021-03-30 南京工业大学 Calibration method for on-machine measurement system of forming gear grinding machine
CN112935420A (en) * 2021-03-23 2021-06-11 重庆大学 Involute gear shaving cutter and three-dimensional modeling method and gear shaving processing method thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王尧等: "正交直齿面齿轮参数化建模方法", 《机械设计与制造》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117029705A (en) * 2023-06-27 2023-11-10 苏州瑞威盛科技有限公司 Gear bar span measuring system and method based on non-contact 3D vision
CN117029705B (en) * 2023-06-27 2024-03-22 苏州瑞威盛科技有限公司 Gear bar span measuring system and method based on non-contact 3D vision

Also Published As

Publication number Publication date
CN113446960B (en) 2023-01-03

Similar Documents

Publication Publication Date Title
CN109063326B (en) Gear accurate modeling method considering microscopic shape correction and actual machining errors
KR101721969B1 (en) Method for the location determination of the involutes in gears
CN110297462B (en) Gear grinding precision prediction modeling method considering influence of machine tool geometric error
US5136522A (en) Stock dividing method and apparatus for gear manufacturing machine
CN112541235B (en) Universal design method for hypoid gear pair
CN112257192B (en) One-tooth tangential integrated deviation calculation method based on actually measured tooth profile deviation
CN109993464B (en) Machining parameter optimization method for reducing mounting error sensitivity of hypoid gear
CN113446960B (en) Tooth surface point cloud theoretical distribution modeling method and measuring method
CN111666643B (en) Method for determining contact performance of complex tooth surface
CN106493439A (en) A kind of spiral bevel gear tooth top chamfering method that is analyzed based on material contact
US20170008106A1 (en) Method of producing a toothed workpiece having a modified surface geometry
CN106980724A (en) A kind of flank of tooth accurate model section setting-out construction method based on complex teeth surfaces grid planning principles
JP2019512684A (en) Worm gear measurement
Özel A study on cutting errors in the tooth profiles of the spur gears manufactured in CNC milling machine
CN114918494B (en) Design method of instantaneous equivalent model for gear generating and grinding
CN109033669B (en) Spiral bevel gear simulation machining modeling method based on universal motion parameter driving
CN110021066A (en) A kind of prolate involute worm gear hobbing three-dimensional modeling method, device and equipment
JPWO2006022336A1 (en) Evaluation method of edge profile of re-sharpened pinion cutter
CN113124800B (en) Archimedes spiral surface worm wheel rolling shaving processing precision detection method
CN109887079A (en) Spiral bevel gear three-dimensional modeling method
CN108648265A (en) Helical gears gear hobbing process flank of tooth three-dimensional modeling method
CN109977530B (en) Three-dimensional modeling method for straight tooth face gear
Połowniak et al. Double enveloping worm gear modelling using CAD environment
CN113536472A (en) Method for calculating involute initial point of gear based on tooth profile normal method
Michalski et al. Modelling the tooth flanks of hobbed gears in the CAD environment

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210928

Assignee: Wenzhou Kehong Robot Co.,Ltd.

Assignor: Wenzhou University

Contract record no.: X2024330000007

Denomination of invention: A modeling method and measurement method for theoretical distribution of tooth surface point clouds

Granted publication date: 20230103

License type: Common License

Record date: 20240124