38
39
40
41#include "implicit_f.inc"
42
43
44
45#include "mvsiz_p.inc"
46
47
48
49
50
51
52 INTEGER IXS(NIXS,*),NEL ,IDDL(*) ,NDOF(*) ,
53 . IADK(*) ,JDIK(*)
54
56 . k_diag(*) ,k_lt(*) ,
57 . k11(9,*),k12(9,*),k13(9,*),k14(9,*),
58 . k15(9,*),k16(9,*),k22(9,*),k23(9,*),
59 . k24(9,*),k25(9,*),k26(9,*),k33(9,*),
60 . k34(9,*),k35(9,*),k36(9,*),k44(9,*),
61 . k45(9,*),k46(9,*),k55(9,*),k56(9,*),
62 . k66(9,*),off(*)
63
64
65
66 INTEGER EP,NN(MVSIZ,6)
67
68 DO ep = 1,nel
69 nn(ep,1)=ixs(2,ep)
70 nn(ep,2)=ixs(3,ep)
71 nn(ep,3)=ixs(4,ep)
72 nn(ep,4)=ixs(6,ep)
73 nn(ep,5)=ixs(7,ep)
74 nn(ep,6)=ixs(8,ep)
75 ENDDO
76
77 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
78 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
79 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
80 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
81 CALL assem_kii(nn(1,5),nel,iddl,iadk,k_diag,k_lt ,k55,3,off)
82 CALL assem_kii(nn(1,6),nel,iddl,iadk,k_diag,k_lt ,k66,3,off)
83
84 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
85 . k_diag,k_lt,k12,3,off)
86 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
87 . k_diag,k_lt,k13,3,off)
88 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
89 . k_diag,k_lt,k14,3,off)
91 . k_diag,k_lt,k15,3,off)
92 CALL assem_kij(nn(1,1),nn(1,6),nel,iddl,iadk,jdik,
93 . k_diag,k_lt,k16,3,off)
94 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
95 . k_diag,k_lt,k23,3,off)
96 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
97 . k_diag,k_lt,k24,3,off)
99 . k_diag,k_lt,k25,3,off)
100 CALL assem_kij(nn(1,2),nn(1,6),nel,iddl,iadk,jdik,
101 . k_diag,k_lt,k26,3,off)
102 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
103 . k_diag,k_lt,k34,3,off)
104 CALL assem_kij(nn(1,3),nn(1,5),nel,iddl,iadk,jdik,
105 . k_diag,k_lt,k35,3,off)
106 CALL assem_kij(nn(1,3),nn(1,6),nel,iddl,iadk,jdik,
107 . k_diag,k_lt,k36,3,off)
108 CALL assem_kij(nn(1,4),nn(1,5),nel,iddl,iadk,jdik,
109 . k_diag,k_lt,k45,3,off)
110 CALL assem_kij(nn(1,4),nn(1,6),nel,iddl,iadk,jdik,
111 . k_diag,k_lt,k46,3,off)
112 CALL assem_kij(nn(1,5),nn(1,6),nel,iddl,iadk,jdik,
113 . k_diag,k_lt,k56,3,off)
114
115
116 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)