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 DIME 2 ;
  60. OPTION ELEM SEG2 MODE PLAN CONT ;
  61. *
  62. SOLAY = 1.1186815E-6 ;
  63. SOLAX = 4.0716677E-7 ;
  64. *
  65. * ---------- DEFINITION DE LA SURFACE TOP DU JOINT ----------
  66. *
  67. A1 = 1.879385242 0.00 ;
  68. B1 = 0.00 0.684040286 ;
  69. *
  70. * ---------- MAILLAGE ----------
  71. *
  72. H1 = A1 DROIT 1 B1 ;
  73. *
  74. L1 = H1 ;
  75. *
  76. * ---------- DEFINITION DE LA SURFACE BOT DU JOINT ----------
  77. *
  78. IA1 = 1.879385242 0.00 ;
  79. IB1 = 0.00 0.684040286 ;
  80. *
  81. * ---------- MAILLAGE ----------
  82. *
  83. IH1 = IA1 DROIT 1 IB1 ;
  84. *
  85. IL1 = IH1 ;
  86. *
  87. * ---------- CREATION DU JOINT JOI2 ----------
  88. *
  89. OPTION ELEM RAC2 ;
  90. VOL = RACCORD 0.01 (INVE L1) (INVE IL1) ;
  91. *
  92. * ---------- DEFINITION DES CONDITIONS LIMITES ----------
  93. *
  94. CL11 = BLOQ IA1 UX ;
  95. CL12 = BLOQ IA1 UY ;
  96. CL1 = CL11 ET CL12 ;
  97. *
  98. CL21 = BLOQ IB1 UX ;
  99. CL22 = BLOQ IB1 UY ;
  100. CL2 = CL21 ET CL22 ;
  101. *
  102. CL = CL1 ET CL2 ;
  103. *
  104. * ---------- DEFINITION DU MODELE DU JOINT ----------
  105. *
  106. MOD1 = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' JOI2 ;
  107. MA1 = MATE MOD1 KS 4.2E08 KN 4.2E10 ;
  108. *
  109. * ---------- MATRICE DE RIGIDITE ----------
  110. *
  111. RI1 = RIGI MOD1 MA1 ;
  112. RI2 = RI1 ET CL ;
  113. *
  114. * ---------- FORCE DE TRACTION ----------
  115. *
  116. FO1 = FORCE ( 34202.01433 93969.26209 ) L1 ;
  117. *
  118. * ---------- RESOLUTION ----------
  119. *
  120. RE = RESO RI2 FO1 ;
  121. MESS ' ' ;
  122. MESS ' ' ;
  123. MESS ' Solution Analytique :' ;
  124. MESS ' ' ;
  125. MESS ' UX =' SOLAX ;
  126. MESS ' UY =' SOLAY ;
  127. MESS ' ' ;
  128. MESS ' ' ;
  129. MESS ' ' ;
  130. MESS ' Solution Calculee :' ;
  131. MESS ' ' ;
  132. LIST RE ;
  133.  
  134. *
  135. * ---------- CODE DE FONCTIONNEMENT ----------
  136. *
  137. DEPA1 = EXTR RE UX A1 ;
  138. DEPA2 = EXTR RE UY A1 ;
  139. RESIX = ABS( (DEPA1-SOLAX)/SOLAX ) ;
  140. RESIY = ABS( (DEPA2-SOLAY)/SOLAY ) ;
  141. *
  142. SI ( (RESIX &lt;EG 1E-4 ) ET (RESIY &lt;EG 1E-4 ) ) ;
  143. ERRE 0 ;
  144. SINO;
  145. ERRE 5 ;
  146. FINSI ;
  147.  
  148. *
  149. * ---------- CALCUL DES DEFORMATIONS ----------
  150. *
  151. EPS1 = EPSI MOD1 RE ;
  152. LIST EPS1 ;
  153. *
  154. * ---------- CALCUL DES CONTRAINTES ----------
  155. *
  156. SIG1 = SIGMA MOD1 MA1 RE ;
  157. LIST SIG1 ;
  158. *
  159. FIN ;
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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