Richard Nakka’s Experimental Rocketry Web Site
________________________________________________
Introduction
to Rocket Design
Appendix K
Parachute Sizing and Loads
Introduction
My webpage Parachute Design and Construction provides a good introduction to the aerodynamics of parachutes,
the basics of parachute design and testing, as well as means of constructing
one’s own parachute. This appendix is intended to serve more as a design aid,
providing details on sizing a parachute based on desired performance and expected
loading conditions.
A bit of personal history. My first parachute equipped EX rocket
was my “B”
rocket that I first launched
in April of 1972. Purchasing a suitable parachute back then was not an option. Consequently,
I designed and made my own parachute using nylon fabric purchased at the local KMart. I fabricated my parachute with individual
gores, alternating red and white fabric (perhaps not surprisingly, my
inspiration was the Apollo
parachute). The gores were sewn
together using my mother’s sewing machine (with her guidance and approval, of
course). Later on, I fortuitously came across some military surplus parachutes
at the local Army
Surplus store which were
perfect for my rockets. Figure 1 illustrates a 33² circular chute (was white, I dyed it red) and a
cross parachute (9²W ´ 28²L), both ripstop nylon. When that surplus supply
dried up, and I needed a parachute for my Frostfire rocket, I decided to make a simple cross parachute. In more recent times, now that commercially
made parachutes for rocketry have become available (and money not being as
tight as it was), I purchase my parachutes. Fruity Chutes, for example, has a large selection of
parachutes of all sizes and of different types, and the quality of their chutes
is impeccable.
Figure 1: Two of my parachutes purchased at a military
surplus store
Overview
Parachutes are, without a doubt, the ideal Recovery Descent system
for rockets. A parachute can return a rocket to the ground more gently than any
other recovery method. Even now in the 21st century, parachutes remain
the means of choice for recovery of commercial spacecraft such as Orion, Dragon
and Starliner. The reasons are simple. Parachutes can be tailored to safely
recover payloads of any size, from eggs to expended Space Shuttle Solid Rocket Boosters,
weighing in at 175,000 pounds each. A parachute is lightweight, can be packed
into a remarkably tiny volume, is generally inexpensive and durable. The desired
rate of descent, or softness of landing, is simply a matter of choosing the
appropriate size. The bigger the chute, the softer the landing. The principle
of operation is simple making a parachute highly reliable. Figure2 illustrates
the main parts of a typical EX parachute.
Figure 2: Parts of a typical parachute
When inflated, a parachute has a huge drag-to-mass ratio. Although
a parachute is always used to descend some kind of payload, it is interesting
to examine the Ballistic
Coefficient of a parachute alone. The
Ballistic Coefficient (see my webpage Aerodynamics and its Role in Experimental Rocket Design ) is given by:
where
m = mass
CD = drag coefficient
A = cross-sectional area (normal to velocity
vector)
As an example, the 42² (1067mm) ellipsoidal parachute I use for
recovery of my Xi rocket has a mass of 120 grams (including shroud lines and
bridle). The nominal cross-sectional, or projected, area is 894,000 mm2.
The CD is 1.5 based on projected
area.
As such, for this parachute:
A comparison of BC for
different objects is given in this table. A ping-pong ball is 55´ as slippery as
this parachute and the Xi rocket is 11,500´ as slippery ! (the more slippery an object, the less drag resistance it offers). As
such, it is clear to see that attaching a parachute to a payload works well to
slow its descent through the atmosphere.
Types of Parachutes
Parachutes, or rather the
parachute canopy, come in a vast assortment of
shapes and configurations. More common shapes are flat circular, conical,
bi-conical, hemispherical, ellipsoidal, cross, annular, toroidal, ringslot,
ringsail and ribbon-type. Ribbon-type or slotted chutes have particularly low
opening forces. Parachutes can be employed singly or clustered. This webpage
will focus on the three most common parachute types used for EX rocketry:
-
Parasheet or flat
circular
-
Elliptical or ellipsoidal
-
Cruciform or cross
I have made and used all
three types over the years. The parasheet is simplest to make, it is basically
a single piece of fabric cut to a circular shape. Shroud lines are sewn onto
the hemmed circumference. This type of parachute is suited to model rockets and
relatively lightweight LoPer class rockets. The Ellipsoidal and Cross type parachute are better
performing (have a higher drag) and can be made to be robust enough to handle
any weight rocket and to withstand greater opening velocity than a simple parasheet.
The three canopy types are illustrated in Figure 3. The cross parachute is a
good choice for deployment at high speed as it has a low opening shock factor.
The cross parachute also performs well at supersonic velocities.
The shape of a parachute
canopy when inflated is determined by a balance of internal pressure forces and
tension in the shroud lines.
Figure 3: Canopy shape for various type parachutes
Notes: [1]
Referenced to canopy surface area
[2]
According to Res.K2, the drag coefficient greatly increases at low velocity
Drag Force of Parachutes
Parachute drag force is a
key parameter, as it determines the rate of descent. Parachute drag will equal
to the weight of the total assembly (rocket, tether and parachute). Parachute
drag force, Fc, is given by:
CD = parachute drag coefficient
S = reference area of canopy associated with CD (ft2
or m2)
q = ½ r v2 =
dynamic pressure due to descent velocity, v and r = local air density (slug/ft3 or kg/m3)
Units of drag, Fc, are either pounds-force (lbf) or Newtons (N).
Unlike streamlined bodies, Reynolds
number does not change the
drag coefficient of parachutes. This is because bodies with separated flow,
such as parachute canopies, are little affected by Reynolds number.
As mentioned earlier, parachutes have a huge drag-to-mass ratio.
Using the same example as earlier, the 42² (1067mm) ellipsoidal parachute that I use for
recovery of my Xi rocket has a mass of 120 grams (including shroud lines and
bridle). Drag coefficient is 1.5 (based on projected area). Descent rate for a
3 kg. Xi rocket is 20 ft/sec. or 6.1 m/sec. The parachute
drag force is (using m.k.s system of
units):
D = 1.5 π/4 (1.067)2 [½ (1.225)
(6.1)2 ] = 30.6 N.
Giving a drag-to-mass ratio of 255 Newtons per kilogram, or in US
units 26 pound-force per pound-mass.
Rate of Descent
A stable parachute in
descent is in equilibrium between the total drag of the parachute and the
weight of the payload. This is illustrated in Figure 4.
Figure 4: Forces acting on parachute in stable
descent
The vertical descent rate (velocity) provided by a parachute in a
stable descent is given by:
where Wt = total weight of body including parachute (kg-m/sec2.
or lbm.), S = canopy
reference area (m2 or ft2), and r = local air density (kg/m3 or slug/ft3).
The descent rate units are metres per
second (m/s) or feet per second (fps).
Things to note about this equation.
1)
Wt is weight, not mass. As such, if using
SI system of units, mass must be multiplied by g,
the acceleration of gravity (9.808 m/sec). i.e. Wt
= mt g
2)
S is the reference area of the canopy upon which the drag
coefficient, CD, is
based. Unfortunately, there exists inconsistency in definition of canopy area.
Three possible reference areas are:
a.
Nominal area. This is the cross-sectional area of a non-inflated
canopy. For a circular parasheet, this is straightforward to calculate.
However, for a shaped canopy (that does not lie flat) such as ellipsoidal,
there is clearly some ambiguity.
b.
Projected area. This is the cross-sectional projected
area of the inflated parachute
canopy. This is similar in concept to the reference area for solid bodies such
as nose cones or rocket bodies. Manufacturers prefer this definition as it
gives the highest numerical value of CD. The problem with this reference area is that
it is not constant. It changes with descent speed and shroud line length. As
such, CD referenced to parachute
projected area should be specified in relation to inflated diameter (or other
defining parameter), shroud line length and descent velocity.
c.
Surface area of the canopy (So).
This is the area of the actual fabric of the canopy. Engineers (such as T.W.Knacke) prefer this definition as the surface area of
the canopy is fixed and as such in
non-ambiguous. Vent holes and slots are included in this area, by definition.
The drag coefficient based on surface (fabric) area indicates how effectively a
parachute canopy produces drag with a given amount of fabric area. As such,
drag coefficient referenced to canopy surface
area (CDo) provides an excellent
measure of parachute efficiency (drag force per canopy mass).
3)
Air density (r) decays with increasing altitude. Therefore
descent velocity gradually lessens as the descending rocket gets closer to the
ground, at a rate inversely proportional to the square root of r. This is one small, but not insigificant, way
that nature is on the rocketeer’s side. The descent velocity (ve) at any altitude is given by:
where veo
is the descent velocity at sea-level, ro is the air density at sea-level and r is the air density at the altitude of interest.
For a given parachute
shape, drag coefficient can be increased by the use of longer shroud lines,
greater number of shroud lines, and lower canopy porosity. Longer shroud lines
increase the inflated diameter of the canopy, thereby increasing drag. This
holds true for both circular cross-section and cross parachutes. My 42” ellipsoidal
Fruity Chutes parachute has twelve
shroud lines of 1.15´ chute diameter. Figure 5 shows the effect of
shroud line length on the drag coefficient for a 40” cross parachute (taken from Res.K2). The
increase in drag coefficient with shroud line length is quite significant. Note
that the drag coefficient (CDo) for this example is referenced to canopy surface area (So = 2 L W – W 2).
Figure 5: Drag
coefficient as function of shroud line length
It is worth noting that a higher drag coefficient does not necessarily mean a parachute with
superior drag characteristics. It is the
drag force that is the key
measurement when comparing parachute shapes. For a given rate of descent,
dynamic pressure is a fixed value. As such, it is the drag area
(CDo So ) that
counts when comparing parachutes. Or more importantly, for weight critical
applications such as a rocket, the drag force per unit mass of the parachute
canopy. My 42” ellipsoidal Fruity Chutes parachute has a published CD
= 1.5. This is presumably based on canopy projected area.
Lower canopy porosity
has the drawback that it can make a parachute less stable with a tendency to
oscillate. Figure 6 illustrates the airflow around a porous and a non-porous
canopy. As airflow cannot freely go through a non-porous canopy, the flow goes
around and separates on the leading edge in alternating vortices, forming a Karman vortex trail. The separation results in alternate pressure areas on opposite
sides of the canopy, generating an alternating force perpendicular to the
airflow velocity vector, thereby generating oscillations.
Figure 6: Comparison between behaviour of porous vs
non-porous canopy
In a similar manner to
increasing porosity, adding a spill hole serves to reduce or eliminate such
oscillation. The hole acts as a vent to release air pressure buildup without
causing uneven air flow on one side of the canopy, and in doing so, stabilizes
the descent by allowing trapped air to escape evenly, preventing the parachute
from rocking or oscillating wildly as it falls, ensuring a more controlled and
straight descent. A spill hole allows for use of a lower-porosity fabric. For
example, my 42” ellipsoidal Fruity Chutes
parachute has a spill hole that is 20% of the chute diameter. This allows for
use of very low-porosity (3 ft3/ft2/min.) 1.2 oz. fabric
(PIA-C 44378 type IV). For cross parachutes, such low canopy porosity results
in rotational instability.
By rule of thumb, a
spill hole area should be 1% of the canopy area. A spill hole is not necessary
for a cross parachute.
Sizing the Parachute
The first step is to choose a suitable descent rate. For my
rockets equipped with a drogue chute, descent from apogee to a lower altitude,
a descent velocity in the range of 70-80 feet per second (21- 24 m/s), based on
ground level air density, was considered appropriate. This has been considered
as a good compromise between getting the rocket down quickly (to minimize
downrange drift due to wind) and minimizing damage to the rocket in case the
main chute fails to fully deploy. For descent under the main parachute, a
suitable descent rate for my rockets has been in the range of 20-30 feet per
second (6 – 9 m/s). The hardness of the terrain that the rocket will touchdown
on is a factor. Less hard terrain such as tall grass or snow allows for a higher
touchdown speed.
Knowing the total weight (Wt)
of our descending rocket, and having settled on a descent speed (v), the surface area of
the canopy (So) can be determined:
In this equation air density (r) should be that at ground level,
whether sizing a drogue or main parachute. We know that the rocket with drogue
chute will descend much faster at high altitudes, as air density drops with
altitude. However, we wish to size the drogue chute based on descent velocity
when either the main deploys (near the ground) or failing that, when the descending
rocket reaches the ground.
It is important to note that the drag coefficient (CDo) in equation 3 must
be related to the canopy surface area as the reference area. If the drag
coefficient is referenced to the canopy projected (or
inflated) area, the equation is similar:
In this case, the canopy area, S, is taken as the
projected area, as shown in Figure 7.
Figure 7: Canopy surface area versus projected area
Once the canopy surface
area (So) has been determined, the
canopy projected diameter (Dp) can
be taken as the following:
Flat Circular diameter
Ellipsoidal base diameter
for b/a = 0.707, where a and b are the major and minor semiaxes, respectively
Cross panel length
where l = W/L
For a circular or
ellipsoidal chute whereby the projected area has been calculated using equation 4:
Note that in all cases units must be consistent. For example, area
in square feet gives diameter (or length) in feet.
Time of Descent
It is useful to know how
long it would take for a rocket to descend from one altitude level (z1) to lower altitude level (z2). The time duration of
descent can, for example, be used to estimate the downrange drift of a rocket,
based on wind speed. Descent time, td,
is given by:
assuming that the rocket descends at local terminal velocity (vd). Using the drag force
equation (Eqn.1), terminal velocity is given by
Atmospheric air density varies significantly with altitude. Using equation 9 and equation 10
together with polynomial curve fitting of the air density data results in the
following expression for descent time between any two altitude points in the
range of zero to 100,000 feet (30 km)
above sea level (ASL):
with constant C1
given by:
The coefficients A,
B, and C
are as follows:
US units of measure (slugs, feet): SI units of measure (m.k.s)
A = 1.057 e-12 2.582e-10
B = -3.756 e-7 -2.798e-5
C = 0.0490 1.112
(based on density lapse rate per International Standard Atmosphere)
Descent time calculated by this method is
accurate to within approximately two percent based on the fidelity of the
second-order curve fit. Details of the method are described here.
Canopy Opening
Parachute canopy inflation is a dynamic and chaotic (non-linear)
process. The canopy expansion is resisted by tension in the shroud lines until
fully inflated. The canopy inflation process is illustrated in Figure 8. The
initial steps (a,b &c) occur relatively slowly, whereas the final steps
(e,f &g) occur in rapid succession.
Figure 8: Parachute canopy inflation process
a)
Opening of canopy mouth
b)
Air mass moves into canopy
c)
Ball of air reaches crown of canopy (less pronounced if spill hole
present)
d)
Influx of air expands the canopy
e)
Expansion of canopy resisted by shroud line tension
f)
Canopy becomes fully inflated
g)
Canopy over-inflates; crown depressed by moment of surrounding air
I have had the good
fortune of witnessing my fair share of parachute deployments and, as such, have
observed this distictive sequence. In particular, how a parachute canopy opens
slowly at first, then rapidly inflates, overinflates, then settles to its
equilibrium shape. One additional behaviour has been observed. Following (g),
it has been noticed that the canopy repeatedly overinflates then underinflates,
doing this three or four times. This is strikingly visible in the video taken
of the Xi-22 parachute (36²) deployment.
Xi-22 Parachute Deployment videoclip ( ¼ actual speed) 15 Mbytes
This cycle of
overinflation/underinflation is certainly due to varying tension in the shroud
lines. The inertia term in equation 12,
acting in collaboration with the elastic nature of the parachute assembly, behaves
as a spring-mass
system.
The force (F) acting (at
the bridle) of a vertically descending parachute plus rocket, as a function of
time, is given by:
where:
CD S = instantaneous drag area of the canopy (ft2 or m2),
which ranges from 0% to 100%
r = air density
(slug/ft3 or
kg/m3)
mr = rocket mass (slug or kg)
mp = parachute mass (slug or kg)
ma = apparent mass (slug or kg)
The apparent mass is that of the air
captured by the canopy. For most EX rockets, this term is much less than the
mass of the rocket, and may be omitted. The middle term in equation 12
is the inertial force term due to the system mass and change in velocity of the
rocket as it decelerates from free-fall (or drogue) descent velocity to parachute
descent velocity .
The drag area (CD S),
as a function of time, can be taken as (see Res.K7):
where b = 6 for a regular (non-slotted) canopy, and tf
is the time to fully open. This relationship is plotted in Figure 9,
showing percent of drag area versus time. As can be seen, the canopy opens at a
slow rate initially, then suddenly blossoms to
fully inflated state. How rapidly this occurs is, of course, dependent upon the
size of the canopy, as well as other factors. When the Xi parachute (36”
ellipsoidal) opens cleanly, this can occur is less than one second. It is this
rapid opening that generates shock or opening force
loading of the parachute.
Figure 9: Canopy
drag area % vs time to fill ratio
It can be useful to know
how long it takes for a canopy to fully inflate. When a rocket is descending,
for example at free-fall velocity such as 85 feet per second (26 m/s.), each
second the rocket gets 85 feet closer to the ground. The fill time needs to be
considered when choosing a deployment altitude. From Res.K1,
the canopy inflation time for a flat circular parachute may be expressed as:
where n = 4.0 for
standard porosity and n = 2.5 for
low-porosity canopies.
Do is the nominal diameter given by:
For example, for the 36² ellipsoidal parachute for Xi-22 rocket, the
canopy area So = 11.5 ft2 , giving Do =
3.82 ft. The free-fall velocity for this flight was 80.3 ft/sec., giving a
canopy fill time of:
tf
= 4.0 (11.5)/(80.3)0.85
= 1.1 second.
This assumes that the fill time for an ellipsoidal parachute is
the same as for a flat circular chute. The value of 1.1 second is close to the
fill time measured from the video, which is slightly less than one second. As
such, the rocket descends approximately 80 feet from the time the chute is
deployed to the time it fully inflates.
The deployment of a parachute generates two distinct and separate
forces to the suspended load. These are known as the snatch force
(Fs) and the opening
force (Fx). Snatch force is the maximum shroud line tension during
deployment prior to canopy inflation. When the parachute is deployed, the
parachute and attached rocket must come to a common velocity. The inertial
force of the rocket generates tension in the bridle when the shroud lines become
taut. The canopy then inflates, followed by a brief overinflation, resulting in
a peak momentary force that is transmitted to the rocket by the parachute. This
is known as the opening force. The attachment force then decays to that of the
equilibrium drag force (Fc). This
sequence is illustrated in Figure 10.
Figure 10: Forces
developed by a deployed and inflating parachute
Based on accelerometer
data from many flights, it is apparent that the more critical of the two is the
opening force (as per Fig.10). The equation for the opening force, taken from Res.K1, is:
Cx is the opening force coefficient assuming an infinite mass payload. This is the case where the mass of
the payload is large compared to the drag force of the parachute. Values for Cx are:
Flat circular 1.7
Ellipsoidal 1.6
Cross 1.2
The force reduction
factor, X1, provides alleviation for
the case where the payload mass is finite, or in
the same order of magnitude (or so) as the parachute mass. This would be the
case for EX rockets. The force reduction factor may be calculated using this graph with ballistic parameter, A, defined as:
For a flat circular,
ellipsoidal, or cross parachute, the value of n = 2 should be assumed when interpreting
the graph. The parameter is dimensionless, and as such, the following units are
applicable:
Wt = total weight of body including
parachute (kg-m/sec2. or lbm.), S = canopy reference area (m2 or ft2), r = local air density (kg/m3 or slug/ft3), g = acceleration of gravity (9.81 m/sec2
or 32.2 ft/sec2), v1 = velocity at chute
deploy (m/sec or ft/sec) and tf = inflation time
(sec.).
Note that the drag area term (CD S) must be consistent. If CD
is referenced to the surface area of the canopy, then S =
So. If CD
is referenced to the projected area of the canopy, then S
is the projected area shown in Figure 7.
For my Xi rocket fitted with 42” ellipsoidal parachute, the force
reduction factor works out to X1
= 0.21. As such, the opening
force comes to 31 lbf. (138 N.).
For design purposes, it is undoubtedly wise to assume the infinite mass condition and as such, conservatively assume X1 = 1.
Resources
Res.K1 NWC
TP 6575 Parachute
Recovery Systems Design Manual, T.W.Knacke, Para Publishing
Res.K2 NOLTR
71-111 Effects of Canopy Geometry
on the Drag Coefficient of a Cross Parachute in the Fully Open and Reefed
Conditions W.P.Ludtke
Res.K3 MIL-C-7020 (Military
Specification) Cloth, Parachute, Nylon-Rip Stop and Twill
Weave
Res.K4 MIL-C-44378 (Military
Specification) Cloth, Parachute, Nylon, Low Permeability
(this spec is same as PIA-C-44378)
Res.K5 Parachute
Canopy Fabric Mil Specs
Res.K6 FDL-TDR-64-155 Drag and Stability of Cross Type Parachutes, R.J.Niccum,
E.L.Haak, & R.Gutenkauf
Res.K7AGARD-AG-295 The Aerodynamics of Parachutes D.J.Cockrell
Last updated January 23,
2025
Originally posted January
23, 2025