摘要：在扭動(dòng)微動(dòng)的應(yīng)力分布解析計(jì)算中,通過接觸表面切應(yīng)力可以推導(dǎo)出接觸體內(nèi)部的應(yīng)力。而轉(zhuǎn)矩卸載過程中,由于接觸表面切向力表達(dá)式復(fù)雜,難以推導(dǎo)出接觸體內(nèi)部的應(yīng)力。為此采用快速傅里葉變換的方法對(duì)球-球接觸模型進(jìn)行數(shù)值模擬,以此獲得卸載時(shí)接觸體內(nèi)部的應(yīng)力分布。首先,計(jì)算出接觸表面加載時(shí)的應(yīng)力分布,并得到轉(zhuǎn)矩與扭轉(zhuǎn)角之間的遲滯曲線。然后,分析法向力和加卸載轉(zhuǎn)矩共同作用下接觸體內(nèi)部的應(yīng)力分布,以及最大應(yīng)力的大小與位置。最后,在遲滯曲線上標(biāo)明最大應(yīng)力點(diǎn)位置,并分析黏著區(qū)半徑c和摩擦因數(shù)f對(duì)最大應(yīng)力位置的影響規(guī)律。研究表明:卸載時(shí)表面切向力qr沿半徑方向,先增大再減小最后持續(xù)增大;法向力在較大范圍內(nèi)影響接觸體應(yīng)力分布,而轉(zhuǎn)矩主要在深度0~0.5a范圍內(nèi)起作用;黏著區(qū)半徑c越大,摩擦因數(shù)f越小,最大應(yīng)力點(diǎn)更可能位于接觸體內(nèi)部。