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 PDF

Info

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
Application number
CN201710057686.1A
Other languages
Chinese (zh)
Other versions
CN106897501B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710057686.1A priority Critical patent/CN106897501B/en
Publication of CN106897501A publication Critical patent/CN106897501A/en
Application granted granted Critical
Publication of CN106897501B publication Critical patent/CN106897501B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

The positioning and optimizing method based on blade parts deformation towards in adaptive machining
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：
C d ( t ) = { c l d ( t ) | c l d ( t ) = Σ i = 1 m 0 N i , 3 ( t ) V i } , l = 1 , 2 ... n
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：
C ^ d ( t ) = { c ^ l d ( t ) | c ^ l d ( t ) = Σ i = 1 m 0 N i , 3 ( t ) ( V i · R l ) } , l = 1 , 2 ... n
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
s ( u , v ) = Σ i = 1 n 1 Σ j = 1 m 1 V i , j N i , 3 ( u ) N j , 3 ( v )
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
f 1 , l U = ( q i 1 , l U - q i 1 , l m ) · n i 1 , l m ≥ 0 f 2 , l U = ( q i 2 , l U - q i 2 , l m ) · n i 1 , l m ≥ 0 f 3 , l U = ( q i 3 , l U - q i 3 , l m ) · n i 1 , l m ≥ 0 f 4 , l U = ( q i 4 , l U - q i 4 , l m ) · n i 1 , l m ≥ 0 , l = 1 , 2 ... n
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
f 1 , l L = ( q i 1 , l L - q i 1 , l m ) · n i 1 , l m ≤ 0 f 2 , l L = ( q i 2 , l L - q i 2 , l m ) · n i 2 , l m ≤ 0 f 3 , l L = ( q i 3 , l L - q i 3 , l m ) · n i 3 , l m ≤ 0 f 4 , l L = ( q i 4 , l L - q i 4 , l m ) · n i 4 , l m ≤ 0 , l = 1 , 2 ... n
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
c l U ( t ) = c l d ( t ) + μ l U · n ( t ) = Σ i = 1 m 0 N i , 3 ( t ) V i + μ l U · n ( t )
c l L ( t ) = c l d ( t ) - ϵ l L · n ( t ) = Σ i = 1 m 0 N i , 3 ( t ) V i - ϵ l L · n ( t )
N (t)=t (t) × ez
t ( t ) = c l ′ d ( t ) | | c l ′ d ( t ) | |
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：
min g 1 , l = ∫ 0 1 c l ′ d ( t ) × c l ′ m ( t ) d t , l = 1 , 2 ... n
min g 2 , l = ∫ 0 1 [ | k l d ( t ) - k l m ( t ) | ] d t = ∫ 0 1 | c l ′ d ( t ) × c l ′ ′ d ( t ) | | c l ′ d ( t ) | 3 - | c l ′ m ( t ) × c l ′ ′ m ( t ) | | c l ′ m ( t ) | 3 | d t , l = 1 , 2 ... n
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.
CN201710057686.1A 2017-01-23 2017-01-23 Towards the positioning and optimizing method based on blade parts deformation in adaptive machining Active CN106897501B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
万能等: "基于等几何法的叶片点解加工过程仿真方法研究", 《航空制造技术》 *

Cited By (11)

* Cited by examiner, † Cited by third party
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