[SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

konyaro · 2017-02-10 17:21:31

Hello,
I'm facing the following problem:

A beam is sliding on a fixed part without friction. The reaction forces should be the same on the moving nodes as on the fixed ones. But I get a difference:

- Y reaction force on the fixed nodes: 46.3 N
- Y reaction force on the moving nodes: 44.3 N
- normal contact force RN: 44.3 N

Where do these 2N on the fixed nodes come from? I guess I'm missing something very simple but I cannot find what.

Thank you for your help,

Konyaro

Last edited by konyaro (2018-11-25 06:54:14)

Attachments:
ReactionForces.png, Size: 90.36 KiB, Downloads: 449

konyaro · 2017-02-10 17:22:37

The comm and med files.

Attachments:
BeamWithContact.zip, Size: 74.2 KiB, Downloads: 426

konyaro · 2017-02-20 17:24:14

Compared with another software, the reaction force of the fixed nodes is similar (cf. attached picture).

Nevertheless I still can't figure out why the reaction force of the moving nodes is not the same in the CA simulation...

Last edited by konyaro (2017-02-21 16:18:11)

Attachments:
Force_comparaison.png, Size: 131.36 KiB, Downloads: 462

bbd · 2017-02-21 15:46:40

there is something I don't understand.
Your second calculation seems relevant.
Then the conclusion is the code_aster simulation is not correct, isn't it ?

konyaro · 2017-02-21 16:17:48

Exactly, I expect to get the same reaction forces on both sides.

bbd · 2017-02-21 16:21:38

there are warnings about the time step (=increment ?) which is not adapted
during the calculation.
The results do not seem correct as the X component oscillates from one step to another, like numerical stick-slip

Last edited by bbd (2017-02-21 16:22:47)

konyaro · 2017-02-21 16:54:20

Thank you for your advice. I tried with a very small time step but the results are still the same. The X force is strange indeed, it should be around 0 as there is no friction and almost no deformation of the contact surface.

bbd · 2017-02-21 17:33:17

newbie question. Isn't there a pb of orientation of edges ?

konyaro · 2017-02-22 08:20:32

maybe, but as I already use the ORIE_PEAU_2D I can't see what else I could change...

bbd · 2017-02-22 08:32:11

the mesh.Refine locally

konyaro · 2017-02-22 15:52:04

I don't think it's due to the mesh. Even with a crazy refined quadratic mesh the reaction forces do not change.

Attachments:
MeshAffineQuad.PNG, Size: 83.84 KiB, Downloads: 347

bbd · 2017-02-22 17:50:17

exactly same values ?

konyaro · 2017-02-23 07:40:46

Not exactly, but close, as depicted on the attached plot.

Attachments:
Forces.png, Size: 145.89 KiB, Downloads: 307

bbd · 2017-02-23 13:55:20

I removed line 103 :
# DEFORMATION='GROT_GDEP',

and it works. I mean the balance of forces is correct.

then I read http://www.code-aster.org/forum2/viewtopic.php?id=18029 and it was written that for non linear calculations (either material and geometry) you must use 3D shells.
Hope it helps.

Last edited by bbd (2017-02-23 19:54:46)

konyaro · 2017-02-24 23:08:45

Thank you for your help bbd,
You are right, the reaction forces are not correct with C_PLAN elements and GROT_GDEP deformations ! I made a simple comparison on a beam in flexure:

1) C_PLAN and PETIT deformations: forces OK
2) C_PLAN and GROT_GDEP deformations: forces KO
3) 3D and PETIT deformations: forces OK
4) 3D and GROT_GDEP deformations: forces OK

That is not a good news as most of my calculations are C_PLAN/GROT_GDEP... I found nothing in the documentation that mentions this limitation.

Last edited by konyaro (2017-02-24 23:12:33)

Attachments:
Image_2D_3D_Forces.png, Size: 240.05 KiB, Downloads: 333

bbd · 2017-02-25 08:03:57

I believe it is written in U2.02.01 section 2.3.2
but I agree the way it is presented is not clear.
The 'available' behaviors and laws are listed in a table but it's worth telling a few words about the risks of using this combination of asumptions (elements and calculation type) as nothing prevents us using such combinations.

konyaro · 2017-02-25 15:02:59

According to U4.51.11 section 4.5.2 it should work...

Attachments:
U4_51_11_P59.png, Size: 39.99 KiB, Downloads: 327

bbd · 2017-03-07 15:57:17

il faudrait peut-être signaler le bug ?

konyaro · 2017-03-08 06:32:27

Oui, je vais le faire.

konyaro · 2018-11-25 06:53:55

I never closed this topic but is is solved.
There is a mandatory workaround since 2017: the following error occurs when using GROT_GDEP with C_PLAN and elastic material :

! On ne peut pas utiliser les contraintes planes en élasticité avec une !
! déformation de type GROT_GDEP. !
! Conseil: utilisez une loi élastoplastique sur laquelle vous aurez mis une !
! limite d'élasticité très grande. !

The trick is to define a non linear material with a very high yield strength, for instance :
mat01 = DEFI_MATERIAU(ECRO_LINE=_F(D_SIGM_EPSI=1.0,
SY=1e+15),
ELAS=_F(E=200000.0,
NU=0.3))

and to use the following relation in COMPORTEMENT :
RELATION='VMIS_ISOT_LINE'

Last edited by konyaro (2018-11-25 12:49:53)

#1 2017-02-10 17:21:31

[SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#2 2017-02-10 17:22:37

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#3 2017-02-20 17:24:14

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#4 2017-02-21 15:46:40

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#5 2017-02-21 16:17:48

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#6 2017-02-21 16:21:38

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#7 2017-02-21 16:54:20

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#8 2017-02-21 17:33:17

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#9 2017-02-22 08:20:32

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#10 2017-02-22 08:32:11

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#11 2017-02-22 15:52:04

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#12 2017-02-22 17:50:17

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#13 2017-02-23 07:40:46

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#14 2017-02-23 13:55:20

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#15 2017-02-24 23:08:45

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#16 2017-02-25 08:03:57

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#17 2017-02-25 15:02:59

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#18 2017-03-07 15:57:17

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#19 2017-03-08 06:32:27

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

#20 2018-11-25 06:53:55

Re: [SOLVED] Reaction forces not equal with C_PLAN/GROT_GDEP

Board footer