外文翻譯--對薄的橫向各同性彈性層軸對稱接觸的漸近解

附錄中文譯文 對 薄 的 橫向各同性彈性層軸對稱接觸 的漸近解 摘要 對剛性球與一個(gè)薄的，橫向各同性彈性層正常接觸的問題進(jìn)行研究，薄的彈性層建立在剛性基礎(chǔ)上。無論是在極限情況下，其中的接口可以是理想的粘結(jié)或無粘結(jié)摩擦，并被認(rèn)為是用于球形壓頭下的接觸壓力和接觸補(bǔ)片的半徑近似解。通過用有限元法預(yù)測和比較這些近似解，來確定解決方案的準(zhǔn)確性。 關(guān)鍵詞 ：固體潤滑劑 ；各向異性涂層 ；有限元法 ；接觸壓力 1 引言 表面涂層已成為廣泛使用在工程領(lǐng)域，由于它們提供了許多潛在的好處。這些涂層可以是各向同性或各向異性的性質(zhì)，通常被用于提高 接觸特性及改進(jìn)的配合部件的摩擦性能。表面涂層最普遍的應(yīng)用是在運(yùn)行于真空或低壓環(huán)境中的一些滾動(dòng)元件或滑塊軸承 1。涂層通常是非常薄層軟金屬（如金，銀， 鋅及鉛 ）層固體（如二硫化鉬和硼氮化物）或聚合物（例如，聚四氟乙烯，聚酰胺復(fù)合材料和酚醛樹脂和環(huán)氧樹脂），其被沉積到襯底上。當(dāng)固體潤滑劑的涂層被施加到基材上時(shí)，接觸應(yīng)力和位移場顯著從赫茲接觸輪廓的半空間內(nèi)偏離。這是設(shè)計(jì)時(shí)的顯著問題，因?yàn)闆]有可用于分析三維涂層表面接觸方法，所以即使簡單的幾何形狀封閉形式也沒很好的解決方案。 回顧文獻(xiàn)，許多研究人員調(diào)查解析與各向 同性彈性層 2-12 軸對稱接觸得問題。各向異性涂料展示更復(fù)雜的接觸特性，不是因?yàn)樗麄兌ㄏ虻囊蕾嚫飨蛲詮椥詫印?Ovaert 13 分析了橫觀各向同性固體潤滑膜的剛性橢球壓頭的打孔問題。 Kuo 和科爾 14 通過漢克爾方式取得橫觀各向同性多層介質(zhì)數(shù)值解變換。 Lovell和洛弗爾以及 Khonsari 16采用有限元方法來研究具有彈性基板上的橫向各向同性彈性層接觸的彈性球法向應(yīng)力和切向的摩擦特性。為了更好地了解層狀介質(zhì)的接觸，本文的重點(diǎn)是剛性的，無摩擦的橫觀各向同性層上剛性基板的軸對稱 縮進(jìn)。彈性層和剛性基板之間的界面上被假定為理想的粘結(jié)或無粘結(jié)摩擦。所采取的方法是延續(xù)約翰遜橫觀各向同性彈性涂料的“平面截面后仍保持平面壓縮”的假設(shè)2。 2 問題描述 在考慮剛性的基礎(chǔ)上，其中旋轉(zhuǎn)對稱層的軸是垂直于所述層和所述襯底之間界面的無限橫觀各向同性彈性層。在所述層和襯底之間的界面將被假定為理想的粘結(jié)或無粘結(jié)摩擦。一個(gè)剛性的球形壓頭穿透層法線方向的自由表面，在表面壓頭與層之間的界面被視為無摩擦。近似解將取決于該層的位移場和應(yīng)力場。這個(gè)問題是軸對稱球形壓頭大約含有多少垂直于橫觀各向同性層和基礎(chǔ)之間界面 的質(zhì)心線。 2.1 控制方程 圓柱坐標(biāo) 系（ r ， ， z）與相關(guān)聯(lián)的正交向量基 zr eee ， ，它定義在該層自由表面的平面，并與對稱軸以及 z 軸一致。如圖所示 1，在其引用的位置，也就是說，之前的壓痕彈性層占據(jù) 0 Z z = 0 處，并定義其自由表面區(qū)域之間的界面層與剛性基層在 z =t 萬噸。該問題的對稱性決定了位移場 u，且具有代表性的是zzrr eueuu ，同時(shí)標(biāo)量分量ru和zu與是相互獨(dú)立的： 假設(shè)變形小， | | u | | 1，并且它遵循從方程（ 1 ），該位移具有代表性）（ rzzrrzzzzzrrrr eeeeeeeeee ，標(biāo)量分量 rr， zz 和 rz，這些是在給定位移場條件的標(biāo)量分量。 因此，應(yīng)力場具有代表性的 ）（rzzrrzzzzzrrrr eeeeeeeeee ，與標(biāo)量分量 rr， zz 的指向和 rz 是由給定應(yīng)變場的標(biāo)量分量的計(jì)算 其中 C11， C12， C13， C33 和 C44 是為橫觀各向同性彈性層的五個(gè)獨(dú)立的彈性常數(shù)。最后，假設(shè)自身重量和物質(zhì)加速度場可以忽略不計(jì)，所以平衡方程 0 只有兩個(gè)非平衡的標(biāo)量分量。 第三個(gè)標(biāo)量平衡方程是平衡滿足時(shí)。 2.2 邊界條件 在變形的結(jié)構(gòu)中，彈性層和剛性球形壓頭之間的接口是一個(gè)球形帽，具有一個(gè)圓形的正交投影，或接觸印痕，（參見圖 1）這個(gè)接口被假定為無摩擦，使邊界條件： 其中 w（ r）是球形帽的公式，這樣，如果 R 是球形壓頭和 w* W（ 0）時(shí)壓痕的深度半徑，則 如果該彈性層是相對于壓頭很薄， 且半徑 R t，瓦特 tw* ,則它遵循 a R 并且 222 2r-1r-1 RR ，當(dāng)所有 r a 時(shí)。因此，對于薄彈性層的壓痕更新方程（ 7）由下式給出 R2r-r 2* ； （ 7） 注意，接觸印痕半徑是關(guān)系到球面壓頭半徑 R 和壓痕深度 W *。要求當(dāng) W（ A） = 0 時(shí)，它對于一個(gè)薄的彈性層的公式如（ 8）。 ,2 r-ar 22 R R2a2* . (9) 該壓痕深度寬 W*比彈性層厚度 t 少 ,可以由不等式條件表示 : t2ta R （ 10） 該彈性層將被假定為從現(xiàn)在開始要薄，并且近似表達(dá)式方程為（ 9），壓痕輪廓和深度將用于接觸補(bǔ)片的外面，在彈性層的自由表面牽引是無邊界。 00zrz0zzz ，ar （ 11） 邊界條件方程（ 6），（ 9）和（ 11）在自由表面處 z = 0 處，都是相似的考慮以所有情況。用于彈性層和剛性基礎(chǔ)之間的界面為 z= t 時(shí)，將要考慮的是該接口在理想化的情況下是理想粘結(jié)，其相應(yīng)的邊界條件是 0uu tzztzr ， ，0r （ 12） 這將被稱為鍵合的情況。其他理想化將要考慮的是這種情況，其中界面彈性層和剛性基礎(chǔ)之間無粘結(jié)和摩 擦，從而使相應(yīng)的邊界條件是 ,0tzrz ， 0u tzz ， ，0r （ 13） 這種將被稱為無粘結(jié)情況。 3，分析解決方案 約翰遜假設(shè)軸對稱變化 2說：“飛機(jī)的部分壓縮后仍然面壓縮”意味著這是 ru 與 z 無關(guān)： ）（ ruu rr ； （ 14） 首 先，注意，該假設(shè)意味著rzrz uc44 ； （ 15） 這意味著壓頭之間的界面在彈性層中，零摩擦邊界條件式（ 6） 1 和正常位移邊界條件方程（ 6） 2 不能同時(shí)滿足。優(yōu)先此處是考慮到正常的位移邊界條件。在另一方面，零摩擦邊界條件式（ 13） 1 是隱含由零正常位移邊界條件式（ 13） 2。 給定了位移場的假定形式后，它不可能找到平衡方程的解，方程（ 4）和（ 5）滿足所有的邊界和內(nèi)部連續(xù)條件。因此，該解決方案在下面衍生，它們近似值可能非平衡或一些邊界以及內(nèi)部連續(xù) 條件可能不會(huì)滿足。不滿足條件的一些誤差程度明確識(shí)別并將指示給予有限元分析比較。 3.1 理想情況下綁定的接口 用于粘合的情況下，徑向位移邊界條件式（ 12） 1 和假設(shè)方程（ 14）意味著徑向位移消失無處不在的彈性層： 0ur ， 0rr ， ，0r ； t0z ， （ 16） 注意的是軸向位移的邊界條件方程（ 6） 2 和（ 12） 2 及平衡方程方程（ 4）和 （ 5）與 此結(jié)果不兼容，也就是說，它們不能同時(shí)滿足。因此，判定為滿足軸向位移的邊界條件為平衡方程。 假定應(yīng)變分量zz是通過均勻彈性層的厚度 . )(zz rzz ； (17) 它遵循由應(yīng)變 - 位移關(guān)系式（ 2） 3， 軸向位移的邊界條件方程（ 6） 2 和 （ 12） 2 和牽引無邊界條件式（ 11），其 因此，該球形壓頭的下方，該組件壓力由下式給出； 然而，牽引 - 自由表面之 下，應(yīng)力場消失： 0 ， ，ar （ 20） 這種近似解滿足所有邊界除摩擦壓條件（ 6） 1 ，該解決方案對應(yīng)的不是一個(gè)理想的結(jié)合壓頭和彈性層之間的界面。另外非平衡的應(yīng)力場和剪切牽引 rz 是在不連續(xù)的整個(gè)表面R = A 上。 壓頭之間的接觸壓力 P（ r）的彈性層為 t2 r-a-r22330zzz RCP （ 21） 并與由下式給出的總壓頭載荷 P； drrr2 a0 PP （ 22） 如下的接觸補(bǔ)片的半徑是與壓頭通過加載； 4133t4a CPR (23) 有趣的是注意到，接觸壓力的分布與接觸面半徑的這個(gè)近似解為獨(dú)立四五個(gè)彈性常數(shù)； 3.2 無粘結(jié)，無摩擦界面 在制定的近似解解析時(shí)，應(yīng)選擇其中的邊界條件，以滿足在某些特別的區(qū)域所得到的值。當(dāng)然，最終的目的，是獲得一個(gè)模型，它提供的預(yù)測，是最接近系統(tǒng)的響應(yīng)。鑒于牽引邊界在方程中給出的條件（ 6） 1 及（ 13） 1，我們將假定為一個(gè)無粘結(jié)的材料，根據(jù)剪應(yīng) 力分量 rz 壓頭是可以忽略不計(jì)的。 0rz a0r ， （ 25） 它遵循由式（ 2） 3 和軸向位移邊界條件方程（ 6） 2，（ 13） 2 是； ）（ z-tt2 r-au 22z R a0r ， （ 26） 結(jié)合方程（ 24）和（ 26），（ 2）和（ 3），把以上結(jié)果代入式（ 4）給出的平衡方程為tcrc-u-drdurdrudr11313rr2 r22R （ 27） 一般的解決這個(gè)歐拉 - 柯西方程為 tc8c-rru 111321r RBB a0r ， （ 28） 那 ru 等于零在 r=0 和 B2=0 的情況下。將方程（ 26）和（ 28）代入（ 2）和（ 3），并要求這里的 zz =0 時(shí)，當(dāng) r a 然后給出： 這種近似解滿足所有邊界條件為 0 r a 時(shí)。然而平衡條件式（ 5）不成立，是因?yàn)槲灰茍龇匠蹋?26）和 （ 28）是在假定方程（ 24）不一致 ，此外， 0 r a 為獨(dú)立方程的一個(gè)近似解，此領(lǐng)域之外的彈性層響應(yīng)的接觸壓力 rP由下式給出； 221121333110zzzr r-atc2c-cc-RP （ 30） 總壓頭的載荷 P 由公式（ 22）給出，相關(guān)的壓頭負(fù)載由該接觸結(jié)合半徑得 41132331111c-ccRtpc4a （ 31） 需要注意的是，在這種情況下，接觸壓力分布和接觸面半徑是獨(dú)立的彈性常數(shù) 44c 。 4 結(jié)果與 討論 對于提出粘結(jié)和無粘結(jié)分析的解決方案在第 3.1 和 3.2，但更重要注意的是有適用性的范圍是有限。這兩種近似解決方案只適用于可壓縮，橫向各向同性材料。為了滿足這些條件，材料的剛度矩陣（見式（ 3）必須是穩(wěn)定和正定： 當(dāng)上述不等式的左邊側(cè)接近為零時(shí)，材料變得不可壓縮，同時(shí)分析解決方案將變得非常不準(zhǔn)確。 4.1 對于降低各向同性材料的解決方案 上面導(dǎo)出的結(jié)果可以通過以下方式獲得 初始 近似 的準(zhǔn)確性，檢查對于各向同性的情況下，降低溶液的涂料。對于各向同性材料，橫觀各向同性中出現(xiàn)的本方程彈性常數(shù) 式（ 3）減少到 其中， E 是楊氏模量，是泊松比。它可以認(rèn)為，對于各向同性的彈性層近似壓力分布公式（ 21）用于粘合的情況下，與壓力分布公式（ 30）用于在粘合的情況下。 這些分別都是由賈法爾 11派生的近似結(jié)果。近似解方程（ 21）和（ 23）用于結(jié)合情況，而方程（ 30）和（ 31）用于在粘合的情況下，兩者可以看作是賈法爾的解決方案擴(kuò)展到橫觀各向同性彈性層。如將要討論下面的部分，方程（ 21）和（ 30）也將與有限元分析預(yù)測結(jié)果。 4.2 有限元分析比較 商用有限元分析軟件 ABAQUS6.3.1 用于預(yù)測的應(yīng)力和位移，用于解決橫 觀各向同性彈性層中的以上問題。二維軸對稱元件被用于彈性層，用來解析球形壓頭剛性表面。該無限彈性層被建模為一個(gè)圓柱體，其半徑為0r= 1 毫米，厚度 t=5 微米。而壓頭的半徑 R=8 毫米。從而 200tr0 ； tR1600 （ 38） 并且從方程（ 10）得出該接觸補(bǔ)丁半徑被限制為 6.56ta （ 39） 盡管人們不會(huì)準(zhǔn)確的 認(rèn)為 預(yù)測前模型的印跡半徑值接近這一限制 值 。網(wǎng)格 5 uM0.5uM中 大約有 2000 個(gè)元素內(nèi)的彈性層。 對于 NbSe2 橫觀各向同性彈性常數(shù)有 C11=106 GPa，C33=54 GPa， C12=14 GPa， C13=31 GPa 和 C44=19.5 GPa。 最大接觸壓力會(huì)發(fā)生在中心 r = 0 的結(jié)合處，在如圖 2 所示，預(yù)測 P（ 0）為 這圖（ 2）模型預(yù)測超過了最大接觸壓力0P,有限元預(yù)測作為無量綱的函數(shù)接觸面半徑為 a/噸。 圖 3，接觸壓力分布在粘結(jié)和無粘結(jié)的情況下，如圖通過分析模型，并在 a/t=10 進(jìn)行有限元分析預(yù)測。 把近似解的最大接觸壓力方程（ 21）和（ 30）進(jìn)行比較，如果有限元分析接觸印痕半徑的值不同。則結(jié)果0P的無量綱接觸面半徑比 a/t 也不同。其中對于任何有限元分析的結(jié)果是粘結(jié)或無粘結(jié)的情況，適當(dāng)?shù)膶τ谶@兩種情況在結(jié)合和未結(jié)合的情況下，把那些好的方案模型進(jìn)行有限元分析時(shí)預(yù)測得到的 a/t 值比 5 更大，從接觸壓力分布的近似解方程（ 21）和（ 30）中用有限元素分析示圖 3 中得出 一個(gè)固定的接觸印痕半徑，其比值 a/ t 為 10。檢查該圖中可以得出的結(jié)論是近似 解解析與對應(yīng)的有限元分析有較好的一致性。 5 總結(jié)與結(jié)論 把剛性球放到一個(gè)薄的，橫向的剛性基片上時(shí)對各向同性彈性層進(jìn)行問題分析。該彈性層和剛性基礎(chǔ)之間的接口被假定為理想的粘結(jié)或無粘結(jié)的摩擦。約翰遜的假設(shè)的推廣使用得出近似解，有些用有限元分析來確定它們的相對精度，并進(jìn)行結(jié)果比較。從分析結(jié)果以及它們相對于有限元分析的結(jié)果來看，得到以下結(jié)論。 粘結(jié)和無粘結(jié)的近似解用于減少到賈法爾的解決方案極限的各向同性彈性層。 在結(jié)合情況下的近似 解逼近有限元分析時(shí)，其預(yù)測作的接觸半徑將增大。 對于無粘合的情況下，基于微不足道的剪應(yīng)力上假設(shè)的近似解與有限元的計(jì)算結(jié)果吻合良好。 參考文獻(xiàn) 1 M.R.洛弗爾， M.M.Khonsari， R.D.馬蘭戈尼 .為二硫化鉬涂層球軸承的摩擦分析：三維有限元分析 .ASME .賽報(bào) ,1997,11(9):754-763. 2 K.L.約翰遜 .接觸力學(xué) .英國劍橋：劍橋大學(xué)出版社， 1985. 3 I.N. Sneddon.傅立葉變換 .紐約：麥格勞希爾， 1964. 4陳 W.T.應(yīng)力和位移的層狀彈 性介質(zhì)詮釋的計(jì)算 J.工程科學(xué)， 1972， 9 ： 775-800. 5 W.T.有機(jī)化學(xué)， P.A.恩格爾 .沖擊和接觸應(yīng)力分析多層介質(zhì) .詮釋固體力學(xué)，1972,8:1257-1281. 6 MJ Matthewson.薄涂層符合軸對稱接觸 J.固體報(bào)， 1981,2（ 9）： 89-113. 7 V.M.亞歷山德羅夫 .關(guān)于彈性和數(shù)學(xué)物理理論的一些積分方程的近似解 .數(shù)學(xué)機(jī)甲 ,1967,