2025.01.10 更新到python3，注意區(qū)分下載！

0.出發(fā)點

數(shù)值仿真數(shù)據(jù)的合理處理對于結果非常重要，特別是細觀模型，包括特殊結構、3D打印材料、顆粒、纖維填充等代表性單元，需要通過提取其平均值才能與試驗進行對比。下面是個簡單的算例，在y方向壓縮極小曲面之Gyroid。

為了對比該極小曲面的應力水平，采用同樣的材料做了單軸壓縮，兩種情況對比如圖所示：

從圖中可以看到，如果僅去極小曲面上的一個點作為其應力應變，其曲線甚至比實心立方體還高，顯然這是不合理的。出現(xiàn)這種現(xiàn)象的主要原因是，此類細觀結構或變形不均勻時，不能取一個點代表整個模型，不然很可能會出現(xiàn)較大的誤差。對此問題，細觀力學有些方法，比如作用力反作用力法，體積平均法等，但也有人認為對于細觀結構，作用力反作用力法 不太合理。體積平均法的簡單表達式如下：

即模型中每一個單元的應力（應變）對單元體積積分后，除以模型整個體積。上述應力應變曲線也證實，采用該方法能夠得到較為真實的數(shù)據(jù)。

如果進提取某個方向的平均應力-應變，也許采用作用力反作用關系更快些；然而，有時候需要不同參量，如von mises應力等的平均值，此時需要通過python進行批量處理。

****基于以上需求，做了各abaqus后處理插件，可以提取任意參量的平均值。******

介紹

插件主界面如下圖所示，功能主要包括：

• 定義任意step；
• 可以選擇每隔n個幀提取，實際上這是為了加快提取速度；
• 可以提取標量：當所輸出為標量時，比如用戶自定義SDV1等，members可以不填寫；同時支持多個變量，比如SDV1,SDV2等；
• 有份量的參量：比如應力S及真實應變LE等都是有方向的，此時不僅需要給出參量名稱，也需要給定方向，比如要提取S11 則在Name填寫S，members填寫11; 如果相提取mises平均，則在Name先填寫S，members 再填mises即可。同樣支持多個量同時提取，比如 Name:S,E 而members：22同時提取S和E的2方向。

• 可以自定義單元集合ElementSet或直接從屏幕選擇：二選一；建議先定義好ElementSet，因為直接選取的話也會消耗時間（單元數(shù)量很多時）；如果前期沒有定義set，后續(xù)將odb以可讀寫方式打開（技巧：直接運行插件一次即可以可讀寫方式打開），如圖所示，再定義set即可：

• 可以對結果轉換坐標輸出：支持前處理或后處理自定義坐標后轉換輸出，如圖所示。

• dVi==1：當前處理沒有輸出EVOL或只想提取算數(shù)平均數(shù)時勾選，此時提取速度快；如果需要提取體積平均如式(1)，則需勾掉，并必須有EVOL。
• 可以動態(tài)顯示：

***前段時間共享源代碼，后來發(fā)現(xiàn)別人做成插件倒賣，所以做成了簡單的插件；希望共同學習，尊重自己、尊重別人*****

***附件是插件，一般可以放到C:\Users\用戶名\abaqus_plugins下即可***

***打開abaqus,菜單欄中點擊plug-ins,里面找出Average即可*****

歡迎咨詢討論。