 Scaling the problem

Hai, Can we model this problem in Fenics?

This looks like the energy functional for plane-strain elasticity. There is a good FEniCS-based tutorial on that problem here, by @bleyerj. That tutorial doesn’t contain the traction term, but it could easily be added to the right-hand side form l as dot(t,u_)*ds, where t is the traction vector and ds integrates over the exterior boundary of the domain.

2 Likes

See the thickness is constant and taken out of the equation and we can do an area integral. can we model this in the fenics? How to model this? any resources

If the thickness is a constant (independent of position in the plane, as needed to move it outside of the integral), then its value won’t affect the displacement solution. The tutorial I linked chooses h^e = 1.

So let’s say if we apply this idea to a plate in the tension of thickness 10 mm, if I multiply with 10 mm for including thickness effect. Will the solution gets affected or remain the same?

Think about it this way: If you find a displacement field \mathbf{u} that minimizes the energy \Pi(\mathbf{u}), then it will also be a minimum of C\Pi(\mathbf{u}), for any constant C. Since every term in your energy functional is scaled by h^e, then you can set it to 1 without loss of generality.

If the superscript e is meant to indicate that the thickness changes from element to element, then it will affect results. The easiest way to implement this would be to represent thickness as a Function in the space FunctionSpace(mesh,"DG",0) (i.e., the space of constants on elements), and simply keep it inside of the dx integral (over the whole domain) used to define variational forms in FEniCS.

1 Like

Thanks for the answer, I have changed my view about the problem. Can you elaborate or show some way to implement this " The easiest way to implement this would be to represent thickness as a Function in the space FunctionSpace(mesh,"DG",0)"

Also, how can one implement area integral, like in the problem state above?