* fichier : rupt24.dgibi ************************************************************************ ************************************************************************ ******************************************************* * Test rupt24.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 3 avec des elements CUB8 sur un maillage complet *symetrique ************************ *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; *POINTS POUR L'AXE DE REVOLUTION p0 = (0. 0. 0.); py = (0. 1. 0.); *COORDONNEE DE LA POINTE DE LA FISSURE pf = (a + ri) 0. 0.; *NOMBRE D'ELEMENTS AUTOUR DE LA POINTE DE LA FISSURE (1 et 2 COUT) nfiss = 10; *TAILLE D'UN ELEMENT DE LA 1ERE ET 2EME COUTURE* tel = 200e-6 ; tel2 = 400e-6 ; *Facteur d'agrandissement de la taille du derafinement ttel2 = 4.*tel2 ; *LONGUEUR DE LA 1ERE ET 2EME COUTURE* lc1 = nfiss * tel ; lc2 = tel2 * nfiss; *NIVEAU DE CHARGEMENT p0T = -400. ; p0P = 400. ; dt0 = 300.; *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * DEBUT DU MAILLAGE *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ************************************************************** ********************** 1ERE COUTURE ************************ ************ (Autour de la pointe de la fissure) *********** ************************************************************** *cout1 = coul jaun cout1 ; ************************************************************** ********************** 2EME COUTURE ************************ ************ (Autour de la pointe de la fissure) *********** ************** Partie au-dessus de la fissure ************** ************************************************************** *cout2 = coul vert cout2 ; 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 ----------------------- couttot = cout1 et cout2 et cout3; ************************************************************** ********************* RESTE DU MAILLAGE ******************** ************** Partie au-dessus de la fissure ************** ************************************************************** *Partie de gauche *---------------- pri = ri 0. 0.; pg = d_partg tran *Partie de droite *---------------- pre = re 0. 0.; pd = dpartd tran *Partie du haut *-------------- dini 1.6e-3 dfin t; *Structure totale EN AXI ************************************* ************************************* *Definition des bords du maillage AXI n1 = 1; deg1 = 0.5 ; *Structure totale EN 3D ************************************************* ************************************************* * Creation des surfaces inferieure et superieure. *Defintion du front de fissure pfx = pfx1 0. pfx3; pfx et frfiss et f1 et f2) 1.e-5; *Trois points sur la surface f2 *============================================================= ************************************************************** *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ * 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 et cl3 et cl4 ; *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' = surlev ; tabJel . 'FRONT_FISSURE' = frfiss ; tabJel . 'COUCHE' = 5 ; g_theta tabJel ; JelT1 = tabJel.resultats.global ; *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.global ; *cas 3 : gradient de temperature tabJel . 'MODELE' = mo0 ; tabJel . 'CARACTERISTIQUES' = ma0 ; tabJel . 'TEMPERATURES' = cht0 ; tabJel . 'SOLUTION_RESO' = utestTH ; g_theta tabJel ; JelTH1 = tabJel.resultats.global ; *Erreurs sur J entre la solution analytique et le MEF 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 ; tabT . 'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ; PASAPAS tabT ; *PROCEDURE G_THETA POUR LES 3 PROBLEMES (UN A CHAQUE PAS DE TEMPS) *ATTENTION, IL FAUT RETIRER LE CHARGERMENT MECA DE PRESSION SUR LES *LEVRES ET UTILISER LE CHARGEMENT PLEV *PROCEDURE G_THETA POUR LES 3 PROBLEMES (UN A CHAQUE PAS DE TEMPS) tabJel . 'SOLUTION_PASAPAS' = tabT ; tabJel . 'LEVRE_SUPERIEURE' = surlev ; tabJel . 'FRONT_FISSURE' = frfiss ; tabJel . 'COUCHE' = 5 ; g_theta tabJel ; JelT2 = tabJel.resultats. 1 . global ; JelP2 = tabJel.resultats. 2 . global ; JelTH2 = tabJel.resultats. 3 . global ; *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.5) ; erre 'Erreur sur le calcul de JelT' ; fins ; si ((abs errP) > 0.5) ; erre 'Erreur sur le calcul de JelP' ; fins ; si ((abs errTH) > 0.2) ; erre 'Erreur sur le calcul de JelTH' ; fins ; FIN ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales