publications
I see no reason to address the - in any case erroneous - comments of your anonymous expert. — Albert Einstein
preprints
2024
- Fractal structure, depinning, and hysteresis of dislocations in high-entropy alloysHoa Thi Le, Wolfram G Nöhring, and Lars PastewkaarXiv:2410.21838 (2024)
High-entropy alloys (HEAs) are complex alloys containing multiple elements in high concentrations. Plasticity in HEAs is carried by dislocations, but the random nature of their composition pins dislocations, effectively hindering their motion. We investigate the resulting complex structure of the dislocation in terms of spatial correlation functions, which allow us to draw conclusions on the fractal geometry of the dislocation. At high temperature, where thermal fluctuations dominate, dislocations adopt the structure of a random walk with Hurst exponent 1/2 or fractal dimension 3/2. At low temperature we find larger Hurst exponents (lower dimensions), with a crossover to an uncorrelated structure beyond a correlation length. These changes in structure are accompanied by an emergence of hysteresis (and hence pinning) in the motion of the dislocation at low temperature. We use a modified Labusch/Edwards-Wilkinson-model to argue that this correlation length must be an intrinsic property of the HEA. This means dislocations in HEAs are an individual pinning limit, where segments of the dislocation are independently pinned by local distortions of the crystal lattice that are induced by chemical heterogeneity.
2022
- Crack-path selection in phase-field models for brittle fractureW Beck Andrews, and Lars PastewkaarXiv:2203.16467 (2022)
This work presents a critical overview of the effects of different aspects of model formulation on crack path selection in quasi-static phase field fracture. We consider different evolution methods, mechanics formulations, fracture dissipation energy formulations, and forms of the irreversibility condition. The different model variants are implemented with common numerical methods based on staggered solution of the phase-field and mechanics sub-problems via FFT-based solvers. These methods mix standard approaches with novel elements, such as the use of bound-constrained conjugate gradients for the phase field sub-problem and a heuristic method for near-equilibrium evolution. We examine differences in crack paths between model variants in simple model systems and microstructures with randomly heterogeneous Young’s modulus. Our results indicate that near-equilibrium evolution methods are preferable for quasi-static fracture of heterogeneous microstructures compared to minimization and time-dependent methods. In examining mechanics formulations, we find distinct effects of crack driving force and the model for contact implicit in phase field fracture. Our results favor the use of a strain-spectral decomposition for the crack driving force but not the contact model. Irreversibility condition and fracture dissipation energy formulation were also found to affect crack path selection, but systematic effects were difficult to deduce due to the overall sensitivity of crack selection within the heterogeneous microstructures. Our findings support the use of the AT1 model over the AT2 model and irreversibility of the phase field within a crack set rather than the entire domain. Sensitivity to these differences in formulation was reduced but not eliminated by reducing the crack width parameter \ell relative to the size scale of the random microstructures.
2021
- Efficient topology optimization using compatibility projection in micromechanical homogenizationIndre Jödicke, Richard J. Leute, Till Junge, and Lars PastewkaarXiv:2107.04123 (2021)
The adjoint method allows efficient calculation of the gradient with respect to the design variables of a topology optimization problem. This method is almost exclusively used in combination with traditional Finite-Element-Analysis, whereas Fourier-based solvers have recently shown large efficiency gains for homogenization problems. In this paper, we derive the discrete adjoint method for Fourier-based solvers that employ compatibility projection. We demonstrate the method on the optimization of composite materials and auxetic metamaterials, where void regions are modelled with zero stiffness.