Télécharger redumode.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : redumode.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * démonstration pour reduire le nombre de modèles
  6. *
  7. debp affpro mai*maillage ia*entier;
  8. mo1= model mai mecanique elastique plastique parfait cons co1;
  9. mo2= model mai mecanique elastique plastique parfait cons co2;
  10. mo3= model mai mecanique elastique plastique parfait cons co3;
  11. mo4= model mai mecanique elastique plastique parfait cons co4;
  12. mo5= model mai mecanique elastique plastique parfait cons co5;
  13. mo6= model mai mecanique elastique plastique parfait cons co6;
  14. mo7= model mai mecanique elastique plastique parfait cons co7;
  15. mo8= model mai mecanique elastique plastique parfait cons co8;
  16. mo9= model mai mecanique elastique plastique parfait cons co9;
  17. mo10= model mai mecanique elastique plastique parfait cons co10;
  18. mot= mo1 et mo2 et mo3 et mo4 et mo5 et mo6 et mo7 et mo8 et
  19. mo9 et mo10;
  20. ib= 20 + ia;
  21. ma1= mo1 mate young 20000 nu 0.3 sigy ib;
  22. ma2= mo2 mate young 20000 nu 0.3 sigy ib;
  23. ma3= mo3 mate young 20000 nu 0.3 sigy ib;
  24. ma4= mo4 mate young 20000 nu 0.3 sigy ib;
  25. ma5= mo5 mate young 20000 nu 0.3 sigy ib;
  26. ma6= mo6 mate young 20000 nu 0.3 sigy ib;
  27. ma7= mo7 mate young 20000 nu 0.3 sigy ib;
  28. ma8= mo8 mate young 20000 nu 0.3 sigy ib;
  29. ma9= mo9 mate young 20000 nu 0.3 sigy ib;
  30. ma10= mo10 mate young 20000 nu 0.3 sigy ib;
  31. mat= ma1 et ma2 et ma3 et ma4 et ma5 et ma6 et ma7 et ma8 et
  32. ma9 et ma10;
  33. finproc mot mat;
  34.  
  35.  
  36. * ajout de option epsilon lineaire pour la precision des test!
  37. OPTION epsilon lineaire;
  38.  
  39. opti dime 2 elem qua4 mode plan defo;
  40. pa= 0 0 ; pb = 10 0;
  41. liab=pa droi 10 pb;
  42. su = liab tran 5 (0 10);
  43. motot= su model mecanique elastique plastique parfait;
  44. matot= mate motot young 200000 nu 0.3 sigy 200.;
  45. cl1= bloq uy liab;
  46. cl2= bloq ux pa;
  47. cltot= cl1 et cl2;
  48.  
  49. ff = force ( cote 3 su) (0 100);
  50.  
  51. ri = rigi motot matot;
  52. de = reso ( ri et cltot) ff;
  53.  
  54. sig= sigma de motot matot;
  55. *trac sig motot matot;
  56.  
  57.  
  58.  
  59. nbe= nbel su;
  60. repe bou nbe;
  61.  
  62. maiel= elem su &bou;
  63. moel mael= affpro maiel &bou;
  64. si ( &bou ega 1) ;
  65. mott = moel; matt= mael;
  66. sinon;
  67. mott = mott et moel;
  68. matt = matt et mael;
  69. finsi;
  70. fin bou;
  71. ritt = rigi mott matt;
  72. dett = reso ( ritt et cltot) ff;
  73.  
  74.  
  75. diffde= de - dett; list ( maxi diffde);
  76.  
  77. temps;
  78. lt= prog 0. 1.;
  79. cf= prog 0. 1.;
  80. ev = evol manu 'temps' lt 'coef' cf;
  81. ch = char 'MECA' ff ev;
  82. ta = table ;
  83. ta.modele= motot;
  84. ta.caracteristiques = matot;
  85. ta . blocages_mecaniques=cltot;
  86. ta.chargement = ch;
  87. ta.temps_calcules = lt;
  88. de1 = ta.deplacements . 1;
  89. temps;
  90. difde1= de - de1;
  91. list ( maxi difde1);
  92. tat = table ;
  93. tat.modele= mott;
  94. tat.caracteristiques = matt;
  95. tat . blocages_mecaniques=cltot;
  96. tat.chargement = ch;
  97. tat.temps_calcules = lt;
  98. pasapas tat;
  99.  
  100. temps;
  101. de1t= tat.deplacements . 1;
  102.  
  103. dide1t= de - de1t; list ( maxi dide1t) ;
  104.  
  105.  
  106.  
  107. *
  108. * on reduit le nombre de modele en s'appuyant sur les constituant
  109. *
  110. debp reduc coi*mot motto*mmodel matto*mchaml;
  111. * petite procedure qui prend tous les models ayant le constituant coi
  112. * et qui les rassemble en un seul. La reduction du champ par element
  113. * doit etre faite simultaneemnt
  114.  
  115. moco1= extr motto cons coi;
  116. maico1= extr moco1 maillage;
  117. motco1= modele maico1 mecanique elastique plastique parfait cons coi;
  118. matco1= redu matto motco1;
  119. finproc motco1 matco1;
  120.  
  121. temps;
  122. mot1 mat1 = reduc co1 mott matt;
  123. mot2 mat2 = reduc co2 mott matt;
  124. mot3 mat3 = reduc co3 mott matt;
  125. mot4 mat4 = reduc co4 mott matt;
  126. mot5 mat5 = reduc co5 mott matt;
  127. mot6 mat6 = reduc co6 mott matt;
  128. mot7 mat7 = reduc co7 mott matt;
  129. mot8 mat8 = reduc co8 mott matt;
  130. mot9 mat9 = reduc co9 mott matt;
  131. mot10 mat10 = reduc co10 mott matt;
  132.  
  133.  
  134. mov= mot1 et mot2 et mot3 et mot4 et mot5 et mot6 et mot7 et
  135. mot8 et mot9 et mot10;
  136.  
  137. mav= mat1 et mat2 et mat3 et mat4 et mat5 et mat6 et mat7 et
  138. mat8 et mat9 et mat10;
  139.  
  140. tav= table;
  141. tav.modele= mov;
  142. tav.caracteristiques = mav;
  143. tav . blocages_mecaniques=cltot;
  144. tav.chargement = ch;
  145. tav.temps_calcules = lt;
  146. pasapas tav;
  147. temps;
  148.  
  149. de1tv= tav.deplacements . 1;
  150. dide1tv= de - de1tv ; list ( maxi dide1tv);
  151.  
  152. aa = prog (maxi dide1tv) (maxi dide1t) (maxi diffde)
  153. ( maxi difde1);
  154. ab = maxi aa;
  155. mess ' erreur max ' ab ' comparee a 1.e-11';
  156. si (( maxi aa) > 1.e-11 ) ;
  157. erreur 5;
  158. finsi;
  159. fin;
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  

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