Accueil | Cast3M

   1 : $$$$ CH_THETA NOTICE  JB251061  21/02/05    21:15:01     10879          
   2 :                                              DATE     21/02/05
   3 : 
   4 :   Procedure CH_THETA                       Voir aussi : G_THETA
   5 :                                                           G_CALCUL  
   6 :                                                           CH_THETX  
   7 :                                                           G_CAS  
   8 :                                                           G_AUX  
   9 :     ------------------  
  10 :     CH_THETA SUPTAB OBJUTI BOOL ;
  11 : 
  12 :                               SUPTAB.'MAILLAGE'
  13 :                                     .'FISSURE'
  14 :                                     .'FRONT_FISSURE'
  15 :                                     .'COUCHE'
  16 :                                     .'POINT_1'
  17 :                                     .'POINT_2'
  18 :                                     .'POINT_3'
  19 :                                     .'PCENTRE'
  20 :                                     .'CHPOINT_TRANSFORMATION'
  21 :                                     .'OPERATEUR'
  22 :                                     .'EPAISSEUR'
  23 : 
  24 : 
  25 :     Objet :
  26 :     _______
  27 : 
  28 :   Cette procedure determine un chpoint de type THETA. c'est-a-dire  
  29 :   un champ/point dont la norme est constante a l'interieur d'une      
  30 :   couronne entourant le front d'une fissure, zero a l'exterieur   
  31 :   cette couronne. Le vecteur represente par ce chpoint indique  
  32 :   la direction de propagation eventuelle de la fissure. Pour un   
  33 :   front de fissure tridimensionnel, la procedure cree autant de
  34 :   champs THETA qu'il y a de points sur le front de fissure. Les
  35 :   champs THETA construits sont des champs locaux qui s'appuient
  36 :   sur les noeuds du plan normal au front de fissure au point
  37 :   considere (le maillage doit etre elabore de maniere a prevoir
  38 :   l'existence de ces plans normaux). En 3D Cette procedure est 
  39 :   applicable uniquement pour les fissures planes. 
  40 : 
  41 :     En entree
  42 :     ---------
  43 : 
  44 :     SUPTAB = Objet de type TABLE dont les indices sont des              
  45 :              objets de type MOT (a ecrire en toutes lettres).
  46 :              
  47 :      
  48 :     OBJUTI = Objet de type TABLE dont les indices sont des              
  49 :              objets de type MOT (a ecrire en toutes lettres).
  50 :              Cette table est definie par les procedures G_THETA
  51 :              et G_CAS avant de faire appel a CH_THETA.
  52 :              
  53 :              
  54 :     BOOL = Objet de type TABLE dont les indices sont des              
  55 :            objets de type MOT (a ecrire en toutes lettres).
  56 :            Cette table est definie par les procedures G_THETA
  57 :            et G_CAS avant de faire appel a CH_THETA.
  58 : 
  59 : 
  60 :     1) ARGUMENTS OBLIGATOIRES DANS TOUS LES CAS
  61 :     -------------------------------------------
  62 :     SUPTAB.'MAILLAGE' = Objet de type MAILLAGE representant soit        
  63 :                         la structure totale etudiee (maillage           
  64 :                         utilise dans l'analyse par elements finis,      
  65 :                         soit, pour reduire le temps de calcul, le       
  66 :                         maillage entourant le plus grand des contours   
  67 :                         qu'on a defini pour calculer le champ THETA.    
  68 :     SUPTAB.'FISSURE'  = Objet de type MAILLAGE donnant les deux  
  69 :                         levres d'une fissure si elle est complete (la   
  70 :                         fissure presente des noeuds doubles), une       
  71 :                         seule levre si l'autre levre n'est pas maillee  
  72 :                         en raison, par example, de la symetrie du       
  73 :                         probleme.                                       
  74 :     SUPTAB.'FRONT_FISSURE' = Objet de type POINT (representant la       
  75 :                              pointe de la fissure) si la fissure est    
  76 :                              une ligne, de type MAILLAGE (representant  
  77 :                              le front de la fissure) si la fissure est  
  78 :                              sur un plan en 3D.                         
  79 :     SUPTAB.'COUCHE'   = Objet de type ENTIER representant le nombre     
  80 :                         de couches d'elements (autour du point de       
  81 :                         fissure) qui se deplacent pour simuler la       
  82 :                         propagation de la fissure.                       
  83 : 
  84 : 
  85 :     2) CAS D'UNE FISSURE CIRCULAIRE DANS UNE GEOMETRIE PLANE               
  86 :     --------------------------------------------------------
  87 :                                                                         
  88 :     SUPTAB.'PCENTRE'  = centre de la fissure circulaire                 
  89 :                                                                         
  90 :                                                                         
  91 :     3) CAS OU L EXTENSION DE FISSURE CORRESPOND A UNE                      
  92 :     SIMPLE TRANSLATION DANS UN TUYAUTERIE DROITE (3D)                   
  93 :     -------------------------------------------------
  94 :                                                                         
  95 :     Dans ce cas on effectue dans la procedure une transformation de     
  96 :     tuyau en plaque en passant au systeme de coordonnees cylindriques.  
  97 :     Il est alors necessaire de fournir :                                
  98 :                                                                         
  99 :     SUPTAB.'POINT_1' = centre du systeme de coordonnees                 
 100 :     SUPTAB.'POINT_2' = POINT tel que l'axe defini par POINT_1           
 101 :                        vers POINT_2 soit l'axe Z poisitif               
 102 :     SUPTAB.'POINT_3' = POINT tel que le plan defini par les 3 points    
 103 :                        POINT_1 POINT_2 POINT_3 donne l'angle theta nul  
 104 :                                                                         
 105 :     4) CAS OU L EXTENSION DE FISSURE NE
 106 :     CORRESPOND PAS A UNE SIMPLE TRANSLATION                                            
 107 :     ---------------------------------------
 108 :                                                                         
 109 :     A) Fissure dans un tuyauterie droite (3D, Rotation)                 
 110 :                                                                         
 111 :     SUPTAB.'POINT_1' = Objet de type POINT                              
 112 :     SUPTAB.'POINT_2' = Objet de type POINT qui, avec le point POINT_1,  
 113 :                        constitue l'axe perpendiculaire a la section     
 114 :                        fissuree.                                        
 115 :                                                                         
 116 :     B) Fissure dans un coude (3D, rotation + transformation)            
 117 :                                                                         
 118 :        Outre les deux points SUPTAB.'POINT_1' et SUPTAB.'POINT_2'       
 119 :        definis en haut on donne encore :                                
 120 :     SUPTAB.'CHPOINT_TRANSFORMATION' = Objet de type CHPOINT utilise     
 121 :                                       pour transformer une coude en     
 122 :                                       un tuyauterie droite.             
 123 :                                                                         
 124 :     SUPTAB.'OPERATEUR' = Objet de type MOT valant 'PLUS' ou 'MOIN'      
 125 :                          pour indiquer l'operateur PLUS ou MOIN a       
 126 :                          utiliser si l'on veut transformer la coude     
 127 :                          en un tuyauterie droite.                       
 128 :                                                                         
 129 :     5) CAS DES ELEMENTS DE COQUE                                        
 130 :     ----------------------------
 131 :                                                                         
 132 :     SUPTAB.'EPAISSEUR' = Objet de type FLOTTANT donnant l'epaisseur     
 133 :                          de la coque a la pointe de la fissure     
 134 :                          
 135 : 
 136 : 
 137 : 
 138 :       En sortie :
 139 :       -----------
 140 : 
 141 :     SUPTAB.'CHAMP_THETA' = Objet de type :                                              
 142 :                                                                        
 143 :         - TABLE indicee par des objets de type POINT contenant
 144 :           des elements de type CHPOINT dans les 3 dimensions.
 145 :           Chaque element contient le champ THETA au noeud du 
 146 :           front, de coordonnees correspondant au point P: TETA.P.
 147 :           Elle est egalement indicee par le mot 'GLOBAL' pour 
 148 :           donner le champ THETA global le long de tout le front
 149 :           de la fissure
 150 :         - CHPOINT contenant le champ THETA en 2 dimensions (ou en
 151 :           3 dimensions avec des elements coque mince) a la pointe 
 152 :           de fissure. 
 153 :           
 154 :     SUPTAB.'CHAMP_PI' = Objet de meme type que SUPTAB.'CHAMP_THETA'.
 155 :                         Calcule dans le cas ou l'objectif est 'DJ/DA'.
 156 :            
 157 :     OBJUTI.'DIRECTION1' (et 2 et 3 eventuellement) :
 158 :         Objet de type POINT donnant la direction de propagation,
 159 :          ainsi que le ou les vecteurs normaux le cas echeant.
 160 : 
 161 : 
 162 :