Source author record

Roger A. Sauer

Roger A. Sauer 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

12works
10topics
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

12 published item(s)

preprint2025arXiv

Efficient snap-to-contact computations for van der Waals interacting fibers

We consider van der Waals interactions between in-plane fibers, where the computational model employs the Lennard-Jones potential and the coarse-grained approach. The involved 6D integral over two interacting fibers is split into a 4D analytical pre-integration over cross sections and the remaining 2D numerical integration along the fibers' axes. Two section-section interaction laws are implemented, refined, and compared. Fibers are modeled using the Bernoulli-Euler beam theory and spatially discretized with isogeometric finite elements. We derive and solve the weak form of both quasi-static and dynamic boundary value problems. Four numerical examples involving highly nonlinear and dynamic snap-to-contact phenomena are scrutinized. We observe that the coarse-graining and pre-integration of interaction potentials enable the efficient modeling of complex phenomena at small length scales.

preprint2022arXiv

A Bayesian regularization-backpropagation neural network model for peeling computations

Bayesian regularization-backpropagation neural network (BR-BPNN) model is employed to predict some aspects of the gecko spatula peeling viz. the variation of the maximum normal and tangential pull-off forces and the resultant force angle at detachment with the peeling angle. K-fold cross validation is used to improve the effectiveness of the model. The input data is taken from finite element (FE) peeling results. The neural network is trained with 75% of the FE dataset. The remaining 25% are utilized to predict the peeling behavior. The training performance is evaluated for every change in the number of hidden layer neurons to determine the optimal network structure. The relative error is calculated to draw a clear comparison between predicted and FE results. It is shown that the BR-BPNN model in conjunction with k-fold technique has significant potential to estimate the peeling behavior.

preprint2022arXiv

An isogeometric finite element formulation for boundary and shell viscoelasticity based on a multiplicative surface deformation split

This work presents a numerical formulation to model isotropic viscoelastic material behavior for membranes and thin shells. The surface and the shell theory are formulated within a curvilinear coordinate system, which allows the representation of general surfaces and deformations. The kinematics follow from Kirchhoff-Love theory and the discretization makes use of isogeometric shape functions. A multiplicative split of the surface deformation gradient is employed, such that an intermediate surface configuration is introduced. The surface metric and curvature of this intermediate configuration follow from the solution of nonlinear evolution laws - ordinary differential equations (ODEs) - that stem from a generalized viscoelastic solid model. The evolution laws are integrated numerically with the implicit Euler scheme and linearized within the Newton-Raphson scheme of the nonlinear finite element framework. The implementation of membrane and bending viscosity is verified with the help of analytical solutions and shows ideal convergence behavior. The chosen numerical examples capture large deformations and typical viscoelasticity behavior, such as creep, relaxation, and strain rate dependence. It is also shown that the proposed formulation can be straightforwardly applied to model boundary viscoelasticity of 3D bodies.

preprint2022arXiv

An isogeometric finite element formulation for frictionless contact of Cosserat rods with unconstrained directors

This paper presents an isogeometric finite element formulation for nonlinear beams with impenetrability constraints, based on the kinematics of Cosserat rods with unconstrained directors. The beam cross-sectional deformation is represented by director vectors of an arbitrary order. For the frictionless lateral beam-to-beam contact, a surface-to-surface contact algorithm combined with an active set strategy and a penalty method is employed. The lateral boundary surface of the beam is parameterized by its axis and cross-sectional boundary curves with NURBS basis functions having at least $C^2$-continuity, which yields a continuous surface metric and curvature for the closest point projection. Three-dimensional constitutive laws of hyperelastic materials are considered. Several numerical examples verify the accuracy and efficiency of the proposed beam contact formulation in comparison to brick element solutions. The lateral contact pressure distribution of the beam formulation is in excellent agreement with the contact pressure of the brick element formulation while requiring much less degrees-of-freedom.

preprint2022arXiv

Nonlinear material identification of heterogeneous isogeometric Kirchhoff-Love shells

This work presents a Finite Element Model Updating inverse methodology for reconstructing heterogeneous material distributions based on an efficient isogeometric shell formulation. It uses nonlinear hyperelastic material models suitable for describing incompressible material behavior as well as initially curved shells. The material distribution is discretized by bilinear elements such that the nodal values are the design variables to be identified. Independent FE analysis and material discretization, as well as flexible incorporation of experimental data, offer high robustness and control. Three elementary test cases and one application example, which exhibit large deformations and different challenges, are considered: uniaxial tension, pure bending, sheet inflation, and abdominal wall pressurization. Experiment-like results are generated from high-resolution simulations with the subsequent addition of up to 4% noise. Local optimization based on the trust-region approach is used. The results show that with a sufficient number of experimental measurements, design variables and analysis elements, the algorithm is capable to reconstruct material distributions with high precision even in the presence of large noise. The proposed formulation is very general, facilitating its extension to other material models, optimization algorithms and meshing approaches. Adapted material discretizations allow for an efficient and accurate reconstruction of material discontinuities by avoiding overfitting due to superfluous design variables. For increased computational efficiency, the analytical sensitivities and Jacobians are provided.

preprint2020arXiv

A modified Coulomb's law for the tangential debonding of osseointegrated implants

Cementless implants are widely used in orthopedic and oral surgery. However, debonding-related failure still occurs at the bone-implant interface. It remains difficult to predict such implant failure since the underlying osseointegration phenomena are still poorly understood. Especially in terms of friction and adhesion at the macro-scale, there is a lack of data and reliable models. The aim of this work is to present a new friction formulation that can model the tangential contact behavior between osseointegrated implants and bone tissue, with focus on debonding. The classical Coulomb's law is combined with a state variable friction law to model a displacement-dependent friction coefficient. A smooth state function, based on the sliding distance, is used to model implant debonding. The formulation is implemented in a 3D nonlinear finite element framework, and it is calibrated with experimental data and compared to an analytical model for mode III cleavage of a coin-shaped, titanium implant (Mathieu et al. 2012). Overall, the results show close agreement with the experimental data, especially the peak and the softening part of the torque curve with a relative error of less than 2.25 %. In addition, better estimates of the bone's shear modulus and the adhesion energy are obtained. The proposed model is particularly suitable to account for partial osseointegration, as is also shown.

preprint2020arXiv

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

We present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff-Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell's mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith's theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires $C^1$-continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-$α$ method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton-Raphson scheme. The interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

preprint2020arXiv

Arbitrary Lagrangian--Eulerian finite element method for curved and deforming surfaces. I. General theory and application to fluid interfaces

An arbitrary Lagrangian--Eulerian (ALE) finite element method for arbitrarily curved and deforming two-dimensional materials and interfaces is presented here. An ALE theory is developed by endowing the surface with a mesh whose in-plane velocity need not depend on the in-plane material velocity, and can be specified arbitrarily. A finite element implementation of the theory is formulated and applied to curved and deforming surfaces with in-plane incompressible flows. Numerical inf--sup instabilities associated with in-plane incompressibility are removed by locally projecting the surface tension onto a discontinuous space of piecewise linear functions. The general isoparametric finite element method, based on an arbitrary surface parametrization with curvilinear coordinates, is tested and validated against several numerical benchmarks. A new physical insight is obtained by applying the ALE developments to cylindrical fluid films, which are computationally and analytically found to be stable to non-axisymmetric perturbations, and unstable with respect to long-wavelength axisymmetric perturbations when their length exceeds their circumference. A Lagrangian scheme is attained as a special case of the ALE formulation. Though unable to model fluid films with sustained shear flows, the Lagrangian scheme is validated by reproducing the cylindrical instability. However, relative to the ALE results, the Lagrangian simulations are found to have spatially unresolved regions with few nodes, and thus larger errors.

preprint2020arXiv

On the presence of a critical detachment angle in gecko spatula peeling -- A numerical investigation using an adhesive friction model

A continuum-based computational contact model is employed to study coupled adhesion and friction in gecko spatulae. Nonlinear finite element analysis is carried out to simulate spatula peeling from a rigid substrate. It is shown that the "frictional adhesion" behavior, until now only observed from seta to toe levels, is also present at the spatula level. It is shown that for sufficiently small spatula pad thickness, the spatula detaches at a constant angle known as the critical detachment angle irrespective of the peeling and shaft angles. The spatula reaches the same energy states at the jump-off contact point, which directly relates to the invariance of the critical detachment angle. This study also reveals that there is an optimum pad thickness associated with the invariance of the critical detachment angle. It is further observed that the sliding of the spatula pad is essential for the invariance of the critical detachment angle.

preprint2016arXiv

A stabilized finite element formulation for liquid shells and its application to lipid bilayers

This paper presents a new finite element (FE) formulation for liquid shells that is based on an explicit, 3D surface discretization using $C^1$-continuous finite elements constructed from NURBS interpolation. Both displacement-based and mixed FE formulations are proposed. The latter is needed for area-incompressible material behavior, where penalty-type regularizations can lead to misleading results. In order to obtain quasi-static solutions, several numerical stabilization schemes are proposed based on either stiffness, viscosity or projection. Several numerical examples are considered in order to illustrate the accuracy and the capabilities of the proposed formulation, and to compare the different stabilization schemes. The presented formulation is capable of simulating non-trivial surface shapes associated with tube formation and protein-induced budding of lipid bilayers. In the latter case, the presented formulation yields non-axisymmetric solutions, which have not been observed in previous simulations. It is shown that those non-axisymmetric shapes are preferred over axisymmetric ones.

preprint2016arXiv

C1-continuous space-time discretization based on Hamilton's law of varying action

We develop a class of C1-continuous time integration methods that are applicable to conservative problems in elastodynamics. These methods are based on Hamilton's law of varying action. From the action of the continuous system we derive a spatially and temporally weak form of the governing equilibrium equations. This expression is first discretized in space, considering standard finite elements. The resulting system is then discretized in time, approximating the displacement by piecewise cubic Hermite shape functions. Within the time domain we thus achieve C1-continuity for the displacement field and C0-continuity for the velocity field. From the discrete virtual action we finally construct a class of one-step schemes. These methods are examined both analytically and numerically. Here, we study both linear and nonlinear systems as well as inherently continuous and discrete structures. In the numerical examples we focus on one-dimensional applications. The provided theory, however, is general and valid also for problems in 2D or 3D. We show that the most favorable candidate -- denoted as p2-scheme -- converges with order four. Thus, especially if high accuracy of the numerical solution is required, this scheme can be more efficient than methods of lower order. It further exhibits, for linear simple problems, properties similar to variational integrators, such as symplecticity. While it remains to be investigated whether symplecticity holds for arbitrary systems, all our numerical results show an excellent long-term energy behavior.

preprint2012arXiv

A computational formulation for constrained solid and liquid membranes considering isogeometric finite elements

A geometrically exact membrane formulation is presented that is based on curvilinear coordinates and isogeometric finite elements, and is suitable for both solid and liquid membranes. The curvilinear coordinate system is used to describe both the theory and the finite element equations of the membrane. In the latter case this avoids the use of local cartesian coordinates at the element level. Consequently, no transformation of derivatives is required. The formulation considers a split of the in-plane and out-of-plane membrane contributions, which allows the construction of a stable formulation for liquid membranes with constant surface tension. The proposed membrane formulation is general, and accounts for dead and live loading, as well as enclosed volume, area, and contact constraints. The new formulation is illustrated by several challenging examples, considering linear and quadratic Lagrange elements, as well as isogeometric elements based on quadratic NURBS and cubic T-splines. It is seen that the isogeometric elements are much more accurate than standard Lagrange elements. The gain is especially large for the liquid membrane formulation since it depends explicitly on the surface curvature.