自定義單元開發的案例
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階頻率信息。
展開 ANSYS梁單元自定義截面
ANSYS梁單元自定義截面
梁單元作為一種簡單且高效的計算單元，在結構分析尤其是建筑結構中得到廣泛的應用。使用梁單元可以避免將結構中梁柱全部轉換為實體單元，從而降低了計算量，且梁單元結構形式簡單，求解精度也相對較高。在ANSYS中，梁單元基本上可以分為線性單元和二次單元，二者之間計算理論不同，經典的二次單元即BEAM189單元的積分點如下圖所示：
在ANSYS中可以為BEAM單元定義截面，其中大部分經典的截面形式都包含在ANSYS的截面庫中，但是經典的梁單元計算時截面方向分為四個單元，這對于一般計算來說是足夠的，但如果需要仔細分析截面方向的內力，可能就略顯的粗糙了。除此之外，鋼管混凝土、組合梁之類也都是異形梁截面，此時標準截面庫中的數據也沒什么用。針對這個問題存在兩種解決方式，一種是使用ASEC自定義截面參數，這個命令不管截面如何，只需要給出截面相關的信息即可，截面的信息輸入如下圖所示：
至于這些截面的參數可以使用簡單的截面計算工具得到，如果是鋼筋混凝土梁這種比較復雜的復合梁，那么需要使用Xtract之類的截面有限元軟件進行計算。將截面信息填入。采用ASEC的截面輸入方式計算效率高，截面信息準確的話，精度也不差，但缺點是不能輸出截面積分點和柵點的數據。
另一種方式就是自定義截面，其基本思路如下：
1.設定MESH200單元，建立截面幾何形狀；
2.用MESH200單元劃分截面，并保存截面數據；
3.建立計算幾何模型，讀取截面數據；
4.賦予模型截面，施加邊界條件計算；
5.后處理。
展開 ANSYS梁單元自定義截面
梁單元作為一種簡單且高效的計算單元，在結構分析尤其是建筑結構中得到廣泛的應用。使用梁單元可以避免將結構中梁柱全部轉換為實體單元，從而降低了計算量，且梁單元結構形式簡單，求解精度也相對較高。在ANSYS中，梁單元基本上可以分為線性單元和二次單元，二者之間計算理論不同，經典的二次單元即BEAM189單元的積分點如下圖所示：
在ANSYS中可以為BEAM單元定義截面，其中大部分經典的截面形式都包含在ANSYS的截面庫中，但是經典的梁單元計算時截面方向分為四個單元，這對于一般計算來說是足夠的，但如果需要仔細分析截面方向的內力，可能就略顯的粗糙了。除此之外，鋼管混凝土、組合梁之類也都是異形梁截面，此時標準截面庫中的數據也沒什么用。針對這個問題存在兩種解決方式，一種是使用ASEC自定義截面參數，這個命令不管截面如何，只需要給出截面相關的信息即可，截面的信息輸入如下圖所示：
至于這些截面的參數可以使用簡單的截面計算工具得到，如果是鋼筋混凝土梁這種比較復雜的復合梁，那么需要使用Xtract之類的截面有限元軟件進行計算。將截面信息填入。采用ASEC的截面輸入方式計算效率高，截面信息準確的話，精度也不差，但缺點是不能輸出截面積分點和柵點的數據。
另一種方式就是自定義截面，其基本思路如下：
1.設定MESH200單元，建立截面幾何形狀；
2.用MESH200單元劃分截面，并保存截面數據；
3.建立計算幾何模型，讀取截面數據；
4.賦予模型截面，施加邊界條件計算；
5.后處理。
展開 ANSYS梁單元自定義截面
ANSYS梁單元自定義截面
梁單元作為一種簡單且高效的計算單元，在結構分析尤其是建筑結構中得到廣泛的應用。使用梁單元可以避免將結構中梁柱全部轉換為實體單元，從而降低了計算量，且梁單元結構形式簡單，求解精度也相對較高。在ANSYS中，梁單元基本上可以分為線性單元和二次單元，二者之間計算理論不同，經典的二次單元即BEAM189單元的積分點如下圖所示：
在ANSYS中可以為BEAM單元定義截面，其中大部分經典的截面形式都包含在ANSYS的截面庫中，但是經典的梁單元計算時截面方向分為四個單元，這對于一般計算來說是足夠的，但如果需要仔細分析截面方向的內力，可能就略顯的粗糙了。除此之外，鋼管混凝土、組合梁之類也都是異形梁截面，此時標準截面庫中的數據也沒什么用。針對這個問題存在兩種解決方式，一種是使用ASEC自定義截面參數，這個命令不管截面如何，只需要給出截面相關的信息即可，截面的信息輸入如下圖所示：
至于這些截面的參數可以使用簡單的截面計算工具得到，如果是鋼筋混凝土梁這種比較復雜的復合梁，那么需要使用Xtract之類的截面有限元軟件進行計算。將截面信息填入。采用ASEC的截面輸入方式計算效率高，截面信息準確的話，精度也不差，但缺點是不能輸出截面積分點和柵點的數據。
另一種方式就是自定義截面，其基本思路如下：
1.設定MESH200單元，建立截面幾何形狀；
2.用MESH200單元劃分截面，并保存截面數據；
3.建立計算幾何模型，讀取截面數據；
4.賦予模型截面，施加邊界條件計算；
5.后處理。
展開 
hypermesh二次開發之界面自定義按鈕的添加 ￥30
<p>hypermesh擁有強大的二次開發功能，用戶可以利用Tcl / Tk語言進行二次開發，涉及自動化腳本、界面按鈕添加、彈窗等，本次主要講述通過在Utility菜單的User頁面添加自定義按鈕，按鈕附加Tcl / Tk腳本，直接進行用戶所需的特殊功能。如下圖所示：</p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202503/attachment/c7171445acb943e7a194e71d2310faac.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202503/attachment/c7171445acb943e7a194e71d2310faac.png" data-mobile-src="https://img.jishulink.com/202503/attachment/c7171445acb943e7a194e71d2310faac.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202503/attachment/c7171445acb943e7a194e71d2310faac.png?
展開 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用戶自定義單元UEL學習資料 ￥49
ABAQUS子程序UEL
一種新思路用于實現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子程序自定義單元例子 ￥3
通過abaqus的uel子程序定義了2節點單元的剛度矩陣，直接運行分析可以得到節點應力、位移等參數，可以快速了解uel的組織架構。
運行方法：abaqus創建job——來源——點擊輸入文件——加載inp文件——編輯作業——通用模塊——加載用戶子程序for文件——運行
APDL自定義截面打開單元形狀后總是對不齊，怎么辦？
ANSYS中的自定義截面功能為用戶定義復雜截面提供了方便，然而部分同學在處理復雜截面時，有遇得到過明明截面外框尺寸都一樣，而網格劃分后打開單元形狀發現截面出現偏離對不齊的情況。
以下面的案例為例，此處自定義了兩個復雜箱梁截面，均采用在CAD建立好面域，導入ANSYS中劃分網格，然后在采用secwrite命令自定義截面。
兩個截面的外觀尺寸完全一致，唯一不同的是內部空心 形狀不一致！
截面一形式如下：
截面二形式如下：
導入ANSYS，劃分網格后，并打開單元形狀，截面如下：
從截圖中可看到，兩者截面交匯處有明顯的截面錯位，放大如下所示：
有的同學會認為兩者截面外部尺寸完全一樣，理論上應該重合才是，可為啥還有這種情況出現呢？
原因如下：
在ANSYS的梁單元截面定義中，ANSYS默認梁單元中心線(單元坐標系X軸)位于截面形心處，當用戶用secplot繪制截面形狀時所出現的Centroid Y和Centroid Z即為截面形心相當于截面定義原點的坐標。如果兩個截面的這兩個坐標不一致，就會出現類似上圖中的錯位現象！
單元坐標系X軸以截面形心位置為準
那么截面原點位于哪個地方呢？自帶截面庫對于雙邊對稱的情況（例如矩形，圓形），原點一般位于截面對稱中心處，而對于單邊對稱的情況（例如槽型截面），原點一般位于左下角原點。如果用戶采用了自定義截面，則原點位于用戶自定義截面時的坐標原點（0,0）處。
例如此處我們在cad中繪制截面時，坐標（0,0）位于頂部中心處，那么這個點就是我們截面的原點。
此處注意坐標系的轉換，雖然我們在cad中是X、Y坐標，但換到ANSYS截面定義中，x坐標代表單元坐標系Y軸，y坐標系代表單元坐標系Z軸。
展開 『原創』請問怎樣自定義殼單元的截面
如題
梁單元截面偏置(用戶自定義位置)計算方法
截面定義
sectype,1,beam,rect
secdata,600,250
secoffset,user,-300,-125
.
.
.
!================
!結構幾何建模
blc4,,,3600,3600
k,5,0,0,-3600
l,5,1
lsel,s,,,5
lgen,2,all,,,3600
lgen,2,all,,,,3600
nummrg,all
!================
!截面屬性賦值
lsel,s,,,1
latt,1,,1,,,4,1
.
.
.
!================
!劃分單元及截面顯示
allsel,all
lesize,all,,,3
lmesh,all
allsel,all
/eshape,1
展開 梁單元截面偏置(用戶自定義位置)計算方法（workbench版本）
經典版見水哥的帖子
鏈接為：
梁單元截面偏置(用戶自定義位置)計算方法
workbench的方法
主要是梁的定義
未打開截面的幾何模型
打開截面的幾何模型具體的對齊方式
網格劃分：
workbench的優點：
操作方便，簡單易懂
workbench的缺點
無法像經典能夠對截面的長寬邊的劃分數目進行控制，只能是1.我找了好久，確實沒發現，貌似workbench計算出面積，慣性矩等參量，直接代入
木網殼雙節點考慮轉動受軸力影響的自定義連接單元-UCEL ￥200
（1）基于ABAQUS UEL子程序接口創建了此單元；
（2）基于WE-P理論分析模型編寫的子程序關系。
『原創』ANSYS中殼單元截面如果能夠自定義該是有限元技術中的一個難點突破!
殼單元是工程實際應用中一種重要的單元形式,能夠解決非常多的實際問題!比如壓力容器,橋梁分析,鋼結構分析,復合材料,汽車,船舶等等!
然而在多年的有限元工程應用中,有一個問題一直都困擾著我,問題描述如下:有一大類薄板結構,其截面是不規則的,如果按照均勻薄板結構來算顯然會有較大出入;若按照梁殼結合,工作量將是非常大,且未必能夠很好的解決!
某突發奇想,如果有限元中能象解決梁截面一樣,在分析中也可以自定義殼截面那改有多好啊!
這個問題我在仿真互動論壇中也發過貼子,在這里希望繼續和大家探討,多交流,看是否還有什么更好的解決辦法!
展開 