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.0125 3. ;
  62.  
  63. pres1 = prog ;
  64. teta = 0. ;
  65. n1 = 19 ;
  66. xc = 90./18 ;
  67. p0bar = r * r * r * 12 / valyoun ;
  68. repeter bloc1 n1 ;
  69. pr1 = (1 + (cos (2*teta))) / p0bar ;
  70. pres1 = pres1 et (prog pr1) ;
  71. teta = teta + xc ;
  72. fin bloc1 ;
  73. *
  74. * Modele Chargement de pression
  75. MOPRE = 'MODE' LIG1 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  76. MAPRE = 'MATE' MOPRE 'PRES' 1 ;
  77. *
  78. * Description du chargement
  79. EVPRE = 'EVOL' 'MANU' 'ESTE' lt 'PRES' lf ;
  80. *
  81. CHPZ = 'MANU' 'CHPO' ('CHAN' 'POI1' ('EXTR' MOPRE 'MAIL')) 1
  82. 'PRES' pres1 ;
  83. CHMZ = 'CHAN' 'CHAM' CHPZ MOPRE ;
  84. CHMZ = 'CHAN' 'TYPE' CHMZ 'CONTRAINTES' ;
  85. *
  86. TTPS = 'TABL' ; TTPS. 0 = 0. ; TTPS. 1 = 1000. ;
  87. TPRE = 'TABL' ; TPRE. 0 = 0.*CHMZ ; TPRE. 1 = 1000.*CHMZ ;
  88. CHAA = 'CHAR' 'PRES' TTPS TPRE ;
  89. *
  90. *------------------------------
  91. * calcul non lineaire
  92. *------------------------------
  93.  
  94. tab2 = table ;
  95. tab2.grands_deplacements = vrai ;
  96. tab2.blocages_mecaniques = cl;
  97. tab2.modele = MO1 'ET' MOPRE ;
  98. tab2.caracteristiques = MA 'ET' MAPRE ;
  99. tab2.chargement = CHAA ;
  100. tab2.temps_calcules = ls;
  101. tab2.'PROCESSEURS' = 'MONO_PROCESSEUR';
  102. tab2.'PRECISION' = 1e-8;
  103. tab2.'K_SIGMA'= vrai;
  104. pasapas tab2 ;
  105.  
  106. *------------------------------
  107. * comparaison avec la reference
  108. *------------------------------
  109.  
  110. xxx = peche tab2 deplacements;
  111. uxb = (extr xxx ux b) ;
  112. uxbth = -86. ;
  113. err = (abs ((uxbth - uxb)/uxbth))*100 ;
  114. mess 'deplacement du point b, p0 = 3 ' ;
  115. mess 'reference 'uxbth 'castem' uxb '%erreur' err ;
  116.  
  117. si (err < 2) ;
  118. erre 0 ;
  119. sinon ;
  120. erre 5 ;
  121. finsi ;
  122.  
  123. si(graph) ;
  124. defo0 = defo lig1 xxx 0. VERT ;
  125. defo1 = defo lig1 xxx 1. ROUG ;
  126. trac (defo0 et defo1) ;
  127. finsi ;
  128.  
  129. fin ;
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  

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