assem_q4.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

assem_q4()

subroutine assem_q4	(	integer, dimension(nixq,*)	ixq,
		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_q4.F.

      use element_mod , only : nixq
C-----------------------------------------------
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXQ(NIXQ,*), NEL, IDDL(*), NDOF(*), IADK(*), JDIK(*)
      my_real
     .   k_diag(*), k_lt(*),
     .   k11(2,2,*), k12(2,2,*), k13(2,2,*), k14(2,2,*), k22(2,2,*), 
     .   k23(2,2,*), k24(2,2,*), k33(2,2,*), k34(2,2,*), k44(2,2,*),
     .   off(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, EP, NN(MVSIZ,4)
      my_real
     .   ktp(3,3,mvsiz)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      DO i=1,4
        DO ep=1,nel
          nn(ep,i) = ixq(i+1,ep)
        ENDDO
      ENDDO
C
      DO ep=1,nel
        DO j=1,3
          ktp(1,j,ep) = zero
        ENDDO
        DO i=2,3
          ktp(i,1,ep) = zero
        ENDDO
      ENDDO
C-----block diagonal kii
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k11(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k22(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k33(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k44(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt,ktp,3,off)
C-----non diag kij
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k12(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k13(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k14(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k23(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k24(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
      DO ep=1,nel
        DO j=1,2
        DO i=1,2
          ktp(i+1,j+1,ep) = k34(i,j,ep)
        ENDDO
        ENDDO
      ENDDO
      CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
     .               k_diag,k_lt,ktp,3,off)
C
      RETURN