abaqus只有動力隱式的案例
abaqus顯示動力學VS隱式動力學 ￥29.99
引言：本文內容綜合參考了《ABAQUS 6.12 有限元分析從入門到精通》、ABAQUS 官方幫助文檔以及《ABAQUS 有限元分析常見問題解答》等資料，同時結合個人在學習與實際應用過程中的體會與思考，旨在幫助讀者對顯示分析步與隱式分析步的差異有更加深入的理解。
需要特別說明的是，文中觀點部分基于作者的學習與實踐經驗，難免存在不足或偏差，誠摯歡迎同行提出寶貴意見與建議，以便相互交流、共同進步。
1、 通用隱式分析步：
圖1為創建“動力，隱式”后的“基本信息”“增量”“其他”三個選項卡。
圖1 隱式動力學分析步
在設置分析步時，“增量”和“其他”兩個選項卡往往容易被忽視。一般來說，選擇自動時間增量時可以通過Half-step residual控制平衡殘差的容差，以兼顧精度與效率；而固定時間增量則可啟用Suppress half-step residual來跳過殘差檢查，加快計算，但可能犧牲穩定性。在“其他”選項卡中，求解技術不涉及接觸迭代，載荷默認按瞬態方式隨時間變化；至于初始加速度，如果是第一個動力學分析步則為零，如果前一步同樣是動力學步則沿用其結束時的加速度，默認情況下ABAQUS會自動計算，但若確認載荷無突變則可關閉以節省運算量。
2、 通用顯示分析步
該分析步用于顯式動力學分析，除了“基本信息”“增量”和“其他”三個選項卡頁面外，其“編輯分析步”對話框還包括一個“質量縮放”選項卡頁面。“基本信息”選項卡頁面中的幾何非線性選項默認為“開”。“增量”選項卡頁面的相關參數如表1所示。
展開 基于ABAQUS顯式動力學和隱式動力學的彎管成型加工分析 ￥50
摘要： 彎管成型分析方法及注意點，分別采用顯式動力學和隱式動力學方法進行彎管成型分析，并對兩種分析方法進行了比較。
素材來源：
本文中所引用的案例素材來源于：錢 峰，潘笑譽，葉小奔. 基于ＡＢＡＱＵＳ的管件彎曲成型的有限元分析[J]. 機械工程與自動化,１６７２－６４１３（２０１７）０４－００１９－０３.
案例中彎管模具以及坯料尺寸皆根據該論文進行建模，論文中未提及的尺寸信息，筆者根據自身過往經驗自行設定。
建 模：
零部件定義：共五個組件，除坯料外，全部設置為離散剛體，解析剛體也可，這里以離散剛體為例。
材料屬性：剛體不需要設置材料，坯料材料屬性需要設置彈性模量，應力應變數據以及密度（采用隱式動力學可不設置密度參數）。
關于應力應變曲線定義需要注意abaqus程序內部對應力應變曲線數據的正則化問題，顯式動力學算法會根據用戶輸入的應力應變數據進行正則化處理，如果數據點過多，或者數據點間存在突變，則可能會產生計算誤差或者直接報錯。
莊茁教授關于顯示動力學中應力應變關系的正則化處理的解釋
分析步：
隱式動力學設置如下：
顯式動力學設置如下：
關于顯式動力學中質量縮放（Mass scaling）的原理請查閱官方幫助文檔。由于一般的材料加工問題都屬于準靜態問題，因此采取一定的質量縮放可以在損失較少計算精度的前提下，大幅提高計算速度，具有非常大的應用價值。當然，如果能夠通過采用更大的網格劃分提高計算速度，那是更好的選擇。在模型調試階段可以采用較大的質量縮放系數，快速計算，debug模型，在模型調整好之后減小質量縮放系數，以獲得精度更高的計算結果。
展開 如何用Abaqus隱式動力學畫出李薩如圖形？ 附Abaqus動力學有限元分析指南張文元下載
李薩如沙擺的模擬可以用Abaqus顯式動力學，考慮到我們的顯式分析應用案例比較多，這期文章我們換一種方法，使用Abaqus隱式動力學來計算這個過程。
Abaqus隱式動力學使用隱式時間積分（Hilber-Hughes-Taylor算法、向后的歐拉算法）來計算系統的瞬態動力學或準靜態響應，首先簡單的看一下Abaqus中隱式動力學的幾種應用方案：
瞬態保真（Transient fidelity）
不含接觸模型的默認選項，涉及最小的系統能量耗散，比如衛星系統的分析，使用較小的時間增量來精確求解結構的振動響應，本文的沙擺采用這種方案。
隱式瞬態保真應用-沙擺振動
中度耗散（Moderate dissipation）
包含接觸模型的默認選項，涉及中度的系統能量耗散，比如動力學系統通過塑性、黏性阻尼或其他效應進行能量耗散，可以用于各種插拔、碰撞和成型分析。
隱式中度耗散應用-棘輪碰撞
準靜態（Quasi-static）
準靜態分析的選項，主要感興趣的是最終的靜態響應，涉及高度的能量耗散，通過引入慣性效應來規范不穩定行為，比如因欠約束導致的剛體位移或“突然跳變”。一個應用場景是指甲刀的捏合分析，首先通過添加和釋放輔助約束的靜力學方法來計算指甲刀的裝配應力，然后將所有部件的相互作改用接觸定義，模型中增加了很多不穩定因素，如果繼續使用靜力學則極其容易發散，改用隱式準靜態可以順利完成計算。
展開 請教abaqus隱式動力求解問題
abaqus做隱式動力學的時候提交文件后一直是running狀態，等了幾個小時還是這樣，結果文件也不更新，sta文件也沒有，一直卻沒有報錯，請問是什么原因？
ABAQUS 靜/動力隱式UEL 二次開發 ￥23456
<p><strong>概述</strong>：結合HHT時程積分法，推導了ABAQUS <strong>靜/動力隱式算法</strong>中的UEL關鍵矩陣表達式，并將公式應用到自編CPS4/CPE4、C3D8 BBAR和C3D20用戶自定義單元中，計算結果均與ABAQUS自帶單元保持一致。其中，靜力計算中，關鍵矩陣AMARTX和RHS等可直接按照剛度矩陣和方程右端不平衡力輸出。動力隱式計算中，這兩者的輸出較為復雜，需要結合HHT時程積分法進行推導，將剛度矩陣、質量矩陣和阻尼矩陣依據LFLAGS數組的數值進行組合，RHS同樣需要進行推導計算，并以合適的方法將解相關的狀態量儲存在SVARS中，供后面的增量步調用。</p><p>關于UEL的程序設計，只支持靜力通用計算分析步的資料有很多，本帖子內容適不但適用于靜力通用，而且適用于動力隱式、頻率分析等分析步計算，采用模塊化程序設計，所有的矩陣求解均被封裝，調用方便，讀者可以針對自己的需求對相應函數進行改編即可，尤其是涉及到動力隱式計算部分，<strong>適用于任何運動方程的動力隱式求解，可直接移植使用</strong>。
展開 ABAQUS 二次開發（UEL+UMAT C3D8 單元動力隱式計算）
概述：采用UEL接口二次開發實現八節點單元，考慮BBAR修正，避免體積自鎖，對標ABAQUS自帶的C3D8單元，計算的剛度矩陣、質量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進行應力和應變數據的可視化，計算的應力應變數據同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動力簡諧荷載為例，采用上述程序，應用動力隱式計算分析步，最終計算的位移、應變等時程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節點均固定。另一面所有的節點施加簡諧荷載，簡諧荷載曲線為（詳細的參數見附件）：
計算的網格圖為：
網格尺寸為1，共計10x10x=1000個單元。總的計算時長為1s，這只固定增量步長為0.01s，所以總增量步數為100。
（）計算結果
以上面網格圖中中間角點為例，提取加載向位移時程曲線如下圖：
0.89s時刻x向應力云圖為：
注：左一為ABAQUS計算結果，中間為umat實現線彈性各向同性本構計算結果，右一為UEL+UMAT計算結果，后面的云圖也是這樣排布。
0.89s時刻y向應力云圖為：
寫文字好麻煩，不寫了！放個動圖算了，有興趣的移步附件。
展開 如何用Abaqus隱式動力學畫出李薩如圖形？
李薩如沙擺的模擬可以用Abaqus顯式動力學，考慮到我們的顯式分析應用案例比較多，這期文章我們換一種方法，使用Abaqus隱式動力學來計算這個過程。
Abaqus隱式動力學使用隱式時間積分（Hilber-Hughes-Taylor算法、向后的歐拉算法）來計算系統的瞬態動力學或準靜態響應，首先簡單的看一下Abaqus中隱式動力學的幾種應用方案：
瞬態保真（Transient fidelity）
不含接觸模型的默認選項，涉及最小的系統能量耗散，比如衛星系統的分析，使用較小的時間增量來精確求解結構的振動響應，本文的沙擺采用這種方案。
隱式瞬態保真應用-沙擺振動
中度耗散（Moderate dissipation）
包含接觸模型的默認選項，涉及中度的系統能量耗散，比如動力學系統通過塑性、黏性阻尼或其他效應進行能量耗散，可以用于各種插拔、碰撞和成型分析。
隱式中度耗散應用-棘輪碰撞
準靜態（Quasi-static）
準靜態分析的選項，主要感興趣的是最終的靜態響應，涉及高度的能量耗散，通過引入慣性效應來規范不穩定行為，比如因欠約束導致的剛體位移或“突然跳變”。一個應用場景是指甲刀的捏合分析，首先通過添加和釋放輔助約束的靜力學方法來計算指甲刀的裝配應力，然后將所有部件的相互作改用接觸定義，模型中增加了很多不穩定因素，如果繼續使用靜力學則極其容易發散，改用隱式準靜態可以順利完成計算。
隱式準靜態應用-指甲刀捏合
沙擺模型中不含接觸，能量耗散比較小，因此宜采用Abaqus的隱式瞬態保真（Transient fidelity）進行計算。
展開 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階頻率信息。
展開 