Problem 4 - Concepts and Applications of Finite Element Analysis

Question

The model shown is a square grillage of uniformly spaced beam elements that lie in the \(x y\) plane and are welded together at nodes. Assume that all elements are identical and that nodal d.o.f. are lateral (z direction) displacements \(w_{i}\) and rotations \(\theta_{d i}\) and \(\theta_{y i}\). Supports impose \(w_{i}=0\) at all nodes \(i\) on the square boundary of the grillage. Internal nodes are not restrained. In each part of this problem, lateral (z direction) forces of magnitude \(P\) act on intemal nodes as described below. State what portion of the grillage constitutes the smallest acceptable model, and what its boundary conditions must be, if: (a) Each node carries a load \(P\), all acting in the same direction. (b) Loads \(P\) act upward for \(y>0\), act downward for \(y<0\), and are omitted on \(y=0\) (c) Loads \(P\) alternate in direction by quadrant. That is, they act upward in the first and third quadrants, downward in the second and fourth quadrants, and are omitted on \(x\) and \(y\) axes. (d) Loads \(P\) alternate in direction by octant.

Step-by-Step Solution

Answer

For case (a), the smallest model is one beam; for case (b), it's half the structure along the \(y\) axis; for case (c), it's one-fourth of the grillage (a single quadrant); for case (d), it's one-eighth of the grillage. The boundary conditions maintain \(\theta_{d i} = 0\) and \(\theta_{y i} = 0\) at the boundary beams.

Step 1: Identifying Symmetry for Case (a)

In this case, the load \(P\) acts in the same direction at each node. The symmetry is uniform across the entire structure. Therefore, the smallest acceptable model will be a single beam with boundary conditions \(\theta_{d i} = 0\) and \(\theta_{y i} = 0\) at its ends.

Step 2: Deciphering Load Distribution for Case (b)

Here, \(P\) acts upward for \(y > 0\) and downward for \(y < 0\), and is not applied on \(y = 0\). This creates a mirror symmetry across the \(y\) axis of the grillage. Therefore, the minimal model here would be half of the structure (along the \(y\) axis) with boundary conditions \(\theta_{d i} = 0\) and \(\theta_{y i} = 0\) at \(y = 0\) and at the outermost beam.

Step 3: Understanding Alternating Load Directions for Case (c)

In this scenario, \(P\) acts upward in first and third quadrants, downward in the second and fourth quadrants, and is omitted on \(x\) and \(y\) axes. There's a rotational symmetry about the origin here. The smallest acceptable model is one-fourth of the grillage (a single quadrant). The boundary conditions remain \(\theta_{d i} = 0\) and \(\theta_{y i} = 0\) at the axes of this quadrant.

Step 4: Analyzing Octant Direction Variation for Case (d)

In the final case, the load \(P\) alternates by octant. The smallest model here would be one-eighth of the grillage, maintaining boundary conditions \(\theta_{d i} = 0\) and \(\theta_{y i} = 0\) at the boundary beams.

Key Concepts

Beam Elements

Beam elements are the fundamental components in finite element analysis used to model structures like beams, frames, and grillages. These elements capture the behavior of beams under various loading conditions, accounting for bending and twisting.

Key characteristics of beam elements include:

• **Geometry:** Defined by nodes at each end that specify the position in a structure.
• **DoF (Degrees of Freedom):** Each node typically has displacements and rotations, e.g., lateral displacement ( w_{i} ) and rotations ( heta_{di} and heta_{yi} ).
• **Uniformity:** In the given model, all beam elements are identical, simplifying the analysis.

Understanding beam elements is crucial for modeling complex structures like the grillage in question.

Nodal Displacements

Nodal displacements refer to how nodes move or rotate under applied loads. In the beam grillage model, nodal displacements are observed in the lateral ( z direction).

Important aspects of nodal displacements include:

• **Lateral Movements:** Displacements in the z direction are denoted as w_{i} , indicating vertical shifts.
• **Rotational Movements:** Rotations about the x and y axes are specified as heta_{di} and heta_{yi} , showing twisting at a node.

Properly calculating nodal displacements helps predict how the structure behaves under various loads.

Boundary Conditions

Boundary conditions are constraints applied to the nodes of a structure to simulate real-world supports or connections. For the grillage model, these are critical in determining the behavior of the system.

Here’s what to note about boundary conditions:

• **Supports:** Imposed at nodes on the boundary, keeping w_{i} = 0 , meaning the grillage's outer edge doesn't move laterally.
• **Symmetry Conditions:** For cases like rotational symmetry, particular nodes have rotations constrained ( heta_{d i} = 0 and heta_{y i} = 0 ), simplifying the model.

Applying the correct boundary conditions ensures that the model accurately reflects the real-world structure.

Structural Symmetry

Structural symmetry refers to using symmetrical patterns to reduce the complexity of a model in finite element analysis. It allows the use of smaller models by exploiting the symmetry in load distribution.

In the exercise, symmetry simplifies each case:

• **Uniform Loads (a):** Full symmetry means the model can be reduced to a single beam.
• **Mirror (b):** Mirror symmetry along the y -axis halves the model.
• **Rotational (c):** Rotational symmetry about the origin allows using a quadrant of the structure.
• **Octant Variation (d):** The model is further reduced to an eighth of the structure.

Recognizing and applying structural symmetry is essential for efficient modeling.