Constitutive modeling of Solargrade silicon materials

Residual stresses in directionally solidified silicon ingots are known to increase the likelihood of material failure, leading to unnecessary costs and material waste. These stresses develop during solidification and cooling, owing to thermal gradients, the multicrystalline nature of the materials, and its impurity content. Silicon multicrystals (mc-Si) are simply aggregates of several, rather large single crystals. Hence, building a model for mc-Si consists in the development of accurate constitutive equations for SoG-Si monocrystals and the construction of multicrystalline aggregates by Finite Element methods. Constitutive models for solar-grade silicon monocrystals (SoG c-Si) of various complexity are derived in this thesis. Their ultimate goal is to allow the quantification of the residual stresses in SoG-Si materials given a thermomechanical loading path. A physical approach to modeling plasticity in semiconductors has been consistently adopted: all equations are based on the introduction of dislocation densities, carrying plastic flow and participating to material hardening with deformation. Constitutive equations are implemented into a rate-dependent crystal plasticity kinematical framework. The benefits of this methodology are twofold: firstly, the models are perfectly useable in the temperature range where silicon behaves as a perfectly elastic material, but they excel at reproducing its particularly temperature-sensitive plastic behavior. Secondly, final dislocation densities can be predicted by the models. These are of great use for industrial applications, since regions of large dislocation contents also exhibit poor electrical performance. Consequently, the constitutive models introduced in this Ph.D. thesis are not limited to stress calculations, but can be coupled to other physical analyses. An exhaustive literature review has been written, covering the aspects of plastic deformation over multiple scales. It lays the foundation for analysis and improvement of existing constitutive models. Systematic identification of constitutive parameters of traditional models calls for the introduction of new internal variables, dislocations stored in dipolar structures. Accounting for these immobile dislocations enables an extension of the model of Alexander & Haasen to correctly reproduce the steady-state of deformation during uniaxial testing of silicon monocrystals oriented for single glide. The influence of dissolved oxygen on the plastic behavior of