Esope source files

   1 : C TUFINC    SOURCE    PV        11/03/07    21:18:39     6885           
   2 : C TUFINC    SOURCE    BROC        94/05/09
   3 : C______________________________________________________________________
   4 :       SUBROUTINE TUFINC(DT,SIG,EPSV,VAR,XMAT,XCAR,EPSVPT,VARPT,
   5 :      &     NSTRS,NVARI,INPLAS,NCOMAT,KERR1)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 : C    Include contenant quelques constantes dont XPI :
   9 : -INC CCREEL
  10 :       DIMENSION SIG(*),EPSV(*),VAR(*),XMAT(*),XCAR(*),
  11 :      &     EPSVPT(*),VARPT(*)
  12 : C______________________________________________________________________
  13 : C CETTE ROUTINE CALCULE A PARTIR DES VALEURS INITIALES SIG VAR & EPSV
  14 : C                ET DE LA VARIATION DU TEMPS DT
  15 : C LES INCREMENTS DE DEFORMATIONS VISQUEUSES ET DE VARIABLES INTERNES
  16 : C____________________________________________________________________
  17 : 
  18 :       THETA = XCAR(9)/2
  19 :       EPAIS = XCAR(2)
  20 :       RAYOM = XCAR(1)-(EPAIS/2)
  21 :       AF1 = XMAT(6)
  22 :       AF2 = XMAT(7)
  23 :       ANG=XPI*THETA/180.0D0
  24 : 
  25 :       EPLU = VAR (1)
  26 :       EMIN = VAR (2)
  27 :       EEFF    = VAR (3)
  28 :       EPSH    = VAR (4)
  29 :       LL      = nint(VAR (5))
  30 :       DO 30 I=1,6
  31 :          EPSVPT(I)=0.0D0
  32 :  30   CONTINUE
  33 :       EPSCRP= EPSV(6)
  34 :       DO 50 I=1,NVARI
  35 :          VARPT(I)=0.0D0
  36 :  50   CONTINUE
  37 :       KERR1=0
  38 : 
  39 : C     
  40 : C  CALCUL DES INCREMENTS DE VARIABLES INTERNES
  41 : C     
  41 :       CALL TUFVAR(SIG,EPSCRP,EPLU,EMIN,EEFF,EPSH,LL,NSTRS)
  42 :       VARPT (1)   =  EPLU
  43 :       VARPT (2)   =  EMIN
  44 :       VARPT (3) =  EEFF
  45 :       VARPT (4) =  EPSH
  46 :       VARPT (5) =  LL
  47 : C     
  48 : C CALCUL DES INCREMENTS DE EPSILON
  49 : C     
  50 : 
  51 : C EVALUATION DE C* MOYENNANT ZAHOOR
  52 : C---- C* CALCULE QUE POUR DES CHARGEMENTS EN OUVERTURE DE FISSURE
  53 : C----ATTENTION SIG(6)<0 SI MOMENT EN OUVERTURE
  54 : 
  55 :       IF (SIG(6).GT.(1.D-6).OR.SIG(1).LT.(-1.D-6)) THEN
  56 :          KERR1=729
  57 :          RETURN
  58 :       ENDIF
  59 :       XM0= 4*RAYOM*RAYOM*EPAIS*(COS(ANG/2)-0.5D0*SIN(ANG))
  60 :       IF (ABS(SIG(6)).GT.(0.99D0*SIG(1)*RAYOM)) THEN
  61 :          XLSLP1 = 1.D0
  62 :       ELSE
  63 :          XLAMDA = ABS(SIG(6))/(SIG(1)*RAYOM)
  64 :          XP0 = 2*RAYOM*EPAIS*(XPI-ANG - 2*ASIN(0.5D0*SIN(ANG)))
  65 :          XPP0=0.5D0*(-1*XLAMDA*RAYOM*XP0*XP0/XM0)
  66 :          XX=0.5D0*(((((XLAMDA*RAYOM*XP0*XP0/XM0)**2) + 4*XP0*XP0)**0.5D0
  67 :      $        ))
  68 :          XPP0=XPP0+XX
  69 :          XLSLP1=XLAMDA/(XLAMDA+1)
  70 :       ENDIF
  71 :       RAP=RAYOM/EPAIS
  72 :       AN=ANG/XPI
  72 :       CALL TUFH1(RAP,AF2,AN,XLSLP1,H1,KERR1)
  73 :       IF (KERR1.NE.0) THEN
  74 :          RETURN
  75 :       ENDIF
  76 :       IF (XLSLP1.GT.(0.99D0)) GOTO 55
  77 :       CETOI=AF1*(XPI-ANG)*RAYOM*(ANG/XPI)
  78 :       CETOI=CETOI*H1*((SIG(1)/XPP0)**(AF2+1))
  79 :       GO TO 60
  80 : 
  81 :  55   CONTINUE
  82 :       CETOI=AF1*XPI*((1-ANG/XPI)**2)*RAYOM
  83 :       CETOI=CETOI*H1*(ABS(SIG(6)/XM0)**(AF2+1))
  84 :  60   CONTINUE
  85 :       VARPT (7) =CETOI/DT
  86 : 
  87 : C RECHERCHE DE LA POSITION ALPHA DE LA FIBRE NEUTRE
  88 : 
  89 :       ALPH1=ANG/2
  90 :       ALPH2=ACOS(0.5D0*SIN(ANG))
  91 :       IF (ABS(SIG(6)).GT.(0.99D0*SIG(1)*RAYOM)) THEN
  92 :          ALPH=ALPH1
  93 :          GOTO 120
  94 :       ELSE
  95 :          X=XPI*ABS(SIG(6))/(SIG(1)*2*RAYOM)
  96 :          IF (X.LT.1.D-6) THEN
  97 :             ALPH=ALPH2
  98 :             GOTO 120
  99 :          ENDIF
 100 :       ENDIF
 101 : 
 102 :       F1=COS(ALPH1)-0.5D0*SIN(ANG)-X*(ALPH1/XPI-ANG/(2*XPI))
 103 :       F2=COS(ALPH2)-0.5D0*SIN(ANG)-X*(ALPH2/XPI-ANG/(2*XPI))
 104 :       ALPH3=ALPH1-((ALPH2-ALPH1)/(F2-F1))*F1
 105 :       F3=COS(ALPH3)-0.5D0*SIN(ANG)-X*(ALPH3/XPI-ANG/(2*XPI))
 106 : 
 107 :       ERR=F3
 108 : 
 109 :  70   CONTINUE
 110 :       IF ((F1*F3).LT.(0.0D0)) THEN
 111 :          ALPH2=ALPH3
 112 :          F2=F3
 113 :       ENDIF
 114 :       IF ((F2*F3).LT.(0.0D0)) THEN
 115 :          ALPH1=ALPH3
 116 :          F1=F3
 117 :       ENDIF
 118 :       ALPH3=ALPH1-((ALPH2-ALPH1)/(F2-F1))*F1
 119 :       F3=COS(ALPH3)-0.5D0*SIN(ANG)-X*(ALPH3/XPI-ANG/(2*XPI))
 120 :       ERR=COS(ALPH3)-0.5D0*SIN(ANG)-X*(ALPH3/XPI-ANG/(2*XPI))
 121 :       IF (ERR.GT.(1.00D-3)) GOTO 70
 122 : 
 123 :       ALPH=ALPH3
 124 : C EVALUATION DE L INCREMENT DE ROTATION DE FLUAGE
 125 : 
 126 :  120  CONTINUE
 127 :       YM0= 4*RAYOM*RAYOM*EPAIS*( COS(ALPH)-0.5D0*SIN(ANG))
 128 :       YP0= 2*RAYOM*EPAIS*(2*ALPH-ANG)
 129 :       IF (YM0.GE.1.0D-1) THEN
 130 :          EPSVPT(6)=CETOI/((SIG(6)/YM0)*RAYOM*(COS(ANG)+SIN(ALPH)))
 131 :       ELSE
 132 :          EPSVPT(6)=CETOI/((SIG(1)/YP0)*RAYOM*(COS(ANG)+SIN(ALPH)))
 133 :       ENDIF
 134 :       VARPT (6) = ABS(EPSVPT(6))
 135 : 
 136 :       RETURN
 137 :       END
 138 : 
 139 : 
 140 : 
 141 : 
 142 :  
 143 :  
 144 :  
 145 :