CN106156409B - Two-parameter special-shaped bolt hole robust design method based on minimum contour change - Google Patents

Two-parameter special-shaped bolt hole robust design method based on minimum contour change Download PDF

Info

Publication number
CN106156409B
CN106156409B CN201610475937.3A CN201610475937A CN106156409B CN 106156409 B CN106156409 B CN 106156409B CN 201610475937 A CN201610475937 A CN 201610475937A CN 106156409 B CN106156409 B CN 106156409B
Authority
CN
China
Prior art keywords
radius
special
hole
shaped hole
dimensionless
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.)
Expired - Fee Related
Application number
CN201610475937.3A
Other languages
Chinese (zh)
Other versions
CN106156409A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201610475937.3A priority Critical patent/CN106156409B/en
Publication of CN106156409A publication Critical patent/CN106156409A/en
Application granted granted Critical
Publication of CN106156409B publication Critical patent/CN106156409B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention discloses a robust design method of two-parameter special-shaped bolt holes based on minimum contour change. In the optimization model, a special-shaped hole comprehensive optimization design model is established by reasonably weighting and summing each optimization target. Meanwhile, a dynamic weight coefficient adjusting method for balancing contribution of different design targets to the comprehensive optimization target is adopted, so that a balanced optimization design effect of the special-shaped hole under a multi-target condition is obtained, the dependence of optimization results of various design variables (namely, the size of the circular arc) of the special-shaped hole on the value range of the design variables is fundamentally eliminated, and the optimal special-shaped hole profile which meets the requirement condition of a strength design task and has the minimum profile size change can be quickly obtained through optimization. The design result has better consistency.

Description

Two-parameter special-shaped bolt hole robust design method based on minimum contour change
Technical Field
The invention belongs to the structural optimization design of turbine rotor parts of an aircraft engine, and particularly relates to a multi-objective optimization design model of the outline shape of a double-shaft symmetrical special-shaped bolt hole of a turbine disc.
Background
The high-pressure turbine disc of the aero-engine has a severe long-term working environment, and the bolt round holes on the front mounting edge of the disc often generate hole edge cracks due to stress concentration, so that the disc becomes an important reason for failure of the disc.
The patent "multi-arc shaped hole" (patent No. ZL 201310287309.9) proposes that an eight-arc shaped hole structure can be adopted. The following disadvantages still remain:
1. the optimal value of the objective function in the existing optimization model has an over-optimization phenomenon (namely, the stress reduction amplitude of the hole edge exceeds the design requirement).
2. The optimization design result excessively depends on the setting of the value range of the design variable, and in the actual design, continuous trial calculation is needed to determine the reasonable value range of the design variable, so that the design efficiency is low.
3. The robustness of the optimization design result is poor, and the convergence speed in the optimization process is low.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects in the prior art, the invention provides a two-parameter special-shaped bolt hole steady design method based on minimum outline change, which can avoid tedious and low-efficiency boundary trial calculation during optimization design, improve the rationality and stability of design and greatly improve the design efficiency.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the technical scheme that:
a robust design method for two-parameter special-shaped bolt holes based on minimum contour change comprises the following steps:
step 1, establishing a special-shaped hole multi-target optimization model according to a task given target value, a main arc radius and a transition arc radius and according to a stress reduction target, a special-shaped hole main arc size change target and a transition arc size change target.
And 2, determining the size constraint of the special-shaped hole multi-target optimization model in the step 1 according to the upper boundary of the radius of the main circular arc and the lower boundary of the radius of the transition circular arc.
And 3, introducing a weighting coefficient, and converting the special-shaped hole multi-target optimization model into a special-shaped hole single comprehensive target optimization model.
And 4, optimizing the single comprehensive target optimization model of the special-shaped hole according to the fixed contour change value and the allowable error to obtain each design variable of the special-shaped hole and complete the design of the special-shaped hole.
The method for establishing the special-shaped hole multi-objective optimization model in the step 1 comprises the following steps:
and 11, establishing a dimensionless stress reduction amplitude target function according to the task given target value, the radius of the main arc and the radius of the transition arc.
And step 12, establishing a contour control objective function of the dimensionless special-shaped hole with the main arc radius according to the radius of the original bolt through hole and the radius of the main arc.
And step 13, establishing a contour control target function of the dimensionless special-shaped hole with the transition arc radius according to the radius of the original bolt through hole and the radius of the transition arc.
And step 14, establishing the special-shaped hole multi-target optimization model according to the dimensionless stress reduction amplitude objective function established in the step 11, the main arc radius dimensionless special-shaped hole contour control objective function established in the step 12 and the transition arc radius dimensionless special-shaped hole contour control objective function established in the step 13.
The expression of the dimensionless stress reduction amplitude objective function in the step 11 is as follows:
Figure GDA0002195969590000021
f1max) For a dimensionless stress-reduced amplitude objective function, σmaxMaximum principal stress at the edge of the hole, σmax=σmax(R1,R2),R1Radius of main circular arc, R2To transition the radius of the arc, σoIndicating a task given target value.
The expression of the control objective function of the dimensionless special-shaped hole profile of the radius of the main circular arc in the step 12 is as follows:
Figure GDA0002195969590000022
wherein f is2(R1) Main circular arc radius dimensionless profile control objective function, R1The radius of the main circular arc is R, and the radius of the original bolt through hole is represented by R.
The expression of the control objective function of the dimensionless profile hole profile of the transition arc radius in the step 13 is as follows:
Figure GDA0002195969590000023
wherein f is3(R2) Control of the target function, R, for the contour of a dimensionless profiled hole of the radius of the transition arc2And (4) representing the radius of the transition circular arc, wherein R represents the radius of the original bolt through hole.
The expression of the special-shaped hole multi-objective optimization model in the step 14 is as follows:
minf(y1=f1max),y2=f2(R1),y3=f3(R2))。
the size constraint established in step 2 is:
Figure GDA0002195969590000031
wherein R is1maxIs the upper bound of the radius of the main circular arc, R2minIs the lower boundary of the transition arc radius, R1Radius of main circular arc, R2And (4) representing the radius of the transition circular arc, wherein R represents the radius of the original bolt through hole.
The abnormal hole single comprehensive target optimization model in the step 3 is as follows:
f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2)。
wherein k is1Weighting coefficients for the dimensional stress reduction amplitude objective function, f1max) Is a dimensionless stress reduction amplitude objective function, k2Controlling the weighting coefficient of the objective function for the dimensionless profile of the hole having the main circular arc radius, f2(R1) Control objective function, k, of main arc radius dimensionless profile hole contour3Controlling the weighting factor of the objective function for the dimensionless profile of the hole in the radius of the transition arc, f3(R2) And controlling an objective function for the profile of the dimensionless special-shaped hole with the radius of the transition arc.
The method for determining the weighting coefficient in the step 3 is as follows:
weighting coefficient k of target function of dimensional stress reduction amplitude1Weighting coefficient k of main circular arc radius dimensionless special-shaped hole outline control objective function2Weighting coefficient k of contour control objective function of dimensionless special-shaped hole with transition circular arc radius3With R1Radius of the main arc and sigmaoThe stress target values are dynamically varied, wherein: radius of main arc R1And transition arcRadius R2When the value range is changed, the influence on the hole edge stress is different.
R1maxWhen the change is carried out:
k2=0.008~0.012
k2and k3Satisfies the following conditions:
k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max
k1and k is2,k3Satisfies the following conditions:
Figure GDA0002195969590000032
wherein σ1rThe maximum main stress value of the edge of the through hole of the original circular bolt is obtained.
When the abnormal-shaped hole single comprehensive target optimization model is optimized in the step 4, only the stress of the optimized abnormal-shaped hole edge meets the set requirement, namely, the stress of the abnormal-shaped hole edge meets the set requirementmaxoAnd | ≦ ε, which is the allowable error. And compared with the original bolt circular through hole, the profile change of the special-shaped bolt hole is smaller than a given profile change value, and the design parameter or variable is an optimized value meeting the requirement.
Has the advantages that: compared with the prior art, the robust design method for the two-parameter special-shaped bolt hole based on the minimum profile change has the following beneficial effects:
1. the optimization target of the optimization problem of the two-variable special-shaped hole is perfected, so that the obtained optimization solution can reduce the hole edge stress according to the expected target.
2. The contour dimension of the special-shaped hole obtained by the optimized design can meet the requirement that the stress is reduced to a specified value, and the deviation from the original contour is not too large, so that the design is more stable.
3. The problem of excessive dependence of the original optimization result on the value range of the design variable is solved, and the size (R) of each circular arc of the special-shaped hole is not required to be adjusted repeatedly in the optimization design1,R2) Thereby improving the stability of the optimized design point, leading the final design scheme to be more reasonable and leading the design efficiency to be more reasonableThe improvement is obvious.
Drawings
FIG. 1 shows the weighting factor k1,k2,k3Dynamically determining a flow chart;
FIG. 2 is a profile of a shaped hole based on a multi-objective two-variable optimization model;
wherein: 1-two-variable special-shaped hole main circular arc R1(ii) a 2-two-variable special-shaped hole transition arc R2(ii) a 3-conventional bolt through hole (base circle).
Detailed Description
The present invention is further illustrated by the following description in conjunction with the accompanying drawings and the specific embodiments, it is to be understood that these examples are given solely for the purpose of illustration and are not intended as a definition of the limits of the invention, since various equivalent modifications will occur to those skilled in the art upon reading the present invention and fall within the limits of the appended claims.
A robust design method of two-parameter special-shaped bolt holes based on minimum profile change aims at the profile optimization problem of special-shaped bolt through holes (hereinafter referred to as special-shaped holes) on a high-pressure turbine disc of an aeroengine, and provides the robust design method. In the optimization model, a special-shaped hole comprehensive optimization design model is established by reasonably weighting and summing each optimization target. In the comprehensive optimization model, because the weight of each target changes along with the change of the value range of the design variable of the irregular hole, the patent inventor designs a dynamic weight coefficient adjusting method capable of keeping the contribution of different design targets to the comprehensive optimization target balanced, so that the weight coefficient (k) related to the stress control target1) Is associated with the value range of the profile design variable of the special-shaped hole, and when the latter changes, k1Can be automatically adjusted, thereby obtaining the balanced optimization design effect of the special-shaped holes under the condition of multiple targets, and fundamentally eliminating the optimization result pair of each design variable (namely the arc size) of the special-shaped holesThe dependency of the value range of the design variable can be quickly optimized to obtain the optimal special-shaped hole profile which meets the requirement condition of the strength design task and has the minimum profile size change. The invention enables the two-variable profiled hole contour design to be optimally designed according to the hole edge stress target value given by a task, meanwhile, the value range of each design variable size of the profiled hole does not need to be repeatedly tentatively adjusted, and the design result has better consistency.
The method specifically comprises the following steps:
step 1, establishing a special-shaped hole multi-target optimization model according to a task given target value, a main arc radius and a transition arc radius and according to a stress reduction target, a special-shaped hole main arc size change target and a transition arc size change target.
The method specifically comprises the following steps:
step 11, establishing a dimensionless stress reduction amplitude target function according to the task given target value, the radius of the main arc and the radius of the transition arc, f1max) For the dimensionless stress reduction amplitude objective function, it is necessary to ensure that the maximum principal stress of the hole edge of the special-shaped hole sufficiently approaches to a certain target value sigma set in advance (task setting)o. The expression is as follows:
f1max)=mind
d represents f1(R) maximum principal stress σ at hole edgemaxWith a predetermined stress target value sigmaoIs measured by the dimensional difference of (a).
Figure GDA0002195969590000051
f1max) For a dimensionless stress-reduced amplitude objective function, σmaxMaximum principal stress at the edge of the hole, σmax=σmax(R1,R2),R1Radius of main circular arc, R2To transition the radius of the arc, σoIndicating a task given target value. In specific implementation,. sigmaoThe stress target value is determined according to the design task target and the working condition
And step 12, establishing a main circular arc radius dimensionless special-shaped hole outline control target function according to the original bolt through hole radius and the main circular arc radius, wherein the target function enables the difference between the obtained special-shaped hole outline and the original outline to be as small as possible, and therefore a more stable design point is obtained. And taking the radius R of the original bolt through hole as a reference. The expression is as follows:
Figure GDA0002195969590000052
wherein f is2(R1) Main circular arc radius dimensionless profile control objective function, R1The radius of the main circular arc is R, and the radius of the original bolt through hole is represented by R.
And the objective function f2(R1),f3(R2) Ensuring that the target f is met1(R) on the premise of obtaining the profile hole contour dimension R (R)1,R2) As close as possible to the conventional design bolt through hole profile, resulting in a more robust design point.
f2(R1)=minδ1,f3(R2)=minδ2
Wherein, delta1,δ2Respectively represent main circular arcs R of the special-shaped holes1And transition arc radius R2The dimensionless difference from the base radius of the circular bolt hole is as follows.
And step 13, establishing a control objective function of the profile of the dimensionless special-shaped hole with the transition arc radius according to the radius of the original bolt through hole and the radius of the transition arc, wherein the objective function enables the difference between the profile of the special-shaped hole obtained by optimization and the original profile to be as small as possible, and therefore a more stable design point is obtained. And taking the radius R of the original bolt through hole as a reference. The expression is as follows:
Figure GDA0002195969590000061
wherein f is3(R2) Control of the target function, R, for the contour of a dimensionless profiled hole of the radius of the transition arc2And (4) representing the radius of the transition circular arc, wherein R represents the radius of the original bolt through hole.
In the optimization design, when the value range of the design variable is determined, R1max(the upper boundary of the radius of the main circular arc) can be 8-20R (R is the radius of the base circle) according to the actual diameter of the bolt hole; r2minThe lower boundary of the transition arc radius depends on the minimum radius under the actual processing condition, and the actual minimum diameter cutter radius can be adopted.
Step 14, establishing a special-shaped hole multi-target optimization model according to the dimensionless stress reduction amplitude target function established in the step 11, the main arc radius dimensionless special-shaped hole contour control target function established in the step 12 and the transition arc radius dimensionless special-shaped hole contour control target function established in the step 13, wherein the expression of the special-shaped hole multi-target optimization model is as follows:
minf(y1=f1max),y2=f2(R1),y3=f3(R2))
all design points R when solving a concrete optimization problemi(R1,R2) Corresponding to f of the objective function1,f2,f3A criterion space, denoted by F, can be formed as follows:
F={fi(f1,f2,f3)∈R3|f1=f1maxi),f2=f2(R1),f3=f3(R2)}
fi=(f1,f2,f3) Is a vector consisting of three objective function values. The optimization process is to all vectors f in the criterion spaceiThe evaluation process of (1).
σmax=σmax(R1,R2),R1,R2Respectively represent a main circular arc and a transition circular arc, which represent a design point R; f. of1max) Is equal to the design point R (R)1,R2) A corresponding stress reduction amplitude objective function; f. of2(R1) Is a main circular arc R1Varying a control objective function; f. of3(R2) Is a transition arc R2The variation controls the objective function. R is the radius of the base circle (i.e., the conventional design bolt through hole), R1maxIs the upper bound of the radius of the main circular arc, R2minThe lower bound of the radius of the transition arc.
The special-shaped hole multi-target optimization model meets the following requirements:
(1) reducing the stress of the hole edge to the target given by the strength design task;
(2) the change of the overall size of the hole is minimal;
(3) the dependence of the optimized profile hole contour dimension on the value range of the design variable is effectively reduced, and the design robustness is improved.
And 2, determining the size constraint of the special-shaped hole multi-target optimization model in the step 1 according to the upper boundary of the radius of the main circular arc and the lower boundary of the radius of the transition circular arc.
Its dimensional constraints are:
Figure GDA0002195969590000071
wherein R is1maxIs the upper bound of the radius of the main circular arc, and can be arbitrarily valued in a reasonable range, R2minIs the lower boundary of the transition arc radius, R1Radius of main circular arc, R2For the transition arc radius, R denotes the base circle, i.e. the radius of the (original) bolt through hole of the conventional design.
Step 3, introducing a weighting coefficient, and converting the abnormal-shaped hole multi-target optimization model into an abnormal-shaped hole single comprehensive target optimization model, namely a target f1max)、f2(R1) And f3(R2) And performing weighted summation processing to simplify the multi-objective optimization problem into a single comprehensive objective.
The optimization model of the single comprehensive target of the special-shaped hole is as follows:
f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2)。
wherein k is1Weighting coefficients for the dimensional stress reduction amplitude objective function, f1max) Is a dimensionless stress reduction amplitude objective function, k2Controlling the objective function for the profile of a dimensionless profiled hole having a radius of the main arcWeighting coefficient, f2(R1) Control objective function, k, of main arc radius dimensionless profile hole contour3Controlling the weighting factor of the objective function for the dimensionless profile of the hole in the radius of the transition arc, f3(R2) And controlling an objective function for the profile of the dimensionless special-shaped hole with the radius of the transition arc.
The multi-objective optimization design model formed by the three new objective functions can simultaneously control the maximum main stress reduction amplitude of the hole edge and the change of the outline size of the special-shaped hole, so that the hole edge stress reduction amplitude meets the set target, the hole pattern design is more stable, and the hole pattern parameters are more stable.
By a weighting coefficient k1,k2,k3The contributions of the three single optimization objectives to the composite optimization objective may be adjusted.
As shown in fig. 1, the determination method of the weighting factor is as follows:
weighting coefficient k of target function of dimensional stress reduction amplitude1Weighting coefficient k of main circular arc radius dimensionless special-shaped hole outline control objective function2Weighting coefficient k of contour control objective function of dimensionless special-shaped hole with transition circular arc radius3With R1Radius of the main arc and sigmaoThe stress target values are dynamically varied, wherein: radius of main arc R1And transition arc radius R2When the value range is changed, the influence on the hole edge stress is different.
R1maxWhen the change is carried out:
k2=0.008~0.012
k2and k3Satisfies the following conditions:
k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max
k1and k is2,k3Satisfies the following conditions:
Figure GDA0002195969590000081
wherein σ1rIs in the shape of a circleThe maximum main stress value of the hole edge of the bolt through hole.
Generated weighting coefficient k1,k2,k3The change of each design parameter can be more balanced, and a steady design point is obtained.
And 4, optimizing the single comprehensive target optimization model of the special-shaped hole according to the fixed contour change value and the allowable error to obtain each design variable of the special-shaped hole and complete the design of the special-shaped hole.
When the optimization model of the single comprehensive target of the special-shaped hole is optimized, the stress of the special-shaped hole edge after optimization reaches the set requirement, namely, the stress of the | sigma & lt/EN & gtmaxoAnd | ≦ ε, which is the allowable error. Compared with the original bolt circular through hole, when the profile change of the special-shaped bolt hole is smaller than a given profile change value, the design parameter or variable is an optimized value meeting the requirement, and three single targets are required to contribute to the comprehensive design result in a given multi-target optimization model in a balanced manner. Due to the weighting coefficient k1,k2,k3In relation to the value range of the design variable, a weighting factor k is given for this purpose1,k2,k3Dynamically generating a solution such that k1,k2,k3The value can follow R1maxAnd adjusting the value according to different value ranges to meet the balance among three single optimization targets, finally obtaining a design point which meets the strength design requirement and is more stable and stable, and ensuring that the three have proper influence factors, thereby ensuring the robustness of the final optimal solution.
The multi-objective optimization design method provided by the invention can obtain a hole pattern design scheme with small hole outline change and more stability on the basis of meeting a stress design target (task value); in addition, in the optimization design, when the allowable variation range of the arc size of the special-shaped hole is obviously changed, a stable optimization solution can still be kept, so that the final optimization design result is free from the dependence on the set range value of the arc design variable, and the design efficiency is improved. The structure of the two-parameter special-shaped hole obtained under the multi-objective optimization model is shown in the attached figure 2.
The invention provides a novel multi-objective optimization design model, which can be used for designing holes according to task specificationsThe edge stress reduction amplitude provides a steady and stable design scheme. The multi-objective optimization model consists of (a) a stress reduction target and (b) a main circular arc size R of a special-shaped hole1Varying the target and (c) the transition arc size R2A change target composition; in the multi-objective optimization model, the contributions of three single objectives (namely (a), (b) and (c)) to the comprehensive objective meet the equilibrium condition; the provided new model and algorithm can avoid complex and inefficient boundary trial calculation during optimization design, improve the reasonability and stability of design and greatly improve the design efficiency.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.

Claims (5)

1. A robust design method for two-parameter special-shaped bolt holes based on minimum profile variation is characterized by comprising the following steps:
step 1, establishing a special-shaped hole multi-target optimization model according to a task given target value, a main arc radius and a transition arc radius and according to a stress reduction target, a special-shaped hole main arc size change target and a transition arc size change target;
the method for establishing the special-shaped hole multi-target optimization model comprises the following steps:
step 11, establishing a dimensionless stress reduction amplitude target function according to a task given target value, a main arc radius and a transition arc radius;
the expression of the dimensionless stress reduction amplitude objective function is as follows:
Figure FDA0002195969580000011
f1max) For a dimensionless stress-reduced amplitude objective function, σmaxMaximum principal stress at the edge of the hole, σmax=σmax(R1,R2),R1Radius of main circular arc, R2Is a transition arc halfDiameter, σoRepresenting a task given target value;
step 12, establishing a main circular arc radius dimensionless special-shaped hole outline control objective function according to the original bolt through hole radius and the main circular arc radius;
the expression of the control objective function of the dimensionless special-shaped hole profile of the radius of the main circular arc in the step 12 is as follows:
Figure FDA0002195969580000012
wherein f is2(R1) Main circular arc radius dimensionless profile control objective function, R1The radius of the main circular arc is R, and the radius of the original bolt through hole is represented by R;
step 13, establishing a transition arc radius dimensionless special-shaped hole profile control target function according to the original bolt through hole radius and the transition arc radius;
the expression of the control objective function of the dimensionless special-shaped hole profile of the radius of the transition circular arc is as follows:
Figure FDA0002195969580000013
wherein f is3(R2) Control of the target function, R, for the contour of a dimensionless profiled hole of the radius of the transition arc2The radius of a transition arc is shown, and R represents the radius of an original bolt through hole;
step 14, establishing a special-shaped hole multi-target optimization model according to the dimensionless stress reduction amplitude objective function established in the step 11, the main arc radius dimensionless special-shaped hole contour control objective function established in the step 12 and the transition arc radius dimensionless special-shaped hole contour control objective function established in the step 13;
the expression of the special-shaped hole multi-objective optimization model is as follows:
minf(y1=f1max),y2=f2(R1),y3=f3(R2));
step 2, determining the size constraint of the special-shaped hole multi-target optimization model in the step 1 according to the upper boundary of the radius of the main arc and the lower boundary of the radius of the transition arc;
step 3, introducing a weighting coefficient, and converting the special-shaped hole multi-target optimization model into a special-shaped hole single comprehensive target optimization model;
and 4, optimizing the single comprehensive target optimization model of the special-shaped hole according to the fixed contour change value and the allowable error to obtain each design variable of the special-shaped hole and complete the design of the special-shaped hole.
2. The two-parameter deformed bolt hole robust design method based on minimum profile variation according to claim 1, characterized in that: the size constraint established in step 2 is:
Figure FDA0002195969580000021
wherein R is1maxIs the upper bound of the radius of the main circular arc, R2minIs the lower boundary of the transition arc radius, R1Radius of main circular arc, R2And (4) representing the radius of the transition circular arc, wherein R represents the radius of the original bolt through hole.
3. The two-parameter deformed bolt hole robust design method based on the minimum profile variation as claimed in claim 2, characterized in that: the abnormal hole single comprehensive target optimization model in the step 3 is as follows:
f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2);
wherein k is1Weighting coefficients for the dimensional stress reduction amplitude objective function, f1max) Is a dimensionless stress reduction amplitude objective function, k2Controlling the weighting coefficient of the objective function for the dimensionless profile of the hole having the main circular arc radius, f2(R1) Control objective function, k, of main arc radius dimensionless profile hole contour3Controlling the weighting factor of the objective function for the dimensionless profile of the hole in the radius of the transition arc, f3(R2) Is a transition arcThe radius dimensionless profile of the hole controls the objective function.
4. The robust design method for two-parameter special-shaped bolt holes based on minimum profile variation according to claim 3, characterized in that: the method for determining the weighting coefficient in the step 3 is as follows:
weighting coefficient k of target function of dimensional stress reduction amplitude1Weighting coefficient k of main circular arc radius dimensionless special-shaped hole outline control objective function2Weighting coefficient k of contour control objective function of dimensionless special-shaped hole with transition circular arc radius3With R1Radius of the main arc and sigmaoThe stress target values are dynamically varied, wherein: radius of main arc R1And transition arc radius R2When the value range is changed, the influence on the hole edge stress is different;
R1maxwhen the change is carried out:
k2=0.008~0.012
k2and k3Satisfies the following conditions:
k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max
k1and k is2,k3Satisfies the following conditions:
Figure FDA0002195969580000031
wherein σ1rThe maximum main stress value of the edge of the through hole of the original circular bolt is obtained.
5. The robust design method for two-parameter special-shaped bolt holes based on minimum profile variation according to claim 4, characterized in that: when the abnormal-shaped hole single comprehensive target optimization model is optimized in the step 4, only the stress of the optimized abnormal-shaped hole edge meets the set requirement, namely, the stress of the abnormal-shaped hole edge meets the set requirementmaxoEpsilon is less than or equal to | and epsilon is an allowable error; and compared with the original bolt circular through hole, the profile change of the special-shaped bolt hole is smaller than the given profile change valueThe design parameter or variable is the optimum value that meets the requirements.
CN201610475937.3A 2016-06-24 2016-06-24 Two-parameter special-shaped bolt hole robust design method based on minimum contour change Expired - Fee Related CN106156409B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610475937.3A CN106156409B (en) 2016-06-24 2016-06-24 Two-parameter special-shaped bolt hole robust design method based on minimum contour change

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610475937.3A CN106156409B (en) 2016-06-24 2016-06-24 Two-parameter special-shaped bolt hole robust design method based on minimum contour change

Publications (2)

Publication Number Publication Date
CN106156409A CN106156409A (en) 2016-11-23
CN106156409B true CN106156409B (en) 2020-01-14

Family

ID=57349242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610475937.3A Expired - Fee Related CN106156409B (en) 2016-06-24 2016-06-24 Two-parameter special-shaped bolt hole robust design method based on minimum contour change

Country Status (1)

Country Link
CN (1) CN106156409B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108664704B (en) * 2018-04-12 2021-05-07 南京航空航天大学 Five-parameter multi-arc special-shaped hole balance optimization method
CN109344494B (en) * 2018-09-28 2021-09-03 中国航发湖南动力机械研究所 Bolt hole of bolt connection type gas turbine baffle and design method thereof

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103334798A (en) * 2013-07-10 2013-10-02 南京航空航天大学 Multi-arc specially-shaped hole
CN105046000A (en) * 2015-07-20 2015-11-11 南京航空航天大学 Engineering discrete variable based two-variable multi-circular arc special-shaped hole balance design method and model

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103334798A (en) * 2013-07-10 2013-10-02 南京航空航天大学 Multi-arc specially-shaped hole
CN105046000A (en) * 2015-07-20 2015-11-11 南京航空航天大学 Engineering discrete variable based two-variable multi-circular arc special-shaped hole balance design method and model

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Evolutionary multi-objective optimization: a historical view of the field;C.A. Coello Coello;《IEEE Computational Intelligence Magazine》;20061231;第1卷(第1期);第28-36页 *
涡轮连接单元异型孔应力场分析及结构优化;张超;《中国优秀硕士学位论文全文数据库工程科技II辑》;20150115(第1期);第40-65页 *

Also Published As

Publication number Publication date
CN106156409A (en) 2016-11-23

Similar Documents

Publication Publication Date Title
CN109409013B (en) Intelligent optimization design method for wind wheel of low-wind-speed wind turbine generator
CN113094840B (en) Axial flow type anti-cavitation impeller design method based on airfoil parameter control
CN106156409B (en) Two-parameter special-shaped bolt hole robust design method based on minimum contour change
CN109145506B (en) Optimal design method for outer side wing profile of wind turbine with high aerodynamic performance and low noise level
CN112329170B (en) Blade airfoil section optimization design method for tidal current energy water turbine
CN105404743B (en) The wind mill airfoil design method that B-spline is combined with the smooth continuity of curvature
CN111027258B (en) Intelligent prediction method for power generation load and heat supply quantity of supercritical unit
CN103488847A (en) Aerodynamic shape optimization method based on neural network integration
CN105700357B (en) Method of Boiler Combustion Control System based on multivariable PID-PFC
CN112502883A (en) Water turbine speed regulator parameter setting method and device considering frequency stability characteristic
Liu et al. Optimization model for rotor blades of horizontal axis wind turbines
CN109783912B (en) Optimized design method for lift type vertical axis wind turbine blade airfoil profile
CN110361974B (en) Water turbine speed regulating system optimization method based on BP-FOA hybrid algorithm
Vasilopoulos et al. CAD-based aerodynamic optimization of a compressor stator using conventional and adjoint-driven approaches
CN110571828A (en) primary frequency modulation optimization method based on target flow optimization
CN110970936B (en) Method for calculating primary frequency modulation performance of deep peak shaving unit
CN106971023B (en) Wheel disc special-shaped hole structure design method based on hyperelliptic curve
CN106294908B (en) Sound serves as a contrast design method
CN112780484B (en) Torsion angle distribution design method and device of wind power blade and electronic equipment
CN109441287A (en) A kind of blower fan tower barrel doorframe Parameters design
CN112906166B (en) Wind turbine blade optimization design method considering aerodynamic efficiency and aerodynamic load
CN111259487A (en) Machining method for changing machining allowance of blisk
CN112943479A (en) Aero-engine acceleration process optimal control method based on improved simplex method
CN113806878B (en) Centrifugal fan blade profile optimization method based on steepest descent method
CN113236583A (en) Design method for arc tenon blade root extending structure

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200114