Difference between revisions of "Contrib:asamit"

From CAELinuxWiki
Jump to: navigation, search
Line 10: Line 10:
  
  
----
+
<code>
 +
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();
 +
</code>

Revision as of 10:52, 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.


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();