abaqus定義線載荷的案例
基于ABAQUS海底滑坡模擬過程中海底載荷（泥線處海水對海底泥線的載荷）如何施加？ ￥3
利用ABAQUS進行海底滑坡或海底沉降或滑坡模擬過程中，如果海底是水平的，則該載荷很容易添加，如果海底存在一定的坡度，則不同位置處海底載荷不相等，那么就需要利用一定的手段進行施加。
本貼內(nèi)容就針對該問題為初學者進行解惑。入門ABAQUS高級使用者請繞路
如果假設(shè)模型模擬參數(shù)如下：
①尺寸：長250m，深125m，最淺處水深200m
那么海底泥線處載荷如何施加呢？
abaqus自定義載荷子程序------Dload使用
abaqus子程序Dload的主要作用：
（1）可用于定義作為位置、時間、單元編號、被加載積分點數(shù)量等的函數(shù)分布載荷大小的變化。
（2）在應力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調(diào)用；
（3）將在每個積分點調(diào)用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；
（4）不能在基于模態(tài)的程序中用于描述負載的時間變化；并且忽略可能與相關(guān)聯(lián)的階躍定義或非均勻分布負載定義一起出現(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)應力分析）的靜態(tài)分析，F(xiàn)必須定義為荷載比例系數(shù)λ的函數(shù)。分布式負載大小不可用于輸出目的。
用于傳遞信息的變量
KSTEP：Step 編號
KINC：增量數(shù)
TIME(1)：當前分析步對應的當前時間
TIME(2)：所有分析步對應的當前時間
NOEL：單元編號
NPT：根據(jù)荷載類型，構(gòu)件內(nèi)或構(gòu)件表面上的荷載積分點編號。
展開 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適用于定義在單元上的載荷（如壓力、密度效應等）；可以利用時間（TIME(1)）、空間坐標（COORD）、元素編號等信息，來定義移動的載荷區(qū)域或強度，主要用于法向載荷的模擬。
展開 Abaqus/CAE中定義與深度相關(guān)的載荷
在使用有限元分析時，有時候會遇到與深度相關(guān)的載荷，如：靜水壓力、地應力等。Abaqus/CAE提供了直接定義這兩種載荷的方式，下面通過簡單的例子介紹如何操作以及使用中應該注意到問題。
靜水壓力載荷的定義
假設(shè)高100M的壓力容器，需在其內(nèi)部定義靜水壓力，如下圖1：
圖1 高100M的壓力容器內(nèi)部施加靜水壓力
靜水壓力在Pressure下定義，分布形式選擇Hydrostatic。根據(jù)壓力的計算公式：p=ρgh，得知最深出的水壓為980000Pa。還需要定義零水位高度坐標值（Zero pressure height）和水位最深高度坐標值（Reference pressure height），這里分別輸入100和0。如下圖2所示。
圖2 靜水壓力的施加
進行簡單的計算后，發(fā)現(xiàn)壓力云圖如下，靜水壓力分布和預期的完全不一樣。
圖3 靜水壓力云圖
問題出在什么地方呢？這是因為Abaqus默認深度在Z軸方向上，而我們在建模時，將深度方向放在了Y軸方向，如圖1。
將深度方向放在Z軸方向上后重新計算，得到的靜水壓力分布如下：
圖4 深度方向放在Z軸（左）以及靜水壓力云圖（右）
注意到筒底的靜水壓力接近于980000Pa，筒頂?shù)撵o水壓力也不是0Pa。這是因為云圖中的數(shù)值來源于單元的積分點，而不是節(jié)點。
靜水壓力也能用于2D模型，只是需要注意到是，此時，深度方向在Y軸方向。
地應力的施加
Abaqus/CAE中施加地應力時，需要在預定義場中定義，如下圖5所示，載荷類型為Geostatic。
圖5 地應力的施加
需要定義的參數(shù)：第一點的深度坐標值與對應的垂直應力；第二點的深度坐標值與對應的垂直應力；兩方向的側(cè)向土壓力系數(shù)。
展開 
abaqus自定義載荷子程序------Dload使用 ￥29.9
abaqus子程序Dload的主要作用：
（1）可用于定義作為位置、時間、單元編號、被加載積分點數(shù)量等的函數(shù)分布載荷大小的變化。
（2）在應力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調(diào)用；
（3）將在每個積分點調(diào)用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；
（4）不能在基于模態(tài)的程序中用于描述負載的時間變化；并且忽略可能與相關(guān)聯(lián)的階躍定義或非均勻分布負載定義一起出現(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)應力分析）的靜態(tài)分析，F(xiàn)必須定義為荷載比例系數(shù)λ的函數(shù)。分布式負載大小不可用于輸出目的。
用于傳遞信息的變量
KSTEP：Step 編號
KINC：增量數(shù)
TIME(1)：當前分析步對應的當前時間
TIME(2)：所有分析步對應的當前時間
NOEL：單元編號
NPT：根據(jù)荷載類型，構(gòu)件內(nèi)或構(gòu)件表面上的荷載積分點編號。
展開 