Abstract:This paper reports a numerical study based on a wake oscillator model, to determine the three-dimensional vortex-induced vibration (VIV) responses on a flexible cylinder with pinned-pinned boundary conditions subjected to a uniform flow. Four different aspect ratios have been selected for the study. The coupling equations of the structural oscillator models and wake oscillator models in both the cross-flow (CF) and in-line (IL) directions have been solved using a standard central difference method of the second order. The structural displacement, structural frequency, response wave pattern, response trajectory, and lift force coefficient, for four aspect ratios, have been compared. The numerical results establish that, for a small aspect ratio, the CF displacements have absolute standing wave behaviors without travelling wave behaviors, and the IL displacements have dominant standing wave behaviors with slight travelling wave behaviors. Further, the VIV trajectory is repeatable and displays figure-eight shapes. However, for a large aspect ratio, the CF displacements display identical characteristics with the IL displacements for a small aspect ratio, and the IL displacements for a large aspect ratio are simultaneously dominated by standing and travelling wave behaviors. Moreover, the VIV trajectory is apparently aperiodic and shows chaotic shapes.
What problem does this paper attempt to address?
This paper aims to solve the problem of response prediction of cylinders in three - dimensional vortex - induced vibration (VIV). Specifically, the research focuses on using the wake oscillator model to predict the response of a flexible cylinder with fixed - fixed boundary conditions in three - dimensional vortex - induced vibration in a uniform flow field. Four different aspect ratios (L/D = 100, 200, 500, 1000) were selected in the study. By numerically solving the coupled equations of the structural oscillator model and the wake oscillator model in the cross - flow (CF) and in - line (IL) directions, the changes in structural displacement, structural frequency, response waveform, response trajectory and lift coefficient under different aspect ratios were explored.
### Main research contents:
1. **Model description**:
- The research object is a flexible cylinder with length \( L \) and diameter \( D \), fixed at both ends and freely vibrating in the cross - flow (CF) and in - line (IL) directions.
- The vibrating of the flexible cylinder is described by using the tensioned beam model, taking into account the tension \( \Theta \) and bending stiffness \( EI \).
- The flow is assumed to be a uniform flow with a velocity of \( U \).
2. **Mathematical model**:
- The structural vibration equations in the cross - flow and in - line directions are respectively expressed as:
\[
m\frac{\partial^2 X(Z,T)}{\partial T^2} + R\frac{\partial X(Z,T)}{\partial T} - \Theta\frac{\partial^2 X(Z,T)}{\partial Z^2} + EI\frac{\partial^4 X(Z,T)}{\partial Z^4} = F_x(Z,T)
\]
\[
m\frac{\partial^2 Y(Z,T)}{\partial T^2} + R\frac{\partial Y(Z,T)}{\partial T} - \Theta\frac{\partial^2 Y(Z,T)}{\partial Z^2} + EI\frac{\partial^4 Y(Z,T)}{\partial Z^4} = F_y(Z,T)
\]
- The wake oscillator model equations are:
\[
\frac{\partial^2 p(Z,T)}{\partial T^2} + \varepsilon_x \Omega_f \left[p^2(Z,T) - 1\right] \frac{\partial p(Z,T)}{\partial T} + 4\Omega_f^2 p(Z,T) = A_x \frac{\partial^2 X(Z,T)}{\partial T^2}
\]
\[
\frac{\partial^2 q(Z,T)}{\partial T^2} + \varepsilon_y \Omega_f \left[q^2(Z,T) - 1\right] \frac{\partial q(Z,T)}{\partial T} + \Omega_f^2 q(Z,T) = A_y \frac{\partial^2 Y(Z,T)}{\partial T^2}
\]
3. **Numerical method**:
- The second - order central difference method is used to discretize the space and time domains to solve the above - mentioned coupled equations.
- The initial conditions are set to zero for displacement and velocity along the entire span.
- The boundary conditions are fixed at both ends, that is, both displacement and bending moment are zero.
4. **Results and discussion**:
- The changes in structural displacement, frequency, waveform, trajectory and lift coefficient under different aspect ratios were studied.
- It was found that for a small aspect ratio (such as L/D = 100), the displacement in the cross - flow direction mainly exhibits standing - wave behavior, while the displacement in the in - line direction exhibits both standing - wave and traveling - wave behaviors.