@Article{cmc.2020.08358,
AUTHOR = {Hamid Ghasemi, Harold S. Park, Xiaoying Zhuang, Timon Rabczuk},
TITLE = {Three-Dimensional Isogeometric Analysis of Flexoelectricity with  MATLAB Implementation},
JOURNAL = {Computers, Materials \& Continua},
VOLUME = {65},
YEAR = {2020},
NUMBER = {2},
PAGES = {1157--1179},
URL = {http://www.techscience.com/cmc/v65n2/39867},
ISSN = {1546-2226},
ABSTRACT = {Flexoelectricity is a general electromechanical phenomenon where the electric 
polarization exhibits a linear dependency to the gradient of mechanical strain and vice versa. 
The truncated pyramid compression test is among the most common setups to estimate the 
flexoelectric effect. We present a three-dimensional isogeometric formulation of 
flexoelectricity with its MATLAB implementation for a truncated pyramid setup. Besides 
educational purposes, this paper presents a precise computational model to illustrate how the 
localization of strain gradients around pyramidal boundary shapes contributes in generation 
of electrical energy. The MATLAB code is supposed to help learners in the Isogeometric 
Analysis and Finite Elements Methods community to learn how to solve a fully coupled 
problem, which requires higher order approximations, numerically. The complete MATLAB 
code which is available as source code distributed under a BSD-style license, is provided in 
the part of Supplementary Materials of the paper.},
DOI = {10.32604/cmc.2020.08358}
}