abaqus 場變量的案例
ABAQUS案例-場變量的應用及材料彈性模量隨場變量而變化 ￥3
ABAQUS中的場變量具有較高的應用價值，可以在一些復雜的工程應用中極大的減輕工作量。本實例即是展示一個場變量應用——材料彈性模量隨場變量而變化，其中它涉及到關鍵字的編輯（關鍵字的具體編輯也在附件中）。本實例在附件的inp文件中。
ABAQUS導入初始場變量（預定義場）多次低速沖擊以及沖擊后壓縮 ￥38
ABAQUS導入初始場變量（預定義場）
通常利用ABAQUS計算時，需要多步驟分析，例如計算多次低速沖擊以及沖擊后壓縮等，下面詳細描述利用數據傳遞方法進行多步驟分析。（建議購買視頻，視頻內包含此帖子）
導入效果圖如下：
導入的損傷云圖
導入的應力場
導入的位移場
分層損傷的導入
1. 計算完成后，新建一個ABAQUS 窗口，切記與上一步計算的ODB文件在同一個文件夾下，導入Part部件
abaqus后處理插件—場變量結果標注 ￥45
abaqus后處理界面中探針功能附帶的標記樣式非常丑陋，基于abaqus的試圖注釋功能進行二次開發，形成了場變量標注插件，方便快速的標注關心區域的應力應變等結果。
插件介紹：
按鈕介紹
從左至右依次是：標記按鈕、隱藏標記按鈕、恢復顯示按鈕、刪除按鈕
示意動畫
使用方法：
1） Probe查看節點結果，并勾選需要標記的節點項；
2）點擊工具欄中的標記按鈕，進行標記。
特點
1） 標記速度快，即使在單元數目達到百萬級及以上的模型中，標記速度仍無明顯延遲；
2）所有標記注釋均在試圖注釋功能界面里，有利于對美觀度有更高要求者進一步修改美化。
展開 Abaqus后處理插件——云圖場變量查詢器 ￥50
基于Python對Abaqus進行后處理的二次開發，形成了場變量查詢插件，寫論文寫報告的實用小工具，話不多說直接上圖。
其實Abaqus后處理模塊本身具有場變量查詢功能（Tools-Query-Probe values），可以滿足使用需求，但其顯示效果非常不友好，以節點（單元）編號+場變量值的形式顯示，理工科氣息十足。不得不說，Abaqus在操作體驗方面而言和Workbench還是有一些差距。
PS：常規的二次開發通過fieldvariable.getsubset方式查詢到的場變量為節點直接插值的計算結果，即節點上的值是從相鄰單元插值后未經過平均化的值，是離散的，而結果顯示云圖默認是按75%平均化處理后的結果，更接近實際。
本插件查詢到的場變量為平均化之后的結果，和云圖的顯示的值一致。
展開 
Abaqus_Python 后處理結點數據提取、場變量云圖輸出 ￥1.99
今天體會了Python腳本操作ODB文件進行數據的輸出、處理等方面的快速高效。
首先介紹一下背景，同事的項目，計算case較多。每個case需要輸出幾個路徑上的部分節點值，而且需要針對不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節點數據，工作量會非常大，費時會遠遠超出計算所用時間；并且，同事還要求我根據case名稱、組件名稱（模型有多個組件組成，且網格是在part上畫的，所以同一個節點編號可能對應多個節點）等保存成不同的文件。
結點溫度數據沒有保存，因此沒有輸出。下圖是自動輸出的png格式云圖。
比較倉促，代碼也比價粗糙。回頭做一個詳細的介紹。
主要代碼見下。
ABAQUS案例—邊坡穩定性分析及場變量在邊坡強度折減中的應用 ￥3
本案例（附件中的inp文件）介紹了如何采用ABAQUS軟件進行邊坡穩定性分析，以及介紹了場變量在邊坡強度折減中的應用。介紹了采用平面應變單元來模擬三維的邊坡穩定問題所需要注意的問題及分析技巧。
一維單元模擬混凝土構件開裂的解決方案
圖1 鋼筋混凝土結構裂縫
本案例使用ABAQUS對一根鋼筋混凝土受拉構件進行裂縫估算分析，根據鋼筋應力狀態計算等效裂縫寬度，并采用Python腳本在ODB結果文件中創建裂縫場變量，實現在ABAQUS中使用一維單元進行快速分析，并在后處理模塊 顯示等效裂縫寬度的目標。
案例涉及的相關技術：
①ABAQUS梁單元Rebar積分點插入；
②利用Python腳本提取ABAQUS場變量數據；
③利用Python腳本創建ABAQUS場變量數據。
計算報告編寫采用操作引導式，希望能為讀者使用ABAUQS場變量創建提供有益參考。操作分析要點為：
①ABAQUSABAQUS梁單元Rebar積分點插入；
②ABAQUS場變量輸出Python腳本getSubset()函數應用；
③ABAQUS場變量編輯Python腳本addData()()函數應用。
二、計算任務
1.模型裝配及接觸連接
計算模型取自《混凝土結構：混凝土結構設計原理》（第六版）習題8-3。
計算模型為鋼筋混凝土屋架下弦按軸心受拉構件，見圖2。模型只包含1個part。截面寬200mm，截面高160mm。因為案例模型較為簡單，混凝土梁采用B21單元模擬，鋼筋通過在與混凝土單元共節點建立鋼筋箱型截面單元實現。混凝土材料為C40，鋼筋為HRB400，為實現材料單軸滯回損傷模擬，混凝土受壓骨架線上升段采用Hognested曲線，下降段為直線；鋼筋均采用《混凝土結構設計規范》（GB 50010-2010）附錄C中鋼筋本構關系，材料本構示意圖見圖3，材料參數見表1。
展開 Matlab“稍作修改”Abaqus-odb結果！！！
; meshplot(parts(2)); title(parts(2).Name);
Matlab繪制部件模型
接下來是重中之重了，也就是如何繪制場變量云圖？
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
場變量云圖，如果我們想更改，那就可以在該區域對應的Excel部分進行修改，即可操作成功。
疲勞分析|Abaqus Goodman方法案例操作 附ABAQUS疲勞分析簡介下載
Abaqus/View結果讀取
讀取分析歷程中的最大交變應力和最小交變應力云圖
新建場變量：Alternating Stress和Mean Stress
根據公式：
在Abaqus后處理新建場變量
輸出場變量值到Excel
針對新建場，輸出單元積分點對應的交變應力和平均應力，并輸出到Excel，與Goodman圖一并繪制。
上圖，
仿真所得單元積分點落到
曲線的上方或下方，
處于上方為疲勞壽命沒達到
臨
界曲值
10
E5
次。
下載地址：ABAQUS疲勞分析簡介
展開 【轉載】ABAQUS之inp文件結構解析
** 在Dload和DSload中的載荷類型的代碼，從ABAQUS Analysis User's Manual的Distributing loads查**找。下例的P代表均布面荷載。**
** LOADS
**
** Name: Load-1 Type: Pressure
*Dsload
_PickedSurf4, P, -100.
**
** OUTPUT REQUESTS
*Restart, write, frequency=0
**不輸出用于重啟動分析的數據。
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**將Abaqus默認的場變量寫入ODB文件。這里可以去field output manage中改。也可以：
**在第一個分析步里，有如下語句：
*Output，field，**variable=PRESELECT 在其后添加如下語句：
*NODE PRINT，NEST=Set-Head-Ref
*RF，
** HISTORY OUTPUT: H-Output-1
*Output, history, variable=PRESELECT
**將Abaqus默認的歷史變量寫入ODB文件。
*End Step
展開 
inp書寫規則（以一個inp文件為例）
**----------------------------------------------------------------------------------------------------------------------------------------------------
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**----------------------------------------------------------------------------------------------------------------------------------------------------
**將Abaqus默認的場變量寫入ODB文件。
展開 【iSolver案例分享64】一對集中力作用下受壓大變形圓環的理論公式、iSolver和Abaqus結果對比
網格收斂性考察表
3 結果與討論
初始破損載荷
首先，依據de Runtz和Hodge提出的理論公式，對本文中的圓環結構進行了估算：
下表展示了理論公式結果、iSolver模擬結果和Abaqus模擬結果的對比。三者結果高度接近，相互印證了計算結果的準確性。值得注意的是，iSolver在模擬中計算出了比Abaqus更接近理論公式的結果，與理論公式之間的誤差僅為0.51 %，表現尤為出色。
初始破損載荷對比表
結構大變形毀傷特征
下圖展示了iSolver和Abaqus在不同場變量下的計算結果對比。通過觀察，可以發現兩者模擬出的毀傷特征和典型位置幾乎一致。iSolver較好地模擬出了薄壁結構的毀傷特征，與Abaqus結果一致，兩者在不同場變量的分布上均表現出極佳的一致性。對于該圓環結構，需要四個塑性鉸來形成破損機構。這一點上，模擬結果也再現出了實驗現象。
文獻中的實驗現象
模擬結果對比圖
接下來，我對一些典型數值進行了統計。可以看出，在各種物理量的計算結果中，iSolver與Abaqus之間的平均誤差僅為0.06%，幾乎完全吻合。
展開 建筑結構動力彈塑性與倒塌分析的參數化建模軟件PA-TRANS
根據所采用的材料參數以及3.2.2節中所選材料本構模型，PA-TRANS自動設定ABAQUS模型中所用材料的材料定義參數。這些材料定義參數可以在轉換后的ABAQUS模型中繼續予以修改。若用戶選擇其他本構模型（非本程序提供或ABAQUS自帶的本構模型），則其材料參數在轉換模型中未給出，必須在轉換后的ABAQUS模型中予以人工補充。
圖3.3 材料參數設置頁面
4. 計算模型設置頁面
此頁面主要包含計算模型的總體轉換與輸出結果信息設置，如圖3.4所示。
圖3.4 計算模型設置頁面
1) 總體轉換信息欄
總體轉換信息主要包括：是否將荷載轉換為質量、是否考慮梁柱桿端的非剛性約束、是否考慮剪力墻的邊緣構件、是否轉換剛性樓板（對于板柱結構，必須轉換樓板）、是否考慮結構的幾何非線性以及是否轉換配筋。
對于帶有剛性構件的結構，需選擇設置剛性桿件剛度的處理方式（剛度增大法或剛體約束法）。程序默認采用剛度增大法，此時需設置剛度增大系數。對于帶有速度型阻尼器的結構，需選擇設置速度型阻尼器的計算模型（Kelvin模型或Maxwell模型）。程序默認采用Kelvin模型，即阻尼與剛度元件并聯模型。
2) 場變量結果輸出信息欄
設置在ABAQUS中的場變量輸出結果選項，包括結構位移與節點內力、速度、加速度、應力應變與材料塑性損傷等結果。
3) 歷史變量結果輸出信息欄
設置在ABAQUS中的歷史變量輸出結果選項，包括結構總能量結果以及各類構件的能量結果。
4) 其他信息
同時，此頁面還包含設置在ABAQUS中的單元網格劃分尺寸、墻板單元截面積分點數、地震分析工況結果輸出的時間間隔以及參與分析的CPU數量等信息。
5.
展開 