Correlation Modeling between Morphology and Compression Behavior of Closed-Cell Al Foams Based on X-ray Computed Tomography Observations

: In the last decades, great attention has been focused on the characterization of cellular foams, because of their morphological peculiarities that allow for obtaining effective combinations of structural properties. A predictive analytical model for the compressive behavior of closed-cell Al foams, based on the correlation between the morphology of the cellular structure and its mechanical response, was developed. The cells’ morphology of cylindrical specimens was investigated at different steps of compression by X-ray computed tomography, in order to detect the collapse evolution. The structure, typically inhomogeneous at local level, was represented by developing a global virtual model consisting of homogeneous cells ordered in space, that was ﬁtted on the experimentally detected structure at each deformation step. As a result, the main parameters characterizing the two-dimensional cells morphology (equivalent diameter, circularity), processed by the model, allowed to simulate the whole compression stress–strain curve by enveloping those obtained for each step. The model, ﬁtted on the previous foam, was validated by comparing the simulated stress–strain curve and the corresponding experimental one, detected for similar foams obtained by different powder compositions. The effectiveness in terms of an accurate prediction of the compression response up to the ﬁnal densiﬁcation regime has been conﬁrmed.


Introduction
Due to their cellular microstructure, metal foams have attracted considerable interest in several industrial applications for their unique morphological characteristics, which allow for joining in a single material an effective combination of structural properties (low density, high capability to absorb energy during deformation), and various functional properties [1]. The field of possible applications has been further expanded by additive manufacturing technologies, which have revealed a marked capability in fabricating structures characterized by almost all types of cell shapes [2].
Regarding the mechanical behavior, in the last decades great attention has been focused on the characterization and modelling of closed-cell aluminum foams. Uniaxial stress-strain behavior has been investigated extensively [3][4][5], also taking strain rate sensitivity into consideration [6,7], and exploiting the advantages of new manufacturing technologies in cellular structures optimization, in order to obtain greater mechanical efficiency without increasing the relative density [8].
The compressive properties of cellular materials, which can be expressed by some key parameters of the stress-strain curve, such as the elastic modulus, plateau stress, and final deformation values, are particularly important for the mechanical design of and curvature distributions on metal closed-cell structures, and to simulate their behavior under compressive loading up to the densification stage [25].
The most complex numerical models used to represent random cells arrangements and geometrical properties variability by means of a large number of elements generated by tessellation algorithms, constitute an advanced approach to the modeling of cellular structures, and requires high computational capacity. They have been recently used for investigating the effects of various morphological properties focusing on some aspects of mechanical behavior. Random Laguerre tessellations was used to generate model foams with strongly varying cell sizes and simulate the effect on elastic stiffness of the structure [26]. The cell size effects on the compressive properties of Al foams were investigated by Voronoi modeling and statistical analysis of the relationship between the micro-structure parameters and the macro-properties of foams [27]. The compression behavior of the Voronoi model was also simulated to investigate the influence of cell irregularity on elastic and plastic behavior of closed-cell foams [28]. The effects of cell-size dispersity on the Young's modulus of foam were investigated through numerical simulations on polydispersed model obtained by the random Laguerre tessellation algorithm [29]. Different magnitudes of pore size and misalignment was simulated by imposing random displacement in irregular Voronoi models, so to assess the effects on the compressive response [30].
As a general observation, the more complex computational approaches based on FE modeling are hardly able to accurately simulate the whole spectrum of regimes that arise during the compression of cellular structures (linear elastic deformation, plastic collapse, densification), generally achieving greater accuracy in the second one, although there are examples in which a good fit of the experimental curve in the densification zone [31] or an acceptable agreement in the linear elastic region [32] is obtained.
Therefore, the predictive simulation of the whole stress-strain curve based on the geometrical properties of the cellular structure remains a challenging task, and the analytical approach has not been pushed to this direction.
Focusing on closed-cell Al foams, the present work proposes an analytical model aimed for the prediction of the compressive behavior, by calculating the main parameters of the stress-strain curve and simulating it up to the final densification, based on the two-dimensional analysis of the morphological characteristics of the cellular structure, and their evolution during the collapse process.
Considering that the precursor conditions to produce the metal foams (powder composition and the parameters of the foaming process) can be correlated to the morphological properties of the final cellular structures [33], the applicative purpose of the model can be addressed to predict with good approximation the mechanical behavior of the structure, by means of a basic two-dimensional detection and analysis of the morphology resulting from the foaming conditions set.
With these purposes, the relationship between the distributions of morphological characteristics of the foam sections longitudinal to the load direction and the mechanical response of the cellular structure, during the compressive progression, was searched for. As a result, an analytical model for morphology-behavior correlation, comprising a procedure for constructing the stress-strain curve, was defined, with the objective of simulating the whole compressive behavior, up to the final densification.
Starting from the consideration that foams generally show a typical non-homogeneous morphology, being characterized by cells having geometry and size different from each other [34,35], the correlation between the morphological properties of the foam and the characteristic parameters of the stress-strain curve (elastic modulus, plastic plateau stress, strain parameters of densification) needs to transform the real cellular structure, nonhomogeneous at local level, into a virtual reference model consisting of cells ordered in space, with global homogeneous geometric properties. The translation of the local morphological parameters into the global parameters of the reference model was obtained by introducing an intermediate model to consider the asymmetry of the cells and their orientation with respect to the direction of compression.
To investigate the potential of the model for the prediction of compressive behavior based on the morphological analysis, it was first fitted on the experimentally detected behavior for a specific foam, following the evolution of morphology during the progression of deformation; subsequently, it was used to simulate the compression behavior of similar foams. For this purpose, closed-cell Al foams were produced by an in-house process based on the compacted powder method, starting from Al and SiC powders with the addition of TiH 2 powder as foaming agent at various contents. By this way it was possible to fit the model on the experimentally detected compression behavior for a foam with a specific powder composition and evaluate its potential in predicting the same behavior for similar foams obtained by different powder compositions.
The morphology of the cellular structure, expressed by the distribution of the values of the equivalent diameters and circularities of the cells through the specimen sections, was investigated by X-ray computed tomography (CT), a technique that, together with basic radioscopy, has already proved particularly suitable for the internal investigation of metal foams, focusing on different aspects: foaming process efficiency [36], morphological analysis, and microstructural characterization of foamed structures [34,35], and their mechanical behavior [37,38]. Specifically, in [39], significant changes in the structure of the internal cells during the deformation process were documented by means of CT images, overcoming previous approaches based on combination of metallographic images and static theory application [13].
For the specimen used for model fitting, the tomographic observations were carried out at different levels of compression in order to study the evolution of cells' morphology during deformation up to the final densification, according to an experimental procedure first presented in [40] for Al foam, and in [41] for a Cu tube filled with an Al alloy foam.
Once fitted on the basis of the CT observations as compression progresses, the analytical model was used to predict the mechanical behavior of the other similar foams, produced by varying the content of the TiH 2 forming agent, and was validated by means of comparison between simulated and experimentally detected stress-strain curves.
The paper has been structured according to three main sections. In Section 2, the materials and methods are detailed. Particularly, the approach to experimental investigation is treated, specifying specimens' fabrication and characteristics, the statement of compression testing and tomographic observations, and the method of conducting the morphology analysis. Section 3 is dedicated to the theoretical development of the morphology-behavior correlation model, reported according to its main steps: the modeling of cellular structure, the simulation of compressive behavior, and the fitting of the correlation model. Section 4 reports the application of the model and the results obtained, focusing on the predictive use of the fitted model. The significance and implications of the main results are discussed.

Specimens Production
The cylindrical specimens, made of Al foam with closed cells, were produced through the compacted powder method (Figure 1), assuming the process parameters developed in previous work [42]. Commercially pure Al powder (44 µm average diameter) was mixed with TiH 2 powder as a foaming agent (5 µm average diameter) and SiC powder (37 µm average diameter), added to increase viscosity and stabilize porosity by acting at the interface metal-bubble. Three different powder compositions were considered, by varying the forming agent content: TiH 2 at 0.1-0.2-0.4 weight %, SiC at 2.8 weight %, and commercially pure Al to balance. The powders were carefully mixed to produce after pressing homogeneous precursors, because, especially for the TiH 2 particles, this is a fundamental requirement to produce foams of good quality [43].
The powders were blended in a rotating jar and compacted at room temperature in a mold, by means of a hydraulic press, at a uniaxial pressure of about 270 MPa to get cylindrical foaming precursors (diameter 15 mm, height 8 mm), ready to be foamed. The precursors were subsequently placed in a copper crucible inside an oven and heated from room temperature to the foaming temperature of 700 • C in 120 s to avoid an initial excessive gas dispersion. After a holding time equal to 5 min, the specimens were water cooled.
The foaming temperature was not much higher than the melting point of pure Al (660 • C). This combination of temperature and holding time inside the mold was chosen because it determines a good compromise between viscosity in the liquid state (melt viscosity would be too low at higher temperatures) and foam growth.
Finally, in order to obtain a cylindrical geometry suitable for the compression test, the foamed specimens were post-processed by slight compression between two planes, in order to flatten their bases (Step 0 state).

Compression Tests
The axial compression tests were performed by an MTS Insight 50 electromechanical testing machine (MTS Systems Corporation, Eden Prairie, MN, USA), under control of deformation, starting from Step 0 state. A crosshead constant speed of 2 mm/min with a data acquisition frequency of 5 Hz was utilized.
For the 0.4% TiH 2 foam the compression test was performed at increasing steps of deformation, in order to detect the evolution of the morphological characteristics during deformation. For each step, the specimen was loaded, unloaded, and observed by X-ray computed tomography (CT) to analyze the morphology of the cellular structure corresponding to the compression state at the end of the step.
The deformation reached at each step was calculated with respect to the height of the specimen in the initial condition (Step 0), and the loading curves of the steps were assembled in sequence in order to obtain the continuous stress-strain curve that represents the overall compression behavior.

Tomographic Observations and Morphology Analysis
The starting cells morphology was observed by means of variable focus computed tomography (CT) for all the specimens, while the evolution during compression was investigated for the 0.4% TiH 2 specimen only. It allows for building a virtual model of the specimen at the starting state (Step 0) and at each compression step with a high density of points (some millions), giving information on the internal volume that otherwise would be not accessible [44]. The specimen to be inspected is put in rotation within the radiation cone produced by an X-ray source. The distribution of the radiation intensity that crosses the specimen is measured by a flat panel detector and digitally acquired.
For the tomographic observations, a YXLON device-model Y-CT Vario (YXLON International GmbH, Hamburg, Germany) was used, working at an acceleration voltage adjustable up to a maximum of 225 kV, in the range microfocus-macrofocus by changing the spot size from 250 to 800 µm, in order to give priority to the spatial resolution (from 30 µm to 200 µm), or to the penetration capacity in materials with high absorption coefficient (up to 10 mm in the case of ferrous alloys), or to an object under examination with great size (up to diameter of 180 mm and height of 250 mm). In our specific case, the observations on the Al foam specimen were carried out at an acceleration voltage of 150 kV, with a spot size equal to 250 µm. The detector system was a flat panel with a resolution of 1920 × 1536 pixel.
The specimen was rotated at increments of 0.5 • /s, until a full rotation of 360 • was achieved, while it was scanned by a conical X-ray beam. The data set returned by the flat panel was processed through the VGStudio Max 2.0 (Volume Graphics GmbH, Heidelberg, Germany) software, in order to obtain a complete 3D reconstruction of the specimen.
Finally, 2D images of the specimen axial sections were considered for our work. Figure 2 shows representative cross and longitudinal sections, and CT rendering, for reference specimen (0.4% TiH 2 foam) at Step 0 state: it can be observed non-homogeneous cells that are far from circular. The 2D images of the longitudinal sections, taken at the various steps of deformation considered, were pre-processed to enhance the contours of the walls and porosity. Subsequently, a quantitative analysis of the areas A and perimeters P of the cells in section were performed by means of ImageJ 1.53c software (Wayne Rasband, National Institute of Health, Bethesda, MD, USA).
From these geometric parameters of the cells, it is possible to characterize the morphology of the cellular structure, expressed through the distribution of the values of the equivalent diameters d eq and circularities C of the cells through the section. These parameters are expressed by: The morphological analysis of the 0.4% TiH 2 foam was performed for the same section at the end of each compression step, to assess the evolution of morphological properties during the overall compression.

Morphology-Behavior Correlation Model Setting
To formulate a prediction of the compression behavior of the cellular structure based on the morphological characteristics, it is necessary to define a model for the correlation between the morphological characteristics (i.e., the distributions of the equivalent diameters and of the circularities of the cells), and the parameters that define the mechanical behavior in compression. As for cell foams, the latter are commonly identified in the following parameters of the stress-strain compression curve: elastic modulus, plastic plateau stress, and densification strain.
In Figure 3, the typical compressive stress-strain curve for foams is reported, with reference to elastic-plastic behavior such as that of metal foams [15]. Three different regimes can be identified:

•
Linear elastic regime at low stresses, for which the value of Young's modulus E f can be defined as the initial slope of the curve. • Plastic collapse regime, associated with a long plateau at a stabilization stress value σ pl , that is due to the combination between different collapse mechanisms at the cell walls (elastic buckling, plastic hinging, brittle crushing), depending on the nature of the material. • Densification regime, with an interruption of the plateau region due to structure crushing and compacting, and rapid rise in stress with further strain; this region is usually characterized by a deformation parameter (densification strain ε d ). The same curve can be assumed as a reference for both open and closed-cell foams, although the latter show more complex behavior particularly in the plateau regime, where the σ pl stress is less stable, and tends to rise with increasing strain. With regard to the densification parameter ε d , its identification for metal foams is conditioned by the ductile elasto-plastic behavior during the walls' collapse, because of which some amounts of pores may still exist in crushed structure [45], so a full densification strain is difficult to be defined. For a univocal identification it can be assumed as densification strain ε d the intersection of the two tangents to the stress-strain curve for the plastic collapse regime and the densification regime [6]. Finally, although in actual closed-cell structures behavior there is no clear-cut transition between the regime of cell walls plastic deformation and the regime of structure compacting and densification, a further parameter can be introduced, that is the strain corresponding to the densification onset ε o , starting from which it is possible to consider the development of the post-plateau compression behavior [15].
Among the models that have been developed to describe the mechanical behavior of cellular structures, the constitutive models are not suitable for correlating the morphological properties of the cellular structure and its mechanical behavior, since they represent the cellular structures as solid materials [14]. The representative models previously introduced in Section 1, instead, are based on abstractions of the cellular structure that simplify and order the geometric characteristics and are the most suitable to correlate the morphological properties and the mechanical behavior. They allow to estimate the relative density ρ f /ρ s (being ρ f the foam density, and ρ s the density of the solid metal of which the foam is made) as a function of significant geometric properties that characterize the cellular structure of the foam. Subsequently, the parameters of the stress-strain curve can be calculated as a function of the relative density ρ f /ρ s .

Cellular Structure Modeling
As highlighted in the introduction (Section 1), the representative models for closed-cell foams vary significantly according to the geometry of the polyhedra used. The skeleton cubic cell model [15] is particularly known because it allowed to define a set of equations frequently used for the estimation of the main compression behavior parameters (elastic modulus E f , plastic plateau stress σ pl , densification strain ε d ) as functions of the relative density ρ f /ρ s . This geometric schematization of the cellular structure is not very representative of the real structure of some types of foams, such as the one in question. The truncated cube model [16], which is obtained by assembling symmetric cubic cells truncated at the vertices, so to build a packed structure of small and large cells, allows for defining analytical formulation to calculate the crushing strength, from which it is possible to derive the plateau stress σ pl as a function of the plastic flow stress σ o .
To obtaining a more realistic geometric representation, the cruciform-hemisphere model [17] improves the truncated cube model by replacing the pyramidal sections at the vertex of the truncated cubes with hemispherical sections.
Starting from the latter model, Hasan's weaker-cruciform-stronger-hemisphere model [19] proposes a further improvement in the accuracy of the geometric representation of the closed-cell structure. This improvement is based on the observation of the material distribution in the cellular structures, which tends to accumulate near the junctions of cells, so that the junction regions of larger cells (where the smaller cells are located) are more massive than the larger cell walls. As a consequence, the walls of the smaller cells tend to be thicker than those of the larger cells. Introducing the thickness differentiation between the hemispherical (smaller) and cruciform (larger) sections in the cruciform-hemisphere model, the model in Figure 4 can be obtained, and assumed as reference model for the closed-cell Al alloy foam in the present study. As is evident from the representation in the figure, this type of model is based on some simplifying assumptions:

•
The cells can be divided into two size classes (large and small cells).

•
The smaller cells have high circularity and are characterized by spheroidal shapes; the larger cells tend to lose circularity and are characterized by prismatic shape.

•
The wall thicknesses of the smaller cells are greater than the larger ones.
Introducing this model, the relative density ρ f /ρ s can be calculated as a function of the diameter d and the wall thickness t s of the hemispherical section of the smaller spheroidal cells, and the width w and the wall thickness t p of the cruciform section of the larger prismatic cells [19]: Ordered and uniform models such as that in Figure 4, commonly used in numerical simulation by finite elements, assume that the real cellular structure, typically not homogeneous at the local level, can be represented by a global virtual model consisting of cells ordered in space, having homogeneous geometric properties. Under this assumption, the geometric parameters w, d, t s , t p in Equation (3) become global parameters, which must be derived from the distributions of the local morphological parameters d eq and C, expressed for each cell by the Equations (1) and (2), respectively.
The translation of the local morphological parameters into global parameters of the virtual reference model is obtained through the following three-stages procedure.

Stage 1
It is assumed that the results of the morphological analysis on a section of the cellular structure (e.g., section XZ in Figure 5a) are the distributions of the equivalent diameter d eq and the circularity C throughout the cells in section. The distribution of d eq can be expressed by the vector V(d eq ) of the ranges of the d eq values defined for the frequency analysis, sorted by increasing values, and the vector of frequency distribution D(d eq ) per range of values: where d eq i is at the same time the maximum value that limit the i-th range, and the minimum value that limit the subsequent (i + 1)-th range, N deq i is the number of cells with d eq within the i-th range, and n is the number of ranges identified for the frequency analysis. If d ts eq is the threshold value to differentiate the values of d eq to be attributed to the smaller spheroidal cells and the larger prismatic cells in the model of Figure 4, the following average values of d eq can be attributed to the two classes of cells: where d s eq is the average value of d eq for the spheroidal cells, d p eq is the average value of d eq for the prismatic cells, and m is the number of ranges sorted by increasing values up to the range that includes the threshold value d ts eq; the latter must be understood as a partition coefficient between the two size classes of cells and can be identified by prefixing the ratio d p eq /d s eq on the basis of experimental observation. The second result of the morphological analysis, that is the distribution of circularity C throughout the cells in section, can be processed analogously. As in Equations (4) and (5), it can be expressed by the vector V(C) of the maximum values that limit each range of C values defined for the frequency analysis, sorted by increasing values, and the vector of frequency distribution D(C) per range of values: By setting the threshold value C ts as partition coefficient of C values between the prismatic and the spheroidal cells, it is possible to obtain the average values of the circularity for the spheroidal cells (C s ) and the prismatic cells (C p ) by means of expressions similar to (5): In this case the threshold value for the partition between the two size classes of cells can be fixed as C ts = 0.785, that is the maximum value of circularity for prismatic cells (when the section is square). Therefore, the circularities C i will be assigned to prismatic or spheroidal cells, depending on whether C i < C ts or C i > C ts respectively. As a consequence, in Equation (7), the expressions for C p and C s are inverted, if compared to the expressions of d p eq and d s eq in Equation (5); m remains the number of ranges sorted by increasing values up to the range that includes the threshold value.

Stage 2
The couples of average values (d p eq ,C p ) and (d s eq ,C s ), for the prismatic and spheroidal cells respectively, obtained by the distributions of the local morphological parameters of the real cellular structure, are used to define the geometric parameters of an intermediate virtual model, where the sections of large and small cells are not yet circular and square such as in the final reference model of Figure 4, but instead they are elliptical and rectangular respectively (Figure 5b). By means of this intermediate model it is possible to partly take into account the asymmetry of the cells, and their orientation with respect to the z axis of the analyzed section, that is parallel to the direction of compression.
For this purpose, in a similar way to how the equivalent diameter of the section of a cell is defined as a function of its area according to Equation (1), the equivalent semiaxes ϕ rct eq and ϕ ell eq are defined for cells with rectangular section and elliptical section respectively: where k rtc and k ell are the ratios between the two semi-axes of the rectangular and elliptical section respectively.
Combining Equation (8) with Equation (1) it is possible to express ϕ rct eq and ϕ ell eq as functions of the equivalent diameter d eq . ϕ rct eq = π/16 k rct · d eq ϕ ell eq = 1/4 k ell · d eq (9) Introducing in Equation (2) the ratios between the two semi-axes of the rectangular and elliptical section k rtc and k ell , the corresponding expressions of circularity C can be obtained: For each of the Equation (10) it is possible to obtain a pair of solutions: Imposing C = C p in the first Equation (11), and C = C s in the second one, being C p and C s calculated by Equation (7), the two pairs of solutions (k rct+ , k rct− ) and (k ell+ , k ell− ) can be obtained. They express k rct and k ell respectively, according to whether they represent the ratio between the major and minor semi-axes, or vice versa.
Substituting in Equation (9) the values obtained by Equation (11) for (k rct+ , k rct− ) and (k ell+ , k ell− ), and the values of d p eq and d s eq calculated by the analysis of the distributions with the Equation (5), it is possible to obtain the corresponding values of major and minor equivalent semi-axes (ϕ rct+ eq , ϕ rct− eq ) and (ϕ ell+ eq , ϕ ell− eq ) for both rectangular section cells and elliptical section cells (Figure 5b).

Stage 3
The intermediate virtual model with rectangular section large cells and elliptical section small cells characterized by equivalent semi-axes (ϕ rct+ eq , ϕ rct− eq ) and (ϕ ell+ eq , ϕ ell− eq ) can be now translated in the reference model (Figure 5c), where square section cells and circular section cells are characterized by the global parameters width w and diameter d calculate as: being ϕ rct* eq and ϕ ell* eq the average values of (ϕ rct+ eq , ϕ rct− eq ) and (ϕ ell+ eq , ϕ ell− eq ), respectively. Coefficients γ p and γ s can assume values in the ranges [ϕ rct− eq /ϕ rct* eq , ϕ rct+ eq /ϕ rct* eq ] and [ϕ ell− eq /ϕ ell* eq , ϕ ell+ eq /ϕ ell* eq ]. By this way, as the coefficients γ p and γ s vary, it is as if the inclination of the rectangular and elliptical section cells varies, placing itself between the orthogonal and parallel positions with respect to the z axis of the section, that is the loading direction (Figure 5b). As a consequence, the global parameters w and d expressed by Equations (12) and (13), are calculated as the projection of the whole axes of the cell section, on the x direction orthogonal to the loading direction.
Finally, assuming that l x is the width of the section of the cellular structure in the x direction, the average values of the thickness of prismatic and spherical cells t p and t s can be estimated: t s = λ t p (15) where n x = int(l x /w) is the number of prismatic cells constituting the virtual structure in the x direction, and λ = t s /t p is the ration between the thickness of spherical and prismatic cells (λ >1 according to the premises on the reference model).
Once the coefficients γ p , γ s , λ have been fixed by fitting on experimental data (as will be shown in detail below), it is possible to define all the global parameters of the reference model w, d, t p , t s , and to calculate the value of the relative density ρ f /ρ s by means of the Equation (3).

Simulation of Compression Behavior
By using representative models, the parameters of the stress-strain curve (Figure 3) can be correlated to the relative density of the foam ρ f /ρ s and calculated, and so the compression behavior can be simulated. The scaling relations by Ashby et al. [9], derived for Al alloy closed-cell foam by previous generalized formulation based on the skeleton cubic cell model [15], can be used to calculate the Young's modulus E f , plateau stress σ pl , and densification strain ε d : where c 1 , c 2 , and c 3 are coefficients to be fitted experimentally (0.1≤ c 1 ≤1.0, 0.1≤ c 2 ≤1.0, 0.9≤ c 3 ≤1.0), E s and σ c are the Young modulus and the compressive strength of the solid metal of which the foam is made. For Al foams characterized by low density (ρ f /ρ s <0.2) the expressions of E f and σ pl proposed by Simone and Gibson and based on a tetrakaidecahedral cell model [46,47] could be preferred: Focusing on the type of geometrical model used, the formulation for σ pl by Santosa and Weirzbicki [16] based on the truncated cube model, and expressed as a function of the plastic flow stress σ o , results to be more reliable, because of the similarity between this model and the cruciform-hemisphere type model by Hasan [19] assumed here as a reference model: where the term inside the second square bracket is the plastic flow stress σ o , expressed as a function of the yield strength σ y , the ultimate strength σ u , and the strain hardening exponent ν of the cell wall material. Finally, the densification onset strain ε o can be estimated considering that the stressstrain curve does not start to rise until the foam is compacted to a relative density of about 0.5 [15]: Starting from ε = ε o it is possible to simulate the post-plateau compression curve by means of the following equation [15]: where ξ and ψ are coefficients to be fitted on experimental curves.
For all formulations, increasing ρ f /ρ s increases the elastic modulus E f and the plateau stress σ pl , and reduces both the strain parameters of densification ε d and ε o .
In this regard, it should be pointed out that the equations for calculating the relative density ρ f /ρ s as a function of the geometric properties of the representative models, such as Equation (3), assume that the density of the walls is equal to that of the solid material. A reduction in its compactness, for example, caused by microporosity, would entail a reduction in the relative density of the cellular structure, with the direct consequences on the values of the compression parameters highlighted above (in this case, the decrease in E f and σ pl , and the increase in ε d and ε o ). To extend the use of Equation (3) by taking this aspect into account, it is possible to introduce experimentally-defined coefficients for the thicknesses t p and t s , which express the reduction of density in the cell walls.

Fitting of the Correlation Model
To use the proposed morphology-behavior correlation model, the following coefficients must be defined by experimental fitting: • γ p and γ s in Equations (12) and (13)  To set these coefficients, the theoretical values of compression behavior parameters E f , σ pl , ε d , calculated by applying the model on morphological properties and using the equations reported in the previous Section 3.3, must be fitted on the corresponding values identified on the experimental compression curves.

Compression Test and Morphology Analysis
Starting from the 0.4% TiH 2 specimen in the initial conditions (Step 0), deformation steps were performed at 10% increasing steps of deformation, from Step 1 to Step 6. The geometric data characterizing the steps are reported in Table 1. Figure 6 summarizes the main results of the whole deformation process, including the loading-unloading curves (expressed in N) corresponding to the six steps, and reporting the final states of compression corresponding to each step, reconstructed by CT to perform the morphological analysis of the cellular structure.
The sequence of the loading curves resulting from each compression step were assembled obtaining the continuous stress-strain curve (expressed in MPa in the same Figure 6). The overall deformation ε reached at the end of each step was calculated as the ratio between the shortening ∆l (mm) and the initial height of the specimen at Step 0. A final step of ultimate compression up to an advanced collapse condition was performed to plot the trend of the curve after densification onset. In this case the morphological analysis corresponding to the final state loses significance due to the high degree of collapse of the cellular structure.  As in the schematic curve of Figure 3, the experimental compression curve in Figure 6 is characterized by an initial region of elastic behavior followed by a substantially stable plateau along which most of the plastic deformation is concentrated. This plateau extends up to an overall deformation of 40% (about the end of Step 4), beyond which the load begins to grow quickly together with the collapse of the cells in the whole sample.
The morphological analysis of the cellular structure was performed on the two longitudinal sections, arranged orthogonally with respect to each other and containing the rotation z axis of the specimen, that is the loading axis, reported in Figure 2.
A preliminary qualitative observation of these sections (corresponding to Step 0 state) allows to verify the morphological preconditions that support the choice of the reference model in Figure 4 for the geometry of the cellular structure:

•
The sections present cells diversified according to different size classes, which can be reduced to small and large cells.

•
The former have greater circularity and tend to spheroidal shape, whereas the latter lose circularity, and tend to a prismatic shape.

•
The wall thicknesses of the smaller cells are generally greater than the larger ones.
The images of these sections were recorded and analyzed not only at Step 0 state, but also at the end state of each compression step, to measure the area A and perimeter P of each cell, and the distributions of the equivalent diameter d eq , and circularity C through the section by means of Equations (1) and (2). This allows to analyze the evolution of morphological properties during overall deformation process.
Taking into consideration here the case of the XZ section (being the modality of analysis for the YZ section strictly similar), its evolution through the end states of the compression steps is shown by the CT images in Figure 7. At the beginning, deformation occurs in the weakest zones (Step 1), where some cells assume oblong shapes. Subsequently compression led to the progressive collapse of larger cells (Steps 4 and 5), giving rise to the specimen densification. The morphological characteristics of the cellular structure are well described by the evolution of equivalent diameter d eq and circularity C distributions during the subsequent compression steps (Figures 8 and 9, where the values on the horizontal axes represent the upper limit of the ranges for frequency distribution).
As the compression advances, in general the cell size is reduced; at each deformation step it can be observed that a certain number of large cells disappear. In particular, the number of cells with the smallest equivalent diameter d eq increases from Step 0 to Step 1, decreases from Step 1 to Step 2 and increases again from Step 3 to Step 5. Before the final Step 6, the distribution is grouped around low values of d eq .
The circularity C shows a trend in agreement to that of the equivalent diameter d eq : the number of cells with C≈1 increases from step 0 to step 1 and remains at a high value during the following steps, allowing to deduce that as densification advances the cells tend to become small, assuming an almost spherical shape. However, at the last observed step the number of spherical cells is shown to not be very high.

Correlation Model Fitting and Mechanical Behavior Simulation
Starting from the morphological characteristics expressed by the evolution of equivalent diameter d eq and circularity C distributions collected in Figures 8 and 9, for each compression step, it is possible to apply the morphology-compression behavior correlation model proposed in Section 3.
The distributions of d eq and C can be expressed by the vectors V(d eq ) and D(d eq ) in Equations (4), and vectors V(C) and D(C) in Equation (6), respectively.
The partition of d eq i values between the two size classes of cells has been performed by prefixing the ratio d p eq /d s eq on the basis of image analysis of the CT sections and using Equation (5) to calculate the values of d p eq and d s eq corresponding to the partition that respect the prefixed ratio.
The values of C s and C p were calculated by Equation (7), by fixing the threshold value C ts = 0.785, as suggested in Section 3.1.
The pairs of average values (d p eq ,C p ) and (d s eq ,C s ) calculated as above for the prismatic and spheroidal cells respectively, were used in Equations (9) and (11), to obtain the values of major and minor equivalent semi-axes (ϕ rct+ eq , ϕ rct− eq ) and (ϕ ell+ eq , ϕ ell− eq ) for both rectangular section cells and elliptical section cells of the intermediate virtual model (Figure 5b), and the corresponding average values ϕ rct* eq and ϕ ell* eq . Finally, the global parameters w and d for the reference model (Figures 4 and 5c) were calculated by Equations (12) and (13), as functions of coefficients γ p and γ s , and the average values of the thickness of prismatic and spheroidal cells t p and t s were calculated by Equations (14) and (15), as functions of their ratio λ = t s /t p . The coefficients γ p , γ s , λ need to be fitted on experimental data.
Defining the global parameters of the reference model w, d, t p , t s , it is possible: • to calculate the value of the relative density ρf/ρs by means of the Equation (3). • to simulate the compression stress-strain curve by calculating the parameters Ef, σpl, εo, εd using the most suitable equations that correlate them to relative density, among those collected in the Section 3.3.
Here, the general scaling relations in Equations (16) and (18) were used to calculate E f and ε d respectively. The alternative Equation (19) for E f was discarded as it was shown to be unsuitable in the present case, being the experimental value of the relative density ρ f /ρ s >0.2 for the 0.4% TiH 2 foam (ρ f /ρ s = 0.220 at Step 0 state). Equation (21) was used to calculate σ pl , being the most appropriate for cruciform-hemisphere type model such as the one adopted as reference. The use of Equation (17) was avoided because it has been found that it overestimates the value of σ pl precisely in the case of Al foams obtained by compacted powder method and TiH 2 as foaming agent [48]. Finally, Equation (22) was used to calculate ε o , and Equation (23) to simulate of the post-plateau compression curve.
Overall, the coefficients γ p and γ s in Equations (12) and (13), λ in Equation (15), c 1 and c 3 in Equations (16) and (18), ξ and ψ in Equation (23), which complete the morphologybehavior correlation model, must be defined by experimental fitting. For this purpose, the theoretical values of compression behavior parameters E f , σ pl , ε d were calculated varying the coefficients γ p , γ s , λ, c 1 , c 3 , and are compared to the corresponding values identified on experimental compression curve in Figure 6, following the indication of Figure 3. In particular, the value of stress at the beginning of the plateau was assumed for σ pl . The fitted values of the coefficients are those who minimize the sum of squared distances between the theoretical and experimental values of compression behavior parameters.
The value of relative density ρ f /ρ s corresponding to the fitted parameters was assumed as relative density of the cellular structure.
Finally, the values of post-plateau compression parameters ξ and ψ in Equation (23) were defined by imposing σ = σ pl at ε = ε o and fitting the slope of the final part of the theoretical densification curve on the experimental one.
In Tables 2 and 3, the values of the fitted coefficients, global parameters of the cellular structure, relative density calculate by Equation (3), and parameters of simulated stressstrain curve are reported for all compression steps performed on 0.4% TiH 2 specimen. As a first validation of the model, the theoretical value of the relative density ρ f /ρ s calculated in correspondence with the morphology of Step 0 (0.217) is well fitted with the experimental value equal to 0.220.

Relative Density Parameters of Simulated Curve
Step 0 0.000 0.217 0.886ε 3 − 0.161ε 2  In Figure 10, the compression curve simulated by the theoretical values of the parameters for the morphology at Step 0 (reported in the first row of Table 3) was superimposed on the experimental curve. Two main observations can be made. The apparent discrepancy of the Young's modulus obtained through the theoretical model, in comparison to the value estimated by assessing the initial slope of the experimental curve, is due to a known uncertainty in the estimate of the modulus E f . In fact, although the initial loading condition in compression test could appear to be linear elastic, in reality the corresponding stretch of the curve has a non-straight trend, as some cells yield at very low loads. Therefore, at the beginning of the real stress-strain curve, at the local level, there is a tendency towards a greater slope than the apparent global one [9]. As a consequence, the value of E f calculated as the slope of the experimental curve is underestimated compared to the true modulus. Anyhow, very large differences between the analytical predictions and the experimental values, particularly for Young's modulus, are usually found, and they are due to the frequent presence of defects in the real foam structure (uneven cell wall curvature, missing cell walls, inhomogeneous density distribution) which are not included in the theoretical models [18].
The second observation concerns the plateau regime, and in particular what was already highlighted when the typical compressive stress-strain curve for closed-cell foams was introduced (Figure 3, Section 3.1): the complex behavior in the plastic collapse regime is characterized by unstable σ pl stress, which tends to rise with increasing strain (as confirmed by the experimental curve). This tendency makes it necessary to fit the correlation model between morphology and mechanical behavior according to several compression steps, as the morphological analysis limited to just one compression state, such as the initial state at Step 0, only would allow to define a strictly constant value plateau, resulting in an overall simulation far from the real behavior, as shown in Figure 10. Figure 11 shows the stress-strain curves simulated by the theoretical values of the parameters obtained applying the correlation model to the morphologies at the end state of each compression step (Table 3 on the right). In this case it can be observed that: • By enveloping these curves, it is possible to define the rise trend of the plateau stress.

•
The sections of the curves in densification regime tend to converge. It is therefore possible to trace the tangents to the plastic collapse plateau and the convergence curve to densification, identifying a densification strain value of simulation convergence ε d '.
Finally, the univocal theoretical stress-strain curve can be simulated by these construction criteria:

•
The elastic regime at low stresses is assumed to be represented by the linear regime simulated for Step 0 morphology.

•
The plastic collapse regime is obtained by enveloping the start points of the plateau of the simulated curves in Figure 11.

•
The start point of the plateau for Step 6 simulated curve is assumed as the end point of plastic collapse regime and start point of post-plateau densification regime, so that the corresponding strain is assumed as densification onset strain ε o ' of the enveloped stress-strain curve.

•
The densification regime is obtained by iteration of densification curve calculation, starting from ε = ε o ' and imposing ε d = ε d ' in Equation (23).
As a result, the simulated curve in Figure 12 is obtained. If compared to the experimental one, a good level of approximation is found, quantified by the congruence between the theoretical and experimental values of the densification strain (0.715 and 0.727, respectively), and by a deviation in stress values of the densification curve within 10%, due to a strain shift within 0.02.

Robustness of the Correlation Model
The results presented above refer to the morphological analysis of the section XZ of the 0.4% TiH 2 specimen, according to the scheme in Figure 2. The calculation of the relative densities for each compression step, on the basis of the corresponding morphological parameters, was applied also to section YZ in the same figure, obtaining negligible deviations, and simulation curves consequently similar to those shown before.
To highlight the effect of morphological variations within the same specimen on the simulation of the compression curve, in the following Figures 13 and 14 the distributions of the morphological properties D eq (Figure 13a,b) and C (Figure 14a,b) are compared for both sections XZ and YZ of the reference specimen (Specimen 1) at the initial condition (Step 0). The comparison is based on the analysis of the main statistical parameters that characterize the distribution of frequencies (i.e., the number of cells per interval of morphology property): the first quartile (Q1), second quartile or median (M), and third quartile (Q3). Since the three quartiles are percentiles of 25%, they divide the set of intervals of each distribution of D eq and C into four subsets of equal extension, i.e., equal incidence on the overall distribution. It can, therefore, be observed that, despite the differences in the distributions of frequencies and in the overall number of cells, the values of the statistical parameters remain substantially stable, both for D eq and C, and thus the distributions are affine, and so are the morphological characteristics of the two sections.  This affinity is reflected in the simulation of the compressive behavior, as confirmed by the corresponding effect on the values of the relative density and the compression curve parameters, calculated for both the XZ and YZ sections at the Step 0 state, using the same procedure described in Section 4.2. As highlighted by the data reported in Table 4, the most significant deviations, detectable for σ pl and ε o , are limited and negligible in the simulation of the compression curve. Table 4. Relative density and parameters of simulated compression curve: comparison between longitudinal sections XZ and YZ of reference Specimen 1, as defined in Figure 2.

Sections
Relative An aspect related to the previous one deals with the reproducibility of results for specimens with the same properties, i.e., made using the same foaming method, with the same composition of the powders, and under the same process conditions. In this regard, the effects of morphological variations on the simulation of the compression curve were analyzed. In Figure 13a,c,d and Figure 14a,c,d the distributions of the morphological properties D eq and C, respectively, corresponding to the XZ section of the reference specimen (Specimen 1) are compared with those characterizing randomly chosen longitudinal sections of other two specimens (Specimen 2 and 3). As regards the distribution of D eq , it is possible to note a close affinity between Specimens 1 and 2 (Figure 13a,c); in Specimen 3 a distribution of D eq with a tendency towards larger cells emerges (in Figure 13d the third quartile Q3 is shifted on the right). The distributions of C are substantially homogeneous among the three specimens; in this case the most marked affinity is that between Specimens 2 and 3 (Figure 14c,d), while the distribution of Specimen 1 shows a tendency towards higher C values (in Figure 14a all the three quartiles are compacted on the right).
Also in this case the similarity between the distributions of the morphological properties is reflected in the simulation of the compressive behavior, as confirmed by the limited variances in the values of the relative density and the compressive parameters calculated for the three specimens and reported in Table 5. This result is in agreement with previous studies on the reproducibility of compressive properties of closed-cell aluminum foam, which is considered stable for values of the ratio between the height of the specimen and the mean size of the cells above 6 [10] (in the cases examined here, this ratio calculated using the third quartile of the D eq distribution instead of the mean value, to take into account the presence of greater size cells, assumes values in the range 10.1 ÷ 14.6). Table 5. Relative density and parameters of simulated compression curve: comparison between longitudinal sections of Specimen 1 (XZ), Specimen 2, and Specimen 3.

Sections
Relative Density Parameters of Simulated Curve These results highlight a limited effect of the differences in the morphological properties of the longitudinal sections in the same specimen, and between different specimens made under the same process conditions, on the estimate of the relative density, and therefore on the theoretical curve of compression behavior. This can be traced back to two aspects:

•
The morphological properties of closed-cell structures of metal foams that affect the mechanical behavior are similar between longitudinal sections, while they change significantly with respect to cross sections [38].
• The process of translating the local morphological parameters of the real sections into geometric parameters of the virtual model of the cellular structure, described in detail in Section 3, attenuates the effects of the local morphological characteristics.
As a result, the proposed model for the correlation between morphology and compression behavior turns out to be robust with respect to the choice of the section for the morphological analysis, and to the variability of morphological properties between different specimens obtained by the same process conditions. The first finding solves the problem of which section to choose, among the countless longitudinal sections that rotate around the loading axis of the same specimen, for the two-dimensional morphological analysis on which the correlation model introduced in Section 3 is based. The second one guarantees the reproducibility of the results as the specimen varies.

Predictive Use of the Fitted Model
The correlation model, once fitted on all compression steps for the specimen with specific properties (0.4% TiH 2 , 2.8% SiC), can be used to predict the behavior of similar cell structures. As reported in Section 2.1, other two types of closed-cell Al foam cylindrical specimens were made and tested, varying the TiH 2 forming agent content in powder composition: 0.1 and 0.2 weight % (with SiC content fixed at 2.8 weight % as for previous specimens). Using only the morphological analysis corresponding to the Step 0 state as input, modeling and predictive simulation of the compression behavior of the two similar foams were performed.
Using the fitting parameters γ p , γ s , λ, obtained previously for 0.4% TiH 2 foam at Step 0 ( Table 2, row 1), the relative density ρ f /ρ s was calculated from the morphological analysis of Step 0 state for both 0.1% and 0.2% TiH 2 foams, obtaining the values 0.294 and 0.251 respectively. These theoretical values are congruent to the corresponding experimental values equal to 0.310 and 0.240.
The following procedure was then performed for predictive simulation of the compression behavior of the two similar foams:

•
The predictive values of the relative densities ρ f /ρ s for the subsequent compression steps were calculated by applying the strain-dependent law determined in the previous case by interpolation (Table 3, column 3), modulated with respect to the value of the relative density at Step 0 of the similar foams.

•
The parameters for the stress-strain curves at each compression step were calculated using the aforementioned Equations (16), (18), (21)- (23), with the predictive values of relative density ρ f /ρ s , and the fitted values in Table 2 for c 1 , c 3 , ξ, ψ.
Finally, the same construction procedure by envelope used in the case of 0.4% TiH 2 foam and described in the Section 3.2 was applied, to obtain the simulated compression curves for 0.1% and 0.2% TiH 2 foams, shown respectively in Figures 15 and 16 in comparison with the corresponding curves obtained experimentally, in analogy with Figure 12.
As a general observation, it is possible to note that in the first case an approximate predictive simulation regarding the elastic and plastic collapse phases, and a rather accurate prediction in the densification regime, are obtained. In the second case, instead, the accuracy of the predictive simulation is higher and more homogeneous.
This result is particularly significant if it is considered that the powder composition to produce the metal foams and the morphological characteristics can be correlated [33]. The predictive potential of the morphology-behavior correlation model proposed here thus becomes an integrative tool of these previous studies, allowing to formulate good approximation predictions on the whole compression behavior of the foams as the powder's composition changes, or vice versa allowing us to proactively identify the precursor composition suitable to obtain a specific mechanical behavior.

Effectiveness of the Model, Limitations, and Field of Application
A quantitative analysis of the key parameters of the compression curves shown in the Figures 12, 15 and 16, and the comparison between their simulated and experimental values, allow to provide a more precise evaluation of the accuracy of the simulations and highlight the improving effect obtained by extending the morphological analysis and the fitting of the model on the compression steps detected during collapse evolution.
For these purposes, the mean plateau stress σ pl and the densification strain ε d were taken into consideration, among the main parameters of the compressive curves. The modulus of elasticity E f has not been considered for the reasons already highlighted in Section 4.2, which exclude the possibility of a theoretical evaluation of Young's modulus comparable to the experimental one. The densification onset ε o , on the other hand, was excluded due to the difficulty of defining univocal values in the experimental curves.
With these premises, Table 6 shows the theoretical values (i.e., the values corresponding to the curves simulated using the proposed model) and experimental values of mean plateau stress σ pl and the densification strain ε d , for the three foams examined, as the TiH 2 forming agent content in powder composition varies. The percentage errors in the theoretical values with respect to the corresponding experimental values are included in brackets. In all cases, the theoretical values of the compression parameters are reported with respect to two conditions:

•
Considering only the starting morphological properties, corresponding to Step 0 pre-compression state (Step 0 curve such as in Figure 10 for 0.4% TiH 2 foam).

•
Constructing the final compression curves according to the procedure defined in Section 4.2, that allowed to obtain the theoretical curves in Figures 12, 15 and 16 (all steps curves, obtained taking into account the evolution of the compression process up to collapse in fitting the model). The main conclusion to highlight by observing the results reported in the table can be summarized as follow:

•
As expected, the most accurate theoretical curve (whose compression parameters have the lowest error values compared to the experimental ones) is the all steps curve corresponding to the 4% TiH 2 foam (Figure 12), which was entirely fitted on the corresponding experimental curve.

•
The two predicted theoretical curves (all steps curves for 1% and 2% TiH 2 foams, in Figures 15 and 16 respectively) show a good accuracy in the estimation of the stress plateau σ pl , and a high precision regarding the simulation of the densification stage (the errors for ε d are very limited).

•
In all cases, the improvement in accuracy due to fit the correlation model between morphology and mechanical behavior according to several compression steps (obtaining the all steps curves), compared to the results obtained by a morphological analysis limited to just one compression state, such as the initial pre-compression one (obtaining Step 0 curves), is clearly represented by the marked reductions in the corresponding errors with respect to the experimental curves (confirming what was anticipated in Section 4.2 regarding the construction of the compression curves).
The same data in Table 6, with particular regard to those obtained for the all steps curves can be used to evaluate the accuracy of the proposed model in comparison with other models for simulating the compression behavior of closed-cell Al foams, which are present in the literature. With this purpose, the comparison was set up according to the following premises:

•
No analytical model in the literature was taken into consideration, due to the limited accuracy in the estimation of the compression parameters [18]; for this reason, the comparison was made with some of the numerical models presented in the introduction (Section 1), and chosen so to diversify the comparison and cover the main modeling types. • Also in this case the compression parameters taken into consideration were the mean plateau stress σ pl and the densification strain ε d ; while for the former it was always possible to extract the data from the literature, regarding the latter it must be pointed out that some models fail to provide a sufficiently clear reconstruction of the densification curve, so the corresponding densification strain cannot be defined (therefore in this cases it is missing in the comparison table).

•
Only the percentage errors between the simulated and experimental values were used as significant indicators for the comparison, because they allow to quantify the accuracy of the models in the simulation of mechanical behavior, regardless of the properties of the cellular structure analyzed; furthermore, since it is clearly evident that the accuracy of the models varies with the relative density of the foam, it was considered appropriate to report this data as well.
With these premises, the results obtained by the proposed analytical model (PAM), reported in previous Table 6 for all steps curves, were compared to the results obtained by the following numerical models from the literature:

•
As representative models based on polyhedra, the cruciform-hemisphere model proposed by Meguid et al. [17] (NM3), and the Hasan's weaker-cruciform-strongerhemisphere model [19] (NM4); the latter was also assumed as reference model to set the preliminary ordered structure of the foam in the present study.

•
As an example of advanced tessellation-based numerical modeling, the irregular Voronoi model implemented by Vengatachalam et al. [30] (NM5). Table 7 collects the values of the percentage errors of the simulations with respect to the experimental data. For some models, the data are diversified by the relative density of the cellular structure and/or by the type of cell shape. As a first consideration, it can be observed that the densification parameter ε d is hardly identified. This is due to a general inaccuracy in the simulation of the densification regime, even in the case of the most sophisticated numerical model (NM5), whose estimated error is indicative. On the contrary, the proposed model PAM allows an accurate simulation of the entire densification curve.
The accuracy of PAM in estimating the mean plateau stress σ pl , excluding the foam used to fit the model and considering only the two predictive simulations (0.1% and 0.2% TiH 2 ), stands at low error levels, as compared to the NM1 model (excluding the case of elliptical cells), and to the uniform density NM3 model, and is similar to the NM5 model. Anyway, it should be noted that the proposed model, while losing precision in the mean value of the plateau stress (if compared to numerical models such as NM2, NM3 with randomized density, and NM4 at 0.200 relative density), simulates in detail the entire plastic plateau, and describes with significant precision the whole compression up to densification (as can be deduced from the qualitative observation of Figures 15 and 16), distinguishing itself from all the numerical models compared.
This result seems to contrast with the simplifying nature of the analytical modeling proposed, in comparison with the higher level of detail obtainable by numerical FE simulation, whose results generally present a great improvement if compared to the analytical predictions [19,21]. As a matter of fact, one of the limitations of the morphology-mechanical behavior correlation model is to gather and represent only some morphological trends of the structure, related to the size and shape differences between the cells, to the thickness differences between the cell walls, to the orientation of the cells with respect to the loading direction. The discussed results, which show a high level of effectiveness in the predictive simulation of the whole stress-strain curve, demonstrate that the proposed analytical model compensates for the approximations and the loss of detail by mean of the accurate fitting on the compression steps analyzed during the evolution of the structure collapse.
Another characteristic of the model, which can represent a limitation of its application, is the anisotropy, due to the fact that it is developed with respect to the load direction. Examining whether a cellular structure is geometrically anisotropic is a critical issue [49]. In this respect, it is known that the effects of the structural anisotropy can be considered by performing compressive tests under different loading directions [50]. This peculiarity of the proposed model does not really differ from the behavior of the type of cellular structures examined here, which are obtained by a highly directional foaming process [51]. Anyway, here the anisotropy of the model is functional to more accurately investigate the compressive behavior along the specimen axis, without taking into account the properties along the radial directions.
Finally, as last remarks on the applicative potentialities, it should be pointed out that the proposed model allows to simulate highly inhomogeneous cells, as shown by the investigation carried out on specimens produced by means of the compacted powder method. Since the model parameters are calculated by fitting on the experimental compressive curve, nothing suggests that the model cannot also be successfully used in the case of more ordered and morphologically homogeneous structures; on the contrary, these structures would involve a minor loss of detail in the geometric modeling.
Among the homogeneous cellular structures, the use of the proposed model in the case of syntactic foams should be excluded, since the simulation of the mechanical behavior of the microspheres that characterize their structure requires specially developed mechanical models [52].

Conclusions
The comparison between the experimental stress-strain curve and the one simulated by the theoretical values of the morphological parameters at the initial state of compression highlights the need to fit the correlation model between morphology and mechanical behavior according to several compression steps. In fact, the morphological analysis limited to just one compression state, would only allow for defining a strictly constant value of the plateau stress, resulting in an overall simulation far from the real behavior, characterized by a plateau stress not well stable, and tending to rise with increasing strain.
Therefore, for one of the produced foams (0.4% TiH 2 , 2.8% SiC), the analysis of the evolution of the cells' morphology during deformation up to the final densification was carried out at different levels of compression by X-ray CT observations. The correlation model was fitted with regard to the morphology of each compression step, and the simulated stress-strain curve was obtained by enveloping the curves obtained for each step. When compared to the experimental one, a good level of approximation of the whole curve was found, with a realistic description of the behavior up to the final densification. In particular, a full congruence between the theoretical and experimental values of the densification strain, with a limited deviation in stress values of the densification curve, was achieved.
The capability of the model (fitted on the morphological evolution during compression for 0.4% TiH 2 foam) to predict the mechanical behavior of other similar foams, by means of the unique morphological analysis of the condition prior to compression, was validated. The comparison between simulated and corresponding experimental stress-strain curves for 0.1% and 0.2% TiH 2 foams, revealed in both cases rather accurate prediction of the whole compression curve, particularly in the densification regime.
Taking into account that previous studies demonstrated how the initial composition of the powder and the morphological characteristics can be correlated, the predictive potential of the morphology-behavior correlation model, proposed here, opens new perspectives on the design of cellular structures with the preliminarily simulated compressive behavior.