Hello, everyone! I would like to ask a question: how can I restart a simulation that was stopped at a specific point? I am writing the script for the .comm file and have a slightly more complex case, but it’s similar to what I'm describing here. Typically, the simulation halts during the modal calculation, which suggests it should restart with the mass, stiffness, and geometric matrices already defined, beginning from that point — specifically, after the 'Assemblage' stage.
I researched on the forum and found a command called POURSUITE(), which seems suitable for this purpose. However, when I tried to restart the simulation, it didn’t recognize the K and Kt variables, which should be loaded from the interrupted case. Therefore, I believe these matrices need to be saved first so they can be loaded later.
If anyone has knowledge of the correct procedure, I would greatly appreciate your assistance. Thank you very much!
`
mesh = LIRE_MAILLAGE(FORMAT='MED',
UNITE=20)
model = AFFE_MODELE(AFFE=(_F(GROUP_MA=('gp_msh_0D_support', ),
MODELISATION='DIS_TR',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('gp_msh_pipe_1', 'gp_msh_pipe_10', 'gp_msh_pipe_11', 'gp_msh_pipe_12', 'gp_msh_pipe_13', 'gp_msh_pipe_14', 'gp_msh_pipe_15', 'gp_msh_pipe_16', 'gp_msh_pipe_17', 'gp_msh_pipe_18', 'gp_msh_pipe_19', 'gp_msh_pipe_2', 'gp_msh_pipe_20', 'gp_msh_pipe_21', 'gp_msh_pipe_3', 'gp_msh_pipe_4', 'gp_msh_pipe_5', 'gp_msh_pipe_6', 'gp_msh_pipe_7', 'gp_msh_pipe_8', 'gp_msh_pipe_9'),
MODELISATION='DKT',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('gp_msh_2D_wall', ),
MODELISATION='DKT',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh)
elemprop = AFFE_CARA_ELEM(COQUE=(_F(EPAIS=0.03,
GROUP_MA=('gp_msh_2D_wall', ),
VECTEUR=(0.0, 1.0, 0.0)),
F(EPAIS=0.0045,
GROUP_MA=('gp_msh_pipe_1', 'gp_msh_pipe_10', 'gp_msh_pipe_11', 'gp_msh_pipe_12', 'gp_msh_pipe_13', 'gp_msh_pipe_14', 'gp_msh_pipe_15', 'gp_msh_pipe_16', 'gp_msh_pipe_17', 'gp_msh_pipe_18', 'gp_msh_pipe_19', 'gp_msh_pipe_2', 'gp_msh_pipe_20', 'gp_msh_pipe_21', 'gp_msh_pipe_3', 'gp_msh_pipe_4', 'gp_msh_pipe_5', 'gp_msh_pipe_6', 'gp_msh_pipe_7', 'gp_msh_pipe_8', 'gp_msh_pipe_9'))),
DISCRET=F(CARA='K_TR_D_N',
GROUP_MA=('gp_msh_0D_support', ),
VALE=(1e+18, 1e+18, 1e+18, 1e+18, 1e+18, 1e+18)),
MODELE=model)
steel_A36_20Cdeg = DEFI_MATERIAU(ELAS=_F(E=200000000000.0,
NU=0.3,
RHO=7800.0))
steel_A36_250Cdeg = DEFI_MATERIAU(ELAS=_F(E=184561557270.0,
NU=0.3,
RHO=7800.0))
steel_A36_500Cdeg = DEFI_MATERIAU(ELAS=_F(E=145873689330.0,
NU=0.3,
RHO=7800.0))
material = AFFE_MATERIAU(AFFE=_F(GROUP_MA=('gp_msh_2D_wall', 'gp_msh_pipe_1', 'gp_msh_pipe_10', 'gp_msh_pipe_11', 'gp_msh_pipe_12', 'gp_msh_pipe_13', 'gp_msh_pipe_14', 'gp_msh_pipe_15', 'gp_msh_pipe_16', 'gp_msh_pipe_17', 'gp_msh_pipe_18', 'gp_msh_pipe_19', 'gp_msh_pipe_2', 'gp_msh_pipe_20', 'gp_msh_pipe_21', 'gp_msh_pipe_3', 'gp_msh_pipe_4', 'gp_msh_pipe_5', 'gp_msh_pipe_6', 'gp_msh_pipe_7', 'gp_msh_pipe_8', 'gp_msh_pipe_9'),
MATER=(steel_A36_500Cdeg, )),
MAILLAGE=mesh)
pressure = AFFE_CHAR_MECA(MODELE=model,
PRES_REP=_F(GROUP_MA='gp_msh_pipe_1',
PRES=3825000.0),
VERI_NORM='NON')
reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
CHAM_MATER=material,
EXCIT=F(CHARGE=pressure),
MODELE=model,
SOLVEUR=F(RESI_RELA=0.0001))
prestress = CREA_CHAMP(NOM_CHAM='SIEF_ELGA',
NUME_ORDRE=1,
OPERATION='EXTR',
RESULTAT=reslin,
TYPE_CHAM='ELGA_SIEF_R')
ASSEMBLAGE(CARA_ELEM=elemprop,
CHAM_MATER=material,
MATR_ASSE=(_F(MATRICE=CO('K'),
OPTION='RIGI_MECA'),
_F(MATRICE=CO('Kg'),
OPTION='RIGI_GEOM',
SIEF_ELGA=prestress),
_F(MATRICE=CO('M'),
OPTION='MASS_MECA')),
MODELE=model,
NUME_DDL=CO('nodfs'))
nodal_force = CALC_CHAMP(FORCE=('FORC_NODA', ),
MODELE=model,
RESULTAT=reslin,
TOUT_ORDRE='OUI')
IMPR_RESU(RESU=_F(RESULTAT=reslin),
UNITE=2)
Kt = COMB_MATR_ASSE(COMB_R=(_F(COEF_R=1.0,
MATR_ASSE=Kg),
_F(COEF_R=1.0,
MATR_ASSE=K)))
modes = CALC_MODES(CALC_FREQ=F(NMAX_FREQ=20),
MATR_MASS=M,
MATR_RIGI=Kt,
OPTION='PLUS_PETITE',
SOLVEUR_MODAL=F(METHODE='TRI_DIAG',
MODE_RIGIDE='NON'))
IMPR_RESU(RESU=_F(RESULTAT=modes),
UNITE=80)
`
