Procedures | Cast3M

   1 : * POSTDDI1  PROCEDUR  CHAT      96/08/05    21:25:21     2238           
   2 : 'DEBPROC' POSTDDI1 var1*MCHAML alpha*FLOTTANT mod0*MMODEL;
   3 : *
   4 : ****************************************************************
   5 : *     Calcul de CHI(Xv) = 1.5*ALPHA*XV1 + (1-ALPHA)*J2(Xv)     *
   6 : *       Avec XV1 = Contrainte principale maximale de Xv        *
   7 : ****************************************************************
   8 : *
   9 : mot1 = 'VALE' 'MODE';
  10 : 'SI' ('EGA' mot1 'AXIS');
  11 : xvrr = ('EXCO' var1 'Y1RR' 'SMRR') + ('EXCO' var1 'Y2RR' 'SMRR');
  12 : xvzz = ('EXCO' var1 'Y1ZZ' 'SMZZ') + ('EXCO' var1 'Y2ZZ' 'SMZZ');
  13 : xvtt = ('EXCO' var1 'Y1TT' 'SMTT') + ('EXCO' var1 'Y2TT' 'SMTT');
  14 : xvrz = ('EXCO' var1 'Y1RZ' 'SMRZ') + ('EXCO' var1 'Y2RZ' 'SMRZ');
  15 : xv = xvrr 'ET' xvzz 'ET' xvtt 'ET' xvrz ;
  16 : 'SINON';
  17 : 'SI' ('EGA' mot1 'TRID');
  18 : xvxx = ('EXCO' var1 'Y1XX' 'SMXX') + ('EXCO' var1 'Y2XX' 'SMXX');
  19 : xvyy = ('EXCO' var1 'Y1YY' 'SMYY') + ('EXCO' var1 'Y2YY' 'SMYY');
  20 : xvzz = ('EXCO' var1 'Y1ZZ' 'SMZZ') + ('EXCO' var1 'Y2ZZ' 'SMZZ');
  21 : xvxy = ('EXCO' var1 'Y1XY' 'SMXY') + ('EXCO' var1 'Y2XY' 'SMXY');
  22 : xvyz = ('EXCO' var1 'Y1YZ' 'SMYZ') + ('EXCO' var1 'Y2YZ' 'SMYZ');
  23 : xvxz = ('EXCO' var1 'Y1ZX' 'SMXZ') + ('EXCO' var1 'Y2ZX' 'SMXZ');
  24 : xv = xvxx 'ET' xvyy 'ET' xvzz 'ET' xvxy 'ET' xvyz 'ET' xvxz;
  25 : 'SINON';
  26 : xvxx = ('EXCO' var1 'Y1XX' 'SMXX') + ('EXCO' var1 'Y2XX' 'SMXX');
  27 : xvyy = ('EXCO' var1 'Y1YY' 'SMYY') + ('EXCO' var1 'Y2YY' 'SMYY');
  28 : xvzz = ('EXCO' var1 'Y1ZZ' 'SMZZ') + ('EXCO' var1 'Y2ZZ' 'SMZZ');
  29 : xvxy = ('EXCO' var1 'Y1XY' 'SMXY') + ('EXCO' var1 'Y2XY' 'SMXY');
  30 : xv = 'REDU' (xvxx 'ET' xvyy 'ET' xvzz 'ET' xvxy) mod0;
  31 : 'FINSI';
  32 : 'FINSI';
  33 : xv = 'CHAN' xv 'TYPE' 'CONTRAINTES';
  34 : j2 = 'VMIS' mod0 xv;
  35 : princxv = 'PRIN' mod0 xv;
  36 : xv1 = 'EXCO' princxv 'SI11' 'SCAL';
  37 : xv2 = 'EXCO' princxv 'SI22' 'SCAL';
  38 : xv3 = 'EXCO' princxv 'SI33' 'SCAL';
  39 : max1 = ( (xv1+xv2) + ('ABS' (xv1-xv2)) ) / 2.;
  40 : max2 = ( (max1+xv3) + ('ABS' (max1-xv3)) ) / 2.;
  41 : chi = (alpha*max2*1.5) + ((1.-alpha)*j2) ;
  42 : 'FINPROC' chi;
  43 : 
  44 :