User:Eml4500.f08.Lulz.js/HW3

Axial displacement calculation for local node 2
During meeting 12, a calculation for the axial displacement of local node 1 for an arbitrary element e was considered. What follows is a similar calculation for local node 2 of the same element.

First, we define the displacement vector at local node 2 (Please note the [#] subscript was used in place of the "box" subscript used in class):

$$ \vec{\mathbf{d}^{(e)}_{[2]}} = d_{3}^{(e)} \vec{i} + d_{4}^{(e)} \vec{j} $$

Next, the axial displacement at node 2 is defined: $$ q_{2}^{(e)} = \vec{\mathbf{d}^{(e)}_{[2]}}\cdot\vec{\tilde{i}} $$

$$ q_{2}^{(e)} = (d_{3}^{(e)}\vec{i} + d_{4}^{(e)}\vec{j})\vec{\tilde{i}} $$

$$ q_{2}^{(e)} = d_{3}^{(e)}(\vec{i}\cdot\vec{\tilde{i}}) + d_{4}^{(e)}(\vec{j}\cdot\vec{\tilde{i}}) $$

From a previous homework question and the class notes, the following relations are valid:

$$\vec{\tilde{i}}\cdot\vec{i}=\cos\theta^{(e)}$$ $$\vec{\tilde{i}}\cdot\vec{j}=\sin\theta^{(e)}$$

Substituting,

$$ q_{2}^{(e)} = d_{3}^{(e)}\cos \theta^{(e)} + d_{4}^{(e)}\sin\theta^{(e)} $$

From the notes, the following relations are valid:

$$\cos\theta^{\left(e\right)} = l^{(e)}$$ $$\sin\theta^{\left(e\right)} = m^{(e)}$$

Substituting,

$$ q_{2}^{(e)} = d_{3}^{(e)}l^{(e)} + d_{4}^{(e)}m^{(e)} $$

Putting the above equation yields the following result, consistent with that of the axial displacement of local node 1:

$$ q_{2}^{(e)} = \begin{bmatrix} l^{(e)} & m^{(e)} \end{bmatrix}\begin{Bmatrix} d_{3}^{(e)}\\ d_{4}^{(e)} \end{Bmatrix} $$

Proof of $$\mathbf{k}^{(e)} = \mathbf{T}^{(e)^{T}} \mathbf{\hat{k}}^{(e)}\mathbf{T}^{(e)}$$
In this question we are asked to prove that the above matrix product generates the generic element stiffness matrix $$\mathbf{k}^{(e)}$$.

First, the three matrices used to generate the product:

$$\mathbf{T}^{(e)} = \begin{bmatrix} l^{(e)}&m^{(e)}&0&0 \\ 0&0&l^{(e)}&m^{(e)} \end{bmatrix}$$

$$\mathbf{T}^{(e)^{T}} = \begin{bmatrix} l^{(e)}&0 \\ m^{(e)}&0 \\ 0&l^{(e)} \\ 0&m^{(e)} \\ \end{bmatrix}$$

$$\mathbf{\hat{k}}^{(e)} = \begin{bmatrix} 1&-1 \\ -1&1 \\ \end{bmatrix}$$

The product of the first two matrices:

$$\mathbf{\hat{k}}^{(e)}\mathbf{T}^{(e)} = \begin{bmatrix} 1&-1 \\ -1&1 \\ \end{bmatrix} \begin{bmatrix} l^{(e)}&m^{(e)}&0&0 \\ 0&0&l^{(e)}&m^{(e)} \end{bmatrix} = \begin{bmatrix} l^{(e)}&m^{(e)}&-l^{(e)}&-m^{(e)} \\ -l^{(e)}&-m^{(e)}&l^{(e)}&m^{(e)} \\ \end{bmatrix}$$

The entire product multiplied out:

$$\mathbf{T}^{(e)^{T}}\mathbf{\hat{k}}^{(e)}\mathbf{T}^{(e)} = \begin{bmatrix} l^{(e)}&0 \\ m^{(e)}&0 \\ 0&l^{(e)} \\ 0&m^{(e)} \\ \end{bmatrix} \begin{bmatrix} l^{(e)}&m^{(e)}&-l^{(e)}&-m^{(e)} \\ -l^{(e)}&-m^{(e)}&l^{(e)}&m^{(e)} \\ \end{bmatrix} $$ $$ = \begin{bmatrix} l^{(e)^{2}} & m^{(e)}l^{(e)} & -l^{(e)^{2}} & -m^{(e)}l^{(e)} \\ m^{(e)}l^{(e)} & m^{(e)^{2}} & -m^{(e)}l^{(e)} & -m^{(e)^{2}}\\ -l^{(e)^{2}} & -m^{(e)}l^{(e)} & l^{(e)^{2}} & m^{(e)}l^{(e)}\\ -m^{(e)}l^{(e)} & -m^{(e)^{2}} & m^{(e)}l^{(e)} & m^{(e)^{2}} \\ \end{bmatrix} $$

As shown in previous homework question, the element stiffness matrix:

$$\mathbf{k}^{(e)} = \begin{bmatrix} l^{(e)^{2}} & m^{(e)}l^{(e)} & -l^{(e)^{2}} & -m^{(e)}l^{(e)} \\ m^{(e)}l^{(e)} & m^{(e)^{2}} & -m^{(e)}l^{(e)} & -m^{(e)^{2}}\\ -l^{(e)^{2}} & -m^{(e)}l^{(e)} & l^{(e)^{2}} & m^{(e)}l^{(e)}\\ -m^{(e)}l^{(e)} & -m^{(e)^{2}} & m^{(e)}l^{(e)} & m^{(e)^{2}} \\ \end{bmatrix}$$

Thus, it is easily shown that the following relation is indeed true:

$$\mathbf{k}^{(e)} = \mathbf{T}^{(e)^{T}} \mathbf{\hat{k}}^{(e)}\mathbf{T}^{(e)}$$