Télécharger la plaquette Cast3M


fr en

Télécharger joi25.dgibi

Retour à la liste

Numérotation des lignes : 

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

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