一种基于雨流计数循环的航空发动机载荷谱相关系数计算方法
【专利摘要】本发明公开了一种基于雨流计数循环的航空发动机载荷谱相关系数计算方法，对不同飞行剖面中各参数的载荷谱采用基于雨流计数法进行统计；滤除幅值小于幅值阈值的载荷循环后，获得各剖面各参数的循环数，同一飞行剖面不同的参数载荷谱得到不同的循环数，将不同参数在各个剖面统计得到的循环数一一对应，每个载荷参数都可得到一列循环数据，进而计算载荷谱各参数之间的相关系数。本发明去除了发动机转速与过载值间的延迟现象，结果与实际更相符。很大程度上代表机动循环的法向过载值与代表气动循环的转速的相关系数接近于1。进一步为航空发动机寿命预测与评估提供依据。
【专利说明】
一种基于雨流计数循环的航空发动机载荷谱相关系数计算 方法
技术领域
[0001] 本发明涉及一种基于雨流计数循环的航空发动机载荷谱相关系数计算方法，属于 航空发动机载荷谱技术领域。
【背景技术】
[0002] 日前我国已经对多种发动机飞行任务剖面进行了比较全面的空测和统计，获得了 大量有用的使用信息。如何利用这些飞行任务剖面编制合理的载荷谱，为发动机的应力分 析和寿命设计提供依据就成为十分突出的问题。以往国内在编制发动机飞行载荷谱时往往 只考虑一个飞行剖面参数的变化，即单参数谱的编制。单参数谱的编制方面国内已经做了 大量深入细致的研究工作。然而，大多数航空发动机构件承受的载荷是一个多分量的载荷 谱，如轴承、封严装置、涡轮盘的轮缘所承受的应力等，不仅与转速有关，而且与温度、压力 等因素有关。而安装节上的应力与发动机的推力和3个方向的飞行过载系数(轴向过载NX、 法向过载NY，侧向过载NZ)有关，是一个承受多参数载荷的零件。单参数编谱对于某些零件 可能十分有效，而对承受复杂载荷的零件则显得明显不足。航空发动机飞行任务剖面的处 理实际上是一个多参数谱问题.在进行编谱时必需考虑飞行任务剖面多参数的作用及其相 互间的影响。但是各参数之间的相关系数求解方法目前尚不成熟有待进一步改进。因此，研 究多参数载荷谱的相关关系对研究多参数载荷谱编制方法具有重要的理论意义和工程应 用价值。
[0003] 航空发动机的飞行载荷谱是随机的，如叶片、转子系统的高频振动载荷谱以及发 动机转速谱、飞行速度谱等等。其中发动机转速谱、发动机重心法向过载系数谱以及飞行速 度谱的谱型与飞机的操作使用密切相关，譬如飞机在做机动动作飞行时发动机转速剖面表 现为从最大到慢车之间的反复变化，形成典型的慢车-最大-慢车循环，发动机重心法向过 载系数谱则从1.0增大到某一极大值后再回到1.0,对应于一个个机动飞行动作。通过对与 操作相关的发动机载荷谱的分析发现(如附图4)，法向过载系数任务段的出现与发动机转 速变化具有比较大的相关性。
[0004] 但是现有的机动载荷与气动载荷的相关性分析是将各载荷参数当作一个完全随 机系列来计算相关系数，结合任务段分析发现，每一个发动机转速任务段的起始时间总是 先于法向过载任务段的起始时间，两者之间的时间差在几秒至几十秒内不等(如附图2)，发 动机转速任务段的峰值时间同样是先于法向过载任务段的峰值时间的（如附图3)，按照上 述方法计算出来的机动载荷与气动载荷的相关系数偏小。因此需要对现有的载荷谱参数相 关系数求解方法做进一步改进。
【发明内容】
[0005] 本发明针对上述问题的不足，提出一种基于雨流计数循环的航空发动机载荷谱相 关系数计算方法，该方法能够弥补现有的航空发动机载荷谱相关系数计算方法计算出的机 动载荷与气动载荷的相关系数偏小的问题，且去除了转速与过载值间的延迟现象。
[0006] 本发明为解决上述技术问题提出的技术方案是：
[0007] -种基于雨流计数循环的航空发动机载荷谱相关系数计算方法，对不同飞行剖面 中各参数的载荷谱采用基于雨流计数法进行统计。给出不同飞行剖面中各参数的幅值阈 值，滤除幅值小于幅值阈值的载荷循环后，获得各剖面各参数的循环数，同一飞行剖面不同 的参数载荷谱得到不同的循环数，将不同参数在各个剖面统计得到的循环数一一对应，每 个载荷参数都可得到一列循环数据，进而计算载荷谱各参数之间的相关系数。
[0008] 所述包括以下步骤：
[0009] 步骤1，获取航空发动机的飞行任务剖面的个数以及每个飞行任务剖面的飞行参 数。
[0010] 步骤2,对步骤1获取的每个飞行任务剖面的飞行参数进行飞行剖面峰谷值的提 取。
[0011] 步骤3,对步骤1中得到的每个飞行任务剖面的飞行参数根据步骤2得到的飞行剖 面峰谷值进行雨流计数法统计。
[0012] 步骤4，将各参数幅值小于[最小幅值+(最大幅值-最小幅值)的10% ]剔除，再将幅 值小于门槛值0.2的法向过载NY、纵向过载NX、侧向过载NZ剔除，获得删除小幅值后的各参 数幅值个数即载荷循环数，完成无效循环幅值的去除。无效幅值的去除可以在雨流计数后 进行，也可在雨流计数中进行。
[0013] 步骤5,根据步骤4去除的无效循环幅值的飞行参数获取各剖面的各个参数经过雨 流计数后的循环数，即幅值或均值的个数。
[0014] 步骤6,将步骤1获取的每个飞行任务剖面的飞行参数在步骤5获得的各剖面的各 个参数一一对应，进而计算载荷谱各参数之间的相关系数。
[0015] 优选的:所述步骤2中飞行剖面峰谷值的提取方法:提取每一个飞行任务剖面的每 一个飞行参数的峰谷值。
[0016] 优选的:所述步骤3中在雨流计数法统计之前，先要对飞行剖面峰谷值进行系列调 整，所述调整方法:先将步骤1中获得的每个飞行任务剖面的飞行参数所形成的波形进行封 闭处理，使其首尾数值相连。再将波形从步骤2获得的飞行剖面峰谷值中的最高峰值处分 开，使其首尾相连构成新的波形。
[0017] 优选的，所述步骤3中雨流计数法的统计方法：
[0018] 步骤31，每一个飞行任务剖面的每一个飞行参数的最高峰值或最低谷值为起点， 重新安排应力-时间历程。
[0019] 步骤32,起始于峰值的雨流将在下一个谷值处落下，或者起始于谷值的雨流将在 下一个峰值处落下，落下的雨流遇下述两种情况则停止：（a)起始于峰值的雨流遇到等于或 高于它的峰值便停止。（b)起始于谷值的雨流遇到等于或低于它的谷值便停止。
[0020] 步骤33,起始于峰值、谷值的雨流遇到自上而下的雨流便停止，此时根据该雨流的 起点和终点记取一个循环。
[0021] 步骤34,取出所有的全循环，并记录下各自的幅值和均值。
[0022] 优选的：所述步骤6中的计算载荷谱各参数之间的相关系数通过以下公式进行计 算：
[0025] 其中，rjk为矩阵X的第j列数据和第k列数据的相关系数，Xlj为各剖面各参数经过 雨流计数后所得的循环数，i = 1,2,...，n;j = 1,2,...，p，p表示载荷参数的个数，η表示剖 面个数，X表示Xlj所组成的矩阵。
[0026] 本发明的基于雨流计数循环的航空发动机载荷谱相关系数计算方法，相比现有技 术，具有以下有益效果：
[0027] 本发明提出的基于雨流计数循环的航空发动机载荷谱相关系数计算方法，是基于 与使用相关的载荷谱分析得到的，且去除了转速与过载值间的延迟现象，并由大量外场数 据进行验证，很大程度上代表机动循环的法向过载值与代表气动循环的转速的相关系数接 近于1，结果与实际更相符，因此弥补了现有的航空发动机载荷谱相关系数计算方法计算出 的机动载荷与气动载荷的相关系数偏小的问题。因此本发明可为航空发动机多参数载荷谱 的编制，进一步为航空发动机寿命预测与评估提供依据。
【附图说明】
[0028] 图1 一段典型的机动飞行载荷谱片断。
[0029]图2航空发动机转速与法向过载起始时间的相关关系。
[0030]图3航空发动机转速与法向过载峰值时间的相关关系。
[0031]图4航空发动机转速与法向过载之间的相关关系。
[0032]图5雨流计数法。
[0033]图6三峰谷值雨流计数的计算机程序框图。
[0034] 图7某型发动机某个飞行剖面的法向过载谱。
[0035] 图8基于雨流计数循环的航空发动机载荷谱相关系数计算方法流程图。
【具体实施方式】
[0036] 附图非限制性地公开了本发明一个优选实施例的结构示意图，以下将结合附图详 细地说明本发明的技术方案。
[0037] 实施例
[0038]航空发动机的机动飞行载荷谱是飞机在完成一次机动飞行任务过程中航空发动 机所受载荷随时间变化的曲线。一次典型的机动飞行任务通常由若干个机动飞行动作组 成，而每一个机动飞行动作必然会使航空发动机所受载荷产生一次或多次明显变化。因此， 航空发动机的机动飞行载荷谱就可以描述成由若干个幅值具有明显变化的载荷片断组成 的载荷一时间曲线。本文将这些载荷片断定义为机动飞行载荷任务段，简称任务段。附图1 是一段典型的机动飞行载荷谱片断。从图1中可以看出，该机动飞行载荷谱片断由四个不同 的任务段组成。第i个任务段从时刻!^开始，经过时间队后结束;第i+1个任务段从时刻T1+1开 始，经过时间D1+1后结束。图1中，Di是第i个任务段的持续时间，1\是第i个任务段的到达时 间，Μ是第i个任务段的峰值，Wi是第i个任务段与第i+Ι个任务段之间的等待时间，T 1+1-T^ 第i+1个任务段到达时间与第i个任务段到达时间之间的时间间隔。
[0039] 结合任务段分析发现，如附图4所示，法向过载系数任务段的出现与发动机转速变 化具有比较大的相关性。但是现有的机动载荷与气动载荷的相关性分析是将各载荷参数当 作一个完全随机系列来计算相关系数，计算出来的机动载荷与气动载荷的相关系数偏小。
[0040] 为了弥补现有的航空发动机载荷谱相关系数计算方法计算出的机动载荷与气动 载荷的相关系数偏小的问题，本发明提出了基于雨流计数循环的航空发动机载荷谱相关系 数计算方法，是基于与使用相关的载荷谱分析得到的，且去除了转速与过载值间的延迟现 象，并由大量外场数据进行验证，结果与实际情况更相符。由此进行机动载荷与气动载荷的 相关关系研究以及航空发动机多参数载荷谱编制，进一步为航空发动机寿命预测与评估提 供依据。
[0041] -种基于雨流计数循环的航空发动机载荷谱相关系数计算方法，如图8所示，对不 同飞行剖面中各参数的载荷谱(如转速谱、重心法向过载系数谱等)采用基于雨流计数法进 行统计。给出不同飞行剖面中各参数的幅值阈值，滤除幅值小于幅值阈值的载荷循环后，获 得各剖面各参数的循环数，同一飞行剖面不同的参数载荷谱得到不同的循环数，将不同参 数在各个剖面统计得到的循环数一一对应(如表3所示），每个载荷参数都可得到一列循环 数据，进而计算载荷谱各参数之间的相关系数。本发明所述的雨流法是一种常用的处理载 荷谱时间历程的双参数循环计数方法，通过统计载荷循环中的两个参量，反映载荷循环几 乎全部的信息。
[0042]具体包括以下步骤：
[0043] 步骤1，获取航空发动机的飞行任务剖面的个数以及每个飞行任务剖面的飞行参 数。
[0044] 设总共有η个飞行任务剖面，每个飞行任务剖面有m个飞行参数，例如航空发动机 法向过载Ny，侧向过载Nz，纵向过载Nx，发动机左发高压转子转速112_1^等。每个飞行剖面各 飞行参数有众多飞行数据。
[0045] 步骤2,飞行剖面峰谷值的提取:对步骤1获取的每个飞行任务剖面的飞行参数进 行飞行剖面峰谷值的提取。
[0046] 其中，飞行剖面峰谷值的提取方法:提取每一个飞行任务剖面的每一个飞行参数 的峰谷值。
[0047] 峰谷值提取的实质就是不计循环载荷的时间因素，将所处理参数的峰(谷)值，按 时间顺序从飞行剖面中提取出来，并形成按时间顺序排列的峰谷值序列，大大压缩了剖面 数据，便于进行参数循环计数。
[0048] 步骤3,参数循环计数:对步骤1中得到的每个飞行任务剖面的飞行参数根据步骤2 得到的飞行剖面峰谷值进行雨流计数法统计。
[0049] 对该剖面该参数数据进行雨流计数。在雨流计数之前，先要进行峰谷值系列调整， 以避免雨流后期出现发散-收敛型的波形，如附图5所示，其方法是先将步骤1中获得的每个 飞行任务剖面的飞行参数所形成的波形进行封闭处理，使其首尾数值相连。再将波形从步 骤2获得的飞行剖面峰谷值中的最高峰值处分开，使其首尾相连构成新的波形（图5所示）， 然后进行雨流计数。
[0050] 雨流计数法的统计方法：
[0051] 步骤31，每一个飞行任务剖面的每一个飞行参数的最高峰值或最低谷值为起点， 重新安排应力-时间历程。
[0052] 步骤32,起始于峰值的雨流将在下一个谷值处落下，或者起始于谷值的雨流将在 下一个峰值处落下，落下的雨流遇下述两种情况则停止：（a)起始于峰值的雨流遇到等于或 高于它的峰值便停止。（b)起始于谷值的雨流遇到等于或低于它的谷值便停止。
[0053]步骤33,起始于峰值、谷值的雨流遇到自上而下的雨流便停止，此时根据该雨流的 起点和终点记取一个循环。
[0054]步骤34,取出所有的全循环，并记录下各自的幅值和均值。
[0055] 应编制特定的程序由计算机自动完成雨流计数过程。对载荷循环的雨流计数有多 种编程方法，（附图6)为采用三峰谷值雨流计数的计算机程序框图。
[0056] 步骤4,无效循环幅值的去除：将各参数幅值小于[最小幅值+(最大幅值-最小幅 值)的10% ]剔除，再将幅值小于门槛值0 · 2的法向过载NY、纵向过载NX、侧向过载NZ剔除，获 得删除小幅值后的各参数幅值个数即载荷循环数，完成无效循环幅值的去除。无效幅值的 去除可以在雨流计数后进行，也可在雨流计数中进行。
[0057] 步骤5,循环的统计计数:根据步骤4去除的无效循环幅值的飞行参数获取各剖面 的各个参数经过雨流计数后的循环数，即幅值或均值的个数。
[0058] 设载荷谱总共有p个载荷参数，剖面数为η个，各剖面各参数经过雨流计数后所得 的循环数为xij(i = l，2,. . .，n; j = l，2,. . .，ρ)，则原始数据矩阵为：
[0060]步骤6,将步骤1获取的每个飞行任务剖面的飞行参数在步骤5获得的各剖面的各 个参数--对应(如表3所示），进而计算载荷谱各参数之间的相关系数。
[0061 ]计算得到各参数之间的相关系数矩阵R:
[0063]其中:rjk为矩阵X的第j列数据和第k列数据的相关系数，
[0066] 其中，rjk为矩阵X的第j列数据和第k列数据的相关系数，Xlj为各剖面各参数经过 雨流计数后所得的循环数，i = 1,2,...，n;j = 1,2,...，p，p表示载荷参数的个数，η表示剖 面个数，X表示Xlj所组成的矩阵。
[0067] ( - )某型发动机数十个飞行剖面的统计结果：
[0068] 共统计了 84个飞行剖面，每个飞行剖面统计了 10个参数，包括高压转子转速N2-L、 法向过载Ny、侧向过载Nz、纵向过载Nx、俯仰角速率FYJSL、横滚角速率HGJSL、航向角速率 HXJSL、X轴加速度Ax、Y轴加速度Ay、Z轴加速度Az;附图7为某个飞行剖面的法向过载Ny谱， 共6167个采样点，采样频率为1HZ。
[0069] (1)飞行剖面峰谷值的提取。提取某个剖面某个参数的峰谷值，峰谷值提取的实质 就是不计循环载荷的时间因素，将所处理参数的峰(谷)值，按时间顺序从飞行剖面中提取 出来，并形成按时间顺序排列的峰谷值序列，大大压缩了剖面数据，便于进行参数循环计 数。提取后各参数数据由原来的6167个降为了 2717个，大大压缩了剖面数据。
[0070] (2)参数循环计数。对该剖面该参数载荷谱进行雨流计数，在雨流计数之前，先要 进行峰谷值系列调整，以避免雨流后期出现发散-收敛型的波形。其方法是先将波形进行封 闭处理，使其首尾数值相连;再将波形从最高峰值处分开，使其首尾相连构成新的波形(如 附图5)，然后进行雨流计数。
[0071 ]雨流法的计数规则为：
[0072] a)重新安排应力-时间历程，以最高峰值或最低谷值为起点；
[0073] b)起始于峰值(谷值)的雨流将在下一个谷值(峰值)处落下，落下的雨流遇下述两 种情况则停止:（a)起始于峰值的雨流遇到等于或高于它的峰值便停止；
[0074] (b)起始于谷值的雨流遇到等于或低于它的谷值便停止；
[0075] c)起始于峰值、谷值的雨流遇到自上而下的雨流便停止，此时可根据该雨流的起 点和终点记取一个循环；
[0076] d)取出所有的全循环，并记录下各自的幅值和均值。
[0077]应编制特定的程序由计算机自动完成雨流计数过程。对载荷循环的雨流计数有多 种编程方法，（附图6)为采用三峰谷值雨流计数的计算机程序框图。对提取的峰谷值进行雨 流计数，最终得到的循环数为1358，最小幅值为0.005，最大幅值为1.3265。
[0078] (3)无效循环幅值的去除。构件的循环疲劳损伤主要取决于循环幅的大小。从峰谷 值序列中去掉一些对疲劳损伤影响不大的小量循环在工程上是非常必要的。将各参数幅值 小于[最小幅值+(最大幅值-最小幅值）的1 0 % ]剔除，再将幅值小于门槛值〇. 2的法向过载 NY、纵向过载NX、侧向过载NZ剔除，获得删除小幅值后的各参数循环数即幅值的个数。无效 幅值的去处可以在雨流计数后进行，也可在雨流计数中进行。
[0079] 表1某个飞行剖面法向过载值雨流计数后结果
[0082] (4)循环的统计计数。获得该剖面该参数载荷谱经过雨流计数后的循环数为37。
[0083] (5)获得该剖面其他参数基于雨流法计数后的循环数。参照上述方法，通过编程可 批量处理该剖面各参数数据，获得该剖面各参数经过雨流计数后的循环数。
[0084]表2某个飞行剖面各参数雨流计数后结果
[0086] (6)获得其他剖面各参数基于雨流法后的循环数。参照上述方法，通过编程可批量 处理其他剖面各参数数据，获得其他剖面各参数经过雨流计数后的循环数。
[0087] 表3多个飞行剖面各参数雨流计数后结果
[0092] (7)计算各参数之间的相关系数。根据上表中基于雨流统计的各剖面各参数的统 计结果，通过公式可以得到各参数之间的相关系数矩阵如下表。
[0094 ]表4基于雨流计数获得的载荷谱参数的相关系数
[0097] (8)比较传统方法计算出的相关系数与基于雨流法计算出的相关系数的区别，进 行分析。传统方法是将所有剖面拼接起来，统计各参数之间的相关系数。
[0098]表5基于传统的相关系数计算方法的载荷谱参数的相关系数
[0100]与传统的相关系数计算方法相比，基于雨流计数循环的相关系数计算方法结果显 示：
[0101 ] (1)有的参数间相关系数有所增大，有的系数间相关系数有所减少。例如：主要代 表机动载荷和气动载荷相关关系的法向过载Ny与转速N2L之间的相关系数由0.28156增大 为0.90617。呢1与啦的相关系数也分别由0.44401增大为0.87458<^与办之间的相关系数 由0.96594减为0.73446。
[0102] (2)代表气体力大小的发动机转速与机动载荷间相关系数都有明显的增大，尤其 是法向过载Ny与转速N2L之间的相关系数由0.28156增大为0.90617。剔除过载值与转速之 间相位差的因素后，得到的结果显然与实际情况更相符。
[0103] (3)部分参数间相关系数依然很小，基本接近0，这表明这些参数之间基本上可以 看作相互独立的。如侧向过载Nz与X轴加速度Αχ之间的相关系数只有-0.01175。
[0104]通过对比可见，本发明所述的基于雨流计数循环的航空发动机载荷谱相关系数计 算方法，是基于与使用相关的载荷谱分析得到的，且去除了转速与过载值间的延迟现象，并 由大量外场数据进行验证，很大程度上代表机动循环的法向过载值与代表气动循环的转速 的相关系数接近于1。结果与实际更相符。。
[0105]上面结合附图所描述的本发明优选具体实施例仅用于说明本发明的实施方式，而 不是作为对前述发明目的和所附权利要求内容和范围的限制，凡是依据本发明的技术实质 对以上实施例所做的任何简单修改、等同变化与修饰，均仍属本发明技术和权利保护范畴。
【主权项】
1. 一种基于雨流计数循环的航空发动机载荷谱相关系数计算方法，其特征在于:对不 同飞行剖面中各参数的载荷谱采用基于雨流计数法进行统计;给出不同飞行剖面中各参数 的幅值阔值，滤除幅值小于幅值阔值的载荷循环后，获得各剖面各参数的循环数，同一飞行 剖面不同的参数载荷谱得到不同的循环数，将不同参数在各个剖面统计得到的循环数一一 对应，每个载荷参数都可得到一列循环数据，进而计算载荷谱各参数之间的相关系数。2. 根据权利要求1所述的基于雨流计数循环的航空发动机载荷谱相关系数计算方法， 其特征在于，包括W下步骤： 步骤1，获取航空发动机的飞行任务剖面的个数W及每个飞行任务剖面的飞行参数； 步骤2,对步骤1获取的每个飞行任务剖面的飞行参数进行飞行剖面峰谷值的提取； 步骤3,对步骤1中得到的每个飞行任务剖面的飞行参数根据步骤2得到的飞行剖面峰 谷值进行雨流计数法统计； 步骤4,将各参数幅值小于[最小幅值+(最大幅值-最小幅值)的10%]剔除，再将幅值小 于口槛值0.2的法向过载NY、纵向过载NX、侧向过载NZ剔除，获得删除小幅值后的各参数幅 值个数即载荷循环数，完成无效循环幅值的去除;无效幅值的去除可W在雨流计数后进行， 也可在雨流计数中进行； 步骤5,根据步骤4去除的无效循环幅值的飞行参数获取各剖面的各个参数经过雨流计 数后的循环数，即幅值或均值的个数； 步骤6,将步骤1获取的每个飞行任务剖面的飞行参数在步骤5获得的各剖面的各个参 数一一对应，进而计算载荷谱各参数之间的相关系数。3. 根据权利要求2所述的基于雨流计数循环的航空发动机载荷谱相关系数计算方法， 其特征在于:所述步骤2中飞行剖面峰谷值的提取方法:提取每一个飞行任务剖面的每一个 飞行参数的峰谷值。4. 根据权利要求2所述的基于雨流计数循环的航空发动机载荷谱相关系数计算方法， 其特征在于:所述步骤3中在雨流计数法统计之前，先要对飞行剖面峰谷值进行系列调整， 所述调整方法:先将步骤1中获得的每个飞行任务剖面的飞行参数所形成的波形进行封闭 处理，使其首尾数值相连;再将波形从步骤2获得的飞行剖面峰谷值中的最高峰值处分开， 使其首尾相连构成新的波形。5. 根据权利要求2所述的基于雨流计数循环的航空发动机载荷谱相关系数计算方法， 其特征在于，所述步骤3中雨流计数法的统计方法： 步骤31，每一个飞行任务剖面的每一个飞行参数的最高峰值或最低谷值为起点，重新 安排应力-时间历程； 步骤32,起始于峰值的雨流将在下一个谷值处落下，或者起始于谷值的雨流将在下一 个峰值处落下，落下的雨流遇下述两种情况则停止：（a)起始于峰值的雨流遇到等于或高于 它的峰值便停止；（b)起始于谷值的雨流遇到等于或低于它的谷值便停止； 步骤33,起始于峰值、谷值的雨流遇到自上而下的雨流便停止，此时根据该雨流的起点 和终点记取一个循环； 步骤34，取出所有的全循环，并记录下各自的幅值和均值。6. 根据权利要求2所述的基于雨流计数循环的航空发动机载荷谱相关系数计算方法， 其特征在于:所述步骤6中的计算载荷谱各参数之间的相关系数通过W下公式进行计算：其中，nk为矩阵X的第j列数据和第k列数据的相关系数，XU为各剖面各参数经过雨流计 数后所得的循环数，i = l，2,. . .，n; j = l，2,. . .，p，p表示载荷参数的个数，η表示剖面个数， X表示义^所组成的矩阵。
【文档编号】G06F19/00GK106096262SQ201610402859
【公开日】2016年11月9日
【申请日】2016年6月8日
【发明人】宋迎东, 孙志刚, 邢广鹏
【申请人】南京航空航天大学