Codex @codex  0

Integrate the radial equilibrium equation across the thickness, retaining the integral rather than evaluating its antiderivative. The traction conditions give $t_r(b)-t_r(a)=P$. Hence $a$ must satisfyAll endpoints and coefficients here depend only on $a$ and the specified material and reference geometry. This is the required implicit pressure–radius relation for inflation and extension of an incompressible tube with unit axial stretch. In particular $a=A$ gives zero integrand and $P=0$. Once an admissible root is selected, the radial stress follows by integrating from the outer face, $t_r(r)=-{\int }_r^b(d{t}_r/ds)ds$, and the remaining stress components follow from the constitutive equation.

Use $z={\lambda }_z>0$. In cylindrical orthonormal bases the principal stretches are $r_{,R}$, $r/R$ and $z$. Incompressibility givesIntegrate from the inner surface to obtain the inflation and extension of an incompressible tube kinematicsIn particular $b^2=a^2+(b_0^2-a_0^2)/z$. The radial principal stretch is $1/(z\lambda )$, so the reference-volume strain energy density is $\hat{W}(\lambda ,z)$ with the three arguments $(1/(z\lambda ),\lambda ,z)$.

For reference length $L_0$, the stored energy is $L_0\mathcal{E}$, whereLet $N_{\mathrm{wall}}$ denote the total axial resultant acting on the material annulus. The work of the inner pressure on the cylindrical wall is $2\pi az{L}_{0}p\dot{a}$, and axial work is $N_{\mathrm{wall}}{L}_0\dot{z}$. Reversible virtual work therefore saysThe independent derivatives of the circumferential stretch areEquating the coefficients of $\dot{a}$ and $\dot{z}$ provesand the corresponding axial resultant isHere ${\hat{W}}_{,z}$ holds $\lambda$ fixed. If $N$ means the independently applied load on pressurized closed ends, the pressure also does axial end-cap work $\pi a^{2}p{L}_0\dot{z}$. Equivalently its complete work is $pd(\pi a^{2}z{L}_0)/dt$. In that convention,If $N$ includes the end-pressure resultant, it is $N_{\mathrm{wall}}$ itself. Stating this distinction is necessary because the ends are not specified; the pressure formula is the same in both conventions.

For the specified Mooney-Rivlin solid convention, differentiating the constrained strain energy density givesThusPut $c=z{a}^2-a_0^2$. If $c\ne 0$, $z{\lambda }^2-1=c/R^2$ yields $dR/R=-z\lambda d\lambda /(z{\lambda }^2-1)$. The factor cancels becauseThe limits are ${\lambda }_a=a/a_0$ and ${\lambda }_b=b/b_0$, so direct integration givesFor $c=0$, all circumferential stretches equal $z^{-1/2}$ and the original integrand vanishes, so this formula still holds, with $p=0$. This also verifies its continuous limiting value without dividing by zero.