Dear Amir, > Why we must calculate reaction forces at a boundary which has external > force? >
You certainly don't have to. In my example I used this as a demonstration that the forces computed from the internal stresses matched the applied (uniformly distributed) load. If you believe this then you could happily apply the same procedure to any surface undergoing displacement control, knowing that it will produce the desired result. > But we usually calculate reaction forces at BCs with prescribed > displacement which it is better to calculate them from residual vector > (which equal to [-1.0 * internal forces] at these BCs). > As the body is in equilibrium, they are balanced by the stresses generated > by the material = internal forces > Ok, I see now where you're coming from. This might be true, but only for constrained degrees-of-freedom for quasi-static problems. But in his initial post Hamed wasn't specifically referring to reaction forces (unless I misunderstood something), but rather the resultant force on an unspecified boundary surface. Furthermore, for time-dependent problems you may have additional contributions to these components of the residual due to the imposition of constraints on the velocity and acceleration field (if these constraints are non-homogeneous). > Please, did you test the attached code for a plastic material model? > No I haven't. Do you anticipate the result will be any different when using a plasticity model? If so, how? Regards, Jean-Paul -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.