Télécharger la plaquette Cast3M


fr en

Télécharger uniakt.eso

Retour à la liste

Numérotation des lignes : 

uniakt

  1. C UNIAKT    SOURCE    CHAT      05/01/13    03:58:26     5004
  2.       SUBROUTINE UNIAKT(XMAT,XVAR,INDIC,XKTAT,KERRE)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. C----------------------------------------------------------------------
  6. C
  7. C     COMPORTEMENT UNIDIRECTIONNEL ACIER_UNI (2D ou coque)
  8. C     (D'APRES FIBSTE DU MODELE A FIBRE)
  9. C     Matrice tangente
  10. C
  11. C----------------------------------------------------------------------
  12. C     Pierre Pegon (ISPRA) Juillet 1994
  13. C     Didier Combescure (ISPRA) Septembre 1995
  14. C----------------------------------------------------------------------
  15. C      XMAT(1) YOUN
  16. C      XMAT(2) V1X
  17. C      XMAT(3) V1Y
  18. C  EPAI (if any) en deformations planes
  19. C----------------------------------------------------------------------
  20. C      XVAR( 1) VARF( 1) TANGENT MODULUS
  21. C----------------------------------------------------------------------
  22. C      XKTAT(4,4)   TANGENT STIFFNESS MATRIX
  23. C----------------------------------------------------------------------
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. *
  28.       DIMENSION XMAT(*),XVAR(*)
  29.       DIMENSION XKTAT(4,4)
  30. *
  31.       KERRE=0
  32. *
  33.       YOUN=XMAT(1)
  34.       V1X=XMAT(2)
  35.       V1Y=XMAT(3)
  36.       TANGE = XVAR(1)
  37.       EPAI = 1.D0
  38.       CALL ZERO(XKTAT,4,4)
  39. C
  40.       IF (INDIC.EQ.2)THEN
  41. C
  42. C  Cas deformations planes
  43. C
  44.         EPAI=XMAT(17)
  45.         XKTAT(1,1)=EPAI*TANGE*V1X**4
  46.         XKTAT(2,2)=EPAI*TANGE*V1Y**4
  47.         XKTAT(2,1)=EPAI*TANGE*(V1X*V1Y)**2
  48.         XKTAT(1,2)=EPAI*TANGE*(V1X*V1Y)**2
  49.         XKTAT(1,3)=EPAI*TANGE*(V1X**3)*V1Y
  50.         XKTAT(2,3)=EPAI*TANGE*V1X*(V1Y**3)
  51.         XKTAT(3,1)=EPAI*TANGE*(V1X**3)*V1Y
  52.         XKTAT(3,2)=EPAI*TANGE*V1X*(V1Y**3)
  53.         XKTAT(3,3)=EPAI*TANGE*(V1X*V1Y)**2
  54.       ELSEIF (INDIC.LE.0)THEN
  55. C
  56. C  Cas contraintes planes
  57. C
  58.        XKTAT(1,1)=TANGE*V1X**4
  59.        XKTAT(2,2)=TANGE*V1Y**4
  60.        XKTAT(2,1)=TANGE*(V1X*V1Y)**2
  61.        XKTAT(1,2)=TANGE*(V1X*V1Y)**2
  62.        XKTAT(1,4)=TANGE*(V1X**3)*V1Y
  63.        XKTAT(2,4)=TANGE*V1X*(V1Y**3)
  64.        XKTAT(4,1)=TANGE*(V1X**3)*V1Y
  65.        XKTAT(4,2)=TANGE*V1X*(V1Y**3)
  66.        XKTAT(4,4)=TANGE*(V1X*V1Y)**2
  67. C       XKTAT(1,1)=TANGE
  68. C       XKTAT(3,3)=TANGE
  69.       ENDIF
  70. C
  71.       RETURN
  72.       END
  73.  
  74.  
  75.  

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