Modal Participation Factor

A forum dedicated to feature requests and the future direction of OpenSees, i.e. what would you like, what do you need, what should we explore

Moderators: silvia, selimgunay, Moderators

Post Reply
ja_abell
Posts: 43
Joined: Tue Nov 03, 2009 1:22 pm
Location: Universidad de los Andes
Contact:

Modal Participation Factor

Post by ja_abell » Tue Dec 15, 2009 10:12 am

Hi OpenSEES developers,

OpenSEES really needs a modal participation factor calculation procedure. It's something very common in dynamic analysis of any type. I really miss it.

It shouldn't be to hard to program, maybe I could do it if I understood how OpenSEES works. But just looking at the codes doesn't help much.

If someone gives me a hint start I could try, and I'd love to collaborate.

Thank you.
Jose Abell

silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Tue Dec 15, 2009 11:18 am

can't this be done at the Tcl level?
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

ja_abell
Posts: 43
Joined: Tue Nov 03, 2009 1:22 pm
Location: Universidad de los Andes
Contact:

Post by ja_abell » Wed Dec 16, 2009 4:06 am

Hi Silvia,

I'm sure it could be done. All I need is the mass matrix, and that I don't really know how to extract. I've seen some threads in the forum and no one has actually answered how to do it. FMK pasted some code in one thread, but I couldn't make anything of it.

You think this is easier?

Thank you,
Jose Abell

silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Wed Dec 16, 2009 5:43 am

i haven't spent enough time researching it.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

pbeng
Posts: 19
Joined: Tue Jun 17, 2008 2:13 pm
Location: West Falls, NY

Post by pbeng » Mon Dec 28, 2009 1:29 pm

For beam problems, lump 1/2 the element mass at each of its nodes. Assuming the eigenvectors are mass normalized, multiply the eigenvector displacement for the direction of interest times the nodal mass. Sum all these products up to yield the participation factor.

For plates and solids a good FEA text will have the equations to calculate the consistent mass matrix for an element, from which the lumped masses at each node can be calculated (there are different ways to do this). Once the lumped mass matrix is found for the element nodes, repeat the procedure above.

This assumes that the influence vector is 1 for all nodes in the direction of interest, and zero elsewhere.

Agreed, it'd be nice to be able to query the mass at a node rather than having to compute it.

PB

Post Reply