Continuum Surface Energy from a Lattice Model

We investigate some connections between the continuum and atomistic descriptions of de- formable crystals, using some interesting results from number theory. The energy of a deformed crystal is calculated in the context of a lattice model with binary interactions in two dimensions. A new bond counting approach is used, which reduces the problem to the lattice point problem of number theory. When the crystal shape is a lattice polygon, we show that the energy equals the bulk elastic energy, plus the boundary integral of a surface energy density, plus the sum over the vertices of a corner energy function. This is an exact result when the interatomic potential has finite range; for infinite-range potentials it is asymptotically valid as the lattice parameter zero. The surface energy density is obtained explicitly as a function of the deformation gradient and boundary normal. The corner energy is found as an explicit function of the deformation gradient and the normals of the two facets meeting at the corner. For more general convex domains with possibly curved boundary, the surface energy density depends on the unit normal in a striking way. It is continuous at irrational directions, discontinuous at rational ones and nowhere differ- entiable. This pathology is alarming since it renders the surface energy minimization problem (under domain variations) ill-posed. An alternative approach of defining the continuum region is introduced, that restores continuity of the surface energy density function

Influence of Interface Scattering on Shock Waves in Heterogeneous Solids

In heterogeneous media, the scattering due to interfaces between dissimilar materials play an important role in shock wave dissipation and dispersion. In this work the influence of interface scattering effect on shock waves was studied by impacting flyer plates onto periodically layered polycarbonate/6061 aluminum, polycarbonate/304 stainless steel and polycarbonate/glass composites. The experimental results (using VISAR and stress gauges) indicate that the rise time of the shock front decreases with increasing shock strength, and increases with increasing mechanical impedance mismatch between layers; the strain rate at the shock front increases by about the square of the shock stress. Experimental and numerical results also show that due to interface scattering effect the shock wave velocity in periodically layered composites decreases. In some cases the shock velocity of a layered heterogeneous composite can be lower than that of either of its components

New Solutions for Slow Moving Kinks in a Forced Frenkel-Kontorova Chain

We construct new traveling wave solutions of moving kink type for a modified, driven, dynamic Frenkel-Kontorova model, representing dislocation motion under stress. Formal solutions known so far are inadmissible for velocities below a thresh- old value. The new solutions fill the gap left by this loss of admissibility. Analytical and numerical evidence is presented for their existence; however, dynamic simula- tions suggest that they are probably unstable

On atomistic-to-continuum couplings without ghost forces in three dimensions

In this paper we construct energy based numerical methods free of ghost forces in three dimen- sional lattices arising in crystalline materials. The analysis hinges on establishing a connection of the coupled system to conforming finite elements. Key ingredients are: (i) a new representation of discrete derivatives related to long range interactions of atoms as volume integrals of gradients of piecewise linear functions over bond volumes, and (ii) the construction of an underlying globally continuous function representing the coupled modeling method

A Model for Compression-Weakening Materials and the Elastic Fields due to Contractile Cells

We construct a homogeneous, nonlinear elastic constitutive law, that models aspects of the mechanical behavior of inhomogeneous fibrin networks. Fibers in such networks buckle when in compression. We model this as a loss of stiffness in compression in the stress-strain relations of the homogeneous constitutive model. Problems that model a contracting biological cell in a finite matrix are solved. It is found that matrix displacements and stresses induced by cell contraction decay slower (with distance from the cell) in a compression weakening material, than linear elasticity would predict. This points toward a mechanism for long-range cell mechanosensing. In contrast, an expanding cell would induce displacements that decay faster than in a linear elastic matrix.Comment: 18 pages, 2 figure

Recent Milestones in Unraveling the Full-Field Structure of Dynamic Shear Cracks and Fault Ruptures in Real-Time: From Photoelasticity to Ultrahigh-Speed Digital Image Correlation

The last few decades have seen great achievements in dynamic fracture mechanics. Yet, it was not possible to experimentally quantify the full-field behavior of dynamic fractures, until very recently. Here, we review our recent work on the full-field quantification of the temporal evolution of dynamic shear ruptures. Our newly developed approach based on digital image correlation combined with ultrahigh-speed photography has revolutionized the capabilities of measuring highly transient phenomena and enabled addressing key ques- tions of rupture dynamics. Recent milestones include the visualization of the complete displacement, particle velocity, strain, stress and strain rate fields near growing ruptures, capturing the evolution of dynamic friction during individual rupture growth, and the detailed study of rupture speed limits. For example, dynamic friction has been the big- gest unknown controlling how frictional ruptures develop but it has been impossible, until now, to measure dynamic friction during spontaneous rupture propagation and to understand its dependence on other quantities. Our recent measurements allow, by simul- taneously tracking tractions and sliding speeds on the rupturing interface, to disentangle its complex dependence on the slip, slip velocity, and on their history. In another application, we have uncovered new phenomena that could not be detected with previous methods, such as the formation of pressure shock fronts associated with “supersonic” propagation of shear ruptures in viscoelastic materials where the wave speeds are shown to depend strongly on the strain rate

Spatiotemporal Properties of Sub‐Rayleigh and Supershear Ruptures Inferred From Full‐Field Dynamic Imaging of Laboratory Experiments

Many earthquakes propagate at sub‐Rayleigh speeds. Earthquakes propagating at supershear speeds, though less common, are by far more destructive. Hence, it is important to quantify the motion characteristics associated with both types of earthquake ruptures. Here we report on the spatiotemporal properties of dynamic ruptures measured in our laboratory experiments using the dynamic digital image correlation technique. Earthquakes are mimicked by the frictional rupture propagating along the interface of two Homalite plates. Digital images of the propagating ruptures are captured by an ultrahigh‐speed camera and processed with digital image correlation in order to produce sequences of evolving displacement and velocity maps. Our measurements reveal the full‐field structure of the velocity components, bridge the gap between previous spatially sparse velocimeter measurements available only at two to three locations, and enable us to quantify the attenuation patterns away from the interface

Full-field optical measurement of curvatures in ultra-thin-film–substrate systems in the range of geometrically nonlinear deformations

This article describes coherent gradient sensing (CGS) as an optical, full-field, real-time, nonintrusive, and noncontact technique for the measurement of curvatures and nonuniform curvature changes in film-substrate systems. The technique is applied to the study of curvature fields in thin Al films (6 mum) deposited on thin circular silicon wafers (105 mum) of "large" in-plane dimensions (50.8 mm in diameter) subjected to thermal loading histories. The loading and geometry is such that the system experiences deformations that are clearly within the nonlinear range. The discussion is focused on investigating the limits of the range of the linear relationship between the thermally induced mismatch strain and the substrate curvature, on the degree to which the substrate curvature becomes spatially nonuniform in the range of geometrically nonlinear deformation, and finally, on the bifurcation of deformation mode from axial symmetry to asymmetry with increasing mismatch strain. Results obtained on the basis of both simple models and more-detailed finite-element simulations are compared with the full-field CGS measurements with the purpose of validating the analytical and numerical models

Observations of transient high temperature vortical microstructures in solids during adiabatic shear banding

By using a unique infrared high-speed camera especially constructed for recording highly transient temperature fields at the microscale, we are able to reveal the spatial and temporal microstructure within dynamically growing shear bands in metals. It is found that this structure is highly nonuniform and possesses a transient, short range periodicity in the direction of shear band growth in the form of an array of intense "hot spots" reminiscent of the well-known, shear-induced hydrodynamic instabilities in fluids. This is contrary to the prevailing classical view that describes the deformations and the temperatures within shear bands as being essentially one-dimensional fields. These observations are also reminiscent of the nonuniform structure of localized shear regions believed to exist, at an entirely different length scale, in the earth's lower crust and upper mantle