Next Article in Journal
Heterogeneous Computing (CPU–GPU) for Pollution Dispersion in an Urban Environment
Previous Article in Journal
An Efficient Two-Layer Non-Hydrostatic Model for Investigating Wave Run-Up Phenomena
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Mesh Sensitivity and Cohesive Properties on Simulation of Typha Fiber/Epoxy Microbond Test

1
Department of Mechanical Engineering, Universitas Syiah Kuala, Darussalam, Banda Aceh 23111, Indonesia
2
School of Industrial Technology, Universiti Sains Malaysia, Penang 11800, Malaysia
*
Authors to whom correspondence should be addressed.
Submission received: 26 November 2019 / Revised: 25 December 2019 / Accepted: 27 December 2019 / Published: 3 January 2020

Abstract

:
The microbond test for natural fibers is difficult to conduct experimentally due to several challenges including controlling the gap distance of the blade, the meniscus shape, and the large data spread. In this study, a finite element simulation was performed to investigate the effects of the bonding characteristics in the interface between the fiber and matrix on the Typha fiber/epoxy microbond test. Our aim was to obtain the accurate mesh and cohesive properties via simulation of the Typha fiber/epoxy microbond test using the cohesive zone model technique. The axisymmetric model was generated to model the microbond test specimen with a cohesive layer between the fiber and matrix. The cohesive parameter and mesh type were varied to determine the appropriate cohesive properties and mesh type. The fine mesh with 61,016 elements and cohesive properties including stiffness coefficients Knn = 2700 N/mm3, Ktt = 2700 N/mm3, and Kss = 2700 N/mm3; fracture energy of 15.15 N/mm; and damage initiation tnn = 270 N/mm2, ttt = 270 N/mm2, and tss = 270 N/mm2 were the most suitable. The cohesive zone model can describe the debonding process in the simulation of the Typha fiber/epoxy microbond test. Therefore, the results of the Typha fiber/epoxy microbond simulation can be used in the simulation of Typha fiber reinforced composites at the macro-scale.

Graphical Abstract

1. Introduction

Natural fiber-reinforced polymer composites have attracted research interest as replacements for synthetic fibers. Synthetic fibers are not environmentally friendly due to being nonbiodegradable, which causes environmental problems. Therefore, natural fibers have rapidly been developed to reinforce polymer composites. Natural fibers are biodegradable and renewable. Several natural fiber composites have been studied including ramie, kenaf, coir, sisal, hemp, and jute [1,2,3]. One of the natural fibers that has the potential to reinforce composites is Typha fiber. Typha is widely available in most countries [4]. Typha grows rapidly on wetlands and is often considered as parasite. Although Typha is abundantly present in nature, its potential is still undeveloped compared to other natural fibers [5]. Typha fiber composites are known to have many advantages, such as good flexibility, light weight, good mechanical strength, low density, and renewability [6,7]. However, several challenges exist to the use of natural fiber-reinforced polymer composite, one of which is the low compatibility between natural fiber and the matrix polymer [8], which causes interlayer gaps and initial cracks along the fiber. The interfacial composite is the contact area between the matrix and the reinforcing fiber that carries the load received by the matrix to the fiber. The structures of the fiber and the matrix considerably influence the bonding of the fiber–matrix interface, which also affects the characteristics of the composite structure [9,10,11].
Several testing methods, such as the fiber push-out test, microbond test, and fragmentation test, have been used to evaluate the shear strength of the interface of the composite. The microbond test was developed by Miller et al. [12] and is one of the most widely used test methods for testing fiber–matrix interfacial shear strength (IFSS).
Microbond tests have been commonly used to study the effects of surface treatments on fibers and the interface behavior between matrices and fibers. The microbond test for natural fibers is hard to conduct experimentally due the several challenges such as difficulty in controlling the gap distance of the blade, the meniscus shape, the plastic yield of the matrix, and the reproducibility of sample preparation [13]. Microbond test results are reported cause large variations in test data and the large data spread becomes an obstacle for determining the shear strength of composite interfaces compared to other tests [14]. Therefore, microbond testing with numerical simulation approaches has attracted researchers. Numerical simulation using the finite element method has been used to simplify the investigation of the bonding characteristics on the interface between the fiber and matrix [15,16].
Several reports have been published about microbond test simulations. Ash et al. developed an axisymmetric model to study microbond tests including realistic droplet shapes on the glass fiber–polymer matrix [17]. Kang et al. modeled microbond tests using two-dimensional (2D) axisymmetric models and considered different droplet shapes, such as cylinders and spheres, on the effect of meniscus on carbon–epoxy fibers [14]. Crack propagation during fiber-on-fiber/matrix push-out was modeled using a cohesive zone approach in the volumetric finite element method (FEM) with an axisymmetric model [18]. The cohesive zone model was used to simulate the interphase between fiber and matrix, providing an ideal tool to study the interfacial debonding progress in composites [19]. However, no study has investigated the behavior of interfacial bonds in the microbond test of Typha fiber/epoxy. The aim of this study was to determine the accurate mesh and cohesive properties on the simulation of the Typha fiber/epoxy microbond test using the cohesive zone model technique.
The Typha fiber/epoxy microbond test was simulated using the FEM. The mesh sensitivity and cohesive properties are two parameters that play important roles in microbond test simulation. Therefore, we wanted to determine the appropriate mesh type and cohesive properties to obtain good agreement between simulation and experimental results. We wanted to analyze the mesh sensitivity and cohesive properties of the Typha fiber/epoxy microbond test as these have not been reported in the literature.

2. Materials and Methods

2.1. FEM

The FEM and geometry used were similar to the microbond test experiment, as shown in Figure 1. The commercial finite element software Abaqus 6.14 (Dassault Systèmes Simulia Corp., Johnston, RI, USA) was used to generate the model and analyze the Typha fiber/epoxy microbond test. The axisymmetric model was generated to model the microbond test specimen with a cohesive layer between the fiber and matrix with 0.01 mm thickness, as shown in Figure 2, where the embedded length, fiber radius, and droplet height were 1.67, 0.125, and 0.3 mm, respectively. The vise blade was 0.78 mm away from the interface. The vise blade was modelled as a rigid body. The hard contact type was selected as the interaction between the blade and resin. Tie constraints were used to tie fibers to the cohesive layer and matrix. A displacement load of 0.466 mm was applied at the end of the fiber. The magnitude of displacement was obtained from experiments. The edge of the fiber and matrix provided the asymmetrical boundary condition, where only the Y plane axis is free (U1 = U2 = UR3 = 0), while the vise blade is in the fully built-in fixed support condition (U1 = U2 = UR1 = UR2 = UR3 = 0). Finite element analysis was performed on an ASUS computer (AsusTek Computer Inc., Taipei, Taiwan) with an Intel core i9 processor (Intel Corp., Santa Clara, CA, USA) and 16 GB RAM.
The mesh was varied to determine the mesh that produces results near to the experimental microbond test results. The element near the interfacial area was adjusted. The quadrilateral linear CAX4 type and linear triangular CAX3 type were the elements used on the epoxy resin and Typha fiber, whereas the elements on the cohesive layer were quadrilateral linear COHAX4 type. The mesh on the Typha fiber/epoxy microbond test model is displayed in Figure 3. The numbers of each element in the microbond testing model are presented in Table 1.

2.2. Material Properties

Dugdale [20] and Barenblatt [21] introduced the cohesive zone model in 1960 and 1962, respectively, when solving the problem of cracking of perfect brittle material. Three modes of fractures, opening mode, sliding mode, and tear mode, can appear at the interface between the fiber and matrix when the loading conditions are applied to the composite. The cohesive zone model specifies three bond conditions for each fracture mode: elastic bond, debonding, and debonded conditions. All the fracture modes can be characterized by a traction–separation (tδ) curve as shown in Figure 4. The condition of the elastic bond is determined by the interfacial stiffness (K), which is the elastic constitutive matrix that relates t and δ, as shown in Equation (1).
t = K δ { t n t s t t } = [ K n n K n s K n t K n s K s s K s t K n t K s t K t t ] { δ n δ s δ t }
where n is normal mode, s is shear mode, and t is tear mode.
The material properties are listed in Table 2. Typha fiber properties were obtained from our previous research [22] and the epoxy resin properties were obtained from Selmi [23].
Adjusting the stiffness parameters and damage criteria for cohesive elements is convenient with finite element software [15]. Therefore, the cohesive properties were varied to obtain agreement of the simulation results with the experimental Typha fiber/epoxy microbond results. Table 3 lists the various cohesive properties. These values, obtained based on several simulations, were accepted as reported using the same method of Dadej [24] and Potukuchi [25].

3. Mesh Sensitivity Analysis

The microbond test has several parameters that must be considered, including fiber diameter, embedded length, meniscus angle, contact angle, and the distance between the vise blade and the specimen. Cohesive parameters must be determined first to produce suitable cohesive properties. When accomplished, the experimental results would be similar to the experimental results [15]. The mesh has an important role in finite element analysis for determining the accuracy of the simulation results to ensure good agreement with the experimental results.
The force–displacement curve was obtained from the mesh sensitivity study in the microbond test model, as shown in Figure 5. A, B, C, and D represent the coarse, medium, fine, and very fine mesh, respectively, corresponding to P, Q, R, and S, respectively, which represent computational time. We found that the higher the number of the elements, the lower the reaction force; this indicated that mesh selection affected the simulation result. Figure 5 shows that the coarse mesh had the highest force, which was 9.7 N; the lowest force on the very fine mesh was 1.5 N. The maximum forces for the medium and fine meshes were 7.5 and 2.9 N, respectively.
The results of testing several cohesive properties on the coarse, medium, fine, and very fine mesh are shown in Figure 6a–d, respectively. The force–displacement curves for the cohesive properties were compared with the experimental results. The experimental curve decreased dramatically after reaching the maximum force, with a maximum force value of 2.6 N. The maximum force in the microbond simulation test with fine mesh had a maximum force resembling the experimental results, as shown in Figure 5.
In the simulation, the damage initiation parameter value was varied from 200 to 270 N/mm2. The variation of the cohesive properties with coarse mesh showed that the maximum force was around 5–10 N. This maximum force exceeded the value obtained from the experiment. For medium mesh, the maximum force was 4–7 N. Based on the image, VR3, VR4, and VR6 show that the curve decreased after reaching maximum force. The simulation of the Typha fiber/epoxy microbond test with various cohesive properties on the fine mesh is shown in Figure 6c. We found that the curve was similar to the experimental curve; almost all curves had a maximum force close to the experimental value except for VR6 and VR7. The VR2 curve was the most similar to the experimental results. Figure 6d depicts the results of various cohesive properties on very fine mesh, showing that the maximum force is below the experimental result. Only the VR4 curve decreased after reaching the maximum force.
Based on finite element theory, element size is a critical aspect related to the result accuracy. The smaller the element size, the more accurate the analysis results. However, in certain cases, a smaller element size produces a higher calculation value compared to actual value obtained from the experiments [26]. In this situation, a limit must be set on the percentage error of the results obtained from calculations to determine the appropriate element size [27]. In our study, the least error occurred with the fine mesh.
The incorrect selection of cohesive properties in the microbond test simulation leads to results incompatible with the experimental values. Dadej and Sadowski analyzed the response of the cohesive zone model parameters on glass/epoxy composite and found that incorrectly choosing the damage initiation parameter causes an unstable growth of crack length [24]. Schellekens and de Borst [28] reported that excessive stiffness coefficient parameter values produce unstable analysis, which affects the numerical computation time. The selection of a suitable mesh is also important. As reported by Soto [29], incorrect selection of the number of elements generates incorrect crack initiation and propagation, causing poor agreement of simulation results with experimental results.
From this result, the interfacial shear strength (IFSS) of Typha fiber/epoxy was obtained using Equation (2) [30].
I F S S = F m a x π d f L e m
where Fmax is the maximum force when debonding occurs between the fiber and matrix, df is fiber diameter, and Lem is embedded length.
The interfacial shear strength obtained from the simulation results of microbond testing with fine mesh was 4.4 MPa; the interfacial shear strength from the experimental results is 3.2 MPa per previous research [22]. The sensitivity mesh study revealed that fine mesh with 61,016 elements is the most suitable for Typha fiber/epoxy microbond test simulation. Therefore, we concluded that the cohesive properties and the type of mesh used in this simulation produced results close to the experimental results.

4. Interfacial Debonding

Figure 7 displays the microbond test simulation results, which show the force and displacement of Typha fiber/epoxy. The force–displacement curve illustrates that the matrix is still attached to the fiber until a maximum force of 2.9 N at 0.31 mm displacement. After reaching the maximum force, the line on the curve decreases sharply, which reveals debonding has occurred. As illustrated in Figure 8, the debonding phenomenon is similar to the results of the simulation conducted by Hao et al. [31] for a carbon fiber/bismaleimide composites microbond test, but with a different magnitude.
Figure 8 depicts a contour plot of shear stress in a simulation of Typha fiber/epoxy composite. When the vise blade hit the matrix first, we observed that the stress at the resin tip is concentrated at the vise blade tip. The shear stress contour spreads from the epoxy resin to the cohesive layer, which was modeled as the interfacial region between the fiber and matrix. Then, the stress was distributed to the fiber, as shown in Figure 8a. In the next step, more shear stress appeared on the fiber and continued to propagate down through the interface area, as shown in Figure 8b. The moment immediately before debonding, depicted in Figure 8c, illustrates the last step, which is the shear stress concentration moving into the bottom of the test specimen before full debonding occurred.
Figure 9 illustrates the maximum traction damage initiation criterion (MAXSCRT), which is a scalar parameter where total failure is equal to one. The blue and red contours indicate lower and higher MAXSCRT values, respectively. At a step time value of 0.5, red and blue contours were valued at 0.697 and 3 × 10−5, respectively, which means failure at the upper cohesive layer has not started. When the step time was 0.689, the upper cohesive layer started to fail due to the value of the red contour reaching 1, while the blue contour was 5 × 10−5. The cohesive layer totally failed after a step time of 0.7, as indicated by both the value of red and blue contours reaching one.
We found that the appropriate mesh type and cohesive properties for the simulation of the Typha fiber/epoxy microbond test. The debonding process of the interfacial area between the fiber and the matrix can be clearly described by the proposed model. The mesh type and cohesive properties that matched the experimental results can facilitate further work on understanding the interfacial behavior of Typha fiber and epoxy at the macro scale.

5. Conclusions

In this study, we conducted Typha fiber/epoxy microbond testing simulations with several different cohesive properties and mesh types. From the results, we concluded:
  • Fine mesh with 61,016 elements is the most suitable type for the Typha fiber/epoxy microbond simulation test based on the comparison of simulation and experimental results.
  • The simulation results agreed well with experimental results of load and displacement curves and interfacial shear strength values when cohesive properties including stiffness coefficients (Knn of 2700 N/mm3, Ktt of 2700 N/mm3, and Kss of 2700 N/mm3), fracture energy of 15.15 N/mm, and damage initiation (tnn = 270 N/mm2, ttt = 270 N/mm2, tss = 270 N/mm2), with fine mesh.
  • The cohesive zone model can describe the debonding process during the simulation of the Typha fiber/epoxy microbond test.
From this research, we expect that future Typha fiber simulation processes on a macro-scale will be able to directly use the mesh parameters and cohesive properties reported here. In future work, the effect of the blade gap distance on the Typha fiber/epoxy microbond test will be tested.

Author Contributions

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

Funding

This work was supported by Ministry Research, Technology and Higher Education of Republic of Indonesia by Program Magister Menuju Doktor untuk Sarjana Unggul (PMDSU) scheme No. 72/UN11.2/PP/SP3/2018.

Acknowledgments

The authors gratefully thankful to the World Class Professor (WCP) Program 2018 and 2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Orue, A.; Jauregi, A.; Unsuain, U.; Labidi, J.; Eceiza, A.; Arbelaiz, A. The effect of alkaline and silane treatments on mechanical properties and breakage of sisal fibers and poly (lactic acid)/sisal fiber composites. Compos. Part A Appl. Sci. Manuf. 2016, 84, 186–195. [Google Scholar] [CrossRef]
  2. Caprino, G.; Carrino, L.; Durante, M.; Langella, A.; Lopresto, V. Low impact behaviour of hemp fibre reinforced epoxy composites. Compos. Struct. 2015, 133, 892–901. [Google Scholar] [CrossRef]
  3. Faruk, O.; Bledzki, A.K.; Fink, H.-P.; Sain, M. Biocomposites reinforced with natural fibers: 2000–2010. Prog. Polym. Sci. 2012, 37, 1552–1596. [Google Scholar] [CrossRef]
  4. Zhou, B.; Tu, T.; Kong, F.; Wen, J.; Xu, X. Revised phylogeny and historical biogeography of the cosmopolitan aquatic plant genus Typha (Typhaceae). Sci. Rep. 2018, 8, 8813. [Google Scholar] [CrossRef]
  5. Ramanaiah, K.; Ratna Prasad, A.V.; Chandra Reddy, K.H. Mechanical Properties and Thermal Conductivity of Typha angustifolia Natural Fiber–Reinforced Polyester Composites. Int. J. Polym. Anal. Charact. 2011, 16, 496–503. [Google Scholar] [CrossRef]
  6. Bajwa, D.S.; Sitz, E.D.; Bajwa, S.G.; Barnick, A.R. Evaluation of cattail (Typha spp.) for manufacturing composite panels. Ind. Crops Prod. 2015, 75, 195–199. [Google Scholar] [CrossRef]
  7. Wuzella, G.; Mahendran, A.R.; Bätge, T.; Jury, S.; Kandelbauer, A. Novel, binder-free fiber reinforced composites based on a renewable resource from the reed-like plant Typha sp. Ind. Crops Prod. 2011, 33, 683–689. [Google Scholar] [CrossRef]
  8. Khalil, H.P.S.A.; Ismail, H.; Rozman, H.D.; Ahmad, M.N. The effect of acetylation on interfacial shear strength between plant fibres and various matrices. Eur. Polym. J. 2001, 37, 1037–1045. [Google Scholar] [CrossRef]
  9. Chen, J.; Hojjati, M. Microdielectric Analysis and Curing Kinetics of an Epoxy Resin System. Polym. Eng. Sci. 2007, 47, 150–158. [Google Scholar] [CrossRef]
  10. Ratnam, M.M.; Lim, J.H.; Khalil, H.P.S.A. Study of three-dimensional deformation of a pallet using phase-shiff shadow moiré and finite-element analysis. Exp. Mech. 2005, 45, 9–17. [Google Scholar] [CrossRef]
  11. Elmatzoglou, M.; Avdelas, A. Numerical Modelling of Double-Steel Plate Composite Shear Walls. Computation 2017, 5, 12. [Google Scholar] [CrossRef] [Green Version]
  12. Miller, B.; Muri, P.; Rebenfeld, L. A microbond method for determination of the shear strength of a fiber/resin interface. Compos. Sci. Technol. 1987, 28, 17–32. [Google Scholar] [CrossRef]
  13. Sockalingam, S.; Nilakantan, G. Fiber-matrix interface characterization through the microbond test: A review. Int. J. Aeronaut. Space Sci. 2012, 13, 282–295. [Google Scholar] [CrossRef] [Green Version]
  14. Kang, S.K.; Lee, D.B.; Choi, N.S. Fiber/epoxy interfacial shear strength measured by the microdroplet test. Compos. Sci. Technol. 2009, 69, 245–251. [Google Scholar] [CrossRef]
  15. Xiong, X.; Shen, S.Z.; Hua, L.; Liu, J.Z.; Li, X.; Wan, X.; Miao, M. Finite element models of natural fibers and their composites: A review. J. Reinf. Plast. Compos. 2018, 37, 617–635. [Google Scholar] [CrossRef]
  16. Huzni, S.; Ilfan, M.; Sulaiman, T.; Fonna, S.; Ridha, M.; Arifin, A.K. Finite element modeling of delamination process on composite laminate using cohesive elements. Int. J. Automot. Mech. Eng. 2013, 7, 1023. [Google Scholar] [CrossRef]
  17. Ash, J.T.; Cross, W.M.; Svalstad, D.; Kellar, J.J.; Kjerengtroen, L. Finite element evaluation of the microbond test: Meniscus effect, interphase region, and vise angle. Compos. Sci. Technol. 2003, 63, 641–651. [Google Scholar] [CrossRef]
  18. Lin, G.; Geubelle, P.H.; Sottos, N.R. Simulation of fiber debonding with friction in a model composite pushout test. Int. J. Solids Struct. 2001, 38, 8547–8562. [Google Scholar] [CrossRef]
  19. Meng, Q.; Wang, Z. Prediction of interfacial strength and failure mechanisms in particle-reinforced metal-matrix composites based on a micromechanical model. Eng. Fract. Mech. 2015, 142, 170–183. [Google Scholar] [CrossRef]
  20. Dugdale, D.S. Yielding of steel sheets containing slits. J. Mech. Phys. Solids 1960, 8, 100–104. [Google Scholar] [CrossRef]
  21. Barenblatt, G.I. The mathematical theory of equilibrium cracks in brittle fracture. In Advances in Applied Mechanics; Elsevier: Amsterdam, The Netherlands, 1962; Volume 7, pp. 55–129. ISSN 0065-2156. [Google Scholar]
  22. Rizal, S.; Gopakumar, D.A.; Huzni, S.; Thalib, S.; Syakir, M.I.; Owolabi, F.A.T.; Aprilla, N.A.S.; Paridah, M.T.; Khalil, H.P.S.A. Tailoring the Effective Properties of Typha Fiber Reinforced Polymer Composite via Alkali Treatment. BioResources 2019, 14, 5630–5645. [Google Scholar]
  23. Selmi, A. Effect of Glass Fiber Content on the Flexural Modulus of Elasticity of Glass-Epoxy Sandwich Composites. Emerg. Trends Eng. Technol. 2014, 2, 187–191. [Google Scholar]
  24. Dadej, K.; Surowska, B. Analysis of cohesive zone model parameters on response of glass-epoxy composite in mode II interlaminar fracture toughness test. Theory Pract. 2016, 16, 180–188. [Google Scholar]
  25. Potukuchi, S.K.S. Fracture Analysis of Carbon Fiber/Epoxy Matrix Interface through Microbond and Cruciform Tests. Master’s Thesis, The University of Illinois, Urbana, IL, USA, 2016. Available online: https://www.ideals.illinois.edu/bitstream/handle/2142/90456/POTUKUCHI-THESIS-2016.pdf?sequence=1&isAllowed=y (accessed on 30 December 2019).
  26. Żmudzki, J.; Walke, W.; Chladek, W. Influence of model discretization density in FEM numerical analysis on the determined stress level in bone surrounding dental implants. In Information Technologies in Biomedicine; Springer: Berlin/Heidelberg, Germany, 2008; pp. 559–567. [Google Scholar]
  27. Liu, Y.; Glass, G. Effects of Mesh Density on Finite Element Analysis; SAE Technical Paper; SAE International: Warrendale, PA, USA, January 2013. [Google Scholar]
  28. Schellekens, J.C.J.; De Borst, R. A non-linear finite element approach for the analysis of mode-I free edge delamination in composites. Int. J. Solids Struct. 1993, 30, 1239–1253. [Google Scholar] [CrossRef] [Green Version]
  29. Soto, A.; González, E.V.; Maimí, P.; Turon, A.; de Aja, J.R.S.; de la Escalera, F.M. Cohesive zone length of orthotropic materials undergoing delamination. Eng. Fract. Mech. 2016, 159, 174–188. [Google Scholar] [CrossRef]
  30. Zhou, Y.; Fan, M.; Chen, L. Interface and bonding mechanisms of plant fibre composites: An overview. Compos. Part B Eng. 2016, 101, 31–45. [Google Scholar] [CrossRef] [Green Version]
  31. Hao, J.; Zhao, Y.; Luo, Y.; Wang, Y.; Chen, D. Experiments and finite element simulation of interfacial properties for monofilament composites. New Carbon Mater. 2014, 29, 176–185. [Google Scholar] [CrossRef]
Figure 1. The photograph of fiber–matrix interface microbond.
Figure 1. The photograph of fiber–matrix interface microbond.
Computation 08 00002 g001
Figure 2. (a) Microbond test model geometry and (b) boundary condition.
Figure 2. (a) Microbond test model geometry and (b) boundary condition.
Computation 08 00002 g002
Figure 3. (a) Coarse, (b) medium, (c) fine, and (d) very fine mesh.
Figure 3. (a) Coarse, (b) medium, (c) fine, and (d) very fine mesh.
Computation 08 00002 g003
Figure 4. Traction–separation curve. K is interfacial stiffness, G c i is interfacial fracture energy, to is interfacial strength, and δ is the separation.
Figure 4. Traction–separation curve. K is interfacial stiffness, G c i is interfacial fracture energy, to is interfacial strength, and δ is the separation.
Computation 08 00002 g004
Figure 5. Maximum forces and computational times vs. number of elements.
Figure 5. Maximum forces and computational times vs. number of elements.
Computation 08 00002 g005
Figure 6. The load displacement curve of properties cohesive variations on (a) coarse, (b) medium, (c) fine, and (d) very fine mesh compared with experimental values.
Figure 6. The load displacement curve of properties cohesive variations on (a) coarse, (b) medium, (c) fine, and (d) very fine mesh compared with experimental values.
Computation 08 00002 g006
Figure 7. Debonding process.
Figure 7. Debonding process.
Computation 08 00002 g007
Figure 8. Contour plot of shear stress (MPa): (a) vise blade first hit matrix, (b) then hit the matrix, and (c) immediately before debonding.
Figure 8. Contour plot of shear stress (MPa): (a) vise blade first hit matrix, (b) then hit the matrix, and (c) immediately before debonding.
Computation 08 00002 g008aComputation 08 00002 g008b
Figure 9. Maximum traction damage initiation criterion.
Figure 9. Maximum traction damage initiation criterion.
Computation 08 00002 g009
Table 1. The number of elements used in the model.
Table 1. The number of elements used in the model.
Mesh TypeSeed MeshTypha Fiber and MatrixCohesive LayerTotal
CAX4CAX3COHAX4
Coarse0.005225618316754114
Medium0.003589623016757801
Fine0.00158,3241017167561,016
Very Fine0.000518,849131451675193,311
Table 2. Material properties Typha fiber and epoxy resin.
Table 2. Material properties Typha fiber and epoxy resin.
MaterialYoung’s Modulus (GPa)Poisson’s Ratio
Typha fiber0.880.35
Epoxy resin3.70.4
Table 3. Variations in cohesive properties.
Table 3. Variations in cohesive properties.
VariationsStiffness Coefficients (N/mm3)Fracture Energy (N/mm)Damage Initiation (N/mm2)
KnnKssKttGIctnntssttt
VR127002700270014.65270270270
VR227002700270015.15270270270
VR330003000300015.00250250250
VR435003500350014.00270270270
VR550005000500010.50267267267
VR68000800080005.50200200200
VR710,00010,00010,0005.00250250250

Share and Cite

MDPI and ACS Style

Ikramullah; Afrizal, A.; Huzni, S.; Thalib, S.; Abdul Khalil, H.P.S.; Rizal, S. Effect of Mesh Sensitivity and Cohesive Properties on Simulation of Typha Fiber/Epoxy Microbond Test. Computation 2020, 8, 2. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8010002

AMA Style

Ikramullah, Afrizal A, Huzni S, Thalib S, Abdul Khalil HPS, Rizal S. Effect of Mesh Sensitivity and Cohesive Properties on Simulation of Typha Fiber/Epoxy Microbond Test. Computation. 2020; 8(1):2. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8010002

Chicago/Turabian Style

Ikramullah, Andri Afrizal, Syifaul Huzni, Sulaiman Thalib, H. P. S. Abdul Khalil, and Samsul Rizal. 2020. "Effect of Mesh Sensitivity and Cohesive Properties on Simulation of Typha Fiber/Epoxy Microbond Test" Computation 8, no. 1: 2. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8010002

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