Télécharger plas4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Plas4.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test plas4.dgibi: Jeux de données *
  10. * --------------------------------- *
  11. *******************************************************
  12.  
  13. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  14. * SINON SI GRAPH DIFFERENT DE N TOUS
  15. * LES GRAPHIQUES SONT AFFICHES
  16.  
  17. GRAPH = 'N' ;
  18.  
  19. SAUT PAGE;
  20. SI (NEG GRAPH 'N') ;
  21. OPTI ECHO 1 ;
  22. OPTI TRAC PSC ;
  23. SINO ;
  24. OPTI ECHO 0 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *******************************************************
  29. * CUBE EN TRACTION UNIAXIALE
  30. *
  31. * TEST ELEMENTAIRE DU GROUPE DE TRAVAIL
  32. * 'STATIQUE NON LINEAIRE' COMMISSION VPCS
  33. *
  34. * LE CALCUL MARCHE ET DONNE EXACTEMENT LES
  35. * RESULTATS THEORIQUES LE 05/05/90
  36. *
  37. * Un cube constitué de 2 éléments CUB8 est soumis à
  38. * une traction uniaxiale (blocages des déplacements
  39. * suivant Y et Z et déplacement imposé suivant X
  40. * aux deux extrémités (surfaces SU1 et SU2))
  41. *
  42. * Les contraintes obtenues sont comparées aux
  43. * résultats théoriques et donnent les résultats
  44. * précisés plus bas.
  45. *
  46. *******************************************************
  47. *
  48.  
  49. SAUT PAGE;
  50. TITRE 'CUBE EN TRACTION UNIAXIALE';
  51. OPTION DIME 3 ECHO 0;
  52. TEMPS ;
  53. DENS 1;
  54. *------------------------------------------------------
  55. ***** DEFINITION DE LA GEOMETRIE
  56. *------------------------------------------------------
  57. A = 1.;
  58. A1 = -1.;
  59. A2 = .5;
  60. B = 1.;
  61. C = 1.;
  62. P1 = A 0 0; Q1 = A2 0 0; R1 = 0 0 0;
  63. P2 = A B 0; Q2 = A2 B 0; R2 = 0 B 0;
  64. P3 = A B C; Q3 = A2 B C; R3 = 0 B C;
  65. P4 = A 0 C; Q4 = A2 0 C; R4 = 0 0 C;
  66. TR = A1 0 0;
  67. OPTION ELEM SEG2;
  68. LI1 = P1 DROIT 1 P2;
  69. LI2 = P2 DROIT 1 P3;
  70. LI3 = P3 DROIT 1 P4;
  71. LI4 = P4 DROIT 1 P1;
  72. CO1 = LI1 ET LI2 ET LI3 ET LI4;
  73. LJ1 = Q1 DROIT 1 Q2;
  74. LJ2 = Q2 DROIT 1 Q3;
  75. LJ3 = Q3 DROIT 1 Q4;
  76. LJ4 = Q4 DROIT 1 Q1;
  77. CO2 = LJ1 ET LJ2 ET LJ3 ET LJ4;
  78. LK1 = R1 DROIT 1 R2;
  79. LK2 = R2 DROIT 1 R3;
  80. LK3 = R3 DROIT 1 R4;
  81. LK4 = R4 DROIT 1 R1;
  82. CO3 = LK1 ET LK2 ET LK3 ET LK4;
  83. OPTION ELEM QUA4;
  84. SU1 = SURFACE CO1 PLANE;
  85. SU2 = SURFACE CO2 PLANE;
  86. SU3 = SURFACE CO3 PLANE;
  87. OPTION ELEM CUB8;
  88. VOL1 = SU1 VOLU 1 SU2;
  89. VOL2 = SU2 VOLU 1 SU3;
  90. VOL = VOL1 ET VOL2;
  91. ELIM .001 VOL;
  92. *------------------------------------------------------
  93. ***** TRACE
  94. *------------------------------------------------------
  95. OEIL = (30 20 20);
  96. SI(NEG GRAPH 'N');
  97. TRAC OEIL 'QUAL' VOL;
  98. FINSI;
  99. *------------------------------------------------------
  100. ***** CONDITIONS AUX LIMITES
  101. *------------------------------------------------------
  102. CL1 = BLOQUE SU3 UX;
  103. CL2 = BLOQUE VOL UY UZ;
  104. CL3 = BLOQUE SU2 UX;
  105. CL4 = BLOQ SU1 UX;
  106. CL = CL1 ET CL2 ET CL3 ET CL4;
  107. *------------------------------------------------------
  108. * VALEURS IMPOSEES DES DEPLACEMENTS
  109. *------------------------------------------------------
  110. EPS0 = .25E-4;
  111. COEF1 = A * EPS0 * .5;
  112. COEF2 = A * EPS0;
  113. D1 = DEPI CL3 COEF1;
  114. D2 = DEPI CL4 COEF2;
  115. DEP = D1 ET D2;
  116. *------------------------------------------------------
  117. ***** CALCUL MODELE CARACTERISTIQUES
  118. *------------------------------------------------------
  119. MO=MODE VOL MECANIQUE ELASTIQUE PLASTIQUE PARFAIT CUB8;
  120. MAT1 = MATE MO YOUN 2.5E11 SIGY 5.0E6 NU .25;
  121. *------------------------------------------------------
  122. ***** CHAMPS RELATIFS A L'OPERATEUR NONLIN
  123. *------------------------------------------------------
  124. LP1 = PROG 0. 1. 2. 3. 4.;
  125. LP2 = PROG 0. 1. 2. 1. 0.;
  126. EV = EVOL MANU T LP1 F(T) LP2;
  127. CHA1 = CHAR DIMP DEP EV;
  128. LIS1 = PROG 1. 2. 3. 4.;
  129. *------------------------------------------------------
  130. ***** CALCUL NON LINEAIRE
  131. *------------------------------------------------------
  132. TAB1 = TABLE;
  133. TAB1.'CARACTERISTIQUES' = MAT1;
  134. TAB1.'CHARGEMENT' = CHA1;
  135. TAB1.'MODELE' = MO;
  136. TAB1.'TEMPS_CALCULES' = LIS1;
  137. TAB1.'BLOCAGES_MECANIQUES' = CL;
  138. PASAPAS TAB1;
  139. *------------------------------------------------------
  140. ***** EXTRACTION DES RESULTATS
  141. *------------------------------------------------------
  142. SI1 = TAB1.CONTRAINTES.1;
  143. SI2 = TAB1.CONTRAINTES.2;
  144. SI3 = TAB1.CONTRAINTES.3;
  145. SI4 = TAB1.CONTRAINTES.4;
  146. SXX1 = EXTR SI1 SMXX 1 1 1;
  147. SYY1 = EXTR SI1 SMYY 1 1 1;
  148. SZZ1 = EXTR SI1 SMZZ 1 1 1;
  149. SXX2 = EXTR SI2 SMXX 1 1 1;
  150. SYY2 = EXTR SI2 SMYY 1 1 1;
  151. SZZ2 = EXTR SI2 SMZZ 1 1 1;
  152. SXX3 = EXTR SI3 SMXX 1 1 1;
  153. SYY3 = EXTR SI3 SMYY 1 1 1;
  154. SZZ3 = EXTR SI3 SMZZ 1 1 1;
  155. SXX4 = EXTR SI4 SMXX 1 1 1;
  156. SYY4 = EXTR SI4 SMYY 1 1 1;
  157. SZZ4 = EXTR SI4 SMZZ 1 1 1;
  158.  
  159.  
  160. SI(NEG GRAPH 'N');
  161. TITR 'Contrainte SMYY increment 4 en MPa';
  162. TRAC OEIL ((EXCO SMYY SI4)/1e6) MO;
  163. FINSI;
  164. *------------------------------------------------------
  165. ***** AFFICHAGE DES RESULTATS
  166. *------------------------------------------------------
  167. SAUT PAGE;
  168. *////////// PARAMETRE D'EVOLUTION T : 1 //////////';
  169. *------------------------------------------------------
  170. MESS ' T = 1 ';
  171. MESS ' RESULTATS THEORIQUES : SIGMAXX = 7.5 E6' ;
  172. MESS ' SIGMAYY = 2.5 E6' ;
  173. MESS ' SIGMAZZ = 2.5 E6' ;
  174. MESS ' ' ;
  175. MESS ' RESULTATS CALCULES : SIGMAXX = 'SXX1 ;
  176. MESS ' SIGMAYY = 'SYY1 ;
  177. MESS ' SIGMAZZ = 'SZZ1 ;
  178. SAUT 1 LIGNE ;
  179.  
  180. RESX1 = 100 * (ABS((SXX1-7.5e6)/7.5e6));
  181. RESY1 = 100 * (ABS((SYY1-2.5e6)/2.5e6));
  182. RESZ1 = 100 * (ABS((SZZ1-2.5e6)/2.5e6));
  183. *
  184. *////////// PARAMETRE D'EVOLUTION T : 2 //////////';
  185. *------------------------------------------------------
  186. MESS ' T = 2 ';
  187. MESS ' RESULTATS THEORIQUES : SIGMAXX = 11.66 E6' ;
  188. MESS ' SIGMAYY = 6.66 E6' ;
  189. MESS ' SIGMAZZ = 6.66 E6' ;
  190. MESS ' ' ;
  191. MESS ' RESULTATS CALCULES : SIGMAXX = 'SXX2 ;
  192. MESS ' SIGMAYY = 'SYY2 ;
  193. MESS ' SIGMAZZ = 'SZZ2 ;
  194. SAUT 1 LIGNE ;
  195. RESX2 = 100 * (ABS((SXX2-11.66e6)/11.66e6));
  196. RESY2 = 100 * (ABS((SYY2-6.66e6) /6.66e6 ));
  197. RESZ2 = 100 * (ABS((SZZ2-6.66e6) /6.66e6 ));
  198. *
  199. *////////// PARAMETRE D'EVOLUTION T : 3 //////////';
  200. *------------------------------------------------------
  201. MESS ' T = 3 ' ;
  202. MESS ' RESULTATS THEORIQUES : SIGMAXX = 4.16 E6' ;
  203. MESS ' SIGMAYY = 4.16 E6' ;
  204. MESS ' SIGMAZZ = 4.16 E6' ;
  205. MESS ' ' ;
  206. MESS ' RESULTATS CALCULES : SIGMAXX = 'SXX3 ;
  207. MESS ' SIGMAYY = 'SYY3 ;
  208. MESS ' SIGMAZZ = 'SZZ3 ;
  209. SAUT 1 LIGNE ;
  210. RESX3 = 100 * (ABS((SXX3-4.16e6) /4.16e6));
  211. RESY3 = 100 * (ABS((SYY3-4.16e6) /4.16e6 ));
  212. RESZ3 = 100 * (ABS((SZZ3-4.16e6) /4.16e6 ));
  213.  
  214. *
  215. *////////// PARAMETRE D'EVOLUTION T : 4 //////////';
  216. *------------------------------------------------------
  217. MESS ' T = 4 ' ;
  218. MESS ' RESULTATS THEORIQUES : SIGMAXX = -3.33 E6' ;
  219. MESS ' SIGMAYY = 1.66 E6' ;
  220. MESS ' SIGMAZZ = 1.66 E6' ;
  221. MESS ' ' ;
  222. MESS ' RESULTATS CALCULES : SIGMAXX = 'SXX4 ;
  223. MESS ' SIGMAYY = 'SYY4 ;
  224. MESS ' SIGMAZZ = 'SZZ4 ;
  225. SAUT 1 LIGNE ;
  226. TEMPS ;
  227. RESX4 = 100 * (ABS((SXX4+3.33e6) /-3.33e6 ));
  228. RESY4 = 100 * (ABS((SYY4-1.66e6) / 1.66e6 ));
  229. RESZ4 = 100 * (ABS((SZZ4-1.66e6) / 1.66e6 ));
  230.  
  231. *------------------------------------------------------
  232. * CODE DE FONCTIONNEMENT
  233. *------------------------------------------------------
  234. RESID=MAXI (PROG RESX1 RESX2 RESX3 RESY1 RESY2 RESY3 RESZ1 RESZ2 RESZ3);
  235.  
  236. SI(RESID <EG 0.5);
  237. ERRE 0;
  238. SINO;
  239. ERRE 5;
  240. FINSI;
  241.  
  242. FIN;
  243.  
  244.  
  245.  
  246.  
  247.  

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