Procedures | Cast3M

   1 : * @SIGREF   PROCEDUR  ALANIER   95/04/03    21:16:09     1603           
   2 : 'DEBP' @SIGREF P*MCHAML M*MCHAML W*MCHAML A*MCHAML MO*MMODEL ;
   3 : 
   4 :   CHP01 = P 'MASQ' 'INFERIEUR'  0.0001 ;
   5 :   CHP02 = P 'MASQ' 'SUPERIEUR' -0.0001 ;
   6 :   CHP0  = CHP01 * CHP02 ;
   7 : 
   8 :   CHM01 = M 'MASQ' 'INFERIEUR'  0.0001 ;
   9 :   CHM02 = M 'MASQ' 'SUPERIEUR' -0.0001 ;
  10 :   CHM0  = CHM01 * CHM02 ;
  11 : 
  12 :   CHMP0 = CHP0 * CHM0 ;
  13 :   CHC = CHMP0 'MASQ' 'INFERIEUR' 0.5 ;
  14 : 
  15 : * CAS GENERAL
  16 : 
  17 :   P_M = P + CHP0 ;M_M = M + CHM0 ;
  18 :  
  19 :   CH1 = 'MANU' 'CHML' MO 'SCAL' 1. 'TYPE' 'SCALAIRE' 'STRESSES' ;
  20 :   PSM = (P_M * (M_M ** -1)) ;
  21 : 
  22 :   B = (PSM * (((((W / 2.) - A) ** 2) / 2.)
  23 :    + ((W ** 2) / 8.))) + A ;
  24 :   D1 = (CH1 + ((CH1 + (PSM * B)) ** 0.5)) * (PSM ** -1) ;
  25 :   D2 = (CH1 - ((CH1 + (PSM * B)) ** 0.5)) * (PSM ** -1) ;
  26 :   S1 = (-1.) * (P_M * ((A + (2. * D1)) ** -1)) ;
  27 :   S2 = (-1.) * (P_M * ((A + (2. * D2)) ** -1)) ;
  28 : 
  29 :   BS = ((W / 2.) - A) ;
  30 : 
  31 :   DD1 = D1 'MASQ' 'INFERIEUR' (BS - ('ABS' (BS * 0.00001))) ;
  32 :   DDD1 = D1 'MASQ' 'SUPERIEUR' ((W / (-2.)) * 0.99999) ;
  33 : 
  34 :   SS1 = S1 * (DDD1 * DD1) ;
  35 : 
  36 :   DD2 = D2 'MASQ' 'INFERIEUR' (BS - ('ABS' (BS * 0.00001))) ;
  37 :   DDD2 = D2 'MASQ' 'SUPERIEUR' ((W / (-2.)) * 0.99999) ;
  38 : 
  39 :   SS2 = S2 * (DDD2 * DD2) ;
  40 : 
  41 :   DDD = (DDD1 * DD1) + (DDD2 * DD2) ;
  42 :   'SI' (('ABS' (('MAXI' DDD) - ('MINI' DDD))) > 0.99) ;
  43 :     'MESS' 'PROBLEME DANS CALCUL DE SIGMA REF' ;
  44 :   'FINSI' ;
  45 : 
  46 :   SSG = SS1 + SS2 ;
  47 : 
  48 : * CAS PARTICULIER OU P EST TRES PETIT
  49 : 
  50 :   SSP0 = (4. * M) * (((W ** 2) - (2. * A * W) + (A ** 2)) ** -1) ;
  51 : 
  52 : * CAS PARTICULIER OU M EST TRES PETIT
  53 : 
  54 :   SSM0 = P * ((W - A) ** -1) ;
  55 : 
  56 : * SOLUTION :
  57 : 
  58 :   SS = (SSP0 * CHP0) + (SSM0 * CHM0) + (SSG * CHC) ;
  59 : 
  60 : 
  61 : 'FINP' SS ;