[U7.03.41] Procédure MACR_ECREVISSE
Contents
1. Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 GROUP MA gma grma GROUP NO ORIG ogno grno GROUP NO EXTR egno grno SECTION ELLIPSE RECTANGLE Q LISTE COTES BL 0 max abs curv DEFAUT lcbl listr8 LISTE VAL BL lvbl listr8 ZETA zeta R RUGOSITE eps R OUVERT REMANENTE ouv rem R Q TORTUOSITE 1 DEFAUT tort R PREF IXE FICHIER FISSUREl DEFAUT prefix Kn DONNEES RELATIVES A L ECOULEMENT ECOULEMENT F Z PR S ENTREE pe R PRES ENTREE FO fpe fonction nappe formule Z PRES SORTIE ps R PRES SORTIE FO fps fonction nappe formule FLUIDE ENTREE 1 2 18 4 or 6 si FLUIDE ENTREE 1 3 4 ou 6 TEMP_ ENTREE te R TEMP ENTREE FO fte fonction nappe formule si FLUIDE ENTREE 2 ou 5 TITR MASS xe R TITR MASS FO fxe fonction nappe formule si FLUIDE ENTREE 4 ou 5 PRES PART pae R PRES PART FO fpae fonction nappe formule CHOIX DES MODELES POUR ECREVISSE MODELE ECRE F Q IVENAC 0 DEFAUT 1 ECOULEMENT SATURATION GELE si ECOULEMENT GELE PRESS EBULLITION peb R FROTTEMENT gt 4 3 ELA 12 Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version de2. chaque instant l autre r capitulant les donn es concernant les fissures En cr ant un r pertoire de sortie dans le profil d tude on pourra y r cup rer toutes les sorties concernant Ecrevisse Il est possible d effectuer des poursuites du calcul Remarque ll est n cessaire de pr ciser le mot cl DEBUG F HIST ETAPE OUl dans DEBUT ou POURSUITE pour utiliser MACR ECREVISSE Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ut Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 2 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Table des Mati res 3 LOp randes TABLE TEMPER D BIT qa 6 3 2 Mot cle ETAT INT uuu ak ns nn un 6 3 3 Op ranoes MODELE MEGA MODELE THBR hu eue unless sense 6 3 4 Op rande LIST INST uu uuu usu a nee usant fabte dep tsapea dan dhanaer data sas abat don ide oh hetenisteto te 6 3 5 Op rande CHAM MATE Rossia TR ia a aaa nn Qh ptsassaathhata kusa 7 3 0 Op rande CARATJELENES 8 sasata taian F S M tle EXC MEGA sna ndadani aaa haana aaa iaaa kninska 7 sS MOTOS CONTACT Rd do ipa punin den nt T 29 Motele EXC TIHEPR uu uuu uuu iaaiiai 7 SO MOC CIE GOME INCR nr nn nn ainiaan al uni iegnamel ne lames tenu 7 S MOr e G
3. l paisseur de la paroi travers e Ce coefficient qui par d faut vaut 1 permet de ne f L prendre en compte ce ph nom ne La longueur de la fissure sera donc gale er avec tort lt 1 or 3 14 9 Op rande PREFIXE FICHIER Q PREFIXE FICHIER prefix Permet de d finir le nom des fichiers de sortie Ecrevisse dans le dossier REPE OUT courbes et fichiers textes et ainsi de distinguer les fichiers associ s chaque fissure Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Po Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 11 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 3 15 Mot cl ECOULEMENT 3 15 1 3 15 2 3 15 3 3 15 4 3 15 5 ECOULEMENT Il s agit dans ce mot cl de d finir la nature du fluide et de ses caract ristiques Op randes PRES ENTREE PRES ENTREE FO PRES ENTREE pe PRES ENTREE FO fpe Cet op rande permet de donner la pression totale Pa ou son volution au cours du temps en amont de la fissure Elle doit tre comprise entre 215 10 Pa et 10 Pa Op randes PRES SORTIE PRES SORTIE FO PRES SORTIE pe PRES SORTIE FO fpe Cet op rande permet de donner la pression totale Pa ou son volution au cours
4. rature en bord de fissure Entr e dans MACR ECRE CALC py Ex cution d ECREVISSE R cup ration du d bit calcul et des tables de r sultats Extraction et copie des donn es sortantes CL thermique r sultat Ecrevisse valuation du crit re de convergence erreur en temp rature et ou pression Red coupage ventuel du pas de temps voire arr t suivant les conditions R cup ration des infos et concepts chargement Le d tail des macro commandes interm diaires CALC ECREVISSE et MACR ECRE CALC est donn en annexe Unit s On peut habituellement se passer des unit s dans Code Aster en restant coh rent Cependant pour les besoins d Ecrevisse il est important ici de respecter le syst me d unit s international SI en particulier pour les longueurs m les pressions Pa et les temp ratures C Quelques conseils d utilisation Le maillage doit tre en 2 dimensions et la fissure maill e explicitement Pour que le calcul soit de bonne qualit il est n cessaire d optimiser le maillage pour qu il soit compatible avec les ph nom nes que l on veut repr senter les chargements m caniques les chargements thermiques mais aussi l coulement du fluide Il est ainsi important de raffiner suffisamment le maillage aux bords de la fissure si on veut voir s chauffer le mat riau En entr e de fissures les ph nom nes
5. 2008 03647 FR Modifications apport es au logiciel ECREVISSE version 3 1 pour obtenir ECREVISSE la version 3 1 1 HERVOUET Chantal 3 H 181 2009 00243 FR Simulation d coulement air eau vapeur travers une fissure logiciel ECREVISSE version 3 1 1 Manuel utilisateur HERVOUET Chantal Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 17 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 7 Annexes 7 1 Syntaxe de la proc dure CALC ECREVISSE CALC ECREVISSE CONCEPT SORTANT CHARGE MECA CO char meca Co CHARGE THERI CO char ther CO CHARGE THER2 CO char ther CO TABLE CO table CO DEBIT CO debit CO MODELES MECANIQUES MODELE MECA m meca modele M
6. T Ti Ter _ Max P Pi ONE s 2 2 e Ver ep Ensuite il peut activer s il le souhaite le red coupage du pas de temps si la valeur du crit re est inf rieure une valeur donn e L exp rience acquise montre toutefois qu il vaut mieux essayer d optimiser la liste d instant que de compter sur le red coupage du pas de temps car les calculs Ecrevisse sont relativement longs 3 17 1 Op rande CRITERE Permet de d finir la nature du crit re utilis pour g rer les pas de temps dans la macro commande Les valeurs possibles sont EXPLICITE pas de red coupage quelque soit la valeur des indicateurs TEMP il y aura red coupage en fonction de la valeur de e PRES il y aura red coupage en fonction de la valeur de ep TEMP_ PRES il y aura red coupage en fonction de la valeur de ec 3 17 2 Op rande TEMP_REF Temp rature de r f rence pour le calcul du crit re en temp rature de la macro commande 3 17 3 Op rande PRES REF Pression de r f rence pour le calcul du crit re en pression de la macro commande 3 17 4 Op rande PREC CRIT Valeur laquelle on compare l erreur obtenue pour activer ou non le red coupage du pas de temps 1 par d faut Par exemple pour la temp rature e lt PREC CRIT 3 17 5 Op rande NUME ORDRE MIN Num ro d ordre partir duquel le crit re d erreur est pris en compte La convergence est forc e sur les num ros d ordre inf rieurs On utilise souvent c
7. sortie sont facultatifs il est n anmoins vivement conseill de renseigner a minima TEMPER sinon il est impossible de contr ler les champs de temp rature dans la structure et de faire une poursuite sile concept est r entrant le mot cl TEMPER ne doit pas tre renseign ici mais sous le mot cl ETAT_INIT op rande EVOL THER 3 2 Motcl ETAT _INIT Q ETAT INIT F EVOL NOLI evol noli EVOL THER evol ther NUME_ORDRE nume ordre Mot cl permettant de d finir l tat initial dans le cadre d une poursuite de la macro commande 3 2 1 Op rande EVOL NOLI Nom du concept de type evol_noli d o sera extrait l tat m canique 3 2 2 Op rande EVOL THER Nom du concept de type evol ther d o sera extrait l tat thermique Ce concept sera enrichi au cours du calcul 3 2 3 Op rande NUME ORDRE L tat initial pour poursuivre le calcul sera d fini partir du num ro d archivage NUME ORDRE pour evol noli etpourevol ther 3 3 Op randes MODELE MECA MODELE THER MODELE MECA m meca MODELE THER m ther Ces mots cl s permettent de renseigner le nom du mod le m meca dont les l ments font l objet du calcul m canique et le nom du mod le m ther dont les l ments font l objet du calcul thermique 3 4 Op rande LIST_INST LIST INST list inst Manuel d utilisation Fascicule u7 03 Chainage avec un autre lo
8. veut ou non prendre en compte le transfert de chaleur par convection entre le fluide et la paroi Les valeurs valides sont 12 11 2 1 0 1 2 11 12 La valeur 0 correspond un calcul sans transfert de chaleur adiabatique Les autres valeurs correspondent diff rentes options pour le calcul du coefficient de convection diff rences uniquement pour le r gime laminaire Se reporter la documentation Ecrevisse 3 16 7 Op randes XMINCH XMAXCH Cet op rateur n a de sens que dans les cas diphasiques si TRANSFERT CHAL lt 0 Il permet de renseigner la valeur des titres massiques en gaz d limitant la zone de transition entre les mod les de calcul du coefficient de convection cf documentation Ecrevisse Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Ma Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 13 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 3 17 Mot cl CONV_CRITERE La macro commande calcule chaque instant trois crit res permettant d estimer l importance des changements intervenus entre deux pas de temps et ainsi la validit du cha nage effectu Pour calculer ces crit res l utilisateur doit d finir des valeurs de r f rence pour la pression et la temp rature permettant de quantifier l cart acceptable _ Max
9. A m meca modele_sdaster MODELE THER m ther modele sdaster LISTE D INSTANTS LIST INST list inst listr8 DONNEES POUR STAT NON LINE ET THER LINEAIRE CHAM MATER chmat cham mater Q CARA ELEM carac cara elem EXCIT MECA ONE CHARGE chi char meca Q FONC MULT fi fonction Q TYPE CHARGE FIXE CSTE DEFAUT SUIV DIDI CONTACT char contact char contact EXCIT THER F CHARGE chi char ther Q FONC MULT fi fonction COMP INCR F voir le document U4 51 11 Q NEWTON F voir le document U4 51 03 CONVERGENCE F voir le document U4 51 03 Q PARM THETA theta DEFAUT 0 57 R DONNEES GEOMETRIQUES RELATIVES A LA FISSURE FISSURE HE au Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Proc dure MACR ECREVISSE Date 22 04 2013 Page 4 20
10. Code Aster Pt Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 1 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Proc dure MACR ECREVISSE R sum La macro commande MACR ECREVISSE a pour but de calculer le d bit de fluide air vapeur liquide travers une ou plusieurs fissure s traversante s dans une structure mod lis e en 2 dimensions Elle r alise pour cela chaque pas de temps le cha nage de deux codes Code Aster qui permet de conna tre l tat thermo m canique de la structure Ecrevisse qui r alise le calcul thermo hydraulique d coulement travers la fissure En pratique la macro commande se charge d effectuer successivement pour tous les pas de temps sun calcul thermique lin aire eun calcul m canique quasi statique d appeler la macro commande CALC ECREVISSE qui elle va tre en charge 1 de r cup rer le profil de la fissure 2 d appeler une troisi me macro commande MACR ECRE CALC qui va g n rer le fichier de commande et lancer Ecrevisse 3 d extraire et copier les r sultats 4 de v rifier les crit res impos s sont v rifi s et si besoin red couper le pas de temps Le concept de sortie principal est la structure de donn es r sultat du calcul m canique evol noli Il est galement possible d obtenir la structure de donn es r sultats du calcul thermique ainsi que 2 tables l une contenant les d bits
11. ISTE VAL TEMP lvt R GRADIENT IMPOSE TEMP1 tml R Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Mo Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 19 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 TEMP2 tm2 R GRADIENT CALCULE EPAISSEUR PAROI epp R CONVECTION AMONT alphe R CONVECTION AVAL alphs R LAMBDA lambd R TEMP FLUIDE AVAL ts R DONNEES RELATIVES A L ECOULEMENT ECOULEMENT idem MACR ECREVISSE ECOULEMENT MODELE ECRE idem MACR ECREVISSE MODELE ECRE DONNEES RELATIVES A LA CONVERGENCE NUMERIQUE Q CONVERGENCE ECREVISSE idem MACR ECREVISSE CONVERGENCE ECREVISSE GENERAL COURBES LOGICIEL VERSION idem MACR ECREVISSE GENERAL 7 4 Op randes sp cifiques MACR ECRE CALC 7 41 Mot cl FISSURE 7 4 1 1 Op rande ANGLE Angle de l coulement par rapport l ascendante verticale degr s 7 4 1 2 Op rande LONGUEUR D finit la longueur r elle de la fissure tenant compte de sa tortuosit 7 4 1 3 Op
12. NEVOIONDQ pussuasaasganaqanaqaqaqapaqashusqahaqaqanuqhhia qhaquna yqappbiqanaspakapanhasiphhuqaqaqanpanapanqhythsqaqi 7 SIA Net ele CONVERGENCE uu uu s 0 eiiaoe aaa 8 3 13 Op rande PARM THETA Uhr Ne en ainia 8 S MOCC FISSURE uu oiia asa daaa a aa aaa nasua pu a Supa matin ennemie nues 8 3 15 Mot ele ECOUEEMENT 2 225 usmashuaaashiykutylaqukakababaqasa kiy caca tenez aimant nat ann none 10 3 16 Mot CE MODELE MECRE EUR Mie tatap upas Gene dt se taie seb 11 SMOE CONS CRITTER Eisra no us nn 11 3 18 M t CIE CONVERGENCE EGPREVISSE ciirsa antaas biaia 12 SA9 Op rande COURBE S mitih iiidid obada ai aaiae hanhia raa eide 13 320 Op ra de LOGICIEL rrna aaa a abaia 13 3 21 Op rande VERSION ioiiomiiiionniniibiiiiaaaimii iaiia iiaa ii 13 4 Foncdionnement du COMPAQ run nine nayasa aa a 13 T Pre Gun u Qasa e a a e aaa eee ne umumunu qas 13 A ZUNO CR Ne 14 4 5 Quelques Conseils d UGISAUOM 2222 22 naaid tete Pa ot ne entend ardent nee an 14 5 Exemple d atiliSat miuu u uuu uuu ieaiaia ida aaa E aaa 14 D Reference ECES SE icenian nianna ahdida aeiaai ki araea h ES 15 ANNONO S eina ta a aa e a E en 16 7 2 Op randes sp cifiques CALC_ECREVISSE ausssssssssssawasasasssasamaqssasssss 17 1 3 Syntaxe d la proc dure MACR ECRE CALC oiran 17 7 4 Op randes sp cifiques MACR_ECRE_CALC ssssssssssssssssssssasssssssssssss
13. ODELE THER m ther modele DONNEES GEOMETRIQUES RELATIVES AUX RESULTATS RESULTAT SFe MECANIQUE rmeca resultat THERMIQUE rther resultat NUME ORDRE nume_ordre I Q INST instant R DONNEES GEOMETRIQUES RELATIVES A LA FISSURE FISSURE ET GROUP MA gma grma GROUP NO ORIG ogno grno GROUP NO EXTR egno grno ZETA zeta R RUGOSITE rug R OUVERT_REMANENTE ouv rem R Q PREF IXE FICHIER FISSURE1 DEFAUT prefix Q TORTUOSITE tort R SECTION ELLIPSE RECTANGLE Q LISTE COTES BL 0 max abs curv DEFAUT Icbl listr8 LISTE VAL BL lvbl listr8 DONNEES RELATIVES A L ECOULEMENT ECOULEMENT idem MACR ECREVISSE ECOULEMENT MODELE ECRE idem MACR ECREVISSE MODELE ECRE DONNEES RELATIVES A LA CONVERGENCE NUMERIQUE Q CONVERGENCE ECREVISSE idem MACR ECREVISSE CONVERGENCE ECREVISSE GENERAL COURBES LOGICIEL VERSION idem MACR ECREVISSE GENERAL INFO Z 1 DEFAUT Z 2 3 Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster nn Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 18 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 72 Op randes sp cifiques CALC ECREVISSE 7 2 1 Op rande CHARGE MECA Concept sortant contenant le chargement m canique fourni pa
14. Si LISTE COTES BL n est pas renseign ISTE VAL BL doit contenir la largeur de la fissure l entr e et la sortie de la fissure 3 14 5 Op rande RUGOSITE RUGOSITE eps Rugosit absolue de la paroi en m tres Cette valeur n influe le r sultat que si le frottement est calcul par Ecrevisse op rande FROTTEMENT gt 0 et si l coulement n est plus laminaire 3 14 6 Op rande ZETA ZETA zeta Coefficient de perte de charge singuli re l entr e sans dimension 3 14 7 Op rande OUVERT REMANENTE OUVERT REMANENTE ouv rem Valeur de l ouverture r manente en m tres pour le calcul Ecrevisse correspondant l ouverture hydraulique r elle de la fissure lorsque les deux l vres sont au contact Ceci traduit le fait que m me lorsque la fissure est referm e une petite quantit de fluide transitera en particulier cause de la rugosit Ce param tre d pend du mat riau par exemple c est de l ordre de 10 u m pour le b ton Ce mot cl est donc pris en compte donc lors du calcul hydraulique Ecrevisse mais non par Code Aster il peut donc bien arriver que le calcul hydraulique soit effectu avec l ouverture r manente est que le contact soit activ dans le calcul m canique 3 14 8 Op rande TORTUOSITE Q TORTUOSITE tort Si la fissure est tortueuse on peut admettre que vue du fluide la longueur de la fissure est sup rieure
15. TION AMONT Coefficient de convection la surface de la paroi cot amont 7 4 2 8 Op rande CONVECTION AVAL Coefficient de convection la surface de la paroi cot aval 7 4 2 9 Op rande LAMBDA Conduction thermique de la paroi 7 4 2 10 Op rande TEMP_FLUIDE_AVAL Temp rature du fluide cot aval Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html
16. cl est identique celui de STAT NON LINE il convient donc de se reporter au document U4 51 03 3 12 Mot cl CONVERGENCE Il s agit ici de renseigner les crit res de convergence pour le calcul m canique La syntaxe correspond celle de STAT_ NON LINE voir le document U4 51 03 Remarque il peut tre n cessaire d utiliser RESI GLOB MAXI plut t que RESI GLOB RELA lorsqu il n y a pas de chargement m canique sinon effort nul donc RESI GLOB RELA ind fini 3 13 Op rande PARM THETA Il s agit ici de renseigner la valeur du param tre caract ristique du theta sch ma utilis pour r soudre le probl me thermique Voir le document U4 54 011 3 14 Mot cl FISSURE FISSURE Ce mot cl permet de d finir tous les param tres ayant trait une fissure Plusieurs occurrences de ce mot cl sont possibles s il y a plusieurs fissures 3 14 1 Op rande GROUP _ MA GROUP MA gma Groupes des mailles d finissant la ou les l vres de la fissure On donne le groupe de la l vre inf rieure et sup rieure 3 14 2 Op randes GROUP NO ORIG et GROUP NO EXTR GROUP NO ORIG ogno GROUP NO EXTR egno Permet de d finir les deux groupes de noeuds qui d finissent les deux extr mit s de la fissures afin de l orienter 3 14 3 Op rande SECTION SECTION ELLIPSE RECTANGLE Il s agit ici de d finir la section de passage de l coulem
17. du temps en sortie de la fissure Elle doit tre comprise entre 215 10 Pa et 10 Pa et tre sup rieure la pression d entr e sinon il n y a pas d coulement et le calcul s arr te Op rande FLUIDE ENTREE FLUIDE ENTREE 1 2 3 4 5 o 6 Cet op rande permet de d finir l tat du fluide l entr e de la fissure 1 Eau sous refroidie ou satur e Fluide diphasique Vapeur satur e ou surchauff e Air vapeur surchauff e Air vapeur satur e Air seul En fonction de la situation du fluide seules certaines des caract ristiques suivantes sont renseigner O O1 N Op randes TEMP_ ENTREE TEMP_ ENTREE FO TEMP ENTREE te TEMP ENTREE FO fte Cet op rande permet de renseigner la temp rature C ou son volution au cours du temps en amont de la fissure Il doit tre renseign si le fluide l entr e est lt Eau sous refroidie ou satur e Vapeur satur e ou surchauff e Air vapeur surchauff e ou Air seul op rande FLUIDE ENTREE 1 3 4 ou 6 Op randes TITR MASS TITR MASS FO TITR MASS xe TITR MASS FO fxe Cet op rande permet de renseigner le titre massique de vapeur ou son volution au cours du temps en amont de la fissure Il correspond au rapport entre la masse de vapeur et la masse d eau liquide et vapeur Il doit tre renseign si le fluide l entr e est lt F
18. ent plan perpendiculaire au plan de mod lisation Celle ci peut tre elliptique ou rectangulaire L une des dimensions de cette section correspond la taille de la fissure dans le plan de la mod lisation Celle ci est d termin e chaque instant par la macro commande L autre dimension correspond la taille de la fissure hors du plan de mod lisation Celle ci est fixe pendant tout le calcul et doit tre renseign gr ce aux op randes LISTE_COTES_BL et LIST_VAL BL 3 14 4 Op randes LISTE COTES BL LISTE VAL BL LISTE COTES BL lcb LISTE VAL BL 1vb Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Cod e A ster Version default Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 10 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 ouverture fissure un gt TEK A p direction de l coulement bi liste de valeurs Om abscisse curviligne liste de cotes Figure 3 14 4 a Repr sentation de la fissure Pour caract riser la taille de la fissure dans le plan non mod lis il faut fournir la valeur dux petit axe de la section elliptique ou la largeur de la section rectangulaire au minimum en 2 points ISTE COTES BL permet de renseigner les abscisses curvilignes o l on donne la dimension de la fissure et LIST VAL BL permet de donner la dimension
19. ette op rande pour d marrer l activation du crit re d erreur partir du deuxi me pas de temps et forcer la convergence sur le premier pas qui est souvent d licat traiter puisqu on impose des chargements brutaux l instant initial 3 17 6 Op rande SUBD NIVEAU Nombre de red coupage du pas de temps autoris Au del de cette valeur on ne red coupe plus et on sort de la macro commande 3 17 7 Op rande SUBD PAS MINI Pas de temps en dessous duquel on ne d coupe plus 3 18 Mot cl CONVERGENCE ECREVISSE 3 18 1 Op rande KGTEST Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster u Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 14 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Param tre de l algorithme it ratif d Ecrevisse Doit tre compris entre 0 et 1 Fix par d faut 0 5 on peut tre amen mettre une valeur l g rement inf rieure si des difficult s de convergence sur le d bit sont observ es 3 18 2 Op rande ITER GLOB MAXI Nombre maximum d it rations autoris pour le calcul du d bit 400 par d faut 3 18 3 Op rande CRIT CONV DEBI Pr cision utilis e pour la convergence du calcul du d bit C est la valeur par rapport laquelle Ecrevisse teste le d bit max et min dans son proc d de calcul 6 Cni l
20. fault Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 5 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 11 4 J ey 2 f l 0 144 19 173 4 O AE Z lt 1 27 Z 14 si FROTTEMENT n gatif REYNOLDS LIM relim R FROTTEMENT LIM frtlim R TRANSFERT CHAL f ST 11 S Z ely 0 1 Z 2 LT Z T2 si TRANSFERT CHAL n gatif XMINCH xminch R XMAXCH xmaxch R CRITERE DE CONVERGENCE DE LA MACRO COMMANDE CONV CRITERE F TEMP REF tref R PRES REF pref R CRITERE TEMP PRESS DEFAUT TEMP PRESS EXPLICITE si CRITERE diff rent de EXPLICITE Q NUME ORDRE MIN n I PREC CRIT prec R Q SUBD NIVEAU nb I Q SUBD PAS MINI nbpas R CONVERGENCE NUMERIQUE D ECREVISSE Q CONVERGENCE ECREVISSE EF Q KGTEST kgtest R Z 20 45 DEFAUT Q ITER GLOB MAX itnmax T 400 DEFAUT Q CRIT CONV DEBI precdb I 1 E 5 DEFAUT GENFRAL Q COURBES AUCUNE DEFAUT POSTSCRIPT Kn INTERACTIF Kn LOGICIEL logiciel Kn VERSION Vas 2 t Kn Manuel d utilisation Fascicule u7 03 Chainage avec un a
21. giciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Fo Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 8 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 3 5 3 6 3 7 3 8 3 9 3 10 3 11 Liste des pas de temps qui correspondent aux instants de calcul impos s En cas de red coupage les instants sont ins r s dans cette liste Op rande CHAM MATER CHAM MATER chmat Nom du champ de mat riau affect au maillage Attention ce champ doit comprendre les donn es associ es au comportement m canique et au comportement thermique mot cl THER de DEFI MATERIAU Il faut galement avoir d fini la temp rature de r f rence sous AFFE VARC Op rande CARA ELEM Q CARA ELEM carac Ce mot cl permet de renseigner les caract ristiques des l ments de poutre barre coque tuyau l ment discret lorsqu ils sont pr sents dans le mod le Attention ce mot cl n est utilis que pour le calcul m canique mais n est pas transmis pour l op rateur de thermique pour viter tout probl me avec les barres Mot cl EXCIT MECA EXCIT MECA F CHARGE chi Q FONC MULT Safi Q TYPE CHARGE s FIXE CSTE f SUIV JSMDIDI II s agit ici de renseigner les conditions aux limites du probl me m canique ainsi que les ventuels chargements
22. luide diphasique gt ou lt Air vapeur satur e op rande FLUIDE ENTREE 2 ou 5 3 15 6 Op randes PRES PART PRES PART FO PRES PART pae PRES PART FO fpae Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Fo Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 12 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Cet op rande permet de renseigner la pression partielle d air Pa ou son volution au cours du temps en amont de la fissure Il doit tre renseign si le fluide l entr e est lt Air vapeur surchauff e gt ou lt Air vapeur satur e gt op rande FLUIDE ENTREE 4 ou 5 3 16 Mot cl MODELE _ECRE MODELE ECRE Ce mot cl permet de renseigner les mod les qui seront utilis s par Ecrevisse pour calculer l coulement du fluide 3 16 1 Op rande IVENAC IVENAC ou 1 Lorsque cet op rande est activ IVENAC 1 le calcul ECREVISSE s effectue avec prise en compte de la vena contracta perte de charge en entr e Sinon cas par d faut le calcul s effectuera sans cette option de mod lisation 3 16 2 Op rande ECOULEMENT ECOULEMENT SATURATION GELE En l absence d air et en pr sence d eau op rande FLUIDE ENTREE 1 ou 2 cet op
23. m caniques qui s appliquent la structure La syntaxe est identique celle du mot cl EXCIT de STAT NON LINE mais il n est pas possible de piloter un chargement Voir le document U4 51 03 Mot cl CONTACT CONTACT char contact On renseigne ici les conditions de contact entre les l vres de la fissure qui peuvent se refermer pendant le calcul Le but de ce mot cl est celui d viter l interp n tration des l vres qui peut se produire sous l effet de l coulement du fluide Mot cl EXCIT _THER EXCIT THER F CHARGE chi Q FONC MULT fi Il s agit ici de renseigner les conditions aux limites du probl me thermique ainsi que les ventuels chargements thermiques La syntaxe est identique celle du mot cl EXCIT de THER LINEAIRE Voir le document U4 54 01 Mot cl COMP_INCR Il s agit ici de renseigner la loi de comportement qui sera utilis pour r soudre le probl me m canique Voir le document U4 51 03 Mot cl NEWTON Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Ma Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 9 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Il s agit ici de renseigner les param tres de l algorithme de Newton pour r soudre le probl me m canique Ce mot
24. r Ecrevisse 7 2 2 Op randes CHARG THER1 et CHARG THER2 Concepts sortants contenant les chargements thermiques fournis par Ecrevisse 7 2 3 Mot cl RESULTAT 7 2 3 1 Op rande MECANIQUE Contient le r sultat de STAT NON LINE 7 2 3 2 Op rande THERMIQUE Contient le r sultat de THER LINEAIRE 7 2 3 3 Op rande NUME ORDRE Num ro d ordre du pas de temps de calcul en cours 7 2 3 4 Op rande INST Valeur de l instant en secondes du pas de temps de calcul en cours 7 3 Syntaxe de la proc dure MACR ECRE CALC MACR ECRE CALC CONCEPT SORTANT TABLE CO table CO DEBIT CO debit CO DONNEES GEOMETRIQUES RELATIVES A LA FISSURE FISSURE E LONGUEUR long R RUGOSITE rug R ANGLE alpha R ZETA zeta R SECTION ELLIPSE RECTANGLE Q LISTE COTES BL 0 max abs curv DEFAUT 7 Lepr listr8 LISTE VAL BL lvbl listr8 DONNEES RELATIVES A L ECOULE ECOULEMENT idem ACR ECREVISSE ECOULEMENT DONNEES RELATIVES AU PROFIL DE TEMPERATURE A TRAVERS LA PAROI Q TEMPERATURE SIERE N GRADIENT FOURNI LISTE COTES TEMP lct R L
25. rande permet de choisir entre le mod le homog ne l quilibre et le mod le d coulement GELE caract ris par une fraction non nulle de liquide m tastable H 3 16 3 Op rande PRESS_EBULLITION Pression d bullition Pa fournir uniquement pour le mod le d coulement gel 3 16 4 Op rande FROTTEMENT Il s agit ici de d finir comment est calcul le frottement Les valeurs valides sont 14 12 11 4 3 2 1 0 1 2 3 4 11 12 14 La valeur 0 correspond un calcul sans frottement Pour les valeurs n gatives l utilisateur fixe la valeur du coefficient de frottement pour les coulements turbulents Pour les valeurs positives le frottement en coulement turbulent est calcul partir du coefficient de rugosit De plus pour les valeurs sup rieurs 10 la loi utilis e pour l coulement laminaire est raccord e avec celle en coulement turbulent Autrement une discontinuit est pr sente Selon la valeur utilis le calcul de la viscosit dynamique se fait diff remment important pour les coulements diphasiques cf documentation Ecrevisse 3 16 5 Op randes REYNOLDS_LIM FROTTEMENT LIM Coefficient de Reynolds limite et coefficient de frottement impos pour un Reynolds sup rieur au Reynolds limite A fournir seulement si FROTTEMENT lt 0 3 16 6 Op rande TRANSFERT_CHAL Cet op rande permet de d finir si on
26. rande LISTE COTES AH Selon la valeur de l op rande SECTION ELLIPSE Liste de cotes des points d finissant le grand axe de la section RECTANGLE Liste de cotes des points d finissant la longueur de la section 7 4 1 4 Op rande LISTE VAL AH Selon la valeur de l op rande SECTION ELLIPSE Liste des valeurs des points d finissant le grand axe de la section RECTANGLE Liste des valeurs des points d finissant la longueur de la section 7 4 2 Mot cl TEMPERATURE 7 4 2 1 Op rande GRADIENT Mod le du gradient de temp rature Fourni Distribution de temp rature fournie Impos Distribution impos e de temp rature Calcul Profil de temp rature calcul 7 4 2 2 Op rande LI STE COTES TEMP Liste des cotes pour les temp ratures 7 4 2 3 Op rande LISTE VAL TEMP Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster os Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 20 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 Liste des valeurs de temp rature 7 4 2 4 Op rande TEMP1 Gradient de temp rature de la paroi le long de l coulement 7 4 2 5 Op rande TEMP2 Temp rature de la paroi l entr e 7 4 2 6 Op rande EPAISSEUR PAROI paisseur de la paroi m 7 4 2 7 Op rande CONVEC
27. ression flux vitesse coefficient de convection 3 24 Op rande INFO INFO 1 presque aucun affichage INFO 2 affiche les informations concernant le crit re de convergence 4 Fonctionnement du couplage 41 Principe Manuel d utilisation Fascicule u7 03 Cha nage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Ps Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 15 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 4 2 4 3 Cette macro commande a t d velopp e dans un cadre pr cis estimer les d bits d une fissure dans un massif en b ton mod lis en 2D contenant ventuellement des aciers mod lis s sous forme de barre non pris en compte dans le calcul thermique L utilisation dans tout autre cadre doit donc se faire avec la plus grande pr caution Toutefois il y a en fait 3 macro commande imbriqu es ce qui peut permettre de n utiliser qu une partie des fonctionnalit s L algorithme est le suivant D but de la macro Boucle sur la liste d instants MACR ECREVISSE py R cup ration chargement thermique et m canique Ex cution de THER LINEAIRE Projection du chargement thermique et d finition des chargements venant d Ecrevisse Ex cutionde STAT NON LINE Entr e dans CALC ECREVISSE py R cup ration de l ouverture de la fissure POST RELEVE T R cup ration de la temp
28. ssss 18 Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Proc dure MACR_ECREVISSE Responsable Marina BOTTONI But Date 22 04 2013 Page 3 20 Cl U7 03 41 R vision 10960 Cette commande permet le cha nage de Code Aster avec le logiciel ECREVISSE Le principe repose pour chaque pas de temps un appel successif TH ER LIN FAIR Remarque ll est n cessaire de pr ciser le mot cl DEBUG F HIST ETAPE OUI STAT NON LIN POURSUITE pour utiliser MACR ECREVISSE E pour obtenir l tat thermique de la structure E pour obtenir l tat m canique de la structure au code ECREVISSE pour conna tre les conditions d coulement du fluide dans DEBUT ou 2 Syntaxe MACR ECREVISSE CONCEPTS SORTANTS TABLE CO table CO Q DEBIT CO debit CO TEMPER CO temp CO ETAT INITIAL 0 ETAT INIT F EVOL NOLI evol _ noli evol_noli EVOL THER evol ther evol_ther NUME ORDRE nume ordre I MODELES THERMO MECANIQUES MODELE MEC
29. t CRIT_CONV_DEBI Par d faut 10 des valeurs jusqu 10 am liore le temps de calcul sans perdre beaucoup de pr cision 3 19 Op rande COURBES Q COURBES AUCUNE DEFAUT POSTSCRIPT INTERACTIF Permet de g n rer ou non les courbes de sortie d Ecrevisse Si les courbes sont g n r es elles sont stock es dans le dossier REPE OUT au format postscript tous les instants sont dans le m me fichier postscript En mode interactif les courbes sont affich es pour chaque pas de temps 3 20 Op rande LOGICIEL Permet de sp cifier en dur le chemin pour trouver l ex cutable d Ecrevisse A pr sent seules les versions 3 2 et 3 2 1 sont compatibles avec la MACR_ECREVISSE II faut renseigner obligatoirement cet op rande ou alternativement l op rande VERSION 3 21 Op rande VERSION Indique la version d Ecrevisse utilis e dans le couplage Il faut renseigner obligatoirement cet op rande ou alternativement l op rande LOGICIEL A pr sent la version 3 2 1 est utilis e Elle est la seule qu on peut renseigner sous l op rande VERSION Si on souhaite utiliser la version 3 2 celle ci doit tre renseign e sous l op rande LOGICIEL 3 22 Op rande ENTETE Sert donner un titre au calcul Ecrevisse 3 23 Op rande IMPRESSION Lorsque l op rande vaut oOUI la macro imprime une table r capitulative des r sultats d Ecrevisse le long de l abscisse curviligne temp rature p
30. thermiques et hydrauliques sont rapides il est donc n cessaire l aussi de raffiner suffisamment le maillage La discr tisation en temps est galement tr s importante et doit tre suffisamment fine sous peine d avoir des carts importants dans les r ponses Mieux vaut essayer d optimiser cette liste d instants plut t que de se reposer uniquement sur les crit res pour red couper car cette strat gie est beaucoup plus co teuse en temps La macro commande peut tre r entrante Les r sultats thermiques et m caniques seront bien enrichis par contre de nouvelles tables de sortie seront cr s pour la nouvelle liste de pas de temps Remarque les chargements sont d truits chaque pas de temps ainsi que le champ de mat riaux les post traitements possibles sont de fait limit s d placement contrainte et variables internes aux points de Gauss Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster a Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 16 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 5 Exemple d utilisation Voir test zzzz218a b et c 6 R f rence Ecrevisse 1 Note H 186 2004 00741 FR Simulation d coulement air eau vapeur travers une fissure logiciel ECREVISSE version 3 0 Note de principe HERVOUET Chantal 2 H 181
31. utre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Proc dure MACR ECREVISSE Responsable Marina BOTTONI ENTETE Q IMPRESSION 9 INFO il 2 titre NON OUI Date 22 04 2013 Page 6 20 Cl U7 03 41 R vision 10960 DEFAUT r DEFAUT r Manuel d utilisation Fascicule u7 03 Chainage avec un autre logiciel Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Aa Titre Proc dure MACR_ECREVISSE Date 22 04 2013 Page 7 20 Responsable Marina BOTTONI Cl U7 03 41 R vision 10960 3 Op randes 3 1 Op randes TABLE TEMPER DEBIT Q TEMPER CO temp TABLE CO table DEBIT CO debit Ces mots cl s permettent de donner le nom des concepts sortant de la macro commande savoir TABLE une table de 7 colonnes donnant pour chaque pas de temps la cote z le flux la pression totale Pa la temp rature du fluide C le coefficient d change par convection W m C7 le nom de la fissure TEMPER la structure de donn es r sultat du calcul thermique DEBIT la valeur du d bit kgs chaque pas de temps et pour chaque fissure Remarques e fous les concepts de
Download Pdf Manuals
Related Search