abaqus提取頻率的案例
針對(duì)ABAQUS掃頻odb結(jié)果各頻率下最大位移快速提取Python程序 ￥2
本帖是針對(duì)ABAQUS掃頻仿真項(xiàng)目中遇到的最大值提取需求而產(chǎn)生的具體應(yīng)用。一般掃頻結(jié)束后有對(duì)各頻率下最大位移結(jié)果進(jìn)行提取并繪制曲線的需求，通常手動(dòng)提取僅可用于較少頻率提取情況，當(dāng)頻率點(diǎn)較多（如500時(shí)）手動(dòng)提取將是災(zāi)難性的操作方法。
這里利用python程序?qū)哳lodb最大值進(jìn)行提取。方法分兩類：遍歷節(jié)點(diǎn)法和Visualization顯示值提取法。前者在《python語言在ABAQUS中的應(yīng)用》一書中有節(jié)點(diǎn)應(yīng)力提取案例描述，優(yōu)點(diǎn)是不需對(duì)ABAQUS界面進(jìn)行python操作，可定位具體節(jié)點(diǎn)信息，缺點(diǎn)是速度慢；后者相反。
后一方法的應(yīng)用也可應(yīng)用到最大Mises應(yīng)力等結(jié)果數(shù)據(jù)的快速提取方面。
如有疏漏，煩請(qǐng)指教。
展開 〔abaqus〕頻率分析與復(fù)數(shù)頻率分析
頻率分析: natural frequency extraction 只能分析對(duì)稱的剛度矩陣和質(zhì)量矩陣， 如果涉及到非對(duì)稱的剛度矩陣，質(zhì)量矩陣，阻尼矩陣。 則必須要使用復(fù)數(shù)頻率分析。 complex eigenvalue extraction 在進(jìn)行復(fù)數(shù)頻率分析之前，必須要先進(jìn)行模態(tài)分析。
ABAQUS螺栓滑移量批量提取工具，可以同時(shí)對(duì)多個(gè)工況多個(gè)面進(jìn)行提取 ￥100
適用所有螺栓滑移量提取，腳本與使用方法付費(fèi)
螺栓滑移量提取1.pptx
ABAQUS彈簧質(zhì)量系統(tǒng)固有頻率求解
今天跟大家聊一聊我們?cè)诮Y(jié)構(gòu)力學(xué)與結(jié)構(gòu)動(dòng)力學(xué)里面常見的一個(gè)計(jì)算公式——彈簧質(zhì)量系統(tǒng)的固有頻率求解：
學(xué)過結(jié)構(gòu)力學(xué)或者結(jié)構(gòu)動(dòng)力學(xué)的同學(xué)都知道我們系統(tǒng)的固有頻率求解，求解公式如下：
式中的f0即為固有頻率，k為系統(tǒng)的剛度（N/m），m為系統(tǒng)質(zhì)量（kg）。
假定我們的模型如下所示：
那么由上我們可以計(jì)算出一個(gè)彈簧質(zhì)量系統(tǒng)的固有頻率，如果我們的k=400N/m，m=10kg，那么通過上式可以計(jì)算得到我們的系統(tǒng)固有頻率為1.00658。由此建立我們的ABAQUS有限元模型如下：
1.建立一個(gè)點(diǎn)部件，坐標(biāo)輸入（0，0，0）
2.鼠標(biāo)左鍵長(zhǎng)按1處圖標(biāo)選擇通過偏移形成參考點(diǎn)，通過參考點(diǎn)RP偏移1000mm生成3處參考點(diǎn)
3.導(dǎo)入點(diǎn)部件進(jìn)行裝配
4.在分析步模塊建立線性攝動(dòng)求解類型，頻率求解分析步
5.采用Lanczos求解，頻率求解值設(shè)為1即可
6.在相互作用模塊對(duì)基準(zhǔn)點(diǎn)建立參考點(diǎn)1，即RP-1
7.在上欄special中的彈簧模塊建立兩點(diǎn)之間的彈簧
8.設(shè)置彈簧剛度，在ABAQUS的mm制單位中剛度設(shè)置為0.4N/mm
9.在上欄special慣性與質(zhì)量中設(shè)置RP-1的質(zhì)量為0.01t
10.設(shè)置兩點(diǎn)的邊界條件，其中RP點(diǎn)6個(gè)自由度完全限制，RP-1點(diǎn)除圖中x方向自由度（即U1）其余自由度完全限制
11.無網(wǎng)格劃分操作，設(shè)置job，求解job得到結(jié)果
由上得到我們的結(jié)果，頻率為1.0066，與我們通過公式計(jì)算所得到的1.00658相差無幾，誤差很小。
以上就是我們今天關(guān)于彈簧質(zhì)量系統(tǒng)的固有頻率求解的討論，謝謝大家！我是食詩吃詞！SSCC!
展開 
UGNX NASTRAN、ALGOR、ABAQUS求解頻率的比較
為了便于比較，對(duì)一個(gè)簡(jiǎn)單的長(zhǎng)立方體進(jìn)行模態(tài)頻率分析。
采用同樣的參數(shù)，差別不大。UG和ALGOR的差別特別小。
模型stp格式：
blockstp.rar
密度：7.829e-09T/mm^3;
楊氏模量：206940MPa；
泊松比：0.288
模型大小：1000*100*100
一端完全固定，網(wǎng)格大小為10。
求前15階模態(tài)頻率。
Abaqus線性動(dòng)力學(xué) – 結(jié)構(gòu)受力後共振頻率改變
前言：
先前的帖子有介紹藉由提取結(jié)構(gòu)共振頻率來探討其動(dòng)態(tài)特性，本次的主題也是跟共振頻率相關(guān)，主要介紹結(jié)構(gòu)受力后共振頻率的改變。回到最基本的關(guān)係式，從單自由度系統(tǒng)我們知道，結(jié)構(gòu)共振頻率不外乎與結(jié)構(gòu)的質(zhì)量及剛度矩陣有關(guān)，只要上述兩個(gè)項(xiàng)目有改變，就會(huì)直接影響結(jié)構(gòu)共振頻率，最常見的就是考慮幾何非線性的情況下，受力后因?yàn)閯偠染仃嚫淖兌绊懝舱?em>頻率。
模型說明：
本文欲探討電路板結(jié)構(gòu)受力后共振頻率之改變，模型如下所示。
分析步設(shè)定：
建立分析步時(shí)，在Frequency前面加上Static, General靜態(tài)分析步，必須在求解共振頻率之前讓結(jié)構(gòu)受力，才能觀察其共振頻率的改變，記得要勾選幾何非線性Nlgeom才能考慮大變形！
邊界條件：
提取共振頻率的分析步將電路板底部bracket完全固定，而在靜力分析不時(shí)，設(shè)定右邊bracket向右位移0.125單位。預(yù)期結(jié)構(gòu)受拉力后共振頻率提高，這個(gè)概念可以聯(lián)想到吉他的弦，當(dāng)弦拉得越緊，來回綁盪的頻率越高，相反則是屈曲，細(xì)長(zhǎng)構(gòu)件受壓到超過容許范圍將發(fā)生屈曲。因此，可以想像結(jié)構(gòu)受拉力剛度提高、受壓力剛度降低。
結(jié)果檢核：
檢查最容易被激發(fā)動(dòng)態(tài)效應(yīng)的第一模態(tài)，結(jié)構(gòu)共振頻率從102.15Hz提高到202.69Hz，符合前面的預(yù)期，結(jié)構(gòu)受到拉力后，因?yàn)閯偠忍岣咴斐晒舱?em>頻率也提高。
后處理小技巧：
同時(shí)有Static, General及Frequency分析步時(shí)，撥放動(dòng)畫會(huì)接連著撥，可以從后處理取消顯示特定的分析步，避免動(dòng)態(tài)顯示跳來跳去。
更多線性動(dòng)力相關(guān)教程請(qǐng)參考：
https://www.yqgqt.org.cn/college/video/c15346
展開 ABAQUS提取面積或體積
也是在站內(nèi)找到別的大佬的評(píng)論看到的，分享給大家。直接使用查詢工具就可以，選擇“質(zhì)量屬性（mass properties）”
點(diǎn)擊后選擇幾何區(qū)域
在上邊選擇“面”
然后直接選擇需要得到面積的面點(diǎn)擊“完成”就可以啦！體積也是同樣的道理，在這就不演示啦！
附圖是站內(nèi)評(píng)論區(qū)看到的，再次感謝大佬！
ABAQUS實(shí)用子程序SPRINC提取主應(yīng)力
實(shí)用子程序SPRINC
在ABAQUS中使用UMAT子程序時(shí)有時(shí)會(huì)使用到最大主應(yīng)力進(jìn)行計(jì)算。通過查閱幫助文檔，ABAQUS實(shí)用子程序SPRINC可以在UMAT中計(jì)算最大主應(yīng)力和最大主應(yīng)變，SPRIND可以計(jì)算最大主應(yīng)力和最大主應(yīng)變的方向。
下面是ABAQUS幫助文檔關(guān)于實(shí)用子程序SPRINC的介紹：
SPRINC (calculate principal values)
Interface
CALL SPRINC(S,PS,LSTR,NDI,NSHR)
Variables to be provided to the utility routine
S
Stress or strain tensor.
LSTR
An identifier. LSTR=1 indicates that S contains stresses; LSTR=2 indicates that S contains strains.
NDI
Number of direct components.
NSHR
Number of shear components.
Variables returned from the utility routine
PS(I), I=1,2,3
The three principal values.
展開 ABAQUS后處理之提取分層損傷面積/分層面積/基體損傷面積（ABAQUS+Photoshop） ￥28
ABAQUS后處理之提取損傷面積（ABAQUS+Photoshop聯(lián)合使用）
為了定量描述損傷程度，提取載荷造成的損傷面積變得尤為重要，下面介紹損傷面積的提取方法。
1. 去除單元網(wǎng)格，以及邊緣
2. 突出顯示損傷區(qū)域，建立損傷與未損區(qū)域色差
3. 導(dǎo)出圖片
Abaqus Python 二次開發(fā)-歷程輸出數(shù)據(jù)提取的python實(shí)現(xiàn) ￥3.49
在Abaqus中，歷程變量的輸出，尤其是反力及位移，經(jīng)常是如下的格式。
此時(shí)，如果想提取數(shù)據(jù)，一些關(guān)于Abaqus 二次開發(fā)的書籍里提供的方式是：odb.step['Step-1'].historyRegions['Node PART-TRIMMED-MESH.288422']. historyOutputs['U1'].data 。這種方式，需要輸入 'Node PART-TRIMMED-MESH.288422' 這種很長(zhǎng)的字符串，繁瑣且容易出錯(cuò)。
在付費(fèi)內(nèi)容中，提供了一種便捷的方法，得到長(zhǎng)字符串，提取相關(guān)的數(shù)據(jù)，并將數(shù)據(jù)保存為csv文件。
展開 利用Python提取ABAQUS的.odb結(jié)果中的輪軌接觸信息 ￥500
ABAQUS 顯式explicit計(jì)算時(shí)所有結(jié)果都將寫入到.odb文件中，在GUI界面進(jìn)行數(shù)據(jù)提取時(shí)效率較低，且有些詳細(xì)結(jié)果不能直接從.odb中獲取。因此，需要利用python讀取.odb，然后進(jìn)行提取關(guān)心的數(shù)據(jù)，以此進(jìn)行下一步計(jì)算。
ABAQUS實(shí)體、殼、梁?jiǎn)卧妮S力、剪力、彎矩的提取方式及準(zhǔn)確性驗(yàn)證 ￥8
在ABAQUS中，對(duì)結(jié)構(gòu)或者構(gòu)件進(jìn)行受力分析除了分析應(yīng)力云圖之外，通常還需要對(duì)部件的軸力、剪力或彎矩的變化趨勢(shì)進(jìn)行分析。本帖基于以下的實(shí)體solid、殼shell、梁/beam（truss）模型，分別提取這三類模型的軸力、剪力、彎矩，并與理論計(jì)算相結(jié)合，驗(yàn)證提取結(jié)果的準(zhǔn)確性，并解釋相應(yīng)有限元的計(jì)算原理。
計(jì)算模型
梁?jiǎn)卧?jì)算結(jié)果
實(shí)體單元計(jì)算結(jié)果
殼單元計(jì)算結(jié)果
帖子內(nèi)容概況
批量提取Abaqus指定Step-x下 Set節(jié)點(diǎn)集的反力RF（Reaction force） ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step-x下的Set節(jié)點(diǎn)集的反力RF（Reaction force）。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會(huì)如何使用該腳本，自動(dòng)化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，RF1（X），<span style="color: rgb(25, 27, 31);">RF2（Y），RF3（Z）</span>）。</p><p><br></p><p>如果還需要提取<strong><em><u>節(jié)點(diǎn)變形量 </u></em></strong>或者<strong><em><u>節(jié)點(diǎn)坐標(biāo) </u></em></strong>可以關(guān)注：<a href="https://www.yqgqt.org.cn/post/1939915" rel="noopener noreferrer" target="_blank">利用Python腳本 批量提取Abaqus的odb文件中Set節(jié)點(diǎn)集 初始節(jié)點(diǎn)坐標(biāo)、指定Step下的變形量、變形后節(jié)點(diǎn)坐標(biāo)_CAE Abaqus提取odb-技術(shù)鄰 (jishulink.com)</a></p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點(diǎn)集的<span style="color: rgb(25, 27, 31);">反力RF（Reaction force）</span>是一項(xiàng)常見任務(wù)。然而，手動(dòng)提取這些數(shù)據(jù)是一項(xiàng)繁瑣且容易出錯(cuò)的工作。
展開 Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對(duì)ODB文件中單元集里所有積分點(diǎn)的應(yīng)力及應(yīng)變進(jìn)行自動(dòng)提取并計(jì)算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個(gè)應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個(gè)應(yīng)變平均值
批量提取Abaqus的節(jié)點(diǎn)坐標(biāo)(初始坐標(biāo)、指定Step下的變形量、變形后節(jié)點(diǎn)坐標(biāo)) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step下的Set節(jié)點(diǎn)集變形量。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會(huì)如何使用該腳本，自動(dòng)化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個(gè)增量下的變形后的節(jié)點(diǎn)坐標(biāo)的話，在提取變形量的基礎(chǔ)上，與初始坐標(biāo)進(jìn)行簡(jiǎn)單的計(jì)算就可以求得坐標(biāo)。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點(diǎn)集變形量是一項(xiàng)常見任務(wù)。然而，手動(dòng)提取這些數(shù)據(jù)是一項(xiàng)繁瑣且容易出錯(cuò)的工作。因此，需要一種自動(dòng)化的方法來批量提取指定步驟下按節(jié)點(diǎn)集組織的變形量數(shù)據(jù)。</p><h2>2. 實(shí)例展示</h2><p>假設(shè)我們有一個(gè)名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節(jié)點(diǎn)集。運(yùn)行以上代碼后，腳本會(huì)自動(dòng)將該步驟下節(jié)點(diǎn)集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 