37 use element_mod , only : nixtg
38
39
40
41#include "implicit_f.inc"
42
43
44
45#include "mvsiz_p.inc"
46
47
48
49
50
51
52 INTEGER
53 . IXTG(NIXTG,*), NEL ,IDDL(*) ,NDOF(*) ,
54 . IADK(*) ,JDIK(*)
55
57 . k_diag(*) ,k_lt(*) ,off(*),
58 . kc11(6,6,*),kc12(6,6,*),kc13(6,6,*),
59 . kc22(6,6,*),kc23(6,6,*),kc33(6,6,*)
60
61
62
63 INTEGER ,EP,NN(MVSIZ,3)
64
65 DO k = 1,3
66 DO ep = 1,nel
67 nn(ep,k)=ixtg(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
78 CALL assemc_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
79 . k_diag,k_lt,kc12,6,off,ndof)
80 CALL assemc_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
81 . k_diag,k_lt,kc13,6,off,ndof)
82 CALL assemc_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
83 . k_diag,k_lt,kc23,6,off,ndof)
84
85 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)