Data-adaptive spline surfaces for non-separable hyperelastic energy functions

Invariant-based models for incompressible isotropic hyperelasticity are typically formulated as functions of the first and second invariants, W=W(Ī1,Ī2). A widely used class of models employs separable representations of the form W(Ī1,Ī2)=W1(Ī1)+W2(Ī2), which enable efficient calibration and straightforward enforcement of modeling constraints. However, this decomposition implicitly restricts the coupling between the invariants and may limit the achievable accuracy for complex material responses. Fully coupled data-driven approaches overcome this limitation but often require nonlinear optimization and large parameter sets. In this contribution, we propose a compact alternative: a bivariate B-spline surface defined directly on the physically admissible invariant domain. By aligning the approximation space with physically realizable states, all model parameters contribute meaningfully to the constitutive response. We utilize homogeneous deformation modes to perform a calibration directly from analytical stress relations, eliminating the need for finite element model updating. Owing to the linear dependence of the spline representation on its coefficients, the resulting parameter identification problem reduces to a constrained linear least-squares problem. This enables fast, robust, and initialization-independent calibration, which makes parameter identification practically instantaneous. The results demonstrate that the proposed model improves accuracy compared to separable approaches while requiring only mild regularization in weakly sampled regions. The combination of computational efficiency and the linear structure of a highly expressive spline surface may make the approach attractive for applications involving repeated calibration procedures, such as uncertainty quantification, Bayesian inference, or interactive material characterization.