CN104573199A - Topological optimization method for constraining and damping blades of aeroengine - Google Patents

Topological optimization method for constraining and damping blades of aeroengine Download PDF

Info

Publication number
CN104573199A
CN104573199A CN201410805697.XA CN201410805697A CN104573199A CN 104573199 A CN104573199 A CN 104573199A CN 201410805697 A CN201410805697 A CN 201410805697A CN 104573199 A CN104573199 A CN 104573199A
Authority
CN
China
Prior art keywords
damping
constraining
blade
born
blades
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.)
Pending
Application number
CN201410805697.XA
Other languages
Chinese (zh)
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.)
Chongqing University
AVIC Shenyang Engine Design and Research Institute
Original Assignee
Chongqing University
AVIC Shenyang Engine Design and Research Institute
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 Chongqing University, AVIC Shenyang Engine Design and Research Institute filed Critical Chongqing University
Priority to CN201410805697.XA priority Critical patent/CN104573199A/en
Publication of CN104573199A publication Critical patent/CN104573199A/en
Pending legal-status Critical Current

Links

Landscapes

  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

The invention relates to a topological optimization method for constraining and damping blades of an aeroengine. The method comprises the following steps: taking lightweight aluminum foil as a constraining material, sequentially sticking a viscoelastic material and the constraining material onto the surface of each blade of the aeroengine, carrying out constraining and damping processing on the blades, and establishing finite element models of the constrained and damped blades; carrying out modal analysis on the constrained and damped blades according to the finite element models of the constrained and damped blades, and filtering sensitivity of constrained and damped cell units by adopting an independent grid filtering technology; solving an optimization model by adopting an asymptotic optimization algorithm, programming an optimization program, searching for the optimal layout of a constraining and damping material on each blade, and cutting and regularizing optimized layout, so as to adapt to engineering application. The method has the advantages that the aim of achieving the maximal aeroengine blade damping vibration reduction effect with the minimal damping material additionally-added amount is realized, so that the vibration stress of the blades of the aeroengine is remarkably lowered, and the fatigue damage resistance of the blades is improved.

Description

A kind of Topology Optimization Method of damping-constraining blade of aviation engine
Technical field
The present invention relates to constrained damping structure Optimization Design, particularly a kind of Topology Optimization Method of damping-constraining blade of aviation engine.
Background technology
The structural intergrity of blade of aviation engine is the key of aeromotor security of operation, and along with the increase of aeromotor thrust-weight ratio, the vibration stress of blade increases, and Fatigue Damage Problem is outstanding day by day.At present, many forms such as employing damping shroud or dry-friction damping etc., increase blade damping, reduce blade vibration stress.Constrained damping structure, as a kind of version of effective increase damping, has the simple and efficient feature of process.
The general all standing of traditional constraints damping, in the surface of pending structure, although effectively can increase structural damping, considerably increases additional mass, is unfavorable for the lightweight of aeromotor.
Summary of the invention
The object of the invention is is the optimization of material collocation method proposing the process of a kind of blade of aviation engine damping-constraining, for the target order optimized and constraint damping material use amount, sets up the Topological optimization model of aeromotor damping-constraining blade.Adopt gradual optimization algorithm, the damping-constraining material layout on blade is optimized.By optimizing, with minimum damping-constraining material additional mass, significantly reducing the vibration stress of blade of aviation engine, reaching best damping vibration attenuation effect.
The invention provides a kind of Topology Optimization Method of damping-constraining blade of aviation engine, it is characterized in that: the Topology Optimization Method of described damping-constraining blade of aviation engine, adopt lightweight aluminium foil as constraint material, viscoelastic material and constraint material is pasted successively on the surface at blade of aviation engine, damping-constraining process is carried out to blade, sets up the finite element model of damping-constraining blade;
Determine target order and the constraint damping material use amount of optimization, set up the Topological optimization model of damping-constraining blade:
find : x i , i = 1,2 . . . n min : η r , r = 1,2,3 . s . t : Σ i = 1 n x i M i - M * ≤ 0 0 ≤ x i ≤ 1 , i = 1 , . . . n
Wherein, x ifor design variable, η rfor the modal loss factor on r rank, M ifor the quality of damping-constraining born of the same parents unit, M *for the use amount of damping-constraining material;
According to the finite element model of damping-constraining blade, model analysis is carried out to damping-constraining blade, adopt Sensitivity Analysis Method, calculate the modal loss factor sensitivity of each damping-constraining born of the same parents unit, adopt the sensitivity of independent grid filtering technique to damping-constraining born of the same parents unit to carry out filtering;
Adopt gradual optimization algorithm, solving-optimizing model, establishment optimizer, searches the optimal location of damping-constraining material on blade, and carries out cutting and regular to optimizing distribution, with applicable engineer applied.
The finite element model of the damping-constraining blade of described foundation, adopts the hybridised units of four nodes.
Described design variable x ifor damping-constraining born of the same parents unit existence whether, its value is that { 0,1} works as x iwhen=0, represent that damping-constraining born of the same parents unit i is deleted; Work as x iwhen=1, represent damping-constraining born of the same parents unit i place, be coated with damping-constraining material, the initial value of design variable is 1.
The modal loss factor Calculation of Sensitivity formula of described damping-constraining born of the same parents unit is:
According to modal strain energy method, the r rank modal loss factor of constrained damping structure is calculated as follows:
η r = η v U vr U r
Wherein, η vfor the dissipation factor of damping material; U vrand U rbe respectively damping material strain energy and the constrained damping structure total strain energy of r rank mode;
Make η rto design variable x iask local derviation, r rank modal loss factor is to design variable x isensitivity expression formula be:
∂ ( η r ) ∂ x i = η v U vri e U r - ( U cri e + U vri e ) U vr U r 2
In formula, with be respectively the modal strain energy of the modal strain energy of the damping unit i of r rank mode and restraint layer unit corresponding to damping unit i.
The present invention adopts topological optimization technology, for the target order optimized and damping-constraining materials'use amount, damping-constraining material layout on blade of aviation engine is optimized, makes that the additional mass of blade is minimum, modal loss factor is maximum, vibration stress is minimum, damping vibration attenuation best results.
Advantage of the present invention:
The Topology Optimization Method of damping-constraining blade of aviation engine of the present invention, adopt topological optimization technology, for the target order optimized and constraint damping material use amount, target is turned to so that modal loss factor is maximum, set up blade of aviation engine damping-constraining material layout Optimized model, adopt gradual optimization algorithm, establishment Optimization Calculation Program, obtain the optimal location of damping-constraining material on blade, achieve with minimum damping material additional mass, reach the target of maximum blade of aviation engine damping vibration attenuation effect, the vibration stress of blade of aviation engine is significantly declined, improve the ability of blade antifatigue damage.
Accompanying drawing explanation
Below in conjunction with drawings and the embodiments, the present invention is further detailed explanation:
Fig. 1 is the damping-constraining blade after optimizing;
Fig. 2 is three-dimensional spectrum, excitation air pressure 0.20MPa, blade;
Fig. 3 is three-dimensional spectrum, the damping-constraining blade after excitation air pressure 0.20MPa, optimization.
Embodiment
Embodiment 1
Present embodiments provide a kind of Topology Optimization Method of damping-constraining blade of aviation engine, it is characterized in that: the Topology Optimization Method of described damping-constraining blade of aviation engine, adopt lightweight aluminium foil as constraint material, viscoelastic material and constraint material is pasted successively on the surface at blade of aviation engine, damping-constraining process is carried out to blade, sets up the finite element model of damping-constraining blade;
Determine target order and the constraint damping material use amount of optimization, set up the Topological optimization model of damping-constraining blade:
find : x i , i = 1,2 . . . n min : η r , r = 1,2,3 . s . t : Σ i = 1 n x i M i - M * ≤ 0 0 ≤ x i ≤ 1 , i = 1 , . . . n
Wherein, x ifor design variable, η rfor the modal loss factor on r rank, M ifor the quality of damping-constraining born of the same parents unit, M *for the use amount of damping-constraining material;
According to the finite element model of damping-constraining blade, model analysis is carried out to damping-constraining blade, adopt Sensitivity Analysis Method, calculate the modal loss factor sensitivity of each damping-constraining born of the same parents unit, adopt the sensitivity of independent grid filtering technique to damping-constraining born of the same parents unit to carry out filtering;
Adopt gradual optimization algorithm, solving-optimizing model, establishment optimizer, searches the optimal location of damping-constraining material on blade, and carries out cutting and regular to optimizing distribution, with applicable engineer applied.
The finite element model of the damping-constraining blade of described foundation, adopts the hybridised units of four nodes.
Described design variable x ifor damping-constraining born of the same parents unit existence whether, its value is that { 0,1} works as x iwhen=0, represent that damping-constraining born of the same parents unit i is deleted; Work as x iwhen=1, represent damping-constraining born of the same parents unit i place, be coated with damping-constraining material, the initial value of design variable is 1.
The modal loss factor Calculation of Sensitivity formula of described damping-constraining born of the same parents unit is:
According to modal strain energy method, the r rank modal loss factor of constrained damping structure is calculated as follows:
η r = η v U vr U r
Wherein, η vfor the dissipation factor of damping material; U vrand U rbe respectively damping material strain energy and the constrained damping structure total strain energy of r rank mode;
Make η rto design variable x iask local derviation, r rank modal loss factor is to design variable x isensitivity expression formula be:
∂ ( η r ) ∂ x i = η v U vri e U r - ( U cri e + U vri e ) U vr U r 2
In formula, with be respectively the modal strain energy of the modal strain energy of the damping unit i of r rank mode and restraint layer unit corresponding to damping unit i.
The present invention adopts topological optimization technology, for the target order optimized and damping-constraining materials'use amount, damping-constraining material layout on blade of aviation engine is optimized, makes that the additional mass of blade is minimum, modal loss factor is maximum, vibration stress is minimum, damping vibration attenuation best results.

Claims (4)

1. the Topology Optimization Method of a damping-constraining blade of aviation engine, it is characterized in that: the Topology Optimization Method of described damping-constraining blade of aviation engine, adopt lightweight aluminium foil as constraint material, viscoelastic material and constraint material is pasted successively on the surface at blade of aviation engine, damping-constraining process is carried out to blade, sets up the finite element model of damping-constraining blade;
Determine target order and the constraint damping material use amount of optimization, set up the Topological optimization model of damping-constraining blade:
find : x i , i = 1,2 . . . n min : η r , r = 1,2,3 . s . t : Σ i = 1 n x i M i - M * ≤ 0 0 ≤ x i ≤ 1 , i = 1 , . . . n
Wherein, x ifor design variable, η rfor the modal loss factor on r rank, M ifor the quality of damping-constraining born of the same parents unit, M *for the use amount of damping-constraining material;
According to the finite element model of damping-constraining blade, model analysis is carried out to damping-constraining blade, adopt Sensitivity Analysis Method, calculate the modal loss factor sensitivity of each damping-constraining born of the same parents unit, adopt the sensitivity of independent grid filtering technique to damping-constraining born of the same parents unit to carry out filtering;
Adopt gradual optimization algorithm, solving-optimizing model, establishment optimizer, searches the optimal location of damping-constraining material on blade, and carries out cutting and regular to optimizing distribution, with applicable engineer applied.
2. according to the Topology Optimization Method of damping-constraining blade of aviation engine according to claim 1, it is characterized in that: the finite element model of the damping-constraining blade of described foundation, adopts the hybridised units of four nodes.
3. according to the Topology Optimization Method of damping-constraining blade of aviation engine according to claim 1, it is characterized in that: described design variable x ifor damping-constraining born of the same parents unit existence whether, its value is that { 0,1} works as x iwhen=0, represent that damping-constraining born of the same parents unit i is deleted; Work as x iwhen=1, represent damping-constraining born of the same parents unit i place, be coated with damping-constraining material, the initial value of design variable is 1.
4. according to the Topology Optimization Method of damping-constraining blade of aviation engine according to claim 1, it is characterized in that: the modal loss factor Calculation of Sensitivity formula of described damping-constraining born of the same parents unit is:
According to modal strain energy method, the r rank modal loss factor of constrained damping structure is calculated as follows:
η r = η v U vr U r
Wherein, η vfor the dissipation factor of damping material; U vrand U rbe respectively damping material strain energy and the constrained damping structure total strain energy of r rank mode;
Make η rto design variable x iask local derviation, r rank modal loss factor is to design variable x isensitivity expression formula be:
∂ ( η r ) ∂ x i = η v U vri e U r - ( U cri e + U vri e ) U vr U r 2
In formula, with be respectively the modal strain energy of the modal strain energy of the damping unit i of r rank mode and restraint layer unit corresponding to damping unit i.
CN201410805697.XA 2014-12-19 2014-12-19 Topological optimization method for constraining and damping blades of aeroengine Pending CN104573199A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410805697.XA CN104573199A (en) 2014-12-19 2014-12-19 Topological optimization method for constraining and damping blades of aeroengine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410805697.XA CN104573199A (en) 2014-12-19 2014-12-19 Topological optimization method for constraining and damping blades of aeroengine

Publications (1)

Publication Number Publication Date
CN104573199A true CN104573199A (en) 2015-04-29

Family

ID=53089253

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410805697.XA Pending CN104573199A (en) 2014-12-19 2014-12-19 Topological optimization method for constraining and damping blades of aeroengine

Country Status (1)

Country Link
CN (1) CN104573199A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017215217A1 (en) * 2016-06-16 2017-12-21 华南理工大学 Topology optimization design method for flexible hinge
CN110688795A (en) * 2019-09-25 2020-01-14 国网湖南省电力有限公司 Transformer box damping vibration attenuation method, system and medium based on topology optimization
CN111382536A (en) * 2020-03-14 2020-07-07 郑州轻工业大学 Collaborative optimization design method of constrained damping structure
CN112287477A (en) * 2020-10-16 2021-01-29 中国航发四川燃气涡轮研究院 Turbine guide device large and small blade layout method based on airflow excitation
CN112765732A (en) * 2021-01-25 2021-05-07 沈阳工业大学 Aviation blade topology optimization design method based on selective laser melting process
CN116541970A (en) * 2023-06-29 2023-08-04 中国航发四川燃气涡轮研究院 Energy method-based compressor blade vibration reduction optimization design evaluation method
CN116976121A (en) * 2023-08-02 2023-10-31 中国长江三峡集团有限公司 Constrained damping wind power tower and parameter optimization method thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103294858A (en) * 2013-05-24 2013-09-11 重庆大学 Constrained damping cylindrical shell topological optimization method based on modal loss factor maximization
CN104182598A (en) * 2014-09-18 2014-12-03 重庆大学 Constraint damping structure optimizing and designing method based on level set method
CN104200048A (en) * 2014-09-18 2014-12-10 重庆大学 Constrained damping plate optimization design method based on bidirectional evolutionary optimization algorithm

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103294858A (en) * 2013-05-24 2013-09-11 重庆大学 Constrained damping cylindrical shell topological optimization method based on modal loss factor maximization
CN104182598A (en) * 2014-09-18 2014-12-03 重庆大学 Constraint damping structure optimizing and designing method based on level set method
CN104200048A (en) * 2014-09-18 2014-12-10 重庆大学 Constrained damping plate optimization design method based on bidirectional evolutionary optimization algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BRET STANFORD.ET AL.: ""Aeroelastic Topology Optimization of Blade-Stiffened Panels"", 《JOURNAL OF AIRCRAFT》 *
刘利军: ""转动ACLD悬臂板的动力学建模及振动控制"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *
房占鹏: ""约束阻尼结构的双向渐进拓扑优化"", 《振动与冲击》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017215217A1 (en) * 2016-06-16 2017-12-21 华南理工大学 Topology optimization design method for flexible hinge
CN110688795A (en) * 2019-09-25 2020-01-14 国网湖南省电力有限公司 Transformer box damping vibration attenuation method, system and medium based on topology optimization
CN111382536A (en) * 2020-03-14 2020-07-07 郑州轻工业大学 Collaborative optimization design method of constrained damping structure
CN111382536B (en) * 2020-03-14 2023-06-09 郑州轻工业大学 Collaborative optimization design method for constraint damping structure
CN112287477A (en) * 2020-10-16 2021-01-29 中国航发四川燃气涡轮研究院 Turbine guide device large and small blade layout method based on airflow excitation
CN112287477B (en) * 2020-10-16 2022-11-22 中国航发四川燃气涡轮研究院 Turbine guide device large and small blade layout method based on airflow excitation
CN112765732A (en) * 2021-01-25 2021-05-07 沈阳工业大学 Aviation blade topology optimization design method based on selective laser melting process
CN112765732B (en) * 2021-01-25 2024-03-01 沈阳工业大学 Aviation blade topology optimization design method based on selective laser melting process
CN116541970A (en) * 2023-06-29 2023-08-04 中国航发四川燃气涡轮研究院 Energy method-based compressor blade vibration reduction optimization design evaluation method
CN116541970B (en) * 2023-06-29 2023-10-03 中国航发四川燃气涡轮研究院 Energy method-based compressor blade vibration reduction optimization design evaluation method
CN116976121A (en) * 2023-08-02 2023-10-31 中国长江三峡集团有限公司 Constrained damping wind power tower and parameter optimization method thereof
CN116976121B (en) * 2023-08-02 2024-05-14 中国长江三峡集团有限公司 Constrained damping wind power tower and parameter optimization method thereof

Similar Documents

Publication Publication Date Title
CN104573199A (en) Topological optimization method for constraining and damping blades of aeroengine
Bak et al. The DTU 10-MW reference wind turbine
Potsdam et al. Unstructured mesh CFD aerodynamic analysis of the NREL Phase VI rotor
Chen et al. Structural optimization study of composite wind turbine blade
Song et al. Optimization design, modeling and dynamic analysis for composite wind turbine blade
Nadarajah Adjoint-based aerodynamic optimization of benchmark problems
Lee et al. Performance improvement of horizontal axis wind turbines by aerodynamic shape optimization including aeroealstic deformation
CN105718621B (en) Engine external support optimization design method
CN103294858A (en) Constrained damping cylindrical shell topological optimization method based on modal loss factor maximization
CN104268326A (en) Optimization criterion method based constraint damping plate topological optimization method
CN103810341A (en) Predicating method for wind turbine blade airfoil fluttering
Luo et al. A study on multidisciplinary optimization of an axial compressor blade based on evolutionary algorithms
CN104408243B (en) A kind of engine mockup installed thrust modification method
CN104200047B (en) Constrained damping drum optimization design method for optimally minimizing square of resonance peak values
Yu et al. Numerical simulation of oscillating-wing based energy harvest mechanism using the high-order spectral difference method
CN103823787A (en) Multi-turning-tool parallel turning stability judgment method based on differential quadrature method
CN104182598A (en) Constraint damping structure optimizing and designing method based on level set method
Song et al. Blade shape optimization of transonic axial flow fan in terms of sectional profiles and stacking line
CN106741924A (en) A kind of high lift-drag ratio Low-torque characteristic 9% thickness rotor-blade airfoil under the conditions of full working scope
刘明周 et al. Dynamic rescheduling optimization of job-shop under uncertain conditions
Razaaly et al. Uncertainty Quantification of an ORC turbine blade under a low quantile constrain
Young et al. Optimization of the skin thickness distribution in the composite wind turbine blade
Iyengar et al. A first-principles based methodology for design of axial compressor Configurations
Chen et al. Nacelle: air intake aerodynamic design and inlet compatibility
Kha et al. A FEA Analysis of Blade of HAWT: A Review

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20150429