Télécharger joi24.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi24.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 0 ;
  5. SAUT PAGE ;
  6. MESS'!==============================================================!';
  7. MESS'! !';
  8. MESS'! TEST JOI24 !';
  9. MESS'! !';
  10. MESS'! ESSAI DE TRACTION SUR UN JOINT 2D !';
  11. MESS'! !';
  12. MESS'! Un joint 2D JOI2 a son bord inferieur encastre. !';
  13. MESS'! Son bord superieur est libre. Un effort de traction est !';
  14. MESS'! exerce sur son bord superieur,perpendiculairement au plan !';
  15. MESS'! du joint. Le joint est incline de 20 degres par rapport !';
  16. MESS'! a l axe X. Les noeuds du joint sont A1, B1, IB1, IA1. !';
  17. MESS'! !';
  18. MESS'! Y !';
  19. MESS'! ! !';
  20. MESS'!B1,IB1!. !';
  21. MESS'! ! . !';
  22. MESS'! ! . !';
  23. MESS'! ! . <---- joint !';
  24. MESS'! ! . !';
  25. MESS'! !________________o__(_._ X !';
  26. MESS'! A1,IA1 !';
  27. MESS'! !';
  28. MESS'! angle o = 20 degres !';
  29. MESS'! longueur du joint = L = 2 !';
  30. MESS'! X(A1) = 1.879385242 !';
  31. MESS'! Y(B1) = 0.684040286 !';
  32. MESS'! !';
  33. MESS'! Force exercee sur le bord superieur = 100000 !';
  34. MESS'!> Deplacement de chacun des noeuds libres : Delta !';
  35. MESS'! Force = Raideur * Delta * Aire !';
  36. MESS'! ou !';
  37. MESS'! Raideur = raideur du joint = 4.2E10 !';
  38. MESS'! Delta = deplacement des noeuds du joint !';
  39. MESS'! Aire = aire du joint !';
  40. MESS'! = Longueur * largeur !';
  41. MESS'! L = 2 !';
  42. MESS'! l = 1 ( contraintes ou deformations planes ) !';
  43. MESS'! !';
  44. MESS'! Delta = 100000 / (4.2E10 * 2 * 1) !';
  45. MESS'! = 1.190476E-6 !';
  46. MESS'! !';
  47. MESS'! 1/ Deplacement suivant X : !';
  48. MESS'! !';
  49. MESS'! Deltax = delta * cos(70) !';
  50. MESS'! = 4.0716677E-7 !';
  51. MESS'! !';
  52. MESS'! !';
  53. MESS'! 2/ Deplacement suivant Y : !';
  54. MESS'! !';
  55. MESS'! Deltay = delta * sin(70) !';
  56. MESS'! = 1.1186815E-6 !';
  57. MESS'! !';
  58. MESS'!==============================================================!';
  59. OPTION ELEM SEG2 MODE PLAN CONT ;
  60. *
  61. SOLAY = 1.1186815E-6 ;
  62. SOLAX = 4.0716677E-7 ;
  63. *
  64. * ---------- DEFINITION DE LA SURFACE TOP DU JOINT ----------
  65. *
  66. A1 = 1.879385242 0.00 ;
  67. B1 = 0.00 0.684040286 ;
  68. *
  69. * ---------- MAILLAGE ----------
  70. *
  71. H1 = A1 DROIT 1 B1 ;
  72. *
  73. L1 = H1 ;
  74. *
  75. * ---------- DEFINITION DE LA SURFACE BOT DU JOINT ----------
  76. *
  77. IA1 = 1.879385242 0.00 ;
  78. IB1 = 0.00 0.684040286 ;
  79. *
  80. * ---------- MAILLAGE ----------
  81. *
  82. IH1 = IA1 DROIT 1 IB1 ;
  83. *
  84. IL1 = IH1 ;
  85. *
  86. * ---------- CREATION DU JOINT JOI2 ----------
  87. *
  88. OPTION ELEM RAC2 ;
  89. VOL = RACCORD 0.01 (INVE L1) (INVE IL1) ;
  90. *
  91. * ---------- DEFINITION DES CONDITIONS LIMITES ----------
  92. *
  93. CL11 = BLOQ IA1 UX ;
  94. CL12 = BLOQ IA1 UY ;
  95. CL1 = CL11 ET CL12 ;
  96. *
  97. CL21 = BLOQ IB1 UX ;
  98. CL22 = BLOQ IB1 UY ;
  99. CL2 = CL21 ET CL22 ;
  100. *
  101. CL = CL1 ET CL2 ;
  102. *
  103. * ---------- DEFINITION DU MODELE DU JOINT ----------
  104. *
  105. MOD1 = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' JOI2 ;
  106. MA1 = MATE MOD1 KS 4.2E08 KN 4.2E10 ;
  107. *
  108. * ---------- MATRICE DE RIGIDITE ----------
  109. *
  110. RI1 = RIGI MOD1 MA1 ;
  111. RI2 = RI1 ET CL ;
  112. *
  113. * ---------- FORCE DE TRACTION ----------
  114. *
  115. FO1 = FORCE ( 34202.01433 93969.26209 ) L1 ;
  116. *
  117. * ---------- RESOLUTION ----------
  118. *
  119. RE = RESO RI2 FO1 ;
  120. MESS ' ' ;
  121. MESS ' ' ;
  122. MESS ' Solution Analytique :' ;
  123. MESS ' ' ;
  124. MESS ' UX =' SOLAX ;
  125. MESS ' UY =' SOLAY ;
  126. MESS ' ' ;
  127. MESS ' ' ;
  128. MESS ' ' ;
  129. MESS ' Solution Calculee :' ;
  130. MESS ' ' ;
  131. LIST RE ;
  132.  
  133. *
  134. * ---------- CODE DE FONCTIONNEMENT ----------
  135. *
  136. DEPA1 = EXTR RE UX A1 ;
  137. DEPA2 = EXTR RE UY A1 ;
  138. RESIX = ABS( (DEPA1-SOLAX)/SOLAX ) ;
  139. RESIY = ABS( (DEPA2-SOLAY)/SOLAY ) ;
  140. *
  141. SI ( (RESIX &lt;EG 1E-4 ) ET (RESIY &lt;EG 1E-4 ) ) ;
  142. ERRE 0 ;
  143. SINO;
  144. ERRE 5 ;
  145. FINSI ;
  146.  
  147. *
  148. * ---------- CALCUL DES DEFORMATIONS ----------
  149. *
  150. EPS1 = EPSI MOD1 RE ;
  151. LIST EPS1 ;
  152. *
  153. * ---------- CALCUL DES CONTRAINTES ----------
  154. *
  155. SIG1 = SIGMA MOD1 MA1 RE ;
  156. LIST SIG1 ;
  157. *
  158. FIN ;
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  

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