uvarm子程序開發(fā)的案例
Abaqus安全系數(shù)的實(shí)現(xiàn)——UVARM子程序 ￥2
上述實(shí)例和子程序見附件。
Abaqus-UVARM子程序入門教學(xué)------案例十七 ￥99
? Abaqus-UVARM子程序入門教學(xué)
UVARM子程序
該子程序主要用于生成單元輸出
Abaqus文檔中關(guān)于該子程序的描述
變量的含義
1， UVAR是待定義變量
2， RIRECT，T是材料方向信息
3， TIME，DTIME代表時間與時間增量
4， CMNAME，ORNAME代表材料名稱和方向名稱
5， NUVARM代表變量個數(shù)
6， NOEL,NPT代表單元數(shù)量和積分點(diǎn)數(shù)量，
7， LAYER，KSPT代表復(fù)合材料層數(shù)，以及截面積分點(diǎn)數(shù)量
8， KSTEP,KINC分別代表，step-number，和increment-number
9， NDI,NSHR代表應(yīng)力分量
10， COORD表示坐標(biāo)
11， 其余為接口變量
利用該子程序可以計(jì)算一些abaqus內(nèi)部未提供的單元輸出（如應(yīng)力三軸度和lode參數(shù)）
Sigma_h為靜水應(yīng)力，Sigma_eq為mises等效應(yīng)力，Sigma_1（2,3）為材料的主應(yīng)力
程序的部分截圖
在材料屬性界面設(shè)置自定義輸出個數(shù)以及step對應(yīng)的場輸出
單向拉伸應(yīng)力三軸度和lode的理論值
應(yīng)力三軸度的計(jì)算值
Lod參數(shù)的計(jì)算值
理論與計(jì)算完全一致
展開 Abaqus內(nèi)置LaRC05失效準(zhǔn)則子程序（附cohesive單元umat子程序開發(fā)教程）
Abaqus從2017版本開始，內(nèi)置了UVARM和DMGINI兩種子程序供用戶調(diào)用。UVARM子程序需要以“ABQ_LARC05_DMGCRT”作為材料名的前綴，可以用來評估是否滿足LaRC05強(qiáng)度指標(biāo)；DMGINI子程序需要以“ABQ_LARC05_DMGINI”作為材料名前綴，可以結(jié)合XFEM來分析裂紋萌生和擴(kuò)展。
在abaqus中可以通過修改關(guān)鍵字來進(jìn)行模型設(shè)置
UVARM設(shè)置
DMGINI設(shè)置
輸出變量的含義
計(jì)算結(jié)果
此外，如果有小伙伴希望學(xué)習(xí)cohesive單元的umat子程序開發(fā)相關(guān)技術(shù)，可以關(guān)注下圖的教程：
點(diǎn)擊鏈接進(jìn)入了解詳情：Abaqus Cohesive單元介紹及其本構(gòu)的Umat子程序開發(fā)教程
最后，有ABAQUS子程序開發(fā)相關(guān)需求歡迎通過微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室。
展開 技術(shù)鄰Abaqus復(fù)合材料高級培訓(xùn)定制（長期有效）
基礎(chǔ)班第二天主要講解復(fù)合材料失效理論、面內(nèi)漸進(jìn)損傷分析、層間/界面漸進(jìn)損傷分析、三維VUMAT子程序開發(fā)等內(nèi)容；通過提高班課程的學(xué)習(xí)，學(xué)員能夠掌握Abaqus初始損傷分析、漸進(jìn)損傷分析、基于Cohesive內(nèi)聚力模型的層間/界面漸進(jìn)損傷分析、基于VCCT虛擬裂紋閉合技術(shù)的層間/界面裂紋擴(kuò)展分析以及VUMAT子程序開發(fā)與應(yīng)用等。
基礎(chǔ)班課程大部分都基于Abaqus/CAE進(jìn)行操作，一小部分需要開發(fā)程序，適合尚不熟悉Abaqus復(fù)合材料分析的在校學(xué)生及工程單位CAE分析人員。
高級班主要講解多種復(fù)合材料子程序的開發(fā)、Puck失效理論子程序實(shí)現(xiàn)、材料非線性表征、率相關(guān)特定表征等方面的內(nèi)容。通過該課程學(xué)習(xí)，能夠掌握復(fù)合材料UVARM子程序、USDFLD子程序、VUSDFLD子程序、UMAT子程序、VUMAT子程序的開發(fā)，能夠基于VUMAT子程序進(jìn)一步針對Puck失效理論、材料非線性、率相關(guān)性進(jìn)行深入開發(fā)。該課程幾乎全部涉及程序開發(fā)，適合已掌握Abaqus基礎(chǔ)復(fù)合材料建模，在失效理論、程序開發(fā)、論文發(fā)表方面有更深入需求的研究生和科研人員。
【培訓(xùn)要求】：
培訓(xùn)期間，學(xué)員應(yīng)自帶筆記本電腦，電腦中預(yù)先安裝好Abaqus以及與Abaqus版本對應(yīng)的Visual Studio和Intel Visual Fortran，并完成Abaqus的子程序驗(yàn)證。
展開 
abaqus子程序開發(fā)：fortran-c-python混合編程開發(fā)
abaqus子程序開發(fā)語言不局限fortran和c了，還能支持python！
這是一套自研的開發(fā)框架，源于客戶需要獲取分析過程的中間數(shù)據(jù)作為python機(jī)器學(xué)習(xí)程序的輸入這樣的需求。這種需求乍一聽確實(shí)非主流、非常規(guī)。
后來仔細(xì)做了一些了解，這種需求有其合理性。第一，python在數(shù)據(jù)分析、機(jī)器學(xué)習(xí)方面有豐富的資源；第二，python程序嵌入abaqus求解器一起運(yùn)行，那開發(fā)就會很方便高效；第三，借助python的pdb可以隨心所欲的設(shè)置斷點(diǎn)，調(diào)試程序。在此之前，python一般作為獨(dú)立程序，依賴數(shù)據(jù)文件的順序傳遞，實(shí)現(xiàn)與abaqus的協(xié)同。直白的講，abaqus算完了，形成數(shù)據(jù)文件再導(dǎo)入python。
于是火力全開研究解決方案，結(jié)果可謂“大快人心”：實(shí)現(xiàn)了abaqus的子程序對python的支持，而且還實(shí)現(xiàn)了fortran、c、python三種語言的混合編程，真正讓不同的編程語言在它最擅長的領(lǐng)域發(fā)揮作用。
以上介紹了abaqus子程序fortran-c-python混合編程，如有這方面的需求，歡迎私信聯(lián)系開展合作。
展開 【12月7日-9日 西安】Abaqus復(fù)合材料分析技術(shù)培訓(xùn)班-三天
基礎(chǔ)班第二天主要講解復(fù)合材料失效理論、面內(nèi)漸進(jìn)損傷分析、層間/界面漸進(jìn)損傷分析、三維VUMAT子程序開發(fā)等內(nèi)容；通過提高班課程的學(xué)習(xí)，學(xué)員能夠掌握Abaqus初始損傷分析、漸進(jìn)損傷分析、基于Cohesive內(nèi)聚力模型的層間/界面漸進(jìn)損傷分析、基于VCCT虛擬裂紋閉合技術(shù)的層間/界面裂紋擴(kuò)展分析以及VUMAT子程序開發(fā)與應(yīng)用等。
高級班安排在第三天，主要講解多種復(fù)合材料子程序的開發(fā)、Puck失效理論子程序實(shí)現(xiàn)、材料非線性表征、率相關(guān)特定表征、子程序調(diào)試等方面的內(nèi)容。通過該課程學(xué)習(xí)，能夠掌握復(fù)合材料UVARM子程序、USDFLD子程序、VUSDFLD子程序、UMAT子程序、VUMAT子程序的開發(fā)，能夠基于VUMAT子程序進(jìn)一步針對Puck失效理論、材料非線性、率相關(guān)性進(jìn)行深入開發(fā)。該課程幾乎全部涉及程序開發(fā)，適合已掌握Abaqus基礎(chǔ)復(fù)合材料建模，在失效理論、程序開發(fā)、論文發(fā)表方面有更深入需求的研究生和科研人員。
學(xué)員可根據(jù)自身的條件和需求選擇適合自己的課程進(jìn)行學(xué)習(xí)。
【基礎(chǔ)班課程大綱】：
【高級班課程大綱】：
【培訓(xùn)講義】：
該培訓(xùn)采用獨(dú)家內(nèi)部講義，該講義分為20個部分，共計(jì)256頁，每一部分形成獨(dú)立章節(jié)，從基本建模操作、理論講解、公式推導(dǎo)到程序編寫全面覆蓋。
展開 Abaqus用戶子程序安裝、開發(fā)、調(diào)試
引言
Abaqus為用戶提供了大量的用戶子程序，用戶可以通過子程序接口開發(fā)Abaqus原本不具備的單元、材料、載荷等。此外，利用子程序能幫助我們?nèi)ダ斫?、推測Abaqus的計(jì)算過程，這對國產(chǎn)自主軟件對標(biāo)Abaqus功能開發(fā)具有一定的參考意義。
UEL子程序開發(fā)——梁單元
圖 3 變形協(xié)調(diào)條件
子程序講解
編制矩陣操作函數(shù)
為了程序的易讀性，我們可以將經(jīng)常使用的的計(jì)算部分編寫成一個子程序，使用時直接call一下即可，比如在UMAT中，我們會常常用到Mises應(yīng)力公式，所以我們可以將Mises的求解編寫為一個子程序，在UMAT中調(diào)用即可，大大減少了代碼量，更能提升易讀性。
本程序中，我編寫了矩陣置零函數(shù)，用于數(shù)組的初始化，大家在自己的程序中編制矩陣相乘函數(shù)和矩陣求逆函數(shù)，大大提升編寫代碼的效率。
Simufact 二次開發(fā)及子程序
軟件適配以及基礎(chǔ)程序理解 目前已做簡單的微觀動態(tài)再結(jié)晶模型 可按要求做代碼
abaqus 材料子程序開發(fā)
小妹最近想學(xué)習(xí)abaqus材料子程序開發(fā)，請問一下高手，需要哪些軟件方面的支持呢，是不是安裝visual studio和fotran就可以了，我的abaqus是6.7的，另外兩個軟件需要什么版本的呢。多謝多謝
vpsc-vumat子程序的開發(fā)與應(yīng)用文章推薦
結(jié)合作者的理論（尤其是分段線性化和應(yīng)力驅(qū)動的求解思路）我們可以把獨(dú)立的vpsc子程序編寫進(jìn)abaqus里面，為了避免復(fù)雜的雅可比推導(dǎo)，以及適用各種復(fù)雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序
以下展示一個使用vpsc-鎂合金本構(gòu)模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數(shù)取自vpsc官方案例,為了減少計(jì)算時間使用高應(yīng)變率進(jìn)行計(jì)算，應(yīng)力高于vpsc模擬）：
等效塑性應(yīng)變：
第一個晶粒的累計(jì)剪切滑移：
發(fā)生孿晶次數(shù);
變形后的形狀演化：
Abaqus 子程序開發(fā)日志01-umat之isotropic
對有限元進(jìn)行了一些理論上的學(xué)習(xí)，現(xiàn)在正式嘗試對abaqus用戶子程序進(jìn)行嘗試書寫，此貼只作為自己學(xué)習(xí)的記錄以及分享。
-------------------------------------------------------------------------------------------------------------------------------
umat子程序只用于通用靜力分析步使用來定義用戶自定義材料，這里首先對較為簡單的各向同性材料本構(gòu)模型進(jìn)行構(gòu)建。
先說明一下umat的工作原理，再靜力分析時，每增加一個（極小的）微分應(yīng)變增量，通過代入微分應(yīng)力關(guān)于微分應(yīng)變的雅可比矩陣B（j）方程，得到應(yīng)力增量，再與原始應(yīng)力累加得到新的應(yīng)力（也可以說更新應(yīng)力）。至此，此umat的工作到此結(jié)束。
對于各向同性的剛度矩陣方程為：
再來就是對矩陣中未知量進(jìn)行求解，對其而言只需兩個獨(dú)立常量就可以表示。對于各向同性材料而言，G=E/2(1+NU),所以我們只要知道G、E、NU中的任意兩個量就可以得到另外一個量，習(xí)慣上我們常用E、NU來表示各向同性材料的材料屬性，所以我們這里就需要將矩陣中的未知量用E與NU表示出來就可以了。
求解的過程不再說明，此為材料力學(xué)的基礎(chǔ)，有需要可以查找相關(guān)資料。
求解完成后將其改寫成微分的形式得到雅可比矩陣后用程序進(jìn)行表達(dá)出來，并使用abaqus提供的接口進(jìn)行嵌套寫入用戶子程序，至此，全部工作完成。
先對其進(jìn)行檢驗(yàn)，abaqus中可以直接定義各向同性材料的材料屬性，在初期子程序的學(xué)習(xí)中檢驗(yàn)umat的正確性提供了很好的對照。
這里就隨便建立一個長寬高為5；5；20的立方體對其進(jìn)行施加一個單位為5的拉伸位移（較長方向拉伸）。
展開 Abaqus 修正劍橋模型的vumat子程序開發(fā)
考慮到莫 爾–庫侖準(zhǔn)則有此特征，故將莫爾–庫侖準(zhǔn)則與劍橋模型相結(jié)合，得
聯(lián)合式(1)，(4)，(5)，有改進(jìn)的屈服函數(shù)：
根據(jù)上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計(jì)算得到得結(jié)果如下圖所示。
[1]袁克闊,陳衛(wèi)忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應(yīng)的修正劍橋模型及數(shù)值實(shí)現(xiàn)[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(08):1574-1579.
LS-DYNA用戶子程序（UMAT）開發(fā)案例
最近在International Journal of Impact Engineering（中科院一區(qū)）上看到一篇文章，題目是”Plasticity based material model for concrete subjected to dynamic loadings（https://doi.org/10.1016/j.ijimpeng.2020.103581）“，作者開發(fā)了一個適用于描述動態(tài)荷載作用下混凝土力學(xué)響應(yīng)的塑性模型，并在LS-DYNA中通過用戶自定義子程序（UMAT）二次開發(fā)的方式實(shí)現(xiàn)了模型。論文中展示的模擬效果圖如下：
模擬鋼筋混凝土板的爆炸響應(yīng)
我和這個作者要到了他開發(fā)的UMAT源代碼，在這里給大家共享一下，對于想學(xué)習(xí)LS-DYNA子程序開發(fā)的絕對是不可多得的資料（免費(fèi)共享，不負(fù)責(zé)答疑哦）
（希望大家下載了關(guān)注我一下就行，有空可以看看我發(fā)的各種教學(xué)視頻和資料貼，985高校在讀博士，長期穩(wěn)定更新各種內(nèi)容）
UMAT.f
展開 內(nèi)聚力單元二次開發(fā)子程序
支持二維，三維(四面體和六面體)。
