Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 1/9 Clé : U2.05.05 Révision : 12708 Calcul de structure en fatigue vibratoire 1 But Ce document a pour but de décrire la mise en œuvre d'un calcul de structure en fatigue vibratoire. Il s'agit plus précisément d'estimer l'amplitude maximale de vibration admissible d'une structure soumise à un chargement statique (connu) et à un chargement dynamique (inconnu). Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 2/9 Clé : U2.05.05 Révision : 12708 Table des Matières 1 But.......................................................................................................................................................1 2 Contexte industriel...............................................................................................................................3 3 Présentation de la méthode de calcul..................................................................................................3 3.1 Cadre du calcul et hypothèses......................................................................................................3 3.2 Calcul des vibrations admissibles.................................................................................................4 4 Mise en œuvre du calcul.....................................................................................................................5 4.1 Calcul des contraintes statiques et dynamiques...........................................................................5 4.2 Définition des propriétés matériau................................................................................................5 4.3 Appel à CALC_FATIGUE..............................................................................................................5 4.4 Interprétation des résultats............................................................................................................5 4.5 Remarques et conseils..................................................................................................................6 5 Exemple...............................................................................................................................................6 5.1 Description du cas de calcul.........................................................................................................6 5.2 Résultats obtenus..........................................................................................................................7 Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN 2 Date : 12/11/2014 Page : 3/9 Clé : U2.05.05 Révision : 12708 Contexte industriel Ce document a pour but de décrire la mise en œuvre d'un calcul de structure en fatigue vibratoire. On s'intéresse plus particulièrement aux études sur les machines tournantes (ailettes de turbines par exemple). Ces composants sont sollicités par un chargement statique (effort centrifuge lié à la rotation de la machine) et par un chargement dynamique (vibrations induites par l’écoulement du fluide). La durée de vie de la structure dépend à la fois de la partie dynamique (amplitude des variations des contraintes) et de la partie statique (contrainte moyenne). Le chargement statique est généralement bien connu ; la partie dynamique est quant à elle difficile tant à mesurer qu’à calculer. On ne peut donc pas mettre directement en œuvre les outils classiques de calcul de la durée de vie, basés sur l’évolution des contraintes au cours du temps (confer remarque ci-dessous). A l’inverse, il peut être intéressant d’estimer l’amplitude de variation maximale admissible par la structure, pour un ou plusieurs modes propres de vibration. Cette amplitude peut ensuite être comparée avec les vibrations mesurées sur site (mesures BVM par exemple). La démarche d’une telle étude est la suivante : •calcul de la contrainte liée au chargement statique ; •calcul de la contrainte associée au mode propre considéré ; •application des critères de fatigue pour calculer, en chaque nœud du maillage, la vibration maximale admissible pour avoir une endurance illimitée de la structure. La dernière étape peut se faire avec l'opérateur CALC_FATIGUE ( TYPE_CALCUL = 'FATIGUE_VIBR' ). Le principe du calcul, sa mise en œuvre dans Code_Aster et un exemple industriel sont décrits dans les chapitres suivants. Remarques : Un certain nombre d'hypothèses sont introduites pour la calcul, cf. paragraphe suivant. Il s'agit ici d'estimer, de manière conservative, un niveau vibratoire acceptable. Si 'évolution temporelle des contraintes dans la structure est connue, il faut utiliser les fonctionnalités classiques des opérateurs CALC_FATIGUE ou POST_FATIGUE (avec des critères uniaxiaux ou multiaxiaux, des méthodes de comptage de cycles, …). En dehors des turbines, la méthodologie présentée ici pourrait s'appliquer également à l’analyse vibratoire de lignes de tuyauterie (la contrainte moyenne correspondant alors à la pression interne), ou pour les éoliennes. Le calcul nécessite de connaître au préalable la contrainte à la rupture et la limite d’endurance du matériau. 3 Présentation de la méthode de calcul 3.1 Cadre du calcul et hypothèses On se place dans le cadre d'une décomposition de la contrainte dans la structure sur une base modale : ∞ total t= stat ∑ i imod cos i ti , i =1 stat la contrainte liée au chargement statique et imod les contraintes associées aux modes propres de la structure. i et i sont respectivement la pulsation (connue) et le déphasage en notant (inconnu) du mode i. Le calcul classique consisterait à évaluer le dommage, les contributions modales et les déphasages étant connus. Ici, réciproquement, on cherche à estimer les contributions modales maximales associées à une endurance illimitée de la structure. Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 4/9 Clé : U2.05.05 Révision : 12708 Il n'est pas possible d'identifier de manière générale la contribution maximale de chacun des modes. On est donc conduit à introduire trois hypothèses simplificatrices : critère de fatigue uniaxial ; imposition a priori du poids relatif des modes ; maximisation de l'amplitude. Critère de fatigue uniaxial : l'utilisation d'un critère de fatigue uniaxial (méthode de Wöhler) revient à supposer que les directions principales du chargement statique et du chargement dynamique sont les mêmes. Cette hypothèse semble licite pour les structures usuelles visées (ailettes, lignes de tuyauterie, …) ; elle induit un conservatisme sans doute excessif dans le cas général. Dans toute la suite, la notation correspondra à la norme de la contrainte (von Mises ou Tresca). A noter qu'une approche multiaxiale (avec des critères de type Crossland ou Dang Van) serait également possible. Ces critères nécessitent toutefois de connaître, en plus de la limite d’endurance en traction du matériau, la limite d’endurance en cisaillement pur. Cette donnée est souvent indisponible pour les matériaux des centrales, et la fonctionnalité aurait été alors difficilement utilisable. Poids relatif des modes : on suppose que le poids relatif des différents modes propres considérés est connu. Ce poids relatif peut par exemple être estimé à partir de mesures sur site. On ne considère également qu'un nombre limité de modes propres N (dans la pratique, on aura le plus souvent N N 4 ou 5). Autrement dit : total t= stat ∑ i imod cos i t i . i =1 Le coefficient est le paramètre que l'on cherche à calculer. Maximisation de l'amplitude : dans l'expression de la contrainte ci-dessus, les déphasages i sont S alt , définie comme la demi-amplitude de variation de la contrainte , inconnus. La contrainte alternée N S alt = ∑ i imod . Cette définition de S alt est conservative. est alors calculée comme suit : i=1 3.2 Calcul des vibrations admissibles Sous les hypothèses introduites ci-dessus, le calcul se résume à identifier le coefficient associé à une endurance illimitée de la structure. La contrainte statique stat correspond ici à une contrainte moyenne, généralement non nulle. Il est nécessaire de prendre en compte cette contrainte dans la courbe de fatigue de Wöhler. Cette prise en compte se fait classiquement à l'aide du diagramme de Haigh [R7.04.01], soit avec la droite de Goodman, soit avec la parabole de Gerber. En notant S l la limite d'endurance et S u la limite à la rupture du matériau, on a : •avec la droite de Goodman : S max alt stat =S l 1− Su , soit stat N i =S l 1− / ∑ i mod . S u i=1 •avec la parabole de Gerber : S Manuel d'utilisation max alt =S l 1− 2stat S 2 u , soit =S l 1− 2stat Su 2 N / ∑ i imod . i=1 Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 5/9 Clé : U2.05.05 Révision : 12708 Salt paraboledeGerber Diagrammede HAIGH droite de Goodman Su σm Figure 3.2-1: Diagramme de Haigh Ce calcul est fait pour tous les nœuds ou points de Gauss de la structure. Les zones où est le plus faible correspondent aux zones qui limitent la durée de vie de la structure. Pour passer du coefficient à l'amplitude de vibration admissible, une opération supplémentaire est à réaliser. Cette opération est décrite dans le § 6 . 4 Mise en œuvre du calcul 4.1 Calcul des contraintes statiques et dynamiques L'estimation du niveau vibratoire admissible nécessite le calcul au préalable des contraintes statiques et dynamiques : •le résultat correspondant au chargement statique (effort centrifuge, pression, …) peut être calculé avec MECA_STATIQUE ou STAT_NON_LINE . Un seul pas de temps doit être présent dans la structure de données résultat. Le calcul des contraintes stat à partir du champ de déplacement se fait à l'aide de la commande CALC_CHAMP (option SIEQ_ELGA ou SIEQ_ELNO selon que l'on souhaite calculer le dommage aux noeuds ou au points de Gauss). •Le calcul des N modes propres considérés peut se faire avec la commande CALC_MODES. Le i calcul des contraintes modales mod se fait à l'aide de la commande CALC_CHAMP (option SIEQ_ELGA ou SIEQ_ELNO selon que l'on souhaite calculer le dommage aux noeuds ou au points de Gauss). 4.2 Définition des propriétés matériau Deux paramètres matériau sont nécessaires pour le calcul : •la valeur de la limite à la rupture du matériau S u . Ce paramètre doit être introduit dans l'opérateur DEFI_MATERIAU [U4.43.01] (mot clé facteur RCCM, opérande Su). •la limite d'endurance du matériau S l . Ce paramètre correspond au premier point de la courbe de Wöhler (opérateur DEFI_MATERIAU , mot clé FATIGUE , opérande WOHLER ). 4.3 Appel à CALC_FATIGUE L'utilisation de CALC_FATIGUE (TYPE_CALCUL = 'FATIGUE_VIBR') nécessite de renseigner les paramètres suivants : •les contraintes statiques et modales préalablement calculées ; •les propriétés matériau ; •la liste des modes à considérer ( NUME_MODE ) et leur poids relatif ( FACT_PARTICI ), correspondant aux coefficients i 1≤i≤ N ; • le choix du mode de prise en compte des contraintes moyennes (Gerber ou Goodman : opérande CORR_SIGM_MOYE ) ; Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 6/9 Clé : U2.05.05 Révision : 12708 •le choix du lieu de calcul du dommage : aux nœuds ( OPTION= 'DOMA_ELNO_SIGM' ) ou aux points de Gauss ( OPTION= 'DOMA_ELGA_SIGM' ). 4.4 Interprétation des résultats En sortie de l'opérateur CALC_FATIGUE , on dispose d'un champ (aux nœuds ou aux points de Gauss) de la valeur admissible de : les zones où est le plus faible correspondent aux zones qui limitent la durée de vie de la structure. Pour passer du coefficient à l'amplitude de vibration admissible, une opération supplémentaire est à réaliser. Supposons par exemple qu'on s'intéresse à l'amplitude des déplacement en un point donné, qu'on notera ∂ u . Ce point peut correspondre par exemple à la position d'un capteur, où à la zone d'amplitude de vibration maximale. i On note u mod le déplacement au point d'intérêt associé au mode i. L'amplitude de vibration admissible au point d'intérêt est alors : N ∂ u =min ∑ i u imod i =1 Ce calcul est illustré dans l'exemple ci-dessous. La valeur minimale de peut être obtenue de trois manières différentes : •soit en post-traitant dans Aster le champ résultat ( POST_RELEVE_T, OPERATION='EXTREMA' ) ; •soit en visualisant le champ résultat ; •soit dans le fichier message. L'information est imprimée comme suit : ------------------------------------------------------Amplitude de vibration maximale admissible par la structure : 1.318438 ------------------------------------------------------------------- 4.5 Remarques et conseils Quelques remarques peuvent être formulées sur cette fonctionnalité : •Le résultat du calcul est nécessairement très sensible à la qualité du calcul des contraintes. L'utilisateur doit donc s'assurer que la finesse de son maillage est suffisante. L'utilisation d'éléments quadratiques est indispensable. •De même, il est conseillé de faire le calcul aux points de Gauss, l'interpolation pour le passage points de Gauss – nœuds étant source d'imprécision. •Les zones critiques correspondent le plus souvent aux singularités géométriques. L'utilisateur devra vérifier si la singularité correspond à une réalité physique ou à une discrétisation trop grossière. A noter que l'utilisation du remaillage (logiciel Homard, opérateur MACR_ADAP_MAIL ) ne permet pas de converger vers une solution stable non nulle si la géométrie initiale est discrétisée trop grossièrement. 5 Exemple 5.1 Description du cas de calcul Le post-traitement en fatigue vibratoire va être illustré sur le cas des ailettes des aubages 5-12 des turbines basse pression des centrales REP. Les ailettes sont fixées au disque par des doigts et des broches. Le diamètre extérieur des aubages est de l'ordre de 2 mètres. Les propriétés matériaux retenues pour cette étude sont les suivantes : - limite d'endurance : - limite à la rupture : Manuel d'utilisation S l =400 MPa S u =800 MPa . Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 7/9 Clé : U2.05.05 Révision : 12708 Le chargement statique correspond à la rotation de la roue à 1500 tours/ min . Trois des premiers modes correspondent respectivement à la flexion tangentielle ( f =100 Hz ), à la flexion axiale ( f =180 Hz ) et à la torsion ( f =360 Hz ) Remarque : Les propriétés matériau et les chargements utilisés ici ne correspondent pas au cas réel. – Figure 5.1-1: Maillage d'un paquet d'ailette de l'aubage 5-12 5.2 Résultats obtenus On s'intéresse ici à l'amplitude maximale de vibration admissible en tête d'ailette (groupe de nœuds 'NO_CAPT' , correspondant à la position sur site du capteur BVM), afin que la durée de vie dans le premier doigt de la première ailette du paquet d'aubage soit infinie (groupe de mailles appelé 'DOI_1_1'). - Résultats issus de CALC_FATIGUE : L'appel à CALC_FATIGUE, pour un calcul du dommage aux points de Gauss pour le premier mode de la structure, se fait comme suit : DMG_MOD1= CALC_FATIGUE( TYPE_CALCUL = 'FATIGUE_VIBR', OPTION = 'DOMA_ELGA_SIGM', HISTOIRE = _F( RESULTAT = PRECONT, MODE_MECA = MODE, NUME_MODE=1, FACT_PARTICI = 1,), DOMMAGE = 'WOHLER', CORR_SIGM_MOYE = 'GOODMAN', MATER = MATE, ) Puis, pour extraire le coefficient minimale sur la zone d'intérêt : Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 8/9 Clé : U2.05.05 Révision : 12708 IMPR_RESU(FORMAT='RESULTAT', RESU=_F(CHAM_GD=DMG_MOD1,GROUP_MA='DOI_1_1'),); Les résultats obtenus, pour les trois premiers modes de la structure et pour une combinaison entre les deux premiers modes sont donnés dans le tableau suivant. Mode 1 (flexion tangentielle) Mode 2 (flexion axiale) Contrainte statique maximale Contrainte modale maximale Mode 3 (torsion) Mode 1 + 0,5 Mode 2 685 MPa 13 MPa 51 MPa 27 MPa min Gerber 1.8 0.8 1.3 1.5 min 3.3 1.5 2.3 2.8 Goodman - Utilisation des résultats : Une opération supplémentaire est à réaliser pour passer du coefficient à l'amplitude de vibration admissible en un nœud donné. Le nœud considéré ici correspond à la position d'un des capteurs de déplacement . Cette opération peut par exemple être fait par des boucles python simples : # Paramètres du calcul (les mêmes que dans CALC_FATIGUE) fact_partici = [1., 0.5] nume_mode = [1,2] nbmode = len(nume_mode) # EXTRACTION DES DEPLACEMENTS DE LA TABLE TDEP = POST_RELEVE_T(ACTION=_F( GROUP_NO='NO_CAPT', INTITULE='Depl. modal', RESULTAT=MODE, NOM_CHAM='DEPL', TOUT_CMP='OUI', OPERATION='EXTRACTION',) ,) ; DX = [None]*nbmode DY = [None]*nbmode DZ = [None]*nbmode tdepl = TDEP.EXTR_TABLE() for i in range(nbmode) : dmod = tdepl.NUME_MODE==nume_mode[i] DX[i] = dmod['DX'].values()['DX'] DY[i] = dmod['DY'].values()['DY'] DZ[i] = dmod['DZ'].values()['DZ'] # EXTRACTION DU COEFFICIENT AMIN TDMG12=POST_RELEVE_T(ACTION=_F( GROUP_NO='ELTS', INTITULE='(MODE 1) + 0,5 * (MODE 2)', CHAM_GD=D_MOD12, TOUT_CMP='OUI', OPERATION='EXTREMA',) ,) ; AMIN = TDMG12.EXTR_TABLE() AMIN = AMIN.EXTREMA=='MIN' AMIN = AMIN['VALE'].values()['VALE'] Manuel d'utilisation Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) Version default Code_Aster Titre : Calcul de structure en fatigue vibratoire Responsable : Van Xuan TRAN Date : 12/11/2014 Page : 9/9 Clé : U2.05.05 Révision : 12708 # CALCUL DU DEPLACEMENT MAXIMAL ADMISSIBLE DXmax = 0. DYmax = 0. DZmax = 0. for i in range(nbmode) : DXmax = DXmax + AMIN[0] * fact_partici[i] * abs(DX[i][0]) DYmax = DYmax + AMIN[0] * fact_partici[i] * abs(DY[i][0]) DZmax = DZmax + AMIN[0] * fact_partici[i] * abs(DZ[i][0]) Dtot = sqrt(DXmax*DXmax + DYmax*DYmax +DZmax*DZmax ) . Cet exemple peut être retrouvé dans le cas test SDLV129A. Les résultats pour la présente étude sont donnés dans le tableau suivant. DX ( mm ) DY ( mm ) DZ ( mm ) Gerber DX max ( mm ) DY max ( mm ) DZ max ( mm ) D max ( mm ) Mode 1 -0.02 0.07 0.99 1.8 0.04 0.11 1.78 1.79 Mode 2 0.07 0.98 0.21 0.8 0.06 0.78 0.77 0.8 Mode 3 0.02 -0.18 -0.82 1.3 0.03 0.23 1.07 1.09 Mode 1 + 0,5 Mode 2 0.05 0.56 1.1 1.5 0.08 0.83 1.64 1.84 Cas Manuel d'utilisation min Fascicule u2.05 : Mécanique de la rupture et de l'endommagement Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html) ">
Public link updated
The public link to your chat has been updated.
