Properties of Materials — Diagnostic Tests
Unit Tests
UT-1: Stress-Strain Curve Analysis
Question:
A metal wire of length 2.00 m 2.00\,\text{m} 2.00 m and cross-sectional area 1.5 × 10 − 6 m 2 1.5 \times 10^{-6}\,\text{m}^2 1.5 × 1 0 − 6 m 2 is tested under tension. The following data are obtained:
Stress / MPa 0 100 200 300 350 400 420 430 Strain / 10 − 3 10^{-3} 1 0 − 3 0 0.50 1.00 1.50 1.80 2.30 3.00 4.50
The wire breaks at a stress of 430 MPa 430\,\text{MPa} 430 MPa .
(a) Calculate Young's modulus from the linear region of the graph.
(b) Calculate the elastic potential energy stored in the wire at a stress of 300 MPa 300\,\text{MPa} 300 MPa using the stress-strain graph.
(c) Calculate the total energy per unit volume absorbed by the wire up to the point of fracture.
Solution:
(a) Young's modulus is the gradient of the linear (proportional) region:
E = ◆ L B ◆ Δ σ ◆ R B ◆◆ L B ◆ Δ ε ◆ R B ◆ = ◆ L B ◆ 300 × 10 6 − 0 ◆ R B ◆◆ L B ◆ 1.50 × 10 − 3 − 0 ◆ R B ◆ = ◆ L B ◆ 300 × 10 6 ◆ R B ◆◆ L B ◆ 1.50 × 10 − 3 ◆ R B ◆ = 2.00 × 10 11 Pa = 200 GPa E = \frac◆LB◆\Delta\sigma◆RB◆◆LB◆\Delta\varepsilon◆RB◆ = \frac◆LB◆300 \times 10^6 - 0◆RB◆◆LB◆1.50 \times 10^{-3} - 0◆RB◆ = \frac◆LB◆300 \times 10^6◆RB◆◆LB◆1.50 \times 10^{-3}◆RB◆ = 2.00 \times 10^{11}\,\text{Pa} = 200\,\text{GPa} E = L ◆ B ◆Δ σ ◆ R B ◆◆ L B ◆Δ ε ◆ R B ◆ = L ◆ B ◆300 × 1 0 6 − 0◆ R B ◆◆ L B ◆1.50 × 1 0 − 3 − 0◆ R B ◆ = L ◆ B ◆300 × 1 0 6 ◆ R B ◆◆ L B ◆1.50 × 1 0 − 3 ◆ R B ◆ = 2.00 × 1 0 11 Pa = 200 GPa
This is consistent with steel.
(b) In the linear region (up to 300 MPa 300\,\text{MPa} 300 MPa ), the energy per unit volume is the area under the stress-strain curve:
u = 1 2 σ ε = 1 2 × 300 × 10 6 × 1.50 × 10 − 3 = 2.25 × 10 5 J m − 3 u = \frac{1}{2}\sigma\varepsilon = \frac{1}{2} \times 300 \times 10^6 \times 1.50 \times 10^{-3} = 2.25 \times 10^5\,\text{J}\,\text{m}^{-3} u = 2 1 σ ε = 2 1 × 300 × 1 0 6 × 1.50 × 1 0 − 3 = 2.25 × 1 0 5 J m − 3
Total elastic energy stored in the wire:
U = u × volume = 2.25 × 10 5 × 2.00 × 1.5 × 10 − 6 = 0.675 J U = u \times \text{volume} = 2.25 \times 10^5 \times 2.00 \times 1.5 \times 10^{-6} = 0.675\,\text{J} U = u × volume = 2.25 × 1 0 5 × 2.00 × 1.5 × 1 0 − 6 = 0.675 J
(c) Total energy per unit volume up to fracture = total area under the stress-strain curve.
Using the trapezium rule on the data:
u total = 1 2 [ ( 0 + 100 ) × 0.50 + ( 100 + 200 ) × 0.50 + ( 200 + 300 ) × 0.50 + ( 300 + 350 ) × 0.30 + ( 350 + 400 ) × 0.50 + ( 400 + 420 ) × 0.70 + ( 420 + 430 ) × 1.50 ] × 10 6 × 10 − 3 u_{\text{total}} = \frac{1}{2}\left[(0 + 100) \times 0.50 + (100 + 200) \times 0.50 + (200 + 300) \times 0.50 + (300 + 350) \times 0.30 + (350 + 400) \times 0.50 + (400 + 420) \times 0.70 + (420 + 430) \times 1.50\right] \times 10^6 \times 10^{-3} u total = 2 1 [ ( 0 + 100 ) × 0.50 + ( 100 + 200 ) × 0.50 + ( 200 + 300 ) × 0.50 + ( 300 + 350 ) × 0.30 + ( 350 + 400 ) × 0.50 + ( 400 + 420 ) × 0.70 + ( 420 + 430 ) × 1.50 ] × 1 0 6 × 1 0 − 3
= 1 2 [ 50 + 150 + 250 + 195 + 375 + 574 + 1275 ] × 10 3 = \frac{1}{2}\left[50 + 150 + 250 + 195 + 375 + 574 + 1275\right] \times 10^3 = 2 1 [ 50 + 150 + 250 + 195 + 375 + 574 + 1275 ] × 1 0 3
= 1 2 × 2869 × 10 3 = 1.435 × 10 6 J m − 3 = \frac{1}{2} \times 2869 \times 10^3 = 1.435 \times 10^6\,\text{J}\,\text{m}^{-3} = 2 1 × 2869 × 1 0 3 = 1.435 × 1 0 6 J m − 3
The elastic energy recoverable (if unloaded from 300 MPa 300\,\text{MPa} 300 MPa ) is only 2.25 × 10 5 J m − 3 2.25 \times 10^5\,\text{J}\,\text{m}^{-3} 2.25 × 1 0 5 J m − 3 , so the majority of the energy is dissipated as heat during plastic deformation.
UT-2: Composite Wire Under Load
Question:
A composite wire consists of a steel section of length 1.5 m 1.5\,\text{m} 1.5 m and cross-sectional area 2.0 × 10 − 6 m 2 2.0 \times 10^{-6}\,\text{m}^2 2.0 × 1 0 − 6 m 2 joined to a copper section of length 1.0 m 1.0\,\text{m} 1.0 m and cross-sectional area 2.0 × 10 − 6 m 2 2.0 \times 10^{-6}\,\text{m}^2 2.0 × 1 0 − 6 m 2 . A tensile force of 400 N 400\,\text{N} 400 N is applied to the composite wire.
Young's modulus of steel = 2.0 × 10 11 Pa = 2.0 \times 10^{11}\,\text{Pa} = 2.0 × 1 0 11 Pa ; Young's modulus of copper = 1.2 × 10 11 Pa = 1.2 \times 10^{11}\,\text{Pa} = 1.2 × 1 0 11 Pa .
(a) Calculate the stress in each section.
(b) Calculate the strain in each section.
(c) Calculate the total extension of the composite wire.
Solution:
(a) Since the wire is in series, the force is the same in both sections:
σ = F A = ◆ L B ◆ 400 ◆ R B ◆◆ L B ◆ 2.0 × 10 − 6 ◆ R B ◆ = 2.0 × 10 8 Pa = 200 MPa \sigma = \frac{F}{A} = \frac◆LB◆400◆RB◆◆LB◆2.0 \times 10^{-6}◆RB◆ = 2.0 \times 10^8\,\text{Pa} = 200\,\text{MPa} σ = A F = L ◆ B ◆400◆ R B ◆◆ L B ◆2.0 × 1 0 − 6 ◆ R B ◆ = 2.0 × 1 0 8 Pa = 200 MPa
Both sections experience the same stress of 200 MPa 200\,\text{MPa} 200 MPa .
(b) Strain in steel: ε s = ◆ L B ◆ σ ◆ R B ◆◆ L B ◆ E s ◆ R B ◆ = ◆ L B ◆ 2.0 × 10 8 ◆ R B ◆◆ L B ◆ 2.0 × 10 11 ◆ R B ◆ = 1.0 × 10 − 3 \varepsilon_s = \frac◆LB◆\sigma◆RB◆◆LB◆E_s◆RB◆ = \frac◆LB◆2.0 \times 10^8◆RB◆◆LB◆2.0 \times 10^{11}◆RB◆ = 1.0 \times 10^{-3} ε s = L ◆ B ◆ σ ◆ R B ◆◆ L B ◆ E s ◆ R B ◆ = L ◆ B ◆2.0 × 1 0 8 ◆ R B ◆◆ L B ◆2.0 × 1 0 11 ◆ R B ◆ = 1.0 × 1 0 − 3
Strain in copper: ε c = ◆ L B ◆ σ ◆ R B ◆◆ L B ◆ E c ◆ R B ◆ = ◆ L B ◆ 2.0 × 10 8 ◆ R B ◆◆ L B ◆ 1.2 × 10 11 ◆ R B ◆ = 1.667 × 10 − 3 \varepsilon_c = \frac◆LB◆\sigma◆RB◆◆LB◆E_c◆RB◆ = \frac◆LB◆2.0 \times 10^8◆RB◆◆LB◆1.2 \times 10^{11}◆RB◆ = 1.667 \times 10^{-3} ε c = L ◆ B ◆ σ ◆ R B ◆◆ L B ◆ E c ◆ R B ◆ = L ◆ B ◆2.0 × 1 0 8 ◆ R B ◆◆ L B ◆1.2 × 1 0 11 ◆ R B ◆ = 1.667 × 1 0 − 3
(c) Extension of steel: Δ l s = ε s × l s = 1.0 × 10 − 3 × 1.5 = 1.50 × 10 − 3 m = 1.50 mm \Delta l_s = \varepsilon_s \times l_s = 1.0 \times 10^{-3} \times 1.5 = 1.50 \times 10^{-3}\,\text{m} = 1.50\,\text{mm} Δ l s = ε s × l s = 1.0 × 1 0 − 3 × 1.5 = 1.50 × 1 0 − 3 m = 1.50 mm
Extension of copper: Δ l c = ε c × l c = 1.667 × 10 − 3 × 1.0 = 1.667 × 10 − 3 m = 1.67 mm \Delta l_c = \varepsilon_c \times l_c = 1.667 \times 10^{-3} \times 1.0 = 1.667 \times 10^{-3}\,\text{m} = 1.67\,\text{mm} Δ l c = ε c × l c = 1.667 × 1 0 − 3 × 1.0 = 1.667 × 1 0 − 3 m = 1.67 mm
Total extension = 1.50 + 1.67 = 3.17 mm = 1.50 + 1.67 = 3.17\,\text{mm} = 1.50 + 1.67 = 3.17 mm
The effective Young's modulus of the composite wire:
E eff = ◆ L B ◆ F ( l s + l c ) ◆ R B ◆◆ L B ◆ A ( Δ l s + Δ l c ) ◆ R B ◆ = ◆ L B ◆ 400 × 2.5 ◆ R B ◆◆ L B ◆ 2.0 × 10 − 6 × 3.17 × 10 − 3 ◆ R B ◆ = ◆ L B ◆ 1000 ◆ R B ◆◆ L B ◆ 6.34 × 10 − 9 ◆ R B ◆ = 1.58 × 10 11 Pa E_{\text{eff}} = \frac◆LB◆F(l_s + l_c)◆RB◆◆LB◆A(\Delta l_s + \Delta l_c)◆RB◆ = \frac◆LB◆400 \times 2.5◆RB◆◆LB◆2.0 \times 10^{-6} \times 3.17 \times 10^{-3}◆RB◆ = \frac◆LB◆1000◆RB◆◆LB◆6.34 \times 10^{-9}◆RB◆ = 1.58 \times 10^{11}\,\text{Pa} E eff = L ◆ B ◆ F ( l s + l c ) ◆ R B ◆◆ L B ◆ A ( Δ l s + Δ l c ) ◆ R B ◆ = L ◆ B ◆400 × 2.5◆ R B ◆◆ L B ◆2.0 × 1 0 − 6 × 3.17 × 1 0 − 3 ◆ R B ◆ = L ◆ B ◆1000◆ R B ◆◆ L B ◆6.34 × 1 0 − 9 ◆ R B ◆ = 1.58 × 1 0 11 Pa
UT-3: Elastic Limit and Proportional Limit
Question:
A rubber band has the following force-extension data:
Force / N 0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 6.0 5.0 4.0 3.0 2.0 1.0 0 Extension / mm 0 15 35 60 90 130 180 250 200 170 150 130 110 80 20
The loading data are the first 8 readings and the unloading data are the last 8 readings.
(a) Determine whether the rubber band obeys Hooke's law at any point in its loading curve.
(b) Calculate the work done on the rubber band during loading and the work done by the rubber band during unloading.
(c) Calculate the energy dissipated as heat per loading-unloading cycle and explain its significance.
Solution:
(a) For Hooke's law to apply, force must be proportional to extension (constant gradient).
Checking ratios: F / x = 1.0 / 15 = 0.067 F/x = 1.0/15 = 0.067 F / x = 1.0/15 = 0.067 , 2.0 / 35 = 0.057 2.0/35 = 0.057 2.0/35 = 0.057 , 3.0 / 60 = 0.050 3.0/60 = 0.050 3.0/60 = 0.050 , 4.0 / 90 = 0.044 4.0/90 = 0.044 4.0/90 = 0.044
The ratio decreases, so the rubber band does not obey Hooke's law at any point during loading. The spring constant (gradient) continuously decreases, meaning the rubber becomes easier to stretch as it extends. This is characteristic of non-Hookean materials like rubber and polymers.
(b) Work done during loading (area under loading curve, using trapezium rule):
W load = 1 2 [ ( 0 + 1 ) × 15 + ( 1 + 2 ) × 20 + ( 2 + 3 ) × 25 + ( 3 + 4 ) × 30 + ( 4 + 5 ) × 40 + ( 5 + 6 ) × 50 + ( 6 + 7 ) × 70 ] × 10 − 3 W_{\text{load}} = \frac{1}{2}\left[(0+1) \times 15 + (1+2) \times 20 + (2+3) \times 25 + (3+4) \times 30 + (4+5) \times 40 + (5+6) \times 50 + (6+7) \times 70\right] \times 10^{-3} W load = 2 1 [ ( 0 + 1 ) × 15 + ( 1 + 2 ) × 20 + ( 2 + 3 ) × 25 + ( 3 + 4 ) × 30 + ( 4 + 5 ) × 40 + ( 5 + 6 ) × 50 + ( 6 + 7 ) × 70 ] × 1 0 − 3
= 1 2 [ 15 + 60 + 125 + 210 + 360 + 550 + 910 ] × 10 − 3 = \frac{1}{2}\left[15 + 60 + 125 + 210 + 360 + 550 + 910\right] \times 10^{-3} = 2 1 [ 15 + 60 + 125 + 210 + 360 + 550 + 910 ] × 1 0 − 3
= 1 2 × 2230 × 10 − 3 = 1.115 J = \frac{1}{2} \times 2230 \times 10^{-3} = 1.115\,\text{J} = 2 1 × 2230 × 1 0 − 3 = 1.115 J
Work done during unloading (area under unloading curve, using the trapezium rule):
W unload = ∑ F i + F i + 1 2 × Δ x W_{\text{unload}} = \sum \frac{F_i + F_{i+1}}{2} \times \Delta x W unload = ∑ 2 F i + F i + 1 × Δ x
= 6 + 5 2 ( 250 − 200 ) + 5 + 4 2 ( 200 − 170 ) + 4 + 3 2 ( 170 − 150 ) + 3 + 2 2 ( 150 − 130 ) + 2 + 1 2 ( 130 − 110 ) + 1 + 0 2 ( 110 − 80 ) + 0 + 0 2 ( 80 − 20 ) = \frac{6+5}{2}(250-200) + \frac{5+4}{2}(200-170) + \frac{4+3}{2}(170-150) + \frac{3+2}{2}(150-130) + \frac{2+1}{2}(130-110) + \frac{1+0}{2}(110-80) + \frac{0+0}{2}(80-20) = 2 6 + 5 ( 250 − 200 ) + 2 5 + 4 ( 200 − 170 ) + 2 4 + 3 ( 170 − 150 ) + 2 3 + 2 ( 150 − 130 ) + 2 2 + 1 ( 130 − 110 ) + 2 1 + 0 ( 110 − 80 ) + 2 0 + 0 ( 80 − 20 )
= 275 + 135 + 70 + 50 + 30 + 15 + 0 = 575 mJ = 0.575 J = 275 + 135 + 70 + 50 + 30 + 15 + 0 = 575\,\text{mJ} = 0.575\,\text{J} = 275 + 135 + 70 + 50 + 30 + 15 + 0 = 575 mJ = 0.575 J
(c) Energy dissipated per cycle = W load − W unload = 1.115 − 0.575 = 0.540 J = W_{\text{load}} - W_{\text{unload}} = 1.115 - 0.575 = 0.540\,\text{J} = W load − W unload = 1.115 − 0.575 = 0.540 J
This energy is dissipated as heat due to internal friction within the rubber (hysteresis). The loading-unloading curve forms a hysteresis loop, and the area enclosed represents the energy lost per cycle. This is why rubber gets warm when repeatedly stretched and released. The fraction of energy dissipated is 0.540 / 1.115 = 48 % 0.540/1.115 = 48\% 0.540/1.115 = 48% , meaning the rubber band returns only about 52 % 52\% 52% of the energy stored in it.
Integration Tests
IT-1: Wire Under Thermal and Mechanical Stress (with Dynamics)
Question:
A steel wire of length 3.0 m 3.0\,\text{m} 3.0 m , cross-sectional area 2.0 × 10 − 6 m 2 2.0 \times 10^{-6}\,\text{m}^2 2.0 × 1 0 − 6 m 2 , and Young's modulus 2.0 × 10 11 Pa 2.0 \times 10^{11}\,\text{Pa} 2.0 × 1 0 11 Pa is fixed between two rigid supports. The linear expansivity of steel is 1.2 × 10 − 5 K − 1 1.2 \times 10^{-5}\,\text{K}^{-1} 1.2 × 1 0 − 5 K − 1 .
(a) The temperature increases by 40 K 40\,\text{K} 40 K . Calculate the stress that develops in the wire if it is not allowed to expand.
(b) The wire is now allowed to expand freely by Δ T = 40 K \Delta T = 40\,\text{K} Δ T = 40 K , and then a mass of 50 kg 50\,\text{kg} 50 kg is hung from it. Calculate the total extension.
(c) Calculate the elastic potential energy stored in the wire in part (b).
Take g = 9.81 m s − 2 g = 9.81\,\text{m}\,\text{s}^{-2} g = 9.81 m s − 2 .
Solution:
(a) If the wire cannot expand, the thermal expansion is fully converted to elastic strain:
Δ l thermal = α l Δ T = 1.2 × 10 − 5 × 3.0 × 40 = 1.44 × 10 − 3 m \Delta l_{\text{thermal}} = \alpha l \Delta T = 1.2 \times 10^{-5} \times 3.0 \times 40 = 1.44 \times 10^{-3}\,\text{m} Δ l thermal = α l Δ T = 1.2 × 1 0 − 5 × 3.0 × 40 = 1.44 × 1 0 − 3 m
The wire must be compressed back by this amount, so the strain is:
ε = ◆ L B ◆ Δ l ◆ R B ◆◆ L B ◆ l ◆ R B ◆ = ◆ L B ◆ 1.44 × 10 − 3 ◆ R B ◆◆ L B ◆ 3.0 ◆ R B ◆ = 4.8 × 10 − 4 \varepsilon = \frac◆LB◆\Delta l◆RB◆◆LB◆l◆RB◆ = \frac◆LB◆1.44 \times 10^{-3}◆RB◆◆LB◆3.0◆RB◆ = 4.8 \times 10^{-4} ε = L ◆ B ◆Δ l ◆ R B ◆◆ L B ◆ l ◆ R B ◆ = L ◆ B ◆1.44 × 1 0 − 3 ◆ R B ◆◆ L B ◆3.0◆ R B ◆ = 4.8 × 1 0 − 4
Stress: σ = E ε = 2.0 × 10 11 × 4.8 × 10 − 4 = 9.6 × 10 7 Pa = 96 MPa \sigma = E\varepsilon = 2.0 \times 10^{11} \times 4.8 \times 10^{-4} = 9.6 \times 10^7\,\text{Pa} = 96\,\text{MPa} σ = E ε = 2.0 × 1 0 11 × 4.8 × 1 0 − 4 = 9.6 × 1 0 7 Pa = 96 MPa
This is a compressive stress. The wire would buckle if not properly supported.
(b) Free thermal expansion produces no stress. The extension due to the hanging mass:
Δ l mechanical = F l A E = ◆ L B ◆ 50 × 9.81 × 3.0 ◆ R B ◆◆ L B ◆ 2.0 × 10 − 6 × 2.0 × 10 11 ◆ R B ◆ = ◆ L B ◆ 1471.5 ◆ R B ◆◆ L B ◆ 4.0 × 10 5 ◆ R B ◆ = 3.68 × 10 − 3 m \Delta l_{\text{mechanical}} = \frac{Fl}{AE} = \frac◆LB◆50 \times 9.81 \times 3.0◆RB◆◆LB◆2.0 \times 10^{-6} \times 2.0 \times 10^{11}◆RB◆ = \frac◆LB◆1471.5◆RB◆◆LB◆4.0 \times 10^5◆RB◆ = 3.68 \times 10^{-3}\,\text{m} Δ l mechanical = A E F l = L ◆ B ◆50 × 9.81 × 3.0◆ R B ◆◆ L B ◆2.0 × 1 0 − 6 × 2.0 × 1 0 11 ◆ R B ◆ = L ◆ B ◆1471.5◆ R B ◆◆ L B ◆4.0 × 1 0 5 ◆ R B ◆ = 3.68 × 1 0 − 3 m
Total extension = Δ l thermal + Δ l mechanical = 1.44 + 3.68 = 5.12 × 10 − 3 m = 5.12 mm = \Delta l_{\text{thermal}} + \Delta l_{\text{mechanical}} = 1.44 + 3.68 = 5.12 \times 10^{-3}\,\text{m} = 5.12\,\text{mm} = Δ l thermal + Δ l mechanical = 1.44 + 3.68 = 5.12 × 1 0 − 3 m = 5.12 mm
(c) Only the mechanical extension stores elastic energy (thermal expansion without constraint does not store elastic energy):
U = 1 2 × stress × strain × volume U = \frac{1}{2} \times \text{stress} \times \text{strain} \times \text{volume} U = 2 1 × stress × strain × volume
= 1 2 × F A × ◆ L B ◆ Δ l mech ◆ R B ◆◆ L B ◆ l ◆ R B ◆ × A l = 1 2 F Δ l mech = \frac{1}{2} \times \frac{F}{A} \times \frac◆LB◆\Delta l_{\text{mech}}◆RB◆◆LB◆l◆RB◆ \times Al = \frac{1}{2}F\Delta l_{\text{mech}} = 2 1 × A F × L ◆ B ◆Δ l mech ◆ R B ◆◆ L B ◆ l ◆ R B ◆ × A l = 2 1 F Δ l mech
= 1 2 × 50 × 9.81 × 3.68 × 10 − 3 = 0.903 J = \frac{1}{2} \times 50 \times 9.81 \times 3.68 \times 10^{-3} = 0.903\,\text{J} = 2 1 × 50 × 9.81 × 3.68 × 1 0 − 3 = 0.903 J
IT-2: Stress in a Rotating Ring (with Circular Motion)
Question:
A thin ring of radius R = 0.50 m R = 0.50\,\text{m} R = 0.50 m is made of steel with density 7800 kg m − 3 7800\,\text{kg}\,\text{m}^{-3} 7800 kg m − 3 and Young's modulus 2.0 × 10 11 Pa 2.0 \times 10^{11}\,\text{Pa} 2.0 × 1 0 11 Pa . The ring has a square cross-section of side 5.0 mm 5.0\,\text{mm} 5.0 mm and rotates about its centre with angular velocity ω \omega ω .
(a) Derive an expression for the tensile stress in the ring due to rotation.
(b) Calculate the angular velocity at which the stress reaches 200 MPa 200\,\text{MPa} 200 MPa .
(c) Calculate the percentage increase in the ring's radius at this angular velocity.
Solution:
(a) Consider a small element of the ring of mass δ m \delta m δ m subtending angle δ θ \delta\theta δ θ at the centre. The element is in circular motion with radius R R R .
Forces on the element: tension T T T on each side, at angle δ θ / 2 \delta\theta/2 δ θ /2 to the tangent.
Net inward force = 2 T sin ( δ θ / 2 ) ≈ T δ θ = 2T\sin(\delta\theta/2) \approx T\delta\theta = 2 T sin ( δ θ /2 ) ≈ T δ θ for small δ θ \delta\theta δ θ .
Centripetal force: δ m × ω 2 R = ( ρ A R δ θ ) ω 2 R \delta m \times \omega^2 R = (\rho A R\delta\theta)\omega^2 R δ m × ω 2 R = ( ρ A R δ θ ) ω 2 R
T δ θ = ρ A R 2 ω 2 δ θ T\delta\theta = \rho A R^2\omega^2\delta\theta T δ θ = ρ A R 2 ω 2 δ θ
T = ρ A R 2 ω 2 T = \rho A R^2\omega^2 T = ρ A R 2 ω 2
Stress: σ = T / A = ρ R 2 ω 2 \sigma = T/A = \rho R^2\omega^2 σ = T / A = ρ R 2 ω 2
(b) σ = ρ R 2 ω 2 \sigma = \rho R^2\omega^2 σ = ρ R 2 ω 2
ω = ◆ L B ◆ ◆ L B ◆ σ ◆ R B ◆◆ L B ◆ ρ R 2 ◆ R B ◆◆ R B ◆ = ◆ L B ◆ ◆ L B ◆ 200 × 10 6 ◆ R B ◆◆ L B ◆ 7800 × 0.25 ◆ R B ◆◆ R B ◆ = ◆ L B ◆ ◆ L B ◆ 200 × 10 6 ◆ R B ◆◆ L B ◆ 1950 ◆ R B ◆◆ R B ◆ = ◆ L B ◆ 1.026 × 10 5 ◆ R B ◆ = 320 rad s − 1 \omega = \sqrt◆LB◆\frac◆LB◆\sigma◆RB◆◆LB◆\rho R^2◆RB◆◆RB◆ = \sqrt◆LB◆\frac◆LB◆200 \times 10^6◆RB◆◆LB◆7800 \times 0.25◆RB◆◆RB◆ = \sqrt◆LB◆\frac◆LB◆200 \times 10^6◆RB◆◆LB◆1950◆RB◆◆RB◆ = \sqrt◆LB◆1.026 \times 10^5◆RB◆ = 320\,\text{rad}\,\text{s}^{-1} ω = ◆ L B ◆ L ◆ B ◆ σ ◆ R B ◆◆ L B ◆ ρ R 2 ◆ R B ◆◆ R B ◆ = ◆ L B ◆ L ◆ B ◆200 × 1 0 6 ◆ R B ◆◆ L B ◆7800 × 0.25◆ R B ◆◆ R B ◆ = ◆ L B ◆ L ◆ B ◆200 × 1 0 6 ◆ R B ◆◆ L B ◆1950◆ R B ◆◆ R B ◆ = ◆ L B ◆1.026 × 1 0 5 ◆ R B ◆ = 320 rad s − 1
This is about 3060 rpm 3060\,\text{rpm} 3060 rpm , or 51 rev s − 1 51\,\text{rev}\,\text{s}^{-1} 51 rev s − 1 .
(c) Strain: ε = σ / E = 200 × 10 6 / ( 2.0 × 10 11 ) = 1.0 × 10 − 3 \varepsilon = \sigma/E = 200 \times 10^6/(2.0 \times 10^{11}) = 1.0 \times 10^{-3} ε = σ / E = 200 × 1 0 6 / ( 2.0 × 1 0 11 ) = 1.0 × 1 0 − 3
Percentage increase in radius = ε × 100 = 0.10 % = \varepsilon \times 100 = 0.10\% = ε × 100 = 0.10%
The circumference also increases by 0.10 % 0.10\% 0.10% .
IT-3: Energy Stored in a Stretched Wire Hanging Under Its Own Weight (with Work-Energy)
Question:
A steel wire of length 10 m 10\,\text{m} 10 m , cross-sectional area 1.0 × 10 − 6 m 2 1.0 \times 10^{-6}\,\text{m}^2 1.0 × 1 0 − 6 m 2 , Young's modulus 2.0 × 10 11 Pa 2.0 \times 10^{11}\,\text{Pa} 2.0 × 1 0 11 Pa , and density 7800 kg m − 3 7800\,\text{kg}\,\text{m}^{-3} 7800 kg m − 3 hangs vertically from a fixed support.
(a) Calculate the total extension of the wire under its own weight.
(b) Calculate the elastic potential energy stored in the wire. Compare this with the naive calculation U = 1 2 × ( total weight ) × ( total extension ) U = \frac{1}{2} \times (\text{total weight}) \times (\text{total extension}) U = 2 1 × ( total weight ) × ( total extension ) .
(c) Explain why the two values differ.
Solution:
(a) The stress varies linearly along the wire. At a distance y y y from the bottom, the stress is due to the weight of wire below:
σ ( y ) = ◆ L B ◆ ρ A g y ◆ R B ◆◆ L B ◆ A ◆ R B ◆ = ρ g y \sigma(y) = \frac◆LB◆\rho A g y◆RB◆◆LB◆A◆RB◆ = \rho g y σ ( y ) = L ◆ B ◆ ρ A g y ◆ R B ◆◆ L B ◆ A ◆ R B ◆ = ρ g y
where y y y is measured from the bottom (0 ≤ y ≤ L 0 \le y \le L 0 ≤ y ≤ L ).
The strain at position y y y : ε ( y ) = σ ( y ) / E = ρ g y / E \varepsilon(y) = \sigma(y)/E = \rho g y/E ε ( y ) = σ ( y ) / E = ρ g y / E
Total extension:
Δ L = ∫ 0 L ε ( y ) d y = ∫ 0 L ◆ L B ◆ ρ g y ◆ R B ◆◆ L B ◆ E ◆ R B ◆ d y = ◆ L B ◆ ρ g L 2 ◆ R B ◆◆ L B ◆ 2 E ◆ R B ◆ \Delta L = \int_0^L \varepsilon(y)\,dy = \int_0^L \frac◆LB◆\rho g y◆RB◆◆LB◆E◆RB◆\,dy = \frac◆LB◆\rho g L^2◆RB◆◆LB◆2E◆RB◆ Δ L = ∫ 0 L ε ( y ) d y = ∫ 0 L L ◆ B ◆ ρ g y ◆ R B ◆◆ L B ◆ E ◆ R B ◆ d y = L ◆ B ◆ ρ g L 2 ◆ R B ◆◆ L B ◆2 E ◆ R B ◆
= ◆ L B ◆ 7800 × 9.81 × 100 ◆ R B ◆◆ L B ◆ 2 × 2.0 × 10 11 ◆ R B ◆ = ◆ L B ◆ 7.651 × 10 6 ◆ R B ◆◆ L B ◆ 4.0 × 10 11 ◆ R B ◆ = 1.913 × 10 − 5 m = 0.0191 mm = \frac◆LB◆7800 \times 9.81 \times 100◆RB◆◆LB◆2 \times 2.0 \times 10^{11}◆RB◆ = \frac◆LB◆7.651 \times 10^6◆RB◆◆LB◆4.0 \times 10^{11}◆RB◆ = 1.913 \times 10^{-5}\,\text{m} = 0.0191\,\text{mm} = L ◆ B ◆7800 × 9.81 × 100◆ R B ◆◆ L B ◆2 × 2.0 × 1 0 11 ◆ R B ◆ = L ◆ B ◆7.651 × 1 0 6 ◆ R B ◆◆ L B ◆4.0 × 1 0 11 ◆ R B ◆ = 1.913 × 1 0 − 5 m = 0.0191 mm
(b) Elastic potential energy per unit volume at position y y y :
u ( y ) = 1 2 σ ( y ) ε ( y ) = 1 2 ◆ L B ◆ ( ρ g y ) 2 ◆ R B ◆◆ L B ◆ E ◆ R B ◆ u(y) = \frac{1}{2}\sigma(y)\varepsilon(y) = \frac{1}{2}\frac◆LB◆(\rho g y)^2◆RB◆◆LB◆E◆RB◆ u ( y ) = 2 1 σ ( y ) ε ( y ) = 2 1 L ◆ B ◆ ( ρ g y ) 2 ◆ R B ◆◆ L B ◆ E ◆ R B ◆
Total energy:
U = ∫ 0 L u ( y ) × A d y = A 2 E ∫ 0 L ( ρ g y ) 2 d y = ◆ L B ◆ A ρ 2 g 2 L 3 ◆ R B ◆◆ L B ◆ 6 E ◆ R B ◆ U = \int_0^L u(y) \times A\,dy = \frac{A}{2E}\int_0^L (\rho g y)^2\,dy = \frac◆LB◆A\rho^2 g^2 L^3◆RB◆◆LB◆6E◆RB◆ U = ∫ 0 L u ( y ) × A d y = 2 E A ∫ 0 L ( ρ g y ) 2 d y = L ◆ B ◆ A ρ 2 g 2 L 3 ◆ R B ◆◆ L B ◆6 E ◆ R B ◆
= ◆ L B ◆ 1.0 × 10 − 6 × ( 7800 ) 2 × ( 9.81 ) 2 × 1000 ◆ R B ◆◆ L B ◆ 6 × 2.0 × 10 11 ◆ R B ◆ = \frac◆LB◆1.0 \times 10^{-6} \times (7800)^2 \times (9.81)^2 \times 1000◆RB◆◆LB◆6 \times 2.0 \times 10^{11}◆RB◆ = L ◆ B ◆1.0 × 1 0 − 6 × ( 7800 ) 2 × ( 9.81 ) 2 × 1000◆ R B ◆◆ L B ◆6 × 2.0 × 1 0 11 ◆ R B ◆
where 7800 2 = 6.084 × 10 7 7800^2 = 6.084 \times 10^7 780 0 2 = 6.084 × 1 0 7 , 9.81 2 = 96.24 9.81^2 = 96.24 9.8 1 2 = 96.24 , L 3 = 1000 L^3 = 1000 L 3 = 1000 :
= ◆ L B ◆ 1.0 × 10 − 6 × 6.084 × 10 7 × 96.24 × 10 3 ◆ R B ◆◆ L B ◆ 6 × 2.0 × 10 11 ◆ R B ◆ = ◆ L B ◆ 5.855 × 10 6 ◆ R B ◆◆ L B ◆ 1.2 × 10 12 ◆ R B ◆ = 4.879 × 10 − 6 J = \frac◆LB◆1.0 \times 10^{-6} \times 6.084 \times 10^7 \times 96.24 \times 10^3◆RB◆◆LB◆6 \times 2.0 \times 10^{11}◆RB◆ = \frac◆LB◆5.855 \times 10^6◆RB◆◆LB◆1.2 \times 10^{12}◆RB◆ = 4.879 \times 10^{-6}\,\text{J} = L ◆ B ◆1.0 × 1 0 − 6 × 6.084 × 1 0 7 × 96.24 × 1 0 3 ◆ R B ◆◆ L B ◆6 × 2.0 × 1 0 11 ◆ R B ◆ = L ◆ B ◆5.855 × 1 0 6 ◆ R B ◆◆ L B ◆1.2 × 1 0 12 ◆ R B ◆ = 4.879 × 1 0 − 6 J
Naive calculation: U naive = 1 2 × ( ρ A L g ) × Δ L U_{\text{naive}} = \frac{1}{2} \times (\rho A L g) \times \Delta L U naive = 2 1 × ( ρ A Lg ) × Δ L
= 1 2 × ( 7800 × 1.0 × 10 − 6 × 10 × 9.81 ) × 1.913 × 10 − 5 = \frac{1}{2} \times (7800 \times 1.0 \times 10^{-6} \times 10 \times 9.81) \times 1.913 \times 10^{-5} = 2 1 × ( 7800 × 1.0 × 1 0 − 6 × 10 × 9.81 ) × 1.913 × 1 0 − 5
= 1 2 × 0.7651 × 1.913 × 10 − 5 = 7.317 × 10 − 6 J = \frac{1}{2} \times 0.7651 \times 1.913 \times 10^{-5} = 7.317 \times 10^{-6}\,\text{J} = 2 1 × 0.7651 × 1.913 × 1 0 − 5 = 7.317 × 1 0 − 6 J
(c) The ratio is U naive / U = 7.317 × 10 − 6 / 4.879 × 10 − 6 = 1.50 U_{\text{naive}}/U = 7.317 \times 10^{-6}/4.879 \times 10^{-6} = 1.50 U naive / U = 7.317 × 1 0 − 6 /4.879 × 1 0 − 6 = 1.50
The naive calculation overestimates by a factor of 1.5. This can be shown analytically:
U naive = 1 2 ( ρ A L g ) Δ L = 1 2 ρ A L g × ◆ L B ◆ ρ g L 2 ◆ R B ◆◆ L B ◆ 2 E ◆ R B ◆ = ◆ L B ◆ ρ 2 A g 2 L 3 ◆ R B ◆◆ L B ◆ 4 E ◆ R B ◆ U_{\text{naive}} = \frac{1}{2}(\rho A L g)\Delta L = \frac{1}{2}\rho A L g \times \frac◆LB◆\rho g L^2◆RB◆◆LB◆2E◆RB◆ = \frac◆LB◆\rho^2 A g^2 L^3◆RB◆◆LB◆4E◆RB◆ U naive = 2 1 ( ρ A Lg ) Δ L = 2 1 ρ A Lg × L ◆ B ◆ ρ g L 2 ◆ R B ◆◆ L B ◆2 E ◆ R B ◆ = L ◆ B ◆ ρ 2 A g 2 L 3 ◆ R B ◆◆ L B ◆4 E ◆ R B ◆
U = ◆ L B ◆ ρ 2 A g 2 L 3 ◆ R B ◆◆ L B ◆ 6 E ◆ R B ◆ U = \frac◆LB◆\rho^2 A g^2 L^3◆RB◆◆LB◆6E◆RB◆ U = L ◆ B ◆ ρ 2 A g 2 L 3 ◆ R B ◆◆ L B ◆6 E ◆ R B ◆
Ratio U naive / U = 6 / 4 = 3 / 2 = 1.5 U_{\text{naive}}/U = 6/4 = 3/2 = 1.5 U naive / U = 6/4 = 3/2 = 1.5
The naive calculation overestimates by a factor of 1.5 because the stress is not uniform along the wire -- it varies linearly from zero at the bottom to a maximum at the top. The naive formula assumes the entire wire experiences the maximum stress (full weight divided by area), whereas the correct integral accounts for the linear variation.