Télécharger Contact3Djeu.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : Contact3Djeu.dgibi
  2. *----------------------------------------------------------------------*
  3. * Cas-test adapté de Contact3D.dgibi avec présence d'un jeu. test du *
  4. * Contact fort en 3D *
  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 cube *
  7. * sur une sol rigide. On impose le deplacement de la face superieure *
  8. * du cube. Sa base entre en contact, le cube est mis en compression. *
  9. * On compare la solution EF a la solution analytique. *
  10. *----------------------------------------------------------------------*
  11. * *
  12. 'OPTI' 'DIME' 3 'ELEM' CUB8 ;
  13. * *
  14. * Si TRACes desires, mettre IG1 a VRAI : *
  15. IG1 = FAUX ;
  16. * *
  17. *----------------------------------------------------------------------*
  18. * MAILLAGE *
  19. *----------------------------------------------------------------------*
  20. * *
  21. * Maillage "Sol" : *
  22. S1 = -10. -5. 0. ;
  23. S2 = +10. -5. 0. ;
  24. NLS1 = 5 ;
  25. LS1 = S1 'DROI' NLS1 S2 ;
  26. SOL1 = LS1 'TRAN' NLS1 (0. 20. 0.) ;
  27. SOL1 = SOL1 'COUL' 'OLIV' ;
  28. * *
  29. * Maillage "Cube" : *
  30. C1 = -5. 0. 0.101 ;
  31. M2 = +5. 0. 0.101 ;
  32. NLC1 = 7 ;
  33. LC1 = C1 'DROI' NLC1 M2 ;
  34. SC1 = LC1 'TRAN' NLC1 (0. 10. 0.) ;
  35. CUB1 = SC1 'VOLU' NLC1 'TRAN' (0. 0. 10.) ;
  36. CUB1 = CUB1 'COUL' 'ORAN' ;
  37. * *
  38. MTOT1 = (SOL1 'ET' CUB1) 'COUL' 'DEFA' ;
  39. * *
  40. * Maillages de contact orientes : *
  41. SOL2 = 'ORIE' ('CHAN' 'TRI3' SOL1) (0. 0. +1.) ;
  42. SCC1 = 'ORIE' ('CHAN' 'TRI3' SC1 ) (0. 0. -1.) ;
  43. **MCONT1 = ('IMPO' 'MAIL' SOL2 SCC1) 'COUL' 'ROUG' ;
  44. * *
  45. * Traces : *
  46. 'SI' IG1 ;
  47. 'TITR' 'Maillages Cube / Sol ' ;
  48. 'TRAC' 'FACE' (CUB1 ET SOL1) ;
  49. 'TITR' 'Maillages Cube / Sol & Maillage de contact (rouge) ' ;
  50. 'TRAC' 'FACE' (MTOT1 ET MCONT1) ;
  51. 'FINS' ;
  52. * *
  53. *----------------------------------------------------------------------*
  54. * MODELES / CARACTERISTIQUES *
  55. *----------------------------------------------------------------------*
  56. * *
  57. MOD1 = 'MODE' CUB1 'MECANIQUE' 'ELASTIQUE' ;
  58. MAT1 = 'MATE' MOD1 'YOUN' 1.E3 'NU' 0.3 ;
  59. * *
  60. MODC1 = 'MODE' SOL2 'CONTACT' 'UNILATERAL' SCC1;
  61. chjeu1 = 'MANU' 'CHML' MODC1 'JEU' 1e-3 ;
  62. MATC1 = 'MATE' MODC1 'JEU' chjeu1 ;
  63. * *
  64. MOD0 = MOD1 'ET' MODC1 ;
  65. MAT0 = MAT1 'ET' MATC1 ;
  66. * *
  67. *----------------------------------------------------------------------*
  68. * C.L. / CHARGEMENT *
  69. *----------------------------------------------------------------------*
  70. * *
  71. * Deplacements imposes : *
  72. CL1 = 'BLOQ' SOL1 'DEPL' ;
  73. * *
  74. SC2 = 'FACE' CUB1 2 ;
  75. SCX0 = (CUB1 'ENVE') 'ELEM' 'APPU' ((CUB1 'COOR' 1) 'POIN' 'MINI') ;
  76. SCY0 = (CUB1 'ENVE') 'ELEM' 'APPU' ((CUB1 'COOR' 2) 'POIN' 'MINI') ;
  77. CLSC2 = 'BLOQ' SC2 'UZ' ;
  78. CLSCX0 = 'BLOQ' SCX0 'UX' ;
  79. CLSCY0 = 'BLOQ' SCY0 'UY' ;
  80. *
  81. CL0 = CL1 'ET' CLSC2 'ET' CLSCX0 'ET' CLSCY0 ;
  82. * *
  83. UZ0 = -0.101 ;
  84. DCLC2 = 'DEPI' CLSC2 UZ0 ;
  85. * *
  86. 'SI' IG1 ;
  87. 'TRAC' ('VECT' ('EXCO' DCLC2 'FLX' 'UZ') 'DEPL' 'ROUG') ('ARET' MTOT1)
  88. 'TITR' 'Deplacement impose a la surface superieure du cube.' ;
  89. 'FINS' ;
  90. * *
  91. * Chargements : *
  92. LTPS1 = 'PROG' 0. 1. ;
  93. EV1 = 'EVOL' 'MANU' 'TEMPS' LTPS1 LTPS1 ;
  94. CHARU1 = 'CHAR' 'DIMP' DCLC2 EV1 ;
  95. CHAR0 = CHARU1 ;
  96. * *
  97. *----------------------------------------------------------------------*
  98. * RESOLUTION *
  99. *----------------------------------------------------------------------*
  100. * *
  101. * Construction de la table PASAPAS : *
  102. * *
  103. TAB1 = 'TABL' ;
  104. TAB1 . 'MODELE' = MOD0 ;
  105. TAB1 . 'CARACTERISTIQUES' = MAT0 ;
  106. TAB1 . 'BLOCAGES_MECANIQUES' = CL0 ;
  107. TAB1 . 'CHARGEMENT' = CHAR0 ;
  108. TAB1 . 'TEMPS_CALCULES' = LTPS1 ;
  109. TAB1 . 'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  110. * *
  111. * Resolution : *
  112. * *
  113. TAB2 = PASAPAS TAB1 ;
  114. * *
  115. *----------------------------------------------------------------------*
  116. * DEPOUILLEMENT *
  117. *----------------------------------------------------------------------*
  118. * *
  119. DEP1 = (TAB2 . 'DEPLACEMENTS' . 1) 'ENLE' 'LX' ;
  120. * *
  121. * Deformee : *
  122. DEFO0 = 'DEFO' (CUB1 'ET' SOL1) DEP1 0. 'VERT' ;
  123. DEFO1 = 'DEFO' (CUB1 'ET' SOL1) DEP1 1. 'ROUG' ;
  124. * *
  125. 'SI' IG1 ;
  126. 'TITR' 'Maillages non deforme (vert) et deforme (rouge).' ;
  127. 'TRAC' (DEFO0 'ET' DEFO1) ;
  128. 'FINS' ;
  129. * *
  130. * On verifie qu'on ne tranverse pas la surface : *
  131. ERR1 = ('MAXI' (SOL1 'COOR' 3))
  132. - ('MINI' ((CUB1 'COOR' 3) + ('EXCO' 'UZ' DEP1))) ;
  133. 'MESS' ' ****** ECART RELATIF POSITION DE CONTACT : ' ERR1 ;
  134. ERR0 = ERR1 ;
  135. * *
  136. * Definition des deplacements solutions et comparaison avec la *
  137. * solution EF : *
  138. 'SI' ('EGA' (TAB1 . 'HYPOTHESE_DEFORMATIONS') 'LINEAIRE') ;
  139. EPXX1 = ((-1. * UZ0) - 0.1) * 0.1 * 0.3 ;
  140. UXCUB1 = (('COOR' 1 CUB1) + 5.) * EPXX1 ;
  141. SOLUX1 = UXCUB1 'NOMC' 'UX' ;
  142. EPYY1 = EPXX1 ;
  143. UYCUB1 = ('COOR' 2 CUB1) * EPYY1 ;
  144. SOLUY1 = UYCUB1 'NOMC' 'UY' ;
  145. EPZZ1 = (UZ0 + 0.1) * 0.1 ;
  146. UZCUB1 = -0.1 + ((('COOR' 3 CUB1) - 0.1) * EPZZ1) ;
  147. SOLUZ1 = UZCUB1 'NOMC' 'UZ' ;
  148. SOLU1 = SOLUX1 'ET' SOLUY1 'ET' SOLUZ1 ;
  149. ERR2 = 'MAXI' ('ABS' ((SOLU1 - DEP1) / ('MAXI' 'ABS' SOLU1))) ;
  150. 'MESS' ' ****** ERREUR RELATIVE EN DEPLACEMENT : ' ERR2 ;
  151. ERR0 = ERR0 + ERR2 ;
  152. 'FINS' ;
  153. * *
  154. 'SI' IG1 ;
  155. 'TITR' 'Champ de deplacements.' ;
  156. 'TRAC' ('CHAN' 'CHAM' DEP1 MOD1 'NOEUDS') MOD1 ;
  157. 'FINS' ;
  158. * *
  159. * Comparaison des champs de contraintes : *
  160. * *
  161. 'SI' ('EGA' (TAB1 . 'HYPOTHESE_DEFORMATIONS') 'LINEAIRE') ;
  162. * Il ne faut pas faire le transport pour comparer a la solution castem
  163. * en lagrangien reactualise car il reviendrait a faire varier
  164. * le module d'young au cours du pas
  165. SIG1 = 'SIGM' MOD1 MAT1 SOLU1 ;
  166. **SIG1 = 'PICA' MOD1 SIG1 SOLU1 ;
  167. SIG2 = TAB2 . 'CONTRAINTES' . 1 ;
  168. ERR3 = ('MAXI' 'ABS' (SIG1 - SIG2)) / ('MAXI' 'ABS' SIG1) ;
  169. 'MESS' ' ****** ERREUR RELATIVE EN CONTRAINTE : ' ERR3 ;
  170. ERR0 = ERR0 + ERR3 ;
  171. 'FINS' ;
  172. * *
  173. * Test calcul critere de Tresca avec modele de contact dans objet *
  174. * modele en argument : *
  175. TRE1 = 'TRES' MOD0 SIG2 ;
  176. 'SI' IG1 ;
  177. 'TITR' 'Champ de contraintes.' ;
  178. 'TRAC' SIG2 MOD0 ;
  179. 'TRAC' TRE1 MOD0 'TITR' 'Critere de Tresca' ;
  180. 'FINS' ;
  181. * *
  182. * Visualisations des reactions : *
  183. * *
  184. 'SI' IG1 ;
  185. REAC1 = TAB2 . 'REACTIONS' . 1 ;
  186. VR1 = 'VECT' REAC1 'FORC' 'ROUG' ;
  187. 'TITR' 'Forces de reaction.' ;
  188. 'TRAC' VR1 ('ARET' MTOT1) ;
  189. 'FINS' ;
  190. * *
  191. ERRMAX1 = 1.E-4 ;
  192. * *
  193. 'OPTI' 'ECHO' 0 ;
  194. 'SAUT' 1 'LIGN' ;
  195. '---------------------- RESULTAT CAS-TEST ------------------------' ;
  196. 'MESS' ' ' ;
  197. 'Ecart relatif a la solution calculee sur les deplacements' ;
  198. 'et les contraintes :' ;
  199. 'MESS' ;
  200. 'SI' (ERR0 '<EG' ERRMAX1) ;
  201. 'MESS' ('CHAI' ' ==> MAXI ERREUR =' ERR0 ' < ' ERRMAX1 ) ;
  202. 'MESS' ' ' ;
  203. 'MESS' ' __________________' ;
  204. 'MESS' ' | |' ;
  205. 'MESS' ' | TEST REUSSI ! |' ;
  206. 'MESS' ' |________________|' ;
  207. 'SAUT' 1 'LIGN' ;
  208. 'SINO' ;
  209. 'MESS' ('CHAI' ' ==> MAXI ERREUR =' ERR0 ' > ' ERRMAX1 ) ;
  210. 'MESS' ' ' ;
  211. 'MESS' ' __________________' ;
  212. 'MESS' ' | |' ;
  213. 'MESS' ' | ERREUR ! |' ;
  214. 'MESS' ' |________________|' ;
  215. 'SAUT' 1 'LIGN' ;
  216. 'ERRE' 5 ;
  217. 'FINS' ;
  218. '------------------------ FIN CAS-TEST ---------------------------' ;
  219. * *
  220. 'FIN' ;
  221. 'OPTI' 'ECHO' 1 ;
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  

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