做nCode DesignLife案例教程這個系列寫到現(xiàn)在已經(jīng)是第十八篇了。經(jīng)常有同行拿著前面文章里的截圖來問我問題我發(fā)現(xiàn)大家卡住的往往不是軟件該點哪里而是不太清楚“這一項為什么要這么設”。今天這篇案例我想認真聊一個我在商用車項目里反復做過很多輪的結構駕駛室懸置支架在隨機振動激勵下的疲勞壽命評估。整個過程覆蓋了DesignLife最核心的幾項能力包括FE結果導入、PSD載荷定義、頻域疲勞計算、焊接結構壽命評估基本就是日常結構耐久分析的主流套路。如果你剛接觸疲勞仿真或者已經(jīng)用DesignLife但只會按教程點不會根據(jù)工況調(diào)參數(shù)這篇應該能幫你少走不少彎路。先說清楚這篇的案例對象。支架是某型商用車駕駛室后懸置上連接車架與空氣彈簧的焊接組合件材料為Q355B主結構板厚6mm焊接方式為MAG焊。這個件在整車耐久試驗里出過裂紋裂紋位置集中在焊縫焊趾處所以需要用DesignLife做一輪壽命排查和結構優(yōu)化對比。隨機振動載荷數(shù)據(jù)來自試驗場實測加速度信號采集點在車架側安裝座上一共三個方向。整個分析鏈路不復雜但里面涉及到的細節(jié)非常多任何一個環(huán)節(jié)參數(shù)填錯算出來的壽命可能差一個數(shù)量級。1. 隨機振動疲勞案例的核心思路1.1 為什么非要選“PSD隨機振動焊接支架”做案例很多初學者一提到疲勞分析第一反應是拿一個正弦掃頻或恒定幅值載荷去算。但真實工程結構尤其是商用車、工程機械上的件受到的激勵幾乎都是寬帶隨機信號。路面不平度、發(fā)動機燃燒波動、輪胎與路面相互作用這些疊加到一起之后時域信號看不出明顯周期。對這類載荷如果強行簡化為等幅正弦要么把壽命算得極其保守要么在某些頻率上嚴重失真。隨機振動疲勞的標準做法是把時域加速度信號轉(zhuǎn)換成功率譜密度PSD在頻域里描述“每個頻率上帶多少能量”然后用頻域疲勞方法估算損傷。Selecting這個案例就是想展示從實測隨機信號到PSD再到DesignLife頻域疲勞計算的全過程。焊接支架又恰好是所有結構里最適合做疲勞案例的對象因為焊接結構存在焊趾應力集中疲勞壽命對局部幾何和殘余應力都極其敏感能充分體現(xiàn)DesignLife在細節(jié)評估上的價值。順帶提一句這個支架在試驗場跑完耐久規(guī)范后焊趾處出現(xiàn)可見裂紋但母材區(qū)域完好。這很典型焊接結構的疲勞壽命往往遠低于母材工程里做優(yōu)化時優(yōu)先查焊縫方向不會錯。1.2 從有限元應力到壽命結果的關鍵鏈路在DesignLife里做隨機振動疲勞底層邏輯并不復雜核心鏈路可以拆成四步用有限元軟件求解結構在單位載荷或模態(tài)空間下的應力響應把應力結果導入DesignLife在DesignLife里讀入PSD載荷譜指定通道與FE結果的對應關系用頻域疲勞模型比如Dirlik法計算各頻點應力循環(huán)的累積損傷換算成壽命。這條鏈路的關鍵在于DesignLife本身并不做結構動力學求解而是把外部FE計算的應力結果和試驗測得的載荷譜結合起來做損傷計算。所以FE結果的精度直接決定后面所有分析的上限很多人在DesignLife里怎么調(diào)都算不準回頭看往往是FE模型輸出就有問題這個后面細說。模態(tài)空間的做法更值得提。對于隨機振動分析一般先在有限元里做模態(tài)分析提取一定階數(shù)的模態(tài)頻率和振型再基于模態(tài)疊加法計算應力響應。這個方法的好處是計算效率高而且天然適配PSD載荷的頻域表達。需要特別注意模態(tài)提取的階數(shù)要覆蓋載荷譜的主要能量頻帶通常建議截止到最高分析頻率的1.2到1.5倍否則高頻響應會被截斷壽命結果會偏樂觀。2. FE模型準備與DesignLife項目搭建2.1 有限元模型里必須先解決的三件事DesignLife的疲勞計算精度一半以上取決于FE結果質(zhì)量。我在這個案例里吃過的虧基本都集中在下面三個地方。第一是網(wǎng)格尺寸。疲勞分析對局部應力梯度很敏感焊趾位置的網(wǎng)格如果太粗熱點應力被平均掉壽命會高得離譜。支架焊縫區(qū)域我一般控制單元尺寸在2mm左右母材區(qū)域可以放到5mm到8mm過渡區(qū)域要平緩。要注意的是DesignLife讀取的是FE節(jié)點應力或單元應力局部網(wǎng)格密度直接決定能否捕捉到真實應力集中。如果網(wǎng)格從8mm直接跳到2mm中間沒有過渡求解出來的應力場會出現(xiàn)不連續(xù)壽命云圖上會出現(xiàn)片狀假熱點。第二是單位體系。這是我在實際項目里遇到最多的問題沒有之一。很多FE模型來自不同工程師有人用mm-N-tonne體系有人用m-N-kg體系還有人用cm。DesignLife不會自動幫你做單位換算載荷通道、材料參數(shù)、應力結果必須心里有一致標準。比如應力單位是MPa還是PaPSD單位是g2/Hz還是(m/s2)2/Hz一個沒對齊壽命結果差10?倍都不奇怪。我在項目的全局設置里會固定好單位建議你也養(yǎng)成這個習慣。第三是模態(tài)數(shù)量。前面說了模態(tài)疊加法需要提取模態(tài)我見過有人圖省事只提了前10階結果載荷譜最大能量在40Hz附近結構第12階模態(tài)正好在那附近激勵起來但沒被算進去壽命結果自然“完美”.正確的做法是先看PSD載荷的主能量頻帶再回頭確認提取的模態(tài)是否覆蓋了這個范圍。一般支架類結構模態(tài)密度不高提取到500Hz基本夠用但前提是載荷譜能量確實衰減到了可忽略的水平。2.2 在DesignLife里搭項目數(shù)據(jù)塊的連接邏輯打開DesignLife這里以nCode 2021及以上版本為例界面邏輯和GlyphWorks一樣是數(shù)據(jù)塊流程化串聯(lián)。新建一個分析頁面左側工具箱里把下面幾個數(shù)據(jù)塊拖到工作區(qū)FE Input負責讀取有限元結果文件如ANSYS的.rst、Nastran的.op2或Abaqus的.odbPSD Input讀取實測加速度PSD譜或時域信號轉(zhuǎn)換后的頻譜Material Database選擇材料SN曲線Vibration Analysis頻域疲勞求解器Display顯示壽命云圖、損傷云圖。連線順序上我的習慣是把FE Input的應力輸出同時連接到Vibration Analysis的FE應力輸入和Material輸入把PSD Input的載荷譜也連接到Vibration Analysis的載荷輸入。這里有一個關鍵點Material Database塊里有材料模型但你要在Vibration Analysis塊內(nèi)部指定“哪個材料賦給哪個零件”如果FE模型里有多個材料分組對應關系要單獨配置。很多教程會讓你把所有塊都連成一條線我的建議是不要。把Material Database單獨連到Vibration Analysis邏輯更清晰排查問題也方便。PSD輸入更是如此不要和FE結果混在一個路徑里后面改載荷譜時容易誤操作。在FE Input設置里要指定結果類型。對于隨機振動分析一般選擇單位載荷下的應力結果stress results這可以是模態(tài)坐標下的單位應力也可以是頻率響應函數(shù)。DesignLife不是直接讀你工況下的絕對應力而是結合PSD載荷去算響應所以FE輸入最好保存成“單位載荷工況”。如果你直接把某個特定工況的應力導進來后面PSD載荷會在這個基礎上放大物理意義就錯了。3. 載荷譜處理與材料參數(shù)設置的細節(jié)3.1 實測加速度信號怎么變成PSD譜通道映射怎么配案例里的實測信號是時域加速度采樣率1024Hz時長120秒三個方向X、Y、Z同步采集。處理時域信號的常規(guī)做法是分段加窗做FFT變換再對多段結果取平均得到PSD。這個轉(zhuǎn)換可以在nCode GlyphWorks里用FFT統(tǒng)計塊完成也可以直接用外部軟件處理好了再導入DesignLife。我這邊更推薦在GlyphWorks里做原因很簡單實測信號里通常有毛刺和漂移尤其是加速度信號低頻漂移會污染PSD的低頻段。你可以在GlyphWorks里先做高通濾波截止頻率設到0.5Hz左右把趨勢項剔除然后再做PSD轉(zhuǎn)換。DesignLife里讀入PSD時支持文本文件、Excel文件以及nCode自己的Data File格式。這里有一個極易踩坑的地方。PSD文件里通道標記必須和FE應力結果的方向定義對應好不是簡單三行數(shù)據(jù)就行。比如PSD定義了Z向加速度譜DesignLife在計算時會把Z向單位載荷應力結果和這個PSD相乘再疊加其他方向。如果PSD里的通道名和FE結果里的方向名對不上軟件會報警告或者干脆不計算。我遇到過一次情況FE結果里通道名是“ACC_Z”PSD里寫成“Z_ACC”軟件直接報錯。所以建議在FE導出和PSD制作時就統(tǒng)一好命名規(guī)則比如全部用“X_Direction”“Y_Direction”“Z_Direction”省得后面來回對。關于PSD的單位實測加速度通常以g9.81m/s2為單位PSD單位是g2/Hz。如果FE模型里用的是mm/s2材料單位也是N/mm2MPa那么PSD的縱軸必須把g2/Hz換算成(mm/s2)2/Hz也就是乘以(9810)2否則應力響應幅值會整體偏差一個系數(shù)。單位換算這步我在項目里會寫成備忘錄每次建模前先確認比到最后再返工強。3.2 材料SN曲線母材隨便選焊縫不能隨便選材料參數(shù)是疲勞分析里最容易被低估的部分。DesignLife自帶材料庫里面有很多常用材料的標準SN曲線比如鋼、鋁、鑄鐵等。對于母材區(qū)域如果手里沒有實測材料數(shù)據(jù)用庫里的Q355B或S355曲線做對比分析是完全可行的。但焊接區(qū)域就不一樣了。焊接接頭的疲勞強度遠低于母材不只是因為材料變了更關鍵的是焊趾處存在幾何不連續(xù)、焊接殘余應力和微缺陷。國際焊接學會IIW的規(guī)范里根據(jù)接頭形式和載荷方向定義了不同的FAT等級疲勞等級即2×10?循環(huán)對應的應力幅值。比如非承橫向?qū)雍缚pFAT等級可能是80MPa而十字角焊縫可能只有63MPa甚至更低。我在這案例里對支架焊縫選擇的是IIW規(guī)范里的FAT80曲線同時開啟了對焊接接頭的厚度修正。為什么選80而不是90或100因為MAG焊角焊縫、焊趾位置受拉壓循環(huán)按照IIW的判定表FAT80是合理的起步值。如果你有試驗數(shù)據(jù)對標可以再用試驗結果校準。還有一點值得說焊接結構的SN曲線在高周疲勞區(qū)斜率不是標準鋼的m3或者m5焊接接頭通常取m3。DesignLife里材料曲線的斜率設置會影響壽命外推如果載荷譜里大部分循環(huán)應力幅都在疲勞極限以下曲線長壽命區(qū)的斜率差異會明顯改變結果。不要圖省事直接用母材曲線代替焊縫曲線那基本等于白算。如果案例里存在點焊焊點DesignLife里還有專門的Spot Weld Analysis數(shù)據(jù)塊需要基于焊點力和節(jié)點力計算壽命和這里講的焊縫評估流程完全不同。這篇案例里支架全部為連續(xù)角焊縫所以用Vibration Analysis加材料庫曲線就夠了點焊那套以后有空單獨寫。3.3 平均應力修正要不要開怎么開頻域疲勞計算里平均應力修正是一個繞不開的選項。隨機振動載荷下的應力循環(huán)均值往往不為零尤其是自重、裝配應力這些恒定載荷疊加在動態(tài)應力上。如果完全不修正會低估損傷修正過度又會把壽命算得過于保守。DesignLife里提供了多種平均應力修正方法Goodman、Gerber、FKM、Soderberg等。我的選擇經(jīng)驗是對于母材用Goodman修正偏保守且工程普遍接受對于焊接結構反而要謹慎因為焊接殘余應力很大焊趾處的實際平均應力受殘余應力主導而不是外在載荷。國際規(guī)范里對焊接接頭很多時候直接不考慮平均應力修正或者用特殊的焊縫修正曲線。這個案例里我采用了較為穩(wěn)妥的方案母材區(qū)域開Goodman修正焊縫區(qū)域關閉平均應力修正。這樣處理的原因很簡單焊接殘余應力通常是屈服強度級別外部載荷帶來的平均應力相對是小量再去用Goodman放大或縮小反而讓結果偏離真實。如果你不確定自己結構里的殘余應力狀態(tài)建議先在焊縫處做一組有修正和無修正的對比觀察壽命差異再決定。4. 求解器參數(shù)設置與后處理結果解讀4.1 頻域疲勞求解器里那些參數(shù)到底怎么填Vibration Analysis數(shù)據(jù)塊的屬性對話框打開之后其實選項不算多但每一項都值得細看。第一個是疲勞計算方法。DesignLife頻域疲勞提供Dirlik、Steinberg、Wirsching等幾種經(jīng)驗公式。Steinberg也叫三區(qū)間法工程上用得很多它基于高斯分布假設把應力幅值分三個區(qū)間估算損傷非?？斓珜Ψ歉咚馆d荷會失真。Dirlik方法通過四個譜矩參數(shù)構造幅值概率密度函數(shù)對寬帶高斯隨機過程擬合效果更好是目前工程界的默認選擇。本案例我選的是Dirlik。第二個是損傷累積規(guī)則。DesignLife默認用Miner線性累積法則把各應力幅值造成的損傷線性疊加損傷總和到1就認為發(fā)生疲勞破壞。雖然Miner法則不考慮載荷順序效應但工程上仍然是最主流、最便于工程判斷的方法。第三個是“置信度”相關參數(shù)。有些模板會提供“Probability of Survival”設置比如99.7%或50%這會影響默認的SN曲線位置設計曲線還是均值曲線。如果你手里只有均值SN曲線要得到偏保守的壽命可以在結果分析時適當用安全系數(shù)或置信度參數(shù)往下修。這里我習慣不額外加置信度系數(shù)因為在對比方案階段我更關心相對壽命差異絕對壽命可以留給試驗對標后再統(tǒng)一修正。第四個是材料分區(qū)指定。在Analysis數(shù)據(jù)塊的“Assignment”頁簽里你可以為不同材料分組指定對應的SN曲線和修正方式。這個案例里材料列表包含母材和焊縫兩組分別在模型里定義成不同的單元屬性DesignLife會根據(jù)材料名自動匹配前提是你把材料數(shù)據(jù)庫里的材料命名和FE模型里的材料名對應好。我習慣把焊縫材料單獨命名為“Weld_Q355B_FAT80”避免混淆。求解時還有一個容易被忽略的選項頻率范圍與步長的設定。DesignLife讀入PSD后會按照PSD定義的頻率點計算應力響應。如果FE結果是基于模態(tài)疊加的頻響函數(shù)需要確保模態(tài)響應文件覆蓋了PSD的最高頻率否則高過文件范圍的頻譜會被截斷。這里沒有太多花活就是檢查頻率范圍是否匹配。4.2 壽命云圖出來之后怎么判斷結果合不合理計算完成后Display塊會顯示壽命云圖。正常焊接支架結構在隨機振動下壽命最小值一般出現(xiàn)在焊趾應力集中區(qū)且壽命數(shù)量級通常在10?到10?循環(huán)之間如果算出來是103或者10?就要回頭看設置。優(yōu)先看三個地方一是壽命最小的位置是不是和試驗裂紋位置吻合。通行的做法是把壽命云圖和真實裂紋照片放在一起對比。本案例里裂紋出現(xiàn)在角焊縫焊趾端部壽命云圖上同樣位置也確實出現(xiàn)了最小壽命點這說明輸入條件基本可信。二是損傷分布是否集中在窄帶頻率附近。隨機振動疲勞的損傷往往由少數(shù)幾個主導模態(tài)貢獻你可以看某個頻率下的響應PSD確認損傷峰值頻率和模態(tài)頻率對得上。如果全頻段都均勻損傷通常說明載荷譜處理有問題。三是最小壽命值是否與試驗耐久時間處于同一量級。商用車駕駛室懸置支架的耐久規(guī)范一般是幾十小時強化試驗換成偽損傷或等效壽命落在合理區(qū)間即可。如果差了一個數(shù)量級以上優(yōu)先懷疑材料曲線FAT等級選錯或者單位換算不對。如果壽命不達標后處理上還可以做“Contribution”分析查看哪個方向載荷對損傷貢獻最大哪個頻率段貢獻最大為結構優(yōu)化提供方向。支架類結構常見優(yōu)化手段包括增加板厚、增加翻邊、改變焊縫位置避開高應力區(qū)。這些優(yōu)化方案可以快速修改FE模型再導入DesignLife對比相對壽命不需要重新做載荷處理。5. 實際操作中最容易翻車的幾個問題5.1 常見報錯與排查速查表寫代碼有bug用DesignLife也有它的“報錯藝術”。下面這張表是我這幾年用nCode處理疲勞問題時整理的高頻問題排查清單每一條都是真實踩過的?，F(xiàn)象可能原因解決辦法計算結果壽命極短如小于1次循環(huán)單位不一致應力被放大1000倍檢查FE模型、材料曲線、PSD載荷單位是否統(tǒng)一計算結果壽命極長如101?SN曲線斜率太大或應力都低于疲勞極限檢查載荷是否施加成功PSD通道是否映射正確報錯“No valid data on input”FE結果里沒有應力數(shù)據(jù)或結果類型錯誤確認FE結果文件是否為應力結果不是位移結果報錯“Channel mismatch”PSD通道名與FE結果通道名不一致統(tǒng)一通道命名按方向?qū)P系重新映射壽命云圖全為均勻色塊網(wǎng)格太粗導致應力梯度丟失或材料全賦成一種加密焊縫局部網(wǎng)格核對材料分組分配某個頻點損傷異常高模態(tài)截斷或PSD頻率點與FE頻率響應不匹配增加模態(tài)階數(shù)保證覆蓋到載荷譜最高頻率的1.2倍焊縫區(qū)域壽命和母材一樣沒有給焊縫指定專用SN曲線用IIW FAT等級曲線單獨賦給焊縫組計算能跑但結果為N/A材料庫曲線缺失或單位格式不合法在材料庫里重新創(chuàng)建SN曲線檢查是否輸入為幅值而非范圍這張表不是萬能的但能覆蓋日常80%的卡殼場景。真遇到解決不了的問題還是要回到數(shù)據(jù)本身用Display的XY繪圖功能把中間過程調(diào)出來看比如看PSD的波形是否正常、應力響應是否奇異等比瞎試參數(shù)有效率得多。5.2 一次焊縫壽命對不上的排查實錄這個案例我曾經(jīng)卡了整整兩天現(xiàn)象是試驗顯示焊趾開裂但DesignLife算出來的壽命云圖熱區(qū)在母材遠離焊縫的位置焊縫處壽命反而很好。最初以為是材料曲線沒賦對反復檢查后發(fā)現(xiàn)焊縫材料確實已指定為FAT80問題不在材料。后來我把焊縫處的應力響應單獨導出來看發(fā)現(xiàn)焊縫焊趾處單元應力值明顯低于周圍母材這不對勁因為焊趾是應力集中區(qū)理論上應力應該更高。再回到FE模型查發(fā)現(xiàn)焊縫是用實體單元模擬的但焊縫與母材之間的接觸定義成了共節(jié)點焊縫倒角位置網(wǎng)格太粗局部應力峰被“磨平”了。重新加密焊趾倒角處網(wǎng)格并把圓角半徑按實際焊腳尺寸修正后壽命熱區(qū)終于回到了焊趾位置數(shù)值也和試驗對上了。這件事給我留下的教訓是DesignLife再強大也只是個后處理計算器決定結果上限的依然是FE模型的應力場質(zhì)量。做焊接疲勞焊趾附近的網(wǎng)格必須細到能分辨幾何應力集中一般建議單元尺寸不超過焊腳高度的一半。如果只是粗網(wǎng)格大模型算出來壽命再漂亮都是自欺欺人。5.3 材料庫曲線缺失時的替代方案工程里經(jīng)常遇到材料庫里沒有的材料比如某些特殊牌號鋼或鋁合金型材。這時候不要直接拿庫里最接近的隨便頂上至少做兩步第一步看材料是軋制態(tài)還是鑄造態(tài)鍛造鋁合金和鑄造鋁合金的疲勞性能差很多第二步如果查標準鋁合金母材可以參考FKM指南或Eurocode 9的細節(jié)等級。對于焊接鋁結構IIW也有對應的FAT等級例如非承載橫向?qū)雍缚pFAT為40MPa左右遠低于鋼的80。把外部SN曲線錄入DesignLife的好處是后續(xù)所有項目都能復用。我習慣把常用的焊接FAT等級曲線統(tǒng)一定義在一個自定義材料庫里命名規(guī)范是“FAT80_Steel_Weld”“FAT40_Al_Weld”這樣導入時一眼就能看懂。材料庫維護看著不起眼但對效率提升非常明顯強烈建議養(yǎng)成習慣。6. 對設計優(yōu)化的反饋和一實戰(zhàn)建議6.1 壽命結果怎么變成可落地的改進建議疲勞分析最終要服務的是設計改進不是出一張漂亮云圖交差。當壽命熱區(qū)鎖定在焊趾后我一般會結合DesignLife里的損傷頻率貢獻和應力分布做三個方向的優(yōu)化建議。第一局部幾何優(yōu)化。如果損傷集中處是焊縫端部可以建議把焊縫端部延出讓應力流更平緩如果是焊根位置考慮增加熔深或改為全熔透焊縫。這個案例里我們最后通過增加焊縫退刀槽、增大焊趾圓弧半徑把局部熱點應力降了大概18%計算壽命提升了約3倍。第二載荷路徑優(yōu)化。支架的裂紋很多時候不是強度不夠而是結構剛度分配不合理導致某個撐腳承受了不成比例的載荷。通過壽命云圖可以看到哪個方向貢獻最大這時可以考慮調(diào)整安裝點位置、增加減震墊剛度等。第三制造工藝反饋。焊接順序、焊接變形、焊腳尺寸偏差都會影響實際疲勞強度。仿真算出剩余壽命不高時工藝上一定要強調(diào)焊腳成型質(zhì)量必要時增加焊后打磨或TIG重熔。6.2 新手上手DesignLife的三條個人建議最后分享一點個人體會。帶過不少新人做nCode我發(fā)現(xiàn)上手最快的那批人并不是軟件操作最熟的而是對“載荷怎么來、應力怎么算、損傷怎么累”這條物理鏈理解最透徹的。軟件本身學起來就那么幾個數(shù)據(jù)塊難點全在于當你看到壽命異常時能不能快速定位是哪一環(huán)出了問題。第一建議是建議每個案例都留一份“設置記錄”把模型單位、PSD處理方式、材料曲線來源、修正方法全部寫清楚。這樣不僅方便自己復查后續(xù)做試驗對標和匯報都能直接用。第二建議是多花時間在FE模型上而不是陷在DesignLife里調(diào)參。疲勞計算對網(wǎng)格、約束、載荷施加方式的敏感度超過大多數(shù)人的直覺。與其反復試材料曲線不如先回頭把焊趾局部網(wǎng)格細化往往一點就通。第三建議是有條件就做試驗對標哪怕只做一個等幅疲勞試樣的數(shù)據(jù)都行。DesignLife的默認材料曲線只能保證趨勢對絕對壽命一定有偏差。通過一次簡單的三點彎曲或拉壓疲勞試驗校準一下材料曲線的FAT等級或斜率后面整個項目的置信度都會不一樣。這個案例系列寫到這里核心的流程和坑點都覆蓋得差不多了。如果照著步驟做一遍相信你對nCode DesignLife的隨機振動疲勞分析會有比較完整的把握。工具只是放大器真正值錢的是你對工程問題的理解。希望這篇能幫你少走一些彎路后面有機會再把點焊疲勞、熱機疲勞這些專題慢慢補上。