源問題

abaqus中，對一圓柱管材沖壓脹形，結果中，怎么提取管材的厚度變化曲線？

問題分析

根據問題描述，所述圓柱管材應該采用實體建模，其厚度變化無法直接獲得，一種可行的思路是在圓柱管材內外壁找到兩個特征節點，以這兩個特征節點之間的距離表示管材的厚度，從而得到厚度變化曲線。

ABAQUS中不支持輸出變形后的節點坐標（當前使用版本6.12，后續更新版本未知），因此本問題的關鍵在于如何獲取變形后的節點坐標。

如何獲取變形后的節點坐標？

目前從網上查找的可行的方法有三種。

方法一

在visualization模塊中，通過query-general queries-Nodes工具，逐個選取查詢變形后的節點坐標（Deformed coordinates (unscaled)），手動記錄后再用其他工具繪制進行處理，得到想要的結果。

此方法適用于較少的節點坐標獲取需求，若需要獲得的節點數量很多的話，工作量不言而喻。

方法二

先完成計算，得到計算結果的odb文件；

新建一個model，在abaqus菜單中 file-import-part，類型選擇odb，導入你的odb，然后在出現的對話框下面選中“import deformed configuration” ，step、frame根據你的需要選擇，此時模型上的節點坐標即為變形后的坐標。

此時輸出的inp文件中的坐標即為變形后的坐標。如果需要的是特定節點的坐標，則可以新建一個節點set，然后插入關鍵字（插入位置為end STEP的前一行）

*Node?print,?nset=nsetname（需要輸出的節點集合名）

coord

這樣就可將所關注的節點坐標寫入到dat文件，使用其他編程工具讀取進行后處理。

此方法整個思路其實相對容易理解，但是操作過于繁雜，而且每次只能提取一個增量步的變形節點坐標，輸出數據還需要編程讀取處理，比較麻煩。

方法三

由于輸出的常變量中沒有變形后的節點坐標，因此可采用Python語言提取節點坐標（original coordinate）和特定分析步-增量步的變形量，通過運算容易得到變形后的節點坐標。

示例：

根據題目描述，建立一個厚度5mm、高度50mm的圓柱管，在上表面施加-10mm的強制位移約束，在圓柱管（part上，非assembly）內外相對應位置分別建立一個單節點的節點集，分別命名為SET-IN、SET-OUT，采用general static分析步，設置分20個增量步求解（便于繪制曲線），其余為常規操作。

本算例腳本已進行升級，此版本不再提供購買下載（原貼為收費改不過來了），請勿購買！！！

新版代碼及使用說明詳見：

http://www.yqgqt.org.cn/content/post/0c63a829-4f12-43c9-970a-070e77481906