隱式-顯式順序求解的案例
LSDYNA-隱式－顯式順序求解
LSDYNA-隱式－顯式順序求解，電腦中存的資料
想學習更多的知識，請聯系我們！
微信公眾號：名稱：“DR有限元”
號碼：“hello_cae”
顯式與隱式計算
這是關于LS-DYNA的顯式與隱式計算地 ppt 講義。
顯式隱式.rar
隱式－顯式順序求解.rar
Ansys vs Abaqus：隱式與顯式求解的終極博弈
在CAE領域，選擇Standard（隱式）還是Explicit（顯式）求解器，本質上是在平衡“計算精度”與“時間尺度”。
1?? 隱式求解 (Implicit/Standard)
核心是求解 $Ku=F$。每一步都需要進行矩陣求逆和牛頓迭代，以確保力平衡。
特點： 絕對收斂。步長可以很大，不受穩定性限制。
擅長： 靜力學、線性振動、緩慢的非線性過程。
痛點： 接觸極度復雜或大變形時，收斂困難，報錯“收斂失敗”是常態。
2?? 顯式求解 (Explicit)
核心是動力學方程 $Ma=F-I$。直接根據當前時刻的狀態推導下一時刻，不求逆陣，不迭代。
特點： 沒有收斂問題。但步長受限于穩定性準則（CFL條件），通常極小（$10^{-7}$s量級）。
擅長： 跌落、碰撞、爆炸、高速切削。
痛點： 適合極短時間內的物理過程。計算長時間問題時，累計誤差大。
3?? 工具選型建議
Abaqus： Standard與Explicit切換極其絲滑，適合處理復雜的非線性接觸（如密封件、橡膠）。
Ansys： 隱式求解器極其高效穩定，配合LS-DYNA插件，在結構靜力和多物理場耦合上具有統治力。
展開 顯式動力學出現必須要用隱式求解的警告
模型設置：使用hinge連接器，一端連接剛體參考點，一端使用運動耦合連接軸孔，在使用顯示求解的時候運動耦合參考點出現了下面的錯誤，請問該錯誤會影響結果嗎，如果會該怎么處理。謝謝
『分享』幾個對于初學者來說還不錯的ansys資料
基礎篇.doc
顯式隱式.ppt
隱式－顯式順序求解.ppt
ls-dyna適合初學者學習
3 Part 定義.pdf
2 單元.pdf
4 材料的定義.pdf
1 概述.pdf
5 加載，剛性體和邊界條件.pdf
7 求解和模擬控制.pdf
9 重啟動.pdf
8 后處理.pdf
10 顯式－隱式順序求解.pdf
6 接觸面.pdf
11 隱式－顯式順序求解.pdf
DYNA從入門到精通
<p>ANSYS LS-DYNA教程</p><p>1 概述</p><p>2 單元</p><p>3 Part 定義</p><p>4 材料的定義</p><p>5 加載，剛性體和邊界條件</p><p>6 接觸面</p><p>7 求解和模擬控制</p><p>8 后處理</p><p>9 重啟動</p><p>10 顯式－隱式順序求解</p><p>11 隱式－顯式順序求解</p><p>ANSYS LS-DYNA實例</p><p>1 風機葉片包容性分析</p><p>2 顯式動力學單元</p><p>3 接觸分析中Part的運用</p><p>4 沖擊分析中的塑性和失效</p><p>5 點焊失效分析</p><p>6 鋁棒沖擊分析</p><p>7 受軸向載荷梁的屈曲分析</p><p>8 用剛性沖頭彎管</p><p>9 小型重啟動分析</p><p>10 回彈分析</p><p>11 施加熱結構預載荷</p><p><br></p><p>整理不易，請點贊自取</p><div contenteditable="false" width="100%"><p><img src="https://img.jishulink.com/static/web/attachment.png" style="display:inline;vertical-align: middle;width: 24px;height:24px;">
<a href="https://img.jishulink.com/202512/attachment/9823d2d3708a4f40b969d0b5f92853c2.7z" target="_blank">ANSYS LS-DYNA教材.7z</a></p></div><p><br></p><p><br></p><p><br></p><figure style="text-align:
展開 ABAQUS顯式與隱式的區別
ABAQUS中動態分析包括兩大類基本方法：
振型疊加法：用于求解線性動態問題；
直接積分法：主要用于求解非線性動態問題。
ABAQUS顯式（explicit）和隱式（standard）算法分別對應著直接積分法中的中心差分法（顯式）和Newmark（隱式）法等。
比較兩種算法，顯式中心差分法非常適合研究波的傳播問題，如碰撞、高速沖擊、爆炸等。顯式中心差分法的M與C矩陣是對角陣，如給定某些有限元節點以初始擾動，在經過一個時間步長后，和它相關的節點進入運動，即U中這些節點對應的分量成為非零量，此特點正好和波的傳播特點相一致。另一方面，研究波傳播的過程需要微小的時間步長，這也正是中心差分法的特點。
而Newmark法更加適合于計算低頻占主導的動力問題，從計算精度考慮，允許采用較大的時間步長以節省計算時間，同時較大的時間步長還可以過濾掉高階不精確特征值對系統響應的影響。隱式方法要轉置剛度矩陣，增量迭代，通過一系列線性逼近（Newton-Raphson）來求解。正因為隱式算法要對剛度矩陣求逆，所以計算時要求整體剛度矩陣不能奇異，對于一些接觸高度非線性問題，有時無法保證收斂。
下面分別介紹這兩種算法
abaqus 顯式與隱式的區別.pdf
展開 本周討論熱點：顯式算法與隱式算法
歡迎大家討論啊，希望各位參與進來，一周后結貼給分
Abaqus隱式轉顯式分析 ￥10
我們在做顯式動態分析時經常會碰到模型中需要考慮螺栓預緊力，重力場，過盈配合，預應力等的情況，此時我們便需要用到abaqus隱式轉顯式的方法。
通過一個預緊力的小例子（隱式加載9000N預緊力，顯式工況為空載）來加以說明。
計算結果
通過調整接觸算法，得到誤差更小的接觸力。(方法2預緊力誤差為0.22)
如何在abaqus實現顯式分析與隱式分析交叉進行
如題，如何在abaqus中如何在abaqus實現顯式分析與隱式分析交叉進行
ABAQUS中隱式和顯式的節點和單元的輸出變量解析
The output variables listed below are available in Abaqus/Explicit.
Mechanical analysis–nodal quantities
CFORCE
Field: yes History: no .fil: no
Contact normal force (CNORMF) and frictional shear force (CSHEARF).
CDISP
Field: yes History: no .fil: no
Contact opening (COPEN) and accumulated tangential motions (CSLIP1, CSLIP2, and
CSLIPEQ) for general contact analyses.
CEDGEACTIVE
Field: yes History: no .fil: no
Status of contact edges for general contact analyses (active as primary, active as
secondary and deactive).
CFRICWORK
Field: yes History: no .fil: no
Contact frictional work for general contact analyses.
CNAREA
Field: yes History: no .fil: no
Contact nodal area for each node with active contact forces in general
展開 【快鏈效應】與自由落體彈簧-隱式結果導入顯式
在之前介紹魔力彈簧Slinky的文章中，我們講到了它有趣的下樓梯動作和詭異的自由落體“懸浮”現象。
魔力彈簧下樓梯
魔力彈簧自由落體“懸浮”現象
魔力彈簧自由落體“懸浮”現象的模擬用到了Abaqus中的一個分析連續性技術，即從Abaqus/Standard分析轉換到Abaqus/Explicit分析，時間比較久了還有人問到這個模型，今天的USim小課堂我把建模關鍵細節說明一下。
Slinky自由落體“懸浮”試驗
→ Standard模型
1. 導入外部的彈簧幾何模型文件，或者通過刪除特征的建模思路，直接在Abaqus中創建一個3D的線體彈簧，具體方法如下：
首先通過Revolution創建一個帶升角的回轉片體結構，然后通過Shell edges生成Wire，并刪除多余的邊，得到一個三維線體的彈簧。
創建彈簧1-旋轉生成3D的片體彈簧
創建彈簧2-生成Wire并刪除多余的邊
2. 設置一個矩形梁的截面屬性（1×4）并賦予彈簧，我們只關心彈簧在自身重力下的最終響應，所以這一個分析采用靜力學分析步，考慮幾何非線性，設置重啟動輸出的間隔為1，Load模塊定義重力、固定彈簧上端點。
彈簧矩形截面梁
彈簧在自重下的變形
靜力分析結果表明，由于彈簧受到重力并且上端固定，其變形和應力自上而下逐漸減小，這個分析最后一個增量步的計算結果要傳遞到下一個分析。
彈簧在自重下的應力
→ Explicit模型
1. 在模型樹中復制Standard模型并重命名（名稱隨意），用作第二個分析，替換Static，General分析步為Dynamic，Explicit，
展開 huang隱式程序修改為顯式及計算案例
然而一些特殊的工況，如切削，軋制，沖壓等隱式存在收斂性問題。因此通常使用顯示程序進行計算。但從頭完成顯式晶體塑性構造對于一般學者顯然難度過高，一個簡單的想法就是直接將現成的黃永剛隱式程序改成顯式。abaqus里這是可以實現的。其基本的步驟是：
1，加入vumat接口程序（見附錄abaqus官網有）
2，對nblock進行循環，計算應力和狀態變量
3，更新應力與狀態變量，重復計算直到增量結束。
值得注意的是，umat與vumat程序里面剪應力分量定義順序與應力不同
umat：12，13，23（工程剪應變）
vumat：12，23，13（2*工程剪應變）
同時采用該方法計算時計算效率顯著高于完全顯式，并允許較大的時間增量。為評估模型計算效率，采用1000個晶粒80000個單元的二維模型進行20%的壓縮模擬。耗時3小時，計算結果與隱式結果類似。
展開 6-有限元之顯式算法和隱式算法
隱式算法（implicit method）（backward Euler method）
考慮同一個方程，在T（n+1）時刻有：
所以在隱示算法中，T(n+1)時刻的值不光由T(n)時刻決定，還由當前時刻T(n+1)決定。也就是說當前時刻的值由上一時刻和當前時刻的值共同決定。隱式算法往往需要求解二次方程。
我們來看看一個具體事例：
設常微分方程：
根據上面的方法，對于顯示算法有：
得出：
對于隱式算法有：
導出二次方程：
求解得:
所以很明顯，在隱式算法中，要求得K+1時刻的值，就需要求解二次方程的根。
關于收斂性
顯式算法不存在收斂性的問題（因為不進行收斂計算），從方程中可以看出來，每個時刻的值由上一時刻所確定，所以一步一步進行下去，當時間步取得較大時，就會偏離真實值。
顯式算法的過程（藍色為真實值）
隱式算法是無條件收斂的，在隱式算法中，在求解二次方程的同時，會通過Newton–Raphson method算法對每一步進行迭代收斂，直至收斂到指定的偏差。如下圖所示：
時間步長（time integration）的依賴性（時間變量只在動力學中涉及）
顯式算法要獲得準確的結果，需要取很小的時間步長
隱式算法對時間步長要求不高，由于是絕對收斂的，往往可以取較大的時間步長。
運用上面的方法，我們以方程為例，通過數值算法求得f(u)。
展開 