fe-safe基礎(chǔ)培訓(xùn)第三天

fe-safe軟件基礎(chǔ)培訓(xùn)資料03第三天培訓(xùn)內(nèi)容彈性有限元結(jié)果的疲勞計(jì)算過(guò)程單一載荷作用FEA結(jié)果的疲勞分析原理a、讀入彈性計(jì)算的FEA節(jié)點(diǎn)的應(yīng)力張量（6個(gè)）；b、用載荷歷程數(shù)據(jù)乘以應(yīng)力張量，得到每個(gè)張量的歷程數(shù)據(jù)；c、計(jì)算表面內(nèi)主應(yīng)力的時(shí)間歷程（對(duì)接觸則考慮面外應(yīng)力）；主應(yīng)力方向不隨歷程改變。d、從應(yīng)力計(jì)算三個(gè)主應(yīng)變的時(shí)間歷程數(shù)據(jù)；e、對(duì)應(yīng)變壽命疲勞分析(如Brown-Miller分析)，用多軸循環(huán)塑性模型將彈性應(yīng)力-應(yīng)變時(shí)間歷程數(shù)據(jù)轉(zhuǎn)化為彈塑性應(yīng)力-應(yīng)變歷程數(shù)據(jù)（對(duì)S-N曲線分析則忽略此步）；f、對(duì)剪應(yīng)變、Brown-Miller分析，在三個(gè)可能的平面上計(jì)算剪應(yīng)變、法向應(yīng)變、法向應(yīng)力的時(shí)間歷程數(shù)據(jù)；對(duì)S-N曲線分析，則定義一個(gè)與表面垂直的平面（主平面之一），計(jì)算與該平面垂直的應(yīng)力的時(shí)間歷程數(shù)據(jù)；g、計(jì)算每個(gè)平面的疲勞損傷；各平面的每種疲勞循環(huán)用雨流計(jì)數(shù)法統(tǒng)計(jì)，總損傷由每種循環(huán)的損傷疊加而成；具有最短壽命的平面即為裂紋發(fā)生平面，其壽命寫入結(jié)果文件。h、修正耐勞極限。如果一個(gè)平面內(nèi)的各種循環(huán)均低于耐勞極限（所對(duì)應(yīng)的循環(huán)），則沒(méi)有損傷產(chǎn)生；如任一循環(huán)有損傷產(chǎn)生，則自動(dòng)將耐勞極限降低25%，損傷曲線基于修正后的耐勞極限進(jìn)行估算。i、對(duì)每個(gè)單元結(jié)點(diǎn)都重復(fù)上述過(guò)程。其它常用功能單元組group缺省地與模型單元的材料編號(hào)一一對(duì)應(yīng)；強(qiáng)度因子（FOS）：定義：假設(shè)規(guī)定設(shè)計(jì)壽命N對(duì)應(yīng)的應(yīng)力值為[σ]，計(jì)算應(yīng)力為σ，則與規(guī)定設(shè)計(jì)壽命對(duì)應(yīng)的強(qiáng)度安全因子為FOC=[σ]/σ；失效率（FailureRate）在一系列規(guī)定目標(biāo)壽命下，考慮材料/載荷波動(dòng)按統(tǒng)計(jì)規(guī)律分布時(shí)的失效率；疲勞安全系數(shù)（FRF）對(duì)無(wú)限壽命的應(yīng)力-壽命算法，按照Goodman圖計(jì)算疲勞安全系數(shù)，可基于平均應(yīng)力為定值時(shí)的應(yīng)力幅、或應(yīng)力幅為定值時(shí)的平均應(yīng)力、或同時(shí)考慮平均應(yīng)力/應(yīng)力幅的影響。Haigh圖：顯示每個(gè)節(jié)點(diǎn)上最危險(xiǎn)的循環(huán)；材料圖顯；近似材料數(shù)據(jù)；載荷敏感性分析：當(dāng)構(gòu)件上作用有多個(gè)方向的載荷時(shí)，可對(duì)關(guān)心位置的疲勞壽命在依次去掉某方向載荷的情況下進(jìn)行重新計(jì)算。用于指導(dǎo)部件實(shí)驗(yàn)、查驗(yàn)重要載荷方向、疲勞位置等。輸出細(xì)節(jié)：應(yīng)力張量、主應(yīng)力、主應(yīng)變的時(shí)間歷程數(shù)據(jù)，損傷參數(shù)（法向應(yīng)力/應(yīng)變、剪應(yīng)變）；具有最大損傷的結(jié)點(diǎn)列表；無(wú)損傷的結(jié)點(diǎn)順序列表；“通過(guò)”、“失效”、“邊沿”的疲勞結(jié)果顯示。標(biāo)準(zhǔn)分析設(shè)置的保存；批處理分析；比例因子等；載荷歷程定義支持的文件工業(yè)標(biāo)準(zhǔn)文件*.dac；解析多通道文件*.amc；ASCII單/多通道文件*.txt、*.asc；ASCII柱狀圖文件*.txt、*.asc;第三方軟件的文件格式：*.sbf、*.sbr、*.rsb、*.tab用戶在操作界面下自行定義實(shí)例圖示為在其中心位置開有一圓孔的薄板，薄板一端固定，另一端承受在1個(gè)工作循環(huán)內(nèi)承受P=50Mpa、-50Mpa、50Mpa為一個(gè)工作循環(huán)交替變化的動(dòng)壓力，試用ANSYS+fe-safe軟件計(jì)算：在10000個(gè)循環(huán)這一規(guī)定壽命下，該薄板的疲勞安全系數(shù)？已知參數(shù)為：薄板長(zhǎng)L=1m、高H=0.5m、厚T=0.005m、圓孔半徑R=0.13m、拉應(yīng)力P=50MPa；薄板為鋼材，其材料牌號(hào)為SAE-950C-Manten，彈性模量E=2.03×1011Pa、泊松比ν=0.3、密度ρ=7850kg/m3。分析步驟1：先在ANSYS環(huán)境，施加50Mpa載荷作一個(gè)靜力分析；見(jiàn)文件fatigue01.txt。分析步驟2：在fe-safe環(huán)境，打開有限元模型（file.rst）、定義載荷歷程比例數(shù)據(jù)為：1、-1、1，指定材料（SAE-950C-Manten，并將ConstAmpEnduranceLimit(2nf)值設(shè)為2e4）、算法（Brown-Miller/Infinitelife/Goodman）等后進(jìn)行相關(guān)計(jì)算。計(jì)算結(jié)果：最大主應(yīng)力為222Mpa；最小安全系數(shù)為3.523（規(guī)定壽命為1000個(gè)循環(huán)時(shí)），見(jiàn)下圖。圖1薄板受力示意圖薄板主應(yīng)力計(jì)算結(jié)果薄板疲勞安全系數(shù)計(jì)算結(jié)果（注意，不是FOS值）注意：后處理中sx、sy、sz、sxy分別為徑向、垂向、橫向、最小安全系數(shù)。多種載荷作用下FEA結(jié)果的疲勞載分析最多可支持4096種同時(shí)作用的載荷歷程。每種載荷有其自己的載荷歷程，每種載荷各有一個(gè)線彈性FEA解。稱為“放樣-合并法”。但不能采用單軸疲勞分析方法，因?yàn)橹鲬?yīng)力方向隨時(shí)要發(fā)生變化。采用臨界平面法確定裂紋發(fā)生面。（參見(jiàn)V2、Section9）原理a、讀入彈性計(jì)算的FEA應(yīng)力張量（6個(gè)）；b、用載荷歷程數(shù)據(jù)乘以應(yīng)力張量，得到每個(gè)張量的歷程數(shù)據(jù)；c、計(jì)算表面主應(yīng)力的時(shí)間歷程（對(duì)接觸則考慮面外應(yīng)力）；d、從應(yīng)力計(jì)算三個(gè)主應(yīng)變的時(shí)間歷程數(shù)據(jù)；e、對(duì)應(yīng)變壽命疲勞分析(如Brown-Miller分析)，用多軸循環(huán)塑性模型將彈性應(yīng)力-應(yīng)變時(shí)間歷程數(shù)據(jù)轉(zhuǎn)化為彈塑性應(yīng)力-應(yīng)變歷程數(shù)據(jù)（對(duì)S-N曲線分析則忽略此步）；f、采用臨界平面法，對(duì)剪應(yīng)變、Brown-Miller分析分析，在三個(gè)可能的平面上計(jì)算剪應(yīng)變、法向應(yīng)變、法向應(yīng)力的時(shí)間歷程數(shù)據(jù)；對(duì)S-N曲線分析，則定義一個(gè)臨界平面，計(jì)算垂直于該平面的時(shí)間歷程應(yīng)力數(shù)據(jù)；g、計(jì)算每個(gè)臨界平面的疲勞損傷；各平面的每種疲勞循環(huán)用雨流計(jì)數(shù)法統(tǒng)計(jì)，總損傷由每種循環(huán)的損傷疊加而成；具有最短壽命的平面即為裂紋發(fā)生平面，其壽命寫入結(jié)果文件。h、修正耐勞極限。如果臨界平面內(nèi)的各種循環(huán)均低于耐勞極限，則沒(méi)有損傷產(chǎn)生；如任一循環(huán)有損傷產(chǎn)生，則自動(dòng)將耐勞極限降低25%，損傷曲線基于修正后的耐勞極限進(jìn)行估算。i、對(duì)每個(gè)單元都重復(fù)上述過(guò)程。疲勞載荷定義簡(jiǎn)單載荷（常幅值載荷）可用一系列數(shù)據(jù)定義或LDF文件定義；有以下兩種方式：有限元解為單位載荷的解，載荷歷程為實(shí)際載荷歷程；有限元解為最大載荷的解，載荷歷程為實(shí)際載荷與最大載荷的比值歷程；用戶界面方式定義一個(gè)BLOCK；增加載荷歷程數(shù)據(jù)或自定義載荷歷程數(shù)據(jù)，程序會(huì)自動(dòng)選擇一個(gè)dataset（計(jì)算結(jié)果），以組成aload*dataset定義。也可直接通過(guò)Add/aload*dataset定義。LDF文件格式略多載荷歷程(疊加載荷)多個(gè)有限元解集，每個(gè)解集對(duì)應(yīng)有一個(gè)載荷歷程；程序計(jì)算時(shí)先完成每個(gè)解集*相應(yīng)載荷歷程的應(yīng)力張量歷程計(jì)算；再完成所有載荷作用下的應(yīng)力歷程計(jì)算；最后完成主應(yīng)力計(jì)算?？稍诮缑嫦露x或LDF文件定義；界面定義方式定義一個(gè)BLOCK；通過(guò)Add/aload*dataset分別定義dataset、loadhistory；LDF定義實(shí)例INITEND#Blocknumber1BLOCKn=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=2,ds=2,scale=1END上述結(jié)果ds1對(duì)應(yīng)信號(hào)文件的1信號(hào)；ds2對(duì)應(yīng)信號(hào)文件的2信號(hào)。解集序列(datasetsequence)表示載荷隨時(shí)間變化的一系列應(yīng)力解，如柴油機(jī)曲軸在不同轉(zhuǎn)角時(shí)的應(yīng)力集序列。其中的每個(gè)結(jié)果集可以定義載荷歷程、循環(huán)次數(shù)、比例因子等參數(shù)?？稍诮缑嫦露x或LDF文件定義；ds=1-n形式可以定義很長(zhǎng)的結(jié)果序列。=1\*GB3①界面定義方式在一個(gè)Block內(nèi)只依次增加dataset，其順序表示先后順序。=2\*GB3②LDF定義方式INITEND#Blocknumber1BLOCKn=1,scale=1ds=1-2,scale=1或（ds=1ds=2）END可在一個(gè)Block內(nèi)定義，也可在多個(gè)Block內(nèi)定義。復(fù)雜loadingblock載荷定義由一系列的loadindblock組成，每個(gè)block含有解集序列或loadhistory*dataset，且可在每個(gè)Block上疊加高頻block。可在界面下定義或LDF文件定義；注：如果一個(gè)疊加載荷的歷程比其它載荷短，則默認(rèn)其長(zhǎng)度與最長(zhǎng)的載荷歷程相同，缺失部分用0填充。如果dataset沒(méi)有指定載荷歷程，則其只施加一次。如下例的ds3僅在信號(hào)的第一個(gè)數(shù)據(jù)點(diǎn)時(shí)被疊加。INITEND#Blocknumber1BLOCKn=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=2,ds=2,scale=1ds=3END載荷定義文件（LDF）含塊定義、有限元數(shù)據(jù)系列定義、載荷歷程定義、高頻載荷定義、溫度變化定義、時(shí)間定義。后兩者僅用于Fe-safe的熱-結(jié)構(gòu)分析。要求：至少有一個(gè)Block定義；每個(gè)塊以Block開始、以END結(jié)束；每個(gè)塊必須從第一列開始定義；每個(gè)塊內(nèi)可含多個(gè)結(jié)果集*載荷歷程；標(biāo)注行第一列為#。具體格式為：=1\*GB3①BLOCK定義BLOCKN=10,SCALE=1.2，dt=3600,TEMP=100上式中：N為BLOCK的重復(fù)數(shù)；SCALE為BLOCK比例，對(duì)BLOCK下所有項(xiàng)均適用；DT為N個(gè)BLOCK對(duì)應(yīng)的時(shí)間長(zhǎng)度；TM為BLOCK對(duì)應(yīng)的溫度。=2\*GB3②結(jié)果集序號(hào)定義DS=1，ES=2，SCALE=1，ESCALE=2.0DS為打開的有限元應(yīng)力結(jié)果集編號(hào)；ES為打開的有限元應(yīng)變結(jié)果集編號(hào)（僅對(duì)彈塑性有限元計(jì)算有效）；SCALE為應(yīng)力比例；ESCALE為應(yīng)變比例（應(yīng)變比例縮放后應(yīng)該用m/m單位，不用με單位）；此式中的scale、escale僅用于應(yīng)力、應(yīng)變的單位轉(zhuǎn)化。=3\*GB3③載荷歷程定義：與定義load*dataset等效lh=<load_history_filename>,signum=1,ds=7,scale=45式中：lh為載荷歷程文件；signum為載荷歷程文件的序列號(hào)；ds為結(jié)果集序號(hào)；scale為dataset的比例因子；=4\*GB3④定義時(shí)間高頻block與低頻block疊加時(shí)，或TMF、Turbolife分析時(shí)，需要定義時(shí)間?？赏ㄟ^(guò)上述的dt參數(shù)定義；也可通過(guò)lhtime參數(shù)定義；后者比前者優(yōu)先；兩者單位必須為秒。用lhtime參數(shù)時(shí)，時(shí)間點(diǎn)數(shù)必須與loadhistory定義取樣點(diǎn)數(shù)相同，若時(shí)間點(diǎn)數(shù)小于loadhistory定義的取樣點(diǎn)數(shù)，則會(huì)以最后一個(gè)時(shí)間點(diǎn)作為隨后的loadhistory取樣點(diǎn)的時(shí)間，出錯(cuò)；dt參數(shù)主要用于loadhistory取樣點(diǎn)的時(shí)間間隔都相同的情況，此時(shí)第一個(gè)取樣點(diǎn)對(duì)應(yīng)0s，最后一個(gè)取樣點(diǎn)對(duì)應(yīng)dt/N。當(dāng)loadhistory取樣點(diǎn)的時(shí)間間隔不相同時(shí)，lhtime的最后一個(gè)值等于dt/N（即block的一個(gè)循環(huán)時(shí)間）。最后一個(gè)取樣點(diǎn)和重復(fù)Block的零時(shí)間點(diǎn)的時(shí)間間隔被假設(shè)為0，因此必須將第一個(gè)取樣點(diǎn)的時(shí)間設(shè)置為非0值，以使lhtime最后一個(gè)時(shí)間點(diǎn)與重復(fù)block的零時(shí)間點(diǎn)間的時(shí)間差不為零。？？？示例1：每個(gè)lh文件定義的樣點(diǎn)數(shù)與lhtime定義樣點(diǎn)數(shù)相同。示例2示例3INITEND#Blocknumber1BLOCKn=10,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=2,ds=2,scale=1lhtime=0，1，2，3，4，6END#usingdtBLOCKn=1,dt=20Ds=1-5END#usinglhtimeBLOCKn=1,dt=20Ds=1-5Lhtime=4，8，12，16，20END=4\*GB3④疊加高頻載荷可疊加到定義的主塊載荷上，每個(gè)主載荷塊可疊加20個(gè)高頻塊。高頻載荷定義時(shí)采用loadhistory*dataset方式，循環(huán)區(qū)域從主block的開始位置到其結(jié)束位置；用dt確定高頻block的時(shí)間長(zhǎng)度，因此主block也要用dt或lhtime定義時(shí)間參數(shù)。示例1：主Block總時(shí)間為100s；HFblock時(shí)間為1s，重復(fù)100次。示例2：主Block總時(shí)間為100s；HFblock時(shí)間為1s，重復(fù)100次。BLOCKn=1,dt=100lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lhtime=1，2，3，4，5、6，7，8，9，10HFBLOCKdt=1Lhtime=0，0.1,0.3,0.8Ds=1-2Ds=7-8hfendENDBLOCKn=1,dt=100lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lhtime=1，2，3，4，5、6，7，8，9，10HFBLOCKdt=1lhtime=0，0.1，0.3，0.8lh=D:\fatigue\test.amc,signum=1,ds=1,scale=1Ds=7-8HFENDEND在hfblock中，lhtime的時(shí)間范圍不要覆蓋dt的時(shí)間范圍。如果主block在t=0時(shí)的幅值無(wú)定義（如通過(guò)lhtime參數(shù)），則程序會(huì)以主block內(nèi)lhtime定義的最后時(shí)間點(diǎn)的幅值代替它，以保證HFBLOCK能疊加在主block的所有時(shí)間范圍，一般地要避免這種情況出現(xiàn)（見(jiàn)207頁(yè)）。如果HFBLOCK內(nèi)的lhtime沒(méi)有定義，則HFBLOCK內(nèi)取樣點(diǎn)的時(shí)間數(shù)值為：0，dt/N，….，dt*（N-1）/N。N為dataset的數(shù)目或loadhistory中的樣點(diǎn)數(shù)目。HFBLOCK的dt值盡量取較小的數(shù)值。如下列例2比例1要節(jié)省計(jì)算時(shí)間。示例1：主Block總時(shí)間為35100s,重復(fù)1次；HFblock時(shí)間為20s，重復(fù)100次。示例2：主Block總時(shí)間為20s,重復(fù)1775次；HFblock時(shí)間為20s，重復(fù)100次。BLOCKn=1,dt=35100Ds=1，0.5Ds=1，0.5HFBLOCKdt=20Ds=2-4hfendENDBLOCKn=1775Ds=1，0.5Ds=1，0.5Lhtime=0，20HFBLOCKdt=20Ds=2-4hfendEND實(shí)例：懸臂軸的拉壓-彎曲-扭轉(zhuǎn)組合疲勞分析（疊加載荷）圖示懸臂軸方型大端外端面固定，小端端面承受載荷為：工況=1\*GB3①拉力F=200N/節(jié)點(diǎn)；工況=2\*GB3②切向壓力P=5Mpa；工況=3\*GB3③周向力Fx=2000N/節(jié)點(diǎn)。若工況=1\*GB3①、工況=2\*GB3②、工況=3\*GB3③的載荷同步作用，且其載荷歷程分別如文件Load01.txt的各列所示，試用ANSYS+fe-safe軟件計(jì)算：該軸能工作多少個(gè)循環(huán)？已知參數(shù)為：軸小端端面至大端內(nèi)側(cè)面距離L=200mm、半徑R=30m；大端厚度T=30mm、截面尺寸為50mm×50mm；軸與大端間過(guò)渡圓弧半徑R01=5mm。懸臂軸為鋼材，其材料牌號(hào)為SAE-0030，彈性模量E1=1.965×105MPa、泊松比ν=0.3；各工況作用載荷如圖1-3示。圖1工況=1\*GB3①載荷及歷程圖2工況=2\*GB3②載荷及歷程圖3工況=3\*GB3③載荷及歷程分析步驟1：在ANSYS環(huán)境采用N-mm單位建模型，作工況=1\*GB3①、=2\*GB3②、=3\*GB3③的應(yīng)力計(jì)算；見(jiàn)文件fatigue04.txt。計(jì)算結(jié)果見(jiàn)圖4-6。圖4工況=1\*GB3①應(yīng)力計(jì)算結(jié)果圖5工況=2\*GB3②應(yīng)力計(jì)算結(jié)果圖6工況=3\*GB3③應(yīng)力計(jì)算結(jié)果分析步驟2：在fe-safe環(huán)境的Opendatafiles窗口中，點(diǎn)擊右鍵打開載荷歷程文件Load01.txt。顯示有三個(gè)信號(hào)。在fe-safe環(huán)境的OpenFEModel窗口中，打開有限元計(jì)算結(jié)果文件（*.rst）、選種各Step的應(yīng)力、應(yīng)變按鈕，設(shè)置應(yīng)力單位為Mpa；在fatiguefromFEA窗口中，點(diǎn)擊loading按鈕，點(diǎn)擊OK；右擊Usinguserprofileloading**按鈕，選擇newloading，將后續(xù)載荷信息保存為一個(gè)*.ldf文件；點(diǎn)擊Add/block，選種OpenFEModel窗口中的應(yīng)力集/應(yīng)變集，在fatiguefromFEA窗口中的block下添加結(jié)果。依次添加Step為1（對(duì)應(yīng)工況1）、Step為2（對(duì)應(yīng)工況2）、Step為3（對(duì)應(yīng)工況3）的結(jié)果集。）點(diǎn)中Opendatafiles窗口中的信號(hào)1（#1），再右擊elasticblock下的stressdataset1，選擇addloadhistory，使應(yīng)力集1與載荷歷程1相結(jié)合。點(diǎn)中Opendatafiles窗口中的信號(hào)2（#2），再右擊elasticblock下的stressdataset2，選擇addloadhistory，使應(yīng)力集2與載荷歷程2相結(jié)合。點(diǎn)中Opendatafiles窗口中的信號(hào)3（#3），再右擊elasticblock下的stressdataset1，選擇addloadhistory，使應(yīng)力集3與載荷歷程3相結(jié)合。點(diǎn)擊OK、YES等。在Opendatabase窗口中，選中SAE-0030材料；在fatiguefromFEA窗口中，點(diǎn)擊material1*欄的表格值，指定材料為SAE-0030材料。點(diǎn)擊Yes。取材料的默認(rèn)算法計(jì)算；在fatiguefromFEA窗口中，點(diǎn)擊designlife按鈕指定規(guī)定壽命，進(jìn)行強(qiáng)度因子計(jì)算；在fatiguefromFEA窗口中，點(diǎn)擊probablity按鈕指定規(guī)定壽命下的存活概率計(jì)算；點(diǎn)擊analysis按鈕進(jìn)行計(jì)算。計(jì)算結(jié)果：最小壽命的對(duì)數(shù)為3.623，應(yīng)力因子為0.5；失效概率為100%。見(jiàn)圖7-9示。（注：fatige05.txt文件可產(chǎn)生帶殼單元的模型，可僅用殼單元作計(jì)算來(lái)比較相關(guān)結(jié)果）圖7壽命對(duì)數(shù)計(jì)算結(jié)果圖8強(qiáng)度因子計(jì)算結(jié)果圖9失效概率計(jì)算結(jié)果特殊應(yīng)用實(shí)例彈塑性有限元解的載荷定義讀入解的應(yīng)力、應(yīng)變結(jié)果。載荷定義要對(duì)每個(gè)計(jì)算結(jié)果集定義一個(gè)應(yīng)力-應(yīng)變對(duì)。Block的scale必須設(shè)置為1，不要縮放；但可用ds、es的scale進(jìn)行單位轉(zhuǎn)換（轉(zhuǎn)換后的單位應(yīng)分別為Pa、m/m）。需要時(shí)用biaxialNeuberRule進(jìn)行彈性-彈塑性修正。主應(yīng)變、BrownMiller應(yīng)變、最大剪應(yīng)變可用于彈塑性有限元結(jié)果的分析；一些實(shí)例：使用ds命令，不進(jìn)行彈性-彈塑性修正實(shí)例1實(shí)例2實(shí)例3BlockN=10，scale=1Ds=1，es=5Ds=2，es=6Ds=3，es=7Ds=4，es=8ENDBlockN=10，scale=11，52，63，74，8ENDBlockN=10，scale=1Ds=1-4，es=5-8END在上述三個(gè)文件中，有限元結(jié)果有4個(gè)，讀入結(jié)果時(shí)，1-4表示4個(gè)應(yīng)力集，5-8表示4個(gè)應(yīng)變集。熱-機(jī)械、高頻載荷定義實(shí)例要用到Fesafe/TMF模塊；載荷定義為時(shí)間、溫度的函數(shù)；因要估算樣點(diǎn)間的應(yīng)變速度變化率，因此，要確保T=0、DT時(shí)的應(yīng)力張量相同，避免出現(xiàn)較大的應(yīng)變速度變化率。時(shí)間(s)應(yīng)力集號(hào)溫度集號(hào)0112022503370449055BLOCKDs=1-5lhtime=0,20,50,70,90dtemp=1-5HFBLOCKdt=1ds=6,lh=\texi.txt,signm=1HFENDEND上例中，主block內(nèi)含有應(yīng)力定義（結(jié)果集1-5）、溫度定義（結(jié)果集1-5），hfblock內(nèi)有高頻定義，應(yīng)力集為6，對(duì)應(yīng)lh定義的載荷歷程；注意下例形式類型含義Blockds=1,lh=D:\fatigue\test_mcg2.amc,signum=1end結(jié)果集1及對(duì)應(yīng)載荷歷程，與下同。lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1BlockDs=1,0.5DS=1,0.5endblocklh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1endBlocklh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1lh=D:\fatigue\test_mcg2.amc,signum=1,ds=1,scale=1Ds=1Ds=2end實(shí)例：彈塑性有限元解的的疲勞分析圖示為在其中心位置開有一圓孔的薄板，薄板一端固定，承受載荷為：工況=1\*GB3①另一端拉壓力P1=30Mpa；工況=2\*GB3②下邊一半拉壓力為30Mpa。若將載荷變化情況工況=1\*GB3①-工況=2\*GB3②-工況=1\*GB3①定義為一個(gè)工作循環(huán)，試用ANSYS+fe-safe軟件計(jì)算：該薄板能工作多少個(gè)循環(huán)？已知參數(shù)為：薄板長(zhǎng)L=1000mm、高H=500m、厚T=5mm、圓孔半徑R=130mm。薄板左半部分為鋼材，其材料牌號(hào)為SAE-0030，彈性模量E1=1.965×105MPa、泊松比ν=0.3、屈服強(qiáng)度為σs=324Mpa，隨動(dòng)強(qiáng)化割線模量E11=5×104Ma；；薄板右半部分為鋁合金，材料牌號(hào)為AL1100-T6，彈性模量E2=6.89×104MPa、泊松比ν=0.3、屈服強(qiáng)度為σs=97MPa，隨動(dòng)強(qiáng)化割線模量E22=1×104MPa；分析步驟1：在ANSYS環(huán)境采用N-mm單位建模型，先作工況=1\*GB3