Next Article in Journal
Architectural Graphics Research: Topics and Trends through Cluster and Map Network Analyses
Previous Article in Journal
General Response Formula for CFD Pseudo-Fractional 2D Continuous Linear Systems Described by the Roesser Model
Previous Article in Special Issue
Dynamics of Nonlocal Rod by Means of Fractional Laplacian
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Axial Porosities on Flexomagnetic Response of In-Plane Compressed Piezomagnetic Nanobeams

1
Department of Mechanics of Materials and Structures, Faculty of Civil and Environmental Engineering, Gdansk University of Technology, 80-233 Gdansk, Poland
2
Laboratory of Mechanics of Biomaterials, Research and Education Center “Materials”, Don State Technical University, Gagarina sq., 1, 344000 Rostov on Don, Russia
3
Department of Civil and Environmental Engineering and Architecture (DICAAR), University of Cagliari, Via Marengo, 2, 09123 Cagliari, Italy
4
Faculty of Mechanical Engineering, Bialystok University of Technology, 15-351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Submission received: 29 October 2020 / Revised: 16 November 2020 / Accepted: 23 November 2020 / Published: 24 November 2020
(This article belongs to the Special Issue Time and Space Nonlocal Operators in Structural Mechanics)

Abstract

:
We investigated the stability of an axially loaded Euler–Bernoulli porous nanobeam considering the flexomagnetic material properties. The flexomagneticity relates to the magnetization with strain gradients. Here we assume both piezomagnetic and flexomagnetic phenomena are coupled simultaneously with elastic relations in an inverse magnetization. Similar to flexoelectricity, the flexomagneticity is a size-dependent property. Therefore, its effect is more pronounced at small scales. We merge the stability equation with a nonlocal model of the strain gradient elasticity. The Navier sinusoidal transverse deflection is employed to attain the critical buckling load. Furthermore, different types of axial symmetric and asymmetric porosity distributions are studied. It was revealed that regardless of the high magnetic field, one can realize the flexomagnetic effect at a small scale. We demonstrate as well that for the larger thicknesses a difference between responses of piezomagnetic and piezo-flexomagnetic nanobeams would not be significant.

1. Introduction

Flexomagneticity arises through elastic strain gradient or magnetic field gradient during electric magnetization in the magneto-elastic coupling in smart structures and actuators [1,2,3]. Such an effect should be significant in nano electro-mechanical systems (NEMS) and other smart sensors and actuators. Similar to this influence, viz. flexoelectricity for centrosymmetric and non-centrosymmetric structures has been estimated widely [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. However, the flexomagneticity effect has been less known. Although the structure and physics of a flexomagneticity phenomenon are very complicated, it has economic implications. Physically, it may be difficult to interpret, but the basic idea satisfactorily shows the importance of flexomagneticity.
In discussing a mechanical response of nanomaterials with magneto-mechanical coupling, the importance of piezomagneticity has been profoundly evaluated by scholars in the contemporary decade [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39]. However, fewer studies are available regarding flexomagneticity [40,41,42,43,44,45]. In the available literature, Sidhardh and Ray [40] worked on the bending response of a thin cantilever nanobeam with flexomagnetic property. They discussed both direct and reverse impacts of magneto-elastic coupling with the presence of the surface elasticity. Zhang et al. [41] conducted the importance of flexomagneticity for a nanoscale size-dependent Euler-Bernoulli beam exposed to transverse static loading. To present the size-dependent mechanical behavior of the structure, they utilized the surface elasticity. Both direct and converse flexomagnetic influences were investigated when the nano-sized beam was kept in ends with fixed, pivot, and free edge conditions. One of their momentous achievements was the conclusion that the flexomagneticity is a scale-dependent property of materials. Recently, Malikan and Eremeyev [42] investigated vibrating nanobeams by taking into account the piezomagnetic, particularly flexomagnetic, properties. The free vibration frequencies were evaluated by linear assumptions of strain and the impact of size-dependency was inspected on the basis of a new nonlocal elasticity theory. The size-dependent behavior of the flexomagnetic effect was affirmed by their results. In another effort, Malikan and Eremeyev [43] explored the flexomagnetic response of a smart nanobeam in a vibrational condition based on large frequency modes. The small scale response of the nanoscale beam was searched by imposing a nonlocal strain gradient elasticity approach into the constitutive equations. Besides these, Malikan and Eremeyev [44] performed research on non-linear static bending of smart nanoscale beams while the material included a remarkable flexomagnetic response. The computational model was solved by coupling between an analytical and numerical solution method. The new finding demonstrated that the presence of the flexomagnetic feature leads to diminishing the deflections. More recently, Malikan et al. [45] examined post-buckling stability of a nanoparticle in which both flexomagnetic and piezomagnetic properties were included. The calculations extracted new findings that are helpful for magnetic nanosensors applications.
To the best of our knowledge, the flexomagnetic studies on the mechanical response of nanostructures referred to above represent the majority of cases in the literature, unless otherwise stated. In these references, the models are restricted to piezomagnetic nanobeams under some mechanical analyses. The present research attempts to demonstrate the flexomagnetic property for the stability problem of a nano-sized beam, while it includes a material imperfection with intentional nonlocality and size-dependent characteristics according to the nonlocal strain gradient constitutive equation. The structural and material imperfection is estimated in the framework of different types of porosities. The nonlocal influences that can affect the flexomagnetic response of the Euler–Bernoulli nanoscale beam are addressed in this paper. A sinusoidal transverse deflection is assumed when applying the Navier approach to beam flexure. Terminally, the nanobeam is considered for variations of key parameters based on the three cases, i.e., a simple nanobeam, piezomagnetic nanobeam, and piezo-flexomagnetic nanobeam.

2. Formulation of the Problem

2.1. Constitutive Relations for Piezo-Flexomagnetic Solids

Pursuing [1,2,3], we briefly introduce constitutive relations for flexomagnetic material. In the following, we assume infinitesimal deformations under isothermal conditions. Therefore, the variables are the displacements u and the magnetic field is H as below
u = u ( x ) ,   H = H ( x ) .
where x is a position vector, and H is a tensor of first-order. Within the flexomagneticity, the free energy density function U has the form
U = U ( ε , η , H ) = 1 2 H a H + 1 2 ε : C : ε + 1 2 η g η + ε : r η . H q : ε H f η
where ε is the strain tensor and its gradient is
ε = 1 2 ( u + u T ) ,   η = ε .
where is the 3D nabla-operator in a general case. We introduced to Equation (2) several tensors of material parameters. q is the third-order piezomagnetic tensor, a is the second-order magnetic permeability tensor, g is the sixth-order gradient elasticity tensor, C is the fourth-order elasticity coefficient tensor, f is the fourth-order flexomagnetic tensor, and the strain and strain-gradient tensors are coupled by r, which itself is a fifth-order tensor. Moreover, “∙”, “:”, and “⋮” stand for scalar (inner) products in spaces of vectors, second-order and third-order tensors, respectively.
In the following, we express H through the magnetic potential ψ [40,41],
H = ψ .
The governing equations of the static flexomagneticity can be derived using the virtual work principle as
δ Π = δ A .
in which Π = V U d V , where V is the volume of domain that occupies the flexomagnetic solid and δ A is the work of external loads.
For simplicity, we assume the simple form of
δ A = V F δ u + V t δ u d s .
where F and t are external mass forces and surface traction.
Using the standard calculus of variations from Equation (5) we get
( σ ξ ) + F = 0 ,
B = 0 .
in which B is the magnetic induction vector. The following constitutive equations are introduced
σ = U ε C : ε + r η H q ,
ξ = U η g : η + ε : r H f ,
B = U H = a H + q : ε + f η .

2.2. The Piezo-Flexomagnetic Beam Model

The piezo-flexomagnetic nanobeam studied in this research is exhibited in Figure 1. In the figure, length and height of the beam are denoted by L and h, respectively.
The constitutive equations (Equation (8)) can be re-written for a beam incorporating the flexomagnetic property as [40,41]
σ x x = C 11 p ε x x q 31 H z ,
ξ x x z = g 31 η x x z f 31 H z ,
B z = a 33 H z + q 31 ε x x + f 31 η x x z .
in which η x x z and ε x x are the gradient of the axial elastic strain and the strain itself, C 11 = C 1111 is the elastic modulus, σ x x is the axial stress, f 31 = f 3311 denotes the component of the fourth-order flexomagnetic coefficients tensor, a 33 represents the component of the second-order magnetic permeability tensor, q 31 = q 311 depicts the component of the third-order piezomagnetic tensor, ξ x x z is the component of the higher-order hyper stress tensor and is an induction of coverse flexomagnetic effect, B z and H z exhibit the magnetic flux and the component of magnetic field, respectively, and g 31 = g 311311 illustrates the influence of the sixth-order gradient elasticity tensor. It is worth mentioning that the piezomagnetic tensor would be non-zero for non-centrosymmetric ferroics only, but the flexomagnetic tensor would be non-zero for both centrosymmetric and non-centrosymmetric materials.
The displacement field with respect to the Euler–Bernoulli beam is available as [46,47,48]
u 1 ( x , z ) = u ( x ) z d w ( x ) d x , u 3 ( x , z ) = w ( x ) .
where ui (i = 1,3) represent the points’ displacements in direction of x and z, u and w are the axial and transverse displacements of the mid-plan, respectively, see Figure 1. To show the thickness coordinate, the z parameter is used.
Due to the linear problem studied in this research, the linear Lagrangian strain can be employed as
ε i j = 1 2 ( u i x j + u j x i ) .
The components of the transverse strain and the strain gradient can be presented on the basis of substituting Equation (12) into Equation (13) as below
ε x x = d u d x z d 2 w d x 2 , η x x z = d ε x x d z = d 2 w d x 2 , η x x x = d ε x x d x = z d 3 w d x 3 .
As the η x x x is small compared to the η x x z , it can be ignored. With respect to the Lagrange’s principle, we have
δ ( Π W + Π U ) = 0 .
in which Π W and Π U depict the performed work by outer loads, and the total internal strain energy (magnetic potential energy and mechanical strain energy). The total strain energy by means of Equation (14) can be demonstrated as
δ Π U = V ( σ x x δ ε x x + ξ x x z δ η x x z B z δ H z ) d V .
One can obtain the governing equation and non-classical boundary conditions as below
δ Π U 1 M e c h = 0 L ( d N x d x δ u + d 2 M x d x 2 δ w + d 2 T x x z d x 2 δ w ) d x ,
δ Π U 2 M a g = 0 L h / 2 h / 2 d B z d z δ Ψ d z d x ,
δ Π U 1 M e c h = ( N x δ u M x d δ w d x T x x z d δ w d x + d M x d x δ w + d T x x z d x δ w ) | 0 L ,
δ Π U 2 M a g = 0 L ( B z δ Ψ ) | h / 2 h / 2 d x .
in which
N x = h / 2 h / 2 σ x x d z ,
M x = h / 2 h / 2 σ x x z d z ,
T x x z = h / 2 h / 2 ξ x x z d z .
The performed work by external factors can be expressed as below [49,50]
Π W = 1 2 0 L N x 0 ( d w d x ) 2 d x .
Its first variational form is
δ Π W = 0 L N x 0 ( d δ w d x d w d x ) d x
in which N x 0 shows the axial membrane load.
Hereafter, the magnetic field’s transverse component can be expressed as
H z + d Ψ d z = 0 .
Assuming the condition of a closed circuit as well as the inverse piezomagnetic effect, one gives the magnetic boundary conditions as
Ψ ( + h 2 ) = ψ ,
Ψ ( h 2 ) = 0 .
where ψ determines the external magnetic potential applied to the upper surface of the beam. With mixing Equations (11), (17b), (18b), (24) and (25), one can derive the magnetic polarization and magnetic field as [40,41]
Ψ = q 31 2 a 33 ( z 2 h 2 4 ) d 2 w d x 2 + ψ h ( z + h 2 ) ,
H z = z q 31 a 33 d 2 w d x 2 ψ h .
The nanoscale atomic interactions can be projected in a continuum space by means of the nonlocal strain gradient elasticity theory (NSGT) given as [51]
( 1 μ d 2 d x 2 ) σ x x N L = ( 1 l 2 d 2 d x 2 ) σ x x L .
In order to allocate the influence of nonlocality, namely stiffness-softening, the nonlocal parameter µ (unit: square nanometers = ( n m ) 2 ) is employed. Note that μ ( n m ) 2 = ( e 0 a ) 2 , in which e 0 and a are two small scale factors that determine the nonlocal parameter. Furthermore, in order to set the effect of the size deduction, namely stiffness-hardening, the l ( n m ) parameter is used, which is named as the strain gradient length scale parameter. It is to be noted that [52,53] confirmed the fact that the aforesaid small scale factors are not material constants and can be varied by respecting a variety of conditions. Additionally, the NL and L indexes respectively express the nonlocal and the local components of stress.
Thus, inserting Equations (26)–(28) into Equations (9)–(11), it is possible to present the stress field component, the higher-order moment stress tensor component, and magnetic induction component in the NSGT form as
( 1 μ d 2 d x 2 ) σ x x = ( 1 l 2 d 2 d x 2 ) [ C 11 p d u d x z ( C 11 p + q 31 2 a 33 ) d 2 w d x 2 + q 31 ψ h ] ,
( 1 μ d 2 d x 2 ) ξ x x z = ( 1 l 2 d 2 d x 2 ) [ ( g 31 + q 31 f 31 z a 33 ) d 2 w d x 2 + f 31 ψ h ] ,
( 1 μ d 2 d x 2 ) B z = ( 1 l 2 d 2 d x 2 ) ( f 31 d 2 w d x 2 a 33 ψ h ) .
Hence, on the basis of Equations (29)–(31), Equations (19)–(21) can be expanded as [54,55,56,57,58,59,60,61]
( 1 μ d 2 d x 2 ) N x = ( 1 l 2 d 2 d x 2 ) ( C 11 p h d u d x + q 31 ψ ) ,
( 1 μ d 2 d x 2 ) M x = ( 1 l 2 d 2 d x 2 ) ( I z ( C 11 p + q 31 2 a 33 ) d 2 w d x 2 ) ,
( 1 μ d 2 d x 2 ) T x x z = ( 1 l 2 d 2 d x 2 ) ( g 31 h d 2 w d x 2 + f 31 ψ ) .
in which I z = A z 2 d A is the moment of inertia.
Due to inevitable variations in the manufacturing processes, the presence of some porosity in nanobeams is unavoidable. Inclusion of this imperfection into mechanical analysis of the piezomagnetic-flexomagnetic nano-sized beam is performed as [62]
C 11 p = C 11 ( 1.21 α λ ( x ) 1.21 ) 2.3
where α denotes the porosity coefficient. Axial porosities are defined mathematically and analytically in Table 1 [62].
On the basis of Equations (17a) and (23), and replacing into Equation (15), one can derive the local governing relations as below
d N x d x = 0 ,
d 2 M x d x 2 + d 2 T x x z d x 2 + N x 0 d 2 w d x 2 = 0 .
Equations (36) and (37) are decoupled and thus, Equation (37) gives the stability equation in order to have values of critical buckling loads.
Thus, we should transfer the local stability equation (Equation (37)) to a size-dependent relation. In so doing, inserting Equation (37) into Equation (33), one gets
M x = μ ( d 2 T x x z d x 2 + N x 0 d 2 w d x 2 ) I z ( C 11 p + q 31 2 a 33 ) ( 1 l 2 d 2 d x 2 ) d 2 w d x 2 .
Then, Equation (37) can be re-derived by mixing Equation (38) and Equation (34) as follows
( 1 μ d 2 d x 2 ) ( B d 4 w d x 4 + N x 0 d 2 w d x 2 ) + D ( 1 l 2 d 2 d x 2 ) d 4 w d x 4 = 0 .
in which B = g 31 h , and D = I z ( C 11 p + q 31 2 a 33 ) .
Here, it is required to define the pre-buckling compression as membrane loads as
N x 0 = N M e c h + N M a g .
where N M e c h and N M a g are the axial membrane magnetic and mechanical loads assumed as follows:
N M e c h = P c r ,
N M a g = q 31 ψ .

3. Solution of the Problem

Regarding the analytical closed-form solution, we apply the following transverse deflection equation
w ( x ) = m = 1 X m ( x ) .
in which X m is an allowable function satisfying the boundary conditions of simply-simply supported (S-S) for two ends.
The dedicated kinematic and nonlocal strain gradient constitutive boundary conditions for the S-S nanobeam are expressed by Table 2 [63,64,65,66]
The pointed conditions in Table 2 may be satisfied by the admissible function given below [57]
X m ( x ) = sin ( α m x ) .
in which α m = m π L .
Substituting Equation (42) into Equation (39), the closed-form size-dependent buckling relation for the piezo-flexomagnetic nanobeam becomes
P c r = 0 L ( B μ + l 2 D ) α m 6 + ( B q 31 ψ μ + D ) α m 4 + q 31 ψ α m 2 μ α m 4 α m 2 Y m d x .
in which Y m is a residue. It is important to bear in mind that all of the results of the present work are given for m = 1.

4. Numerical Results

4.1. Validation of Results

The literature survey clearly and obviously showed that the static bifurcation buckling of a piezo-flexomagnetic nanobeam has not been studied thus far. Therefore, the only path to validate correctness of the formulation is neglecting piezomagneticity, flexomagneticity, and the strain gradient model to compare the results with a nano-sized circular beam [67,68]. Both parts of the literature employed a classical beam, however, ref. [67] used the differential transform solution method and ref. [68] used an explicit solution method. As can be shown from Table 3, our results are entirely matched with the literature.

4.2. Stability Analysis

Investigating the flexomagneticity effect on the stability of a nanostructure is the main goal of this work. Table 4 presents applied material properties [40,41]. To take a rational amount for the nonlocal parameter, 0.5 nm < e0a < 0.8 nm [69], and 0 < e0a ≤ 2 nm [70,71], are utilized.
Given Figure 2a, the nonlocal coefficient variations are plotted for four nanobeam states. That is, first, we have just the usual nanobeam by eliminating magnetic effects; the latter is a state where we merely have the piezomagnetic effect and finally the third and fourth states are when we have both the piezo- and flexomagnetic effects but with two different values of flexomagnetic property. As is clear from the figure, increasing the numerical value of the nonlocal parameter reduces the critical loads in all four states. Thereby, one of the important results of this graph is that when we consider a positive magnetic field and we have the piezomagnetic as well as flexomagnetic effect, the nanobeam has greater stability against the axial membrane forces. Furthermore, while purely having the piezomagnetic effect, the least stability is observed for the magneto-mechanical nanobeams. Further investigation is needed and is shown by the next figures. On the other hand, Figure 2b shows a comparison of both piezomagnetic and piezo-flexomagnetic nanoscale beams while the strain gradient parameter is the changeable factor of the horizontal axis of the figure. As it is clarified by Figure 2b, it is noteworthy to say that increasing the strain gradient parameter results in a decrease in difference between the results of the piezo-flexomagnetic nanobeam with the piezomagnetic nanobeam. It is worth underlining that this proximity pertains to the flexomagnetic feature. This means that the increase of the value of the length scale parameter leads to a stiffening effect; therefore, by increasing the values of this parameter, the critical load’s results tend to each other in magnetic beams. It can be argued that the length scale parameter makes the flexomagnetic effect ineffective. Eventually, if the length scale parameter is a large numerical value, it can be stated that the flexomagnetic effect is nothing and all the magnetic nanobeams will respond similarly.
In Figure 3, we investigate the effects of given patterns of the porosities for two cases of the nano-sized magnetic beam, the first one with the flexomagnetic property (PFM) and the later one ignoring this physical feature. We observed in the previous figure that in attending the positive magnetic field, if the piezomagnetic nanobeam has a flexomagnetic property (piezo-flexomagnetic), it is the most stable case in terms of magneto-mechanical nanobeams. It can now be seen from this figure that for all PFM cases, the further in-plane resistance can be observed. Hence, one can prove that this effect makes material stiffer. Furthermore, increase of the value of the porosity parameter (α) leads to reduction of material stability in all cases. However, this decline is more noticeable for X ¯ and O ¯ porosities and insignificant for O and X samples. It is also notable that in terms of X ¯ and O ¯ porosity types, an increase of value of the porosity parameter leads to more gaps between PFM and PM. This means that these kinds of porosities make the flexomagnetic property more important.
Figure 4 is presented correspondingly to reveal the incremental variations of the external magnetic potential. In this figure, as in the previous ones, the nanobeam is investigated in different cases. As can be seen, the increased magnetization potential leads to greater stability of the magnetic nanobeams. Of importance in the diagram is the magnetic nanobeams containing lesser critical loads than the conventional nanobeam. As a matter of fact, the in-plane static stability for such beams in negative amounts of magnetic potential is smaller than that of the conventional nanobeam. This downward/upward trend is linear in the form of a steep slope, indicating that in very strong positive magnetic fields, the piezomagnetic effects and, especially the flexomagnetic ones, will be more and more pronounced. However, these effects are also significant and undeniable in a weak magnetic environment.
Figure 5 displays the aforementioned states of nanobeams concerning changes in their thicknesses. It can be inferred from the figure that at very small thicknesses, and also positive magnetic potential, the flexomagnetic property plays a vital role in in-plane stability. With the thickening of the nanobeams, the difference between the results of the piezomagnetic and piezo-flexomagnetic nanobeams decreases. Thereupon, for larger thicknesses, the flexomagnetic effect seems to be unimportant.

5. Conclusions

This study presented the stability capacity of a porous nanobeam involving piezomagnetic as well as flexomagnetic impacts. To date, it is known that the mechanism of action of nanostructures is based on two principles of hardening and softening. This research applied these actions to a piezo-flexomagnetic nanobeam. Substituting Lagrangian, and nonlocal theory of strain gradient elasticity, the stability relation of the piezo-flexomagnetic nanobeam was gained. Thereafter, the Navier method gave a closed-form solution to reach numerical amounts of the in-plane static stability.
The flexomagnetic effect as a complex physical phenomenon into the magneto-mechanical coupling is known as a size-dependent property that was also here affected by small-scale parameters. Furthermore, the variations of thickness of the nanobeam affected the flexomagneticity, and this property is further noticeable for lower thicknesses of nanobeams. Moreover, comparing a piezomagnetic nanobeam with a piezo-flexomagnetic one showed that this feature presents more stable material. In addition, based on our observations, it was proved that despite the flexomagnetic effect being important in the lowest positive external magnetic potentials, if the potential is sufficiently large, the capacity of static stability for piezo-flexomagnetic nanobeams will be enhanced markedly. Another significant point obtained in this research work was the influence of porosity on the flexomagnetic response of the piezomagnetic nano-sized beam. The results showed that in some patterns of porosity, this imperfection can affect the flexomagnetic behavior of the material.

Author Contributions

Conceptualization, M.M., V.A.E. and K.K.Ż.; methodology, M.M., V.A.E. and K.K.Ż.; software, M.M.; validation, M.M.; formal analysis, M.M.; investigation, M.M. and V.A.E.; resources, M.M., V.A.E. and K.K.Ż.; data curation, M.M., V.A.E. and K.K.Ż.; writing—original draft preparation, M.M.; writing—review and editing, V.A.E. and K.K.Ż.; visualization, M.M.; supervision, V.A.E.; project administration, V.A.E.; funding acquisition, V.A.E. and K.K.Ż. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

V.A.E. acknowledges the support of the Government of the Russian Federation (contract No. 14.Z50.31.0046). K.K.Ż. acknowledges the support of the Ministry of Science and Higher Education of Poland (project No. W/WM-IIM/3/2020).

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Kabychenkov, A.F.; Lisovskii, F.V. Flexomagnetic and Flexoantiferromagnetic Effects in Centrosymmetric Antiferromagnetic Materials. Tech. Phys. 2019, 64, 980–983. [Google Scholar] [CrossRef]
  2. Eliseev, E.A.; Morozovska, A.N.; Glinchuk, M.D.; Blinc, R. Spontaneous flexoelectric/flexomagnetic effect in nanoferroics. Phys. Rev. B 2009, 79, 165433. [Google Scholar] [CrossRef]
  3. Lukashev, P.; Sabirianov, R.F. Flexomagnetic effect in frustrated triangular magnetic structures. Phys. Rev. B 2010, 82, 094417. [Google Scholar] [CrossRef]
  4. Ma, W. Flexoelectricity: Strain gradient effects in ferroelectrics. Phys. Scr. 2007, T129, 180–183. [Google Scholar] [CrossRef]
  5. Lee, D.; Yoon, A.; Jang, S.Y.; Yoon, J.-G.; Chung, J.-S.; Kim, M.; Scott, J.F.; Noh, T.W. Giant Flexoelectric Effect in Ferroelectric Epitaxial Thin Films. Phys. Rev. Lett. 2011, 107, 057602. [Google Scholar] [CrossRef]
  6. Nguyen, T.D.; Mao, S.; Yeh, Y.-W.; Purohit, P.K.; McAlpine, M.C. Nanoscale Flexoelectricity. Adv. Mater. 2013, 25, 946–974. [Google Scholar] [CrossRef]
  7. Zubko, P.; Catalan, G.; Tagantsev, A.K. Flexoelectric Effect in Solids. Annu. Rev. Mater. Res. 2013, 43, 387–421. [Google Scholar] [CrossRef] [Green Version]
  8. Yudin, P.V.; Tagantsev, A.K. Fundamentals of flexoelectricity in solids. Nanotechnology 2013, 24, 432001. [Google Scholar] [CrossRef]
  9. Jiang, X.; Huang, W.; Zhang, S.; Tagantsev, A.K. Flexoelectric nano-generators: Materials, structures and devices. Nano Energy 2013, 2, 1079–1092. [Google Scholar] [CrossRef]
  10. Yurkov, A.S.; Tagantsev, A.K. Strong surface effect on direct bulk flexoelectric response in solids. Appl. Phys. Lett. 2016, 108, 022904. [Google Scholar] [CrossRef]
  11. Wang, B.; Gu, Y.; Zhang, S.; Chen, L.-Q. Flexoelectricity in solids: Progress, challenges, and perspectives. Prog. Mater. Sci. 2019, 106, 100570. [Google Scholar] [CrossRef]
  12. Cross, L.E. Flexoelectric effects: Charge separation in insulating solids subjected to elastic strain gradients. J. Mater. Sci. 2006, 41, 53–63. [Google Scholar] [CrossRef]
  13. Ma, W.; Cross, L.E. Observation of the flexoelectric effect in relaxor Pb(Mg1/3Nb2/3)O3 ceramics. Appl. Phys. Lett. 2001, 78, 2920–2921. [Google Scholar] [CrossRef]
  14. Ma, W.; Cross, L.E. Flexoelectricity of barium titanate. Appl. Phys. Lett. 2006, 88, 232902. [Google Scholar] [CrossRef]
  15. Zubko, P.; Catalan, G.; Buckley, A.; Welche, P.R.L.; Scott, J.F. Strain-gradient-induced polarization in SrTiO3 single crystals. Phys. Rev. Lett. 2007, 99, 167601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Eremeyev, V.A.; Ganghoffer, J.-F.; Konopińska-Zmysłowska, V.; Uglov, N.S. Flexoelectricity and apparent piezoelectricity of a pantographic micro-bar. Int. J. Eng. Sci. 2020, 149, 103213. [Google Scholar] [CrossRef]
  17. Samani, M.S.E.; Beni, Y.T. Size dependent thermo-mechanical buckling of the flexoelectric nanobeam. Mater. Res. Express 2018, 5, 085018. [Google Scholar] [CrossRef]
  18. Malikan, M.; Eremeyev, V.A. On the Dynamics of a Visco–Piezo–Flexoelectric Nanobeam. Symmetry 2020, 12, 643. [Google Scholar] [CrossRef] [Green Version]
  19. Singhal, A.; Sedighi, H.M.; Ebrahimi, F.; Kuznetsova, I. Comparative study of the flexoelectricity effect with a highly/weakly interface in distinct piezoelectric materials (PZT-2, PZT-4, PZT-5H, LiNbO3, BaTiO3). Waves Random Complex Media 2019, 1–19. [Google Scholar] [CrossRef]
  20. Chen, W.; Yan, Z.; Wang, L. On mechanics of functionally graded hard-magnetic soft beams. Int. J. Eng. Sci. 2020, 157, 103391. [Google Scholar] [CrossRef]
  21. Ke, L.-L.; Wang, Y.-S. Free vibration of size-dependent magneto-electro-elastic nanobeams based on the nonlocal theory. Phys. E Low-Dimens. Syst. Nanostruct. 2014, 63, 52–61. [Google Scholar] [CrossRef]
  22. Ping, T. Nonlinear free vibration analysis of nanobeams under magnetic field based on nonlocal elasticity theory. J. Vibroeng. 2016, 18, 1912–1919. [Google Scholar] [CrossRef]
  23. Baghani, M.; Mohammadi, M.; Farajpour, A. Dynamic and Stability Analysis of the Rotating Nanobeam in a Nonuniform Magnetic Field Considering the Surface Energy. Int. J. Appl. Mech. 2016, 8, 1650048. [Google Scholar] [CrossRef]
  24. Ebrahimi, F.; Barati, M.R. Porosity-dependent vibration analysis of piezo-magnetically actuated heterogeneous nanobeams. Mech. Syst. Signal Process. 2017, 93, 445–459. [Google Scholar] [CrossRef]
  25. Zenkour, A.M.; Arefi, M.; AlShehri, N.A. Size-dependent analysis of a sandwich curved nanobeam integrated with piezomagnetic face-sheets. Results Phys. 2017, 7, 2172–2182. [Google Scholar] [CrossRef]
  26. Sun, X.-P.; Hong, Y.-Z.; Dai, H.; Wang, L. Nonlinear frequency analysis of buckled nanobeams in the presence of longitudinal magnetic field. Acta Mech. Solida Sin. 2017, 30, 465–473. [Google Scholar] [CrossRef]
  27. Arefi, M.; Zenkour, A.M. Transient sinusoidal shear deformation formulation of a size-dependent three-layer piezo-magnetic curved nanobeam. Acta Mech. 2017, 228, 3657–3674. [Google Scholar] [CrossRef]
  28. Liu, H.; Yang, J. Vibration of FG magneto-electro-viscoelastic porous nanobeams on visco-Pasternak foundation. Compos. Part B Eng. 2018, 155, 244–256. [Google Scholar] [CrossRef]
  29. Karami, B.; Shahsavari, D.; Li, L. Hygrothermal wave propagation in viscoelastic graphene under in-plane magnetic field based on nonlocal strain gradient theory. Phys. E Low-Dimens. Syst. Nanostruct. 2018, 97, 317–327. [Google Scholar] [CrossRef]
  30. Arefi, M.; Arani, A.H.S. Higher order shear deformation bending results of a magnetoelectrothermoelastic functionally graded nanobeam in thermal, mechanical, electrical, and magnetic environments. Mech. Based Des. Struct. Mach. 2018, 46, 669–692. [Google Scholar] [CrossRef]
  31. Alibeigi, B.; Beni, Y.T. On the size-dependent magneto/electromechanical buckling of nanobeams. Eur. Phys. J. Plus 2018, 133, 398. [Google Scholar] [CrossRef]
  32. Azrar, A.; Ben Said, M.; Azrar, L.; Aljinaidi, A. Dynamic instability analysis of magneto-electro-elastic beams with uncertain parameters under static and parametric electric and magnetic fields. Compos. Struct. 2019, 226, 111185. [Google Scholar] [CrossRef]
  33. Zhen, Y.-X.; Wen, S.-L.; Tang, Y. Free vibration analysis of viscoelastic nanotubes under longitudinal magnetic field based on nonlocal strain gradient Timoshenko beam model. Phys. E Low-Dimens. Syst. Nanostructures 2019, 105, 116–124. [Google Scholar] [CrossRef]
  34. Arefi, M.; Kiani, M.; Civalek, Ö. 3-D magneto-electro-thermal analysis of layered nanoplate including porous core nanoplate and piezomagnetic face-sheets. Appl. Phys. A 2020, 126, 76. [Google Scholar] [CrossRef]
  35. Mirjavadi, S.S.; Forsat, M.; Nikookar, M.; Barati, M.R.; Hamouda, A. Nonlinear forced vibrations of sandwich smart nanobeams with two-phase piezo-magnetic face sheets. Eur. Phys. J. Plus 2019, 134, 508. [Google Scholar] [CrossRef]
  36. Mirjavadi, S.S.; Forsat, M.; Barati, M.R.; Abdella, G.M.; Hamouda, A.M.S.; Afshari, B.M.; Rabby, S. Post-buckling analysis of piezo-magnetic nanobeams with geometrical imperfection and different piezoelectric contents. Microsyst. Technol. 2019, 25, 3477–3488. [Google Scholar] [CrossRef]
  37. Jalaei, M.; CivalekÖ. On dynamic instability of magnetically embedded viscoelastic porous FG nanobeam. Int. J. Eng. Sci. 2019, 143, 14–32. [Google Scholar] [CrossRef]
  38. Ghane, M.; Saidi, A.R.; Bahaadini, R. Vibration of fluid-conveying nanotubes subjected to magnetic field based on the thin-walled Timoshenko beam theory. Appl. Math. Model. 2020, 80, 65–83. [Google Scholar] [CrossRef]
  39. Żur, K.; Arefi, M.; Kim, J.; Reddy, J. Free vibration and buckling analyses of magneto-electro-elastic FGM nanoplates based on nonlocal modified higher-order sinusoidal shear deformation theory. Compos. Part B Eng. 2020, 182, 107601. [Google Scholar] [CrossRef]
  40. Sidhardh, S.; Ray, M. Flexomagnetic response of nanostructures. J. Appl. Phys. 2018, 124, 244101. [Google Scholar] [CrossRef]
  41. Zhang, N.; Zheng, S.; Chen, D. Size-dependent static bending of flexomagnetic nanobeams. J. Appl. Phys. 2019, 126, 223901. [Google Scholar] [CrossRef]
  42. Malikan, M.; Eremeyev, V.A. Free Vibration of Flexomagnetic Nanostructured Tubes Based on Stress-driven Nonlocal Elasticity. In Engineering Design Applications; Springer Science and Business Media LLC.: Berlin/Heidelberg, Germany, 2020; Volume 134, pp. 215–226. [Google Scholar]
  43. Malikan, M.; Eremeyev, V.A. On the geometrically nonlinear vibration of a piezo-flexomagnetic nanotube. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  44. Malikan, M.; Eremeyev, V.A. On Nonlinear Bending Study of a Piezo-Flexomagnetic Nanobeam Based on an Analytical-Numerical Solution. Nanomaterials 2020, 10, 1762. [Google Scholar] [CrossRef] [PubMed]
  45. Malikan, M.; Uglov, N.S.; Eremeyev, V.A. On instabilities and post-buckling of piezomagnetic and flexomagnetic nanostructures. Int. J. Eng. Sci. 2020, 157, 103395. [Google Scholar] [CrossRef]
  46. Song, X.; Li, S.-R. Thermal buckling and post-buckling of pinned–fixed Euler-Bernoulli beams on an elastic foundation. Mech. Res. Commun. 2007, 34, 164–171. [Google Scholar] [CrossRef]
  47. Reddy, J. Nonlocal nonlinear formulations for bending of classical and shear deformation theories of beams and plates. Int. J. Eng. Sci. 2010, 48, 1507–1518. [Google Scholar] [CrossRef]
  48. Barretta, R.; Feo, L.; Luciano, R.; De Sciarra, F.M. Variational formulations for functionally graded nonlocal Bernoulli-Euler nanobeams. Compos. Struct. 2015, 129, 80–89. [Google Scholar] [CrossRef]
  49. Akgöz, B.; Civalek, Ö. Buckling analysis of functionally graded microbeams based on the strain gradient theory. Acta Mech. 2013, 224, 2185–2201. [Google Scholar] [CrossRef]
  50. Sarparast, H.; Ebrahimi-Mamaghani, A.; Safarpour, M.; Ouakad, H.M.; Dimitri, R.; Tornabene, F. Nonlocal study of the vibration and stability response of small-scale axially moving supported beams on viscoelastic-Pasternak foundation in a hygro-thermal environment. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  51. Lim, C.; Zhang, G.; Reddy, J. A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation. J. Mech. Phys. Solids 2015, 78, 298–313. [Google Scholar] [CrossRef]
  52. Ansari, R.; Sahmani, S.; Arash, B. Nonlocal plate model for free vibrations of single-layered graphene sheets. Phys. Lett. A 2010, 375, 53–62. [Google Scholar] [CrossRef]
  53. Khorshidi, M.A. The material length scale parameter used in couple stress theories is not a material constant. Int. J. Eng. Sci. 2018, 133, 15–25. [Google Scholar] [CrossRef]
  54. She, G.L. Wave propagation of FG polymer composite nanoplates reinforced with GNPs. Steel Compos. Struct. 2020, 37, 27–35. [Google Scholar]
  55. She, G.L.; Liu, H.B.; Karami, B. On resonance behavior of porous FG curved nanobeams. Steel Compos. Struct. 2020, 36, 179–186. [Google Scholar]
  56. Malikan, M.; Dimitri, R.; Tornabene, F. Transient response of oscillated carbon nanotubes with an internal and external damping. Compos. Part B Eng. 2019, 158, 198–205. [Google Scholar] [CrossRef]
  57. Malikan, M.; Krasheninnikov, M.; Eremeyev, V.A. Torsional stability capacity of a nano-composite shell based on a nonlocal strain gradient shell model under a three-dimensional magnetic field. Int. J. Eng. Sci. 2020, 148, 103210. [Google Scholar] [CrossRef]
  58. Malikan, M. On the plastic buckling of curved carbon nanotubes. Theor. Appl. Mech. Lett. 2020, 10, 46–56. [Google Scholar] [CrossRef]
  59. Malikan, M.; Eremeyev, V.A. Post-critical buckling of truncated conical carbon nanotubes considering surface effects embedding in a nonlinear Winkler substrate using the Rayleigh-Ritz method. Mater. Res. Express 2020, 7, 025005. [Google Scholar] [CrossRef]
  60. Karami, B.; Janghorban, M.; Rabczuk, T. Dynamics of two-dimensional functionally graded tapered Timoshenko nanobeam in thermal environment using nonlocal strain gradient theory. Compos. Part B Eng. 2020, 182, 107622. [Google Scholar] [CrossRef]
  61. Apuzzo, A.; Barretta, R.; Faghidian, S.A.; Luciano, R.; De Sciarra, F.M. Nonlocal strain gradient exact solutions for functionally graded inflected nano-beams. Compos. Part B Eng. 2019, 164, 667–674. [Google Scholar] [CrossRef]
  62. Tang, H.; Li, L.; Hu, Y. Buckling analysis of two-directionally porous beam. Aerosp. Sci. Technol. 2018, 78, 471–479. [Google Scholar] [CrossRef]
  63. Barretta, R.; De Sciarra, F.M. Constitutive boundary conditions for nonlocal strain gradient elastic nano-beams. Int. J. Eng. Sci. 2018, 130, 187–198. [Google Scholar] [CrossRef]
  64. Barretta, R.; De Sciarra, F.M. Variational nonlocal gradient elasticity for nano-beams. Int. J. Eng. Sci. 2019, 143, 73–91. [Google Scholar] [CrossRef]
  65. De Sciarra, F.M.; Barretta, R. A new nonlocal bending model for Euler-Bernoulli nanobeams. Mech. Res. Commun. 2014, 62, 25–30. [Google Scholar] [CrossRef]
  66. Romano, G.; Barretta, R.; Diaco, M.; De Sciarra, F.M. Constitutive boundary conditions and paradoxes in nonlocal elastic nanobeams. Int. J. Mech. Sci. 2017, 121, 151–156. [Google Scholar] [CrossRef]
  67. Pradhan, S.; Reddy, G. Buckling analysis of single walled carbon nanotube on Winkler foundation using nonlocal elasticity theory and DTM. Comput. Mater. Sci. 2011, 50, 1052–1056. [Google Scholar] [CrossRef]
  68. Wang, C.M.; Zhang, Y.Y.; Ramesh, S.S.; Kitipornchai, S. Buckling analysis of micro- and nano-rods/tubes based on nonlocal Timoshenko beam theory. J. Phys. D Appl. Phys. 2006, 39, 3904–3909. [Google Scholar] [CrossRef]
  69. Ansari, R.; Sahmani, S.; Rouhi, H. Rayleigh-Ritz axial buckling analysis of single-walled carbon nanotubes with different boundary conditions. Phys. Lett. A 2011, 375, 1255–1263. [Google Scholar] [CrossRef]
  70. Duan, W.H.; Wang, C.M. Exact solutions for axisymmetric bending of micro/nanoscale circular plates based on nonlocal plate theory. Nanotechnology 2007, 18, 385704. [Google Scholar] [CrossRef]
  71. Duan, W.H.; Wang, C.M.; Zhang, Y.Y. Calibration of nonlocal scaling effect parameter for free vibration of carbon nanotubes by molecular dynamics. J. Appl. Phys. 2007, 101, 024305. [Google Scholar] [CrossRef]
Figure 1. Geometry and description of a continuum nanobeam as a square actuator installed on simple end conditions.
Figure 1. Geometry and description of a continuum nanobeam as a square actuator installed on simple end conditions.
Symmetry 12 01935 g001
Figure 2. (a) Nonlocal parameter vs. four cases of non-porous nanobeams (l = 0.5 nm, L = 10 h, ψ = 1 mA). (b) The length scale strain gradient parameter vs. different cases of non-porous nanobeams (e0a = 0.5 nm, L = 10 h, ψ = 1 mA).
Figure 2. (a) Nonlocal parameter vs. four cases of non-porous nanobeams (l = 0.5 nm, L = 10 h, ψ = 1 mA). (b) The length scale strain gradient parameter vs. different cases of non-porous nanobeams (e0a = 0.5 nm, L = 10 h, ψ = 1 mA).
Symmetry 12 01935 g002
Figure 3. Types of porosity vs. two cases of nanobeams (e0a = 0.5 nm, l = 1 nm, L = 20 h, ψ = 1 mA).
Figure 3. Types of porosity vs. two cases of nanobeams (e0a = 0.5 nm, l = 1 nm, L = 20 h, ψ = 1 mA).
Symmetry 12 01935 g003
Figure 4. Magnetic potential vs. two cases of non-porous nanobeams (e0a = 0.5 nm, l = 1 nm, L = 10 h).
Figure 4. Magnetic potential vs. two cases of non-porous nanobeams (e0a = 0.5 nm, l = 1 nm, L = 10 h).
Symmetry 12 01935 g004
Figure 5. Thickness vs. two cases of non-porous nanobeams (e0a = 0.5 nm, l = 1 nm, ψ = 1 mA, L = 20 h).
Figure 5. Thickness vs. two cases of non-porous nanobeams (e0a = 0.5 nm, l = 1 nm, ψ = 1 mA, L = 20 h).
Symmetry 12 01935 g005
Table 1. Axial porosity distribution patterns.
Table 1. Axial porosity distribution patterns.
Porosity Type λ ( x ) Ranges   of   α
“O” type distribution η 1 2 sin ( π L x ) 0 α < 0.344
O ¯ ” type distribution η 2 2 [ 1 sin ( π L x ) ] 0 α < 0.112
“X” type distribution η 1 η 2 sin ( π L x ) 0 α < 0.197
X ¯ ” type distribution η 1 η 2 [ 1 sin ( π L x ) ] 0 α < 0.197
Uniform type distribution1 0 α < 0.85
η 1 = π 2 , η 2 = π π 2
Table 2. Constitutive boundary conditions for S-S.
Table 2. Constitutive boundary conditions for S-S.
Nonlocal Strain Gradient Conditions at
(0, L)
Local Conditions
( l = μ = 0 )   at   ( 0 ,   L )
w = 0
M n l = ( 1 + l 2 d 2 d x 2 ) M l + μ d d x ( d 2 T x x z d x 2 + N x 0 d 2 w d x 2 ) = 0 *
T x x z = B d 2 w d x 2 + f 31 ψ = 0
w = 0
M l = D d 2 w d x 2 = 0 *
T x x z = B d 2 w d x 2 + f 31 ψ = 0
* Sub-indexes (nl and l) are nonlocal and local phases, respectively.
Table 3. Critical loads’ validation with literature (E = 1 TPa, ν = 0.19, d = 1 nm).
Table 3. Critical loads’ validation with literature (E = 1 TPa, ν = 0.19, d = 1 nm).
PCr (nN)
L (nm)µ = 0 nm2µ = 1 nm2µ = 4 nm2
[67][68]Present[67][68]Present[67][68]Present
104.84474.84474.84474.40954.40954.40953.47353.47353.4735
123.36443.36443.36443.14863.14863.14862.64052.64052.6405
142.47182.47182.47182.35332.35332.35332.05742.05742.0574
161.89251.89251.89251.82221.82221.82221.63961.63961.6396
181.49531.49531.49531.45111.45111.45111.33291.33291.3329
201.21121.21121.21121.1821.1821.1821.10241.10241.1024
Table 4. Material parameters of the piezo-flexomagnetic nanobeam.
Table 4. Material parameters of the piezo-flexomagnetic nanobeam.
CoFe2O4
C11 = 286 GPa
q31 = 580.3 N/A.m
a33 = 1.57 × 10−4 N/A2
(A = Ampere)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Malikan, M.; Eremeyev, V.A.; Żur, K.K. Effect of Axial Porosities on Flexomagnetic Response of In-Plane Compressed Piezomagnetic Nanobeams. Symmetry 2020, 12, 1935. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121935

AMA Style

Malikan M, Eremeyev VA, Żur KK. Effect of Axial Porosities on Flexomagnetic Response of In-Plane Compressed Piezomagnetic Nanobeams. Symmetry. 2020; 12(12):1935. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121935

Chicago/Turabian Style

Malikan, Mohammad, Victor A. Eremeyev, and Krzysztof Kamil Żur. 2020. "Effect of Axial Porosities on Flexomagnetic Response of In-Plane Compressed Piezomagnetic Nanobeams" Symmetry 12, no. 12: 1935. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121935

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop