Télécharger la plaquette Cast3M


fr en

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