Esope source files

   1 : C MATHOK    SOURCE    CHAT      05/01/13    01:34:13     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 : C Matriz de Hook:
   6 : C     n   : dimension de la matriz
   7 : C     nsig: dimension de sigma
   8 : C     ifl : +1 matriz de Hook
   9 : C           -1 inversa
  10 : C Crisfield 1, pp107 (+1) (sigma, tau)=E (epsil, gamma=2 epsil)
  11 : CCC
  12 :       SUBROUTINE MATHOK (E,N,NSIG,IFL)
  13 :       IMPLICIT INTEGER(I-N)
  14 :       integer n,nsig,i,j,ifl
  15 :       real*8 E(n,n),aux0,aux1,aux2,aux3,you,xnu
  16 :       common /elasdata/ you,xnu
  17 :       if (ifl.eq.1) then
  18 :          aux0=you/((1.D0+xnu)*(1.D0-2.D0*xnu))
  19 :          aux1=(1.D0-xnu)*aux0
  20 :          aux2=xnu*aux0
  21 :          aux3=you/(2.D0*(1.D0+xnu))
  22 :       else if (ifl.eq.-1) then
  23 :          aux1=1.D0/you
  24 :          aux2=-xnu/you
  25 :          aux3=2.D0*(1.D0+xnu)/you      
  26 :       else
  27 :          stop ' Error peticion a MatHook'
  28 :       endif
  29 :       do i=1,n
  30 :       do j=1,n
  31 :         E(i,j)=0.D0
  32 :       enddo
  33 :       enddo
  34 :       do i=1,3
  35 :          E(i,i)=aux1
  36 :       enddo
  37 :       E(1,2)=aux2
  38 :       E(1,3)=aux2
  39 :       E(2,1)=aux2
  40 :       E(2,3)=aux2
  41 :       E(3,1)=aux2
  42 :       E(3,2)=aux2
  43 :       if (n.lt.4) return
  44 :       do i=4,nsig
  45 :         E(i,i)=aux3
  46 :       enddo      
  47 :       return
  48 :       end
  49 :  
  50 :