Next Article in Journal
A Review of Current Research on the Use of Geopolymer Recycled Aggregate Concrete for Structural Members
Next Article in Special Issue
Crashworthiness Analysis of Thin-Walled Square Columns with a Hole Trigger
Previous Article in Journal
Synthesis of ((CeO2)0.8(Sm2O3)0.2)@NiO Core-Shell Type Nanostructures and Microextrusion Printing of a Composite Anode Based on Them
Previous Article in Special Issue
Load Eccentricity of Compressed Composite Z-Columns in Non-Linear State
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of Thin-Walled Perforated Composite Columns Using Finite Element Method

by
Katarzyna Falkowicz
Department of Machine Design and Mechatronics, Faculty of Mechanical Engineering, Lublin University of Technology, Nadbystrzycka 36, 20-618 Lublin, Poland
Materials 2022, 15(24), 8919; https://doi.org/10.3390/ma15248919
Submission received: 3 November 2022 / Revised: 24 November 2022 / Accepted: 12 December 2022 / Published: 13 December 2022
(This article belongs to the Special Issue Compressive Behavior of Materials and Structures)

Abstract

:
Open holes or cut-outs have been commonly used in composite structures for various engineering purposes. Those elements often demand perforation especially for weight reduction and to ease maintenance and servicing operations, for example, in aircraft wing ribs. This work presents a numerical study of the stability behavior of composite perforated columns subjected to a compressive load. Profiles were made of CFRP laminate and weakened by three types of cut-out. Four parameters, spacing ratio S/D0, opening ratio D/D0, hole shape and arrangement of layers, were selected to check their effect on the buckling load and postbuckling behavior of the tested channel profiles. To carry out the numerical analysis, the Abaqus software was used. The results obtained during the analysis helped to identify the best combination of tested parameters to obtain the highest critical load. The performed analysis show that the columns’ behavior is sensitive to configuration of composite, opening ratio and hole shape.

1. Introduction

Thin-walled structures belong to the load-carrying elements characterized by high stiffness and strength while maintaining low weight [1,2,3,4] and that allows designers to have freedom in shaping the construction form. Thanks to the mentioned properties, these kinds of structures have wide application, where the low weight is important, i.e., aerospace and automotive sectors and civil engineering [1,5,6]. Unfortunately, as in any thin-walled elements, they have important disadvantages such as loss of stability when they are compressed or sheared [7,8,9,10]. However, they have one of the most important features which is the possibility of working in the postcritical state when the buckling of a thin-walled element is local and elastic [11,12,13,14,15,16]. In the worldwide available literature, more studies about linear and nonlinear problems of thin-walled structures’ stability can be found [17,18,19]. Moreover, thin-walled elements often demand perforation for weight reduction and to ease servicing and maintenance operations, for example, in aircraft wing ribs. These perforations cause a redistribution of stresses in the structure that may change the strength of the structural element and the elastic stiffness. The buckling behavior of structural perforated profiles is significantly influenced by the size, shape, location or number of perforations. Previous investigations on composite plates [20,21,22,23,24] have shown the significant effect of the geometric parameters and shape of holes, as well as the composite configuration, on the stability behavior. Therefore, it seems imperative to investigate the influence of the abovementioned parameters on thin-walled perforated profiles’ buckling behavior.
Thin-walled channel profiles are currently widely used and often weakened by holes to reduce their volume. An example of such elements are perforated profiles used, among others, for shelves, balustrades, etc. In Figure 1, an example of popular design solutions based on the use of openwork beams as load-carrying structures is shown.
In the literature, there are a lot of investigations about the design and behavior of columns with different kinds of cross-sections made of composite [4,25,26] or traditional materials [27,28]. Moreover, in the available literature there are some examinations where the authors tested elements with different kinds of hole shapes [29,30,31]. However, the obtained results showed that most studies of the literature concern structures without holes, or with holes, but made of traditional materials. However, the knowledge about buckling behavior of perforated composite columns is insufficient. Furthermore, there are no specific standards for composite perforated profiles such as in case of perforated profilers made of steel or aluminum.
The load-bearing capacity of composite profiles, due to wide usage of composite elements in a thin plate form, is very important. Furthermore, different kinds of cut-outs are very often made as a part of the project. This makes the information about the buckling and postbuckling behavior of composite profiles weakened by holes necessary.
The idea of this work is based on the study carried out by Khazaal et al. [32], who tested aluminum alloy profiles subjected to compression. To optimize the results and to check which parameter has the most influence on the buckling behavior of the aluminum thin-walled elements, they investigated three shapes of hole: rectangular, circular and hexagonal. A very similar investigation was described in paper [33], where the authors took under consideration profiles made from GFRP composite. It is worth adding here that, in previous articles [34,35], the authors tested composite profiles with different kinds of cross-section but without perforation. Therefore, the abovementioned studies together with [32,33] have ignited motivation for the presented study, hence contributing towards the knowledge on thin-walled composite structures with holes.
This paper includes the linear and nonlinear analysis of compressed, thin-walled, perforated profiles with a channel cross-section. Numerical calculations were carried out with the commercial Abaqus program using the finite element method (FEM). The effect of parameters such as hole shapes, spacing ratio, opening ratio and arrangement of layers’ layout on the buckling and postbuckling behavior of the perforated channel cross-section profiles was investigated. Moreover, obtained results were compared with unperforated profiles which were tested in previous works [9,35,36]. It should be added here that the tested columns were made of CFRP composite material, and previous research has been carried out on unperforated profiles or on perforated profiles but made from classic materials, including paper [31], where the authors tested perforated profiles made of GFRP composite. It should also be mentioned here that there are no specific standards for the perforation of composite profiles. Furthermore, columns made of CFRP composite have not been taken under this kind of consideration. Therefore, it was a good motivation for conducting the analysis presented in this paper.

2. Research Subject and Methodology

The research subjects were thin-walled composite profiles with a C cross-section subjected to axial compression. The geometric parameters of tested profiles are shown in Figure 2 (length—250 mm, the web—60 mm, the profile wall h = 30 mm). The tested columns consisted of 8 plies symmetric to the laminate mid-plane. Each ply of profiles had a thickness of 0.105 mm and the total thickness of the wall was 0.84 mm. The profiles were made of CFRP composite, whose mechanical properties are presented in Table 1.
According to the previous research [7,36,37], four configurations of composite were chosen for the analysis: P1: [0/45/−45/90]s, P2: [45/−45/90/0]s, P3: [90/−45/45/0]s, P4: [90/0/90/0]s. The profiles were compared with profiles without cut-outs, which were also tested experimentally and numerically in previous works [34,37,38]. The exemplary obtained buckling and postcritical deformation forms are shown in Figure 3.
The results presented above became a reference point for the research performed and presented in this paper. The current paper focuses on the behavior of channel cross-section profiles weakened by holes and subjected to axial compression. This work does not concern the in-depth presentation of the experimental results, it focuses on the numerical analysis, but this is the start for further experimental research.
In this work, according to [32,33], three shapes of holes (circular, square, hexagonal) made in the profile web were designed (Figure 4). The overall geometric parameters were constant. In [32], the authors tested aluminum profiles where the geometry of web holes must be within the given ranges of Eurocode (Equations (1) and (2)) in order to prevent any unwanted failures such as cracks between holes, and to achieve a maximum possible reduction in weight.
1.25 < D/D0 < 1.75
1.08 < S/D0 < 1.5
where: D—width of profile web; D0—diameter/width of hole; S—distance between holes; D/D0—opening ratio; S/D0—spacing ratio.
For composite profiles, this kind of Eurocode does not exist but it was a good reason to take those equations under consideration and perform numerical analysis for composite columns with a channel section. The above equations were used in [33] to analyze perforated profiles made of GFRP composite.
Three parameters, opening ratio (D/D0), spacing ratio (S/D0) and hole shape, were taken under consideration, and their influence on buckling and postbuckling behavior was measured. The opening ratio (D/D0) and spacing ratio (S/D0) were selected according to Equations (1) and (2). Parameters which were outside the bounds of the above equations were also taken into consideration. The parameters and levels are presented in Table 2.

3. Numerical Analysis

The numerical calculations were carried out in the Abaqus program with the finite element method, which has a very wide range of applications [39,40]. Analysis was performed in two stages. The first stage was the linear stability analysis, based on the solution of the following generalized eigenvalue problem:
|[K] + λi [H]| = 0
where: [K]—structural stiffness matrix, λi—the i-th eigenvalue, [H]—stress stiffness matrix.
The critical state of the structure was described by a linear model using the minimum total potential energy principle. Equation (3) presents mathematical notation of the loss of stability phenomenon. The scope of numerical simulations involved describing the critical state of the structure and influence of shape and parameters of cut-outs on critical force.
The second stage of analysis tackled the nonlinear stability problem, based on the incremental iterative Newton–Raphson method. This analysis allowed us to determine the equilibrium path for the structure. The calculations were performed until failure initiation of first layer according to the Tsai–Wu criterion [41].
The discrete model was prepared by using 8-node shell elements (S8R) of second order with reduced integration. The structure of the composite material was defined depending the thickness of the finite element. Mesh density, which was used for the composite profiles, equaled 4 mm and was selected as the optimal density according to the previous analyses on unperforated profiles with channel cross-sections [32]. The value of buckling load obtained for mesh size 4 mm was the closest to the buckling load obtained by experiment and the percentage error was around 2%. However, additionally, the convergence was studied for three different mesh sizes: 3 mm, 4 mm and 5 mm. The results of tests for profiles with P2 lay-up are presented in Figure 5. Obtained results showed that the mesh size has no significant influence on buckling behavior. Furthermore, denser mesh lengthens the processing time. Similar conclusions were presented in [34,40,42].
The boundary conditions of the tested model are presented in Figure 6. The model consisted of two rigid plates as a support. The upper plate had the ability to move along the Z-axis (compressive force direction) and the lower plate was fully fixed. The boundary conditions were defined in a created reference point (RP), connected to plates. Between the plates and composite profile, a contact interaction was introduced with a friction coefficient of 0.2.

4. Results and Discussion

In this section, the results from the performed numerical analysis are presented. Table 3 shows the values of critical forces for all tested cases with different kinds of holes shape. The below results present the buckling load for the P1 configuration of the composite. Ratios which exceed ranges of Equations (2) and (3) are marked in red. In the last column is the percentage difference between critical force for unperforated profiles and for perforated tested profiles. In Figure 7, for better visualization, the effect of spacing and opening ratio on the buckling load for three different hole shapes is presented.
The presented results from FE analysis show that the introduction of perforation leads to the load decrease. The circular holes caused a smaller decrease of critical load and for all tested cases it was in the range of 17 ÷ 28.48%. For square holes, this range was 21.30 ÷ 35.06% and for the hexagonal holes it was 18.76 ÷ 31.54%. It can be observed that when increasing the opening ratio, the buckling load also increased (Figure 7a), whereas the effect of the spacing ratio (Figure 7b) did not show any significant influence on the value of critical load. Those observations confirm the results obtained in [33] for GFRP composite profiles.
It can be observed that among all tested samples with different kinds of hole shapes, the circular holes had the greatest impact on the buckling load with opening ratio D/D0 = 2 and with spacing ratio S/D0 = 1.67. In [33], the best ultimate strength was also seen for the circular shape, while in [32] the best results were for the hexagonal shape.
In this section, the results of the lay-up effect on the buckling behavior of the tested profiles are presented. Table 4 shows the obtained values of critical loads for perforated and unperforated profiles, for four different composite lay-ups. For this analysis, the profile with a circular hole shape with the best results of parameters, which were tested in the first part of the analysis, was chosen. Figure 8 additionally presents the same results but in the form of a column chart.
It can be observed that the laminate lay-up has a significant effect on buckling load. Furthermore, the weakest configuration is P4 and the difference between buckling load for unperforated and perforated profiles for this configuration is more than 20%, compared to the other three configurations where the difference was less than 20%. The highest critical load was obtained for the P2 lay-up. This is ascribable to the fact that, for P2, most outer layers are at an angle of 45 degrees, whereas for P4 the lay-up consists of a few layers at an angle of 90 degrees, which reduces the overall structural stiffness. In all cases, the shape of the buckling mode maintained its symmetrical character with respect to the symmetry plane of the C-profile (Figure 9). However, for almost all lay-ups, the local buckling of the web and shelf was characterized by different numbers of half-waves. The number of half-waves for the perforated profiles was lower by one or two compared to the unperforated profiles.
The place of maximum deformation was also different: for the P1 and P2 configurations, it was on the shelf closest to the rigid plate and for the P3 and P4 configuration it was in the middle of the profile shelf.
In the second part of the research, the nonlinear stability of the structure in the postbuckling range was measured. The aim of the analysis was to investigate the profile behavior in the postcritical state with the implemented lower buckling mode, which was determined in the first part of the tests. During the loading of the structure, a deepening form of laminate buckling was observed. The load-carrying capacity of tested profiles was determined by the Tsai–Wu failure criterion. This criterion was successfully used to determine the failure load of a compressed thin-walled composite plate element weakened by holes and columns with complex cross-sectional shapes in [21,43,44]. Figure 10 shows a designated example map of the Tsai–Wu criterion, where areas in which damage initiation of the first composite layer may occur, corresponding to the achievement of the value of 1 critical parameter, determined according to the above criterion, are presented. Moreover, in Figure 11, failure maps for perforated and for unperforated profiles for the chosen P2 configuration are compared.
It can be observed that the most critical area is the corner near to the end of the cross-section. Moreover, the introduction of perforation did not change this localization. Figure 12 additionally presents the obtained equilibrium paths for the perforated profiles for the four considered lay-ups. Moreover, in Figure 13, the equilibrium paths are compared for the perforated and unperforated profiles with P2 lay-up.
The obtained results show that for all of the tested profiles, the postbuckling equilibrium paths are stable. In addition, the postbuckling behavior of perforated and unperforated profiles is similar at the beginning. However, with the deepening of the deflection, the differences between the equilibrium paths increase. The highest deflection at failure was observed for the P4 lay-up and the lowest one for the P2 lay-up. Furthermore, the equilibrium path for the P2 laminate lay-up seems to be the one with the highest stiffness.

5. Conclusions

In this study, the linear and nonlinear analysis of perforated composite profiles with a C cross-section subjected to axial compression was performed. The effect of three hole shapes, opening ratio D/D0, the spacing ratio S/D0 and composite arrangement of layout on buckling behavior of perforated profiles was studied. The obtained results were compared with profiles without holes. On the basis of the obtained results, some conclusions can be drawn:
The introduction of perforations caused a decrease in the buckling loads: 21.30 ÷ 35.06% for the square shape, 17 ÷ 28.48% for the circular shape and 18.76 ÷ 31.54% for the hexagonal shape.
The obtained results show that the shape of holes, opening ratio and arrangement of laminate layers have the greatest impact on the value of buckling load, whereas the spacing ratio has no significant influence on the buckling load.
The circular holes with D/D0 = 2 and S/D0 = 1.67 gave the highest value of critical force.
The highest critical load was obtained for the P2 lay-up, while the lowest one for the P4 lay-up.
The introduced perforation caused not only a decrease in the critical load value, but also a change in the buckling form. The decrease in the critical load for all tested configurations was: 17.4% for P1, 15.27% for P2, 18.08% for P3 and 29.32% for P4.
The presented results allow us to expand our knowledge about the design of thin-walled composite structures weakened by holes with potential significance for practical applications. However, to confirm the obtained numerical results, a deeper analysis of nonlinear stability conditions and experimental validation are necessary.

Funding

The research was financed in the framework of the project of Lublin University of Technology—Regional Excellence Initiative, funded by the Polish Ministry of Science and Higher Education (contract no. 030/RID/2018/19).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Lublin University of Technology, Faculty of Mechanical Engineering.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Kopecki, T.; Mazurek, P.; Lis, T.; Chodorowska, D. Post-Buckling Deformation States of Semi-Monocoque Cylindrical Structures with Large Cut-Outs under Operating Load Conditions. Numerical Analysis and Experimental Tests. EiN 2016, 18, 16–24. [Google Scholar] [CrossRef]
  2. Kopecki, T.; Mazurek, P.; Lis, T. Experimental and Numerical Analysis of a Composite Thin-Walled Cylindrical Structures with Different Variants of Stiffeners, Subjected to Torsion. Materials 2019, 12, 3230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Baker, A.A.; Dutton, S.; Kelly, D. (Eds.) Composite Materials for Aircraft Structures, 2nd ed.; Education Series; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2004; ISBN 978-1-56347-540-5. [Google Scholar]
  4. Debski, H.; Rozylo, P.; Wysmulski, P.; Falkowicz, K.; Ferdynus, M. Experimental Study on the Effect of Eccentric Compressive Load on the Stability and Load-Carrying Capacity of Thin-Walled Composite Profiles. Compos. Part B Eng. 2021, 226, 109346. [Google Scholar] [CrossRef]
  5. Chróścielewski, J.; Miśkiewicz, M.; Pyrzowski, Ł.; Rucka, M.; Sobczyk, B.; Wilde, K. Modal Properties Identification of a Novel Sandwich Footbridge–Comparison of Measured Dynamic Response and FEA. Compos. Part B Eng. 2018, 151, 245–255. [Google Scholar] [CrossRef]
  6. Chróścielewski, J.; Miśkiewicz, M.; Pyrzowski, Ł.; Sobczyk, B.; Wilde, K. A Novel Sandwich Footbridge-Practical Application of Laminated Composites in Bridge Design and in Situ Measurements of Static Response. Compos. Part B Eng. 2017, 126, 153–161. [Google Scholar] [CrossRef]
  7. Falkowicz, K. Experimental and Numerical Analysis of Compression Thin-Walled Composite Plates Weakened by Cut-Outs. Arch. Civ. Eng. 2017, 63, 161–172. [Google Scholar] [CrossRef]
  8. Falkowicz, K.; Debski, H.; Teter, A. Design Solutions for Improving the Lowest Buckling Loads of a Thin Laminate Plate with Notch. In AIP Conference Proceedings; AIP Publishing LLC: Lublin, Poland, 2018; p. 080004. [Google Scholar]
  9. Wysmulski, P.; Debski, H. Post-Buckling and Limit States of Composite Channel-Section Profiles under Eccentric Compression. Compos. Struct. 2020, 245, 112356. [Google Scholar] [CrossRef]
  10. Rozylo, P.; Wysmulski, P. Failure Analysis of Thin-Walled Composite Profiles Subjected to Axial Compression Using Progressive Failure Analysis (PFA) and Cohesive Zone Model (CZM). Compos. Struct. 2021, 262, 113597. [Google Scholar] [CrossRef]
  11. Banat, D.; Mania, R.J. Damage Analysis of Thin-Walled GLARE Members under Axial Compression–Numerical and Experiment Investigations. Compos. Struct. 2020, 241, 112102. [Google Scholar] [CrossRef]
  12. Falkowicz, K.; Debski, H. The Post-Critical Behaviour of Compressed Plate with Non-Standard Play Orientation. Compos. Struct. 2020, 252, 112701. [Google Scholar] [CrossRef]
  13. Kubiak, T.; Samborski, S.; Teter, A. Experimental Investigation of Failure Process in Compressed Channel-Section GFRP Laminate Columns Assisted with the Acoustic Emission Method. Compos. Struct. 2015, 133, 921–929. [Google Scholar] [CrossRef]
  14. Madukauwa-David, I.D.; Drissi-Habti, M. Numerical Simulation of the Mechanical Behavior of a Large Smart Composite Platform under Static Loads. Compos. Part B Eng. 2016, 88, 19–25. [Google Scholar] [CrossRef]
  15. Rozylo, P. Experimental-Numerical Study into the Stability and Failure of Compressed Thin-Walled Composite Profiles Using Progressive Failure Analysis and Cohesive Zone Model. Compos. Struct. 2021, 257, 113303. [Google Scholar] [CrossRef]
  16. Ribeiro, M.L.; Vandepitte, D.; Tita, V. Damage Model and Progressive Failure Analyses for Filament Wound Composite Laminates. Appl. Compos. Mater. 2013, 20, 975–992. [Google Scholar] [CrossRef]
  17. Ovesy, H.R.; Loughlan, J.; GhannadPour, S.A.M. Geometric Non-Linear Analysis of Channel Sections under End Shortening, Using Different Versions of the Finite Strip Method. Comput. Struct. 2006, 84, 855–872. [Google Scholar] [CrossRef]
  18. Einafshar, N.; Lezgy-Nazargah, M.; Beheshti-Aval, S.B. Buckling, Post-Buckling and Geometrically Nonlinear Analysis of Thin-Walled Beams Using a Hypothetical Layered Composite Cross-Sectional Model. Acta Mech. 2021, 232, 2733–2750. [Google Scholar] [CrossRef]
  19. Rozylo, P.; Falkowicz, K. Stability and Failure Analysis of Compressed Thin-Walled Composite Structures with Central Cut-out, Using Three Advanced Independent Damage Models. Compos. Struct. 2021, 273, 114298. [Google Scholar] [CrossRef]
  20. Falkowicz, K.; Debski, H. The Work of a Compressed, Composite Plate in Asymmetrical Arrangement of Layers. In AIP Conference Proceedings; AIP Publishing LLC: Depok, Indonesia, 2019; p. 020005. [Google Scholar]
  21. Falkowicz, K.; Debski, H. Stability Analysis of Thin-Walled Composite Plate in Unsymmetrical Configuration Subjected to Axial Load. Thin-Walled Struct. 2021, 158, 107203. [Google Scholar] [CrossRef]
  22. Falkowicz, K. Composite Plate Analysis Made in an Unsymmetric Configuartion. J. Phys. Conf. Ser. 2021, 2130, 012014. [Google Scholar] [CrossRef]
  23. Falkowicz, K.; Ferdynus, M.; Rozylo, P. Experimental and Numerical Analysis of Stability and Failure of Compressed Composite Plates. Compos. Struct. 2021, 263, 113657. [Google Scholar] [CrossRef]
  24. Falkowicz, K. Experimental and Numerical Failure Analysis of Thin-Walled Composite Plates Using Progressive Failure Analysis. Compos. Struct. 2023, 305, 116474. [Google Scholar] [CrossRef]
  25. Falkowicz, K. Numerical Investigations of Perforated CFRP Z-Cross-Section Profiles, under Axial Compression. Materials 2022, 15, 6874. [Google Scholar] [CrossRef] [PubMed]
  26. Wysmulski, P. The Effect of Load Eccentricity on the Compressed CFRP Z-Shaped Columns in the Weak Post-Critical State. Compos. Struct. 2022, 301, 116184. [Google Scholar] [CrossRef]
  27. Yuan, H.X.; Wang, Y.Q.; Chang, T.; Du, X.X.; Bu, Y.D.; Shi, Y.J. Local Buckling and Postbuckling Strength of Extruded Aluminium Alloy Stub Columns with Slender I-Sections. Thin-Walled Struct. 2015, 90, 140–149. [Google Scholar] [CrossRef]
  28. Zhao, Y.; Zhai, X.; Sun, L. Test and Design Method for the Buckling Behaviors of 6082-T6 Aluminum Alloy Columns with Box-Type and L-Type Sections under Eccentric Compression. Thin-Walled Struct. 2016, 100, 62–80. [Google Scholar] [CrossRef]
  29. Akbulut, H.; Sayman, O. An Investigation on Buckling of Laminated Plates with Central Square Hole. J. Reinf. Plast. Compos. 2001, 20, 1112–1124. [Google Scholar] [CrossRef]
  30. Aydin Komur, M.; Sen, F.; Ataş, A.; Arslan, N. Buckling Analysis of Laminated Composite Plates with an Elliptical/Circular Cutout Using FEM. Adv. Eng. Softw. 2010, 41, 161–164. [Google Scholar] [CrossRef]
  31. Ouinas, D.; Achour, B. Buckling Analysis of Laminated Composite Plates [(θ/−θ)] Containing an Elliptical Notch. Compos. Part B Eng. 2013, 55, 575–579. [Google Scholar] [CrossRef]
  32. Khazaal, D.S.; AL-khafaji, H.M.; Abdulsahib, I.A. Buckling Behavior of Aluminum Alloy Thin-Walled Beam with Holes under Compression Loading. Jcoeng 2020, 26, 137–154. [Google Scholar] [CrossRef]
  33. Bin Kamarudin, M.N.; Mohamed Ali, J.S.; Aabid, A.; Ibrahim, Y.E. Buckling Analysis of a Thin-Walled Structure Using Finite Element Method and Design of Experiments. Aerospace 2022, 9, 541. [Google Scholar] [CrossRef]
  34. Wysmulski, P.; Debski, H.; Falkowicz, K. Sensitivity of Compressed Composite Channel Columns to Eccentric Loading. Materials 2022, 15, 6938. [Google Scholar] [CrossRef] [PubMed]
  35. Rozylo, P.; Falkowicz, K.; Wysmulski, P.; Debski, H.; Pasnik, J.; Kral, J. Experimental-Numerical Failure Analysis of Thin-Walled Composite Columns Using Advanced Damage Models. Materials 2021, 14, 1506. [Google Scholar] [CrossRef] [PubMed]
  36. Wysmulski, P.; Falkowicz, K.; Filipek, P. Buckling State Analysis of Compressed Composite Plates with Cut-Out. Compos. Struct. 2021, 274, 114345. [Google Scholar] [CrossRef]
  37. Wysmulski, P.; Debski, H.; Falkowicz, K. Stability Analysis of Laminate Profiles under Eccentric Load. Compos. Struct. 2020, 238, 111944. [Google Scholar] [CrossRef]
  38. Wysmulski, P.; Debski, H.; Falkowicz, K.; Rozylo, P. The Influence of Load Eccentricity on the Behavior of Thin-Walled Compressed Composite Structures. Compos. Struct. 2019, 213, 98–107. [Google Scholar] [CrossRef]
  39. Jonak, J.; Karpiński, R.; Wójcik, A. Influence of the Undercut Anchor Head Angle on the Propagation of the Failure Zone of the Rock Medium. Materials 2021, 14, 2371. [Google Scholar] [CrossRef] [PubMed]
  40. Jonak, J.; Karpiński, R.; Siegmund, M.; Wójcik, A.; Jonak, K. Analysis of the Rock Failure Cone Size Relative to the Group Effect from a Triangular Anchorage System. Materials 2020, 13, 4657. [Google Scholar] [CrossRef]
  41. Tsai, S.W.; Wu, E.M. A General Theory of Strength for Anisotropic Materials. J. Compos. Mater. 1971, 5, 58–80. [Google Scholar] [CrossRef]
  42. Jonak, J.; Karpiński, R.; Wójcik, A.; Siegmund, M. The Influence of the Physical-Mechanical Parameters of Rock on the Extent of the Initial Failure Zone under the Action of an Undercut Anchor. Materials 2021, 14, 1841. [Google Scholar] [CrossRef]
  43. Falkowicz, K.; Debski, H.; Wysmulski, P. Effect of Extension-Twisting and Extension-Bending Coupling on a Compressed Plate with a Cut-Out. Compos. Struct. 2020, 238, 111941. [Google Scholar] [CrossRef]
  44. Wysmulski, P. Load Eccentricity of Compressed Composite Z-Columns in Non-Linear State. Materials 2022, 15, 7631. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Example of application of perforated elements.
Figure 1. Example of application of perforated elements.
Materials 15 08919 g001
Figure 2. Geometric parameters of the tested profile.
Figure 2. Geometric parameters of the tested profile.
Materials 15 08919 g002
Figure 3. (a) Form of buckling, (b) form of postcritical deformation.
Figure 3. (a) Form of buckling, (b) form of postcritical deformation.
Materials 15 08919 g003
Figure 4. Shapes and dimensions of holes: (a) circle, (b) square, (c) hexagonal.
Figure 4. Shapes and dimensions of holes: (a) circle, (b) square, (c) hexagonal.
Materials 15 08919 g004
Figure 5. The results from convergence test: (a) mesh size = 3 mm, (b) mesh size = 4 mm, (c) mesh size = 5 mm.
Figure 5. The results from convergence test: (a) mesh size = 3 mm, (b) mesh size = 4 mm, (c) mesh size = 5 mm.
Materials 15 08919 g005
Figure 6. Boundary conditions of tested profile.
Figure 6. Boundary conditions of tested profile.
Materials 15 08919 g006
Figure 7. The effect of opening ratio (a) and spacing ratio (b) on buckling load for three different hole shapes.
Figure 7. The effect of opening ratio (a) and spacing ratio (b) on buckling load for three different hole shapes.
Materials 15 08919 g007
Figure 8. Effect of opening ratio and space ratio on critical buckling load for different hole shapes.
Figure 8. Effect of opening ratio and space ratio on critical buckling load for different hole shapes.
Materials 15 08919 g008
Figure 9. Buckling simulation results for the P2 lay-up: (a) unperforated profile, (b) profile with circular holes.
Figure 9. Buckling simulation results for the P2 lay-up: (a) unperforated profile, (b) profile with circular holes.
Materials 15 08919 g009
Figure 10. Tsai–Wu critical parameter map with the obtained form of loss of stability for the configurations: (a) P1, (b) P2, (c) P3, (d) P4.
Figure 10. Tsai–Wu critical parameter map with the obtained form of loss of stability for the configurations: (a) P1, (b) P2, (c) P3, (d) P4.
Materials 15 08919 g010
Figure 11. Tsai–Wu critical parameter map with the obtained form of loss of stability for the P2 configuration for: (a) unperforated, (b) perforated profile.
Figure 11. Tsai–Wu critical parameter map with the obtained form of loss of stability for the P2 configuration for: (a) unperforated, (b) perforated profile.
Materials 15 08919 g011
Figure 12. Postbuckling equilibrium paths for perforated profiles with all tested lay-ups.
Figure 12. Postbuckling equilibrium paths for perforated profiles with all tested lay-ups.
Materials 15 08919 g012
Figure 13. Postbuckling equilibrium paths for perforated profiles with P2 lay-up and for unperforated profiles with P2 lay-up (results from numerical and experimental analysis).
Figure 13. Postbuckling equilibrium paths for perforated profiles with P2 lay-up and for unperforated profiles with P2 lay-up (results from numerical and experimental analysis).
Materials 15 08919 g013
Table 1. Mechanical properties of a CFRP lamina.
Table 1. Mechanical properties of a CFRP lamina.
Young’s Modulus [MPa]Shear Modulus
[MPa]
Poisson’s RatioTensile Strength
[MPa]
Shear Strength
[MPa]
Compression Strength
[MPa]
E1
E2
90°
G1,2V12FTU1
FTU2
90°
FSU
45°
FCU1
FCU2
90°
143,530582638450.3622214983.5641114
Table 2. Parameters and levels.
Table 2. Parameters and levels.
ParameterLevels
123
Shapes of holesCircularSquareHexagonal
D/D0(S = const)21.671.5
S/D0(D0 = const)1.251.451.6
Table 3. Buckling load value for all tested cases for P1 composite configuration.
Table 3. Buckling load value for all tested cases for P1 composite configuration.
No.Hole ShapeOpening Ratio
D/D0
Spacing Ratio
S/D0
Pcr
[N]
Difference
[%]
1Without holes------2172.4-
2Circular21.671802.317.04
3Square1709.621.30
4Hexagonal1764.818.76
5Circular1.671.391667.323.25
6Square154129.06
7Hexagonal1612.925.75
8Circular1.51.251561.328.13
9Square1415.234.85
10Hexagonal1492.131.32
11Circular1.51.451551.328.59
12Square1410.835.06
13Hexagonal1487.331.54
14Circular1.51.61553.828.48
15Square1422.934,50
16Hexagonal1493.231.26
Table 4. Values of critical load for perforated and unperforated profiles, for four different composite lay-ups for chosen case.
Table 4. Values of critical load for perforated and unperforated profiles, for four different composite lay-ups for chosen case.
Laminate Lay-UpSymbolPcr [N]Difference [%]
Unperforated ProfilePerforated Profile
[0/45/−45/90/0]sP12172.41802.317.04
[45/−45/90/0]sP23259.32761.515.27
[90/−45/45/0]sP32183.31788.518.08
[90/0/90/0]sP41590.61124.229.32
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Falkowicz, K. Stability Analysis of Thin-Walled Perforated Composite Columns Using Finite Element Method. Materials 2022, 15, 8919. https://doi.org/10.3390/ma15248919

AMA Style

Falkowicz K. Stability Analysis of Thin-Walled Perforated Composite Columns Using Finite Element Method. Materials. 2022; 15(24):8919. https://doi.org/10.3390/ma15248919

Chicago/Turabian Style

Falkowicz, Katarzyna. 2022. "Stability Analysis of Thin-Walled Perforated Composite Columns Using Finite Element Method" Materials 15, no. 24: 8919. https://doi.org/10.3390/ma15248919

APA Style

Falkowicz, K. (2022). Stability Analysis of Thin-Walled Perforated Composite Columns Using Finite Element Method. Materials, 15(24), 8919. https://doi.org/10.3390/ma15248919

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