Télécharger la plaquette Cast3M


fr en

Télécharger tufi_relax.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  tufi_relax.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTION ECHO 0;
  5. opti dime 3 mode TRID elem seg2 ;
  6. rap=10.;T1=10.;RMOY=T1 * rap;REXT1=RMOY + (T1 / 2.);
  7. P1=0. 0. 0.;P2=0. 0. 0.;                                                        
  8. TUFSS=MANU SEG2 P1 P2;                                                          
  9. OBJTUFI= MODE TUFSS  MECANIQUE ELASTIQUE FLUAGE NORTON TUFI;
  10. CARTOT= mate OBJTUFI YOUN 2.E5 NU 0.3
  11. AF1 3.e-15 AF2 5. AF3 1. SMAX 200. RAYO REXT1 EPAI T1
  12. VX 1. VY 0. VZ 0. VXF 0. VYF 0. VZF 1. ANGL 120;
  13. CDL=BLOQ depla rota p1;
  14. RIG1= RIGI OBJTUFI CARTOT;RIG = RIG1 ET CDL;
  15. DD = RESO RIG (MOME 100.E5 'MY' P2);
  16. cdli=bloq 'RY' P2;CDL=CDL et cdli;
  17. RIG = RIG1 ET CDL;
  18. vdepli=extr DD 'RY' P2;
  19. depli=depi cdli vdepli;
  20. DD = RESO RIG depli;
  21. M0=react DD cdli;vM0=extr M0 'MY' P2;
  22. CMM=vdepli / vM0;
  23. SS = SIGMA CARTOT OBJTUFI DD;
  24. EVT = EVOL MANU 'T' (PROG 0. 1.E10) 'F(T)' (PROG 1. 1.);
  25. dt=10.;LIS = PROG 0. pas dt 110. ;
  26. FDT=CHAR 'DIMP' depli EVT;                                                                
  27. TENTR=TABLE;
  28. TENTR.CARACTERISTIQUES=CARTOT;
  29. TENTR.MODELE=OBJTUFI;
  30. TENTR.CHARGEMENT=FDT;
  31. TENTR.'CONTRAINTES'=tabl;
  32. TENTR.'DEPLACEMENTS'=tabl;
  33. TENTR.'CONTRAINTES' . 0 =SS;
  34. TENTR.'DEPLACEMENTS' . 0 = DD;
  35. TENTR.BLOCAGES_MECANIQUES=CDL;
  36. TENTR.TEMPS_CALCULES= LIS;
  37. pasapas TENTR;
  38. t10=TENTR.TEMPS . 10;
  39. U10=TENTR.DEPLACEMENTS . 10;
  40. VAR10=TENTR.VARIABLES_INTERNES . 10;
  41. M10=react U10 cdli;vM10=extr M10 'MY' P2;phi_el=CMM * vM10;
  42. phi_f=vdepli - phi_el;
  43. vC10=extr VAR10 'JPOI' 1 1 1;
  44. phi_fb=extr VAR10 'EPSE' 1 1 1;
  45. U9=TENTR.DEPLACEMENTS . 9;
  46. Mi=react U9 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
  47. phi_f9=vdepli - phi_el;
  48. U10=TENTR.DEPLACEMENTS . 10;
  49. Mi=react U10 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
  50. phi_f10=vdepli - phi_el;
  51. U11=TENTR.DEPLACEMENTS . 11;
  52. Mi=react U11 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
  53. phi_f11=vdepli - phi_el;
  54. phi_fpt1=(phi_f10 - phi_f9) / dt;
  55. phi_fpt2=(phi_f11 - phi_f10) / dt;
  56. phi_fpt=(phi_fpt1 + phi_fpt2) / 2.;
  57.  
  58. * COMPARAISON A LA SOLUTION ANALYTIQUE
  59. Mth10=7.747E6;phi_f10=1.850E-4;vCth10=4.576E-3;
  60. ph_fpt10=1.020E-6;
  61.  
  62. errM=abs ((vM10 - Mth10) / Mth10);
  63. errphif=abs ((phi_f - phi_f10) / phi_f10);
  64. errphifb=abs ((phi_fb - phi_f10) / phi_f10);
  65. errC=abs ((vC10 - vCth10) / vCth10);
  66. errphfp=abs ((phi_fpt - ph_fpt10) / ph_fpt10);                    
  67. err=errM * errphif * errphifb * errC * errphfp;
  68. SI (ERR < (1.E-3 ** 5))ERRE 0;SINO;ERRE 5;FINSI;
  69. fin;
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  

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