Next Article in Journal
Hydrogen Recovery from Waste Gas Streams to Feed (High-Temperature PEM) Fuel Cells: Environmental Performance under a Life-Cycle Thinking Approach
Previous Article in Journal
Optimal Robot Motion for Accelerated Life Testing of a 6-DoF Industrial Robot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability and Failure of the Edge-Closed Honeycomb Sandwich Panels with Face/Core Debonding

Department of Aircraft Design, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(21), 7457; https://doi.org/10.3390/app10217457
Submission received: 19 September 2020 / Revised: 20 October 2020 / Accepted: 21 October 2020 / Published: 23 October 2020

Abstract

:
This study investigated the mechanical performance of the edge-closed honeycomb sandwich structure with face/core debonding under compressive load by experimental and numerical methods. Uniaxial compression tests of asymmetric sandwich structures with various debonding sizes between the carbon/epoxy face sheets and the honeycomb core were conducted. The experimental results showed that the failure of debonding specimens was mainly caused by the local buckling of face sheets at the debonding area. The failure zone of sandwich structures gradually translated from the edge-closed beveled area to the debonding area with the increase of debonding sizes correspondingly, Meanwhile, the stability and the load carrying capacity of sandwich structures were in a downtrend. Compared with the nondestructive specimens, the residual strength of two kinds of defective specimens decreased by 9.9% and 22.1%, respectively. The simulation of the developed numerical model based on the linear buckling theory and the continuum damage mechanics agreed well with the experimental data. The study has guide to assess on the damage tolerance of edge-closed honeycomb sandwich panels with face/core debonding.

1. Introduction

Honeycomb sandwich composite is one of the most effective material structures applied to improve the bending stiffness and reduce the weight of the structure. In recent years, it has been widely used in aerospace aircraft structures such as the fuselage, the wings, the rudders and so on [1]. For instance, the proportions of honeycomb sandwich used in B-58 high-altitude bombers and F-111 fighters as high as 85% and 90%, respectively [2]. However, defects of honeycomb sandwich structures, such as cracks, debonding, core crush, and laminate delamination, etc., which generally caused by alternating loads, external impacts, lightning strike and environment during manufacturing and using, greatly affects the mechanical properties. Among them, face/core debonding is the most universal and dangerous defect in honeycomb sandwich structures which undermines the structural integrity seriously [3].
Studies on characterizing face/core debonding of sandwich structures have been conducted. In order to measure the mode I critical fracture toughness, the most common test method adopted is single cantilever beam (SCB) sandwich specimen which has been selected as the candidate of ASTM International test standard [4]. The SCB test method is easy to perform; Saseendran et al. [5] used this test method to measure mode I critical fracture toughness of various sandwich configurations, and concluded that the fracture toughness increased with the rise of core density, face-sheet thickness, and core cell size. Saseendran et al. [6] also adopted this method to examine the effect of face and core moduli and thickness on crack rotation in SCB sandwich specimen. For the measurement of mode II critical fracture toughness, the Cracked Sandwich Beam (CSB) test [7] and End-Notched Flexure (ENF) test [8] are the common test methods. As for the test method under mix-mode condition, the Mixed Mode Bending (MMB) test is used [9,10], but the possible range of mode mixity phase angle in a MMB test is limited [11]. Therefore, Saseendran and Berggreen [11] adopted the Double Cantilever Beam–Uneven Bending Moments (DCB-UBM) test method to overcome the limitation of MMB test method.
Compressive strength is seriously affected by face/core debonding of sandwich structures [12]. The comprehensive experimental study on compression failure of sandwich structures with face/core debonding was firstly carried out by Vadakke and Carlsson [13]. The results showed that failure of sandwich structures was mainly caused by buckling of debonding face sheet, and followed by rapid propagation of the debonding. Test results given by Aviles and Carlsson [14] showed that panels with square debonding failed at lower loads than those with circular debonding under in-plane compression, and the interface (core/resin) failure increased with the core density increased. Wenzhi et al. [15] studied the edgewise compressive strength of Carbon Fiber Reinforced Plastic (CFRP) sandwich structure with debonding area of different shapes and sizes, which showed that the size of the rectangular debonding defect had a slight effect on edgewise compressive strength, and the larger area of the elliptical debonding caused edgewise compressive strength to decrease by 18%. Nieh et al. [16] studied the effects of panel thickness and debonding length on the compressive strength and failure mechanism of composite sandwich panels through numerical simulation and experimental methods. The results showed that specimens of thicker sheets and smaller debonding area had higher buckling loads and failure strength. Pietrek and Horst [17] demonstrated that the failure load and debonding behavior of asymmetric sandwich structures with aluminum sheets were related to the size and depth of delamination by in-plane compressive tests. In terms of simulation, the buckling behavior of sandwich structures was developed by Frostig and Sokolinsky [18] based on high-order sandwich theory. The propagation behavior of interface debonding and core crack was predicted by El-Sayed and Sridharan [19] by a model with cohesive elements. Berggreen and Simonsen [20] developed a 3D model based on the finite element method to predict the residual strength of debonded sandwich panels. Based on finite element analysis and linear elastic fracture mechanics, Moslemian et al. [21] studied the initial buckling load and compressive strength of a delaminated sandwich panel with a circular core. Sayyidmousavi et al. [22] used a new 3D finite element modeling method to study the effects of face/core debonding size, shape, aspect ratio, fiber direction, and core stiffness on the buckling behavior of honeycomb panels under different boundary conditions.
The objects of the above studies are mostly for sandwich structures of the element level and constant thickness. However, the edges of the honeycomb sandwich structure are usually beveled and closed [23], which is easy to connect and transfer load. Castanie et al. [24] developed a new technology, which was well known for its asymmetric sandwich structure, and it was used to design lightweight structures. Yu [2] experimentally studied the effect of debonding defects in the beveled area on the failure process when edges-closed honeycomb sandwich structure was subject to compressive loads. The results showed that the debonding defect in the beveled zone had a small effect on the compressive strength of the structure, but had a large effect on the strain of the inner panel. Pan et al. [25] studied the influence of plate/core debonding defects on the buckling load and instability mode of edge-closed honeycomb panels by finite element method. Deng et al. [26] studied the stability and failure behavior of asymmetric sandwich structures under asymmetric compressive load through experiments and finite element analysis.
However, there was less research on the edge-closed honeycomb sandwich structure with debonding defects in the public literature. In this paper, experimental and simulation studies on the mechanical properties under axial compression load of the edge-closed honeycomb sandwich structure were developed. Two sizes of face/core debonding area were designed according to variation law of buckling load calculated by finite element analysis, and corresponding compression tests were designed. Then, the buckling load, ultimate load and failure mode of edge-closed honeycomb structures with different debonding sizes under axial compressive load were obtained. The results of experiments were in good agreement with the results of the progressive failure model of face/core debonding composite sandwich structures established by combining cohesive zone model (CZM) technique and progressive failure analysis. The present study had guiding significance for the assessment on damage tolerance of edge-closed honeycomb sandwich panels with face/core debonding.

2. Experimental Procedures

2.1. Subsection

The specimens designed in this paper were divided into non-defective samples and samples with face/core debonding, which were consist of skin, honeycomb core and adhesive. The sheets were stacked by plain carbon fabric prepreg Cycom970/T300, provided by Cytec Industries Inc. (Paterson in United States). The thickness of prepreg was 0.216 mm and mechanical parameters are given in Table 1. Nomex honeycomb JY1-3.2–48 manufactured by Jiangsu Junyuan New Material Co., Ltd. (Nantong, China). was used in this study, and its equivalent mechanical properties [27] were given in Table 2. TS-JM-200 was used to adhere face-sheet onto core, and its mechanical properties are given in Table 3. All specimens were produced in autoclave.
The basic configuration of the two type of test specimens is shown in Figure 1. The flat panel of specimens was the outer skin, and the convex panel was the inner skin. The square area between the inner and outer skins was the honeycomb core. The design bases [28] of specimens are summarized in Table 4. The stacking sequence of sheets in honeycomb area from outside to inside was [(45/0)2/45/C10/45/0/45], and that of edge-closed area from outside to inside was [(45/0)2/45/(0/45)S/45/0/45]. More layers and reinforced plates were laid in the closed area to prevent the edge-closed area from destroy in advance.
Compared with non-defective specimens, specimens with debonding defects form debonding areas by prefabricating a polytetrafluoroethylene film with a diameter D between the outer skin and the honeycomb core. For selecting the representative debonding size D, ABAQUS was used to analyze the influence of eight different debonding sizes between 10 mm and 80 mm on the buckling load of honeycomb structure, as shown in Figure 2. According to the change law of the first-order buckling load of the structure with the debonding size, the debonding diameter with a large load change was finally selected, 40 mm and 80 mm respectively.
The test matrix in this paper is shown in Table 5. For research purposes, all specimens were numbered as CDX-Y. C represented the compressive load, DX represented the debonding size, and Y represented the serial number of various samples. For example, the number CD0-1 represents the first non-defective compressive sample.

2.2. Test Procedure

Back-to-back resistance strain gages (gage numbering was BE 120-3BC which was provided by Zhonghang Electronic Measuring Instruments Co., Ltd.) were arranged and grouped on both sides of the compressive specimens. Only longitudinal strain was required in result processing, so the longitudinal strain gages are drawn in Figure 3. For specimens with defects, all gages were divided into four group. (1) Group one: the resistance strain gages stuck at the center of the beveled area near the edge-closed area longitudinally (numbered 8 and 14 on inner skin, 1 and 7 on outer skin); (2) group two: the resistance strain gages set on the intersection points of the horizontal and vertical 1/4 lines of the honeycomb area (numbered 2, 3, 5 and 6 on inner skin, 9, 10, 12 and 13 on outer skin); (3) group three: the resistance strain gages at the center of the specimens (numbered 4 and 11 on inner and outer skin respectively); (4) group four: the resistance strain gages arranged on the intersection of the debonding boundary and the centerline of the honeycomb zone (numbered 15 and 17 on inner skin, 16 and 18 on outer skin). The fourth group of strain gages were not arranged in non-defective compressive specimens.
The INSTRON static hydraulic testing system was used as the loading equipment for this compressive test. The maximum applying force of this machine was up to 1500 KN. The JM3813 data acquisition system was used to measure strain of specimens. The size of the loading ends of specimens exceeded the width of loading heads of testing machine, so it was necessary to design a suitable fixture to support the ends of the specimens. The structure of the fixture is shown in Figure 4. The U-shaped loading ends of fixture were connected to the loading end and the fixed end of the specimen by bolts. In order to avoid the direct contact between the bolts and specimen, which could cause the surface of sample to be damaged, cushion blocks were added between the U-shaped loading ends and the sample. The columns on both sides were used to prevent specimens from being destabilized during loading and causing specimens to fail in advance. The displacement control method was used during test, and the loading rate was set to 0.5 mm/min. The strain was measured every 5 kN until the specimen was failed. Finally, load-displacement curve was output, and the maximum load value was recorded.

3. Test Results and Discussion

3.1. Strain Distribution

As the quantity and variety of specimens were too large, only strain measurement results of typical specimens of which experimental phenomena were easy and clear to be observed were selected for analysis in this section, and the trends of the strain-load curves of other samples were basically consistent. As shown in Figure 5, Figure 6 and Figure 7, the strain values of the inner and outer skin remained the same only at the initial stage of loading, this was attributed to the beveled area closed to the edge of the specimen which led to the initial loading eccentricity. Thereafter, the strain values of the back-to-back resistance strain gages were different, which indicates that the structure was subjected to bending load due to the existence of eccentric loads. When the slope of the strain-load curves showed a reverse trend (vertical red dotted line in the figure), it indicated that the structure was buckled.
The strain-load curves of the non-defective specimens (CD0-3) are shown in Figure 5. When the load reached 65 kN, the slope of strain-load curves of outside skins in the beveled zone turned to a positive value (such as 8 in Figure 5). Subsequently, the slope of the curves at other positions also showed a reverse trend (such as 3, 4 in Figure 5). The beveled area was the weak area of the edge-closed honeycomb structure, so the beveled area was generally damaged before other locations, which indicated that the structure had been damaged before the overall buckling occurred.
The strain-load curve of the test piece (CD40-3) with debonding defects is shown in Figure 6. When the load reached 60 kN, the slope of the strain-load curves of inner sheets in the debonding area turned to a positive value (such as 4, 17). The slope of the curve of the inner sheets around the honeycomb area changed from negative value to positive one (such as 3, 5) when the load continued to increase to 65 kN. Finally, the curve slope of the outer sheets in the beveled area showed a reverse trend (such as 14) when the load was 75 kN. Firstly, the sandwich panels with a defective size of 40 mm locally buckled in the debonding area, and then local debonding of the panel and the honeycomb core occurred in the beveled area near the loading end.
The strain-load curves of the test specimens (CD80-3) with debonding defects are shown in Figure 7. When the load reached 35 kN, the curve slope of the inner sheets in the center of the debonding area showed a reverse trend (such as 4). When the load continued to 50 kN, the slope of the curve of the outside sheets in the beveled area changed to positive value (such as 8). Meanwhile, the strain curves in the debonding area changed randomly (such as 15, 17). Therefore, the sandwich structure with a defect size of 80 mm firstly suffered from local buckling in the defect area at a lower load, and subsequently damage occurred in the debonding area and beveled area.
The strain-load curves of gages 4 and gages 11 in the center of the inner and outer sheets of the different type specimens are selected to drawn in Figure 8, and are contrasted to discuss. For edge-closed honeycomb sandwich panels, the phenomenon of strain-load curves’ slope changing from negative to positive occurred on the inner sheets (gages 4), therefore, the inner sheets with beveled section was more prone to buckling or damage. In addition, the larger the debonding size was, the more likely the curve slope was to have inflection points of positive and negative value under low load.

3.2. Fracture Load

Owing to the initial gaps in the fixture, specimens and loading heads of testing machine, there was a rigid body displacement at the beginning of the test, which resulted in “unreal” initial non-linearity in the load-displacement curves of all the test samples in Figure 9. Table 6 shows the ultimate loads and their average values for all test samples. It can be seen from the Table 6 that the data had little discreteness (the coefficient of variation of the results are 11.25%, 1.39%, 5.46% respectively), indicating that the manufacturing process of different specimens had good consistency. Comparing the results of three types tests, the ultimate loads of all specimen, especially that of CD80, decreased significantly with the increase of the debonding sizes.

3.3. Ultimate Failure Modes

The failure mode of the specimen CD0-3 is shown in Figure 10. Due to the additional bending moment caused by the eccentric load, the structure of the CD0-3 bulged out along the inner sheet. The honeycomb core in the beveled area was crushed and sheared, which caused the sheets to lose support and form a large area of depression. Meanwhile, fiber crushing and delamination occurred in the reinforced area. Finally, the structure failed and could not continue to carry.
The failure mode of the specimen CD40-3 is shown in Figure 11. The honeycomb core between debonding area and beveled area was crushed and sheared, so the inner and outer skins lost support and a large area bulged out. Meanwhile, fiber crushing and delamination occurred on the both sides of laminated area. As a result, the structure could not continue to carry.
The failure mode of the CD80-3 specimen is shown in Figure 12. The honeycomb core in the defect area and the beveled area was crushed and sheared due to the local buckling of the specimens. The inner skin of the defect area was dented and the outer skin bulged. Damage continued to accumulate and the structure could not continue to carry. Unlike the failure modes of CD0 and CD40 specimens, no fiber breakage and delamination were found in the inner and outer skins on both sides of the lamination zone during the entire test.

4. Progressive Damage Model

4.1. Failure Theory of Honeycomb Sandwich Panels

Progressive Failure Analysis (PFA) is a computational method based on the assumption that damaged materials can continue to bear the load according to their degraded properties [29]. It is widely used in strength prediction and damage propagation analysis of composite structures. For sandwich structures under compression, damage may initiate in face-sheets, in the core and at the interface between sheets and core. Therefore, the progressive failure model of face/core debonding composite sandwich structures was established to predict the strength and observe the failure process. The process of progressive failure analysis is shown in Figure 13.

4.1.1. Intra-Laminar Progressive Failure

The material properties of outside the plane could be ignored because of the relatively thin thickness of composite face-sheet, and the structure was considered to be in a state of plane stress. Therefore, the two-dimensional hashin criterion [30,31] was used to simulate the panel failure in this paper. Four distinct failure modes were modeled, including fiber compressive failure (FCF), matrix compressive failure (MCF), fiber tensile failure (FTF) and matrix tensile failure (MTF). Considering the brevity of this paper, the abbreviations of four failure modes were used to introduce failure process in following sections.

4.1.2. Debonding

The cohesive zone model (CZM) technique based on damage mechanics and fracture mechanics has an irreplaceable advantage in predicting the initiation and expansion of cracks. This technique was originally proposed by Dugdale [32] and Barenblatt [33] to deal with the fracture process area at crack tips. Heretofore, CZM technique has been widely used to simulate the inter-lamina of composite laminates [34] and the adhesive of sandwich composites [5]. The conclusions in literature [5] show that the CZM approach is robust, and is able to predict debonding phenomenon of sandwich structure accurately.
Several cohesive laws exist charactering the relationship between displacement and traction of the crack, among them, the bilinear law is widely used, and it can take into account both the precision and efficiency of calculation [35]. Therefore, bilinear constitutive model is adopted in this paper to model the adhesive.
Damage onset occurs when the following quadratic stress failure criterion is satisfied [36]
( σ n N ) 2 + ( τ s S ) 2 + ( τ t T ) 2 = 1
where 〈 〉 is the Macaulay bracket, σ n , τ S and τ t represent the normal and two shear tractions, N, S and T are the normal and two shear strengths.
A fracture energy-based linear softing law is adopted to address the damage evolution of the adhesive, which gives
σ d δ = G C
where σ and δ are the effective traction and displacements at interfaces. Gc is the mixed-mode fracture energy release rate which determined by Benzeggagh and Kenane criterion [37]. As shown by Equation (3).
G C = G I C + ( G I I C G I C ) ( G I I G I + G I I ) η
where GIC and GIIC are the critical fracture energy release rate of mode I and mode II. η is an interaction parameter that describe the effect of mixed-mode ratio GII/(GI + GII) of total fracture toughness Gc.

4.1.3. Honeycomb Core Progressive Failure

The honeycomb core is equivalent to an orthotropic structure by using the sandwich board theory, and the Besant criterion [38] is adopted as the damage criterion of the honeycomb core. Its expression is as follows:
F c o r e 2 = ( σ T T X T T ) n + ( τ L T S L T ) n + ( τ W T S W T ) n
where, σ T T , σ L T , σ W T are the normal stress of the honeycomb core in the T direction and the shear stress in the LT and WT direction respectively. XTT, XLT, XWT are the compressive strength and shear strength of the honeycomb core in the corresponding direction respectively. The power index n = 2. When the damage variable Fcore = 1, it is considered that the honeycomb core is damaged.
The direct degradation method can efficiently and accurately simulate the evolution of the honeycomb core at various stages. The stiffness reduction method of the honeycomb core in the User-defined Material Mechanical Behavior (UMAT) subroutine is as follows:
E 11 = λ 11 E 11 , E 22 = λ 22 E 22 , E 33 = λ 33 E 33 G 12 = λ 12 G 12 , G 13 = λ 13 G 13 , G 23 = λ 23 G 23
The superscript in the formula represents the stiffness after reduction. λ i j ( i , j = 1 , 2 , 3 ) is the reduction factor in each direction, and the factors [39] in different damaged status are given in Table 7.

4.2. Finite Element Model

The progressive failure model of face/core debonding composite sandwich structures is modelled in the well-known finite element (FE) software ABAQUS supported by Dassault SIMULIA (Dassault Systèmes, Vélizy-Villacoublay, France). An UMAT subroutine contained the damage initiate and propagation criterion of honeycomb core is coded to simulate the progressive failure of Honeycomb core in ABAQUS/Standard platform. Figure 14 illustrates the grid discretization and boundary conditions of the edge closed honeycomb panel with debonding defects. The 0° layering direction of composite material follows the red arrows direction in the Figure 14 (A) of the cross-sectional view A-A. The continuous shell element SC8R (green elements) was used to simulate each layer of the inner and outer skins, and the total number of elements is 7810. The three-dimensional solid element C3D8R (white elements) was used to simulate the honeycomb core, and the total number of elements is 16426. In order to accurately express the face/core debond, the cohesive element COH3D8 (red elements) was used to simulate the adhesive connecting the skin with the honeycomb core. The thickness of adhesive was set to 0.02 mm, and the total number of such elements was 7382. The failure methods in Section 4.1 were used to judge the failure of the three materials. To ensure the continuity of the displacement of the element nodes, the skin, the adhesive and the honeycomb core were connected by common nodes. No elements were established in the center of the specimen ((b) in A-A section view) which represented the debonding area between carbon panel and honeycomb, and general contact was generally established to prevent the inner and outer skins from embedding each other. The bottom of the honeycomb structure was fixed, and constraints in the z direction were set in the remaining three sides. Based on the buckle module of the ABAQUS/Standard solver, the concentrated force in the Y-direction was applied to the top of the honeycomb panel to analyze the structural buckling. At the same time, keywords in the input file were modified to output the node information of all modes. Based on the Static General module, the concentrated load was replaced by the displacement load, and the keywords of the input file were modified again. The node information obtained from the buckling analysis, as an initial defect, was introduced to the post-buckling analysis model.

5. Numerical Results and Discussion

5.1. Effect of Imperfections on Stability

According to the method described by Tafreshi and Oswald [40], the instability modes of contained global buckling, mixed buckling and local buckling. As shown in Figure 15, the out-of-plane displacement of the red region was larger than that of the blue region. When D was 0, the instability mode of the defect-free model was global buckling. When D is 40 mm or 80 mm, the local buckling of the panel occurred in the debonding area. The larger the proportion of the structure occupied by the defect was, the more obvious the local instability mode of the structure was. It is worth noting that the small damages or defects occurring in specimens during the molding, manufacturing, and transportation processes may lead to asymmetry in the structure during the test process, and these factors were not considered in the linear buckling analysis.
As given in Table 8, the experimental results of the first-order buckling load of the two types of defective honeycomb panels were basically consistent with the simulated values. However, the simulated value of the defect-free honeycomb panel was higher than its failure load, this is because the material failure was not considered during the linear buckling analysis. The sloped area was just a weak link of the edge closed structure, so before the global buckling of the structure occurred, the sloped area was damaged, which resulted in losing bearing capacity in advance. Overall, further nonlinear post-buckling analysis is required in simulation in order to accurately assess the final bearing capacity of the structure.

5.2. Effect of Imperfections on Stability

The particularity of the material molding process made the structure itself have certain initial defects. After the buckling deformation, the structure showed obvious nonlinearity. Therefore, the nonlinear analysis of sandwich panels needed to introduce initial geometric defects to accurately obtain the final failure load of the structure. Considering that the first-order buckling mode and buckling load obtained by simulation were close to the test results, the node displacement of first-order buckling was introduced into the post-buckling model of the honeycomb panel. As given in Table 9, the ultimate loads of the honeycomb panel obtained by the numerical analysis method in this paper were in good agreement with the experimental values, and the maximum relative error did not exceed 10%. Compared with CD0, the carrying capacities of the CD40 and CD80 structures decreased by 9.9% and 22.1%, respectively.

5.3. Effect of Imperfections on Honeycomb Failure

As shown in Figure 16, the red area represents failed elements, and the blue represents non-failed elements. In the model CD0, the additional bending moment generated by the asymmetric load caused the structure to buckle along the inner skin. Observed from the sides of inner skin or outer skin, the honeycomb damage of the sloped area that was the weak link appeared first. As the buckling displacement of the structure increased, the damage expanded from the periphery to the center of the structure. Compared to the inner skin, the outer skin was closer to the loading plane of the edge closed structure, which resulted in more serious damage.
As shown in Figure 17, the sloped area and the debonding area were the weak areas of CD40 model. Because the diameter of the debonding area was still small, the honeycomb damage in the defect area did not appear before that of the sloped area. Observed from the sides of inner skin or outer skin, the honeycomb in the sloped area bore carrying more stress relatively, so damage appeared there first. As the degree of buckling increased, the damage expanded toward the center of the structure.
Owing to the large debonding diameter in the model CD80, local buckling appeared obviously in the debonding area of the skin, as shown in Figure 18. Observed from the inner skin, the honeycomb at the edge of the defect area carried more stress relatively, so damage appeared there first. As the degree of buckling increased, the damage extends to the sloped area in the vertical loading direction. Observed from the outer skin, the honeycomb of the inclined area appeared damaged first, and then the damage expanded along the vertical loading direction.

5.4. Effect of Imperfections on Skin Failure

As shown in Figure 19, Figure 20 and Figure 21, the clouds of damage modes and damage states of the skins were obtained by different angles and various types of models. In the model CD0 and CD40, the honeycomb in bevel area was subject to more stress under the compressive load, which resulted in its damage at first, but at this moment, the composite skin was not damaged due to its high strength. As the degree of buckling increased, the damage area of the honeycomb core gradually expanded. After the support of the honeycomb core was lost between inner and outer skins, the bending stiffness of the structure suddenly decreased, and the local stress and strain were too large, which caused the damage of fibers and matrix. As the damage of the honeycomb core and the composite skin expanded, the structure lost its bearing capacity. Because the skin of the model CD80 did not fail during loading, there was no failure cloud.

6. Conclusions

Finite element analysis and experimental verification are applied in this paper to study the failure behavior of edge-closed honeycomb panels with non-defect and debonding defects. Significant conclusions from the study are three-fold:
(1)
For specimens with face/core debond, the buckling load, ultimate load and failure mode obtained by finite element analysis are in accordance with the test results, which verifies the reliability of the finite element model and failure theory of honeycomb panels.
(2)
Edge-closed honeycomb panels with non-defect usually occur panel damage or honeycomb damage in the sloped area and closed area. However, for the ones with face/core debond (e.g., CD40), the debonding area becomes the main damaged area. When the debonding diameter is larger (e.g., CD80), the serious local buckling causes that the structure loses its bearing capacity without obvious damage to the panel.
(3)
With the increase of the defective size, the stability and bearing capacity of the structure show a downward trend. It can be seen that the interface debonding defect will seriously affect the mechanical properties of edge-closed composite honeycomb panels, so appropriate repair methods should be adopted in time to repair the defective structure.

Author Contributions

Conceptualization, B.Y. and S.P.; methodology, B.Y.; software, S.P. and J.Y.; validation, B.Y., S.P. and J.Y.; formal analysis, B.Y.; investigation, B.Y. and S.P.; resources, B.Y. and F.L.; data curation, S.P.; writing—original draft preparation, B.Y. and S.P.; writing—review and editing B.Y., J.Y., F.L., X.W. and M.T.; visualization, S.P.; supervision, M.T.; project administration, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to gratefully acknowledge Nanjing University of Aeronautics and Astronautics Unmanned Air Vehicle Institute as the supplier of the test specimens.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Harris, C.E.; Starnes, J.H., Jr.; Shuart, M.J. Design and manufacturing of aerospace composite structures, state-of-the-art assessment. J. Aircr. 2002, 39, 545–560. [Google Scholar] [CrossRef]
  2. Yu, K. Investigation on Compression Strength of Honeycomb Sandwich Structure with Closing Frame. Master’s Thesis, Nanchang University, Nanchang, China, 2019. [Google Scholar]
  3. Liu, Z.; Cheng, X.; Hu, R.; Wang, J. Compression properties of honeycomb sandwich plates with different shape delamination defects. Acta Mater. Compos. Sin. 2009, 26, 153–158. (In Chinese) [Google Scholar]
  4. Ratcliffe, J.; Krueger, R. Face sheet/core disbonding in sandwich composite components: A road map to standardization: Test method development. In Proceedings of the 11th International Conference on Sandwich Structures, Dania Beach, FL, USA, 20–22 May 2016. [Google Scholar]
  5. Saseendran, V.; Varatharaj, P.; Perera, S.; Seneviratne, W. Damage initiation and fracture analysis of honeycomb core single cantilever beam sandwich specimens. J. Sandw. Struct. Mater. 2020, in press. [Google Scholar] [CrossRef]
  6. Saseendran, V.; Carlsson, L.A.; Berggreen, C.; Seneviratne, W. Crack Length Correction and Root Rotation Angle in a Sandwich Single Cantilever Beam (SCB) Fracture Specimen. Int. J. Lightwght Mater. Manuf. 2020, 3, 426–434. [Google Scholar]
  7. Carlsson, L.A. On the design of the cracked sandwich beam (CSB) specimen. J. Reinf. Plast. Compos. 1991, 10, 434–444. [Google Scholar] [CrossRef]
  8. Zenkert, D. Poly (vinyl chloride) sandwich core materials: Fracture behaviour under mode II loading and mixed-mode conditions. Mater. Sci. Eng. A 1989, 108, 233–240. [Google Scholar] [CrossRef]
  9. Quispitupa, A.; Berggreen, C.; Carlsson, L.A. On the analysis of a mixed mode bending sandwich specimen for debond fracture characterization. Eng. Fract. Mech. 2009, 76, 594–613. [Google Scholar] [CrossRef] [Green Version]
  10. Quispitupa, A.; Berggreen, C.; Carlsson, L.A. Design analysis of the mixed mode bending sandwich specimen. J. Sandw. Struct. Mater. 2010, 12, 253–272. [Google Scholar] [CrossRef]
  11. Saseendran, V.; Berggreen, C. Mixed-mode fracture evaluation of aerospace grade honeycomb core sandwich specimens using the Double Cantilever Beam–Uneven Bending Moment test method. J. Sandw. Struct. Mater. 2020, 22, 991–1018. [Google Scholar] [CrossRef] [Green Version]
  12. Goswami, S.; Becker, W. The effect of facesheet/core delamination in sandwich structures under transverse loading. Compos. Struct. 2001, 54, 515–521. [Google Scholar] [CrossRef]
  13. Vadakke, V.; Carlsson, L.A. Experimental investigation of compression failure of sandwich specimens with face/core debond. Compos. Part. B Eng. 2004, 35, 583–590. [Google Scholar] [CrossRef]
  14. Aviles, F.; Carlsson, L.A. Experimental study of debonded sandwich panels under compressive loading. J. Sandw. Struct. Mater. 2006, 8, 7–31. [Google Scholar] [CrossRef]
  15. Hui, W.; Hang, J.; Xie, Z. Study on lateral compressive strength of honeycomb sandwich structure with face/core debonding defects. In Proceedings of the 3rd China International Conference on Composite Materials Science and Technology, Hangzhou, China, 21 October 2017. [Google Scholar]
  16. Nieh, Y.-H.; Lai, Y.-K.; Tsai, J.-L. Characterizing failure behaviors of debonded composite sandwich structures under compressive loading. J. Sandw. Struct. Mater. 2015, 17, 578–594. [Google Scholar] [CrossRef]
  17. Pietrek, M.; Horst, P. Analysis and numerical prediction of the delamination behavior of debonded asymmetric sandwich shells with a thin-walled skin considering plastic deformation. Compos. Struct. 2018, 188, 220–232. [Google Scholar] [CrossRef]
  18. Frostig, Y.; Sokolinsky, V. Higher-order buckling of debonded (delaminated) sandwich panels with soft core. AIAA J. 2000, 38, 2147–2159. [Google Scholar] [CrossRef]
  19. El-Sayed, S.; Sridharan, S. Cohesive layer models for predicting delamination growth and crack kinking in sandwich structures. Int. J. Fract. 2002, 117, 63–84. [Google Scholar] [CrossRef]
  20. Berggreen, C.; Simonsen, B.C. Non-uniform compressive strength of debonded sandwich panels–II. Fracture mechanics investigation. J. Sandw. Struct. Mater. 2005, 7, 483–517. [Google Scholar] [CrossRef]
  21. Moslemian, R.; Quispitupa, A.; Berggreen, C.; Hayman, B. Failure of uniformly compression loaded debond damaged sandwich panels—An experimental and numerical study. J. Sandw. Struct. Mater. 2012, 14, 297–324. [Google Scholar] [CrossRef]
  22. Sayyidmousavi, A.; Malekzadeh, K.; Bougharara, H. Finite element buckling analysis of laminated composite sandwich panels with transversely flexible core containing a face/core debond. J. Compos. Mater. 2012, 46, 193–202. [Google Scholar] [CrossRef]
  23. Feng, C.; Zhao, Q. Analysis of typical failure modes in beveled area of honeycomb sandwich structure. China Sci. Technol. Inf. 2017, 10, 33–34. (In Chinese) [Google Scholar]
  24. Castanié, B.; Barrau, J.J.; Jaouen, J.P.; Rivallant, S. Combined shear/compression structural testing of asymmetric sandwich structures. Exp. Mech. 2004, 44, 461–472. [Google Scholar] [CrossRef] [Green Version]
  25. Pan, S.; Wang, X.; Chen, X. Stability Research of Edge Closed Honeycomb Sandwich Panels with Face/ Core Debond Under Compression. J. Nanjing Univ. Aeronaut. Astronaut. 2019, 51, 35–40. (In Chinese) [Google Scholar]
  26. Deng, J.; Peng, A.; Chen, W.; Zhou, G.; Wang, X. On stability and damage behavior of asymmetric sandwich panels under uniaxial compression. J. Sandw. Struct. Mater. 2020, in press. [Google Scholar] [CrossRef]
  27. Zhang, T.; Ding, Y.; Jin, H. Comparative analysis of equivalent models for honeycomb sandwich plates. Chin. J. Appl. Mech. 2011, 28, 275–282. (In Chinese) [Google Scholar]
  28. Niu, M.C.Y. Composite Airframe Structures; Hong Kong Conmilit Press Ltd.: Hong Kong, China, 2005. [Google Scholar]
  29. Wang, Y.; Zhu, S.; Tong, M. Three Dimensional Modelling of Progressive Failure Analysis in Fibre-Reinforced Composite Laminates. In Proceedings of the Presented in the 50th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Palm Springs, CA, USA, 4–7 May 2009. [Google Scholar]
  30. Dassault SIMULIA. ABAQUS 6.14 Documentation. Providence, R. 2014. Available online: http://ivt-abaqusdoc.ivt.ntnu.no:2080/v6.14/index.html (accessed on 4 September 2020).
  31. Hashin, Z. Failure Criteria for Unidirectional Fiber Composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  32. Dugdale, D.S. Yielding of steel sheets containing slits. J. Mech. Phys. Solids 1960, 8, 100–104. [Google Scholar] [CrossRef]
  33. Barenblatt, G.I. The mathematical theory of equilibrium cracks in brittle fracture. Adv. Appl. Mech. 1962, 7, 55–129. [Google Scholar]
  34. Bak, B.V.; Turon, A.; Lindgaard, E.; Lund, E. A simulation method for high-cycle fatiguedriven delamination using a cohesive zone model. Int. J. Numer. Methods Eng. 2016, 106, 163–191. [Google Scholar] [CrossRef]
  35. Deng, J.; Zhou, G.; Yin, Q.; Xiang, C.; Cai, D. Progressive damage analysis of double-strap bonding repaired laminates under buckling compression. Acta Aeronaut. Astronaut. Sin. 2016, 37, 1526–1535. (In Chinese) [Google Scholar]
  36. Camanho, P.P.; Davila, C.; De Moura, M. Numerical simulation of mixed-mode progressive delamination in composite materials. J. Compos. Mater. 2003, 37, 1415–1438. [Google Scholar] [CrossRef]
  37. Benzeggagh, M.L.; Kenane, M. Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus. Compos. Sci. Technol. 1996, 56, 439–449. [Google Scholar] [CrossRef]
  38. Besant, T.; Davies, G.; Hitchings, D. Finite element modelling of low velocity impact of composite sandwich panels. Compos. Part. A Appl. Sci. 2001, 32, 1189–1196. [Google Scholar] [CrossRef]
  39. Zhou, H. Numerical Analysis on Damage and Repair of Honeycomb Sandwich Composite Panels. Master’s Thesis, Harbin Institute of Technology, Harbin, China, 2017. [Google Scholar]
  40. Tafreshi, A.; Oswald, T. Global buckling behaviour and local damage propagation in composite plates with embedded delaminations. Int. J. Pres. Ves. Pip. 2003, 80, 9–20. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The configuration of specimens.
Figure 1. The configuration of specimens.
Applsci 10 07457 g001
Figure 2. The influence of debonding sizes on buckling load.
Figure 2. The influence of debonding sizes on buckling load.
Applsci 10 07457 g002
Figure 3. Strain arrangement on both sides of specimens.
Figure 3. Strain arrangement on both sides of specimens.
Applsci 10 07457 g003
Figure 4. Fixture of tests.
Figure 4. Fixture of tests.
Applsci 10 07457 g004
Figure 5. Typical strain-load curves of the third non-defective compressive sample (CD0-3).
Figure 5. Typical strain-load curves of the third non-defective compressive sample (CD0-3).
Applsci 10 07457 g005
Figure 6. Typical strain-load curves of the third compressive sample with defective diameter 40 mm (CD40-3).
Figure 6. Typical strain-load curves of the third compressive sample with defective diameter 40 mm (CD40-3).
Applsci 10 07457 g006
Figure 7. Typical strain-load curves of the third compressive sample with defective diameter 80 mm (CD80-3).
Figure 7. Typical strain-load curves of the third compressive sample with defective diameter 80 mm (CD80-3).
Applsci 10 07457 g007
Figure 8. Typical strain-load curves of specimens with different defective.
Figure 8. Typical strain-load curves of specimens with different defective.
Applsci 10 07457 g008
Figure 9. The load-displacement curves of all types of specimens.
Figure 9. The load-displacement curves of all types of specimens.
Applsci 10 07457 g009
Figure 10. Compressive failure mode of CD0-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Figure 10. Compressive failure mode of CD0-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Applsci 10 07457 g010
Figure 11. Compressive failure mode of CD40-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Figure 11. Compressive failure mode of CD40-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Applsci 10 07457 g011
Figure 12. Compressive failure mode of CD80-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Figure 12. Compressive failure mode of CD80-3. (a) Failure mode of inner skin; (b) Failure mode of outer skin.
Applsci 10 07457 g012
Figure 13. The progressive failure process for composite honeycomb panels.
Figure 13. The progressive failure process for composite honeycomb panels.
Applsci 10 07457 g013
Figure 14. The grid discretization and boundary conditions of the debonding model.
Figure 14. The grid discretization and boundary conditions of the debonding model.
Applsci 10 07457 g014
Figure 15. Buckling modes of different honeycomb panels. (a) the first-order and second-order buckling modes of non-defective specimens; (b) the first-order and second-order buckling modes of specimens with defective diameter 40 mm; (c) the first-order and second-order buckling modes of specimens with defective diameter 80 mm.
Figure 15. Buckling modes of different honeycomb panels. (a) the first-order and second-order buckling modes of non-defective specimens; (b) the first-order and second-order buckling modes of specimens with defective diameter 40 mm; (c) the first-order and second-order buckling modes of specimens with defective diameter 80 mm.
Applsci 10 07457 g015
Figure 16. Honeycomb damage propagation in the CD0 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Figure 16. Honeycomb damage propagation in the CD0 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Applsci 10 07457 g016
Figure 17. Honeycomb damage propagation in the CD40 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Figure 17. Honeycomb damage propagation in the CD40 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Applsci 10 07457 g017
Figure 18. Honeycomb damage propagation in the CD80 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Figure 18. Honeycomb damage propagation in the CD80 Model. (a) the damage propagation of inner skin; (b) the damage propagation of outer skin.
Applsci 10 07457 g018
Figure 19. Failure modes of inner skin in the CD0 Model. (a) Matrix compressive failure (MCF); (b) Fiber compressive failure (FCF); (c) Matrix tensile failure (MTF); (d) Fiber tensile failure (FTF).
Figure 19. Failure modes of inner skin in the CD0 Model. (a) Matrix compressive failure (MCF); (b) Fiber compressive failure (FCF); (c) Matrix tensile failure (MTF); (d) Fiber tensile failure (FTF).
Applsci 10 07457 g019
Figure 20. Failure modes of outer skin in the CD0 Model. (a) Matrix compressive failure (MCF); (b) Fiber compressive failure (FCF); (c) Fiber tensile failure (FTF); (d) Matrix tensile failure (MTF).
Figure 20. Failure modes of outer skin in the CD0 Model. (a) Matrix compressive failure (MCF); (b) Fiber compressive failure (FCF); (c) Fiber tensile failure (FTF); (d) Matrix tensile failure (MTF).
Applsci 10 07457 g020
Figure 21. Failure modes of outer skin in the CD40 Model. (a) Fiber compressive failure (FCF); (b) Matrix compressive failure (MCF).
Figure 21. Failure modes of outer skin in the CD40 Model. (a) Fiber compressive failure (FCF); (b) Matrix compressive failure (MCF).
Applsci 10 07457 g021
Table 1. The mechanical parameters of plain carbon fabric prepreg.
Table 1. The mechanical parameters of plain carbon fabric prepreg.
E1 (GPa)E2 (GPa)G12 (GPa)υ12XT = XC (MPa)YT = YC (MPa)S12 (MPa)
49493.60.044676.2615.5261.5
Table 2. Equivalent mechanical properties of Nomex honeycomb.
Table 2. Equivalent mechanical properties of Nomex honeycomb.
E1 (MPa)E2 (MPa)E3 (MPa)G12 (MPa)G13 (MPa)G23 (MPa)υ12υ13 = υ23
0.19530.19531360.07340230.940.0003
XTT (MPa)SLT (MPa)SWT(MPa)
1.931.260.62
Table 3. Mechanical properties of adhesive.
Table 3. Mechanical properties of adhesive.
E (GPa)G (GPa)υtn (MPa)ts = tt (MPa)Gn (N·mm−1)Gs = Gt (N·mm−1)
3.21.230.358650.7443.816
Table 4. The design basis of specimens.
Table 4. The design basis of specimens.
Design ParametersReference
The minimum number of skin layers is 2Prevent moisture penetration and meet the minimum structural strength requirements
The number of skin layers reaches 5Avoid damage caused by low energy impact
The empirical height value of Non-metal core is about 12.7 mmProvides appropriate shear strength and stiffness
±45 plyEnhanced in-plane shear
0/90 plyStrengthening structural bending
Beveling angle is generally 20~30°Minimize the bonding pressure on the inclined surface
Table 5. Test matrix.
Table 5. Test matrix.
TypeDefective Diameter (mm)Number of Specimens
Non-defective specimens03
Debonding specimens403
803
Table 6. The ultimate loads of various specimens.
Table 6. The ultimate loads of various specimens.
NumberUltimate Loads (kN)Mean Value (kN)CV
123
CD097.698.392.396.011.25%
CD4086.585.487.686.51.39%
CD8074.476.972.974.75.46%
Table 7. Equivalent stiffness reduction factor of honeycomb.
Table 7. Equivalent stiffness reduction factor of honeycomb.
Damaged Status λ 11 λ 22 λ 33 λ 12 λ 13 λ 23
Intact111111
Damaged0.40.40.450.40.40.4
Fold0.20.20.20.20.20.2
Table 8. Buckling loads of various honeycomb panels.
Table 8. Buckling loads of various honeycomb panels.
NumberMean Value of Experimental Results (kN)Simulated Value (kN)Error/%
CD06510866.2
CD40606813.3
CD803532−8.6
Table 9. Fracture loads of various honeycomb panels.
Table 9. Fracture loads of various honeycomb panels.
NumberMean Value of Experimental Results (kN)Simulated Value (kN)Error/%
CD096.090.7−5.5
CD4086.582.2−4.9
CD8074.769.3−7.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yan, B.; Wang, X.; Pan, S.; Tong, M.; Yu, J.; Liu, F. Stability and Failure of the Edge-Closed Honeycomb Sandwich Panels with Face/Core Debonding. Appl. Sci. 2020, 10, 7457. https://doi.org/10.3390/app10217457

AMA Style

Yan B, Wang X, Pan S, Tong M, Yu J, Liu F. Stability and Failure of the Edge-Closed Honeycomb Sandwich Panels with Face/Core Debonding. Applied Sciences. 2020; 10(21):7457. https://doi.org/10.3390/app10217457

Chicago/Turabian Style

Yan, Bing, Xinfeng Wang, Song Pan, Mingbo Tong, Jian Yu, and Fan Liu. 2020. "Stability and Failure of the Edge-Closed Honeycomb Sandwich Panels with Face/Core Debonding" Applied Sciences 10, no. 21: 7457. https://doi.org/10.3390/app10217457

APA Style

Yan, B., Wang, X., Pan, S., Tong, M., Yu, J., & Liu, F. (2020). Stability and Failure of the Edge-Closed Honeycomb Sandwich Panels with Face/Core Debonding. Applied Sciences, 10(21), 7457. https://doi.org/10.3390/app10217457

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