Next Article in Journal
Low Temperature Transient Liquid Phase Bonding of Alumina Ceramics with the Bi2O3-ZnO Interlayer
Next Article in Special Issue
Investigation of an Effective Anchoring Length of CFRP Tapes Used to Strengthen Steel Thin-Walled Beams with a Rectangular Cross-Section Subjected to Four-Point Bending
Previous Article in Journal
Adsorption Study of Continuous Heavy Metal Ions (Pb2+, Cd2+, Ni2+) Removal Using Cocoa (Theobroma cacao L.) Pod Husks
Previous Article in Special Issue
Numerical Investigations of Perforated CFRP Z-Cross-Section Profiles, under Axial Compression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensitivity of Compressed Composite Channel Columns to Eccentric Loading

Department of Machine Design and Mechatronics, Faculty of Mechanical Engineering, Lublin University of Technology, 36 Nadbystrzycka Street, 20-618 Lublin, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(19), 6938; https://doi.org/10.3390/ma15196938
Submission received: 5 August 2022 / Revised: 28 September 2022 / Accepted: 4 October 2022 / Published: 6 October 2022
(This article belongs to the Special Issue Computational and Experimental Mechanics of Thin-Walled Structures)

Abstract

:
This study investigated short thin-walled channel columns made of carbon/epoxy laminate. Columns with two multi-ply composite layups [0/45/−45/90]s and [90/−45/45/0]s were tested, with each layup having eight plies symmetric to the midplane. The columns were subjected to compressive loads, including an eccentric compressive load applied relative to the center of gravity of their cross-section. Simple support boundary conditions were applied to the ends of the columns. The scope of the study included analyzing the effect of load eccentricity on the buckling mode, bifurcation load (idealized structure), and critical load (structure with initial imperfections). The critical load for the actual structure was determined with the use of approximation methods, based on experimental postbuckling equilibrium paths. In parallel with the experiments, a numerical analysis was conducted using the finite element method and Abaqus® software (Dassault Systèmes, Vélizy-Villacoublay, France). The first stage of the numerical analysis consisted of solving an eigenproblem, in order to determine the mode of the loss of structural stability and to calculate the bifurcation loads for structures under axial and eccentric compression. The second stage of the numerical analysis involved examining the non-linear state of pre-deflected structures. Numerical postbuckling equilibrium paths were used to estimate the critical loads with an approximation method. The experimental results were used to validate the numerical models. This made it possible to determine the effect of compressive load eccentricity on the buckling mode and critical load of the tested structures. The results confirmed that compressive load eccentricity had a significant impact on the load-carrying capacity in the postbuckling range. This may potentially lead to premature damage to composite materials and, ultimately, to a reduced load-carrying capacity of structures.

1. Introduction

Bearing elements in modern thin-walled structures are predominantly made of advanced composite materials characterized by a low specific weight and high mechanical performance [1,2,3,4,5,6,7,8,9]. This primarily concerns thin-walled stiffeners, which are usually designed as open and closed profiles, with complex cross-sectional shapes. The main role of these structural members is to transfer axial and bending loads, which means that they are exposed to buckling, particularly under compressive loads. Although a loss of stability by these members may cause local weakening in a load-bearing structure itself, this poses no direct threat to the safe performance of such structures [10,11,12,13]. Studies have shown [13,14,15,16,17,18] that thin-walled composite structural members can continue carrying loads in the buckled state, provided that the buckling is elastic and that the postbuckling equilibrium path is stable [19,20,21,22,23,24,25,26]. Considering the above, in addition to investigating the buckling of structural members, it is also necessary to investigate the postbuckling behavior of entire structures and their load-carrying capacity in a buckled state. Stability loss in thin-walled stiffeners usually takes the form of local buckling of the walls and web of the thin-walled column, and this is manifested by the occurrence of a certain number of half-waves along the column axis [27,28,29,30,31,32,33,34]. The load that causes columns to buckle is known as the critical load, its value being dependent on many factors. These factors include the following: the column’s cross-sectional shape and dimensions, its length, boundary conditions, and—for composite materials—the ply layup in the column’s walls and web [35,36,37,38,39]. Knowledge of the critical load value is of vital importance when designing thin-walled structures, because it allows taking measures to prevent these structures from buckling within the service load range. The determination of critical loads for actual (not idealized) structures is complicated, which significantly impedes the design of such structures. The critical load is determined by approximation methods [39,40,41], which make it possible to determine an approximate value of the critical force value based on the postbuckling equilibrium path of a given structure; this path is usually plotted as the load versus displacement perpendicular to the column’s walls. The critical load also depends on the extent of geometrical imperfections, which usually result from the manufacturing inaccuracy of an actual structure. One other important aspect of the design of thin-walled structures concerns the real loads, which in most cases differ from the idealized theoretical loading conditions. One of the key problems in this respect is the load eccentricity resulting from, for example, assembly inaccuracies or unstable boundary conditions. The behavior of a structure under eccentric loading significantly differs from the idealized axial loading response. Load eccentricity may lead to a premature loss of stability, with buckling occurring within the service load range. Changes may also occur in the way the walls and web of a column are loaded, because the eccentricity of the force application induces additional flexural loads that may be a direct cause of premature buckling. If a thin-walled structure enters into the postbuckling state too early, this may result in a significant reduction in its strength parameters and thus lead to a faster failure of its load-carrying members. This is all the more dangerous for actual structure performance, because the vast majority of design solutions do not take the above fact into consideration at the design stage. Most structures are usually designed on the assumption that their structural members are under idealized loading conditions. This research undertook to determine the qualitative and quantitative effect of compressive load eccentricity on the buckling and postbuckling behavior of thin-walled composite columns. In the study, different load eccentricity values were tested using experimental and numerical methods. It is also worth mentioning that this study involved using the popular finite element method, which is widely used in many fields [42,43,44,45,46,47,48,49].

2. Object of the Study

The object of this study were short thin-walled columns with a channel cross-section under axial and eccentric compression. The columns were standard thin-walled structures with perpendicular walls; they were flat plate elements joined by their longer ends [33,38,39,50,51,52,53,54]. The columns were fabricated from epoxy/carbon laminate by autoclaving (Figure 1a). The composite channel columns had the following layups: C1 [0/45/−45/90/90/−45/45/0]T and C2 [90/−45/45/0/0/45/−45/90]T. The laminate consisted of eight plies arranged symmetrically relative to the midplane. The columns had the following overall dimensions: 60 × 30 × 250 × 0.84 mm—Figure 1b.
A schematic design of a channel column under an eccentric compressive load is shown in Figure 1b. Eccentricity was induced by moving the point of compressive force application relative to the longitudinal axis of the column. This point was moved from the center of gravity of the column’s cross-section to the 0° axis by a value of e = 10 mm, and then the column was gradually rotated by 0° < e = 10 mm < 90° every 15° relative to its axis (the point of compressive force application was changed in an arc with a radius of R = 10 mm, with its midpoint located in the center of gravity of the column—Figure 1b). The mechanical properties of a single composite ply, which were determined experimentally in compliance with the relevant ISO standard, were as follows: Young’s modulus along the fiber direction—E1 = 143.5 GPa, Young’s modulus perpendicular to the fiber direction—E2 = 5.83 GPa, Poisson’s ratio on the ply plane—ν12 = 0.36, and Kirchhoff’s modulus—G12 = 3.85 GPa.

3. Methodology of the Study

In this study, experimental and numerical methods were employed to investigate the buckling and early postbuckling responses of thin-walled composite structures under axial and eccentric compression. Experiments were performed up to approx. 200% of the numerical bifurcation load value. The objective of the study was to investigate the buckling behavior of these columns and to determine their postbuckling characteristics, describing the relationship between compressive load and deflection perpendicular to the column wall in the early postbuckling regime. This made it possible to assess the effect of compressive load eccentricity on the stability of the channel section columns. This paper is an extension of the team’s research, which was presented in a previous work [55].

3.1. Experimental

Experiments were conducted at ambient temperature on a static testing machine, Zwick-100 (Zwick Roell, Ulm, Germany), with the upper cross-bar travel speed maintained constant at 2 mm/min. During the experimental tests, the deformations of the column web were recorded using resistance strain gauges, taped on both sides of the web at the location of the expected largest deflections. In addition, the deflections of the column walls in the direction perpendicular to the wall plane were measured using a laser sensor. These measurements made it possible to determine the postbuckling characteristics of the columns and thus describe their postbuckling behavior for the tested loading conditions. The column was simply-supported at both ends, with ball-and-socket joint heads mounted on the pins of the testing machine—Figure 2b.
Specially designed sliding tables were mounted on the heads, in order to ensure that the eccentricity values were entered accurately via positioning the column relative to the static axis of the testing machine (Figure 2a). The centering of the channel column with respect to the pins of the testing machine was done by using plates with holes drilled at the center of gravity of the cross-section (Figure 3b). In order to ensure control of the load eccentricity realization with an accuracy of to 0.01 mm, an additional measuring system using electronic sensors was used (Figure 3a). The experiments were conducted for one axial compressive load case and two cases of eccentric loading: e10mm_0° and e10mm_90°. As a result, it was possible to validate the numerical FEM models for the extreme eccentricity values. Critical loads for the actual structure were calculated based on postbuckling equilibrium paths using approximation methods. Figure 3b shows the physical model of a channel column mounted in the heads of the testing machine, as well as the plates for positioning the column’s center of gravity relative to the axis of the testing machine.
During the tests, the course of the compression force, the deformation of the column web, and the deflection of the column walls in the perpendicular direction were registered using a laser sensor. The resulting equilibrium post-critical paths, which determine the dependence of the compressive load on the deflection of the P-w column walls, made it possible to determine the value of the critical load and to assess the performance of the structure in the in the post-critical range.

3.2. Numerical Analysis

In parallel with the experiments, numerical simulations were performed, to develop FEM models for analyzing the effect of the axial compressive load and eccentric load applied in a range of 0° < e = 10 mm < 90° every 15° on the tested structures. The columns were discretized, with the use of shell elements, and each element having three translational and three rotational degrees of freedom in every node. The numerical models of the columns were constructed with S8R finite elements, which are eight-node elements with a quadratic shape function and reduced integration [30,33,38,39]. The plates reproducing the simple support on the column ends were modelled using four-node rigid finite elements, R3D4.
To investigate the effect of meshing size, this study tested five different element sizes: from 1 × 1 mm to 5 × 5 mm, for the chosen channel column. The mesh convergence was established by increasing the mesh density in each part of the model. It was observed that there were no considerable changes in load response between a 1-mm and 5-mm element size (Figure 4), but the processing time was considerable and an increment in mesh density was unnecessary, as shown in Figure 4.
The designed numerical models had a uniform finite element mesh, with each finite element having a size of 4 × 4 mm. In effect, the individual walls of the column could be divided in a uniform manner by generating a constant-density mesh. This approach made it possible to examine the strains and stresses occurring in the compressed thin-walled structures.
The FEM model attempted to represent the actual support of the column at both ends of the testing machine (Figure 5a). Reference points corresponding to the position of the centers of gravity of the ball joints of the mounting heads were used, where simple support boundary conditions were defined. Two translational degrees of freedom (Ux = Uy = 0) and rotation relative to the axis of the column (URz = 0) were constrained at the point reproducing the upper ball-and-socket joint, whereas three translational degrees of freedom (Ux = Uy = Uz = 0) and rotation (URz = 0) were constrained at the reference point for the lower ball-and-socket joint. The reference points were linked with the rigid plates supporting the columns ends, by coupling all degrees of freedom of the points and the plates. The column was loaded with a compressive force concentrated at the upper reference point. A given value of the eccentricity parameter “e” was applied by changing the position of the column model relative to the rigid plates and displacing it perpendicular to the web (Figure 1b) by a value of 10 mm, and then rotating it with a fixed value every 15°. Contact relations were defined between the column ends and the plates. The structure of the laminate structure was modelled using the layup–ply technique, which made it possible to reproduce the layup of the composite material and ensure a uniform thickness of the finite elements by defining the ply thickness, material type, and fiber orientation (C1—Figure 5b, C2—Figure 5c).

3.3. Analysis of the Early Postbuckling Range

FEM numerical calculations were performed in two stages. The first involved solving an eigenproblem, to determine both the bifurcation load of an idealized structure using the minimum potential energy principle [56] and the corresponding buckling mode—Figure 4. For such a system, an eigenvalue problem was solved using the following equation [57]:
( [ K ] + λ i [ H ] ) { ψ } i = 0
where [K] is the structural stiffness matrix, [H] is the stress stiffness matrix, λi is the i-th eigenvalue, and ψ is the i-th eigenvector of displacement.
Equation (1) is satisfied if the eigenvector of displacement is equal to zero, or if the determinant of the term in brackets is equal to zero. When {ψ}i = 0, this is a trivial solution, and it is out of interest. This means that the structure remains in the initial state of equilibrium. The term in brackets in Equation (1) gives the following solution:
| [ K ] + λ i [ H ] | = 0
Equation (2) represents the eigenvalue problem that allows finding the n values of the buckling load multiplier λ and the corresponding buckling mode shape.
The other stage of the numerical calculations was a nonlinear analysis of the early postbuckling response up to 200% of the bifurcation load. A nonlinear analysis was performed, with the use of a model of the structure with initial geometric imperfections reflecting the lowest buckling mode; the initial amplitude of these imperfections being w0 = 0.1 mm—Figure 6. The value of the initial imperfections was determined by experimental validation of the numerical postbuckling characteristics of the tested structures. The nonlinear stability problem was solved using the Newton–Raphson method. The other stage of the calculations involved determining the postbuckling equilibrium paths, which were used as a basis for determining the critical loads of the structures with initial geometric imperfections.

3.4. Method of Determining the Critical Load of a Structure with Initial Imperfections

Various types of inaccuracies could be observed in the experiments, resulting from different factors, such as geometrical imperfections, the load application, and boundary conditions. In effect, it was difficult to determine the exact critical load values for the thin-walled structures. In cases such as these, approximation methods must be employed to determine critical loads based on the obtained experimental and numerical postbuckling equilibrium paths. In this study, critical loads were determined using an approximation method developed by Koiter.
The Koiter method was used to approximate the experimental and numerical postbuckling equilibrium paths describing the relationship between the load and deflection perpendicular to the column wall. The postbuckling equilibrium paths P-w for the early postbuckled state were approximated using a quadratic function, having the form [40]:
P = P c r a 2 a 0 w 2 + P c r a 1 a 0 + P c r
where a0, a1, and a2 are the parameters of the function, P is the load, Pcr is the value of an unknown critical load, and w is the increase in deflection measured perpendicular to the column’s wall.
According to Koiter’s method, the critical force is defined as a point of intersection between the quadratic function and the vertical axis of the coordinate system in the P-w plot (Figure 7). The accuracy of the critical load value depends on the employed approximation range, which means that for a stable postbuckling equilibrium path, the leading coefficient of a quadratic polynomial must be positive. Approximation accuracy is estimated with a correlation coefficient R2. The correlation coefficient value shows whether the approximation function agrees with a selected range of the experimental approximated curve. The higher the correlation coefficient value, the more accurate the approximation process becomes. With the approximation method employed in this study, the minimal value of the correlation coefficient was equal to R2 ≥ 95%.

4. Results and Discussion

The experiments and numerical analysis of thin-walled composite channel columns under axial and eccentric compression provided data about the effect of compressive load eccentricity on the buckling and postbuckling behavior of the tested columns. The results allowed for an examination of the buckling state based on the buckling modes of the actual and numerical structures. Figure 8 shows the experimental and numerical buckling modes obtained for C1 [0/45/−45/90/90/−45/45/0]T, while Figure 9 shows the experimental and numerical buckling modes for C2 [90/−45/45/0/0/45/−45/90]T.
For the axially compressed column C1 (Figure 8a,b), local buckling of the web and walls occurs, having the form of two longitudinal half-waves symmetric to the longitudinal symmetry plane of the column. There is no change in the buckling mode of the C1 column when the eccentricity value of e10mm_0° is applied; however, one can observe increased loading of the column’s walls. This results from the fact that the point of compressive force application was moved, which led to unloading of the web of the column (Figure 8c,d). The eccentricity value of 15° < e = 10 mm < 75° causes a total change in the buckling mode, from two to three half-waves (Figure 8e–h). For the extreme C1 eccentricity case of e10mm_90°, local buckling returns to the initial mode, i.e., two half-waves.
The axial compression of C2 (Figure 9a,b) is characterized by local buckling of the web and walls, with the buckling having the form of five half-waves symmetric to the symmetry plane of the column. For the eccentricity applied from the web of the C2 column, i.e., e10mm_0° (Figure 9c,d), the number of half-waves on the column walls and web remains the same, but—as in the case of C1—unloading of the web and increased loading of the walls occurs at the same time. For the eccentricity value of 15° < e = 10 mm < 90° there is a change in the buckling mode; the half-waves on the unloaded wall and the web disappear, while the half-waves in proximity of the bottom end of the column become more visible (Figure 9e–k). The experimental and numerical buckling modes showed a qualitative agreement.
To determine the quantitative relationship between load eccentricity and buckling, the obtained bifurcation loads describing the lowest buckling modes of the compressed composite columns were analyzed numerically. Table 1 shows the FEM eigenvalue buckling loads obtained for these structures. The difference between the eigenvalue buckling loads of C1 and C2 for the axial compression case was insignificant and only amounts to 0.5%, whereas the greatest difference of 4.6% can be observed for the eccentricity value of e10mm_30°. That means that both structures, which differ with respect to the position of the 0° and 90° plies, have a similar stiffness and thus undergo buckling under similar critical loads.
The numerical bifurcation loads obtained for the tested eccentric compression loads are shown graphically in Figure 10 and Figure 11. The C1 column exhibits a clear decrease in the bifurcation load under eccentric loading—Figure 10. The highest drop in the bifurcation load can be observed for C1 under e10mm_0° and amounts to 55% compared to the axial compression case, while the minimum difference is the 21% observed for C1 under e10mm_90°. The C2 column [90/−45/45/0]s is slightly more sensitive to load eccentricity, ranging <33% ÷ 57% >—Figure 11. The highest decrease in the bifurcation load of C2 can be observed for the eccentricity value of e10mm_15°. The results demonstrate that C1 [0/45/−45/90]s and C2 [90/−45/45/0]s had similar bifurcation loads under the applied eccentric load conditions.
The critical loads of the pre-deflected structures were calculated based on the experimental and numerical postbuckling equilibrium paths describing a load-deflection relationship measured perpendicularly to the channel column wall. A selection of the experimental approximated characteristics are shown in Figure 12, while the numerical ones are shown in Figure 13.
The postbuckling equilibrium paths describing the relationship between a compressive load P and a deflection w served as the basis for determining critical loads using Koiter’s approximation method [38,39,40]. For all tested cases the approximated experimental curve range included some portion of the early postbuckling path, from a clear inflection point of the load-deflection line, until the end of the experimental curve. For all analyzed cases, the coefficient of correlation was maintained at a very high level of R2 ˃ 95%. The numerical and experimental critical loads obtained with this method are given in Table 2.
The approximated values of the experimental and numerical critical loads for the structures with initial imperfections show a good agreement. The greatest difference between the critical load values of C1 [0/45/−45/90]s can be observed for the eccentricity of e10mm_90°, with this difference amounting to about 1%, whereas for C2 [90/−45/45/0]s the difference is slightly higher, yet does not exceed 3% (C2 e10mm_0°). It should be emphasized that the experimental critical loads are the lower estimation of the FEM critical load values, which results from the structural idealization of the numerical model.
To estimate the effect of load eccentricity on the critical load of the analyzed columns (C1 and C2), diagrams illustrating the relationship between the critical load and eccentricity were created (Figure 14 and Figure 15). The diagrams show a very good agreement between the numerical results and the experimental findings, which confirms that the numerical models were correct.
The results of the postbuckling behavior of the compressed channel columns demonstrate that the deflection reflecting the lowest buckling modes increases when increasing the compressive load. The tests were conducted up to 200% of the bifurcation load. Figure 16 and Figure 17 show the buckling modes for the two compressed channel columns: C1 [0/45/−45/90]s and C2 [90/−45/45/0]s.
The postbuckling behavior of the tested columns was analyzed by comparing the experimental and numerical postbuckling equilibrium paths obtained for these structures. The experimental equilibrium paths describing the relationship between the load and deflection were created based on the load results and measurements made with a laser beam pointed perpendicularly at a half-wave on the channel wall with the largest deflection. The numerical equilibrium paths were generated by measuring the column wall’s deflection perpendicularly, with the measurements made at the same point as in the experiments. Figure 18 compares the experimental and numerical paths obtained for C1 and C2.
The experimental and numerical results show a good agreement. Similarly to buckling, the experimental postbuckling equilibrium paths show a slight reduction in the column’s stiffness compared to the numerical paths. Nevertheless, the experimental and numerical paths indicate a stable behavior of the structures, which confirms that the structures can carry loads in the postbuckling range. Figure 19 shows the numerical postbuckling equilibrium paths obtained for C1 and C2, and considering the eccentricity.
The results demonstrate that the postbuckling stiffness of C1 and C2 decreased when eccentricity was applied. Regarding C1 (Figure 19a), the maximum decrease in the stiffness of this structure under eccentric load was observed for e10mm_30° and amounted to 57% (compared to the axial compression case). As for C2 (Figure 19b), the eccentricity of e10mm_0° caused the highest decrease in the column’s stiffness, of 52%. An analysis of the other cases tested demonstrated a 10% increase in the stiffness of C2 for the eccentricity value of e10mm_90°. The smallest decrease in the stiffness of eccentrically loaded C1 and C2 was observed for the eccentricity value of e10mm_90°.
The high costs of producing carbon–epoxy composite profiles make tests of this type expensive. This is a major limitation for scientists who research this type of structure. It should be added that FEM modelling may be interesting for considering the progressive failure and failure processes in composites [58,59,60,61,62,63,64,65]. For this reason, future research directions will focus on the analysis of the post-critical state in the full range of loads leading to the damage and failure of the composite structure.

5. Conclusions

This study investigated the effect of compressive load eccentricity on the stability and postbuckling behavior of thin-walled composite channel columns. A combined qualitative and quantitative analysis of the effect of the compressive load eccentricity on critical loads and postbuckling equilibrium paths of these structures was performed. The results demonstrated that the eccentricity applied perpendicular and parallel to the web of the column did not change its buckling mode, when compared to the axial compression case. The buckling mode would only change if the eccentricity was applied simultaneously in two directions (in an angle range of 15° ÷ 75°). The weakening of the column cross-section and the change in the form of buckling was due to an increase in the additional bending state, resulting from the eccentricity of the load.
The results showed that the tested ply layups had an insignificant effect on the buckling behavior of the columns. Fiber orientation had a negligible effect on the critical load. Despite changing the fiber orientations of the 0° and 90° plies, no significant differences were observed between the critical loads obtained for C1 and C2, irrespective of the eccentricity value. Nevertheless, a clear relationship was observed between the eccentricity and critical load. Compared to the axial compression case, the C1 column showed a significant decrease of 55% in the bifurcation load for the eccentricity value of e10mm_0°. The C2 column was even more sensitive to load eccentricity; the bifurcation load decrease was as high as 57% for the eccentricity value of e10mm_15°. Therefore, for the tested cases, load eccentricity may have caused premature buckling of the structures and thus considerably reduced their postbuckling load-carrying capacity. It should also be highlighted that the approach adopted in this study for determining the critical load based on postbuckling equlibrium paths was correct, as proven by the very high agreement between the experimental and numerical critical loads for two extreme eccentricity cases (axial and eccentric compression, e10mm_0° and e10mm_90°). This fact also confirms that the numerical models of the tested structures were correct.
The results also showed the good agreement of the experimental and numerical postbuckling equilibrium paths describing the relationship between the load and deflection measured perpendicular to the plane of the column wall. The results demonstrated that load eccentricity had a significant impact on structural stiffness in the postbuckling range. Compared to the axial compression case, the application of the eccentricity value of e10mm_30° caused a 57% reduction in the stiffness of C1, whereas the eccentricity of e10mm_0° resulted in a 52% decrease in the stiffness of C2. These results confirmed that the compressive load eccentricity had a significant impact on load-carrying capacity in the postbuckling range. This may potentially lead to premature damage of composite materials and, ultimately, to a reduced load-carrying capacity of structures. The results are of practical significance, particularly for the design of thin-walled composite structures that are exposed to unexpected inaccuracies, which may considerably affect their performance.

Author Contributions

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

Funding

The project/research was financed within the framework of the project 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

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carlsson, L.A.; Adams, D.F.; Pipes, R.B. Experimental Characterization of Advanced Composite Materials, 4th ed.; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  2. Wu, M.Q.; Zhang, W.; Niu, Y. Experimental and numerical studies on nonlinear vibrations and dynamic snap-through phenomena of bistable asymmetric composite laminated shallow shell under center foundation excitation. Eur. J. Mech. A Solid. 2021, 89, 104303. [Google Scholar] [CrossRef]
  3. Yang, S.W.; Hao, Y.X.; Zhang, W.; Yang, L.; Liu, L.T. Free vibration and buckling of eccentric rotating FG-GPLRC cylindrical shell using first-order shear deformation theory. Compos. Struct. 2021, 263, 113728. [Google Scholar] [CrossRef]
  4. Siriguleng, B.; Zhang, W.; Liu, T.; Liu, Y.Z. Vibration modal experiments and modal interactions of a large space deployable antenna with carbon fiber material and ring-truss structure. Eng. Struct. 2020, 207, 109932. [Google Scholar] [CrossRef]
  5. Yang, S.W.; Zhang, W.; Hao, Y.X.; Niu, Y. Nonlinear vibrations of FGM truncated conical shell under aerodynamics and in-plane force along meridian near internal resonances. Thin-Walled Struct. 2019, 142, 369–391. [Google Scholar] [CrossRef]
  6. Liu, T.; Zhang, W.; Wang, J.F. Nonlinear dynamics of composite laminated circular cylindrical shell clamped along a generatrix and with membranes at both ends. Nonlinear Dyn. 2017, 90, 1393–1417. [Google Scholar] [CrossRef]
  7. Chen, J.E.; Zhang, W.; Guo, X.Y.; Sun, M. Theoretical and experimental studies on nonlinear oscillations of symmetric cross-ply composite laminated plates. Nonlinear Dyn. 2013, 73, 1697–1714. [Google Scholar] [CrossRef]
  8. Bui, T.Q.; Hu, X. A review of phase-field models, fundamentals and their applications to composite laminates. Eng. Fract. Mech. 2021, 248, 107705. [Google Scholar] [CrossRef]
  9. Zhang, P.; Hu, X.; Bui, T.Q.; Yao, W. Phase field modeling of fracture in fiber reinforced composite laminate. Int. J. Mech. Sci. 2019, 161–162, 105008. [Google Scholar] [CrossRef]
  10. Bazant, Z.P.; Cedolin, L. Stability of Structures. In Elastic, Inelastic, Fracture and Damage Theories; Oxford University Press: Oxford, UK, 2010. [Google Scholar]
  11. Singer, J.; Arbocz, J.; Weller, T. Buckling Experiments. In Experimental Methods in Buckling of Thin-Walled Structure: Basic Concepts, Columns, Beams, and Plates; John Wiley & Sons Inc.: New York, NY, USA, 1998; Volume 2, p. 1. [Google Scholar]
  12. Kopecki, T.; Bakunowicz, J.; Lis, T. Post-critical deformation states of composite thin-walled aircraft load-bearing structures. J. Theor. Appl. Mech. 2016, 54, 195–204. [Google Scholar] [CrossRef]
  13. 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. Eksploat. Niezawodn. Maint. Reliab. 2016, 18, 16–24. [Google Scholar] [CrossRef]
  14. Fascetti, A.; Feo, L.; Nistico, N.; Penna, R. Web-flange behavior of pultruded GFRP Ibeams: A lattice model for the interpretation of experimental results. Compos. B Eng. 2016, 100, 257–269. [Google Scholar] [CrossRef]
  15. Feo, L.; Latour, M.; Penna, R.; Rizzano, G. Pilot study on the experimental behavior of GFRP-steel slip-critical connections. Compos. B Eng. 2017, 115, 209–222. [Google Scholar] [CrossRef]
  16. Berardi, V.P.; Perrella, M.; Feo, L.; Cricrì, G. Creep behavior of GFRP laminates and their phases: Experimental investigation and analytical modeling. Compos. B Eng. 2017, 122, 136–144. [Google Scholar] [CrossRef]
  17. Rozylo, P. Failure analysis of thin-walled composite structures using independent advanced damage models. Compos. Struct. 2021, 262, 113598. [Google Scholar] [CrossRef]
  18. Kubiak, T.; Kolakowski, Z.; Swiniarski, J.; Urbaniak, M.; Gliszczynski, A. Local buckling and post-buckling of composite channel-section beams—Numerical and experimental investigations. Compos. B Eng. 2016, 91, 176–188. [Google Scholar] [CrossRef]
  19. Wysmulski, P.; Dȩbski, H. The analysis of sensitivity to eccentric load of compressed thin-walled laminate columns. AIP Conf. Proc. 2019, 2078, 020006. [Google Scholar] [CrossRef]
  20. Debski, H.; Teter, A.; Kubiak, T.; Samborski, S. Local buckling, post-buckling and collapse of thin-walled channel section composite columns subjected to quasi-static compression. Compos. Struct. 2016, 136, 593–601. [Google Scholar] [CrossRef]
  21. Madukauwa-David, I.D.; Drissi-Habti, M. Numerical simulation of the mechanical behavior of a large smart composite platform under static loads. Compos. B Eng. 2016, 88, 19–25. [Google Scholar] [CrossRef]
  22. 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]
  23. Banat, D.; Mania, R.J. Failure assessment of thin-walled FML profiles during buckling and postbuckling response. Compos. B Eng. 2017, 112, 278–289. [Google Scholar] [CrossRef]
  24. Chroscielewski, J.; Miskiewicz, M.; Pyrzowski, L.; Sobczyk, B.; Wilde, K. A novel sandwich footbridge—Practical application of laminated composites in bridge design and in situ measurements of static response. Compos. B Eng. 2017, 126, 153–161. [Google Scholar] [CrossRef]
  25. Wysmulski, P. The analysis of buckling and post buckling in the compressed composite columns. Arch. Mater. Sci. Eng. 2017, 85, 35–41. [Google Scholar] [CrossRef]
  26. 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]
  27. Rozylo, P. Stability and failure of compressed thin-walled composite columns using experimental tests and advanced numerical damage models. Int. J. Numer. Methods Eng. 2021, 122, 5076–5099. [Google Scholar] [CrossRef]
  28. Debski, H.; Kubiak, T.; Teter, A. Buckling and postbuckling behaviour of thin-walled composite channel section column. Compos. Struct. 2013, 100, 195–204. [Google Scholar] [CrossRef]
  29. Debski, H.; Kubiak, T.; Teter, A. Experimental investigation of channel-section composite profiles behaviour with various sequences of plies subjected to static compression. Thin-Walled Struct. 2013, 71, 147–154. [Google Scholar] [CrossRef]
  30. Wysmulski, P.; Debski, H. The Effect of Eccentricity of Load on the Behavior of Compressed Composite Columns in Critical State. Polym. Compos. 2019, 40, 70–77. [Google Scholar] [CrossRef] [Green Version]
  31. Wysmulski, P.; Debski, H. Stability Analysis of Composite Columns under Eccentric Load. Appl. Compos. Mater. 2019, 26, 683–692. [Google Scholar] [CrossRef] [Green Version]
  32. Debski, H.; Rozylo, P.; Wysmulski, P. Stability and load-carrying capacity of short open-section composite columns under eccentric compression loading. Compos. Struct. 2020, 252, 112716. [Google Scholar] [CrossRef]
  33. Wysmulski, P.; Debski, H.; Falkowicz, K. Stability analysis of laminate profiles under eccentric load. Compos. Struct. 2020, 238, 111944. [Google Scholar] [CrossRef]
  34. Banat, D.; Kołakowski, Z.; Mania, R. Investigations of FML profile buckling and post-buckling behaviour under axial compression. Thin-Walled Struct. 2016, 107, 335–344. [Google Scholar] [CrossRef]
  35. Ascione, F. Influence of initial geometric imperfections in the lateral buckling problem of thin walled pultruded GFRP I-profiles. Compos. Struct. 2014, 112, 85–99. [Google Scholar] [CrossRef]
  36. Li, Z.M.; Qiao, P. Buckling and postbuckling behavior of shear deformable anisotropic laminated beams with initial geometric imperfections subjected to axial compression. Eng. Struct. 2015, 85, 277–292. [Google Scholar] [CrossRef]
  37. Urbaniak, M.; Teter, A.; Kubiak, T. Influence of boundary conditions on the critical and failure load in the GFPR channel cross-section columns subjected to compression. Compos. Struct. 2015, 134, 199–208. [Google Scholar] [CrossRef]
  38. 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. B Eng. 2021, 226, 109346. [Google Scholar] [CrossRef]
  39. 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]
  40. Paszkiewicz, M.; Kubiak, T. Selected problems concerning determination of the buckling load of channel section beams and columns. Thin-Walled Struct. 2015, 93, 112–121. [Google Scholar] [CrossRef]
  41. Wysmulski, P.; Falkowicz, K.; Filipek, P. Buckling state analysis of compressed composite plates with cut-out. Compos. Struct. 2021, 274, 114345. [Google Scholar] [CrossRef]
  42. Jonak, J.; Karpinski, R.; Wojcik, A.; Siegmund, M.; Kalita, M. Determining the Effect of Rock Strength Parameters on the Breakout Area Utilizing the New Design of the Undercut/Breakout Anchor. Materials 2022, 15, 851. [Google Scholar] [CrossRef] [PubMed]
  43. Gliszczynski, A.; Samborski, S.; Wiacek, N.; Rzeczkowski, J. Mode I interlaminar fracture of glass/epoxy unidirectional laminates. Part II: Numerical analysis. Materials 2019, 12, 1604. [Google Scholar] [CrossRef]
  44. Jonak, J.; Karpinski, R.; Siegmund, M.; Machrowska, A.; Prostanski, D. Experimental Verification of Standard Recommendations for Estimating the Load-Carrying Capacity of Undercut Anchors in Rock Material. Adv. Sci. Technol. Res. J. 2021, 15, 230–244. [Google Scholar] [CrossRef]
  45. Grzejda, R. FE-modelling of a contact layer between elements joined in preloaded bolted connections for the operational condition. Adv. Sci. Technol. Res. J. 2014, 8, 19–23. [Google Scholar] [CrossRef]
  46. Jonak, J.; Karpinski, R.; Wojcik, 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]
  47. Grzejda, R. Finite element modeling of the contact of elements preloaded with a bolt and externally loaded with any force. J. Comput. Appl. Math. 2021, 393, 113534. [Google Scholar] [CrossRef]
  48. Jonak, J.; Karpinski, R.; Wojcik, A. Influence of the Undercut Anchor Head Angle on the Propagation of the Failure Zone of the Rock Medium—Part II. Materials 2021, 14, 3880. [Google Scholar] [CrossRef]
  49. Jonak, J.; Karpinski, R.; Wojcik, 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] [PubMed]
  50. 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]
  51. Falkowicz, K.; Debski, H. The post-critical behaviour of compressed plate with non-standard play orientation. Compos. Struct. 2020, 252, 112701. [Google Scholar] [CrossRef]
  52. Falkowicz, K.; Dębski, H.; Wysmulski, P.; Różyło, P. The behaviour of compressed plate with a central cut-out, made of composite in an asymmetrical arrangement of layers. Compos. Struct. 2019, 214, 406–413. [Google Scholar] [CrossRef]
  53. 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]
  54. 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]
  55. 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]
  56. Abaqus HTML Documentation; Abaqus: Johnston, RI, USA, 2016.
  57. Kubiak, T. Static and Dynamic Buckling of Thin-Walled Plate Structures; Springer: Cham, Switzerland; Heidelberg, Germany; New York, NY, USA; Dordrecht, The Netherlands; London, UK, 2013. [Google Scholar] [CrossRef]
  58. Szewczak, I.; Rozylo, P.; Snela, M.; Rzeszut, K. Impact of Adhesive Layer Thickness on the Behavior of Reinforcing Thin-Walled Sigma-Type Steel Beams with CFRP Tapes. Materials 2022, 15, 1250. [Google Scholar] [CrossRef]
  59. Sohn, M.S.; Hu, X.Z.; Kim, J.K.; Walker, L. Impact damage characterization of carbon fibre/epoxy composites with multi-layer reinforcement. Compos. B Eng. 2000, 31, 681–691. [Google Scholar] [CrossRef]
  60. Rozylo, P. Comparison of failure for thin-walled composite columns. Materials 2022, 15, 167. [Google Scholar] [CrossRef]
  61. Batra, R.C.; Gopinath, G.; Zheng, J.Q. Damage and failure in low energy impact of fiber-reinforced polymeric composite laminates. Compos. Struct. 2012, 94, 540–547. [Google Scholar] [CrossRef]
  62. Rozylo, P.; Debski, H. Failure study of compressed thin-walled composite columns with top-hat cross-section. Thin-Walled Struct. 2022, 180, 109869. [Google Scholar] [CrossRef]
  63. Hu, H.; Niu, F.; Dou, T.; Zhang, H. Rehabilitation effect evaluation of CFRP-lined prestressed concrete cylinder pipe under combined loads using numerical simulation. Math. Probl. Eng. 2018, 2018, 3268962. [Google Scholar] [CrossRef]
  64. Jonak, J.; Karpinski, R.; Wojcik, A. Numerical analysis of undercut anchor effect on rock. J. Phys. Conf. Ser. 2021, 2130, 012011. [Google Scholar] [CrossRef]
  65. Jonak, J.; Karpinski, R.; Wojcik, A. Numerical analysis of the effect of embedment depth on the geometry of the cone failure. J. Phys. Conf. Ser. 2021, 2130, 012012. [Google Scholar] [CrossRef]
Figure 1. Channel column: (a) actual specimen, (b) geometrical model of a column under eccentric compression.
Figure 1. Channel column: (a) actual specimen, (b) geometrical model of a column under eccentric compression.
Materials 15 06938 g001
Figure 2. Experimental setup: (a) jointed head, with attached table for precision eccentricity of load, (b) cross-section model of the ball-and-socket joint head.
Figure 2. Experimental setup: (a) jointed head, with attached table for precision eccentricity of load, (b) cross-section model of the ball-and-socket joint head.
Materials 15 06938 g002
Figure 3. Experimental setup: (a) calibration of ball-and-socket joint head position, (b) position of a specimen on the sliding tables.
Figure 3. Experimental setup: (a) calibration of ball-and-socket joint head position, (b) position of a specimen on the sliding tables.
Materials 15 06938 g003
Figure 4. The convergence test results.
Figure 4. The convergence test results.
Materials 15 06938 g004
Figure 5. FEM analysis: (a) discrete model of a column and the boundary conditions for the eccentricity value of e10mm_45, (b) layup C1, (c) layup C2.
Figure 5. FEM analysis: (a) discrete model of a column and the boundary conditions for the eccentricity value of e10mm_45, (b) layup C1, (c) layup C2.
Materials 15 06938 g005
Figure 6. Postbuckling equilibrium paths of two structures: idealized, and actual (with initial geometric imperfections described by an amplitude w0).
Figure 6. Postbuckling equilibrium paths of two structures: idealized, and actual (with initial geometric imperfections described by an amplitude w0).
Materials 15 06938 g006
Figure 7. Equilibrium path approximation using Koiter’s method.
Figure 7. Equilibrium path approximation using Koiter’s method.
Materials 15 06938 g007
Figure 8. Buckling modes for C1 [0/45/−45/90]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Figure 8. Buckling modes for C1 [0/45/−45/90]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Materials 15 06938 g008aMaterials 15 06938 g008b
Figure 9. Buckling modes for C2 [90/−45/45/0]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Figure 9. Buckling modes for C2 [90/−45/45/0]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Materials 15 06938 g009aMaterials 15 06938 g009b
Figure 10. Eigenvalue buckling loads for C1 [0/45/−45/90]s.
Figure 10. Eigenvalue buckling loads for C1 [0/45/−45/90]s.
Materials 15 06938 g010
Figure 11. Eigenvalue buckling loads for C2 [90/−45/45/0]s.
Figure 11. Eigenvalue buckling loads for C2 [90/−45/45/0]s.
Materials 15 06938 g011
Figure 12. Experimental critical loads for structures with initial imperfections: (a) C1 axis, (b) C1 e10_0°, (c) C2 axis, (d) C2 e10_0°.
Figure 12. Experimental critical loads for structures with initial imperfections: (a) C1 axis, (b) C1 e10_0°, (c) C2 axis, (d) C2 e10_0°.
Materials 15 06938 g012
Figure 13. Numerical critical loads for structures with initial imperfections: (a) C1 e10_45°, (b) C1 e10_75°, (c) C2 e10_0°, (d) C2 e10_15°.
Figure 13. Numerical critical loads for structures with initial imperfections: (a) C1 e10_45°, (b) C1 e10_75°, (c) C2 e10_0°, (d) C2 e10_15°.
Materials 15 06938 g013
Figure 14. Eccentricity versus critical load for C1 [0/45/−45/90]s.
Figure 14. Eccentricity versus critical load for C1 [0/45/−45/90]s.
Materials 15 06938 g014
Figure 15. Eccentricity versus critical load for C2 [90/−45/45/0]s.
Figure 15. Eccentricity versus critical load for C2 [90/−45/45/0]s.
Materials 15 06938 g015
Figure 16. Early postbuckling modes for C1 [0/45/−45/90]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Figure 16. Early postbuckling modes for C1 [0/45/−45/90]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Materials 15 06938 g016
Figure 17. Early postbuckling modes for C2 [90/−45/45/0]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Figure 17. Early postbuckling modes for C2 [90/−45/45/0]s: (a) axis EXP, (b) axis FEM, (c) e10_0° EXP, (d) e10_0° FEM, (e) e10_15° FEM, (f) e10_30° FEM, (g) e10_45° FEM, (h) e10_60° FEM, (i) e10_75° FEM, (j) e10_90° EXP, (k) e10_90° FEM.
Materials 15 06938 g017
Figure 18. Experimental and numerical postbuckling equilibrium paths for: (a) C1 e10mm_0°, (b) C1 e10mm_90°, (c) C2 e10mm_0°, (d) C2 e10mm_90°.
Figure 18. Experimental and numerical postbuckling equilibrium paths for: (a) C1 e10mm_0°, (b) C1 e10mm_90°, (c) C2 e10mm_0°, (d) C2 e10mm_90°.
Materials 15 06938 g018
Figure 19. Load eccentricity versus postbuckling behavior for: (a) C1, (b) C2.
Figure 19. Load eccentricity versus postbuckling behavior for: (a) C1, (b) C2.
Materials 15 06938 g019
Table 1. Eigenvalue buckling loads of the tested structures.
Table 1. Eigenvalue buckling loads of the tested structures.
Eigenvalue Buckling Loads [N]
Eccentricity
Axise10_0°e10_15°e10_30°e10_45°e10_60°e10_75°e10_90°
C1 [0/45/−45/90]s217098799310111072119013901714
C2 [90/−45/45/0]s21829629489641023113513261641
Diference−0.55%2.53%4.53%4.65%4.57%4.62%4.60%4.26%
Table 2. Critical loads for structures with initial imperfections.
Table 2. Critical loads for structures with initial imperfections.
Critical Loads (Koiter Method) [N]
Eccentricity
Axise10_0°e10_15°e10_30°e10_45°e10_60°e10_75°e10_90°
C1_FEM21419549649861057119813271665
C1_EXP2123951-----1649
Diference0.84%0.31%-----0.96%
C2_FEM2145912871934979109513011611
C2_EXP2126887-----1576
Diference0.89%2.74%-----2.17%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wysmulski, P.; Debski, H.; Falkowicz, K. Sensitivity of Compressed Composite Channel Columns to Eccentric Loading. Materials 2022, 15, 6938. https://doi.org/10.3390/ma15196938

AMA Style

Wysmulski P, Debski H, Falkowicz K. Sensitivity of Compressed Composite Channel Columns to Eccentric Loading. Materials. 2022; 15(19):6938. https://doi.org/10.3390/ma15196938

Chicago/Turabian Style

Wysmulski, Pawel, Hubert Debski, and Katarzyna Falkowicz. 2022. "Sensitivity of Compressed Composite Channel Columns to Eccentric Loading" Materials 15, no. 19: 6938. https://doi.org/10.3390/ma15196938

APA Style

Wysmulski, P., Debski, H., & Falkowicz, K. (2022). Sensitivity of Compressed Composite Channel Columns to Eccentric Loading. Materials, 15(19), 6938. https://doi.org/10.3390/ma15196938

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