Télécharger la plaquette Cast3M


fr en

Télécharger simpl2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  simpl2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph='N';
  5. saut page;
  6. mess 'test de SIMPLE : resistance limite d"un treillis de 3 barres';
  7. opti dime 2 elem seg2;
  8. *****************************************************************
  9. * test du simplex sur un treillis de 3 barres (avec utilisation 
  10. * de la procedure ANLIMTRE)
  11. *
  12. * on considere 3 barres pouvant suportees une contrainte limite
  13. * egale a condlim. L'arrangement des trois barres est le suivant:
  14. *
  15. *            p3                            1           2
  16. *             +                   
  17. *            ***                           ^
  18. *           * * *                          |
  19. *          *  *  *                         |
  20. *         *   *   *                        |         ---->
  21. *        *    *    *
  22. *       +     +     +          
  23. *      p1    p4    p2                direction de sollicitation           
  24. *
  25. * on cherche l'intensite de la charge limite suportee par ce
  26. * treillis, pour les 2 directions de sollicitation indiquees
  27. * dessus. Le resultat attendu est evidemment:
  28. *
  29. * cas 1: intensite=(1+2*sqrt(2))*conlim
  30. * cas 2: intensite=   2*sqrt(2) *conlim
  31. *
  32. * lire les commentaires associes a ANLIMTRE
  33. *
  34. * PP 1/9/92
  35. *****************************************************************
  36. conlim=1.;
  37. reacmax=5.;
  38. *
  39. * 1) preparation du maillage
  40. *
  41. p1=0 0; p2=2 0; p3=1 1; p4=1 0;
  42. mesh=(p1 d 1 p3) et (p3 d 1 p2) et (p3 d 1 p4);
  43. *
  44. * WARNING: on tasse pour avoir les noeuds du maillage
  45. *          avec une numerotation ABSOLUE correcte
  46. tass mesh;
  47. *
  48. * 2) point(s) bloque(s) (du maillage)
  49. *
  50. pbloq=p1 et p2 et p4;
  51. si (ega (type pbloq) 'POINT   ');
  52.   pbloq=pbloq et pbloq;
  53.   pbloq=pbloq elem 1;
  54. finsi;
  55. *
  56. * 3) point/composante sollicite (du maillage)
  57. *
  58. psoll=p3; vsoll1=0. 1.; vsoll2=1. 0.;
  59. *
  60. * 4) preparation de la table d'entree de ANLIMTRE
  61. *
  62. tt=table;
  63. tt.'MESH'=mesh;
  64. tt.'PBLOQ'=pbloq;
  65. tt.'PSOLL'=psoll;
  66. tt.'CONLIM'=conlim;
  67. tt.'REACMAX'=reacmax;
  68. *
  69. * 5) premier cas - impression et erreur -
  70. *
  71. tt.'VSOLL'=vsoll1;
  72. iok=ANLIMTRE tt;
  73. si (iok ega 0);
  74.   solut=tt.'CHARLIM';
  75.   theorie=conlim * (1 + (2. ** .5));
  76.   mess 'premiere direction de sollicitation';
  77.   mess '--->theorie=' theorie ' calcul=' solut;
  78.   si (abs (theorie - solut) > 1.d-6); erre 5;
  79.   sinon;                               erre 0;  finsi;
  80. finsi;
  81. *
  82. * 6) deuxieme cas - impression et erreur -
  83. *
  84. tt.'VSOLL'=vsoll2;
  85. iok=ANLIMTRE tt;
  86. si (iok ega 0);
  87.   solut=tt.'CHARLIM';        
  88.   theorie=conlim * (2. ** .5);
  89.   mess 'deuxieme direction de sollicitation';
  90.   mess '--->theorie=' theorie ' calcul=' solut;
  91.   si (abs (theorie - solut) > 1.d-6); erre 5;
  92.   sinon;                               erre 0;  finsi;
  93. finsi;
  94. *
  95. fin;
  96.  
  97.  
  98.  
  99.  

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