Télécharger la plaquette Cast3M


fr en

Télécharger super1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  super1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. opti echo 0;
  6. *
  7. * definition of the mesh points and of the element
  8. *
  9. p1=1 0; p2=2 0; p3=3 0;
  10. mes1=manu 'SUPE' p1 p2 p3;
  11. *
  12. *                               |2 1 0|
  13. * definition of the K matrix  K=|1 2 1|
  14. *                               |0 1 2|
  15. *
  16. lis1=prog 2 1 0 1 2 1 0 1 2;
  17. rig1=manu 'RIGIDITE' 'TYPE' 'RIGIDITE' mes1 (mots 'UY') lis1;
  18. *
  19. * use of the 'MASSE' option
  20. *
  21. masbl1=mass 'UY' 1 p3;
  22. super1=super 'RIGIDITE' rig1 masbl1;
  23. bb2 mas2=super 'MASSE' super1 masbl1 'LCHP';
  24. u2=extr bb2 1;
  25. *
  26. *mess 'Solution with the "MASSE" option';
  27. *mess (extr u2 'UY' p1) (extr u2 'UY' p2) (extr u2 'UY' p3);
  28. *
  29. * solution of the dirichlet problem
  30. *
  31. bl1=bloq p3 'UY';
  32. f1=depi 1 bl1;
  33. u1=reso (rig1 et bl1) f1;
  34. u1=redu u1 (extr u2 'MAIL');
  35. *
  36. * error check
  37. *
  38. er1=u2-u1;
  39. si ((xtx er1) > 1.D-15); erre 5;
  40. sinon;                 ; erre 0; finsi;
  41. FIN;
  42.  
  43.  
  44.  
  45.  
  46.  

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