CN111797554B - 一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 - Google Patents
一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 Download PDFInfo
- Publication number
- CN111797554B CN111797554B CN202010610897.5A CN202010610897A CN111797554B CN 111797554 B CN111797554 B CN 111797554B CN 202010610897 A CN202010610897 A CN 202010610897A CN 111797554 B CN111797554 B CN 111797554B
- Authority
- CN
- China
- Prior art keywords
- mortise
- model
- shot
- dislocation
- nozzle
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种涡轮榫槽喷丸离散元‑有限元耦合多尺度模拟方法,步骤为:(1)根据榫槽结构对称性对欲研究的榫槽结构进行切割,赋予高应变率本构模型;(2)根据实际工艺情况设置喷口与榫槽的相对位置,定义喷口尺寸、喷口角度,完成榫槽与喷口建模;(3)基于高应变率下应力应变数据获取位错演化模型参数,建立榫槽材料的位错演化模型,关联晶粒尺寸、位错密度与宏观应力应变;(4)利用DEM模块定义粒子生成器*Particle,设置弹丸参数,将榫槽与喷口进行组合,设置弹丸与榫槽之间接触特性;(5)进行数值仿真,获取榫槽结构强化之后表面完整性参数。
Description
技术领域
本发明是一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法,它是一种能够高效准确模拟出涡轮榫槽结构喷丸强化效果的方法,属于航空航天发动机技术领域。
背景技术
喷丸强化技术是目前航空工业应用最为成熟广泛的表面强化技术,利用高速弹丸流喷射结构表面发生碰撞,引入可观的残余应力场,改变表面形貌并产生可控的塑性变形。涡轮榫槽部位长期服役在高温、高压、高转速的环境中,服役环境恶劣,极易诱发产生疲劳问题,因此喷丸强化在涡轮榫槽上的应用具有极大的工程价值。目前针对喷丸强化的研究仍以试验为主,喷丸强化有限元数值模拟的研究主要聚焦在平面结构,而对于涡轮榫槽这类几何特征复杂、曲面特征明显结构的数值模拟,单纯的有限元法存在计算效率低下、弹丸建模困难的问题,难以有效模拟涡轮榫接的喷丸强化过程。离散元方法作为一种无网格方法,特别适用于喷丸强化这种颗粒繁多的力学过程的仿真,因此将离散元方法与有限元相结合,发展一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法,是解决涡轮榫槽结构喷丸强化模拟的有效思路。
目前关于涡轮榫槽喷丸强化的研究基本以试验为主,但是如果针对涡轮榫槽这种几何复杂的结构进行残余应力、粗糙度及微观组织的检测,试验成本较高。而且试验存在次数限制,无法针对所有喷丸强化参数进行试验。因此发展有效的涡轮榫槽数值模拟方法是必须的。
发明内容
本发明技术解决方案:克服现有技术的不足,提供一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法,有效反映涡轮榫槽喷丸强化后残余应力场分布,服务与支撑喷丸强化在涡轮榫槽部位的应用与后续疲劳寿命评估。
本发明技术解决方案:一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法,结合位错演化模型,利用大型商用有限元软件ABAQUS中定义关键字的方式,实现弹丸不断的产生并冲击靶材,有效提高计算效率,实现步骤如下:
第一步,根据榫槽模型对称性,对欲研究的榫槽结构进行切割,赋予高应变率本构模型。真实涡轮盘榫槽结构在涡轮盘周向具有周期对称性,同时同一个榫槽沿轴线对称面具有对称性,因此涡轮榫槽的喷丸强化模拟仅仅需要针对半个榫槽进行。在UG等建模软件中对模型进行切分,采用能够反应高应变率下应力应变关系的J-C本构模型,其表达式如下:
第二步,根据实际涡轮榫槽进行喷丸强化工艺处理时喷口的真实情况,包括喷口尺寸、喷口角度、喷口与榫槽结构的相对位置等,利用ABAQUS有限元软件建立面单元,在装配界面对喷口角度、喷口与榫槽结构的相对位置等进行调整,保证与真实工艺的一致性。
第三步,基于高应变率下应力应变数据获取位错演化模型参数,建立待研究材料的位错演化模型,关联晶粒尺寸、位错密度与宏观应力应变;开展高应变率下霍普金森压杆试验或者根据已经建立的待研究材料高应变率本构模型,获取应力应变数据,通过拟合得到材料的位错演化模型参数,用以计算喷丸过程材料表层的晶粒细化;设置位错演化模型需要拟合的参数的初始值,代入位错演化模型中计算得到应力。然后将计算得到的高应变率下的应力结果与试验测定或者模型计算得到的结果进行对比,利用遗传算法将得到的结果误差减小到最小值,从而得到材料的位错演化模型参数;位错演化模型表达式如下:
式中,ρc、ρw分别为位错胞内、胞壁中的位错密度,分别表示位错胞内、胞壁中的位错密度变化率;α*、β*、k0和n*分别为控制位错密度演化率的参数;b为伯氏矢量;分别为位错胞内、胞壁中的分切剪应变率,由于胞壁、胞内的边界需保证应变协调,通常认为这两个参量相同,即 为工艺的参考分切剪应变率;d为晶粒尺寸;f为位错胞壁所占的体积分数,通常可通过透射电子显微镜试验检测获得。
将上述方程利用ABAQUS有限元软件的VUSDFLD子程序进行编程,方便后续计算。
第四步:在ABAQUS的DEM模块定义中粒子生成器*Particle,在第二步生成的inp文件的基础上,设置第二步建立的喷口为粒子生成器,根据实际工艺情况,设置弹丸的尺寸大小、弹丸尺寸满足的分布方式,设置弹丸的弹性参数以及从喷口处产生的流量与速度,定义弹丸与弹丸之间满足赫兹接触。将榫槽与喷口进行组合,设置弹丸与榫槽之间切向罚函数接触与法向硬接触。上述操作利用Python脚本进行,方便快捷。
第五步,基于ABAQUS软件建立有限元模型进行喷丸强化数值仿真,根据实际榫槽设置边界条件,利用ABAQUS的Python脚本语言进行编程,对靶材进行网格划分,冲击区域进行网格细化。计算时同时提交Fortran程序,完成常规计算的同时计算出残余应力、晶粒尺寸的结果。
本发明与现有技术相比的优点在于:
本发明将离散元-有限元耦合的方法引入涡轮榫槽结构的喷丸强化多尺度模拟中,通过位错演化模型实现宏观参数与微观参数的关联,目前相关技术未见报道,本发明弥补了相关研究不足。
附图说明
图1为本发明的一种涡轮榫槽结构喷丸强化离散元-有限元耦合多尺度模拟方法实施流程;
图2为本次举例的涡轮榫槽喷丸强化模拟有限元模型示意图;
图3为本次举例的涡轮榫槽喷丸强化模拟的残余应力结果。
具体实施方式
下面结合附图,通过举例的方式,对本发明一种涡轮榫槽结构喷丸强化离散元-有限元耦合多尺度模拟方法的技术方案做进一步说明。本实施例的预研究材料为镍基高温合金GH4169。如图1所示,本发明一种涡轮榫槽结构喷丸强化离散元-有限元耦合多尺度模拟方法具体实现过程如下:
第一步根据榫槽模型对称性,对欲研究的榫槽结构进行切割,赋予高应变率本构模型。所述高应变率本构模型是指喷丸强化过程应变率达到103~105/s,所用本构模型需要反应该应变率下应力应变行为真实涡轮盘榫槽结构在涡轮盘周向具有周期对称性,同时同一个榫槽沿轴线对称面具有对称性,因此涡轮榫槽的喷丸强化模拟仅仅需要针对半个榫槽进行。本次模拟的榫槽在UG建模软件中对模型进行切分。采用能够反应高应变率下应力应变关系的J-C本构模型,其余高应变率本构模型的应用类似。J-C模型的表达式为:
预研究材料GH4169的J-C模型参数为A=1241MPa,B=622,C=0.0134,m=1.05,n=0.6522。给定应变、应变率与温度项,可以得到材料的流动应力。
第二步,根据实际涡轮榫槽进行喷丸强化工艺处理时喷口的真实情况,包括喷口尺寸、喷口角度、喷口与榫槽结构的相对位置等,利用ABAQUS有限元软件建立面单元,在装配界面对喷口角度、喷口与榫槽结构的相对位置等进行调整,保证与真实工艺的一致性。本次模拟的榫槽冲击主要是榫齿面,根据真实工艺情况,设置喷口与轴向夹角为45°,喷口与榫槽距离为15mm,喷口尺寸为10mm,建模完成之后,生成inp文件。
第三步,基于高应变率下应力应变数据获取位错演化模型参数,建立待研究材料的位错演化模型,关联晶粒尺寸、位错密度与宏观应力应变;开展高应变率下霍普金森压杆试验或者根据已经建立的待研究材料高应变率本构模型,获取应力应变数据,通过拟合得到材料的位错演化模型参数,用以计算喷丸过程材料表层的晶粒细化;利用SEM、TEM等技术对材料的晶粒尺寸与位错密度进行测量,将测量结果作为输入,根据位错胞内部位错密度的演化规律。设置位错演化模型需要拟合的参数的初始值,按体积分数叠加计算出整体分解剪应力,代入模型中计算得到应力。然后将计算得到的高应变率下的应力结果与J-C模型计算得到的结果进行对比,设置最大允许的拟合误差,利用遗传算法不断更新位错演化模型的参数,将两种模型得到的结果误差减小到允许值,从而得到材料的位错演化模型参数。位错演化模型表达式如下:
式中,ρc、ρw分别为位错胞内、胞壁中的位错密度;α*、β*、k0和n*分别为控制位错密度演化率的参数;b为伯氏矢量;分别为位错胞内、胞壁中的分切剪应变率,由于胞壁、胞内的边界需保证应变协调,通常认为这两个参量相同,即 为工艺的参考分切剪应变率;d为晶粒尺寸;f为位错胞壁所占的体积分数,通常可通过透射电子显微镜试验检测获得。
将上述方程利用ABAQUS有限元软件的VUSDFLD子程序进行编程,方便后续计算。
第四步:在ABAQUS的DEM模块定义中粒子生成器*Particle,在第二步生成的inp文件的基础上,设置第二步建立的喷口为粒子生成器,根据实际工艺情况,设置弹丸的尺寸大小、弹丸尺寸满足的分布方式,设置弹丸的弹性参数以及从喷口处产生的流量与速度,定义弹丸与弹丸之间满足赫兹接触。将榫槽与喷口进行组合,设置弹丸与榫槽之间切向罚函数接触与法向硬接触。在本次模拟中,设计喷口正向为弹丸产生方向,弹丸直径为0.4mm,尺寸均布,弹丸流量为5kg/s,弹丸流速为62m/s,弹丸流量与流速通过定义幅值曲线获得。弹丸与弹丸之间接触特征为赫兹接触,弹丸与榫槽设置为切向罚函数接触,摩擦系数为0.3,法向硬接触。弹丸产生的示意图如图2。
第五步,基于ABAQUS软件建立有限元模型进行喷丸强化数值仿真,根据实际榫槽设置边界条件,本次模拟设置的边界条件为底面全约束。利用ABAQUS的Python脚本语言进行编程,对靶材进行网格划分,冲击区域进行网格细化。计算时同时提交Fortran程序,完成常规计算的同时计算出残余应力、晶粒尺寸的分布情况,完成多尺度模拟。残余应力分布如图3所示,
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (1)
1.一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法,其特征在于,包括步骤如下:
步骤(1):根据榫槽模型对称性,对欲研究的榫槽结构进行切割,赋予高应变率本构模型;所述榫槽模型对称性是指榫槽沿涡轮盘周向存在周期对称性,同时同一个榫槽关于径向存在对称性;所述切割是指在建模软件或者有限元软件中根据模型对称性切割榫槽模型,减小后续计算量;所述高应变率本构模型是指喷丸强化过程应变率达到103~105/s,所述本构模型反应该应变率下应力应变行为;
步骤(2):有限元建模生成喷口,根据榫槽模型,考虑实际工艺情况设置喷口与榫槽的相对位置与榫槽喷丸区域,定义喷口尺寸、喷口角度;所述实际工艺情况是指真实榫槽结构在进行喷丸处理的工艺过程情况,包括喷口与榫槽的相对位置、喷口尺寸、喷口角度;
步骤(3):基于高应变率下应力应变数据获取位错演化模型参数,建立榫槽材料的位错演化模型,关联晶粒尺寸、位错密度与宏观应力应变;所述位错演化模型参数是指基于位错理论建立的描述位错演化过程的模型中未知的材料参数;
步骤(4):利用离散元模块定义粒子生成器*Particle,设置弹丸参数,将榫槽与喷口进行组合,设置弹丸与榫槽之间接触特性,建立离散元-有限元耦合模型;所述离散元模块是指利用离散元方法处理问题的模块;所述粒子生成器*Particle是指在用于生成粒子的关键字*Particle;所述弹丸参数是指弹丸弹性模量、泊松比、流量、速度、弹丸间接触特性;
步骤(5):进行数值仿真,利用步骤(4)建立的离散元-有限元耦合模型,获取榫槽结构强化之后表面完整性参数;所述表面完整性参数包括宏观参数残余应力以及微观参数晶粒尺寸;所述第(1)步中,在建模软件中对模型进行切分,采用能够反应高应变率下应力应变关系的J-C本构模型,其表达式如下:
所述步骤 (3)中位错演化模型表达式如下:
式中,ρc、ρw分别为位错胞内、胞壁中的位错密度,分别表示位错胞内、胞壁中的位错密度变化率;α*、β*、k0和n*分别为控制位错密度演化率的参数;b为伯氏矢量;分别为位错胞内、胞壁中的分切剪应变率,由于胞壁、胞内的边界需保证应变协调,这两个参量相同,即 为工艺的参考分切剪应变率;d为晶粒尺寸;f为位错胞壁所占的体积分数,通过透射电子显微镜试验检测获得;步骤(4)中设置的弹丸相互间接触特性满足赫兹接触,相应的设置方法是通过修改inp文件进行的;所述步骤(4)中:
在ABAQUS的DEM模块定义中粒子生成器*Particle,在第二步生成的inp文件的基础上,设置第二步建立的喷口为粒子生成器,根据实际工艺情况,设置弹丸的尺寸大小、弹丸尺寸满足的分布方式,设置弹丸的弹性参数以及从喷口处产生的流量与速度,定义弹丸与弹丸之间满足赫兹接触;设置的弹丸与榫槽之间接触特性包括切向接触与法向接触特性两类;将榫槽与喷口进行组合,设置弹丸与榫槽之间切向罚函数接触与法向硬接触;步骤(4)基于ABAQUS软件建立有限元模型进行喷丸强化数值仿真,根据实际榫槽设置边界条件,设置的边界条件为底面全约束,对靶材进行网格划分,冲击区域进行网格细化;完成常规计算的同时计算出残余应力、晶粒尺寸的分布情况,完成多尺度模拟。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010610897.5A CN111797554B (zh) | 2020-06-30 | 2020-06-30 | 一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010610897.5A CN111797554B (zh) | 2020-06-30 | 2020-06-30 | 一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111797554A CN111797554A (zh) | 2020-10-20 |
CN111797554B true CN111797554B (zh) | 2022-03-15 |
Family
ID=72809863
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010610897.5A Active CN111797554B (zh) | 2020-06-30 | 2020-06-30 | 一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111797554B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112338814B (zh) * | 2020-10-29 | 2022-06-07 | 中国航发南方工业有限公司 | 涡轮盘的复合喷丸加工方法 |
CN112560259A (zh) * | 2020-12-11 | 2021-03-26 | 北京航空航天大学 | 一种基于弹塑性碰撞的传动轴喷丸表面残余应力快速预测方法 |
CN115828472B (zh) * | 2023-02-15 | 2023-05-26 | 太原理工大学 | 一种用于模拟滚磨光整加工工件表面残余应力的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102862124B (zh) * | 2012-07-26 | 2015-09-09 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种弧形榫槽的喷丸方法 |
CN103898313B (zh) * | 2014-04-10 | 2016-01-13 | 西安航空动力股份有限公司 | 一种涡轮盘榫槽结构的激光冲击强化方法 |
CN104484538A (zh) * | 2014-12-29 | 2015-04-01 | 滁州汽车与家电技术及装备研究院 | 一种基于有限元分析的喷丸强化处理工艺参数确定方法 |
CN108531714B (zh) * | 2018-07-04 | 2019-05-10 | 北京航空航天大学 | 一种榫接结构激光冲击强化的多重精度优化方法 |
CN109946006B (zh) * | 2019-01-17 | 2019-09-24 | 中国石油大学(华东) | 基于混合硬化模型的微粒流冲击金属材料力学行为预测方法 |
CN109725123B (zh) * | 2019-02-28 | 2020-05-22 | 北京航空航天大学 | 一种考虑喷丸强化表层晶粒细化的裂纹扩展寿命确定方法 |
-
2020
- 2020-06-30 CN CN202010610897.5A patent/CN111797554B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111797554A (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111797554B (zh) | 一种涡轮榫槽喷丸离散元-有限元耦合多尺度模拟方法 | |
Liu et al. | Aerospace applications of laser additive manufacturing | |
CN112100885B (zh) | 一种高能喷丸表面硬度数值模拟方法 | |
CN112036071A (zh) | 一种考虑表面粗糙度的有限元与离散元相耦合的喷丸强化数值模拟方法 | |
CN112765732A (zh) | 一种基于选区激光熔化工艺的航空叶片拓扑优化设计方法 | |
CN112580228B (zh) | 一种涡扇发动机混合结构风扇叶片结构设计优化方法 | |
Mao et al. | Finite element modeling of a bird striking an engine fan blade | |
CN109583037A (zh) | 一种航空发动机叶片喷丸加工变形的参数控制方法 | |
CN111985131B (zh) | 一种考虑曲面形状的涡轮榫槽激光冲击强化数值模拟方法 | |
Zhang et al. | Notch-based probabilistic fatigue analysis of automobile transmission gear considering size effect | |
Taylor | Industrial perspectives on geometry handling for aerodynamics | |
EP3742246A1 (en) | Computer implemented method and system for generating slope synchronized tool paths for incremental sheet forming | |
Sudhir Sastry et al. | Impact analysis of compressor rotor blades of an aircraft engine | |
Xu et al. | Control of mistuning level using hard coating for blisks based on mistuning identification technique | |
Kulakov et al. | The assessment of impact of the crack size on the fracture load of a cylindrical element | |
Cao et al. | DNN-based surrogate modeling-based feasible performance reliability design methodology for aircraft engine | |
Briasco et al. | A NURBS-based optimization tool for axial compressor cascades at design and off-design conditions | |
Caivano et al. | Defect-Driven Topology Optimisation: TopFat algorithm validation via 3D components re-design for real industrial applications | |
Ou et al. | Reduction in post forging errors for aerofoil forging using finite element simulation and optimization | |
Wang et al. | Optimization of laser cladding process parameters based on genetic algorithm and neural networks | |
CN109902349B (zh) | 自由曲面薄壁叶片分层分区域加载残余应力的建模方法 | |
Duo et al. | Assessment of foreign object damage (FOD) to aero engine blades | |
Brown et al. | Effect of Additive Manufacturing Induced Material Stiffness Anisotropy on Topology Optimisation | |
Abrusci | Preliminary optimization of an hollow Low Pressure Turbine blade | |
Chaouch et al. | Check for updates Performance Prediction of Abrasive Water Jet Machining Composite Laminates Using Artificial Neural Networks and Regression |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |