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.  
  74.  
  75.  

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