# Finite Element Analysis of Temperature Distribution and Stress Behavior of Squeeze Pressure Composites

Authors: P. Gurusamy; T. Sathish; V. Mohanavel; Alagar Karthick; M. Ravichandran; Omaima Nasif; Saleh Alfarraj; Velu Manikandan; S. Prasath
Journal: Advances in Materials Science and Engineering
 (2021)

---

## Abstract

Aluminium-reinforced composites play a vital role in the engineering industry because of their better strength and stiffness. The properties are directly related to the solidification phenomenon of the cast alloy. The design engineer should understand the importance of the solidification behavior of base alloy and its reinforcement. Composites’ solidification study is rare, and the reviews are limited. The solidification process is analyzed using the finite element method (FEM), and this would fetch a lot of information about the cooling rate of the composites and also helps to reduce the time in experimentation. This paper reports and plots the cooling curves of Al/SiCp composites using simulation software. Cylindrical-shaped composites were developed using the squeeze casting method, and the experimental cooling curves were plotted using a K-type thermocouple. Composites samples were prepared at the following squeeze pressures: 0, 30, 50, 70, 100, and 130 MPa; melt and die temperature was kept constant at 800 and 400°C, respectively. The experimental and FEA cooling curves were compared, and it was agreed that the increase in the squeeze pressure increases the cooling rate of the developed composite. Furthermore, the effect of temperature distribution from the inner region of the melt and die material which causes the radial and tangential stress of components has also been examined.


---

## Body

## 1. Introduction

Aluminium alloys, reinforced with SiC, have become an essential material that is used in a wide range of applications in the aviation, marine, electronic packaging, defence, and automobile industries [1]. The nature of the matrix metal is generally soft which contains good mechanical and thermal properties. The most common metals are aluminium, copper, titanium, magnesium, etc. The standard properties offered by the metals are malleability, ductility, fatigue strength, machinability, and thermal and electrical conductivity. The reinforcing material has to offer strong bonding to the metal in order to make the composite tensile, and the stiffness to weight ratio and strength to weight ratio is higher than the metal matrices. The hard and robust nature is obtained when the correct amount of the reinforcements is added to the matrix. The particle-reinforced metal matrix composites exhibit better physical and mechanical properties; it can be modified to reach metallurgist requirements. Also, this material can show excellent thermophysical properties due to which it has been used in all the structural applications. It can pose a tremendous amount of strength at high-temperature applications; the capability to reveal low thermal shock and excellent wear resistance at less weight is the unique characteristics. There are many manufacturing methods to produce MMCs on which solid- and liquid-state processing is the promising route for developing the composites [2]. In solid-state processing, the interfacial reaction between base alloy and ceramic particles is not right in agreement. Hence, liquid-state processing can resolve this issue through stir casting. Further improving the better properties of the composites, the process adopts the squeeze casting method. It is a method to remove shrinkage and porosity during the liquid-state processing of composites. Hence, near net shape components are fabricated by this method [3].Ghomashchi and Vikhrov [4] reported that the essential parameters for fabricating any alloys using squeeze casting routes are melt temperature, die temperature, and application of pressure, the main premise for this research. The base alloy’s temperature, heat transfer coefficient, thermal conductivity, die temperature, and pressure during solidification would change the phase relationship of any alloy system. Thus, the application of pressure helps to decrease the interatomic movement during the process of solidification, which allows changes in the temperature of the melt. Many other factors would influence increase of the pressure, which are impurities, trace elements, and intersolubility of constituents’ elements. Yue and Chadwick [5] have observed that the application of pressure would reduce the gap between the mould and melt, which leads to a more significant cooling rate. It influences the cooling rate. In permanent mould casting, due to the surface contact between the mould and casting, the cooling rate and solidification time must be considered for better results [6]. The mechanical properties of the casting are derived from the solidification time of the casting. Sukumaran et al. [7] investigated the squeeze cast A2124 alloy with 10% ceramic particles with increasing pressure during solidification of alloy and the developed composite. He reported that the application of pressure directly relates to the enhancement of composites. The applied pressure of about 100 MPa has created uniform particle distribution and other factors such as porosity and shrinkages and obtained fine refinement of microstructure. Furthermore, he reported that the squeeze pressure coupled with die temperature is increased; the heat transfer rate in composites is higher compared with that of the base alloy. The heat transfer rate is directly related to the amount of particles present in the alloys, which decides the cooling rate of the composite. Hanumanth and Irons [8] have reported that particle distribution can increase the cooling rate up to a certain percentage of reinforcing particles; after that, the cooling rate is decreased. Increasing the particles reduces the heat transfer rate in the melt, which causes the thermal diffusivity of the composite to reduce. Studies by Rajan et al. [9] have shown varying volume fractions of graphite and SiC particles with AA6082 alloy in the sand, steel, and graphite moulds, and heat transfer at casting and mould interface is also investigated. The heat transfer rate is enhanced due to the addition of one percent of magnesium particles in the aluminium alloy; this increases the wettability between the particle and the base alloy. The magnesium particles can be useful for hunting the oxide layer as well as increasing the interfacial reaction between the melt and the mould. Many researchers found, during the study of interfacial reaction, that the effect of intermetallic compounds plays a vital role. Zhou [10] reveals that the intermetallic compounds (Al4C3) would decrease when the solidification process is disturbed, which means it affects the fluidity of the materials. The compunds are varied with respect to the base alloy. Studies on Al and Graphene-bonding strength and wettability literature are not enough so far, these two intermetallic compounds, such as Al4C3 and Al2O3, are helping the base alloy in terms of transferring loads [10–12]. Li et al. [13] claimed that no evidence had noticed the formation of Al4C3 between Al and graphene nanoplatelets (GNPs). Nonetheless, stacking faults were observed near the Al matrix, which could be the reason for the strengthening of composites. Incorporating metal reinforcement (Cu-GNPs and Ni-GNPs) along with GNPs would enhance the strength of the composite; perhaps, the metal particles could have stopped the formation of Al4C3 [14]. Liu et al. [15] reported that the application of pressure would enhance the interface strength of Al/Mg-based composites. The intermetallic compounds are directly related to the cooling rate of the composites, and the application of squeeze pressure enhances the cooling rate of the composite [16]. The effect of intermetallic compound (Al4C3) is modeled with the help of Ansys software, and the impact of that is compared with the rule of mixture predictions. No significant results have been found in the Al-Si/CNT composites [17].There are many studies based on fracture behaviour, interpenetrating with some structure models, simulation of moulds [18, 19], addition of particles with base alloy [20–22], and mechanical properties of particles [23–25] and using the finite element method. The study based on the solidification and cooling rate of the composites is highly complicated. Because, during the solidification process, the complex shapes’ grain refinement and the features of microstructures are changed. Also, there are many studies which are focused only on the base alloy cooling rate and solidification process. This study gives importance of two different materials’ solidification process because of its different behaviours. Hence, the squeeze casting process would give better results and integrate the cooling rate and microstructural features. There is a lack of information regarding the application pressure on the solidification study of cast composites, which is an essential feature in the high-pressure die casting industry. Therefore, the present work focuses on both experimental and theoretical research on the solidification of AA6082 composites with temperature behaviour and stress behaviour using the finite element method.


## 2. Boundary Conditions for the Cylindrical Model

The cooling curves have been plotted by the simulation software-ANSYS (version 19.2). Structural analysis, followed by thermal analysis, is carried out by the coupled field analysis. Hence, the structural and thermal capability of the element (PLANE 13) is used. Thus, each node has up to four-degree degrees of freedom. In this study, first, nonlinear thermal analysis was carried out, and then, static stress analysis was carried out. The half portion of a cylinder with the melt model was plotted for this symmetrical problem. Figure1 shows the two-dimensional model of casting with the melt with its dimensions. The thermophysical properties of ceramic materials and AA6082 are shown in Table 1. The processing conditions for the preparation of cast samples are listed in Table 2. Tables 3 and 4 show die and ceramic material properties [9, 26]. Table 5 shows the boundary conditions of the FEA model. Convection properties were initially given to the model, and the initial room temperature is also applied to the die wall. Figure 2(a) shows the model for melt and die without boundary conditions, and Figure 2(b) shows the convection boundary condition for the die material. The rule of mixture was used to calculate the thermal conductivity of the composite. In contrast, the value of enthalpy is considered as the same value of base alloy because there is no phase change or ceramic material [22]. A step-by-step process of the solidification process concerning time is established by solution control. Finally, the cooling curve is plotted with the help of a time-history processor.Figure 1
The cylindrical model with BCs for the die temperature of 400°C (melt temperature 800°C and pressure 130 MPa).Table 1
Chemical composition of AA6082 aluminium alloy.
ElementSiFeCuMnMgCrZnTiAlWeight (%)1.00.600.100.50.70.250.100.20BalTable 2
Process conditions for composites.
Die temperature (°C)Matrix temperature (°C)Pressure (MPa)4008000, 30, 50, 70, 100, 130Table 3
Properties of SiC [9].
MaterialThermal conductivity (W/m°C)Density (kg/m3)Specific heat (J/kg·k)Silicon carbide10032001300Table 4
Properties of H13 [26].
MaterialThermal conductivity (W/m°C)Density (kg/m3)Specific heat (J/kg·k)AISI H1328.67761461Table 5
Boundary conditions for the applied pressure of 130 MPa.
Key pointsVariables (mm)Parameters (°C)1x = 0; y = 0Td = 4002x = 60; y = 0Td = 4003x = 60; y = 175Td = 4004x = 25; y = 175Tm = 8005x = 25; y = 25Tm = 8006x = 0; y = 25Tm = 8007x = 0; y = 175Tm = 800Figure 2
(a) FEA model of melt and die without BCs. (b) Convection boundary conditions of the FEA model.
(a)(b)


## 3. Experimental Work

The base plate width was 120 mm, the die inner diameter was 50 mm, thickness was 35 mm, and height was 150 mm. Figure3(a) shows the melt stirring with reinforcement setup, and Figure 3(b) shows the punch and die have been preheated. Figure 4(a) shows the baseplate and die, and Figure 4(b) shows the appropriate locations of thermocouples. Figure 4(c) shows the fabrication of the squeeze cast setup. The squeeze casting technique was adopted for fabricating the composites, as discussed in [3]. The temperature values were being recorded by K-type thermocouple, which was recorded every one second. Figure 4(d) shows the thermocouples in its exact locations. The data acquisition system was incorporated in the experimental setup to record all temperature values, precisely. The recorded values were transferred to the computer to draw solidification curves. Aluminium composites were fabricated with increasing pressures when the die and melt temperatures were 400 and 800°C, respectively.Figure 3
(a) Melt stirring setup. (b) Die and punch preheating for the preparation of composite.Figure 4
(a) Die and baseplate. (b) Appropriate locations of thermocouples for the complete setup. (c) Experimental squeeze cast setup. (d) Exact locations of thermocouples in the baseplate and die material.
(a)(b)(c)(d)


## 4. Results and Discussion

### 4.1. Cooling Rate and Time

The temperature readings measured from the baseplate readings (T0 to T2) are used to plot cooling curves and solidification time. In contrast, the die materials’ temperature readings (T3 to T5) are used to measure the temperature distribution of the die. The temperature readings are used to determine the solidification time for any type of cast alloy [26].Increasing the pressure decreases the solidification time from 60 seconds to 42 seconds [2]. Figure 5 shows the various squeeze pressures cooling curves for both experiment and FEA study of Al/SiCp composite. The cooling rate is observed from the experimental cooling curves, and the values are 2.28°C/s 2.52°C/s, 2.68°C/s, 2.87°C/s, 3.5°C/s, and 3.75°C/s for 0, 30, 50, 70, 100, and 130 MPa, respectively. The flattening region is the freezing zone, and it is clearly shown in Figure 5.Figure 5
Comparison of experimental and FEA curves for Al AA6082/SiCp composites.The melting range of AA6082 alloy in the solid state is 577°C and the liquid state is 615°C, and the experimental freezing range is closely agreed with that of Hanumanth and Irons [8]. Also, in squeeze casting, the application of pressure brings the melt to solidify first, and then, the cooling rate is decided by the location of the particle [27]. Maleki et al. [28] reported that the application of pressure enhances the die and melt contact surface, which promotes high heat transfer and cooling rate. Hence, increasing the pressure increases the cooling rate of the composite. The predicted cooling rates are 1.23°C/s, 1.31°C/s, 1.47°C/s, 1.81°C/s, 1.98°C/s, and 2.24°C/s for the squeeze pressure of 0, 30, 50, 70, 100, and 130 MPa, respectively. Particle and base alloy reaction, formation of air inside the mould and melt, aspirations, and flow velocity were not considered in this analysis. The mentioned factors have to be considered for the FEA study [22], even though the cooling rate is almost close to the experimental cooling rate.


### 4.2. Temperature Distribution and Thermal Stresses

To avoid multiple trial-and-error casting production or to avoid recreating the system repeatedly, the FEA-based analysis is recommended to fulfill the technical problems [13, 14]. Temperature distribution within the melt and distribution in the die are essential factors to identify the rate of the solidification process, mainly when the melt is in the application of pressure. Figure 6 shows the temperature dispersion in the composites at various squeeze pressures up to 800 seconds. From the FEA analysis, the observed values are 314°C, 318°C, 319°C, 321°C, 322°C, and 323°C for different squeeze pressures 0, 30, 50, 70, 100, and 130 MPa, respectively. The amount of heat transfer is observed in the melt and reinforcement based on the applied pressure [29]. Hence, the temperature is varied based on the applied pressure.Figure 6
Temperature dispersion in the composites at various squeeze pressures up to 800 seconds: (a) 0 MPa; (b) 30 MPa; (c) 50 MPa; (d) 70 MPa; (e) 100 MPa; (f) 130 MPa.
(a)(b)(c)(d)(e)(f)It is important to focus complex shape engineering application-oriented castings, using design and solidification process using the finite element method. FEA is significant for casting industry because components’ rejection leads to increase in cost of production. The finite element method is better in analyzing heat transfer of components within the system and outside the system. In this analysis, the initial thermal stresses are kept at 0°C, and it was considered that the composites’ mechanical property values are taken as per the experimental work because the mechanical property values may change at higher temperature. With these constraints, the temperature (linear) distribution and tangential and radial stresses are plotted for the squeeze-cast specimen. The melt and reinforcements have different thermal expansions; hence, the study of radial and tangential stresses is essential in this analysis. Figure7 shows the radial and tangential stresses obtained with increasing squeeze pressures. The application of pressure is in the axial direction; hence, the radial stresses are high at the bottom of the die. At the melt and die interface location, the tangential stress is maximum. Table 6 shows the stress dispersion of the squeeze cast composite for 800 seconds. Tangential stresses are comparatively higher than the radial stresses, and in the interfacial region, the tangential stresses are decreased [30]. It is due to the application pressure, which leads to a thermal mismatch within the composite [31]. Hence, the stresses are closely related to the orthotropic behavior of composite material.Figure 7
The radial and tangential stresses obtained with increasing squeeze pressures: (a) 0 MPa; (b) 30 MPa; (c) 50 MPa; (d) 70 MPa; (e) 100 MPa; (f) 130 MPa.
(a)(b)(c)(d)(e)(f)Table 6
Stress dispersion of the squeeze cast composite for 800 seconds.
Squeeze pressure (MPa)Radial stress (MPa)Tangential stress (MPa)MaxMinMaxMin0252−468639−47830256−475635−48350257−477629−48670258−479623−488100259−481613−490130260−483603−492


## 4.1. Cooling Rate and Time

The temperature readings measured from the baseplate readings (T0 to T2) are used to plot cooling curves and solidification time. In contrast, the die materials’ temperature readings (T3 to T5) are used to measure the temperature distribution of the die. The temperature readings are used to determine the solidification time for any type of cast alloy [26].Increasing the pressure decreases the solidification time from 60 seconds to 42 seconds [2]. Figure 5 shows the various squeeze pressures cooling curves for both experiment and FEA study of Al/SiCp composite. The cooling rate is observed from the experimental cooling curves, and the values are 2.28°C/s 2.52°C/s, 2.68°C/s, 2.87°C/s, 3.5°C/s, and 3.75°C/s for 0, 30, 50, 70, 100, and 130 MPa, respectively. The flattening region is the freezing zone, and it is clearly shown in Figure 5.Figure 5
Comparison of experimental and FEA curves for Al AA6082/SiCp composites.The melting range of AA6082 alloy in the solid state is 577°C and the liquid state is 615°C, and the experimental freezing range is closely agreed with that of Hanumanth and Irons [8]. Also, in squeeze casting, the application of pressure brings the melt to solidify first, and then, the cooling rate is decided by the location of the particle [27]. Maleki et al. [28] reported that the application of pressure enhances the die and melt contact surface, which promotes high heat transfer and cooling rate. Hence, increasing the pressure increases the cooling rate of the composite. The predicted cooling rates are 1.23°C/s, 1.31°C/s, 1.47°C/s, 1.81°C/s, 1.98°C/s, and 2.24°C/s for the squeeze pressure of 0, 30, 50, 70, 100, and 130 MPa, respectively. Particle and base alloy reaction, formation of air inside the mould and melt, aspirations, and flow velocity were not considered in this analysis. The mentioned factors have to be considered for the FEA study [22], even though the cooling rate is almost close to the experimental cooling rate.


## 4.2. Temperature Distribution and Thermal Stresses

To avoid multiple trial-and-error casting production or to avoid recreating the system repeatedly, the FEA-based analysis is recommended to fulfill the technical problems [13, 14]. Temperature distribution within the melt and distribution in the die are essential factors to identify the rate of the solidification process, mainly when the melt is in the application of pressure. Figure 6 shows the temperature dispersion in the composites at various squeeze pressures up to 800 seconds. From the FEA analysis, the observed values are 314°C, 318°C, 319°C, 321°C, 322°C, and 323°C for different squeeze pressures 0, 30, 50, 70, 100, and 130 MPa, respectively. The amount of heat transfer is observed in the melt and reinforcement based on the applied pressure [29]. Hence, the temperature is varied based on the applied pressure.Figure 6
Temperature dispersion in the composites at various squeeze pressures up to 800 seconds: (a) 0 MPa; (b) 30 MPa; (c) 50 MPa; (d) 70 MPa; (e) 100 MPa; (f) 130 MPa.
(a)(b)(c)(d)(e)(f)It is important to focus complex shape engineering application-oriented castings, using design and solidification process using the finite element method. FEA is significant for casting industry because components’ rejection leads to increase in cost of production. The finite element method is better in analyzing heat transfer of components within the system and outside the system. In this analysis, the initial thermal stresses are kept at 0°C, and it was considered that the composites’ mechanical property values are taken as per the experimental work because the mechanical property values may change at higher temperature. With these constraints, the temperature (linear) distribution and tangential and radial stresses are plotted for the squeeze-cast specimen. The melt and reinforcements have different thermal expansions; hence, the study of radial and tangential stresses is essential in this analysis. Figure7 shows the radial and tangential stresses obtained with increasing squeeze pressures. The application of pressure is in the axial direction; hence, the radial stresses are high at the bottom of the die. At the melt and die interface location, the tangential stress is maximum. Table 6 shows the stress dispersion of the squeeze cast composite for 800 seconds. Tangential stresses are comparatively higher than the radial stresses, and in the interfacial region, the tangential stresses are decreased [30]. It is due to the application pressure, which leads to a thermal mismatch within the composite [31]. Hence, the stresses are closely related to the orthotropic behavior of composite material.Figure 7
The radial and tangential stresses obtained with increasing squeeze pressures: (a) 0 MPa; (b) 30 MPa; (c) 50 MPa; (d) 70 MPa; (e) 100 MPa; (f) 130 MPa.
(a)(b)(c)(d)(e)(f)Table 6
Stress dispersion of the squeeze cast composite for 800 seconds.
Squeeze pressure (MPa)Radial stress (MPa)Tangential stress (MPa)MaxMinMaxMin0252−468639−47830256−475635−48350257−477629−48670258−479623−488100259−481613−490130260−483603−492


## 5. Conclusion

The following observations are identified from the experimental and FEA analysis of an aluminium-based metal matrix composite.(1)
Experimental cooling rates are increased by increasing the squeeze pressure. The effect of pressure enhances the heat transfer rate of the melt and the die material.(2)
Predicted cooling rates are increased by increasing the pressure, but the values are almost close to experimental cooling rates.(3)
The amount of heat transfer is observed in the melt and reinforcement based on the applied pressure. Hence, the temperature is varied based on the applied pressure.(4)
Tangential stresses are comparatively higher than the radial stresses, and in the interfacial region, the tangential stresses are decreased. It is due to the application pressure, which leads to a thermal mismatch within the composite. Hence, the stresses are closely related to the orthotropic behaviour of composite material.(5)
Particle and base alloy reaction, the formation of air inside the mould, and flow velocity have to be considered for the FEA study.



---

*Source: 8665674-2021-09-27.xml*