Mechanical Engineering Fundamentals for Mechatronics

Chapter 4
Mechanics and Strength of Materials

Bachelor Program in Mechatronics Engineering
Mechanical Engineering Fundamentals for Mechatronics   •   Bachelor

Mechanics and Strength of Materials

Introduction

Mechanics and strength of materials provide the analytical foundation for designing safe mechanical, robotic, and mechatronic systems. A mechatronic product may contain motors, shafts, gears, brackets, frames, pressure vessels, sensors, and actuators. Each component transmits force or torque and must resist deformation and failure while remaining sufficiently accurate for its function.

This chapter develops the principal models used to evaluate stress, strain, elastic behavior, thermal stress, torsion, bending, shear, combined loading, stress concentration, beam deflection, and column buckling. These models are idealizations. Their value depends on using them within their assumptions and comparing predictions with testing, standards, and, where useful, finite-element analysis (FEA).

Learning Objectives

By the end of this chapter, a student should be able to:

  1. Distinguish force, stress, strain, stiffness, strength, and failure.
  2. Calculate normal, shear, bending, torsional, and bearing stresses.
  3. Apply Hooke’s law and generalized elastic relationships.
  4. Evaluate deformation caused by mechanical and thermal loading.
  5. Analyze shafts, beams, combined loading, and stress concentrations.
  6. Estimate beam deflection and column buckling loads.
  7. Select suitable analytical, experimental, and finite-element methods.
  8. Recognize modeling assumptions, uncertainty, and safe engineering practice.

1. Loads, Equilibrium, and Internal Actions

A load is an externally applied force, moment, pressure, or temperature change. Before calculating stress, isolate the component with a free-body diagram. For static equilibrium in Cartesian coordinates,

∑Fx = 0, ∑Fy = 0, ∑Fz = 0, ∑M = 0

The internal actions exposed by a section cut are axial force N, shear forces V, bending moments M, and torque T. Their signs depend on the chosen convention. A correct free-body diagram and consistent sign convention are more important than memorizing isolated formulas.

Typical loading modes are:

Loading modePrincipal internal actionTypical component
Tension or compressionAxial force NTie rod, actuator link
Direct shearShear force VPin or bolt
BendingMoment MRobot arm or bracket
TorsionTorque TMotor shaft
Stability loadingCompressive force PColumn or slender link

2. Stress and Strain

Stress describes internal force intensity. Average normal stress over an area A is

σ = P/A

where σ is normal stress in pascals (Pa) or megapascals (MPa), and P is an axial force in newtons (N). The expression assumes a prismatic member, a centrally applied load, and a sufficiently uniform stress distribution away from load introduction and geometric discontinuities.

Average direct shear stress is

τavg = V/A

where τ_avg is shear stress and V is transverse force. In beams, the actual shear stress is generally nonuniform and is calculated from

τ = VQ/(It)

Here V is shear force, Q is the first moment of the area outside the point of interest about the neutral axis, I is the second moment of area about the neutral axis, and t is the local width measured parallel to the shear direction. For a rectangular section, the maximum shear stress is 1.5 times the average value; for a circular section under transverse shear, the distribution is also nonuniform.

Normal strain is change in length divided by original length:

ε = ΔL/L₀

Strain is dimensionless, although it is often reported in microstrain (με), where 1 με = 10⁻⁶.

Shear strain γ is the angular distortion of a material element and is measured in radians, which are dimensionless in SI units. Tensile stress and strain are normally positive; compressive values are negative under the usual convention.

Worked Example 1: Axial Stress and Strain

A steel tie rod of diameter 12 mm carries a tensile load of 18 kN. Take E = 200 GPa. Determine its average stress and elastic strain.

The cross-sectional area is

A = πd²/4 = π(0.012 m)²/4 = 1.131 × 10⁻⁴ m²

Therefore,

σ = P/A = 18,000/(1.131 × 10⁻⁴) = 159.2 MPa

Using Hooke’s law,

ε = σ/E = (159.2 × 10⁶)/(200 × 10⁹) = 7.96 × 10⁻⁴

The rod experiences approximately 159 MPa tensile stress and 796 με elastic strain. This result is valid only if the material remains elastic and the load is sufficiently concentric.

3. Elastic Behavior and Material Properties

The engineering stress–strain curve identifies several important regions. In the elastic region, unloading returns the specimen approximately to its original dimensions. The proportional limit is the end of an approximately linear relationship. The yield strength marks the onset of significant plastic deformation, while the ultimate tensile strength is the maximum engineering stress. Fracture occurs later, often after necking in a ductile metal.

For uniaxial linear elasticity,

σ = Eε

where E is Young’s modulus. E measures stiffness, not strength. A material can be very stiff but brittle, or relatively compliant but tough.

For isotropic three-dimensional elasticity, the normal strain in the x direction is

εx = [σx − ν(σy + σz)]/E

where ν is Poisson’s ratio. The corresponding relations apply in the other directions. For pure shear in the elastic range,

τ = Gγ

where G is the shear modulus. For an isotropic material,

G = E/[2(1 + ν)]

These equations assume small strains, linear elasticity, isotropic material behavior, and reversible deformation. Composites, polymers, biological materials, and additively manufactured parts may require anisotropic, viscoelastic, plastic, or nonlinear models.

The total elongation of a uniform axial member is

δ = PL/(AE)

For a variable member or variable load, the more general small-deformation expression is

δ = ∫[N(x)/(A(x)E(x))] dx

4. Thermal Stress and Strain

A freely expanding uniform member subjected to a temperature change ΔT undergoes thermal strain

εT = αΔT

where α is the coefficient of thermal expansion in K⁻¹. If expansion is completely prevented, the induced stress is

σT = −EαΔT

for a uniform, linear-elastic, uniaxial member. Heating a restrained member therefore produces compression; cooling produces tension. Partial restraint requires compatibility and stiffness analysis rather than the fully restrained formula.

Worked Example 2: Fully Restrained Thermal Member

An aluminum bar has E = 70 GPa and α = 23 × 10⁻⁶ K⁻¹. It is fully restrained and heated by 40 °C. The thermal stress is

σT = −EαΔT = −(70 × 10⁹)(23 × 10⁻⁶)(40) = −64.4 MPa

The stress is 64.4 MPa in compression. This calculation neglects plasticity, temperature-dependent properties, and restraint flexibility.

Thermal loading is especially important in precision mechanisms, motor housings, brake systems, welded structures, and assemblies made from dissimilar materials. Design may use expansion joints, compliant mounts, matched materials, or active thermal compensation.

5. Torsion of Circular Shafts

For a circular shaft in Saint-Venant elastic torsion,

τ(r) = Tr/J

where T is applied torque, r is radial distance from the shaft center, and J is the polar second moment of area. The maximum stress occurs at the outer radius c:

τmax = Tc/J

For a solid circular shaft,

J = πd⁴/32

For a hollow circular shaft,

J = π(D⁴ − d⁴)/32

where D and d are outside and inside diameters. The angle of twist over length L is

φ = TL/(JG)

with φ in radians. These equations assume a circular, prismatic shaft, elastic material, uniform torque, and Saint-Venant torsion without severe warping or stress concentration.

Worked Example 3: Motor Shaft

A solid steel shaft of diameter 20 mm transmits 120 N·m over 0.5 m. Take G = 80 GPa. Since

J = π(0.020)⁴/32 = 1.571 × 10⁻⁷ m⁴

then

τmax = (120)(0.010)/(1.571 × 10⁻⁷) = 7.64 × 10⁶ Pa = 7.64 MPa

and

φ = (120)(0.5)/[(1.571 × 10⁻⁷)(80 × 10⁹)] = 0.0478 rad = 2.74°

A shaft design must check both strength and torsional stiffness. Excessive twist can cause encoder errors or gear misalignment even when stress is below the material limit.

6. Bending and Shear in Beams

For a slender beam in elastic bending, the flexure formula is

σx = −My/I

where M is the internal bending moment, y is the distance from the neutral axis, and I is the second moment of area. The maximum magnitude occurs at the farthest fiber, y = c:

σmax = Mc/I = M/S

where S = I/c is the elastic section modulus.

For common sections:

Irect = bh³/12

for a rectangle of width b and height h, and

Icircle = πd⁴/64

for a solid circular section. The strong dependence on h³ or d⁴ explains why depth and diameter are highly effective ways to increase bending resistance.

The relationship between distributed load w(x), shear V(x), and bending moment M(x) is

dV/dx = −w(x), dM/dx = V(x)

subject to the selected sign convention. In a simply supported beam with a central point load P and span L, the maximum bending moment is PL/4. In a simply supported beam with uniform load w over the entire span, it is wL²/8.

Beam shear is important near supports and load applications. The shear formula τ = VQ/(It) is preferable to V/A for nonuniform distributions. Thin webs may be governed by shear stress or shear buckling even when the flanges govern bending strength.

7. Beam Deflection

Strength alone does not ensure functional performance. Deflection can affect positioning accuracy, vibration, sealing, and fatigue life. For a prismatic Euler–Bernoulli beam,

d²v/dx² = M(x)/(EI)

with sign depending on convention. E is Young’s modulus and I is the second moment of area. This model assumes small deflection, slender beams, linear elasticity, and negligible shear deformation.

Useful cases include:

Beam and loadingMaximum deflection
Simply supported, central point load PPL³/(48EI)
Cantilever, end load PPL³/(3EI)
Simply supported, full uniform load w5wL⁴/(384EI)
Cantilever, full uniform load wwL⁴/(8EI)

Worked Example 4: Cantilever Positioning Arm

A cantilever arm has L = 0.4 m, carries an end load of 100 N, and has E = 70 GPa and I = 2.0 × 10⁻⁸ m⁴. Its end deflection is

δ = PL³/(3EI) = (100)(0.4)³/[3(70 × 10⁹)(2.0 × 10⁻⁸)] = 0.0152 m

The deflection is 15.2 mm. If the mechanism permits only 1 mm, increasing strength alone is not sufficient; the section geometry, span, material, or support arrangement must be changed. Because deflection varies with L³, reducing the unsupported length is especially effective.

For short, deep beams, Timoshenko beam theory includes shear deformation. For large deflections, nonlinear geometry must be considered. Boundary conditions must represent the actual joint stiffness; an assumed perfectly fixed support can substantially underpredict displacement.

8. Combined Loading and Principal Stress

Components commonly experience axial force, bending, torsion, and shear simultaneously. At a point, the normal stress from axial force and bending may be written

σ = P/A ± Mc/I

The sign depends on whether axial and bending stresses reinforce or oppose one another. Torsion adds shear stress τ. For a plane stress state with σ_x, σ_y, and τ_xy, the principal stresses are

σ₁,₂ = (σx + σy)/2 ± √{[(σx − σy)/2]² + τxy²}

The maximum in-plane shear stress is

τmax = √{[(σx − σy)/2]² + τxy²}

For a common shaft case with σ_y = 0,

σ₁,₂ = σx/2 ± √{(σx/2)² + τxy²}

For ductile metals, the von Mises equivalent stress under plane stress is

σvm = √(σx² − σxσy + σy² + 3τxy²)

A preliminary yield check compares σ_vm with the yield strength, with a factor of safety selected according to uncertainty, consequence of failure, loading variation, and applicable standards. This is not a complete fracture, fatigue, creep, or stability assessment.

9. Stress Concentration

Abrupt changes in geometry—holes, keyways, grooves, threads, sharp corners, and shoulders—disturb the stress field. The theoretical stress concentration factor is

Kt = σmax,local/σnominal

where σ_nominal is calculated using a specified nominal area or section convention. K_t depends on geometry and loading mode, not directly on material strength.

In fatigue design, the fatigue notch factor K_f is often lower than K_t because material sensitivity and notch-root plasticity affect fatigue behavior. A commonly used relationship is

Kf = 1 + q(Kt − 1)

where q is notch sensitivity, 0 ≤ q ≤ 1. This is an empirical design model and must be used with appropriate material and notch data.

Good practice includes generous fillet radii, smooth transitions, avoidance of machining marks at high-stress locations, correct thread reliefs, and inspection of critical surfaces. Linear-elastic FEA can identify stress concentration, but a sharp idealized corner may produce a mathematically unbounded peak; the result must then be interpreted using a realistic radius or stress-averaging method.

10. Column Buckling

A slender compression member may fail by instability before its material reaches compressive yield. For an ideal pin-ended Euler column,

Pcr = π²EI/(KL)²

where P_cr is the critical load, L is the unsupported length, E is Young’s modulus, I is the least second moment of area, and K is the effective-length factor. For ideal pin-ended, fixed-fixed, fixed-pinned, and fixed-free conditions, K is approximately 1.0, 0.5, 0.7, and 2.0 respectively, though actual joint behavior may differ.

The radius of gyration is

r = √(I/A)

and slenderness ratio is

λ = KL/r

Euler theory applies to sufficiently slender, straight, elastic columns with idealized loading and supports. Short and intermediate columns require inelastic or empirical relations, such as a code-specified column curve. Initial crookedness, eccentricity, residual stress, and imperfections reduce practical capacity.

11. Testing and Finite-Element Analysis

Analytical equations should be treated as models to be validated, not as substitutes for engineering judgment. A tensile test supplies an elastic modulus, yield strength, ultimate strength, and ductility. Strain gauges measure local strain; load cells measure force; torque transducers measure torque; displacement sensors measure deflection. Calibration, alignment, sampling rate, temperature, and uncertainty must be documented.

A responsible FEA workflow is:

  1. Define the design question and failure criterion.
  2. Idealize geometry while retaining load paths and critical features.
  3. Select material behavior and units consistently.
  4. Apply realistic constraints and loads.
  5. Mesh the model, refining areas with gradients or concentrations.
  6. Check reactions, equilibrium, deformation shape, and convergence.
  7. Compare selected results with hand calculations or test data.
  8. Perform sensitivity studies and report assumptions.

FEA does not automatically correct incorrect boundary conditions. A visually detailed model can still be physically wrong. Singular stresses at point loads, perfectly sharp corners, and fully constrained surfaces require special interpretation. For mechatronic systems, coupled thermal, structural, contact, modal, fatigue, and nonlinear analyses may be needed, but complexity should be added only when justified by the design question.

Illustrative Case: Aluminum Robotic Arm Bracket

Background

Consider a fictional preliminary design of an aluminum bracket supporting a robotic actuator. The bracket behaves approximately as a cantilever of length 0.25 m, with a 600 N end load. Its rectangular cross-section is 40 mm wide and 20 mm deep. Take E = 70 GPa and an allowable normal stress of 120 MPa for this preliminary check.

Problem and Analysis

The end moment is

M = PL = (600)(0.25) = 150 N·m

The second moment of area, with the 20 mm dimension vertical, is

I = bh³/12 = (0.040)(0.020)³/12 = 2.667 × 10⁻⁸ m⁴

The extreme-fiber distance is c = 0.010 m. Thus,

σmax = Mc/I = (150)(0.010)/(2.667 × 10⁻⁸) = 56.3 MPa

The cantilever end deflection is

δ = PL³/(3EI) = (600)(0.25)³/[3(70 × 10⁹)(2.667 × 10⁻⁸)] = 0.0167 m

The predicted bending stress is below the stated allowable value, but the deflection is 16.7 mm. If the actuator requires positioning accuracy within 2 mm, the design fails its serviceability requirement.

Approach and Solution

Increasing the vertical depth from 20 mm to 30 mm changes I in proportion to h³, giving a factor of (30/20)³ = 3.375 increase. The estimated deflection becomes approximately 4.94 mm, still above the target. Shortening the arm, adding a rib, changing the support arrangement, or using a boxed section may be more effective than simply selecting a stronger alloy. A refined FEA model should include bolt holes, fillets, contact regions, and the actual actuator load path. A physical load test should then verify stiffness and permanent deformation.

Outcome and Learning Points

The preliminary stress check alone would approve the bracket, whereas the deflection check rejects it. The case illustrates that stiffness, strength, fatigue, joint behavior, and manufacturability are separate design requirements. It also demonstrates why analytical calculations are valuable for checking FEA and for identifying the dominant design variable before computational refinement.

Limitations, Safety, and Professional Practice

Strength calculations can be invalidated by fatigue, corrosion, wear, residual stress, temperature-dependent properties, creep, impact loading, manufacturing defects, and inaccurate loads. Safety factors should not be used to conceal poor modeling. Applicable machinery, pressure, lifting, structural, welding, and occupational-safety standards must be identified for the product and jurisdiction.

Engineers have a duty to communicate assumptions, uncertainty, inspection requirements, and consequences of failure. Critical components require traceable material data, controlled fabrication, appropriate nondestructive testing, and independent review. AI-assisted optimization and automated FEA can accelerate design-space exploration, but generated geometry and predictions require human verification, physical constraints, version control, and documented validation. A model that is statistically plausible but mechanically inconsistent must not be accepted.

Best-Practice Design Procedure

  1. Draw the free-body diagram and identify all load cases.
  2. Establish units, coordinates, sign conventions, and material properties.
  3. Calculate reactions and internal forces.
  4. Check axial, shear, bending, torsional, thermal, and stability effects.
  5. Evaluate both strength and stiffness.
  6. Account for concentrations, fatigue, joints, imperfections, and environment.
  7. Compare stresses with the appropriate failure criterion and design code.
  8. Validate critical predictions with test or verified FEA.
  9. Document assumptions, margins, mesh convergence, and limitations.
  10. Review the design for safety, manufacturability, maintainability, and sustainable material use.

Chapter Summary

Stress is internal force per unit area, while strain measures deformation relative to original dimensions. Linear elasticity relates normal stress to strain through Young’s modulus and shear stress to shear strain through the shear modulus. Axial members, shafts, beams, and columns require different models: axial stress uses P/A, torsion uses Tc/J, bending uses Mc/I, beam shear uses VQ/(It), and elastic buckling uses π²EI/(KL)².

Thermal expansion becomes thermal stress when restrained. Combined loading requires superposition only within linear-elastic assumptions and may be assessed using principal stresses or von Mises stress. Geometric discontinuities amplify local stress and must be managed in static and fatigue design. Beam deflection and column instability may govern before material strength is reached.

Analytical models provide transparent first estimates and essential checks for testing and FEA. Reliable engineering practice combines equations, realistic boundary conditions, material data, verification, physical testing, applicable standards, and professional judgment.

Key Terms

Stress: Internal force intensity, commonly force divided by area.

Strain: Deformation divided by original dimension.

Young’s modulus: Elastic modulus relating uniaxial normal stress to strain.

Shear modulus: Elastic modulus relating shear stress to shear strain.

Poisson’s ratio: Ratio of transverse strain to axial strain in uniaxial loading, with sign convention considered.

Elasticity: Reversible deformation after unloading.

Plasticity: Permanent deformation after unloading.

Normal stress: Stress acting perpendicular to a surface.

Shear stress: Stress acting tangentially to a surface.

Bending moment: Internal moment that produces curvature in a beam.

Torsion: Twisting produced by an applied torque.

Second moment of area: Geometric measure governing bending resistance.

Polar second moment: Geometric measure governing circular-shaft torsion.

Stress concentration: Local stress amplification caused by geometric discontinuity.

Principal stress: Normal stress on a plane where shear stress is zero.

Von Mises stress: Distortion-energy equivalent stress used commonly for ductile-metal yielding.

Deflection: Translational displacement of a structural member.

Buckling: Instability involving rapid lateral deformation under compression.

Finite-element analysis: Numerical method that approximates field behavior by subdividing a body into elements.

Factor of safety: Ratio representing design capacity relative to the specified demand or allowable criterion.

Review Questions

Short Answer Questions

  1. Define normal stress and normal strain, and state their SI units.
  2. Why can a component satisfy a strength requirement but fail a stiffness requirement?
  3. State the assumptions behind the circular-shaft torsion equation.
  4. What physical features commonly produce stress concentrations?
  5. Why is the least second moment of area used in column buckling calculations?

Long Answer Questions

  1. Explain the engineering stress–strain curve for a ductile metal and distinguish stiffness, strength, and toughness.
  2. Derive the relationship between torque, shear stress, and angle of twist for a uniform circular shaft.
  3. Explain how bending stress and transverse shear stress vary across a rectangular beam section.
  4. Discuss the assumptions and limitations of Euler column theory.
  5. Describe a defensible workflow for validating an FEA model using analytical calculations and physical testing.

Analytical / Application Questions

  1. A 16 mm diameter steel rod carries 25 kN in tension. Using E = 200 GPa, calculate its average stress and elastic strain. State whether the result can be accepted without a yield-strength value.
  2. A simply supported beam of span 1.2 m carries a central 800 N load. Explain which section properties control maximum bending stress and maximum deflection, and identify design changes that would reduce both.
  3. A motor shaft experiences simultaneous bending and torsion. Develop a calculation sequence using axial/bending normal stress, torsional shear stress, principal stress, and von Mises stress. State which additional information is required for a design decision.

Multiple Choice Questions with Answers

Question 1

For a prismatic member carrying a concentric axial load P over area A, the average normal stress is:

A. P/A

B. PA

C. A/P

D. P²/A

Correct Answer: A. P/A

Question 2

A fully restrained elastic bar is heated by ΔT. Which expression gives the induced thermal stress for a uniform uniaxial member?

A. EΔT/α

B. −EαΔT

C. αΔT/E

D. −α/(EΔT)

Correct Answer: B. −EαΔT

Question 3

For a solid circular shaft, the polar second moment of area is:

A. πd²/4

B. πd⁴/64

C. πd⁴/32

D. bd³/12

Correct Answer: C. πd⁴/32

Question 4

In elastic beam bending, the maximum normal stress occurs at the location where:

A. The shear force is zero in every case

B. The distance from the neutral axis is greatest

C. The second moment of area is zero

D. The beam width is greatest in every case

Correct Answer: B. The distance from the neutral axis is greatest

Question 5

The maximum deflection of a cantilever with an end point load P is:

A. PL/(3EI)

B. PL²/(2EI)

C. PL³/(3EI)

D. PL⁴/(8EI)

Correct Answer: C. PL³/(3EI)

Question 6

Which quantity is included in the beam shear formula τ = VQ/(It)?

A. The first moment of area Q

B. The polar radius c only

C. The thermal expansion coefficient α

D. The column effective-length factor K only

Correct Answer: A. The first moment of area Q

Question 7

For plane stress, the von Mises equivalent stress is used primarily to assess:

A. Ductile-metal yielding

B. Thermal expansion alone

C. Beam support reactions only

D. Fluid pressure without material data

Correct Answer: A. Ductile-metal yielding

Question 8

The Euler critical load of an ideal column varies with unsupported length L as:

A. L²

B. L

C. 1/L

D. 1/L²

Correct Answer: D. 1/L²

Question 9

The theoretical stress concentration factor K_t is defined as:

A. Nominal stress divided by local maximum stress

B. Local maximum stress divided by nominal stress

C. Yield stress divided by elastic modulus

D. Shear strain divided by normal strain

Correct Answer: B. Local maximum stress divided by nominal stress

Question 10

A linear-elastic FEA result at a mathematically sharp corner shows an extremely high peak stress. The most appropriate first interpretation is that:

A. The material must have infinite strength

B. The peak may be a stress singularity requiring realistic geometry or stress interpretation

C. The mesh is always sufficiently converged

D. The load should automatically be deleted

Correct Answer: B. The peak may be a stress singularity requiring realistic geometry or stress interpretation

Mechanical Engineering Fundamentals for Mechatronics
Chapter 4 of 15