Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


Options
Go to last post Go to first unread
Offline oscampo  
#1 Posted : 14 September 2023 02:41:44(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 247
Man
Colombia
Location: Cali, Colombia

Was thanked: 87 time(s) in 66 post(s)
Here is a problem from a Matrix Analysis of Structures. I'm trying to solve by using the file provided here by Alvaro Gavilan, but for this exercise, the calculated stiffness matrix is singular.
Any ideas on how to solve and found the correct answers?

a Frame.png

A Frame with internal releases.sm (249kb) downloaded 17 time(s).

Edited by user 14 September 2023 20:03:11(UTC)  | Reason: Not specified

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline Razonar  
#2 Posted : 14 September 2023 06:59:28(UTC)
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,356
Uruguay

Was thanked: 815 time(s) in 516 post(s)
Hola Oscar. Hope this helps: Please check example 7.1, page 353, from this book: Aslam_Kassimali_Matrix_Analysis_of_Structure. There are also an SI version of that book, but I don't have it.

Best regards.
Alvaro.
Offline oscampo  
#3 Posted : 14 September 2023 12:19:45(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 247
Man
Colombia
Location: Cali, Colombia

Was thanked: 87 time(s) in 66 post(s)
Originally Posted by: Razonar Go to Quoted Post
Hola Oscar. Hope this helps: Please check example 7.1, page 353, from this book: Aslam_Kassimali_Matrix_Analysis_of_Structure. There are also an SI version of that book, but I don't have it.

Best regards.
Alvaro.


Gracias Alvaro,

I have precisely extracted this example from that book (SI Version) and I am trying to solve it using the file provided by Alvaro Gavilán, which he has successfully tested with other examples, but in this one, even though the loading conditions and the model for each bar are correct, it is not possible to calculate the system...
thanks 1 user thanked oscampo for this useful post.
on 15/09/2023(UTC)
Offline oscampo  
#4 Posted : 14 September 2023 20:02:48(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 247
Man
Colombia
Location: Cali, Colombia

Was thanked: 87 time(s) in 66 post(s)
Solved...
Frame solved.PNG

Custom FEM implementation.
I hope it would be useful for others:

A Frame with internal releases-SOLVED.sm (345kb) downloaded 24 time(s).

Best,
Oscar
thanks 2 users thanked oscampo for this useful post.
on 15/09/2023(UTC),  on 15/09/2023(UTC)
Offline Alvaro Gavilán  
#5 Posted : 16 September 2023 01:05:30(UTC)
Alvaro Gavilán


Rank: Advanced Member

Groups: Registered
Joined: 16/04/2020(UTC)
Posts: 65
Man
Paraguay
Location: France

Was thanked: 25 time(s) in 16 post(s)
Hola oscampo!

When singular matrices appear its probably due to a problem with hinged supports or connections. I verified the FE matrices and they are correct. I just changed some loads (that are corrected in your solved file) and then I changed two things:

1) Supports are always "rigid" when considering the member type. If we have a simple support for a frame, the beam that arrives to that support is a "rigid-rigid" member, we will just free the rotation in the constraints matrix (phi_z = 0).

2) When two members of a frame are connected with a hinge, one will be "rigid-hinge" and the other "rigid-rigid" (not "hinge-rigid"!). In these situations you have two rotations for the two beams that arrive at the hinge. To know the rotation corresponding to a certain beam, this has to be the "rigid-rigid" one.

A way of seeing this criteria is: we always prescribe rotations as unknowns, someone has to "withstand" those rotations (you can later "free" them with the constraints matrix), but if none of the elements can withstand that rotation, you will have a row (and column) full of zeros (zero stiffness) for that displacement in your K matrix.

Hope it is more clear!

Saludos,

Alvaro
A Frame with internal releases - v2.sm (257kb) downloaded 39 time(s).
Users browsing this topic
Guest
Forum Jump  
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.