Télécharger test_pres_cham.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test1_pres_cham.dgibi
  2. ************************************************************************
  3. * Utilisation de l'operateur PRES avec un MCHAML *
  4. * *
  5. * Une comparaison est faite entre les utilisation dans PRES *
  6. * - d'un CHPOINT (CHAR1) *
  7. * - d'un MCHAML (CHAR2) *
  8. * *
  9. * Calcul mecanique avec d'une poutre avec section en T : *
  10. * - encastrement de la face gauche *
  11. * - pression variable repartie sur la face superieure *
  12. ************************************************************************
  13. 'OPTI' 'DIME' 3 'ELEM' 'CUB8' 'TRAC' 'X' 'ECHO' 0 ;
  14.  
  15. * Indicateur pour traces
  16. ITRAC = FAUX ;
  17. *
  18. *-----------------------------------------------------------------------*
  19. * Maillage *
  20. *-----------------------------------------------------------------------*
  21. * Dimensions section
  22. LA = 0.004 ; LB = 0.02 ; H1 = 0.03 ;
  23. * Longueur de la poutre
  24. LONG = 0.2 ;
  25. * Taille de maille cible
  26. OPTI 'DENS' LA ;
  27. *
  28. P0 = 0. 0. 0. ;
  29. P1 = P0 PLUS (0. LA 0.) ;
  30. P2 = P1 PLUS (0. 0. (H1 - (2*LA))) ;
  31. P3 = P2 PLUS (0. (LB - LA) 0.) ;
  32. P4 = P3 PLUS (0. 0. (2*LA)) ;
  33. P5 = P0 PLUS (0. 0. H1) ;
  34. LIG1 = D P0 P1 ;
  35. LIG1 = D LIG1 P2 ;
  36. LIG1 = D LIG1 P3 ;
  37. LIG1 = D LIG1 P4 ;
  38. LIG1 = D LIG1 P5 ;
  39. LIG1 = D LIG1 P0 ;
  40. SUR1 = SURF LIG1 'PLAN' ;
  41. SUR1 = SUR1 ET (SUR1 SYME 'DROI' P0 P5) ;
  42. ELIM SUR1 1.E-9 ;
  43. *
  44. VOL1 = SUR1 VOLU 'TRAN' 50 (LONG 0. 0.) ;
  45. 'SI' ITRAC ;
  46. TRAC 'CACH' VOL1 'TITR' 'Maillage' ;
  47. 'FINSI' ;
  48. *
  49. ENV1 = ENVE VOL1 ;
  50. PHAUT = ENV1 POIN 'PLAN' P4 P5 (P4 PLUS (1. 0. 0.)) ;
  51. SHAUT = ENV1 ELEM 'APPU' 'STRI' PHAUT ;
  52. *
  53. *-----------------------------------------------------------------------*
  54. * Modeles *
  55. *-----------------------------------------------------------------------*
  56. * Modele mecanique
  57. MOMEC = 'MODE' VOL1 'MECANIQUE' 'ELASTIQUE' ;
  58. MAMEC = 'MATE' MOMEC 'YOUN' 200.E9 'NU' 0.3 ;
  59. *
  60. *-----------------------------------------------------------------------*
  61. * Conditions aux limites *
  62. *-----------------------------------------------------------------------*
  63. * Encastrement
  64. BL1 = BLOQ 'DEPL' SUR1 ;
  65. *
  66. *-----------------------------------------------------------------------*
  67. * Chargements *
  68. *-----------------------------------------------------------------------*
  69. * Pression repartie sur la face superieure
  70. PMAX = 2.E6 ;
  71. CHX = 'COOR' SHAUT 1 ;
  72. CHP = PMAX * ((CHX - ('MINI' CHX)) / (('MAXI' CHX) - ('MINI' CHX))) ;
  73. LTPS = 'PROG' 0. 5. 10. ; LCOE = 'PROG' 0. 1. 2. ;
  74. EV1 = EVOL 'MANU' 'TEMP' LTPS 'COEF' LCOE ;
  75. *
  76. * Chargement PRES avec un CHPOINT
  77. F1 = 'PRES' 'MASS' MOMEC CHP ;
  78. CHAR1 = 'CHAR' 'MECA' F1 EV1 ;
  79. *
  80. * Chargement PRES avec un MCHAML
  81. F2 = 'PRES' 'MASS' MOMEC ('CHAN' 'CHAM' CHP SHAUT) ;
  82. CHAR2 = 'CHAR' 'MECA' F2 EV1 ;
  83. *
  84. * Trace
  85. V1 = 'VECT' F1 'FORC' 'ROUG' ;
  86. 'SI' ITRAC ;
  87. TITR_Z = 'CHAI' 'Pression repartie' ;
  88. TRAC 'CACH' V1 (SHAUT ET ('ARETE' ENV1)) 'TITR' TITR_Z ;
  89. 'FINSI' ;
  90. *
  91. *-----------------------------------------------------------------------*
  92. * Resolution *
  93. *-----------------------------------------------------------------------*
  94. * Resolution avec PASAPAS
  95. TAB1 = TABL ;
  96. TAB1.'MODELE' = MOMEC ;
  97. TAB1.'CARACTERISTIQUES' = MAMEC ;
  98. TAB1.'BLOCAGES_MECANIQUES' = BL1 ;
  99. TAB1.'CHARGEMENT' = CHAR1 ;
  100. TAB1.'TEMPS_CALCULES' = LTPS ;
  101. PASAPAS TAB1 ;
  102. *
  103. TAB2 = TABL ;
  104. TAB2.'MODELE' = MOMEC ;
  105. TAB2.'CARACTERISTIQUES' = MAMEC ;
  106. TAB2.'BLOCAGES_MECANIQUES' = BL1 ;
  107. TAB2.'CHARGEMENT' = CHAR2 ;
  108. TAB2.'TEMPS_CALCULES' = LTPS ;
  109. PASAPAS TAB2 ;
  110. *
  111. *-----------------------------------------------------------------------*
  112. * Visualisation *
  113. *-----------------------------------------------------------------------*
  114. * Deformee de la poutre
  115. DEP1 = TAB1.'DEPLACEMENTS' ;DEP2 = TAB2.'DEPLACEMENTS' ;
  116. DEF_T0 = 'DEFO' VOL1 DEP1. 0 1 VERT ;
  117. 'REPE' BDEP ('DIME' DEP1) ;
  118. DEF1_TI = 'DEFO' VOL1 DEP1. (&BDEP - 1) 10. ROUGE ;
  119. DEF2_TI = 'DEFO' VOL1 DEP2. (&BDEP - 1) 10. BLEU ;
  120. 'SI' ITRAC ;
  121. TITR_Z = 'CHAI' 'Instant ' TAB1.'TEMPS'. (&BDEP - 1) ;
  122. TRAC (DEF_T0 'ET' DEF1_TI 'ET' DEF2_TI) 'TITR' TITR_Z 'CACH' ;
  123. 'FINSI' ;
  124. 'FIN' BDEP ;
  125. *
  126. *-----------------------------------------------------------------------*
  127. * Comparaison des resultats *
  128. *-----------------------------------------------------------------------*
  129. PT = VOL1 POIN PROC (P5 PLUS (LONG 0. 0.)) ;
  130. DEP1_PT = 'EXTR' (TAB1.'DEPLACEMENTS'. (('DIME' DEP1) - 1)) 'UZ' PT ;
  131. DEP2_PT = 'EXTR' (TAB2.'DEPLACEMENTS'. (('DIME' DEP2) - 1)) 'UZ' PT ;
  132. ERRMAX = 'ABS' (DEP1_PT - DEP2_PT) ;
  133. 'SI' (ERRMAX '>' 1.E-3) ;
  134. 'ERRE' 5 ;
  135. 'FINSI' ;
  136. *
  137. 'FIN' ;
  138.  
  139.  
  140.  
  141.  
  142.  

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