Dear all,
I downloaded salome_meca a few days ago and tried my first calculations, mostly according to youtube tutorials, to learn how to use the program. I found the following interesting tutorial which I tried to reproduce: Contact non linear analysis using code aster|salome meca tutorial|paraview tuto|tut-49. If I do it with exactly the provided instructions in the tutorial, it works. If I do it with my own files I always get errors.
I have now some initial questions:
- In the Tutorial, at 2:14, the STEP-import was closed with "No", which means, the units from the file are set to the system settings (which is normally "m" for meters") is this correct? So, if the model was built in millimeters, the model dimensions are converted to meter, but the size is unchanged. Is this correct?
- In the Tutorial, at 11:42 the Young's Module is set to 2.1e5 which means 210'000N/mm2. If the system is set to meters, should there be the value in Pascal, means 210'000'000'000?
If I import my own model (
), mesh it and add the following codes, the solver ends in errors:
DEBUT()
mesh = LIRE_MAILLAGE(FORMAT='MED',
UNITE=20)
model = AFFE_MODELE(AFFE=_F(MODELISATION='3D',
PHENOMENE='MECANIQUE',
TOUT='OUI'),
MAILLAGE=mesh)
mater = DEFI_MATERIAU(ELAS=_F(E=210000000000.0,
NU=0.3))
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
TOUT='OUI'),
MAILLAGE=mesh)
Load_Function = DEFI_FONCTION(NOM_PARA='INST',
VALE=(0.0, 0.0, 1.0, 1.0))
listr = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1.0,
PAS=0.1))
times = DEFI_LIST_INST(DEFI_LIST=_F(LIST_INST=listr),
METHODE='AUTO')
load = AFFE_CHAR_MECA(DDL_IMPO=_F(BLOCAGE=('DEPLACEMENT', ),
GROUP_MA=('fixation', )),
MODELE=model)
load0 = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=-0.009,
DY=0.0,
DZ=0.0,
GROUP_MA=('displacement', )),
MODELE=model)
contact = DEFI_CONTACT(FORMULATION='CONTINUE',
MODELE=model,
ZONE=_F(GROUP_MA_ESCL=('slave', ),
GROUP_MA_MAIT=('master', )))
resnonl = STAT_NON_LINE(CHAM_MATER=fieldmat,
COMPORTEMENT=_F(DEFORMATION='PETIT',
RELATION='ELAS'),
CONTACT=contact,
EXCIT=(_F(CHARGE=load),
_F(CHARGE=load0,
FONC_MULT=Load_Function)),
INCREMENT=_F(LIST_INST=listr),
MODELE=model)
IMPR_RESU(FORMAT='MED',
RESU=_F(MAILLAGE=mesh,
RESULTAT=resnonl,
TOUT_CHAM='OUI'),
UNITE=2)
FIN()
Does anyone have any idea what could be wrong?
