Télécharger la plaquette Cast3M


fr en

Télécharger 1ddl.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :1ddl.dgibi 
  2. ************************************************************************
  3. ************************************************************************
  4. ***************************************
  5. *
  6. * Etude d'un système 1 ddl
  7. * Exemple d'utilisation de OSCI, SPO et SPON
  8. * 
  9. * D. Combescure aout 2006
  10. *
  11. ***************************************
  12. *GRAPH = VRAI;
  13. GRAPH = FAUX;
  14. *COMPLET = VRAI;
  15. COMPLET = FAUX;
  16. *
  17. opti dime 2 elem seg2;
  18. *
  19. fr1 = 5.;
  20. *
  21. dt = 0.005;
  22. prtime = prog 0. pas dt 100.;
  23. *
  24. * Durée de la rampe de montée et de descente du chargement
  25. dtramp = 0.001;
  26. *
  27. * Durée de l'impulsion
  28. *
  29. dtchar = 0.05;
  30. prtim1 = prog 0. 0.1 (0.1 + dtramp) (0.1 + dtramp + dtchar) 
  31.                      (0.1 + dtchar + (2.*dtramp)) 100.;
  32. prload1 = ((2.*pi*fr1)**2)*(prog 0 0. 1. 1. 0. 0.);
  33. evload1 = evol manu prtime (ipol prtime prtim1 prload1);
  34. prtimec = prog 0. pas dt 1.;
  35. *
  36. SI (EGA GRAPH VRAI);
  37.  dess evload1 xbord 0. 1. ;
  38. FINSI;
  39. *
  40. * Réponse 1ddl elastique
  41. *
  42. ev1 = osci evload1 amor 0.05 freq fr1 temps prtimec depl 0.;
  43. SI (EGA GRAPH VRAI);
  44.  dess ev1 xbord 0. 1.;
  45.  @excel1 ev1 'ev1.excel'; 
  46. FINSI;
  47. *
  48. * 
  49. *
  50. * Calcul spectre oscillateur
  51. *
  52. SI COMPLET; 
  53.  lfreq = (prog 0.01 pas 0.1 20.);
  54. SINON;
  55.  lfreq  = prog 1. pas 1. 20.;
  56. FINSI;
  57. spev1 = spo evload1 amor (prog 0.05)
  58.           freq lfreq
  59.           acce coul rouge;
  60. *
  61. * Calcul spectre inélastique - oscilateur élastoplastique parfait
  62. *
  63. * Option DEMA => on trace la pseudoaccélération correspondant 
  64. *                à la limite élastique (pour le dimensionnement)
  65. *
  66. * Ductilité 1
  67. spnlev1 = spon DEMA SIGN evload1 SPEL spev1 ACCE
  68.           AMOR (prog 0.05) CINE (prog 1. 0.00) ACCE; 
  69. * Ductilité 2
  70. spnlev2 = spon DEMA SIGN evload1 SPEL spev1 ACCE
  71.           AMOR (prog 0.05) CINE (prog 2. 0.00) ACCE;   
  72. * Ductilité 5
  73. spnlev5 = spon DEMA SIGN evload1 SPEL spev1 ACCE
  74.           AMOR (prog 0.05) CINE (prog 5. 0.00) ACCE;   
  75. *
  76. SI (EGA GRAPH VRAI);
  77.  dess (spnlev1 et spnlev2 et spnlev5 et spev1);
  78. FINSI;    
  79. *
  80. ERR = maxi ((extr spnlev1 ordo) - (extr spev1 ordo)) abs;
  81. *
  82. SI (ERR > 1.D-3);
  83.  ERRE 4;
  84. FINSI;
  85. *
  86. FIN;
  87. *
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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