Télécharger frocable.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : frocable.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * fichier: frocable.dgibi
  6. *
  7.  
  8. opti dime 2 mode plan cont elem qua4;
  9.  
  10. *=== GEOMETRIES===
  11.  
  12. P1 = 0 0;
  13. P2 = 5 0;
  14. liab = P1 droi 5 P2;
  15. su = liab trans 15 (0 3);
  16.  
  17. P3 = 0 1.5;
  18. P4 = 5 1.5;
  19. *P5 = 0 2.5;
  20. *P6 = 5 2.5;
  21.  
  22. cab1 = p3 droi 11 P4;
  23. *cab2 = p5 droi 10 p6;
  24. *cabt = cab1 et cab2;
  25.  
  26. bordy = cote 4 su;
  27.  
  28. cab1= coul rouge cab1;
  29. su = coul vert su;
  30. toto= cab1 et su;
  31.  
  32. *trac toto;
  33.  
  34. *opti donn 5;
  35. *=== MODELES===
  36.  
  37. MODCAB1 = model cab1 mecanique elastique barre;
  38.  
  39. *list modcab1;
  40. *list modcab2;
  41. cabfr= impf cab1;
  42. *MODFRO1 = model cabfr contact frottant frocable modcab1 su;
  43. modfro1= model cab1 contact frottant frocable su;
  44. *list modfro1;
  45.  
  46. * juste pour voir azppel a rfco
  47. *cdep crr = rfco modfro1 vrai ;
  48. *list bb;
  49. *opti donn 5;
  50.  
  51. *list ( frig modfro1);
  52. MODBET = model su mecanique elastique;
  53.  
  54. *list modfro1;
  55.  
  56. *=== MATERIAUX===
  57.  
  58. Ebet=34567e6;
  59. nubet=0.2;
  60. rhobet=2240;
  61.  
  62. matbet = mate MODBET YOUN Ebet NU nubet RHO rhobet;
  63.  
  64. Ecab=200000e6;
  65. nucab=0.3;
  66. rhocab=8000;
  67.  
  68. frotcab1=0.18;
  69.  
  70. phicab1 =0.002;
  71. * on multiplie phicab1 par 10 pour avoir une influence
  72. * plus importante
  73. phicab1=phicab1*10;
  74.  
  75. matcab1 = mate MODCAB1 YOUN Ecab NU nucab RHO rhocab sect 0.15 ;
  76. matfro1 = mate MODFRO1 FF frotcab1 PHIF phicab1;
  77.  
  78. mattot= matcab1 et matfro1 et matbet;
  79.  
  80. modtot= modfro1 et modcab1 et modbet;
  81.  
  82.  
  83. *list matfro1;
  84.  
  85.  
  86. *===RIGIDITES===
  87.  
  88. ribet= RIGI MODBET matbet;
  89. ricab= RIGI MODCAB1 matcab1;
  90.  
  91. *list riac1;
  92. *list riac2;
  93.  
  94. *===PRECONTRAINTES===
  95.  
  96. TAB = table;
  97. TAB.'FF '=frotcab1;
  98. TAB.'PHIF'=phicab1;
  99. TAB.'GANC'=0;
  100. TAB.'RMU0'=0.43;
  101. TAB.'FPRG'=1.7d9;
  102. TAB.'RHL0'=2.5;
  103.  
  104. *===CALCUL DES FORCES===
  105.  
  106. PRE1 = PREC MODCAB1 MATCAB1 TAB 3.e5;
  107. PRE2 = PREC MODCAB1 MATCAB1 TAB 6.e5;
  108. PRE3 = PREC MODCAB1 MATCAB1 TAB 9.e5;
  109. PRE4 = PREC MODCAB1 MATCAB1 TAB 1.2e6;
  110. PRE5 = PREC MODCAB1 MATCAB1 TAB 1.5e6;
  111. PRE6 = PREC MODCAB1 MATCAB1 TAB 1.8e6;
  112. PRE7 = PREC MODCAB1 MATCAB1 TAB 2.1e6;
  113. PRE8 = PREC MODCAB1 MATCAB1 TAB 2.4e6;
  114. PRE9 = PREC MODCAB1 MATCAB1 TAB 2.7e6;
  115. PRE10 = PREC MODCAB1 MATCAB1 TAB 3.e6;
  116. *list pre10;
  117. el1=cab1 elem 1;
  118. el2=cab1 elem 4;
  119.  
  120. FORC1 = FORC FX 3.e6 p4;
  121. *list forc1;
  122.  
  123. FORC2 = FORC FX -3.e6 p3;
  124.  
  125. depenc= bloque depl ( bordy );
  126.  
  127. deppoi= bloque maxi UX P3;
  128. *deppoi= bloque mini UX P4;
  129. *riac1 = rela glissant modfro1 su 0.001 ;
  130. *list riac1;
  131. *opti donn 5;
  132.  
  133. *cl= riac1 et depenc et deppoi;
  134. cl= depenc et deppoi;
  135.  
  136. *rig= rigi modtot mattot;
  137. *rigtot=rig et cl;
  138. *dep=reso rigtot forc1;
  139.  
  140. *VEC1 = VECT forc1 1.E-6 FX FY turq;
  141. *trac (vec1) toto;
  142.  
  143. *def1=defo toto dep 0. blanc;
  144. *def2=defo toto dep rouge;
  145.  
  146. *trac (def1 et def2);
  147.  
  148.  
  149. *sig1=sigma modcab1 matcab1 dep;
  150. *sxx1=extr sig1 comp;
  151. *xx =exco 'EFFX' sig1;list xx;
  152.  
  153. *===PAS A PAS===
  154.  
  155. *list1=prog 0 pas 0.5 1. pas 1. 3.;
  156. *list2=prog 0 1 1 1 1;
  157. *list3=prog 0 0 1 1.01 1.02;
  158.  
  159. list1=prog 0 pas 0.5 1 2 3 ;
  160. list2=prog 0 1 1 1 1 ;
  161. list3=prog 0 0 0.99 pas 0.005 1. ;
  162.  
  163. evt1=evol manu list1 list2;
  164. evt2=evol manu list1 list3;
  165.  
  166. *dess (evt1 et evt2);
  167.  
  168.  
  169. cha1=char meca forc1 evt1;
  170. cha2=char meca forc2 evt2;
  171.  
  172. chat=cha1 et cha2;
  173.  
  174. tab1 = table;
  175. tab1.'BLOCAGES_MECANIQUES'=cl;
  176. tab1.'MODELE'=modtot;
  177. tab1.'CHARGEMENT'=chat;
  178. list mattot;
  179. tab1.'CARACTERISTIQUES'=mattot;
  180. tab1.'TEMPS_CALCULES'=list1;
  181. tab1.'PRECISION'=1.e-7;
  182.  
  183. pasapas tab1;
  184.  
  185. *opti sauv 'donnees.sauv';
  186. *sauv tab1;
  187. debproc ecri ta*table mo*mmodel ge*maillage;
  188. co= ta.contraintes ;
  189. de= ta.deplacements ;
  190. re = ta. reactions ;
  191. cab= extr mo maillage;
  192. na = (dime co ) - 1;
  193. repe bo na;
  194. mess ' contraintes au pas ' &bo;
  195. list (( co . &bo redu mo) exco EFFX);
  196. fin bo;
  197. *repe ba na;
  198. *mess ' deplacements au pas ' &ba;
  199. *list ( de . &ba redu ge);
  200. **fin ba;
  201. repe bi na;
  202. mess ' reactions au pas ' &bi;
  203. list ( re . &bi redu ge) ;
  204. fin bi;
  205. finp ;
  206. *ecri tab1 modcab1 cab1;
  207. ffr1=(tab1.contraintes . 1) exco EFFX ;
  208. list ffr1;
  209. valref= ffr1 extr EFFX 1 6 2;
  210. ffr4= (tab1.contraintes . 4) exco EFFX ;
  211. minffr4= mini ffr4;
  212. * teste de la valeur min trouvée au temps 1 par rapport a
  213. * solution fournie par precontrainte( pre10)
  214. xpre= mini pre10;
  215. xsol= mini ffr1;
  216. mess ' xsol ' xsol ' xpre ' xpre;
  217. ersol= abs ( (xsol - xpre) / xpre);
  218. mess ' erreur sur la precontrainte ' ersol ;
  219.  
  220. si ( ersol > 1.e-2) ;
  221. erreur 5 ;
  222. finsi;
  223.  
  224. * on compare la valeur trouvé au temps 1
  225. *pour le point central (valref) et celle trouvée
  226. * pour le temps 4
  227.  
  228. er=abs((minffr4 - valref)/valref);
  229. mess ' erreur en % ' er;
  230. si ( er > 0.01);
  231. erreur 5;
  232. finsi;
  233. *list pre10;
  234. *list ffr1;
  235. fin;
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  

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