Accueil | Cast3M

   1 : $$$$ @ANA_LIM NOTICE  CHAT      11/09/12    21:15:06     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :    Procedure @ANA_LIM                        Voir aussi :                     
   5 :     ------------------  
   6 :     @ANA_LIM TAB1;
   7 : 
   8 :    TAB1.  BLOCAGES_MECANIQUES           MAXITERATION
   9 :           BRIDE                         MECARUINE
  10 :           CADIM                         MODELE
  11 :           CALCUL_COQUE                  PAS 
  12 :           CHARGEMENT                    PRECISION
  13 :           CONTRAINTES                   PROCEDURE_PERSO1
  14 :           CRITERE                       REPRISE
  15 :           DEPLA                         RESULTAT
  16 :           EVOLCL                        SIGLIM
  17 :           LISTCL                        
  18 :           LISTITER 
  19 :          
  20 : 
  21 :     Objet :
  22 :     _______
  23 : 
  24 : Cette procedure permet de calculer le chargement limite d'une structure
  25 : sous la forme d'une suite monotone decroissante convergeant vers la 
  26 : solution. 
  27 : La modelisation est equivalente a un calcul plastique parfait.
  28 : Version 04/96 
  29 : Elle s'utilise comme PASAPAS
  30 : 
  31 : Commentaires
  32 : ------------
  33 : 
  34 : TAB1   : Objet de type TABLE.                                        
  35 : 
  36 : En entree, TAB1 sert a definir les options et les  parametres du calcul.
  37 : Les indices de l'objet TAB1 sont des mots dont voici la liste :                 
  38 : *********************************************************************           
  39 :         
  40 : MODELE          MMODEL          objet modele (MECANIQUE ELASTIQUE)
  41 :                                 Indispensable
  42 : 
  43 : CHARGEMENT      CHPOINT         Forces nodales equivalentes 
  44 :                                 Indispensable
  45 :         
  46 : BLOCAGES_MECANIQUES RIGIDITE    blocages mecaniques
  47 :                                 Indispensable
  48 :                                         
  49 : CRITERE         MOT             choix du critere de plasticite
  50 :                                 VONMISES OU tresca (taper TRSCA)
  51 :                                 defaut: VONMISES
  52 :                                 
  53 : SIGLIM          FLOTTANT        limite d'ecoulement 
  54 :                 ou MCHAML       pour un MCHAML il faut un type SCALAIRE
  55 :                                 appuye aux points de calcul des 
  56 :                                 contraintes
  57 :                                 (voir operateur CHAN -> type "SCALAIRE"
  58 :                                 et "STRESSES")
  59 :                                  valeur par defaut: 1.
  60 : 
  61 : CALCUL_COQUE    MCHAML          active un calcul en coque
  62 :         stype CARACTERISTIQUES  Caracteristiques de la coque (EPAI,
  63 :                                 ALFA) pour un veritable calcul 
  64 :                                  plastique ALFA=2/3
  65 : 
  66 : PRECISION       FLOTTANT        valeur de la precision           
  67 :                                         valeur par defaut: 1e-3
  68 :         
  69 : MAXITERATION     ENTIER         nombre maximum d'iterations      
  70 :                                         valeur par defaut: 50
  71 :         
  72 : PAS              ENTIER         pas d'affichage des resultats      
  73 :                                         valeur par defaut: 1
  74 :        
  75 : 
  76 : REPRISE         LOGIQUE         si REPRISE=VRAI il s'agit de la
  77 :                                 reprise d'un calcul
  78 :                                 (exemple a la fin)
  79 :                                 
  80 : PROCEDURE_PERSO1 LOGIQUE        VRAI si procedure perso a la fin de 
  81 :                                 chaque iteration
  82 :                                 
  83 :         
  84 : En sortie, TAB1 permet de retrouver les resultats:
  85 : 
  86 :  indice         type objet
  87 : 
  88 : MECARUINE       MCHAML          VMises a convergence            
  89 :           sstype CONTRAINTES    
  90 :         
  91 : CONTRAINTES     MCHAML          contraintes a convergence
  92 : 
  93 : 
  94 : YOUNG           MCHAML          young a convergence
  95 :                 
  96 : 
  97 : LISTITER        LISTREEL        liste des iterations
  98 :                                         
  99 :                                         
 100 : LISTCL          LISTREEL        liste des charges limites
 101 :                                         
 102 :                                         
 103 : EVOLCL          EVOLUTION       evolution des charges limites
 104 :                                         en fonction des iterations
 105 :                                         
 106 : CADIM           FLOTTANT        charge limite resultat de 
 107 :                                 chargement limite/chargement elastique
 108 :                                     
 109 : DEPLA           CHPO            deplacements a convergence
 110 :         
 111 :        
 112 : 
 113 : le barometre permet de verifier que les Von Mises (ou les Tresca )
 114 : s'appuient bien sur la surface de charge.
 115 : 
 116 : liste des avertissements:
 117 : 
 118 : AVERT1:vous etes totalement elastique (passage au dessous de la limite
 119 : elastique)
 120 : AVERT2:la suite des charges limites est croissante -> invalidite de la 
 121 : theorie
 122 : AVERT3:pas assez d'iterations
 123 : 
 124 : remarques:
 125 : -@ANA_LIM fonctionne soit en depplacements imposes soit en forces
 126 : imposes,mais pas avec les deux conditions pour un meme calcul.
 127 : Lorsqu'un calcul est fait en 'forces imposees', le parametre "charge 
 128 : limite" represente la valeur a multiplier au chargement initial imposes
 129 : pour obtenir le chargement limite.
 130 : - Lorsqu'un calcul est fait en 'deplacements imposes', il convient de
 131 : depouiller le chargement a la fin de chaque iteration par l'intermediaire
 132 : d'une procedure PERSO1 qui traitera le resultat pour le transformer en 
 133 : reactions aux endroits d'application des deplacements (REDU BSIG)
 134 : - Pour une reprise, a la fin du premier calcul, ajouter 
 135 : table.REPRISE=vrai et relancer @ANA_LIM 
 136 : - Le "ratio elastique" correspond au Vmises maxi apres un calcul 
 137 : elastique incompressible  divise par la limite d'ecoulement.
 138 : Pour plus de precisions voir David Plancq (42257492)