Télécharger boobj.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : boobj.dgibi
  2. SAUT PAGE ;
  3. **********************************************************************
  4. * Utilisation des opérateurs CHI1 et CHI2
  5. * test avec échange
  6. * Ce test est identique à bo2.dgibi mais les entrées sont des OBJETS
  7. **********************************************************************
  8. *
  9. * DEFINITION DU MAILLAGE
  10. *
  11. *
  12. n1 = 1 ;
  13. *n2 = 160 ;
  14. *n2 = 80 ;
  15. n2 = 1 ;
  16. *
  17. * POINT SERVANT A DEFINIR LE CONTOUR
  18. *
  19. a = 0.0 0.0 ;
  20. b = 1. 0.0 ;
  21. c = 1. 8. ;
  22. d = 0. 8. ;
  23. *
  24. option elem qua4 ;
  25. *
  26. ab = a droit n1 b ;
  27. bc = b droit n2 c ;
  28. cd = c droit n1 d ;
  29. da = d droit n2 a ;
  30. *
  31. * DEFINITION DU MAILLAGE
  32. *
  33. *
  34. GP = AB BC CD DA DALL 'PLAN' ;
  35. ELIM 0.001 GP ;
  36. AB= CHANG POI1 AB ;
  37. G1= CHANG POI1 GP ;
  38. G0= DIFF G1 AB ;
  39. *
  40. TABDON=OBJET DONCHI1 ;
  41. TABDON%GIDEN (LECT 1 22 5 103 50);
  42. COMP1=OBJET LINVCOMP ;
  43. COMP1%COM_IDEN 162 ;
  44. COMP1%COM_NOM X162 ;
  45. COMP1%COM_CHAR -1 ;
  46. TABDON%GNVCOMP 1 COMP1 ;
  47. TABESP1= OBJET LIESPECE ;
  48. TABESP1%ESP_IDEN 163 ;
  49. TABESP1%ESP_LOGK 7. ;
  50. TABESP1%ESP_ITYP 2 ;
  51. TABESP1%ESP_COMP (LECT 162 50) ;
  52. TABESP1%ESP_STOE (PROG 1. 1.) ;
  53. TABESP2= OBJET LIESPECE;
  54. TABESP2%ESP_IDEN 164 ;
  55. TABESP2%ESP_LOGK 1.5 ;
  56. TABESP2%ESP_ITYP 2 ;
  57. TABESP2%ESP_COMP ( LECT 162 5) ;
  58. TABESP2%ESP_STOE ( PROG 1. 1.) ;
  59. *TABDON.NVESP= TABLE ;
  60. TABDON%GNVESP 1 TABESP1 ;
  61. TABDON%GNVESP 2 TABESP2 ;
  62. TABDON%GECHANGE ( LECT 162) ;
  63. *
  64. TB1=CHI1 TABDON COMP '/u/castem/divers/COMPOM'
  65. LOGK '/u/castem/divers/COMPOM' ;
  66. *
  67. * Table de données de CHI2
  68. *
  69. *
  70. *TBPAR2= TABLE ;
  71. *TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  72. TBPAR2= OBJET DONCHI2 ;
  73. TBPAR2%GLOGC (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%GTOT ( TOTCA + TOTLI + TOTNA + TOTCL + TOTSF +TOTH) ;
  92. TBPAR2%GFIONI ( MANU CHPO G1 1 SCAL 0.001) ;
  93. TBPARM= OBJET PARMCHI2 ;
  94. TBPARM%GITMAX 80;
  95. TBPARM%GEPS 1.D-8 ;
  96. TBPARM%GNFI 4 ;
  97. TBPARM%GITERSOL 15 ;
  98. TBPARM%GSORTIE ( 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.  
  134.  
  135.  

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