abaqus讀取頻率的案例
Abaqus二次開發讀取變形后的節點坐標并輸出到txt文件中。讀取模態頻率到txt 文件中。 ￥10
有時候在abaqus中，我們需要知道某一個集合中的節點變形后的坐標。以此為輸入來進行一些研究。這里我用一個自編函數，將節點集合變形后的坐標寫入到一個txt文件中。格式化的寫入文件，方便用此坐標來進行計算研究。
同樣有時候，需要將計算的模態頻率值提取出來。同樣用一個函數將模態頻率提取出來，放進txt文件中，方便后續研究。
讀取的節點結果如下圖所示：
〔abaqus〕頻率分析與復數頻率分析
頻率分析: natural frequency extraction 只能分析對稱的剛度矩陣和質量矩陣， 如果涉及到非對稱的剛度矩陣，質量矩陣，阻尼矩陣。 則必須要使用復數頻率分析。 complex eigenvalue extraction 在進行復數頻率分析之前，必須要先進行模態分析。
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程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節點的各種計算結果，或者簡要的數據處理，臨時用一用，可以采用如下方式，不必費心寫復雜的程序，不值當浪費時間，
青春易逝，頭發寶貴
。當然了，通過阿信這里給出的基本程序改一改，能省幾根頭發。
限于時間關系，這里只列出基本的應力提取程序，寫入文本比較簡單，隨意百度即可實現，不再給出。
展開 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
展開 Abaqus利用腳本讀取結果并創建新的輸出場
Abaqus利用腳本讀取結果并創建新的輸出場
長安CAE
Abaqus的后處理有時候采用Python腳本會比較方便，此次實例通過腳本語言讀取已經存在的Odb結果文件，并且通過程序實現不同載荷步下的結果場數據運算，以得到新的場數據，并顯示其云圖。
Abaqus二次開發之Python腳本讀取場輸出位移
操作腳本如下：
結果如下：
Node label: 1000
Displacement in X direction: -7.891572e-34
Displacement in Y direction: -76.45553
Displacement: 76.4555282593
注：本文中涉及的模型文件和完整腳本文件可在Up的公眾號：易木木響叮當，中回復“讀取場輸出數據”自動獲取。
讀取ABAQUS結果文件中的單元應變能并輸出至excel文件的處理腳本 ￥1
腳本內容如題。
ABAQUS彈簧質量系統固有頻率求解
今天跟大家聊一聊我們在結構力學與結構動力學里面常見的一個計算公式——彈簧質量系統的固有頻率求解：
學過結構力學或者結構動力學的同學都知道我們系統的固有頻率求解，求解公式如下：
式中的f0即為固有頻率，k為系統的剛度（N/m），m為系統質量（kg）。
假定我們的模型如下所示：
那么由上我們可以計算出一個彈簧質量系統的固有頻率，如果我們的k=400N/m，m=10kg，那么通過上式可以計算得到我們的系統固有頻率為1.00658。由此建立我們的ABAQUS有限元模型如下：
1.建立一個點部件，坐標輸入（0，0，0）
2.鼠標左鍵長按1處圖標選擇通過偏移形成參考點，通過參考點RP偏移1000mm生成3處參考點
3.導入點部件進行裝配
4.在分析步模塊建立線性攝動求解類型，頻率求解分析步
5.采用Lanczos求解，頻率求解值設為1即可
6.在相互作用模塊對基準點建立參考點1，即RP-1
7.在上欄special中的彈簧模塊建立兩點之間的彈簧
8.設置彈簧剛度，在ABAQUS的mm制單位中剛度設置為0.4N/mm
9.在上欄special慣性與質量中設置RP-1的質量為0.01t
10.設置兩點的邊界條件，其中RP點6個自由度完全限制，RP-1點除圖中x方向自由度（即U1）其余自由度完全限制
11.無網格劃分操作，設置job，求解job得到結果
由上得到我們的結果，頻率為1.0066，與我們通過公式計算所得到的1.00658相差無幾，誤差很小。
以上就是我們今天關于彈簧質量系統的固有頻率求解的討論，謝謝大家！我是食詩吃詞！SSCC!
展開 UGNX NASTRAN、ALGOR、ABAQUS求解頻率的比較
為了便于比較，對一個簡單的長立方體進行模態頻率分析。
采用同樣的參數，差別不大。UG和ALGOR的差別特別小。
模型stp格式：
blockstp.rar
密度：7.829e-09T/mm^3;
楊氏模量：206940MPa；
泊松比：0.288
模型大小：1000*100*100
一端完全固定，網格大小為10。
求前15階模態頻率。
Abaqus線性動力學 – 結構受力後共振頻率改變
前言：
先前的帖子有介紹藉由提取結構共振頻率來探討其動態特性，本次的主題也是跟共振頻率相關，主要介紹結構受力后共振頻率的改變?；氐阶罨镜年P係式，從單自由度系統我們知道，結構共振頻率不外乎與結構的質量及剛度矩陣有關，只要上述兩個項目有改變，就會直接影響結構共振頻率，最常見的就是考慮幾何非線性的情況下，受力后因為剛度矩陣改變而影響共振頻率。
模型說明：
本文欲探討電路板結構受力后共振頻率之改變，模型如下所示。
分析步設定：
建立分析步時，在Frequency前面加上Static, General靜態分析步，必須在求解共振頻率之前讓結構受力，才能觀察其共振頻率的改變，記得要勾選幾何非線性Nlgeom才能考慮大變形！
邊界條件：
提取共振頻率的分析步將電路板底部bracket完全固定，而在靜力分析不時，設定右邊bracket向右位移0.125單位。預期結構受拉力后共振頻率提高，這個概念可以聯想到吉他的弦，當弦拉得越緊，來回綁盪的頻率越高，相反則是屈曲，細長構件受壓到超過容許范圍將發生屈曲。因此，可以想像結構受拉力剛度提高、受壓力剛度降低。
結果檢核：
檢查最容易被激發動態效應的第一模態，結構共振頻率從102.15Hz提高到202.69Hz，符合前面的預期，結構受到拉力后，因為剛度提高造成共振頻率也提高。
后處理小技巧：
同時有Static, General及Frequency分析步時，撥放動畫會接連著撥，可以從后處理取消顯示特定的分析步，避免動態顯示跳來跳去。
更多線性動力相關教程請參考：
https://www.yqgqt.org.cn/college/video/c15346
展開 針對ABAQUS掃頻odb結果各頻率下最大位移快速提取Python程序 ￥2
本帖是針對ABAQUS掃頻仿真項目中遇到的最大值提取需求而產生的具體應用。一般掃頻結束后有對各頻率下最大位移結果進行提取并繪制曲線的需求，通常手動提取僅可用于較少頻率提取情況，當頻率點較多（如500時）手動提取將是災難性的操作方法。
這里利用python程序對掃頻odb最大值進行提取。方法分兩類：遍歷節點法和Visualization顯示值提取法。前者在《python語言在ABAQUS中的應用》一書中有節點應力提取案例描述，優點是不需對ABAQUS界面進行python操作，可定位具體節點信息，缺點是速度慢；后者相反。
后一方法的應用也可應用到最大Mises應力等結果數據的快速提取方面。
如有疏漏，煩請指教。
展開 