Research ArticleOn the Effect of Un

Research Article
On the Effect of Unit-Cell Parameters in Predicting the Elastic
Response of Wood-Plastic Composites
Fatemeh Alavi,1 Amir Hossein Behravesh,1 Abbas S.Milani,2 and Davoud Karimi1
1 Department of Mechanical Engineering, Tarbiat Modares University, Tehran 14115-111, Iran
2 School of Engineering, University of British Columbia, Kelowna, BC, Canada V1V 1V7
Correspondence should be addressed to Amir Hossein Behravesh; [e-mail]
Received 27 November 2012; Accepted 8 February 2013
Academic Editor: Toshio Hattori
Copyright © 2013 Fatemeh Alavi et al.This is an open access article distributed under the Creative Commons Attribution License,
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper presents a study on the effect of unit-cell geometrical parameters in predicting elastic properties of a typical wood
plastic composite (WPC). The ultimate goal was obtaining the optimal values of representative volume element (RVE) parameters
to accurately predict the mechanical behavior of the WPC. For each unit cell, defined by a given combination of the above
geometrical parameters, finite element simulation in ABAQUS was carried out, and the corresponding stress-strain curve was
obtained. A uniaxial test according to ASTM D638-02a type V was performed on the composite specimen. Modulus of elasticity
was determined using hyperbolic tangent function, and the results were compared to the sets of finite element analyses.Main effects
of RVE parameters and their interactions were demonstrated and discussed, specially regarding the inclusion of two adjacent wood
particles within one unit cell of the material. Regression analysis was performed tomathematicallymodel the RVE parameter effects
and their interactions over themodulus of elasticity response.Themodel was finally employed in an optimization analysis to arrive
at an optimal set of RVE parameters that minimizes the difference between the predicted and experimental moduli of elasticity.
1. Introduction
It is well known that particles arrangement in a matrix
affects the local stress/strain field in the ensuing composite,
and in turn influences the macrolevel behavior of the
material. Fortunately, effectivemechanical properties of some
heterogeneous materials rely on the average response of
their microstructures and properties of their individual
constituents [1]; hence, basicmicromechanics theories can be
sufficiently used in analyzing such materials.The connection
between the macro- and mesolevel studies is traditionally
viewed via the concept of the representative volume element
(RVE). However, today little quantitative knowledge is
available about minimum RVE sizes for various engineering
materials. Several attempts have been made to determine the
optimum size of an RVE [2]. Numerical-statistical, analytical
approaches, and experimental observations are three methods
applied by researchers to determine this size.
In the numerical-statistical approach, multiple realizations,
finite element simulations of materials unit cells,
and appropriate statistical procedures are employed. Kanit
et al. [3] proposed a quantitative definition of RVE size
which was associated with a given precision in estimating
the desired overall properties and the number of realizations
for a given volume of microstructure. Eventually, they
demonstrated how a minimal volume size for the computation
of effective properties can be determined depending
on the chosen precision and number of realizations.
Trias et al. [4] analyzed the random distribution of fibers
by means of optical microscopy, and the obtained images
were used to generate realizations of statistical representative
volume elements (SRVEs) at microscale. They solved
finite element models of real-microstructure SRVEs with
arbitrary boundary conditions to obtain probability density
functions of stress, strain, and dilatational energy density
and then related the results to stress tensor at any point in
the macroscale by means of a two-scale approach. Other
numerical-statistical approaches based on setting a tolerance
for the scatter in results are given by Vinogradov [5].
Monte-Carlo simulation is a related approach which was
successfully applied by Ostoja-Starzewski [6] and Gusev [7]
2 Journal of Engineering
to predict the overall elastic constants of the studied periodic
composite.
Among analytical approaches, Drugan and Willis [8]
employed an explicit nonlocal constitutive equation by considering
averaged strain fields varying with the position of
material points. A micromechanics model was also proposed
by Jiang et al. [1] for studying the effective elasticmodulus of
composites containing regularly distributed sphere particles.
Three typical particle arrangements in the forms of simple
cubic lattice, body-centered cubic lattice, and face-centered
cubic lattice were investigated. In the case of irregular-shaped
particles, Li and Wongsto [9] derived new unit cells capable
of dealing with problems involving reinforcing particles of
irregular geometries and local imperfections such as debonding
and microcracks in the matrix. Boundary conditions
for their proposed unit cells were derived from appropriate
considerations of the conditions resulting from translational
symmetry transformations. Giraud et al. [10] investigated
the arbitrarily oriented ellipsoidal inhomogeneities to determine
the macroscopic poroelastic properties of transversely
isotropic geomaterials or rock-like composites. They mainly
dealt with separating the effect of matrix anisotropy and of
inhomogeneities in fiber orientation distribution and shape.
In the RVEs based on experimental observations, experimental
analysis often involves the selection of particular
sample geometries for mechanical testing and subjecting
specimens to image analysis after testing. Graham and Yang
[11] and Shan and Gokhale [12] employed this approach for
an instance. In the case of wood-plastic composites, analytical
macromodels have been employed by Hugot and Cazaurang
[13] to predict theWPC effective properties.
Based on the above review, RVE modeling has been
proven as an efficient approach to represent mechanical properties
of composites including WPCs, as micromechanical
modeling of a whole material structure would be computationally
costly and often infeasible. On the other hand,
the determination of optimum RVE parameters including
its dimensions, particle size, and orientation is of major
concern during such an analysis. Having appropriate RVE
parameters, the analysis can further proceed to predict a
composite’s effective properties.Majority of themicromodels
in the current literature of WPCs encompass only one particle,
which precludes the interaction effect of particles at the
microscale. In this investigation, for the first time to address
this issue, two particles of different geometrical specifications
are included in a WPC RVE. Particles are modeled using
elliptical geometries, the sizes ofwhich together with the RVE
size, location, and the angles of the ellipses center lines are
considered as study parameters. Main and interaction effects
of the parameters on the ensuing macrolevel modulus of
elasticity (MOE) are investigated using finite element analysis.
The Taguchi method of DOE (design of experiments) is
employed to determine effects of RVE parameters. Regression
analysis is performed to mathematically model the MOE as
a function of RVE parameters. Finally, sequential quadratic
programming is employed to optimize the RVE parameters to
bestmatch themodeledMOEwiththe preparedWPCsample
with experimentally determined mechanical properties.
Figure 1: A typicalwood plastic composite sample.

ผลลัพธ์ (ไทย) 1:

บทความวิจัยในผลของพารามิเตอร์หน่วยเซลล์ในการทำนายการยืดหยุ่นของไม้พลาสติกคอมโพสิตFatemeh Alavi, Behravesh Hossein Amir 1 ระเบิด 1 S.Milani,2 และ Davoud Karimi11 ภาควิชาวิศวกรรมเครื่องกล มหาวิทยาลัย Tarbiat Modares เตหราน 14115-111 อิหร่าน2 โรงเรียนวิศวกรรม มหาวิทยาลัยของรัฐบริติชโคลัมเบีย เคโลว์น่า BC แคนาดา V1V 1V7ควรส่งจดหมายถึง Amir Hossein Behravesh [e-mail]ได้รับ 27 2555 พฤศจิกายน ยอมรับ 8 2013 กุมภาพันธ์บรรณาธิการวิชาการ: Toshio Hattoriสงวนลิขสิทธิ์ © 2013 Fatemeh Alavi et alนี่คือบทความเข้าเปิดเผยแพร่ภายใต้ลิขสิทธิ์แสดงความคิดสร้างสรรค์ของคอมมอนส์ซึ่งอนุญาตให้ใช้จำกัด กระจาย และการทำซ้ำในสื่อใด ๆ ให้ถูกต้องมีอ้างงานต้นฉบับเอกสารนี้นำเสนอการศึกษาในลักษณะของหน่วยเซลล์ geometrical พารามิเตอร์ในการทำนายคุณสมบัติยืดหยุ่นของไม้ทั่วไปประกอบพลาสติก (WPC) เป้าหมายสูงสุดได้รับค่าสูงสุดของไดรฟ์ข้อมูลตัวแทนพารามิเตอร์องค์ประกอบ (RVE)เพื่อสามารถทำนายพฤติกรรมเชิงกลของ WPC สำหรับแต่ละเซลล์หน่วย กำหนดไว้ โดยการกำหนดข้างต้นพารามิเตอร์ geometrical ทำออกอัตราไนต์ ABAQUS และโค้งต้องใช้ความเครียดที่เกี่ยวข้องได้ได้รับการ ทดสอบตามมาตรฐาน ASTM D638 02a ชนิด V uniaxial ทำตามตัวอย่างประกอบ โมดูลัสของความยืดหยุ่นถูกกำหนดโดยใช้ฟังก์ชันไฮเพอร์โบลิสัมผัส และผลลัพธ์ได้เมื่อเทียบกับชุดของการวิเคราะห์องค์ประกอบจำกัดผลหลักRVE พารามิเตอร์และการโต้ตอบที่แสดง และกล่าว ถึง เป็นพิเศษเกี่ยวกับการรวมสองติดไม้อนุภาคภายในหนึ่งหน่วยเซลล์ของวัสดุ วิเคราะห์การถดถอยได้ทำ tomathematicallymodel RVE พารามิเตอร์ผลและการโต้ตอบผ่าน themodulus ของการตอบสนองความยืดหยุ่นสุดท้ายถูกจ้างในการวิเคราะห์ปรับ Themodel ถึงในชุด RVE พารามิเตอร์เหมาะสม ที่ช่วยลดความแตกต่างระหว่าง moduli คาดการณ์ และการทดลองของความยืดหยุ่น1. บทนำมันเป็นที่รู้จักเรียงอนุภาคในเมทริกซ์ส่งผลต่อฟิลด์ความเครียด/ต้องใช้เครื่องคอมโพสิตเพราะและจะมีผลต่อการทำงาน macrolevel ของการวัสดุ โชคดี effectivemechanical คุณสมบัติของใช้วัสดุต่างชนิดการตอบสนองเฉลี่ยmicrostructures และคุณสมบัติของแต่ละของพวกเขาของพวกเขาconstituents [1]; ดังนั้น ทฤษฎี basicmicromechanics สามารถพอใช้ในการวิเคราะห์วัสดุดังกล่าวการเชื่อมต่อระหว่างศึกษาแมโครและ mesolevel เป็นประเพณีดูผ่านแนวคิดขององค์ประกอบปริมาณพนักงาน(RVE) อย่างไรก็ตาม ความรู้เชิงปริมาณน้อยวันนี้เป็นมีเกี่ยวกับขนาด RVE ขั้นต่ำวิศวกรรมต่าง ๆวัสดุ ได้ทำความพยายามหลายครั้งเพื่อกำหนดขนาดสูงสุดของการ RVE [2] ตัวเลขสถิติ วิเคราะห์วิธี และสังเกตทดลองอยู่สามวิธีใช้ โดยนักวิจัยจะกำหนดขนาดนี้In the numerical-statistical approach, multiple realizations,finite element simulations of materials unit cells,and appropriate statistical procedures are employed. Kanitet al. [3] proposed a quantitative definition of RVE sizewhich was associated with a given precision in estimatingthe desired overall properties and the number of realizationsfor a given volume of microstructure. Eventually, theydemonstrated how a minimal volume size for the computationof effective properties can be determined dependingon the chosen precision and number of realizations.Trias et al. [4] analyzed the random distribution of fibersby means of optical microscopy, and the obtained imageswere used to generate realizations of statistical representativevolume elements (SRVEs) at microscale. They solvedfinite element models of real-microstructure SRVEs witharbitrary boundary conditions to obtain probability densityfunctions of stress, strain, and dilatational energy densityand then related the results to stress tensor at any point inthe macroscale by means of a two-scale approach. Othernumerical-statistical approaches based on setting a tolerancefor the scatter in results are given by Vinogradov [5].Monte-Carlo simulation is a related approach which wassuccessfully applied by Ostoja-Starzewski [6] and Gusev [7]2 Journal of Engineeringto predict the overall elastic constants of the studied periodiccomposite.ระหว่างวิธีวิเคราะห์ Drugan และวิลส์ [8]จ้างชัดเจน nonlocal ขึ้นสมการ โดยพิจารณาaveraged ต้องใช้เขตข้อมูลที่แตกต่างกันกับตำแหน่งของสถานที่วัสดุ ได้เสนอแบบจำลอง micromechanicsโดยเจียง et al. [1] สำหรับการ elasticmodulus ผลของการศึกษาคอมโพสิตที่ประกอบด้วยการกระจายอย่างสม่ำเสมอทรงกลมอนุภาคสามอนุภาคทั่วไปจัดในรูปแบบเรียบง่ายโครงตาข่ายประกอบลูกบาศก์ ร่างกายแปลกลูกบาศก์โครงตาข่ายประกอบ และ แปลกหน้าโครงตาข่ายประกอบลูกบาศก์ถูกสอบสวน ในกรณีที่รูปทรงไม่สม่ำเสมอเซลล์หน่วยใหม่สามารถได้รับอนุภาค Li และ Wongsto [9]จัดการกับปัญหาที่เกี่ยวข้องกับอนุภาคของภาคเอกชนรูปทรงเรขาคณิตที่ไม่สม่ำเสมอและความไม่สมบูรณ์ในท้องถิ่นเช่น debondingและ microcracks ในเมตริกซ์ เงื่อนไขขอบเขตสำหรับการเสนอหน่วย เซลล์ซึ่งที่เหมาะสมพิจารณาเงื่อนไขที่เกิดจาก translationalสมมาตรแปลง สอบสวน Giraud et al. [10]inhomogeneities ellipsoidal โดยมุ่งเน้นการตรวจสอบคุณสมบัติ macroscopic poroelastic ของ transverselyisotropic geomaterials หรือหินเหมือนคอมโพสิต พวกเขาส่วนใหญ่แยกผล ของ anisotropy เมตริกซ์ และinhomogeneities ไฟเบอร์แนวกระจายและรูปร่างใน RVEs ตามสังเกตทดลอง ทดลองวิเคราะห์มักจะเกี่ยวข้องกับการเลือกเฉพาะตัวอย่างรูปทรงเรขาคณิตการทดสอบทางกลแล้วก็กดชัตเตอร์specimens to image analysis after testing. Graham and Yang[11] and Shan and Gokhale [12] employed this approach foran instance. In the case of wood-plastic composites, analyticalmacromodels have been employed by Hugot and Cazaurang[13] to predict theWPC effective properties.Based on the above review, RVE modeling has beenproven as an efficient approach to represent mechanical propertiesof composites including WPCs, as micromechanicalmodeling of a whole material structure would be computationallycostly and often infeasible. On the other hand,the determination of optimum RVE parameters includingits dimensions, particle size, and orientation is of majorconcern during such an analysis. Having appropriate RVEparameters, the analysis can further proceed to predict acomposite’s effective properties.Majority of themicromodelsin the current literature of WPCs encompass only one particle,which precludes the interaction effect of particles at themicroscale. In this investigation, for the first time to addressthis issue, two particles of different geometrical specificationsare included in a WPC RVE. Particles are modeled usingelliptical geometries, the sizes ofwhich together with the RVEsize, location, and the angles of the ellipses center lines areconsidered as study parameters. Main and interaction effectsof the parameters on the ensuing macrolevel modulus ofelasticity (MOE) are investigated using finite element analysis.The Taguchi method of DOE (design of experiments) isemployed to determine effects of RVE parameters. Regressionanalysis is performed to mathematically model the MOE asa function of RVE parameters. Finally, sequential quadraticprogramming is employed to optimize the RVE parameters tobestmatch themodeledMOEwiththe preparedWPCsamplewith experimentally determined mechanical properties.Figure 1: A typicalwood plastic composite sample.