czfintce.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	czfintce (jft, jlt, thk, c1, a_i, x13, x24, y13, y24, z1, mx23, mx13, mx34, my13, my23, my34, vstre, mstre, vf, vm, nel)
subroutine	czfintcrz (jft, jlt, thk, vol, area, x13, x24, y13, y24, z1, mx23, mx13, mx34, my13, my23, my34, vstre, vsrz, vf, vmz, bm0rz, krz, vrlz, dt1c, eint, off, nel)

Function/Subroutine Documentation

czfintce()

subroutine czfintce	(	integer	jft,
		integer	jlt,
			thk,
			c1,
			a_i,
			x13,
			x24,
			y13,
			y24,
			z1,
			mx23,
			mx13,
			mx34,
			my13,
			my23,
			my34,
			vstre,
			mstre,
			vf,
			vm,
		integer	nel )

Definition at line 29 of file czfintce.F.

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
C        CALCUL LA PARTIE CONSTANTE
C         VF(3,NODE),VM(2,NODE) : FORCES INTERNES GENERALISEES
C        ENTREES  : THK,AREA,X13,X24,Y13,Y24,MX13,MX34,MY13,MY34,Z1,...
C        SORTIES  : VF,VM
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 JFT ,JLT,NEL
        my_real 
     .     thk(*)   ,c1(*)     ,a_i(*)     ,z1(*)     ,
     .     x13(*)   ,x24(*)    ,y13(*)     ,y24(*)    ,
     .     mx13(*)  ,mx23(*)   ,mx34(*)    ,my13(*)   ,my23(*)  ,
     .     my34(*)  ,vstre(nel,5),mstre(nel,3) ,vf(mvsiz,3,4) ,vm(mvsiz,2,4)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
        INTEGER I ,J
        my_real
     .     s1, c2, y13s8, y24s8, s42s, s52s, x13s8, x24s8
C
        DO i=jft,jlt
          x13s8=x13(i)*mstre(i,3)
          x24s8=x24(i)*mstre(i,3)
          y13s8=y13(i)*mstre(i,3)
          y24s8=y24(i)*mstre(i,3)
C
          s1=(my34(i)*mx23(i)-my23(i)*mx34(i))*thk(i)
          s42s=s1*vstre(i,5)
          s52s=s1*vstre(i,4)
C-------------------------------------------------
C        PARTIE ANTISYM STOCKEE DANS NOEUDS 1,2
C        PARTIE SYM NOEUDS STOCKEE DANS NOEUDS 3,4
C-------------------------------------------------
          vf(i,1,1)=thk(i)*( y24(i)*vstre(i,1)-x24(i)*vstre(i,3))
          vf(i,2,1)=thk(i)*(-x24(i)*vstre(i,2)+y24(i)*vstre(i,3))
          vf(i,3,1)=thk(i)*(-x24(i)*vstre(i,4)+y24(i)*vstre(i,5))
C
          vm(i,1,1)=c1(i)*(x24(i)*mstre(i,2)-y24s8)-my13(i)*vf(i,3,1)
          vm(i,2,1)=c1(i)*(y24(i)*mstre(i,1)-x24s8)+mx13(i)*vf(i,3,1)
C
          vm(i,1,3)=-s52s
          vm(i,2,3)= s42s
C
          vf(i,1,2)=thk(i)*(-y13(i)*vstre(i,1)+x13(i)*vstre(i,3))
          vf(i,2,2)=thk(i)*( x13(i)*vstre(i,2)-y13(i)*vstre(i,3))
          vf(i,3,2)=thk(i)*( x13(i)*vstre(i,4)-y13(i)*vstre(i,5))
C
          vm(i,1,2)=c1(i)*(-x13(i)*mstre(i,2)+y13s8)+my13(i)*vf(i,3,2)
          vm(i,2,2)=c1(i)*(-y13(i)*mstre(i,1)+x13s8)-mx13(i)*vf(i,3,2)
C
          vm(i,1,4)=vm(i,1,3)
          vm(i,2,4)=vm(i,2,3)
C
          c2=c1(i)*z1(i)*4.*a_i(i)
C
          vf(i,1,1)=vf(i,1,1)+c2*(x13(i)*mstre(i,1)+y13s8)
          vf(i,2,1)=vf(i,2,1)+c2*(y13(i)*mstre(i,2)+x13s8)
          vf(i,1,2)=vf(i,1,2)-c2*(x24(i)*mstre(i,1)+y24s8)
          vf(i,2,2)=vf(i,2,2)-c2*(y24(i)*mstre(i,2)+x24s8)
       ENDDO
       RETURN

czfintcrz()

subroutine czfintcrz	(	integer	jft,
		integer	jlt,
			thk,
			vol,
			area,
			x13,
			x24,
			y13,
			y24,
			z1,
			mx23,
			mx13,
			mx34,
			my13,
			my23,
			my34,
			vstre,
			vsrz,
			vf,
			vmz,
			bm0rz,
			krz,
			vrlz,
			dt1c,
			eint,
			off,
		integer	nel )

Definition at line 110 of file czfintce.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,NEL
        my_real 
     .     thk(*)   ,vol(*)     ,area(*)     ,z1(*)     ,
     .     x13(*)   ,x24(*)    ,y13(*)     ,y24(*)    ,
     .     mx13(*)  ,mx23(*)   ,mx34(*)    ,my13(*)   ,my23(*) ,
     .     my34(*)  ,vstre(nel,5),vsrz(nel,19) ,vf(mvsiz,3,4) ,vmz(mvsiz,4),
     .     bm0rz(mvsiz,4,4),krz(*),vrlz(*),dt1c,eint(nel,2),off(*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
        INTEGER I ,J
        my_real
     .     s1, c2, a05,erz,off05
C
C      --------const sig_rz added here---------
       DO i=jft,jlt
        erz= vrlz(i)*dt1c
        off05=half*off(i)*vol(i)
        eint(i,1) = eint(i,1)+ vsrz(i,19)*erz*off05
        vsrz(i,19)= vsrz(i,19)+krz(i)*erz*off(i)
        eint(i,1) = eint(i,1)+ vsrz(i,19)*erz*off05
       ENDDO
       DO j=1,4
       DO i=jft,jlt
        vmz(i,j)=(bm0rz(i,1,j)*vstre(i,1)+bm0rz(i,2,j)*vstre(i,2)
     .           +bm0rz(i,3,j)*vstre(i,3))*thk(i)
       ENDDO
       ENDDO
C        
        DO i=jft,jlt
         c2=thk(i)*vsrz(i,19)*half
         DO j=1,4
         vmz(i,j)= vmz(i,j)+bm0rz(i,4,j)*c2
         ENDDO
        ENDDO
C-------------------------------------------------
C        PARTIE ANTISYM STOCKEE DANS NOEUDS 1,2
C        PARTIE SYM NOEUDS STOCKEE DANS NOEUDS 3,4
C-------------------------------------------------
        DO i=jft,jlt
          c2=thk(i)*vsrz(i,19)*half
C
          vf(i,1,1)=vf(i,1,1)+c2*x24(i)
          vf(i,2,1)=vf(i,2,1)+c2*y24(i)
          vf(i,1,2)=vf(i,1,2)-c2*x13(i)
          vf(i,2,2)=vf(i,2,2)-c2*y13(i)
        ENDDO
C
       RETURN