Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Eigen analysis of deflected beam? #101

Open
RibeiroAndre opened this issue Mar 7, 2023 · 4 comments
Open

Eigen analysis of deflected beam? #101

RibeiroAndre opened this issue Mar 7, 2023 · 4 comments

Comments

@RibeiroAndre
Copy link

Hi, Taylor.
Can you advise on how to calculate the modal properties of a wing during a time_domain_analysis?
Say a beam is vibrating and I'd like to do eigenvalue_analysis to see its eigenmodes and frequencies at different deflections. I'm not sure this is possible at the moment.
Thanks a lot!

@taylormcd
Copy link
Collaborator

taylormcd commented Mar 8, 2023

A linear model of a nonlinear system is a valid qualitative description of the nonlinear system (in a limited area) when the system is linearized about a steady state operating point. That's why the default is to perform eigenvalue analyses about steady state operating conditions. There is a steady keyword argument which you can use to skip the steady solve (the states in system will be used instead).

@RibeiroAndre
Copy link
Author

Hi, Taylor.
I tried doing this:

system,history,converged = time_domain_analysis(assembly, t;
    prescribed_conditions = prescribed_conditions,
    point_masses = structure.pointmass,
    structural_damping = false);
system, λ, V, converged = GXBeam.eigenvalue_analysis!(system, assembly, prescribed_conditions = prescribed_conditions, nev = 10, point_masses=pointmass,structural_damping=false,steady=false);
state = AssemblyState(system, assembly; prescribed_conditions);
eigenstates = [ AssemblyState(system, assembly, V[:,k]; prescribed_conditions)
    for k = 1:10 ];
write_vtk("eigenmodes", assembly, state, λ[1], eigenstates[1]);

But the eigenmodes are getting are identical to the undeformed beam and when I look at the VTK it is also undeformed. I'm deforming the beam quite aggressively, so I'd expect some effect on the eigenfrequencies and to see the eigenmodes on the deformed shape. Am I doing something wrong?
Thanks a lot!

@taylormcd
Copy link
Collaborator

Looks like you might need to set linearization_state=system.x

@RibeiroAndre
Copy link
Author

RibeiroAndre commented Mar 14, 2023

In the eigenvalue_analysis! calculation? It didn't seem to change anything. Also, after the eigenvalue_analysis, system.x is full of zeros.
Thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants