Figure 9.34 shows the members and idealised support conditions for a roof truss. All three members of which are steel and have the same cross-sectional area such that AE = 12 MN throughout. Using the displacement based finite element method, treating the truss as a pin-jointed plane frame and each

Question

Figure 9.34 shows the members and idealised support conditions for a roof truss. All three members of which are steel and have the same cross-sectional area such that AE = 12 MN throughout. Using the displacement based finite element method, treating the truss as a pin-jointed plane frame and each member as a rod:

(a) assemble the necessary terms in the structural stiffness matrix;
(b) hence, determine the nodal displacements with respect to the global coordinates, for the condition shown in Fig. 9.34.
(c) If, under load, the left support sinks by 5 mm, determine the resulting new nodal displacements, with respect to the global coordinates.

Accepted Answer

(a) Figure 9.35 shows suitable node, dof. and element labelling. Lack of symmetry prevents any advantage being taken to reduce the calculations. However, since both ends of the horizontal member are fixed it is redundant therefore and does not need to be considered further.

All three members have the same AE, hence

[latex]\left(AE ight)^{a} =\left(AE ight)^{b}=AE[/latex]

With reference to §9.7, the element stiffness matrix with respect to global coordinates is given by

[latex][k^{(e)}]=(\frac{AE}{L} )^{(e)}\begin{bmatrix} \cos ^{2} \alpha & & & \ \sin \alpha \cos \alpha & \sin ^{2}\alpha & symmetric & \ -\cos ^{2} \alpha & -\sin \alpha \cos \alpha& \cos ^{2} \alpha & \ -\sin \alpha \cos \alpha & -\sin ^{2}\alpha & \sin \alpha \cos \alpha &\sin ^{2}\alpha \end{bmatrix} ^{(e)}[/latex]

Evaluating the stiffness matrix for both elements:

Element a

[latex]L^{\left(a ight) } =2m , \alpha ^{\left(a ight) }=60° , \cos \alpha ^{\left(a ight)}=1/2 , \sin \alpha ^{\left(a ight)}=\sqrt{3} /2[/latex]

[latex][k^{\left(a ight) }]=\frac{AE}{8}\begin{bmatrix} 1 & \sqrt{3}&-1&-\sqrt{3} \ 3 & 3 &-\sqrt{3} &-3 \-1&-\sqrt{3}&1&\sqrt{3}\-\sqrt{3}&-3& \sqrt{3}&3 \end{bmatrix}[/latex]

Element b

[latex]L^{\left(b ight) } =2\sqrt{3} m , \alpha ^{\left(b ight) }=150° , \cos \alpha ^{\left(b ight)}=-\sqrt{3} /2 , \sin \alpha ^{\left(a ight)}=1 /2[/latex]

[latex][k^{\left(b ight) }]=\frac{AE}{8}\begin{bmatrix} \sqrt{3} &-1&-\sqrt{3} &1 \ -1 & 1/\sqrt{3} &1 &-1/\sqrt{3}&\-\sqrt{3}&1& \sqrt{3}&-1 \1 &-1/\sqrt{3} & -1& 1/\sqrt{3} \end{bmatrix}[/latex]

The structural stiffness matrix can now be multi-assembled using a dof. correspondence table, (ref. §9.10), and will be of order 6 × 6. Only the upper sub-matrices need to be completed, i.e. [latex][K_{αα}] ext{ and } [K_{αβ}][/latex] , since the reactions are not required in this example. The necessary structural governing equations and hence the required structural stiffness matrix are therefore given as

Row/column in [latex][K^{(e)}][/latex]	Row/column in [K]
Row/column in [latex][K^{(e)}][/latex]	a	b
1	3	5
2	4	6
3	1	1
4	2	2

[latex] \begin{bmatrix} X_{1} \\Y_{1} \\X_{2} \\Y_{2} \\X_{3} \\Y_{3} \\end{bmatrix}=\frac{AE}{8}\begin{bmatrix} 1 & &\sqrt{3} & &-1 & &-\sqrt{3}& & & & & \ & \sqrt{3} & &-1 & & && & &-\sqrt{3}& &1 \\sqrt{3} & &3 & &-\sqrt{3} & &-3 & & & & & \ & -1 & & 1/\sqrt{3} & & & & & & 1 &-1/\sqrt{3} \ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\ & & & & & & & & & &\end{bmatrix} \begin{bmatrix} u_{1} \\\upsilon _{1} \\u_{2} \\\upsilon _{2}\ \u_{3} \\\upsilon _{3} \\end{bmatrix}[/latex]

(b) Corresponding to [latex]u_{2}= \upsilon _{2}=u_{3}=\upsilon _{3}=0[/latex] , the partitioned equations reduce to:

[latex]\begin{bmatrix} X _{1} \ Y_{2} \end{bmatrix}= \begin{bmatrix}-20 imes 10^{3} \ -40 imes 10^{3} \end{bmatrix}=\frac{AE}{8}\begin{bmatrix} 2.7321 & 0.7321 \0.7321 & 3.5774 \end{bmatrix} \begin{bmatrix} \upsilon _{1} \ \upsilon _{2} \end{bmatrix}[/latex] i.e. [latex]\left\{p_{α} ight\} = [K_{αα}]\left\{p_{α} ight\}[/latex]

Inverting [latex] [K_{αα}] [/latex] to enable a solution for the displacements from [latex]\left\{p_{α} ight\} = [K_{αα}]^{-1}\left\{p_{α} ight\}[/latex]

[latex]adj [K_{\alpha \alpha }]=\begin{bmatrix} 3.5774& -0.7321 \-0.7321 & 2.7321 \end{bmatrix}[/latex] and [latex]det[K_{\alpha \alpha }]=(AE/8)^{2} 9.2378[/latex]

Then [latex][K_{αα}]^{-1}[/latex]

[latex]=\frac{8}{AE} \begin{bmatrix} 0.3873 &-0,079251 \-0.07925 &0.2958 \end{bmatrix} [/latex] Check: [latex]\frac{8}{AE} \begin{bmatrix} 0.3873 &-0.079251 \-0.07925 &0.2958 \end{bmatrix}\frac{8}{AE} \begin{bmatrix} 2.7321 &0.7321 \0.7321 &3.57741 \end{bmatrix} =[I][/latex]

Hence, the required displacements are given by

[latex]\left\{P_{α } ight\} =[K_{\alpha \alpha }]^{-1}\left\{P_{\alpha } ight\}[/latex]

Substituting,[latex] \begin{bmatrix} u_{1} \ \upsilon _{1} \end{bmatrix} =\frac{8}{12×10^{6}}\begin{bmatrix} 0.3873 & -0.079251 \ -0.07925&0.2958 \end{bmatrix}\begin{bmatrix}-20.10^{3} \ -40.10^{3} \end{bmatrix}=\begin{bmatrix} -3.051 \ -6.83\end{bmatrix}_{mm}[/latex]

The required nodal displacements are therefore [latex]u_{1} = -3.05 \ mm[/latex] and [latex]\upsilon _{1} = -6.83 \ mm[/latex].

(c) With reference to §9.12, for non-zero prescribed displacements, i.e. [latex] \left\{p_{β} ight\} ≠ \left\{0 ight\},[/latex] the full partition of the governing equation is required, namely,

[latex]\left\{P_{\alpha } ight\} =[K_{\alpha \alpha }]\left\{P_{\alpha } ight\}+[K_{\alpha \beta }]\left\{P_{\beta } ight\}[/latex]

Rearranging for the unknown displacements

[latex]\left\{P_{\alpha } ight\} =[K_{\alpha \alpha }]^{-1}\left\{P_{\alpha } ight\}-[K_{\alpha \alpha }]^{-1}[K_{\alpha \beta }]\left\{P_{\beta } ight\}[/latex]

Evaluating

[latex][K_{\alpha \alpha }]^{-1}[K_{\alpha \beta }]=\begin{bmatrix} 0.3873& -0.079251 \ -0.07925&0.2958 \end{bmatrix}\begin{bmatrix} -1 & -1.7321 &-1.7321 &1\ -1.7321 & -3 &1 & -0.5773 \end{bmatrix} = \begin{bmatrix} -0.25 & -0.4331 &-0.75 &0.4331\ -0.4331 & -0.75&0.4331& -0.25 \end{bmatrix}[/latex]

and [latex][K_{\alpha \alpha }]^{-1}[K_{\alpha \beta }]\left\{P_{\beta } ight\}=\begin{bmatrix} -0.25 & -0.4331 &-0.75 &0.4331\ -0.4331 & -0.75&0.4331& -0.25 \end{bmatrix} \begin{bmatrix} 0 \-5 imes 10^{-3}\0\0 \end{bmatrix} =\begin{bmatrix} 2.1655\3.75 \end{bmatrix} _{mm}[/latex]

Recalling from part (b) that

[latex][K_{\alpha \alpha }]^{-1}\left\{P_{\alpha } ight\}=\begin{bmatrix} -3.05\-6.83 \end{bmatrix}_{mm}[/latex]

and substituting into the above rearranged governing equation,

i.e [latex]\begin{bmatrix} u_{1} \ \upsilon _{1} \end{bmatrix}=\begin{bmatrix} -3.05\-6.83 \end{bmatrix}-\begin{bmatrix} 2.1655\3.75\end{bmatrix}=\begin{bmatrix} -5.22\-10.58\end{bmatrix}_{mm}[/latex]

yields the required new nodal displacements, namely, [latex]u_{1} = -5.22 \ mm[/latex] and [latex]\upsilon _{1} = - 10.58 \ mm[/latex].

Question 9.2: Figure 9.34 shows the members and idealised support conditio...

Related Answered Questions

Verified Answer:

Verified Answer:

Verified Answer:

Derive the stiffness matrix in global coordinates for a three-node triangular membrane element for plane stress analysis. Assume that the elastic modulus, E, and thickness, t, are constant throughout, and that the displacement functions are u(x, y) =α1 + α2x + α3y υ(x, y) = α4 + α5x + α6y ...

Verified Answer:

Verified Answer:

Figure 9.31 shows a planar steel support structure, all three members of which have the same axial stiffness, such that AE/L = 20 MN/m throughout. Using the displacement based finite element method and treating each member as a rod: (a) assemble the necessary terms in the structural stiffness ...

Verified Answer: