
發(fā)布
注冊
/
登錄ABAQUS角度輸出的案例
基于ABAQUS的python不同噴丸角度的隨機彈丸噴丸插件
Python實現(xiàn)可變噴丸角度的隨機彈丸分布
噴丸是一種典型的強化手段,采用高速彈丸反復(fù)撞擊材料表面,使表層發(fā)生塑性變形,從而引入顯著殘余壓應(yīng)力,增大零件表面硬度,細化晶粒,有利于提高材料的抗疲勞性能。
噴丸覆蓋率和噴丸角度作為噴丸工藝中最重要的工藝參數(shù)之一,直接影響噴丸工藝的加工質(zhì)量。覆蓋率C表示為被噴零部件表面上的彈痕面積與零部件總面積的比值,可通過阿夫拉米公式計算:
單彈丸模擬結(jié)果
r為單個彈丸撞擊后留下的彈坑半徑;N 為彈丸個數(shù)。
ABAQUS軟件作為一款強大的非線性軟件,在沖擊領(lǐng)域具有重要應(yīng)用。本帖主要通過ABAQUS的python二次開發(fā)實現(xiàn)機械噴丸強化過程,并根據(jù)單彈丸噴丸后的結(jié)果結(jié)合阿夫拉米公式近似得到一定覆蓋率下的彈丸個數(shù)。
噴丸速度與噴丸強度有關(guān),實際實驗中以ALMEN試片的弧高值表征噴丸強度,在模擬中以100%覆蓋率下的噴丸速度表征噴丸強度,下面是ALMEN試片噴丸后的變形特征。
ALMEN試片噴丸后的位移變化
下圖為通用機械噴丸插件,可以根據(jù)覆蓋率輸入彈丸參數(shù):彈丸個數(shù)、彈丸大小,彈丸本插件采用的是可變形體,需要材料參數(shù),噴丸角度以90度最佳,小角度噴丸在相同噴丸強度下容易出現(xiàn)損傷;靶材參數(shù):靶材大小與靶材的材料參數(shù)。其次還有分析步時間以及網(wǎng)格大小,分析步時間其實可以通過噴丸速度與最高彈丸位置點計算,本插件仍然可以優(yōu)化,彈丸采用解析剛體減少計算量。
機械噴丸插件
下面是一些插件的案例。
展開 abaqus中旋轉(zhuǎn)角度?
abaqus中我想給個長方體兩端面固定,然后給長方體除了兩端面添加旋轉(zhuǎn)角度,這個長方體會發(fā)生變形,最后想看變形后的力?該如何做?
abaqus中旋轉(zhuǎn)角度?
abaqus中我想給長方體兩端面施加固定約束,除了長方體兩端面添加旋轉(zhuǎn)角度,看長方體發(fā)生變形,最后看變形中的力?該如何設(shè)置?
從編程角度闡述有限元最佳入門方法:以Abaqus 和OpenSees 為例
而我們常見的商業(yè)有限元軟件Abaqus和ansys 的求解模塊均由Fortran 語言編寫的,這也是為什么,我們經(jīng)常看到很多學(xué)者基于自己課題研究的需要,通過Fortran語言編寫了很多材料的umat 或者vumat的子程序,實際上就是對abaqus 的材料庫進行了擴充。但是Fortran 語言相比C在調(diào)用計算機計算資源這方面,卻遠遠不及,這也是為什么,我們發(fā)現(xiàn)相同的自由度規(guī)模的分析,Opensees 的速度較Abaqus 快很多。更不用說,Opensees 中擁有最為強大的材料庫和單元庫,每年都會有很多來自全球的學(xué)者開發(fā)相當(dāng)?shù)拇a來擴充OS 的求解模塊。
而對于有限元軟件的前后模塊,則對計算性能要求不高,所以通常用其他的編程語言進行GUI的設(shè)計,比如我們常用的abaqus 就是采用python 語言進行前后處理器的設(shè)計的。前處理模塊,主要功能就是實現(xiàn)和使用者的交互,你在窗口上的點擊交互,會激發(fā)它自動生成某些代碼,比如你在用鼠標(biāo)操作abaqus的每一步,都會在對應(yīng)的abaqus.py文件中輸出相應(yīng)的代碼如圖1:
圖1 Abaqus.py中對應(yīng)GUI 操作代碼
最終abaqus 會自動把這些由GUI出發(fā)產(chǎn)生的代碼轉(zhuǎn)換為通用格式ASCII .inp 文件,由這些命令流去觸發(fā)abaqus 求解模塊,使其采用相應(yīng)的材料,單元,作用,算法,荷載,來形成有限元模型,并對其求解。
上述分析,也就說明,我們在建立abaqus 模型時,可以采用三種方式建立模型,一是GUI操作,二是基于python 寫模型的script, 三是直接寫.inp 文件。 但是對于上述過程,對于初學(xué)者而言,或者大多數(shù)已經(jīng)使用abaqus 的用戶,如果不做二次開發(fā)的話,想必對后兩者的方法都相當(dāng)陌生。
展開 
abaqus二次開發(fā):后處理批量提取場輸出和歷程輸出結(jié)果(源碼帶注釋) ¥198
abaqus二次開發(fā):后處理批量提取場輸出和歷程輸出結(jié)果(購買后需要加V私聊獲取使用教程)
購買后,請加V:wzd_1021_
更多插件內(nèi)容,請私聊
同時提供定制hypermesh/hyperview/abaqus前后處理插件開發(fā)。
【ABAQUS模態(tài)動力學(xué)】Composite&abaqus 預(yù)應(yīng)力模態(tài)分析&輸出單元剛度矩陣
從上面這個理解出發(fā),ABAQUS預(yù)應(yīng)力模態(tài)只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣:
【ABAQUS 二次開發(fā)筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解,不一定完全正確,本文僅作為個人學(xué)習(xí)記錄之用,其他概不負責(zé)。
Abaqus-高清圖片輸出 ¥4.99
<p>在使用 Abaqus 的過程中,輸出高清圖片是一個常見的需求。</p><p>高清圖片的本質(zhì)是 DPI(<span style="background-color: rgba(18, 18, 18, 0);">Dots Per Inch,圖像每英寸長度內(nèi)的像素點數(shù))達到一定數(shù)值。</span>網(wǎng)上的很多資料,僅僅是調(diào)整輸出圖片的格式,或者調(diào)整 Abaqus 內(nèi)的顯示設(shè)置,并沒有從根本上改善圖片的 DPI。</p><p>本帖提供了操作簡便的處理思路,可以精準(zhǔn)的控制圖片的清晰度及尺寸,適用于 原始文件(cae文件)的幾何模型、裝配模型以及結(jié)果文件(odb文件)的應(yīng)力云圖、變形圖等圖片的高清輸出。</p><p>高清圖片導(dǎo)出--圖文流程,見付費部分。</p>
展開 【Abaqus】輸出任意參量平均值的插件
</p><p class="ql-align-justify">****基于以上需求,做了各abaqus后處理插件,可以提取任意參量的平均值。******</p><h2 class="ql-align-center"><strong>介紹</strong></h2><p class="ql-align-justify">插件主界面如下圖所示,功能主要包括:</p><p class="ql-align-justify"><br></p><ul><li class="ql-align-justify">定義任意step;</li><li class="ql-align-justify">可以選擇每隔n個幀提取,實際上這是為了加快提取速度;</li><li class="ql-align-justify">可以提取標(biāo)量:當(dāng)所輸出為標(biāo)量時,比如用戶自定義SDV1等,members可以不填寫;同時支持多個變量,比如<span style="color: rgb(18, 18, 18);">SDV1,SDV2等;</span></li><li class="ql-align-justify">有份量的參量:比如應(yīng)力S及真實應(yīng)變LE等都是有方向的,此時不僅需要給出參量名稱,也需要給定方向,比如要提取S11 則在Name填寫S,members填寫11; 如果相提取mises平均,則在Name先填寫S,members 再填mises即可。同樣支持多個量同時提取,比如 Name:S,E 而members:22同時提取S和E的2方向。
展開 ABAQUS 結(jié)果文件輸出
區(qū)別是
a.odb文件默認(rèn)生成,是二進制文件,abaqus才能訪問,python后處理也是通過abaqus訪問。
b.fil文件不是默認(rèn)生成,需要在inp中添加關(guān)鍵字(*FILE FORMAT, ASCII),添加的位置可以是(*Restart)命令之后。是文本文件,文本編輯器可以訪問。
要輸出的量可以通過cae界面選擇,也可以通過inp添加。
輸出到fil文件屬于歷史輸出范圍,輸出的是所有增量步的值。fil文件的大小似乎是固定的,在cae界面中,歷史輸出可以設(shè)置只輸出最后一個increment的值,對應(yīng)在inp文件中表示為(*Output, history, frequency=99999),輸出頻率非常大,故只輸出最后一幀。這個命令針對fil的輸出無用,這里姑且認(rèn)為fil只能輸出的所有增量步的值。
*NODE FILE 輸出節(jié)點變量,位移(UCOORD-101),節(jié)點坐標(biāo)(COORD-107)
*EL FILE 輸出單元變量,MISES應(yīng)力(S),應(yīng)力分量(SINV-12)
要使用URDFIL子程序,首先必須先創(chuàng)建.fil文件,子程序讀取.fil文件,F(xiàn)ortran代碼實現(xiàn)結(jié)果的格式化輸出。
展開 關(guān)于Abaqus圖片輸出的總結(jié) 附ABAQUS 6.14.1安裝包下載和安裝教程下載
關(guān)于Abaqus圖片輸出的總結(jié)
1.改變abaqus界面背景顏色
Aabaqus默認(rèn)的背景為漸變藍色,如果想把背景變?yōu)榘咨奖銓?dǎo)出圖片,則可以進行以下操作。
View-Graphics Options-Viewport Background-Solid(點擊圓圈)-選擇顏色(白色)即可。
2.設(shè)定界面文字是否顯示
1Abaqus界面上有默認(rèn)輸出的文字和坐標(biāo)和圖例等,可以通過以下操作選擇是否顯示這些內(nèi)容。
2Viewport-Viewport Annotation Options-General-Visibility
3取消勾選其中的內(nèi)容即可
3.圖例大小
為了美觀需要顯示較小的圖例,通過以下操作即可
Viewport-Viewport Annotation Options-Legend-Set Font-Size
修改Size的大小即可
4.將云圖中的網(wǎng)格除去
有時為了輸出圖片的美觀,需要將劃分好的網(wǎng)格去除,通過以下操作即可。
Option-Common Plot Options-No edges
點選No edges即可
5.輸出高清png圖片
File-Print
需要高清圖片,點掉256色,進一步選擇PNG Formal Options即可進入分辨率界面。
以上即為Abaqus圖片輸出相關(guān)的操作過程
下載地址:ABAQUS 6.14.1安裝包下載和安裝教程
展開 ABAQUS—輸出骨架曲線
ABAQUS構(gòu)件擬靜力模擬時,輸出滯回曲線的同時,也可以輸出骨架曲線。
可以點下面鏈接去看視頻。
https://mp.weixin.qq.com/s/RvIABJS3VVD8USUCVnpWWw
教程:
1、將加載制度里每級加載的最大位移對應(yīng)的時間找到,寫成一列,復(fù)制。
2、在ABAQUS的Step模塊中,點擊Create Histoty Output,Domain里選擇Set,Set-4(Set-4為輸出骨架曲線的點的集合,自己定義即可)。Frequency里選From times points,點擊Create Time Point,將復(fù)制的時間列輸入進去,如下:
選擇輸出的RF和U。
3、計算完成后,在Create XY Data里點ODB histoty output(第一個),就有骨架曲線力和位移的數(shù)據(jù)了。
點擊Save XY Data,復(fù)制到Excel中,升序,畫圖,和滯回曲線對比,基本吻合,僅黑框處有所差異,原因在于你輸出的是最大位移及其對應(yīng)的力(一般情況下即為骨架曲線點),但是試驗和模擬中,也會存在每級加載中最大的力是在最大位移之前出現(xiàn)的,即產(chǎn)生了不同。
有用請點贊!關(guān)注我,關(guān)注公眾號 結(jié)構(gòu)工程師 ,可以學(xué)習(xí)更多!
展開 
ABAQUS輸出單元積分點坐標(biāo)
方法
在ABAQUS CAE的場輸出中選擇的坐標(biāo)點是節(jié)點的坐標(biāo),而節(jié)點是從積分點插值出來的,單元積分點的信息相對真實。所以最好是獲取積分點的信息,其中積分點的坐標(biāo)無法在CAE中獲取,需要在關(guān)鍵字中添加。具體在每個分析步的單元輸出下面添加COORD,如果需要輸出節(jié)點的坐標(biāo)也可以在節(jié)點場輸出下面添加COORD(這和CAE中場輸出選擇節(jié)點坐標(biāo)的效果是一致的)。具體如下圖:
2.注意
在ODB結(jié)果中創(chuàng)建場輸出時會附帶著一份XYZ坐標(biāo),這個應(yīng)該也可以當(dāng)做單元的坐標(biāo),,但是我比較過這個附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo),有時候有點差別,可能是數(shù)據(jù)精度的問題。
展開 Abaqus 輸出矩陣的方法
在比較深入的CAE分析過程中,我們可能會關(guān)注模型中某些部分的剛度分布情況,因此需要提取剛度矩陣參數(shù),因為使用相對較少,目前的Abaqus前處理模塊CAE中還無法直接實現(xiàn)該設(shè)置,因此需要對其inp文件進行局部修改來定義。
Abaqusn輸出剛度矩陣主要分為兩大類,一類是輸出某一部分的單元的剛度矩陣,一類是輸出整體模型的剛度矩陣(可含單元、彈簧、質(zhì)量、阻尼、力等)。
下面我們分別對這兩中情況分別討論:
1、輸出單元矩陣
輸出單元剛度矩陣的方式相對容易理解,當(dāng)然你也可以從自帶的英文幫助中需找答案,你只需要在inp文件中添加下述關(guān)鍵字即可:
*File Format,ASCII
*Element Matrix Output,Elset=xxx,
File Name=xxx,Frequency=1,Output File=User Defined,Stiffness=Yes
其中第一行的目的是將軟件默認(rèn)生成的二進制文件轉(zhuǎn)化為ASCII格式文件,以方便直接查看和讀取避免亂碼。Elset后面的xxx對應(yīng)的是之前已經(jīng)定義好的單元集的名稱,也就是你想要獲取的那部分單元。File Name后面的xxx是abaqus計算后所提取出的單元矩陣所保存的文件名稱,后綴為.mtx ,也就是說最終在工作目錄中會生成一個xxx.mtx的文件,可以通過記事打開的剛度矩陣文件。
展開 Abaqus Python 二次開發(fā)-歷程輸出數(shù)據(jù)提取的python實現(xiàn) ¥3.49
在Abaqus中,歷程變量的輸出,尤其是反力及位移,經(jīng)常是如下的格式。
此時,如果想提取數(shù)據(jù),一些關(guān)于Abaqus 二次開發(fā)的書籍里提供的方式是:odb.step['Step-1'].historyRegions['Node PART-TRIMMED-MESH.288422']. historyOutputs['U1'].data 。這種方式,需要輸入 'Node PART-TRIMMED-MESH.288422' 這種很長的字符串,繁瑣且容易出錯。
在付費內(nèi)容中,提供了一種便捷的方法,得到長字符串,提取相關(guān)的數(shù)據(jù),并將數(shù)據(jù)保存為csv文件。
展開 abaqus鋼管壓縮及能量輸出
abaqus鋼管壓縮及能量輸出