Procedures | Cast3M

   1 : * AUTOPILO  PROCEDUR  CHAT      11/03/16    21:15:35     6902           
   2 : 'DEBPROC' AUTOPILO DELT*'CHPOINT' DELA*'CHPOINT' ZMOD*'MMODEL'
   3 :  ZMAT*'MCHAML' TTT*TABLE;
   4 : *
   5 : *  cette procedure est utile en cas de limitation d'increment de
   6 : *  charge pour NONLIN option AUTOMATIQUE
   7 : *  En entree DELT est un champ de deplacements qui est l'increment de
   8 : *  deplacement que propose INCREME
   9 : *            DELA est la resolution elastqiue tangente du probleme
  10 : *  c'est a dire K-1.(F1-F0)
  11 : *            ZMOD et ZMAT contiennent le model et le champ par element
  12 : *  de materiau (et caracteristique) .
  13 : *            TTT est la table fournie en entree de NONLIN qui contient
  14 : *   notamment  TTT.'AUTOCRIT' qui est le critere de pilotage
  15 : *
  16 : *  Le but de cette procedure est de fournir un coefficient pour
  17 : *  normer DELT afin que ce dernier soit compatible avec le critere
  18 : *  de pilotage. Ce coefficient est une norme de DELT au sens du critere
  19 : *  En principe il est positif toutefois lui donner un signe permet
  20 : *  de forcer l'algorithme vers un direction donnée.
  21 : *
  22 : * EXEMPLE 1: Si on veut limiter a  TTT.'AUTOCRIT' l'increment de
  23 : *            deplacement pour un  point  P1  dans la direction
  24 : *            UY il faut faire :
  25 : *                     U1MA = ABS DELT 'EXTR' 'UR' A;
  26 : *
  27 : *
  28 : *  EXEMPLE 2 : on veut que le travail d'une force contenue dans
  29 : *              TTT.'TOTO' soit limite a TTT.'AUTOCRIT'
  30 : *              Il faut :
  31 : *                MODEP = MOTS 'UX' 'UY' 'UZ' 'RX' 'RY' 'RZ';
  32 : *                MOFOR = MOTS 'FX' 'FY' 'FZ' 'MX' 'MY' 'MZ';
  33 : *                U1MA = 'XTY' DELT  TTT.'TOTO'  MODEP MOFOR;
  34 : *
  35 : *
  36 : *  EXEMPLE 3 : On veut limiter a TTT.'AUTOCRIT' la norme de
  37 : *              l'increment de deplacement Max de epsilon
  38 : *               U1 = 'EPSI' ZMOD  DELT  ZMAT;
  39 : *               U1MA= 'MAXI' U1 'ABS';
  40 : *
  41 : *  en standard limitation sur le max de l'increment de deformation totale 
  42 : *  et de l'increment de deformation "plastique" (total - elastique) 
  43 : *
  44 : 
  45 :                U1 = 'EPSI' ZMOD  DELT  ZMAT;
  46 :                U2 = 'EPSI' ZMOD (DELT-DELA) ZMAT ;
  47 : 
  48 :             SI ( NON (TTT.PLASTIQUE OU TTT.VISCOPLASTIQUE OU
  49 :                      TTT.ENDOMMAGEMENT OU TTT.VISCODOMMAGE) ) ;
  50 : *            si on est en elastique pur on oublie les def. plas.
  51 :                   U2 = U1 ;
  52 :             FINSI ;
  53 :             U1MAX = 'MAXI' U1 'ABS';
  54 :             U2MAX = 'MAXI' U2 'ABS';
  55 :             NORM1 = MAXI ( PROG U1MAX U2MAX) ;
  56 :             
  57 : FINPROC   NORM1;
  58 : 
  59 : 
  60 : 
  61 : 
  62 : 
  63 :