不穩定性的案例
CFD理論|流動不穩定性
導讀：介紹兩種流動不穩定現象：開爾文-亥姆霍茲不穩定、瑞利-泰勒不穩定。
流動穩定性
流動穩定性（hydrodynamic stability） 流動受初始擾動后恢復原先運動狀態的能力。外界的擾動如果會自動衰減，原先的流動便是穩定的；外界的擾動如果會發展，并轉變為新的流動狀態，這就是 流動失穩現象。
流動穩定性理論研究流體運動穩定的條件和失穩后流動的發展變化，包括轉捩為湍流的過程。
層流到湍流的轉捩，一般始于失穩。但隨著某流動參數（如雷諾數）的逐漸增大，流動失穩后也有可能過渡為另一種更為復雜的層流，最后再失去層流的規律性，轉捩為湍流。
本文介紹的兩種不穩定現象是屬于有一個明確界面的穩定性問題，
開爾文-亥姆霍茲不穩定
開爾文-亥姆霍茲不穩定性（英語：Kelvin–Helmholtz instability，名稱來自開爾文男爵和赫爾曼·馮·亥姆霍茲）是在有剪力速度的連續流體內部或有速度差的兩個不同流體的界面之間發生的不穩定現象。
KH不穩定性廣泛存在于高能量密度物理、地球和天體物理、慣性約束聚變、燃燒、玻色-愛因斯坦凝聚、石墨烯等領域。充分發展的KH不穩定性導致了星際颶風、星系旋臂、太陽風與地球磁層相互作用中大規模旋渦結構的形成；另一方面，被顯著抑制的KH不穩定性有助于高準直、高長寬比、高穩定性超聲速天體射流的形成。
比如說風吹過水面時，在水面上表面的波的不穩定。而這種不穩定狀況更常見于云、海洋、土星的云帶、木星的大紅斑、太陽的日冕中
瑞利-泰勒不穩定
當重流體處于輕流體上方時，如果界面無限平整且不存在擾動，則該流體系統處于不穩定平衡狀態。由于在自然界中擾動的不可避免性，即便是原本無限平整的界面在重力作用下也會發生失穩。
展開 ：通過構筑熱力學不穩定性和動力學穩定性實現彈性體修復后變強
總結：作者利用離聚物以及聚電解質類材料的熱力學不穩定的特質，從分子設計上引入大位阻，制備了在常溫下動力學穩定但熱力學不穩定的材料。當動力學穩定性被熱或者力刺激破壞后，熱力學不穩定性使得材料中未配對的離子進行配對進而形成更多更大的聚集體，這些聚集體作為更強的物理交聯點賦予材料更強的力學性能。從而真正實現了像生物材料一樣的超量恢復行為。
該工作被發表在Materials Horizons雜志上（Materials Horizons, 2021, DOI: 10.1039/D1MH00638J），第一作者為博士生彭燕，通訊作者為吳錦榮教授。該工作由國家自然科學基金（51873110）和四川省科技計劃項目（2021JDJQ0018）。
原文鏈接：
https://pubs.rsc.org/en/content/articlelanding/2021/mh/d1mh00638j#!divAbstract
展開 Nat.Commun：利用納米尺度的力學不穩定性實現超低剛度金屬
【小結】
這篇文章揭示了在納米尺度內力學不穩定性對超低剛度材料性能的影響。通過材料在相變過程中經歷的力學不穩定性和對超材料的分子動力學模擬的手段，在準靜態條件下得到了剛度低于2GPa的超金屬。研究結果表明，在準靜態條件下，超材料通過一系列的手段得到了完全致密、全面強度納米級的超低剛度。
文獻鏈接：Harnessing mechanical instabilities at the nanoscale to achieve ultra-low stiffness metals（Nat.Commun, 2017，DOI:10.1038/s41467-017-01260-6）
本文轉自材料牛
展開 【CAE案例】細長結構在軸向流體作用下的變形
這種不穩定性，是在某些條件下，系統在擾動的影響下自發地離開其平衡狀態，包括靜態不穩定性和動態不穩定性兩種類型。
求解這類不穩定性問題的一般方法，通常采用瞬態方法，而這里將采用靜態穩定性的直接研究。采用簡化的不可壓縮、無湍流、無粘度的勢流模型來描述流動的不穩定性，并在一個簡單案例中通過實驗和數值結果進行驗證。實驗裝置的實物圖和實驗模型示意圖如圖2和圖3所示。
圖1 流體狀態的變化
圖2 實驗裝置實物圖
圖3 實驗模型的示意圖
02 解決方案
求解不穩定性問題主要分為兩步，第一步如圖4所示，循環迭代1到N階模態，按照各階模態陣型定義幾何模型，生成對應的網格，計算幾何模型周圍的勢流，進而獲得流體對板的壓力；第二步如圖5所示，基于固有模態投影壓力載荷，建立靜態平衡，尋找系統的特定值，從而確定不穩定性的閾值，即臨界失穩速度Ucrit。
圖4 不穩定性問題第一步求解流程
圖5不穩定性問題第二步求解流程
本研究中針對拉普拉斯類問題進行有效求解，基于Python內置命令實現自動計算與連續網格，采用固體力學仿真軟件求解細長結構在軸向流體作用下的變形問題。
具體而言，先進行流體邊界的計算，如圖6所示，在一組節點或網格上施加邊界條件。采用軟件內置的operator THER_LINEAIRE進行勢流計算，完成溫度場T與勢能場φ的轉換，將流體速度場等效為流體分量場FLUX，FLUY，FLUZ，再使用operator FORMULE定義壓力公式，使用operator CALC_CHAMP計算流場各處的壓力，使用operator CREA_CHAMP存儲各節點上的壓力。
展開 
本構模型(Constitutive Models)選擇
第一個問題是物理不穩定性。如果巖土材料破壞時出現軟化行為，就會出現物理不穩定性。當物理不穩定性發生時，部分材料儲存的能量會以動能形式釋放，這會導致解答可能無法收斂。
第二個特點是非線性材料的路徑依賴性。大多數地質力學系統有無數的解決方案可以滿足描述系統的平衡、兼容性和本構關系。為了找到一個"正確"的解決方案，必須指定一條路徑。例如通過爆炸的開挖是突然進行的，那么解答可能會受到慣性效應的影響，從而引入材料的附加破壞；但如果開挖是逐漸進行的，那么這種情況就不會出現。數值求解方法應該能夠適應不同的加載路徑以正確應用本構模型。
第三個特點是應力-應變響應的非線性。這包括彈性剛度和強度包絡對約束應力的非線性依賴，包括極限破壞后的行為，它根據應力水平而改變特性(例如在拉伸、非約束和約束條件下的不同破壞后響應)。數值求解方法需要能夠適應這些不同形式的非線性。
上述在地質力學數值模擬中所面臨的困難---物理不穩定性、路徑依賴性以及極端非線性本構模型的模擬可以通過使用顯式動態的求解方法來解決，例如FLAC3D和3DEC。顯式動態的求解方法允許數值分析以真實的方式跟蹤地質力學系統的演變，而不必擔心數值不穩定的問題。顯式動態求解方法建立了完整的動態運動方程，即使模擬的物理系統不穩定，例如邊坡的突然倒塌，但數值解是穩定的。當系統中的一些應變能量轉化為動能時，因為包括了慣性項，所以能夠耗散動能而達到靜力平衡。顯式動態求解方法直接模擬了這一過程。(zone mechanical energy active/clear)
相反，如果求解方法中不包括慣性項，那么必須使用一些數值程序來處理物理不穩定性。即使該過程能成功地防止數值不穩定，但所應用的路徑也可能不是真實的。系統的物理演化方式會影響到解的結果。顯示動態的求解方案可以遵循物理路徑。
展開 氣動聲學建模指南
后兩種類型的波依賴于背景流場的速度進行對流，而不以聲速傳播。聲波在傳播過程中可以（通過反應項）與流體相互作用，聲能量可以傳遞給聲學模式，并從聲學模式傳遞給渦旋模式和熵模式?？刂品匠讨械姆磻検且痤愃频牧鲃勇晫W耦合的原因。這是因為渦旋波和熵波對背景流場的解產生了非聲學（類似于 CFD）擾動，因此在某種程度上，反應項模擬的是 CFD 和聲學之間的線性相互作用。
許多氣動聲學公式忽略了反應項，因為它們也是 開爾文-亥姆霍茲不穩定性 的產生過程的背后原因。這些過程很難進行數值處理。另一方面，如果忽略這些項，則不能對聲音衰減和放大進行精確建模。反應項完全包含在線性納維-斯托克斯接口中。
不穩定性增長在 COMSOL Multiphysics 中有兩種處理方式。我們可以通過選擇頻域而非時域公式來解決隨時間增長的不穩定性。如果不能正確地對渦旋模式進行解析，可能出現空間不穩定性，這時可以利用伽遼金最小二乘法穩定方案有效地進行處理。
根據利用線性納維-斯托克斯方程模擬的不同應用，我們可能需要解析聲學、粘性和熱邊界層。如果存在無滑移和等溫邊界條件，則需要在固體表面上創建上述邊界條件，從而對振蕩流進行解析。通常情況下，在大型模型（同邊界層厚度相比）中，沒有必要考慮邊界層的損耗。在液體中，我們通常也可以忽略熱邊界層，但是氣體中一定要添加。通過在壁邊界條件下勾選滑移或絕熱選項，可以忽略這兩種效應。
值得一提的是，我們還能在背景流場和聲場之間創建另一個間接耦合。當聲波在涉及湍流背景流場的區域中傳播時會衰減。只要將 CFD RANS 模型的湍流粘度耦合到聲學模型中，即可將衰減效應引入模型。比如說，在分析存在流動現象的消聲器系統的傳輸損耗時，衰減效應是一個重要的考慮因素。
展開 光伏發電消納是什么意思？如何消納？
1.光伏電池組成的系統
這種技術主要是針對光伏發電的不穩定性和間歇性，使用多組光伏電池并聯的方式來提高光伏發電的穩定性和可靠性。通過并聯多個光伏電池組，可以降低單個光伏電池組產生的電壓和電流的不穩定性，從而提高整個光伏發電系統的穩定性和可靠性。
2.電池儲能系統
這種技術主要是利用電池儲能系統對光伏發電進行平衡控制。在光伏發電的高峰期，將多余的電能存儲到電池中，以備低谷期間使用，從而實現光伏發電的平衡消納。
3.智能控制系統
這種技術主要是利用人工智能、大數據和云計算等技術手段，對光伏發電進行智能化控制和管理。通過對光伏發電系統進行實時監測和預測，以及對電力系統負荷的實時調度，實現光伏發電的高效消納和穩定運行。
二、光伏消納的影響
l 電力系統穩定性：由于光伏發電的不穩定性和間歇性，當光伏發電容量達到一定規模時，會對電力系統的穩定性造成一定的影響，容易引起電壓波動和頻率變化等問題。
l 電力系統安全性：光伏發電在電力系統中的接入需要遵守一定的技術規范和標準，以保證電力系統的安全運行。如果不合理接入，可能會對電力系統的設備和人員安全造成威脅。
l 能源結構優化：光伏消納的推廣可以促進能源結構的優化和調整，降低對傳統能源的依賴，實現能源的可持續發展和環保。
l 經濟效益：光伏消納可以降低電力系統的運行成本，提高電力系統的效益，同時也可以為光伏發電提供更好的發展環境和市場前景。
以上就是小編分享的全部內容了，如果還想了解更多內容，或者對光伏電站開發感興趣，可以繼續關注鷓鴣云，也可以私信評論小編。
展開 轉動部件和靜止部件之間相對運動引起的不穩定相互影響：
1. 潛在相互影響：從上游和下游傳播來的壓力波動引起的流動不穩定性
2. 尾跡效應的影響：從上游向下游傳遞的尾跡引起的流動不穩定性
3. 沖擊波影響：對于亞音速/超音速流動，沖擊波沖擊下游扇葉引起的不穩定性
MRF和MPM都忽略了這些相互作用，因此只限于用在這些交互作用很弱的流動問題上。
而Slide Mesh模型對這些交互作用的預測比較準確。
III. FLUENT處理旋轉動力機械問題的4種模式：
1. SRF (Single Rotating Frame)：整個流場都以同一旋轉參考系為參考（最簡單的模式，不再贅述）
2. MRF(Multiple Rotating Frame)：有一個以上的參考系，忽略流域之間的相互影響。在流動穿過轉動區域（流進然后又流出轉動區域的外邊界）的地方可能得出誤導性的結果。
3. MPM (Mixing Plane)：對于有多組扇葉而且相鄰兩組扇葉的扇葉數量不一樣的旋轉動力機械，不能用MRF模型建立周期性流場，此時可用Mixing Plane模式建立周期性流場。
4.SMM (Slide Mesh) : Slide Mesh模型考慮了所有流域間相對運動引起的交互影響，可以得到比較精確的結果。但只能用于瞬態計算，耗費計算資源較大。
Ⅳ. MPM (Mixing Plane)模型
1. 為什么要用MPM模型？
渦輪機械，不管是軸向式還是離心式，通常包括1個或多個stage，每個stage由定子葉片組和轉子葉片組組成。并且相鄰葉片組的葉片數量也是不相同的。每組葉片的periodic angle不同的時候，無法用MRF模型建立周期性流場。這種情況可以用Mixing Plane模型來簡化計算。
展開 轉動部件和靜止部件之間相對運動引起的不穩定相互影響：
1. 潛在相互影響：從上游和下游傳播來的壓力波動引起的流動不穩定性
2. 尾跡效應的影響：從上游向下游傳遞的尾跡引起的流動不穩定性
3. 沖擊波影響：對于亞音速/超音速流動，沖擊波沖擊下游扇葉引起的不穩定性
MRF和MPM都忽略了這些相互作用，因此只限于用在這些交互作用很弱的流動問題上。
而Slide Mesh模型對這些交互作用的預測比較準確。
III. FLUENT處理旋轉動力機械問題的4種模式：
1. SRF (Single Rotating Frame)：整個流場都以同一旋轉參考系為參考（最簡單的模式，不再贅述）
2. MRF(Multiple Rotating Frame)：有一個以上的參考系，忽略流域之間的相互影響。在流動穿過轉動區域（流進然后又流出轉動區域的外邊界）的地方可能得出誤導性的結果。
3. MPM (Mixing Plane)：對于有多組扇葉而且相鄰兩組扇葉的扇葉數量不一樣的旋轉動力機械，不能用MRF模型建立周期性流場，此時可用Mixing Plane模式建立周期性流場。
4.SMM (Slide Mesh) : Slide Mesh模型考慮了所有流域間相對運動引起的交互影響，可以得到比較精確的結果。但只能用于瞬態計算，耗費計算資源較大。
Ⅳ. MPM (Mixing Plane)模型
1. 為什么要用MPM模型？
渦輪機械，不管是軸向式還是離心式，通常包括1個或多個stage，每個stage由定子葉片組和轉子葉片組組成。并且相鄰葉片組的葉片數量也是不相同的。每組葉片的periodic angle不同的時候，無法用MRF模型建立周期性流場。這種情況可以用Mixing Plane模型來簡化計算。
展開 案例21-環肋圓柱的非線性穩定屈曲與后屈曲分析
施加非線性穩定
有兩種方法可用于控制穩定力：
• 應用阻尼系數（STABILIZE，，damping，，，）
• 應用能量比（STABILIZE，，energy，，）
選擇使用哪種方法，然后為能量或阻尼因子選擇正確的值，都是非常重要的決定。要使用的最佳選項和值取決于不穩定性的類型、模型中使用的單元類型和大小以及加載步的結束時間和子步數。在大多數情況下，決策是基于試錯糾正過程。目標是以最小的穩定力實現收斂，該穩定力可以通過能量比或阻尼系數進行控制。
阻尼法
第一次試驗通常使用能量選項進行，因為：
• 它具有可使用的特定值范圍（介于0和1之間），以及
• 運行后，在下一子步開始時計算阻尼系數，為可使用的阻尼系數提供參考值（如果選擇阻尼作為穩定方法）。
嘗試的第一個能量值提供了0.1e-2的阻尼因子，這有助于收斂，但不會產生顯著的屈曲。然后，阻尼值逐漸減小到0.1e-5，即使在每種情況下都實現了收斂，但由于屈曲導致的坍塌也沒有發生。0.1e-6的值不會導致收斂，應用阻尼系數法的試驗停止。
觀察到的是由于該問題固有的局部化現象。當應用阻尼時，所有單元都使用指定值。當施加的阻尼值太大時，結構上施加的穩定力太大，因此系統剛度太大，容易收斂，不會產生太大變形。相反，當施加的阻尼太小時，不穩定單元不會從足夠的穩定力中受益，因此解會發散。
結論是，阻尼選項對于涉及顯著局部屈曲的問題并不理想。
能量法
由于該問題的特點是局部不穩定性而非全局不穩定性，因此穩定化的能量方法更為有用。能量法對不同的單元使用不同的阻尼值，因此不穩定性大的單元具有較高的阻尼值而不穩定性小的單元具有較小的阻尼值。因此，可以補償不穩定性而不會使系統剛度變得過大。
展開 幾篇論文amesim
基于AMEsim 的液力變矩器進出口
定壓閥動態仿真研究
AMESim_用于不穩定性分析的重型汽車動力傳動建模（英文）
基于AMESim的液力變距器進出口定壓閥動態仿真研究.pdf
AMESim_用于不穩定性分析的重型汽車動力傳動建模.pdf
核聚變箍縮反應堆
總部位于西雅圖的Zap能源公司的FuZE-Q反應堆計劃于2022年年中完工，這個反應堆不需要昂貴復雜的磁線圈，而是由機器沿著高導電等離子體柱發送電流脈沖，產生的磁場可同時限制、壓縮和加熱電離氣體。這種方法被稱為Z箍縮（Z-pinch），這樣命名是因為電流沿著三維網格的第三軸或Z軸箍縮等離子體，其產生能量的裝置有可能會比笨重的托卡馬克裝置以及當前正在開發的激光聚變設備更簡單、更小巧、更便宜。
Z-箍縮等離子體過去一直受到不穩定性的困擾。在沒有完全均勻擠壓的情況下，等離子體會在幾十納秒內起皺、蜷縮和倒塌，時間太短，無法產生有效的電量。
Zap能源公司的方法叫做“剪切流穩定”，其通過調整等離子體沿立柱的流動來克服不穩定性。這種設計用速度更快的等離子體流將等離子體包裹在立柱中心軸附近。離子體就像在高速公路中間車道上行駛的一輛輛汽車，由于兩側車道呼嘯而過的車輛高速而密集，因此無法變換車道。與之前的Z箍縮結構相比，這種設置保持聚變反應等離子體被集結和被壓縮的時間更長。
Zap能源公司的研發總監本?萊維特（Ben Levitt）說：“我們認為我們的反應堆是最便宜、最緊湊、最可擴展的解決方案，也是實現商業聚變發電的最短路徑?！比R維特預測，到2023年中期，Zap將達到Q=1的科學盈虧平衡點，即聚變原子釋放的能量等于創造聚變條件所需的能量。它將成為第一個達到該點的聚變項目。
由于聚變能源研究長期以來都未能兌現諾言，這類說法值得懷疑。但Zap能源公司在令人生畏的陡峭技術曲線上迅速攀登，令人印象深刻。這家初創公司成立于2017年，由華盛頓大學聚變Z箍縮實驗（FuZE）研究團隊派生而來。第二年，該公司就生產了第一次批聚變反應產品。在公司成立之前，該大學團隊曾與美國勞倫斯?利弗莫爾國家實驗室的研究人員合作。
展開 清潔能源 | 如何利用仿真技術應對氫燃料挑戰
與典型碳氫化合物相比，氫燃料具有更高的火焰速度（高8倍）、更低的點火能量要求（低15倍）以及更大的可燃性限值（4%-70%）。氫的這些特征，為設計基于氫燃料及氫混合燃料的能源轉換系統提供了機遇，但同時也帶來了挑戰。
比如，氫的特征有助于提高效率和燃燒穩定性。然而，氫更高的火焰速度和更大的可燃性限值為回火及其它安全相關問題帶來了關鍵挑戰；氫火焰更高的火焰溫度，則為氮氧化物和金屬保護帶來了挑戰。由于氫的路易斯數（熱擴散率與質量擴散率之比）較低，導致其存在顯著的差異擴散效應，而這是引起燃燒不穩定性的主要因素。差異擴散效應將導致局部等效比變化，從而導致沿火焰前緣的反應速率發生變化。因此，大規模采用氫作為更清潔的燃料的進程，取決于解決與回火、氮氧化物排放和燃燒不穩定有關問題的速度。
一些研究小組正在研究如何利用實驗室測試和仿真來緩解這些挑戰。大田韓國科學技術院和Ansys正在制定計算流體力學（CFD）方法和最佳實踐，以利用大渦模擬仿真（LES）預測氫甲烷混合火焰的火焰結構。
韓國科學技術院燃燒動力學與診斷實驗室開展的研究
KAIST CDDL正在研究重型燃氣輪機燃燒室、飛行器發動機加力燃燒室及雙推進劑液體火箭發動機的低頻及高頻燃燒不穩定性。其目前的研究工作側重于：
高頻燃燒不穩定性的觸發，以及相關復雜模態動力學和多千赫茲橫向熱聲波動的機理細節。
氨/氫基無碳燃氣輪機燃燒的激光診斷測量和數值仿真。
詳細表征先進軸向燃料分級燃燒系統中，橫向反應射流的自激勵動力學所產生的拓撲特征。
氫燃料及氫混合燃料的Ansys仿真方法
過去的燃燒模型和最佳實踐，大多是多年來基于碳氫化合物燃料發展而來的，并且有大量實驗數據作為支撐。不過，這些模型和最佳實踐還需要基于氫燃料和氫混合燃料進行研究驗證。
展開 消費后回收材料(PCR)的高附加值應用開發
圖1：評估PCR材料應用可行性的技術路線圖
基于這樣的背景，誠模精密的材料應用研究中心開發了PCR材料應用可行性評估方法（圖1），利用此方法篩選出合適的PCR材料，可以在不升高模溫的條件下，有效改善射出件外觀（凹痕、表面浮纖等），縮短成型周期，降低生產成本；利用玻璃纖維和PCR材料的協同效果，可以在不降低材料力學強度的前提下改善射出制品外觀，進一步提升PCR的創新應用附加價值，并成功應用于薄壁產品的生產（圖2）。結合Moldex3D模流仿真分析技術，不同批次PCR原料的結構和性能變化的不穩定因素將得到有效管理和控制。模流仿真過程會適應PCR的結構變化，快速得到適合每一批次PCR的成型工藝條件，進而指導現場射出，有效降低現場射出成型的調機成本，大大提升PCR產品的質量。
圖2：PCR加纖用于薄壁產品開發
從材料分子結構考慮，熔融指數是影響制品外觀的重要因素，因此，我們嘗試利用消費后回收材料(PCR)來改善射出件外觀。塑膠材料經過一段時間的使用，可能會發生材料的降解、老化，因此平均分子量可能會減小，從而使得PCR材料再加工的時候流動性變好，進而有可能改善外觀，但同時力學強度也會受到影響。因此，合理、有效、快速地評估PCR材料的結構和性能變化，是找到PCR材料高附加值應用的關鍵。
目前PCR材料的原料來源極不穩定，盡管各家原料商都有提供PCR材料的物性表，但依然無法保持原料中所使用的PCR材料結構和性能一致。因此，我們需要借助實驗室分析資源，對于PCR材料的不穩定性做出有效評估，并將此不穩定性納入模具設計和成型工藝的考量之中，通過Moldex3D模流仿真分析，搭建材料不穩定性與射出成型的關系。
結語
如何評估和使用PCR材料，提升PCR材料應用的附加價值，是實現產業綠色發展的重要一步。
展開 abaqus顯示動力學應用-鉆頭切削巖石
abaqus擁有強大的顯示動力學求解能力，應用abaqus的Explicit做了個牙輪鉆頭切削巖石的案例，總結以下幾個遇到的問題：
1.不開多線程可以正常求解，一開多線程就報錯？
在顯示動力學接觸設置中，abaqus軟件默認的是動態接觸算法（Kinematic contact algorithm）,當開啟多線程時就會由于求解速度過高而產生計算的不穩定性，而該算法的接觸約束嚴格性很高，因此當遇到求解不穩定時就會產生報錯從而導致計算終止。由于罰函數法的接觸約束嚴格性要低于動態接觸算法，因此改為罰函數法（penalty contact method）即可。
2.切削中鉆頭和巖石發生穿透？
在切削仿真中鉆頭和巖石間的接觸壓力、接觸剛度和許用穿透量之間的平衡被打破?？梢约毣佑|區域網格; 修改接觸剛度;用軟接觸代替硬接觸
3.仿真中的求解不穩定性問題？
由于abaqus explicit的接觸算法對接觸面的類型有較嚴格的限制,而切削仿真又是一個高度非線性求解過程，這些都會導致求解的不穩定?？梢圆捎眉毣W格、調節增量步長、采用ALE技術、在接觸中引入阻尼等來完成的。
4.巖石的切削形態如何控制？
建議失效選用位移方式，合理選擇失效數值
5.求解時間較長？
顯式動力學是采用顯式算法進行動力學方程的求解，顯式算法最大優點是有較好的穩定性，不存在隱式算法中的收斂性問題。顯式動力學最適合發生在短時間，幾毫秒內的事件或更小時間。持續1秒以上的事件可以模擬但是需要較長的時間，通過諸如質量縮放和動態松弛之類的技術可用于提高模擬效率減少計算時長。
求解效果圖如下：
展開 