Using CAble Element to simulate the swinging of a simple pendulum
Hello there, I am trying to simulate a pendulum using a cable element and it is giving me this error:
════════════════════════════════════════════════════════════════════════════════════════════════╗
║ <EXCEPTION> <DVP_1> ║
║ ║
║ Erreur de programmation. ║
║ ║
║ Condition non respectée: ║
║ .false. ║
║ Fichier ║
║ /local01/home/es4c769n/Constructions/meca/V2022_scibian_univ/tools/src/Code_aster_testing-1620 ║
║ /bibfor/debug/contex_param.F90, ligne 147 ║
║ ║
║ -------------------------------------------- ║
║ Contexte du message : ║
║ Option : CHAR_MECA_PESA_R ║
║ Type d'élément : MECABL2 ║
║ Maillage : 00000001 ║
║ Maille : M1 ║
║ Type de maille : SEG2 ║
║ Cette maille appartient aux groupes de mailles suivants : ║
║ Group_1 ║
║ Position du centre de gravité de la maille : ║
║ x=0.000000 y=-0.000000 z=0.207500 ║
║ ║
║ ║
║ Il y a probablement une erreur dans la programmation. ║
║ Veuillez contacter votre assistance technique. ║
╚══════════════════════════════════════════════
════════════════════════════════════════════════════════════════════════════════════════════════╗
║ <EXCEPTION> <DVP_1> ║
║ ║
║ Erreur de programmation. ║
║ ║
║ Condition non respectée: ║
║ .false. ║
║ Fichier ║
║ /local01/home/es4c769n/Constructions/meca/V2022_scibian_univ/tools/src/Code_aster_testing-1620 ║
║ /bibfor/debug/contex_param.F90, ligne 147 ║
║ ║
║ ║
║ Il y a probablement une erreur dans la programmation. ║
║ Veuillez contacter votre assistance technique.
please can anyone help me solve this problem
This is my input code:
DEBUT(LANG='EN')
mesh = LIRE_MAILLAGE(UNITE=3)
model = AFFE_MODELE(AFFE=_F(MODELISATION='CABLE',
PHENOMENE='MECANIQUE',
TOUT='OUI'),
MAILLAGE=mesh)
elemprop = AFFE_CARA_ELEM(CABLE=_F(GROUP_MA=('Group_1', ),
N_INIT=10.0,
SECTION=0.002828571),
MODELE=model)
mater = DEFI_MATERIAU(CABLE=F(EC_SUR_E=0.0001),
ELAS=F(E=90000000000.0,
NU=0.36,
RHO=8e-09))
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
TOUT='OUI'),
MAILLAGE=mesh,
MODELE=model)
#comment: times0 = DEFI_LIST_INST(DEFI_LIST=_F(VALE=(0.0, 2.0)))
listr = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=10.0,
NOMBRE=1))
load = AFFE_CHAR_MECA(DDL_IMPO=F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('Fix', )),
FORCE_NODALE=F(FZ=-20.0,
GROUP_NO=('force', )),
MODELE=model)
load0 = AFFE_CHAR_MECA(MODELE=model,
PESANTEUR=_F(DIRECTION=(0.0, 0.0, -1.0),
GRAVITE=9.81))
#comment: reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
#comment: CHAM_MATER=fieldmat,
#comment: EXCIT=_F(CHARGE=load),
#comment: MODELE=model)
resnonl = STAT_NON_LINE(CHAM_MATER=fieldmat,
COMPORTEMENT=F(DEFORMATION='GROT_GDEP',
RELATION='CABLE',
TOUT='OUI'),
CONVERGENCE=F(ITER_GLOB_MAXI=300),
EXCIT=(_F(CHARGE=load),
F(CHARGE=load0)),
INCREMENT=F(LIST_INST=listr),
MODELE=model,
NEWTON=_F(MATRICE='TANGENTE',
PREDICTION='TANGENTE',
REAC_ITER=1))
resnonl = CALC_CHAMP(reuse=resnonl,
CONTRAINTE=('SIGM_NOEU', 'SIGM_ELNO'),
CRITERES=('SIEQ_NOEU', ),
MODELE=model,
RESULTAT=resnonl)
IMPR_RESU(RESU=_F(MAILLAGE=mesh,
RESULTAT=resnonl),
UNITE=80)
FIN()
