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

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