CN105046000B - The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable - Google Patents
The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable Download PDFInfo
- Publication number
- CN105046000B CN105046000B CN201510427969.1A CN201510427969A CN105046000B CN 105046000 B CN105046000 B CN 105046000B CN 201510427969 A CN201510427969 A CN 201510427969A CN 105046000 B CN105046000 B CN 105046000B
- Authority
- CN
- China
- Prior art keywords
- design
- variable
- discrete
- value
- standard
- 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
Links
- 238000013461 design Methods 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000012546 transfer Methods 0.000 claims abstract description 9
- 238000013507 mapping Methods 0.000 claims abstract description 6
- 238000004519 manufacturing process Methods 0.000 claims abstract 2
- 230000006870 function Effects 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 6
- 238000009510 drug design Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
Landscapes
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
The invention discloses a kind of more circular-arc abnormal-shape hole Equalization Design methods of two variables based on engineering discrete variable and model, this method establishes the standard alternate pond of each design variable according to the arc radius standard series of national regulations first, to store arc radius standard value；Then, rounding is carried out to former continuous variable with self-defined transfer function；The design space combined using marginal check module to discrete variable is controlled, and for crossing the border, parameter is rejected；Last Optimized Iterative, before optimization, discrete design presses size order numbering；In Optimized Iterative, when selecting new design point, select new design point to number by numbering computing, then determine new design point by the numbering of discrete design point and the mapping relations of actual design point.The present invention can make two variable profiled holes profile design radiuses meet industrial processes standard, and to meet the principle of equal strength during the design of aero-engine high-pressure turbine disk connection unit, improve design efficiency.
Description
Technical field
The invention belongs to the Optimal Structure Designing of aero-turbine rotor part, more particularly to turbine disk bolt hole are different
The contour shape Optimization Design and model of type hole.
Background technology
Side upper bolt hole is installed before aero-engine high-pressure turbine disk and mainly uses borehole structure, long-term work high pressure,
Under hot environment, easily there is stress concentration phenomenon in circular hole hole side, produces hole-edge crack, has a strong impact on the turbine disk or even entirely start
The service life of machine.
For this problem, turbine disk upper bolt hole is improved to a kind of special-shaped pore structure by the engines of CFM56- III first, south
Capital Aero-Space university computer MSR Information system is made that numerous studies to the contour structure feature and optimization method of the profiled holes, provides
The Derivation of Mathematical Model of profiled holes, and propose Optimization Design.Current optimization method, which can provide, to be designed profiled holes
Recommendation arc radius value with reference value, around-the-hole stress can be reduced to a certain extent, improve around-the-hole stress distribution feelings
Condition.But still suffer from following deficiency：
1. the profiled holes profile of design gained excessively relies on the span of the design variable of setting, it is unfavorable for engineering and sets
Meter.
It is difficult comparatively fast to converge to optimal solution when being designed using former method 2. being difficult to the design point stablized, takes longer.
3. former design method carries out optimizing for continuous variable, resulting optimum radius size does not meet actual industrial and added
Work standard, come to actual processing and detection band difficult.
The content of the invention
The shortcomings that for above-mentioned prior art, it is an object of the invention to provide a kind of two variables based on engineering discrete variable
More circular-arc abnormal-shape hole Equalization Design methods, enable to profiled holes arc radius design result to meet industrial standard.
To achieve the above object, the technical solution adopted by the present invention is：
A kind of more circular-arc abnormal-shape hole Equalization Design methods of two variables based on engineering discrete variable, comprise the following steps：
(1) sliding-model control of design variable
The standard alternate pond of each design variable is established according to the arc radius standard series of national Specification first, to
Store arc radius standard value；Then, write self-defined transfer function with VB language and rounding is carried out to former continuous variable；Ensure
Each design variable sample is the centrifugal pump randomly generated in radii criteria sequence pond；
(2) control of discrete variable design space
The design space combined using marginal check module to discrete variable is controlled, and for crossing the border, parameter is picked
Remove；
As a further improvement, step (3) Optimized Iterative process is also included after step (2), it is specially：
Before optimization, discrete design presses size order numbering；In Optimized Iterative, when selecting new design point, pass through numbering
Computing selects new design point to number, and is then determined by the numbering of discrete design point with the mapping relations of actual design point
New design point, by the comparison of design point, pick out more preferable new design point.
In step (1), rounding is carried out to former continuous variable with following self-defined transfer functions, the self-defined transfer function
Mathematical modeling is：
Xt=min { Xs|Xs≥Xi}
Wherein, XtRepresent the discrete heat sources value after rounding, XsRepresent the standard value in standard alternate pond, XiRepresent with
Continuous variable value caused by machine.
When carrying out sliding-model control with the function pair design variable, the continuous variable value and standard alternate that will randomly generate
Data are compared one by one in pond, and until approaching standard value that it is not more than and immediate, it is random continuous to take the standard value
Centrifugal pump corresponding to variable.
In step (2), the screening rule of marginal check module is as follows：
If Xt>XupOr Xt<Xlow
Then Xt=Xnew
Wherein, XtRepresent discrete variable combination to be tested, XupIt is the arc radius upper bound artificially set, XlowIt is artificial
Arc radius lower bound, X are setnewFor a variate-value in newly-generated standard design space.
The screening module is by the way that the discrete variable regenerated value is compared with the radius variable bound pre-set
Compared with carrying out judgment variable value whether in rational design space, randomly selected for getting over bound variable in the standard design space
A new discrete variable value substitute.
It is a further object to provide a kind of turbine disk profiled holes Optimized model, with solve optimum results excessively according to
The problem of relying design variable codomain, more sane design knot can be obtained in the case where meeting Structural Strength Design requirement
Fruit.
To achieve the above object, the technical solution adopted by the present invention is：
Equilibrium considers the design feature and the variable profiled holes of the turbine disk two of aircraft engine high pressure rotor linkage unit
Around-the-hole stress changing rule, it is proposed that the new Optimized model of the variable profiled holes of the turbine disk two is as follows
A kind of variable profiled holes Optimized model of the turbine disk two, the model are：
min f(x)
Wherein, R1For main arc radius, R2For transition arc radius, f (x) is hole side principal stress value, σ1maxFor hole side first
Maximum principal stress, σ1stFor restraint stress, δ is a small amount of, Rmin1For main arc radius lower bound, Rmax1For the main arc radius upper bound,
Rmin2For transition arc radius lower bound, Rmax2For the transition arc radius upper bound.
σ1stLateral opening side maximum stress is installed by reference to other components in matched turbine assembly to determine, or to set
Target call is counted to determine；R1And R2Span press Patent No. ZL201310287309.9 Chinese patent " more circular-arc abnormal-shapes
Hole " and actual processing request determine.
The beneficial effects of the invention are as follows：
1) sliding-model control is carried out to design variable, the optimization knot for meeting processing criterion can be directly obtained in design
Fruit, improve the quality of optimal solution.Meanwhile the scale in the alternative pond of sample can be effectively reduced using discrete variable, reduce design
Total operand；
2) Optimized model considered using Equalization Design so that design process more meets actual requirement, not only can effectively press down
Radial design result processed too relies on the situation on variable-value border, and can substantially accelerate the convergence rate of Optimization Progress, carries
High design efficiency；
3) introducing of restraint stress so that around-the-hole stress drops to reasonable value optimization process and just stopped, and coordinates appropriate
Circular arc border parameter setting, can obtain with circular hole overall size similar in optimal solution, and the requirement of stress decrease rate can be met, this
It ensure that the harmony of design.
Brief description of the drawings
Fig. 1 is two variable profiled holes design profile figures；
Fig. 2 is the profiled holes profile for not considering to obtain under Equalization Design thought；
Fig. 3 is the profiled holes profile obtained under the new Optimized model considered based on Equalization Design；
Fig. 4 is the concrete numerical value example of rounding mapping；
Wherein：The main circular arc Rs of variable profiled holes of 1- two1；The variable profiled holes transition arc R of 2- two2；3- bolt circular hole basic circles；4-
Two variable profiled holes transition arcs of Equalization Design are not considered；5- does not consider the two main circular arcs of variable profiled holes of Equalization Design；6-
Consider two variable profiled holes transition arcs of Equalization Design；7- considers the two main circular arcs of variable profiled holes of Equalization Design.
Embodiment
The present invention is further described with reference to embodiment.
1. the implementation of design variable discretization
1) species Discrete is initialized：Each design is established according to the arc radius standard series of national Specification first to become
The standard alternate pond of amount, to store arc radius standard value；Then, self-defined transfer function is write with VB language to connect original
Continuous variable carries out rounding.Principle is：More than initial value, and away from its nearest arc radius preferred value.It ensure that each design
Variable sample is the centrifugal pump randomly generated in radii criteria sequence pond.
The mathematical modeling of self-defined transfer function is：
Xt=min { Xs|Xs≥Xi}
Wherein, XtRepresent the discrete heat sources value after rounding, XsRepresent the standard value in standard alternate pond, XiRepresent with
Continuous variable value caused by machine.
When carrying out sliding-model control with the function pair design variable, the continuous variable value and standard alternate that will randomly generate
Data are compared one by one in pond, and until approaching standard value that it is not more than and immediate, it is random continuous to take the standard value
Centrifugal pump corresponding to variable.
Reference picture 4, by taking concrete numerical value as an example, explain the realization of rounding mapping.
2) discrete variable design space controls：
Marginal check module is introduced, inspection of crossing the border is carried out to the combined sample of discrete variable, it is full to ensure each design point
The discrete variable of sufficient prime design requirement, and without departing from the alternative pond scope of design variable.The screening rule of marginal check module
It is as follows：
If Xt>XupOr Xt<Xlow
Then Xt=Xnew
Wherein, XtRepresent discrete variable combination to be tested, XupIt is the arc radius upper bound artificially set, XlowIt is artificial
Arc radius lower bound, X are setnewFor a variate-value in newly-generated standard design space.
The screening module is by the way that the discrete variable regenerated value is compared with the radius variable bound pre-set
Compared with carrying out judgment variable value whether in rational design space, randomly selected for getting over bound variable in the standard design space
A new discrete variable value substitute.
3) Optimized Iterative process
Before optimization, discrete design presses size order numbering；In Optimized Iterative, when selecting new design point, pass through numbering
Computing selects new design point to number, and is then determined by the numbering of discrete design point with the mapping relations of actual design point
New design point, by the comparison of design point, more preferable new design point is picked out, this method can accelerate to promote Optimized Iterative to enter
Journey.
2. the Optimized model construction method that Equalization Design considers：
Two new variable profiled holes Optimized models of structure are considered based on Equalization Design, introduce new stress constraint condition.Should
The focus of invention is the introducing and determination of restraint stress, and the principal element that it need to consider when setting is as follows：
1) consider for uniform strength design principle, turbine disk installation side bolt hole around-the-hole stress should be not higher than associated with it
Installation side bolt hole around-the-hole stress value (restraint stress), turbine disk bolt hole around-the-hole stress is reduced in restraint stress during design
Nearby.
2) determination of restraint stress value also needs to consider influence of the CONSTRUCTED SPECIFICATION size (chamfering etc.) to profiled holes around-the-hole stress
Suitably adjusted.
The variable profiled holes Optimized model of the turbine disk two is：
minf(x)
Wherein, R1For main arc radius, R2For transition arc radius, f (x) is hole side principal stress value, σ1maxFor hole side first
Maximum principal stress, σ1stFor the restraint stress newly introduced, σ1stIt is lateral opening with reference to other component installations in matched turbine assembly
Side maximum stress determines, or requires to determine with design object；δ is a small amount of, Rmin1For main arc radius lower bound, Rmax1For main circular arc
The radius upper bound, Rmin2For transition arc radius lower bound, Rmax2For the transition arc radius upper bound.R1And R2Span press patent
" more circular-arc abnormal-shape holes " (patent No. ZL2013 1 0287309.9) and actual processing request determine.
Two variable profiled holes optimum results under Different Optimization model are as shown in table 1.
The Different Optimization model of table 1 is to two variable profiled holes optimum results contrast tables
Note：Design variable span is：6mm≤R1≤ 40mm, 2mm≤R2≤5mm；What model 1 and model 2 obtained
Optimum results can meet design requirement.
Its result is respectively as shown in accompanying drawing 2, accompanying drawing 3, it is seen that is not introduced into the profiled holes for optimizing before Equalization Design theory to obtain
Profile is not easy to process and detected closer to a square hole.Latter scheme computing sample number is few, and desin speed can substantially add
It hurry up.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, some improvement can also be done under the premise without departing from the principles of the invention, and these improvement also should be regarded as the guarantor of the present invention
Protect scope.
Claims (2)
1. a kind of more circular-arc abnormal-shape hole Equalization Design methods of two variables based on engineering discrete variable, it is characterised in that：Including with
Lower step：
（1）The sliding-model control of design variable
The standard alternate pond of each design variable is established according to the arc radius standard series of national regulations first, to store circular arc
Radii criteria value；Then, rounding is carried out to former continuous variable；
Wherein, rounding is carried out to former continuous variable with following self-defined transfer functions, the mathematical modeling of the self-defined transfer function
For：
X t =min{X s |X s ≥X i }
Wherein, X t Represent the discrete heat sources value after rounding, X s Represent the standard value in standard alternate pond, X i Represent random production
Raw continuous variable value；
When carrying out rounding to former continuous variable with self-defined transfer function, the continuous variable value and standard alternate that will randomly generate
Data are compared one by one in pond, until approaching standard value that it is not more than and immediate, are taken the standard value to be used as and are connected at random
Centrifugal pump corresponding to continuous variable；
（2）The control of discrete variable design space
The design space combined using marginal check module to discrete variable is controlled, and for crossing the border, parameter is rejected；
The screening rule of the marginal check module is as follows：
If X t > X up Or X t < X low
Then X t = X new
Wherein, X t Represent discrete variable combination to be tested, X up It is the arc radius upper bound artificially set, X low It is artificially to set
Arc radius lower bound, X new For a variate-value in newly-generated standard design space；
The screening module is by the way that the discrete variable regenerated value compared with the radius variable bound pre-set, is come
Judgment variable value is used in one randomly selected in standard design space whether in rational design space, for getting over bound variable
New discrete variable value substitutes.
2. the more circular-arc abnormal-shape hole Equalization Design methods of two variables as claimed in claim 1 based on engineering discrete variable, it is special
Sign is：In step（2）Also include step afterwards（3）Optimized Iterative process, it is specially：Before optimization, discrete design point is by size
Serial number；In Optimized Iterative, select new design point to number by numbering computing, then pass through the numbering of discrete design point
New design point is determined with the mapping relations of actual design point, by the comparison of design point, picks out more preferable new design point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510427969.1A CN105046000B (en) | 2015-07-20 | 2015-07-20 | The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510427969.1A CN105046000B (en) | 2015-07-20 | 2015-07-20 | The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105046000A CN105046000A (en) | 2015-11-11 |
| CN105046000B true CN105046000B (en) | 2017-12-08 |
Family
ID=54452541
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510427969.1A Expired - Fee Related CN105046000B (en) | 2015-07-20 | 2015-07-20 | The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105046000B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106156409B (en) * | 2016-06-24 | 2020-01-14 | 南京航空航天大学 | Two-parameter special-shaped bolt hole robust design method based on minimum contour change |
| CN106971023B (en) * | 2017-03-01 | 2020-05-19 | 南京航空航天大学 | Wheel disc special-shaped hole structure design method based on hyperelliptic curve |
| 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 (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2639189Y (en) * | 2003-09-03 | 2004-09-08 | 瑞安市金泰包装机械厂 | Punching device for special shape hole of continuous roll-cutting type soft packing bag |
| CN101592176A (en) * | 2008-05-28 | 2009-12-02 | 宽仕工业股份有限公司 | Hexagonal irregularly-shaped hole screw |
| JP2010174956A (en) * | 2009-01-28 | 2010-08-12 | Nichirin Co Ltd | Intermediately holding fitting of hose |
| CN103334798A (en) * | 2013-07-10 | 2013-10-02 | 南京航空航天大学 | Multi-arc specially-shaped hole |
-
2015
- 2015-07-20 CN CN201510427969.1A patent/CN105046000B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2639189Y (en) * | 2003-09-03 | 2004-09-08 | 瑞安市金泰包装机械厂 | Punching device for special shape hole of continuous roll-cutting type soft packing bag |
| CN101592176A (en) * | 2008-05-28 | 2009-12-02 | 宽仕工业股份有限公司 | Hexagonal irregularly-shaped hole screw |
| JP2010174956A (en) * | 2009-01-28 | 2010-08-12 | Nichirin Co Ltd | Intermediately holding fitting of hose |
| CN103334798A (en) * | 2013-07-10 | 2013-10-02 | 南京航空航天大学 | Multi-arc specially-shaped hole |
Non-Patent Citations (4)
| Title |
|---|
| Stress analysis of single-bolt,single-lap,countersunk composite joints with variable bolt-hole clearance;B.Egan et al.;《Composite Structures》;20120228;第94卷(第3期);全文 * |
| 涡轮盘双轴对称异形孔结构建模与优化;陈秋任;《航空动力学报》;20130630;第28卷(第6期);全文 * |
| 涡轮盘异型孔结构设计与优化;陈秋任;《中国优秀硕士学位论文全文数据库工程科技II辑》;20150215(第2期);全文 * |
| 涡轮连接单元异型孔应力场分析及结构优化;张超;《中国优秀硕士学位论文全文数据库工程科技II辑》;20150115(第1期);论文第40-43页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105046000A (en) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105046000B (en) | The more circular-arc abnormal-shape hole Equalization Design methods of two variables and model based on engineering discrete variable | |
| CN107103175B (en) | A kind of wind power generating set disorder data recognition method and device | |
| GB2584060A8 (en) | Multi-core sensor system, and isolation and recovery method therefor | |
| EP3905069A3 (en) | Resource recommendation method and apparatus, parameter determination method and apparatus, device, and medium | |
| CN110685651B (en) | Yield splitting method and system for multilayer commingled production gas well | |
| CN105373852A (en) | BIM-based personnel evacuation prompting system and method | |
| WO2020239910A3 (en) | An intelligent computer aided decision support system | |
| CN102305595A (en) | Method for automatically detecting spring pitch by using mechanical vision system | |
| CN106611021B (en) | Data processing method and equipment | |
| CN104346442B (en) | A kind of Rules extraction method of Process-Oriented object data | |
| CN103324888B (en) | Based on virus characteristic extraction method and the system of family's sample | |
| CN104037787B (en) | Power system stabilizer coordination optimization method for improving sending capacity of multiple small hydropower stations | |
| CN106056168B (en) | The determination method of gas-steam combined circulating generation unit operating condition optimal value | |
| CN104268640A (en) | Pilot designating method | |
| CN115034422A (en) | Wind power short-term power prediction method and system based on fluctuation identification and error correction | |
| CN111082473A (en) | Station level reactive power coordination control method and system | |
| CN107886174B (en) | Generator set maintenance arrangement method and device | |
| CN115936232A (en) | Method and system for predicting density and encryption adjustment time of reasonable well pattern of tight gas reservoir | |
| CN105335467A (en) | Novelty judgment method and system for hot topics in publishing industry | |
| CN106156409B (en) | Two-parameter special-shaped bolt hole robust design method based on minimum contour change | |
| CN105183013A (en) | Multiple steam pipeline systems and steam balance method thereof | |
| CN105138753B (en) | The determination method of high Pressure Gas Pipeline node actual flow under accident conditions | |
| CN110957742B (en) | Electric power system small-interference stable auxiliary decision-making method and system based on deep learning | |
| Yong et al. | The simulation of 3D plant leaves wilting | |
| CN102747212A (en) | System and method for optimizing heat buckling in strip steel furnace |
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: 20171208 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |