Hello,

I am running a DYNA_NON_LINE calculation with a model consisting of:

1d beam (mesh SEG2) -> mesh of 21 elements SEG2
2d section (mesh TRIA3), affecting the beam with MACR_CARA_POUTRE -> mesh of 400 elements TRIA3
3d brackets (supports), that hold the beam -> mesh of 42000 elements TETRA4

The objective is to calculate deformation and stress at the brackets when an object is dropped from 2 meters and left hanging from the beam. I would also know which is the force on the beam caused by this drop.

To do so, I have split the case in 2 stages, first one to setup the mesh, model, material, etc. I apply gravity (to have a load, because if not the calculation fails) and a speed with CREA_CAMP (equivalent to the 2 m fall, according to formula v=sqr(2gH) ):

speed0 = CREA_CHAMP(identifier='21:1',
AFFE=_F(GROUP_NO=('load', ),
NOM_CMP=('DZ', ),
VALE=(-6.261, )),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_DEPL_R')

At stage 2, I take the result of stage 1 as initial state:

dynanl1 = DYNA_NON_LINE(identifier='7:2',
...
ETAT_INIT=_F(EVOL_NOLI=statnl0,
VITE=speed0),
...
NEWTON=_F(MATR_RIGI_SYME='OUI'),
SCHEMA_TEMPS=_F(FORMULATION='DEPLACEMENT',
SCHEMA='HHT'))

But I have some doubts:

Is correct to impose the speed as CREA_CHAMP at the end of stage 1 and run the calculation at stage 2?
Is "HHT - deplacement" the correct schema? (It took it from jeanpierreaubry's book)
In my case, shouldn't it be "HHT - vitesse"? I tried but it says that is not compatible...
What is the difference between schemas Newmark and HHT ? I read U4.53.01, chapter 3.19 but it is not clear about cases of use...
The calculation is running for 3 days now, which leads me to think there is something wrong with the setup...
Any help or idea about the physical phenonmenon will be welcome.

Thanks.

P.S: Full code:

#stage 1: setup and application of drop
DEBUT(LANG='EN')
mesh1 = LIRE_MAILLAGE(identifier='0:1',
UNITE=7)
mesh2 = LIRE_MAILLAGE(identifier='1:1',
UNITE=8)
mesh3 = LIRE_MAILLAGE(identifier='2:1',
UNITE=2)
msh_complete = ASSE_MAILLAGE(identifier='5:1',
MAILLAGE_1=mesh1,
MAILLAGE_2=mesh2,
OPERATION='SUPERPOSE')
model = AFFE_MODELE(identifier='6:1',
AFFE=(_F(GROUP_MA=('beam', ),
MODELISATION='POU_D_T_GD',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('3dA', '3dB'),
MODELISATION='3D',
PHENOMENE='MECANIQUE')),
MAILLAGE=msh_complete)
table = MACR_CARA_POUTRE(identifier='7:1',
GROUP_MA_BORD=('border_ext', ),
GROUP_MA_INTE=('border_int', ),
GROUP_NO=('nodeX', ),
MAILLAGE=mesh3,
NOM='profile',
TABLE_CARA='OUI')
elemprop = AFFE_CARA_ELEM(identifier='8:1',
MODELE=model,
ORIENTATION=_F(CARA='VECT_Y',
GROUP_MA=('beam', ),
VALE=(0.0, -1.0, 0.0)),
POUTRE=_F(GROUP_MA=('beam', ),
NOM_SEC='profile',
SECTION='GENERALE',
TABLE_CARA=table))
alum = DEFI_MATERIAU(identifier='9:1',
ELAS=_F(E=67000000000.0,
NU=0.33,
RHO=2700.0))
steel = DEFI_MATERIAU(identifier='10:1',
ECRO_LINE=_F(D_SIGM_EPSI=2100000000.0,
SY=275000000.0),
ELAS=_F(E=210000000000.0,
NU=0.3,
RHO=7500.0))
fieldmat = AFFE_MATERIAU(identifier='11:1',
AFFE=(_F(GROUP_MA=('beam', ),
MATER=(alum, )),
_F(GROUP_MA=('3dA', '3dB'),
MATER=(steel, ))),
MODELE=model)
fload = DEFI_FONCTION(identifier='12:1',
NOM_PARA='INST',
VALE=(0.0, 0.0, 1.0, 1.0))
listr0 = DEFI_LIST_REEL(identifier='13:1',
DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1,
NOMBRE=10))
times0 = DEFI_LIST_INST(identifier='14:1',
ADAPTATION=_F(VALE_I=8),
DEFI_LIST=_F(LIST_INST=listr0),
ECHEC=_F(SUBD_PAS=5),
METHODE='AUTO')
BC = AFFE_CHAR_MECA(identifier='15:1',
DDL_IMPO=(_F(DRY=0.0,
DRZ=0.0,
DX=0.0,
DZ=0.0,
GROUP_NO=('supportC', )),
_F(DRY=0.0,
DRZ=0.0,
DX=0.0,
DZ=0.0,
GROUP_NO=('supportD', )),
_F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('fixA', 'fixB'))),
MODELE=model)
binding = AFFE_CHAR_MECA(identifier='17:1',
LIAISON_UNIF=(_F(DDL=('DX', 'DY', 'DZ'),
GROUP_NO=('bracketA1', 'supportA1')),
_F(DDL=('DX', 'DY', 'DZ'),
GROUP_NO=('bracketA2', 'supportA2')),
_F(DDL=('DX', 'DY', 'DZ'),
GROUP_NO=('bracketB1', 'supportB1')),
_F(DDL=('DX', 'DY', 'DZ'),
GROUP_NO=('bracketB2', 'supportB2'))),
MODELE=model)
gravity = AFFE_CHAR_MECA(identifier='18:1',
MODELE=model,
PESANTEUR=_F(DIRECTION=(0.0, 0.0, -1.0),
GRAVITE=9.8))
statnl0 = STAT_NON_LINE(identifier='19:1',
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
COMPORTEMENT=(_F(DEFORMATION='GROT_GDEP',
GROUP_MA=('beam', ),
RELATION='ELAS_POUTRE_GR'),
_F(DEFORMATION='SIMO_MIEHE',
GROUP_MA=('3dA', '3dB'),
RELATION='VMIS_ISOT_LINE')),
EXCIT=(_F(CHARGE=BC),
_F(CHARGE=binding),
_F(CHARGE=gravity)),
INCREMENT=_F(LIST_INST=times0),
MODELE=model,
NEWTON=_F(MATR_RIGI_SYME='OUI'))
speed0 = CREA_CHAMP(identifier='21:1',
AFFE=_F(GROUP_NO=('load', ),
NOM_CMP=('DZ', ),
VALE=(-6.261, )),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_DEPL_R')
IMPR_RESU(identifier='22:1',
FORMAT='RESULTAT',
RESU=_F(CARA_ELEM=elemprop,
GROUP_MA=('beam', ),
INST=(1.0, ),
NOM_CHAM=('DEPL', 'SIEF_ELGA'),
RESULTAT=statnl0),
UNITE=3)
FIN()

#stage 2: dynamic calculation
POURSUITE(LANG='EN')
listr2 = DEFI_LIST_REEL(identifier='3:2',
DEBUT=2.0,
INTERVALLE=_F(JUSQU_A=3.0,
NOMBRE=10))
times2 = DEFI_LIST_INST(identifier='4:2',
ADAPTATION=_F(VALE_I=8),
DEFI_LIST=_F(LIST_INST=listr2),
ECHEC=_F(SUBD_PAS=5),
METHODE='AUTO')
dynanl1 = DYNA_NON_LINE(identifier='7:2',
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
COMPORTEMENT=(_F(DEFORMATION='GROT_GDEP',
GROUP_MA=('beam', ),
RELATION='ELAS_POUTRE_GR'),
_F(DEFORMATION='SIMO_MIEHE',
GROUP_MA=('3dA', '3dB'),
RELATION='VMIS_ISOT_LINE')),
ETAT_INIT=_F(EVOL_NOLI=statnl0,
VITE=speed0),
EXCIT=(_F(CHARGE=BC),
_F(CHARGE=binding)),
INCREMENT=_F(LIST_INST=times2),
MODELE=model,
NEWTON=_F(MATR_RIGI_SYME='OUI'),
SCHEMA_TEMPS=_F(FORMULATION='DEPLACEMENT',
SCHEMA='HHT'))
dynanl1 = CALC_CHAMP(identifier='8:2',
reuse=dynanl1,
CRITERES=('SIEQ_ELGA', ),
INST=(3.0, ),
RESULTAT=dynanl1)
IMPR_RESU(identifier='9:2',
FORMAT='RESULTAT',
RESU=(_F(INST=(3.0, ),
NOM_CHAM=('DEPL', 'SIEF_ELGA'),
RESULTAT=dynanl1),
_F(GROUP_MA=('3dA', ),
INST=(3.0, ),
NOM_CHAM=('SIEQ_ELGA', ),
NOM_CMP=('VMIS', ),
RESULTAT=dynanl1,
VALE_MAX='OUI',
VALE_MIN='OUI')),
UNITE=5)
IMPR_RESU(identifier='10:2',
FORMAT='MED',
RESU=_F(CARA_ELEM=elemprop,
LIST_INST=listr2,
RESULTAT=dynanl1),
UNITE=80)
FIN()