ABSTRACT Flexoelectricity, the coupling between strain gradients and electric polarization, poses significant computational challenges due to its governing fourth‐order partial differential equations that require C 1 ‐continuous solutions. To address these issues, we propose a physics‐informed neural network (PINN) framework grounded in an energy‐based formulation that treats both forward and inverse problems within a unified architecture. The forward problem is recast as a saddle‐point optimization of the total potential energy, solved via the deep energy method (DEM), which circumvents the direct computation of high‐order derivatives. For the inverse problem of identifying unknown flexoelectric coefficients from sparse measurements, we introduce an additional variational loss that enforces stationarity with respect to the electric potential, ensuring robust and stable parameter inference. The framework integrates finite element‐based numerical quadrature for stable energy evaluation and employs hard constraints to rigorously enforce boundary conditions. Numerical results for both direct and converse flexoelectric effects show excellent agreement with mixed‐FEM solutions, and the inverse model accurately recovers material parameters from limited data. This study establishes a unified, mesh‐compatible, and scalable PINN approach for high‐order electromechanical problems, offering a promising alternative to traditional simulation techniques.
Moon et al. (2026) studied this question.