Télécharger la plaquette Cast3M


fr en

Télécharger nafems-le3-ortho.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  nafems-le3-ortho.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *** Options ...
  6.  
  7.     opti dime 3 mode trid elem seg2 echo 1 ;
  8.     graph = 'N' ;
  9.  
  10. *** Paramètres ...
  11.  
  12.     rayon = 10. ;
  13.  
  14.     nb_div = 2 ;
  15.  
  16.     valyoun = 6.825e+7 ;
  17.     valnu   = 0.3 ;
  18.     valcis  = valyoun / (2 * (1.0+valnu)) ;
  19.     lepaiss = 0.04 ;
  20.  
  21. *** Points ...
  22.  
  23.     dens (pi * rayon / 4 / nb_div) ;
  24.  
  25.     x1 = rayon ;
  26.     x2 = (rayon * rayon / 2) ** 0.5 ;
  27.     x3 = (rayon * rayon / 3) ** 0.5 ;
  28.  
  29.     p0 = 0. 0. 0. ;
  30.     px = x1 0. 0. ;
  31.     py = 0. x1 0. ;
  32.     pz = 0. 0. x1 ;
  33.  
  34.     pxy = x2 x2 0. ;
  35.     pxz = x2 0. x2 ;
  36.     pyz = 0. x2 x2 ;
  37.  
  38.     pxyz = x3 x3 x3 ;
  39.  
  40. *** Lignes ...
  41.  
  42.     lixy = cerc nb_div px p0 pxy ;
  43.     lixz = cerc nb_div px p0 pxz ;
  44.  
  45.     liyx = cerc nb_div py p0 pxy ;
  46.     liyz = cerc nb_div py p0 pyz ;
  47.  
  48.     lizx = cerc nb_div pz p0 pxz ;
  49.     lizy = cerc nb_div pz p0 pyz ;
  50.  
  51.     licxy = cerc nb_div pxyz p0 pxy ;
  52.     licxz = cerc nb_div pxyz p0 pxz ;
  53.     licyz = cerc nb_div pxyz p0 pyz ;
  54.  
  55. *** Surfaces ...
  56.  
  57.     opti elem qua4 ;
  58.  
  59.     quartx = dall lixy (inve licxy) licxz (inve lixz) sphe p0 ;
  60.     quarty = dall liyz (inve licyz) licxy (inve liyx) sphe p0 ;
  61.     quartz = dall lizx (inve licxz) licyz (inve lizy) sphe p0 ;
  62.  
  63.     su1 = quartx et quarty et quartz ;
  64.  
  65.     si ( EGA   graph  'O') ;
  66.        titr 'Le maillage (un quart de la hemisphere)' ;
  67.        trac su1 ;
  68.     finsi ;
  69.  
  70. *** Modèle ...
  71.  
  72.     mo1 = mode su1 mecanique elastique orthotrope coq4 ;
  73.     ma1 = mate mo1 DIRECTION (1 -1 1) PARALLELE
  74.               YG1 valyoun YG2 valyoun NU12 valnu 
  75.               G12 valcis G23 valcis G13 valcis
  76.               EPAI lepaiss ;
  77.  
  78. *** Rigidité et CL ...
  79.  
  80.     ri1 = rigi mo1 ma1 ;
  81.     cl1 = symt depl rota px p0 pz su1 ;
  82.     cl2 = symt depl rota py p0 pz su1 ;
  83.     cl3 = bloq uz pz ;
  84.  
  85.     rielas = ri1 et cl1 et cl2 et cl3 ;
  86.  
  87.     si(ega 'O' graph) ;
  88.        titr 'La structure avec ses conditions aux limites' ;
  89.        trac rielas ;
  90.     finsi ;
  91.  
  92. *** Chargement ...
  93.  
  94.     val_forc = 1. ;
  95.  
  96.     fo1 = forc fx     val_forc  px ;
  97.     fo2 = forc fy (-1*val_forc) py ;
  98.  
  99.     fortot = fo1 et fo2 ;
  100.  
  101. *** Résolution ...
  102.  
  103.     de1 = reso rielas fortot ;
  104.  
  105.     si(ega 'O' graph) ;
  106.        titr 'La structure deformee' ;
  107.        vf1 = vect fortot FX FY FZ (rayon / val_forc / 5.) JAUN ;
  108.        defo0 = defo su1 de1 0.     VERT ;
  109.        defo1 = defo su1 de1    vf1 ROUG ;
  110.        trac (defo0 et defo1) ;
  111.     finsi ;
  112.  
  113. *** Test ...
  114.  
  115.     val_ref = 0.0925 ;
  116.  
  117.     dx =       extr de1 ux px  ;
  118.     dy = -1 * (extr de1 uy py) ;
  119.  
  120.     errx = 100 * (dx - val_ref) / val_ref ;
  121.     erry = 100 * (dy - val_ref) / val_ref ;
  122.  
  123.     opti echo 0 ;
  124.     mess 'Déplacement de référence = ' val_ref ;
  125.     mess 'Déplacement en X =' dx ' soit' errx ' % d erreur' ;
  126.     mess 'Déplacement en Y =' dy ' soit' erry ' % d erreur' ;
  127.     opti echo 1 ;
  128.  
  129. *** Bye ...
  130.  
  131.     fin ;
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales