Past exam of the mathematics course of the University of Cambridge / 2007 / iii / Paper 76 (Nonlinear Continuum Mechanics)

Use both material isotropy and material frame indifference. The polar decomposition in continuum mechanics writes $F=VR$, where $V$ is the left stretch tensor. Isotropy removes the reference rotation: $\sigma (F)=\sigma (V)$. Choose an orthonormal eigenbasis of $V$. A half-turn about any principal axis leaves $V$ fixed but reverses appropriate off-diagonal entries of $\sigma (V)$. Frame indifference therefore makes those entries zero. If a principal stretch is repeated, rotations in its eigenspace additionally force the stress to be isotropic there. Thus the principal stresses share the stretch axes: this proves coaxiality of isotropic elastic stress with $V$ and with $B=V^2=F{F}^T$.

For simple shear, the left Cauchy-Green deformation tensor isIts principal axes make an angle $\vartheta$ with the coordinate axes satisfying $2\cot (2\vartheta )=\gamma$ when $\gamma \ne 0$. If the two principal stresses are $s_1,s_2$, rotating their diagonal matrix back givesConsequently the universal simple-shear normal-stress identity isEquivalently, it follows immediately from the off-diagonal entry of $\sigma B-B\sigma =0$. At $\gamma =0$, the undeformed material has isotropic stress, so the identity holds there too. No particular elastic constitutive equation has been used.

The inextensible fibre constraint gives $|{\mathbf{x}}_{,1}|=1$. Write the fibre direction and its transverse unit direction asThen $x_{1,1}=\cos \theta$ and $x_{2,1}=\sin \theta$. In plane strain, incompressibility gives $\det F=1$. Decomposing the second column of the deformation gradient along $e,n$ therefore gives ${\mathbf{x}}_{,2}=\gamma e+n$, where $\gamma ={\mathbf{x}}_{,2}\cdot e$. HenceThis is a local rotation followed by simple shear.

The compatibility of the deformation map requires ${\mathbf{x}}_{,12}={\mathbf{x}}_{,21}$. Since $e_{,j}={\theta }_{,j}n$ and $n_{,j}=-{\theta }_{,j}e$, this isTaking components gives the compatibility relations for incompressible inextensible-fibre plane strain:Along a reference curve,Thus $\theta$ is constant when $d{X}_2/d{X}_1=-1/\gamma$. The equivalent condition $d{X}_1+\gamma d{X}_2=0$ also covers $\gamma =0$. Its current tangent is proportional to $F(-\gamma ,1{)}^T=n$, so these transverse curves are straight wherever the deformation is regular.

To construct the stress with zero body force, specify the local shear stress $h$ as the material's function or functional of $\gamma$ under simple shear. The normal stresses are constraint reactions: incompressibility supplies pressure, and the inextensible fibre constraint supplies fibre tension. WriteThe compatible geometry has $\nabla \theta =\kappa e$ with $\kappa ={\theta }_{,1}$, so $\nabla \cdot e=0$, $\nabla \cdot n=-\kappa$, $(e\cdot \nabla )e=\kappa n$ and $(n\cdot \nabla )e=0$. Resolving $\nabla \cdot \sigma =0$ along the two directions givesThe triangular equilibrium construction for inextensible-fibre plane strain now determines $s$ by integration along fibres and then $t$ along transverse curves. Thus any regular compatible deformation has a local zero-body-force stress construction once the shear law and suitable boundary integration data are supplied. There is no additional constitutive normal-stress law to specify: the two normal components are fixed by force balance and boundary tractions. This does not assert that arbitrary separately prescribed boundary tractions can accompany an arbitrary deformation.

Put $c=a^2-a_0^2$ and use reference cylindrical coordinates $(\rho ,\Theta ,Z)$ and current coordinates $(r,\theta ,z)$. The deformation map is $r=\sqrt{{\rho }^2+c}$, $\theta =\Theta +\alpha Z$, $z=Z$. In the associated orthonormal bases its deformation gradient isIts determinant is one, as required by incompressibility. The radial principal stretch is ${\lambda }_0=\rho /r$. The other two squared principal stretches are the eigenvalues of the tangential-axial block of $F{F}^T$:Writing $K=(r/\rho {)}^2+{\alpha }^2{r}^2+1$, they areThe positive square roots give ${\lambda }_{\pm }$, and ${\lambda }_0{\lambda }_{+}{\lambda }_{-}=1$.

The strain energy density is per reference volume, so the energy per height isThe work rates per height are $M\dot{\alpha }$ from the end torque and $2\pi ap\dot{a}$ from the inner pressure. The outer surface contributes no work because it is traction free, and the fixed height has zero axial work. Thus the virtual work identities for twist and inflation of an incompressible tube are

For a neo-Hookean solid, an irrelevant energy constant can be omitted:Differentiating under the finite integral givesAn antiderivative divided by $\mu$ is $\log \rho -\frac{1}{2}\log ({\rho }^2+c)-c/[2({\rho }^2+c)]+{\alpha }^2{\rho }^2/2$. With $b^2=b_0^2+c$, the neo-Hookean tube pressure under twist and inflation is thereforeIt vanishes for an untwisted undeformed tube, and the twist contributes the last term at fixed radii.

Use the material-first nominal stress tensor convention $P_{Ii}$, with $F_{iI}$ the deformation gradient. For any fixed reference material subvolume $V_0$ with outward unit normal $N$, conservation of energy isHere $b$ is body force per mass, $r$ is heat supply per mass and $q$ is nominal heat flux. Using momentum conservation to remove the kinetic energy balance gives the integral Lagrangian internal energy balance:The integral Lagrangian entropy inequality iswhere $\theta >0$ is temperature and $q/\theta$ is the thermodynamic entropy flux.

Since the subvolume is arbitrary, the divergence theorem givesExpand $\dot{U}=U_{,F_{iI}}{\dot{F}}_{iI}+U_{,\eta }\dot{\eta }+U_{,{\xi }_r}{\dot{\xi }}_r$ and eliminate the heat supply between these relations. The Clausius-Duhem inequality becomesUnder the usual local constitutive independence assumptions, the reversible mechanical and thermal variations can have either sign. The Coleman-Noll procedure therefore imposesDefine the thermodynamic force conjugate to an internal variable by $f_r=-{\rho }_0{U}_{,{\xi }_r}$. The energy balance and the remaining entropy constraint now give, respectively,

The specific Helmholtz free energy $\psi =U-\theta \eta$ satisfies $d\psi ={\rho }_0^{-1}{P}_{Ii}d{F}_{iI}-\eta d\theta -{\rho }_0^{-1}{f}_{r}d{\xi }_r$. For the multiplicative decomposition of the deformation gradient $F=F_{\ast }A$, adopt $P_{Ji}^{\ast }={\rho }_0\partial \psi /\partial (F_{\ast }{)}_{iJ}$, consistently with the material-first nominal stress tensor. Differentiating $F_{\ast }=F{A}^{-1}$ givesConsequentlyorThese relations follow from the chain rule; the transpose convention matters.

Differentiating the dissipation potential, with the Frobenius norm of $Q$, givesAn integrating factor yields the exponential memory from tensor relaxation:For $n>0$, the driving term is continuously zero at $Q=0$.

Taking the product $T=PF$ explicitly defined in the question, $F=F_{\ast }A$ gives $T=QA$, and hence its requested history functional isThere is a naming error in the PDF: with its $P_{Ii}$ convention, the actual second Piola-Kirchhoff stress tensor is $S=P{F}^{-T}$, while $PF=SC$ is the transpose of the Mandel stress tensor, where $C=F^{T}F$. Thus $S=QA{C}^{-1}$ if that stress measure is required. This agrees with the stress push-forward used in Question 5; the printed product is fully determined by the history above, whereas recovering $S$ also requires the deformation or an invertible elastic constitutive relation.

The multiplicative decomposition of the deformation gradient is used on histories where $A$ stays invertible. Thermodynamic admissibility additionally requires $\mathrm{tr}(Q\dot{A})=\alpha \Vert Q{\Vert }^{n+1}-\mathrm{tr}(AQ)/\tau \ge 0$: positivity of the constants alone does not ensure this for arbitrary prescribed histories.

Write $\mathbf{\tau }$ for the Kirchhoff stress tensor and $S$ for the second Piola-Kirchhoff stress tensor. Since $\mathbf{\tau }=FS{F}^T$ and $\dot{F}=LF$ for the velocity gradient,Therefore the upper-convected rate as a stress push-forward is

For the simple shear flow, put $\gamma =f_{,2}$ and $g=\dot{\gamma }$. The velocity gradient and rate-of-strain tensor areTheir upper-convected derivative givesDenote the symmetric extra stress by $E={\sigma }^d$, and the relaxation time by ${\tau }_r$ to distinguish it from $\mathbf{\tau }$. The Oldroyd-B model in the question becomesFor the steady constitutive response, $E_{22}=E_{33}=E_{13}=E_{23}=0$, $E_{12}=\mu g$ and $E_{11}=2(\mu -{\mu }_r){\tau }_r{g}^2$. HenceA constant imposed shear rate can still have startup transients; these values describe the steady response.

For circular Couette flow $v(r)e_{\theta }$, the velocity gradient in the orthonormal polar basis $(e_r,e_{\theta })$ isFor completeness, in fixed Cartesian coordinates let $\omega (r)=v(r)/r$ and $J=\left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)$. Then $v=\omega Jx$ and $L=\omega J+({\omega }^{\prime }/r)(Jx)\otimes x$. Along a particle, the polar basis rotates at angular speed $\omega$. Its contribution to the material derivative cancels the rigid-rotation part of $L$ in the upper-convected derivative. The rotating-frame reduction of circular viscoelastic shear therefore leavesThe flow direction is now $e_{\theta }$ and the gradient direction is $e_r$. With $g=v^{\prime }-v/r$, the same simple-shear constitutive calculation gives

Angular force balance implies $r^2\mu g$ is constant. Integrating $v^{\prime }-v/r\propto r^{-2}$ gives $v=Ar+B/r$. With $v(a)=V_a$ and $v(b)=0$,The radial force balance is $-p^{\prime }-2(\mu -{\mu }_r){\tau }_r{g}^2/r=0$. Thus the inertialess Oldroyd-B circular Couette pressure isThe constant $p_0$ is fixed by a pressure or normal-traction boundary condition.

Write $c=\cos 2\varphi$, $s=\sin 2\varphi$ and $\sigma =({\sigma }_{11}+{\sigma }_{22})/2$. On the fixed yield surface, ${\xi }_{,l}=-c$ and ${\tau }_{,l}=-s$, where $l$ is stress-space arc length. Substituting ${\sigma }_{11}=\sigma +\xi$, ${\sigma }_{22}=\sigma -\xi$, ${\sigma }_{12}=\tau$ into planar force balance givesThus $\nabla \sigma =M\nabla l$, with $M=\left(\begin{array}{cc}c& s\\ s& -c\end{array}\right)$.

For a spatial curve $x(s_0)$ and a function $\mathcal{F}(\sigma ,l)$, the chain rule givesHere primes refer to the spatial parameter $s_0$, not the yield-curve parameter $l$. The displayed row-matrix product being zero is therefore sufficient for $\mathcal{F}$ to remain constant.

The eigenvectors of $M$ are $e_{\alpha }=(\cos \varphi ,\sin \varphi )$ and $e_{\beta }=(-\sin \varphi ,\cos \varphi )$, with eigenvalues $+1$ and $-1$. Along these slip lines, respectively, $d\sigma =dl$ and $d\sigma =-dl$. This proves the stress characteristics for an anisotropic yield curve:

For the velocity relations, the outward normal to the reduced yield curve is proportional to $(-\sin 2\varphi ,\cos 2\varphi )$. The associated flow rule, using the symmetric stress work pairing $D:\delta \sigma$, therefore gives a rate-of-strain tensor proportional toThe off-diagonal stress variation contributes $2D_{12}\delta {\sigma }_{12}$, which ensures this normality calculation uses the correct factor. Direct multiplication gives $e_{\alpha }^{T}D{e}_{\alpha }=e_{\beta }^{T}D{e}_{\beta }=0$: neither slip line direction extends instantaneously.

Resolve the velocity as $w=u{e}_{\alpha }+v{e}_{\beta }$. Since $d{e}_{\alpha }=e_{\beta }d\varphi$ and $d{e}_{\beta }=-e_{\alpha }d\varphi$, its longitudinal derivative along an $\alpha$ curve is $du-vd\varphi$, and along a $\beta$ curve it is $dv+ud\varphi$. Setting these extensions to zero gives the Geiringer velocity relations: