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