Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

COMSOL, structural 3D, shell elements, livelink for Matlab, mphmatrix, constrained mass and stiffness matrices

Please login with a confirmed email address before reporting spam

Hello,

I wonder how to obtain the mass and stiffness matrices of the problem I am working with. The structural problem is modelled using 3D shell elements, with a total of 624 nodes of which 12 nodes are constrained, i.e. with 612*6=3672 degrees of freedom. My intention is to obtain the mass and stiffness matrices of the constrained structure, i.e. the matrices with size (3672x3672) and the description of the degrees of freedom. The livelink for matlab function mphmatrix gives me either:

Unconstrained mass and stiffness matrices E and K with size (14832x14832)

Constrained mass and stiffness matrices Ec and Kc with size (12240x12240)

I don't understand why the constrained matrices have more degrees of freedom than: number of degrees of freedom per node x number of nodes

Kind regards,

Lucas



1 Reply Last Post 2020年10月10日 GMT-4 19:54
COMSOL Moderator

Hello Lucas Silveira

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.


Please login with a confirmed email address before reporting spam

Posted: 4 years ago 2020年10月10日 GMT-4 19:54

Hej Lucas,

After checking with Livelink the nodes coordinates the number for different nodes is 2472. If you multiply with the number of DoFs you get your full matrix size. An eliminated matrix has removed the DoFs related to boundary conditions. A way to verify the matrices is eig(A,B) and remembering that the eigenfrequencies that you ask from Comsol should be calculated from the mass and stiffness matrix.

https://www.comsol.com/multiphysics/eigenfrequency-analysis

Hej då, Theodore

Hej Lucas, After checking with Livelink the nodes coordinates the number for different nodes is 2472. If you multiply with the number of DoFs you get your full matrix size. An eliminated matrix has removed the DoFs related to boundary conditions. A way to verify the matrices is eig(A,B) and remembering that the eigenfrequencies that you ask from Comsol should be calculated from the mass and stiffness matrix. https://www.comsol.com/multiphysics/eigenfrequency-analysis Hej då, Theodore

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.