Télécharger acti3.procedur

Retour à la liste

Numérotation des lignes :

  1. * ACTI3 PROCEDUR NOBODY 92/01/01 00:00:00 001
  2. * appele par INCREME pour extrapoler en utilisant une estimation
  3. * de la matrice tangente
  4.  
  5. 'DEBPROC' ACTI3 CHD3*CHPOINT CHD2*CHPOINT CHD1*CHPOINT CHD0*CHPOINT
  6. FO3*CHPOINT FO2*CHPOINT FO1*CHPOINT FO0*CHPOINT ;
  7.  
  8. L = 0 ; M = 0 ; N = 0 ;
  9.  
  10. 'REPETER' LABEUR 1;
  11.  
  12. LMT = 'MOTS' 'FX' 'FY' 'FZ' 'FR' 'FT' 'MX' 'MY' 'MZ' 'MR' 'MT';
  13.  
  14. F1_F0 = FO1 - FO0 ; C0_C1 = CHD0 - CHD1 ;
  15. F2_F0 = FO2 - FO0 ; C0_C2 = CHD0 - CHD2 ;
  16. F3_F0 = FO3 - FO0 ; C0_C3 = CHD0 - CHD3 ;
  17.  
  18. F12 = 'XTX' F1_F0 ; F22 = 'XTX' F2_F0; F32 = 'XTX' F3_F0;
  19.  
  20. F1F2= F1_F0 'XTY' F2_F0 LMT LMT ;
  21. F1F3= F1_F0 'XTY' F3_F0 LMT LMT ;
  22. F2F3= F2_F0 'XTY' F3_F0 LMT LMT ;
  23. F0F1= FO0 'XTY' F1_F0 LMT LMT ;
  24. F0F2= FO0 'XTY' F2_F0 LMT LMT ;
  25. F0F3= FO0 'XTY' F3_F0 LMT LMT ;
  26.  
  27. A1=F12 ; B1=F1F2 ; C1=F1F3 ; J1 = (-1)*F0F1;
  28. D1=F1F2; E1=F22 ; F1=F2F3 ; K1 = (-1)*F0F2;
  29. G1=F1F3; H1=F2F3 ; I1=F32 ; L1 = (-1)*F0F3;
  30.  
  31. 'SI' (ABS(I1) < 1E-30) ; 'QUITTER' LABEUR; 'FINSI' ;
  32.  
  33. A2=A1-(G1*C1/I1); B2=B1-(H1*C1/I1); J2=J1-(L1*C1/I1);
  34. D2=D1-(G1*F1/I1); E2=E1-(H1*F1/I1); K2=K1-(L1*F1/I1);
  35.  
  36. 'SI' (ABS(E2) < 1E-30) ; 'QUITTER' LABEUR; 'FINSI' ;
  37.  
  38. A3=A2-(D2*B2/E2); J3=J2-(K2*B2/E2);
  39.  
  40. 'SI' (ABS(A3) < 1E-30) ; 'QUITTER' LABEUR; 'FINSI' ;
  41.  
  42. L=J3/A3;
  43.  
  44. B4=B2; J4=J2-(A2*L);
  45.  
  46. 'SI' (ABS(B4) < 1E-30) ; 'QUITTER' LABEUR; 'FINSI' ;
  47.  
  48. M=J4/B4;
  49.  
  50. C5=C1; J5=J1-(A1*L)-(B1*M);
  51.  
  52. 'SI' (ABS(C5) < 1E-30) ; 'QUITTER' LABEUR; 'FINSI' ;
  53.  
  54. N=J5/C5;
  55.  
  56. 'FIN' LABEUR;
  57.  
  58. 'MESS' 'coefficients acceleration ' (L*-1) (M*-1) (N*-1) ;
  59.  
  60. TEMP1 = L * C0_C1 ;
  61. TEMP2 = M * C0_C2 ;
  62. TEMP3 = N * C0_C3 ;
  63. TEMP4 = TEMP1 + TEMP2 ;
  64. RES = TEMP4 + TEMP3 ;
  65.  
  66. 'DETRUIRE' TEMP1 ; 'DETRUIRE' TEMP2 ; 'DETRUIRE' TEMP3 ;
  67. 'DETRUIRE' TEMP4 ;
  68. 'DETRUIRE' F1_F0 ; 'DETRUIRE' F2_F0 ; 'DETRUIRE' F3_F0 ;
  69. 'DETRUIRE' C0_C1 ; 'DETRUIRE' C0_C2 ; 'DETRUIRE' C0_C3 ;
  70.  
  71. 'FINPROC' RES ;
  72.  

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