Télécharger kp2_test.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : kp2_test.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************
  5. *
  6. * test de la matrice de rigidité associée à un champ
  7. * de pression linéaire.
  8. * Il s'agit de calculer la variation de la composante
  9. * verticale des forces de pression d'un reservoir
  10. * semispherique, due à un mouvement verticale dans un champ
  11. * de pression hydrostatique.
  12. *
  13. *
  14. * reference : rapport DMT/95/
  15. *
  16. ******************************************************
  17. r1 = 1. ; comm 'rayon de la sphere' ;
  18. d1 = 1. ; comm 'deplacement verticale' ;
  19. rg = 1. ; comm 'module du gradient de la pression' ;
  20. ******************************************************
  21. opti dime 3 elem qua4 ;
  22. p1 = 0.001 0 0 ;
  23. p2 = r1 0. r1 ;
  24. c1 = 0 0 r1 ;
  25. l1 = cerc p1 c1 p2 dini .05 dfin .15 ;
  26. *
  27. ** on maille un secteur de 2 degrees
  28. surf1 = rota l1 2. 1 (0 0 0) c1 ;
  29. elim surf1 .0015;
  30. surf1 = rege surf1 ;
  31. modl1 = mode surf1 mecanique coq4 dkt ;
  32. *
  33. rig1 = kp modl1 rg (0 0 -1) ;
  34. dep1 = manu chpo surf1 1 uz d1 ;
  35. f1 = resu (rig1 * dep1) ;
  36. fz1 = 180 * (maxi (exco fz f1));
  37. mess 'analytique :'(-1*pi) 'numerique :'fz1 ;
  38. err1 = ((fz1 + pi)/fz1)*100 ;
  39. mess ' erreur :'err1'%';
  40.  
  41. si ((abs err1) > .03) ;
  42. erre 5 ;
  43. sinon ;
  44. erre 0 ;
  45. finsi ;
  46.  
  47. fin ;
  48.  
  49.  
  50.  
  51.  
  52.  

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