Télécharger joi43.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi43.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 0 ;
  5. SAUT PAGE ;
  6. *===================================================================*
  7. * *
  8. * TEST JOI43 *
  9. * *
  10. * ESSAI DE CISAILLEMENT SUR UN JOINT 3D *
  11. * *
  12. * Un joint 3D JOI4 a sa surface inferieure encastree. Sa *
  13. * surface superieure est libre. Un effort de cisaillement est *
  14. * exercee sur sa surface superieure. *
  15. * *
  16. * Force totale exercee sur la surface superieure = 100000 *
  17. * ===> Force repartie sur les quatre noeuds de la surface = *
  18. * 100000 / 4 = 25000 *
  19. * *
  20. * Solution analytique : *
  21. * *
  22. * En ecrivant le principe des travaux virtuels, on montre *
  23. * que le deplacement delta en chaque noeud est egal, dans le *
  24. * cas d'un essai de traction uniaxial, a : *
  25. * F *
  26. * delta = ---------------- *
  27. * k * L * l *
  28. * ou *
  29. * F = force totale exercee sur la structure *
  30. * k = raideur *
  31. * L = longueur de l'element joint *
  32. * l = largeur de l'element joint *
  33. * *
  34. * Application numerique : *
  35. * *
  36. * F = 100000.0 *
  37. * k = 4.2E10 *
  38. * L = 2.0 *
  39. * l = 2.0 *
  40. * *
  41. * ===> Delta = 100000 / (4.2E10 * 2.0 * 2.0) *
  42. * = 5.95238E-7 *
  43. * *
  44. *===================================================================*
  45. OPTION DIME 3 ;
  46. OPTION ELEM QUA4 MODE TRID ;
  47. *
  48. SOLANA = 5.95238E-7 ;
  49. SOLNUL = 0.0 ;
  50. *
  51. * ---------- DEFINITION DE LA SURFACE TOP DU JOINT ----------
  52. *
  53. A1 = 0.00 0.00 0.00 ;
  54. B1 = 2.00 0.00 0.00 ;
  55. C1 = 2.00 2.00 0.00 ;
  56. D1 = 0.00 2.00 0.00 ;
  57. *
  58. * ---------- MAILLAGE ----------
  59. *
  60. O1 = 0.00 0.00 10.00 ;
  61. H1 = A1 DROIT 1 B1 ;
  62. H2 = B1 DROIT 1 C1 ;
  63. H3 = C1 DROIT 1 D1 ;
  64. H4 = D1 DROIT 1 A1 ;
  65. *
  66. L1 = H1 ET H2 ET H3 ET H4 ;
  67. S1 = SURF L1 PLAN ;
  68. *
  69. * ---------- DEFINITION DE LA SURFACE BOT DU JOINT ----------
  70. *
  71. IA1 = 0.00 0.00 0.00 ;
  72. IB1 = 2.00 0.00 0.00 ;
  73. IC1 = 2.00 2.00 0.00 ;
  74. ID1 = 0.00 2.00 0.00 ;
  75. *
  76. * ---------- MAILLAGE ----------
  77. *
  78. IH1 = IA1 DROIT 1 IB1 ;
  79. IH2 = IB1 DROIT 1 IC1 ;
  80. IH3 = IC1 DROIT 1 ID1 ;
  81. IH4 = ID1 DROIT 1 IA1 ;
  82. *
  83. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  84. IS1 = SURF IL1 PLAN ;
  85. *
  86. * ---------- CREATION DU JOINT JOI4 ----------
  87. *
  88. OPTION ELEM LIA4 ;
  89. VOL = LIAISON 0.00001 S1 IS1 ;
  90. *
  91. * ---------- DEFINITION DES CONDITIONS LIMITES ----------
  92. *
  93. CL11 = BLOQ IA1 UX ;
  94. CL12 = BLOQ IA1 UY ;
  95. CL13 = BLOQ IA1 UZ ;
  96. CL1 = CL11 ET CL12 ET CL13 ;
  97. *
  98. CL21 = BLOQ IB1 UX ;
  99. CL22 = BLOQ IB1 UY ;
  100. CL23 = BLOQ IB1 UZ ;
  101. CL2 = CL21 ET CL22 ET CL23 ;
  102. *
  103. CL31 = BLOQ IC1 UX ;
  104. CL32 = BLOQ IC1 UY ;
  105. CL33 = BLOQ IC1 UZ ;
  106. CL3 = CL31 ET CL32 ET CL33 ;
  107. *
  108. CL41 = BLOQ ID1 UX ;
  109. CL42 = BLOQ ID1 UY ;
  110. CL43 = BLOQ ID1 UZ ;
  111. CL4 = CL41 ET CL42 ET CL43 ;
  112. *
  113. CL = CL1 ET CL2 ET CL3 ET CL4 ;
  114. *
  115. * ---------- DEFINITION DU MODELE DU JOINT ----------
  116. *
  117. MOD1 = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' JOI4 ;
  118. MA1 = MATE MOD1 KS 4.2E10 KN 4.2E10 ;
  119. *
  120. * ---------- MATRICE DE RIGIDITE ----------
  121. *
  122. RI1 = RIGI MOD1 MA1 ;
  123. RI2 = RI1 ET CL ;
  124. *
  125. * ---------- FORCE DE TRACTION ----------
  126. *
  127. FO1 = FORCE ( 0.0 100000.0 0.0 ) S1 ;
  128. *
  129. * ---------- RESOLUTION ----------
  130. *
  131. RE = RESO RI2 FO1 ;
  132. MESS ' ' ;
  133. MESS ' ___________________________________________________________ ' ;
  134. MESS '! !' ;
  135. MESS '! JOINT 3D JOI4 !' ;
  136. MESS '! !' ;
  137. MESS '! - Surface inferieure du joint encastree. !' ;
  138. MESS '! - Surface superieure du joint libre . !' ;
  139. MESS '! - Force de traction exercee sur la surface superieure !' ;
  140. MESS '! dans le plan du joint dans la direction y !' ;
  141. MESS '! !' ;
  142. MESS '! !' ;
  143. MESS '! Solution analytique : !' ;
  144. MESS '! !' ;
  145. MESS '! En ecrivant le principe des travaux virtuels, on !' ;
  146. MESS '!montre, sous certaines hypotheses, que le deplacement delta!' ;
  147. MESS '!en chaque noeud, dans la direction de l essai uniaxial, est!' ;
  148. MESS '!egal a : !' ;
  149. MESS '! F !' ;
  150. MESS '! delta = ---------------- !' ;
  151. MESS '! k * L * l !' ;
  152. MESS '! ou !' ;
  153. MESS '! F = force totale exercee sur la structure !' ;
  154. MESS '! k = raideur !' ;
  155. MESS '! L = longueur de l element joint !' ;
  156. MESS '! l = largeur de l element joint !' ;
  157. MESS '! !' ;
  158. MESS '! Application numerique : !' ;
  159. MESS '! !' ;
  160. MESS '! F = 100000.0 !' ;
  161. MESS '! k = 4.2E10 !' ;
  162. MESS '! L = 2.0 !' ;
  163. MESS '! l = 2.0 !' ;
  164. MESS '! !' ;
  165. MESS '! delta = 100000 / (4.2E10 * 2.0 * 2.0) !' ;
  166. MESS '! = 5.95238E-7 !' ;
  167. MESS '!___________________________________________________________!' ;
  168. MESS ' ' ;
  169. MESS ' ' ;
  170. MESS ' Solution Analytique :' ;
  171. MESS ' ' ;
  172. MESS ' UX =' SOLNUL ;
  173. MESS ' UY =' SOLANA ;
  174. MESS ' UZ =' SOLNUL ;
  175. MESS ' ' ;
  176. MESS ' ' ;
  177. MESS ' ' ;
  178. MESS ' Solution Calculee :' ;
  179. MESS ' ' ;
  180. LIST RE ;
  181. *
  182. * ---------- CODE DE FONCTIONNEMENT ----------
  183. *
  184. DEPA1 = EXTR RE UY A1 ;
  185. RESI = ABS( (DEPA1-SOLANA)/SOLANA ) ;
  186. *
  187. SI (RESI <EG 1E-4 ) ;
  188. ERRE 0 ;
  189. SINO;
  190. ERRE 5 ;
  191. FINSI ;
  192. *
  193. * ---------- CALCUL DES DEFORMATIONS ----------
  194. *
  195. EPS1 = EPSI MOD1 RE ;
  196. LIST EPS1 ;
  197. *
  198. * ---------- CALCUL DES CONTRAINTES ----------
  199. *
  200. SIG1 = SIGMA MOD1 MA1 RE ;
  201. LIST SIG1 ;
  202. *
  203. * ---------- CALCUL DE BSIGMA ----------
  204. *
  205. BSI1 = BSIGMA MOD1 SIG1 ;
  206. LIST BSI1 ;
  207. *
  208. FIN ;
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  

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