Constant Strain Triangle (CST) — Solved Questions

Chapter 4.6 — IOE Past Questions with Worked FEM Setups and Results
Collected & Prepared By: Er. Sanjay Kumar Sah  |  Dept. of Civil Engineering, Pulchowk Campus

Question Index

Reference Formulas Used Repeatedly
$$2A = \begin{vmatrix}1 & x_1 & y_1 \\ 1 & x_2 & y_2 \\ 1 & x_3 & y_3\end{vmatrix}, \qquad [B] = \frac{1}{2A}\begin{bmatrix}b_1&b_2&b_3&0&0&0\\0&0&0&c_1&c_2&c_3\\c_1&c_2&c_3&b_1&b_2&b_3\end{bmatrix}$$
$$[D] = \frac{E}{1-\nu^2}\begin{bmatrix}1&\nu&0\\\nu&1&0\\0&0&\frac{1-\nu}{2}\end{bmatrix}, \qquad [K] = tA[B]^T[D][B]$$
$$\{f_t\}_{ij} = \frac{tl_{ij}}{2}\begin{Bmatrix}t_x\\t_y\\t_x\\t_y\end{Bmatrix}, \qquad \{f_b\} = \frac{\gamma At}{3}\begin{Bmatrix}0\\-1\\0\\-1\\0\\-1\end{Bmatrix}$$
Working note: The sketches in the scanned questions are slightly ambiguous. The numerical solutions below follow the support and coordinate interpretation consistent with the figure captions and the CST notes page already in this workspace. For the unconstrained element-only question, the singularity is stated explicitly instead of forcing an artificial displacement solution.

Question 1 — Equilateral CST with self-weight and edge traction

2080 Chaitra, Q10, Group B
Plane stress, $E = 210 \times 10^3$ MPa, $\nu = 0.3$, $t = 1$ cm, $\gamma = 78.5$ kN/m$^3$, equilateral triangle with side $100$ mm, and edge traction $10$ kN/m$^2$ on face 2-3. Nodes 1 and 2 are restrained.
Q1 — Equilateral CST | side = 100 mm x y O 1 (0, 0) 2 (100, 0) 3 (50, 86.6) 10 kN/m² (on edge 2–3) γ = 78.5 kN/m³ 100 mm Material Data E = 210 × 10³ MPa ν = 0.3 t = 10 mm γ = 78.5 kN/m³ Plane stress

Figure Q1 — Equilateral CST (100 mm side), pin-supported at nodes 1 & 2, with 10 kN/m² traction on edge 2–3

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(-0.05,\;0)\;\text{m}, \quad \text{Node 2: }(0.05,\;0)\;\text{m}, \quad \text{Node 3: }(0,\;0.09)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & -0.05 & 0 \\ 1 & 0.05 & 0 \\ 1 & 0 & 0.09\end{vmatrix}$$ $$= -0.05(0-0.09) + 0.05(0.09-0) + 0(0-0)$$ $$= 0.00866\;\text{m}^2$$
$$\boxed{A = 0.00433\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.09 = -0.09$$ $$b_2 = y_3 - y_1 = 0.09 - 0 = 0.09$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0 - 0.05 = -0.05$$ $$c_2 = x_1 - x_3 = -0.05 - 0 = -0.05$$ $$c_3 = x_2 - x_1 = 0.05 - -0.05 = 0.1$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-10.00 & 0.00 & 10.00 & 0.00 & 0.00 & 0.00 \\ 0.00 & -5.77 & 0.00 & -5.77 & 0.00 & 11.55 \\ -5.77 & -10.00 & -5.77 & 10.00 & 11.55 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}230769230.77 & 69230769.23 & 0.00 \\ 69230769.23 & 230769230.77 & 0.00 \\ 0.00 & 0.00 & 80769230.77\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}1115814.53 & 375000.00 & -882647.01 & -28846.15 & -233167.53 & -346153.85 \\ 375000.00 & 682827.23 & 28846.15 & -16634.30 & -403846.15 & -666192.93 \\ -882647.01 & 28846.15 & 1115814.53 & -375000.00 & -233167.53 & 346153.85 \\ -28846.15 & -16634.30 & -375000.00 & 682827.23 & 403846.15 & -666192.93 \\ -233167.53 & -403846.15 & -233167.53 & 403846.15 & 466335.05 & 0.00 \\ -346153.85 & -666192.93 & 346153.85 & -666192.93 & 0.00 & 1332385.86\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

(a) Body force (self-weight in $-y$ direction):

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \frac{78.5 \times 0.00433 \times 0.01}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -1.1330 \times 10^{-3} \\ 0 \\ -1.1330 \times 10^{-3} \\ 0 \\ -1.1330 \times 10^{-3}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 2–3: The edge length is $l_{23} = 0.1$ m. The outward normal to edge 2–3 has direction $(0.87, 0.5)$, so the traction components are:

$$t_x = 10 \times 0.87 = 8.66\;\text{kN/m}^2, \quad t_y = 10 \times 0.5 = 5\;\text{kN/m}^2$$
$$\{f_t\}_{23} = \frac{t \cdot l_{23}}{2} \begin{Bmatrix} 0\\0\\t_x\\t_y\\t_x\\t_y \end{Bmatrix} = \begin{Bmatrix}0 \\ 0 \\ 4.3300 \times 10^{-3} \\ 2.5000 \times 10^{-3} \\ 4.3300 \times 10^{-3} \\ 2.5000 \times 10^{-3}\end{Bmatrix}\;\text{kN}$$

(c) Combined load vector:

$$\{F\} = \{f_b\} + \{f_t\} = \begin{Bmatrix}0 \\ -1.1330 \times 10^{-3} \\ 4.3300 \times 10^{-3} \\ 1.3670 \times 10^{-3} \\ 4.3300 \times 10^{-3} \\ 1.3670 \times 10^{-3}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Nodes 1 and 2 are restrained (pinned): $u_1 = v_1 = u_2 = v_2 = 0$.

Free DOFs: $u_3$ (DOF 5) and $v_3$ (DOF 6). Extract the $2 \times 2$ sub-system.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}466335.05 & 0.00 \\ 0.00 & 1332385.86\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}4.3300 \times 10^{-3} \\ 1.3670 \times 10^{-3}\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_3 = 9.2852 \times 10^{-9}\;\text{m}, \qquad v_3 = 1.0260 \times 10^{-9}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ 1.1847 \times 10^{-8} \\ 1.0722 \times 10^{-7} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = 0.82\;\text{kN/m²}, \quad \sigma_y = 2.73\;\text{kN/m²}, \quad \tau_{xy} = 8.66\;\text{kN/m²}}$$
For a CST element, strain and stress are constant throughout the element, so the centroid stress is identical to the stress computed above.

Question 5 — Element stiffness matrix and traction load vector

2078 Chaitra, Q6, Group B
Steel plate of thickness $10$ mm with a uniform traction of $5$ kN/m$^2$ on the left edge. Determine the CST element stiffness matrix, load vector, and nodal displacements.
Q5 — Right-Triangle CST | 500 × 500 mm | Free element x y O 1 (0, 0) 2 (0, 500) 3 (500, 500) 5 kN/m² 500 mm 500 mm Material Data E = 2 × 10⁵ MPa ν = 0.3 t = 10 mm Plane stress No supports (free element)

Figure Q5 — Right-triangle CST (free element) with 5 kN/m² traction on edge 1–2 (left vertical edge)

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0,\;0.5)\;\text{m}, \quad \text{Node 3: }(0.5,\;0.5)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0 & 0.5 \\ 1 & 0.5 & 0.5\end{vmatrix}$$ $$= 0(0.5-0.5) + 0(0.5-0) + 0.5(0-0.5)$$ $$= -0.25\;\text{m}^2$$
$$\boxed{A = 0.12\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0.5 - 0.5 = 0$$ $$b_2 = y_3 - y_1 = 0.5 - 0 = 0.5$$ $$b_3 = y_1 - y_2 = 0 - 0.5 = -0.5$$
$$c_1 = x_3 - x_2 = 0.5 - 0 = 0.5$$ $$c_2 = x_1 - x_3 = 0 - 0.5 = -0.5$$ $$c_3 = x_2 - x_1 = 0 - 0 = 0$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-0.00 & -0.00 & -2.00 & -0.00 & 2.00 & -0.00 \\ -0.00 & -2.00 & -0.00 & 2.00 & -0.00 & -0.00 \\ -2.00 & -0.00 & 2.00 & -2.00 & -0.00 & 2.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}219780219.78 & 65934065.93 & 0.00 \\ 65934065.93 & 219780219.78 & 0.00 \\ 0.00 & 0.00 & 76923076.92\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}384615.38 & 0.00 & -384615.38 & 384615.38 & 0.00 & -384615.38 \\ 0.00 & 1098901.10 & 329670.33 & -1098901.10 & -329670.33 & 0.00 \\ -384615.38 & 329670.33 & 1483516.48 & -714285.71 & -1098901.10 & 384615.38 \\ 384615.38 & -1098901.10 & -714285.71 & 1483516.48 & 329670.33 & -384615.38 \\ 0.00 & -329670.33 & -1098901.10 & 329670.33 & 1098901.10 & 0.00 \\ -384615.38 & 0.00 & 384615.38 & -384615.38 & 0.00 & 384615.38\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vector

Traction of $5$ kN/m² acts in the $+x$ direction on edge 1–2 (left vertical edge). Edge length $l_{12} = 0.5$ m.

$$\{F\} = \frac{t \cdot l_{12}}{2} \begin{Bmatrix} t_x \\ t_y \\ t_x \\ t_y \\ 0 \\ 0 \end{Bmatrix} = \frac{0.01 \times 0.5}{2} \begin{Bmatrix} 5\\0\\5\\0\\0\\0 \end{Bmatrix} = \begin{Bmatrix}1.2500 \times 10^{-2} \\ 0 \\ 1.2500 \times 10^{-2} \\ 0 \\ 0 \\ 0\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

No support conditions are given — this is a free (unsupported) element.

Step 9 — Solve for Unknown Displacements

Since no boundary conditions are imposed, the stiffness matrix is singular (rank 3). A unique displacement solution cannot be obtained from the free element alone.

Important: The element stiffness matrix [K] and load vector {F} are well-defined. However, without support conditions, the system has 3 rigid-body modes, making [K] singular. To obtain displacements, support conditions must first be specified.

Step 10 — Compute Element Strains and Stresses

Without support conditions, strains and stresses cannot be computed for the free element.

Important: A free single CST element has rank 3 (3 rigid-body modes: 2 translations + 1 rotation). The $6 \times 6$ stiffness matrix is singular. To solve for displacements, at least 3 independent support conditions must be imposed. The stiffness matrix $[K]$ and load vector $\{F\}$ computed above are the required answers.

Question 6 — Loaded right-triangle CST with supports

2078 Kartik, Q6, Group B
Plane stress, $t = 10$ mm, $\gamma = 78.5$ kN/m$^3$, $E = 2 \times 10^5$ MPa, $G = 105 \times 10^3$ MPa. Determine the stiffness matrix, nodal load vector and nodal displacement.
Q6 — Right-Triangle CST | 800 × 600 mm x y O 1 (0, 0) 2 (800, 0) 3 (0, 600) 5 kN/m² 800 mm 600 mm Material Data E = 2 × 10⁵ MPa G = 105 × 10³ MPa ν (from G) t = 10 mm Plane stress

Figure Q6 — Right-triangle CST with nodes 1 & 2 pin-supported, horizontal traction 5 kN/m² on edge 1–3

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0.8,\;0)\;\text{m}, \quad \text{Node 3: }(0,\;0.6)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0.8 & 0 \\ 1 & 0 & 0.6\end{vmatrix}$$ $$= 0(0-0.6) + 0.8(0.6-0) + 0(0-0)$$ $$= 0.48\;\text{m}^2$$
$$\boxed{A = 0.24\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.6 = -0.6$$ $$b_2 = y_3 - y_1 = 0.6 - 0 = 0.6$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0 - 0.8 = -0.8$$ $$c_2 = x_1 - x_3 = 0 - 0 = 0$$ $$c_3 = x_2 - x_1 = 0.8 - 0 = 0.8$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-1.25 & 0.00 & 1.25 & 0.00 & 0.00 & 0.00 \\ 0.00 & -1.67 & 0.00 & 0.00 & 0.00 & 1.67 \\ -1.67 & -1.25 & 0.00 & 1.25 & 1.67 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}200454545.45 & -9545454.55 & 0.00 \\ -9545454.55 & 200454545.45 & 0.00 \\ 0.00 & 0.00 & 105000000.00\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}1451704.55 & 477272.73 & -751704.55 & -525000.00 & -700000.00 & 47727.27 \\ 477272.73 & 1730113.64 & 47727.27 & -393750.00 & -525000.00 & -1336363.64 \\ -751704.55 & 47727.27 & 751704.55 & 0.00 & 0.00 & -47727.27 \\ -525000.00 & -393750.00 & 0.00 & 393750.00 & 525000.00 & 0.00 \\ -700000.00 & -525000.00 & 0.00 & 525000.00 & 700000.00 & 0.00 \\ 47727.27 & -1336363.64 & -47727.27 & 0.00 & 0.00 & 1336363.64\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

Poisson's ratio from given data: $G = E/[2(1+\nu)]$, so $\nu = E/(2G) - 1 = 200000/(2 \times 105000) - 1 = -0.05$.

This gives a near-zero Poisson's ratio. Some exam versions may intend $\nu \approx 0.3$ for steel — use the value your question sheet specifies.

(a) Body force:

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -6.2800 \times 10^{-2} \\ 0 \\ -6.2800 \times 10^{-2} \\ 0 \\ -6.2800 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 1–3 ($l_{13} = 0.6$ m, acting in $+x$):

$$\{f_t\} = \frac{t \cdot l_{13}}{2} \begin{Bmatrix} 5\\0\\0\\0\\5\\0 \end{Bmatrix} = \begin{Bmatrix}1.5000 \times 10^{-2} \\ 0 \\ 0 \\ 0 \\ 1.5000 \times 10^{-2} \\ 0\end{Bmatrix}\;\text{kN}$$
$$\{F\} = \{f_b\} + \{f_t\} = \begin{Bmatrix}1.5000 \times 10^{-2} \\ -6.2800 \times 10^{-2} \\ 0 \\ -6.2800 \times 10^{-2} \\ 1.5000 \times 10^{-2} \\ -6.2800 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Nodes 1 and 2 are restrained: $u_1 = v_1 = u_2 = v_2 = 0$. Free DOFs: $u_3, v_3$.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}700000.00 & 0.00 \\ 0.00 & 1336363.64\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}1.5000 \times 10^{-2} \\ -6.2800 \times 10^{-2}\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_3 = 2.1429 \times 10^{-8}\;\text{m}, \qquad v_3 = -4.6993 \times 10^{-8}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ -7.8322 \times 10^{-8} \\ 3.5714 \times 10^{-8} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = 0.75\;\text{kN/m²}, \quad \sigma_y = -15.7\;\text{kN/m²}, \quad \tau_{xy} = 3.75\;\text{kN/m²}}$$

Question 7 — Triangular plate with body force and side traction

2076 Bhadra, Q8, Group B
Plane stress, $E = 210$ GPa, thickness $15$ mm, unit weight $80$ kN/m$^3$, $\nu = 0.3$, and side traction $15$ kN/m$^2$ on the right edge. Determine displacement and stress.
Q7 — Triangular Plate | 750 × 500 mm x y O 1 (0, 0) 2 (750, 0) 3 (750, 500) 15 kN/m² γ = 80 kN/m³ 750 mm 500 mm Material Data E = 210 GPa ν = 0.3 t = 15 mm γ = 80 kN/m³ Plane stress

Figure Q7 — Triangular plate with body force and 15 kN/m² horizontal traction on edge 2–3

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0.75,\;0)\;\text{m}, \quad \text{Node 3: }(0.75,\;0.5)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0.75 & 0 \\ 1 & 0.75 & 0.5\end{vmatrix}$$ $$= 0(0-0.5) + 0.75(0.5-0) + 0.75(0-0)$$ $$= 0.38\;\text{m}^2$$
$$\boxed{A = 0.19\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.5 = -0.5$$ $$b_2 = y_3 - y_1 = 0.5 - 0 = 0.5$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0.75 - 0.75 = 0$$ $$c_2 = x_1 - x_3 = 0 - 0.75 = -0.75$$ $$c_3 = x_2 - x_1 = 0.75 - 0 = 0.75$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-1.33 & 0.00 & 1.33 & 0.00 & 0.00 & 0.00 \\ 0.00 & 0.00 & 0.00 & -2.00 & 0.00 & 2.00 \\ 0.00 & -1.33 & -2.00 & 1.33 & 2.00 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}230769230.77 & 69230769.23 & 0.00 \\ 69230769.23 & 230769230.77 & 0.00 \\ 0.00 & 0.00 & 80769230.77\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}1153846.15 & 0.00 & -1153846.15 & 519230.77 & 0.00 & -519230.77 \\ 0.00 & 403846.15 & 605769.23 & -403846.15 & -605769.23 & 0.00 \\ -1153846.15 & 605769.23 & 2062500.00 & -1125000.00 & -908653.85 & 519230.77 \\ 519230.77 & -403846.15 & -1125000.00 & 3000000.00 & 605769.23 & -2596153.85 \\ 0.00 & -605769.23 & -908653.85 & 605769.23 & 908653.85 & 0.00 \\ -519230.77 & 0.00 & 519230.77 & -2596153.85 & 0.00 & 2596153.85\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

(a) Body force ($\gamma = 80$ kN/m³):

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -7.5000 \times 10^{-2} \\ 0 \\ -7.5000 \times 10^{-2} \\ 0 \\ -7.5000 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 2–3 (right vertical edge, $l_{23} = 0.5$ m). The traction acts in the $-x$ direction with magnitude $15$ MPa $= 15000$ kN/m².

$$\{f_t\}_{23} = \frac{t \cdot l_{23}}{2} \begin{Bmatrix} 0\\0\\t_x\\0\\t_x\\0 \end{Bmatrix} = \begin{Bmatrix}0 \\ 0 \\ -56.25 \\ 0 \\ -56.25 \\ 0\end{Bmatrix}\;\text{kN}$$

(c) Combined:

$$\{F\} = \begin{Bmatrix}0 \\ -7.5000 \times 10^{-2} \\ -56.25 \\ -7.5000 \times 10^{-2} \\ -56.25 \\ -7.5000 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Nodes 1 and 2 are restrained: $u_1 = v_1 = u_2 = v_2 = 0$. Free DOFs: $u_3, v_3$.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}908653.85 & 0.00 \\ 0.00 & 2596153.85\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}-56.25 \\ -7.5000 \times 10^{-2}\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_3 = -6.1905 \times 10^{-5}\;\text{m}, \qquad v_3 = -2.8889 \times 10^{-8}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ -5.7778 \times 10^{-8} \\ -1.2381 \times 10^{-4} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = -4\;\text{kN/m²}, \quad \sigma_y = -13.33\;\text{kN/m²}, \quad \tau_{xy} = -10000\;\text{kN/m²}}$$

Question 10 — Stress at centroid of fixed triangular plate

2074 Magh, Q8, Group B
Steel plate of thickness $10$ mm, $E = 200$ GPa, $\nu = 0.3$, self-weight $78.5$ kN/m$^3$, and uniform load $5$ kN/m$^2$ on the bottom edge. Find the stress at the centroid using CST.
Q10 — Fixed-Edge Triangular Plate | 800 × 600 mm x y O 1 (0, 0) 2 (800, 0) 3 (0, 600) 5 kN/m² 800 mm 600 mm Material Data E = 200 GPa ν = 0.3 t = 10 mm γ = 78.5 kN/m³ Plane stress fixed edge

Figure Q10 — Fixed left-edge triangular plate with 5 kN/m² downward traction on edge 1–2

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0.8,\;0)\;\text{m}, \quad \text{Node 3: }(0,\;0.6)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0.8 & 0 \\ 1 & 0 & 0.6\end{vmatrix}$$ $$= 0(0-0.6) + 0.8(0.6-0) + 0(0-0)$$ $$= 0.48\;\text{m}^2$$
$$\boxed{A = 0.24\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.6 = -0.6$$ $$b_2 = y_3 - y_1 = 0.6 - 0 = 0.6$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0 - 0.8 = -0.8$$ $$c_2 = x_1 - x_3 = 0 - 0 = 0$$ $$c_3 = x_2 - x_1 = 0.8 - 0 = 0.8$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-1.25 & 0.00 & 1.25 & 0.00 & 0.00 & 0.00 \\ 0.00 & -1.67 & 0.00 & 0.00 & 0.00 & 1.67 \\ -1.67 & -1.25 & 0.00 & 1.25 & 1.67 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}219780219.78 & 65934065.93 & 0.00 \\ 65934065.93 & 219780219.78 & 0.00 \\ 0.00 & 0.00 & 76923076.92\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}1336996.34 & 714285.71 & -824175.82 & -384615.38 & -512820.51 & -329670.33 \\ 714285.71 & 1753663.00 & -329670.33 & -288461.54 & -384615.38 & -1465201.47 \\ -824175.82 & -329670.33 & 824175.82 & 0.00 & 0.00 & 329670.33 \\ -384615.38 & -288461.54 & 0.00 & 288461.54 & 384615.38 & 0.00 \\ -512820.51 & -384615.38 & 0.00 & 384615.38 & 512820.51 & 0.00 \\ -329670.33 & -1465201.47 & 329670.33 & 0.00 & 0.00 & 1465201.47\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

(a) Body force:

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -6.2800 \times 10^{-2} \\ 0 \\ -6.2800 \times 10^{-2} \\ 0 \\ -6.2800 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 1–2 (bottom edge, $l_{12} = 0.8$ m). Uniform pressure $5$ MPa $= 5000$ kN/m² acting downward ($-y$).

$$\{f_t\}_{12} = \frac{t \cdot l_{12}}{2} \begin{Bmatrix} 0\\-5000\\0\\-5000\\0\\0 \end{Bmatrix} = \begin{Bmatrix}0 \\ -20 \\ 0 \\ -20 \\ 0 \\ 0\end{Bmatrix}\;\text{kN}$$

(c) Combined:

$$\{F\} = \begin{Bmatrix}0 \\ -20.0628 \\ 0 \\ -20.0628 \\ 0 \\ -6.2800 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Left edge is fixed (wall), meaning nodes 1 and 3 (both at $x = 0$) are fully restrained: $u_1 = v_1 = u_3 = v_3 = 0$.

Free DOFs: $u_2, v_2$ (node 2 only).

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}824175.82 & 0.00 \\ 0.00 & 288461.54\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}0 \\ -20.0628\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_2 = 0\;\text{m}, \qquad v_2 = -6.9551 \times 10^{-5}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ 0 \\ -8.6939 \times 10^{-5} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = 0\;\text{kN/m²}, \quad \sigma_y = 0\;\text{kN/m²}, \quad \tau_{xy} = -6687.6\;\text{kN/m²}}$$
For a CST element the stress is constant everywhere, so the centroid stress is simply the element stress computed above.

Question 11 — Nodal displacements and stresses under 25 kN/m^2 side load

2073 Bhadra, Q10, Group B
Thickness $8$ mm, plane stress, $E = 210 \times 10^5$ MPa, $\nu = 0.3$, self-weight $78.5$ kN/m$^3$, base $400$ mm, height $400$ mm, and horizontal traction $25$ kN/m$^2$ on the right edge. Determine nodal displacements and stresses.
Q11 — Triangular Plate | 400 × 400 mm x y O 1 (0, 0) 2 (400, 0) 3 (400, 400) 25 kN/m² γ = 78.5 kN/m³ 400 mm 400 mm Material Data E = 2.1 × 10⁵ MPa ν = 0.3 t = 8 mm γ = 78.5 kN/m³ Plane stress

Figure Q11 — Right-triangle CST with 25 kN/m² horizontal traction on edge 2–3 and self-weight

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0.4,\;0)\;\text{m}, \quad \text{Node 3: }(0.4,\;0.4)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0.4 & 0 \\ 1 & 0.4 & 0.4\end{vmatrix}$$ $$= 0(0-0.4) + 0.4(0.4-0) + 0.4(0-0)$$ $$= 0.16\;\text{m}^2$$
$$\boxed{A = 0.08\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.4 = -0.4$$ $$b_2 = y_3 - y_1 = 0.4 - 0 = 0.4$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0.4 - 0.4 = 0$$ $$c_2 = x_1 - x_3 = 0 - 0.4 = -0.4$$ $$c_3 = x_2 - x_1 = 0.4 - 0 = 0.4$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-2.50 & 0.00 & 2.50 & 0.00 & 0.00 & 0.00 \\ 0.00 & 0.00 & 0.00 & -2.50 & 0.00 & 2.50 \\ 0.00 & -2.50 & -2.50 & 2.50 & 2.50 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}230769230.77 & 69230769.23 & 0.00 \\ 69230769.23 & 230769230.77 & 0.00 \\ 0.00 & 0.00 & 80769230.77\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}923076.92 & 0.00 & -923076.92 & 276923.08 & 0.00 & -276923.08 \\ 0.00 & 323076.92 & 323076.92 & -323076.92 & -323076.92 & 0.00 \\ -923076.92 & 323076.92 & 1246153.85 & -600000.00 & -323076.92 & 276923.08 \\ 276923.08 & -323076.92 & -600000.00 & 1246153.85 & 323076.92 & -923076.92 \\ 0.00 & -323076.92 & -323076.92 & 323076.92 & 323076.92 & 0.00 \\ -276923.08 & 0.00 & 276923.08 & -923076.92 & 0.00 & 923076.92\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

(a) Body force:

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -1.6747 \times 10^{-2} \\ 0 \\ -1.6747 \times 10^{-2} \\ 0 \\ -1.6747 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 2–3 (right vertical edge, $l_{23} = 0.4$ m). Uniform $25$ MPa acting in the $-x$ direction:

$$\{f_t\}_{23} = \frac{t \cdot l_{23}}{2} \begin{Bmatrix} 0\\0\\-25000\\0\\-25000\\0 \end{Bmatrix} = \begin{Bmatrix}0 \\ 0 \\ -40 \\ 0 \\ -40 \\ 0\end{Bmatrix}\;\text{kN}$$

(c) Combined:

$$\{F\} = \begin{Bmatrix}0 \\ -1.6747 \times 10^{-2} \\ -40 \\ -1.6747 \times 10^{-2} \\ -40 \\ -1.6747 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Nodes 1 and 2 are restrained (pinned): $u_1 = v_1 = u_2 = v_2 = 0$. Free DOFs: $u_3, v_3$.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}323076.92 & 0.00 \\ 0.00 & 923076.92\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}-40 \\ -1.6747 \times 10^{-2}\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_3 = -1.2381 \times 10^{-4}\;\text{m}, \qquad v_3 = -1.8142 \times 10^{-8}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ -4.5356 \times 10^{-8} \\ -3.0952 \times 10^{-4} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = -3.14\;\text{kN/m²}, \quad \sigma_y = -10.47\;\text{kN/m²}, \quad \tau_{xy} = -25000\;\text{kN/m²}}$$

Question 15 — Thin plate with nodal point loads

2071 Bhadra, Q4c, Group B
Plate thickness $0.3$ in, Poisson's ratio $0.3$, modulus $30 \times 10^6$ psi, nodes 1 and 2 fixed, and node 3 loaded by $+2000$ lb in $x$ and $-2000$ lb in $y$. Determine nodal displacement and element stress.
Q15 — Point-Loaded Triangular Plate (inches) x y O 1 (0, 10) 2 (0, 0) 3 (20, 10) 2000 lb 2000 lb 20 in 10 in Material Data E = 30 × 10⁶ psi ν = 0.3 t = 0.3 in Plane stress fixed edge

Figure Q15 — Thin triangular plate fixed at nodes 1 & 2, with point loads at node 3

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (in inches):

$$\text{Node 1: }(0,\;10)\;\text{in}, \quad \text{Node 2: }(0,\;0)\;\text{in}, \quad \text{Node 3: }(20,\;10)\;\text{in}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 10 \\ 1 & 0 & 0 \\ 1 & 20 & 10\end{vmatrix}$$ $$= 0(0-10) + 0(10-10) + 20(10-0)$$ $$= 200\;\text{in}^2$$
$$\boxed{A = 100\;\text{in}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 10 = -10$$ $$b_2 = y_3 - y_1 = 10 - 10 = 0$$ $$b_3 = y_1 - y_2 = 10 - 0 = 10$$
$$c_1 = x_3 - x_2 = 20 - 0 = 20$$ $$c_2 = x_1 - x_3 = 0 - 20 = -20$$ $$c_3 = x_2 - x_1 = 0 - 0 = 0$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-0.05 & 0.00 & 0.00 & 0.00 & 0.05 & 0.00 \\ 0.00 & 0.10 & 0.00 & -0.10 & 0.00 & 0.00 \\ 0.10 & -0.05 & -0.10 & 0.00 & 0.00 & 0.05\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}32967032.97 & 9890109.89 & 0.00 \\ 9890109.89 & 32967032.97 & 0.00 \\ 0.00 & 0.00 & 11538461.54\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}5934065.93 & -3214285.71 & -3461538.46 & 1483516.48 & -2472527.47 & 1730769.23 \\ -3214285.71 & 10755494.51 & 1730769.23 & -9890109.89 & 1483516.48 & -865384.62 \\ -3461538.46 & 1730769.23 & 3461538.46 & 0.00 & 0.00 & -1730769.23 \\ 1483516.48 & -9890109.89 & 0.00 & 9890109.89 & -1483516.48 & 0.00 \\ -2472527.47 & 1483516.48 & 0.00 & -1483516.48 & 2472527.47 & 0.00 \\ 1730769.23 & -865384.62 & -1730769.23 & 0.00 & 0.00 & 865384.62\end{bmatrix}$$

(Units: lb/in — stiffness values in the consistent unit system)

Step 7 — Load Vector

Concentrated point loads applied directly at node 3:

$$\{F\} = \begin{Bmatrix} 0\\0\\0\\0\\+2000\\-2000 \end{Bmatrix}\;\text{lb}$$

(No body force or traction — just point loads at node 3.)

Step 8 — Apply Boundary Conditions

Nodes 1 and 2 are fixed (left edge): $u_1 = v_1 = u_2 = v_2 = 0$. Free DOFs: $u_3, v_3$.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}2472527.47 & 0.00 \\ 0.00 & 865384.62\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}2000 \\ -2000\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_3 = 8.0889 \times 10^{-4}\;\text{in}, \qquad v_3 = -0.0023\;\text{in}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 4.0444 \times 10^{-5} \\ 0 \\ -1.1556 \times 10^{-4} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = 1333.33\;\text{lb/in²}, \quad \sigma_y = 400\;\text{lb/in²}, \quad \tau_{xy} = -1333.33\;\text{lb/in²}}$$

Question 16 — Centroid stress in smaller fixed triangular plate

2071 Magh, Q5b, Group B
Steel plate thickness $10$ mm, plane stress, $E = 200$ GPa, $\nu = 0.3$, bottom edge load $5$ kN/m$^2$, geometry $400$ mm by $300$ mm. Determine centroid stress using CST.
Q16 — Fixed-Edge Triangular Plate | 400 × 300 mm x y O 1 (0, 0) 2 (400, 0) 3 (0, 300) 5 kN/m² 400 mm 300 mm Material Data E = 200 GPa ν = 0.3 t = 10 mm Plane stress fixed edge

Figure Q16 — Smaller fixed-edge triangular plate with 5 kN/m² downward traction on edge 1–2

Solution

Step 1 — Node Coordinates & Given Data

Assign coordinates (convert mm to metres for consistent units):

$$\text{Node 1: }(0,\;0)\;\text{m}, \quad \text{Node 2: }(0.4,\;0)\;\text{m}, \quad \text{Node 3: }(0,\;0.3)\;\text{m}$$

Step 2 — Area of the Triangular Element

Using the determinant formula:

$$2A = \begin{vmatrix}1 & 0 & 0 \\ 1 & 0.4 & 0 \\ 1 & 0 & 0.3\end{vmatrix}$$ $$= 0(0-0.3) + 0.4(0.3-0) + 0(0-0)$$ $$= 0.12\;\text{m}^2$$
$$\boxed{A = 0.06\;\text{m}^2}$$

Step 3 — Coefficients $b_i$ and $c_i$

These are computed from cyclic differences of node coordinates:

$$b_1 = y_2 - y_3 = 0 - 0.3 = -0.3$$ $$b_2 = y_3 - y_1 = 0.3 - 0 = 0.3$$ $$b_3 = y_1 - y_2 = 0 - 0 = 0$$
$$c_1 = x_3 - x_2 = 0 - 0.4 = -0.4$$ $$c_2 = x_1 - x_3 = 0 - 0 = 0$$ $$c_3 = x_2 - x_1 = 0.4 - 0 = 0.4$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (must always be zero)

Step 4 — Strain-Displacement Matrix $[B]$

Using the standard CST formula:

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Substituting the values:

$$[B] = \begin{bmatrix}-2.50 & 0.00 & 2.50 & 0.00 & 0.00 & 0.00 \\ 0.00 & -3.33 & 0.00 & 0.00 & 0.00 & 3.33 \\ -3.33 & -2.50 & 0.00 & 2.50 & 3.33 & 0.00\end{bmatrix}$$

Step 5 — Constitutive (Material) Matrix $[D]$

For plane stress:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix}$$
$$[D] = \begin{bmatrix}219780219.78 & 65934065.93 & 0.00 \\ 65934065.93 & 219780219.78 & 0.00 \\ 0.00 & 0.00 & 76923076.92\end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

Using the formula $[K] = t \cdot A \cdot [B]^T [D] [B]$:

$$[K] = \begin{bmatrix}1336996.34 & 714285.71 & -824175.82 & -384615.38 & -512820.51 & -329670.33 \\ 714285.71 & 1753663.00 & -329670.33 & -288461.54 & -384615.38 & -1465201.47 \\ -824175.82 & -329670.33 & 824175.82 & 0.00 & 0.00 & 329670.33 \\ -384615.38 & -288461.54 & 0.00 & 288461.54 & 384615.38 & 0.00 \\ -512820.51 & -384615.38 & 0.00 & 384615.38 & 512820.51 & 0.00 \\ -329670.33 & -1465201.47 & 329670.33 & 0.00 & 0.00 & 1465201.47\end{bmatrix}$$

(Units: force/m — stiffness values in the consistent unit system)

Step 7 — Load Vectors

(a) Body force:

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \begin{Bmatrix}0 \\ -1.5700 \times 10^{-2} \\ 0 \\ -1.5700 \times 10^{-2} \\ 0 \\ -1.5700 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

(b) Traction on edge 1–2 (bottom edge, $l_{12} = 0.4$ m). Uniform $5$ MPa downward:

$$\{f_t\}_{12} = \frac{t \cdot l_{12}}{2} \begin{Bmatrix} 0\\-5000\\0\\-5000\\0\\0 \end{Bmatrix} = \begin{Bmatrix}0 \\ -10 \\ 0 \\ -10 \\ 0 \\ 0\end{Bmatrix}\;\text{kN}$$

(c) Combined:

$$\{F\} = \begin{Bmatrix}0 \\ -10.0157 \\ 0 \\ -10.0157 \\ 0 \\ -1.5700 \times 10^{-2}\end{Bmatrix}\;\text{kN}$$

Step 8 — Apply Boundary Conditions

Left edge fixed: $u_1 = v_1 = u_3 = v_3 = 0$. Free DOFs: $u_2, v_2$.

Step 9 — Solve for Unknown Displacements

Extract the reduced $2 \times 2$ system:

$$[K]_{ff} = \begin{bmatrix}824175.82 & 0.00 \\ 0.00 & 288461.54\end{bmatrix}$$
$$\{F\}_{f} = \begin{Bmatrix}0 \\ -10.0157\end{Bmatrix}$$

Solving $[K]_{ff}\{u\}_f = \{F\}_f$:

$$\boxed{u_2 = 0\;\text{m}, \qquad v_2 = -3.4721 \times 10^{-5}\;\text{m}}$$

Step 10 — Compute Element Strains and Stresses

Strain (constant throughout the element):

$$\{\varepsilon\} = [B]\{u\} = \begin{Bmatrix} \varepsilon_x \\ \varepsilon_y \\ \gamma_{xy} \end{Bmatrix} = \begin{Bmatrix} 0 \\ 0 \\ -8.6803 \times 10^{-5} \end{Bmatrix}$$

Stress (constant throughout the element — and equals centroid stress for CST):

$$\boxed{\sigma_x = 0\;\text{kN/m²}, \quad \sigma_y = 0\;\text{kN/m²}, \quad \tau_{xy} = -6677.13\;\text{kN/m²}}$$
Since CST gives constant strain, the centroid stress equals the element stress above.

Question 17 — CST with usual notation, symbolic setup

2070 Bhadra, Q3b, Group A
Plane stress CST with $E = 30 \times 10^5$ psi, $t = 0.3$ in, $\gamma = 460$ lb/in$^3$, $\nu = 0.3$, and side traction $T = 360$ psi. Determine nodal displacement and stresses using usual notation.
Q17 — CST with Symbolic Setup (usual notation) x y O 1 (x₁, y₁) 2 (x₂, y₂) 3 (x₃, y₃) T = 360 psi (on edge 2–3) γ = 460 lb/in³ Edge 1–2 Edge 1–3 Material Data E = 30 × 10⁵ psi ν = 0.3 t = 0.3 in γ = 460 lb/in³ Plane stress

Figure Q17 — Symbolic CST setup with usual notation, traction T on edge 2–3

Solution

Step 1 — Node Coordinates & Given Data

Let the three nodes have coordinates:

$$\text{Node 1: }(x_1, y_1), \quad \text{Node 2: }(x_2, y_2), \quad \text{Node 3: }(x_3, y_3)$$

Given: $E = 30 \times 10^5$ psi, $\nu = 0.3$, $t = 0.3$ in, $\gamma = 460$ lb/in³, traction $T = 360$ psi on edge 2–3.

Step 2 — Area of the Triangular Element

$$2A = \begin{vmatrix}1 & x_1 & y_1 \\ 1 & x_2 & y_2 \\ 1 & x_3 & y_3\end{vmatrix} = x_1(y_2 - y_3) + x_2(y_3 - y_1) + x_3(y_1 - y_2)$$

Step 3 — Coefficients $b_i$ and $c_i$

$$b_1 = y_2 - y_3, \quad b_2 = y_3 - y_1, \quad b_3 = y_1 - y_2$$ $$c_1 = x_3 - x_2, \quad c_2 = x_1 - x_3, \quad c_3 = x_2 - x_1$$
Check: $b_1+b_2+b_3 = 0$ and $c_1+c_2+c_3 = 0$ (always).

Step 4 — Strain-Displacement Matrix $[B]$

$$[B] = \frac{1}{2A} \begin{bmatrix} b_1 & 0 & b_2 & 0 & b_3 & 0 \\ 0 & c_1 & 0 & c_2 & 0 & c_3 \\ c_1 & b_1 & c_2 & b_2 & c_3 & b_3 \end{bmatrix}$$

Step 5 — Constitutive Matrix $[D]$ (Plane Stress)

$$[D] = \frac{E}{1 - \nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix} = \frac{30 \times 10^5}{1 - 0.09} \begin{bmatrix} 1 & 0.3 & 0 \\ 0.3 & 1 & 0 \\ 0 & 0 & 0.35 \end{bmatrix}$$
$$= 3296703.30 \begin{bmatrix} 1 & 0.3 & 0 \\ 0.3 & 1 & 0 \\ 0 & 0 & 0.35 \end{bmatrix}$$

Step 6 — Element Stiffness Matrix $[K]$

$$[K] = t \cdot A \cdot [B]^T [D] [B]$$

Substitute $t = 0.3$ in, $A$ from Step 2, $[B]$ from Step 4, and $[D]$ from Step 5 to get the $6 \times 6$ stiffness matrix.

Step 7 — Load Vectors

(a) Body force (self-weight in $-y$):

$$\{f_b\} = \frac{\gamma A t}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix} = \frac{460 \times A \times 0.3}{3} \begin{Bmatrix} 0\\-1\\0\\-1\\0\\-1 \end{Bmatrix}$$

(b) Traction on edge 2–3 (length $l_{23}$). Resolve $T = 360$ psi into components $t_x, t_y$ perpendicular to edge 2–3:

$$\{f_t\}_{23} = \frac{t \cdot l_{23}}{2} \begin{Bmatrix} 0\\0\\t_x\\t_y\\t_x\\t_y \end{Bmatrix}$$

(c) Total: $\{F\} = \{f_b\} + \{f_t\}_{23}$.

Step 8 — Apply Boundary Conditions

If nodes 1 and 2 are restrained: $u_1 = v_1 = u_2 = v_2 = 0$.

Partition the system and extract the $2 \times 2$ reduced system for the free DOFs $(u_3, v_3)$.

Step 9 — Solve for Displacements

$$[K]_{ff} \begin{Bmatrix} u_3 \\ v_3 \end{Bmatrix} = \begin{Bmatrix} F_{3x} \\ F_{3y} \end{Bmatrix}$$
$$\begin{Bmatrix} u_3 \\ v_3 \end{Bmatrix} = [K]_{ff}^{-1} \begin{Bmatrix} F_{3x} \\ F_{3y} \end{Bmatrix}$$

Step 10 — Compute Strains and Stresses

$$\{\varepsilon\} = [B]\{u\}, \qquad \{\sigma\} = [D]\{\varepsilon\} = [D][B]\{u\}$$
$$\boxed{\text{Substitute the computed } u_3, v_3 \text{ into the full displacement vector and use the equations above.}}$$
Note: The full geometry dimensions were not clearly visible in the scanned question. Once specific $(x_i, y_i)$ coordinates are known, substitute them into each step above to get numerical answers.
← 4.6 CST Element Notes