Sandeep P. Patil

Rate-Dependent Behavior of the Amorphous Phase of Spider Dragline Silk

The time-dependent stress-strain behavior of spider dragline silk was already observed decades ago, and has been attributed to the disordered sequences in silk proteins, which compose the soft amorphous matrix. However, the actual molecular origin and magnitude of internal friction within the amorphous matrix has remained inaccessible, because experimentally decomposing the mechanical response of the amorphous matrix from the embedded crystalline units is challenging. Here, we used atomistic molecular dynamics simulations to obtain friction forces for the relative sliding of peptide chains of Araneus diadematus spider silk within bundles of these chains as a representative unit of the amorphous matrix in silk fibers. We computed the friction coefficient and coefficient of viscosity of the amorphous phase to be in the order of 10(-6) Ns/m and 10(4) Ns/m(2), respectively, by extrapolating our simulation data to the viscous limit. Finally, we used a finite element method for the amorphous phase, solely based on parameters derived from molecular dynamics simulations including the newly determined coefficient of viscosity. With this model the time scales of stress relaxation, creep, and hysteresis were assessed, and found to be in line with the macroscopic time-dependent response of silk fibers. Our results suggest the amorphous phase to be the primary source of viscosity in silk and open up the avenue for finite element method studies of silk fiber mechanics including viscous effects.

Force distribution analysis of mechanochemically reactive dimethylcyclobutene.

Internal molecular forces can guide chemical reactions, yet are not straightforwardly accessible within a quantum mechanical description of the reacting molecules. Here, we present a force-matching force distribution analysis (FM-FDA) to analyze internal forces in molecules. We simulated the ring opening of trans-3,4-dimethylcyclobutene (tDCB) with on-the-fly semiempirical molecular dynamics. The self-consistent density functional tight binding (SCC-DFTB) method accurately described the force-dependent ring-opening kinetics of tDCB, showing quantitative agreement with both experimental and computational data at higher levels. Mechanical force was applied in two different ways, namely, externally by a constant pulling force and internally by embedding tDCB within a strained macrocycle-containing stiff stilbene. We analyzed the distribution of tDCB internal forces in the two different cases by FM-FDA and found that external force gave rise to a symmetric force distribution in the cyclobutene ring, which also scaled linearly with the external force, indicating that the force distribution was uniquely determined by the symmetric architecture of tDCB. In contrast, internal forces due to stiff stilbene resulted in an asymmetric force distribution within tDCB, which indicated a different geometry of force application and supported the important role of linkers in the mechanochemical reactivity of tDCB. In addition, three coordinates were identified through which the distributed forces contributed most to rate acceleration. These coordinates are mostly parallel to the coordinate connecting the two CH3 termini of tDCB. Our results confirm previous observations that the linker outside of the reactive moiety, such as a stretched polymer or a macrocycle, affects its mechanochemical reactivity. We expect FM-FDA to be of wide use to understand and quantitatively predict mechanochemical reactivity, including the challenging cases of systems within strained macrocycles.

Mechanics of Nanostructured Porous Silica Aerogel Resulting from Molecular Dynamics Simulations

Silica aerogels are nanostructured, highly porous solids which have, compared to other soft materials, special mechanical properties, such as extremely low densities. In the present work, the mechanical properties of silica aerogels have been studied with molecular dynamics (MD) simulations. The aerogel model of 192 000 atoms was created with different densities by direct expansion of β-cristobalite and subjected to series of thermal treatments. Because of the high number of atoms and improved modeling procedure, the proposed model was more stable and showed significant improvement in the smoothness of the resulting stress-strain curves in comparison to previous models. Resulting Poissons ratio values for silica aerogels lie between 0.18 and 0.21. The elasticity moduli display a power law dependence on the density, with the exponent estimated to be 3.25 ± 0.1. These results are in excellent agreement with reported experimental as well as computational values. Two different deformation scenarios have been discussed. Under tension, the low-density aerogels were more ductile while the denser ones behaved rather brittle. In the compression simulations of low-density aerogels, deformation occurred without significant increase in stress. However, for high densities, atoms offer a higher resistance to the deformation, resulting in a more stiff response and an early densification. The relationship between different mechanical parameters has been found in the cyclic loading simulations of silica aerogels with different densities. The residual strain grows linearly with the applied strain (≥0.16) and can be approximated by a phenomenological relation ϵp = 1.09ϵmax - 0.12. The dissipation energy also varies with the compressive strain according to a power law with an exponent of 2.31 ± 0.07. Moreover, the tangent modulus under cyclic loading varies exponentially with the compressive strain. The results of the study pave the way toward multiscale modeling of silica as well as reinforced silica aerogels.

Viscous friction between crystalline and amorphous phase of dragline silk.

The hierarchical structure of spider dragline silk is composed of two major constituents, the amorphous phase and crystalline units, and its mechanical response has been attributed to these prime constituents. Silk mechanics, however, might also be influenced by the resistance against sliding of these two phases relative to each other under load. We here used atomistic molecular dynamics (MD) simulations to obtain friction forces for the relative sliding of the amorphous phase and crystalline units of Araneus diadematus spider silk. We computed the coefficient of viscosity of this interface to be in the order of 102 Ns/m2 by extrapolating our simulation data to the viscous limit. Interestingly, this value is two orders of magnitude smaller than the coefficient of viscosity within the amorphous phase. This suggests that sliding along a planar and homogeneous surface of straight polyalanine chains is much less hindered than within entangled disordered chains. Finally, in a simple finite element model, which is based on parameters determined from MD simulations including the newly deduced coefficient of viscosity, we assessed the frictional behavior between these two components for the experimental range of relative pulling velocities. We found that a perfectly relative horizontal motion has no significant resistance against sliding, however, slightly inclined loading causes measurable resistance. Our analysis paves the way towards a finite element model of silk fibers in which crystalline units can slide, move and rearrange themselves in the fiber during loading.

Numerical modelling of the gas detonation process of sheet metal forming

Gas detonation forming is an unconventional technique, which has the potential to form complex geometries, including sharp angles and undercuts in a very short process time. To date, most of the numerical studies on detonation forming neglect the highly dynamic pressure profile of the detonation obtained from experiments. In the present work, it is emphasised that the consideration of the actual detonation pressure as measured in the experiment is crucial. The thickness distribution and radial strain are studied using a strain-rate dependent Johnson-Cook material model. The obtained results vary significantly with change in loading rate. Moreover, the model is capable of predicting extremely sharp edges.

Joining of tubes by gas detonation forming

For many applications, such as in structural components, it is required to join two tubes - sometimes with dissimilar material properties. Only few research studies have investigated the joining of tubular metallic components by means of high-velocity forming processes. In this paper, we present the novel process of joining of two tubes by a gas detonation pressure wave. In particular, the joining of a copper and a steel tube is discussed by means of a finite element study and a conducted experiment.

A NANO-MACRO BOTTOM-UP APPROACH TOWARDS BRITTLE FRACTURE

A novel combined method for highly brittle materials, which provides an efficient and accurate insight into multi-scale fracture modeling, is proposed. In particular, physicallymotivated molecular dynamics simulations are performed to predict crack propagation, in the nanoscale, and therewith determine material and other parameters required for the macroscale modeling under a phase-field continuum approach. The proposed computational approach, which does not require any empirical parameters, contributes towards an improved understanding of mechanics at all length-scale levels.