🔬 physics

Size-dependent elastostatic behavior of non-centrosymmetric solids via a mixed variational finite element method

This paper presents a novel mixed variational finite element method based on consistent couple stress theory to investigate the size-dependent elastostatic behavior of non-centrosymmetric solids, offering new insights into micro- and nano-scale responses and identifying potential experiments for measuring couple stress parameters.

Original authors: G. F. Dargush, B. T. Darrall, T. K. Scully, O. Lavan

Published 2026-08-10

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: G. F. Dargush, B. T. Darrall, T. K. Scully, O. Lavan

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

Imagine the world of materials science as a giant, invisible Lego set. For a long time, scientists built their models of how these blocks behave using "classical elasticity," a set of rules that treats every piece of matter like a smooth, uniform blob. In this old view, if you squish a block of clay, it just squishes; if you pull it, it stretches. The rules say that size doesn't matter—a tiny speck of clay should behave exactly like a giant boulder of the same stuff, just scaled down.

But there's a catch. Some materials, like the crystals inside your phone or certain types of sand, aren't perfectly symmetrical. They have a "handedness" or a specific direction, much like your left hand is different from your right. In the old rules, these asymmetrical materials were mostly famous for one trick: piezoelectricity. This is the magic where squeezing a crystal creates electricity (like in a lighter or a microphone). Scientists knew these materials were special, but they were missing a second, hidden trick. They suspected that when you squeeze these asymmetrical crystals, they might also twist or bend in ways that depend entirely on how small they are. It's like realizing that a tiny, asymmetrical Lego brick doesn't just squish; it might also curl up like a leaf when you push it, but only if the brick is microscopic. Understanding this "size-dependent" behavior is crucial because as we build smaller and smaller machines (nanotechnology), these tiny quirks could make or break our designs.

This paper dives into that missing trick using a new mathematical lens called "consistent couple stress theory." Think of this theory as upgrading the rules of the game. Instead of just tracking how much a material stretches (strain), the new rules also track how much it twists or curves (mean curvature). The authors, a team of engineers and mathematicians, realized that for asymmetrical crystals, there's a secret handshake between stretching and twisting. Just as squeezing a piezoelectric crystal creates electricity, squeezing an asymmetrical crystal might create a "couple stress"—a tiny internal torque that tries to rotate the material.

The team didn't just guess at this; they built a sophisticated computer simulation, a "mixed variational finite element method," to act as a virtual laboratory. Imagine this as a high-tech video game where they can build tiny beams, bars, and curved bridges out of different crystals (like Silicon Carbide and Quartz) and push, pull, and twist them to see what happens. They created a special set of math tools that allowed them to solve these problems without getting bogged down by the usual computational headaches, focusing specifically on how the material's internal "handedness" changes its stiffness.

What they found is fascinating. First, they confirmed that as these materials get smaller, they get stiffer, just like previous theories suggested. But the real surprise came from the asymmetry. In their simulations, when they pushed on a tiny bar made of hexagonal Silicon Carbide, it didn't just get shorter; it actually bent sideways! This is a purely linear effect, meaning it happens instantly under load without needing a critical "tipping point" like a buckling column. It's as if pushing down on a tiny, asymmetrical spring makes it curl to the left or right depending on the direction of the push.

The size of this effect depends on two invisible "length scales" within the material. One scale ( $l$ ) is like the material's natural thickness, and the other ( $l_p$ ) is a new, mysterious parameter specific to asymmetrical crystals. The simulations showed that as the ratio of these lengths changes, the material's stiffness can jump dramatically—sometimes becoming over 40 times stiffer than classical physics would predict. However, this isn't a universal rule; it depends heavily on the shape of the object and the direction of the force. For example, a deep beam being sheared (slid sideways) barely noticed these effects, while a thin, bending beam went crazy with stiffness.

The authors are careful to note that these are computer simulations, not physical measurements yet. They haven't measured the exact values of these new "couple stress" parameters in a real lab because the equipment to do so is incredibly difficult to build. Instead, they used their simulations to design a roadmap for future experiments. They suggest that if scientists build tiny cantilever beams (like a diving board) out of these crystals and push them, they should see this sideways bending and massive stiffening. If they do, it will prove that these asymmetrical crystals have a hidden, size-dependent personality that classical physics missed. The paper concludes that while the math is complex, the physical reality it predicts is a new frontier for understanding how the microscopic world twists and turns under pressure.

Technical Summary: Size-dependent elastostatic behavior of non-centrosymmetric solids via a mixed variational finite element method

Problem Statement
While classical elasticity and consistent couple stress theory (C-CST) have been extensively applied to centrosymmetric materials, the size-dependent elastostatic response of non-centrosymmetric solids remains largely unexplored within the C-CST framework. Non-centrosymmetric crystal structures, known for phenomena like piezoelectricity, possess a third-order constitutive tensor that couples force-stress with strain and couple-stress with mean curvature. This paper addresses the lack of computational tools and theoretical formulations necessary to analyze these materials at micro- and nano-scales, where size-dependent effects become significant. The primary challenge lies in developing a formulation that fully accounts for the thermodynamic and continuum mechanical interactions between force-stress, couple-stress, strain, and engineering mean curvature without requiring high-order continuity in the displacement field.

Methodology
The authors develop a novel mixed variational finite element method (FEM) based on consistent couple stress theory (C-CST). The methodology proceeds through the following steps:

Theoretical Formulation: The governing equations are derived from C-CST, assuming quasistatic conditions and infinitesimal deformations. The constitutive relations are established using a generalized Heckmann diagram approach, drawing an analogy to piezoelectricity. This introduces a third-order constitutive tensor ( $Q$ or $P$ ) that couples strain and mean curvature to force-stress and couple-stress, respectively. The formulation distinguishes between stiffness-based and compliance-based tensors, defining their relationships under different experimental constraints (e.g., constant strain vs. constant stress).
Mixed Variational Principle: To avoid the computational difficulty of $C^1$ continuity required by standard displacement-based formulations (which involve second derivatives of displacement), the authors derive a mixed variational principle. This principle treats displacements ( $u_i$ ) and couple-stresses ( $\mu_m$ ) as primary variables. By introducing Lagrange multipliers to satisfy kinematic constraints, the authors derive a "Principle of Stationary Mixed Potential Energy." This reduces the continuity requirement to $C^0$ for both displacement and couple-stress fields.
Finite Element Implementation: A two-dimensional plane strain discretization is implemented using low-order triangular elements.
- Displacements: Standard three-node linear shape functions are used.
- Couple-Stresses: To represent the couple-stress vector, the authors employ edge-based shape functions. Unlike centrosymmetric cases where the couple-stress vector is divergence-free, non-centrosymmetric materials generally require both tangential and normal components along element edges. This results in two degrees of freedom per edge.
- Numerical Solution: The weak form is discretized into a system of linear algebraic equations involving symmetric and skew-symmetric stress components. The system is solved using sparse matrix techniques in MATLAB.

Key Contributions

First Application to Non-Centrosymmetric Solids: This work represents the first computational exploration of size-dependent elastostatic behavior in non-centrosymmetric solids using C-CST.
Novel Mixed Variational Framework: The paper extends the mixed variational method (previously used for centrosymmetric materials) to the bi-anisotropic (non-centrosymmetric) case, successfully handling the full coupling between force/couple-stresses and strain/curvature.
Constitutive Modeling: The authors explicitly define the thermodynamic relationships and constitutive tensors for non-centrosymmetric materials, clarifying how material parameters depend on the variables held constant during testing.
Element Development: The formulation introduces a specific triangular edge element capable of capturing both tangential and normal components of couple-stress, essential for general non-centrosymmetric symmetries.

Results
The methodology was validated through four computational examples involving three non-centrosymmetric materials: Cubic Silicon Carbide (SiC), Hexagonal SiC, and Trigonal Quartz. The problems included a prismatic cantilever under bending, a prismatic bar under axial load, a deep beam under shear, and a curved beam under combined loading.

Convergence: The mixed FEM solutions demonstrated excellent convergence. The relative energy norm converged with a rate slightly less than unity for bending problems and nearly unity for axial loading, consistent with the use of linear shape functions.
Size-Dependent Stiffness: As observed in centrosymmetric materials, structures became stiffer as the non-dimensional intrinsic length scale ( $l/h$ ) increased, eventually reaching a saturation regime. Bending-dominated problems (e.g., cantilevers) showed significant stiffness amplification (up to ~43x), whereas shear-dominated problems (deep beams) showed minimal amplification (<2x).
Non-Centrosymmetric Effects:
- Coupling Parameter ( $l_p$ ): The second length scale parameter ( $l_p$ ), associated with the non-centrosymmetric coupling tensor, was found to further increase stiffness, particularly as $l_p$ approached $2l$ .
- Anisotropy: The response was highly dependent on the crystal orientation and loading plane. For example, in Hexagonal SiC, the coupling between axial load and transverse curvature was significant in the $x-y$ plane but negligible in the $x-z$ plane.
- Linear Transverse Deflection: A critical finding was that under axial compressive loading, non-centrosymmetric materials exhibit a linear transverse deflection (bending) without a critical buckling load. This is a purely linear mechanical effect analogous to the converse piezoelectric effect, where the sign of the deflection depends on the sign of the coupling parameter.
- Stress Asymmetry: The skew-symmetric part of the force-stress tensor was non-zero, leading to unequal shear stress components ( $\sigma_{xy} \neq \sigma_{yx}$ ) and distinct stress distributions compared to classical elasticity.

Significance and Claims
The paper claims to provide the first computational framework for analyzing size-dependent behavior in non-centrosymmetric elastic solids. The authors emphasize that the observed linear coupling between axial loading and transverse curvature offers a potential mechanism for physical experiments to measure the elusive couple-stress constitutive parameters ( $l_p$ and $P$ tensor components). The study suggests that shallow cantilever bending tests are effective for estimating the intrinsic length scale $l$ , while axial loading tests could isolate the non-centrosymmetric coupling parameters. The authors conclude that while the current work focuses on linear elastostatics, the developed framework lays the groundwork for future investigations into vibration, transient dynamics, three-dimensional problems, and finite deformations, including potential large-deformation induced symmetry breaking in centrosymmetric crystals.

More like this