Representative volume element to estimate buckling behavior of graphene/polymer nanocomposite.  
Jump to Full Text  
MedLine Citation:

PMID: 22994951 Owner: NLM Status: PubMednotMEDLINE 
Abstract/OtherAbstract:

The aim of the research article is to develop a representative volume element using finite elements to study the buckling stability of graphene/polymer nanocomposites. Research work exploring the full potential of graphene as filler for nanocomposites is limited in part due to the complex processes associated with the mixing of graphene in polymer. To overcome some of these issues, a multiscale modeling technique has been proposed in this numerical work. Graphene was herein modeled in the atomistic scale, whereas the polymer deformation was analyzed as a continuum. Separate representative volume element models were developed for investigating buckling in neat polymer and graphene/polymer nanocomposites. Significant improvements in buckling strength were observed under applied compressive loading when compared with the buckling stability of neat polymer. 
Authors:

Avinash Parashar; Pierre Mertiny 
Publication Detail:

Type: Journal Article Date: 20120920 
Journal Detail:

Title: Nanoscale research letters Volume: 7 ISSN: 1556276X ISO Abbreviation: Nanoscale Res Lett Publication Date: 2012 
Date Detail:

Created Date: 20130201 Completed Date: 20130204 Revised Date: 20130411 
Medline Journal Info:

Nlm Unique ID: 101279750 Medline TA: Nanoscale Res Lett Country: United States 
Other Details:

Languages: eng Pagination: 515 Citation Subset:  
Affiliation:

University of Alberta, 49 Mechanical Engineering Building, Edmonton, Alberta, T6G 2G8, Canada. pmertiny@ualberta.ca. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:


Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Nanoscale Res Lett Journal ID (isoabbrev): Nanoscale Res Lett ISSN: 19317573 ISSN: 1556276X Publisher: Springer 
Article Information Download PDF Copyright ©2012 Parashar and Mertiny; licensee Springer. openaccess: Received Day: 26 Month: 6 Year: 2012 Accepted Day: 4 Month: 8 Year: 2012 collection publication date: Year: 2012 Electronic publication date: Day: 20 Month: 9 Year: 2012 Volume: 7 Issue: 1 First Page: 515 Last Page: 515 PubMed Id: 22994951 ID: 3561215 Publisher Id: 1556276X7515 DOI: 10.1186/1556276X7515 
Representative volume element to estimate buckling behavior of graphene/polymer nanocomposite  
Avinash Parashar1  Email: draparashar@yahoo.com 
Pierre Mertiny1  Email: pmertiny@ualberta.ca 
1University of Alberta, 49 Mechanical Engineering Building, Edmonton, Alberta, T6G 2G8, Canada 
In the recent past, graphene has emerged as a potential candidate for developing nanocomposites with improved properties [^{1},^{2}]. The experimental characterization of graphene/polymer nanocomposites is a challenging process, and hence, computational approaches for predicting the behavior of such materials have also extensively been employed. Various multiscale models are available in the literature for predicting the properties of carbon nanotube (CNT)based nanocomposites [^{3}^{}^{5}], but very few models have been presented to study graphene nanocomposites. For example, Cho et al. [^{6}] developed a numerical model in conjunction with MoriTanaka approach to study the elastic constants of randomly distributed graphene in polymer. Awasthi and his team [^{7}] investigated the load transfer mechanism between polyethylene and graphene sheets. Montazeri and Tabar [^{8}] developed a finite element (FE)based multiscale model to investigate the elastic constants of graphenebased nanocomposites.
Buckling in isolated graphene sheets was modeled by several researchers [^{9}^{}^{11}]. However, buckling stability of graphene/polymer nanocomposites was only reported by Rafiee et al. [^{12}]. Using an experimental and analytical approach, up to 50% and 32% improvement in the buckling stability of nanocomposites was reported respectively. In the analytical approach, an Euler buckling formulation was employed, and elastic properties required in the Euler equation were estimated by experimental means. The discrepancies between the two buckling stabilities were attributed to scaling issues.
It is well established that the reinforcement of polymer with graphene increases the elastic modulus of the material which further improves buckling stability. The aim of this study is to propose a numerical model which can estimate the increase in buckling stability with different volume fractions of graphene and can further be extended to complex shapes and structures.
It has been reported that achieving a uniform dispersion of twodimensional graphene sheets in polymer is more challenging compared to the mixing of onedimensional CNT. Moreover, the application of nanocomposites is not limited to simple structures, and the comprehension of material behavior in complex structures is restricted when employing experimental and analytical methods. Consequently, research efforts are increasingly focused on numerical approaches. To overcome some of the limitations that exist in experimental and analytical work, a multiscale representative volume element (RVE) is proposed in this paper to investigate buckling phenomena in graphene/polymer nanocomposites under the assumption that graphene is uniformly distributed in the polymer. To the knowledge of the present authors, no numerical model has been reported yet to study the effect of graphene on the buckling strength of nanocomposites. In the proposed technique, graphene was modeled in the atomistic scale, whereas polymer deformation was analyzed as a continuum.
In this paper a finite element technique was employed in conjunction with molecular and continuum mechanics to simulate buckling in graphene/polymer nanocomposites. In the proposed RVE, the polymer, epoxy in this case, was modeled as a continuum phase whereas the deformation in graphene was evaluated using an atomistic description. Nonbonded interactions were considered as the load transfer mechanism or interphase between polymer and graphene. FE modeling was performed in the ANSYS (Version 13) software environment (ANSYS Inc., Canonsburg, PA, USA).
To model graphene in the proposed RVE, it was assumed that graphene behaves like a space frame structure in which the covalent bonding between CC atoms constitutes the loadbearing element while atoms act like a joint. A molecularmechanicsbased approach was employed to estimate properties for graphene in the atomistic scale. In the molecular form, potential energy of graphene can be represented as
(1)
∑Vgraphene=∑Vnonbonded+∑Vbonded 
In Equation 1V_{bonded} represents the bonded interaction between atoms and is considered primarily responsible for maintaining the structural integrity. The potential energy term ‘V_{nonbonded}’ includes electrostatic forces, e.g., van der Waals forces, which are considered weak in nature and can generally be neglected in comparison to bonded interactions. Li and Chou [^{13}] established a correlation between interatomic molecular potential energies and corresponding strain energies of a beam element as described by Equations 2, 3, and 4.
(2)
EAL=KR 
(3)
EIL=Kθ 
(4)
GJL=Kτ 
where K_{r} (938 kcalmol^{−1}Å^{−2}), K_{θ} (126 kcalmol^{−1}rad^{−2}), and K_{τ} (40 kcalmol^{−1}rad^{−2}) are the bond stretching, bond bending, and bond torsional resistance force constants. E represents the Young's modulus; A, the cross sectional area; I, the moment of inertia; G, the shear modulus; J, the polar moment of inertia of the beam element, respectively. BEAM4 elements were used in the ANSYS software environment to model graphene by connecting nodes, and material properties for those elements were estimated with the help of Equations 2, 3, and 4.
The volume fraction of graphene in polymer ranges commonly up to 10%. Most of the material volume is therefore occupied by polymer. Simulating the polymer phase on the atomistic scale would require large efforts in dealing with large numbers of degrees of freedom as well as substantial computational cost. Therefore, as a reasonable compromise, the polymer phase was modeled as a continuum, and threedimensional SOLID45 elements were used for meshing the geometry. Epoxy with a Young's modulus of 3.4 GPa and a Poisson’s ratio of 0.42 was considered as the polymer material in the present work.
In this research paper, nonbonded interactions were considered as load transfer mechanisms between the polymer and graphene. Naturally, load transfer between nonfunctionalized graphene sheets and polymer takes place through van der Waals interactions. The Lennard Jones ‘612’ potential given in Equation 5 was employed to estimate the properties for the interface region, which is depicted in the schematic in Figure 1.
(5)
Ur=4γϕr12+ϕr6 
where r is the atomistic distance between atoms, φ (0.34 nm) is the hard sphere radius, and γ (0.556 kcal/mole) is the potential well depth. In the finite element environment, these nonbonding interactions were modeled with the help of LINK8 truss elements. As shown in Equation 6, properties for those truss elements were estimated by comparing the classical continuum strain energy of truss elements with the Lennard Jones ‘612’ potential. Note that the truss model described herein was earlier employed by Li and Chou [^{14}] to simulate the interaction between CNT and epoxy polymer.
(6)
Er=8γREqAr−REq2ϕr12+ϕr6 
where R_{Eq} is the initial undeformed length, r is the deformed length, A is the cross sectional area, and E(r) is the Young's modulus of the truss element.
A linear analysis for mode one buckling was performed in this work, which is associated with the computation of a bifurcation load and corresponding buckling mode. The analysis in the finite element environment was divided into two sections, i.e., a prebuckling and a postbuckling analysis [^{15}].
The prebuckling analysis was performed in the finite element environment to compute the reference stresses S* within each element. These computed reference values are then used for estimating the geometric stiffness matrix K_{G}, which is needed in the postbuckling analysis. To initiate the analysis in the finite element model, the general matrix form equation provided in Equation 7 was employed to estimate nodal displacements [Q] when the structure is subjected to the reference unit load P_{R}. In Equation 7 the total stiffness matrix [K_{0}] is the assembled global stiffness matrix for the multiscale structure composed of beam, truss, and solid elements. Conventional assembly techniques can be used to obtain the global stiffness matrix.
(7)
K0Q=PR 
Equation 8 represents the general solution that is obtained in the finite element analysis after solving the assembled matrices defined in Equation 7.
(8)
Q*=K0−1PR 
The nodal displacement vectors {Q*} obtained from Equation 8 are further postprocessed to compute strains and corresponding stress values S* within each element using Equation 9, in which matrix [D] represents the elasticity matrix (material property) while matrix [B] is a function of approximating polynomial or shape functions.
(9)
S*=DBQ* 
The postbuckling analysis corresponds to a general eigenvalue problem as defined by Equation 10, in which λ and v are the load factor and the eigenvector of displacements. The geometric stiffness matrix K_{G} used in Equation 10 is defined with the help of Equation 11, where Γ is composed of derivatives of shape functions and S is a function of the stress S* estimated in the prebuckling analysis.
(10)
K0+λKGν=0 
(11)
KG=∫VΓTSΓdV 
By solving Equation 10, the lowest possible value for the load factor λ_{min} can be determined, which is then employed to estimate the buckling force given by Equation 12.
(12)
Bucklingforce=λmin×Referenceload 
The proposed RVE was employed to understand the buckling behavior of graphene in a polymer matrix when graphene was assumed to be uniformly distributed. A second separate RVE structure with equal dimensions was developed using finite element modeling to study the buckling in neat (homogenous) polymer.
To validate the results of the proposed RVE, an alternative method (indirect method) was also employed in this paper. The proposed RVE model with boundary conditions as shown in Figure 2 was first used to estimate the Young's modulus of the RVE as a whole for various volume fractions of graphene. The estimated material properties (E) were then used in defining the homogenous rectangular plate with dimensions of the proposed RVE model. The buckling loads were calculated for this homogenous rectangular plate shown in Figure 3a and were compared with the direct approach defined in Figure 3b.
The analytical formulation to estimate the buckling load for a rectangular plate is given in Equation 13, where P_{x} is the applied unidirectional compressive load, w is the displacement in the outward normal direction, and D is flexural rigidity (i.e., a function of the Young's modulus). Boundary conditions were kept identical in all RVE models. Only the value of flexure rigidity was varied, which afforded the proposed RVE model the capability of yielding E for different graphene volume fractions.
(13)
D∇4wx,y+Px∂2w∂2x=0 
The boundary conditions along with dimensions for the proposed RVE model are shown in Figure 3. The thickness of graphene in the atomistic scale and for epoxy as a continuum phase was kept constant at 0.344 nm, whereas the thickness of the interphase was kept at 0.172 nm according to [^{14}]. The graphene volume fraction in the proposed RVE model was varied by changing the size of the graphene sheet, leading to graphene volume fractions ranging from 2% to 6%. Filler volume fractions of reasonable and practical magnitude were thus studied, omitting the agglomeration effects in graphene nanocomposites with high filler content.
The results obtained from the developed RVE structures were plotted in Figure 4. These data show a significant improvement in the buckling performance of nanocomposites under compressive loading. The buckling strength of neat epoxy was herein considered as the reference level. Buckling strengths that are calculated from the direct and indirect approach are in good agreement, which validates the proposed numerical technique. In the current work, up to 26% improvement in the buckling strength of epoxy was estimated for only 6% volume fraction of graphene.
In this study a representative volume element method was successfully employed to investigate the buckling phenomenon in graphene/polymer nanocomposites, where graphene was assumed to be uniformly distributed. Graphene was modeled in the atomistic scale and polymer as a continuum. A significantly enhanced buckling strength of graphene reinforced polymers was observed as compared to neat polymer, i.e., buckling strength of graphene/polymer nanocomposite improved by 26% with only 6% filler volume fraction.
The authors declare that they have no competing interests.
Both authors made equally valuable contributions to this paper. Both authors read and approved the final manuscript.
The authors would like to acknowledge the financial support by the Alberta Ingenuity Graduate Scholarship in Nanotechnology.
References
Rafiee MA,Rafiee J,Srivastave I,Wang Z,Song H,Yu ZZ,Koratkar N,Fracture and fatigue in graphene nanocompositesSmallYear: 2010617918310.1002/smll.20090148019924737  
Rafiee MA,Rafiee J,Wang Z,Song H,Yu ZZ,Koratkar N,Enhanced mechanical properties of nanocomposites at low graphene contentACS NanoYear: 200333884389019957928  
Ayatollahi MR,Shadlou S,Shokrieh MM,Multiscale modeling for mechanical properties of carbon nanotube reinforced nanocomposites subjected to different types of loadingCompos StructYear: 2011932250225910.1016/j.compstruct.2011.03.013  
Meguid SA,Wernik JM,Cheng ZQ,Atomistic based continuum representation of the effective properties of nanoreinforced epoxiesInt J Solids StructYear: 2010471723173610.1016/j.ijsolstr.2010.03.009  
Odegard GM,Gates TS,Wise KE,Park C,Siochi E,Constitutive modeling of nanotubereinforced polymer compositesCompos Sci TechnolYear: 2003316711687  
Cho J,Luo JJ,Daniel IM,Mechanical characterization of graphite/epoxy nanocomposites by multiscale analysisCompos Sci TechnolYear: 2007672399240710.1016/j.compscitech.2007.01.006  
Awasthi AP,Lagoudas DC,Hammerand DC,Modeling of graphenepolymer interfacial mechanical behavior using molecular dynamicsModell Simul Mater Sci EngYear: 2009173137  
Montazeri A,Tabar HR,Multiscale modeling of graphene and nanotubebased reinforced polymer nanocompositesPhys Lett AYear: 20113754034404010.1016/j.physleta.2011.08.073  
Ru CQ,Effect of Van der Waals forces on axial buckling of doublewalled carbon nanotubeJ Appl PhysYear: 2000877227723110.1063/1.372973  
Pradhan SC,Buckling of single layer graphene sheet based on nonlocal elasticity and higher order shear deformation theoryPhys Lett AYear: 20093734182418810.1016/j.physleta.2009.09.021  
SakhaeePour A,Elastic buckling of singlelayered graphene sheetComput Mater SciYear: 20094526627010.1016/j.commatsci.2008.09.024  
Rafiee MA,Rafiee J,Yu ZZ,Koratkar N,Buckling resistant graphene nanocompositesAppl Phys LettYear: 20099522310310.1063/1.3269637  
Li C,Chou TW,A structural mechanics approach for the analysis of carbon nanotubesInt J Solids StructYear: 20034024842499  
Li C,Chou TW,Multiscale modeling of compressive behavior of carbon nanotube/polymer compositesCompos Sci TechnolYear: 2006662409241410.1016/j.compscitech.2006.01.013  
Waszczyszyn Z,Cichon C,Radwanska M,Stability of Structures by Finite Element MethodsYear: 1994Amsterdam: Elsevier 
Figures
Article Categories:
Keywords: Multiscale model, Nanocomposite, Buckling, Finite element method. 
Previous Document: Clinical Evaluation Of A Dioic AcidBased Formulation On Facial Skin In An Indian Population.
Next Document: Nanostructuring molecular materials as particles and vesicles for drug delivery, using compressed an...