CN104200003A - Self-locking damping vane gas flow bending stress design method based on fluid-solid interaction - Google Patents
Self-locking damping vane gas flow bending stress design method based on fluid-solid interaction Download PDFInfo
- Publication number
- CN104200003A CN104200003A CN201410358976.6A CN201410358976A CN104200003A CN 104200003 A CN104200003 A CN 104200003A CN 201410358976 A CN201410358976 A CN 201410358976A CN 104200003 A CN104200003 A CN 104200003A
- Authority
- CN
- China
- Prior art keywords
- blade
- bending stress
- fluid
- solid
- vane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
The invention discloses a self-locking damping vane gas flow bending stress design method based on fluid-solid interaction. The method is characterized in that a one-to-one iteration method for node coordinates is designed with a self-locking damping vane and surrounding fluid serving as a three-dimensional coupling system, and the gas flow bending stress of the self-locking damping vane is obtained by solving the deformation and the fluid field of the vane through coupling. The method includes the steps of firstly, establishing a fluid calculation area through Turbogrid software, and inputting a grid into the CFX to calculate the fluid field distribution around the vane; secondly, establishing a solid calculation area through solidworks software, mapping the vane surface pressure calculated through the CFX to a solid calculation grid through the node coordinates to be calculated, obtaining the gas flow bending stress through calculation in the ANSYS, and judging whether the requirement is met through safety magnification calibration and verification. By means of the self-locking damping vane gas flow bending stress design method based on fluid-solid interaction, the gas flow bending stress distribution condition can be more truly reflected, vane design optimization is facilitated, vane test expenses are reduced, and the method is suitable for wide-range popularization.
Description
Technical field
The present invention relates to turbomachine field, especially a kind of interlocked blade gas-flow bending stress method for designing based on the solid coupling of stream.
Background technology
At present, the design of large flow variable rotating speed industrial steam turbine low pressure stage group exhaust stage blade is that industry drags one of the difficult point in field always.Meeting under high efficiency designing requirement, guarantee the safe operation of unit, will carry out Comprehensive Designing to the Safety and Service Life of blade.Because last stage movable vane newly developed is frequency-modulating blade, its blade is thinner, the factor of paying close attention to when gas-flow bending stress is design.At present, the conventional line integral method that most of research blade gas-flow bending stress is used, do not consider the impact of coupled system, also do not consider the impact of centrifugal prestress and blade boss chamfering, blade root chamfering, therefore when the distribution situation of reaction blade gas-flow bending stress, data are not accurate enough true.
Summary of the invention
The object of the invention is to provide in order to solve the deficiency of above-mentioned technology the interlocked blade gas-flow bending stress method for designing based on the solid coupling of stream of the distribution situation of the accurate actual response blade of a kind of energy gas-flow bending stress.
In order to achieve the above object, a kind of interlocked blade gas-flow bending stress method for designing based on the solid coupling of stream that the present invention is designed, it is characterized in that: these computing method are: (1) is using interlocked blade and surrounding fluid as a three-dimensional coupled system, design a set of nodes oriented coordinate alternative manner one to one, and solve by coupling the gas-flow bending stress that deformable blade and flow field obtain blade; (2) according to the physical size requirement of interlocked blade, adopt Turbogrid software to set up fluid calculation region, grid is input to CFX and calculates this blade Flow Field Distribution around; (3) adopt Solidworks software to set up solid zoning, the blade surface pressure that CFX is calculated is mapped on solid computing grid one by one by node coordinate, flow solid coupling and calculate, and in ANSYS, calculate and be with prestressed gas-flow bending stress; (4) by safe multiplying power, check criterion, judge whether gas-flow bending stress distribution meets design requirement, and then adjusts vane foil parameter, and cycle calculations obtains the interlocked blade that will design.
As preferably, a kind of based on the solid interlocked blade gas-flow bending stress method for designing being coupled of stream, it is characterized in that: this calculating detailed step comprises the following steps:
Step 1, the foundation of how much physical models; According to the interlocked blade parameter of primary design, first adopt Turbogrid software to set up fluid domain, then adopt Solidworks software to set up solid domain, and blade shroud band and boss are rotated periodically and are cut apart.
Step 2, grid is divided; Convection cell territory adopts Turbogrid to carry out hexahedron structure grid and divides, and to solid domain, adopts Ansysmeshing to carry out hexahedron tetrahedron mixed mesh.
Step 3, Flow Field Calculation; Fluid grid is imported to CFX, and arrange with downstream condition: moving territory adopts multiple system method to process, and blade adopts periodic boundary, the given stagnation pressure of entrance, stagnation temperature, mass dryness fraction, outlet is to constant static-pressure, and fluid adopts IF97 water vapor model, and turbulence model is set
, fluid governing equation adopts Finite Volume Method to solve, and selects second order spatial spreading form.
Step 4, solid static stress is calculated; Solid grid is imported to ANSYS, adopts following Cyclic Symmetry periodic boundary condition: fixedly blade root surface of contact, given rotating speed load, and then it is discrete to adopt implicit expression method of dynamic analysis, governing equation to solve employing Finite Element Method, the distribution of calculating centrifugal force.
Step 5, the solid Coupled Numerical of stream calculates; At step 4 blade, bear on centrifugal prestressed basis, the blade surface pressure of fluid calculation is mapped on the blade surface of solid domain one by one by node coordinate, adopt iteration coupling analytical method, governing equation adopts complete Newton iteration method, and the blade stress calculating with gas-flow bending stress distributes.
Step 6, gas-flow bending stress aftertreatment; The blade stress that step 4 and step 5 are calculated distributes and outputs to Excel and do difference computing, and its result is blade and distributes at the gas-flow bending stress bearing under pre-stressed state.
Step 7, safe multiplying power is checked; According to blade material composite fatigue intensity curve, find resistance to the shake intensity of blade under mean stress, and the computing formula of safe multiplying power is
, according to the performance of vibration strength of blade safety criterion judge blade; If do not met, remodify profile parameters of vane, return to the first step, loop iteration, until obtain the interlocked blade meeting design requirement
The resulting a kind of interlocked blade gas-flow bending stress method for designing based on the solid coupling of stream of the present invention, can react more really the distribution situation of blade gas-flow bending stress compared with prior art, and these computing method mainly adopt is to take the fatigue resistance criterion that the static strength criterion that permissible stress is benchmark and the safe multiplying power of take be benchmark, like this can the blade design stage just obtain blade various factors in service for the impact of leaf longevity loss and blade on the distribution of various piece life consumption, thereby find out the weak link on blade, contribute to the testing expenses of optimizing blade design and reducing the blade design stage, the effect size that can also quantitative test affects the various factors of blade safe operation, thereby not only for Turbine Blade Design provides a new appraisal procedure of evaluating blade security from life consumption aspect, and new approaches and utility are provided for carrying out leaf longevity management aspect the steam turbine operation mode day by day complicated, therefore be suitable for promoting on a large scale and using.
Accompanying drawing explanation
Fig. 1 is Technology Roadmap of the present invention.
Fig. 2 is the present invention's certain interlocked blade geometric model used.
Fig. 3 is fluid calculation of the present invention territory geometric model and motion pattern.
Fig. 4 is interlocked blade centrifugal force of the present invention and the gas-flow bending stress distribution plan of making a concerted effort.
Fig. 5 is the interlocked blade gas-flow bending stress of the present invention cloud charts of making a concerted effort.
Fig. 6 is that interlocked blade of the present invention is along the high cross section of leaf maximum airflow bending stress distribution plan.
Fig. 7 is the safe multiplying power check of interlocked blade of the present invention figure.
Embodiment
Below by embodiment, the invention will be further described by reference to the accompanying drawings.
Embodiment 1.
As shown in Fig. 1-7, a kind of interlocked blade gas-flow bending stress method for designing based on the solid coupling of stream that the present embodiment is described, it is characterized in that: these computing method are: (1) is using interlocked blade and surrounding fluid as a three-dimensional coupled system, design a set of nodes oriented coordinate alternative manner one to one, and solve by coupling the gas-flow bending stress that deformable blade and flow field obtain blade; (2) according to the physical size requirement of interlocked blade, adopt Turbogrid software to set up fluid calculation region, grid is input to CFX and calculates this blade Flow Field Distribution around; (3) adopt Solidworks software to set up solid zoning, the blade surface pressure that CFX is calculated is mapped on solid computing grid one by one by node coordinate, flow solid coupling and calculate, and in ANSYS, calculate and be with prestressed gas-flow bending stress; (4) by safe multiplying power, check criterion, judge whether gas-flow bending stress distribution meets design requirement, and then adjusts vane foil parameter, and cycle calculations obtains the interlocked blade that will design, and this calculating detailed step comprises the following steps:
Step 1, the foundation of how much physical models; According to the interlocked blade parameter of primary design, first adopt Turbogrid software to set up fluid domain, then adopt Solidworks software to set up solid domain, and blade shroud band and boss are rotated periodically and are cut apart;
Step 2, grid is divided; Convection cell territory adopts Turbogrid to carry out hexahedron structure grid and divides, and to solid domain, adopts Ansysmeshing to carry out hexahedron tetrahedron mixed mesh;
Step 3, Flow Field Calculation; Fluid grid is imported to CFX, and arrange with downstream condition: moving territory adopts multiple system method to process, and blade adopts periodic boundary, the given stagnation pressure of entrance, stagnation temperature, mass dryness fraction, outlet is to constant static-pressure, and fluid adopts IF97 water vapor model, and turbulence model is set
, fluid governing equation adopts Finite Volume Method to solve, and selects second order spatial spreading form;
Step 4, solid static stress is calculated; Solid grid is imported to ANSYS, adopts following Cyclic Symmetry periodic boundary condition: fixedly blade root surface of contact, given rotating speed load, and then it is discrete to adopt implicit expression method of dynamic analysis, governing equation to solve employing Finite Element Method, the distribution of calculating centrifugal force;
Step 5, the solid Coupled Numerical of stream calculates; At step 4 blade, bear on centrifugal prestressed basis, the blade surface pressure of fluid calculation is mapped on the blade surface of solid domain one by one by node coordinate, adopt iteration coupling analytical method, governing equation adopts complete Newton iteration method, and the blade stress calculating with gas-flow bending stress distributes;
Step 6, gas-flow bending stress aftertreatment; The blade stress that step 4 and step 5 are calculated distributes and outputs to Excel and do difference computing, and its result is blade and distributes at the gas-flow bending stress bearing under pre-stressed state;
Step 7, safe multiplying power is checked; According to blade material composite fatigue intensity curve, find resistance to the shake intensity of blade under mean stress, and the computing formula of safe multiplying power is
wherein k is safety coefficient, passes judgment on the performance of blade according to vibration strength of blade safety criterion; If do not met, remodify profile parameters of vane, return to the first step, loop iteration, until obtain the interlocked blade meeting design requirement.
Claims (2)
1. based on stream, consolidate the interlocked blade gas-flow bending stress method for designing being coupled for one kind, it is characterized in that: (1) is using interlocked blade and surrounding fluid as a three-dimensional coupled system, design a set of nodes oriented coordinate alternative manner one to one, and solve by coupling the gas-flow bending stress that deformable blade and flow field obtain blade; (2) according to the physical size requirement of interlocked blade, adopt Turbogrid software to set up fluid calculation region, grid is input to CFX and calculates this blade Flow Field Distribution around; (3) adopt Solidworks software to set up solid zoning, the blade surface pressure that CFX is calculated is mapped on solid computing grid one by one by node coordinate, flow solid coupling and calculate, and in ANSYS, calculate and be with prestressed gas-flow bending stress; (4) by safe multiplying power, check criterion, judge whether gas-flow bending stress distribution meets design requirement, and then adjusts vane foil parameter, and cycle calculations obtains the interlocked blade that will design.
2. according to claim 1 a kind of based on the solid interlocked blade gas-flow bending stress method for designing being coupled of stream, it is characterized in that: this calculating detailed step comprises the following steps:
Step 1, the foundation of how much physical models; According to the interlocked blade parameter of primary design, first adopt Turbogrid software to set up fluid domain, then adopt Solidworks software to set up solid domain, and blade shroud band and boss are rotated periodically and are cut apart;
Step 2, grid is divided; Convection cell territory adopts Turbogrid to carry out hexahedron structure grid and divides, and to solid domain, adopts Ansysmeshing to carry out hexahedron tetrahedron mixed mesh;
Step 3, Flow Field Calculation; Fluid grid is imported to CFX, and arrange with downstream condition: moving territory adopts multiple system method to process, and blade adopts periodic boundary, the given stagnation pressure of entrance, stagnation temperature, mass dryness fraction, outlet is to constant static-pressure, and fluid adopts IF97 water vapor model, and turbulence model is set
, fluid governing equation adopts Finite Volume Method to solve, and selects second order spatial spreading form;
Step 4, solid static stress is calculated; Solid grid is imported to ANSYS, adopts following Cyclic Symmetry periodic boundary condition: fixedly blade root surface of contact, given rotating speed load, and then it is discrete to adopt implicit expression method of dynamic analysis, governing equation to solve employing Finite Element Method, the distribution of calculating centrifugal force;
Step 5, the solid Coupled Numerical of stream calculates; At step 4 blade, bear on centrifugal prestressed basis, the blade surface pressure of fluid calculation is mapped on the blade surface of solid domain one by one by node coordinate, adopt iteration coupling analytical method, governing equation adopts complete Newton iteration method, and the blade stress calculating with gas-flow bending stress distributes;
Step 6, gas-flow bending stress aftertreatment; The blade stress that step 4 and step 5 are calculated distributes and outputs to Excel and do difference computing, and its result is blade and distributes at the gas-flow bending stress bearing under pre-stressed state;
Step 7, safe multiplying power is checked; According to blade material composite fatigue intensity curve, find resistance to the shake intensity of blade under mean stress, and the computing formula of safe multiplying power is
, according to the performance of vibration strength of blade safety criterion judge blade; If do not met, remodify profile parameters of vane, return to the first step, loop iteration, until obtain the interlocked blade meeting design requirement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410358976.6A CN104200003B (en) | 2014-07-28 | 2014-07-28 | Interlocked blade gas-flow bending stress design method based on fluid structurecoupling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410358976.6A CN104200003B (en) | 2014-07-28 | 2014-07-28 | Interlocked blade gas-flow bending stress design method based on fluid structurecoupling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104200003A true CN104200003A (en) | 2014-12-10 |
| CN104200003B CN104200003B (en) | 2017-09-19 |
Family
ID=52085296
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410358976.6A Active CN104200003B (en) | 2014-07-28 | 2014-07-28 | Interlocked blade gas-flow bending stress design method based on fluid structurecoupling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104200003B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108829917A (en) * | 2018-04-25 | 2018-11-16 | 内蒙古工业大学 | A kind of stalk throws impeller Prediction method for fatigue life |
| CN109753716A (en) * | 2018-12-28 | 2019-05-14 | 山东大学 | Core based on flow field simulation/thermal power steam turbine group fluid excitation numerical computation method and system |
| CN110020489A (en) * | 2019-04-15 | 2019-07-16 | 哈尔滨汽轮机厂有限责任公司 | The method for determining turbine blade erosion protection sheild geomery is analyzed based on CFD |
| CN110929419A (en) * | 2018-12-29 | 2020-03-27 | 山东大学 | Method for quickly predicting instability limit of steam turbine rotor system based on shroud zero damping |
| CN112861288A (en) * | 2021-03-08 | 2021-05-28 | 浙江水泵总厂有限公司 | Fluid machinery checking method |
| CN116720265A (en) * | 2023-08-09 | 2023-09-08 | 北京航空航天大学 | Unified parameterization method and device for variable-dimension FFD-based turbine |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101158990A (en) * | 2007-11-29 | 2008-04-09 | 北京航空航天大学 | Big and small blade integral leaf disc structural integrity fluid-solid coupling integrated design method |
-
2014
- 2014-07-28 CN CN201410358976.6A patent/CN104200003B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101158990A (en) * | 2007-11-29 | 2008-04-09 | 北京航空航天大学 | Big and small blade integral leaf disc structural integrity fluid-solid coupling integrated design method |
Non-Patent Citations (4)
| Title |
|---|
| FRIEDRICH-KARL BENRA 等: "COMPARISON OF PUMP IMPELLER ORBIT CURVES OBTAINED BY MEASUREMENT AND FSI SIMULATION", 《2007 ASME PRESSURE VESSELS AND PIPING DIVISION CONFERENCE》 * |
| 张晓辉: "考虑流固耦合的压气机某级叶片静强度分析", 《中国优秀硕士学位论文全文数据库-工程科技Ⅱ辑》 * |
| 潘旭 等: "轴流泵叶片流固耦合强度分析", 《水力发电学报》 * |
| 隋永枫 等: "工业汽轮机CAE技术研发与应用", 《CAD/CAM与制造业信息化》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108829917A (en) * | 2018-04-25 | 2018-11-16 | 内蒙古工业大学 | A kind of stalk throws impeller Prediction method for fatigue life |
| CN108829917B (en) * | 2018-04-25 | 2022-04-01 | 内蒙古工业大学 | Fatigue life prediction method for straw throwing impeller |
| CN109753716A (en) * | 2018-12-28 | 2019-05-14 | 山东大学 | Core based on flow field simulation/thermal power steam turbine group fluid excitation numerical computation method and system |
| CN109753716B (en) * | 2018-12-28 | 2020-12-11 | 山东大学 | Nuclear/thermal power turboset fluid excitation numerical calculation method and system based on flow field simulation |
| CN110929419A (en) * | 2018-12-29 | 2020-03-27 | 山东大学 | Method for quickly predicting instability limit of steam turbine rotor system based on shroud zero damping |
| CN110020489A (en) * | 2019-04-15 | 2019-07-16 | 哈尔滨汽轮机厂有限责任公司 | The method for determining turbine blade erosion protection sheild geomery is analyzed based on CFD |
| CN112861288A (en) * | 2021-03-08 | 2021-05-28 | 浙江水泵总厂有限公司 | Fluid machinery checking method |
| CN112861288B (en) * | 2021-03-08 | 2021-11-23 | 浙江水泵总厂有限公司 | Fluid machinery checking method |
| WO2022188394A1 (en) * | 2021-03-08 | 2022-09-15 | 浙江水泵总厂有限公司 | Fluid machinery checking method |
| CN116720265A (en) * | 2023-08-09 | 2023-09-08 | 北京航空航天大学 | Unified parameterization method and device for variable-dimension FFD-based turbine |
| CN116720265B (en) * | 2023-08-09 | 2023-12-08 | 北京航空航天大学 | Unified parameterization method and device for variable-dimension FFD-based turbine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104200003B (en) | 2017-09-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Toja-Silva et al. | A review of computational fluid dynamics (CFD) simulations of the wind flow around buildings for urban wind energy exploitation | |
| Subramanian et al. | Effect of airfoil and solidity on performance of small scale vertical axis wind turbine using three dimensional CFD model | |
| CN104200003A (en) | Self-locking damping vane gas flow bending stress design method based on fluid-solid interaction | |
| Chen et al. | Performance improvement of a vertical axis wind turbine by comprehensive assessment of an airfoil family | |
| Gong et al. | Application of entropy production theory to hydro-turbine hydraulic analysis | |
| Li et al. | Aerodynamic optimization of wind turbine airfoils using response surface techniques | |
| Shehata et al. | Enhancement of performance of wave turbine during stall using passive flow control: First and second law analysis | |
| Chandrala et al. | Aerodynamic analysis of horizontal axis wind turbine blade | |
| Zhou et al. | Influence of inflow conditions on turbine loading and wake structures predicted by large eddy simulations using exact geometry | |
| Rogowski et al. | Scale adaptive simulation model for the darrieus wind turbine | |
| Tabib et al. | A full-scale 3D Vs 2.5 D Vs 2D analysis of flow pattern and forces for an industrial-scale 5MW NREL reference wind-turbine. | |
| Oggiano | CFD simulations on the NTNU wind turbine rotor and comparison with experiments | |
| Kecskemety et al. | Influence of wake effects and inflow turbulence on wind turbine loads | |
| Manca et al. | Numerical investigation of thermal and fluid dynamic behavior of solar chimney building systems | |
| Suffer et al. | Modeling and numerical simulation of a vertical axis wind turbine having cavity vanes | |
| Zhang et al. | Multiobject reliability analysis of turbine blisk with multidiscipline under multiphysical field interaction | |
| Buonomo et al. | Thermal and fluid dynamic analysis of solar chimney building systems | |
| Sun et al. | Numerical simulation of the unsteady flow and power of horizontal axis wind turbine using sliding mesh | |
| Balduzzi et al. | A hybrid BEM-CFD model for effective numerical siting analyses of wind turbines in the urban environment | |
| Chen et al. | Application of numercal simulation to obtain the optimization pitch angle for VAWT | |
| Aksenov et al. | A validation study using nrel phase VI experiments, Part I: Low computational resource scenario | |
| Balasubramanian et al. | Novel curvature-based airfoil parameterization for wind turbine application and optimization | |
| Micallef et al. | Validating BEM, Direct and Inverse Free Wake Models with the MEXICO Experiment. | |
| Ozdemir et al. | An advanced method for wind turbine wake modeling | |
| CN110489887A (en) | Modeling method that a kind of turbine blade based on CFD is through-flow |
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 | ||
| CP03 | Change of name, title or address |
Address after: 310022 Building 1, No. 608, Kangxin Road, Linping District, Hangzhou, Zhejiang Patentee after: Hangzhou Turbine Power Group Co.,Ltd. Address before: No. 357, Shiqiao Road, Jianggan District, Hangzhou, Zhejiang 310022 Patentee before: HANGZHOU STEAM TURBINE Co.,Ltd. |
|
| CP03 | Change of name, title or address |