ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/mgoperation/src/main.cpp
Revision: 1050
Committed: Thu Sep 17 22:50:46 2020 UTC (4 years, 11 months ago) by francois
File size: 94527 byte(s)
Log Message:
correction  de l'aide SiDolo + verification de l'existance des fichiers dans le calcul magic 

File Contents

# Content
1
2 //------------------------------------------------------------
3 //------------------------------------------------------------
4 // MAGiC
5 // Jean Christophe Cuilli�e et Vincent FRANCOIS
6 // D�artement de G�ie M�anique - UQTR
7 //------------------------------------------------------------
8 // Le projet MAGIC est un projet de recherche du d�artement
9 // de g�ie m�anique de l'Universit�du Qu�ec �
10 // Trois Rivi�es
11 // Les librairies ne peuvent �re utilis�s sans l'accord
12 // des auteurs (contact : francois@uqtr.ca)
13 //------------------------------------------------------------
14 //------------------------------------------------------------
15 //
16 // main.cpp
17 //
18 //------------------------------------------------------------
19 //------------------------------------------------------------
20 // COPYRIGHT 2000
21 // Version du 02/03/2006 �11H25
22 //------------------------------------------------------------
23 //------------------------------------------------------------
24 #include "gestionversion.h"
25 #ifdef WINDOWS_VERSION
26 #include "fenetre.h"
27 #endif
28
29
30 #include <string.h>
31 #include "mg_file.h"
32 #include "mg_export.h"
33 #include "fct_export.h"
34 #include "mg_import.h"
35 #include "step_import.h"
36 #include "acis_import.h"
37 #include "step_import.h"
38 #include "stl_import.h"
39 #include "occ_import.h"
40 #include "mailleur_fem.h"
41 #include "mc_gestionnaire.h"
42 #include "vct_multi_modele.h"
43 #include "mgaster.h"
44 #include <ctype.h>
45 #include <math.h>
46 #include "mailleur_analyse.h"
47 #include "magic_application.h"
48 #include "fct_taille_fem_solution.h"
49 #include "fem_maillage_quadratique_outils.h"
50 #include "fem_maillage_outils.h"
51 #include "mg_table.h"
52 #include "occ_fonction_v2017.h"
53 #include "mg_cg_modele.h"
54 #include "mg_cg_assemblage.h"
55 //---------------------------------------------------------------------------
56
57
58 void minuscule(char * mess)
59 {
60 int nb=strlen(mess);
61 for (int i=0;i<nb;i++)
62 mess[i]=tolower(mess[i]);
63 }
64
65
66 int main(int argc,char **argv)
67 {
68 std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
69 MAGIC_PARAMETRE_APPLICATION p1(2,(char*)"-gmsh",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en gmsh");
70 MAGIC_PARAMETRE_APPLICATION p2(3,(char*)"-nodeforme",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ne convertit pas les deformes des maillages FEM");
71 MAGIC_PARAMETRE_APPLICATION p3(4,(char*)"-copie",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Copie un fichier MAGiC en un autre");
72 MAGIC_PARAMETRE_APPLICATION p4(5,(char*)"-out",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de sortie");
73 MAGIC_PARAMETRE_APPLICATION p5(6,(char*)"-aster",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion MAGiC en format ASTER (uniquement calcul mecanique)");
74 MAGIC_PARAMETRE_APPLICATION p6(7,(char*)"-coderesu",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Choix des valeurs à calculer par EF - Defaut 11111111 - |Surface moyenne|Surface inferieure|Surface superieur|Deplacement|Contraintes|Deformations|Contraintes equivalentes|Energie|");
75 MAGIC_PARAMETRE_APPLICATION p7(8,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion carte de taille ancien format en format gmsh");
76 MAGIC_PARAMETRE_APPLICATION p8(9,(char*)"-step",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier step en format MAGiC");
77 MAGIC_PARAMETRE_APPLICATION p9(10,(char*)"-unite",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur d'unité pour obtenir des métres - Defaut 1.");
78 MAGIC_PARAMETRE_APPLICATION p10(11,(char*)"-stepocc",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier step en format MAGiC via opencascade");
79 MAGIC_PARAMETRE_APPLICATION p11(12,(char*)"-sat",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier sat en format MAGiC");
80 MAGIC_PARAMETRE_APPLICATION p12(13,(char*)"-eps",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Tolérance pour un double - Defaut 0.000001");
81 MAGIC_PARAMETRE_APPLICATION p13(14,(char*)"-epsgeo",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Distance max entre la triangulation et la géométrie - Defaut 1.");
82 MAGIC_PARAMETRE_APPLICATION p14(15,(char*)"-importstl",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Importation de la triangulation STL");
83 MAGIC_PARAMETRE_APPLICATION p15(16,(char*)"-importtriangulation",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Importation de la triangulation STL avec noeuds fusionnés");
84 MAGIC_PARAMETRE_APPLICATION p16(17,(char*)"-occ",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier opencascade en format MAGiC");
85 MAGIC_PARAMETRE_APPLICATION p17(18,(char*)"-fem",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Création d'un maillage fem");
86 MAGIC_PARAMETRE_APPLICATION p18(19,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de maillage - Defaut le premier");
87 MAGIC_PARAMETRE_APPLICATION p19(20,(char*)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Degre du maillage à créer");
88 MAGIC_PARAMETRE_APPLICATION p20(21,(char*)"-astermagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation des resultats ASTER dans un maillage");
89 MAGIC_PARAMETRE_APPLICATION p21(22,(char*)"-fichieraster",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier aster");
90 MAGIC_PARAMETRE_APPLICATION p22(23,(char*)"-optis",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier magic en format OPTIS - Fonction esactivée");
91 MAGIC_PARAMETRE_APPLICATION p23(24,(char*)"-optismagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation des résultats OPTIS");
92 MAGIC_PARAMETRE_APPLICATION p24(25,(char*)"-fichieroptis",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier optis");
93 MAGIC_PARAMETRE_APPLICATION p25(26,(char*)"-fichierflottant",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier flottant - Optionnel");
94 MAGIC_PARAMETRE_APPLICATION p26(27,(char*)"-seuil",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Seuil de densité pour conserver les éléments");
95 MAGIC_PARAMETRE_APPLICATION p27(28,(char*)"-niveau",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur pour le seuil (parametre incertain)");
96 MAGIC_PARAMETRE_APPLICATION p28(29,(char*)"-selectgr",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Selection de groupe topologique");
97 MAGIC_PARAMETRE_APPLICATION p29(30,(char*)"-deselectgr",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Deselection de groupe topologique");
98 MAGIC_PARAMETRE_APPLICATION p30(31,(char*)"-numgeo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de geometrie - Defaut la première");
99 MAGIC_PARAMETRE_APPLICATION p31(32,(char*)"-numgt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de groupe topologique - Defaut le premier");
100 MAGIC_PARAMETRE_APPLICATION p32(33,(char*)"-numentite",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro d'entité à selectionner");
101 MAGIC_PARAMETRE_APPLICATION p33(34,(char*)"-expand",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Prend en compte la topologie sous-jacente dans la selection");
102 MAGIC_PARAMETRE_APPLICATION p34(35,(char*)"-examine",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Contenu d'un fichier MAGiC. Si un numero de maillage est présent l'action analyse ce maillage");
103 MAGIC_PARAMETRE_APPLICATION p35(36,(char*)"-angle",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Angle limite pour l'analyse de retournement de triangle - Defaut 0.03");
104 MAGIC_PARAMETRE_APPLICATION p36(37,(char*)"-fichierlog",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier log");
105 MAGIC_PARAMETRE_APPLICATION p37(38,(char*)"-deforme",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Calcul un maillage déforme");
106 MAGIC_PARAMETRE_APPLICATION p38(39,(char*)"-fusiongeo",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Fusionne deux géométries");
107 MAGIC_PARAMETRE_APPLICATION p39(40,(char*)"-numsol1",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 1");
108 MAGIC_PARAMETRE_APPLICATION p40(41,(char*)"-numsol2",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 2");
109 MAGIC_PARAMETRE_APPLICATION p41(42,(char*)"-numsol3",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 3");
110 MAGIC_PARAMETRE_APPLICATION p42(43,(char*)"-numchamp1",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 1");
111 MAGIC_PARAMETRE_APPLICATION p43(44,(char*)"-numchamp2",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 2");
112 MAGIC_PARAMETRE_APPLICATION p44(45,(char*)"-numchamp3",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 3");
113 MAGIC_PARAMETRE_APPLICATION p45(46,(char*)"-ccf",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie");
114 MAGIC_PARAMETRE_APPLICATION p46(47,(char*)"-ccft",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous jacente");
115 MAGIC_PARAMETRE_APPLICATION p47(48,(char*)"-ccfhelp",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Liste des conditions aux limites disponibles");
116 MAGIC_PARAMETRE_APPLICATION p48(49,(char*)"-valeur",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Condition aux limites à appliquer sous la forme id:ccf:valeur(:NS|SU) ou id:ccf:formule:var1,var2(:NS|SU)");
117 MAGIC_PARAMETRE_APPLICATION p49(50,(char*)"-oricoque",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Orientation d'une coque");
118 MAGIC_PARAMETRE_APPLICATION p50(51,(char*)"-numcoque",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de la coque - Defaut la premiere");
119 MAGIC_PARAMETRE_APPLICATION p51(52,(char*)"-stlmagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation de fichier STL");
120 MAGIC_PARAMETRE_APPLICATION p52(53,(char*)"-epsf",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Tolérance pour un float - Defaut 0.0001");
121 MAGIC_PARAMETRE_APPLICATION p53(54,(char*)"-calculaster",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Calcul ASTER");
122 MAGIC_PARAMETRE_APPLICATION p54(55,(char*)"-correspondance",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'un fichier de correspondance entre GMSH et MAGIC");
123 MAGIC_PARAMETRE_APPLICATION p55(56,(char*)"-carte2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion carte de taille nouvelle version en format gmsh");
124 MAGIC_PARAMETRE_APPLICATION p56(57,(char*)"-stlout",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en stl");
125 MAGIC_PARAMETRE_APPLICATION p57(58,(char*)"-comparetaille",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compare la taille d'un maillage avec la taille d'une carte de taille isotrope ancienne version");
126 MAGIC_PARAMETRE_APPLICATION p58(59,(char*)"-cartetaille",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom de la carte de taille a fournir");
127 MAGIC_PARAMETRE_APPLICATION p59(60,(char*)"-comparetaille2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compare la taille d'un maillage avec la taille d'une carte de taille isotrope nouvelle version");
128 MAGIC_PARAMETRE_APPLICATION p60(61,(char*)"-coefdef",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Coefficient de déformé - par defaut 0. = pas de creation de maillage");
129 MAGIC_PARAMETRE_APPLICATION p61(62,(char*)"-detachemaillage",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Détache un maillage des ses liens avec la géométrie");
130 MAGIC_PARAMETRE_APPLICATION p62(63,(char*)"-elastique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en elasticité");
131 MAGIC_PARAMETRE_APPLICATION p63(64,(char*)"-thermique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en thermique");
132 MAGIC_PARAMETRE_APPLICATION p64(65,(char*)"-contrainteplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en contrainte plane");
133 MAGIC_PARAMETRE_APPLICATION p65(66,(char*)"-deformationplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en deformation plane");
134 MAGIC_PARAMETRE_APPLICATION p66(67,(char*)"-sonde",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Obtenir un resultat sur une arete dans un fichier csv");
135 MAGIC_PARAMETRE_APPLICATION p67(68,(char*)"-numsol",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de la solution (la premiere par defaut)");
136 MAGIC_PARAMETRE_APPLICATION p68(69,(char*)"-numchamp",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero du champs de la solution (le premier par defaut)");
137 MAGIC_PARAMETRE_APPLICATION p69(70,(char*)"-csvout",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier csv");
138 MAGIC_PARAMETRE_APPLICATION p70(71,(char*)"-numare",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de l'arete à sonder (la premiere par default)");
139 MAGIC_PARAMETRE_APPLICATION p71(72,(char*)"-projsolution",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Projeter une solution sur un autre maillage");
140 MAGIC_PARAMETRE_APPLICATION p72(73,(char*)"-fichiersolution",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier contenant la solution");
141 MAGIC_PARAMETRE_APPLICATION p73(74,(char*)"-comparesolution",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Comparer deux champs de solution. La reference est la deuxième solution.");
142 MAGIC_PARAMETRE_APPLICATION p74(75,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Analyse le maillage quadratique. Defaut non.");
143 MAGIC_PARAMETRE_APPLICATION p75(76,(char*)"-table",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Extraire des tables de valeurs d'un maillage fem");
144 MAGIC_PARAMETRE_APPLICATION p76(77,(char*)"-noeud",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les noeuds");
145 MAGIC_PARAMETRE_APPLICATION p77(78,(char*)"-element1",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element1D");
146 MAGIC_PARAMETRE_APPLICATION p78(79,(char*)"-element2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element2D");
147 MAGIC_PARAMETRE_APPLICATION p79(80,(char*)"-element3",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element3D");
148 MAGIC_PARAMETRE_APPLICATION p80(81,(char*)"-numtopo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Action uniquement sur le numero de topologie mentionnée");
149 MAGIC_PARAMETRE_APPLICATION p81(82,(char*)"-maigmsh",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation d'un maillage de GMSH (-in peut être omis -out peut etre omis si -in n'est pas omis -numgeo le defaut correspond à un non attechement du maillage)");
150 MAGIC_PARAMETRE_APPLICATION p82(83,(char*)"-fichiergmsh",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier gmsh");
151 MAGIC_PARAMETRE_APPLICATION p83(84,(char*)"-param",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier parametre");
152 MAGIC_PARAMETRE_APPLICATION p84(85,(char*)"-creeparam",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Création d'un fichier parametre");
153 MAGIC_PARAMETRE_APPLICATION p85(86,(char*)"-importefonction",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation d'une fonction dans une géométrie");
154 MAGIC_PARAMETRE_APPLICATION p86(87,(char*)"-fichier",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier à importer");
155 MAGIC_PARAMETRE_APPLICATION p87(88,(char*)"-elastique_gd",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en elasticite grandes deformations");
156 MAGIC_PARAMETRE_APPLICATION p88(89,(char*)"-plastique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en plasticite");
157 MAGIC_PARAMETRE_APPLICATION p89(90,(char*)"-base",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Base de numerotation des entités code aster - Defaut 10");
158 MAGIC_PARAMETRE_APPLICATION p90(91,(char*)"-element_gauss1",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element1D");
159 MAGIC_PARAMETRE_APPLICATION p91(92,(char*)"-element_gauss2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element2D");
160 MAGIC_PARAMETRE_APPLICATION p92(93,(char*)"-element_gauss3",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element3D");
161 MAGIC_PARAMETRE_APPLICATION p93(94,(char*)"-ccfi",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie par identificateur d'une autre entité MAGiC");
162 MAGIC_PARAMETRE_APPLICATION p94(95,(char*)"-ccfit",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous jacente par identificateur d'une autre entité MAGiC");
163 MAGIC_PARAMETRE_APPLICATION p95(96,(char*)"-ccff",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie sous forme de formule");
164 MAGIC_PARAMETRE_APPLICATION p96(97,(char*)"-ccfft",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous forme de formule");
165 MAGIC_PARAMETRE_APPLICATION p97(98,(char*)"-opt",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Optimisation de la numerotation des noeuds");
166 MAGIC_PARAMETRE_APPLICATION p98(99,(char*)"-rigidite",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite elastique par code aster");
167 MAGIC_PARAMETRE_APPLICATION p99(100,(char*)"-rigiditecontrainteplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite en contrainte plane par code aster");
168 MAGIC_PARAMETRE_APPLICATION p100(101,(char*)"-rigiditedeformationplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite en deformation plane par code aster");
169 MAGIC_PARAMETRE_APPLICATION p101(102,(char*)"-elementnd",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les tous les elements (solution de calcul mixte)");
170 MAGIC_PARAMETRE_APPLICATION p102(103,(char*)"-changefonction",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Modifie une fonction dans une géométrie");
171 MAGIC_PARAMETRE_APPLICATION p103(104,(char*)"-numfonction",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero d'une fonction dans la géométrie (defaut la premiere)");
172 MAGIC_PARAMETRE_APPLICATION p104(105,(char*)"-offout",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en off");
173 MAGIC_PARAMETRE_APPLICATION p105(106,(char*)"-mince",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion d'une coquille mince");
174 MAGIC_PARAMETRE_APPLICATION p106(107,(char*)"-volcor",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Ajout d'une face de base et d'une face en vis a vis pour un volume mince sous la forme idgeo:idvol:idfacebase:idfacevisavis");
175 MAGIC_PARAMETRE_APPLICATION p107(108,(char*)"-elastique_nl",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en elasticité non lineaire");
176 MAGIC_PARAMETRE_APPLICATION p108(109,(char*)"-fusion_entite_similaire",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Fusionner les entites similaires par recherche iterative de correspondance");
177 MAGIC_PARAMETRE_APPLICATION p109(110,(char*)"-brep_occ",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier BREP OpenCASCADE en format MAGiC V2017");
178 MAGIC_PARAMETRE_APPLICATION p110(111,(char*)"-step_occ",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier STEP en format MAGiC via OpenCASCADE V2017");
179 MAGIC_PARAMETRE_APPLICATION p111(112,(char*)"-compressmagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compression et rassemblement d'un fichier magic");
180 MAGIC_PARAMETRE_APPLICATION p112(113,(char*)"-compresscarte",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compression et rassemblement d'un fichier carte ancienne version");
181 MAGIC_PARAMETRE_APPLICATION p113(114,(char*)"-compresscarte2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compression et rassemblement d'un fichier carte nouvelle version");
182 MAGIC_PARAMETRE_APPLICATION p114(115,(char*)"-tristl",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation de la triangulation STL via OpenCASCADE V2017");
183 MAGIC_PARAMETRE_APPLICATION p115(116,(char*)"-brep_out",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Exportation du fichier .brep via OpenCASCADE V2017");
184 MAGIC_PARAMETRE_APPLICATION p116(117,(char*)"-abaqus",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en abaqus");
185
186 p1.ajouter_dependance(1);p1.ajouter_dependance(3);p1.ajouter_dependance(55);
187 p3.ajouter_dependance(1);p3.ajouter_dependance(5);
188 p5.ajouter_dependance(1);p5.ajouter_dependance(7);p5.ajouter_dependance(84);p5.ajouter_dependance(85);
189 p7.ajouter_dependance(1);
190 p8.ajouter_dependance(1);p8.ajouter_dependance(10);
191 p10.ajouter_dependance(1);p10.ajouter_dependance(13);p10.ajouter_dependance(14);p10.ajouter_dependance(15);;p10.ajouter_dependance(16);
192 p11.ajouter_dependance(1);p11.ajouter_dependance(10);
193 p16.ajouter_dependance(1);p16.ajouter_dependance(13);p16.ajouter_dependance(14);p16.ajouter_dependance(15);;p16.ajouter_dependance(16);
194 p17.ajouter_dependance(1);p17.ajouter_dependance(5);p17.ajouter_dependance(19);p17.ajouter_dependance(20);p17.ajouter_dependance(75);p17.ajouter_dependance(98);
195 p20.ajouter_dependance(1);p20.ajouter_dependance(19);p20.ajouter_dependance(22);p20.ajouter_dependance(90);
196 //p22.ajouter_dependance(1);p22.ajouter_dependance(7);
197 p23.ajouter_dependance(1);p23.ajouter_dependance(19);p23.ajouter_dependance(25);p23.ajouter_dependance(26);p23.ajouter_dependance(27);p23.ajouter_dependance(28);
198 p28.ajouter_dependance(1);p28.ajouter_dependance(31);p28.ajouter_dependance(32);p28.ajouter_dependance(33);p28.ajouter_dependance(34);
199 p29.ajouter_dependance(1);p29.ajouter_dependance(31);p29.ajouter_dependance(32);p29.ajouter_dependance(33);p29.ajouter_dependance(34);
200 p34.ajouter_dependance(1);p34.ajouter_dependance(19);p34.ajouter_dependance(36);p34.ajouter_dependance(37);
201 p37.ajouter_dependance(1);p37.ajouter_dependance(19);p37.ajouter_dependance(40);p37.ajouter_dependance(41);p37.ajouter_dependance(42);p37.ajouter_dependance(43);p37.ajouter_dependance(44);p37.ajouter_dependance(45);p37.ajouter_dependance(61);
202 p38.ajouter_dependance(1);p38.ajouter_dependance(31);
203 p45.ajouter_dependance(1);p45.ajouter_dependance(31);p45.ajouter_dependance(49);
204 p46.ajouter_dependance(1);p46.ajouter_dependance(31);p46.ajouter_dependance(49);
205 p49.ajouter_dependance(1);p49.ajouter_dependance(19);p49.ajouter_dependance(51);p49.ajouter_dependance(5);
206 p51.ajouter_dependance(1);p51.ajouter_dependance(10);p51.ajouter_dependance(53);
207 p53.ajouter_dependance(1);p53.ajouter_dependance(5);p53.ajouter_dependance(7);p53.ajouter_dependance(19);p53.ajouter_dependance(63);p53.ajouter_dependance(64);p53.ajouter_dependance(65);p53.ajouter_dependance(66);p53.ajouter_dependance(84);p53.ajouter_dependance(85);p53.ajouter_dependance(88);p53.ajouter_dependance(89);p53.ajouter_dependance(99);p53.ajouter_dependance(100);p53.ajouter_dependance(101);p53.ajouter_dependance(108);
208 p55.ajouter_dependance(1);
209 p56.ajouter_dependance(1);p56.ajouter_dependance(19);
210 p57.ajouter_dependance(1);p57.ajouter_dependance(19);p57.ajouter_dependance(59);
211 p59.ajouter_dependance(1);p59.ajouter_dependance(19);p59.ajouter_dependance(59);
212 p61.ajouter_dependance(1);p61.ajouter_dependance(5);p61.ajouter_dependance(19);
213 p66.ajouter_dependance(1);p66.ajouter_dependance(31);p66.ajouter_dependance(68);p66.ajouter_dependance(69);p66.ajouter_dependance(70);p66.ajouter_dependance(71);
214 p71.ajouter_dependance(1);p71.ajouter_dependance(73);p71.ajouter_dependance(68);p71.ajouter_dependance(69);p71.ajouter_dependance(19);p71.ajouter_dependance(5);
215 p73.ajouter_dependance(1);p73.ajouter_dependance(5);p73.ajouter_dependance(40);p73.ajouter_dependance(41);p73.ajouter_dependance(43);p73.ajouter_dependance(44);
216 p75.ajouter_dependance(1);p75.ajouter_dependance(5);p75.ajouter_dependance(19);p75.ajouter_dependance(77);p75.ajouter_dependance(78);p75.ajouter_dependance(79);p75.ajouter_dependance(80);p75.ajouter_dependance(81);p75.ajouter_dependance(91);p75.ajouter_dependance(92);p75.ajouter_dependance(93);p75.ajouter_dependance(102);
217 p81.ajouter_dependance(1);p81.ajouter_dependance(5);p81.ajouter_dependance(31);p81.ajouter_dependance(83);
218 p85.ajouter_dependance(1);p85.ajouter_dependance(31);p85.ajouter_dependance(87);
219 p93.ajouter_dependance(1);p93.ajouter_dependance(31);p93.ajouter_dependance(49);
220 p94.ajouter_dependance(1);p94.ajouter_dependance(31);p94.ajouter_dependance(49);
221 p95.ajouter_dependance(1);p95.ajouter_dependance(31);p95.ajouter_dependance(49);
222 p96.ajouter_dependance(1);p96.ajouter_dependance(31);p96.ajouter_dependance(49);
223 p102.ajouter_dependance(1);p102.ajouter_dependance(31);p102.ajouter_dependance(87);p102.ajouter_dependance(104);
224 p104.ajouter_dependance(1);p104.ajouter_dependance(19);
225 p105.ajouter_dependance(1);p105.ajouter_dependance(107);
226 p109.ajouter_dependance(1);p109.ajouter_dependance(5);p109.ajouter_dependance(13);p109.ajouter_dependance(14);p109.ajouter_dependance(15);p109.ajouter_dependance(16);p109.ajouter_dependance(109);
227 p110.ajouter_dependance(1);p110.ajouter_dependance(5);p110.ajouter_dependance(13);p110.ajouter_dependance(14);p110.ajouter_dependance(15);p110.ajouter_dependance(16);p110.ajouter_dependance(109);
228 p111.ajouter_dependance(1);
229 p112.ajouter_dependance(1);
230 p113.ajouter_dependance(1);
231 p114.ajouter_dependance(1);p114.ajouter_dependance(5);p114.ajouter_dependance(14);p114.ajouter_dependance(31);
232 p115.ajouter_dependance(1);p115.ajouter_dependance(5);p115.ajouter_dependance(31);
233 p116.ajouter_dependance(1);p116.ajouter_dependance(3);
234 lst.push_back(p1);
235 lst.push_back(p2);
236 lst.push_back(p3);
237 lst.push_back(p4);
238 lst.push_back(p5);
239 lst.push_back(p6);
240 lst.push_back(p7);
241 lst.push_back(p8);
242 lst.push_back(p9);
243 lst.push_back(p10);
244 lst.push_back(p11);
245 lst.push_back(p12);
246 lst.push_back(p13);
247 lst.push_back(p14);
248 lst.push_back(p15);
249 lst.push_back(p16);
250 lst.push_back(p17);
251 lst.push_back(p18);
252 lst.push_back(p19);
253 lst.push_back(p20);
254 lst.push_back(p21);
255 lst.push_back(p22);
256 lst.push_back(p23);
257 lst.push_back(p24);
258 lst.push_back(p25);
259 lst.push_back(p26);
260 lst.push_back(p27);
261 lst.push_back(p28);
262 lst.push_back(p29);
263 lst.push_back(p30);
264 lst.push_back(p31);
265 lst.push_back(p32);
266 lst.push_back(p33);
267 lst.push_back(p34);
268 lst.push_back(p35);
269 lst.push_back(p36);
270 lst.push_back(p37);
271 lst.push_back(p38);
272 lst.push_back(p39);
273 lst.push_back(p40);
274 lst.push_back(p41);
275 lst.push_back(p42);
276 lst.push_back(p43);
277 lst.push_back(p44);
278 lst.push_back(p45);
279 lst.push_back(p46);
280 lst.push_back(p47);
281 lst.push_back(p48);
282 lst.push_back(p49);
283 lst.push_back(p50);
284 lst.push_back(p51);
285 lst.push_back(p52);
286 lst.push_back(p53);
287 lst.push_back(p54);
288 lst.push_back(p55);
289 lst.push_back(p56);
290 lst.push_back(p57);
291 lst.push_back(p58);
292 lst.push_back(p59);
293 lst.push_back(p60);
294 lst.push_back(p61);
295 lst.push_back(p62);
296 lst.push_back(p63);
297 lst.push_back(p64);
298 lst.push_back(p65);
299 lst.push_back(p66);
300 lst.push_back(p67);
301 lst.push_back(p68);
302 lst.push_back(p69);
303 lst.push_back(p70);
304 lst.push_back(p71);
305 lst.push_back(p72);
306 lst.push_back(p73);
307 lst.push_back(p74);
308 lst.push_back(p75);
309 lst.push_back(p76);
310 lst.push_back(p77);
311 lst.push_back(p78);
312 lst.push_back(p79);
313 lst.push_back(p80);
314 lst.push_back(p81);
315 lst.push_back(p82);
316 lst.push_back(p83);
317 lst.push_back(p84);
318 lst.push_back(p85);
319 lst.push_back(p86);
320 lst.push_back(p87);
321 lst.push_back(p88);
322 lst.push_back(p89);
323 lst.push_back(p90);
324 lst.push_back(p91);
325 lst.push_back(p92);
326 lst.push_back(p93);
327 lst.push_back(p94);
328 lst.push_back(p95);
329 lst.push_back(p96);
330 lst.push_back(p97);
331 lst.push_back(p98);
332 lst.push_back(p99);
333 lst.push_back(p100);
334 lst.push_back(p101);
335 lst.push_back(p102);
336 lst.push_back(p103);
337 lst.push_back(p104);
338 lst.push_back(p105);
339 lst.push_back(p106);
340 lst.push_back(p107);
341 lst.push_back(p108);
342 lst.push_back(p109);
343 lst.push_back(p110);
344 lst.push_back(p111);
345 lst.push_back(p112);
346 lst.push_back(p113);
347 lst.push_back(p114);
348 lst.push_back(p115);
349 lst.push_back(p116);
350 MAGIC_APPLICATION app((char*)"Operation dans l'environnement MAGiC",argc,argv,lst,true,false);
351 char nomfichier[3000];
352 char* p;
353 if (app.get_action()!=48)
354 {
355 if (app.get_erreur()==true) return 0;
356 if ((app.get_action()!=82) && (app.get_action()!=54)&& (app.get_action()!=6))
357 {
358 app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
359 p=strrchr(nomfichier,'.');
360 }
361 }
362 char chaine[1000];for (int i=0;i<1000;i++) chaine[i]=0;
363 if (app.get_action()==2)
364 {
365 bool desactivedeforme;
366 app.recupere_parametre_bool_avec_defaut(4,desactivedeforme,false);
367 bool corres;
368 app.recupere_parametre_bool_avec_defaut(55,corres,false);
369 if (app.get_erreur()==true) return 0;
370 char mess[2000];
371 sprintf(mess,"Conversion MAGiC vers GMSH");
372 app.affiche(mess);
373 MG_FILE gest(nomfichier);
374 int nb=gest.get_nb_mg_maillage();
375 for (int i=0;i<nb;i++)
376 {
377 MG_MAILLAGE* mai=gest.get_mg_maillage(i);
378 unsigned long id=mai->get_id();
379 strncpy(chaine,nomfichier,p-nomfichier);
380 char chaine2[3000];
381 sprintf(chaine2,"%s%lu",chaine,id);
382 std::string namefic=chaine2;
383 MG_EXPORT exp;
384 sprintf(mess," Maillage geometrique : creation de %s.msh",namefic.c_str());
385 app.affiche(mess);
386 if (corres==false) exp.gmsh(mai,namefic); else exp.gmsh_avec_correspondance(mai,namefic);
387 }
388 nb=gest.get_nb_fem_maillage();
389 for (int i=0;i<nb;i++)
390 {
391 FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
392 if (desactivedeforme) mai->desactive_deforme();
393 unsigned long id=mai->get_id();
394 strncpy(chaine,nomfichier,p-nomfichier);
395 char chaine2[3000];
396 sprintf(chaine2,"%s%lu",chaine,id);
397 std::string namefic=chaine2;
398 MG_EXPORT exp;
399 sprintf(mess," Maillage EF : Creation de %s.msh",namefic.c_str());
400 app.affiche(mess);
401 if (corres==false) exp.gmsh(mai,namefic); else exp.gmsh_avec_correspondance(mai,namefic);
402 }
403 }
404 if (app.get_action()==4)
405 {
406 char fichierout[2000];
407 app.recupere_parametre_string(5,fichierout,(char*)"Fichier sortant manquant");
408 if (app.get_erreur()==true) return 0;
409 char mess[2000];
410 sprintf(mess,"Copie de fichier MAGiC");
411 app.affiche(mess);
412 MG_FILE gest(nomfichier);
413 gest.enregistrer(fichierout);
414 }
415 if (app.get_action()==6)
416 {
417 char mess[500];
418 char coderesu[10];
419 char nomparam[500];
420 bool creeparam;
421 app.recupere_parametre_string_avec_defaut(7,coderesu,(char*)"11111111");
422 app.recupere_parametre_bool_avec_defaut(85,creeparam,false);
423 if (!creeparam) app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant");
424 if (!creeparam) app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
425 //app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant");
426 if (app.get_erreur()==true) return 0;
427 if (creeparam)
428 {
429 char mess[500];
430 sprintf(mess," Calcul ASTER\n\nCréation d'un fichier parametre\n");
431 app.affiche(mess);
432 MG_EXPORT exp;
433 exp.ecrire_params_aster(nomparam);
434 }
435 else
436 {
437 sprintf(mess,"Conversion MAGiC vers code ASTER");
438 app.affiche(mess);
439 MG_FILE gest(nomfichier);
440 int nb=gest.get_nb_fem_maillage();
441 for (int i=0;i<nb;i++)
442 {
443 FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
444 MG_GEOMETRIE* geo=mai->get_mg_geometrie();
445 unsigned long id=mai->get_id();
446 p=strrchr(nomfichier,'.');
447 strncpy(chaine,nomfichier,p-nomfichier);
448 char chaine2[3000];
449 sprintf(chaine2,"%s%lu",chaine,id);
450 std::string namefic=chaine2;
451 MG_EXPORT exp;
452 exp.lire_params_aster(nomparam);
453 exp.aster(geo,mai,namefic,MAGIC::CALCUL_ASTER::ELASTIQUE,coderesu);
454 }
455
456 }
457 }
458 if (app.get_action()==8)
459 {
460 if (app.get_erreur()==true) return 0;
461 char mess[500];
462 sprintf(mess," Conversion carte de taille ancienne version vers GMSH");
463 app.affiche(mess);
464 FCT_GENERATEUR_3D<4> carte;
465 carte.lire(nomfichier);
466 strncpy(chaine,nomfichier,p-nomfichier);
467 std::string namefic=chaine;
468 FCT_EXPORT exp;
469 exp.gmsh(&carte,namefic);
470 }
471 #ifdef BREP_STEP
472 if (app.get_action()==9)
473 {
474 double unite;
475 app.recupere_parametre_double_avec_defaut(10,unite,1.);
476 if (app.get_erreur()==true) return 0;
477 char mess[500];
478 sprintf(mess," Conversion STEP vers MAGiC");
479 app.affiche(mess);
480 STEP_IMPORT stepimport;
481 MG_GESTIONNAIRE gest;
482 stepimport.importer(gest,nomfichier) ;
483 strncpy(chaine,nomfichier,p-nomfichier);
484 std::string namefic=chaine;
485 namefic=namefic + ".magic";
486 int nb=gest.get_nb_mg_geometrie();
487 for (int i=0;i<nb;i++)
488 gest.get_mg_geometrie(i)->change_valeur_unite(unite);
489 gest.enregistrer(namefic.c_str());
490 }
491 #endif
492 #ifdef BREP_OCC
493 if (app.get_action()==11)
494 {
495 double eps,eps2;
496 bool importstl,importtriangulation;
497 app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
498 app.recupere_parametre_double_avec_defaut(14,eps2,1.);
499 app.recupere_parametre_bool_avec_defaut(15,importstl,false);
500 app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
501 if (app.get_erreur()==true) return 0;
502 int importtri=0;
503 if (importstl) importtri=1;
504 if (importtriangulation) importtri=2;
505 char mess[500];
506 sprintf(mess," Conversion STEP via OpenCascade vers MAGiC avec epsilon=%lf",eps);
507 app.affiche(mess);
508 MG_GESTIONNAIRE gest;
509 OCC_IMPORT occimport;
510 MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIERSTEP,eps) ;
511 if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
512 strncpy(chaine,nomfichier,p-nomfichier);
513 std::string namefic=chaine;
514 namefic=namefic + ".magic";
515 gest.enregistrer(namefic.c_str());
516 }
517 #endif
518 #ifdef BREP_SAT
519 if (app.get_action()==12)
520 {
521 double unite;
522 app.recupere_parametre_double_avec_defaut(10,unite,1.);
523 if (app.get_erreur()==true) return 0;
524 char mess[500];
525 sprintf(mess," Conversion SAT vers MAGiC");
526 app.affiche(mess);
527 MG_GESTIONNAIRE gest;
528 ACIS_IMPORT acisimport;
529 acisimport.importer(gest,nomfichier) ;
530 strncpy(chaine,nomfichier,p-nomfichier);
531 std::string namefic=chaine;
532 namefic=namefic + ".magic";
533 int nb=gest.get_nb_mg_geometrie();
534 for (int i=0;i<nb;i++)
535 gest.get_mg_geometrie(i)->change_valeur_unite(unite);
536 gest.enregistrer(namefic.c_str());
537 }
538 #endif
539 #ifdef BREP_OCC
540 if (app.get_action()==17)
541 {
542 double eps,eps2;
543 bool importstl,importtriangulation;
544 app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
545 app.recupere_parametre_double_avec_defaut(14,eps2,1.);
546 app.recupere_parametre_bool_avec_defaut(15,importstl,false);
547 app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
548 if (app.get_erreur()==true) return 0;
549 int importtri=0;
550 if (importstl) importtri=1;
551 if (importtriangulation) importtri=2;
552 char mess[500];
553 sprintf(mess," Conversion OpenCascade vers MAGiC");
554 app.affiche(mess);
555 MG_GESTIONNAIRE gest;
556 OCC_IMPORT occimport;
557 MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIEROCC,eps) ;
558 if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
559 strncpy(chaine,nomfichier,p-nomfichier);
560 std::string namefic=chaine;
561 namefic=namefic + ".magic";
562 gest.enregistrer(namefic.c_str());
563 }
564 #endif
565 if (app.get_action()==18)
566 {
567 int degre,nummai;
568 bool analyse;
569 bool optno;
570 char fichierout[2000];
571 app.recupere_parametre_string_avec_defaut(5,fichierout,nomfichier);
572 app.recupere_parametre_int_avec_defaut(19,nummai,0);
573 app.recupere_parametre_int(20,degre,(char*)"Degré du maillage manquant");
574 app.recupere_parametre_bool_avec_defaut(75,analyse,false);
575 app.recupere_parametre_bool_avec_defaut(98,optno,false);
576 if (app.get_erreur()==true) return 0;
577 char mess[500];
578 sprintf(mess," Creation d'un maillage FEM");
579 app.affiche(mess);
580 MG_FILE gest(nomfichier);
581 MG_MAILLAGE* mai;
582 if (nummai==0) mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
583 FEM_MAILLAGE* fem=new FEM_MAILLAGE(mai->get_mg_geometrie(),mai,degre);
584 gest.ajouter_fem_maillage(fem);
585 MAILLEUR_FEM m;
586 m.active_affichage(app.affiche);
587 if (!optno) m.maille(fem,0); else m.maille(fem,1);
588 gest.enregistrer(fichierout);
589 if (analyse==true)
590 {
591 MAILLEUR_ANALYSE m3d(fem);
592 m3d.active_affichage(app.affiche);
593 m3d.analyse();
594 }
595 /*if (degre==2)
596 {
597 //char message [5000];
598 FEM_MAILLAGE_QUADRATIQUE_OUTILS ou ;
599 double moy=ou.get_distortion_moy(fem);
600 double min=ou.get_distortion_min (fem);
601 int nb=ou.get_nombre_elements_invalides(fem);
602 sprintf(mess," maillage quadratique :");
603 app.affiche(mess);
604 sprintf(mess," ");
605 app.affiche(mess);
606 sprintf(mess," distorsion moy %f ",moy);
607 app.affiche(mess);
608 sprintf(mess," distorsion min %f ",min);
609 app.affiche(mess);
610 sprintf(mess," nombre elements invalides %d ",nb);
611 app.affiche(mess);
612 }*/
613
614
615 }
616 if (app.get_action()==21)
617 {
618 int nummai;
619 char fichieraster[500];
620 int base;
621 app.recupere_parametre_int_avec_defaut(19,nummai,0);
622 app.recupere_parametre_string(22,fichieraster,(char*)"Fichier aster manquant");
623 app.recupere_parametre_int_avec_defaut(90,base,10);
624 if (app.get_erreur()==true) return 0;
625 char mess[500];
626 sprintf(mess," Importation des résultats ASTER vers MAGIC");
627 app.affiche(mess);
628 MG_FILE gest(nomfichier);
629 FEM_MAILLAGE* mai;
630 if (nummai==0) mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
631 MG_IMPORT imp;
632 std::string res=imp.aster(base,mai,fichieraster,nomfichier);
633 app.affiche((char*)("\033[1;31m"+res+"\033[1;32m").c_str());
634 gest.enregistrer(nomfichier);
635 }
636
637 if (app.get_action()==23)
638 {
639
640 }
641 if (app.get_action()==24)
642 {
643 int nummai;
644 char fichieroptis[500];
645 char fichierflottant[500];
646 double niveau,seuil;
647 app.recupere_parametre_int_avec_defaut(19,nummai,0);
648 app.recupere_parametre_string(25,fichieroptis,(char*)"Fichier OPTIS manquant");
649 app.recupere_parametre_string_avec_defaut(26,fichieroptis,(char*)"");
650 app.recupere_parametre_double(27,seuil,(char*)"Seuil manquant");
651 app.recupere_parametre_double(28,niveau,(char*)"Niveau manquant");
652 if (app.get_erreur()==true) return 0;
653 char mess[500];
654 sprintf(mess," Importation des résultats d'optimisation vers MAGIC");
655 app.affiche(mess);
656 MG_FILE gest(nomfichier);
657 FEM_MAILLAGE* mai;
658 if (nummai==0) mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
659 MG_IMPORT imp;
660 int flottant=0;
661 if (strlen(fichierflottant)>0) flottant=1;
662 imp.optis(mai,fichieroptis,fichierflottant,seuil,niveau,flottant);
663 gest.enregistrer(nomfichier);
664 }
665
666 if (app.get_action()==29)
667 {
668 int numgeo,numgt,numentite;
669 bool expand;
670 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
671 app.recupere_parametre_int_avec_defaut(32,numgt,0);
672 app.recupere_parametre_int(33,numentite,(char*)"Numéro d'entité manquant");
673 app.recupere_parametre_bool_avec_defaut(34,expand,false);
674 if (app.get_erreur()==true) return 0;
675 char mess[500];
676 sprintf(mess," Selection de groupe topologique");
677 app.affiche(mess);
678 MG_FILE gest(nomfichier);
679 MG_GEOMETRIE* mggeo;
680 if (numgeo==0) mggeo=gest.get_mg_geometrie(numgeo); else mggeo=gest.get_mg_geometrieid(numgeo);
681 MG_GROUPE_TOPOLOGIQUE* mggt;
682 if (numgt==0)
683 {
684 mggt=new MG_GROUPE_TOPOLOGIQUE;
685 mggeo->ajouter_mg_groupe_topologique(mggt);
686 }
687 else mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
688 MG_ELEMENT_TOPOLOGIQUE* ele;
689 ele=mggeo->get_mg_sommetid(numentite);
690 if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
691 if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
692 if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
693 if (ele!=NULL)
694 {
695 mggt->ajouter(ele);
696 if (expand)
697 {
698 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
699 ele->get_topologie_sousjacente(&lst);
700 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
701 for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
702 mggt->ajouter(ele2);
703 }
704 }
705 gest.enregistrer(nomfichier);
706 }
707 if (app.get_action()==30)
708 {
709 int numgeo,numgt,numentite;
710 bool expand;
711 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
712 app.recupere_parametre_int_avec_defaut(32,numgt,0);
713 app.recupere_parametre_int(33,numentite,(char*)"Numéro d'entité manquant");
714 app.recupere_parametre_bool_avec_defaut(34,expand,false);
715 if (app.get_erreur()==true) return 0;
716 char mess[500];
717 sprintf(mess," Deselection de groupe topologique");
718 app.affiche(mess);
719 MG_FILE gest(nomfichier);
720 MG_GEOMETRIE* mggeo=gest.get_mg_geometrieid(numgeo);
721 MG_GROUPE_TOPOLOGIQUE* mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
722 MG_ELEMENT_TOPOLOGIQUE* ele;
723 ele=mggeo->get_mg_sommetid(numentite);
724 if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
725 if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
726 if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
727 if (ele!=NULL)
728 {
729 mggt->supprimer(ele);
730 if (expand)
731 {
732 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
733 ele->get_topologie_sousjacente(&lst);
734 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
735 for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
736 mggt->supprimer(ele2);
737 }
738 }
739 gest.enregistrer(nomfichier);
740 }
741
742 if (app.get_action()==35)
743 {
744 int nummai;
745 double angle;
746 char fichierlog[500];
747 strcpy(fichierlog,"");
748 app.recupere_parametre_int_avec_defaut(19,nummai,-1);
749 app.recupere_parametre_double_avec_defaut(36,angle,0.03);
750 app.recupere_parametre_string_avec_defaut(37,fichierlog,(char*)fichierlog);
751 if (app.get_erreur()==true) return 0;
752 char mess[500];
753 sprintf(mess," Contenu d'un fichier MAGIC");
754 app.affiche(mess);
755 MG_FILE gest(nomfichier);
756 if (nummai!=-1)
757 {
758 MG_MAILLAGE* mai=gest.get_mg_maillageid(nummai);
759 if (mai!=NULL)
760 {
761 MAILLEUR_ANALYSE m3d(mai);
762 m3d.active_affichage(app.affiche);
763 m3d.change_eps_angle_retourne(angle);
764 if (strlen(fichierlog)==0)
765 m3d.analyse();
766 else
767 m3d.analyse(fichierlog);
768 }
769 else
770 {
771 FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
772 if (mai!=NULL)
773 {
774 MAILLEUR_ANALYSE m3d(mai);
775 m3d.active_affichage(app.affiche);
776 m3d.change_eps_angle_retourne(angle);
777 if (strlen(fichierlog)==0)
778 m3d.analyse();
779 else
780 m3d.analyse(fichierlog);
781 }
782
783 }
784
785 }
786
787 else
788 {
789 char message [50000];
790 int nb_mgcg_modele=gest.get_nb_mgcg_modele();
791 sprintf(message," \033[1;32m%d\033[1;33m mgcg_modele :",nb_mgcg_modele);
792 for (int i=0;i<nb_mgcg_modele;i++)
793 {
794 MG_CG_MODELE* mgcg_modele=gest.get_mgcg_modele(i);
795 sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,mgcg_modele->get_id());
796 }
797 app.affiche(message);
798 for (int i=0;i<nb_mgcg_modele;i++)
799 {
800 MG_CG_MODELE* mgcg_modele=gest.get_mgcg_modele(i);
801 sprintf(message," mgcg_modele \033[1;31m%lu\033[1;33m",mgcg_modele->get_id());
802 app.affiche(message);
803 int nb_forme=mgcg_modele->get_nb_mgcg_forme();
804 sprintf(message," \033[1;32m%d\033[1;33m mgcg_forme",nb_forme);
805 app.affiche(message);
806 int nb_operateur=mgcg_modele->get_nb_mgcg_operateur();
807 sprintf(message," \033[1;32m%d\033[1;33m mgcg_operateur",nb_operateur);
808 app.affiche(message);
809 int nb_assemblage=mgcg_modele->get_nb_mgcg_assemblage();
810 sprintf(message," \033[1;32m%d\033[1;33m mgcg_assemblage : ",nb_assemblage);
811 std::map<long,MG_CG_ASSEMBLAGE*>::iterator it_assemblage;
812 for(MG_CG_ASSEMBLAGE* mgcg_assemblage=mgcg_modele->get_premier_mgcg_assemblage(it_assemblage);mgcg_assemblage!=NULL;mgcg_assemblage=mgcg_modele->get_suivant_mgcg_assemblage(it_assemblage))
813 {
814 sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,mgcg_assemblage->get_id());
815 }
816 app.affiche(message);
817 }
818 int nbarbre=gest.get_nb_mg_arbre();
819 sprintf(message," \033[1;32m%d\033[1;33m arbres :",nbarbre);
820 for (int i=0;i<nbarbre;i++)
821 {
822 MG_ARBRE* arb=gest.get_mg_arbre(i);
823 sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,arb->get_id());
824 }
825 app.affiche(message);
826 for (int i=0;i<nbarbre;i++)
827 {
828 MG_ARBRE* arb=gest.get_mg_arbre(i);
829 sprintf(message," arbre \033[1;31m%lu\033[1;33m",arb->get_id());
830 app.affiche(message);
831 int nbass=arb->get_nb_mg_assemblage();
832 sprintf(message," \033[1;32m%d\033[1;33m assemblages",nbass);
833 app.affiche(message);
834 int nbprim=arb->get_nb_mg_primitive();
835 sprintf(message," \033[1;32m%d\033[1;33m primitives",nbprim);
836 //for (int i=0;i<nbprim;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,arb->get_mg_primitive(i)->get_id());
837 app.affiche(message);
838 int nbope=arb->get_nb_mg_operateur_boolean();
839 sprintf(message," \033[1;32m%d\033[1;33m operateurs booleans",nbope);
840 app.affiche(message);
841 }
842 int nbgeo=gest.get_nb_mg_geometrie();
843 sprintf(message," \033[1;32m%d\033[1;33m geometries :",nbgeo);
844 for (int i=0;i<nbgeo;i++)
845 {
846 MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
847 sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,geo->get_id());
848 }
849 app.affiche(message);
850 for (int i=0;i<nbgeo;i++)
851 {
852 MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
853 sprintf(message," geometrie \033[1;31m%lu\033[1;33m",geo->get_id());
854 app.affiche(message);
855 int nbmg=geo->get_nb_mg_groupe_topologique();
856 sprintf(message," \033[1;32m%d\033[1;33m groupe topologique",nbmg);
857 for (int j=0;j<nbmg;j++)
858 {
859 MG_GROUPE_TOPOLOGIQUE* mggp=geo->get_mg_groupe_topologique(j);
860 sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,mggp->get_id());
861 }
862 app.affiche(message);
863 int nbvolume=geo->get_nb_mg_volume();
864 if (nbvolume==0) sprintf(message," \033[1;32m%d\033[1;33m volumes",geo->get_nb_mg_volume());
865 else sprintf(message," \033[1;32m%d\033[1;33m volumes : ",geo->get_nb_mg_volume());
866 for (int i=0;i<nbvolume;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_volume(i)->get_id());
867 app.affiche(message);
868 int nbcoque=geo->get_nb_mg_coque();
869 if (nbcoque==0) sprintf(message," \033[1;32m%d\033[1;33m coques",geo->get_nb_mg_coque());
870 else sprintf(message," \033[1;32m%d\033[1;33m coques : ",geo->get_nb_mg_coque());
871 for (int i=0;i<nbcoque;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_coque(i)->get_id());
872 app.affiche(message);
873 int nbpoutre=geo->get_nb_mg_poutre();
874 if (nbpoutre==0) sprintf(message," \033[1;32m%d\033[1;33m poutres",geo->get_nb_mg_poutre());
875 else sprintf(message," \033[1;32m%d\033[1;33m poutres : ",geo->get_nb_mg_poutre());
876 for (int i=0;i<nbpoutre;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_poutre(i)->get_id());
877 app.affiche(message);
878 sprintf(message," \033[1;32m%d\033[1;33m coquilles",geo->get_nb_mg_coquille());
879 app.affiche(message);
880 sprintf(message," \033[1;32m%d\033[1;33m faces",geo->get_nb_mg_face());
881 app.affiche(message);
882 sprintf(message," \033[1;32m%d\033[1;33m aretes",geo->get_nb_mg_arete());
883 app.affiche(message);
884 sprintf(message," \033[1;32m%d\033[1;33m sommets",geo->get_nb_mg_sommet());
885 app.affiche(message);
886 int nbgef=geo->get_nb_mg_geom_fonction();
887 if (nbgef==0) sprintf(message," \033[1;32m%d\033[1;33m fonctions geometriques",nbgef);
888 else sprintf(message," \033[1;32m%d\033[1;33m fonctions geometriques : ",nbgef);
889 for (int i=0;i<nbgef;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_geom_fonction(i)->get_id());
890 app.affiche(message);
891 }
892 int nbmai=gest.get_nb_mg_maillage();
893 sprintf(message," \033[1;32m%d\033[1;33m maillage geometrique :",nbmai);
894 app.affiche(message);
895 for (int i=0;i<nbmai;i++)
896 {
897 MG_MAILLAGE* mai=gest.get_mg_maillage(i);
898 sprintf(message," maillage geometrique \033[1;31m%lu\033[1;33m ",mai->get_id());
899 app.affiche(message);
900 sprintf(message," \033[1;32m%d\033[1;33m noeuds",mai->get_nb_mg_noeud());
901 app.affiche(message);
902 sprintf(message," \033[1;32m%d\033[1;33m segments",mai->get_nb_mg_segment());
903 app.affiche(message);
904 sprintf(message," \033[1;32m%d\033[1;33m triangles",mai->get_nb_mg_triangle());
905 app.affiche(message);
906 sprintf(message," \033[1;32m%d\033[1;33m quadrangles",mai->get_nb_mg_quadrangle());
907 app.affiche(message);
908 sprintf(message," \033[1;32m%d\033[1;33m tetras",mai->get_nb_mg_tetra());
909 app.affiche(message);
910 sprintf(message," \033[1;32m%d\033[1;33m hexas",mai->get_nb_mg_hexa());
911 app.affiche(message);
912 sprintf(message," \033[1;32m%d\033[1;33m pentas",mai->get_nb_mg_penta());
913 app.affiche(message);
914 }
915 int nbmaif=gest.get_nb_fem_maillage();
916 sprintf(message," \033[1;32m%d\033[1;33m maillage FEM :",nbmaif);
917 app.affiche(message);
918 for (int i=0;i<nbmaif;i++)
919 {
920 FEM_MAILLAGE* maif=gest.get_fem_maillage(i);
921 sprintf(message," maillage FEM \033[1;31m%lu\033[1;33m degre \033[1;31m%d\033[1;33m ",maif->get_id(),maif->get_degre());
922 app.affiche(message);
923 sprintf(message," \033[1;35m%lu\033[1;33m maillage géométrique de référence",maif->get_mg_maillage()->get_id());
924 app.affiche(message);
925 sprintf(message," \033[1;32m%d\033[1;33m noeuds",maif->get_nb_fem_noeud());
926 app.affiche(message);
927 sprintf(message," \033[1;32m%d\033[1;33m elements 1D",maif->get_nb_fem_element1());
928 app.affiche(message);
929 sprintf(message," \033[1;32m%d\033[1;33m elements 2D",maif->get_nb_fem_element2());
930 app.affiche(message);
931 sprintf(message," \033[1;32m%d\033[1;33m elements 3D",maif->get_nb_fem_element3());
932 app.affiche(message);
933 double rx,ry,rz;
934 maif->calcul_somme_reaction(rx,ry,rz);
935 sprintf(message," Reactions : \n Rx=\033[1;31m%le\033[1;33m\n Ry=\033[1;31m%le\033[1;33m\n Rz=\033[1;31m%le\033[1;33m",rx,ry,rz);
936 app.affiche(message);
937
938 }
939 int nbsol=gest.get_nb_mg_solution();
940 sprintf(message," \033[1;32m%d\033[1;33m solution maillage geometrique :",nbsol);
941 app.affiche(message);
942 message[0]=0;
943 for (int i=0;i<nbsol;i++)
944 {
945 MG_SOLUTION* sol=gest.get_mg_solution(i);
946 char legende[500];
947 legende[0]=0;
948 message[0]=0;
949 for (int j=0;j<sol->get_nb_champ();j++)
950 sprintf(legende,"%s%s_%s ",legende,(char*)sol->get_nom().c_str(),(char*)sol->get_legende(j).c_str());
951 sprintf(message," \033[1;31m%lu(%s)\033[1;33m ",sol->get_id(),legende);
952 app.affiche(message);
953 }
954 int nbsolf=gest.get_nb_fem_solution();
955 sprintf(message," \033[1;32m%d\033[1;33m solution maillage FEM :",nbsolf);
956 app.affiche(message);
957 message[0]=0;
958 for (int i=0;i<nbsolf;i++)
959 {
960 FEM_SOLUTION* sol=gest.get_fem_solution(i);
961 char legende[500];
962 legende[0]=0;
963 message[0]=0;
964 for (int j=0;j<sol->get_nb_champ();j++)
965 sprintf(legende,"%s%s_%s ",legende,(char*)sol->get_nom().c_str(),(char*)sol->get_legende(j).c_str());
966 sprintf(message," \033[1;31m%lu(%s)\033[1;33m ",sol->get_id(),legende);
967 app.affiche(message);
968 }
969 }
970 }
971 if (app.get_action()==38)
972 {
973 int nummai;
974 int numsol1,numsol2,numsol3;
975 int numchamp1,numchamp2,numchamp3;
976 double coef;
977 app.recupere_parametre_int_avec_defaut(19,nummai,0);
978 app.recupere_parametre_int(40,numsol1,(char*)"Numéro solution 1 manquant");
979 app.recupere_parametre_int(41,numsol2,(char*)"Numéro solution 2 manquant");
980 app.recupere_parametre_int(42,numsol3,(char*)"Numéro solution 3 manquant");
981 app.recupere_parametre_int(43,numchamp1,(char*)"Numéro champs 1 manquant");
982 app.recupere_parametre_int(44,numchamp2,(char*)"Numéro champs 2 manquant");
983 app.recupere_parametre_int(45,numchamp3,(char*)"Numéro champs 3 manquant");
984 app.recupere_parametre_double_avec_defaut(61,coef,0.);
985 if (app.get_erreur()==true) return 0;
986 char mess[500];
987 sprintf(mess," Calcul de deforme");
988 app.affiche(mess);
989 MG_FILE gest(nomfichier);
990 FEM_MAILLAGE* mai;
991 if (nummai==0) mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
992 FEM_SOLUTION* sol1=gest.get_fem_solutionid(numsol1);
993 FEM_SOLUTION* sol2=gest.get_fem_solutionid(numsol2);
994 FEM_SOLUTION* sol3=gest.get_fem_solutionid(numsol3);
995 mai->calcul_deforme(sol1,numchamp1,sol2,numchamp2,sol3,numchamp3);
996 if (coef>1e-16)
997 {
998 sprintf(mess," Création du maillage deforme");
999 app.affiche(mess);
1000 MG_MAILLAGE* mgmai=new MG_MAILLAGE(mai,coef);
1001 }
1002 sprintf(mess," Enregistrement");
1003 app.affiche(mess);
1004 gest.enregistrer(nomfichier);
1005 }
1006 if (app.get_action()==39)
1007 {
1008 int numgeo;
1009 app.recupere_parametre_int_avec_defaut (31,numgeo,0);
1010 if (app.get_erreur()==true) return 0;
1011 char mess[500];
1012 sprintf(mess," Fusion de géométrie");
1013 app.affiche(mess);
1014 MG_FILE gest(nomfichier);
1015 MG_GEOMETRIE* mggeo;
1016 if (numgeo==0) mggeo=gest.get_mg_geometrie(numgeo); else mggeo=gest.get_mg_geometrieid(numgeo);
1017 VCT_MULTI_MODELE fusion(mggeo);
1018 fusion.recherche_identite();
1019 gest.enregistrer(nomfichier);
1020 }
1021 if ((app.get_action()==46) || (app.get_action()==47))
1022 {
1023 int numgeo;
1024 char ccf[500];
1025 app.recupere_parametre_int_avec_defaut (31,numgeo,0);
1026 app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
1027 if (app.get_erreur()==true) return 0;
1028 char mess[500];
1029 sprintf(mess," Selection de condition aux limites");
1030 unsigned long identity;
1031 char typeccf[3];
1032 char effortsuiv[3]={'N','S',0};
1033 typeccf[2]=0;
1034 double valccf=0.5; // Valeur par défaut
1035 sscanf(ccf,"%lu:%c%c:%le:%c%c\n",&identity,&typeccf[0],&typeccf[1],&valccf,&effortsuiv[0],&effortsuiv[1]);
1036 MC_GESTIONNAIRE mcgest;
1037 int res=mcgest.existe(typeccf);
1038 if (res==0)
1039 {
1040 sprintf(mess," ERREUR!!!! Code conditions aux limites inconnu");
1041 app.affiche(mess);
1042 return 0;
1043 }
1044 std::string effortsuivstr=effortsuiv;
1045 if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
1046 {
1047 sprintf(mess," ERREUR!!!! type de suivi des efforts inconnus");
1048 app.affiche(mess);
1049 return 0;
1050 }
1051 sprintf(mess," Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
1052 app.affiche(mess);
1053 MG_FILE gest(nomfichier);
1054 MG_GEOMETRIE* geo;
1055 if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
1056 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
1057 MG_VOLUME* ele;
1058 ele=geo->get_mg_volumeid(identity);
1059 if (ele!=NULL) liste.ajouter(ele);
1060 else
1061 {
1062 MG_COQUE* ele;
1063 ele=geo->get_mg_coqueid(identity);
1064 if (ele!=NULL) liste.ajouter(ele);
1065 else
1066 {
1067 MG_POUTRE* ele;
1068 ele=geo->get_mg_poutreid(identity);
1069 if (ele!=NULL) liste.ajouter(ele);
1070 else
1071 {
1072 MG_FACE* ele;
1073 ele=geo->get_mg_faceid(identity);
1074 if (ele!=NULL) liste.ajouter(ele);
1075 else
1076 {
1077 MG_ARETE* ele;
1078 ele=geo->get_mg_areteid(identity);
1079 if (ele!=NULL) liste.ajouter(ele);
1080 else
1081 {
1082 MG_SOMMET* ele;
1083 ele=geo->get_mg_sommetid(identity);
1084 if (ele!=NULL) liste.ajouter(ele);
1085 }
1086 }
1087 }
1088 }
1089 }
1090 if (app.get_action()==47)
1091 {
1092 int nb=liste.get_nb();
1093 for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
1094 }
1095 int nb=liste.get_nb();
1096 for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,effortsuiv);
1097 gest.enregistrer(nomfichier);
1098 }
1099 if ((app.get_action()==94) || (app.get_action()==95))
1100 {
1101 int numgeo;
1102 char ccf[500];
1103 app.recupere_parametre_int_avec_defaut (31,numgeo,0);
1104 app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
1105 if (app.get_erreur()==true) return 0;
1106 char mess[500];
1107 sprintf(mess," Selection de condition aux limites");
1108 unsigned long identity;
1109 char typeccf[3];
1110 char effortsuiv[3]={'N','S',0};
1111 typeccf[2]=0;
1112 unsigned long valccf=0; // Valeur par défaut
1113 sscanf(ccf,"%lu:%c%c:%lu:%c%c\n",&identity,&typeccf[0],&typeccf[1],&valccf,&effortsuiv[0],&effortsuiv[1]);
1114 MC_GESTIONNAIRE mcgest;
1115 int res=mcgest.existe(typeccf);
1116 if (res==0)
1117 {
1118 sprintf(mess," ERREUR!!!! Code conditions aux limites inconnu");
1119 app.affiche(mess);
1120 return 0;
1121 }
1122 std::string effortsuivstr=effortsuiv;
1123 if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
1124 {
1125 sprintf(mess," ERREUR!!!! type de suivi des efforts inconnus");
1126 app.affiche(mess);
1127 return 0;
1128 }
1129 sprintf(mess," Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
1130 app.affiche(mess);
1131 MG_FILE gest(nomfichier);
1132 MG_GEOMETRIE* geo;
1133 if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
1134 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
1135 MG_VOLUME* ele;
1136 ele=geo->get_mg_volumeid(identity);
1137 if (ele!=NULL) liste.ajouter(ele);
1138 else
1139 {
1140 MG_COQUE* ele;
1141 ele=geo->get_mg_coqueid(identity);
1142 if (ele!=NULL) liste.ajouter(ele);
1143 else
1144 {
1145 MG_POUTRE* ele;
1146 ele=geo->get_mg_poutreid(identity);
1147 if (ele!=NULL) liste.ajouter(ele);
1148 else
1149 {
1150 MG_FACE* ele;
1151 ele=geo->get_mg_faceid(identity);
1152 if (ele!=NULL) liste.ajouter(ele);
1153 else
1154 {
1155 MG_ARETE* ele;
1156 ele=geo->get_mg_areteid(identity);
1157 if (ele!=NULL) liste.ajouter(ele);
1158 else
1159 {
1160 MG_SOMMET* ele;
1161 ele=geo->get_mg_sommetid(identity);
1162 if (ele!=NULL) liste.ajouter(ele);
1163 }
1164 }
1165 }
1166 }
1167 }
1168 if (app.get_action()==95)
1169 {
1170 int nb=liste.get_nb();
1171 for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
1172 }
1173 int nb=liste.get_nb();
1174 for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,effortsuiv);
1175 gest.enregistrer(nomfichier);
1176 }
1177 if ((app.get_action()==96) || (app.get_action()==97))
1178 {
1179 int numgeo;
1180 char ccf[500];
1181 app.recupere_parametre_int_avec_defaut (31,numgeo,0);
1182 app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
1183 if (app.get_erreur()==true) return 0;
1184 char mess[500];
1185 sprintf(mess," Selection de condition aux limites");
1186 unsigned long identity;
1187 char typeccf[3];
1188 char effortsuiv[3]={'N','S',0};
1189 typeccf[2]=0;
1190 char valccf[500];
1191 char formule[500];
1192 sscanf(ccf,"%lu:%c%c:%s:%c%c\n",&identity,&typeccf[0],&typeccf[1],formule,&effortsuiv[0],&effortsuiv[1]);
1193 OT_CHAINE ot;
1194 std::vector<std::string> variable=ot.split(formule,':');
1195 strcpy(valccf,variable[0].c_str());
1196 std::vector<std::string> listvariable=ot.split(variable[1].c_str(),',');
1197 MC_GESTIONNAIRE mcgest;
1198 int res=mcgest.existe(typeccf);
1199 if (res==0)
1200 {
1201 sprintf(mess," ERREUR!!!! Code conditions aux limites inconnu");
1202 app.affiche(mess);
1203 return 0;
1204 }
1205 std::string effortsuivstr=effortsuiv;
1206 if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
1207 {
1208 sprintf(mess," ERREUR!!!! type de suivi des efforts inconnus");
1209 app.affiche(mess);
1210 return 0;
1211 }
1212 sprintf(mess," Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
1213 app.affiche(mess);
1214 MG_FILE gest(nomfichier);
1215 MG_GEOMETRIE* geo;
1216 if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
1217 TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
1218 MG_VOLUME* ele;
1219 ele=geo->get_mg_volumeid(identity);
1220 if (ele!=NULL) liste.ajouter(ele);
1221 else
1222 {
1223 MG_COQUE* ele;
1224 ele=geo->get_mg_coqueid(identity);
1225 if (ele!=NULL) liste.ajouter(ele);
1226 else
1227 {
1228 MG_POUTRE* ele;
1229 ele=geo->get_mg_poutreid(identity);
1230 if (ele!=NULL) liste.ajouter(ele);
1231 else
1232 {
1233 MG_FACE* ele;
1234 ele=geo->get_mg_faceid(identity);
1235 if (ele!=NULL) liste.ajouter(ele);
1236 else
1237 {
1238 MG_ARETE* ele;
1239 ele=geo->get_mg_areteid(identity);
1240 if (ele!=NULL) liste.ajouter(ele);
1241 else
1242 {
1243 MG_SOMMET* ele;
1244 ele=geo->get_mg_sommetid(identity);
1245 if (ele!=NULL) liste.ajouter(ele);
1246 }
1247 }
1248 }
1249 }
1250 }
1251 if (app.get_action()==97)
1252 {
1253 int nb=liste.get_nb();
1254 for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
1255 }
1256 int nb=liste.get_nb();
1257 for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,listvariable,effortsuiv);
1258 gest.enregistrer(nomfichier);
1259 }
1260 if (app.get_action()==48)
1261 {
1262 char mess[500];
1263 sprintf(mess," Liste des conditions aux limites disponibles");
1264 app.affiche(mess);
1265 MC_GESTIONNAIRE mcgest;
1266 int ok=0;
1267 do
1268 {
1269 std::string code;
1270 std::string description;
1271 int res=mcgest.get_description(code,description);
1272 if (res)
1273 {
1274 sprintf(mess," %s : %s",code.c_str(),description.c_str());
1275 app.affiche(mess);
1276 }
1277 else ok=1;
1278 }
1279 while (ok==0);
1280 }
1281 if (app.get_action()==50)
1282 {
1283 int numgeo,nummai,numcoque;
1284 char fichierout[500];
1285 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
1286 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1287 app.recupere_parametre_int_avec_defaut(51,numcoque,0);
1288 app.recupere_parametre_string_avec_defaut(5,fichierout,nomfichier);
1289 if (app.get_erreur()==true) return 0;
1290 char mess[500];
1291 sprintf(mess," Orientation d'une coque");
1292 app.affiche(mess);
1293 MG_FILE gest(nomfichier);
1294 MG_GEOMETRIE* geo;
1295 if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
1296 MG_MAILLAGE* mai;
1297 if (nummai==0) mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
1298 MG_COQUE* coque;
1299 if (numcoque==0) coque=mai->get_mg_geometrie()->get_mg_coque(numcoque); else coque=mai->get_mg_geometrie()->get_mg_coqueid(numcoque);
1300 coque->orienter(mai);
1301 gest.enregistrer(fichierout);
1302 }
1303 if (app.get_action()==52)
1304 {
1305 double epsf,unite;
1306 app.recupere_parametre_double_avec_defaut(10,unite,1.);
1307 app.recupere_parametre_double_avec_defaut(53,epsf,0.0001);
1308 if (app.get_erreur()==true) return 0;
1309 char mess[500];
1310 sprintf(mess," Conversion STL vers MAGiC");
1311 app.affiche(mess);
1312 MG_GESTIONNAIRE gest;
1313 STL_IMPORT stlimport;
1314 stlimport.change_eps(epsf);
1315 stlimport.importer(gest,nomfichier) ;
1316 strncpy(chaine,nomfichier,p-nomfichier);
1317 std::string namefic=chaine;
1318 namefic=namefic + ".magic";
1319 int nb=gest.get_nb_mg_geometrie();
1320 for (int i=0;i<nb;i++)
1321 gest.get_mg_geometrie(i)->change_valeur_unite(unite);
1322 gest.enregistrer(namefic.c_str());
1323 }
1324 if (app.get_action()==54)
1325 {
1326 int nummai;
1327 char coderesu[10];
1328 bool elastique;
1329 bool elastique_nl;
1330 bool elastique_gd;
1331 bool plastique;
1332 bool thermique;
1333 bool conplane;
1334 bool defplane;
1335 bool rigiditeconplane;
1336 bool rigiditedefplane;
1337 bool creeparam;
1338 bool rigidite;
1339 char nomparam[500];
1340 char nomfichierout[500];
1341 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1342 app.recupere_parametre_string_avec_defaut(7,coderesu,(char*)"11111111");
1343 app.recupere_parametre_bool_avec_defaut(63,elastique,false);
1344 app.recupere_parametre_bool_avec_defaut(108,elastique_nl,false);
1345 app.recupere_parametre_bool_avec_defaut(64,thermique,false);
1346 app.recupere_parametre_bool_avec_defaut(65,conplane,false);
1347 app.recupere_parametre_bool_avec_defaut(66,defplane,false);
1348 app.recupere_parametre_bool_avec_defaut(85,creeparam,false);
1349 app.recupere_parametre_bool_avec_defaut(88,elastique_gd,false);
1350 app.recupere_parametre_bool_avec_defaut(89,plastique,false);
1351 app.recupere_parametre_bool_avec_defaut(99,rigidite,false);
1352 app.recupere_parametre_bool_avec_defaut(100,rigiditeconplane,false);
1353 app.recupere_parametre_bool_avec_defaut(101,rigiditedefplane,false);
1354 if (!creeparam) app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant",2);
1355 app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant",2);
1356 app.recupere_parametre_string_avec_defaut(5,nomfichierout,nomfichier);
1357 if (app.get_erreur()==true) return 0;
1358 if (creeparam)
1359 {
1360 char mess[500];
1361 sprintf(mess," Calcul ASTER\n\nCréation d'un fichier parametre\n");
1362 app.affiche(mess);
1363 MG_EXPORT exp;
1364 exp.ecrire_params_aster(nomparam);
1365 }
1366 else
1367 {
1368 char mess[500];
1369 sprintf(mess," Calcul ASTER");
1370 app.affiche(mess);
1371 char nometude[1000];
1372 char *p=strchr(nomfichierout,'.');
1373 strncpy(nometude,nomfichierout,p-nomfichierout);
1374 nometude[p-nomfichierout]=0;
1375 MG_FILE gest(nomfichier);
1376 FEM_MAILLAGE* mai;
1377 if (nummai==0) mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
1378 MGASTER mgaster;
1379 mgaster.active_affichage(app.affiche);
1380 if (elastique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,coderesu);
1381 if (elastique_nl) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_NON_LINEAIRE_3D,coderesu);
1382 if (rigidite) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE,coderesu);
1383 if (conplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE,coderesu);
1384 if (defplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::DEFORMATION_PLANE,coderesu);
1385 if (rigiditeconplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE_CONTRAINTE_PLANE,coderesu);
1386 if (rigiditedefplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE_DEFORMATION_PLANE,coderesu);
1387 if (thermique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)"");
1388 if (elastique_gd) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF,coderesu);
1389 if (plastique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::PLASTIQUE,coderesu);
1390
1391 }
1392 }
1393 if (app.get_action()==56)
1394 {
1395 if (app.get_erreur()==true) return 0;
1396 char mess[500];
1397 sprintf(mess," Conversion carte de taille nouvelle version vers GMSH");
1398 app.affiche(mess);
1399 FCT_TAILLE_FEM_SOLUTION carte(nomfichier);
1400 strncpy(chaine,nomfichier,p-nomfichier);
1401 std::string namefic=chaine;
1402 MG_EXPORT exp;
1403 exp.gmsh(&carte,namefic);
1404 sprintf(mess," Carte de taille FEM : Creation de %s.msh",namefic.c_str());
1405 app.affiche(mess);
1406 }
1407 if ((app.get_action()==58) || (app.get_action()==60))
1408 {
1409 int nummai;
1410 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1411 char nomcarte[500];
1412 if (app.get_action()==58)
1413 app.recupere_parametre_string(59,nomcarte,(char*)"Fichier carte manquant");
1414 if (app.get_action()==60)
1415 app.recupere_parametre_string(59,nomcarte,(char*)"Fichier carte manquant");
1416 if (app.get_erreur()==true) return 0;
1417 char mess[500];
1418 sprintf(mess," Comparaison de la taille reelle d'un maillage avec la carte de taille");
1419 app.affiche(mess);
1420 sprintf(mess," Lecture fichier MAGiC");
1421 app.affiche(mess);
1422 MG_FILE gest(nomfichier);
1423 MG_MAILLAGE* mai;
1424 if (nummai==0) mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
1425 sprintf(mess," Lecture fichier carte de taille");
1426 app.affiche(mess);
1427 FCT_TAILLE *carte;
1428 if (app.get_action()==58)
1429 carte=new FCT_GENERATEUR_3D<4>(nomcarte);
1430 if (app.get_action()==60)
1431 carte=new FCT_TAILLE_FEM_SOLUTION(nomcarte);
1432 char nometude[1000];
1433 char *p=strchr(nomfichier,'.');
1434 strncpy(nometude,nomfichier,p-nomfichier);
1435 nometude[p-nomfichier]=0;
1436 char nomsolution[500];
1437 sprintf(nomsolution,"%s_comp%lu.sol",nometude,mai->get_id());
1438 sprintf(mess," Comparaison");
1439 app.affiche(mess);
1440 MAILLEUR_ANALYSE m3d(mai);
1441 sprintf(mess," Respect global de la taille demandee");
1442 app.affiche(mess);
1443 double tabtaille[3];
1444 m3d.compare_maillage_carte_isotrope(carte,nomsolution,tabtaille);
1445 sprintf(mess," Nombre de maille théorique : %lf - Nombre de maille reelle : %lf",tabtaille[0],tabtaille[1]);
1446 app.affiche(mess);
1447 sprintf(mess," Erreur : %lf%%",tabtaille[2]);
1448 app.affiche(mess);
1449 MG_SOLUTION *sol=gest.get_mg_solution(gest.get_nb_mg_solution()-1);
1450 sol->active_solution(2);
1451 LISTE_MG_NOEUD::iterator it;
1452 int tab[201];
1453 for (int i=0;i<201;i++) tab[i]=0;
1454 for (MG_NOEUD* no=mai->get_premier_noeud(it);no!=NULL;no=mai->get_suivant_noeud(it))
1455 {
1456 double val=no->get_solution();
1457 int valint=(int)val+100;
1458 if (valint<0) valint=0;
1459 if (valint>200) valint=200;
1460 tab[valint]++;
1461 }
1462 sprintf(mess," Repartition du respect de la taille demandee");
1463 app.affiche(mess);
1464 char nomanalyse[500];
1465 sprintf(nomanalyse,"%s_comp%lu.csv",nometude,mai->get_id());
1466 FILE *in=fopen(nomanalyse,"wt");
1467 for (int i=0;i<201;i++)
1468 if (i<101) fprintf(in,"%d;%d;\n",i-100,tab[i]);
1469 else fprintf(in,"%d;%d;%d;\n",i-100,tab[i],tab[i]+tab[200-i]);
1470 fclose(in);
1471 sprintf(mess," Enregistrement");
1472 app.affiche(mess);
1473 sprintf(nomsolution,"%s_comp%lu.magic",nometude,mai->get_id());
1474 gest.enregistrer(nomsolution);
1475 delete carte;
1476 }
1477 if (app.get_action()==57)
1478 {
1479 int nummai;
1480 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1481 if (app.get_erreur()==true) return 0;
1482 char mess[500];
1483 sprintf(mess," Conversion mg maillage en stl");
1484 app.affiche(mess);
1485 MG_FILE gest(nomfichier);
1486 strncpy(chaine,nomfichier,p-nomfichier);
1487 MG_MAILLAGE* mai;
1488 if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
1489 sprintf(chaine,"%s%lu",chaine,mai->get_id());
1490 std::string namefic=chaine;
1491 MG_EXPORT exp;
1492 exp.stl(mai,namefic);
1493 sprintf(mess," creation de %s.stl",namefic.c_str());
1494 app.affiche(mess);
1495
1496 }
1497 if (app.get_action()==105)
1498 {
1499 int nummai;
1500 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1501 if (app.get_erreur()==true) return 0;
1502 char mess[500];
1503 sprintf(mess," Conversion mg maillage en off");
1504 app.affiche(mess);
1505 MG_FILE gest(nomfichier);
1506 strncpy(chaine,nomfichier,p-nomfichier);
1507 MG_MAILLAGE* mai;
1508 if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
1509 sprintf(chaine,"%s%lu",chaine,mai->get_id());
1510 std::string namefic=chaine;
1511 MG_EXPORT exp;
1512 exp.off(mai,namefic);
1513 sprintf(mess," creation de %s.off",namefic.c_str());
1514 app.affiche(mess);
1515
1516 }
1517 if (app.get_action()==62)
1518 {
1519 int nummai;
1520 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1521 char outfile[500];
1522 app.recupere_parametre_string(5,outfile,(char*)"Fichier de sortie manquant");
1523 if (app.get_erreur()==true) return 0;
1524 char mess[500];
1525 sprintf(mess," Détachement d'un mg maillage de sa geometrie");
1526 app.affiche(mess);
1527 MG_FILE gest(nomfichier);
1528 MG_MAILLAGE* mai;
1529 if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
1530 MG_GESTIONNAIRE gest2;
1531 MG_MAILLAGE *mai2=mai->detacher(&gest2);
1532 sprintf(mess," Enregistrement");
1533 app.affiche(mess);
1534 gest2.enregistrer(outfile);
1535 }
1536 if (app.get_action()==67)
1537 {
1538 int numsol;
1539 int numchamp;
1540 int numare;
1541 int numgeo;
1542 char outfile[500];
1543 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
1544 app.recupere_parametre_int_avec_defaut(68,numsol,0);
1545 app.recupere_parametre_int_avec_defaut(69,numchamp,0);
1546 app.recupere_parametre_int_avec_defaut(71,numare,0);
1547 app.recupere_parametre_string(70,outfile,(char*)"Fichier de sortie manquant");
1548 if (app.get_erreur()==true) return 0;
1549 MG_FILE gest(nomfichier);
1550 FEM_SOLUTION* sol;
1551 if (numsol==0) sol=gest.get_fem_solution(0); else sol=gest.get_fem_solutionid(numsol);
1552 MG_GEOMETRIE* geo;
1553 if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrie(numgeo);
1554 MG_ARETE *are;
1555 if (numare==0) are=geo->get_mg_arete(0); else are=geo->get_mg_areteid(numare);
1556 char mess[500];
1557 sprintf(mess," Sonde de %s_%s sur l'arète %lu",(char*)sol->get_nom().c_str(),(char*)sol->get_legende(numchamp).c_str(),are->get_id());
1558 app.affiche(mess);
1559 are->recupere_resultat(sol,numchamp,outfile);
1560
1561 }
1562 if (app.get_action()==72)
1563 {
1564 int numsol;
1565 int numchamp;
1566 int nummai;
1567 char outfile[500];
1568 char infile[500];
1569 app.recupere_parametre_int_avec_defaut(68,numsol,0);
1570 app.recupere_parametre_int_avec_defaut(69,numchamp,0);
1571 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1572 app.recupere_parametre_string(73,infile,(char*)"Fichier de solution manquant");
1573 app.recupere_parametre_string_avec_defaut(5,outfile,infile);
1574 if (app.get_erreur()==true) return 0;
1575 app.affiche((char*)"Lecture du fichier MAGIC");
1576 MG_FILE gest(nomfichier);
1577 app.affiche((char*)"Lecture de la solution a projeter");
1578 MG_FILE gestsol(infile);
1579 FEM_SOLUTION* sol;
1580 if (numsol==0) sol=gestsol.get_fem_solution(0); else sol=gestsol.get_fem_solutionid(numsol);
1581 FEM_MAILLAGE* fem;
1582 if (nummai==0) fem=gest.get_fem_maillage(0); else fem=gest.get_fem_maillageid(nummai);
1583 char mess[500];
1584 sprintf(mess," Projection de la solution %s_%s du maillage %lu sur le maillage %lu",(char*)sol->get_nom().c_str(),(char*)sol->get_legende(numchamp).c_str(),sol->get_maillage()->get_id(),fem->get_id());
1585 app.affiche(mess);
1586 FEM_MAILLAGE_OUTILS otfem;
1587 int res=otfem.projete_solution_maillage(sol,numchamp,fem);
1588 if (res==MAGIC::PROJECTION::SUCCES) gest.enregistrer(outfile);
1589 if (res==MAGIC::PROJECTION::DIMENSION_MAILLAGE_INCOMPATIBLE) app.affiche((char*)" Erreur : Dimension des maillages incompatibles");
1590 if (res==MAGIC::PROJECTION::TYPE_SOLUTION_INCORRECT) app.affiche((char*)" Erreur : Solution a projeter n'est pas exprimer aux noeuds");
1591 if (res==MAGIC::PROJECTION::NON_IMPLEMENTER) app.affiche((char*)" Erreur : Cas de projection non implementer");
1592 }
1593 if (app.get_action()==74)
1594 {
1595 int numsol1;
1596 int numchamp1;
1597 int numsol2;
1598 int numchamp2;
1599 char outfile[500];
1600 app.recupere_parametre_int_avec_defaut(40,numsol1,0);
1601 app.recupere_parametre_int_avec_defaut(41,numsol2,0);
1602 app.recupere_parametre_int_avec_defaut(43,numchamp1,0);
1603 app.recupere_parametre_int_avec_defaut(44,numchamp2,0);
1604 app.recupere_parametre_string_avec_defaut(5,outfile,nomfichier);
1605 if (app.get_erreur()==true) return 0;
1606 MG_FILE gest(nomfichier);
1607 FEM_SOLUTION* sol1;
1608 if (numsol1==0) sol1=gest.get_fem_solution(0); else sol1=gest.get_fem_solutionid(numsol1);
1609 FEM_SOLUTION* sol2;
1610 if (numsol2==0) sol2=gest.get_fem_solution(0); else sol2=gest.get_fem_solutionid(numsol2);
1611 char mess[500];
1612 sprintf(mess," Comparaison de %s_%s et de %s_%s",(char*)sol1->get_nom().c_str(),(char*)sol1->get_legende(numchamp1).c_str(),(char*)sol2->get_nom().c_str(),(char*)sol2->get_legende(numchamp2).c_str());
1613 app.affiche(mess);
1614 FEM_MAILLAGE_OUTILS otfem;
1615 otfem.compare_champs_solution(sol1,numchamp1,sol2,numchamp2);
1616 gest.enregistrer(outfile);
1617 }
1618 if (app.get_action()==76)
1619 {
1620 bool bnoeud,belement1,belement2,belement3,bgelement1,bgelement2,bgelement3,belementnd;
1621 int nummai,numtopo;
1622 char out[500];
1623 app.recupere_parametre_string(5,out,(char*)"Fichier sortant manquant");
1624 app.recupere_parametre_int_avec_defaut(19,nummai,0);
1625 app.recupere_parametre_bool_avec_defaut(77,bnoeud,false);
1626 app.recupere_parametre_bool_avec_defaut(78,belement1,false);
1627 app.recupere_parametre_bool_avec_defaut(79,belement2,false);
1628 app.recupere_parametre_bool_avec_defaut(80,belement3,false);
1629 app.recupere_parametre_bool_avec_defaut(91,bgelement3,false);
1630 app.recupere_parametre_bool_avec_defaut(92,bgelement3,false);
1631 app.recupere_parametre_bool_avec_defaut(93,bgelement3,false);
1632 app.recupere_parametre_bool_avec_defaut(102,belementnd,false);
1633 app.recupere_parametre_int_avec_defaut(81,numtopo,0);
1634 if (app.get_erreur()==true) return 0;
1635 char mess[500];
1636 sprintf(mess," Extraction de tables de valeurs");
1637 app.affiche(mess);
1638 MG_FILE gest(nomfichier);
1639 FEM_MAILLAGE *fem;
1640 MG_ELEMENT_TOPOLOGIQUE *topo;
1641 if (nummai==0) fem=gest.get_fem_maillage(0); else fem=gest.get_fem_maillageid(nummai);
1642 if (numtopo==0) topo=NULL;
1643 else
1644 {
1645 MG_GEOMETRIE* geo=fem->get_mg_geometrie();
1646 topo=geo->get_mg_faceid(numtopo);
1647 if (topo==NULL) topo=geo->get_mg_volumeid(numtopo);
1648 if (topo==NULL) topo=geo->get_mg_areteid(numtopo);
1649 if (topo==NULL) topo=geo->get_mg_sommetid(numtopo);
1650 }
1651 MG_TABLE tab;
1652 if (bnoeud)
1653 {
1654 app.affiche((char*)" Création de la table aux noeuds");
1655 char nomout[500];
1656 sprintf(nomout,"%s_noeud.csv",out);
1657 tab.info_noeud(fem,nomout,topo);
1658 }
1659 if (belement1)
1660 {
1661 app.affiche((char*)" Création de la table aux elements 1D");
1662 char nomout[500];
1663 sprintf(nomout,"%s_element1.csv",out);
1664 tab.info_element1(fem,nomout,topo);
1665 }
1666 if (belement2)
1667 {
1668 app.affiche((char*)" Création de la table aux elements 2D");
1669 char nomout[500];
1670 sprintf(nomout,"%s_element2.csv",out);
1671 tab.info_element2(fem,nomout,topo);
1672 }
1673 if (belement3)
1674 {
1675 app.affiche((char*)" Création de la table aux elements 3D");
1676 char nomout[500];
1677 sprintf(nomout,"%s_element3.csv",out);
1678 tab.info_element3(fem,nomout,topo);
1679 }
1680 if (belementnd)
1681 {
1682 app.affiche((char*)" Création de la table aux elements de solutions mixtes");
1683 char nomout[500];
1684 sprintf(nomout,"%s_elementnd.csv",out);
1685 tab.info_element_nd(fem,nomout,topo);
1686 }
1687 if (bgelement1)
1688 {
1689 app.affiche((char*)" Création de la table aux point de gauss des elements 1D");
1690 char nomout[500];
1691 sprintf(nomout,"%s_elementgauss1.csv",out);
1692 tab.info_element_gauss1(fem,nomout,topo);
1693 }
1694 if (bgelement2)
1695 {
1696 app.affiche((char*)" Création de la table aux point de gauss des elements 2D");
1697 char nomout[500];
1698 sprintf(nomout,"%s_elementgauss2.csv",out);
1699 tab.info_element_gauss2(fem,nomout,topo);
1700 }
1701 if (bgelement3)
1702 {
1703 app.affiche((char*)" Création de la table aux point de gauss des elements 3D");
1704 char nomout[500];
1705 sprintf(nomout,"%s_elementgauss3.csv",out);
1706 tab.info_element_gauss3(fem,nomout,topo);
1707 }
1708 }
1709 if (app.get_action()==82)
1710 {
1711 int numgeo;
1712 char out[500];
1713 char in[500];
1714 char fgmsh[500];
1715 app.recupere_parametre_string_avec_defaut(1,in,(char*)"");
1716 if (strcmp(in,"")==0)
1717 app.recupere_parametre_string(5,out,(char*)"Fichier sortant manquant");
1718 else
1719 app.recupere_parametre_string_avec_defaut(5,out,in);
1720 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
1721 app.recupere_parametre_string(83,fgmsh,(char*)"Fichier GMSH manquant");
1722 if (app.get_erreur()==true) return 0;
1723 char mess[2000];
1724 sprintf(mess,"Importation d'un maillage GMSH");
1725 app.affiche(mess);
1726 MG_GESTIONNAIRE *gest;
1727 if (strcmp(in,"")==0) gest=new MG_GESTIONNAIRE; else gest=new MG_FILE(in);
1728 MG_GEOMETRIE *geo;
1729 if ((numgeo!=0) && (gest->get_nb_mg_geometrie()>0)) geo=gest->get_mg_geometrieid(numgeo);
1730 MG_IMPORT imp;
1731 imp.maillage_gmsh(gest,fgmsh,geo);
1732 sprintf(mess,"Enregistrement du fichier");
1733 app.affiche(mess);
1734 gest->enregistrer(out);
1735 delete gest;
1736 }
1737 if ((app.get_action()==86)||(app.get_action()==103))
1738 {
1739 int numgeo,numfonc;
1740 char chemin[1000];
1741 app.recupere_parametre_int_avec_defaut(31,numgeo,0);
1742 app.recupere_parametre_string(87,chemin,(char*)"Fichier de fonction manquant");
1743 if (app.get_action()==103) app.recupere_parametre_int_avec_defaut(104,numfonc,0);
1744 if (app.get_erreur()==true) return 0;
1745 char mess[2000];
1746 if (app.get_action()==86) sprintf(mess,"Importation d'une fonction geometrique");
1747 if (app.get_action()==103) sprintf(mess,"Modification d'une fonction geometrique");
1748 app.affiche(mess);
1749 sprintf(mess,"\nLecture du fichier MAGiC");
1750 app.affiche(mess);
1751 MG_FILE gest(nomfichier);
1752 MG_GEOMETRIE *geo;
1753 if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
1754 if (app.get_action()==86)
1755 {
1756 sprintf(mess,"Importation de la courbe");
1757 app.affiche(mess);
1758 MG_GEOM_FONCTION *gf=new MG_GEOM_FONCTION(chemin);
1759 geo->ajouter_mg_geom_fonction(gf);
1760 gest.enregistrer(nomfichier);
1761 }
1762 if (app.get_action()==103)
1763 {
1764 sprintf(mess,"Modification de la courbe");
1765 app.affiche(mess);
1766 MG_GEOM_FONCTION *gf;
1767 if (numfonc==0) gf=geo->get_mg_geom_fonction(0); else gf=geo->get_mg_geom_fonctionid(numfonc);
1768 gf->change(chemin);
1769 gest.enregistrer(nomfichier);
1770 }
1771 }
1772 if (app.get_action()==106)
1773 {
1774 char mess[255];
1775 sprintf(mess,"Définition d'une coque mince");
1776 app.affiche(mess);
1777 char valeur[500];
1778 app.recupere_parametre_string(107,valeur,(char*)"Faces en vis a vis manquantes");
1779 if (app.get_erreur()==true) return 0;
1780 MG_FILE gest(nomfichier);
1781 unsigned long id,id1,id2,idgeo;
1782 sscanf(valeur,"%lu:%lu:%lu:%lu\n",&idgeo,&id,&id1,&id2);
1783 MG_GEOMETRIE* geo=gest.get_mg_geometrieid(idgeo);
1784 if (geo!=NULL)
1785 {
1786 MG_VOLUME* vol=geo->get_mg_volumeid(id);
1787 if (vol!=NULL)
1788 {
1789 MG_FACE* face1=geo->get_mg_faceid(id1);
1790 MG_FACE* face2=geo->get_mg_faceid(id2);
1791 if (face1!=NULL)
1792 if (face2!=NULL)
1793 {
1794 vol->ajouter_face_correspondante(face1,face2);
1795 sprintf(mess," Ajout de la face de base %lu et de la face correspondante %lu dans le volume mince %lu",id1,id2,id);
1796 app.affiche(mess);
1797 }
1798 }
1799 }
1800 sprintf(mess,"Enregistrement du fichier");
1801 gest.enregistrer(nomfichier);
1802 }
1803
1804 #ifdef ALL_OCC
1805 if (app.get_action()==110)
1806 {
1807 double eps,eps2;
1808 bool importstl,importtriangulation,fusion_entite_similaire;
1809 char fichier_sortie[1000];
1810 char fichier_entree[1000];
1811 app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
1812 app.recupere_parametre_double_avec_defaut(14,eps2,1.);
1813 app.recupere_parametre_bool_avec_defaut(15,importstl,false);
1814 app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
1815 app.recupere_parametre_bool_avec_defaut(108,fusion_entite_similaire,false);
1816 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier brep manquant");
1817 app.recupere_parametre_string(5,fichier_sortie,(char*)"Nom du fichier MAGiC de sortie manquant");
1818 std::string str_fichier_sortie(fichier_sortie);
1819 std::string str_fichier_sortie_caf;
1820 std::string str_nom_fichier_caf_mggeo;
1821 if(str_fichier_sortie.size()>6)
1822 {
1823 std::size_t find_magic = str_fichier_sortie.find(".magic");
1824 if(find_magic!=std::string::npos)
1825 {
1826 str_fichier_sortie_caf = str_fichier_sortie;
1827 str_fichier_sortie_caf.erase(find_magic,6);
1828 str_nom_fichier_caf_mggeo = str_fichier_sortie_caf;
1829 str_fichier_sortie_caf = str_fichier_sortie_caf + ".ocaf";
1830 }
1831 }
1832 if (app.get_erreur()==true) return 0;
1833 bool importer_triangulation=false;
1834 int importtri=0;
1835 if (importstl) importtri=1;
1836 if (importtriangulation) importtri=2;
1837 if(importtri>0) importer_triangulation=true;
1838 char mess[500];
1839 sprintf(mess,"Conversion BREP OpenCascade vers MAGiC V2017");
1840 app.affiche(mess);
1841 MG_GESTIONNAIRE gest;
1842 OCC_IMPORT occ_import;
1843 MG_GEOMETRIE* mggeo = occ_import.importer_fichier_brep_V2017(gest,
1844 fichier_entree,
1845 (char*)str_fichier_sortie_caf.c_str(),
1846 FICHIEROCC,
1847 eps,
1848 fusion_entite_similaire);
1849 OCC_FONCTION_V2017* occ_fonction_v2017 = (OCC_FONCTION_V2017*)&mggeo->get_occ_fonction();
1850 sprintf(mess,"________________________________________________\n");
1851 app.affiche(mess);
1852 sprintf(mess,"-> NB MG_VOLUME : %i",mggeo->get_nb_mg_volume());
1853 app.affiche(mess);
1854 sprintf(mess,"-> NB MG_COQUILLE : %i",mggeo->get_nb_mg_coquille());
1855 app.affiche(mess);
1856 sprintf(mess,"-> NB MG_FACE : %i",mggeo->get_nb_mg_face());
1857 app.affiche(mess);
1858 sprintf(mess,"-> NB MG_COFACE : %i",mggeo->get_nb_coface());
1859 app.affiche(mess);
1860 sprintf(mess,"-> NB MG_BOUCLE : %i",mggeo->get_nb_mg_boucle());
1861 app.affiche(mess);
1862 sprintf(mess,"-> NB MG_ARETE : %i",mggeo->get_nb_mg_arete());
1863 app.affiche(mess);
1864 sprintf(mess,"-> NB MG_COARETE : %i",mggeo->get_nb_coarete());
1865 app.affiche(mess);
1866 sprintf(mess,"-> NB MG_SOMMET : %i",mggeo->get_nb_mg_sommet());
1867 app.affiche(mess);
1868 sprintf(mess,"-> NB MG_COSOMMET : %i\n",mggeo->get_nb_cosommet());
1869 app.affiche(mess);
1870 if(importer_triangulation)
1871 {
1872 occ_import.importer_triangulation_V2017(gest,mggeo,eps2,importtri);
1873 sprintf(mess,"Importation de la triangulation : OK");
1874 app.affiche(mess);
1875 }
1876 gest.enregistrer((char*)fichier_sortie);
1877 sprintf(mess,"\nFichier OpenCascade sortie : %s",str_fichier_sortie_caf.c_str());
1878 app.affiche(mess);
1879 sprintf(mess,"Fichier MAGiC sortie : %s",str_fichier_sortie.c_str());
1880 app.affiche(mess);
1881 sprintf(mess,"________________________________________________\n");
1882 app.affiche(mess);
1883 }
1884 if (app.get_action()==111)
1885 {
1886 double eps,eps2;
1887 bool importstl,importtriangulation,fusion_entite_similaire;
1888 char fichier_sortie[1000];
1889 char fichier_entree[1000];
1890 app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
1891 app.recupere_parametre_double_avec_defaut(14,eps2,1.);
1892 app.recupere_parametre_bool_avec_defaut(15,importstl,false);
1893 app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
1894 app.recupere_parametre_bool_avec_defaut(108,fusion_entite_similaire,false);
1895 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier STEP manquant");
1896 app.recupere_parametre_string(5,fichier_sortie,(char*)"Nom du fichier MAGiC de sortie manquant");
1897 std::string str_fichier_sortie(fichier_sortie);
1898 std::string str_fichier_sortie_caf;
1899 std::string str_nom_fichier_caf_mggeo;
1900 if(str_fichier_sortie.size()>6)
1901 {
1902 std::size_t find_magic = str_fichier_sortie.find(".magic");
1903 if(find_magic!=std::string::npos)
1904 {
1905 str_fichier_sortie_caf = str_fichier_sortie;
1906 str_fichier_sortie_caf.erase(find_magic,6);
1907 str_nom_fichier_caf_mggeo = str_fichier_sortie_caf;
1908 str_fichier_sortie_caf = str_fichier_sortie_caf + ".ocaf";
1909 }
1910 }
1911 if (app.get_erreur()==true) return 0;
1912 bool importer_triangulation=false;
1913 int importtri=0;
1914 if (importstl) importtri=1;
1915 if (importtriangulation) importtri=2;
1916 if(importtri>0) importer_triangulation=true;
1917 char mess[500];
1918 sprintf(mess,"Conversion STEP vers MAGiC via OpenCASCADE V2017");
1919 app.affiche(mess);
1920 MG_GESTIONNAIRE gest;
1921 OCC_IMPORT occ_import;
1922 MG_GEOMETRIE* mggeo = occ_import.importer_fichier_STEP_V2017(gest,
1923 fichier_entree,
1924 (char*)str_fichier_sortie_caf.c_str(),
1925 FICHIERSTEP,
1926 eps,
1927 fusion_entite_similaire);
1928 OCC_FONCTION_V2017* occ_fonction_v2017 = (OCC_FONCTION_V2017*)&mggeo->get_occ_fonction();
1929 sprintf(mess,"________________________________________________\n");
1930 app.affiche(mess);
1931 sprintf(mess,"-> NB MG_VOLUME : %i",mggeo->get_nb_mg_volume());
1932 app.affiche(mess);
1933 sprintf(mess,"-> NB MG_COQUILLE : %i",mggeo->get_nb_mg_coquille());
1934 app.affiche(mess);
1935 sprintf(mess,"-> NB MG_FACE : %i",mggeo->get_nb_mg_face());
1936 app.affiche(mess);
1937 sprintf(mess,"-> NB MG_COFACE : %i",mggeo->get_nb_coface());
1938 app.affiche(mess);
1939 sprintf(mess,"-> NB MG_BOUCLE : %i",mggeo->get_nb_mg_boucle());
1940 app.affiche(mess);
1941 sprintf(mess,"-> NB MG_ARETE : %i",mggeo->get_nb_mg_arete());
1942 app.affiche(mess);
1943 sprintf(mess,"-> NB MG_COARETE : %i",mggeo->get_nb_coarete());
1944 app.affiche(mess);
1945 sprintf(mess,"-> NB MG_SOMMET : %i",mggeo->get_nb_mg_sommet());
1946 app.affiche(mess);
1947 sprintf(mess,"-> NB MG_COSOMMET : %i\n",mggeo->get_nb_cosommet());
1948 app.affiche(mess);
1949 if(importer_triangulation)
1950 {
1951 occ_import.importer_triangulation_V2017(gest,mggeo,eps2,importtri);
1952 sprintf(mess,"Importation de la triangulation : OK");
1953 app.affiche(mess);
1954 }
1955 gest.enregistrer((char*)fichier_sortie);
1956 sprintf(mess,"\nFichier OpenCascade sortie : %s",str_fichier_sortie_caf.c_str());
1957 app.affiche(mess);
1958 sprintf(mess,"Fichier MAGiC sortie : %s",str_fichier_sortie.c_str());
1959 app.affiche(mess);
1960 sprintf(mess,"________________________________________________\n");
1961 app.affiche(mess);
1962 }
1963 if (app.get_action()==112)
1964 {
1965 char fichier_entree[1000];
1966 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier MAGiC manquant");
1967 if (app.get_erreur()==true) return 0;
1968 char mess[1000];
1969 sprintf(mess,"Compression et rassemblement d'un fichier MAGiC");
1970 app.affiche(mess);
1971 MG_FILE gest(fichier_entree);
1972 std::vector<std::string> liste;
1973 gest.get_fichier_dependant(liste);
1974 char message[40000];
1975 sprintf(message,"tar -zcvf %s.tar.gz %s",fichier_entree,fichier_entree);
1976 for (int i=0;i<liste.size();i++)
1977 sprintf(message,"%s %s",message,liste[i].c_str());
1978 system(message);
1979 }
1980 if ((app.get_action()==113) || (app.get_action()==114))
1981 {
1982 char fichier_entree[1000];
1983 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier carte manquant");
1984 if (app.get_erreur()==true) return 0;
1985 char mess[1000];
1986 sprintf(mess,"Compression et rassemblement d'un fichier carte de taille");
1987 app.affiche(mess);
1988
1989 FCT_TAILLE* taille;
1990 if (app.get_action()==113) taille=new FCT_GENERATEUR_3D<4>(fichier_entree);
1991 if (app.get_action()==114) taille=new FCT_TAILLE_FEM_SOLUTION(fichier_entree);
1992 std::vector<std::string> liste;
1993 taille->get_fichier_dependant(liste);
1994 char message[40000];
1995 sprintf(message,"tar -zcvf %s.tar.gz %s",fichier_entree,fichier_entree);
1996 for (int i=0;i<liste.size();i++)
1997 sprintf(message,"%s %s",message,liste[i].c_str());
1998 system(message);
1999 delete taille;
2000 }
2001
2002 if (app.get_action()==115)
2003 {
2004 #ifdef ALL_OCC
2005 char fichier_entree[1000];
2006 char fichier_sortie[1000];
2007 int numgeo;
2008 double eps;
2009 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier MAGiC manquant");
2010 app.recupere_parametre_string(5,fichier_sortie,(char*)"Nom du fichier MAGiC manquant");
2011 app.recupere_parametre_double_avec_defaut(14,eps,0.01);
2012 app.recupere_parametre_int_avec_defaut(31,numgeo,-1);
2013 if (app.get_erreur()==true) return 0;
2014 char mess[1000];
2015 sprintf(mess,"Importation de la triangulation STL via OpenCASCADE V2017");
2016 app.affiche(mess);
2017 MG_FILE *gest = new MG_FILE(fichier_entree);
2018 MG_GEOMETRIE *mggeo;
2019 if(numgeo!=-1) mggeo = gest->get_mg_geometrieid(numgeo);
2020 else mggeo = gest->get_mg_geometrie(0);
2021 OCC_IMPORT occ_import;
2022 MG_MAILLAGE* tristl = occ_import.importer_triangulation_V2017(*(MG_GESTIONNAIRE*)gest,mggeo,eps,1);
2023 gest->enregistrer(fichier_sortie);
2024 delete gest;
2025 #endif
2026 }
2027
2028 if (app.get_action()==116)
2029 {
2030 #ifdef ALL_OCC
2031 char fichier_entree[1000];
2032 char fichier_sortie[1000];
2033 int numgeo;
2034 double eps;
2035 app.recupere_parametre_string(1,fichier_entree,(char*)"Nom du fichier MAGiC manquant");
2036 app.recupere_parametre_string(5,fichier_sortie,(char*)"Nom du fichier .brep manquant");
2037 app.recupere_parametre_double_avec_defaut(14,eps,0.01);
2038 app.recupere_parametre_int_avec_defaut(31,numgeo,-1);
2039 if (app.get_erreur()==true) return 0;
2040 char mess[1000];
2041 sprintf(mess,"Exportation du fichier .brep via OpenCASCADE V2017");
2042 app.affiche(mess);
2043 MG_FILE *gest = new MG_FILE(fichier_entree);
2044 MG_GEOMETRIE *mggeo;
2045 if(numgeo!=-1) mggeo = gest->get_mg_geometrieid(numgeo);
2046 else mggeo = gest->get_mg_geometrie(0);
2047 OCC_FONCTION_V2017* occ_fonction_v2017 = (OCC_FONCTION_V2017*)&mggeo->get_occ_fonction();
2048 occ_fonction_v2017->exporter_fichier_brep(fichier_sortie,mggeo);
2049 char message[5000];
2050 std::sprintf(message,"Creation du fichier : %s",fichier_sortie);
2051 app.affiche(message);
2052 delete gest;
2053 #endif
2054 }
2055 if (app.get_action()==117)
2056 {
2057 if (app.get_erreur()==true) return 0;
2058 char mess[2000];
2059 sprintf(mess,"Conversion MAGiC vers abaqus");
2060 app.affiche(mess);
2061 MG_FILE gest(nomfichier);
2062 int nb=gest.get_nb_fem_maillage();
2063 for (int i=0;i<nb;i++)
2064 {
2065 FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
2066 unsigned long id=mai->get_id();
2067 strncpy(chaine,nomfichier,p-nomfichier);
2068 char chaine2[3000];
2069 sprintf(chaine2,"%s%lu",chaine,id);
2070 std::string namefic=chaine2;
2071 MG_EXPORT exp;
2072 sprintf(mess," Maillage EF : Creation de %s.inp",namefic.c_str());
2073 app.affiche(mess);
2074 exp.abaqus(mai,namefic);
2075 }
2076 }
2077
2078
2079 #endif
2080 char mess[255];
2081 sprintf(mess,"Fin");
2082 app.affiche(mess);
2083 }
2084
2085
2086
2087
2088
2089
2090
2091
2092
2093
2094
2095
2096
2097
2098