Télécharger la plaquette Cast3M


fr en

Télécharger vsur3.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  vsur3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *********************************************************
  5. *  vsur3.dgibi : Test de VSUR en coq6 et coq8           *
  6. *********************************************************
  7. opti echo 0 dime 3 elem qua8;
  8. *
  9. *  Définition de la géométrie
  10. *
  11. p0 = 0. 0. 0.;
  12. p1 = 9. 0. 6.;
  13. p2 = 7. 6. 10.;
  14. p3 = -2. 6. 4.;
  15. p4 = 8. 3. 8.;
  16. p5 = -1. 3. 2.;
  17. l1 = D 3 p0 p1;
  18. l2 = D 2 p1 p2;
  19. l3 = D 5 p2 p3;
  20. l4 = D 4 p3 p0;
  21. *
  22. cn1 = l1 et l2 et l3 et l4;
  23. obj = SURF cn1 'PLANE';
  24. *
  25. *  Définition des modèles
  26. *
  27. obj6 = obj elem 'TRI6';
  28. obj8 = obj elem 'QUA8';
  29. mo6 =  modeli obj6  mecanique elastique isotrope coq6;
  30. mo8 =  modeli obj8  mecanique elastique isotrope coq8;
  31. mod1 = mo6 et mo8;
  32. *
  33. *  Comparaison des champs obtenus par VSUR  
  34. *
  35. * Premier modèle : coq6 et coq8
  36. che1 = VSUR mod1;
  37. chn1 = VSUR mod1 'NORM';
  38. chj1 = JACO mod1;
  39. cnx1 = exco 'VX' che1 'SCAL';
  40. cny1 = exco 'VY' che1 'SCAL';
  41. cnz1 = exco 'VZ' che1 'SCAL';
  42. cnn1 = (cnx1*cnx1) + (cny1*cny1) + (cnz1*cnz1);
  43. cjn1 = chj1*chj1;
  44. cze1 = cnn1 - cjn1;
  45. *list cze1;
  46. cnx1 = exco 'VX' chn1 'SCAL';
  47. cny1 = exco 'VY' chn1 'SCAL';
  48. cnz1 = exco 'VZ' chn1 'SCAL';
  49. cnn1 = (cnx1*cnx1) + (cny1*cny1) + (cnz1*cnz1);
  50. *list cnn1;
  51. CPZ1 = CHAN CHPO MOD1 CZE1;
  52. CPN1 = CHAN CHPO MOD1 CNN1;
  53. RE11 = MAXI CPZ1 'ABS';
  54. RE12 = (MAXI CPN1 'ABS') - 1.;
  55. RE13 = (MAXI CPN1 'ABS') - (MINI CPN1 'ABS');
  56. RES1 = RE11 + RE12 + RE13;
  57. SI (res1 < 1.e-6) ;
  58.   MESS 'OPERATEUR <VSUR> ERR  0';
  59.   ERRE  0;
  60. SINON;
  61.   MESS 'OPERATEUR <VSUR> ERR  5';
  62.   ERRE  5;
  63. FINSI; 
  64. fin;
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  

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