Commande POST_GP
Contents
1. 5 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 3 Calcul du param tre Gp Le calcul du param tre G n cessite un maillage bien particulier La zone situ e en aval de l entaille ou zone de propagation virtuelle de lentaille doit tre maill e en copeaux Chaque copeau doit tre identifi par un groupe de mailles et d crit dans la liste des copeaux en progressant partir du premier copeau n 1 situ en fond d entaille jusqu au dernier copeau TATANTA A A aA TAT a ATA TRA NA AAYA 5 PA VAN VAT AA AV ETAT PAPA PATATE TATA ra DN EEE ira Far NA A A TAVA Y LA VATA VA VAT AT A TA VA VA AT TA VV VA Pa A Va P Ex VA Aara A PA ratata YAYA VAVA YAYATA TAYA TAVAYA YATATA YATA TATATATA ta OTA ATA A PE VAVA VA TATAN VAVAVAVAVAYAVAVAVAVAVAVAVAYAVAVAYAVAVAVAVA VA AVAT Aa CA Z VA VAVAVA a A 7 o SR RA KW a EL eee EEK y TARER ERP OC 2RRERPEERIS Nu D A Pa A VAVA A A A A VAYA A A A AA A P RERCRERE KI SOOOCOCR CPC TAVAYA VA YAYA TAVAYA VA VAYAVAV TATAYAYATATATATAYA TAVAYA A Ya a C Pa ry AES A oN SAT a PATA VA YAVA VAYA A A YA TAVAYA YAYATA YAYATA TAYANA a a NA Aa A PA Va a A Va VA Va AA A YAYA A VA AAA AV TA A VA YAYA Va s A 4 A a SE AN A AA NN AAA DS A D OS A AAAA YAVAS Soa SAV TA A VA VA VAT AVATA TA VA VA VA TA PAPA A AA VA Va VA VA AA AP A A O TT A EN NW AW AY CRIER ER RTE OS K SPAVAVA YATA TAVAYA TAYA TA VAATA YA TAVA AVA YATA TAYA TATATA YATA YA2. comportement utilis e lors du calcul m canique Utilis par CALC _G 4 6 Op rande TYPE DEF 0 TYPE DEF YPETIT DEFAUT GRAND Indique si l utilisateur est en petites d formations ou en grandes d formations Dans le cas de petites d formations le code ne prend pas en compte les changement de g om trie et la surface du copeau dans le plan de sym trie n volue pas Par contre dans la seconde option les surfaces sont recalcul es chaque instant en tenant compte de la d formation des l ments 4 7 Op rande SYME CHAR SYME CHAR SANS DEFAUT SYME ANTI Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Ho Titre Commande POST_GP Date 06 07 2011 Page 8 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 Indique si le chargement est sym trique ou antisym trique dans le cas o on ne mod lise que la moiti du solide par rapport la fissure Dans ce cas un facteur multiplicatif de 2 appara t dans le calcul de G Utilis galement par CALC G 4 8 Op randes THETA 2D THETA 3D FOND FISS et DIRECTION On utilise THETA 2D ou THETA 3D suivant qu on travaille en 2D ou en 3D THETA 2D F GROUP NO fond entaille THETA 3D F GROUP MA fo
3. AX GPMAX GP CRIT PREDICTION EDICTION NUMI E TRANCHE La colonne PREDICTION vaut 0 tant que G est inf rieur au G critique 1 quand G atteint ou d passe G critique En 3D le r sultat DELTALMAX est une distance associ e au copeau Il s agit de la moyenne arithm tique des distances des noeuds situ s la fronti re interne du copeau au fond d entaille Autrement dit pour le copeau n la valeur DELTALMAX est la distance des noeuds communs entre le copeau n et le copeau n 1 au fond d entaille Pour le premier copeau la valeur DELTALMAX est nulle 5 Exemple d utilisation On trouvera des exemples et conseils d utilisation dans le document U2 05 08 dans le test 2D ssnp131 V6 03 131 et dans le test 3D ssnv207 V6 04 207 Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html
4. Code Aster Pt Titre Commande POST_GP Date 06 07 2011 Page 1 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 Macro commande POST _GP 1 But L objet de cette macro commande est de calculer le crit re nerg tique G l issue d un calcul thermo m canique dans les deux situations suivantes e identifier les valeurs critiques du param tre G en fonction de t nacit s critiques donn es une temp rature fix e e pr dire les instants de rupture sur un transitoire thermom canique partir de valeurs critiques de G pr c demment identifi s pour chaque temp rature Les diff rentes tapes accomplies par POST_GP sont e cr ation des champs theta avec CALC _ THETA e calcul de G avec CALC G e calcul de l nergie lastique avec POST_ELEM calcul de G f G E w La macro commande peut retourner trois tables e la premi re celle gauche du signe contenant les volutions de K G K noya Gus Pa e la deuxi me au mot cl TABL GPMAX contenant les r sultats de l identification ou de la pr diction La troisi me au mot cl TABL GP contenant les valeurs de G calcul s aux instants pr cis s G en fonction du temps Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster u Titre Commande POST_GP Date 06 07 2011
5. Le calcul est plus long lorsque cette option est activ e 416 Mot cl IDENTIFICATION KJ CRIT TEMP kj crit temp On renseigne avec ses deux mots cl s N couples kj crit temp valeur de K critique la temp rature donn e Pour chaque couple ki crit temp e on r cup re le module d Young et le coefficient de Poisson dans le mat riau fourni au mot cl MATER dont les valeurs permettent de calculer les valeurs des K i tant l indice de la couronne du champ theta partir des G e on calcule les moyennes K o CATE e d termination de G critique e en 2D on cherche quel instant K a atteint la valeur kj crit ce qui permet moy de trouver la valeur G critique G cet instant P max e en 3D on cherche quel instant et quel n ud sommet K a atteint la valeur moy kj_crit ce qui permet de trouver la valeur G critique G cet instant P max La table des r sultats est compos e de ces colonnes KJ CRIT INST GPMAX KGPMAX DELTALMAX en 3D LTALMAX NUME TRANCHI Auxquelles s ajoute la colonne NUME TRANCH KJ_CRIT INST GPMAX KGPMAX D nm m eai En 3D le r sultat DELTALMAX est une distance associ e au copeau ll s agit de la moyenne arithm tique des distances des noeuds situ s la fronti re interne du copeau au fond d entaille Autrement dit pour le copeau n la valeur DELTALMAX est la distance des noeuds communs entre le copea
6. Page 2 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 La macro commande fonctionne en 2D ou en 3D En 2D on utilise le mot cl THETA 2D en 3D THETA _ 3D Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Commande POST_GP Date 06 07 2011 Page 3 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 2 Syntaxe tab table POST GP RESULTAT resumeca resultat RESU THER resuther evol_ther MODELE modele modele MATER mater materiau COMP ELAS F 0 RELATION z ELAS DEFAUT ELAS VMIS LINE ELAS VMIS TRAC DEFORMATION A PETIT DEFAUT YGREEN TYPE DEF YBETIT DEFAUT GRAND 0 EXCIT a E CHARGE charg char_meca char_cine_meca FONC MULT fonet A fonction formule nappe Q TYPE CHARGE FIXE SYME CHAR SANS DEFAUT YSYME CANTI THETA 2D HU GROUP NO fond entaille 1 group no R INF LE inf R R SUP r sup R THETA 3D _F GROUP MA fond entaille 1 group ma R INF AT R R SUP x sup R FOND FISS fond entaille fond f
7. TA YAVAN S a KI KL EEE nn ne en ar AN AS ee A EAN ZA X V A V4 PAL AA DEN Re neue tend SA Ta VAA VATAYATAYA I MA copeau n 1 En 3D on introduit la notion de tranche chaque tranche contenant plusieurs copeaux O Te On range les copeaux dans une liste de groupes de mailles de la fa on suivante le 1er groupe correspond au copeau le plus proche du fond d entaille le 2 me groupe correspond au 1er copeau et au copeau adjacent le i me groupe correspond aux i copeaux les plus proches du fond d entaille Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl htmi Code Aster Fo Titre Commande POST_GP Date 06 07 2011 Page 6 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 _ En 3D e on d finit les copeaux de la 1 re tranche de la m me mani re qu en 2D la diff rence que les mailles sont volumiques et sont hexa driques e on poursuit la liste en ajoutant les copeaux de la 2 me tranche de la m me mani re e on obtient au final une liste de nb X ND panches groupes de mailles G est calcul partir de l nergie lastique sur la zone des copeaux situ e entre le copeau n 1 et le copeau n i et ceci pour i variant de 1 N N tant le num ro du dernier copeau cf U2 05 01 totale Am R o Am est la mesure du copeau au f
8. iss TRANCHE EX GROUP MA copeaux l group ma 0 DIRECTION dir LER Q CRIT_MAXI_GP YABSOLU DEFAUT RELATIF RAYON AXIS R R 1 DEFAUT TRAC COMP OUI DEFAUT M thode de discr tisation de th ta en fond de fissure 3D local 0 LISSAGI E EX LEGENDRE LAGRANGE LISSAGE THETA DEFAUT Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Commande POST _GP Date 06 07 2011 Page 4 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 LAGRANGE REGU LISSAGE G LEGENDRE DEFAUT _ LAGRANGE LAGRANGE NO NO LAGRANGE REGU DEGRE 0 Ho Je 2 T 3 4 53 DEFAUT 6 17 r IDENTIFICATION EX KJ CRIT kj crit R TEMP temp 1_R Jy PREDICTION 2 GP CRIT gp_crit R TEMP temp 1 _R Jg TABL GPMAX CO tablresu CO 0 TABL GP CO tabqgp CO LIST INST linst 1 R INFO 1 DEFAUT f 2 Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster ou Titre Commande POST_GP Date 06 07 2011 Page
9. les GROUP MA que de tranches Sous le mot cl GROUP MA les groupes de mailles sont d finis dans l ordre de Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster M Titre Commande POST_GP Date 06 07 2011 Page 9 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 proximit du fond d entaille Ainsi le premier groupe de mailles d fini dans GROUP MA sera le groupe de mailles le plus proche du fond d entaille Exemple dans le cas de 2 tranches avec Ti Cj correspondant la tranche et au j i me copeau tel que d fini au 3 et cf figure ci apr s TRANCHE F GROUP MA COPS 1 1 COPS 1 2 COPS 1 n GROUP MA COPS 2 1 COPS 2 2 COPS 2 n Frontd entaille COPS 1 1 COPS 2_1 COPS 1 2 C0PS 2 2 cops 1 3 COPS 2 3 4 10 Op rande CRIT MAXI GP CRIT MAXI GP YABSOLU DEFAUT YRELATIF D signe le crit re utilis pour extraire la valeur maximale de G e ABSOLU G valeur maximale de long de entaille e RELATIF G dernier pic atteint le long de l entaille 4 11 Op rande RAYON AXIS RAYON AXIS R D signe la valeur du rayon utilis dans la formule d Irwin en d formations planes _ GE NRI 4 12 Op rande TABL GPMAX TABL GPMAX CO tablresu tablresu est le nom dans le fichier de commande
10. nd entaille En 2D les mots cl s DIRECTION THETA 2D et FOND FISS sont utilis s par CALC THETA pour cr er le champ theta Le groupe fond entaille d signe le n ud du fond d entaille dont on r cup rera l volution de temp rature En 3D La d finition du fond d entaille dans DEFI FOND FISS est indispensable Le groupe de mailles d finissant le fond d entaille doit appartenir un plan et tre orient Celui ci doit tre d fini l aide de la commande NORMALE de DEFI FOND FISS Dans POST Gp e on utilise l op rande THETA 3D e la direction de propagation de l entaille est calcul e par la commande CALC G comme le produit vectoriel de l orientation donn e au front d entaille et de la normale d finie par le mot cl NORMALE telle que pr sent sur la figure ci dessous Attention donc la bonne orientation du fond d entaille et de la normale sous peine de r sultats faux Directions calul es YY US NZ Noeud final Normale Noeud initial 49 Op rande TRANCHE Ce mot cl d finit les groupes de maille composant les copeaux Exemple En 2D TRANCHE GROUP MA COPS 1l COPS 2 COPS 3 COPS 4 COPS 5 En 3D TRANCHE GROUP MA COPS 1 1 COPS 1 2 GROUP MA COPS 2 l1 COPS 2 2 r L ordre des copeaux est important le premier ensemble de groupes de mailles constituera la premi re tranche et ainsi de suite A fortiori il Y aura donc autant de groupes de mail
11. o commande et donc appelle en interne d autres commandes de Code Aster La plupart des mots cl s sont transmis tels quels aux autres commandes On indiquera par la suite dans quelle s commande s sont utilis s les mots cl s 4 1 Op rande RESULTAT D signe le r sultat du calcul thermo m canique pour lequel on calcule le param tre G Utilis par CALC_G et POST_ELEM 4 2 Op rande RESU THER D signe le r sultat du calcul thermique le cas ch ant Il est utilis pour extraire la temp rature en fond d entaille En l absence de ce mot cl on utilise la temp rature fournie au mot cl TEMP de IDENTIFICATION ou PREDICTION En IDENTIFICATION on boucle sur les diff rentes temp ratures fournies derri re TEMP et pour chacun des associ on en d duit un G critique En PREDICTION la liste GP_CRIT TEMP et la temp rature en fond d entaille issue de RESU THER permettent de d terminer le bon instant critique Remarque le lien entre RESU_THER et RESU MECA se fait sur NUME ORDRE et pas sur INST 4 3 Op rande MODELE D signe le mod le utilis lors du calcul m canique Utilis par CALC_THETA et par POST_ELEM 4 4 Op rande EXCIT D signe le chargement utilis lors du calcul m canique Utilis par CALC_G et la liste des charges est fournie POST_ELE z 4 5 Op rande COMP_ELAS D signe la relation de
12. ond d entaille e en 2D la mesure calcul e est la distance du copeau au fond d entaille e en 3D la mesure calcul e est la surface du copeau au fond d entaille R est la valeur fournie sous RAYON AXIS en 2D fact_syme vaut 2 en cas de sym trie 1 sinon G fact symex En 3D on calcule G sur chaque tranche puis on d termine quelle tranche maximise G La table tab produite par POST GP contient l volution des K G j 1 au nombre de couronnes du champ theta K G et G en fonction du temps moy moy En 3D G est calcul aux n uds du fond de fissure alors que G est calcul par tranche Pour obtenir une table r sultat contenant l ensemble des grandeurs on utilise la convention suivante pour chaque l ment on transf re G au 1er n ud sommet de l l ment et au n ud milieu dans le cas d un maillage quadratique sauf pour la derni re tranche o on transf re galement au dernier n ud sommet On ajoute une colonne NUME TRANCHE pour indiquer sur quelle tranche a t calcul le Gpmax Les r sultats de l identification ou de la pr diction sont fournis via le mot cl TABL GPMAX Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Hs Titre Commande POST_GP Date 06 07 2011 Page 7 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 4 Op randes POST GP est une macr
13. s de la table des r sultats de l identification ou de la pr diction Le contenu de la table est d taill ci apr s en fonction du type d analyse 4 13 Op rande TABL GP 0 TABL GP CO tabgp tabgp contient l volution de G en fonction du temps pour chaque copeau 4 14 Op randes LIST_INST et INST LIST INST linst Cet op rande permet de d finir une liste d instants auxquels les post traitements de la commande seront effectu s Pour ce faire la macro commande POST_GP utilise la commande EXTR RESU pour extraire des r sultats aux instants d finis L instant correspondant au num ro d ordre 0 est interdit Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Fo Titre Commande POST_GP Date 06 07 2011 Page 10 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 4 15 Mot cl TRAC COMP 0 TRAC _ COMP YOUL Cette option permet de prendre en compte la compression dans le calcul de l nergie pour chaque copeau l nergie dans chaque l ment est sign e par la trace de la contrainte si traction si compression et la sommation dans le copeau se fait avec les nergies sign es Si le r sultat est positif il est inchang s il est n gatif il est mis z ro En l absence de cette option on utilise POST ELEM pour calculer l nergie lastique de la liste des copeaux
14. u n et le copeau n 1 au fond d entaille Pour le premier copeau la valeur DELTALMAX est nulle L instant est obtenu par interpolation si celui ci ne correspond pas un instant de la liste d finie dans LIST_INST avec une pr cision relative de 0 01 Une alarme est alors mise afin que l utilisateur puisse juger de la pertinence de l interpolation le calcul est non lin aire 4 17 Mot cl PREDICTION gp_crit temp GP CRIT TEMP Manuel d utilisation Fascicule u4 82 M canique de la rupture Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Commande POST_GP Date 06 07 2011 Page 11 11 Responsable Renaud BARGELLINI Cl U4 82 31 R vision 6693 On renseigne avec ses deux mots cl s M couples gp crit temp Variation de G critique en fonction de la temp rature En 2D pour chaque instant du transitoire on peut ainsi valuer G critique la temp rature du fond d entaille et le comparer au G obtenu En 3D pour chaque instant du transitoire et pour chaque tranche du front d entaille on peut valuer G critique la temp rature du fond d entaille et le comparer au Ge obtenu La table des r sultats est compos e de ces colonnes NUME ORDRE INST TEMP DELTALMAX GPMAX GP CRIT PRI Auxquelles s ajoute la colonne NUME TRANCHE en 3D NUME ORDRE INST TEMP DELTALM
Download Pdf Manuals
Related Search