航空发动机螺栓连接结构薄层单元静力学非线性模型修正_刘玉.pdf
《航空发动机螺栓连接结构薄层单元静力学非线性模型修正_刘玉.pdf》由会员分享,可在线阅读,更多相关《航空发动机螺栓连接结构薄层单元静力学非线性模型修正_刘玉.pdf(9页珍藏版)》请在咨信网上搜索。
1、文章编号:1000-8055(2023)06-1480-09doi:10.13224/ki.jasp.20220899航空发动机螺栓连接结构薄层单元静力学非线性模型修正刘玉1,赵迪文1,艾延廷1,付鹏哲2,刘仕运2(1.沈阳航空航天大学辽宁省航空推进系统先进测试技术重点实验室,沈阳110136;2.中国航发沈阳发动机研究所,沈阳110066)摘要:针对于航空发动机机匣带有止口的螺栓连接结构的非线性刚度特征提出一种新的非线性模拟方法以及模型修正方法。首先提出了一种非线性薄层单元方法对螺栓连接结构的非线性刚度特征进行模拟,并对模拟结果进行评估;然后,提出一种新的双加权响应面修正方法;最后,使用该修
2、正方法对螺栓连接结构静力学非线性模型修正:结果表明:本文建立的非线性薄层单元简化模型可以预测实际螺栓连接结构考虑止口接触面粗糙度的非线性刚度特征;提出双加权响应面修正方法修正前后的刚度误差维持在 3.5%以内,验证了建模以及修正方法的可行性。关键词:螺栓连接结构;有限元建模;非线性刚度模拟;非线性模型修正;刚度测试中图分类号:V232文献标志码:AModificationofstaticnonlinearmodelofthinlayerelementforaeroengineboltedstructureLIUYu1,ZHAODiwen1,AIYanting1,FUPengzhe2,LIUSh
3、iyun2(1.LiaoningKeyLaboratoryofAdvancedMeasurementandTestTechnologyforAviationPropulsionSystem,ShenyangAerospaceUniversity,Shenyang110136,China;2.ShenyangEngineResearchInstitute,AeroEngineCorporationofChina,Shenyang110066,China)Abstract:Anewnonlinearsimulationmethodandmodelmodificationmethodwereprop
4、osedforthenonlinearstiffnesscharacteristicsofboltedstructureswithnotchesinaeroenginecasings.Anonlinearthin layer element method was proposed to simulate the nonlinear stiffness characteristics of boltedstructures,andthesimulationresultswereevaluated.Anewdualweightedresponsesurfacecorrectionmethodwas
5、alsoproposed.Themodifiedmethodwasusedtomodifythestaticnonlinearmodelofboltedstructures.The results showed that the simplified nonlinear thin layer element model can predict thenonlinearstiffnesscharacteristicsofactualboltedstructuresconsideringtheroughnessoftheseamcontactsurface;adualweightedrespons
6、esurfacecorrectionmethodwasproposedtomaintainthestiffnesserrorwithin3.5%beforeandaftercorrection,verifyingthefeasibilityofthemodelingandcorrectionmethod.Keywords:boltedstructure;finiteelementmodeling;nonlinear-stiffnesssimulation;nonlinearmodelmodification;stiffness-test收稿日期:2022-11-23基金项目:辽宁省教育厅面上项
7、目(JYT2020053);沈阳航空航天大学博士启动基金(19YB39)作者简介:刘玉(1989),男,讲师,博士,主要从事航空发动机结构强度、力学振动研究。引用格式:刘玉,赵迪文,艾延廷,等.航空发动机螺栓连接结构薄层单元静力学非线性模型修正J.航空动力学报,2023,38(6):1480-1488.LIUYu,ZHAODiwen,AIYanting,etal.ModificationofstaticnonlinearmodelofthinlayerelementforaeroengineboltedstructureJ.JournalofAerospacePower,2023,38(6):
8、1480-1488.第38卷第6期航空动力学报Vol.38No.62023年6月JournalofAerospacePowerJune2023有限元分析方法因其理论基础简明、物理概念清晰、模拟精度高等优点,被广泛应用于许多航空发动机性能研究中。准确模拟参数对航空发动机力学特征分析至关重要。然而,在有限元建模过程中必然存在一定的阶次误差和参数误差1-2,所以需要进行模型修正3-5。在过去的 30 年里,有限元模型修正技术在航空发动机性能研究方面发展迅速,各种新理论在该领域得到广泛应用。根据修正目标的不同,可以分为直接修正技术和参数修正技术6-7。与传统的直接修正技术相比,参数修正技术由于修正后物
9、理意义明确,在工程实践中普遍采用。Stentson8确定了模态参数以及变体积结构参数,最后计算出修正后的有限元模型的刚度矩阵和质量矩阵。Kabe9做了进一步的研究,保留了修正后的计算模型刚度矩阵。Kabe 的研究是一个很大的进步,但这种方法只适用于自由度较少的系统。Mottershead 和 Friswell等10-11对参数模型修正技术进行了更加细致的研究,极大地推动了设计参数修正技术的发展。任伟新等12研究了基于 2 阶响应法的结构静力学有限元模型修正方法;邓苗毅、马双超等13-15以两跨连续梁结构为研究背景,验证了模型修正的计算效率,获得了高精度的仿真计算模型。总之,基于响应法的有限元模
10、型修正是在修正理论的基础上结合统计理论发展起来的,具有较高的修正精度,并且提高了计算效率。但是在修正优化迭代的过程中,需要不断地调用仿真模型,每进行一次优化,就要调用一次仿真模型,当涉及多参数目标进行动力学优化时,就会导致所需时间长,计算量大,动力学特征不明显等问题。为解决上述问题,基于实测数据提高模型修正的效率受到极大关注。频率特征值、振动模态振型、频响函数等动力学实测数据在试验测试过程中易受噪声影响,且对系统结构的局部损伤不太敏感,相比之下,静力学实测数据具有测试准确、受噪声干扰小等优点,使得基于静力测量数据的模型修正技术具有广阔的应用发展前景16。2003 年向天宇等17以实测形变响应值
11、为基准,结合灵敏度分析方法,对两跨连续梁结构进行有限元模型修正和损伤识别,取得了很好的识别效果。2004 年邱春图等18以静力学试验测试响应值为依据,对飞机整体有限元模型进行静力学修正研究,并提出了实测数据的可靠性评估和确认方法。2007 年 Lee 等19基于静力测量数据对连接翼结构有限元初始模型进行迭代优化计算,使修正后的有限元模型与试验测试结果具有较好的一致性。2009 年邓苗毅等20以实际桥梁结构为背景,利用实测静力载荷数据对桥梁结构进行静力学模型修正研究,研究结果为同类桥梁结构的实际静力性能的评价工作提供了依据。2015 年翟学、王建军等21-22基于实测静刚度试验数据对某型航空发动
12、机静子机匣结构进行修正研究,使修正后的静子系统响应值与试验测试值十分接近,取得了较好的效果。现有研究中,普通的响应面修正方法只能对线性问题进行修正,针对于静力学非线性刚度的修正方法仍有待完善。准确建立螺栓连接结构结合面的等效模型一直以来都是国内外学者研究的热点。Song 等23提出一种调整的 Iwan 梁单元来模拟结合面,并通过实时机器学习力场(MLFF)神经网络识别出该单元中的参数。徐超等24通过改进的 Iwan 模型来建立非线性连接单元,分析了螺栓连接梁的非线性动力学特征。Liao 等25提出了一种使用虚拟梯度材料模拟结合面的方法,并且基于胡克定律推导出材料的弹性模量和切变模量。Beaud
13、oin等26提出一种非线性的块状模型,该模型将连接件的几何参数纳入考虑范畴,并且表征出了接触面的间隙和摩擦行为。薄层单元由于原理简单,易在有限元软件中建模计算,因此也是目前最常用的结合部等效模型之一。Zhao 等27探讨了螺栓预紧力对固有频率的影响,以薄层单元表示界面接触,对比试验和有限元分析,得到了接触刚度与螺栓预紧力的关系。姜东等28基于薄层单元理论对螺栓搭接结构进行有限元建模,且对薄层单元参数进行识别。姚星宇、刘玉等29-31针对航空发动机中的螺栓连接结构,研究了薄层单元材料参数对连接刚度的影响规律,并推导出薄层单元材料参数的计算公式。本文以某航空发动机机匣带有止口的螺栓连接结构为对象,
14、基于静力学模型修正的基本原理,利用薄层单元对刚度矩阵的非线性特征进行模拟,对模型进行验证后,依据实测静力学响应信息,提出一种新的双加权响应面修正方法,并完成静力学非线性模型修正。本文建立的薄层单元简化模型结构简单、计算方便,具有良好的静力学特第6期刘玉等:航空发动机螺栓连接结构薄层单元静力学非线性模型修正1481征,形成的参数型修正方法可以避免对模型进行多次调用,提高了模型修正的效率,具有很强的工程适用性。1基于薄层单元法的螺栓连接结构非线性建模方法1.1薄层单元基本原理对螺栓连接处进行薄层单元法建模。假设分块的薄层单元为 8 节点六面体实体单元,单元中任一点形变可表示为|x=8i=1Nixi
15、y=8i=1Niyiz=8i=1Nizi(1)Nixiyizii式中是 8 节点六面体单元第 i 个节点的形函数,、和 是第 个节点坐标。根据弹性力学基本方程,单元应变向量 为=BaT(2)Baa式中 为单元几何矩阵,为单元中的节点形变向量。单元应力向量 与节点形变向量 的关系为=D=DBaT(3)D式中为弹性矩阵。对于螺栓连接结构对接面,忽略各个方向间的耦合刚度,则弹性矩阵变为D=|c11c22c33c44c55c66|(4)c11c22c33c44c55c66式中、代表三个方向的弹性模量,、代表三个平面的切变模量,弹性矩阵所在坐标系与节点坐标系一致。根据虚功原理,单元刚度矩阵为Ke=V0B
16、TDBdV0(5)V0KeFeKF式中是单元体积。最后,集合单元刚度矩阵和节点力向量形成螺栓连接机匣法兰结构的刚度矩阵和螺栓连接机匣法兰结构的节点力向量,从而静力学有限元方程可表示为KaT=FT(6)|Nix=AixNiy=AiyNiz=Aiz(7)AixAiyAiz式中、为节点各个方向形函数Ni的偏导数。根据薄层单元法的有限元方程,需要明确航空发动机螺栓连接结构薄层单元法的力学参数,该力学参数可以表征螺栓连接结构的刚度特性。如图 1 所示,螺栓连接机匣法兰结构的轴向刚度定义为kn=Fnl(8)Fnl式中是结构所受的轴向力,是整个结构的轴向形变。Fn图1螺栓连接结构轴向刚度定义Fig.1Def
17、initionofaxialstiffnessofboltedstructures当螺栓连接机匣法兰结构受到轴向力作用时,整个结构的横向形变几乎为 0,可以忽略不计,所以在只考虑轴向刚度问题时,根据式(8)螺栓连接机匣法兰结构的轴向刚度可以进一步表示为kn=Fnl=ki=1V0(A1xc11)dV0Xil+ki=1V0(Aixc11)dV0Xil+ki=1V0(A(i+1)xc11)dV0X(i+1)l+ki=1V0(Aixc11)dV0Xil(9)kk8Xiki=1V0(A1xc11)dV0Xil+ki=1V0(Aixc11)dV0Xil式中 为某个节点被单元共有的单元数,为 对 应 单 元
18、 的 轴 向 形 变。部分代表不同扇区1482航空动力学报第38卷ki=1V0(A(i+1)xc11)dV0X(i+1)l+ki=1V0(Aixc11)dV0Xilknc11的 薄 层 单 元(圆 形 区 域),部分代表整个结构剩下部分的非薄层单元。可以看出,螺栓连接机匣法兰结构的轴向刚度与轴向弹性模量、单元形函数和节点的轴向形变有关。对于一个具体的结构,横截面尺寸和单元类型是确定的,因此,当节点坐标系为笛卡儿坐标时,螺栓连接机匣法兰结构的轴向刚度由不同区域薄层单元的轴向弹性模量确定。1.2考虑止口接触面的非线性刚度薄层单元建模与模拟在航空发动机转子与静子机匣结构中,连接件与被连接件之间一般都
19、存在止口,如图 2 所示,主要起到增加连接件与被连接件之间的径向紧度和密封等作用,在一定程度上还能起到定位的作用。根据实际机匣结构,基于接触单元建立实体模型如图 3 所示。设置为轴向方向一端施加约束,另一端施加轴向载荷 02kN 的边界条件,计算得到非线性刚度特征如图 4 所示。由于接触单元网格复杂、计算缓慢,因此利用薄层单元法对螺栓连接结构进行简化,去掉螺栓和螺栓孔,用上下法兰接触面间的薄层单元模拟螺栓连接,采用壳单元对圆筒部件进行等效简化来等效模拟螺栓连接结构力学行为,在薄层单元与对应接触的壳单元之间建立绑定接触,螺栓的质量采用 Mass21单元集中模拟,建立的简化模型如图 5 所示。由于
20、止口过盈量导致摩擦力存在,螺栓连接结构的轴向刚度会表现出一定的非线性特征,如图 6 所示。在 ANSYS 软件中使用一种双线性的随动强化模型,通过对关键项进行设置,可以实现非线性刚度模拟。对于大多数受应变载荷的金属而言,遵循随动强化规律。随动强化满足包辛格效应,即刚度屈服前后,体现出明显的非线性。本文对薄层单元的材料赋予双线性的随动强化模型,分别针对屈服前后的刚度进行模拟,通过设置屈服前后的应力和弹性模量模拟螺栓连接结构的非线性轴向刚度特征,得到如图 7 所示的非线性特征。图2带有止口的螺栓连接结构Fig.2Boltedstructurewithseam(a)圆筒结构(b)局部放大图3螺栓连接
21、接触单元圆筒结构有限元模型Fig.3Finiteelementmodelofcylindricalstructurewithboltedcontactelement2.42.01.61.20.80.405101520轴向载荷/kN轴向形变/105 m图4接触单元非线性刚度特征Fig.4Nonlinearstiffnesscharacteristicsofcontactelements图5螺栓连接薄层单元圆筒结构有限元模型Fig.5Finiteelementmodelofcylinderstructurewithboltedthinlayerelement第6期刘玉等:航空发动机螺栓连接结构薄层
22、单元静力学非线性模型修正1483由图 7 可知,在载荷较小阶段,接触面处于黏滞状态,只发生弹性形变。随着横向载荷的增大,止口接触面出现滑移。在滑移阶段止口上、下两侧的接触状态发生改变使得接触刚度减小,轴向刚度减小,发生刚度屈服。利用接触单元和薄层单元都可以对螺栓连接结构的非线性刚度特征进行模拟,但是两种不同模拟方法的预测结果之间具有一定的差异性,与真实结果的准确性无法保证,因此需要针对于薄层单元模型进行模型修正。2非线性刚度薄层单元模型修正2.1双加权目标函数构建以静力学实测数据为基准的模型修正方法假定结构为线弹性及小形变,其修正的基本原理是:建立结构初始有限元模型后,选择合理的模型修正参数,
23、基于静力实测数据,通过优化迭代计算来不断调整有限元模型的参数,使得仿真模型响应计算值与结构实测响应值差异最小,进而实现模型参数修正。这个过程可转化为数学优化求解问题,根据不同的实际工程求解问题,构造相应形式的目标函数,如式(10)所示:f=fSFEM(P,c)Stest(P)(10)SFEM(P,c)式中是有限元仿真模型计算的响应值;Stest(P)Pc是试验数据的响应值;是施加的轴向载荷;是对应不同阶段的轴向弹性模量。本文将响应面方法引入到非线性模型修正中,并采用二次不含交叉项的双加权响应面多项式模型来建立响应面模型,对屈服前后分别构建具体的目标函数如下:J1(c1)=a1(c1)a1ta1
24、t2(11)J2(c2)=a2(c2)a2ta2t2(12)J1J2c1c2a1(c1)a2(c2)a1ta2tm1m2式中是发生屈服前模型预测轴向形变和试验轴向形变的残差;是发生屈服后模型预测轴向形变和试验轴向形变的残差;是发生屈服前待修正的薄层单元弹性模量;是发生屈服后待修正的薄层单元弹性模量;是发生屈服前的有限元轴向形变;是发生屈服后的有限元轴向形变;是发生屈服前的试验轴向形变;是发生屈服后的试验轴向形变。由式(11)式(12)可知,屈服前后薄层单元的弹性模量和螺栓连接结构发生刚度屈服前后的不同刚度特征之间存在着两个阶段的对应关系。但是对于薄层单元的弹性模量进行修正只能反映出屈服前后的刚
25、度特征,无法体现出发生屈服时的轴向载荷的临界值,因此,本文引入加权系数、,构造出一种双加权的响应面目标函数,如式(13)所示:J(c)=a1J1(c1)+a2J2(c2)(13)a1a2将修正后不同阶段的弹性模量带入式中,以发生屈服时轴向载荷的临界值为新的修正目标对加权系数、进一步修正。J1J2本文对屈服前后的目标函数、进行整合,构造双加权响应面目标函数,具体为J(c)=m1a1(c1)a1ta1t2+m2a2(c2)a2ta2t2(14)m1m2式中、为双加权响应面系数。2.2优化算法及误差指标E目标函数建立后,非线性有限元模型修正问题就转化为结构参数的优化问题。本文采用软件 Matlab
- 配套讲稿:
如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。