abaqus進行計算的案例
#ABAQUS使用vcct進行裂紋擴展計算插件
<p>虛擬裂紋閉合方法‘vcct’，廣泛應用于斷裂力學中裂紋擴展的模擬分析中，虛擬裂紋閉合法是經過對斷裂力學的不斷深入研究而提出的能夠有效計算裂紋尖端能量釋放率的一種極為有效的方法，ABAQUS軟件中在早期并未提供虛擬裂紋閉合法vcct的功能，因此必須自己二次開發程序才能進行順利進行能量的提取計算，否則，需要經過手動處理大量的數據文件才能獲得想要的結果，這個過程是極為復雜的，經過ABAQUS版本的不斷改進，在6.10版本之后的ABAQUS中都加入了vcct的計算，這使得裂紋尖端能量的計算及裂紋擴展的分析又向前邁進了一步，此帖子主要對比使用二次開發的插件及較新版本中自帶的vcct計算之間的結果對比。</p><p>模型建立：</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/201812/a5231bea6dd34249ac610ea6ca0698e3.jpg" title="1.jpg" alt="1.jpg" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/201812/a5231bea6dd34249ac610ea6ca0698e3.jpg?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201812/a5231bea6dd34249ac610ea6ca0698e3.jpg?
展開 自動調用abaqus 進行計算的插件 ￥30
在使用abaqus進行仿真過程中，如果是通過別的軟件進行前處理導出得到inp，常常使用的方法是新建一個bat 文件，在文件中輸入
call abaqus job=XXX
XXX是inp的名字，那么問題來了，在調試過程中每次都得輸入inp的名字，如果在調試過程中，在一個文件夾下可能出現多個inp，每次修改起來依然不夠方便。因而，筆者編寫一個插件，每次能自動調用該文件夾下最新的inp 文件計算。
使用者每次在文件夾中右鍵，即可彈出該按鈕，單擊該按鈕即可自動調用，如下圖所示
使用方法：
解壓附錄文件，將
附錄：
如何在windows 添加右鍵
https://jingyan.baidu.com/article/3065b3b6455d6dbecff8a4b2.html
展開 abaqus批量提交inp文件進行計算
此時在abaqus插入腳本即可使用了
當然還有別的方法，但這中腳本方法親測我覺得最靠譜。
abaqus批處理提交計算3個完成后自動退出不能繼續進行的解決辦法
abaqus批處理提交計算3個完成后自動退出不能繼續進行的解決辦法.pdf
一個奇怪的問題（很多人沒有遇到過）：
Abaqus進行批處理時，批處理文件中有超過3個job時，第三個完成后自動退出，后續的不再進行了，必須設置多個批處理文件才能完成更多作業（每個批處理文件中不超過3個job），最初的設置如下（大家通常使用的）：
只能計算出A1/A2/A3， A4/A5不再進行計算。
解決辦法：
ABAQUS斷裂模擬收徒 ，保證快速學會各種ABAQUS斷裂模擬方法 1200/人（將享有各種插件以及程序，價值3000+、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）
展開 
換熱器流固熱耦合計算，四面體網格多面體網格分開畫好后組裝再進行計算設置(含fluent計算設置視頻) ￥30
外部氣流和內部水流
組裝后的網格
Abaqus、Tosca和Fe-safe聯合仿真進行疲勞優化
1綜述
運用Abaqus、Tosca、Fe-safe聯合仿真，實現產品的疲勞優化。Abaqus進行有限元計算，結果ODB文件導入到Fe-safe中進行疲勞分析，疲勞分析的損傷值作為Tosca形狀優化的優化目標，Tosca對表面節點進行擾動，更新后的inp文件導入給Abaqus，如此循環實現疲勞優化。
圖1疲勞優化流程
2模型準備
2.1ABAQUS模型
有限元分析中采用線性分析，有2個LOADCASE，載荷分別為150MPa、70MPa（如圖2）。由于TOSCA中不支持*Part、*Instance、*Assemble等關鍵字，輸出inp文件時需進行設置，Model>Editattribute>Model>Donotusepartsandassembliesininputfile，如圖3所示。
圖2載荷
圖3輸出設置
導出inp后，寫批處理命令運行inp文件。
callabaqusjob=holeplate_damcpus=4int
2.2FE-SAFE模型
FE-SAFE中疲勞分析設置過程如下圖，導入FEA模型、設置分析集合材料、設置載荷工況，然后進行疲勞分析計算。
圖4FE-SAFE疲勞分析設置過程
疲勞計算完成后，在.\jobs\job_01\fe-results文件下生成holeplate_damResults.odb，last_run.stlx等文件。
展開 利用Python對Abaqus進行后處理結果輸出
利用Python對Abaqus進行后處理結果輸出
1 概述
在Abaqus的二次開發過程中，通常需要采用Python腳本語言將Abaqus的計算結果進行輸出，然后再進行處理。Python使Abaqus的內核語言，使用較為方便，Abaqus運行Python語言的方式有多種，可以直接命令窗口，也可以讀入腳本，還可以采用類似批處理的方式。
本次以一個例子細說Python語言在Abaqus后處理中的應用，模型的計算結果云圖如圖1所示。
圖1 計算結果
2 輸出所有節點的Mises應力
直接上Python代碼：
import os
myodb=openOdb(path='Job-1.odb')
cpFile=open('artlcF1.txt','w')
RF=myodb.steps['Step-1'].frames[1].fieldOutputs['S'].values
for i in range(len(RF)) :
cpFile.write('%10.3F\n' % (RF[i].mises))
else:
cpFile.close()
#引入模塊，因為需要打開結果文件
#打開結果文件，并復制給變量myodb
#打開一個txt文件
#將輸出場賦值給RF
#循環語句，向txt文件逐行寫入mises應力
Abaqus的結構層次分的很細，比如結果文件下分如下：
圖2 Model data
使用過Abaqus的都知道step表示載荷步，frame表示載荷子步，因而在讀取Mises應力時需要詳細地指定輸出哪一步的應力，而應力結果是輸出場數據（fieldOutput）的中一種，需要指定是何種應力，程序才知道怎么讀取并寫入。
展開 使用C＃模塊進行自定義計算
VirtualLab Fusion中的C＃模塊為用戶提供了完全的自由度，可以用于不同的計算目的。 它們可以基于給定的公式快速實現相對簡單的任務，同時，C＃模塊保留與VirtualLab Fusion內部典型文檔交互的能力。 例如，我們演示了如何實現C＃模塊中光導耦合的光柵周期計算（基于給定公式），以及VirtualLab Fusion內兩個場之間偏差的計算。
測量給定結果的準確性是科學和工程的基礎。 在本用例中，我們將向您展示如何對模塊進行編程以計算兩個場之間的標準偏差。
在波前計算中OpticStudio如何進行采樣
這和FFT的計算方式以及OpticStudio中方向的定義有關。當網格的中心點在某個域（例如空間域）的坐標為n/2+1時，在另一個域中（例如頻域）的中心點坐標則為n/2。仔細觀察下圖你會發現最左一列和最下一行的數據為空白。
在MTF計算中，MTF為波前的自相關函數，并且通常來講其像素數量為波前圖的兩倍（不考慮坐標軸變化）。因此為了MTF，OpticStudio首先會用數據0將32x32個數據點補充為64x64個數據點，然后再進行自相關計算。對于三維FFT MTF（Surface FFT MTF），OpticStudio會對波前的FFT平方，然后再對其進行FFT計算，換句話說MTF是PSF的傅里葉變換。將三維FFT MTF參數設置如下：
我們得到如下結果：
可以看到峰值點在坐標（32,32）處，或者為（n/2,n/2）。
OpticStudio通過自相關函數的邊界1/(lambda*F/#)來確定三維FFT MTF的頻率間隔，其中lambda為系統中最短的波長（如果我們計算的是多波長結果）。OpticStudio通過實際計算所有波長的截止頻率和F數的乘積，并依據他們中最大的結果來縮放整個圖表。其他波長的數據則會在光瞳空間中縮放，來使所有的PSF以同樣的間距進行采樣。
展開 【simufact.welding】無法進行求解計算
如圖，simufact.welding，在提交計算后，求解器不運行。無論是軟件自帶的模型還是自己導入的模型都無法求解，且系統檢查無警告，無錯誤。
請有相關經驗的前輩指教如何解決。
最大實體要求如何在尺寸鏈中進行計算？
下面我們通過孔軸裝配案例來進行分析，如下圖所示，來計算基準面B處的間隙大小，分析孔軸零件是否能順利裝配。
孔、軸的被測要素均有最大實體要求，對于這種情形如何來進行尺寸鏈計算？首先我們分別對孔、軸零件的圖樣標注進行解析。
軸：
圖4a 圖示解釋
圖4b 動態公差圖
孔：
圖5a 圖示解釋
圖5b 動態公差圖
所以由最大實體要求的定義和動態公差圖可知，最大實體要求應用在被測要素上，實際上是在一定范圍內放松了對幾何公差的要求，當被測要素偏離最大實體尺寸時，偏離的值可以補償給幾何公差。
在尺寸鏈計算時，由于幾何公差有最大實體要求，會使得孔軸的邊界尺寸發生變化，所以將幾何公差疊加至尺寸公差進行尺寸鏈計算，下面我們利用專業尺寸鏈計算及公差分析工具（DCC）進行計算。
展開 
如何用概率法進行尺寸鏈手工計算？
在GB/T5847 -2004中規定了極值法和概率法兩種計算方法，大多數工程師做手工計算時常采用極值法進行計算，概率法相對復雜一些，很多工程師不太熟悉，下面這個案例我們通過GB/T5847-2004中的概率法來計算。
某機械結構如圖1所示，A/B零件采用人工裝配，裝配過程中保持A/B零件側面平行，分析間隙C1大小，是否可以順利裝配？
圖1
已知：零件尺寸如圖2和3所示，所有尺寸的分布規律為正態3西格瑪。
圖2
圖3
這是一個典型的孔軸裝配案例，A、B兩板通過右側孔軸配合，當A板向右移動至極限位置，即可使左側間隙C1達到最大值，此時C1若大于零，即可保證順利裝配。根據此裝配關系繪制尺寸鏈圖如圖4所示。
展開 FLUENT中進行共軛傳熱計算
由于流體求解器同時具備流體與固體傳熱計算的能力，因此可以直接采用流體求解器進行求解，無需使用流固耦合計算。流體求解器能夠求解流體對流、傳導、輻射傳熱，對于固體傳熱計算，只能求解熱傳導方程。
本例演示共軛傳熱問題在FLUENT中的求解方法。
1、問題描述
如圖1所示的計算區域，既包含流體區域也包含固體區域。在初始狀態下，流體域與固體與溫度均為293K，然后給固體域底部施加恒定溫度434K，計算分析計算域內溫度隨時間分布規律。邊界條件如圖中所示。
圖1計算域描述
2、建立幾何模型并劃分網格
利用DM建立如圖1所示2D平面幾何。采用全四邊形網格劃分，如圖2所示。
為所有邊界命名，尤其是流體和固體區域交界面，后面需要在求解器中進行設置。
3、進入Fluent求解設置
本例為瞬態計算。
涉及到熱量傳遞，因此需要激活能量方程。
流體介質為理想氣體，考慮其在溫度影響下密度變化。
考慮重力影響，設置重力加速度向量[0,-9.81,0]，設置操作密度為0。如圖3所示。
壓力-速度耦合方程采用PISO求解方式，對流項計算采用QUICK算法，其他項采用二階迎風格式。
圖2網格模型
圖3 操作項設置面板
設置流體域介質為air，固體域介質為默認的AL。
按圖1所示邊界條件設置計算域邊界。
創建交界面，如圖4所示進行設置。
圖4 設置交界面
4、初始化計算
設置初始化溫度293K，如圖5所示。
圖5初始化面板
設置自動保存選項與動畫錄制項。
設置時間步長0.1s，時間步數100，內迭代次數20。
進行求解計算即可。
展開 [NEWSLETTER] 使用C＃模塊進行自定義計算
VirtualLab Fusion中的C＃模塊為用戶提供了完全的自由度，可以用于不同的計算目的。 它們可以基于給定的公式快速實現相對簡單的任務，同時，C＃模塊保留與VirtualLab Fusion內部典型文檔交互的能力。 例如，我們演示了如何實現C＃模塊中光導耦合的光柵周期計算（基于給定公式），以及VirtualLab Fusion內兩個場之間偏差的計算。
用于光導耦合的光柵周期計算的模塊
生成VirtualLab Fusion中的模塊，計算耦合光柵周期的范圍，以滿足光導的傳導條件。
編程計算標準偏差的模塊
測量給定結果的準確性是科學和工程的基礎。 在本用例中，我們將向您展示如何對模塊進行編程以計算兩個場之間的標準偏差。
了解更多信息，請發送郵件至：support@infotek.com.cn / support@infocrops.com
展開 workbench自動調用CFX進行求解計算
在workbench框架下使用命令行自動調用CFX進行建模、劃分網格和流體域設置，設置好之后 如何提交計算？（包括如何設置并行計算以及結果文件的保存目錄）
