Micromechanics of composites/Average stress in a RVE

Average Stress in a RVE

Let the average stress in the RVE be defined as

${\langle {\boldsymbol {\sigma }}\rangle :={\cfrac {1}{V}}~\int _{\Omega }{\boldsymbol {\sigma }}(\mathbf {x} )~{\text{dV}}}$

where $V$ is the volume of $\Omega$.

We would like to find out the relation between the average stress in a RVE and the applied tractions on the boundary of the RVE. To do that, recall the relation (see Appendix)

$\int _{\partial {\Omega }}\mathbf {v} \otimes ({\boldsymbol {S}}^{T}\bullet \mathbf {n} )~{\text{dA}}=\int _{\Omega }[{\boldsymbol {\nabla }}\mathbf {v} \cdot {\boldsymbol {S}}+\mathbf {v} \otimes ({\boldsymbol {\nabla }}\bullet {\boldsymbol {S}}^{T})]~{\text{dV}}~.$

If we choose $\mathbf {v}$ such that ${\boldsymbol {\nabla }}\mathbf {v} ={\boldsymbol {\mathit {1}}}$, we have

$\int _{\partial {\Omega }}\mathbf {v} \otimes ({\boldsymbol {S}}^{T}\bullet \mathbf {n} )~{\text{dA}}=\int _{\Omega }[{\boldsymbol {\mathit {1}}}\cdot {\boldsymbol {S}}+\mathbf {v} \otimes ({\boldsymbol {\nabla }}\bullet {\boldsymbol {S}}^{T})]~{\text{dV}}=\int _{\Omega }{\boldsymbol {S}}~{\text{dV}}+\int _{\Omega }\mathbf {v} \otimes ({\boldsymbol {\nabla }}\bullet {\boldsymbol {S}}^{T})~{\text{dV}}~.$

Therefore,

$\int _{\Omega }{\boldsymbol {S}}~{\text{dV}}=-\int _{\Omega }\mathbf {v} \otimes ({\boldsymbol {\nabla }}\bullet {\boldsymbol {S}}^{T})~{\text{dV}}+\int _{\partial {\Omega }}\mathbf {v} \otimes ({\boldsymbol {S}}^{T}\bullet \mathbf {n} )~{\text{dA}}~.$

If we choose ${\boldsymbol {S}}$ to be the stress tensor ${\boldsymbol {\sigma }}$, and involve the symmetry of the stress tensor, we get

$\int _{\Omega }{\boldsymbol {\sigma }}~{\text{dV}}=-\int _{\Omega }\mathbf {v} \otimes ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})~{\text{dV}}+\int _{\partial {\Omega }}\mathbf {v} \otimes ({\boldsymbol {\sigma }}\bullet \mathbf {n} )~{\text{dA}}~.$

Now, the divergence of the stress is zero (from the conservation of linear momentum). Therefore,

$\int _{\Omega }{\boldsymbol {\sigma }}~{\text{dV}}=\int _{\partial {\Omega }}\mathbf {v} \otimes ({\boldsymbol {\sigma }}\bullet \mathbf {n} )~{\text{dA}}~.$

Using the traction boundary condition, we have

$\int _{\Omega }{\boldsymbol {\sigma }}~{\text{dV}}=\int _{\partial {\Omega }}\mathbf {v} \otimes {\bar {\mathbf {t} }}~{\text{dA}}~.$

Now ${\boldsymbol {\nabla }}\mathbf {v} ={\boldsymbol {\mathit {1}}}$ if $\mathbf {v} =\mathbf {x}$. Therefore, we have

$\int _{\Omega }{\boldsymbol {\sigma }}~{\text{dV}}=\int _{\partial {\Omega }}\mathbf {x} \otimes {\bar {\mathbf {t} }}~{\text{dA}}~.$

Hence the average stress is given by

${\langle {\boldsymbol {\sigma }}\rangle :={\cfrac {1}{V}}~\int _{\partial {\Omega }}\mathbf {x} \otimes {\bar {\mathbf {t} }}~{\text{dA}}~.}$

This implies that the average stress is completely determined by the applied tractions!

Symmetry of the average stress and the effect of rigid body translation

Let us now assume that the applied tractions are self equilibrating. Then the resultant forces and moments due to the applied tractions vanish and we have

$\int _{\partial {\Omega }}{\bar {\mathbf {t} }}~{\text{dA}}=\mathbf {0} \qquad {\text{and}}\qquad \int _{\partial {\Omega }}\mathbf {x} \times {\bar {\mathbf {t} }}~{\text{dA}}=\mathbf {0} ~.$

From the moment balance equation above we can show that (see Appendix)

$\int _{\partial {\Omega }}\mathbf {x} \otimes {\bar {\mathbf {t} }}~{\text{dA}}=\int _{\partial {\Omega }}{\bar {\mathbf {t} }}\otimes \mathbf {x} ~{\text{dA}}~.$

Therefore the average stress tensor $\langle {\boldsymbol {\sigma }}\rangle$ is symmetric if the applied tractions are self equilibrated.

Now, if we translate the body by a constant amount $\mathbf {u} _{0}$ (rigid body translation), we get

${\bar {\langle {\boldsymbol {\sigma }}\rangle }}={\cfrac {1}{V}}~\int _{\partial {\Omega }}(\mathbf {x} +\mathbf {u} _{0})\otimes {\bar {\mathbf {t} }}~{\text{dA}}={\cfrac {1}{V}}~\int _{\partial {\Omega }}[\mathbf {x} \otimes {\bar {\mathbf {t} }}+\mathbf {u} _{0}\otimes {\bar {\mathbf {t} }}]~{\text{dA}}~.$

or

${\bar {\langle {\boldsymbol {\sigma }}\rangle }}={\cfrac {1}{V}}\left[\int _{\partial {\Omega }}\mathbf {x} \otimes {\bar {\mathbf {t} }}~{\text{dA}}+\mathbf {u} _{0}\otimes \int _{\partial {\Omega }}{\bar {\mathbf {t} }}~{\text{dA}}\right]=\langle {\boldsymbol {\sigma }}\rangle$

Therefore, the average stress is not affected by a rigid body translation only if the applied tractions are self equilibrated.

We can conclude that the average stress $\langle {\boldsymbol {\sigma }}\rangle$ is an acceptable measure of stress in a RVE if the applied tractions are self equilibrated.

← Average displacement in a RVE

Micromechanics of composites

Average stress power in a RVE →