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