pke3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	pke3 (jft, jlt, nel, mtn, ismstr, pm, ncc, x, elbuf_tab, geo, offset, ikgeo, etag, iddl, ndof, k_diag, k_lt, iadk, jdik)
subroutine	pkeg3 (jft, jlt, fx, al, ke11, ke12, ke22)

Function/Subroutine Documentation

pke3()

subroutine pke3	(	integer	jft,
		integer	jlt,
		integer	nel,
		integer	mtn,
		integer	ismstr,
			pm,
		integer, dimension(nixp,*)	ncc,
			x,
		type (elbuf_struct_), target	elbuf_tab,
			geo,
		integer	offset,
		integer	ikgeo,
		integer, dimension(*)	etag,
		integer, dimension(*)	iddl,
		integer, dimension(*)	ndof,
			k_diag,
			k_lt,
		integer, dimension(*)	iadk,
		integer, dimension(*)	jdik )

Definition at line 39 of file pke3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod
      use element_mod , only : nixp
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"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NCC(NIXP,*),JFT,JLT,OFFSET,NEL ,MTN,ISMSTR,IKGEO,
     .        ETAG(*),IDDL(*)  ,NDOF(*)  ,IADK(*) ,JDIK(*)     
      my_real 
     .   pm(*), x(*),geo(*),
     .   ke11(36,mvsiz),ke12(36,mvsiz),ke22(36,mvsiz),off(mvsiz),
     .   k_diag(*) ,k_lt(*)
C
      TYPE (ELBUF_STRUCT_), TARGET:: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER         
     .           I 
      INTEGER  MAT(MVSIZ),PID(MVSIZ),NGL(MVSIZ)
      my_real
     .    ym(mvsiz),g(mvsiz),al(mvsiz),ali(mvsiz),
     .    sh1(mvsiz),sh2(mvsiz),
     .    area(mvsiz),b1(mvsiz),b2(mvsiz),b3(mvsiz),
     .    r11(mvsiz),r12(mvsiz),r13(mvsiz),
     .    r21(mvsiz),r22(mvsiz),r23(mvsiz),
     .    r31(mvsiz),r32(mvsiz),r33(mvsiz),fx(mvsiz)
      TYPE(G_BUFEL_)  ,POINTER :: GBUF
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      gbuf => elbuf_tab%GBUF
 
      CALL pcoork3(jft         ,jlt      ,x      ,ncc   ,gbuf%OFF,
     2             pid         ,mat      ,ngl    ,al    ,ali     , 
     3             gbuf%LENGTH,ismstr    ,r11    ,r12   ,r13     ,
     4             r21        ,r22       ,r23    ,r31   ,r32     ,
     5             r33        ,gbuf%SKEW ,off    )
      CALL pmat3(jft    ,jlt    ,pm     ,geo     ,ym      ,
     2           g      ,area   ,b1     ,b2     ,b3       ,
     3           al     ,ali    ,sh1    ,sh2    ,mtn      ,
     4           mat    ,pid    )
C-------------------------
C     STIFFNESS MATRIX
C-------------------------
      CALL psumg3(jft    ,jlt    ,ym      ,g      ,area   ,
     2           b1     ,b2     ,b3      ,al     ,ali    ,
     3           sh1    ,sh2    ,r11     ,r12    ,r13    ,
     4           r21    ,r22    ,r23     ,r31    ,r32    ,
     5           r33    ,ke11   ,ke12    ,ke22   )
      IF (ikgeo>0) THEN
       DO i =jft ,jlt
         fx(i)=gbuf%FOR(2*nel+i)
       ENDDO 
       CALL pkeg3(jft    ,jlt    ,
     1            fx     ,al     ,ke11    ,ke12    ,ke22   )
      ENDIF 
      
             IF (neig>0) CALL peoff(
     1                   nixp, jft, jlt, ncc, etag,
     2                   off )
             CALL assem_p(nixp,6 ,
     1         ncc       ,nel       ,iddl      ,ndof      ,k_diag    ,
     2         k_lt      ,iadk      ,jdik      ,ke11      ,ke12      ,
     3         ke22      ,off       )
C
      RETURN

pkeg3()

subroutine pkeg3	(	integer	jft,
		integer	jlt,
			fx,
			al,
			ke11,
			ke12,
			ke22 )

Definition at line 129 of file pke3.F.

C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT
      my_real 
     .    fx(*),al(*),ke11(6,6,*),ke12(6,6,*),ke22(6,6,*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER EP,I
      my_real 
     .    h(mvsiz)
C-----------
       DO ep=jft,jlt 
        h(ep)=fx(ep)/al(ep)
       ENDDO
       DO i=1,3
        DO ep=jft,jlt
         ke11(i,i,ep) = ke11(i,i,ep)+h(ep)
         ke12(i,i,ep) = ke12(i,i,ep)-h(ep)
         ke22(i,i,ep) = ke22(i,i,ep)+h(ep)
        ENDDO
       ENDDO
C
      RETURN