abaqus 能量解讀的案例
Abaqus節(jié)點(diǎn)力(NFORC)如何解讀
ABAQUS中對(duì)于節(jié)點(diǎn)力NFORC的定義是：Nodal force due to element stresses，可以理解成節(jié)點(diǎn)力是由節(jié)點(diǎn)所在的單元上的應(yīng)力按照一定規(guī)則等效到節(jié)點(diǎn)上的力。
對(duì)于下面的模型，頂部右半邊施加均布?jí)毫d荷P=20MPa，受壓面積為50*100。講道理，頂部所有節(jié)點(diǎn)的節(jié)點(diǎn)力之和應(yīng)該是等于100000N。
經(jīng)過計(jì)算，并查詢頂部節(jié)點(diǎn)的Z向節(jié)點(diǎn)力，可以看到，第一列節(jié)點(diǎn)的節(jié)點(diǎn)力數(shù)值為1666N，第二列為833N，第三列為0。
然后，通過Creat XY Data的相關(guān)操作，將頂部所有節(jié)點(diǎn)的Z向節(jié)點(diǎn)力求和，其結(jié)果為35000N（具體操作過程不作說明），與正確的數(shù)值100000N相差甚遠(yuǎn)，說明節(jié)點(diǎn)力是有問題的，那么問題出在哪兒？
根源還是在NFORC的理解上。
Abaqus默認(rèn)計(jì)算是采用了Avg:75%這個(gè)評(píng)價(jià)準(zhǔn)則，這一塊可以在網(wǎng)上查到專門的解釋。
當(dāng)不采用評(píng)價(jià)準(zhǔn)則，通過Result-Options-Computation，設(shè)置Average threshold(%)為0時(shí)，重新顯示結(jié)果，再次查詢相同節(jié)點(diǎn)的節(jié)點(diǎn)力。
相比之前，最大的區(qū)別在于，現(xiàn)在該節(jié)點(diǎn)的Z向節(jié)點(diǎn)力有4個(gè)數(shù)值，而非之前的一個(gè)。
這四個(gè)數(shù)值其實(shí)代表的是這個(gè)節(jié)點(diǎn)所屬的四個(gè)單元（節(jié)點(diǎn)屬于四個(gè)單元的公共節(jié)點(diǎn)）分別等效的節(jié)點(diǎn)力，該節(jié)點(diǎn)的實(shí)際節(jié)點(diǎn)力應(yīng)該是這四個(gè)數(shù)值相加之和。
再回到前面Avg:75%的結(jié)果，評(píng)價(jià)后僅有一個(gè)數(shù)值，將該數(shù)值乘以4才應(yīng)該是該節(jié)點(diǎn)的節(jié)點(diǎn)力。
所以：為了得到節(jié)點(diǎn)的準(zhǔn)確節(jié)點(diǎn)力，需要在Avg:75%的結(jié)果基礎(chǔ)上，依據(jù)節(jié)點(diǎn)所共有的單元數(shù)目，將結(jié)果乘以該單元數(shù)目才是準(zhǔn)確的節(jié)點(diǎn)力。（注意：有的節(jié)點(diǎn)是4個(gè)單元的公共節(jié)點(diǎn)，有的節(jié)點(diǎn)只有兩個(gè)單元，例如邊界上的節(jié)點(diǎn)，還有的節(jié)點(diǎn)沒有公共單元，例如角上的頂點(diǎn)）。
展開 ABAQUS變量解讀：教你讀懂應(yīng)力/應(yīng)變/損傷
結(jié)語
ABAQUS中的變量雖多，但歸根結(jié)底都是為描述材料行為和判斷結(jié)構(gòu)狀態(tài)服務(wù)的。理解每個(gè)變量的物理意義，結(jié)合實(shí)際分析目標(biāo)選擇合適的輸出，才能讓仿真結(jié)果真正“說話”。
如果你對(duì)某個(gè)變量還有疑問，或者想了解更深入的應(yīng)用場景，歡迎在評(píng)論區(qū)留言！
Abaqus應(yīng)用之能量篇
該文檔主要介紹了Abaqus/Explicit分析中的能量相關(guān)內(nèi)容，包括能量平衡表述、輸出及相關(guān)特點(diǎn)，具體如下：
一、能量平衡表述
1．能量平衡公式：
其中， EI 為內(nèi)能，EV為粘性耗散能，EFD是摩擦耗散能，EKE是動(dòng)能，EW是外加載荷所做的功。這些能量分量的總和為，它必須是個(gè)常數(shù)。在數(shù)值模型中， 只是近似的常數(shù)，一般有小于1％的誤差。
2．內(nèi)能組成
內(nèi)能是能量的總和，它包括可恢復(fù)的彈性應(yīng)變能EE；非彈性過程的能量耗散（例如塑性）EP；粘彈性或者蠕變過程的能量耗散ECD；和偽應(yīng)變能EA：
l 彈性應(yīng)變能EE：可恢復(fù)的能量。
l 非彈性耗散能EP：如塑性變形等過程的能量耗散。
l 粘彈性耗散能ECD：粘彈性或蠕變過程的能量耗散。
l 偽應(yīng)變能EA：包括沙漏阻力及殼和梁單元橫向剪切中的能量，大量偽應(yīng)變能表明需對(duì)網(wǎng)格進(jìn)行細(xì)劃或修改，單元偽應(yīng)變能密度可查看各單元偽應(yīng)變能情況，偽應(yīng)變能ALLSE<5%時(shí)沙漏可控制。
3．粘性能：
由阻尼機(jī)制引起，包括體粘性阻尼和材料阻尼，與粘彈性或非彈性過程耗散能量不同。
4．外力功
由節(jié)點(diǎn)力（力矩）和位移（轉(zhuǎn)角）定義，指定邊界條件也有貢獻(xiàn)。
展開 Abaqus仿真告訴你網(wǎng)球撿球神器為何這么神 附abaqus手冊(cè)線性粘彈性UMAT詳細(xì)解讀下載
Abaqus Model
“撿球神器”Abaqus分析模型
撿球器Abaqus模型的主要部分是沿圓周均布的40根輻條，輻條材料為鋼，模擬采用線彈性材料本構(gòu)，單元類型選用beam；輻條可整體繞中心軸轉(zhuǎn)動(dòng)，接觸到地面、網(wǎng)球時(shí)會(huì)發(fā)生彈性變形。
邊界條件設(shè)置撿球器手柄的下壓、前推、上提過程。
撿球中...
看一下?lián)烨蚱魇侨绾伟亚虺赃M(jìn)去的：
輻條“搭”到網(wǎng)球上
輻條的彈性變形
汽車輪胎花紋的橫紋縫隙里經(jīng)常會(huì)卡到一些比縫隙大的石子，這個(gè)過程某種程度上和上面撿球器的力學(xué)過程是類似的。
月壤樣本采集系統(tǒng)
撿球器是利用輻條的被動(dòng)變形來工作的，下面這個(gè)有點(diǎn)像打蛋器的裝置，利用主動(dòng)變形，調(diào)節(jié)輻條開口縫隙，實(shí)現(xiàn)拾取不同尺寸月壤樣本的功能。
月壤樣本采集系統(tǒng)
采樣爪參數(shù)調(diào)節(jié)
通過改變參數(shù)Ang1在93°到103°之間變化，采樣爪會(huì)進(jìn)入三種不同模式，可以拾取不同尺寸的月壤樣本。
采樣爪的三種模式
下載地址：abaqus手冊(cè)線性粘彈性UMAT詳細(xì)解讀
展開 
ABAQUS能量平衡輸出變量
Total energy output quantities
ALLAE
“Artificial” strain energy associated with constraints used to remove singular modes (such as hourglass control), and with constraints used to make the drill rotation follow the in-plane rotation of the shell elements.
ALLCD
Energy dissipated by creep, swelling, and viscoelasticity.
ALLEE
Electrostatic energy.
ALLFD
Total energy dissipated through frictional effects. (Available only for the whole model.)
ALLIE
Total strain energy. (ALLIE = ALLSE + ALLPD + ALLCD + ALLAE + ALLQB + ALLEE + ALLDMD.)
ALLJD
Electrical energy dissipated due to flow of electrical current.
ALLKE
Kinetic energy.
ALLKL
Loss of kinetic energy at impact. (Available only for the whole model.)
ALLPD
Energy dissipated by rate-independent and rate-dependent
展開 ABAQUS中的能量平衡
2、我分析的是結(jié)構(gòu)在地震波作用下的反應(yīng)，通過查看資料發(fā)現(xiàn)，能量方程有兩種，即相對(duì)能量方程和絕對(duì)能量方程，相對(duì)能量方程是動(dòng)力方程兩邊對(duì)“相對(duì)位移”積分得到的，絕對(duì)能量方程式動(dòng)力方程兩邊對(duì)“絕對(duì)位移”積分得到的，我在ABAQUS幫助文件中沒有找到這些能量項(xiàng)計(jì)算的公式說明，請(qǐng)問這些能量項(xiàng)的計(jì)算公式在哪本資料上能夠查的到？還想再問下，如果求地震波輸入結(jié)構(gòu)中的能量，是否就是外力功ALLWK這項(xiàng)呢？
3、根據(jù)您博客中的內(nèi)容，內(nèi)能=彈性應(yīng)變能（Es）+塑性能（Ep）+阻尼耗能（Ev）+蠕變耗能（Ec），而彈性應(yīng)變能（Es）=損傷耗能（Ed）+可恢復(fù)能（Ee），這里的損傷耗能Ed是ABAQUS中變量ALLDMD（Damage dissipation energy）嗎？它代表什么意思呢？
比對(duì)公式2，如果ALLSE指的是可恢復(fù)的彈性應(yīng)變能，那么是否公式2中少了一項(xiàng)損傷耗能（Es），不知該怎么理解？
展開 abaqus鋼管壓縮及能量輸出
abaqus鋼管壓縮及能量輸出
求abaqus用戶手冊(cè)粘彈性例子詳細(xì)解讀
有沒有大佬懂abaqus用戶手冊(cè)中給出的線性粘彈性的子程序及本構(gòu)方程的詳細(xì)解讀啊，孩子確實(shí)連第一個(gè)公式都看不懂怎么來的。
ABAQUS網(wǎng)格大小對(duì)混凝土本構(gòu)模型影響的案例分析 附Abaqus混凝土材料模型解讀與參數(shù)設(shè)置 V2
本文就以ABAQUS模擬棱柱體混凝土試塊為例，混凝土強(qiáng)度等級(jí)為C110，棱柱體尺寸為100mm*100mm*300mm。（就是我們平常做高強(qiáng)混凝土軸心抗壓強(qiáng)度試塊的尺寸）
模擬數(shù)據(jù)
本文采用受壓本構(gòu)數(shù)據(jù)如下：
本文采用受拉本構(gòu)數(shù)據(jù)如下：
模擬時(shí)網(wǎng)格分別設(shè)為10mm、30mm、50mm和90mm。
加載方式采用在參考點(diǎn)處施加位移的方式，設(shè)置參考點(diǎn)與棱柱體頂面耦合。
邊界條件設(shè)置為與實(shí)際試塊加載的約束條件相同。
模擬結(jié)果
模擬得到的力和位移數(shù)據(jù)經(jīng)過處理，可以得到應(yīng)力和應(yīng)變關(guān)系曲線，如下圖。
從模擬結(jié)果來看，網(wǎng)格大小確實(shí)對(duì)混凝土本構(gòu)有影響。
1，整體趨勢來看，網(wǎng)格越小，混凝土模型表現(xiàn)出的抗壓強(qiáng)度越大，峰值應(yīng)變?cè)叫。_(dá)到峰值后承載力下降越快，相當(dāng)于混凝土越脆。
2，網(wǎng)格10mm和網(wǎng)格30mm的本構(gòu)基本完全相同，但10mm網(wǎng)格的計(jì)算時(shí)間是30mm的8倍。因此采用10mm的網(wǎng)格不太經(jīng)濟(jì)。
3，網(wǎng)格10mm和網(wǎng)格30mm的本構(gòu)峰值強(qiáng)度比原始本構(gòu)下降6.6%，網(wǎng)格50mm的下降了10.5%，網(wǎng)格90mm的下降了11.7%。下降幅度倒是差別不大。
所以網(wǎng)格的大小確實(shí)會(huì)影響模型的響應(yīng)，導(dǎo)致其表現(xiàn)出的本構(gòu)與實(shí)際不同。
下載地址：Abaqus混凝土材料模型解讀與參數(shù)設(shè)置 V2
展開 ABAQUS動(dòng)態(tài)分析中的能量平衡、沙漏及結(jié)果評(píng)估
ABAQUS動(dòng)態(tài)分析中的能量平衡、沙漏及結(jié)果評(píng)估
ABAQUS歷史輸出中，各能量變量（ALLAE、ALLIE、ETOTAL等）意義
ABAQUS中，對(duì)于很多動(dòng)態(tài)問題，尤其像高速?zèng)_擊模擬中，對(duì)結(jié)果評(píng)價(jià)很重要的一點(diǎn)就是要保證模型能量守恒，這就涉及到ABAQUS歷史輸出中各能量變量的意義，下面最各簡單整理：
ALLAE：人工偽應(yīng)變能，六面體、殼網(wǎng)格中沙漏發(fā)生情況指標(biāo)
ALLCD：蠕變、膨脹以及粘彈性消耗的能量
ALLFD：摩擦消耗的能量
ALLIE：總的內(nèi)能，ALLIE=ALLSE+ALLPD+ALLCD+ALLAE+ALLQB+ALLEE+ALLDMD
ALLKE：動(dòng)能
ALLKL：碰撞消耗的能量
ALLVD：粘性消耗的能量
ALLDMD：破壞消耗的能量
ETOTAL：所有能量的總和
NOTE：ALLAE（偽應(yīng)變能）的理解：偽應(yīng)變能就是控制沙漏變形所消耗的主要能量。如果偽應(yīng)變能過高，說明過多的應(yīng)變能被用來控制沙漏變形了。一般通過比較偽應(yīng)變能和其他內(nèi)部能量的值來判斷偽應(yīng)變能是否過高，以及判斷過高的偽應(yīng)變能的來源。
展開 
Abaqus準(zhǔn)靜態(tài)小例子: 能量平衡、質(zhì)量放大
Abaqus準(zhǔn)靜態(tài)問題小例子
1，定義：將顯示問題用于準(zhǔn)靜態(tài)問題需要一些特殊的考慮，由于準(zhǔn)靜態(tài)問題一般是較長時(shí)間的求解過程，它將需要大量的時(shí)間增量，為了獲得較經(jīng)濟(jì)的解答，采取一些措施是必要的。但是帶的問題是隨著加載速度的增加，慣性力起著越來越多的作用。因此必須保證慣性力的影響不至于影響到結(jié)果的準(zhǔn)確度。
那么加載的速度應(yīng)該是多大呢？ 我們知道，結(jié)構(gòu)的最低階模態(tài)通常控制著結(jié)構(gòu)的響應(yīng)。如果知道結(jié)構(gòu)的最低階模態(tài)，我們就可以估計(jì)出適當(dāng)?shù)撵o態(tài)響應(yīng)所需要的時(shí)間，通常理想的加載時(shí)間的最低階模態(tài)周期的10倍，以確保解答是真正的準(zhǔn)靜態(tài)。在加載過程中，保證施加的載荷光滑性是非常重要的。光滑幅值曲線為我們提供了一個(gè)好的方法。
2 質(zhì)量放大技術(shù)，
利用中心差分法求解時(shí)，解的穩(wěn)定性是時(shí)間步長t必須小于該問題求解方程性質(zhì)中某個(gè)臨近值t1—聲波通過該單元的時(shí)間.
網(wǎng)格中的最小單元尺寸將決定t1的大小，網(wǎng)格尺寸越小，t1越小從而使計(jì)算量越大。t1還可以表示成，t1=L/C ，其中C=sqrt（E/P）—聲波的傳輸速度
P是材料的密度，E是彈性模量，泊松比假設(shè)為0.從這個(gè)方程上我們可以看出，如果網(wǎng)格中某個(gè)單元的尺寸過小，計(jì)算量將會(huì)產(chǎn)生不合理的增加。如果我們將材料密度認(rèn)為的增加f^2倍，波速C就會(huì)降低f倍，臨界穩(wěn)定時(shí)間t1就會(huì)增加f倍。所需要的時(shí)間就會(huì)相應(yīng)的減少。對(duì)質(zhì)量縮放的部分，最好是集中在單元尺寸小的網(wǎng)格上。
3，能量平衡
評(píng)估模擬是否產(chǎn)生了正確的準(zhǔn)靜態(tài)響應(yīng)，當(dāng)模型太復(fù)雜時(shí)，單純的通過應(yīng)力作用不明顯，通常的作法是通過對(duì)比能量的角度，能量平衡方程式：
E1+EV+EKE+EFD-EW=Etotal=constant
其中E1是內(nèi)能增量（包括彈性和塑性應(yīng)變能），EV是粘性耗散吸收的能量，EKE是系統(tǒng)的動(dòng)能，EFD是摩擦耗散的能量，EW是外力作的功，Etotal系統(tǒng)中的總能量。
展開 【公益帖】abaqus中的ETOTAL 能量澄清，解釋為何有時(shí)為負(fù)？
顯式計(jì)算中，觀察ETOTOL發(fā)現(xiàn)有時(shí)能量為負(fù)，覺得不太可能，主要原因是大家對(duì)ETOTAL的物理含義有誤解，在abaqus的幫助文檔中 4.2.1 Abaqus/Explicit output variable identifiers給出了ETOTAL物理含義，代表的是總的平衡能，具體代表什么呢？ 如下
ETOTAL 總的能量平衡
ETOTAL = ALLKE + ALLVD + ALLSD + ALLKL + ALLFD + ALLJD + ALLIE – ALLWK
ALLKE 動(dòng)能
ALLVD 耗散能（粘性引起，不包括 ALLSD 和 ALLCD）
ALLSD 耗散能（自動(dòng)穩(wěn)定引起，如接觸）
ALLKL 動(dòng)能損失（沖擊引起）
ALLFD 耗散能（摩擦引起）
ALLJD 耗散能（電流引起）
ALLIE 總應(yīng)變能
ALLIE = ALLSE + ALLPD + ALLCD + ALLAE + ALLQB + ALLEE + ALLDMD
ALLSE 應(yīng)變能(可恢復(fù))
ALLPD 耗散能（塑性變形引起）
ALLCD 耗散能（粘彈性、蠕變、膨脹引起）
ALLAE 偽應(yīng)變能
ALLQB 耗散能（無限單元引起，如無反射邊界）
ALLEE 靜電能
ALLDMD 耗散能（裂紋引起）
ALLWK 外力功
那么出現(xiàn)了ETOTAL能量為負(fù)，并不是代表計(jì)算就有問題，總能來說ETOTAL月接近零越好。但是有時(shí)可能是一個(gè)比較大的值，這是否就說明計(jì)算有問題嗎？ 也不一定要看ETOTAL與 總動(dòng)能（或者總勢能）的比值是否可以比擬，如果在1%以內(nèi)我們認(rèn)為也是可以接受的。
展開 ABAQUS輸出能量說明（Total energy output）——轉(zhuǎn)自公眾號(hào)CAE仿真空間
拿沙漏現(xiàn)象來講，它是零能模式，需要人為控制，增加一定的剛度，由此產(chǎn)生的能量并不是模型自身存在的，所以屬于偽應(yīng)變能；而殼單元的面內(nèi)轉(zhuǎn)動(dòng)自由度屬于“虛假自由度”，是為了防止整體剛度矩陣奇異而引入的，所以由此而帶來的能量也不屬于系統(tǒng)自身產(chǎn)生的，也屬于為應(yīng)變能。
我們都知道，當(dāng)發(fā)現(xiàn)偽應(yīng)變能過大時(shí)候，一般超過5%則需要考慮細(xì)化網(wǎng)格或做其它相應(yīng)處理，否則計(jì)算結(jié)果可能存在問題。
ALLCD—粘彈性耗散能
材料粘滯性和彈性綜合作用耗散的能量，不包括線性粘彈性規(guī)律的超彈性等材料。
ALLFD—摩擦耗散能
整個(gè)模型因摩擦作用耗散的能量，僅對(duì)整個(gè)模型而言。
ALLKE—?jiǎng)幽?ALLPD—非彈性耗散能
幫助里寫的是率無關(guān)或率相關(guān)的塑性變形而耗散的能量
ALLSE—可恢復(fù)的應(yīng)變能
ALLVD—粘性耗散能
主要是因系統(tǒng)粘性阻尼和材料阻尼耗散的能量
ALLWK—外力做的功
僅針對(duì)整個(gè)系統(tǒng)模型而言
ALLIHE—內(nèi)部熱能
ALLHF—外部對(duì)流的熱能
ALLDMD—裂紋損傷導(dǎo)致的能量損失
ALLDC—單元扭曲控制耗散的能量
ALLFC—流體腔能
針對(duì)整個(gè)模型而言，是所有流體腔做的負(fù)功（？？什么意思？？懂的朋友可以留言講講）
ALLPW—罰函數(shù)接觸產(chǎn)生的能量
包含通用接觸對(duì)和罰函數(shù)/運(yùn)動(dòng)學(xué)接觸對(duì)，僅對(duì)整個(gè)模型有效。
展開 