Processing math: 100%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Mathematics LibreTexts

4.00: Physical Problems for Simultaneous Linear Equations

( \newcommand{\kernel}{\mathrm{null}\,}\)

Lesson 1: General Engineering

Summary

To find the velocity profile of a rocket requires a solution of simultaneous linear equations.

Modeling

The upward velocity of a rocket is given at three different times in the following table

Time, t (s) Velocity, v (m/s)
5 106.8
8 177.2
12 279.2
A rocket launching into space.
Figure 4.0.1.1. A rocket launching into space.

The velocity data is approximated by a polynomial as

v(t)=at2+bt+c, 5t12.(4.0.1.1)

Set up the equations in matrix form to find the coefficients a, b, c of the velocity profile.

The polynomial is going through three data points (t1,v1),(t2,v2), and (t3,v3) where from the above table

t1=5,v1=106.8

t2=8,v2=177.2

t3=12,v3=279.2

Requiring that v(t)=at2+bt+c passes through the three data points gives

v(t1)=v1=at21+bt1+cv(t2)=v2=at22+bt2+cv(t3)=v3=at23+bt3+c

Substituting the data (t1, v1), (t2, v2), (t3, v3) gives

a(52)+b(5)+c=106.8a(82)+b(8)+c=177.2a(122)+b(12)+c=279.2

or

25a+5b+c=106.864a+8b+c=177.2144a+12b+c=279.2

This set of equations can be rewritten in the matrix form as

[25a+5b+c64a+8b+c144a+12b+c]=[106.8177.2279.2]

The above equation can be written as a linear combination as follows

a[2564144]+b[5812]+c[111]=[106.8177.2279.2]

and further using matrix multiplications gives

[25516481144121] [abc]=[106.8177.2279.2]

The solution of the above three simultaneous linear equations will give the value of a,b,c.

Questions

(1) Solve for the values of a,b,c.

(2) Verify if you get back the value of the velocity data at t=5 s.

(3) Estimate the velocity of the rocket at t=7.5 s?

(4) Estimate the acceleration of the rocket at t=7.5 s?

(5) Estimate the distance covered by the rocket between t=5.5 s and 8.9 s.

(6) If the following data is given for the velocity of the rocket as a function of time, and you are asked to use a quadratic polynomial to approximate the velocity profile to find the velocity at t=16 s, what data points would you choose and why?

\(t \ (\text{s})\) \(v(t) \ (\text{m/s})\)
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67

Lesson 2: Civil Engineering

Summary

A pressure vessel made of a compounded cylinder has a higher capability to handle internal pressure over a single cylinder. To find this capability, one has to solve a set of simultaneous linear equations.

Modeling

A pressure vessel can only be subjected to an amount of internal pressure that is limited by the strength of the material used. For example, take a pressure vessel of an internal radius of a=5 and outer radius b=8, made of ASTM 36 steel (yield strength of ASTM 36 steel is 36 ksi). How much internal pressure can this pressure vessel take before it is considered to have failed?

Top view of a single-cylinder pressure vessel, showing the internal radius a and the outer radius b.
Figure 4.0.2.1. A single-cylinder pressure vessel with internal radius a and outer radius b.

The hoop and radial stress in a cylindrical pressure vessel is given by Equations (4.0.2.1) and (4.0.2.2).

σr=a2pib2a2(1b2r2)(4.0.2.1)

σθ=a2pib2a2(1+b2r2)(4.0.2.2)

The maximum normal stress anywhere in the cylinder is the hoop stress at the inner radius, a:

σθ|max=pi(b2+a2b2a2)(4.0.2.3)

Assuming a factor of safety of 2, while the yield strength is given as 36 ksi,

36×1032=pi(82+528252)

pi=7.887 ksi

You can see from Equation (4.0.2.3) that even for b>>a, the maximum internal pressure one can apply is only pi=18ksi. Therefore, what can an engineer do to maximize the internal pressure while keeping the material and radial dimensions the same? They can use a compounded cylinder. They can create a compounded cylinder by shrink-fitting one cylinder into another, and hence creating pre-existing favorable stresses to allow more internal pressure. Let us see how that would work.

Top view of a compounded cylinder pressure vessel showing the internal radius a, the outer radius b, and the interface at r=c.
Figure 4.0.2.2. A compounded cylinder pressure vessel with internal radius a, outer radius b, and interface at r=c.

Let us make the compounded cylinder of two cylinders (Figure 4.0.2.2). Cylinder 1 has an internal radius of a=5, and outer radius c=6.5, while Cylinder 2 has an internal radius of c=6.5 and outer radius, b=8. Assume that that radial interference, δ=0.007 occurs at the interface of a compounded cylinder at r=c=6.5. How does then one find the pressure that can be applied to the compounded cylinder of internal radius, a=5 and outer radius, b=8?

For cylinder 1, the radial displacement u1 is given by

u1=c1r+c2r(4.0.2.4)

and the radial stress σ1r and hoop stress σ1θ are given by

σ1r=E1ν2[c1(1+ν)c2(1νr2)](4.0.2.5)

σ1θ=E1ν2[c1(1+ν)+c2(1νr2)](4.0.2.6)

where

E=Young's modulus of steel,

ν=Poisson's ratio of steel.

For cylinder 2, the radial displacement u2 is given by

u2=c3r+c4r(4.0.2.7)

and the radial stress σ2r and hoop stress σ2θ are given by

σ2r=E1ν2[c3(1+ν)c4(1νr2)](4.0.2.8)

σ2θ=E1ν2[c3(1+ν)+c4(1νr2)](4.0.2.9)

So if one is able to find the four constants, c1,c2,c3 and c4, one can find the stresses in the compounded cylinder to be able to find what internal pressure can be applied. So how do we find the four unknown constants?

The boundary and interface conditions are the following.

The radial stress at the inner radius, r=a, is the applied internal pressure

σ1r(r=a)=pi(4.0.2.10)

The radial stress is continuous at the interface, r=c:

σ1r(r=c)=σ2r(r=c)(4.0.2.11)

The radial displacement at the interface r=c has a jump of the radial interference, δ

u2(r=c)u1(r=c)=δ(4.0.2.12)

The radial stress at the outer radius r=b is

σ2r(r=b)=0(4.0.2.13)

This will set up four equations and four unknowns, if we know what internal pressure we are applying. Assume we are applying the same pressure as the single cylinder can take, that is, pi=7.887 ksi and let us see later what stresses it creates in the compounded cylinder.

Assuming E=30×106 psi, ν=0.3, Equations (4.0.2.10) through (4.0.2.13) become

30×10610.32[c1(1+0.3)c2(10.352)]=7.887×103

30×10610.32[c1(1+0.3)c2(10.36.52)]=30×10610.32[c3(1+0.3)c4(10.36.52)]

c3(6.5)+c46.5c1(6.5)c26.5=0.007

30×10610.32[c3(1+0.3)c4(10.382)]=0(4.0.2.14ad)

Writing the Equations (4.0.2.14ad) in matrix form, we get

[4.2857×1079.2307×105004.2857×1075.4619×1054.2857×1075.4619×1056.50.153846.50.15384004.2857×1073.6057×105][c1c2c3c4]=[7.887×10300.0070](4.0.2.15)

Solving these four simultaneous linear equations, we can find the four constants.

References

(1) A.C. Ugural, S.K. Fenster, Advanced strength and applied elasticity, Third Edition, Prentice Hall, New York, 1995.

(2) J.E. Shigley, C.R. Mischke, Chapter 19 - Limits and fits, Standard handbook of machine design, McGraw-Hill, New York, 1986.

Questions

(1) Find the unknown constants of Equation (4.0.2.14) using different numerical methods.

(2) Knowing that the critical points in the compounded cylinder are r=a,c,c+, and b, find the maximum hoop stress in the compounded cylinder. What is its value compared to the maximum hoop stress allowable of 18 ksi?

(3) Find the maximum internal pressure you can apply to the compounded cylinder. Compare it with the maximum possible internal pressure for a single cylinder of the same dimensions.

(4) The radial interference at the interface is created by making the inner cylinder 1 to have a larger outer radius than the inner radius of cylinder 2. Standard interference fits dictate the limits of these dimensions. If cylinder 2 is fit into cylinder 1, there is an upper and lower limit by which the nominal diameter of each cylinder varies at the interface. This limit L, in thousands of an inch, is given by

L=CD1/3 where D (nominal diameter) is in inches, and the coefficient C, based on the type of fit, is given in Table 1 below.

Table 4.0.2.1. Limits of cylindrical interference fits.
Cylinder Limit Class of fit
    FN2 FN3
1 Lower 0.000 0.000
  Upper 0.907 0.907
2 Lower 2.717 3.739
  Upper 3.288 4.310

Assuming FN2 fit at the interface, find the maximum internal pressure you would recommend.

Lesson 3: Mechanical Engineering

Summary

To find the diametrical contraction during shrink-fitting in a trunnion requires one to first find how the coefficient of linear thermal expansion of steel is related to temperature. The relation is given through a second-order polynomial, and calculating the coefficients of the polynomial requires solving a set of simultaneous linear equations.

Modeling

To make the fulcrum (Figure 4.0.3.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 shrunk-fit into the girder of the bridge.

A bascule bridge.
A trunnion fitted into a hub, which is in turn fitted into a bridge girder.
Figure 4.0.3.1. Trunnion-Hub-Girder (THG) assembly.

This is done by first immersing the trunnion in a cold medium such as 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 4.0.3.2).

A trunnion being slid through a hub after contracting.
Figure 4.0.2. Trunnion is slid through the hub after contracting.

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 be ordered at a 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 was because the trunnion had not contracted enough to slide through the hole. Can you find out why?

A hollow trunnion of an outside diameter 12.363 is to be fitted in a hub of an inner diameter 12.358. The trunnion was put in a dry-ice/alcohol mixture (temperature of the fluid - dry ice/alcohol mixture is 108F) 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 is required between the trunnion and the hub. Assuming the room temperature is 80F, is immersing it in dry-ice/alcohol mixture a correct decision?

Solution

To calculate the contraction in the diameter of the trunnion, the coefficient of linear thermal expansion at room temperature is used. In that case, the reduction ΔD in the outer diameter of the trunnion is

ΔD=DαΔT(4.0.3.1)

where

D=outer diameter of the trunnion,

α=coefficient of linear thermal expansion at room temperature, and

ΔT=change in temperature,

Given

D=12.363

α=6.817×106in/in/F at 80F

ΔT=TfluidTroom=10880=188F where

Tfluid=temperature of dry-ice/alcohol mixture

Troom=room temperature

the reduction in the trunnion outer diameter is given by

ΔD=(12.363)(6.47×106)(188)=0.01504

So, the trunnion is predicted to reduce in diameter by 0.01504. But is this enough reduction in diameter? As per specifications, we need the trunnion to contract by

=trunnion outside diameterhub inner diameter+diametral clearance=12.36312.358+0.01=0.015

So, according to these calculations, immersing the steel trunnion in dry-ice/alcohol mixture gives the desired contraction of 0.015 as we predict a contraction of 0.01504.

But as shown in Figure 4.0.3.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, Equation (4.0.3.1) would overestimate the thermal contraction.

Graph of the varying coefficient of linear thermal expansion as a function of temperature for cast steel.
Figure 4.0.3.3. Varying coefficient of linear thermal expansion as a function of temperature for cast steel.

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

ΔD=DTfluidTroomαdT(4.0.3.2)

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 polynomial through the data given in Table 4.0.3.1.

Table 4.0.3.1. Coefficient of linear thermal expansion as a function of temperature.
Temperature Coefficient of linear thermal expansion
F μin/in/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
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

Assuming that the coefficient of linear thermal expansion is related to temperature by a second-order polynomial,

α=a0+a1T+a2T2(4.0.3.3)

Given the data points (α1,T1), (α2,T2), …..,(αn,Tn) as in Figure 4.0.3.3 and Table 4.0.3.1, the sum of the square of the residuals (sum of the square of the differences between the observed and predicted values) is Sr=ni=1(αi{a0+a1Ti+a2T2i})2=ni=1(αia0a1Tia2T2i)2(4.0.3.4)

To minimize the value of the sum of the square of the residuals, we take the derivative with respect to each of the three unknown coefficients, a0, a1, and a2, to give

Sra0=ni=12(αia0a1Tia2T2i)(1)=2[ni=1αi+na0+a1ni=1Ti+a2ni=1T2i]

Sra1=ni=12(αia0a1Tia2T2i)(Ti)=2[ni=1αiTi+a0ni=1Ti+a1ni=1Ti2+a2ni=1T3i]

Sra2=ni=12(αia0a1Tia2T2i)(T2i)=2[ni=1αiTi2+a0ni=1Ti2+a1ni=1T3i+a2ni=1T4i]   (4.0.3.5)

Setting the three partial derivatives in Equation (4.0.3.5) equal to zero gives,

na0+a1ni=1Ti+a2ni=1T2i=ni=1αi(4.0.3.6a)

a0ni=1Ti+a1ni=1T2i+a2ni=1T3i=ni=1αiTi(4.0.3.6b)

a0ni=1T2i+a1ni=1T3i+a2ni=1T4i=ni=1αiT2i(4.0.3.6c)

The set of equations given by equations 4.0.3.6a-c are simultaneous linear equations. The number of data points in Figure 4.0.3.3 is 24, as given in Table 4.0.3.1. Hence

n=24

24i=1Ti=2860

24i=1T2i=7.26×105

24i=1T3i=1.86472×108

24i=1T4i=5.24357×1010

24i=1αi=1.057×104

24i=1αiTi=1.04162×102

24i=1αiT2i=2.56799

24a02860a1+7.26×105a2=1.057×1042860a0+7.26×105a11.8647×108a2=1.04162×1027.26×105a01.86472×108a1+5.24357×1010a2=2.56799(4.0.3.7)

In matrix form, the three simultaneous linear equations can be written as

[2428607.26×10528607.26×1051.86472×1087.26×1051.86472×1085.24357×1010][a0a1a2]=[1.057×1041.04162×1022.56799](4.0.3.8)

Questions

(1) Can you now find the contraction in the trunnion outer diameter?

(2) Is the magnitude of contraction more than 0.015 as required?

(3) If that is not the case, what if the trunnion were immersed in liquid nitrogen (boiling temperature=321F)? Will that give enough contraction in the trunnion?

(4) Redo problem #1 using a third-order polynomial as the regression model. How much different is the estimate of contraction using the third-order polynomial?

(5) Find the optimum polynomial order to use for the regression model.

(6) Find the effect of the number of significant digits used in solving the set of the equations for problem #4, as you must have noticed that the order of the numbers in the coefficient matrix varies quite a bit.


This page titled 4.00: Physical Problems for Simultaneous Linear Equations is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Autar Kaw via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?