Datasets examples

* fichier :  elas18.dgibi
saut page;
********************************************************************
* Etude d'une poutre encastree a une extremite et chargee a l'autre
*    on néglige l'énergie associée à l'effort tranchant
*    La flexion se fait dans le plan xoy
*
*  encastrement                       extremite chargee
*
*     |-------------------------------------*
*                                           ^
*     p1                                    p2
********************************************************************
* Notations:
*
* h =hauteur
* b = largeur
* L =longueur totale
* I = b * h**3 / 12
* P =charge appliquee en P2
* E =module d'Young
* nu=Poisson
********************************************************************
* Objectif du test:
* on cherche à verifier la solution theorique en flexion pure
* donnant la flèche en bout
********************************************************************
* 
* 1) option de calcul
opti dime 3 elem seg2;
*
* 2) quelques donnees et calcul de quelques constantes
*     pour la solution theorique
*
h=1.; E=1.e10; nu=.3; L=10.; b = 0.5 ;
P = 50000. ;
h3 = h ** 3 ;
L3 = l ** 3 ;
inerz = ( b * h3 ) /12. ;
inery = ( (b ** 3 ) * h ) /12. ;
*
*
* 3) maillage
*
nel=5;
p1 = 0 0 0; p2 = L 0 0;
d1= p1 d nel p2;
*
* 4) modele et materiaux
*    et  caracteristiques geometriques
*
modti=modl d1 mecanique elastique pout;
mater=matr modti 'YOUN' E 'NU' nu;
carac=carb modti 'SECT' ( b * h )  'INRY' inery 'INRZ' inerz
          'TORS' (inery + inerz)   VECT ( 0 1 0 );
mc=mater et carac;
*
* 5) chargement
*
chrti = force 'FY' P  p2 ;
*
* 6) bloquage
*
bloqtot=bloq depl rota p1;
*
*
* 7) calcul de la rigidite
*
  rigti=rigi modti mc;
  rigtti=rigti et bloqtot ;
*
* 8) solution
*
  depti=reso rigtti chrti;
*
* 9) calcul de la solution analytique
*
  theorie=  ( P * L3 ) / ( ( 3. * E ) * inerz ) ;
  calcul= extr depti 'UY' p2 ;
  erreur=abs ((calcul - theorie) / theorie);
*
* 10) comparaison et message
*
  mess 'calcul,theorie,erreur:  ' calcul theorie erreur;
  mess '***********************';
  si (erreur > 5.d-4);
       erre 5;
  sinon;
       erre 0;
  finsi;
fin;