Skip to content

Commit 198450f

Browse files
authored
Update README.md
1 parent 05c0b47 commit 198450f

File tree

1 file changed

+0
-33
lines changed

1 file changed

+0
-33
lines changed

README.md

-33
Original file line numberDiff line numberDiff line change
@@ -22,36 +22,3 @@ Learn more about the *Python for Mechanical and Aerospace Engineering* book [her
2222

2323
<center><a rel="license" href="http://creativecommons.org/licenses/by-nc-sa/4.0/"><img alt="Creative Commons License" style="border-width:0" src="https://i.creativecommons.org/l/by-nc-sa/4.0/88x31.png" /></a></center><br />
2424
This work is licensed under a <a rel="license" href="http://creativecommons.org/licenses/by-nc-sa/4.0/">Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License</a>.
25-
26-
## Correction
27-
28-
Turns out, I was right to be skeptical of 1 billion psf as the maximum dynamic pressure in Chapter 3! I mis-interpreted the NASA atmospheric model and treated pressure as density.
29-
This corrected density function should give you the original Max Q psf average acceleration of 1,395 psf:
30-
31-
def density_corrected(height: float) -> float:
32-
"""
33-
Returns the air density in slug/ft^3 based on altitude
34-
Equations from https://www.grc.nasa.gov/www/k-12/rocket/atmos.html
35-
:param height: Altitude in feet
36-
:return: Density in slugs/ft^3
37-
"""
38-
if height < 36152.0:
39-
temp = 59 - 0.00356 * height
40-
p = 2116 * ((temp + 459.7)/518.6)**5.256
41-
elif 36152 <= height < 82345:
42-
temp = -70
43-
p = 473.1*np.exp(1.73 - 0.000048*height)
44-
else:
45-
temp = -205.05 + 0.00164 * height
46-
p = 51.97*((temp + 459.7)/389.98)**-11.388
47-
48-
rho = p/(1718*(temp + 459.7))
49-
return rho
50-
51-
You will have to adjust the plotting offsets in `chap3.py` to correct for the changed scale and magnitude:
52-
53-
> xytext=(x_values[ind] + 15, ~~max_val + 1E5~~),
54-
55-
should be
56-
57-
> xytext=(x_values[ind] + 15, **max_val + 15**),

0 commit comments

Comments
 (0)