abaqus 網(wǎng)格信息的案例
Abaqus采用Python腳本提取ODB信息，輸出為VTK的UnstructureGrid文件格式 ￥49.99
在有些場景下，需要將Abaqus計算結(jié)果提取出來保存為VTK格式，用于一些可視化和數(shù)據(jù)處理（VTK在可視化這方面還是強大）。最近寫了兩個腳本，一個是保存VTK非結(jié)構(gòu)網(wǎng)格格式的文件，一個是提取Abaqus計算結(jié)果中的網(wǎng)格信息、field信息，將兩個腳本結(jié)合起來就可以實現(xiàn)提取ODB中的信息，并保存成為VTK非結(jié)構(gòu)網(wǎng)格文件，當然，還是需要有Abaqus環(huán)境。
注意：提取ODB信息的那個腳本目前僅僅是針對我自己提供的那個算例，有能力的朋友完全可以修改一下用于自己的情況。寫VTK非結(jié)構(gòu)網(wǎng)格文件的那個腳本目前應該沒有太大的bug。ODB信息的提取還是需要有Abaqus環(huán)境，并不是解析二進制文件。請看清楚再確定是否需要。
之所以做這個，還是希望慢慢完善，從而積累出一定數(shù)量的腳本，用于ODB數(shù)據(jù)的提取。
以下是abaqus中的網(wǎng)格和云圖，和用paraview查看的vtk文件的網(wǎng)格和云圖：
展開 ICEM CFD拉伸網(wǎng)格的一些信息
有些CFD軟件在計算二維時，仍然需要具有一層網(wǎng)格的假三維網(wǎng)格，如CFX、OpenFOAM等，如果在ICEM CFD中直接采用三維網(wǎng)格劃分，則網(wǎng)格數(shù)會不可接受，因此其拉伸網(wǎng)格功能就非常有用了，下面介紹幾個在生成拉伸網(wǎng)格過程中的幾個信息。
(1) 先生成二維網(wǎng)格。在這個例子中，由于管束周圍的網(wǎng)格非常密，而在遠離管束處，網(wǎng)格不需要那么密，以減小網(wǎng)格數(shù)，因此需要用到比例劃分邊的功能，為了設置合理的比例，最好打開線的方向（勾選curve direction）。
(2) 拉伸網(wǎng)格。將生成的二維網(wǎng)格選上，最好在New side part name中輸入新的part名，如newside，也在New top part name中輸入新part名，如newtop，這樣的目的是為了方便后面定義邊界條件，不然和以前的邊界名攪和在一起，在定義邊界時，難以選擇網(wǎng)格單元。注意：vector不像建模時的需要用“{}”，而是直接輸入方向矢量，如0 0 -1。
Apply之后的網(wǎng)格如下面所示，這里我們生成了一層網(wǎng)格，當然也可以根據(jù)需要生成多層。
(3) 關于網(wǎng)格數(shù)量。生成網(wǎng)格后，我們采用info輸出網(wǎng)格信息，可以看到生成的fluid體積域中包含了193735個單元。但這里總的單元數(shù)Total elements顯示為592608，其包含了面單元、點、線等所有信息，所以不要搞混了，然后將網(wǎng)格導入到fluent中，查看，也可以看到網(wǎng)格總數(shù)為193735。
展開 基于Matlab的有限元網(wǎng)格自動生成算法 | Q4、Q8、Abaqus單元網(wǎng)格
今日給大家?guī)淼闹饕獌?nèi)容是二維問題下四邊形單元有限元網(wǎng)格如何自動生成？
單元網(wǎng)格的形成實際上屬于有限元計算中的前處理部分，即確定單元節(jié)點信息，當模型較為復雜時，用戶可在Abaqus、Ansys等大型商業(yè)有限元軟件中進行建模，導出網(wǎng)格信息。
當模型較為簡單時，如二維平面板模型，用戶可基于一些較為基礎的網(wǎng)格生成算法，在自己的程序中通過控制模型長、寬等信息，即可生成有限元網(wǎng)格。
看似應用有限，但是在一些比較復雜的領域內(nèi)，往往需要先在簡單的模型中得到理論驗證，如此以來，有利于自編程代碼的完整性，即前處理、內(nèi)核計算、后處理于一體。
本篇推文，木木就帶著大家學習一下Q4、Q8單元網(wǎng)格的自動生成以及Abaqus網(wǎng)格節(jié)點順序解讀。
代碼獲取：
基于Matlab的有限元網(wǎng)格自動生成算法 | Q4、Q8、Abaqus單元網(wǎng)格
Q4單元網(wǎng)格
單元自動網(wǎng)格劃分
如下圖所示，為4節(jié)點四邊形單元網(wǎng)格生成示意圖，圖中NXE和NYE分別是模型橫向和縱向單元個數(shù)，dhx和dhy分別是單元的橫向、縱向長度。
展開 用hypermesh劃分網(wǎng)格，用dynaform進行拉延設置總是提示相同的錯誤信息
用hypermesh劃分網(wǎng)格，用dynaform進行拉延設置，提交文件，再用LS_DYNS971分析，不論怎么設置，總是提示相同的錯誤信息，但是
在dynaform中還能運算。
有圖，高手幫幫忙吧。信息提示的好像與內(nèi)存設置和單位有關系。
我的qq：373437611
ABAQUS網(wǎng)格控制屬性詳解（三種網(wǎng)格劃分技術） ￥12
><span style="font-family:'Calibri';font-size:12.0pt;white-space:pre-wrap;">，它ABAQUS是決定采用何種策略劃分網(wǎng)格的選項</span><span style="font-family:'Calibri';font-size:12.0pt;white-space:pre-wrap;">。
用多面體網(wǎng)格劃分壓縮機渦輪，再通過fluent進行求解（case文件中包含fluent所有設置信息） ￥15
渦輪網(wǎng)格
壓力云圖
速度云圖
利用Python提取ABAQUS的.odb結(jié)果中的輪軌接觸信息 ￥500
ABAQUS 顯式explicit計算時所有結(jié)果都將寫入到.odb文件中，在GUI界面進行數(shù)據(jù)提取時效率較低，且有些詳細結(jié)果不能直接從.odb中獲取。因此，需要利用python讀取.odb，然后進行提取關心的數(shù)據(jù)，以此進行下一步計算。
Abaqus視圖界面顯示信息設置
就是部分場景用圖時，部分信息用不上需要隱藏時，具體如何操作：
1.在結(jié)果視圖模塊下。
2.菜單欄Viewport——Viewport Annotation Options...
3.然后彈出下圖中的窗口
4.每條控制界面顯示的地方如下圖所示。
DS.SIMULIA.Abaqus 2022-2026補丁信息
2022最新補丁：2540 大概率不再更新了
2023最新補丁2613：
2024 最新補丁2615
2025最新補丁2610
2026最新補丁2614
DS的Abaqus每年大的新版本網(wǎng)絡公開和諧，但是每個大版本的補丁基本沒人發(fā)布，使用過程有問題的，建議升級個大版本的補丁；
Abaqus cae noGUI 運行 python 腳本，命令行窗口輸出提示信息 ￥5
<p> 通過 <a href="https://www.yqgqt.org.cn/major/abaqus" rel="noopener noreferrer" target="_blank">Abaqus</a> Command 窗口 或者 Windows 的 CMD 窗口，以 <a href="https://www.yqgqt.org.cn/major/abaqus" rel="noopener noreferrer" target="_blank">abaqus</a> <a href="https://www.yqgqt.org.cn/major/cae" rel="noopener noreferrer" target="_blank">cae</a> noGUI=myscript.py 的方式，運行 <a href="https://www.yqgqt.org.cn/major/abaqus" rel="noopener noreferrer" target="_blank">Abaqus</a> 的 Python 腳本，是一個非常方便且高效的操作。但是，這種操作有一個明顯的缺點，那就是 Python 腳本中的 print( ) 語句無法在命令行窗口輸出提示內(nèi)容，著對于 Python 腳本運行情況的監(jiān)控是極其不利且不方便的。</p><p> <a href="https://www.yqgqt.org.cn/major/abaqus" rel="noopener noreferrer" target="_blank">Abaqus</a> 內(nèi)置的 prettyPrintToTerm( ) 方法，雖然可以實現(xiàn)在命令行窗口輸出提示內(nèi)容，但是樣式很丑。
展開 利用Python提取ABAQUS的計算結(jié)果（ODB）信息
圖6在ABAQUS命令行接口（基于Python的Abaqus靜力分析操作實例）表示的應該是最后一次循環(huán)后各變量的值，可以看出RefU表示的nodeSets['U-XYZ']的位移值，圖6中顯示為一個索引；而RefUValues則表示該節(jié)點集中每個節(jié)點的位移值，體現(xiàn)在values中；而values[0]表示第一個節(jié)點的位移值，由于存在三個方向的位移，U1,U2,U3,則后續(xù)通過data[0]，data[1]，以及data[2]分別表示三個方向的位移值。
圖6 運行相關結(jié)果查看（通過ABAQUS命令行接口）
特別的，由于U-XYZ只有一個點存在于Set中，如果我們在這里用values[1]，看會出現(xiàn)什么樣的效果？如圖7。可以明顯看到Sequence index out of range表示超出了范圍，因為RefUValues[1]表示提取第二個節(jié)點，但我們這里只有一個節(jié)點存在于U-XYZ點Set中，所以出現(xiàn)了范圍超過的報錯提示。
圖7 可以存在RefUValues[1]嗎？（通過ABAQUS命令行接口）
需要說一下的，如果大家覺得ABAQUS命令行接口的查看范圍太窄，可以往上拉動窗口，這樣就可以看到更多的代碼，如圖8所示，但與之相反的就是模型就會顯得很小，與我們一般的常規(guī)abaqus操作有點相反。
圖8 ABAQUS命令行接口顯示更多的內(nèi)容
再對代碼中的Python round函數(shù)做個說明，通過下面的例子我們可以看到round函數(shù)是四舍五入函數(shù)，一般為四舍五入后為整數(shù)，round(a,b)中的a為需要進行四舍五入的數(shù)，而b為需要保留的小數(shù)位數(shù)，如圖9所示。
展開 
Abaqus錯誤與警告信息匯總（適合初學者）（轉(zhuǎn)載內(nèi)容）
（四）
比較有價值的的信息考察。比如：
1）Numerical sigularity solver problem. numerical sigularity when processing node105 instance 表示：數(shù)值奇異：剛體位移（欠約束）
2）Zero pivot 表示：過約束
3）對于TIME INCREMENT REQUIRED IS LESS THAN THE MINIMUM SPECIFIED，Too many attamps have been made
4）對于“網(wǎng)格扭曲”的警告： excessively distorted elements 前面有提到。
第一步：采用display查看“ ErrElemExcessDistortion-Step1 ”在模型的哪些部位，做到心中有數(shù)。
第二步：檢查模型的網(wǎng)格質(zhì)量： mesh步---verify----Analysis Check選取模型。這種情況，一開始計算即出現(xiàn)“distorted element”的信息。除此之外，很多其他問題也會網(wǎng)格扭曲警告。比如，幾何模型導入有誤需要修補、單元類型選取錯誤、邊界條件有誤、材料屬性錯誤、接觸設置不合理、子程序錯誤等。
第三步：即使你的網(wǎng)格劃分很好，如果變形過大，也會導致網(wǎng)格扭曲。然后修改網(wǎng)格劃分，不要出現(xiàn)紅色，關鍵區(qū)域不要出現(xiàn)黃色。（當然最好是所有的網(wǎng)格都用structure劃分，且都沒有紅色、黃色出現(xiàn)。網(wǎng)格質(zhì)量就比較好。這種情況，警告信息往往是在計算到一定步驟之后才出現(xiàn)“distorted element”。 這種情況建議采用ALE等方式。
展開 切削仿真（abaqus軟件）中的典型報錯信息及對應解決方案(一)
一、材料問題
報錯提示1：
報錯原因：沒有輸入密度參數(shù)
解決方案：輸入密度參數(shù)
報錯提示2：
報錯原因：密度參數(shù)為零
解決方案：修改密度參數(shù)，使之不為零
報錯提示3：刀具單元顯示為紅色（正常為綠色，如右圖所示），同時提示如下圖所示
▲報錯
▼正常
報錯原因：沒有創(chuàng)建刀具截面
解決方案：創(chuàng)建刀具截面
二、裝配問題
報錯提示4：提交后無報錯信息，但是提交后會出現(xiàn)刀具整體“穿透”工件的現(xiàn)象，如下圖所示：
報錯原因：在assembly模塊刀具和工件分別裝配了兩次。
解決方案：在左邊的目錄樹刪掉一組即可，如下圖所示：
三、單元問題
報錯提示5：
報錯原因：分析步選擇了熱位移耦合step ,但是刀具單元沒有選擇熱位移耦合單元。
解決方案：刀具單有選擇熱位移耦合單元即可，如下圖所示：
四、場變量問題
結(jié)果提示6：
產(chǎn)生原因：沒有勾選狀態(tài)變量
解決方案：在step的場變量中中勾選狀態(tài)變量即可，如下圖所示：
結(jié)果提示7：ODB文件過大
產(chǎn)生原因：場變量輸出幀數(shù)過多，如下圖所示：
解決方案：在step的場變量中，將輸出幀數(shù)適當減少即可，默認為20，推薦值200。
展開 abaqus系列技巧2：如何在abaqus中用掃掠的方法畫六面體網(wǎng)格
在abaqus中畫網(wǎng)格并不是一件快樂的事情，很多時候回比較苦惱，尤其是我們需要一個六面體網(wǎng)格的時候。作者對待網(wǎng)格的策略是，不太復雜的網(wǎng)格選擇在ab中完成，復雜的在hypermesh中完成。當然這個復雜的邊界時很模糊的，每個人都不一樣。
在六面體的劃分選項里面，其實還是略為簡單的。如下圖所示：
兩種方法。由于structured無可調(diào)節(jié)項，反正我是基本不同，sweep由于給予了一定的調(diào)節(jié)空間，通過合理的選擇參數(shù)，還是能滿足需要的。
這里面主要講下掃掠的幾個要素：
一個掃掠需要三個主要要素，源面，目標面和掃掠路徑，缺一不可。如常見的正方體，圓柱體都是。但是圓錐體不是，因為在圓錐體上你找不到源面和目標面。當然源面和目標面并不一定要求一樣大，但一定要“相對”。其次關于掃掠路徑，一定要連續(xù)，光滑，不能有折線的情況。
如下圖，就不可以直接用掃掠完成劃分。
這時候我們需要對其進行切分，把這三個要素都湊齊
在abaqus中，掃掠的三個要素一般只需要指定路徑，制定的方法如下：
到這里，基本就可以了。上面的algorithm，挨個試下，哪個漂亮用哪個，沒必要太在意。
我的視頻課程中有一些具體的劃分案例，有興趣的可以看一看
想獲得幻想飛翔最新CAE技術文章，請關注幻想飛翔公眾賬號：幻想飛翔CAE。
想獲得幻想飛翔最新CAE技術文章，請關注幻想飛翔公眾賬號：幻想飛翔CAE。
也歡迎加入abaqus交流群516073058進行討論研究
展開 輸出abaqus inp如何去除instance信息，而且單元節(jié)點編號也不會有重復
解決的方法：在command line輸入 mdb.models['modelName'].setValues(noPartsInputFile=ON) 重新生成inp文件，這樣在寫出的inp文件里面就沒有instance信息，而且單元節(jié)點編號也不會有重復了。
