Plane Trusses Finite Elements 2 - Assembled Matrix

AI Thread Summary
The discussion centers on the assembly of a global stiffness matrix for plane trusses using finite element analysis. Participants express confusion regarding the naming convention and the assembly process, particularly how global and element numbers are represented in the matrix. The matrix provided illustrates the contributions from different elements, emphasizing the symmetry and the degrees of freedom at each node. A method for creating an Element Freedom Table (EFT) is shared, demonstrating how to derive global node numbers based on degrees of freedom. Overall, the conversation seeks clarity on the assembly technique and its underlying principles.
bugatti79
Messages
786
Reaction score
4
Folks,

I am having difficulty understanding how this global matrix is assembled with the naming convention used as shown in attached.

The numbers in the corners such as 1(1,2) etc in figure 4.6.3 (b) denote the global and element numbers respectively.

Can anyone shed light on how this assembled...? Thanks
 

Attachments

  • Assembled Matrix.jpg
    Assembled Matrix.jpg
    11.9 KB · Views: 501
Engineering news on Phys.org
bugatti79 said:
Folks,

I am having difficulty understanding how this global matrix is assembled with the naming convention used as shown in attached.

The numbers in the corners such as 1(1,2) etc in figure 4.6.3 (b) denote the global and element numbers respectively.

Can anyone shed light on how this assembled...? Thanks

Maybe if I write out the matrix attached for easier read...the blanks indicate symmetry.

##\begin{bmatrix}
k^1_{11}+k^3_{11} &k^1_{12}+k^3_{12} &k^1_{13} &k^1_{14} &k^3_{13} &k^3_{14} \\
&k^1_{22}+k^3_{22} &k^1_{23} &k^1_{24} &k^3_{23} &k^3_{24} \\
& & k^1_{33}+k^2_{11} &k^1_{34}+k^2_{12} &k^2_{13} &k^2_{14} \\
& & &k^1_{44}+k^2_{22} &k^2_{23} &k^2_{24} \\
& & & & k^2_{33}+k^3_{33} &k^2_{24}+k^3_{34} \\
& & & & & k^2_{44}+k^3_{44}
\end{bmatrix}##

The above is the matrix I am trying to understand how it was assembled based on the attached picture..

The numbers in the corners such as 1(1,2) etc in figure 4.6.3 (b) denote the global and element numbers respectively.

There are 2 displacement degrees of freedom (horizontal and vertical) at each node of the element...thanks
 

Attachments

  • DSC_0892.jpg
    DSC_0892.jpg
    30.2 KB · Views: 474
bugatti79 said:
Maybe if I write out the matrix attached for easier read...the blanks indicate symmetry.

##\begin{bmatrix}
k^1_{11}+k^3_{11} &k^1_{12}+k^3_{12} &k^1_{13} &k^1_{14} &k^3_{13} &k^3_{14} \\
&k^1_{22}+k^3_{22} &k^1_{23} &k^1_{24} &k^3_{23} &k^3_{24} \\
& & k^1_{33}+k^2_{11} &k^1_{34}+k^2_{12} &k^2_{13} &k^2_{14} \\
& & &k^1_{44}+k^2_{22} &k^2_{23} &k^2_{24} \\
& & & & k^2_{33}+k^3_{33} &k^2_{24}+k^3_{34} \\
& & & & & k^2_{44}+k^3_{44}
\end{bmatrix}##

The above is the matrix I am trying to understand how it was assembled based on the attached picture..

The numbers in the corners such as 1(1,2) etc in figure 4.6.3 (b) denote the global and element numbers respectively.

There are 2 displacement degrees of freedom (horizontal and vertical) at each node of the element...thanks

After some searching online I have a found an easy way of assembling the global matrix for this problem.

If we focus on element 3 which has global nodes 1 and 3. we can create the element freedom table 'EFT' for this element by the following

2 dof's times global number 1 minus 1=1
2 dof's times the global number 1 =2

2 dof's times global number 3 minus 1=5
2 dof's times the global number 3 =6

( I am interested to know what the above technique is based on)

Thus the EFT is {1,2,5,6}. Similarly for the other 2 elements.

Then one combines the EFT for each element into the global matrx (2 dof's times number elements 3= 6 gives a 6 matrix.)
 
How did you find PF?: Via Google search Hi, I have a vessel I 3D printed to investigate single bubble rise. The vessel has a 4 mm gap separated by acrylic panels. This is essentially my viewing chamber where I can record the bubble motion. The vessel is open to atmosphere. The bubble generation mechanism is composed of a syringe pump and glass capillary tube (Internal Diameter of 0.45 mm). I connect a 1/4” air line hose from the syringe to the capillary The bubble is formed at the tip...
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'
Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...
I'd like to create a thread with links to 3-D Printer resources, including printers and software package suggestions. My motivations are selfish, as I have a 3-D printed project that I'm working on, and I'd like to buy a simple printer and use low cost software to make the first prototype. There are some previous threads about 3-D printing like this: https://www.physicsforums.com/threads/are-3d-printers-easy-to-use-yet.917489/ but none that address the overall topic (unless I've missed...
Back
Top