Télécharger la plaquette Cast3M


fr en

Télécharger stru2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  stru2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Stru2.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. * Pour visualiser les traces, mettre GRAPH a 'O' :
  8. GRAPH = 'MOT' 'N' ;
  9. OPTION  ECHO 0;
  11. *                                                                               
  12. MESS ' **********************************************';                 
  13. MESS ' *                                            *';                  
  14. MESS ' * COMMISSION V.P.C.S.   TEST STRU2.DGIBI     *';
  15. MESS ' *                       TEST SSLS 33/91      *';
  16. MESS ' * GROUPE : STRUCTURE ASSEMBLEE               *';                  
  17. MESS ' *          STATIQUE LINEAIRE                 *';                  
  18. MESS ' *                                            *';                  
  19. MESS ' * TEST : PLAQUE CARREE ORTHOTROPE            *';                  
  20. MESS ' *        ENCASTREE                           *';                  
  21. MESS ' *                                            *';                  
  22. MESS ' * REMARQUE :  CAS TEST FAIT AVEC ELEMENTS    *';
  23. MESS ' *             QUA4 EN CONTRAINTES PLANES     *';
  24. MESS ' *                                            *';
  25. MESS ' * LES VALEURS DE REFERENCES SONT LES         *';
  26. MESS ' * MOYENNES DE RESULTATS OBTENUS PAR SIX      *';
  27. MESS ' * CODES DIFFERENTS                           *';
  28. MESS ' *                                            *';
  29. MESS ' * TEST MODIFIE LE 3/07/92                    *';
  30. MESS ' **********************************************';                 
  31. *
  32. OPTION DIME 2 ELEM QUA4;                                                        
  33. OPTION MODE PLAN CONT;
  34. *                                                          
  35. *-----------------------------------------------------*       
  36. *   GEOMETRIE                                         *            
  37. *-----------------------------------------------------*       
  38. *
  39. PA =  0.5  0.5;                                                                 
  40. PB = -0.5  0.5;                                                                 
  41. PC = -0.5 -0.5;                                                                 
  42. PD =  0.5 -0.5;                                                                 
  43. PE =  0.   0.;                                                                  
  44. *                                                                               
  45. D1 = PA DROIT 10 PB;                                                            
  46. D2 = PB DROIT 10 PC;                                                            
  47. D3 = PC DROIT 10 PD;                                                            
  48. D4 = PD DROIT 10 PA;                                                            
  49. *                                                                               
  50. SUR1 = DALLAGE D1 D2 D3 D4 PLAN;                                                
  51. SUR = SUR1 ET PE;                                                               
  52. *                                                                               
  53. ELIM SUR 0.001;                                                                 
  54. *                                                                               
  55. SI (NEG GRAPH 'N');                                                             
  56.   TITRE 'PLAQUE CARREE ORTHOTROPE';                                             
  57.   TRACE SUR ;                                                                   
  58. FINSI ;                                                                         
  59. *                                                                               
  60. NBNOSUR = NBNO SUR1;                                                            
  61. NBELSUR = NBEL SUR1;                                                            
  62. *                                                                               
  63. *-----------------------------------------------------*       
  64. *         MODELE                                      *        
  65. *-----------------------------------------------------*       
  66. *                                                                               
  67. MOD = MODE SUR1 MECANIQUE ELASTIQUE ORTHOTROPE QUA4;                            
  68. *                                                                               
  69. *-----------------------------------------------------*       
  70. *       CARACTERISTIQUES - MATERIAUX                  *       
  71. *-----------------------------------------------------*       
  72. *                                                                               
  73. DIR1 = (COS 30.) (SIN 30.);                                                     
  74. *                                                                               
  75. MAT = MATE MOD DIRE DIR1 YG1 2.5737E10 YG2 7.377E9  
  76. YG3 7.377E9 NU12 0.31 NU23 0.31 NU13 0.31 G12 2.319E9;                                                             
  77. *                                                                               
  78. *-----------------------------------------------------*       
  79. *           CONDITIONS AUX LIMITES  -  RIGIDITES      *        
  80. *-----------------------------------------------------*       
  81. *                                                                               
  82. CL = BLOQUE DEPL ROTA D2;                                                       
  83. RI = RIGI MOD MAT ;                                                             
  84. RIG = RI ET CL;                                                                 
  85. *                                                                               
  86. *-----------------------------------------------------*       
  87. *                      CHARGEMENTS                    *       
  88. *-----------------------------------------------------*       
  89. * IL S'AGIT D'UN CHARGEMENT LINEIQUE DE 1e4 NEWTONS/M 
  90. * SUR AD. L'EPAISSEUR EST DE 0.01M. ON A DONC UNE  
  91. * FORCE PAR UNITE D'EPAISSEUR DE 1E6 N/M
  92.  
  93. FTOT = 1.E6;                                                                    
  94. NBNOD4 = NBNO D4;                                                               
  95. F5 = FTOT / ((2 * NBNOD4) - 2);                                                 
  96. F50 = F5 0.;                                                                    
  97. I = 1;                                                                          
  98. NBI = NBNOD4 - 1;                                                               
  99. CHA = FORC (0. 0.) (D4 POIN 1);                                                 
  100. REPETER BLOC1 NBI;
  101.   CHA1 = FORC F50 (D4 POIN I);                                                  
  102.   CHA2 = FORC F50 (D4 POIN (I+1));                                              
  103.   CHA = CHA ET CHA1 ET CHA2;                                                    
  104.   I = I + 1;                                                                    
  105. FIN BLOC1;                                                                      
  106.  
  107. *-----------------------------------------------------*       
  108. *RESULTATS : DEPLACEMENTS,  CONTRAINTES               *
  109. *-----------------------------------------------------*       
  110. *                                                                               
  111. DE = RESO RIG CHA;                                                              
  112. *                                                                               
  113. SI (NEG GRAPH 'N');                                                             
  114.   TITR ' PLAQUE CARREE ORTHOTROPE : DEFORMEE .';                                
  115.   DEF0 = DEFO SUR  DE 0.;                                                       
  116.   DEF1 = DEFO SUR  DE 10. ROUGE;                                                
  117.   TRAC (DEF0 ET DEF1);                                                          
  118. FINSI;                                                                          
  119. *                                                                               
  120. UA = EXTR DE UX PA; UA1 = 1.098e-4;
  121. ERGUA = 100 * (ABS ((UA1 - UA) / UA1));
  122. *
  123. UD = EXTR DE UX PD; UD1 = 1.003e-4;
  124. ERGUD = 100 * (ABS ((UD1 - UD) / UD1));
  125. *
  126. UE = EXTR DE UX PE; UE1 = 4.768E-5;
  127. ERGUE = 100 * (ABS ((UE1 - UE) / UE1));
  128. *
  129. VA = EXTR DE UY PA; VA1 = -1.309E-4;
  130. ERGVA = 100 * (ABS ((VA1 - VA) / VA1));
  131. *
  132. VD = EXTR DE UY PD; VD1 = -7.134E-5;
  133. ERGVD = 100 * (ABS ((VD1 - VD) / VD1));
  134. *
  135. VE = EXTR DE UY PE; VE1 = -4.965E-5;
  136. ERGVE = 100 * (ABS ((VE1 - VE) / VE1));
  137. *                                                                                                                                                             
  139. MESS '    Resultats : test stru2.dgibi';
  140. MESS '    ----------------------------';                                          
  142. MESS ' Deplacements';
  143. MESS ' ------------';
  145. *
  146. MESS '   Point A :';
  147. MESS '   ---------';
  148. MESS '  Deplacement theorique horizontal :' UA1 'm';
  149. MESS '  Deplacement calcule horizontal   :' UA 'm';
  150. MESS '      Soit un ecart de :' ERGUA '%';
  152. MESS '  Deplacement theorique vertical   :' VA1 'm';
  153. MESS '  Deplacement calcule vertical     :' VA 'm';
  154. MESS '      Soit un ecart de :' ERGVA '%';
  156. *
  157. MESS '   Point D :';
  158. MESS '   ---------';
  159. MESS '  Deplacement theorique horizontal :' UD1 'm';
  160. MESS '  Deplacement calcule horizontal   :' UD 'm';
  161. MESS '      Soit un ecart de :' ERGUD '%';
  163. MESS '  Deplacement theorique vertical   :' VD1 'm';
  164. MESS '  Deplacement calcule vertical     :' VD 'm';
  165. MESS '      Soit un ecart de :' ERGVD '%';
  167. *
  168. MESS '   Point E :';
  169. MESS '   ---------';
  170. MESS '  Deplacement theorique horizontal :' UE1 'm';
  171. MESS '  Deplacement calcule horizontal   :' UE 'm';
  172. MESS '      Soit un ecart de :' ERGUE '%';
  174. MESS '  Deplacement theorique vertical   :' VE1 'm';
  175. MESS '  Deplacement calcule vertical     :' VE 'm';
  176. MESS '      Soit un ecart de :' ERGVE '%';
  178. *
  179. *SAUT 1 LIGN;
  180. *MESS ' Contraintes au point E';
  181. *MESS ' ----------------------';
  182. *SAUT 1 LIGN;
  183. *                                                                               
  184. CHAM1 = SIGM MOD MAT DE;                                                        
  185. CHAM2 = CHAN NOEUD MOD CHAM1;                                                   
  186. CHPO1 = CHAN CHPO MOD CHAM2;                                                    
  187. *                                                                               
  188. SMXXE = EXTR CHPO1 SMXX PE;                                                     
  189. SMXXE = SMXXE / 1.E6;
  190. SMXXE1 = 1.004;
  191. ERGXX = 100 * (ABS ((SMXXE1 - SMXXE) / SMXXE1));                                                          
  192. *
  193. SMYYE = EXTR CHPO1 SMYY PE;                                                     
  194. SMYYE = SMYYE / 1.E6;                                                           
  195. SMXYE = EXTR CHPO1 SMXY PE;                                                     
  196. SMXYE = SMXYE / 1.E6;                                                           
  197. *  
  198. *MESS ' Contrainte theorique horizontale :' SMXXE1 
  199. *'MPa';
  200. *MESS ' Contrainte calculee horizontale  :' SMXXE 
  201. *'MPa';
  202. *MESS '      Soit un ecart de :' ERGXX '%';
  204. *MESS ' Contrainte calculee verticale    :' SMYYE 
  205. *'MPa';
  207. *MESS ' Contrainte calculee de cisaillement :' SMXYE 
  208. *'MPa';
  210. *                                                                               
  211. MESS ' Nombre de noeuds  :' NBNOSUR;                                            
  212. MESS ' Nombre d elements :' NBELSUR;                                            
  213. *                                                                               
  214. * UTILISATION DE LA PROCEDURE AFFICHE                                           
  215. * 
  216. SI (NEG GRAPH 'N');
  217.   AFFICHE RIG CHA SUR 1.E-6 1.E-6; 
  218. FINSI;                                               
  219. *                                                                               
  220. SI (NEG GRAPH 'N');                                                             
  221.   CHAM10 = EXCO SMXX CHAM1;                                                     
  222.   CHAM10 = 'CHAN' 'TYPE' CHAM10 'SCALAIRE' ;
  223.   TITRE 'CONTRAINTE SMXX';                                                      
  224.   TRACE MOD CHAM10 SUR;                                                         
  225. *                                                                               
  226.   CHAM11 = EXCO SMYY CHAM1;                                                     
  227.   CHAM11 = 'CHAN' 'TYPE' CHAM11 'SCALAIRE' ;
  228.   TITRE 'CONTRAINTE SMYY';                                                      
  229.   TRACE MOD CHAM11 SUR;                                                         
  230. *                                                                               
  231.   CHAM12 = EXCO SMXY CHAM1;                                                     
  232.   CHAM12 = 'CHAN' 'TYPE' CHAM12 'SCALAIRE' ;
  233.   TITRE 'CONTRAINTE SMXY';                                                      
  234.   TRACE MOD CHAM12 SUR;                                                         
  235. *                                                                               
  236.   CHAM13 = VMIS MOD CHAM1;                                                      
  237.   CHAM13 = 'CHAN' 'TYPE' CHAM13 'SCALAIRE' ;
  238.   TITRE 'CONTRAINTE DE VON-MISES';                                              
  239.   TRACE MOD CHAM13 SUR;                                                         
  240. FINSI ;                                                                         
  241. *                                                                               
  242. *=====================================================*
  243. *                code fonctionnement                  *
  244. *=====================================================*
  245. *
  246. ERGMAX=MAXI (PROG ERGUA ERGUD ERGUE ERGVA ERGVD ERGVD);
  247. *
  248. * Il y a erreur si un des pourcentages ERG* 
  249. * est superieur a 5%
  250. *
  252. SI (ERGMAX <EG 5.);
  253.     ERRE  0 ;                                                                   
  254. SINO;                                                                           
  255.     ERRE  5 ;                                                                   
  256. FINSI;                                                                          
  257. *
  259. TEMPS;
  261. *
  262. FIN;                                                                            
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  

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