Télécharger test_coupe.dgibi

Retour à la liste

Numérotation des lignes :

  1. ** Un petit cas test de l'operateur COUP
  2.  
  3. ** Options generales
  4. OPTI 'DIME' 3 'ELEM' 'CUB8' ;
  5. itrac = FAUX ;
  6.  
  7. ** Maillage
  8. p1 = 0. 0. 0. ;
  9. p2 = 1. 0. 0. ;
  10. p3 = 1. 1. 0. ;
  11. p4 = 0. 1. 0. ;
  12. d1 = DROI 2 p1 p2 ;
  13. d2 = DROI 2 p2 p3 ;
  14. d3 = DROI 2 p3 p4 ;
  15. d4 = DROI 2 p4 p1 ;
  16. s = SURF (d1 ET d2 ET d3 ET d4) 'PLAN' ;
  17. v = s VOLU TRAN 10 (0. 0. 10.) ;
  18.  
  19. ** Plan de coupe et coupe du maillage
  20. pc1 = 0. 0. 3.1 ;
  21. pc2 = 1. 0. 3.1 ;
  22. pc3 = 0. 1. 4.3 ;
  23. pc = pc1 ET pc2 ET pc3 ;
  24. s1 = COUP v pc1 pc2 pc3 ;
  25.  
  26. ** Controle du maillage de la section
  27. * verification avec VERM
  28. VERM s1 ;
  29. * les points sont ils bien dans le plan ?
  30. ps1 = s1 POIN 'PLAN' pc1 pc2 pc3 1.E-10 ;
  31. SI ((NBNO ps1) NEG (NBNO s1)) ;
  32. ERRE 'Les points du maillage de la coupe ne sont pas dans le plan !' ;
  33. FINSI ;
  34. * trace du maillage
  35. SI itrac ;
  36. TRAC ((ARET v) ET (pc COUL 'ROUG') ET (s1 COUL 'VERT')) ;
  37. FINSI ;
  38.  
  39.  
  40. FIN ;
  41.  
  42.  
  43.  

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