Optimization of Precision Investment Casting for High-Temperature Alloy Components via Numerical Simulation and Robust Design of Experiments

The manufacturing of thin-walled, high-temperature alloy components with intricate internal cavities for aerospace applications presents a significant challenge. Traditional methods such as mechanical machining or forging are often inadequate for forming these complex geometries. Consequently, precision investment casting has become the predominant manufacturing route. This technique offers distinct advantages for the rapid forming of complex parts with high dimensional accuracy, excellent surface finish, and precise elemental control, ultimately leading to castings with superior high-temperature performance. Furthermore, it enhances material utilization and reduces overall production costs.

The precision investment casting process is extensively employed in producing complex thin-walled components for aircraft engines and gas turbines. The use of integrally cast, large-scale intricate parts drastically reduces machining requirements, lowers costs, and improves both service performance and time-between-overhauls, thereby enabling more advanced designs.

Numerical simulation technology stands as the optimal pathway for designing and validating casting processes. With advancements in computational power, simulations that once took months can now be completed in a week or less. This rapid development has led to the widespread application of simulation techniques, including in precision investment casting. Researchers have utilized these tools to optimize parameters like pouring temperature and gating system design for alloys like K4169, achieving refined microstructures and reduced defects. Increasingly, methodologies like the Taguchi method are being combined with high-throughput numerical simulations to efficiently optimize critical process parameters and enhance casting performance.

K4169, a Fe-Ni based superalloy similar to IN718C, is strengthened primarily by the γ” phase and secondarily by the γ’ phase. It exhibits excellent comprehensive properties at temperatures up to 650°C, including high tensile and yield strength, as well as superior stress rupture strength and plasticity. This alloy is widely used for integrally cast casings and associated components in aero-engines. However, research on the pouring parameters for complex thin-walled precision investment casting of K4169 is still evolving. This article employs ProCAST software to simulate the precision investment casting process of an aerospace K4169 component. It analyzes filling, solidification, and defect formation, and integrates the Taguchi experimental design method to optimize process parameters with the goal of minimizing shrinkage porosity, thereby establishing a robust and high-quality casting process.

Component, Methodology, and Initial Analysis

The subject is an annular disk-class aerospace component made from K4169 superalloy, demanding high surface finish, dimensional precision, and freedom from defects like shrinkage cavities, porosity, or cracks. The component features a complex structure with significant variation in wall thickness, regularly arranged peripheral holes, and a narrowed upper section, making it prone to mistrun, gas entrapment, and shrinkage defects during casting. The relatively poor fluidity of molten K4169 further exacerbates the risk of defects, with shrinkage being particularly critical.

The initial casting design, designated Scheme L, employed a traditional bottom-gating system. Numerical simulation under baseline parameters (Pouring Temperature: 1450°C, Shell Preheat Temperature: 950°C, Pouring Time: 5s) revealed a significant volume of shrinkage porosity concentrated in the upper region of the casting. Analysis identified two root causes: 1) The overly vertical sprue and undersized, few ingates led to turbulent filling, potential oxide inclusion, and insufficient metal feed. 2) The upper section, being the last to fill and solidify, suffered from inadequate metal feeding and poor venting, leading to defect formation.

Based on this analysis, an improved gating system, Scheme M, was designed. The key modifications included increasing the number and size of the ingates and adding strategically placed feeder heads (risers) to provide the necessary molten metal for feeding the upper section during solidification.

To systematically optimize the process, three key factors were identified: Pouring Temperature (A), Mold Shell Preheat Temperature (B), and Pouring Time (C). The target for optimization was the volume of shrinkage porosity within the casting body. The melting range for K4169 is 1243–1359°C, and pouring temperature is typically 100–150°C above the liquidus. Therefore, a range of 1400–1500°C was selected for Factor A. To avoid excessive thermal shock, the shell preheat temperature (Factor B) was varied between 950°C and 1050°C. Given the component’s relatively small volume, the pouring time (Factor C) was tested between 4 and 8 seconds. Each factor was studied at three levels, as detailed in the orthogonal array.

Factor Level 1 Level 2 Level 3
A: Pouring Temperature (°C) 1400 1450 1500
B: Shell Temperature (°C) 950 1000 1050
C: Pouring Time (s) 4 6 8

A standard L9 (3^4) Taguchi orthogonal array was used to design the experiments, with the three factors assigned to three columns. The response variable was the shrinkage porosity volume (Y) obtained from ProCAST simulations for each run using the improved Scheme M gating design.

Exp. No. A: Pouring Temp. (°C) B: Shell Temp. (°C) C: Pouring Time (s) Shrinkage Volume, Y (mm³) Signal-to-Noise Ratio, S/N (dB)
1 1400 950 4 29.733 -29.465
2 1400 1000 6 113.509 -41.101
3 1400 1050 8 130.530 -42.314
4 1450 950 6 34.418 -30.736
5 1450 1000 8 48.107 -33.644
6 1450 1050 4 8.347 -18.431
7 1500 950 8 13.133 -22.367
8 1500 1000 4 3.496 -10.871
9 1500 1050 6 10.661 -20.556

Table 1: Taguchi L9 Orthogonal Array Design and Simulation Results for Shrinkage Porosity.

The Signal-to-Noise (S/N) ratio was calculated for each experiment using the “smaller-the-better” characteristic, as the goal is to minimize shrinkage. The formula is:

$$ S/N = \eta = -10 \log_{10}\left(\frac{1}{n}\sum_{i=1}^{n} y_i^2\right) $$

where \( n \) is the number of repetitions (here, n=1 per simulation run) and \( y_i \) is the measured shrinkage volume. A higher (less negative) S/N ratio indicates better performance and greater robustness against noise factors.

Analysis of Results and Determination of Optimal Parameters

Direct observation of Table 1 shows that Experiment 8 (A3B2C1: 1500°C, 1000°C, 4s) yielded the minimum shrinkage volume (3.496 mm³) and the highest S/N ratio (-10.871 dB), indicating it as the best parameter combination within the designed experimental space.

To determine the statistical significance and influence of each factor, Analysis of Variance (ANOVA) and mean response analysis were performed on the S/N ratio data.

Factor Degrees of Freedom (f) Sum of Squares (S) Variance (V) F-Ratio Contribution (%)
A: Pouring Temperature 2 581.94 290.97 9.65 64.5%
B: Shell Temperature 2 46.57 23.28 0.77 5.2%
C: Pouring Time 2 303.42 151.71 5.03 33.6%
Error 2 30.16 15.08 -3.3%*
Total 8 901.77 100%

Table 2: Analysis of Variance (ANOVA) for S/N Ratios. (*Pooled error contribution is often recalculated; the negative value here indicates the factors explain more than 100% of the variation due to the small error term).

The ANOVA table (Table 2) clearly shows that Pouring Temperature (Factor A) has the highest F-ratio, confirming it as the most statistically significant factor affecting shrinkage porosity. Pouring Time (Factor C) is also significant, while Shell Preheat Temperature (Factor B) has a relatively minor influence within the tested range for this specific gating design.

The mean response for each factor level was calculated by averaging the S/N ratios for all experiments containing that level.

$$ M_{A1} = \frac{\eta_1 + \eta_2 + \eta_3}{3} = \frac{-29.465 -41.101 -42.314}{3} = -37.627 $$

$$ M_{A2} = \frac{\eta_4 + \eta_5 + \eta_6}{3} = \frac{-30.736 -33.644 -18.431}{3} = -27.604 $$

$$ M_{A3} = \frac{\eta_7 + \eta_8 + \eta_9}{3} = \frac{-22.367 -10.871 -20.556}{3} = -17.931 $$

Similarly, the means for Factors B and C (\(M_{Bj}\), \(M_{Ck}\)) are calculated. The range \(R\) for each factor, which is the difference between the maximum and minimum mean S/N, quantifies the factor’s effect magnitude.

$$ R_A = \max(M_{Ai}) – \min(M_{Ai}) = (-17.931) – (-37.627) = 19.696 $$
$$ R_B = 1.439 $$
$$ R_C = 13.186 $$

The order of influence based on the range values is: Pouring Temperature (A) > Pouring Time (C) > Shell Temperature (B). This aligns with the ANOVA conclusions. The optimal level for each factor is the one corresponding to the highest mean S/N ratio: A3 (1500°C), B2 (1000°C), C1 (4s).

Therefore, the optimized precision investment casting parameters for the K4169 component with the Scheme M gating system are definitively established as: Pouring Temperature = 1500°C, Shell Preheat Temperature = 1000°C, Pouring Time = 4 seconds.

Validation and Discussion of Single-Factor Effects

A final simulation was run using the optimized parameters (A3B2C1) to validate the results. The filling process was smooth and complete. The solidification sequence showed a clear directional progression from the casting extremities towards the feeder heads. The feeder heads themselves solidified last, successfully acting as reservoirs to feed the shrinking casting. The shrinkage porosity defect was entirely isolated within the feeder heads, resulting in a sound casting body free from internal shrinkage, as predicted.

To further understand the underlying physical mechanisms, single-factor effect studies were conducted while keeping the other two factors at their optimal levels.

Effect of Pouring Temperature: With Shell Temperature at 1000°C and Pouring Time at 4s, the pouring temperature was varied from 1400°C to 1525°C. Shrinkage volume followed a distinct “V” trend, reaching a minimum at 1500°C. This can be modeled by a quadratic relationship:
$$ V_{sh} \approx \alpha (T – T_{opt})^2 + V_{min} $$
where \(V_{sh}\) is shrinkage volume, \(T\) is pouring temperature, \(T_{opt}\) is the optimal temperature (1500°C), \(V_{min}\) is the minimum shrinkage, and \(\alpha\) is a positive coefficient. Temperatures too low impair fluidity, leading to mistrun and poor feeding. Temperatures too high increase total heat content and may promote excessive grain growth and localized shrinkage due to prolonged liquid phase existence and greater contraction upon solidification. The optimal temperature balances fluidity with controlled solidification characteristics.

Effect of Shell Preheat Temperature: With Pouring Temperature at 1500°C and Pouring Time at 4s, the shell temperature was varied from 800°C to 1050°C. The effect was less pronounced, with shrinkage volume showing a shallow minimum around 1000°C. The relationship is milder and can be approximated linearly near the optimum:
$$ V_{sh} \approx \beta |B – B_{opt}| + V_{min} $$
where \(B\) is shell temperature and \(B_{opt}\) is 1000°C. A shell temperature too low creates a steep thermal gradient, potentially causing surface chilling and rapid solidification that hinders interdendritic feeding. A shell temperature too high reduces the gradient excessively, potentially leading to a wider mushy zone and equiaxed grain formation, which can also impede feeding. The optimal value promotes a favorable thermal gradient for directional solidification towards the feeders.

Effect of Pouring Time: With Pouring Temperature at 1500°C and Shell Temperature at 1000°C, the pouring time was varied from 3s to 8s. Shrinkage volume increased markedly with longer pouring times, showing an approximately exponential rise:
$$ V_{sh} \approx \gamma e^{\lambda t} + C $$
where \(t\) is pouring time and \(\gamma, \lambda, C\) are constants. A very short pour time (e.g., 3-4s) provides sufficient momentum for complete filling and establishes a beneficial thermal field. Excessively long pour times lead to significant heat loss in the early-poured metal, creating temperature stratification within the mold. This disrupts the desired progressive solidification front, causing isolated liquid pools and severe shrinkage. Therefore, a rapid but controlled fill is essential in precision investment casting of thin-walled sections.

The confluence of these single-factor studies confirms the optimal point identified by the Taguchi method. It also provides a scientific rationale: the optimized set (1500°C, 1000°C, 4s) achieves an ideal balance. The high pouring temperature ensures excellent fluidity. The matched, high shell temperature prevents premature chilling and maintains a controllable solidification rate. The short pouring time prevents thermal degradation and establishes a coherent thermal field for directional solidification, which is then effectively managed by the properly designed Scheme M gating and feeding system.

Gating System Comparison and Concluding Framework

The critical role of gating design cannot be overstated. The initial Scheme L, despite potential parameter adjustments, was fundamentally flawed due to its inability to establish effective feeding. Its solidification simulation revealed isolated liquid pools in the upper casting region, which inevitably transformed into macro-porosity. In contrast, Scheme M, by providing adequate feed metal reservoirs (risers) and improved filling characteristics, enabled a sound solidification pattern. This highlights a fundamental principle in precision investment casting: a robust gating and feeding design is a prerequisite for successful process parameter optimization. Parameters can only be fine-tuned within the window of feasibility created by a sound geometrical design.

In conclusion, this work demonstrates a comprehensive, simulation-driven framework for optimizing precision investment casting processes for complex high-temperature alloy components. The methodology integrates:

  1. Defect Diagnosis via Simulation: Using numerical tools like ProCAST to identify defect locations and root causes in initial designs.
  2. Design of Experiments for Efficiency: Applying the Taguchi method to systematically and efficiently explore the multi-parameter space (Pouring Temperature, Shell Temperature, Pouring Time) and identify their relative influence and optimal combination.
  3. Physical Validation and Insight: Using single-factor simulation studies to validate the optimum and understand the underlying thermal and fluid flow principles governing defect formation.

The key findings for this specific K4169 annular component are:

  1. The order of influence on shrinkage porosity is: Pouring Temperature > Pouring Time > Shell Preheat Temperature.
  2. The optimal precision investment casting parameters are: 1500°C Pouring Temperature, 1000°C Shell Temperature, and 4s Pouring Time.
  3. The gating system design is the most critical factor; an effective feeding system (like Scheme M) must be in place for parameter optimization to yield a defect-free casting.

This framework, combining numerical simulation, statistical experimental design, and fundamental process analysis, provides a powerful and generalizable approach for enhancing quality, yield, and performance in the demanding field of aerospace precision investment casting.

Scroll to Top