Seyedeh Mohadeseh Taheri Mousavi

Anisotropy governs strain stiffening in nanotwinned-materials

In their papers, Li et al.1,2 proposed an indentation strain stiffening mechanism to explain the experimentally reported ultra-high hardness of nanotwinned (nt-) cBN3 and nt-diamond4 at extremely small twin boundary (TB) thicknesses (λ ≤ 5 nm). Here, however, we show that the strain stiffening mechanism proposed by these authors is not exclusive to nt-covalent-bonding materials and also exists in nt-metals which are known to exhibit a softening behavior below a critical twin spacing5,6. To demonstrate this, molecular dynamics (MD) simulations were first conducted on a single crystalline Cu sample with dimensions of 40 × 20 × 56 nm3. Periodic boundary conditions were imposed at the boundaries of the sample. The embedded atom method potential for Cu was adopted to describe the interatomic interactions7, and the Nose–Hoover thermostat was used to maintain an NPT ensemble. The time step is set at 1 fs. The sample was relaxed at 1 K for 200 ps. A total shear strain of 90% was applied along the (111)1⁄2112 (hard shearing) and (111) 1⁄2112 (weak shearing) directions at a constant strain rate of 109 s −1. The simulated stress-strain responses displayed in the left inset of Fig. 1 shows that continuous loading along the weak direction causes all atomic layers to be transformed into the hard direction, and the sample exhibits strain stiffening. Next, similar MD simulations were performed on nt-Cu samples, with the same dimensions of 40 × 21 × 50 nm3, containing TBs with thicknesses λ= 0.63 nm, 0.83 nm, 1.25 nm, and 2.54 nm. After 200 ps relaxation, the samples were subjected to both shear and compression loading with ratio fixed at a constant value of tan (68°), which is identical to the loading method applied in the original paper by Li et al.1,2 (see the right inset of Fig. 1). The simulated stress-strain curves in Fig. 1 for nt-Cu exhibit peak strengths which are insensitive to the TB spacing. The stress drops on the curves are caused by the movement of partial dislocations on TBs leading to TB migrations. The saw-tooth pattern continues until the nt-samples become twin-free and the strain stiffening takes place. The slope of the strain stiffening part of the curves is equivalent for all the samples as they share an identical twin-free single crystalline structure. It is noteworthy that the model in Fig. 1, which is identical to that used in refs. 1 and 2, cannot capture the TB-thickness dependence of strength as no grain boundaries have been considered. The above simulations indicate that the strain stiffening mechanism proposed by Li et al.1,2 is not limited to nt-covalent bonding ceramics and exists also in nt-metals. Since it is known that nt-Cu exhibits softening behavior below a critical TB spacing5,6, this finding places significant doubt on the validity of the model of Li et al.1,2 in explaining the observed stiffening behaviors of nt-CBN and nt-Diamond. DOI: 10.1038/s41467-018-03972-9 OPEN

Modeling Abrasion Resistant Materials by Modeling Large Sliding Frictional Contact

— In this paper, our goal is to simulate abrasion resistance material. We therefore need a robust algorithm to model this phenomenon which is a kind of large frictional contact problem. In order to reach to our aim, we have proposed a new method to impose contact constraints in eXtended Finite Element Method (XFEM) framework. In this algorithm, we have modeled large sliding contact problems by using the Node To Segment (NTS) concept. Furthermore, friction between two sliding interface has been modeled based on the Coulomb friction law. In addition, the penalty method which is the most convenient way of imposing non-penetration constraints has been employed. In our algorithm, new Lagrangian shape functions have been used to solve the problems of the conventional Heaviside enrichment function. Finally, a numerical simulation has been delivered to prove the accuracy and capability of our new algorithm.