Ultimate and Yield Plastic Bending Allowables Using Cozzone Stresses

The following analysis calculates the Cozzone stresses and plastic bending allowables given material modulus, ultimate strength, yield, ultimate strain, and section shape factor.

  • Use of Cozzone method for plastic bending in short transverse (ST) direction is not recommended

  • Elevated temperature values of fm and fo must be based on a corrected full range stress-strain curve at temperature

  • Use smaller of tension or compression stress-strain curves and assume other is identical - small conservative error

Material and Section Properties Input


Unit System: /

Material Allowables





Bending Section Properties



×

Cozzone Method Equation Derivation

Background

The Cozzone method provides a practical approach for calculating plastic bending allowables based on the material's stress-strain curve. This method accounts for the redistribution of stress across a section as the outer fibers yield and enter the plastic range.

1. Elongation Strain (εu')

The elongation strain represents the plastic component of the total ultimate strain:

εu' = εu - Ftu/E

Where:

  • εu = total ultimate strain (elastic + plastic)
  • Ftu/E = elastic component of strain
  • εu' = plastic component (elongation)

2. Material Shape Factor in Plastic Range (m)

The parameter m characterizes the shape of the stress-strain curve in the plastic region using the Ramberg-Osgood formulation:

m = log(εu' / 0.002) / log(Ftu / Fty)

This relates the plastic strain at ultimate (εu') to the 0.2% offset yield criterion and the ratio of ultimate to yield strength.

3. Ultimate Intercept Stress (fou)

The intercept stress is derived by integrating the stress distribution across a beam section in bending, assuming:

  • Plane sections remain plane (Bernoulli-Euler assumption)
  • Material follows Ramberg-Osgood stress-strain relationship
  • Maximum fiber stress equals Ftu

Derivation Steps:

  1. The strain distribution across the section is linear: ε = (y/c) × εmax
  2. The stress at any fiber is obtained from the Ramberg-Osgood equation
  3. The moment is found by integrating: M = ∫ σ × y × dA
  4. For a rectangular section, this leads to the bending stress equation: Fb = fm + fo × (k - 1)

The intercept stress formula accounts for the nonlinear stress distribution:

fou = ((6/εu2) × [(1/3) × (Ftu/E)2 + εu' × ((m+1)/(m+2)) × (Ftu/E) + (m/(2m+1)) × (εu')2] - 2) × Ftu

The three terms in brackets represent:

  • (1/3) × (Ftu/E)2 - Contribution from elastic region
  • εu' × ((m+1)/(m+2)) × (Ftu/E) - Transition region contribution
  • (m/(2m+1)) × (εu')2 - Fully plastic region contribution

4. Yield Intercept Stress (foy)

Similarly, the yield intercept stress uses the yield properties:

εy = 0.002 + Fty/E

foy = ((6/εy2) × [(1/3) × (Fty/E)2 + (εy - Fty/E) × ((m+1)/(m+2)) × (Fty/E) + (m/(2m+1)) × (εy - Fty/E)2] - 2) × Fty

5. Section Shape Factor (k)

The shape factor k accounts for different cross-sectional geometries:

  • k = 1.000 - Two-noded beam (pure bending, linear stress distribution)
  • k = 1.273 - Thin round tube
  • k = 1.333 - Hour-glass section
  • k = 1.500 - Rectangular section (fully plastic)
  • k = 1.698 - Round section (solid circular)
  • k = 2.000 - Diamond section

6. Plastic Bending Allowables

The final bending allowables combine the maximum fiber stress with the intercept stress contribution:

Fbu = fmu + fou × (k - 1) (ultimate)

Fby = fmy + foy × (k - 1) (yield)

Where:

  • fmu = Ftu (maximum fiber stress at ultimate)
  • fmy = Fty (maximum fiber stress at yield)
  • (k - 1) factor adjusts for section geometry

Physical Interpretation

The term (k - 1) represents the benefit gained from stress redistribution as the outer fibers yield. A higher k value indicates greater capacity for plastic stress redistribution, resulting in higher bending strength compared to elastic theory.

×

Ramberg-Osgood Shape Factors: m vs n

Overview

The Ramberg-Osgood equation describes nonlinear stress-strain behavior in the plastic range. There are two common parameterizations using different shape factors:

Original Form (using n)

The classic Ramberg-Osgood equation is:

ε/ε₀ = σ/σ₀ + (σ/σ₀)n

Where:

  • n = strain hardening exponent (material constant)
  • σ₀ = reference stress (often Fty, the 0.2% offset yield stress)
  • ε₀ = reference strain (often 0.002 for 0.2% offset)
  • n typically ranges from ~5 to 30+ depending on material

Physical Meaning of n:

  • Higher n → sharper yield transition, material behaves more like elastic-perfectly-plastic
  • Lower n → gradual yielding, more rounded stress-strain curve
  • n → ∞ → approaches perfectly elastic-plastic (no strain hardening)

Cozzone Form (using m)

The Cozzone method uses a related parameter m that characterizes the plastic region specifically:

m = log(ε′u / 0.002) / log(Ftu / Fty)

Where:

  • ε′u = plastic strain at ultimate (total strain minus elastic strain)
  • 0.002 = reference plastic strain (0.2%)
  • Ftu/Fty = ratio of ultimate to yield strength
  • m is derived directly from material test data

Relationship Between m and n

The parameters are related but not identical:

  • n is a general material constant valid across the entire plastic range
  • m is specifically derived for the region between yield and ultimate stress
  • Both describe "roundness" of the stress-strain curve in plastic range
  • Approximate relationship: m ≈ 1/(n-1) for many materials
  • They're inversely related - high n corresponds to low m

Typical Values

Aluminum Alloys:

  • n ≈ 15-25
  • m ≈ 20-25

Steel:

  • n ≈ 5-10 (sharper yield transition)
  • m ≈ higher values

Copper:

  • n ≈ 20-30
  • m ≈ varies widely

Use in Cozzone Method

The parameter m appears in the intercept stress calculations with terms like:

  • (m+1)/(m+2) - weight factor for transition region
  • m/(2m+1) - weight factor for fully plastic region

These factors come from integrating the stress distribution across the bending section using the Ramberg-Osgood stress-strain relationship.

Key Advantage

The Cozzone method is convenient because m can be directly calculated from standard material properties (Ftu, Fty, εu, E) without needing to fit the entire stress-strain curve!




References

Cozzone, F. P., Melcon, M. A., & Hoblit, F. M. (1943). Bending Strength in the Plastic Range. Journal of the Aeronautical Sciences, Vol. 10, No. 5, pp. 137-151.

Bruhn, E. F. (1973). Analysis and Design of Flight Vehicle Structures. Tri-State Offset Company, Cincinnati, OH. Section C3.3.

Ramberg, W., & Osgood, W. R. (1943). Description of Stress-Strain Curves by Three Parameters. NACA Technical Note 902. National Bureau of Standards, Washington, D.C.


DISCLAIMER: The methods and calculation on aNOWysis.com are provided as means to educate the reader. There is no warranty, expessed or implied, as to the relevance or accuracy of the methods or calculations. All calculations should be checked by a qualified analyst for appropriateness and accuracy. We will not be held liable for any damages, injuries or deaths caused by using methods or calculations from this site.