Preprint
Article
Free-edge Effect on the Tensile Properties of 3D Woven Composites
Altmetrics
Downloads
105
Views
29
Comments
0
This version is not peer-reviewed
Submitted:
05 May 2023
Posted:
05 May 2023
You are already at the latest version
Alerts
Abstract
Free-edge effect is one of the factors affecting the mechanical properties of three-dimensional woven composites under tensile load. However, current research is relatively poorly understood regarding the effect of free-edge on the stiffness and strength of the material. This paper aims at examining the influence of free-edge effect on the mechanical properties of 3D woven composites under tension through experimental and simulation mathods. The three-dimensional DIC technique is used to collect the full-field strains on the specimen surface during the test, and the stress-strain differences in different regions in the width direction are analyzed, and the overlap of the curves in each region is found to be high. Experimental studies are conducted on specimens of different widths (within the range of 15-20 mm), and the results indicate that the differences in mechanical properties of 3D woven composites under tension loading in this width range are not significant. A progressive damage finite element model is developed for calculation and compared with experimental results.
Keywords:
Subject: Engineering - Mechanical Engineering
1. Introduction
In the fields of aerospace and shipbuilding, fiber-reinforced composites are extensively used due to their exceptional mechanical characteristics, including high specific stiffness and specific strength. These composites are frequently subjected to tensile forces, and scientists have investigated and examined their tensile properties.
Tensile property is a significant mechanical feature of fiber-reinforced composites, which are excellent engineering materials, and scholars studied this property in depth from various aspects[1,2,3].Hang et al[4] researched the damage and deformation of 2.5D woven composites under tensile loading through a micro/mesoscale analytical framework, analyzed and predicted the strength and elastic properties of fiber bundles, and verified them by comparing with experimental results.
Researchers have conducted numerous early studies on the free-edge effect of composite laminates and found that in many cases the difference in Poisson's ratio between adjacent plies can lead to laminate free-edge damage [5]. The study of the free edge effect in woven composites is further complicated by the presence of inhomogeneous deformation within the single cell [6].
Chao et al [7] focused on the effect of free-edge effect on the feedback of triaxial woven composites under tensile loading and the effect of different specimen widths on the tangential modulus and breaking strength of the material in transverse tension. Finally, the free-edge periodic distribution law was obtained and the relationship between the effective modulus and the specimen width was quantified.
Researchers investigated ways to enhance the mechanical properties of composites by dealing with the free-edge effect. Brauning et al [8] adopted a method of treating the edges of composite sheets with a high-performance adhesive to eliminate the free-edge effect and demonstrated that the method was effective in improving the tensile properties by about 11%.
The influence of the geometry of the specimen on the properties of the composite is also not negligible [9,10]. Li et al [11]found through experimental studies that different preparation angles have significant effects on the tensile properties and damage mechanisms of 3D braided composites. Tl et al [12] used a new modified fixture for combined loading to investigate the difference in the performance of three laminates with the same size and elastic properties and different layup size and layup orientation angles under combined loading. The study by Cai et al [13] on the free-edge effect of triaxially woven composites found that the damage mechanism changed from fiber bundle pull-out to fiber bundle shear as the specimen width increased.
Relatively little research was reported on the free-edge effect on the tensile properties of 3D woven composites, because of the complex and diverse internal weaving forms of three-dimensional woven composites, diverse weaving patterns lead to different material properties [14,15], there is a necessity to study more comprehensively and extensively on the free-edge effect of 3D woven composites.
In this paper, the free-edge effect of three-dimensional woven composites under tensile loading is investigated. The three-dimensional digital image correlation (3D DIC) technique is used to observe the process in order to collect the distribution of full-field strain on the surface of material. The strain variation in different regions along the width of the specimen during the tensile process is analyzed, and the effect of width on the tensile properties of the material is investigated ensuring that the width of the specimen is greater than two single cell dimensions, and the variation of tensile strength, tensile modulus and failure strain are analyzed. The affect of the free-edge effect on the tensile properties of the material is predicted by using a progressive damage finite element model.
2. Materials and Methods
2.1 Materials
Three-dimensional woven composite is one of the typical non-uniform anisotropic material composed of fiber reinforcement and resin matrix, whose internal fiber bundles include warp fiber bundles, weft fiber bundles and Z-directional fiber bundles, which form a fiber reinforcement by complex and diverse winding forms between them. 3D woven composites are classified into three categories: orthogonal woven, angle-interlocked woven and inter-layered interlocked woven composites, depending on their Z-directional fiber bundle winding methods.
The material used in this study is an interlocking three-dimensional woven composite material provided by the Institute of Composites Structure of Tianjin University of Technology. The thickness of the material sheet is 5 mm, in which there are 8 layers of warp fiber bundles and 9 layers of weft fiber bundles, and Z-directional fiber bundles are interspersed between the warp fiber bundles and link the layers of weft fiber bundles as a whole. The structure schematic diagram of the material is shown in Figure 1.
2.2 Test Methods and Equipment
As an advanced composite material, the American Society for Testing and Materials (ASTM) developed some test standards for fiber-reinforced composites under quasi-static loading in the early years, and later the ASTM D 30 Composites Committee revised the original performance test methods several times according to the characteristics of polymer matrix composites, forming such as ASTM D 3039, ASTM D 3410, ASTM D 3518, ASTM D 5379, ASTM D 6641 and other common test standards. The reference standard for the quasi-static uniaxial tension study of 3D woven composites in this paper is ASTM D 3039/D 3039M-08[16].
In order to ensure that there are at least two complete cells in the specimen width, the specimen width is set to 20mm when studying the boundary effect and 15-20 mm when studying the dimensional effect of the specimen, and the specimen size and shape are shown in Figure 2. To reduce the early damage of the specimen caused by the stress concentration in the clamping part at both ends of the specimen, and to increase the clamping friction to prevent the slippage of the specimen during the test, the clamping part at both ends of the specimen is polished smooth and a reinforced aluminum sheet with a thickness of 2mm and a length of 65mm is pasted, and the effective length of the specimen is 120mm.
The uniaxial tensile experiment is conducted at room temperature, and loaded by Lishi – 250 kN Servo-Hydraulic Fatigue Testing Machine. The loading method is displacement-controlled loading with a rate of 0.5 mm/min.
2.3 Three-Dimensional Digital Image Correlation (3D DIC) Technology
A non-contact, full-field strain measurement system developed using digital image correlation (DIC) technology can be used to capture the optical scatter on the surface of the object to be measured, track the displacement data at each point on the surface, and calculate the strain and other data. Compared with the traditional contact strain measurement method, the DIC system has the advantages of no impact on the experimental process due to the actual contact with the specimen, high accuracy, and ease of use, etc. The 3D DIC system establishes a 3D spatial coordinate system in the field of view by obtaining spatial information between cameras through the calibration process. The 3D DIC system can calculate the displacement and strain in one more direction than the 2D DIC system, and the strain obtained by the 3D DIC system can eliminate the fake strain caused by the off-surface movement of the specimen during the experiment when calculating the plane strain.
To obtain the full-field strain on the surface of the test piece during the tensile experiment, the experimental setup shown in Figure 3 is built in this study. The surface of the specimen is photographed with two mutually angled (5°~10° is appropriate) fixed-focus camera lenses, and the two cameras are connected to the image acquisition system, and the shooting frequency is set to 2 Hz, i.e., each camera takes two pictures per second. Clamp the specimen in the upper and lower chuck of the test machine, according to the range of the camera, adjust the distance between the camera and the specimen, the angle between the two cameras to focus, to ensure that the image acquisition system to obtain a clear picture, if necessary, LED lights can be used for lighting compensation. Then, with the help of calibration board, several sets of images with different postures and azimuth angles are taken by two cameras at the same time, which are used as the parameters of the VIC-3D analysis software to calibrate the relative positions of the cameras, so that the real strain field on the surface of the specimen can be obtained.
3. Results
3.1 Boundary Effect Analysis
In order to study the boundary effect of the material, the DIC technique is used to obtain the full-field strain information on the specimen surface during the uniaxial tension of the 3D woven composite. The full-field strains on the specimen surface at different moments are shown in Figure 4, and the numbers on the right side of Figure 4(1–16) are the average strains in the effective area of the specimen surface corresponding to different moments. It can be seen that during the change from the average strain of 1557 to 16481, the strain along the latitudinal direction of the specimen surface is basically the same at each moment, and the overall strain is relatively uniform, which indicates to a certain extent that the boundary effect of the specimen is not obvious.
With the purpose of further studying the influence of the specimen boundary on the overall performance. The strain information of six different regions on the specimen surface are extracted separately, and the average stress-strain curves of different regions are obtained, as shown in Figure 5. The colored square areas in Figure 5a represent the selected five small areas at different locations, and the black area represents the middle large area with the boundary removed. The average stress-strain curves of these six regions are given in Figure 5b. It can be found that the stress-strain curves of these six regions basically overlap before the specimens fracture, indicating that the boundary effect of the test piece does not have much influence on the macroscopic properties of the material. Therefore, in order to save materials and ensure more reliable experimental results, the width of the test piece is selected as 20 mm.
3.2 Size Effect Analysis
In this study, tensile experimental studies are conducted on three different widths (20 mm, 17.5 mm and 15 mm, respectively) to investigate the effect of specimen size on the experimental results, and macroscopic stress-strain curves are obtained for the three widths, as shown in Figure 6.
As can be seen, the dispersion of the stress-strain curves of the three sizes of test pieces is quite minor. The tension experiment results of the three sizes of test pieces are given in Table 1, and it can be observed that the average elastic modulus, average tensile strength and average failure strain of the three sizes are generally consistent, and the dispersion coefficient of the tension results of the three sizes of test pieces is less than 0.03, which indicates that the size effect of the three-dimensional woven composite specimens is not obvious between 15 mm and 20 mm in width.
4. Simulation Verification
4.1. Finite Element Model
In this paper, a progressive damage finite element model is also developed to investigate the effect of free-edge effect on the tensile properties of 3D woven composites [17]. The VUMAT subroutine based on ABAQUS software is written to use different material properties by determining the component material type at the location of the integration point, thus avoiding the difficulty of meshing the fiber bundle and the matrix separately. In case of tensile or compressive damage, the Puck criterion is used as the failure criterion for the fiber bundle and the parabolic criterion for the matrix, while in case of shear damage, the maximum shear stress criterion is used as the failure criterion for both the fiber bundle and the matrix. An exponential model is used to discount the material properties during the damage evolution. The macroscopic full-size model is established according to the test dimensions, and the schematic diagram of the model is shown in Figure 7. The mesh is divided into positive hexahedral cells with a mesh size of 0.5 mm, and the cell type is an eight-node reduced integral cell (C3D8R), and the boundary conditions are imposed according to the real force conditions during the test to simulate the damage failure behavior of 3D woven composite materials under tensile loading.
4.2 Model Validation
The damage and failure behaviors of the 3D woven composite under uniaxial tensile loading are obtained by the above finite element model calculation, and the final damage schematic is shown in Figure 8. The final damage and failure of the finite element model shows that the final failure is severe at a certain cross-section of the model mid-section and a large number of units fail, which is consistent with the actual material failure.
Different regions are divided according to the width direction of the specimen, and the stress-strain curves are extracted separately, as shown in Figure 9. The stress-strain curves in the model in the figure for different regions along the width direction show that the curves in each region largely overlap before the material fails and damages. This indicates that the mechanical properties of the edge and middle regions of the material under tensile loading are basically the same.
It can be seen that the final damage situation obtained from the simulation is consistent with the damage situation obtained from the test, and the test and simulation curves basically match, indicating that the model can effectively predict the failure and surface strain distribution of the 3D woven composite under uniaxial tensile loading.
4.3 Calculation Results Analysis
The tensile properties of the specimens with different widths are obtained and the data extracted by changing the model width through the above model calculations, and the results are shown in Figure 10.
From Figure 10, in the width range of 15-20 mm, the failure strain, elastic modulus and tensile strength of the material vary somewhat but are within a reasonable range. When the specimen width is less than 15 mm, the mechanical properties of the material under tensile load are somewhat weakened with the decrease of the specimen width.
The reason for this phenomenon is that the three-dimensional woven composites are assembled by the arrangement of single cells, and when the width of the specimen is smaller than the size of two single cells, the mechanical properties of the material under tensile loading are reduced.
5. Conclusion
In order to study the free-edge effect of three-dimensional woven composites under uniaxial tension, the 3D DIC system is used in this paper to collect and analyze the full-field strain of the specimens during the test, and the tensile experiments are studied for specimens of different widths on the basis of ensuring that the width of the specimens is larger than two single cell sizes, and the following conclusions are obtained:
- Under the uniaxial tension, the surface strain changes of the three-dimensional woven composites are generally consistent and uniformly distributed along the width direction, and the stress-strain curves of the specimens in different areas along the width direction are in good agreement, indicating that the boundary effect has a small effect on the tensile properties of the three-dimensional woven composites.
- The stress-strain curves of the three-dimensional woven composite specimens with different widths (15-20 mm) under uniaxial tensile loading overlap well, and the average elastic modulus, average tensile strength and average failure strain dispersion coefficients of the three types of specimens are less than 0.03, demonstrating that the size effect of the specimens in this width range is not obvious.
- By using the progressive damage finite element model, it is found that the stress-strain curves in different regions of the same specimen width direction under tensile loading of three-dimensional woven composites basically the same, and are not greatly affected by boundary effects. and the mechanical properties of different width specimens have not much change. When the width of the specimen is less than twice the size of the single cell, the tensile mechanical properties of the material are decreased.
Author Contributions
Conceptualization, Z.W.; methodology, G.L.; software, C.L.; validation, Z.W. and F.L.; formal analysis, S.Z.; investigation, C.L.; resources, G.L.; data curation, F.L.; writing—original draft preparation, C.L.; writing—review and editing, G.L.; visualization, F.L.; supervision, Z.W.; project administration, S.Z.; funding acquisition, Z.W.; G.L. and S.Z. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the National Natural Science Foundation of China (No. 12102215), the Natural Science Foundation of Zhejiang Province (No. LQ21A020001), the Jiangxi Provincial Natural Science Foundation (Grant No. 20212BAB211007 and Grant No. 20202BAB211005), Department of Education of Jiangxi Province (Grant No. GJJ190328).
Institutional Review Board Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
Not applicable.
Conflicts of Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
References
- Deng, L.; Hao, Z.; Zhang, L.; Liu, L.; Experimental and numerical investigation of progressive damage and failure behavior for 2.5D woven alumina fiber/silica matrix composites under a complex in-plane stress state. Composite structures. 2021, 270: 114032. [CrossRef]
- Kohlman, L.W.; Bail, J.L.; Roberts, G.D.; Salem, J.A.; Martin, R.E.; Binienda, W.K. A notched coupon approach for tensile testing of braided composites. Compos. Part A: Appl. Sci. Manuf. 2012, 43, 1680–1688. [Google Scholar] [CrossRef]
- Liu, Y.; Zhang, L.; Li, Z.; Chen, Z.; Huang, K.; Guo, L. Investigation on damage evolution of open-hole plain woven composites under tensile load by acoustic emission signal analysis. Compos. Struct. 2023, 305, 116481. [Google Scholar] [CrossRef]
- Hang, C.; Cui, H.; Liu, H.; Suo, T. Micro/meso-scale damage analysis of a 2.5D woven composite including fiber undulation and in-situ effect. Compos. Struct. 2021, 256, 113067. [Google Scholar] [CrossRef]
- Mittelstedt, C.; Becker, W. Interlaminar Stress Concentrations in Layered Structures: Part I - A Selective Literature Survey on the Free-Edge Effect since 1967. J. Compos. Mater. 2004, 38, 1037–1062. [Google Scholar] [CrossRef]
- Kohlman, L. W. Evaluation of test methods for triaxial braid composites and the development of a large multiaxial test frame for validation using braided tube specimens. The University of Akron; 2012.
- Zhang, C.; Binienda, W.K.; Goldberg, R.K. Free-edge effect on the effective stiffness of single-layer triaxially braided composite. Compos. Sci. Technol. 2015, 107, 145–153. [Google Scholar] [CrossRef]
- A Brauning, K.; Kunza, A.; Alarifi, I.M.; Asmatulu, R. Mitigations of machine-damaged free-edge effects on fiber-reinforced composites. J. Compos. Mater. 2020, 55, 1621–1633. [Google Scholar] [CrossRef]
- Kueh. ; Ahmad, B. H. Size-influenced mechanical isotropy of singly-plied triaxially woven fabric composites. Compos. Part A Appl. Sci. Manuf. 2014, 57, 76–87. [Google Scholar] [CrossRef]
- Flauder, S.; Bombarda, I.; D'Ambrosio, R.; Langhof, N.; Lazzeri, A.; Krenkel, W.; Schaffoner, S. Size effect of carbon fiber-reinforced silicon carbide composites (C/C-SiC): Part 2—tensile testing with alignment device. J. Eur. Ceram. Soc. 2021, 42, 1227–1237. [Google Scholar] [CrossRef]
- Li, D.-S.; Han, W.-F.; Jiang, L. On the tensile properties and failure mechanisms of 3D six-directional braided composites at elevated temperatures. Compos. Commun. 2021, 28, 100884. [Google Scholar] [CrossRef]
- Laux, T.; Gan, KW.; Dulieu-Barton, JM.; Thomsen, OT. Ply thickness and fiber orientation effects in multidirectional composite laminates subjected to combined tension/compression and shear. Compos. Part A: Appl. Sci. Manuf. 2020, 133, 105864. [Google Scholar] [CrossRef]
- Cai, Y.; Zhao, Z.; Tie, Y.; Cao, Y.; Chen, J.; Binienda, W.K.; Zhang, C. Size-dependency of the transverse-tensile failure behavior for triaxially braided composites. Compos. Sci. Technol. 2021, 206. [Google Scholar] [CrossRef]
- Zhou, G.; Sun, Q.; Meng, Z.; Li, D.; Peng, Y.; Zeng, D.; Su, X. Experimental investigation on the effects of fabric architectures on mechanical and damage behaviors of carbon/epoxy woven composites. Compos. Struct. 2021, 257, 113366. [Google Scholar] [CrossRef] [PubMed]
- Jiao, W.; Xie, J.; Yang, Z.; Fang, J.; Chen, L. Effect of weaving structures on the geometry variations and mechanical properties of 3D LTL woven composites. Compos. Struct. 2020, 252, 112756. [Google Scholar] [CrossRef]
- ASTM_International A. D3039: Standard test method for tensile properties of polymer matrix composite materials. West Conshohocken (PA): ASTM International; 2008.
- Liu, G.; Zhang, L.; Guo, L.; Liao, F.; Zheng, T.; Zhong, S. Multi-scale progressive failure simulation of 3D woven composites under uniaxial tension. Compos. Struct. 2019, 208, 233–243. [Google Scholar] [CrossRef]
Figure 1.
The structure diagram of 3D woven composites.
Figure 2.
Uniaxial tensile specimen of 3D woven composites.
Figure 3.
Setup of 3D DIC Technology.
Figure 4.
Surface strain nephogram of specimen during uniaxial tension.
Figure 5.
Boundary effect analysis of uniaxial tension (a) schematic diagram for selecting different areas on the surface of a specimen (b) stress-strain curves of different area.
Figure 5.
Boundary effect analysis of uniaxial tension (a) schematic diagram for selecting different areas on the surface of a specimen (b) stress-strain curves of different area.
Figure 6.
Stress-strain curves of specimens with different widths.
Figure 7.
Schematic diagram of finite element model (a) load condition (b) mesh
Figure 8.
Final damage diagram of finite element model
Figure 9.
Comparison of stress-strain curves in different areas
Figure 10.
Finite element simulation results of tensile properties for different width models
Table 1.
Uniaxial tensile experimental data of specimens with different widths
| Specimen width | Specimen number | Geometrical dimension（mm） | Elastic modulus（GPa） | Tensile strength（MPa） | Failure strain |
|---|---|---|---|---|---|
| 20mm | UT-20-1 | 250×20×5 | 53.583 | 897.41 | 1.69% |
| UT-20-2 | 53.103 | 909.976 | 1.68% | ||
| UT-20-3 | 51.416 | 894.933 | 1.74% | ||
| Average value | 52.701 | 900.773 | 1.70% | ||
| 17.5mm | UT-17.5-1 | 250×17.5×5 | 55.692 | 932.001 | 1.71% |
| UT-17.5-2 | 52.337 | 873.523 | 1.68% | ||
| UT-17.5-3 | 54.585 | 894.821 | 1.69% | ||
| Average value | 54.205 | 900.115 | 1.69% | ||
| 15mm | UT-15-1 | 250×15×5 | 54.323 | 874.693 | 1.67% |
| UT-15-2 | 52.625 | 921.112 | 1.69% | ||
| UT-15-3 | 55.474 | 936.882 | 1.72% | ||
| Average value | 54.141 | 910.896 | 1.69% | ||
| Dispersion coefficient | 0.027 | 0.025 | 0.013 | ||
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
MDPI Initiatives
Important Links
© 2024 MDPI (Basel, Switzerland) unless otherwise stated