* fichier : rupt19.dgibi ************************************************************************ ************************************************************************ ******************************************************* * Test rupt19.dgibi: Jeux de données * * --------------------------------- * ******************************************************* * CAS TEST DU 15/12/15 PROVENANCE : TEST *Cas test de validation pour le calcul de J sous plusieurs chargement *avec les procedures g_theta.procedur et g_calcul.procedur * *- chargement en traction *- chargement avec pression sur levres *- chargement thermique *Calcul en dimension 2 avec des elements QUA8 sur un maillage complet *non symetrique * ig1 : mettre a vrai pour activer traces ig1 = faux ; ************************ *Données paramètriques : ************************ * a : profondeur de la fissure * * t : epaisseur du tube * * ri, re : rayon interne/externe * * h : hauteur du tube * h = 1. ; t = 60.e-3 ; a = t/5 ; ri = t*5; re = ri+t; *COORDONNEE DE LA POINTE DE LA FISSURE pf = (a + ri) 0. ; *NOMBRE D'ELEMENTS AUTOUR DE LA POINTE DE LA FISSURE (1 et 2 COUT) n_fiss = 10 ; *TAILLE D'UN ELEMENT DE LA 1ERE ET 2EME COUTURE* t_el = 200e-6 ; t_el2 = 400e-6 ; *Facteur d'agrandissement de la taille du derafinement tt_el2 = 4.*t_el2 ; *LONGUEUR DE LA 1ERE ET 2EME COUTURE* lc1 = n_fiss * t_el ; lc2 = t_el2 * n_fiss; *NIVEAU DE CHARGEMENT p0T = -400. ; p0P = 400. ; dt0 = 300.; *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * DEBUT DU MAILLAGE *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ************************************************************** ********************** 1ERE COUTURE ************************ ************ (Autour de la pointe de la fissure) *********** ************************************************************** ************************************************************** ********************** 2EME COUTURE ************************ ************ (Autour de la pointe de la fissure) *********** ************** Partie au-dessus de la fissure ************** ************************************************************** cout1et2 = cout1 et cout2; ************************************************************** ***************** DERAFINEMENT DES COUTURES **************** ************** Partie au-dessus de la fissure ************** ************************************************************** *------------------( DERAF A 4 ELEMENT )------------------- *------------------( DERAF A 3 ELEMENT )------------------- saa1 = sa1 ; repe i0 1 ; fin i0 ; sa1 = sa1 et ssa1 ; *---------------------- PARTIE GAUCHE ----------------------- *---------------------- PARTIE DROITE ----------------------- *---------------------- PARTIE HAUTE ----------------------- *lignes diagonales pour la symetrie *---------------------- PARTIE COIN ----------------------- cout_tot = cout1 et cout2 et cout3; ************************************************************** ********************* RESTE DU MAILLAGE ******************** ************** Partie au-dessus de la fissure ************** ************************************************************** *Partie de gauche *---------------- p_ri = ri 0.; pg = d_partg tran *Partie de droite *---------------- p_re = re 0.; pd = d_partd tran *Partie du haut *-------------- dini 1.6e-3 dfin t ; *Structure haute *--------------- *Structure basse *--------------- *Structure totale *************************************************************** *************************************************************** *DEFINTION DES LEVRES *lvsup = lvsup coul roug ; lvinf = lvinf coul bleu; *Definition des bords *bord haut *bord gauche *bord droit *bord bas si ig1 ; fins ; *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * FIN DU MAILLAGE *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * PARTIE CALCULS *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * PROPRIETE MATERIAUX A 300°C E0 = 185e3; nu0 = 0.3; alfa0 = 13.08e-6; *Elevation de temperature *CONDITIONS AUX LIMITES *Blocages cl0 = cl1 et cl2; *Traction uniaxiale (via un modele de pression) *Pression sur les levres (via un modele de pression) *********************** * SOLUTIONS ANALYTIQUES *********************** *Fonctions d'influence i0 = 1.211 ; i1 = 0.718 ; *Contraintes imposées pour le gradient de temperature sig0 = ((E0*alfa0*dt0)/(1-nu0)) * (ri/(3*t)) * ((2*(re**2))/(ri*(re+ri)) - 1); sig1 = -1. * ((E0*alfa0*dt0)/(1-nu0)); * J analytiques JT = (1-(nu0**2)) * ((i0*(-1.*p0T)*((pi*a)**(1./2.)))**2) / E0; JP = (1-(nu0**2)) * ((i0*(-1.*p0P)*((pi*a)**(1./2.)))**2) / E0; JTH = (1-(nu0**2)) * ((((i0*sig0)+ (i1*sig1*(a/t)))*((pi*a)**(1./2.))) **2)/E0; **************************************************** * CALCUL ELASTIQUE AVEC RESO - CALCUL DE J ELASTIQUE **************************************************** * Construction des second membres * RESOLUTION ELASTIQUE DES 3 PROBLEMES mena ; *PROCEDURE G_THETA *cas 1 : traction seule tabJel = table ; tabJel . 'MODELE' = mo0 et moph ; tabJel . 'CARACTERISTIQUES' = ma0 ; tabJel . 'PRESSION' = maph ; tabJel . 'BLOCAGES_MECANIQUES' = cl0 ; tabJel . 'SOLUTION_RESO' = utestT ; tabJel . 'LEVRE_SUPERIEURE' = lvsup ; tabJel . 'LEVRE_INFERIEURE' = lvinf ; tabJel . 'FRONT_FISSURE' = pf ; tabJel . 'COUCHE' = 5 ; g_theta tabJel ; JelT1 = tabJel.resultats ; *cas 2 : pression sur les levres tabJel . 'MODELE' = mo0 ET mopl; tabJel . 'CARACTERISTIQUES' = ma0 ; tabJel . 'PRESSION' = mapl ; tabJel . 'SOLUTION_RESO' = utestP ; g_theta tabJel ; JelP1 = tabJel.resultats ; *cas 3 : gradient de temperature tabJel . 'MODELE' = mo0 ; tabJel . 'CARACTERISTIQUES' = ma0 ; tabJel . 'TEMPERATURES' = cht0 ; tabJel . 'SOLUTION_RESO' = utestTH ; g_theta tabJel ; JelTH1 = tabJel.resultats ; *Erreurs sur J : solution analytique VS calcul RESO + G_THETA errT1 = ((JelT1-JT)/JT)*100.; errP1 = ((JelP1-JP)/JP)*100.; errTH1 = ((JelTH1-JTH)/JTH)*100.; ******************************************************* * CALCUL ELASTIQUE AVEC PASAPAS - CALCUL DE J ELASTIQUE ******************************************************* * Chargements de pression (obligatoires si modele de pression) * Chargement thermique *RESOLUTION AVEC PASAPAS DES 3 PROBLEMES (UN A CHAQUE PAS DE TEMPS) *AU PAS 1 : Traction seule *AU PAS 2 : Pression sur les levres *AU PAS 3 : Gradient de temperature tabT . 'MODELE' = mo0 ET moph ET mopl ; tabT . 'CARACTERISTIQUES' = ma0 ; tabT . 'BLOCAGES_MECANIQUES' = cl0 ; tabT . 'CHARGEMENT' = chaph ET chapl ET chath ; PASAPAS tabT ; *PROCEDURE G_THETA POUR LES 3 PROBLEMES (UN A CHAQUE PAS DE TEMPS) tabJel . 'SOLUTION_PASAPAS' = tabT ; tabJel . 'LEVRE_SUPERIEURE' = lvsup ; tabJel . 'LEVRE_INFERIEURE' = lvinf ; tabJel . 'FRONT_FISSURE' = pf ; tabJel . 'COUCHE' = 5 ; g_theta tabJel ; JelT2 = tabJel.resultats.1 ; JelP2 = tabJel.resultats.2 ; JelTH2 = tabJel.resultats.3 ; *Erreurs sur J : solution analytique VS calcul PASAPAS + G_THETA errT2 = ((JelT2-JT)/JT)*100.; errP2 = ((JelP2-JP)/JP)*100.; errTH2 = ((JelTH2-JTH)/JTH)*100.; **************************************** * AFFICHAGE DES RESULTATS ET DES ERREURS **************************************** mess 'Solution Theorique : ' JT JP JTH ; mess ; mess 'Erreur en % : ' errT1 errP1 errTH1 ; mess ; mess 'Erreur en % : ' errT2 errP2 errTH2 ; * Test sur les erreurs si ((abs errT) > 0.3) ; erre 'Erreur sur le calcul de JelT' ; fins ; si ((abs errP) > 0.3) ; erre 'Erreur sur le calcul de JelP' ; fins ; si ((abs errTH) > 0.4) ; erre 'Erreur sur le calcul de JelTH' ; fins ; FIN ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales