Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-11-03 Thread Rossi, Simone
I solved the P2 problem in Ex4 and Ex5. Everything was fine, but for nu = 0.49995, the linear solver was not converging. I did not realize that as I was not printing the output of the linear solver on screen. Changing the preconditioner, I was able to solve the nearly-incompressible P2

Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread David Knezevic
On Thu, Oct 13, 2016 at 12:34 PM, Rossi, Simone wrote: > Dear David, > thanks for your answer. In this case (nu=0.49995), first order elements > typically lock, but second order elements typically do not lock. > In fact many use second order lagrangian elements for nearly >

Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread Rossi, Simone
Dear David, thanks for your answer. In this case (nu=0.49995), first order elements typically lock, but second order elements typically do not lock. In fact many use second order lagrangian elements for nearly incompressible materials. I wanted to use this example just to show that second order

Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread David Knezevic
On Thu, Oct 13, 2016 at 11:57 AM, Rossi, Simone wrote: > Dear all, > I’m playing around with the elasticity tests in the system of equations > examples (more specifically Ex4 and Ex5). > In particular I’m trying to set the poisson ratio to nu = 0.49995. > With this choice

[Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread Rossi, Simone
Dear all, I’m playing around with the elasticity tests in the system of equations examples (more specifically Ex4 and Ex5). In particular I’m trying to set the poisson ratio to nu = 0.49995. With this choice the solution I get using second order lagrangian elements does not make any sense. For