Télécharger la plaquette Cast3M


fr en

Télécharger bo2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  bo2.dgibi
  2.  SAUT PAGE ;
  3. **********************************************************************
  4. *                Utilisation des opérateurs CHI1 et CHI2
  5. *               test avec échange
  6. **********************************************************************
  7. * repertoire des fichiers "divers"
  8. DIVERS = VENV 'CASTEM_DIVERS';
  9. *
  10. OPTION DIME 2 ;                                                         
  11. *
  12. *                DEFINITION DU MAILLAGE
  13. *
  14. *
  15. n1 = 1 ;
  16. *n2 = 160 ;
  17. *n2 = 80  ;
  18. n2 = 1   ;
  19. *
  20. *                 POINT SERVANT A DEFINIR LE CONTOUR
  21. *
  22. a = 0.0 0.0 ;
  23. b = 1. 0.0 ;
  24. c = 1. 8. ;
  25. d = 0. 8. ;
  26. *
  27. option elem qua4 ;
  28. *
  29. ab = a droit n1 b ; 
  30. bc = b droit n2 c ; 
  31. cd = c droit n1 d ; 
  32. da = d droit n2 a ; 
  33. *
  34. *                         DEFINITION DU MAILLAGE
  35. *
  36. *
  37. GP = AB BC CD DA DALL 'PLAN' ;
  38. ELIM 0.001 GP ;
  39. AB= CHANG POI1 AB ;
  40. G1= CHANG POI1 GP ;
  41. G0= DIFF G1 AB ;
  42. *
  43. TABDON=TABLE ;                                                          
  44. TABDON.IDEN= LECT 1 22 5 103 50 ;
  45. TABDON.NVCOMP= TABLE ;
  46.  COMP1=TABLE ;
  47.  COMP1.IDEN= 162 ;
  48.  COMP1.NOM= X162 ;
  49.  COMP1.CHARGE= -1 ;
  50. TABDON.NVCOMP.1= COMP1 ;
  51.    TABESP1= TABLE   ;
  52.     TABESP1.IDEN= 163 ;
  53.     TABESP1.LOGK= 7. ;
  54.     TABESP1.ITYP= 2 ;
  55.     TABESP1.COMP=  LECT 162 50 ;
  56.     TABESP1.STOECH=  PROG  1.  1.  ;
  57.    TABESP2= TABLE   ;
  58.     TABESP2.IDEN= 164 ;
  59.     TABESP2.LOGK= 1.5 ;
  60.     TABESP2.ITYP= 2 ;
  61.     TABESP2.COMP=  LECT 162 5 ;
  62.     TABESP2.STOECH=  PROG  1. 1.   ;
  63.  TABDON.NVESP= TABLE ;
  64.  TABDON.NVESP.1= TABESP1 ;
  65.  TABDON.NVESP.2= TABESP2 ;
  66.  TABDON.ECHANGE= LECT 162 ;
  67. *
  68. *  LOGK  '/u2/anne/mch/LOGK2' ;
  69.  TB1=CHI1 TABDON  COMP  ('CHAINE' DIVERS '/COMPOM')
  70.    LOGK  ('CHAINE' DIVERS '/COMPOM') ;
  71. *
  72. *                        Table de données de CHI2
  73. *
  74. *
  75.  TBPAR2= TABLE ;
  76.  TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  77.  TBPAR2.LOGC= MANU CHPO G1 6 X001 -6. X022 -6. X005 -6.
  78.                X050 -6. X162 -2. X103 -6. ;
  79.   TOTCA= MANU CHPO G1 1 X001 1.D-15 ;
  80.   TOTLI1= MANU CHPO G1 1 X022 1.248D-9 ;
  81.   TOTLI2= MANU CHPO AB 1 X022 (1.D-3 - 1.248D-9) ;
  82.   TOTLI = TOTLI1 + TOTLI2 ;
  83.   TOTNA1= MANU CHPO G1 1 X005 1.100131015D-1 ;
  84.   TOTNA2= MANU CHPO AB 1 X005 (1.D-1 - 1.100131015D-1) ;
  85.   TOTNA= TOTNA1 + TOTNA2 ;
  86.   TOTCL1= MANU CHPO G1 1 X103 1.D-1 ;
  87.   TOTCL2= MANU CHPO AB 1 X103 (1.02D-1 - 1.D-1)  ;
  88.    TOTCL= TOTCL1+ TOTCL2 ;
  89.   TOTSF1= MANU CHPO G1 1 X162  1.D-2 ;
  90.   TOTSF2= MANU CHPO AB 1 X162  (1.D-15 - 1.D-2);
  91.   TOTSF= TOTSF1 + TOTSF2 ;
  92.   TOTH1= MANU CHPO G1 1 X050 -1.3012D-5 ;
  93.   TOTH2= MANU CHPO AB 1 X050 (1.0001D-3 + 1.3012D-5) ;
  94.    TOTH = TOTH1 + TOTH2 ;
  95.  TBPAR2.TOT= TOTCA + TOTLI + TOTNA + TOTCL + TOTSF +TOTH ;
  96.  TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.001 ;
  97.  TBPARM= TABLE ;
  98.  TBPARM.ITMAX = 80;
  99.  TBPARM.EPS= 1.D-8 ;
  100.  TBPARM.NFI= 4 ;
  101.  TBPARM.ITERSOLI=15 ;
  102.  TBPARM.SORTIE= MOTS 'FION' 'TYP5' 'SURF' ;
  103. *
  104.  TB3= CHI2 TB1 TBPARM TBPAR2 ;
  105. *
  106. *                                   controle des résultats
  107. *
  108.  FIONTE1=MANU CHPO G1 1 SCAL 1.00030E-01 'NATURE' DISCRET ;
  109.  FIONTE2=MANU CHPO AB 1 SCAL (1.02000E-01 - 1.00030E-01)
  110.  'NATURE' DISCRET ;
  111.  FIONTES= FIONTE1 + FIONTE2 ;
  112.  VERR1= ( ABS ( FIONTES - TB3.FION )) MASQUE SUPERIEUR 5.D-6 SOMME ;
  113.  SURFTE1= MANU CHPO G0 2  W006  9.98262E-03  W010  1.73813E-05
  114.    'NATURE' DISCRET ;
  115.  SURFTE2= MANU CHPO AB 2 W006  3.15844E-19  W010  9.98885E-16
  116.  'NATURE' DISCRET ;
  117.  SURFTES= SURFTE1+SURFTE2 ;
  118.  SURFD=SURFTES /  50. ;
  119.  VERR2= ( ABS ( SURFTES - TB3.SURF )) MASQUE SUPERIEUR SURFD SOMME ;
  120.  TY5TE1=  MANU CHPO G0 1 W011  3.12037E-20 'NATURE' DISCRET ;
  121.  TY5TE2=  MANU CHPO AB 1 W011   9.43503E-33  'NATURE' DISCRET ;
  122.  TY5TES= TY5TE1+TY5TE2 ;
  123.  TY5D= TY5TES/ 50. ;
  124.  VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  125. *
  126.  VERR= VERR1+VERR2+VERR3 ;
  127.  SI (VERR EGA 0 ) ;
  128.   ERRE 0 ;
  129.  SINO  ;
  130.   ERRE 5 ;
  131.  FINSI ;
  132. FIN ;
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  

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