abaqus 載荷 順序的案例
ABAQUS順序熱力耦合分析實例
ABAQUS順序熱力耦合分析實例.pdf
Abaqus 熱-力順序耦合與 DFLUX 詳解 ￥59.9
——科研到工程：Abaqus Goldak 雙橢球 + FROM FILE 實現(xiàn)可復(fù)現(xiàn)實驗結(jié)果（含 Goldak 熱源 DFLUX ）
適用人群：做焊接/鍵合殘余應(yīng)力/變形預(yù)測、增材制造熱-力場分析的工程師與研究生
代碼環(huán)境：Abaqus/CAE 2019（Python 2.7），Abaqus/Standard（DFLUX Fortran 子程序）
本文提供 兩個腳本（Abaqus/CAE Python 自動網(wǎng)格建模腳本 + Fortran DFLUX 熱源子程序）梳理成一套可復(fù)用的 有限元計算流程：
從物理到實現(xiàn) 的清晰鏈路：能量輸入 → 傳熱 → 溫度–時間歷程 → FROM FILE 映射 → 彈塑性力學(xué)響應(yīng)；
建模與求解流程：幾何、分區(qū)、網(wǎng)格、邊界、步長、輸出與文件命名；
Goldak 雙橢球熱源與熱力耦合理論： 在 DFLUX 中的實現(xiàn)原理與關(guān)鍵參數(shù)；
目錄
- 用 Abaqus 做焊接/鍵合熱-力耦合的“一鍵批量建模與計算”
- 目錄
- 1. 為什么要做焊接/鍵合熱-力耦合？
- 2. Goldak 雙橢球熱源與能量守恒
- 3. 總體流程與工程目錄
- 4. 熱分析建模要點（Thermal）
- 5. 力學(xué)分析建模要點（Mechanical）
- 6. 自動化批量建模腳本（Python，最終版）
- 7. DFLUX：Goldak 體熱源子程序（Fortran）
- 8. 模型驗證
- 9. 參考參數(shù)與推薦文獻(xiàn)
1. 為什么要做焊接/鍵合熱-力耦合？
焊接/鍵合是強(qiáng)非線性、強(qiáng)非穩(wěn)態(tài)的多物理場過程：移動熱源瞬時把能量輸入到極小體積，熱擴(kuò)散與對流/輻射把能量帶走，材料在不同溫度區(qū)間內(nèi)經(jīng)歷彈性–塑性–循環(huán)硬化乃至回復(fù)。
展開 ABAQUS批量inp順序作業(yè)
work-more.txt
ABAQUS批量inp順序作業(yè).docx
ABAQUS 單元-節(jié)點排布順序解析（重點講解分析方法）
進(jìn)行ABAQUS UEL二次開發(fā)、或者研究界面問題的時候，比如cohsive單元界面問題，會涉及到單元-節(jié)點的排布順序。ABAQUS inp文件中的單元-節(jié)點順序排布很有講究，不能搞錯，這是因為節(jié)點的排布順序與內(nèi)部程序有關(guān)聯(lián)，內(nèi)部的程序我們不得而知，但是節(jié)點順序的規(guī)律可以從inp文件中看到，再對比cae節(jié)點中的節(jié)點編號可以總結(jié)規(guī)律。以下內(nèi)容介紹一些經(jīng)典的實體單元-節(jié)點排布順序。
這部分內(nèi)容其實挺好玩兒的！！！
（一）三角形單元
以一個10x10（無單位）的挖孔板為例，劃分的網(wǎng)格為：（注意到被單元為30、138和23包圍的三角形沒有編號，這是因為那個三角形不是單元，我在那里挖了個孔，因為背景顯示的問題，網(wǎng)格顏色和背景一樣，故沒有區(qū)分。作此說明）
這里我們首先查看節(jié)點的排布方式。分別以單元1、12、123為例。
編號為“1”的單元，在inp文件中，單元-節(jié)點的編號為：
1, 61, 1, 2
在cae界面中，我們只顯示編號為“1”的單元為：
發(fā)現(xiàn)節(jié)點排布的方式是“逆時針”順序。
編號為“12”的單元，在inp文件中，單元-節(jié)點的編號為：
12, 30, 9, 10
在cae界面中，我們只顯示編號為“12”的單元為：
發(fā)現(xiàn)節(jié)點排布的方式是“逆時針”順序。
編號為“123”的單元，在inp文件中，單元-節(jié)點的編號為：
123, 83, 34, 33
在cae界面中，我們只顯示編號為“12”的單元為：
發(fā)現(xiàn)節(jié)點排布的方式是“逆時針”順序。
至此我們可以得出結(jié)論，在一階三角形單元中，單元-節(jié)點的索引方式遵循“逆時針”順序，即按照節(jié)點順序彎曲右手四個手指，大拇指指向屏幕外。
展開 
Abaqus中平面應(yīng)力單元高斯積分點的順序
可以輸出umat接口中的變量coords進(jìn)行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結(jié)果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應(yīng)力單元高斯積分點的順序為：
展開 ABAQUS之熱應(yīng)力 順序耦合分析-終章
圖6 溫度及應(yīng)力分布云圖
文章來源：ABAQUS仿真世界
ABAQUS激光切割（熱力順序耦合DFLUX+VUSDFLD）仿真案例講解
abaqus熱力耦合---順序（間接）耦合和完全（直接）完全耦合的結(jié)果對比 ￥200
<p> 前言</p><p>使用abaqus分析熱力學(xué)計算的例子很多，但是并沒有見有人發(fā)過順序耦合和直接完全耦合的對比，而且網(wǎng)上關(guān)于熱力耦合分析的教程又很少，而相關(guān)書籍上一般都用預(yù)定義場分析均勻溫度場，但是對于梯度載荷需要用到順序耦合或完全耦合。
展開 基于ABAQUS海底滑坡模擬過程中海底載荷（泥線處海水對海底泥線的載荷）如何施加？ ￥3
利用ABAQUS進(jìn)行海底滑坡或海底沉降或滑坡模擬過程中，如果海底是水平的，則該載荷很容易添加，如果海底存在一定的坡度，則不同位置處海底載荷不相等，那么就需要利用一定的手段進(jìn)行施加。
本貼內(nèi)容就針對該問題為初學(xué)者進(jìn)行解惑。入門ABAQUS高級使用者請繞路
如果假設(shè)模型模擬參數(shù)如下：
①尺寸：長250m，深125m，最淺處水深200m
那么海底泥線處載荷如何施加呢？
abaqus順序熱力耦合odb文件提取溫度數(shù)據(jù)加入后提示錯誤，求大佬解答
***ERROR: System Error in Opening DDB file. Please report the problem to your
local office.
THE PROGRAM HAS DISCOVERED 1 FATAL ERRORS
** EXECUTION IS TERMINATED **
Abaqus移動載荷 ￥25
ABAQUS——DLOAD和VDLOAD子程序應(yīng)用（移動載荷隱式和顯示）
abaqus自定義載荷子程序------Dload使用 ￥29.9
abaqus子程序Dload的主要作用：
（1）可用于定義作為位置、時間、單元編號、被加載積分點數(shù)量等的函數(shù)分布載荷大小的變化。
（2）在應(yīng)力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調(diào)用；
（3）將在每個積分點調(diào)用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；
（4）不能在基于模態(tài)的程序中用于描述負(fù)載的時間變化；并且忽略可能與相關(guān)聯(lián)的階躍定義或非均勻分布負(fù)載定義一起出現(xiàn)的任何幅度參考。
子程序接口界面
SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
1 COORDS,JLTYP,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION TIME(2), COORDS (3)
CHARACTER*80 SNAME
user coding to define F
RETURN
END
待定義變量
F：分布載荷的大小。表面載荷的單位為FL?2，體力的單位為FL?3。F將作為基于單元或基于表面的分布式載荷定義的一部分指定的載荷大小傳遞到例程中。如果未定義大小，F(xiàn)將作為零傳入。對于使用修正Riks法（靜態(tài)應(yīng)力分析）的靜態(tài)分析，F(xiàn)必須定義為荷載比例系數(shù)λ的函數(shù)。分布式負(fù)載大小不可用于輸出目的。
用于傳遞信息的變量
KSTEP：Step 編號
KINC：增量數(shù)
TIME(1)：當(dāng)前分析步對應(yīng)的當(dāng)前時間
TIME(2)：所有分析步對應(yīng)的當(dāng)前時間
NOEL：單元編號
NPT：根據(jù)荷載類型，構(gòu)件內(nèi)或構(gòu)件表面上的荷載積分點編號。
展開 Abaqus 波浪載荷計算 Step by Step ￥3
Abaqus波浪載荷計算-01-15.pdf
ABAQUS中橢圓形移動載荷DLOAD和UTRACLOAD子程序詳解：從定義到實現(xiàn) ￥288
圖5 切向載荷分布
3、子程序DLOAD和UTRACLOAD編程實現(xiàn)
確定好法向和切向移動載荷分布后，利用FORTRAN編寫DLOAD和UTRACLOAD用戶子程序，實現(xiàn)法向和切向移動載荷施加，載荷施加在滾動接觸體表面（以圖6所示的鋼軌踏面施加移動載荷為例）。
圖6 鋼軌表面施加法向和切向移動載荷
3.1 法向移動載荷實現(xiàn)DLOAD
DLOAD子程序是ABAQUS中定義體載荷、面載荷、線載荷等的一種接口，通過Fortran代碼自定義每個積分點上的載荷值。DLOAD適用于定義在單元上的載荷（如壓力、密度效應(yīng)等）；可以利用時間（TIME(1)）、空間坐標(biāo)（COORD）、元素編號等信息，來定義移動的載荷區(qū)域或強(qiáng)度，主要用于法向載荷的模擬。
展開 Abaqus邊界和載荷op參數(shù)詳解
Abaqus中施加邊界和載荷時，op參數(shù)設(shè)置不同必然帶來結(jié)果差異，說明如下
a.載荷
Cload施加集中載荷時，op選項默認(rèn)設(shè)為mod，即保留當(dāng)前載荷（前一分析步定義載荷）并在此基礎(chǔ)上修改已有載荷，或增加新集中載荷（在未定義自由度上）；如op設(shè)為new，則所有當(dāng)前施加在模型上的載荷均被移除（移除時在當(dāng)前分析步內(nèi)線性減至零或突降至零，取決于分析步內(nèi)定義的Amplitude卡片。靜力分析中線性移除，動力分析則即刻移除），并重新定義新載荷。
以如下模型為例，第一個分析步在1號螺栓上施加500N力，第二個分析步在2號螺栓上施加1000N力。在第二個分析步中對1號螺栓外載分別按mod和new選項進(jìn)行定義。
圖1 載荷模型
分別輸出耦合面參考點的支反力和兩加載點的外載變化歷程。支反力和兩點外載變化曲線如下圖
圖2 支反力
圖3 外載變化
mod設(shè)定，1號螺栓外載在step 1中線性增加至指定的500N，在step 2 中保持500N定值；2號螺栓外載在step 1中未定義，在step 2開始線性加載至1000N。總的支反力在step 1結(jié)束時與1號螺栓外載平衡，step 2結(jié)束時與兩個螺栓外載平衡。
New設(shè)定，從step 2開始，1號螺栓外載在分析時段內(nèi)線性減少至零（該點未定義新載荷，默認(rèn)新值為零）；2號螺栓載荷因盡在step 2定義，兩條曲線完全重合。總的支反力與螺栓外載疊加曲線相平衡。
b.邊界
Boundary卡片，op為mod時（默認(rèn)值），更改已存邊界條件或為之前未被約束自由度添加邊界條件；op為new時，所有當(dāng)前生效的邊界條件均被移除。僅移除選定邊界條件時，使用new選項并重新指定所有需要保留的邊界條件。如果在standard的盈利/位移分析中移除一個邊界條件，它會被一個與前一分析步終止時該約束自由度處支反力相等的集中力代替。
展開 