Accueil | Cast3M

   1 : $$$$ RAFF     NOTICE  PASCAL    20/06/08    21:15:10     10623          
   2 :                                              DATE     20/06/08
   3 :  
   4 :    Operateur RAFF                           Voir aussi :
   5 :     
   6 :     1ere fonction : raffinement d'un maillage
   7 :     _______________                           
   8 : 
   9 :       GEO2  =  RAFF  GEO1  CHPO1 ; 
  10 : 
  11 : 
  12 :     Objet :
  13 :     _______
  14 : 
  15 :        L'operateur RAFF part d'un maillage existant (GEO1) pour le raffiner
  16 :     en respectant un champ de densite  (CHPO1). Tant que la densite n'est 
  17 :     pas atteinte un element est divise en sous elements etc... 
  18 :     Le maillage genere contient le resultat de la division des elements 
  19 :     plus des elements de types relations (itypel=22) qui permettent de 
  20 :     realiser au mieux la conformite en deplacement des elements. 
  21 : 
  22 : 
  23 :     Commentaire :
  24 :     _____________
  25 : 
  26 :     Entrees :
  27 : 
  28 :       - GEO1  :  Maillage initial
  29 :                                                    
  30 :       - CHPO1 :  Objet CHPOINT de densite                                         
  31 :     
  32 : 
  33 :     Sortie :
  34 : 
  35 :       - GEO2  : Maillage final, contenant les relations a imposer.                
  36 : 
  37 :     
  38 :     Remarque : 
  39 :     __________
  40 : 
  41 :     Apres avoir raffiné un maillage il est nécéssaire de faire appel à 
  42 :     l'opérateur RELA pour contruire les relations de conformité 
  43 :     entre les différentes zones de raffinement.
  44 : 
  45 : 
  46 :     2e fonction : raffinement d'une liste de reels
  47 :     -------------                           
  48 : 
  49 :       LRE2  =  RAFF  LRE1  | ENT1  ;
  50 :                            [ FLOT1
  51 : 
  52 :     Objet
  53 :     _____
  54 : 
  55 :        L'operateur RAFF permet egalement de raffiner un LISTREEL.
  56 : 
  57 : 
  58 :     Commentaire
  59 :     ___________
  60 : 
  61 :     Entrees :
  62 :       - LRE1  :  objet LISTREEL, liste initiale.
  63 :                                                      
  64 :       - ENT1  :  objet ENTIER, nombre de sous-decoupages de chaque
  65 :                  intervalle de la liste. ENT1 peut etre negatif
  66 :                  (voir remarque).
  67 : 
  68 :       - FLOT1 :  objet FLOTTANT, valeur cible de la taille
  69 :                  des intervalles de la liste.
  70 : 
  71 :     Sortie :
  72 :       - LRE2  :  objet LISTREEL, liste raffinee.
  73 : 
  74 : 
  75 :     Remarque 1
  76 :     __________
  77 : 
  78 :     Si ENT1 est negatif, la taille du decoupage de deux intervalles
  79 :     successifs de tailles differentes suit une progression geometrique.
  80 :     Le nombre d'intervalles peut alors etre different de celui attendu
  81 :     (abs(ENT1)).
  82 : 
  83 :     Remarque 2
  84 :     __________
  85 : 
  86 :     Si le LISREEL contient une succession de valeurs identiques,
  87 :     les intervalles entre ces veleurs sont egalement decoupes.
  88 :     Par exemple, la liste {1 1} raffinee d'un facteur 3 est
  89 :     la liste {1 1 1 1}.
  90 :     Ce comportement est souhaitable pour obtenir des LISTREEL
  91 :     de meme dimension en sortie de RAFF, quelles que soient
  92 :     les valeurs du LISTREEL en entree.
  93 : 
  94 : 
  95 : 
  96 :     Exemple d'utilisation de la 1ere fonction
  97 :     -----------------------------------------
  98 : 
  99 :     opti elem qua4 mode plan defo dime 2; 
 100 :     dens 2.;
 101 : * mesh 10x6 
 102 :     pa= 0 0; pb= 10 0;pc= 20 0;
 103 :     liab= pa droi pb;libc= pb droi pc;
 104 :     su = (liab et libc)  trans  ( 0 12);
 105 :     trac su;
 106 : * definition of density
 107 :     x y = coor su;
 108 :     distance = ((x - 10 ) * ( x- 10) + ( y * y)) ** 0.5;
 109 :     den = 0.3 + (0.18*distance);
 110 :     trac su den;
 111 : * new mesh 
 112 :     su2= raff su den;
 113 :     hh = elem su2 SURE ;
 114 :     sureal = su2 diff hh ;
 115 : * use of this mesh
 116 : * definition of model and caracteristic
 117 :     mo= mode su2 mecanique elastique isotrope ;    
 118 :     ma = mate mo  YOUN 2.e5  NU 0.3 ; 
 119 : * conformity relations
 120 :     rel = rela mo;
 121 : * loads
 122 :     psupe = su2 poin droite ( 0 12) ( 10 12) 0.1;
 123 :     lisupe = elem ( contour su2) appu stric psupe;
 124 :     ff = pres ( redu mo sureal) massif lisupe  -1.;
 125 : * displacements conditions
 126 :     py0= point su2 droit pa pc 0.01;
 127 :     liy0= elem ( contou su2) appu strict py0;
 128 :     li2bc = liy0 elem compris pb pc;
 129 :     cl1= bloqu li2bc UY;
 130 :     cl2= bloq UX pb;
 131 :     cltot= cl1 et cl2;
 132 : * compute elastic solution
 133 :     ri = rigi mo ma;
 134 :     displa =  reso ( ri et cltot et rel) ff;
 135 :     stre = sigma displa mo ma;
 136 :     vm = vmis stre mo ma;
 137 :     trac  su2 vm mo ma;    
 138 : * compute stress intensity factor
 139 :     gt = table;
 140 :     gt.'OBJECTIF' = MOT 'J';
 141 :     lifis = liy0 elem compris pa pb;   
 142 :     gt.'LEVRE_SUPERIEURE' = lifis;
 143 :     gt.'FRONT_FISSURE' = Pb;
 144 :     gt.'CARACTERISTIQUES' = ma;
 145 :     gt.'MODELE' = mo;
 146 :     gt.'SOLUTION_RESO' = displa; 
 147 :     rea = reaction  (ri et rel) displa;
 148 :     gt.'CHARGEMENTS_MECANIQUES'=rea;  
 149 :     naa = 5;opti veri 1;
 150 :     repe no naa;
 151 :       gt.'COUCHE' = &no;                
 152 :       G_THETA gt; 
 153 :       si ( &no ega 1) ; g2=prog gt.resultats;sinon;
 154 :       g2 = g2 et ( prog gt.resultats); finsi;
 155 :     fin no; 
 156 :     xx = prog 1 pas 1 naa;
 157 :     ev= evol manu  'nb of rows' xx 'G ' g2;
 158 :     ttt=table;
 159 :     ttt.1 = mot 'MARQ CROI';
 160 :     tt2=table;
 161 :     tt2.1= ' G '; 
 162 :     ttt.'TITRE'=tt2;
 163 :     dess ev lege ttt;
 164 : $$$$
 165 :  
 166 :  
 167 :  
 168 :