One of the reasons for the success of the finite element method in Solid Mechanics, among other Applied Sciences, is its versatility to deal with bodies of arbitrary shape. In case the problem at hand is modeled by second‐order partial differential equations with Dirichlet conditions prescribed on a curvilinear boundary, method's isoparametric version for meshes consisting of curved triangles or tetrahedra has been mostly employed to recover the optimal approximation properties known to hold for polygonal or polyhedral domains and methods of order greater than one based on standard straight‐edged elements. However, besides algebraic and geometric inconveniences, the isoparametric technique is limited in scope, since its extension to degrees of freedom other than function values is not straightforward. In previous work the author introduced a simple alternative, which bypasses the above drawbacks without eroding qualitative approximation properties. Among other advantages, it can do without curved elements and is based only on polynomial Algebra. The purpose of this paper is two‐fold: On the one hand the mathematical foundations of this method are established, in the framework of a model Poisson problem solved with Lagrange elements; in particular optimal‐rate error estimates both in energy norm and in the mean‐square sense are demonstrated and a numerical comparison with the isoparametric technique is carried out. On the other hand this study is extended to Hermite elements with normal‐derivative degrees of freedom for solving the equation modeling the bending of clamped thin plates with a curvilinear edge.