1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

buckling of thin composite plates reinforced with randomly oriented, straight singlewalled carbon nanotubes using B3-Spline finite strip method
H. Foroughia, H.R. Askariyehb and M. Azharia, 1
a

Department of Civil Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran
Department of Civil Engineering, Faculty of Engineering, Yazd University, Yazd,P.O.BOX 89195-741, Iran

b

Abstract
This paper is devoted to the buckling analysis of thin composite plates under straight single-walled
carbon nanotubes reinforcement with uniform distribution and random orientations. In order to
develop the fundamental equations, the B3-Spline finite strip method along with the classical plate
theory (CPT) is employed and the total potential energy is minimized which leads to an eigenvalue
problem. For deriving the effective modulus of thin composite plates reinforced with carbon
nanotubes, the Mori-Tanaka method is used in which each straight carbon nanotube is modeled as a
fiber with transversely isotropic elastic properties. The numerical results including the critical
buckling loads for rectangular thin composite plates reinforced by carbon nanotubes with various
boundary conditions and different volume fractions of nanotubes are provided and the positive effect
of using carbon nanotubes reinforcement in buckling of thin plates is illustrated.
Keywords: mechanical buckling; thin composite plates; single-walled carbon nanotubes, B3-Spline
finite strip method

1

Corresponding Author:

Mojtaba Azhari, Department of Civil Engineering, Isfahan University of Technology, Isfahan, Iran.
Phone: +98 3113913804; Fax: ++98 3113912700; Email: mojtaba@cc.iut.ac.ir

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

1.Introduction
In recent years, there has been a growing trend towards using composite materials in many
engineering fields. Many researchers have attempted to reinforce composite structures by
high strength reinforcements, among which the carbon nanotubes are widely used.
In this regard, many researchers have elaborated on the material properties of the
composites reinforced with CNTs. Among these investigations, it can be referred to Fukuda
and Kawatawho evaluated the Young modulus of composite structures reinforced by short
fibers [1]. Also Ajayan et al. studied the polymer composites reinforced by aligned CNT
arrays [2]. Griebel and Hamaekersused the molecular dynamics simulations to evaluate the
elastic moduli of CNT reinforced composite structures [3]. Esawi and Faragshowed that
reinforcing composite structures by CNTs leads to a considerably high stiffness while
economical [4].
In order to examine the behavior of CNT reinforced composites in bending and buckling,
several researches have been conducted. Vodenitcharova and Zhang studied the bending and
local buckling of a nanocomposite beam reinforced by a single-walled CNT by employing
the Airy stress-function method to analyze the deformation of the matrix [5]. Shen studied the
postbuckling of nanotube-reinforced composite cylindrical shells in thermal environments
[6]. Zhu et. Al discussed the static and free vibration of composite plates under CNT
reinforcement by employing the finite element method [7].
Finite strip method is a modified finite element method introduced by Cheung and has been
obtained by discretizing the plate into some finite strips [8]. Bradford and Azhari investigated
the buckling of plates with different end conditions using the finite strip method [9].Azhari et
al. applied the spline finite strip method augmented by the bubble functions for elastic
analysis of rectangular plates to show that efficiently improves the convergence of the spline
finite strip method with respect to strip subdivision [10]. Lotfi et al. studied the inelastic local
buckling of skew thin thickness-tapered plates using Spline finite strip method [11].
In the current study, the buckling analysis of thin composite plates reinforced with singlewalled CNTs is conducted. The distribution of nanotubes in the matrix is uniform and they
are randomly oriented. In order to derive the fundamental equations to be applied in the
buckling relations, the classical plate theory and B3-Spline finite strip method is used. It is
noteworthy that the characterization of the CNTs properties is obtained through the MoriTanaka approach.
2. Eshelby-Mori-Tanaka approach
The Mori-Tanaka approach is introduced as an analytical method to characterize a linear
elastic polymer matrix reinforced by a large number of dispersed CNTs that are aligned or
randomly oriented, straight and of infinite length.
This method assumes that each nanotube is embedded in an infinite matrix subjected to an
effective stress ( σ m ) or an effective strain ( ε m ).According to Benveniste the effective elastic
tensor is revised as follows [12]

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

C  C m V CN ((CCN  Cm ).A ).[V m I V CN A ]1

(1)

where I is the fourth-order identity tensor, the subscripts m and CN denote the quantities of
the matrix and the nanotubes, respectively, V m and V CN stand for the volume fractions, and
Cm and CCN denote the elastic moduli tensors and A is the fourth-order tensor depending on
the characteristics of the carbon nanotube.
in which
A  [I  S.(Cm ) 1.(CCN  Cm )]1

(2)

where S is the Eshelby tensor derived by Mura [13].
In this study, the CNTs are assumed to be randomly oriented in the matrix. The orientation of
a straight CNT is described by two Euler angles  and  . In order to transform in
global(x,y,z) coordinates to local (x',y',z') coordinates, the transformation matrix g is used as
follows
 cos 
g   sin 
 0

 cos  sin 
cos  cos 
sin 

sin  sin  
 sin  cos  
cos 


(3)

Moreover, the orientation distribution of CNTs in a composite is characterized by a
probability density function  ( ,  ) satisfying the normalization condition as below
2

 
0



2
0

 ( ,  ) sin  d  d   1

(4)

In the above relation, in case that the orientations of CNTs are completely random, the
density function is defined as  ( ,  )   / 2 .
Considering the constitutive equations of the matrix phase, the strain and the stress vectors
εCN ( ,  ) and σCN ( ,  ) are related to the stress matrix σ m as below
εCN ( ,  )  A ( ,  ).Cm1.σ m , σCN ( ,  )  [C r .A ( ,  ).C m1 ].σ m

(5)

The average strain and stress in all randomly oriented CNTs can be written as

 2 

avg
εCN
    2  ( ,  ) A( ,  )sin  d  d   .ε m ,
0
0


 2 

avg
σCN
    2  ( ,  )[CCN .A( ,  ).Cm1 ]sin  d  d   .σ m
(6)
0
0


where, the superscript "avg" denotes the average over special orientations. According the
average theorems and substituting Eq. (9) in to the constitutive relation σ  C.ε , the effective
elastic modulus tensor for randomly oriented CNTs in the matrix will be derived, which is not
given here for brevity.
According to the Mori-Tanaka approach, the composite is considered isotropic, and so its
bulk modulus K and shear modulus G are defined as
V (  3K m CN )
V (  2G m CN )
(7)
K  K m  CN CN
, G  G m  CN CN
2(V m V CN CN )
3(V m V CN CN )

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

where K m and G m are the bulk and shear moduli of the matrix phase, respectively, and

CN 

3( K m  G m )  k CN  lCN
,
3(G m  k CN )

2 G m (3K m  G m )  G m (3K m  7G m ) 
1  4G m  2k CN  lCN
4G m


,
5  3(G m  k CN )
G m  pCN
G m (3K m  G m )  mCN (3K m  7G m ) 

CN  
1

CN   nCN  2lCN 
3


1 2

(2k CN  lCN )(3K m  2G m  lCN ) 
,
G m  k CN

8G p

8m G (3K  4G )

CN m
m
m

CN   (nCN  lCN )  m CN 
5 3
Gm  pCN 3K m (mCN Gm ) Gm (7mCN Gm )

2(k CN  lCN )(2Gm  lCN ) 

3(Gm  k CN )


(8)

where k CN , lCN , mCN , nCN and pCN are taken from analytical results of Popov et. al (2000).
Finally, the effective Young’s modulus E and Poisson’s ratio  of the composite in terms of
the bulk modulus K and shear modulus G is defined as
9KG
3K  2G
(9)
E 
,
3K  G
6K  2G
The constitutive stress-strain relations are then introduced as

 xx  Q11 Q12
0   xx 
E (z )
  
 
,Q12   Q11 ,Q 66  G
 yy   Q12 Q 22 0   yy  , Q11  Q 22 
2

1

   0
0 Q 66   xy 
 xy  

(10)

3. Mechanical buckling
3.1. The displacement field

In this section, mechanical buckling of thin composite plates reinforced with single-walled
CNTs is discussed. Based on the classical plate theory (CPT), the displacement field is
defined as
w
w
u ( x, y , z )  u0 ( x, y )  z 0 v ( x, y , z )  v0 ( x, y )  z 0 w( x, y , z )  w0 ( x, y )
(11)
x
y
where u 0 , v 0 and w 0 are the mid-plane displacements .
3.2. B3-Spline Finite strip approach

According to the classical plate theory for analysis thin plates, the ensuring displacement
fields for one strip is
 win 
T

  in  δ in
 w 
 w  
 n   wi , 
 , wj,
    w jn   δ jn
 x  i
 x  j  n 


  jn 

 

(12)

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

nodal lines

h
h

b
h

m sections

h
Fig. 1.nodal lines for one strip in Spline finite strip method

A generic function f ( y ) , say a generalized, can be represented as a linear combination of
local B3- Spline functions  i ( y ) and parameters  i for the jth node, according to the
expression:
m 1

f ( y )    i i ( y )

(13)

i 1

where
 (y  y i 2 )3
2
3
 3
2
 1   h 3h  y  y i 1  3h  y  y i 1  3 y  y i 1 
 i (y )   3   3
2
3
 6h  h 3h 2  y i 1  y  3h  y i 1  y  3 y i 1  y 
  y i 2  y 

y i 2  y  y i 1
y i 2  y  y i 1,

(14)

y  y  y i 1,
y i 1  y  y i 2 ,

The shape function,w ij could be introduced as
w

ij

 N i ( x ) j ( y )

(15)

In addition, N i  x  are the following Hermitian cubic polynomials.
2
3
2
2
3

 x 2  x
x 
x 
x  x  
x 
x 
  2   , N 2 x  x 1  2       , N 3 x  3    2   , N 4 x  x    

 b   b
b 
b 
 b   b  
b 
b 

N 1  x   1  3

 

 

 



 

(16)

The out of plane buckling displacement of the plate,  w  can be expressed as
4

m 1

w (x , y )     N i (x ). j ( y ). ij



(17)

i 1 j 1

    w 1, 1 ,w 1,0 ,w 1,1 ,...,w 1, m 1,  1, 1 , 1,0 , 1,1 ,..., 1, m 1, w 1, 1 ,w 1,0 ,w 1,1 ,...,w 1,m 1 ,  2, 1 , 2,0 , 2,1 ,..., 2,m 1 T (18)
In Eq.(19) and Eq.(20),    , is the nodal displacement vector. Furthermore, the following
strain- displacement relations are employed.

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

 ij 

1
 ui, j  u j ,i  uk ,iuk , j 
2

(19)

Assuming that the plate is thin based on Kirchhorf hypothesis [14], the vectors of internal
moments, and infinitesimal buckling curvatures, may be given by well-known expressions:
T
T
  f    M x , M y , M xy  ,   f     x ,  y ,  xy 

(20)

In which

 2w
 2w
 2w
x 
,y 
,  xy 
x 2
y 2
x y

(21)

4. The derivation of the critical buckling load
On the basis of the present approach and according to the total potential theorem, the total
potential energy of a plate can be expressed as

U total  U  Vp

(22)

in which U is the strain energy of the plate and Vp is potential energy of the external loads.
According the procedure of minimizing the total potential energy in the standard finite strip
method, the strain energy and stiffness matrix can obtain as
U

1
1
T
ε f T σ f dA   δ K δ, K    B f



2
2

T  Q  B f  dA

(23)

in which
  2   N 

2
 2 x
    N 
Bf  
2
 2 y
    N 

x y


 



 


 

 3c

(24)

where K is the stiffness matrix of the plate.
Following the procedure typically used in B3-Spline finite strip method, the geometry
stiffness matrix of the strip for a plate subjected to  x0 ,  y0 and  xy0 in-plane stresses
Vp 

 0 w 2
1
1
0 w 2
0 w w 
T
 x ( )   y ( ) 2 xy
 dV   δ K g δ

2
x
y
x
y
2







where K g is the geometric stiffness matrix and may be obtained as

(25)

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

 K g   t   B g    g  B g 
T

(26)

In which
    N 

 B g       Nx


y


 



 

 2c

,

  x0  xy0 
0
0 
  xy  y 

 g   

(27)

All of the equations are developed for one strip. After assembling the stiffness matrix of the
plate, and the geometric stiffness matrix over all the strips of the plate, followed by imposing
the boundary conditions. The critical buckling load of a composite plate reinforced by CNTs,
defined as  cr , can be determined by solving the following eigenvalue problem using
standard computer subroutines.
K   cr K g  0

(28)

4. Numerical results
In this section, several numerical examples are presented and discussed to illustrate the
efficiency of the B3-Spline finite strip method in mechanical buckling analysis of thin
rectangular plates reinforced with single-walled CNTs .The matrix in which the CNTs are
embedded is PmPV (Poly{(m-phenylenevinylene)-co-[(2.5-dioctoxy-p-phenylene)vinylene]},
the material properties of which are  m  0.34 and E m  2.1 GPa at room temperature of
T  300 K .
In order to validate the results obtained by the present method, the normalized critical
buckling load (  ) for composite plates with no CNT (V CN  0 ) are presented in Table 1 and
compared with those obtained by Bradford and Azhari [13]. In Table 1, the results are
provided for various boundary conditions. As it can be observed from Table 1, the results
obtained from the present method have excellent accuracy.
Table 2. Normalized critical buckling factor (  

b

  cr 12(1 2 )b 2 
) for different boundary conditions
  2 Eh 3 
Boundary conditions

h

Methods

SSSS

SCSC

SCSS

CCCC

present

4.000

7.721

5.979

10.072

Bradford et al[13]

4.000

7.72

5.97

10.08

0.01

Table 3illustrates the values of  for CNT reinforced composite plates of Fig. 1 for different
nanotube volumes (V CN ). Here again the results for various boundary are discussed. As it was
expected, by increasing the volume fraction of CNTs, the normalized critical buckling load
increases as a result. Moreover, by comparing the results for different boundary conditions, it
can be seen that the normalized critical buckling load for the free boundary conditions,
showed with F in the tables, are considerably lower than those for the clamped boundary
conditions which is quite realistic. Furthermore, as the volume fraction of CNTs increases,

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

the normalized critical buckling loads increases which shows the influence of CNTs in
increasing the stiffness of the plate.
  cr 12(1 m 2 )b 2 
Table 3. Normalized critical buckling loads (  
) for various carbon nanotubes volume
  2E m h 3 
fraction and different boundary conditions

V CN

b /h

SCSC

SCSS

SSSS

SFSF

100

10.512

7.768

5.390

1.271

50

10.320

7.711

5.375

1.269

0.01
0.05
0.1

100

21.351

15.774

10.944

2.607

50

20.960

15.658

10.913

2.604

100

35.167

25.978

18.023

2.607

50

34.522

25.787

17.972

2.604

In Fig 2, the critical buckling load ( N cr ) versus aspect ratio ( a / b ) for composite plates
reinforced with three CNTs volume fractions of 0.01, 0.05 and 0.1 is depicted. It is observed
that the normalized critical buckling load increases with the volume fraction of CNTs.
Moreover, as proved in the previous studies, by the increase in the aspect ratio, the critical
buckling coefficient converges.

Ncr (KN/M)

60

CNTs %10

50

CNTs %5

40

CNTs %1

30
20
10
0
0.5

1

1.5

a/b

2

2.5

3

Fig 2. Critical buckling load of composite plate reinforced with various volume fraction of carbon nanotubes
versus aspect ratio

5. Conclusions
In this paper, the B3-Spline finite strip method is successfully applied to the mechanical
buckling analysis of thin composite plates reinforced with single-walled carbon nanotubes.
The mechanical properties of carbon nanotubes embedded in the matrix are derived through
the Mori-Tanaka approach. Comparison studies were conducted on the critical buckling load
of composites without nanotubes to verify the present method and the results were found to
be in good agreement with those available in literature. Several numerical examples for
composites under carbon nanotubes reinforcement were provided and it was concluded that
the change in carbon nanotubes volume fraction, boundary conditions, aspect ratio, thickness
to width strongly affect the critical buckling load of the plate.

1st National Conference on New Materials and Systems in Civil Engineering, Graduate University of Advanced Technology, Feb 2013

References
[1] Fukuda, H., Kawata, K. (1974).
composites.”Fibre.Sci.Technol.,7, 207-222.

“On

Young’s

modulus

of

short

fibre

[2] Ajayan, P.M., Stephan, O., Colliex, C., and Trauth, D. (1994).“Aligned carbon nanotube arrays
formed by cutting a polymer resin—nanotube composite.”Science.,265(5176), 1212-1214.
[3] Griebel, M., Hamaekers, J. (2004) “Molecular dynamics simulations of the elastic moduli of
polymer–carbon nanotube composites.”Mech. Eng., 193(17-20), 1773-1788.
[4] Esawi , M. K., Farag, M. (2007). “Carbon Nanotube Reinforced Composites: Potential and
Current Challenges.” Mater.Des.,28(9), 2394-2401.
[5] Vodenitcharova, T., Zhang, L. C. (2006).“Bending and local buckling of a nanocomposite beam
reinforced by a single-walled carbon nanotube.”Int. J. Solids. Struct.,43(10), 3006-3024.
[6] Shen, H. S. (2011). “Postbuckling of nanotube-reinforced composite cylindrical shells inthermal
environments – Part I: Axially-loaded shells.”Compos.Struct., 93(8), 2096-2108
[7] Zhu, P., Lei, Z.X., Liew, K.M. (2012). “Static and free vibration analyses of carbon
nanotubereinforced composite plates using finite element method with first order shear deformation
plate theory.” Compos.Struct., 94(4), 1450-60.
[8] Cheung, Y.K. (1976). “Finite strip method in structural analysis.” Pergamon Press, Oxford-New
York.
[9] Bradford, M.A., Azhari, M. (1995). “Buckling of plates with different end conditions using the
finite strip method.” Comput.Struct., 56(1), 75-83.
[10]Azhari, M., Hoshdar, S., Bradford, M.A. (2000).“On the use of bubble functions in the local
buckling analysis of plate structures by the spline finite strip method.”Int J Numer Methods Eng., 48,
583-593.
[11] .Lotfi , S., Azhari , M., Heidarpour , A. (2011). “Inelastic initial local buckling of skew thin
thickness-tapered plates with and without intermediate supports using the isoparametric spline finite
strip method.”Thin-Walled Structures., 49, 1475-1482.
[12] Benveniste, Y. (1987). “A new approach to the application of Mori–Tanaka’s theory in
composite materials.”Mech Mater., 6(2), 147-157.
[13] Mura, T.(1987). “Micromechanics of Defects in Solids.”MartinusNijhoff Publishers, Dordrecht.
[14] Vrcelj, Z., Bradford, M.A. (2008). “A simple method for the inclusion of external and internal
supports in the spline finite strip method (SFSM) of buckling analysis.”ComputStruct., 86, 529-544.

[15] Foroughi, H., Askariyeh, H.R., Azhari, M. (2013).“Mechanical Buckling of Thick Composite Plates
Reinforced with Randomly Oriented, Straight, Single-Walled Carbon Nanotubes Resting on an Elastic Foundation
using the Finite Strip Method.” Journal of Nanomechanics and Micromechanics., 3, 49-58.

