Eigen Command: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(11 intermediate revisions by 5 users not shown)
Line 4: Line 4:


{|  
{|  
| style="background:yellow; color:black; width:800px" | '''analyze <$type> <$solver> $numEigenvalues'''
| style="background:lightgreen; color:black; width:800px" | '''eigen  <$solver> $numEigenvalues'''
|}
|}


Line 12: Line 12:
|  style="width:150px" | '''$numEigenvalues''' || number of eigenvalues required
|  style="width:150px" | '''$numEigenvalues''' || number of eigenvalues required
|-
|-
|  '''$type''' || optional string detailing type of eigen analysis, -standard of -generalized (default generalized)
| '''$solver'''|| optional string detailing type of solver: -genBandArpack, -fullGenLapack (default: -genBandArpack)
|-
| '''$solver'''|| optional string detailing type of solver, -UmfPack, -genBandArpack, -symmBandLapack, -SuperLU (default genBandArpack)
|}
|}


Line 25: Line 23:
NOTES:
NOTES:


#The eigenvectors are stored at the nodes and can be printed out using a Node Recorder, the nodeEigenvector command,  or the Print command.
# The eigenvectors are stored at the nodes and can be printed out using a Node Recorder, the nodeEigenvector command,  or the Print command.
# The default eigensolver is able to solve only for N-1 eigenvalues, where N is the number of inertial DOFs. When running into this limitation the -fullGenLapack solver can be used instead of the default Arpack solver.
# The -fullGenLapack option is VERY SLOW for moderate to large models




THEORY:
THEORY:


Standard Problem:
:<math> \left (K - \lambda I \right ) \Phi = 0 </math>
Generalized Problem:


:<math> \left (K - \lambda M \right ) \Phi = 0 </math>
:<math> \left (K - \lambda M \right ) \Phi = 0 </math>
Line 43: Line 38:
EXAMPLE:
EXAMPLE:


set eigenvalues [eigenvalue 10];  
set eigenvalues [eigen 10];  





Latest revision as of 19:38, 19 August 2022




This command is used to perform the analysis.

eigen <$solver> $numEigenvalues

$numEigenvalues number of eigenvalues required
$solver optional string detailing type of solver: -genBandArpack, -fullGenLapack (default: -genBandArpack)


RETURNS:

a tcl string containg eigenvalues.


NOTES:

  1. The eigenvectors are stored at the nodes and can be printed out using a Node Recorder, the nodeEigenvector command, or the Print command.
  2. The default eigensolver is able to solve only for N-1 eigenvalues, where N is the number of inertial DOFs. When running into this limitation the -fullGenLapack solver can be used instead of the default Arpack solver.
  3. The -fullGenLapack option is VERY SLOW for moderate to large models


THEORY:


<math> \left (K - \lambda M \right ) \Phi = 0 </math>



EXAMPLE:

set eigenvalues [eigen 10];



Code Developed by: fmk