Materials Theory最新文献

Crack initiation emerges due to a combination of elasticity, plasticity, and disorder, and it displays strong dependence on the material’s microstructural details. The characterization of the structural uncertainty in the original microstructure is typically empirical and systematic characterization protocols are lacking. In this paper, we propose an investigational tool in the form of the curvature of an ellipsoidal notch: As the radius of curvature at the notch increases, there is a dynamic phase transition from notch-induced crack initiation to disorder-induced crack nucleation. We argue that the this transition may unveil the characteristic length scale of structural disorder in the material. We investigate brittle but elastoplastic metals with continuum, microstructural disorder that could originate in a manufacturing process, such as alloying. We perform extensive and realistic simulations, using a phase-field approach coupled to crystal plasticity, where microstructural disorder and notch width are systematically varied. We identify the brittle-to-quasi-brittle transition for various disorder strengths in terms of the damage and stress evolution. Moreover, we investigate precursors to crack initiation that we quantify in terms of the expected stress drops during displacement control loading.

Macroscopic properties of structural materials are strongly dependent on their microstructure. However, the modeling of their evolution is a complex task because of the mechanisms involved such as plasticity, recrystallization, and phase transformations, which are common processes taking place in metallic alloys. This complexity led to a growing interest in atomistic simulations formulated without any auxiliary hypotheses beyond the choice of interatomic potential. In this context, we propose here a model based on an overdamped stochastic evolution of particles interacting through inter-atomic forces. The model settles to the correct thermal equilibrium distribution in canonical and grand-canonical ensembles and is used to study the grain boundary migration. Finally, a comparison of our results with those obtained by molecular dynamics shows that our approach reproduces the complex atomic-scale dynamics of grain boundary migration correctly.

We derive the Green tensor of Mindlin’s anisotropic first strain gradient elasticity. The Green tensor is valid for arbitrary anisotropic materials, with up to 21 elastic constants and 171 gradient elastic constants in the general case of triclinic media. In contrast to its classical counterpart, the Green tensor is non-singular at the origin, and it converges to the classical tensor a few characteristic lengths away from the origin. Therefore, the Green tensor of Mindlin’s first strain gradient elasticity can be regarded as a physical regularization of the classical anisotropic Green tensor. The isotropic Green tensor and other special cases are recovered as particular instances of the general anisotropic result. The Green tensor is implemented numerically and applied to the Kelvin problem with elastic constants determined from interatomic potentials. Results are compared to molecular statics calculations carried out with the same potentials.

A mesoscale model is introduced to study the dynamics of material defects lying at interface junctions. The proposed framework couples the dynamics of discrete dislocation and disclination lines. Disclinations are expected to be natural defects at interface junctions; their presence serving the purpose of accommodating discontinuities in rotation fields at material interface junctions. Crystallography-based rules are proposed to describe the kinematics of disclination motion. A discrete-continuous couple-stress framework, in which discrete defect lines are introduced as plastic eigenstrains and eigencurvatures, is proposed to explicitly follow the dynamics of interfacial defects. The framework is then applied to study (left (10bar {1}2right)) twin transverse propagation and thickening in magnesium. Focusing first on the case of a twin domain, It is shown that a disclination based representation of twin domains allows for an appropriate mechanistic description of the kinematics of shear transformations. In what concerns twin thickening, the stability of defects at twin interfaces is further studied. To this end, a 3D crater lying on a twin interface is described as a dipole of disclination loops. Upon self-relaxation, it is found that out of plane motion of disclinations followed by the nucleation of twinning dislocations can be activated; thereby showing that conservative non-planar motion of disclinations can be thermodynamically favorable; mechanism that had been postulated some 50 years ago.

An electrochemical cell contains three open thermodynamic systems that, in dynamic equilibrium, equalize their electrochemical potentials with that of their surrounding by forming an electric-double-layer-capacitor at the interface of the electrolyte with each of the two electrodes. Since the electrode/electrolyte interfaces are heterojunctions, the electrochemical potentials or Fermi levels of the two materials that contact the electrolyte at the two electrodes determine the voltage of a cell. The voltage is the sum of the voltages of the two interfacial electric-double-layer capacitors at the two electrode/electrolyte interfaces. A theoretical analysis of the thermodynamics that gives a quantitative prediction of the observed voltages in an asymmetric cell with an S₈ relay at the positive electrode is provided. In addition, new discharge data and an X-ray photoelectron spectroscopy analysis of the lithium plated on the positive electrode of a discharged cell is presented. Ab initio, DFT methods were used to calculate the band structure and surface-state energies of the crystalline S₈ solid sulfur relay. The theoretical exposition of the thermodynamics of the operative driving force of the chemical reactions in an electrochemical cell demonstrate that our initial experimental data and conclusions are valid. Other reported observations of lithium plating on the positive electrode, observations that were neither exploited nor their origins specified, are also cited.

We explore new ways of regulating defect behavior in systems of conservation laws. Contrary to usual regularization schemes (such as a vanishing viscosity limit), which attempt to control defects by making them smoother, our schemes result in defects which are more singular, and we thus refer to such schemes as “irregularizations”. In particular, we seek to produce delta shock defects which satisfy a condition of stationarity. We are motivated to pursue such exotic defects by a physical example arising from dislocation dynamics in materials physics, which we describe.

An exact analytical solution of the Fick’s second law was developed and applied to the mixed-mode growth of a multicomponent ellipsoidal precipitate growing with constant eccentricities in the quasi-stationary regime. The solution is exact if the nominal composition, equilibrium concentrations and material properties are assumed constant, and can be applied to compounds having no limitations in the number of components. The solution was compared to the solution calculated by a diffusion-controlled application software and it was found that the solute concentrations at the interface can be determined knowing only the nominal composition, the full equilibrium concentrations and the coefficients of diffusion. The thermodynamic calculations owing to find alternative tie-lines are proven to be useless in the mixed-mode model. From this, it appears that the search of alternative tie-lines is computationally counterproductive, even when the interface has a very high mobility. A more efficient computational scheme is possible by considering that a moving interface is not at equilibrium.

Continuum dislocation dynamics (CDD) aims at representing the evolution of systems of curved and connected dislocation lines in terms of density-like field variables. Here we discuss how the processes of dislocation multiplication and annihilation can be described within such a framework. We show that both processes are associated with changes in the volume density of dislocation loops: dislocation annihilation needs to be envisaged in terms of the merging of dislocation loops, while conversely dislocation multiplication is associated with the generation of new loops. Both findings point towards the importance of including the volume density of loops (or ’curvature density’) as an additional field variable into continuum models of dislocation density evolution. We explicitly show how this density is affected by loop mergers and loop generation. The equations which result for the lowest order CDD theory allow us, after spatial averaging and under the assumption of unidirectional deformation, to recover the classical theory of Kocks and Mecking for the early stages of work hardening.

A formal asymptotic analysis of two classes of phase field models for void growth and coarsening in irradiated solids has been performed to assess their sharp-interface kinetics. It was found that the sharp interface limit of type B models, which include only point defect concentrations as order parameters governed by Cahn-Hilliard equations, captures diffusion-controlled kinetics. It was also found that a type B model reduces to a generalized one-sided classical Stefan problem in the case of a high driving thermodynamic force associated with the void growth stage, while it reduces to a generalized one-sided Mullins-Sekerka problem when the driving force is low in the case of void coarsening. The latter case corresponds to the famous rate theory description of void growth. Type C models, which include point defect concentrations and a non-conserved order parameter to distinguish between the void and solid phases and employ coupled Cahn-Hilliard and Allen-Cahn equations, are shown to represent mixed diffusion and interfacial kinetics. In particular, the Allen-Cahn equation of model C reduces to an interfacial constitutive law representing the attachment and emission kinetics of point defects at the void surface. In the limit of a high driving force associated with the void growth stage, a type C model reduces to a generalized one-sided Stefan problem with kinetic drag. In the limit of low driving forces characterizing the void coarsening stage, however, the model reduces to a generalized one-sided Mullins-Sekerka problem with kinetic drag. The analysis presented here paves the way for constructing quantitative phase field models for the irradiation-driven nucleation and growth of voids in crystalline solids by matching these models to a recently developed sharp interface theory.

We present results of a two-scale model of disordered cellular materials where we describe the microstructure in an idealized manner using a beam network model and then make a transition to a Cosserat-type continuum model describing the same material on the macroscopic scale. In such scale transitions, normally either bottom-up homogenization approaches or top-down reverse modeling strategies are used in order to match the macro-scale Cosserat continuum to the micro-scale beam network. Here we use a different approach that is based on an energetically consistent continuization scheme that uses data from the beam network model in order to determine continuous stress and strain variables in a set of control volumes defined on the scale of the individual microstructure elements (cells) in such a manner that they form a continuous tessellation of the material domain. Stresses and strains are determined independently in all control volumes, and constitutive parameters are obtained from the ensemble of control volume data using a least-square error criterion. We show that this approach yields material parameters that are for regular honeycomb structures in close agreement with analytical results. For strongly disordered cellular structures, the thus parametrized Cosserat continuum produces results that reproduce the behavior of the micro-scale beam models both in view of the observed strain patterns and in view of the macroscopic response, including its size dependence.