Source author record

Paul Steinmann

Paul Steinmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2025arXiv

Cutting soft materials: how material differences shape the response

Cutting soft materials is a complex process governed by the interplay of bulk large deformation, interfacial soft fracture, and contact forces with the cutting tool. Existing experimental characterizations and numerical models often fail to capture the variety of observed cutting behaviors, especially the transition from indentation to cutting and the roles of dissipative mechanisms. Here, we combine novel experimental cutting tests on three representative materials-a soft hydrogel, elastomer, and food-based materials-with a coupled computational model that integrates soft fracture, adhesion, and frictional interactions. Our experiments reveal material-dependent cutting behaviors, with abrupt or smooth transitions from indentation to crack initiation, followed by distinct steady cutting regimes. The computational model captures these behaviors and shows that adhesion and viscous cohesive forces dominate tangential stresses, while Coulomb friction plays a negligible role due to low contact pressures. Together, these results provide new mechanistic insights into the physics of soft cutting and offer a unified framework to guide the design of soft materials, cutting tools, and cutting protocols, with direct relevance to surgical applications and the engineering of food textures optimized for mastication.

preprint2022arXiv

LEFM is agnostic to geometrical nonlinearities arising at atomistic crack tips

Various fields such as mechanical engineering, materials science, etc., have seen a widespread use of linear elastic fracture mechanics (LEFM) at the continuum scale. LEFM is also routinely applied to the atomic scale. However, its applicability at this scale remains less well studied, with most studies focusing on non-linear elastic effects. Using a harmonic (snapping spring) nearest-neighbor potential which provides the closest match to LEFM on a discrete lattice, we show that the discrete nature of an atomic lattice leads to deviations from the LEFM displacement field during energy minimization. We propose that these deviations can be ascribed to geometrical nonlinearities since the material does not have a nonlinear elastic response prior to bond breaking. We demonstrate that crack advance and the critical stress intensity factor in an incremental loading scenario is governed by the collectively loaded region, and can not be determined analytically from the properties (max. elongation, max. sustained force, etc.) of the stressed crack tip bond alone.

preprint2022arXiv

On continuum modeling of cell aggregation phenomena

Cellular aggregates play a significant role in the evolution of biological systems such as tumor growth, tissue spreading, wound healing, and biofilm formation. Analysis of such biological systems, in principle, includes examining the interplay of cell-cell interactions together with the cell-matrix interaction. These two interaction types mainly drive the dynamics of cellular aggregates which is intrinsically out of equilibrium. Here we propose a non-linear continuum mechanics formulation and the corresponding finite element simulation framework to model the physics of cellular aggregate formation. As an example, we focus in particular on the process of bacterial colony formation as recently studied by Kuan et al. Thereby we describe the aggregation process as an active phase separation phenomenon. We develop a Lagrangian continuum description of the problem which yields a substantial simplification to the formulations of the governing equations. Due to the presence of spatial Hessian and Laplacian operators, a gradient-enhanced approach is required to incorporate C1 continuity. In addition, a robust and efficient finite element formulation of the problem is provided. Taylor-Hood finite elements are utilized for the implementation to avoid instabilities related to the LBB condition. Finally, through a set of numerical examples, the influence of various parameters on the dynamics of the cellular aggregate formation is investigated. Our proposed methodology furnishes a general framework for the investigation of the rheology and non-equilibrium dynamics of cellular aggregates.

preprint2021arXiv

A p-adaptive, implicit-explicit mixed finite element method for reaction-diffusion problems

A new class of implicit-explicit (IMEX) methods combined with a p-adaptive mixed finite element formulation is proposed to simulate the diffusion of reacting species. Hierarchical polynomial functions are used to construct an $H(\mathrm{Div})$-conforming base for the flux vectors, and a non-conforming $L^2$ base for the mass concentration of the species. The mixed formulation captures the distinct nonlinearities associated with the constitutive flux equations and the reaction terms. The IMEX method conveniently treats these two sources of nonlinearity implicitly and explicitly, respectively, within a single time-stepping framework. The combination of the p-adaptive mixed formulation and the IMEX method delivers a robust and efficient algorithm. The proposed methods eliminate the coupled effect of mesh size and time step on the algorithmic stability. A residual based a posteriori error estimate that provides an upper bound of the natural error norm is derived. The availability of such estimate which can be obtained with minimal computational effort and the hierarchical construction of the finite element spaces allow for the formulation of an efficient p-adaptive algorithm. A series of numerical examples demonstrate the performance of the approach. It is shown that the method with the p-adaptive strategy accurately solves problems involving travelling waves, and those with discontinuities and singularities. The flexibility of the formulation is also illustrated via selected applications in pattern formation and electrophysiology.

preprint2021arXiv

A Variational Integrator for the Discrete Element Method

A novel implicit integration scheme for the Discrete Element Method (DEM) based on the variational integrator approach is presented. The numerical solver provides a fully dynamical description that, notably, reduces to an energy minimisation scheme in the quasi-static limit. A detailed derivation of the numerical method is presented for the Hookean contact model and tested against an established open source DEM package that uses the velocity-Verlet integration scheme. These tests compare results for a single collision, long-term stability and statistical quantities of ensembles of particles. Numerically, the proposed integration method demonstrates equivalent accuracy to the velocity-Verlet method.

preprint2020arXiv

Assessment of an Isogeometric Approach with Catmull-Clark Subdivision Surfaces using the Laplace-Beltrami Problems

An isogeometric approach for solving the Laplace-Beltrami equation on a two-dimensional manifold embedded in three-dimensional space using a Galerkin method based on Catmull-Clark subdivision surfaces is presented and assessed. The scalar-valued Laplace-Beltrami equation requires only C0 continuity and is adopted to elucidate key features and properties of the isogeometric method using Catmull-Clark subdivision surfaces. Catmull-Clark subdivision bases are used to discretise both the geometry and the physical field. A fitting method generates control meshes to approximate any given geometry with Catmull-Clark subdivision surfaces. The performance of the Catmull-Clark subdivision method is compared to the conventional finite element method. Subdivision surfaces without extraordinary vertices show the optimal convergence rate. However, extraordinary vertices introduce error, which decreases the convergence rate. A comparative study shows the effect of the number and valences of the extraordinary vertices on accuracy and convergence. An adaptive quadrature scheme is shown to reduce the error.

preprint2020arXiv

Concurrent consideration of cortical and cancellous bone within continuum bone remodelling

Continuum bone remodelling is an important tool for predicting the effects of mechanical stimuli on bone density evolution. While the modelling of only cancellous bone is considered in many studies based on continuum bone remodelling, this work presents an approach of modelling also cortical bone and the interaction of both bone types. The distinction between bone types is made by introducing an initial volume fraction. A simple point-wise example is used to study the behaviour of novel model options, as well as a proximal femur example, where the interaction of both bone types is demonstrated using initial density distributions. The results of the proposed model options indicate that the consideration of cortical bone remarkably changes the density evolution of cancellous bone, and should therefore not be neglected.

preprint2015arXiv

A finite element implementation of surface elasticity at finite strains using the deal.II library

The potentially significant role of the surface of an elastic body in the overall response of the continuum can be described using the mature theory of surface elasticity. The objective of this contribution is to detail the finite element approximation of the underlying governing equations (both in the volume and on its surface) and their solution using the open-source finite element library deal.II. The fully-nonlinear (geometric and material) setting is considered. The nonlinear problem is solved using a Newton--Raphson procedure wherein the tangent contributions from the volume and surface are computed exactly. The finite element formulation is implemented within the total Lagrangian framework and a Bubnov-Galerkin spatial discretization of the volume and the surface employed. The surface is assumed material. A map between the degrees of freedom on the surface and on the boundary of the volume is used to allocate the contribution from the surface to the global system matrix and residual vector. The deal.II library greatly facilitates the computation of the various surface operators, allowing the numerical implementation to closely match the theory developed in a companion paper. Key features of the theory and the numerical implementation are elucidated using a series of benchmark example problems. The full, documented source code is provided.

preprint2015arXiv

A theory of finite deformation magneto-viscoelasticity

This paper deals with the mathematical modelling of large strain magneto-viscoelastic deformations. Energy dissipation is assumed to occur both due to the mechanical viscoelastic effects as well as the resistance offered by the material to magnetisation. Existence of internal damping mechanisms in the body is considered by decomposing the deformation gradient and the magnetic induction into `elastic' and `viscous' parts. Constitutive laws for material behaviour and evolution equations for the non-equilibrium fields are derived that agree with the laws of thermodynamics. To illustrate the theory the problems of stress relaxation, magnetic field relaxation, time dependent magnetic induction and strain are formulated and solved for a specific form of the constitutive law. The results, that show the effect of several modelling parameters on the deformation and magnetisation process, are illustrated graphically.

preprint2015arXiv

On rate-dependent dissipation effects in electro-elasticity

This paper deals with the mathematical modelling of large strain electro-viscoelastic deformations in electro-active polymers. Energy dissipation is assumed to occur due to mechanical viscoelasticity of the polymer as well as due to time-dependent effective polarisation of the material. Additive decomposition of the electric field $\mathbb{E} = \mathbb{E}_e + \mathbb{E}_v$ and multiplicative decomposition of the deformation gradient $\mathbf{F} = \mathbf{F}_e \mathbf{F}_v$ are proposed to model the internal dissipation mechanisms. The theory is illustrated with some numerical examples in the end.