基于PANDA框架的非線性靜力學有限元

1、基于PANDA框架的非線性靜力學有限元論文導讀：基于PANDA框架。能夠分析千萬自由度規(guī)模的彈塑性靜力學問題。非線性求解策略。形成了面向對象有限元并行計算框架PANDA。并行計算，基于PANDA框架的非線性靜力學有限元。關鍵詞：PANDA，靜力學，非線性，有限元，并行計算1 引言特種武器結構復雜，在整個庫存到靶序列（Stockpile to TargetSequence，STS）全壽命周期內要經歷復雜嚴酷的載荷和環(huán)境條件，結構響應呈現出高度的材料非線性、邊界非線性和幾何非線性。為提高特種武器的設計、試驗和庫存維護水平，對武器結構在各種條件下響應的精細建模和分析至關重要，需要充分考慮結構的幾何細

2、節(jié)和物理內涵，所建立的有限元模型可達上千萬自由度規(guī)模乃至更高，而傳統(tǒng)的商用有限元程序由于國外對我國的出口限制，非線性有限元模型的分析規(guī)模被限制在幾百萬自由度以下，且計算周期較長，無法快速響應設計和維護的需要。為了提升特種武器的工程數值模擬能力，適應不斷提高的武器工程數值模擬需求，迎接和加速由現階段小規(guī)模低效率計算向大規(guī)模高效并行計算的轉變，2007年中國工程物理研究院啟動了院預研重大項目“武器工程大規(guī)模并行計算框架研究及基礎平臺開發(fā)”。該項目在已有源碼程序的基礎上，通過在有限元并行計算方法方面開展研究與軟件開發(fā)，初步形成了面向對象有限元并行計算框架PANDA，并基于PANDA框架初步開發(fā)了可應

3、用于部分靜力、振動、沖擊和傳熱武器工程問題求解的大規(guī)模有限元并行計算模擬程序。針對特種武器研制中的非線性靜力學有限元大規(guī)模精細分析需求，充分消化吸收開放源代碼的程序設計思想和技巧，基于PANDA框架，開發(fā)非線性靜力學有限元分析所需的單元類型、材料模型、非線性并行求解策略，集成大規(guī)模線性方程組并行求解算法，初步形成了可求解小應變、有限應變線彈性和彈塑性靜力學問題的非線性靜力學程序。懸臂梁彈塑性有限元分析模型達到了千萬自由度規(guī)模，并行求解時間低于一小時。本文介紹了基于PANDA框架的單元類型、材料模型、非線性求解策略設計，并初步驗證了非線性靜力學有限元并行計算程序的計算精度和千萬自由度規(guī)模分析能力

4、。2 基于PANDA框架的非線性靜力學有限元并行計算程序設計通過中國工程物理研究院的預研重大項目，采用面向對象、層次化、組件化的設計思想，對工程結構非結構網格有限元分析程序的基本數據結構、并行通信、求解控制等方面的共性和可重用部分進行抽象和程序實現，并集成了區(qū)域分割、解法器等服務組件，形成了面向對象有限元并行計算框架PANDA，提供經過系統(tǒng)規(guī)劃設計的應用程序開發(fā)接口，以提供服務的形式引導應用程序的設計和實現，初步建立了結構分析有限元并行計算應用程序的集成開發(fā)環(huán)境?？萍颊撐?，并行計算?；赑ANDA框架，結構分析有限元并行計算應用程序的開發(fā)工作變得較為簡單和高效，程序開發(fā)工作量大為減少。在PAN

5、DA框架既設的應用軟件架構下，應用程序開發(fā)者可以將精力集中到本應用程序獨有的個性部分，并充分利用框架中集成的經過充分驗證的高效解法器等服務組件，場和網格數據的組織、存儲和管理由框架負責，應用程序開發(fā)者無需關心其底層數據結構等實現細節(jié)。對非線性靜力學有限元并行計算程序的開發(fā)而言，通過繼承PANDA框架中場、節(jié)點、單元類型、材料模型、空間積分器、求解控制等共性部分，開發(fā)適用于非線性靜力學有限元計算的場、節(jié)點、單元類型、材料模型和非線性求解策略，使用由框架提供的輸入參數解析、并行通信、區(qū)域分割、線性方程組解法器等組件，就可較快速形成可求解大規(guī)模非線性靜力學問題的高性能有限元并行計算程序。2.1 單元

6、類型和材料模型PANDA框架中的ElementBaseT基類抽象出了各種單元類型的共性(屬性和操作)，在應用程序中，具體單元類型都可由該基類逐級派生并添加自身特有的屬性和操作得到，添加新單元類型較為方便。目前，在PANDA非線性靜力學有限元并行計算程序中已基于PANDA框架開發(fā)了結構分析中常用的八節(jié)點六面體小應變單元和更新拉格朗日有限應變單元，后者在每一載荷增量步對幾何構形進行更新?；赑ANDA框架，應用程序中具體材料模型在各個類繼承層次上定義，在程序結構設計時便于添加各種具體的材料模型。目前已實現了適用于小應變單元和更新拉格朗日單元的二維、三維各向同性線彈性和彈塑性材料模型。其中的彈塑性材

7、料模型是線彈性模型與J2各向同性屈服條件的組合，塑性變形階段的當前屈服應力由下述四種各向同性塑性硬化函數之一進行描述。（1）線性各向同性硬化函數(1)上式中為初始屈服應力，為硬化模量K（EET）/（E-ET）(2)其中，E為彈性模量，ET為割線模量?？萍颊撐?，并行計算。當硬化模量K為0時，采用線性各向同性硬化函數的小應變各向同性彈塑性材料模型就成為小應變理想彈塑性材料模型。（2）linear-exponentialsaturation各向同性硬化函數(3)上式中為初始屈服應力，為線性硬化模量。為saturation硬化模量，稱為saturation應變，是一特征等效塑性應變量，用于描述等效塑性

8、應變對saturation硬化的影響，的值越大，相同等效塑性應變時的當前屈服應力越小。隨著等效塑性應變值增大，linear-exponentialsaturation各向同性硬化函數逐步趨近于由表示的線性各向同性硬化函數（其初始屈服應力為）。在ANSYS中，Nonlinear Isotropic Hardening (NLISO)材料模型對塑性硬化規(guī)律的描述與此相似。（3）general cubic spline（立方樣條）各向同性硬化函數(4)上式中，n為用于定義樣條的點數（對應于樣條的n-1個子區(qū)間和2個端點外區(qū)域）。在每個子區(qū)間上，樣條具有下面的形式：(5)上式中，。從樣條點計算系數a(

9、i)的條件為：樣條通過這些數據，且樣條的一階導數在子區(qū)間之間連續(xù)。此外，樣條的端條件被規(guī)定為：，以滿足樣條點數據范圍之外的線性延伸?？萍颊撐?，并行計算。盡管不是必須和強制性的，人們一般選擇，以便用定義初始屈服應力。（4）power law（冪函數）各向同性硬化函數(6)由上式可知，初始屈服應力。2.2 非線性求解策略PANDA非線性靜力學有限元并行計算程序實現了標準牛頓法、帶線性搜索的非線性牛頓法、帶線性搜索的非線性預處理共軛梯度法等三種策略。它將整個載荷劃分為多個載荷增量（可根據迭代求解過程的情況進行自動增減），在每個載荷增量內進行反復迭代求解直至滿足某個收斂準則，而每一迭代過程中通過調用PANDA框架中集成的線性方程組解法器組件實現線性方程組的求解（線性方程組的求解可選用直接解法或迭代解法）。三種求解策略使用最多17個參數對非線性求解的迭代過程、精度、自動載荷增量步、線性搜索、預處理過程等進行控制，都使用基于殘余力范數的兩個收斂評估準則（絕對容差和相對容差，后者是兩次迭代的殘余力范數的比值），只要滿足其中一個準則，一次非線性迭代結束。3 測試算例為了初步驗證基于P