Authors:
Ikram Feddal Department of Industrial and Civil Sciences and Technologies, National School of Applied Sciences, University Abdelmalek Essaadi, 93030, Tetouan, Morocco

Search for other papers by Ikram Feddal in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0003-2186-8119
and
Abdelatif Khamlichi Department of Industrial and Civil Sciences and Technologies, National School of Applied Sciences, University Abdelmalek Essaadi, 93030, Tetouan, Morocco

Search for other papers by Abdelatif Khamlichi in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Laminated composite shell panels take part in several engineering structures. Due to their complex nature, failure modes in composites are highly dependent on the geometry, direction of loading and orientation of the fibers. However, the design of composite parts is still a delicate task because of these fiber failure modes, which includes matrix failure modes or other so-called interlaminar interface failure such as delamination, that corresponds to the separation of adjacent layers of the laminate as a consequence of the weakening of interface layer between them. In this work, impact-induced delamination represented as a circular single delamination is investigated, as it can reduce greatly the structural integrity without getting detected. Furthermore, attention is focused on its effect upon the post-buckling response and the compressive strength of a composite panel. The delamination buckling was modelled using the cohesive element technique under Abaqus software, in order to predict delamination growth and damage propagation while observing their effects on the critical buckling load.

Abstract

Laminated composite shell panels take part in several engineering structures. Due to their complex nature, failure modes in composites are highly dependent on the geometry, direction of loading and orientation of the fibers. However, the design of composite parts is still a delicate task because of these fiber failure modes, which includes matrix failure modes or other so-called interlaminar interface failure such as delamination, that corresponds to the separation of adjacent layers of the laminate as a consequence of the weakening of interface layer between them. In this work, impact-induced delamination represented as a circular single delamination is investigated, as it can reduce greatly the structural integrity without getting detected. Furthermore, attention is focused on its effect upon the post-buckling response and the compressive strength of a composite panel. The delamination buckling was modelled using the cohesive element technique under Abaqus software, in order to predict delamination growth and damage propagation while observing their effects on the critical buckling load.

1 Introduction

Aerospace industry are increasingly moving towards the massive use of composite materials in aircraft structures such as in Airbus 350, Airbus 380 and Boeing 787 [1]. Composites have attractive mechanical properties like the increased strength and high specific stiffness combined with weight reduction, compared to conventional materials. Hence, the extensive use of composite materials instead of aluminum alloys, contributes to reduce operating costs significantly. This makes it possible to meet the prediction announced in the aeronautics industry to reduce costs by 20% over the short period and to consider a 50% reduction in the long term [2]. These advantages should not cover fragility of composites to undergo damage even under low-speed impacts. This is due to the low thicknesses used and the risk of inter-ply delamination which corresponds to the decohesion of two layers of the laminate. This damage can be caused by falling tools during manufacturing and maintenance operations, bird strike, and can seriously degrade the laminate’s compressive strength and buckling stability [3, 4].

Due to their complex nature, failure modes in composites are highly dependent on the geometry, direction of loading and orientation of the plies [5]. There is damage at the level of the fibers, at the level of the matrix or other so-called inter-laminar as delamination, which is a particular concern. It is one of the predominant forms of damage caused by manufacturing defects and high stress concentrations due to geometric discontinuity [6]. In buckling conditions, delamination can expand and further reduce the strength of composite structures. It then poses a serious threat to the safety of the structure and the consequences can be catastrophic.

Delamination is a failure in a composite material, which leads to the separation of plies. It is due to various causes such as [7]:

  1. - Geometric discontinuities: Stresses that occur between stiffeners and thin plates, free edges, joints, and holes promote the initiation of delamination and trigger mechanisms of delamination between plies. In a laminate composed of several plies with various orientations, the plies mutually limit the deformations of Poisson by developing interline shear stresses, also as normal stresses, counting on the direction of the thickness. A transfer of interlaminar loads occurs particularly at the sides of the structure and may cause delamination.
  2. - Curved sections: In the case of curved segments, tubes, cylinders and spheres under pressure, normal and shear stresses develop at the interface of two adjacent layers. they’ll cause an interlaminar crack to develop and cause loss of adhesion.
  3. - Hygrothermal effects: The difference within the thermal expansion coefficients of the matrix and therefore the reinforcement leads to mechanical stresses under thermal gradients which may be a source of delamination also as mechanical changes. Furthermore, the anisotropic dimensional response of the laminates because of moisture absorption can also cause interline cracks. Plasticization of the polymer matrix and chemical deterioration of the constituent materials due to interaction with penetrating water may end in additional damage like cracks or cracks within the matrix.
  4. - Poor manufacturing process: Delamination can arise from the manufacturing stage on account of non-uniform resin distribution or because of the presence of voids resulting from improper practices when applying the layers. Residual stresses also can be induced by differential shrinkage of the constituents upon cooling from the cure temperature to ambient temperature.
  5. - Low speed impacts: Transverse concentrated loads caused by low-energy impacts, like a tool dropped during maintenance or propelled runway debris during take-off or landing, may cause interlaminar detachment between adjacent plies having different orientations. Impact-induced delamination is initiated because of the interaction of the matrix cracks and therefore the resin-rich zone along the ply interface. This is often likely to end in complex damage with multiple delamination, fiber rupture and matrix cracks. These damages are classified as BVID (Barely Visible Impact Damage), they arise below the surface of the laminate and are not easily detected during maintenance tasks. This explains why impact damage is so of concern when designing a composite structure, because undetected hidden delamination can cause ruptures with none external warning signs [8].

Juhász and Szekrényes [9] have studied the effect of delamination on the critical buckling force of composite plates where orthotropic rectangular plates with through-the-width delamination are modelled using special sorts of Mindlin plate finite elements. It had been found that the presence of delamination with different sizes is affecting the buckling loads. Also, for the generalization of the results, they have used an equivalent boundary conditions and layup with different plate sizes, only the ratio of the delamination and also the plate length matters.

Chen et al. [10] have developed an analytical model with a generalized Rayleigh–Ritz approach so as to review the characteristics of the buckling response of VAT (Variable angle tow) composite plates with the width, or an embedded rectangular delamination under axial compressive loading. Numerical results of VAT composite plates with one delamination were obtained with reference to various delamination sizes and positions. It had been found that the in-plane deformation of delaminated portions within the delamination region is vitally important for the delamination buckling analysis. The buckling loads decrease with a rise in delamination size.

Monda and Ramachandra [11] have investigated numerically the nonlinear dynamic pulse buckling of imperfect composite plate with embedded delamination, the dynamic buckling load is calculated using Tsai-Wu quadratic interaction criterion. The results have showed that the varied sort of impulsive loading, plate condition affects the pulse global buckling load of the structure. Moreover, the delamination may arise reducing the stiffness of the plate which results in the buckling of the structure then its collapse.

The shape of the various delamination resulting from an impact is complex and depends on the geometry of the structure, the properties of the fabric and therefore the impact energy. The active damage mechanisms are more complex than within the case of simple delamination, but the geometry observed during experimental tests is often incorporated or simplified during finite element simulation. The delamination zone is often taken as a square, circular or elliptical domain so as to confirm a satisfactory compromise between the important realistic representation of the geometry of the real delamination and also the simple insertion of the artificial damage. The placement of non-adhesive inserts of known shape and position in the structures allows the effect of delamination geometry on mechanical properties to be studied in a more controlled manner. In the study presented in this paper, we have considered only the case of damage induced by an impact and represented by a circular delamination.

2 Delamination buckling

Delamination may result during a significant reduction within the compressive strength of a carbon fiber reinforced polymer (CFRP) composite structure. A drastic reduction in bending stiffness is additionally observed. In the presence of compressive loads [12], interline cracks can trigger the local buckling of the thinnest sub-laminate. Once buckling occurs, interline delamination can expand and further reduce the strength of the structure. Fuselage panels and upper wing skins are samples of composite components that are particularly sensitive to in-service compression buckling and shear loads [13]. Due to the complexity of the analysis of those structures, which response is extremely nonlinear, engineers are reluctant to explore all the alternatives offered and thus limit the chances of design optimization. Especially, analysis within the field of post-buckling, during which nonlinearities play a dominant role, is extremely laborious [14]. What characterizes delamination buckling is the complex interaction between integrated delamination and therefore the response to buckling and post-buckling. This results in the reduction in compressive strength as has been demonstrated experimentally in damaged laminates. Therefore, the stresses in these layers are greater than those that would exist in an undamaged panel, which greatly reduces the breaking load [13].

When a delaminated composite panel is subjected to uniaxial compression within the plane, different buckling modes may develop counting on the dimensions and position of delamination. These modes are illustrated within Fig. 1.

Fig. 1.
Fig. 1.

Buckling modes of a delaminated composite panel

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

Local buckling, as shown in Fig. 1(a), may occur for instance, when the upper sub-laminate is thin and therefore the delaminated area is large. In this case, the stiffness of the upper sub-laminate is low compared to that of the lower sub-laminate and, therefore, it undergoes buckling. the appliance of a further load increment to the local buckling load may trigger a mode change to the mixed buckling mode as shown in Fig. 1(b). This mode may be a combination of local and global buckling modes and should occur before the critical load.

Depending on the case, the overall instability of the structure also can be observed consistent with the overall buckling mode, Fig. 1(c). This mode is characterized by the two layers buckling within the same direction and with an equivalent out-of-plane movement. It is going to be the primary to occur when the initial delamination features a small area and is found near the median surface of the panel. In some cases, the failure load is often reached without spreading the built-in damage.

2.1 Numerical model of delamination: cohesive elements

Several numerical techniques have been proposed to assess the problem of delamination in composite structures. Traditional numerical tools are formulated as part of linear fracture mechanics (LMR). This technique is predicated on Griffith’s fracture theory [15] which, consistent with the first law of thermodynamics, postulates that the reduction in strain energy because of crack propagation is employed to form new crack surfaces. This hypothesis is valid for brittle materials, during which the dissipation of the energy derived from the plastic deformation, during the rupture, are often neglected.

Since the nonlinear crack-tip processes happen in a very small plastic area, compared to the tiniest characteristic dimension of the structure, the numerical approaches supported the LMR (Linear Mechanics of Rupture) assume that the fracture mechanisms are often associated with the simple propagation of the delamination front. The growth of the discontinuity occurs when a mixture of the components of the rate of energy release G is such that is equal to or greater than a critical value G C .

Techniques like J-integral, EDT (Energy Derivative Technique), tangent stiffness, or VCCT (Virtual Crack Closure Technique) are employed to calculate the components of the critical rate of energy release G C using the principles of LMR. A more modern approach to delamination analysis consists of modelling, in which nonlinear crack-tip processes are represented explicitly instead of being considered to be infinitely localized on the discontinuity front [16, 17]. The formation of a cohesive zone on a spread surface of displacement discontinuity characterizes the progressive degradation of stiffness associated with the irreversible damage that occurs between the laminates.

This approach seems to permit an improved description of the physical mechanisms that develop during cracking by delamination and overcoming some important drawbacks related to models supported the LMR. Nevertheless, the hypothesis of the self-similar spread of delamination, namely that the delamination front does not change its shape throughout the loading history, requires prior knowledge of the location of the crack and therefore the direction of propagation. This requirement prevents VCCT from getting used for several categories of delamination problems like delamination caused by low-speed impact, because it cannot accurately predict the onset of damage [18, 19].

The formulation of the cohesive elements is based on the CZM (Cohesive Zone Model) model proposed by Dugdale and Barenblatt [20] to simulate complex fracture mechanisms at the crack front. The main advantage of cohesive zone models is the ability to predict the onset and spread of delamination without first knowing the location and direction of defect propagation. Therefore, unlike VCCT, problems such as the study of the compressive behavior of composite plates containing several built-in artificial delamination as well as the analysis of fracture of composite joints can be investigated numerically. For this reason, the formulation of the cohesive zone has become a very useful tool in the design of damage tolerant composite structures.

The cohesive zone approach models an extended cohesive zone, or process zone, at interfaces where delamination may occur, in which tractions or cohesive forces resist interfacial separations, often referred to as relative displacements in the literature. Indeed, the cohesive damage zone is the part of the cohesive layer closest to the delamination front in which any irreversible degradation of the properties of the interface takes place, Fig. 2. The elements within this zone are characterized by meeting the specified damage initiation criterion which governs the start of the progressive damage process. Physically, the zone of cohesive damage represents the way in which the rigidity of the material degrades locally due to the combination of the cracks around the crack tip (Fig. 3).

Fig. 2.
Fig. 2.

Cohesive zone model schematic

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

Fig. 3.
Fig. 3.

Geometry of the 8-node cohesive element

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

According to Camanho and Davila [21], the cohesive zone approaches can be related to Griffith’s fracture theory if the size of the cohesive zone is small compared to the characteristic dimensions of the structure. This condition is met by materials that exhibit near-brittle fracture behavior, such as Polyether-ether-ketone (PEEK) or epoxy-based composites used in aerospace structures. The connection between the two theories can be established by setting the area under the tensile-separation curve equal to the corresponding fracture toughness, regardless of its shape. In case of delamination under single-mode loading, we have:
0 δ 3 f τ 3 ( δ 3 ) d δ 3 = G I C ; 0 δ 2 f τ 2 ( δ 2 ) d δ 2 = G I I C ; 0 δ 1 f τ 1 ( δ 1 ) d δ 1 = G I I I C
Where GIC, GIIC et GIIIC are respectively mode I, mode II and mode III fracture toughness.

δ 3 f , δ 2 f and δ 1 f are the corresponding final separations relative to the crack propagation under single-mode loading.

In general, we can write the expression related to delamination under mixed mode loading as follow:
0 δ m f τ m ( δ m ) d δ m = G C
  • G C is the mixed-mode fracture toughness

  • δ m is the equivalent separation

  • τ m is the corresponding equivalent traction

  • δ m f is the final equivalent separation associated to complete debonding.

The so-called zero thickness elements are suitable in situations where the intermediate glue material is extremely thin and for all practical purposes are often considered zero thickness, like the bonded composite laminates investigated during this work. In this case, the macroscopic properties of the fabric are not directly relevant, and that we must use concepts derived from fracture mechanics.

2.2 Constitutive law

The constitutive law which relates the tractions of elements τ to the separations of elements δ expressed in local isoparametric coordinates is:
τ s = D s r δ r
The following assumptions are taken into account when defining the constitutive operator D sr :
  1. - Linear elastic response prior to damage onset.
τ s = K p δ s
Where Kp is the penalty stiffness.
  1. - Failure of the cohesive element is characterized by progressive degradation of the material stiffness, which is driven by a damage process.
  2. - The cohesive layer does not undergo damage under pure normal compressive stresses or strains.
A single variable constitutive law for delamination in mixed mode is able to follow the evolution of interface damage. The irreversible softening behavior of the cohesive zone and contact problems resulting from compression can be observed by defining the constitutive operator D sr as follows:
D s r = { δ s r ¯ K p δ s r ¯ K p [ ( 1 d ) K p + d K p δ 3 δ 3 δ s 3 ¯ ] δ s 3 ¯ δ s 3 ¯ δ 3 δ 3 K p
where δ s r ¯ is the Kronecker Delta. Three different areas of material behavior delimited by equivalent separations corresponding to the appearance of damage δ m 0 and total decohesion δ m f can be identified:
  1. - Linear elastic response before the onset of damage (5.a)
  2. - Gradual softening during additional loading after the onset of damage (5.b)
  3. - Control of the propagation of damage following the ultimate failure of the element (5.c)

It is noted that the constituent equation is coupled with the law of specified damage evolution through the values of critical separations δ m 0 and δ m f , which can be calculated by defining respectively the appropriate criteria for occurrence and propagation of damage, and through the variable of global damage d, which depends on the form of the law of softening as well as on δ m 0 and δ m f themselves. The next part discusses an appropriate damage law defined in terms of equivalent separations δ m and corresponding equivalent tractions τ m as well as its effect on the critical buckling load.

  1. A. Damage evolution law:

The law of evolution of the damages to be coupled with the constitutive equation in order to describe the behaviour at fracture of the cohesive layer is defined by:
  1. - A criterion for the appearance of damage: The onset of damage refers to the onset of degradation of the stiffness of the cohesive member [22]. The softening process begins when the stresses and/or strains meet a predefined damage initiation criterion. The corresponding values of equivalent separation and equivalent traction are respectively δ m 0 and τ m 0 .
  2. - In mixed mode loading, the onset of softening behaviour can occur before any of the tensile components reach their own unique allowable mode, namely interlaminar tensile strength N and interlaminar shear forces S and T. Therefore, the criterion for the occurrence of damage must take into account the interaction between normal loads and shear loads. Hence, the criterion for the appearance of damage must take into account the interaction between normal loads and shear loads. In this work, the quadratic constraint criterion proposed by [23] was used because it has been shown to provide reasonable predictions for composite structures [24]:
( τ 3 N ) 2 + ( τ 2 S ) 2 + ( τ 1 T ) 2 = 1
Among the different ways to measure the mode ratio under the mixed mode loading conditions, the mix mode ratio β has been used. For an opening separation δ 3 greater than zero, β is defined as:
β = δ c i s a i l l e m e n t δ 3
Using the quadratic stress criterion, the equivalent mixed-mode separation corresponding to the initiation of damage is given by
δ m 0 = { δ 3 0 δ 1 0 1 + β 2 ( δ 1 0 ) 2 + ( β δ 3 0 ) 2 δ 3 > 0 (8a) δ c i s a i l l e m e n t 0 δ 3 0 (8b)
Using the same penalty stiffness for modes I, II and III and assuming isotropic shear behaviour, i.e., S = T, equivalent separations δ 3 0 , δ 1 0 and δ c i s a i l l e m e n t 0 corresponding to the onset of softening under single-mode loading are equal to:
δ 3 0 = N K p , δ 1 0 = δ 2 0 = δ c i s a i l l e m e n t 0 = S K p
Finally, the equivalent traction τ m 0 corresponding to the damage initiation can be calculated as follows:
τ m 0 = K p δ m 0
  1. B. Softening law:

The formulation of the cohesive finite elements is predicated on the CZM approach with the bilinear traction-separation law. Their shape and corresponding parameters are usually determined empirically according to the expected behaviour of the adhesive. Indeed, it is widespread to use a bilinear softening (i.e. triangular) law shape for brittle adhesive and an elasto-plastic (i.e. trapezoidal) shape for ductile adhesive. Among the numerous softening models commonly used (bilinear, exponential, perfectly plastic), the bilinear law is utilized to simulate the various interlaminar toughness tests was adopted within the present study, Fig. 4.

Fig. 4.
Fig. 4.

Linear damage evolution

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

The line AB represents the linear softening envelope of the bilinear constitutive equation. It is assumed that unloading after the appearance of damage occurs linearly towards the origin of tension-separation plane. Reloading after unloading is additionnally performed along an equivalent linear path until the softening envelope is reached. A further reload follows this envelope as indicated by the arrow until it is reached.

Assuming a linear softening law, the global damage variable d is reduced to the expression proposed by Camanho [21]:
d = δ m f ( δ m m a x m 0 ( ) ) δ m m a x m m f 0 , d [ 0 , 1 ]
  1. C. Damage propagation criteria:

The definition of an appropriate criterion of propagation of the damages makes it possible to determine the equivalent separation δ m corresponding to the reduction to zero of the tractions which hold together the faces of the cohesive element. According to the principles of linear elastic fracture mechanics, the growth of mixed-mode delamination is predicted when G T G C .

Several laws have been implemented in finite element codes to account for the mixing of modes when calculating the toughness of a material, such as the power law’s criterion [25], the B-K’s criterion [26] and the Reeder’s criterion. Indeed, the predictions provided by the B-K criterion are well in line with the experimental results, as this criterion makes it possible to predict the toughness in mixed mode of materials with the same critical breaking energies for both shear modes. For the linear softening law, the explicit B-K criterion in terms of the mixture ratio of modes is given by:
G I C + ( G I I C G I C ) ( β 2 1 + β 2 ) η = G C

The parameter η is found by the least squares fitting of the experimental data obtained from a mixed-mode bending test performed at different mode ratios.

3 Finite element modelling of the delamination buckling of a composite panel

In this work, we will discuss finite element modeling of delamination buckling of composite panels using Abaqus software. A 3D model with an 8-node composite shell member is used. The panel is divided into two sub-laminates by a plane containing the delamination. The two sublayers are modeled separately using an 8-node composite shell member. Appropriate stress conditions are added for the nodes in the safe region. A set of finite element models has been implemented in Abaqus to predict delamination growth and damage evolution while observing their effects on critical buckling load.

The delamination analysis was performed by finite elements using Abaqus software (version 6.19). There are several ways to model the panel for delamination analysis. For the present study, the technique of cohesive elements was adopted. The panel is divided into two regions, the first is the intact region and the other is the cohesive region.

The reduced integration 4 node S4R shell element was used. The reference surface of the sub-laminates was moved from the median surface by using the option OFFSET of the menu of the properties of the shells that makes it possible to model the cohesive contact.

Two levels of mesh refinement were considered in the panel, as shown in Fig. 5.

  • Intact region: where the propagation of delamination is not expected. The upper and lower layers of the shell elements are linked together by beam-type multipoint stresses making it possible to simulate the contact connection.

  • Cohesive zone: where the propagation of delamination can occur. A layer of cohesive elements was placed between the top and bottom of the shell elements to simulate the evolution of damage in the interlaminar interface.

Fig. 5.
Fig. 5.

Mesh patterns of the studied panel

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

The solution obtained by using cohesive elements in the simulation of delamination mechanisms is highly dependent on the refinement of the mesh in the cohesive zone region. Indeed, the size of the element in this boundary area determines how the delamination front can extend. In addition, in order to provide realistic predictions of the evolution of delamination, the discretization of the cohesive zone must be fine enough to ensure an accurate representation of the interlaminar stress field in front of the crack tip. For these reasons, research on the dependence of the simulated response to the density of the mesh is necessary to select the characteristic size of the cohesive elements (Fig. 6).

Fig. 6.
Fig. 6.

Finite element model of the considered panel

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

A square panel 100 mm × 100 mm with a central circular delamination is considered. The delamination radius R is 40 mm as shown in Fig. 7. Due to panel symmetry, only a quarter of the panel is modeled for finite element analysis. The position of the delamination through the thickness is considered at the midplane of the laminate (Fig. 8).

Fig. 7.
Fig. 7.

Geometry of the studied panel

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

Fig. 8.
Fig. 8.

Position of delamination through the thickness in the composite panel

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

The material used for the numerical simulations is a carbon fiber reinforced Epoxy prepreg composite with a 45-degree ply orientation and, is provided by ILSB (Institute of Lightweight Design and Structural Biomechanics (E317))[27], E is the Young’s modulus, G is the shear modulus,
E 1 = E 2 = 57280.65 M P a , ν 12 = 0.0452 and G 12 = G 13 = G 23 = 4205.64 M P a

As well as the material for the cohesive zone as indicated in the subroutine presented above. The load and boundary conditions corresponding to the uniaxial loading as shown in Fig. 9. We design by u the displacements and θ the rotations.

Fig. 9.
Fig. 9.

Boundary conditions and loads

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

The aim of this study is to see the effect of behavioral parameters of cohesive elements on the ultimate buckling load in the presence of delamination. For this reason, we have varied a number of parameters such as the Damage initiation criteria and the Quads damage propagation criteria. For each criterion, we took three levels as shown in Table 1.

Table 1.

Levels of the considered criteria

Initiation damage criteria Propagation damage criteria
0.9 0.9
1 1
1.1 1.1

4 Results and discussion

In order to perform a full analysis of the structure, a first linear analysis is considered to determine the Eulerian buckling modes, Fig. 10. Subsequently, a nonlinear analysis of the structure is performed using the Riks algorithm. Initial geometric imperfections can also be introduced into the model.

Fig. 10.
Fig. 10.

The first four buckling modes of the composite panel in the presence of a delamination defect

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

The analysis makes it possible to deduce the critical buckling load for a perfect composite plate, i.e. delaminated but without having considered initial geometric imperfections, but also for a real plate suffering from imperfections.

Three levels of imperfections were considered: 10%, 20% and 30%, all modal according to the first Eulerian buckling mode of the delaminated composite plate. By recording the value of the limit load, the results obtained are shown in Table 2. This table shows that the critical buckling load decreases as a function of the magnitude of the initial geometric imperfection, which is quite expected. Figure 11 shows an example of the evolution of the curve giving the load as a function of the shortening for the three levels of initial geometric imperfection which were considered.

Table 2.

Levels of the considered criteria

Imperfection's level Buckling load (N)
10% 11,553
20% 11,136
30% 10,699
Fig. 11.
Fig. 11.

Force - displacement curve for a composite plate for three different levels of initial imperfection

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

In order to examine the influence of the cohesive behavior on the buckling resistance, we performed a parametric study in which we varied the damage criteria defined in Table 1 according to a full factorial design containing 32 = 9 combinations. For each case, we calculated the buckling load in the presence of delamination. The results obtained are given in Table 3.

Table 3.

Delamination buckling simulation results

Quads Damage Criterion (Qdc) Damage Evolution (De) Buckling load (N)
0.9 0.9 43,165
0.9 1 45,109
0.9 1.1 50,184
1 0.9 47,099
1 1 48,554
1 1.1 54,396
1.1 0.9 47,795
1.1 1 57,449
1.1 1.1 58,498

From Table 3, we can notice that the critical load varies with the two criteria that control the behaviour of cohesive elements. As expected, the critical load increases with the damage initiation threshold and the stiffness of damage propagation. Figure 12 shows the delamination buckling of the composite plate where the onset of delamination growth occurs at a slightly lower load than that associated with the overall buckling. The failure of the delaminated plate occurs during the propagation of the delamination crack.

Fig. 12.
Fig. 12.

Buckling of the plate in presence of delamination

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00354

From Table 3, it is possible to build a polynomial response surface allowing the buckling load to be expressed as a function of the parameters of the cohesive behavior. Firstly, we’ve performed an analysis of variance on the results obtained. We have started with a linear model (Table 4), we can notice that it is the initiation threshold factor that has the greatest influence on the variability of the buckling resistance.

Table 4.

Delamination buckling simulation results

Source Sum Sq. d.f. Mean Sq. F Prob>F
Qdc 1.07E+08 2 5.34E+07 10.02 2.77E-02
De 1.04E+08 2 5.22E+07 9.79 0.288
Error 2.13E+07 4 5.33E+06
Total 2.33E+08 8
The linear model obtained by the “regstats” command of Matlab is written in the form:
P c r = 10 4 × ( 0.35905 + 4.21408 × Q d c + 4.16995 × D e )
R 2 = 90.7 %

We can also derive with an interacting model to improve the representation of the buckling load and increase the coefficient of determination R 2. The results of the analysis of variance are given in Table 5.

Table 5.

Delamination buckling simulation results

Source Sum Sq. d.f. Mean Sq. F Prob>F
Qdc 1.07E+08 2 5.34E+07 Inf NaN
De 1.04E+08 2 5.22E+07 Inf NaN
Qdc*De 2.13E+07 4 5.33E+06 Inf NaN
Error −1.49E−07 0 0
Total 2.33E+08 8
The obtained interacting model is given by:
P c r = 10 4 × ( 5.85295 4.99792 × Q d c 5.04205 × D e + 9.212 × Q d c × D e )
R 2 = 92.15 %
Finally, we can interpolate the results obtained using a quadratic polynomial. We then find:
P c r = 10 4 × ( 7.53268 11.9989 × Q d c 1.42305 × D e + 9.212 × Q d c × D e + 3.5005 × Q d c 2 1.8095 × D e 2 )
R 2 = 92.3 %

We note that the quadratic model is complex but that it does not provide more information on the correlation than the interacting model. The coefficient of determination R2 does not improve significantly compared to the interacting model.

We also notice that the linear model has the merit of showing how the buckling resistance increases with the material performances of the cohesive zone. This is not as obvious in the case of the other two models.

5 Conclusion

In the present work, the problem of buckling in interaction or not with delamination occurring in a composite panel subjected to axial compression was simulated. The method of cohesive elements was used. We performed a parametric study to quantify the influence of parameters describing the cohesive behavior on buckling resistance in the presence of circular delamination. In the studied field of parameters, the results obtained showed a certain sensitivity with respect to the amplitude of the initial displacements out of planes associated with the imperfection. Buckling occurs more easily as the magnitude of the defect increases. We have also found that the two criteria of cohesive behavior play an important role. A significant increase in the maximum load is observed as the initiation threshold or the stiffness of the cohesive members increases. Further studies are needed to objectively quantify the delamination buckling interaction and in particular to analyze the influence of delamination geometry and its position on strength.

Nomenclature

Abbreviations

BVID

Barely Visible Impact Damage

CFRP

carbon fiber reinforced polymer

CZM

Cohesive Zone Model

Df

degrees of freedom

F

F-statistic is the ratio of the mean squared errors

EDT

Energy Derivative Technique

FEM

Finite Element Method

LMR

linear mechanics of rupture

Mean Sq

Mean square of the error term

PEEK

Polyether-ether-ketone

Sum Sq

Sum of Squared Elements

VAT

Variable angle tow

VCCT

Virtual Crack Closure Technique

Symbols

d

Global damage variable

Dsr

Constitutive operator

G

Energy release rate

GI

Mode I energy release rate

GII

Mode II energy release rate

GIII

Mode III energy release rate

GC

Mixed-mode fracture toughness

GIC

Mode I critical energy release rate

GIIC

Mode II critical energy release rate

GIIIC

Mode III critical energy release rate

Gcisaillement

Shear mode energy release rate

GT

Total energy release rate

Kp

Penalty stiffness

N

Interlaminar tensile strength

S and T

Interlaminar shear forces

β

Mix mode ratio

δ 3 f , δ 2 f , δ 1 f

Final separations relative to the crack propagation under single-mode loading

δ m

Equivalent separation

τ m

Corresponding equivalent traction

δ m f

Final equivalent separation associated to complete debonding

δ sr ¯

Kronecker Delta

References

  • [1]

    M. Chaves-Vargas , a. Dafnis , H. G. Reimerdes , and K. U. Schröder , “Modal parameter identification of a compression-loaded CFRP stiffened plate and correlation with its buckling behaviour,” Prog. Aerosp. Sci., vol. 78, pp. 3949, 2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [2]

    E. S. Greenhalgh , “Delamination-dominated failures in polymer composites,” Fail. Anal. Fractography Polym. Compos., pp. 164237, 2009.

    • Search Google Scholar
    • Export Citation
  • [3]

    N. J. Pagano and G. A. Schoeppner , “Delamination of polymer matrix composites: problems and assessment,” Compr. Compos. Mater., pp. 433528, 2000.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [4]

    I. Feddal , A. Khamlichi , and K. Ameziane , “Resistance to buckling of a stiffened panel under axial compression; Effect of imperfections resulting from welding process,” Proced. Manufacturing, vol. 32, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [5]

    A. Shahrjerdi , and B. Bahramibabamiri , “The effect of different geometrical imperfection of buckling of composite cylindrical shells subjected to axial loading,” Int. J. Mech. Mater. Eng., vol. 10, no. 1, 2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [6]

    R. Degenhardt , K. Rohwer , D. Wilckens , and J. Tessmer , “Improved material exploitation of composite airframe structures by accurate simulation of collapse - the COCOMAT project,” Tech. Isr. Inst. Technol. - 48th Isr. Annu. Conf. Aerosp. Sci. 2008, vol. 1, no. March, pp. 458463, 2008.

    • Search Google Scholar
    • Export Citation
  • [7]

    N. P. , M. T. H. Ku , and H. Wang , A Review on the Tensile Properties of Natural Fibre Reinforced Polymer Composites. University of Southern Queensland, 2014.

    • Search Google Scholar
    • Export Citation
  • [8]

    M. J. Pavier , and M. P. Clarke , “Experimental techniques for the investigation of the effects of impact damage on carbon-fibre composites,” Compos. Sci. Technol., vol. 55, no. 2, pp. 157169, 1995.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [9]

    Z. Juhász and A. Szekrényes , “The effect of delamination on the critical buckling force of composite plates: experiment and simulation,” Compos. Struct., vol. 168, pp. 456464, 2017.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [10]

    X. Chen , Z. Wu , G. Nie , and P. Weaver , “Buckling analysis of variable angle tow composite plates with a through-the-width or an embedded rectangular delamination,” Int. J. Sol. Struct., vol. 138, pp. 166180, 2018.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [11]

    S. Mondal , and L. S. Ramachandra , “Nonlinear dynamic pulse buckling of imperfect laminated composite plate with delamination,” Int. J. Sol. Struct., vol. 198, pp. 170182, 2020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [12]

    A. A. Griffits , “The phenomena of rupture and flow in solids,” Masinovedenie, no. 1, pp. 914, 1995.

  • [13]

    A. A. Griffith , “The Phenomena of Rupture and Flow in Solids VI . The Phenomena of Rupture and Flow in Solids . larger question of the nature of intermolecular cohesion . The original object of the work , which was carried out at the Royal Aircraft Estab- lishment, was t,” Philos. Trans. R. Soc. Lond. Ser. A, Contain. Pap. a Math. or Phys. Character, vol. 221, no. 1921, pp. 163198, 1921.

    • Search Google Scholar
    • Export Citation
  • [14]

    D. Albiol , Buckling Analyses of Composite Laminated Panels with Delamination. Facolta di Ingegneria Industriale - Politecnico Di Milano, 2010.

    • Search Google Scholar
    • Export Citation
  • [15]

    A. C. Orifici , R. S. Thomson , R. Degenhardt , C. Bisagni , and J. Bayandor , “Development of a degradation model for the collapse analysis of composite aerospace structures,” Proc. 3rd Eur. Conf. Comput. Mech. Solids, Struct. Coupled Probl. Eng., no. July, 2006.

    • Search Google Scholar
    • Export Citation
  • [16]

    A. D. Dimarogonas , “Vibration of cracked structures: a state of the art review,” Eng. Fract. Mech., vol. 55, no. 5, pp. 831857, 1996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [17]

    D. S. Dugdale , “Yielding of steel sheets containing slits,” J. Mech. Phys. Sol., vol. 8, no. 2, pp. 100104, 1960.

  • [18]

    P. P. Camanho , C. G. Davila , and S. S. Pinho , “Fracture analysis of composite co-cured structural joints using decohesion elements,” Fatigue Fract. Eng. Mater. Struct., vol. 27, no. 9, pp. 745757, 2004.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [19]

    F. Marulo , M. Guida , L. Maio , and F. Ricci , “Numerical simulations and experimental experiences of impact on composite structures,” in Dynamic Response and Failure of Composite Materials and Struc tures, V. Lopresto , L. Antionio , and A. Serge , Eds, Woodhead Publishing, 2017, pp. 85125.

    • Search Google Scholar
    • Export Citation
  • [20]

    M. R. Wisnom , Z. J. Petrossian , and M. I. Jones , “Interlaminar failure of unidirectional glass/epoxy due to combined through thickness shear and tension,” Compos. Part A. Appl. Sci. Manuf., vol. 27, no. 10, pp. 921929, 1996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [21]

    P. Camanho , and C. G. Davila , “Mixed-mode decohesion finite elements in for the simulation composite of delamination materials,” Nasa, vols TM-2002-21, no. June, pp. 137, 2002.

    • Search Google Scholar
    • Export Citation
  • [22]

    M. Kenane , and M. L. Benzeggagh , “Mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites under fatigue loading,” Compos. Sci. Technol., vol. 57, no. 5, pp. 597605, 1997.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [23]

    J. Edward , M. Wu , and R. C. Router , Crack Extension in Fiberglass Reinforced Plastics. Urbana, Illinois, 1965, Contract N° 64-0178-d.

    • Search Google Scholar
    • Export Citation
  • [24]

    G. Marom , Environmental Effects on Fracture Mechanical Properties of Polymer Composites, vol. 6, no. C. Elsevier Science Publishers B.V., 1989.

    • Search Google Scholar
    • Export Citation
  • [25]

    A. J. Brunner , B. R. K. Blackman , and P. Davies , “A status report on delamination resistance testing of polymer-matrix composites,” Eng. Fract. Mech., vol. 75, no. 9, pp. 27792794, 2008.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [26]

    N. F. Knight , “User-defined material model for thermo- mechanical progressive failure analysis,” no. December, 2008.

  • [27]

    T. Ceglar , “Finite element method simulations of the delamination in laminated composite structures,” no. February, 2017.

  • [1]

    M. Chaves-Vargas , a. Dafnis , H. G. Reimerdes , and K. U. Schröder , “Modal parameter identification of a compression-loaded CFRP stiffened plate and correlation with its buckling behaviour,” Prog. Aerosp. Sci., vol. 78, pp. 3949, 2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [2]

    E. S. Greenhalgh , “Delamination-dominated failures in polymer composites,” Fail. Anal. Fractography Polym. Compos., pp. 164237, 2009.

    • Search Google Scholar
    • Export Citation
  • [3]

    N. J. Pagano and G. A. Schoeppner , “Delamination of polymer matrix composites: problems and assessment,” Compr. Compos. Mater., pp. 433528, 2000.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [4]

    I. Feddal , A. Khamlichi , and K. Ameziane , “Resistance to buckling of a stiffened panel under axial compression; Effect of imperfections resulting from welding process,” Proced. Manufacturing, vol. 32, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [5]

    A. Shahrjerdi , and B. Bahramibabamiri , “The effect of different geometrical imperfection of buckling of composite cylindrical shells subjected to axial loading,” Int. J. Mech. Mater. Eng., vol. 10, no. 1, 2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [6]

    R. Degenhardt , K. Rohwer , D. Wilckens , and J. Tessmer , “Improved material exploitation of composite airframe structures by accurate simulation of collapse - the COCOMAT project,” Tech. Isr. Inst. Technol. - 48th Isr. Annu. Conf. Aerosp. Sci. 2008, vol. 1, no. March, pp. 458463, 2008.

    • Search Google Scholar
    • Export Citation
  • [7]

    N. P. , M. T. H. Ku , and H. Wang , A Review on the Tensile Properties of Natural Fibre Reinforced Polymer Composites. University of Southern Queensland, 2014.

    • Search Google Scholar
    • Export Citation
  • [8]

    M. J. Pavier , and M. P. Clarke , “Experimental techniques for the investigation of the effects of impact damage on carbon-fibre composites,” Compos. Sci. Technol., vol. 55, no. 2, pp. 157169, 1995.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [9]

    Z. Juhász and A. Szekrényes , “The effect of delamination on the critical buckling force of composite plates: experiment and simulation,” Compos. Struct., vol. 168, pp. 456464, 2017.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [10]

    X. Chen , Z. Wu , G. Nie , and P. Weaver , “Buckling analysis of variable angle tow composite plates with a through-the-width or an embedded rectangular delamination,” Int. J. Sol. Struct., vol. 138, pp. 166180, 2018.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [11]

    S. Mondal , and L. S. Ramachandra , “Nonlinear dynamic pulse buckling of imperfect laminated composite plate with delamination,” Int. J. Sol. Struct., vol. 198, pp. 170182, 2020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [12]

    A. A. Griffits , “The phenomena of rupture and flow in solids,” Masinovedenie, no. 1, pp. 914, 1995.

  • [13]

    A. A. Griffith , “The Phenomena of Rupture and Flow in Solids VI . The Phenomena of Rupture and Flow in Solids . larger question of the nature of intermolecular cohesion . The original object of the work , which was carried out at the Royal Aircraft Estab- lishment, was t,” Philos. Trans. R. Soc. Lond. Ser. A, Contain. Pap. a Math. or Phys. Character, vol. 221, no. 1921, pp. 163198, 1921.

    • Search Google Scholar
    • Export Citation
  • [14]

    D. Albiol , Buckling Analyses of Composite Laminated Panels with Delamination. Facolta di Ingegneria Industriale - Politecnico Di Milano, 2010.

    • Search Google Scholar
    • Export Citation
  • [15]

    A. C. Orifici , R. S. Thomson , R. Degenhardt , C. Bisagni , and J. Bayandor , “Development of a degradation model for the collapse analysis of composite aerospace structures,” Proc. 3rd Eur. Conf. Comput. Mech. Solids, Struct. Coupled Probl. Eng., no. July, 2006.

    • Search Google Scholar
    • Export Citation
  • [16]

    A. D. Dimarogonas , “Vibration of cracked structures: a state of the art review,” Eng. Fract. Mech., vol. 55, no. 5, pp. 831857, 1996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [17]

    D. S. Dugdale , “Yielding of steel sheets containing slits,” J. Mech. Phys. Sol., vol. 8, no. 2, pp. 100104, 1960.

  • [18]

    P. P. Camanho , C. G. Davila , and S. S. Pinho , “Fracture analysis of composite co-cured structural joints using decohesion elements,” Fatigue Fract. Eng. Mater. Struct., vol. 27, no. 9, pp. 745757, 2004.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [19]

    F. Marulo , M. Guida , L. Maio , and F. Ricci , “Numerical simulations and experimental experiences of impact on composite structures,” in Dynamic Response and Failure of Composite Materials and Struc tures, V. Lopresto , L. Antionio , and A. Serge , Eds, Woodhead Publishing, 2017, pp. 85125.

    • Search Google Scholar
    • Export Citation
  • [20]

    M. R. Wisnom , Z. J. Petrossian , and M. I. Jones , “Interlaminar failure of unidirectional glass/epoxy due to combined through thickness shear and tension,” Compos. Part A. Appl. Sci. Manuf., vol. 27, no. 10, pp. 921929, 1996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [21]

    P. Camanho , and C. G. Davila , “Mixed-mode decohesion finite elements in for the simulation composite of delamination materials,” Nasa, vols TM-2002-21, no. June, pp. 137, 2002.

    • Search Google Scholar
    • Export Citation
  • [22]

    M. Kenane , and M. L. Benzeggagh , “Mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites under fatigue loading,” Compos. Sci. Technol., vol. 57, no. 5, pp. 597605, 1997.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [23]

    J. Edward , M. Wu , and R. C. Router , Crack Extension in Fiberglass Reinforced Plastics. Urbana, Illinois, 1965, Contract N° 64-0178-d.

    • Search Google Scholar
    • Export Citation
  • [24]

    G. Marom , Environmental Effects on Fracture Mechanical Properties of Polymer Composites, vol. 6, no. C. Elsevier Science Publishers B.V., 1989.

    • Search Google Scholar
    • Export Citation
  • [25]

    A. J. Brunner , B. R. K. Blackman , and P. Davies , “A status report on delamination resistance testing of polymer-matrix composites,” Eng. Fract. Mech., vol. 75, no. 9, pp. 27792794, 2008.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [26]

    N. F. Knight , “User-defined material model for thermo- mechanical progressive failure analysis,” no. December, 2008.

  • [27]

    T. Ceglar , “Finite element method simulations of the delamination in laminated composite structures,” no. February, 2017.

  • Collapse
  • Expand

Senior editors

Editor-in-Chief: Ákos, LakatosUniversity of Debrecen, Hungary

Founder, former Editor-in-Chief (2011-2020): Ferenc Kalmár, University of Debrecen, Hungary

Founding Editor: György Csomós, University of Debrecen, Hungary

Associate Editor: Derek Clements Croome, University of Reading, UK

Associate Editor: Dezső Beke, University of Debrecen, Hungary

Editorial Board

  • Mohammad Nazir AHMAD, Institute of Visual Informatics, Universiti Kebangsaan Malaysia, Malaysia

    Murat BAKIROV, Center for Materials and Lifetime Management Ltd., Moscow, Russia

    Nicolae BALC, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Umberto BERARDI, Toronto Metropolitan University, Toronto, Canada

    Ildikó BODNÁR, University of Debrecen, Debrecen, Hungary

    Sándor BODZÁS, University of Debrecen, Debrecen, Hungary

    Fatih Mehmet BOTSALI, Selçuk University, Konya, Turkey

    Samuel BRUNNER, Empa Swiss Federal Laboratories for Materials Science and Technology, Dübendorf, Switzerland

    István BUDAI, University of Debrecen, Debrecen, Hungary

    Constantin BUNGAU, University of Oradea, Oradea, Romania

    Shanshan CAI, Huazhong University of Science and Technology, Wuhan, China

    Michele De CARLI, University of Padua, Padua, Italy

    Robert CERNY, Czech Technical University in Prague, Prague, Czech Republic

    Erdem CUCE, Recep Tayyip Erdogan University, Rize, Turkey

    György CSOMÓS, University of Debrecen, Debrecen, Hungary

    Tamás CSOKNYAI, Budapest University of Technology and Economics, Budapest, Hungary

    Anna FORMICA, IASI National Research Council, Rome, Italy

    Alexandru GACSADI, University of Oradea, Oradea, Romania

    Eugen Ioan GERGELY, University of Oradea, Oradea, Romania

    Janez GRUM, University of Ljubljana, Ljubljana, Slovenia

    Géza HUSI, University of Debrecen, Debrecen, Hungary

    Ghaleb A. HUSSEINI, American University of Sharjah, Sharjah, United Arab Emirates

    Nikolay IVANOV, Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia

    Antal JÁRAI, Eötvös Loránd University, Budapest, Hungary

    Gudni JÓHANNESSON, The National Energy Authority of Iceland, Reykjavik, Iceland

    László KAJTÁR, Budapest University of Technology and Economics, Budapest, Hungary

    Ferenc KALMÁR, University of Debrecen, Debrecen, Hungary

    Tünde KALMÁR, University of Debrecen, Debrecen, Hungary

    Milos KALOUSEK, Brno University of Technology, Brno, Czech Republik

    Jan KOCI, Czech Technical University in Prague, Prague, Czech Republic

    Vaclav KOCI, Czech Technical University in Prague, Prague, Czech Republic

    Imre KOCSIS, University of Debrecen, Debrecen, Hungary

    Imre KOVÁCS, University of Debrecen, Debrecen, Hungary

    Angela Daniela LA ROSA, Norwegian University of Science and Technology, Trondheim, Norway

    Éva LOVRA, Univeqrsity of Debrecen, Debrecen, Hungary

    Elena LUCCHI, Eurac Research, Institute for Renewable Energy, Bolzano, Italy

    Tamás MANKOVITS, University of Debrecen, Debrecen, Hungary

    Igor MEDVED, Slovak Technical University in Bratislava, Bratislava, Slovakia

    Ligia MOGA, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Marco MOLINARI, Royal Institute of Technology, Stockholm, Sweden

    Henrieta MORAVCIKOVA, Slovak Academy of Sciences, Bratislava, Slovakia

    Phalguni MUKHOPHADYAYA, University of Victoria, Victoria, Canada

    Balázs NAGY, Budapest University of Technology and Economics, Budapest, Hungary

    Husam S. NAJM, Rutgers University, New Brunswick, USA

    Jozsef NYERS, Subotica Tech College of Applied Sciences, Subotica, Serbia

    Bjarne W. OLESEN, Technical University of Denmark, Lyngby, Denmark

    Stefan ONIGA, North University of Baia Mare, Baia Mare, Romania

    Joaquim Norberto PIRES, Universidade de Coimbra, Coimbra, Portugal

    László POKORÁDI, Óbuda University, Budapest, Hungary

    Roman RABENSEIFER, Slovak University of Technology in Bratislava, Bratislava, Slovak Republik

    Mohammad H. A. SALAH, Hashemite University, Zarqua, Jordan

    Dietrich SCHMIDT, Fraunhofer Institute for Wind Energy and Energy System Technology IWES, Kassel, Germany

    Lorand SZABÓ, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Csaba SZÁSZ, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Ioan SZÁVA, Transylvania University of Brasov, Brasov, Romania

    Péter SZEMES, University of Debrecen, Debrecen, Hungary

    Edit SZŰCS, University of Debrecen, Debrecen, Hungary

    Radu TARCA, University of Oradea, Oradea, Romania

    Zsolt TIBA, University of Debrecen, Debrecen, Hungary

    László TÓTH, University of Debrecen, Debrecen, Hungary

    László TÖRÖK, University of Debrecen, Debrecen, Hungary

    Anton TRNIK, Constantine the Philosopher University in Nitra, Nitra, Slovakia

    Ibrahim UZMAY, Erciyes University, Kayseri, Turkey

    Andrea VALLATI, Sapienza University, Rome, Italy

    Tibor VESSELÉNYI, University of Oradea, Oradea, Romania

    Nalinaksh S. VYAS, Indian Institute of Technology, Kanpur, India

    Deborah WHITE, The University of Adelaide, Adelaide, Australia

International Review of Applied Sciences and Engineering
Address of the institute: Faculty of Engineering, University of Debrecen
H-4028 Debrecen, Ótemető u. 2-4. Hungary
Email: irase@eng.unideb.hu

Indexing and Abstracting Services:

  • DOAJ
  • ERIH PLUS
  • Google Scholar
  • ProQuest
  • SCOPUS
  • Ulrich's Periodicals Directory

 

2023  
Scimago  
Scimago
H-index
11
Scimago
Journal Rank
0.249
Scimago Quartile Score Architecture (Q2)
Engineering (miscellaneous) (Q3)
Environmental Engineering (Q3)
Information Systems (Q4)
Management Science and Operations Research (Q4)
Materials Science (miscellaneous) (Q3)
Scopus  
Scopus
Cite Score
2.3
Scopus
CIte Score Rank
Architecture (Q1)
General Engineering (Q2)
Materials Science (miscellaneous) (Q3)
Environmental Engineering (Q3)
Management Science and Operations Research (Q3)
Information Systems (Q3)
 
Scopus
SNIP
0.751


International Review of Applied Sciences and Engineering
Publication Model Gold Open Access
Online only
Submission Fee none
Article Processing Charge 1100 EUR/article
Regional discounts on country of the funding agency World Bank Lower-middle-income economies: 50%
World Bank Low-income economies: 100%
Further Discounts Limited number of full waivers available. Editorial Board / Advisory Board members: 50%
Corresponding authors, affiliated to an EISZ member institution subscribing to the journal package of Akadémiai Kiadó: 100%
Subscription Information Gold Open Access

International Review of Applied Sciences and Engineering
Language English
Size A4
Year of
Foundation
2010
Volumes
per Year
1
Issues
per Year
3
Founder Debreceni Egyetem
Founder's
Address
H-4032 Debrecen, Hungary Egyetem tér 1
Publisher Akadémiai Kiadó
Publisher's
Address
H-1117 Budapest, Hungary 1516 Budapest, PO Box 245.
Responsible
Publisher
Chief Executive Officer, Akadémiai Kiadó
ISSN 2062-0810 (Print)
ISSN 2063-4269 (Online)

Monthly Content Usage

Abstract Views Full Text Views PDF Downloads
Sep 2024 0 100 16
Oct 2024 0 323 16
Nov 2024 0 384 22
Dec 2024 0 425 22
Jan 2025 0 239 16
Feb 2025 0 119 5
Mar 2025 0 0 0