Accueil | Cast3M

   1 : $$$$ RAFF     NOTICE  PASCAL    19/03/08    21:15:04     10153          
   2 :                                              DATE     19/03/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 : 
  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 coupes dans chaque intervalle
  65 :                  de valeurs.
  66 :                  ENT1 peut etre negatif (voir remarque ci-dessous).
  67 : 
  68 :     Sortie :
  69 :       - LRE2  :  objet LISTREEL, liste raffinee.
  70 : 
  71 : 
  72 :     Remarque 1
  73 :     __________
  74 : 
  75 :     Si ENT1 est negatif, la taille du decoupage de deux intervalles
  76 :     successifs de tailles differentes suit une progression geometrique.
  77 :     Le nombre d'intervalles peut alors etre different de celui attendu
  78 :     (abs(ENT1)).
  79 : 
  80 :     Remarque 2
  81 :     __________
  82 : 
  83 :     Si le LISREEL contient une succession de valeurs identiques,
  84 :     les intervalles entre ces veleurs sont egalement decoupes.
  85 :     Par exemple, la liste {1 1} raffinee d'un facteur 3 est
  86 :     la liste {1 1 1 1}.
  87 :     Ce comportement est souhaitable pour obtenir des LISTREEL
  88 :     de meme dimension en sortie de RAFF, quelles que soient
  89 :     les valeurs du LISTREEL en entree.
  90 : 
  91 : 
  92 : 
  93 :     Exemple d'utilisation de la 1ere fonction
  94 :     -----------------------------------------
  95 : 
  96 :     opti elem qua4 mode plan defo dime 2; 
  97 :     dens 2.;
  98 : * mesh 10x6 
  99 :     pa= 0 0; pb= 10 0;pc= 20 0;
 100 :     liab= pa droi pb;libc= pb droi pc;
 101 :     su = (liab et libc)  trans  ( 0 12);
 102 :     trac su;
 103 : * definition of density
 104 :     x y = coor su;
 105 :     distance = ((x - 10 ) * ( x- 10) + ( y * y)) ** 0.5;
 106 :     den = 0.3 + (0.18*distance);
 107 :     trac su den;
 108 : * new mesh 
 109 :     su2= raff su den;
 110 :     hh = elem su2 SURE ;
 111 :     sureal = su2 diff hh ;
 112 : * use of this mesh
 113 : * definition of model and caracteristic
 114 :     mo= mode su2 mecanique elastique isotrope ;    
 115 :     ma = mate mo  YOUN 2.e5  NU 0.3 ; 
 116 : * conformity relations
 117 :     rel = rela mo;
 118 : * loads
 119 :     psupe = su2 poin droite ( 0 12) ( 10 12) 0.1;
 120 :     lisupe = elem ( contour su2) appu stric psupe;
 121 :     ff = pres ( redu mo sureal) massif lisupe  -1.;
 122 : * displacements conditions
 123 :     py0= point su2 droit pa pc 0.01;
 124 :     liy0= elem ( contou su2) appu strict py0;
 125 :     li2bc = liy0 elem compris pb pc;
 126 :     cl1= bloqu li2bc UY;
 127 :     cl2= bloq UX pb;
 128 :     cltot= cl1 et cl2;
 129 : * compute elastic solution
 130 :     ri = rigi mo ma;
 131 :     displa =  reso ( ri et cltot et rel) ff;
 132 :     stre = sigma displa mo ma;
 133 :     vm = vmis stre mo ma;
 134 :     trac  su2 vm mo ma;    
 135 : * compute stress intensity factor
 136 :     gt = table;
 137 :     gt.'OBJECTIF' = MOT 'J';
 138 :     lifis = liy0 elem compris pa pb;   
 139 :     gt.'LEVRE_SUPERIEURE' = lifis;
 140 :     gt.'FRONT_FISSURE' = Pb;
 141 :     gt.'CARACTERISTIQUES' = ma;
 142 :     gt.'MODELE' = mo;
 143 :     gt.'SOLUTION_RESO' = displa; 
 144 :     rea = reaction  (ri et rel) displa;
 145 :     gt.'CHARGEMENTS_MECANIQUES'=rea;  
 146 :     naa = 5;opti veri 1;
 147 :     repe no naa;
 148 :       gt.'COUCHE' = &no;                
 149 :       G_THETA gt; 
 150 :       si ( &no ega 1) ; g2=prog gt.resultats;sinon;
 151 :       g2 = g2 et ( prog gt.resultats); finsi;
 152 :     fin no; 
 153 :     xx = prog 1 pas 1 naa;
 154 :     ev= evol manu  'nb of rows' xx 'G ' g2;
 155 :     ttt=table;
 156 :     ttt.1 = mot 'MARQ CROI';
 157 :     tt2=table;
 158 :     tt2.1= ' G '; 
 159 :     ttt.'TITRE'=tt2;
 160 :     dess ev lege ttt;
 161 : $$$$
 162 :  
 163 :  
 164 :