Télécharger la plaquette Cast3M


fr en

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