Télécharger la plaquette Cast3M


fr en

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