Procedures | Cast3M

   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 ;