Difference between revisions of "Contrib:asamit"
From CAELinuxWiki
| Line 12: | Line 12: | ||
---- | ---- | ||
| − | 1.[[http://www.caelinux.org/wiki/images/e/eb/Asamit-shaft.med.tar.gz ''' | + | 1.[[http://www.caelinux.org/wiki/images/e/eb/Asamit-shaft.med.tar.gz '''Download MED''']] [[http://www.caelinux.org/wiki/images/1/16/Asamit-rota.export.tar.gz '''Download export''']] |
<code> | <code> | ||
DEBUT(); | DEBUT(); | ||
Revision as of 14:13, 8 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 MED] [Download export]
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();
