Third medium contact method

Last updated June 17, 2024

The third medium contact(TMC) is an implicit formulation for contact mechanics. Contacting bodies are embedded in a highly compliant medium (the third medium), which becomes increasingly stiff under compression. The stiffening of the third medium allows tractions to be transferred between the contacting bodies when the third medium between the bodies is compressed. In itself, the method is inexact; however, in contrast to most other contact methods, the third medium approach is continuous and differentiable, which makes it applicable to applications such as topology optimization.^[1]^[2]

The method was first proposed by Peter Wriggers et al. where a St. Venant-Kirchhoff material was used to model the third medium.^[3] This approach requires explicit treatment of surface normals. A simplification to the method was offered by Bog et al. by applying a Hencky material with the inherent property of becoming rigid under ultimate compression.^[4] This property has made the explicit treatment of surface normals redundant, thereby transforming the third medium contact method into a fully implicit method. The addition of a new regularization by Bluhm et al. to stabilize the third medium further extended the method to applications involving moderate sliding, rendering it practically applicable.^[1]

Methodology

A material with the property that it becomes increasingly stiff under compression is augmented by a regularization term. In terms of strain energy density, this may be expressed as

$\Psi (u)=W(u)+\mathbb {H} (u)\,{\boldsymbol {\scriptstyle {\vdots }}}\,\mathbb {H} (u)$ ,

where ${\textstyle \Psi (u)}$ represents the augmented strain energy density in the third medium, ${\textstyle \mathbb {H} (u)\,{\boldsymbol {\scriptstyle {\vdots }}}\,\mathbb {H} (u)}$ is the regularization term representing the inner product of the spatial Hessian by itself, and ${\textstyle W(u)}$ is the underlying strain energy density of the third medium, e.g. a Neo-Hookean solid or another hyperelastic material. The term ${\textstyle \mathbb {H} (u)\,{\boldsymbol {\scriptstyle {\vdots }}}\,\mathbb {H} (u)}$ is commonly referred to as HuHu-regularization.

Related Research Articles

Bra–ket notation, also called Dirac notation, is a notation for linear algebra and linear operators on complex vector spaces together with their dual space both in the finite-dimensional and infinite-dimensional case. It is specifically designed to ease the types of calculations that frequently come up in quantum mechanics. Its use in quantum mechanics is quite widespread.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

A wavelet is a wave-like oscillation with an amplitude that begins at zero, increases or decreases, and then returns to zero one or more times. Wavelets are termed a "brief oscillation". A taxonomy of wavelets has been established, based on the number and direction of its pulses. Wavelets are imbued with specific properties that make them useful for signal processing.

In continuum mechanics, stress is a physical quantity that describes forces present during deformation. For example, an object being pulled apart, such as a stretched elastic band, is subject to tensile stress and may undergo elongation. An object being pushed together, such as a crumpled sponge, is subject to compressive stress and may undergo shortening. The greater the force and the smaller the cross-sectional area of the body on which it acts, the greater the stress. Stress has dimension of force per area, with SI units of newtons per square meter (N/m²) or pascal (Pa).

A Newtonian fluid is a fluid in which the viscous stresses arising from its flow are at every point linearly correlated to the local strain rate — the rate of change of its deformation over time. Stresses are proportional to the rate of change of the fluid's velocity vector.

Quantum statistical mechanics is statistical mechanics applied to quantum mechanical systems. In quantum mechanics a statistical ensemble is described by a density operator S, which is a non-negative, self-adjoint, trace-class operator of trace 1 on the Hilbert space H describing the quantum system. This can be shown under various mathematical formalisms for quantum mechanics.

In mathematics, a norm is a function from a real or complex vector space to the non-negative real numbers that behaves in certain ways like the distance from the origin: it commutes with scaling, obeys a form of the triangle inequality, and is zero only at the origin. In particular, the Euclidean distance in an Euclidean space is defined by a norm on the associated Euclidean vector space, called the Euclidean norm, the 2-norm, or, sometimes, the magnitude of the vector. This norm can be defined as the square root of the inner product of a vector with itself.

In quantum information theory, a quantum channel is a communication channel which can transmit quantum information, as well as classical information. An example of quantum information is the state of a qubit. An example of classical information is a text document transmitted over the Internet.

Topology optimization is a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions and constraints with the goal of maximizing the performance of the system. Topology optimization is different from shape optimization and sizing optimization in the sense that the design can attain any shape within the design space, instead of dealing with predefined configurations.

<span class="mw-page-title-main">Level-set method</span> Conceptual framework used in numerical analysis of surfaces and shapes

The Level-set method (LSM) is a conceptual framework for using level sets as a tool for numerical analysis of surfaces and shapes. LSM can perform numerical computations involving curves and surfaces on a fixed Cartesian grid without having to parameterize these objects. LSM makes it easier to perform computations on shapes with sharp corners and shapes that change topology. These characteristics make LSM effective for modeling objects that vary in time, such as an airbag inflating or a drop of oil floating in water.

<span class="mw-page-title-main">S wave</span> Type of elastic body wave

In seismology and other areas involving elastic waves, S waves, secondary waves, or shear waves are a type of elastic wave and are one of the two main types of elastic body waves, so named because they move through the body of an object, unlike surface waves.

A strain energy density function or stored energy density function is a scalar-valued function that relates the strain energy density of a material to the deformation gradient.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

In applied mathematics, topological data analysis (TDA) is an approach to the analysis of datasets using techniques from topology. Extraction of information from datasets that are high-dimensional, incomplete and noisy is generally challenging. TDA provides a general framework to analyze such data in a manner that is insensitive to the particular metric chosen and provides dimensionality reduction and robustness to noise. Beyond this, it inherits functoriality, a fundamental concept of modern mathematics, from its topological nature, which allows it to adapt to new mathematical tools.

The structured support-vector machine is a machine learning algorithm that generalizes the Support-Vector Machine (SVM) classifier. Whereas the SVM classifier supports binary classification, multiclass classification and regression, the structured SVM allows training of a classifier for general structured output labels.

In continuum mechanics, a hypoelastic material is an elastic material that has a constitutive model independent of finite strain measures except in the linearized case. Hypoelastic material models are distinct from hyperelastic material models in that, except under special circumstances, they cannot be derived from a strain energy density function.

<span class="mw-page-title-main">Gradient discretisation method</span>

In numerical mathematics, the gradient discretisation method (GDM) is a framework which contains classical and recent numerical schemes for diffusion problems of various kinds: linear or non-linear, steady-state or time-dependent. The schemes may be conforming or non-conforming, and may rely on very general polygonal or polyhedral meshes.

Quantum optimization algorithms are quantum algorithms that are used to solve optimization problems. Mathematical optimization deals with finding the best solution to a problem from a set of possible solutions. Mostly, the optimization problem is formulated as a minimization problem, where one tries to minimize an error which depends on the solution: the optimal solution has the minimal error. Different optimization techniques are applied in various fields such as mechanics, economics and engineering, and as the complexity and amount of data involved rise, more efficient ways of solving optimization problems are needed. Quantum computing may allow problems which are not practically feasible on classical computers to be solved, or suggest a considerable speed up with respect to the best known classical algorithm.

The variational multiscale method (VMS) is a technique used for deriving models and numerical methods for multiscale phenomena. The VMS framework has been mainly applied to design stabilized finite element methods in which stability of the standard Galerkin method is not ensured both in terms of singular perturbation and of compatibility conditions with the finite element spaces.

The streamline upwind Petrov–Galerkin pressure-stabilizing Petrov–Galerkin formulation for incompressible Navier–Stokes equations can be used for finite element computations of high Reynolds number incompressible flow using equal order of finite element space by introducing additional stabilization terms in the Navier–Stokes Galerkin formulation.

References

1 2 Bluhm, Gore Lukas; Sigmund, Ole; Poulios, Konstantinos (2021-03-04). "Internal contact modeling for finite strain topology optimization". Computational Mechanics. 67 (4): 1099–1114. arXiv: 2010.14277 . Bibcode:2021CompM..67.1099B. doi:10.1007/s00466-021-01974-x. ISSN 0178-7675. S2CID 225076340.
↑ Frederiksen, Andreas Henrik; Sigmund, Ole; Poulios, Konstantinos (2023-10-07). "Topology optimization of self-contacting structures". Computational Mechanics. arXiv: 2305.06750 . doi: 10.1007/s00466-023-02396-7 . ISSN 1432-0924.
↑ Wriggers, P.; Schröder, J.; Schwarz, A. (2013-03-30). "A finite element method for contact using a third medium". Computational Mechanics. 52 (4): 837–847. Bibcode:2013CompM..52..837W. doi:10.1007/s00466-013-0848-5. ISSN 0178-7675. S2CID 254032357.
↑ Bog, Tino; Zander, Nils; Kollmannsberger, Stefan; Rank, Ernst (October 2015). "Normal contact with high order finite elements and a fictitious contact material". Computers & Mathematics with Applications. 70 (7): 1370–1390. doi:10.1016/j.camwa.2015.04.020. ISSN 0898-1221.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[:0-1] 1 2 Bluhm, Gore Lukas; Sigmund, Ole; Poulios, Konstantinos (2021-03-04). "Internal contact modeling for finite strain topology optimization". Computational Mechanics. 67 (4): 1099–1114. arXiv: 2010.14277 . Bibcode:2021CompM..67.1099B. doi:10.1007/s00466-021-01974-x. ISSN 0178-7675. S2CID 225076340.

[2] Frederiksen, Andreas Henrik; Sigmund, Ole; Poulios, Konstantinos (2023-10-07). "Topology optimization of self-contacting structures". Computational Mechanics. arXiv: 2305.06750 . doi: 10.1007/s00466-023-02396-7 . ISSN 1432-0924.

[3] Wriggers, P.; Schröder, J.; Schwarz, A. (2013-03-30). "A finite element method for contact using a third medium". Computational Mechanics. 52 (4): 837–847. Bibcode:2013CompM..52..837W. doi:10.1007/s00466-013-0848-5. ISSN 0178-7675. S2CID 254032357.

[4] Bog, Tino; Zander, Nils; Kollmannsberger, Stefan; Rank, Ernst (October 2015). "Normal contact with high order finite elements and a fictitious contact material". Computers & Mathematics with Applications. 70 (7): 1370–1390. doi:10.1016/j.camwa.2015.04.020. ISSN 0898-1221.

[1]

[2]

[3]

[4]