# RD-E: 5502 Rotating Fan Blade Ice Impact and Failure

Two hailstone ice balls are added to the model from RD-E: 5501 Fan Blade Rotation Initialization to study the deformation and failure of the blades when impacted.

The stress in the blades, due to rotation is accounted for by using Radioss state files, /STATE. A Johnson-Cook Failure Model, /FAIL/JOHNSON, is applied to the blade elements to take into account material failure.

## Options and Keywords Used

- Centrifugal force pre-load in rotating structures
- State files used to apply pre-load stress in a new simulation (/STATE)
- Implicit simulation options (Implicit Solution)
- Centrifugal force field (/LOAD/CENTRI)
- Rotational velocity about an axis (/INIV/AXIS/Z/1)
- Boundary Condition remove in Engine file (/BCSR)
- Johnson-Cook failure model (/FAIL/JOHNSON)

```
# To save the full stress use these 3 STATE SHELL options
/STATE/SHELL/STRESS/FULL
/STATE/SHELL/AUX/FULL
/STATE/SHELL/STRAIN/FULL
# time to write state file and which part
/STATE/DT
0.1 0.1
1
```

Now the hailstone ice ball model and the stress State file from the implicit simulation are included in the full simulation. The initial rotational velocity on the blades is included in the Engine file (/INIV/AXIS/Z/1) but could have also been be added Starter file using /INIVEL/AXIS. The ice impact simulation is ran for 0.08 seconds.

## Input Files

- <install_directory>/hwsolvers/demos/radioss/example/55_Fan_Blade_Out_FBO/2_Rotation_and_1_ice/*

## Model Description

Units: mm, s, Mg, N, and MPa

^{1}

**Blade Titanium Material Properties**- Density
- $4.43e-9\frac{Mg}{m{m}^{3}}$
- Young's modulus
- 113400 $\left[\mathrm{MPa}\right]$
- Poisson's ratio
- 0.342
- Yield stress
- 1098 $\left[\mathrm{MPa}\right]$
- Plastic hardening parameter
- 1092 $\left[\mathrm{MPa}\right]$
- Plastic hardening exponent
- 0.93

**Failure Material Properties**`D`_{1}- -0.09
`D`_{2}- 0.25
`D`_{3}- -0.05
`D`_{4}- 0.014
`I`_{fail_sh}- 1

**Case Steel Material Properties**- Density
- $7.9e-9\frac{Mg}{m{m}^{3}}$
- Young's modulus
- 210000 $\left[\mathrm{MPa}\right]$
- Poisson's ratio
- 0.3
- Yield stress
- 200 $\left[\mathrm{MPa}\right]$
- Plastic hardening parameter
- 450 $\left[\mathrm{MPa}\right]$
- Plastic hardening exponent
- 0.5
- Maximum stress
- 425 $\left[\mathrm{MPa}\right]$

**Ice Ball Material Properties**- Density
- $9.167e-10\frac{Mg}{m{m}^{3}}$
- Young's modulus
- 8304 $\left[\mathrm{MPa}\right]$
- Poisson's ratio
- 0.2
- Yield stress
- 10.3 $\left[\mathrm{MPa}\right]$

**Ice Ball Material Properties**- Plastic hardening parameter
- 2500 $\left[\mathrm{MPa}\right]$
- Plastic hardening exponent
- ${C}_{0}={C}_{2}={C}_{3}={C}_{4}={C}_{5}=0$

- Blade Center constrained all directions, except Rz
- Imposed Rotational Speed = 1000 = 104.72 $$\left[\frac{\text{rad}}{\text{s}}\right]$$
- Edges of case are fully constrained in X, Y, Z directions

### Model Method

This simulation builds on the results in RD-E: 5501 Fan Blade Rotation Initialization. The blades are impacted by two 3.07 kg hailstone ice balls with a radius of ~ 85 mm. A Johnson-Cook Failure model is added to the blade material to demonstrate how material failure can be modeled. The Johnson-Cook failure model relates plastic failure strain as a function of stress triaxiality (normalized mean stress). Additional details can be found in Materials of the Radioss User Guide, /FAIL/JOHNSON and /FAIL/TAB1.

## Results

### Conclusion

The stress in the blades from steady state rotation can be correctly accounted by including stress state files from a separate simulation. This saves time since it is not necessary to rerun the pre-load simulation for every design change of the case. The Johnson-Cook failure model accounts for the material failure of the blades due to their impact with hailstone ice balls.

^{1}Don Lesuer, Experimental Investigations of Material Models for Ti-6AL4V and 2024-T3, Lawrence Livermore National Laboratory, May 3, 1999