Télécharger la plaquette Cast3M


fr en

Télécharger A1DDL.dgibi

Retour à la liste

Numérotation des lignes : 

  1. **********************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * EXEMPLE A1DDL.dgibi
  5. **********************************************
  6. * Entrée :      Chargement sismique
  7.  
  8. **********************************************
  9. * Sortie :      Sans objet
  10. *               
  11. **********************************************
  12. * Commentaire : Test de la procedure  
  13. *               @A1DDL.PROCEDUR
  14. **********************************************
  15. * Developpeur : Benjamin Richard
  16. *               CEA, DEN, DANS, DM2S, SEMT, EMSI
  17. *               benjamin.richard@cea.fr
  18. **********************************************
  19. **********************************************
  20. * OPTIONS DE CALCUL
  21. **********************************************
  22. * repertoire des fichiers "divers"
  23. DIVERS = VENV 'CASTEM_DIVERS';
  24. *
  25.  
  26. RUN = 1;
  27.  
  28. **********************************************
  29. * ACQUISITION RUN I
  30. **********************************************
  31. * RUN 1
  32. SI (RUN EGA 1);
  33. *OPTI ACQU './exp/21time.txt';
  34. OPTI ACQU ('CHAINE' DIVERS '/21time.txt');
  35. *ACQU TIME*LISTREEL 8192;
  36. ACQU TIME*LISTREEL 550;
  37. *OPTI ACQU './exp/21axtab.txt';
  38. OPTI ACQU ('CHAINE' DIVERS '/21axtab.txt');
  39. *ACQU AXTAB*LISTREEL 8192;
  40. ACQU AXTAB*LISTREEL 550;
  41. FINSI;
  42.  
  43. * RUN 2
  44. SI (RUN EGA 2);
  45. OPTI ACQU '.\exp\22time.txt';
  46. ACQU TIME*LISTREEL 8192;
  47. OPTI ACQU '.\exp\22axtab.txt';
  48. ACQU AXTAB*LISTREEL 8192;
  49. FINSI;
  50.  
  51. * RUN 3
  52. SI (RUN EGA 3);
  53. OPTI ACQU '.\exp\23time.txt';
  54. ACQU TIME*LISTREEL 8192;
  55. OPTI ACQU '.\exp\23axtab.txt';
  56. ACQU AXTAB*LISTREEL 8192;
  57. FINSI;
  58.  
  59. * RUN 4
  60. SI (RUN EGA 4);
  61. OPTI ACQU '.\exp\24time.txt';
  62. ACQU TIME*LISTREEL 8192;
  63. OPTI ACQU '.\exp\24axtab.txt';
  64. ACQU AXTAB*LISTREEL 8192;
  65. FINSI;
  66.  
  67. **********************************************
  68. * REMPLISSAGE DE LA TABLE D ENTREE A @A1DDL
  69. **********************************************
  70. * Parametres dynamiques
  71. *TABDYN . 1  = M;     -- masse de la structure
  72. *TABDYN . 2  = BETA;  -- coefficient beta
  73. *TABDYN . 3  = GAMMA; -- coefficient gamma
  74. *TABDYN . 4  = FPLAS; -- effort de plastification
  75. *TABDYN . 5  = KA;    -- raideur de l'acier
  76. *TABDYN . 6  = KB;    -- raideur du beton
  77. *TABDYN . 7  = ACTU   -- type d actualisation
  78. *                        --> 1 : 
  79. *                        --> 2 : 
  80. *                        --> 3 : 
  81. *TABDYN . 8  = TIME   -- liste de temps
  82. *TABDYN . 9  = AXTAB  -- liste d acceleration en pied de maquette
  83. *TABDYN . 10 = XI0    -- taux d amortissement initial
  84. *TABDYN . 11 = XMIN   -- taux d amrotissement minimal
  85. *TABDYN . 12 = AMMAX  -- taux d amortissement maximal 
  86. *TABDYN . 13 = NC     -- indicateur sur le type d actualisation
  87. *                        --> 0 : 
  88. *                        --> 1 : 
  89. *TABDYN . 14 = DPLUS  -- endommagement positif
  90. *TABDYN . 15 = DMOIN  -- endommagement negatif
  91. *TABDYN . 16 = MAXDP  -- maximum deplacement positif
  92. *TABDYN . 17 = MAXDM  -- maximum deplacement negatif
  93. *TABDYN . 18 = AOLD   -- taux d amortissement au premier pas
  94.  
  95. TABDYN = TABLE;
  96. TABDYN . 1  = 2980.0;
  97. TABDYN . 2  = 0.25;
  98. TABDYN . 3  = 0.50;
  99. TABDYN . 4  = 35000.0;
  100. TABDYN . 5  = 729.7855;
  101. TABDYN . 6  = 2.8213E6;
  102. TABDYN . 7  = 1;
  103. TABDYN . 8  = TIME;
  104. TABDYN . 9  = AXTAB;
  105. TABDYN . 10 = 0.05;
  106. TABDYN . 11 = 0.02;
  107. TABDYN . 12 = 0.30;
  108. TABDYN . 13 = 0.0;
  109. TABDYN . 14 = 0.0;
  110. TABDYN . 15 = 0.0;
  111. TABDYN . 16 = 0.0;
  112. TABDYN . 17 = 0.0;
  113. TABDYN . 18 = 0.0;
  114.  
  115. **********************************************
  116. * APPEL A @A1DDL
  117. **********************************************
  118. EVDR EVVR EVAR EVXI TAB3 = @A1DDL TABDYN;
  119.  
  120. **********************************************
  121. * ECRITURES EN SORTIES
  122. **********************************************
  123. *@EXCEL1 EVDR (CHAI '.\evdr' RUN '.txt');
  124. *@EXCEL1 EVVR (CHAI '.\evvr' RUN '.txt');
  125. *@EXCEL1 EVAR (CHAI '.\evar' RUN '.txt');
  126. *@EXCEL1 EVXI (CHAI '.\evxi' RUN '.txt');
  127.  
  128. LVDRA = extr (extr EVDR ordo) 500;
  129. LVVRA = extr (extr EVVR ordo) 500;
  130. LVARA = extr (extr EVAR ordo) 500;
  131. LVXIA = extr (extr EVXI ordo) 500;
  132.  
  133. err1 = abs (1.56839E-05 - LVDRA);
  134. err2 = abs (-7.66832E-03 - LVVRA);
  135. err3 = abs (-2.50373E-02 - LVARA);
  136. err4 = abs (2.10603E-02 - LVXIA);
  137.  
  138. si (> err1 1.0E-5);
  139. erreur(5);
  140. finsi;
  141.  
  142. si (> err2 1.0E-5);
  143. erreur(5);
  144. finsi;
  145.  
  146. si (> err3 1.0E-5);
  147. erreur(5);
  148. finsi;
  149.  
  150. si (> err4 1.0E-5);
  151. erreur(5);
  152. finsi;
  153.  
  154. FIN;
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  

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