The Role of Interfacial Adhesion on the Mechanical Behavior of Thin Metal/Polymer Laminate with Surface Roughness

Metal/polymer laminate has versatile applications in industry due to the essential roles of its constituents in controlling its mechanical behavior. Therefore, efforts to enhance a laminate’s performance should target its mechanical behavior. One of the most influencing features of the mechanical behavior of this type of thin laminate is the interface layer properties. This study concentrates on the mechanical response of thin aluminum (Al) foil deposited on a polymer substrate by calibrating its interfacial layer properties based on available uniaxial tensile tests performed on thin Al/polymer laminate. Then, taking into account the calibrated parameters for the interface layer, which leads to mimicking the real conditions of the laminate, one type of imperfection is introduced as a wavy roughness on the surface of each layer with different amplitudes to investigate its influence on the overall mechanical behavior of the laminate and its failure mode. The results highlighted that the existence of the roughness on the surface of the polymer layer reduces the maximum engineering stress of the laminate more severely compared to other conditions. As the roughness amplitude increases, the maximum stress reduces a lot. The distribution of equivalent plastic strains represents the appearance of the shear bands in the Al layer and an almost uniform distribution for the polymer layer. In the case of existing roughness on each layer, a higher amount of plastic strain accumulation occurs in the middle of the polymer layer and top corners of the thin Al layer. Due to the significant effect of interfacial layer properties to improve the maximum strength of the laminate and its final elongation, a parametric study is performed, taking into account different interfacial properties. The results indicate that laminate behavior with weaker separation properties in the interface layer is mostly unaffected by adopting higher tractions, and no change happens in the case of high separation considering weak tractions.


Introduction
Metal composites have attracted much attention due to their vast applications. One of the applications of metal composites is flexible electronics, which are produced by assembling electronic circuits and mounting electronic gadgets on flexible plastic substrates made of polyimide or transparent conductive polyester [1]. Applications concern flexible (rollable) displays, printable thin film solar cells, paper-like displays and electronic (artificial) skins [2]. To make the structure stretchable, all the components must comply with some degree of bending or a specific percentage of stretch without losing their capability to function. The functional appropriateness of flexible electronics strongly depends on the adhesion of the metal/polymer interface as well as on the geometry of the metal and polymer layers separately. Using numerical modeling, two competing failure mechanisms (metal film necking versus grain boundary cracking) of thin metal films on polymer substrates were studied [3]. They showed that when the grain boundaries of the metal films are weak, the thin metal films tend to rupture by grain boundary cracking. On the other experimental test behavior using a numerical method. The results of this study are crucial for a better understanding of the mechanical behavior of Al foil in Al/polymer laminate.
In laminated structures, the influences of cracks are significant, and it is essential to know the fracture characteristics of a crack along or perpendicular to the laminate interface. To study these characteristics, analytical methods were used to derive basic equations for cracks perpendicular to the interface of a finite and infinite solid [15]. The dislocation simulation approach as a series of the first Chebyshev polynomial showed that for the crack situated in weaker material, the stress intensity factor is smaller than that in the homogeneous material and larger for cracks in a stiffer material. They observed that for the finite and infinite solids of biomaterial, despite the existence of different stress intensity factors, the stress distributions were similar. Another research studied the micromechanisms of fracture in a laminate composed of an Al foil and a polymer film [16]. A micro-mechanical approach utilizing scanning electron microscope (SEM)-micrographs, micro-mechanisms and an analytical expression motivated the derivation of an equation suitable to calculate the work of failure of freestanding and laminated thin Al foil. The authors applied a tensile test on the thin laminate of Al foil and LDPE to evaluate their stress-strain relationships. With the assumption of perfect adhesion between these two layers, LDPE could increase the rupture strain of Al foil to a reasonable extent.
The mechanical and fracture behavior of a thin Al foil and polymer laminate was studied by [17]. The possibility of controlling the path of the growing crack propagation by adjusting the adhesion level and the property of the polymer layer was investigated. First, the fracture process of the Al foil was studied experimentally. Secondly, the mechanical and fracture behavior of the laminate were evaluated. A theory for the mechanics of the composite material was used to evaluate a series of experiments. Then, the behavior of a crack in one of the layers, perpendicular to the bi-material interface in a finite solid was evaluated by formulating a dislocation superposition method. The results showed that the analytical methods for an asymptotically small crack extension could also be applied for a fairly large amount of crack growth. As adopting a proper set of material parameters for each layer in the laminate leads to realistic results, the fracture properties together with damage parameters of the laminate and the adhesion between two layers were defined based on experimental and theoretical methods by [18].
The opening procedure in the packaging system was also studied in [19], which simulated a numerical model including realistic material behavior from experimental tests and considering large deformations and fracture mechanic concepts. The result of the centercracked tension specimen under fracture mode I was compared with the experimental test and showed good agreement. Various models were developed [20] to recover the adhesion properties of the Al/polymer laminate and highlighted the difficulties in characterizing the interface properties.
In the above studies, different investigations indicated the important role of interfacial adhesion properties on the mechanical behavior of the laminate, which makes it necessary to obtain those properties comparable with experimental data to be applicable in real conditions. In the current study, first, the interface properties are calibrated considering the bi-linear T-S law and based on the available uniaxial tensile test on the laminate. Then, the influence of surface roughness as a kind of imperfection on the mechanical behavior and failure of laminates is studied. The surface roughness is assumed to be sinusoidal, and the effects of its wavelength and amplitude are investigated independently. The details of the interfacial layer properties are investigated in a complete parametric analysis due to their considerable significance. The development of the numerical model of the adhesion layer with its properties like the real laminate condition is an advancement to study this kind of thin laminates' mechanical behavior numerically. With the knowledge of interface layer properties and the adoption of possible imperfection in numerical models, the results of the current study are useful in real applications.

Materials and Numerical Methods
The materials used in the current study consist of thin Al foil deposited on an LDPE layer, and their mechanical behavior together with material properties are described in the following.

Kinematics and Constitutive Law
Based on the continuum mechanic theory for elastoplastic materials behavior [21], logarithmic strain defined in the numerical model is formulated based on large deformation theory as: where λ i is the principle stretch in principle direction of N i . The logarithmic strain can be related to the rate of deformation by: .
where F is the gradient of deformation. With the assumption of small elastic strain and large plastic strain, the total strain rate can be decomposed additively as where . ε el and . ε p are elastic strain rate and plastic strain rate, respectively, in which the response with small elastic strain is formulated as where C is the elastic stiffness matrix. The elastic domain is restricted by a yield surface in the form of f (σ, σ) = 0, where σ corresponds to the evolved yield limit considering exponential hardening according to relationship where σ 0 is the initial yield surface, and n is the hardening exponent. ε p represents the equivalent plastic strain during plastic deformation and with the associated material response, plastic strains in relationship 7 are defined by the following flow rule: where .
Λ is the plastic multiplier defining the loading/unloading criterion through the complementary conditions:

Finite Element Model
To investigate the mechanical behavior of the laminate numerically, finite element models are developed that implements both material and geometrical nonlinearity. An elastic-plastic constitutive law based on the classical Huber-Hencky-Mises yield criterion is assumed for Al and LDPE. True stress versus the logarithmic strain curve of Al and LDPE, extracted from the experimental tensile tests performed by [16] are used as the material's  Table 1. The problem domain is discretized by four-node and eight-node elements with randomly distributed, fine meshes considering either two-dimensional (2D) or threedimensional (3D) quasi-static analyses, respectively. The plane stress element is used for the 2D problem, and solid elements are adopted for 3D analyses. The clamped end boundary conditions are applied on one end, and the loading is applied on the other end in the strain-controlled condition.

Material Parameters Calibration
To model the softening behavior of the laminate, an interface layer is considered between laminate layers. Since the focus here is to calibrate the interface properties, the laminate is modeled numerically considering its average mechanical behavior [16] as a homogeneous solid structure and by using plain stress conditions and neglecting its composition of Al foil and LDPE layers separately. The average mechanical properties of the laminate are listed in Table 2. The geometry of the plain and notched laminates is the same as the tested samples in experiments performed by [16] and is listed in Table 3. The discretized 2D numerical model of the laminate visualizing its boundary condition is depicted in Figure 1. The details of experimental tests on the thin laminate are discussed completely in [16].   Table 3. Geometry details of the laminate.

Length (mm) Width (mm) Thickness (mm)
The interface layer is modeled by zero-thickness cohesive elements, and its behavior is characterized by a bilinear T-S law (visualized in Figure 2), whose parameters are determined by calibrating based on the behavior of the plain laminate samples under uniaxial load in the experimental test.
The failure mechanism represented in Figure 2 is characterized by three parameters; the peak traction stress σ p , the corresponding separation displacement δ p and the ultimate separation displacement or displacement jump when decohesion occurs δ u . The same normal and shear tractions are considered in the following study, and separations are assumed to be normal to the interface and in the shear direction, respectively. The initial stiffness of the cohesive zone model should be high enough to keep a stiff connection between two neighboring bulk elements and to obtain realistic results [22]. The initial stiffness is considered to be 1000 times that of Al foil Young's modulus [23]. The maximum normal/shear strength is chosen from the maximum stress experienced in the nominal stress-strain curve of the laminate under tensile loading in experiments. The ultimate separation parameter, which describes the strain localization and the softening branch of the laminate, is determined iteratively to capture the strain localization of the plain sample acquired by a tensile machine. The calibrated interface properties are listed in Table 4.
The details of experimental tests on the thin laminate are discussed completely in [16].   Table 3. Geometry details of the laminate. The interface layer is modeled by zero-thickness cohesive elements, and its behavior is characterized by a bilinear T-S law (visualized in Figure 2), whose parameters are determined by calibrating based on the behavior of the plain laminate samples under uniaxial load in the experimental test.

Length (mm) Width (mm) Thickness (mm)
The failure mechanism represented in Figure 2 is characterized by three parameters; the peak traction stress , the corresponding separation displacement and the ultimate separation displacement or displacement jump when decohesion occurs . The same normal and shear tractions are considered in the following study, and separations are assumed to be normal to the interface and in the shear direction, respectively. The initial stiffness of the cohesive zone model should be high enough to keep a stiff connection between two neighboring bulk elements and to obtain realistic results [23]. The initial stiffness is considered to be 1000 times that of Al foil Young's modulus [24]. The maximum normal/shear strength is chosen from the maximum stress experienced in the nominal stress-strain curve of the laminate under tensile loading in experiments. The ultimate separation parameter, which describes the strain localization and the softening branch of the laminate, is determined iteratively to capture the strain localization of the plain sample acquired by a tensile machine. The calibrated interface properties are listed in Table 4.   Figure 3, in the experimental tests data, after maximum stress, the laminate first displays a softening behavior due to the degradation of the Al foil layer and then a plateau (in the case of the plain sample) and another softening branch with a different slope, which is supported seen by only a polymer layer. While in the numerical models of plain and notched laminates, after maximum stress, just one softening branch exists that is due to the modeling of the laminate based on its average mechanical behavior from experimental tests and neglecting the Al foil and LDPE layers behavior separately.   Figure 3 shows the engineering stress-strain behavior of the plain and notched laminates (with central crack lengths of 15 and 45 mm) under tensile loading by calibrating the interfacial properties based on the plain samples from the experimental test. The predicted curves of two notched laminates highlight that the numerical model could capture both the maximum nominal stress and the softening behavior. As visualized in Figure 3, in the experimental tests data, after maximum stress, the laminate first displays a softening behavior due to the degradation of the Al foil layer and then a plateau (in the case of the plain sample) and another softening branch with a different slope, which is supported seen by only a polymer layer. While in the numerical models of plain and notched laminates, after maximum stress, just one softening branch exists that is due to the modeling of the

Maximum Traction (MPa)
Ultimate Separation Displacement (mm) 17.9×10 6 17.5 36 Figure 3. Comparison of the engineering stress-strain curve of plain and notched laminates obtained by experimental test and numerical methods using interface layer.

The Influence of Surface Roughness
In the current study, as a kind of imperfection in the laminate, surface roughness is considered on the surface of the laminate layer and its impact is studied numerically. Here, we focus on the effect of the surface roughness of Al and the LFPE on the laminate's overall behavior, taking into account the influence of interfacial damage using calibrated interfacial layer properties.
For simplicity, it is assumed that the wavy surface has the following idealized sinusoidal form [25]: where a and λ denote the amplitude and wavelength, with the x-axis pointing to the right horizontally (see Figure 4a). Depending on how the thin Al film is deposited, three possible cases are envisioned to exist: (a) the Al surface is flat, as sketched in Figure 4a; (b) the Al surface is wavy, identical to that of the LDPE surface ( Figure 4b) and (c) the Al surface is wavy, but the LDPE surface is flat (Figure 4c).

The Influence of Surface Roughness
In the current study, as a kind of imperfection in the laminate, surface roughness is considered on the surface of the laminate layer and its impact is studied numerically. Here, we focus on the effect of the surface roughness of Al and the LFPE on the laminate's overall behavior, taking into account the influence of interfacial damage using calibrated interfacial layer properties.
For simplicity, it is assumed that the wavy surface has the following idealized sinusoidal form [24]: y = a sin(π/λ)x, where a and λ denote the amplitude and wavelength, with the x-axis pointing to the right horizontally (see Figure 4a). Depending on how the thin Al film is deposited, three possible cases are envisioned to exist: (a) the Al surface is flat, as sketched in Figure 4a; (b) the Al surface is wavy, identical to that of the LDPE surface ( Figure 4b) and (c) the Al surface is wavy, but the LDPE surface is flat (Figure 4c).  To investigate the effect of surface roughness on the overall behavior of the laminate, three different values for amplitude ratio (a/λ) are considered with three wavelengths (λ) equal to h/2, ℎ and 2h, in which h is the Al foil thickness. Therefore, the ratio of a/λ in the numerical models is changed to 0, 0.025 and 0.05. The behavior of the laminate reveals that upon applying the load, necking occurs in the Al foil layer and the thickness of Al foil plays an important role in the appearance of this necking and final failure mode of the  To investigate the effect of surface roughness on the overall behavior of the laminate, three different values for amplitude ratio (a/λ) are considered with three wavelengths (λ) equal to h/2, h and 2h, in which h is the Al foil thickness. Therefore, the ratio of a/λ in the numerical models is changed to 0, 0.025 and 0.05. The behavior of the laminate reveals that upon applying the load, necking occurs in the Al foil layer and the thickness of Al foil plays an important role in the appearance of this necking and final failure mode of the laminate. In this regard, the wavelength of the roughness was chosen to be proportional to the Al foil thickness. Figure 5a represents the influence of surface roughness on Al and LDPE separately and on both surfaces with different amplitudes and λ equal to h/2. The results highlight that maximum stress occurs for the case of roughness only on the Al foil surface, and roughness only on the LDPE results in the lowest level of stress. Roughness also affects the softening branch such that compared to the without-surface roughness, the laminate degradation happens earlier. Figure 5b depicts the effect of surface roughness with a wavelength equal to the Al foil thickness. The same behavior in the λ equal to h/2 is also observed here, but the rate of decrease in the stress level by changing the roughness from the Al foil surface to the LDPE surface is higher, and it will be magnified by higher amplitudes.
Increasing the wavelength leads to a significant difference in the overall behavior of the laminate with roughness on its layers. The results in Figure 5c indicate that a higher value of wavelength weakens the overall mechanical behavior of the laminate, and strain localization triggers at a much lower strain level. It is worth noting that in spite of the lower level of stress for the case of roughness on the LDPE layer, the anticipation of softening branch happens for the roughness on both Al foil and LDPE layers.
The distribution of stresses parallel to the loading direction (along the x-axis) at maximum nominal stress for different wavelengths is shown in Figure 6. In Figure 6a, there is no roughness on the surface of laminate layers, the stress in the Al foil layer is very high and its distribution is uniform. On the other hand, the stress value is very low at the polymer layer but still with a uniform distribution. The distribution of stresses in the laminate with roughness on Al foil, LDPE and both with the amplitude equal to 0.05 is visualized in Figure 6b-d. Surface roughness causes a change in the distribution of stress in the Al foil layer only, thus reducing stress, particularly at the surface but increasing stress at the edges. In addition, at the Al foil layer, stress concentration occurs at specific locations along the Al foil thickness. Applying higher elongation to the laminate can reveal the stress concentration more clearly by forming cross bands.  To investigate the effect of surface roughness on the overall behavior of the laminate, three different values for amplitude ratio (a/λ) are considered with three wavelengths (λ) equal to h/2, ℎ and 2h, in which h is the Al foil thickness. Therefore, the ratio of a/λ in the numerical models is changed to 0, 0.025 and 0.05. The behavior of the laminate reveals that upon applying the load, necking occurs in the Al foil layer and the thickness of Al foil plays an important role in the appearance of this necking and final failure mode of the laminate. In this regard, the wavelength of the roughness was chosen to be proportional to the Al foil thickness. Figure 5a represents the influence of surface roughness on Al and LDPE separately and on both surfaces with different amplitudes and λ equal to h/2. The results highlight that maximum stress occurs for the case of roughness only on the Al foil surface, and roughness only on the LDPE results in the lowest level of stress. Roughness also affects the softening branch such that compared to the without-surface roughness, the laminate degradation happens earlier. Figure 5b depicts the effect of surface roughness with a wavelength equal to the Al foil thickness. The same behavior in the λ equal to h/2 is also observed here, but the rate of decrease in the stress level by changing the roughness from the Al foil surface to the LDPE surface is higher, and it will be magnified by higher amplitudes.
Increasing the wavelength leads to a significant difference in the overall behavior of the laminate with roughness on its layers. The results in Figure 5c indicate that a higher value of wavelength weakens the overall mechanical behavior of the laminate, and strain localization triggers at a much lower strain level. It is worth noting that in spite of the lower level of stress for the case of roughness on the LDPE layer, the anticipation of softening branch happens for the roughness on both Al foil and LDPE layers.
(a) The distribution of stresses parallel to the loading direction (along the x-axis) at maximum nominal stress for different wavelengths is shown in Figure 6. In Figure 6a, there is no roughness on the surface of laminate layers, the stress in the Al foil layer is very high and its distribution is uniform. On the other hand, the stress value is very low at the polymer layer but still with a uniform distribution. The distribution of stresses in the laminate with roughness on Al foil, LDPE and both with the amplitude equal to 0.05 is visualized in Figure 6b-d. Surface roughness causes a change in the distribution of stress in the Al foil layer only, thus reducing stress, particularly at the surface but increasing stress at the edges. In addition, at the Al foil layer, stress concentration occurs at specific locations along the Al foil thickness. Applying higher elongation to the laminate can reveal the stress concentration more clearly by forming cross bands. The distribution of equivalent plastic strain at maximum stress is visualized in Figure 7 for a perfect laminate and laminate with different wavelengths. When there is roughness on the LDPE surface, the distribution of equivalent plastic strain will not remain constant at the polymer layer compared to the flat Al surface. Hence, in the presence of surface irregularities, the polymer layer accompanies Al foil to redistribute the applied strain. The results also highlight that in the presence of surface roughness, the strain will localize to create shear bands, as represented in Figure 7; by reducing the wavelength, the number of shear bands increases. In general, increasing the wavelength results in a higher amount of plastic strain in the laminate and comparing three different surface roughness conditions demonstrates that the laminate with roughness on both layers experiences severe plastic strain localization.
The results in Figure 8 depict the influence of surface roughness on the failure mode of the laminate at the end of loading with two different wavelengths. As visualized in Figure 8a-d, the failure mode of the laminate with a wavelength of h/2 is almost the same for the case of roughness only on the Al foil surface (Figure 8b), on the LDPE surface ( Figure 8c) and on both Al and LDPE surfaces (Figure 8d) in which there will be one shear band due to necking in the Al foil layer. On the other hand, in the case of no surface roughness, the laminate failure mode changes by increasing the wavelength, resulting in the emergence of necking in the Al foil layer (see Figure 8a,e). As visualized in Figure 8f-h, two neckings will be on both sides of the Al foil due to roughness on the Al surface, one necking in the middle of the Al foil layer with roughness on the LDPE surface (Figure 8g), and one side necking with multi-debonding along the length of the laminate interface when there is roughness on both Al foil and LDPE surfaces (Figure 8h).  The distribution of equivalent plastic strain at maximum stress is visualized in Figure  7 for a perfect laminate and laminate with different wavelengths. When there is roughness on the LDPE surface, the distribution of equivalent plastic strain will not remain constant at the polymer layer compared to the flat Al surface. Hence, in the presence of surface irregularities, the polymer layer accompanies Al foil to redistribute the applied strain. The results also highlight that in the presence of surface roughness, the strain will localize to create shear bands, as represented in Figure 7; by reducing the wavelength, the number of shear bands increases. In general, increasing the wavelength results in a higher amount of plastic strain in the laminate and comparing three different surface roughness conditions demonstrates that the laminate with roughness on both layers experiences severe plastic strain localization.  The results in Figure 8 depict the influence of surface roughness on the failure mode of the laminate at the end of loading with two different wavelengths. As visualized in Figure 8a-d, the failure mode of the laminate with a wavelength of h/2 is almost the same for the case of roughness only on the Al foil surface (Figure 8b), on the LDPE surface (Figure 8c) and on both Al and LDPE surfaces (Figure 8d) in which there will be one shear band due to necking in the Al foil layer. On the other hand, in the case of no surface roughness, the laminate failure mode changes by increasing the wavelength, resulting in the emergence of necking in the Al foil layer (see Figure 8a,e). As visualized in Figure 8f-h, two neckings will be on both sides of the Al foil due to roughness on the Al surface, one necking in the middle of the Al foil layer with roughness on the LDPE surface (Figure 8g), and one side necking with multi-debonding along the length of the laminate interface when there is roughness on both Al foil and LDPE surfaces (Figure 8h).

The Influence of Different Interfacial Properties
This section investigates the influence of interfacial properties without imperfections (surface roughness). To study the effect of interface properties on the laminate behavior, different values for traction and separation are defined as listed in Table 5. Logically, the initial stiffness will change accordingly by changing the traction and separation values. The overall behavior of the material with different traction and separation values is provided in Figure 9. The results emphasize that a change in traction leads to different behavior, particularly in the softening branch, such that the higher traction retards the Al degradation, and by reducing the traction, the maximum stress reduces as well.
By increasing the separation quantities ( /ℎ) from 0.01 to 100, the results highlight that if the separation is low, adopting high quantities for traction leads to the same behavior, while for higher separation, the laminate behavior is almost unaffected by lower tractions.
The primary function of LDPE connected to a thin Al foil is to improve the Al foil's final elongation. In order to attain this goal in a real-world application, the interfacial layer properties should be adjusted to have weak separation and medium traction properties.

The Influence of Different Interfacial Properties
This section investigates the influence of interfacial properties without imperfections (surface roughness). To study the effect of interface properties on the laminate behavior, different values for traction and separation are defined as listed in Table 5. Logically, the initial stiffness will change accordingly by changing the traction and separation values. The overall behavior of the material with different traction and separation values is provided in Figure 9. The results emphasize that a change in traction leads to different behavior, particularly in the softening branch, such that the higher traction retards the Al degradation, and by reducing the traction, the maximum stress reduces as well.

Closing Remarks
This research concerns material systems made of thin metal foils and polymer ply coupling. In the present work, computational approaches were employed to investigate the mechanical behavior of thin laminated structures, considering a kind of imperfection By increasing the separation quantities (δ p /h) from 0.01 to 100, the results highlight that if the separation is low, adopting high quantities for traction leads to the same behavior, while for higher separation, the laminate behavior is almost unaffected by lower tractions.
The primary function of LDPE connected to a thin Al foil is to improve the Al foil's final elongation. In order to attain this goal in a real-world application, the interfacial layer properties should be adjusted to have weak separation and medium traction properties.

Closing Remarks
This research concerns material systems made of thin metal foils and polymer ply coupling. In the present work, computational approaches were employed to investigate the mechanical behavior of thin laminated structures, considering a kind of imperfection and the interfacial layer effect. First, the interfacial layer properties were calibrated based on available uniaxial experimental tests on thin laminates, with and without crack. Calibrated parameters could reflect the nominal stress-strain behavior of plain and notched laminates closely comparable to experimental behavior. They could capture both maximum stress and the softening behavior after the degradation of the thin Al foil.
Then, a wavy surface with a sinusoidal shape was defined on both the Al foil and polymer surfaces to investigate its influence on the overall behavior of the thin laminate and its failure mode. A numerical model consisting of Al foil coupled to the polymer under plane strain conditions was developed, and an interface layer between Al and polymer with zero-thickness cohesive elements characterized by T-S law (using calibrated parameters) was defined. The results based on the developed numerical models can be highlighted as:

•
Among three cases (roughness on Al, LDPE and both), the wavy surface on only the LDPE layer led to the most critical condition such that the maximum level of stress decreased significantly when increasing the amplitude of the wavy interface.

•
The failure mode of the laminate changes by varying the imperfection location and different failure modes observed by defining surface roughness on each and both layers.

•
As the existence of roughness on the surface of each layer or both in the laminate is inevitable during manufacturing, the results of this study indicate the significant influence of this type of imperfection on the functionality of thin metal/polymer laminates, which is worthwhile in real-world applications.
As the interface layer properties play an important role in the laminate's overall behavior, a parametric study was also performed, changing the traction and separation quantities in the T-S law. The results indicated that: • At a constant separation, increasing tractions leads to higher stress and retarding the strain localization. On the other hand, at low separation, adopting higher traction results in almost the same mechanical behavior, and at high separation, the behavior of the laminate is unchanged by reducing the traction.

•
Regarding the mentioned applications, these results can emphasize that adopting very high traction or very low separation does not necessarily lead to a strong interface behavior, but an optimum combination should be adopted.