Problem 17 - Concepts and Applications of Finite Element Analysis

Question

Consider a straight, uniform shaft of solid circular cross section, with a node at each end. (a) Let nodal d.o.f. be angular rotation vectors parallel to the bar, one at each end. Nodal loads are axially directed torque vectors. What is \([\mathrm{k}]\), in terms of the length, shear modulus, and radius of the cross section? (b) Let the bar be inclined at angle \(\beta\) to the \(x\) axis, with \(\theta_{x}\) and \(\theta_{y}\) as d.o.f. at each node (rotations about the \(x\) and \(y\) coordinate axes). What is \([\mathrm{k}]\) ? As in part (a), consider torsional stiffness only.

Step-by-Step Solution

Answer

The stiffness [k] in part (a) can be found by using the formula for the stiffness of a shaft, while considering the given parameters. In part (b), the stiffness [k] would need to be calculated after transformation to the global coordinates, taking into account the incline and new d.o.fs.

Step 1 - Calculate [k] in Part (a)

Stiffness of a shaft is given by the formula \[ \mathrm{k} = \frac{\pi G r^4}{2 L} \] where, G is shear modulus, r is radius, and L is length of the shaft. Substitute the values of G, r, L in above formula to get [k] for part (a).

Step 2 - Determine the new d.o.fs. in Part (b)

The shaft is now inclined at an angle \(\beta\). The rotation d.o.fs. \(\theta_{x}\) and \(\theta_{y}\) are given, instead of \(\theta\). Write the relations for \(\theta_{x}\) and \(\theta_{y}\) in terms of \(\theta\) and \(\beta\). Create the transformation matrix T based on these relations. It will be a diagonal matrix with cosine values of the angle \(\beta\).

Step 3 - Calculate [k] in Part (b)

Transform the stiffness matrix from part (a) using transformation matrix T. The transformed stiffness matrix in global coordinates is given by the expression \[ \mathrm{k}' = T^T \mathrm{k} T \] In the above expression, replace the values of T and [k], and calculate [k'] giving the stiffness in global coordinates for part (b).

Key Concepts

Torsional Stiffness

When dealing with structures like shafts in engineering, understanding torsional stiffness is fundamental. Torsional stiffness is essentially the resistance a shaft offers to twisting. It's a crucial property when determining how much a material will deform under a torsional load.
Torsional stiffness depends on several factors:

• Shear Modulus (G): This is a material property that describes its rigidity. A higher shear modulus means the material is less likely to deform under shear stress.
• Radius (r) of the cross-section: Since the resistance to twist increases with the fourth power of the radius, even small changes in radius can significantly affect torsional stiffness.
• Length (L) of the shaft: Longer shafts are more flexible under torsional loads.

The formula for calculating torsional stiffness, as given by \[ \mathrm{k} = \frac{\pi G r^4}{2 L} \] shows how these factors interact. Given this relationship, engineers can design shafts that are adequately stiff by selecting appropriate materials and dimensions.

Shear Modulus

The shear modulus (G) is a fundamental mechanical property that reflects how a material responds to shear stress. This modulus is critical in understanding materials used in construction and mechanical components.
The shear modulus is defined as the ratio of shear stress to the shear strain in a material:\[G = \frac{\text{shear stress}}{\text{shear strain}}\] Key points about shear modulus include:

• Materials with a high shear modulus are more rigid and less likely to deform when subjected to shear forces.
• Understanding shear modulus helps in predicting how much a material will twist under a given torque.
• It is proportional to the torsional stiffness of a shaft, as indicated in the stiffness equation \(k = \frac{\pi G r^4}{2 L}\).

Shear modulus is a vital parameter for calculating torsional stiffness in finite element analysis, allowing engineers to understand better and design systems subjected to torsional stress.

Transformation Matrix

In part (b) of engineering problems, often, geometric considerations like angles play a role in transforming coordinate systems. A transformation matrix helps convert elements, like stiffness matrices, from a local to a global coordinate system.
Transformation matrices are valuable when:

• The element is oriented at an angle with respect to a primary axis (like the shaft inclined at angle \(\beta\) to the \(x\) axis).
• We need to connect the rotations \(\theta_{x}\) and \(\theta_{y}\) to the principal rotation \(\theta\).
• Ensuring calculations are aligned with the physical orientation of components in the real world.

The transformation matrix \(T\) for aligning the rotations, especially when angles of inclination like \(\beta\) are involved, uses trigonometric relationships. Typically, for a simple angle of inclination:\[ T = \begin{bmatrix} \cos(\beta) & 0 \ 0 & \cos(\beta) \end{bmatrix} \]This diagonal form reflects the simplicity when aligned aspects correspond directly to global x and y axes. Applying this transformation ultimately enables the translation of local stiffness to global coordinates using:\[ \mathrm{k}' = T^T \mathrm{k} T \]This transformation ensures that analysis adheres to the realistic conditions under which structures operate.