Numerical Homogenization of Orthotropic Functionally Graded Periodic Cellular Materials: Method Development and Implementation

This study advances the state of the art by computing the macroscopic elastic properties of 2D periodic functionally graded microcellular materials, incorporating both isotropic and orthotropic solid phases, as seen in additively manufactured components. This is achieved through numerical homogenization and several novel MATLAB implementations (known in this study as <i>Cellular_Solid</i>, <i>Homogenize_test</i>, <i>homogenize_ortho</i>, and <i>Homogenize_test_ortho_principal</i>). The developed codes in the current work treat each cell as a material point, compute the corresponding cell elasticity tensor using numerical homogenization, and assign it to that specific point. This is conducted based on the principle of scale separation, which is a fundamental concept in homogenization theory. Then, by deriving a fit function that maps the entire material domain, the homogenized material properties are predicted at any desired point. It is shown that this method is very capable of capturing the effects of orthotropy during the solid phase of the material and that it effectively accounts for the influence of void geometry on the macroscopic anisotropies, since the obtained elasticity tensor has different E1 and E2 values. Also, it is revealed that the complexity of the void patterns and the intensity of the void size changes from one cell to another can significantly affect the overall error in terms of the predicted material properties. As the stochasticity in the void sizes increases, the error also tends to increase, since it becomes more challenging to interpolate the data accurately. Therefore, utilizing advanced computational techniques, such as more sophisticated fitting methods like the Fourier series, and implementing machine learning algorithms can significantly improve the overall accuracy of the results. Furthermore, the developed codes can easily be extended to accommodate the homogenization of composite materials incorporating multiple orthotropic phases. This implementation is limited to periodic void distributions and currently supports circular, rectangular, square, and hexagonal void shapes.