jeanpierreaubry
mesh = LIRE_MAILLAGE(UNITE=20)
model = AFFE_MODELE(AFFE=(_F(MODELISATION='3D',
PHENOMENE='MECANIQUE',
TOUT='OUI'),
_F(GROUP_MA=('Pt_pilote2', ),
MODELISATION='DIS_TR',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh)
elemprop = AFFE_CARA_ELEM(DISCRET=(_F(CARA='K_TR_D_N',
GROUP_MA=('Pt_pilote2', ),
VALE=(1.0, 1.0, 1.0, 1.0, 1.0, 1.0)),
_F(CARA='M_TR_D_N',
GROUP_MA=('Pt_pilote2', ),
VALE=(1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0))),
MODELE=model)
spep = DEFI_MATERIAU(ELAS=_F(E=195000000000.0,
NU=0.3,
RHO=7850.0))
beton = DEFI_MATERIAU(ELAS=_F(E=35077000000.0,
NU=0.2,
RHO=2500.0))
fieldmat = AFFE_MATERIAU(AFFE=(_F(GROUP_MA=('solid_beton', ),
MATER=(beton, )),
_F(GROUP_MA=('Grspep_Volumes', ),
MATER=(spep, ))),
MODELE=model)
load = AFFE_CHAR_MECA(FACE_IMPO=F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('face_fix', )),
FORCE_NODALE=F(FX=101000.0,
FY=1000.0,
FZ=15000.0,
GROUP_NO=('pt_pilote', )),
LIAISON_RBE3=F(COEF_ESCL=(1.0, ),
DDL_ESCL=('DX-DY-DZ', ),
DDL_MAIT=('DX', 'DY', 'DZ'),
GROUP_NO_ESCL=('force_noeuds', ),
GROUP_NO_MAIT=('pt_pilote', )),
MODELE=model,
PESANTEUR=F(DIRECTION=(-1.0, 0.0, 0.0),
GRAVITE=9.81,
GROUP_MA=('solid_beton', 'Grspep_Volumes')))
contact = AFFE_CHAR_MECA(LIAISON_MAIL=_F(GROUP_MA_ESCL=('conct_spep_1_rotated', ),
GROUP_MA_MAIT=('solid_beton', ),
TYPE_RACCORD='MASSIF'),
MODELE=model)
reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
EXCIT=(_F(CHARGE=load),
_F(CHARGE=contact)),
MODELE=model)
IMPR_RESU(RESU=_F(CARA_ELEM=elemprop,
RESULTAT=reslin,
TOUT_CHAM='OUI'),
UNITE=2)
