Stiffness Matrix Method: Symmetry vs Introducing a new node / joint

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
Master1022
Messages
590
Reaction score
116
TL;DR
This question is about comparing two methods in a stiffness matrix analysis question - using symmetry on half of the structure vs. introducing a new joint at the load
Hi,

In the question outlined in the images (apologies for the poor quality of the scans), the chosen solution has opted to use a symmetry argument and proceed from there.

ScannableDocumentStiffnessMatrix.jpg
StiffnessMatrix2.jpg


Question is from "Structures: theory and analysis" by Williams & Todd

My question is: How could we approach the same problem by introducing an extra node at the point where the load [itex]P[/itex] acts?

Method: I believe that all we would need to do if we were considering an extra joint would be to effectively double the entry in the bottom right corner of the global stiffness matrix: (i.e.

[tex]K =<br /> \begin{pmatrix}<br /> \frac{4EI}{l} & \frac{-6EI}{l^2} \\<br /> \frac{-6EI}{l^2} & (\frac{24EI}{l^3} + k)\\<br /> \end{pmatrix}[/tex]

and changing the force vector to: [tex]F = \begin{bmatrix} 0 \\ -P \\ \end{bmatrix}[/tex]

I thought of doubling the bottom right entry of K as we would be adding it to itself again in another 2 x 2 global stiffness matrix and that would be the only parameter that we care about within [itex]K_{BC}[/itex]

However, solving this yields different answers (letting EI = 1 and l = 1 for ease of typing):
[tex]\begin{bmatrix} \theta_{z_A} \\ u_{y_b} \\ \end{bmatrix} = \frac{-P}{15 + k} \begin{bmatrix} \frac{3}{2} \\ 1 \\ \end{bmatrix}[/tex].

Can anyone see where I have gone wrong in my thinking here?

Thanks in advance
 
Physics news on Phys.org
Update: I realized that I added up the stiffness matrices of the right and left half incorrectly.

The global stiffness matrix (for variables [itex]\theta_{z_A}, u_{y_B}, \theta_{z_C}[/itex] respectively) should have been:
[tex] \begin{pmatrix}<br /> \frac{4EI}{l^3} & \frac{-6EI}{l^2} & 0 \\<br /> \frac{-6EI}{l^2} & \frac{24EI}{l^3} + k & \frac{6EI}{l^2} \\<br /> 0 & \frac{6EI}{l^2} & \frac{4EI}{l^3}<br /> \end{pmatrix}[/tex]
Using this matrix leads to the same result as in the book.