abaqus應(yīng)變?cè)茍D的案例
ABAQUS SpectrumColor云圖配色插件 ￥98
插件介紹
AbyssFish SpectrumColor插件可實(shí)現(xiàn)一鍵式修改ABAQUS可視化中的云圖繪制顏色譜。
插件內(nèi)置150余種不同顏色映射模式，可實(shí)現(xiàn)多種樣式的云圖繪圖風(fēng)格。
可通過(guò)插件中的顏色預(yù)覽卡，查看需要設(shè)置的云圖顏色譜風(fēng)格及名稱(chēng)。
也可在插件內(nèi)勾選離散，繪制顏色離散分布的云圖。
如需自定義顏色映射模式，可通過(guò)修改插件文件夾內(nèi)的Color文件夾下的txt顏色表實(shí)現(xiàn)，顏色格式設(shè)置如下。
說(shuō)明提醒
該插件可在 Windows 10 和 Windows 11 系統(tǒng)上運(yùn)行，支持 Abaqus 2024、2025 及更高版本。
插件需要注冊(cè)，注冊(cè)完成后即可永久使用。該插件為單機(jī)許可銷(xiāo)售模式，購(gòu)買(mǎi)后請(qǐng)通過(guò) QQ：1135122921 或微信：AbyssFish_LJR 聯(lián)系我們以獲取許可證。
展開(kāi) abaqus子程序VUSDFLD——考慮應(yīng)變率與應(yīng)變軟化效應(yīng)的軟土模型 ￥25
<p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫(xiě)的由Einav與Randolph提出的西澳模型，用于求解軟黏土體劇烈變形后的強(qiáng)度變化，可應(yīng)用于的大變形計(jì)算。</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png" title="8$U(VZ82]O{OEMQB}[P(ZMB.png" alt="8$U(VZ82]O{OEMQB}[P(ZMB.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?
展開(kāi) abaqus云圖
step1. step2跑完了，云圖只顯示到了1.5，咋回事呀，大佬們
Abaqus后處理-云圖變量含義（部分）
結(jié)果中分量說(shuō)明：
S11、S22、S33指各軸正應(yīng)力；
S12指作用于XZ平面(與“2”，即Y軸垂直的平面)內(nèi)，沿1方向剪應(yīng)力；
S13指作用于XY平面內(nèi)，沿1方向剪應(yīng)力；
S23指作用于XY平面內(nèi)，沿2方向剪應(yīng)力。
若為柱坐標(biāo)，S12、S13、S23分別指：由徑向向環(huán)向的剪應(yīng)力、由徑向向軸向的剪應(yīng)力、由環(huán)向向軸向的剪應(yīng)力。
ABAQUS中求解某部分單元的平均應(yīng)力或平均應(yīng)變 ￥10
1、參考模型：?jiǎn)蜗蚶w維的RVE模型；
2、腳本功能：針對(duì)指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計(jì)算方法。對(duì)于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過(guò)對(duì) RVE 內(nèi)每一個(gè)單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計(jì)算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對(duì)每個(gè)單元進(jìn)行應(yīng)力（應(yīng)變）的輸出和計(jì)算。
Abaqus使用Python腳本自動(dòng)保存云圖
示例代碼見(jiàn)文末，功能室自動(dòng)化保存每一幀的Mises應(yīng)力云圖。
在Abaqus/CAE中截圖的時(shí)候有一些參數(shù)可選，在使用腳本截圖的時(shí)候這些參數(shù)在printOptions中，可以調(diào)用函數(shù)
session.printOptions.setValues()
來(lái)設(shè)置。比如要取消保存的截圖上的標(biāo)題欄，可以將參數(shù)vpDecorations置為OFF。
如何在ABAQUS對(duì)可視化后，生成的變形圖云圖，顯示等值線(xiàn)？
<p>1、點(diǎn)擊：視圖切面管理器</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png" style="text-align: center">
<img src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png" data-mobile-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png?image_process=/format,webp" data-initial-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png">
</figure>
</div><p>2、點(diǎn)擊：create創(chuàng)建</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com
展開(kāi) Abaqus后處理插件——云圖場(chǎng)變量查詢(xún)器 ￥50
基于Python對(duì)Abaqus進(jìn)行后處理的二次開(kāi)發(fā)，形成了場(chǎng)變量查詢(xún)插件，寫(xiě)論文寫(xiě)報(bào)告的實(shí)用小工具，話(huà)不多說(shuō)直接上圖。
其實(shí)Abaqus后處理模塊本身具有場(chǎng)變量查詢(xún)功能（Tools-Query-Probe values），可以滿(mǎn)足使用需求，但其顯示效果非常不友好，以節(jié)點(diǎn)（單元）編號(hào)+場(chǎng)變量值的形式顯示，理工科氣息十足。不得不說(shuō)，Abaqus在操作體驗(yàn)方面而言和Workbench還是有一些差距。
PS：常規(guī)的二次開(kāi)發(fā)通過(guò)fieldvariable.getsubset方式查詢(xún)到的場(chǎng)變量為節(jié)點(diǎn)直接插值的計(jì)算結(jié)果，即節(jié)點(diǎn)上的值是從相鄰單元插值后未經(jīng)過(guò)平均化的值，是離散的，而結(jié)果顯示云圖默認(rèn)是按75%平均化處理后的結(jié)果，更接近實(shí)際。
本插件查詢(xún)到的場(chǎng)變量為平均化之后的結(jié)果，和云圖的顯示的值一致。
展開(kāi) Abaqus_Python 后處理結(jié)點(diǎn)數(shù)據(jù)提取、場(chǎng)變量云圖輸出 ￥1.99
每個(gè)case需要輸出幾個(gè)路徑上的部分節(jié)點(diǎn)值，而且需要針對(duì)不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節(jié)點(diǎn)數(shù)據(jù)，工作量會(huì)非常大，費(fèi)時(shí)會(huì)遠(yuǎn)遠(yuǎn)超出計(jì)算所用時(shí)間；并且，同事還要求我根據(jù)case名稱(chēng)、組件名稱(chēng)（模型有多個(gè)組件組成，且網(wǎng)格是在part上畫(huà)的，所以同一個(gè)節(jié)點(diǎn)編號(hào)可能對(duì)應(yīng)多個(gè)節(jié)點(diǎn)）等保存成不同的文件。
結(jié)點(diǎn)溫度數(shù)據(jù)沒(méi)有保存，因此沒(méi)有輸出。下圖是自動(dòng)輸出的png格式云圖。
比較倉(cāng)促，代碼也比價(jià)粗糙。回頭做一個(gè)詳細(xì)的介紹。
主要代碼見(jiàn)下。
Abaqus結(jié)果云圖中無(wú)增量步顯示（ODB加載問(wèn)題）
<p>我做的是掃頻分析</p><p>掃頻分析步<strong>已添加場(chǎng)輸出（默認(rèn)是沒(méi)有場(chǎng)輸出的）</strong>，提交計(jì)算完成后，監(jiān)視器中有增量步，結(jié)果云圖中無(wú)增量步顯示（ODB加載問(wèn)題）解決方法：</p><p><br></p><ol><li>上方菜單欄File-Close ODB-全部清除</li><li>Job模塊-重新打開(kāi)計(jì)算文件的Results</li><li>查看云圖，增量步結(jié)果可查看</li></ol><p><br></p><p>總結(jié)：關(guān)閉ODB，重新打開(kāi)</p><p><br></p>
展開(kāi) abaqus系列技巧15：說(shuō)一說(shuō)后處理中云圖的變形縮放倍數(shù)
開(kāi)篇聲明：對(duì)象為剛接觸abaqus的同學(xué)，力求深入淺出，不求嚴(yán)謹(jǐn)
什么是云圖的變形縮放倍數(shù)呢？英文叫做 deformation scale factor。請(qǐng)對(duì)號(hào)入座。
當(dāng)我們打開(kāi)結(jié)果云圖是，常常遇到這樣的情況：
此視不明真相的觀(guān)眾通常的操作是，打開(kāi)手機(jī)QQ，拍個(gè)照，然后發(fā)到各種交流群：大神，看看我這個(gè)怎么了？
其實(shí)完全不需要，我們關(guān)注最后一行
具體我就不翻譯了，應(yīng)該都看的懂，縮放倍數(shù)明顯太大了。
那好，我們來(lái)改正他，具體步驟如下：
是不是風(fēng)平浪靜，一切安好啊。
不著急，我們?cè)诳纯磏onuniform是做什么的。這叫做舉一反三。它是一個(gè)自定義的變形比例，如下圖設(shè)置
我們只定義Y方向的變形比例，其他兩個(gè)方向的為0，那么只有Y方向的變形被放大，這在分方向確定變形趨勢(shì)的時(shí)候很有用！
想獲得幻想飛翔最新CAE技術(shù)文章，請(qǐng)關(guān)注幻想飛翔公眾賬號(hào)：幻想飛翔CAE。
想獲得幻想飛翔最新CAE技術(shù)文章，請(qǐng)關(guān)注幻想飛翔公眾賬號(hào)：幻想飛翔CAE。
也歡迎加入群516073058進(jìn)行交流
展開(kāi) 
有限元理論基礎(chǔ)及Abaqus內(nèi)部實(shí)現(xiàn)方式研究系列27： Abaqus內(nèi)部計(jì)算和顯示的應(yīng)變
（1）顯示應(yīng)變：Abaqus計(jì)算完畢后得到導(dǎo)入結(jié)果，在后處理中查看，應(yīng)變E11=8.528e-1,E22=-5.173e-1如下：
（2）計(jì)算應(yīng)變：Abaqus中采用UMAT子程序，利用我們的子程序調(diào)試插件DUS調(diào)試UMAT，在Visual Studio中查看dStran的值，發(fā)現(xiàn)在計(jì)算完應(yīng)變后，進(jìn)入U(xiǎn)MAT時(shí)，E11=8.528e-1,E22=-5.173e-1，調(diào)試如下：
可以發(fā)現(xiàn)殼單元Abaqus的計(jì)算應(yīng)變和顯示應(yīng)變一樣，猜測(cè)都是對(duì)數(shù)應(yīng)變。
1.5.3 iSolver的應(yīng)變
iSolver中采用自帶材料進(jìn)行計(jì)算，材料參數(shù)和UMAT的輸入完全一致。
為了計(jì)算和Abaqus完全一致，iSolver也采用對(duì)數(shù)應(yīng)變計(jì)算方式，得到的應(yīng)變顯示如下，可發(fā)現(xiàn)和Abaqus完全一致。
==總結(jié)==
由上可以看到，在實(shí)際計(jì)算中，對(duì)體單元，Abaqus和iSolver都采用變形率積分方式來(lái)計(jì)算應(yīng)變，對(duì)殼單元，Abaqus和iSolver都采用對(duì)數(shù)應(yīng)變。一般理論書(shū)都認(rèn)為Abaqus是因?yàn)閷?duì)數(shù)應(yīng)變計(jì)算復(fù)雜才采用別的應(yīng)變，但個(gè)人認(rèn)為應(yīng)該不是這個(gè)原因，因?yàn)?em>Abaqus對(duì)體單元為了顯示對(duì)數(shù)應(yīng)變，依然重新計(jì)算了一遍，說(shuō)明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒(méi)研究清楚，歡迎探討。
如果有任何其它疑問(wèn)或者項(xiàng)目合作意向，也歡迎聯(lián)系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元?jiǎng)偠染仃囇芯俊?http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質(zhì)量矩陣研究。
展開(kāi) ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
批量提交的核心是需要等待當(dāng)前計(jì)算任務(wù)的結(jié)束，上圖中若干個(gè)job的提交代碼如下：
# coding: utf-8
#微信公眾號(hào)：ABAQUS二次開(kāi)發(fā)
#作者：阿信老師CAE
#email：axin_cae@163.com
#2022.03.17
from abaqus import *
from abaqusConstants import *
jobList = mdb.jobs.keys() #獲取所有計(jì)算任務(wù)的jobName
for jobName in jobList:
mdb.jobs[jobName].submit(consistencyChecking=OFF) #提交計(jì)算
mdb.jobs[jobName].waitForCompletion() #等待計(jì)算完成
print jobName , "is completed"
不過(guò)在實(shí)際的操作過(guò)程中，我們可能會(huì)需要避開(kāi)一些job不提交，或者job窗口事先并沒(méi)有建立job，等等，總之實(shí)際問(wèn)題永遠(yuǎn)比任何教程都復(fù)雜，不過(guò)只需要靈活面對(duì)就可以了，處理起來(lái)并不難。
展開(kāi) 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提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個(gè)odb）生成csv（excel）文件。提供源文件，注釋詳細(xì)，可根據(jù)需要進(jìn)行修改。
