Télécharger la plaquette Cast3M


fr en

Télécharger soravs.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  soravs.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *********************************************************
  5. *
  6. *  Test du fonctionnement de la directive SORT 'AVS' et
  7. *  de l'opérateur LIRE 'AVS'. On effectue les deux 
  8. *  opérations et on vérifie que le CHPOINT lu est égal 
  9. *  à celui qui a été écrit.
  10. *
  11. *  Auteur : M. Bulik
  12. *
  13. *  Date : Septembre '95
  14. *
  15. *********************************************************
  16. *** options générales ...
  17.  
  18.     opti dime 3 echo 1 elem seg2 mode trid ;
  19.  
  20. *** noeuds ...
  21.  
  22.     dens 0.1 ;
  23.  
  24.     p1 = 0 0 0 ;
  25.     p2 = 1 0 0 ;
  26.  
  27. *** lignes ...
  28.  
  29.     li1 = p1 d 1 p2 ;
  30.     li2 = li1 plus (0.1 1. 0.) ;
  31.  
  32. *** surfaces ...
  33.  
  34.     opti elem qua4 ;
  35.  
  36.     su1 = li1 regl 1 li2 ;
  37.  
  38. *** volumes ...
  39.  
  40.     opti elem cub8 ;
  41.  
  42.     vo1 = su1 volu tran 2 (-0.1 0.1 1.) ;
  43.  
  44. *** Modèle ...
  45.  
  46.     mo1 = mode vo1 mecanique elastique ;
  47.  
  48. *** Matériau ...
  49.  
  50.     ma1 = mate mo1 YOUN 2.e+11 NU 0.3 ;
  51.  
  52. *** Rigidité et CL ...
  53.  
  54.     ri1 = rigi mo1 ma1 ;
  55.     cl1 = bloq depl su1 ;
  56.     ritot = ri1 et cl1 ;
  57.  
  58. *** Chargement ...
  59.  
  60.     fo1 = forc fx 10000 (vo1 poin proc (1 1 1)) ;
  61.  
  62. *** Résolution ...
  63.  
  64.     de1 = reso ritot fo1 ;
  65.  
  66. *** Les contraintes ...
  67.  
  68.     si1 = sigm mo1 ma1 de1 ;
  69.     si1g = chan GRAVITE si1 mo1 ;
  70.  
  71. *** la sortie AVS ...
  72.  
  73.     opti sort 'test-soravs.inp' ;
  74.  
  75.     valtemp = 1.23456789 ;
  76.     sort 'AVS' de1 si1g 'SUIT' 'TEMP' valtemp ;
  77.  
  78.     de2 = 2 * de1 ;
  79.     si2g = 2 * si1g ;
  80.     sort 'AVS' de2 si2g 'SUIT' ;
  81.  
  82. *** La lecture AVS ...
  83.  
  84.     opti sort 'toto' ;
  85.     opti lect 'test-soravs.inp' ;
  86.  
  87.     tblect1 = lire AVS ;
  88.  
  89.     tblect2 = lire AVS ;
  90.  
  91. *** La comparaison ...
  92.  
  93.     list (nbno vo1) ;
  94.     elim (vo1 et (tblect1 . LEMAILLA) et (tblect2 . LEMAILLA)) ;
  95.  
  96.     de1redu = redu de1 vo1 ;
  97.     chpodiff = de1redu - (tblect1 . LECHPOIN) ;
  98.     list (abs (maxi chpodiff)) ;
  99.     si ((abs (maxi chpodiff)) > 1.e-14) ;
  100.        erreur 5 ;
  101.     finsi ;
  102.  
  103. *    chamdiff = si1g + (-1 * (tblect . LEMCHAML)) ;
  104. *    list (abs (maxi chamdiff)) ;
  105. *    si ((abs (maxi chamdiff)) > 1.e-10) ;
  106. *       erreur 5 ;
  107. *    finsi ;
  108.  
  109.     valdiff = abs (valtemp - (tblect1 . 'time'));
  110.     list valdiff ;
  111.     si(valdiff > 1.e-7) ;
  112.        erreur 5 ;
  113.     finsi ;
  114.  
  115.     de2redu = redu de2 vo1 ;
  116.     chpodiff = de2redu - (tblect2 . LECHPOIN) ;
  117.     list (abs (maxi chpodiff)) ;
  118.     si ((abs (maxi chpodiff)) > 1.e-14) ;
  119.        erreur 5 ;
  120.     finsi ;
  121.  
  122. *** Bye ...
  123.  
  124.     fin ;
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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