Hello everyone, I am having a problem assembling elements:
Initially, I define three nodes for a beam element (Pou_D_E)
Node1: initial node
Node2: interface node
Node3: final node
I'm generating two submeshes for each macro element with a different circular section. Each macro element have 10 subelements, that is, between nodes 1 and 2 there are 10 elements, and the same apply between nodes 2 and 3.
For some reason, I have a rigid movement (I understand this to be a false assembly between elements), but when I check the mesh, it appears to have continuity:
I am leaving you the .comm and the mesh in .dat in case you can give me a hand.
.comm
DEBUT(LANG='FR')
malla = LIRE_MAILLAGE(FORMAT='MED',
UNITE=20)
model = AFFE_MODELE(AFFE=F(MODELISATION='POU_D_E',
PHENOMENE='MECANIQUE',
TOUT='OUI'),
MAILLAGE=malla)
modeloDiametroMenorYMayor = AFFE_CARA_ELEM(MODELE=model,
POUTRE=(F(CARA=('EP', 'R'),
GROUP_MA=('mallaIzquierdaEdge', ),
SECTION='CERCLE',
VALE=(2.5, 7.5)),
F(CARA=('EP', 'R'),
GROUP_MA=('mallaDerechaEdge', ),
SECTION='CERCLE',
VALE=(7.0, 12.0))))
Acero = DEFI_MATERIAU(ELAS=F(E=200000.0,
NU=0.3))
AceroAplicado = AFFE_MATERIAU(AFFE=F(GROUP_MA=('mallaIzquierdaEdge', 'mallaDerechaEdge'),
MATER=(Acero, )),
MODELE=model)
empotramiento = AFFE_CHAR_MECA(DDL_IMPO=F(DRX=0.0,
DRY=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('nodoEmpotrado', )),
MODELE=model)
carga = AFFE_CHAR_MECA(FORCE_NODALE=F(FZ=-30.0,
GROUP_NO=('nodoCargado', ),
MX=-6200.0),
MODELE=model)
reslin = MECA_STATIQUE(CARA_ELEM=modeloDiametroMenorYMayor,
CHAM_MATER=AceroAplicado,
EXCIT=(F(CHARGE=empotramiento),
F(CHARGE=carga)),
MODELE=model)
ReaccionesYEsfuerzos = CALC_CHAMP(CONTRAINTE=('EFGE_NOEU', ),
FORCE=('FORC_NODA', 'REAC_NODA'),
RESULTAT=reslin)
IMPR_RESU(FORMAT='MED',
RESU=F(RESULTAT=reslin),
UNITE=80)
IMPR_RESU(FORMAT='MED',
RESU=_F(RESULTAT=ReaccionesYEsfuerzos),
UNITE=2)
FIN()
mesh
21 20
1 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00
2 4.20000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
3 4.60000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
4 4.20000000000000e+00 0.00000000000000e+00 0.00000000000000e+00
5 8.40000000000000e+00 0.00000000000000e+00 0.00000000000000e+00
6 1.26000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
7 1.68000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
8 2.10000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
9 2.52000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
10 2.94000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
11 3.36000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
12 3.78000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
13 4.24000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
14 4.28000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
15 4.32000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
16 4.36000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
17 4.40000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
18 4.44000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
19 4.48000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
20 4.52000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
21 4.56000000000000e+01 0.00000000000000e+00 0.00000000000000e+00
1 102 1 4
2 102 4 5
3 102 5 6
4 102 6 7
5 102 7 8
6 102 8 9
7 102 9 10
8 102 10 11
9 102 11 12
10 102 12 2
11 102 2 13
12 102 13 14
13 102 14 15
14 102 15 16
15 102 16 17
16 102 17 18
17 102 18 19
18 102 19 20
19 102 20 21
20 102 21 3
