Télécharger la plaquette Cast3M


fr en

Télécharger varihc.procedur

Retour à la liste

Numérotation des lignes : 

  1. * VARIHC    PROCEDUR  FANDEUR   22/06/02    21:15:07     11372          
  2. ****************************************************
  3. * Procedure de tracking par resolution d'un probleme
  4. * de conduction thermique stationnaire
  5. * Application au modele E-FEM
  6. * Version  : 
  7. *     Variable online global tracking 
  8. * Creation : 
  9. *     Francesco RICCARDI 22/08/2016
  10. *     Benjamin RICHARD
  11. * Contact :  
  12. *     Francesco.Riccardi[at]cea.fr
  13. *     Benjamin.Richard[at]cea.fr 
  14. * Institution : 
  15. *     CEA\DEN\DANS\DM2S\SEMT\EMSI
  16. ****************************************************
  17. * Voir notice TRACKING pour plus d'informations
  18. ****************************************************
  19. *
  20. *******----DEBUT Procedure VARIHC----*******
  21. *
  22. DEBP VARIHC TAB2*'TABLE';
  23.  
  24.  
  25. * Indice d'evolution
  26. NPAS = dime tab2. deplacements;
  27. NPAS = NPAS - 1;
  28. *LIST NPAS;
  29.  
  30. * Modeles
  31. MODEL1 = TAB2. 'MODELE';
  32. MODELS = EXTR MODEL1 'ZONE';
  33.  
  34. * Modele et materiau EFEM
  35. MODEL11 = MODELS. 1;
  36. SURFACE1 = MODELS. 2;
  37. MATr1 = TAB2. 'WTABLE'. 'CARACTERISTIQUES';
  38. MATr11 = REDU MATr1 SURFACE1;
  39. E11 = EXCO MATr11 YOUN;
  40. NU11 = EXCO MATr11 NU;
  41. P11 = EXCO MATr11 FT;
  42. P12 = EXCO MATR11 XNX;
  43. P13 = EXCO MATR11 XNY;
  44. P14 = EXCO MATr11 IND1;
  45. P15 = EXCO MATR11 XIE;
  46. P16 = EXCO MATR11 YIE;
  47.  
  48. * Coordonnees de la geometrie
  49. X0 = COOR 1 SURFACE1;
  50. Y0 = COOR 2 SURFACE1;
  51. X0S = CHAN 'CHAM' X0 MODEL11 'STRESSES';
  52. X0S = CHAN 'TYPE' X0S 'CONTRAINTES PRINCIPALES';
  53. X0S = CHAN 'COMP' 'SI11' X0S;
  54. Y0S = CHAN 'CHAM' Y0 MODEL11 'STRESSES';
  55. Y0S = CHAN 'TYPE' Y0S 'CONTRAINTES PRINCIPALES';
  56. Y0S = CHAN 'COMP' 'SI11' Y0S;
  57.  
  58. * Nombre de root elements
  59. NROOTS = TAB2. 'NROOTS'; 
  60. *LIST NROOTS;
  61.  
  62. * Nombre maximum de fissures
  63. MAXR = TAB2. 'NRMAX'; 
  64.  
  65. * Variables internes : NFLA
  66. FLAG = EXCO (TAB2. 'VARIABLES_INTERNES'. NPAS) 'NFLA';
  67. *TRAC FLAG MODEL11;
  68. VFLAG = MAXI FLAG;
  69. ELACT1 = FLAG ELEM 'EGALE' 1.;
  70.  
  71. * Maximum des contraintes principales
  72. SIG1 = TAB2. 'CONTINUATION'. 'CONTRAINTES';
  73. PRIN1 = PRIN SIG1 MODEL11;
  74. SIGMA11 = EXCO 'SI11' PRIN1;
  75. COS1X = NOMC 'SI11' (EXCO PRIN1 'COX1');
  76. COS1Y = NOMC 'SI11' (EXCO PRIN1 'COY1');
  77. SIGMAX = MAXI SIGMA11;
  78.  
  79. * 
  80. *----- Debut algorithme de tracking --------------------------- 
  81. * 
  82. SI ((SIGMAX > 0.) ET (NPAS < (TAB2. 'NPAS_TRACKING')));
  83.  
  84.  
  85. * Elements qui violent le critere de resistance
  86. ***dans la phase elastique lineaire VFLAG=0***
  87. SI (VFLAG EGA 0.);
  88.     ELSP11 = SIGMA11 ELEM 'MAXI';
  89. SINON;
  90.     MESH0 = TAB2. 'MESH'. (NPAS - 1);
  91.     SURFRED = SURFACE1 DIFF MESH0;
  92.     SIGRED = REDU SIGMA11 SURFRED;
  93. ***pour assurer que les points test n'appartiennent pas a MESH0***  
  94. ***cela peut se passer si on a des bloccages***   
  95.     ELSP11 = SIGRED ELEM 'EGSUPE' (MAXI P11);
  96. FINSI;
  97. NTEST = NBEL ELSP11;
  98. SI (NTEST > 0);
  99.     SP11 = REDU SIGMA11 ELSP11;
  100. FINSI;
  101.  
  102. TABISO0 = TABLE;
  103. TABROOT = TABLE;
  104. TABEROOT = TABLE;
  105. SI (VFLAG EGA 0);
  106.     TABEROOT. 1 = SIGMA11 ELEM 'MAXI';    
  107.     TABROOT. 1 = BARY (TABEROOT. 1);
  108.     SI (NROOTS EGA 1);
  109.         TABROOT. 1 = TAB2. 'PROOTS'. 1;
  110.         TABEROOT. 1 = TAB2. 'EROOTS'. 1;
  111.     FINSI;
  112. SINON;
  113.     i = 0;
  114. ***NROOTS est surement >0***
  115.     REPETER BOUCLE1 NROOTS;
  116.         indice = i * 4;
  117.         LXIN0 = TAB2. 'TABISO'. (indice + 1);
  118.         LYIN0 = TAB2. 'TABISO'. (indice + 2);
  119.         LXOUT0 = TAB2. 'TABISO'. (indice + 3);
  120.         LYOUT0 = TAB2. 'TABISO'. (indice + 4);
  121.         MESH0 = TAB2. 'MESH_ISO'. (i + 1);
  122. * Nombre de points qui composent la fissure  
  123.         DIMM0 = DIME LXIN0;     
  124. ***si tous les elements sont actives on touche rien*** 
  125.         FLAGISO = REDU FLAG MESH0;
  126.         MESHACT = FLAGISO ELEM 'SUPERIEUR' 0.;
  127.         NACT = NBEL MESHACT;   
  128.         SI ((NACT > 0) ET (NACT < DIMM0));
  129.             SI ((NACT + 1) EGA DIMM0);
  130.                 LISTEN0 = DIMM0;
  131.             SINON;
  132.                 LISTEN0 = ENTI (PROG (NACT + 1) PAS 1 DIMM0);
  133.             FINSI;
  134. * Partie activee de la fissure            
  135.             LXIN0 = ENLE LXIN0 LISTEN0;
  136.             LYIN0 = ENLE LYIN0 LISTEN0;
  137.             LXOUT0 = ENLE LXOUT0 LISTEN0;
  138.             LYOUT0 = ENLE LYOUT0 LISTEN0;
  139.         FINSI;
  140.         LXG0 = (LXIN0 + LXOUT0) / 2.; 
  141.         LYG0 = (LYIN0 + LYOUT0) / 2.;  
  142.         TABROOT. (i + 1) = TAB2. 'PROOTS'. (i + 1);
  143.         TABEROOT. (i + 1) = TAB2. 'EROOTS'. (i + 1);
  144. ***le point root est le point output du dernier element active***
  145. ***si NACT=0 il coincide avec celui du pas precedent***           
  146.         SI (NACT > 0);
  147.             TABROOT. (i + 1) = (EXTR LXOUT0 NACT) (EXTR LYOUT0 NACT);
  148.             TABEROOT. (i + 1) = MESH0 ELEM 'CONTENANT'
  149.             ((EXTR LXG0 NACT) (EXTR LYG0 NACT));
  150.         FINSI;
  151.         TABISO0. (indice + 1) = LXIN0;
  152.         TABISO0. (indice + 2) = LYIN0;
  153.         TABISO0. (indice + 3) = LXOUT0;
  154.         TABISO0. (indice + 4) = LYOUT0;
  155.         i = i + 1;
  156.     FIN BOUCLE1;
  157. FINSI;    
  158.  
  159. SI (NTEST > 0);
  160. * Points root potentiels ordonnes par contrainte decroissante
  161.     TABTEST = TABLE;
  162.     i = 0;
  163.     LSP11  = 'PROG' ;
  164.     LPTEST = 'PROG' ;
  165.     REPETER BOUCLE2 NTEST;
  166.         i = i + 1;
  167.         STEST = EXTR SP11 'SI11' 1 i 1;
  168.         LSP11 = LSP11 ET STEST ;
  169.         LPTEST = LPTEST ET i ;
  170.     FIN BOUCLE2;
  171.     TABTEST. 1 = LPTEST;
  172.     TABTEST. 2 = LSP11;
  173.     TABTEST = ORDO TABTEST 2 'DECR';
  174. FINSI;
  175.  
  176. *
  177. *\\\\ Resolution du probleme de conduction ////*
  178. *
  179. * Cosinus des tangentes
  180. COS2X = EXCO PRIN1 'COX2';
  181. COS2Y = EXCO PRIN1 'COY2';
  182. COS2Z = EXCO PRIN1 'COZ2';
  183.  
  184. * Champ vectoriel des tangentes
  185. *SIGMA22 = MANU 'CHML' MODEL11 'SI22' 1. 'TYPE' 
  186. *'CONTRAINTES PRINCIPALES' 'STRESSES';
  187. *PRIN12 = (SIGMA22 ET COS2X ET COS2Y ET COS2Z);
  188. *VECTAN = VECT PRIN12 MODEL11 0.001 'SI22';    
  189.  
  190. * Modele thermique
  191. MOTH = MODE SURFACE1 THERMIQUE ANISOTROPE;
  192. * Perturbation conductivite numerique
  193. TTPETIT = 1.e-6;
  194. * Tenseur conductivite anisotrope
  195. TTX = NOMC 'UX' COS2X;
  196. TTY = NOMC 'UX' COS2Y;
  197. KXX = (TTX ** 2) + ((TTX ** 0) * TTPETIT);
  198. KYY = (TTY ** 2) + ((TTY ** 0) * TTPETIT);
  199. KXY = TTX * TTY;
  200. PX = 1.0 0.;
  201. MATH = MATE MOTH 'DIRECTION' PX K11 KXX K21 KXY K22 KYY;
  202. * Matrice de rigidite
  203. KTH = COND MOTH MATH;
  204. * Conditions aux limites
  205. PCL1 = CHAN 'POI1' (TAB2. 'LBC'); 
  206. CLT1 = BLOQ PCL1 'T';
  207. CLTH = CLT1;
  208. * Temperature imposee
  209. VBC = TAB2. 'BCSTH';
  210. DE1 = DEPI CLT1 VBC;
  211. DETH = DE1;
  212. * Assemblage
  213. KTHTOT = KTH ET CLTH;
  214. * Resolution
  215. RESTH = RESO KTHTOT DETH;
  216. T0 = CHAN 'CHAM' RESTH MOTH 'NOEUD';
  217. T0 = EXCO T0 'T';
  218. TAB2. 'RESTH'. NPAS = RESTH;   
  219.  
  220. *    
  221. *\\\\ Extraction des isovaleurs ////*   
  222. *
  223. SI (NROOTS EGA MAXR);
  224.     NTEST = 0;
  225. FINSI;
  226. TISOPAS = TABLE;
  227. TABMESH1 = TABLE;
  228. MESH1 = VIDE 'MAILLAGE';
  229. ***on choisit un point du domaine en faisant l'hypothese***
  230. ***que aucune fissure passe par ceci***
  231. P0 = MANU 'POI1' ((X0 * Y0) POIN 'MINI');
  232. i = 0;
  233. REPETER BOUCLE3 (NROOTS + NTEST);          
  234.     i = i + 1;
  235. * Fissure activees aux pas precedents   
  236.     SI (i &lt;EG NROOTS);    
  237.         PPA = TABROOT. i;
  238.         MPPA = MANU 'POI1' PPA;
  239.         T0P = PROI T0 MPPA; 
  240.         TG = EXTR T0P 'T' PPA;
  241.         ISOVG = ISOV T0 TG;
  242. *        TRAC (ISOVG ET (CONT SURFACE1));
  243. ***cette phase est critique***
  244. ***si on perd l'unicite de la solution ORDO ne marche pas***
  245. ***arreter la procedure de tracking***
  246.         ISOVG = ORDO ISOVG;         
  247.                 T0G = MANU 'CHML' MOTH 'T' TG 'NOEUD';
  248.                 T01 = T0 - T0G;
  249.                 T02 = T01 MASQ 'SUPERIEUR' 0.;
  250.                 MAILA = T02 ELEM 'EGALE' 1.;
  251.                 MAILB = T02 ELEM 'EGALE' 0.; 
  252.                 MESHISO = SURFACE1 DIFF (MAILA ET MAILB);
  253. * Extraction du maillage  
  254.                 MM1 = (CHAN 'POI1' MESHISO) ET P0;
  255.                 MMT = (CHAN 'POI1' MESH1) ET P0;        
  256.                 CHECK = MMT INTE MM1;
  257.                 SI ((NBEL CHECK) > 1);
  258.             ISOVG = TAB2. 'ISOTEST'. i;
  259.                 FINSI;
  260.         MESH1 = MESH1 ET MESHISO;
  261.         TISOPAS. i = ISOVG;
  262.         TABMESH1. i = MESHISO;
  263.         TAB2. 'ISOTEST'. i = ISOVG;       
  264. * Nouvelles fissures
  265. ***on fait l'hypothese d'avoir une seule nouvelle fissure***
  266. ***par pas de temps***        
  267.     SINON;
  268.         iTEST = ENTI (EXTR (TABTEST. 1) (i - NROOTS));
  269.         ETEST = ELSP11 ELEM iTEST;
  270.         PPP = BARY ETEST;
  271.         MPPP = MANU 'POI1' PPP;
  272.         T0P = PROI T0 MPPP; 
  273.         TG = EXTR T0P 'T' PPP;
  274.         ISOVG = ISOV T0 TG;        
  275. * Extraction du maillage 
  276.                 T0G = MANU 'CHML' MOTH 'T' TG 'NOEUD';
  277.                 T01 = T0 - T0G;
  278.                 T02 = T01 MASQ 'SUPERIEUR' 0.;
  279.                 MAILA = T02 ELEM 'EGALE' 1.;
  280.                 MAILB = T02 ELEM 'EGALE' 0.;
  281.                 MESHISO = SURFACE1 DIFF (MAILA ET MAILB);        
  282.                 MM1 = (CHAN 'POI1' MESHISO) ET P0;
  283.                 MMT = (CHAN 'POI1' MESH1) ET P0;        
  284.                 CHECK = MMT INTE MM1;        
  285. * On verifie s'il y a des intersection entre maillages
  286. ***dans ce cas la on ecarte l'isovaleur***        
  287.         SI ((NBEL CHECK) EGA 1);
  288.             ISOVG = ORDO ISOVG;
  289. *            TRAC (ISOVG ET (CONT SURFACE1));
  290.             T0RED = REDU T0 ETEST;
  291.             ISOROOT = ISOV T0RED TG;
  292.             PROOT = ISOROOT POIN 1;
  293.             CHECK1 = SURFACE1 ELEM 'CONTENANT' PROOT 'TOUS';
  294. ***on fait l'hypothese que le root element appartient au contour***            
  295.             SI ((NBEL CHECK1) > 1);
  296.                 PROOT = ISOROOT POIN 2;
  297.             FINSI;
  298.             TABROOT. (NROOTS + 1) = PROOT;
  299.                         TABEROOT. (NROOTS + 1) = ETEST; 
  300.             TISOPAS. (NROOTS + 1) = ISOVG;
  301.             TAB2. 'ISOTEST'. (NROOTS + 1) = ISOVG;
  302.             TABMESH1. (NROOTS + 1) = MESHISO;
  303.             TAB2. 'NROOTS' = (TAB2. 'NROOTS') + 1;
  304.             QUIT BOUCLE3;
  305.                 FINSI;
  306.     FINSI;    
  307. FIN BOUCLE3;    
  308. TAB2. 'PROOTS' = TABROOT;
  309. TAB2. 'EROOTS' = TABEROOT;
  310.  
  311. ***cette partie sert a creer un systeme d'abscisses curvilignes***
  312. ***avec sense positif selon celui de propagation des fissures***
  313. ***l'origine est en pointe de fissure***
  314. TABISO1 = TABLE;
  315. NN = DIME TABROOT;
  316. *LIST NN;
  317. i = 0;
  318. REPETER BOUCLE4 NN;
  319.     i = i + 1;
  320.     indice = (i - 1) * 4;
  321.     PROOT = TABROOT. i;    
  322.     EROOT = TABEROOT. i;
  323.     ISOVG = TISOPAS. i;
  324.     MESHISO = TABMESH1. i;
  325.     NUMEL = (NBNO ISOVG) - 1;
  326.     NACT = 0;
  327.     SI ((VFLAG > 0.) ET (i &lt;EG NROOTS));
  328.         FLAGISO = REDU FLAG (TAB2. 'MESH_ISO'. i);
  329.         MESHACT = FLAGISO ELEM 'SUPERIEUR' 0.;
  330.         NACT = NBEL MESHACT;
  331.     FINSI;
  332.     XISOVG = COOR 1 ISOVG;
  333.     YISOVG = COOR 2 ISOVG;
  334.     EVOLX = EVOL 'CHPO' XISOVG ISOVG;
  335.     EVOLY = EVOL 'CHPO' YISOVG ISOVG;
  336.     LISX = EXTR EVOLX 'SCAL';
  337.     LISY = EXTR EVOLY 'SCAL';
  338.     LISABSC = EXTR EVOLX 'ABSC';
  339.     NUMP = NBNO ISOVG;
  340.     LISP = PROG 1 PAS 1 NUMP;    
  341.     j = 0;
  342.     REPETER B_INT NUMP;
  343.         j = j + 1;
  344.         Pj = ISOVG POIN j;
  345.         SI ('&lt;EG' ('DIST' Pj PROOT) 1.E-10);
  346.             jroot = j;
  347. *            LIST jroot;
  348. *            LIST NUMP;
  349.             QUIT B_INT;
  350.         FINSI;
  351.     FIN B_INT;  
  352.     SROOT = EXTR LISABSC jroot;
  353.     LISABSC1 = LISABSC - ((LISP / LISP) * SROOT);
  354.     MM = MAXI LISABSC1;
  355.     TABABS = TABLE;
  356.     TABABS. 1 = LISABSC1;
  357.     TABABS. 2 = LISP;
  358.     TABABS. 3 = LISX;
  359.     TABABS. 4 = LISY;
  360.     SI (NACT EGA 0);
  361. ***si MM=0 il faut changer le signe***    
  362.         SI ((MM EGA 0.));
  363.             TABABS. 1 = -1. * LISABSC1;
  364.             TABABS = ORDO TABABS 1 'CROI';
  365.         FINSI;
  366. ***le dernier point est un point de sortie***        
  367.         LXIN = ENLE (TABABS. 3) (DIME LISX);
  368.         LYIN = ENLE (TABABS. 4) (DIME LISY);
  369. ***le premier point est un point d'entree***        
  370.         LXOUT = ENLE (TABABS. 3) 1;
  371.         LYOUT = ENLE (TABABS. 4) 1;
  372.         TABISO1. (indice + 1) = LXIN;
  373.         TABISO1. (indice + 2) = LYIN;
  374.         TABISO1. (indice + 3) = LXOUT;
  375.         TABISO1. (indice + 4) = LYOUT;        
  376.     SINON;
  377.         TABISO1. (indice + 1) = TABISO0. (indice + 1);
  378.         TABISO1. (indice + 2) = TABISO0. (indice + 2); 
  379.         TABISO1. (indice + 3) = TABISO0. (indice + 3);
  380.         TABISO1. (indice + 4) = TABISO0. (indice + 4);   
  381.         SI (NACT < NUMEL);
  382. ***on cherche le point correspondant a la pointe de fissure***        
  383.             PTEST = ISOVG POIN (jroot + 1);   
  384.             TEST = NBEL (EROOT ELEM 'CONTENANT' PTEST 'NOVERIF');
  385.             SI (TEST > 0);   
  386.                 TABABS. 1 = -1. * LISABSC1;
  387.                 jroot = (NUMP - jroot) + 1; 
  388.             FINSI;  
  389.             TABABS = ORDO TABABS 1 'CROI';
  390.             SI (jroot EGA 2);
  391.                 LISTEN = 1;
  392.             SINON;
  393.                 LISTEN = ENTI (PROG 1 PAS 1 (jroot - 1));
  394.             FINSI;
  395.             LISX1 = ENLE (TABABS. 3) LISTEN;
  396.             LISY1 = ENLE (TABABS. 4) LISTEN;
  397. * Partie potentielle de la fissure       
  398.             LXIN1 = ENLE LISX1 (DIME LISX1);
  399.             LYIN1 = ENLE LISY1 (DIME LISY1);  
  400.             LXOUT1 = ENLE LISX1 1;
  401.             LYOUT1 = ENLE LISY1 1;
  402. * On integre avec la partie activee         
  403.             TABISO1. (indice + 1) = 
  404.             TABISO0. (indice + 1) ET LXIN1; 
  405.             TABISO1. (indice + 2) = 
  406.             TABISO0. (indice + 2) ET LYIN1;
  407.             TABISO1. (indice + 3) = 
  408.             TABISO0. (indice + 3) ET LXOUT1;
  409.             TABISO1. (indice + 4) = 
  410.             TABISO0. (indice + 4) ET LYOUT1;
  411.         FINSI;
  412.     FINSI;        
  413. FIN BOUCLE4;
  414. TAB2. 'TABISO' = TABISO1;
  415.  
  416. *
  417. *\\\\ Creation des chamelem ////*       
  418. *
  419. ***CHM vaut sur les isovaleurs et 0 ailleurs***
  420. CHM = SIGMA11 * 0.;
  421. * Abscisses des points d'entree des fissures
  422. CHIX = CHM;
  423. * Ordonnees des points d'entree des fissures
  424. CHIY = CHM;
  425. * Composante x des normales
  426. CHNX = CHM;
  427. * Composante y des normales
  428. CHNY = CHM;
  429. TAB2. 'MESH'. NPAS = VIDE 'MAILLAGE';
  430. i = 0;
  431. REPETER BOUCLE5 NN;
  432.     i = i + 1;
  433.     indice = (i - 1) * 4;
  434.     MESHISO = TABMESH1. i;
  435.     NACT = 0;
  436.     SI ((VFLAG > 0.) ET (i &lt;EG NROOTS));
  437.         FLAGISO = REDU FLAG (TAB2. 'MESH_ISO'. i);
  438.         MESHACT = FLAGISO ELEM 'SUPERIEUR' 0.;
  439.         NACT = NBEL MESHACT;
  440.     FINSI;
  441.     LXIN = TABISO1. (indice + 1);
  442.     LYIN = TABISO1. (indice + 2);
  443.     LXOUT = TABISO1. (indice + 3);
  444.     LYOUT = TABISO1. (indice + 4);
  445.     LXG = (LXIN + LXOUT) / 2.; 
  446.     LYG = (LYIN + LYOUT) / 2.;
  447.     SI (NACT EGA 0);
  448.         MTEST = MESHISO;
  449.     SINON;
  450.         MDIFF = MESHISO DIFF MESHACT;
  451.         MTEST = MDIFF ET MESHACT; 
  452.     FINSI;
  453.     j = 0;
  454.     REPETER B_INT (DIME LXIN);
  455.         j = j + 1;
  456. * Coordonnees des points entree et sortie
  457.         XIj = EXTR LXIN j;
  458.         YIj = EXTR LYIN j;
  459.         XOj = EXTR LXOUT j;
  460.         YOj = EXTR LYOUT j;
  461. * Coordonnees des points moyens   
  462.         XGj = EXTR LXG j;
  463.         YGj = EXTR LYG j;
  464.         PGj = XGj YGj;    
  465. * Calcul de la tangente et de la normale a la fissure   
  466.         LL = (((XIj - XOj) ** 2) + ((YIj - YOj) ** 2)) ** 0.5;
  467.         SI (YIj > YOj);
  468.             TY = (YIj - YOj) / LL;
  469.             TX = (XIj - XOj) / LL;
  470.         SINON;
  471.             TY = (YOj - YIj) / LL;
  472.             TX = (XOj - XIj) / LL;
  473.         FINSI;
  474.         SI (YIj EGA YOj);
  475.             TY = 0.;
  476.             TX = -1.;
  477.         FINSI;
  478.         NX = TY;
  479.         NY = -1. * TX;    
  480. * Element contenant le point moyen
  481.         ELj = MTEST ELEM 'CONTENANT' PGj;
  482.         VALXj = MAXI (REDU X0S ELj);
  483.         VALYj = MAXI (REDU Y0S ELj);
  484. * Champ nul sauf dans ELi 
  485.         CHPX = X0S MASQ 'EGALE' VALXj;
  486.         CHPY = Y0S MASQ 'EGALE' VALYj;
  487.         CHP = CHPX * CHPY;
  488.  
  489. * Champs par elements pour le model EFEM
  490.         CHM = CHM + CHP;    
  491.         CHIX = CHIX + (CHP * XIj);
  492.         CHIY = CHIY + (CHP * YIj);
  493.         CHNX = CHNX + (CHP * NX);
  494.         CHNY = CHNY + (CHP * NY);
  495.  
  496. * Extraction de la mesh de la fissure potentielle
  497.         SI(j EGA 1);
  498.             LISO = DROI 1 (XIj YIj) (XOj YOj); 
  499.             MESHi = ELj;
  500.         SINON;
  501.             LISO = LISO ET (DROI 1 (XIj YIj) (XOj YOj));
  502.             MESHi = MESHi ET ELj;
  503.         FINSI;        
  504.     FIN B_INT; 
  505.     TAB2. 'ISOTOT'. i = LISO;
  506.     TAB2. 'MESH_ISO'. i = MESHi; 
  507.     TAB2. 'MESH'. NPAS = (TAB2. 'MESH'. NPAS) ET MESHi;     
  508.  
  509. FIN BOUCLE5;
  510. *TRAC CHM MODEL11;
  511.  
  512. *
  513. *\\\\ Mise a jour des parametres materiaux ////*        
  514. *
  515. MATr2 = MATE MODEL11 YOUN E11 NU NU11 FT P11 XNX CHNX
  516. XNY CHNY IND1 CHM XIE CHIX YIE CHIY;
  517. MATr1 = MATr2;
  518. TAB2. 'WTABLE'. 'CARACTERISTIQUES' = MATR1;
  519.  
  520.  
  521. FINSI;
  522.  
  523.  
  524. FINP MATR1 ;
  525.  
  526.  
  527.  
  528.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales