Télécharger la plaquette Cast3M


fr en

Télécharger stru1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  stru1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/06/13     PROVENANCE : TEST
  5. SAUT PAGE;
  6.  
  7. *********************************************************************
  8. *                                                                   *
  9. *   TEST STRU1                                                      *
  10. *                                                                   *
  11. *   PLAQUE RAIDIE SUR APPUIS SIMPLES, SOUS CHARGE UNIFORMEMENT      *
  12. *     REPARTIE                                                      *
  13. *                                                                   *
  14. *     La plaque rectangulaire de 2 metre de long, de 1 metre de     *
  15. *   large et 0.06 metre d'epaisseur est raidie a l'aide de 6        *
  16. *   raidisseurs de 0.01 metre d'epaisseur.                          *
  17. *                                                                   *
  18. *     La plaque superieure (uniquement) repose simplement sur les   *
  19. *   quatre cotes et est uniformement chargee.                       *
  20. *                                                                   *
  21. *     Pour des raisons de symetrie on ne maille que le quart de la  *
  22. *   structure.                                                      *
  23. *                                                                   *
  24. *     On se propose de chercher :                                   *
  25. *       - d'une part le deplacement vertical du milieu de la plaque,*
  26. *       - d'autre part le deplacement vertical maximal.             *
  27. *                                                                   *
  28. *   Ref. : JL. BATOZ,J.CHATELAIN,JP. JAMEUX                         *
  29. *         ETUDE D UNE PLAQUE RAIDIE                                 *
  30. *         BULLETIN DU CLUB S.A.P - VOLUME VIII, N1                 *
  31. *                                                                   *
  32. *********************************************************************
  33.  
  34. GRAPH = 'N';
  35.  
  36. SAUT PAGE;
  37. SI (NEG GRAPH 'N') ;
  38.   OPTI ECHO 1   ;
  39.   OPTI TRAC PSC ;
  40. SINO ;
  41.   OPTI ECHO 1   ;
  42. FINSI ;
  43.  
  44. SAUT PAGE;
  45. *
  46. OPTION DIME 3;
  47. OPTION ELEM QUA4;
  48.  
  49. *
  50. *  QUELQUES DONNEES
  51. *
  52. OEIL = 10000 -10000 10000;
  53. *
  54. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  55. *                             GEOMETRIE                              *
  56. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  57. *
  58. *            PLAQUE SUPERIEURE
  59. *
  60. PA = 0 .5 0; PB =1. .5 0; PD = 1. 0 0; PC = 0 0 0;
  61. PF = .6667 0 0; PE = .3333 0 0;
  62. *
  63. DAB = PA DROIT 18 PB; DBD = PB DROIT 10 PD;
  64. DDC = PD DROIT 18 PC; DCA = PC DROIT 10 PA;
  65. SU1 = DALLER DAB DBD DDC DCA PLAN;
  66. *
  67. *            RAIDISSEUR CENTRAL
  68. *
  69. PA2 = 0 .5 -.1; PC2 = 0 0 -.1;
  70. D1 = PA DROIT 10 PC;
  71. D3 = PC2 DROIT 10 PA2;
  72. D2 = PC DROIT 2 PC2;
  73. D4 = PA2 DROIT 2 PA;
  74. RA1 = DALL D1 D2 D3 D4 PLAN;
  75. *
  76. * SEGMENTS DE DROITES POUR CONDITIONS DE SYMETRIE
  77. D5 = D2 PLUS (.33333 0. 0.);
  78. D6 = D2 PLUS (.66667 0. 0.);
  79. *
  80. *            RAIDISSEURS 2 ET 3
  81. *
  82. RA2 = RA1 PLUS (.33333 0 0) ;
  83. RA3 = RA1 PLUS (.66667 0 0) ;
  84. *
  85. SURF = SU1 ET RA1 ET RA2 ET RA3;
  86. *
  87. ELIM SURF .0001;
  88. ELIM (SURF ET D1 ET D2 ET D3 ET D4 ET D5 ET D6) .001;
  89. *
  90. SI (NEG GRAPH 'N');
  91.   TITR ' TEST STRU1 : MAILLAGE';
  92.   TRAC CACH OEIL QUAL SURF;
  93. FINSI;
  94. *
  95. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  96. *      MODELE                                                   *
  97. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  98. *
  99. MODL1 = MODE RA1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  100. MODL2 = MODE RA2 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  101. MODL3 = MODE RA3 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  102. MODL4 = MODE SU1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  103. MODL5 = MODL1 ET MODL2 ET MODL3 ET MODL4;
  104. *
  105. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  106. *     MATERIAUX - CARACTERISTIQUES                              *
  107. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  108. *
  109. MATR1 = MATE MODL1 YOUN 2.E11 NU 0.3 EPAI 0.005 EXCEN 0.0025;
  110. MATR2 = MATE MODL2 YOUN 2.E11 NU 0.3 EPAI 0.01  EXCEN 0.0;
  111. MATR3 = MATE MODL3 YOUN 2.E11 NU 0.3 EPAI 0.01  EXCEN 0.0;
  112. MATR4 = MATE MODL4 YOUN 2.E11 NU 0.3 EPAI 0.006  EXCEN 0.0;
  113. MATR5 = MATR1 ET MATR2 ET MATR3 ET MATR4;
  114. *
  115. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  116. *             CONDITIONS AUX LIMITES                                *
  117. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  118. *
  119. CL1 = BLOQUE DAB UZ;
  120. CL2 = BLOQUE DAB RY;
  121. CL3 = BLOQUE DBD UZ;
  122. CL4 = BLOQUE DBD RX;
  123. CL5 = SYMT DEPL ROTA (0 0 0) (0 1 0) (0 0 1) SURF 0.0001;
  124. CL6 = SYMT DEPL ROTA (0 0 0) (1 0 0) (0 0 1) SURF 0.0001;
  125. CL = CL1 ET CL2 ET CL3 ET CL4 ET CL5 ET CL6;
  126. *
  127. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  128. *                 CHARGEMENT
  129. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  130. *
  131. MOP = 'MODE' SU1 'CHARGEMENT' 'PRESSION' 'COQ4' ;
  132. MAP = 'PRES' MOP 'PRES' 5E4 ;
  133. CHARG = 'BSIG' MOP MAP ;
  134. *
  135. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  136. *           RESOLUTION - CALCUL ELASTIQUE
  137. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  138. *
  139. RIG = RIGI MODL5 MATR5;
  140. RIG = RIG ET CL;
  141. *
  142. DEP = RESO RIG CHARG;
  143. *
  144. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  145. * RESULTATS
  146. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  147. *
  148. NOSURF = NBNO SURF ;
  149. ELSURF = NBEL SURF ;
  150. *
  151. SI (NEG GRAPH 'N');
  152.   AMPL1 = 1.E-5;
  153.   REA1 = REAC DEP RIG;
  154.   VE = VECT REA1 AMPL1 FX FY FZ BLEU ;
  155.   VF = VECT CHARG (AMPL1 * 20) FX FY FZ ROUGE;
  156.   VV = VE ET VF;
  157.   DEF0 = DEFO SURF DEP 0.;
  158.   DEF1 = DEFO SURF DEP 25. VV VERT;
  159.   TITR ' STRU1 : DEFORMEE, CHARGEMENT, REACTIONS' ;
  160.   TRACE OEIL (DEF0 ET DEF1);
  161.  
  162.   TITR 'STRU1 : Deplacements UZ';
  163.   TRAC CACH SURF (EXCO DEP UZ);
  164.  
  165.   DEF1 = DEFO SURF DEP 25. ROUGE;
  166.   OEIL = 0.25 -1.E6 0.;
  167.   TITR ' STRU1 : DEFORMEE';
  168.   TRACE OEIL (DEF0 ET DEF1);
  169. FINSI;
  170. *
  171. UZC = EXTR DEP UZ PC ;
  172. UZC1 = -0.475E-3;
  173. ERG1 = ABS ((UZC - UZC1) / UZC1) * 100 ;
  174. *
  175. DZ = EXCO DEP UZ;
  176. ZMIN = MAXI (ABS (DZ)); ZMIN = -1. * ZMIN;
  177. ZMIN1 = -0.126E-2;
  178. ERG2 = ABS ((ZMIN - ZMIN1) / ZMIN1) * 100 ;
  179. *
  180. SAUT PAGE ;
  182. *
  183. MESS ' RESULTATS ';
  184. MESS ' --------- ';
  185. *
  187. MESS ' DEPLACEMENT THEORIQUE AU POINT PC : ' UZC1 'METRES';
  188. MESS ' DEPLACEMENT CALCULE AU POINT PC   : ' UZC  'METRES' ;
  189. MESS '     SOIT UN ECART DE              : ' ERG1 '%' ;
  190. *
  192. MESS ' DEPLACEMENT THEORIQUE MAXIMUM (ABS) : ' ZMIN1 'METRES';
  193. MESS ' DEPLACEMENT CALCULE MAXIMUM (ABS)   : ' ZMIN 'METRES';
  194. MESS '     SOIT UN ECART DE                : ' ERG2 '%' ;
  195. *
  196. *SAUT 1 LIGNE;
  197. *MESS ' ABSICE DU POINT OU SE TROUVE LE MAXIMUM : ' ;
  198. *MESS '         - ABSICE THEORIQUE : ' XX1 'METRES';
  199. *MESS '         - ABSICE CALCULEE  : ' XX  'METRES';
  200. *MESS '     SOIT UN ECART DE       : '  ERG3  '%' ;
  201. *
  202. SAUT 1 LIGNE;
  203. MESS ' Nombre de noeuds  : ' NOSURF ;
  204. MESS ' Nombre d elements : ' ELSURF ;
  205. MESS '  Soit ' (NOSURF * 6) 'd.d.l.';
  206. SAUT 1 LIGNE;
  207. *
  208. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  209. *            CODE FONCTIONNEMENT                                     *
  210. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  211. *
  212. ERGMAX = MAXI (PROG ERG1 ERG2);
  213. *
  214. * Il y a erreur si un des pourcentages ERG* est suprieur  5%
  215. *
  216. SI (ERGMAX <EG 5.);
  217.     ERRE  0 ;
  218. SINO;
  219.     ERRE  5 ;
  220. FINSI;
  221. *
  223. TEMPS;
  225.  
  226. FIN;
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  

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