Reproducing kernel collocation method applied to the non-linear dynamics of pipe whip in a plane

Amit Shaw, D. Roy, Stephen R Reid, Majid Aleyaasin

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

A windowed collocation method, based on a moving least squares reproducing kernel particle approximation of functions, is explored for spatial discretization of the strongly non-linear system of partial differential equations governing large, planar whipping motion of a cantilever pipe Subjected to a follower force pulse (the blow-down force) normal to the deflected centreline at its tip. This problem was discussed by Reid et al. [An elastic plastic hardening-softening cantilever beam subjected to a force pulse at its tip: a model for pipe whip. Proc R Soc London A 1998;454:997-1029] where a space-time finite difference discretization was employed to solve the governing partial differential equation of motion. It was shown that, despite the deflected shape predictions being accurate, numerical solutions of these equations might exhibit problematic (possibly spurious) steep localized gradients. The resolution of this problem in the context of structural mechanics is novel and is the subject of this paper. In particular, it is demonstrated that it is possible to reduce significantly such spurious and localized numerical instabilities through a windowed collocation approach with a suitable choice of the window size. The collocation procedure presently adopted is based on the moving least squares reproducing kernel particle method. Material and structural non-linearity in the beam (pipe) model is incorporated via an elastic-plastic-hardening-softening moment-curvature relationship. The projected ordinary differential equations are then integrated in time through a fifth order, explicit Runge-Kutta method with adaptive step sizes.

The whipping motion of a pipe is often associated with the formation of one or more so-called moving plastic 'hinges' (regions). Depending on the magnitude and duration of the follower load pulse, plastic hinges may form at intermediate positions along the pipe. In such cases numerical solutions, especially those for curvatures and their derivatives, obtained through central difference or classical finite element methods may exhibit considerable fluctuations in space and time and may even 'blow up' following the formation of a kink at the location of the plastic hinge. However, it is shown that the presently adopted windowed collocation strategy has a far superior numerical performance, arresting these fluctuations to a great extent. Indeed, as demonstrated in this paper, there is a band of window sizes for which such fluctuations may nearly vanish.

Original languageEnglish
Pages (from-to)1637-1654
Number of pages18
JournalInternational Journal of Impact Engineering
Volume34
Issue number10
Early online date2 Jan 2007
DOIs
Publication statusPublished - Oct 2007

Keywords

  • pipe whip
  • collocation methods
  • travelling plastic hinges
  • large deflections
  • elastic-plastic-hardening-softening beams
  • large-deformation analysis
  • particle methods
  • element
  • formulation
  • cantilever
  • partition
  • force
  • space
  • model
  • tip

Cite this

Reproducing kernel collocation method applied to the non-linear dynamics of pipe whip in a plane. / Shaw, Amit; Roy, D.; Reid, Stephen R; Aleyaasin, Majid.

In: International Journal of Impact Engineering, Vol. 34, No. 10, 10.2007, p. 1637-1654.

Research output: Contribution to journalArticle

@article{29b6c1f734d041a3bdcf9cfa128f8290,
title = "Reproducing kernel collocation method applied to the non-linear dynamics of pipe whip in a plane",
abstract = "A windowed collocation method, based on a moving least squares reproducing kernel particle approximation of functions, is explored for spatial discretization of the strongly non-linear system of partial differential equations governing large, planar whipping motion of a cantilever pipe Subjected to a follower force pulse (the blow-down force) normal to the deflected centreline at its tip. This problem was discussed by Reid et al. [An elastic plastic hardening-softening cantilever beam subjected to a force pulse at its tip: a model for pipe whip. Proc R Soc London A 1998;454:997-1029] where a space-time finite difference discretization was employed to solve the governing partial differential equation of motion. It was shown that, despite the deflected shape predictions being accurate, numerical solutions of these equations might exhibit problematic (possibly spurious) steep localized gradients. The resolution of this problem in the context of structural mechanics is novel and is the subject of this paper. In particular, it is demonstrated that it is possible to reduce significantly such spurious and localized numerical instabilities through a windowed collocation approach with a suitable choice of the window size. The collocation procedure presently adopted is based on the moving least squares reproducing kernel particle method. Material and structural non-linearity in the beam (pipe) model is incorporated via an elastic-plastic-hardening-softening moment-curvature relationship. The projected ordinary differential equations are then integrated in time through a fifth order, explicit Runge-Kutta method with adaptive step sizes. The whipping motion of a pipe is often associated with the formation of one or more so-called moving plastic 'hinges' (regions). Depending on the magnitude and duration of the follower load pulse, plastic hinges may form at intermediate positions along the pipe. In such cases numerical solutions, especially those for curvatures and their derivatives, obtained through central difference or classical finite element methods may exhibit considerable fluctuations in space and time and may even 'blow up' following the formation of a kink at the location of the plastic hinge. However, it is shown that the presently adopted windowed collocation strategy has a far superior numerical performance, arresting these fluctuations to a great extent. Indeed, as demonstrated in this paper, there is a band of window sizes for which such fluctuations may nearly vanish.",
keywords = "pipe whip, collocation methods, travelling plastic hinges, large deflections, elastic-plastic-hardening-softening beams, large-deformation analysis, particle methods, element, formulation, cantilever, partition, force, space, model, tip",
author = "Amit Shaw and D. Roy and Reid, {Stephen R} and Majid Aleyaasin",
year = "2007",
month = "10",
doi = "10.1016/J.IJIMPENG.2006.09.004",
language = "English",
volume = "34",
pages = "1637--1654",
journal = "International Journal of Impact Engineering",
issn = "0734-743X",
publisher = "Elsevier Limited",
number = "10",

}

TY - JOUR

T1 - Reproducing kernel collocation method applied to the non-linear dynamics of pipe whip in a plane

AU - Shaw, Amit

AU - Roy, D.

AU - Reid, Stephen R

AU - Aleyaasin, Majid

PY - 2007/10

Y1 - 2007/10

N2 - A windowed collocation method, based on a moving least squares reproducing kernel particle approximation of functions, is explored for spatial discretization of the strongly non-linear system of partial differential equations governing large, planar whipping motion of a cantilever pipe Subjected to a follower force pulse (the blow-down force) normal to the deflected centreline at its tip. This problem was discussed by Reid et al. [An elastic plastic hardening-softening cantilever beam subjected to a force pulse at its tip: a model for pipe whip. Proc R Soc London A 1998;454:997-1029] where a space-time finite difference discretization was employed to solve the governing partial differential equation of motion. It was shown that, despite the deflected shape predictions being accurate, numerical solutions of these equations might exhibit problematic (possibly spurious) steep localized gradients. The resolution of this problem in the context of structural mechanics is novel and is the subject of this paper. In particular, it is demonstrated that it is possible to reduce significantly such spurious and localized numerical instabilities through a windowed collocation approach with a suitable choice of the window size. The collocation procedure presently adopted is based on the moving least squares reproducing kernel particle method. Material and structural non-linearity in the beam (pipe) model is incorporated via an elastic-plastic-hardening-softening moment-curvature relationship. The projected ordinary differential equations are then integrated in time through a fifth order, explicit Runge-Kutta method with adaptive step sizes. The whipping motion of a pipe is often associated with the formation of one or more so-called moving plastic 'hinges' (regions). Depending on the magnitude and duration of the follower load pulse, plastic hinges may form at intermediate positions along the pipe. In such cases numerical solutions, especially those for curvatures and their derivatives, obtained through central difference or classical finite element methods may exhibit considerable fluctuations in space and time and may even 'blow up' following the formation of a kink at the location of the plastic hinge. However, it is shown that the presently adopted windowed collocation strategy has a far superior numerical performance, arresting these fluctuations to a great extent. Indeed, as demonstrated in this paper, there is a band of window sizes for which such fluctuations may nearly vanish.

AB - A windowed collocation method, based on a moving least squares reproducing kernel particle approximation of functions, is explored for spatial discretization of the strongly non-linear system of partial differential equations governing large, planar whipping motion of a cantilever pipe Subjected to a follower force pulse (the blow-down force) normal to the deflected centreline at its tip. This problem was discussed by Reid et al. [An elastic plastic hardening-softening cantilever beam subjected to a force pulse at its tip: a model for pipe whip. Proc R Soc London A 1998;454:997-1029] where a space-time finite difference discretization was employed to solve the governing partial differential equation of motion. It was shown that, despite the deflected shape predictions being accurate, numerical solutions of these equations might exhibit problematic (possibly spurious) steep localized gradients. The resolution of this problem in the context of structural mechanics is novel and is the subject of this paper. In particular, it is demonstrated that it is possible to reduce significantly such spurious and localized numerical instabilities through a windowed collocation approach with a suitable choice of the window size. The collocation procedure presently adopted is based on the moving least squares reproducing kernel particle method. Material and structural non-linearity in the beam (pipe) model is incorporated via an elastic-plastic-hardening-softening moment-curvature relationship. The projected ordinary differential equations are then integrated in time through a fifth order, explicit Runge-Kutta method with adaptive step sizes. The whipping motion of a pipe is often associated with the formation of one or more so-called moving plastic 'hinges' (regions). Depending on the magnitude and duration of the follower load pulse, plastic hinges may form at intermediate positions along the pipe. In such cases numerical solutions, especially those for curvatures and their derivatives, obtained through central difference or classical finite element methods may exhibit considerable fluctuations in space and time and may even 'blow up' following the formation of a kink at the location of the plastic hinge. However, it is shown that the presently adopted windowed collocation strategy has a far superior numerical performance, arresting these fluctuations to a great extent. Indeed, as demonstrated in this paper, there is a band of window sizes for which such fluctuations may nearly vanish.

KW - pipe whip

KW - collocation methods

KW - travelling plastic hinges

KW - large deflections

KW - elastic-plastic-hardening-softening beams

KW - large-deformation analysis

KW - particle methods

KW - element

KW - formulation

KW - cantilever

KW - partition

KW - force

KW - space

KW - model

KW - tip

U2 - 10.1016/J.IJIMPENG.2006.09.004

DO - 10.1016/J.IJIMPENG.2006.09.004

M3 - Article

VL - 34

SP - 1637

EP - 1654

JO - International Journal of Impact Engineering

JF - International Journal of Impact Engineering

SN - 0734-743X

IS - 10

ER -