1 : C UO2CE SOURCE CHAT 05/01/13 03:58:49 5004 2 : SUBROUTINE UO2CE(MC,MM,SIGMA,DXV1,W,WMAX,SMAX,WRUPT,XLTR,XINVL, 3 : & BTR,NFISSU,NVF,FCRIT,VF,YOUN,PRECIZ,JRESU,KERRE) 4 : C---------------------------------------------------------------------- 5 : C MODELE UO2 (OTTOSEN ET GATT_MONERIE) 6 : C---------------------------------------------------------------------- 7 : C 8 : C SIGNIFICATION DU CONTENU DU TABLEAU MM 9 : C -------------------------------------- 10 : C 1 NOUVELLE FISSURATION SELON LA DIRECTION 1 11 : C 2 NOUVELLE FISSURATION SELON LA DIRECTION 2 12 : C 3 NOUVELLE FISSURATION SELON LA DIRECTION 3 13 : C 4 POURSUITE DE LA FISSURATION SELON LA DIRECTION 1 14 : C 5 POURSUITE DE LA FISSURATION SELON LA DIRECTION 2 15 : C 6 POURSUITE DE LA FISSURATION SELON LA DIRECTION 3 16 : C 7 RUPTURE SELON LA DIRECTION 1 17 : C 8 RUPTURE SELON LA DIRECTION 2 18 : C 9 RUPTURE SELON LA DIRECTION 3 19 : C 10 REOUVERTURE SELON LA DIRECTION 1 20 : C 11 REOUVERTURE SELON LA DIRECTION 2 21 : C 12 REOUVERTURE SELON LA DIRECTION 3 22 : C 13 REFERMETURE SELON LA DIRECTION 1 23 : C 14 REFERMETURE SELON LA DIRECTION 2 24 : C 15 REFERMETURE SELON LA DIRECTION 3 25 : C 26 : C ENTREES 27 : C ------- 28 : C MM(20) = TYPES DES CRIT. DE FISS. SUSCEPTIBLES D ETRE ATTEINTS 29 : C MC = NBR. DE CRITERES DE FISS. SUSCEPTIBLES D ETRE ATTEINTS 30 : C SIGMA(6) = CONTRAINTES 31 : C DXV1(3) = INCREMENT DES DEFORMATIONS DE FISSURATION (OUV.) 32 : C W(3) = OUVERTURES DES FISSURES 33 : C WMAX(3) = OUVERTURES MAXIMALES DES FISSURES 34 : C SMAX(3) = CONTR. CORRESPONDANT A WMAX 35 : C WRUPT(3) = OUVERTURES CONDITIONNANT LA RUPTURE 36 : C XLTR(3) = LIMITES EN TRACTION POUR LA FISSURATION 37 : C XINVL(3) = PARAMETRES DE TAILLE 38 : C BTR = PARAMETRE DE FERMETURE 39 : C NFISSU = NOMBRE DE FISSURES 40 : C NVF = NBR DE DIRECTIONS IMPOSEES POUVANT DEVENIR 41 : C DES DIRECTIONS DE FISSURATION 42 : C YOUN = MODULE D YOUNG 43 : C PRECIZ = PRECISION POUR TESTS SUR CONTRAINTES 44 : C 45 : C SORTIES 46 : C ------- 47 : C FCRIT(20) = CRIT. DE FISS. SUSCEPTIBLES D ETRE ATTEINTS 48 : C VF(3,3) = VECTEURS DES DIRECTIONS DE FISSURATION 49 : C NON NULS SEULEMENT SI NOUVELLES FISSURATIONS POSSIBLES 50 : C ET CRITERES CORRESPONDANTS ATTEINTS 51 : C JRESU(3) = 1 2 OU 3 SELON QUE LA DIRECTION POTENTIELLE DE FISS. 52 : C EST CONNUE, INCONNUE EN PLAN OU INCONNUE EN 3D 53 : C KERRE = GESTION DES ERREURS 54 : C---------------------------------------------------------------------- 55 : C 56 : IMPLICIT INTEGER(I-N) 57 : IMPLICIT REAL*8(A-H,O-Z) 58 : -INC CCOPTIO 59 : C 60 : DIMENSION MM(*),SIGMA(*),DXV1(*) 61 : DIMENSION W(*),WMAX(*),SMAX(*),WRUPT(*),XLTR(*),XINVL(*) 62 : DIMENSION FCRIT(*),VF(3,*),JRESU(*) 63 : DIMENSION FCF(3) 64 : C 65 : C ----------------------- 66 : C INITIALISATIONS 67 : C ----------------------- 68 : KERRE=0 69 : DO IC=1,15 70 : FCRIT(IC)=-1.D4*PRECIZ 71 : ENDDO 71 : CALL IANUL(JRESU,3) 72 : C 73 : C 74 : LAPPEL=0 75 : DO IC=1,MC 76 : JC=MM(IC) 77 : IF(JC.GE.1.AND.JC.LE.3) THEN 78 : LAPPEL=1 79 : GO TO 200 80 : ENDIF 81 : ENDDO 82 : C 83 : 200 CONTINUE 84 : C 85 : C --- pas d appel a OTTOFI en cas de pre-fissuration 86 : C 87 : IF(XINVL(1)*XINVL(2)*XINVL(3).EQ.0.D0.AND.LAPPEL.EQ.1) THEN 87 : CALL OTTOFI(SIGMA,NFISSU,NVF,XLTR,XINVL, 88 : & VF,JRESU,FCF,PRECIZ,KERRE) 89 : IF(KERRE.NE.0) RETURN 90 : IF(IIMPI.EQ.42) THEN 91 : WRITE(IOIMP,76600) 92 : 76600 FORMAT(// 2X, ' UO2CE - APRES OTTOFI ' /) 93 : WRITE(IOIMP,76601) (JRESU(IC),IC=1,3) 94 : 76601 FORMAT(2X, ' JRESU ',2X,3I3/) 95 : ENDIF 96 : ENDIF 97 : C 98 : C 99 : DO IC=1,MC 100 : JC=MM(IC) 101 : C 102 : GO TO (1,1,1,4,4,4,7,7,7,10,10,10,13,13,13),JC 103 : C 104 : KERRE=99 105 : RETURN 106 : C 107 : 1 CONTINUE 108 : C 109 : C nouvelle fissure selon une direction 110 : C ------------------------------------ 111 : IF(JRESU(JC).NE.0) THEN 112 : FCRIT(JC)=FCF(JC) 113 : ENDIF 114 : GO TO 100 115 : C 116 : 4 CONTINUE 117 : KC=JC-3 118 : C 119 : C poursuite de la fissuration selon la direction KC ( 1, 2 OU 3 ) 120 : C --------------------------------------------------------------- 121 : FCRIT(JC)= SIGMA(KC) - SMAX(KC) 122 : GO TO 100 123 : C 124 : 7 CONTINUE 125 : KC=JC-6 126 : C 127 : C rupture selon la direction KC ( 1, 2 OU 3 ) 128 : C ------------------------------------------- 129 : FCRIT(JC)= -SIGMA(KC) 130 : GO TO 100 131 : C 132 : 10 CONTINUE 133 : KC=JC-9 134 : C 135 : C reouverture selon la direction KC ( 1, 2 OU 3 ) 136 : C ----------------------------------------------- 137 : FCRIT(JC)= SIGMA(KC) 138 : GO TO 100 139 : C 140 : 13 CONTINUE 141 : KC=JC-12 142 : C 143 : C refermeture selon la direction KC ( 1, 2 OU 3 ) 144 : C ----------------------------------------------- 145 : WREOUV=BTR*MIN(WMAX(KC),WRUPT(KC)) 146 : WAUX=W(KC)+DXV1(KC)/XINVL(KC) 147 : FCRIT(JC)=YOUN*(WREOUV-WAUX)*XINVL(KC) 148 : GO TO 100 149 : C 150 : C 151 : 100 CONTINUE 152 : C 153 : ENDDO 154 : C 155 : C IF(IIMPI.EQ.42) THEN 156 : C WRITE(IOIMP,77000) 157 : C77000 FORMAT(// 2X, ' EN SORTIE DE UO2CE ' /) 158 : C DO IC=1,MC 159 : C JC=MM(IC) 160 : C WRITE(IOIMP,77001) IC,JC,FCRIT(JC) 161 : C77001 FORMAT( 2X, ' UO2CE - IC=',I3,2X,'JC=',I3,2X, 162 : C & 'FCRIT=',1PE12.5/) 163 : C ENDDO 164 : C ENDIF 165 : C 166 : RETURN 167 : END 168 : 169 :
© Cast3M 2003 - All rights reserved.
Disclaimer