Difference between revisions of "Contrib:asamit"
From CAELinuxWiki
Line 14: | Line 14: | ||
1.[[http://www.caelinux.org/wiki/images/2/26/Asamit-rota-mesh.tar.gz '''Download Mesh''']] | 1.[[http://www.caelinux.org/wiki/images/2/26/Asamit-rota-mesh.tar.gz '''Download Mesh''']] | ||
<code> | <code> | ||
− | DEBUT(); | + | DEBUT(); |
− | + | MA1=LIRE_MAILLAGE(UNITE=20, | |
− | MA1=LIRE_MAILLAGE(UNITE=20, | + | |
FORMAT='MED',); | FORMAT='MED',); | ||
− | + | MA2=CREA_MAILLAGE(MAILLAGE=MA1, | |
− | MA2=CREA_MAILLAGE(MAILLAGE=MA1, | + | |
MODI_MAILLE=_F(TOUT='OUI', | MODI_MAILLE=_F(TOUT='OUI', | ||
OPTION='QUAD8_9',),); | OPTION='QUAD8_9',),); | ||
− | + | MA3=CREA_MAILLAGE(MAILLAGE=MA2, | |
− | MA3=CREA_MAILLAGE(MAILLAGE=MA2, | + | |
CREA_POI1=_F(NOM_GROUP_MA='poi1', | CREA_POI1=_F(NOM_GROUP_MA='poi1', | ||
GROUP_NO='a_shaft',),); | GROUP_NO='a_shaft',),); | ||
− | + | MO=AFFE_MODELE(MAILLAGE=MA3, | |
− | MO=AFFE_MODELE(MAILLAGE=MA3, | + | |
AFFE=(_F(GROUP_MA='shaft', | AFFE=(_F(GROUP_MA='shaft', | ||
PHENOMENE='MECANIQUE', | PHENOMENE='MECANIQUE', | ||
Line 34: | Line 30: | ||
PHENOMENE='MECANIQUE', | PHENOMENE='MECANIQUE', | ||
MODELISATION='DIS_TR',),),); | MODELISATION='DIS_TR',),),); | ||
− | + | CARA_ELE=AFFE_CARA_ELEM(MODELE=MO, | |
− | CARA_ELE=AFFE_CARA_ELEM(MODELE=MO, | + | |
COQUE=_F(GROUP_MA='shaft', | COQUE=_F(GROUP_MA='shaft', | ||
EPAIS=1.0,), | EPAIS=1.0,), | ||
Line 42: | Line 37: | ||
GROUP_MA='poi1', | GROUP_MA='poi1', | ||
VALE=(1,1,1,1,1,1,),),); | VALE=(1,1,1,1,1,1,),),); | ||
− | + | STEEL=DEFI_MATERIAU(ELAS=_F(E=200000.0, | |
− | STEEL=DEFI_MATERIAU(ELAS=_F(E=200000.0, | + | |
NU=.3, | NU=.3, | ||
RHO=7.8e-9,),); | RHO=7.8e-9,),); | ||
− | + | MAT1=AFFE_MATERIAU(MAILLAGE=MA3, | |
− | MAT1=AFFE_MATERIAU(MAILLAGE=MA3, | + | |
AFFE=_F(TOUT='OUI', | AFFE=_F(TOUT='OUI', | ||
MATER=STEEL,),); | MATER=STEEL,),); | ||
− | + | MECA=AFFE_CHAR_MECA(MODELE=MO, | |
− | MECA=AFFE_CHAR_MECA(MODELE=MO, | + | |
DDL_IMPO=_F(GROUP_MA='poi1', | DDL_IMPO=_F(GROUP_MA='poi1', | ||
DX=0.0, | DX=0.0, | ||
Line 60: | Line 52: | ||
DRZ=0.0,), | DRZ=0.0,), | ||
LIAISON_SOLIDE=_F(GROUP_MA=('ax_shaft','poi1',),),); | LIAISON_SOLIDE=_F(GROUP_MA=('ax_shaft','poi1',),),); | ||
− | + | L_INST=DEFI_LIST_REEL(DEBUT=0.0, | |
− | L_INST=DEFI_LIST_REEL(DEBUT=0.0, | + | |
INTERVALLE=_F(JUSQU_A=1.0, | INTERVALLE=_F(JUSQU_A=1.0, | ||
NOMBRE=10,),); | NOMBRE=10,),); | ||
− | + | inst=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO', | |
− | inst=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO', | + | |
LIST_INST=L_INST, | LIST_INST=L_INST, | ||
PAS_MINI=1e-8,),); | PAS_MINI=1e-8,),); | ||
− | + | RAMP=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0, | |
− | RAMP=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0, | + | |
1.0,1.0, | 1.0,1.0, | ||
),); | ),); | ||
− | + | STA=STAT_NON_LINE(MODELE=MO, | |
− | STA=STAT_NON_LINE(MODELE=MO, | + | |
CHAM_MATER=MAT1, | CHAM_MATER=MAT1, | ||
CARA_ELEM=CARA_ELE, | CARA_ELEM=CARA_ELE, | ||
Line 86: | Line 74: | ||
SOLVEUR=_F(METHODE='MULT_FRONT', | SOLVEUR=_F(METHODE='MULT_FRONT', | ||
STOP_SINGULIER='NON',),); | STOP_SINGULIER='NON',),); | ||
− | + | STA=CALC_ELEM(reuse =STA, | |
− | STA=CALC_ELEM(reuse =STA, | + | |
MODELE=MO, | MODELE=MO, | ||
CHAM_MATER=MAT1, | CHAM_MATER=MAT1, | ||
RESULTAT=STA, | RESULTAT=STA, | ||
OPTION='EQUI_ELNO_SIGM',); | OPTION='EQUI_ELNO_SIGM',); | ||
− | + | IMPR_RESU(FORMAT='MED', | |
− | IMPR_RESU(FORMAT='MED', | + | |
UNITE=70, | UNITE=70, | ||
RESU=_F(MAILLAGE=MA3, | RESU=_F(MAILLAGE=MA3, | ||
RESULTAT=STA, | RESULTAT=STA, | ||
TOUT_CHAM='OUI',),); | TOUT_CHAM='OUI',),); | ||
− | + | FIN(); | |
− | FIN(); | + | |
</code> | </code> |
Latest revision as of 11:17, 10 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.
1.[Download Mesh]
DEBUT(); MA1=LIRE_MAILLAGE(UNITE=20, FORMAT='MED',); MA2=CREA_MAILLAGE(MAILLAGE=MA1, MODI_MAILLE=_F(TOUT='OUI', OPTION='QUAD8_9',),); MA3=CREA_MAILLAGE(MAILLAGE=MA2, CREA_POI1=_F(NOM_GROUP_MA='poi1', GROUP_NO='a_shaft',),); MO=AFFE_MODELE(MAILLAGE=MA3, AFFE=(_F(GROUP_MA='shaft', PHENOMENE='MECANIQUE', MODELISATION='COQUE_3D',), _F(GROUP_MA='poi1', PHENOMENE='MECANIQUE', MODELISATION='DIS_TR',),),); CARA_ELE=AFFE_CARA_ELEM(MODELE=MO, COQUE=_F(GROUP_MA='shaft', EPAIS=1.0,), DISCRET=_F(REPERE='GLOBAL', CARA='K_TR_D_N', GROUP_MA='poi1', VALE=(1,1,1,1,1,1,),),); STEEL=DEFI_MATERIAU(ELAS=_F(E=200000.0, NU=.3, RHO=7.8e-9,),); MAT1=AFFE_MATERIAU(MAILLAGE=MA3, AFFE=_F(TOUT='OUI', MATER=STEEL,),); MECA=AFFE_CHAR_MECA(MODELE=MO, DDL_IMPO=_F(GROUP_MA='poi1', DX=0.0, DY=0.0, DZ=0.0, DRX=6.28, DRY=0.0, DRZ=0.0,), LIAISON_SOLIDE=_F(GROUP_MA=('ax_shaft','poi1',),),); L_INST=DEFI_LIST_REEL(DEBUT=0.0, INTERVALLE=_F(JUSQU_A=1.0, 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,), SOLVEUR=_F(METHODE='MULT_FRONT', STOP_SINGULIER='NON',),); 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=MA3, RESULTAT=STA, TOUT_CHAM='OUI',),); FIN();