Télécharger stru2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : stru2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Stru2.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. * Pour visualiser les traces, mettre GRAPH a 'O' :
  8. GRAPH = 'MOT' 'N' ;
  9. OPTION ECHO 0;
  10. *
  11. MESS ' **********************************************';
  12. MESS ' * *';
  13. MESS ' * COMMISSION V.P.C.S. TEST STRU2.DGIBI *';
  14. MESS ' * TEST SSLS 33/91 *';
  15. MESS ' * GROUPE : STRUCTURE ASSEMBLEE *';
  16. MESS ' * STATIQUE LINEAIRE *';
  17. MESS ' * *';
  18. MESS ' * TEST : PLAQUE CARREE ORTHOTROPE *';
  19. MESS ' * ENCASTREE *';
  20. MESS ' * *';
  21. MESS ' * REMARQUE : CAS TEST FAIT AVEC ELEMENTS *';
  22. MESS ' * QUA4 EN CONTRAINTES PLANES *';
  23. MESS ' * *';
  24. MESS ' * LES VALEURS DE REFERENCES SONT LES *';
  25. MESS ' * MOYENNES DE RESULTATS OBTENUS PAR SIX *';
  26. MESS ' * CODES DIFFERENTS *';
  27. MESS ' * *';
  28. MESS ' * TEST MODIFIE LE 3/07/92 *';
  29. MESS ' **********************************************';
  30. *
  31. OPTION DIME 2 ELEM QUA4;
  32. *
  33. *-----------------------------------------------------*
  34. * GEOMETRIE *
  35. *-----------------------------------------------------*
  36. *
  37. PA = 0.5 0.5;
  38. PB = -0.5 0.5;
  39. PC = -0.5 -0.5;
  40. PD = 0.5 -0.5;
  41. PE = 0. 0.;
  42. *
  43. D1 = PA DROIT 10 PB;
  44. D2 = PB DROIT 10 PC;
  45. D3 = PC DROIT 10 PD;
  46. D4 = PD DROIT 10 PA;
  47. *
  48. SUR1 = DALLAGE D1 D2 D3 D4 PLAN;
  49. SUR = SUR1 ET PE;
  50. *
  51. ELIM SUR 0.001;
  52. *
  53. SI (NEG GRAPH 'N');
  54. TITRE 'PLAQUE CARREE ORTHOTROPE';
  55. TRACE SUR ;
  56. FINSI ;
  57. *
  58. NBNOSUR = NBNO SUR1;
  59. NBELSUR = NBEL SUR1;
  60. *
  61. *-----------------------------------------------------*
  62. * MODELE *
  63. *-----------------------------------------------------*
  64. *
  65. MOD = MODE SUR1 MECANIQUE ELASTIQUE ORTHOTROPE QUA4;
  66. *
  67. *-----------------------------------------------------*
  68. * CARACTERISTIQUES - MATERIAUX *
  69. *-----------------------------------------------------*
  70. *
  71. DIR1 = (COS 30.) (SIN 30.);
  72. *
  73. MAT = MATE MOD DIRE DIR1 YG1 2.5737E10 YG2 7.377E9
  74. YG3 7.377E9 NU12 0.31 NU23 0.31 NU13 0.31 G12 2.319E9;
  75. *
  76. *-----------------------------------------------------*
  77. * CONDITIONS AUX LIMITES - RIGIDITES *
  78. *-----------------------------------------------------*
  79. *
  80. CL = BLOQUE DEPL ROTA D2;
  81. RI = RIGI MOD MAT ;
  82. RIG = RI ET CL;
  83. *
  84. *-----------------------------------------------------*
  85. * CHARGEMENTS *
  86. *-----------------------------------------------------*
  87. * IL S'AGIT D'UN CHARGEMENT LINEIQUE DE 1e4 NEWTONS/M
  88. * SUR AD. L'EPAISSEUR EST DE 0.01M. ON A DONC UNE
  89. * FORCE PAR UNITE D'EPAISSEUR DE 1E6 N/M
  90.  
  91. FTOT = 1.E6;
  92. NBNOD4 = NBNO D4;
  93. F5 = FTOT / ((2 * NBNOD4) - 2);
  94. F50 = F5 0.;
  95. I = 1;
  96. NBI = NBNOD4 - 1;
  97. CHA = FORC (0. 0.) (D4 POIN 1);
  98. REPETER BLOC1 NBI;
  99. CHA1 = FORC F50 (D4 POIN I);
  100. CHA2 = FORC F50 (D4 POIN (I+1));
  101. CHA = CHA ET CHA1 ET CHA2;
  102. I = I + 1;
  103. FIN BLOC1;
  104.  
  105. *-----------------------------------------------------*
  106. *RESULTATS : DEPLACEMENTS, CONTRAINTES *
  107. *-----------------------------------------------------*
  108. *
  109. DE = RESO RIG CHA;
  110. *
  111. SI (NEG GRAPH 'N');
  112. TITR ' PLAQUE CARREE ORTHOTROPE : DEFORMEE .';
  113. DEF0 = DEFO SUR DE 0.;
  114. DEF1 = DEFO SUR DE 10. ROUGE;
  115. TRAC (DEF0 ET DEF1);
  116. FINSI;
  117. *
  118. UA = EXTR DE UX PA; UA1 = 1.098e-4;
  119. ERGUA = 100 * (ABS ((UA1 - UA) / UA1));
  120. *
  121. UD = EXTR DE UX PD; UD1 = 1.003e-4;
  122. ERGUD = 100 * (ABS ((UD1 - UD) / UD1));
  123. *
  124. UE = EXTR DE UX PE; UE1 = 4.768E-5;
  125. ERGUE = 100 * (ABS ((UE1 - UE) / UE1));
  126. *
  127. VA = EXTR DE UY PA; VA1 = -1.309E-4;
  128. ERGVA = 100 * (ABS ((VA1 - VA) / VA1));
  129. *
  130. VD = EXTR DE UY PD; VD1 = -7.134E-5;
  131. ERGVD = 100 * (ABS ((VD1 - VD) / VD1));
  132. *
  133. VE = EXTR DE UY PE; VE1 = -4.965E-5;
  134. ERGVE = 100 * (ABS ((VE1 - VE) / VE1));
  135. *
  136. MESS ' Resultats : test stru2.dgibi';
  137. MESS ' ----------------------------';
  138. MESS ' Deplacements';
  139. MESS ' ------------';
  140. *
  141. MESS ' Point A :';
  142. MESS ' ---------';
  143. MESS ' Deplacement theorique horizontal :' UA1 'm';
  144. MESS ' Deplacement calcule horizontal :' UA 'm';
  145. MESS ' Soit un ecart de :' ERGUA '%';
  146. MESS ' Deplacement theorique vertical :' VA1 'm';
  147. MESS ' Deplacement calcule vertical :' VA 'm';
  148. MESS ' Soit un ecart de :' ERGVA '%';
  149. *
  150. MESS ' Point D :';
  151. MESS ' ---------';
  152. MESS ' Deplacement theorique horizontal :' UD1 'm';
  153. MESS ' Deplacement calcule horizontal :' UD 'm';
  154. MESS ' Soit un ecart de :' ERGUD '%';
  155. MESS ' Deplacement theorique vertical :' VD1 'm';
  156. MESS ' Deplacement calcule vertical :' VD 'm';
  157. MESS ' Soit un ecart de :' ERGVD '%';
  158. *
  159. MESS ' Point E :';
  160. MESS ' ---------';
  161. MESS ' Deplacement theorique horizontal :' UE1 'm';
  162. MESS ' Deplacement calcule horizontal :' UE 'm';
  163. MESS ' Soit un ecart de :' ERGUE '%';
  164. MESS ' Deplacement theorique vertical :' VE1 'm';
  165. MESS ' Deplacement calcule vertical :' VE 'm';
  166. MESS ' Soit un ecart de :' ERGVE '%';
  167. *
  168. *SAUT 1 LIGN;
  169. *MESS ' Contraintes au point E';
  170. *MESS ' ----------------------';
  171. *SAUT 1 LIGN;
  172. *
  173. CHAM1 = SIGM MOD MAT DE;
  174. CHAM2 = CHAN NOEUD MOD CHAM1;
  175. CHPO1 = CHAN CHPO MOD CHAM2;
  176. *
  177. SMXXE = EXTR CHPO1 SMXX PE;
  178. SMXXE = SMXXE / 1.E6;
  179. SMXXE1 = 1.004;
  180. ERGXX = 100 * (ABS ((SMXXE1 - SMXXE) / SMXXE1));
  181. *
  182. SMYYE = EXTR CHPO1 SMYY PE;
  183. SMYYE = SMYYE / 1.E6;
  184. SMXYE = EXTR CHPO1 SMXY PE;
  185. SMXYE = SMXYE / 1.E6;
  186. *
  187. *MESS ' Contrainte theorique horizontale :' SMXXE1
  188. *'MPa';
  189. *MESS ' Contrainte calculee horizontale :' SMXXE
  190. *'MPa';
  191. *MESS ' Soit un ecart de :' ERGXX '%';
  192. *MESS ' Contrainte calculee verticale :' SMYYE
  193. *'MPa';
  194. *MESS ' Contrainte calculee de cisaillement :' SMXYE
  195. *'MPa';
  196. *
  197. MESS ' Nombre de noeuds :' NBNOSUR;
  198. MESS ' Nombre d elements :' NBELSUR;
  199. *
  200. * UTILISATION DE LA PROCEDURE AFFICHE
  201. *
  202. SI (NEG GRAPH 'N');
  203. AFFICHE RIG CHA SUR 1.E-6 1.E-6;
  204. FINSI;
  205. *
  206. SI (NEG GRAPH 'N');
  207. CHAM10 = EXCO SMXX CHAM1;
  208. CHAM10 = 'CHAN' 'TYPE' CHAM10 'SCALAIRE' ;
  209. TITRE 'CONTRAINTE SMXX';
  210. TRACE MOD CHAM10 SUR;
  211. *
  212. CHAM11 = EXCO SMYY CHAM1;
  213. CHAM11 = 'CHAN' 'TYPE' CHAM11 'SCALAIRE' ;
  214. TITRE 'CONTRAINTE SMYY';
  215. TRACE MOD CHAM11 SUR;
  216. *
  217. CHAM12 = EXCO SMXY CHAM1;
  218. CHAM12 = 'CHAN' 'TYPE' CHAM12 'SCALAIRE' ;
  219. TITRE 'CONTRAINTE SMXY';
  220. TRACE MOD CHAM12 SUR;
  221. *
  222. CHAM13 = VMIS MOD CHAM1;
  223. CHAM13 = 'CHAN' 'TYPE' CHAM13 'SCALAIRE' ;
  224. TITRE 'CONTRAINTE DE VON-MISES';
  225. TRACE MOD CHAM13 SUR;
  226. FINSI ;
  227. *
  228. *=====================================================*
  229. * code fonctionnement *
  230. *=====================================================*
  231. *
  232. ERGMAX=MAXI (PROG ERGUA ERGUD ERGUE ERGVA ERGVD ERGVD);
  233. *
  234. * Il y a erreur si un des pourcentages ERG*
  235. * est superieur a 5%
  236. *
  237. SI (ERGMAX <EG 5.);
  238. ERRE 0 ;
  239. SINO;
  240. ERRE 5 ;
  241. FINSI;
  242. *
  243. TEMPS;
  244. *
  245. FIN;
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  

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