Dload子程序的案例
ABAQUS DLOAD (子程序3) 自定義荷載 (附Fortran和input) ￥10
該例子可用于校核 ABAQUS關聯Visual Studio和Fortran compiler時 該關聯是否成功；通過簡單的Fortran代碼理解用Fortran自定義荷載時，Fortran代碼與ABAQUS主程序的接口問題；該例子適合新手入門，如何用Fortran自定義荷載
所編DLOAD子程序
對應用于校核該DLOAD子程序的Abaqus inp文件
ABAQUS——DLOAD子程序應用（移動載荷） ￥9.99
DLOAD子程序可對上述過程進行有效模擬。
DLOAD子程序介紹：
在網上對于單移動載荷的實現辦法介紹較為全面，實現的效果如下：
同時附上我編寫的子程序主體片段，關鍵是需要對移動載荷生效區域和失效區域的準確描述（通過坐標），本例中，車輪與地面間的接觸區域簡化為長方形（寬度即為輪寬）：
但若要實現多個載荷在不同位置同時移動呢？這就需要花費一點心思，觀察子程序的可用參數：對COORDS或者SNAME的加強判斷即可實現多載荷移動，具體效果如下：
收費內容為上述模型的cae文件和上述兩個案例的子程序文件示例
基于Abaqus的DLOAD子程序實現移動載荷 ￥12
<p>有限元模型中，當載荷比較復雜時，難以通過ABAQUS/CAE界面直接進行設置，這時候就需要使用DLOAD子程序。在工程實際應用中，經常會遇到移動載荷的例子，如車轍實驗，汽車過橋等。本案例介紹基于ABAQUS子程序的DLOAD實現移動載荷。
Abaqus用Dload子程序實現移動載荷
今天就介紹一下用abaqus子程序Dload實現移動載荷——“CAE仿真實驗室”出品。
abaqus自定義載荷子程序------Dload使用 ￥29.9
abaqus子程序Dload的主要作用：
（1）可用于定義作為位置、時間、單元編號、被加載積分點數量等的函數分布載荷大小的變化。
（2）在應力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調用；
（3）將在每個積分點調用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；
（4）不能在基于模態的程序中用于描述負載的時間變化；并且忽略可能與相關聯的階躍定義或非均勻分布負載定義一起出現的任何幅度參考。
子程序接口界面
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將作為零傳入。對于使用修正Riks法（靜態應力分析）的靜態分析，F必須定義為荷載比例系數λ的函數。分布式負載大小不可用于輸出目的。
用于傳遞信息的變量
KSTEP：Step 編號
KINC：增量數
TIME(1)：當前分析步對應的當前時間
TIME(2)：所有分析步對應的當前時間
NOEL：單元編號
NPT：根據荷載類型，構件內或構件表面上的荷載積分點編號。
展開 abaqus自定義載荷子程序------Dload使用
abaqus子程序Dload的主要作用：
（1）可用于定義作為位置、時間、單元編號、被加載積分點數量等的函數分布載荷大小的變化。
（2）在應力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調用；
（3）將在每個積分點調用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；
（4）不能在基于模態的程序中用于描述負載的時間變化；并且忽略可能與相關聯的階躍定義或非均勻分布負載定義一起出現的任何幅度參考。
子程序接口界面
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將作為零傳入。對于使用修正Riks法（靜態應力分析）的靜態分析，F必須定義為荷載比例系數λ的函數。分布式負載大小不可用于輸出目的。
用于傳遞信息的變量
KSTEP：Step 編號
KINC：增量數
TIME(1)：當前分析步對應的當前時間
TIME(2)：所有分析步對應的當前時間
NOEL：單元編號
NPT：根據荷載類型，構件內或構件表面上的荷載積分點編號。
展開 有限元理論基礎及Abaqus內部實現方式研究系列36：DLOAD用戶子程序開發步驟
通用結構有限元軟件iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第36篇： DLOAD用戶子程序開發步驟==
前段時間很多朋友都問我們iSolver怎么實現DLOAD子程序，相比前面第八、九篇介紹的UMAT用戶自定義材料的開發，第二十、二十三篇介紹的UEL用戶自定義單元的開發，DLOAD要簡單很多，可以說，沒有任何復雜的算法，所以我們在此也僅僅介紹一下DLOAD的開發步驟，也順便給iSolver用戶提供一個簡單的算例幫助大家更好的使用iSolver的DLOAD功能。
實際生活中很多時候載荷是隨著某些物理量變化的，最常見的是隨著時間或者位置的變化，如果僅僅是隨著時間變化，那么在Abaqus中可以用Amplitude設置一個載荷和時間的對應表就行了；如果僅僅是隨著加載位置變化，那么可以直接把載荷位置獨立劃分出來，譬如設成一個單獨Set，在整個時間段內都只在該位置加載就行。那如果載荷既隨時間變化，同時在時間推進過程中加載位置也在變怎么辦呢？
譬如大壩的水位隨著潮汐的作用一天內會不同，導致不同時間段對大壩的水壓也不同。再譬如汽車在大橋上行駛，每一個時刻車子所在位置不同，對橋的荷載自然也不同了。
上述這些情況就必須要用到子程序了，如果是分布載荷，那么就是DLOAD子程序，DLOAD主要是將用戶特定的分布載荷隨時間、單元、積分點等變化算法編寫為計算機語言表示的公式，并實現和求解器之間的交互迭代。
展開 ABAQUS中橢圓形移動載荷DLOAD和UTRACLOAD子程序詳解：從定義到實現 ￥288
圖5 切向載荷分布
3、子程序DLOAD和UTRACLOAD編程實現
確定好法向和切向移動載荷分布后，利用FORTRAN編寫DLOAD和UTRACLOAD用戶子程序，實現法向和切向移動載荷施加，載荷施加在滾動接觸體表面（以圖6所示的鋼軌踏面施加移動載荷為例）。
圖6 鋼軌表面施加法向和切向移動載荷
3.1 法向移動載荷實現DLOAD
DLOAD子程序是ABAQUS中定義體載荷、面載荷、線載荷等的一種接口，通過Fortran代碼自定義每個積分點上的載荷值。DLOAD適用于定義在單元上的載荷（如壓力、密度效應等）；可以利用時間（TIME(1)）、空間坐標（COORD）、元素編號等信息，來定義移動的載荷區域或強度，主要用于法向載荷的模擬。
展開 ABAQUS二次開發（DLOAD實現車輪移動荷載）
算例采用DLOAD子程序實現，涉及到DLOAD子程序編寫。
（）幾何模型信息
模型為四輪車，每個輪子與地面的接觸區域為0.1x0.3（不考慮實際情況），車輛在平板居中放置，從右邊界啟動，向左移動，示意圖如下：
（）計算模型
網格整體如下圖：
網格橫斷面如下圖：
車輪經過的區域網格設置加密，車輛中間區域與外部區域網格尺寸設置漸變。其中，車輪經過區域的網格如下圖：
網格信息詳見附件。
（）效果
不想費勁寫介紹了，直接看視頻，有需要的就直奔附件吧！
！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
！視頻聲音忘記去除了，是個很吵的DJ，請提前關閉聲音觀看！！
！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
（）附件
ABAQUS
JOB-1.INP
BRIDGE.CAE
BRIDGE.JNL
DYNAMIC_LOAD.FOR
RUN.BAT
ABAQUS.zip
展開 VDLOAD與VUAMP聯合使用——實現變位置變幅值復雜加載 ￥15
測試案例展示效果如下：
ABAQUS——DLOAD子程序應用（移動載荷）一文，我們討論了DLOAD子程序在車橋通過問題的應用。
https://www.yqgqt.org.cn/content/post/1257703
在本例中，我們按照同樣的思路，在VDLOAD中實現了相同的過程，VDLOAD中載荷的移動速度和時間具有真實的物理含義，但如何實現變化的移動載荷呢？
由幫助文檔可知，添加 “*AMPLITUDE” 即可在VDLOAD中傳遞在CAE中定義的幅值，該幅值通過TAB表實現，并通過線性插值擬合。
在CAE界面定義好VDLOAD子程序和幅值表后，需要手動修改inp文件，將VDLOAD和幅值相對應，便可將表中的幅值傳遞給子程序，定義格式如下。
采用上述方法具有一定的局限性，比如遇到復雜周期問題、幅值需要速度位移參與計算、依賴于上步幅值解、復雜函數（三角函數疊加）等情況變難以實現。
周期幅值：隨機幅值
實現隨機幅值的加載僅需調用隨機數函數：random_number()。
故本文介紹另一種解決辦法：將VUAMP子程序（幅值自定義子程序）和VDLOAD子程序聯合使用，實現位置可變、幅值自定義的復雜壓力加載過程。
測試案例中，一邊采用VUAMP子程序方式定義幅值，另一邊采用幅值表進行定義，均設置為隨時間的線性遞增加載。
得到的仿真效果如下，左右兩側效果相同，驗證了上述思路的有效性。
為進一步驗證實現方案（VUAMP+VDLOAD），兩個幅值不同時的加載效果當一端遞增，一端為三角函數周期幅值時
同一個面下的不同幅值的移動載荷
展開 基于ABAQUS單點顯式VDLOAD/隱式DLOAD激光沖擊加載(圓形光斑和方形光斑) ￥50
ABAQUS軟件可以通過顯式VDLOAD或隱式DLOAD子程序二次開發進行單點/多點激光沖擊模擬，效率高。不同之處在于隱式相對計算時間長但是可以在第一步沖擊后直接在第二步進行回彈分析，無需進行數據傳遞來計算回彈過程。
本帖基于ABAQUS的VDLOAD/DLOAD子程序對比顯式/隱式算法下不同光斑形狀的應力和塑性應變。首先進行模型構建。
材料采用Ti-6Al-4V鈦合金，有限單元區賦予彈塑性參數，無限單元賦予彈性參數。具體參數如下：
密度：4.5e-9；彈性：1.2e5,0.34；塑性：A:1098 B:1092 C:0.014 n:0.93 參考應變率：1
裝配：全局坐標原點與有限單元頂點重合。
分析步創建：隱式/顯式分析步，增量采用固定增量步
接口設定：
網格劃分，無限單元部分定義網格節點方向排布
建立模型輸出inp文件，將無限單元部分的單元類型改為CIN3D8
如果沖擊光斑為圓形光斑，網格細化至50微米，如果為方向光斑，網格100微米
圓形光斑在空間表現為高斯分布，表現為從中心區域到邊緣沖擊載荷逐漸變小。
方向光斑在空間表現為均布載荷，其峰值載荷為圓形光斑的0.618倍，一些研究表明相同激光參數下方形光斑搭接沖擊材料疲勞性能較高。
展開 
Dload循環荷載
大佬們，請教一個問題，模型z方向是80m，循環四次 ，分析步時長也改了4s，用DLoad子程序跑出來怎么只循環一次在3.5s到4s。這個是子程序的代碼，請問是不是do這里沒寫對，應該怎么改啊。
某移動罩下軌道梁，在移動罩運動時，產生較大變形，通過有限元分析，使用動載荷分析 ￥20
有限元瞬態分析步驟：
幾何建模：細化載荷移動路徑網格（尺寸≤1/10波長）；
接觸定義：采用面-面接觸模擬輪軌/車橋相互作用；
載荷施加：通過APDL命令流或用戶子程序實現移動載荷；
求解設置：時間步長滿足 Δt≤Tmin?/10?為最小振動周期）。
將各載荷添加于模型，其中移動罩載荷使用ABAQUS中DLOAD子程序實現，如圖1所示。
（a）高軌軌道梁尺寸
（b）高軌軌道梁模型及載荷
abaqus凍土路基的溫度-水分-變形多場耦合分析
本文中，溫度場分析通過film子程序和dflux子程序定義溫度邊界，通過hetval子程序定義相變熱。變形場分析通過dload子程序定義車輛載荷，通過uexpan子程序引入凍脹影響。
abaqus凍土路基的溫度-水分-變形多場耦合分析
本文中，溫度場分析通過film子程序和dflux子程序定義溫度邊界，通過hetval子程序定義相變熱。變形場分析通過dload子程序定義車輛載荷，通過uexpan子程序引入凍脹影響。
原文鏈接：
https://www.yqgqt.org.cn/content/post/1822630
歡迎大家前去投票打call.
最后，有需要歡迎通過微信公眾號聯系我們。
