Esope source files

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