“…This extra kinematic descriptor is constrained to be equal to the gradient of displacement field. We are aware of the fact that such a numerical code can be optimized for lower computational burdens, for instance, by using an isogeometric approach (as those presented in Fischer et al [74], Cuomo et al [75], Cazzani et al [76][77][78][79], Niiranen et al [80][81][82], Balobanov and Niiranen [83], Khakalo and Niiranen [84,85], Yaghoubi et al [86], Yildizdag et al [87,88], Capobianco et al [89], Greco et al [90], and Balobanov et al [91]).…”