接觸厚度的案例
LS-DYNA中的接觸問題（五）（接觸參數）
8.接觸厚度的相關建議
有時可以通過修改SST和MST來減小接觸厚度，以消除初始穿透，但與進一步精細化網格相比，這只是一個不太好的替代方案。如果需要調整這兩個參數，我們建議根據實際情況來修改，若設置的數值太小，例如0.1mm，那么可能會出現問題，這是因為接觸厚度用來確定在釋放穿透節點之前所允許的最大穿透深度。通常情況下，增加接觸厚度可以避免接觸部件太薄所導致的計算失敗的情況。根據經驗，SST和MST不能小于0.6-0.7mm。
由于修改SST和MST的值會應用到整個模型中的所有接觸中，所以需要更謹慎的情況下，我們最好在關鍵字*PART_CONTACT中通過OPTT 或 SFT參數單獨對某一部件的接觸進行定義，尤其是在many parts of widely ranging thickness are included in a single contact這種情況下。
9.接觸厚度縮放系數
參數選項為：SFST 和 SFMT, card 3, *CONTACT_option)
除了直接修改接觸厚度以外，我們還可以利用SFST 和（或） SFMT作為縮放系數來修改接觸厚度。在關鍵字*SECTION_SHELL中，縮放系數和殼單元厚度相乘可得接觸厚度，默認值均為1.0。
10.接觸厚度縮放系數的相關建議
相關建議內容已經上述第8節，但還是要特別強調設置縮放系數時不要使接觸厚度小于0.6-0.7mm。
11. 粘性阻尼
參數選項為：V DC (Card 2,*CONTACT_option)
該參數默認為0。最初，接觸粘性阻尼用來抑制金屬成型模擬中接觸面法向震蕩，目前我們發現粘性阻尼對于碰撞和沖擊問題中的高頻震蕩也有很好的抑制作用。
展開 LS-DYNA中的接觸問題（一）（工作原理，接觸類型，單向接觸）
所以這種情況下我們建議使用自動接觸，因為這種接觸無需給定取向，也就是可以檢測來自殼單元每一側的接觸。在LS-DYNA中，如果*CONTACT關鍵字中存在AUTOMATIC字樣，這說明該接觸為自動接觸。與舊版本的接觸相比，自動接觸中使用的算法使其更適合于不連貫的網格。對于殼單元，自動接觸給定的接觸面平行于殼單元的中面，距離為接觸厚度的一半；同時在殼單元面的邊界上，接觸面以接觸面厚度的一半為半徑繞邊界翻轉，從而形成一個連續的接觸面。我們一般以殼單元的厚度為參考選定接觸面到殼中面的距離（此句存疑）。接觸定義過程中有許多可選參數供使用，用戶可以據此直接定義或者通過縮放來定義接觸厚度。若用戶沒有定義接觸厚度，則接觸厚度等于殼厚度（在單面接觸中，接觸厚度等于橋厚度和單元邊長中的最小值）。同樣的，在考慮梁的接觸是，梁單元的接觸面與梁中軸線的距離為梁橫截面的等效半徑。正是因為接觸面是由殼的中面和梁的中性軸偏置而來，所以在有限元建模中，為了表示殼的厚度和梁的橫截面尺寸而留下的殼部件之間與梁部件之間的間隙尺寸就顯得至關重要。不合適的間隙尺寸則會導致接觸面上出現初始穿透的現象。LS-DYNA將通過移動穿透的從節點到主面的方式消除初始穿透。并不是所有的初始穿透都需要消除，否則可能會導致不符合物理規則的接觸行為的發生。設置一個精確地初始幾何模型總是會花費很大時間。
在從節點開始被釋放，其接觸力設為0之前，LS-DYNA中的大部分接觸類型會針對最大穿透深度給定一個上限。這主要發生于自動接觸中，用于避免從節點穿過殼的中面時在反方向產生較大的接觸力。不同接觸類型對應的最大接觸深度以表格的形式列在了用戶手冊v.960的表6.1中。有時自動接觸并未起到作用，這是因為早在計算的開始就達到了接觸的起始點，這主要發生于接觸面中存在有非常薄的殼單元的情況下。
展開 abaqus顯式接觸厚度問題
我做了一個落錘沖擊梁的模型，落錘用的離散剛體，加工成殼，用動力隱式可以計算，但用動力顯式就報錯。
動力隱式容易不收斂，所以現在不知道在顯式中怎么改。
For two-dimensional models, if a non-unity thickness is specified for two-dimensional solid elements and these elements are involved in an interaction such as contact, the same thickness should be specified for the out-of-plane thickness of the corresponding surface under *surface interaction.
展開 LS-DYNA中的接觸問題：工作原理，接觸類型，單向接觸
所以這種情況下我們建議使用自動接觸，因為這種接觸無需給定取向，也就是可以檢測來自殼單元每一側的接觸。在LS-DYNA中，如果*CONTACT關鍵字中存在AUTOMATIC字樣，這說明該接觸為自動接觸。與舊版本的接觸相比，自動接觸中使用的算法使其更適合于不連貫的網格。對于殼單元，自動接觸給定的接觸面平行于殼單元的中面，距離為接觸厚
度的一半；同時在殼單元面的邊界上，接觸面以接觸面厚度的一半為半徑繞邊界翻轉，從而形成一個連續的接觸面。我們一般以殼單元的厚度為參考選定接觸面到殼中面的距離（此句存疑）。接觸定義過程中有許多可選參數供使用，用戶可以據此直接定義或者通過縮放來定義接觸厚度。若用戶沒有定義接觸厚度，則接觸厚度等于殼厚度（在單面接觸中，接觸厚度等于橋厚度和單元邊長中的最小值）。同樣的，在考慮梁的接觸是，梁單元的接觸面與梁中軸線的距離為梁橫截面的等效半徑。正是因為接觸面是由殼的中面和梁的中性軸偏置而來，所以在有限元建模中，為了表示殼的厚度和梁的橫截面尺寸而留下的殼部件之間與梁部件之間的間隙尺寸就顯得至關重要。不合適的間隙尺寸則會導致接觸面上出現初始穿透的現象。LS-DYNA將通過移動穿透的從節點到主面的方式消除初始穿透。并不是所有的初始穿透都需要消除，否則可能會導致不符合物理規則的接觸行為的發生。設置一個精確地初始幾何模型總是會花費很大時間。
在從節點開始被釋放，其接觸力設為0之前，LS-DYNA中的大部分接觸類型會針對最大穿透深度給定一個上限。這主要發生于自動接觸中，用于避免從節點穿過殼的中面時在反方向產生較大的接觸力。不同接觸類型對應的最大接觸深度以表格的形式列在了用戶手冊v.960的表6.1中。有時自動接觸并未起到作用，這是因為早在計算的開始就達到了接觸的起始點，這主要發生于接觸面中存在有非常薄的殼單元的情況下。
展開 
abaqus里的非線性薄層單元，零厚度cohesive單元，goodman接觸單元等的基本形式是什么？如何構建與應用？
在使用Abaqus，Comsol等軟件進行薄層區域的力學分析過程中，例如在研究水壓致裂、裂縫擴展，接觸粘結滑移的這類薄層力學性質時，我們經常需要采用應力-相對位移(σ-u)關系，而不是傳統本構描述的應力-應變(σ-ε)關系來描述，例如Abaqus里面的Cohesive單元，Goodman單元，以及Comsol里的彈性薄層（在后面我把這類單元統稱為增量非線性力學薄層）。這類單元厚度非常小甚至為0，薄層兩側的節點（單元）用一組力（應力）與相對位移的關系方程聯系起來，例如給出一個形式最為簡單的典型應力-位移方程
此方程描述了1,2,3方向（通常是法向和兩個切向）上相對位移與應力的關系，應力與相對位移呈線性關系，類似于“線性彈簧”。但是對于土-結構接觸、裂縫的張開閉合這類問題，線性方程已經不足以準確描述這些物理量之間的關系，這時就需要引入增量非線性方程來構建薄層單元。
引入增量非線性薄層的概念之前，首先介紹一下全量非線性薄層以理解非線性的概念，首先給出以下公式
這是一個全量非線性薄層，其非線性的表現可以用下面幾個例子體現，
對比①和②項，可以發現僅存在3方向上的位移變化的情況下，1，2方向上的力也會發生改變，體現了彈簧三個方向力學性質的非獨立性，對比①和③項，可以發現力的大小并不和位移大小成正比，也就是非線性特征。
所以對于增量非線性方程，就是把應力-位移關系方程寫成應力增量-位移增量的關系方程，例如
寫成微分形式的好處是，可以體現出應力路徑對位移結果的影響，也就是類似于“塑性”特征（所以所有的彈塑性本構也都是增量方程）。但是對于此類微分方程的求解，必須給定一個力的初始值。
展開 LS-Dyna的穿透問題解決建議
在LS-Dyna中初始穿透（initial penetration）會導致負的滑動界面能產生，同時由于并不是所有的從節點都會移到主表面上，存在的穿透節點將會導致不切實際的接觸行為。一般有以下三種方法可以解決初始穿透的問題：
（1）在建立模型時應當花費時間和精力避免有初始穿透，盡量保持接觸對中的接觸空隙（考慮殼單元的厚度），但對于復雜的模型，不可避免會出現初始穿透，此時可根據第一次遞交后程序給出的穿透信息（在 MESSAGE 和 D3HSP 文件中有詳細記錄），按照提示移動相關節點，調整計算模型，消除穿透；
（2）對于比較小的初始穿透問題，可以通過減小接觸厚度來解決，對應于 *CONTACT關鍵字中的控制參數SFST 和 SFMT。但由于縮小了接觸厚度，為保持接觸力的穩定，應相應增大罰函數剛度（控制參數SFS和SFM）。該方法只對很小的初始穿透效果好，對于大的初始穿透，可能會導致錯誤的結果；
（3）對于初始穿透問題，LSTC 公司在 LS-DYNA960 中增加了相關控制參數來處理該問題，在關鍵字
*control_contact 中有一參數 IGNORE，有多個選項可以控制、消除初始穿透，對應*CONTACT 關鍵字中也有相同的參數可以對單個的接觸對進行初始穿透處理。
展開 LS-DYNA中的接觸問題（六）（整車模型中的接觸問題）
定義面段為一個集合 by *SET SEGMENT
除此之外，*DEFINE_BOX可以定義一個三維的立方體box，限制在這一立方體之內部件或面段所參與的接觸可以設置為計算的起點，這么做可以縮小處理接觸問題的范圍，從未節省計算時間。
5.摩擦
當包含有多個部件的整車模型中定義了一個整體的基礎，那么設置一個同一的摩擦系數（可以為0）是可以接受的。但是當摩擦行為起到較大作用時，我們建議使用*PART_CONTACT來對各個部件單獨設置摩擦屬性。當且僅當*CONTACT中的FS設置為-1時，*PART_CONTACT中設置的摩擦系數會覆蓋掉*CONTACT中設置的摩擦系數。請注意，只有當衰減系數DC不為0時，動摩擦系數FD才會生效。
6.接觸厚度
可以通過修改*CONTACT中的全局SST和MST選項為適當值來減少初始穿透的影響。在單獨設置部件參數時，*PART_CONTACT中的 OPTT選項可以覆蓋上述兩個參數。用戶不應將接觸厚度設置為極小的值，因為這會導致計算錯誤；實際上對于厚度小于1mm的很薄的殼單元，有必要增加接觸厚度來避免計算失敗的發生。
（未完待續）
展開 參賽：LS-DYNA中的接觸界面模擬
2 如否，不作處理，搜索結束；如是，則在從節點與主面上的接觸點間附加
一法向接觸力 Fn。
3 處理摩擦力。
4 將接觸力 Fn和摩擦力 投影到總體坐標，組集到總體載荷向量中。
4 接觸類型
在具體介紹各種類型的接觸前，先闡述幾個基本的概念。
在殼單元中，自動接觸通過法向投影中面的1/2“Contact Thickness”來確定接觸面。這就是“shell thickness offsets”。接觸厚度可以在接觸的定義中明確指定。如果接觸厚度沒有指定，則等于殼的厚度（在單面接觸中，為殼厚度或單元邊長的最小值）。相同的，在梁的接觸中，接觸面從梁的基線偏置梁截面等效半徑距離。因此，在有限元幾何建模時，為考慮殼厚、梁截面尺寸必須在殼、梁的part間有適當的間隙，否則會有初始穿透現象發生（即發生不真實的接觸現象）。雖然LS-DYNA可以通過移動穿透的從節點到主面上來消除初始穿透，但是并不是所有的初始穿透都能檢查出。
DYNA中大多數的接觸有一個“極限穿透深度”，如侵徹超過這個深度則從節點被釋放，接觸力置為0。這主要用在自動接觸中，防止過大接觸力的產生而引起數值不穩定性。然而在有些情況下，因為這個閾值過早達到而使接觸失效（常發生在非常薄的殼單元中）。此時應采取的措施是放大接觸厚度因子或設置接觸厚度為大于殼厚度的一個值，或者改變接觸剛度的計算方法（如改為Soft=1）。
LS-DYNA中的接觸允許從節點與主段間壓縮載荷的傳遞。如接觸摩擦激活，也允許切向載荷的傳遞。Coulomb摩擦列式用來處理從靜到動摩擦的轉換，這種轉換要求一個衰減系數、靜摩擦系數大于動摩擦系數。
關于接觸搜索方法，這里僅給出幾個簡單的要點，詳細描述見Theoretical Manual of LS-DYNA。
展開 ANSYS LS-DYNA常見問題之Sliding Energy解答
圖6在message文件中查找初始穿透
l 調整關鍵字參數
■ 調整*CONTACT中的SFST和SFMT，即調整該接觸對從面和主面的接觸厚度的比例因子。該方法對微小的初始穿透效果很明顯，但是對于初始穿透值很大，效果不好，容易導致計算結果錯誤。
圖7 接觸卡片中的SFST和SFMT
■ 設置*CONTACT中的IGNORE參數，該參數默認值為0，即采用*CONT ROL_CONTACT中IGNORE的值。可以對每個接觸對進行修改，建議修改值為1或者2。
圖8 接觸卡片中的IGNORE
第二種情況的解決方法
眾所周知，ANSYS LS-DYNA存在2種方式設置接觸，即有Contact Region和Body Interaction。
■Body Interaction接觸即LS-DYNA接觸關鍵字*CONTACT_AUTOMATIC_SINGLE_SURFACE，主要用于整個模型的接觸或自接觸。
■ Contact Region 用于建立接觸對，當幾何模型導入到ANSYS以后，Mechanical默認對每個零件建立bond接觸。其輸出關鍵字為*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET ，其set類型為segment set。
對于接觸面和目標面這兩個segment進行接觸搜索時，在交界處的節點無法檢查出穿透，然而當節點滑到接觸厚度中時，就會發現穿透節點的存在，因此系統必定會給該節點施加一個力，把它拉回到接觸面上，在這個過程中系統對該節點做功，消耗其接觸勢能，能量曲線輸出表現為負的滑移能。
其解決方法之一是擴大主段的接觸面。
由*CONTACT中的關鍵字MAXPAR進行調整。
展開 LS-DYNA 中的接觸界面模擬(2)
6 接觸控制參數
LS-DYNA 提供了多個與接觸相關的控制參數。根據不同接觸問題的具體特點，設置不
同的控制參數，對提高“接觸模型”的精確性是非常必要的。
LS-DYNA 中的接觸控制參數可以在*Control_Contact、*Contact 或*Part_Contact 中設置，
而有些參數也可以同時在多個命令中設定。如一個參數在多個命令中設置，則這樣的設置有
一定的優先次序。*Control_Contact 對整個模型中的接觸提供一種“全局性”的“缺省”參
數設置；*Contact 對每個具體的接觸提供“局部”的參數設置，優先權較高；*Part_Contact
則為某個具體的Part 涉及的接觸提供最高級別的參數控制。
6.1 Thickness offset: Automatic, SLTHK(Card 1,*Control_Contact, Option Card A)
LS-DYNA 中非自動接觸類型：
1 *CONTACT_SURFACE_TO_SURFACE
2 *CONTACT_NODES_TO_SURFACE
3 *CONTACT_ONE_WAY_SURFACE_TO_SURFACE
利用參數SHLTHK 確定是否考慮“ 厚度偏置” （ 見下圖） ， 該參數可以在
*CONTROL_CONTACT 中全局定義，也可以在Optional Card B 中局部定義。如果
SHLTHK=0，不考慮厚度偏置，采用incremental search 方法來確定從節點最接近的主段；如
SHLTHK=1，考慮變形體的厚度偏置，但不考慮剛體厚度偏置；如SHLTHK=2，變形體、
剛體的厚度偏置都考慮。如SHLTHK 為1 或2，程序采用global bucket search 來確定接觸對。
展開 在LS-PrePost中查找并刪除初始穿透
3、穿透檢查
在 LS-PREPOST 中，您可以選擇通過切換“Select Parts”來檢查所選零件的穿透，也可以通過切換“*CONTACT”來檢查關鍵字文件中的現有接觸定義 ，然后從接觸列表中選擇要檢查的接觸。此方法可以檢查殼單元、梁單元和固體單元零件的穿透。
通過“Select Parts”進行穿透檢查
如果您尚未創建任何接觸定義，則此選項很有用。
■ 除非激活“Thickness”切換，否則此檢查的觸點厚度取自剖面卡數據。
■ 如果切換“Thickness”，所有殼單元和beam單元將得到用戶指定的接觸厚度，實體元件的厚度為零。
■ 執行的檢查遵循與 *CONTACT_GENERAL相同的規則，即將檢查node-to-surface、free_edge-to-free_edge 以及所有零件（包括自接觸）的所有部件邊緣穿透的節點。
按接觸進行穿透檢查交叉邊緣/ Crossed edges
從模型中定義的接觸列表中選擇一個接觸，將顯示接觸中所包含的單元，并檢查所有影響接觸厚度的參數（SST、MST、SFST、SFMT、SHLTHK、SLDTHK、SSTHK、OPTT、SFT、TH、TH_SF），以及要檢查哪些節點/單元的穿透度。
C2500車型發動機罩的穿透檢查
選擇發動機罩上的兩個部分并激活“Penetration”切換，將顯示與穿透距離成正比的白色方塊和箭頭的穿透節點。
如圖-14所示，檢查顯示出25個節點到曲面穿透和12個邊緣到邊緣穿透的報告以及最大穿透距離。
邊緣到邊緣穿透的節點報告方式與節點到曲面穿透的節點報告方式相同。通過切換“顯示：Penet”，可以直觀地關閉顯示的穿透。
展開 
LS-DYNA 中的接觸界面模擬(1)
2 如否，不作處理，搜索結束；如是，則在從節點與主面上的接觸點間附加一法向接觸力 Fn。
3 處理摩擦力。
4 將接觸力 Fn 和摩擦力 投影到總體坐標，組集到總體載荷向量中。
4 接觸類型
在具體介紹各種類型的接觸前，先闡述幾個基本的概念。
在殼單元中，自動接觸通過法向投影中面的1/2“Contact Thickness”來確定接觸面。這
就是“shell thickness offsets”。接觸厚度可以在接觸的定義中明確指定。如果接觸厚度沒有
指定，則等于殼的厚度（在單面接觸中，為殼厚度或單元邊長的最小值）。相同的，在梁的
接觸中，接觸面從梁的基線偏置梁截面等效半徑距離。因此，在有限元幾何建模時，為考慮
殼厚、梁截面尺寸必須在殼、梁的part 間有適當的間隙，否則會有初始穿透現象發生（即發
生不真實的接觸現象）。雖然LS-DYNA 可以通過移動穿透的從節點到主面上來消除初始穿
透，但是并不是所有的初始穿透都能檢查出。
DYNA 中大多數的接觸有一個“極限穿透深度”，如侵徹超過這個深度則從節點被釋放，
接觸力置為0。這主要用在自動接觸中，防止過大接觸力的產生而引起數值不穩定性。然而
在有些情況下，因為這個閾值過早達到而使接觸失效（常發生在非常薄的殼單元中）。此時
應采取的措施是放大接觸厚度因子或設置接觸厚度為大于殼厚度的一個值，或者改變接觸剛
度的計算方法（如改為Soft=1）。
LS-DYNA 中的接觸允許從節點與主段間壓縮載荷的傳遞。如接觸摩擦激活，也允許切
向載荷的傳遞。Coulomb 摩擦列式用來處理從靜到動摩擦的轉換，這種轉換要求一個衰減系
數、靜摩擦系數大于動摩擦系數。
關于接觸搜索方法，這里僅給出幾個簡單的要點，詳細描述見Theoretical Manual of
LS-DYNA。
展開 ANSYS LS-DYNA常見問題之Sliding Energy問題
對于接觸面和目標面這兩個segment進行接觸搜索時，在交界處的節點無法檢查出穿透，然而當節點滑到接觸厚度中時，就會發現穿透節點的存在，因此系統必定會給該節點施加一個力，把它拉回到接觸面上，在這個過程中系統對該節點做功，消耗其接觸勢能，能量曲線輸出表現為負的滑移能。
① 其解決方法之一是擴大主段的接觸面。
由*CONTACT中的關鍵字MAXPAR進行調整。MAXPAR默認值為0，即表示設置值為1.025，其LS-DYNA推薦值為1.2。
▲ 圖-9 接觸卡片中的MAXPAR
② 其解決方法之二是調換接觸算法。
ANSYS LS-DYNA在Mechanical中設定的bond接觸是*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET。在LS-DYNA接觸選項中帶OFFSET，其目的是把基于約束類型的接觸算法切換到基于罰函數的接觸算法。
因此，當出現滑移能異常，可以把*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET切換為*CONTACT_TIED_SURFACE_TO_SURFACE。
下文是一個簡單的對比，使用bond接觸的模型如圖-10所示，其采用不同的接觸算法的能量圖如圖-11所示。
展開 在LS-PrePost中查找并刪除初始穿透
穿透檢查
在 LS-PREPOST 中，您可以選擇通過切換“Select Parts”來檢查所選零件的穿透，也可以通過切換“*CONTACT”來檢查關鍵字文件中的現有接觸定義 ，然后從接觸列表中選擇要檢查的接觸。此方法可以檢查殼單元、梁單元和固體單元零件的穿透。
通過“Select Parts”進行穿透檢查
如果您尚未創建任何接觸定義，則此選項很有用。
■ 除非激活“Thickness”切換，否則此檢查的觸點厚度取自剖面卡數據。
■ 如果切換“Thickness”，所有殼單元和beam單元將得到用戶指定的接觸厚度，實體元件的厚度為零。
■ 執行的檢查遵循與 *CONTACT_GENERAL相同的規則，即將檢查node-to-surface、free_edge-to-free_edge 以及所有零件（包括自接觸）的所有部件邊緣穿透的節點。
按接觸進行穿透檢查
從模型中定義的接觸列表中選擇一個接觸，將顯示接觸中所包含的單元，并檢查所有影響接觸厚度的參數（SST、MST、SFST、SFMT、SHLTHK、SLDTHK、SSTHK、OPTT、SFT、TH、TH_SF），以及要檢查哪些節點/單元的穿透度。
展開 ABAQUS中的光滑粒子流體動力學 ( SPH ) 方法
相互作用
使用粒子定義的部件可以與其他有限單元定義的部件建立接觸關系。接觸關系類似于基于節點的面（與粒子相關聯）和基于單元的面（或解析剛體）之間的接觸關系。SPH支持通用接觸和接觸對。任何基于節點的面定義的接觸類型和準則都可用，包括cohesive behavior。支持定義不同的接觸特性。默認情況下，粒子不是通用接觸域中的一部分，這與一維單元（1-node elements）類似（例如point masses）。默認的粒子接觸厚度與在截面定義中指定的特征長度值相同；因此，為了實現接觸，粒子的接觸行為等效成一個球體，其半徑等于小立方體（與粒子相關聯的）的內接球的半徑。
對于PC3D單元，接觸厚度不能為0，否則計算難以收斂。推薦使用默認的或指定合理的接觸厚度。
允許在使用PC3D單元定義的不同部件間建立接觸關系。但是，只有當相互碰撞的SPH部件是由相同的流體材料組成時，接觸才有意義。例如水滴落入裝水的桶中。在與固體相關的應用中，比如子彈穿透裝甲板，建模時一個部件必須使用常規有限元定義。
粒子和歐拉體之間不能定義接觸。
Input 輸入文件使用：
使用下列語句來定義網格單元或解析剛體面與粒子定義的面之間的接觸:
*CONTACT
*CONTACT INCLUSIONS
node-based particle surface, element-based/analytical_surface
輸出
PC3D單元能夠輸出連續介質單元所有的與力學相關的參數：應力，應變，能量，狀態變量，場，用戶自定義變量。Abaqus/Explicit分析中的所有節點變量均可輸出。
限制
SPH法有以下限制:
· 該方法在變形不太嚴重,單元不扭曲時不如一般的拉格朗日有限元分析準確，在較大的變形時不如耦合歐拉-拉格朗日分析準確。
展開 