Procedures | Cast3M

   1 : * IDBHT     PROCEDUR  JEANV     94/10/10    21:31:41     1285           
   2 : 'DEBPROC' IDBHT  M1*'MCHAML'  MO1*'MMODEL'  BHTMX*'CHPOINT'
   3 :                  T2*'CHPOINT' ;
   4 : 
   5 : * Identificazione dei modelli di cemento ad alta temperatura
   6 : IFER= FAUX;
   7 : 
   8 : *------------------- modèle drucker-parfait ------------------------
   9 : 
  10 : 'SI' ('EXISTE' MO1 'MATE' 'DRUCKER_PARFAIT');
  11 :    LCO1 = 'EXTR' MO1 'COMP';
  12 :    DIM1 = 'DIME' LCO1;
  13 :    ISU1 = FAUX;
  14 :    ISU2 = FAUX;
  15 :    'REPE' B2BHT DIM1;
  16 :       ICOMP=&B2BHT;
  17 :       CC=LCO1.ICOMP;
  18 :      'SI' ('OU'('OU'('OU'('EGA' CC 'LTR'  )
  19 :                          ('EGA' CC 'LCS' ))
  20 :                          ('EGA' CC 'NU'  ))
  21 :                          ('EGA' CC 'YOUN'));
  22 :         'SI' ISU1;
  23 :            MA1 = MA1 'ET' ('EXCO' M1 CC CC);
  24 :         'SINON';
  25 :            MA1 = 'EXCO' M1 CC CC;
  26 :            ISU1 = VRAI;
  27 :         'FINSI';
  28 :      'SINON';
  29 :         'SI' ISU2;
  30 :            MA2 = MA2 'ET' ('EXCO' M1 CC CC);
  31 :         'SINON';
  32 :            MA2 = 'EXCO' M1 CC CC;
  33 :            ISU2 = VRAI;
  34 :         'FINSI';
  35 :      'FINSI';
  36 :    'FIN' B2BHT;
  37 :    'SI' ISU1;
  38 :       CCC1 = 'VARI' MO1 MA1 BHTMX;
  39 :    'FINSI';
  40 :    'SI' ISU2;
  41 :       CCC2 = 'VARI' MO1 MA2 T2;
  42 :    'FINSI';
  43 :    'SI' (ISU1 'ET' ISU2);
  44 :       CCC = CCC1 'ET' CCC2;
  45 :    'SINON';
  46 :       'SI' ISU1;
  47 :          CCC = CCC1;
  48 :       'FINSI';
  49 :       'SI' ISU2;
  50 :          CCC = CCC2;
  51 :       'FINSI';
  52 :    'FINSI';
  53 :    MI  = 'MATR' MO1 CC CCC ;
  54 : 
  55 : *------------------- modele drucker_prager ---------------------
  56 : 
  57 : 'SINON';
  58 :    'SI'( 'EXISTE' MO1 'MATE' 'DRUCKER_PRAGER');
  59 :      'ERREUR' 'Le modele DRUCKER_PRAGER H. T. est indisponible.';
  60 :       IFER =  VRAI;
  61 :      'QUITTER' IDBHT;
  62 :    'SINON';
  63 : 
  64 : *------------------------ modèle mazars -------------------------
  65 : 
  66 :      'SI'('EXISTE' MO1 'MATE' 'MAZARS');
  67 :         LCO1 = 'EXTR' MO1 'COMP';
  68 :         DIM1 = 'DIME' LCO1;
  69 :         ISU1 = FAUX;
  70 :         ISU2 = FAUX;
  71 :         TMAZ=TABLE;
  72 :        'REPE' B3BHT DIM1;
  73 :           ICOMP = &B3BHT;
  74 :           CC=LCO1.ICOMP;
  75 :           'SI' ('OU'('OU'('OU'('OU'('OU'('OU'('EGA' CC 'LTR'  )
  76 :                                             ('EGA' CC 'LCS' ))
  77 :                                             ('EGA' CC 'NU'  ))
  78 :                                             ('EGA' CC 'YOUN'))
  79 :                                             ('EGA' CC 'EMAX'))
  80 :                                             ('EGA' CC 'ATRA'))
  81 :                                             ('EGA' CC 'BTRA'));
  82 : 
  83 :              MA1 = 'EXCO' M1 CC CC;
  84 :              TMAZ.CC = 'VARI' MO1 MA1 BHTMX;
  85 :            'SINON';
  86 :              MA1 = 'EXCO' M1 CC CC;
  87 :              TMAZ.CC = 'VARI' MO1 MA1   T2;
  88 :            'FINSI';
  89 :        'FIN' B3BHT;
  90 : *
  91 :         ED0MAZ = TMAZ.'LTR' * ( TMAZ.'YOUN' **(-1.) ) ;
  92 :         BCMAZ=(-1./(2.**.5))*((TMAZ.'NU' * TMAZ.'EMAX')**(-1.));
  93 :         BCED0MAZ= BCMAZ * ED0MAZ;
  94 :         ACMAZ1=((2.**.5)* TMAZ.'LCS' * TMAZ.'NU' *
  95 :                (TMAZ.'YOUN' **(-1.)  )) + ED0MAZ;
  96 :         ACMAZ2= BCMAZ * ( BCED0MAZ - ('EXP'(
  97 :           ('MANU' 'CHML' MO1 'SCAL' -1. 'TYPE' 'SCALAIRE' 'RIGIDITE')
  98 :           +  BCED0MAZ)));
  99 :         ACMAZ= ACMAZ1 * ACMAZ2;
 100 :         MI1 = 'MATR' MO1
 101 :               'YOUN' ('EXCO' TMAZ.'YOUN' 'SCAL' 'NOID' 'YOUN')
 102 :               'NU'   ('EXCO' TMAZ.'NU'   'SCAL' 'NOID' 'NU'  )
 103 :               'ALPH' ('EXCO' TMAZ.'ALPH' 'SCAL' 'NOID' 'ALPH');
 104 :         MI2 = 'MATR' MO1
 105 :               'ATRA' ('EXCO' TMAZ.'ATRA'  'SCAL' 'NOID' 'ATRA' )
 106 :               'BTRA' ('EXCO' TMAZ.'BTRA'  'SCAL' 'NOID' 'BTRA' )
 107 :               'BCOM' ('EXCO' BCMAZ        'SCAL' 'NOID' 'BCOM' )
 108 :               'ACOM' ('EXCO' ACMAZ        'SCAL' 'NOID' 'ACOM' )
 109 :               'KTR0' ('EXCO' ED0MAZ       'SCAL' 'NOID' 'KTR0' )
 110 :               'BETA' ('EXCO' TMAZ.'BETA'  'SCAL' 'NOID' 'BETA' );
 111 : *
 112 :         MI = MI1 + MI2 ;
 113 :         'SI' ('EXIST' TMAZ 'RHO' );
 114 :            MI= MI + ('MATR' MO1  'RHO'
 115 :                    ('EXCO' TMAZ.'RHO'   'SCAL' 'NOID' 'RHO'  ));
 116 :         'FINSI';
 117 : *
 118 :      'SINON';
 119 :         IFER=VRAI;
 120 :      'ERREUR' 'Le modele donne ne correspond pas a ceux attendus';
 121 :      'FINSI';
 122 :    'FINSI';
 123 : 'FINSI';
 124 : 'FINPROC' MI IFER;