このサブルーチンは PxNUy 形式の分布荷重ラベルで特徴づけられる非一様な圧力のために使用されます。*DLOADを参照してください。荷重ラベルの上限は20文字です。特に y は非一様なフラックスのパターンを判別するために使用できます。ヘッダーと変数の意味は以下の通りです。
subroutine dload(f,kstep,kinc,time,noel,npt,layer,kspt,coords,jltyp,loadtype,vold,co,lakonl,konl,ipompc,nodempc,coefmpc,nmpc,ikmpc,ilmpc,iscale,veold,rho,amat,mi)
!
! ユーザーサブルーチン dload
!
!
! 入力：
!
! kstep ステップ番号
! kinc インクリメント番号
! time(1) 現在のステップ時間
! time(2) 現在の総時間
! noel 要素番号
! npt 積分点番号
! layer 現在は不使用
! kspt 現在は不使用
! coords(1..3) 積分点のグローバル座標
! jltyp 荷重面コード：
! 21 = 面1
! 22 = 面2
! 23 = 面3
! 24 = 面4
! 25 = 面5
! 26 = 面6
! loadtype 荷重タイプラベル
! vold(0..4,1..nk) 全節点で解かれる場
! 0: 温度
! 1: グローバル X 方向の変位
! 2: グローバル Y 方向の変位
! 3: グローバル Z 方向の変位
! 4: 静圧
! veold(0..3,1..nk) 全節点での解かれる場の時間微分
! 0: 温度変化
! 1: グローバル X 方向の速さ
! 2: グローバル Y 方向の速さ
! 3: グローバル Z 方向の速さ
! co(3,1..nk) 全節点の座標
! 1: グローバル X 方向の座標
! 2: グローバル Y 方向の座標
! 3: グローバル Z 方向の座標
! lakonl 要素ラベル
! konl(1..20) 要素に属する節点
! ipompc(1..nmpc) ipompc(i) は場 nodempc の MPC i の最初の項を指す
! nodempc(1,*) MPC の項の節点番号
! nodempc(2,*) MPC の項の座標方向
! nodempc(3,*) 0以外の場合：場 nodempc の MPC の次の項を指す
! 0の場合：MPC 定義終了
! coefmpc(*) MPC の項の係数
! nmpc MPC の数
! ikmpc(1..nmpc) 順序付けられたMPCのグローバルな自由度
! グローバルな自由度は MPC の依存項の 8*(node-1)+direction（direction = 0：温度, 1-3：変位, 4：静圧, 5-7：回転）
! ilmpc(1..nmpc) ilmpc(i) は ikmpc(i) での参照番号に対応する MPC 番号
! rho ローカルな密度
! amat 材料名
! mi(1) 要素あたりの積分点の最大数（全要素にわたる最大値）
! mi(2) 場での節点あたりの最大自由度（全節点にわたる最大値）。例 v(0:mi(2))。
!
! 出力：
!
! f 分布荷重の大きさ
! iscale 静計算でのステップ時間より小さなインクリメントでフラックスをスケールするかどうか
! 0：スケールしない
! 1：スケールする（デフォルト）
!