Search

In this paper, a novel chemo-mechanical model is proposed for the description of the stretch-dependent chemical processes known as Bayliss effect and their impact on the active contraction in vascular smooth muscle. These processes are responsible for the adaptive reaction of arterial walls to changing blood pressure by which the blood vessels actively support the heart in providing sufficient blood supply for varying demands in the supplied tissues. The model is designed to describe two different stretch-dependent mechanisms observed in smooth muscle cells (SMCs): a calcium-dependent and a calcium-independent contraction. For the first one, stretch of the SMCs leads to an inlet of calcium ions which activates the myosin light chain kinase (MLCK). The increased activity of MLCK triggers the contractile units of the cells resulting in the contraction on a comparatively short time scale. For the calcium-independent contraction mechanism, stretch-dependent receptors of the cell membrane stimulate an intracellular reaction leading to an inhibition of the antagonist of MLCK, the myosin light chain phosphatase resulting in a contraction on a comparatively long time scale. An algorithmic framework for the implementation of the model in finite element programs is derived. Based thereon, it is shown that the proposed approach agrees well with experimental data. Furthermore, the individual aspects of the model are analyzed in numerical simulations of idealized arteries subject to internal pressure waves with changing intensities. The simulations show that the proposed model is able to describe the experimentally observed contraction of the artery as a reaction to increased internal pressure, which can be considered a crucial aspect of the regulatory mechanism of muscular arteries.

The brief history of relaxation in continuum mechanics ranges from early application of non-convex plasticity and phase transition formulations to small and large strain continuum damage mechanics. However, relaxed continuum damage mechanics formulations are still limited in the following sense that their material response lack to model strain softening and the convexification of the non-convex incremental stress potential is computationally costly. This paper presents a reduced model for relaxed continuum damage mechanics at finite strains which includes strain softening by a fiber-specific damage in the microsphere approach. Computational efficiency is achieved by novel adaptive algorithms for the fast convexification of the one-dimensional fiber material model. The algorithms are benchmarked against state-of-the-art methods, and the choice of quadrature schemes for the microsphere approach is discussed. This contribution is finalized by a mesh independence test.

In this paper, a framework for the simulation of crack propagation in brittle and ductile materials is proposed. The framework is derived by extending the eigenerosion approach of Pandolfi and Ortiz (Int J Numer Methods Eng 92(8):694–714, 2012. 10.1002/nme.4352) to finite strains and by connecting it with a generalized energy-based, Griffith-type failure criterion for ductile fracture. To model the elasto-plastic response, a classical finite strain formulation is extended by viscous regularization to account for the shear band localization prior to fracture. The compression–tension asymmetry, which becomes particularly important during crack propagation under cyclic loading, is incorporated by splitting the strain energy density into a tensile and compression part. In a comparative study based on benchmark problems, it is shown that the unified approach is indeed able to represent brittle and ductile fracture at finite strains and to ensure converging, mesh-independent solutions. Furthermore, the proposed approach is analyzed for cyclic loading, and it is shown that classical Wöhler curves can be represented.

The modeling of damage processes in materials constitutes an ill-posed mathematical problem which manifests in mesh-dependent finite element results. The loss of ellipticity of the discrete system of equations is counteracted by regularization schemes of which the gradient enhancement of the strain energy density is often used. In this contribution, we present an extension of the efficient numerical treatment, which has been proposed in [1], to materials that are subjected to large deformations. Along with the model derivation, we present a technique for element erosion in the case of severely damaged materials. Efficiency and robustness of our approach is demonstrated by two numerical examples.

The mechanized tunnel construction is carried out by tunnel boring machines, in which the soil in front of the working face is removed, and the tunnel lining is carried out with shotcrete or the setting of segments and their back injection. Advancements in this field aim towards increase of the excavation efficiency and increase of the tool lifetime, especially in rock-dominated grounds. The latter is achieved by understanding the wear mechanisms abrasion and surface-fatigue, and by knowledge of the microstructure-property relation of the utilized materials. Improvements for tool concepts are derived, based on experiments and simulations. A key parameter towards efficient rock excavation is the shape of the cutting edge of the utilized disc cutters. Sharp cutting edges have proven to generate higher rock excavation rates compared to blunt ones. The compressive strength of the utilized steel has to be high, to inhibit plastic deformation and thereby to maintain sharp cutting edges. This requirement competes with the demand for toughness, which is necessary to avoid crack-growth in the case of cyclic loading. Solutions for this contradiction lie in specially designed multiphase microstructures, containing both hard particles and ductile microstructural constituents. Besides adapting the alloying concept, these required microstructures and the associated properties can be adjusted by specific heat-treatments.

An established strategy for material modeling is provided by energy-based principles such that evolution equations in terms of ordinary differential equations can be derived. However, there exist a variety of material models that also need to take into account non-local effects to capture microstructure evolution. In this case, the evolution of microstructure is described by a partial differential equation. In this contribution, we present how Hamilton’s principle provides a physically sound strategy for the derivation of transient field equations for all state variables. Therefore, we begin with a demonstration how Hamilton’s principle generalizes the principle of stationary action for rigid bodies. Furthermore, we show that the basic idea behind Hamilton’s principle is not restricted to isothermal mechanical processes. In contrast, we propose an extended Hamilton principle which is applicable to coupled problems and dissipative microstructure evolution. As example, we demonstrate how the field equations for all state variables for thermo-mechanically coupled problems, i.e., displacements, temperature, and internal variables, result from the stationarity of the extended Hamilton functional. The relation to other principles, as the principle of virtual work and Onsager’s principle, is given. Finally, exemplary material models demonstrate how to use the extended Hamilton principle for thermo-mechanically coupled elastic, gradient-enhanced, rate-dependent, and rate-independent materials.

In this paper we present a fully-coupled, two-scale homogenization method for dynamic loading in the spirit of FE$$^2$$ 2 methods. The framework considers the balance of linear momentum including inertia at the microscale to capture possible dynamic effects arising from micro heterogeneities. A finite-strain formulation is adapted to account for geometrical nonlinearities enabling the study of e.g. plasticity or fiber pullout, which may be associated with large deformations. A consistent kinematic scale link is established as displacement constraint on the whole representative volume element. The consistent macroscopic material tangent moduli are derived including micro inertia in closed form. These can easily be calculated with a loop over all microscopic finite elements, only applying existing assembly and solving procedures. Thus, making it suitable for standard finite element program architectures. Numerical examples of a layered periodic material are presented and compared to direct numerical simulations to demonstrate the capability of the proposed framework. In addition, a simulation of a split Hopkinson tension test showcases the applicability of the framework to engineering problems.

Im Rahmen des DFG‐Graduiertenkollegs GRK 2250/1 „Impaktsicherheit von Baukonstruktionen durch mineralisch gebundene Komposite“ werden neuartige Verstarkungsmaterialien mit unterschiedlichen Faserbewehrungen entwickelt, die die Widerstandsfahigkeit bestehender Konstruktionen gegen kurzzeitdynamische Belastungen deutlich erhohen sollen. Des Weiteren werden adaquate Prufkonzepte und Messverfahren fur eine grundlegende Charakterisierung und Beschreibung des dehnratenabhangigen Materialverhaltens geschaffen. Die experimentell erzielten Ergebnisse stellen eine wichtige Grundlage fur die Formulierung numerischer Modelle zur Simulation von impaktbeanspruchten Kompositen sowie von verstarkten Betonstrukturen mit Kopplung unterschiedlicher Raum‐ und Zeitskalen dar. Schlieslich werden im Rahmen des Projekts Methoden fur die Analyse und Bewertung der Nachhaltigkeit und Resilienz von Verstarkungsmasnahmen mit neuen Verbundwerkstoffen erarbeitet.

Tissue degradation plays a crucial role in vascular diseases such as atherosclerosis and aneurysms. Computational modeling of vascular hemodynamics incorporating both arterial wall mechanics and tissue degradation has been a challenging task. In this study, we propose a novel finite element method-based approach to model the microscopic degradation of arterial walls and its interaction with blood flow. The model is applied to study the combined effects of pulsatile flow and tissue degradation on the deformation and intra-aneurysm hemodynamics. Our computational analysis reveals that tissue degradation leads to a weakening of the aneurysmal wall, which manifests itself in a larger deformation and a smaller von Mises stress. Moreover, simulation results for different heart rates, blood pressures and aneurysm geometries indicate consistently that, upon tissue degradation, wall shear stress increases near the flow-impingement region and decreases away from it. These findings are discussed in the context of recent reports regarding the role of both high and low wall shear stress for the progression and rupture of aneurysms.

In this paper, we propose a novel, semi-analytic approach for the two-scale, computational modeling of concentration transport in packed bed reactors. Within the reactor, catalytic pellets are stacked, which alter the concentration evolution. Firstly, the considered experimental setup is discussed and a naive one-scale approach is presented. This one-scale model motivates, due to unphysical fitted values, to enrich the computational procedure by another scale. The computations on the second scale, here referred to as microscale, are based on a proper investigation of the diffusion process in the catalytic pellets from which, after continuum-consistent considerations, a sink term for the macroscopic advection–diffusion–reaction process can be identified. For the special case of a spherical catalyst pellet, the parabolic partial differential equation at the microscale can be reduced to a single ordinary differential equation in time through a semi-analytic approach. After the presentation of our model, we show results for its calibration against the macroscopic response of a simple standard mass transport experiment. Based thereon, the effective diffusion parameters of the catalyst pellets can be identified.

A method to quantify uncertain macroscopic material properties resulting from variations of a material’s microstructure morphology is proposed. Basis is the computational homogenization of virtual experiments as part of a Monte-Carlo simulation to obtain the associated uncertain macroscopic material properties. A new general approach is presented to construct a set of artificial microstructures, which exhibits a statistically similar variation of the morphology as the real material’s microstructure. The individual artificial microstructures are directly constructed in a way that a lower discretization effort is required compared to real microstructures. The costs to perform the computational homogenization for all considered SSVEs are reduced by an adapted form of the Finite Cell concept and by applying the multilevel Monte-Carlo method. As an illustrative example, the proposed method is applied to a real Dual-Phase steel microstructure.

This paper proposes a new method for in vivo and almost real-time identification of biomechanical properties of the human cornea based on non-contact tonometer data. Further goal is to demonstrate the method’s functionality based on synthetic data serving as reference. For this purpose, a finite element model of the human eye is constructed to synthetically generate full-field displacements from different data sets with keratoconus-like degradations. Then, a new approach based on the equilibrium gap method combined with a mechanical morphing approach is proposed and used to identify the material parameters from virtual test data sets. In a further step, random absolute noise is added to the virtual test data to investigate the sensitivity of the new approach to noise. As a result, the proposed method shows a relevant accuracy in identifying material parameters based on full-field displacements. At the same time, the method turns out to work almost in real time (order of a few minutes on a regular workstation) and is thus much faster than inverse problems solved by typical forward approaches. On the other hand, the method shows a noticeable sensitivity to rather small noise amplitudes rendering the method not accurate enough for the precise identification of individual parameter values. However, analysis show that the accuracy is sufficient for the identification of property ranges which might be related to diseased tissues. Thereby, the proposed approach turns out promising with view to diagnostic purposes.