【正文】 出發(fā)，利用半逆解法，假設材料的所有電彈性常數(shù)沿厚度方向按同一函數(shù)規(guī)律變化，研 究功能梯度壓電懸臂梁在上表面受均布荷載、自由端受集中力和集中力矩聯(lián)合作用問題的力電耦合場，求得了相應的應力和位移。 給出數(shù)值算例，假設梯度函數(shù)為指數(shù)形式，驗證功能梯度壓電材料懸臂梁簡化理論，分析不同的梯度變化對功能梯度壓電梁靜力響應的影響。 畢業(yè)論文 第 7 頁 共 35 頁 2 功能梯度材料梁、板、殼結構的分析方法 梁、板、殼以及它們的組合是目前功能梯度材料的常見結構形式， 分析方法主要以簡化的結構理論為主，也有一些是從三維方程出發(fā)進行直接求解 。 由于以往針對均勻材料建立的經典梁、板、殼結構理論并不能完全適用于功能梯度材 料，因此必須發(fā)展針對功能梯度材料特點的梁、板、殼結構分析方法，尋找功能梯度材料梁、板、殼結構的解析解、半解析解以及簡化理論解，下面將分別闡述 。 解析解 對于功能梯度材料梁、板、殼結構，由于控制方程通常為變系數(shù)的微分方程，其邊值問題的解析求解難度非常大，因此在現(xiàn)有文獻中，解析解非常少 。 解析解的價值在于它們能為其他理論模型或計算方案提供假設的依據(jù)和驗證的考題 。 Sankar 針對材料常數(shù)沿厚度為指數(shù)分布的正交各向異性功能梯度梁受任意垂直載荷作用的情形推導得到了彈性力學精確解 。 稍后 Sanka 等給出了考慮溫度效 應的一個精確解 。 這些工作都假設材料參數(shù)沿厚度以同一指數(shù)變化 。 對于沿厚度非指數(shù)變化的情況，還沒有彈性力學精確解的報道 。 在解析解方面， Lekhnitskii 在其名著中給出了正交各向異性懸臂梁端部受剪 、彎的解，但其后近 40 年幾 乎沒有任何進展 。 最近 ，仲政等和丁皓江等系統(tǒng)發(fā)展了考慮力、熱、電、磁等作用的各向異性功能梯度材料直梁二維問題的廣義應力函數(shù)解法，獲得的解析解適用于各材料參數(shù)沿厚度方向的任意梯度分布情況，可在 SaintVenant 意義上考慮簡支、固支等不同的位移約束邊界條件和集中力、彎矩等端部載荷作用 。 這一方法將 Lekhnitskii 在各向異性平面彈性問題上的研究推進了一大步，豐富了非均勻材料力學的求解手段，具有重要的科學意義 。Huang 等進一步研究了功能梯度壓電執(zhí)行器的響應特性，給出了解析表達式并與一維梁理論、有限元和實驗結果等進行了比較，說明了解析解的有效性及一維梁理論解的不足之處 。 Jiang 等則發(fā)展了基于調和函數(shù)表示的通解解法，求得了密度功能梯度材料懸臂梁的兩個解析解 。 智能功能梯度結構的響應特性近 年得到了很多研究者的關注 。 已有的研究大都假設壓電層和功能梯度材料主體結構之問粘結完好 。 計算發(fā)現(xiàn)，由于界 面的弱化，使得結構系統(tǒng)的靜動力特性發(fā)生了一定程度的改變，從而可能導致作動器或者傳感器性能的變化，必須加以重視 [6]。 陳偉球等將辛空間 的彈性力學解法推廣應用到功能梯度材料結構的解析研究中，提畢業(yè)論文 第 8 頁 共 35 頁 出了移位 Hamilton 矩陣的概念， 并且 發(fā)現(xiàn)材料的功能梯度屬性將改變有關特征函數(shù)的形式，對不同的 特征函數(shù)建立了 辛正交關系，并提出了穩(wěn)定的數(shù)值計算方法 。 最近，他們又將工作進一步推廣到功能梯度壓電材料，詳細研究了非均勻參數(shù)對 SaintVenant 解的影響 。 在軸對稱載荷作用下的圓 /環(huán)板的彎曲問題是彈性力學的一個經典問題 。 Han 等假設彈性模量沿徑向按冪律分布，考察了功能梯度旋轉圓盤的一維軸對稱變形問題 。 Li等將 Ding 等提出的直梁的分析方法進一步推廣應用于功能梯度彈性圓板 。 對于橫觀各向同性功能梯度圓板的純彎曲問題，首次給出了解析形式的彈性力學解 。 與直梁問題類似，這一方法可以考察材料常數(shù)沿厚度方向的任意變化以及圓周邊界上的簡支或固支邊界條件 。 Li 等通過在位移表達式中引入適當?shù)膶?shù)函數(shù)項，獲得了內外圓周邊界條件任意組合時在均布載荷作用下橫觀各向同性功能梯度圓環(huán)板的解析解，該解可以考慮材料常數(shù)沿厚度方向的任意分布 。 對于均勻 電勢作用下的功能梯度圓板和均布載荷作用下的功能梯度壓電壓磁板，同樣通過假設位移、電勢和 (或 )磁勢的適當形式，導出了相應的軸對稱解析解 [24]。 針對沿厚度方向材料常數(shù)任意變化的情形， Kaprielian 等及 Mian 等提出了一種非常有效的求解方法，即通過假設位移的適當形式將三維問題化為二維問題以簡化求解過程 。 但與經典的板殼理論不同，在位移模式中含有很多與厚度坐標有關的未知函數(shù)，因此可以考慮材料特性沿板厚的任意變化 。 但他們的研究僅局限于各向同性材料，且板的上下表面沒有載荷作用 。 Yang 等推廣了非均勻材料 板的 Spence 分析方法，并求解了均布載荷作用下橫觀各向同性材料環(huán)板的位移和應力場，從二維簡化解構造出相應的材料沿厚度方向任意變化的圓環(huán)板的三維解 。 對于對邊簡支矩形板的彎曲問題，通過在級數(shù)項外引入適當?shù)亩囗検巾?，一方面簡化了求解，另一方面也加快了解的收?[37]。 張曉日等假設材料的力學和電學性質沿板厚方向按統(tǒng)一的指數(shù)函數(shù)形式梯度分布，獲得了周 邊為廣義剛性滑動和廣義簡支兩種邊界條件下軸對稱功能梯度壓電圓板自 由振動問題的 精確頻率方程，數(shù)值求解不同板厚和不同梯度變化情況下的軸對稱圓板自 由振動的固有頻率 。 Chen 等通過假設徑向位移和軸向位移的適當形式，獲得了勻速旋轉的橫觀各向同性功能梯度圓板和圓環(huán)板問題的解析解，該解可考慮材料參數(shù)沿厚度方向的任意分布情形 [14]。 Zhong 等假設材料的力學、電學和熱學性質沿板厚方向按統(tǒng)一的指數(shù)函數(shù)形式梯度畢業(yè)論文 第 9 頁 共 35 頁 分布，獲得了四邊簡支、接地和等溫的正交各向異性功能梯度壓電材料矩形板，在 上下表面作用機械載荷、電載荷和熱載荷情況下的三維靜力精確解。 Zhang 等對于任意的材料梯度分布形式，利用 Haar 小波級數(shù)展開法，得到四邊簡支功能梯度材料矩形板和功能梯度壓電材料矩形板彎曲問題的三維 Haar 小波級數(shù)解。 Huang 等考察了雙參數(shù)地基對功能梯度材料板彎曲特性的影響，發(fā)現(xiàn)地基參數(shù)的變化對結構的應力分布影響較大，尤其是對于較厚的梁板，地基參數(shù)越大，中性面的位置逐漸向加載的表面移動，應力分布也逐漸趨于緩和，即應力水平逐漸減小 [46]。 高立名等采用分層法、 Frobenius 法和同倫分析法研究功能梯度平板表面波的傳播問題，發(fā)現(xiàn)表面波的兩種變化模式在功能梯度材料中的效應是不同的 。 分層法雖然概念簡單，但收斂速度較慢 ， Frobenius 法可以求出精確解，但求解過程繁瑣，而且對數(shù)值計算過程的精度要求比較高 。 同倫法可以求得包括 Frobenius 解的一般解，而且收斂速度和精度可以通過選擇適當?shù)膮?shù)來調整 。 潘永東等利用 Peano 級數(shù)展開法研究功能梯度材料中波的傳播特性，求解了圓桿和圓管表面上功能梯度涂層的表面波頻散曲線，通過數(shù)值計算說明了表面波頻散曲線對涂層材料聲速梯度分布的依賴性，為激光超聲實驗測量表征材料梯度分布特性提供了理論依據(jù) 。 針對功能梯度材料厚板中波的傳播問題，Chen 等提出了回傳射線矩陣分析方法 ， 采用均勻層合模型將各向同性功能梯度材料板劃分成若干層，在每層建立互反的兩個局部坐標系并獲得層內的相位關系， 根據(jù)界面處連續(xù)條件獲得散射關系，最終建立體系的回傳矩陣整體代數(shù)方程 。 由于采用了兩個局部坐標系，該方法摒除了具有正實部指數(shù)的指數(shù)函數(shù)，從而在計算中避免了傳統(tǒng)狀態(tài)空問法中普遍存在的大數(shù)相減，得到了絕對穩(wěn)定的數(shù)值計算結果，表明回傳射線矩陣法對于結構高頻振動計算具有不可替代的優(yōu)越性 。 從上述最一般的理論出發(fā)，文獻對非均勻彈性桿、各向異性功能梯度板 (包括層合板 )和功能梯度壓電薄膜及其體波諧振器內的波傳播問題進行了系統(tǒng)的研究 [13]。 半解析半數(shù)值解 傳統(tǒng)的狀態(tài)空 間 法只能適用于具有特殊邊界條件 (如簡支、滑支 )的 結構的求解 。 為了克服這一限制，范家讓提出用疊加原理和廣義函數(shù)的 Fourier 級數(shù)展開來求解具有固支和 自 由邊界的板殼結構靜動力響應。 Chen 等針對復合材料梁板提出了一種有效的半解析半數(shù)值方法，該方法將狀態(tài)空 間 法和微分求積法結合在一起，在面內采用微分求積法離散可以處理任意側面邊界條件，而在厚度方向采用狀態(tài)空 間 法進行精確求解 。 該方畢業(yè)論文 第 10 頁 共 35 頁 法拓寬了彈性力學傳統(tǒng)狀態(tài)空問法的求解范圍，可以用于分析層合和厚度方向功能梯度梁、板、殼結構靜動力問題 。 Li 等利用該方法獲得了功能梯度層合梁以及壓電層合板在簡支、固支、 自 由等邊界條件下的 自 由振動與靜力彎曲問題的解答 。 在計算中發(fā)現(xiàn)，當微分求積法中點數(shù)取值較大時，會產生一定的數(shù)值不穩(wěn)定問題，為此 Li 等采用 Nagem和 Williams 針對空 間 結構分析提出的 耦 合節(jié)點矩陣法，進行功能梯度梁的熱應力分析和多跨板的振動分析，克服了數(shù)值不穩(wěn)定的現(xiàn)象 [6]。 Nie 等對狀態(tài)空 間 法和微分求積法相結合的半解析半數(shù)值方法進行了改進，提出了將位移和位移的一階導數(shù)作為狀態(tài)變量的求解方法，與應力和位移混合作為狀態(tài)變量的狀態(tài)空問法相比較，更易處理某些邊界條件和圓板中心的正則性條件 。 利用該方法，研究了功能梯度圓板、圓環(huán)板 、扇形板在簡支、固支、 自 由等各種邊界條件下的靜力響應與 自 由振動問題 。 基于狀態(tài)空 間 的微分求積法還被推廣應用于求解雙向功能梯度材料梁和雙向功能梯度材料圓板的分析 。 簡化理論解 曹志遠等人采用經典板殼理論，獲得了各類功能梯度矩形板固有頻率與振型的解析解和各類功能梯度復合材料圓柱殼固有頻率的解析解 。 為了對更復雜條件下功能梯度梁、板、殼結構進行分析，必須發(fā)展適合于功能梯度材料特點的梁、板、殼簡化理論 。 在各種高階板理論中， Soldatos 等提出的含有隨厚度坐標變化的翹曲函數(shù)理論獨具優(yōu)勢 。 Bian 等將 Soldatos 層合板理論推廣用于功能梯度板的分析，構造了狀態(tài)空 間 列式來確定翹曲函數(shù)的厚度方向分布，然后針對 Soldatos 層合板理論構造了狀態(tài)空 間 列式對單跨以及多跨柱形板的彎曲進行求解 。 Li等針對功能梯度材料薄膜，基于經典薄板假設，提出了考慮表面彈性的連續(xù)介質模型 。 研究表明，對于給定的表面材料性質，當薄膜厚度減小到微米尺度時，其抗彎剛度和固有頻率表現(xiàn)出明顯的尺度效應，并隨厚度的進一步減小，表面效應的影響越來越明顯 。 采用 Soldatos 板理論的位移假設，以表面力學基本方程代替經典理論中表面剪應力為零的約束 條件，確定 自 適應形函數(shù)的具體形式，獲得了無限長功能梯度材料薄膜的解析解 。 李堯臣等基于若干基本假設，提出了指數(shù)型功能梯度壓電材料圓板在軸對稱載荷作用下的簡化理論與解析解，獲得了板的周邊固定或簡支并接地情況下中性層法線轉角的解和用 FourierBessel 級數(shù)表示的電勢解，這個解有足夠的精度，在形式上