Problem 1 - Concepts and Applications of Finite Element Analysis

Question

For a straight bar on the \(x\) axis that may stretch and bend in the \(x y\) plane, strain and strain energy for small displacements are given by $$ \varepsilon_{x}=u_{x_{x}}+\frac{1}{2} v_{x}^{2}-y v_{x x} \quad U=\int \frac{1}{2} E \varepsilon_{x}^{2} d V $$ Show how these expressions yield familiar conventional element stiffness matrices and a stress stiffness matrix.

Step-by-Step Solution

Answer

The derivation process ultimately yields the conventional element stiffness matrices and a stress stiffness matrix from the given expressions for strain and strain energy, which are applicable to small displacements in the \(x y\) plane of a straight bar along the \(x\) axis.

Step 1: Identify the Given Strain and Strain Energy Equations

The exercise provides the equation for strain \(\varepsilon_{x}\) as \(u_{x_{x}}+\frac{1}{2} v_{x}^{2}-y v_{x x}\) and the equation for strain energy \(U\) as \( \int \frac{1}{2} E \varepsilon_{x}^{2} d V\). It should be noted that these equations are non-linear. This suggests that the derivation process may involve linearization.

Step 2: Linearize the Strain Equation

Since we are dealing with small displacements, we can linearize the strain equation by ignoring the non-linear terms. Therefore, the strain equation becomes \(\varepsilon_{x}=u_{x_{x}}-y v_{x x}\).

Step 3: Derive the Strain Energy Integral

Substitute the linearized strain equation into the strain energy equation to obtain an expression for strain energy \(U\) in terms of the displacement gradients \(u_{x_{x}}\) and \(v_{x x}\). This results in the integral \(U=\int \frac{1}{2} E (u_{x_{x}}-y v_{x x})^{2} d V\).

Step 4: Develop the Element Stiffness Matrices

The element stiffness matrices can be obtained by taking the derivative of the strain energy integral with respect to the displacement gradients and setting it equal to zero, as this minimizes the strain energy. This process will yield two element stiffness matrices, one for the displacement \(u\) and the other for the displacement \(v\).

Step 5: Derive the Stress Stiffness Matrix

Similarly, the stress stiffness matrix can be derived by differentiating the strain energy integral with respect to the displacements and then setting the derivative equal to zero. This will result in a stress stiffness matrix that relates the stresses to the small displacements in the plane.

Key Concepts

Strain Energy in Finite Elements

Understanding the concept of strain energy in finite elements is crucial for engineering applications that deal with structural analysis. Strain energy refers to the energy stored within a material due to deformation caused by applied forces or displacements. In the context of finite element analysis (FEA), we examine how this strain energy is calculated and distributed throughout a structure.

According to the provided textbook exercise, for a straight bar that may stretch and bend in the xy plane, the strain energy (\(U\)) is expressed as an integral: \[U = \int \frac{1}{2} E \varepsilon_{x}^{2} d V\] where \(E\) represents the material's modulus of elasticity, and \(\varepsilon_{x}\) is the strain component in the x-direction. The equation highlights that strain energy is directly proportional to the square of the strain, indicating that even small strains can lead to significant amounts of energy being stored in the material, especially if the material is very stiff (i.e., has a high modulus of elasticity).

To help students better understand, it's helpful to visualize the bar as a spring. When a force stretches or compresses it, the bar stores energy just like a spring does when it's stretched or compressed. The strain energy is an integral measure of all the 'spring-like' responses throughout the element volume. Imagine pushing a soft mattress versus a stiff one; the mattress (or material) resistance and the deformation pattern directly affect how the strain energy is distributed. This analogy helps simplify the concept of strain energy for better comprehension.

Element Stiffness Matrices

Moving on to the concept of element stiffness matrices, these are fundamental to finite element analysis. The stiffness matrix of an element measures how resistant it is to deformation when subject to loading. It mathematically characterizes the relationship between nodal forces and the corresponding displacements.

In the step-by-step solution, once the strain equation is linearized, we derive the element stiffness matrices by setting the partial derivative of the strain energy with respect to displacement gradients to zero. This is done because, physically, a stable system will 'seek' a state where potential energy (in this case, strain energy) is minimized. It is akin to letting go of the mattress we discussed earlier; after displacing it, it returns to its original shape, representing a state of minimal energy.

The process involves calibration of the material's response to various types of loads and binds each element's reaction to deformation within the larger structure. The matrices contain values that correlate to the bar's ability to resist bending and stretching, which can be visualized as the sturdiness and rigidity of individual parts of a framework. Here's how we can make this concept digestible: think of building a model bridge with different materials for each section - some might be made of strong metal, while others of balsa wood. Each piece's resistance to bending under weight can be described by its own stiffness matrix and impacts how the entire bridge would respond to a load.

Stress Stiffness Matrix

Finally, the stress stiffness matrix is another critical component in FEA. While the element stiffness matrices relate to how the structure deforms under load, the stress stiffness matrix associates the induced stresses to the displacements in the material. It essentially captures the additional stiffness provided by the pre-stress state of the structure.

The textbook solution reveals that this matrix is also derived from the strain energy integral. By differentiating the strain energy with respect to the displacements and setting it to zero, we obtain a relationship that reflects how stresses within the bar relate to its small movements or deformations.

To paint a clearer picture, imagine a tightly stretched rubber band. Even before you apply any new forces, the pre-stretched rubber band already has internal stresses resisting deformation. This pre-loaded state lends the rubber band a certain stiffness on its own. Similarly, in engineering structures, initial stresses (such as those due to pre-tensioning in beams) can significantly influence how the structure will respond to additional loading and deformations. Therefore, the stress stiffness matrix isn’t just academic exercise; it's a representation of a real-world phenomenon that's vital for predicting the performance of pre-stressed structures and components.