非線性子程序開發(fā)
關(guān)注創(chuàng)建者：匿名 創(chuàng)建時(shí)間：2026-01-05
非線性子程序開發(fā)的視頻教程
非線性各向同性強(qiáng)化彈塑性umat子程序教程
本課程針對最一般的非線性強(qiáng)化彈塑性本構(gòu)，講述了材料在復(fù)雜應(yīng)力情況下何時(shí)屈服，如何強(qiáng)化，怎樣塑性流動(dòng)等彈塑性力學(xué)知識，在此基礎(chǔ)上又講述了如何將理論變成程序，同時(shí)通過一個(gè)簡單的算例來驗(yàn)證本umat子程序的正確性。為了幫助大家編程和調(diào)試，真實(shí)再現(xiàn)了作者在編寫子程序所遇到的問題以及如何調(diào)試解決的方案，希望本次課程能給有些編程基礎(chǔ)的umat子程序用戶提供一些幫助！
￥30 1小時(shí)33分鐘 5645播放
查看
采用非線性彈簧模擬鋼管混凝土柱粘結(jié)滑移（不用子程序）
本人碩士與博士課題都是與鋼-混結(jié)構(gòu)相關(guān)的，此次采用非線性彈簧來做了一個(gè)鋼管混凝土粘結(jié)滑移的案例，案例基于本人的一次試驗(yàn)，有限元模擬與試驗(yàn)擬合較好。如果你的課題與粘結(jié)滑移相關(guān)（不僅僅是鋼管混凝土），這個(gè)視頻全部搞定一步到位！此次視頻只需兩杯奶茶錢，肯定會對你有所幫助，解開心中的疑惑，視頻適合新手入門以及老手進(jìn)階，歡迎大家購買觀看，有問題可以留言或者私信。
￥99 42分鐘 1505播放
查看
【A01】Abaqus 子程序開發(fā)入門=安裝教程+子程序關(guān)聯(lián)+子程序開發(fā)工具+示例
介紹Abaqus子程序開發(fā)的入門級教程。 完整視頻包括：1、軟件安裝；2、Abaqus子程序使用環(huán)境配置；3、子程序開發(fā)環(huán)境配置；4、子程序開發(fā)工具（Abaci）教程；5、復(fù)刻一個(gè)Abaqus UMAT子程序示例。 子程序語言為Fortran（自由格式）。課程附件是視頻中用到的文件。 最后1個(gè)大視頻里有5個(gè)章節(jié)，是上面5個(gè)視頻的融合，并對其中一些細(xì)節(jié)做了補(bǔ)充更新。
￥28.88 1小時(shí)22分鐘 495播放
查看
非線性子程序開發(fā)的實(shí)例教程
問題描述：假設(shè)某種材料，彈性模量E=200GPa，泊松比u=0.3，非線性強(qiáng)化特性如下曲線所示，在一個(gè)直徑d=10mm，長度L0=50mm的試樣的一端約束Z向和周向變形（不約束斷面的徑向，這樣可更好的模擬拉伸實(shí)驗(yàn)），另一端拉伸伸長DL=3.2mm，求試樣Z方向的應(yīng)力。
其硬化曲線的方程為
自定義材料模型為
加載方式為
計(jì)算結(jié)果為
對于這個(gè)問題，通過簡單計(jì)算可以發(fā)現(xiàn)試樣已經(jīng)發(fā)生塑性變形，通過自編的Umat子程序計(jì)算最后試樣應(yīng)力為300MPa。我們知道這個(gè)問題是有理論解的，下面我們來求理論解。
假設(shè)z向應(yīng)力為，總的應(yīng)變?yōu)?聯(lián)立后，得方程
解方程，并取大于200的解為試樣的軸向應(yīng)力
基于Umat子程序的計(jì)算結(jié)果與理論值完全一致。
感興趣的朋友，請觀看非線性強(qiáng)化彈塑性umat子程序視頻教程，也歡迎大家下載本次的Abaqus模型文件和*.for文件。
視頻教程鏈接：
http://www.yqgqt.org.cn/college/video/c14070
展開 非線性彈簧參數(shù)與泊松比
! 4個(gè)Maxwell體參數(shù)。props(5)=E1,props(6)=o1, props(7)=E2,props(8)=o2,
! props(9)=E3,props(10)=o3, props(11)=E4,props(12)=o4
REAL*8 TONGSHI
在有限元分析中，我們經(jīng)常會和非線性打交道，如材料非線性、幾何非線性、邊界非線性。非線性有限元一直是有限元中較為困難的一部分，在非線性有限元中我們經(jīng)常碰到諸如Newton-Raphson迭代法，切線剛度陣等概念，今天我們就單的介紹一下非線性吧。
1.簡單實(shí)例
首先看一個(gè)簡單的彈簧桿件結(jié)構(gòu)，如圖所示，中間節(jié)點(diǎn)作用一個(gè)F的力，會產(chǎn)生一個(gè)位移v
由靜力平衡關(guān)系可得到
該方程為典型的非線性方程，對于這個(gè)方程，如果給定一個(gè)位移v就能求得F，如下圖所示，從圖中曲線可以看到非線性的含義了。圖中不同k對應(yīng)的曲線，可以看到k比較小時(shí)，桿內(nèi)力起主要作用，呈現(xiàn)出幾何非線性，K較大時(shí)，彈簧起主要作用，呈現(xiàn)出彈簧的線彈性。
2.牛頓迭代法
但是在實(shí)際中，我們往往是不知道位移v的，而是知道F，那么給定一個(gè)F，怎么求v呢？這時(shí)候牛頓迭代法就要上場了。牛頓迭代法的思想是將非線性方程線性化，以線性方程的解逼近非線性方程的解，具體操作如下：
牛頓迭代法圖形解釋
對于非線性方程f(x)=的迭代解法有如下格式
3.非線性有限元迭代法
雖然上文只是簡單的一維問題，但是我們可以把它當(dāng)做位移法有限元的原型，對于一般有限元，離散平衡方程一般具有如下形式：
對于試探解、一般有
該方程的求解有如下形式
（1）直接迭代法
直接迭代法中要求K矩陣為u的顯式函數(shù)，只適用于和變形歷史無關(guān)的非線性問題。
展開 問題描述：假設(shè)某種材料，彈性模量E=200GPa，泊松比u=0.3，線性強(qiáng)化特性如下曲線所示，在一個(gè)直徑d=10mm，長度L0=50mm的試樣的一端約束Z向和周向變形（不約束斷面的徑向，這樣可更好的模擬拉伸實(shí)驗(yàn)），另一端拉伸伸長DL=0.9mm，求試樣Z方向的應(yīng)力。
對于這個(gè)問題，通過簡單計(jì)算可以發(fā)現(xiàn)試樣已經(jīng)發(fā)生塑性變形，通過自編的Umat子程序計(jì)算最后試樣應(yīng)力為509.1MPa。我們知道這個(gè)問題是有理論解的，下面我們來求理論解。
先求硬化系數(shù)H
總的應(yīng)變?yōu)?聯(lián)立后，試樣的軸向應(yīng)力為
基于Umat子程序的計(jì)算結(jié)果與理論值完全一致。
接下來請大家觀看Umat子程序逐句編寫視頻，也歡迎大家下載本次的Abaqus模型文件和*.for文件。
或者觀看視頻教程。
http://www.yqgqt.org.cn/college/video/c14014
展開 2026年清明節(jié)前完成了非線性rbe3也就是Abaqus考慮有限變形下分布耦合的開發(fā)工作。從計(jì)算結(jié)果和迭代效率來看，做到了和Abaqus的完全一致性。大變形rbe3的開發(fā)難度要遠(yuǎn)大于大變形rbe2，加上連接單元的開發(fā)，現(xiàn)在已初步具備了有限轉(zhuǎn)動(dòng)下一階二階變分求導(dǎo)的能力。
非線性子程序開發(fā)的相關(guān)專題、標(biāo)簽、搜索
非線性子程序開發(fā)的最新內(nèi)容
原始文獻(xiàn)：《A multiscale simulation framework of the accumulative roll bonding process accounting for texture evolution》
DOI：10.1016/j.msea.2015.02.005
在金屬材料研究領(lǐng)域，粘塑性自洽模型（Visco-Plastic Self-Consistent, VPSC
非線性分布耦合的開發(fā)1個(gè)月前
2026年清明節(jié)前完成了非線性rbe3也就是Abaqus考慮有限變形下分布耦合的開發(fā)工作。從計(jì)算結(jié)果和迭代效率來看，做到了和Abaqus的完全一致性。大變形rbe3的開發(fā)難度要遠(yuǎn)大于大變形rbe2，加上連接單元的開發(fā)，現(xiàn)在已初步具備了有限轉(zhuǎn)動(dòng)下一階二階變分求導(dǎo)的能力。
C3D8單元幾何非線性算法研究及UEL開發(fā)4個(gè)月前
因科研需要，一直在研究一些單元算法，看著網(wǎng)上相關(guān)資料很多，但是和商軟對標(biāo)的非線性單元技術(shù)相對較少。非線性這方面ABAQUS比較受人認(rèn)可，所以打算用空余時(shí)間研究一下ABAQUS的單元技術(shù)，推導(dǎo)編寫一下相關(guān)程序供大家討論。本人水平十分有限，主要是學(xué)習(xí)ABAQUS的文檔，力學(xué)理論和代碼方面的問題請大家不吝賜教。
本文主要推導(dǎo)ABAQUS在幾何非線性（大變形）有限元分析中，用于計(jì)算單元切線剛度矩陣的算法
[圖片]
[圖片]
單向復(fù)合材料在層內(nèi)只有一個(gè)方向的纖維，然而對于平紋復(fù)合材料在層內(nèi)具有經(jīng)向和緯向兩個(gè)方向交織的纖維束。因此在建立均質(zhì)化模型時(shí)，平紋復(fù)合材料的剛度矩陣，損傷起始準(zhǔn)則，損傷演化方法以及退化的剛度矩陣與單向復(fù)合材料具有明顯的差異。主要體現(xiàn)為平紋復(fù)合材料在面內(nèi)的兩個(gè)方向均有纖維，因?yàn)閷τ谄郊y復(fù)合材料的失效模式主要有：經(jīng)向拉伸/壓縮損傷，緯向拉伸/壓縮損傷以及厚度方向上的拉伸/壓縮損傷,此外還可以通過在層間插入
剪切修正模型的數(shù)值實(shí)現(xiàn)------《Nielsen KL, Tvergaard V. Ductile shear failure or plug failure of spot welds modelled by modified Gurson model. Engineering Fracture Mechanics 2010;77:1031–47.》
GTN模型是韌性斷裂的一個(gè)廣為人知的微觀力學(xué)模型
正交各向異性（工程常數(shù)）：
EngineeringConstants-xplD.obj
EngineeringConstants-xpl.obj
共享庫：
explicitU.dll
explicitU-D.dll
參數(shù)輸入依次為：
E1,E2,E3,Nu23,Nu13,Nu12,G23,G13,G12
引言
Abaqus為用戶提供了大量的用戶子程序，用戶可以通過子程序接口開發(fā)Abaqus原本不具備的單元、材料、載荷等。此外，利用子程序能幫助我們?nèi)ダ斫狻⑼茰yAbaqus的計(jì)算過程，這對國產(chǎn)自主軟件對標(biāo)Abaqus功能開發(fā)具有一定的參考意義。
一、安裝Abaqus，并配置使用用戶子程序
Win10(x64)成功安裝VisualStudio2012+Intel Visual Fortran
此次不做測試。
完全各向異性材料的本構(gòu)方程為：
由矩陣可以知道，完全各向異性材料需要21個(gè)獨(dú)立常數(shù)。
此umat的材料參數(shù)設(shè)置參考abaqus，參數(shù)自上而下分別為：
使用方法和之前一致。
若選用縮減積分單元時(shí)要開啟沙漏控制。
完全各向異性umat：
fullyanisotropic-std.obj
