Procedures | Cast3M

   1 : * @INITIA   PROCEDUR  KICH      99/05/31    21:15:08     3569           
   2 : DEBP @INITIA ttt*TABLE ;
   3 : 
   4 : *-----------------------------------------------------------------------
   5 : * @INITIA : initialise le tableau de calcul en
   6 : *          - stockant les conditions aux limites (blocages + rigidite)
   7 : *            et les contraintes initiales (mises a 0)
   8 : *          - initialisant a 0 les deformations plastiques, variables
   9 : *            internes ...
  10 : *          - resolvant le probleme elastique a l'iteration 0
  11 : *-----------------------------------------------------------------------
  12 : 
  13 :   ttt . 'DEPLACEMENTS'            = tabl ;
  14 :   ttt . 'DEFORMATIONS'            = tabl ;
  15 :   ttt . 'CONTRAINTES'             = tabl ;
  16 :   ttt . 'DEFORMATIONS_PLASTIQUES' = tabl ;
  17 :   ttt . 'VARIABLES_INTERNES'      = tabl ;
  18 :   ttt . 'CONTRAINTES_PLASTIQUES'  = tabl ;
  19 :   ttt . 'FORCES_PLASTIQUES'       = tabl ;
  20 :   ttt . 'CRITERE_PLASTICITE'      = tabl ;
  21 :   ttt . 'CONVERGENCE'             = tabl ;
  22 :   si ((vale DIME) ega 2) ;
  23 :     ttt . EP2D = tabl ;
  24 :     ttt . VA2D = tabl ;
  25 :   finsi ;
  26 :   si ((vale DIME) ega 3) ;
  27 :     ttt . EP3D = tabl ;
  28 :     ttt . VA3D = tabl ;
  29 :   finsi ;
  30 : 
  31 : *** stockage des donnees
  32 : 
  33 :   ttt . 'BLOCAGES_MECANIQUES'   = (ttt . 'BLOCAGES_MECANIQUES') et 
  34 :                                   (rigi (ttt . 'MODELE')
  35 :                                   (ttt . 'CARACTERISTIQUES')) ;
  36 :   ttt . 'CONTRAINTES_INITIALES' = zero (ttt . 'MODELE') CONTRAIN ;
  37 : 
  38 : *** initialisation des tableaux de calcul
  39 : 
  40 :   ttt . 'CONVERGENCE'             . 0 = FAUX ;
  41 : *  ttt . 'DEFORMATIONS_PLASTIQUES' . 0 = zero (ttt . 'MODELE') DEFORMAT ;
  42 :   ttt . 'DEFORMATIONS_PLASTIQUES' . 0 = zero (ttt . 'MODELE') DEFINELA ;
  43 :   ttt . 'VARIABLES_INTERNES'      . 0 = zero (ttt . 'MODELE') VARINTER ;
  44 :  lm2 = extr ttt . modele deformat ;
  45 :  lm1 = extr ttt . 'DEFORMATIONS_PLASTIQUES' . 0 comp ;
  46 :  def_ut = changer comp ttt . 'DEFORMATIONS_PLASTIQUES' . 0 lm1 lm2 ;
  47 :  def_ut = changer def_ut type 'DEFORMATIONS' ;
  48 : *  ttt . 'CONTRAINTES_PLASTIQUES'  . 0 = (ttt . 'MODELE')
  49 : *                   (hook (ttt . 'MODELE') (ttt . 'CARACTERISTIQUES')) * 
  50 : *                   (ttt . 'DEFORMATIONS_PLASTIQUES' . 0) ;
  51 :   ttt . 'CONTRAINTES_PLASTIQUES'  . 0 = (ttt . 'MODELE')
  52 :                    (hook (ttt . 'MODELE') (ttt . 'CARACTERISTIQUES')) * 
  53 :                   def_ut  ;
  54 :   ttt . 'FORCES_PLASTIQUES'       . 0 = bsig (ttt . 'MODELE')
  55 :                    (ttt . 'CONTRAINTES_PLASTIQUES' . 0) ;
  56 : 
  57 :   si ((vale DIME) ega 2) ;
  58 : *    ttt . EP2D . 0 = zero (MOD . 1) DEFORMAT ;
  59 :     ttt . EP2D . 0 = zero (MOD . 1) DEFINELA ;
  60 :     ttt . VA2D . 0 = zero (MOD . 1) VARINTER ;
  61 :   finsi ;
  62 :   si ((vale DIME) ega 3) ;
  63 : *    ttt . EP3D . 0 = zero (MOD . 1) DEFORMAT ;
  64 :     ttt . EP3D . 0 = zero (MOD . 1) DEFINELA ;
  65 :     ttt . VA3D . 0 = zero (MOD . 1) VARINTER ;
  66 :   finsi ;
  67 : 
  68 : *** solution elastique a l'iteration 0
  69 : 
  70 :   ttt . 'DEPLACEMENTS' . 0 = reso (ttt . 'BLOCAGES_MECANIQUES')
  71 :                                   (ttt . 'CHARGEMENT') ;
  72 :   ttt . 'DEFORMATIONS' . 0 = epsi (ttt . 'MODELE')
  73 :                                   (ttt . 'DEPLACEMENTS' . 0) ;
  74 :   ttt . 'CONTRAINTES'  . 0 = sigm (ttt . 'MODELE')
  75 :                                   (ttt . 'CARACTERISTIQUES')
  76 :                                   (ttt . 'DEPLACEMENTS' . 0) ;
  77 : 
  78 : FINP ttt ;
  79 :