Télécharger la plaquette Cast3M


fr en

Télécharger stru3.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  stru3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 1;
  5. GRAPH = 'N';
  6. * 
  7. *********************************************************************
  8. *  TEST : SSLS34/91 DE LA COMMISSION "STRUCTURES ASSEMBLEES"        *
  9. *  ----                                                             *
  10. *                                                                   *
  11. *  DESCRIPTION : PLAQUE ORTHOTROPE RAIDIE SUR APPUIS SIMPLES,       *
  12. *                SOUS UNE CHARGE UNIFORMEMENT REPARTIE.             *
  13. *                Elements DKT                                       *
  14. *                                                                   *
  15. *                                                                   *
  16. *  REF. : JL. BATOZ,J.CHATELAIN,JP. JAMEUX                          *
  17. *         ETUDE D'UNE PLAQUE RAIDIE - BULLETIN DU CLUB S.A.P -      *
  18. *         VOLUME VIII, N1.                                          *
  19. *  Comparaison avec les codes CA.ST.OR, ASELF, ELFINI, PERMAS       *
  20. *  Test envoye pour evolution le 17/09/92                           *
  21. *                                                                   *
  22. *********************************************************************
  23. OPTI ECHO 1;
  25. OPTI ELEM tri3;
  26. * 
  27. *  QUELQUES DONNEES
  28. *
  29. OEIL = 10000 -10000 10000;
  30. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  31. *                             GEOMETRIE                              * 
  32. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  33. * 
  34. *            PLAQUE SUPERIEURE 
  35. *
  36. PA = 0 .5 0; PB =1. .5 0; PD = 1. 0 0; PC = 0 0 0;
  37. PF = .6667 0 0; PE = .3333 0 0;
  38. *
  39. DAB = PA DROIT 36 PB; DBD = PB DROIT 20 PD;
  40. DDC = PD DROIT 36 PC; DCA = PC DROIT 20 PA;
  41. SU1 = DALLER DAB DBD DDC DCA PLAN;
  42. *
  43. *            RAIDISSEUR CENTRAL
  44. * 
  45. PA2 = 0 .5 -.1; PC2 = 0 0 -.1;
  46. D1 = PA DROIT 20 PC;
  47. D3 = PC2 DROIT 20 PA2;
  48. D2 = PC DROIT 4 PC2;
  49. D4 = PA2 DROIT 4 PA;
  50. RA1 = DALL D1 D2 D3 D4 PLAN;
  51. *                                                                               
  52. * SEGMENTS DE DROITES POUR CONDITIONS DE SYMETRIE                               
  53. *
  54. D5 = D2 PLUS (.33333 0. 0.);
  55. D6 = D2 PLUS (.66667 0. 0.);
  56. *                                                                               
  57. *            RAIDISSEURS 2 ET 3                                                 
  58. *                                                                               
  59. RA2 = RA1 PLUS (.33333 0 0);
  60. RA3 = RA1 PLUS (.66667 0 0);
  61. SURF = SU1 ET RA1 ET RA2 ET RA3;
  62. ELIM (SURF ET D1 ET D2 ET D3 ET D4 ET D5 ET D6) .001;
  63. *                                                                               
  64. SI (NEG GRAPH 'N');
  65.   TRAC OEIL SURF;
  66. FINSI;
  67. *                                                                               
  68. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*               
  69. *              MODELE
  70. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*               
  71. *
  72. MODL1 = MODE RA1 MECANIQUE ELASTIQUE ORTHOTROPE dkt;
  73. MODL2 = MODE RA2 MECANIQUE ELASTIQUE ORTHOTROPE dkt;
  74. MODL3 = MODE RA3 MECANIQUE ELASTIQUE ORTHOTROPE dkt;
  75. MODL4 = MODE SU1 MECANIQUE ELASTIQUE ORTHOTROPE dkt;
  76. MODL5 = MODL1 ET MODL2 ET MODL3 ET MODL4;
  77.  
  78. *                                             
  79. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  80. *          MATERIAUX - CARACTERISTIQUES                         *
  81. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  82. *                                                                
  83. DIR1 = 0. 1. 0.;
  84. *                                                                               
  85. MAT1 = MATE MODL1 DIRECTION DIR1 PARALLELE YG1 2.E11 YG2 0.74E11 NU12 .3 G12 .222E11 EPAI 0.005;
  86. *                                                                               
  87. MAT2 = MATE MODL2 DIRECTION DIR1 PARALLELE YG1 2.E11 YG2 0.74E11 NU12 .3 G12 .222E11 EPAI 0.01 ;
  88. *                                                                               
  89. MAT3 = MATE MODL3 DIRECTION DIR1 PARALLELE YG1 2.E11 YG2 0.74E11 NU12 .3 G12 .222E11 EPAI 0.01 ;
  90. *                                                                               
  91. MAT4 = MATE MODL4 DIRECTION DIR1 PARALLELE YG1 2.E11 YG2 0.74E11 NU12 .3 G12 .222E11 EPAI 0.006 ;
  92.  
  93. MAT = MAT1 ET MAT2 ET MAT3 ET MAT4;
  94. *                                                                               
  95.  
  96. *                                                                               
  97. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*           
  98. *             CONDITIONS AUX LIMITES                                *           
  99. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*           
  100. *                                                                               
  101. CL1 = BLOQUE DAB UZ;
  102. CL2 = BLOQUE DAB RY;
  103. CL3 = BLOQUE DBD UZ;
  104. CL4 = BLOQUE DBD RX;
  105. CL5 = SYMT DEPL ROTA (0 0 0) (0 1 0) (0 0 1) SURF 0.0001;
  106. CL6 = SYMT DEPL ROTA (0 0 0) (1 0 0) (0 0 1) SURF 0.0001;
  107. CL = CL1 ET CL2 ET CL3 ET CL4 ET CL5 ET CL6;
  108. *                                                                               
  109. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  110. *                    CHARGEMENT                                     *
  111. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  112. *                                                                               
  113. MOP = 'MODE' SU1 'CHARGEMENT' 'PRESSION' 'DKT' ;
  114. MAP = 'PRES' MOP 'PRES' 5.E4 ;
  115. CH  = 'BSIG' MOP MAP ;
  116. *                                                                               
  117. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  118. *           RESOLUTION - CALCUL ELASTIQUE                           *
  119. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  120. *                                                                               
  121. RI = RIGI MODL5 MAT;
  122. RIG = RI ET CL ;
  123. DEP = RESO RIG CH ;
  124. *                                                                               
  125. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  126. * RESULTATS                                                         *
  127. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
  128. *                                                                               
  129. NOSURF = NBNO SURF ;
  130. ELSURF = NBEL SURF ;
  131. *                                                                               
  132. SI (NEG GRAPH 'N');
  133.   TITR ' PLAQUE RAIDIE : DEFORMEE ' ;
  134.   DEF0 = DEFO SURF DEP 0. ;
  135.   DEF1 = DEFO SURF DEP 25. ROUGE;
  136.   TRACE OEIL (DEF0 ET DEF1) ;
  137.   OEIL =0.25 -1.E6 0.;
  138.   TRACE OEIL (DEF0 ET DEF1) ;
  139. FINSI ;
  140. *                                                                               
  141. UZC = EXTR DEP UZ PC ;
  142. UZC1 = -.603e-3;
  143. ERG1 = ABS ((UZC - UZC1) / UZC1) * 100 ;
  144. *                                                                               
  145. DZ = EXCO DEP UZ;
  146. ZMIN = MINI DZ;
  147. ZMIN1 = -.319e-2;
  148. ERG2 = ABS ((ZMIN - ZMIN1) / ZMIN1) * 100 ;
  149.  
  150.  
  151. *                                                                               
  152. X1  = DZ POIN 'MINI';
  153. XX  = MAXI ( COOR 1 X1);
  154. XX1 = 0.833;
  155. ERG3 = ABS ((XX - XX1) / XX1) * 100 ;
  156.  
  157. SAUT PAGE ;
  159. OPTI ECHO 1;
  160. *                                                                               
  161. MESS ' TEST : SSLS34/90 (FICHIER stru3.DGIBI)';
  162. MESS ' ----';
  163. MESS '  PLAQUE ORTHOTROPE RAIDIE                                                 CHARGEE UNIFORMEMENT SUR APPUIS SIMPLES '                                        *                                                                                SAUT 2 LIGN ;
  164. MESS ' RESULTATS ';
  165. MESS ' --------- ';
  166. *                                                                               
  168. MESS ' DEPLACEMENT THEORIQUE AU POINT PC : ' UZC1 'METRES';
  169. MESS ' DEPLACEMENT CALCULE AU POINT PC   : ' UZC  'METRES' ;
  170. MESS '     SOIT UN ECART DE              : ' ERG1 '%' ;
  171. *                                                                               
  173. MESS ' DEPLACEMENT THEORIQUE MAXIMUM (ABS) : ' ZMIN1 'METRES';
  174. MESS ' DEPLACEMENT CALCULE MAXIMUM (ABS)   : ' ZMIN 'METRES';
  175. MESS '     SOIT UN ECART DE                : ' ERG2 '%' ;
  176. *                                                                               
  177. SAUT 1 LIGNE;
  178. MESS ' ABSICE DU POINT OU SE TROUVE LE MAXIMUM : ' ;
  179. MESS '         - ABSICE THEORIQUE : ' XX1 'METRES';
  180. MESS '         - ABSICE CALCULEE  : ' XX  'METRES';
  181. MESS '     SOIT UN ECART DE       : '  ERG3  '%' ;
  182. *                                                                               
  183. SAUT 1 LIGNE;
  184. MESS ' NOMBRE DE NOEUDS  : ' NOSURF ;
  185. MESS ' NOMBRE D ELEMENTS : ' ELSURF ;
  186. *
  187. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  188. *            CODE FONCTIONNEMENT                                     *          
  189. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  190. *                                                                               
  191. ERGMAX = MAXI (PROG ERG1 ERG2);
  192. *                                                                               
  193. * Il y a erreur si un des pourcentages ERG* est suprieur  5%
  194. *
  196. *                                                                               
  197. SI (ERGMAX <EG 5.);
  198.     ERRE  0 ;
  199. SINON;
  200.     ERRE  5 ;
  201. FINSI;
  202. *
  204. TEMPS;
  206. *
  207. FIN;
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  

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