abaqus uel單元的案例
abaqus中uel單元如何施加重力
1000是最模型的高度
最后的結(jié)果比真實(shí)值小一半，求助大佬解疑
ABAQUS用戶自定義單元UEL學(xué)習(xí)資料 ￥49
ABAQUS子程序UEL
ABAQUS UEL - 損傷材料本構(gòu)簡(jiǎn)單應(yīng)用于4節(jié)點(diǎn)平面單元 ￥300
利用ABAQUS自定義單元子程序，既可以開(kāi)發(fā)新的單元，同時(shí)也可以定義新的材料本構(gòu)模型。本文以損傷模型簡(jiǎn)單應(yīng)用于4節(jié)點(diǎn)平面單元為案例，介紹ABAQUS UEL的開(kāi)發(fā)和使用。
如上圖所示，該單元包含4個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)有兩個(gè)自由度，分別在水平（X）和垂直（Y）方向運(yùn)動(dòng)。節(jié)點(diǎn)1的兩個(gè)自由度被固定，節(jié)點(diǎn)4的水平自由度被固定，節(jié)點(diǎn)2的垂直自由度被固定。節(jié)點(diǎn)3和節(jié)點(diǎn)4在垂直方向上向上運(yùn)動(dòng)，位移為0.1mm。該正方形單元的邊長(zhǎng)為100mm。在input文件里，坐標(biāo)表示為，
定義節(jié)點(diǎn)組合與邊界條件為，
為了讓模型收斂性更好，采用quasi-newton 求解器。時(shí)間步設(shè)置為，
在文件夾中通過(guò)Powershell提交job和子程序，
單個(gè)單元的變形為，
采用不同的 ??
，在后處理中得到損傷因子的變化，
相對(duì)應(yīng)的力-時(shí)間關(guān)系為，
對(duì)于多個(gè)單元的情況，比如9單元組成的模型，
具體介紹見(jiàn)知乎：ABAQUS UEL - 損傷材料本構(gòu)簡(jiǎn)單應(yīng)用于4節(jié)點(diǎn)平面單元 - 知乎 (zhihu.com)
相對(duì)應(yīng)的UEL代碼和input文件在付費(fèi)內(nèi)容中，
展開(kāi) ABAQUS用戶定義單元UEL與VUEL從入門(mén)到放棄系列1 ￥20
最初來(lái)技術(shù)鄰，講實(shí)話其實(shí)就是來(lái)學(xué)UEL的，后來(lái)發(fā)現(xiàn)沒(méi)帖子，然后就各種找論壇，simweABAQUS版塊的帖子二零零幾年的帖子都被我挖干凈了，做UEL的人確實(shí)比做UMAT的少多了，但技術(shù)鄰這個(gè)平臺(tái)依舊讓我受益良多，諸如藍(lán)牙老師，isolver團(tuán)隊(duì)snowwave02等大佬，他們的分享，甚至很多是免費(fèi)分享確實(shí)高屋建瓴，知識(shí)的學(xué)習(xí)必是知其所以然，我想只有這樣的學(xué)術(shù)氛圍越來(lái)越重，這樣樂(lè)于分享的大佬越來(lái)越多，中國(guó)的自研工業(yè)軟件才能走向世界，我們的工業(yè)4.0才有可能實(shí)現(xiàn)。
學(xué)于此，也該傳于此，作為一只科研狗，我想我比誰(shuí)都清楚什么叫從“入門(mén)”到“放棄”，未知的知識(shí)就是這樣，不懂的時(shí)候它比金子都貴，等你學(xué)會(huì)了它好像又像水一樣普通，亦如水一樣珍貴。
之后也會(huì)逐漸在技術(shù)鄰更一些帖子，可能也會(huì)有視頻吧，內(nèi)容大概會(huì)圍繞有限元理論，用戶子程序開(kāi)發(fā)，ABAQUS python開(kāi)發(fā)，XFEM，lamb波的SHM等，因?yàn)椴恢老雽W(xué)UEL\VUEL的人群多大，也不知道大家都想了解些什么，不知道應(yīng)該從何說(shuō)起，可能有些想學(xué)UEL的小白基礎(chǔ)很差連有限元理論都沒(méi)學(xué)全，再加上我確實(shí)是個(gè)懶癌晚期，其實(shí)很多東西都做好了，懶的發(fā)，也希望能治療一下自己的懶癌吧。
今天的正事，發(fā)兩個(gè)UEL的Fortran程序和對(duì)應(yīng)的inp算例，我發(fā)現(xiàn)好像很多人不會(huì)寫(xiě)UEL的inp，下次有機(jī)會(huì)給你們講講吧（懶癌警告）。
一個(gè)是平面應(yīng)力\平面應(yīng)變單元，2d4node程序
一個(gè)是3D stress單元，C3D8程序
展開(kāi) 
ABAQUS UEL - Embedded crack model 在CST單元中的簡(jiǎn)單實(shí)現(xiàn) ￥500
為了采用有限元模擬材料的開(kāi)裂，cohesive單元常被預(yù)設(shè)在連續(xù)體單元之間，通過(guò)cohesive塑性或者損傷本構(gòu)模擬出材料強(qiáng)度下降的過(guò)程。但是，這種方式使得裂縫只能在預(yù)設(shè)cohesive單元的位置開(kāi)裂。為了更好的模擬連續(xù)體單元的開(kāi)裂，Embedded crack model假設(shè)在單元內(nèi)部存在開(kāi)裂面，并通過(guò)cohesive zone model定義開(kāi)裂面的本構(gòu)關(guān)系。開(kāi)裂面的角度通過(guò)最大主應(yīng)力確定。這種方式極大地方便了網(wǎng)格劃分過(guò)程。 本文通過(guò)3節(jié)點(diǎn)CST單元，介紹如何使用UEL子程序在ABAQUS中實(shí)現(xiàn)簡(jiǎn)單的Embedded crack model。
1. Embedded crack model 的概念
假設(shè)有一個(gè)三角形單元如下，
在節(jié)點(diǎn)1的水平位置施加位移0.05，單元的變形為，
對(duì)于延展性良好的材料，上述變形發(fā)生時(shí)材料發(fā)生塑性變形；但是對(duì)于脆性材料，單元發(fā)生開(kāi)裂，開(kāi)裂的效果為下圖，
開(kāi)裂面上的坐標(biāo)系存在法向（n）和切向（s）行為。材料發(fā)生開(kāi)裂時(shí)，開(kāi)裂面的本構(gòu)是非線性的，而其他連續(xù)體部分依然是線彈性的。
相比起將整個(gè)單元都定義為非線性，這種將非線性開(kāi)裂面內(nèi)嵌在線彈性連續(xù)體內(nèi)的方法更接近脆性材料的性質(zhì)。
2. 開(kāi)裂面本構(gòu)關(guān)系
連續(xù)體線彈性的本構(gòu)采用Plane strain relation，在大部分有限元參考書(shū)中都有介紹，其剛度矩陣為，
開(kāi)裂面的本構(gòu)建立在塑性力學(xué)框架中，采用了簡(jiǎn)單的 tension cut-off 屈服準(zhǔn)則。
tn 表示開(kāi)裂面上法向應(yīng)力，ft 則是材料的受拉強(qiáng)度，可以表達(dá)為，
3.
展開(kāi) ABAQUS UEL-梯度損傷模型應(yīng)用于4節(jié)點(diǎn)平面單元 ￥600
傳統(tǒng)損傷模型對(duì)于單元的尺寸十分敏感，不同單元尺寸會(huì)導(dǎo)致有限元模型精度出現(xiàn)明顯偏差。針對(duì)該問(wèn)題，梯度損傷（Gradient-damage）模型的概念被提了出來(lái)。
本文詳細(xì)介紹了如何將梯度損傷模型應(yīng)用于4節(jié)點(diǎn)平面單元，并在有限元模型中進(jìn)行模擬。
ABAQUS提供了UEL（user defined element）給使用者進(jìn)行開(kāi)發(fā)。筆者利用UEL開(kāi)發(fā)4節(jié)點(diǎn)平面單元，其邊界條件如下圖所示。其中，節(jié)點(diǎn)1的X、Y方向被限制住，節(jié)點(diǎn)2的Y方向被限制，節(jié)點(diǎn)4的X方向被限制，節(jié)點(diǎn)3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個(gè)節(jié)點(diǎn)除了X和Y方向的位移，還帶有非局部應(yīng)變（nonlocal strain）。
單個(gè)單元模型，
多個(gè)單元模型，
具體內(nèi)容可參見(jiàn)知乎文章：
ABAQUS UEL-梯度損傷模型應(yīng)用于4節(jié)點(diǎn)平面單元 - 知乎 (zhihu.com)
相應(yīng)的input文件和uel代碼付費(fèi)可見(jiàn)，
展開(kāi) abaqus最簡(jiǎn)單的uel子程序自定義單元例子 ￥3
通過(guò)abaqus的uel子程序定義了2節(jié)點(diǎn)單元的剛度矩陣，直接運(yùn)行分析可以得到節(jié)點(diǎn)應(yīng)力、位移等參數(shù)，可以快速了解uel的組織架構(gòu)。
運(yùn)行方法：abaqus創(chuàng)建job——來(lái)源——點(diǎn)擊輸入文件——加載inp文件——編輯作業(yè)——通用模塊——加載用戶子程序for文件——運(yùn)行
ABAQUS用戶定義單元UEL與VUEL從入門(mén)到放棄系列2 ￥50
最初來(lái)技術(shù)鄰，講實(shí)話其實(shí)就是來(lái)學(xué)UEL的，后來(lái)發(fā)現(xiàn)沒(méi)帖子，然后就各種找論壇，simweABAQUS版塊的帖子二零零幾年的帖子都被我挖干凈了，做UEL的人確實(shí)比做UMAT的少多了，但技術(shù)鄰這個(gè)平臺(tái)依舊讓我受益良多，諸如藍(lán)牙老師，isolver團(tuán)隊(duì)snowwave02等大佬，他們的分享，甚至很多是免費(fèi)分享確實(shí)高屋建瓴，知識(shí)的學(xué)習(xí)必是知其所以然，我想只有這樣的學(xué)術(shù)氛圍越來(lái)越重，這樣樂(lè)于分享的大佬越來(lái)越多，中國(guó)的自研工業(yè)軟件才能走向世界，我們的工業(yè)4.0才有可能實(shí)現(xiàn)。
學(xué)于此，也該傳于此，作為一只科研狗，我想我比誰(shuí)都清楚什么叫從“入門(mén)”到“放棄”，未知的知識(shí)就是這樣，不懂的時(shí)候它比金子都貴，等你學(xué)會(huì)了它好像又像水一樣普通，亦如水一樣珍貴。
之后也會(huì)逐漸在技術(shù)鄰更一些帖子，可能也會(huì)有視頻吧，內(nèi)容大概會(huì)圍繞有限元理論，用戶子程序開(kāi)發(fā)，ABAQUS python開(kāi)發(fā)，XFEM，lamb波的SHM等，因?yàn)椴恢老雽W(xué)UEL\VUEL的人群多大，也不知道大家都想了解些什么，不知道應(yīng)該從何說(shuō)起，可能有些想學(xué)UEL的小白基礎(chǔ)很差連有限元理論都沒(méi)學(xué)全，再加上我確實(shí)是個(gè)懶癌晚期，其實(shí)很多東西都做好了，懶的發(fā)，也希望能治療一下自己的懶癌吧。
今日份正事，給大家簡(jiǎn)單講一講ABAQUS中的壓電耦合單元C3D8E。
那么什么叫壓電耦合單元呢，簡(jiǎn)單的說(shuō)就是你給它加載電壓（電勢(shì)的差，ABAQUS中為位移加載），那么就會(huì)引起單元力學(xué)場(chǎng)的變化，比如位移、應(yīng)力、應(yīng)變等等；同樣的你給它加載力或者位移，亦會(huì)引起單元電場(chǎng)的變化。
單元壓電耦合場(chǎng)的廣義本構(gòu)方程表示如下：
即為：廣義應(yīng)力=廣義彈性矩陣·廣義應(yīng)變。其中廣義力中D為電位移，廣義應(yīng)變中E為電場(chǎng)強(qiáng)度。
在廣義彈性矩陣中，C矩陣為力學(xué)場(chǎng)的彈性矩陣，e矩陣為壓電常數(shù)矩陣，右下角為介電常數(shù)矩陣。
展開(kāi) ABAQUS 二次開(kāi)發(fā)（UEL+UMAT C3D8 單元動(dòng)力隱式計(jì)算）
概述：采用UEL接口二次開(kāi)發(fā)實(shí)現(xiàn)八節(jié)點(diǎn)單元，考慮BBAR修正，避免體積自鎖，對(duì)標(biāo)ABAQUS自帶的C3D8單元，計(jì)算的剛度矩陣、質(zhì)量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進(jìn)行應(yīng)力和應(yīng)變數(shù)據(jù)的可視化，計(jì)算的應(yīng)力應(yīng)變數(shù)據(jù)同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動(dòng)力簡(jiǎn)諧荷載為例，采用上述程序，應(yīng)用動(dòng)力隱式計(jì)算分析步，最終計(jì)算的位移、應(yīng)變等時(shí)程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節(jié)點(diǎn)均固定。另一面所有的節(jié)點(diǎn)施加簡(jiǎn)諧荷載，簡(jiǎn)諧荷載曲線為（詳細(xì)的參數(shù)見(jiàn)附件）：
計(jì)算的網(wǎng)格圖為：
網(wǎng)格尺寸為1，共計(jì)10x10x=1000個(gè)單元。總的計(jì)算時(shí)長(zhǎng)為1s，這只固定增量步長(zhǎng)為0.01s，所以總增量步數(shù)為100。
（）計(jì)算結(jié)果
以上面網(wǎng)格圖中中間角點(diǎn)為例，提取加載向位移時(shí)程曲線如下圖：
0.89s時(shí)刻x向應(yīng)力云圖為：
注：左一為ABAQUS計(jì)算結(jié)果，中間為umat實(shí)現(xiàn)線彈性各向同性本構(gòu)計(jì)算結(jié)果，右一為UEL+UMAT計(jì)算結(jié)果，后面的云圖也是這樣排布。
0.89s時(shí)刻y向應(yīng)力云圖為：
寫(xiě)文字好麻煩，不寫(xiě)了！放個(gè)動(dòng)圖算了，有興趣的移步附件。
展開(kāi) ABAQUS UEL二次開(kāi)發(fā)（動(dòng)力隱式C3D8自定義單元）
概述：開(kāi)發(fā)了適用于靜力通用、頻率分析和動(dòng)力隱式（固定增量步長(zhǎng)和自適應(yīng)增量步長(zhǎng)均可）的三維八節(jié)點(diǎn)線性UEL，即ABAQUS自帶的C3D8單元，該UEL考慮了B-BAR修正，避免體積鎖死。采用編寫(xiě)的UEL，分別設(shè)置了靜力通用分析步、頻率分析和動(dòng)力隱式分析步，將計(jì)算結(jié)果與ABAQUS對(duì)比，位移、速度和加速度與ABAQUS均保持一致，說(shuō)明該UEL復(fù)現(xiàn)了一小部分C3D8單元的計(jì)算功能。
（一）模型信息
懸臂梁尺寸：10x10x100，密度1e10，密度200，泊松比0.25。不設(shè)置單位，純驗(yàn)證。
網(wǎng)格如下圖：，每個(gè)單元尺寸為：2x2x5。
（二）靜力分析
邊界條件如下圖：
荷載大小為：1e6，采用固定增量步長(zhǎng)，計(jì)算總時(shí)長(zhǎng)為10（靜力計(jì)算中計(jì)算時(shí)長(zhǎng)無(wú)意義，僅為驗(yàn)證設(shè)置），增量步長(zhǎng)為0.01，總增量步數(shù)為1000。
總位移云圖如下圖：
加載向（U2、Y向）位移云圖如下圖：
梁向（U3，z向）位移云圖如下圖：
U1，x向位移云圖如下圖：
懸臂端角點(diǎn)加載向位移-荷載歷程如下圖：
注意這里的時(shí)間并沒(méi)有物理意義，在靜力通用分析步中僅僅是增量步的計(jì)算意義。
（三）頻率分析
邊界條件，梁一端固定。計(jì)算前100階頻率信息。
展開(kāi) 關(guān)于Abaqus UEL中RHS數(shù)組長(zhǎng)度大于單元總自由度數(shù)的一些解釋
在Abaqus中，用戶自定義元素子程序（UEL）的開(kāi)發(fā)需要遵循一些特定的規(guī)則和約定。其中一個(gè)關(guān)鍵約定是關(guān)于子程序中的RHS（右手邊）向量的維數(shù)，我最近在嘗試用UEL做一些二次開(kāi)發(fā)，也發(fā)現(xiàn)了RHS向量的維數(shù)比單元的總自由度數(shù)多了4個(gè)這一現(xiàn)象，結(jié)合在站內(nèi)一些同行的猜測(cè)，我認(rèn)為這可能是由于Abaqus的內(nèi)部工作方式所導(dǎo)致的。
首先，有同行懷疑是因?yàn)橛昧怂墓?jié)點(diǎn)單元，所以多了四個(gè)，我開(kāi)發(fā)的單元是12個(gè)節(jié)點(diǎn)的，依然多了四個(gè)維度，因此排除是單元內(nèi)節(jié)點(diǎn)個(gè)數(shù)導(dǎo)致的。
在Abaqus中，RHS向量的維數(shù)實(shí)際上包括了除了單元的位移自由度外的其他項(xiàng)，這些項(xiàng)用于處理多種情況，例如：
體積力和表面力的計(jì)算：RHS向量可能包括用于計(jì)算體積力和表面力的額外自由度。這些自由度用于存儲(chǔ)單元內(nèi)的體積力和表面力的貢獻(xiàn)。
約束和邊界條件：Abaqus可能需要額外的自由度來(lái)處理約束條件和邊界條件，以確保數(shù)值穩(wěn)定性和正確的求解結(jié)果。
內(nèi)部狀態(tài)變量：某些材料模型和非線性分析可能需要存儲(chǔ)和更新一些內(nèi)部狀態(tài)變量，這些變量也可以占用RHS向量中的額外位置。
因此，RHS向量的維數(shù)不僅僅包括單元的位移自由度，還包括其他與分析和模型特性相關(guān)的項(xiàng)。這是Abaqus設(shè)計(jì)的一部分，旨在確保通用性和可擴(kuò)展性，以處理各種復(fù)雜的問(wèn)題。其中我認(rèn)為可擴(kuò)展性是一項(xiàng)比較重要的應(yīng)用，我目前所做的工作可能會(huì)利用到這一點(diǎn)，等有結(jié)果了會(huì)繼續(xù)更新。
展開(kāi) 
基于ABAQUS的UEL子程序定義4節(jié)點(diǎn)平面應(yīng)變等參單元的剛度問(wèn)題
摘要：
采用基于ABAQUS的UEL子程序開(kāi)發(fā)4節(jié)點(diǎn)平面應(yīng)變等參單元，采用雙線性形函數(shù)，4點(diǎn)高斯積分，本構(gòu)關(guān)系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對(duì)比ANSYS的單元剛度矩陣，結(jié)果顯示兩者也不相同。問(wèn)題出在哪里呢？本文檔將對(duì)此問(wèn)題進(jìn)行回答。
本文可以作為ABAQUS高級(jí)子程序UEL的入門(mén)級(jí)教程，做UEL的應(yīng)該關(guān)注下！
基于ABAQUS的UEL子程序定義4節(jié)點(diǎn)平面應(yīng)變等參單元的剛度問(wèn)題（技術(shù)鄰 藍(lán)牙）.pdf
展開(kāi) ABAQUS用戶定義單元UEL與VUEL從入門(mén)到放棄系列3————講一講的RHS和NBLOCK ￥50
我們?cè)賮?lái)看VUEL中關(guān)于RHS的說(shuō)明
依舊定義為單元對(duì)系統(tǒng)方程右端的貢獻(xiàn)，但此時(shí)不再是外力-內(nèi)力，而是定義為包含單元的內(nèi)力或是由分布式載荷計(jì)算而來(lái)的外部載荷，故在沒(méi)有分布式載荷時(shí)，RHS=內(nèi)力KU（不考慮非線性）；因此可再次猜測(cè)，ABAQUS UEL和VUEL所說(shuō)的外載荷就是分布式載荷，其它INP中定義的外載荷不做考慮。
接下來(lái)討論一個(gè)VUEL中特有的數(shù)組NBLOCK，譯為塊，也正是這個(gè)數(shù)組標(biāo)志著VUEL與UEL運(yùn)行的不同；眾所周知，ABAQUS調(diào)用UEL計(jì)算時(shí)，是一次進(jìn)一個(gè)單元，然后計(jì)算這個(gè)單元的RHS、AMATRX等數(shù)組并回傳給ABAQUS，然后再進(jìn)第二個(gè)單元，如此進(jìn)行...；
而ABAQUS調(diào)用VUEL計(jì)算時(shí)，是一次進(jìn)入一批單元，這個(gè)一批單元個(gè)數(shù)即為NBLOCK的大小（之前一位名為Xujianqing的作者曾經(jīng)發(fā)過(guò)一個(gè)帖子說(shuō)明說(shuō)明過(guò)NBLOCK最大值為136，但我后來(lái)測(cè)試發(fā)現(xiàn)，在低版本如6.14中，最大NBLCOK=136，在高版本如2020中，最大NBLOCK=144）。
這是個(gè)什么意思呢？
展開(kāi) 最新項(xiàng)目需求丨在線派單，先到先得！
點(diǎn)擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8281
項(xiàng)目五：
預(yù)算范圍：3000-5000搶單中
發(fā)布時(shí)間：2023-02-28
工期：3.10號(hào)
使用的軟件： abaqus中UEL單元開(kāi)發(fā)
需求描述：
abaqus uel子結(jié)構(gòu)實(shí)現(xiàn)混凝土水管冷卻溫度場(chǎng)計(jì)算，程序功能包括子結(jié)構(gòu)的凝聚、水溫迭代計(jì)算等
點(diǎn)擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8284
項(xiàng)目六：
預(yù)算范圍：100搶單中
發(fā)布時(shí)間：2023-02-27
使用的軟件： ansys
需求描述：
*CONSTRAINED_BEAM_IN_SOLID模型怎么施加粘結(jié)滑移關(guān)系
點(diǎn)擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8279
項(xiàng)目七：
預(yù)算范圍：當(dāng)面咨詢搶單中
發(fā)布時(shí)間：2023-02-26
工期：一天
使用的軟件： Abaqus聲學(xué)穩(wěn)態(tài)動(dòng)力學(xué)分析，聲固耦合
需求描述：
我有一個(gè)笛子的內(nèi)部空氣柱模型，要求對(duì)笛子內(nèi)部空氣柱進(jìn)行諧振分析(不用直接法)，我不知道內(nèi)部參數(shù)應(yīng)該如何定義，主要目的一是得出諧振，二是希望可以告知子空間法以及模態(tài)法的載荷定義以及阻尼的相關(guān)設(shè)置。
點(diǎn)擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8274
項(xiàng)目八：
預(yù)算范圍：3000搶單中
發(fā)布時(shí)間：2023-02-26
工期：可探討
使用的軟件： Fluent
需求描述：
使用UDF實(shí)現(xiàn)顆粒擴(kuò)散泳熱泳仿真，出不同粒徑粒子擴(kuò)散泳沉積量隨時(shí)間變化的曲線圖。
展開(kāi) 【JY】 ABAQUS子程序UEL的有限元原理與應(yīng)用
本文對(duì)ABAQUS/Standard中UEL子程序的有限元原理及編寫(xiě)原則進(jìn)行初步梳理，并編寫(xiě)了平面三角形單元靜力分析的UEL子程序，以加深對(duì)UEL的理解與認(rèn)識(shí)。
【ABAQUS的UEL單元之有限元基本原理簡(jiǎn)述】
當(dāng)用戶需要用ABAQUS/Standard單元庫(kù)中沒(méi)有的單元進(jìn)行分析時(shí)，可通過(guò)ABAQUS提供的UEL子程序接口進(jìn)行二次開(kāi)發(fā)，編寫(xiě)適用于特定分析的單元。下面以平面三角形單元為例，給出編寫(xiě)UEL單元的基本流程。
系統(tǒng)的勢(shì)能Π由下式給出
其中為單元的應(yīng)變能。
展開(kāi) 