Home Analysis on the impact response of fiber-reinforced composite laminates: an emphasis on the FEM simulation
Article Open Access

Analysis on the impact response of fiber-reinforced composite laminates: an emphasis on the FEM simulation

  • Jian He , Liang He EMAIL logo and Bin Yang EMAIL logo
Published/Copyright: August 8, 2018

Abstract

The effects of units, material parameters, and constitutive relationships on the dynamic mechanical response of composite laminates subjected to high- and low-velocity impacts were investigated. Additionally, the role of impact or shape, including hemispherical, flat, and conical, on the damage area of the adhesive layer and displacement of the center of the laminated plates was investigated. The results show that the energy absorption of composite laminates increases with impact velocity, and specific energy absorption changes with the density of the contact surface, which is affected by ply thickness. Moreover, the target energy absorption decreases with increasing layer angle. Under a low-velocity impact, the maximum contact force, damage area of the adhesive layer, and displacement of the center of the laminated plate increase as the impact energy increases, thus showing that impact energy is not directly related to contact duration and energy absorption of composite laminates. The results of different geometric shapes show that the damage area of the adhesive layer and the displacement of the center of the laminated plates are largest for a conical impactor and smallest for a flat impactor.

1 Introduction

As a kind of lightweight material, composite materials are widely used in bulletproof [1, 2] and aerospace fields [3, 4]. Research on composite materials damage under high-velocity impact has theoretical and practical significance for military preparedness [5]. Composite materials will inevitably face unexpected collisions during use. Obvious damage may not be found on the surface under low-velocity impacts, although the composite laminates may suffer internal delamination damage. This internal damage will greatly influence the material’s mechanical properties and lead to the decline of the material bearing capacity [6]. Therefore, it is valuable to study the impact response of fiber-reinforced composite laminates.

Considerable effort has been made in analyzing the damage modes and the impact response of composite materials. Yang et al. [7] built a 3-D finite element model (FEM) to investigate the high-velocity impact response of fiber-metal laminates (FMLs) using ABAQUS/Explicit platform; the results showed that the effect of fiber stacking sequence on the impact performance of FMLs is very limited under critical penetration velocity, and the damage pattern depends on the projectile incident angle. In reference [8], they also verified their mode of interply hybrid composites based on woven fabrics and polymerized cyclic butylene terephthalate resin subjected to low-velocity impact. Long et al. [9] established a damage analysis process of composite laminates subjected to low-velocity impact by building numerical models on a damage model where cohesive contact method was involved according to the experimental observation. Based on the results obtained, they presented a prediction process of delamination shape for composite laminates under low-velocity impact. Zhang et al. [10] proposed a FEM to investigate the dynamic mechanical response and damage modes in cross-ply composite laminates under transverse low-velocity impact. By adopting Hashin criterion and a gradual degradation scheme, the simulation results agree well with the available experimental data. Using the proposed model, they also investigated the effect of interface friction on the delamination response of cross-ply composite laminates under impact. Zhang and Zhang [11] built a simplified and efficient FEM with the delamination damage adjacent to the mid-thickness plane and the lower face taken into consideration to study the impact response and corresponding damage mode of composite laminates subjected to low-velocity impact. The proposed model can deliver acceptable prediction of the critical impact force and damage shape and size. In addition, foam-filled sandwich panels with six types of facesheets were manufactured by Yang et al. [12], and the low-velocity impact test and compression after impact test of the prepared panels were performed. They found that the stacking sequence of fabric clothes has some effects on the impact properties of the hybrid panels.

The paper is organized as follows: section 1 introduces the research background. In section 2, the dynamic response of laminated composites under high-velocity impact is analyzed using the finite element software LS-DYNA, which simulated the laminated plate with a SHELL163 unit without considering the bonding layer, using CHANG-CHANG criterion as the constitutive model. The effects on the energy absorption of the composite laminates, which were produced by the impact velocity, target thickness, and ply angle, were analyzed in detail. In section 3, the low-velocity impact dynamic response of composite laminates was studied. The laminated board layer was simulated with a SOLID164 unit, and the adhesive layer was simulated with a SOLID ELEMENT No. 19 unit, using CHANG-CHANG criterion as the constitutive model. The effects on laminates bonded layer damage and energy absorption, produced by the bonding layer damage development, sub-layer stress change of different ply angles, impact energy, and punch shape, were analyzed in detail. Section 4 concludes the paper.

2 Dynamic response of composite laminates under high-velocity impact

2.1 Damage criterion of laminates

Compared with isotropic materials, the failure modes of composite materials under impact loads are much more complex because of not only the complexity of composite strength itself but also the characteristics of multiphase combinations. From a microscopic angle, laminate damage is mainly divided into three categories. First is the fiber fracture which is dominated by fiber damage, including fiber breakage and micro buckling. Second is the matrix cracking, which is dominated by matrix damage, including meso-cracks and holes. Third is the spatial delamination fracture, including the interface debonding and crack propagation. Fiber breakage and matrix cracking occur mainly inside the composite single plate, while delamination damage occurs mainly within the composite sub-layers. In this paper, the CHANG-CHANG failure criterion is selected based on the LS-DYNA.

(1){fiber fracture failure (σ11>0)(σ11/XT)2+τ¯=1matrix crack failure (σ22>0)(σ22/YT)2+τ¯=1matrix compression failure (σ22<0)(σ222S)2+[(Yc2S)21](σ22YC)+τ¯1

where τ¯=(τ12S)21+32αG12σ1221+32αG12S2,XT, YT, and YC represent the longitudinal tensile strength, the transverse tensile strength, and the transverse compression strength, respectively. α is a constant determined by experiment, and α=0 when considering linear elasticity only.

2.2 The finite element model

In the finite element simulation, a 7.62-mm handgun bullet is used, and the bullet head is simplified as a hemisphere. The SOLID164 unit and the MAT_RIGID model are employed. SOLID164 is an eight-node unit with three degrees of freedom, and the bullet head contains 500 elements. In the process of impact, the bullet is regarded as an ideal rigid material without deformations. The material parameters are shown in Table 1.

Table 1:

Material parameters of bullet.

Density (g/cm3)Young modulus (Gpa)Poisson’s ratioMass (g)
7.852070.35.8

The laminated plate size is 60 mm×60 mm, simulated by the specific SHELL163 unit of ANSYS/LS-DYNA. SHELL163 is a four-node unit with six degrees of freedom. The model is chosen from the material library, and the material parameters are shown in Table 2 . Considering both calculation precision and run time, only one fourth of the complete model is put forward in this paper since the model is symmetric, and the plate contains 800 elements. During the simulation, the processing time varies from 0.5 h to 1 h and depends on the incident velocities.

Table 2:

Material parameters of composite.

Young modulus (GPa)Poisson’s ratioShear modulus (GPa)Bulk modulus (GPa)
EXEYEZNUXYNUXZNUYZGABGBCGCA
23.3990.310.310.391.653.021.6520.4
Shear strength (GPa)Compressive strength (GPa)Tensilestrength (GPa)Density (g/cm3)
SBASCASCBXXCYYCZZCXXTYYTρ
0.390.390.390.650.650.651.3050.131.44

The target plate is subdivided into the bullet shooting region, and the remaining region is relatively rougher. The grid size in the bullet contact center region is 1 mm, and the remainder of the area has a grid size of 2 mm, as shown in Figure 1.

Figure 1: Finite element model of the laminated plates under bullet impact. (A) Bullet model, (B) ply situation, (C) global grid, and (D) local grid.
Figure 1:

Finite element model of the laminated plates under bullet impact. (A) Bullet model, (B) ply situation, (C) global grid, and (D) local grid.

2.3 Analysis of the numerical simulation results

In order to verify the proposed simulation method, we consider the work of Aymerich et al. [13]. According to the experiments, the laminate is simulated with the angle ply sequence [0/90/0/90/0/90/0/90/0/90]. The cylindrical bullet diameter is 7.62 mm, with an aspect ratio of 2. The target size is 20 cm*90 cm, and the thickness is 5 mm. In our simulation, the impact energy is handled by assigning different velocities to the impactor. Simulations are performed with the velocities of 246.6 m/s, 300.2 m/s, 388.8 m/s, 471.2 m/s, and 551.5 m/s, respectively.

The comparison of bullet residual velocity from the numerical simulation and experimental results are listed in Table 3. It can be seen that the user-defined material numerical results match the experimental results, which illustrates that the numerical simulation method in this paper is effective.

Table 3:

Comparison of bullet residual velocity numerical simulation results and experimental results.

Bullet velocity (m/s)246.6300.2388.8471.2551.5
Bullet residual velocity (m/s)62.6190.7302.4416.9498.1
Numerical simulation results (m/s)50.1178.6287.3392.6475.2

2.3.1 Anti-ballistic properties of composite materials against different impact velocities

The anti-ballistic performance of composite materials depends greatly on the velocity effect, which means that in a certain range of velocity, the energy-absorbing capacity may increase with the velocity of the bullet [14]. In the simulation, the one-fourth bullet model penetrates the target at speeds of 200 m/s, 250 m/s, 300 m/s, 350 m/s, 400 m/s, and 450m/s, while the other conditions remain unchanged. Figure 2 shows the kinetic energy variation curve of the bullet at different initial velocities.

Figure 2: Bullet kinetic energy-time curves with different impact velocities.
Figure 2:

Bullet kinetic energy-time curves with different impact velocities.

In order to analyze the relationship between bullet velocity and target plate absorption energy, the relationship between bullet velocity and kinetic energy loss is shown in Figure 3. It can be observed in Figure 3 that in the velocity ranging from 200 m/s to 450 m/s, kinetic energy loss increases with increasing bullet velocity. This finding illustrates that target energy absorption will increase as the initial velocity increases in a certain range, demonstrating the energy absorption speed effect mentioned previously.

Figure 3: Bullet velocity-kinetic energy reduction curve.
Figure 3:

Bullet velocity-kinetic energy reduction curve.

2.3.2 Relationship analysis between target plate thickness and energy absorption effect

It is generally recognized that the area density of the fiber composite laminate is directly related to its energy-absorbing capability. The area density is the quality per unit area in the penetration direction. The area density is related to both the material characteristics and the target board thickness [15]. Since the laminated composite material remains constant in this paper, when target thickness changes, the reduction of kinetic energy under different plate thicknesses can be obtained, which is approximately equal to the energy absorbed by the target board.

Figure 4 shows the bullet kinetic energy loss curves for the target thicknesses of 1 mm, 2 mm, 3 mm, 4 mm, and 5 mm, respectively. It can be observed that, with all other conditions kept constant, bullet residual velocity reduces from 293.01 m/s to 180 m/s with the increase in plate thickness. With the increase in area density, the loss of bullet kinetic energy, which is approximately equal to the target energy absorption, increases from 12 J to 164.32 J. This response occurs mainly due to the fact that the increased target thickness leads to higher energy consumption in the penetration process. Relevant literature indicates that the specific energy absorption is an important standard to measure the material energy absorption ability, which measures the quality of materials’ anti-ballistic properties using the ratio of target energy absorption to area density [16].

Figure 4: Bullet kinetic energy reduction-time curves for target plates with different thicknesses.
Figure 4:

Bullet kinetic energy reduction-time curves for target plates with different thicknesses.

Figure 5 shows the change of specific energy absorption for the target plate thicknesses ranging from 1 mm to 5 mm. The specific energy absorption increases as the target thickness increases. The specific energy absorption increases sharply when the target is thinner and tends to flatten when the target is thicker. This can be explained by the fact that the size of the target is close to the bullet size, the energy produced in penetration cannot be ignored, and it does not belong in the category of thin plate.

Figure 5: Target plate thickness-specific energy absorption curve.
Figure 5:

Target plate thickness-specific energy absorption curve.

2.3.3 Relationship analysis between the ply angle and target plate energy absorption

For composite materials, different ply directions may lead to different mechanical properties. In order to investigate how the degree of anisotropy influences the composite energy absorption, the 3-mm composite plate is subjected to different penetration angles under the same speed of 300 m/s. Figure 6 shows the one-fourth bullet kinetic energy reduction trend under ply angles [0]8, [30]8, [45]8, and [60]8.

Figure 6: Bullet kinetic energy-time curves with different layer angles.
Figure 6:

Bullet kinetic energy-time curves with different layer angles.

At ply angle [0]8, the bullet loses the most kinetic energy, signifying the best anti-ballistic performance. On the contrary, at ply angle [60]8, the bullet loses the least kinetic energy. The laminate energy absorption decreases as the ply angle increases with other conditions unchanged, as shown in Figure 7. This is due to the fact that the increase in ply angle reduces the strain rate strengthening effect; thus, the dynamic mechanical properties are deteriorating [17]. In addition, as the ply angle increases, the reflection of the shock wave in the crossing point increases and the decreased impact response area is not conducive to the spread of energy, thus reducing the anti-ballistic properties of the laminates. To summarize, the laminate ply angle has substantial influence on the absorbed energy.

Figure 7: Ply angle-bullet kinetic energy damage curve.
Figure 7:

Ply angle-bullet kinetic energy damage curve.

By comparing the damage results under different ply angles, it can be found that the breakage shape is approximately an ellipse when ply angle is 0, and the fiber direction is along the short axis because of the orthotropic characteristics of the composite materials. It demonstrates that the mechanical property along the fiber direction is higher than that along its vertical direction. According to the damage pattern obtained from the numerical results, the damage shape of composite laminates changes from an ellipse to a circle as the ply angle changes.

3 Damage response of composite laminates subjected to low-velocity impact

3.1 Damage model of composite laminates subjected to low-velocity impact

In simulation, LS-DYNA is used to simulate and compute the impact. CHANG-CHANG failure criterion is selected as the sub-layer failure criterion, and the E material model for the new version of LS-DYAN program is used in this paper, which is extended from the CHANG-CHANG criterion. The formula expression is as follows:

(2){fiber fracture (σ11>0)(σ11/Xt)2+(σ12/S12)2=1fiber compression (σ11<0)(σ11/Xc)2=1matrix cracking (σ22>0)(σ22/Yt)2+(σ12/S12)2=1matrix compression failure (σ22<0)14(σ22/S12)2+Ycσ22/(4S122)σ22/Yc+(σ12/S12)2=1

A number of experimental studies have illustrated that the delamination damage in composite laminates generally arises between the sub-layers that undergo ply angle changes. To address this concept, a thin layer of interface element layer can be added between the angle-changed sub-layers in the numerical simulation. Then, the delamination damage of composite laminates can be simulated through the interfacial layer damage failure. The interface layer is shown in Figure 8.

Figure 8: Interface element model.
Figure 8:

Interface element model.

The selected SOLID ELEMENT 19 unit is based on COHESIVE ZONE theory, mainly considering the interfacial force Pi (i=1, 2, 3) and the displacement relationship between the top and bottom surfaces [Pi=ki (δi) and δi (i=1, 2, 3)]. The bilinear constitutive model is adopted, as shown in Figure 9.

Figure 9: Bonding layer simplified bilinear model. (A) Shear model. (B) Normal model.
Figure 9:

Bonding layer simplified bilinear model. (A) Shear model. (B) Normal model.

In practical applications, re-coupling or the linear coupling criteria are used as the basis of judgment for a bonding layer crack, as expressed in Equation (3), in which the coefficients GiC(i=a, b, and c) represent the critical energy release rates of three types of cracks.

(3)fg=(Ga/GaC)2+(Gb/GbC)2+(Gc/GcC)2fg=Ga/GaC+Gb/GbC+Gc/GcC

3.2 Numerical simulation of composite materials subjected to low-velocity impact

In Ref. [13], the experiment of composite materials subjected to low-velocity impact was simulated. The ply angle was designated as [0/90/0/90/0/90/0/90/0/90], and the laminate was hinge-constrained with a size of 87.5 mm×65 mm×2 mm and effective center area with size of 67.5 mm×45 mm. The experimental rigid hemispherical punch has the radius of 0.625 cm and a density of 7.8 g/cm3. In order to reduce the finite element running time, the simulation used a one-fourth model of the punch with a SOLID164 solid element.

A SOLID164 element is also used to simulate laminated sub-layers. Three single plates with 0 ply angle were modeled as one sub-layer. A six-layer single plate with a 90 ply angle layer in the middle was modeled as one sub-layer. A three-layer single plate with a 0 ply angle was modeled as one sub layer. The material parameters to model 0 ply angle layers are shown in Table 4, and they can be multiplied by a stiffness transformation matrix to obtain parameters to model for the 90 ply angle layers.

Table 4:

Material parameters of sub layer with 0 ply angle

Young modulus (GPa)Poisson’s ratioShear modulus (GPa)Bulk modulus (GPa)
EXEYEZNUXYNUXZNUYZGABGBCGCA
93.77.457.450.260.260.263.973.973.9720.4
Shear strength (GPa)Compressive strength (GPa)Tensilestrength (GPa)Density (g/cm3)
SBASCASCBXXCYYCZZCXXTYYTρ
0.360.360.360.650.650.651.3060.131.44

Improved CHANG-CHANG failure criterion is used as sub-layer damage failure criterion. An interface layer is inserted between the sub-layers with changing ply angles, forming a laminated structure. Furthermore, an adhesive layer is simulated with specific material unit COHESIVSOLID ELEMENT 19 element in the LS-DYNA971 version. Constraints between top and bottom layers are established to ensure the continuity of the whole model movement. The new model is applied to build a bonding layer model, using the bilinear material model. In order to simplify the calculation, it is assumed that the shear stresses in the two directions on the interface are the same due to the same constitutive model. Moreover, both the tearing type and the sliding type of cracking mode have the same critical release rate, which can be calculated from Equation (3). The following equation can be derived.

(4)k2=k3,GbC=GcC,fg=(Ga/GaC)2+(Gb/GbC)2

In order to precisely obtain the dynamic response of laminates subjected to low-velocity impact, the meshes of the impact region are refined, particularly in the rectangular region of X[0 mm, 15 mm], Y[0 mm,10 mm], as shown in Figure 10.

Figure 10: Geometric modeling diagram.
Figure 10:

Geometric modeling diagram.

3.3 Analysis of the simulation result

To compare the simulation results with the experimental results, the delamination damage chart of the bottom layer and the impact contact force time curve are shown in Figure 11.

Figure 11: The time curve of the impact contact force of the laminate plate for experiment and simulation. (A) The layered damage for experiment in Ref. [13]. (B) The damage of adhesive layer for simulation.
Figure 11:

The time curve of the impact contact force of the laminate plate for experiment and simulation. (A) The layered damage for experiment in Ref. [13]. (B) The damage of adhesive layer for simulation.

3.3.1 Damage and extension on bonding layer

The bonding layer can be damaged easily under low-velocity impact when the direction of the fiber changes. At the same time, the area and shape of the damage as well as the strain should be approximately the same. It is convenient and clear to observe the process of layered production, extension, and evolution of damage by the post-processing method in simulation, as shown in Figure 12.

Figure 12: The damage of the adhesive layer.
Figure 12:

The damage of the adhesive layer.

3.3.2 Plate stress under low-velocity impact

It can be found that the stress nephogram exhibits a peanut-shaped profile for both conditions (layer angle 0° or 90°), and the main axis’s direction of the peanut-shape nephogram is perpendicular to the layer angle of the sub-layer. This observation is consistent with the experimental phenomenon. As time elapses, the strain transfers towards both sides, and the damage increases, as shown in Figure 13.

Figure 13: The Von Mises stress at different times.
Figure 13:

The Von Mises stress at different times.

3.3.3 Damage analysis under different impact energies

In this section, five groups of punches with different initial energies are set up to impact the laminate to investigate the corresponding damages. The results are shown in Table 5.

Table 5:

Conditions of calculation.

Impact energy (J)Punch mass (kg)Impact velocity (m/s)
1.9831.91.444
2.1912.11.444
2.42.31.414
2.72.31.532
32.31.615

As shown in Figure 14, the peak of the impact load increases as the impact energy increases. Within a certain range, the increased punch energy caused by the mass effect affects the contact time between the punch and laminated plate. Moreover, this finding also shows that the duration time of the impact does not change after the energy reaches a certain value.

Figure 14: Time curves of contact force under different impact energies.
Figure 14:

Time curves of contact force under different impact energies.

3.3.4 Effect of impact energy on the damage area and plate center displacement

As shown in Figures 15 and 16, within a certain range, when impact energy increases, the damage area of the adhesive layer and the displacement of the plate center increase in an approximately linear manner.

Figure 15: The variations of delamination area as a function of energy.
Figure 15:

The variations of delamination area as a function of energy.

Figure 16: The curve between the impact energy and damage area of the adhesive layer and displacement at the plate center.
Figure 16:

The curve between the impact energy and damage area of the adhesive layer and displacement at the plate center.

3.3.5 Effect of impact energy on energy absorption of laminated plate

Compared with high-velocity impact, the energy damage of the punch under low-velocity impact is relatively small, as shown in Figure 17. Besides, the damage energy of the punch and the energy absorption of the target plate exhibit small changes for different impact energies. The results show that the initial energy is not the main factor that affects the laminated plate during the low-velocity impact.

Figure 17: The variations of punch kinetic energy as a function of time.
Figure 17:

The variations of punch kinetic energy as a function of time.

3.3.6 Effect of punch shapes on the response of the plate under low velocity

In this paper, three kinds of punches with different shapes are used to analyze the response of the plate under low-velocity impact, as shown in Figure 18. All of the punches have the same mass, radius, and initial impact energy.

Figure 18: The schematic diagram of the impactor geometric shapes (flat, hemispherical, and conical).
Figure 18:

The schematic diagram of the impactor geometric shapes (flat, hemispherical, and conical).

Figure 19 shows the variations of contact force as a function of time for different types of impactors. As seen, the laminated plate impacted by a flat punch is the first to reach the maximum contact force, while the plate impacted by a conical punch is damaged last. In contrast, the laminated plate impacted by a flat punch bears the highest maximum contact force, while the conical punch bears the lowest. Because it takes the longest time to achieve maximum contact force, it gets the most serious damage among these three conditions. The following conclusion can be drawn: Incisive punch causes severe damage to the plate.

Figure 19: The time curve on the impact contact force with different impactor geometric shapes.
Figure 19:

The time curve on the impact contact force with different impactor geometric shapes.

3.3.7 Effect of geometric shapes of the punch on adhesive layer damage and plate center displacement

Table 6 shows that the punch geometric shape does not significantly affect the adhesive layer damage and the plate center displacement. Moreover, the table also shows that the more incisive the punch, the more serious the damage to the laminate plate.

Table 6:

The damaged area of the adhesive layer and the maximum displacement of the plate center caused by the punch with different geometric shapes.

The shape of the punchThe damaged area of the adhesive layer (mm2)The maximum displacement of the plate’s center (mm)
Flat3452.2
Hemispherical3512.5
Conical3623.02

4 Conclusion

In this paper, the dynamic response of composite laminated plates was studied and simulated under both high- and low-velocity impacts. Under high-velocity impact, different types of bullets, acceleration, and energy are analyzed. The results show that the energy absorption of the studied composite laminates increases with higher-velocity impacts, specific energy absorption linearly increases with plate thickness, and the target energy absorption decreases as the layer angle increases.

The contact impact force of laminated plates and the damage evolution of cohesive layer and sub-layers with different angles and under low-velocity impact are discussed. The results show that the maximum contact force, the damage area of adhesive layer, and the displacement of the center of laminated plates increase as the impact energy increases. It also shows that the impact energy is not directly related to impact duration and energy absorption of composite laminates.

The impactors with hemispherical, flat, and conical punches are also discussed. Under the impact of conical punch, the laminated plate experiences the longest contact time, the largest damage area, and the most center displacement. Under the impact of flat punch, the laminated plate experiences the shortest contact time, the smallest damage area, and the least center displacement.

Award Identifier / Grant number: 11702097, 11532013, and 11302054

Funding statement: This work was supported by the National Natural Science Foundation of China (grant numbers 11702097, 11532013, and 11302054, Funder Id: 10.13039/501100001809) and the Natural Science Foundation of Heilongjiang Province of China (grant number A2015012).

References

[1] Shim GI, Kim SH, Ahn DL, Park JK, Jin DH, Chung DT, Choi SY. Compos. Part B-Eng 2016, 97, 150–161.10.1016/j.compositesb.2016.04.078Search in Google Scholar

[2] Zhao L, Liu ML. Mater. Res. 2015, 1090, 22–25.10.4028/www.scientific.net/AMR.1090.22Search in Google Scholar

[3] Alderliesten RC. Mater. Design 2015, 66, 421–428.10.1016/j.matdes.2014.06.068Search in Google Scholar

[4] Wulfsberg J, Herrmann A, Ziegmann G, Lonsdorfer G, Stöß N, Fette M. Procedia Engineer. 2014, 81, 1601–1607.10.1016/j.proeng.2014.10.197Search in Google Scholar

[5] Shao D, Yu M, Lian J, Sawyer S. Nanotechnology 2013, 24, 295701.10.1088/0957-4484/24/29/295701Search in Google Scholar PubMed

[6] Kim EH, Rim MS, Lee I, Hwang TK. Compos. Struct. 2013, 95, 123–134.10.1016/j.compstruct.2012.07.002Search in Google Scholar

[7] Yang B, He L, Gao Y. Sci. Eng. Compos. Mater. 2017, 25, 621–631.10.1515/secm-2016-0226Search in Google Scholar

[8] Yang B, Wang ZQ, Zhou LM, Zhang JF, Liang WY. Compos. Struct. 2015, 132, 464–476.10.1016/j.compstruct.2015.05.069Search in Google Scholar

[9] Long S, Yao X, Zhang X. Compos. Struct. 2015, 132, 290–298.10.1016/j.compstruct.2015.05.037Search in Google Scholar

[10] Zhang C, Duodu EA, Gu J. Compos. Struct. 2017, 173, 219–227.10.1016/j.compstruct.2017.04.017Search in Google Scholar

[11] Zhang J, Zhang X. Compos. Struct. 2015, 130, 85–94.10.1016/j.compstruct.2015.04.023Search in Google Scholar

[12] Yang B, Wang ZQ, Zhou LM, Tong LL, Zhang JF, Liang WY. Compos. Struct. 2015, 132, 1129–1140.10.1016/j.compstruct.2015.07.058Search in Google Scholar

[13] Aymerich F, Dore F, Priolo P. Compos. Sci. Technol. 2009, 69, 1699–1709.10.1016/j.compscitech.2008.10.025Search in Google Scholar

[14] Wang YB, Wang XJ, Hu XZ, Wang F. Engineer. Mechan. 2005, 3, 014.Search in Google Scholar

[15] Aymerich F, Dore F, Priolo P. Compos. Sci. Technol. 2008, 68, 2383–2390.10.1016/j.compscitech.2007.06.015Search in Google Scholar

[16] Tao L, Xue JY, Feng MY, Xiu LZ, Ai L. Eng. Plast. Appl. 2010, 10, 003.Search in Google Scholar

[17] Findik F, Tarim N. Compos. Struct. 2003, 61, 187–192.10.1016/S0263-8223(03)00065-5Search in Google Scholar

Received: 2017-06-30
Accepted: 2018-06-02
Published Online: 2018-08-08
Published in Print: 2019-01-28

©2019 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Analysis on the impact response of fiber-reinforced composite laminates: an emphasis on the FEM simulation
  2. Artificial neural network for predicting the flexural bond strength of FRP bars in concrete
  3. Cyclic behavior of GFRP strengthened infilled RC frames with low and normal strength concrete
  4. Durability of basalt fiber-reinforced polymer bars in wet-dry cycles alkali-salt corrosion
  5. Effect of B4C particle size on the mechanical properties of B4C reinforced aluminum matrix layered composite
  6. Enhanced dielectric properties of BaTiO3 ceramics with cerium doping, manganese doping and Ce-Mn co-doping
  7. Free and forced vibration analysis of rectangular/circular/annular plates made of carbon fiber-carbon nanotube-polymer hybrid composites
  8. Influence of nano-SiO2 on the bonding strength and wear resistance properties of polyurethane coating
  9. Investigation of wear behavior of nanoalumina and marble dust-reinforced dental composites
  10. Negative effect of clay fillers on the polyvinyl alcohol biodegradation: technical note
  11. Photocatalytic activity of Cu2O/ZnO nanocomposite for the decomposition of methyl orange under visible light irradiation
  12. Sub-surface mechanical properties and sub-surface creep behavior of wood-plastic composites reinforced by organoclay
  13. Surface integrity in wire-EDM tangential turning of in situ hybrid metal matrix composite A359/B4C/Al2O3
  14. The influence of the WC-Co composite microstructure model on stress field heterogeneity at the microstructure level: FEM based study
  15. Vibration-damping characterization of the basalt/epoxy composite laminates containing graphene nanopellets
  16. A review on nanocomposite hydrogels and their biomedical applications
  17. Optimization and simulation analysis of structure parameters of OPCM ultrasonic longitudinal wave actuating element
  18. Research Article
  19. Preparation of POSS-triol/wollastonite composite particles by Liquid phase mechanochemical method and its application in UV curable coatings
  20. Research on preload relaxation for composite pre-tightened tooth connections
  21. Dough moulding compound reinforced silicone rubber insulating composites using polymerized styrene butadiene rubber as a compatibilizer
  22. Hydration And Microstructure Of Astm Type I Cement Paste
  23. Effects of NiO content on the microstructure and mechanical properties of AgSnO2NiO composites
  24. Overall buckling behaviour of laminated CFRP tubes with off-axis ply orientation in axial compression
  25. UV sensing optode for composite materials environment monitoring
  26. On crushing characteristics of hybrid sandwich aluminum-cardboard panels reinforced with glass fiber composite rods
  27. Preparation and characterization of Ni-Cu composite nanoparticles for conductive paints
  28. A research on the preparation of oil-adsorbing hydrophobic porous resins by high internal phase emulsions (HIPEs) template
  29. Material characteristics of random glass-mat-reinforced thermoplastic under cryogenic thermal cycles
  30. Differentiation of non-black fillers in rubber composites using linear discriminant analysis of principal components
  31. Research Article
  32. Efficiency of TiO2 catalyst supported by modified waste fly ash during photodegradation of RR45 dye
  33. Synthesis and performance of polyurethane/silicon oxide nano-composite coatings
  34. Study on preparation of magnesium-rich composite coating and performance enhancement by graft modification of epoxy resin
  35. Research Article
  36. Mechanical and wear properties of polyetheretherketone composites filled with basalt fibres
  37. Mechanical Properties of Al 25 wt.% Cu Functionally Graded Material
  38. Research Article
  39. Weight reduction of a carbon fibre composite wheel
  40. Synthesis, electrical properties, and kinetic thermal analysis of polyaniline/ polyvinyl alcohol - magnetite nanocomposites film
  41. Seismic Behaviour of TRC-Strengthened RC Columns under Different Constraint Conditions
  42. Characterization of neat and modified asphalt binders and mixtures in relation to permanent deformation
  43. Microstructures, interface structure and room temperature tensile properties of magnesium materials reinforced by high content submicron SiCp
  44. Research Article
  45. Effect of Cutting Temperature on Bending Properties of Carbon Fibre Reinforced Plastics
  46. Mechanical and tribological properties of B-C-N coatings sliding against different wood balls
  47. Thermal conductivity of unidirectional composites consisting of randomly dispersed glass fibers and temperature-dependent polyethylene matrix
  48. Effects of Waste Eggshells addition on Microstructures, Mechanical and Tribological Properties of Green Metal Matrix Composite
  49. Investigation of porosity effect on flexural analysis of doubly curved FGM conoids
  50. Review Article
  51. Utilization of tailings in cement and concrete: A review
  52. Research Article
  53. Equivalent stiffness prediction and global buckling analysis using refined analytical model of composite laminated box beam
  54. Mechanochemical synthesis of zincite doped with cadmium in various amounts
  55. Size-dependent vibration analysis of graphene-PMMA lamina based on non-classical continuum theory
  56. Automated, Quality Assured and High Volume Oriented Production of Fiber Metal Laminates (FML) for the Next Generation of Passenger Aircraft Fuselage Shells
  57. Research Article
  58. An investigation of the stitching effect on single lap shear joints in laminated composites
  59. The low-velocity impact and compression after impact (CAI) behavior of foam core sandwich panels with shape memory alloy hybrid face-sheets
  60. Effect of granulometric distribution on electromagnetic shielding effectiveness for polymeric composite based on natural graphite
  61. The enhancement of filament winding in marine launching rubber gasbag
  62. Research on ELID Grinding Mechanism and Process Parameter Optimization of Aluminum-Based Diamond Composites for Electronic Packaging
Downloaded on 13.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/secm-2017-0222/html
Scroll to top button