 Solution

In the Voigt elastic approximation, every local strain is $E$. The corresponding local stress is ${\sigma }_V(x)=D_{\epsilon }W(E,x)$. Differentiating the averaged strain energy density under the integral givesThe uniform strain need not generate an equilibrated local stress; this calculation concerns the stated approximation, not an assertion that it solves the heterogeneous boundary problem.

For the Reuss elastic approximation, use the convex conjugateUnder differentiable strictly convex duality, the local strain produced by uniform stress $S$ is $e_R(x)=D_S{W}^{\ast }(S,x)$. Set $F(S)=⟨W^{\ast }(S,x)⟩$. Averaging the local strains gives $E=DF(S)$. At the maximizing stress in $W_R(E)=F^{\ast }(E)$, the stationarity condition is precisely $DF(S)=E$. ConsequentlyThus the derivative of each approximate energy is its own mean stress, includingFor positive linear elasticity, $W=e:C(x)e/2$ and $W^{\ast }=S:C(x{)}^{-1}S/2$, giving the useful check $C_V=⟨C⟩$ and $C_R=⟨C^{-1}{⟩}^{-1}$. More general convex duality gives the same conjugacy relation with subgradients when differentiability or strict convexity is absent.