[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Getfem-users] Computation of reaction forces associated to a dirichlet

From: David Danan
Subject: [Getfem-users] Computation of reaction forces associated to a dirichlet condition
Date: Fri, 30 Jul 2021 12:41:30 +0200

Dear GetFem community,

i am trying to compute the reaction forces associated to a dirichlet condition, in order to do so, i have followed the advices given in this thread

You can also find enclosed my test case, using the python interface. The steps are the following:
-Solve a simple 2D elasticity problem on a rectangle (dirichlet on the left, neumann on the top)
-Compute the reaction forces arising from the residual without the dirichlet condition, by using the solution in displacement, the tangent matrix and the right hand side term
-Retrieve the value of the multiplier associated to the dirichlet condition

In fact, i expected both quantity to be the same but they were not (2 order of magnitude different, i didn't expect such large values for the multipliers) and i was wondering which one was correct.
If it helps, you can find both quantities in Mult.png (Multipliers is the true multiplier and MultipliersResidual is the residual without the dirichlet condition).

Next, i built a second model similar to the first one but without a dirichlet condition. I added an explicit rhs using the MultipliersResidual  ( Multipliers gave strange results) and tried to solve the problem, you can see the original solution "Displacement" and newsolution "Displacement_test" result in Sol.png.

It seems i am missing something there, i expected the solutions to be the same, not just close enough.
Can you enlighten me?

Best regards,

Description: Text document

Attachment: Sol.PNG
Description: PNG image

Attachment: Mult.PNG
Description: PNG image

reply via email to

[Prev in Thread] Current Thread [Next in Thread]