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 - GEN 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, using the example of UMATs for viscoplastic models, we demonstrate how available user subroutines can be incorporated into the interface while maintaining computational performance of FEniCSx comparable to that of Abaqus. KW - Finite element method KW - Constitutive models KW - FEniCSx KW - UMAT KW - Rust KW - Python KW - C++ PY - 2024 DO - https://doi.org/10.5281/zenodo.13980988 PB - Zenodo CY - Geneva AN - OPUS4-62112 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 - 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 - CONF A1 - Kindrachuk, Vitaliy T1 - Computational investigation of DGG kinetics by phase-field method N2 - Non-equilibrium microstructure evolution in additive manufacturing (AM) is a major barrier for establishing a safe and sustainable application of AM in industrial processes. The constant heat source in the AM continuously affects the just-solidified grains beneath the melting pool, leading to directional grain growth (DGG). While real-time measurements of the non-equilibrium microstructure evolution is challenging, here developing a computational framework to systematically explore DGG becomes imperative. We have advanced a comprehensive approach, integrating mean-field modeling and phase-field simulations, to elucidate the dynamics of DGG under an external driving force. Our simulations unveil a steady-state power-law grain growth kinetics during DGG, characterized by the interplay between curvature-driven dynamics at grain boundary junctions and directional driving forces. T2 - Tagung DGM Additive Fertigung CY - Bremen, Germany DA - 12.06.2024 KW - Directional grain growth KW - Phase-field simulation KW - Additive manufacturing PY - 2024 AN - OPUS4-60750 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 - GEN A1 - Jentzsch, Stefan A1 - Stock, Daniel A1 - Häcker, Ralf A1 - Skrotzki, Birgit A1 - Darvishi Kamachali, Reza A1 - Klingbeil, Dietmar A1 - Kindrachuk, Vitaliy T1 - Split Hopkinson Pressure Bar (SHPB) investigations of steel S355 specimens and complementary characterization methods N2 - The workbook SHPB_S355.xlsx contains in the main spreadsheet S355_TestOverview an overview on the S355 specimens, which were tested in Split Hopkinson Pressure Bar (SHPB) and complementary quasi-static (QS) tests. The tests were conducted with notched hat-shaped specimens without a notch offset (S355_TestOverview/column geometry: “no offset”) and a small notch offset of x=0.35 mm (S355_TestOverview/column geometry: “offset”), cf. Fig. 1 [1]. Furthermore, for selected specimens, links are provided to DIC and bar strain measurement files as well as to evaluations from further characterization methods (microhardness, EBSD). For the boundary conditions at the SHPB projectile impact, the pressures of the compressor p, driving the projectile, and the associated projectile impact velocities are provided in SHPB_S355.xlsx. The DIC displacement measurements are provided in FurtherMeasurements/DIC with frame output times in the file labels, which are associated with an imaginary trigger at the left end of the shortened incident bar (length 300 mm), which is considered within the SHPB simulation setup, see [1]. Furthermore, the DIC reference coordinate systems are provided as COS.jpg files in the respective DIC folders. Starting from the strain signals at the bars, captured by strain gauges at the incident (file name BC_Inc) and transmission bar (BC_Trans), displacement boundary conditions (which are provided for the tests in FurtherMeasurements/BarDisplBCs) are calculated by eq. (12) in [1], incorporating the acoustic velocity equal to 4639 m/s at the bars and a correction factor. Fig. 1 shows the shear specimen geometry (lengths in mm) with the offset of the notches x. For the quasi-static tests force(displacements)-values are directly provided in FurtherMeasurements/Fu_curves, which are considered from the relative displacements of the specimens, evaluated by DIC. The Vickers microhardness (HV 0.01) distribution across the shear localization zone was assessed by QNESS 60A+ EVO (DIN EN ISO 6507-1) for quasi-statically and dynamically tested specimens, applying the small notch offset, such that the specimens did not fail in the localization region. Therefore, raw data is provided in FurtherMeasurements/Microhardness. For fitting the microhardness distribution perpendicularly to the shear bands (as i.e. provided for the SHPB specimen in [1]), representative microhardness profiles were considered. Similarly for specimens with the small notch offset, EBSD data on dynamic and quasi-static tests is provided in FurtherMeasurements/EBSD. The grain size distributions for positions at the notch and the undeformed region are included in the subdirectory ./GrainLists_Shortened and EBSD images in the PowerPoint Presentations. Further information as the sizes and misorientations of the single grains, is incorporated in the subfolders ./GrainLists_FurtherInf. KW - Adiabatic shear bands KW - Split Hopkinson pressure bar KW - Digital image correlation PY - 2025 DO - https://doi.org/10.5281/zenodo.17591439 PB - Zenodo CY - Geneva AN - OPUS4-64842 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 -