Accueil | Cast3M

   1 : $$$$ MEC3     NOTICE  CHAT      11/09/12    21:17:13     7124           
   2 :                                              DATE     11/09/12
   3 :    Voir aussi :
   4 : *--------------------------------------------------------------------*
   5 : *                                                                    *
   6 : *           Exemple simple de calcul mecanique elastique             *
   7 : *    fonctionnant avec les nouvelles structures MMODEL et MCHAML     *
   8 : *                                                                    *
   9 : *--------------------------------------------------------------------*
  10 : * Definition des options
  11 : *
  12 :            OPTI DIME 3 ELEM CU20 MODEL TRIDIM ; DENS 10 ;
  13 : *
  14 : * Maillage d'une barre de 5 elements CU20
  15 : *
  16 :   P1 = 0 0 0 ; P2 = 10 0 0 ; LI= P1 D 1 P2 ; SUB=LI TRAN (0 10 0) ;
  17 :    TOTAL = SUB VOLU TRAN (0 0 50) COUL JAUN ; SUH = TOTAL FACE 2 ;
  18 :      LIH = SUH COTE 2 ; TRAC TOTAL CACH (1000 -2000 1000) QUAL ;
  19 : *
  20 : * Formulation, materiau defini a l'aide de
  21 : *
  22 :    MODL1 = MODL TOTAL MECANIQUE ELASTIQUE CU20 ;
  23 :     CARM = MATR MODL1 'YOUN' 2.E11 'NU' 0.3 'RHO' 7800. 'ALPH' 12.E-6 ;
  24 : *
  25 : * Conditions de blocages
  26 : *
  27 :                         ENC1 = BLOQ UZ SUB ;
  28 :                ENC2 = BLOQ UY LI ; ENC3 = BLOQ UX P1 ;
  29 : *
  30 : * Blocage pour deplacements imposes
  31 : *
  32 :                         ENC4 = BLOQ UZ SUH ;
  33 :                     ENC = ENC1 ET ENC2 ET ENC3 ;
  34 : *
  35 : * Raideur totale
  36 : *
  37 :                  RIG = (RIGI MODL1 CARM) ET ENC ET ENC4 ;
  38 : *
  39 : * Valeur des deplacements imposes
  40 : *
  41 :                       FEXT = DEPI ENC4 1E-9  ;
  42 : *
  43 : * Force ponctuelle et force totale
  44 : *
  45 :          FNOD = FORCE LIH (0 1000 0) ; FTOT = FEXT ET FNOD ;
  46 : *
  47 : * Resolution  pour obtenir les deplacements
  48 : *
  49 :                         DEP = RESOU RIG FTOT ;
  50 : *
  51 : * Vecteur force et reaction
  52 : *
  53 :                VEF = VECTEUR FTOT 0.03 FX FY FZ ROUG ;
  54 :                    REA = REAC (ENC4 ET ENC) DEP ;
  55 :                VER = VECTEUR REA 0.01 FX FY FZ BLEU ;
  56 : *
  57 : * Contraintes et vonmises
  58 : *
  59 :                SSI = SIGMA MODL1 DEP CARM ; VM = VMIS MODL1 SSI ;
  60 : *
  61 : * Deformee et trace
  62 : *
  63 :              DEF = DEFO TOTAL DEP 50 (VEF ET VER) VERT ;
  64 :                    TRAC DEF MODL1 VM (1000 -2000 1000) ;
  65 :                                 FIN ;
  66 :