Let me try to simplify the description of the error.
Following the suggestion of Sound_Spinning, this is the model:
https://drive.google.com/file/d/1CwmqcsLH7PtwAEx2R36OADuJkd-KPPFD/view?usp=drive_link
In SSNL118 a result is read and, since I don't have one for my case, I used the commands below to create one:
vvent0 = CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='AFFE',
INFO=2,
MODELE=model,
AFFE=(_F(TOUT='OUI',
NOM_CMP=('DX', 'DY', 'DZ'),
VALE=(0.0, 0.0, 0.0),),
F(GROUP_NO=('body', ),
NOM_CMP=('DX', 'DY', 'DZ'),
VALE=(0.0, 0.0, 0.0),),
),
)
vvent1 = CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='AFFE',
INFO=2,
MODELE=model,
AFFE=(F(TOUT='OUI',
NOM_CMP=('DX', 'DY', 'DZ'),
VALE=(0.0, 0.0, 0.0),),
F(GROUP_NO=('body', ),
NOM_CMP=('DX', 'DY', 'DZ'),
VALE=(0.0, -5.0, 0.0),),
),
)
evvent = CREA_RESU(OPERATION='AFFE',
TYPE_RESU='EVOL_CHAR',
AFFE=(F(CHAM_GD=vvent0,
INST=(0.0, ),
PRECISION=0.0,
CRITERE='RELATIF',),
_F(CHAM_GD=vvent1,
INST=(1.0, ),
PRECISION=0.0,
CRITERE='RELATIF',),
),
NOM_CHAM='VITE_VENT',)
vent = AFFE_CHAR_MECA(MODELE=model,
EVOL_CHAR=evvent,
INFO=2,
VERI_NORM='OUI',
DOUBLE_LAGRANGE='OUI',)
This field, if I understood correctly, will generate the normal relative velocity that will be used to evaluate the force defined in
fcxv = FORMULE(VALE='VENT**2',
NOM_PARA=('VENT', ),)
and applied in the 'body' (a simple POU_D_E with springs in the extremities).
elprop = AFFE_CARA_ELEM(MODELE=model,
POUTRE=F(SECTION='CERCLE',
FCX=fcxv,
GROUP_MA=('body', ),
CARA=('R', ),
VALE=(0.1, ),
VARI_SECT='CONSTANT',
MODI_METRIQUE='NON',
TUYAU_NSEC=16,
TUYAU_NCOU=3,),
DISCRET=F(REPERE='GLOBAL',
CARA='K_TR_D_N',
GROUP_MA=('a', 'b'),
VALE=(0.0, 1.E5, 0.0, 1.E6, 1.E6, 1.E6),
SYME='OUI',),
INFO=1,)
All of this is finally put to
resnonl0 = DYNA_NON_LINE(MODELE=model,
CHAM_MATER=fdmat,
CARA_ELEM=elprop,
INFO=2,
EXCIT=(_F(CHARGE=fixs,
TYPE_CHARGE='FIXE_CSTE',
MULT_APPUI='NON',),
F(TYPE_CHARGE='SUIV',
CHARGE=vent,
MULT_APPUI='NON',),
),
COMPORTEMENT=F(GROUP_MA=('body', ),
RELATION='ELAS',
DEFORMATION='PETIT_REAC',
RESI_CPLAN_RELA=1.E-06,
RESI_INTE_RELA=1.E-06,
ITER_INTE_MAXI=20,
ITER_INTE_PAS=0,
ITER_CPLAN_MAXI=1,
SYME_MATR_TANG='OUI',
PARM_THETA=1.0,),
INCREMENT=F(LIST_INST=times,
PRECISION=1.E-06,),
SCHEMA_TEMPS=F(SCHEMA='NEWMARK',
BETA=0.25,
GAMMA=0.5,
FORMULATION='DEPLACEMENT',
COEF_MASS_SHIFT=0.0,),
INFO_DBG=('APPARIEMENT', 'MECA_NON_LINE'),
AMOR_RAYL_RIGI='TANGENTE',
METHODE='NEWTON',
NEWTON=F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,
REAC_ITER_ELAS=0,
MATR_RIGI_SYME='NON',),
CONVERGENCE=F(ITER_GLOB_MAXI=10,
ITER_GLOB_ELAS=25,
ARRET='OUI',),
SOLVEUR=F(RENUM='AUTO',
NPREC=8,
ELIM_LAGR='LAGR2',
STOP_SINGULIER='OUI',
TYPE_RESOL='AUTO',
ACCELERATION='AUTO',
LOW_RANK_SEUIL=0.0,
PRETRAITEMENTS='AUTO',
POSTTRAITEMENTS='AUTO',
PCENT_PIVOT=20,
RESI_RELA=-1.0,
GESTION_MEMOIRE='AUTO',
FILTRAGE_MATRICE=-1.0,
MIXER_PRECISION='NON',
MATR_DISTRIBUEE='NON',
METHODE='MUMPS',),
MESURE=F(TABLE='NON',),
ARCHIVAGE=_F(PRECISION=1.E-06,
CRITERE='RELATIF',),
)
and the final message, the fatal one, is this:
!--------------------------------------------------------------------------------------!
<EXCEPTION> <CHARGES3_10> !
Les composantes dans le champ de vent &&NMVGME.FNOE_CALC doivent être exactement DX, DY et DZ. !
!---------------------------------------------------------------------------------------!
What do Aster mean by this? Should I use another names in the variables? I put all in 'D' terms and I can't find or imagine any other error.
Could anyone help me?
I thank you very much for the attention.
Helio.
