subroutine fem3d(...) |
10 call addmon(...) |
if(numelh.ne0) call solide |
subroutine solide |
do 20 i=1,nelt |
do 20 j=1,nelg |
call unpki |
call strain |
call force |
20 continue |
if(...) return |
goto 10 |
subroutine force(...) |
do 10 i=lft,llt |
sgv(i) = sig1(i)-qp(i)*vol(i) |
10 continue |
do 50 n=1,nnc |
i0=ia(n) |
i1=ia(n+1)-1 |
do 50 i=i0,i1 |
e(1,ix(i))=e(1,ix(i))+ep11(i) |
50 continue |