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

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