Télécharger Contact2Djeu.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : Contact2Djeu.dgibi
  2. *---------------------------------------------------------------------*
  3. * Ce cas-test est une adaptation du cas-test Contact2D avec l'ajout*
  4. * d'un jeu entre les deux objets. Il permet de tester la gestion du *
  5. * contact fort 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.9998 ;
  25. S2 = 10. 0.9998 ;
  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. *chjeu1 = 'MANU' 'CHML' MCONT1 'JEU' 1e-4 ;
  39. *list MCONT1;
  40. * *
  41. * Traces : *
  42. * *
  43. 'SI' IG1 ;
  44. 'TITR' 'Maillages' ;
  45. 'TRAC' 'FACE' (LS1 'ET' SM1 'ET' MCONT1) ;
  46. 'FINS' ;
  47. * *
  48. *---------------------------------------------------------------------*
  49. * MODELES / CARACTERISTIQUES *
  50. *---------------------------------------------------------------------*
  51. * *
  52. MODM1 = 'MODE' SM1 'MECANIQUE' 'ELASTIQUE' ;
  53. MATM1 = 'MATE' MODM1 'YOUN' 1.E3 'NU' 0.3 ;
  54. * *
  55. *---------------------------------------------------------------------*
  56. * C.L. / CHARGEMENT *
  57. *---------------------------------------------------------------------*
  58. * *
  59. * Deplacements imposes : *
  60. * *
  61. LM3 = SM1 'COTE' 3 ;
  62. LM4 = SM1 'COTE' 4 ;
  63. CLLM3 = 'BLOQ' LM3 'UY' ;
  64. CLLM4 = 'BLOQ' LM4 'UX' ;
  65. CLS1 = 'BLOQ' LS1 'DEPL' ;
  66. CL0 = CLLM3 'ET' CLLM4 'ET' CLS1 ;
  67. * *
  68. UY0 = -0.0002 ;
  69. DCLLM3 = 'DEPI' CLLM3 UY0 ;
  70. * *
  71. 'SI' IG1 ;
  72. 'TITR' 'Deplacement impose au bord superieur du carre.' ;
  73. 'TRAC' ('VECT' (DCLLM3 'NOMC' 'UY') 1. 'UX' 'UY' 'VERT')
  74. (LS1 'ET' SM1) ;
  75. 'FINS' ;
  76. * *
  77. * Chargements : *
  78. * *
  79. LTPS1 = 'PROG' 0. 1. ;
  80. EV1 = 'EVOL' 'MANU' 'TEMPS' LTPS1 ('PROG' 0. 1.) ;
  81. * *
  82. CHARU1 = 'CHAR' 'DIMP' DCLLM3 EV1 ;
  83. CHAR0 = CHARU1 ;
  84. MODCONTA= model lm1 contact unilateral (inve ls1);
  85. chjeu1 = 'MANU' 'CHML' MODCONTA 'JEU' 1e-4 ;
  86. list modconta;
  87. MATCONTA= 'MATE' MODCONTA 'JEU' chjeu1;
  88. * *
  89. *---------------------------------------------------------------------*
  90. * RESOLUTION *
  91. *---------------------------------------------------------------------*
  92. * *
  93. * Construction de la table PASAPAS : *
  94. * *
  95. TAB1 = 'TABL' ;
  96. TAB1 . 'TEMPS_CALCULES' = LTPS1 ;
  97. TAB1 . 'MODELE' = MODM1 et modconta;
  98. TAB1 . 'CARACTERISTIQUES' = MATM1 'ET' MATCONTA ;
  99. TAB1 . 'BLOCAGES_MECANIQUES' = CL0 ;
  100. TAB1 . 'CHARGEMENT' = CHAR0 ;
  101. *TAB1 . 'CONTACT' = MCONT1 ;
  102. *TAB1 . 'GRANDS_DEPLACEMENTS' = FAUX ;
  103. * *
  104. * Resolution : *
  105. * *
  106. TAB2 = PASAPAS TAB1 ;
  107. * *
  108. *---------------------------------------------------------------------*
  109. * DEPOUILLEMENT *
  110. *---------------------------------------------------------------------*
  111. * *
  112. DEP1 = (TAB2 . 'DEPLACEMENTS' . 1) 'ENLE' 'LX' ;
  113. * *
  114. * Deformee : *
  115. * *
  116. DEFO0 = 'DEFO' (SM1 'ET' LS1) DEP1 0. 'VERT' ;
  117. DEFO1 = 'DEFO' (SM1 'ET' LS1) DEP1 1. 'ROUG' ;
  118. * *
  119. 'SI' IG1 ;
  120. 'TITR' 'Maillages non deforme (vert) et deforme (rouge).' ;
  121. 'TRAC' (DEFO0 'ET' DEFO1) ;
  122. 'FINS' ;
  123. * *
  124. * Definition des deplacements solutions et comparaison avec la *
  125. * solution EF : *
  126. * *
  127. EPXX1 = ((-1. * UY0) - 0.0001) * 0.1 * 0.3;
  128. UXSM1 = (('COOR' 1 SM1) + 5.) * EPXX1 ;
  129. SOLUX1 = UXSM1 'NOMC' 'UX' ;
  130. EPYY1 = ((-1. * UY0) - 0.0001) * 0.1 ;
  131. UYSM1 = ((('COOR' 2 SM1) - 1.) * EPYY1) + 0.0001;
  132. SOLUY1 = (-1. * UYSM1) 'NOMC' 'UY' ;
  133. SOLU1 = SOLUX1 'ET' SOLUY1 ;
  134. ERR1 = 'MAXI' ('ABS' ((SOLU1 - DEP1) / ('MAXI' 'ABS' (SOLU1)))) ;
  135. * *
  136. 'SI' IG1 ;
  137. 'TITR' 'Champ de deplacements.' ;
  138. 'TRAC' ('CHAN' 'CHAM' DEP1 MODM1 'NOEUDS') MODM1 ;
  139. 'FINS' ;
  140. * *
  141. * Comparaison des champs de contraintes : *
  142. * *
  143. EPS1 = 'EPSI' MODM1 SOLU1 ;
  144. SIG1 = 'ELAS' MODM1 EPS1 MATM1 ;
  145. SIG2 = TAB2 . 'CONTRAINTES' . 1 ;
  146. ERR2 = 'MAXI' ('ABS' ((SIG1 - SIG2) / ('MAXI' 'ABS' (SIG1)))) ;
  147. * *
  148. 'SI' IG1 ;
  149. 'TITR' 'Champ de contraintes.' ;
  150. 'TRAC' SIG2 MODM1 ;
  151. 'FINS' ;
  152. * *
  153. * Visualisations des reactions : *
  154. * *
  155. 'SI' IG1 ;
  156. REAC1 = TAB2 . 'REACTIONS' . 1 ;
  157. VR1 = 'VECT' REAC1 0.8E-2 'FX' 'FY' 'ROUG' ;
  158. 'TITR' 'Forces de reaction.' ;
  159. 'TRAC' VR1 (LS1 'ET' SM1) ;
  160. 'FINS' ;
  161. * *
  162. LERR0 = 'PROG' ERR1 ERR2 ;
  163. ERR0 = 'MAXI' LERR0 ;
  164. ERRMAX1 = 1.E-3 ;
  165. * *
  166. 'OPTI' 'ECHO' 0 ;
  167. 'SAUT' 1 'LIGN' ;
  168. '------------------------ RESULTAT CAS-TEST ------------------------' ;
  169. 'SAUT' 1 'LIGN' ;
  170. 'Ecart relatif a la solution calculee sur les deplacements' ;
  171. 'et les contraintes :' ;
  172. 'MESS'
  173. '-----------------------------------------------------------' ;
  174. 'MESS' ;
  175. 'MESS' ' MAX. ERREUR =' ERR0 ;
  176. 'SAUT' 1 'LIGN' ;
  177. 'SI' (ERR0 '<EG' ERRMAX1) ;
  178. 'MESS' '==> Erreur relative inferieure a' ERRMAX1 ':' ;
  179. 'MESS' '' ;
  180. 'MESS' ' __________________' ;
  181. 'MESS' ' | |' ;
  182. 'MESS' ' | TEST REUSSI ! |' ;
  183. 'MESS' ' |________________|' ;
  184. 'SAUT' 1 'LIGN' ;
  185. 'SINO' ;
  186. 'MESS' '==> Erreur relative superrieure a' ERRMAX1 ':' ;
  187. 'MESS' '' ;
  188. 'MESS' ' __________________' ;
  189. 'MESS' ' | |' ;
  190. 'MESS' ' | ERREUR ! |' ;
  191. 'MESS' ' |________________|' ;
  192. 'SAUT' 1 'LIGN' ;
  193. 'ERRE' 5 ;
  194. 'FINS' ;
  195. '-------------------------- FIN CAS-TEST ---------------------------' ;
  196. * *
  197. 'FIN' ;
  198. 'OPTI' 'ECHO' 1 ;
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  

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