Problem 1 - Concepts and Applications of Finite Element Analysis

Question

Let a uniform bar element of axial stiffness \(k=A E / L\) be allowed only \(x\)-direction displacements (along its axis). Its stiffness matrix [k] operates on nodal displacements \(u_{1}\) and \(u_{2}\). Transform [k] so that it operates on nodal d.o.f. \(u_{1}\) and \(u_{r}\), where \(u_{r}\) is the displacement of node 2 relative to node \(1 .\)

Step-by-Step Solution

Answer

The transformed stiffness matrix to operate on \(u_{1}\) and \(u_{r}\) can be obtained by pre-multiplying and post-multiplying the original stiffness matrix by the transformation matrix and its transpose, respectively.

Step 1: Define Constants and Variables

The constants in this problem are the axial stiffness \(k = AE/L\), while variables include nodal displacements \(u_{1}\) and \(u_{2}\). Transform [k] to operate on \(u_{1}\) and \(u_{r}\). hence, \(u_{r}\) is the displacement of node 2 relative to node 1 (\(u_{r} = u_{2} - u_{1}\))

Step 2: Original Stiffness Matrix

The original stiffness matrix [k] that operates on \(u_{1}\) and \(u_{2}\) can be represented as: \[k] = \(k * [ [1, -1], [-1, 1] ]\)

Step 3: Transformation Matrix

We need a transformation matrix to transform the original stiffness matrix to operate on \(u_{1}\) and \(u_{r}\). The transformation matrix [T] is given by: \[T] = \( [ [1, 0], [-1, 1] ] \)

Step 4: Transform the Stiffness Matrix

We can now transform the original stiffness matrix [k] by pre-multiplying and post-multiplying it by the transformation matrix [T] and its transpose. This yields the transformed stiffness matrix: \[k'] = \ (T^T * [k] * T)\

Key Concepts

Stiffness Matrix Transformation

A foundational concept in finite element analysis is the stiffness matrix transformation. It's essential when the nodal degrees of freedom (d.o.f) do not align with the global coordinates, or when we wish to express the stiffness relating to different d.o.f. In our exercise, we need to make the transformation so that the stiffness matrix, originally operating on nodal displacements along the x-axis, can operate on one original d.o.f ((u_1)) and a relative displacement ((u_r)).

The transformation matrix ((T]), acts as a bridge, converting the influence of one set of displacements ((u_1) and (u_2)) to another ((u_1) and (u_r)). The process involves two operations using the transformation matrix: pre-multiplication and post-multiplication of the original stiffness matrix. The resulting matrix expresses the same physical properties but in the context of the new nodal dof framework, streamlining the analysis for systems where relative displacements are more meaningful or intuitive.

Uniform Bar Element

A uniform bar element is a basic structural element within finite element analysis characterized by a constant cross-sectional area (A), uniform material properties such as Young's modulus (E), and a constant length (L). Its behavior is one-dimensional with stress and strain only occurring along its length, making it a perfect subject for studying axial deformation. Axial stiffness is a key property, derived from the material and geometric properties, dictating how the bar will deform under axial forces.

In practical terms, it means that the uniform bar will only allow for stretching or compression along its length, simplifying the calculations and aiding in building an intuitive understanding of how the structure responds to forces.

Axial Stiffness

Axial stiffness ((k)) of a uniform bar element is defined as the ratio of axial force to the displacement along the axis ((F/u)). It is directly calculated from the bar's cross-sectional area (A), its length (L), and the material's Young modulus (E). The formula (k = AE/L) succinctly captures the relationship between these factors.

High axial stiffness indicates a material's resistance to being deformed axially. In our exercise, this axial stiffness underpins the kinematic relationships captured within the stiffness matrix, allowing us to evaluate the bar's reaction to applied nodal displacements. Understanding axial stiffness is crucial for designing and analyzing structures, as it helps predict their behavior under load.

Nodal Displacements

Nodal displacements refer to the movement of nodes, points of interest within the finite element mesh, in response to forces or other physical effects. These displacements are a pivotal part of finite element analysis given they serve as key inputs and outputs of the model. From the calculated displacements, stresses and strains can be deduced, revealing how a structure might behave under certain conditions.

In the context of our example, we're concerned with (u_1) and (u_2), which represent the displacements of two nodes along the x-axis of the bar. We then define a relative displacement (u_r) as the difference between (u_2) and (u_1), which simplifies the analysis of the element's deformation. Monitoring nodal displacements is crucial for ensuring a structure's integrity, aiding in predicting potential failure points and the overall structural performance.