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