簡介面向工程分析人員和初學者的Workbench LS-DYNA技術(shù)培訓PDF系統(tǒng)講解ANSYS LS-DYNA顯式動力學分析的整體流程與核心功能。包內(nèi)為單個PDF文件大小3.74MB內(nèi)容依次介紹Workbench LS-DYNA前后處理器、Mechanical操作界面、集成LS-DYNA前后處理器選項并重點拆解前處理六大要素工程材料庫支持彈性、彈塑性、超彈性、狀態(tài)方程、失效模型與熱學模型等本構(gòu)定義ANSYS DesignModeler可完成二維草圖、三維參數(shù)化建模以及第三方CAD模型導入和修改接觸定義則涵蓋綁定、不分離、有摩擦、無摩擦等常見類型確保動態(tài)分析中部件相互作用準確模擬。教程還配有模型演示與操作演示直觀呈現(xiàn)金屬成型、沖擊碰撞、爆炸及流固耦合等典型場景的建模與求解過程。已有82人學習下載。讀者可借此掌握從CAD模型處理、材料參數(shù)設(shè)置、接觸配置到結(jié)果分析的一整套顯式動力學分析方法既適合新手快速入門也能幫助有經(jīng)驗的工程師系統(tǒng)查漏補缺。 翻開《Workbench LS-DYNA技術(shù)培訓》這份資料之前我以為它和ANSYS官方那些操作手冊差不多打開Workbench、拖一個LS-DYNA模塊、設(shè)幾項參數(shù)、點Solve完事。但真按這個思路走一遍才發(fā)現(xiàn)顯式動力學和平時做的靜力學完全是兩套邏輯。資料里講得輕描淡寫的地方恰恰是最容易把人卡死的環(huán)節(jié)。這篇文章就是圍繞培訓主線把PPT里一筆帶過、但真正決定你能不能算出來的那些東西展開聊透。1. 翻開培訓資料之前Workbench里的LS-DYNA和經(jīng)典版差在哪1.1 它到底是什么先把這個最基本的問題說清楚Workbench里的LS-DYNA本質(zhì)上是經(jīng)典LS-DYNA求解器外套了一層Workbench圖形界面。求解器還是那個求解器計算核心沒有縮水但前處理、后處理的方式和經(jīng)典版截然不同。經(jīng)典LS-DYNA的標準工作流是編輯K文件文本格式的關(guān)鍵字卡片→ 提交給求解器 → 用LS-PrePost看結(jié)果。Workbench的路子是在圖形界面里建模、劃網(wǎng)格、設(shè)接觸系統(tǒng)自動生成K文件算完之后再用Mechanical或者LS-PrePost做后處理。這意味著一個關(guān)鍵事實Workbench界面能操作的功能只是LS-DYNA全部能力的一個子集。LS-DYNA有上千個關(guān)鍵字卡片Workbench界面只暴露了最常用的那些。真正的高手在Workbench里建模最后也要回到K文件層面做微調(diào)。培訓資料里那句部分高級功能需通過關(guān)鍵字卡片實現(xiàn)說的就是這么回事。1.2 和Workbench里其他顯式模塊怎么選型Workbench里做沖擊、碰撞、爆炸這類問題會碰到三個容易混淆的模塊Explicit Dynamics、Autodyn、LS-DYNA。很多新手上來就蒙了不知道選哪個。我按自己的使用經(jīng)驗給個參照模塊核心求解器適合場景特點Explicit DynamicsAUTODYN或ANSYS顯式求解器通用顯式問題、入門學習和Workbench集成最緊密操作最簡單AutodynAUTODYN爆炸、沖擊波、流體-結(jié)構(gòu)耦合網(wǎng)格和算法選項豐富火工品領(lǐng)域常用LS-DYNA經(jīng)典LS-DYNA碰撞、跌落、金屬成形、多物理場材料庫最龐大汽車、航天工業(yè)認可度最高三者求解原理接近但材料庫、接觸算法、關(guān)鍵字體系差異很大。如果你面對的是汽車碰撞安全這類行業(yè)標準問題或者要和上下游交換K文件直接用LS-DYNA模塊。如果只是想快速驗證一個沖擊響應(yīng)的量級Explicit Dynamics更省事。一句話總結(jié)培訓資料選擇LS-DYNA作為主線是有道理的它是顯式分析這個圈子里通用性最強的工具多花點時間入門不虧。1.3 定位清楚了再動手你打算做什么類型的分析LS-DYNA能做的事遠比東西摔地上碎沒碎廣得多。按我自己的項目經(jīng)驗它的大頭應(yīng)用集中在幾類跌落與碰撞手機跌落、包裝防護、汽車碰撞、保險杠沖擊金屬成形沖壓、鍛造、輥彎LS-DYNA在這個領(lǐng)域有專門的求解功能爆炸與防護爆轟波傳播、防護結(jié)構(gòu)響應(yīng)但這里用Autodyn的也不少多物理場耦合流固耦合、熱力耦合比如油箱晃動、鳥撞發(fā)動機動手前先明確你做的是哪一類后續(xù)材料模型、接觸設(shè)置、控制卡片的選擇都會跟著變。培訓資料里一般以跌落或碰撞案例為主線因為它最能體現(xiàn)顯式分析的核心邏輯。2. 從幾何到結(jié)果跑通一次跌落分析的標準流程2.1 第一步就翻車單位制別等算完才發(fā)現(xiàn)錯了我在指導新人時發(fā)現(xiàn)跌落分析跑完發(fā)現(xiàn)結(jié)果差了6個數(shù)量級十有八九是單位制搞錯。LS-DYNA求解器本身不關(guān)心你用什么單位它只要求所有輸入保持一致的物理單位體系。K文件頭部有一段*UNIT_SYSTEM可以聲明單位制但很多人根本沒注意。工業(yè)界做顯式分析最常用的一套是毫米、千克、毫秒。對應(yīng)力和材料參數(shù)的單位就是牛頓、兆帕MPa、毫米/毫秒。這和Workbench默認的m-kg-s體系不一樣所以新建分析系統(tǒng)后第一件事就是把單位切到mm-kg-ms。很多人習慣把靜力學里算的速度直接填進去結(jié)果差了一千倍。比如一個部件以5 m/s的速度撞向地面在mm-ms單位制下應(yīng)該填5 mm/ms也就是5乘以1000除以1000數(shù)值上剛好是5。但如果是50 km/h的汽車碰撞50 km/h換算過來是13.889 mm/ms不是13889。數(shù)值對不上結(jié)果就是荒謬的。我自己的習慣是在工作表里先做一張單位對照表把長度、質(zhì)量、時間、速度、壓力全部換算好再動手建模。這個習慣救過我很多次你照做不會吃虧。2.2 材料定義Workbench面板只是冰山一角跌落分析里最常見的材料是金屬和塑料。Workbench的Engineering Data里LS-DYNA材料庫會以Density、Isotropic Elasticity這些熟悉的字段出現(xiàn)但對顯式分析來說靜力學設(shè)置里那些參數(shù)遠遠不夠。以彈塑性金屬材料為例LS-DYNA經(jīng)典的*MAT_PLASTIC_KINEMATIC也就是MAT_003需要密度、彈性模量、泊松比、屈服應(yīng)力、切線模量、失效應(yīng)變。Workbench界面里能填其中一部分但有些關(guān)鍵參數(shù)——尤其應(yīng)變率相關(guān)的Cowper-Symonds系數(shù)C和P——默認是不顯示的得切換到關(guān)鍵字視圖或者直接改K文件。做碰撞和跌落應(yīng)變率效應(yīng)不可忽略。鋼材在高速沖擊下的屈服應(yīng)力可能比靜力下高出50%以上。如果培訓資料里沒強調(diào)這一點等你算出來變形量偏大、和實驗對不上時回來查大概率是這個原因。碳纖維、泡沫、橡膠這類材料Workbench面板覆蓋得更少基本都要去K文件里寫材料卡片。所以培訓資料教你的核心能力不是在面板里填參數(shù)而是讀懂材料卡片里每個字段的含義。2.3 網(wǎng)格劃分顯式分析的網(wǎng)格和靜力學思路完全相反靜力學里網(wǎng)格局部加密天經(jīng)地義應(yīng)力集中區(qū)細劃遠處粗劃。但顯式分析最忌諱的恰恰是局部單元過小。原因在于顯式求解的時間步長受到最小單元尺寸限制。時間步長近似等于最小單元的特征長度除以材料聲速。鋼材聲速約5.9 mm/ms如果你在倒角處劃了一個0.2 mm的單元時間步長就大約只有0.03 ms。一個10 ms的跌落過程需要三十多萬步迭代每步都要完成接觸檢測這計算量足夠讓人懷疑人生。所以在LS-DYNA里網(wǎng)格規(guī)劃的目標是盡可能均勻而不是盡可能細。平均4 mm的網(wǎng)格局部最小能做到2 mm已經(jīng)是比較健康的狀態(tài)。想捕捉應(yīng)力集中靠的不是局部加密而是靠材料模型和接觸算法。培訓和實戰(zhàn)中我總結(jié)的網(wǎng)格檢查流程是先用自動網(wǎng)格跑一遍然后在求解信息里查看實際計算的時間步長再回到網(wǎng)格里查看最小單元出現(xiàn)在哪個位置。如果是尖角、倒角這種無關(guān)區(qū)域直接用虛擬拓撲或幾何清理把那個特征抹掉收益立竿見影。2.4 接觸、邊界和求解設(shè)置里的隱藏細節(jié)接觸設(shè)置是顯式分析和靜力學差異最大的地方之一。靜力學里接觸只要有就行顯式分析里接觸質(zhì)量直接決定結(jié)果可信度。最常見的選擇是*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE自動面面接觸。這個自動的含義是接觸方向由程序自動判斷可以在大變形過程中自動更新法向極其適合跌落碰撞這類拓撲劇烈變化的問題。什么情況要換TIEBREAK接觸焊接、膠接這種需要帶失效的粘接或者螺栓連接在碰撞中要模擬斷裂用TIEBREAK更合適。平面/殼體單元默認的接觸類型是自動單面接觸*CONTACT_AUTOMATIC_SINGLE_SURFACE它允許所有外表面互相接觸連哪兩個面可能接觸都不用指定。第一次做跌落時我選了它省掉了很多麻煩。接觸剛度也有講究默認的罰函數(shù)剛度夠用但接觸穿透太明顯時可以適當增大罰剛度因子不過調(diào)太大又會導致時間步長進一步縮小、計算變慢。這是典型的矛盾權(quán)衡。分析設(shè)置里最常改的三項終止時間、時間步長因子、結(jié)果輸出頻率。終止時間別設(shè)太長跌落分析前10毫秒內(nèi)沖擊響應(yīng)基本都出來了。時間步長因子默認0.9這是給計算留的安全余量正常情況下不需要動。結(jié)果輸出頻率決定d3plot文件的間隔輸出太密磁盤會爆太稀后處理又看不清細節(jié)一般按每毫秒輸出一個結(jié)果來控制。2.5 后處理別繞彎子LS-PrePost才是主力Workbench自帶的后處理能看應(yīng)力應(yīng)變、位移動畫但做LS-DYNA分析的人最后基本都去用LS-PrePost。原因很簡單它直接讀d3plot文件能看的東西比Workbench多得多。比如能量曲線——動能、內(nèi)能、沙漏能Workbench面板里不是不能看但要操作半天。LS-PrePost里加載History文件后兩條命令就能畫出來。做顯式分析能量曲線是判斷計算結(jié)果是否正常的第一道關(guān)口沙漏能占比超過總能量的5%到10%結(jié)果就值得懷疑了。實操中我也經(jīng)常在Mechanical里看應(yīng)力云圖在LS-PrePost里看動畫和能量曲線兩者互補。但如果你只肯用其中一個我建議你學LS-PrePost。3. K文件是LS-DYNA的靈魂關(guān)鍵字卡片沒那么玄3.1 一張K文件在說什么K文件用純文本描述整個模型。格式看起來神秘拆開其實就那么幾類段節(jié)點與單元*NODE、*ELEMENT_SOLID/SHELL材料與截面*MAT_xxx、*SECTION_SOLID/SHELL部件聚合*PART把材料、截面、單元關(guān)聯(lián)起來接觸*CONTACT_xxx控制參數(shù)*CONTROL_TERMINATION、*CONTROL_TIMESTEP、*CONTROL_HOURGLASS輸出*DATABASE_BINARY_D3PLOT決定結(jié)果文件寫什么邊界與載荷*BOUNDARY_xxx、*LOAD_xxx、*INITIAL_VELOCITYWorkbench的圖形界面操作最終就是在拼裝這些卡片。你設(shè)了接觸對它就生成對應(yīng)的CONTACT卡片你設(shè)了固定約束就生成BOUNDARY_SPC_SET。理解K文件結(jié)構(gòu)等于拿到了排查問題的底層視角。報錯信息報的是關(guān)鍵字名你不認識它就沒法排錯。3.2 關(guān)鍵字卡片的格式列對齊是刻在骨子里的約定K文件的卡片格式有個古老傳統(tǒng)每個字段占據(jù)固定的列位置字段之間用逗號分隔也行但老派做法是嚴格列對齊。LS-DYNA求解器解析時按列讀字段所以漏寫一個逗號或者錯位一個空格可能導致某個參數(shù)被錯誤地讀到另一個字段去求解器不一定報錯但結(jié)果會非常詭異。我遇到過一次詭異透頂?shù)那闆r*MAT_PLASTIC_KINEMATIC里失效應(yīng)變的大字段填錯了位置讀到了密度欄里結(jié)果材料密度變成幾十萬模型重得離譜跌落分析里那個零件紋絲不動。排查了兩天才發(fā)現(xiàn)是K文件里逗號數(shù)量不對字段整體錯位。經(jīng)驗教訓直接從Workbench導出K文件后除非你有十足把握盡量別手改卡片字段的位置。3.3 在Workbench里改K文件的正確姿勢培訓資料一般會告訴你兩種方式第一在Workbench的LS-DYNA模塊里插入Commands或使用Keyword Editor直接在界面里追加關(guān)鍵字。這種方式適合往里加獨立的關(guān)鍵字卡片比如*CONTROL_HOURGLASS、*DATABASE_GLSTAT。操作路徑一般是右鍵模型樹中的對應(yīng)環(huán)節(jié)選擇Insert Commands然后在編輯窗口里寫關(guān)鍵字。Workbench會在生成K文件時把你寫的卡片原樣帶上。第二導出K文件到本地用文本編輯器改完再提交。這種方式最靈活適合大規(guī)模修改但改完的文件脫離了Workbench的數(shù)據(jù)聯(lián)動以后再回到Workbench改幾何、改網(wǎng)格K文件里的手改內(nèi)容就會被覆蓋掉。這一點必須牢記。我自己的習慣是幾何和網(wǎng)格調(diào)整階段一直在Workbench里反復迭代不做手改模型最終凍結(jié)后再導出K文件做微調(diào)和精修。手改的內(nèi)容用版本管理起來每條改動都寫注釋防止過兩天自己都看不懂。4. 求解報錯不用慌時間步長、沙漏與負體積排查鏈路4.1 時間步長過小先查單元再查材料最后查接觸第一次提交求解最常見的報錯信息類似time step size is too small或者干脆卡在某一步計算極慢。排查順序按概率來排查網(wǎng)格最小單元找出全模型最小的那個單元在哪里。遇到過一個小伙伴因為幾何里留了一段0.01 mm的裝飾圓角時間步長直接給干到幾千分之一后面所有工作都白搭。清理幾何特征后計算速度立馬回到正常。查材料參數(shù)材料剪切模量或彈性模量填錯也能導致聲速異常。金屬材料還好說泡沫、橡膠類超彈性材料如果參數(shù)單位錯了時間步長會瞬間崩掉。查接觸剛度接觸罰剛度太高等效剛度也會拉低時間步長。有一種職業(yè)選手的做法叫質(zhì)量縮放允許求解器人為增加某些極小單元的質(zhì)量換取可接受的時間步長。設(shè)置*CONTROL_TIMESTEP里的質(zhì)量縮放開關(guān)和最大時間步長可以穩(wěn)定求解速度但代價是質(zhì)量增加會帶來慣性誤差。跌落問題里質(zhì)量縮放比例控制在百分之幾以內(nèi)工程上可接受但具體多少需要算完對比能量曲線來判斷。4.2 沙漏能超標單點積分單元的老毛病LS-DYNA默認使用單點積分單元一點積分優(yōu)點是計算快、能承受大變形缺點是會產(chǎn)生沙漏模式——一種非物理的零能變形模式看起來像網(wǎng)格在顫抖或扭曲但對應(yīng)單元的應(yīng)變能是零。判斷標準是看沙漏能和總能量的比值。如果沙漏能占到總能量5%到10%以上計算結(jié)果基本就不能用了。控制沙漏的手段按優(yōu)先級排列細化網(wǎng)格最直接有效但代價是計算量增加改用全積分單元*SECTION_SOLID里的ELFORM2對八面體單元有效但計算更貴增加沙漏控制剛度*CONTROL_HOURGLASS里的IHQ參數(shù)默認1是粘性沙漏控制還有基于剛度的4、5、6選項工程上的實用建議是先細化局部區(qū)域網(wǎng)格再看沙漏能曲線是否有回落。不要在細網(wǎng)格都沒試過的情況下就硬加沙漏控制剛度那相當于靠數(shù)值手段掩蓋模型缺陷。4.3 負體積大變形分析最讓人頭禿的錯誤負體積negative volume指單元在變形過程中節(jié)點順序翻轉(zhuǎn)體積變成負值通常出現(xiàn)在固體單元極度壓縮或扭曲時。報錯信息直接告訴你negative volume in solid element #XXXXX。排查鏈路通常是先定位是哪個部件、哪個區(qū)域的單元出問題然后在后處理里回放動畫看那個位置在出錯前經(jīng)歷了什么變形。多數(shù)原因是局部網(wǎng)格被壓扁或過度剪切本質(zhì)是單元已經(jīng)無法承受這種變形。治本的辦法是重新劃分網(wǎng)格讓該區(qū)域單元形態(tài)更勻稱或者檢查材料模型是否合理——比如不考慮失效的材料在超高壓下會無限壓縮最終必然負體積。我一向建議做跌落分析時把可能接觸地面區(qū)域的網(wǎng)格先細化好不要在報錯后再補。還有一類負體積是接觸穿透嚴重導致的solve里接觸檢測和網(wǎng)格運動互相牽扯局部單元被擠爆。這種時候同時調(diào)整接觸罰剛度和網(wǎng)格尺寸往往能解決。4.4 排查方法論一次只動一個變量這是我?guī)ы椖孔顖猿值囊粭l紀律。顯式分析的報錯原因經(jīng)常互相耦合時間步長小了可能是網(wǎng)格問題也可能是接觸問題沙漏超標可能是單元積分問題也可能是載荷太快。新人常見做法是看到報錯同時把網(wǎng)格加密、把沙漏控制打開、把接觸剛度調(diào)大然后求解跑通了皆大歡喜但根本不知道是哪個變量起作用了。模型稍微一改又崩了一切從頭再來。正確做法是先停一步只改一個變量記錄結(jié)果再改下一個。雖然看起來慢但每一步都幫你積累對模型行為的理解。算顯式分析調(diào)試時間占整個項目周期的一半以上很正常這部分投入不能省。5. 培訓資料之外的實戰(zhàn)提醒與工具鏈5.1 環(huán)境配置的坑裝好不等于能算Workbench安裝是很多人的第一道坎?？吹接信笥严螺d安裝包后卡在許可證配置或者裝完打開就閃退。以我的經(jīng)驗這類問題多數(shù)出在兩方面一是許可證服務(wù)沒起來或者和版本不匹配二是顯卡驅(qū)動太老導致圖形界面渲染異常。Workbench的圖形界面和臺式機獨顯的兼容性總體沒問題但筆記本雙顯卡集顯獨顯環(huán)境下偶爾會出現(xiàn)UI卡頓甚至閃退更新的顯卡驅(qū)動能解決大部分。還有就是跑LS-DYNA求解時別把核心數(shù)設(shè)置成物理核心數(shù)的兩倍給超線程求解器是科學計算負載超線程不但不加速反而可能由于內(nèi)存帶寬爭用略微降速。5.2 方法論層面的三個建議如果培訓結(jié)束你要獨立做分析下面三條個人經(jīng)驗供參考第一從小模型起步。別上來就仿整車碰撞。先做一個單零件的跌落驗證材料、接觸、邊界都正常再逐步加復雜度。一個能在幾分鐘內(nèi)算完的模型是你調(diào)試參數(shù)的最佳試驗場。第二用實驗數(shù)據(jù)校準模型參數(shù)。仿真算出來的結(jié)果和實驗一定存在偏差。材料參數(shù)、摩擦系數(shù)、接觸剛度不同數(shù)值組合出來結(jié)果差異很大。靠譜的做法是拿到試驗數(shù)據(jù)后反向調(diào)校材料模型的關(guān)鍵參數(shù)直到仿真結(jié)果和實驗在關(guān)鍵響應(yīng)點上吻合。第三做好模型版本管理。顯式分析調(diào)試過程中K文件、配置、網(wǎng)格版本經(jīng)常反復迭代沒有版本管理改著改著就不知道自己改了什么。不需要多復雜的工具哪怕每天壓縮存檔一次也比沒有強。我在踩了不少坑后才養(yǎng)成的習慣希望你能早一點用上。本文還有配套的精品資源點擊獲取