39 use element_mod , only : nixc
40
41
42
43#include "implicit_f.inc"
44
45
46
47#include "mvsiz_p.inc"
48
49
50
51
52
53
54 INTEGER
55 . IXC(NIXC,*), NEL ,IDDL(*) ,
56 . IADK(*) ,JDIK(*) ,NDOF(*)
57
59 . off(*),k_diag(*) ,k_lt(*) ,
60 . kc11(6,6,*),kc12(6,6,*),kc13(6,6,*),
61 . kc22(6,6,*),kc23(6,6,*),kc24(6,6,*),
62 . kc33(6,6,*),kc34(6,6,*),kc44(6,6,*),kc14(6,6,*)
63
64
65
66 INTEGER K,EP,NN(MVSIZ,4)
67
68 DO k = 1,4
69 DO ep = 1,nel
70 nn(ep,k)=ixc(k+1,ep)
71 ENDDO
72 ENDDO
73
74 CALL assemc_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,kc11,6,off,
75 1 ndof)
76 CALL assemc_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,kc22,6,off,
77 1 ndof)
78 CALL assemc_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,kc33,6,off,
79 1 ndof)
80 CALL assemc_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,kc44,6,off,
81 1 ndof)
82
83 CALL assemc_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
84 . k_diag,k_lt,kc12,6,off,ndof)
85 CALL assemc_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
86 . k_diag,k_lt,kc13,6,off,ndof)
87 CALL assemc_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
88 . k_diag,k_lt,kc14,6,off,ndof)
89 CALL assemc_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
90 . k_diag,k_lt,kc23,6,off,ndof)
91 CALL assemc_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
92 . k_diag,k_lt,kc24,6,off,ndof)
93 CALL assemc_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
94 . k_diag,k_lt,kc34,6,off,ndof)
95
96 RETURN
subroutine assemc_kij(ni, nj, nel, iddl, iadk, jdik, k_diag, k_lt, kij, nd, off, ndof)
subroutine assemc_kii(ni, nel, iddl, iadk, k_diag, k_lt, kii, nd, off, ndof)