Laminate Calculator — Calculation Documentation

Oxess / C-Sink Sports  ·  Composite laminate analysis for snowboards and skateboards

Contents

  1. Coordinate System & Conventions
  2. Material Database
  3. Cross-Section Analysis
    1. Neutral Axis
    2. Bending Stiffness EI
    3. Torsional Stiffness GJ
    4. Stiffness Balance
    5. Strength Balance
    6. Strain Balance
  4. Board-Level Analysis
  5. Standard Hardness (EI equiv)
  6. Standard Torsion Stiffness (GJ equiv)
  7. Balance Uniformity
  8. 3-Point Flex Test
  9. All Displayed Metrics — Quick Reference

1 · Coordinate System & Conventions

The cross-section coordinate $y$ is measured upward from the base (ski/board riding surface). The wood core sits between the below-core (tension) and above-core (compression) composite layers.

All thickness quantities are in mm. Moduli are in GPa. Stiffness quantities (EI, GJ) are computed per unit width (i.e. for a strip of width $b = 1\,\text{mm}$), so the unit is $\text{GPa·mm}^3/\text{mm}$. To obtain the stiffness of the full board, multiply by the board width in mm.

2 · Material Database

Each composite ply is characterised by the following properties:

PropertySymbolUnitDescription
Tensile modulus$E_t$GPaYoung's modulus on the tension side of NA
Compressive modulus$E_c$GPaYoung's modulus on the compression side of NA
Tensile strength$\sigma_t$MPaFailure stress in tension
Compressive strength$\sigma_c$MPaFailure stress in compression
In-plane shear modulus$G_{xy}$GPaDominates torsion; highest for ±45° plies
Density$\rho$g/cm³Layer density
Default thickness$t$mmNominal laminate thickness

Wood core types (paulownia, poplar, ash, beech) use the same property set with separate tension and compression moduli representing the anisotropy of wood along its grain direction.

±45° layers have low $E_t$ (≈ 14–18 GPa) but very high $G_{xy}$ (10–40 GPa), so they contribute little to bending stiffness but dominate torsional stiffness. 0/90° and UD layers are the opposite — high $E$, low $G_{xy}$.

3 · Cross-Section Analysis

Function analyzeSection(belowLayers, coreThickness, aboveLayers, coreType) analyses a single cross-section at a given position along the board. Layers are stacked from $y=0$ upward in the order: below-core plies → wood core → above-core plies.

3.1 · Neutral Axis

The neutral axis $y_\text{NA}$ is the $y$-position where bending stress is zero. Because the tensile and compressive moduli of composite plies may differ ($E_t \neq E_c$), the neutral axis is found iteratively:

$$y_\text{NA} = \frac{\sum_k E_k \cdot A_k \cdot y_{m,k}}{\sum_k E_k \cdot A_k}$$

where for each layer $k$: $A_k = t_k \cdot b$ is the cross-sectional area (per unit width $b=1$), $y_{m,k}$ is the centroidal $y$-coordinate, and $E_k$ is assigned as $E_t$ for layers fully below $y_\text{NA}$ and $E_c$ for layers fully above it. Layers split by the neutral axis contribute both moduli proportionally. Iteration continues until $|y_\text{NA}^\text{new} - y_\text{NA}^\text{old}| < 10^{-8}\,\text{mm}$ (typically converges in < 10 steps).

3.2 · Bending Stiffness EI

The flexural rigidity per unit width is computed by summing the parallel-axis contributions of each layer about $y_\text{NA}$:

$$EI = \sum_k E_k \left[ \frac{b\,t_k^3}{12} + b\,t_k \left(\bar{y}_k - y_\text{NA}\right)^2 \right]$$

where $\bar{y}_k = (y_{\text{bot},k} + y_{\text{top},k})/2$ is the layer centroid. For layers split by $y_\text{NA}$, the tension and compression sub-portions are treated separately with their respective moduli.

Unit: $\text{GPa·mm}^3/\text{mm}$ (per unit width). Displayed scaled by $10^{-3}$ as $\times 10^3\,\text{GPa·mm}^3/\text{mm}$.

3.3 · Torsional Stiffness GJ

Torsional stiffness is computed using Classical Laminate Theory (CLT). The plate torsion stiffness $D_{66}$ per unit width is:

$$D_{66} = \frac{1}{3}\sum_k G_{xy,k}\left(z_{\text{top},k}^3 - z_{\text{bot},k}^3\right)$$

where $z = y - y_\text{geom}$ is measured from the geometric mid-plane $y_\text{geom} = t_\text{total}/2$ (not from the neutral axis). The factor-of-2 relationship between $D_{66}$ and the effective torsional stiffness for a plate gives:

$$GJ = 2\,D_{66} = \frac{2}{3}\sum_k G_{xy,k}\left(z_{\text{top},k}^3 - z_{\text{bot},k}^3\right)$$

Unit: $\text{GPa·mm}^3/\text{mm}$ (per unit width). To get the absolute torsional stiffness of the board, multiply by the board width: $GJ_\text{board} = GJ \times w$.

The $z^3$ dependence means layers far from the geometric mid-plane contribute much more to torsion. ±45° layers with high $G_{xy}$ placed at the board surfaces (furthest from mid-plane) are most effective for torsional rigidity.

3.4 · Stiffness Balance

Stiffness balance measures the symmetry of bending stiffness between the tension and compression faces:

$$\text{Balance} = \frac{\sum_{k \in \text{tension}} E_{t,k}\,t_k}{\sum_{k \in \text{compression}} E_{c,k}\,t_k} = \frac{(E \times t)_\text{tension}}{(E \times t)_\text{compression}}$$

where the sums run over the portions of each layer on the respective side of $y_\text{NA}$. A value of 1.0 means both faces are equally stiff.

3.5 · Strength Balance

Analogous to stiffness balance, using strength instead of modulus:

$$\text{Strength Balance} = \frac{\sum_{k \in \text{tension}} \sigma_{t,k}\,t_k}{\sum_{k \in \text{compression}} \sigma_{c,k}\,t_k}$$

A value of 1.0 means both faces accumulate the same "strength-thickness product". Values far from 1.0 indicate that one face will yield before the other under increasing bending load.

3.6 · Strain Balance

Identifies which side of the laminate has the layer that will fail first (the weakest link). The failure strain of each layer is:

$$\varepsilon_{\text{fail},k} = \frac{\sigma_k}{E_k}$$

The minimum failure strain on each side is found:

$$\varepsilon_\text{tension}^\text{min} = \min_{k \in \text{tension}} \frac{\sigma_{t,k}}{E_{t,k}}, \qquad \varepsilon_\text{compression}^\text{min} = \min_{k \in \text{compression}} \frac{\sigma_{c,k}}{E_{c,k}}$$
$$\text{Strain Balance} = \frac{\varepsilon_\text{tension}^\text{min}}{\varepsilon_\text{compression}^\text{min}}$$

A value of 1.0 means both sides reach their failure strain simultaneously. Values significantly different from 1.0 indicate the lower-side face will fail first.

4 · Board-Level Analysis

Function analyzeBoard(belowLayers, aboveLayers, coreProfile, coreType) applies analyzeSection at every point in the core thickness profile. The core profile is a list of $(x_i, t_{\text{core},i})$ pairs giving the core thickness at each position $x_i$ along the board length.

The result is a 1-D array of section analyses indexed by position. Charts show how EI, GJ, balance, and neutral axis vary from tip to tip. The midpoint section (at the position of maximum core thickness) is used for the headline metric values.

5 · Standard Hardness (EI equiv)

The single-number "Standard Hardness" answers: what constant EI would give the same mid-span deflection as the real, varying-stiffness board under a centred 3-point load?

For a simply-supported beam of length $L$ with a central point load $P$, the virtual-work deflection integral is:

$$\frac{\delta}{P} = \int_0^L \frac{m(x)^2}{EI(x)}\,dx$$

where $m(x)$ is the bending moment diagram for a unit load ($P=1$):

$$m(x) = \begin{cases} x/2 & 0 \le x \le L/2 \\ (L-x)/2 & L/2 \le x \le L \end{cases}$$

The equivalent uniform stiffness that produces the same deflection is then:

$$EI_\text{equiv} = \frac{L^3}{48 \displaystyle\int_0^L \frac{m(x)^2}{EI(x)}\,dx}$$

The integral is evaluated numerically using the trapezoidal rule over the core profile points. The $m(x)^2/EI$ weighting means the stiff centre contributes more than the thin, flexible tips — soft tips pull the equivalent stiffness well below the midpoint value.

Compare Standard Hardness to EI midpoint: the ratio reveals how much the tip taper softens the overall board response. For a board with no taper, both values are equal.

6 · Standard Torsion Stiffness (GJ equiv)

The analogous question for torsion: what constant GJ would give the same total twist angle under an applied end torque as the real board?

For a shaft of varying torsional stiffness $GJ(x)$ subjected to a torque $T$, torsional compliance adds in series. The total twist angle is:

$$\varphi_\text{total} = T \int_0^L \frac{1}{GJ(x)}\,dx$$

The equivalent uniform stiffness that gives the same total twist is:

$$GJ_\text{equiv} = \frac{L}{\displaystyle\int_0^L \frac{1}{GJ(x)}\,dx}$$

This is the harmonic-mean-like average of $GJ$ along the board. Because compliance sums in series, the softest sections dominate — the thin, low-$GJ$ tips have a disproportionately large influence. This is the opposite weighting from the bending case, where the stiff centre dominates.

Compare Standard Torsion Stiffness to GJ midpoint: a large difference means the tips are torsionally much softer than the centre and significantly limit the whole-board torsional rigidity.

Both Standard Hardness and Standard Torsion Stiffness are per unit width (GPa·mm³/mm). Multiply by the board width in mm to obtain absolute board stiffness values for comparison with test data.

7 · Balance Uniformity

Uniformity measures how consistently the Stiffness Balance is maintained along the entire board length:

$$\text{Uniformity} = \left| \text{Balance}_\text{centre} - \text{Balance}_\text{tips} \right|$$

where $\text{Balance}_\text{centre}$ is the balance value at the 50th percentile position along the board, and $\text{Balance}_\text{tips}$ is the average of the first and last three profile points.

At the thin tips, the wood core contributes very little stiffness, so the composite layers on both faces dominate the balance ratio. Any asymmetry between below-core and above-core layups that is masked by the thick core in the centre becomes fully exposed at the tips.

8 · 3-Point Flex Test Simulation

The flex test panel simulates a symmetric 3-point bend with supports placed symmetrically about the board centre. Given the input parameters (applied moment $M$, board width $w$, and distance between supports $s$):

Support positions

$$x_L = \frac{L}{2} - \frac{s}{2}, \qquad x_R = \frac{L}{2} + \frac{s}{2}$$

where $L$ is the board length from the core profile.

Equivalent EI over the test span

The same virtual-work integral as for Standard Hardness is applied, but restricted to the span $[x_L, x_R]$:

$$EI_\text{span} = \frac{s^3}{48 \displaystyle\int_{x_L}^{x_R} \frac{m_s(x)^2}{EI(x)}\,dx}$$

where $m_s(x)$ is the moment diagram for a unit load at the centre of the span.

Predicted deflection

The moment $M$ is applied at the span centre. The equivalent point force and mid-span deflection are:

$$P = \frac{4M}{s}, \qquad \delta = \frac{M \cdot s^2}{12 \cdot EI_\text{span} \cdot w \cdot 10^3}$$

where $M$ is in N·mm, $s$ in mm, $EI_\text{span}$ in GPa·mm³/mm, and $w$ is the board width in mm. The factor $10^3$ converts GPa to N/mm².

Spring rate

$$k = \frac{P}{\delta} \quad [\text{N/mm}]$$

Back-calculation from measured deflection

If a measured deflection $\delta_\text{meas}$ is entered, the actual EI is back-calculated by inverting the deflection formula:

$$EI_\text{measured} = \frac{M \cdot s^2}{12 \cdot \delta_\text{meas} \cdot w \cdot 10^3}$$

The calibration factor $EI_\text{measured} / EI_\text{span}$ indicates whether the real board is stiffer (> 1) or softer (< 1) than the model predicts.

9 · All Displayed Metrics — Quick Reference

Metric Formula / Definition Unit Target / Good range
Stiffness Balance $(E{\times}t)_\text{tension}\ /\ (E{\times}t)_\text{compression}$ 0.85 – 1.15
Uniformity $|$Balance$_\text{centre}$ − Balance$_\text{tips}|$ < 0.05
Strength Balance $(\sigma{\times}t)_\text{tension}\ /\ (\sigma{\times}t)_\text{compression}$ ≈ 1.0
Strain Balance $\varepsilon^\text{min}_\text{tension}\ /\ \varepsilon^\text{min}_\text{compression}$ ≈ 1.0
Neutral Axis $y_\text{NA}$ from iterative $\sum E_k A_k y_k / \sum E_k A_k$ mm from base
E×t Tension $\sum_{k\,\in\,\text{tension}} E_{t,k}\,t_k$ GPa·mm
E×t Compression $\sum_{k\,\in\,\text{compression}} E_{c,k}\,t_k$ GPa·mm
EI midpoint $\sum_k E_k [b t_k^3/12 + b t_k d_k^2]$ at thickest section ×10³ GPa·mm³/mm
EI mean $\int_0^L EI(x)\,dx\ /\ L$ ×10³ GPa·mm³/mm
Standard Hardness $L^3 \Big/ \Big(48 \int_0^L m(x)^2/EI(x)\,dx\Big)$ ×10³ GPa·mm³/mm Compare boards at equal width
GJ midpoint $2D_{66} = \frac{2}{3}\sum_k G_{xy,k}(z_\text{top}^3 - z_\text{bot}^3)$ at thickest section ×10³ GPa·mm³/mm
Standard Torsion Stiffness $L \Big/ \int_0^L 1/GJ(x)\,dx$ ×10³ GPa·mm³/mm Compare boards at equal width
Flex Test EI span $s^3 \Big/ \Big(48 \int_{x_L}^{x_R} m_s^2/EI\,dx\Big)$ ×10³ GPa·mm³/mm
Predicted deflection δ $M s^2 / (12 \cdot EI_\text{span} \cdot w \cdot 10^3)$ mm
Spring rate k $P / \delta = 4M / (s \cdot \delta)$ N/mm
Note on per-unit-width values: EI and GJ quantities are all computed for a strip of width $b = 1\,\text{mm}$. For an absolute comparison with physical measurements or for computing actual board deflection, multiply by the board width $w$ [mm]. The ratio between two boards of the same nominal width is unaffected by this scaling.

Laminate Calculator v2  ·  Oxess / C-Sink Sports  ·  Documentation generated from source code