Skip to main content

Compute D matrix for plane problems

Short answer

For plane stress, D = [[219780,65934,0],[65934,219780,0],[0,0,76923]] MPa. For plane strain, D = [[269231,115385,0],[115385,269231,0],[0,0,76923]] MPa. With epsilon_x=0.001, plane stress gives sigma_x=219.8, sigma_y=65.9, sigma_z=0 MPa; plane strain gives sigma_x=269.2, sigma_y=115.4, sigma_z=115.4 MPa.

For a steel with E=200E = 200 GPa and ν=0.3\nu = 0.3, compute the elasticity matrix [D][D] for (a) plane stress and (b) plane strain. For a strain state εx=0.001\varepsilon_x = 0.001, εy=0\varepsilon_y = 0, γxy=0\gamma_{xy} = 0, find the stresses σx\sigma_x, σy\sigma_y and σz\sigma_z in each case.

(a) Plane stress

E1−ν2=200×1031−0.09=219 780 MPa\frac{E}{1 - \nu^2} = \frac{200\times10^3}{1 - 0.09} = 219\,780\ \text{MPa} [D]=219 780[10.300.310000.35]=[219 78065 934065 934219 78000076 923] MPa[D] = 219\,780\begin{bmatrix} 1 & 0.3 & 0 \\ 0.3 & 1 & 0 \\ 0 & 0 & 0.35 \end{bmatrix} = \begin{bmatrix} 219\,780 & 65\,934 & 0 \\ 65\,934 & 219\,780 & 0 \\ 0 & 0 & 76\,923 \end{bmatrix}\ \text{MPa}

Note that 1−ν2=0.35\dfrac{1 - \nu}{2} = 0.35 and 219 780(0.35)=76 923=G=E2(1+ν)219\,780(0.35) = 76\,923 = G = \dfrac{E}{2(1 + \nu)}.

(b) Plane strain

E(1+ν)(1−2ν)=200×1031.3(0.4)=384 615 MPa\frac{E}{(1 + \nu)(1 - 2\nu)} = \frac{200\times10^3}{1.3(0.4)} = 384\,615\ \text{MPa} [D]=384 615[0.70.300.30.70000.2]=[269 231115 3850115 385269 23100076 923] MPa[D] = 384\,615\begin{bmatrix} 0.7 & 0.3 & 0 \\ 0.3 & 0.7 & 0 \\ 0 & 0 & 0.2 \end{bmatrix} = \begin{bmatrix} 269\,231 & 115\,385 & 0 \\ 115\,385 & 269\,231 & 0 \\ 0 & 0 & 76\,923 \end{bmatrix}\ \text{MPa}

Stresses for εx=0.001\varepsilon_x = 0.001

Caseσx\sigma_x (MPa)σy\sigma_y (MPa)σz\sigma_z (MPa)
Plane stress219 780(0.001)=219.8219\,780(0.001) = 219.865 934(0.001)=65.965\,934(0.001) = 65.90
Plane strain269 231(0.001)=269.2269\,231(0.001) = 269.2115 385(0.001)=115.4115\,385(0.001) = 115.4ν(σx+σy)=0.3(384.6)=115.4\nu(\sigma_x + \sigma_y) = 0.3(384.6) = 115.4

The shear modulus term is the same in both cases. Plane strain gives higher stresses because εz=0\varepsilon_z = 0 is enforced.

Answer: plane stress σx=219.8\sigma_x = 219.8 MPa, σy=65.9\sigma_y = 65.9 MPa, σz=0\sigma_z = 0; plane strain σx=269.2\sigma_x = 269.2 MPa, σy=115.4\sigma_y = 115.4 MPa, σz=115.4\sigma_z = 115.4 MPa.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Finite Element Method chapter 6 questions →