发明名称 一种基于细节疲劳额定值的高机动飞机寿命分析方法
摘要 本发明属于飞机疲劳学领域,涉及一种基于细节疲劳额定值的高机动飞机寿命分析方法,其特征在于,包括如下步骤:第一,编制高机动飞机DFR分析的载荷谱,第二,DFR基准值分析,第三,计算DFR<sub>0</sub>值,第四,计算构件疲劳额定值系数R<sub>C</sub>,第五,确定结构DFR许用值[DFR],第六,许用应力[σ<sub>max</sub>]计算,第七,疲劳强度裕度计算与评估,第八,结构可靠性寿命计算,第九,建立应力水平控制曲线,给出目标寿命对应的应力许用值。本发明的优点是:应用DFR法进行疲劳设计和寿命评估的过程比较简单,计算工作量也相对较少,能大幅度地提高整体部门的工作效率,对缩短新机研制周期、节约研制成本具有重要意义,经济效益显著。
申请公布号 CN104809321A 申请公布日期 2015.07.29
申请号 CN201310744491.6 申请日期 2013.12.31
申请人 中国航空工业集团公司沈阳飞机设计研究所 发明人 隋福成;周丽君;陈亮
分类号 G06F19/00(2011.01)I;G01N3/00(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京慕达星云知识产权代理事务所(特殊普通合伙) 11465 代理人 高原
主权项 一种基于细节疲劳额定值的高机动飞机寿命分析方法,其特征在于,包括如下步骤:第一,编制高机动飞机DFR分析的载荷谱针对机型的疲劳载荷谱,通过雨流计数法抽取完整的载荷(应力)循环。按当量载荷循环的峰值进行排序,初步选取峰值载荷,使疲劳载荷谱中在选取的峰值载荷上下的载荷循环对应的损伤相当。按照当量等幅载荷谱与疲劳载荷谱下达到疲劳寿命时有相同临界损伤的原则,导出疲劳载荷谱单位飞行小时对应的当量等幅载荷循环数,用过载表示的公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>N</mi><mi>fh</mi></msub><mo>=</mo><mfrac><mrow><mfrac><mn>1</mn><msub><mi>T</mi><mn>0</mn></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>yi</mi></msub><mo>)</mo></mrow><mi>dl</mi><mi>m</mi></msubsup></mrow><msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>df</mi></msub><mo>)</mo></mrow><mi>dl</mi><mi>m</mi></msubsup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000011.GIF" wi="1290" he="229" /></maths>当疲劳载荷谱为应力谱时,公式中的过载用应力取代。根据等损伤原则将选取的峰值载荷(对应脉动循环)按应力比R=0.1进行等损伤折算,最终确定当量等幅载荷循环的峰值。第二,DFR基准值分析选取飞机结构典型结构部位制成试验件,按上述载荷谱进行疲劳试验,并记录其均值寿命N<sub>50</sub>。高机动飞机结构的可靠性指标为99.9%的可靠度和90%的置信度,求出与其对应的试验可靠性寿命:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mn>99.9</mn><mo>/</mo><mn>90</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>N</mi><mn>50</mn></msub><mrow><msub><mi>S</mi><mi>R</mi></msub><msub><mi>S</mi><mi>C</mi></msub><msub><mi>S</mi><mi>T</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000012.GIF" wi="1266" he="137" /></maths>式中:S<sub>R</sub>‑可靠系数S<sub>C</sub>‑置信系数S<sub>T</sub>‑试件系数DFR基准值按下式计算:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>DFR</mi><mi>base</mi></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>R</mi><mo>)</mo></mrow><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><msub><mi>&sigma;</mi><mi>max</mi></msub></mrow><mrow><mn>0.9</mn><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mn>0.1</mn><mo>-</mo><mi>R</mi><mo>)</mo></mrow><msub><mi>&sigma;</mi><mi>max</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>N</mi><mrow><mn>99.9</mn><mo>/</mo><mn>90</mn></mrow></msub><msub><mi>N</mi><mn>50</mn></msub></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>m</mi></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000021.GIF" wi="1505" he="183" /></maths>式中:σ<sub>m0</sub>‑等寿命曲线参数,对于铝合金、钛合金、中强钢(σ<sub>b</sub>≤1600MPa)和高强钢(σ<sub>b</sub>&gt;1600MPa),σ<sub>m0</sub>分别取为310、620、930和1240MPa。σ<sub>max</sub>,R‑试验的最大应力与应力比。m‑S‑N曲线的斜率,一般取4.0。第三,计算DFR<sub>0</sub>值DFR<sub>0</sub>是实际结构的基本DFR值,通过各种修正系数以及应力集中的修正与典型结构的DFR基准值(DFR<sub>base</sub>)建立联系,其关系式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>DFR</mi><mn>0</mn></msub><mo>=</mo><msub><mi>DFR</mi><mi>base</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>A</mi><mo>&CenterDot;</mo><mi>B</mi><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><msub><mrow><mo>(</mo><msub><mi>K</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>base</mi></msub><msub><mi>K</mi><mi>t</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000022.GIF" wi="1416" he="153" /></maths>式中A、B……即为修正系数;(K<sub>t</sub>)<sub>base</sub>为典型结构细节处的应力集中系数,K<sub>t</sub>为实际结构细节处的应力集中系数。第四,计算构件疲劳额定值系数R<sub>C</sub>R<sub>c</sub>定义为实际结构(当量严重细节数为n)的细节疲劳额定值DFR<sub>n</sub>与基本结构(单细节)的细节疲劳额定值DFR<sub>1</sub>的比值,其公式如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>R</mi><mi>c</mi></msub><mo>=</mo><msup><mn>10</mn><mrow><mfrac><mrow><mo>[</mo><msup><mi>&phi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mn>0.999</mn><mo>)</mo></mrow><mo>-</mo><msup><mi>&phi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msup><mn>0.999</mn><mrow><mn>1</mn><mo>/</mo><mi>n</mi></mrow></msup><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mn>0</mn></msub></mrow><mi>m</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></msup></mrow>]]></math><img file="FSA0000099694700000023.GIF" wi="1264" he="131" /></maths>式中:σ<sub>0</sub>为对数寿命标准差,取为0.14;φ<sup>‑1</sup>()为标准正态分布函数的逆函数。第五,确定结构DFR许用值[DFR][DFR]按下式计算:[DFR]=DFR<sub>0</sub>·R<sub>c</sub>    (6)第六,许用应力[σ<sub>max</sub>]计算[σ<sub>max</sub>]按下式计算:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>[</mo><msub><mi>&sigma;</mi><mi>max</mi></msub><mo>]</mo><mo>=</mo><mfrac><mrow><mn>0.9</mn><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&CenterDot;</mo><mo>[</mo><mi>DFR</mi><mo>]</mo></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>R</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>N</mi><mi>e</mi></msub><mrow><mn>5</mn><mo>&times;</mo><msup><mn>10</mn><mn>4</mn></msup></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>m</mi></mfrac></msup><mo>+</mo><mrow><mo>(</mo><mi>R</mi><mo>-</mo><mn>0.1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>[</mo><mi>DFR</mi><mo>]</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000031.GIF" wi="1306" he="221" /></maths>式中:N<sub>e</sub>‑目标寿命。第七,疲劳强度裕度计算与评估根据上述计算结果可以得出疲劳强度裕度为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfrac><mrow><mo>[</mo><msub><mi>&sigma;</mi><mi>max</mi></msub><mo>]</mo></mrow><msub><mi>&sigma;</mi><mi>max</mi></msub></mfrac><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000032.GIF" wi="1230" he="151" /></maths>疲劳裕度为负,意指必须减小使用应力,同时(或者)通过改进细节设计,增加其细节疲劳额定值。要使在给定使用应力下满足设计目标寿命,改进设计的指标是使要求的DFR所对应的疲劳裕度等于或大于零。第八,结构可靠性寿命计算可靠性寿命的计算公式如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>N</mi><mo>=</mo><msup><mrow><mo>{</mo><mfrac><mrow><mo>[</mo><mi>DFR</mi><mo>]</mo></mrow><mi>Y</mi></mfrac><mo>}</mo></mrow><mi>m</mi></msup><mo>&CenterDot;</mo><mfrac><mrow><mn>5</mn><mo>&times;</mo><msup><mn>10</mn><mn>4</mn></msup></mrow><mi>FRF</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000033.GIF" wi="1304" he="162" /></maths>式中,FRF‑可靠性系数,Y是任意应力比下的载荷循环通过等寿命线折算到指定应力比0.1下的最大应力的转换公式,其表达式如下:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>R</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>max</mi></msub></mrow><mrow><mn>0.9</mn><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mn>0.1</mn><mo>-</mo><mi>R</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>max</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000099694700000034.GIF" wi="1385" he="154" /></maths>第九,建立应力水平控制曲线,给出目标寿命对应的应力许用值应力水平控制曲线即特征应力(1g应力或单位载荷应力)‑可靠性寿命曲线。由多个特征应力‑可靠性寿命数据点,即可建立应力水平控制曲线。由应力水平控制曲线即可给出目标寿命对应的应力许用值。
地址 110035 辽宁省沈阳皇姑区塔湾街40号
您可能感兴趣的专利