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