CN111859729B - 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 - Google Patents
考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 Download PDFInfo
- Publication number
- CN111859729B CN111859729B CN202010500736.0A CN202010500736A CN111859729B CN 111859729 B CN111859729 B CN 111859729B CN 202010500736 A CN202010500736 A CN 202010500736A CN 111859729 B CN111859729 B CN 111859729B
- Authority
- CN
- China
- Prior art keywords
- shot blasting
- blasting process
- stress
- shot
- wheel disc
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/20—Recycling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑多弹丸随机分布的喷丸模型对轮盘寿命的快速计算方法,属于航空发动机轮盘寿命管理领域。首先,本发明基于所需的喷丸工艺参数,判断其是否包含在已建立的数据库中,若是,则从数据库中快速识别提取出残余应力,否则通过二次开发程序建立多弹丸随机分布的喷丸模型并进行有限元分析,获得不同深度处的残余应力。然后,将获得的残余应力结合考虑轮盘工作应力,利用通用权函数法计算得到应力强度因子。最后,利用裂纹扩展Nasgro公式计算得到轮盘寿命。本发明为考虑实际喷丸过程的轮盘寿命评估奠定基础,通过考虑多弹丸随机分布模型规避了少弹丸对称分布模型带来的误差,并建立数据库,实现轮盘寿命高效准确计算。
Description
技术领域
本发明属于航空发动机轮盘的寿命管理领域,涉及一种考虑多弹丸随机分布的喷丸模型对轮盘寿命的快速计算方法。
背景技术
航空发动机的安全性是保证飞行器安全运行的关键,安全性保障需通过安全性分析技术实现;美国航空工业界针对旋转轮盘等关键部件的安全性保障,是通过一套完整的寿命限制件寿命评估流程实现的,其核心即概率失效风险评估技术。该技术在航空发动机轮盘的寿命管理中展示出很大的潜力,因而这些年来获得了广泛的应用。它可以处理载荷、工作环境和材料等参数的随机性对部件安全性的影响,适航规章中规定航空发动机在进行适航取证时需要使用此方法表明限制寿命件在寿命期内的失效概率小于设计目标风险(DTR),发动机才能获得最终的型号合格证。
应力分布及寿命计算是概率风险评估方法的重要输入参数,其中,寿命计算主要受到应力分布的影响。作为典型限制寿命件的轮盘受到的应力不仅包括离心应力、热应力等工作应力,还包括加工过程产生的轮盘内部及表面的残余应力。轮盘类部件在进行喷丸处理过程中,大量高速运动的弹丸不断撞击金属靶件表面,在靶件表面产生不均匀塑性变形,使靶件表面产生残余压应力,此残余压应力场会影响裂纹的裂纹扩展速率,进而会影响工件疲劳寿命和工件表面完整性。
然而,喷丸过程中弹丸的形态是大量随机的,在商用仿真软件中无法实现多弹丸快速生成并随机分布的要求,故目前针对喷丸加工的研究多是基于少弹丸对称分布的模型,但这与实际喷丸过程不符。并且由于仿真模拟计算量庞大,无法根据给定喷丸工艺参数快速得到对应残余应力分布,这就给轮盘寿命计算提出了难题。
发明内容
针对现有技术中,喷丸加工的模拟基于少弹丸对称分布实现,无法体现实际喷丸加工过程的影响,本发明提出一种考虑多弹丸随机分布的喷丸模型对轮盘寿命的快速计算方法,利用本方法可更加准确快速、自动地模拟真实喷丸情况下残余应力分布情况,并得到相应的寿命增长值。
本发明提供了一种考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法,包括如下步骤:
步骤一:设定喷丸工艺参数,利用二次开发前处理程序建立多丸粒初始位置服从随机分布的喷丸强化有限元模型,并进行有限元分析;
步骤二:利用二次开发后处理程序对步骤一中有限元分析获得的结果进行后处理分析,提取所需求解区域不同深度处的所有应力分布,并对相同深度处所有节点应力计算平均值,作为该深度处的残余应力结果输出;
步骤三:将步骤一中的二次开发前处理程序和步骤二中的二次开发后处理程序整合开发为参数建模二次开发程序,并利用此参数建模二次开发程序实现参数化建模,根据常用的喷丸工艺参数范围值,建立以喷丸工艺参数为输入条件、模拟得到的应力分布为输出的数据库;
步骤四:根据生产所需的喷丸工艺参数值调用步骤三中所建立的数据库,判断生产所需的喷丸工艺参数值是否超出步骤三中所建立的数据库中的喷丸工艺参数阈值;如果超出,执行步骤五,否则执行步骤六;
步骤五:将生产所需的喷丸工艺参数值输入到步骤三中整合开发的参数建模二次开发程序中,通过有限元分析获得随深度变化的残余应力分布,之后采用通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力的 应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,得到设定喷丸工艺参数下的轮盘寿命增加值;
步骤六:将生产所需的喷丸工艺参数值,与步骤三中所建立的数据库中的喷丸工艺参数范围值作对比,从数据库中查找到对应的残余应力分布,根据通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力的 应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,从而得到设定喷丸工艺参数下的轮盘寿命增加值。
进一步,设定的喷丸工艺参数可以包括弹丸直径、弹丸速度和弹丸个数。
本发明的有益效果:
1)本发明考虑了实际喷丸过程中更加准确的分析模型及残余应力分布,从而为考虑喷丸过程的轮盘寿命评估奠定基础;
2)本发明基于随机的原则自动生成随机分布的弹丸,规避了人为生成对称分布的弹丸带来的不准确性与复杂性;
3)本发明通过建立常用不同喷丸工艺参数下残余应力分布数据库,可避免重复的有限元仿真工作,从而使轮盘寿命计算更加简单高效。
附图说明
图1为本发明的考虑多弹丸随机分布的喷丸模型对轮盘寿命的快速计算方法流程图;
图2为本发明考虑多弹丸随机分布的喷丸模型示意图;
图3为本发明考虑喷丸情况下裂纹扩展寿命增长示意图。
具体实施方式
如图1所示,本发明提供的考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法,包括如下步骤:
步骤一:设定喷丸工艺参数,利用二次开发前处理程序,建立多丸粒初始位置服从随机分布的喷丸强化有限元模型,并提交作业进行有限元分析。优选地,喷丸工艺参数包括弹丸直径、弹丸速度和弹丸个数等参数。
步骤二:利用二次开发后处理程序对步骤一中有限元分析获得的结果进行后处理分析,可以实现提取所需求解区域不同深度处的所有应力分布,并对相同深度处所有节点应力计算平均值,作为该深度处的残余应力结果输出(txt文件);
步骤三:将上述步骤一中的二次开发前处理程序和步骤二中的二次开发后处理程序整合开发为一套完整的参数建模二次开发程序,并利用此参数建模二次开发程序实现参数化建模,根据常用的喷丸工艺参数范围值,建立以喷丸工艺参数为输入条件、模拟得到的应力分布为输出的数据库,实现根据喷丸工艺参数快速获取应力并进行后续轮盘寿命分析;
步骤四:根据实际生产所需的喷丸工艺参数值调用步骤三中所建立的数据库,首先判断实际生产所需的喷丸工艺参数值是否超出步骤三中所建立的数据库中的步骤一中设定的喷丸工艺参数的阈值;如果超出,执行步骤五,否则执行步骤六;
步骤五:将实际生产所需的喷丸工艺参数值输入到步骤三中整合开发的参数建模二次开发程序中,进行有限元分析并获得随深度变化的残余应力分布,之后采用通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,得到所需的喷丸工艺参数值下的轮盘寿命增加值;
步骤六:将实际生产所需喷丸工艺参数值,与所建立的应力分布数据库中输入的喷丸工艺参数范围值做对比,找到范围中包含实际生产所需的喷丸工艺参数值的的子数据库,将喷丸工艺参数值输入到子数据库中,快速查找到生产所需的喷丸工艺参数对应的残余应力分布,根据通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,从而得到实际生产所需的喷丸工艺参数值下的轮盘寿命增加值。应该理解,本发明的工作应力是为了计算轮盘寿命人为给定的。
下面结合附图和具体实施例来进一步说明本发明。
选择一组常用的喷丸工艺参数:弹丸直径d=0.28mm、弹丸速度v=80m/s、弹丸个数n=75个,利用二次开发前处理程序,建立多丸粒初始位置服从随机分布的喷丸强化有限元模型,如图2所示,将实际生产所需的喷丸工艺参数值输入到数据库中得到对应的残余应力分布,并选定常用的工作应力 500Mpa进行后续计算。利用权函数法计算残余应力与工作应力应力强度因子K:
其中,a为裂纹半长;σ为不含裂纹时加载的应力分布;F为含有裂纹体几何修正系数;W是带有裂纹体几何结构的宽度;L是带有裂纹体几何结构的长度。
然后将由于施加的载荷和残余应力而产生的应力强度因子K相加,得出总应力强度因子。再利用Nasgro公式计算裂纹扩展速率,计算表达式如下:
其中,N为疲劳载荷的循环次数;a是裂纹长度;R为应力比;ΔK为应力强度因子幅值;Kmax是最大应力强度因子;ΔKth是应力强度因子阈值;KC是临界应力强度因子;C,m,p,q是通过实验确定的材料相关经验参数; f是裂纹张开公式。按照上述本发明轮盘寿命计算方法,获得最终考虑喷丸下的裂纹扩展寿命增长示意图,如图3所示,裂纹扩展寿命增长率约 22.2%。
对于本领域的普通技术人员来说,在不脱离本发明创造构思的前提下,还可以对本发明的实施例做出若干变型和改进,这些都属于本发明的保护范围。
Claims (2)
1.一种考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法,其特征在于,包括如下步骤:
步骤一:设定喷丸工艺参数,利用二次开发前处理程序建立多丸粒初始位置服从随机分布的喷丸强化有限元模型,并进行有限元分析;
步骤二:利用二次开发后处理程序对步骤一中有限元分析获得的结果进行后处理分析,提取所需求解区域不同深度处的所有应力分布,并对相同深度处所有节点应力计算平均值,作为该深度处的残余应力结果输出;
步骤三:将步骤一中的二次开发前处理程序和步骤二中的二次开发后处理程序整合开发为参数建模二次开发程序,并利用此参数建模二次开发程序实现参数化建模,根据常用的喷丸工艺参数范围值,建立以喷丸工艺参数为输入条件、模拟得到的应力分布为输出的数据库;
步骤四:根据生产所需的喷丸工艺参数值调用步骤三中所建立的数据库,判断生产所需的喷丸工艺参数值是否超出步骤三中所建立的数据库中的喷丸工艺参数阈值;如果超出,执行步骤五,否则执行步骤六;
步骤五:将生产所需的喷丸工艺参数值输入到步骤三中整合开发的参数建模二次开发程序中,通过有限元分析获得随深度变化的残余应力分布,之后采用通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力的 应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,得到设定喷丸工艺参数下的轮盘寿命增加值;
步骤六:将生产所需的喷丸工艺参数值,与步骤三中所建立的数据库中的喷丸工艺参数范围值作对比,从数据库中查找到对应的残余应力分布,根据通用权函数法计算残余应力的应力强度因子,再根据应力叠加法将残余应力的 应力强度因子与工作应力强度因子叠加得到总的应力强度因子,利用考虑裂纹闭合效应的NASGRO公式计算裂纹扩展速率,得到设定喷丸工艺参数下的轮盘寿命增加值。
2.根据权利要求1所述的计算方法,其特征在于,设定的喷丸工艺参数包括弹丸直径、弹丸速度和弹丸个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010500736.0A CN111859729B (zh) | 2020-06-04 | 2020-06-04 | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010500736.0A CN111859729B (zh) | 2020-06-04 | 2020-06-04 | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111859729A CN111859729A (zh) | 2020-10-30 |
CN111859729B true CN111859729B (zh) | 2022-07-12 |
Family
ID=72985073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010500736.0A Active CN111859729B (zh) | 2020-06-04 | 2020-06-04 | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111859729B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112329242B (zh) * | 2020-11-09 | 2022-04-08 | 南通大学 | 一种激光喷丸强化处理后材料的疲劳寿命直观评估方法 |
CN115879247B (zh) * | 2023-03-02 | 2023-06-30 | 中国航发四川燃气涡轮研究院 | 一种基于系统辨识的轮盘关键部位应力计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484538A (zh) * | 2014-12-29 | 2015-04-01 | 滁州汽车与家电技术及装备研究院 | 一种基于有限元分析的喷丸强化处理工艺参数确定方法 |
CN106649994A (zh) * | 2016-11-16 | 2017-05-10 | 西北工业大学 | 一种钛合金喷丸强化残余应力场的预测方法 |
CN109359365A (zh) * | 2018-09-30 | 2019-02-19 | 北京航空航天大学 | 一种考虑弹丸随机效应的喷丸工艺数值模拟方法 |
WO2019051616A1 (en) * | 2017-09-18 | 2019-03-21 | Polyvalor, Limited Partnership | METHOD AND SYSTEM FOR REALIZING TRAINING SIMULATION |
CN109725123A (zh) * | 2019-02-28 | 2019-05-07 | 北京航空航天大学 | 一种考虑喷丸强化表层晶粒细化的裂纹扩展寿命确定方法 |
-
2020
- 2020-06-04 CN CN202010500736.0A patent/CN111859729B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484538A (zh) * | 2014-12-29 | 2015-04-01 | 滁州汽车与家电技术及装备研究院 | 一种基于有限元分析的喷丸强化处理工艺参数确定方法 |
CN106649994A (zh) * | 2016-11-16 | 2017-05-10 | 西北工业大学 | 一种钛合金喷丸强化残余应力场的预测方法 |
WO2019051616A1 (en) * | 2017-09-18 | 2019-03-21 | Polyvalor, Limited Partnership | METHOD AND SYSTEM FOR REALIZING TRAINING SIMULATION |
CN109359365A (zh) * | 2018-09-30 | 2019-02-19 | 北京航空航天大学 | 一种考虑弹丸随机效应的喷丸工艺数值模拟方法 |
CN109725123A (zh) * | 2019-02-28 | 2019-05-07 | 北京航空航天大学 | 一种考虑喷丸强化表层晶粒细化的裂纹扩展寿命确定方法 |
Non-Patent Citations (1)
Title |
---|
王帅等.喷丸表面粗糙度对疲劳寿命影响研究.《组合机床与自动化加工技术》.2018,(第2期), * |
Also Published As
Publication number | Publication date |
---|---|
CN111859729A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111859729B (zh) | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 | |
CN107194050A (zh) | 随机载荷作用下涡轮盘结构的概率疲劳寿命预测方法 | |
CN108531714B (zh) | 一种榫接结构激光冲击强化的多重精度优化方法 | |
Hu et al. | Probabilistic design for turbine disk at high temperature | |
Bu¨ che et al. | Automated design optimization of compressor blades for stationary, large-scale turbomachinery | |
CN109583147B (zh) | 一种离心叶轮预旋转盘心梯度模拟件设计方法 | |
CN110990948A (zh) | 一种航空发动机叶片外物损伤疲劳强度预测方法 | |
Frija et al. | Optimization of peened-surface laser shock conditions by method of finite element and technique of design of experiments | |
Chen et al. | Analysis and optimization of process parameter intervals for surface quality in polishing Ti-6Al-4V blisk blade | |
CN104535324A (zh) | 一种发动机低循环寿命控制方法 | |
Liu et al. | Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine | |
CN112528535B (zh) | 基于热-力-流多场耦合的榫槽拉削工艺仿真分析方法 | |
Shlyannikov et al. | Method for assessment of the residual life of turbine disks | |
Paspulati et al. | LCF Initiated-HCF propagated crack life estimation of gas turbine bolts | |
CN110987388B (zh) | 一种基于缺口疲劳强度的等效机械加工缺口的方法 | |
Corran et al. | The development of anomaly distributions for machined holes in aircraft engine rotors | |
Arkhipov et al. | Probabilistic assessment of life for gas turbine engine parts considering manufacture tolerances | |
Kaminskas | SNS mercury target design optimization | |
CN109977462A (zh) | 考虑表面初始加工刀痕的轮盘低循环疲劳寿命预测方法 | |
Whitney-Rawls et al. | Comparison of aero engine component lifing methods | |
Nikhamkin et al. | Effect of blade geometry and foreign object kinetic energy on blades damage | |
Merculov et al. | Improving the accuracy of the behaviour simulation of the material of the turbojet aircraft engine fan rotor blades in the event of a bird strike by using adapted finite element computational models | |
Garrison | High cycle fatigue (HCF) science and technology program 2000 annual report | |
Vasilyev et al. | Modeling of Residual Stresses After Shot Peening and the Effect of Accounting Their Influence on the Stress State of Turbine Blade Fir Tree Roots | |
Kravchenko et al. | Evaluation of optimal grain hardening modes to increase fatigue resistance of parts made of high-strength steels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Li Guo Inventor after: Ding Shuiting Inventor after: Teng Yida Inventor after: Liu Xiaojing Inventor before: Li Guo Inventor before: Ding Shuiting Inventor before: Yida Inventor before: Liu Xiaojing |
|
GR01 | Patent grant | ||
GR01 | Patent grant |