45.6. Boundary Conditions and Loading

Fixed-support conditions are applied near the disk portion of the cyclic-sector fan blade model, as shown in the following figure:

Figure 45.9: Fixed Support Condition Near Disk

Fixed Support Condition Near Disk

The following loads are considered for the cyclic symmetry analyses:

  • Centrifugal loads caused by rotational velocity

  • Thermal loads due to the difference in reference and operating temperature

  • Unsteady flow pressure applied on the fan blade (imported from CFX)

The rotational velocity (CGOMGA,0,0,1680) is applied along the global Z-axis. The bladed disk is modeled with a coefficient of thermal expansion of 1.2 x 10-5 °C-1. The reference temperature is held at 22 °C and an operating temperature of 50 °C is applied on all nodes of the model to generate the thermal load vector.

The thermal load vector generated from the base static solve can be ignored in the subsequent analysis (THEXPAND).

The unsteady flow pressures, which are imported from Ansys CFX, are generated for EO = 2 engine order excitation at a rotational frequency of 534.76 Hz. These pressure data are then mapped to the structural rotor fan blade model in Ansys Mechanical APDL using its mapping capabilities via the /MAP processor.

The following example input shows the steps involved in this mapping procedure:

/map                                   ! Mapping processor
target, nbf1                           ! Specifies target nodes for mapping pressures onto 
                                       ! surface effect elements
FTYPE, cfxtbr,1                        ! Specifies file and pressure type for the 
                                       ! subsequent import of source points and pressures
READ,CFX_ExportResults_FT_10P_EO2.csv' ! Read coordinate and pressure data from a file

allsel,all,all
csys,0

/show,jpeg,rev
plgeom                                 ! Plots source and target geometries
MAP,,2,,1,                             ! Map pressures from source points to target surface elements
plmap,target                           ! Plot target and source pressure		
plmap,target,,,1
/show,close
WRITEMAP,'mapped_on_cyclic_model.dat'  ! Write interpolated pressure 			
                                       ! data to a file
finish

To illustrate the step-by-step mapping procedure, the complete workflow in Mechanical APDL is shown below:

Using the PLGEOM command, the target and source geometries can be plotted as shown:

Figure 45.10: Target and Source Geometries of the Rotor Fan Blade

Target and Source Geometries of the Rotor Fan Blade

The real and imaginary components of mapped target pressure are plotted using the PLMAP command, as shown:

Figure 45.11: Mapped Target Pressure Plots on the Rotor Fan Blade Surface

Mapped Target Pressure Plots on the Rotor Fan Blade Surface

The procedure for computing aerodynamic coefficients (stiffness and damping terms) for each interblade phase angle is described in Aero Coupling. The AEROCOEFF command uses unsteady pressures obtained from a CFD flutter analysis to compute an aerodynamic coefficient array. This aerodynamic coefficient array is compatible with the CYCFREQ,AERO command. If all interblade phase angles are not available from the CFD computation, the coefficients can be interpolated between nodal diameters and inserted into the aerodynamic coefficient array.

In the case of a mistuned bladed disk, blades are subjected to engine order (EO) excitation. Engine order excitation is the effective traveling wave excitation that a bladed disk experiences as it passes through disturbances in the flow-field for each revolution. For example, one disturbance in the flow field results in an EO = 1 excitation; two disturbances in the flow field result in an EO = 2 excitation.

In Mechanical APDL, the engine order excitation is applied using the CYCFREQ command with Option = EO. Internally, the program computes the aliased engine order, including its sign. The aliased engine order is determined from the input engine order as follows:

Table 45.1: Aliased Engine Order (Excited Harmonic Index)

Engine Order CAliased Engine Order
N EvenN Odd
C ≤ N / 2C ≤ (N - 1) / 2C
N / 2 < C ≤ N - 12(N + 1) ≤ C ≤ NN - C
N ≤ C ≤ 3N / 2N ≤ C ≤ (3N - 1) / 2C - N
3N / 2 ≤ C ≤ 2N(3N +1) / 2 ≤ C ≤ 2N2N - C

where N = number of sectors.

For the case of mistuned forced-response analysis, the forcing/excitation frequency expressed in terms of EO is as follows:

where is in Hz and is in RPM.

The excitation frequency is specified with the HARFRQ command.

In general, mistuning in bladed disks is typically modeled as small, random perturbations to the stiffness of the blade DOF. The deviation in stiffness for blade n is expressed as follows:

where is the mistuning parameter for blade n.

The mistuning parameters are provided in an array parameter of size N × 1, where N is the number of blades. These mistuning parameters are defined using the CYCFREQ command with Option = MIST and Value1 = K.

For the perturbed mode-superposition harmonic analysis, the unsteady flow pressures applied as per the engine order excitation on the fan blade are treated as a harmonically varying load.