UMESHMOTION子程序應用的案例
Abaqus子程序Umeshmotion案例 附umeshmotion翻譯下載
下載地址：umeshmotion翻譯
umeshmotion子程序
PDF
UMESHMOTION子程序，磨損分析 ￥10
UMESHMOTION子程序可用于磨損和燒蝕，網上關于UMESHMOTION子程序的資料較少，筆者找了一些資料整理如下。由于本人并非專業人士，無法對質量進行過多評價，希望能對需要的人員有所幫助：
1. ABAQUS用戶手冊——輪胎磨損經典例子，考慮輪胎紋理邊角的磨損，涉及到不同磨損區域的判定和局部磨損方向的定義。（本人對其進行了逐行的詳細解釋）
http://wufengyun.com:888/books/exa/default.htm
2. ABAQUS子程序講義，內含該程序的基本概念解釋。
https://fdocuments.in/document/user-sub.html
3. 一本介紹磨損的書籍，內含一些ABAQUS磨損案例，可通過這些經典案例的曲線形狀判斷子程序結果準確與否。
書名：Numerical Modeling of the Effect of Fretting Wear on Fretting Fatigue
http://eprints.nottingham.ac.uk/10681/1/ThesisFinal.pdf
4. 一個外文資料，基于Archard理論和實驗進行對比，并對相關理論公式進行說明。
https://www.scirp.org/html/1-2190164_96936.htm#f3
5. 百度文庫中的一個例子，簡單二維磨損案例有配圖，沒有基于磨損理論進行子程序編寫，也沒有說明磨損過程。
https://wenku.baidu.com/view/305572e86aec0975f46527d3240c844769eaa0c4.html
6. 一個簡單的材料燒蝕帖子及啟發（附案例）。
https://www.smart3de.com/forum.php?
展開 ABAQUS用戶自定義子程序UMESHMOTION仿真案例講解
[圖片]
Abaqus輪-軌接觸摩擦磨損(UMESHMOTION子程序)仿真案例講解 ￥600
[圖片]
abaqus通過umeshmotion子程序模擬沖壓過程中模具的磨損行為
abaqus通過umeshmotion子程序模擬沖壓過程中模具的磨損行為
沖壓過程中，模具磨損是最常見的一種現象。模具磨損不僅會影響磨具的壽命，也會影響沖壓件的成型質量。因此沖壓磨具的磨損分析具有重要意義。
目前磨損分析中使用最為廣泛的理論為Archard理論。Archard模型的一般公式為
式中：dV為磨損體積，dP為接觸面的法向壓力，dL為切向相對滑移，H為模具硬度，K為磨損因子。根據式（1）可以得到模具磨損深度的計算公式
Abaqus中可以通過umeshmotion子程序進行結構的磨損分析。本文編寫了基于archard模型的umeshmotion子程序，并結合model change對沖壓過程中磨具的磨損行為進行了分析,有限元模型如下。
模擬結果如下，為了簡便，這里只模擬了三次沖壓的磨損情況。
磨損前后的輪廓對比
磨損量和沖壓次數的關系
展開 ABAQUS子程序umeshmotion報錯，磨損深度能超過網格單元高度嗎
1 elements in the adaptive mesh domain are distorting so much that
they turn inside out. The elements have been identified in element
set WarnElemAdaptMeshDistortStep3Inc69.
ADAPTIVE MESH CONSTRAINTS CANNOT BE APPLIED IN THIS INCREMENT.
CONSIDER RERUNNING THE ANALYSIS WITH SMALLER TIME INCREMENTATION IN
THIS STEP.
展開 vpsc-vumat子程序的開發與應用文章推薦
結合作者的理論（尤其是分段線性化和應力驅動的求解思路）我們可以把獨立的vpsc子程序編寫進abaqus里面，為了避免復雜的雅可比推導，以及適用各種復雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序
以下展示一個使用vpsc-鎂合金本構模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數取自vpsc官方案例,為了減少計算時間使用高應變率進行計算，應力高于vpsc模擬）：
等效塑性應變：
第一個晶粒的累計剪切滑移：
發生孿晶次數;
變形后的形狀演化：
ABAQUS用戶子程序USDFLD應用實例
一、引言
為滿足用戶多樣化的分析需求，ABAQUS提供了眾多子程序接口。其中USDFLD可用于復雜材料模型的建模分析，例如粉末冶金、損傷分析等。本文以某藥物粉末的模壓為例，介紹了ABAQUS用戶子程序USDFLD的使用方法。
二、用戶子程序USDFLD簡介
1.USDFLD即Use Defined Field，通過將材料屬性定義為場變量fi的函數，來實現復雜材料行為的建模
①ABAQUS/Standard中的大多數材料屬性都可定義為場變量的函數。
②USDFLD允許用戶在單元的每個積分點處定義fi。
設計仿真 | Adams線性化分析子程序應用
02
線性化分析子程序功能
新版添加的c_get_linear_mat_states函數，可以使得工程師進行更加方便的用戶子程序綜合應用。之前，為了獲取系統的狀態空間矩陣，比如標準的ABCD矩陣形式，只能通過Linear/statemat命令實現，沒有其它的功能可以獲取該類數據。因此，當工程師需要在自己編寫的用戶子程序進行狀態空間矩陣數據調用時，只能通過先將這些數據求解并存儲于指定文件中，然后再通過編程讀取這些數據，非常的不方便。
有了c_get_linear_mat_states功能后，工程師可以通過子程序CONSUB進行狀態空間矩陣的求解后，直接進行其數據的調用，也就是狀態空間矩陣求解以及狀態向量獲取可以直接在用戶子程序中完成，不再需要額外的輸入輸出工作。
上圖仍舊是單自由度彈簧振子模型，只不過在本模型中，剛度和阻尼會隨時間變化，并且將會在0s，1s，5s通過子程序CONSUB分別對系統進行線性化處理，當然在子程序中有c_get_linear_mat_states函數的調用。另外，模型上施加的外部載荷、彈簧力都是通過Gforce實現的，通過函數的綜合應用實現更靈活的設置。
為了實現上述所要求的仿真工況，必須借助腳本命令來完成，在初始階段首先對模型進行初始化求解，接著完成靜平衡計算，然后才進行第一次線性化求解，進而進行動態求解，在1秒完成時，進行第二次線性化求解，以此類推。
展開 設計仿真 | Adams線性化分析子程序應用
02
線性化分析子程序功能
新版添加的c_get_linear_mat_states函數，可以使得工程師進行更加方便的用戶子程序綜合應用。之前，為了獲取系統的狀態空間矩陣，比如標準的ABCD矩陣形式，只能通過Linear/statemat命令實現，沒有其它的功能可以獲取該類數據。因此，當工程師需要在自己編寫的用戶子程序進行狀態空間矩陣數據調用時，只能通過先將這些數據求解并存儲于指定文件中，然后再通過編程讀取這些數據，非常的不方便。
有了c_get_linear_mat_states功能后，工程師可以通過子程序CONSUB進行狀態空間矩陣的求解后，直接進行其數據的調用，也就是狀態空間矩陣求解以及狀態向量獲取可以直接在用戶子程序中完成，不再需要額外的輸入輸出工作。
上圖仍舊是單自由度彈簧振子模型，只不過在本模型中，剛度和阻尼會隨時間變化，并且將會在0s，1s，5s通過子程序CONSUB分別對系統進行線性化處理，當然在子程序中有c_get_linear_mat_states函數的調用。另外，模型上施加的外部載荷、彈簧力都是通過Gforce實現的，通過函數的綜合應用實現更靈活的設置。
展開 
ABAQUS中使用USDFLD子程序的應用例子 ￥10
今天來解讀一個在ABAQUS中使用USDFLD子程序實現損傷計算的例子。
USDFLD子程序的界面如下：
需要定義的變量
數組FIELD（NFIELD）包含當前材料（積分）點的場變量fi。
——在當前增量步的末尾，這些傳入了由節點插值得到的值，由*INITIAL CONDITIONS選項或*FIELD選項定義。
——更新的fi用于計算材料屬性值，其作為場變量的函數。更新的fi傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。
可能需要定義的變量
• 包含依賴于解的狀態變量的數組STATEV（NSTATV），能夠在USDFLD中定義。
——這些在增量步的開始作為值傳入。
——在該子程序中STATEV在所有情況下都能夠更新，所有更新的值被傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。
——與該材料點相關的狀態變量的數量由*DEPVAR選項定義。
• 建議的新時間步長與現在正在使用的時間步長（DTIME）的比例PNEWDT，可以給定。
——該變量允許用戶在ABAQUS中的自動時間增量步長算法中提供輸入。
展開 ABAQUS——DLOAD子程序應用（移動載荷） ￥9.99
DLOAD子程序可對上述過程進行有效模擬。
DLOAD子程序介紹：
在網上對于單移動載荷的實現辦法介紹較為全面，實現的效果如下：
同時附上我編寫的子程序主體片段，關鍵是需要對移動載荷生效區域和失效區域的準確描述（通過坐標），本例中，車輪與地面間的接觸區域簡化為長方形（寬度即為輪寬）：
但若要實現多個載荷在不同位置同時移動呢？這就需要花費一點心思，觀察子程序的可用參數：對COORDS或者SNAME的加強判斷即可實現多載荷移動，具體效果如下：
收費內容為上述模型的cae文件和上述兩個案例的子程序文件示例
【JY】 ABAQUS子程序UEL的有限元原理與應用
采用FORTRAN77 輸出變量的方式：
（假設變量為雅各比矩陣Jacobi）PRINT *, 'Jacobi=' ,Jacobi
4.編寫完成，開始debug
編寫完UEL子程序之后，先將UEL子程序應用于單個單元進行測試（從靜力、線性攝動、動力時程等各個方面逐一測試），在保證單個單元的計算結果正確無誤的前提下，再將UEL子程序應用于較多單元的分析案例中，可保證分析結果的準確性。
【在ABAQUS中的UEL單元研究】
根據以上平面三角形單元的有限元分析思路，編寫對應的平面二維三角形單元UEL子程序，并通過兩個不同的有限元分析算例驗證該子程序的有效性。
展開 考慮溫度影響的vumat子程序在木材受火后強度分析中的應用
本文編寫了考慮溫度影響的木材本構vumat子程序，并對子程序的正確性進行了驗證。
木材是一種復雜的各向異性復合材料，可以分為三個方向，即縱向順紋(L)、徑向橫紋(R)、弦向橫紋(T)。其復雜的本構關系主要體現為在拉或剪力作用下發生 脆性破壞，而在壓力作用下發生塑性變形，且在橫紋壓力作用下變形較大，同時拉壓強度不相等。
屈服準則
木材是各向異性材料，且L、R、T三個方向的拉壓屈服強度不一樣，屬于拉壓非對稱材料。為了準確地預測木材的失效需要選擇合適的各向異性屈服準則，目前常用的各向異性屈服準則有：Hill準則，Hosford準則，Yamada-Sun屈服準則等。本采用Hashin準則作為木材的屈服準則
損傷演化準則
本文木材本構關系模型定義了兩種不同的損傷演化模型，受壓延性破壞采用理想彈塑性模型，受拉脆性破壞采用線彈性軟化模型.引入損傷變量D來描述木材的受損狀態。
溫度影響
木材隨著溫度的升高，發生不同程度的炭化，其強度、彈性模量、斷裂能也隨之發生變化。本文考慮了溫度對木材的模量、強度的影響，并且認為溫度對拉壓性能產生的影響不同。
根據上述相關理論編寫了abaqus vumat子程序，并通過單胞模型對子程序進行驗證。
下圖為不同溫度下單向拉壓結果
下圖為三點彎曲載荷下的破壞行為
最后，歡迎通過公眾號“320科技工作室”聯系我們
展開 