JNTU B.Tech IV FINITE ELEMENT METHODS FOR CIVIL ENGINEERING Previous Question Paper (July/August - 2022 - R18)

The Finite Element Method (FEM) is a numerical technique for finding approximate solutions to boundary value problems for partial differential equations. The core concept involves discretizing a complex physical domain into a collection of smaller, simpler sub-domains called finite elements. Within each element, the unknown field variable (such as displacement, temperature, or pressure) is approximated by simple interpolation functions, typically polynomials, defined at specific points called nodes.

The formulation process generally follows these steps:

1. Discretization: Dividing the continuum into a finite element mesh.
2. Selection of Interpolation Functions: Choosing shape functions to represent the variation of the field variable within the element.
3. Element Characteristic Equations: Deriving the element stiffness matrix and load vector using variational principles (e.g., Principle of Minimum Potential Energy) or the Method of Weighted Residuals (e.g., Galerkin method).
4. Assembly: Combining individual element equations into a global system of equations, typically represented as [K]{u} = {F}, where [K] is the global stiffness matrix, {u} is the nodal displacement vector, and {F} is the global force vector.
5. Application of Boundary Conditions: Incorporating known constraints to solve the system.
6. Solution and Post-processing: Solving the algebraic equations for nodal values and deriving secondary quantities like stresses or strains.

FEM has emerged as a powerful engineering tool due to its versatility and robustness. Unlike analytical methods limited to simple geometries and boundary conditions, FEM can handle complex, irregular geometries, non-homogeneous material properties, and diverse loading conditions. The rapid advancement in computational power has allowed engineers to simulate large-scale, non-linear structural behaviors, fluid-structure interactions, and thermal analyses with high precision. By reducing the reliance on expensive physical prototyping and enabling iterative design optimization, FEM has become indispensable in modern structural, geotechnical, and mechanical engineering.

To solve for a simply supported beam of length L, flexural rigidity EI, under a uniformly distributed load w, we use the Rayleigh-Ritz method. The potential energy functional Π is given by:

We assume a deflection function v(x) that satisfies the essential boundary conditions (v=0 at x=0 and x=L). Let v(x) = a sin(πx/L). Calculating the derivatives:

Substituting into the energy functional:

Evaluating the integrals over 0 to L:

To minimize Π, set dΠ/da = 0:

The mid-span deflection is v(L/2) = a sin(π/2) = 4wL⁴ / (EIπ⁵) ≈ 0.01307 wL⁴/EI. The exact solution is 5wL⁴ / 384EI ≈ 0.01302 wL⁴/EI. The bending moment M = -EI(d²v/dx²) at mid-span is:

The exact mid-span bending moment is wL²/8 = 0.125 wL². The Rayleigh-Ritz method provides a close approximation with a single-term trigonometric function.