Numerical prediction of orthotropic elastic properties of 3D-printed materials using micro-CT and representative volume element

• Published in: Acta mechanica 2020-02, Vol.231 (2), p.503-516
• Main Authors: Biswas, P., Guessasma, S., Li, J.
• Format: Article
• Language: English
• Subjects: Angles (geometry); Anisotropy; Boundary conditions; Classical and Continuum Physics; Computed tomography; Computer-aided design; Condensed Matter; Control; CT imaging; Diagnostic imaging; Dynamical Systems; Elastic anisotropy; Elastic properties; Engineering; Engineering Fluid Dynamics; Engineering Thermodynamics; Filaments; Finite element method; Heat and Mass Transfer; Mathematical models; Mechanical properties; Micromechanics; Numerical prediction; Original Paper; Physics; Porosity; Printed materials; Process parameters; Solid Mechanics; Statistical Mechanics; Stiffness matrix; Theoretical and Applied Mechanics; Three dimensional models; Three dimensional printing; Vibration
• ISSN: 0001-5970; 1619-6937
• DOI: 10.1007/s00707-019-02544-2

The mechanical property of 3D-printed components often exhibits anisotropic behaviors and a strong dependence on printing orientations and process parameters. In this study, computational models based on microstructures of 3D-printed ABS polymers are developed using micromechanics of a representative volume element (RVE) to investigate the orthotropic elastic properties. Two finite element models, based on Micro-CT or CAD geometry, with different raster angles: 0 ∘ , 30 ∘ , 45 ∘ , and 60 ∘ (corresponding to the layups of 0 ∘ / 90 ∘ , 30 ∘ / - 60 ∘ , 45 ∘ / - 45 ∘ , and 60 ∘ / - 30 ∘ ), are considered. The Micro-CT model useds the realistic geometry of a 3D-printed cube reconstructed from Micro-CT scans. The CAD model useds periodic layers and filaments with sizes specified according to the dimensional statistics from the Micro-CT model, including bond width, layer height, filament width, and total porosity. All models are subjected to six independent load cases of macroscopically uniform boundary conditions (BCs) admitted by the Hill–Mandel condition to obtain a fully orthotropic elastic stiffness matrix. Scale-dependent bounds from those BCs are used to establish the convergence of RVE responses. The theoretical method of Mori–Tanaka homogenization is also employed to verify numerical predictions. The results of both Micro-CT and CAD models are close to each other and agree well with experimental results in the literature. Parametric studies are further performed in CAD models by varying layer height, filament width, and bond width to investigate their effects on the orthotropic elastic properties. It is found that the porosity plays a significant role, and more porosity leads to larger elastic anisotropy in 3D-printed materials.