還原央企高管薪酬真相Restorethetruthofexecutivecompensationofcentralenterprises

《還原央企高管薪酬真相Restorethetruthofexecutivecompensationofcentralenterprises》由會(huì)員分享，可在線閱讀，更多相關(guān)《還原央企高管薪酬真相Restorethetruthofexecutivecompensationofcentralenterprises（99頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、 分類號(hào)：O34 U D C ：080102 密級(jí)： 學(xué)號(hào)：356006511003 南 昌 大 學(xué) 博 士 研 究 生 學(xué) 位 論 文 非橢圓夾雜 Eshelby問題的擴(kuò)展研究 - 光滑夾雜的外場(chǎng)問題和多邊形夾雜的多項(xiàng)式本征應(yīng)變問題 Extended Studies on Eshelby's Problems of Non-elliptical Inclusions - Exterior Fields of Smooth Inclusions and Polynomial Eigenstrain Problems of Polygonal Inclusions

2、 李永剛 培養(yǎng)單位（院、系）： 建筑工程學(xué)院 指導(dǎo)教師姓名、職稱： 鄒文楠 教授 申請(qǐng)學(xué)位的學(xué)科門類： 工學(xué) 學(xué)科專業(yè)名稱： 論文答辯日期： 固體力學(xué) 2016 年 12 月 18 日 答辯委員會(huì)主席： 評(píng)閱人： 2016 年 12 月 22 日 一、學(xué)位論文獨(dú)創(chuàng)性聲明 本人聲明所呈交的學(xué)位論文是本人在導(dǎo)師指導(dǎo)下進(jìn)行的研究工作及取得的 研究成果。據(jù)我所知，除了文中特別加以標(biāo)注和致謝的地方外，論文中不包含 其他人已經(jīng)發(fā)表或撰寫過的研究成果，也不包含為獲得南昌大學(xué)或其他教育機(jī) 構(gòu)的學(xué)位或證書而使用過的材料。與我一同工作的同志對(duì)本研究所做的任何貢 獻(xiàn)均

3、已在論文中作了明確的說明并表示謝意。 學(xué)位論文作者簽名（手寫）： 簽字日期： 年 月 日 二、學(xué)位論文版權(quán)使用授權(quán)書 本學(xué)位論文作者完全了解南昌大學(xué)有關(guān)保留、使用學(xué)位論文的規(guī)定，同意 學(xué)校有權(quán)保留并向國(guó)家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版，允許論文 被查閱和借閱。本人授權(quán)南昌大學(xué)可以將學(xué)位論文的全部或部分內(nèi)容編入有關(guān) 數(shù)據(jù)庫(kù)進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存、匯編本學(xué)位論 文。同時(shí)授權(quán)北京股份有限公司和中國(guó)學(xué)術(shù)期刊（光盤版）電子雜志 社將本學(xué)位論文收錄到《中國(guó)學(xué)位論文全文數(shù)據(jù)庫(kù)》和《中國(guó)優(yōu)秀博碩士學(xué)位 論文全文數(shù)據(jù)庫(kù)》中全文發(fā)表，并通過網(wǎng)絡(luò)向社會(huì)公眾

4、提供信息服務(wù)，同意按 “章程”規(guī)定享受相關(guān)權(quán)益。 學(xué)位論文作者簽名（手寫）： 導(dǎo)師簽名（手寫）： 簽字日期： 年 月 日 簽字日期： 年 月 日 論文題目 非橢圓夾雜 Eshelby 問題的擴(kuò)展研究 姓 名 李永剛 學(xué)號(hào) 356006511003 論文級(jí)別 博士R 碩士£ 院/系/所 建筑工程學(xué)院 專業(yè) 固體力學(xué) E_mail lyg1056@ 備注： □公開 □保密（向校學(xué)位辦申請(qǐng)獲批準(zhǔn)為“保密”， 年 月后公開） 摘要 摘 要 Eshelby 問題的研究歷史已經(jīng)超過半個(gè)世紀(jì)，建立在橢球夾雜 Eshelby 張量

5、 均勻性基礎(chǔ)上的等效夾雜法已經(jīng)成為復(fù)合材料細(xì)觀力學(xué)的基石，由此發(fā)展出多 種用于估計(jì)非均勻材料有效性質(zhì)的方法，如自恰法、Mori-Tanaka 方法、IDD 法 等。然而對(duì)于“均勻性僅限于橢球形狀”的 Eshelby 猜想的證實(shí)或證偽一直懸而 未決，直到 2008 年才有了徹底的證明。在對(duì) Eshelby 猜想反復(fù)辨析的過程中， 非橢球（圓）夾雜問題的研究吸引了眾多研究者的興趣。另一方面，真實(shí)的夾 雜往往是非橢球（圓）的，這一物理背景也進(jìn)一步推動(dòng)了非橢球（圓）夾雜問 題的研究。 三維非橢球夾雜問題的解析研究除了多面體外，由于三維形狀幾何表述的 復(fù)雜性，鮮有涉及其他的形狀；而對(duì)于

6、二維問題，多邊形夾雜和洛朗多項(xiàng)式型 光滑曲線夾雜作為兩類典型的非橢圓夾雜，相應(yīng)的研究雖然都有工作開展，但 仍存在一些問題沒有解決，如洛朗多項(xiàng)式型光滑夾雜 Eshelby 張量場(chǎng)的內(nèi)外完整 性問題，考慮多項(xiàng)式本征應(yīng)變時(shí)解析解的一般性問題等。 本論文從兩個(gè)方面對(duì)非橢圓夾雜問題進(jìn)行擴(kuò)展研究，一是洛朗多項(xiàng)式型光 滑夾雜的外場(chǎng)解問題，二是多邊形夾雜的多項(xiàng)式本征應(yīng)變問題，具體的工作和 得到的結(jié)論如下： （1）從 Eshelby 張量的邊界積分公式出發(fā)，發(fā)展了一種導(dǎo)出洛朗多項(xiàng)式型 光滑夾雜外場(chǎng)解的通用方法，給出了任意旋輪線形和準(zhǔn)平行四邊形夾雜外場(chǎng)的 顯式解，通過數(shù)值計(jì)算并與橢圓等標(biāo)準(zhǔn)模型的

7、對(duì)比研究，得到了夾雜形狀對(duì)外 場(chǎng)的影響范圍。 （2）通過本征位移構(gòu)造了本征應(yīng)變的任意階次多項(xiàng)式形式，基于各向同性 材料的平面彈性復(fù)變函數(shù)方法將任意形狀?yuàn)A雜發(fā)生多項(xiàng)式本征應(yīng)變時(shí)的彈性場(chǎng) 問題轉(zhuǎn)化為基本函數(shù)的邊界積分問題。對(duì)于任意多邊形夾雜，導(dǎo)出了基本函數(shù) 所包含邊界積分的顯式結(jié)果，從而得到了彈性場(chǎng)的解析解。通過數(shù)值計(jì)算得到 了三角形、正方形和多邊形逼近的橢圓夾雜在均勻、線性和二次形式的本征應(yīng) 變作用下的應(yīng)力場(chǎng)和位移場(chǎng)，分析了夾雜的幾何形狀和本征應(yīng)變多項(xiàng)式次數(shù)對(duì) 彈性場(chǎng)的影響。 （3）對(duì)于各向異性磁電彈材料的多項(xiàng)式本征應(yīng)變問題，通過廣義本征位移 摘要 實(shí)現(xiàn)任意階

8、次多項(xiàng)式廣義本征應(yīng)變的構(gòu)造，基于廣義的 Stroh 理論將擾動(dòng)物理場(chǎng) 的求解問題歸結(jié)為兩組基本函數(shù)的邊界積分問題。對(duì)于任意多邊形夾雜，導(dǎo)出 了基本函數(shù)所包含邊界積分的顯式結(jié)果，從而得到了擾動(dòng)物理場(chǎng)的解析解。通 過數(shù)值計(jì)算證實(shí)了“Eshelby 多項(xiàng)式守恒定理”對(duì)于各向異性磁電彈材料橢圓夾 雜的適用性，分析了多邊形頂點(diǎn)處的場(chǎng)量集中和奇性特征，并利用基本函數(shù)的 解析公式進(jìn)行了闡釋。 關(guān)鍵詞：Eshelby 問題；非橢圓夾雜；洛朗多項(xiàng)式型光滑夾雜；外場(chǎng)解；多 項(xiàng)式本征應(yīng)變；多邊形夾雜 Abstract ABSTRACT Studies on the Eshelby’

9、s inclusion problem have continued for over half a century. During this time, the equivalent inclusion method (EIM) based on uniformity of the interior Eshelby tensor of an ellipsoidal inclusion has proved to be the cornerstone of composite micromechanics, and delivered many standard schemes to

10、estimate the efficient properties of heterogeneous materials, such as the self-consistent method, the Mori-Tanaka method, the IDD method, etc. However, until recently, a fundamental problem has been up in the air, that is to verify or falsify the Eshelby’s conjecture, which claims that an ellipso

11、idal inclusion is the only one to possess uniform interior fields. In dedications about this issue, inclusions with non-ellipsoidal/elliptical shapes attract many researchers’ interests. On the other hand, real inclusions in application are usually non-ellipsoidal/elliptical, which is another pr

12、omotion of studying non-ellipsoidal/elliptical inclusion problems. In three-dimensional (3D) problems, except polyhedrons, analytical solutions of non-ellipsoidal inclusions are rarely reported for the complexity of descripting inclusions’ shapes. In two-dimensional (2D) cases, two kinds of typic

13、al non-elliptical inclusions, polygonal ones and those characterized by Laurent polynomials, are touched widely. However, some associated topics are still open, such as the integrity of the Eshelby tensor fields of smooth inclusions characterized by Laurent polynomials, the general analytical so

14、lutions of non-uniform eigenstrain problems in the aforementioned two kinds of 2D inclusions. In this thesis, we dedicated to do some extended research about the non-elliptical inclusion problem from two aspects: (1) the exterior elastic fields of smooth inclusions of Laurent polynomial type; (2

15、) the polynomial eigenstrain problems of polygonal inclusions. The contributions of this dissertation mainly include the following three parts: (1) Starting from the boundary integral formulation of the Eshelby tensor, we developed a general method to explicitly derive the exterior elastic field

16、s of smooth inclusions of Laurent polynomial type, and achieved the close-form solutions of the Abstract (N+1)-gonal hypocycloidal and the quasi-parallelogram shaped inclusions. By comparison with those of some classical and benchmarking models, the circle, the ellipse, etc., the influen

17、ce ranges of the exterior fields from different inclusion’s geometries are analyzed quantitatively. (2) Based on the complex variable method of isotropic elasticity, by constructing the prescribed eigenstrains in polynomial form of arbitrary order through the associated eigendisplacements, the e

18、lastic fields of arbitrary inclusions induced by polynomial eigenstrains are attributed to calculate the boundary integrals involved in a set of basic functions. For arbitrary polygonal inclusions, the involved boundary integrals are explicitly carried out, and the stress and displacement fields

19、of some specific shaped inclusions, the triangle, the square, and the ellipse approximated by N-sided polygons, are numerically achieved. The effects of the shapes of inclusions and the forms of eigenstrains are analyzed by figures. (3) For the polynomial eigenstrain problem in anisotropi

20、c magneto-electro-elastic (MEE) materials, the prescribed extended eigenstains are also constructed through the associated extended eigendisplacements. From the extended Stroh formulism, some tedious derivations finally lead that the induced physical fields depend on two sets of basic functions,

21、 which contain boundary integrals over the transformed inclusion domain. For arbitrary polygonal inclusions, the involved boundary integrals are worked out explicitly. Numerical calculations conforms that the Eshelby’s polynomial conservation theorem is valid for anisotropic MEE materials. Field

22、 concentrations and singularities at vertexes of polygons are analyzed thoroughly, and demonstrated by the explicit forms of basic functions. Keywords: Eshelby problems; non-elliptical inclusions; smooth inclusions of Laurent polynomial type; exterior elastic fields; polynomial eigenstrains; arbi

23、trary polygonal inclusions. 目錄 目 錄 摘 要 ·········································

24、3;················································IV ABSTRA

25、CT ·················································

26、83;·································VI 目 錄 ···············&

27、#183;·················································&

28、#183;······················VIII 第 1 章 緒論·························

29、;··················································

30、;········· 1 1.1 Eshelby 問題概述 ·····································&#

31、183;···························· 1 1.2 國(guó)內(nèi)外研究現(xiàn)狀···················

32、;················································· 2 1

33、.2.1 非橢圓夾雜的外場(chǎng)解···············································

34、3;······ 4 1.2.2 非均勻本征應(yīng)變問題········································

35、83;············· 5 1.3 本文研究的目的和意義·································

36、83;························· 8 1.4 本文的主要內(nèi)容······················

37、·············································· 8 第 2 章 洛朗多項(xiàng)式型光滑夾雜的外彈性

38、場(chǎng)············································10 2.1 Eshelby 張量場(chǎng)的邊界積分表示 ·&#

39、183;··············································10 2.2 洛朗多項(xiàng)式型光滑夾雜的外

40、場(chǎng)解··············································12 2.3 顯式的結(jié)果與算例·

41、··················································

42、·············14 2.3.1 旋輪線形夾雜··································&#

43、183;···························14 2.3.2 準(zhǔn)平行四邊形夾雜···················&#

44、183;····································18 2.3.3 外場(chǎng)解的影響范圍分析··········

45、········································19 2.3.4 G 和G 在夾雜邊界上跳躍性分析 ·····

46、83;·······························24 2 4 2.4 本章小結(jié)···············

47、83;·················································

48、83;··········25 第 3 章 多邊形夾雜的多項(xiàng)式本征應(yīng)變問題： 各向同性材料 ·····················26 3.1 平面彈性問題的復(fù)變函數(shù)表示理論 ·········

49、;··································26 3.2 多項(xiàng)式本征應(yīng)變問題的一般解············&#

50、183;····································27 3.3 多邊形夾雜的顯式解··········

51、83;·················································

52、83;30 3.4 數(shù)值結(jié)果與討論···············································

53、····················32 3.4.1 三角形夾雜···························

54、83;·····································32 3.4.2 正方形夾雜··········

55、··················································

56、·····35 VIII 目錄 3.4.3 多邊形逼近的橢圓夾雜·······································

57、;···········38 3.5 本章小結(jié)·····································

58、;·······································42 第 4 章 多邊形夾雜的多項(xiàng)式本征應(yīng)變問題： 各向異性磁電彈材料 ····

59、3;·······43 4.1 預(yù)備知識(shí)········································

60、3;···································44 4.1.1 各向異性磁電彈問題的控制方程··········

61、3;···························44 4.1.2 二維問題的廣義 Stroh 表示··················&#

62、183;··························46 4.1.3 廣義應(yīng)力和廣義應(yīng)變的對(duì)稱表示···················&#

63、183;··················47 4.2 多項(xiàng)式本征應(yīng)變問題的一般解····························

64、;·····················48 4.3 與格林函數(shù)解之間的關(guān)系·························

65、3;·····························53 4.4 任意多邊形夾雜的顯式解·················

66、83;·····································55 4.5 算例與討論··········&#

67、183;·················································&#

68、183;············58 4.5.1 與現(xiàn)有結(jié)果的比較··································&#

69、183;·····················59 4.5.2 三角形夾雜··························

70、;·······································60 4.5.3 正方形夾雜········&#

71、183;·················································&#

72、183;······63 4.5.4 多邊形逼近的圓形夾雜········································

73、··········66 4.5.5 正多邊形夾雜發(fā)生二次本征應(yīng)變時(shí)的部分場(chǎng)值····················70 4.5.6 多邊形夾雜頂點(diǎn)處的奇性分析···········&

74、#183;·····························73 4.6 本章小結(jié)··················&

75、#183;·················································&

76、#183;·······75 第 5 章 結(jié)論與展望·······································

77、3;··································76 致 謝 ··············&#

78、183;·················································&#

79、183;·························77 參考文獻(xiàn) ·······················

80、;··················································

81、;··············78 附錄 A ··································&#

82、183;·················································&#

83、183;·····87 A.1 從 K-M 雙勢(shì)函數(shù)導(dǎo)出 Eshelby 張量·······································&

84、#183;··87 A.2 Eshenlby 張量在夾雜域內(nèi)的顯式解 ··········································&#

85、183;88 攻讀學(xué)位期間的研究成果 ··············································

86、3;···················91 IX 第 1 章 緒論 第1章 緒論 1.1 Eshelby 問題概述 Eshelby[1, 2]于 1957、1959 年在英國(guó)皇家學(xué)會(huì)會(huì)刊上分別發(fā)表了兩篇論文討 論各向同性無限域中橢球核的彈性場(chǎng)問題，其中第一篇利用格林函數(shù)方法得到 了用橢圓積分表達(dá)的擾動(dòng)彈性場(chǎng)，并得到了一個(gè)重要的結(jié)論，即均勻本征應(yīng)變 引起的橢球內(nèi)部彈性場(chǎng)也是均

87、勻的，同時(shí)發(fā)現(xiàn)如果橢球域內(nèi)和域外的材料性質(zhì) 不同，這時(shí)在均勻的遠(yuǎn)場(chǎng)加載下橢球內(nèi)的擾動(dòng)彈性場(chǎng)同樣為常值，可以通過“等 效夾雜”的方法來實(shí)現(xiàn)求解。Eshelby 的第二篇論文更加詳細(xì)地研究了橢球核的 外部彈性場(chǎng)，構(gòu)造了完全用橢球簡(jiǎn)諧勢(shì)表達(dá)的外彈性場(chǎng)封閉解。兩年后，Eshelby [3] 再次發(fā)文并斷言 ，在均勻本征應(yīng)變的作用下，只有橢球才具有內(nèi)部彈性場(chǎng)的均 勻性，該結(jié)論被后人稱為 Eshelby 猜想，同時(shí)分析指出，如果與應(yīng)力無關(guān)的相變 x 、x 、x 3 N 的 次多項(xiàng)式形式給出，則由此引起的橢球內(nèi) 位移（本征位移）以坐標(biāo) 1 2 N 部的擾動(dòng)位移場(chǎng)也為 次多項(xiàng)式形

88、式，這一結(jié)論被后人稱為 Eshelby 多項(xiàng)式守恒 定理。Eshelby 關(guān)于橢球核的一系列工作奠定了復(fù)合材料細(xì)觀力學(xué)的基礎(chǔ)，后來 由此發(fā)展了多種細(xì)觀力學(xué)方法（自恰法、廣義自恰法、Mori-Tanaka 方法、IDD 法等）用于估計(jì)復(fù)合材料、多晶材料等各種非均勻材料的有效性質(zhì)[4, 5]，給橢球 夾雜的解帶來了廣泛的應(yīng)用，比如復(fù)合材料有效性質(zhì)的估計(jì)、壓電/壓磁器件的 性能設(shè)計(jì)、應(yīng)變誘導(dǎo)的量子點(diǎn)/線的生長(zhǎng)等[6-9] 。 Eshelby 的工作無疑是開創(chuàng)性的，但在物理概念的使用上略顯局限，稱夾雜 內(nèi)部發(fā)生的這種與應(yīng)力無關(guān)的應(yīng)變?yōu)橄嘧儜?yīng)變（stress-free transform

89、ation strain）。 [10] “本征應(yīng)變”的概念是后來 Mura 提出的，用來泛指在彈性體內(nèi)的一個(gè)局部區(qū) 域由于某種原因在無約束情況下產(chǎn)生的一類永久的非彈性變形，比如熱膨脹應(yīng) 變（thermal dilatation strains）、濕熱膨脹應(yīng)變（hydrothermal dilatation strains）、 晶格參數(shù)錯(cuò)配引起的應(yīng)變（misfit strains）、相變應(yīng)變（phase transformation strains）、 塑性應(yīng)變（plastic strains）、磁-力或電-力耦合引起的應(yīng)變（magneto-elastic or electro-e

90、lastic strains）。Mura 同時(shí)發(fā)展了一套方法來處理與橢球核有關(guān)的問題。 [1] 按照 Eshelby 的描述 ，在無限大空間中有一橢球子域，先假設(shè)橢球子域從 周圍基體中切割出來，然后子域內(nèi)部發(fā)生均勻本征應(yīng)變從而產(chǎn)生變形，接著在 子域邊界上施加面力使其回到原來的形狀，再將其放回到原來的基體中去，同 1 第 1 章 緒論 時(shí)將施加的面力去掉，這時(shí)由于基體材料的存在，在橢球子域內(nèi)部必將產(chǎn)生一 部分附加變形，即擾動(dòng)應(yīng)變，在基體中也將產(chǎn)生相應(yīng)的變形，最終達(dá)到彈性平 衡狀態(tài)。如果子域與基體的材料相同，子域內(nèi)發(fā)生本征應(yīng)變，這時(shí)稱為夾雜問 題（inclusio

91、n problem），也就是通常所講的 Eshelby 問題，或第一類 Eshelby 問 題；如果子域內(nèi)材料與基體材料不同，稱為異質(zhì)問題（inhomogeneity problem）， [11] 又稱為第二類 Eshelby 問題 。 Eshelby 問題已經(jīng)成為線彈性材料的一個(gè)經(jīng)典問題，在過去的幾十年里，眾 多的研究者開展了大量相關(guān)的研究工作，從各向同性到各向異性，從空間問題 到平面問題，從無限域到有限域[12-16]，從純彈性到彈性與熱、電、磁的耦合， 從靜態(tài)到動(dòng)態(tài)[17-19]，從橢球（圓）到非橢球（圓），從均勻本征應(yīng)變問題到考 [10] [4] 慮非均勻的本征應(yīng)變

92、等等，相應(yīng)的文獻(xiàn)可參考 Mura 、Nemat-Nasser & Hori 、 [7] [20] [21] [22] Li & Wang 、Nomura 、張研 & 韓林 和 Li & Gao 等人的專著，以及 Huang [23] [24] [25] [6] [26] [8] & Mura 、Mura 、Mura 、Gutkin 、Ovid'ko & Sheinerman 和 Zhou et al. 等人的綜述性文章。近二十年來，人們發(fā)現(xiàn) Eshelby 夾雜問題在各種納米材料和 納米結(jié)構(gòu)（如量子點(diǎn)（QDs）

93、、量子線（QWRs）等）的建模和分析中也具有廣 泛的應(yīng)用，這方面的內(nèi)容可以參考 Maranganti & Sharma[27]和 Duan et al. 的文 章及他們?cè)谖闹兴械膮⒖嘉墨I(xiàn)。Eshelby 問題中所涉及的 Eshelby 張量、Hill 張量、Moment 張量和 Concentration 張量在牛頓勢(shì)問題和線性彈性靜力學(xué)問題中 [28] [29] 具有廣泛的應(yīng)用，最近 Parnell 就此方面對(duì)橢球異質(zhì)問題進(jìn)行了綜述。 1.2 國(guó)內(nèi)外研究現(xiàn)狀 Eshelby 關(guān)于夾雜問題的工作有兩個(gè)先決條件，一是夾雜形狀要為橢球，二 是夾雜內(nèi)的本征應(yīng)變要是均勻的。

94、如前所述，橢球夾雜的解已成為復(fù)合材料細(xì) 觀力學(xué)的重要基礎(chǔ)[4, 10, 30]，并在斷裂力學(xué)等其他領(lǐng)域帶來了廣泛的應(yīng)用，盡管如 此，在過去幾十年里，人們還是嘗試著剔除掉 Eshelby 的兩個(gè)先決條件從而將問 題推廣到更一般的情況，即：（1）夾雜形狀為非橢球；（2）夾雜內(nèi)的本征應(yīng) 變?yōu)榉蔷鶆虻?，圍繞這兩個(gè)方面開展了大量的工作[31-36] 。 另一方面，對(duì) Eshelby 猜想的證實(shí)或證偽問題數(shù)十年來一直處于懸而未決的 困境，直到 2008 年才得到徹底的證明[37-39]。在對(duì) Eshelby 猜想的反復(fù)辨析中， 大量的研究工作集中在各種非橢球（圓）夾雜的彈性場(chǎng)問題上[40-

95、47]。在實(shí)際應(yīng) 2 第 1 章 緒論 用方面，不論是天然的非均勻材料，還是各種人造的復(fù)合材料，它們所含顆粒 的真實(shí)形狀通常都是非橢球（圓）的，這一背景成為對(duì)非橢球（圓）夾雜問題 進(jìn)行研究的另一個(gè)重要驅(qū)動(dòng)[48, 49] 。 在各類 Eshelby 問題的研究中，格林函數(shù)方法發(fā)揮了重要的作用。Eshelby[1] 通過對(duì)點(diǎn)力引起的位移格林函數(shù)G 在橢球域內(nèi)的積分得到了橢球夾雜的 ij Eshelby 張量，并得到了“橢球域內(nèi)的應(yīng)變場(chǎng)是均勻的”這一重要結(jié)論。Wang[50] 利用格林函數(shù)方法研究了三維壓電材料橢球夾雜的 Eshelby 問題。除橢球外，格 [

96、51] 林函數(shù)方法也被用來求解其他形狀?yuàn)A雜的彈性場(chǎng)問題，如 Chiu 分析了長(zhǎng)方體 夾雜問題，Wu & Du[52], [53]研究了圓柱形夾雜的應(yīng)力場(chǎng)問題，Wu & Du 得到了 半球夾雜引起的彈性場(chǎng)。在二維問題中，Ma[55] 導(dǎo)出了相變應(yīng)變核引起的 Muskhelishvili 勢(shì)函數(shù)的基本解，并利用格林函數(shù)方法構(gòu)造了平面夾雜問題的積 分解，分析了由馬氏體相變和/或鐵彈相變引起的相變?cè)鲰g現(xiàn)象。在另一篇文章 [54] [56] 中，Ma et al. 按照同樣的思路并利用 Stroh 理論分析了各向異性材料的相變?cè)? 韌現(xiàn)象。 在非橢球（圓）夾雜問題的眾多

97、研究工作中，Zou et al.的工作是比較系統(tǒng)性 的。Zou et al.[47], [57]基于張量的不可約分解系統(tǒng)地研究了非橢圓夾雜第一類 Eshelby 問題，得到了二維彈性任意非橢圓夾雜的解析解，并得到如下重要結(jié)論： 對(duì)凸形夾雜采用橢圓夾雜 Eshelby 張量來近似平均 Eshelby 張量是可接受的，但 對(duì)非凸形夾雜則未必；而對(duì)在細(xì)觀力學(xué)的非橢圓（無論是凸的還是非凸的）顆 [58] 粒問題中引入的廣義 Eshelby 張量 ，無論是取為橢圓夾雜 Eshelby 張量還是平 均 Eshelby 張量都是不可接受的。二維任意形狀熱夾雜問題、二維各向異性壓電 材料任意形

98、狀?yuàn)A雜問題、二維多鐵雙材料（multiferroic bimaterials）任意形狀?yuàn)A 雜問題也都得到解決[59-61]，都得到了相應(yīng)問題的解析解。針對(duì)有限域問題，Zou [62] et al. 提出了一種基于疊加原理的分析方法，并將其應(yīng)用于平面有限域夾雜問 題的求解，得到了圓形基體內(nèi)含有正方形及偏心圓形夾雜時(shí)的解析解，同時(shí)將 這一方法推廣到二維有限域內(nèi)熱夾雜問題和反平面問題的求解[63, 64]。Zou & Zheng[65]指出求解含有異質(zhì)的材料在遠(yuǎn)場(chǎng)加載下的彈性平衡問題和第二類 Eshelby 問題是等價(jià)的。值得一提的是，Zou et al.關(guān)于非橢圓夾雜的一系列

99、工作 都最終實(shí)現(xiàn)了問題的解析求解，而且在導(dǎo)出過程中沒有利用一般的格林函數(shù)方 法，從而避免了在夾雜內(nèi)部可能會(huì)出現(xiàn)的奇異性。 3 第 1 章 緒論 1.2.1 非橢圓夾雜的外場(chǎng)解 對(duì)于橢球夾雜，Eshelby 早在 1959 年就得到了夾雜的外場(chǎng)解，后來 Mura 等 人[66, 67]對(duì)這一問題進(jìn)行了進(jìn)一步研究，剔除了 Eshelby 假設(shè)本征應(yīng)變?cè)跈E球夾雜 內(nèi)均勻分布的先決條件，得到了用單位球域內(nèi)的積分表達(dá)的外場(chǎng)解，并將結(jié)果 [68] 推廣到各向異性的情況。然而，多年后 Ju & Sun 仍然認(rèn)為前人得到的橢球夾 雜的外場(chǎng)解是不完善的，他們?cè)跈E球夾雜

100、的外部引入一個(gè)虛擬橢球，并利用該 橢球的外法線單位矢量導(dǎo)出了另一種外場(chǎng) Eshelby 張量的完全解，相比 Mura 等 [68] 人的結(jié)果，這個(gè)解更簡(jiǎn)單、更完備，而且物理意義更明確。利用 Ju & Sun 導(dǎo) [69] 出的關(guān)于橢球夾雜外場(chǎng) Eshelby 張量的一般公式，Jin et al. 得到了橢圓柱夾雜 [70] 的外場(chǎng)顯式解，其中還指出了 Kim & Lee 解答的錯(cuò)誤在于不滿足內(nèi)對(duì)稱性。 對(duì)于涉及非橢球夾雜的 Eshelby 問題，對(duì)夾雜邊界進(jìn)行恰當(dāng)?shù)拿枋鲆恢笔乔? 解與分析的一個(gè)關(guān)鍵問題。實(shí)際的情況是，夾雜的邊界可以完全由光滑的曲線 （面）構(gòu)成

101、，可以都是直線段（平面）（多邊形/多面體），也可以既包含光滑 曲線（面）又包含直線段（平面）。由于幾何描述的復(fù)雜性，三維夾雜問題中 除了多面體和橢球外其他形狀的夾雜鮮有報(bào)道，對(duì)于二維問題，除了多邊形和 橢圓以外，用洛朗多項(xiàng)式表達(dá)的光滑曲線可以描述一大類非橢圓夾雜，而且數(shù) 學(xué)處理上也相對(duì)簡(jiǎn)單，因此吸引了不少研究者的興趣[46, 47]。對(duì)于多邊形（多面 體）夾雜，Eshelby 張量的解析解通常會(huì)表達(dá)為各段（片）邊界求和的形式，而 且在夾雜內(nèi)部和夾雜外部具有統(tǒng)一的形式[41, 71-73]。但是，對(duì)于光滑曲線（面） 夾雜，外場(chǎng)與內(nèi)場(chǎng)的推導(dǎo)存在著較大的差異，而且外場(chǎng)的分析更為復(fù)雜[2

102、, 10, 46, 68] 實(shí)際上就我們所知，由光滑曲線圍成的非橢圓夾雜 Eshelby 問題只能在文獻(xiàn)中找 到屈指可數(shù)的結(jié)果，而且有的不是完全解析的，有的給出的解答并不完備。除 了前面提到的橢圓和橢球夾雜，Onaka et al.[49, 74-77]研究了甜甜圈形和超球形的夾 。 [46] 雜問題，通過數(shù)值的方法得到了 Eshelby 張量在夾雜域內(nèi)的平均值。Ru 基于 復(fù)變函數(shù)方法并利用解析延拓及保形映射等技巧發(fā)展了一種適用于求解任意形 狀?yuàn)A雜彈性場(chǎng)的解析方法，可惜的是他僅僅給出了橢圓夾雜的解。Zou et al.[47] 基于張量的不可約分解導(dǎo)出了全新的 Eshelb

103、y 張量場(chǎng)的邊界積分公式，并且顯式 地得到了若干洛朗多項(xiàng)式型非橢圓夾雜的內(nèi)部彈性場(chǎng)的解析解以及它們的域內(nèi) 平均值，然而遺憾的是它們所對(duì)應(yīng)的外場(chǎng)解并沒有相應(yīng)的給出。 [2] 就如 Eshelby 所指出的 ，“夾雜內(nèi)的彈性場(chǎng)可以單獨(dú)地計(jì)算，而不必先找 到夾雜外部的彈性場(chǎng)”，而且“大量的信息可以由內(nèi)部彈性場(chǎng)的解答中單獨(dú)得 4 第 1 章 緒論 到”。然而，當(dāng)我們關(guān)注這幾個(gè)方面：（1）Eshelby 張量場(chǎng)的完整性；（2）夾 雜之間的相互作用；（3）夾雜特性對(duì)基體中的影響范圍；（4）夾雜與夾雜之 間、夾雜與有限域邊界之間以及夾雜與其他類型的材料缺陷之間的相互作用，

104、 這時(shí)，夾雜的外部彈性場(chǎng)就變得至關(guān)重要了。因此，作為非橢圓夾雜 Eshelby 問 題擴(kuò)展研究的第一個(gè)內(nèi)容，在本論文中我們首先對(duì)洛朗多項(xiàng)式型光滑夾雜的外 部彈性場(chǎng)問題進(jìn)行研究。 1.2.2 非均勻本征應(yīng)變問題 一直以來，Eshelby 夾雜問題的研究主要考慮的是均勻本征應(yīng)變，涉及非均 勻本征應(yīng)變的 Eshelby 問題，由于其數(shù)學(xué)推演的復(fù)雜性，無論是各向同性材料， 還是各向異性材料，相關(guān)的研究都不多。然而在實(shí)際問題中，非均勻的本征應(yīng) 變很多情況下更接近物理真實(shí)，用均勻本征應(yīng)變?nèi)ツM只能算作一定程度上的 近似，如電子芯片常用的各種場(chǎng)效應(yīng)管（FETs），在受到高度局部化的點(diǎn)狀熱

105、源的瞬時(shí)加熱時(shí)，這時(shí)的溫度場(chǎng)是高度不均勻的，可以用高斯形式或指數(shù)形式 的非均勻本征應(yīng)變場(chǎng)來描述[78]。如果本征應(yīng)變來源于材料的擴(kuò)散以及材料組分 的改變等，這時(shí)本征應(yīng)變應(yīng)該遵循擴(kuò)散微分方程，也將導(dǎo)致其呈非均勻分布。 如果考慮本征應(yīng)變的形成過程，即便最終的狀態(tài)是均勻的，但在某一瞬時(shí)則可 能是非均勻的。另外在處理非橢球的異質(zhì)問題時(shí)，非均勻的本征應(yīng)變，特別是 [41] 多項(xiàng)式形式的本征應(yīng)變，可以發(fā)揮重要的作用 。 [3] Eshelby 指出如果橢球夾雜內(nèi)的本征應(yīng)變以任意階次多項(xiàng)式給出，則夾雜 域內(nèi)的擾動(dòng)應(yīng)變場(chǎng)也必將是同階多項(xiàng)式形式，然而這一多項(xiàng)式具體的顯式表達(dá) [79] 直到

106、 2002 年 Rahman 才最終給出，并將這一結(jié)論稱為 Eshelby 多項(xiàng)式守恒定理。 [80] 繼 Eshelby 之后，Moschovidis 在他的博士論文中研究了橢球核內(nèi)多項(xiàng)式形式 的本征應(yīng)變問題，Mura[10]提出了一個(gè)基于多極展開的一般性方法，用以量化在 多項(xiàng)式本征應(yīng)變下橢球夾雜內(nèi)部的應(yīng)變場(chǎng)，然而實(shí)際上除了初始的少數(shù)幾項(xiàng)外， [79] 徹底完成所有項(xiàng)的推演有相當(dāng)難度。Rahman 通過勢(shì)積分來表達(dá)這一多項(xiàng)式， 而這些勢(shì)積分能夠進(jìn)一步轉(zhuǎn)換成用任意階多項(xiàng)式表達(dá)的代數(shù)形式，而且特別適 合用計(jì)算機(jī)做符號(hào)運(yùn)算。橢球夾雜多項(xiàng)式本征應(yīng)變問題的研究意義非凡，不僅 具有數(shù)學(xué)

107、演繹上的自然之美，也具有相當(dāng)重要的應(yīng)用價(jià)值，拓寬了人們處理橢 球夾雜非均勻本征應(yīng)變問題的范圍[79]。在許多實(shí)際的問題中，本征應(yīng)變可能不 會(huì)顯式地以多項(xiàng)式形式出現(xiàn)，但是只要表征本征應(yīng)變的函數(shù)是連續(xù)的，就可以 5 第 1 章 緒論 根據(jù) Bernstein’s 定理，在橢球域內(nèi)用多項(xiàng)式來做近似。在這種近似的具體實(shí)現(xiàn) 過程中，可以借助最小二乘法、Marquardt-Levenberg 方法等一些完善的算法。 橢球夾雜多項(xiàng)式本征應(yīng)變問題的解也被用來處理動(dòng)態(tài)的橢球夾雜 Eshelby 問題 [17-19, 81-84]。 最近幾十年，納米科技的發(fā)展為非均勻本征應(yīng)變問題的

108、研究提供了更加廣 闊的物理背景。例如，當(dāng)靜電場(chǎng)、靜磁場(chǎng)或純彈性場(chǎng)中的橢球顆粒受到非均勻 的本征場(chǎng)時(shí)，如何計(jì)算作用在其上的力或力矩[85]？這一問題有大量的應(yīng)用，比 如設(shè)計(jì)磁性納米鑷子（magnetic nanotweezers）時(shí)，為了能夠精確地控制和操縱 納米顆粒，關(guān)鍵的問題就是如何確定納米顆粒所受的非均勻本征場(chǎng)與作用在納 米顆粒上的力和力矩之間的關(guān)系[86, 87]。另外，在研究某種溶液中顆粒積聚或分 離的機(jī)制時(shí)，為了精確考量相鄰顆粒之間的相互作用力，也必須假設(shè)本征物理 場(chǎng)是非均勻的。 [88] [89] Shodja & Shokrolahi-Zadeh 和

109、Avazmohammadi et al. 分別研究了雙材料 平面內(nèi)包含任意取向橢圓夾雜的問題，假設(shè)夾雜域內(nèi)的本征應(yīng)變?yōu)槎囗?xiàng)式形式。 后來發(fā)現(xiàn)多項(xiàng)式本征應(yīng)變可以用來處理異質(zhì)之間的相互作用[90, 91]。Shodja & [92] Sarvestani 發(fā)現(xiàn)在處理含有涂層的異質(zhì)問題時(shí)也可以用多項(xiàng)式形式的本征應(yīng)變 [93] 來處理。Rodin & Hwang 系統(tǒng)地研究了多項(xiàng)式形式的本征應(yīng)變?cè)谔幚矶鄠€(gè)橢球 [94] [85] 異質(zhì)相互作用時(shí)的應(yīng)用。Taya 和 Liu 利用數(shù)學(xué)工具嚴(yán)格地證明了橢圓夾雜內(nèi) 多項(xiàng)式本征應(yīng)變引起的擾動(dòng)應(yīng)變場(chǎng)也是同階多項(xiàng)式形式，同時(shí)將擾

110、動(dòng)應(yīng)變場(chǎng)多 項(xiàng)式表達(dá)式中的系數(shù)通過橢圓積分顯式地給出，嚴(yán)格證明了 Eshelby 多項(xiàng)式守恒 定理。Chen[95]2014 年應(yīng)用復(fù)變函數(shù)方法和保形映射技術(shù)求解了平面橢圓夾雜發(fā) 生多項(xiàng)式形式本征應(yīng)變時(shí)的彈性平衡問題。 關(guān)于橢球（圓）夾雜非均勻本征應(yīng)變問題的有限工作大都集中于多項(xiàng)式形 式的本征應(yīng)變，除此之外，其他形式分布的非均勻本征應(yīng)變也有涉及。Sharma & Sharma[78]2003 年研究了三維橢球夾雜（ W : x2 / a2 + x2 / a2 + x2 / a2 £ 1 ）發(fā)生指 1 1 2 2 3 3 數(shù)形式或高斯形式本征應(yīng)變時(shí)

111、的彈性平衡問題，但僅限于均勻膨脹型的本征應(yīng) e0 = e0d 變 （這種本征應(yīng)變常見于電子芯片中高度局部化的點(diǎn)狀熱源等）。 ij ij 對(duì)各向異性材料的非均勻本征應(yīng)變問題只有少量的工作開展。Asaro & Barnett[96]利用一般各向異性無限域的格林函數(shù)分析了橢球夾雜內(nèi)發(fā)生M 次多項(xiàng)式 形式的本征應(yīng)變時(shí)，在橢球內(nèi)部引起的擾動(dòng)應(yīng)變場(chǎng)和應(yīng)力場(chǎng)也是同樣次數(shù)的多項(xiàng)式 [97] [98] 形式。Kinoshita & Mura 和 Mura & Kinoshita 研究了同樣的問題，并給出了 6 第 1 章 緒論 [96] 相應(yīng)的

112、位移場(chǎng)，進(jìn)一步完善了 Asaro & Barnett 的結(jié)果。在研究異質(zhì)問題時(shí)，多 項(xiàng)式本征應(yīng)變有重要應(yīng)用。Nie et al.[33, 99, 100]求得了正交各向異性材料中橢圓夾 雜發(fā)生多項(xiàng)式形式的本征應(yīng)變時(shí)的解析解，其中假設(shè)夾雜域內(nèi)的彈性場(chǎng)也為多 項(xiàng)式形式。Huang et al.[101]分析了一般各向異性材料橢圓夾雜的非均勻本征應(yīng)變 問題。 [65] Zou & Zheng 指出求解非均勻材料有效性質(zhì)的遠(yuǎn)場(chǎng)問題與第二類 Eshelby 問題是等價(jià)的，對(duì)于非橢球（圓）的異質(zhì)顆粒，如果將本征應(yīng)變假設(shè)為多項(xiàng)式 形式，就可以實(shí)現(xiàn)與第一類 Eshelby 問題

113、的等價(jià)變換，如圖 1.1 所示，通過積分 方程解出待定系數(shù)，實(shí)現(xiàn)非橢球（圓）異質(zhì)問題的求解。簡(jiǎn)言之，考慮多項(xiàng)式 本征應(yīng)變的 Eshelby 夾雜問題是非橢球（圓）異質(zhì)問題求解的一條重要途徑。 圖 1.1 異質(zhì)問題與夾雜問題的轉(zhuǎn)換關(guān)系。 以上關(guān)于橢球（圓）夾雜多項(xiàng)式本征應(yīng)變的工作從不同的角度證實(shí)了 Eshelby 多項(xiàng)式守恒定理對(duì)各向同性材料和各向異性材料都是適用的。而對(duì)于非橢球 （圓）夾雜，多項(xiàng)式本征應(yīng)變問題的研究雖然具有重要的理論和應(yīng)用價(jià)值，但 相關(guān)的研究并不很多。Cheng et al.[102]利用格林函數(shù)研究了各向同性矩形夾雜內(nèi) 發(fā)生四次多項(xiàng)式形式本征應(yīng)變時(shí)的彈性場(chǎng)問題

114、，發(fā)現(xiàn)在夾雜內(nèi)部的彈性場(chǎng)存在 奇性，而在夾雜外部不存在奇性。各向異性材料非橢圓夾雜多項(xiàng)式本征應(yīng)變問 題的研究最近幾年才有報(bào)道。Sun et al.[103]研究了包含任意多邊形夾雜的二維各 向異性壓電材料的線性本征應(yīng)變問題，給出了彈性和壓電場(chǎng)的封閉解答。Chen et al.[104]將問題推廣到半平面情況，Yue et al.[105]最近進(jìn)一步研究了二次多項(xiàng)式形式 的本征應(yīng)變問題。以上研究代表了非均勻本征應(yīng)變問題的最新成果，然而在研 究中廣泛利用了平面問題（包括各向同性材料、壓電材料、磁電彈材料）的格 林函數(shù)[106-109]，由于格林函數(shù)本身在夾雜域內(nèi)的奇性特征，隨著本征應(yīng)

115、變多項(xiàng)式 次數(shù)的增加，問題解答的復(fù)雜程度會(huì)極大增加。因此，發(fā)展一種一般性的理論 和方法來處理任意形狀?yuàn)A雜的任意階次多項(xiàng)式本征應(yīng)變問題變得重要而迫切。 鑒于以上分析，作為非橢圓夾雜 Eshelby 問題擴(kuò)展研究的第二個(gè)內(nèi)容，本論文將 7 第 1 章 緒論 分別以各向同性材料和各向異性磁電彈材料為背景，就任意多邊形夾雜的多項(xiàng) 式本征應(yīng)變問題進(jìn)行探索研究。 1.3 本文研究的目的和意義 如前所述，非橢圓夾雜 Eshelby 問題的研究仍有很多問題是開放的，本文將 從兩個(gè)方面對(duì)非橢圓夾雜 Eshelby 問題進(jìn)行擴(kuò)展研究，一是洛朗多項(xiàng)式型光滑夾 雜的外場(chǎng)解問題，二是

116、多邊形夾雜的多項(xiàng)式本征應(yīng)變問題，具體的研究目的有 三： 1、 發(fā)展一種求解洛朗多項(xiàng)式型光滑夾雜外場(chǎng)解析解的通用方法，得到 若干該類型夾雜外場(chǎng)的顯式解，分析并得到夾雜形狀對(duì)外場(chǎng)的影響 范圍； 2、 3、 以現(xiàn)有二維非橢圓夾雜均勻本征應(yīng)變問題的理論和方法為基礎(chǔ)，發(fā) 展一套適用于多項(xiàng)式本征應(yīng)變問題的新的求解方法； 對(duì)于各向同性材料和各向異性磁電彈材料，分別得到二維非橢圓（多 邊形）夾雜發(fā)生任意多項(xiàng)式形式的本征應(yīng)變時(shí)擾動(dòng)物理場(chǎng)的解析 解，分析得到夾雜形狀和本征應(yīng)變多項(xiàng)式次數(shù)對(duì)擾動(dòng)物理場(chǎng)的影響 機(jī)制。 本文研究洛朗多項(xiàng)式型光滑夾雜的外場(chǎng)解析解，將使該類非橢圓夾雜的彈 性場(chǎng)變得完整，拓寬了非橢圓夾雜相關(guān)問題的研究領(lǐng)域；通過對(duì)多邊形夾雜在 多項(xiàng)式本征應(yīng)變條件下擾動(dòng)物理場(chǎng)的解析研究，探討微結(jié)構(gòu)幾何參數(shù)（大小、 形狀、取向等）對(duì)擾動(dòng)物理場(chǎng)的影響機(jī)制，為非橢圓異質(zhì)問題的解析求解，基 體-顆粒型復(fù)合材料有效性質(zhì)的合理估計(jì)和性能設(shè)計(jì)等細(xì)觀力學(xué)問題奠定更堅(jiān)實(shí) 的理論基礎(chǔ)。 1.4 本文的主要內(nèi)容 本文圍繞二維非橢圓夾雜 Eshelby 問題進(jìn)行擴(kuò)展研究，以解析研究為特色。 各章的研究?jī)?nèi)容如下： 第 2 章從 Eshelby 張量的邊界積分公式出發(fā)，發(fā)展了一種導(dǎo)出洛朗多項(xiàng)式型 光滑夾雜外場(chǎng)解的通用方法，給出了任意旋輪線形和準(zhǔn)平行四邊形夾雜外場(chǎng)的 顯式解，通