abaqus 膜單元定義的案例
ABAQUS任意單元表面加入膜單元或加入復合材料纖維層
以上內容來自360百科
本期是教大家如何在ABAQUS有限元模型中在任意實體單元表面加入殼單元作為纖維增強材料來模擬復合材料：
孔眼壁上的膜單元來模擬壁面加固材料
內加入纖維增強材料
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
Abaqus運用Vusdfld實現自定義單元刪除
Abaqus運用Vusdfld實現自定義單元刪除
———公眾號‘CAE仿真實驗室’出品
Abaqus提供了幾種單元刪除的方法，如下表所示，單元刪除的本質是讓單元喪失承受應力的能力，使其對模型剛度沒有貢獻。在顯式分析中，可以利用*depvar，delete實現單元刪除，該方法是通過給單元設定一個狀態變量來控制單元是否刪除，當狀態變量為0時，單元刪除，狀態變量為1時，單元不刪除。今天就來介紹一些通過VUSDFLD子程序來自定義單元刪除準則的例子
1、例子
單軸拉伸樣條，左端固定，右端給定一個速度拉伸，假設當應變超過0.02時，單元就失效，刪除單元。為了讓中間某些單元更快達到失效應變，事先在式樣中間刪除幾個單元，使中間應變更大。
定義材料模型時，設置Depvar和控制element deletion的狀態變量并設置User Defined Field.Vusdfld關鍵代碼如下，關于vusdfld用法參考之前有關usdfld的帖子。
2、結果
應變結果如下，開始時應變均勻增大，達到失效應變后，中間單元應變極劇增大，說明該單元基本沒有‘剛度’。
應變云圖和應變曲線
單元刪除后，單元仍然在，但是單元的應力為0，右端拉力在失效前均勻增大，失效后，急劇下降，之后結果波動，此時結果已經沒有意義~
應力云圖和應力曲線
3、結論
由上面結果可以看出，基本實現了單元刪除的功能，但是該例子較為簡單，用于復雜模型還沒有驗證。
展開 一種新思路用于實現ABAQUS用戶自定義單元
</span></blockquote><p><br></p><blockquote><strong style="color: rgb(25, 27, 31);">這種自定義單元完美解決了ABAUS傳統的用于自定義單元無法完成可視化的問題</strong> <span style="color: rgb(25, 27, 31);">。搞過UEL二次開發的人都被自定義單元的可視化問題困擾，要么是將數據導出到第三方軟件進行處理，要么是采用UMAT套一層單元進行可視化，這兩種方法都需要大量的編程工作。</span></blockquote><p><br></p><blockquote><strong>這種自定義單元完美的解決了ABAQUS傳統的用戶自定義單元無法施加面力、體力等復雜力的短板</strong>。</blockquote><p><br></p><blockquote>更多的功能還在探索之中...
展開 ABAQUS用戶自定義單元UEL學習資料 ￥49
ABAQUS子程序UEL
abaqus最簡單的uel子程序自定義單元例子 ￥3
通過abaqus的uel子程序定義了2節點單元的剛度矩陣，直接運行分析可以得到節點應力、位移等參數，可以快速了解uel的組織架構。
運行方法：abaqus創建job——來源——點擊輸入文件——加載inp文件——編輯作業——通用模塊——加載用戶子程序for文件——運行
淺談abaqus針對不同單元類型定義初始溫度場
在進行熱-應力分析時，初始溫度場的定義為最常見的。針對不同的單元類型（Solid單元、Shell單元、Beam單元），Abaqus提供了多種不同的定義初始溫度場的方法，可以根據實際情況靈活的選擇不同的定義方式，從而更加精確的實現仿真分析。下面簡單的介紹一下在Abaqus中以上三種單元定義初始溫度場的方法。
l
Solid單元初始溫度場定義
l
Shell單元初始溫度場定義
l
Beam單元的初始溫度場定義
這三部分單元的初始溫度場定義詳見附件：
淺談abaqus針對不同單元類型的初始溫度場定義.pdf
展開 ABAQUS用戶定義單元UEL與VUEL從入門到放棄系列2 ￥50
本構關系的張量表達式為：
其中，廣義應變的有限元格式可表示為：
其中，電場強度E為負的電勢梯度：
則廣義應變列陣記為：
那么，單元的剛度矩陣可以表示為：
其中Kuu為C3D8原本的剛度矩陣，Kuv與Kvu為壓電耦合剛度矩陣，Kvv為電場的廣義剛度矩陣。
ABAQUS中的C3D8和C3D8E都是做了一些剛度修正的，比如C3D8為了防止單元自鎖，采用了B-Bar方法，得出的剛度矩陣是介于C3D8和C3D8R之間的值，同樣的C3D8E也有一些類似的修正，以下我將提供一個不包含修正的版本，對ABAQUS剛度修正方法感興趣的朋友可以去拿去跟ABAQUS CAE對比。
展開 ABAQUS UEL二次開發（動力隱式C3D8自定義單元）
概述：開發了適用于靜力通用、頻率分析和動力隱式（固定增量步長和自適應增量步長均可）的三維八節點線性UEL，即ABAQUS自帶的C3D8單元，該UEL考慮了B-BAR修正，避免體積鎖死。采用編寫的UEL，分別設置了靜力通用分析步、頻率分析和動力隱式分析步，將計算結果與ABAQUS對比，位移、速度和加速度與ABAQUS均保持一致，說明該UEL復現了一小部分C3D8單元的計算功能。
（一）模型信息
懸臂梁尺寸：10x10x100，密度1e10，密度200，泊松比0.25。不設置單位，純驗證。
網格如下圖：，每個單元尺寸為：2x2x5。
（二）靜力分析
邊界條件如下圖：
荷載大小為：1e6，采用固定增量步長，計算總時長為10（靜力計算中計算時長無意義，僅為驗證設置），增量步長為0.01，總增量步數為1000。
總位移云圖如下圖：
加載向（U2、Y向）位移云圖如下圖：
梁向（U3，z向）位移云圖如下圖：
U1，x向位移云圖如下圖：
懸臂端角點加載向位移-荷載歷程如下圖：
注意這里的時間并沒有物理意義，在靜力通用分析步中僅僅是增量步的計算意義。
（三）頻率分析
邊界條件，梁一端固定。計算前100階頻率信息。
展開 ABAQUS用戶定義單元UEL與VUEL從入門到放棄系列1 ￥20
最初來技術鄰，講實話其實就是來學UEL的，后來發現沒帖子，然后就各種找論壇，simweABAQUS版塊的帖子二零零幾年的帖子都被我挖干凈了，做UEL的人確實比做UMAT的少多了，但技術鄰這個平臺依舊讓我受益良多，諸如藍牙老師，isolver團隊snowwave02等大佬，他們的分享，甚至很多是免費分享確實高屋建瓴，知識的學習必是知其所以然，我想只有這樣的學術氛圍越來越重，這樣樂于分享的大佬越來越多，中國的自研工業軟件才能走向世界，我們的工業4.0才有可能實現。
學于此，也該傳于此，作為一只科研狗，我想我比誰都清楚什么叫從“入門”到“放棄”，未知的知識就是這樣，不懂的時候它比金子都貴，等你學會了它好像又像水一樣普通，亦如水一樣珍貴。
之后也會逐漸在技術鄰更一些帖子，可能也會有視頻吧，內容大概會圍繞有限元理論，用戶子程序開發，ABAQUS python開發，XFEM，lamb波的SHM等，因為不知道想學UEL\VUEL的人群多大，也不知道大家都想了解些什么，不知道應該從何說起，可能有些想學UEL的小白基礎很差連有限元理論都沒學全，再加上我確實是個懶癌晚期，其實很多東西都做好了，懶的發，也希望能治療一下自己的懶癌吧。
今天的正事，發兩個UEL的Fortran程序和對應的inp算例，我發現好像很多人不會寫UEL的inp，下次有機會給你們講講吧（懶癌警告）。
一個是平面應力\平面應變單元，2d4node程序
一個是3D stress單元，C3D8程序
展開 基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 ABAQUS用戶定義單元UEL與VUEL從入門到放棄系列3————講一講的RHS和NBLOCK ￥50
我們再來看VUEL中關于RHS的說明
依舊定義為單元對系統方程右端的貢獻，但此時不再是外力-內力，而是定義為包含單元的內力或是由分布式載荷計算而來的外部載荷，故在沒有分布式載荷時，RHS=內力KU（不考慮非線性）；因此可再次猜測，ABAQUS UEL和VUEL所說的外載荷就是分布式載荷，其它INP中定義的外載荷不做考慮。
接下來討論一個VUEL中特有的數組NBLOCK，譯為塊，也正是這個數組標志著VUEL與UEL運行的不同；眾所周知，ABAQUS調用UEL計算時，是一次進一個單元，然后計算這個單元的RHS、AMATRX等數組并回傳給ABAQUS，然后再進第二個單元，如此進行...；
而ABAQUS調用VUEL計算時，是一次進入一批單元，這個一批單元個數即為NBLOCK的大小（之前一位名為Xujianqing的作者曾經發過一個帖子說明說明過NBLOCK最大值為136，但我后來測試發現，在低版本如6.14中，最大NBLCOK=136，在高版本如2020中，最大NBLOCK=144）。
這是個什么意思呢？
展開 
