Télécharger mathok.eso

Retour à la liste

Numérotation des lignes :

  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.  
  51.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales