CN106897501A - The positioning and optimizing method based on blade parts deformation towards in adaptive machining - Google Patents
The positioning and optimizing method based on blade parts deformation towards in adaptive machining Download PDFInfo
- Publication number
- CN106897501A CN106897501A CN201710057686.1A CN201710057686A CN106897501A CN 106897501 A CN106897501 A CN 106897501A CN 201710057686 A CN201710057686 A CN 201710057686A CN 106897501 A CN106897501 A CN 106897501A
- Authority
- CN
- China
- Prior art keywords
- section line
- design section
- line
- sectional plane
- blade
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Numerical Control (AREA)
Abstract
Present invention proposition is a kind of towards the positioning and optimizing method based on blade parts deformation in adaptive machining, and the n bars first in blade design model obtain detection sectional plane line along the design section line of the parallel distribution in long-pending folded axle directionParameter expression, next sets up the multiple objective function of positioning and optimizing, and constraint when carrying out multiple-objection optimization solution is set up, and solution is finally optimized, find the purpose of target processing curve in preform to reach by the shape of mathematical optimization models section line.The present invention finds target processing curve with the non-uniform profile degree tolerance in blade type face as geometry constraint conditions in constraint, so as to meet follow-up blade profile detection require, improves the qualification rate of product.And positioning and optimizing method proposed by the present invention allows design section line to be deformed in constraint, while ensure that detection sectional plane wire shaped is closest with design section line, improve search out can processing curve ability.
Description
Technical field
It is specifically a kind of based on non-equal the present invention relates to the positioning and optimizing field of adaptive machining Leaf class part
The blade profile line deformation positioning and optimizing method of even profile tolerance tolerance constraints.
Background technology
In recent years, near-net-shape+accurate digital control technique realizes that the adaptive machining of the preforming blade of small surplus turns into main
The measurement point that preforming blade after the casting of trend, i.e. finish forge/essence is obtained by on-machine measurement with design a model that to do optimization registering, according to
Nc program is compensated according to the positioning posture of optimization rear blade.But because blade space is complex-shaped, near-net-shape technique
Blade afterwards there may be mass defect causes adaptive machining position fixing process to can not find target processing curve and turn into waste product, leaf
Piece, as the expensive high precision part produced in enormous quantities and the cycle is more long, is a huge loss for enterprise.Therefore near
The target processing curve as programming foundation is found in preform after net forming technology as far as possible, the conjunction to improving blade
Lattice rate has very important meaning.
Document " blade processing surplus adaptive optimization method [J] computer technologies and development based on on-line checking,
2014,24(11):226-229. " establishes allowance adaptive optimization for the uneven problem of blade blank allowance balance
Uniform mathematical model, realizes being accurately positioned for blade.This method thinks to change weight with machining benchmark when the design basis of blade
After conjunction, design a model and the European conversion in space is done in preform to meet allowance for finish, but finding best orientation
The profile tolerance tolerance of blade profile line is not considered during attitude, causes the blade after finishing to exceed when blade profile is detected
Form tolerance and turn into waste product;Furthermore, current positioning and optimizing method thinks that workpiece is rigid, does not recognize cutting for blade
Upper thread can deform in its profile tolerance tolerance constraints region, can cause to can not find target processing song for the blade of small surplus
Face and be mistakenly considered " waste product ", therefore current positioning and optimizing method not fully adds suitable for the self adaptation of near-net-shape class blade
Work.
The content of the invention
Machined surface profile degree tolerance is not considered during blade adaptive machining positioning and optimizing about in order to solve existing method
Beam and searching low the two problems of target processing curve ability.The present invention propose it is a kind of towards in adaptive machining based on leaf
The positioning and optimizing method of piece class part deformation, can not only meet the form tolerance requirement of subsequent detection blade profile, but also
The blade that current localization method is considered " scrapping " can be saved.
Blade is that by family's design section line traffic control, design section line can be divided into leaf basin, blade back, preceding according to function difference
Edge and the subregion of trailing edge four, the present invention with the non-uniform profile degree tolerance of design section line as geometry constraint conditions, by excellent
The shape for changing the section line that designs a model finds the purpose of target processing curve to reach in preform.
The technical solution adopted for the present invention to solve the technical problems is:
It is described a kind of towards the positioning and optimizing method based on blade parts deformation in adaptive machining, it is characterised in that:
Comprise the following steps:
Step 1:N bars in blade design model along the parallel distribution in long-pending folded axle direction design section lineDetection sectional plane line is obtained using following stepsParameter expression, l=1,2 ... n;
Step 1.1:Design section line is represented using cubic B-spline parametric equation:
Wherein Ni,3T () is the basic function of B-spline Curve, ViIt is the coordinate of the control vertex of B-spline curves, m0Represent
The quantity of control vertex;
Step 1.2:Design section line is around reference axis x, the postrotational parameter expression of y, z:
Wherein Rl=Rx(αl)·Ry(θl)·Rz(βl) it is the l articles spin matrix of design section line, Rx(αl) it is the l articles
Design section line is around x-axis rotation alphalThe spin matrix at angle, Ry(θl) θ is rotated around y-axis for the l bars design section linelThe spin moment at angle
Battle array, Rz(βl) for the l bars design section line around z-axis rotation βlThe spin matrix at angle;
Step 1.3:Postrotational design section line base vector is obtained according to the spin matrix of design section line;Then adjust
CurvePosition of the control vertex on base vector direction obtain deformed sections line
The wherein the l articles shift in position amount of i-th control vertex of design section line is { Δ x 'i,l,Δy′i,l};
Step 1.4:By deformed sections lineSetting-out generates target processing curve s (u, v);Target is processed
The parametric equation of curved surface s (u, v) is
Wherein Ni,3(u),Nj,3V () is respectively basic function of the Bicubic B-Spline Surfaces on u and v directions, n1It is u direction
The number of upper control vertex, m1It is the number of control vertex on v directions, Vi,jIt is the control vertex coordinate of target processing curve;
Step 1.5:For the l bars design section line, according to its z directions height zl, calculate s (u, v)=zlWhen friendship
Point parameterN0Represent the number of intersection point parameter;WillSubstitute into mesh
Mark processing curve parametric equation, obtains the coordinate of discrete pointAccording to discrete point coordinates, utilize
Interpolation method obtains the l bars detection sectional plane lineParameter expression;
Step 2:According to the detection sectional plane line that step 1 is obtainedParameter expression, set up positioning it is excellent
The multiple objective function of change, optimization aim is that detection sectional plane line is closest with corresponding design section wire shaped;
Step 3:Foundation carries out constraint when multiple-objection optimization is solved:Including
Constraint 1, detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four are within upper profile tolerance tolerance range
Constraint function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four
In upper profile tolerance tolerance rangePoint;
Constraint 2, detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four are outside bottom profiled degree tolerance range
Constraint function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four
In bottom profiled degree tolerance rangePoint;Represent in the l articles profile tolerance tolerance range o'clock the l articles inspection
Survey corresponding closest approach on section line;Represent the list of corresponding closest approach on the l bars detection sectional plane line
The outer method arrow in position;I1=1,2 ... N1, i2=1,2 ... N2, i3=1,2 ... N3, i4=1,2 ... N4, N1,N2,N3,N4Respectively section line
The quantity of the point in leading edge, trailing edge, leaf basin and blade back region;
Constraint 3, constraint function F of the blade blank measurement point on the outside of target processing curve2:
F2=(pi-s(ui,vi))·ni>=0 i=1,2 ... N
Wherein piRepresent that on-machine measurement obtains the ith measurement point of the blade blank after near-net-shape technique, N is measurement point
Quantity;s(ui,vi) it is closest approach of the ith measurement point on target processing curve, niRepresent closest approach s (ui,vi) in target
The outer method arrow of unit on processing curve;
Without flex point constraint in the line segment of constraint 4, detection sectional plane line;
Step 4:With design section lineAround x, y, the z-axis anglec of rotation { αl,θl,βl, l=1,2 ... n
And variation { the Δ x ' of control vertexi,l,Δy′i,l, i=1,2 ... m0It is optimized variable, the multiple target set up to step 2
Function optimizes solution, the design section line after being optimized, and obtains target processing curve s (u, v).
Further preferred scheme, it is described a kind of towards the positioning and optimizing based on blade parts deformation in adaptive machining
Method, it is characterised in that:The profile tolerance tolerance range of the l bars design section line is
N (t)=t (t) × ez
WhereinIt is the l bars upper profile tolerance band of design section line,It is upper profile variation amount, n (t) is design
Section lineIn the outer method arrow of the unit of parameter t,It is the l articles bottom profiled tolerance range of design section line,It is bottom profiled
Departure;T (t) is design section lineArrow, e are cut in the unit of parameter tzRepresent the unit vector in z directions;To set
Meter section lineArrow is cut in the single order of parameter t.
Further preferred scheme, it is described a kind of towards the positioning and optimizing based on blade parts deformation in adaptive machining
Method, it is characterised in that:It is constrained to without flex point in the line segment of detection sectional plane line:The controlling polygon convex-concave of detection sectional plane line line segment
Property is consistent.
Further preferred scheme, it is described a kind of towards the positioning and optimizing based on blade parts deformation in adaptive machining
Method, it is characterised in that:Multiple objective function is:
WhereinRepresent that the single order of the l bars design section line cuts arrow,Represent the l articles single order of detection sectional plane line
Cut arrow;The l articles design section curvature of a curve is represented,The l articles detection sectional plane curvature of a curve is represented,Represent that the l articles sets
The second order for counting section line cuts arrow,Represent that the second order of the l bars detection sectional plane line cuts arrow.
Beneficial effect
The present invention has the following advantages that compared to traditional surplus optimization method:
1. the present invention finds target with the non-uniform profile degree tolerance in blade type face as geometry constraint conditions in constraint
Processing curve, so as to meet follow-up blade profile detection require, improves the qualification rate of product.
2., for the blade of small surplus after near-net-shape technique, current positioning and optimizing method can not find target processing curve
And using blade as " waste product ", positioning and optimizing method proposed by the present invention allows design section line to be deformed in constraint,
Simultaneously ensure detection sectional plane wire shaped it is closest with design section line, improve search out can processing curve ability.
Additional aspect of the invention and advantage will be set forth in part in the description, and will partly become from the following description
Obtain substantially, or recognized by practice of the invention.
Brief description of the drawings
Of the invention above-mentioned and/or additional aspect and advantage will become from description of the accompanying drawings below to embodiment is combined
Substantially and be readily appreciated that, wherein:
Fig. 1 blade positioning and optimizing schematic flow sheets.
Fig. 2 blades design a model and type face section line design requirement schematic diagram.
Section line setting-out after Fig. 3 rotational deformations is curved surface schematic diagram.
The detection sectional plane line schematic diagram of Fig. 4 blades.
Fig. 5 B-spline curves concavity and convexities judge schematic diagram.
Fig. 6 is certain section Line contour degree tolerance range schematic diagram.
The relation schematic diagram of preform and the blade that designs a model before Fig. 7 positioning and optimizings.
The relation schematic diagram of preform and the blade that designs a model after Fig. 8 positioning and optimizings.
The relation schematic diagram of detection sectional plane line and profile tolerance tolerance after Fig. 9 positioning and optimizings.
Specific embodiment
Embodiments of the invention are described below in detail, the embodiment is exemplary, it is intended to for explaining the present invention, and
It is not considered as limiting the invention.
By taking the blade shown in Fig. 2 as an example, according to the flow in the content of the invention, 3-9 makes into one to the present invention referring to the drawings
The explanation of step, it is comprised the following steps that:
Step 1:5 design section lines along the parallel distribution in long-pending folded axle direction in blade design model(implication of unified representation design in subscript d texts), detection sectional plane line is obtained using following stepsThe parameter expression of (implication of unified representation detection in subscript m texts), l=1,2 ... 5.
Step 1.1:Design section line is represented using cubic B-spline parametric equation:
Wherein Ni,3T () is the basic function of B-spline Curve, ViIt is the coordinate of the control vertex of B-spline curves, m0Represent
The quantity of control vertex.
Step 1.2:Design section line is around reference axis x, the postrotational parameter expression of y, z:
Wherein Rl=Rx(αl)·Ry(θl)·Rz(βl) it is the l articles spin matrix of design section line, Rx(αl) it is the l articles
Design section line is around x-axis rotation alphalThe spin matrix at angle, Ry(θl) θ is rotated around y-axis for the l bars design section linelThe spin moment at angle
Battle array, Rz(βl) for the l bars design section line around z-axis rotation βlThe spin matrix at angle.
Step 1.3:Postrotational design section line base vector is obtained according to the spin matrix of design section line;Then adjust
CurvePosition of the control vertex on base vector direction obtain deformed sections line
The wherein the l articles shift in position amount of i-th control vertex of design section line is { Δ x 'i,l,Δy′i,l}。
Blade profile line be deformed by adjust B-spline curvesControl vertex along base vector
Position realize, section lineNot parallel xoy planes, therefore the base vector of plane after rotation need to be found.
The base vector of design section line is { (el,x,el,y) | l=1,2 ... 5 }, the base vector of plane is { (e ' after rotationl,x,e′l,y)=
(el,x·Rl,el,y·Rl) | l=1,2 ... 5 }, wherein (el,x,el,y) it is basal orientation of the l bars design section line on x, y directions
Amount, (e 'l,x,e′l,y) it is the l articles postrotational section line base vector planar.Therefore the expression formula of section line after deforming
For:
ΔViIt is i-th variation of control vertex, above formula is also referred to as:
WhereinCoordinate components of the section line on x, y, z direction after respectively deforming, { Δ
x′i,l,Δy′i,lIt is respectively section line after the l articles rotationI-th control vertex along base vector (e 'l,x,e′l,y)
The amount of movement in direction.
Step 1.4:By deformed sections lineSetting-out generates target processing curve s (u, v);Target is processed
The parametric equation of curved surface s (u, v) is
I.e.
Wherein sx(u,v),sy(u,v),sz(u, v) is respectively the target processing curve x at parameter (u, v) place, y, z value,
Vx,Vy,VzThe respectively coordinate V of control vertexi,jIn x, the component of y, z, Ni,3(u),Nj,3V () is respectively bi-cubic B spine bent
Basic function of the face on u and v directions, n1It is the number of control vertex on u direction, m1It is the number of control vertex on v directions, Vi,j
It is the control vertex coordinate of target processing curve.
Step 1.5:Solve detection sectional plane lineSubstantially solve target processing curve s (u, v) with
The intersection of the xoy planes of specified altitude assignment.The equation of the l articles plane is:Z=zl, wherein zlWith the l bars design section lineIn z
The height in direction is the same.For the l bars design section line, according to its z directions height zl, calculate s (u, v)=zlWhen
Intersection point parameterN0Represent the number of intersection point parameter.
Solved using Newton iteration method, it is contemplated that design section line anglec of rotation very little, use design section line's
ParameterDetection sectional plane line is calculated as Newton iteration methodParameterInitial value, and parameterBe
[0,1] it is discrete severalObtain.
WillTarget processing curve parametric equation is substituted into, the coordinate of discrete point is obtainedAccording to discrete point coordinates, the l bars detection sectional plane line is obtained using cubic spline interpolation methodParameter expression.
If the l bars detection sectional plane lineParametric equation be:
WhereinIt is the control vertex coordinate of detection sectional plane line, p (t) is discrete point coordinates, and parametric equation can be expressed as
Solve the control vertex V that above-mentioned matrix equation obtains curvei m, inspection is determined by the control vertex and knot vector of curve
Survey section lineParameter expression.
Step 2:According to the detection sectional plane line that step 1 is obtainedParameter expression, set up positioning it is excellent
The multiple objective function of change, optimization aim is that detection sectional plane line is closest with corresponding design section wire shaped.
Multiple objective function is:
WhereinRepresent that the single order of the l bars design section line cuts arrow,Represent the l articles single order of detection sectional plane line
Cut arrow;The l articles design section curvature of a curve is represented,The l articles detection sectional plane curvature of a curve is represented,Represent that the l articles sets
The second order for counting section line cuts arrow,Represent that the second order of the l bars detection sectional plane line cuts arrow.
Step 3:Foundation carries out constraint when multiple-objection optimization is solved:Including
Constraint 1, detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four are within upper profile tolerance tolerance range
Constraint function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four
In upper profile tolerance tolerance rangePoint;
Constraint 2, detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four are outside bottom profiled degree tolerance range
Constraint function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four
In bottom profiled degree tolerance rangePoint;Represent in the l articles profile tolerance tolerance range o'clock the l articles inspection
Survey corresponding closest approach on section line;Represent the list of corresponding closest approach on the l bars detection sectional plane line
The outer method arrow in position;I1=1,2 ... N1, i2=1,2 ... N2, i3=1,2 ... N3, i4=1,2 ... N4, N1,N2,N3,N4Respectively section line
The quantity of the point in leading edge, trailing edge, leaf basin and blade back region;
Constrain 3, to ensure that target processing curve has allowance, then must assure that each measurement point in target processing curve
Outside, so setting up constraint function F of the blade blank measurement point on the outside of target processing curve2:
F2=(pi-s(ui,vi))·ni>=0 i=1,2 ... N
Wherein piRepresent that on-machine measurement obtains the ith measurement point of the blade blank after near-net-shape technique, N is measurement point
Quantity;s(ui,vi) it is closest approach of the ith measurement point on target processing curve, niRepresent closest approach s (ui,vi) in target
The outer method arrow of unit on processing curve;
Without flex point constraint in the line segment of constraint 4, detection sectional plane line;
Detection sectional plane lineMathematic(al) representation be B-spline Curve.By the convex of B-spline curves
Bag property understands, when its controlling polygon is convex polygon, then its curve is convex curve;When controlling polygon is concave polygon
When, then curve is sag vertical curve.Ensure in curved section without flex point, then the controlling polygon convex-concave property of substantially section inner curve is
Consistent,
Set up the constraint f without flex point in convex detection sectional plane line segment1For
Set up the constraint f without flex point in recessed detection sectional plane line segment2For
viIt is i-th control vertex of curve, vi,x,vi,yIt is i-th x of control vertex, y-coordinate component, vkFor curve is recessed
The boundary control vertex of convexity.
And profile tolerance tolerance range is set up by procedure below in constraining 1 and constraint 2:
.The profile tolerance tolerance of blade profile line is considered design section lineAlong the offset line of normal distance, blade
Section line can be divided into leaf basin, blade back, leading edge, the subregion of trailing edge four according to action scope difference, and the profile tolerance of this four parts is public
Difference requires difference, and the mathematic(al) representation of the l articles design section line tolerance range is:
N (t)=t (t) × ez
WhereinIt is the l bars upper profile tolerance band of design section line,For upper profile variation amount (is united in subscript U texts
One meaning for representing, the meaning in subscript L texts under unified representation), n (t) is design section lineOutside the unit of parameter t
Method arrow,It is the l articles bottom profiled tolerance range of design section line,It is bottom profiled departure;T (t) is design section lineArrow, e are cut in the unit of parameter tzRepresent the unit vector in z directions;It is design section lineThe one of parameter t
Rank cuts arrow.
Step 4:With design section lineAround x, y, the z-axis anglec of rotation { αl,θl,βl, l=1,2 ... 5
And variation { the Δ x ' of control vertexi,l,Δ′yi,l, i=1,2 ... m0It is optimized variable, the multiple target set up to step 2
Function optimizes solution, obtains detection sectional plane line.
The design section line of solution is as shown in table 1 around x, y, z-axis rotation and along the result of the translational movement of base vector, excellent
Preform before change is with the relation for designing a model as shown in fig. 7, preform and the relation for designing a model after optimization
As shown in figure 8, the detection sectional plane line after optimization is as shown in Figure 9 with the relation of profile tolerance tolerance.See the present invention from optimum results
Positioning and optimizing method preform there is no allowance in the case of have found target processing curve, while meeting blade type
Face form tolerance requirement.
The optimized variable result of table 1
Although embodiments of the invention have been shown and described above, it is to be understood that above-described embodiment is example
Property, it is impossible to limitation of the present invention is interpreted as, one of ordinary skill in the art is not departing from principle of the invention and objective
In the case of above-described embodiment can be changed within the scope of the invention, change, replace and modification.
Claims (4)
1. a kind of towards the positioning and optimizing method based on blade parts deformation in adaptive machining, it is characterised in that:Including with
Lower step:
Step 1:N bars in blade design model along the parallel distribution in long-pending folded axle direction design section line
Detection sectional plane line is obtained using following stepsParameter expression, l=1,2 ... n;
Step 1.1:Design section line is represented using cubic B-spline parametric equation:
Wherein Ni,3T () is the basic function of B-spline Curve, ViIt is the coordinate of the control vertex of B-spline curves, m0Represent control
The quantity on summit;
Step 1.2:Design section line is around reference axis x, the postrotational parameter expression of y, z:
Wherein Rl=Rx(αl)·Ry(θl)·Rz(βl) it is the l articles spin matrix of design section line, Rx(αl) it is the l articles design
Section line is around x-axis rotation alphalThe spin matrix at angle, Ry(θl) θ is rotated around y-axis for the l bars design section linelThe spin matrix at angle, Rz
(βl) for the l bars design section line around z-axis rotation βlThe spin matrix at angle;
Step 1.3:Postrotational design section line base vector is obtained according to the spin matrix of design section line;Then curve is adjustedPosition of the control vertex on base vector direction obtain deformed sections lineIts
In the l articles shift in position amount of i-th control vertex of design section line be { Δ x 'i,l,Δy′i,l};
Step 1.4:By deformed sections lineSetting-out generates target processing curve s (u, v);Target processing curve
The parametric equation of s (u, v) is
Wherein Ni,3(u),Nj,3V () is respectively basic function of the Bicubic B-Spline Surfaces on u and v directions, n1To be controlled on u direction
The number on summit, m1It is the number of control vertex on v directions, Vi,jIt is the control vertex coordinate of target processing curve;
Step 1.5:For the l bars design section line, according to its z directions height zl, calculate s (u, v)=zlWhen intersection point ginseng
NumberN0Represent the number of intersection point parameter;WillTarget is substituted into add
Work surface parameter equation, obtains the coordinate of discrete pointAccording to discrete point coordinates, using interpolation
Method obtains the l bars detection sectional plane lineParameter expression;
Step 2:According to the detection sectional plane line that step 1 is obtainedParameter expression, set up positioning and optimizing
Multiple objective function, optimization aim is that detection sectional plane line is closest with corresponding design section wire shaped;
Step 3:Foundation carries out constraint when multiple-objection optimization is solved:Including
The constraint 1, constraint of detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four within upper profile tolerance tolerance range
Function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four upper
Profile tolerance tolerance rangePoint;
The constraint 2, constraint of detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four outside bottom profiled degree tolerance range
Function
WhereinRepresent the l articles detection sectional plane line leading edge, trailing edge, leaf basin and the subregion of blade back four under
Profile tolerance tolerance rangePoint;Represent in the l articles profile tolerance tolerance range o'clock the l articles detection cut
Corresponding closest approach on upper thread;Represent on the l bars detection sectional plane line outside the unit of corresponding closest approach
Method is sweared;I1=1,2 ... N1, i2=1,2 ... N2, i3=1,2 ... N3, i4=1,2 ... N4, N1,N2,N3,N4Respectively before section line
The quantity of the point in edge, trailing edge, leaf basin and blade back region;
Constraint 3, constraint function F of the blade blank measurement point on the outside of target processing curve2:
F2=(pi-s(ui,vi))·ni>=0 i=1,2 ... N
Wherein piRepresent that on-machine measurement obtains the ith measurement point of the blade blank after near-net-shape technique, N is the number of measurement point
Amount;s(ui,vi) it is closest approach of the ith measurement point on target processing curve, niRepresent closest approach s (ui,vi) in target processing
The outer method arrow of unit on curved surface;
Without flex point constraint in the line segment of constraint 4, detection sectional plane line;
Step 4:With design section lineAround x, y, the z-axis anglec of rotation { αl,θl,βl, l=1,2 ... n and
Variation { the Δ x ' of control vertexi,l,Δy′i,l, i=1,2 ... m0It is optimized variable, the multiple objective function set up to step 2
Solution is optimized, the design section line after being optimized, and obtain target processing curve s (u, v).
2. a kind of towards the positioning and optimizing method based on blade parts deformation in adaptive machining according to claim 1,
It is characterized in that:The profile tolerance tolerance range of the l bars design section line is
N (t)=t (t) × ez
WhereinIt is the l bars upper profile tolerance band of design section line,It is upper profile variation amount, n (t) is design section
LineIn the outer method arrow of the unit of parameter t,It is the l articles bottom profiled tolerance range of design section line,It is bottom profiled deviation
Amount;T (t) is design section lineArrow, e are cut in the unit of parameter tzRepresent the unit vector in z directions;It is design section
LineArrow is cut in the single order of parameter t.
3. according to claim 1 or claim 2 a kind of towards the positioning and optimizing side based on blade parts deformation in adaptive machining
Method, it is characterised in that:It is constrained to without flex point in the line segment of detection sectional plane line:The controlling polygon concavity of detection sectional plane line line segment
Matter is consistent.
4. a kind of towards the positioning and optimizing method based on blade parts deformation in adaptive machining according to claim 3,
It is characterized in that:Multiple objective function is:
WhereinRepresent that the single order of the l bars design section line cuts arrow,Represent that the single order of the l bars detection sectional plane line cuts arrow;The l articles design section curvature of a curve is represented,The l articles detection sectional plane curvature of a curve is represented,Represent that the l articles design cuts
The second order of upper thread cuts arrow,Represent that the second order of the l bars detection sectional plane line cuts arrow.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710057686.1A CN106897501B (en) | 2017-01-23 | 2017-01-23 | Towards the positioning and optimizing method based on blade parts deformation in adaptive machining |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710057686.1A CN106897501B (en) | 2017-01-23 | 2017-01-23 | Towards the positioning and optimizing method based on blade parts deformation in adaptive machining |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106897501A true CN106897501A (en) | 2017-06-27 |
CN106897501B CN106897501B (en) | 2019-07-05 |
Family
ID=59198663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710057686.1A Active CN106897501B (en) | 2017-01-23 | 2017-01-23 | Towards the positioning and optimizing method based on blade parts deformation in adaptive machining |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106897501B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536932A (en) * | 2018-03-26 | 2018-09-14 | 华中科技大学 | Based on the mutual aerial blade long-pending folded axle verticality computational methods turned round under constraints |
CN109332772A (en) * | 2018-11-22 | 2019-02-15 | 中国航发沈阳黎明航空发动机有限责任公司 | A kind of aero-engine stator blade front and rear edge adaptive machining method |
CN111259487A (en) * | 2019-12-31 | 2020-06-09 | 苏州千机智能技术有限公司 | Machining method for changing machining allowance of blisk |
CN111428334A (en) * | 2020-01-17 | 2020-07-17 | 西北工业大学 | Robot station planning method in laser radar measurement |
CN111498083A (en) * | 2020-04-15 | 2020-08-07 | 成都飞机工业(集团)有限责任公司 | Laminar flow wing aircraft aerodynamic outer edge tolerance control method |
CN113868802A (en) * | 2021-10-04 | 2021-12-31 | 西北工业大学 | Hollow blade margin optimization model building and solving method under variable wall thickness constraint |
CN114676526A (en) * | 2022-03-31 | 2022-06-28 | 中国科学院西安光学精密机械研究所 | Blade integral flexible deformation method based on mean camber line |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011153520A1 (en) * | 2010-06-04 | 2011-12-08 | The Gleason Works | Adaptive control of a machining process |
CN104750892A (en) * | 2013-12-31 | 2015-07-01 | 中国航空工业集团公司沈阳飞机设计研究所 | Three-dimensional modeling method for thickness-variable curved-surface part inner shape surface |
CN105261065A (en) * | 2015-11-02 | 2016-01-20 | 南京航空航天大学 | Front and rear edge profile redesign method for adaptive processing |
CN105739440A (en) * | 2016-04-29 | 2016-07-06 | 南京航空航天大学 | Adaptive machining method of wide-chord hollow fan blade |
CN106021782A (en) * | 2016-05-31 | 2016-10-12 | 西北工业大学 | Mean-line-based blade front and back edge fitting and section line smooth reconstruction method |
-
2017
- 2017-01-23 CN CN201710057686.1A patent/CN106897501B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011153520A1 (en) * | 2010-06-04 | 2011-12-08 | The Gleason Works | Adaptive control of a machining process |
CN104750892A (en) * | 2013-12-31 | 2015-07-01 | 中国航空工业集团公司沈阳飞机设计研究所 | Three-dimensional modeling method for thickness-variable curved-surface part inner shape surface |
CN105261065A (en) * | 2015-11-02 | 2016-01-20 | 南京航空航天大学 | Front and rear edge profile redesign method for adaptive processing |
CN105739440A (en) * | 2016-04-29 | 2016-07-06 | 南京航空航天大学 | Adaptive machining method of wide-chord hollow fan blade |
CN106021782A (en) * | 2016-05-31 | 2016-10-12 | 西北工业大学 | Mean-line-based blade front and back edge fitting and section line smooth reconstruction method |
Non-Patent Citations (1)
Title |
---|
万能等: "基于等几何法的叶片点解加工过程仿真方法研究", 《航空制造技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536932A (en) * | 2018-03-26 | 2018-09-14 | 华中科技大学 | Based on the mutual aerial blade long-pending folded axle verticality computational methods turned round under constraints |
CN108536932B (en) * | 2018-03-26 | 2020-07-28 | 华中科技大学 | Aviation blade stacking axis verticality calculation method based on mutual torsion constraint condition |
CN109332772A (en) * | 2018-11-22 | 2019-02-15 | 中国航发沈阳黎明航空发动机有限责任公司 | A kind of aero-engine stator blade front and rear edge adaptive machining method |
CN111259487A (en) * | 2019-12-31 | 2020-06-09 | 苏州千机智能技术有限公司 | Machining method for changing machining allowance of blisk |
CN111259487B (en) * | 2019-12-31 | 2023-10-13 | 苏州千机智能技术有限公司 | Processing method for changing machining allowance of blisk |
CN111428334A (en) * | 2020-01-17 | 2020-07-17 | 西北工业大学 | Robot station planning method in laser radar measurement |
CN111498083A (en) * | 2020-04-15 | 2020-08-07 | 成都飞机工业(集团)有限责任公司 | Laminar flow wing aircraft aerodynamic outer edge tolerance control method |
CN113868802A (en) * | 2021-10-04 | 2021-12-31 | 西北工业大学 | Hollow blade margin optimization model building and solving method under variable wall thickness constraint |
CN113868802B (en) * | 2021-10-04 | 2022-09-06 | 西北工业大学 | Hollow blade margin optimization model building and solving method under variable wall thickness constraint |
CN114676526A (en) * | 2022-03-31 | 2022-06-28 | 中国科学院西安光学精密机械研究所 | Blade integral flexible deformation method based on mean camber line |
CN114676526B (en) * | 2022-03-31 | 2023-01-06 | 中国科学院西安光学精密机械研究所 | Blade overall flexible deformation method based on mean camber line |
Also Published As
Publication number | Publication date |
---|---|
CN106897501B (en) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106897501A (en) | The positioning and optimizing method based on blade parts deformation towards in adaptive machining | |
CN103488832B (en) | A kind of geometry restorative procedure of complex curved surface parts damaged area | |
CN102637216B (en) | Method for generating numerical-control side milling machining tool path for complicated curved surfaces | |
CN107077126B (en) | The generation method and lathe of cutter path | |
CN106354927B (en) | A kind of construction method of precision forged blade front and rear edge adaptive machining Optimized model | |
CN104899345B (en) | Method for determining complex-curved shape workpiece laser shot forming technological parameter | |
CN105739440A (en) | Adaptive machining method of wide-chord hollow fan blade | |
CN110069041A (en) | A kind of Work piece processing method and system based on on-machine measurement | |
US20100095526A1 (en) | Method of repairing machined components such as turbomachine blades or blades of blisks | |
CN108262648A (en) | Axial ultrasonic vibration-assisted grinding workpiece surface appearance simulated prediction method | |
CN106994483B (en) | A kind of method of the accurate type face processing of Automobile Cover Drawing Die | |
CN109343468B (en) | Projection offset-based blade multi-axis track generation method | |
CN105069271A (en) | Bevel gear tooth surface machining error correction method | |
CN104697462A (en) | Method for extracting profiling feature parameters of aeronautical blade based on central axis | |
CN106125666A (en) | The Machining of Curved Surface cutter path planing method being constraint with cutting force fluctuation | |
CN108544181B (en) | Repair method for damaged blades of blisk | |
CN109213083A (en) | A kind of ruled surface processing path generation method, device and equipment | |
Koo et al. | Progress in aerodynamic shape optimization based on the reynolds-averaged navier-stokes equations | |
Dong et al. | An improved calculation method for cutting contact point and tool orientation analysis according to the CC points | |
CN103136426A (en) | Aviation blade circular arc leading-trailing edge process model generation method | |
Khaimovich et al. | Practice in the application of specialized technologic process patterns of 5-axis machining operation of blade rings of blisks | |
Neng et al. | A new localization theory of adaptive machining of near-net-shape blades | |
CN111474899A (en) | Triangular-based complex cavity high-speed numerical control milling spiral path generation method | |
CN113714731B (en) | Tolerance constraint-based blade surface vibration mark repairing method | |
CN109725595A (en) | Compensation method, processing method and the workpiece of the machining path of workpiece |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |