Télécharger joi22.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi22.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Joi22.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTI ECHO 0 ;
  8. SAUT PAGE ;
  9. MESS'!==============================================!';
  10. MESS'! !';
  11. MESS'! TEST JOI22 !';
  12. MESS'! !';
  13. MESS'! ESSAI DE TRACTION SUR UN JOINT 2D !';
  14. MESS'! !';
  15. MESS'! Un joint 2D JOI2 a son bord inferieur !';
  16. MESS'! encastre. !';
  17. MESS'! Son bord superieur est libre. Un effort de !';
  18. MESS'! traction est exercé sur son bord superieur. !';
  19. MESS'! !';
  20. MESS'! !';
  21. MESS'! Solution analytique : !';
  22. MESS'! !';
  23. MESS'! En ecrivant le principe des travaux !';
  24. MESS'! virtuels,on montre que le deplacement delta !';
  25. MESS'! en chaque noeud est egal, dans le cas d un !';
  26. MESS'! essai de traction uniaxial, a : !';
  27. MESS'! F !';
  28. MESS'! delta = ---------------- !';
  29. MESS'! k * L * l !';
  30. MESS'! ou !';
  31. MESS'! F = force totale exercee sur la structure !';
  32. MESS'! k = raideur !';
  33. MESS'! L = longueur de l element joint !';
  34. MESS'! l = largeur de l element joint !';
  35. MESS'! !';
  36. MESS'! Application numerique : !';
  37. MESS'! !';
  38. MESS'! F = 100000.0 !';
  39. MESS'! k = 4.2E10 !';
  40. MESS'! L = 2.0 !';
  41. MESS'! l = 1.0 !';
  42. MESS'! !';
  43. MESS'! Delta = 100000 / (4.2E10 * 2.0 * 1.0) !';
  44. MESS'! = 1.19047619E-6 !';
  45. MESS'! !';
  46. MESS'!==============================================!';
  47. OPTION ELEM SEG2 MODE PLAN CONT ;
  48. *
  49. SOLANA = 1.190476E-6 ;
  50. SOLNUL = 0.0 ;
  51. *
  52. * ------- DEFINITION DE LA SURFACE TOP DU JOINT -------
  53. *
  54. A1 = 0.00 0.00 ;
  55. B1 = 2.00 0.00 ;
  56. *
  57. * ---------- MAILLAGE ----------
  58. *
  59. H1 = A1 DROIT 1 B1 ;
  60. *
  61. L1 = H1 ;
  62. *
  63. * ------ DEFINITION DE LA SURFACE BOT DU JOINT --------
  64. *
  65. IA1 = 0.00 0.00 ;
  66. IB1 = 2.00 0.00 ;
  67. *
  68. * ---------- MAILLAGE ----------
  69. *
  70. IH1 = IA1 DROIT 1 IB1 ;
  71. *
  72. IL1 = IH1 ;
  73. *
  74. * ---------- CREATION DU JOINT JOI2 ----------
  75. *
  76. OPTION ELEM RAC2 ;
  77. VOL = RACCORD 0.00001 L1 IL1 ;
  78. *
  79. * --------- DEFINITION DES CONDITIONS LIMITES ---------
  80. *
  81. CL11 = BLOQ IA1 UX ;
  82. CL12 = BLOQ IA1 UY ;
  83. CL1 = CL11 ET CL12 ;
  84. *
  85. CL21 = BLOQ IB1 UX ;
  86. CL22 = BLOQ IB1 UY ;
  87. CL2 = CL21 ET CL22 ;
  88. *
  89. CL = CL1 ET CL2 ;
  90. *
  91. * ----- DEFINITION DU MODELE DU JOINT ---------
  92. *
  93. MOD1 =MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' JOI2;
  94. MA1 = MATE MOD1 KS 4.2E08 KN 4.2E10 ;
  95. *
  96. * ------------- MATRICE DE HOOK ------------
  97. *
  98. HOO1 = HOOKE MOD1 MA1 ;
  99. ******LIST HOO1 ;
  100. *
  101. * ---------- MATRICE DE RIGIDITE ---------
  102. *
  103. ******RI1 = RIGI MOD1 MA1 ;
  104. RI1 = RIGI MOD1 HOO1 ;
  105. RI2 = RI1 ET CL ;
  106. *
  107. * ---------- FORCE DE TRACTION ---------
  108. *
  109. FO1 = FORCE ( 0. 100000.0 ) L1 ;
  110. *
  111. * --------- RESOLUTION ----------
  112. *
  113. RE = RESO RI2 FO1 ;
  114. MESS ' ' ;
  115. MESS ' ' ;
  116. MESS ' Solution Analytique :' ;
  117. MESS ' ' ;
  118. MESS ' UX =' SOLNUL ;
  119. MESS ' UY =' SOLANA ;
  120. MESS ' ' ;
  121. MESS ' ' ;
  122. MESS ' ' ;
  123. MESS ' Solution Calculee :' ;
  124. MESS ' ' ;
  125. LIST RE ;
  126. *
  127. * ---------- CODE DE FONCTIONNEMENT ----------
  128. *
  129. DEPA1 = EXTR RE UY A1 ;
  130. RESI = ABS( (DEPA1-SOLANA)/SOLANA ) ;
  131. *
  132. SI (RESI <EG 1E-4 ) ;
  133. ERRE 0 ;
  134. SINO;
  135. ERRE 5 ;
  136. FINSI ;
  137. *
  138. * ---------- CALCUL DES DEFORMATIONS ----------
  139. *
  140. EPS1 = EPSI MOD1 RE ;
  141. LIST EPS1 ;
  142. *
  143. * ---------- CALCUL DES CONTRAINTES ----------
  144. *
  145. SIG1 = SIGMA MOD1 MA1 RE ;
  146. LIST SIG1 ;
  147. *
  148. FIN ;
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  

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