Télécharger la plaquette Cast3M


fr en

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

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