User login

Navigation

You are here

Eigenvalues of the incompressibile elasticity problem

WaiChing Sun's picture

Dear all,

    I have a question about the eigenvalue of the incompressible elasticity problem. Assume that we have a incompressible elasticity continuum body whose constitutive response holds normality. Obviously, we can treat the incompressiblity limit as a Lagrangian constraint and express the governing equation via the Lagrange-Multiplier method. If mixed finite element method is used, one can use Courant-Fischer theorem to show that the tangential stiffness matrix of  constrained variational form would have its maximum eigenvalue larger than its unconstrained counterpart and its minimum eigenvalue smaller than its uncsontrained counterpart.

On the other hand, if penalty method is used, then Wely's monotonicity theorem shows that the smallest eigenvale of the stiffness matrix K_ij + alpha I_i _j would be always larger than that of K_ij.

Now, suppose we choose not to discretize the problem. Instead, assume that the elasticity problem can be expressed with a compact, self-adjointed linear operators A on a Hibert space such that

A(u) = f

where A(u) = div (C: epsilon(u))  

Would the incompressible constraint impose a simiar effect on the linear operator A? Is there any published results that show how the eigenvalue of the operator A changes to the constraint which is somehow "equvalient" to the monotonicity therorem/Courant-Fischer theorem applied to the discrete Hessian matrix? Any comment/suggestion is appreciated. Thank you.

Regards,

WaiChing Sun 

 

 

Comments

Not sure exactly what you're looking for.  I assume you've already looked at you  Auchmuty, G., Variational principles for self-adjoint elliptic eigenproblems and references in Nonsmooth/nonconvex mechanics by Gao, Ogden and Stavroulakis?

-- Biswajit

WaiChing Sun's picture

Hi Biswajit,

      Thank you for pointing out the reference by Auchmuty. Actually, I am interested in how the smallest eigenvalue will shift if a certain Lagrangian constraint is imposed to an unconstrained system. Since solutions of finite element are in finite dimensional space, it is not hard to determine how  the smallest eigenvalue increase/decrease based on the property of the Hessian matrix. However, I am not sure how the eigenvalue of the unconstrained system would changes when a f(u)=0 constraint is imposed (assme that C does not depend on u for now and no discrietization take place). My gut feeling is that the smallest eigenvalue would shift in a similiar fashion regardless of whether a discrietization take place or not, but I am not sure how to prove/disprove such thing.  I have not yet finished reading the paper you kindly pointed out to me but it seems quite useful.  

Best Regards,

WaiChing Sun 

 

Hi Dear All

I am a beginner in using abaqus, I am trying to model buckling of a plate but I recieve an error "TOO MANY ITERATIONS NEEDED TO SOLVE THE EIGENVALUE PROBLEM" , does any one has any comments? I really need help?

Saba

Subscribe to Comments for "Eigenvalues of the incompressibile elasticity problem"

Recent comments

More comments

Syndicate

Subscribe to Syndicate