Skip to main content

Plane stress, plane strain and axisymmetric

The basic relationship between stress and strain is \sigma = C \epsilon$ and the stiffness matrix CC is of the form

C=[c11c12c13000c21c22c23000c31c32c33000000c44000000c55000000c66]C = \begin{bmatrix} c_{11} & c_{12} & c_{13} & 0 & 0 & 0 \\ c_{21} & c_{22} & c_{23} & 0 & 0 & 0 \\ c_{31} & c_{32} & c_{33} & 0 & 0 & 0 \\ 0 & 0 & 0 & c_{44} & 0 & 0 \\ 0 & 0 & 0 & 0 & c_{55} & 0 \\ 0 & 0 & 0 & 0 & 0 & c_{66} \\ \end{bmatrix}

For plane strain the strains ϵzz\epsilon_{zz}, ϵyz\epsilon_{yz} and ϵzz\epsilon_{zz} are zero, and we are not directly interested in the corresponding stresses so we can reduce the stiffness matrix to a 3 ×\times 3 matrix.

C=[c11c120c21c22000c44]C = \begin{bmatrix} c_{11} & c_{12} & 0 \\ c_{21} & c_{22} & 0 \\ 0 & 0 & c_{44} \\ \end{bmatrix}

Axisymmetric is similar to plane strain, reducing the problem to two dimensions, radial (xx) and axial (yy), but in this instance the strain in the third (hoop zz) direction is related to the radial strain. In this case the strains ϵyz\epsilon_{yz} and ϵzx\epsilon_{zx} are zero, and we are not directly interested in the corresponding stresses so we can reduce the stiffness matrix to a 4 ×\times 4 matrix.

C=[c11c12c130c21c22c230c31c32c330000c44]C = \begin{bmatrix} c_{11} & c_{12} & c_{13} & 0 \\ c_{21} & c_{22} & c_{23} & 0 \\ c_{31} & c_{32} & c_{33} & 0 \\ 0 & 0 & 0 & c_{44} \\ \end{bmatrix}

For plane stress the stresses σzz\sigma_{zz}, σyz\sigma_{yz} and σzx\sigma_{zx} are zero, so we can partition the stiffness matrix

C=[c11c120c1300c21c220c230000c44000c31c320c33000000c55000000c66]C = \begin{bmatrix} c_{11} & c_{12} & 0 & c_{13} & 0 & 0 \\ c_{21} & c_{22} & 0 & c_{23} & 0 & 0 \\ 0 & 0 & c_{44} & 0 & 0 & 0 \\ c_{31} & c_{32} & 0 & c_{33} & 0 & 0 \\ 0 & 0 & 0 & 0 & c_{55} & 0 \\ 0 & 0 & 0 & 0 & 0 & c_{66}\\ \end{bmatrix}

so that the stress strain relationship can focus on the terms of interest ii and the terms that can be removed rr.

[σi0]=[CiiCirCriCrr][ϵiϵr]\begin{bmatrix} \sigma_i \\ 0 \\ \end{bmatrix} = \begin{bmatrix} C_{ii} & C_{ir} \\ C_{ri} & C_{rr} \\ \end{bmatrix} \begin{bmatrix} \epsilon_i \\ \epsilon_r \\ \end{bmatrix}

This allows the strain corresponding to unstressed term jj to be removed giving the updated equation

σi=(CiiCirCrr1Cri)ϵi\sigma_i = (C_{ii} - C_{ir} C_{rr}^{-1} C_{ri}) \epsilon_i

where CiiCirCrr1CriC_{ii} - C_{ir} C_{rr}^{-1} C_{ri} is the compliance matrix for plane stress.