Télécharger la plaquette Cast3M


fr en

Télécharger hotan.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : hotan.dgibi
  2. *Option de calcul
  3. OPTI 'DIME' 2 ELEM TRI3 MODE PLAN DEFO;
  4. *OPTI TRAC OPEN;
  5. OPTI 'SAUV' 'FORMAT' 'XDR' SAVING;
  6. OPTI SORT 'RES' ;
  7.  
  8. ** Parametres geometriques
  9. L1 = 1. ;
  10. R1 = 0.15 ;
  11. ** Maillage
  12. *DENmail0 = R1 / 15. ;
  13. *DENS2 = L1 / 15. ;
  14. P0 = 0. 0. ;
  15. P1 = R1 0. ;
  16. P2 = L1 0. ;
  17. P3 = L1 L1 ;
  18. P4 = 0. L1 ;
  19. P5 = 0. R1 ;
  20. CER = CERC 15 P5 P0 P1 ;
  21. L1 = DROI 20 P1 P2 ;
  22. L2 = DROI 20 P2 P3 ;
  23. L3 = DROI 20 P3 P4 ;
  24. L4 = DROI 20 P4 P5 ; 
  25. C1 = CER ET L1 ET L2 ET L3 ET L4;
  26. mail0 = SURF PLAN C1;
  27. *Tracer le maillage 
  28. *TRAC mail0;
  29. *Nombre de noeuds et d'éléments
  30. nbno0 = nbno mail0 ; mess 'Nombre de noeuds :' nbno0 ;
  31. nbel0 = nbel mail0 ; mess 'Nombre d éléments :' nbel0 ;
  32.  
  33. ** Modèle élasto-plastique Von Mises à l'écrouissage isotrope
  34. * Courbe de traction à déclarer. Limite d'élasticité 237 MPa
  35. x = prog  0. 0.00112860 0.024 0.031 0.039 0.047 0.054
  36.           0.062 0.070 0.077 0.086 0.095 0.104 0.112 ;
  37. *0.00112860
  38. y = prog   0. 237. 266.  282.  296.  308.  318. 
  39.           327.  335.  342.  350.  355.  360.  364.;
  40.  
  41. *opti donn 5 ;
  42.  
  43. * Courbe d'ecrouissage : a priori, le 1er point donne Sy
  44. epsy1 = extr x 2 ;
  45. lep   = (x enle 1) - epsy1 ;
  46. lsig1 = y enle 1 ;
  47. x     = lep ; y = lsig1 ;
  48.  
  49. *Conversion de l'unité de contraintes (MPa)
  50. y = y * 1.E6;
  51.  
  52. *Créer et tracer la courbe de traction
  53. xy = evol roug manu x y ;
  54. *dess xy titx 'Déformation' tity 'Contrainte [MPa]' ;
  55.  
  56. mod0 = MODELE mail0 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'ISOTROPE';
  57. mat0 = MATE mod0 'YOUN' 210E9 'NU' 0.30 'ECRO' xy;
  58.  
  59. *Modèle J2 pour contourner la déclaration d'une courbe de traction
  60. *mod0 = MODELE mail0 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'J2';
  61. *mat0 = MATE mod0 'YOUN' 210E9 'NU' 0.30 'SIG0' 250E6 'KISO' 21E9 'SIGI' 400E6 'VELO' 0;
  62.  
  63. *mod0 = MODELE mail0 MECANIQUE 'ELASTIQUE' 'NON_LINEAIRE' 'EQUIPLAS'; ce truc est pour l'élasticité non linéaire
  64.  
  65.  
  66. ** Conditions aux limites
  67. CLB = BLOQ L1 'UY';
  68. CLG = BLOQ L4 'UX';
  69. CL = CLB ET CLG;
  70.  
  71. ** Chargement évolutive
  72. *On se donne une sollicitation pendant 10s. Une charge/décharge qui monte linéairement 
  73. *jusqu'à 250MPa et puis descend à 0. 
  74. timestep = 0.5;
  75. time1 = PROG 0.0 PAS timestep 5.0 ;
  76. time2 = PROG 5.0 PAS timestep 10.0 ;
  77. time = PROG 0.0 PAS timestep 10.0;
  78. Pext1 = prog 'LINE' 'A' 45E6 'B' 0 time1;
  79. Pext2 = prog 'LINE' 'A' -45E6 'B' 450E6 time2;
  80.  
  81. ev1 = evol manu 'Temps' time1 'Pression' Pext1;
  82. ev2 = evol manu 'Temps' time2 'Pression' Pext2;
  83. *dess ev1;
  84. *dess ev2;
  85. ev = CONCAT ev1 ev2;
  86. *dess ev;
  87. PRES1 = PRES 'MASS' mod0 -1.0 L2;
  88. Charge1 = CHAR MECA PRES1 ev;
  89. Vec = VECT PRES1 'FORC';
  90. *TRAC Vec mail0;
  91.  
  92. RIG1 = RIGI mod0 mat0;
  93. RIG1 = RIG1 ET CL;
  94. * Initialisation de la table de calcul
  95.  SAVING_TBL = TABLE ;
  96.  SAVING_TBL . 'DEFIN' = VRAI;
  97.  SAVING_TBL . 'DEFTO' = VRAI;
  98.  SAVING_TBL . 'DEFLO' = FAUX;
  99.  SAVING_TBL . 'ROTAF' = FAUX;
  100.  
  101.  
  102. tab1 = TABLE ;
  103. tab1.BLOCAGES_MECANIQUES = CL ;
  104. tab1.MODELE = mod0 ;
  105. tab1.RIGI=RIG1;
  106. tab1.CARACTERISTIQUES = mat0 ;
  107. tab1.CHARGEMENT = Charge1 ;
  108. tab1.TEMPS_CALCULES = time ;
  109. tab1.TEMPS_SAUVES = time;
  110. tab1.TEMPS_SAUVEGARDES = time;
  111. tab1.MES_SAUVEGARDES=SAVING_TBL;
  112. *tab1.K_TANGENT = VRAI;
  113. *tab1.K_TANGENT_SYME = VRAI;
  114.  
  115. ** Calcul
  116. PASAPAS tab1;
  117.  
  118. ** Post-traitement (juste pour tester)
  119. dd1 = tab1 . DEPLACEMENTS ;
  120.  
  121. *
  122. td1 = dd1 . 0 ;
  123. td2 = dd1 . 9 ;
  124. td3 = dd1 . 20 ;
  125. *
  126. Def0 = DEFO mail0 td1 100. BLANC ;
  127. Def1 = DEFO mail0 td2 100. ROUGE;
  128. def2 = DEFO mail0 td3 100. BLEU ;
  129. TITRE 'Déformée Blanc = ini , Rouge = milieu , Bleu = fin ' ;
  130.  
  131. *TRAC ( Def0 ET Def1 ET Def2 ) QUAL ;
  132.  
  133.  
  134. *Export des résultats vers un fichier inp.
  135.  
  136. Pasmabou = 1 ;
  137.  
  138. REPETER mabou 20  ;
  139.  
  140. SIG_G = tab1 . CONTRAINTES . Pasmabou ;
  141. EPS_G = tab1 . DEFORMATIONS . Pasmabou;
  142. EPSP_G = tab1 . VARIABLES_INTERNES . Pasmabou;
  143. *vmis1 = VMIS mod0 SIG_G ;
  144.  
  145. SIGXX = EXCO SIG_G SMXX;
  146. SIGYY = EXCO SIG_G SMYY;
  147. SIGXY = EXCO SIG_G SMXY;
  148. EPSXX = EXCO EPS_G EPXX;
  149. EPSXY = EXCO EPS_G EPYY;
  150. GAMXY = EXCO EPS_G GAXY;
  151. *EPSXY = GAMXY/2;
  152.  
  153. MAT_TAN = HOTA mod0 SIG_G EPSP_G mat0;
  154.  
  155. SORT 'AVS' mail0 SIGXX 'SUIT' 'TEMP' Pasmabou;
  156. SORT 'AVS' mail0 SIGYY 'SUIT' 'TEMP' Pasmabou;
  157. SORT 'AVS' mail0 SIGXY 'SUIT' 'TEMP' Pasmabou;
  158.  
  159.  
  160. Pasmabou = Pasmabou + 1 ;
  161. FIN mabou;
  162.  
  163. FIN;
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  

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