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