2009年6月第26卷第3期沈阳航卒工业学院学报JournalofShenyangInstituteofAeronauticalEngineeringJun.2009V01.26No.3文章编号:1007—1385(2009)03—0010—04基于雨流计数法的随机声疲劳寿命估算方法研究郭小鹏沙云东张军(沈阳航空工业学院动力与能源T:程学院,辽宁沈阳110136)摘要:针对飞机结构中的声疲劳破坏问题,以雨流计数法的基本原理和计数规则为基础,研究了雨流计数法在Matlab语言中的具体实现方法并得到载荷谱,结合材料性能S—N曲线在Miner线性累积损伤理论准则下估算结构疲劳寿命。并以某航空薄壁柱壳结构为例进行疲劳寿命估算,取得了良好的效果,从而为结构的声疲劳寿命估算提供了一种可行的方法。关键词:声疲劳;疲劳寿命;雨流计数法;薄壁结构中图分类号:V215.5;V232.6文献标识码:A声疲劳是结构在高频低幅值随机载荷下,因结构中快速交变应力的作用而使结构产生声疲劳破坏的现象,这种快速交变应力是因为噪声激励引起结构振动而产生的…。随机载荷计数法有十多种,其中最常用的有雨流法、幅度计数法和峰值计数法口。3J。航空发动机噪声及相关声疲劳寿命问题属于航空结构动强度分析和寿命预测中的学科前沿问题。本文通过雨流计数法对某航空薄(3)当雨流遇到来自上面屋顶流下的雨流时,即行停止,取出所有的全循环,并记录下各自的幅值和均值。壁柱壳结构进行疲劳寿命估算,其结果将对我国航空发动机结构完整性设计、提高结构动强度和使用寿命提供重要参考依据。1雨流计数法基本原理雨流计数法是以双参数法为基础的一种计数葺F坞;蕊2r・‘t皇>-7O拶。X—I{4:上』’吣=≥枷:9ll法,考虑了动强度(幅值)和静强度(均值)两个变量,符合疲劳载荷固有的特性,其主要功能是把实测载荷历程简化为若干个载荷循环,供疲劳寿命估算和编制疲劳试验载荷谱使用HJ。萝1{3‰≮逡≤:ll’▲●’17L●f地常用的雨流计数法的计数规则为:(1)重新安排应力一时间历程,以最高峰值或最低谷值(视两者的绝对值哪一个更大)为起点。(2)雨流依次从每个峰(或谷)的内侧向下流,在下一个峰(或谷)处落下,直到对面有一个图1雨流计数法计数原理雨流计数法的基本原理如图1所示,第一个雨流自0点处第一个谷的内侧流下,从1点落至1’后流至5,然后下落。第二个雨流从峰1点内侧流至2点落下,由于1点的峰值低于5点的峰值,故停止。第三个雨流自谷2点的内侧流到3,自3点落下至37,流到1’处碰上上面屋顶流下的雨流而停止。如此下去,可以得到如下的计数循环块:比起点更高的峰值(或更低的谷值)停止。收稿日期:2008—12一04基金项目:航空基础科学基金(项目编号02C54007)作者简介:郭小鹏(1982一),男,河北张家口人,在读研究牛,主要研究方向:航空发动机强度振动及噪声。E—mail:guoxiaopen9312@163.corn;沙云东(1966一),男,黑龙江阿城人。教授。3—4—3’、1—2—17、6—7—67、8—9—8’、11—12—1l、13—14一137和12—15一12’。2雨流计数法在程序中具体实现过程2.1数据压缩万方数据第2期郭小鹏等:基于雨流计数法的随机疲劳寿命估箅方法研究数据压缩就是把原始的计算数据处理成便于计取循环数的数组。为避免在峰谷值提取时产生影响,应先将相邻的等值数压缩(即把相同数值的保留一个),然后再提取峰谷值。如图2所示,当遇到相等数时,取的是其中的第一个数值,其余的点则不予考虑。在峰谷值检测时,如图3所示,若为真值时就说明该点同时大于(小于)两边的数字,则可知该点就是峰(谷)值。否则,则认为该点不是峰值也不是谷值,可以不予考虑。图2雨流计数法等值数压缩检测流程图图3雨流计数法峰谷值检测的流程图2.2循环块的提取(1)一次雨流计数一次雨流计数是从压缩处理过的数据中提取万方数据循环,并记录其特性值,如峰值、谷值、幅值以及均值等。如图4所示,在数组中能提取出循环的形式有两种,若满a大于等于b并且b小于C时就可以认为提取一个循环(图中构成三角形的部分),同时去除e(i)和e(i+1)两点。这样对这个数组反复判断就可以取出所有的循环和记录所剩的各点。这种判断适用于四个点组成的各种情形,提取循环的流程图参见图5所示。,,●,图4雨流计数法循环的判断图图5雨流计数法提取循环的流程图(2)对接一次雨流计数剩下的点构成的波形是一个标准的发散收敛型波形,如图6所示。这时按雨流计数法则无法再形成整循环,只能将其在绝对值最大的点处截开再进行首尾对接。(3)二次雨流计数二次雨流计数是将完成对接的波形继续提取循环直到剩下三个点(即是数组中最值构成的整循环)为止,程序中实现也很简单,只需将对接完成的数组放入一次雨流计数中就行了。(4)应力谱雨流计数的结果可以用一张均、幅值矩阵来12沈阳航空工业学院学报第26卷裁荷时问图6发散一收敛波表示,这就是载荷谱[5】。对薄壁柱壳构件的三个测点测得的应力时间历程结果进行雨流计数,将得到的计数结果填人载荷谱表,就可以得到应力谱。3Miner线性累计损伤理论Miner线性累积损伤理论是工程上广泛采用的一种疲劳寿命计算方法"]。如果是单级加载,循环比等于l时即出现破坏。如果是多级加载,总损伤等于各循环比(或损伤比)的总和,且当循环比总和等于1时发生破坏。在多个应力水平及随机载荷作用下的损伤计算公式D=兰主袭(,)i=1J=1”口式中:nii表示构件在第i级应力幅值和歹应力均值载荷下的工作循环次数由前面提取出的均值、幅值载荷谱得到;Ⅳii表示在存活率为P第i级应力幅值和_『应力均值载荷下构件的疲劳寿命值,数据来源于S一Ⅳ曲线。根据Miner线性疲劳累积损伤理论的破坏准则:ⅣrD=1(2)式中:Ⅳ,——所测得的载荷谱作用的次数。根据式(2),可求得薄壳的寿命1T=Nr・At=去・At(3)上,式中:缸——所测得的载荷谱作用的时间。根据式3,再加上所需的材料性能Js一Ⅳ曲线,就可以进行疲劳寿命的估算了一】。应用线性累积损伤理论来计算疲劳寿命时,主要是将载荷谱中各级载荷循环作用的顺序和次数用雨流计数法筛选出来,各级载荷之间的迟滞效应则不予考虑。万方数据4应用实例4.1声激励载荷有限带宽高斯白噪声在其频带宽度范围内,各频率分量是均匀分布的。薄壁柱壳结构在有限带宽高斯白噪声载荷作用下的响应更能反映出结构声疲劳问题的本质,同时还可以比较同一结构相同频带宽度范围内不同声压级噪声作用下的响应问题。G。∽=4/Af×10‘奇-1∞(4)式中:G。(.厂)表示有限带宽高斯白噪声功率谱密度;SPL表示带宽声压级;Af表示频带宽度。取频带宽度为0—1024Hz,频率分辨率为8Hz,根据公式4计算带宽声压级分别为120dB、140dB、160dB的有限带宽高斯白噪声功率谱密度,数据如表1所示。表l有限带宽高斯白噪声率谱密度4.2算例分析以开孑L圆柱壳为算例模型,如图7所示。通过Ansys有限元软件的计算,获得薄壁柱壳上三个测点的应力一时间历程。节点A(自由边上)、B(周向中点)和c(固定端)三点分别在噪声载荷120dB、140dB、160dB下应力均方值分别如表2、表3和表4所示。对比三张表可知在开孔边处的应力值是所有应力之中最大的,这说明孑L边容易}}{现应力集中;同样说明了声压级每增加20dB,结构的均方应力值相应增加约两个数量级。图7开孔圆柱壳节位置示意图节点B在120dB、140dB、160dB下的应力响应功率谱密度如图8所示。结合模态频率值,从应力频谱图中发现,在中低频带内有多个模态被激起了响应峰值,声压级每增加20dB,应力功率谱密度就增加的两个数量级,这与前面均方应力第2期郭小鹏等:基于雨流计数法的随机疲劳寿命估箅方法研究13罢‘窑舄I瑟皇峥童真避图8节点B处的VonMises应力相应增加的数量级是一致的。在获得应力一时间历程的前提下,应用雨流计数法计数,结合S—N曲线和疲劳累积损伤理论分别在噪声载荷120dB、140dB、160dB下对点A、点B和点c进行疲劳寿命寿命估算,其结果如表5所示。表2声压级120dB时圆柱壳结构应力均方值(Pa2)表3声压级140dB时圆柱壳结构应力均方值(Pa2)表4声压级160dB时圆柱壳结构应力均方值(Pa2)表5开孔圆柱壳疲劳寿命估算值万方数据通过开孔圆柱壳疲劳寿命估算值可知在同一应力水平下。节点A、B和C的寿命估算值排列顺序从大到小依次是:c、A、B。即自由端的疲劳寿命最长,固定端次之,周向中点的疲劳寿命最短。并且声压级每增大lOdB寿命值就降低约一个数量级。5结论(1)本文论述了雨流计数法基本原理和计数规则,雨流法可有效地进行数据压缩,加快疲劳寿命估算速度,从而更好地实现构件的疲劳寿命预测的实时陛,它的计数结果具有较高的准确性,由此进行的疲劳寿命计数比较接近实际;(2)通过开孔圆柱壳算例可知雨流计数法估算疲劳寿命方法是可行的,其结果对于航空航天飞行器疲劳寿命估算具有重要的理论价值和工程应用价值;(3)该方法对航空发动机燃烧室火焰简结构完整性设计、提高结构使用寿命具有一定的工程应用价值。参考文献:[1]金奕山,李琳.关于航窄发动机结构声疲劳寿命估算方法的探讨[J].航空动力学报,2003,6(3):374—377.[2]姚卫星.结构振动疲劳寿命估算方法研究[D].南京:南京航空航天大学,2005.[3]阎楚良,王公权.峰值计数法疲劳载荷谱的编制与统计处理程序[J].农业机械学报,1984,3(1):51—57.[4]王宏伟,刑波。骆红云.雨流计数法及其在疲劳寿命估算中的应用[J].矿山机械,2006,3(34):95—9r7.[5]刘惟信.机械可靠性设计[M].北京:清华大学出版社,2006.[6]A.Fmemi,L.Yang.Cumulativefatiguedamageandlifepredie—tiontheories:asurveyofthestateoftheartforhomogeneousma-terials[J].InternationaljoumaloffatigIle,1998,20.9—34.[7]王东锋,汪定江,王新坤.构件安全疲劳寿命估算中的P—s—N曲线修正[A].第二届中国航空学会青年科技论坛文集[C]。2006.[8]张军,沙云东,倪绍华.基于功率谱密度的结构声疲劳寿命估算方法研究[J].沈阳航搴工业学院学报,2008,25(1):1l—14.(下转第9页)第3期苏春峰等:接触非线性仿真中接触刚度因子选取的方法研究9能低估了接触刚度的值,如果总体的收敛困难是究[J].北京航空航天大学学报。1998,13(4):l一5.由于不平衡力和位移增量达到收敛值需要过多的【4]白润波,曹平周,曹茂森,陈建锋.基于优化一反分析法的接触迭代次数,而不是由于过大的渗透量,那么接触刚刚度因子的确定[J].建筑科学,2008,24(1):46—47.[5]F.J.uu,Analysisfailureonattenuatingvibrationofsqueeze度的值可能被高估;filmdamper[J],JournalofBeijingUniversityofAeronauticsand(3)按需要调查接触刚度值,重新分析,直至Astronautics.1998,24(3):259—300结果达到要求的收敛精度。[6】T.B.Yu,z.JiaIIg,Y.D.Gong,etal,Experimentofaqueezefilmdamperfordeepholemachining[J].JournalofNortheastern参考文献:university(NaturalScience),2002,23(6):577—580.[7]娄小宝,沙云东,苏春峰.航空发动机燃烧室噪声产生机理及[1]G.M.L.Gladwell,经典弹性问题理论中的接触问题[M].北其主要影响因素[J].沈阳航空T业学院学报,2008,25(5):京:北京理工大学H{版社.1991.1—2:118一155.28—31.[2]曹磊,高德平,江和甫.弹性环式挤压油膜阻尼器设计因素研[8]王英,沙云东.航空发动机故障诊断技术综述[J].沈阳航空究[J].燃气涡轮试验与研究,2006:19(4):l一5.工业学院学报,2007,24(2):11—14.[3]周明,李其汉,晏砺堂.弹性环式挤压油膜阻尼器减振机理研Studyoncontactstiffnessfactorselectionmethodincontactnon—.1inearsimulationSUChun—-fengAIYan—-tingLOUXiao——bao(ShenyangInstituteofAeronauticalEngineering,LiaoningKeyLaboratoryofDigitalProcessandTesting,Ha-oningShenyang110136)Abstract:Itisdifficultandcrucialtogetthecontactstiffnessvalueincontactanalysis.And,thereisnosatis—factoryselectingmethodforcontactstiffnessselectionatpresent.ThefiniteelementsoftwareANSYSiSusedtoanalyzetheoreticallyelasticringsqueezeoilfilmdamper(SFD)contactbehavior,includingitscontactdynam—ics、conditionsandparameters,especiallythecontactstiffnessfactorsselection.Thecorrespondingcalcula—tionresultsofdifferentcontactstiffnessvaluesarelistedandtheconvergenceiSanalyzed。andthenumericalrecursivemethodofcontactstiffnessfactorisputforward.WhichcanbeappliedtotheFEcalculationofsimi—larcontactcase.Keywords:damper;contactstiffness;numericalrecursivemethod(上接第13页)Randomfatiguelifeestimationbasedonrain・--flowcountingmethodGUOXiao—pengSHAYun—dongZHANG—jan(AircraftPowerandEnergyCollege,ShenyangInstituteofAeronauticalEngineering,LiaoningShenyang110136)Abstract:Interviewofairplanestructuresoundfatiguefailurequestion,basedontheprincipleandregulationofrain—flowcountingmethod,thispaperestablishesanelaboratedprogrammewithmatlab,takeacertainairofthin—walledcylindricalshellstructureasanexample,estimatingthefatiguelifebyloadingspectrumandS—Ncurveanalysing.Theresultsprovethatthemethodhasasuperiorskillforfatiguelifeestimating.Keywords:soundfatigue;fatiguelife;rain—flowcountingmethod万方数据基于雨流计数法的随机声疲劳寿命估算方法研究
作者:作者单位:刊名:英文刊名:年,卷(期):引用次数:
郭小鹏, 沙云东, 张军, GUO Xiao-peng, SHA Yun-dong, ZHANG jun沈阳航空工业学院动力与能源工程学院,辽宁,沈阳,110136沈阳航空工业学院学报
JOURNAL OF SHENYANG INSTITUTE OF AERONAUTICAL ENGINEERING2009,26(3)0次
参考文献(8条)
1.金奕山.李琳 关于航窄发动机结构声疲劳寿命估算方法的探讨 2003(3)2.姚卫星 结构振动疲劳寿命估算方法研究 2005
3.阎楚良.王公权 峰值计数法疲劳载荷谱的编制与统计处理程序 1984(1)
4.王宏伟.邢波.骆红云 雨流计数法及其在疲劳寿命估算中的应用[期刊论文]-矿山机械 2006(3)5.刘惟信 机械可靠性设计 2006
6.A Faterni.L Y ang Cumulative fatigue damage and life prediction theories:a survey of the state ofthe art for homogeneous materials 1998
7.王东锋.汪定江.王新坤 构件安全疲劳寿命估算中的p-S-N曲线修正 2006
8.张军.沙云东.倪绍华 基于功率谱密度的结构声疲劳寿命估算方法研究[期刊论文]-沈阳航空工业学院学报2008(1)
相似文献(10条)
1.期刊论文 杜三虎.黎在良.杨嘉陵 复合材料结构声疲劳的一种数值模拟方法 -北京航空航天大学学报2003,29(1)
复合材料结构声疲劳问题是许多工程技术人员关心的课题,经常采用数值方法对其进行分析.利用边界元法对其动力响应情况进行了数值模拟,并发展了一种可以进行复合材料结构声疲劳完整数值模拟的方法.所得结果可以为预测复合材料结构的声疲劳寿命提供参考.
2.学位论文 杜三虎 复合材料工程结构体声疲劳的边界元数值模拟 2001
该文以现代飞行器设计中涉及的结构声疲劳损伤破坏问题为背景,利用边界元方法(BoundaryElementMethod(BEM))对纤维增强类复合材料层合结构动力响应及声疲劳问题进行了数值模拟.理论和计算的模型侧重于二维平面应变及三维情形的BEM公式推导及其编程的实现,为了考察理论模型的可靠性,还进行了复合材料结构的振动实验.
3.期刊论文 张军.沙云东.倪绍华.ZHANG Jun.SHA Yun-dong.NI Shao-hua 基于功率谱密度的结构声疲劳寿命估算方法研究 -沈阳航空工业学院学报2008,25(1)
针对航空结构的声疲劳破坏问题,提出了结构声疲劳寿命估算的两种方法:时域统计法和功率谱密度法,其主要采用功率谱密度方法,应用振动理论、强度理论、Miner线性累积损伤理论和S-N曲线关系式等相关理论,建立了结构受窄带和宽带随机噪声载荷作用下的疲劳寿命估算公式,为工程上估算结构疲劳寿命提供了一定的参考价值.
4.会议论文 金奕山.李琳 航空发动机结构声疲劳寿命估算方法的研究 2002
航空发动机是一种功率强大且分布复杂的噪声源,处于这种宽频带高能级声激励环境中的发动机薄壁机匣结构极易发生疲劳破坏,其声疲劳分析在发动机强度设计中占有很重要的地位.有关发动机机匣结构声载荷的预测及声疲劳应力的估算人们已经做了大量研究,但到目前为止,还没有一种合适的方法用来进行此种结构的疲劳寿命估算.本文介绍了一种由疲劳应力功率谱密度函数直接计算雨流循环概率密度的理论方法,该法在估算发动机机匣结构的声疲劳寿命中将起到很大作用.
5.学位论文 孙龙生 飞机结构声疲劳寿命预测及分析平台的构筑 2003
该文针对现阶段声疲劳试验的不足,以及工程中常用的诺谟图方法的局限,整理了目前比较成熟的典型蒙皮壁板声疲劳寿命研究成果,提出设计并编制了一套功能强、易于使用,可扩展性好的声疲劳寿命预测及分析平台.通过对声疲劳问题研究的回顾,根据典型飞机蒙皮壁板模型,整理了四边固支,四边简支两种边界条件下多阶模态固有频率,非平稳及平稳随机白噪声载荷下的位移和应力响应,以及如何从等幅S-N曲线求得随机载荷应力下的疲劳寿命.在该文的后半部分,为了代替目前常用的列线图求解声疲劳响应方法,着重运用面向对象程序设计技术,完成声疲劳寿命预测及分析平台的设计和实现.通过已有的算例进行验证,结果表明计算迅速可靠,可以代替以列线图为主要手段的求解方法,广泛的应用于工程领域.
6.学位论文 葛森 飞机壁板高温声疲劳特性分析 1999
飞机壁板结构的高温声疲劳研究是声疲劳领域内的一项前沿性研究课题,也是航空工程领域中所面临的急待解决的问题。所谓高温声疲劳是指温度与噪声载荷联合作用于飞机结构而且量值都大到足以在结构中产生疲劳损伤或疲劳破坏。与常温下的声疲劳不同,高温声疲劳分析时所采用的是考虑温度影响的运动方程,由于横向位移的量值与板厚度尺寸相比不现进一个小量,根据大挠度非线性理论,该运动方程为高阶的偏微分方程。该文采用了成熟的等价线性化方法,将该高阶偏微分方程的求解问题化为二同分方程的求解。为了便于理解各种影响参数的物理意义,该文利用解析分析的办法,采用考虑基频响应为主的单自由度模型,探讨出一整套完整的求解分析的方法与步骤,最终求得了结构内各种体现声疲劳特性的响应统计特征。所得到的结果是对飞机壁板结构进行高温声疲劳分析、计算或设计的重要结果。一旦有条件,能够通过试验获得材料与结构的高温与强噪声联合作用条件下的基础数据曲线,用所得结构便可以直接查图表曲线,得到疲劳寿命的有关数据。为了验证这些结构的可入性与不关问献中的一些试验结果进行了对比,证明了结果是合理的,并通过比较明确了今后进一步深入研究方向。
7.期刊论文 金奕山.李琳 关于航空发动机结构声疲劳寿命估算方法的探讨 -航空动力学报2003,18(3)
噪声载荷是一种高频、宽带随机载荷,在低循环疲劳寿命预测中被广泛采用的时域计数法很难直接利用.根据现有的累积疲劳损伤理论,噪声载荷下结构疲劳寿命分析的关键在于如何在频域中获得雨流循环的概率密度函数.本文对从疲劳应力功率谱密度出发获得寿命估算所需要的应力峰值或应力雨流循环幅值的概率密度函数的方法进行了分析与比较,并对发动机结构声疲劳分析中值得深入研究及发展的方法给出了建议.
8.会议论文 张秀义 歼八Ⅱ型飞机进气道结构声疲劳研究 19949.学位论文 屈伸 航空薄壁板结构在噪声载荷作用下响应的研究 2006
本文结合航空科学基金课题(课题号:02C54007)完成,主要研究飞行器薄壁结构在高声级宽频带随机噪声载荷作用下的响应问题,这是研究结构声疲劳问题的基础。声疲劳问题的研究主要包括三方面:1.建立随机噪声载荷模型;2.计算结构在随机噪声载荷作用下的响应;3.预测结构的疲劳寿命。本文重点研究结构在随机噪声载荷作用下位移响应和应力响应的计算方法。在基本理论研究基础上,建立了结构响应的计算模型。选取具有固支边界和简支边界的矩形金属薄壁板单元作为研究对象,在时域和频域中研究该薄壁板在噪声载荷作用下响应的计算方法,并完成了位移响应和应力响应的计算。 1.采用Miles单自由度理论,假定噪声载荷作用下结构响应以基频为主,选取试验测得的飞机地面开车喷流噪声载荷作为输入,计算四边固支和四边简支矩形薄壁板上各点的位移响应、速度响应、加速度响应和应力响应的均方值,并将计算结果与试验数据进行比较。计算结果给出了该矩形薄壁板上各点位移响应均方值和应力响应均方值的分布,应力响应最大的位置是结构上容易出现疲劳破坏的位置。 2.根据弹性力学中薄板的振动理论,计算出结构的各阶振动模态。从薄板的振动微分方程出发,引入结合受纳函数,以有限带宽高斯白噪声载荷和试验测得的噪声载荷作为输入,采用功率谱密度方法计算结构上各点位移响应的功率谱密度和应力响应的功率谱密度。该方法在频域揭示了声激励、结构固有特性和响应的关系,为薄壁板结构设计中避免出现声疲劳破坏提供理论依据。 3.根据弹性力学中的能量原理,采用变分法求解薄壁板结构在随机噪声载荷作用下的大挠度弯曲问题,计算出薄壁板结构中各点位移响应和应力响应的时间历程。Miles单自由度理论和功率谱密度方法均基于薄板的小挠度弯曲假设,考虑到在高声级随机噪声载荷作用下结构响应有时是大挠度的,而且,目前大多数疲劳分析都是在时域内进行的,本文采用变分法求解薄壁板结构在随机噪声载荷作用下的大挠度弯曲振动响应,具有创新性。 本文研究的航空薄壁板结构在随机噪声载荷作用下响应的计算方法和获得的计算结果对分析和解决航空薄壁板结构的声疲劳问题、提高结构动强度具有重要意义。
10.会议论文 邵汉林.马建忠 声疲劳实验中的声致振动响应分析 1997
该文介绍了在声疲劳试验方法研究中对钢板所做的声致振动实验,并根据实验结果对声致振动响应进行了分析和讨论,将实验装置的声学系统元件与试件的耦合作用归结为对试件非耦合时固有频率产生偏移的结果,按照非耦合时的情况进行近似处理,这将便于对试件的疲劳寿命进行分析和对实验采取合适的监控措施。
本文链接:http://d.g.wanfangdata.com.cn/Periodical_syhkgyxyxb200903003.aspx
下载时间:2010年3月8日