Finite element method has been widely applied in modeling natural fibers and natural fiber reinforced composites. This paper is a comprehensive review of finite element models of natural fibers and natural fiber reinforced composites, focusing on the micromechanical properties (strength, deformation, failure, and damage), thermal properties (thermal conductivity), and macro shape deformation (stress–strain and fracture). Representative volume element model is the most popular homogenization-based multi-scale constitutive method used in the finite element method to investigate the effect of microstructures on the mechanical and thermal properties of natural fibers and natural fiber reinforced composites. The representative volume element models of natural fibers and natural fiber reinforced composites at various length scales are discussed, including two types of geometrical modeling methods, the computer-based modeling method and the image-based modeling method. Their modeling efficiency and accuracy are also discussed.
XueDMiaoM andHuH.Permeability anisotropy of flax nonwoven mats in vacuum‐assisted resin transfer molding. J Text I2011;
102: 612–620.
5.
LiuKTakagiH andYangZ.Evaluation of transverse thermal conductivity of Manila hemp fiber in solid region using theoretical method and finite element method. Mater Des2011;
32: 4586–4458.
6.
WangHLeiYPWangJSet al.
Theoretical and computational modeling of clustering effect on effective thermal conductivity of cement composites filled with natural hemp fibers. J Compos Mater2015;
50: 1509–1521.
7.
WangHQinQ-H andXiaoY.Special n-sided Voronoi fiber/matrix elements for clustering thermal effect in natural-hemp-fiber-filled cement composites. Int J Heat Mass Transf2016;
92: 228–235.
Dan-MallamYAbdullahMZ andMegat YusoffPSM.Influence of recycling frequency on mechanical and physical properties of kenaf fiber reinforced polyoxymethylene composite. J Nat Fibers2016;
13: 532–546.
10.
Dan-MallamYAbdullahMZ andYusoffPSMM.The effect of hybridization on mechanical properties of woven kenaf fiber reinforced polyoxymethylene composite. Polym Compos2014;
35: 1900–1910.
11.
ChenHMiaoM andDingX.Chemical treatments of bamboo to modify its moisture absorption and adhesion to vinyl ester resin in humid environment. J Compos Mater2011;
45: 1533–1542.
12.
RegazziALégerRCornSet al.
Modeling of hydrothermal aging of short flax fiber reinforced composites. Compos Part A Appl S2016;
90: 559–566.
13.
ThuaultABazinJEveSet al.
Numerical study of the influence of structural and mechanical parameters on the tensile mechanical behaviour of flax fibres. J Ind Text2013;
44: 22–39.
14.
Saavedra FloresEI andHaldarS.Micro–macro mechanical relations in Palmetto wood by numerical homogenisation. Compos Struct2016;
154: 1–10.
15.
KernWTKimWArgentoAet al.
Finite element analysis and microscopy of natural fiber composites containing microcellular voids. Mater Des2016;
106: 285–294.
16.
SunCT andVaidyaRS.Prediction of composite properties from a representative volume element. Compos Sci Technol1996;
56: 171–179.
17.
BaleyC.Analysis of the flax fibres tensile behaviour and analysis of the tensile stiffness increase. Compos Part A Appl Sci2002;
33: 939–948.
18.
NilssonT andGustafssonPJ.Influence of dislocations and plasticity on the tensile behaviour of flax and hemp fibres. Compos Part A Appl Sci2007;
38: 1722–1728.
19.
IlczyszynFCherouatAMontayG.Effect of hemp fibre morphology on the mechanical properties of vegetal fibre composite material. AMR2014;
875–877: 485–489.
20.
KentP.Micromechanical modelling of wood and fibre properties. PhD Thesis, Lund University, Sweden, 2000.
21.
TrivaudeyFPlacetVGuicheret-RetelVet al.
Nonlinear tensile behaviour of elementary hemp fibres. Part II: modelling using an anisotropic viscoelastic constitutive law in a material rotating frame. Compos Part A Appl Sci2015;
68: 346–355.
22.
ThygesenLG andGierlingerN.The molecular structure within dislocations in Cannabis sativa fibres studied by polarised Raman microspectroscopy. J Struct Biol2013;
182: 219–225.
23.
PlacetVCisséO andLamine BoubakarM.Nonlinear tensile behaviour of elementary hemp fibres. Part I: investigation of the possible origins using repeated progressive loading with in situ microscopic observations. Compos Part A Appl Sci2014;
56: 319–327.
24.
Del MastoATrivaudeyFGuicheret-RetelVet al.
Nonlinear tensile behaviour of elementary hemp fibres: a numerical investigation of the relationships between 3D geometry and tensile behaviour. J Mater Sci2017;
52: 6591–6610.
25.
DanielsHE.The statistical theory of the strength of bundles of threads. I. Proc R Soc Lond Soc Math Phys Eng Sci1945;
183: 405–435.
ZwebenC andRosenBW.A statistical theory of material strength with application to composite materials. J Mech Phys Solids1970;
18: 189–206.
28.
HarlowDG andPhoenixSL.The chain-of-bundles probability model for the strength of fibrous materials I: analysis and conjectures. J Compos Mater1978;
12: 195–214.
29.
HarlowDG andPhoenixSL.The chain-of-bundles probability model for the strength of fibrous materials II: a numerical study of convergence. J Compos Mater1978;
12: 314–334.
30.
BeyerleinIJ andLandisCM.Shear-lag model for failure simulations of unidirectional fiber composites including matrix stiffness. Mech Mater1999;
31: 331–350.
BeakouA andCharletK.Mechanical properties of interfaces within a flax bundle – part II: numerical analysis. Int J Adhes Adhes2013;
43: 54–59.
33.
CharletK andBéakouA.Mechanical properties of interfaces within a flax bundle – part I: experimental analysis. Int J Adhes Adhes2011;
31: 875–881.
34.
Buet-GautierK andBoisseP.Experimental analysis and modeling of biaxial mechanical behavior of woven composite reinforcements. Exp Mech2001;
41: 260–269.
35.
BoissePGasserA andHivetG.Analyses of fabric tensile behaviour: determination of the biaxial tension–strain surfaces and their use in forming simulations. Compos Part A Appl Sci2001;
32: 1395–1414.
36.
TephanyCSoulatDGillibertJet al.
Influence of the non-linearity of fabric tensile behavior for preforming modeling of a woven flax fabric. Text Res J2016;
86: 604–617.
37.
NajjarWLegrandXPupinCet al.
A simple discrete method for the simulation of the preforming of woven fabric reinforcement. KEM2012;
504–506: 213–218.
38.
NajjarWLegrandXDal SantoPet al.
Analysis of the blank holder force effect on the preforming process using a simple discrete approach. KEM2013;
554–557: 441–446.
39.
HaldarSGheewalaNGrande-AllenKJet al.
Multi-scale mechanical characterization of Palmetto Wood using digital image correlation to develop a template for biologically-inspired polymer composites. Exp Mech2010;
51: 575–589.
40.
GershonALBruckHAXuSet al.
Multiscale mechanical and structural characterizations of Palmetto wood for bio-inspired hierarchically structured polymer composites. Mater Sci Eng C2010;
30: 235–244.
41.
HaldarS andBruckHA.Characterization of dynamic damage mechanisms in Palmetto wood as biological inspiration for impact resistant polymer composites. Mech Mater2013;
57: 97–108.
42.
HofstetterKHellmichC andEberhardsteinerJ.Development and experimental validation of a continuum micromechanics model for the elasticity of wood. Eur J Mech A Solid2005;
24: 1030–1053.
43.
MishnaevskyL andQingH.Micromechanical modelling of mechanical behaviour and strength of wood: state-of-the-art review. Comp Mater Sci2008;
44: 363–370.
44.
QingH andMishnaevskyL.3D hierarchical computational model of wood as a cellular material with fibril reinforced, heterogeneous multiple layers. Mech Mater2009;
41: 1034–1049.
45.
QingH andMishnaevskyL.3D multiscale micromechanical model of wood: from annual rings to microfibrils. Int J Solids Struct2010;
47: 1253–1267.
46.
MurataKNagaiH andNakanoT.Estimation of width of fracture process zone in spruce wood by radial tensile test. Mech Mater2011;
43: 389–396.
47.
Saavedra FloresEIde Souza NetoEA andPearceC.A large strain computational multi-scale model for the dissipative behaviour of wood cell-wall. Comp Mater Sci2011;
50: 1202–1211.
48.
RitschelFBrunnerAJ andNiemzP.Nondestructive evaluation of damage accumulation in tensile test specimens made from solid wood and layered wood materials. Compos Struct2013;
95: 44–52.
49.
LoTYCuiHZ andLeungHC.The effect of fiber density on strength capacity of bamboo. Mater Lett2004;
58: 2595–2598.
50.
LoTYCuiHZTangPWCet al.
Strength analysis of bamboo by microscopic investigation of bamboo fibre. Constr Build Mater2008;
22: 1532–1535.
51.
DixonPGAhvenainenPAijaziANet al.
Comparison of the structure and flexural properties of Moso, Guadua and Tre Gai bamboo. Constr Build Mater2015;
90: 11–17.
52.
MaJChenWZhaoLet al.
Elastic buckling of bionic cylindrical shells based on bamboo. J Bionic Eng2008;
5: 231–238.
53.
ZhaoLMaJWangTet al.
Lightweight design of mechanical structures based on structural bionic methodology. J Bionic Eng2010;
7: S224–S231.
54.
GibsonLJ.The hierarchical structure and mechanics of plant materials. J R Soc Interface2012;
9: 2749–2766.
55.
EstradaMLineroDL andRamírezF.Constitutive relationship of the fiber cluster of bamboo Guadua angustifolia, determined by means of a Weibull probability function and a model of progressive failure. Mech Mater2013;
63: 12–20.
56.
ZouMXuSWeiCet al.
A bionic method for the crashworthiness design of thin-walled structures inspired by bamboo. Thin Wall Struct2016;
101: 222–230.
57.
JiaoHZhangY andChenW.The lightweight design of low RCS pylon based on structural bionics. J Bionic Eng2010;
7: 182–190.
58.
SilvaECNWaltersMC andPaulinoGH.Modeling bamboo as a functionally graded material: lessons for the analysis of affordable materials. J Mater Sci2006;
41: 6991–7004.
59.
ChandNShuklaM andSharmaMK.Analysis of mechanical behaviour of bamboo (dendrocalamus strictus) by using FEM. J Nat Fibers2008;
5: 127–137.
60.
WSFuZZHanWZhouJBResearch on finite element model for parallel to bamboo culms axial shear. AMM2013;
477–478: 986–989.
61.
PalombiniFLKindleinWde OliveiraBFet al.
Bionics and design: 3D microstructural characterization and numerical analysis of bamboo based on X-ray microtomography. Mater Charact2016;
120: 357–336.
62.
GuH andMiaoM.Optimising fibre alignment in twisted yarns for natural fibre composites. J Compos Mater2013;
48: 2993–3002.
MuralidharK.Equivalent conductivity of a heterogeneous medium. Int J Heat Mass Tranasf1990;
33: 1759–1766.
65.
MingqingZBomingY andDuanmingZ.An analytical solution for transverse thermal conductivities of unidirectional fibre composites with thermal barrier. J Phys D Appl Phys2002;
35: 1867–1873.
66.
WangMHeJYuJet al.
Lattice Boltzmann modeling of the effective thermal conductivity for fibrous materials. Int J Therm Sci2007;
46: 848–855.
67.
ZhengG.Numerical investigation of characteristic of anisotropic thermal conductivity of natural fiber bundle with numbered lumens. Math Probl Eng2014;
2014: 1–8.
68.
WangHXiaoY andQinQ.2D hierarchical heat transfer computational model of natural fiber bundle reinforced composite. Sci Iran2016;
23: 268–276.
69.
IslamMR andPramilaA.Thermal conductivity of fiber reinforced composites by the FEM. J Compos Mater1999;
33: 1699–1715.
70.
HasselmanDPH andJohnsonLF.Effective thermal conductivity of composites with interfacial thermal barrier resistance. J Compos Mater1987;
21: 508–515.
71.
WangHZhaoX andWangJ.Interaction analysis of multiple coated fibers in cement composites by special n-sided interphase/fiber elements. Compos Sci Technol2015;
118: 117–126.
72.
DevireddySBR andBiswasS.Physical and thermal properties of unidirectional banana-jute hybrid fiber-reinforced epoxy composites. J Reinf Plast Compos2016;
35: 1157–1172.
73.
Tucker IiiCL andLiangE.Stiffness predictions for unidirectional short-fiber composites: review and evaluation. Compos Sci Technol1999;
59: 655–671.
74.
ZhangDMilanovicNRZhangYet al.
Effects of humidity conditions at fabrication on the interfacial shear strength of flax/unsaturated polyester composites. Compos Part B Eng2014;
60: 186–192.
75.
SliserisJYanL andKasalB.Numerical modelling of flax short fibre reinforced and flax fibre fabric reinforced polymer composites. Compos Part B Eng2016;
89: 143–154.
76.
ModniksJ andAndersonsJ.Modeling elastic properties of short flax fiber-reinforced composites by orientation averaging. Comp Mater Sci2010;
50: 595–599.
77.
ModniksJ andAndersonsJ.Modeling the non-linear deformation of a short-flax-fiber-reinforced polymer composite by orientation averaging. Compos Part B Eng2013;
54: 188–193.
78.
KossaA andSzabóL.Exact integration of the von Mises elastoplasticity model with combined linear isotropic-kinematic hardening. Int J Plast2009;
25: 1083–1106.
79.
BourmaudAAusiasGLebrunGet al.
Observation of the structure of a composite polypropylene/flax and damage mechanisms under stress. Ind Crop Prod2013;
43: 225–236.
80.
DrayDGilorminiP andRegnierG.Comparison of several closure approximations for evaluating the thermoelastic properties of an injection molded short-fiber composite. Compos Sci Technol2007;
67: 1601–1610.
ZhongYTranLQNKureemunUet al.
Prediction of the mechanical behavior of flax polypropylene composites based on multi-scale finite element analysis. J Mater Sci2017;
52: 4957–4967.
83.
LubineauG andLadevèzeP.Construction of a micromechanics-based intralaminar mesomodel, and illustrations in ABAQUS/Standard. Comp Mater Sci2008;
43: 137–145.
84.
LapczykI andHurtadoJA.Progressive damage modeling in fiber-reinforced materials. Compos Part A Appl Sci2007;
38: 2333–2341.
85.
LiuPFGuZPYangYHet al.
A nonlocal finite element model for progressive failure analysis of composite laminates. Compos Part B Eng2016;
86: 178–196.
86.
ReisPNBFerreiraJAMAntunesFVet al.
Flexural behaviour of hybrid laminated composites. Compos Part A Appl Sci2007;
38: 1612–1620.
87.
AssararMZouariWSabhiHet al.
Evaluation of the damping of hybrid carbon–flax reinforced composites. Compos Struct2015;
132: 148–154.
88.
YanLChouwN andYuanX.Improving the mechanical properties of natural fibre fabric reinforced epoxy composites by alkali treatment. J Reinf Plast Compos2012;
31: 425–437.
89.
GhasemnejadHSoroushVRMasonPJet al.
To improve impact damage response of single and multi-delaminated FRP composites using natural Flax yarn. Mater Des 1980–20152012;
36: 865–873.
90.
BrighentiRCarpinteriA andScorzaD.Mechanics of interface debonding in fibre-reinforced materials. J Compos Mater2015;
50: 2699–2718.
91.
XuZLiJWuXet al.
Effect of kidney-type and circular cross sections on carbon fiber surface and composite interface. Compos Part A Appl Sci2008;
39: 301–307.
92.
BarbieriGBiolziLBocciarelliMet al.
Size and shape effect in the pull-out of FRP reinforcement from concrete. Compos Struct2016;
143: 395–417.
93.
JägerJSauseMGRBurkertFet al.
Influence of plastic deformation on single-fiber push-out tests of carbon fiber reinforced epoxy resin. Compos Part A Appl Sci2015;
71: 157–167.
94.
ChenHMiaoM andDingX.Influence of moisture absorption on the interfacial strength of bamboo/vinyl ester composites. Compos Part A Appl Sci2009;
40: 2013–2019.
95.
MokhtarKAlexandreVAlainBet al.
Damage analysis of composites reinforced with Alfa fibers: viscoelastic behavior and debonding at the fiber/matrix interface. J Appl Polym Sci2016; 133.
96.
KeckS andFullandM.Effect of fibre volume fraction and fibre direction on crack paths in flax fibre-reinforced composites. Eng Fract Mech2016;
167: 201–209.
97.
WangWLoweADaveySet al.
Establishing a new Forming Limit Curve for a flax fibre reinforced polypropylene composite through stretch forming experiments. Compos Part A Appl Sci2015;
77: 114–123.
98.
PanamoottilSMDasR andJayaramanK.Anisotropic continuum damage model for prediction of failure in flax/polypropylene fabric composites. Polym Compos2016;
37: 2588–2597.
99.
VoyiadjisGZ andVensonAR.Experimental damage investigation of a SiC-Ti aluminide metal matrix composite. Int J Damage Mech1995;
4: 338–361.
100.
Jacob FishKS.Multiscale analysis of composite materials and structures. Compos Sci Technol2000;
60: 2547–2556.
101.
Heydari-MeybodiMSaber-SamandariSSadighiM.A new approach for prediction of elastic modulus of polymer/nanoclay composites by considering interfacial debonding: experimental and numerical investigations. Compos Sci Technol2015;
117: 379–385.
102.
ZhiCLongH andMiaoM.Influence of microbond test parameters on interfacial shear strength of fiber reinforced polymer-matrix composites. Compos Part A Appl Sci2017;
100: 55–63.
103.
Guillebaud-BonnafousCVasconcellosDTouchardFet al.
Experimental and numerical investigation of the interface between epoxy matrix and hemp yarn. Compos Part A Appl Sci2012;
43: 2046–2058.
104.
VaughanTJ andMcCarthyCT.Micromechanical modelling of the transverse damage behaviour in fibre reinforced composites. Compos Sci Technol2011;
71: 388–396.
105.
YangLYanYLiuYet al.
Microscopic failure mechanisms of fiber-reinforced polymer composites under transverse tension and compression. Compos Sci Technol2012;
72: 1818–1825.
106.
ShokrianMDShelesh-NezhadK andSoudmandBH.3D FE analysis of tensile behavior for co-PP/SGF composite by considering interfacial debonding using CZM. J Reinf Plast Compos2015;
35: 365–374.
107.
AlamPFagerlundPHägerstrandPet al.
L-Lysine templated CaCO3 precipitated to flax develops flowery crystal structures that improve the mechanical properties of natural fibre reinforced composites. Compos Part A Appl Sci2015;
75: 84–88.
108.
PerrierAEcaultRTouchardFet al.
Towards the development of laser shock test for mechanical characterisation of fibre/matrix interface in eco-composites. Polym Test2015;
44: 125–134.
109.
HbibMGuessasmaSBassirDet al.
Interfacial damage in biopolymer composites reinforced using hemp fibres: finite element simulation and experimental investigation. Compos Sci Technol2011;
71: 1419–1426.
110.
ZhandarovS andMäderE.An alternative method of determining the local interfacial shear strength from force–displacement curves in the pull-out and microbond tests. Int J Adhes Adhes2014;
55: 37–42.
111.
ScarponiC.Hemp fiber composites for the design of a Naca cowling for ultra-light aviation. Compos Part B Eng2015;
81: 53–63.
112.
ScarponiC andMessanoM.Comparative evaluation between E-Glass and hemp fiber composites application in rotorcraft interiors. Compos Part B Eng2015;
69: 542–549.
113.
KongCLeeH andParkH.Design and manufacturing of automobile hood using natural composite structure. Compos Part B Eng2016;
91: 18–26.
