Esope source files

   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 : -INC CCOPTIO
  25 : *
  26 :       DIMENSION XMAT(*),XVAR(*)
  27 :       DIMENSION XKTAT(4,4) 
  28 : *
  29 :       KERRE=0
  30 : *
  31 :       YOUN=XMAT(1)
  32 :       V1X=XMAT(2)
  33 :       V1Y=XMAT(3)
  34 :       TANGE = XVAR(1)
  35 :       EPAI = 1.D0
  35 :       CALL ZERO(XKTAT,4,4)
  36 : C
  37 :       IF (INDIC.EQ.2)THEN
  38 : C
  39 : C  Cas deformations planes
  40 : C
  41 :         EPAI=XMAT(17)
  42 :         XKTAT(1,1)=EPAI*TANGE*V1X**4
  43 :         XKTAT(2,2)=EPAI*TANGE*V1Y**4
  44 :         XKTAT(2,1)=EPAI*TANGE*(V1X*V1Y)**2
  45 :         XKTAT(1,2)=EPAI*TANGE*(V1X*V1Y)**2
  46 :         XKTAT(1,3)=EPAI*TANGE*(V1X**3)*V1Y
  47 :         XKTAT(2,3)=EPAI*TANGE*V1X*(V1Y**3)
  48 :         XKTAT(3,1)=EPAI*TANGE*(V1X**3)*V1Y
  49 :         XKTAT(3,2)=EPAI*TANGE*V1X*(V1Y**3)
  50 :         XKTAT(3,3)=EPAI*TANGE*(V1X*V1Y)**2
  51 :       ELSEIF (INDIC.LE.0)THEN
  52 : C
  53 : C  Cas contraintes planes
  54 : C
  55 :        XKTAT(1,1)=TANGE*V1X**4
  56 :        XKTAT(2,2)=TANGE*V1Y**4
  57 :        XKTAT(2,1)=TANGE*(V1X*V1Y)**2
  58 :        XKTAT(1,2)=TANGE*(V1X*V1Y)**2
  59 :        XKTAT(1,4)=TANGE*(V1X**3)*V1Y
  60 :        XKTAT(2,4)=TANGE*V1X*(V1Y**3)
  61 :        XKTAT(4,1)=TANGE*(V1X**3)*V1Y
  62 :        XKTAT(4,2)=TANGE*V1X*(V1Y**3)
  63 :        XKTAT(4,4)=TANGE*(V1X*V1Y)**2
  64 : C       XKTAT(1,1)=TANGE
  65 : C       XKTAT(3,3)=TANGE
  66 :       ENDIF
  67 : C
  68 :       RETURN
  69 :       END
  70 :  
  71 :