Exact Admissibility and Spectral Ordering of Unequal-Order Four-Field Polynomial Discretizations for Free-Free Timoshenko Beams

Imagine a long, slender bridge or a diving board. When engineers design these structures, they often rely on a simplified model that assumes the material is so thin and stiff that it bends without shearing, like a perfectly rigid ruler. This works well for very thin beams, but when the beam is shorter, thicker, or made of a material that flexes easily, that simple model breaks down. It fails to account for how the material slides against itself as it bends, a phenomenon known as shear deformation. To handle these real-world scenarios, scientists use a more complex model called the Timoshenko beam theory. This model treats the beam as having two distinct behaviors: the up-and-down movement of the beam itself, and the rotation of its cross-section. It also tracks the internal forces, specifically the shear force that tries to slide the layers past one another and the bending moment that causes the curve.

The challenge arises when scientists try to solve the equations for these beams using computers. They must break the continuous beam down into a finite set of mathematical pieces, a process called discretization. For the Timoshenko model, this is tricky because the different parts of the equation—the movement, the rotation, and the internal forces—must be approximated with just the right amount of mathematical detail. If the computer model uses too little detail for the forces, or too much for the movement, the simulation can produce "ghost" results: vibrations that shouldn't exist, or a failure to recognize that the beam can simply float or slide without bending. These errors are not just small mistakes; they can make the entire simulation useless. For decades, engineers have known that the choice of mathematical detail matters, but they have lacked a precise rulebook for how to mix and match these details when using a specific type of high-precision calculation based on global polynomials, which are smooth curves used to approximate the beam's shape across its entire length.

In a recent study, researchers Ming-Hung Shu and You-Pei Li tackled this specific puzzle. They investigated a four-field approach where the displacement, rotation, shear force, and bending moment are all represented by independent mathematical curves. The central question was simple yet profound: if you choose a certain level of complexity for the movement and rotation, exactly how much complexity do you need for the forces to get a correct answer? They discovered that the answer is not a simple rule of thumb, such as "the forces must always be calculated with one degree more detail than the movement." Instead, the relationship is a complex, piecewise map where the requirements for the bending force depend on the specific choices made for the shear force and the rotation.

The researchers found that the system has a hidden mechanism of compensation. If you increase the mathematical detail used to calculate the shear force, you can actually get away with using a simpler, lower-order curve for the bending moment. This cross-field compensation means that the four variables are deeply interconnected; you cannot judge the stability of the model by looking at the shear and bending pairs in isolation. The team derived an exact rule that tells you, for any combination of four different levels of complexity, whether the model will work or fail. They proved that if the model is set up correctly, it will capture exactly two ways the beam can move without bending—sliding sideways and rotating as a whole—and no other "ghost" movements. If the setup is wrong, the model creates extra, fake zero-energy modes that look like vibrations but are actually mathematical artifacts.

To verify these findings, the authors did not rely on standard computer approximations that can be prone to rounding errors. Instead, they used exact algebraic methods, treating the numbers as precise fractions rather than floating-point decimals. They tested thousands of different combinations of polynomial orders. In every case, their exact rule predicted the outcome perfectly. They also examined what happens when you add more detail to the force calculations while keeping the movement calculation fixed. They showed that the calculated frequencies of the beam's vibration will always increase or stay the same as you add more detail to the forces, but they will never exceed a specific upper limit defined by a different, simpler type of calculation. Crucially, they demonstrated that getting closer to this upper limit does not always mean getting closer to the true physical reality. Sometimes, a less detailed model might accidentally be closer to the true answer because two different types of errors cancel each other out, while a more detailed model moves away from the true answer as it corrects one error but exposes another.

The study also highlighted a practical danger in how these calculations are performed on computers. The mathematical equations used to describe the beam have a peculiar structure where some parts are singular, meaning they can produce infinite values that are not physical. When researchers try to solve the full set of equations directly, the computer's choice of mathematical coordinates can drastically affect the result. The researchers showed that using a standard set of polynomial coordinates can cause the computer to misclassify these infinite values, leading it to report a high-frequency vibration that does not actually exist. However, by switching to a different, more stable set of coordinates known as shifted-Legendre polynomials, the computer correctly separates the real vibrations from the mathematical noise. This finding serves as a warning: even if the underlying math is correct, the way the computer is asked to solve it can introduce false results.

Ultimately, this work provides a rigorous, exact map for engineers and scientists who use these advanced polynomial methods to model beams. It replaces vague guidelines with a precise boundary that defines when a model is trustworthy. The researchers proved that the stability of the simulation depends on the specific interplay between all four fields, not just a simple comparison of their individual complexities. By establishing these exact conditions, they ensure that future simulations of free-floating beams will correctly identify the physical modes of vibration and avoid the trap of spurious, non-physical results. The study confirms that while the mathematics of these beams is intricate, the rules governing their stability are exact and discoverable, provided one looks at the system as a whole rather than in isolated parts.