An analysis was attempted based on previously reported analysis examples.
The following tools were used.
OS used: Windows 11
Soft SALOME-MECA 2023 for Windows
Request:
As described below, the calculation is not performed, so please advise how to resolve it.
Model
See entity column
Importing model
The following step file was imported test-contact-1.stp
Creating an entity
Creating a mesh
An error occurred when creating a mesh for each object as a single unit, but it was resolved by using a submesh.
Creating the analysis code
I tried to convert the analysis code to 2023 using the past example (2009) as a reference. Shown below.
Calculation results
An example of the situation is shown below.
The program will be aborted due to an error reading the output file.
I tried linear analysis with a similar model, but the output file was aborted when loading, just like the above calculation results.
Analysis code
<1>
Test contact-1 by Fujii mod. ソースコード 2025-03-27
FIN()DEBUT(LANG=‘EN’)
MAIL = LIRE_MAILLAGE(identifier=‘0:1’,
FORMAT=‘MED’,
UNITE=20)
MAIL = MODI_MAILLAGE(identifier=‘1:1’,
reuse=MAIL,
MAILLAGE=MAIL,
ORIE_PEAU=_F(GROUP_MA_PEAU=(‘contPL’, ‘contBASE’, ‘PRESS’)))
MODE = AFFE_MODELE(identifier=‘2:1’,
AFFE=_F(MODELISATION=‘3D’,
PHENOMENE=‘MECANIQUE’,
TOUT=‘OUI’),
MAILLAGE=MAIL)
MAT = DEFI_MATERIAU(identifier=‘3:1’,
ELAS=_F(E=130300.0,
NU=0.343))
MATE = AFFE_MATERIAU(identifier=‘4:1’,
AFFE=_F(MATER=(MAT, ),
TOUT=‘OUI’),
MODELE=MODE)
INST = DEFI_LIST_REEL(identifier=‘5:1’,
DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1.0,
PAS=0.2))
RAMP = DEFI_FONCTION(identifier=‘6:1’,
NOM_PARA=‘INST’,
VALE=(0.0, 0.0, 1.0, 1.0))
times = DEFI_LIST_INST(identifier=‘7:1’,
DEFI_LIST=_F(LIST_INST=INST))
CHAR = AFFE_CHAR_MECA(identifier=‘8:1’,
DDL_IMPO=(_F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=(‘FIX’, )),
_F(DX=0.0,
DZ=0.0,
GROUP_MA=(‘PRESS’, ))),
MODELE=MODE)
loadP = AFFE_CHAR_MECA(identifier=‘9:1’,
DDL_IMPO=_F(DY=-0.5,
GROUP_MA=(‘PRESS’, )),
MODELE=MODE)
contact = DEFI_CONTACT(identifier=‘10:1’,
ALGO_RESO_CONT=‘NEWTON’,
FORMULATION=‘CONTINUE’,
MODELE=MODE,
ZONE=_F(ALGO_CONT=‘STANDARD’,
APPARIEMENT=‘MAIT_ESCL’,
COEF_CONT=100.0,
CONTACT_INIT=‘OUI’,
GROUP_MA_ESCL=(‘contBASE’, ),
GROUP_MA_MAIT=(‘contPL’, )))
result = STAT_NON_LINE(identifier=‘11:1’,
ARCHIVAGE=_F(PAS_ARCH=1),
CHAM_MATER=MATE,
COMPORTEMENT=_F(DEFORMATION=‘GREEN_LAGRANGE’,
RELATION=‘ELAS’,
TOUT=‘OUI’),
CONTACT=contact,
CONVERGENCE=_F(ITER_GLOB_MAXI=30),
EXCIT=(_F(CHARGE=CHAR),
_F(CHARGE=loadP,
FONC_MULT=RAMP)),
INCREMENT=_F(LIST_INST=INST),
METHODE=‘NEWTON’,
MODELE=MODE,
NEWTON=_F(MATRICE=‘TANGENTE’,
REAC_INCR=1,
REAC_ITER=1))
unnamed = CALC_CHAMP(identifier=‘12:1’,
CONTRAINTE=(‘SIGM_ELNO’, ‘SIGM_NOEU’),
CRITERES=(‘SIEQ_ELNO’, ‘SIEQ_NOEU’),
DEFORMATION=(‘EPSI_NOEU’, ),
FORCE=(‘REAC_NODA’, ),
MODELE=MODE,
RESULTAT=result)
IMPR_RESU(identifier=‘13:1’,
RESU=_F(MAILLAGE=MAIL,
NOM_CHAM=(‘DEPL’, ‘SIEQ_NOEU’, ‘SIGM_NOEU’),
RESULTAT=result),
UNITE=2)
FIN()
References
2025-03-30T15:00:00Z