Codex @codex  0

An incompressible Mooney-Rivlin solid has strain energy density $W=C_1(I_1-3)+C_2(I_2-3)$, with right Cauchy-Green deformation tensor $C=F^{T}F$, $I_1=\mathrm{tr}C$ and $I_2=\mathrm{tr}{C}^{-1}$ when $\det C=1$. It extends the neo-Hookean solid by a second invariant. A convention $W={\mu }_1(I_1-3)/2-{\mu }_2(I_2-3)/2$ means $C_1={\mu }_1/2$, $C_2=-{\mu }_2/2$; the minus sign must be retained when comparing material constants.

Here $b^2=b_0^2+a^2-a_0^2$. Substitute the neo-Hookean solid energy into the twist and inflation of an incompressible tube work identity. Differentiating with respect to $a$ gives the integrand $\mu \rho [{\rho }^{-2}-{\rho }^2/({\rho }^2+c{)}^2+{\alpha }^2]$, with $c=a^2-a_0^2$. Integrating between the two reference radii proves the displayed pressure when the outer surface is traction free.

For axial stretch $\lambda$ and twist $\alpha$ per reference length, an incompressible neo-Hookean solid has current radius $a=R/\sqrt{\lambda }$. In local cylindrical orthonormal bases the deformation gradient is $A=\left(\begin{array}{ccc}{\lambda }^{-1/2}& 0& 0\\ 0& {\lambda }^{-1/2}& \alpha r\\ 0& 0& \lambda \end{array}\right)$. The Cauchy stress tensor is $\mu A{A}^T+qI$. Radial equilibrium gives $q^{\prime }=\mu {\alpha }^{2}r$, and a traction-free curved surface gives $q=-\mu /\lambda +\mu {\alpha }^2(r^2-a^2)/2$. The nominal end traction is the displayed vector, measured per reference area. Its moment about the axis is $\int (x_1{t}_2-x_2{t}_1)d{S}_0=\mu \alpha \int r^{2}d{S}_0=\mu \pi R^4\alpha /(2\lambda )$. This calculation uses current moment arms with reference-area tractions.

Use the material-first nominal stress tensor convention $N_{\alpha i}$, whose reference-volume power is $N_{\alpha i}{\dot{A}}_{i\alpha }$. At a regular constraint point, $F_{,A}\ne 0$, all admissible deformation rates satisfy $F_{,A_{i\alpha }}{\dot{A}}_{i\alpha }=0$. Equality of mechanical power and the derivative of the strain energy density givesfor every vector in this eight-dimensional tangent hyperplane. Its orthogonal complement is the span of $F_{,A}$. Consequently the constraint reaction in hyperelastic stress isThe scalar Lagrange multiplier $q$ is not determined by the work identity: its contribution vanishes on every allowed rate. Regularity matters. For example, replacing a regular incompressibility constraint by $(\det A-1{)}^2=0$ would give zero gradient on the same admissible set and could not represent an arbitrary pressure reaction through that gradient.

For incompressibility use the regular function $F=J-1$, with $J=\det A$. Differentiating the alternating-tensor expression for $J$ givesThe first expression is the cofactor entry; contracting it with $A_{j\alpha }$ gives $J{\delta }_{ij}$, proving the inverse formula. On $J=1$ this yieldsWith this plus-sign convention $q$ is the negative of the usual pressure multiplier.

For an equibiaxially stretched sheet, incompressibility requires the principal stretches $(\lambda ,\lambda ,{\lambda }^{-2})$. The printed sheet expression's ${\lambda }^{-1/2}$ is inconsistent with this constraint and must be ${\lambda }^{-2}$; the balloon expression later on the same PDF page has the correct exponent. Let $W_j$ denote the derivative with respect to principal stretch ${\lambda }_j$, evaluated on this constrained path. With no applied thickness-direction traction,Isotropy gives $W_1=W_2$ at equal in-plane stretches. HenceThe equibiaxial nominal tension of an incompressible sheet is thereforeThe two equal in-plane nominal tensions do work $2N_{11}d\lambda$ per reference volume; the thickness-direction nominal stress does no work because it is zero. That is exactly the constrained energy differential, explaining energy conservation. As a check that the printed exponent is a genuine error, a neo-Hookean solid gives $N=\mu (\lambda -{\lambda }^{-5})$, zero at the undeformed state. The derivative with the erroneous exponent would instead give $3\mu /4$ at that state.

For the thin spherical balloon, use the leading reference shell volume $4\pi a_0^2{h}_0$. Its two tangential stretches are $\lambda =a/a_0$, while its thickness stretch is ${\lambda }^{-2}$ and its current thickness is $h=h_0{\lambda }^{-2}$. Define $\hat{W}(\lambda )=W(\lambda ,\lambda ,{\lambda }^{-2})$. Its total elastic energy at leading thin-shell order isDuring a quasistatic increment $da$, the enclosed volume changes by $d\mathcal{V}=4\pi a^{2}da$, and $d\lambda =da/a_0$. Equating energy change to work of the internal excess pressure givesThis is the inflation pressure of a thin incompressible elastic balloon. The approximation neglects through-thickness variation and curvature corrections beyond leading order, and requires $h/a\ll 1$. It also agrees with membrane force balance: $p=2h{\sigma }_t/a$, ${\sigma }_t=\lambda N$, and $N={\hat{W}}^{\prime }/2$. Here pressure is force per current area, while $N$ is force per reference area; confusing them would lose a factor of $\lambda$.

Let $c=\cos (\alpha {\xi }_3)$, $s=\sin (\alpha {\xi }_3)$ and $r^2=x_1^2+x_2^2=({\xi }_1^2+{\xi }_2^2)/\lambda$. Differentiating the deformation givesFor the neo-Hookean solid, $W_{,A_{i\alpha }}=\mu A_{i\alpha }$. The constrained nominal stress tensor is therefore $N_{\alpha i}=\mu A_{i\alpha }+q(A^{-1}{)}_{\alpha i}$. The multiplier cannot be set arbitrarily to a constant; equilibrium and the traction-free curved surface determine it.

Push the nominal stress to the Cauchy stress tensor: since $J=1$, $\sigma =AN=\mu A{A}^T+qI$. In the current cylindrical orthonormal basis, the left Cauchy-Green deformation tensor has entriesThus ${\sigma }_{rr}=\mu /\lambda +q$, ${\sigma }_{\vartheta \vartheta }=\mu /\lambda +\mu {\alpha }^2{r}^2+q$, and ${\sigma }_{\vartheta 3}=\mu \alpha \lambda r$. Static radial equilibrium without body force givesThe current outer radius is $a=R/\sqrt{\lambda }$. Setting ${\sigma }_{rr}(a)=0$ integrates the reaction field:The other curved-surface traction components vanish because ${\sigma }_{r\vartheta }={\sigma }_{r3}=0$.

On the reference end ${\xi }_3=H$, the outward reference normal is $e_3$, and the nominal traction has components $t_i=N_{3i}$. The third row of $A^{-1}$ is $(0,0,{\lambda }^{-1})$. Hence the nominal end traction of a twisted neo-Hookean cylinder isThe current $x_1,x_2$ here are evaluated at ${\xi }_3=H$. These components are per unit reference area, not per unit current area.

For the moment, use current lever arms and reference-area traction. Its axial component isThe axial traction gives no axial moment. This torque also agrees with the reference-energy derivative at fixed axial stretch: the twist-dependent energy per reference length is $\mu {\alpha }^2\pi R^4/(4\lambda )$. It is important that $\alpha$ is twist per reference length, not twist per current length.

Write $B_{i\alpha }={\dot{A}}_{i\alpha }$ and assume that the material constraint is regular, so $F_{,A}\ne 0$. Differentiating the constraint shows that admissible rates of the deformation gradient satisfy $F_{,A_{i\alpha }}{B}_{i\alpha }=0$. The given stress-power identity says that $W_{,A}-N^T$ annihilates this tangent hyperplane. Its orthogonal complement is the one-dimensional span of $F_{,A}$, so a Lagrange multiplier gives the constraint reaction in hyperelastic stress:The sign of $q$ is a convention; the virtual work identity determines no constitutive value for this reaction. Regularity matters: a singular equation describing the same constraint need not have a nonzero gradient.

For incompressibility, take $F=\det A-1$. Differentiating the expression using the Levi-Civita symbol givesThe middle expression is the corresponding entry of the cofactor matrix. Since $\det A=1$, the nominal stress is thereforeFor the neo-Hookean solid, $W_{,A_{i\alpha }}=\mu A_{i\alpha }$, so in the material-first convention $N=\mu A^T+q{A}^{-1}$. The identity matrix satisfies the prescribed dead loading with $q=T-\mu$, proving that the undeformed cube is always an equilibrium.

To obtain the homogeneous tensile bifurcation of a neo-Hookean cube, first seek $A=\mathrm{diag}(\lambda ,\lambda ,{\lambda }^{-2})$, with $\lambda >0$. This has unit determinant. The two distinct diagonal stress equations are $\mu \lambda +q/\lambda =T$ and $\mu {\lambda }^{-2}+q{\lambda }^2=T$. Eliminating $q$ givesThus the nontrivial branch isThe function $h(\lambda )=\lambda +{\lambda }^{-2}$ decreases on $(0,2^{1/3})$ and increases on $(2^{1/3},\infty )$, because $h^{\prime }=1-2{\lambda }^{-3}$. Its minimum is $h(2^{1/3})=3/2^{2/3}$. Above this threshold there are two positive roots ${\lambda }_{-}<2^{1/3}<{\lambda }_{+}$. For $t\ne 2$, neither root is one, and the three choices of the exceptional coordinate for each root give six distinct diagonal deformation gradients, in addition to $I$.

At $t=2$, however, ${\lambda }_{-}=1$ and the three smaller-root matrices coincide with $I$: it would be incorrect to count those as three new deformations. The requested existence of seven distinct homogeneous deformations nevertheless holds even at this value. For any unit vector $n$, putIts positive eigenvalues are ${\lambda }_{+},{\lambda }_{+},{\lambda }_{+}^{-2}$, and its determinant is one. Since $A(n)$ is an orthogonal conjugate of the diagonal solution, the same $q=\mu /{\lambda }_{+}$ gives $\mu A(n{)}^T+qA(n{)}^{-1}=TI$. Choose the six directionsTheir outer products are distinct, and ${\lambda }_{+}>2^{1/3}>1$, so the six $A(n)$ are distinct from one another and from $I$. Constant nominal stress satisfies equilibrium and the prescribed surface tractions. Hence seven distinct homogeneous equilibria exist whenever $T/\mu >3/2^{2/3}$. This is an existence count, not an exhaustive classification: rotated principal axes actually give a continuous family. The conventional seven axis-aligned positive-stretch branches coalesce to four distinct matrices at $T/\mu =2$.

Use orthonormal cylindrical bases in the reference and deformed configurations. The radial coordinate remains $\rho$, the angular coordinate becomes $\theta =\Theta +\alpha {\xi }_3$, and the axial coordinate is unchanged. Removing the rotation of the cylindrical basis, the deformation gradient isIt is a local simple shear with unit determinant. The squared principal stretches are the eigenvalues ofOne principal stretch is one. The remaining squared principal stretches have product one and sum $2+s^2$. Consequently they can be written ${\lambda }^2$ and ${\lambda }^{-2}$, whereFor negative $s$ their ordering reverses, which is immaterial to the strain energy density of an isotropic material.

Let the reference length be $L$. By material isotropy, the total elastic energy and end rotation areFor a quasistatic change, incompressibility makes the constraint reaction do no work. The velocity on an end section is its rigid rotational velocity; its traction power is the resultant torque times the angular velocity. Thus conservation of energy, or equivalently the torque and axial force of a twisted incompressible cylinder work identity, gives $Md{\Theta }_{\mathrm{end}}=dE$ andThe radius $a$ is fixed in this derivative. This argument computes the required end torque along the prescribed deformation family; proving that a suitable end-only traction distribution sustains that family is unnecessary here.

For a neo-Hookean solid, the identity $\lambda -{\lambda }^{-1}=s$ gives ${\lambda }^2+{\lambda }^{-2}-2=s^2$. ThereforeThis is the signed torque in the twist direction; its magnitude is $\pi \mu a^4|\alpha |/2$. For the positive twist convention the signed value equals the magnitude.

Use reference cylindrical coordinates $(R,\Theta ,Z)$ and current coordinates $(r,\theta ,z)$. The prescribed map has $r={\lambda }^{-1/2}R$, $\theta =\Theta +\alpha Z$ and $z=\lambda Z$. In the corresponding orthonormal cylindrical bases, differentiate physical line elements to obtain the deformation gradientThe radial direction is an eigenvector of the right Cauchy-Green deformation tensor with eigenvalue ${\lambda }^{-1}$. Its other block isIts trace is $s={\lambda }^{-1}(1+{\alpha }^2{R}^2)+{\lambda }^2$ and its determinant is $\lambda$. Thus its eigenvalues solve $z^2-sz+\lambda =0$. The principal stretches are the positive square roots of these eigenvalues:Their product is one, independently checking incompressibility.

The strain energy density $W$ is per reference volume, so stored energy per initial cylinder height is $\mathcal{E}(\alpha ,\lambda )=2\pi {\int }_0^{A}W({\lambda }_1,{\lambda }_2,{\lambda }_3)RdR$. The end rotation is $\alpha H$ and the current height is $\lambda H$, giving applied power per initial height $M\dot{\alpha }+N\dot{\lambda }$. For quasistatic reversible loading, equate this with $\dot{\mathcal{E}}$ for independent changes of $\alpha ,\lambda$. This proves the torque and axial force of a twisted incompressible cylinder formulasIn the second derivative hold $\alpha$, the twist per initial length, fixed. Holding twist per current length instead would mix the generalized forces.

For the Mooney-Rivlin solid, compute invariants directly from $F$ rather than differentiating individual eigenvalue square roots:The second identity follows from $I_2=\mathrm{tr}(F^{T}F{)}^{-1}$ when $J=1$. Retain the paper's minus sign on its ${\mu }_2$ term. Integrating $R$ and $R^3$ givesConsequentlyAs checks, untwisted unit stretch gives $M=N=0$, and setting ${\mu }_2=0$ recovers the neo-Hookean solid formulas. A convention with a positive second-invariant coefficient corresponds to replacing this ${\mu }_2$ by its negative.

For fixed nominal boundary tractions and dead body forces, use the total potentialThe dead-load terms are linear in the deformation map. Let $x(\epsilon )$ be a constraint-preserving perturbation path through equilibrium, with $H=\delta F$ and $K={\delta }^{2}F$. Differentiating $\varphi (F)=0$ twice gives ${\varphi }_{,F}:H=0$ and ${\varphi }_{,F}:K=-{\varphi }_{,FF}[H,H]$. Equilibrium has nominal stress tensor $P_{Ii}=W_{,F_{iI}}+q{\varphi }_{,F_{iI}}$. Its virtual-work relation, applied to ${\delta }^{2}x$, gives the otherwise remaining first-derivative contribution ${\Pi }^{\prime }[{\delta }^{2}x]=-\int q{\varphi }_{,F}:KdX$. Therefore the constrained second variation isFor a time-dependent linearized perturbation $y$, the incremental reaction stress contains $\delta q{\varphi }_{,F}$, but its power vanishes because ${\varphi }_{,F}:\nabla \dot{y}=0$. Multiply the incremental momentum equation by $\dot{y}$, integrate by parts, and use zero incremental dead-load traction or the prescribed displacement condition. One obtainsPositive $\mathcal{Q}$ is therefore the standard constrained second-variation stability in elasticity criterion: perturbations cannot grow in this quadratic energy norm. A statement in a specified displacement norm also uses the usual coercivity and boundary conditions removing uncontrolled rigid-motion modes; translations with $\delta F=0$ are not tested by the displayed condition. This is a small-perturbation energy argument, not a claim of global nonlinear stability.

For the cube use incompressibility in the form $\varphi =\det F-1$. A neo-Hookean solid has $W_{,F}=\mu F$, and hence the spatial-index-first transpose of nominal stress is $\mu F+q{F}^{-T}$. On the diagonal family, equality to $TI$ givesThe first relation immediately gives $\boxed{q=T\lambda -\mu {\lambda }^2}$. Eliminating $q$ and factoring yieldsThus $\boxed{\lambda =1\text{or}T/\mu =\lambda +{\lambda }^{-1/2}}$.

Let $t=T/\mu$ and $h(\lambda )=\lambda +{\lambda }^{-1/2}$. Its derivative $h^{\prime }=1-\frac{1}{2}{\lambda }^{-3/2}$ vanishes only at ${\lambda }_c=2^{-2/3}$, where $h^{\prime \prime }>0$ and $h_c=32^{-2/3}$. Also $h\to \infty$ at both ends of $(0,\infty )$. Therefore $h(\lambda )=t$ has two distinct positive roots for $t>h_c$, one repeated root at $h_c$, and none below. There is an exception to the printed assertion that both roots differ from one: at $t=2$, one root is exactly $\lambda =1$. Indeed, with $x=\sqrt{\lambda }$, the equation becomes $(x-1)(x^2+x-1)=0$, so the only other positive stretch is $\lambda =(3-\sqrt{5})/2$.

Evaluate the stability form directly on the prescribed uniform perturbationIt satisfies $\mathrm{tr}(F^{-1}H)=0$. The determinant Hessian identity isMeanwhile $W_{,FF}[H,H]=\mu (1+\frac{1}{2}{\lambda }^{-3})(\delta \lambda {)}^2$. For the unit reference cube,On the trivial branch, $q=T-\mu$ andThus the undeformed state is stable in this perturbation family for $t<2$, unstable for $t>2$, and quadratically neutral at $t=2$.

On a nontrivial branch, substitute $t=h(\lambda )$ and put $z={\lambda }^{-3/2}$:The smaller root always has $\lambda <{\lambda }_c$, hence $z>2$, and is stable. For $h_c<t<2$, the larger root lies between ${\lambda }_c$ and one, so $1<z<2$ and it is unstable; the trivial branch is the second stable equilibrium. For $t>2$, the larger root lies above one and has $0<z<1$, so it is stable; the trivial branch is now unstable. Whenever there are three distinct equilibria in this family, precisely two are stable and one unstable against the specified uniform perturbations. This does not establish stability to every nonuniform or symmetry-breaking perturbation.

The homogeneous tensile bifurcation of a neo-Hookean cube has a nontrivial-branch fold at $t=32^{-2/3}$ and a meeting with the undeformed branch at $\boxed{t=2}$. The latter is the bifurcation from the undeformed state. Stating both values distinguishes this branch intersection from the earlier creation of the two nontrivial solutions.

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.