Télécharger gdef1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gdef1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. * *
  6. * TEST GDEF1 *
  7. * *
  8. * CISAILLEMENT PUR EN GRANDES DEFORMATIONS ELASTIQUES *
  9. * *
  10. * On compare avec les valeurs obtenues a la solution *
  11. * analytique *
  12. * *
  13. ********************************************************************
  14. *
  15. OPTI ECHO 0;
  16. TITRE ' TEST GDEF1 ' ;
  17. GRAPH = 'N';
  18. OPTI DIME 2 ELEM QUA4 MODE PLAN CONT;
  19. * ajout de option epsilon lineaire pour la precision des test!
  20. OPTION epsilon JAUMANN;
  21.  
  22. *
  23. *==================================================================*
  24. * geometrie : maillage *
  25. *==================================================================*
  26. *
  27. P1 = 0 0 ;
  28. P2 = 1. 0. ;
  29. P3 = 1. 1. ;
  30. P4 = 0. 1. ;
  31. L1 = P1 D 1 P2 ;
  32. L2 = P2 D 1 P3;
  33. L3 = P3 D 1 P4 ;
  34. L4 = P4 D 1 P1 ;
  35. SU = DALLER L1 L2 L3 L4 ;
  36. SI (NEG GRAPH 'N');
  37. TRAC SU ;
  38. FINSI;
  39. *
  40. *==================================================================*
  41. * modele - materiau - caracteristiques en Pa *
  42. * conditions aux limites *
  43. *==================================================================*
  44. *
  45. MO = MODE SU MECANIQUE ELASTIQUE ;
  46. YU = 1000. ; XNU = 0.3 ;
  47. DEUG = YU / ( 1. + XNU ) ;
  48. MA = MATE MO YOUN YU NU XNU ;
  49. BL1 = BLOQ DEPL P1 ;
  50. BL2 = BLOQ DEPL P2 ;
  51. BL3 = BLOQ UX L3 ;
  52. BL4 = BLOQ UY L3 ;
  53. FF = DEPI BL3 1. ;
  54. *
  55. *==================================================================*
  56. * Initialisation de la table pour appel a PASAPAS *
  57. *==================================================================*
  58. *
  59. TAB1 = TABLE;
  60. tab1.grands_deplacements=vrai;
  61.  
  62. TAB1 . BLOCAGES_MECANIQUES= BL1 ET BL2 ET BL3 ET BL4 ;
  63. TAB1 . MODELE = MO;
  64. TAB1 . CARACTERISTIQUES = MA ;
  65. EV1 = EVOL MANU X (PROG 0 5 10 15 20 )
  66. Y (PROG 0 5 10 5 0 );
  67. TAB1 . CHARGEMENT = CHAR MECA FF EV1 ;
  68. TAB1 . TEMPS_CALCULES = PROG 0.2 PAS 0.2 8;
  69. TMASAU=table;
  70. tab1 . 'MES_SAUVEGARDES'=TMASAU;
  71. TMASAU .'DEFTO'=VRAI;
  72. TMASAU .'DEFIN'=VRAI;
  73. *TMASAU .'DEFLO'=VRAI;
  74. PASAPAS TAB1 ;
  75.  
  76. *dep=tab1.deplacements; co = tab1.contraintes ;
  77. *list co.1;
  78. *list co.2;
  79. *fin;
  80.  
  81. *
  82. *==================================================================*
  83. * affichage des resultats *
  84. *==================================================================*
  85. *
  86. MESS ' RESULTATS TEST GDEF1 ';
  87. MESS ' -------------------- ';
  88. MESS ' On cherche les déformations';
  89. MESS ' Elles sont comparees a des valeurs de reference ';
  90. *
  91. *----------- r{sultats --------------------------------------------------
  92. *
  93. opti echo 1;
  94. AA = TAB1 . DEPLACEMENTS ;
  95. BB = TAB1 . CONTRAINTES ;
  96. CC = TAB1 . DEFORMATIONS ;
  97. N = DIME AA ;
  98. SP11= PROG ;
  99. SP12= PROG ;
  100. SP22= PROG ;
  101. EP11= PROG ;
  102. EP12= PROG ;
  103. EP22= PROG ;
  104. EPTH11= PROG ;
  105. EPTH12= PROG ;
  106. EPTH22= PROG ;
  107. EP = PROG ;
  108. EP = TAB1 . TEMPS_CALCULES ;
  109. REPE BOU (N - 1 ) ;
  110. CSX = CHAN CHPO MO ( BB . &BOU ) ;
  111. ESX = CHAN CHPO MO ( CC . 0 ) ;
  112. S11 = EXTR CSX SMXX P1 ;
  113. SP11 = SP11 ET ( PROG S11 ) ;
  114. S22 = EXTR CSX SMYY P1 ;
  115. SP22 = SP22 ET ( PROG S22 ) ;
  116. S12 = EXTR CSX SMXY P1 ;
  117. SP12 = SP12 ET ( PROG S12 ) ;
  118. E11 = EXTR ESX EPXX P1 ;
  119. EP11 = EP11 ET ( PROG E11 ) ;
  120. E22 = EXTR ESX EPYY P1 ;
  121. EP22 = EP22 ET ( PROG E22 ) ;
  122. E12 = EXTR ESX GAXY P1 ;
  123. EP12 = EP12 ET ( PROG (E12 * 0.5) ) ;
  124. DX = ( EXTR ( AA . &BOU ) UX P3) * 57.29577951 ;
  125. CDX = 0.5 * ( 1. - ( COS DX ) ) ;
  126. SDX = 0.5 * (SIN DX) ;
  127. EPTH11= EPTH11 ET ( PROG CDX ) ;
  128. EPTH12= EPTH12 ET ( PROG SDX ) ;
  129. EPTH22= EPTH22 ET ( PROG ( -1. * CDX ) ) ;
  130. *EX = LOG ( 1. + DX) ;
  131. *EP = EP ET ( PROG EX ) ;
  132. FIN BOU ;
  133. SP11 = SP11 / DEUG ;
  134. SP12 = SP12 / DEUG ;
  135. SP22 = SP22 / DEUG ;
  136. EV11 = EVOL ROSE MANU 'DEPLA' EP 'SIG11' SP11 ;
  137. EV12 = EVOL ROSE MANU 'DEPLA' EP 'SIG12' SP12 ;
  138. EV22 = EVOL ROSE MANU 'DEPLA' EP 'SIG22' SP22 ;
  139. FV11 = EVOL ROSE MANU 'DEPLA' EP 'EPS11' EP11 ;
  140. FV12 = EVOL ROSE MANU 'DEPLA' EP 'EPS12' EP12 ;
  141. FV22 = EVOL ROSE MANU 'DEPLA' EP 'EPS22' EP22 ;
  142. FVTH11 = EVOL TURQ MANU 'DEPLA' EP 'EPS11' EPTH11 ;
  143. FVTH12 = EVOL TURQ MANU 'DEPLA' EP 'EPS12' EPTH12 ;
  144. FVTH22 = EVOL TURQ MANU 'DEPLA' EP 'EPS22' EPTH22 ;
  145. SI (NEG GRAPH 'N');
  146. * DESS ( FVTH11 ET FV11) ;
  147. * DESS ( FVTH12 ET FV12) ;
  148. * DESS ( FVTH22 ET FV22) ;
  149. DESS ( EV11 ET FVTH11 ) ;
  150. DESS ( EV12 ET FVTH12 ) ;
  151. DESS ( EV22 ET FVTH22 ) ;
  152. FINSI;
  153. *
  154. * COMPARAISON SOL ANALYTIQUE
  155. * ecarts
  156. *
  157. ER11 =(MAXI ABS ( sp11 - EPTH11 ))/(MAXI ABS EPTH11);
  158. ER12 =(MAXI ABS ( sp12 - EPTH12 ))/(MAXI ABS EPTH12);
  159. ER22 =(MAXI ABS ( sp22 - EPTH22 ))/(MAXI ABS EPTH22);
  160. *
  161. MESS ' Ecart relatif en déformations eps11 :' ER11 ;
  162. MESS ' Ecart relatif en déformations eps12 :' ER12 ;
  163. MESS ' Ecart relatif en déformations eps22 :' ER22 ;
  164. *
  165. *
  166. *==================================================================*
  167. * code fonctionnement *
  168. *==================================================================*
  169. *
  170. * L'ecart maximum entre valeurs de reference et calculee doit etre
  171. * inferieur a 5.d-2
  172. *
  173. SAUT 2 LIGNE;
  174. SI (ER11 <EG 0.05 ) ;
  175. SI (ER12 <EG 0.05 ) ;
  176. SI (ER22 <EG 0.05 ) ;
  177. ERRE 0;
  178. SINON;
  179. ERRE 5;
  180. FINSI;
  181. SINON;
  182. ERRE 5;
  183. FINSI;
  184. SINON;
  185. ERRE 5;
  186. FINSI;
  187. *
  188. *
  189. FIN;
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  

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