次へ 上へ 前へ 目次へ
次へ:Heat convection (film.f) 上へ:User-defined loading 前へ:Distributed flux (dflux.f) 目次へ


分布荷重(dload.f)

このサブルーチンは 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:スケールする(デフォルト)
!



guido dhondt 2016-03-08