Télécharger la plaquette Cast3M


fr en

Télécharger test_CHAB_SINH_X.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : test_CHAB_SINH_X.dgibiq
  2. ************************************************************************
  3. ************************************************************************
  4. **--------------------------------------
  5. **Nouvelle loi de comportement mécanique 
  6. **--------------------------------------
  7.  
  8.  
  9. **Données
  10. **-------
  11. *longueur de la plaque
  12. Lo=100.;
  13. *largeur de la demi plaque
  14. La=20.;
  15.  
  16. **Points du mailage
  17. **-----------------
  18. OPTI 'DIME' 2 'ELEM' QUA4 MODE AXIS;
  19.  
  20. A1=0. 0.;
  21. A2=la 0.;
  22. A3=la Lo;
  23. A4=0. Lo;
  24. L1=DROITE 1 A4 A1;
  25. L2=DROITE 1 A1 A2;
  26. L3=DROITE 1 A2 A3;
  27. L4=DROITE 1 A3 A4;
  28. cube=DALL L1 L2 L3 L4;
  29. *TRAC cube;
  30.  
  31. **Définition modèle et matériau
  32. **-----------------------------
  33. *modèle mécanique
  34. MODM=MODEL cube MECANIQUE ELASTIQUE VISCOPLASTIQUE 'CHAB_SINH_X';
  35. MATM=MATER MODM 'YOUN' 180014.4716 'NU' 0.3 'SIGY' 359.8729
  36. 'K' 2.2 'N' 0.6 'B' 23.1521 'RI' 40.3074 'epso' 6.E-8 'C1' 16535.60 
  37. 'gam1' 165.5549 'C2' 784.2466;
  38.  
  39. **Blocages mécaniques
  40. **-------------------
  41. CLM1=BLOQ 'UZ' L2;
  42. CLM2=BLOQ 'UZ' L4;
  43.  
  44. CL_M=CLM1 ET CLM2;
  45.  
  46. **Déplacements imposés
  47. **--------------------
  48. def1=2.;
  49. def2=4.;
  50. dep1=(def1/100) * Lo;
  51. dep2=(def2/100) * Lo;
  52. V=0.0001;
  53. T_def1=(def1/100)/V;
  54. T_def2=(def2/100)/V;
  55.  
  56. *-------------------------
  57. *boucles à la première def
  58. *-------------------------
  59. EVT=PROG 0.;
  60. EVDEF=PROG 0;
  61. EVT=EVT ET (PROG (T_def1 / 2.));
  62. EVDEF=EVDEF ET (PROG dep1);
  63. REPE bou 8;
  64.   ite=&bou;
  65.   EVT=EVT ET (PROG ((EXTR EVT (DIME EVT)) + T_def1));
  66.   EVDEF=EVDEF ET (PROG ((-1) ** (ite) * dep1));
  67. FIN bou;
  68. REPE bou 6;
  69.   ite=&bou;
  70.   EVT=EVT ET (PROG ((EXTR EVT (DIME EVT)) + T_def2));
  71.   EVDEF=EVDEF ET (PROG ((-1) ** (ite) * dep2));
  72. FIN bou;
  73.  
  74. DI1=DEPI CLM2 1.;
  75. EV2=EVOL MANU 'temps' EVT 'deplacement' EVDEF;
  76. CHAM1=CHAR 'DIMP' DI1 EV2;
  77.  
  78.  
  79. **Mise en place de la table PASAPAS
  80. **---------------------------------
  81. tm=table;
  82. tm.MES_SAUVEGARDES=TABLE;
  83. tm.MES_SAUVEGARDES.'DEFTO'=VRAI;
  84. tm.MES_SAUVEGARDES.'DEFIN'=VRAI;
  85. tm.precision=1.E-3;
  86. tm.modele=MODM;
  87. tm.caracteristiques=MATM;
  88. tm.blocages_mecaniques=CL_M;
  89. tm.chargement=CHAM1;
  90. tm.temps_calcules=prog 0. pas (T_def1/10.) (EXTR EVT (DIME EVT));
  91.  
  93.  
  94. **Analyse des champs résultants
  95. **-----------------------------
  96. IDIM=(DIME (tm.temps)) - 1;
  97. *déformée
  98. DEPF=tm.deplacements.IDIM;
  99. DEFI=DEFO DEPF cube 0. VERT;
  100. DEFF=DEFO DEPF cube 1. ROUGE;
  101. *TRAC (DEFI ET DEFF) TITRE 'déformée';
  102. *Contraintes de V.Mises
  103. CO=tm.contraintes.IDIM;
  104. COP2=VMIS (MODM) CO;
  105. *TRAC DEFF COP2 (MODM) TITRE 'Contrainte de Von Mises';
  106. *tracé des isovaleurs des déplacements
  107. ZZ1=EXCO DEPF 'UZ';
  108. OPTI ISOV 'SURFACE';
  109. *TRAC cube ZZ1 TITRE 'déplacement Uz';
  110. XX1=EXCO DEPF 'UR';
  111. *TRAC cube XX1 TITRE 'déplacement Ur';
  112. *tracé de la contrainte SMZZ
  113. TCON=EXCO tm.contraintes.IDIM 'SMZZ';
  114. *TRAC TCON (MODM) TITRE 'contrainte suivant zz';
  115. *écrouissage isotrope
  116. EI=EXCO tm.variables_internes.IDIM 'r';
  117. *TRAC cube MODM EI TITRE 'écrouissage isotrope';
  118. *déformation plastique cumulée
  119. DPC=EXCO tm.variables_internes.IDIM 'EPSE';
  120. *TRAC cube MODM DPC TITRE 'déformation plastique cumulée';
  121.  
  122.  
  123. **caractéristiques en fonction du temps
  124. **-------------------------------------
  125. *initilisation
  126. DEPP5=PROG 0.;
  127. CONPC=PROG 0.;
  128. DEFTPC=PROG 0.;
  129. SEC=PROG 0.;
  130. RPC=PROG 0.;
  131. PPC=PROG 0.;
  132. DEI=PROG 0.;
  133. X1=PROG 0.;
  134. X2=PROG 0.;
  135. *boucle
  136. REPETER bou IDIM;
  137. ite=&bou;
  138. TE=tm.temps.ite;
  139. DEP=tm.deplacements.ite;
  140. CON=tm.contraintes.ite;
  141. DEFT=tm.deformations.ite;
  142. VARI2=tm.variables_internes.ite;
  143. DEFI=tm.deformations_inelastiques.ite;
  144. *echelle de temps
  145. SEC=SEC ET (PROG TE);
  146. *déplacement du bout de l'éprouvette
  147. DEPP5n=EXTR DEP 'UZ' A3;
  148. DEPP5=DEPP5 ET (PROG DEPP5n);
  149. *contrainte au point de contrôle
  150. CONPCn=EXTR CON 'SMZZ' 1 1 1;
  151. CONPC=CONPC ET (PROG CONPCn);
  152. *déformation totale au point de contrôle
  153. DEFTPCn=EXTR DEFT 'EPZZ' 1 1 1;
  154. DEFTPC=DEFTPC ET (PROG DEFTPCn);
  155. *écrouissage au point de contrôle
  156. RPCn=(EXTR VARI2 'r' 1 1 1) * 40.3074 * 23.1521;
  157. RPC=RPC ET (PROG RPCn);
  158. *p au point de contrôle
  159. PPCn=EXTR VARI2 'EPSE' 1 1 1;
  160. PPC=PPC ET (PROG PPCn);
  161. *deformation inelastique suivant Uz
  162. DEIn=EXTR DEFI 'EIZZ' 1 1 1;
  163. DEI=DEI ET (PROG DEIn);
  164. *écrouissage cinématique non linéaire suivant z:
  165. X1n=(EXTR VARI2 'A1ZZ' 1 1 1) * 16535.60;
  166. X1=X1 ET (PROG X1n);
  167. *écrouissage cinématique linéaire suivant z:
  168. X2n=(EXTR VARI2 'A2ZZ' 1 1 1) * 784.2466;
  169. X2=X2 ET (PROG X2n);
  170. FIN bou;
  171.  
  172. **graphiques
  173. **----------
  174. *deplacement(t)
  175. EVDEPB=EVOL 'MANU' 'temps' SEC 'déplacement' DEPP5;
  176. *DESS EVDEPB TITRE 'deplacement du bout de la poutre (t)';
  177. *contrainte suivant z au point de contrôle (t)
  178. EVCONPC=EVOL 'MANU' 'temps' SEC 'contrainte suivant z' CONPC;
  179. *DESS EVCONPC TITRE 'contrainte suivant z au point de contrôle';
  180. *tracé de la déformation totale au point de contrôle (t)
  181. EVDEFTPC=EVOL 'MANU' 'temps' SEC 'déformation totale' DEFTPC;
  182. *DESS EVDEFTPC TITRE 'déformation totale au point de contrôle';
  183. *tracé de la contrainte en fonction de la déformation
  184. EVCD=EVOL 'MANU' 'epsilon' DEFTPC 'contrainte' CONPC;
  185. *DESS EVCD TITRE 'contrainte/deformation totale';
  186. *tracé de l'écrouissage en fonction de p
  187. EVRD=EVOL 'MANU' 'p' PPC 'écrouissage' RPC;
  188. *DESS EVRD TITRE 'ecrouissage/p';
  189. *tracé de p en fonction de la déformation
  190. EVPD=EVOL 'MANU' 'epsilon' DEFTPC 'p' PPC;
  191. *DESS EVPD TITRE 'p/deformation totale';
  192. *tracé de la deformation inelastique
  193. EVDEFI=EVOL 'MANU' 'epsilon' DEFTPC 'deformation inelastique' DEI;
  194. *DESS EVDEFI TITRE 'deformation inelastique/deformation totale';
  195. *tracé de X1
  196. EVX1=EVOL 'MANU' 'p' PPC 'X1' X1;
  197. *DESS EVX1 TITRE 'X1/p';
  198. *tracé de X2
  199. EVX2=EVOL 'MANU' 'p' PPC 'X2' X2;
  200. *DESS EVX2 TITRE 'X2/p';
  201.  
  202. **sauvegarde des résultats
  203. **------------------------
  204. @EXCEL1 EVCD 'et-contrainte.xls';
  205. @EXCEL1 EVRD 'p-R.xls';
  206. @EXCEL1 EVX1 'p-X1.xls';
  207. @EXCEL1 EVX2 'p-X2.xls';
  208.  
  209. FIN; 
  210.  
  211.  
  212.  
  213.  

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