CN110132195A - The probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method - Google Patents
The probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method Download PDFInfo
- Publication number
- CN110132195A CN110132195A CN201910461858.0A CN201910461858A CN110132195A CN 110132195 A CN110132195 A CN 110132195A CN 201910461858 A CN201910461858 A CN 201910461858A CN 110132195 A CN110132195 A CN 110132195A
- Authority
- CN
- China
- Prior art keywords
- vector
- point
- blade profile
- nurbs curve
- discrete
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
Abstract
The probe that the present invention provides a kind of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method, include the following steps: using each survey ball's dead center M of three-dimensional coordinates measurement as the data point P of nurbs curve, the knot vector K of nurbs curve is calculated using accumulation chord length method, and calculates the basic function matrix N of nurbs curvei,m;According to offset dot matrix control of reverse computing vertex matrix D, the nurbs curve of fitting is made to cross data point;Equidistant discrete nurbs curve calculates the normal vector N that discrete point is parallel to X/Y planexy;Discrete point O is O' along the subpoint of Z-direction plane where blade profile, calculates vectorEnable vectorIt is equal to N in the direction of X/Y planexyDirection, vectorIt is 0 in the component of Z-direction, thenWithFormed compensation vectorAs compensate three-dimensional radius compensation vector.The present invention solves blade profile and scans the cosine error occurred in radius compensation.
Description
Technical field
The present invention relates to three coordinate measuring engine measurement technical field, especially a kind of blade profile contact type scanning measurement
Probe surveys ball three-dimensional radius compensation method.
Background technique
Blade has important application, while processing and inspection to blade in aero-engine, turbine, field of wind power generation
Required precision is surveyed increasingly to increase.Contact and contactless two kinds of detection methods are mainly used in current vane detection, the former
Detection accuracy has tens microns of error compared with the latter, can only or big workpiece application low to some accuracy request classes.To high-precision
Crop leaf measuring, three coordinate measuring machine are first-selected instruments, during the scanning process to reduce stylus abrasion, are increased in stylus tip anti-
The ruby ball of abrasion as surveying ball, due to three coordinate measuring machine read be ruby ball center coordinate value, ground reducing
The problem of compensation of the ruby radius of a ball is brought while damage.
Current gauge head measuring needle radius compensation method is broadly divided into curved surface compensation and curve compensation.Curved surface compensation needs to measure leaf
Body numerical value, then the inwardly biasing ruby radius of a ball, this is just needed using five axis gauge heads, and sweeping measurement, this measurement method is not only
Data volume is big, and gauge head itself is expensive, so crop leaf measuring majority uses section line Scanning Detction at present, and carries out curve
Compensation.Curve compensation method is a kind of two dimension compensation, does not account for the cosine angle bring cosine error of blade twist generation.Needle
To this problem, patent CN 104330068A proposes a kind of use cross-sectional scans mode, measures multiple section lines and is fitted measurement song
Face according still further to the method for surface method arrow biasing, then intercepts section line after biasing.The method theoretically can reduce cosine error,
But need to measure multiple sections, measurement efficiency is difficult to ensure.
Summary of the invention
It is an object of the present invention to overcome the shortcomings of the prior art and provide a kind of surveys of blade profile contact type scanning
The probe of amount surveys ball three-dimensional radius compensation method, scans the cosine error occurred in radius compensation to solve blade profile.This
Invention the technical solution adopted is that:
A kind of probe survey ball three-dimensional radius compensation method of blade profile contact type scanning measurement, includes the following steps:
Step S101, coordinate system XYZ is used by this method of blade coordinate system to be measured, by each survey ball of three-dimensional coordinates measurement
Data point P of the central point M as nurbs curve calculates the knot vector K of nurbs curve using accumulation chord length method, and calculates
The basic function matrix N of nurbs curvei,m;According to offset dot matrix control of reverse computing vertex matrix D, make the nurbs curve mistake of fitting
Data point;
Step S102, equidistant discrete nurbs curve, each discrete point are expressed as point O, calculate discrete point and are parallel to X/Y plane
Normal vector Nxy;
Step S103, discrete point O is O' along the subpoint of Z-direction plane where blade profile, calculates vector
A point is to survey ball and the practical contact point of blade, enables vectorIt is equal to N in the direction of X/Y planexyDirection, vectorIn Z axis
The component in direction is 0, thenWithFormed compensation vectorAs compensate three-dimensional radius compensation vector.
Further, step S101 is specifically included:
Coordinate system XYZ is used by this method of blade coordinate system to be measured, blade profile to be measured is parallel to plane where X/Y axis
That is X/Y plane, plane is perpendicular to Z axis where blade profile to be measured;
Using each survey ball's dead center M of three-dimensional coordinates measurement as data point P, n+1 data point is shared, nurbs curve
Number m;
The equation C of nurbs curve is expressed are as follows:
Wherein, k is variable, diFor control vertex, wiFor weight, K=[k0,…,ki,…,kn+m+1] it is nurbs curve
Knot vector can be expressed as:
PiFor data point coordinate;
In formula (1), Ni,mFor the basic function of nurbs curve, it is expressed as the coefficient matrix of formula (3) or formula (5);
Control vertex matrix D is expressed by the control vertex in formula (1) are as follows:
D=[d0 … dn]T. (4)
Formula (4) is solved by formula (5), and the coefficient matrix of formula (5) is a kind of easy expression of formula (3);
Coefficient matrix in formula (5):
For deformation matrix of the offset dot matrix under knot vector control;
Wherein,。
Further, step S102 is specifically included:
Equidistant discrete nurbs curve, l point of covariance, discrete distance determine that discrete point is denoted as C by discrete point quantityi
(k), [1, l] i ∈, each discrete point are expressed as point O, and all discrete points project to blade profile place plane along Z-direction, in blade
Section planar calculates the tangent vector of each discrete point:
Ti(k)=Ci'(k). (6)
Wherein, Ci' (k) be Ci(k) derivative is then parallel to the normal vector N of X/Y planexyExpression are as follows:
Normal vector NxyIt is expressed as Nxy=(nx,ny), nxIndicate component of the normal vector in X-direction, nyIndicate normal direction arrow
Measure the component in Y direction.
Further, in step S103, compensation vectorExpression are as follows:
Wherein, z is the Z-direction numerical value of discrete point, z0By survey section Z-direction numerical value.
Advantages of the present invention: the probe of blade profile contact type scanning measurement proposed by the present invention surveys ball three-dimensional radius compensation
Method is compensated using ball sphere center position and section a difference in height existing for Z-direction is surveyed when probe scanning section in two dimension
On the basis of increase this elevation information, improve blade profile scanning in radius compensation precision.
Detailed description of the invention
Fig. 1 is flow chart of the invention.
Fig. 2 is compensation vector schematic diagram of the invention.
Specific embodiment
Below with reference to specific drawings and examples, the invention will be further described.
The present invention proposes that a kind of probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method, including such as
Lower step:
Step S101, coordinate system XYZ is used by this method of blade coordinate system to be measured, by each survey ball of three-dimensional coordinates measurement
Central point M is as NURBS (non-uniform rational B-spline, Non-Uniform Rational B-Splines)) data point of curve
P calculates the knot vector K of nurbs curve using accumulation chord length method, and calculates the basic function matrix N of nurbs curvei,m;According to
Offset dot matrix control of reverse computing vertex matrix D makes the nurbs curve of fitting cross data point;
Specifically,
Coordinate system XYZ is used by this method of blade coordinate system to be measured, blade profile to be measured is parallel to plane where X/Y axis
That is X/Y plane, plane is perpendicular to Z axis where blade profile to be measured;
Using each survey ball's dead center M of three-dimensional coordinates measurement as data point P, n+1 data point is shared, the present invention uses 3
The number m of secondary nurbs curve, i.e. nurbs curve takes 3;
The equation C of nurbs curve is expressed are as follows:
Wherein, k is variable, diFor control vertex, wiFor weight, weight of the present invention takes 1, K=[k0,…,ki,…,kn+m+1]
For the knot vector of nurbs curve, can be expressed as:
PiFor data point coordinate;
In formula (1), Ni,mFor the basic function of nurbs curve, it can be expressed as formula (3), easier expression way is
The coefficient matrix of formula (5);
Control vertex matrix D is expressed by the control vertex in formula (1) are as follows:
D=[d0 … dn]T. (4)
Formula (4) is solved by formula (5), and the coefficient matrix of formula (5) is a kind of easy expression of formula (3);
Coefficient matrix in formula (5):
For deformation matrix of the offset dot matrix under knot vector control;
Wherein,
Step S102, equidistant discrete nurbs curve, each discrete point are expressed as point O, calculate discrete point and are parallel to X/Y plane
Normal vector Nxy;
It specifically includes:
Equidistant discrete nurbs curve, l point of covariance, discrete distance determine that discrete point is denoted as C by discrete point quantityi
(k), [1, l] i ∈, each discrete point are expressed as point O, and all discrete points project to blade profile place plane along Z-direction, in blade
Section planar calculates the tangent vector of each discrete point:
Ti(k)=Ci'(k). (6)
Wherein, Ci' (k) be Ci(k) derivative is then parallel to the normal vector N of X/Y planexyIt can be expressed as:
Normal vector NxyIt is expressed as Nxy=(nx,ny), nxIndicate component of the normal vector in X-direction, nyIndicate normal direction arrow
Measure the component in Y direction;
Step S103, point O is O' in the subpoint of X/Y plane along Z-direction, calculates vectorA point is to survey ball and leaf
The practical contact point of piece, enables vectorIt is equal to N in the direction of X/Y planexyDirection, vectorIt is in the component of Z-direction
0, thenWithFormed compensation vectorAs compensate three-dimensional radius compensation vector.As shown in Figure 2.
Step S103, discrete point O is O' along the subpoint of Z-direction plane where blade profile, calculates vector
A point is to survey ball and the practical contact point of blade, enables vectorIt is equal to N in the direction of X/Y planexyDirection, vectorIn Z axis
The component in direction is 0, thenWithFormed compensation vectorAs compensate three-dimensional radius compensation vector.Such as Fig. 2 institute
Show.
It specifically includes:
Compensation vectorExpression are as follows:
Wherein, z is the Z-direction numerical value of discrete point, z0By survey section Z-direction numerical value.
It should be noted last that the above specific embodiment is only used to illustrate the technical scheme of the present invention and not to limit it,
Although being described the invention in detail referring to example, those skilled in the art should understand that, it can be to the present invention
Technical solution be modified or replaced equivalently, without departing from the spirit and scope of the technical solution of the present invention, should all cover
In the scope of the claims of the present invention.
Claims (4)
1. a kind of probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method, which is characterized in that including such as
Lower step:
Step S101, coordinate system XYZ is used by this method of blade coordinate system to be measured, by each survey ball center of three-dimensional coordinates measurement
Data point P of the point M as nurbs curve calculates the knot vector K of nurbs curve using accumulation chord length method, and calculates NURBS
The basic function matrix N of curvei,m;According to offset dot matrix control of reverse computing vertex matrix D, the nurbs curve of fitting is made to cross offset
Point;
Step S102, equidistant discrete nurbs curve, each discrete point are expressed as point O, calculate the method that discrete point is parallel to X/Y plane
To vector Nxy;
Step S103, discrete point O is O' along the subpoint of Z-direction plane where blade profile, calculates vectorA point is
Ball and the practical contact point of blade are surveyed, vector is enabledIt is equal to N in the direction of X/Y planexyDirection, vectorIn Z-direction
Component is 0, thenWithFormed compensation vectorAs compensate three-dimensional radius compensation vector.
2. the probe of blade profile contact type scanning measurement as described in claim 1 surveys ball three-dimensional radius compensation method, special
Sign is,
Step S101 is specifically included:
Coordinate system XYZ is used by this method of blade coordinate system to be measured, blade profile to be measured is parallel to plane i.e. XY where X/Y axis
Plane, plane is perpendicular to Z axis where blade profile to be measured;
Using each survey ball's dead center M of three-dimensional coordinates measurement as data point P, n+1 data point, the number of nurbs curve are shared
m;
The equation C of nurbs curve is expressed are as follows:
Wherein, k is variable, diFor control vertex, wiFor weight, K=[k0,…,ki,…,kn+m+1] be nurbs curve node
Vector can be expressed as:
PiFor data point coordinate;
In formula (1), Ni,mFor the basic function of nurbs curve, it is expressed as the coefficient matrix of formula (3) or formula (5);
Control vertex matrix D is expressed by the control vertex in formula (1) are as follows:
D=[d0 … dn]T. (4)
Formula (4) is solved by formula (5), and the coefficient matrix of formula (5) is a kind of easy expression of formula (3);
Coefficient matrix in formula (5):
For deformation matrix of the offset dot matrix under knot vector control;
Wherein,
3. the probe of blade profile contact type scanning measurement as claimed in claim 2 surveys ball three-dimensional radius compensation method, special
Sign is,
Step S102 is specifically included:
Equidistant discrete nurbs curve, l point of covariance, discrete distance determine that discrete point is denoted as C by discrete point quantityi(k), i
∈ [1, l], each discrete point are expressed as point O, and all discrete points project to blade profile place plane along Z-direction, in blade profile
Planar calculate the tangent vector of each discrete point:
Ti(k)=Ci'(k). (6)
Wherein, Ci' (k) be Ci(k) derivative is then parallel to the normal vector N of X/Y planexyExpression are as follows:
Normal vector NxyIt is expressed as Nxy=(nx,ny), nxIndicate component of the normal vector in X-direction, nyIndicate that normal vector exists
The component of Y direction.
4. the probe of blade profile contact type scanning measurement as claimed in claim 3 surveys ball three-dimensional radius compensation method, special
Sign is,
In step S103, compensation vectorExpression are as follows:
Wherein, z is the Z-direction numerical value of discrete point, z0By survey section Z-direction numerical value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910461858.0A CN110132195B (en) | 2019-05-30 | 2019-05-30 | Three-dimensional radius compensation method for probe ball measurement of blade section contact type scanning measurement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910461858.0A CN110132195B (en) | 2019-05-30 | 2019-05-30 | Three-dimensional radius compensation method for probe ball measurement of blade section contact type scanning measurement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110132195A true CN110132195A (en) | 2019-08-16 |
CN110132195B CN110132195B (en) | 2021-04-20 |
Family
ID=67582874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910461858.0A Active CN110132195B (en) | 2019-05-30 | 2019-05-30 | Three-dimensional radius compensation method for probe ball measurement of blade section contact type scanning measurement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110132195B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112033338A (en) * | 2020-09-14 | 2020-12-04 | 中国航空工业集团公司北京长城计量测试技术研究所 | Blade curved surface contact type scanning measurement probe radius surface compensation method |
CN117606327A (en) * | 2024-01-24 | 2024-02-27 | 天津工业大学 | Aeroengine blade profile detection equipment |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4612622A (en) * | 1983-07-12 | 1986-09-16 | United Technologies Corporation | Probe for coordinate measuring machine |
CN103777570A (en) * | 2014-01-07 | 2014-05-07 | 浙江大学 | Machining error rapid detection and compensation method based on NURBS curved surface |
CN106202822A (en) * | 2016-07-27 | 2016-12-07 | 西北工业大学 | B-spline surface model reconstruction method towards blade adaptive machining |
CN106393106A (en) * | 2016-09-13 | 2017-02-15 | 东南大学 | Parameter adapting and calibrating robot NURBS curvilinear motion interpolation method |
CN107238364A (en) * | 2017-06-30 | 2017-10-10 | 四川大学 | Contact type measurement chaining pin Probe-radius fine compensation method |
CN108106522A (en) * | 2017-11-29 | 2018-06-01 | 中国航发沈阳黎明航空发动机有限责任公司 | A kind of method for three-dimensional measurement of irregular surface |
CN108692644A (en) * | 2018-03-26 | 2018-10-23 | 华中科技大学 | A kind of complex-curved Three-coordinate measurer and error compensating method |
CN109341632A (en) * | 2018-11-29 | 2019-02-15 | 株洲中航动力精密铸造有限公司 | Gauge head and three-coordinates measuring machine for turbo blade blade profile measurement |
-
2019
- 2019-05-30 CN CN201910461858.0A patent/CN110132195B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4612622A (en) * | 1983-07-12 | 1986-09-16 | United Technologies Corporation | Probe for coordinate measuring machine |
CN103777570A (en) * | 2014-01-07 | 2014-05-07 | 浙江大学 | Machining error rapid detection and compensation method based on NURBS curved surface |
CN106202822A (en) * | 2016-07-27 | 2016-12-07 | 西北工业大学 | B-spline surface model reconstruction method towards blade adaptive machining |
CN106393106A (en) * | 2016-09-13 | 2017-02-15 | 东南大学 | Parameter adapting and calibrating robot NURBS curvilinear motion interpolation method |
CN107238364A (en) * | 2017-06-30 | 2017-10-10 | 四川大学 | Contact type measurement chaining pin Probe-radius fine compensation method |
CN108106522A (en) * | 2017-11-29 | 2018-06-01 | 中国航发沈阳黎明航空发动机有限责任公司 | A kind of method for three-dimensional measurement of irregular surface |
CN108692644A (en) * | 2018-03-26 | 2018-10-23 | 华中科技大学 | A kind of complex-curved Three-coordinate measurer and error compensating method |
CN109341632A (en) * | 2018-11-29 | 2019-02-15 | 株洲中航动力精密铸造有限公司 | Gauge head and three-coordinates measuring machine for turbo blade blade profile measurement |
Non-Patent Citations (2)
Title |
---|
刘荷辉等: "自由曲面的二维自适应测量及测球半径的三维补偿", 《机械工程学报》 * |
王建华等: "复杂型面测量中测头中心的轨迹曲面及测头半径的三维补偿", 《计量学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112033338A (en) * | 2020-09-14 | 2020-12-04 | 中国航空工业集团公司北京长城计量测试技术研究所 | Blade curved surface contact type scanning measurement probe radius surface compensation method |
CN117606327A (en) * | 2024-01-24 | 2024-02-27 | 天津工业大学 | Aeroengine blade profile detection equipment |
CN117606327B (en) * | 2024-01-24 | 2024-04-05 | 天津工业大学 | Aeroengine blade profile detection equipment |
Also Published As
Publication number | Publication date |
---|---|
CN110132195B (en) | 2021-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105627923B (en) | A kind of planning parameters of scanning paths method that blade dense point cloud based on laser ranging is obtained | |
CN109115095B (en) | Structural parameter optimization method of non-contact R-test measuring instrument | |
CN103292732B (en) | A kind of large-scale free form surface on-machine measurement device of extension type | |
CN101387494A (en) | Geometrical dimensional measurement apparatus and method for large-sized tunnel tunnel segment component | |
CN102873586B (en) | Fast on-line measuring device for curvature radius of workpiece processed in numerically controlled manner | |
CN104070418B (en) | Axial symmetry surface shape of optical aspheric surface On-line Measuring Method | |
CN103438828B (en) | A kind of method of laser detection of screw rotor end sectional shape | |
CN110132195A (en) | The probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method | |
CN107990856B (en) | Spatial position error detection method for over-range workpiece | |
CN102944204A (en) | Method for detecting cross-section profile by utilizing two-dimensional measuring machine | |
CN103009194B (en) | Non-contact inner parallel plane distance measuring method for large-sized workpiece | |
CN108344361A (en) | A kind of planar process vector quantity measurement method based on laser tracker | |
CN111536877A (en) | Method for calibrating attitude of line laser sensor on three-coordinate measuring machine | |
CN105519264B (en) | The method of contactless composite measurement large flat leaky antenna three-dimensional surface shape and gap size and system thereof | |
CN108528755A (en) | A kind of digitlization installation method of acceleration transducer | |
Liu et al. | A coaxiality measurement model based on multi-systematic errors separation for turbine shaft | |
CN102749028A (en) | Method for measuring intermediate diameter of taper thread gauge of universal tool microscope in contact way | |
CN201885667U (en) | Skirt sleeve type line measuring tool for pistons | |
CN203274694U (en) | Telescopic large free-form surface on machine measurement device | |
CN208269868U (en) | The measuring device of high-precision industrial robot repetitive positioning accuracy | |
CN110186405B (en) | Blade section contact type scanning probe sphere measuring three-dimensional radius compensation and cross compensation point correction method | |
CN104880156A (en) | Method for measuring large R value by using CAD technology | |
CN209386967U (en) | A kind of device of auto-measuring structural elements geometry initial imperfection | |
CN206756128U (en) | A kind of joint arm measuring force error identification device | |
CN112595280A (en) | Method for measuring angled complex surface shape |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190816 Assignee: Jiangsu Jihui Huake Intelligent Equipment Technology Co., Ltd. Assignor: Wuxi research institute of the Central China University of Science and Technology Contract record no.: X2019980000288 Denomination of invention: The probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method License type: Common License Record date: 20191023 |
|
GR01 | Patent grant | ||
GR01 | Patent grant |