0
$\begingroup$

I have a large scale system with the following equations of motions

$M\ddot z(t)+K z(t)=F(t)$

$M$ is the mass matrix

$K$ is the stiffness matrix

And I want to solve the equation for displacement $z(t)$, the eigenvalues and the eigenvectors with load dependent Ritz Analysis ( the algorithm is described succinctly in here).

enter image description here

I know that I can use APRACK to for large scale eigenvalue problems, but what about load dependent Ritz Analysis? Is there already a solver ( commercial, open source, all doesn't matter) out there that does this kind of thing?

Just to be clear, I am looking for a programming library that I can use, not a commercial software that I can't call programmatically .

$\endgroup$
  • $\begingroup$ Sorry for answering without any prior knowledge on the subject, but to me the algorithm described in your scan looks like it could be implemented very easily using any library that can do the two things: 1) solution of a linear system (in particular, steps (3) and (5)) and 2) matrix-vector products (steps (4), (6) and (7)). This library could be, e.g., SciPy, cholmod, SuperLU, etc. depending on the properties of $K$. $\endgroup$ – knl Jan 17 '17 at 11:31
  • $\begingroup$ @knl, would you like to post your comment as an answer? Because it is $\endgroup$ – Graviton Feb 24 '17 at 6:35
2
$\begingroup$

The algorithm described in your scan looks like it could be implemented using a library that can do the two things: 1) solution of a linear system (in particular, steps (3) and (5)) and 2) matrix-vector products (steps (4), (6) and (7)). This library could be, e.g., SciPy, cholmod, SuperLU, etc. depending on the properties of the matrix $K$.

$\endgroup$
0
$\begingroup$

The commercial program, SAP 2000, http://www.comp-engineering.com/SAPManE.htm has an option for dynamic analysis with Ritz vectors. Take a look at page 74 in the Basic Analysis Reference Manual.

$\endgroup$
  • $\begingroup$ This is not a solution that I can integrate into my programming development, see the updated question $\endgroup$ – Graviton Jun 17 '16 at 12:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.