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
special
arc radius
shaped hole
hole
objective function
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

本发明公开了一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,本发明在优化设计中,采用(a)减小孔边应力至任务规定值和(b)控制异型孔轮廓变化幅度(即避免异型孔轮廓与原圆形轮廓差异过大)等优化目标。在优化模型中,通过对每个优化目标进行合理的加权求和处理,建立了异型孔综合优化设计模型。同时采用不同设计目标对综合优化目标贡献均衡的权重系数动态调整方法,获得了多目标条件下的异型孔均衡优化设计效果,从根本上消除了异型孔各设计变量(即圆弧尺寸)优化结果对设计变量取值范围的依赖性,可快速优化得到满足强度设计任务要求条件且的轮廓尺寸变化最小的最佳异型孔轮廓。设计结果具有较好的一致性。

Figure 201610475937

The invention discloses a robust design method for two-parameter special-shaped bolt holes based on the minimum contour change. In the optimized design, the invention adopts (a) reducing the stress on the edge of the hole to a specified value and (b) controlling the change range of the special-shaped hole contour (that is, to avoid the excessive difference between the contour of the special-shaped hole and the original circular contour) and other optimization goals. In the optimization model, a comprehensive optimization design model for special-shaped holes is established by performing a reasonable weighted summation process for each optimization objective. At the same time, the dynamic adjustment method of the weight coefficient in which the contribution of different design objectives to the comprehensive optimization objective is balanced is adopted, and the balanced optimization design effect of the special-shaped hole under the condition of multiple objectives is obtained, which fundamentally eliminates the optimization results of each design variable (ie arc size) of the special-shaped hole. The dependence on the value range of the design variables can be quickly optimized to obtain the best special-shaped hole contour that meets the requirements of the strength design task and has the smallest contour size change. The design results have good consistency.

Figure 201610475937

Description

一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法A Robust Design Method for Two-Parameter Special-shaped Bolt Holes Based on Minimum Profile Variation

技术领域technical field

本发明属于航空发动机涡轮转子部件的结构优化设计,尤其涉及涡轮盘双轴对称异型螺栓孔的轮廓形状的多目标优化设计模型。The invention belongs to the structural optimization design of aero-engine turbine rotor components, in particular to a multi-objective optimization design model of the contour shape of a turbine disk biaxially symmetrical special-shaped bolt hole.

背景技术Background technique

航空发动机高压涡轮盘长期工作环境恶劣,其前安装边上螺栓圆孔常因应力集中现象,产生孔边裂纹,成为轮盘失效的一个重要原因。The long-term working environment of aero-engine high-pressure turbine disks is harsh, and the bolt holes on the front mounting side often cause cracks on the side of the holes due to stress concentration, which is an important reason for the failure of the disk.

专利“多圆弧异型孔”(专利号:ZL2013 1 0287309.9)提出了可采用一种八圆弧异型孔结构。但仍存在如下不足:The patent "multi-arc special-shaped hole" (patent number: ZL2013 1 0287309.9) proposes an eight-arc special-shaped hole structure. But there are still the following shortcomings:

1.已有的优化模型中的目标函数的最优值存在过度优化现象(即孔边应力下降幅超过设计要求)。1. The optimal value of the objective function in the existing optimization model has an over-optimization phenomenon (that is, the drop of the hole edge stress exceeds the design requirement).

2.优化设计结果过分依赖于设计变量取值范围的设定,实际设计中,需不断试算以确定合理的设计变量取值范围,设计效率低。2. The results of the optimal design are overly dependent on the setting of the value range of the design variables. In the actual design, continuous trial calculations are required to determine the reasonable range of the design variables, and the design efficiency is low.

3.优化设计结果鲁棒性差,优化过程收敛速度较慢。3. The robustness of the optimization design results is poor, and the convergence speed of the optimization process is slow.

发明内容SUMMARY OF THE INVENTION

发明目的:为了克服现有技术中存在的不足,本发明提供一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,该方法可避免优化设计时繁琐、低效的边界试算,提高设计的合理性及稳定性,同时大大提高设计效率。Purpose of the invention: In order to overcome the deficiencies in the prior art, the present invention provides a robust design method for two-parameter special-shaped bolt holes based on minimum contour changes, which can avoid cumbersome and inefficient boundary trial calculations during optimal design, and improve design. the rationality and stability, and greatly improve the design efficiency.

技术方案:为实现上述目的,本发明采用的技术方案为:Technical scheme: In order to realize the above-mentioned purpose, the technical scheme adopted in the present invention is:

一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,包括以下步骤:A robust design method for two-parameter special-shaped bolt holes based on minimal contour variation, including the following steps:

步骤1.根据任务给定目标值、主圆弧半径、过渡圆弧半径、以应力下降目标、异型孔主圆弧尺寸变化目标以及过渡圆弧尺寸变化目标建立异型孔多目标优化模型。Step 1. Establish a multi-objective optimization model for special-shaped holes according to the given target value of the task, the radius of the main arc, the radius of the transition arc, the target of stress reduction, the target of the size of the main arc of the special-shaped hole and the change of the size of the transition arc.

步骤2.根据主圆弧半径上界、过渡圆弧半径下界确定步骤1中的异型孔多目标优化模型的尺寸约束。Step 2. Determine the size constraints of the multi-objective optimization model of the special-shaped hole in step 1 according to the upper bound of the main arc radius and the lower bound of the transition arc radius.

步骤3.引用加权系数,将异型孔多目标优化模型转化为异型孔单一综合目标优化模型。Step 3. Referring to the weighting coefficient, the multi-objective optimization model of the special-shaped hole is converted into a single comprehensive objective optimization model of the special-shaped hole.

步骤4.根据定轮廓变化值和允许的误差对异型孔单一综合目标优化模型进行优化,获得异型孔各设计变量,完成异型孔的设计。Step 4. Optimize the single comprehensive target optimization model of the special-shaped hole according to the change value of the fixed contour and the allowable error, obtain each design variable of the special-shaped hole, and complete the design of the special-shaped hole.

所述步骤1中异型孔多目标优化模型建立方法包括以下步骤:The method for establishing a multi-objective optimization model for the special-shaped hole in the step 1 includes the following steps:

步骤11,根据任务给定目标值、主圆弧半径、过渡圆弧半径建立无量纲应力下降幅目标函数。Step 11: Establish a dimensionless stress reduction objective function according to the given target value of the task, the radius of the main arc, and the radius of the transition arc.

步骤12,根据原螺栓通孔半径、主圆弧半径建立主圆弧半径无量纲异型孔轮廓控制目标函数。Step 12: Establish the main arc radius dimensionless special-shaped hole contour control objective function according to the original bolt through hole radius and the main arc radius.

步骤13,根据原螺栓通孔半径、过渡圆弧半径建立过渡圆弧半径无量纲异型孔轮廓控制目标函数。Step 13 , according to the original bolt through hole radius and the transition arc radius, establish the transition arc radius dimensionless special-shaped hole contour control objective function.

步骤14,根据步骤11建立的无量纲应力下降幅目标函数、步骤12建立的主圆弧半径无量纲异型孔轮廓控制目标函数、步骤13建立的过渡圆弧半径无量纲异型孔轮廓控制目标函数建立异型孔多目标优化模型。Step 14, establish according to the objective function of the dimensionless stress reduction amplitude established in step 11, the main arc radius dimensionless special-shaped hole contour control objective function established in step 12, and the transition arc radius dimensionless special-shaped hole contour control objective function established in step 13. A multi-objective optimization model for shaped holes.

所述步骤11中无量纲应力下降幅目标函数的表达式如下:In the step 11, the expression of the objective function of the dimensionless stress reduction amplitude is as follows:

Figure GDA0002195969590000021
Figure GDA0002195969590000021

f1max)为无量纲应力下降幅目标函数,σmax为孔边最大主应力,σmax=σmax(R1,R2),R1为主圆弧半径,R2为过渡圆弧半径,σo表示任务给定目标值。f 1max ) is the objective function of the dimensionless stress reduction range, σ max is the maximum principal stress at the edge of the hole, σ maxmax (R 1 , R 2 ), R 1 is the main arc radius, and R 2 is the transition circle Arc radius, σ o represents the task given target value.

所述步骤12中主圆弧半径无量纲异型孔轮廓控制目标函数的表达式如下:In described step 12, the expression of the main arc radius dimensionless special-shaped hole contour control objective function is as follows:

Figure GDA0002195969590000022
Figure GDA0002195969590000022

其中,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,R1为主圆弧半径,R*表示原螺栓通孔半径。Among them, f 2 (R 1 ) is the main arc radius of the dimensionless special-shaped hole contour control objective function, R 1 is the main arc radius, and R* represents the original bolt through hole radius.

所述步骤13中过渡圆弧半径无量纲异型孔轮廓控制目标函数的表达式如下:In described step 13, the expression of the control objective function of the non-dimensionally shaped hole contour of the transition arc radius is as follows:

Figure GDA0002195969590000023
Figure GDA0002195969590000023

其中,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数,R2为过渡圆弧半径,R*表示原螺栓通孔半径。Among them, f 3 (R 2 ) is the transition arc radius dimensionless hole contour control objective function, R 2 is the transition arc radius, and R* represents the original bolt through hole radius.

所述步骤14中异型孔多目标优化模型的表达式如下:The expression of the multi-objective optimization model for special-shaped holes in the step 14 is as follows:

minf(y1=f1max),y2=f2(R1),y3=f3(R2))。minf(y 1 =f 1max ), y 2 =f 2 (R 1 ), y 3 =f 3 (R 2 )).

所述步骤2中建立的尺寸约束为:The size constraints established in step 2 are:

Figure GDA0002195969590000031
Figure GDA0002195969590000031

其中,R1max为主圆弧半径上界,R2min为过渡圆弧半径下界,R1为主圆弧半径,R2为过渡圆弧半径,R*表示原螺栓通孔半径。Among them, R 1max is the upper bound of the main arc radius, R 2min is the lower bound of the transition arc radius, R 1 is the main arc radius, R 2 is the transition arc radius, and R* represents the original bolt through hole radius.

所述步骤3中的异型孔单一综合目标优化模型为:The single comprehensive objective optimization model of the special-shaped hole in the step 3 is:

f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2)。f(f 1 , f 2 , f 3 )=k 1 ·f 1max )+k 2 ·f 2 (R 1 )+k 3 ·f 3 (R 2 ).

其中,k1为量纲应力下降幅目标函数的加权系数,f1max)为无量纲应力下降幅目标函数,k2为主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,k3为过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数。Among them, k 1 is the weighting coefficient of the objective function of dimensionless stress reduction amplitude, f 1max ) is the objective function of dimensionless stress reduction amplitude, k 2 is the weighting coefficient of the main arc radius dimensionless special-shaped hole contour control objective function, f 2 (R 1 ) the main arc radius dimensionless special-shaped hole contour control objective function, k 3 is the weighting coefficient of the transition arc radius dimensionless special-shaped hole contour control objective function, f 3 (R 2 ) is the transition arc radius infinite The contour control objective function of the special-shaped hole.

所述步骤3中加权系数的确定方法如下:The method for determining the weighting coefficient in the step 3 is as follows:

量纲应力下降幅目标函数的加权系数k1、主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k2、过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k3随R1主圆弧半径的取值范围以及σo应力目标值的不同而动态变化,其中:主圆弧半径R1和过渡圆弧半径R2取值范围变化时,对孔边应力的影响不同。The weighting coefficient k 1 of the objective function of the dimensional stress reduction amplitude, the weighting coefficient k 2 of the main arc radius dimensionless special-shaped hole contour control objective function, the weighting coefficient k 3 of the transition arc radius dimensionless special-shaped hole contour control objective function varies with R 1 The value range of the main arc radius and the target value of σ o stress vary dynamically, among which: when the value range of the main arc radius R 1 and the transition arc radius R 2 changes, the impact on the hole edge stress is different.

R1max变化时:When R 1max changes:

k2=0.008~0.012k 2 =0.008~0.012

k2和k3之间满足:Between k 2 and k 3 :

k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max k 3 =g 1 (k 2 ,R 1max )=(1.25 to 2.1)·k 2 ·R 1max

k1与k2,k3之间满足:Between k 1 and k 2 , k 3 satisfies:

Figure GDA0002195969590000032
Figure GDA0002195969590000032

其中,σ1r为原圆形螺栓通孔孔边最大主应力值。Among them, σ 1r is the maximum principal stress value of the original circular bolt through hole edge.

所述步骤4中对异型孔单一综合目标优化模型进行优化时,只有在优化后的异型孔孔边应力达到既定要求,即|σmaxo|≤ε,ε为允许误差。且异型螺栓孔与原螺栓圆形通孔相比,轮廓变化小于给定轮廓变化值时,设计参数或变量才是满足要求的优化值。When optimizing the single comprehensive target optimization model of the special-shaped hole in the step 4, only when the optimized special-shaped hole hole edge stress meets the predetermined requirements, that is, |σ maxo |≤ε, ε is the allowable error. And compared with the original bolt circular through hole, when the contour change of the special-shaped bolt hole is less than the given contour change value, the design parameters or variables are the optimized values that meet the requirements.

有益效果:本发明提供的一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,相比现有技术,具有以下有益效果:Beneficial effects: Compared with the prior art, a robust design method for two-parameter special-shaped bolt holes based on minimum profile changes provided by the present invention has the following beneficial effects:

1.完善了两变量异型孔优化问题的优化目标,使得到的优化解可按预期目标降低孔边应力。1. The optimization objective of the two-variable special-shaped hole optimization problem is improved, so that the obtained optimal solution can reduce the hole edge stress according to the expected objective.

2.优化设计获得的异型孔轮廓尺寸既可满足应力下降到规定值,又不至偏离原有轮廓过大,因而设计更加稳健。2. The contour size of the special-shaped hole obtained by the optimized design can not only meet the stress reduction to the specified value, but also not deviate too much from the original contour, so the design is more robust.

3.克服了原优化结果对设计变量取值范围的过分依赖的问题,无需再优化设计中反复调整异型孔各圆弧尺寸(R1,R2)的设计范围,从而提高了优化设计点的稳定性,使得最终设计方案更加合理,设计效率也明显提高。3. Overcome the problem that the original optimization result relies too much on the value range of design variables, and it is not necessary to repeatedly adjust the design range of each arc size (R 1 , R 2 ) of the special-shaped hole in the optimization design, thereby improving the optimization design point. The stability makes the final design scheme more reasonable, and the design efficiency is also significantly improved.

附图说明Description of drawings

图1是加权系数k1,k2,k3的动态确定流程图;Fig. 1 is a flow chart of dynamic determination of weighting coefficients k 1 , k 2 , k 3 ;

图2是基于多目标两变量优化模型下得到的异型孔轮廓;Figure 2 is the contour of the special-shaped hole obtained based on the multi-objective two-variable optimization model;

其中:1-两变量异型孔主圆弧R1;2-两变量异型孔过渡圆弧R2;3-传统螺栓通孔(基圆)。Among them: 1- the main arc R 1 of the two-variable special-shaped hole; 2- the transition arc R 2 of the two-variable special-shaped hole; 3- the traditional bolt through hole (base circle).

具体实施方式Detailed ways

下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。Below in conjunction with the accompanying drawings and specific embodiments, the present invention will be further clarified. It should be understood that these examples are only used to illustrate the present invention and are not used to limit the scope of the present invention. Modifications in the form of valence all fall within the scope defined by the appended claims of the present application.

一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,该模型针对航空发动机高压涡轮盘上的异型螺栓通孔(以下简称异型孔)轮廓优化问题,提出了一种稳健设计方法,即在优化设计中,采用(a)减小孔边应力至任务规定值和(b)控制异型孔轮廓变化幅度(即避免异型孔轮廓与原圆形轮廓差异过大)等优化目标,以此来获得合理、稳健的设计方案。在优化模型中,通过对每个优化目标进行合理的加权求和处理,建立了异型孔综合优化设计模型。在综合优化模型中,由于每个目标的权重将随异型孔设计变量取值范围的变化而变化,为此专利发明者设计了一种可保持不同设计目标对综合优化目标贡献均衡的权重系数动态调整方法,使得与应力控制目标有关的权重系数(k1)与异型孔轮廓设计变量取值范围相关联,在后者发生变化时,k1可自动调整,因而获得了多目标条件下的异型孔均衡优化设计效果,从根本上消除了异型孔各设计变量(即圆弧尺寸)优化结果对设计变量取值范围的依赖性,可快速优化得到满足强度设计任务要求条件且的轮廓尺寸变化最小的最佳异型孔轮廓。本发明使得两变量异型孔轮廓设计可按任务给出的孔边应力目标值进行优化设计,同时无需对异型孔各设计变量尺寸取值范围反复试探调整,设计结果具有较好的一致性。A robust design method for two-parameter special-shaped bolt holes based on minimum contour changes. This model proposes a robust design method for the contour optimization problem of special-shaped bolt through holes (hereinafter referred to as special-shaped holes) on aero-engine high-pressure turbine disks. In the optimization design, optimization objectives such as (a) reducing the stress on the edge of the hole to the specified value of the task and (b) controlling the variation range of the contour of the special-shaped hole (that is, to avoid the excessive difference between the contour of the special-shaped hole and the original circular contour) are adopted to obtain Reasonable and robust design. In the optimization model, a comprehensive optimization design model for special-shaped holes is established by performing a reasonable weighted summation process for each optimization objective. In the comprehensive optimization model, since the weight of each objective will change with the value range of the special-shaped hole design variables, the patent inventor designed a dynamic weight coefficient that can keep the contribution of different design objectives to the comprehensive optimization objective balanced. The adjustment method makes the weight coefficient (k 1 ) related to the stress control target associated with the value range of the contour design variables of the special-shaped hole. When the latter changes, k 1 can be automatically adjusted, thus obtaining the special-shaped hole under the multi-objective condition. The effect of hole balance optimization design fundamentally eliminates the dependence of the optimization results of each design variable (ie arc size) of the special-shaped hole on the value range of the design variables, and can be quickly optimized to meet the requirements of the strength design task with the smallest change in the contour size. optimum hole profile. The invention enables the two-variable special-shaped hole contour design to be optimally designed according to the target value of the hole edge stress given by the task, and at the same time, it is not necessary to repeatedly test and adjust the size range of each design variable of the special-shaped hole, and the design results have good consistency.

具体包括以下步骤:Specifically include the following steps:

步骤1.根据任务给定目标值、主圆弧半径、过渡圆弧半径、以应力下降目标、异型孔主圆弧尺寸变化目标以及过渡圆弧尺寸变化目标建立异型孔多目标优化模型。Step 1. Establish a multi-objective optimization model for special-shaped holes according to the given target value of the task, the radius of the main arc, the radius of the transition arc, the target of stress reduction, the target of the size of the main arc of the special-shaped hole and the change of the size of the transition arc.

具体包括以下步骤:Specifically include the following steps:

步骤11,根据任务给定目标值、主圆弧半径、过渡圆弧半径建立无量纲应力下降幅目标函数,f1max)为无量纲应力下降幅目标函数,需保证异型孔孔边最大主应力足够逼近于事先设定(任务给定)的某一目标值σo。其表达式如下:Step 11: According to the given target value of the task, the radius of the main arc, and the radius of the transition arc, a dimensionless stress reduction range objective function is established, f 1max ) is the dimensionless stress reduction range objective function, and it is necessary to ensure that the maximum hole edge The principal stress is sufficiently close to a certain target value σ o set in advance (task given). Its expression is as follows:

f1max)=mindf 1max )=mind

d表示f1(R)所描述的孔边最大主应力σmax与事先设定的应力目标值σo的无量纲差值。d represents the dimensionless difference between the maximum principal stress σ max described by f 1 (R) and the pre-set stress target value σ o .

Figure GDA0002195969590000051
Figure GDA0002195969590000051

f1max)为无量纲应力下降幅目标函数,σmax为孔边最大主应力,σmax=σmax(R1,R2),R1为主圆弧半径,R2为过渡圆弧半径,σo表示任务给定目标值。在具体实施时,σo应力目标值需根据设计任务目标及工作条件确定f 1max ) is the objective function of the dimensionless stress reduction range, σ max is the maximum principal stress at the edge of the hole, σ maxmax (R 1 , R 2 ), R 1 is the main arc radius, and R 2 is the transition circle Arc radius, σ o represents the task given target value. In the specific implementation, the σ o stress target value should be determined according to the design task objectives and working conditions

步骤12,根据原螺栓通孔半径、主圆弧半径建立主圆弧半径无量纲异型孔轮廓控制目标函数,该目标函数将使得优化获得的异型孔轮廓与原轮廓差异尽可能小,从而得到更加稳健的设计点。以原螺栓通孔半径R*为基准。其表达式如下:Step 12, establish the main arc radius dimensionless special-shaped hole contour control objective function according to the original bolt through-hole radius and the main arc radius. Robust design point. Based on the original bolt through hole radius R*. Its expression is as follows:

Figure GDA0002195969590000052
Figure GDA0002195969590000052

其中,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,R1为主圆弧半径,R*表示原螺栓通孔半径。Among them, f 2 (R 1 ) is the main arc radius of the dimensionless special-shaped hole contour control objective function, R 1 is the main arc radius, and R* represents the original bolt through hole radius.

而目标函数f2(R1),f3(R2)保证满足目标f1(R)的前提下,获得的异型孔轮廓尺寸R(R1,R2)尽可能接近传统设计螺栓通孔轮廓,从而获得更加稳健的设计点。On the premise that the objective functions f 2 (R 1 ) and f 3 (R 2 ) can satisfy the objective f 1 (R), the obtained profile size R(R 1 , R 2 ) of the special-shaped hole is as close as possible to the traditionally designed bolt through hole contour for a more robust design point.

f2(R1)=minδ1,f3(R2)=minδ2f 2 (R 1 )=minδ 1 , f 3 (R 2 )=minδ 2 ;

其中,δ1,δ2分别表征异型孔主圆弧R1和过渡圆弧半径R2与圆螺栓孔基圆半径的无量纲差,具体如下。Among them, δ 1 and δ 2 respectively represent the dimensionless difference between the main arc R 1 and transition arc radius R 2 of the special-shaped hole and the base circle radius of the bolt hole, as follows.

步骤13,根据原螺栓通孔半径、过渡圆弧半径建立过渡圆弧半径无量纲异型孔轮廓控制目标函数,该目标函数将使得优化获得的异型孔轮廓与原轮廓差异尽可能小,从而得到更加稳健的设计点。以原螺栓通孔半径R*为基准。其表达式如下:Step 13: According to the original bolt through hole radius and the transition arc radius, establish the transition arc radius dimensionless special-shaped hole contour control objective function, which will make the difference between the optimized special-shaped hole contour and the original contour as small as possible, so as to obtain more Robust design point. Based on the original bolt through hole radius R*. Its expression is as follows:

Figure GDA0002195969590000061
Figure GDA0002195969590000061

其中,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数,R2为过渡圆弧半径,R*表示原螺栓通孔半径。Among them, f 3 (R 2 ) is the transition arc radius dimensionless hole contour control objective function, R 2 is the transition arc radius, and R* represents the original bolt through hole radius.

优化设计中,确定设计变量取值范围时,R1max(主圆弧半径上界)可根据实际螺栓孔直径取为8~20R*(R*为基圆半径);R2min为过渡圆弧半径下界,取决于实际加工条件下的最小半径,可取实际用最小直径刀具半径。In the optimization design, when determining the value range of the design variables, R 1max (the upper bound of the main arc radius) can be taken as 8~20R* (R* is the base circle radius) according to the actual bolt hole diameter; R 2min is the transition arc radius The lower bound depends on the minimum radius under actual machining conditions, and the actual minimum diameter tool radius can be taken.

步骤14,根据步骤11建立的无量纲应力下降幅目标函数、步骤12建立的主圆弧半径无量纲异型孔轮廓控制目标函数、步骤13建立的过渡圆弧半径无量纲异型孔轮廓控制目标函数建立异型孔多目标优化模型,其表达式如下:Step 14, establish according to the objective function of the dimensionless stress reduction amplitude established in step 11, the main arc radius dimensionless special-shaped hole contour control objective function established in step 12, and the transition arc radius dimensionless special-shaped hole contour control objective function established in step 13. The multi-objective optimization model for special-shaped holes is expressed as follows:

minf(y1=f1max),y2=f2(R1),y3=f3(R2))minf(y 1 =f 1max ),y 2 =f 2 (R 1 ),y 3 =f 3 (R 2 ))

求解具体优化问题时,所有设计点Ri(R1,R2)对应目标函数的f1,f2,f3可构成一个判据空间,用F表示,如下:When solving a specific optimization problem, all design points R i (R 1 , R 2 ) correspond to f 1 , f 2 , f 3 of the objective function to form a criterion space, which is represented by F, as follows:

F={fi(f1,f2,f3)∈R3|f1=f1maxi),f2=f2(R1),f3=f3(R2)}F={f i (f 1 , f 2 , f 3 )∈R 3 |f 1 =f 1maxi ),f 2 =f 2 (R 1 ),f 3 =f 3 (R 2 )}

fi=(f1,f2,f3)为三个目标函数值组成的向量。寻优过程即为对判据空间中所有向量fi的评估过程。f i =(f 1 , f 2 , f 3 ) is a vector composed of three objective function values. The optimization process is the evaluation process of all vectors f i in the criterion space.

σmax=σmax(R1,R2),R1,R2分别代表主圆弧和过渡圆弧,它们代表了设计点R;f1max)为与设计点R(R1,R2)对应的应力下降幅目标函数;f2(R1)为主圆弧R1变化控制目标函数;f3(R2)为过渡圆弧R2变化控制目标函数。R*为基圆(即传统设计螺栓通孔)半径,R1max为主圆弧半径上界,R2min为过渡圆弧半径下界。σ maxmax (R 1 , R 2 ), R 1 , R 2 represent the main arc and the transition arc respectively, which represent the design point R; f 1max ) is the same as the design point R(R 1 , R 2 ) corresponds to the objective function of the stress drop amplitude; f 2 (R 1 ) is the objective function of the change control of the main arc R 1 ; f 3 (R 2 ) is the objective function of the change of the transition arc R 2 . R* is the radius of the base circle (that is, the traditional design bolt through hole), R 1max is the upper bound of the main arc radius, and R 2min is the lower bound of the transition arc radius.

异型孔多目标优化模型满足:The multi-objective optimization model of special-shaped holes satisfies:

(1)降低孔边应力至强度设计任务给定目标;(1) Reduce the stress on the edge of the hole to the given target of the strength design task;

(2)孔轮廓尺寸改变最小;(2) The size of the hole contour changes the smallest;

(3)有效降低优化所得的异型孔轮廓尺寸对设计变量取值范围的依赖性,提高设计的鲁棒性。(3) Effectively reduce the dependence of the contour size of the special-shaped hole obtained by optimization on the value range of the design variables, and improve the robustness of the design.

步骤2.根据主圆弧半径上界、过渡圆弧半径下界确定步骤1中的异型孔多目标优化模型的尺寸约束。Step 2. Determine the size constraints of the multi-objective optimization model of the special-shaped hole in step 1 according to the upper bound of the main arc radius and the lower bound of the transition arc radius.

其尺寸约束为:Its size constraints are:

Figure GDA0002195969590000071
Figure GDA0002195969590000071

其中,R1max为主圆弧半径上界,可在合理范围内任意取值,R2min为过渡圆弧半径下界,R1为主圆弧半径,R2为过渡圆弧半径,R*表示基圆,即传统设计(原)螺栓通孔半径。Among them, R 1max is the upper bound of the main arc radius, which can be any value within a reasonable range, R 2min is the lower bound of the transition arc radius, R 1 is the main arc radius, R 2 is the transition arc radius, and R* represents the base Circle, i.e. traditional design (original) bolt through hole radius.

步骤3.引用加权系数,将异型孔多目标优化模型转化为异型孔单一综合目标优化模型,即对目标f1max)、f2(R1)和f3(R2)进行加权求和处理,将多目标优化问题简化为单一综合目标。Step 3. Convert the multi-objective optimization model of special-shaped holes into a single comprehensive objective optimization model of special-shaped holes by referring to the weighting coefficient, that is, weighted calculation of the objectives f 1max ), f 2 (R 1 ) and f 3 (R 2 ) and processing, reducing the multi-objective optimization problem into a single comprehensive objective.

异型孔单一综合目标优化模型为:The single comprehensive objective optimization model for special-shaped holes is:

f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2)。f(f 1 , f 2 , f 3 )=k 1 ·f 1max )+k 2 ·f 2 (R 1 )+k 3 ·f 3 (R 2 ).

其中,k1为量纲应力下降幅目标函数的加权系数,f1max)为无量纲应力下降幅目标函数,k2为主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,k3为过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数。Among them, k 1 is the weighting coefficient of the objective function of dimensionless stress reduction amplitude, f 1max ) is the objective function of dimensionless stress reduction amplitude, k 2 is the weighting coefficient of the main arc radius dimensionless special-shaped hole contour control objective function, f 2 (R 1 ) the main arc radius dimensionless special-shaped hole contour control objective function, k 3 is the weighting coefficient of the transition arc radius dimensionless special-shaped hole contour control objective function, f 3 (R 2 ) is the transition arc radius infinite The contour control objective function of the special-shaped hole.

本发明由三个新目标函数构成的多目标优化设计模型,可以同时对孔边最大主应力下降幅和异型孔轮廓尺寸变化进行控制,使得孔边应力降低幅满足既定目标,孔型设计更为稳健,孔型参数更为稳定。The multi-objective optimization design model composed of three new objective functions in the present invention can control the maximum principal stress reduction range of the hole edge and the profile size change of the special-shaped hole at the same time, so that the reduction range of the hole edge stress can meet the set target, and the hole shape design is more efficient. Robust, the pass parameters are more stable.

通过加权系数k1,k2,k3可调整三个单一优化目标对综合优化目标的贡献。The contribution of the three single optimization objectives to the comprehensive optimization objective can be adjusted by the weighting coefficients k 1 , k 2 , and k 3 .

如附图1所示,加权系数的确定方法如下:As shown in accompanying drawing 1, the determination method of weighting coefficient is as follows:

量纲应力下降幅目标函数的加权系数k1、主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k2、过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k3随R1主圆弧半径的取值范围以及σo应力目标值的不同而动态变化,其中:主圆弧半径R1和过渡圆弧半径R2取值范围变化时,对孔边应力的影响不同。The weighting coefficient k 1 of the objective function of the dimensional stress reduction amplitude, the weighting coefficient k 2 of the main arc radius dimensionless special-shaped hole contour control objective function, the weighting coefficient k 3 of the transition arc radius dimensionless special-shaped hole contour control objective function varies with R 1 The value range of the main arc radius and the target value of σ o stress vary dynamically, among which: when the value range of the main arc radius R 1 and the transition arc radius R 2 changes, the impact on the hole edge stress is different.

R1max变化时:When R 1max changes:

k2=0.008~0.012k 2 =0.008~0.012

k2和k3之间满足:Between k 2 and k 3 :

k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max k 3 =g 1 (k 2 ,R 1max )=(1.25 to 2.1)·k 2 ·R 1max

k1与k2,k3之间满足:Between k 1 and k 2 , k 3 satisfies:

Figure GDA0002195969590000081
Figure GDA0002195969590000081

其中,σ1r为原圆形螺栓通孔孔边最大主应力值。Among them, σ 1r is the maximum principal stress value of the original circular bolt through hole edge.

生成的加权系数k1,k2,k3可使得各设计参量变化更为均衡,获得稳健的设计点。The generated weighting coefficients k 1 , k 2 , and k 3 can make the variation of each design parameter more balanced, and obtain a robust design point.

步骤4.根据定轮廓变化值和允许的误差对异型孔单一综合目标优化模型进行优化,获得异型孔各设计变量,完成异型孔的设计。Step 4. Optimize the single comprehensive target optimization model of the special-shaped hole according to the change value of the fixed contour and the allowable error, obtain each design variable of the special-shaped hole, and complete the design of the special-shaped hole.

对异型孔单一综合目标优化模型进行优化时,只有在优化后的异型孔孔边应力达到既定要求,即|σmaxo|≤ε,ε为允许误差。且异型螺栓孔与原螺栓圆形通孔相比,轮廓变化小于给定轮廓变化值时,设计参数或变量才是满足要求的优化值,在所给的多目标优化模型中,要求三个单一目标对综合设计结果贡献均衡。由于加权系数k1,k2,k3与设计变量的取值范围相关,为此给出了加权系数k1,k2,k3动态生成办法,使得k1,k2,k3值可随R1max取值范围的不同,调整取值,以满足三个单一优化目标之间的均衡,最终得到满足强度设计要求、且更为稳健、稳定的设计点,保证三者具有恰当的影响因数,从而保证了最终最优解的鲁棒性。When optimizing the single comprehensive objective optimization model of the special-shaped hole, only when the optimized special-shaped hole hole edge stress meets the established requirements, that is, |σ maxo |≤ε, ε is the allowable error. And compared with the circular through hole of the original bolt, when the contour change of the special-shaped bolt hole is smaller than the given contour change value, the design parameters or variables are the optimized values that meet the requirements. In the given multi-objective optimization model, three single parameters are required. Objectives contribute equally to the overall design results. Since the weighting coefficients k 1 , k 2 , k 3 are related to the value range of the design variables, a dynamic generation method for the weighting coefficients k 1 , k 2 , k 3 is given, so that the values of k 1 , k 2 , k 3 can be With the different value ranges of R 1max , adjust the value to meet the balance between the three single optimization objectives, and finally obtain a more robust and stable design point that meets the strength design requirements, ensuring that the three have appropriate influencing factors , thus ensuring the robustness of the final optimal solution.

本发明中提出的多目标优化设计方法可在满足应力设计目标(任务值)的基础上,获得孔轮廓改变更小、更加稳健的孔型设计方案;此外,优化设计中,在异型孔圆弧尺寸允许变化范围明显变化时,仍能保持稳定的优化解,使得最终优化设计结果摆脱了对圆弧设计变量设定范围值的依赖,提高了设计效率。该多目标优化模型下得到的两参量异型孔结构如附图2所示。The multi-objective optimization design method proposed in the present invention can obtain a hole shape design scheme with smaller hole profile changes and more robust on the basis of satisfying the stress design objective (task value); in addition, in the optimization design, the special-shaped hole arc When the allowable variation range of the size changes significantly, the stable optimal solution can still be maintained, so that the final optimal design result can get rid of the dependence on the setting range value of the arc design variable and improve the design efficiency. The two-parameter shaped pore structure obtained under the multi-objective optimization model is shown in FIG. 2 .

本发明提供了一种新的多目标优化设计模型,可按任务规定的孔边应力下降幅给出稳健、稳定的设计方案。该多目标优化模型由(a)应力下降目标、(b)异型孔主圆弧尺寸R1变化目标和(c)过渡圆弧尺寸R2变化目标构成;多目标优化模型中,三个单一目标(即(a)、(b)和(c))对综合目标的贡献应满足均衡条件;所提供的新模型及算法可避免优化设计时繁琐、低效的边界试算,提高设计的合理性及稳定性,同时大大提高设计效率。The invention provides a new multi-objective optimization design model, which can give a robust and stable design scheme according to the drop range of the hole edge stress specified by the task. The multi-objective optimization model consists of (a) the stress reduction objective, (b) the change objective of the main arc size R1 of the special-shaped hole and (c) the change objective of the transition arc size R2 ; in the multi-objective optimization model, three single objectives (i.e. (a), (b) and (c)) the contribution to the comprehensive goal should satisfy the equilibrium condition; the new model and algorithm provided can avoid the tedious and inefficient boundary trial calculation when optimizing the design, and improve the rationality of the design and stability, while greatly improving design efficiency.

以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only the preferred embodiment of the present invention, it should be pointed out that: for those skilled in the art, without departing from the principle of the present invention, several improvements and modifications can also be made, and these improvements and modifications are also It should be regarded as the protection scope of the present invention.

Claims (5)

1.一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,其特征在于,包括以下步骤:1. a robust design method for two-parameter special-shaped bolt holes based on minimum profile changes, is characterized in that, comprises the following steps: 步骤1.根据任务给定目标值、主圆弧半径、过渡圆弧半径、以应力下降目标、异型孔主圆弧尺寸变化目标以及过渡圆弧尺寸变化目标建立异型孔多目标优化模型;Step 1. According to the given target value of the task, the radius of the main arc, the radius of the transition arc, the target of stress reduction, the target of the size of the main arc of the special-shaped hole and the change of the size of the transition arc, the multi-objective optimization model of the special-shaped hole is established; 异型孔多目标优化模型建立方法包括以下步骤:The method for establishing a multi-objective optimization model for special-shaped holes includes the following steps: 步骤11,根据任务给定目标值、主圆弧半径、过渡圆弧半径建立无量纲应力下降幅目标函数;Step 11, establish a dimensionless stress drop range objective function according to the given target value of the task, the radius of the main arc, and the radius of the transition arc; 无量纲应力下降幅目标函数的表达式如下:The expression of the objective function of dimensionless stress drop amplitude is as follows:
Figure FDA0002195969580000011
Figure FDA0002195969580000011
f1max)为无量纲应力下降幅目标函数,σmax为孔边最大主应力,σmax=σmax(R1,R2),R1为主圆弧半径,R2为过渡圆弧半径,σo表示任务给定目标值;f 1max ) is the objective function of the dimensionless stress reduction range, σ max is the maximum principal stress at the edge of the hole, σ maxmax (R 1 , R 2 ), R 1 is the main arc radius, and R 2 is the transition circle Arc radius, σ o represents the given target value of the task; 步骤12,根据原螺栓通孔半径、主圆弧半径建立主圆弧半径无量纲异型孔轮廓控制目标函数;Step 12, establish the main arc radius dimensionless special-shaped hole contour control objective function according to the original bolt through hole radius and the main arc radius; 所述步骤12中主圆弧半径无量纲异型孔轮廓控制目标函数的表达式如下:In described step 12, the expression of the main arc radius dimensionless special-shaped hole contour control objective function is as follows:
Figure FDA0002195969580000012
Figure FDA0002195969580000012
其中,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,R1为主圆弧半径,R*表示原螺栓通孔半径;Among them, f 2 (R 1 ) is the main arc radius of the dimensionless special-shaped hole contour control objective function, R 1 is the main arc radius, and R* represents the original bolt through hole radius; 步骤13,根据原螺栓通孔半径、过渡圆弧半径建立过渡圆弧半径无量纲异型孔轮廓控制目标函数;Step 13, according to the original bolt through hole radius and the transition arc radius, establish the transition arc radius dimensionless special-shaped hole contour control objective function; 过渡圆弧半径无量纲异型孔轮廓控制目标函数的表达式如下:The expression of the objective function of the transition arc radius dimensionless hole contour control is as follows:
Figure FDA0002195969580000013
Figure FDA0002195969580000013
其中,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数,R2为过渡圆弧半径,R*表示原螺栓通孔半径;Among them, f 3 (R 2 ) is the transition arc radius dimensionless hole profile control objective function, R 2 is the transition arc radius, and R* represents the original bolt through hole radius; 步骤14,根据步骤11建立的无量纲应力下降幅目标函数、步骤12建立的主圆弧半径无量纲异型孔轮廓控制目标函数、步骤13建立的过渡圆弧半径无量纲异型孔轮廓控制目标函数建立异型孔多目标优化模型;Step 14, establish according to the objective function of the dimensionless stress reduction amplitude established in step 11, the main arc radius dimensionless special-shaped hole contour control objective function established in step 12, and the transition arc radius dimensionless special-shaped hole contour control objective function established in step 13. Special-shaped hole multi-objective optimization model; 异型孔多目标优化模型的表达式如下:The expression of the multi-objective optimization model for special-shaped holes is as follows: minf(y1=f1max),y2=f2(R1),y3=f3(R2));minf(y 1 =f 1max ),y 2 =f 2 (R 1 ),y 3 =f 3 (R 2 )); 步骤2.根据主圆弧半径上界、过渡圆弧半径下界确定步骤1中的异型孔多目标优化模型的尺寸约束;Step 2. Determine the size constraints of the multi-objective optimization model of the special-shaped hole in step 1 according to the upper bound of the main arc radius and the lower bound of the transition arc radius; 步骤3.引用加权系数,将异型孔多目标优化模型转化为异型孔单一综合目标优化模型;Step 3. Quoting the weighting coefficient, convert the multi-objective optimization model of the special-shaped hole into a single comprehensive target optimization model of the special-shaped hole; 步骤4.根据定轮廓变化值和允许的误差对异型孔单一综合目标优化模型进行优化,获得异型孔各设计变量,完成异型孔的设计。Step 4. Optimize the single comprehensive target optimization model of the special-shaped hole according to the change value of the fixed contour and the allowable error, obtain each design variable of the special-shaped hole, and complete the design of the special-shaped hole.
2.根据权利要求1所述的基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,其特征在于:所述步骤2中建立的尺寸约束为:2. The robust design method for two-parameter special-shaped bolt holes based on minimum profile changes according to claim 1, wherein the size constraints established in the step 2 are:
Figure FDA0002195969580000021
Figure FDA0002195969580000021
其中,R1max为主圆弧半径上界,R2min为过渡圆弧半径下界,R1为主圆弧半径,R2为过渡圆弧半径,R*表示原螺栓通孔半径。Among them, R 1max is the upper bound of the main arc radius, R 2min is the lower bound of the transition arc radius, R 1 is the main arc radius, R 2 is the transition arc radius, and R* represents the original bolt through hole radius.
3.根据权利要求2所述的基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,其特征在于:所述步骤3中的异型孔单一综合目标优化模型为:3. The robust design method for two-parameter special-shaped bolt holes based on minimum profile changes according to claim 2, wherein: the special-shaped hole single comprehensive target optimization model in the step 3 is: f(f1,f2,f3)=k1·f1max)+k2·f2(R1)+k3·f3(R2);f(f 1 , f 2 , f 3 )=k 1 ·f 1max )+k 2 ·f 2 (R 1 )+k 3 ·f 3 (R 2 ); 其中,k1为量纲应力下降幅目标函数的加权系数,f1max)为无量纲应力下降幅目标函数,k2为主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f2(R1)主圆弧半径无量纲异型孔轮廓控制目标函数,k3为过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数,f3(R2)为过渡圆弧半径无量纲异型孔轮廓控制目标函数。Among them, k 1 is the weighting coefficient of the objective function of dimensionless stress reduction amplitude, f 1max ) is the objective function of dimensionless stress reduction amplitude, k 2 is the weighting coefficient of the main arc radius dimensionless special-shaped hole contour control objective function, f 2 (R 1 ) the main arc radius dimensionless special-shaped hole contour control objective function, k 3 is the weighting coefficient of the transition arc radius dimensionless special-shaped hole contour control objective function, f 3 (R 2 ) is the transition arc radius infinite The contour control objective function of the special-shaped hole. 4.根据权利要求3所述的基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,其特征在于:所述步骤3中加权系数的确定方法如下:4. The robust design method for two-parameter special-shaped bolt holes based on minimum profile changes according to claim 3, wherein the method for determining the weighting coefficient in the step 3 is as follows: 量纲应力下降幅目标函数的加权系数k1、主圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k2、过渡圆弧半径无量纲异型孔轮廓控制目标函数的加权系数k3随R1主圆弧半径的取值范围以及σo应力目标值的不同而动态变化,其中:主圆弧半径R1和过渡圆弧半径R2取值范围变化时,对孔边应力的影响不同;The weighting coefficient k 1 of the objective function of the dimensional stress reduction amplitude, the weighting coefficient k 2 of the main arc radius dimensionless special-shaped hole contour control objective function, the weighting coefficient k 3 of the transition arc radius dimensionless special-shaped hole contour control objective function varies with R 1 The value range of the main arc radius and the target value of σ o stress vary dynamically, among which: when the value ranges of the main arc radius R 1 and the transition arc radius R 2 change, the impact on the stress on the edge of the hole is different; R1max变化时:When R 1max changes: k2=0.008~0.012k 2 =0.008~0.012 k2和k3之间满足:Between k 2 and k 3 : k3=g1(k2,R1max)=(1.25~2.1)·k2·R1max k 3 =g 1 (k 2 ,R 1max )=(1.25 to 2.1)·k 2 ·R 1max k1与k2,k3之间满足:Between k 1 and k 2 , k 3 satisfies:
Figure FDA0002195969580000031
Figure FDA0002195969580000031
其中,σ1r为原圆形螺栓通孔孔边最大主应力值。Among them, σ 1r is the maximum principal stress value of the original circular bolt through hole edge.
5.根据权利要求4所述的基于最小轮廓变化的两参量异型螺栓孔稳健设计方法,其特征在于:所述步骤4中对异型孔单一综合目标优化模型进行优化时,只有在优化后的异型孔孔边应力达到既定要求,即|σmaxo|≤ε,ε为允许误差;且异型螺栓孔与原螺栓圆形通孔相比,轮廓变化小于给定轮廓变化值时,设计参数或变量才是满足要求的优化值。5. The method for robust design of two-parameter special-shaped bolt holes based on minimum contour changes according to claim 4, characterized in that: in the step 4, when optimizing a single comprehensive target optimization model for special-shaped holes, only the optimized special-shaped The hole edge stress meets the established requirements, that is, |σ maxo |≤ε, ε is the allowable error; and when the contour change of the special-shaped bolt hole is smaller than the given contour change value compared with the original bolt circular through hole, the design parameters Or variable is the optimized value that satisfies the requirement.
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 南京航空航天大学 A five-parameter multi-arc shaped hole balanced 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
CN109376418B (en) An aerodynamic robust optimization design method for small wind turbine airfoils suitable for turbulent flow conditions
CN103277245B (en) Large-thickness blunt-trailing-edge wind-power airfoil profiles and a design method thereof
CN103136422B (en) The integrated intermediate gauge Airfoil Design method be combined with B-spline of aerofoil profile
CN105673322B (en) Realize the variable element Nonlinear Feedback Control Method of wind energy conversion system MPPT controls
CN105404743B (en) The wind mill airfoil design method that B-spline is combined with the smooth continuity of curvature
CN106156409B (en) Two-parameter special-shaped bolt hole robust design method based on minimum contour change
CN109145506A (en) The optimum design method of aerofoil profile on the outside of a kind of high aeroperformance low noise level wind energy conversion system
CN105425612A (en) Preferred method of water turbine adjustment system control parameter
CN107479365A (en) The fuzzy synthetic appraisement method of PID suspension controllers based on genetic algorithm
CN112329170B (en) Blade airfoil section optimization design method for tidal current energy water turbine
CN111859801B (en) Method for designing stall-controlled wind turbine airfoil
CN106611951B (en) A kind of power distribution network multi-objective reactive optimization method based on two stages particle group optimizing
CN103047079B (en) Efficient low-load wing section special for horizontal shaft wind turbine blade and designing method of wind section
CN110361974B (en) Optimization method of turbine speed control system based on BP-FOA hybrid algorithm
CN108708823B (en) On-line optimization method and system for optimal gain parameters of wind turbines
CN105279579B (en) A kind of preferred method of turbine-generator units excitation system pid control parameter
CN109783912B (en) Optimized design method for lift type vertical axis wind turbine blade airfoil profile
CN107194122B (en) Improvement method for multi-working-condition pneumatic optimization of variable speed fan blade
CN106971023B (en) Design method of special-shaped hole structure of roulette based on hyperelliptic curve
CN105046000B (en) The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable
CN112054549A (en) Method for improving stability of direct-current power transmission system
CN113654143B (en) Air-conditioning water system pump valve linkage control method based on neural network
CN117034753A (en) Decision variable grouping multi-objective optimization method based on extreme point pareto front edge
CN109802440B (en) Offshore wind farm equivalence method, system and device based on wake effect factor
CN110970936A (en) Method for calculating primary frequency modulation performance of deep peak shaving unit

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

Granted publication date: 20200114

CF01 Termination of patent right due to non-payment of annual fee