abaqus返回方法的案例
Abaqus用戶子程序USDFLD調(diào)用GETVRM 返回值為0
SDEG云圖是有數(shù)值的但是調(diào)用一直為零麻煩各位大佬幫忙看看 以下是源代碼 SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER, 2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO, 3 LACCFLA) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 CMNAME,ORNAME CHARACTER*3 FLGRAY(15) DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3), 1 T(3,3),TIME(2) DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*), 1 COORD(*) C real alfa,K0,D,DAMAGED_K,K PARAMETER(TOLER=1.0E-6) alfa=4.0 K0=1.0E-6 C Absolute value of current strain: CALL GETVRM('SDEG',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) D=ARRAY(1) IF (D.LE.TOLER) THEN D=0 END IF DAMAGED_K=EXP(alfa*D) K=K0*DAMAGED_K C field variable FIELD(1)=D C state variable STATEV(1)=D STATEV(2)=K C If error, write comment to .DAT file: IF(JRCD.NE
展開(kāi) SHPB可控多脈沖加載技術(shù)與Abaqus仿真方法 ￥15
1、問(wèn)題介紹
SHPB多脈沖加載方法一般有兩種：多次反射加載法、多級(jí)撞擊桿法。多次反射加載法，利用入射桿的反射波在端面二次反射形成加載波，實(shí)際上常規(guī)的SHPB試驗(yàn)都是多次反射加載，只不過(guò)在處理數(shù)據(jù)時(shí)只截取了第一次加載的數(shù)據(jù)，其特點(diǎn)是相鄰加載時(shí)間間隔是固定值（入射桿桿長(zhǎng)的兩倍與桿彈性波速的比值）；多級(jí)撞擊桿法，是基于撞擊桿或者加載結(jié)構(gòu)設(shè)計(jì)，將撞擊桿設(shè)計(jì)成可實(shí)現(xiàn)多次撞擊的結(jié)構(gòu)，撞擊間隔可調(diào)可控，多級(jí)撞擊桿一般有串聯(lián)結(jié)構(gòu)、夾心結(jié)構(gòu)等形式。
本案例主要介紹SHPB夾心結(jié)構(gòu)的多級(jí)撞擊桿技術(shù)與仿真方法。
2、內(nèi)容
2.1 基于夾心撞擊桿的多脈沖加載SHPB結(jié)構(gòu)
夾心撞擊桿形式的多脈沖加載SHPB結(jié)構(gòu)如下：
夾心形式的撞擊桿主要由外桿和內(nèi)桿組成，內(nèi)桿與外桿端面間隔d。實(shí)際試驗(yàn)中，內(nèi)桿是圓柱體，尺寸與外桿內(nèi)徑相同（留有公差），內(nèi)桿與外桿可以滑動(dòng)，外桿自由端封閉，靠近撞擊端的端面裝配有端蓋。
進(jìn)行實(shí)驗(yàn)時(shí)，內(nèi)桿、外桿以相同的初速度運(yùn)動(dòng)，由于間隔d的存在，外桿先撞擊入射桿，然后經(jīng)過(guò)一定的時(shí)間間隔后內(nèi)桿再撞擊入射桿，因此通過(guò)調(diào)節(jié)間隔d的大小可以控制多脈沖加載的時(shí)間間隔。
2.2 時(shí)間間隔計(jì)算
根據(jù)一維應(yīng)力波理論，可知：
（1）加載脈寬：
第一次加載（加載波1）：
第二次加載（加載波2）：
（2）兩次沖擊時(shí)間間隔：
其中，初始撞擊速度，撞擊外桿長(zhǎng)度，撞擊內(nèi)桿長(zhǎng)度，波速，間隔長(zhǎng)度，為波阻抗比值。
展開(kāi) 在Abaqus中Python通過(guò)findAt方法定位建立region區(qū)域的方法匯總
在Abaqus中Python通過(guò)findAt方法定位建立region區(qū)域的方法匯總
（轉(zhuǎn)帖自 <）" _src="http://blog.sina.com.cn/s/blog_4fff04900102vff3.html>）" style="font-family:宋體, SimSun;font-size:14px;">http://blog.sina.com.cn/s/blog_4fff04900102vff3.html>）
Abaqus腳本接口中許多命令都包含region參數(shù)，包括：
1）Load命令：使用region指定施加荷載的區(qū)域。例如集中力施加在頂點(diǎn)上，壓力施加在邊或面上。
2）Mesh命令：使用region指定單元的類型、網(wǎng)格種子的定義區(qū)域。
3）Set命令：使用region參數(shù)指定集合的區(qū)域。例如：節(jié)點(diǎn)集，單元集等。
4）Coupling命令：使用region參數(shù)指定需要耦合的區(qū)域或者參考點(diǎn)，節(jié)點(diǎn)等。
如果在Abaqus/CAE中，Abaqus將自動(dòng)為模型的特征進(jìn)行編號(hào)，例如節(jié)點(diǎn)編號(hào)，單元編號(hào)，表面中各條邊的編號(hào)，參考點(diǎn)的編號(hào)。但在腳本中無(wú)法獲知這些模型特征的ID，而且對(duì)某個(gè)模型特征編輯后ID將發(fā)生變化，新舊版本ID不同不兼容。
鑒于定位尋找模型特征（頂點(diǎn)（vertex）、邊（edge）、面（face）、體（cell）、參考點(diǎn)（referencePoint））的重要性——可以用來(lái)施加荷載邊界條件，劃分網(wǎng)格、建立集合、建立耦合，有必要通過(guò)坐標(biāo)的方法尋找定位這些模型特征！
findAt方法就可以達(dá)到這一點(diǎn)。用法參考《Abaqus Scripting Reference Manual》。
展開(kāi) 使用abaqus中CEL方法模擬氣囊充氣過(guò)程 ￥49.9
image_process=/format,webp/resize,w_760" data-initial-src="https://img.jishulink.com/202406/attachment/1b75fafef7124891a8a37aecf85b8b23.png">
</figure>
</div><p>2、設(shè)置材料</p><p>氣囊材料選用常規(guī)線彈性材料，可根據(jù)實(shí)際需要選取其他材料模型（本案例教程只為說(shuō)明仿真方法，材料參數(shù)并不準(zhǔn)確）</p><p>歐拉計(jì)算域材料為氣體</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202406/attachment/f7543c5d42164a548fd8f8de9635817a.png" style="text-align: center">
<img src="https://img.jishulink.com/202406/attachment/f7543c5d42164a548fd8f8de9635817a.png" data-mobile-src="https://img.jishulink.com/202406/attachment/f7543c5d42164a548fd8f8de9635817a.png?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/202406/attachment/f7543c5d42164a548fd8f8de9635817a.png?
展開(kāi) 
ABAQUS案例—ABAQUS中聲固耦合、聲輻射分析方法 ￥4
本案例（附件中inp）講述了ABAQUS中的聲固耦合分析、聲輻射分析方法。ABAQUS中有一套完整的聲固耦合分析方法。
噪聲輻射分析中，需要模擬附著在結(jié)構(gòu)上的外部空氣，而且它是向外無(wú)限延伸的，因此直接用聲學(xué)有限單元去模擬無(wú)限的空氣區(qū)域是不合理的。在Abaqus中可以通過(guò)兩種方式來(lái)模擬無(wú)限聲學(xué)介質(zhì)的影響：一,使用聲學(xué)無(wú)限單元；二，用阻抗邊界來(lái)模擬。
在對(duì)外部的噪聲輻射問(wèn)題進(jìn)行仿真分析時(shí)，無(wú)限單元法的應(yīng)用已經(jīng)越來(lái)越廣泛。無(wú)限單元可以直接在結(jié)構(gòu)上定義，或者也可以在聲學(xué)有限單元區(qū)域的終面上定義。
對(duì)于邊界阻抗技術(shù)，實(shí)質(zhì)上屬于無(wú)反射邊界條件。然而當(dāng)用此來(lái)模擬結(jié)構(gòu)外部的區(qū)域時(shí)，結(jié)構(gòu)與輻射表面的距離必須足夠大（通常取聲波波長(zhǎng)的1/3）。
聲學(xué)無(wú)限單元計(jì)算公式與聲輻射阻抗邊界的計(jì)算有幾個(gè)關(guān)鍵的區(qū)別：無(wú)限單元采用更高階的差值函數(shù)，而聲輻射邊界則采用一階差值函數(shù)。雖然無(wú)限元計(jì)算每個(gè)單元的花費(fèi)更高，但是無(wú)限單元的要比阻抗邊界精確很多，因此通過(guò)減小無(wú)限元的單元規(guī)模，從而可以大大的降低結(jié)構(gòu)總的計(jì)算時(shí)間；本案例即是講解無(wú)限元單元法在模擬噪聲分析中的應(yīng)用。
展開(kāi) Abaqus中Fortran子程序的調(diào)試方法 附ABAQUS2020關(guān)聯(lián)Fortran下載
近兩年，收到許多同學(xué)關(guān)于Abaqus軟件使用過(guò)程中Fortran子程序調(diào)試的問(wèn)題，在的錯(cuò)誤面前許多初學(xué)者無(wú)從下手，不知如何去查找問(wèn)題所在。本文將對(duì)Abaqus中子程序的調(diào)試方法以及如何對(duì)錯(cuò)誤進(jìn)行預(yù)警做一個(gè)簡(jiǎn)單的介紹。
通過(guò)總結(jié)各類錯(cuò)誤問(wèn)題，大致可以把子程序編寫(xiě)過(guò)程中的常見(jiàn)錯(cuò)誤歸為三類：模型錯(cuò)誤、語(yǔ)法錯(cuò)誤和邏輯錯(cuò)誤。
（一）模型錯(cuò)誤：
有限元模型中自身材料、截面屬性、網(wǎng)格類型、網(wǎng)格劃分方式、接觸設(shè)置、邊界/載荷設(shè)置、輸出設(shè)置、軟件子程序配置等引起的錯(cuò)誤，非子程序本身的錯(cuò)誤。此類錯(cuò)誤一般會(huì)在任務(wù)監(jiān)控界面下的Error或Warning欄中有具體提示。
模型錯(cuò)誤屬于有限元模型常規(guī)建模錯(cuò)誤，可以先使用非子程序模型予以測(cè)試，排除建模錯(cuò)誤之后，再改為子程序計(jì)算模型。該類錯(cuò)誤屬于最容易解決的一類。
（二）語(yǔ)法錯(cuò)誤：
在使用FORTRAN編寫(xiě)子程序時(shí)，使用了非法的語(yǔ)句或者非法的格式，而引起的錯(cuò)誤。該類錯(cuò)誤會(huì)直接在任務(wù)監(jiān)控界面提示Error:Problem during compilation編譯出錯(cuò)。如下圖所示。
具體的語(yǔ)法錯(cuò)誤內(nèi)容將會(huì)顯示在與任務(wù)名同名的XXX.log文件中，該文件位于Abaqus當(dāng)前的工作目錄，語(yǔ)法錯(cuò)誤點(diǎn)及所在行數(shù)都有明確的提示。
例如，將以下代碼中第29行注釋掉，第28行的語(yǔ)句被打斷，變得不完整，將會(huì)有哪些提示信息呢？提示信息又將在哪里顯示呢？
使用文本編輯軟件（notepad++或UE等）將log文件打開(kāi)，可以考單提示信息如下，具體的錯(cuò)誤類型和錯(cuò)誤位置都有明確的提示：
語(yǔ)法錯(cuò)誤是許多子程序初學(xué)者經(jīng)常會(huì)犯的錯(cuò)誤，按照上述方式就可以快速定位錯(cuò)誤位置了，至于基本的Fortran語(yǔ)言語(yǔ)法，可以參照響應(yīng)的語(yǔ)法書(shū)。
展開(kāi) Abaqus 非線性屈曲分析方法 附ABAQUS分析手冊(cè)分析卷下載
當(dāng)然，對(duì)于方筒這類實(shí)際上是通過(guò)顯示方法實(shí)現(xiàn)的，更準(zhǔn)確的講是動(dòng)力屈曲分析，所以我們還得判斷動(dòng)能、塑形耗散等能量參數(shù)，才能使結(jié)果更加準(zhǔn)確。
下載地址：ABAQUS分析手冊(cè)分析卷
abaqus系列技巧2：如何在abaqus中用掃掠的方法畫(huà)六面體網(wǎng)格
在abaqus中畫(huà)網(wǎng)格并不是一件快樂(lè)的事情，很多時(shí)候回比較苦惱，尤其是我們需要一個(gè)六面體網(wǎng)格的時(shí)候。作者對(duì)待網(wǎng)格的策略是，不太復(fù)雜的網(wǎng)格選擇在ab中完成，復(fù)雜的在hypermesh中完成。當(dāng)然這個(gè)復(fù)雜的邊界時(shí)很模糊的，每個(gè)人都不一樣。
在六面體的劃分選項(xiàng)里面，其實(shí)還是略為簡(jiǎn)單的。如下圖所示：
兩種方法。由于structured無(wú)可調(diào)節(jié)項(xiàng)，反正我是基本不同，sweep由于給予了一定的調(diào)節(jié)空間，通過(guò)合理的選擇參數(shù)，還是能滿足需要的。
這里面主要講下掃掠的幾個(gè)要素：
一個(gè)掃掠需要三個(gè)主要要素，源面，目標(biāo)面和掃掠路徑，缺一不可。如常見(jiàn)的正方體，圓柱體都是。但是圓錐體不是，因?yàn)樵趫A錐體上你找不到源面和目標(biāo)面。當(dāng)然源面和目標(biāo)面并不一定要求一樣大，但一定要“相對(duì)”。其次關(guān)于掃掠路徑，一定要連續(xù)，光滑，不能有折線的情況。
如下圖，就不可以直接用掃掠完成劃分。
這時(shí)候我們需要對(duì)其進(jìn)行切分，把這三個(gè)要素都湊齊
在abaqus中，掃掠的三個(gè)要素一般只需要指定路徑，制定的方法如下：
到這里，基本就可以了。上面的algorithm，挨個(gè)試下，哪個(gè)漂亮用哪個(gè)，沒(méi)必要太在意。
我的視頻課程中有一些具體的劃分案例，有興趣的可以看一看
想獲得幻想飛翔最新CAE技術(shù)文章，請(qǐng)關(guān)注幻想飛翔公眾賬號(hào)：幻想飛翔CAE。
想獲得幻想飛翔最新CAE技術(shù)文章，請(qǐng)關(guān)注幻想飛翔公眾賬號(hào)：幻想飛翔CAE。
也歡迎加入abaqus交流群516073058進(jìn)行討論研究
展開(kāi) 【Abaqus爆炸】 SPH方法分析彈頭侵徹爆破 ￥69
</span></h2><p>前段時(shí)間分享了用Abaqus CEL方法分析侵徹爆破問(wèn)題，詳情見(jiàn)：</p><p><br></p><p><br></p><div contenteditable="false" width="100%">
<figure class="figure-link" data-title="【Abaqus爆炸】 CEL方法分析彈頭侵徹爆破" data-link="https://www.yqgqt.org.cn/post/1790179" data-regular="true">
<a href="https://www.yqgqt.org.cn/post/1790179" target="_blank">
<div class="link-card">
<span class="link-title">【Abaqus爆炸】 CEL方法分析彈頭侵徹爆破</span>
<div class="link-url">
https://www.yqgqt.org.cn/post/1790179
</div>
</div></a>
</figure>
</div><p>然而有人反饋說(shuō)歐拉區(qū)域填充比較麻煩，而且網(wǎng)格尺寸非常影響材料填充精度。另外，對(duì)于復(fù)雜結(jié)構(gòu)，仍要求網(wǎng)格是六面體網(wǎng)格；也不支持質(zhì)量縮放（歐拉單元），這可能是CEL方法的不足。實(shí)際上上述侵徹爆破問(wèn)題可以用SPH方法進(jìn)行分析，而且SPH方法可以跟普通的拉格朗日單元分析類似進(jìn)行操作，下面展示SPH方法的使用及注意事項(xiàng)。
展開(kāi) Abaqus python二次開(kāi)發(fā)方法 附Abaqus Python Reader v1.9.4.1
復(fù)制你新生成的python文件，并對(duì)該文件中的參數(shù)進(jìn)行修改，在提交給abaqus--cae就可以完成重復(fù)建模了，如此可以無(wú)限重復(fù)，其實(shí)python語(yǔ)言都是大白話，你能看懂的
需要指出的是：
1） 可以結(jié)合其它編程語(yǔ)言如 VB、VC 配合修改參數(shù)并生成py文件，使用更為靈活
2） 生成py文件可以直接在cae中選擇file-run script，選擇你生成的python文件
3） 可以用python文件直接生成cae模型文件，可在py文件最后添加
"mdb.saveAs(pathName='" *** "')"
4） 可以通過(guò)cmd命令直接將py文件提交個(gè)abaqus內(nèi)核，讓abaqus進(jìn)行運(yùn)算，cmd命令為“Shell "C:\Windows\SysWOW64\cmd.exe /k abaqus cae noGUI=" **** ".py ", vbHide
等待abaqus運(yùn)算的py語(yǔ)言"myJob.submit(consistencyChecking=OFF, datacheckJob=True)"
下載地址：Abaqus Python Reader v1.9.4.101028
展開(kāi) 疲勞分析|Abaqus Goodman方法案例操作 附ABAQUS疲勞分析簡(jiǎn)介下載
Abaqus/View結(jié)果讀取
讀取分析歷程中的最大交變應(yīng)力和最小交變應(yīng)力云圖
新建場(chǎng)變量：Alternating Stress和Mean Stress
根據(jù)公式：
在Abaqus后處理新建場(chǎng)變量
輸出場(chǎng)變量值到Excel
針對(duì)新建場(chǎng)，輸出單元積分點(diǎn)對(duì)應(yīng)的交變應(yīng)力和平均應(yīng)力，并輸出到Excel，與Goodman圖一并繪制。
上圖，
仿真所得單元積分點(diǎn)落到
曲線的上方或下方，
處于上方為疲勞壽命沒(méi)達(dá)到
臨
界曲值
10
E5
次。
下載地址：ABAQUS疲勞分析簡(jiǎn)介
Abaqus流固耦合仿真方法 附ABAQUS初學(xué)者用戶子程序小例子下載
這種方法是通過(guò)粒子間的物理距離、平滑距離h來(lái)確定“內(nèi)核方程”，進(jìn)而確定周?chē)Ｗ訉?duì)中心粒子自由度的影響。
應(yīng)用：
液體晃動(dòng)、波浪、噴漆、水壓切割、氣流、填塞、破碎后的二次沖擊，鳥(niǎo)撞分析、射流爆破等。
方法：
把有限元網(wǎng)格轉(zhuǎn)化為SPH粒子，可以按照應(yīng)力、應(yīng)變或時(shí)間閾值來(lái)確定SPH粒子的轉(zhuǎn)化。
Tips：上圖是Abaqus自帶的流固耦合經(jīng)典案例，幫助文檔搜索“Impact of a water-filled bottle”即可找到水壺跌落CEL/SPH兩種方法的inp文件。當(dāng)初第一次看到CEL方法很震驚，下載完模型怎么也找不到里面的水在哪！（奧秘就藏在VFT工具里）~為了搞懂，我把這個(gè)案例說(shuō)明翻譯了一遍，很有收獲。
三大類方法在流固耦合方面的表現(xiàn)：
根據(jù)具體計(jì)算需求，結(jié)合各種方法的特點(diǎn)，選擇最合適的手段進(jìn)行流固耦合仿真分析。
下載地址：ABAQUS初學(xué)者用戶子程序小例子
展開(kāi) Abaqus材料庫(kù)及其導(dǎo)入方法
abaqus材料庫(kù).zip
abaqus新建材料庫(kù)
進(jìn)入Property模塊時(shí)，CAE界面左側(cè)常用于顯示模型樹(shù)的一欄中會(huì)出現(xiàn)材料庫(kù)Material Library。
點(diǎn)擊Creat，定義新的材料庫(kù)。
輸入材料庫(kù)名稱，選擇存儲(chǔ)位置(Home位于ABAQUS的數(shù)據(jù)定義的目錄下，例如在我電腦的位置是：C:\SIMULIA\CAE\2019\win_b64\lib\abaqus_plugin中，Current是當(dāng)前工作目錄下，如：D:\temp\abaqus_plugins)。
如圖所示可以任意調(diào)用材料庫(kù)中的材料。
PS. 大家有沒(méi)有做abaqus增材制造模擬的可以一起討論交流
展開(kāi) ABAQUS批處理方法 ￥5
ABAQUS也可以像ANSYS一樣采用批處理方法，該方法可以打開(kāi)ABAQUS界面，也可以直接在后臺(tái)進(jìn)行求解。這樣的好處是可通過(guò)其他軟件工具創(chuàng)建界面程序，然后調(diào)用批處理文件進(jìn)行后臺(tái)求解，也是ABAQUS非常靈活的二次開(kāi)發(fā)的一種形式。
1. 首先，編寫(xiě)Python的腳本語(yǔ)言程序，本歷程中采用腳本語(yǔ)言錄制的方法來(lái)生成一段腳本語(yǔ)言，如下所示。在最后一行添加運(yùn)行該腳本的程序。
ABAQUS中沖擊動(dòng)力學(xué)問(wèn)題的求解方法
在顯式方法中可以很容易的模擬接觸條件和其他一些極度不連續(xù)的情況，并且能夠一個(gè)節(jié)點(diǎn)、一個(gè)節(jié)點(diǎn)地求解而不必迭代。為了平衡接觸時(shí)的內(nèi)力和外力，可以調(diào)整節(jié)點(diǎn)加速度。
顯式方法最顯著的特點(diǎn)是沒(méi)有在隱式方法中所需要的整體切向剛度矩陣。由于是顯式地前推模型的狀態(tài)，所以不需要迭代和收斂準(zhǔn)則。
2.4 顯式方法的條件穩(wěn)定性
穩(wěn)定性限制了ABAQUS/Explicit求解器所能采用的最大時(shí)間步長(zhǎng)，這是ABAQUS/Explicit 進(jìn)行計(jì)算的一個(gè)重要因素。應(yīng)用顯式方法，基于在增量步開(kāi)始時(shí)刻
的模型狀態(tài)，通過(guò)時(shí)間增量
前推到當(dāng)前時(shí)刻的模型狀態(tài)。使得狀態(tài)能夠前推并仍能保持對(duì)問(wèn)題的精確描述的時(shí)間非常短。如果時(shí)間增量大于這個(gè)最大的時(shí)間步長(zhǎng)，則時(shí)間增量已經(jīng)超出了穩(wěn)定性限制（stability limit）。超過(guò)穩(wěn)定性限制的后果是數(shù)值不穩(wěn)定，可能導(dǎo)致計(jì)算不收斂。由于一般不可能精確地確定穩(wěn)定性限制，因而采用保守的估計(jì)值。為了提高計(jì)算效率，ABAQUS/Explicit選擇時(shí)間增量，使其盡可能地接近而且又不超過(guò)穩(wěn)定性限制。
在系統(tǒng)中，以最高頻率（
）的形式定義穩(wěn)定性限制。
展開(kāi) 