TY - JOUR

T1 - Three-dimensional isogeometric analysis of flexoelectricity with MATLAB implementation

AU - Ghasemi, Hamid

AU - Park, Harold S.

AU - Zhuang, Xiaoying

AU - Rabczuk, Timon

N1 - Funding Information: Acknowledgment: Hamid Ghasemi acknowledge the support of the Mechanical Engineering department at Arak University of Technology. Xiaoying Zhuang gratefully acknowledge the financial support by European Research Council for COTOFLEXI project (802205). Harold Park acknowledges the support of the Mechanical Engineering department at Boston University. Timon Rabczuk gratefully acknowledge financial support by the 2019 Foreign Experts Plan of Hebei Province.

PY - 2020/8/20

Y1 - 2020/8/20

N2 - 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.

AB - 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.

KW - B-spline elements

KW - Coupled electromechanical problem

KW - Finite elements

KW - Flexoelectricity

KW - Isogeometric Analysis (IGA)

KW - MATLAB

UR - http://www.scopus.com/inward/record.url?scp=85090947321&partnerID=8YFLogxK

U2 - 10.32604/cmc.2020.08358

DO - 10.32604/cmc.2020.08358

M3 - Article

AN - SCOPUS:85090947321

VL - 65

SP - 1157

EP - 1179

JO - Computers, Materials and Continua

JF - Computers, Materials and Continua

SN - 1546-2218

IS - 2

ER -