abaqus節點應力輸出的案例
ABAQUS 輸出節點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節點坐標
展開 通過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中隱式和顯式的節點和單元的輸出變量解析
FSLIP
Field: yes History: no .fil: no
Length of contact slip path at secondary nodes during contact (FSLIPEQ) and in some
cases (see About contact pairs in Abaqus/Explicit) components of net contact slip in local
tangent directions (FSLIP1 and FSLIP2). These variables remain constant while a
secondary node is not in contact.
FSLIPR
Field: yes History: no .fil: no
Magnitude of contact slip rate at secondary nodes during contact (FSLIPR) and in some
cases (see About contact pairs in Abaqus/Explicit) components of contact slip rate in local
tangent directions (FSLIPR1 and FSLIPR2). These variables are set to zero while a
secondary node is not in contact.
BONDSTAT
Field: no History: yes .fil: no
Spot weld bond status.
展開 輸出abaqus inp如何去除instance信息，而且單元節點編號也不會有重復
解決的方法：在command line輸入 mdb.models['modelName'].setValues(noPartsInputFile=ON) 重新生成inp文件，這樣在寫出的inp文件里面就沒有instance信息，而且單元節點編號也不會有重復了。
abaqus拉伸后處理（應力-應變，位移-力的輸出）
abaqus拉伸后處理（應力-應變，位移-力的輸出）
lashen.zip
Abaqus二次開發讀取變形后的節點坐標并輸出到txt文件中。讀取模態頻率到txt 文件中。 ￥10
有時候在abaqus中，我們需要知道某一個集合中的節點變形后的坐標。以此為輸入來進行一些研究。這里我用一個自編函數，將節點集合變形后的坐標寫入到一個txt文件中。格式化的寫入文件，方便用此坐標來進行計算研究。
同樣有時候，需要將計算的模態頻率值提取出來。同樣用一個函數將模態頻率提取出來，放進txt文件中，方便后續研究。
讀取的節點結果如下圖所示：
ABAQUS Fortran基于提供的代碼 P53，完善 8 節點單元程序，并增加對應力結果的處理。
1.改寫輸入數據格式，使之能適應任意幾何（可利用節點坐標輸 入節點，利用單元-節點關系輸入單元）； 2. 計算節點應力，給出并實現至少一種應力處理方案，提供處理 前后的應力結果（可用表格和云圖表示），可與其它軟件對比； 3.提交總結報告（包括方法/方案描述、帶詳細注釋的代碼、程序框圖、算例描述、結果比較分析等）、可編譯源代碼、可執行文件、 數據文件、結果文件
program p53
!-----------------------------------------------------------------------------
! program 5.3 plane strain of an elastic solid using uniform
! 8-node quadrilateral elements numbered in the x direction
!-----------------------------------------------------------------------------
use new_library ; use geometry_lib ; implicit none
integer::nels,nxe,neq,nband,nn,nr,nip,nodof=2,nod=8,nst=3,ndof,loaded_nodes,&
i,k,iel,ndim=2
real::aa,bb,e,v,det ; character(len=15) :: element = 'quadrilateral'
!
展開 ABAQUS-如何求結構的節點位移單元應力分量和支反力
這兩個輸出的是節點位移與支反力。
圖15
圖16
在位置
下拉框中選擇積分點，在新窗口選擇S下拉菜單中的S11。這輸出的是單元應力分量。
圖17
點擊窗口的設置按鈕，更改名稱為link.rpt。
圖18
在保存目錄中用記事本打開
link
.
rpt
，即可得到相應的結果。
圖19
微信公眾號：CAE小花生
有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 ABAQUS-利用python輸出單元應力S11,S12,S13,S22,S23,S33和位移 ￥20
ABAQUS可以通過讀取rpy文件進行一些簡單的輸出圖像和數據后處理的工作，如出圖時，可以找到對應的旋轉命令和圖例編輯命令
但若要批量輸出odb內集合的應力、變形結果，還需要編寫.py文件，并利用.bat調用這個py文件。
【問題描述】：在C3D4單元的inp中已創建需要輸出結點應力結果的集合'SET-6'，且利用abaqus已計算得到了一個odb，為g.odb，需要輸出結點集合的編號和應力結果？
下面進行展示：
