Télécharger la plaquette Cast3M


fr en

Télécharger cou22.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  cou22.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 1 ;
  5. SAUT PAGE ;
  6. MESS '                                                       ';
  7. MESS '                                                       ';
  8. MESS '        TEST DE CISAILLEMENT SUR UN JOINT              ';
  9. MESS '                                                       ';
  10. MESS '                                                       ';
  11. MESS '               C_________________D                     ';
  12. MESS '               A                 B                     ';
  13. MESS '                                                       ';
  14. MESS '  JOINT FORME DES NOEUDS CDAB                          ';
  15. MESS '  CONDITIONS AUX LIMITES  :                            ';
  16. MESS '       - AB BLOQUE EN UX ET EN UY                      ';
  17. MESS '  CHARGEMENT              :                            ';
  18. MESS '       - FORCE DE CISAILLEMENT SUIVANT CD              ';
  19. MESS '                                                       ';
  20. MESS '  LE CHARGEMENT EST TEL QUE L ON RESTE DANS LE DOMAINE ';
  21. MESS '  ELASTIQUE                                            ';
  22. MESS '                                                       ';
  23. ********************************************************************
  24. *
  25. OPTION DIME 2 ;
  26. OPTION ELEM SEG2 MODE PLAN DEFO ;
  27. *
  28. COJANA = 1.5 ;
  29. DEPANA = 0.015 ;
  30. *
  31. * ----------DEFINITION DE LA GEOMETRIE DU JOINT ----------
  32. *
  33. A = 0.00 0.00 ;
  34. B = 1.00 0.00;
  35. C = 0.00 0.00 ;
  36. D = 1.00 0.00 ;
  37. *
  38. * ----------                  LIGNE                  ----------
  39. *
  40. AB =  A DROIT 1 B ;
  41. CD =  C DROIT 1 D ;
  42. *
  43. * ----------         MAILLAGE                ----------
  44. *
  45. OPTION ELEM RAC2 ;
  46. *
  47. JOINT = RACC 0.001 CD AB ;
  48. *
  49. * ----------     DEFINITION DES CONDITIONS AUX LIMITES   ----------
  50. *
  51. *
  52. *                        JOINT
  53. *
  54. CL1 = BLOQ AB UX ;
  55. CL2 = BLOQ AB UY ;
  56. CLJOINT =  CL1 ET CL2;
  57. *
  58. * ----------         DEFINITION DU MODELE DU JOINT      ----------
  59. *
  60. MODJOI = MODE JOINT MECANIQUE ELASTIQUE PLASTIQUE COULOMB JOI2 ;
  61. MAJOI  = MATE MODJOI KS   100.0  KN   100.0   EF  150.0   ECN  0.05
  62.                      COHE 2.0    FRIC  45.0   ;
  63. *
  64. * ----------            DEFINITION DU CHARGEMENT       ----------
  65. *
  66. CHPFO = FORCE (1.5 00.00) CD ;
  67. LIST CHPFO;
  68. *
  69. * ----------   DEFINITION DE LA PROGRESSION DU CHARGEMENT ----------
  70. *
  71.    LI1 = PROG 0. 2.0 ; LI2 = PROG 1.0 1.0 ;                                   
  72.    EV = EVOL MANU X LI1 Y LI2 ;     
  73.    CAR1 = CHAR 'MECA' CHPFO EV ;                                                                             
  74. *
  75. * ----------            RESOLUTION            ----------
  76. *
  77.    TAB = TABLE ;
  78.    LIS1 = PROG 0.0 1.0 2.0 ;
  79.    TAB.'BLOCAGES_MECANIQUES' = CLJOINT;
  80.    TAB.'CARACTERISTIQUES' = MAJOI;
  81.    TAB.'MODELE' = MODJOI;
  82.    TAB.'TEMPS_CALCULES' = LIS1;
  83.    TAB.'CHARGEMENT' = CAR1;
  84. TMASAU=table;
  85. tab . 'MES_SAUVEGARDES'=TMASAU;
  86. TMASAU .'DEFTO'=VRAI;
  87. TMASAU .'DEFIN'=VRAI;
  88.    PASAPAS TAB ;
  89. *
  90. * ---------- POST TRAITEMENT ----------
  91. *
  92. *
  93. * -----RECUPERATION DE LA TABLE DES DEPLACEMENTS
  94. *
  95. dtab1=tab.temps;
  96. *
  97. ndime= (dime dtab1) - 1;
  98. *
  99. * ----- EXTRACTION DES DIFFERENTES TABLES
  100. *
  101. *                                    DEPLACEMENTS 
  102. *
  103. i=0 ;
  104. repeter bloc1 ndime ;
  105.    i=i+1 ;
  106.    dep1=tab.deplacements.i ;
  107.    sig1=tab.contraintes.i ;
  108.    var1=tab.variables_internes.i ;
  109.    def1=tab.deformations_inelastiques.i ;
  110.    eps1 = epsi modjoi dep1 ;
  111. *
  112. * ----- REDUCTION DES CHPOINTS A L'ELEMENT JOINT 
  113. *
  114. SIGJOI = REDU SIG1 JOINT ;
  115. DEFJOI = REDU EPS1 JOINT ;
  116. DPLJOI = REDU DEP1 JOINT ;
  117. DEPJOI = REDU DEF1 JOINT ;
  118. VARJOI = REDU VAR1 JOINT ;
  119. *
  120. mess '              '   ;
  121. mess ' PAS DE TEMPS ' i ;
  122. mess ' ____________ '   ;
  123. mess '              '   ;
  124. mess ' CONTRAINTES  '   ;
  125. list sigjoi ;
  126. mess ' DEFO ELASTI  '   ;
  127. list defjoi ;
  128. mess ' DEFO INELASTI'   ;
  129. list depjoi ;
  130. mess ' VARIABLES INTERNES '   ;
  131. list varjoi ;
  132. mess ' DEPLACEMENTS '   ;
  133. list dpljoi ;
  134. *
  135. SIGJOI1 = EXTR SIGJOI SMSN 1 1 1 ;
  136. SIGJOI2 = EXTR SIGJOI SMSN 1 1 2 ;
  137. DEFJOI1 = EXTR DEFJOI DRSN 1 1 1 ;
  138. DEFJOI2 = EXTR DEFJOI DRSN 1 1 2 ;
  139. *
  140. RESI2 = ABS ( ( SIGJOI1 - COJANA)/COJANA ) ;
  141. RESI3 = ABS ( ( SIGJOI2 - COJANA)/COJANA ) ;
  142. RESI4 = ABS ( ( DEFJOI1 - DEPANA)/DEPANA ) ;
  143. RESI5 = ABS ( ( DEFJOI2 - DEPANA)/DEPANA ) ;
  144. *
  145. SI (  (RESI2 <EG 1E-4 ) ET (RESI3 <EG 1E-4 ) ET
  146.       (RESI4 <EG 1E-4 ) ET (RESI5 <EG 1E-4 ) ) ;
  147.    ERRE  0 ;
  148. SINON;
  149.    ERRE  5 ;
  150. FINSI ;
  151. *
  152. MESS '                                                       ' ;
  153. MESS '   SOLUTION THEORIQUE :                                ' ;
  154. MESS '     CONTRAINTES DANS LE JOINT' COJANA                   ;
  155. MESS '     DEFORMATION DANS LE JOINT' DEPANA                   ;
  156. MESS '                                                       ' ;
  157. MESS '   SOLUTION CALCULEE  :                                ' ;
  158. MESS '     CONTRAINTES DANS LE JOINT' SIGJOI1                  ;
  159. MESS '     DEFORMATION DANS LE JOINT' DEFJOI1                  ;
  160. fin bloc1 ;
  161. FIN ;
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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