無限單元的案例
Abaqus無限單元的建立方法
什么是無限單元
無限單元是有限元的一種單元，它的單元形式如下圖所示。
2. 無限單元的作用
使用無限單元作為反射邊界，應(yīng)力波將無反射，防止在邊界上產(chǎn)生的應(yīng)力波反射重新進(jìn)入模型，從而確保結(jié)果的正確性。
3. 怎么建立無限單元
（1）先建立有限元模型，然后將需要設(shè)置為無限單元的部分分割，在mesh模塊中設(shè)置該部分單元類型，這樣在inp文件中需要設(shè)置為無限單元的部分就會(huì)集中在一起，方便修改其節(jié)點(diǎn)的編號(hào)。
（2）在inp文件里對(duì)無限單元進(jìn)行單元屬性及編號(hào)的改變，需要注意無限單元的方向，在二維無限單元中，前兩個(gè)節(jié)點(diǎn)所組成直線中點(diǎn)A與后兩個(gè)節(jié)點(diǎn)所組成直線中點(diǎn)B，無限單元的方向就是A指向B的方向，如圖1，二維無限元的方向是朝下；在三維無限單元中，前四個(gè)節(jié)點(diǎn)所組成平面的中心點(diǎn)C與后四個(gè)節(jié)點(diǎn)所組成平面的中心點(diǎn)D，無限單元的方向就是C指向D的方向，如圖1，三維無限單元的方向朝右。
（3）將修改后的inp導(dǎo)入，建立job提交就可以了。
實(shí)例講解-鋼丸撞擊金屬板
機(jī)械噴丸的模擬，其中設(shè)置金屬板邊界部分為無限單元。
（1）建立有限元模型，將金屬板的邊界partition切割，全部設(shè)置為C3D8R單元類型，修改邊界部分為C3D8I（邊界部分的單元類型和其他部分不一樣，方便下步統(tǒng)一更改為無限單元），如圖，建立job-write input 生成inp文件。
（2）修改inp文件的單元屬性及編號(hào)
inp文件以記事本的方式打開，將Element, type（單元類型）C3D8I修改為CIN3D8（無限單元），
下一步修改無限單元節(jié)點(diǎn)編號(hào)的排列順序，在inp的data lines中第一列是單元編號(hào)，后面的是節(jié)點(diǎn)編號(hào)，如圖7。在本文中，無限單元的方向是從金屬板的中心向外。
展開 ABAQUS中無限單元的使用（一）——2D操作案例
ABAQUS中無限單元的使用（一）
——2D操作案例
作者：Nick_Liu
使用無限單元可以用消除應(yīng)力波的反射，從而以較小的模型得到和整體模型一致的計(jì)算結(jié)果，在節(jié)省計(jì)算資源和提高計(jì)算精度方面具有很大優(yōu)勢。一般局部加載對(duì)整體影響不大或者問題為無邊界時(shí)可以使用無限單元模型。
ABAQUS 6.14提供了17種無限單元。非聲學(xué)問題中常用的無限單元有：CINAX4、CIN3D8、CINPE4和CINPS4。使用時(shí)只需在input文件中修改單元類型的關(guān)鍵字即可。
使用無限單元的關(guān)鍵在于模型區(qū)域的劃分。本文通過隨時(shí)空變化壓力載荷作用于圓柱體的案例，演示二維無限單元模型的建立過程。
1）建模過程
Step-1：建立軸對(duì)稱模型，并將模型分成3個(gè)區(qū)域，如圖1所示：
圖1 建模和分塊
Step-2：設(shè)置邊界區(qū)域厚度方向的種子密度為1，寬度方向種子數(shù)量相同；內(nèi)部區(qū)域單元類型為CAX4R，邊界區(qū)域單元類型為CAX4；劃分網(wǎng)格如圖2所示：
圖2 模型網(wǎng)格
Step-3：孤立網(wǎng)格，調(diào)整堆疊方向，如圖3所示：
圖3 調(diào)整堆疊方向
Step-4：賦予材料屬性，建立裝配體和動(dòng)態(tài)分析步（Dynamic, Explicit）。
Step-5：在模型頂部通過VDLoad子程序施加大小隨時(shí)間和空間變化的壓力；同時(shí)對(duì)模型施加對(duì)稱和固定約束，如圖4所示：
圖4 載荷和約束
Step-5：新建Job，輸出Input文件，將Input文件中的單元類型CAX4改為CINAX4。
Step-6：通過修改后的Input文件新建Job并提交計(jì)算。
2）計(jì)算結(jié)果
計(jì)算結(jié)果如圖5所示：
圖5 計(jì)算結(jié)果
可以看到，應(yīng)力傳遞至白色的無限單元邊界時(shí)沒有進(jìn)行反射，其結(jié)果與實(shí)際情況是一致的。
3）注意事項(xiàng)
a. 厚度方向只能有一個(gè)單元；
b.
展開 ABAQUS案例—ABAQUS中聲固耦合、聲輻射分析方法 ￥4
噪聲輻射分析中，需要模擬附著在結(jié)構(gòu)上的外部空氣，而且它是向外無限延伸的，因此直接用聲學(xué)有限單元去模擬無限的空氣區(qū)域是不合理的。在Abaqus中可以通過兩種方式來模擬無限聲學(xué)介質(zhì)的影響：一,使用聲學(xué)無限單元；二，用阻抗邊界來模擬。
在對(duì)外部的噪聲輻射問題進(jìn)行仿真分析時(shí)，無限單元法的應(yīng)用已經(jīng)越來越廣泛。無限單元可以直接在結(jié)構(gòu)上定義，或者也可以在聲學(xué)有限單元區(qū)域的終面上定義。
對(duì)于邊界阻抗技術(shù)，實(shí)質(zhì)上屬于無反射邊界條件。然而當(dāng)用此來模擬結(jié)構(gòu)外部的區(qū)域時(shí)，結(jié)構(gòu)與輻射表面的距離必須足夠大（通常取聲波波長的1/3）。
聲學(xué)無限單元計(jì)算公式與聲輻射阻抗邊界的計(jì)算有幾個(gè)關(guān)鍵的區(qū)別：無限單元采用更高階的差值函數(shù)，而聲輻射邊界則采用一階差值函數(shù)。雖然無限元計(jì)算每個(gè)單元的花費(fèi)更高，但是無限單元的要比阻抗邊界精確很多，因此通過減小無限元的單元規(guī)模，從而可以大大的降低結(jié)構(gòu)總的計(jì)算時(shí)間；本案例即是講解無限元單元法在模擬噪聲分析中的應(yīng)用。
展開 基于ABAQUS單點(diǎn)顯式VDLOAD/隱式DLOAD激光沖擊加載(圓形光斑和方形光斑) ￥50
材料采用Ti-6Al-4V鈦合金，有限單元區(qū)賦予彈塑性參數(shù)，無限單元賦予彈性參數(shù)。具體參數(shù)如下：
密度：4.5e-9；彈性：1.2e5,0.34；塑性：A:1098 B:1092 C:0.014 n:0.93 參考應(yīng)變率：1
裝配：全局坐標(biāo)原點(diǎn)與有限單元頂點(diǎn)重合。
分析步創(chuàng)建：隱式/顯式分析步，增量采用固定增量步
接口設(shè)定：
網(wǎng)格劃分，無限單元部分定義網(wǎng)格節(jié)點(diǎn)方向排布
建立模型輸出inp文件，將無限單元部分的單元類型改為CIN3D8
如果沖擊光斑為圓形光斑，網(wǎng)格細(xì)化至50微米，如果為方向光斑，網(wǎng)格100微米
圓形光斑在空間表現(xiàn)為高斯分布，表現(xiàn)為從中心區(qū)域到邊緣沖擊載荷逐漸變小。
方向光斑在空間表現(xiàn)為均布載荷，其峰值載荷為圓形光斑的0.618倍，一些研究表明相同激光參數(shù)下方形光斑搭接沖擊材料疲勞性能較高。
展開 
無限單元要保證單元堆棧方向都朝外側(cè)。需要自己重新排列節(jié)點(diǎn)順序
【公益帖】abaqus中的ETOTAL 能量澄清，解釋為何有時(shí)為負(fù)？
如下
ETOTAL 總的能量平衡
ETOTAL = ALLKE + ALLVD + ALLSD + ALLKL + ALLFD + ALLJD + ALLIE – ALLWK
ALLKE 動(dòng)能
ALLVD 耗散能（粘性引起，不包括 ALLSD 和 ALLCD）
ALLSD 耗散能（自動(dòng)穩(wěn)定引起，如接觸）
ALLKL 動(dòng)能損失（沖擊引起）
ALLFD 耗散能（摩擦引起）
ALLJD 耗散能（電流引起）
ALLIE 總應(yīng)變能
ALLIE = ALLSE + ALLPD + ALLCD + ALLAE + ALLQB + ALLEE + ALLDMD
ALLSE 應(yīng)變能(可恢復(fù))
ALLPD 耗散能（塑性變形引起）
ALLCD 耗散能（粘彈性、蠕變、膨脹引起）
ALLAE 偽應(yīng)變能
ALLQB 耗散能（無限單元引起，如無反射邊界）
ALLEE 靜電能
ALLDMD 耗散能（裂紋引起）
ALLWK 外力功
那么出現(xiàn)了ETOTAL能量為負(fù)，并不是代表計(jì)算就有問題，總能來說ETOTAL月接近零越好。但是有時(shí)可能是一個(gè)比較大的值，這是否就說明計(jì)算有問題嗎？ 也不一定要看ETOTAL與 總動(dòng)能（或者總勢能）的比值是否可以比擬，如果在1%以內(nèi)我們認(rèn)為也是可以接受的。
展開 有限元經(jīng)典之二《有限單元法基本原理與應(yīng)用》（第二版）
2-3 單元應(yīng)變
2-4 初應(yīng)變
2-5 單元應(yīng)力
2-6 等效結(jié)點(diǎn)力力與單元剛度矩陣
2-7 結(jié)點(diǎn)荷載
2-8 結(jié)點(diǎn)平衡方程與整體剛度矩陣
2-9 用編碼法建立整體剛度矩陣
2-10 計(jì)算實(shí)例
參考文獻(xiàn)
第3章 單元分析
3-1 虛位移原理
3-2 單元位移
3-3 單元應(yīng)變與應(yīng)力
3-4 結(jié)點(diǎn)力與單元剛度矩陣
3-5 結(jié)點(diǎn)荷載
3-6 虛位移原理應(yīng)用實(shí)例——梁單元
3-7 應(yīng)變能和余應(yīng)變能
3-8 最小勢能原理
3-9 最小余能原理
3-10 雜交單元
3-11 雜交單元實(shí)例——平面矩形單元
3-12 混合能量原理
3-13 復(fù)合單元
參考文獻(xiàn)
第4章 整體分析
第5章 平面問題高次單元
第6章 彈性力學(xué)軸對(duì)稱問題
第7章 彈性力學(xué)空間問題
第8章 形函數(shù)、坐標(biāo)變換、等參數(shù)單元與無限單元
第9章 各種平面與空間單元的比較、應(yīng)用實(shí)例
第10章 彈性薄板
第11章 彈性薄殼
第12章 軸對(duì)稱殼
第13章 彈性厚板和厚殼
第14章 流體力學(xué)問題
第15章 熱傳導(dǎo)問題
第16章 非線性有限元分析方法
第17章 塑性力學(xué)問題
第18章 混凝土徐變、一般粘彈性及粘塑性問題
第19章 彈性穩(wěn)定問題
第20章 大位移問題
第21章 斷裂力學(xué)問題
第22章 結(jié)構(gòu)動(dòng)力學(xué)問題
第23章 巖石力學(xué)問題
第24章 土力學(xué)問題
第25章 混凝土與鋼筋混凝土結(jié)構(gòu)
第26章 工程反分析
第27章 網(wǎng)格自動(dòng)生成、誤差估計(jì)與自適應(yīng)技術(shù)
附錄
第一部分
有限單元法原理與應(yīng)用 第2版(朱伯芳)[1].part01.rar
展開 Actran助力汽車全頻率段聲學(xué)響應(yīng)預(yù)測與優(yōu)化
使用非連續(xù)伽遼金法求解線性化的歐拉方程，支持單元階次的自適應(yīng)，能夠解決大規(guī)模工程中的聲學(xué)仿真問題
渦輪機(jī)械噪聲分析模塊—Actran TM. 采用并行求解技術(shù)，能夠模擬大馬赫數(shù)、大幾何尺寸及形狀復(fù)雜的大模型，如航空發(fā)動(dòng)機(jī)渦輪機(jī)械聲輻射問題
2. 軟件優(yōu)勢
基于通用的有限單元和無限單元的建模環(huán)境，具有完備的結(jié)構(gòu)單元庫、材料庫、邊界條件和求解方法
支持頻率自適應(yīng)網(wǎng)格自動(dòng)生成
支持并行計(jì)算提高仿真效率
支持與上游結(jié)構(gòu)有限元軟件，如Nastran、Abaqus、Ansys等進(jìn)行鏈?zhǔn)椒抡妫行蚣啥鄠€(gè)學(xué)科，提高整體仿真的準(zhǔn)確性
支持求解高速、強(qiáng)對(duì)流區(qū)域的聲傳播、超大規(guī)模聲學(xué)域的仿真問題
采用能量分析方法作為有限元分析的補(bǔ)充，擴(kuò)展振動(dòng)聲學(xué)分析的頻率范圍
結(jié)果后處理功能豐富，包括瀑布圖、階次分析、頻響函數(shù)查看、聲音文件恢復(fù)
通過腳本及API編程實(shí)現(xiàn)優(yōu)化
典型應(yīng)用
??新能源電驅(qū)動(dòng)系統(tǒng)噪聲仿真模擬
針對(duì)新能源汽車中變頻電機(jī)的電流諧波和電磁諧波引起的高頻電磁噪聲，利用Actran對(duì)電機(jī)輻射噪聲進(jìn)行仿真、評(píng)估及優(yōu)化，為電驅(qū)動(dòng)總成聲學(xué)性能的提高提供依據(jù)。
展開 無限元在Abaqus靜力分析中的應(yīng)用
因此，一種定義無限元的常用方法為首先確定極點(diǎn)的位置，然后根據(jù)極點(diǎn)和第一個(gè)單元面（或單元邊）的節(jié)點(diǎn)確定延伸方向，最后根據(jù)延伸方向和距離確定非第一個(gè)單元面的剩余節(jié)點(diǎn)，這樣的方法可以保證所有的無限元均圍繞同一個(gè)極點(diǎn)。
圖7 可接受和不可接受的二維無限元
2.3
在平面應(yīng)力/應(yīng)變分析中使用二維無限元
下面本文簡要介紹如何在平面應(yīng)力/應(yīng)變分析中建立包含二維無限元的有限元模型，在三維應(yīng)力分析中的建立方法可參考該方法進(jìn)行。以4節(jié)點(diǎn)平面應(yīng)力無限元CINPS4為例，它與常規(guī)平面應(yīng)力單元CPS4具有相同的節(jié)點(diǎn)數(shù)量，唯一的不同體現(xiàn)在節(jié)點(diǎn)順序的定義上，因?yàn)锳baqus要求只能是無限元的第一條單元邊與常規(guī)單元相連。因此，我們可以在劃分網(wǎng)格時(shí)將無限元按常規(guī)單元處理，而后在inp文件中通過修改單元類型和相應(yīng)的節(jié)點(diǎn)定義來建立無限元。需要注意的是，在Abaqus中單元節(jié)點(diǎn)都是以逆時(shí)針順序編號(hào)的，例如圖8給出了一個(gè)常規(guī)4節(jié)點(diǎn)平面單元的節(jié)點(diǎn)定義。
圖8 4節(jié)點(diǎn)平面單元節(jié)點(diǎn)定義
假定該單元的單元號(hào)為1，則該單元在inp文件中一種可能的節(jié)點(diǎn)定義順序?yàn)椋?1, 6, 8, 1, 3
如果將圖8中的單元看作是一個(gè)無限元，并且節(jié)點(diǎn)1和節(jié)點(diǎn)3與常規(guī)單元連接，則上面的節(jié)點(diǎn)定義順序顯然有誤，因?yàn)榇藭r(shí)節(jié)點(diǎn)6和節(jié)點(diǎn)8構(gòu)成的單元邊為第一條單元邊，正確的節(jié)點(diǎn)定義順序應(yīng)為：
1, 1, 3, 6, 8
顯然定義無限元的關(guān)鍵在于確定無限元與常規(guī)單元相連接的節(jié)點(diǎn)號(hào)（下面稱為近場邊界節(jié)點(diǎn)），然后根據(jù)這些節(jié)點(diǎn)確定節(jié)點(diǎn)定義信息中的起始節(jié)點(diǎn)（即節(jié)點(diǎn)1），由于節(jié)點(diǎn)總是按逆時(shí)針排序的，因此只要起始節(jié)點(diǎn)確定，則其余節(jié)點(diǎn)依次排列即可。
一旦確定了常規(guī)單元和無限元所在區(qū)域，則近場邊界節(jié)點(diǎn)的編號(hào)是已知的，問題轉(zhuǎn)換為如何判斷起始節(jié)點(diǎn)，即判斷圖8中節(jié)點(diǎn)1和節(jié)點(diǎn)3哪一個(gè)為起始節(jié)點(diǎn)。
展開 ANSYS有限元網(wǎng)格介紹
結(jié)構(gòu)分析中有結(jié)構(gòu)點(diǎn)單元（如MASS21）、結(jié)構(gòu)線單元（如LINK10、LINK180）、結(jié)構(gòu)梁單元（如BEAM23、BEAM188）、結(jié)構(gòu)實(shí)體單元（如PLANE25、SOLID45、SOLID185、SOLID186）、結(jié)構(gòu)殼單元（如SHELL51、SHELL181）、結(jié)構(gòu)管單元（如PIPE59）、接觸單元（如CONTAC174、TARGE170）、顯示動(dòng)力單元（如LINK160、SOLID164、BEAM161）等。
熱力學(xué)分析中有熱點(diǎn)單元（如MASS71）、熱線單元（如LINK31、LINK32）、熱實(shí)體單元（如PLANE35、SOLID70、SOLID87）、熱殼單元（如SHELL57、SHELL131、SHELL132）、熱電單元（如PLANE67、SOLID69）等。
流體結(jié)構(gòu)分析中還包括一些特有的單元，如FLUID29、FLUID30、FLUID116、FLUID129、FLUID136、FLUID141等。
電子分析中包括電磁單元（如PLANE53、PLANE230、SOLID96、SOLID123、INTER115、HF118、HF121）、電路單元（如SOURC36、CIRCU94、CIRCU124）、機(jī)電單元（如TRANS109、TRANS126）等。
除此以外ANSYS中還包括一些其他單元，如耦合場單元PLANE13、矩陣單元MATRIX27、無限單元INFIN9、網(wǎng)分單元MESH200、表面單元SURF151等等。
展開 基于ABAQUS蠕變儲(chǔ)層稠油蒸汽吞吐開發(fā)過程數(shù)值模擬
圖1 地層幾何模型
網(wǎng)格劃分
整個(gè)地層劃分為11個(gè)不同的層位，其中具有孔隙壓力的軸對(duì)稱縮減積分單元CAX8RP用于模擬井附近的巖層，當(dāng)使用二階單元時(shí)，一般采用縮減積分，因?yàn)樗ǔＬ峁└鼫?zhǔn)確的結(jié)果并且比完全積分具有更小的計(jì)算成本；遠(yuǎn)場區(qū)域則使用軸對(duì)稱無限單元CINAX5R建模，以提供橫向剛度。網(wǎng)格劃分后的有限元模型如圖2所示。
圖2 地層有限元模型
模擬參數(shù)
地表土層與泥巖層
對(duì)地表土層S1、T1與深部泥巖層U1和L1使用Drucker-Prager塑性模型建模，其彈性和非彈性材料屬性均列于表1中。本例使用沒有中間主應(yīng)力效應(yīng)的Drucker-Prager模型的線性形式，因此流應(yīng)比，即三軸拉伸強(qiáng)度與三軸壓縮強(qiáng)度之比k=1；該模型假設(shè)為流動(dòng)法則非關(guān)聯(lián)，即在式（1）塑性本鉤矩陣中屈服面、加載面（后續(xù)加載面）與塑性勢面不同，即材料剛度矩陣不對(duì)稱，因此，使用非對(duì)稱矩陣存儲(chǔ)和非對(duì)稱求解器解可以顯著改善該非線性解的收斂性；硬化/軟化行為是Drucker-Prager塑性材料定義的擴(kuò)展，其數(shù)據(jù)列于表1中。因?yàn)檫@些層位遠(yuǎn)離頻繁加載，所以可不考慮蠕變。
式（1）
表1 地表土層與泥巖層Drucker-Prager模型參數(shù)
蠕變特性顯著的中間層
D1到D7使用改進(jìn)的Drucker-Prager Cap塑性模型模擬，材料屬性數(shù)據(jù)列于表2中。
展開 
ABAQUS---輪軌瞬態(tài)滾動(dòng)接觸有限元模型(直線半輪對(duì)) ￥888
</p><p><span style="color: rgb(25, 27, 31);"> </span>③流體單元、無限單元、彈簧單元、質(zhì)量單元、阻尼單元以及旋轉(zhuǎn)慣性單元不能進(jìn)行數(shù)據(jù)傳遞。</p><p><span style="color: rgb(25, 27, 31);"> </span>④已經(jīng)定義了初始狀態(tài)場的部件實(shí)體、網(wǎng)格及其相應(yīng)的部件，都不能再在CAE中進(jìn)行修改，不允許為其賦予新的截面屬性、材料方向、定義質(zhì)量和慣性矩，在后續(xù)分析中會(huì)自動(dòng)替換為原始分析中的材料。</p><p><span style="color: rgb(25, 27, 31);"> </span>⑤在standard和explicit之間進(jìn)行數(shù)據(jù)傳遞時(shí)，要保證沙漏的計(jì)算方法一致。</p><h3>（1）Assembly</h3><p><span style="color: rgb(25, 27, 31);"> </span>復(fù)制隱式CAE模型，<strong>模型中部件名字千萬不能更改替換</strong>。</p><p><span style="color: rgb(25, 27, 31);"> </span>注意：取掉隱式模型中在edit keywords中加入的*contact controls, stabilize，該關(guān)鍵詞不適用于顯式，否則報(bào)錯(cuò)。
展開 設(shè)計(jì)仿真 | 2023年海克斯康工業(yè)軟件中國區(qū)高校合作計(jì)劃招募公告
Actran提供了一個(gè)豐富的材料模型庫、一個(gè)包含無限單元的完整單元庫、解決大型問題的高性能解決方案和一個(gè)用戶友好的GUI。Actran具有模塊化的特性，可以根據(jù)您的應(yīng)用程序和需求進(jìn)行部署。
11
多尺度材料平臺(tái)—Digimat
Digimat提供了設(shè)計(jì)工具，通過對(duì)復(fù)合材料的局部性能準(zhǔn)確描述，使用戶對(duì)復(fù)合材料產(chǎn)品有100%的信心。精確的材料建模可以減少“安全系數(shù)”——允許復(fù)合材料100%地發(fā)揮它們的潛力，最大限度地提高它們對(duì)金屬的競爭力，并進(jìn)行產(chǎn)品的輕量化設(shè)計(jì)。
在 COMSOL 中簡化地下水流建模
半徑為 20m、高度為 3m 的水庫中，半徑為 0.5m 的井的幾何模型，其周圍是一個(gè)無限的單元域。
使用無限元是為了使我們可以在離井很遠(yuǎn)的地方施加壓力而不增加建模域。這里顯示的幾何圖形將井解析為一個(gè)圓柱形的表面。為了能夠應(yīng)用邊界條件，必須將井的圓柱體從儲(chǔ)層中切割出來。另外，也可以使用質(zhì)量通量的邊條件，但前提是我們要應(yīng)用質(zhì)量通量而不是壓力。我們可以使用井邊界條件，它適用于壓力和質(zhì)量通量條件。
我們用完全相同的網(wǎng)格設(shè)置來比較這兩種情況下的網(wǎng)格。在這個(gè)案例中，我們劃分了 65,674 個(gè)域單元，而使用井邊界條件，僅劃分了 28,728 個(gè)域單元。這還不到網(wǎng)格單元數(shù)量的一半。
使用相同設(shè)置的在完全解析井時(shí)和使用井邊界條件時(shí)的網(wǎng)格比較。
這個(gè)優(yōu)勢只有在我們得到一個(gè)準(zhǔn)確的解時(shí)才有用。繼續(xù)使用這個(gè)測試案例，我們在井口施加一個(gè) 1 kg/s 的質(zhì)量流速，M0。這相當(dāng)于在面積為 A 的邊界處的質(zhì)量通量為
。在長度為 l 的邊處的質(zhì)量通量為
。壓力在外部無限元的邊界是固定的。
一維繪圖顯示沿中心線的壓力與井外的方法幾乎完全一致。
沿截線的壓力比較。
與指定邊的質(zhì)量通量相反，井功能考慮了井半徑，即使沒有明確解析也考慮了。質(zhì)量通量 邊特征不合適計(jì)算井中的壓力，因?yàn)樗豢紤]徑向的膨脹。井邊界條件提供了一個(gè)變量 dl.well1.p，該變量給出了井壓。
傳熱與多孔介質(zhì)流的耦合應(yīng)用：地?zé)峄毓? 井邊界條件可用于前面提到的地?zé)峄毓嗍纠Ｎ覀冊谥暗奈恼轮薪榻B的模型的略微修改版本。在這種情況下，地?zé)岬叵滤?150 l/s 的速度通過生產(chǎn)井生產(chǎn)。在用于產(chǎn)生熱量后，水再以相同的速率重新注入，此時(shí)水溫為 5°C。在外部邊界施加 2 mm/m 的水平水力梯度。
地?zé)峄毓嗄Ｐ偷哪Ｐ驮O(shè)置(左)和網(wǎng)格(右)。
展開 abaqus檢驗(yàn)總結(jié)1-論壇整理
網(wǎng)格質(zhì)量檢查在Mesh 功能模塊中，點(diǎn)擊主菜單Mesh→Verify，可以選擇部件、實(shí)體、幾何區(qū)域或單元，檢查其網(wǎng)格的質(zhì)量，獲得節(jié)點(diǎn)和單元信息。在Verify Mesh 對(duì)話框，選擇Statistical Checks（統(tǒng)計(jì)檢查）可以檢查單元的幾何形狀，選擇Analysis Checks（分析檢查）可以檢查分析過程中會(huì)導(dǎo)致錯(cuò)誤或警告信息的單元。單擊Highlight 按鈕，符合檢查判據(jù)的單元就會(huì)以高亮度顯示出來。
單元選擇
ABAQUS 擁有433 種單元，分8 大類：連續(xù)體單元（continuum element，即實(shí)體單元solid element）、殼單元、薄膜單元、梁單元、桿單元、剛體單元、連接單元和無限單元
線性單元（即一階單元）；二次單元（即二階單元）；修正的二次單元（只有Tri 或Tet 才有此類型）。 ABAQUS/Explicit 中沒有二次完全積分的連續(xù)體單元。
線性完全積分單元的缺點(diǎn)：承受彎曲載荷時(shí)，會(huì)出現(xiàn)剪切自鎖，造成單元過于剛硬，即使劃分很細(xì)的網(wǎng)格，計(jì)算精度仍然很差。
二次完全積分單元的優(yōu)點(diǎn)：
應(yīng)力計(jì)算結(jié)果很精確，適合模擬應(yīng)力集中問題；
一般情況下，沒有剪切自鎖問題。
二次完全積分單元的缺點(diǎn)：
不能用于接觸分析；
對(duì)于彈塑性分析，如果材料不可壓縮（例如金屬材料），則容易產(chǎn)生體積自鎖；
當(dāng)單元發(fā)生扭曲或彎曲應(yīng)力有梯度時(shí)，有可能出現(xiàn)某種程度的自鎖。
線性減縮積分單元在單元中心只有一個(gè)積分點(diǎn)，存在沙漏數(shù)值問題而過于柔軟。采用這種單元模擬承受彎曲載荷的結(jié)構(gòu)時(shí)，沿厚度方向上至少應(yīng)劃分四個(gè)單元。
展開 