動態松弛的案例
關于lsdyna動態松弛卡片的一些使用建議
在一些分析中我們經常會考慮到預應力或者其他一些preload的影響，或者使用動態松弛來消除動能進行多次沖擊，這時候選擇使用動態松弛是一個不錯的選擇。
關鍵字：*CONTROL_DYNAMIC_RELAXTION。
1.考慮預應力的影響
一般來說，卡片使用默認設置就行，其中有兩個設置選項需要注意：
1：DRTOL：該卡片為收斂容差的選擇，默認為0.001，你可以設置自己想要的收斂容差，想要快速收斂，可以設置成0.1。
對于螺栓預緊力來說，默認設置都能很快收斂，對于重力影響，收斂會比較慢，這時候可以稍微調大一些容差加快收斂。
2：IDRFLG：這是控制動態松弛關閉和開啟的選項，選項很多，但我們基本只用0和1。0代表關閉，1代表開啟。
當你在進行顯示動力學分析時，開啟和關閉都對計算沒有影響，只要在載荷曲線（Curve）中SIDR一欄設置為1，就能強制開啟動態松弛。
當進行隱式分析時候，建議打開將IDRFLG設置為1（當然SIDR也要設置為1），當該選項為0的時候會出現一些奇怪的結果，大家可以去試一下。
2.消除動能
首先我更建議大家使用阻尼法來消除動能進行多次沖擊，如果大家想使用動態松弛的話，要注意幾點：
在進行第一次沖擊分析完成后，需要將*CONTROL_DYNAMIC_RELAXTION中DRTERM設置好結束時間，IDRFLG設置1開啟，
*CONTROL_TERMINATION中的ENDTIM設置為0，即不設置終止時間，這就代表著模型將開始動態松弛來消除動能。
展開 拉伸試驗CAE分析對比（涉及殘余應力映射、動態松弛） ￥15
本文主要講述：
1、拉伸試驗的CAE建模及分析，涉及樣片拉伸試驗仿真的約束和加載等；
2、通過關鍵字輸出拉伸試驗后樣片的殘余應力應變厚度變化等信息；
3、通過映射和動態松弛，將殘余應力應變引入試片拉伸分析，驗證加工硬化的影響。
拉伸試驗樣片基礎尺寸如下：
拉伸試驗CAE建模：
1、網格基本尺寸2mm，試片厚度1.2mm，材料B250P1。
2、左端對兩排單元的節點進行全約束（*BOUNDARY_SPC_option），右端對兩排單元的節點施加強迫運動（*BOUNDARY_PRESCRIBED_MOTION_RIGID）。
3、在試片中間建立彈簧單元來模擬標距，可以通過彈簧的變化量來計算應變。
工況一：加載端強迫位移15mm。
工況二：加載端強迫位移3mm，輸出dynain文件（包含殘余應力應變等）。
工況三：對拉伸試片映射工況二的殘余應力應變后，采用動態松弛，最后加載端強迫位移15mm。
以上僅作為學習研究的方法，涉及具體拉伸試驗對標等工作，需要做一定的調整。
展開 MATLAB/FORTRAN | 鍵基近場動力學（BBPD）動態松弛法實現準靜態單軸壓縮模擬（含預制裂隙），反力計算應力應變曲線 ￥119
程序采用經典的動態松弛算法（Dynamic Relaxation），將動力學方程轉化為解決準靜態問題的工具，模擬二維材料在單軸壓縮載荷下的響應及裂紋擴展過程。
準靜態模擬方案：利用動態松弛代碼，通過人為阻尼迭代，穩定求解準靜態單軸壓縮過程。
預制裂隙建模：代碼內置預制裂隙邏輯，用戶可根據需求自定義裂隙的位置、角度和長度，觀察裂隙對材料強度的影響。
鍵基 PD 理論基礎：嚴格遵循 BBPD 理論，涵蓋近場半徑（Horizon）確定、微模量計算及斷裂準則。
單軸壓縮工況：預設標準的單軸壓縮邊界條件，模擬材料在受壓狀態下的損傷演化。
應力應變曲線計算：通過反力計算試樣的應力應變曲線。
MATLAB/Fortran 編寫：代碼結構清晰，算法邏輯直觀，無須配置復雜的第三方環境，適合學習與二次開發。
損傷演化可視化：程序包含后處理模塊，可生成裂紋擴展路徑、損傷場分布圖。
參數可調：材料參數、幾何尺寸、離散間距及迭代終止條件均可靈活修改。
展開 LS-DYNA動態松弛實現應力初始化設置重力效應時，關于關鍵字設置的一些思考
k文件來自于論壇大神，原作者冰刀，Email： yj152052520@163.com QQ395550334
技術鄰原帖：https://www.yqgqt.org.cn/content/post/284766
drelax1test和drelax2test分別為動力松弛過程應力初始化和后續載荷施加過程，第一步重力加速度加載實現應力初始化，第二步重力加速度繼續加載；
drelax3test為動力松弛過程應力初始化和后續加載一步過程，下面我們看看它們的關鍵字設置具體區別在哪：
下圖為drelax1test計算得到的，用時37s，可以看到時間顯示是0，只有兩步
這里設置了IDRFLG=1，ENDTIM=0.0，*DEFINE_CURVE的SIDR=1，意味著該曲線只能用來應力初始化
運行drelax2test進行重啟動時，需要選取Implicit-to-explicit Sequential Solution，
然后計算時候會提示要求輸入重啟動文件，在命令框輸入m=drdisp.sif點擊回車
就會繼續計算，實現后面的計算
注意這里設置了IDRFLG=2 ENDTIM=0.03，*DEFINE_CURVE的SIDR=1，意味著該曲線只用于瞬態分析或其他應用。
下圖是Drelax3test計算得到的，用時49s，發現有后續的計算，也就是后續的重力加載，導致結果稍有差別，可以看到時間顯示是0.03，一共32步
注意這里設置了IDRFLG=1 ENDTIM=0.03，*DEFINE_CURVE的SIDR=2，意味著該曲線同時用于初始化和瞬態分析
注意，約束的是端面
使用上述的動態松弛法進行土壤重力的施加，得到的土壤重力分布如下圖所示
展開 
基于Ls-dyna電池包密封圈壓縮仿真過程詳解（附直播推薦）
三、仿真思路簡述
1）密封圈壓縮仿真需要考慮密封螺栓預緊過程，基于Ls-dyna的動態松弛關鍵字實現這一過程；
圖1 創建螺栓截面關鍵字
圖2 施加螺栓預緊力關鍵字
圖3 動態松弛關鍵字
2）密封圈類型分兩種，一種是完全不可壓縮材料（如硅橡膠），另一種完全可壓縮材料（如發泡材料），需要使用不同的材料本構進行模擬；
圖4 完全不可壓縮材料本構MATL27
圖5 完全可壓縮材料本構MATL57
3）密封圈壓縮仿真會出現密封圈網格畸變導致仿真報錯終止、主從剛度差異較大導致接觸穿透等棘手問題，本人總結了許多操作技巧解決了網格畸變、接觸穿透等問題：
a.根據密封圈的設計形式（開孔型、閉合型）分別采取不同的網格離散方式；
b.對于和密封圈接觸的零部件不能進行過度的前處理簡化；
c.密封圈材料不推薦使用全積分單元；
d.密封圈網格邊界應小于其它施壓零部件的網格邊界；
e.采用SPH方法模擬密封圈壓縮過程。
4）上述仿真思路將會在2024年11月26日技術鄰直播中做詳細展示，敬請期待，點擊下方圖片報名預約??
四、結果
圖6 壓塊重力預緊密封墊被壓縮(動態松弛法)
圖7壓塊受到推力與密封墊脫離接觸
圖8 動態松弛預緊轉推力仿真過程
五、重要說明
上述仿真思路是本人基于Ls-dyna官方學習資料和試用版軟件總結的，注意僅限于學習交流，請勿傳播，請勿商用，違者必究。
-end-
技術鄰簡介:
技術鄰專注于工科技術社區，從最早的CAE技術社區（中國CAE聯盟）發展而來，在CAE領域有20年的教學和咨詢服務經驗。
展開 FEA的核心思想-仿真時間步-隱式算法顯示算法
其中沖壓成型就是動態松弛法的主要例子。
c)在求解靜力學問題時，有時候將其看作動力學問題來處理而采用動態松弛法，這是顯式算法。Flac就是主要例子。
顯式算法
隱式算法
每步求解方法
矩陣乘法
線性方程組
時步穩定性
有條件
無條件
適用問題
動力中心差分法
動力動態松弛法
靜力動態松弛法
動力Newmark
動力靜態松弛法
基于LS-DYNA的過盈裝配的小案例 ￥20
以軸承和軸的熱裝配為例，過盈量為2mm，進行裝配
基本思路： 采用*CONTACT_SURFACE_TO_SURFACE_INTERFERENCE_ID進行動態松弛預應力裝配
動態松弛計算結果
動態松弛過程.gif
關鍵計算關鍵字 卡片已經列下下面，有一定dyna基礎的就應該自行實現！
在Dyna中對螺栓施加預緊力
就會生成關鍵字DEFINE_CURVE，但需要將SIDR改為1，表示激活動態松弛分析。
創建截面：
打開關鍵字DATABASE_CROSS_SECTION_PLANE，定義平面法向向量頭尾坐標，這個平面可以是矩形、圓形，也可以是無限大，這里定義一個圓形，輸入半徑為100.
這個關鍵字可以獲得該平面經過的Part的截面力或者截面扭矩。
創建Part集合：
點擊CreEnt>Set Data>*SET_PART，點擊Cre，框選螺栓，點擊Apply。最后在關鍵字中INITIAL_STRESS_SECTION，將這三個填入即可。
控制設置：
打control_DYNAMIC_RELAXATION，填寫收斂步數（NRCYCK）為100，收斂容差（DRTOL）默認，點擊Accept。
輸出設置：
為了創建動態松弛二進制輸出內容，點擊DATABASE_BINARY_D3DRLF，將CYCL改為1，其他都默認，最后保存。
后處理：
打開d3drlf，從動畫中可以看到截面周圍的單元有壓縮的趨勢，點擊Post > History > Element，制定其中某個單元，可以查看單元的應力。
展開 流固耦合方法進行隧道爆破開挖，怎么實現添加初始地應力？
在進行動態松弛方法模擬地應力時，為什么每次獲取動態松弛文件進行顯式計算都無法正常計算，實現這一操作的規范步驟是怎樣的，有無大佬指導一下，痛苦死了。
【12月19日項目懸賞】
【單號5850】
預算范圍：500以上
使用軟件：ABAQUS
需求描述：高阻尼橡膠支座的模擬，要求高阻尼橡膠的本構輸入，希望能有相應的本構論文做依據
立即搶單
【單號5860】
預算范圍：800
使用軟件：Hypermesh+LS-DYNA
需求描述：模型是滾石（剛體）撞擊鋼筋混凝土橋墩，重力是通過動態松弛方式施加的，現在的問題是自己把動態松弛添加上去了，計算中途停止，可能是動態松弛設置有問題，求助大神幫忙調一下；計算結果的能量平衡曲線也有一些問題（沙漏偏大），也需要調整一下。
立即搶單
【單號5861】
預算范圍：300-400
使用軟件： procast2016
需求描述：
任務：奧迪轎車鋁合金減震塔高真空壓鑄成型工藝過程模擬，目前有減震塔模型，需要設計出澆注系統和溢流槽，模擬過程主要考慮溫度，壓射比壓，壓射速度，真空度等的影響。縮松縮孔、冷隔現象改善以及澆注系統、溢流槽的優化。軟件：procast2016，時間：3-4周，預算可以詳談。
立即搶單
【單號5871】
預算范圍：2000
使用軟件：Fluent,SCDM,ICEM
需求描述：軸流風機氣動噪聲Fluent仿真
立即搶單
【單號5877】
預算范圍：1000-2000
需求描述：大斷面球墨鑄件中間心部出現球墨衰退，需要模擬鑄件的整個鑄造過程，分析中間部位液體的流動和凝固。
展開 【12月12日項目懸賞】
立即搶單
【單號5860】
預算范圍：800
使用軟件：Hypermesh+LS-DYNA
需求描述：模型是滾石（剛體）撞擊鋼筋混凝土橋墩，重力是通過動態松弛方式施加的，現在的問題是自己把動態松弛添加上去了，計算中途停止，可能是動態松弛設置有問題，求助大神幫忙調一下；計算結果的能量平衡曲線也有一些問題（沙漏偏大），也需要調整一下。
立即搶單
【單號5850】
預算范圍：500以上
使用軟件：ABAQUS
需求描述：高阻尼橡膠支座的模擬，要求高阻尼橡膠的本構輸入，希望能有相應的本構論文做依據
立即搶單
【單號5890】
預算范圍：500
使用軟件：ABAQUS
需求描述：ABAQUS在模擬抗浮錨桿分級加載的抗拔過程中，當分五級的時候施加荷載是750k可以收斂，當分七級施加荷載到1400的時候就不收斂了，是什么原因？應變3mm左右。變形縮放系數為1
立即搶單
下載“技術鄰”APP，或點擊鏈接查看所有派單： http://www.yqgqt.org.cn/requirement/more
注：目前手機不支持搶單，如需搶單請用電腦瀏覽器打開該網址，或打開技術鄰首頁，在首頁右側點擊“懸賞”進行搶單。
客服微信：jishulink888
展開 
近場動力學（PD）態基模型零能模式抑制算法對比：Silling、Li Pan、Wan Ji 多方法復現 ￥49
動態松弛法 (Dynamic Relaxation)：采用 Madenci 專著中的動態松弛策略，確保靜力學問題的準靜態求解穩定性。可視化后處理：內置 3D 散點云圖顯示、實時能量曲線監控（Energy Balance Check），數據可靠性高。結構清晰：包含 CommonFiles 庫調用、形狀張量（Shape Tensor）計算、變形梯度（Deformation Gradient）提取等核心 PD 算子。</p>
靜水壓力調試經驗貼
*EOS_LINEAR_POLYNOMIAL
$ EOSID C0 C1 C2 C3 C4 C5 C6
1 -1E5 0.00 0.00 0.00 0.40 0.40 0.00
$ E0 V0
2.5E5 1
$
調試結果
1.重力+無初始化+無質量阻尼
壓力云圖
壓力時程曲線
波形震蕩非常厲害
2.重力+無初始化+質量阻尼（采用4*pi/基頻周期）
質量阻尼采用*DAMPING_PART_MASS關鍵字
壓力時程曲線
波形有明顯改善，附k文件
hydrostastic_float.rar
3.重力+動力松弛靜水壓力初始化+質量阻尼（采用4*pi/基頻周期）
*DEFINE_CURVE開啟動態松弛，見k文件
float11.zip
，底部壓力正好是pgh
壓力時程曲線
4.重力+LOAD_DENSITY_DEPTH+質量阻尼（采用4*pi/基頻周期)
壓力云圖
壓力時程
展開 有限元與離散元。
顯式解法
顯式解法用于動力問題的求解或動態松弛法的靜力求解顯式算法無須建立像有限元法那樣的大型剛度矩陣只需將單元的運動分別求出計算比較簡單數據量較少并且允許單元發生很大的平移和轉動可以用來求解一些含有復雜物理力學模型的非線性問題時間積分采用中心差分法由于條件收斂的限制使得計算步長不能太大因而增加了計算時間。
隱式解法
隱式解法用于求解靜力問題的靜態松弛法隱式解法的動態松弛法式直接找導塊體失去平衡后達到再平衡的力位移關系建立隱式方法解聯立方程組并通過迭代求解以完全消除塊體的殘余力和力矩。
有限元方法的基本概念
將介質復雜幾何區域離散為具有簡單幾何形狀的單元而單元內的材料性質和控制方程通過單元節點的未知量來進行表達再通 過單元集成、外載和約束條件的處理得到方程組求解該方程組就可以得到該介質行為的近似表達。
有限元法的優點
可以用有限的、相互關聯的單元模擬無限的復雜體，無論多么復雜的幾何體都能用相應的單元簡化，從而建模分析計算出結果。使復雜的、感覺無處下手的工程問題簡單化，這是最大的優點。 有限元法采用矩陣形式表達，編程性高
對于線彈性問題當實際結構位移場函數連續光滑時能夠得到收斂解。
對于任意復雜結構理論上總是可以通過細分單元的方法獲得足夠近似的模擬。 長期大量工程應用積累了豐富的經驗。
有限元的插值是基于網格的，所以需要人為做好單元，這很耗時間，但是單元就好像人們修了路一樣，計算的時候可以節省很多時間，效率比較高。同時，這也是有限元法的一個缺點，大變形問題中的網格畸變問題，本質在于單元插值造成的
缺點：精確度浮動性比較大。基于建模的水平和邊界條件、載荷工況的模擬是否真實等等。
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
展開 離散元與有限元優缺點（附inp文件）
顯式解法
顯式解法用于動力問題的求解或動態松弛法的靜力求解顯式算法無須建立像有限元法那樣的大型剛度矩陣只需將單元的運動分別求出計算比較簡單數據量較少并且允許單元發生很大的平移和轉動可以用來求解一些含有復雜物理力學模型的非線性問題時間積分采用中心差分法由于條件收斂的限制使得計算步長不能太大因而增加了計算時間。
隱式解法
隱式解法用于求解靜力問題的靜態松弛法隱式解法的動態松弛法式直接找導塊體失去平衡后達到再平衡的力位移關系建立隱式方法解聯立方程組并通過迭代求解以完全消除塊體的殘余力和力矩。
有限元方法的基本概念
將介質復雜幾何區域離散為具有簡單幾何形狀的單元而單元內的材料性質和控制方程通過單元節點的未知量來進行表達再通 過單元集成、外載和約束條件的處理得到方程組求解該方程組就可以得到該介質行為的近似表達。
有限元法的優點
可以用有限的、相互關聯的單元模擬無限的復雜體，無論多么復雜的幾何體都能用相應的單元簡化，從而建模分析計算出結果。使復雜的、感覺無處下手的工程問題簡單化，這是最大的優點。 有限元法采用矩陣形式表達，編程性高
對于線彈性問題當實際結構位移場函數連續光滑時能夠得到收斂解。
對于任意復雜結構理論上總是可以通過細分單元的方法獲得足夠近似的模擬。 長期大量工程應用積累了豐富的經驗。
有限元的插值是基于網格的，所以需要人為做好單元，這很耗時間，但是單元就好像人們修了路一樣，計算的時候可以節省很多時間，效率比較高。
展開 