.fil文件的案例
FE-SAFE使用Abaqus的fil文件進行疲勞運算（一）
FE-SAFE使用Abaqus的fil文件進行疲勞運算（一）
FE-SAFE是一款高級疲勞耐久性分析和信號處理的軟件，它是多軸疲勞分析解決方案的領導者，算法先進，功能全面細致，是世界公認精度最高的疲勞分析軟件之一。
本文將通過一個實例（Tutorial 105），介紹使用Abaqus的fil文件進行疲勞運算的過程。
準備階段：啟動fe-safe之后，首先需要設置工作目錄。同時，我們需要對上次所留下的數據進行清理，點擊Tools>>Clear Data and Settings…，選中所有的選項，點擊OK，數據清理完成。。。。。。
接下來的步驟附件中會有詳細說明，大家可以先看看，相互交流，相互學習，有什么好資料好想法希望大家都可以拿出來分享下!
FE-SAFE使用Abaqus的fil文件進行疲勞運算.pdf
展開 [分析處理] ABAQUS生成fil文件詳細步驟 inp轉換fil
ABAQUS結合疲勞分析軟件fe-safe或者fatigue進行疲勞分析，需要導入ABAQUS分析結果生成的文件fil，由于ABAQUS軟件默認生成odb格式的計算結果文件，而軟件fe-safe雖然能導入odb但是效果很可能不太理想，于是如何生成fil格式顯得很重要。
論壇上以有許多人給出解釋說明，但無奈不夠完整，講述十分模糊。
本人在大概學會轉換過程后，進行了總結：大致只需要在inp文件末尾添加一段代碼即可將inp轉換為二進制fil格式。
記事本打開分析結果生成的inp文件，加入代碼如下：
*Output, history, variable=PRESELECT
--------------以下為添加代碼
*El Print, freq=1
*Node Print, freq=1
*NODE FILE
CF,
RF,
U,
*EL FILE, POSITION=NODES
S,
SINV,
---------------
*End Step
打開ABAQUS軟件，新建一個job，你在新建任務的時候，Souce 要選Input File（輸入文件），然后找到修改過的那個INP文件，提交運算，運算結束會自動生成fil文件。
如果不這樣，按照原來的方法點提交運算，這樣ABAQUS又自動生成一個跟原來一樣的INP文件，把修改的覆蓋掉了。
一般的應力問題都可以成功轉換。
如果要
用 patran 來讀取的話，可能需要在step 前加入*FILE FORMAT,ASCII
皇者謝霆鋒 2015年01月06日
展開 ABAQUS 結果文件輸出
結果文件有output database（.odb）文件和result file（.fil）文件。區別是
a.odb文件默認生成，是二進制文件，abaqus才能訪問，python后處理也是通過abaqus訪問。
b.fil文件不是默認生成，需要在inp中添加關鍵字（*FILE FORMAT, ASCII），添加的位置可以是（*Restart）命令之后。是文本文件，文本編輯器可以訪問。
要輸出的量可以通過cae界面選擇，也可以通過inp添加。
輸出到fil文件屬于歷史輸出范圍，輸出的是所有增量步的值。fil文件的大小似乎是固定的，在cae界面中，歷史輸出可以設置只輸出最后一個increment的值，對應在inp文件中表示為（*Output, history, frequency=99999），輸出頻率非常大，故只輸出最后一幀。這個命令針對fil的輸出無用，這里姑且認為fil只能輸出的所有增量步的值。
*NODE FILE 輸出節點變量，位移（UCOORD-101），節點坐標（COORD-107）
*EL FILE 輸出單元變量，MISES應力（S），應力分量（SINV-12）
要使用URDFIL子程序，首先必須先創建.fil文件，子程序讀取.fil文件，Fortran代碼實現結果的格式化輸出。
展開 ABAQUS fil 文件如何閱讀
fil 文件是一條一條的，各條記錄之間用 *I 隔開，到下一個 *I 是一條完整的記錄
每條文件內部用大寫的 I 隔開，每個數據第一個數表示位數，緊跟著后面表示該數據，如211I，表示該數據有兩位，是11；如41900I，表示4位數的1900
每條數據第一個數據段表示該條數據的數據量，如
*I 211I 11I 11I 11I 10I 10A I 13I 13I 10I 10，第一個211I表示該條數據有11個數據；
*I 18I 211D 5.143470386335691D+00D 3.097019956138536D-03D-4.272948703247792D-04D-4.131188576751562D-03D-7.064966453604931D-04D-3.626273143709050D-04，第一個數據18I表示該條數據有8個數據
每條數據第二個數據段表示該條數據對應的輸出量，如
*I 211I 11I 11I 11I 10I 10A I 13I 13I 10I 10，第二個11I表示這條數據對應的RecordKey為1的輸出，http://www.wufengyun.com:888/v6.14/books/usb/default.htm?
展開 
abaqus與fatigue結合疲勞分析
圖2 應力云圖
2、*.Fil文件說明
*.fil文件是ABAQUS的一種二進制輸出文件，供其他軟件（如Patran）后處理使用，如生成X-Y曲線，制作二維表格等，可以輸出的項目包括：單元、節點、接觸面、能量、模態、梁截面等的輸出信息，輸出的方法是在INP文件中增加輸出指令，
生成*.fil文件的步驟如下
對ABAQUS/Standard，可以直接輸出.fil文件，步驟如下：
在inp文件中，step步驟之后， end step步驟之前，加上以下內容：
*NODE FILE
RF,U,V
**輸出節點的作用力(RF)，位移(U,V)到*.fil中
*EL FILE
S,E
**輸出單元應力(S)，應變(E)到*.fil中
在abaqus的job界面重新運行inp文件，即可得到對應的fil文件
3、疲勞壽命估算
疲勞壽命估算需用到MSC.Patran軟件中的MSC.Fatigue模塊。如圖3所示，MSC.Fatigue位于MSC.Patran的Tools菜單下，點擊Main Interface即可進入MSC.Fatigue模塊主界面。
圖3 在MSC.Patran中進入MSC.Fatigue界面
對結構施加的疲勞載荷譜見表1。
表1 名義應力譜
級數
Smax
Smin
循環次數
1
318
-121
2
2
176
41
982
將載荷譜導入MSC.Fatigue后顯示如圖4所示。
圖4載荷譜塊譜示意圖
將模型的結果文件（.fil文件）導入MSC.Patran中，點擊輸入材料和載荷譜信息，進行壽命估算，得到模型的對數壽命云圖，如圖5所示。
展開 Abaqus 非線性屈曲分析方法
具體做法如下（注意關鍵字的插入位置和書寫格式）：
再提交計算后會生成相應的.fil文件，該.fil文件將在后屈曲分析中進行引用！
2
后屈曲分析
后屈曲分析通常在線性屈曲分析后，通常的做法是將原線性屈曲模型復制生成新的模型，并調整和修改分析步和載荷工況、接觸等，比如方筒的壓潰，就需要修改為顯示動力學分析，添加自接觸關系等。當然最終要的是要引入缺陷擾動，即調用之前生產的.fil文件。具體做法如下：
其中關鍵字中file后面的名稱即調用的前面的.fil文件的名稱，第二行第一個數字1表示引入的是第一階線性屈曲的擾動結果，0.5e-3表示引入的擾動量的大小。
擾動量（或稱為缺陷因子）的大小準確的做法是進行試驗矯正，一般是按照經驗的做法，取殼厚或桿長的 1‰ ~ 2‰ 。
實際上一定程度上，你也可以調整擾動量的大小進行試算并比對后屈曲狀態和其他計算結果參數，如果結果變化不大，說明結構對缺陷的敏感性較低，反之表示相對敏感。當然這個觀點還有待深入考究。
上圖中左圖沒有考慮引入擾動量，右圖是考慮的擾動量的結果。大家可以仔細觀察下結果的區別。
實際上右圖在壓潰的過程中，相對光滑，結果相對更加合理。
當然對于方筒這類實際上是通過顯示方法實現的，更準確的講是動力屈曲分析，所以我們還得判斷動能、塑形耗散等能量參數，才能是結果更加準確，限于偏于以后單獨在顯示動力學分析中再單獨討論。
本文轉自 極速有限元 旨在分享知識，侵刪。
展開 Abaqus 非線性屈曲分析方法 附ABAQUS分析手冊分析卷下載
具體做法如下（注意關鍵字的插入位置和書寫格式）：
再提交計算后，會生成相應的.fil 文件，該.fil 文件將在后屈曲分析中進行引用！
02
后屈曲分析
后屈曲分析通常在線性屈曲分析后，通常的做法是將原線性屈曲模型復制生成新的模型，并調整和修改分析步和載荷工況、接觸等。比如方筒的壓潰，就需要修改為顯示動力學分析，添加自接觸關系等。當然最重要的是要引入缺陷擾動，即調用之前生成的.fil 文件。具體做法如下：
其中，關鍵字中file 后面的名稱即調用的前面的.fil 文件的名稱，第二行第一個數字1表示引入的是第一階線性屈曲的擾動結果，0.5e-3表示引入的擾動量的大小。
擾動量（或稱為缺陷因子）的大小準確的做法是進行試驗矯正，一般是按照經驗的做法，取殼厚或桿長的1‰~2‰。
實際上，在一定程度上你也可以調整擾動量的大小進行試算并比對后屈曲狀態和其他計算結果參數，如果結果變化不大，說明結構對缺陷的敏感性較低，反之表示相對敏感。當然，這個觀點還有待深入考究。
上圖中左圖沒有考慮引入擾動量，右圖是考慮擾動量的結果。大家可以仔細觀察下結果的區別。
實際上右圖在壓潰的過程中，相對光滑，結果相對更加合理。
當然，對于方筒這類實際上是通過顯示方法實現的，更準確的講是動力屈曲分析，所以我們還得判斷動能、塑形耗散等能量參數，才能使結果更加準確。
下載地址：ABAQUS分析手冊分析卷
展開 hypermesh-abaqus視頻教程
6.如何把abaqus生成的.fil文件翻譯成Hypermesh/hyperview可以識別的.hmres格式；
7.如何使用hypermesh和hyperview做后處理。
視頻中的模型描述：
模型是由實際工程問題加以簡化和抽象得到的。經過簡化后，模型比較簡單，但是又能夠說明問題。
模型由兩個tower組成,在法蘭處用四個螺栓聯結起來。螺栓上施加有預緊力。tower的底面約束住。在tower的頂端施加了一個彎矩Mx.
視頻教程中所用的軟件版本：
Hypermesh v8.0
Hyperview v8.0
abaqus v6.8.1(你使用低版本也可以的)
視頻基本信息：
1.視頻文件格式：可直接播放的.exe文件，無需任何播放器。
2.文件大小：83M.
3.播放時間：大約40分鐘。
Tower.part2.rar
Tower.part1.rar
展開 fe-safe中用s-n曲線預測壽命實例
fe-safe預測疲勞壽命的原理是對比材料的S-N曲線，一般情況下，fe-safe的材料屬性里會包含一個S-N曲線的選項，如下所示：
在使用fe-safe進行分析時，將有限元分析的應力結果（比如odb或者fil文件）導入到fe-safe中，進行處理后與S-N曲線進行比對，從而得出構件的疲勞壽命，接下來就詳細的介紹下這個流程。
一、對構建進行有限元分析，求出靜力結果
二、打開fe-safe，清除所有痕跡
三、載入有限元分析結果(這里載入的是odb文件)
四、創建一個荷載譜（這里以標準正弦波為例）
五、在Ladling settings里設置荷載譜
六、在Analysis setting里配置材料等計算所需項目，然后點擊analysis進行計算
七、將分析后的結果文件導入有限元軟件里進行后處理
此時即可顯示構件的壽命云圖
fe-safe中用s-n曲線預測壽命實例.pdf
展開 加筋板屈曲riks法分析實例
算例來源如下圖所示：（對于這種沿著壓載方向剛度均勻的結構來說，不能直接進行riks分析，必須先進行擾動分析）
先進行線性屈曲特征值分析，inp文件請查看附件，僅計算第一階模態，模態的構型如下圖所示：
求得特征值為5.9655，施加的線壓力載荷為100，故屈曲線壓力載荷為596.55，總的屈曲力為32214，第一圖文獻1的給出的屈曲載荷為29160，可見誤差不大，可以接受。
在線性屈曲分析的基礎上，進行riks后屈曲分析，inp文件請查看附件，分別引入初始缺陷為shell厚度的1%、10%、50%、100%進行比較，所得到的整理后的載荷比例因子—位移如下圖所示：
可知：1）初始缺陷的引入對后屈曲行為有較大的影響；2）此壁板結構在初次屈曲失穩后仍可繼續加載（如1%情況），且未出現負剛度。
額外之言：如果僅進行線性特征值分析，所有操作均可在CAE中完成；若想在其基礎上繼續進行后屈曲分析，則需要手動更改inp文件添加*NODE FILE語句，以生成后屈曲分析所需要的fil文件；后屈曲分析也需要在inp文件中手動添加*imperfection語句，詳細情況請參看分析手冊。所以本帖附件沒有給出CAE文件，僅給出inp文件。希望此小例子能給大家帶來幫助！
buckling.rar
BUCKLING-CAE.rar
展開 薄殼圓柱屈曲分析
2 操作步驟
a 建立模型
這里用到薄殼模型，建模步驟我就省略了
賦予材料的話，要給于塑形，各種參數見上圖
b 建立相關連接
先建立兩個參考點（0，0，0）和（0，0，500）
再將頂部的所有點、底部的所有點與相應的參考點進行耦合（剛性）
c 建立分析步
d 加載邊界條件
由于與力值大小無關，因此在頂端參考點給-1的力，底部參考點完全約束
e 修改inp活著關鍵字，將節點的位移量進行輸出
f 計算求解，確保你的結果里面有fil文件（網格自行劃分，種子數給的10）
注意到屈曲力大概在1.2e6這么大
g 復制model-1到model-2
h 刪除原來的分析步，建立riks分析步
i 加載邊界條件
這里的力值參考最大屈曲力1.2e6，但是實際上這個值需要嘗試，特別是加載imperfection關鍵字之后，力值往往會急劇見效，這里設置為1.2e5
j 加入自定義關鍵字imperfection
在step之前加入
這里的imperfection是將特征值屈曲時的模態按照組合進行初始化，更接近實際結果
k 計算
展開 
ABAQUS非線性屈曲分析
我們此次課程中采用屈曲分析方式，先計算屈曲模態，也就是先做特征值屈曲分析，此分析為線性屈曲分析，在小變形的情況下進行，得出臨界載荷（一般取一階模態的eigenvalue乘以加載的單位載荷1），且需要在inp文件中輸入如下圖字符，輸入次字符的目的是將初始缺陷的節點輸出為.fil文件；然后將1階屈曲模態做為初始缺陷引入極限載荷后屈曲分析，后屈曲分析可以定義非線性材料及幾何非線性，所以risk屈曲分析也成為非線性屈曲分析.
分析模型，為一塊薄長板
特征值屈曲分析
后屈曲分析
結果對比
（1）特征值屈曲分析
（2）引入初始缺陷的后屈曲分析
（3）直接加載極值載荷做屈曲分析，也就是我們開始介紹的第一種方法，對比結果看差異很大
動圖結果
更多仿真分析案例及相關知識技巧，歡迎訪問“元王仿真云”（www.featechweb.cn）——集仿真項目咨詢，CAE學習培訓，設計工具應用，仿真人才招聘、仿真新聞資訊于一體的CAE仿真行業垂直服務平臺。服務仿真，助力設計，專業為你！
深圳市有限元科技有限公司（簡稱元王）成立于2007年，是一家以計算機輔助工程CAE(Computer Aided Engineering)為主業，以工程仿真軟件開發為核心，集CAE咨詢、CAE培訓、CAE軟件研發與銷售為一體的高科技企業。 依托深厚的CAE技術背景和工程經驗，為客戶提供高水平的CAE工程咨詢服務，為企業創造價值。聯系方式13632683051（微信同號）
展開 常用結構疲勞分析軟件的比較
疲勞分析軟件一般作為有限元軟件的后處理來進行結構的疲勞分析，需要把有限元的結果文件導入進行分析，和目前廣泛使用的用限元軟件的良好匹配是很重要的；疲勞分析軟件是否能為工程師提供滿意的解決實際問題的工具是一個十分關心的問題。所以在這里主要對比的內容是：接口和功能部分。操作和可視化也有所提及。至于價錢不在分析的范圍之內。
本部分是關于上述軟件與其它有限元軟件接口的比較。
1.1 Nsoft軟件的接口
在Nsoft軟件中，可以直接將有限元元軟件的計算結果直接倒入到軟件中，并能直接讀出模型的材料和單元分組、以及應力的信息。并可以將導入的文件顯示出來
Nsoft軟件中，可以直接導入的有限元軟件及文件格式如下
Patran/nastran的文件形式.op2
I-DEAS的文件形式.unv
ANSYS的文件形式.rst
Pro/Mechanica的文件形式.neu
ABAQUS的文件形式.fil
1.2 FE-safe 軟件的接口
FE-safe 軟件可以直接讀入有限元軟件的分析數據，例如應力、應變和溫度等。
可以直接讀入的軟件的結果文件格式如下：
ABAQUS軟件的格式為：.fil文件、
ANSYS軟件的格式為：.rst文件
NASTRAN軟件的格式為：.f06文件、.op2文件 `
Pro/Engineer軟件的格式為：.s01文件、s02文件等文件（應力）.d01文件 （溫度）
SDRC I-DEAS軟件的格式為：.unv文件
1.3 WinLIFE 軟件的接口
WinLIFE 軟件的合作者是NASTRAN(MAS.visual NASTRAN4 Windows)和NE/NASTRAN。這兩個軟件是采用FEMAP作為前處理和后處理的。這就決定了WinLIFE所能識別的數據的結構。
展開 基于Hypermesh和Abaqus的子模型分析
在輸出inp文件之前，使用Hypermesh的Organize功能將網格放在同一個Component中。
注2：Component的名稱不能出現數字，否則導入到Abaqus中將可能出現沒有網格信息的情況。
2）子模型屬性設置和子模型邊界條件設置
子模型屬性設置中需輸入全局模型的結果文件名。需要讀入的全局模型結果文件包括.odb文件（或.fil文件）和.prt文件，這些文件都應該在該計算模型所對應的工作路徑下。
上述子模型邊界條件設置的含義為：子模型邊界上的驅動變量為位移U1、U2和U3，讀入全局模型中第一個分析步的位移結果。對于同一個分析步，全局模型和子模型的增量步長可以不一樣，Abaqus會自動對其進行插值處理（對于大變形分析也沒有問題）。
建立第一個子模型時仍然使用全局模型在此部位原有的網格，將子模型的分析結果和全局模型的分析結果相比較，如果二者幾乎完全相同，就可以初步驗證子模型的各項參數是正確的。
然后再對子模型進行網格細化，得到更精確的分析結果。
子模型分析包括以下幾個基本步驟。
（1）完成對全局模型的分析，并保存子模型邊界附近的分析結果；
（2）創建子模型，定義子模型邊界；
（3）設置各個分析步中的驅動變量；
（4）設置子模型的邊界條件、載荷、接觸和約束；
（5）提交對子模型的分析，檢查分析結果。
展開 HyperMesh 與abaqus的接口問題
定義Load Steps
首先為了在HyperView和HyperMesh中的后處理，需要Abaqus計算出.fil文件。在BCs里面的output block中創建一個輸出控制，要求輸出節點位移、單元應力、約束反力和接觸應力等相關結果，各參數的含義請參考Abaqus手冊。在本例中，創建一個名為standard的output block。
在load steps中創建一個step1，在loadcols中選擇constraints和forces，在outputblock中選擇剛才創建的standard。
在編輯step1的時候，需要為其指定如下的step parameters。
Name——step的名稱
Nlgeom——幾何非線性
在AnalysisProcedure中指定這次分析為Static。
其它選項可以視情況而定。
輸出*.inp文件
在HyperMesh的file面板中輸出*.inp文件。該文件可以直接用于Abaqus計算。
1如果不做動態分析，不需要輸入材料密度。
展開 