Computation of Frequencies and Nodal Diameters

Blade Row Frequencies Computation

For a case with N blade rows, the complete set of blade passing frequencies for blade row j is:

Figure 1. EQUATION_DISPLAY
ω j = i = 1 i j N [ m i B i ( Ω j - Ω i ) + f i ]
(5050)

where, for blade row or wake i , m i is the blade row index and can take all values between - M i and M i . M i is the specified number of modes.

When using the Consider All Blade Rows option, i represents any blade row, wake, time varying boundary, or flutter condition, except blade row j . The relative rotation rate between any two blade rows is kept non-zero by adding a small fictitious rotation rate to the blade row if necessary; see Fictitious Rotation Rate Computation. When using the Consider Only Neighbors option, the only blade rows that are considered are the rows next to j that share a common indirect interface with j , along with any wake, time varying boundary, or flutter conditions that are directly attached to blade row j .

The variables B i , Ω i , and f i have different meanings, depending on what kind of blade row they refer to. These meanings are described in the following table.

Type of Blade Row i Meaning of B i Meaning of Ω i Meaning of f i
Modeled blade row The number of blades within the given pitch as specified in the Number of Blades property of the Blade Row n > Blades Per Pitch node. See Blades Per Pitch Properties. Rotation rate of the modeled blade row, which is inferred from the motion specification of the blade row region. See Defining Region Motion. Always zero. f i is not used for a modeled blade row.
Wake The number of blades within the given pitch as specified in the Number of Blades property of the Wake > Blades Per Pitch node. See Blades Per Pitch Properties. Rotation rate that is specified in the Rotation Rate property of the Wake node under Regions > Boundary > Physics Values. See Understanding Wake Options. Always zero. f i is not used for a wake blade row.
Time-varying option Set equal to ( B r σ ) / ( 2 π ) , where B r is the number of blades in the modeled blade row to which this boundary condition applies, and σ is specified in the Phase Angle property of the Time Varying Option Parameters node. See Time Varying Option. Set equal to the rotation rate of the modeled blade row to which this boundary condition applies. The fundamental frequency of variation in the boundary conditions, in radians/s. Specified in the Frequency property of the Time Varying Option Parameters node. See Time Varying Option.
Blade flutter Set equal to ( B r σ i b ) / ( 2 π ) , where B r is the number of blades in the modeled blade row to which this flutter condition applies, and σ i b is the inter-blade phase angle that is specified in the Flutter Phase Angle property of the Flutter Parameters node. See Flutter Parameters Properties. Set equal to the rotation rate of the modeled blade row to which this flutter condition applies. The frequency of unsteadiness ω used in Eqn. (5059). Specified in the Flutter Frequency property of the Flutter Parameters node. See Flutter Parameters Properties.

The maximum value for ω j is set with the Modes Table Filter Option. When this option is set to Cut-Off, the model computes the maximum value ω max using the cut-off factor γ :

Figure 2. EQUATION_DISPLAY
ωmax=MAXi{γMAXj[Mj[Bj(ΩiΩj)+fj]]}
(5051)

When the option is set to Maximum Limits, the model sets ω max to the value of the Maximum Frequency property.

Blade Row Nodal Diameter Computation

For a case with N blade rows, the complete set of nodal diameters that correspond to each frequency computed above (in Eqn. (5050)) for blade row j is:

Figure 3. EQUATION_DISPLAY
N j = i = 1 i j N m i B i
(5052)

where, for blade row or wake i :

  • m i is the modal index and can take all values between - M i and M i , where M i is the number of modes.
  • B i is the number of blades.

When using the Consider All Blade Rows option, i represents all blade rows except blade row j . When using the Consider Only Neighbors option, only the adjacent blade rows that share a common indirect interface with blade row j are considered.

The maximum value for N j is set with the Modes Table Filter Option. When this option is set to Cut-Off, the model computes the maximum value N max using the cut-off factor γ :

Figure 4. EQUATION_DISPLAY
N max = M A X i { γ M A X j [ M j B j ] }
(5053)

When the option is set to Maximum Limits, the model sets N max to the value of the Maximum Nodal Diameter property.

Fictitious Rotation Rate Computation

When you use Consider All Blade Rows as your Frequency Compute Option, every blade must be assigned a different frequency for the purposes of computation, even if some have identical frequencies in the User Specified Mode Table. The Fictitious Rotation Rate Factor ε determines the offset from the originally specified rotation rate.

The offset is an integer multiple of f where:

Figure 5. EQUATION_DISPLAY
f = ε Ω max
(5054)

Ω max is the maximum user-specified rotation rate in the blade row. The default value of ε is 0.001. You can set the value of ε in the Fictitious Rotation Rate Factor property of the Blade Rows manager node, but it is best to make the value small, so that the departure from the original value stays small.