Télécharger dyna_nl4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * FICHIER : DYNA_NL4.DGIBI
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************
  5.  
  6. *******************************************************
  7. * *
  8. * CALCUL DU REBOND D'UNE BARRE *
  9. * AVEC PASAPAS ET RELA *
  10. * *
  11. * D. COMBESCURE AOUT 2006 *
  12. *******************************************************
  13. *GRAPH = VRAI;
  14. *COMPLET = VRAI;
  15. GRAPH = FAUX;
  16. COMPLET = FAUX;
  17. *
  18. OPTION ECHO 1 ;
  19. OPTIO DIME 2 ELEM QUA8 MODE PLAN CONT;
  20. *
  21. * MAILLAGE
  22. *
  23. NL = 10;
  24. H1 = 1.D-1;
  25. L1 = 0.5;
  26. P1 = (0.) (0.5*H1);
  27. P2 = (0.) ((-0.5)*H1);
  28. VH1 = L1 0.;
  29. P3 = L1 (0.5*H1);
  30. P4 = L1 ((-0.5)*H1);
  31. P3B = L1 (0.5 * H1);
  32. P4B = L1 ((-0.5)*H1);
  33. P5 = ((-1.)*L1) (0.5 * H1);
  34. P6 = ((-1.)*L1) ((-0.5)*H1);
  35. LIG1 = D 1 P1 P2;
  36. LIG2 = D 1 P3 P4;
  37. LIG2B = D 1 P3B P4B;
  38. LIG3 = D 1 P5 P6;
  39. S1 = (LIG1 TRAN NL VH1) ET (LIG1 TRAN NL ((-1.)*VH1));
  40.  
  41. * S1 = (LIG1 TRAN NL VH1);
  42. LIGH = (D NL P5 P1) ET (D NL P1 P3);
  43. ELIM 0.000001 (LIG1 ET LIG2 ET S1 ET LIG3 ET LIGH);
  44. *
  45. * MODELE
  46. *
  47. MOD1 = MODELE S1 MECANIQUE ELASTIQUE ;
  48. *
  49. * PROPRIETES DE MATERIAUX
  50. *
  51. RO = 7800.;
  52. MAT1 = MATER MOD1 YOUN 200.E6 NU 0.0 RHO RO ;
  53. *
  54. * CONDITIONS AUX LIMITES
  55. *
  56. CL1 = (BLOQ S1 UY) ;
  57. CL2 = RELA 'MAXI' UX LIG2 - 'UX' LIG2B;
  58. CL3 = BLOQ DEPL LIG2B;
  59. DEPIM2 =DEPI CL2 0.01;
  60. *
  61. * MATRICE DE RIGIDITE
  62. *
  63. RIG1 = RIGI MOD1 MAT1 ;
  64. *
  65. * MATRICE DE MASSE
  66. *
  67. MAS1 = MASSE MOD1 MAT1 ;
  68. *
  69. * CHARGEMENT
  70. *
  71. DDT = 2.0D-5;
  72. LI1 = (PROG 0. DDT (2.*DDT) 10.);
  73. LI2 = (PROG 1. 1. 1. 1.);
  74. EV = EVOL MANU 'T' LI1 'F(T)' LI2 ;
  75. *DESS EV YBORD 0. 2. XBORD 0. 1.0;
  76. F1 = MAS1*(MANU CHPO S1 1 UX (9.81));
  77. F2 = MAS1*(MANU CHPO S1 1 UX (9.81));
  78. CHA1 = CHAR (F1 ET DEPIM2) EV ;
  79. CHA2 = CHAR (F2 ET DEPIM2) EV ;
  80. *CHA1 = CHAR (DEPIM2) EV ;
  81. *CHA2 = CHAR (DEPIM2) EV ;
  82. *
  83. * PREPARATION DE LA TABLE DE DONNEES POUR PASAPAS
  84. *
  85.  
  86. SI (EGA COMPLET VRAI) ;
  87. TIME = PROG 0. PAS 1.E-3 50.D-2 ;
  88. SINON ;
  89. TIME = PROG 0. PAS 1.E-3 1.D-2 ;
  90. FINSI ;
  91.  
  92.  
  93.  
  94.  
  95. *
  96. TAB1 = TABLE;
  97. TAB1 . 'MODELE' = MOD1;
  98. TAB1 . 'CARACTERISTIQUES' = MAT1;
  99. TAB1 . 'BLOCAGES_MECANIQUES' = CL1 ET CL2 ET CL3;
  100. TAB1 . 'CHARGEMENT' = CHA1;
  101. *TAB1 . 'AMORTISSEMENT' = 0.01*(RIGI MOD1 MAT1);
  102. TAB1 . 'TEMPS_CALCULES' = TIME;
  103. TAB1 . 'DYNAMIQUE' = VRAI;
  104. TAB1 . 'CONVERGENCE_FORCEE' = FAUX;
  105. PASAPAS TAB1;
  106. *
  107. * EXTRAIRE LES OBJETS D'EVOLUTION DE DEPLACEMENTS
  108. *
  109. LT = VIDE 'LISTREEL' ;
  110. LUX1 = VIDE 'LISTREEL' ;
  111. LUX2 = VIDE 'LISTREEL' ;
  112. LUX3 = VIDE 'LISTREEL' ;
  113. LVX1 = VIDE 'LISTREEL' ;
  114. LVX2 = VIDE 'LISTREEL' ;
  115. LVX3 = VIDE 'LISTREEL' ;
  116. LEPOT = VIDE 'LISTREEL' ;
  117. LENEC = VIDE 'LISTREEL' ;
  118. LEELA = VIDE 'LISTREEL' ;
  119. LEREA = VIDE 'LISTREEL' ;
  120. I = 1 ;
  121. M1 = P5 ;M2 = P1 ; M3 = P3 ;
  122. TEMPS = INDEX ( TAB1.DEPLACEMENTS);
  123. NBB = (DIME TEMPS) - 2;
  124.  
  125. REPETER BB1 NBB ;
  126. I = I + 1 ;
  127. LT = LT ET (PROG (TAB1.'TEMPS'. (TEMPS. I))) ;
  128. DEP1 = TAB1 . 'DEPLACEMENTS' . (TEMPS. I) ;
  129. EPS1 = EPSI DEP1 MOD1 MAT1;
  130. VIT1 = TAB1 . 'VITESSES' . (TEMPS. I) ;
  131. CON1 = TAB1 . 'CONTRAINTES' . (TEMPS. I) ;
  132. REA1 = REAC DEP1 CL3;
  133. ENE1 = ENER MOD1 EPS1 CON1;
  134. LUX1 = LUX1 ET( PROG (EXTR DEP1 UX P1) ) ;
  135. LUX2 = LUX2 ET( PROG (EXTR DEP1 UX P2) ) ;
  136. LUX3 = LUX3 ET( PROG (EXTR DEP1 UX P3) ) ;
  137. LVX1 = LVX1 ET( PROG (EXTR VIT1 UX P1) ) ;
  138. LVX2 = LVX2 ET( PROG (EXTR VIT1 UX P2) ) ;
  139. LVX3 = LVX3 ET( PROG (EXTR VIT1 UX P3) ) ;
  140.  
  141. * LEELA = LEELA ET (PROG (INTG ENE1 MOD1));
  142. LEELA = LEELA ET
  143. (PROG ((0.5*('XTY' DEP1 (RIG1*DEP1)
  144. ('MOTS' UX UY)
  145. ('MOTS' FX FY))) ));
  146. LENEC = LENEC ET
  147. (PROG ((0.5*(XTY VIT1 (MAS1*VIT1)
  148. (MOTS UX UY) (MOTS FX FY))) ));
  149. LEPOT = LEPOT ET
  150. (PROG (XTY DEP1 (MAS1*(MANU CHPO S1 1 UX (9.81)))
  151. (MOTS UX UY) (MOTS FX FY)));
  152. LEREA = LEREA ET (PROG (EXTR (RESU REA1) P1 FX));
  153. FIN BB1 ;
  154.  
  155. EVU1 = EVOL MANU 'TIME' LT 'UX ' ((-1.)*LUX1) ;
  156. EVU2 = EVOL ROUGE MANU 'TIME' LT 'UX ' ((-1.)*LUX2);
  157. EVU3 = EVOL BLEU MANU 'TIME' LT 'UX ' ((-1.)*LUX3) ;
  158. EVREA = EVOL BLEU MANU 'TIME' LT 'UX ' (LEREA) ;
  159.  
  160. EVEC = EVOL BLEU MANU 'TIME' LT 'EC' LENEC ;
  161. EVPOT = EVOL ROUGE MANU 'TIME' LT 'EP' LEPOT ;
  162. EVEL = EVOL MANU 'TIME' LT 'EP' LEELA ;
  163. EVBILAN = EVOL MANU 'TIME' LT (((-1.)*LEPOT) + LEELA + LENEC);
  164.  
  165. EVV1 = EVOL MANU 'TIME' LT 'UX ' ((-1.)*LVX1) ;
  166. EVV2 = EVOL ROUGE MANU 'TIME' LT 'UX ' ((-1.)*LVX2) ;
  167. EVV3 = EVOL BLEU MANU 'TIME' LT 'UX ' ((-1.)*LVX3) ;
  168.  
  169. *
  170. * --- TRACER LES COURBES
  171. *
  172. SI GRAPH ;
  173. TA = TABLE;
  174. TA.1 = 'MARQ CARR';
  175. DESS (EVU1 ET EVU2 ET EVU3) TITRE 'DEPLACEMENT';
  176. DESS (EVV1 ET EVV2 ET EVV3) TITRE 'VITESSE';
  177. DESS (EVU3) XBORD 0. 0.20 TITRE 'DEPLACEMENT';
  178. DESS (EVV3) XBORD 0. 0.20 TITRE 'VITESSE';
  179. TA = TABLE;
  180. TA.1 = 'TIRR';
  181. TA.2= 'TIRC';
  182. DESS (EVPOT ET EVEC ET EVEL ET EVBILAN) TA XBORD 0. 0.20
  183. YBORD 0. 90. TITRE 'ENERGIE';
  184. FINSI ;
  185.  
  186. *
  187. TEST = (MAXI (EXTR EVBILAN ORDO) ABS)
  188. / (MAXI (EXTR EVPOT ORDO) ABS);
  189. *
  190. SI (TEST > 1.D-6) ;
  191. ERR 4;
  192. FINSI ;
  193.  
  194. FIN;
  195.  
  196.  
  197.  
  198.  
  199.  

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