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