Updated Lagrange法で応力積分式のトレス項の符号
Updated Lagrange法における応力積分のサブルーチン(fistr1/src/lib/static_LIB_3d.f90)Update_Stress3D(…)の653行目に右辺の最後項dumstress*trDの符号は'-'はずではないでしょうか?
dum(:,:) = dumstress + matmul( rot,dumstress ) - matmul( dumstress, rot ) + dumstress*trD
添付のリンク先(http://www.cee.northwestern.edu/people/bazant/PDFs/Papers/522.pdf)参考文献の式(41)より、右辺後ろの3項を左辺に移動させると、ちょうど上のコードの右辺後ろの3項に該当しますが、最後項の符号は'-'になるべきと思います。 ご確認いただければ幸いです。