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