Télécharger gdep2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gdep2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ******************************************************************
  6. * anneau sous pression non uniforme
  7. ******************************************************************
  8. *
  9. * position du probleme
  10. * --------------------
  11. * il s agit de determiner la deformee d un anneau
  12. * sous pression non uniforme suiveuse en grands
  13. * deplacements
  14. * reference:
  15. * ----------
  16. * Argyris J., Mlejnek H. P., Dynamics of structure, North Holland
  17. * 1991, p545
  18. ******************************************************************
  19.  
  20. option echo 1;
  21. option dime 2 elem seg2 ;
  22. opti mode 'PLAN' 'CONT' ;
  23.  
  24. graph = faux ;
  25.  
  26. *-------------------------------
  27. * definition de la geometrie
  28. *-------------------------------
  29.  
  30. r = 100 ;
  31. a = 0 0 ;
  32. b = 100 0 ;
  33. c = 0 100 ;
  34. lig1 = b cerc 18 a c ;
  35.  
  36. *--------------------------------------
  37. *definition des conditions aux limites
  38. *--------------------------------------
  39. cl1 = bloq rota (b et c) ;
  40. cl2 = bloq uy b ;
  41. cl3 = bloq ux c ;
  42. cl = cl1 et cl2 et cl3 ;
  43.  
  44. *---------------------------------------------------------
  45. *definition du modele, du materiau et des caracteristiques
  46. *---------------------------------------------------------
  47.  
  48. valyoun = 2.1e7 ;
  49. mo1 = modele lig1 mecanique elastique isotrope coq2 ;
  50. ma = mate mo1 youn valyoun nu 0. rho 7.8e3 ;
  51. ep1 = 1. ;
  52. ca = cara mo1 epai ep1 ;
  53. ma = ma et ca;
  54.  
  55. *-----------------------------
  56. * calcul de l objet chargement
  57. *-----------------------------
  58.  
  59. lf = prog 0 1000. ;
  60. lt = prog 0 1000. ;
  61. ls = prog 0 pas 0.025 1. pas 0.025 3. ;
  62.  
  63. p0bar = r * r * r * 12 / valyoun ;
  64. chx1 = lig1 coor 1 ;
  65. chp1 = chx1 ** 2 / r / r * 2. / p0bar ;
  66. chp1 = chp1 nomc pr ;
  67.  
  68. * Modele Chargement de pression
  69. MOPRE = 'MODE' LIG1 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  70. CHMZ = pres MOPRE (chan cham MOPRE chp1) ;
  71. *
  72. TTPS = 'TABL' ; TTPS. 0 = 0. ; TTPS. 1 = 1000. ;
  73. TPRE = 'TABL' ; TPRE. 0 = 0.*CHMZ ; TPRE. 1 = 1000.*CHMZ ;
  74. CHAA = 'CHAR' 'PRES' TTPS TPRE ;
  75. *
  76. *------------------------------
  77. * calcul non lineaire
  78. *------------------------------
  79.  
  80. tab2 = table ;
  81. tab2.grands_deplacements = vrai ;
  82. tab2.blocages_mecaniques = cl;
  83. tab2.modele = MO1 'ET' MOPRE ;
  84. tab2.caracteristiques = MA ;
  85. tab2.chargement = CHAA ;
  86. tab2.temps_calcules = ls;
  87. tab2.'PROCESSEURS' = 'MONO_PROCESSEUR';
  88. tab2.'PRECISION' = 1e-8;
  89. tab2.'K_SIGMA'= vrai;
  90. pasapas tab2 ;
  91.  
  92. *------------------------------
  93. * comparaison avec la reference
  94. *------------------------------
  95.  
  96. xxx = peche tab2 deplacements;
  97. uxb = (extr xxx ux b) ;
  98. uxbth = -86. ;
  99. err = (abs ((uxbth - uxb)/uxbth))*100 ;
  100. mess 'deplacement du point b, p0 = 3 ' ;
  101. mess 'reference 'uxbth 'castem' uxb '%erreur' err ;
  102.  
  103. si (err < 2) ;
  104. erre 0 ;
  105. sinon ;
  106. erre 5 ;
  107. finsi ;
  108.  
  109. si(graph) ;
  110. defo0 = defo lig1 xxx 0. VERT ;
  111. defo1 = defo lig1 xxx 1. ROUG ;
  112. trac (defo0 et defo1) ;
  113. finsi ;
  114.  
  115. fin ;
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  

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