A linear complementarity problem formulation for periodic solutions to unilateral contact problems

Markus B. Meingast, Mathias Legrand, Christophe Pierre

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Presented is an approach for finding periodic responses of structural systems subject to unilateral contact conditions. No other non-linear terms, e.g. large displacements or strains, hyper-elasticity, plasticity, etc. are considered. The excitation period due to various forcing conditions - from harmonic external or contact forcing due to a moving contact interface - is discretized in time, such that the quantities of interest - displacement, velocity, acceleration as well as contact force - can be approximated through time-domain schemes such as backward difference, Galerkin, and Fourier. The solution is assumed to exist and is defined on a circle with circumference T to directly enforce its periodicity. The strategy for approximating time derivative terms within the discretized period, i.e. velocity and acceleration, is hence circulant in nature. This results in a global circulant algebraic system of equations with inequalities that can be translated into a unique Linear Complementarity Problem (LCP). The LCP is then solved by dedicated and established methods such as Lemke's Algorithm. This allows for the computation of approximate periodic solutions exactly satisfying unilateral contact constraints on a discrete time set. The implementation efficiency and accuracy are discussed in comparison to classical time marching techniques for initial value problems combined with a Lagrange Multiplier contact treatment. The LCP algorithm is validated using a simple academic problem. The extension to large-scale systems is made possible through the implementation of a Craig-Bampton type Modal Component Synthesis. The method shows applicability to industrial rotor/casing contact set-ups as shown by studying a compressor blade. A good agreement to time marching simulations is found, suggesting a viable alternative to time marching or Fourier-based harmonic balance simulations.

Original languageEnglish (US)
Pages (from-to)18-27
Number of pages10
JournalInternational Journal of Non-Linear Mechanics
Volume66
DOIs
StatePublished - Jan 1 2014

Fingerprint

Unilateral Problems
Unilateral Contact
Linear Complementarity Problem
Contact Problem
Periodic Solution
Formulation
Lagrange multipliers
Initial value problems
Contact
Turbomachine blades
Plasticity
Compressors
Large scale systems
Elasticity
Rotors
Forcing
Derivatives
Hyperelasticity
Harmonic Balance
Large Displacements

Cite this

@article{1fc648f2fedd4e8883c3051aa8a0378e,
title = "A linear complementarity problem formulation for periodic solutions to unilateral contact problems",
abstract = "Presented is an approach for finding periodic responses of structural systems subject to unilateral contact conditions. No other non-linear terms, e.g. large displacements or strains, hyper-elasticity, plasticity, etc. are considered. The excitation period due to various forcing conditions - from harmonic external or contact forcing due to a moving contact interface - is discretized in time, such that the quantities of interest - displacement, velocity, acceleration as well as contact force - can be approximated through time-domain schemes such as backward difference, Galerkin, and Fourier. The solution is assumed to exist and is defined on a circle with circumference T to directly enforce its periodicity. The strategy for approximating time derivative terms within the discretized period, i.e. velocity and acceleration, is hence circulant in nature. This results in a global circulant algebraic system of equations with inequalities that can be translated into a unique Linear Complementarity Problem (LCP). The LCP is then solved by dedicated and established methods such as Lemke's Algorithm. This allows for the computation of approximate periodic solutions exactly satisfying unilateral contact constraints on a discrete time set. The implementation efficiency and accuracy are discussed in comparison to classical time marching techniques for initial value problems combined with a Lagrange Multiplier contact treatment. The LCP algorithm is validated using a simple academic problem. The extension to large-scale systems is made possible through the implementation of a Craig-Bampton type Modal Component Synthesis. The method shows applicability to industrial rotor/casing contact set-ups as shown by studying a compressor blade. A good agreement to time marching simulations is found, suggesting a viable alternative to time marching or Fourier-based harmonic balance simulations.",
author = "Meingast, {Markus B.} and Mathias Legrand and Christophe Pierre",
year = "2014",
month = "1",
day = "1",
doi = "https://doi.org/10.1016/j.ijnonlinmec.2014.01.007",
language = "English (US)",
volume = "66",
pages = "18--27",
journal = "International Journal of Non-Linear Mechanics",
issn = "0020-7462",
publisher = "Elsevier Ltd",

}

A linear complementarity problem formulation for periodic solutions to unilateral contact problems. / Meingast, Markus B.; Legrand, Mathias; Pierre, Christophe.

In: International Journal of Non-Linear Mechanics, Vol. 66, 01.01.2014, p. 18-27.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A linear complementarity problem formulation for periodic solutions to unilateral contact problems

AU - Meingast, Markus B.

AU - Legrand, Mathias

AU - Pierre, Christophe

PY - 2014/1/1

Y1 - 2014/1/1

N2 - Presented is an approach for finding periodic responses of structural systems subject to unilateral contact conditions. No other non-linear terms, e.g. large displacements or strains, hyper-elasticity, plasticity, etc. are considered. The excitation period due to various forcing conditions - from harmonic external or contact forcing due to a moving contact interface - is discretized in time, such that the quantities of interest - displacement, velocity, acceleration as well as contact force - can be approximated through time-domain schemes such as backward difference, Galerkin, and Fourier. The solution is assumed to exist and is defined on a circle with circumference T to directly enforce its periodicity. The strategy for approximating time derivative terms within the discretized period, i.e. velocity and acceleration, is hence circulant in nature. This results in a global circulant algebraic system of equations with inequalities that can be translated into a unique Linear Complementarity Problem (LCP). The LCP is then solved by dedicated and established methods such as Lemke's Algorithm. This allows for the computation of approximate periodic solutions exactly satisfying unilateral contact constraints on a discrete time set. The implementation efficiency and accuracy are discussed in comparison to classical time marching techniques for initial value problems combined with a Lagrange Multiplier contact treatment. The LCP algorithm is validated using a simple academic problem. The extension to large-scale systems is made possible through the implementation of a Craig-Bampton type Modal Component Synthesis. The method shows applicability to industrial rotor/casing contact set-ups as shown by studying a compressor blade. A good agreement to time marching simulations is found, suggesting a viable alternative to time marching or Fourier-based harmonic balance simulations.

AB - Presented is an approach for finding periodic responses of structural systems subject to unilateral contact conditions. No other non-linear terms, e.g. large displacements or strains, hyper-elasticity, plasticity, etc. are considered. The excitation period due to various forcing conditions - from harmonic external or contact forcing due to a moving contact interface - is discretized in time, such that the quantities of interest - displacement, velocity, acceleration as well as contact force - can be approximated through time-domain schemes such as backward difference, Galerkin, and Fourier. The solution is assumed to exist and is defined on a circle with circumference T to directly enforce its periodicity. The strategy for approximating time derivative terms within the discretized period, i.e. velocity and acceleration, is hence circulant in nature. This results in a global circulant algebraic system of equations with inequalities that can be translated into a unique Linear Complementarity Problem (LCP). The LCP is then solved by dedicated and established methods such as Lemke's Algorithm. This allows for the computation of approximate periodic solutions exactly satisfying unilateral contact constraints on a discrete time set. The implementation efficiency and accuracy are discussed in comparison to classical time marching techniques for initial value problems combined with a Lagrange Multiplier contact treatment. The LCP algorithm is validated using a simple academic problem. The extension to large-scale systems is made possible through the implementation of a Craig-Bampton type Modal Component Synthesis. The method shows applicability to industrial rotor/casing contact set-ups as shown by studying a compressor blade. A good agreement to time marching simulations is found, suggesting a viable alternative to time marching or Fourier-based harmonic balance simulations.

UR - http://www.scopus.com/inward/record.url?scp=84907056833&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84907056833&partnerID=8YFLogxK

U2 - https://doi.org/10.1016/j.ijnonlinmec.2014.01.007

DO - https://doi.org/10.1016/j.ijnonlinmec.2014.01.007

M3 - Article

VL - 66

SP - 18

EP - 27

JO - International Journal of Non-Linear Mechanics

JF - International Journal of Non-Linear Mechanics

SN - 0020-7462

ER -