ALE方法及其附屬的流固耦合方法，旨在模擬一系列流體與固體間具有較大動量和能量轉換特點的瞬態工程問題。 LS-DYNA ALE 多材料單元模式允許同一網格中多種流體共存。 進而它所帶的流固耦合算法可分析固體結構與各單個流體之間的相互作用。 這種優點使得它被廣泛用于分析多種工程領域的問題。

ALE/FSI 組件可很好地解決攜帶較大動量或能量密度的流體撞擊，侵入結構這一類工程問題。 例如，爆炸，油箱液體晃動， 容器跌落，飛鳥撞擊， 彈藥撞擊，飛行器濺落，輪胎打滑等。

新近開發的ALE 本質邊界條件（ALE ESSENTIAL BOUNDARY） 功能可大降低在處理流體與剛體間耦合的模擬時間。 這一功能將在包裝，石油，化工，制造行業內得到利用，來模擬管道流， 樹脂成型等問題。

ALE網格的任意性 (Flexible ALE Mesh Motion)

ALE網格與通常有限元網格不同，它既非隨物體本身變形（拉格朗日方法）， 也非在空間內保持不動（歐拉方法）。它可隨特定物理問題，采用自己特別的特殊移動方式。 這種任意性使得ALE網格可以根據問題特質而靈活改變，從而減少所需單元和計算時間。例如，在模擬彈頭飛行并穿透裝甲時，ALE網格可跟隨彈頭移動；又如，模擬高能炸藥起爆時，網格可隨炸藥膨脹而擴張。

ALE方法的簡便性 (Compatible ALE Solver) 

LS-DYNA ALE 采用算子分裂法（Operator Splitting）分別處理擴散項與遷移項。 這樣，一個時間步長內，單元會經歷一個正常的拉格朗日時間步，以及一個額外的平流時間步（Advection timestep）。

在拉格朗日時間步內，我們求解動量守恒方程并更新加速度，速度和位移；在平流時間步內，我們計算質量，速度，歷史變量的遷移項。 由于拉格朗日時間步的求解利用LS-DYNA求解器（Solver），ALE求解器與通常的結構求解器（Structure Solver）直接兼容。 相對于傳統計算流體力學方法，LS-DYNA 有著所解方程數少，能直接兼容LS-DYNA材料庫，求解簡便的優點。

ALE多流體模擬 (ALE Multi-material Capability)

LS-DYNA ALE 使用界面重建法（Interface Reconstruction）來構造不同流體間的物質界面。  這使得同一網格中多流體的計算成為可能。 例如，求解油箱晃動問題時，ALE網格中可分別定義汽油和油箱內和油箱外空氣。這樣，我們不僅可以模擬汽油對油箱的沖擊，還可以研究油箱內空氣負壓對結構的損傷。 在爆炸防護問題中，ALE網格中可定義炸藥，土壤和空氣，這樣土壤和炸藥對結構的沖擊也可以被加以考慮，而不象傳統的CONWEP方法僅僅對結構施加空氣中的爆炸波。相比其它界面描述方式如LEVELSET方法等，LS-DYNA ALE的界面重建法有著即無需求解額外方程，也無需儲存額外變量的優點。 LS-DYNA ALE多材料單元類型11（SOLID ELEMENT TYPE=11)  支持多流體計算。

緊密的ALE流固耦合 (Tightly-coupled ALE FSI) 

使用罰函數接觸法（Penalty Coupling Method）。 ALE/FSI 流固耦合算法專門針對固體與ALE多流體之間的接觸，相對于其它Loosely Coupled Methods，流體與固體間的信息交換簡便易行，也無需多步迭代，而且易于并行計算。

ALE/FSI 已成功解決多種工程領域內的流固耦合問題。在這些流固耦合問題中，流體一般都帶有大的動量或能量密度。這些大動量或能量密度會在很短時間內被傳遞到固體結構上，造成結構的大變形或損傷。ALE/FSI在處理以瞬態，高能，大動量，結構大變形為特征的，例如爆炸，油箱液體晃動， 容器跌落，飛鳥撞擊， 彈藥撞擊，飛行器濺落等工程問題上，有著較強的優勢。

經濟的ALE 投影功能(Cost-Efficient ALE Mapping) 

處理固體大變形的ALE單材料單元方法 (ALE Single Material Element to Cure Element Distortion) 

ALE方法在處理固體大變形方面，有著很大的優點。通常的拉格朗日有限元方法，可以用加高網格密度或采用高階單元的方法處理網格畸變（Mesh Distortion)。 但這種做法只能減輕網格畸變而無法有效阻止。 ALE方法可采用網格舒緩法（Mesh Smoothing）降低畸變處的網格密度而有效防止網格畸變，從而保證有限元計算的進行。 LS-DYNA ALE單材料單元類型5（SOLID ELEMENT TYPE=5) 采用網格舒緩法。與其它大變形處理方法相比，ALE單材料單元所需計算時間較少，內存較小，處理很大和超大變形的能力更強。

Structured ALE: 針對規則網格的快速求解器 

        在處理這些問題時，我們發現，絕大多數情形下，ALE模型采用規則正交網格(rectilinear)，也稱為IJK網格。網格不一定非要均勻，因為在關鍵區域通常我們會細化網格；但它一定是正交的，規則的立方體網格。這種特殊網格的幾何信息其簡單，在善加利用的情況下，可以很大程度上降低算法的復雜度，從而達到減少運算時間，降低內存需求的目的。

        同時，舊有的ALE求解器開發之初是用來解決固體大變形問題的。這類問題中，網格隨物質邊界變形而移動，而固體也只是用單材料單元來模擬。雖然LSTC的開發者不斷擴展原程序來支持多材料，支持網格移動等等。但舊有的算法和邏輯遠非很好。

        還有，最近十年以來，ALE模型的單元個數增長，由百萬到現在的千萬量級。輸入文件變得龐大，編輯文件耗時很長也非常不便。而對規則網格而言，我們完全可以根據用戶提供的簡單幾何信息，由程序本身自行創建網格。而省去用戶創建網格->程序讀入的麻煩，同時我們也節省了大量讀寫操作帶來的運算時間和內存需要。

輪胎土壤耦合

S-ALE與ALE: 相同與不同

        理論方面, S-ALE求解器與ALE完全相同，但S-ALE在實現上是獨立于ALE而全新開發的。網格的簡單與單一性也使得程序變得簡潔與易維護。簡潔的程序又反過來提高了運行的效率。綜合目前為止的大部分算例結果以及用戶的反饋，運行時間大致減少20%到40%之間。并且，借助于MPP算法的全新設計，MPP的scalability極大提高，大型算例一般可以保有0.9的加速率到大約400核。

        輸入卡片方面，S-ALE使用*ALE_STRUCTURED_MESH卡片創建網格并進行運算。其它卡片，包括流固耦合，使用原有的ALE卡片。這里額外說明，S-ALE中網格運動由*ALE_STRUCTURED_MESH中的原點平動以及局部坐標系轉動控制，不再使用*ALE_REFERENCE_SYSTEM卡片。

        S-ALE,引入兩個新概念：材料PART和網格PART，并且施加一個特殊的限制：ALE的多材料必須與材料PART一一對應。這種做法可以使得建模過程成為標準化的三步操作，清晰而簡單，避免舊有ALE方法中PART的雙重屬性帶來的概念混淆與困擾。具體作法請參閱http://ftp.lstc.com/anonymous/outgoing/hao/sale/

S-ALE: 不斷進步

        自問世以來，S-ALE不斷加入新的功能。例如漸變網格密度；多網格支持；ALE輸入文件->S-ALE自動轉換；運行過程中刪除網格；多網格合并等等。我們也不斷添加新的算例(http://ftp.lstc.com/anonymous/outgoing/hao/sale/models)，來幫助用戶建模與參考。S-ALE將不斷進步，以其能得到中國用戶的認可。

*CONSTRAINED_BEAM_IN_SOLID：嵌入結構的約束

        加筋混凝土是一種典型的嵌入結構，鋼筋嵌入在混凝土中來提供加強作用。進行有限元離散時，一般對鋼筋使用梁單元，對混凝土使用固體單元。然后一種作法是創建單元網格時，強迫固體單元與梁單元共享節點。這樣對建模有很高的要求，而且對于復雜一點的加筋情況就無法處理了。所以目前一般的作法是在獨立建模的鋼筋和混凝土之間施加約束。這一約束是施加在速度與加速度兩者之上的，它保持系統的動量守恒。這里請注意，約束法一般的情況下是不能滿足動能守恒的。時間步長越小，約束越頻繁，系統損失的動能就會越多。所以對于瞬時效應強烈的問題，例如沖擊載荷下的加筋混凝土，我們要采用懲罰法。

軸向滑移和脫粘(debonding)模擬　*CONSTRAINED_BEAM_IN_SOLID一個特別的功能是模擬debonding現象。這需要我們只在垂直于梁的方向施加約束，而對軸向約束加以放松(CDIR=1)。Debonding過程中的切向滑移力可以用User function(AXFOR=-N)或者User subroutine(AXFOR=N,N>1000)來施加。值得一提的是，這個功能在模擬往復作用例如地震載荷中比以往用*CONTACT_1D的作法更為準確。

梁單元端點間額外耦合點 理想情況下，在梁單元穿過的每個固體單元中，我們希望都至少有一個耦合點。這樣每個固體單元都可以和梁單元有相互作用，而不會出現“漏掛”的情形從而導致固體局部應變過大或過小。單用梁單元端點（節點）不能保證這一點，所以我們用NCOUP=N在梁單元中間施加額外的耦合點(coupling points)。這里需要特別加以說明的是，舊有的*CONSTRAINED_LAGRANGE_IN_SOLID, CTYPE=2中有類似功能(NQUAD)，但它的算法和實現上都有錯誤，會導致系統動能人為增加。

四面體和五面體單元 四面體和五面體單元在舊有的*CONSTRAINED_LAGRANGE_IN_SOLID中是當作degenerated 六面體來處理的。這樣會造成四面體和五面體的各個節點權重不同，結果有很大偏差，并且會隨網格的改變而劇烈變化。CBIS對四面體和五面體分別采用了它們自己的形函數，各節點擁有正確的權重，從而保證結果的正確性。

本質邊界條件 一般而言，約束法不一定能保證速度場重新分布后，還能滿足原來的速度邊界條件。例如，如果梁端點在固體單元中心，而這個固體單元的右側節點有速度為零。普通算法下，梁端點在這些節點處的動量投影不為零，造成這些節點的速度邊條被破壞。*CONSTRAINED_BEAM_IN_SOLID對于固體上的速度邊界條件做了特殊處理，可以在保證動量守恒的情況下，滿足速度邊界條件。

支持自適應網格(r-adaptivity) 這個功能主要是為模擬纖維復合材料的制造過程－－在可塑性強的材料中添加纖維，然后在模具中施壓后最終成型。纖維由梁單元模擬，在用固體單元離散的基體中沿切向滑動并產生滑移力?；屏突凭嚯x的關系使用AXFOR選項給定。但是現在的問題是固體單元在模壓過程中的變形很大，如果不采用某些特殊處理的話，很快問題就會因為單元畸變而無法繼續求解。一種解決畸變的方式就是在固體變形較大時重新構筑網格，稱為自適應網格(remeshingadaptivity, r-adaptivity)。這種方法在實現中，有一個停止運行－》寫出網格和歷史變量－》重構網格并投影歷史變量－》繼續運行的過程。CBIS在這一過程中需要記錄并重投影與切向滑移有關的歷史變量。

懲罰法　約束法對于瞬時效應較小的工程問題表現很好，但對于應變梯度大的問題會有動能丟失的毛病。這是因為在動量在梁單元和固體單元來回分配的時候會造成平滑效應(smoothing)。其中一個直觀表現就是采用不同的時間步長，問題的結果相差很大。檢查能量平衡時會發現時間步長越小，約束次數越多，動能丟失越大。為解決這種問題，我們提供了耦合的懲罰法選項－－*CONSTRAINED_BEAM_IN_SOLID_PENALTY。簡單而言，就是在梁和固體間構造彈簧，它們之間的相對運動會造成彈簧拉伸從而產生懲罰力。這樣，動能會轉化為彈性能儲存在彈簧中并最終釋放，達到既保持能量平衡，又達成運動一致的目的。

纖維強化塑料成型

Curtsey Shingo Hayashi, JSOL Corp.