含界面裂紋圓形夾雜的十二次對稱二維準晶熱應力分析

翟 婷, 趙雪芬, 丁生虎

(1. 寧夏大學數學統計學院，銀川 750021; 2. 寧夏大學新華學院，銀川 750021)

0 引言

由于熱應力問題在新材料，核工程以及航空等領域的重要性，含夾雜、孔洞或裂紋的復合材料熱應力問題受到了廣泛的關注。Chao 和Shen[1–2]借助Hilbert 公式以及復應力函數法研究了粘結異種材料曲線裂紋的熱彈性問題并且給出了圓形夾雜熱彈性問題的一般解答。Kattis 和Meguid[3]采用二相復勢的方法研究了曲線邊界夾雜的平面熱彈性問題，獲得了橢圓形夾雜在均勻溫變下的熱彈性場。Ding 等[4–5]根據奇異積分方程方法計算了混合邊值問題并獲得了在均勻熱流作用下兩種粘結正交各向異性介質中界面裂紋問題的平面彈性解；運用疊加原理和傅里葉積分變換，對在熱載荷作用下，正交異性涂層-基底結構內埋裂紋的熱彈性問題進行了分析。

準晶作為一種新型智能材料，與金屬合金相比，它形成了一系列具有奇特物理和機械性能的特殊結構。由于準晶體脆性高，其工作環境多為高壓高溫等更苛刻的條件，容易產生夾雜、裂紋或孔洞等缺陷，這會嚴重降低器件的性能，故對含缺陷準晶材料的研究是十分有必要的。Hu 等[6]、Fan 和Mai[7]對準晶體的物理性質及其對材料科學的意義進行了全面的概述。近年來，含缺陷準晶的力學分析得到了充分的研究[8–14]。劉媛等[15]采用Muskhelishvil 復變函數法研究了十次對稱二維準晶材料中橢圓孔邊部分裂紋面受到均布應力作用問題。Dang 等[16]運用位移不連續法研究了一維六方準晶涂層在反平面載荷作用下界面裂紋的斷裂行為，用傅里葉變換方法導出了反平面位移不連續的基本解，分析了裂紋前緣應力的奇異性。

在考慮熱效應時，準晶材料含缺陷問題的研究相對較少。Yang 等[17]借助非局部理論，得到了以厚度方向為準周期方向的二維均勻十面準晶納米板的精確解，給出了二維準晶簡支納米板在上表面溫度變化下的三維熱彈性解析解。王旭等[18–19]推導了考慮熱效應時十次對稱二維準晶體平面應變問題的通解表示，將穩恒線熱源作用下變形問題劃為平面應變問題進行求解，給出了均勻熱流作用下十次對稱二維準晶體中一絕緣橢圓孔洞引起的熱應力問題的解答。Wang 和Schiavone[20]討論了含橢圓夾雜的二維十次對稱準晶在熱力耦合作用下的平面熱彈性問題。Wang 和Schiavone[21]結合復變函數和傳遞矩陣方法研究了含橢圓夾雜的N層十次準晶復合材料在三種不同熱載荷作用下的熱力學響應解。Li 和Guo[22]通過復變函數方法研究了八面準晶含橢圓孔的熱應力分析。Fan 等[23]基于傅里葉變換，運用擴展位移不連續法對一維六方準晶周期裂紋的平面熱效應進行了分析。Li 等[24]給出了均布熱流作用下一維六方準晶含幣形裂紋問題的解析解。文獻[25]借助廣義勢理論法研究了二維六方準晶體無限空間中含半無限平面裂紋的三維基本熱彈性問題。本文利用復變函數方法研究了十二次對稱二維準晶在點熱源作用下的共弧圓形夾雜界面裂紋的平面問題，本問題的研究未見相關報道。

1 問題描述與基本公式

如圖1 所示，無限大十二次對稱二維準晶基體S?中包含一個半徑為R的圓形彈性夾雜S+。在無限大準晶基體S?內任意一點z0作用著一點熱源，其強度為q0。同時，在準晶基體與彈性夾雜界面上有n條互不相交的裂紋L′1,L′2,··· ,L′n，裂紋的尖端分別用aj和bj表示，記L′=L′1+L′2+···+L′n。除裂紋L′外的界面剩余部分記為L。不失一般性，假設在L′上無面力作用。當準晶基體與彈性夾雜粘結完好時，位移和應力在L上連續。由于十二次對稱二維準晶體的聲子場和相位子場不耦合，為了界面L上滿足位移和應力連續則相位子場上的位移和應力全為0。

圖1 點熱源作用在圓形夾雜界面裂紋附近

設圓形夾雜的中心在復平面z=x1+ix2的坐標原點，則夾雜與準晶基體界面上的點可表示為t=Reiθ。此時，界面上位移和應力的邊界條件可表示為

其中σr和σrθ表示極坐標中彈性夾雜和準晶基體聲子場的應力分量。Hr和Hrθ表示極坐標中準晶相位子場的應力分量。ur和urθ是極坐標系中彈性夾雜和準晶聲子場的位移分量。最后一個下標1 和2 分別表示區域S+和S?，上標“+”和“?”表示當z分別從S+和S?趨向于界面時函數所取的邊界值。

由二維彈性復勢理論及熱彈性問題的基本公式[26]，圓形夾雜內的彈性場可用兩個復勢應力函數Φ(z)和Ψ(z)再加一個溫度復函數g′1(z)來確定，即

由文獻[27]知，十二次對稱二維準晶基體的彈性場可用兩個復勢函數?(z)、Θ(z)和一個溫度函數g′2(z)來確定，即

由二維熱傳導理論知，熱流Q，熱流強度q和溫度T可以由單個復溫度函數表示為

其中κt表示熱傳導系數，q1、q2和qρ、qφ分別表示熱流強度的直角坐標分量和極坐標方向的分量，Re 和Im 分別表示對復數函數取實部和虛部。

2 溫度場

準晶基體中的復溫度函數可以表示為

考慮到界面裂紋是非滲透的，邊界上熱流強度和溫度連續邊界條件可表示為

其中qρ表示熱流強度徑向極坐標分量。

設Fi(z) =g′i(z)，利用邊界條件(12)和(13)，由解析延拓技術和奇性主部分析方法可求得兩個區域復勢函數的一般解答為

其中G0(z)、G∞(z)、Gz0(z)和Gz?(z)表示被積函數N1G(z)/X0(z)在z=0, ∞, z0, z?點處的奇性主部，并且

對式(15)，還需要確定n個常數c1,c2,··· ,cn。把式(15)在無窮遠域展開后與F′2(z)比較1/z前的系數，求得1/z前的系數。為了確保L上溫度的連續性，補充n個端點溫度相等條件。先假定端點a1處溫度確定，則可得到n ?1 個方程

式(16)是n ?1 個方程，可確定n ?1 個未知數，至此F′1(z)和F′2(z)可完全確定，這樣表示整個溫度場也完全確定。

3 聲子場的熱應力分析

如問題所述，當準晶基體S?內任意一點z0作用一強度為q0的點熱源時，可設夾雜S+內和準晶基體聲子場的復勢函數為

同理，在S?內引入函數

這樣，Φ(z)可經過圓弧L′解析到S?內。

考慮邊界條件(2)～(4)，把式(20)和式(23)代入，由廣義Liouville 定理[26]，可得

其中?1(z)=?4Mi(L+M)?(z)且D0為常數。

對式(7)關于θ求導，由邊界條件式(1)及式(10)、式(24)得

把式(21)和式(24)代入，可得一個R-H 邊值問題，它的解為

其中

X1(z)是沿L割開平面上的一單值分支，滿足limz→∞[znX1(z)]=1。

計算式(26)的積分后，再代入到式(21)和式(24)后，即可得到Φ(z)和?(z)。由式(20)和式(23)，可知

這樣，彈性夾雜和準晶基體聲子場的復勢函數形式解已經求得。下面來確定τ0,τ1,··· ,τn和D0這n+2 個常數。由無窮遠處應力為0，得?(∞) = 0，則由式(21)和式(26)可得τ0=0。

由于在無窮遠鄰域內

與式(21)在無窮遠鄰域內的洛朗級數展開式比較可得一方程。

又考慮到Φ(∞)=D0，再代入式(23)，即得

為了保證位移連續條件還需補充在n個端點位移相等的條件。考慮相對剛體位移，可假設其中b1重合，則有n ?1 個獨立條件，由此可得

把式(21)、式(24)代入整理得

其中

至此，Φ(z)和?(z)已完全確定。

4 典型問題

4.1 一條裂紋

如圖2 所示，考慮在準晶基體內任一點z0作用強度為q0的點熱源，在圓形夾雜界面上只含有一條裂紋，不失一般性，假設圓弧裂紋關于x軸對稱，在圓形界面|z| =R上有端點a=Re?iα, b=Reiα。

圖2 一條界面裂紋

4.1.1 溫度場

當n=1, c0=0 時，有

將上述溫度場的解答退化到只含有圓形夾雜時，即α= 0, a=b=R，由式(32)和式(33)得

該結果與文獻[2]一致。

4.1.2 聲子場熱應力

當n=1, τ0=0 時，進而

再把在X1(z)無窮遠鄰域內展開有

與式(29)作比較，即得τ1=A2。

這樣聲子場復勢函數中的待定常數都已求出，聲子場受到的熱應力也可確定。

4.2 兩條裂紋

如圖3 所示，考慮在準晶基體內任意一點z0作用強度為q0的點熱源，在圓形夾雜界面上含有兩條裂紋，不失一般性，假設兩條圓弧裂紋分別關于x軸對稱，在圓形界面|z|=R上有端點a=Re?iα, b=Reiα, a′=Rei(π?α), b′=Rei(π+α)。

圖3 兩條界面裂紋

4.2.1 溫度場

此時c0=0, n=2，得到

4.2.2 聲子場熱應力

當n=2, τ1=0 時，D0=0，可得

4.3 裂紋尖端的熱應力強度因子

為了計算圓弧裂紋端點b的熱應力強度因子，首先做如下坐標變化[26]

這樣在Z平面X軸與裂紋相切與端點b1(Z平面b的對應點)，且

在b1點領域的奇性主項可寫為

聲子場熱應力強度因子定義如下

從式(39)中可以看出，裂紋尖端的聲子場熱應力強度因子與裂紋的大小，熱源的位置及準晶材料的彈性常數有關。

5 數值算例

本節將給出一些數值結果，主要分析夾雜材料為彈性材料，基體材料為十二次對稱二維準晶在點熱源作用下界面弧形裂紋尖端的熱應力及熱應力強度因子的變化趨勢，點熱源位于第四象限。本文考慮一組滿足正交條件的準晶材料參數[28]

取C12= 100 GPa。由于對稱性，只需要計算一個裂紋尖端的熱應力及熱應力強度因子。定義無量綱熱應力及熱應力強度因子為

取α= 30?，得到復平面上圓形夾雜界面上(t=Reiθ(θ= 0,1,··· ,2π))裂紋尖端的熱應力分布曲線圖。

觀察圖4 和圖5 可以看出，裂紋尖端熱應力的大小與外載荷熱源強度呈線性關系，同樣觀察到熱應力在θ= 30?和θ= 330?附近集中。圖6 和圖7 描述了含兩條界面裂紋時，在不同強度的點熱源作用下熱應力隨極坐標轉換過程中θ變化的情形，發現熱源強度越大，裂紋尖端的熱應力越大，觀察到熱應力在θ= 30?, θ= 145?, θ= 215?以及θ=330?附近應力集中。圖8 和圖9 分別給出了只含有一條和兩條界面裂紋時，當點熱源位置位于x軸上時，界面熱應力在不同熱源強度下的變化情形。通過觀察我們發現此時應力峰值相等。

圖4 不同q0 下，邊界上含一條裂紋的徑向應力σr0

圖5 不同q0 下，邊界上含一條裂紋的剪應力σrθ0

圖6 不同q0 下，邊界上含兩條裂紋的徑向應力σr0

圖7 不同q0 下，邊界上含兩條裂紋的剪應力σrθ0

圖8 含一條裂紋邊界上的徑向應力σr0

圖9 含兩條裂紋邊界上的徑向應力σr0

圖10 和圖11 給出了在圓形夾雜大小不同的情況下，無量綱熱應力強度因子K10及K20隨裂紋角度α的變化情形。我們可以觀察到夾雜半徑R的大小對裂紋尖端熱應力強度因子K10有很明顯的影響，R越大熱應力強度因子變化趨勢越明顯并且峰值越高，對于熱應力強度因子K20而言，夾雜半徑越小，熱應力強度因子越小。我們發現熱應力強度因子不會隨裂紋角度α增大而單調增加。如圖10 所示復合材料的熱應力強度因子K10開始隨裂紋角度α增大而增大，在α= 78?附近達到最大值，然后隨裂紋角度進一步增大而減小。而熱應力強度因子K20開始隨裂紋角度α增大而減小，當界面裂紋沿彎曲延伸至57?時產生的熱應力強度因子K20最小，然后又隨裂紋角度α增大而增大。

圖10 不同R 下K10 隨裂紋角度α 的變化

圖11 不同R 下K20 隨裂紋角度α 的變化

6 總結

本文研究了十二次對稱二維準晶在點熱源作用下的共弧圓形夾雜界面裂紋的平面問題，利用Muskhelishvil 復變函數方法，給出了夾雜內外溫度場，聲子場的熱應力分析。當含有一條裂紋以及含有兩條裂紋時，獲得了分區復勢函數和應力場的封閉形式解答。當只含有一條裂紋時求出了裂紋尖端熱應力強度因子的解析表達式，還給出了裂紋尖端熱應力強度因子隨裂紋幾何條件變化而產生的變化情形。觀察到熱應力以及熱應力強度因子的變化受夾雜半徑、熱源強度以及裂紋角度等影響。結果表明在點熱源作用下含界面裂紋圓形夾雜裂紋尖端附近應力集中，應力場呈現振蕩奇異性。本文的解答為準晶材料斷裂性能的研究提供了理論依據。