abaqus塑性變形計(jì)算的案例
彈塑性材料分析-殘余變形計(jì)算
對(duì)于塑性材料，當(dāng)結(jié)構(gòu)屈服之后不能恢復(fù)原形，如果沒有設(shè)置塑性參數(shù)，其與彈性材料比較變形和應(yīng)力結(jié)果都有一定差異。
點(diǎn)擊下方超鏈接觀看視頻
彈塑性材料分析-殘余變形計(jì)算
abaqus彈塑性粒狀材料的有限變形
Finite deformation of an elastic-plastic granular material.rar
ABAQUS 在計(jì)算到step3后計(jì)算中斷，變形過大
我在原本可以計(jì)算的模型的基礎(chǔ)上修改的，只是刪了幾個(gè)樓板與梁連接的栓釘
ABAQUS混凝土損傷塑性模型-C30EXCEL計(jì)算表格（含計(jì)算結(jié)果） ￥3.7
計(jì)算表格（如下）中標(biāo)黃部分的參數(shù)可自行設(shè)定后，EXCEL程序會(huì)自動(dòng)計(jì)算“抗拉強(qiáng)度、非彈性應(yīng)變、受拉損傷因子”。
基于Abaqus的高純鋁不同應(yīng)變率下單晶塑性變形的取向依賴性研究
（原因是位錯(cuò)密度模型提供了位錯(cuò)之間更合理的相互作用形式）
作者的研究思路
一，通過實(shí)驗(yàn)獲得兩種取向的單晶在不同應(yīng)變率下的流動(dòng)應(yīng)力，并發(fā)現(xiàn)了單晶變形的流動(dòng)應(yīng)力與應(yīng)變率和取向是強(qiáng)相互影響的
二，為了捕捉這種應(yīng)變率響應(yīng)，作者在huang經(jīng)典程序的技術(shù)上進(jìn)行了修改，提出了兩類新的流動(dòng)和硬化的晶體塑性模型，三類模型分別如下：
類型一：經(jīng)典單晶唯象本構(gòu)方程
類型二：應(yīng)變率修正的單晶唯象本構(gòu)方程（類似于JC的修改形式，考慮溫度和應(yīng)變率）
類型三：考慮熱激活和位錯(cuò)密度的物理本構(gòu)模型
作者詳細(xì)探討了三種模型的差異性和預(yù)測(cè)能力，并對(duì)物理模型的參數(shù)給出了大概的區(qū)間以及參數(shù)的影響，對(duì)于使用位錯(cuò)密度模型提供了很好的范例
三類模型的預(yù)測(cè)結(jié)果如下所示：
模型一的預(yù)測(cè)結(jié)果
模型二的預(yù)測(cè)結(jié)果
模型三的預(yù)測(cè)結(jié)果
作者認(rèn)為模型三對(duì)單晶變形的預(yù)測(cè)能力最好，因?yàn)槠洳蹲搅烁嗟奈锢硖卣鳌?最后，有需要的小伙伴歡迎聯(lián)系。
展開 ABAQUS混凝土損傷塑性模型-C30EXCEL計(jì)算表格（含計(jì)算結(jié)果）
這是根據(jù)GB50010-2010中混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范中的混凝土本構(gòu)模型，結(jié)合文獻(xiàn)所述的損傷因子定義，編制的計(jì)算C30混凝土非彈性應(yīng)變和損傷因子的EXCEL表格。也是邊學(xué)變做，希望能和大家多交流。
C30砼本構(gòu)（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計(jì)算依據(jù).rar
Abaqus+PyQt+Python平面變形歐拉角計(jì)算
下面以簡(jiǎn)單例子介紹平面變形、指向歐拉角的計(jì)算，包括絕對(duì)歐拉角、相對(duì)歐拉角。
1 簡(jiǎn)化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個(gè)面的厚度不同，在上側(cè)3個(gè)面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計(jì)算要求
計(jì)算六面體上面3個(gè)面的變形歐拉角，包括3個(gè)面的絕對(duì)歐拉角，平面2、3相對(duì)與平面1的相對(duì)歐拉角。平面1、2、3如下圖所示。
3 數(shù)據(jù)處理
使用平面節(jié)點(diǎn)坐標(biāo)、位移數(shù)據(jù)計(jì)算平面變形歐拉角?？梢允褂肞ython腳本輸出平面節(jié)點(diǎn)編號(hào)、節(jié)點(diǎn)坐標(biāo)（X、Y、Z）、節(jié)點(diǎn)位移（U1、U2、U3），如下圖所示。下圖為平面1的10個(gè)工況的數(shù)據(jù)文件，打開的文本文件中7列數(shù)據(jù)為節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移。
三個(gè)平面10個(gè)工況的節(jié)點(diǎn)數(shù)據(jù)文件如下圖所示。每個(gè)文件中包含一個(gè)工況一個(gè)平面的節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移數(shù)據(jù)。
4 絕對(duì)歐拉角計(jì)算
使用PyQt+Python開發(fā)了一個(gè)簡(jiǎn)單的小軟件，計(jì)算絕對(duì)歐拉角、相對(duì)歐拉角。
首先計(jì)算各平面的絕對(duì)歐拉角。
計(jì)算平面1的10個(gè)工況的絕對(duì)歐拉角。
平面1變形的絕對(duì)歐拉角計(jì)算結(jié)果如下圖所示。
伴隨絕對(duì)歐拉角計(jì)算結(jié)果，軟件同時(shí)寫出了平面變形前后的坐標(biāo)系數(shù)據(jù)，如下圖。每行18個(gè)數(shù)據(jù)，每3個(gè)數(shù)據(jù)為一個(gè)坐標(biāo)軸向量，變形前后2個(gè)坐標(biāo)系，6個(gè)坐標(biāo)軸，18個(gè)數(shù)據(jù)。
5 相對(duì)歐拉角計(jì)算
利用計(jì)算絕對(duì)歐拉角時(shí)得到的坐標(biāo)系文件，計(jì)算平面變形相對(duì)歐拉角，如下圖所示，計(jì)算平面2相當(dāng)于平面1、平面3相對(duì)與平面1的相對(duì)歐拉角。
計(jì)算結(jié)果如下圖所示。
6 小結(jié)
上述軟件用的算法申請(qǐng)了發(fā)明專利，軟件申請(qǐng)了軟著。
展開 ABAQUS網(wǎng)格劇烈變形導(dǎo)致計(jì)算終止
模型就是封面那樣，是一個(gè)螺紋連接二維對(duì)稱模型，我要算多孔金屬塑性，就是GTN模型，然后我在1這邊加了一個(gè)位移載荷，2這邊加一個(gè)U1的位移轉(zhuǎn)角，3這里是完全約束。
各相互接觸的螺紋牙之間都設(shè)置的面與面接觸，具體參數(shù)沒有設(shè)置，都是用的系統(tǒng)默認(rèn)。
但是就算我載荷設(shè)置的再小，都會(huì)出現(xiàn)這種錯(cuò)誤
The ratio of deformation speed to wave speed exceeds 1.0000 in at least one element. This usually indicates an error with the model definition. Additional diagnostic information may be found in the message file.
展開 Abaqus如何保存變形后的模型并在此基礎(chǔ)上計(jì)算
Abaqus如何保存變形后的模型并在此基礎(chǔ)上計(jì)算
在Abaqus中計(jì)算完成后會(huì)生成odb文件，也就是結(jié)果文件。如果想重新采用變形后的模型進(jìn)行相關(guān)分析，可以通過Abaqus的導(dǎo)入操作完成。
新建一個(gè)model，通過File下拉的import—Part，之后選擇相應(yīng)的odb模型，可以自己選擇導(dǎo)入哪個(gè)載荷步或者哪個(gè)子步的變形結(jié)果。
圖1
很明顯，導(dǎo)入后的模型如圖2所示。而之前的模型如圖3所示。
圖2
圖3
導(dǎo)進(jìn)來之后的模型需要重新定義材料、加載組件、定義載荷步、加載，之后可以選擇在編寫基礎(chǔ)之上繼續(xù)計(jì)算。
如果只是想得到變形后的單元節(jié)點(diǎn)信息，則可以不考慮在Abaqus中定義材料和加載等，而是直接在Job模塊直接創(chuàng)建一個(gè)Job，通過Write input生成記錄了該變形的inp文件。該inp文件里面的信息是單元和節(jié)點(diǎn)信息。
Inp文件可以作為數(shù)據(jù)文件進(jìn)行傳遞，對(duì)于HyperMesh，可以與Abaqus進(jìn)行有效的連接，可以讀取Abaqus的inp文件生成相關(guān)信息。
展開 線性強(qiáng)化彈塑性umat子程序-umat在abaqus計(jì)算流程中的意義及調(diào)試方法
abaqus子程序用fortran語言編寫，且子程序又通過abaqus調(diào)用，因而涉及到很多的力學(xué)及有限元理論知識(shí)，使得調(diào)試難度陡然增加。那么abaqus子程序常用的調(diào)試方法是什么呢？其實(shí)說出來也很簡(jiǎn)單，就是打印變量法，即把程序中的變量打印出來（一般打印結(jié)果在當(dāng)前工作目錄下的job-name.log的文件里），然后考察其合理性。
如在上述umat子程序中定義了一個(gè)save在內(nèi)存的變量，用來記錄進(jìn)入umat的次數(shù)。
直接打印變量的方法很有用，但是當(dāng)單元數(shù)增加后，眾多的打印信息形成很多的干擾因素，這時(shí)候我們需要控制打印的頻率，常用的方法為針對(duì)某個(gè)單元的某個(gè)積分點(diǎn)打印，如下面的一段程序所示，紅色圓圈里內(nèi)容的含義是當(dāng)單元編號(hào)為1且積分點(diǎn)編號(hào)為1時(shí)才打印相關(guān)信息，這樣調(diào)試更具針對(duì)性。
我們通過打印信息發(fā)現(xiàn)，一次迭代二次進(jìn)入umat，第一次進(jìn)入umat僅是給abaqus返回雅可比矩陣，第二次進(jìn)入umat的目的是為了更新應(yīng)力等信息。
視頻教程有這個(gè)帖子的更詳細(xì)解釋，感興趣的點(diǎn)擊下面的鏈接觀看
http://www.yqgqt.org.cn/college/video/c14014
展開 ABAQUS混凝土損傷塑性模型-2010混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范中C50混凝土-彈模34400Mpa-損傷因子計(jì)算及EXCEL
這是我自己計(jì)算的2010規(guī)范中ABAQUS混凝土損傷塑性模型-2010混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范中C50混凝土-彈模34400Mpa-損傷因子計(jì)算及EXCEL
首先用自己的數(shù)據(jù)計(jì)算2010規(guī)范中規(guī)定的混凝土本構(gòu)關(guān)系
然后借助文件夾中02版規(guī)范的方法，計(jì)算損傷因子。
以后還會(huì)有詳細(xì)計(jì)算方法，此數(shù)據(jù)僅供參考。
2010規(guī)范用C50混凝土損傷塑性本構(gòu)關(guān)系數(shù)據(jù)-彈模34400MPa-帶損傷因子-自己數(shù)據(jù)計(jì)算得出.rar
