r12sumg3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	r12sumg3 (jft, jlt, kx, dkx, r11, r21, r31, r12, r22, r32, ke11, ke12, ke13, ke22, ke23, ke33)

Function/Subroutine Documentation

r12sumg3()

subroutine r12sumg3	(	integer	jft,
		integer	jlt,
			kx,
			dkx,
			r11,
			r21,
			r31,
			r12,
			r22,
			r32,
			ke11,
			ke12,
			ke13,
			ke22,
			ke23,
			ke33 )

Definition at line 28 of file r12sumg3.F.

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
#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 
     .    kx(*),dkx(*),r11(*),r21(*),r31(*),r12(*),r22(*),r32(*)
       my_real
     .    ke11(3,3,*),ke22(3,3,*),ke33(3,3,*),ke12(3,3,*),
     .    ke13(3,3,*),ke23(3,3,*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I, J, EP
       my_real
     .    k1(3,3,mvsiz),k2(3,3,mvsiz),k3(3,3,mvsiz)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
C---------------------------------------
C   TRANS LOCAL-->GLOBAL 
C---------------------------------------
       DO ep=jft,jlt
        k1(1,1,ep)=r11(ep)*r11(ep)
        k1(1,2,ep)=r11(ep)*r21(ep)
        k1(1,3,ep)=r11(ep)*r31(ep)
        k1(2,2,ep)=r21(ep)*r21(ep)
        k1(2,3,ep)=r21(ep)*r31(ep)
        k1(3,3,ep)=r31(ep)*r31(ep)
        k1(2,1,ep)=k1(1,2,ep)
        k1(3,2,ep)=k1(2,3,ep)
        k1(3,1,ep)=k1(1,3,ep)
       ENDDO 
       DO ep=jft,jlt
        k2(1,1,ep)=r12(ep)*r12(ep)
        k2(1,2,ep)=r12(ep)*r22(ep)
        k2(1,3,ep)=r12(ep)*r32(ep)
        k2(2,2,ep)=r22(ep)*r22(ep)
        k2(2,3,ep)=r22(ep)*r32(ep)
        k2(3,3,ep)=r32(ep)*r32(ep)
        k2(2,1,ep)=k2(1,2,ep)
        k2(3,2,ep)=k2(2,3,ep)
        k2(3,1,ep)=k2(1,3,ep)
       ENDDO 
       DO ep=jft,jlt
        k3(1,1,ep)=r11(ep)*r12(ep)
        k3(1,2,ep)=r11(ep)*r22(ep)
        k3(1,3,ep)=r11(ep)*r32(ep)
        k3(2,1,ep)=r21(ep)*r12(ep)
        k3(2,2,ep)=r21(ep)*r22(ep)
        k3(2,3,ep)=r21(ep)*r32(ep)
        k3(3,1,ep)=r31(ep)*r12(ep)
        k3(3,2,ep)=r31(ep)*r22(ep)
        k3(3,3,ep)=r31(ep)*r32(ep)
       ENDDO 
C
       DO i=1,3 
        DO j=i,3 
         DO ep=jft,jlt
          ke11(i,j,ep)=kx(ep)*k1(i,j,ep)
          ke22(i,j,ep)=kx(ep)*(k1(i,j,ep)+k2(i,j,ep)+
     .                         k3(i,j,ep)+k3(j,i,ep))
          ke33(i,j,ep)=kx(ep)*k2(i,j,ep)
         ENDDO
        ENDDO
       ENDDO
       DO i=1,3 
        DO j=1,3 
         DO ep=jft,jlt
          ke12(i,j,ep)=-kx(ep)*(k1(i,j,ep)+k3(i,j,ep))
          ke13(i,j,ep)=kx(ep)*k3(i,j,ep)
          ke23(i,j,ep)=-kx(ep)*(k2(i,j,ep)+k3(i,j,ep))
         ENDDO
        ENDDO
       ENDDO
C-----frotement---
       DO ep=jft,jlt
        IF (dkx(ep)>zero) THEN
         DO i=1,3 
          DO j=i,3 
           ke11(i,j,ep)=ke11(i,j,ep)+dkx(ep)*k1(i,j,ep)
           ke22(i,j,ep)=ke22(i,j,ep)+dkx(ep)*(k1(i,j,ep)-k3(i,j,ep))
           ke33(i,j,ep)=ke33(i,j,ep)+dkx(ep)*k2(i,j,ep)
           ke22(i,j,ep)=ke22(i,j,ep)+dkx(ep)*(k2(i,j,ep)-k3(j,i,ep))
          ENDDO
         ENDDO
         DO i=1,3 
          DO j=1,3 
           ke12(i,j,ep)=ke12(i,j,ep)-dkx(ep)*(k1(i,j,ep)-k3(i,j,ep))
           ke13(i,j,ep)=ke13(i,j,ep)-dkx(ep)*k3(i,j,ep)
           ke23(i,j,ep)=ke23(i,j,ep)-dkx(ep)*(k2(i,j,ep)-k3(i,j,ep))
          ENDDO
         ENDDO
        ENDIF 
       ENDDO
C
       DO i=1,3 
        DO j=i+1,3 
         DO ep=jft,jlt
          ke11(j,i,ep)=ke11(i,j,ep)
          ke22(j,i,ep)=ke22(i,j,ep)
          ke33(j,i,ep)=ke33(i,j,ep)
         ENDDO
        ENDDO
       ENDDO
C
      RETURN