本發(fā)明公開了一種基于ABAQUS二次開發(fā)的花鍵磨損與疲勞壽命仿真方法，先建立內(nèi)外花鍵模型，在ABAQUS中完成前處理設(shè)置，求解當(dāng)前循環(huán)步得到花鍵旋轉(zhuǎn)一周各接觸節(jié)點(diǎn)的接觸壓力、應(yīng)力、應(yīng)變和相對(duì)滑移距離；采用Fortran編寫UMESHMOTION子程序，基于Archard方程計(jì)算各個(gè)節(jié)點(diǎn)的磨損深度，通過(guò)自適應(yīng)網(wǎng)格技術(shù)更新網(wǎng)格；利用宏命令提取各節(jié)點(diǎn)應(yīng)力和應(yīng)變，利用三種方法計(jì)算當(dāng)前循環(huán)步各節(jié)點(diǎn)疲勞壽命，計(jì)算出各節(jié)點(diǎn)在所有循環(huán)步的累積損傷值D并繪制磨損分析圖像。本發(fā)明利用ABAQUS建立花鍵磨損和疲勞壽命分析模型，可準(zhǔn)確判斷花鍵失效形式，為花鍵壽命提升提供指導(dǎo)，多種臨界平面法可同時(shí)考慮鍵齒表面的壓縮應(yīng)力狀態(tài)和齒根區(qū)域的拉伸應(yīng)力狀態(tài)，使分析結(jié)果更準(zhǔn)確。