@article{ELLMER2026118349,

title = {Stretch-based hyperelastic electromechanical constitutive metamodels via gradient enhanced Gaussian predictors using hierarchical structure discovery},

author = {Nathan Ellmer and Rogelio Ortigosa and Jesús Martínez-Frutos and Roman Poya and Johann Sienz and Antonio J. Gil},

url = {https://www.sciencedirect.com/science/article/pii/S0045782525006218},

doi = {https://doi.org/10.1016/j.cma.2025.118349},

issn = {0045-7825},

year  = {2026},

date = {2026-01-01},

journal = {Computer Methods in Applied Mechanics and Engineering},

volume = {448},

pages = {118349},

abstract = {This paper introduces a new approach to developing electromechanical constitutive metamodels via the use of Gradient Enhanced Gaussian Predictors (Kriging). The formulation uses principal stretches for the isotropic mechanics, invariants for the electrostatics and coupling terms, and accounts for anisotropy through the relevant inclusion of anisotropic invariants associated with a respective symmetry integrity basis. Three novelties are presented in this paper. The first is the use of orthogonal projections to identify the most appropriate set of inputs - related to material anisotropy - for use in the metamodel. By projecting the stress and electric field data into several derivative bases - defined for each anisotropic class - and then reconstructing the quantities, the errors in reconstruction can be assessed thus inferring the most appropriate class of anisotropy. Furthermore, the procedure forms a pre-processing stage and is particularly useful when an underlying model is completely unknown as seen when modelling Relative Volume Elements. The second novelty arises from the use of a hybrid formulation, namely the principal stretches for isotropic mechanics and the electromechanical anisotropic invariants. This is beneficial during the projection procedure in reducing the cases where the projection matrix becomes singular but requires careful development of the correlation function to maintain physical symmetry conditions. Thirdly, the electromechanical metamodels are calibrated upon the concentric styled data before being integrated within a Finite Element framework and tested upon a range of challenging simulations including bending actuators with induced torsion, frilling due to bending with selected electrode placement, as well as buckling plates tested with three rank-one laminate materials with increasing levels of anisotropy due to physical contrasts. The successful calibration and implementation of the metamodels can be witnessed amongst the wide range of presented numerical examples.},

keywords = {Anisotropy, Constitutive modelling, electromechanics, Gradient Kriging, Principal stretches},

pubstate = {published},

tppubtype = {article}

}

@article{castanar2025stabilized,

title = {Stabilized Mixed Formulations for Incompressible Finite Strain Electromechanics Including Stress Accurate Analysis},

author = {Inocencio Castañar and Jesús Martínez-Frutos and Rogelio Ortigosa and Ramon Codina},

url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/nme.70089},

doi = {10.1002/nme.70089},

year  = {2025},

date = {2025-01-01},

urldate = {2025-01-01},

journal = {International Journal for Numerical Methods in Engineering},

volume = {126},

number = {22},

pages = {e70089},

keywords = {electromechanics, incompressible hyperelasticity, mixed formulations, orthogonal subgrid scales, PID2022-141957OA-C22, stabilization methods},

pubstate = {published},

tppubtype = {article}

}