Télécharger la plaquette Cast3M


fr en

Télécharger elas4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  elas4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. *           Test elas4.dgibi: Jeux de données         *
  7. *           ---------------------------------         *
  8. *******************************************************
  9.  
  10. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  11. * SINON SI GRAPH DIFFERENT DE N TOUS
  12. * LES GRAPHIQUES SONT AFFICHES
  13.  
  14. GRAPH = 'N' ;
  15.  
  16. SAUT PAGE;
  17. SI (NEG GRAPH 'N') ;
  18.   OPTI ECHO 1   ;
  19.   OPTI TRAC PSC ;
  20. SINO ;
  21.   OPTI ECHO 0   ;
  22. FINSI ;
  23.  
  24. SAUT PAGE;
  25. *______________________________________________________
  26. *  TEST ELAS4                                         *
  27. *                                                     *
  28. *  POUTRE A SECTION EN FORME DE "Z" CHARGEE EN TORSION*
  29. *                                                     *
  30. *  Cas test NAFEMS numero LE5 (Z-section cantilever)  *
  31. *                                                     *
  32. *      La poutre est maillee a l'aide d'elements de   *
  33. *  coques COQ4.                                       *
  34. *  Le maillage a ete affine car on ne peut utiliser   *
  35. *  l'operateur CALP avec l'element COQ8.              *
  36. *                                                     *
  37. *      La poutre est soumis a un couple de torsion    *
  38. *  de 1200 m.N, obtenue par deux forces de 600 N      *
  39. *  uniformement reparties sur chaque flanc.           *
  40. *                                                     *
  41. *      On se propose de caculer la contrainte axiale  *
  42. *  de la surface au point A.                          *
  43. *                                                     *
  44. *______________________________________________________
  46. OPTI ELEM QUA4;
  47. *======================================================*
  48. *             geometrie : maillage                     *
  49. *======================================================*
  50. *
  51. * Dimension en metres
  52. *
  53. * Points
  54. *
  55. OEIL = 1000 -1000 1000;
  56. *
  57. P1 = 0 0 0; P21 = (P1 PLUS (10 0 0));
  58. P2 = 0 0 1; P22 = (P2 PLUS (10 0 0));
  59. P3 = 0 2 1; P23 = (P3 PLUS (10 0 0));
  60. P4 = 0 2 2; P24 = (P4 PLUS (10 0 0));
  61. *
  62. P1M = 10 0 .5; P2M = 10 2 1.5;
  63. *
  64. * Droites et lignes
  65. *
  66. * Remarque : ne pouvant utiliser d'element coq8,
  67. * et seulement des coq4 le maillage a ete affine.
  68. *
  69. N1 = 2;
  70. LIG1 = P1 D N1 P2 D N1 P3 D N1 P4;
  71. *
  72. * Maillage
  73. *
  74. N2 = 16;
  75. SUR1 = LIG1 TRANS N2 (10 0 0);
  76.  
  77. ELIM (SUR1 ET P21 ET P22 ET P23 ET P24 ET P1M ET  P2M) 0.001;
  78.  
  79. SI (NEG GRAPH 'N');
  80.    TITR ' TEST ELAS4 : MAILLAGE ';
  81.    TRAC OEIL QUAL SUR1;
  82. FINSI;
  83.  
  84. *====================================================*
  85. *           modele - materiau - caracteristique      *
  86. *           rigidite - conditions aux limites        *
  87. *====================================================*
  88. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  89. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3 EPAI 0.1;
  90. RI1   = RIGI MATR1 MODL1;
  91. CL1   = BLOQ DEPL LIG1;
  92. RIG1  = RI1 ET CL1;
  93.  
  94. *====================================================*
  95. *                      Chargement                    *
  96. *====================================================*
  97. *   Couple de torsion de 1.2e6 m.N applique dans la
  98. * section x = 10 m  par deux forces uniformement
  99. * reparties egales de 600 N sur chaque flanc.
  100.  
  101. F1  = 0 0 -150E3;
  102. F1M = 0 0 -300E3;
  103. F2  = 0 0  150E3;
  104. F2M = 0 0  300E3;
  105.  
  106. FORC1  = FORC F1  P21;
  107. FORC2  = FORC F1  P22;
  108. FORC1M = FORC F1M P1M;
  109. FORC3  = FORC F2  P23;
  110. FORC4  = FORC F2  P24;
  111. FORC2M = FORC F2M P2M;
  112.  
  113. CHA1 = FORC1 ET FORC2 ET FORC3 ET FORC4 ET FORC1M ET
  114. FORC2M;
  115.  
  116. *=====================================================*
  117. *            resolution : champs de deplacements      *
  118. *                         champs de contraintes       *
  119. *=====================================================*
  120. DEP1 = RESO RIG1 CHA1;
  121.  
  122. *
  123. * Trace facultatif de la deformee et des reactions
  124. SI (NEG GRAPH 'N');
  125.   AMPVEC = .5E-6;
  126.   REA1 = REAC DEP1 RIG1;
  127.   VECT1 = VECT CHA1 (AMPVEC * 20) FX FY FZ ROUGE;
  128.   VECT2 = VECT REA1 AMPVEC FX FY FZ BLEU;
  129.   DEF0 = DEFO 0.  DEP1 SUR1;
  130.   DEF1 = DEFO 10. DEP1 SUR1 (VECT1 ET VECT2) VERT;
  131.   TITR ' ELAS4 : Deformee, chargement, reactions';
  132.   TRAC OEIL SUR1 (DEF0 ET DEF1);
  133. FINSI;
  134.  
  135. *======================================================*
  136. *                  extraction d'une contrainte         *
  137. *            et comparaison avec solution analytique   *
  138. *======================================================*
  139. *     On cherche la contrainte axiale au point A sur la
  140. * surface moyenne.
  141.  
  142. * Calcul des contraintes generalisees dans le repere
  143. * local des elements
  144. CHAM1 = SIGM MODL1 MATR1 DEP1;
  145.  
  146. * Calcul du champs de contraintes au sens des milieux
  147. * continus.
  148. *   Pour cela on utilise l'operateur CALP.
  149. CHAM2 = CALP CHAM1 MATR1 MODL1 MOYEN;
  150. CHAM3 = CHAN NOEUD CHAM2 MODL1;
  151. CHPO3 = CHAN CHPO CHAM3 MODL1;
  152. PA = SUR1 POIN PROC (2.5 0 0);
  153. *
  154. * Trace facultatif des contraintes aux noeuds.
  155. *
  156. SI (NEG GRAPH 'N');
  157.    TITR 'ELAS4 : Maillage';
  158.    TRAC CACH SUR1 'QUAL';
  159.    TITR 'ELAS4 : Contraintes Nodales';
  160.    TRAC CACH SUR1 ((EXCO CHPO3 SMYY)*1e-6);
  161. FINSI;
  162.  
  163. SMXXA2 = EXTR CHPO3 SMYY PA;
  164. SMXXA2 = SMXXA2 / 1.E6;
  165. SMXXA1 = -108;
  166. ERGXXA = 100 * (ABS ((SMXXA2 - SMXXA1) / SMXXA1));
  167.  
  168. *=======================================================*
  169. *                   affichage des resultats             *
  170. *=======================================================*
  171. MESS ' RESULTATS TEST ELAS4 ';
  172. MESS ' -------------------- ';
  174. MESS '   On cherche la contrainte axiale sur la surface moyenne';
  175. MESS ' au point A.';
  176. MESS '   Elle est comparee a une valeur theorique obtenue';
  177. MESS ' analytiquement.';
  179. MESS ' Contrainte axiale theorique en A :' SMXXA1 'MPa';
  180. MESS ' Contrainte axiale calculee en A  :' SMXXA2 'MPa';
  181. MESS '    Soit un ecart de :' ERGXXA '%';
  183.  
  184. ELSUR1 = NBEL SUR1;
  185. NOSUR1 = NBNO SUR1;
  186.  
  187. MESS '    Nombre d elements :' ELSUR1;
  188. MESS '    Nombre de noeuds  :' NOSUR1;
  189.  
  190. *========================================================*
  191. *                  code fonctionnement                   *
  192. *========================================================*
  193. SAUT 2 LIGNE;
  194. SI (ERGXXA <EG 5);
  195.     ERRE 0;
  196. SINON;
  197.     ERRE 5;
  198. FINSI;
  199.  
  200. TEMPS;
  201. FIN;
  202.  
  203.  
  204.  
  205.  
  206.  

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