Télécharger la plaquette Cast3M


fr en

Télécharger usadac.procedur

Retour à la liste

Numérotation des lignes : 

  1. * USADAC    PROCEDUR  MB234859  24/02/06    21:15:02     11833          
  2. ************************************************************************
  3. *                         PROCEDURE USADAC
  4. *
  5. * Procedure qui approxime la largeur de la zone de contact L en fonction 
  6. * du nombre de cycle N via une fonction de la forme :
  7. *
  8. *                    L(N) = L0 + [L1 '*' N '**' L2]
  9. *                dL(N)/dN = L1 '*' L2 '*' N '**' (L2 '-' 1)
  10. *
  11. * Cette procedure est celle a fournir a LEVM pour determiner les valeurs 
  12. * des parametres L1 et L2 au cours de la resolution.
  13. *
  14. * Lorsque DlCyc  = 0, elle retourne dL(N)/dN
  15. *
  16. * Lorsque DlCyc /= 0, elle calcule  : 
  17. *                  L(N) '-' L0 = L1 '*' N '**' L2
  18. *                      dL/dL1  = N '**' L2
  19. *                      dL/dL2  = L1 '*' LOG N '*' N '**' L2
  20. *
  21. * Il s'agit du schema developpe par Q.Caradec lors de sa these.
  22. *
  23. * Entrees :
  24. *   T1    : Table de PASAPAS
  25. *
  26. * Appelee par USACCE et utilisee par LEVM
  27. *
  28. ************************************************************************
  29. 'DEBP' USADAC lCyc*'LISTREEL' lLambda*'LISTREEL';
  30. *
  31.   DlCyc   = 'DIME' lCyc ;
  32.   lambda1 = 'EXTR' lLambda 1;
  33.   lambda2 = 'EXTR' lLambda 2;
  34. *-----------------------------------------------------------------------
  35.   'SI' ('EGA' DlCyc 0) ;
  36.     Llcont  = 0. ;
  37.     Ncycle  = 'EXTR' lLambda 3;
  38.     Ldlcont = lambda1 '*' lambda2 '*' (NCycle '**' (lambda2 '-' 1.));
  39.   'SINON' ;
  40.     Llcont = 'PROG' ; Ldlcont = 'PROG' ;
  41.     'REPE' i DlCyc ;
  42.       Cyc_i = 'EXTR' lCyc &i;
  43.       dlcl1 = Cyc_i '**' lambda2 ;
  44.       lcont = lambda1 '*' dlcl1 ;
  45.       dlcl2 = lambda1 '*' (LOG Cyc_i) '*' dlcl1 ;
  46. *     
  47.       Llcont  = Llcont  'ET' lcont ;
  48.       Ldlcont = Ldlcont 'ET' dlcl1 'ET' dlcl2 ;
  49.     'FIN' i;
  50.   'FINSI' ;
  51. *-----------------------------------------------------------------------
  52. 'FINP' Llcont Ldlcont ;
  53.  
  54.  

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