36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "mvsiz_p.inc"
44
45
46
47
48
49
50 INTEGER IXS(NIXS,*),NEL ,IDDL(*) ,NDOF(*) ,
51 . IADK(*) ,JDIK(*)
52
54 . k_diag(*) ,k_lt(*) ,
55 . k11(9,*) ,k12(9,*) ,k13(9,*) ,k14(9,*) ,k22(9,*) ,
56 . k23(9,*) ,k24(9,*) ,k33(9,*) ,k34(9,*) ,k44(9,*),
57 . off(*)
58
59
60
61 INTEGER EP,NN(MVSIZ,4)
62
63 DO ep = 1,nel
64 nn(ep,1)=ixs(2,ep)
65 nn(ep,2)=ixs(4,ep)
66 nn(ep,3)=ixs(7,ep)
67 nn(ep,4)=ixs(6,ep)
68 ENDDO
69
70 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
71 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
72 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
73 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
74
75 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
76 . k_diag,k_lt,k12,3,off)
77 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
78 . k_diag,k_lt,k13,3,off)
79 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
80 . k_diag,k_lt,k14,3,off)
81 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
82 . k_diag,k_lt,k23,3,off)
83 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
84 . k_diag,k_lt,k24,3,off)
85 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
86 . k_diag,k_lt,k34,3,off)
87
88
89 RETURN
subroutine assem_kii(ni, nel, iddl, iadk, k_diag, k_lt, kii, nd, off)
subroutine assem_kij(ni, nj, nel, iddl, iadk, jdik, k_diag, k_lt, kij, nd, off)