Télécharger la plaquette Cast3M


fr en

Télécharger @initia.procedur

Retour à la liste

Numérotation des lignes : 

  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.  
  80.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales