Accueil | Cast3M

   1 : $$$$ RAFF     NOTICE  CHAT      11/09/12    21:17:49     7124           
   2 :                                              DATE     11/09/12
   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 : 
  34 :    Exemple :
  35 :    _________
  36 : 
  37 :     opti elem qua4 mode plan defo dime 2; 
  38 :     dens 2.;
  39 : * mesh 10x6 
  40 :     pa= 0 0; pb= 10 0;pc= 20 0;
  41 :     liab= pa droi pb;libc= pb droi pc;
  42 :     su = (liab et libc)  trans  ( 0 12);
  43 :     trac su;
  44 : * definition of density
  45 :     x y = coor su;
  46 :     distance = ((x - 10 ) * ( x- 10) + ( y * y)) ** 0.5;
  47 :     den = 0.3 + (0.18*distance);
  48 :     trac su den;
  49 : * new mesh 
  50 :     su2= raff su den;
  51 :     hh = elem su2 SURE ;
  52 :     sureal = su2 diff hh ;
  53 : * use of this mesh
  54 : * definition of model and caracteristic
  55 :     mo= model su2 mecanique elastique isotrope ;    
  56 :     ma = matr mo  YOUN 2.e5  NU 0.3 ; 
  57 : * loads
  58 :     psupe = su2 point droite ( 0 12) ( 10 12) 0.1;
  59 :     lisupe = elem ( contour su2) appu stric psupe;
  60 :     ff = pres ( redu mo sureal) massif lisupe  -1.;
  61 : * displacements conditions
  62 :     py0= point su2 droit pa pc 0.01;
  63 :     liy0= elem ( contou su2) appu strict py0;
  64 :     li2bc = liy0 elem compris pb pc;
  65 :     cl1= bloqu li2bc UY;
  66 :     cl2= bloq UX pb;
  67 :     cltot= cl1 et cl2;
  68 : * compute elastic solution
  69 :     ri = rigi mo ma;
  70 :     displa =  reso ( ri et cltot) ff;
  71 :     stre = sigma displa mo ma;
  72 :     vm = vmis stre mo ma;
  73 :     trac  su2 vm mo ma;    
  74 : * compute stress intensity factor
  75 :     gt = table;
  76 :     gt.'OBJECTIF' = MOT 'J';
  77 :     lifis = liy0 elem compris pa pb;   
  78 :     gt.'LEVRE_SUPERIEURE' = lifis;
  79 :     gt.'FRONT_FISSURE' = Pb;
  80 :     gt.'CARACTERISTIQUES' = ma;
  81 :     gt.'MODELE' = mo;
  82 :     gt.'SOLUTION_RESO' = displa; 
  83 :     rea = reaction  ri displa;
  84 :     gt.'CHARGEMENTS_MECANIQUES'=rea;  
  85 :     naa = 5;opti veri 1;
  86 :     repe no naa;
  87 :       gt.'COUCHE' = &no;                
  88 :       G_THETA gt; 
  89 :       si ( &no ega 1) ; g2=prog gt.resultats;sinon;
  90 :       g2 = g2 et ( prog gt.resultats); finsi;
  91 :     fin no; 
  92 :     xx = prog 1 pas 1 naa;
  93 :     ev= evol manu  'nb of rows' xx 'G ' g2;
  94 :     ttt=table;
  95 :     ttt.1 = mot 'MARQ CROI';
  96 :     tt2=table;
  97 :     tt2.1= ' G '; 
  98 :     ttt.'TITRE'=tt2;
  99 :     dess ev lege ttt;
 100 : $$$$
 101 :  
 102 :