ABAQUS分析中ALE

1、ALE方法詳解及各選項(xiàng)卡參數(shù)意義與設(shè)置(ABAQUS/Explicit )為了方便理解，先整體介紹一下ALE網(wǎng)格自適應(yīng)方法的基本過程一個完整的 ALE過程可以分為若干個網(wǎng)格remesh子過程，而每一次remesh的過程可以分為一下兩步：1生成一個新的網(wǎng)格( create a new mesh)利用各種算法以及控制策略生成一個良好的網(wǎng)格2環(huán)境變量的轉(zhuǎn)換(advection variales)也就是將舊網(wǎng)格中的變量信息利用remapping技術(shù)轉(zhuǎn)換到新網(wǎng)格中，也有不同算法，其中包括靜變量(應(yīng)力場，應(yīng)變場等)的轉(zhuǎn)換與動變量(速度場，加速度場等)的轉(zhuǎn)換生成的新網(wǎng)格成功與否則在于對這兩步的控制首先來看第

2、1個大部分-對整個adaptive remesh過程控制：包括兩個方面：一個是對adaptive remesh過程的算法控制，另一個對adaptive mesh過程強(qiáng)度的控制1. 對adaptive remesh的算法控制算法控制包括兩部分算法控制，其一為網(wǎng)格算法控制，其二為變量轉(zhuǎn)換算法控制(1) 網(wǎng)格重畫算法控制：在ABAQUS中是如何生成新網(wǎng)格的呢？即使用網(wǎng)格掃掠技術(shù)(mesh sweep technique)，每sweep 一次，生成一套新的網(wǎng)格。但 是當(dāng)你使用的算法不同時，sweep出來的網(wǎng)格也是不同的，打個不是很恰當(dāng)?shù)谋确剑河貌煌墓ぞ咦鐾患|西，做出來的質(zhì)量與精度會不一樣，同樣，

3、用不同的算法來sweep網(wǎng)格，得到的網(wǎng)格質(zhì)量也會不一樣。在ABAQUS顯示模塊中，sweep算法用英語來說就是 mesh smoothing method,有三種算 法來sweep網(wǎng)格，如下所示：1) 體積算法(volume smoothing )該算法十分健壯，為默認(rèn)算法，再絕大多數(shù)情況下適用2) 拉普拉斯算法(laplacian smoothing )耗費(fèi)資源最少的算法，能力一般，作用與體積算法類似(一階算法，類似于求平均值)，對于曲率比較高的曲線曲面邊界時，效果不是很理想3) .等位算法(equipotential smoothing )比較復(fù)雜的算法，是基于拉普拉斯算法的解之上的算法，

4、對曲率較大的曲線曲面邊 界效果較好，在節(jié)點(diǎn)被非結(jié)構(gòu)化網(wǎng)格包圍時，次算法為推薦算法，若節(jié)點(diǎn)被結(jié)構(gòu)化 網(wǎng)格包圍，其效果與體積算法類似。三種算法可以結(jié)合適用，利用權(quán)重值來定義，需要記住的是，三種算法各占的權(quán)值加 起來必須等于1。在ABAQUS CAE 中，在step模塊中定義選項(xiàng)卡如下所示：TB tdrt ME Adaptive Me*h Controls(2)變量轉(zhuǎn)換算法控制:在ABAQUS中是如何將舊網(wǎng)格中的環(huán)境變量轉(zhuǎn)換到新網(wǎng)格中的呢？即使用remapping技術(shù)，對于靜變量(應(yīng)力場，應(yīng)變場，位移場等)的轉(zhuǎn)換(advection), 有兩種算法即為一階算法(first order)與(secon

5、d order)算法，second order算法適用于所 有問題，為推薦算法，一階算法比較簡單，占資源少，速度快；對于動變量(速度，加速度 等)轉(zhuǎn)換(momentum advection ),也有兩種算法, element center projection method 與 half-index shift method，前者為推薦算法，選擇前者就ok 了，如果想仔細(xì)研究，查查 ABAqus文檔就可以了，里面寫的很清楚。在ABAQUS CAE 中，在step模塊中定義選項(xiàng)卡如下所示：Nwne: Ada-LMeshng and SmoothinoPtwrty; * JmprDVB aspec

6、t ratio： Presetve 心舊 fnesh 3必中SfnootH-ig gorthm: DstAnrihed by tMhnii product |Meshing prlctor: Q Current defofrrwcf pcstunPofUicrt fromprewus 鮑 E 4d4ptiw mesh ncremeixCurvature r#fnment： 1BoLndary Reon SwcthngIrvbal f«4tur*3QTr jrrtttxi fey we v-. 30NMh nntraat angfej 60 AdwcbonFr$t p Setefid

7、orderMornentun jdvedion: * Ekmsnt ctntei prcjettion Hat-rwjejr shrft2. 對adaptive remesh過程強(qiáng)度控制在ABAQUS中如何對adaptive remesh過程的強(qiáng)度進(jìn)行控制呢？ ？這就涉及到一下幾個（1） Frequency控制（頻率控制），主要是對整個 step time中網(wǎng)格remesh的次數(shù)進(jìn)行 控制。Remesh次數(shù) n 可以由 n=Increment number /Frequency 來表達(dá)其意義，當(dāng) frequency 的 值為i時，表示每i個增量步進(jìn)行一次 remesh；一個典型的 ALE過程，

8、在每5-100個增量步就需要一次 remesh,對于拉格朗日問題， 改參數(shù)默認(rèn)值為 10,若變形實(shí)在太大，可適當(dāng)調(diào)高，以增加網(wǎng)格重畫的強(qiáng)度，對于爆炸， 碰撞等變形時間極短的問題求解，則在每一個增量步都需要一次 remesh,這時Frequency的值需要設(shè)置得很小，比如設(shè)為1,當(dāng)然，adaptive remesh過程的強(qiáng)度也很高，也會很廢時。對于其他變形不是很劇烈的問題求解，該參數(shù)值可以適當(dāng)調(diào)高。對于歐拉問題，默認(rèn)值為1。（2） Mesh Sweeps參數(shù)控制，當(dāng)該參數(shù)的值為 n時，每一個remesh過程將對網(wǎng)格進(jìn)行n次sweep,其實(shí)這個參數(shù)可以理解為對整個adaptive remesh過程

9、的每一個子過程（remesh過程）的強(qiáng)度進(jìn)行控制。那么，我們先來理解一下sweep的概念，每sweep 一次，abaqus將利用我們設(shè)置好的算法（體積算法，拉普拉斯算法或等位算法）生成一套新的網(wǎng)格，但這個網(wǎng)格不一定是符合 要求的，因此，需要在生成的新網(wǎng)格的基礎(chǔ)上用同樣的方式再進(jìn)行sweep,就像我們求解方程時迭代的概念是一樣的。就這樣一直sweep下去直到sweep的次數(shù)達(dá)到 mesh sweeps參數(shù)的值，這樣就完成了一個 remesh過程中的新網(wǎng)格的生成。 同樣，mesh sweeps參數(shù)的值越高， adaptive remesh過程強(qiáng)度越高，網(wǎng)格優(yōu)化的狀況良好的機(jī)率也就越大。在ABAQU

10、S CAE 中，在step模塊中定義選項(xiàng)卡如下所示：T1 E dit A1.E Adaptive Mesh Domain接下來看第二個部分：對 adaptive remesh其他方面控制主要包括以下幾個參數(shù)設(shè)置：(1) smoothing algorithm 參數(shù)設(shè)置：選項(xiàng): Use enhanced algorithm based on evolving element geometry主要是在幾何學(xué)的方面對我們定義的網(wǎng)格sweep算法(前面提到的三種算法)進(jìn)行增強(qiáng)，目的是為了保證adaptive remesh過程的健壯性，為推薦選項(xiàng)，選它就行了選項(xiàng)：conventional smoothi

11、ng利用我們定義好的算法進(jìn)行計(jì)算，無幾何增強(qiáng)。n* E dll ALE Adaptive Meth LonlroKMame: Ada-1Meshrwj and SmoorfwigPrtofity： ® Jnnprot，村點(diǎn)Prewvc lolhM wsh gr*dhng5modtfwig akgorthni; Wsniwied by dnlrlysisi produdMershrig predKtcr: b Current defevmed potionPositkri from previouK ME adaptive mesh ncrementCurvafcje refwerwt

12、 1(2)Bcwr4r SmwthTfliwsd f«lbjr« tfNg|»= |aaMiKbciMUtrMnt anqte; 60|AdwciionFirst order * Second orderHw5umOemert 5:就 prgrtEM Fde# MOK r:priority參數(shù)設(shè)置：也就是指網(wǎng)格梯度控制（是否保持初始網(wǎng)格梯度，若需要保持初始網(wǎng)格梯度，則對 網(wǎng)格的質(zhì)量將會有影響）。對于拉格朗日問題選擇Improve aspect ratio ,在計(jì)算過程中將考慮到網(wǎng)格單元高寬比的改善，不考慮對初始網(wǎng)格梯度的保持。對于歐拉問題選擇：Preserve i

13、nitial mesh grading ,在計(jì)算過程中保證初始的網(wǎng)格梯 度，但不會考慮到網(wǎng)格寬高比的改善。(3)Meshing predictor 參數(shù)設(shè)置也就是網(wǎng)格節(jié)點(diǎn)位置控制（理想的網(wǎng)格節(jié)點(diǎn)位置控制，將會減少需要的網(wǎng)格 次數(shù)，減少資源浪費(fèi)）sweeps對于拉格朗 日問題選擇Current deformed position ,對于歐拉問題選擇 Position from previous adaptive mesh increment這里就不細(xì)說了，有興趣的話可以看看文檔(4) Curvature refinement 參數(shù)設(shè)置也就是曲率較大的曲線曲面邊界的網(wǎng)格密度控制，默認(rèn)為 的網(wǎng)格密度

14、也就會越大，比較簡單。1,該值越大，則圓角區(qū)(5)Initial remeshing sweeps 參數(shù)設(shè)置也就是ALE過程開始之前對網(wǎng)格的一個優(yōu)化，概念與mesh sweeps類似，因?yàn)槲覀冇锌赡芾靡呀?jīng)變形的很厲害的網(wǎng)格進(jìn)行分析，這時，在分析開始之前，就需要對 網(wǎng)格進(jìn)行重畫。在 ABAQUS CAE中的ale mesh control選項(xiàng)卡剩下的參數(shù)就三個了，即 boundary region smoothing下面的三個參數(shù)：(1) initial feature angle即初始檢測角度的設(shè)置，當(dāng)兩個相鄰的面的法向量大于該角度值的時候，這兩個相 鄰面形成的corner將被檢測出來，在

15、sweep時，網(wǎng)格不允許通過這個corner小于的話就說明，該 corner足夠圓滑，網(wǎng)格可以通過，當(dāng)然，該 corner應(yīng)該是具有 活性的，對 corner活性的控制由下面一個參數(shù)(Transition feature angle)控制，否則也不會被考慮。(2) Transition feature angle控制被檢測出的corner的活性的，如果被檢測處的corner的兩面法線夾角大于該值則該corner在ale過程中是會被考慮的，否則就不會考慮。(3) Mesh constraint angle控制分析過程的一個角度參數(shù)，一般大于45度，設(shè)為默認(rèn)值就可以，在分析過程中，當(dāng)網(wǎng)格內(nèi)某一個角度大于該參數(shù)值時，分析終止，文檔有詳細(xì)介紹。第三部分 ALE adaptive mesh constraint網(wǎng)格約束方面的控制該部分比較復(fù)雜，內(nèi)容很多，我理解得也不清楚，望有高手指教。5個增量步進(jìn)行一次 remesh,ALE具體過程可由下圖表示：mesh