孔隙对3D打印功能梯度材料验证轮盘破裂转速影响分析.pdf
《孔隙对3D打印功能梯度材料验证轮盘破裂转速影响分析.pdf》由会员分享,可在线阅读,更多相关《孔隙对3D打印功能梯度材料验证轮盘破裂转速影响分析.pdf(8页珍藏版)》请在咨信网上搜索。
1、研究论文RESEARCH引文格式:杨雅丽,赖希宁,赵伟,等.孔隙对3D打印功能梯度材料验证轮盘破裂转速影响分析.航空制造技术,2 0 2 3,6 6(13):9 4-10 1.YANG Yali,LAI Xining,ZHAO Wei,et al.Effect of pores on burst speed of functionally graded material validation turbinedisk manufactured by 3D printingJ.Aeronautical Manufacturing Technology,2023,66(13):94-101.孔隙对3
2、D打印功能梯度材料验证轮盘破裂转速影响分析杨雅丽,赖希宁,赵伟”,叶超,陈(1.厦门大学,厦门36 10 0 5;2.中国航发四川燃气涡轮研究院,绵阳6 2 10 0 0)【摘要目前多材料融合3D打印技术对于短寿命、低成本小型发动机的研发具有很好的应用前景,如3D打印功能梯度材料在新型涡轮盘典型热端关键部件上的应用。为探究3D打印技术产生的孔隙缺陷对涡轮盘破裂转速的影响,基于极限应变法开展了50 0 测试均匀温度场及真实温度场下功能梯度材料验证轮盘的破裂转速分析。研究中主要考虑孔隙率、大孔隙所处区域、大孔隙个数、孔隙间距及孔隙与起裂位置的距离等孔隙缺陷表征参量及相关因素对验证轮盘破裂转速的影响
3、。研究结果表明,对于3D打印轮盘的性能分析,不能只考虑孔隙缺陷的随机分布,分布于高应变区(危险截面)的大孔隙将导致验证轮盘破裂转速的显著下降,在3D打印时应严格控制距离预测起裂位置较近的高应变区域内的缺陷。关键词:破裂转速;孔隙;3D打印;功能梯度材料;有限元方法Effect of Pores on Burst Speed of Functionally Graded Material Validation Turbine DiskYANG Yali,LAI Xining,ZHAO Wei?,YE Chao,CHEN Yun,CHEN Lijie2.AECC Sichuan Gas Turbi
4、ne Establishment,Mianyang 621000,China)ABSTRACT At present,the multi-material fusion 3D printing technology has a good application prospect for thedevelopment of short-life and low-cost small engines,such as the application of 3D printing functionally graded materialsin the typical hot section parts
5、 of new-type turbine disk.In order to explore the influence of pore defects generated by3D printing technology on the burst speed of turbine disk in aero-engine,the burst speed analysis of functionally gradedmaterial validation turbine disk was carried out based on the strain-based fracture criteria
6、 under the testing uniformtemperature field of 500 C and the real temperature field respectively.Our studies mainly focus on the influences of poredefect characterization parameters and related factors,such as porosity,the position of large pores located,the numberof large pores,large-pore spacing a
7、nd the distance between large pore and bursts initiation position on the burst speed ofvalidation turbine disk.The results show that the random distribution of pore defects is not the only consideration for theperformance analysis of 3D printed disk,The large pores distributed in the high-strain reg
8、ion(hazard section)will lead toa significant decrease in the burst speed of validation turbine disk,and the defects in the high-strain region,close to thepredicted bursts initiation position,should be strictly controlled during 3D printing.Keywords:Burst speed;Pore;3D printing;Functionally graded ma
9、terial;Finite element methodD0I:10.16080/j.issn1671-833x.2023.13.094*基金项目:中国航发自主创新专项资金(ZZCX-2018-017)。94航空制造技术2023年第6 6 卷第13期云,陈立杰Manufactured by 3D Printing(1.Xiamen University,Xiamen 361005,China;RESEARCH研究论文涡轮盘是航空发动机关键热端零部件之一,涡轮盘破裂将极大影响飞机的飞行安全,造成灾难性飞行事故,所有轮盘和绝大多数轮缘的破裂属于非包容性故障。因此,为保障航空涡轮盘安全性,国内外设计
10、准则通常要求轮盘满足相应的适航要求2-4,当轮盘承受最大温度梯度和最高工作温度时,轮盘的设计破裂转速不低于最高允许稳态转速的12 2%,对涡轮盘展开破裂转速分析是轮盘设计关键内容之一。功能梯度材料涡轮盘是利用功能梯度3D打印技术5将两种材料连续梯度混合制成的整体叶盘,结合了组成材料的性能优点,与传统连接方式焊接的双合金涡轮盘相比,避免了不同材料间明显界面造成的性能突变问题。然而,3D打印成形件内部会不可避免地存在多种缺陷,如孔隙、未熔合缺陷和裂纹等。其中孔隙缺陷是普遍存在的一种缺陷8 ,会降低材料致密度,有时对成形件的力学性能会造成显著影响9-10 ,且孔隙缺陷会成为应力集中及裂纹源,在外部载
11、荷作用下可能产生裂纹并逐步扩展直至断裂失效,造成严重后果。但目前,梯度材料零件3D打印技术普遍存在工艺可控性差、成形效率低及在材料选择及其路线精准受控方面有局限性等缺点,打印工件质量难以把控,成形技术与系统还未十分完善12-13。关于3D打印功能梯度材料涡轮盘的制备还鲜有相关研究,在孔隙缺陷对涡轮盘破裂影响方面更缺少系统的评估,不足以用于确定3D打印轮盘的缺陷容限。本文针对某航空发动机功能梯度材料验证轮盘,探索增材制造孔隙缺陷对验证轮盘超转破裂的影响,为3D打印轮盘的缺陷容限制定提供一定理论依据。1含孔隙验证轮盘建模及破裂转速计算方法1.1含孔隙验证轮盘有限元模型建立为更好地模拟实际3D打印样
12、件内部孔隙缺陷特征,建立随机孔隙验证轮盘模型,利用激光近净成形技术(LENS技术),采用环向扫描策略逐层打印了双金属梯度材料样件,对其进行了内部缺陷的无损检测,通过(a)正视图(b)侧视图图13D打印双金属功能梯度材料样件内部缺陷三维透视图Fig.13D perspective of internal defects of 3D printed bimetallic functionally graded material samples工业CT扫描,获取了内部孔隙缺陷相关数据,样件内部缺陷三维透视图如图1所示。通过三维透视图可看出大孔隙基本分布在样件底部(靠近基板)及样品靠近表面处。实际增材
13、制造会留出一定加工余量,可在后续对增材制造样件打印质量较差的部位进行减材加工,降低样件孔隙率及表面粗糙度。针对本文研究的验证轮盘结构,实际打印将在迎风面及背风面留出余量做减材加工处理,但验证轮盘内腔周围无法进一步机加处理。结合样件缺陷特征CT扫描统计结果,验证轮盘内腔表面随机分布大孔隙,大孔隙半径在0.38 0.9 mm,占孔隙率比例为9 6.2 1%,其余部分设置均匀分布小孔隙,小孔隙半径为0.0 2 5 0.2 8 mm。由于增材制造工艺过程中产生的缺陷具有随机性和弥散性,为了达到产生随机孔隙的目的,采用Python脚本通过随机函数控制在一定范围内生成孔隙,孔隙形状为圆形,且孔隙之间两两互
14、不相交。为保证计算精度,首先对验证轮盘进行网格收敛性测试,考虑到轮盘为周期旋转对称结构,取1/47 的轮盘进行仿真分析,调整全局网格和孔隙网格密度进行含随机孔隙轮盘应力分析,得到验证轮盘最大Mises应力随网格密度的变化情况,如图2 所示。根据计算结果选取全局尺寸为0.6 mm,孔隙处单元数设置32 个,最小尺寸占全局尺寸比例为0.0 1,单元类型为C3D10,网格划分情况如图3所示。1.2功能梯度材料参数设置本文功能梯度材料由镍基高温合金GH4169和GH3625混合而成,两种材料在不同温度下的力学性能见表 1 14-1 和表 2 14,1 考虑到功能梯度材料力学性能沿径向的梯度变化,采用分
15、层法定义功能梯度材料,即将功能梯度材料沿径向等距划分为多个区域,分别进行材料性能赋值。本文将功能梯度材料分为11层,自盘心向叶片方向从GH3625逐渐过渡至GH4169,每个区域两种材料组分变化10%的比例。即整个轮盘的GH3625和GH4169两种材料过渡体积/mm30.2470.1970.1480.0990.0400(c)俯视图2023年第6 6 卷第13期航空制造技术9 5研究论文RESEARCH740720700680(a)最大Mises应力与全局网格密度的关系Fig.2 Relationship between different mesh densities and maximum
16、 Mises stress区域材料混合比逐层变化,为盘心10 0%GH362590%GH3625+10%GH416980%GH3625+20%GH416910%GH3625+90%GH4169 叶片10 0%GH4169,共含有11种不同的材料组分。不同组分比下的材料性能由两种单材料力学性能线性插值得到,编写Python脚本实现功能梯度材料相关参数的赋予,效果如图4所示。本构模型采用简化弹塑性模型,如图5所示。1.3石破裂转速计算方法涡轮盘破裂转速预测的常用方法有平均应力法17-18 、极限应力法、极限应变法、残余变形法等。平均应力法认为当轮盘某一截面上的平均周向或径向应力达到轮盘材料抗拉强度
17、极限时,轮盘发生破裂。该方法虽简便,但计算结果精确度依赖于公式中修正因子的选取,需要温度/弹性模量/GPa2020410020030040050060070080090095096航空制造技术2023年第6 6 卷第13期9559409259108958808651.00.8全局尺寸/mm图2 不同网格密度与最大Mises应力关系表 1 GH4169 材料力学性能 14-15Table1 Mechanical properties of GH4169l-+1泊松比屈服强度/MPa极限强度/MPa0.30120011001760.311600.321500.321410.33最小尺寸占全局尺寸比
18、0.1一0.0 1一0.0010.60.4108010809106101500.212(b)最大Mises应力与孔隙网格密度的关系开展大量试验建立修正因子的确定方法119 。基于弹塑性有限元分析方法,采用局部破裂准则预测轮盘破裂转速的极限应变法的破裂准则为:在单调增孔隙网格划分盘心图31/47 含随机孔隙验证轮盘网格划分模型Fig.3 Meshing of 1/47 validation turbine disk with random pores延伸率线膨胀系数(10-)14400.2213300.1613000.2013100.1910000.266150.261800.8020孔隙网格单
19、元数/个腹板11.811.813.013.514.114.414.815.417.018.42836叶片密度/(gcm)8.2444RESEARCH研究论文表 2 GH3625 材料力学性能 14,16 1Table 2 Mechanical properties of GH362514,1温度/弹性模量/GPa202051002003004005006007008009001000加的载荷作用下,当轮盘内任一局部点的等效塑性应变达到材料的延伸率时,该局部点就会有裂纹萌生,一旦裂纹形成就会快速扩展,使轮盘失去承载能力,轮盘发生破裂,此时的转速为轮盘的破裂转速。当极限应变法计算中采用大应变有限元
20、分析方法时,轮盘的破裂准则可以描述为:在单调增加的转速作用下,轮盘内任一局部点的等效塑性应变&n达到ln(1+8)时(8 为材料延伸率),该点发生破坏2 0 。极限应变法与平均应力法相比,其主要通过轮盘材料的弹塑性本构模型来预测轮盘的破裂转速,极限应变法不包含任何修正因子,不需要像平均应力法通过开展大量的实际轮盘破裂试验来确定轮盘的修正系数,方便应用于试验前期阶段的轮盘设计与优化。与此同时,极限应变法不仅可以较精确地描述轮盘破裂前的变形行为,清晰地呈现出轮盘破裂的最薄弱部位,破裂转速预测结果与实际测试结果更接近2 0-2 1。因此,本文将采用极限应变法开展含随机孔隙验证轮盘的超转破裂分析。1.
21、4载荷及边界条件设置航空发动机涡轮盘一般要承受离心载荷及迎面而来的气流对叶片的冲击产生的气动载荷、温度载荷、振动载荷和装配应力等多种载荷的作用。实际工况下涡轮盘的离心载荷和温度载荷占涡轮盘所承受载荷的主要部分,其他气动载荷、振动载荷及装配应力的影响相对较小,故本文对验证轮盘的破裂转速分析主要考虑离心载荷和温度载荷,其他载荷忽略不计。顺序耦合法是目前涡轮盘多物理场耦合分析的常用方法2-2 4,该方法是将前一个物理场的计算结果当成后一个物理场计算的初始载荷条件,按照顺序单独计算泊松比屈服强度/MPa0.315985561850.311790.321730.321680.33极限强度/MPa9719
22、4043484547382948472239240398104Fig.4 Schematic diagram of multi-material distribution of00.200.2%图5简化弹塑性模型应力-应变曲线Fig.5 Stress-strain curves of simplified elastic-plastic model各物理场作用结果的耦合计算方法。本文进行顺序耦合分析时,把预先计算完成的温度场结果作为预定义场加入初始条件中,再设置离心载荷进行热力耦合计算。真实温度场根据150 0 0 个完整叶轮盘实测点的温度数据值拟合得到(图6),拟合结果显示轮心温度最低,其值
23、为2023年第6 6 卷第13期航空制造技术9 7延伸率线膨胀系数(10-)0.4012.712.70.5313.013.50.6713.60.6114.014.30.6715.00.6115.616.10.78IGH36251GH4169多材料区域图4验证轮盘多材料分布示意图validationturbinedisk8/2应变密度/(gcm)8.44等效塑性应变(Avg:75%)研究论文RESEARCH281.8,叶片温度最高,为9 8 9.1。固定验证轮盘的轴向和周向位移,并在对称面施加循环对称边界条件。2破裂转速计算结果分析2.1不同温度场下孔隙率对验证轮盘破裂转速影响由于条件限制,验证
24、轮盘的超转试验较难模拟真实温度环境,实际测试时通常在均匀温度场下进行,并根据与实际温度场下破裂比相当的原则来确定测试转速。温度/+9.890e+02+9.300e+02+8.711e+02+8.122e+02+7.532e+02+6.943e+02+6.354e+02+5.764e+02+5.175e+02+4.586e+02+3.996e+02+3.407e+02+2.818e+02Fig.6 Contour of real temperature field of validation turbine diskMax:+5.856e-01GH4169GH4169(a)无孔隙(b)孔隙率0.
25、1%图7 真实温度场下等效塑性应变云图Fig.7 Contour of equivalent plastic strain in real temperature fieldMax:+2.183e-01Max:+3.030e-0190%GH3625十10%GH4169(a)无孔隙Fig.8 Contour of equivalent plastic strain at the testing uniform temperature field of 500 C98航空制造技术2023年第6 6 卷第13期因此,为探究不同工况下孔隙对验证轮盘性能的影响,本文分别进行了真实温度场与50 0 测试均
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 孔隙 打印 功能 梯度 材料 验证 轮盘 破裂 转速 影响 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。