The main motivation behind this study is to develop a methodology to generate surrogate model coefficients for parametrized Spheroidal Graphite Iron (SGI) for thermal expansion and Young’s modulus response. To accomplish this objective, a 3D model was created by Abaqus, and Python was used for developing the script. The design of experiments (DoE) was utilized to compute two objectives parametrically. The inputs in this thesis were Poisson’s ratio, Young's modulus, radius, and volume fraction of Graphite, Ferrite, and Pearlite. Then DoE is integrated with Abaqus to get homogenized properties of the SGI in terms of Coefficient of thermal expansion (CTE) and Young’s modulus. Using these homogenized properties, the surrogate model coefficients were obtained, and Young’s modulus-CTE graph was visualized.