TY - JOUR A1 - Rosenbusch, Sjard Mathis A1 - Diercks, Philipp A1 - Kindrachuk, Vitaliy A1 - Unger, Jörg F. T1 - Integrating custom constitutive models into FEniCSx: A versatile approach and case studies N2 - The development and integration of user-defined constitutive relationships into finite element (FE) tools using standardized interfaces play a pivotal role in advancing the capabilities of FE solvers for structural mechanics applications. While commercial FE solvers like Abaqus and Ansys have designed their interfaces to provide custom stresses, tangents, and updated history variables, the open-source solver FEniCSx remains efficient only when the constitutive update has an analytical representation. This restricts the application of FEniCSx for non-linear structural mechanics. Since FEniCSx has become a powerful and popular open-source tool for solving partial differential equations, particularly due to its automatic computation of Hessians, we aim to develop a generalized interface to enhance its capability for constitutive modeling. This approach will address complex constitutive equations that require iterative solutions at the quadrature point level. Specific implementation challenges, such as using return-mapping procedures, can then be managed commonly. The provided interface for custom constitutive models offers a versatile way to implement them in various languages, including C++, Python, Rust, and Fortran. Finally, with UMATs for viscoplastic models as an example, we demonstrate how existing user subroutines can be incorporated into the interface and utilized within the FEniCSx framework. KW - Finite element method KW - Constitutive models KW - FEniCSx KW - UMAT KW - Rust KW - Python KW - C++ PY - 2025 UR - https://nbn-resolving.org/urn:nbn:de:kobv:b43-630439 DO - https://doi.org/10.1016/j.advengsoft.2025.103922 SN - 0965-9978 VL - 206 SP - 1 EP - 11 PB - Elsevier CY - Amsterdam AN - OPUS4-63043 LA - eng AD - Bundesanstalt fuer Materialforschung und -pruefung (BAM), Berlin, Germany ER - TY - JOUR A1 - Kindrachuk, Vitaliy A1 - Klunker, Andre T1 - Phase field modeling of Hertzian cone cracks under spherica indentation N2 - A phase field model of brittle fracture has been developed to simulate the Hertzian crack induced by penetration of a rigid sphere to an isotropic linear-elastic half-space. The fracture formation is regarded as a diffusive field variable, which is zero for the intact material and unity if there is a crack. Crack growth is assumed to be driven by a strain invariant. The numerical implementation is performed with the finite element method and an implicit time integration scheme. The mechanical equilibrium and the phase field equations are solved in a staggered manner, sequentially updating the displacement field and the phase field variable. Numerical examples demonstrate the capability of the model to reproduce the nucleation and growth of the Hertzian cone crack. KW - Hertzian cracks KW - Phase field model KW - Contact mechanics PY - 2021 DO - https://doi.org/10.1007/s11223-021-00251-9 VL - 52 IS - 6 SP - 967 EP - 974 AN - OPUS4-52276 LA - eng AD - Bundesanstalt fuer Materialforschung und -pruefung (BAM), Berlin, Germany ER - TY - JOUR A1 - Kindrachuk, Vitaliy A1 - Darvishi Kamachali, Reza T1 - Mean-field modeling and phase-field simulation of grain growth under directional driving forces N2 - Directional grain growth is a common phenomenon in the synthetic and natural evolution of various polycrystals. It occurs in the presence of an external driving force, such as a temperature gradient, along which grains show a preferred, yet competitive, growth. Novel additive manufacturing processes, with intense, localized energy deposition, are prominent examples of when directional grain growth can occur, beneath the melting pool. In this work, we derive a phenomenological mean-field model and perform 3D phase-field simulations to investigate the directional grain growth and its underlying physical mechanisms. The effect of the intensity of driving force is simulated and systematically analyzed at the evolving growth front as well as various cross-sections perpendicular to the direction of the driving force. We found that although the directional growth significantly deviates from normal grain growth, it is still governed by a power law relation α tⁿ with an exponent n ~ 0.6–0.7. The exponent n exhibits a nontrivial dependence on the magnitude of the directional driving force, such that the lowest growth exponent is observed for intermediate driving forces. We elaborate that this can originate from the fact that the forces at grain boundary junctions evolve out of balance under the influence of the directional driving force. With increasing the driving forces, the growth exponent asymptotically approaches a value of n ≈ 0.63, imposed by the largest possible grain aspect ratio for given grain boundary energies. The current combined mean-field and phase-field framework pave the way for future exploration in broader contexts such as the evolution of complex additively manufactured microstructures. KW - Additive manufacturing KW - Phase-field simulation KW - Grain growth KW - Mean-field modelling KW - Directional grain growth PY - 2024 UR - https://nbn-resolving.org/urn:nbn:de:kobv:b43-593210 DO - https://doi.org/10.1016/j.mtla.2023.101989 SN - 2589-1529 VL - 33 SP - 1 EP - 10 PB - Elsevier AN - OPUS4-59321 LA - eng AD - Bundesanstalt fuer Materialforschung und -pruefung (BAM), Berlin, Germany ER - TY - JOUR A1 - Jentzsch, Stefan A1 - Stock, Daniel A1 - Häcker, Ralf A1 - Skrotzki, Birgit A1 - Darvishi Kamachali, Reza A1 - Klingbeil, Dietmar A1 - Kindrachuk, Vitaliy T1 - Shear Band Formation with Split Hopkinson Bar Experiments N2 - The essence of dynamic failure is closely linked to dramatic shear deformations which often lead to the formation of adiabatic shear bands (ASB). Under high loading velocities and the subsequent rapid temperature increase, the localization of shear strain is crucial in view of safety issues of systems in mechanical and aircraft engineering, especially with respect to fast rotating components and diverse crash scenarios. In this research, we perform high speed impact tests at the split Hopkinson pressure bar (SHPB) setup and use particular hat-shaped specimen geometries that resemble the stresses and failure conditions at the component level. In the first step, we specify a notched specimen geometry using finite element (FE) simulations to ensure pure shear. Further, quasi-static compressive tests and a series of impact tests at high strain rates of 10^3-10^4 s^-1 are conducted on specimens manufactured from a fine-grain structural steel with the properties of S355. Optical microscopy and electron backscatter diffraction (EBSD) of the sheared zones unveil significant localization to maximal shear strains of about 0.9 accompanied by grain refinement by factors 5 to 14. The displacements across the surface of the specimens are captured with subset-based local digital image correlation (DIC) during the impact time, and serve as an objective to validate a viscoplastic constitutive relationship. More precisely, the deformation distribution is accurately reproduced by the widely recognized Johnson-Cook (JC) model, which features an enhanced description of damage evolution. Thus, combining experimental and characterization techniques, continuum mechanics and reasonable optimization strategies for the identification of model parameters provides an efficient approach for comprehensive insights into the strain localization behaviour and its impact on the mechanical performance of S355 under extreme strain rates and deformations. KW - Adiabatic shear bands KW - Finite element analysis KW - Viscoplastic material modelling KW - Split Hopkinson pressure bar KW - Hat-shaped specimen KW - Johnson–Cook parameter identification PY - 2024 UR - https://nbn-resolving.org/urn:nbn:de:kobv:b43-613339 DO - https://doi.org/10.1016/j.ijmecsci.2024.109749 VL - 284 SP - 1 EP - 14 PB - Elsevier BV AN - OPUS4-61333 LA - eng AD - Bundesanstalt fuer Materialforschung und -pruefung (BAM), Berlin, Germany ER - TY - JOUR A1 - Dub, S. A1 - Haftaoglu, Cetin A1 - Kindrachuk, Vitaliy T1 - Estimate of theoretical shear strength of C60 single crystal by nanoindentation N2 - The onset of plasticity in a single crystal C60 fullerite was investigated by nanoindentation on the (111) crystallographic plane. The transition from elastic to plastic deformation in a contact was observed as pop-in events on loading curves. The respective resolved shear stresses were computed for the octahedral slip systems ⟨011¯¯¯⟩{111}, supposing that their activation resulted in the onset of plasticity. A finite element analysis was applied, which reproduced the elastic loading until the first pop-in, using a realistic geometry of the Berkovich indenter blunt tip. The obtained estimate of the C60 theoretical shear strength was about 1/11 of the shear modulus on {111} planes. KW - Finite element analysis KW - Fullerite KW - Nanoindentation PY - 2021 UR - https://nbn-resolving.org/urn:nbn:de:kobv:b43-523208 DO - https://doi.org/10.1007/s10853-021-05991-2 VL - 56 IS - 18 SP - 10905 EP - 10914 PB - Springer Nature AN - OPUS4-52320 LA - eng AD - Bundesanstalt fuer Materialforschung und -pruefung (BAM), Berlin, Germany ER -