Finite Element Analysis of Stress Evolution in Al-Si Alloy Sudha Joseph and S. Kumar (Submitted July 10, 2014; in revised form September 28, 2014; published online November 7, 2014) A 2D multi-particle model is carried out to understand the effect of microstructural variations and loading conditions on the stress evolution in Al-Si alloy under compression. A total of six parameters are varied to create 26 idealized microstructures: particle size, shape, orientation, matrix temper, strain rate, and tem- perature. The effect of these parameters is investigated to understand the fracture of Si particles and the yielding of Al matrix. The Si particles are modeled as a linear elastic solid and the Al matrix is modeled as an elasto-plastic solid. The results of the study demonstrate that the increase in particle size decreases the yield strength of the alloy. The particles with high aspect ratio and oriented at 0° and 90° to the loading axis show higher stress values. This implies that the particle shape and orientation are dominant factors in controlling particle fracture. The heat treatment of the alloy is found to increase the stress levels of both particles and matrix. Stress calculations also show that higher particle fracture and matrix yielding is expected at higher strain rate deformation. Particle fracture decreases with increase in temperature and the Al matrix plays an important role in controlling the properties of the alloy at higher temperatures. Further, this strain rate and temperature dependence is more pronounced in the heat-treated microstructure. These predictions are consistent with the experimentally observed Si particle fracture in real microstructure. Keywords Al-Si alloy, heat treatment, finite element modeling, fracture, Si particle morphology, stress analysis 1. Introduction The ductility of Al-Si-based cast alloy mainly depends on the fracture characteristics of eutectic Si particles, since the primary source of damage is cracking of these particles at relatively low values of strain (Ref 1). Under monotonic loading, the fracture and debonding of Si particles leads to the nucleation, growth, and coalescence of voids (Ref 1, 2). Consequently, the failure of Al-Si cast alloy is very sensitive to the Si particle morphology. In general, coarse and irregular shaped particles promote rapid damage evolution and result in a cast Al alloy with low ductility and strength. In addition, the orientation of the particles with the loading axis and the heat treatment given to the alloy also play an important role. Thus, in order to predict the failure, it is necessary to understand the stress states of Si particles and Al matrix. It is well known that the properties of the alloy can be enhanced by reducing the size and aspect ratio of the Si particles, by increasing the cooling rate and/or modifying by Sr addition. It has been reported that the size of the largest Si particles controls the ductility of the alloy (Ref 3). It has also been suggested that the ductility of unmodified alloys is controlled by the mean size of the Si particles whilst that of modified alloy is determined by their distribution (Ref 4). McLellan (Ref 5) reported that ductility is determined by Si particle aspect ratio. It is also reported that the heat-treated Al-Si-based alloys exhibit lower ductility though strength is enhanced (Ref 6, 7). The detailed microstructural examinations reveal that the particles with high aspect ratios consistently show an increased propensity for fracture and the particlesÕ orientation also has strong effect on particle fracture (Ref 8–11). We have carried out a comprehensive experimental investigation to understand the particle fracture characteristics in a near-eutectic Al-Si- based cast alloy with different microstructures under compres- sion and the results have been published elsewhere (Ref 12, 13). It was observed that Si particles with higher aspect ratio and oriented nearly perpendicular to the loading axis fracture more. The heat treatment of the alloy also increases particle fracture. Further, Si particle fracture is found to increase with increase in strain and strain rate and decrease with increase in temperature. All these effects on particle fracture can be explained clearly if the stress state in Si particles is understood, which is possible only by numerical calculations using materials models. There are various models to explain the particle fracture and they are mainly based on dislocation-based theory or contin- uum theory. Yeh and Liu (Ref 8) proposed a micromechanical model based on dislocation pile-ups to explain the particle fracture of the Al-Si-based alloy. The authors concluded that the particle fracture occurs by the impingement of dislocations on the particles. However, their model does not include the influence of geometrical aspects, such as aspect ratio and orientation, of the particles on fracture. These effects can be captured well by models based on continuum theory. In the present manuscript, we have used the continuum based multi-particle model, in which the stress state of Si particles and matrix is quantified by varying parameters such as, particle size, aspect ratio, and orientation. The stress evolution in Al-Si alloys with different microstructures under fatigue loading was carried out by Gall et al. (Ref 14) using Sudha Joseph and S. Kumar, Department of Materials Engineering, Indian Institute of Science, Bangalore 560 012, India. Contact e-mails: sudhajos.iisc@gmail.com and skumar@materials.iisc.ernet.in. JMEPEG (2015) 24:253–260 ÓASM International DOI: 10.1007/s11665-014-1288-x 1059-9495/$19.00 Journal of Materials Engineering and Performance Volume 24(1) January 2015—253