Télécharger la plaquette Cast3M


fr en

Télécharger super2.dgibi

Retour à la liste

Numérotation des lignes : 

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

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