投稿者 スレッド: 【要素関係】直交異方弾性体を20節点ソリッド要素で解析したときに出力される応力値が零になるバグ  (参照数 1144 回)

ghashimoto

  • Okulab
  • Newbie
  • *****
  • 投稿: 17
    • プロフィールを見る
直交異方弾性体を20節点ソリッド要素で解析したとき,
出力される応力値が零になるバグがあります.

原因は,直交異方性を決める局所座標系を取得できていないためです.

src/lib/static_LIB_3d.f90の
サブルーチンUpdateST_C3の
      IF( matlaniso ) THEN
        CALL set_localcoordsys( coords, g_LocalCoordSys(cdsys_ID), coordsys, serr )
        IF( serr == -1 ) STOP "Fail to setup local coordinate"
        IF( serr == -2 ) THEN
          WRITE(*, *) "WARNING! Cannot setup local coordinate, it is modified automatically"
        END IF
      END IF

      IF( cdsys_ID > 0 ) THEN
        CALL set_localcoordsys( coords, g_LocalCoordSys(cdsys_ID), coordsys, serr )
        IF( serr == -1 ) STOP "Fail to setup local coordinate"
        IF( serr == -2 ) THEN
          WRITE(*, *) "WARNING! Cannot setup local coordinate, it is modified automatically"
        END IF
      END IF
に修正する必要があります.
« 最終編集: 5月 29, 2015, 07:25:48 午後 by ghashimoto »