Télécharger la plaquette Cast3M


fr en

Télécharger crit_pplan.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : crit_pplan.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *---------------------------------------------------------------------*
  5. * Test du critere d'ecart a la planeite de l'operateur 'SURF', option *
  6. * 'PLAN' : (S1 = SURF GEO1 PLAN CRIT ;)                               *
  7. *---------------------------------------------------------------------*
  8. 'OPTI' 'DIME' 3 'ELEM' 'TET4' ;
  9. *                                                                     *
  10. * Mettre IG1 a VRAI pour trace :                                      *
  11. IG1      = FAUX ;
  12. *                                                                     *
  13. * Contour "test" :                                                    *
  14. P0       = 0. 0. 0. ;
  15. P1       = 1. 0. 0. ;
  16. P2       = 1. 1. 0. ;
  17. P3       = 0. 1. 0. ;
  18. N1       = 10 ;
  19. CNT1     = (P0 'DROI' N1 P1) 'ET' (P1 'DROI' N1 P2)
  20.   'ET' (P2 'DROI' N1 P3) 'ET' (P3 'DROI' N1 P0) ;
  21. 'DEPL' P2 'PLUS' (0. 0. 2.5) ;
  22. *                                                                     *
  23. * Calcul de la normale au contour selon meme methode que 'SURF' :     *
  24. NB1      = 'NBNO' CNT1 ;
  25. PN1      = P0 ;
  26. PI1      = CNT1 'POIN' NB1 ;
  27. 'REPE' BOU1 NB1 ;
  28.   PI2      = CNT1 'POIN' &BOU1 ;
  29.   PN1      = PN1 'PLUS' (PI1 'PVEC' PI2) ;
  30.   PI1      = PI2 ;
  31. 'FIN' BOU1 ;
  32. PN1      = PN1 / ('NORM' PN1) ;
  33. *                                                                     *
  34. * Calcul du critere de planeite :                                     *
  35. PG1      = 'BARY' CNT1 ;
  36. XTEST1   = ('PSCA' (P2 'MOIN' PG1) PN1) / ('NORM' (P2 'MOIN' PG1)) ;
  37. *                                                                     *
  38. S1       = CNT1 'SURF' 'PLAN' (1.0001 * XTEST1) ;
  39. *                                                                     *
  40. 'SI' IG1 ;
  41.   'TRAC' S1 'TITR' 'Critere d"ecart a la planeite de "SURF" OK !' ;
  42. 'FINS' ;
  43. *                                                                     *
  44. 'FIN' ;
  45.  
  46.  
  47.  
  48.  
  49.  

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