分布フラックス(dflux.f)
次へ 上へ 前へ 目次へ
次へ:Distribruted load (dload.f) 上へ:User-defined loading 前へ:Concentrated load (cload.f) 目次へ


分布フラックス(dflux.f)

このサブルーチンは SxNUy 形式の分布荷重ラベルで特徴づけられる非一様な熱フラックスのために使用されます。*DFLUXを参照してください。荷重ラベルの上限は20文字です。特に y は非一様なフラックスのパターンを判別するために使用できます。ヘッダーと変数の意味は以下の通りです。

      subroutine dflux(flux,sol,kstep,kinc,time,noel,npt,coords,jltyp,temp,press,loadtype,area,vold,co,lakonl,konl,ipompc,nodempc,coefmpc,nmpc,ikmpc,ilmpc,iscale,mi)
!
!     ユーザーサブルーチン dflux
!
!
!     入力:
!
!     sol                現在の温度値
!     kstep              ステップ番号
!     kinc               インクリメント番号
!     time(1)            現在のステップ時間
!     time(2)            現在の総時間
!     noel               要素番号
!     npt                積分点番号
!     coords(1..3)       積分点のグローバル座標
!     jltyp              荷重面コード:
!                        1  = 体積フラックス
!                        11 = 面1 
!                        12 = 面2 
!                        13 = 面3 
!                        14 = 面4 
!                        15 = 面5 
!                        16 = 面6
!     temp               現在は不使用
!     press              現在は不使用
!     loadtype           荷重タイプラベル
!     area               表面フラックスの場合:積分点によってカバーされる面積
!                        体積フラックスの場合:積分点によってカバーされる体積
!     vold(0..4,1..nk)   全節点で解かれる場
!                        0: 温度
!                        1: グローバル X 方向の変位
!                        2: グローバル Y 方向の変位
!                        3: グローバル Z 方向の変位
!                        4: 静圧
!     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 番号
!     mi(1)              要素あたりの積分点の最大数(全要素にわたる最大値)
!     mi(2)              場での節点あたりの最大自由度(全節点にわたる最大値)。例 v(0:mi(2))。
!
!     出力:
!
!     flux(1)            フラックスの大きさ
!     flux(2)            不使用:値を代入しないでください
!     iscale             静計算でのステップ時間より小さなインクリメントでフラックスをスケールするかどうか
!                        0:スケールしない
!                        1:スケールする(デフォルト)
!



guido dhondt 2016-03-08