Télécharger la plaquette Cast3M


fr en

Télécharger Rccmtest.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  rccmtest.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************************
  5. *          TEST DES ROUTINES INTERNES DE @RCCM                    *
  6. *                                                                 *  
  7. *      CALCUL EN AXISYMETRIQUE  ET EN 3D   SUR UN SEGMENT D APPUI *
  8. *      COMPARAISON DES RESULTATS                                  *  
  9. *******************************************************************
  10. OPTION DIME 2 ELEM QUA8 ;
  11. GRAPH = 'N' ;
  12. OO = 0 0  ; OX = 1 0 ; OY = 0 1 ;
  13. DENS 1. ;
  14. P1 = 10 0 ; P2 = 15 0 ;HT = 50. ;
  15. LIG  = D DINI 3 DFIN 1. OO P1 D P2 ;
  16. TOUT = D P1 P2  TRAN 5 ( 0 HT ) ; 
  17. CCONT = TOUT CONTOUR ;
  18. LIINT = TOUT COTE 4 ;
  19. *-------------------------------------------------------------
  20. OPTION MODE AXIS ;
  21. *-------------------------------------------------------------
  22. MOD1 = MODE TOUT MECANIQUE ELASTIQUE ISOTROPE ;
  23. MAT1 = MATE MOD1 YOUN  2.E5  NU .3 ;
  24. *
  25. FPRE = PRESSION MASS  5. MOD1  LIINT  ;
  26. RIG1 = RIGI MOD1 MAT1 ;
  27. C1 = (TOUT COTE 1 ) ET ( TOUT COTE 3) ;
  28. *  -------  conditions limites -----------------------
  29. CLM = BLOQUER DEPLA  C1  ;
  30. SOL1 = RESO ( RIG1 ET CLM) FPRE ;
  31. DEF0 = DEFO CCONT SOL1 0. BLAN ;
  32. DE1 = MAXI SOL1 AVEC ( MOTS UR) ;
  33. AMP = 1./DE1 ;
  34. DEF1 = DEFO CCONT SOL1 AMP ROUG ;
  35.  
  36. 'SI' ( NEG GRAPH N) ;
  37. TRAC ( DEF0 ET DEF1) ;
  38. 'FINSI' ;
  39. SIG1 = SIGMA SOL1 MOD1 MAT1 ;
  40. *
  41. H = 25. ;
  42. PP1 = P1 PLUS ( 0 H ) ; PP2 = P2 PLUS (0 H ) ;
  43. NSEG = 12 ;
  44. **********  NI TRACES NI ECRITURES DE FICHIERS EN INTERNE *****
  45. GRAPL*LOGIQUE =FAUX ;
  46. ECRI*LOGIQUE  =FAUX ;NOMFICH = 'BLANC' ; 
  47. ETAT = 1 ; ICOU = 1 ; 
  48. TITRE ' RCCMCO2   '  ;  
  49. TABB1 = @RCCMCO2 ETAT  ICOU  PP1 PP2  NSEG SIG1  MOD1
  50.                        .1  GRAPL  ECRI  NOMFICH  (0 0 ) ;  
  51. *--------------------------------------------------------------
  52. *          3D 
  53. *--------------------------------------------------------------
  54. OPTION DIME 3 ELEM CU20 ;
  55. OZ = 0 0 1 ;
  56. VOL1 = TOUT VOLU ROTA  3  20  OO OY ;
  57. ELIM .001 VOL1 ;
  58. VOL1 = REGE VOL1 ;
  59. DEPLACER VOL1 TOURNER -10 OO OY ;
  60.  
  61. MOD2  =  MODE VOL1 MECANIQUE ELASTIQUE  ;
  62. MAT2  =  MATE MOD2 YOUN  2.E5  NU .3 ;
  63. VOENV =  VOL1 ENVEL ;
  64. SUI1  =  VOENV POINT CYLI OO OY P1 .01 ;
  65. SUINT = VOENV ELEM APPU STRIC   SUI1 ; 
  66. FPRE = PRESSION MASS  5. MOD2   SUINT   ;
  67. RIG2 = RIGI MOD2 MAT2 ;
  68. P11 = P1 TOURNER 20  OO OY ;P3 = P1 TOURNET 20  OO OY ;
  69. OP = OO PLUS ( 0  HT 0) ;OP1 = P1 PLUS ( 0  HT 0) ;
  70. OP3 = P3 PLUS ( 0  HT 0) ;
  71. *------------------  conditions limites ---------------------
  72. CL1  = (BLOQUER  DEPLA ( VOL1 POINT PLAN  OO P1  P3   .001 ))
  73.    ET  (BLOQUER  DEPLA ( VOL1 POINT PLAN  OP OP1 OP3   .001)) ;
  74. CL2  = (SYMT DEPLA   OO OP  P1   VOL1 .001 )
  75.   ET   (SYMT DEPLA   OO OP  P3   VOL1 .001 )  ;
  76. SOL2 = RESO ( RIG2 ET CL1 ET CL2 ) FPRE ;
  77. SIG2 = SIGMA SOL2 MOD2 MAT2 ;
  78. *
  79. VPER = 0  0  1 ;
  80. 'SI' ( NEG GRAPH N) ;
  81. TRAC ( VOENV ET ( D 2 PP1 PP2 COUL ROUG) ) ;
  82. 'FINSI' ;
  83.  
  84. TABB2 = @RCCMCO2 ETAT ICOU  PP1 PP2  NSEG SIG2 MOD2 
  85.                     .1  GRAPL  ECRI  NOMFICH  VPER ; 
  86.  
  87.  
  88. * ------------ COMPARAISON DES RESULTATS SIGNIFICATIF -------
  89. *   ------------  ENTRE 2D  ET 3D --------------------
  90. *    ATTENTION SMTT  2D   CORRESPOND A SMYY EN 3D 
  91. OPTION ECHO 0 ;  
  92. III = INDEX     TABB2;
  93. REPETER BAB 3  ;
  94. LRE1 = TABB1.&BAB ;   
  95. LRE2 = TABB2.&BAB ;
  96. EPS1 = ABS(((EXTR LRE1 1 ) - (EXTR LRE2 1))/(EXTR LRE2 1)) ;
  97. EPS2 = ABS(((EXTR LRE1 2 ) - (EXTR LRE2 3))/(EXTR LRE2 3)) ;
  98. EPS3 = ABS(((EXTR LRE1 3 ) - (EXTR LRE2 2))/(EXTR LRE2 2)) ;
  99.  
  100.  MESS  ' ERREURS  ' EPS1 EPS2 EPS3 ;
  101. *
  102. 'SI' ((EPS1 > 0.01) OU (EPS2 > 0.01) OU (EPS3 > 0.01)) ;
  103.  ERRE  5 ;
  104.  'SINON' ;
  105.  ERRE  0 ;
  106. 'FINSI' ;               
  107. 'FIN' BAB ;
  108. 'FIN' ;
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  

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