Model change的案例
*Model change 應用技巧進階
在assembly 中設置“需要鈍化”的單元集，以便在model change進行相關操作。如果不在此處設置單元集，則在model change調用的時候需要指定part屬性。
3.添加“完全彈性”材料
要求完全彈性，即相對于非鈍化單元的彈性模量小，最好能小5個數量級；要求質輕，以避免自重引起變形過大，可近似取自重密度比其彈性模型模量小5個數量級。
4.在適當階段鈍化單元集
** Remove all elements in element set Original_Set
*MODEL CHANGE, REMOVE
Original_Set
需要注意的是鈍化掉的是原始單元集，而不是“追蹤單元集”，否則無法實現追蹤功能。
5.在需要激活的階段激活已鈍化單元集
** Reactivate elements in element set Original_Set
*MODEL CHANGE, ADD=STRAIN FREE
Original_Set
在此步中也可以將“追蹤單元集”remove掉。
附件為修改后的模型，共學習參考。該模型在model change命令中采用的是直接調用instance內部的單元集。也可以在assembly中定義單元集后直接調用單元集。
testModelchage-3.rar
展開 ABAQUS生死單元（Model change）入門
幫助文檔這樣描述“Amodel change interaction allows you to deactivate and reactivate elements to simulate removal of part of the model, either temporarily or for the remainder of the analysis. You can create a model change interaction in all
Abaqus/Standard steps except for static, Riks steps and linear perturbationsteps. For a detailed discussion, see Element and contact pair removal andreactivation.”
通過這段說明我們可以知道，生死單元的使用條件是：in all Abaqus/Standard steps except for static, Riks steps andlinear perturbation steps。
生死單元在CAE界面實現方法為： create interaction》選擇需要移除或者激活的分析步 step并選擇model change，如下圖中選擇step-3使用 model change。
點擊繼續之后，可選取相應的區域來施加model change以及相關選項。
下面這個案例簡要的演示了model change的使用，該案例為一個簡單的拉伸模型，僅為演示用法，并無特殊的工程意義，讀者在使用生死單元的時候可根據自己的實際工況加以靈活選用。
展開 Model change 應用技巧進階（隧道開挖應用）（轉）
在assembly 中設置“需要鈍化”的單元集，以便在model change進行相關操作。如果不在此處設置單元集，則在model change調用的時候需要指定part屬性。
3.添加“完全彈性”材料
要求完全彈性，即相對于非鈍化單元的彈性模量小，最好能小5個數量級；要求質輕，以避免自重引起變形過大，可近似取自重密度比其彈性模型模量小5個數量級。
4.在適當階段鈍化單元集
** Remove all elements in element set Original_Set
*MODEL CHANGE, REMOVE
Original_Set
需要注意的是鈍化掉的是原始單元集，而不是“追蹤單元集”，否則無法實現追蹤功能。
5.在需要激活的階段激活已鈍化單元集
** Reactivate elements in element set Original_Set
*MODEL CHANGE, ADD=STRAIN FREE
Original_Set
在此步中也可以將“追蹤單元集”remove掉。
testModelchage-3.rar
展開 *Model change 的應用技巧初步
為此，要求完全彈性，相對于非鈍化單元彈性模量小，最好能小5個數量級（有待進一步研究），相對質量要輕（因為避免自身變形過大，就得考慮自重密度與彈模之間的比例關系，過自重密度與彈模比較大，則變形必定很大，暫定5個數量級）
3>assembly
a）按照正常程序進行assembly
b）在此階段需設置多個單元集和節點集（至少有“鈍化單元節點集”，“非鈍化單元節點集”，“備份鈍化單元節點集”）
c）對備份鈍化單元集的要求
與鈍化單元節點號完全一致；與鈍化單元數量完全一致；單元總數會因此發生變化，增量數即為備份鈍化單元數；
注意，經測試，鈍化單元集只能在assembly階段設置，否則無法鈍化和激活
4>step
按一般方法設置
a）在第一階段鈍化單元集
** Remove all elements in element set EL-SET
*MODEL CHANGE, REMOVE
EL_SET
b）在需要激活的階段激活鈍化單元集
** Reactivate elements in element set EL_Set
*MODEL CHANGE, ADD=STRAIN FREE
EL_SET
c）“完全彈性單元集”不能鈍化，否則失去追逐坐標的功能
5>load
注意，各工況的計算結果是否需要累計
6>mesh
一般方法
7>job
在job管理器中，create —source file(inp) 提交inp文件
備注，目前是適于于ADD=STRAIN FREE的情況，其他情況尚未測試。“完全彈性材料”彈性模量的取值以及材料非線性以及與時間有關的因素也有待進一步探討。
附錄為一個T形梁的測試文件，即在第一個階段鈍化翼緣板，在第三個階段激活該翼緣板。
Modelchage_test.rar
展開 
Model change 的應用技巧初步
為此，要求完全彈性，相對于非鈍化單元彈性模量小，最好能小5個數量級（有待進一步研究），相對質量要輕（因為避免自身變形過大，就得考慮自重密度與彈模之間的比例關系，過自重密度與彈模比較大，則變形必定很大，暫定5個數量級）
3>assembly
a）按照正常程序進行assembly
b）在此階段需設置多個單元集和節點集（至少有“鈍化單元節點集”，“非鈍化單元節點集”，“備份鈍化單元節點集”）
c）對備份鈍化單元集的要求
與鈍化單元節點號完全一致；與鈍化單元數量完全一致；單元總數會因此發生變化，增量數即為備份鈍化單元數；
注意，經測試，鈍化單元集只能在assembly階段設置，否則無法鈍化和激活
4>step
按一般方法設置
a）在第一階段鈍化單元集
** Remove all elements in element set EL-SET
*MODEL CHANGE, REMOVE
EL_SET
b）在需要激活的階段激活鈍化單元集
** Reactivate elements in element set EL_Set
*MODEL CHANGE, ADD=STRAIN FREE
EL_SET
c）“完全彈性單元集”不能鈍化，否則失去追逐坐標的功能
5>load
注意，各工況的計算結果是否需要累計
6>mesh
一般方法
7>job
在job管理器中，create —source file(inp) 提交inp文件
備注，目前是適于于ADD=STRAIN FREE的情況，其他情況尚未測試。“完全彈性材料”彈性模量的取值以及材料非線性以及與時間有關的因素也有待進一步探討。
附錄為一個T形梁的測試文件，即在第一個階段鈍化翼緣板，在第三個階段激活該翼緣板
Modelchage_test.rar
T梁的測試文件inp格式
展開 abaqus中的漸進單元激活與model change
漸進單元激活相對于model change，可以對單元的激活進行更精確的控制。
1.Model change（單元以及接觸對的移除與重新激活）
· 可用于模擬從模型中刪除part，可以是暫時的或在剩下的分析中永久刪除;
· 允許單元在無應變或有應變的情況下重新激活;
· 當不需要接觸對時，可用于節省計算時間;
· 只能用于通用分析步;
· 只有在原始分析中使用或激活時，才能在重新啟動分析中使用。
您可以在通用分析步中從模型中刪除指定的單元。就在刪除的分析步之前，Abaqus/Standard存儲了被刪除區域在它們之間的邊界節點處對模型其余部分施加的力/通量。在刪除對應的分析步中，這些力降為零；因此，只有在刪除分析步結束時，刪除區域對模型其余部分的影響才完全消失。力逐漸下降，以確保刪除對模型的影響平滑。
2.Progressive element activation（漸進單元激活）
· 可用于在分析過程中激活單元;
· 可用于在單元激活過程中應用本征應變;
· 只能與實體和殼單元一起使用;
· 只能用于熱傳導或靜力分析。
可以在分析步的每個增量步中激活單元。必須首先定義在分析過程中可以激活的單元，然后在每個可以激活它們的分析步中引用它們。激活特性在一個分析步中被打開的單元可以通過在每個增量的開始為一個單元分配一個材料的體積分數來激活。
支持單元的完全和部分激活。對于完全激活，添加的材料體積分數必須等于0或1(也就是說，單元的狀態只能從不激活變為完全激活)。對于部分激活，加入的材料體積分數可以是任意的;然而，在實際應用中，一個單元的體積分數不能太小，從而不能避免數值奇異性問題。
在應力-位移分析中，假定加在單元上的材料是無應力的。
展開 單元生死在ABA/CAE的實現
Model Change以前只能在隱式分析中通過inp操作（顯式分析不支持Model Change功能），現在v6.10版可以ABA/CAE的實現，近來很多人問起如何操作。工作過程中尚未用到這類工況，今天隨便假想一個工況，做一個小實例。
Model Change與其他類型分析的唯一區別是在interaction management---create---選擇相應的step----model change--按照圖中步驟選取你需要刪除或者添加的幾何（或者集合）。
模型簡介：
梁兩端固支，在中間劃分出一小段用于remove。建立2個step：
step1 在前2段加載均布壓力p
step2 將中間一段remove，并將均布壓力增大。
結果顯示：
step1的結果和一個兩端固支梁加載無異；
step2的結果顯示右端沒有應力。
附：
1）inp文件和6.10版的cae文件：
6.10版的cae文件
model-change.rar
inp文件
change.rar
展開 生死單元（Model change）中網格副本的建立
生死單元中網格副本的設置.pdf
一邊督促自己，一邊幫助別人，我開心，你快樂，世界多美好！
ABAQUS焊接分析- Part 1:手動定義
激活/抑制單位
在模型最初的情況下焊料需要被去除，然后在分析過程中使用Model change命令添加焊料。
焊接過程中焊料熔化，在仿真的過程中可以理解為焊料從無到有。此過程時間非常短，僅為1E-7s。這樣，有無焊料都假定溫度場這個短時間內不會發生明顯改變。
在Create Interaction中找到Model Change（圖2）。
圖2：Model change窗口
在分析步1去除焊料（Deactivated）。（圖3）
圖3：在分析步1去除焊料
在分析步4，修改Interaction，重新激活單元，以添加焊料。
添加初始溫度
使用預定義場定義幾何初始溫度。工件溫度為室溫。焊料溫度為高溫（1500℃）。
添加邊界條件
在分析步3中，在工件和焊料之間的邊緣處施加一個漸增的溫度邊界條件，在分析步3中進行，在分析步4中移除。
圖4：添加溫度邊界條件
添加薄膜條件
在結構的外表面不同區域手動選擇施加薄膜條件。
展開 新手模擬隧道開挖的常用幾種方法
采用model change 控制圍巖開挖和襯砌加入。但由于TIE約束本身的缺陷，該方法劃分網格后，節點調整后，對于大的3D模型，求解和精度都有一些問題。
2、隨著認知的增加，漸漸你會發現，其實一般隧道的模擬不是這么做的，一般是在只建立一個part，partition出襯砌和開挖塊。對于單獨襯砌建立集合和小集合后。劃分網格，在job模塊，write input，然后在inp文件中用elcopy對襯砌集合生成另一個節點號相同但是單元號不同的集合。對該集合賦予襯砌性質。最后再采用生死單元model change。
3、采用field和temperature來變化材料性質。這個方法很明了，就是只能采用的本構模型比較簡單。實際模擬不好用。inp文件下次再傳。
4、采用input方法，直接導入模型節點，和各材料的單元。可以借鑒《ABAQUS在隧道與地下工程中的應用》一書。其實該書還是不錯的。雖然很多童鞋覺得不好、但是abaqus雖然作為CAE軟件，input文件的修改和書寫還是很關鍵的能力。還是要像ANSYS一樣啃命令流呀。查閱下自帶的幫助文件。
展開 基于Abaqus的生死單元技術的焊接仿真分析 ￥11.99
編輯模型屬性（Model-->Edit Attributes）
設置絕對零度(-273.15)和波爾茨曼常數(5.67E-8)。
8. 生死單元實現
（1）查看焊錫部分兩端的單元編號，以及單元編號的分布規律。本例中焊錫部分每層有8個單元，兩端的單元編號分別為1~8和793~800。
（2）采用Python腳本（參考附錄：第二部分（循環編輯關鍵字Model change））實現關鍵字Model change的自動添加。
9. Job模塊
建立Job，在General選項卡添加User subroutine file（Fortran文件），并提交。
展開 
隧道開挖模擬-2
.,-1.
（2）修改inp語句，在geostatic分析步中移除all2，在襯砌生成步中重新激活
*Step, name=Step-1, nlgeom=NO
*Geostatic
**
*Model Change, remove
part-2-1.all2
*Step, name=Step-3, nlgeom=NO
*Static
1., 1., 1e-05, 1.
*Model Change, add
part-2-1.all2
圖3和圖4是開挖后和激活襯砌及接觸面后的豎向位移圖，激活襯砌后土體豎向變形沒有變化，符合實際。
圖3
圖4
關注點：
（1） 隧道開挖后的土體的上抬變形，有可能是abaqus調整接觸面上的過盈量造成的。
（2） 建立接觸面時，注意土體內部面的方向，否則襯砌和土體之間不會有接觸約束。
來源：ABAQUS在巖土工程中的應用
作者：費康
展開 基于Abaqus的光固化3D打印結構仿真
另外一個需要考慮的是樹脂固化過程中額體積收縮，體積收縮與固化度的關系采用三線性模型
用到的子程序：umat(材料本構),uexpan(熱應變和化學收縮應變)
建模方法
3D打印時，材料是層層疊加的，因此我們這里需要用到Abaqus中的Model Change功能，對材料進行分層激活（新版本中多了UEPACTIVATIONSETUP子程序，也可以實現與Model Change相同的功能）。由于3D打印層數非常多，直接人工進行建模需要耗費大量時間，我們通過Python腳本對建模過程進行簡化。腳本的建模思路為：按層批量分割實體，批量建立熱力耦合分析步以及按分析步激活實體。
模擬結果
正向打印
反向打印
固化度和溫度變化趨勢
展開 重磅發布！Abaqus焊接仿真指南V2.0：從DFLUX子程序到FSW全流程詳解 ￥29.9
3?? 生死單元技術（Model Change）
想模擬真實的材料填充過程？必須掌握生死單元。文檔詳細演示了如何在Abaqus中設置 Model Change，以及如何通過 Python 腳本 自動創建大量的Set集和分析步，告別機械重復的體力活。
4?? 進階：攪拌摩擦焊（FSW）CEL法仿真
針對復雜的固相連接工藝——攪拌摩擦焊，指南中詳細講解了基于 CEL（耦合歐拉-拉格朗日） 方法的建模全流程。 從歐拉域的網格劃分、體積分數填充，到攪拌頭的剛體設置、下壓/旋轉/移動的邊界條件加載，再到使用Meta進行后處理，全流程無死角覆蓋。
?? 為什么你需要這份南？
實戰導向：不是枯燥的幫助文檔翻譯，而是基于真實案例的操作手冊。
圖文并茂：關鍵步驟均有軟件截圖，參數設置一目了然。
源碼交付：配套的 INP文件、Fortran子程序文件、Python腳本，我都準備好了。
?? 如何獲取？
歡迎加入我的 知識星球。在星球內，你不僅可以下載到：
《焊接模擬工程師指南 V2.0》高清PDF版 (見下方付費內容)
TIG焊接完整INP文件 + DFLUX子程序源碼 (.f)
Python自動化腳本 (.py)
FSW攪拌摩擦焊完整模型文件
更重要的是，你在練習過程中遇到的報錯、環境配置（VS+IVF+Abaqus）等問題，都可以在星球內向我提問，我會一對一解答。
?? 掃碼加入星球，開啟你的焊接仿真進階之路！
寫在最后：仿真之路，道阻且長。這份V2.0版本的指南，是我對近期焊接仿真經驗的一個總結。希望能幫大家少走彎路，快速上手。如果你覺得有用，請點個“在看”支持一下！
展開 ABAQUS在產品開發中的應用及示例
STEP1中，接觸被去掉（*MODEL CHANGE， REMOVE， TYPE = CONTACT PAIR）；STEP2中，接觸被加上（*MODEL CHANGE， ADD， TYPE = CONTACT PAIR）；STEP3中，接觸仍然存在。
5．邊界條件的處理：
在STEP1中，可以只定義彈簧夾片下部的接點（DOF=1，6），將自由端取一排結點，通過結點位移方法使彈簧夾片打開（DOF3=-5）；在STEP2中，將彈簧夾片卸載，即取消自由端結點的邊界條件（通過定義其它邊界條件時，令OP=NEW），讓其自由回彈，夾緊接觸部位。此時，附加合理約束來保證鏡盒本體平衡，在STEP3中再將它去掉；在STEP3中，在轉動點處向下加載位移（DOF3=65）來保證鏡盒本體向下饒轉軸轉動90度。
6．結果輸出要求： 輸出：應力、應變、位移，接觸力的時間曲線（用于判斷接觸的情況）
7．結果判斷準則：考慮到疲勞問題，要求一次操作時頸部應力不應超過15MPa。
8．前處理、求解及后處理： 用ABAQUS/CAE完成前后處理，用ABAQUS/Sandards來求解。
9．結果如下：
A：接觸反力曲線：
由圖可見，在T=1.37左右，彈簧夾片與本體開始接觸，在T=2.0時，接觸力達到最大。在STEP3中(T：2.0~3.0)，接觸力逐漸下降，在T=2.5左右出現接觸波動(不連續)，在接近T=3.0左右時，接觸力降為0，說明彈簧夾片與鏡盒本體脫開。
B：鏡盒本體應力云圖：
如圖所示，雖然有一個單元的應力較高（VON MISS STRESS=14.5MPa），但是不能作為評價值，而頸部面積較大的綠色區域為真實的應力情況，VON MISS STRESS=8.5MPa。該值滿足設計要求。
展開 