Télécharger gdep2_boucle.dgibi

Retour à la liste

Numérotation des lignes :

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

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