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 - 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 - 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 -