The in-plane equilibrium equations are and . On a simply connected patch, each divergence-free row has a streamfunction: write , , . Symmetry gives , so another potential exists with , . Hence
If two such Airy stress functions give the same stress, their difference has zero Hessian and is therefore an arbitrary affine function on a connected domain. This is the full nonuniqueness for a fixed stress. The existence argument is local, or global on a simply connected domain; domains that are multiply connected regions require the associated force and moment period conditions for a single-valued global potential. Equilibrium alone should not silently erase that topological qualification.
For plane strain put . The constitutive equation gives . Consequently
The displacement Saint-Venant compatibility equations then becomes
For stable elastic moduli the coefficient is nonzero, so . If is bounded and smooth and on its boundary, integrating twice by parts gives . Thus is harmonic with zero boundary value and by uniqueness of the Dirichlet problem. The normal-derivative condition is what removes the extra fourth-order boundary terms.
For the pressurized annulus, substituting the radial constitutive equations into equilibrium gives , so . The stresses are
Using gauge stresses relative to ambient, the boundaries are , . Therefore
and
The maximum occurs at the bore. Once the chosen yield criterion is reached, Hooke's law cannot describe the total strain alone: write total strain as elastic plus plastic strain and supplement the elastic stress law with a yield condition, a plastic flow rule and any required work hardening law. A plastic region starts near the bore and must be matched to the remaining elastic region. The exact conversion between a principal-stress difference and the quoted yield stress depends on whether the material criterion is Tresca or von Mises. This is the elastic stress concentration in a pressurized thick cylinder.