abaqus溫度調(diào)用的案例
ABAQUS UMAT調(diào)用后為什么損傷變量出現(xiàn)負(fù)值？
ABAQUS UMAT調(diào)用后為什么損傷變量和單元?jiǎng)h除出現(xiàn)負(fù)值和大于1得值？
調(diào)用子程序后，計(jì)算復(fù)合材料損傷過程，損傷變量和單元?jiǎng)h除出現(xiàn)負(fù)值和大于1得值，真誠求助。附上子程序，請(qǐng)老師指教
20251203.txt
fortran 調(diào)用 abaqus
matlab和python都可以方便的調(diào)用<a href="/major/abaqus,為算法優(yōu)化帶來了極大的便利。我在優(yōu)化過程中需要用到Fortran調(diào)用abaqus,更準(zhǔn)確的說是在abaqus的子程序中調(diào)用abaqus，目前遇到的一個(gè)問題是log中會(huì)有一個(gè)warning,暫不清楚會(huì)有什么影響，下面開始寫如何做到的。
我的warning：*Warning: it took 18 seconds to delete C:\Users\Bidong\AppData\Local\Temp\Bidong_Job-4_13268
1.Fortran調(diào)用abaqus代碼：
program XXX
use dflib
implicit none
logical(4) cmdout
cmdout=system('abaqus job=Job9_A2M inp=D:\abaqus\software\temp\Job9_A2M.inp')
end program XXX
這段代碼可以運(yùn)行指定文件夾的inp,最后的結(jié)果文件在F文件項(xiàng)目所在的文件夾中。
2.在abaqus job1的子程序中運(yùn)行新的job2，將上述代碼寫成一個(gè)子程序，用call的形式調(diào)用。生成的job2結(jié)果文件會(huì)在job1運(yùn)行結(jié)束后在內(nèi)存中釋放，job2的結(jié)果可以在job2的程序或子程序中進(jìn)行提取。
3.在遇到已經(jīng)存在的任務(wù)時(shí)，abaqus會(huì)詢問是否覆蓋先前的 job，如果不想被詢問直接覆蓋先前的任務(wù)，可以在上述命令結(jié)尾加上 ask=off
展開 調(diào)用ABAQUS內(nèi)置JH2模型模擬沖擊損傷-ABAQUS例子
結(jié)果如下：
impactsiliconcarbide_jh2.txt
把附件的txt后綴直接改為inp文件即可運(yùn)行
ABAQUS斷裂模擬收徒 ，快速學(xué)會(huì)各種ABAQUS斷裂模擬方法 **/人（將有機(jī)會(huì)享有各種插件以及程序，價(jià)值**、專門定制視頻、全程親自教學(xué)、各種模型調(diào)試及解答問題等等，傾囊相教）
Abaqus圓形激光溫度-位移耦合案例教學(xué) ￥19.98
1、 引言
本案例通過力 - 熱耦合分析方法，探究圓形激光載荷作用下玻璃板的溫度分布及應(yīng)力響應(yīng)特性。通過開發(fā)定制化子程序生成激光熱源，并結(jié)合溫度 - 位移耦合分析步，建立高精度有限元模型，最終實(shí)現(xiàn)對(duì)溫度場與應(yīng)力場的多物理場耦合求解與結(jié)果分析。
2、 幾何模型與材料參數(shù)
（1） 模型構(gòu)建：建立三維實(shí)體模型模擬玻璃板，尺寸為178×127×0.3（需根據(jù)實(shí)際場景設(shè)定具體參數(shù)），
圖1模型構(gòu)建
（2） 材料屬性：定義玻璃板的熱物理參數(shù)（如導(dǎo)熱系數(shù)、比熱容、熱膨脹系數(shù)）與力學(xué)參數(shù)（如彈性模量、泊松比），考慮材料屬性隨溫度的非線性變化（如需）。
圖2 材料屬性構(gòu)建
3、 激光熱源子程序開發(fā)
（1） 熱源特性：采用高斯分布模擬圓形激光束，功率密度函數(shù)為：
其中，P 為激光功率，r0為光斑半徑，r 為徑向坐標(biāo)
（2） 子程序?qū)崿F(xiàn)：基于ABAQUS的用戶子程序接口（如DFLUX或HETVAL），編寫 Fortran/Python 程序生成動(dòng)態(tài)加載的圓形激光熱源，通過時(shí)間 - 空間函數(shù)控制熱源移動(dòng)軌跡（如需模擬掃描過程）。
圖3 使用荷載子程序
5、 計(jì)算結(jié)果與分析
（1） 溫度場分布特征
1． 云圖可視化：通過后處理軟件顯示不同時(shí)刻的溫度場云圖，典型結(jié)果包括：激光光斑中心區(qū)域出現(xiàn)局部高溫峰值，溫度梯度沿徑向快速衰減；隨時(shí)間延長，熱擴(kuò)散導(dǎo)致高溫區(qū)域擴(kuò)大，穩(wěn)態(tài)時(shí)形成穩(wěn)定溫度分布。
2． 數(shù)據(jù)提取：提取特征點(diǎn)（如光斑中心、邊緣）的溫度 - 時(shí)間曲線，分析升溫速率與峰值溫度隨激光功率 / 作用時(shí)間的變化規(guī)律。
圖7 溫度云圖可視化
（2） 應(yīng)力場響應(yīng)規(guī)律
1．
展開 
自動(dòng)調(diào)用abaqus 進(jìn)行計(jì)算的插件 ￥30
在使用abaqus進(jìn)行仿真過程中，如果是通過別的軟件進(jìn)行前處理導(dǎo)出得到inp，常常使用的方法是新建一個(gè)bat 文件，在文件中輸入
call abaqus job=XXX
XXX是inp的名字，那么問題來了，在調(diào)試過程中每次都得輸入inp的名字，如果在調(diào)試過程中，在一個(gè)文件夾下可能出現(xiàn)多個(gè)inp，每次修改起來依然不夠方便。因而，筆者編寫一個(gè)插件，每次能自動(dòng)調(diào)用該文件夾下最新的inp 文件計(jì)算。
使用者每次在文件夾中右鍵，即可彈出該按鈕，單擊該按鈕即可自動(dòng)調(diào)用，如下圖所示
使用方法：
解壓附錄文件，將
附錄：
如何在windows 添加右鍵
https://jingyan.baidu.com/article/3065b3b6455d6dbecff8a4b2.html
展開 ABAQUS中調(diào)用MATLAB的函數(shù)
實(shí)例測試
在運(yùn)行程序時(shí)，需要注意的是，必須將用到的自定義的m文件的目錄都包括進(jìn)去
from mlab.releases import latest_release as matlab
# 添加matlab程序的位置，需要引入所有需要依賴的m文件的目錄
x = matlab.path(matlab.path(),r'C:\Users\15321\abaqus_plugins\path')
# 運(yùn)行函數(shù)
number = matlab.plot(x)
ABAQUS顯式子程序調(diào)用規(guī)則及nblock變量解釋
ABAQUS在調(diào)用VUHARD子程序時(shí)，每次向子程序提供136個(gè)單元（單元數(shù)＜136則提供所有單元）進(jìn)行計(jì)算，使用(do k=1,nblock → end do)計(jì)算每一個(gè)單元的相關(guān)變量。因此對(duì)于1000單元來說，一共調(diào)用8次子程序。
設(shè)置全局變量commom /globals/ kdtest，在ABAQUS每一次調(diào)用子程序之后，給其加一，統(tǒng)計(jì)模擬過程中的總循環(huán)數(shù)，在(do k=1,nblock → end do) 循環(huán)內(nèi)部輸出變量kdtest的值。（全局變量可以不跟著k的循環(huán)而變化，用戶可以根據(jù)需求設(shè)置其在代碼中的功能），結(jié)果截圖如下：
共1000個(gè)數(shù)據(jù)，1~7各重復(fù)了136次，8重復(fù)了48次，與上述分析一致。
因此在nblock實(shí)際代表的是ABAQUS提供給子程序的材料點(diǎn)塊，這個(gè)塊區(qū)包含的單元數(shù)與模型單元數(shù)有關(guān)，而k則是對(duì)該材料點(diǎn)塊實(shí)現(xiàn)一個(gè)遍歷，確保每個(gè)單元都被考慮到。本研究僅針對(duì)于單核計(jì)算來講，多核模擬將在后續(xù)展開介紹。
展開 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
展開 Abaqus調(diào)用內(nèi)置子程序模擬形狀記憶合金 ￥19.89
形狀記憶合金在外力下產(chǎn)生塑性變形，去掉外力后變形不能完全恢復(fù)，但將合金加熱到一定的溫度后，其變形消失，恢復(fù)到原始形狀。
形狀記憶合金最典型的特征包括形狀記憶效應(yīng)和超彈性。形狀記憶效應(yīng)是指通過加熱使材料溫度達(dá)到Af以上時(shí)，馬氏體相變?yōu)閵W氏體，材料最終恢復(fù)原始形狀，如圖1所示
圖 1 形狀記憶效應(yīng)示意圖
超彈性則是指，當(dāng)材料所處的環(huán)境溫度高于奧氏體相變結(jié)束的臨界溫度 Af，材料處于奧氏體相的熱力學(xué)穩(wěn)定狀態(tài)，馬氏體處于不穩(wěn)定狀態(tài)。在此溫度下，應(yīng)力誘發(fā)的馬氏體相變與溫度誘發(fā)的馬氏體相變有所不同，應(yīng)力誘發(fā)相變所產(chǎn)生的馬氏體并不是自相適應(yīng)的，材料會(huì)產(chǎn)生很大的變形(大于 5%)；卸載以后材料產(chǎn)生的變形將會(huì)完全消失。超彈性示意圖如圖2所示
圖2 超彈性示意圖
為了在Abaqus中模擬形狀記憶合金的形狀記憶效應(yīng)以及超彈性行為，我們可以通過編寫Umat/Vumat子程序來實(shí)現(xiàn)。但是由于編寫子程序需要很高的門檻，同時(shí)也需要花費(fèi)大量時(shí)間精力，因此本文向大家介紹了一種直接調(diào)用Abaqus內(nèi)部SMA材料本構(gòu)的方法。
SMA內(nèi)置本構(gòu)的調(diào)用方法與自編子程序相比更加便捷，無需安裝Fortran開發(fā)環(huán)境。同時(shí)Abaqus內(nèi)置的SMA子程序適用于隱式分析和顯示分析。
通過Abaqus模擬得到的SMA單向拉伸載荷位移曲線如下所示
展開 abaqus材料庫插件、材料庫修改、材料庫的調(diào)用
1概要
ABAQUS擁有強(qiáng)大的非線性處理能力，但是不提供材料庫，每次都需要去查找、并重復(fù)輸入，大大降低了工作效率，考慮到這一點(diǎn)，ABAQUS提供了材料庫接口，以*.lib文件形式進(jìn)行存儲(chǔ)。
POLARIS_MAT_BASE是星辰-北極星團(tuán)隊(duì)開發(fā)的一款基礎(chǔ)材料庫插件，共包含318種材料，主要涉及材料密度、彈性模量、塑性、熱膨脹系數(shù)、比熱、熱傳導(dǎo)率。參數(shù)由網(wǎng)絡(luò)資源轉(zhuǎn)換而來，避免不了可能存在的錯(cuò)誤，還請(qǐng)查證后使用。如您發(fā)現(xiàn)錯(cuò)誤，請(qǐng)及時(shí)提醒作者，避免錯(cuò)誤進(jìn)一步傳播。
2 ABAQUS材料庫使用
打開軟件后，進(jìn)入Property模塊，左側(cè)將增加ABAQUS材料庫使用界面，如下圖所示：
3 POLARIS基礎(chǔ)材料庫
3.1 POLARIS_MAT_BASE基礎(chǔ)材料庫下載
ABAQUS材料庫插件_POLARIS_MAT_BASE.zip
3.2 文件說明
壓縮包共包含兩個(gè)lib文件，分布為：POLARIS_MAT_BASE_SI_m.lib和POLARIS_MAT_BASE_SI_mm.lib，分別表示國際_米制（Kg-m-s）和（T-mm-s）兩種單位制的材料，相互之間的轉(zhuǎn)換關(guān)系請(qǐng)查看：《有限元的單位》。
展開 abaqus調(diào)用damask實(shí)現(xiàn)FCC，BCC，HCP多晶織構(gòu)演化和應(yīng)力應(yīng)變場分布模擬
FCC------以鋁為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
RD拉伸織構(gòu)：
RD壓縮織構(gòu)：
ND平面應(yīng)變壓縮織構(gòu)：
BCC------以鐵素體為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
拉伸織構(gòu)：
壓縮織構(gòu)：
平面應(yīng)變壓縮織構(gòu)：
HCP------以鎂為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
拉伸織構(gòu)：
壓縮織構(gòu)：
平面應(yīng)變壓縮織構(gòu)：
多晶局部應(yīng)力應(yīng)變場分布模擬與宏觀應(yīng)力應(yīng)變響應(yīng)。以FCC－鋁為例子。BCC與HCP同理。
展開 
abaqus調(diào)用damask實(shí)現(xiàn)FCC，BCC，HCP多晶織構(gòu)演化和應(yīng)力應(yīng)變場分布模擬
FCC------以鋁為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
RD拉伸織構(gòu)：
RD壓縮織構(gòu)：
ND平面應(yīng)變壓縮織構(gòu)：
BCC------以鐵素體為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
拉伸織構(gòu)：
壓縮織構(gòu)：
平面應(yīng)變壓縮織構(gòu)：
HCP------以鎂為代表，參數(shù)使用原始abaqus提供的參數(shù)
織構(gòu)演化模擬模型使用包含1000個(gè)單元的1*1*1mm立方體，其中每個(gè)單元表示一個(gè)特定取向的單晶，初始織構(gòu)使用軟件生成1000組隨機(jī)取向，并分配給不同的單元，模型和初始織構(gòu)如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進(jìn)行平面應(yīng)變壓縮時(shí)的織構(gòu)
拉伸織構(gòu)：
壓縮織構(gòu)：
平面應(yīng)變壓縮織構(gòu)：
多晶局部應(yīng)力應(yīng)變場分布模擬與宏觀應(yīng)力應(yīng)變響應(yīng)。以FCC－鋁為例子。BCC與HCP同理。
展開 Abaqus調(diào)用damask實(shí)現(xiàn)軋制變形中FCC，BCC織構(gòu)演化分析------案例六
Abaqus調(diào)用damask實(shí)現(xiàn)軋制變形中FCC，BCC織構(gòu)演化分析
案例實(shí)操一
1，使用abaqus建立20*20*20（mm）的立方塊
2，對(duì)立方塊進(jìn)行單元?jiǎng)澐止舶?000個(gè)單元
3，假設(shè)每個(gè)單元代表一個(gè)單獨(dú)的晶粒，通過腳本隨機(jī)賦予每個(gè)單元材料屬性
4，施加對(duì)應(yīng)的邊界提交（60%的下壓量）
5，提交與后處理材料數(shù)據(jù)
包含1000個(gè)晶粒的有限元模型
材料的初始取向分布
FCC軋制后的取向分布情況
BCC軋制后的取向分布情況
Abaqus中溫度輸出的規(guī)定 ￥10
Abaqus中溫度輸出的規(guī)定，詳細(xì)解釋了abaqus中實(shí)體單元，殼單元，梁單元溫度輸出的規(guī)定，并用實(shí)例進(jìn)行了展示。
基于abaqus溫度法多地層隧道開挖
1、自動(dòng)地應(yīng)力平衡方法在多地層模型中的實(shí)現(xiàn)方法；
2、溫度法以模擬真實(shí)工況下的軟化模量；
3、每次開挖前都進(jìn)行模量軟化
。
