Difference between revisions of "Contrib:asamit"
From CAELinuxWiki
| 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 11:52, 7 April 2010
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();
