styファイルを用いた初期値設定に関して
（下記についてすでに英語でも投稿したのですが、未解決のままなので日本語でも投稿させていただきます）
現在、事前計算から得られる*.styファイルを用いて、四面体要素に初期応力（initial stress）を、ノードに初速度（initial velocity）を与えて計算を行うことを試みています。使用しているソルバーはRadioss、モデル作成にはHyperMeshを用いています。
*.styファイル内の結果は、事前計算のEngine file内で以下のKeywordを記述して出力しています。
/OUTP/BRIC/STRES/FULL
/OUTP/VECT/VEL
事前計算で得られた初期値が、本計算で正しく読み込まれているかを確認するために、
各計算で出力される*.styファイル内の応力値および速度を比較しました。
その結果、以下に示すように同じ四面体要素で異なる応力テンソルの値が格納されていることがわかりました。
[事前計算]
/SOLID /TENSOR /STR_FUL
Full stress tensor + plastic strain
#FORMAT:(NPT, ISOLNOD (2I10/2E20.13), EINT(I),RHO(I),,I=1,NUMSOL
#FORMAT: (1P6E20.13/E20.13) ((TX(I,J),TY(I,J),TZ(I,J),TXY(I,J),TYZ(I,J),TZX(I,J),
#EPSP(J,I),J=1,NPT),I=1,NUMSOL,NPT=1)
1 4
-6.7621997110738E-07 1.4999960395455E-03
7.5204702163451E-02-3.8801280336932E-02 1.1782171715919E-02 3.1060056356859E-02-1.6245572187241E-02-2.6996609790200E-02
0.0000000000000E+00
[本計算]
/SOLID /TENSOR /STR_FUL
Full stress tensor + plastic strain
#FORMAT:(NPT, ISOLNOD (2I10/2E20.13), EINT(I),RHO(I),,I=1,NUMSOL
#FORMAT: (1P6E20.13/E20.13) ((TX(I,J),TY(I,J),TZ(I,J),TXY(I,J),TYZ(I,J),TZX(I,J),
#EPSP(J,I),J=1,NPT),I=1,NUMSOL,NPT=1)
1 4
-6.7621997110738E-07 1.4999960395455E-03
-2.8076925403528E-07-2.1057678031223E-17 2.8076925396508E-07 0.0000000000000E+00 0.0000000000000E+00 0.0000000000000E+00
0.0000000000000E+00
一方で、初速度は正しく反映されていました。
すでにいただいたアドバイスで、*.styファイルではなく*.staファイルを用いることを推奨されましたが、
上記と同様に、事前計算の応力テンソルの値が、本計算で正しく反映されていませんでした。
stress initializationに関する投稿もいくつか確認させていただきましたが、
いまだに解決策が見つけられていません。
なにかアドバイスがありましたら教えていただけたら幸いです。
Answers
-
本件、弊社技術サポートにて対応済みです。
0
Altair Employee