6.00: Physical Problems for Regression
- Page ID
- 126415
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Lesson 1: General Engineering
Summary
A student wants to pick the best mousetrap out of many for use in a mousetrap car. To do so, one needs to regress the torsion in the mousetrap spring to the angle through which it is opened.
Modeling
In the summer of 2003, a high school student worked with me to study the physics of a mousetrap car to build the car that would go the fastest and/or the farthest.
So what is a mousetrap car?
A mousetrap car (Figure \(\PageIndex{1.1}\)) is a vehicle that is powered by the energy stored in the spring of a mousetrap. One of the basic ways most mousetrap vehicles are set into motion is by connecting the lever of the mousetrap bar through a string to the car’s axle. As the mousetrap lever is released, the tension that was built up in the spring is released, and the car sets into motion.
So to get the car to go the farthest and fastest, we need to store as much potential energy in the mousetrap spring. Since there is no specification on the mousetrap other than using a certain brand, we could use a mousetrap with the largest spring constant as that would translate to the largest potential energy stored.
Therefore, we bought several mousetrap springs and conducted a simple experiment to determine the spring constant of the mousetrap spring. These values of the spring constant would then allow us to find the potential energy stored in the mousetraps.
To find the spring constant of the mousetrap spring, we conducted the following experiment. We pulled the end of the mousetrap lever by a force probe (a fishhook) and measured the angle through which it is pulled by a protractor (Figure \(\PageIndex{1.2}\)). We pulled on the end of the spring via a string so that the force on the lever is applied at a right angle to the lever. This way, the torque applied is simply the product of the force, as measured by the fishhook (Figure \(\PageIndex{1.3}\)) and the length of the lever arm.
We measured the force needed to pull the lever for several different angles. The data for one of the mousetraps is given in Table \(\PageIndex{1.1}\).
The relationship between the torque applied and the angle of the rotation of the spring rotation is assumed to be a straight line
\[T = k_{0} + k_{1}{\theta}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.1}) \nonumber\]
where
\[T = \text{Torque}\ \text{(N-m)} \nonumber\]
\[\theta = \text{Angle through which the spring is rotated} \nonumber\]
Angle (degrees) | Force (lbs) |
---|---|
\(40\) | \(0.9\) |
\(55\) | \(1.0\) |
\(65\) | \(1.1\) |
\(90\) | \(1.2\) |
\(110\) | \(1.5\) |
The data in Table \(\PageIndex{1.1}\) is converted to the SI system of units, and the torque is calculated using the measured lever moment arm of \(47\ \text{mm}\), that is,
\[T = FL\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.2}) \nonumber\]
where
\[T =\text{ Torque} \ \text{(N-m)} \nonumber\]
\[F =\text{Force applied}\ \text{(N)} \nonumber\]
\[L =\text{ Moment arm}\ \text{(m)} \nonumber\]
Other mousetrap springs were tested accordingly, and the data is given in Tables \(\PageIndex{1.2}\) and \(\PageIndex{1.3}\).
From these tables, we can find the \(T = k_{0} + k_{1}{\theta}\) relationship for each of three mousetraps. The potential energy, \(U\) stored in the spring would be given by
\[\begin{split} U &= \int_{\theta_{\text{low}}}^{\theta_{\text{high}}} T \ d \theta \\ &= \int_{\theta_{\text{low}}}^{\theta_{\text{high}}} \left(k_{0} + k_{1}\ \theta \right) \ d\theta \end{split}\]
Knowing that in our case, \(\theta_{\text{low}} = 0\) and \(\theta_{\text{high}} = \pi\), the maximum potential energy stored
\[\begin{split} U_{\max} &= \int_{0}^{\pi}(k_{0} + k_{1}\ \theta)d\theta\\ &= k_{0}\pi + k_{1}\frac{\pi^{2}}{2} \end{split}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.3})\]
Angle (degrees) | Force (lbs) |
---|---|
\(38\) | \(0.9\) |
\(58\) | \(1.0\) |
\(65\) | \(1.1\) |
\(90\) | \(1.2\) |
\(120\) | \(1.6\) |
Angle (degrees) | Force (lbs) |
---|---|
\(40\) | \(0.9\) |
\(57\) | \(1.0\) |
\(65\) | \(1.1\) |
\(90\) | \(1.2\) |
\(135\) | \(1.7\) |
Questions
(1) Find the constants of the regression model for the spring for the three cases (Tables 1, 2, and 3).
(2) Based on maximum potential energy stored, which of the three springs would you choose?
Lesson 2: Mechanical Engineering
Summary
To find the diametrical contraction during shrink-fitting in a trunnion requires one to first find how the thermal expansion coefficient of steel is related to temperature. The relation is given through a second-order polynomial regression model.
Modeling
To make the fulcrum (Figure \(\PageIndex{2.1}\)) of a bascule bridge, a long hollow steel shaft called the trunnion is shrunk-fit into a steel hub. The resulting steel trunnion-hub assembly is then shrink fit into the girder of the bridge.
This shrink-fitting is done by first immersing the trunnion in a cold medium such as a dry-ice/alcohol mixture. After the trunnion reaches the steady-state temperature, that is, the temperature of the cold medium, the trunnion outer diameter contracts. The trunnion is taken out of the medium and slid through the hole of the hub (Figure \(\PageIndex{1.2}\)).
When the trunnion heats up, it expands and creates an interference fit with the hub. In 1995, on one of the bridges in Florida, this assembly procedure did not work as designed. Before the trunnion could be inserted fully into the hub, the trunnion got stuck. Luckily, the trunnion was taken out before it got stuck permanently. Otherwise, a new trunnion and hub would need to have been ordered at the cost of $50,000. Coupled with construction delays, the total loss could have been more than a hundred thousand dollars.
Why did the trunnion get stuck? This failure was because the trunnion had not contracted enough to slide through the hole. Can you find out why?
A hollow trunnion of outside diameter \({12.36}{3}{^{\prime\prime} }\) is to be fitted in a hub of inner diameter \(12.358{^{\prime\prime} }\). The trunnion was put in a dry ice/alcohol mixture (temperature of the fluid - dry ice/alcohol mixture is \(- 108{^\circ}\text{F}\)) to contract the trunnion so that it can be slid through the hole of the hub. To slide the trunnion without sticking, a diametrical clearance of at least \(0.01{^{\prime\prime} }\) is required between the trunnion and the hub. Assuming the room temperature is \(80{^\circ}\text{F}\), is immersing it in a dry-ice/alcohol mixture a correct decision?
Solution
To calculate the contraction in the diameter of the trunnion, coefficient of linear thermal expansion at room temperature is used. In that case the reduction, \({\Delta D}\) in the outer diameter of the trunnion is
\[\Delta D = D\alpha\Delta T\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{2.1}) \nonumber\]
where
\[D = \text{outer diameter of the trunnion,} \nonumber\]
\[\alpha = \text{coefficient of linear thermal expansion at room temperature, and} \nonumber\]
\[\Delta T = \text{change in temperature.} \nonumber\]
Given
\[D = 12.363{^{\prime\prime} } \nonumber\]
\[\alpha = 6.817 \times 10^{- 6}\ \text{in/in/}{^\circ}\text{F}\text{ at }80{^\circ}\text{F} \nonumber\]
\[\begin{split} \Delta T&=T_{\text{fluid}} - T_{\text{room}}\\ &= - 188{^\circ}\text{F} \end{split}\]
where
\[T_{\text{fluid}} = \text{temperature of dry-ice/alcohol mixture,} \nonumber\]
\[T_{\text{room}} = \text{room temperature,} \nonumber\]
the reduction in the trunnion outer diameter (Equation \(\PageIndex{2.1}\)) is given by
\[\begin{split} \Delta D &= (12.363)\left( 6.47 \times 10^{- 6} \right)\left( - 188 \right)\\ &= - 0.01504^{\prime\prime} \end{split}\]
So the trunnion is predicted to reduce in diameter by \(0.01504{^{\prime\prime} }\). But is this enough reduction in diameter? As per specifications, we need the trunnion to contract by
\[\begin{split} &= \text{trunnion outside diameter - hub inner diameter + diametric clearance}\\ &= 12.363^{\prime\prime} - 12.358^{\prime\prime} + 0.01^{\prime\prime} \\ &= 0.015^{\prime\prime} \end{split}\]
So, according to these calculations, immersing the steel trunnion in a dry-ice/alcohol mixture gives the desired contraction of \(0.015^{\prime\prime}\) as we predict a contraction of \(0.01504{^{\prime\prime} }\).
But as shown in Figure \(\PageIndex{1.3}\), the coefficient of linear thermal expansion of steel decreases with temperature and is not constant over the range of temperature the trunnion goes through. Hence, the above formula (Equation \(\PageIndex{2.1}\)) would overestimate the thermal contraction.
The contraction in the diameter for the trunnion for which the coefficient of linear thermal expansion varies as a function of temperature is given by
\[\Delta D = D\int_{T_{\text{room}}}^{T_{\text{fluid}}} \alpha \ dT \;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{2.2}) \nonumber\]
So one needs to find the curve to find the coefficient of linear thermal expansion as a function of temperature. This curve is found by regression where we best fit a function through the data given in Table 1. In this case, we may fit a second-order polynomial
\[\alpha = a_{0} + a_{1} \times T + a_{2} \times T^{2}.\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{2.3}) \nonumber\]
The values of the coefficients in the above equations will be found by polynomial regression. Knowing the values of \(a_{0}\), \(a_{1}\),and \(a_{2}\), we can then find the contraction in the trunnion diameter from Equation (2) and (3) as
\[\begin{split} \Delta D &= D\int_{T_{\text{room}}}^{T_{\text{fluid}}}{(a_{0} + a_{1} \times T + a_{2} \times T^{2}}) \ dT\\ &= D\lbrack a_{0} \times (T_{\text{fluid}} - T_{\text{room}}) + a_{1} \times \frac{({T_{\text{fluid}}}^{2} - {T_{\text{room}}}^{2})}{2} + a_{2} \times \frac{({T_{\text{fluid}}}^{3} - {T_{\text{room}}}^{3})}{3}\rbrack\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{2.4}) \end{split}\]
Questions
(1) Find the contraction in the outer diameter of the trunnion.
(2) As required, is the magnitude of contraction more than \(0.015{^{\prime\prime} }\)?
(3) If the magnitude of the contraction is less than \(0.015{^{\prime\prime} }\), what if the trunnion were immersed in liquid nitrogen (boiling temperature\(= - 321{^\circ}F\))? Will that give enough contraction in the trunnion?
(4) Rather than regressing the data to a second-order polynomial so that one can find the contraction in the outer diameter of the trunnion, how would you use the trapezoidal rule of integration for unequal segments to find the contraction? What is the relative difference between the two results? The data for the coefficient of linear thermal expansion as a function of temperature is given in Table \(\PageIndex{2.1}\).
(5) We chose a second-order polynomial for regression. What would make a different order polynomial a better choice for regression? Is there an optimum order of polynomial you can find?
Temperature | Coefficient of Thermal Expansion |
---|---|
\({^\circ}F\) | \(\mu \text{in} / \text{in}^{\circ}\text{F}\) |
\(80\) | \(6.47\) |
\(60\) | \(6.36\) |
\(40\) | \(6.24\) |
\(20\) | \(6.12\) |
\(0\) | \(6.00\) |
\(-20\) | \(5.86\) |
\(-40\) | \(5.72\) |
\(-60\) | \(5.58\) |
\(-80\) | \(5.43\) |
\(-100\) | \(5.28\) |
\(-120\) | \(5.09\) |
\(-140\) | \(4.91\) |
\(-160\) | \(4.72\) |
\(-180\) | \(4.52\) |
\(-200\) | \(4.30\) |
\(-220\) | \(4.08\) |
\(-240\) | \(3.83\) |
\(-260\) | \(3.58\) |
\(-280\) | \(3.33\) |
\(-300\) | \(3.07\) |
\(-320\) | \(2.76\) |
\(-340\) | \(2.45\) |