⚡ electrical engineering

Application of Gröbner basis methodology to geometrically nonlinear axisymmetric circular isotropic plates

This paper introduces a novel framework that combines weighted-residual approximations with Gröbner basis algebraic elimination to derive closed-form analytical solutions for the geometrically nonlinear behavior of axisymmetric circular isotropic plates, offering a computationally efficient alternative to traditional numerical methods while maintaining high accuracy.

Original authors: T. M. Harrell, Y. J. Liu, J. Peddieson

Published 2026-09-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: T. M. Harrell, Y. J. Liu, J. Peddieson

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Engineers have long relied on a set of rules to predict how solid objects behave under pressure. When a bridge carries a heavy truck or a floor supports a crowd, the materials bend. For small bends, the math is straightforward: if you double the weight, the bend doubles. This linear relationship is easy to calculate and has served the construction industry for centuries. However, nature does not always follow simple rules. When a structure bends significantly—when the curve becomes deep and the material stretches as much as it bends—the old math breaks down. The forces interact in complex, tangled ways that make the equations difficult to solve. For decades, the only way to handle these large, messy deformations was to use powerful computers to crunch numbers step-by-step, a process that gives an answer but hides the underlying logic.

A team of researchers has found a way to untangle these complex problems without relying solely on brute-force computing. They applied a sophisticated branch of mathematics, originally developed for pure algebra, to the physical world of bending plates. By treating the equations of physics as a puzzle of polynomials, they were able to strip away the complexity and reveal a clean, direct formula. This approach allows engineers to see exactly how a structure will react to a load, with every variable—such as thickness, material strength, and load size—written out clearly in a single expression. It is a shift from getting a number to understanding the relationship itself.

The study focused on circular plates, which are common in engineering, ranging from manhole covers to the floors of large tanks. The researchers examined three specific scenarios: a plate held tightly around its edges, a plate resting on a rim that cannot move sideways, and a plate that extends beyond its support like a balcony. In each case, the plate was subjected to a steady, even pressure. The goal was to predict how much the plate would bend and what kind of stress it would feel, especially when the bending was large enough to make the material stretch.

To solve this, the team first used a standard engineering technique to turn the difficult, continuous problem of bending into a set of algebraic equations. Think of this as translating a flowing river into a series of connected buckets; the water still flows, but it is now represented by discrete amounts. This translation resulted in a system of equations that were tightly linked together, where changing one value affected all the others. Traditionally, solving such a system requires a computer to guess and check until it finds the right answer. Instead, the researchers applied a method called a Gröbner basis. This is a systematic way of rearranging a set of equations so that they fall into a neat, step-by-step order. The first equation contains only one unknown, the next contains that one unknown plus a second, and so on. This transformation turns a tangled knot into a straight line that can be solved one step at a time.

The result was a set of closed-form solutions. These are complete, symbolic formulas that describe the behavior of the plate using the language of algebra. Unlike computer simulations that output a single number for a specific set of conditions, these formulas retain all the variables. An engineer can look at the formula and immediately see how changing the thickness of the plate or the strength of the material will change the outcome, without running a new simulation. The researchers tested these new formulas against detailed computer simulations and found them to be highly accurate. In cases where the plate bent significantly, the new formulas matched the computer results with an error of less than four percent.

The study also highlighted a crucial limit of this method. When the researchers tried to make the model even more precise by adding more variables to the equations, the resulting algebraic puzzle became too complex to solve with a simple formula. The math showed that once the equations reached a certain level of complexity, a clean, symbolic answer was no longer possible. However, the researchers found that the simpler version was already accurate enough for practical use. The improvement gained by adding extra complexity was tiny compared to the massive error introduced by using the old, linear theories, which could be off by more than eighty percent in these large-bending scenarios.

This work demonstrates that a powerful tool from abstract mathematics can be brought to bear on real-world structural problems. It offers a new way to derive analytical solutions for nonlinear mechanics, a field that has long been dominated by numerical approximations. By proving that these complex, large-deformation problems can be solved with explicit formulas, the researchers have provided engineers with a clearer window into the behavior of structures. The method is not limited to circular plates; it can be applied to any structural problem where the physics can be described by polynomial equations. This opens the door to faster, more insightful design processes across civil, mechanical, and aerospace engineering, where understanding the precise relationship between load and deformation is critical for safety and efficiency.

Technical Summary: Application of Gröbner Basis Methodology to Geometrically Nonlinear Axisymmetric Circular Isotropic Plates

Problem Statement
The paper addresses the challenge of obtaining analytical, closed-form solutions for geometrically nonlinear mechanics problems, specifically focusing on axisymmetric circular isotropic plates under static distributed loads. Traditional approaches to such problems, governed by Kirchhoff–von Kármán plate theory, typically rely on numerical methods (e.g., Finite Element Analysis) because the governing equations, when reduced via weighted-residual approximations like Rayleigh–Ritz or Galerkin methods, result in highly coupled systems of multivariate polynomial algebraic equations. While linear theories fail to accurately predict behavior at large deflections (specifically when transverse displacement exceeds half the plate thickness), existing nonlinear solutions often lack explicit symbolic forms, hindering efficient parametric studies and design investigations.

Methodology
The authors propose a systematic framework that combines weighted-residual approximations with Gröbner basis algebraic elimination to transform coupled nonlinear polynomial systems into decoupled forms. The methodology proceeds as follows:

Formulation: The total potential energy of the plate is defined using Kirchhoff–von Kármán theory. Displacement fields (in-plane $u_r$ and transverse $w$ ) are approximated using trial functions that satisfy boundary conditions.
Reduction: Substituting these trial functions into the energy functional and minimizing with respect to unknown coefficients yields a system of coupled multivariate polynomial equations.
Decoupling: The Gröbner basis methodology (specifically using an elimination order) is applied to this system. This transforms the coupled equations into an equivalent triangular (sequential) system where the first equation contains only one unknown.
Solution: If the primary polynomial in the decoupled system is of degree four or lower, it can be solved analytically. Subsequent equations, typically of similar or lower degree, are then solved sequentially to obtain closed-form expressions for all coefficients.
Implementation: The symbolic computation software Maple was used to perform the Gröbner basis transformations. The results were validated against Finite Element Method (FEM) simulations (using ANSYS) and compared against exact linear theory solutions.

Case Studies
The methodology was demonstrated on three distinct boundary condition configurations:

Case 1: Fully clamped circular plate.
Case 2: Simply supported immovable circular plate.
Case 3: Simply supported immovable plate with an overhang (a more complex configuration involving piecewise displacement fields).

For each case, the authors utilized a "Current Study" (CS) formulation with a single unknown coefficient for transverse displacement and two for in-plane displacement. A "Two-Coefficient Study" (TCS) was also tested to evaluate accuracy limits; however, the TCS resulted in a primary Gröbner basis equation of ninth order, which precluded an analytical solution due to the Abel–Ruffini theorem.

Key Results

Closed-Form Solutions: The study successfully derived explicit symbolic closed-form solutions for transverse displacement, moments, shear forces, and stresses for all three cases. These solutions retain geometric ( $a, t$ ), material ( $E, \nu$ ), and loading ( $q$ ) parameters in explicit form.
Accuracy: The closed-form solutions (CS) showed excellent agreement with FEM results, with errors remaining below 4% even for large deflections (up to $1.5t$ ). In contrast, linear theory solutions deviated significantly, showing errors exceeding 80% at high displacement levels.
Parametric Insight: Because the solutions are symbolic, the influence of parameters (such as thickness-to-radius ratio) on displacement and stress can be examined directly without re-running numerical simulations.
Limitations: Increasing the complexity of the trial functions (TCS) improved accuracy slightly (approx. 2%) but rendered the system unsolvable analytically, highlighting a trade-off between the degree of the polynomial system and the feasibility of obtaining a closed-form solution.

Significance and Claims
The paper claims that the primary novelty lies in exploiting Gröbner basis techniques as an enabling tool to derive closed-form solutions from nonlinear mechanics formulations that are traditionally solved numerically. The authors assert that this approach offers a practical pathway for:

Deriving analytical solutions for classes of nonlinear structural mechanics problems.
Facilitating rapid parametric studies and design investigations by retaining explicit parameter dependence.
Providing a general framework applicable to any nonlinear mechanics formulation that produces polynomial algebraic systems through weighted-residual approximations.

The study concludes that while the method is limited to cases where the resulting primary polynomial is of degree four or less, it successfully demonstrates the potential of a new general framework for analytical solutions in nonlinear engineering mechanics, offering a significant accuracy improvement over linear theory while maintaining the efficiency of symbolic analysis.

Technical Summary: Application of Gröbner Basis Methodology to Geometrically Nonlinear Axisymmetric Circular Isotropic Plates

More like this