abaqus提取彎矩值的案例
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實(shí)用技巧：彎矩、軸力、剪力提取方法詳解
1
前言
最近在一個(gè)交流群里有人問(wèn)到Abaqus中彎矩、軸力、剪力等內(nèi)力的提取方法，我就好奇的查閱資料并實(shí)踐了一下。基本上通過(guò)Abaqus的gui界面操作可以獲取任意截面在任意時(shí)刻的內(nèi)力，總的來(lái)看方法都是一樣的，小差異在于軟件中的某些選項(xiàng)的設(shè)置。在使用Abaqus進(jìn)行結(jié)構(gòu)分析時(shí)，提取彎矩、軸力和剪力是非常常見(jiàn)的需求。下面我將詳細(xì)介紹一些在Abaqus中提取這些結(jié)果的實(shí)用技巧。
2
操作步驟
1.進(jìn)入后處理，點(diǎn)擊view cut,如下圖
2.選擇截面位置和截面數(shù)量，具體見(jiàn)下圖
3.進(jìn)入report，點(diǎn)擊free body cut,然后進(jìn)行下圖的設(shè)置
4.把導(dǎo)出的數(shù)據(jù)復(fù)制到excel中即可，分析和繪圖等
3
結(jié)語(yǔ)
最近在學(xué)習(xí)和對(duì)比一些界面本構(gòu)的差異和效果，發(fā)現(xiàn)網(wǎng)上一些cohesive單元插入插件使用起來(lái)不僅方便，而且可以實(shí)現(xiàn)一些特殊的效果，如果大家有這方面的插件希望分享下，我也是剛開(kāi)始接觸，這些對(duì)我也很重要，后續(xù)我也會(huì)繼續(xù)分享使用所得的技巧和收集的相關(guān)資料。
文章來(lái)源：力學(xué)混子愛(ài)AI
展開(kāi) 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寫小程序進(jìn)行后處理的腳本處理(abaqus)，近期嘗試用python寫了個(gè)讀取文件/判斷最值/寫出文件的一個(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) 