abaqus單元VUEL的案例
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與VUEL從入門到放棄系列2 ￥50
最初來技術鄰，講實話其實就是來學UEL的，后來發現沒帖子，然后就各種找論壇，simweABAQUS版塊的帖子二零零幾年的帖子都被我挖干凈了，做UEL的人確實比做UMAT的少多了，但技術鄰這個平臺依舊讓我受益良多，諸如藍牙老師，isolver團隊snowwave02等大佬，他們的分享，甚至很多是免費分享確實高屋建瓴，知識的學習必是知其所以然，我想只有這樣的學術氛圍越來越重，這樣樂于分享的大佬越來越多，中國的自研工業軟件才能走向世界，我們的工業4.0才有可能實現。
學于此，也該傳于此，作為一只科研狗，我想我比誰都清楚什么叫從“入門”到“放棄”，未知的知識就是這樣，不懂的時候它比金子都貴，等你學會了它好像又像水一樣普通，亦如水一樣珍貴。
之后也會逐漸在技術鄰更一些帖子，可能也會有視頻吧，內容大概會圍繞有限元理論，用戶子程序開發，ABAQUS python開發，XFEM，lamb波的SHM等，因為不知道想學UEL\VUEL的人群多大，也不知道大家都想了解些什么，不知道應該從何說起，可能有些想學UEL的小白基礎很差連有限元理論都沒學全，再加上我確實是個懶癌晚期，其實很多東西都做好了，懶的發，也希望能治療一下自己的懶癌吧。
今日份正事，給大家簡單講一講ABAQUS中的壓電耦合單元C3D8E。
那么什么叫壓電耦合單元呢，簡單的說就是你給它加載電壓（電勢的差，ABAQUS中為位移加載），那么就會引起單元力學場的變化，比如位移、應力、應變等等；同樣的你給它加載力或者位移，亦會引起單元電場的變化。
單元壓電耦合場的廣義本構方程表示如下：
即為：廣義應力=廣義彈性矩陣·廣義應變。其中廣義力中D為電位移，廣義應變中E為電場強度。
在廣義彈性矩陣中，C矩陣為力學場的彈性矩陣，e矩陣為壓電常數矩陣，右下角為介電常數矩陣。
展開 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）。
這是個什么意思呢？
展開 