第一篇:ABAQUS單元小結
ABAQUS單元小結
1、單元表征
單元族:單元名字里開始的字母標志著這種單元屬于哪一個單元族。C3D8I是實體單元; S4R是殼單元; CINPE4是無限元; 梁單元; 剛體單元; 膜單元;
特殊目的單元,例如彈簧,粘壺和質量; 桁架單元。
自由度dof(和單元族直接相關):每一節點處的平動和轉動 1 1方向的平動 2 2方向的平動 3 3方向的平動 4 繞1軸的轉動 5 繞2軸的轉動 6 繞3軸的轉動 開口截面梁單元的翹曲 8 聲壓或孔隙壓力 9 電勢 度(或物質擴散分析中歸一化濃度)12+梁和殼厚度上其它點的溫度 軸對稱單元 r方向的平動 2 z方向的平動 6 r-z方向的轉動
節點數:決定單元插值的階數
數學描述:定義單元行為的數學理論
積分:應用數值方法在每一單元的體積上對不同的變量進行積分。大部分單元采用高斯積分方法計算單元內每一高斯點處的材料響應。單元末尾用字母“R”識別減縮積分單元,否則是全積分單元。
ABAQUS擁有廣泛適用于結構應用的龐大單元庫。單元類型的選擇對模擬計算的精度和效率有重大的影響;
節點的有效自由度依賴于此節點所在的單元類型;
單元的名字完整地標明了單元族、單元的數學描述、節點數及積分類型; 所用的單元都必須指定單元性質選項。單元性質選項不僅用來提供定義單元幾何形狀的附加數據,而且用來識別相關的材料性質定義;
對于實體單元,ABAQUS參考整體笛卡爾坐標系來定義單元的輸出變量,如應力和應變。可以用*ORIENTATION選項將整體坐標系改為局部坐標系;
對于三維殼單元,ABAQUS參考建立在殼表面上的一個坐標系來定義單元的輸出變量。可以用*ORIENTATION選項更改這個參考坐標系。2.實體單元(C)實體單元可在其任何表面與其他單元連接起來。C3D:三維單元
CAX:無扭曲軸對稱單元,模擬3600的環,用于分析受軸對稱載荷作用,具有軸對稱幾何形狀的結構;
CPE:平面應變單元,假定離面應變ε33為零,用力模擬厚結構; CPS:平面應力單元,假定離面應力σ33為零,用力模擬薄結構; 廣義平面應變單元包括附加的推廣:離面應變可以隨著模型平面內的位置線性變化。這種數學描述特別適合于厚截面的熱應力分析。
可以扭曲的軸對稱單元:用來模擬初始時為軸對稱的幾何形狀,且能沿對稱軸發生扭曲。這些單元對于模擬圓柱形結構,例如軸對稱橡膠套管的扭轉很有用。
反對稱單元的軸對稱單元:用來模擬初始為軸對稱幾何形狀的反對稱變形。適合于模擬像承受剪切載荷作用的軸對稱橡膠支座一類的問題。
如果不需要模擬非常大的應變或進行一個復雜的,改變接觸條件的問題,則應采用二次減縮積分單元(CAX8R,CPE8R,CPS8R,C3D20R)
如果存在應力集中,則應在局部采用二次完全積分單元(CAX8,CPE8,CPS8,C3D20等)。
對含有非常大的網格扭曲模擬(大應變分析),采用細網格劃分的線性減縮積分單元(CAX4R,CPE4R,CPS4R,C3D8R等)
對接觸問題采用線性減縮積分單元或非協調元(CAX4I,CPE4I,CPS4I, C3D8I)的細網格劃分。
如果在模型中采用非協調元應使網格扭曲減至最小。三維情況應盡可能采用塊狀單元(六面體)。當幾何形狀復雜時,完全采用塊體單元構造網格會很困難,因此可能有必要采用稧形和四面體單元,但盡量少用,并遠離需要精確求解的區域。
一些前處理程序包括網格劃分方法,它們可用四面體單元構造任意形狀的網格。只要采用二次四面體單元(C3D10),其結果對小位移問題應該是合理的。
小結:
在實體單元中所用的數學公式和積分階數對分析的精度和花費有顯著的影響;
使用完全積分單元,尤其是一階(線性)單元,容易形成自鎖現象,正常情況不用;
一階減縮積分單元容易出現沙漏現象;充分的單元細化可減小這種問題; 在分析中如有彎曲位移,且采用一階減縮積分單元時,應在厚度方向至少用4個單元;
沙漏現象在二階減縮積分單元中較少見,一般問題應考慮應用這些單元; 非協調單元的精度依賴于單元扭曲的量值;
結果的數值精度依賴于所用的網格,應進行網格細化研究以確保該網格對問題提供了唯一的解答。但是應記住使用一個收斂網格不能保證計算結果與問題的實際行為相匹配:它還依賴于模型其他方面的近似化和理想化程度;
通常只在想要得到精確結果的區域細劃網格; ABAQUS具有一些先進特點如子模型,它可以幫助對復雜模擬得到有用的結果。
3.殼單元(S)可以模擬有一維尺寸(厚度)遠小于另外兩維尺寸,且垂直于厚度方向的應力可以忽略的結構。
一般殼單元:S4R,S3R,SAX1,SAX2,SAX2T。對于薄殼和厚殼問題的應用均有效,且考慮了有限薄膜應變;
薄殼單元:STRI3,STRI35,STRI65,S4R5,S8R5,S9R5,SAXA。強化了基爾霍夫條件,即:垂直于殼中截面的平面保持垂直于中截面;
厚殼單元:S8R,S8RT。二階四邊形單元,在小應變和載荷使計算結果沿殼的跨度方向上平緩變化的情況下,比普通單元產生的結果更精確;
對于給定的應用,判斷是屬于薄殼還是厚殼問題,一般:如果單一材料制造的各向同性殼體的厚度和跨度之比在1/20-1/10之間,認為是厚殼問題;如果比值小于1/30,則認為是薄殼問題;若介于1/30-1/20之間,則不能明確劃分。由于橫向剪切柔度在復合材料層合殼結構中作用顯著,故比值(厚跨比)將遠小于“薄”殼理論中采用的比值。具有高柔韌中間層的復合材料(“三明治”復合材料)有很低的橫向剪切剛度并且幾乎總是被用來模擬“厚”殼;
橫向剪切力和剪切應變存在于普通殼單元和厚殼單元中。對于三維單元,提供了可估計的橫向剪切應力。計算這些應力時忽略了彎曲和扭轉變形的耦合作用,并假定材料性質和彎曲力矩的空間梯度很小;
殼單元可以使用每個單元的局部材料方向,各項異型材料的數據,如纖維增強復合材料,以及單元輸出變量,如應力和應變,都按局部材料方向而定義。在大位移分析中,殼單元上的局部材料軸隨著材料各積分點上的平均運動而轉動;
線性、有限薄膜應變、四邊形殼單元(S4R)是較完備的而且適合于普通范圍的應用;
線性、有限薄膜應變、三角形殼單元(S3R)可作為通用的殼單元來應用。由于在單元內部近似為應變場,精細的網格劃分可用于求解彎曲變形和高應變梯度;
考慮到在復合材料層合殼模型中剪切柔度的影響,將采用“厚”殼單元(S4R,S3R,S8R)
四邊形或三角形的二次殼單元,用于一般的小變形薄殼是很有效的。它們對剪力自鎖和薄膜鎖死是不敏感的;
在接觸模擬中不用選用二階三角形殼單元(STRI65),要采用9節點的四邊形殼單元(S9R5);對于僅經歷幾何線性行為的非常大的模型,線性、薄殼單元(S4R5)一般將比通用殼單元花費更少;
小結:
殼單元的橫截面特性可以由沿厚度方向的數值積分確定(*SHELL SECTION),或在分析開始時應用計算的橫截面剛度(*SHELL GENERAL SECTION);
*SHELL GENERAL SECTION是非常有效的,但僅用于線性材料,*SHELL SECTION可用于線性和非線性材料;
數值積分在沿殼厚度方向的一系列積分點上進行。這些積分點就是單元變量可以被輸出的位置。最外層的積分點位于殼單元的表面。
殼單元法線方向決定了單元的正和負表面,為了正確地定義接觸和解釋輸出數據,必須知道其對應的是哪個面。殼法線還定義了施加在單元上正壓力載荷的方向,并可以在ABAQUS/Post中畫出;
殼單元利用材料方向局部化到每個單元。在大位移分析中,局部材料軸隨單元而轉動。*ORIENTATION被用來定義非默認的局部坐標系統。單元的變量,如應力和應變,在局部方向輸出;
*TRANSFORM定義節點的局部坐標系,集中載荷和邊界條件被應用在局部坐標系中。所用節點的輸出,如位移,也默認為基于局部的坐標系;
矢量圖可以使模擬結果可視化,特別是用來觀察結構的運動和載荷路徑。
4.梁單元(B)
模擬一維尺寸(長度)遠大于另外二維尺寸的構件,且只有長度方向上的應力比較顯著。
對于包含接觸的任何模擬,應使用一階、剪切變形的梁單元(B21,B31)如果結構剛度非常大或者非常柔軟,在幾何非線性模擬中應當使用雜交梁單元(B21H,B32H,等)
使用歐拉-伯努利(三次)梁單元(B23,B33)精度很高,可模擬承受分布載荷作用的梁,例如動態振動分析。如果橫向剪切變形也很重要,要使用鐵摩辛柯(二次型)梁單元(B22,B32)
模擬有開口薄壁橫截面的結構,應當使用考慮了開口截面翹曲理論的梁單元(B31OS,B32OS)
小結:
梁單元的性質由截面(*BEAM SECTION或*BEAM GENERAL SECTION)的數值積分決定,或直接給出截面積、慣性矩和扭轉常數(*BEAM GENERAL SECTION);
當使用*BEAM GENERAL SECTION選項時,模擬開始時進行一次數值積分,并且假定材料是彈性的;
ABAQUS包括大量的標準橫截面形狀。其它形狀可以通過定義SECTION=ARBITRARY來模擬;
必須定義橫截面取向,方法是通過給出第三個節點,或者在單元性質定義中定義一個矢量。截面取向在ABAQUS后處理中可以畫出;
當梁作為殼的加強構件使用時,梁的橫截面可能偏離節點;
線性和二次型包含剪切變形的影響,三次型梁不考慮剪切柔度。開口截面梁準確地模擬了扭轉和薄壁開口截面翹曲(包括翹曲約束)的影響;
多點約束和約束方程可以用來連接模型中鉸接、剛性連接等節點的自由度;
“彎矩”型圖使得像梁這樣的一維單元的結果很清楚地表示出來; ABAQUS后處理圖的硬拷貝可以得到PostScript和HPGL的格式。
5.桁架單元(T)
只能承受拉伸和壓縮載荷的桿,不能承受彎曲,模擬鉸接框架結構,近似模擬線纜和彈簧。
6.剛體單元(R)
沒有獨立的自由度。7.非線性分析
小結:
結構問題中存在著三種非線性來源:材料、幾何和邊界(接觸)。這些因素的任意組合都可以出現在ABAQUS的分析中;
幾何非線性發生在位移量值影響結構響應的情況下。這包括大位移和轉動效應、突然翻轉和載荷硬化;
非線性問題是利用牛頓-拉弗森方法來進行迭代求解的。非線性問題比線性問題所需要的計算機資源要高許多倍;
非線性分析步被分為許多增量步。ABAQUS通過迭代,在新的載荷增量結束時近似地達到靜力學平衡。ABAQUS在整個模擬計算中完全控制載荷的增量和收斂性;
狀態文件允許在分析運行時監控分析過程的進展。信息文件包含了載荷增量和迭代過程的詳細信息;
在每個增量步結束時可以保存計算結果,這樣結構響應的演化就可以用ABAQUS/Post顯示出來。計算結果也可以用x-y圖的形式繪出。
8.材料
小結:
ABAQUS包含一個廣泛的材料庫,可模擬各種工程材料的性質。其中包括金屬塑性和橡膠彈性模型;
金屬塑性模型的應力-應變數據必須用真實應變定義; 金屬塑性模型假定材料具有一旦屈服即不可壓縮的性質。這將對應用于彈-塑性模擬的單元類型帶來某些限制;
多項式和奧根應變能函數可應用于橡膠材料的彈性(超彈性)。兩種模型均允許直接用實驗數據來確定材料的系數。實驗數據必須是名義應力和名義應變的值;
在超彈性材料模型中的穩定性警告,說明所要分析的應變范圍不合適; 存在對稱性時,可以只考慮部分模型從而減小模擬的尺寸。可通過施加適當的邊界條件來反映結構其余部分的效應;
大畸變問題的網格設計比小位移問題更加困難。在分析的任何階段,網格中的單元務必不能過于畸變;
ABAQUS/Post中的*DEFINE CURVE命令允許處理曲線上的數據以生成新的曲線。兩條曲線或一條曲線與一個常數可以加、減、乘、除。曲線還可以求導、積分和合并。
9.動態問題
具有下列特征的問題適于采用線性瞬態動力學分析:
系統是線性的:線性材料行為,無接觸條件,無非線性的幾何效應; 響應只受較少的頻率支配。當響應中各頻率成分增加時,例如撞擊和沖擊情況,振型疊加方法的效果將大大降低;
載荷的主要頻率在可得到的固有頻率范圍內,以確保對載荷的描述足夠精確;
由于任何突然加載所產生的初始加速度能用特征模型精確描述; 對系統的阻尼不能過大。小結:
動態分析包括結構的慣性效應;
*FREQUENCY可以計算結構的固有頻率和振型;
通過振型疊加,可以確定線性系統的動態響應。這一方法盡管有效,但是不能用于非線性問題;
線性動態過程可以計算瞬態載荷的瞬態響應、諧振動下的穩態響應、支座移動造成的響應峰值和隨機載荷的響應;
為了準確表示結構的動態行為,必須選擇足夠多的振型。總的等效模型質量應占可動質量的90%以上;
用戶可以給定直接模態阻尼、瑞利阻尼和復合模態阻尼。但是由于固有頻率和振型的計算都是基于無阻尼的結構,所以此法只適用于低阻尼結構;
模態技術不適用于非線性的動態模擬。在這種情況下必須采用自己的時間積分方法(*DYNAMIC)
*AMPLITUDE選項可以描述隨時間任意變化的載荷,以及給定的邊界條件; 振型和瞬態結果可以在ABAQUS/Post中用動畫顯示。這對于理解動態響應和非線性靜態分析十分有幫助。
10.多步驟分析
小結:
一個ABAQUS模擬過程可以包含任意數目的步驟;
一個分析步驟就是一段“時間”,在這段時間里ABAQUS計算模型對一套指定載荷和邊界條件的響應。這一步驟中所用的特殊分析過程確定了這個響應的特征;
在一個一般分析步驟中,結構的響應可能是線性的,也可能是非線性的; 每一個一般步驟的開始狀態是上一個一般步驟的結束狀態。這樣,在一個模擬中模型的響應隨一系列一般步驟而演化;
線性擾動步驟計算結構對擾動載荷的線性響應。這個響應的基本狀態是相對于最后一個一般步驟結束時模型的狀態所定義的;
在一般步驟中任何載荷選項里的OP參數(例如*BOUNDARY,*CLOAD和*DLOAD中)控制著這些選項中所指定的數值是如何與前面步驟中定義的數值相互作用的;
只要存儲了一個重新啟動文件就可以進行重新啟動分析。重新啟動文件可以用來繼續一個中斷的分析或者給模擬添加附加的載荷過程。11.接觸
小結:
接觸分析需要一個謹慎的邏輯方法。如果必要,將分析分解成幾步執行,并緩慢地施加荷載,以保證很好地建立接觸條件;
一般地,對分析的每一步最好采用分離步驟進行,即使只是因為載荷而改變邊界條件。您幾乎肯定要比預期情況應用更多的步驟,但模型則收斂得更容易。如果想一步就將所有的載荷加上,接觸分析是難以完成的;
在對結構施加工作載荷之前,要在所有部件之間取得穩定的接觸條件。如果必要,采用臨時的邊界條件,在以后階段再消除這些約束。只要所提供的約束不產生永久的變形,對最后的結果應該毫無影響;
不用對接觸面上的節點施加邊界條件,即在接觸方向上限制節點。如果有摩擦,不要在任何自由度上約束這些節點:可能導致零主元信息;
對于接觸模擬,總要試圖使用一階單元。
第二篇:Abaqus遇到的問題小結
Abaqus遇到的問題小結
1.Abaqus量綱系統(SI):m、N、Kg、s、Pa、J、Kg/m。
2.有限元單元類型(Element Type):實體單元、殼單元、桿件單元(梁Beam、桁架Truss)3.Abaqus怎樣將部件的不同部分設定為不同的子集?
3tools >partition cell/plane(如果是要分層使用,則進行分割,分別賦予不同的屬性;或者對其中的單元進行設置Set;如果要在運算中替換則需要在關鍵詞里添加語句定義)
4.建模的過程中輸入的尺寸被自動降低精度(四舍五入),請問如何設置尺寸精度(小數點位數)?
在 sketch option里打開sketch模塊(左下角最下面)打開sketchoption對話框,其中decimal places 表示小數點位數,默認為2,最高能調到6。
5.abaqus中出現,Dependent partinstancesconnot be edited,怎么辦?
mesh on part , mesh on instance一個非獨立實體只是原始部件的一個指針,可以對原始部件劃分網格,但是不能對一個非獨立實體劃分網格,即mesh on part。
方法:左邊模型樹里面,Assembly->instances里面的子選項右擊,點Make independent;或者點擊model旁邊的Part,對part進行網格劃分。
6.材料屬性中拉伸強度和壓縮強度怎么定義?
計算結果有每個計算節點的應力,計算中用不到強度值(彈性力學三組基本方程不涉及強度,求解過程用不到)。(強度值是用來校核的:斷裂,屈服或超過線性變形)
7.如何選擇創建獨立實體還是非獨立實體?
如果集合中包含許多具有相同性質的部件,則創建非獨立實體(Dependent)。反之,如果集合中包含許多不同性質的部件,則創建獨立實體(Independent)要有優勢些。
8.地基土的分層問題?
方法1:定義為一個part,然后partition分層賦材料屬性,不需要設置接觸。(超級贊!!)方法2:使用多個part,每個part單獨定義屬性。一般不用設置一些tie,各層間只有彈性模量相差較大時再在各層間定義tie。
9.tie捆綁約束。
tie功能很強大,可以模擬實體與實體,實體與殼體,殼體和殼體的連接,而且允許連接對象之間有距離。
在tie面板中,Position Tolerance位置容差是一個重要的參數。一般選擇默認設置usecomputed default,abaqus會根據你要tie的對象的網格尺寸,自動選取個容差值,一般為網格基本尺寸5%。于是,連接對象之間,距離小于這5%容差值的節點便施加tie約束,也就是對應的兩節點的自由度擁有一樣的數值。
如果連接對象之間有距離怎么辦?若還是用usecomputed default,則是無法起到tie約束效果的。這就要手動設置specify distance為一個適合的參數。這個值要大于連接對象對應節點之間的距離。若不是很確定,可以先設一個較大的值,計算完后在Visualization調出Node sets,觀察tie實際起作用的節點區域。再對specify distance值進行調整。
在使用specify distance,需將Adjust slave surface initial position選項取消。否則abaqus會將兩連接對象拉近,網格產生畸變,產生錯誤,中止計算!
10.abaqus怎么將模型分割畫網格?
劃分模型的時候有全局劃分網格按鈕,旁邊有個局部劃分網格按鈕,(經常要用到partition技術),點它給需要另外劃分的部分部種子,然后下面有按區域劃分網格。全局劃分網格和局部劃分網格分別是mesh part和mesh region。
網格劃分完畢之后,最好在verify中對其中質量不好的單元進行編輯處理,避免在job模塊里面出現錯誤或警告信息。
11.邊界條件設置在初始步還是后續的分析步中?
設置在initial。
12.interaction面選擇不上?
在此模塊里面,一般先定義相應的surfaces,并分別賦予有意義的容易識別的名字,如slavesurface、master surface,并選擇接觸面的正確方向,如果方向選擇錯誤的話,經常在job模塊中會出現無法收斂的錯誤信息。
13.通過ABAQUS計算得出的應力值不連續是什么原因造成的?
這主要是因為ABAQUS的應力和應變等單元變量值(相對位移等節點變量)的計算方式所引起的 這些單元變量值主要是通過單元內的高斯積分點計算得出,然后再通過外插值法計算出節點上的應力應變值,因此對于兩個相鄰單元共同占用的節點,從兩個單元外插得出的值有時是不一樣的。對于大多情況.如果網格足夠密,高斯積分點上的值不會相差很大.所以可得到一個連續性比較好的結果另外可以通過設定平均值的臨界值來決定contour的連續性(即legend下默認的75%)。比如說,如果同一節點扶不同單元得出的外插值相差75%以下(即較大的外插值減去較小的外插值的差乘以大值),那么contour輸出這兩個值的平均值。
同時需要注意的是,這個默認的75% 可以修改,如果設為0%.那么整個contour的單元間幾乎都不連續,如果設在100%,那么整個圖肯定是連續的。此外,還可以通過設定legend下的limit方式(即外插值或者是平均值)來控制最后的結果。
14.平衡地應力是出現An initial condition has been specified on element 0 but thiselement has not been defined?
是因為修改inp文件時,文件中有未刪除的空行。
15.土體初始地應力平衡方法。
地應力平衡原理:平衡原理和屈服條件。目的:為了使數值模擬獲得一個存在初始應力,而無初始應變的狀態。
平衡地應力方法:
1.Autobanlance 2.*Initial Conditions,type=stress,Geostatic;(關鍵字定義初始地應力,只適合土體表面水平的土體,該方法需給出不同材料區域的最高和最低點的自重應力及其相應坐標)
3.*Initial Conditions,type=stress,input=xxx.csv;(初始地應力提取法,此方法最為通用,可用于不同材料、不同地形,適應性強。)
方法3的具體步驟: i.ii.iii.iv.v.在Abaqus中建立完整的模型,定義部件、截面屬性、裝配件、網格、邊界條件等。單元類型使用的默認的線性縮減積分單元C3D8R。
在Step功能模塊中創建分析步,將類型設為Geostatic。在Load功能模塊中,為整個土體定義重力載荷(類型為Geostatic),并在Component3中輸入-9.8。
在命令行輸入:mdb.models['Model-1'].setValues(noPartsInputFile=ON);在Job模塊中創建名為Job-NoInitialCondition的分析作業,提交分析。
將分析得到的應力場保存為一個文本文件。具體的方法為:打開分析得到的ODB文件Job-NoInitialCondition.odb,選擇菜單Report---FieldOutput,選中Variable對話框中積分點上的各個應力分量(S11、S22、S33、S12、S13、S23);在Setup標簽頁,創建inp輸出文件,并命名為XX.inp,不選Append to file和Column total 和column min/max。
用Excel修改輸出的XX.Inp文件,刪除inp文件中的所有敘述性文字和空格(例如若土體有N層,需找出N-1處文字性敘述的地方并刪除),只保留數字結果,將數字間在順序不變的前提下用逗號隔開,保存成csv以逗號隔開的文件形式XXX.csv。vi.vii.viii.為模型定義初始應力場。因在Abaqus中無法直接定義初始應力,只能手工添加關鍵詞,具體做法:選中菜單Model---Edit Keywords,在*Boundary語句之后,*Step語句之前添加語句*Initial Conditions,type=stress,input=XXX.csv 在Job功能模塊中定義Job-WithInitialCondition作業,提交分析,即得到所需要的結果。
16.Part模塊建立剛體時一定要指定參考點referencepoint,在后面分析過程中的位移、載荷的施加都在施加在此參考點上的。剛體只有平動和轉動,不會有變形,因此參考點的選擇即可以在建立的rigid part上,也可以建立在其他你覺得方便的位置。
17.load模塊,是施加載荷和位移邊界條件的,先對要施加載荷和邊界條件的面、線、節點等建立set,再對相應sets分別設置。在施加載荷的時候,注意載荷的方向性,另外此模塊也可以施加自定義載荷、predifined field等。
18.abaqus結果顯示時數據太小看不清楚怎么辦?
Viewport-->Viewport Annotation Options..在出來的對話框中選擇“Legend”標簽,然后點擊“Set Front...” 之后選擇字體大小,確認。
BY
HIT-WGJ
第三篇:ABAQUS與FORTRAN論壇問題小結
ABAQUS與FORTRAN論壇問題小結
Q:索了一下論壇發現以前的問題都是不同類型的子程序,如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上面是區分不同umat的方法,區分uel應該類似吧
Q:行abaqus的時候出現這樣的提示“Problem during compilation-df.exe not found in PATH.”請問這是什么原因.A: 就是說找不到用于編譯的df.exe文件(好像是Fortran的編譯運行文件),在Windows控制面板的“系統”里,打開高級屬性,將系統的PATH添加上Fortran的路徑,最好也一并把INCLUDE和LIB也都添加上對應的FORTRAN的路徑。
Q: 如何在ABAQUS中調用用戶子程序UMAT
A:abaqus job=*** user=*.for or *.objA: 裝好fortran編譯器的前提下,有兩種方式:
1、命令行:abaqus job=*.inp user=*.for或 abaqus job=*.inp user=*.obj(可以由*.f90編譯后得到)
2、菜單格式:先運行CAE,在general設置中可以指定Umat用戶子程序A: job=jobname user=yourfortran.for(其它類型子程序)A:如果是本構的二次開發:*material, user
Q: abaqus job=job-name user={source-file | object-file}source-file和object-file各是什么? A: fortran的源文件(.for)或者目標文件(.o)在command下面輸入以上命令行就可以正常運行.用戶子程序再調用別的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我也學了沒多久,按照這個過程成功了一個!你們試一下吧!Q: 想用FORTRAN實現:...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窗口里手工輸入文件的絕對路徑則計算正常。但我的目的是實現Fortran程序對ABAQUS的計算程序的自動調用。
A: 說個投機的花招,有沒有試過在前面加上result = SYSTEMQQ('path=e:/abaqu***amples')這樣的命令~~方法簡單,不過能解決問題最好
A:把要算的工況都集中在一個文件夾下,fortran生成.exe文件,放在該文件夾下,然后雙擊運行。計算結束之后再把各種工況分類。笨招,但好用。
A: 我想,或者把systemQQ改一下,先cd e:...,然后再systemQQ(abaqus job=jobname interactive)。
Q: 偶想在FORTRAN里設定一個計時器,其功能就是每隔一段時間(如10s)就去讀一次ABAQUS的(.fil)文件。哪個函數可以實現這個功能?怎么實現?? A: 通過一個while循環和下面的語句:call sleepqq(時間間隔)就可以了。
Q: 本人嘗試使用ABAQUS的子程序SUBROUTINE,根據幫助文件需要編譯通過生成*.for和*.lib文件,可是編譯的時候出現錯誤:Error: Cannot open include file 'ABA_PARAM.INC'不知道是什么原因,怎么解決?
A: 1.確認Visual Fortran及VC++正確安裝,并在安裝ABAQUS之前安裝。2.首先將用戶子程序*.f文件后綴名改為*.for。3.到控制面板-》系統-》高級-》環境變量中的 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函數作批處理遇到問題,在fortran中使用result = systemqq(abaqus job=job1 interactive)作批處理,總是在計算兩個job之后就不繼續往下算了,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如何調用ABAQUS?
A: use the function “system” or “systemqq” USE DFLIB LOGICAL(4)result result = SYSTEMQQ('abaqus job=xx')QA: 首先我用ABAQUS的VARIFY命令對系統進行修正和檢測,這個命令會自動運行,對ABAQUS進行調試等操作,其中涉及到對用戶子程序的模擬計算,如果計算失敗,就會在VARIFY.LOG文件中說明,確實發現“USER SUBROUTINE...和SINGLE...和DOUBLE...”失敗,說明用戶子程序安裝不成功,肯定與VISUAL FORTRAN6.5有關!然后,找到FORTRAN_MAKE文本文件,發現錯誤,可能是與FORTRAN安裝有關!于是卸載FORTRAN,重新安裝在C盤!安裝完畢后,改變環境變量如下: 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我們只需要關注MICROSOFT VISUAL STUDIO文件夾!最后,再次執行VARIFY命令,等運行完畢后,檢查VARIFY文本文件,發現: 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://www.tmdps.cn/bbs/redirect.php?fid=10&tid=9031&goto=nextnewset
第四篇:abaqus在基坑工程中使用小結
abaqus在基坑工程中使用小結
本人在使用abaqus解決基坑支護問題的過程中,得到了版上aba_aba,以及ustc的life大俠等人的幫助,寫些自己使用abaqus的一些感受吧,(其實大多數東西都是從版上得到的相關參考)希望對其他新手能有所幫助。**首先是建模方式
對于初學者比較好掌握的是CAE界面操作方式。但是由于ABAQUS的CAE模塊和輸入文件兩種建模方式是由兩家不同的公司研制的,CAE模塊推出比較晚,所以功能還不是很完全,對于CAE和INP交互操作需要注意幾點:
1、一些功能只能通過編輯INP輸入文件來實現,比如模擬開挖建造時需要用到殺死或者激活單元的語句model change,此時需要直接在CAE下編輯關鍵字,或者生成INP文件后通過記事本等工具打開INP文件后加入關鍵字。
2、CAE對于INP輸入文件中的一些關鍵字不支持,如果在CAE下編輯,則系統會給出不支持關鍵字的提示。對于這樣的情況則只能通過文本編輯工具打開INP來操作。
3、對INP輸入文件進行編輯的方式能支持所有關鍵字功能,因此應用更加廣泛。在許多情況下都需要對INP文件進行操作。**然后是初始應力場的施加。
土工有限元分析中初始應力場的施加是個關鍵。在ABAQUS中,施加初始應力場很方便,可以使用以下語句來完成。
initial conditions,type=stress,geostatic 可以有兩種實現的辦法:
1、如果地表水平,且初始應力就是自重應力時,可以在數據行中按照如下順序填寫數據 施加初始應力的集合,應力值1,應力值1對應的坐標值1,應力值2,應力值1對應的坐標值2,x方向側壓力系數,y方向側壓力系數(在二維問題中可以省略)
這樣的順序定義初始應力實際上默認了初始應力就是自重應力,只需要給定集合中2點的應力值及其對應的坐標值,其他部位的應力值就可以通過1,2兩點的連線采用差值來得到,而水平向的地應力則乘一個側壓力系數得到。(這些內容在abaqus手冊中都有)
2、采用input文件的方式,即先通過計算一個重力荷載作用下沒有初始應力的模型而得到各點的應力值,然后將這些值導入到新的模型中,在*initial conditions,type=stress,后面加入input=XXX.XXX,讓模型接受存有初始應力值的input文件而施加。這樣的方法通用性非常強,但是比1稍微復雜些,具體可以見相關手冊。
另外,對于施加初始應力,對于二維的情況,系統默認地應力方向為Y方向,對于三維情況,默認為Z方向,這一點在初學者中容易弄錯。**模擬開挖建造
ABAQUS中模擬開挖建造通過關鍵字model change,type=element來實現,如果是殺死單元則在后面加上remove,如果激活單元則加上add。
事先需要將需要殺死或激活的土體單元定義為一個集合set。整個應用格式如下: *Model change, type=element, remove Set name1,set name2
其意義如下:首先Model change說明是模型發生變化,后面是種類,因為除了殺死激活單元外,還常用到殺死激活接觸對(type=contact),此處需要指明是殺死激活單元導致的模型變化,后面remove表示殺死單元。重啟一行寫出需要殺死激活的單元集合列表。
在模擬建造過程中,通常一個單元的材料屬性會發生變化,比如土釘支護中面層處單元,在沒有開挖時是土體材料,當開挖建造之后就變成了面層的屬性了。在ABAQUS中模擬材料屬性隨開挖步進行可以使用場量*Field關鍵字來實現。具體過程如下:首先是材料定義上要指定材料性質隨著場變量而變化,如下面語句中,定義彈性時,利用dependencies=1表示材料根據1個場變量變化而變化,在接下來的屬性數據行中每行后面的數字1,2表示了該行屬性對應的這個場變量的值,如第一行的彈性模量和泊松比值是場變量為1的值,而第二行的彈性模量和泊松比值是場變量為2的值。*Elastic, dependencies=1
4e+07, 0.35, , 1.1.5e+10, 0.2, , 2.接下來就是在每個分析步中利用*field,variable=1這樣語句指定分析步中使用的場變量數量,本例中值為1個場變量。接下來的數據行則是表示使用場變量的集合名稱及其在該分析步中的場變量值。如果這個值是1,根據在定義材料時的場變量定義,值為1時對應于材料屬性第一行值,這樣就實現了材料屬性隨分析步變化的功能。
**模擬地連墻與土的接觸ABAQUS中對于接觸的模擬可以通過定義接觸對來實現。系統可以允許面對面、點對面的接觸等多種接觸形式。在模擬開挖過程中,開挖土體與墻體的接觸屬性會喪失,根據經驗,最好能將接觸按照開挖步驟劃分為許多對接觸對(所謂的細化),根據開挖的進行,將接觸對先殺死,這樣系統比較容易找到準確的接觸狀態,使接觸模擬順利進行。模擬接觸時,定義接觸對要分清主控面和受控面。主控面和受控面的選擇一般遵循以下原則:
一、材質較軟的為受控面,材質較硬的為主控面。
二、網格較密的為受控面,網格稀疏的為主控面。應用到網格劃分的時候就要先根據材質的相對軟硬程度確定主控面,以便在網格劃分中按照網格疏密關系正確劃分網格。主控面和受控面定義一定要正確,這是影響接觸模擬的關鍵。如果一個節點上同時定義了兩個接觸面,或者同時定義了一個接觸面和邊界條件,系統則會產生過約束(over constraint)的提示,很容易影響到接觸的模擬。
分享<隧道施工應力釋放法的ABAQUS實現【inp及代碼見2#】> 應力釋放法模擬隧道開挖主要有兩種實現途徑,大家可以參考帖子: http://forum.simwe.com/thread-765338-1-1.html
簡單的說,主要就是參數弱化法和施加節點反力法。這兩種方法各有優劣。參數弱化法:
優點:
1、便于實現;
2、便于控制應力釋放過程;
缺點:
1、無法給出明確的應力釋放率;
2、無法用在劍橋模型中,因為找不到一個合適的參數進行弱化。施加節點反力法:
優點:
1、能夠明確給出應力釋放率;
2、物理意義較為明確; 缺點:實現過程不是很方便
下面跟大家匯報一下我的做法,希望有經驗的童鞋給出建議,各種板磚,放心拍吧~
一、參數弱化法
參數弱化方法其實沒有什么難度,版內已經討論了無數次了。各種牛們也是不厭其煩的跟新手們說了一遍又一遍,可惜有些童鞋天生喜歡求助,上來就發帖問,完全置搜索功能于不顧。所以我在這要再次呼吁,搜索處要更加醒目,加紅字,加粗,加浮動,另外加旁注:親,如果你想快速入門aba,搜索吧,這是悶騷技術人新一代的學習方法,包會哦!
參數弱化,主要就是弱化彈性模量,這個可以通過關鍵字*field來實現,具體腫么實現,搜索吧。下面給出計算后的云圖,拋磚引玉。地應力平衡,勉強可以吧:
待開挖區域彈性模量折減20%:
如何添加襯砌以及如何設置追蹤單元我想再單開一貼和大家進行總結探討。上面位移矢量圖隱去了待開挖部分,是為了矢量顯示清晰,實際上這一步待開挖(中心土體)部分實際還是存在的,只不過彈性模量折減過了。從圖中可以看出,這種方法計算出的位移場并不是很符合實際情況。隧道開挖土體位移場可以參見帖子:http://forum.simwe.com/thread-985421-1-1.html
在這里也順便給這個帖子做個個人的解答吧。關于該貼提出的強制位移法,我想在后面專開一貼和大家討論總結。
有經驗的童鞋可否告知你們采用折減法做出的位移場是否也和上面一樣的?
二、施加節點反力法
這一方法的技術難點在于,如何獲得地應力平衡后隧道周邊土體節點反力。通過搜索和查看manual,基本可以鎖定兩種輸出量:NFORC & RF。其中,RF只能在有邊界約束處輸出,NFORC根據manual我猜測應該是由節點所涉及單元的應力進行外插平均后得到的節點力,參見帖子:http://forum.simwe.com/thread-883192-1-1.html10#。不管怎樣,總的思路是,通過對開挖后隧道相鄰土體節點施加非均布節點力,使得此時的模型(土體)保持初始地應力和位移準零狀態,然后施加原節點荷載*(1-應力釋放率),以此來模擬不同的應力釋放率下隧道及土體的力學響應。
以前的帖子對這個問題討論了很多,有代表性的列舉如下(話說論壇的收藏是個好功能):
http://forum.simwe.com/thread-872620-1-1.html 該貼lz給出了四種獲取節點力的方法,都很有參考價值。最后采用了NFORC,但得到的效果不是很理想。我做了一下,發現的確如該lz所說,NFORC輸出的節點力是不能直接施加到原節點上的,即使乘以2倍,也只能是接近實際的節點力。所得到的位移場在初始狀態就有很大的位移,不符合前面提出的總體思路。該貼5#和7#都提出采用約束洞周位移獲取RF。事實證明,這個方法是可行的,但不同的操作方式得到的結果千差萬別。該貼7#提出的思路是對的,但給出的文件中將洞周位移約束貫穿了整個施工過程,這樣造成的后果是,其計算起始點已經不再是green field的情況,得到的應力場自然不符合實際情況,那么在錯誤的應力場上做的計算怎么能對呢?
http://forum.simwe.com/thread-854591-1-1.html
該貼給出了施加洞周約束獲取反力的正確思路,即施加洞周位移約束要在地應力平衡后進行。這樣能保證施加位移約束時的地應力初始狀態。下面就我做的例子詳細說明。
模型:
材料參數:
修正劍橋模型
*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.劍橋模型參數的設置參見: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...",如果你確定有地應力平衡,那么就請你檢查是否是用的直接加重力平衡的吧。多孔介質彈性模型是需要在initial conditions中定義初始地應力的。
地應力平衡:
平衡效果一般,主要是做實例,只平衡了一次,到了10E-4就沒再管了。
施加洞周邊界約束,同時remove掉待開挖部分
可以看出,這一步進行完后,剩余土體的應力基本沒有發生改變,也就是基本保持了初始地應力,在這種情況下導出節點力再倒入折減后支護力才是正確的。可以看出施加洞周位移約束后土體位移分布發生了變化,但數量級還是很小,基本沒變,可以忽略。最后一張圖是洞周節點反力云圖,在這種情況下就可以輸出節點反力了。下面介紹一下我的做法。
第五篇:abaqus自我總結
如何輸入inp文件
inelastic heat fraction:非彈性功(一般是塑性功)轉換為熱量的比例,比例越大,產熱越高,溫升越高
rate dependent:本構率相關參數,表示材料的屈服流動應力不僅和等效塑性應變相關,而且還與應變率有關,一般是正相關
如何將二維圖表的背景變成白色
我的怎么就改變了呢 說了雙擊
舉例子說說我的理解吧:history output->盯住一個節點,這個節點在整個加載時間段內應力變化情況,用xy圖表示就是一條TIME-STRESS曲線;field output->盯住一個時刻,所有節點在這個時刻時的應力值,直觀表示是云圖。當然在xyplot中你可以選擇盯住一個節點按時間輸出該節點的力(即ODB field output),不過最多只能在各個幀對應的時刻輸出,xy圖中點的數量較少。如果我錯了的話請大家指正哈
view---overlayplot--然后每一個時刻創建一個新圖層,然后plotoverlay
怎么樣才能使圓柱體繞另一固定軸轉動(不是圓柱體中心軸),類似于公轉之類的?我嘗試了很多辦法,包括在邊界條件上施加UR3 或VR3的值,貌似物體都是固定不動的,或者是沿著直線走的,根本就不是預想的那樣(公轉)只加UR3肯定是不行的。可以試試下面的方法:
1、取圓柱上一個節點(最好是端面圓心位置),定義其U1、U2位移,通過控制U1、U2使圓柱體到軸的距離不變,這樣就有了在X-Y平面內的剛體位移,要沿圓弧運動就看你怎么控制U1、U2了。
2、要不要自轉,隨便定義UR3即可。
注意:圓周運動的軌跡實際是位移,不是自轉。另外,提供一個簡單方法。
在軸上定義一個rp,rp與圓柱端面圓心建立MPC或其他約束,約束rp的位移自由度,給rp加UR3位移,這樣圓柱也能跟著rp繞rp所在的軸轉動。
運動問題從幾何角度考慮一下,很簡單。
簡化方法:我記得以前學數學的時候,圓周可以看成是無限多邊形(現在很多軟件其實也是這么實現的,比如CAD的圓就是多邊形簡化的)
既然這樣,可以把繞軸公轉 圓弧軌跡 簡化成 多邊形軌跡
也就是
描點了 設置多個分析步 多個U1 U2(U1 和 U2取圓周上的點)
方法二: ABAQUS 里面有connecter 其中有一種連接方式 是hinge 可以實現繞軸旋轉 建立兩個參考點 在兩個參考點之間hinge連接
一個點做軸 一個點和圓柱coupling 在做軸的點上加UR3 即可實現公轉
最近碰到一個計算時間比較長的Explicit模型,分析步4個,每個分析步計算時間都超過8小時,由于用筆記本計算,下班就得中斷計算關機走人,這樣就不得不考慮重啟動的問題。
經過這幾天的摸索,發現對于這種計算時間長的Explicit模型,要想是實現意外中斷以后繼續計算,最好的方法就是用recover。
方法其實很簡單:
要想在中斷以后實現無縫連接(接著上一個increment繼續計算),必須在計算之前做一個Restart的設定,CAE:step模塊,output->restart requests->intervals,默認為1,將這個數據改大一點,我改為100,與field output的inteval對應起來。有了這個前提,不管什么時候中斷,只需把該job type改為recover,然后提交就會接著中斷的部分繼續算。再中斷再提交就行了,不需要做任何改動,直至計算完成。Command命令:abaqus job=jobname recover int。
解釋:
restart前提條件的設定,Explicit默認的會寫出各種重啟動需要的文件,只不過默認的interval=1,就是說只有在每個分析步的開始和結尾才有重啟動點,這樣不能實現哪中斷哪重啟。要實現無縫連接,理論上這個數據越大越好,只要不超過一個分析步increment總數。但是大了會增加數據寫入量,也沒必要。這個前提設定和restart一樣。
好處:
只需要在第一次重啟動的時候改job type為recover,提交就ok了,當然要記得刪除lck文件。后面中斷了就只需提交就行了,不需做任何改動,方便。而且odb,sta等文件都是在原文件的基礎上續寫,不存在合并odb文件的問題。用這個方法,甚至可以人為中斷后再接著計算。
一點體會,供參考。
1、首先需要重啟動的case計算時有設置重啟動輸出,即inp中以下類似代碼關鍵有write: *Restart, write, number interval=8, time marks=YES
2、查看計算結果后綴為.sta的文件,觀察最后restart number數值,如Restart Number 8 at 3.0600;
3、合并原計算結果文件,jobname為斷電前使用的文件名(需先刪除原文件后綴為.lck的文件)
abq6101 job=jobname convert=all int
4、設置重啟動讀入值,interval等于最后restart number減去1,例 *RESTART, READ, STEP=4, END STEP, INTERVAL=7
5、設置重啟動inp文件,新建一個AMPLITUDE文件,時間設置為計算總時間,*Amplitude, name=Amp-3,TIME=total time
3.0, 0.,3.06,1.6、改變相應邊界條件; 斷電前使用
*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
注意:
默認使用的AMPLITUDE中time=step time,跟AMPLITUDE相關的邊界條件都需要改變;
對時間增量步理解
abaqus的step里有maximum number of increment、initial increment、minimum increment、maximum increment四個量許多網友不知怎樣設置合理,合理設置是建立在深刻理解基礎上的。
要理解這個問題,首先需要了解abaqus的計算過程和有限元計算收斂性問題,abaqus首先用initial值輸入進行疊代計算,如果計算結果收斂,則繼續以這個值代入計算下一步,如果不收斂,則自動減小時間步長(time increment)重新計算直到收斂然后計算下一步。
但是如果時間步長減小到最小值minimum時計算結果還是不收斂,則abaqus將停止計算,由此可知maximum值和minimum值分別是abaqus在收斂計算時時間步長的上下限,同時total time=求和(time increment*number),當時間步長很小時,需要計算的步數number相應增大(電腦計算花的時間也隨之增大),因此number一般要設置較大值。minimum并不是越小越好,因為1)number即計算時間增大2)abaqus計算精度約在10^(-5),當時間步長小于這個值,計算結果已經沒什么意義了。
有限元計算收斂性與(最小空間步長/時間步長)值有關,若minimum設為10^(-5),還是不收斂,可適當減小空間步長(即把網格畫細點),當然還有一些其他辦法,如果實在計算不了,也許是模型本身有點問題,或改為顯示explicit計算
總而言之,maximun number要適當設置較大值,initial可適當改小(如-2,-3量級),minimum(-5量級)不要修改,maximum值影響不大,可不改
顯示分析的求解時間與以下幾個方面有關: 材料密度。加快求解有一種方法就是增大密度,密度增加100倍,時間降低為原來的十分之一。當然,這種情況下必須是密度這一屬性對你要求的結果影響不大的前提下。最小的網格尺寸.這決定計算機每步求解的最大步長。網格尺寸越小,最大時間步長越大。即使只有一個網格尺寸很小,其他都很大,最大時間步長也提不高。所以,避免產生個別小的網格尺寸是十分必要的。網格數量。網格數量越大,每求解一步所需的時間越長。因此,可以適當優化網格,采用局部加密的方法,減少網格數量。4 模型復雜程度,主要是指模型是否有接觸、大變形、多物理場耦合等。越復雜,時間越長。
以上,希望對你有用。
發現當轉動模型時模型就會變成靜態線框下的透明體,如何才能把他改成轉動時不是透明的呢?