次へ 上へ 前へ 目次へ
次へ:Hardening (uhardening.f) 上へ:User subroutines. 前へ:User subroutines. 目次へ


クリープ(creep.f)

ユーザーサブルーチン「creep.f」を使用すると、入力デッキでキーワードシーケンス *CREEP,LAW=USER を選択することでユーザーは独自のクリープ則を組み込むことができるようになります。入力/出力は材料の種類に依存します。材料の弾性プロパティーが等方性の場合にはフォン・ミーゼス応力が入力され、等価偏差クリープひずみ増分とそれをフォン・ミーゼス応力で微分したものが対応物として出力されます。材料の弾性プロパティーが異方性の場合には等価偏差クリープひずみ増分が入力され、フォン・ミーゼス応力、等価偏差クリープひずみ増分をフォン・ミーゼス応力で微分したものが出力されます。ただしクリープ状態は常に等方性です。弾性状態が等方性となるか異方性となるかは変数の値によります。ヘッダーと入出力変数の意味は以下の通りです。

      subroutine creep(decra,deswa,statev,serd,ec,esw,p,qtild,temp,dtemp,predef,dpred,time,dtime,cmname,leximp,lend,coords,nstatv,noel,npt,layer,kspt,kstep,kinc)
!
!     ユーザークリープルーチン
!
!     入力(一般):
!
!     statev(1..nstatv)  内部変数
!     serd               不使用
!     ec(1)              インクリメント開始時の等価クリープひずみ
!     ec(2)              不使用
!     esw(1..2)          不使用
!     p                  不使用
!     temp               インクリメント終了時の温度
!     dtemp              不使用
!     predef             不使用
!     dpred              不使用
!     time(1)            インクリメント終了時のステップ時間の値
!     time(2)            インクリメント終了時の総時間の値
!     dtime              時間増分
!     cmname             材料名
!     leximp             不使用
!     lend               2の場合: 等方性クリープ
!                        3の場合: 異方性クリープ
!     coords(1..3)       現在の積分点の座標
!     nstatv             内部変数番号
!     noel               要素番号
!     npt                積分点番号
!     layer              不使用
!     kspt               不使用
!     kstep              不使用
!     kinc               不使用
!
!    入力(等方性弾性材料のみ):
!     qtild              フォン・ミーゼス応力
!
!    入力(異方性弾性材料のみ):
!     decra(1)           等価偏差クリープひずみ増分
!
!
!     出力(一般):
!
!     decra(1)           等価偏差クリープひずみ増分
!     decra(2..4)        不使用
!     decra(5)           等価偏差クリープひずみ増分をフォン・ミーゼス応力で微分したもの
!     deswa(1..5)        不使用
!
!     出力(等方性弾性材料のみ):
!     decra(1)           等価偏差クリープひずみ増分
!
!     出力(異方性弾性材料のみ):
!     qtild              フォン・ミーゼス応力
!



guido dhondt 2016-03-08