Accueil | Cast3M

   1 : $$$$ RAFF     NOTICE  GG250959  18/01/09    21:15:25     9688           
   2 :                                              DATE     18/01/09
   3 :  
   4 :    Operateur RAFF                           Voir aussi :
   5 :     --------------  
   6 :       GEO2  =  RAFF  GEO1  CHPO1 ; 
   7 :                                 
   8 : 
   9 : 
  10 : 
  11 :     Objet :
  12 :     _______
  13 : 
  14 :     L'operateur RAFF part d'un maillage existant (GEO1) pour le raffiner
  15 :  en respectant un champ de densite  (CHPO1). Tant que la densite n'est 
  16 :  pas atteinte un element est divise en sous elements etc... 
  17 :  Le maillage genere contient le resultat de la division des elements 
  18 :  plus des elements de types relations (itypel=22) qui permettent de 
  19 :  realiser au mieux la conformite en deplacement des elements. 
  20 : 
  21 : 
  22 :       Commentaire :
  23 :       _____________
  24 : 
  25 :       GEO1  :  Maillage initial
  26 :                                                      
  27 :       CHPO1 :  Objet CHPOINT de densite                                         
  28 :                                                      
  29 :       
  30 :       GEO2  : Maillage final, contenant les relations a imposer.                
  31 : 
  32 :       
  33 :       Remarque : 
  34 :       __________
  35 : 
  36 :       Apres avoir raffinÃ© un maillage il est nÃ©cÃ©ssaire de faire appel Ã  
  37 :       l'opÃ©rateur RELA pour contruire les relations de conformitÃ© 
  38 :       entre les diffÃ©rentes zones de raffinement.
  39 : 
  40 :    Exemple :
  41 :    _________
  42 : 
  43 :     opti elem qua4 mode plan defo dime 2; 
  44 :     dens 2.;
  45 : * mesh 10x6 
  46 :     pa= 0 0; pb= 10 0;pc= 20 0;
  47 :     liab= pa droi pb;libc= pb droi pc;
  48 :     su = (liab et libc)  trans  ( 0 12);
  49 :     trac su;
  50 : * definition of density
  51 :     x y = coor su;
  52 :     distance = ((x - 10 ) * ( x- 10) + ( y * y)) ** 0.5;
  53 :     den = 0.3 + (0.18*distance);
  54 :     trac su den;
  55 : * new mesh 
  56 :     su2= raff su den;
  57 :     hh = elem su2 SURE ;
  58 :     sureal = su2 diff hh ;
  59 : * use of this mesh
  60 : * definition of model and caracteristic
  61 :     mo= mode su2 mecanique elastique isotrope ;    
  62 :     ma = mate mo  YOUN 2.e5  NU 0.3 ; 
  63 : * conformity relations
  64 :     rel = rela mo;
  65 : * loads
  66 :     psupe = su2 poin droite ( 0 12) ( 10 12) 0.1;
  67 :     lisupe = elem ( contour su2) appu stric psupe;
  68 :     ff = pres ( redu mo sureal) massif lisupe  -1.;
  69 : * displacements conditions
  70 :     py0= point su2 droit pa pc 0.01;
  71 :     liy0= elem ( contou su2) appu strict py0;
  72 :     li2bc = liy0 elem compris pb pc;
  73 :     cl1= bloqu li2bc UY;
  74 :     cl2= bloq UX pb;
  75 :     cltot= cl1 et cl2;
  76 : * compute elastic solution
  77 :     ri = rigi mo ma;
  78 :     displa =  reso ( ri et cltot et rel) ff;
  79 :     stre = sigma displa mo ma;
  80 :     vm = vmis stre mo ma;
  81 :     trac  su2 vm mo ma;    
  82 : * compute stress intensity factor
  83 :     gt = table;
  84 :     gt.'OBJECTIF' = MOT 'J';
  85 :     lifis = liy0 elem compris pa pb;   
  86 :     gt.'LEVRE_SUPERIEURE' = lifis;
  87 :     gt.'FRONT_FISSURE' = Pb;
  88 :     gt.'CARACTERISTIQUES' = ma;
  89 :     gt.'MODELE' = mo;
  90 :     gt.'SOLUTION_RESO' = displa; 
  91 :     rea = reaction  (ri et rel) displa;
  92 :     gt.'CHARGEMENTS_MECANIQUES'=rea;  
  93 :     naa = 5;opti veri 1;
  94 :     repe no naa;
  95 :       gt.'COUCHE' = &no;                
  96 :       G_THETA gt; 
  97 :       si ( &no ega 1) ; g2=prog gt.resultats;sinon;
  98 :       g2 = g2 et ( prog gt.resultats); finsi;
  99 :     fin no; 
 100 :     xx = prog 1 pas 1 naa;
 101 :     ev= evol manu  'nb of rows' xx 'G ' g2;
 102 :     ttt=table;
 103 :     ttt.1 = mot 'MARQ CROI';
 104 :     tt2=table;
 105 :     tt2.1= ' G '; 
 106 :     ttt.'TITRE'=tt2;
 107 :     dess ev lege ttt;
 108 : $$$$
 109 :  
 110 :  
 111 :