43
44
45
46#include "implicit_f.inc"
47
48
49
50#include "mvsiz_p.inc"
51
52
53
54
55
56
57 INTEGER IXS(NIXS,*),NEL ,IDDL(*) ,NDOF(*) ,
58 . IADK(*) ,JDIK(*)
59
61 . k_diag(*) ,k_lt(*) ,
62 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*),
63 . k15(3,3,*),k16(3,3,*),k17(3,3,*),k18(3,3,*),
64 . k22(3,3,*),k23(3,3,*),k24(3,3,*),k25(3,3,*),
65 . k26(3,3,*),k27(3,3,*),k28(3,3,*),
66 . k33(3,3,*),k34(3,3,*),k35(3,3,*),k36(3,3,*),
67 . k37(3,3,*),k38(3,3,*),k44(3,3,*),k45(3,3,*),
68 . k46(3,3,*),k47(3,3,*),k48(3,3,*),
69 . k55(3,3,*),k56(3,3,*),k57(3,3,*),k58(3,3,*),
70 . k66(3,3,*),k67(3,3,*),k68(3,3,*),
71 . k77(3,3,*),k78(3,3,*),k88(3,3,*),off(*)
72
73
74
75 INTEGER K,EP,NN(MVSIZ,8)
76
77 DO k = 1,8
78 DO ep = 1,nel
79 nn(ep,k)=ixs(k+1,ep)
80 ENDDO
81 ENDDO
82
83 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
84 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
85 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
86 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
87 CALL assem_kii(nn(1,5),nel,iddl,iadk,k_diag,k_lt ,k55,3,off)
88 CALL assem_kii(nn(1,6),nel,iddl,iadk,k_diag,k_lt ,k66,3,off)
89 CALL assem_kii(nn(1,7),nel,iddl,iadk,k_diag,k_lt ,k77,3,off)
90 CALL assem_kii(nn(1,8),nel,iddl,iadk,k_diag,k_lt ,k88,3,off)
91
92 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
93 . k_diag,k_lt,k12,3,off)
94 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
95 . k_diag,k_lt,k13,3,off)
96 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
97 . k_diag,k_lt,k14,3,off)
98 CALL assem_kij(nn(1,1),nn(1,5),nel,iddl,iadk,jdik,
99 . k_diag,k_lt,k15,3,off)
100 CALL assem_kij(nn(1,1),nn(1,6),nel,iddl,iadk,jdik,
101 . k_diag,k_lt,k16,3,off)
102 CALL assem_kij(nn(1,1),nn(1,7),nel,iddl,iadk,jdik,
103 . k_diag,k_lt,k17,3,off)
104 CALL assem_kij(nn(1,1),nn(1,8),nel,iddl,iadk,jdik,
105 . k_diag,k_lt,k18,3,off)
106 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
107 . k_diag,k_lt,k23,3,off)
108 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
109 . k_diag,k_lt,k24,3,off)
110 CALL assem_kij(nn(1,2),nn(1,5),nel,iddl,iadk,jdik,
111 . k_diag,k_lt,k25,3,off)
112 CALL assem_kij(nn(1,2),nn(1,6),nel,iddl,iadk,jdik,
113 . k_diag,k_lt,k26,3,off)
114 CALL assem_kij(nn(1,2),nn(1,7),nel,iddl,iadk,jdik,
115 . k_diag,k_lt,k27,3,off)
116 CALL assem_kij(nn(1,2),nn(1,8),nel,iddl,iadk,jdik,
117 . k_diag,k_lt,k28,3,off)
118 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
119 . k_diag,k_lt,k34,3,off)
120 CALL assem_kij(nn(1,3),nn(1,5),nel,iddl,iadk,jdik,
121 . k_diag,k_lt,k35,3,off)
122 CALL assem_kij(nn(1,3),nn(1,6),nel,iddl,iadk,jdik,
123 . k_diag,k_lt,k36,3,off)
124 CALL assem_kij(nn(1,3),nn(1,7),nel,iddl,iadk,jdik,
125 . k_diag,k_lt,k37,3,off)
126 CALL assem_kij(nn(1,3),nn(1,8),nel,iddl,iadk,jdik,
127 . k_diag,k_lt,k38,3,off)
128 CALL assem_kij(nn(1,4),nn(1,5),nel,iddl,iadk,jdik,
129 . k_diag,k_lt,k45,3,off)
130 CALL assem_kij(nn(1,4),nn(1,6),nel,iddl,iadk,jdik,
131 . k_diag,k_lt,k46,3,off)
132 CALL assem_kij(nn(1,4),nn(1,7),nel,iddl,iadk,jdik,
133 . k_diag,k_lt,k47,3,off)
134 CALL assem_kij(nn(1,4),nn(1,8),nel,iddl,iadk,jdik,
135 . k_diag,k_lt,k48,3,off)
136 CALL assem_kij(nn(1,5),nn(1,6),nel,iddl,iadk,jdik,
137 . k_diag,k_lt,k56,3,off)
138 CALL assem_kij(nn(1,5),nn(1,7),nel,iddl,iadk,jdik,
139 . k_diag,k_lt,k57,3,off)
140 CALL assem_kij(nn(1,5),nn(1,8),nel,iddl,iadk,jdik,
141 . k_diag,k_lt,k58,3,off)
142 CALL assem_kij(nn(1,6),nn(1,7),nel,iddl,iadk
143 . k_diag,k_lt,k67,3,off)
144 CALL assem_kij(nn(1,6),nn(1,8),nel,iddl,iadk,jdik,
145 . k_diag,k_lt,k68,3,off)
146 CALL assem_kij(nn(1,7),nn(1,8),nel,iddl,iadk,jdik,
147 . k_diag, k_lt,k78,3,off)
148
149
150 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)