Télécharger la plaquette Cast3M


fr en

Télécharger cou31.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  cou31.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 0;                                                                    
  5. SAUT PAGE;   
  6. MESS '                                                       ';
  7. MESS '                                                       ';
  8. MESS '                                                       ';
  9. MESS '  TEST AVEC 2 CUBES ET UN JOINT (3D)                   ';
  10. MESS '                                                       ';
  11. MESS '  LE CUBE DU HAUT EST SOUMIS A UN DEPLACEMENT IMPOSE   ';
  12. MESS '  DE 0.1 DANS LA DIRECTION DE CISAILLEMENT X.          ';
  13. MESS '  LE CUBE DU BAS EST BLOQUE EN X, EN Y ET EN Z.        ';
  14. MESS '  LE JOINT PLASTIFIE EN CISAILLEMENT                   ';
  15. MESS '  LA CONTRAINTE DE CISAILLEMENT DANS LA                ';
  16. MESS '  DIRECTION X ATTEINT LA COHESION, ALORS QUE LA        ';
  17. MESS '  CONTRAINTE DE CISAILLEMENT DANS LA DIRECTION Y       ';
  18. MESS '  ET LA CONTRAINTE NORMALE VALENT 0.                   ';
  19. MESS '  LE DEPLACEMENT DU JOINT DANS LA DIRECTION X EST DE   ';
  20. MESS '  0.1                                                  ';
  21. MESS '                                                       ';
  22. MESS '                                                       ';
  23. MESS '                                                       ';
  24. MESS '               B _______________ C                     ';
  25. MESS '                |               |                      ';
  26. MESS '                |               |                      ';
  27. MESS '                |    CUBE 1     |                      ';
  28. MESS '                |               |                      ';
  29. MESS '                |               |                      ';
  30. MESS '                |_______________|                      ';
  31. MESS '               A|    JOINT      |H                     ';
  32. MESS '               G|_______________|D                     ';
  33. MESS '                |               |                      ';
  34. MESS '                |               |                      ';
  35. MESS '                |    CUBE 2     |                      ';
  36. MESS '                |               |                      ';
  37. MESS '   Z            |_______________|                      ';
  38. MESS '   |   Y       K                 F                     ';
  39. MESS '   |  /                                                ';
  40. MESS '   | /                                                 ';
  41. MESS '   |/_______ X                                         ';
  42. ********************************************************************                                                                   
  43. *
  44. OPTION DIME 3 ;
  45. OPTION ELEM QUA4 MODE TRID ;
  46. *
  47. OEIL = -1.E1 -1.E6 1.E1 ;
  48. *
  49. * ----------DEFINITION DE LA GEOMETRIE DU CUBE SUPERIEUR----------
  50. *
  51. SA1 = 0.00 0.00 0.00 ;
  52. SB1 = 1.00 0.00 0.00 ;
  53. SC1 = 1.00 1.00 0.00 ;
  54. SD1 = 0.00 1.00 0.00 ;
  55. *
  56. * ----------                  MAILLAGE                  ----------
  57. *
  58. SH1 =  SA1 DROIT 1 SB1 ;
  59. SH2 =  SB1 DROIT 1 SC1 ;
  60. SH3 =  SC1 DROIT 1 SD1 ;
  61. SH4 =  SD1 DROIT 1 SA1 ;
  62. *
  63. SL1 = SH1 ET SH2 ET SH3 ET SH4 ;
  64. SS1 = SURF SL1 PLAN ;
  65. *
  66. * ----------DEFINITION DE LA GEOMETRIE DU CUBE INFERIEUR----------
  67. *
  68. IA1 = 0.00 0.00 0.00 ;
  69. IB1 = 1.00 0.00 0.00 ;
  70. IC1 = 1.00 1.00 0.00 ;
  71. ID1 = 0.00 1.00 0.00 ;
  72. *
  73. * ----------                  MAILLAGE                  ----------
  74. *
  75. IH1 =  IA1 DROIT 1 IB1 ;
  76. IH2 =  IB1 DROIT 1 IC1 ;
  77. IH3 =  IC1 DROIT 1 ID1 ;
  78. IH4 =  ID1 DROIT 1 IA1 ;
  79. *
  80. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  81. IS1 = SURF IL1 PLAN ;
  82. *
  83. OPTION ELEM CUB8 ;
  84. *
  85. VOLS = SS1 VOLU 1 TRAN ( 0.00 0.00   1.00 ) ;
  86. VOLI = IS1 VOLU 1 TRAN ( 0.00 0.00  -1.00 ) ;
  87. FACS = FACE 2 VOLS ;
  88. FACI = FACE 2 VOLI ;
  89. VT1  = VOLS ET VOLI ;
  90. *   definition pour la verif finale
  91. P3 = NOEUD 3 ;
  92. *
  93. * ----------DEFINITION DE LA GEOMETRIE DU JOINT----------
  94. *
  95. * JOINT FORME DES FACES SS1 ET IS1
  96. *
  97. * ----------         MAILLAGE                ----------
  98. *
  99. OPTION ELEM LIA4 ;
  100. *
  101. JOINT = LIAISON 0.001 SS1 IS1 ;
  102. *
  103. *
  104. * ----------      MAILLAGE TOTAL               ----------
  105. *
  106. GEOTOT = VT1 ET JOINT  ;
  107. *
  108. * ---------- DEFINITION DES CONDITIONS AUX LIMITES ----------
  109. *
  110. CL11 = BLOQ VOLI UX ;
  111. CL12 = BLOQ VOLI UY ;
  112. CL13 = BLOQ VOLI UZ ;
  113. CL1  = CL11 ET CL12 ET CL13 ;
  114. *
  115. CL21 = BLOQ FACI UX ;
  116. CL22 = BLOQ FACI UY ;
  117. CL23 = BLOQ FACI UZ ;
  118. CL2  = CL21 ET CL22 ET CL23 ;
  119. *
  120. CL3  = BLOQ SS1 UX  ;
  121. DEP  = DEPI CL3 0.1 ;
  122. *
  123. CL   = CL1 ET CL3 ;
  124. *
  125. * ----------        DEFINITION DU MODELE DES CUBES       ----------
  126. *
  127. MODCUB = MODE VT1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  128. MACUB  = MATE MODCUB YOUN 100.0 NU 0.0000001  ;
  129. *
  130. * ----------         DEFINITION DU MODELE DU JOINT      ----------
  131. *
  132. MODJOI = MODE JOINT MECANIQUE ELASTIQUE PLASTIQUE COULOMB JOI4 ;
  133. MAJOI  = MATE MODJOI KS   100.0  KN   100.0   EF  150.0   ECN  0.05
  134.                      COHE 2.0     FRIC 45.0    ;
  135. *
  136. * ---------- MODELE ET MATERIAU DE TOUTE LA STRUCTURE ----------
  137. *
  138. MODTOT = MODCUB ET MODJOI;
  139. MATTOT = MACUB  ET MAJOI ;
  140. *
  141. * ----------        DEFINITION DU CHARGEMENT         ----------
  142. *                                                                             
  143.    LI1 = PROG 0.0  1.0  2.0 ;
  144.    LI2 = PROG 0.0  1.0  1.0 ;                                   
  145.    EV = EVOL MANU X LI1 Y LI2 ;     
  146. *
  147.    CAR1 = CHAR 'DIMP' DEP EV ; 
  148. *
  149. * ----------            RESOLUTION            ----------
  150. *
  151.    TAB = TABLE ;
  152.    LIS1 = PROG 0.0 PAS 1.0 2.0 ; 
  153.    TAB.'BLOCAGES_MECANIQUES' = CL;
  154.    TAB.'CARACTERISTIQUES' = MATTOT;
  155.    TAB.'CHARGEMENT' = CAR1;
  156.    TAB.'TEMPS_CALCULES' = LIS1;
  157.    TAB.'MODELE' = MODTOT;
  158.  
  159.    TMASAU=table;
  160. tab . 'MES_SAUVEGARDES'=TMASAU;
  161. TMASAU .'DEFTO'=VRAI;
  162. TMASAU .'DEFIN'=VRAI;
  163. PASAPAS TAB ;
  164. *
  165. * ---------- POST TRAITEMENT ----------
  166. *
  167. *
  168. * -----RECUPERATION DE LA TABLE DES DEPLACEMENTS
  169. *
  170. dtab1= tab.temps;
  171. *
  172. ndime= (dime dtab1) - 1 ;
  173. *
  174. * ----- EXTRACTION DES DIFFERENTES TABLES
  175. *
  176. *                                    DEPLACEMENTS 
  177. *
  178. i=0 ;
  179. repeter bloc1 ndime ;
  180.    i=i+1 ;
  181.    dep1=tab.deplacements.i ;
  182.    deft=defo geotot dep1 1. ;
  183. **   trac      deft ;
  184. **   trac OEIL deft ;
  185.    sig1=tab.contraintes.i ;
  186.    var1=tab.variables_internes.i ;
  187.    defp=tab.deformations_inelastiques.i ;
  188.    eps1 = epsi modjoi dep1 ;
  189.    epscu = epsi modcub dep1 ;
  190. *
  191. * ----- REDUCTION DES CHPOINTS A L'ELEMENT JOINT 
  192. *
  193. *
  194. *
  195. *
  196. DPTCU1 = REDU DEP1  VOLS  ;
  197. SIGCU1 = REDU SIG1  VOLS  ;
  198. DEFCU1 = REDU EPSCU VOLS ;
  199. DPTCU2 = REDU DEP1  VOLI  ;
  200. SIGCU2 = REDU SIG1  VOLI  ;
  201. DEFCU2 = REDU EPSCU VOLI ;
  202. DPTJOI = REDU DEP1 JOINT ;
  203. SIGJOI = REDU SIG1 JOINT ;
  204. DEFJOI = REDU EPS1 JOINT ;
  205. DEPJOI = REDU DEFP JOINT ;
  206. VARJOI = REDU VAR1 JOINT ;
  207. *
  208. mess '               '   ;
  209. mess ' TEMPS  ' D ;
  210. mess ' ____________  '   ;
  211. mess '               '   ;
  212. mess ' deplacements du cube 1  '   ; 
  213. list dptcu1 ;
  214. mess '               '   ;
  215. mess ' deplacements du cube 2  '   ; 
  216. list dptcu2 ;
  217. mess '               '   ;
  218. mess ' deplacements du joint   '   ; 
  219. list dptjoi ;
  220. mess '               '   ;
  221. mess ' contraintes du cube 1    '   ; 
  222. list sigcu1 ;
  223. mess '               '   ;
  224. mess ' contraintes du cube 2    '   ; 
  225. list sigcu2 ;
  226. mess '               '   ;
  227. mess ' contraintes du joint  '   ; 
  228. list sigjoi ;
  229. mess '               '   ;
  230. mess ' def totales du joint'   ;
  231. list defjoi ;
  232. mess '               '   ;
  233. mess ' def inelasti. du joint '   ;
  234. list depjoi ;
  235. mess ' variables internes du joint'   ;
  236. list varjoi ;
  237. mess ' defo. elasti du cube du haut '   ;
  238. list defcu1 ;
  239. mess ' defo. elasti du cube du bas  '   ;
  240. list defcu2 ;
  241. *
  242. *
  243. DEP13 = EXTR DPTCU1 UX P3 ;
  244. DIF13 = DEP13-0.1;
  245. XDIF1 = ABS(DIF13);
  246.     SI (XDIF1 <EG 1.E-4) ;
  247.       ERRE 0 ;
  248.     SINON ;
  249.       ERRE 5 ;
  250.     FINSI ;
  251. *
  252. DEP23 = EXTR DPTCU1 UZ P3 ;
  253. DIF23 = DEP23-0.08;
  254. ZDIF1 = ABS(DIF13);
  255.     SI (ZDIF1 <EG 1.E-4) ;
  256.       ERRE 0 ;
  257.     SINON ;
  258.       ERRE 5 ;
  259.     FINSI ;
  260. *
  261. SMNJ  = EXTR SIGJOI SMN 1 1 1   ;
  262. ASMN  = ABS(SMNJ) ;
  263.     SI (ASMN <EG 1.E-4) ;
  264.       ERRE 0 ;
  265.     SINON ;
  266.       ERRE 5 ;
  267.     FINSI ;
  268. *
  269. SMZC  = EXTR SIGCU1 SMZZ 1 1 1   ;
  270. DIFZ  = SMZC ;
  271. ASZC  = ABS(DIFZ) ;
  272.     SI (ASZC <EG 1.E-4) ;
  273.       ERRE 0 ;
  274.     SINON ;
  275.       ERRE 5 ;
  276.     FINSI ;
  277. *
  278. fin bloc1 ;
  279. *
  280. FIN;
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  
  316.  
  317.  
  318.  

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