Formula sheet

Polymer Composite and Ceramic Materials Formula Sheet

Materials engineering formula sheet for polymers, composites, and ceramics covering mixture rules, anisotropy, creep, moisture, voids, thermal stress, fracture, proof screening, and validation.

This formula sheet collects first-pass calculations for polymer, composite, and ceramic material systems. Use it for hand checks, design reviews, test planning, and interpretation of material data. It is not a substitute for qualified allowables, material standards, supplier data, coupon testing, environmental conditioning, or project-specific design rules.

Polymers, composites, and ceramics are strongly controlled by process history, defects, interfaces, temperature, moisture, loading rate, and inspection capability. A formula is useful only when the material system behind the formula matches the real part.

Conventions and Units

State units before comparing values:

  • modulus, stress, and strength: Pa, MPa, or GPa;
  • density: kg/m^3 or g/cm^3;
  • thickness, flaw size, and bond length: m or mm;
  • temperature change: K or deg C difference;
  • moisture content and void content: percent or fraction, stated clearly;
  • creep time: s, h, or years with temperature and stress level;
  • probability quantities: dimensionless, with sample volume or area basis.

Do not mix mass fraction and volume fraction. Composite stiffness formulas usually use volume fraction. Purchasing and batch records often use mass fraction.

Volume Fraction from Mass and Density

For constituents with masses m_i and densities \rho_i, volume fraction is:

\displaystyle V_i=\frac{m_i/\rho_i}{\sum_j m_j/\rho_j}

For a fiber and matrix composite:

\displaystyle V_f=\frac{m_f/\rho_f}{m_f/\rho_f+m_m/\rho_m}

and:

V_m=1-V_f

where V_f is fiber volume fraction and V_m is matrix volume fraction.

Worked Check

A coupon contains m_f=58\ \text{g} carbon fiber and m_m=42\ \text{g} epoxy. Use:

\rho_f=1.78\ \text{g/cm}^3
\rho_m=1.20\ \text{g/cm}^3

Then:

\displaystyle V_f=\frac{58/1.78}{58/1.78+42/1.20}
\displaystyle V_f=\frac{32.58}{32.58+35.00}=0.482

The fiber volume fraction is about 48.2\%.

Engineering Comment

This result is lower than the mass fraction because carbon fiber is denser than the matrix. If a stiffness calculation accidentally uses mass fraction as volume fraction, the predicted modulus can be badly biased.

Composite Density

For a void-free composite:

\rho_c=V_f\rho_f+V_m\rho_m

When mass fractions w_i are known:

\displaystyle \rho_c=\frac{1}{\sum_i w_i/\rho_i}

These equations assume no voids and no significant density change during cure, sintering, crystallization, or chemical reaction.

Longitudinal Rule of Mixtures

For a unidirectional composite loaded along the fiber direction:

E_1\approx V_fE_f+V_mE_m

where:

  • E_1 is longitudinal modulus;
  • E_f is fiber modulus;
  • E_m is matrix modulus.

Strength is sometimes screened similarly:

\sigma_{1}\approx V_f\sigma_f+V_m\sigma_m

but this is much less reliable than the modulus relation because failure depends on fiber defects, matrix strain limit, interface quality, compression behavior, waviness, stress concentration, and test method.

Worked Check

Use:

V_f=0.55
E_f=230\ \text{GPa}
E_m=3.0\ \text{GPa}

Then:

E_1=0.55(230)+0.45(3.0)=127.9\ \text{GPa}

The composite is stiff in the fiber direction, but this value says nothing about transverse modulus, interlaminar tension, bearing, impact damage, or delamination.

Transverse Modulus Screen

A simple inverse rule for transverse modulus is:

\displaystyle \frac{1}{E_2}\approx\frac{V_f}{E_f}+\frac{V_m}{E_m}

so:

\displaystyle E_2\approx\left(\frac{V_f}{E_f}+\frac{V_m}{E_m}\right)^{-1}

This is a screening model. More accurate composite micromechanics depends on fiber packing, matrix constraint, interface quality, voids, temperature, moisture, and test direction.

Worked Check

Using the same material system:

\displaystyle \frac{1}{E_2}=\frac{0.55}{230}+\frac{0.45}{3.0}=0.1524\ \text{GPa}^{-1}

Therefore:

E_2=6.56\ \text{GPa}

The transverse modulus is about 5\% of the longitudinal modulus. This is why ply direction and load path matter.

Shear Modulus and Poisson Ratio

For an isotropic phase:

\displaystyle G=\frac{E}{2(1+\nu)}

where G is shear modulus and \nu is Poisson ratio.

For a composite lamina, do not assume isotropy. A unidirectional ply needs direction-specific properties such as E_1, E_2, G_{12} and \nu_{12}. If only isotropic formulas are used, bearing, torsion, shear transfer, and free-edge behavior can be misread.

Laminate Modulus Screen

For a rough axial modulus screen of a laminate loaded in one direction:

\displaystyle E_{x,lam}\approx\sum_k E_{x,k}\frac{t_k}{t}

where:

  • E_{x,k} is the effective modulus contribution of ply k in the load direction;
  • t_k is ply thickness;
  • t=\sum_k t_k is total laminate thickness.

This formula is only a screening estimate. Real laminate analysis transforms orthotropic ply stiffness, assembles laminate stiffness matrices, and checks coupling, symmetry, balance, damage tolerance, open holes, compression after impact, and free-edge effects.

Specific Stiffness and Specific Strength

For weight-sensitive designs:

\displaystyle \text{specific stiffness}=\frac{E}{\rho}
\displaystyle \text{specific strength}=\frac{\sigma_{allow}}{\rho}

These ratios are useful for comparison, but they do not include joining, inspection, fire, damage tolerance, environmental knockdowns, cost, process yield, or repair.

Average Adhesive Shear Stress

For a simple bonded overlap carrying load P:

\displaystyle \tau_{avg}=\frac{P}{bL_b}

where b is bond width and L_b is effective bond length.

Worked Check

If:

P=8.0\ \text{kN}
b=25\ \text{mm}
L_b=80\ \text{mm}

then:

\displaystyle \tau_{avg}=\frac{8000}{0.025(0.080)}=4.0\times10^6\ \text{Pa}=4.0\ \text{MPa}

Engineering Comment

Average shear is not peak shear. Real joints have peel stress, edge stress, adherend bending, surface preparation sensitivity, cure history, temperature, moisture, fatigue, and inspection limits. A low average value can still fail at the edge.

Void Content from Density

If theoretical void-free density is \rho_{theory} and measured composite density is \rho_{meas}:

\displaystyle V_v=1-\frac{\rho_{meas}}{\rho_{theory}}

where V_v is void volume fraction.

Worked Check

For:

\rho_{theory}=1.54\ \text{g/cm}^3
\rho_{meas}=1.50\ \text{g/cm}^3

the void fraction is:

\displaystyle V_v=1-\frac{1.50}{1.54}=0.0260

So:

V_v=2.6\%

Engineering Comment

Void content should be interpreted with process route and inspection evidence. A small void fraction may be acceptable in one laminate and unacceptable near a highly loaded bond, pressure boundary, fatigue detail, or compression-critical zone.

Moisture Uptake

Moisture content by mass is:

\displaystyle M_t=\frac{m_t-m_d}{m_d}\times100\%

where m_t is conditioned mass and m_d is dry mass.

For early-time one-dimensional diffusion in a thin plate, a screening relation is:

\displaystyle \frac{M_t}{M_{\infty}}\approx\frac{4}{h}\sqrt{\frac{Dt}{\pi}}

where:

  • M_{\infty} is saturation moisture content;
  • h is plate thickness;
  • D is diffusion coefficient;
  • t is exposure time.

This approximation is valid only at early uptake and for the assumed geometry and boundary condition.

Polymer Creep Compliance

For linear viscoelastic screening under constant stress:

\epsilon(t)=J(t)\sigma

where:

  • \epsilon(t) is strain at time t;
  • J(t) is creep compliance;
  • \sigma is applied stress.

If an initial elastic strain and creep strain are separated:

\displaystyle \epsilon(t)=\frac{\sigma}{E_0}+\epsilon_{creep}(t)

Worked Check

A polymer clip carries:

\sigma=12\ \text{MPa}

At service temperature, the creep compliance at one year is:

J(1\ \text{year})=0.090\ \text{GPa}^{-1}

Convert stress to GPa:

12\ \text{MPa}=0.012\ \text{GPa}

Then:

\epsilon=0.090(0.012)=0.00108

So:

\epsilon=0.108\%

Engineering Comment

The value is meaningful only for the stated temperature, stress, material condition, and time. Polymers may show nonlinear creep near high stress, elevated temperature, chemical exposure, or long service time.

Creep Retention Factor

A simple retention factor can compare long-term property to short-term property:

\displaystyle R_E(t)=\frac{E(t)}{E_0}

or for allowable stress:

\displaystyle R_{\sigma}(t)=\frac{\sigma_{allow}(t)}{\sigma_{allow,0}}

Use the factor only with the same test temperature, environment, strain rate, and failure criterion.

Thermal Strain and Thermal Stress

Free thermal strain is:

\epsilon_{th}=\alpha\Delta T

If expansion is fully restrained in a simple linear elastic member:

\sigma_{th}=E\alpha\Delta T

For a plane-strain-like restraint, a conservative screen may use:

\displaystyle \sigma_{th}\approx\frac{E\alpha\Delta T}{1-\nu}

where \alpha is coefficient of thermal expansion.

Worked Check

A ceramic layer with:

E=210\ \text{GPa}
\alpha=8.0\times10^{-6}\ \text{K}^{-1}

is fully restrained during:

\Delta T=120\ \text{K}

Simple restrained stress:

\sigma_{th}=210\times10^9(8.0\times10^{-6})(120)
\sigma_{th}=202\ \text{MPa}

Engineering Comment

This stress is a warning, not a full thermal-stress solution. Coating thickness, substrate stiffness, temperature gradient, creep, cracking, residual stress, edge geometry, and interface toughness can dominate the real failure mode.

Fracture Toughness and Critical Flaw Size

For a brittle material or ceramic under tensile stress:

K_I=Y\sigma\sqrt{\pi a}

Failure is screened when:

K_I\geq K_{IC}

Critical flaw size is:

\displaystyle a_c=\frac{1}{\pi}\left(\frac{K_{IC}}{Y\sigma}\right)^2

where:

  • K_I is stress intensity factor;
  • K_{IC} is fracture toughness;
  • Y is geometry factor;
  • a is flaw size;
  • \sigma is tensile stress.

Worked Check

A ceramic has:

K_{IC}=4.0\ \text{MPa}\sqrt{\text{m}}

with:

Y=1.1

and tensile stress:

\sigma=120\ \text{MPa}

Then:

\displaystyle a_c=\frac{1}{\pi}\left(\frac{4.0}{1.1(120)}\right)^2
a_c=2.92\times10^{-4}\ \text{m}=0.292\ \text{mm}

Engineering Comment

If inspection cannot reliably detect flaws below this size in the critical region, the design needs lower stress, tougher material, proof screening, compressive residual stress, geometry change, or a different inspection basis.

Weibull Failure Probability

For brittle strength scatter, a common Weibull screening form is:

\displaystyle P_f=1-\exp\left[-\left(\frac{\sigma}{\sigma_0}\right)^m\frac{V}{V_0}\right]

where:

  • P_f is probability of failure;
  • \sigma is applied tensile stress;
  • \sigma_0 is scale stress for reference volume V_0;
  • m is Weibull modulus;
  • V is stressed volume.

This relation is sensitive to stress distribution, flaw population, volume basis, surface condition, proof testing, and statistical evidence.

Proof Stress Screening

If a brittle component survives proof stress \sigma_p, and service tensile stress is \sigma_s, a simple proof ratio is:

\displaystyle R_p=\frac{\sigma_p}{\sigma_s}

This is not a complete reliability proof. Proof testing can damage some materials, miss time-dependent degradation, or fail to represent multiaxial stress and contact damage.

Permeability Flux

For a thin polymer barrier with steady one-dimensional transport:

\displaystyle J=\frac{P\Delta p}{h}

where:

  • J is molar or mass flux, depending on the units of permeability P;
  • \Delta p is pressure or concentration driving difference;
  • h is barrier thickness.

Permeability depends on penetrant, polymer chemistry, crystallinity, temperature, humidity, aging, fillers, defects, and seams.

Rheology and Viscosity

For simple Newtonian shear:

\tau=\mu\dot{\gamma}

where:

Many polymer melts, filled resins, adhesives, and slurries are non-Newtonian. Their apparent viscosity can change with shear rate, temperature, filler loading, cure state, moisture, and residence time.

Damping and Loss Factor

Loss factor can be screened as:

\displaystyle \eta=\frac{E''}{E'}

where E' is storage modulus and E'' is loss modulus.

High damping may help vibration and acoustic response, but the same viscoelastic behavior may increase heat buildup, creep, hysteresis, and temperature sensitivity.

Measurement and Validation Checks

For a material property used in release:

\displaystyle \text{margin}=\frac{\text{allowable}-\text{demand}}{\text{demand}}

or:

\displaystyle MS=\frac{\text{allowable}}{\text{demand}}-1

where MS is margin of safety.

When the property comes from test data, record:

  1. specimen orientation and geometry;
  2. process route, lot, cure, sinter, heat treatment, or conditioning;
  3. test temperature, humidity, strain rate, and environment;
  4. measurement uncertainty and calibration status;
  5. defect population and inspection method;
  6. statistical basis and outlier treatment;
  7. whether the coupon represents the real part.

Common Mistakes

Common errors include:

  • using mass fraction where volume fraction is required;
  • applying longitudinal composite modulus to transverse, shear, bearing, or interlaminar loading;
  • using short-term polymer strength for long-term creep-controlled service;
  • ignoring moisture, temperature, ultraviolet exposure, chemicals, or strain rate;
  • treating void content as harmless without checking location and load path;
  • using ceramic average strength without flaw size, surface finish, proof test, or volume effect;
  • accepting a material change without repeating process, inspection, and validation evidence.

Validity Limits

These formulas are screening tools. They are most useful for early sizing, plausibility checks, test planning, and independent review. They become unsafe when used beyond their assumptions: nonlinear polymer behavior, damaged composites, thick laminates, impact damage, complex joints, multiaxial ceramic stress, thermal gradients, residual stress, manufacturing defects, aging, or unvalidated material substitutions.

For release decisions, connect every formula to material lot, process route, environmental conditioning, inspection evidence, uncertainty, and the failure mode that actually controls the design.

REF

See also