abaqus提取應(yīng)力值的案例
ABAQUS中mises應(yīng)力云圖顯示的最大值還不到屈服應(yīng)力值為啥還有PEEQ值
ABAQUS中mises應(yīng)力云圖顯示的最大值還不到屈服應(yīng)力值為啥還有PEEQ值，PEEQ云圖有變形值
ANSYS瞬態(tài)分析全時(shí)程結(jié)構(gòu)響應(yīng)最大值的提取方法（變形、應(yīng)力、應(yīng)變、能量） ￥100
</p><p>同樣的方法，可以提取全時(shí)程最大的位移、應(yīng)力、應(yīng)變、能量等結(jié)果。</p><p><img src="https://img.jishulink.com/upload/202302/623025b5c0d646b9973cd2adc6c6037f.png" alt="1.png"></p><p>收費(fèi)內(nèi)容為相關(guān)命令流。</p>
Abaqus中的特征值提取
在ABAQUS中tigong了三種提取特征值的求解器：AMS (Automatic Multi-level Substructuring) 求解器、Lanczos特征值求解器和子空間迭代求解器。Lanczos特征值求解器是默認(rèn)的方法。
AMS特征值求解器是一種高效的，針對(duì)大規(guī)模問(wèn)題的能提取大量特征值的方法，主要適用于1百萬(wàn)自由度以上的模型及500階模態(tài)以上。
它包含3個(gè)求解步驟：
（1）生成子結(jié)構(gòu)；
（2）獲得特征值；
（3）從縮減的向量中獲得全部特征向量。在Abaqus中采用AMS特征值求解器的句法（提取100階頻率）
*STEP
*FREQUENCY,ELGENSOLVER=AMS
,,100
對(duì)于傳統(tǒng)的結(jié)構(gòu)，Lanczos是默認(rèn)的特征值提取方法。然而對(duì)于大規(guī)模問(wèn)題，相比AMS方法，Lanczos方法效率較低。使用Lanczos方法，可以指定最大頻率的提取或提取的頻率數(shù)量，也可以指定最小頻率的提取。Lanczos特征值求解器允許計(jì)算到特征值真正的誤差限制時(shí)才終止，可以滿足正常的終止原則。對(duì)于多數(shù)問(wèn)題，相對(duì)誤差為1.E-12數(shù)量級(jí)，因此Lanczos求解器的計(jì)算結(jié)果精度一般要比子空間迭代法高。而子空間迭代法的終止條件是通過(guò)判斷從這一次迭代到下一次迭代過(guò)程中特征值的相對(duì)變化來(lái)實(shí)現(xiàn)的，如果相對(duì)變化小于1.E-5則認(rèn)為已經(jīng)收斂，結(jié)束計(jì)算。
abaqus中的特征值提取.pdf
展開(kāi) ABAQUS 提取接觸對(duì)主面節(jié)點(diǎn)滑移值出現(xiàn)的問(wèn)題
the selected primary variable is not available in the current frame for any elements in the current display group
場(chǎng)輸出選擇了整個(gè)模型的CDISP,就是兩個(gè)簡(jiǎn)單的方塊受力接觸，結(jié)構(gòu)如封面所示，單獨(dú)查看接觸對(duì)主面所在的單元滑移值出現(xiàn)這個(gè)錯(cuò)誤，無(wú)法顯示云圖，單獨(dú)查看應(yīng)力值，變形值都可以，唯獨(dú)滑移值不可以，單獨(dú)查看從面單元的滑移值也可以，只有主面的滑移無(wú)法單獨(dú)查看，提取主面節(jié)點(diǎn)的滑移，會(huì)說(shuō)沒(méi)有相對(duì)應(yīng)的xy數(shù)據(jù)
展開(kāi) 
基于python 提取整個(gè)模型 節(jié)點(diǎn)力最大最小值的腳本(abaqus)
以前用fortran寫(xiě)小程序進(jìn)行后處理的腳本處理(abaqus)，近期嘗試用python寫(xiě)了個(gè)讀取文件/判斷最值/寫(xiě)出文件的一個(gè)腳本
也許對(duì)一些朋友有借鑒作用
import re
Data = []
f = open('data.txt')
line = f.readline()
Num_time=0
Num_nd=0
Num_line=1
Num=0
time=[]
column=1
while line:
searchObj= re.search( r'n o d a l f o r c e g r o u p o u t p u t t=', line)
searchObj2= re.search( r'nd#', line)
if searchObj:
Data.append([])
Num_time+=1
Num_nd=0
s=line
parts= (re.split( '\s*=\s*' ,s ))
TT=float(parts[1])
time.append(TT)
Data[Num_time-1].append([])
Data[Num_time-1].append([])
Data[Num_time-1].append([])
if searchObj2:
Num_nd+=1
# print (Num_nd)
s=line
parts= (re.split( '\s*' ,s ))
Num_node
展開(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實(shí)用子程序SPRINC提取主應(yīng)力
調(diào)用的時(shí)候直接在UMAT子程序中輸入：CALL SPRINC(S,PS,LSTR,NDI,NSHR)
需要提供給該實(shí)用子程序的變量包括：
實(shí)用子程序返回UMAT的變量
實(shí)例：
CALL SPRINC(S,PS,LSTR,NDI,NSHR) 可以直接寫(xiě)成
CALL SPRINC(STRESS,PS,1,3,3)或CALL SPRINC(STRESS,PS,1, NDI,NSHR)
這行代碼的意思是將UMAT中的應(yīng)力分量數(shù)組STRESS提供給實(shí)用子程序SPRINC計(jì)算；1表示LSTR=1，表明提供給實(shí)用子程序SPRINC的為應(yīng)力分量數(shù)組；后面的兩個(gè)“3”分別表示正應(yīng)力分量個(gè)數(shù)為3和切應(yīng)力分量個(gè)數(shù)為3，也可以不用改變直接使用NDI和NSHR，因?yàn)樵谌S模型的UMAT中NDI和NSHR已經(jīng)被賦予了值為3。PS數(shù)組儲(chǔ)存返回的計(jì)算得到的主應(yīng)力，這個(gè)需要自己在umat程序中定義數(shù)組PS.
下面是關(guān)于返回主應(yīng)力順序的驗(yàn)證：
CALL SPRINC(STRESS,PS,1,3,3)
STATEV(4)=PS(1)
STATEV(5)=PS(2)
STATEV(6)=PS(3)
需要注意以下點(diǎn)：
1.PS數(shù)組一般需要在程序中自定義數(shù)組變量；
2.S數(shù)組一般也需要自定義數(shù)組變量，且需要為S數(shù)組賦予應(yīng)力分量或應(yīng)變分量的值；
3.如果SPRINC或SPRIND返回的主應(yīng)力為零，檢查輸入的變量或數(shù)組是否未定義或者未賦值，其次檢查輸出數(shù)組是否未定義。
展開(kāi) Abaqus插件——平均應(yīng)力應(yīng)變提取 ￥60
通過(guò)該插件可實(shí)現(xiàn):
1）提取所有幀的任意單元集合的平均應(yīng)力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區(qū)域的x、y、z方向的平均應(yīng)變（事先定義節(jié)點(diǎn)集合，如圖中的SET-2）
3）將以上數(shù)據(jù)保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項(xiàng)******************************
1、插件使用過(guò)程中，如有任何問(wèn)題請(qǐng)發(fā)郵件至shenz1hao@126.com
2、插件僅做學(xué)習(xí)交流使用，尊重原創(chuàng)者，切勿以營(yíng)利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車(chē)
2、將Stress-Strain文件夾解壓至當(dāng)前目錄下
3、打開(kāi)abaqus，菜單欄中點(diǎn)擊plug-ins，里面找出Stress-Strain
4、輸入相應(yīng)參數(shù)（hx、hy、hz表示x、y、z方向模型長(zhǎng)度，當(dāng)以上三參數(shù)取1時(shí)輸出的為該方向位移）
*****************************************************************
展開(kāi) 在abaqus中采用python提取結(jié)果-平均應(yīng)力 ￥5
Python腳本操作ODB文件進(jìn)行計(jì)算結(jié)果數(shù)據(jù)的提取。
因?qū)Y(jié)果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應(yīng)力，目前網(wǎng)上雖有很多代碼是關(guān)于單元集的平均應(yīng)力的提取，但并未有針對(duì)每一子步都需要提取結(jié)果的代碼。故針對(duì)此需要編寫(xiě)了python代碼。
ABAQUS提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個(gè)odb）生成csv（excel）文件。提供源文件，注釋詳細(xì)，可根據(jù)需要進(jìn)行修改。
怎么在abaqus提取xfem下，動(dòng)態(tài)三維應(yīng)力強(qiáng)度因子？
求解答
Matlab“稍作修改”Abaqus-odb結(jié)果 | 提取三維裝配體Mises應(yīng)力
數(shù)據(jù)及函數(shù)文件獲取方式：
Matlab“稍作修改”Abaqus-odb結(jié)果 | 提取三維裝配體Mises應(yīng)力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
ABAQUS怎么提取扣件支反力和結(jié)構(gòu)層壓應(yīng)力
ABAQUS怎么提取扣件支反力？軌道結(jié)構(gòu)各結(jié)構(gòu)層壓應(yīng)力是提取S22還是最大主拉應(yīng)力
Abaqus CAE Python后處理提取每一幀最大等效應(yīng)力
使用Python語(yǔ)言對(duì)Abaqus CAE后處理結(jié)果進(jìn)行分析，并提取一個(gè)分析步中每一幀的最大等效應(yīng)力，其中Python代碼如下：
from abaqus import *
from abaqusConstants import *
from odbAccess import *
import visualization
myFile=open('DATA.txt','w')
print('********************************\n')
myFile.write('********************************\n')
myOdb=openOdb(path='viewer_tutorial.odb')
myStepValue=myOdb.steps.values()
for step in myStepValue:
print('The current step is: %s.\n'%step.name)
myFile.write('The current step is: %s.\n'%step.name)
frameID=0
for frame in step.frames:
print('The current frame is: %d.\n'%frameID)
myFile.write('The current frame is: %d.
展開(kāi) 