The Maugis analysis is applied to adhesive contact between a cylinder with various wave profiles and a semi-infinite, elastic half-plane. We extend the analysis of Waters, Lee and Guduru, who consider the adhesive contact of a Hertzian indenter on a semi-infinite, elastic half-space with axi-symmetric, wave profiles. This work gives the closed-form contact mechanical solution for continuous, line contact without the need for any approximation. The resulting semi-analytical model serves to complement existing (numerical) models of adhesive line contact with the static load-area response as a reference. Herewith we analyse adhesion-induced loading-unloading hysteresis and contrast semi-analytical and numerical result to assess the limit of the former analysis. We confirm that roughness-induced dissipation vanishes with increasing wave roughness and decreasing Maugis parameter due to an increase in the range of adhesion and cavitation. Instability and cavitation are mutually exclusive at a given load-area locus yet occur successively in the same contact. An interesting result is that the Johnson parameter, that is known to govern the amplification of adhesion in the JKR-limit, bounds the load-area envelope irrespective of Maugis parameter. However, the Johnson parameter does not control the occurrence of roughness-induced dissipation and thus interface toughening.

Viscoelasticity and roughness are among the possible causes of the adhesive hysteresis displayed by soft contacts. Viscoelasticity causes an increased effective work of adhesion due to stiffening of the contact, while roughness is responsible for elastic instabilities. Herein, we explore the interplay between viscoelasticity and roughness by simulating in two dimensions the retraction of a rigid cylinder, with wavy surface profile, from a viscoelastic half-space. The wave amplitude and length are varied to induce instabilities in the load–to–area response, while the retraction velocity is increased to promote viscoelasticity. Results show that, in the regime where viscoelasticity is confined to the edges of the wavy contact, the contributions of viscoelasticity and waviness to adhesive hysteresis are nearly independent and additive. At low retraction rates, the instabilities in the load-area curve typical of rough elastic contacts are suppressed by viscoelasticity: the contact stiffens to promote a stable decrease of the contact area with load. This occurs with a minimal change in work of adhesion. However, when the instantaneous limit is met at high retraction rates, mechanical instabilities appear.

Modelling dislocation glide over the initial part of a stress–strain curve of metals received little attention up to now. However, dislocation glide is essential to ones understanding of the fundamental relationship between inelastic deformation and the evolution of the dislocation network structure. Therefore, we present a model of dislocation-driven deformation under static loading conditions. We reproduce repeated cyclic uniaxial tensile tests on Interstitial-Free and Low-Alloy steels. The elastic mechanical behaviour is described by isotropic linear elasticity, pre-yield anelastic mechanical behaviour by a dislocation bow-out model with dissipation, and the post-yield evolution of dislocation network structure by a statistical storage model. We hypothesise that when the local anelastic compliance is lower than the global plastic compliance, deformation is mechanically recoverable, and vice versa. This hypothesis is corroborated with the classical Taylor relation. We report the relation between stable and unstable dislocation glide using this prototypical modelling framework. We find four structural variables, that are based on dislocation physics, to describe the stress–strain curve: total dislocation density, average dislocation segment length, dislocation junction formation rate, and average dislocation junction length. Firstly, we quantify the dislocation network evolution during uniaxial monotonic loading, and verify work-hardening by dislocation junction formation and a Taylor-type equation for flow. Finally, we present a semi-empirical relation for the evolution of the dislocation network structure. Which allows us to: refine the physical interpretation of the Taylor relationship, and rationalise experimental observations on apparent modulus degradation by thermomechanical processing. Both these findings circumvent the limitations of current, physics-based hardening models.

The volume increase and shape change during austenite to martensite transformation in dual-phase (DP) steels are largely accommodated in the microstructure by the deformation of the surrounding ferrite matrix. Accurate estimation of transformation-induced deformation of ferrite via experiments and modeling is essential for predicting the subsequent mechanical behavior of DP steels. This study aims to illustrate the disadvantages of simplifying the anisotropic transformation deformation of martensite to isotropic dilatation for modeling the transformation-induced deformation of ferrite. A novel methodology is developed comprising sequential experimental and numerical research on DP steels to quantify transformation-induced strains in ferrite. This methodology combines the results of prior austenite grain reconstruction, phenomenological theory of martensite crystallography and electron backscatter diffraction (EBSD) orientation data to estimate variant-specific transformation deformation. Subsequently, by comparison of full-field micromechanical calculation results on a virtual DP steel microstructure with experimental EBSD kernel average misorientation and geometrically necessary dislocation measurement results it is shown that neglecting the shear deformation associated with the martensitic transformation leads to significant underestimation in the prediction of transformation-induced strains in ferrite.

Predicting microstructure and (micro-)texture evolution during thermo-mechanical processing requires the combined simulation of plastic deformation and recrystallization. Here, a simulation approach based on the coupling of a full-field dislocation density based crystal plasticity model and a cellular automaton model is presented. A regridding/remeshing procedure is used to transfer data between the deformed mesh of the large-strain crystal plasticity model and the regular grid of the cellular automaton. Moreover, a physics based nucleation criterion has been developed based on dislocation density difference and changes in orientation due to deformation. The developed framework is used to study meta-dynamic recrystallization during double-hit compression tests and multi-stand rolling in high-resolution representative volume elements. These simulations reveal a good agreement with experimental results in terms of texture evolution, mechanical behaviour and growth kinetics, while enabling insights regarding the effect of nucleation on kinetics and crystallographic texture evolution.

Intricate knowledge of dislocation networks in metals has proven paramount in understanding the constitutive behaviour of these materials but current experimental methods yield limited information on the characteristics of these networks. Recently, the isotropic anelastic response of metals has been used to investigate complex dislocation networks through the well-known phenomenon that the observed elastic constants are influenced by dislocations. Considering the dependence of the behaviour of a Frank-Read (FR) source on its initial dislocation character and using discerning characteristics of dislocations, i.e. Burgers vector, line sense and slip system, the present paper takes dislocation character, crystal structure and dislocation network geometry into account and obtains the anisotropic mechanical response for a generic Poisson's ratio. In this work, the tensile test tangent moduli and yield points are presented for spatially uniform and nonuniform dislocation distributions across slip systems. First, the reversible shear strain of the FR source is derived as a function of initial dislocation character. The area swept by a mobile and initially straight dislocation segment pinned at both ends is given as an explicit function of the line stress. Secondly, the anisotropic anelastic strain contribution of FR sources to the total pre- and at-yield strain in single crystallites is calculated. For a given normal stress and superposition of the principal infinitesimal linear elastic lattice strain and anelastic dislocation strain, the tangent moduli are presented. The moduli and the inception of plastic flow have a notable dependence on initial dislocation character, spatial dislocation distribution and loading direction.

The retraction of a cylindrical rigid indenter from a viscoelastic substrate is studied by means of an efficient Green's function method. Hysteresis is observed in the load to area relationship in accordance with experimental results. Although our model relaxes many assumptions posed by LEFM-based analytical theories, the results fall between the limits, at high and low retraction velocities, predicted by the theories. Approaching the high velocity instantaneous limit requires, however, very high velocities or Maugis parameter. The work of adhesion is found to change during retraction. A non-dimensional parameter is proposed to estimate the effect of viscoelasticity in adhesive hysteresis.

It is well established that, at small loads, a linear relation exists between contact area and reduced pressure for elastic bodies with non-adhesive rough surfaces. In the case of adhesive contacts, however, there is not yet a general consensus on whether or not linearity still holds. In this work evidence is provided, through numerical simulations, that the relation is non-linear. The simulations here presented can accurately describe contact between self-affine adhesive rough surfaces, since they rely on Green's function molecular dynamics to describe elastic deformation and on coupled phenomenological traction-separation laws for the interfacial interactions. The analysis is performed for two-dimensional compressible and incompressible bodies under plane strain conditions. Interfaces with various roughness parameters and work of adhesion are considered.

The relative contact area of rough surface contacts is known to increase linearly with reduced pressure, with proportionality factor κ. In its common definition, the reduced pressure contains the root-mean-square gradient (RMSG) of the surface. Although easy to measure, the RMSG of the entire surface does not coincide, at small loads, with the RMSG over the actual contact area g¯ _r, which gives a better description of the contact between rough surfaces. It was recently shown that, for Hertzian contacts, linearity between area and load is indeed obtained only if the RMSG is determined over the actual contact area. Similar to surface contacts, in line contacts, numerical data are often studied using theories that predict linearity by design. In this work, we revisit line contact problems and examine whether or not the assumption of linearity for line contacts holds true. We demonstrate, using Green’s function molecular dynamics simulations, that κ for line contacts is not a constant: It depends on both the reduced pressure and the Hurst exponent. However, linearity holds when the RMSG is measured over the actual contact area. In that case, we could compare κ for line and surface contacts and found that their ratio is approximately 0.9. Finally, by analytically deriving the proportionality factor using g¯ _r in the original model of Greenwood and Williamson, a value is obtained that is surprisingly in good agreement with our numerical results for rough surface contacts.