Télécharger test-coller1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test1_coller1.dgibi
  2. *---------------------------------------------------------------------*
  3. * CAS TEST DE VERIFICATION DE COLLER1
  4. *
  5. OPTI DIME 3 ELEM QUA4 ;
  6.  
  7. *----------------------------- MAILLAGE ------------------------------*
  8. * Maillage poutre
  9. P1 = 0. 0. 0. ;
  10. PZ1 = 0. 0. 1. ;
  11. POUT1 = DROI 1 P1 PZ1 ;
  12.  
  13. * Maillage Coque
  14. VX = 2. ;
  15. PCOQ1 = P1 MOIN (VX 1. 0.) ;
  16. PCOQ2 = P1 PLUS (VX -1 0.) ;
  17. LCOQ1 = DROI 2 PCOQ1 PCOQ2 ;
  18. SCOQ1 = LCOQ1 TRAN 2 (0. 2. 0.) ;
  19. *
  20. P2 = P1 'PLUS' (VX 0. 0.) ;
  21. P2 = SCOQ1 'POIN' 'PROC' P2 ;
  22. V0 = P2 'MOINS' P1 ;
  23. N0 = 'NORME' V0 ;
  24. *
  25. ELIM POUT1 SCOQ1 1.E-9 ;
  26. REP1 = @REPERE P1 (PROG 3 * 0.2) VRAI TURQ ;
  27. *
  28. MTOT = POUT1 et SCOQ1 ;
  29. *
  30. *---------------------- CONDITIONS AUX LIMITES -----------------------*
  31. * Conditions aux limites
  32. BL1 = BLOQ UX UY UZ RX RY PZ1 ;
  33. BL2 = BLOQ RZ PZ1 ;
  34. BL3 = BLOQ UX UY P1 ;
  35. COL1 = COLLER1 SCOQ1 POUT1 ;
  36. BLOT = BL1 ET BL2 ET BL3 ET COL1 ;
  37. *
  38. *---------------------------- CHARGEMENT -----------------------------*
  39. * Chargement
  40. ANG_Z = 15. ;
  41. ROT1 = DEPI BL2 (ANG_Z*pi/180.) ;
  42. *
  43. SCOQ2 = SCOQ1 'TOUR' ANG_Z P1 PZ1 ;
  44. *
  45. *------------------------------ MODELES -------------------------------*
  46. * Modele
  47. e = 10. ;
  48. SECT5 = e*e ;
  49. INRY5 = SECT5*SECT5/12. ;
  50. INRZ5 = INRY5;
  51. TORS5 = INRY5 + INRZ5;
  52. Y_LOCAL = 0. 0. 0.;
  53.  
  54. MOP = MODE POUT1 MECANIQUE ELASTIQUE POUT ;
  55. MAP = MATE MOP YOUN 200.E9 NU 0.3 ;
  56. MAP = MAP ET (CARA MOP 'INRY' INRY5 'INRZ' INRZ5 'SECT' SECT5
  57. 'TORS' TORS5 'VECT' Y_LOCAL);
  58.  
  59. MOC = MODE SCOQ1 MECANIQUE ELASTIQUE COQ4 ;
  60. MAC = MATE MOC YOUN 200.E9 NU 0.3 EPAI 1.;
  61.  
  62. MODTOT = MOC ET MOP ;
  63. MATTOT = MAC ET MAP ;
  64. RIGTOT = RIGI MODTOT MATTOT ;
  65. BLOTOT = RIGTOT ET BLOT ;
  66. *
  67. *---------------------------- RESOLUTION -----------------------------*
  68. *
  69. NBF = 10 ;
  70. L_TPS = PROG 0. PAS (ANG_Z/NBF) ANG_Z ;
  71. N_PAS = DIME L_TPS ;
  72. L_ANG = PROG 0. PAS (1./(N_PAS-1)) 1. ;
  73. EVOLC = EVOL MANU L_TPS L_ANG ;
  74. *
  75. TAB1 = TABLE ;
  76. TAB1 . 'MODELE' = MODTOT ;
  77. TAB1 . 'CARACTERISTIQUES' = MATTOT ;
  78. TAB1 . 'BLOCAGES_MECANIQUES' = BLOT ;
  79. TAB1 . 'CHARGEMENT' = 'CHAR' 'DIMP' ROT1 EVOLC ;
  80. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  81. TAB1 . 'TEMPS_CALCULES' = L_TPS ;
  82. TAB1 . 'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  83. PASAPAS TAB1 ;
  84. *
  85. CONF0 = 'FORM' ;
  86. FORM (TAB1.'DEPLACEMENTS'. NBF) ;
  87. V1 = P2 'MOINS' P1 ;
  88. FORM CONF0 ;
  89. N1 = 'NORME' V1 ;
  90. PSV01 = 'PSCA' V0 V1 ;
  91. COS01 = PSV01 / (N0*N1) ;
  92. ANG01 = ACOS COS01 ;
  93. MESS 'ANGLE IMPOSE ' ANG_Z ;
  94. MESS 'ANGLE PASAPAS ' ANG01 ;
  95. ERRANG = (ABS (ANG_Z-ANG01)) / ANG_Z ;
  96. 'SI' ('>' ERRANG 0.01) ;
  97. 'ERRE' 5 ;
  98. 'FINSI' ;
  99.  
  100. FIN ;
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  

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