×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Which Solver can do eigenvalues with Lagrande multipliers

Which Solver can do eigenvalues with Lagrande multipliers

Which Solver can do eigenvalues with Lagrande multipliers

(OP)
I am trying to figure out which solver (in Fortran) will be able to calculate eigenvalues for a matrix that has been augmented with Lagrange multipliers.

Basically I have a FEM model with 20node solids and 4node plates that are coupled with Lagrange multies in a inhouse FE program. I have been using the Lapack DSPGVX subroutine in Fortran 90 with my other models but now with this coupling I am not able because of singular matrix.

Anyone with some words if advice would be of great help...


RE: Which Solver can do eigenvalues with Lagrande multipliers

If i understand well your matrix is singular.
This means some zero value eignevalues
This means rigid body motion in structures or ill conditioning.
Am i right?

Dr. Costas J. Tsaprounis

RE: Which Solver can do eigenvalues with Lagrande multipliers

Nastran has the option to use Lagrange rigid elements (they have stiffness matrices instead of MPS's which are linear rigid elements), if thats any help.

Your singualarity is caused by a degree of freedom without any structural stiffness (or very low), and as such cannot be inverted.
Typical usual examples of this are 1) forgetting to eliminate rigid body modes, 2) connecting dissimilar element types, 3) cross connecting wrong degrees of freedom.

RE: Which Solver can do eigenvalues with Lagrande multipliers

Hi,
due to the mixed-elem type formulation, I'd guess that 40818's point 2) is the responsible. You wouldn't like to solve a singular system, I presume.

Regards

RE: Which Solver can do eigenvalues with Lagrande multipliers

(OP)
Well, Costas... you are right with the rigidbody mode. I am looking at the structures behavior in water so I need the free state.

Ansys, Hypermesh  and other programs are able to come about eigen results even if the matrix is sigular or closed to it.

But my in house program with can not compute a structure that should have the first 6  eigenvalues =0 . I read about some solvers that can handle these situations but the Fortran Lapack DSPGVX subroutine can't.

So HOW do Ansys, Abaqus and ther FE programs deal with this sort of matrix?

Thank you for your time


J.Blair M.Sc.A,ing.

RE: Which Solver can do eigenvalues with Lagrande multipliers

I you run a free normal modes analysis, you should have 6 rigid body modes at 0Hz, and zero strain energy associated with it. It is a standard FE model check to ensuer you dont have any rogue elements flying of into infinity.
I will have to delve back into some theory guides to tell you how the FE solves a singular matrix.

RE: Which Solver can do eigenvalues with Lagrande multipliers

"I will have to delve back into some theory guides to tell you how the FE solves a singular matrix."

This is done by applying an "eigen shift" , that is a factored amount of the mass matrix is added to the stiffness matrix to make the stiffness matrix non-singular, solution is then possible, the mode shapes are not affected by this procedure but all the eigen values are "shifted" to become greater than zero, the solver will then automatically remove this shift from the eigen values before presenting results to the user.

RE: Which Solver can do eigenvalues with Lagrande multipliers

You can use dggev although it won't be as efficient as a solver which is coded specifically for symmetric semi-definite problems like yours. Many of the sparse solvers (MUMPS for example) so you may be able to massage an eigenvalue computation out of them....

RE: Which Solver can do eigenvalues with Lagrande multipliers

I think that Jonhors an pja are correct.
If you want to see how eigne shift is done see
the "Numerical Recipies".
Good luck

Dr. Costas J. Tsaprounis

RE: Which Solver can do eigenvalues with Lagrande multipliers

(OP)
Thank you all, I will be looking in to the eigen shift method.

Kind regards
J Blair

RE: Which Solver can do eigenvalues with Lagrande multipliers

jabla71, did you manage to implement an eigen shift into your model ?

RE: Which Solver can do eigenvalues with Lagrande multipliers

(OP)
Well finally after getting back to this problem.

The DGGEVX routine did the job quite well.

Thank you all for your input

Cheers

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources