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 ;
  181. *
  182. MESS ' RESULTATS ';
  183. MESS ' --------- ';
  184. *
  185. MESS ' DEPLACEMENT THEORIQUE AU POINT PC : ' UZC1 'METRES';
  186. MESS ' DEPLACEMENT CALCULE AU POINT PC : ' UZC 'METRES' ;
  187. MESS ' SOIT UN ECART DE : ' ERG1 '%' ;
  188. *
  189. MESS ' DEPLACEMENT THEORIQUE MAXIMUM (ABS) : ' ZMIN1 'METRES';
  190. MESS ' DEPLACEMENT CALCULE MAXIMUM (ABS) : ' ZMIN 'METRES';
  191. MESS ' SOIT UN ECART DE : ' ERG2 '%' ;
  192. *
  193. *SAUT 1 LIGNE;
  194. *MESS ' ABSICE DU POINT OU SE TROUVE LE MAXIMUM : ' ;
  195. *MESS ' - ABSICE THEORIQUE : ' XX1 'METRES';
  196. *MESS ' - ABSICE CALCULEE : ' XX 'METRES';
  197. *MESS ' SOIT UN ECART DE : ' ERG3 '%' ;
  198. *
  199. SAUT 1 LIGNE;
  200. MESS ' Nombre de noeuds : ' NOSURF ;
  201. MESS ' Nombre d elements : ' ELSURF ;
  202. MESS ' Soit ' (NOSURF * 6) 'd.d.l.';
  203. SAUT 1 LIGNE;
  204. *
  205. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  206. * CODE FONCTIONNEMENT *
  207. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  208. *
  209. ERGMAX = MAXI (PROG ERG1 ERG2);
  210. *
  211. * Il y a erreur si un des pourcentages ERG* est suprieur 5%
  212. *
  213. SI (ERGMAX <EG 5.);
  214. ERRE 0 ;
  215. SINO;
  216. ERRE 5 ;
  217. FINSI;
  218. *
  219. TEMPS;
  220.  
  221. FIN;
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  

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