Esope source files

   1 : C CUNIAC    SOURCE    KICH      09/09/11    21:17:54     6505           
   2 : C UNIACI    SOURCE    AM        96/12/24    21:33:08     2448           
   3 :       SUBROUTINE CUNIAC(wrk52,wrk53,wrk54,NSTRS1,MFR1)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 : C----------------------------------------------------------------------
   7 : C
   8 : C     COMPORTEMENT UNIDIRECTIONNEL ACIER_UNI (2D ou coque)
   9 : C     (D'APRES FIBSTE DU MODELE A FIBRE)
  10 : C
  11 : C----------------------------------------------------------------------
  12 : C     Pierre Pegon (ISPRA) Juillet 1994
  13 : C----------------------------------------------------------------------
  14 : C      XMAT( 1) YOUNG
  15 : C      XMAT( 2) NU
  16 : C      XMAT( 3) V1X
  17 : C      XMAT( 4) V1Y
  18 : C      XMAT( 5) MAX. ELAST. STRESS
  19 : C      XMAT( 6) ULTIMATE STRESS
  20 : C      XMAT( 7) HARDENING STRAIN
  21 : C      XMAT( 8) ULTIMATE STRAIN
  22 : C      XMAT( 9) Ro FACTOR
  23 : C      XMAT(10) b  FACTOR
  24 : C      XMAT(11) a1 FACTOR
  25 : C      XMAT(12) a2 FACTOR
  26 : C      XMAT(13) FACTOR L/D
  27 : C      XMAT(14) a6 FACTOR
  28 : C      XMAT(15) c FACTOR
  29 : C      XMAT(16) a FACTOR
  30 : C----------------------------------------------------------------------
  31 : C      VAR0( 1) VARF( 1) FLAG
  32 : C      VAR0( 2) VARF( 2) EPSO
  33 : C      VAR0( 3) VARF( 3) R FACTOR
  34 : C      VAR0( 4) VARF( 4) TANGENT MODULUS
  35 : C      VAR0( 5) VARF( 5) LAST REV. STRAIN
  36 : C      VAR0( 6) VARF( 6) INTERSEPTION POINT
  37 : C      VAR0( 7) VARF( 7) MAX. REV. POST. STRAIN
  38 : C      VAR0( 8) VARF( 8) MIN. REV. POST. STRAIN
  39 : C      VAR0( 9) VARF( 9) LAST REV. STRESS
  40 : C      VAR0(10) VARF(10) INTERSEPTION POINT
  41 : C      VAR0(11) VARF(11) MAX. REV. POST. STRESS
  42 : C      VAR0(12) VARF(12) MIN. REV. POST. STRESS
  43 : C----------------------------------------------------------------------
  44 : C      xcarb( 1) EPAIS (if any)
  45 : C----------------------------------------------------------------------
  46 : -INC CCOPTIO
  47 : -INC DECHE
  48 : 
  49 :       DIMENSION EPSON(5),SIGMA(5)
  50 :       DIMENSION COBMA(10)
  51 : *
  52 :       KERRE=0
  53 : *
  54 : C
  55 : C Calcul des contraintes unidirectionnelles effectives
  56 : C
  57 :       YOUN=XMAT(1)
  58 :       V1X=XMAT(3)
  59 :       V1Y=XMAT(4)
  60 :       IF (IFOUR.EQ.2.AND.(MFR1.EQ.3.OR.MFR1.EQ.9))THEN
  61 :         EPAI=xcarb(1)
  62 :         S0IUNI=(SIG0(1)*V1X**2+SIG0(2)*V1Y**2+V1X*V1Y*SIG0(3))/EPAI
  63 :         DSIUNI=(DEPST(1)*V1X**2+DEPST(2)*V1Y**2+V1X*V1Y*DEPST(3))
  64 :         DSIUNI=YOUN*DSIUNI 
  65 :       ELSEIF (IFOUR.LE.0.AND.MFR1.EQ.1)THEN
  66 :         S0IUNI=SIG0(1)*V1X**2+SIG0(2)*V1Y**2+V1X*V1Y*SIG0(4)
  67 :         DSIUNI=DEPST(1)*V1X**2+DEPST(2)*V1Y**2+V1X*V1Y*DEPST(4)
  68 :         DSIUNI=YOUN*DSIUNI
  69 :       ELSE
  70 :         KERRE=99
  71 :         RETURN
  72 :       ENDIF
  73 : C
  74 :       EYOUN = YOUN 
  75 :       FSY__ = XMAT( 5)
  76 :       FSU__ = XMAT( 6)
  77 :       EPSH_ = XMAT( 7)
  78 :       EPSU_ = XMAT( 8)
  79 :       ROFAC = XMAT( 9)
  80 :       BFACT = XMAT(10)
  81 :       A1FAC = XMAT(11)
  82 :       A2FAC = XMAT(12)
  83 :       FACLD = XMAT(13)
  84 :       A6FAC = XMAT(14)
  85 :       CFACT = XMAT(15)
  86 :       AFACT = XMAT(16)
  87 : C
  88 :       IFLAG = NINT(VAR0(1))
  89 :       STRAI = VAR0(2) + DSIUNI/EYOUN
  90 :       RFACT = VAR0(3)
  91 :       TANGE = VAR0(4)
  92 : C
  93 :       DO IFACT = 1, 2
  94 :          EPSON(IFACT) = VAR0(4+IFACT)
  95 :          SIGMA(IFACT) = VAR0(8+IFACT)
  96 :       END DO
  97 : C
  98 :       EPSON(3) = VAR0(2)
  99 :       SIGMA(3) = S0IUNI
 100 : C
 101 :       DO IFACT = 4, 5
 102 :          EPSON(IFACT) = VAR0(3+IFACT)
 103 :          SIGMA(IFACT) = VAR0(7+IFACT)
 104 :       END DO
 105 : C
 105 :       CALL STEEL1
 106 :      .       ( EYOUN, FSY__, FSU__, EPSH_, EPSU_, ROFAC, BFACT,
 107 :      .         A1FAC, A2FAC, FACLD, A6FAC, CFACT, AFACT, IFLAG,
 108 :      .         STRAI, RFACT, EPSON, SIGMA, TANGE, STRES )
 109 : C
 110 :       VARF(1) = DBLE(IFLAG)
 111 :       VARF(2) = STRAI
 112 :       VARF(3) = RFACT
 113 :       VARF(4) = TANGE
 114 : C
 115 :       DO IFACT = 1, 2
 116 :          VARF(4+IFACT) = EPSON(IFACT)
 117 :          VARF(8+IFACT) = SIGMA(IFACT)
 118 :       END DO
 119 : C
 120 :       DO IFACT = 4, 5
 121 :          VARF(3+IFACT) = EPSON(IFACT)
 122 :          VARF(7+IFACT) = SIGMA(IFACT)
 123 :       END DO
 124 : C
 125 :       IF (IFOUR.EQ.2.AND.(MFR1.EQ.3.OR.MFR1.EQ.9))THEN
 126 :         EPAI=xcarb(1)
 127 :         SIGF(1)=EPAI*STRES*V1X**2
 128 :         SIGF(2)=EPAI*STRES*V1Y**2
 129 :         SIGF(3)=EPAI*STRES*V1X*V1Y
 130 :         DO ISTRS=4,NSTRS1
 131 : C PPDC
 132 : C          SIGF(ISTRS)=SIG0(ISTRS)+DSIGT(ISTRS)
 133 :           SIGF(ISTRS)=SIG0(ISTRS)
 134 : C PPDC
 135 :         ENDDO
 136 :       ELSEIF (IFOUR.LE.0.AND.MFR1.EQ.1)THEN
 137 :         SIGF(1)=STRES*V1X**2
 138 :         SIGF(2)=STRES*V1Y**2
 139 :         SIGF(3)=0.D0
 140 :         SIGF(4)=STRES*V1X*V1Y
 141 :       ENDIF
 142 : C
 143 :       DO ISTRS=1,NSTRS1
 144 :         DEFP(ISTRS)=DEPST(ISTRS)-((SIGF(ISTRS)-SIG0(ISTRS)))/EYOUN
 145 :       END DO
 146 : C
 147 :  1000 RETURN
 148 :       END
 149 :  
 150 :  
 151 :  
 152 :  
 153 :  
 154 :  
 155 :