3d.comm

DEBUT();
mesh=LIRE_MAILLAGE( INFO=1,UNITE=20,FORMAT='MED',);
mesh=DEFI_GROUP(
 MAILLAGE=mesh,
 reuse =mesh,
 CREA_GROUP_MA=_F(NOM='TOUT',TOUT='OUI',),
 CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI',),
);
model=AFFE_MODELE(
MAILLAGE=mesh,
 AFFE=(
 _F(
 GROUP_MA=('poutre',),
 PHENOMENE='MECANIQUE',
 MODELISATION='3D',
 ),
 ),
);
steel=DEFI_MATERIAU(
 ELAS=_F(E=210000.,NU=0.3,RHO=8e-9,ALPHA=12e-6,),
);
material=AFFE_MATERIAU(
 MAILLAGE=mesh,
 AFFE=_F(GROUP_MA=('TOUT',),MATER=steel,),
);
encas=AFFE_CHAR_MECA(
 MODELE=model,
 DDL_IMPO=
 _F(GROUP_NO=('encastr',),DX=0,DY=0,DZ=0,),
);
force=AFFE_CHAR_MECA(
 MODELE=model,
 FORCE_NODALE=_F(GROUP_NO=('force',),FZ=-500,),
);
stat=MECA_STATIQUE(
 MODELE=model,
 CHAM_MATER=material,
 INFO=2,
 EXCIT=(
 _F(CHARGE=encas,),
 _F(CHARGE=force,),
 ),
);
stat=CALC_CHAMP(
 reuse =stat,
 RESULTAT=stat,
 CONTRAINTE=(
 'SIEF_ELNO',
 'SIPM_ELNO',
 'SIGM_ELNO',
 ),
 FORCE=('REAC_NODA'),
 CRITERES=('SIEQ_NOEU'),
);
IMPR_RESU(
 FORMAT='MED', 
 UNITE=80,
 RESU=(
 _F(
 GROUP_MA=('poutre'),
 RESULTAT=stat,NOM_CHAM=('DEPL',),
 ),
 _F(
 GROUP_MA=('poutre',),
 RESULTAT=stat,
 NOM_CHAM='SIEQ_NOEU',NOM_CMP='VMIS',
 NOM_CHAM_MED='vmis',
 ),
 ),
);
FIN();