Difference between revisions of "Contrib:asamit"

From CAELinuxWiki
Jump to: navigation, search
Line 13: Line 13:
  
 
1.
 
1.
[[http://www.example.com link title '''Dounload Mesh''']]
 
 
<code>
 
<code>
 
  DEBUT();
 
  DEBUT();

Revision as of 14:18, 7 April 2010

Asamit-roller-2.jpg Asamit-roller.jpg Asamit-roller-3.jpg

Problems

1. How to rotate the magenta object 360 degrees.

2. How to connect the magenta and the yellow.

3. How to add rigid bars(spokes) to the yellow.



1.

DEBUT();
MA1=LIRE_MAILLAGE(UNITE=20,
                 FORMAT='MED',);
MA1=DEFI_GROUP(reuse =MA1,
              MAILLAGE=MA1,
              CREA_GROUP_MA=_F(NOM='TOUT',
                               TOUT='OUI',),);
MA2=CREA_MAILLAGE(MAILLAGE=MA1,
                 MODI_MAILLE=_F(TOUT='OUI',
                                OPTION='QUAD8_9',),);
MO=AFFE_MODELE(MAILLAGE=MA2,
              AFFE=_F(TOUT='OUI',
                      PHENOMENE='MECANIQUE',
                      MODELISATION='COQUE_3D',),);
CARA_ELE=AFFE_CARA_ELEM(MODELE=MO,
                       COQUE=_F(GROUP_MA='TOUT',
                                EPAIS=1.0,),);
STEEL=DEFI_MATERIAU(ELAS=_F(E=200000.0,
                           NU=.3,
                           RHO=7.8e-9,),);
MAT1=AFFE_MATERIAU(MAILLAGE=MA2,
                  AFFE=_F(TOUT='OUI',
                          MATER=STEEL,),);
MECA=AFFE_CHAR_MECA(MODELE=MO,
                   DDL_IMPO=_F(GROUP_MA='ax_shaft',
                               DX=0.0,
                               DY=0.0,
                               DZ=0.0,
                               DRX=6.6,
                               DRY=0.0,
                               DRZ=0.0,),);
L_INST=DEFI_LIST_REEL(DEBUT=0.0,
                     INTERVALLE=_F(JUSQU_A=1,
                                   NOMBRE=10,),);
inst=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO',
                                LIST_INST=L_INST,
                                PAS_MINI=1e-8,),);
RAMP=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0,
                        1.0,1.0,
                        ),);
STA=STAT_NON_LINE(MODELE=MO,
                 CHAM_MATER=MAT1,
                 CARA_ELEM=CARA_ELE,
                 EXCIT=_F(CHARGE=MECA,
                          FONC_MULT=RAMP,),
                 COMP_INCR=_F(RELATION='ELAS',
                              DEFORMATION='GROT_GDEP',),
                 INCREMENT=_F(LIST_INST=inst,),
                 NEWTON=_F(PREDICTION='ELASTIQUE',
                           MATRICE='TANGENTE',
                           REAC_ITER=1,),);
STA=CALC_ELEM(reuse =STA,
             MODELE=MO,
             CHAM_MATER=MAT1,
             RESULTAT=STA,
             OPTION='EQUI_ELNO_SIGM',);
IMPR_RESU(FORMAT='MED',
         UNITE=70,
         RESU=_F(MAILLAGE=MA2,
                 RESULTAT=STA,
                 TOUT_CHAM='OUI',),);
FIN();