abaqus計算線程的案例
多線程會影響Abaqus計算精度嗎？
Abaqus域級并行
一般情況下，多線程運行Abaqus/Explicit時，線程數(shù)并不會影響計算結(jié)果。
然而，有時運行同一個inp文件，當(dāng)域級并行的拓撲域數(shù)目設(shè)置的不一樣時，我們發(fā)現(xiàn)，會存在計算結(jié)果不一致的問題，這是如何發(fā)生的呢？
先來看一下拓撲域的設(shè)置，顯式計算時，并行設(shè)置有Domain和Loop兩種方式，Domain即拓撲域，適用于大多數(shù)顯式計算，創(chuàng)建任務(wù)時，可以把它設(shè)置為線程數(shù)的整數(shù)倍，后處理時，也可以通過Color Code顯示這些拓撲區(qū)域。
域級并行設(shè)置和后處理顯示
Abaqus將劃分好的拓撲域均勻地分配給每個線程，所有域內(nèi)單獨進行計算，每一個時間增量內(nèi)，有共同邊界的域之間會進行通信，合并數(shù)據(jù)。
設(shè)置了8個拓撲域的金屬點陣沖擊模型
Domain如何影響計算結(jié)果？
在分析過程中，不同的分解形式會造成數(shù)值浮點運算次序的改變，由此產(chǎn)生不同的截斷誤差，也就會導(dǎo)致不同的計算結(jié)果。
展開 AMD EPYC 128核心256線程 CPU計算服務(wù)器/GPU服務(wù)器仿真計算、HPC計算、大數(shù)據(jù)分析、
適用場景:
CAE/仿真計算： 如Fluent, Abaqus, ANSYS等，能極大縮短求解時間。
大數(shù)據(jù)與數(shù)據(jù)分析： 海量內(nèi)存和多核心能輕松處理TB級數(shù)據(jù)集。
人工智能與機器學(xué)習(xí)： 適合模型訓(xùn)練和推理，尤其適合中等規(guī)模或作為大型集群的一個計算節(jié)點。
科研計算： 在物理、化學(xué)、生物、氣象等領(lǐng)域進行復(fù)雜的數(shù)值模擬。
虛擬化與云計算： 可以創(chuàng)建大量的虛擬機，作為私有云或虛擬桌面的主機。
媒體與娛樂： 用于三維渲染、視頻編碼等任務(wù)。
SimuFact.Forming V9 多線程計算設(shè)置說明
目前進行鍛壓、擠壓分析的4大主流的商用軟件：Deform、Simufact.Forming（原MSC.SuperForge）、Forge、Qform中，都支持多核運算；現(xiàn)在主流的CPU一般也都在2個核心以上了，所以計算時最大程度的利用CPU的計算資源，能夠大幅度的節(jié)省計算時間。下面就SimuFact.Forming的FE求解器的多線程(多核/多CPU)求解設(shè)置過程；
默認情況下，SimuFact.Forming只支持一個核心或者說一個線程進行運算（比如雙核四線程CPU），在這里需要首先了解一下概念：
A：SimuFact.Forming的FE求解器為MARC求解器，F(xiàn)V求解器為 Dytran；
B：Marc求解器支持多電腦、多線程(多核/多CPU)運算，設(shè)置時需要控制到每一個點；
C：進行多核運算需要License授權(quán)；
了解了以上概念后，我們對SimuFact.Forming進行設(shè)置以支持多CPU運算：
加圖太麻煩，看附件吧
20110128 Simufact.Forming多核運算詳解.pdf
展開 ABAQUS 在計算到step3后計算中斷，變形過大
我在原本可以計算的模型的基礎(chǔ)上修改的，只是刪了幾個樓板與梁連接的栓釘
Mesh Free-眼鏡剛強度校核計算，附Abaqus計算結(jié)果對比
采用Mesh Free對某品牌眼鏡整體剛度、強度進行校核，如下圖所示，Mesh Free支持在不用做幾何清理的前提下進行計算分析，導(dǎo)入模型部件可以包含細節(jié)特征，比如螺釘上的倒角。
眼鏡定義了5種線彈性材料：鏡架主體采用鈦合金；眼鏡片采用樹脂；螺釘?shù)冗B接件采用鋼；鼻托和鏡片扎線采用兩種不同的尼龍材料。
眼鏡腿一只固定，另一只向上掰，加力1N，模擬分析此種工況下眼鏡整體結(jié)構(gòu)的剛度、強度。
Mesh Free所有接觸面定義為完全剛性連接，Abaqus作同樣處理，不考慮非線性因素，對比二者的線性計算結(jié)果。
Mesh Free給出的眼鏡最大變形為23.92mm，Abaqus的結(jié)果為23.46mm。
Mesh Free給出的眼鏡最大應(yīng)力為303.4MPa，Abaqus的結(jié)果為308.3MPa。
談?wù)凪esh Free使用感受：
雖然我常用ABQ，但是不得不說，對于包含細節(jié)幾何特征的復(fù)雜裝配結(jié)構(gòu)建模分析，Mesh Free真的要比Abaqus高效的多。
據(jù)我了解Mesh Free的非線性也在大力的開發(fā)之中，目前已經(jīng)支持經(jīng)典塑性材料非線性、邊界條件非線性也可以設(shè)置滑動和一般的摩擦接觸。
對不熟悉常規(guī)有限元操作的結(jié)構(gòu)設(shè)計人員來說，不用幾何清理、不用劃網(wǎng)格是極好的體驗。
關(guān)鍵是Mesh Free的結(jié)果也確實很準，目前的CAE無非是追求更準的基礎(chǔ)上算的更快，這兩點Mesh Free無疑是滿足的。
Mesh Free
Abaqus
展開 ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結(jié)果） ￥3.7
計算表格（如下）中標黃部分的參數(shù)可自行設(shè)定后，EXCEL程序會自動計算“抗拉強度、非彈性應(yīng)變、受拉損傷因子”。
ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結(jié)果）
這是根據(jù)GB50010-2010中混凝土結(jié)構(gòu)設(shè)計規(guī)范中的混凝土本構(gòu)模型，結(jié)合文獻所述的損傷因子定義，編制的計算C30混凝土非彈性應(yīng)變和損傷因子的EXCEL表格。也是邊學(xué)變做，希望能和大家多交流。
C30砼本構(gòu)（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計算依據(jù).rar
隨機振動分析-abaqus（附一個電池包計算案例） ￥20
四、如何將時域隨機振動曲線轉(zhuǎn)換得到功率譜密度曲線
五、 隨機振動分析理論
附.常見功率譜密度曲線給出形式
附.以dB/oct形式給出的功率譜密度曲線如何計算
附.國標中定義的PSD譜總均方根加速度值是如何計算的？
六. 隨機振動分析案例-abaqus
第一步：計算結(jié)構(gòu)模態(tài)，輸出位移和應(yīng)力。
第二步：隨機振動分析
2.1 定義輸出頻率上下限和模態(tài)阻尼
2.2 定義PSD載荷及加載
2.3 定義輸出
2.4 隨機振動計算頭文件設(shè)置
2.5 隨機振動分析結(jié)果
2.6 隨機振動σ應(yīng)力結(jié)果評價
Abaqus接觸非線性在有限元計算分析中的應(yīng)用 附莊茁ABAQUS非線性有限元分析與實例下載
來源：有限元在線
ABAQUS的非線性主要在有三種：幾何非線性，材料非線性以及接觸非線性。接觸非線性在ABAQUS的有限元計算分析中應(yīng)用非常廣泛，特別是動態(tài)顯式的求解，只要模型中包含兩個以上相互接觸的部件，就要用到接觸非線性。
ABAQUS接觸非線性的設(shè)置主要在Interation模塊中完成，設(shè)置接觸的屬性時，可以設(shè)置摩擦系數(shù)，阻尼系數(shù)，損壞，失效準則等非線性參數(shù)，如圖1所示。
如圖2所示，在接觸定義界面，可以選擇通用接觸、面-面接觸、自接觸等各種非線性接觸方式。
在接觸編輯界面，可以選擇機械約束方式為運動學(xué)接觸算法，或是懲罰接觸方式，還可選擇滑移方式為有限滑移或小滑移，如圖3所示。
這是對模型定義非線性接觸后得到的分析結(jié)果，以供參考。
下載地址：莊茁ABAQUS非線性有限元分析與實例
展開 自動調(diào)用abaqus 進行計算的插件 ￥30
在使用abaqus進行仿真過程中，如果是通過別的軟件進行前處理導(dǎo)出得到inp，常常使用的方法是新建一個bat 文件，在文件中輸入
call abaqus job=XXX
XXX是inp的名字，那么問題來了，在調(diào)試過程中每次都得輸入inp的名字，如果在調(diào)試過程中，在一個文件夾下可能出現(xiàn)多個inp，每次修改起來依然不夠方便。因而，筆者編寫一個插件，每次能自動調(diào)用該文件夾下最新的inp 文件計算。
使用者每次在文件夾中右鍵，即可彈出該按鈕，單擊該按鈕即可自動調(diào)用，如下圖所示
使用方法：
解壓附錄文件，將
附錄：
如何在windows 添加右鍵
https://jingyan.baidu.com/article/3065b3b6455d6dbecff8a4b2.html
展開 abaqus簡單立方體胞元周期性邊界條件施加計算腳本源代碼 ￥39.9
</p><p class="ql-align-justify">在abaqus中通過file→run script選擇該腳本運行，腳本會輸出X，Y，Z三個方向的計算結(jié)果ODB文件。</p><p><br></p>
ABAQUS流體計算
誰有ABAQUS的流體計算資料，越詳細越好， 初學(xué)者。
Abaqus仿真計算中的單元選擇
目前第一、二期直播已結(jié)束（聯(lián)系文末客服看回放），第三期直播<Abaqus仿真計算中的單元選擇>，已經(jīng)開啟報名，歡迎參加~
”
對于有限元分析的網(wǎng)格模型，不僅需要較高的網(wǎng)格質(zhì)量，還需要擁有合適的單元類型。ABAQUS為用戶提供了豐富的單元庫，幾乎可以模擬實際工程中任意幾何形狀的有限元模型，在對一個問題進行分析時，可以根據(jù)情況選擇使用。
如何才能選取出適合于分析的單元類型呢？
本次分享首先介紹ABAQUS中對于單元的分類，每種單元特定的使用范圍，各種單元類型的節(jié)點數(shù)目、單元形狀、插值函數(shù)階次以及單元構(gòu)造的方式。
展開 Abaqus計算內(nèi)存的問題：aba_sint_cap
因此，查閱了網(wǎng)絡(luò)上相關(guān)帖子和幫助文檔，加上親身經(jīng)歷，整理如下：
方法一：可以直接在abaqus_v6.env中設(shè)置memory='80%'或其他方式，詳見幫助文檔
方法二：可通過高級系統(tǒng)設(shè)置--->環(huán)境變量--->變量名：ABA_SINT_CAP 值：16384
（說明：16384這個值可以更大，但似乎沒有意義！）
基于上述方法設(shè)置了足夠大的內(nèi)存后，但仍然出現(xiàn)該錯誤。具體原因：
Abaqus中每8GB的物理內(nèi)存，使用的內(nèi)存將不會超過2GB
在滿足上述條件的基礎(chǔ)上，最多使用內(nèi)存限制為16GB。
遇到這個問題后，建議重新調(diào)整模型，減小inp文件的體量。目前尚未找到特別好的方法！
展開 Abaqus+PyQt+Python平面變形歐拉角計算
下面以簡單例子介紹平面變形、指向歐拉角的計算，包括絕對歐拉角、相對歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側(cè)3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計算要求
計算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。
3 數(shù)據(jù)處理
使用平面節(jié)點坐標、位移數(shù)據(jù)計算平面變形歐拉角。可以使用Python腳本輸出平面節(jié)點編號、節(jié)點坐標（X、Y、Z）、節(jié)點位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數(shù)據(jù)文件，打開的文本文件中7列數(shù)據(jù)為節(jié)點編號、坐標、位移。
三個平面10個工況的節(jié)點數(shù)據(jù)文件如下圖所示。每個文件中包含一個工況一個平面的節(jié)點編號、坐標、位移數(shù)據(jù)。
4 絕對歐拉角計算
使用PyQt+Python開發(fā)了一個簡單的小軟件，計算絕對歐拉角、相對歐拉角。
首先計算各平面的絕對歐拉角。
計算平面1的10個工況的絕對歐拉角。
平面1變形的絕對歐拉角計算結(jié)果如下圖所示。
伴隨絕對歐拉角計算結(jié)果，軟件同時寫出了平面變形前后的坐標系數(shù)據(jù)，如下圖。每行18個數(shù)據(jù)，每3個數(shù)據(jù)為一個坐標軸向量，變形前后2個坐標系，6個坐標軸，18個數(shù)據(jù)。
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當(dāng)于平面1、平面3相對與平面1的相對歐拉角。
計算結(jié)果如下圖所示。
6 小結(jié)
上述軟件用的算法申請了發(fā)明專利，軟件申請了軟著。
展開 