abaqus 節點位移的案例
ABAQUS-如何求結構的節點位移單元應力分量和支反力
一.問題
如圖1所示，顯示了四根桁架結構的尺寸與約束情況，材料為鋼，彈性模量設置為2.96Gpa，橫截面積為100，求該模型的
節點位移、單元應力分量、支反力。
圖1
二.部件與材料
首先按照圖1創建部件，選擇二維平面，特征為
線，繪制相應的草圖，并生成實體，命名為link。
圖2
如圖3所示，在屬性模塊創建材料屬性，選擇力學-彈性，在彈框中填寫楊氏模量2960，泊松比0。
圖3
如圖4、5、6所示，創建截面，選擇類型梁-桁架，并賦予材料屬性，填寫截面面積為100。
圖4
圖5
圖6
三.裝配與分析步
如圖7、8所示，將部件進行裝配，創建靜力通用分析步。
圖7
圖8
四.邊界條件與載荷
如圖9所示，創建邊界條件，選擇轉角/位移，約束點1與點2U1和U2的位移；約束點3U2的位移。
圖9
如圖10 所示，創建載荷，選擇點4，給定集中力CF2=-150N。
圖10
五.網格與作業提交
選擇單元族為桁架，單元類型為T2D2，提交作業。
圖11
圖12
六.結果展示與后處理
圖13是放大之后的變形圖，可以看到點4向右下方位移。
圖13
對結果進行處理，提取模型節點位移，單元應力分量和支反力。首先點擊進入可視化模塊，依次點擊選項-通用，在彈框中選擇標簽，勾選顯示單元編號和顯示節點編號，在右側可以自行選擇顏色。
圖14
依次點擊報告-場輸出，在彈框中位置下拉框中選擇唯一節點的，在新窗口選擇RF下拉菜單中的RF1與RF2；選擇U下拉菜單中的U1與U2。
展開 接觸分析節點位移過大
從hypermesh里面前處理好之后，導入ansys（采用mpc184加轉矩），每次都會提示總有節點位移過大。我把相應節點選出來，相關單元選出來一看，都是面單元，是怎么回事啊？我想著把這些單元刪除，刪了很久，沒刪多少，求大神指導……很糾結。 發現那些面單元都在接觸面上。
基于節點位移的應力強度因子外推法
今天木木給大家分享的是基于節點位移求解應力強度因子，相比于上一期出的基于單元應力求解應力強度因子得出的結果更加接近解析解。這一期包括以下內容：（1）簡要講述INP文件（2）運用最小二乘法進行線性擬合（3）對裂尖數據進行特殊處理。
Abaqus圓形激光溫度-位移耦合案例教學 ￥19.98
7、 附件：本案例中的abaqus模型文件（包括cae和激光子程序）
批量提取Abaqus的節點坐標(初始坐標、指定Step下的變形量、變形后節點坐標) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發來批量提取ABAQUS輸出數據庫（ODB）文件中指定Step下的Set節點集變形量。通過詳細的步驟說明、代碼示例和圖片展示，您將學會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個增量下的變形后的節點坐標的話，在提取變形量的基礎上，與初始坐標進行簡單的計算就可以求得坐標。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領域，提取ABAQUS輸出數據庫（ODB）文件中的節點集變形量是一項常見任務。然而，手動提取這些數據是一項繁瑣且容易出錯的工作。因此，需要一種自動化的方法來批量提取指定步驟下按節點集組織的變形量數據。</p><h2>2. 實例展示</h2><p>假設我們有一個名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節點集。運行以上代碼后，腳本會自動將該步驟下節點集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 PhiPsi-高性能有限元程序
最近，PhiPsi結合了PyQt5，形成了一個軟件的界面形式PPView：
軟件的使用，相比于前期版本的使用方法要方便很多，只需鼠標點點點，即可完成數據讀取、內核計算、后處理自定義顯示，作為有限元寶藏up主，怎么可能不整點花活兒呢~
模型測試
在體驗完XFEM的一些功能后，我就又回到了有限元求解上，網上都說Fortran計算效率高，我倒要看看到底有多高，拿一個164萬個C3D8單元有限元模型作測試，對比一下Abaqus的計算速度與PhiPsi的計算速度！
底部約束所有的自由度，頂部施加向上的壓力，簡單的拉伸模型：
求解設置
Abaqus采用迭代解6線程并行計算
PhiPsi采共享內存并行計算技術：OpenMP，6線程并行計算，在剛度矩陣的求解上采用了EBE-PCG預處理迭代求解技術，相關理論可參考《有限元方法編程（第五版）》P55。
求解時間對比
場變量均只輸出節點位移，PhiPsi的結果文件按照二進制格式輸出。
Abaqus共使用5分42秒
PhiPsi共使用298秒，也就是4分58秒，比Abaqus快了44秒
求解精度對比
Abaqus的節點位移最大為:
0.0162334625623317
PhiPsi節點位移最大為:
0.016423374759488234
吻合度98.83%，精度可以接受。
番外篇
對于查看Abaqus模型中場邊量的最大值，可以在cae中進行查看：
小數點后面最大可以保留九位，對于上文提及的0.0162334625623317是怎么得到的呢？
有限元UP主又可以開始整活兒了！直接來一個python腳本，自定義程度極高，想查看什么就查看什么，相對數據怎么操作就怎么操作，想轉換什么數據格式就轉換...好了進入正題！
展開 ABAQUS批量提取部件節點集節點編號及坐標腳本 ￥30
基于python的ABAQUS批量提取部件節點集節點編號及坐標二次開發腳本
abques鋼管混凝土柱-型鋼梁節點，在梁端用位移加載，柱頂施加軸力，為什么軸力變大就不收斂呢？怎么解決？
abques鋼管混凝土柱-型鋼梁節點，在梁端用位移加載，柱頂施加軸力，為什么軸力變大就不收斂呢？怎么解決？
單元積分點應力如何外插至節點上 | 數值實現篇
繼上次的推文：有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】，本次分享單元積分點應力外插至節點處的數值實現過程。
數值實現
借助以上理論，我們可以基于matlab平臺編制以下代碼段：
% 將積分點應力外插至單元節點上，這里只列舉了Q4的情況
for i = 1:3
StressElem(e,:,i) = [1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3) -0.5;
-0.5 1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3);
1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3) -0.5;
-0.5 1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3)]*...
[stress(e,1,i);stress(e,2,i);stress(e,3,i);stress(e,4,i)];
end
對標Abaqus
模型材料參數為普通的線彈性材料，單元類型選擇CPS4，網格劃分及邊界條件設置如下：
在結果對標過程中，可以先對比自研程序與Abaqus的節點位移場：
Abaqus位移場結果
自研程序位移場結果
在位移場一致的前提下，我們再來對標應力結果。以常見的mises應力為例：
Abaqus位移應力場結果
自研程序應力場結果
結果是一致的，說明了程序的正確性。
展開 適用于ABAQUS的粘彈性人工邊界及其等效節點力的施加程序 ￥150
程序適用于二維多土層粘彈性邊界和地震波等效節點力的加載；可以實現P波和SV波的斜入射。程序用MATLAB編寫
注意:本程序用MATLAB編寫；本程序僅限于模型網格是規則的，請參考圖片；由于本物品并非實體，因此賣出概不退換，因此購買前請詢問清楚。
編輯
基于hyperworks/abaqus位移加載-02 ￥12
本案例是基于hyperworks/abaqus簡單的模擬位移加載分析，重點在于說明如何在hyperworks中完成前處理（部件建立、網格劃分、材料創建、屬性定義、位移加載設置、約束設置、接觸設置、分析步設置等），接著導出inp模型文件并在abaqus中進行求解計算，abaqus只是扮演一個求解器的角色，hyperview中進行后處理。
本案例模型文件前處理全部在hyperworks中完成，要查看前處理如何設置，只需要在hyperworks的abaqus操作界面，導入inp模型便可查看。凡購買本案例的朋友在操作上有什么疑問，可以私信我！如果你只在hyperworks中完成部件建立、網格劃分、材料創建、屬性定義，連接關系的創建，然后在abaqus中完成加載、約束、接觸等設置并提交計算的話，遇到一些常見的問題可以關注我之前發的帖子《Hyperworks其它模塊轉到ABAQUS模塊中常會遇到的問題及解決方法匯總》。
展開 
通過Abaqus python腳本批量獲取節點的應力 ￥25
背景
有限單元法計算單元積分點的應力應變值，而對于節點的應力應變值是通過外插得到的，Abaqus中云圖顯示的就是經過插值和平均后的節點的值。通過工具欄的Query-Probe values可以查看單元或節點的應力應變等結果。
對于自動化的后處理場景，通常需要自動批量地獲取單元/節點的結果，通常都需要通過python腳本來實現。通過類似odb.steps['Step-1'].frames[-1].fieldOutputs['S']的場輸出可以比較方便地直接獲得單元的積分點應力，但沒有直接的API可以獲取節點的應力應變等結果。
如果需要獲取部件表面節點應力，可以通過創建路徑+XYData的方式實現，但想要獲得最大節點應力，則該方式不便實現。
2. 通過python腳本獲取節點應力結果
本文通過fieldOutput.getSub()函數獲取所有單元的節點結果，并對每一節點關聯的多個單元的節點值進行平均后得到節點的結果。以下以某個簡單的odb結果進行驗證。
（1）批量獲得節點的mises應力值
（2）批量獲得節點的X方向正應力值
（3）批量獲得節點的最大主應力值
（4）獲取節點的最大mises應力及編號
3. 獲取節點應變等結果
只需將腳本程序中的應力場改為應變成E等即可，此處不再演示。
以下為本文的python腳本代碼（代碼中作了必要的簡單注釋）。
展開 如何在abaqus中得到荷載-位移的數據
一般需要一個參考點（就是想得到某處的曲線，就在這定義個參考點），在step設置輸出變量field out 時，單獨對這個參考點輸出位移和反力兩個變量
1.在后處理時（visualization模塊下） 有一個按鈕（上邊是XY下面幾行是空白 鼠標放上去會顯示Create XY Data）點擊
2. 在彈出的對話框中選第四個 operate on XY data 然后 continue
3. 在彈出的操作框中最底下一行 頭一個按鈕 create XY data ，在彈出的對話框中選第二個odb field output然后continue
4. 在variables選項卡中的position下拉框里選擇unique nodal 在下面的變量里勾選RF或RT（反力）、U（位移）一般只選某個方向的（如2方向）；在elements/nodes選項卡中的method選擇Node sets，右邊選擇你定義的參考點 點擊Save
5.這時在操作框里XY Data欄下會有兩個數據，他們是參考點處的反力和位移隨時間的變化，在右邊的operators里有一個函數combine(x,x),點一下這個函數會出現在expression欄里，將兩個數據位移和反力用add to expression添加到combine函數的括號里，注意位移在前，反力在后，中間的逗號是英文的“,”
6.將expression另存為（save as按鈕）一個新的名字，可以用plot expression查看曲線，也可以在主窗口的XY Data manager用plot查看，用edit讀取數值
如果覺得位移和反力的符號是相反的，可以在第5步combine之前將兩個數據反號另存為新的數據之后combine
展開 基于hyperworks/abaqus位移加載-01 ￥5
本案例是基于hyperworks/abaqus簡單的模擬位移加載，重點在于說明如何在hyperworks中完成前處理（部件建立、網格劃分、材料創建、屬性定義、約束設置、位移加載設置、分析步設置等），接著導出inp模型文件并在abaqus中進行求解計算，abaqus只是扮演一個求解器的角色，hyperview中進行后處理。
本案例inp模型文件前處理全部在hyperworks中完成，要查看前處理如何設置，只需要在hyperworks的abaqus操作界面，導入inp模型便可查看。如果你只在hyperworks中完成部件建立、網格劃分、材料創建、屬性定義，連接關系的創建，然后在abaqus中完成加載、約束、接觸等設置并提交計算的話，遇到一些常見的問題可以關注我之前發的帖子《Hyperworks其它模塊轉到ABAQUS模塊中常會遇到的問題及解決方法匯總》。
展開 abaqus拉伸后處理（應力-應變，位移-力的輸出）
abaqus拉伸后處理（應力-應變，位移-力的輸出）
lashen.zip
