abaqus 讀取坐標的案例
Abaqus二次開發讀取變形后的節點坐標并輸出到txt文件中。讀取模態頻率到txt 文件中。 ￥10
有時候在abaqus中，我們需要知道某一個集合中的節點變形后的坐標。以此為輸入來進行一些研究。這里我用一個自編函數，將節點集合變形后的坐標寫入到一個txt文件中。格式化的寫入文件，方便用此坐標來進行計算研究。
同樣有時候，需要將計算的模態頻率值提取出來。同樣用一個函數將模態頻率提取出來，放進txt文件中，方便后續研究。
讀取的節點結果如下圖所示：
ANSYS與Abaqus球坐標系下的結果讀取
ANSYS與Abaqus球坐標系下的結果讀取
1 概述
采用ANSYS和Abaqus軟件計算的結果通常默認的結果是在總體笛卡爾坐標系下產生的結果，這對于應力或者應變等分量的分析有時候不方便，比如對于一個圓筒體，比較關心其徑向應力和環向應力，而這個結果直接讀取使不可能的，需要一定的轉換。
這就是結果坐標系轉換。
在軟件里，應力分量表示為sx,xy,xz（ANSYS），s11,s22,s22（Abaqus），當其轉換到柱坐標或者球坐標時，對應的應力分量就發生變化，sx和s11均表示徑向應力。
2 ANSYS
建立一個球體模型，如圖1，加載求解，得到其總體坐標系下的sx應力分量。
圖1
在后處理器中，將結果坐標系轉換為球坐標系，采用的命令為：RSYS。查詢ANSYS幫助文檔，如圖2：
圖2 RSYS
0,1,2分別代表笛卡爾坐標系，柱坐標系，球坐標系。
輸入命令：RSYS，2
顯式結果sx為圖3，此時的sx應力分量為徑向應力。
圖3
3 Abaqus
建立模型加載求解，得到s11應力分量如圖4.
圖4
轉換結果坐標系，Visualization模塊下選擇 Tools--Create Coordinate Aystem，按指定方法建立局部坐標系，然后選擇Result-Option，選擇Transformation標簽，User-specified，就可以看到新建立的坐標系，選擇新建的坐標系即可完成坐標轉換。
如圖5，圖6
圖5
建立球坐標系的時候根據Abaqus窗口下方的提示進行操作。
圖6
最終轉換為徑向應力的顯式結果，如圖7
圖7
展開 通過讀取Z Buffer深度值將屏幕2D坐標轉化為場景3D坐標（含程序和源碼）
主要有2種方法來實現屏幕2D坐標轉化為場景3D坐標：一種是通過揀選射線，我的《APRG Demo》就是這樣做的，有興趣的可以查看源代碼：http://blog.csdn.net/skyman_2001/archive/2005/10/06/495938.aspx；另一種是通過讀取Z Buffer深度值。相比較而言后者要簡單些。這里專門講如何用后者來實現。
1. 使用雙緩存,打開深度測試：
glClearDepth(1.0f); // 深度緩存設置
glEnable(GL_DEPTH_TEST); // 打開深度測試
glDepthFunc(GL_LEQUAL); // 設置深度測試類型
2. 繪制場景：
3.
展開 基于tcl語言實現CSV文件中的坐標讀取并創建節點 ￥35
本案例在于如何使用tcl語言實現在hyperworks中實現批量讀取節點坐標及應用該節點坐標實現批量創建節點。凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！再次說明下本案例相比于《tcl實現文件讀取及節點創建》這篇案例，本案例是另外一種方法，現將CSV中的每一列坐標值付給三個列表，再讀取三個列表的每一行創建節點。
基于TCL實現文件的讀取以及節點坐標的輸出、節點編號進行renumber ￥50
本案例在于如何使用tcl語言實現在hyperworks中實現批量讀取節點坐標及應用該節點坐標實現批量創建節點。同時對這些節點的節點編號進行renumber，以及如何讀取這些節點的坐標值并將其寫入文件。凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！
批量提取Abaqus的節點坐標(初始坐標、指定Step下的變形量、變形后節點坐標) ￥40
如果需要提取變形后的節點坐標的話，我們還要繼續進行。</p><h3>3.4提取Set組中初始坐標</h3><p>方法一：提取初始節點坐標可以通過Abaqus后處理查詢結果，在CSDN上可以查詢到，此處不再贅述。</p><p>方法二：通過Python腳本，按節點編號提取初始節點坐標。</p><p><span style="color: rgb(18, 18, 18);">代碼為附件：2 </span>Python腳本-初始節點提取</p><p><img src="https://img.jishulink.com/202405/attachment/90a3bc50ff224f1887815b8d0943ed64.jpg"></p><h2>4.數據處理</h2><p>通過數據處理，將對應值求和即可求得變形后的節點坐標。以編號1號節點為例，節點初始坐標（X=19.42577，Y=2.96978），變形量（U1=0.54588，U2=0.48286），可求得編號1號節點變形后的坐標（X+U1，Y+U2），即為（19.97165，3.45264）。</p><p>驗證坐標如下圖所示，使用Hyperview后處理，可以看出提取節點坐標與Python腳本后處理的節點坐標一致。綜上所述可以看出該腳本可滿足需求。
展開 在ABAQUS中基于圓柱坐標系設置關于坐標函數的表面力（keyword 曲面加載,圓柱坐標,面力）
例如下圖所示，受Y方向某拉力作用，各點應力狀態為：
在圓孔中心位置建立圓柱坐標系，該應力狀態在圓柱坐標系下的公式為：
在這種情況下反推物理量，需要對曲面施加基于圓柱坐標系的面力。
案例如下：在圓弧面基于圓柱坐標系施加等效于單向應力狀態的面力。
加載前先建立圓柱坐標系（注意R軸方向為0度位置，T軸方向為角度增大方向，示意圖見文后的加載圖）
具體設置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設置徑向力，按以下參數設置：
Distribution：應力分配，點擊后面的f(x)創建一個基于圓柱坐標系的表達式，Local system 要選擇圓柱坐標系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點擊選擇圖標后，依次選擇（0，0，0） （-1，0，0） ，坐標選擇建立的圓柱坐標系。
注：面力方向矢量是基于所選坐標系，（-1，0，0）就是沿圓柱坐標系下的R軸反向。
Magnitude：選擇應力大小為1。
然后在創建一個Load,設置切向力，如下圖所示，也是基于圓柱坐標系。
再創建一個Load，在整體坐標系下對兩側的平面施加Y方向的面力，大小為1，同時對后面的面施加全約束。
最后加載形式為下圖所示：
求解結果如下圖：
大部分位置應力在0.99~1.01之間，為單向應力狀態，加載方式正確。
本問題的關鍵是面力的方向問題，在選擇面力的方向矢量時，是基于所選坐標系。對于圓柱坐標系，切向力矢量為（0，-1，0）時，即力的方向只沿著theta的反方向。
展開 ABAQUS 輸出節點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節點坐標
展開 abaqus怎么提取極坐標一下的節點坐標？
如題
ABAQUS二次開發-參考點坐標自動識別，與更新坐標 ￥80
ABAQUS二次開發-參考點坐標自動識別，與更新坐標
Abaqus利用region讀取結果數據
Abaqus利用region讀取結果數據
長安CAE
Abaqus計算完成后，為輸出某部分模型的結果數據，需要用到region，并且配套使用的是getSubset函數。
示例一創建一個變量center，將一個叫做PUNCH的節點set賦值給這個變量，通過整個模型的位移場displacement得到該節點set的輸出數據，代碼如下：
center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']
centerDisplacement = displacement.getSubset(region=center)
centerValues = centerDisplacement.values
for v in centerValues:
print v.nodeLabel, v.data
其中的displacement是已經創建好的位移場數據，centerDisplacement變量得到該節點集的所有數據，centerValues變量得到該節點集的數值數據，通過循環輸出節點號和當前數據。
展開 
Abaqus利用region讀取結果數據
單元質心點
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
問題2
如何從ODB結果文件中讀取計算結果數據
第二個問題是大家經常要面對的，常規的數據提取采用ABAQUS自帶的后處理工具已經能滿足需求了，但是如果我們需要進行一些寫出文本，數據篩選，數據的二次計算處理。。。，我們就需要用程序語言的幫助了。
針對今天的問題——
提取結果并寫出到文本
，如果用自帶的工具，可以導出abaqus.rpt文件，不過這種方式并不能滿足我們較為苛刻的要求，用python程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節點的各種計算結果，或者簡要的數據處理，臨時用一用，可以采用如下方式，不必費心寫復雜的程序，不值當浪費時間，
青春易逝，頭發寶貴
。當然了，通過阿信這里給出的基本程序改一改，能省幾根頭發。
限于時間關系，這里只列出基本的應力提取程序，寫入文本比較簡單，隨意百度即可實現，不再給出。
展開 Abaqus二次開發之Python腳本讀取場輸出位移
操作腳本如下：
結果如下：
Node label: 1000
Displacement in X direction: -7.891572e-34
Displacement in Y direction: -76.45553
Displacement: 76.4555282593
注：本文中涉及的模型文件和完整腳本文件可在Up的公眾號：易木木響叮當，中回復“讀取場輸出數據”自動獲取。
Abaqus利用腳本讀取結果并創建新的輸出場
Abaqus利用腳本讀取結果并創建新的輸出場
長安CAE
Abaqus的后處理有時候采用Python腳本會比較方便，此次實例通過腳本語言讀取已經存在的Odb結果文件，并且通過程序實現不同載荷步下的結果場數據運算，以得到新的場數據，并顯示其云圖。
