Static structural analysis for rotating wheel

DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=2.1e+11,
                         NU=0.3,
                         RHO=7800,),
                      );    

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse=MAIL,
                   MAILLAGE=MAIL,
                   ORIE_PEAU_3D=_F(GROUP_MA=('fix',),),
                   );

# create node groups fom element groups 'SLAVE' and 'fix' with the same names
MAIL=DEFI_GROUP(reuse=MAIL,
		MAILLAGE=MAIL,
		CREA_GROUP_NO=_F(GROUP_MA=('SLAVE','fix')),
		);

# create a node group which is the difference of SLAVE with fix named SLAVE_OK
MAIL=DEFI_GROUP(reuse=MAIL,
		MAILLAGE=MAIL,
		CREA_GROUP_NO=_F(DIFFE=('SLAVE','fix'), NOM='SLAVE_OK'),
		);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
                 AFFE=_F(TOUT='OUI',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=_F(TOUT='OUI',
                           MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    ROTATION=_F(GROUP_MA='rot',
                                VITESSE=104.7,
                                AXE=(0,0,1,),),
                    DDL_IMPO=_F(GROUP_MA=('fix',),
                                DX=0.0,
                                DY=0.0,
                                DZ=0.0,),);

CYCL=AFFE_CHAR_MECA(MODELE=MODE,
                    LIAISON_MAIL=_F(GROUP_MA_MAIT='rot',
                                    GROUP_NO_ESCL='SLAVE_OK',
                                    ANGL_NAUT=(18.0,0.0,0.0,),
                                    CENTRE=(0.0,0.0,0.0,),),
                    );	

RESU=MECA_STATIQUE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=(_F(CHARGE=CHAR,),
			  _F(CHARGE=CYCL,),),);

RESU=CALC_CHAMP(reuse =RESU,
               MODELE=MODE,
               CHAM_MATER=MATE,
               RESULTAT=RESU,
               CONTRAINTE=('SIGM_NOEU','SIGM_ELGA','SIGM_ELNO',),
               CRITERES=('SIEQ_NOEU','SIEQ_ELGA','SIEQ_ELNO',),
               );

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL','SIGM_ELGA','SIGM_ELNO',),),);

FIN();