CN103292655B - A kind of computing method of the cylindrical acts having benchmark to retrain - Google Patents

A kind of computing method of the cylindrical acts having benchmark to retrain Download PDF

Info

Publication number
CN103292655B
CN103292655B CN201310250956.2A CN201310250956A CN103292655B CN 103292655 B CN103292655 B CN 103292655B CN 201310250956 A CN201310250956 A CN 201310250956A CN 103292655 B CN103292655 B CN 103292655B
Authority
CN
China
Prior art keywords
axis
cylinder
contact point
plane
jump
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.)
Active
Application number
CN201310250956.2A
Other languages
Chinese (zh)
Other versions
CN103292655A (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.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic 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 Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN201310250956.2A priority Critical patent/CN103292655B/en
Publication of CN103292655A publication Critical patent/CN103292655A/en
Application granted granted Critical
Publication of CN103292655B publication Critical patent/CN103292655B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

Computing method for the cylindrical acts having benchmark to retrain, first the method is measured and is obtained the measuring point coordinate on tested cylinder and benchmark; Then based on Minimum Area matching benchmark, provide the initial parameter of tested cylinder, and measuring point is projected in the plane perpendicular to measured circle mast axis, the method matching subpoint of the action face of computation circle; Inquire about the measuring point that tested cylinder contacts with action face, and according to the coordinate of contact point, carry out rotating variation or dimensional variations accordingly to action face; And then inquire about contact point, carry out the variation of action face, successively iterative computation, until meet criterion, export the parameter of tested cylinder.The present invention accurately can calculate the optimal value of measured circle column parameter.

Description

A kind of computing method of the cylindrical acts having benchmark to retrain
Technical field
The present invention relates to a kind of computing method of the cylindrical acts having benchmark to retrain, belong to delicate metering and computer application field, the cylindrical acts that can be used for having benchmark to retrain in geometric product in various situation, the qualification of any direction straight line degree of tilt detect, and provide guidance for the improvement of process and processing technology.
Background technology
Cylinder is one of modal element in mechanical component, and its precision is on quality, the performance of product and have important impact serviceable life, and cylindrical acts is the important technological parameters of cylindrical component.Maximum inscribed cylinder and Minimum circumscribed cylinder are referred to as cylindrical action face, and the radius of maximum inscribed cylinder and the radius of Minimum circumscribed cylinder are referred to as cylindrical acts.
The computing method of maximum inscribed cylinder and Minimum circumscribed cylinder are the computing method set up based on the Cleaning Principle of smooth cylindrical ring gauge.Maximum inscribed cylinder embody measured hole by maximum fitted shaft, thus obtained deviation from cylindrical form can be considered the maximal clearance between measured hole and maximum fitted shaft; And Minimum circumscribed cylinder embody measured axis by minimum mating holes, thus obtained deviation from cylindrical form can be considered the maximal clearance between measured axis and minimum mating holes.
The present invention mainly solves the calculating of the right cylinder acts of benchmark constraint.The right cylinder acts having benchmark to retrain, except the above using value introduced, the computation process of the minimum circumscribed circle column dimension having benchmark to retrain can be applied to the bank error calculating of any direction straight line.
The calculating of the right cylinder acts having benchmark to retrain, belongs to non-differentiability complex optimization problem, and at present, Chinese scholars mainly adopts traditional optimization, intelligent algorithm, computational geometry method etc.These methods owing to there is computational stability difference, counting yield is low, to adopting, a quantity is restricted, result of calculation is difficult to reach the defects such as exact solution, is difficult to application of putting into production.Ripe least square method is generally all adopted to calculate the right cylinder acts having benchmark to retrain approx in the market.
Summary of the invention
In order to overcome above-mentioned technical disadvantages, the object of this invention is to provide a kind of computing method of the cylindrical acts having benchmark to retrain.The evaluation of the cylindrical acts retrained for there being benchmark, any direction straight line degree of tilt, this method not only increases the accuracy of detection of surveying instrument, and algorithm stability is good, counting yield is high.
For solving the problem, the computing method of a kind of cylindrical acts having benchmark to retrain of the present invention, during for calculating cylindrical minimum external size, mainly comprise the following steps:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system.
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector .
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for minimum circumscribed circle, calculate the center of circle of minimum circumscribed circle o, by the center of circle oas the axis of tested cylinder on a bit .
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done.
Step 5: calculate each measuring point of tested cylinder successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to ultimate range, measuring point corresponding to the ultimate range recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of ultimate range.
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards overlapping, and determines to rotate variation initial value.
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7.
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely reduce the diameter of cylinder action face, determine dimensional variations amount initial value.
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9.
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations , ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step.
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.
In said method, described calculation process can be applied to the straight line bank error evaluation of any direction.
The computing method of cylindrical acts that the present invention's another kind has benchmark to retrain, for calculate most imperial palace connect size time, mainly comprise the following steps:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system.
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector .
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for maximum inscribed circle, calculate the center of circle of maximum inscribed circle o, by the center of circle oas the axis of tested cylinder on a bit .
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done.
Step 5: to calculate on tested cylinder each measuring point successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to minor increment, measuring point corresponding to the minor increment recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of minor increment.
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes cylinder axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards vertical, determines to rotate variation initial value.
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7.
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely increase the diameter of cylinder action face, determine dimensional variations amount initial value.
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9.
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations with ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step.
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.
The present invention, according to the cylindrical feature having benchmark to retrain, on the direction optimized, carries out corresponding rotation variation or dimensional variations to cylinder action face, finally calculates the optimal value of measured circle column parameter.
Beneficial effect of the present invention is: the relative position of the contact point that this method inquires meets criterion, the right cylinder acts having benchmark to retrain calculated or any direction straight line bank error are optimal value, and computational stability is good, counting yield is high.
To those skilled in the art, according to application design disclosed by the invention, easily can carry out various distortion and change to the present invention program, it should be noted that all these distortion and change all should belong to scope of the present invention.
Accompanying drawing explanation
Fig. 1 is the calculation flow chart of the minimum external size of the right cylinder having benchmark to retrain of the present invention.
Embodiment
Embodiment 1:
Computing method for the cylindrical acts having benchmark to retrain, during for calculating right cylinder minimum external size, comprise the steps, as shown in Figure 1:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system.
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector .
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for minimum circumscribed circle, calculate the center of circle of minimum circumscribed circle o, by the center of circle oas the axis of tested cylinder on a bit .
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done.
Step 5: calculate each measuring point of tested cylinder successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to ultimate range, measuring point corresponding to the ultimate range recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of ultimate range.
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards overlapping, and determines to rotate variation initial value.
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7.
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely reduce the diameter of cylinder action face, determine dimensional variations amount initial value.
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9.
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations , ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step.
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.
Embodiment 2:
A kind of computing method of the cylindrical acts having benchmark to retrain, during for calculating minimum external size, need to change embodiment 1, in embodiment 2, the rotation variation of action face and dimensional variations are the reverse variations rotating variation and dimensional variations in embodiment 1, and concrete steps are as follows:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system.
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector .
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for maximum inscribed circle, calculate the center of circle of maximum inscribed circle o, by the center of circle oas the axis of tested cylinder on a bit .
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done.
Step 5: to calculate on tested cylinder each measuring point successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to minor increment, measuring point corresponding to the minor increment recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of minor increment.
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes cylinder axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards vertical, determines to rotate variation initial value.
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7.
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely increase the diameter of cylinder action face, determine dimensional variations amount initial value.
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9.
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations with ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step.
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.

Claims (2)

1. computing method for the cylindrical acts having benchmark to retrain, for calculating minimum external size, it is characterized in that, concrete steps are as follows:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system;
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector ;
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for minimum circumscribed circle, calculate the center of circle of minimum circumscribed circle o, by the center of circle oas the axis of tested cylinder on a bit ;
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done;
Step 5: calculate each measuring point of tested cylinder successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to ultimate range, measuring point corresponding to the ultimate range recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of ultimate range;
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards overlapping, and determines to rotate variation initial value;
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7;
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely reduce the diameter of cylinder action face, determine dimensional variations amount initial value;
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9;
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations , ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step;
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.
2. computing method for the cylindrical acts having benchmark to retrain, connect size for calculating most imperial palace, it is characterized in that, concrete steps are as follows:
Step 1: be placed on measuring table by the part comprising tested cylinder, obtains the measuring point coordinate of reference characteristic and tested cylinder respectively in measurement space rectangular coordinate system;
Step 2: judgment standard feature is plane or straight line;
If reference characteristic is plane, then applies minimum area method matching reference plane, obtain the normal vector of reference plane; If reference characteristic is straight line, then applies minimum area method matching reference line, obtain the direction vector of reference line; The direction vector unification of the normal vector of reference plane or reference line is designated as base vector ;
Step 3: setup algorithm error permissible value ; Provide measured circle mast axis at random direction vector initial value, and to ensure with become the desired angle that drawing specifies ; The measuring point of tested cylinder is projected in perpendicular to plane on, in projection plane, by the method for maximum inscribed circle, calculate the center of circle of maximum inscribed circle o, by the center of circle oas the axis of tested cylinder on a bit ;
Step 4: carry out coordinate transform, make base vector be parallel to coordinate system z-axis, simultaneously measuring point, the axis of tested cylinder direction vector and axis on a bit also corresponding coordinate transform is done;
Step 5: to calculate on tested cylinder each measuring point successively to cylinder axis distance, and record each measuring point to cylinder axis the measuring point corresponding to minor increment, measuring point corresponding to the minor increment recorded forms the test points set contacted with cylinder action face, the now diameter of cylinder action face equal each measuring point to arrive 2 times of minor increment;
Step 6: judge whether the quantity of contact point in contact point set equals 2;
If the quantity of contact point is not equal to 2, jump to step 8;
If the quantity of contact point equals 2, cylinder axis need around base vector rotate variation, the direction of variation makes cylinder axis ? xoyprojection in plane and 2 contact points exist xoyprojection join line in plane is tending towards vertical, determines to rotate variation initial value;
Step 7: cylinder axis around base vector rotate axis is obtained after angle , calculate axis direction vector; Axis on a bit it is the middle point coordinate of 2 contact points; Calculating contact point arrives distance , and calculate each noncontact measuring point according to this and arrive distance ; If , represent that rotation amount is excessive, become ; If , represent that rotation amount is too small, become ;
If value be less than or equal to the error permissible value of setting , by axis direction vector assignment to axis direction vector , will point coordinate assignment is given point, and judge whether to meet criterion, jumps to step 5 if do not met, satisfied then jump to step 11; If value be greater than the error permissible value of setting , repeat step 7;
Step 8: judge in contact point set, whether contact point quantity equals 3;
If contact point quantity is not equal to 3, jump to step 5; If contact point quantity equals 3, carry out dimensional variations, namely increase the diameter of cylinder action face, determine dimensional variations amount initial value;
Step 9: set contact point coordinate as if the axis that after dimensional variations, tested cylinder is corresponding is , direction vector be , on be a bit , ? xoyplane inner projection with xthe angle of axle is ; Each contact point all meets following formula,
Thereby establish ternary Nonlinear System of Equations, wherein unknown quantity is , ; By solving Nonlinear System of Equations, can obtain direction vector ; Calculate each noncontact measuring point according to this to axis distance ; If , represent that variation is excessive, become ; If , represent that variation is too small, become ;
If value be less than or equal to the error permissible value of setting , will assignment is given , will assignment is given , and judge whether to meet criterion, if met, jump to step 11, if do not met, jump to step 5; If value be greater than the error permissible value of setting , repeat step 9;
Step 10: judge whether contact point quantity is more than or equal to 4;
If contact point quantity is less than 4, then jump to step 5;
If contact point quantity is more than or equal to 4, be 1 combination with 3 contact points, with wherein 1 combination in 3 contact points for calculating object, carry out dimensional variations according to the method for step 9, driven dimension variation , and upgrade after dimensional variations with ;
If dimensional variations amount be greater than the error permissible value of setting , jump to step 5; If dimensional variations amount be less than or equal to the error permissible value of setting , then judge whether to meet criterion, if do not meet criterion, then change a combination, recalculate dimensional variations amount, repeat step 10, if meet criterion, then jump to next step;
Step 11: calculate and export deviation from cylindrical form, cylinder acts tand cylinder axis parameter.
CN201310250956.2A 2013-06-24 2013-06-24 A kind of computing method of the cylindrical acts having benchmark to retrain Active CN103292655B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310250956.2A CN103292655B (en) 2013-06-24 2013-06-24 A kind of computing method of the cylindrical acts having benchmark to retrain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310250956.2A CN103292655B (en) 2013-06-24 2013-06-24 A kind of computing method of the cylindrical acts having benchmark to retrain

Publications (2)

Publication Number Publication Date
CN103292655A CN103292655A (en) 2013-09-11
CN103292655B true CN103292655B (en) 2015-12-02

Family

ID=49093926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310250956.2A Active CN103292655B (en) 2013-06-24 2013-06-24 A kind of computing method of the cylindrical acts having benchmark to retrain

Country Status (1)

Country Link
CN (1) CN103292655B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106289145B (en) * 2016-10-27 2018-10-09 北京理工大学 A kind of cylindricity On-machine Test method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1305088A (en) * 1999-12-14 2001-07-25 大地公司 Method of monitoring the diameter for cylindrical body made by injection process
RU2190829C2 (en) * 1999-12-22 2002-10-10 Самарский государственный технический университет Process of ball control
US6729936B1 (en) * 1999-06-25 2004-05-04 Toyoda Koki Kabushiki Kaisha Apparatus for measuring dimensional errors of eccentric cylinder by utilizing movement of measuring member held in contact with such eccentric cylinder
CN102445137A (en) * 2011-09-20 2012-05-09 飞翼股份有限公司 Device for measuring radius of spherical crown and measuring method thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6729936B1 (en) * 1999-06-25 2004-05-04 Toyoda Koki Kabushiki Kaisha Apparatus for measuring dimensional errors of eccentric cylinder by utilizing movement of measuring member held in contact with such eccentric cylinder
CN1305088A (en) * 1999-12-14 2001-07-25 大地公司 Method of monitoring the diameter for cylindrical body made by injection process
RU2190829C2 (en) * 1999-12-22 2002-10-10 Самарский государственный технический университет Process of ball control
CN102445137A (en) * 2011-09-20 2012-05-09 飞翼股份有限公司 Device for measuring radius of spherical crown and measuring method thereof

Also Published As

Publication number Publication date
CN103292655A (en) 2013-09-11

Similar Documents

Publication Publication Date Title
CN103256916B (en) Evaluation method of part flatness error based on minimum area
CN108871256B (en) Roundness error evaluation algorithm
CN102661699A (en) Precision measurement method of large-radius short arc part
CN103292773B (en) A kind of symmetry error evaluation method based on Minimum Area
CN103486984B (en) The detection method of profile right alignment in a kind of wind-tunnel
CN103278126B (en) A kind of sphericity error assessment method for part based on Minimum Area
CN107330142B (en) On-machine detection measuring point normal vector estimation method based on STL model
CN103292769B (en) A kind of plane inclination error detection method based on Minimum Area
CN102944163A (en) Device and method for measuring profile tolerance of annular dovetail groove of any axial section
CN103294920B (en) A kind of part Cylindricity error evaluation based on Minimum Area
WO2021128614A1 (en) Method for measuring and evaluating error of feature line-based arc cam profile
CN104197823B (en) The measuring method of radio frequency quadrupole field plate cartridge three-dimensional space curved surface
CN103471572A (en) Total station networking measurement method of large-scale structural component
CN103292654B (en) A kind of method calculating function size of cylindrical part
CN103292655B (en) A kind of computing method of the cylindrical acts having benchmark to retrain
CN111765902A (en) Laser point cloud precision evaluation method based on polygonal pyramid target
CN103310106B (en) The computing method of a kind of hole in piece part system acts
US10132623B2 (en) Method for measuring slant wall thickness dimension of hub
CN102564383B (en) Method for determining interaxial angles of triaxial magnetic sensor based on axial surface angles
CN102095366A (en) Method for measuring large-gradient aspheric surface profile
CN106646413B (en) Radar networking vertical line cross fusion positioning method and error resolving method
CN103292770B (en) A kind of method calculating function size of cone part
CN205373636U (en) A mark piece for image and laser survey sensor data fusion
CN111366125B (en) Method for measuring rotation angle of crawler crane
CN114018202A (en) Novel algorithm for rapidly evaluating roundness

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Chen Leilei

Inventor after: Huang Meifa

Inventor after: Zhong Yanru

Inventor after: Xu Bensheng

Inventor after: Fan Kai

Inventor after: Sun Yingchao

Inventor after: Bao Jiading

Inventor before: Chen Leilei

Inventor before: Zhong Yanru

Inventor before: Sun Yonghou

Inventor before: Huang Meifa

Inventor before: Kuang Bing

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20151029

Address after: 541004 East District School of mechanical and electrical engineering, Guilin University of Electronic Technology, 1 Jinji Road, Qixing District, Guilin, the Guangxi Zhuang Autonomous Region

Applicant after: Guilin University of Electronic Technology

Address before: 541004 graduate school, Guilin University of Electronic Technology, 1 Jinji Road, Qixing District, Guilin, the Guangxi Zhuang Autonomous Region

Applicant before: Chen Leilei

C14 Grant of patent or utility model
GR01 Patent grant