Télécharger Contact2D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : Contact2D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *---------------------------------------------------------------------*
  5. * Ce cas-test permet de tester la gestion du contact par PASAPAS. *
  6. * Il simule la mise en contact, en deplacements imposes, d'un carre *
  7. * sur une surface rigide. Le probleme est traite en 2D, contraintes *
  8. * planes. On impose le deplacement de l'arete sup. du carre. Sa base *
  9. * entre en contact, le carre est mis en compression. On compare la *
  10. * solution EF a la solution analytique. *
  11. *---------------------------------------------------------------------*
  12. * *
  13. **'OPTI' ECHO 0 ;
  14. * *
  15. 'OPTI' 'DIME' 2 'ELEM' QUA4 'MODE' 'PLAN' 'CONT' ;
  16. * *
  17. * Si TRACes desires, mettre IG1 a VRAI : *
  18. IG1 = FAUX ;
  19. * *
  20. *---------------------------------------------------------------------*
  21. * MAILLAGE *
  22. *---------------------------------------------------------------------*
  23. * *
  24. S1 = -10. 0.9999 ;
  25. S2 = 10. 0.9999 ;
  26. NLS1 = 5 ;
  27. LS1 = S1 'DROI' NLS1 S2 ;
  28. * *
  29. M1 = -5. 1. ;
  30. M2 = 5. 1. ;
  31. NLM1 = 7 ;
  32. LM1 = M1 'DROI' NLM1 M2 ;
  33. SM1 = LM1 'TRAN' NLM1 (0. 10.) ;
  34. * *
  35. * Maillages de contact : *
  36. * *
  37. MCONT1 = ('IMPO' 'MAIL' LM1 ('INVE' LS1)) 'COUL' 'JAUN' ;
  38. *list MCONT1;
  39. * *
  40. * Traces : *
  41. * *
  42. 'SI' IG1 ;
  43. 'TITR' 'Maillages' ;
  44. 'TRAC' 'FACE' (LS1 'ET' SM1 'ET' MCONT1) ;
  45. 'FINS' ;
  46. * *
  47. *---------------------------------------------------------------------*
  48. * MODELES / CARACTERISTIQUES *
  49. *---------------------------------------------------------------------*
  50. * *
  51. MODM1 = 'MODE' SM1 'MECANIQUE' 'ELASTIQUE' ;
  52. MATM1 = 'MATE' MODM1 'YOUN' 1.E3 'NU' 0.3 ;
  53. * *
  54. *---------------------------------------------------------------------*
  55. * C.L. / CHARGEMENT *
  56. *---------------------------------------------------------------------*
  57. * *
  58. * Deplacements imposes : *
  59. * *
  60. LM3 = SM1 'COTE' 3 ;
  61. LM4 = SM1 'COTE' 4 ;
  62. CLLM3 = 'BLOQ' LM3 'UY' ;
  63. CLLM4 = 'BLOQ' LM4 'UX' ;
  64. CLS1 = 'BLOQ' LS1 'DEPL' ;
  65. CL0 = CLLM3 'ET' CLLM4 'ET' CLS1 ;
  66. * *
  67. UY0 = -0.0003 ;
  68. DCLLM3 = 'DEPI' CLLM3 UY0 ;
  69. * *
  70. 'SI' IG1 ;
  71. 'TITR' 'Deplacement impose au bord superieur du carre.' ;
  72. 'TRAC' ('VECT' (DCLLM3 'NOMC' 'UY') 1. 'UX' 'UY' 'VERT')
  73. (LS1 'ET' SM1) ;
  74. 'FINS' ;
  75. * *
  76. * Chargements : *
  77. * *
  78. LTPS1 = 'PROG' 0. 1. ;
  79. EV1 = 'EVOL' 'MANU' 'TEMPS' LTPS1 ('PROG' 0. 1.) ;
  80. * *
  81. CHARU1 = 'CHAR' 'DIMP' DCLLM3 EV1 ;
  82. CHAR0 = CHARU1 ;
  83. MODCONTA= model MCONT1 contact unilateral;
  84. list modconta;
  85. * *
  86. *---------------------------------------------------------------------*
  87. * RESOLUTION *
  88. *---------------------------------------------------------------------*
  89. * *
  90. * Construction de la table PASAPAS : *
  91. * *
  92. TAB1 = 'TABL' ;
  93. TAB1 . 'TEMPS_CALCULES' = LTPS1 ;
  94. TAB1 . 'MODELE' = MODM1 et modconta;
  95. TAB1 . 'CARACTERISTIQUES' = MATM1 ;
  96. TAB1 . 'BLOCAGES_MECANIQUES' = CL0 ;
  97. TAB1 . 'CHARGEMENT' = CHAR0 ;
  98. *TAB1 . 'CONTACT' = MCONT1 ;
  99. *TAB1 . 'GRANDS_DEPLACEMENTS' = FAUX ;
  100. * *
  101. * Resolution : *
  102. * *
  103. TAB2 = PASAPAS TAB1 ;
  104. * *
  105. *---------------------------------------------------------------------*
  106. * DEPOUILLEMENT *
  107. *---------------------------------------------------------------------*
  108. * *
  109. DEP1 = (TAB2 . 'DEPLACEMENTS' . 1) 'ENLE' 'LX' ;
  110. * *
  111. * Deformee : *
  112. * *
  113. DEFO0 = 'DEFO' (SM1 'ET' LS1) DEP1 0. 'VERT' ;
  114. DEFO1 = 'DEFO' (SM1 'ET' LS1) DEP1 1. 'ROUG' ;
  115. * *
  116. 'SI' IG1 ;
  117. 'TITR' 'Maillages non deforme (vert) et deforme (rouge).' ;
  118. 'TRAC' (DEFO0 'ET' DEFO1) ;
  119. 'FINS' ;
  120. * *
  121. * Definition des deplacements solutions et comparaison avec la *
  122. * solution EF : *
  123. * *
  124. EPXX1 = ((-1. * UY0) - 0.0001) * 0.1 * 0.3;
  125. UXSM1 = (('COOR' 1 SM1) + 5.) * EPXX1 ;
  126. SOLUX1 = UXSM1 'NOMC' 'UX' ;
  127. EPYY1 = ((-1. * UY0) - 0.0001) * 0.1 ;
  128. UYSM1 = ((('COOR' 2 SM1) - 1.) * EPYY1) + 0.0001;
  129. SOLUY1 = (-1. * UYSM1) 'NOMC' 'UY' ;
  130. SOLU1 = SOLUX1 'ET' SOLUY1 ;
  131. ERR1 = 'MAXI' ('ABS' ((SOLU1 - DEP1) / ('MAXI' 'ABS' (SOLU1)))) ;
  132. * *
  133. 'SI' IG1 ;
  134. 'TITR' 'Champ de deplacements.' ;
  135. 'TRAC' ('CHAN' 'CHAM' DEP1 MODM1 'NOEUDS') MODM1 ;
  136. 'FINS' ;
  137. * *
  138. * Comparaison des champs de contraintes : *
  139. * *
  140. EPS1 = 'EPSI' MODM1 SOLU1 ;
  141. SIG1 = 'ELAS' MODM1 EPS1 MATM1 ;
  142. SIG2 = TAB2 . 'CONTRAINTES' . 1 ;
  143. ERR2 = 'MAXI' ('ABS' ((SIG1 - SIG2) / ('MAXI' 'ABS' (SIG1)))) ;
  144. * *
  145. 'SI' IG1 ;
  146. 'TITR' 'Champ de contraintes.' ;
  147. 'TRAC' SIG2 MODM1 ;
  148. 'FINS' ;
  149. * *
  150. * Visualisations des reactions : *
  151. * *
  152. 'SI' IG1 ;
  153. REAC1 = TAB2 . 'REACTIONS' . 1 ;
  154. VR1 = 'VECT' REAC1 0.8E-2 'FX' 'FY' 'ROUG' ;
  155. 'TITR' 'Forces de reaction.' ;
  156. 'TRAC' VR1 (LS1 'ET' SM1) ;
  157. 'FINS' ;
  158. * *
  159. LERR0 = 'PROG' ERR1 ERR2 ;
  160. ERR0 = 'MAXI' LERR0 ;
  161. ERRMAX1 = 1.E-3 ;
  162. * *
  163. 'OPTI' 'ECHO' 0 ;
  164. 'SAUT' 1 'LIGN' ;
  165. '------------------------ RESULTAT CAS-TEST ------------------------' ;
  166. 'SAUT' 1 'LIGN' ;
  167. 'Ecart relatif a la solution calculee sur les deplacements' ;
  168. 'et les contraintes :' ;
  169. 'MESS'
  170. '-----------------------------------------------------------' ;
  171. 'MESS' ;
  172. 'MESS' ' MAX. ERREUR =' ERR0 ;
  173. 'SAUT' 1 'LIGN' ;
  174. 'SI' (ERR0 '<EG' ERRMAX1) ;
  175. 'MESS' '==> Erreur relative inferieure a' ERRMAX1 ':' ;
  176. 'MESS' '' ;
  177. 'MESS' ' __________________' ;
  178. 'MESS' ' | |' ;
  179. 'MESS' ' | TEST REUSSI ! |' ;
  180. 'MESS' ' |________________|' ;
  181. 'SAUT' 1 'LIGN' ;
  182. 'SINO' ;
  183. 'MESS' '==> Erreur relative superrieure a' ERRMAX1 ':' ;
  184. 'MESS' '' ;
  185. 'MESS' ' __________________' ;
  186. 'MESS' ' | |' ;
  187. 'MESS' ' | ERREUR ! |' ;
  188. 'MESS' ' |________________|' ;
  189. 'SAUT' 1 'LIGN' ;
  190. 'ERRE' 5 ;
  191. 'FINS' ;
  192. '-------------------------- FIN CAS-TEST ---------------------------' ;
  193. * *
  194. 'FIN' ;
  195. 'OPTI' 'ECHO' 1 ;
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  

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