Eur. Phys. J. E manuscript No. (will be inserted by the editor) Memory effect of external oscillation on residual stress in a paste Jumpei Moritaa,1, Michio Otsukib,2 1Department of Physics and Material Science, Shimane University, Matsue, Shimane 690-8504, Japan. 2Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan. Received: date / Accepted: date Abstract We numerically investigate the stress distribution of a paste when an external oscillation is applied. The paste memorizes the oscillation through plastic deformation. Due to the plastic deformation, the residual stress remains after the oscillation, where the residual stress distribution depends on the number of cycles in the oscillation. As this number increases, the symmetry of the stress distribution is enhanced, which is consistent with the crack patterns observed in the experiments using a drying paste. 1 Introduction criterion. In this regard, recent experiments to mea- sure stresses in drying paste with memory effect [18, When a paste containing powder and water is dried, 19] deserve special attention. According to these ex- cracks are formed [1]. Crack patterns on the surface periments, the fracture criterion is unchanged, but the of the paste are usually random and isotropic in shal- external oscillation induces the residual stress in the di- low containers [2]. However, recent experiments have rection of the oscillation, which increases due to the des- revealed that the crack patterns become anisotropic iccation and leads to the perpendicular crack pattern. on application of external fields [3,4,5,6,7,8,9,10,11, However, the mechanism by which plastic deformation 12,13]. In particular, Nakahara and Matsuo reported causes residual stress is not clear in the experiments. that subjecting the container to horizontal oscillation The stress-measuring experiments [18,19] had been before desiccation imprints a “memory” in the paste. motivated by the theoretical conjecture that anisotropy This memory induces lamellar crack patterns perpen- in residual stresses caused by plastic deformation is re- dicular to the direction of the oscillation after several sponsible for the memory effect [20,21,22]. However, days [3]. theoretical predictions in the previous studies seem to They also discovered that such crack patterns ap- be incoherent. On the one hand, a quasi-linear analy- pear when the stress induced by the oscillation exceeds sis of an elastoplastic model based on the infinitesimal the yield stress of the paste [3]. This experimental re- strain theory was performed in Ref. [20]. In this anal- sult indicates that the plastic deformation caused by the ysis, a non-uniform plastic deformation remains after arXiv:2012.14675v2 [cond-mat.soft] 2 Sep 2021 oscillation is related to the memory effect. The visual- the external oscillation owing to the boundary condi- ization of the plastic deformation in the drying paste tion at the walls of the container. The plastic defor- showing the memory effect supports this conjecture [1]. mation causes an asymmetric stress distribution in the The cracks are formed to release the tension caused direction of the oscillation, and the tension increases in by the desiccation when a fracture criterion is satisfied some parts of the paste. The desiccation after the os- [14,15,16,17]. Therefore, the memory effect can be ex- cillation enhances the asymmetric stress, which creates plained if plastic deformation affects one of the factors cracks perpendicular to the direction of the external in crack formation, namely the tension or the fracture oscillation. Therefore, in this analysis, the crack pat- ae-mail: [email protected] (Current affiliation is Ne- terns are asymmetric. On the other hand, in Ref. [21, turen Co., Ltd.) 22], a nonlinear analysis based on finite strain theory be-mail: [email protected] is performed. This analysis shows that a uniform plas- 2 tic deformation under periodic boundary conditions in- t. Therefore, ∂t represents the derivative with respect creases the tension over the entire area of the paste. to t, with X being constant. This uniform increase of the tension leads to uniform The line element dr between two points labeled as and symmetric crack patterns along the direction of the X and X +dX with an infinitesimal vector dX is given oscillation. by However, the seemingly incoherent predictions may i dr = (∂ir)dX , (3) represent two different limiting cases, as is suggested i by a recent experiment [23], where the symmetric and where ∂i = ∂/∂X , and we use the Einstein notation for asymmetric crack patterns occur depending on the num- the summation of repeated indices. The square of the ber of cycles in the external oscillation. For a few cycles, distance between the points in the current configuration the asymmetric crack pattern predicted in the quasi- is linear analysis was formed. For a sufficiently large num- 2 i 2 i j (ds) = (∂ir)dX = gij dX dX , (4) ber of cycles, the symmetric pattern predicted in the | | nonlinear analysis was formed. However, the present where gij is a component of the Euclidean metric tensor theories cannot explain the dependence of the crack g given by pattern on the number of cycles. Therefore, a unified gij = (∂ir) (∂j r). (5) theoretical model is required. · In this study, we formulate such a unified model, Assuming the existence of a local stress-free “natu- with which we numerically examine the residual stress ral state”, the line element between the two points in in a paste after external oscillation. In Sec. 2, we re- this state is denoted as dr♮ [21]. The square of the dis- view a Lagrangian description of kinematics. In Sec. tance between the points in the “natural state” is given 3, we derive the time evolution equations for a paste. by The memory effect of the external oscillation is numeri- cally investigated in Sec. 4. The setup of our simulation (ds♮)2 = dr♮ 2 = g♮ dXidXj , (6) | | ij is shown in Sec. 4.1. We demonstrate the relation be- ♮ ♮ tween the plastic deformation and the residual stress where gij is the natural metric tensor satisfying gij = ♮ in Sec. 4.2. In Secs. 4.3 and 4.4, we present the resid- gji. The change in the natural metric tensor represents ual stress after the oscillation for different numbers of the plastic deformation. cycles, respectively. The dependence of the stress dis- The elastic deformation is represented by the trans- tribution on the number of cycles is examined in Sec. formation from dr♮ to dr: 4.5. We discuss our results and conclude our paper in r F r♮ Sec. 5. In Appendix A, the left Cauchy–Green tensor in d = d , (7) our setup is derived. The stretching tensor is obtained where F is the deformation gradient tensor. The left in Appendix B. The constitutive relation of the paste is Cauchy–Green Tensor is defined as [24,25] derived in Appendix C. In Appendix D, we discuss the T shear stress under steady shear in our model. B = FF . (8) B is given by 2 Kinematics ij B = g (∂ir) (∂j r), (9) ♮ ⊗ The configuration of a 3-dimensional paste is repre- ij where (g♮ ) denotes the inverse of the component matrix sented by the current coordinates r mapped from the ♮ of the metric tensor (g ) satisfying reference coordinates X as ij ♮ jk k x(X,Y,Z,t) gij g♮ = δi . (10) (X,t) r(X,t)= y(X,Y,Z,t) , (1) 7→ See Appendix A for the derivation of Eq. (9). z(X,Y,Z,t) C The difference in velocities between two points la- where [ ]C denotes the representation in terms of Carte- beled as X and X + dX is given by · sian components. The velocity v is given by i dv = (∂iv) dX . (11) v(X,t)= ∂tr(X,t), (2) The time rate of the deformation is described by the 2r and the acceleration is given by ∂t . Here, ∂t stands velocity gradient tensor W, which satisfies [24,25] for the time derivative in the Lagrangian description, and the independent variables of the paste are X and dv = Wdr. (12) 3 The stretching tensor is defined as the Hadamard strain energy. Substituting Eq. (19) with Eqs. (20) and (21) into Eq. (18), we obtain 1 T D = W + W . (13) (EP) 2 σ = Σ + I3Ψ ′(I3)+ µ I + µ (B I) . (22) { } − As shown in Appendix B, the stretching tensor is rep- Here, I is the unit tensor. See Appendix C for the deriva- resented as tion of Eq. (22). For σ(V), we adopt the linear viscous 1 ij stress tensor as D = ∂tg (∂ir) (∂j r), (14) −2 ⊗ σ(V) =2ηD (23) where (gij ) is the inverse of the component matrix of with the viscosity η [24]. the metric tensor (gij ) satisfying Following Ref. [21,22], the plastic deformation is de- jk k gij g = δi . (15) scribed by the temporal evolution of the natural metric ij g♮ as ij ij 3 Equation of motion and constitutive relation (1 + τ∂t)g♮ = Kg . (24) The equation of motion for r(X,t) is given by Here, K is given by 3 3 ρ∂2r = divσ + f , (16) K = = (25) t e ♮ ij I gij g 1 where ρ, σ, and fe are the density in the current config- ♮ g♮ uration, the Cauchy stress tensor, and the body force, to satisfy the incompressibility condition J = det = 1 in the natural state. It should be noted that we assume respectively. The Cauchy stress tensor consists of an p elastoplastic part σ(EP) and a viscous part σ(V) as effective incompressibility for g by adopting Poisson’s ratio ν λ/(λ + µ)/2 1/2 in our numerical simula- ≡ ≃ σ = σ(EP) + σ(V).