Télécharger la plaquette Cast3M


fr en

Télécharger piloindi.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PILOINDI  PROCEDUR  HO261890  20/12/22    21:15:08     10832          
  2. 'DEBP' PILOINDI PRECED*'TABLE' U*'CHPOINT' DU*'CHPOINT' DUI*'CHPOINT' DUII*'CHPOINT' DTAU*'FLOTTANT' ;
  3.  
  4. WTAB=PRECED.'WTABLE';
  5.  
  6. REGION = WTAB . 'ZONE_DE_PILOTAGE' ;
  7. 'SI' ('EGA' REGION FAUX) ; 
  8.    MODTOT =  WTAB.'MO_TOT'  ;
  9. 'SINON' ;
  10.    MODTOT = REDU WTAB.'MO_TOT' REGION ;
  11. 'FINSI' ;
  12.  
  13. lis1 = Extr Modtot 'DEFO';
  14. lis2 = Extr Modtot 'CONT';
  15.  
  16. *   --- equation de pilotage indirecte --------
  17. CHATAU = 'MANU' 'CHML' MODTOT 'SCAL' DTAU 'STRESSES';
  18.  
  19. EPS_TOT = 'EPSI' 'LINE' MODTOT dU;
  20. def_I   = 'EPSI' 'LINE' MODTOT duI ;
  21. CH_A_0  = EPS_TOT + def_I;
  22. CH_A_1  = 'EPSI' 'LINE' MODTOT duII ;
  23.  
  24. EPS_TOT1 = 'EPSI' 'LINE' MODTOT u;
  25. dir1  = CHAN 'TYPE' EPS_TOT1 'CONTRAINTES';
  26. dir1  = CHAN 'COMP' lis1 lis2 dir1;
  27. norm1 = ENER  MODTOT dir1  EPS_TOT1  ;
  28. norm1 = norm1**0.5;
  29.  
  30. petit1 = 1.E-20;
  31. normnul = ((maxi norm1) < petit1);
  32.  
  33. SI (normnul);
  34.   REPE BC1 (DIME(lis2));
  35.     SI (&bc1 EGA 1);
  36.       chm =         'MANU' 'CHML' MODTOT (EXTR lis2 (&BC1)) 1.0 'STRESSES';
  37.     SINON;
  38.       chm = chm et ('MANU' 'CHML' MODTOT (EXTR lis2 (&BC1)) 1.0 'STRESSES');
  39.     FINSI;
  40.   FIN BC1;
  41. SINON;
  42.   REPE BC1 (DIME(lis2));
  43.     SI (&bc1 EGA 1);
  44.       aux1 = 'CHAN' 'COMP' (EXTR lis2 (&BC1)) norm1;
  45.       nm1 = aux1;
  46.     SINON;
  47.       aux1 = 'CHAN' 'COMP' (EXTR lis2 (&BC1)) norm1;
  48.       nm1 = nm1 et aux1;
  49.     FINSI;
  50.   FIN BC1;
  51.   chm = dir1*(nm1**(-1.0));
  52. FINSI;
  53.  
  54. CH_A_1  = ENER  MODTOT chm CH_A_1 ;
  55. CH_A_0  = ENER  MODTOT chm CH_A_0 ;
  56.  
  57. CH_ETA  = (CHATAU - CH_A_0)*(CH_A_1**(-1.0));
  58.  
  59. MA1 = 'MASQ' CH_A_1 'SUPERIEUR' 0. ;
  60. CH_LAMB_MAX = (MA1*CH_ETA);
  61. MA2 = 'MASQ' CH_A_1 'EGINFE' 0. ;
  62. CH_LAMB_MIN = (MA2*CH_ETA);
  63.  
  64. big_value = 1.0E20;
  65. CH_LAMB_MAX = CH_LAMB_MAX + (big_value*MA2);
  66. CH_LAMB_MIN = CH_LAMB_MIN - (big_value*MA1);
  67. LIM_INF = 'MAXI' CH_LAMB_MIN 'AVEC' (MOTS 'SCAL');
  68. LIM_SUP = 'MINI' CH_LAMB_MAX 'AVEC' (MOTS 'SCAL');
  69.  
  70. D_ETA = LIM_SUP;
  71.  
  72. 'FINPROC' D_ETA ;
  73.  

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