Error in solve command while analysing a cantilever beam under prescribed displacement

[dokken]

There are many posts regarding this on the forum, including

Issue in calculating reaction force by integration of stress over a subdomain

Hello everyone, I am trying to calculate the reaction force through integration of stresses over a subdomain as done on the following link: https://comet-fenics.readthedocs.io/en/latest/demo/tips_and_tricks/computing_reactions.html I am attaching a minimal working code for 1D bar in tension. Left end of this bar is fixed and a displacement is applied on the right end. The program is providing the zero reaction force. The same method worked for me in dolfin but in dolfinx it is giving zero rea…

and

Computation of reaction force

Because all external forces come from BCs, and the solution is in static equilibrium, you expect the net reaction force to all BCs to be zero. If you want to get the reaction force from one BC, you need to apply the others to the residual. For example, to get the reaction to the force on the right end, you could do something like the following: for n in range(num_steps): t += dt solve(F == 0, w, bc) (u,p)=w.split(True) p_n.assign(p) f_int = assemble(inn…