第一篇：ABAQUS單元小結(jié)

ABAQUS單元小結(jié)

1、單元表征

單元族：單元名字里開始的字母標(biāo)志著這種單元屬于哪一個單元族。C3D8I是實體單元； S4R是殼單元； CINPE4是無限元； 梁單元； 剛體單元； 膜單元；

特殊目的單元，例如彈簧，粘壺和質(zhì)量； 桁架單元。

自由度dof（和單元族直接相關(guān)）：每一節(jié)點處的平動和轉(zhuǎn)動 1 1方向的平動 2 2方向的平動 3 3方向的平動 4 繞1軸的轉(zhuǎn)動 5 繞2軸的轉(zhuǎn)動 6 繞3軸的轉(zhuǎn)動 開口截面梁單元的翹曲 8 聲壓或孔隙壓力 9 電勢 度（或物質(zhì)擴散分析中歸一化濃度）12＋梁和殼厚度上其它點的溫度 軸對稱單元 r方向的平動 2 z方向的平動 6 r-z方向的轉(zhuǎn)動

節(jié)點數(shù)：決定單元插值的階數(shù)

數(shù)學(xué)描述：定義單元行為的數(shù)學(xué)理論

積分：應(yīng)用數(shù)值方法在每一單元的體積上對不同的變量進(jìn)行積分。大部分單元采用高斯積分方法計算單元內(nèi)每一高斯點處的材料響應(yīng)。單元末尾用字母“R”識別減縮積分單元，否則是全積分單元。

ABAQUS擁有廣泛適用于結(jié)構(gòu)應(yīng)用的龐大單元庫。單元類型的選擇對模擬計算的精度和效率有重大的影響；

節(jié)點的有效自由度依賴于此節(jié)點所在的單元類型；

單元的名字完整地標(biāo)明了單元族、單元的數(shù)學(xué)描述、節(jié)點數(shù)及積分類型； 所用的單元都必須指定單元性質(zhì)選項。單元性質(zhì)選項不僅用來提供定義單元幾何形狀的附加數(shù)據(jù)，而且用來識別相關(guān)的材料性質(zhì)定義；

對于實體單元，ABAQUS參考整體笛卡爾坐標(biāo)系來定義單元的輸出變量，如應(yīng)力和應(yīng)變?？梢杂?ORIENTATION選項將整體坐標(biāo)系改為局部坐標(biāo)系；

對于三維殼單元，ABAQUS參考建立在殼表面上的一個坐標(biāo)系來定義單元的輸出變量。可以用*ORIENTATION選項更改這個參考坐標(biāo)系。2．實體單元(C)實體單元可在其任何表面與其他單元連接起來。C3D:三維單元

CAX:無扭曲軸對稱單元，模擬3600的環(huán)，用于分析受軸對稱載荷作用，具有軸對稱幾何形狀的結(jié)構(gòu)；

CPE:平面應(yīng)變單元，假定離面應(yīng)變ε33為零，用力模擬厚結(jié)構(gòu)； CPS:平面應(yīng)力單元，假定離面應(yīng)力σ33為零，用力模擬薄結(jié)構(gòu)； 廣義平面應(yīng)變單元包括附加的推廣：離面應(yīng)變可以隨著模型平面內(nèi)的位置線性變化。這種數(shù)學(xué)描述特別適合于厚截面的熱應(yīng)力分析。

可以扭曲的軸對稱單元：用來模擬初始時為軸對稱的幾何形狀，且能沿對稱軸發(fā)生扭曲。這些單元對于模擬圓柱形結(jié)構(gòu)，例如軸對稱橡膠套管的扭轉(zhuǎn)很有用。

反對稱單元的軸對稱單元：用來模擬初始為軸對稱幾何形狀的反對稱變形。適合于模擬像承受剪切載荷作用的軸對稱橡膠支座一類的問題。

如果不需要模擬非常大的應(yīng)變或進(jìn)行一個復(fù)雜的，改變接觸條件的問題，則應(yīng)采用二次減縮積分單元（CAX8R,CPE8R,CPS8R,C3D20R）

如果存在應(yīng)力集中，則應(yīng)在局部采用二次完全積分單元（CAX8,CPE8,CPS8,C3D20等）。

對含有非常大的網(wǎng)格扭曲模擬（大應(yīng)變分析），采用細(xì)網(wǎng)格劃分的線性減縮積分單元（CAX4R,CPE4R,CPS4R,C3D8R等）

對接觸問題采用線性減縮積分單元或非協(xié)調(diào)元（CAX4I,CPE4I,CPS4I, C3D8I）的細(xì)網(wǎng)格劃分。

如果在模型中采用非協(xié)調(diào)元應(yīng)使網(wǎng)格扭曲減至最小。三維情況應(yīng)盡可能采用塊狀單元（六面體）。當(dāng)幾何形狀復(fù)雜時，完全采用塊體單元構(gòu)造網(wǎng)格會很困難，因此可能有必要采用稧形和四面體單元，但盡量少用，并遠(yuǎn)離需要精確求解的區(qū)域。

一些前處理程序包括網(wǎng)格劃分方法，它們可用四面體單元構(gòu)造任意形狀的網(wǎng)格。只要采用二次四面體單元（C3D10），其結(jié)果對小位移問題應(yīng)該是合理的。

小結(jié)：

在實體單元中所用的數(shù)學(xué)公式和積分階數(shù)對分析的精度和花費有顯著的影響；

使用完全積分單元，尤其是一階（線性）單元，容易形成自鎖現(xiàn)象，正常情況不用；

一階減縮積分單元容易出現(xiàn)沙漏現(xiàn)象；充分的單元細(xì)化可減小這種問題； 在分析中如有彎曲位移，且采用一階減縮積分單元時，應(yīng)在厚度方向至少用4個單元；

沙漏現(xiàn)象在二階減縮積分單元中較少見，一般問題應(yīng)考慮應(yīng)用這些單元； 非協(xié)調(diào)單元的精度依賴于單元扭曲的量值；

結(jié)果的數(shù)值精度依賴于所用的網(wǎng)格，應(yīng)進(jìn)行網(wǎng)格細(xì)化研究以確保該網(wǎng)格對問題提供了唯一的解答。但是應(yīng)記住使用一個收斂網(wǎng)格不能保證計算結(jié)果與問題的實際行為相匹配：它還依賴于模型其他方面的近似化和理想化程度；

通常只在想要得到精確結(jié)果的區(qū)域細(xì)劃網(wǎng)格； ABAQUS具有一些先進(jìn)特點如子模型，它可以幫助對復(fù)雜模擬得到有用的結(jié)果。

3．殼單元(S)可以模擬有一維尺寸（厚度）遠(yuǎn)小于另外兩維尺寸，且垂直于厚度方向的應(yīng)力可以忽略的結(jié)構(gòu)。

一般殼單元：S4R,S3R,SAX1,SAX2,SAX2T。對于薄殼和厚殼問題的應(yīng)用均有效，且考慮了有限薄膜應(yīng)變；

薄殼單元：STRI3,STRI35,STRI65,S4R5,S8R5,S9R5,SAXA。強化了基爾霍夫條件，即：垂直于殼中截面的平面保持垂直于中截面；

厚殼單元：S8R,S8RT。二階四邊形單元，在小應(yīng)變和載荷使計算結(jié)果沿殼的跨度方向上平緩變化的情況下，比普通單元產(chǎn)生的結(jié)果更精確；

對于給定的應(yīng)用，判斷是屬于薄殼還是厚殼問題，一般：如果單一材料制造的各向同性殼體的厚度和跨度之比在1/20－1/10之間，認(rèn)為是厚殼問題；如果比值小于1/30，則認(rèn)為是薄殼問題；若介于1/30－1/20之間，則不能明確劃分。由于橫向剪切柔度在復(fù)合材料層合殼結(jié)構(gòu)中作用顯著，故比值（厚跨比）將遠(yuǎn)小于“薄”殼理論中采用的比值。具有高柔韌中間層的復(fù)合材料（“三明治”復(fù)合材料）有很低的橫向剪切剛度并且?guī)缀蹩偸潜挥脕砟M“厚”殼；

橫向剪切力和剪切應(yīng)變存在于普通殼單元和厚殼單元中。對于三維單元，提供了可估計的橫向剪切應(yīng)力。計算這些應(yīng)力時忽略了彎曲和扭轉(zhuǎn)變形的耦合作用，并假定材料性質(zhì)和彎曲力矩的空間梯度很??；

殼單元可以使用每個單元的局部材料方向，各項異型材料的數(shù)據(jù)，如纖維增強復(fù)合材料，以及單元輸出變量，如應(yīng)力和應(yīng)變，都按局部材料方向而定義。在大位移分析中，殼單元上的局部材料軸隨著材料各積分點上的平均運動而轉(zhuǎn)動；

線性、有限薄膜應(yīng)變、四邊形殼單元（S4R）是較完備的而且適合于普通范圍的應(yīng)用；

線性、有限薄膜應(yīng)變、三角形殼單元（S3R）可作為通用的殼單元來應(yīng)用。由于在單元內(nèi)部近似為應(yīng)變場，精細(xì)的網(wǎng)格劃分可用于求解彎曲變形和高應(yīng)變梯度；

考慮到在復(fù)合材料層合殼模型中剪切柔度的影響，將采用“厚”殼單元（S4R,S3R,S8R）

四邊形或三角形的二次殼單元，用于一般的小變形薄殼是很有效的。它們對剪力自鎖和薄膜鎖死是不敏感的；

在接觸模擬中不用選用二階三角形殼單元（STRI65），要采用9節(jié)點的四邊形殼單元（S9R5）;對于僅經(jīng)歷幾何線性行為的非常大的模型，線性、薄殼單元（S4R5）一般將比通用殼單元花費更少；

小結(jié)：

殼單元的橫截面特性可以由沿厚度方向的數(shù)值積分確定（*SHELL SECTION），或在分析開始時應(yīng)用計算的橫截面剛度（*SHELL GENERAL SECTION）；

*SHELL GENERAL SECTION是非常有效的，但僅用于線性材料，*SHELL SECTION可用于線性和非線性材料；

數(shù)值積分在沿殼厚度方向的一系列積分點上進(jìn)行。這些積分點就是單元變量可以被輸出的位置。最外層的積分點位于殼單元的表面。

殼單元法線方向決定了單元的正和負(fù)表面，為了正確地定義接觸和解釋輸出數(shù)據(jù)，必須知道其對應(yīng)的是哪個面。殼法線還定義了施加在單元上正壓力載荷的方向，并可以在ABAQUS/Post中畫出；

殼單元利用材料方向局部化到每個單元。在大位移分析中，局部材料軸隨單元而轉(zhuǎn)動。*ORIENTATION被用來定義非默認(rèn)的局部坐標(biāo)系統(tǒng)。單元的變量，如應(yīng)力和應(yīng)變，在局部方向輸出；

*TRANSFORM定義節(jié)點的局部坐標(biāo)系，集中載荷和邊界條件被應(yīng)用在局部坐標(biāo)系中。所用節(jié)點的輸出，如位移，也默認(rèn)為基于局部的坐標(biāo)系；

矢量圖可以使模擬結(jié)果可視化，特別是用來觀察結(jié)構(gòu)的運動和載荷路徑。

4.梁單元（B）

模擬一維尺寸（長度）遠(yuǎn)大于另外二維尺寸的構(gòu)件，且只有長度方向上的應(yīng)力比較顯著。

對于包含接觸的任何模擬，應(yīng)使用一階、剪切變形的梁單元（B21,B31）如果結(jié)構(gòu)剛度非常大或者非常柔軟，在幾何非線性模擬中應(yīng)當(dāng)使用雜交梁單元（B21H,B32H,等）

使用歐拉－伯努利（三次）梁單元（B23,B33）精度很高，可模擬承受分布載荷作用的梁，例如動態(tài)振動分析。如果橫向剪切變形也很重要，要使用鐵摩辛柯（二次型）梁單元（B22,B32）

模擬有開口薄壁橫截面的結(jié)構(gòu)，應(yīng)當(dāng)使用考慮了開口截面翹曲理論的梁單元（B31OS,B32OS）

小結(jié)：

梁單元的性質(zhì)由截面（*BEAM SECTION或*BEAM GENERAL SECTION）的數(shù)值積分決定，或直接給出截面積、慣性矩和扭轉(zhuǎn)常數(shù)（*BEAM GENERAL SECTION）；

當(dāng)使用*BEAM GENERAL SECTION選項時，模擬開始時進(jìn)行一次數(shù)值積分，并且假定材料是彈性的；

ABAQUS包括大量的標(biāo)準(zhǔn)橫截面形狀。其它形狀可以通過定義SECTION=ARBITRARY來模擬；

必須定義橫截面取向，方法是通過給出第三個節(jié)點，或者在單元性質(zhì)定義中定義一個矢量。截面取向在ABAQUS后處理中可以畫出；

當(dāng)梁作為殼的加強構(gòu)件使用時，梁的橫截面可能偏離節(jié)點；

線性和二次型包含剪切變形的影響，三次型梁不考慮剪切柔度。開口截面梁準(zhǔn)確地模擬了扭轉(zhuǎn)和薄壁開口截面翹曲（包括翹曲約束）的影響；

多點約束和約束方程可以用來連接模型中鉸接、剛性連接等節(jié)點的自由度；

“彎矩”型圖使得像梁這樣的一維單元的結(jié)果很清楚地表示出來； ABAQUS后處理圖的硬拷貝可以得到PostScript和HPGL的格式。

5．桁架單元（T）

只能承受拉伸和壓縮載荷的桿，不能承受彎曲，模擬鉸接框架結(jié)構(gòu)，近似模擬線纜和彈簧。

6．剛體單元（R）

沒有獨立的自由度。7．非線性分析

小結(jié)：

結(jié)構(gòu)問題中存在著三種非線性來源：材料、幾何和邊界（接觸）。這些因素的任意組合都可以出現(xiàn)在ABAQUS的分析中；

幾何非線性發(fā)生在位移量值影響結(jié)構(gòu)響應(yīng)的情況下。這包括大位移和轉(zhuǎn)動效應(yīng)、突然翻轉(zhuǎn)和載荷硬化；

非線性問題是利用牛頓－拉弗森方法來進(jìn)行迭代求解的。非線性問題比線性問題所需要的計算機資源要高許多倍；

非線性分析步被分為許多增量步。ABAQUS通過迭代，在新的載荷增量結(jié)束時近似地達(dá)到靜力學(xué)平衡。ABAQUS在整個模擬計算中完全控制載荷的增量和收斂性；

狀態(tài)文件允許在分析運行時監(jiān)控分析過程的進(jìn)展。信息文件包含了載荷增量和迭代過程的詳細(xì)信息；

在每個增量步結(jié)束時可以保存計算結(jié)果，這樣結(jié)構(gòu)響應(yīng)的演化就可以用ABAQUS/Post顯示出來。計算結(jié)果也可以用x-y圖的形式繪出。

8．材料

小結(jié)：

ABAQUS包含一個廣泛的材料庫，可模擬各種工程材料的性質(zhì)。其中包括金屬塑性和橡膠彈性模型；

金屬塑性模型的應(yīng)力－應(yīng)變數(shù)據(jù)必須用真實應(yīng)變定義； 金屬塑性模型假定材料具有一旦屈服即不可壓縮的性質(zhì)。這將對應(yīng)用于彈－塑性模擬的單元類型帶來某些限制；

多項式和奧根應(yīng)變能函數(shù)可應(yīng)用于橡膠材料的彈性（超彈性）。兩種模型均允許直接用實驗數(shù)據(jù)來確定材料的系數(shù)。實驗數(shù)據(jù)必須是名義應(yīng)力和名義應(yīng)變的值；

在超彈性材料模型中的穩(wěn)定性警告，說明所要分析的應(yīng)變范圍不合適； 存在對稱性時，可以只考慮部分模型從而減小模擬的尺寸?？赏ㄟ^施加適當(dāng)?shù)倪吔鐥l件來反映結(jié)構(gòu)其余部分的效應(yīng)；

大畸變問題的網(wǎng)格設(shè)計比小位移問題更加困難。在分析的任何階段，網(wǎng)格中的單元務(wù)必不能過于畸變；

ABAQUS/Post中的*DEFINE CURVE命令允許處理曲線上的數(shù)據(jù)以生成新的曲線。兩條曲線或一條曲線與一個常數(shù)可以加、減、乘、除。曲線還可以求導(dǎo)、積分和合并。

9．動態(tài)問題

具有下列特征的問題適于采用線性瞬態(tài)動力學(xué)分析：

系統(tǒng)是線性的：線性材料行為，無接觸條件，無非線性的幾何效應(yīng)； 響應(yīng)只受較少的頻率支配。當(dāng)響應(yīng)中各頻率成分增加時，例如撞擊和沖擊情況，振型疊加方法的效果將大大降低；

載荷的主要頻率在可得到的固有頻率范圍內(nèi)，以確保對載荷的描述足夠精確；

由于任何突然加載所產(chǎn)生的初始加速度能用特征模型精確描述； 對系統(tǒng)的阻尼不能過大。小結(jié)：

動態(tài)分析包括結(jié)構(gòu)的慣性效應(yīng)；

*FREQUENCY可以計算結(jié)構(gòu)的固有頻率和振型；

通過振型疊加，可以確定線性系統(tǒng)的動態(tài)響應(yīng)。這一方法盡管有效，但是不能用于非線性問題；

線性動態(tài)過程可以計算瞬態(tài)載荷的瞬態(tài)響應(yīng)、諧振動下的穩(wěn)態(tài)響應(yīng)、支座移動造成的響應(yīng)峰值和隨機載荷的響應(yīng)；

為了準(zhǔn)確表示結(jié)構(gòu)的動態(tài)行為，必須選擇足夠多的振型?？偟牡刃Ｐ唾|(zhì)量應(yīng)占可動質(zhì)量的90%以上；

用戶可以給定直接模態(tài)阻尼、瑞利阻尼和復(fù)合模態(tài)阻尼。但是由于固有頻率和振型的計算都是基于無阻尼的結(jié)構(gòu)，所以此法只適用于低阻尼結(jié)構(gòu)；

模態(tài)技術(shù)不適用于非線性的動態(tài)模擬。在這種情況下必須采用自己的時間積分方法（*DYNAMIC）

*AMPLITUDE選項可以描述隨時間任意變化的載荷，以及給定的邊界條件； 振型和瞬態(tài)結(jié)果可以在ABAQUS/Post中用動畫顯示。這對于理解動態(tài)響應(yīng)和非線性靜態(tài)分析十分有幫助。

10．多步驟分析

小結(jié)：

一個ABAQUS模擬過程可以包含任意數(shù)目的步驟；

一個分析步驟就是一段“時間”，在這段時間里ABAQUS計算模型對一套指定載荷和邊界條件的響應(yīng)。這一步驟中所用的特殊分析過程確定了這個響應(yīng)的特征；

在一個一般分析步驟中，結(jié)構(gòu)的響應(yīng)可能是線性的，也可能是非線性的； 每一個一般步驟的開始狀態(tài)是上一個一般步驟的結(jié)束狀態(tài)。這樣，在一個模擬中模型的響應(yīng)隨一系列一般步驟而演化；

線性擾動步驟計算結(jié)構(gòu)對擾動載荷的線性響應(yīng)。這個響應(yīng)的基本狀態(tài)是相對于最后一個一般步驟結(jié)束時模型的狀態(tài)所定義的；

在一般步驟中任何載荷選項里的OP參數(shù)（例如*BOUNDARY,*CLOAD和*DLOAD中）控制著這些選項中所指定的數(shù)值是如何與前面步驟中定義的數(shù)值相互作用的；

只要存儲了一個重新啟動文件就可以進(jìn)行重新啟動分析。重新啟動文件可以用來繼續(xù)一個中斷的分析或者給模擬添加附加的載荷過程。11．接觸

小結(jié)：

接觸分析需要一個謹(jǐn)慎的邏輯方法。如果必要，將分析分解成幾步執(zhí)行，并緩慢地施加荷載，以保證很好地建立接觸條件；

一般地，對分析的每一步最好采用分離步驟進(jìn)行，即使只是因為載荷而改變邊界條件。您幾乎肯定要比預(yù)期情況應(yīng)用更多的步驟，但模型則收斂得更容易。如果想一步就將所有的載荷加上，接觸分析是難以完成的；

在對結(jié)構(gòu)施加工作載荷之前，要在所有部件之間取得穩(wěn)定的接觸條件。如果必要，采用臨時的邊界條件，在以后階段再消除這些約束。只要所提供的約束不產(chǎn)生永久的變形，對最后的結(jié)果應(yīng)該毫無影響；

不用對接觸面上的節(jié)點施加邊界條件，即在接觸方向上限制節(jié)點。如果有摩擦，不要在任何自由度上約束這些節(jié)點：可能導(dǎo)致零主元信息；

對于接觸模擬，總要試圖使用一階單元。

第二篇：Abaqus遇到的問題小結(jié)

Abaqus遇到的問題小結(jié)

1.Abaqus量綱系統(tǒng)（SI）：m、N、Kg、s、Pa、J、Kg/m。

2.有限元單元類型（Element Type）：實體單元、殼單元、桿件單元（梁Beam、桁架Truss）3.Abaqus怎樣將部件的不同部分設(shè)定為不同的子集？

3tools >partition cell/plane（如果是要分層使用，則進(jìn)行分割，分別賦予不同的屬性;或者對其中的單元進(jìn)行設(shè)置Set；如果要在運算中替換則需要在關(guān)鍵詞里添加語句定義）

4.建模的過程中輸入的尺寸被自動降低精度（四舍五入），請問如何設(shè)置尺寸精度（小數(shù)點位數(shù)）？

在 sketch option里打開sketch模塊(左下角最下面)打開sketchoption對話框,其中decimal places 表示小數(shù)點位數(shù),默認(rèn)為2,最高能調(diào)到6。

5.abaqus中出現(xiàn)，Dependent partinstancesconnot be edited，怎么辦？

mesh on part , mesh on instance一個非獨立實體只是原始部件的一個指針,可以對原始部件劃分網(wǎng)格，但是不能對一個非獨立實體劃分網(wǎng)格，即mesh on part。

方法：左邊模型樹里面，Assembly->instances里面的子選項右擊，點Make independent;或者點擊model旁邊的Part,對part進(jìn)行網(wǎng)格劃分。

6.材料屬性中拉伸強度和壓縮強度怎么定義？

計算結(jié)果有每個計算節(jié)點的應(yīng)力，計算中用不到強度值（彈性力學(xué)三組基本方程不涉及強度，求解過程用不到）。（強度值是用來校核的：斷裂，屈服或超過線性變形）

7.如何選擇創(chuàng)建獨立實體還是非獨立實體？

如果集合中包含許多具有相同性質(zhì)的部件，則創(chuàng)建非獨立實體（Dependent）。反之，如果集合中包含許多不同性質(zhì)的部件，則創(chuàng)建獨立實體（Independent）要有優(yōu)勢些。

8.地基土的分層問題？

方法1：定義為一個part，然后partition分層賦材料屬性，不需要設(shè)置接觸。（超級贊?。。┓椒?：使用多個part，每個part單獨定義屬性。一般不用設(shè)置一些tie，各層間只有彈性模量相差較大時再在各層間定義tie。

9.tie捆綁約束。

tie功能很強大，可以模擬實體與實體，實體與殼體，殼體和殼體的連接，而且允許連接對象之間有距離。

在tie面板中，Position Tolerance位置容差是一個重要的參數(shù)。一般選擇默認(rèn)設(shè)置usecomputed default，abaqus會根據(jù)你要tie的對象的網(wǎng)格尺寸，自動選取個容差值，一般為網(wǎng)格基本尺寸5%。于是，連接對象之間，距離小于這5%容差值的節(jié)點便施加tie約束，也就是對應(yīng)的兩節(jié)點的自由度擁有一樣的數(shù)值。

如果連接對象之間有距離怎么辦？若還是用usecomputed default，則是無法起到tie約束效果的。這就要手動設(shè)置specify distance為一個適合的參數(shù)。這個值要大于連接對象對應(yīng)節(jié)點之間的距離。若不是很確定，可以先設(shè)一個較大的值，計算完后在Visualization調(diào)出Node sets,觀察tie實際起作用的節(jié)點區(qū)域。再對specify distance值進(jìn)行調(diào)整。

在使用specify distance,需將Adjust slave surface initial position選項取消。否則abaqus會將兩連接對象拉近，網(wǎng)格產(chǎn)生畸變，產(chǎn)生錯誤，中止計算！

10.abaqus怎么將模型分割畫網(wǎng)格？

劃分模型的時候有全局劃分網(wǎng)格按鈕，旁邊有個局部劃分網(wǎng)格按鈕，（經(jīng)常要用到partition技術(shù)），點它給需要另外劃分的部分部種子，然后下面有按區(qū)域劃分網(wǎng)格。全局劃分網(wǎng)格和局部劃分網(wǎng)格分別是mesh part和mesh region。

網(wǎng)格劃分完畢之后，最好在verify中對其中質(zhì)量不好的單元進(jìn)行編輯處理，避免在job模塊里面出現(xiàn)錯誤或警告信息。

11.邊界條件設(shè)置在初始步還是后續(xù)的分析步中？

設(shè)置在initial。

12.interaction面選擇不上？

在此模塊里面，一般先定義相應(yīng)的surfaces，并分別賦予有意義的容易識別的名字，如slavesurface、master surface，并選擇接觸面的正確方向，如果方向選擇錯誤的話，經(jīng)常在job模塊中會出現(xiàn)無法收斂的錯誤信息。

13.通過ABAQUS計算得出的應(yīng)力值不連續(xù)是什么原因造成的?

這主要是因為ABAQUS的應(yīng)力和應(yīng)變等單元變量值(相對位移等節(jié)點變量)的計算方式所引起的 這些單元變量值主要是通過單元內(nèi)的高斯積分點計算得出，然后再通過外插值法計算出節(jié)點上的應(yīng)力應(yīng)變值，因此對于兩個相鄰單元共同占用的節(jié)點，從兩個單元外插得出的值有時是不一樣的。對于大多情況．如果網(wǎng)格足夠密，高斯積分點上的值不會相差很大．所以可得到一個連續(xù)性比較好的結(jié)果另外可以通過設(shè)定平均值的臨界值來決定contour的連續(xù)性(即legend下默認(rèn)的75％)。比如說，如果同一節(jié)點扶不同單元得出的外插值相差75％以下(即較大的外插值減去較小的外插值的差乘以大值)，那么contour輸出這兩個值的平均值。

同時需要注意的是，這個默認(rèn)的75％ 可以修改，如果設(shè)為0％．那么整個contour的單元間幾乎都不連續(xù)，如果設(shè)在100％，那么整個圖肯定是連續(xù)的。此外，還可以通過設(shè)定legend下的limit方式(即外插值或者是平均值)來控制最后的結(jié)果。

14.平衡地應(yīng)力是出現(xiàn)An initial condition has been specified on element 0 but thiselement has not been defined？

是因為修改inp文件時，文件中有未刪除的空行。

15.土體初始地應(yīng)力平衡方法。

地應(yīng)力平衡原理：平衡原理和屈服條件。目的：為了使數(shù)值模擬獲得一個存在初始應(yīng)力，而無初始應(yīng)變的狀態(tài)。

平衡地應(yīng)力方法：

1.Autobanlance 2.*Initial Conditions,type=stress,Geostatic;（關(guān)鍵字定義初始地應(yīng)力，只適合土體表面水平的土體，該方法需給出不同材料區(qū)域的最高和最低點的自重應(yīng)力及其相應(yīng)坐標(biāo)）

3.*Initial Conditions,type=stress,input=xxx.csv;（初始地應(yīng)力提取法，此方法最為通用，可用于不同材料、不同地形，適應(yīng)性強。）

方法3的具體步驟： i.ii.iii.iv.v.在Abaqus中建立完整的模型，定義部件、截面屬性、裝配件、網(wǎng)格、邊界條件等。單元類型使用的默認(rèn)的線性縮減積分單元C3D8R。

在Step功能模塊中創(chuàng)建分析步，將類型設(shè)為Geostatic。在Load功能模塊中，為整個土體定義重力載荷（類型為Geostatic），并在Component3中輸入-9.8。

在命令行輸入：mdb.models['Model-1'].setValues(noPartsInputFile=ON);在Job模塊中創(chuàng)建名為Job-NoInitialCondition的分析作業(yè)，提交分析。

將分析得到的應(yīng)力場保存為一個文本文件。具體的方法為：打開分析得到的ODB文件Job-NoInitialCondition.odb，選擇菜單Report---FieldOutput，選中Variable對話框中積分點上的各個應(yīng)力分量（S11、S22、S33、S12、S13、S23）；在Setup標(biāo)簽頁，創(chuàng)建inp輸出文件，并命名為XX.inp，不選Append to file和Column total 和column min/max。

用Excel修改輸出的XX.Inp文件，刪除inp文件中的所有敘述性文字和空格（例如若土體有N層，需找出N-1處文字性敘述的地方并刪除），只保留數(shù)字結(jié)果，將數(shù)字間在順序不變的前提下用逗號隔開，保存成csv以逗號隔開的文件形式XXX.csv。vi.vii.viii.為模型定義初始應(yīng)力場。因在Abaqus中無法直接定義初始應(yīng)力，只能手工添加關(guān)鍵詞，具體做法：選中菜單Model---Edit Keywords，在*Boundary語句之后，*Step語句之前添加語句*Initial Conditions,type=stress,input=XXX.csv 在Job功能模塊中定義Job-WithInitialCondition作業(yè)，提交分析，即得到所需要的結(jié)果。

16.Part模塊建立剛體時一定要指定參考點referencepoint，在后面分析過程中的位移、載荷的施加都在施加在此參考點上的。剛體只有平動和轉(zhuǎn)動，不會有變形，因此參考點的選擇即可以在建立的rigid part上，也可以建立在其他你覺得方便的位置。

17.load模塊，是施加載荷和位移邊界條件的，先對要施加載荷和邊界條件的面、線、節(jié)點等建立set，再對相應(yīng)sets分別設(shè)置。在施加載荷的時候，注意載荷的方向性，另外此模塊也可以施加自定義載荷、predifined field等。

18.abaqus結(jié)果顯示時數(shù)據(jù)太小看不清楚怎么辦？

Viewport-->Viewport Annotation Options..在出來的對話框中選擇“Legend”標(biāo)簽，然后點擊“Set Front...” 之后選擇字體大小，確認(rèn)。

BY

HIT-WGJ

第三篇：ABAQUS與FORTRAN論壇問題小結(jié)

ABAQUS與FORTRAN論壇問題小結(jié)

Q:索了一下論壇發(fā)現(xiàn)以前的問題都是不同類型的子程序，如UEL, UMAT...這樣把它們放在一個.for文件里不會有誤會，但是如果都是UEL的話，該怎么識別呢？ A：subroutine Umat()IF(CMNAME(1:4).EQ.'MAT1')THENCALL UMAT_MAT1(argument_list)ELSE IF(CMNAME(1:4).EQ.'MAT2')THENCALL UMAT_MAT2(argument_list)END IF上面是區(qū)分不同umat的方法，區(qū)分uel應(yīng)該類似吧

Q:行abaqus的時候出現(xiàn)這樣的提示“Problem during compilation-df.exe not found in PATH.”請問這是什么原因.A: 就是說找不到用于編譯的df.exe文件（好像是Fortran的編譯運行文件），在Windows控制面板的“系統(tǒng)”里，打開高級屬性，將系統(tǒng)的PATH添加上Fortran的路徑，最好也一并把INCLUDE和LIB也都添加上對應(yīng)的FORTRAN的路徑。

Q: 如何在ABAQUS中調(diào)用用戶子程序UMAT

A：abaqus job=*** user=*.for or *.objA: 裝好fortran編譯器的前提下，有兩種方式：

1、命令行:abaqus job=*.inp user=*.for或 abaqus job=*.inp user=*.obj(可以由*.f90編譯后得到)

2、菜單格式：先運行CAE，在general設(shè)置中可以指定Umat用戶子程序A: job=jobname user=yourfortran.for(其它類型子程序)A：如果是本構(gòu)的二次開發(fā)：*material, user

Q: abaqus job=job-name user={source-file | object-file}source-file和object-file各是什么？ A: fortran的源文件（.for）或者目標(biāo)文件（.o）在command下面輸入以上命令行就可以正常運行.用戶子程序再調(diào)用別的fortran程序時，要加入abaqus 子程序中的include 'ABA_PARAM.INC'。

Q: 自己編寫子程序需要哪些工作步驟 A: 1.建立工作目錄

2.將abaqus安裝目錄6.4-pr11site下的aba_param_dp.inc 或 aba_param_sp.inc拷貝到工作目錄，并改名為aba_param.inc ； 3.將編好的fortran程序拷貝到工作目錄； 4.編譯fortran程序，將.obj文件拷貝到工作目錄； 5.建立好輸入文件.inp;6.運行abaqus job=inp file name user=fortan name即可。我的要運行abq64pr11 job=inp file name user=fortan name我也學(xué)了沒多久，按照這個過程成功了一個！你們試一下吧！Q: 想用FORTRAN實現(xiàn)：...result = SYSTEMQQ('abq job=e:abaqu***amplesbeam interactive')...出錯提示：ABAQUS ERROR：Unable to create a temporary directory C:DOCUME~1ADMINI~1LOCALS~1TEMPAdministrato_ e:abaqu***amplesbeam_3628 to be used for scratch files.ABAQUS/Analysis exited with errors.但如果是 result = SYSTEMQQ('abq job=beam interactive')，但需要在FORTRAN彈出的DOS窗口里手工輸入文件的絕對路徑則計算正常。但我的目的是實現(xiàn)Fortran程序?qū)BAQUS的計算程序的自動調(diào)用。

A: 說個投機的花招，有沒有試過在前面加上result = SYSTEMQQ('path=e:/abaqu***amples')這樣的命令~~方法簡單，不過能解決問題最好

A：把要算的工況都集中在一個文件夾下，fortran生成.exe文件，放在該文件夾下，然后雙擊運行。計算結(jié)束之后再把各種工況分類。笨招，但好用。

A: 我想，或者把systemQQ改一下，先cd e:...,然后再systemQQ(abaqus job=jobname interactive)。

Q: 偶想在FORTRAN里設(shè)定一個計時器，其功能就是每隔一段時間（如10s）就去讀一次ABAQUS的(.fil)文件。哪個函數(shù)可以實現(xiàn)這個功能？怎么實現(xiàn)？？ A: 通過一個while循環(huán)和下面的語句：call sleepqq（時間間隔）就可以了。

Q: 本人嘗試使用ABAQUS的子程序SUBROUTINE，根據(jù)幫助文件需要編譯通過生成*.for和*.lib文件，可是編譯的時候出現(xiàn)錯誤：Error: Cannot open include file 'ABA_PARAM.INC'不知道是什么原因，怎么解決？

A: 1.確認(rèn)Visual Fortran及VC++正確安裝，并在安裝ABAQUS之前安裝。2.首先將用戶子程序*.f文件后綴名改為*.for。3.到控制面板－》系統(tǒng)－》高級－》環(huán)境變量中的 Administrator用戶變量:修改—— 變量名：lib變量值:C:Program FilesMicrosoft Visual StudioDF98LIB;C:Program FilesMicrosoft Visual StudioVC98LIB（我的問題就是因為只加了一個LIB，必須兩個都加！）完成以上步驟以后，就可以運行了。

Q: 怎樣把自己編寫的子程序倒入到*.INP文件中。

A: Using the following command : abaqus job=job-id As can be seen this is no different from running a standard abaqus job.The user subroutine itself can be embedded in the abaqus input file.Here it is illustrated with the umat subroutine.<....part of the abaqus input file....>

........*END STEP*USER SUBROUTINES SUBROUTINE UMAT(........)

........END

........Alternatively the user subroutine can be in a separate file(say my_material.f)and the INPUT parameter is set to that file name.<....part of the abaqus input file....>

........*USER SUBROUTINES, INPUT=my_material.f

........Q: 用Fortran SYSTEMQQ函數(shù)作批處理遇到問題，在fortran中使用result = systemqq(abaqus job=job1 interactive)作批處理，總是在計算兩個job之后就不繼續(xù)往下算了，cmd明令行提示：run time error:floating-point...請問這是怎么回事？

A: 個方法就是做一個批處理文件，abaqus job=beam interactiveabaqus job=beam1 interactiveabaqus job=beam2 interactiveabaqus job=beam n? interactive存為.bat文件即可。python好像一般是做CAE前后處理界面的。

A：前面加call.call abaqus job=1 interactivecall abaqus job=2 interactivecall abaqus job=3 interactive

Q:FORTRAN如何調(diào)用ABAQUS？

A: use the function “system” or “systemqq” USE DFLIB LOGICAL(4)result result = SYSTEMQQ('abaqus job=xx')QA: 首先我用ABAQUS的VARIFY命令對系統(tǒng)進(jìn)行修正和檢測,這個命令會自動運行,對ABAQUS進(jìn)行調(diào)試等操作,其中涉及到對用戶子程序的模擬計算,如果計算失敗,就會在VARIFY.LOG文件中說明,確實發(fā)現(xiàn)“USER SUBROUTINE...和SINGLE...和DOUBLE...”失敗,說明用戶子程序安裝不成功,肯定與VISUAL FORTRAN6.5有關(guān)!然后,找到FORTRAN_MAKE文本文件,發(fā)現(xiàn)錯誤,可能是與FORTRAN安裝有關(guān)!于是卸載FORTRAN,重新安裝在C盤!安裝完畢后,改變環(huán)境變量如下: INCLUDE C:Program FilesMicrosoft Visual StudioDF98IMSLINCLUDE;C:Program FilesMicrosoft Visual StudioVC98INCLUDE;D:程序文件VC98INCLUDE;LIB C:Program FilesMicrosoft Visual StudioDF98IMSLLIB;C:Program FilesMicrosoft Visual StudioDF98LIB;C:Program FilesMicrosoft Visual StudioVC98LIB;D:程序文件VC98LIB;PATH C:Program FilesMicrosoft Visual StudioCommonTools;D:程序文件VC98BIN;%SystemRoot%system32;%SystemRoot%;%SystemRoot%System32Wbem;C:KAV6;D:程序文件KV6;c:program filesmicrosoft visual studiodf98bin;c:program filesmicrosoft visual studiovc98bin;c:program filesmicrosoft visual studiocommonmsdev98bin;D:ABAQUSCommands;C:Program FilesMicrosoft Visual StudioDF98BIN我們只需要關(guān)注MICROSOFT VISUAL STUDIO文件夾!最后,再次執(zhí)行VARIFY命令,等運行完畢后,檢查VARIFY文本文件,發(fā)現(xiàn): ABAQUS/Standard with User Subroutines...PASSContinuing...ABAQUS/Explicit with User Subroutines(single precision)...PASSContinuing...ABAQUS/Explicit with User Subroutines(double precision)...PASSContinuing...OK?。。。。。。。?！搞定

Q: 因為ABAQUS是用PYTHON編寫的，那么我想要PYTHON 的原代碼，因為我想 在其中加一個橢圓程序，不知可行么？（訪問原代碼問題）

A: 所有的PY都是差不多一樣的,把ABAQUS的PY拿到MARC試試看

文檔來源：http://004km.cn/bbs/redirect.php?fid=10&tid=9031&goto=nextnewset

第四篇：abaqus在基坑工程中使用小結(jié)

abaqus在基坑工程中使用小結(jié)

本人在使用abaqus解決基坑支護問題的過程中，得到了版上aba_aba，以及ustc的life大俠等人的幫助，寫些自己使用abaqus的一些感受吧，（其實大多數(shù)東西都是從版上得到的相關(guān)參考）希望對其他新手能有所幫助。**首先是建模方式

對于初學(xué)者比較好掌握的是CAE界面操作方式。但是由于ABAQUS的CAE模塊和輸入文件兩種建模方式是由兩家不同的公司研制的，CAE模塊推出比較晚，所以功能還不是很完全，對于CAE和INP交互操作需要注意幾點：

1、一些功能只能通過編輯INP輸入文件來實現(xiàn)，比如模擬開挖建造時需要用到殺死或者激活單元的語句model change，此時需要直接在CAE下編輯關(guān)鍵字，或者生成INP文件后通過記事本等工具打開INP文件后加入關(guān)鍵字。

2、CAE對于INP輸入文件中的一些關(guān)鍵字不支持，如果在CAE下編輯，則系統(tǒng)會給出不支持關(guān)鍵字的提示。對于這樣的情況則只能通過文本編輯工具打開INP來操作。

3、對INP輸入文件進(jìn)行編輯的方式能支持所有關(guān)鍵字功能，因此應(yīng)用更加廣泛。在許多情況下都需要對INP文件進(jìn)行操作。**然后是初始應(yīng)力場的施加。

土工有限元分析中初始應(yīng)力場的施加是個關(guān)鍵。在ABAQUS中，施加初始應(yīng)力場很方便，可以使用以下語句來完成。

initial conditions,type=stress,geostatic 可以有兩種實現(xiàn)的辦法：

1、如果地表水平，且初始應(yīng)力就是自重應(yīng)力時，可以在數(shù)據(jù)行中按照如下順序填寫數(shù)據(jù) 施加初始應(yīng)力的集合，應(yīng)力值1，應(yīng)力值1對應(yīng)的坐標(biāo)值1，應(yīng)力值2，應(yīng)力值1對應(yīng)的坐標(biāo)值2，x方向側(cè)壓力系數(shù)，y方向側(cè)壓力系數(shù)（在二維問題中可以省略）

這樣的順序定義初始應(yīng)力實際上默認(rèn)了初始應(yīng)力就是自重應(yīng)力，只需要給定集合中2點的應(yīng)力值及其對應(yīng)的坐標(biāo)值，其他部位的應(yīng)力值就可以通過1，2兩點的連線采用差值來得到，而水平向的地應(yīng)力則乘一個側(cè)壓力系數(shù)得到。（這些內(nèi)容在abaqus手冊中都有）

2、采用input文件的方式，即先通過計算一個重力荷載作用下沒有初始應(yīng)力的模型而得到各點的應(yīng)力值，然后將這些值導(dǎo)入到新的模型中，在*initial conditions,type=stress,后面加入input=XXX.XXX，讓模型接受存有初始應(yīng)力值的input文件而施加。這樣的方法通用性非常強，但是比1稍微復(fù)雜些，具體可以見相關(guān)手冊。

另外，對于施加初始應(yīng)力，對于二維的情況，系統(tǒng)默認(rèn)地應(yīng)力方向為Y方向，對于三維情況，默認(rèn)為Z方向，這一點在初學(xué)者中容易弄錯。**模擬開挖建造

ABAQUS中模擬開挖建造通過關(guān)鍵字model change,type=element來實現(xiàn)，如果是殺死單元則在后面加上remove,如果激活單元則加上add。

事先需要將需要殺死或激活的土體單元定義為一個集合set。整個應(yīng)用格式如下： *Model change, type=element, remove Set name1，set name2

其意義如下：首先Model change說明是模型發(fā)生變化，后面是種類，因為除了殺死激活單元外，還常用到殺死激活接觸對(type=contact)，此處需要指明是殺死激活單元導(dǎo)致的模型變化，后面remove表示殺死單元。重啟一行寫出需要殺死激活的單元集合列表。

在模擬建造過程中，通常一個單元的材料屬性會發(fā)生變化，比如土釘支護中面層處單元，在沒有開挖時是土體材料，當(dāng)開挖建造之后就變成了面層的屬性了。在ABAQUS中模擬材料屬性隨開挖步進(jìn)行可以使用場量*Field關(guān)鍵字來實現(xiàn)。具體過程如下：首先是材料定義上要指定材料性質(zhì)隨著場變量而變化，如下面語句中，定義彈性時，利用dependencies=1表示材料根據(jù)1個場變量變化而變化，在接下來的屬性數(shù)據(jù)行中每行后面的數(shù)字1，2表示了該行屬性對應(yīng)的這個場變量的值，如第一行的彈性模量和泊松比值是場變量為1的值，而第二行的彈性模量和泊松比值是場變量為2的值。*Elastic, dependencies=1

4e+07, 0.35, , 1.1.5e+10, 0.2, , 2.接下來就是在每個分析步中利用*field,variable=1這樣語句指定分析步中使用的場變量數(shù)量，本例中值為1個場變量。接下來的數(shù)據(jù)行則是表示使用場變量的集合名稱及其在該分析步中的場變量值。如果這個值是1，根據(jù)在定義材料時的場變量定義，值為1時對應(yīng)于材料屬性第一行值，這樣就實現(xiàn)了材料屬性隨分析步變化的功能。

**模擬地連墻與土的接觸ABAQUS中對于接觸的模擬可以通過定義接觸對來實現(xiàn)。系統(tǒng)可以允許面對面、點對面的接觸等多種接觸形式。在模擬開挖過程中，開挖土體與墻體的接觸屬性會喪失，根據(jù)經(jīng)驗，最好能將接觸按照開挖步驟劃分為許多對接觸對（所謂的細(xì)化），根據(jù)開挖的進(jìn)行，將接觸對先殺死，這樣系統(tǒng)比較容易找到準(zhǔn)確的接觸狀態(tài)，使接觸模擬順利進(jìn)行。模擬接觸時，定義接觸對要分清主控面和受控面。主控面和受控面的選擇一般遵循以下原則：

一、材質(zhì)較軟的為受控面，材質(zhì)較硬的為主控面。

二、網(wǎng)格較密的為受控面，網(wǎng)格稀疏的為主控面。應(yīng)用到網(wǎng)格劃分的時候就要先根據(jù)材質(zhì)的相對軟硬程度確定主控面，以便在網(wǎng)格劃分中按照網(wǎng)格疏密關(guān)系正確劃分網(wǎng)格。主控面和受控面定義一定要正確，這是影響接觸模擬的關(guān)鍵。如果一個節(jié)點上同時定義了兩個接觸面，或者同時定義了一個接觸面和邊界條件，系統(tǒng)則會產(chǎn)生過約束(over constraint)的提示，很容易影響到接觸的模擬。

分享<隧道施工應(yīng)力釋放法的ABAQUS實現(xiàn)【inp及代碼見2#】> 應(yīng)力釋放法模擬隧道開挖主要有兩種實現(xiàn)途徑，大家可以參考帖子： http://forum.simwe.com/thread-765338-1-1.html

簡單的說，主要就是參數(shù)弱化法和施加節(jié)點反力法。這兩種方法各有優(yōu)劣。參數(shù)弱化法：

優(yōu)點：

1、便于實現(xiàn)；

2、便于控制應(yīng)力釋放過程；

缺點：

1、無法給出明確的應(yīng)力釋放率；

2、無法用在劍橋模型中，因為找不到一個合適的參數(shù)進(jìn)行弱化。施加節(jié)點反力法：

優(yōu)點：

1、能夠明確給出應(yīng)力釋放率；

2、物理意義較為明確； 缺點：實現(xiàn)過程不是很方便

下面跟大家匯報一下我的做法，希望有經(jīng)驗的童鞋給出建議，各種板磚，放心拍吧~

一、參數(shù)弱化法

參數(shù)弱化方法其實沒有什么難度，版內(nèi)已經(jīng)討論了無數(shù)次了。各種牛們也是不厭其煩的跟新手們說了一遍又一遍，可惜有些童鞋天生喜歡求助，上來就發(fā)帖問，完全置搜索功能于不顧。所以我在這要再次呼吁，搜索處要更加醒目，加紅字，加粗，加浮動，另外加旁注：親，如果你想快速入門aba，搜索吧，這是悶騷技術(shù)人新一代的學(xué)習(xí)方法，包會哦！

參數(shù)弱化，主要就是弱化彈性模量，這個可以通過關(guān)鍵字*field來實現(xiàn)，具體腫么實現(xiàn)，搜索吧。下面給出計算后的云圖，拋磚引玉。地應(yīng)力平衡，勉強可以吧：

待開挖區(qū)域彈性模量折減20%：

如何添加襯砌以及如何設(shè)置追蹤單元我想再單開一貼和大家進(jìn)行總結(jié)探討。上面位移矢量圖隱去了待開挖部分，是為了矢量顯示清晰，實際上這一步待開挖（中心土體）部分實際還是存在的，只不過彈性模量折減過了。從圖中可以看出，這種方法計算出的位移場并不是很符合實際情況。隧道開挖土體位移場可以參見帖子：http://forum.simwe.com/thread-985421-1-1.html

在這里也順便給這個帖子做個個人的解答吧。關(guān)于該貼提出的強制位移法，我想在后面專開一貼和大家討論總結(jié)。

有經(jīng)驗的童鞋可否告知你們采用折減法做出的位移場是否也和上面一樣的？

二、施加節(jié)點反力法

這一方法的技術(shù)難點在于，如何獲得地應(yīng)力平衡后隧道周邊土體節(jié)點反力。通過搜索和查看manual，基本可以鎖定兩種輸出量：NFORC & RF。其中，RF只能在有邊界約束處輸出，NFORC根據(jù)manual我猜測應(yīng)該是由節(jié)點所涉及單元的應(yīng)力進(jìn)行外插平均后得到的節(jié)點力，參見帖子：http://forum.simwe.com/thread-883192-1-1.html10#。不管怎樣，總的思路是，通過對開挖后隧道相鄰?fù)馏w節(jié)點施加非均布節(jié)點力，使得此時的模型（土體）保持初始地應(yīng)力和位移準(zhǔn)零狀態(tài)，然后施加原節(jié)點荷載*（1-應(yīng)力釋放率），以此來模擬不同的應(yīng)力釋放率下隧道及土體的力學(xué)響應(yīng)。

以前的帖子對這個問題討論了很多，有代表性的列舉如下（話說論壇的收藏是個好功能）：

http://forum.simwe.com/thread-872620-1-1.html 該貼lz給出了四種獲取節(jié)點力的方法，都很有參考價值。最后采用了NFORC，但得到的效果不是很理想。我做了一下，發(fā)現(xiàn)的確如該lz所說，NFORC輸出的節(jié)點力是不能直接施加到原節(jié)點上的，即使乘以2倍，也只能是接近實際的節(jié)點力。所得到的位移場在初始狀態(tài)就有很大的位移，不符合前面提出的總體思路。該貼5#和7#都提出采用約束洞周位移獲取RF。事實證明，這個方法是可行的，但不同的操作方式得到的結(jié)果千差萬別。該貼7#提出的思路是對的，但給出的文件中將洞周位移約束貫穿了整個施工過程，這樣造成的后果是，其計算起始點已經(jīng)不再是green field的情況，得到的應(yīng)力場自然不符合實際情況，那么在錯誤的應(yīng)力場上做的計算怎么能對呢？

http://forum.simwe.com/thread-854591-1-1.html

該貼給出了施加洞周約束獲取反力的正確思路，即施加洞周位移約束要在地應(yīng)力平衡后進(jìn)行。這樣能保證施加位移約束時的地應(yīng)力初始狀態(tài)。下面就我做的例子詳細(xì)說明。

模型：

材料參數(shù)：

修正劍橋模型

*Material, name=soil-3 *Clay Plasticity, intercept=1.45 0.11, 1.27, 0., 1., 1., , 1.*Density 1770., *Porous Elastic 0.009, 0.32, 0.劍橋模型參數(shù)的設(shè)置參見：http://forum.simwe.com/thread-968882-1-1.html

這里再啰嗦幾句，最近很多人在問為什么用劍橋模型的時候會提示“THE SUM OF INITIAL PRESSUE AND ELASTIC TENSILE STRENGTH FOR THE POROUS ELASTIC MATERIAL IN ELEMENT..."，如果你確定有地應(yīng)力平衡，那么就請你檢查是否是用的直接加重力平衡的吧。多孔介質(zhì)彈性模型是需要在initial conditions中定義初始地應(yīng)力的。

地應(yīng)力平衡：

平衡效果一般，主要是做實例，只平衡了一次，到了10E-4就沒再管了。

施加洞周邊界約束，同時remove掉待開挖部分

可以看出，這一步進(jìn)行完后，剩余土體的應(yīng)力基本沒有發(fā)生改變，也就是基本保持了初始地應(yīng)力，在這種情況下導(dǎo)出節(jié)點力再倒入折減后支護力才是正確的?？梢钥闯鍪┘佣粗芪灰萍s束后土體位移分布發(fā)生了變化，但數(shù)量級還是很小，基本沒變，可以忽略。最后一張圖是洞周節(jié)點反力云圖，在這種情況下就可以輸出節(jié)點反力了。下面介紹一下我的做法。

第五篇：abaqus自我總結(jié)

如何輸入inp文件

inelastic heat fraction：非彈性功（一般是塑性功）轉(zhuǎn)換為熱量的比例，比例越大，產(chǎn)熱越高，溫升越高

rate dependent：本構(gòu)率相關(guān)參數(shù)，表示材料的屈服流動應(yīng)力不僅和等效塑性應(yīng)變相關(guān)，而且還與應(yīng)變率有關(guān)，一般是正相關(guān)

如何將二維圖表的背景變成白色

我的怎么就改變了呢 說了雙擊

舉例子說說我的理解吧：history output->盯住一個節(jié)點，這個節(jié)點在整個加載時間段內(nèi)應(yīng)力變化情況，用xy圖表示就是一條TIME-STRESS曲線；field output->盯住一個時刻，所有節(jié)點在這個時刻時的應(yīng)力值，直觀表示是云圖。當(dāng)然在xyplot中你可以選擇盯住一個節(jié)點按時間輸出該節(jié)點的力（即ODB field output），不過最多只能在各個幀對應(yīng)的時刻輸出，xy圖中點的數(shù)量較少。如果我錯了的話請大家指正哈

view---overlayplot--然后每一個時刻創(chuàng)建一個新圖層，然后plotoverlay

怎么樣才能使圓柱體繞另一固定軸轉(zhuǎn)動（不是圓柱體中心軸），類似于公轉(zhuǎn)之類的？我嘗試了很多辦法，包括在邊界條件上施加UR3 或VR3的值，貌似物體都是固定不動的，或者是沿著直線走的，根本就不是預(yù)想的那樣（公轉(zhuǎn)）只加UR3肯定是不行的?？梢栽囋囅旅娴姆椒ǎ?/p>

1、取圓柱上一個節(jié)點（最好是端面圓心位置），定義其U1、U2位移，通過控制U1、U2使圓柱體到軸的距離不變，這樣就有了在X-Y平面內(nèi)的剛體位移，要沿圓弧運動就看你怎么控制U1、U2了。

2、要不要自轉(zhuǎn)，隨便定義UR3即可。

注意：圓周運動的軌跡實際是位移，不是自轉(zhuǎn)。另外，提供一個簡單方法。

在軸上定義一個rp，rp與圓柱端面圓心建立MPC或其他約束，約束rp的位移自由度，給rp加UR3位移，這樣圓柱也能跟著rp繞rp所在的軸轉(zhuǎn)動。

運動問題從幾何角度考慮一下，很簡單。

簡化方法：我記得以前學(xué)數(shù)學(xué)的時候，圓周可以看成是無限多邊形（現(xiàn)在很多軟件其實也是這么實現(xiàn)的，比如CAD的圓就是多邊形簡化的）

既然這樣，可以把繞軸公轉(zhuǎn) 圓弧軌跡 簡化成 多邊形軌跡

也就是

描點了 設(shè)置多個分析步 多個U1 U2（U1 和 U2取圓周上的點）

方法二: ABAQUS 里面有connecter 其中有一種連接方式 是hinge 可以實現(xiàn)繞軸旋轉(zhuǎn) 建立兩個參考點 在兩個參考點之間hinge連接

一個點做軸 一個點和圓柱coupling 在做軸的點上加UR3 即可實現(xiàn)公轉(zhuǎn)

最近碰到一個計算時間比較長的Explicit模型，分析步4個，每個分析步計算時間都超過8小時，由于用筆記本計算，下班就得中斷計算關(guān)機走人，這樣就不得不考慮重啟動的問題。

經(jīng)過這幾天的摸索，發(fā)現(xiàn)對于這種計算時間長的Explicit模型，要想是實現(xiàn)意外中斷以后繼續(xù)計算，最好的方法就是用recover。

方法其實很簡單：

要想在中斷以后實現(xiàn)無縫連接（接著上一個increment繼續(xù)計算），必須在計算之前做一個Restart的設(shè)定，CAE：step模塊，output->restart requests->intervals，默認(rèn)為1，將這個數(shù)據(jù)改大一點，我改為100，與field output的inteval對應(yīng)起來。有了這個前提，不管什么時候中斷，只需把該job type改為recover，然后提交就會接著中斷的部分繼續(xù)算。再中斷再提交就行了，不需要做任何改動，直至計算完成。Command命令：abaqus job=jobname recover int。

解釋：

restart前提條件的設(shè)定，Explicit默認(rèn)的會寫出各種重啟動需要的文件，只不過默認(rèn)的interval=1，就是說只有在每個分析步的開始和結(jié)尾才有重啟動點，這樣不能實現(xiàn)哪中斷哪重啟。要實現(xiàn)無縫連接，理論上這個數(shù)據(jù)越大越好，只要不超過一個分析步increment總數(shù)。但是大了會增加數(shù)據(jù)寫入量，也沒必要。這個前提設(shè)定和restart一樣。

好處：

只需要在第一次重啟動的時候改job type為recover，提交就ok了，當(dāng)然要記得刪除lck文件。后面中斷了就只需提交就行了，不需做任何改動，方便。而且odb，sta等文件都是在原文件的基礎(chǔ)上續(xù)寫，不存在合并odb文件的問題。用這個方法，甚至可以人為中斷后再接著計算。

一點體會，供參考。

1、首先需要重啟動的case計算時有設(shè)置重啟動輸出，即inp中以下類似代碼關(guān)鍵有write： *Restart, write, number interval=8, time marks=YES

2、查看計算結(jié)果后綴為.sta的文件，觀察最后restart number數(shù)值，如Restart Number 8 at 3.0600；

3、合并原計算結(jié)果文件，jobname為斷電前使用的文件名（需先刪除原文件后綴為.lck的文件）

abq6101 job=jobname convert=all int

4、設(shè)置重啟動讀入值，interval等于最后restart number減去1，例 *RESTART, READ, STEP=4, END STEP, INTERVAL=7

5、設(shè)置重啟動inp文件，新建一個AMPLITUDE文件，時間設(shè)置為計算總時間，*Amplitude, name=Amp-3,TIME=total time

3.0, 0.，3.06，1.6、改變相應(yīng)邊界條件； 斷電前使用

*Boundary, amplitude=Amp-2, type=VELOCITY 改為

*Boundary, amplitude=Amp-3, type=VELOCITY

7、分析計算，restartname為新建的重啟動inp文件名稱，jobname為斷電前使用的文件名 abq6101 job=restartname oldjob=jobname cpus=12 int mp_mode=threads double

注意：

默認(rèn)使用的AMPLITUDE中time=step time，跟AMPLITUDE相關(guān)的邊界條件都需要改變；

對時間增量步理解

abaqus的step里有maximum number of increment、initial increment、minimum increment、maximum increment四個量許多網(wǎng)友不知怎樣設(shè)置合理，合理設(shè)置是建立在深刻理解基礎(chǔ)上的。

要理解這個問題，首先需要了解abaqus的計算過程和有限元計算收斂性問題，abaqus首先用initial值輸入進(jìn)行疊代計算，如果計算結(jié)果收斂，則繼續(xù)以這個值代入計算下一步，如果不收斂，則自動減小時間步長（time increment）重新計算直到收斂然后計算下一步。

但是如果時間步長減小到最小值minimum時計算結(jié)果還是不收斂，則abaqus將停止計算，由此可知maximum值和minimum值分別是abaqus在收斂計算時時間步長的上下限，同時total time=求和(time increment*number)，當(dāng)時間步長很小時，需要計算的步數(shù)number相應(yīng)增大(電腦計算花的時間也隨之增大)，因此number一般要設(shè)置較大值。minimum并不是越小越好，因為1)number即計算時間增大2)abaqus計算精度約在10^(-5),當(dāng)時間步長小于這個值，計算結(jié)果已經(jīng)沒什么意義了。

有限元計算收斂性與（最小空間步長/時間步長）值有關(guān)，若minimum設(shè)為10^(-5),還是不收斂，可適當(dāng)減小空間步長（即把網(wǎng)格畫細(xì)點），當(dāng)然還有一些其他辦法，如果實在計算不了，也許是模型本身有點問題，或改為顯示explicit計算

總而言之，maximun number要適當(dāng)設(shè)置較大值，initial可適當(dāng)改?。ㄈ?2，-3量級），minimum（-5量級）不要修改，maximum值影響不大，可不改

顯示分析的求解時間與以下幾個方面有關(guān)： 材料密度。加快求解有一種方法就是增大密度，密度增加100倍，時間降低為原來的十分之一。當(dāng)然，這種情況下必須是密度這一屬性對你要求的結(jié)果影響不大的前提下。最小的網(wǎng)格尺寸.這決定計算機每步求解的最大步長。網(wǎng)格尺寸越小，最大時間步長越大。即使只有一個網(wǎng)格尺寸很小，其他都很大，最大時間步長也提不高。所以，避免產(chǎn)生個別小的網(wǎng)格尺寸是十分必要的。網(wǎng)格數(shù)量。網(wǎng)格數(shù)量越大，每求解一步所需的時間越長。因此，可以適當(dāng)優(yōu)化網(wǎng)格，采用局部加密的方法，減少網(wǎng)格數(shù)量。4 模型復(fù)雜程度，主要是指模型是否有接觸、大變形、多物理場耦合等。越復(fù)雜，時間越長。

以上，希望對你有用。

發(fā)現(xiàn)當(dāng)轉(zhuǎn)動模型時模型就會變成靜態(tài)線框下的透明體，如何才能把他改成轉(zhuǎn)動時不是透明的呢？