牛牛天天人人综合影院 每列表示接触线某一单元的应力时程接触线的实现过程

       对接触线上每个单元的应力时程分别进行雨流计数，将各单元的应力时程储存于一个二维数组中，二维数组的每行表示同一时刻接触线各单元的应力，每列表示接触线某一单元的应力时程，具体实现过程如下:
(1)相邻等值数压缩。分别对二维数组的每列判断应力时程前后两个元素是否相等，若相等则只保留一个。
(2)提取谷峰值。判断经压缩后的二维数组每列相邻数值的差的乘积是否大于零，若大于零则提取峰(或谷)值。
(3)谷峰值调整。对二维数组的每列判断峰谷值点总数是奇数还是偶数，若是奇数，则对经以上两步处理过的应力时程不做变化；若是偶数，则去掉最后一个点。判断首尾点是峰值还是谷值，若是峰值，则使首尾两点的值都取二者中的较大者；若是谷值，则使首尾两点的值都取二者中的较小者。
(4)对接。寻找调整过的应力时程二维数组的每列最高波峰(低波谷)点，将应力时程的每列从该点处截断，将左段的起点与右段的末点对接，使新的应力时程首尾皆为最高波峰(低波谷)点。
(5 )用四峰谷值法进行一次雨流计数，4个点的选取遵循单向选取、循环进行的原则提取循环应力幅值和应力均值，并记录相应的应力循环数。
6)雨流计数循环。对一次雨流计数后剩下的点编制循环命令流，继续提取应力循环直至剩下3个点(即是数组中最值构成的整循环)为止。对提取的应力循环、每次雨流计数剩下的点、用于删除计数点的数组 等分别存放于不同的二维数组；用一维数组存放雨流计数剩下的点和应力循环每列的个数，并在每次雨流计数后更新；用于删除计数点的数组和雨流计数剩下的点的数组初始值为应力时程经前四步处理后的数组，在每次雨流计数后更新。
(7 )无效幅值去除。用变程阈值公式去除幅值很小的可以忽略不计的应力循环。
8)按1MPa的应力幅分级范围统计不同应力幅下的应力循环数目。 经过以上雨流计数过程，得到接触线单元的应力幅值集合、应力均值集合、相应的应力循环数集合。