【正文】 50E6MPa， xyG = xzG =50E6MPa， yzG =20E6MPa， xyPR 畢業(yè)論文 25 = yzPR =， xzPR =。從圖 56 可以看出每個點的累積概率，說明可以看出當 FAILURE 值小于零的概率，也就是失效概率。損傷檢測和修復可以使得飛機結(jié)構(gòu)在修復過后，結(jié)構(gòu)強度恢復到一定的水平。 （ 2） 繼續(xù)用 ANSYS 模擬整個損傷容限的過程， 對 所有的過程的變量進行分析和計算。 感謝他們在 我畢業(yè)設計期間 在學習上對我的幫助！ 最后，我還要特別感謝我的父母，他們對我無私的愛和關心，對我的養(yǎng)育之恩、不懈的鼓勵，使我能夠以順利地完成學業(yè)！ 感謝所有幫助我，關心我的老師和同學！ 還有感謝一些有關有限元、 MATLAB 和蒙特卡洛的論壇網(wǎng)站。 （或者是 %壓縮下剩余強度 Lth=[0,15,20,30,40,50,60,100,125,150]。 %服役損傷時刻 (均勻分布 ) %t=[0:3000]； %損傷尺寸檢測概率 Livh=[0,10,20,25,30,40,50,60,70,80,90,100,110,120]。 %修復 %Krhf=Krkf=。 %溫度和強度修正系數(shù) Tk=[20,60,100,125,140,160]。 Rm222=0。 Rc=0。 Iivk=0。 畢業(yè)論文 34 end else Rm11=0。 %產(chǎn)生孔和裂紋的損傷次數(shù) if (N3=1) Lsh33=25+exprnd(20,N3,1)。 %均勻分布產(chǎn)生損傷時刻 end N=N3+N4。amp。 end Rc(i)=Rsk44(m11)。 elseif(L120) Pivh(m11)=1。 end Iiih(m11)=geornd(Piih(m11),1)。L0) Pivk(m11)=spline(Livk,Plivk,L)。 end 畢業(yè)論文 37 Iiik(m11)=geornd(Piik(m11),1)。 if(L=40amp。 m11=m11+1。 if(L=70amp。 %修復完成時間 %產(chǎn)生每個尺寸的修復次數(shù)（時間） if(N=1) for m11=1:N 畢業(yè)論文 36 if(N3=1) for m22=1:N3 if(t(m11)==t33(m22)) L=Lsh33(m22)。 else t=0。 elseif (N3=1amp。 %產(chǎn)生孔和裂紋損傷尺寸 for m11=1:N4 Rsk44(m11)=spline(Ltk,Rltk,0)spline(Ltk,Rltk,Lsk44(m11))。 end else Rm22=0。 N2=poissrnd(,1)。 Piih=0。 Rsk44=0。 for J=1:100000 %初始化各個變量 Lm11=0。 %亞音速飛行應力超越數(shù)曲線 %S11=[0,37,55,70,85,102,120,137,154]。 Livk=[0,10,20,25,30,40,50,60]。） %服役損傷 %Lsh=[0,25,30,35,40,60,100]。 %拉伸下剩余強度 Lth=[0,15,20,30,40,50,60,100,125,150]。 [3] Sheldonross主編，趙選民等譯，概率論基礎教程 [M]， 北京： 機械工業(yè)出版社 ，南加州大學， ，85~159. [4] M. W. Long and J. D. Narciso， PROBABILISTIC DESIGN METHODOLOGY FOR COMPOSITE AIRCRAFT STRUCTURES[R]， FAA Technical report， June 1999. [5] Ushakov, A. Stewart, I. Mishulin, and A. Pankov， PROBABILISTIC DESIGN OF DAMAGE TOLERANT COMPOSITE AIRCRAFT STRUCTURES[R]， FAA Technical report ,January 2020. [6] 沈觀林 , 胡更開 ,復合材料力學 [M], 北京 ： 清華大學出版社 , 2020. [7] 飛思科技產(chǎn)品研發(fā)中心 編著， MATLAB7 基礎與提高 [M]， 北京： 電子工業(yè)出版社， . [8] 王文健主編，試驗數(shù)據(jù)分析處理與軟件應用 [M]， 北京： 電子工業(yè)出版社， . [9] 王明強 ,朱彤 , 應用 ANSYS 概率有限元法的連桿結(jié)構(gòu)強度可靠性分析 [J]. 現(xiàn)代制造工程 2020 年第 3期 .54~57. [10] 趙國華 , 用 ANSYS 進行復合材料層合板的沖擊應力場模擬 [J]. 《 機床與液壓》 ~161. [11] 高娟 , 羅奇峰 , 車偉 .蒙特卡羅法理論及其在 ANSYS 中的實現(xiàn) [J], 青島理工大學學報第 29 卷第 4 期 ~22. [12] 章慧健 , 工程結(jié)構(gòu)可靠度計算的 MATLAB 實現(xiàn) [J], 計算機應用 , 第 27 卷 1 期 .154~156. [13] 何鳳霞 ,張翠蓮 . 蒙特卡羅方法的應用及算例 [J]. 華北電力大學學報 , 第 32 卷第 3 期 ,110~113. [14] 馮曉波 ,楊樺 , 用 MATLAB 實現(xiàn)蒙特卡羅法計算結(jié)構(gòu)可靠度 [J], 中國農(nóng)村水利水電 ,2020 年第 8期 .50~51. 畢業(yè)論文 31 致 謝 在 大 三下時候 ，我很幸運 找到 一位好導師 —— 教授。 （ 2） 由于很難確定小失效概率和數(shù)據(jù)精確性的置信水平，僅僅在有大量數(shù)據(jù)的情況下使用結(jié)構(gòu)概率分析方法可以使結(jié)果更加精確； （ 3） 對結(jié)構(gòu)失效概率的估算沒有考慮到載荷在飛機結(jié)構(gòu)上的傳遞和再分配的問題，這樣的假設比較保守，會使結(jié)構(gòu)的重量有所增加。從而可以用數(shù)學表達式來建立變量之間的模型。選擇 MCM 直接法抽樣方法對本例題進 行模擬分析，取抽樣次數(shù) 100。所以簡化了這個模型，部分不是很合理。記憶意味著避免了重復抽取已經(jīng)計算過的樣本。 圖 51 可靠性分析流程 概率設計分析模塊 PDS 簡介 ANSYS 中可靠性分析數(shù)據(jù)流向如圖 52 所示。我們假設從飛機結(jié)構(gòu)上取一塊復合材料層合板，將損傷和缺陷等效為板中帶圓孔，這樣就建立一個帶孔的復合材料層合板模型。 步驟 4 服役期間，通過一二三步驟之后，我們可以計算出每個飛行小時的 剩余強度。 圖 41 隨機數(shù)的產(chǎn)生 這里要經(jīng)過三個步驟：首先是均勻分布隨機數(shù)的產(chǎn)生，產(chǎn)生所謂的偽隨機數(shù)，常用的方法有迭代取中法、移位指令加法和同余法等；其次對產(chǎn)生的隨機數(shù)進行檢驗，看其是否滿足我們希望的統(tǒng)計性質(zhì) (隨機性檢驗主要有四個：均勻性檢驗、獨立性檢驗 (或不相關性檢驗 )、組合規(guī)律性檢驗、無連貫性檢驗 )；最后進行隨機變量的抽樣，目的是為了得到具有一定分布的隨機數(shù)。所以通常情況下對連續(xù)型分布采用直接抽樣是有一定困難的。 模擬法這里主要指蒙特卡羅法。對于拉伸載荷，是比較保守的?？毯蹟?shù)據(jù)來自參考文獻 [5]，如表 315 所示。 樣條曲線擬合 損傷表面的強度退化函數(shù) 考慮到?jīng)_擊損傷的兩種類型：一是刻痕損傷類型；二是孔 +裂紋的損傷類型。 表 35服役損傷刻痕的超越數(shù) 尺寸 (mm) 每平方每小時 超越數(shù) 10^5 每 平方米每個壽命 超越數(shù) 0 10 25 10 30 35 6 40 60 100 畢業(yè)論文 11 表 36服役損傷刻痕的累積概率 刻痕尺寸 (mm) 每個壽命每 平方米 累積概率 0 0 25 0 30 35 40 60 100 圖 33 刻痕損傷的累積概率分布 表 37 服役損傷孔 +裂紋的超越數(shù) 尺寸 (mm) 每平方米每小時 超越數(shù) 10^5 每 平方米每壽命 超越數(shù) 0 5 10 5 25 2 30 35 1 40 60 100 畢業(yè)論文 12 表 38服役損傷孔 +裂紋的累積概率 孔 +裂紋 mm 每個壽命每 平方米 累積概率 0 0 10 0 25 30 35 40 60 100 圖 34服役損傷孔 +裂紋的累積概率分布 整個 過程中，損傷尺寸發(fā)生次數(shù)，我們可以把在每個壽命每 平方米中最大超越數(shù)看作超出期望值，那樣的話，我們就可以通過泊松分布產(chǎn)生每個壽命每 平方米的發(fā)生次數(shù)。 按前面的處理方法，結(jié)果如表 34 所示。 概率分布擬合 載荷應力超越數(shù)以及分析擬合 在亞音速設計載荷發(fā)生的情況下，我們判斷載荷大小的發(fā)生概率。對于每個時間段內(nèi)的失效概率的計算，都是通過計算每次飛行中的失效概率來確定的。每次的損傷，檢測和修復的時間都是隨機產(chǎn)生的，這里假設檢測和修復的時間是相同的。 模擬過程 整個程序的模擬過程如圖 22 所示： 圖 22 模擬過程 畢業(yè)論文 6 a) 對于每種設計載荷情況，每次迭代都開始于初始的強度值（ 0 0t ? ）?？梢詮闹械玫揭淮物w行的失效概率。失效概率通過確定應力大于強度的概率計算出來。經(jīng)驗表明在復合材料飛機結(jié)構(gòu)產(chǎn)生災難性破壞的最大的原因是在制造過程和服役期間的力學沖擊損傷中產(chǎn)生的缺陷，這些缺陷目視不易發(fā)現(xiàn)。得到各個變量的概率分布圖和其他變量的樣條曲線圖。因此就需要應用使用中的損傷容限準則預測復合材料結(jié)構(gòu)在損傷發(fā)生后 的情況。各種方法有各自的優(yōu)缺點。設計變量還包括了結(jié)構(gòu)尺寸、飛行的載荷情況以及飛行的環(huán)境。以 Mig29 垂直尾翼為例進行詳細的分析和編程計算。接著在通過使用商業(yè)有限元軟件 ANSYS， 對一個帶孔的復合材料層合板，使用參數(shù)化設計語言（ APDL）和概率設計模塊（ PDS），在隨機輸入變量的情況下，計算此板的失效概率。如果僅由傳統(tǒng)的安全系數(shù)來考慮結(jié)構(gòu)設計是很危險的，安全系數(shù)可能太大，或某些情況下卻太小，無法確定結(jié)構(gòu)的可靠性，最嚴重極端的情況就是整個方法可能導致混合及低效的設計。我們要合理將它們應用在工程上。 本文的研究工作 本文是建立在國內(nèi)外的參考文獻的數(shù)據(jù)基礎之上，對四大概率設計方法的之一的概率損傷容限設計方法（ TsAGI 方法） 進行研究。 第四章主要用 MATLAB 語言對整個損傷過程進行蒙特卡洛模擬。因此就需要應用使用中的損傷容限準則預測復合材料結(jié)構(gòu)在損傷發(fā)生后的情況。 在每一個時間間隔，將應力和強度相比較，這種模擬依賴于初始強度狀態(tài)和操作過程中的隨機變量。損傷和未損傷結(jié)構(gòu)的失效概率的計算方法不同。結(jié)構(gòu)的強度由很多設計變量確定，最后的強度值是一個概率統(tǒng)計的變量。事實上，不同類型的損傷的檢測間隔都是根據(jù)結(jié)構(gòu)的構(gòu)型，查找統(tǒng)一的檢測間隔時間表得到的。 畢業(yè)論文 7 第三章 數(shù)據(jù)處理 引言 復合材料結(jié)構(gòu)損傷損傷容限設計的變量有很多，這些變量主要是隨機變量。載荷大小從 037 的載荷應力是不發(fā)生的，因為載荷應力太小了。對這些樣本數(shù)據(jù)進行擬合分布曲線，通過描點法，從曲線形狀可以看出是正態(tài)分布， MATLAB 語言擬合后的結(jié)果如圖 32 所示。 （ 2）制造缺陷 根據(jù)制造中垂直尾翼的檢測數(shù)據(jù)，從而獲得兩類損傷超越數(shù)：一是有關分層或刻痕的；二是有關孔 +裂紋的。在不同損傷類型的情況下，強度退化的函數(shù)是不一樣的?？毯蹣訔l曲線擬合后，如圖 39 所示。強度修復系數(shù)和溫度關系的參數(shù)來自參考文 獻 [5]，經(jīng)過樣條曲線擬合后，關系如圖 311 所示： 表 318 強度修復系數(shù)和溫度關系 溫度 修正系數(shù) 20 60 100 125 140 160 20 40 60 80 100 120 140 160 18000 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 91溫度 T強度修正系數(shù) k k v s . T溫度 修正系數(shù) 圖 311 溫度和強度修正系數(shù)的關系 畢業(yè)論文 19 第四章 用 MATLAB編程實現(xiàn)計算失效概率 引言 MATLAB是一種功能極其強大的科學和工程計算數(shù)學軟件系統(tǒng) ,匯集了大量數(shù)學、統(tǒng)計、科學和工程所需的函 數(shù)。蒙特卡羅法在目前可靠度計算中 ， 被認為是一種相對精確的方法 ， 它不受極限狀態(tài)方程非線性的限制。利用 MATLAB的強大數(shù)值計算功能 ,實現(xiàn)了在