對接觸線上每個單元的應力時程分別進行雨流計數，將各單元的應力時程儲存于一個二維數組中，二維數組的每行表示同一時刻接觸線各單元的應力，每列表示接觸線某一單元的應力時程，具體實現過程如下:
(1)相鄰等值數壓縮。分別對二維數組的每列判斷應力時程前后兩個元素是否相等，若相等則只保留一個。
(2)提取谷峰值。判斷經壓縮后的二維數組每列相鄰數值的差的乘積是否大于零，若大于零則提取峰(或谷)值。
(3)谷峰值調整。對二維數組的每列判斷峰谷值點總數是奇數還是偶數，若是奇數，則對經以上兩步處理過的應力時程不做變化；若是偶數，則去掉最后一個點。判斷首尾點是峰值還是谷值，若是峰值，則使首尾兩點的值都取二者中的較大者；若是谷值，則使首尾兩點的值都取二者中的較小者。
(4)對接。尋找調整過的應力時程二維數組的每列最高波峰(低波谷)點，將應力時程的每列從該點處截斷，將左段的起點與右段的末點對接，使新的應力時程首尾皆為最高波峰(低波谷)點。
(5 )用四峰谷值法進行一次雨流計數，4個點的選取遵循單向選取、循環進行的原則提取循環應力幅值和應力均值，并記錄相應的應力循環數。
6)雨流計數循環。對一次雨流計數后剩下的點編制循環命令流，繼續提取應力循環直至剩下3個點(即是數組中最值構成的整循環)為止。對提取的應力循環、每次雨流計數剩下的點、用于刪除計數點的數組 等分別存放于不同的二維數組；用一維數組存放雨流計數剩下的點和應力循環每列的個數，并在每次雨流計數后更新；用于刪除計數點的數組和雨流計數剩下的點的數組初始值為應力時程經前四步處理后的數組，在每次雨流計數后更新。
(7 )無效幅值去除。用變程閾值公式去除幅值很小的可以忽略不計的應力循環。
8)按1MPa的應力幅分級范圍統計不同應力幅下的應力循環數目。 經過以上雨流計數過程，得到接觸線單元的應力幅值集合、應力均值集合、相應的應力循環數集合。