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