Télécharger test-coller1.dgibi

Retour à la liste

Numérotation des lignes :

  1. *---------------------------------------------------------------------*
  2. * CAS TEST DE VERIFICATION DE COLLER1
  3. *
  4. OPTI DIME 3 ELEM QUA4 ;
  5.  
  6. *----------------------------- MAILLAGE ------------------------------*
  7. * Maillage poutre
  8. P1 = 0. 0. 0. ;
  9. PZ1 = 0. 0. 1. ;
  10. POUT1 = DROI 1 P1 PZ1 ;
  11.  
  12. * Maillage Coque
  13. VX = 2. ;
  14. PCOQ1 = P1 MOIN (VX 1. 0.) ;
  15. PCOQ2 = P1 PLUS (VX -1 0.) ;
  16. LCOQ1 = DROI 2 PCOQ1 PCOQ2 ;
  17. SCOQ1 = LCOQ1 TRAN 2 (0. 2. 0.) ;
  18. *
  19. P2 = P1 'PLUS' (VX 0. 0.) ;
  20. P2 = SCOQ1 'POIN' 'PROC' P2 ;
  21. V0 = P2 'MOINS' P1 ;
  22. N0 = 'NORME' V0 ;
  23. *
  24. ELIM POUT1 SCOQ1 1.E-9 ;
  25. REP1 = @REPERE P1 (PROG 3 * 0.2) VRAI TURQ ;
  26. *
  27. MTOT = POUT1 et SCOQ1 ;
  28. *
  29. *---------------------- CONDITIONS AUX LIMITES -----------------------*
  30. * Conditions aux limites
  31. BL1 = BLOQ UX UY UZ RX RY PZ1 ;
  32. BL2 = BLOQ RZ PZ1 ;
  33. BL3 = BLOQ UX UY P1 ;
  34. COL1 = COLLER1 SCOQ1 POUT1 ;
  35. BLOT = BL1 ET BL2 ET BL3 ET COL1 ;
  36. *
  37. *---------------------------- CHARGEMENT -----------------------------*
  38. * Chargement
  39. ANG_Z = 15. ;
  40. ROT1 = DEPI BL2 (ANG_Z*pi/180.) ;
  41. *
  42. SCOQ2 = SCOQ1 'TOUR' ANG_Z P1 PZ1 ;
  43. *
  44. *------------------------------ MODELES -------------------------------*
  45. * Modele
  46. e = 10. ;
  47. SECT5 = e*e ;
  48. INRY5 = SECT5*SECT5/12. ;
  49. INRZ5 = INRY5;
  50. TORS5 = INRY5 + INRZ5;
  51. Y_LOCAL = 0. 0. 0.;
  52.  
  53. MOP = MODE POUT1 MECANIQUE ELASTIQUE POUT ;
  54. MAP = MATE MOP YOUN 200.E9 NU 0.3 ;
  55. MAP = MAP ET (CARA MOP 'INRY' INRY5 'INRZ' INRZ5 'SECT' SECT5
  56. 'TORS' TORS5 'VECT' Y_LOCAL);
  57.  
  58. MOC = MODE SCOQ1 MECANIQUE ELASTIQUE COQ4 ;
  59. MAC = MATE MOC YOUN 200.E9 NU 0.3 EPAI 1.;
  60.  
  61. MODTOT = MOC ET MOP ;
  62. MATTOT = MAC ET MAP ;
  63. RIGTOT = RIGI MODTOT MATTOT ;
  64. BLOTOT = RIGTOT ET BLOT ;
  65. *
  66. *---------------------------- RESOLUTION -----------------------------*
  67. *
  68. NBF = 10 ;
  69. L_TPS = PROG 0. PAS (ANG_Z/NBF) ANG_Z ;
  70. N_PAS = DIME L_TPS ;
  71. L_ANG = PROG 0. PAS (1./(N_PAS-1)) 1. ;
  72. EVOLC = EVOL MANU L_TPS L_ANG ;
  73. *
  74. TAB1 = TABLE ;
  75. TAB1 . 'MODELE' = MODTOT ;
  76. TAB1 . 'CARACTERISTIQUES' = MATTOT ;
  77. TAB1 . 'BLOCAGES_MECANIQUES' = BLOT ;
  78. TAB1 . 'CHARGEMENT' = 'CHAR' 'DIMP' ROT1 EVOLC ;
  79. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  80. TAB1 . 'TEMPS_CALCULES' = L_TPS ;
  81. *
  82. OPTI EPSI 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.  

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