assem_s4.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	assem_s4 (ixs, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, off)

Function/Subroutine Documentation

assem_s4()

subroutine assem_s4	(	integer, dimension(nixs,*)	ixs,
		integer	nel,
		integer, dimension(*)	iddl,
		integer, dimension(*)	ndof,
			k_diag,
			k_lt,
		integer, dimension(*)	iadk,
		integer, dimension(*)	jdik,
			k11,
			k12,
			k13,
			k14,
			k22,
			k23,
			k24,
			k33,
			k34,
			k44,
			off )

Definition at line 33 of file assem_s4.F.

      use element_mod , only : nixs
C----6---------------------------------------------------------------7---------8
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXS(NIXS,*),NEL   ,IDDL(*)  ,NDOF(*)  ,
     .          IADK(*)  ,JDIK(*)        
C     REAL
      my_real
     .   k_diag(*) ,k_lt(*)   ,
     .   k11(9,*) ,k12(9,*)  ,k13(9,*)  ,k14(9,*)  ,k22(9,*)  ,
     .   k23(9,*) ,k24(9,*)  ,k33(9,*)  ,k34(9,*)  ,k44(9,*), 
     .   off(*)  
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER EP,NN(MVSIZ,4)
C----6---------------------------------------------------------------7---------8
        DO ep = 1,nel
         nn(ep,1)=ixs(2,ep)
         nn(ep,2)=ixs(4,ep)
         nn(ep,3)=ixs(7,ep)
         nn(ep,4)=ixs(6,ep)
        ENDDO 
C-------block diagonal k11,k22...k44-----
       CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off) 
       CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off) 
       CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off) 
       CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off) 
C-------non diag kij-----
       CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k12,3,off) 
       CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k13,3,off) 
       CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k14,3,off) 
       CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k23,3,off) 
       CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k24,3,off) 
       CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
     .                 k_diag,k_lt,k34,3,off) 
C
C----6---------------------------------------------------------------7---------8
      RETURN