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 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales