1 : $$$$ @LISPA16 NOTICE MB234859 17/10/02 21:15:01 9577 2 : DATE 17/10/02 3 : 4 : Procedure @LISPA16 Voir aussi : 5 : ------------------ 6 : @LISPA16 TAB1 TAB2 ; 7 : 8 : TAB1.'FAT' .'CC' .'ROO2' .'FLU' .'CT' 9 : .'EPS' .'A316' 10 : .'MFIS' .'MTOT' .'OBJ' .'OBJT' 11 : .'MAT' .'MATT' . 'CAR' .'CART' 12 : .'BLO' 13 : 14 : TAB2.'DP' .'DC' .'PA' .'PB' .'CA' 15 : .'CB' 16 : .'DPTH' .'DCTH' .'TPSM' 17 : 18 : 19 : 20 : 21 : Objet : 22 : ------- 23 : 24 : Cette procedure evalue la propagation en fatigue-fluage le long d'un 25 : front de fissure modelise par des elements LISP, en fonction du niveau 26 : de sollicitation, en appliquant le procedure de l'Annexe A16 du RCC-MR 27 : (Rapport DMT 94/043). 28 : Pour l'acier 316, les caracteristiques du materiau a 525°C sont deja 29 : rentrees dans la procedure. 30 : Cette procedure a ete developpee dans le cadre d'un travail presente 31 : dans le rapport 94/612. 32 : 33 : 34 : 35 : Description des arguments d'entres et de sortie : 36 : ------------------------------------------------- 37 : 38 : 1) Arguments d'entree : 39 : ----------------------- 40 : 41 : TAB1 (type TABLE) : table contenant les donnees du materiau ( a noter 42 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ que les unites sont a respecter). 43 : 44 : 45 : * Pour un materiau autre que A316 46 : 47 : -> donnees en fatigue 48 : 49 : TAB1.FAT.'COEFFICIENT' : coefficient C de la loi de propagation 50 : n 51 : de Paris da/dN = C (DKeff). 52 : 53 : TAB1.FAT.'EXPOSANT' : exposant n de la loi de propagation 54 : n 55 : de Paris da/dN = C (DKeff) . 56 : 57 : C et n sont entres pour un DKeff exprime en MPaVm et da/dN en 58 : m/cycle. 59 : TAB1.CC.'EPS' : abscisses de la courbe de traction 60 : cyclique du materiau. 61 : 62 : TAB1.CC.'SIG' : ordonnees de la courbe de traction 63 : cyclique du materiau; 64 : 65 : TAB1.'ROO2' : limite d'elasticite du materiau a 0,2%. 66 : 67 : 68 : -> donnees en fluage 69 : 70 : TAB1.FLU.'COEFFICIENT' : coefficient A de la loi de propagation en 71 : fluage 72 : *q 73 : da/dt = A C . 74 : 75 : TAB1.FLU.'EXPOSANT' : exposant q de la loi de propagation en 76 : *q 77 : fluage da/dt = A C . 78 : 79 : TAB1.CT.'EPS' : abscisses de la courbe de traction 80 : monotone du materiau. 81 : 82 : TAB1.CT.'SIG' : ordonnees de la courbe de traction 83 : monotone du materiau. 84 : 85 : TAB1.EPS.'COEFFSEC' : coefficient C de la loi de fluage 86 : n 87 : secondaireEpsf = 100 C sig . 88 : 89 : TAB1.EPS.'EXPSEC' : exposant n de la loi de fluage secondaire 90 : n 91 : Epsf = 100 C sig . 92 : 93 : * 94 : A et q sont entres pour un C exprime en MPam/h et un da/dt en m/h. 95 : C est entre pour les contraintes exprimees en MPa. 96 : 97 : * Si le materiau est l'acier 316 98 : 99 : TAB1.'A316' : VRAI (objet de type LOGIQUE) 100 : 101 : Cette option contient les lois de comportement necessaires au 102 : calcul en fatigue-fluage pour l'acier 316 a 525°C. On a tenu 103 : compte du fluage primaire et du fluage secondaire. 104 : 105 : 106 : 107 : TAB2 (type TABLE) : table contenant les donnees du chargement 108 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 109 : 110 : TAB2.'DP' : objet de type CHPOINT de forces appliquees 111 : a la structure pour la contribution primaire 112 : du chargement en fatigue. 113 : 114 : TAB2.'DC' : objet de type CHPOINT de forces appliquees a 115 : la structure pour le chargement complet de 116 : fatigue. 117 : 118 : TAB2.'PA' : objet de type CHPOINT de forces appliquees a 119 : la structure pour la contribution primaire du 120 : chargement A. 121 : 122 : TAB2.'PB' : objet de type CHPOINT de forces appliquees a 123 : la structure pour la contribution primaire du 124 : chargement B. 125 : 126 : TAB2.'CA' : objet de type CHPOINT de forces appliquees a 127 : la structure pour le chargement A complet. 128 : 129 : TAB2.'CB' : objet de type CHPOINT de forces appliquees a 130 : la structure pour le chargement B complet. 131 : 132 : TAB2.'DPTH' : objet de type CHPOINT de trmperatures 133 : appliquees a la structure pour la 134 : contribution primaire du chargement thermique. 135 : 136 : TAB2.'DCTH' : objet de type CHPOINT de temperatures 137 : appliquees a la structure pour le chargement 138 : thermique complet. 139 : 140 : TAB2.'TPSM' : objet de type FLOTTANT. Temps de maintien du 141 : chargement (B) en fluage. 142 : Il est entre en heure. 143 : 144 : Remarque : Selon l'Annexe A16 version2, le chargement de fatigue est 145 : defini par DP = PB - PA (DC = CB - CA). S'il existe un transitoire 146 : thermique T entre les etats permanents A et B alors 147 : DP = PB - PA + PT (DC = CB - CA + CT) (PT et CT sont les champs 148 : de forces associes respectivement aux champs DPTH et DCTH). 149 : 150 : 151 : AUTRES ARGUMENTS OBLIGATOIRES : 152 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 153 : 154 : TAB1.'MFIS' : objet de type MAILLAGE constitue par les 155 : elements LISP. 156 : 157 : TAB1.'MTOT' : objet de type MAILLAGE constitue par les 158 : autres elements. 159 : 160 : TAB1.'OBJ' : objet de MMODEL, modele associe au maillage 161 : constitue d'elements LISP. 162 : 163 : TAB1.'OBJT' : objet de MMODEL, modele associe au reste du 164 : maillage. 165 : 166 : TAB1.'MAT' : objet de MCHAML, materiau associe au maillage 167 : constitue d'elements LISP. 168 : 169 : TAB1.'MATT' : objet de MCHAML, materiau associe au reste du 170 : maillage. 171 : 172 : TAB1.'CAR' : objet de MCHAML, caracteristiques associees au 173 : maillage constitue d'elements LISP. 174 : 175 : TAB1.'CART' : objet de MCHAML, caracteristiques associees au 176 : reste du maillage. 177 : 178 : Remarque : TAB1.'CAR' contient la definition geometrique de la 179 : fissure au niveau de chaque element LISP. Son contenu evolue au 180 : cours du calcul en fatigue-fluage. 181 : Le maillage doit etre en mm. 182 : 183 : TAB1.'BLO' : objet de RIGIDITE contenant toutes les 184 : conditions limites de blocage pour effectuer 185 : une resolution. 186 : 187 : 188 : ARGUMENTS OPTIONNELS : 189 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 190 : 191 : TAB1.'REACTUA' : objet de type FLOTTANT, s'exprimant en 192 : pourcentage de la profondeur initialisee a la 193 : reactualisation precedente ( < 10% valeur par 194 : defaut). 195 : 196 : TAB1.'ALIMIT' : objet de type FLOTTANT, rapport a/w 197 : (a profondeur, w epaisseur de coque) 198 : maximum marquant la fin de calcul (0,7 par 199 : defaut). 200 : 201 : TAB1.'FATUNI' : objet de type LOGIQUE. Si vrai on ne calcule 202 : que la contribution de la propagation par 203 : fatigue (fatigue et fluage pris en compte par 204 : defaut). 205 : 206 : TAB1.'NBREA' : objet de type ENTIER representant le nombre 207 : de reactualisations. 208 : 209 : TAB1.'PROPA' : objet de type ENTIER permettant 210 : eventuellement en fluage primaire de changer 211 : le nombre de pas de calcul durant le temps de 212 : maintien ( par defaut 500). C'est le nombre 213 : maximum entre deux reactualisations. 214 : 215 : TAB1.'RELAX' : objet de type LOGIQUE. Si vrai, ceci permet de 216 : tenir compte de la relaxation des contraintes 217 : au cours du temps de maintien. 218 : 219 : Remarque : Pour le calcul de la contrainte de reference, on fait 220 : appel a la procedure @SIGREF. 221 : 222 : 223 : 224 : 2) Arguments de sortie : 225 : ------------------------ 226 : 227 : Ces arguments sont rajoutes dans TAB1. 228 : Les tables sont indicees par le numero de reactualisation. 229 : 230 : 231 : Nombre de cycles : 232 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 233 : 234 : TAB1.'EVNCY' : table contenant le nombre total N de cycles 235 : atteint a la nieme reactualisation 236 : (N = TAB1.'EVNCY'.n). 237 : 238 : TAB1.EVCAR.'TOTAL' : table contenant les objets caracteristiques 239 : associes a chaque reactualisation cumulant 240 : fatigue et fluage. 241 : 242 : TAB1.EVCAR.'FATIGUE' : table contenant seulement la contribution 243 : de fatigue. 244 : 245 : TAB1.EVCAR.'FLUAGE' : table contenant seulement la contribution 246 : de fluage. 247 : 248 : Remarque : Pour ces tables, n = 0 repere l'objet caracteristique 249 : initiale. 250 : 251 : 252 : Champs des efforts generalises issus des resolutions: 253 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 254 : 255 : TAB1.EVCON.'DELTSIG' : table contenant le champ des efforts issus 256 : du chargement complet de fatigue. 257 : 258 : TAB1.EVCON.'DELTSIG1': table contenant le champ des efforts issus 259 : du chargement primaire de fatigue. 260 : 261 : TAB1.EVCON.'CHSIG-A' : table contenant le champ des efforts issus 262 : du chargement complet A. 263 : 264 : TAB1.EVCON.'CHSIG-A1': table contenant le champ des efforts issus 265 : du chargement primaire A. 266 : 267 : TAB1.EVCON.'CHSIG-B' : table contenant le champ des efforts issus 268 : du chargement complet B (fluage). 269 : 270 : TAB1.EVCON.'CHSIG-B1': table contenant le champ des efforts issus 271 : du chargement primaire B (fluage). 272 : 273 : Remarque : Ces champs contiennent les efforts fournis par les s 274 : resolution successives avec les elements LISP definis par unite 275 : de longueur et par rapport a la fibre neutre du barreau sain 276 : (Nzz = sigm t ; Mxx = sigb t²/6 ; t est l'epaisseur du barreau), 277 : ainsi que le facteur d'intensite de contrainte KI. 278 : 279 : 280 : 281 : Champs des parametres principaux de l'Annexe A16 ( elements LISP): 282 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 283 : 284 : TAB1.EVCON.'SDEFRAS' : table des champs de contrainte reelle pres 285 : du defaut sous chargement A. 286 : 287 : TAB1.EVCON.'SDEFRBS' : table des champs de contrainte reelle pres 288 : du defaut sous chargement B. 289 : 290 : TAB1.EVCON.'SDEFR' : table des champs de contrainte reelle pres 291 : du defaut sous chargement de fatigue. 292 : 293 : TAB1.EVCON.'VALR' : table des champs de la valeur R, rapport des 294 : valeurs extremes que voit le defaut en chacun 295 : de ces points. 296 : 297 : TAB1.EVCON.'VALQ' : table des champs de la valeur Q : 298 : si R < 0 Q = (1 - 0,5 R) / (1 - R), 299 : si R > 0 Q = 1 / (1 - 0,5 R). 300 : 301 : TAB1.EVCON.'VALPSI' : table des champs de la valeur psi (correction 302 : de plasticite). 303 : 304 : TAB1.EVCON.'VALKDJ' : table des champs de la valeur kDJ(methode Js). 305 : 306 : TAB1.EVCON.'VALDJ' : table des champs de la valeur DJ (prenant en 307 : compte la plasticite). 308 : 309 : TAB1.EVCON.'VALDKEFF': table des champs de la valeur DKeff 310 : (variation de la valeur efficace du facteur 311 : d'intensite de contrainte corrige). 312 : 313 : TAB1.EVCON.'EPSILON' : table des champs de la valeur epsfluage 314 : cumulee en % (methode 4 du rapport 94/612). 315 : 316 : 317 : Listes de reels des valeurs maxi au cours du temps: 318 : ¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨ 319 : 320 : TAB1.'LR-NBCYT' : liste de reels des nombres de cycles totaux. 321 : 322 : 323 : TAB1.'LR-DAFAT' : liste de reels des Da de fatigue maxi. 324 : 325 : TAB1.'LR-DAFLU' : liste de reels des Da de fluage maxi. 326 : 327 : 328 : TAB1.'LR-ATOT1' : liste de reels de la profondeur maxi sur 329 : tous les elements LISP due a la contribution 330 : de fatigue. 331 : 332 : TAB1.'LR-ATOT2' : liste de reels de la profondeur maxi due a la 333 : contribution de fluage. 334 : 335 : TAB1.'LR-ATOT' : liste de reels de la profondeur maxi due aux 336 : contributions combinees de la fatigue et du 337 : fluage. 338 : 339 : TAB1.'LR-EPFLU' : liste de reels des epsfluage cumule maxi. 340 : 341 : TAB1.'LR-SEFB' : liste de reels de la contrainte sig(B)efficace 342 : maxi. 343 : 344 : 345 : 346 : Traces: 347 : ¨¨¨¨¨¨¨ 348 : 349 : EVCHSIGB : objet de type EVOLUTION, trace des contraintes 350 : efficaces maxi (en B). 351 : 352 : EVDAFLU : objet de type evolution, trace des da de 353 : fluage maxi en fonction des nombres de cycles 354 : totaux. 355 : 356 : EVDAATOT : objet de type evolution, trace des profondeurs 357 : maxi due aux contributions combinees de la 358 : fatigue et du fluage en fonction des nombre 359 : de cycles totaux. 360 : 361 : EVEPFLU : objet de type evolution, trace des epsfluage 362 : cumule maxi en fonction des nombres de cycles 363 : totaux. 364 : 365 : Ce sont des traces fournis par la procedure juste apes l'execution 366 : (obtenus en appelant DESS avec les listes precedentes). 367 : 368 : 369 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales