## Abstract

This study developed microstructure-based finite element (FE) models to investigate the behavior of cold-sprayed aluminum–alumina (Al-Al_{2}O_{3}) metal matrix composite (MMC) coatings subject to indentation and quasi-static compression loading. Based on microstructural features (i.e., particle weight fraction, particle size, and porosity) of the MMC coatings, 3D representative volume elements (RVEs) were generated by using Digimat software and then imported into ABAQUS/Explicit. State-of-the-art physics-based modeling approaches were incorporated into the model to account for particle cracking, interface debonding, and ductile failure of the matrix. This allowed for analysis and informing on the deformation and failure responses. The model was validated with experimental results for cold-sprayed Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} metal matrix composite coatings under quasi-static compression by comparing the stress versus strain histories and observed failure mechanisms (e.g., matrix ductile failure). The results showed that the computational framework is able to capture the response of this cold-sprayed material system under compression and indentation, both qualitatively and quantitatively. The outcomes of this work have implications for extending the model to materials design and for applications involving different types of loading in real-world application (e.g., erosion and fatigue).

## Introduction

Aluminum (Al)-based metal matrix composites (Al-MMCs) have been widely used to enhance surface performance (Ref 1,2,3,4) due to their superior mechanical properties such as high stiffness and strength (Ref 5), low coefficient of thermal expansion (Ref 6), improved wear (Ref 7) and corrosion resistance (Ref 8), and better creep–fatigue performance (Ref 9). To fabricate Al-MMCs, a wide range of reinforcing ceramic particles are added to an Al matrix (e.g., SiC (Ref 10), Al_{2}O_{3} (Ref 1), TiC (Ref 11), TiB_{2} (Ref 12), and B_{4}C (Ref 13)). Among these reinforcing particles, Al_{2}O_{3} has been frequently used due to improved corrosion resistance and chemical stability (Ref 14). Recently, the Al-Al_{2}O_{3} composites fabricated by the cold spray additive manufacturing process have been used as protective coatings against wear, erosion, corrosion, and high temperature degradation in aerospace and other industrial sectors (Ref 15, 16). Among many deposition routes for producing coatings made of pure Al, cold spray stands out due to minimal heat loading of the substrate during the deposition process (Ref 17). In addition, this additive manufacturing method provides the possibility of producing multi-phase coatings via mixing feedstock powders by which the hard phases such as SiC and Al_{2}O_{3} can be incorporated effectively in the Al-MMCs coatings (Ref 18). These hard particles play a critical role in lowering the wear rate of the ceramic–metal coatings (Ref 19).

In research studies on experimental mechanics of ceramic–metal composite coatings, the behavior of cold-sprayed MMC coatings has been extensively addressed in terms of dry sliding wear (Ref 20), indentation (Ref 15), flexural properties (Ref 21), and erosion (Ref 16). However, there are a limited number of articles that investigate the response of MMC coatings through computational approaches, which is addressed in the current study. For example, Bolelli et al. (Ref 22) generated 2D RVEs based on the SEM images of WC–CoCr and WC–FeCrAl coatings to simulate ball-on-disk test and compression. The matrix and hard particles were modeled as elastic–plastic and pure elastic materials, respectively. The numerical results of Bolelli et al. (Ref 22) were compared with experiments based on the morphology of the worn surfaces and Young’s modulus. In a separate study by Balokhonov et al. (Ref 23), 2D models at micro-, meso-, and macroscales have been implemented by using the finite difference method for MMC coatings under tension and compression. It was found that curvilinear interfaces lead to stress concentration giving rise to the formation of shear bands in the Al matrix locally as well as cracking in the ceramic particles.

From a computational perspective, the majority of previous studies on Al-MMCs have focused on axial tensile loading, by which the microstructure has been modeled via 2D (Ref 24) and 3D (Ref 25) RVEs. In a large number of studies, the occurrence of the three competitive damage mechanisms in particulate-reinforced Al matrix composites (PRAMCs) subject to tension, such as matrix ductile failure, matrix/particle debonding, and particle cracking, has been mainly explored by using phenomenological ductile failure criteria (Ref 25), cohesive zone models (Ref 26), and a conventional brittle cracking model (Ref 13) (i.e., elastic cracking behavior which employs the Rankine criterion for failure initiation (Ref 48)), respectively. For example, Zhang et al. (Ref 25) investigated the behavior of a 7vol.% SiCp/Al composite made by a stir casting technique and incorporated the three damage mechanisms in a real microstructure-based 3D RVE. The numerical results revealed that particle fracture and interfacial debonding emerged as the initial failure mechanisms in the composite under tensile loading.

In contrast to the numerous numerical studies that explore the tensile behavior of Al-MMCs, limited efforts have been made to address the indentation and compressive behavior computationally, particularly with emphasis on damage mechanisms. For example, Park et al. (Ref 27) investigated Al-SiC MMCs under quasi-static compression up to 1% strain using a statistical synthetic RVE made by the DREAM.3D software—no damage mechanisms were incorporated into the model. The indentation behavior of MMCs has been frequently explored via 2D models (Ref 28). More recently, Shedbale et al. (Ref 29) employed homogeneous and heterogeneous 3D FE models to study the ball indentation response of particulate-reinforced MMCs. The results showed that the heterogeneous model tends to overestimate the hardness compared to the experiments due to the local concentration of particles under the indenter.

Motivated by previous numerical studies, which primarily focus on the tensile response of MMCs, this work aims to investigate the compressive and indentation behavior of Al-Al_{2}O_{3} coatings by using 3D RVEs produced by Digimat software based on microstructural characteristics obtained using scanning electron microscope images. To account for the damage mechanisms, the Gurson–Tvergaard–Needleman (GTN) model (Ref 30) was applied to the Al matrix and the matrix/particle debonding was modeled by the CZM method (Ref 31). For the ceramic particle phase, the Johnson–Holmquist II (JH2) model (Ref 32) was used to incorporate particle cracking into the FE model. The model was validated with the experimental data for Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} coatings in terms of stress–strain histories, failure mechanisms, and Vickers hardness. The results show that the model has the potential to be employed for future parametric studies for material design and optimization, which tailors concentration of reinforcing particles to balance strength and density for weight-sensitive applications.

## Experimental Procedures

### Material and Specimen Preparation

In this study, pure Al and Al-Al_{2}O_{3} composite coatings were fabricated using a low-pressure cold spray system (SST series P, CenterLine Ltd., Windsor, ON, Canada), as shown schematically in Fig. 1(a), which was connected to a volumetric powder feeder (5MPE, Oerlikon Metco, Westbury, NY, USA). Based on previous studies (Ref 2, 33, 34), the air temperature and pressure were set to 375 °C and 620 kPa, respectively. The nozzle was manipulated by a robot (Motoman HP-20, Yaskawa Electric Corp., Waukegan, IL, USA). The cold spray nozzle traversed across the Al substrate at a speed of 15 mm/s to transfer the feedstock powder to the substrate, and the deposition process produced five layers of the coating. The feedstock powder blend was developed through a three-step process: gas atomization, sieving, and mixing (see Fig. 1a). Aluminum (99.0%) powder (CenterLine Ltd., Windsor, ON, Canada) and α-Al_{2}O_{3} with a purity of 99.5% (Amdry 6060, Oerlikon Metco Inc., Westbury, NY, USA) were used in this study. The Al and Al_{2}O_{3} powders were sieved to obtain a size distribution of 40–60 and 30–45 *µ*m, respectively. The Al and Al_{2}O_{3} powders were admixed to produce powder blends containing 0, 60, and 90 wt.% Al_{2}O_{3}. This process was conducted using a cylinder with a 20 mm diameter whose angular velocity and operating time were set to 20 RPM and 30 minutes, respectively. As shown by Shao et al. (Ref 33), when deposited into coatings, the mixed powder blend with 60 wt.% Al_{2}O_{3} produced coatings that were Al-34 wt.% Al_{2}O_{3} and the mixed powder blend with 90 wt.% Al_{2}O_{3} produced coatings that were Al-46 wt.% Al_{2}O_{3}.

### Mechanical Testing and Characterization

As shown in Fig. 1(b), the cold-sprayed coating deposits were cut into cuboidal specimens with dimensions of 2.3 mm in length, 2.7 mm in width, and 3.5 mm in height using wire electrical discharge machining. The samples were used for quasi-static compression testing, where the loading was applied in the direction of the 3.5 mm dimension. The experiments were conducted using the displacement control technique up to a maximum displacement of 1 mm at a constant rate of 1 × 10^{−3} s^{−1} using an Instron 3365 testing machine (Instron, Norwood, Massachusetts, USA). To visualize the features of the macroscopic deformation of the specimen surface, the machine was equipped with an AOS PROMON U750 high-speed camera with a full resolution of 1280 × 1024 and a VIC 900170WOF LED laser light guide for illumination. This assembly was coupled with digital image correlation (DIC) capabilities using VIC-2D (v6 2018) software (Ref 35) (Correlated Solutions Irmo, SC, USA) to monitor the strain fields, which is detailed in Shao et al. (Ref 33). The specimen was aligned between the loading platens made from M2-graded high-speed steel with a diameter of 1 inch (see Fig. 1b), and extreme pressure grease was applied on the interfaces to eliminate the effect of friction and allow free lateral expansion. The compression tests were carried out as per ASTM Standard C1424-15 (Ref 36) at room temperature and repeated four times for each coating with different reinforcing particle content. To inform the microstructure-based models as related to the reinforcing particle content, porosity, particle shape, size, and distribution, microstructural characterization was done using a field-emission SEM coupled with energy-dispersive x-ray spectroscopy (EDS) operated at 20 kV (Zeiss Sigma, Oberkochen, Baden-Württemberg, Germany), as shown in Fig. 1(c) for the Al-46 wt.% Al_{2}O_{3} composite. The porosity of the samples was estimated using ImagePro software coupled with the SEM images, and the porosity was found to be 2.84 ± 0.31 vol.% in pure Al, 0.23 ± 0.04 vol.% in Al-34 wt.% Al_{2}O_{3}, and 0.17 ± 0.03 vol.% in Al-46 wt.% Al_{2}O_{3}. In addition, the EDS analysis revealed that the feedstock powders with 60 and 90 wt.% Al_{2}O_{3} led to depositions with 34 ± 2.56 and 46 ± 2.04 wt.% of ceramic particles, respectively.

## Numerical Methodology

A 3D FE model based on the coating microstructural features is presented to explore the behavior of Al-Al_{2}O_{3} MMC coatings under quasi-static compressive loading. The 3D representative volume elements (RVEs) were generated by the Digimat software for Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} MMCs. The RVEs were imported into the ABAQUS/Explicit solver (release 6.14). For the micro-indentation test, the homogenization approach (Ref 29) was applied and experimental compression data were used to extract the effective mechanical properties for each MMC coating sample.

### Geometry and Model Description

RVEs with different sizes have been considered in previous studies (Ref 13, 37). For example, Ma et al. (Ref 37) found no significant difference in the tensile stress–strain responses by varying the RVE size from 20 µm to 50 µm. In this study, an RVE length of 100 µm was chosen based on the microstructural features (e.g., average reinforcing particle size of 15 µm). The SEM images (see Fig. 1c and d) were first used to extract the distribution of particle size in the composites (e.g., the size of the alumina particles range from 1 to 30 µm in the Al-46 wt. % Al_{2}O_{3} composites (see Fig. 1c). Next, the measured range of distribution was incorporated into the RVE using a uniform distribution in the Digimat software, which accounts for the variation and uncertainty in the particle size distribution that is likely to be observed through the SEM images from different locations of the sample. Additionally, the reinforcing ceramic particles represent an irregular shape in the SEM micrographs. Here, an icosahedron geometry was used to account for the shape of the particles, which has also been used to represent the irregular-shaped particles in previous studies involving MMC coatings (Ref 37,38,39). The time step was set at 2 µs, which was found to meet the quasi-static loading condition (Ref 40). As shown in Fig. 2(a), the top and bottom boundary faces have been fully restricted to the reference points (RPs) by the kinematic coupling constraints to facilitate application of load/boundary conditions and obtain the stress–strain response of the RVE. A corner of the RVE was fully restricted to prevent the material from rigid boundary motion. The degrees of freedom of the bottom boundary face were fixed within all degrees of freedom, except for the in-plane displacements (i.e., Y- and Z-directions in Fig. 2b). The same boundary conditions for application of compressive loading were also used in other studies (Ref 41, 42). All of the constituents were discretized by 3D linear tetrahedral C3D4 elements. Following a mesh quality assessment to decrease the likelihood of element distortion at high strains, an average element size of 1.5 *μ*m was used. On this basis, the RVE for Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} MMCs was meshed by 741,222 and 1,275,358 elements, respectively. Micro-indentation Vickers testing of the composite coatings was simulated via the homogenization approach (Ref 29). The effective mechanical properties, including the Young’s modulus and flow stress, were extracted from the experimental compressive stress–strain histories for each particle concentration in the coatings and were an input into the approach. Figure 3 shows the FE model of the Vickers test. Due to symmetry, only one quarter of the homogenized block with symmetric constraints was modeled. The size of the block was determined as per the work of Shedbale et al. (Ref 29) to achieve convergence in the indentation response.

As shown by Shedbale et al. (Ref 29), the Vickers indenter was considered as a discrete rigid body and fully confined, except for the vertical direction Y (see Fig. 3). The bottom surface of the block was fixed, and the lateral surfaces were free to deform. A frictionless contact was defined between the indenter and the top surface of the block, which was implemented by the standard surface-to-surface contact algorithm (Ref 29). By using an average element size of 5 µm, the block and the indenter were meshed by using 30,276 first-order eight-node 3D elements with reduced integration (C3D8R) and 2104 quadrilateral four-node 3D rigid elements (R3D4), respectively.

### Material Models

In this study, the micromechanical-based Gurson–Tvergaard–Needleman (GTN) model (Ref 30) was applied to capture the matrix failure. Given experimental evidence of particle cracking and interfacial debonding under compression, the Johnson–Holmquist II (JH2) model (Ref 32) for the failure of ceramic particles and the cohesive zone model (CZM) approach (Ref 31) for the matrix/particle debonding failure were incorporated into the model.

#### The GTN Model for Matrix Failure

Ductile failure of metals occurs as a result of a three-stage mechanism of nucleation, growth, and coalescence of voids. Voids are nucleated from an inclusion or as a consequence of either cracking or decohesion of second-phase particles and then grow due to the localization of plastic strain under a high stress triaxiality state (Ref 43). To analyze this phenomenon theoretically, the porous metal plasticity model proposed by Tvergaard and Needleman (i.e., the GTN model) (Ref 30), as a modified version of the Gurson model (Ref 44), has been widely used (Ref 45,46,47). The yield function of the GTN model is expressed as follows:

where \(\emptyset\) denotes the non-dilatational strain energy and \(q_{1}\), \(q_{2}\), and \(q_{3}\) are the constants proposed by Tvergaard (Ref 48) to account for the effects induced by void interaction due to multiple-void arrays and to provide better consistency with experimental data. Here, \(\sigma_{q}\) and \(\sigma_{y}\) represent von Mises stress and the flow stress of the undamaged material. To model the rapid deterioration of stress carrying capacity caused by void coalescence, the parameter, \(f^{*}\), known as the effective porosity was first introduced by Tvergaard–Needleman (Ref 30). The function is specified as follows:

where \(f_{c}\) is the critical void volume fraction (VVF) at the onset of the coalescence, \(f_{u}^{*} = \frac{1}{{q_{1} }}\) corresponds to zero stress carrying capacity, and \(f_{F}\) denotes the VVF when the material has completely failed, which governs the element deletion process. The increase in the VVF is deemed as the summation of the increment owing to void nucleation and the growth of existing voids. The function can be written as

Assuming plastically incompressible behavior for the material, the void growth rate (i.e., \({\text{d}}f_{{{\text{growth}}}}\)) can be expressed as a function of the plastic volume change as follows:

where \({\text{d}}\varepsilon_{ii}^{p}\) denotes the trace of the plastic strain rate tensor. The nucleation of voids is considered to be exclusively dependent on the plastic strain, and it was assumed that occurrence of void nucleation occurred only under hydrostatic tension (Ref 49, 50). On this basis, the function is written as

where \(p\) represents the hydrostatic stress, \(f_{N}\) is the void volume fraction of the nucleated void, \(\varepsilon_{N}\) is the mean equivalent plastic strain for void nucleation, and \(s_{N}\) is the standard deviation of the distribution. Here, the rate of equivalent plastic strain,\(d\overline{\varepsilon }_{ }^{p}\), is obtained by enforcing equality between the matrix plastic dissipation and the rate of macroscopic plastic work as follows:

Table 1 summarizes the constants of the GTN model used in this study. The proposed values by Tvergaard (Ref 48) were utilized for the \(q_{1}\), \(q_{2}\), and \(q_{3}\) constants. From SEM characterization, the initial porosity, \(f_{0}\), was assumed to be an average quantity of 0.0017 and 0.0023 for Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} composites, respectively. The parameters \(\varepsilon_{N}\), \(s_{N}\), and \(f_{N}\) were obtained from previous studies for pure Al (Ref 39, 51).

#### The JH2 Model for Particle Cracking

For damage modeling of ceramics, the phenomenological Johnson–Holmquist models have been commonly used (Ref 38, 52,53,54) to depict the behavior of these materials, including pressure-dependent strength, strain rate dependency, and dilation or bulking effects (Ref 32). The strength and damage are expressed as analytical functions of pressure and other parameters as (Ref 32)

where \(\sigma^{*}\) and \(\sigma_{i}^{*}\) represent the normalized intact equivalent stress, \(\sigma_{f}^{*}\) denotes the normalized fracture stress, and *D* is the damage variable, varying from 0 to 1. Here, \(\sigma^{*} = \sigma /\sigma_{{{\text{HEL}}}}\), \(P^{*} = P/P_{{{\text{HEL}}}}\), \(T^{*} = T/P_{{{\text{HEL}}}}\), and \(\dot{\varepsilon }^{*} = \dot{\varepsilon }/\dot{\varepsilon }_{0}\), where \(\sigma\) is the actual equivalent stress, \(\sigma_{{{\text{HEL}}}}\) is the equivalent stress at the Hugoniot elastic limit (HEL), *P* is the actual pressure, *P*_{HEL} is the pressure at the HEL, *T* is the maximum tensile hydrostatic pressure tolerated by the material, \(\dot{\varepsilon }\) is the actual strain rate, and \(\dot{\varepsilon }_{0} = 1\) is the reference strain rate. *A*, *B*, *C*, *M*, and *N* are material constants which need to be calibrated for each material. The maximum value of \(\sigma_{f}^{*}\) is defined by SFMAX (i.e., the maximum limitation of the normalized fractured strength). Once the yield function is met as per Eq 11, the damage begins to accumulate based on the incremental equivalent plastic strain defined as

where \(\sigma_{q}\), *D*_{1}, and *D*_{2} are deemed as material constants. To calculate pressure, *P*, a polynomial equation of state (EOS) is employed, which is defined as

where \(K_{1}\) denotes the bulk modulus, \(K_{2}\) and \(K_{3}\) are material constants, \(\eta\) is the specific volume, *U* represents the internal energy which is related to the equivalent plastic flow stress \(\sigma_{y}\) by a quadratic expression, \(\beta\) is the fraction of the elastic energy loss converted to potential hydrostatic energy, and the shear modulus is shown by *G*. The 21 constants of the JH2 model for Al_{2}O_{3} were obtained from previous studies (Ref 54, 55) and are summarized in Table 2.

#### The Cohesive Zone Model for the Interfacial Debonding

The cohesive zone model (CZM) was first proposed by Barrenblatt (Ref 31) and Dugdale (Ref 56) and is now widely used as an effective approach for modeling the fracture process in materials such as polymers, metals, ceramics, concretes, and laminated composites (Ref 57). It was assumed that matrix/particle interface debonding occurs when the quadratic interaction function involving the nominal stress ratios attains unity as follows (Ref 58):

where \(T_{n}\), \(T_{s}\), and \(T_{t}\) represent the tractions acting on the interface at a load increment in normal and in two in-plane shear directions, respectively. Likewise, \(T_{n}^{0}\), \(T_{s}^{0}\), and \(T_{t}^{0}\) denote the tractions at the onset of damage initiation in normal and in two in-plane shear directions. Note that the normal traction is tensile and pure compressive stress does not lead to decohesion. The components of the traction–separation law are written as follows:

where \(T_{n}^{*}\), \(T_{s}^{*}\), and \(T_{t}^{*}\) are the stress components calculated by the elastic traction–separation behavior for the current strain prior to the damage initiation. *D* denotes the damage variable which begins to gradually increase from 0 to 1 with further loading once the debonding initiation criterion expressed by Eq 18 is met. The damage variable is defined as (Ref 58)

where \(\delta_{m}^{0}\) and \(\delta_{m}^{f}\) represent the effective separations at damage initiation and complete failure, respectively. The maximum value of the effective displacement during the loading process at each increment is shown by \(\delta_{m}^{\max }\). The effective separation at each load increment \(\delta_{m}\) is calculated as (Ref 58)

where \(\delta_{n}\), \(\delta_{s}\), and \(\delta_{t}\) are the nominal separations in normal and in two in-plane shear directions, respectively. To obtain the effective separation at complete decohesion, one can use the fracture energy \(G_{c}\), which is given as

The CZM constants used in this study are presented in Table 3. For Al-MMC composites, the CZM constants have been reported in different ranges in the previous studies (Ref 13, 59,60,61,62). For example, the reported values for the interface strength vary from a quantity on the order of MPa (Ref 61) to GPa (Ref 62). The values in Table 3 were selected to establish the best match between the experimental and numerical outcomes in this research study.

## Results

### Quasi-Static Compression

The numerically predicted results for Al-34 wt.% Al_{2}O_{3} and Al-46 wt.% Al_{2}O_{3} coatings under quasi-static compression are compared with those of the experiments in terms of the stress versus strain histories and observed failure mechanisms. This can provide insights for establishing an accurate computational framework to further explore the behavior of the material that can eventually give rise to a tool for material design and optimization. Figure 4 shows the predicted stress–strain responses in comparison with those measured by experiments. The pure Al and MMC samples were all experimentally tested in different directions, namely the nozzle traverse (travel) direction, the second in-plane direction perpendicular to the nozzle travel direction, and the deposition direction represented by *X*, *Y*, and *Z*, respectively, in Fig. 4.

For the pure Al matrix, the results based on the data for the *Z*-direction were compared to experiments, as shown in the red solid curve in Fig. 4. The predicted curve for the pure Al matrix aligns with the experimentally measured one, which shows the accuracy of the approach for modeling the pure Al coating. Regarding the predicted responses for the MMCs, the model can reasonably capture the stiffness (i.e., the Young’s modulus) and the maximum load bearing capacity of the composite coatings with different reinforcing particle concentrations. In addition, the experimental trend toward decreasing ductility with an increase in reinforcing particle content from 34 to 46 wt.% is reasonably reflected in the numerically predicted curves. Namely, plastic deformation in the Al-34 wt.% Al_{2}O_{3} composite coating begins to take place in the model at a strain of 0.6%, which leads to a yield stress of 141 MPa. The experimental yield strain ranges from 0.53% to 0.65% and the yield stress ranges from 135 MPa to 178 MPa based on the different coordinate directions that were explored. The model tends to slightly underestimate the yield stress of the MMC with 34 wt.% of the reinforcing inclusions. This is likely due to cold working, which introduces a hardening and strengthening effect in the matrix during the cold spray process. This is a result of the high-velocity hard ceramic particles impacting on the Al grains (Ref 18, 63), which is not yet considered in the model. In addition to cold working, there are also other mechanisms that are likely to cause the deviations observed between the predicted and experimental stress versus strain histories as it relates to strain hardening behavior and elongation at failure. Among these mechanisms, grain size effects (Ref 64, 65) (e.g., Hall–Petch relationship (Ref 66)), which require implementation of strain gradient plasticity models (Ref 67) in the model of this study, may also play a role.

Figure 5 illustrates the contour of the plastic deformation (Ref 68), interfacial failure (Ref 26), and particle cracking (Ref 24) at two different axial strains during the loading process for the Al-34 wt.% Al_{2}O_{3} composite. With the further application of load, the plastic strain begins to accumulate in the ligament between the particles closely aligned together (see Fig. 5a), which has been reported in previous studies of 2D RVEs (Ref 68). The localization of plastic strain leads to the growth of the void volume fraction in the matrix resulting in a gradual decrease in matrix flow stress. (See the solid black curve in Fig. 4.) In addition, as the strain exceeds 4%, matrix/particle debonding and a particle cracking failure mechanism manifest and develop with the increase in applied load within the RVE (see Fig. 5b and c), particularly at the sharp corner and concaves of the particles (Ref 68). This behavior can be attributed to the mismatch between the mechanical constants as well as the stress concentration at the curvilinear interfaces (Ref 23). Consequently, the flow stress remains almost constant from 5% strain to a strain of 15% and then rises slightly with a further increase in load up to the end of the loading cycle. This behavior is also reflected in the experimental curves, as matrix failure happens locally and does not lead to fracture of the sample (Fig. 4—the dashed curves for the 34 wt.% Al_{2}O_{3} MMCs). From simulation, it can be implied that the spatial distribution of the particles affected by the number of particles and the mean free path parameter (i.e., inter-particle spacing of the reinforcing phase) (Ref 69)—which is determined by the weight percentage of the inclusions—does not lead to fracture of the Al-34 wt.% Al_{2}O_{3} MMC sample, since the number of thin ligaments is not critical enough to develop a 45° shear band (Ref 13, 68), fracturing the specimen. This numerical implication (i.e., the microcracks in the matrix do not evolve to fracture the sample) is also corroborated by the experimental observation of the deformed sample at the end of the loading and the corresponding SEM micrograph shown in Fig. 5(d). As shown, damage mechanisms such as ductile matrix failure and interfacial decohesion emerge locally, leading to the formation of dispersed microcracks that do not coalesce to cause failure in the sample at macro-length scales. For the Al-46 wt.% Al_{2}O_{3} composite, the predicted yield strain and yield strength are 0.7% and 253 MPa, respectively, while the measured yield strain varies between 0.79% and 0.83%, and the yield strength varies between 298 MPa and 317 MPa for the different coordinate directions. The larger discrepancy in the predicted yield stress of the Al-46 wt.% Al_{2}O_{3} MMC compared to that of the Al-34 wt.% Al_{2}O_{3} MMC can be another indication of the importance of including the cold spray-induced hardening and strengthening effect in the model in order to produce more accurate numerical results, especially for MMCs with a high percentage of particle reinforcement. Additionally, the earlier onset of the debonding failure mechanism in the model compared to the experiments can also play a role in the loss of stiffness prior to achieving peak load.

Figure 6 and 7 shows a comparison between the predicted and experimentally observed failure mechanisms in the Al-46 wt.% Al_{2}O_{3} composite. As shown, the FE model can capture the occurrence of the three competitive damage mechanisms (i.e., matrix ductile failure, matrix/particle debonding, and particle cracking) of metal–ceramic coatings under compressive loading. In comparison with the MMC coating with 34 wt.% alumina reinforcement, the plastic strain is severely localized in the thin ligaments between the particles (see Fig. 6a) in the model. The model predicts the formation of 45° shear cracks passing through the Al matrix between the particles as observed in SEM images, as shown in Fig. 7. Once the strain exceeds 4%, particle cracking and matrix/particle debonding are initiated at the sharp corners of the particles (see Fig 6b and c)—which have been experimentally observed, as shown in Fig. 7—and then propagate within the RVE. These failure mechanisms are accompanied by matrix failure due to void growth in the thin ligaments (see Fig. 6a). Once the strain exceeds 5%, the stress bearing capacity starts to decrease slightly and then remains constant up to a strain of 12% (Fig. 4—the solid blue curve). This is a consequence of the development of the damage modes. The elements of the Al matrix in which the porosity has exceeded the critical value are removed from the mesh, leading to an abrupt decrease in load sustaining capacity (see Fig. 4—the solid blue curve at a strain of ~12%). This behavior is in agreement with the experimental trend that the material’s load sustaining ability decreases after a given strain is reached. Overall, the reasonable agreement between the numerical and experimental findings in terms of stress–strain behaviors and failure mechanisms reveals the applicability of the model to conduct parametric studies that translate into tailoring particle and concentration size to control competition between failure mechanisms toward improving strength–density trade-offs.

### Vickers Micro-indentation

The numerical outcomes of the homogenization approach were validated by Vickers hardness experiments. In the experiments, Vickers micro-indentation was applied to the samples with a load of up to 10 N as per ASTM Standard E384 (Ref 70). Figure 8(a) shows the plastically deformed area after complete unloading for the homogenized model of each MMC composite. From the figure, the diagonal length of the indented area decreased as the particle concentration increased, which results in a higher hardness. In other words, at a critical load, the composite material with a higher particle content is less deformed due to the enhanced stiffness and flow stress induced by the hard alumina particles.

Figure 8(b) illustrates the residual von Mises stress distribution for a load of 10 N after unloading in the homogenized model of Al-34 wt.% Al_{2}O_{3} MMC. The distribution pattern for different particle concentrations is the same, and the magnitude of the residual stress rises with an increase in the percentage of alumina particles. From Fig. 8(b), the residual stress follows a continuous distribution as the model does not explicitly account for the microstructure, while the observations reported after using 2D heterogeneous models (Ref 28, 71) showed that the residual stress is localized between the particles. This implies that 3D microstructure-based models are needed to study how the residual stress is developed in the material more realistically. The Vickers hardness of the numerical data was computed using (Ref 72):

where *F* represents the applied load and *d* is the diagonal length of the indented area. Figure 9 shows the numerically predicted Vickers hardness results in comparison with the experimental ones measured in three different directions, namely the nozzle travel direction, the deposition direction, and the third perpendicular direction represented by *X*, *Z*, and *Y*, respectively. As shown, the homogenized model predicts the hardness of MMCs with an error of 10 ± 3% compared to the measured values. The heterogeneous modeling approach was reported to overestimate the hardness of MMCs in previous studies (Ref 28, 29, 71, 73) due to particle consolidation (Ref 74), which required calibration of a significant number of parameters to reach acceptable agreement with the experiments, when compared to the homogenization approach. However, a 3D microstructure-based model allows flexibility to investigate the effect of the size, randomness, and morphology of the particles on the indentation behavior and hardness of cold-sprayed MMC coatings.

## Discussion

This study is the first of its kind to develop 3D finite element models to explore the compressive behavior and hardness of Al-Al_{2}O_{3} composite coatings fabricated by low-pressure cold spray. The presented microstructure-based models for the MMC coatings were built on previous studies of 3D modeling of tensile behavior (Ref 10, 11, 13, 37, 75,76,77,78) and 2D modeling of the compressive response (Ref 54, 68) of particulate-reinforced metal matrix composites. Our model culminated in an acceptable consistency between the numerical and experimental findings, both quantitatively (i.e., based on yield strain, yield strength, and stiffness) and qualitatively (i.e., based on initiation and propagation of damage mechanisms as well as the trends of stress versus strain histories), which lays the foundation to fill the gap in our computational knowledge of ceramic–metal composite coatings. The extended 3D models of Al-Al_{2}O_{3} MMC coatings were generated by Digimat software, where microstructural characteristics obtained by SEM micrographs and EDS analysis (Ref 33) were incorporated into the model, including the distribution of size, shape, and weight fraction (i.e., ~ 34wt. % and ~ 46wt. %) of reinforcing ceramic particles. Once the RVEs were generated, the models were validated using the experimental outcomes in terms of quantitative (i.e., stress versus strain histories) and qualitative (i.e., failure mechanisms) comparisons. This method of validation has also been used in previous studies (Ref 41, 42, 59, 79). The experimental measurements linked with the DIC technique showed that the compressive strength of the material was between 135 and 178 MPa for the composite with 34wt. % of alumina and between 298 and 317 MPa for the composite with 46wt. % of alumina. These are among the highest values reported in the literature (Ref 80,81,82,83) due to the fabrication strategy in this research that employs both matrix strengthening and dispersion strengthening mechanisms (Ref 33). In addition, the measured stiffness of our coatings varied from 28 GPa to 62 GPa for the Al-34 wt. % Al_{2}O_{3} composites and from 48 GPa to 63 GPa for the Al-46 wt.% Al_{2}O_{3} composites. The developed model predicted the strength of the material with an error based on the average experimental quantities of 9.6% and 17.6% for the 34wt. % alumina and the 46wt. % alumina samples, respectively. As well, the numerical prediction of the stiffness yielded a value of 38 GPa and 53 GPa for the 34wt. % alumina and the 46wt. % alumina MMCs, respectively, resulting in an error of 15.5 % and 4.5% with respect to the average measured values. With the maximum error of 17.6%, the model is reasonably in quantitative agreement with the experimental outcomes. Qualitatively, for the first time to the best of our knowledge, the manifestation and evolution of experimentally observed failure mechanisms in ceramic–metal coatings under compression (i.e., matrix ductile failure, matrix/particle debonding, and particle cracking) were all numerically captured through a 3D micromechanical finite element framework. This improves upon previous studies using 2D models (Ref 23, 54, 68) or single-particle 3D models (Ref 11). Additionally, the necessity of developing 3D models in this study was illustrated by Böhm et al. (Ref 84) and Soppa et al. (Ref 85) to adequately capture the plastic strain distribution in two-phase materials.

Finally, the deviations between the simulation and experimental results can stem from the differences between the assumed and real boundary conditions (Ref 86), the complexities of the real microstructures such as reinforcing particle clustering (Ref 79) that are not yet incorporated into the generated RVEs, the fracture of reinforcing particles during the cold spray deposition process (Ref 18) leading to damage accumulation, and the increase in porosity due to interface decohesion and particle cracking (Ref 87), which are not considered in the GTN model implemented in the present study. In addition, the work hardening effect (Ref 88) induced in the Al matrix by the high-velocity impact of hard ceramic particles gives rise to a significant increase in dislocation density (Ref 89), which results in higher strength and hardness in the experimental samples when compared to the model material systems. To account for these crystallographic orientation effects (Ref 90) in the model, electron backscatter diffraction (EBSD)-based RVEs (Ref 91, 92) can be employed in the future as a promising future direction to numerically explore the micromechanical behavior of MMC coatings. Altogether, the presented model established a reasonable match between the predicted and measured outcomes which allows for further exploration of the microstructure–property relationships of the material (i.e., the effect of matrix porosity (Ref 93) as well as particle size (Ref 2, 94), shape (Ref 95), and distribution (Ref 96) of particles on the macroscale behavior). This paves the way to create a 3D computational tool for the design and optimization of ceramic–metal cold-sprayed composite coatings via tailoring the microstructure.

## Conclusion

This study explored the behavior of Al-Al_{2}O_{3} composite coatings under quasi-static compression and indentation loading via FE analysis, both quantitatively (i.e., stress versus strain response) and qualitatively (i.e., the manifestation of damage mechanisms, including matrix ductile failure, interfacial debonding, and particle cracking). For the FE models, 3D RVEs were generated by Digimat software based on the microstructural features of the MMC coating samples with different particle concentrations, and the homogenization approach was employed for modeling the Vickers micro-indentation test. To account for the matrix ductile failure and the matrix/particle decohesion, the GTN model and the CZM approach were employed, respectively. The ceramic particles were modeled using the phenomenological JH2 model to incorporate particle damage accumulation. The FE model was validated by stress–strain histories, Vickers hardness, and damage mechanisms obtained experimentally, and a reasonable agreement was observed between the results. Altogether, the outcomes of this study confirm the applicability of the model to be used as a computational tool for spatially tailoring matrix and particle properties and geometries to develop high-performing gradient coating structures.

## References

- 1.
K.J. Hodder, H. Izadi, A.G. McDonald and A.P. Gerlich, Fabrication of Aluminum-Alumina Metal Matrix Composites Via Cold Gas Dynamic Spraying at Low Pressure Followed by Friction Stir Processing,

*Mater. Sci. Eng. A*, 2012,**556**, p 114–121. - 2.
K.J. Hodder, J.A. Nychka and A.G. McDonald, Comparison of 10 μm and 20 nm Al-Al

_{2}O_{3}Metal Matrix Composite Coatings Fabricated by Low-Pressure Cold Gas Dynamic Spraying,*J. Therm. Spray Technol.*, 2014,**23**(5), p 839–848. - 3.
Y. Tao, T. Xiong, C. Sun, H. Jin, H. Du and T. Li, Effect of α-Al

_{2}O_{3}on the Properties of Cold Sprayed Al/α-Al_{2}O_{3}Composite Coatings on AZ91D Magnesium Alloy,*Appl. Surf. Sci.*, 2009,**256**(1), p 261–266. - 4.
Q. Wang, Q. Sun, M.-X. Zhang, W.-J. Niu, C.-B. Tang, K.-S. Wang, X. Rui, L. Zhai and L. Wang, The Influence of Cold and Detonation Thermal Spraying Processes on the Microstructure and Properties of Al-Based Composite Coatings on Mg Alloy,

*Surf. Coat. Technol.*, 2018,**352**, p 627–633. - 5.
C. San Marchi, F. Cao, M. Kouzeli and A. Mortensen, Quasistatic and Dynamic Compression of Aluminum-Oxide Particle Reinforced Pure Aluminum,

*Mater. Sci. Eng. A*, 2002,**337**(1), p 202–211. - 6.
Yashpal, Sumankant, C.S. Jawalkar, A.S. Verma, N.M. Suri, Fabrication of Aluminium Metal Matrix Composites with Particulate Reinforcement: A Review.

*Mater. Today: Proc.*2017,**4**(2, Part A), p. 2927-2936 - 7.
J.M. Shockley, S. Descartes, P. Vo, E. Irissou and R.R. Chromik, The Influence of Al

_{2}O_{3}Particle Morphology on the Coating Formation and Dry Sliding Wear Behavior of Cold Sprayed Al-Al_{2}O_{3}Composites,*Surf. Coat. Technol.*, 2015,**270**, p 324–333. - 8.
O. Meydanoglu, B. Jodoin and E.S. Kayali, Microstructure, Mechanical Properties and Corrosion Performance of 7075 Al Matrix Ceramic Particle Reinforced Composite Coatings Produced by the Cold Gas Dynamic Spraying Process,

*Surf. Coat. Technol.*, 2013,**235**, p 108–116. - 9.
D. Giugliano, D. Barbera, H. Chen, N.-K. Cho and Y. Liu, Creep-Fatigue and Cyclically Enhanced Cree p. Mechanisms in Aluminium Based Metal Matrix Composites,

*Eur. J. Mech. A Solids*, 2019,**74**, p 66–80. - 10.
Q. Shen, Z. Yuan, H. Liu, X. Zhang, Q. Fu and Q. Wang, The Damage Mechanism of 17vol.%SiCp/Al Composite Under Uniaxial Tensile Stress,

*Mater. Sci. Eng. A*, 2020,**782**, p 139274. - 11.
R. Balokhonov, V. Romanova and A. Kulkov, Microstructure-Based Analysis of Deformation and Fracture in Metal-Matrix Composite Materials,

*Eng. Fail. Anal.*, 2020,**110**, p 104412. - 12.
Y. Ma, Z. Chen, M. Wang, D. Chen, N. Ma and H. Wang, High Cycle Fatigue Behavior of the In-Situ TiB2/7050 Composite,

*Mater. Sci. Eng. A*, 2015,**640**, p 350–356. - 13.
P. Peng, M. Gao, E. Guo, H. Kang, H. Xie, Z. Chen and T. Wang, Deformation Behavior and Damage in B4Cp/6061Al Composites: An Actual 3D Microstructure-Based Modeling,

*Mater. Sci. Eng. A*, 2020,**781**, p 139169. - 14.
Z. Zhang, F. Liu, E.-H. Han, L. Xu and P.C. Uzoma, Effects of Al

_{2}O_{3}on the Microstructures and Corrosion Behavior of Low-Pressure Cold Gas Sprayed Al 2024-Al_{2}O_{3}Composite Coatings on AA 2024–T3 Substrate,*Surf. Coat. Technol.*, 2019,**370**, p 53–68. - 15.
R. Fernandez and B. Jodoin, Cold Spray Aluminum-Alumina Cermet Coatings: Effect of Alumina Content,

*J. Therm. Spray Technol.*, 2018,**27**(4), p 603–623. - 16.
T. Peat, A. Galloway, A. Toumpis, P. McNutt and N. Iqbal, The Erosion Performance of Particle Reinforced Metal Matrix Composite Coatings Produced by Co-Deposition Cold Gas Dynamic Spraying,

*Appl. Surf. Sci.*, 2017,**396**, p 1623–1634. - 17.
H. Assadi, H. Kreye, F. Gärtner and T. Klassen, Cold Spraying—A Materials Perspective,

*Acta Mater.*, 2016,**116**, p 382–407. - 18.
L. He and M. Hassani, A Review of the Mechanical and Tribological Behavior of Cold Spray Metal Matrix Composites,

*J. Therm. Spray Technol.*, 2020,**29**(7), p 1565–1608. - 19.
J.M. Shockley, E.F. Rauch, R.R. Chromik and S. Descartes, TEM Microanalysis of Interfacial Structures After Dry Sliding of Cold Sprayed Al-Al

_{2}O_{3},*Wear*, 2017,**376–377**, p 1411–1417. - 20.
Y.T.R. Lee, H. Ashrafizadeh, G. Fisher and A. McDonald, Effect of Type of Reinforcing Particles on the Deposition Efficiency and Wear Resistance of Low-Pressure Cold-Sprayed Metal Matrix Composite Coatings,

*Surf. Coat. Technol.*, 2017,**324**, p 190–200. - 21.
L. Gyansah, N.H. Tariq, J.R. Tang, X. Qiu, B. Feng, J. Huang, H. Du, J.Q. Wang and T.Y. Xiong, Cold Spraying SiC/Al Metal Matrix Composites: Effects of SiC Contents and Heat Treatment on Microstructure, Thermophysical and Flexural Properties,

*Mater. Res. Express*, 2018,**5**(2), p 026523. - 22.
G. Bolelli, A. Candeli, H. Koivuluoto, L. Lusvarghi, T. Manfredini and P. Vuoristo, Microstructure-Based Thermo-Mechanical Modelling of Thermal Spray Coatings,

*Mater. Des.*, 2015,**73**, p 20–34. - 23.
R.R. Balokhonov, V.A. Romanova, S. Schmauder and E.S. Emelianova, A Numerical Study of Plastic Strain Localization and Fracture Across Multiple Spatial Scales in Materials with Metal-Matrix Composite Coatings,

*Theor. Appl. Fract. Mech.*, 2019,**101**, p 342–355. - 24.
S.I. Gad, M.A. Attia, M.A. Hassan and A.G. El-Shafei, A Random Microstructure-Based Model to Study the Effect of the Shape of Reinforcement Particles on the Damage of Elastoplastic Particulate Metal Matrix Composites,

*Ceram. Int.*, 2020 https://doi.org/10.1016/j.ceramint.2020.09.189 - 25.
J. Zhang, Q. Ouyang, Q. Guo, Z. Li, G. Fan, Y. Su, L. Jiang, E.J. Lavernia, J.M. Schoenung and D. Zhang, 3D Microstructure-Based Finite Element Modeling of Deformation and Fracture of SiCp/Al Composites,

*Compos. Sci. Technol.*, 2016,**123**, p 1–9. - 26.
M.A. Eltaher and A. Wagih, Micromechanical Modeling of Damage in Elasto-Plastic Nanocomposites Using Unit Cell Representative Volume Element and Cohesive Zone Model,

*Ceram. Int.*, 2020 https://doi.org/10.1016/j.ceramint.2020.01.046 - 27.
H.K. Park, J. Jung and H.S. Kim, Three-Dimensional Microstructure Modeling of Particulate Composites Using Statistical Synthetic Structure and Its Thermo-Mechanical Finite Element Analysis,

*Comput. Mater. Sci.*, 2017,**126**, p 265–271. - 28.
R. Ekici, M. Kemal Apalak, M. Yildirim and F. Nair, Simulated and Actual Micro-Structure Models on the Indentation Behaviors of Particle Reinforced Metal Matrix Composites,

*Mater. Sci. Eng. A*, 2014,**606**, p 290–298. - 29.
A.S. Shedbale, I.V. Singh and B.K. Mishra, Heterogeneous and Homogenized Models for Predicting the Indentation Response of Particle Reinforced Metal Matrix Composites,

*Int. J. Mech. Mater. Des.*, 2017,**13**(4), p 531–552. - 30.
V. Tvergaard and A. Needleman, Analysis of the Cup-Cone Fracture in a Round Tensile Bar,

*Acta Metall.*, 1984,**32**(1), p 157–169. - 31.
G.I. Barenblatt, The mathematical theory of equilibrium cracks in brittle fracture,

*Advances in Applied Mechanicsed*. H.L. Dryden, T. von Kármán, G. Kuerti, F.H. van den Dungen, L. Howarth Ed., Elsevier, Amsterdam, 1962, p 55–129 - 32.
G.R. Johnson and T.J. Holmquist, An Improved Computational Constitutive Model for Brittle Materials,

*AIP. Conf. Proc.*, 1994,**309**(1), p 981–984. - 33.
C. Shao, C. Lo, K. Bhagavathula, A. McDonald and J. Hogan, High Strength Particulate Aluminum Matrix Composite Design: Synergistic Strengthening Strategy,

*Compos. Commun.*, 2021,**25**, p 100697. - 34.
M.P. Dewar, A.G. McDonald and A.P. Gerlich, Interfacial Heating During Low-Pressure Cold-Gas Dynamic Spraying of Aluminum Coatings,

*J. Mater. Sci.*, 2012,**47**(1), p 184–198. - 35.
Vic-2D, Correlated Solutions Inc., Irmo, South Carolina. URL http:// correlatedsolutions.com/vic-2d/.

- 36.
"STM C1425-15, Standard Test Method for Monotonic Compressive Strength of Advanced Ceramics at Ambient Temperature", West Conshohocken (PA)

- 37.
S. Ma, X. Zhuang and X. Wang, 3D Micromechanical Simulation of the Mechanical Behavior of an In-Situ Al

_{3}Ti/A356 Composite,*Compos. Part. B-ENG.*, 2019,**176**, p 107115. - 38.
M. Parsazadeh, G. Fisher, A. McDonald and J. Hogan, Computational Investigation of the Effect of Microstructure on the Scratch Resistance of Tungsten-Carbide Nickel Composite Coatings,

*Wear*, 2021,**478–479**, p 203888. - 39.
N.K. Sharma, R.K. Mishra and S. Sharma, 3D Micromechanical Analysis of Thermo-Mechanical Behavior of Al

_{2}O_{3}/Al Metal Matrix Composites,*Comput. Mater. Sci.*, 2016,**115**, p 192–201. - 40.
S. Sayahlatifi, G.H. Rahimi and A. Bokaei, The Quasi-Static Behavior of Hybrid Corrugated Composite/Balsa Core Sandwich Structures in Four-Point Bending: Experimental Study and Numerical Simulation,

*Eng. Struct.*, 2020,**210**, p 110361. - 41.
B. Amirian, H.Y. Li and J.D. Hogan, An Experimental and Numerical Study of Novel Nano-Grained (γ+α

_{2})-TiAl/Al_{2}O_{3}Cermets,*Mater. Sci. Eng. A*, 2019,**744**, p 570–580. - 42.
B. Amirian, H.Y. Li and J.D. Hogan, The Mechanical Response of a α

_{2}(Ti_{3}Al) + γ(TiAl)-Submicron Grained Al_{2}O_{3}Cermet Under Dynamic Compression: Modeling and Experiment,*Acta Mater.*, 2019,**181**, p 291–308. - 43.
S. Gatea, H. Ou, B. Lu and G. McCartney, Modelling of Ductile Fracture in Single Point Incremental Forming Using a Modified GTN Model,

*Eng. Fract. Mech.*, 2017,**186**, p 59–79. - 44.
A.L. Gurson, Continuum Theory of Ductile Rupture by Void Nucleation and Growth: Part I—Yield Criteria and Flow Rules for Porous Ductile Media,

*J. Eng. Mater. Technol.*, 1977,**99**(1), p 2–15. - 45.
G. Li and S. Cui, Meso-Mechanics and Damage Evolution of AA5182-O Aluminum Alloy Sheet Based on the GTN Model,

*Eng. Fract. Mech.*, 2020,**235**, p 107162. - 46.
A. Kami, B.M. Dariani, A. Sadough Vanini, D.S. Comsa and D. Banabic, Numerical Determination of the Forming Limit Curves of Anisotropic Sheet Metals Using GTN Damage Model,

*J. Mater. Process. Technol.*, 2015,**216**, p 472–483. - 47.
F. Abbassi, S. Mistou and A. Zghal, Failure Analysis Based on Microvoid Growth for Sheet Metal During Uniaxial and Biaxial Tensile Tests,

*Mater. Des.*, 2013,**49**, p 638–646. - 48.
V. Tvergaard, Influence of Void Nucleation on Ductile Shear Fracture at a Free Surface,

*J. Mech. Phys. Solids*, 1982,**30**(6), p 399–425. - 49.
C.C. Chu and A. Needleman, Void Nucleation Effects in Biaxially Stretched Sheets,

*J. Eng. Mater. Technol.*, 1980,**102**(3), p 249–256. - 50.
Abaqus, User’s manual, Version 6.14, Dassault Systèmes, 2014.

- 51.
N. Benseddiq and A. Imad, A Ductile Fracture Analysis Using a Local Damage Model,

*Int. J. Press. Vessel. Pip.*, 2008,**85**(4), p 219–227. - 52.
J. Venkatesan, M.A. Iqbal and V. Madhu, Experimental and Numerical Study of the Dynamic Response of B

_{4}C Ceramic Under Uniaxial Compression,*Thin-Walled Struct.*, 2020,**154**, p 106785. - 53.
M.R.I. Islam, J.Q. Zheng and R.C. Batra, Ballistic Performance of Ceramic and Ceramic-Metal Composite Plates with JH1, JH2 and JHB Material Models,

*Int. J. Impact Eng*, 2020,**137**, p 103469. - 54.
J.-T. Zhang, L.-S. Liu, P.-C. Zhai, Z.-Y. Fu and Q.-J. Zhang, The Prediction of the Dynamic Responses of Ceramic Particle Reinforced MMCs by Using Multi-Particle Computational Micro-Mechanical Method,

*Compos. Sci. Technol.*, 2007,**67**(13), p 2775–2785. - 55.
Y. Zhu, G.R. Liu, Y. Wen, C. Xu, W. Niu and G. Wang, Back-Spalling Process of an Al

_{2}O_{3}Ceramic Plate Subjected to an Impact of Steel Ball,*Int. J. Impact Eng*, 2018,**122**, p 451–471. - 56.
D.S. Dugdale, Yielding of Steel Sheets Containing Slits,

*J. Mech. Phys. Solids*, 1960,**8**(2), p 100–104. - 57.
D.V. Kubair, P.H. Geubelle and Y.Y. Huang, Analysis of a Rate-Dependent Cohesive Model for Dynamic Crack Propagation,

*Eng. Fract. Mech.*, 2002,**70**(5), p 685–704. - 58.
P. Camanho, C. Dávila, Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials, 2002

- 59.
S.I. Gad, M.A. Attia, M.A. Hassan and A.G. El-Shafei, Predictive Computational Model for Damage Behavior of Metal-Matrix Composites Emphasizing the Effect of Particle Size and Volume Fraction,

*Mater.*, 2021,**14**(9), p 2143. - 60.
G.V. Jagadeesh and S. Gangi Setti, A Review on Micromechanical Methods for Evaluation of Mechanical Behavior of Particulate Reinforced Metal Matrix Composites,

*J. Mater. Sci.*, 2020,**55**(23), p 9848–9882. - 61.
Q. Meng and Z. Wang, Prediction of Interfacial Strength and Failure Mechanisms in Particle-Reinforced Metal-Matrix Composites Based on a Micromechanical Model,

*Eng. Fract. Mech.*, 2015,**142**, p 170–183. - 62.
A. Sazgar, M.R. Movahhedy, M. Mahnama and S. Sohrabpour, Development of a Molecular Dynamic Based Cohesive Zone Model for Prediction of an Equivalent Material Behavior for Al/Al

_{2}O_{3}Composite,*Mater. Sci. Eng. A*, 2017,**679**, p 116–122. - 63.
X.-T. Luo, G.-J. Yang and C.-J. Li, Multiple Strengthening Mechanisms of Cold-Sprayed cBNp/NiCrAl Composite Coating,

*Surf. Coat. Technol.*, 2011,**205**(20), p 4808–4813. - 64.
J. Wang, C. Li, Y. Wan, C. Zhang, J. Ran and M.W. Fu, Size Effect on the Shear Damage Under Low Stress Triaxiality in Micro-Scaled Plastic Deformation of Metallic Materials,

*Mater. Des.*, 2020,**196**, p 109107. - 65.
K.S. Tun, V. Tungala, Q.B. Nguyen, J. Chan, R. Kwok, J.V.M. Kuma and M. Gupta, Enhancing Tensile and Compressive Strengths of Magnesium Using Nanosize (Al

_{2}O_{3}+ Cu) Hybrid Reinforcements,*J. Compos. Mater.*, 2011,**46**(15), p 1879–1887. - 66.
C.S. Pande and K.P. Cooper, Nanomechanics of Hall-Petch Relationshi p. in Nanocrystalline Materials,

*Prog. Mater Sci.*, 2009,**54**(6), p 689–706. - 67.
H. Lim, C.C. Battaile, J.E. Bishop and J.W. Foulk, Investigating Mesh Sensitivity and Polycrystalline RVEs in Crystal Plasticity Finite Element Simulations,

*Int. J. Plast.*, 2019,**121**, p 101–115. - 68.
J.-T. Zhang, H.-J. Shi, L.-S. Liu and P.-C. Zhai, Numerical Analysis on the Effects of Particle Configuration on the Damage and Mechanical Properties of Particle Reinforced MMCs Under Dynamic Compression,

*Comput. Mater. Sci.*, 2010,**50**(2), p 496–502. - 69.
G. Munday, J. Hogan and A. McDonald, On the Microstructure-Dependency of Mechanical Properties and Failure of Low-Pressure Cold-Sprayed Tungsten Carbide-Nickel Metal Matrix Composite Coatings,

*Surf. Coat. Technol.*, 2020,**396**, p 125947. - 70.
"Standard Test Method for Knoo p. and Vickers Hardness of Materials", E 384, Annual Book of ASTM Standards, ASTM, 2010, p. 1-42

- 71.
R. Ekici, M. Kemal Apalak, M. Yıldırım and F. Nair, Effects of Random Particle Dispersion and Size on the Indentation Behavior of SiC Particle Reinforced Metal Matrix Composites,

*Mater. Des.*, 2010,**31**(6), p 2818–2833. - 72.
ISO 6507-1: Metallic materials—Vickers hardness test—Part 1: Test method.

- 73.
Y.L. Shen and Y.L. Guo, Indentation Modelling of Heterogeneous Materials,

*Modell. Simul. Mater. Sci. Eng.*, 2001,**9**(5), p 391–398. - 74.
R. Pereyra and Y.L. Shen, Characterization of Indentation-induced ‘Particle Crowding’ in Metal Matrix Composites,

*Int. J. Damage Mech.*, 2005,**14**(3), p 197–213. - 75.
J.F. Zhang, H. Andrä, X.X. Zhang, Q.Z. Wang, B.L. Xiao and Z.Y. Ma, An Enhanced Finite Element Model Considering Multi Strengthening and Damage Mechanisms in Particle Reinforced Metal Matrix Composites,

*Compos. Struct.*, 2019,**226**, p 111281. - 76.
L. Weng, T. Fan, M. Wen and Y. Shen, Three-Dimensional Multi-Particle FE Model and Effects of Interface Damage, Particle Size and Morphology on Tensile Behavior of Particle Reinforced Composites,

*Compos. Struct.*, 2019,**209**, p 590–605. - 77.
X. Gao, X. Zhang, M. Qian and L. Geng, Effect of Reinforcement Shape on Fracture Behaviour of SiC/Al Composites with Network Architecture,

*Compos. Struct.*, 2019,**215**, p 411–420. - 78.
J.F. Zhang, X.X. Zhang, Q.Z. Wang, B.L. Xiao and Z.Y. Ma, Simulations of Deformation and Damage Processes of SiCp/Al Composites During Tension,

*J. Mater. Sci. Technol*, 2018,**34**(4), p 627–634. - 79.
S. Ma, X. Zhuang and X. Wang, Particle Distribution-Dependent Micromechanical Simulation on Mechanical Properties and Damage Behaviors of Particle Reinforced Metal Matrix Composites,

*J. Mater. Sci.*, 2021,**56**(11), p 6780–6798. - 80.
S. Amirkhanlou, M.R. Rezaei, B. Niroumand and M.R. Toroghinejad, Refinement of Microstructure and Improvement of Mechanical Properties of Al/Al

_{2}O_{3}Cast Composite by Accumulative Roll Bonding Process,*Mater. Sci. Eng. A*, 2011,**528**(6), p 2548–2553. - 81.
D. Garbiec, M. Jurczyk, N. Levintant-Zayonts and T. Mościcki, Properties of Al-Al

_{2}O_{3}Composites Synthesized by Spark Plasma Sintering Method,*Arch. Civ. Mech. Eng.*, 2015,**15**(4), p 933–939. - 82.
M. Rahimian, N. Parvin and N. Ehsani, Investigation of Particle Size and Amount of Alumina on Microstructure and Mechanical Properties of Al Matrix Composite Made by Powder Metallurgy,

*Mater. Sci. Eng. A*, 2010,**527**(4), p 1031–1038. - 83.
M.P. Reddy, F. Ubaid, R.A. Shakoor, G. Parande, V. Manakari, A.M.A. Mohamed and M. Gupta, Effect of Reinforcement Concentration on the Properties of Hot Extruded Al-Al

_{2}O_{3}Composites Synthesized Through Microwave Sintering Process,*Mater. Sci. Eng. A*, 2017,**696**, p 60–69. - 84.
H.J. Böhm and W. Han, Comparisons Between Three-Dimensional and Two-Dimensional Multi-Particle Unit Cell Models For Particle Reinforced Metal Matrix Composites,

*Modell. Simul. Mater. Sci. Eng.*, 2001,**9**(2), p 47–65. - 85.
E. Soppa, J. Nellesen, V. Romanova, G. Fischer, H.-A. Crostack and F. Beckmann, Impact of 3D-Model Thickness on FE-Simulations of Microstructure,

*Mater. Sci. Eng. A*, 2010,**527**(3), p 802–811. - 86.
H. Qing, Automatic Generation of 2D Micromechanical Finite Element Model of Silicon-Carbide/Aluminum Metal Matrix Composites: Effects of the Boundary Conditions,

*Mater. Des.*, 2013,**44**, p 446–453. - 87.
K. Dorhmi, L. Morin, K. Derrien, Z. Hadjem-Hamouche and J.-P. Chevalier, A Homogenization-Based Damage Model for Stiffness Loss in Ductile Metal-Matrix Composites,

*J. Mech. Phys. Solids*, 2020,**137**, p 103812. - 88.
A. Góral, W. Żórawski and M. Makrenek, The Effect of the Standoff Distance on the Microstructure and Mechanical Properties of Cold Sprayed Cr3C2-25(Ni20Cr) Coatings,

*Surf. Coat. Technol.*, 2019,**361**, p 9–18. - 89.
B. Liu, W. Huang, L. Huang and H. Wang, Size-Dependent Compression Deformation Behaviors of High Particle Content B4C/Al Composites,

*Mater. Sci. Eng. A*, 2012,**534**, p 530–535. - 90.
M. Groeber, S. Ghosh, M.D. Uchic and D.M. Dimiduk, A Framework for Automated Analysis and Simulation of 3D Polycrystalline Microstructures.: Part 1: Statistical Characterization,

*Acta Mater.*, 2008,**56**(6), p 1257–1273. - 91.
J. Genée, S. Berbenni, N. Gey, R.A. Lebensohn and F. Bonnet, Particle Interspacing Effects on the Mechanical Behavior of a Fe-TiB

_{2}Metal Matrix Composite Using FFT-Based Mesoscopic Field Dislocation Mechanics,*Adv. Model. Simul. Eng. Sci.*, 2020,**7**(1), p 6. - 92.
W. Kayser, A. Bezold and C. Broeckmann, EBSD-Based FEM Simulation Of Residual Stresses in a WC6wt.-%Co Hardmetal,

*Int. J. Refract. Met. Hard Mater.*, 2018,**73**, p 139–145. - 93.
W. Jibran, J. Hogan and A. McDonald, Towards Optimization of Thickness, Hardness, and Porosity of Low-Pressure Cold Sprayed WC-Ni Coatings,

*Int. J. Adv. Manuf. Tech.*, 2021 https://doi.org/10.1007/s00170-021-07500-w - 94.
J. Liu, X. Huang, K. Zhao, Z. Zhu, X. Zhu and L. An, Effect of Reinforcement Particle Size on Quasistatic and Dynamic Mechanical Properties of Al-Al

_{2}O_{3}Composites,*J. Alloys Compd.*, 2019,**797**, p 1367–1371. - 95.
A. Rasool and H.J. Böhm, Effects of Particle Shape on the Macroscopic and Microscopic Linear Behaviors of Particle Reinforced Composites,

*Int. J. Eng. Sci.*, 2012,**58**, p 21–34. - 96.
Q. Liu, F. Qi, Q. Wang, H. Ding, K. Chu, Y. Liu and C. Li, The Influence of Particles Size and Its Distribution on the Degree of Stress Concentration in Particulate Reinforced Metal Matrix Composites,

*Mater. Sci. Eng. A*, 2018,**731**, p 351–359.

## Acknowledgments

The authors gratefully acknowledge funding support from Imperial Oil (Esso), the Natural Science and Engineering Research Council Canada, the Canada Foundation for Innovation, and the Province of Alberta Ministry of Jobs, Economy, and Innovation.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

## About this article

### Cite this article

Sayahlatifi, S., Shao, C., McDonald, A. *et al.* 3D Microstructure-Based Finite Element Simulation of Cold-Sprayed Al-Al_{2}O_{3} Composite Coatings Under Quasi-Static Compression and Indentation Loading.
*J Therm Spray Tech* (2021). https://doi.org/10.1007/s11666-021-01260-5

Received:

Revised:

Accepted:

Published:

### Keywords

- Al-Al
_{2}O_{3}MMC coating - aluminum
- compression test
- damage mechanisms
- finite element simulation
- microstructure-based model