Télécharger vsur1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : vsur1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Vsur1.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. *******************************************************
  8. * vsur1.dgibi : Test de l'opérateur VSUR en coq2 *
  9. *******************************************************
  10. opti echo 0 dime 2 elem seg2 mode axis;
  11. *
  12. * Définition de la géométrie
  13. *
  14. p0 = 0. 0.;
  15. p1 = 10. 10.;
  16. pc = 0. 10.;
  17. obj1 = cerc 10 p0 pc p1;
  18. cl1 = bloq 'DEPL' 'ROTA' p1;
  19. cl2 = bloq 'UZ' p0;
  20. *
  21. * Création du modèle
  22. *
  23. mod = modeli obj1 mecanique elastique isotrope coq2;
  24. *
  25. * Comparaison des champs obtenus par VSUR
  26. *
  27. * che1 : champ de vecteurs surfaciques
  28. che1 = VSUR mod;
  29. * chn1 : champ de normaes
  30. chn1 = VSUR mod 'NORM';
  31. * chj1 : champ des jacobiens
  32. chj1 = JACO mod;
  33. cnr1 = exco 'VR' che1 'SCAL';
  34. cnz1 = exco 'VZ' che1 'SCAL';
  35. * champ des normes
  36. cnn1 = (cnr1*cnr1) + (cnz1*cnz1);
  37. cjn1 = chj1*chj1;
  38. cze1 = cnn1 - cjn1;
  39. *list cze1;
  40. cpz1 = chan 'CHPO' mod cze1;
  41. res1 = maxi cpz1 'ABS';
  42. *list res1;
  43. cnr1 = exco 'VR' chn1 'SCAL';
  44. cnz1 = exco 'VZ' chn1 'SCAL';
  45. cnn1 = (cnr1*cnr1) + (cnz1*cnz1);
  46. *list cnn1;
  47. cpn1 = chan 'CHPO' mod cnn1;
  48. res2 = mini cpn1 'ABS';
  49. res3 = maxi cpn1 'ABS';
  50. *list res2;
  51. *list res3;
  52. EC1 = RES3 - RES2;
  53. SI ((res1 < 1.e-6) ET (ABS(res2 - 1.) < 1.E-6) ET (EC1 < 1.E-6)) ;
  54. MESS 'OPERATEUR <VSUR> ERR 0';
  55. ERRE 0;
  56. SINON;
  57. MESS 'OPERATEUR <VSUR> ERR 5';
  58. ERRE 5;
  59. FINSI;
  60. fin;
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  

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