1 | subroutine fem3d(...) |
2 | 10 call addmon(...) |
3 | if(numelh.ne0) call solide |
4 | subroutine solide |
5 | do 20 i=1,nelt |
6 | do 20 j=1,nelg |
7 | call unpki |
8 | call strain |
9 | call force |
10 | 20 continue |
11 | if(...) return |
12 | goto 10 |
13 | subroutine force(...) |
14 | do 10 i=lft,llt |
15 | sgv(i) = sig1(i)-qp(i)*vol(i) |
16 | 10 continue |
17 | do 50 n=1,nnc |
18 | i0=ia(n) |
19 | i1=ia(n+1)-1 |
20 | do 50 i=i0,i1 |
21 | e(1,ix(i))=e(1,ix(i))+ep11(i) |
22 | 50 continue |