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

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