Télécharger Contact2Djeufaible.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : Contact2Djeufaible.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 faible 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' 'FAIB' 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) 'FAIBLE';
  85. chjeu1 = 'MANU' 'CHML' ModCONTa 'JEU' 1e-4 ;
  86. MATCONTA= 'MATE' MODCONTA 'JEU' chjeu1;
  87. * *
  88. *---------------------------------------------------------------------*
  89. * RESOLUTION *
  90. *---------------------------------------------------------------------*
  91. * *
  92. * Construction de la table PASAPAS : *
  93. * *
  94. TAB1 = 'TABL' ;
  95. TAB1 . 'TEMPS_CALCULES' = LTPS1 ;
  96. TAB1 . 'MODELE' = MODM1 et modconta;
  97. TAB1 . 'CARACTERISTIQUES' = MATM1 'ET' MATCONTA ;
  98. TAB1 . 'BLOCAGES_MECANIQUES' = CL0 ;
  99. TAB1 . 'CHARGEMENT' = CHAR0 ;
  100. *TAB1 . 'CONTACT' = MCONT1 ;
  101. *TAB1 . 'GRANDS_DEPLACEMENTS' = FAUX ;
  102. * *
  103. * Resolution : *
  104. * *
  105. TAB2 = PASAPAS TAB1 ;
  106. * *
  107. *---------------------------------------------------------------------*
  108. * DEPOUILLEMENT *
  109. *---------------------------------------------------------------------*
  110. * *
  111. DEP1 = (TAB2 . 'DEPLACEMENTS' . 1) 'ENLE' 'LX' ;
  112. * *
  113. * Deformee : *
  114. * *
  115. DEFO0 = 'DEFO' (SM1 'ET' LS1) DEP1 0. 'VERT' ;
  116. DEFO1 = 'DEFO' (SM1 'ET' LS1) DEP1 1. 'ROUG' ;
  117. * *
  118. 'SI' IG1 ;
  119. 'TITR' 'Maillages non deforme (vert) et deforme (rouge).' ;
  120. 'TRAC' (DEFO0 'ET' DEFO1) ;
  121. 'FINS' ;
  122. * *
  123. * Definition des deplacements solutions et comparaison avec la *
  124. * solution EF : *
  125. * *
  126. EPXX1 = ((-1. * UY0) - 0.0001) * 0.1 * 0.3;
  127. UXSM1 = (('COOR' 1 SM1) + 5.) * EPXX1 ;
  128. SOLUX1 = UXSM1 'NOMC' 'UX' ;
  129. EPYY1 = ((-1. * UY0) - 0.0001) * 0.1 ;
  130. UYSM1 = ((('COOR' 2 SM1) - 1.) * EPYY1) + 0.0001;
  131. SOLUY1 = (-1. * UYSM1) 'NOMC' 'UY' ;
  132. SOLU1 = SOLUX1 'ET' SOLUY1 ;
  133. ERR1 = 'MAXI' ('ABS' ((SOLU1 - DEP1) / ('MAXI' 'ABS' (SOLU1)))) ;
  134. * *
  135. 'SI' IG1 ;
  136. 'TITR' 'Champ de deplacements.' ;
  137. 'TRAC' ('CHAN' 'CHAM' DEP1 MODM1 'NOEUDS') MODM1 ;
  138. 'FINS' ;
  139. * *
  140. * Comparaison des champs de contraintes : *
  141. * *
  142. EPS1 = 'EPSI' MODM1 SOLU1 ;
  143. SIG1 = 'ELAS' MODM1 EPS1 MATM1 ;
  144. SIG2 = TAB2 . 'CONTRAINTES' . 1 ;
  145. ERR2 = 'MAXI' ('ABS' ((SIG1 - SIG2) / ('MAXI' 'ABS' (SIG1)))) ;
  146. * *
  147. 'SI' IG1 ;
  148. 'TITR' 'Champ de contraintes.' ;
  149. 'TRAC' SIG2 MODM1 ;
  150. 'FINS' ;
  151. * *
  152. * Visualisations des reactions : *
  153. * *
  154. 'SI' IG1 ;
  155. REAC1 = TAB2 . 'REACTIONS' . 1 ;
  156. VR1 = 'VECT' REAC1 0.8E-2 'FX' 'FY' 'ROUG' ;
  157. 'TITR' 'Forces de reaction.' ;
  158. 'TRAC' VR1 (LS1 'ET' SM1) ;
  159. 'FINS' ;
  160. * *
  161. LERR0 = 'PROG' ERR1 ERR2 ;
  162. ERR0 = 'MAXI' LERR0 ;
  163. ERRMAX1 = 1.E-3 ;
  164. * *
  165. 'OPTI' 'ECHO' 0 ;
  166. 'SAUT' 1 'LIGN' ;
  167. '------------------------ RESULTAT CAS-TEST ------------------------' ;
  168. 'SAUT' 1 'LIGN' ;
  169. 'Ecart relatif a la solution calculee sur les deplacements' ;
  170. 'et les contraintes :' ;
  171. 'MESS'
  172. '-----------------------------------------------------------' ;
  173. 'MESS' ;
  174. 'MESS' ' MAX. ERREUR =' ERR0 ;
  175. 'SAUT' 1 'LIGN' ;
  176. 'SI' (ERR0 '<EG' ERRMAX1) ;
  177. 'MESS' '==> Erreur relative inferieure a' ERRMAX1 ':' ;
  178. 'MESS' '' ;
  179. 'MESS' ' __________________' ;
  180. 'MESS' ' | |' ;
  181. 'MESS' ' | TEST REUSSI ! |' ;
  182. 'MESS' ' |________________|' ;
  183. 'SAUT' 1 'LIGN' ;
  184. 'SINO' ;
  185. 'MESS' '==> Erreur relative superrieure a' ERRMAX1 ':' ;
  186. 'MESS' '' ;
  187. 'MESS' ' __________________' ;
  188. 'MESS' ' | |' ;
  189. 'MESS' ' | ERREUR ! |' ;
  190. 'MESS' ' |________________|' ;
  191. 'SAUT' 1 'LIGN' ;
  192. 'ERRE' 5 ;
  193. 'FINS' ;
  194. '-------------------------- FIN CAS-TEST ---------------------------' ;
  195. * *
  196. 'FIN' ;
  197. 'OPTI' 'ECHO' 1 ;
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  

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