# Compaction of quasi-one-dimensional elastoplastic materials

M. Reza Shaebani, Javad Najafi, Ali Farnudi, Daniel Bonn, Mehdi Habibi

► To cite this version:

M. Reza Shaebani, Javad Najafi, Ali Farnudi, Daniel Bonn, Mehdi Habibi. Compaction of quasi-one-dimensional elastoplastic materials. Nature Communications, 2017, 8 (1), pp.15568. 10.1038/ncomms15568. hal-04304636

HAL Id: hal-04304636

https://hal.science/hal-04304636v1

Submitted on 24 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
ARTICLE

Received 7 Jul 2016 | Accepted 7 Apr 2017 | Published 6 Jun 2017

DOI: 10.1038/ncomms15568

OPEN

# Compaction of quasi-one-dimensional elastoplastic materials

M. Reza Shaebani<fnref n="1" />, Javad Najafi²,³, Ali Farnudi³, Daniel Bonn⁴ & Mehdi Habibi³,⁴,†

Insight into crumpling or compaction of one-dimensional objects is important for understanding biopolymer packaging and designing innovative technological devices. By compacting various types of wires in rigid confinements and characterizing the morphology of the resulting crumpled structures, here, we report how friction, plasticity and torsion enhance disorder, leading to a transition from coiled to folded morphologies. In the latter case, where folding dominates the crumpling process, we find that reducing the relative wire thickness counter-intuitively causes the maximum packing density to decrease. The segment size distribution gradually becomes more asymmetric during compaction, reflecting an increase of spatial correlations. We introduce a self-avoiding random walk model and verify that the cumulative injected wire length follows a universal dependence on segment size, allowing for the prediction of the efficiency of compaction as a function of material properties, container size and injection force.

<footnote n="1">¹Department of Theoretical Physics, Saarland University, 66041 Saarbrücken, Germany. ²Department of Experimental Physics, Saarland University, 66041 Saarbrücken, Germany. ³Department of Physics, Institute for Advanced Studies in Basic Sciences, Zanjan 45195, Iran. ⁴Van der Waals-Zeeman Institute, University of Amsterdam, 1098 XH Amsterdam, The Netherlands. †Present address: Laboratory of Physics and Physical Chemistry of Foods, Wageningen University, Wageningen, The Netherlands. Correspondence and requests for materials should be addressed to M.R.S. (email: shaebani@lusi.uni-sb.de) or to D.B. (email: D.Bonn@uva.nl) or to M.H. (email: M.Habibi@uva.nl).</footnote>
Compaction of slender objects in confined geometries is ubiquitous in nature. Perhaps the most important example is DNA packaging in viral and bacteriophage capsids and cell nuclei¹⁻⁵. Other pertinent examples are the folding of insect wings in cocoons⁶, and flower or plant leaves in buds⁷. The process of compaction may result in complex morphologies depending on the applied forces and constraints⁸. Recent numerical studies⁹⁻¹³ showed that the space-filling properties of two-dimensional (2D) crumpled sheets are influenced by parameters such as self-avoidance and plasticity-ingredients that are difficult to disentangle in experiments. Self-avoidance alters the hierarchical nature of the compaction process and induces stronger self-correlations as the compression increases¹³,¹⁴. Thus, considering the structural evolution is key for understanding the efficiency of *in vitro* compaction.

Aiming to provide quantitative insights into the role of self-avoidance, we turn to one-dimensional (1D) wires. Due to its very nonlinear nature it is easier to study 1D systems than the more complicated crumpling process in 2D sheets¹³,¹⁵⁻¹⁷. For 1D-compaction, how the morphology of crumpled objects develops is of particular importance in technological and biological applications as, for example, in endovascular coiling treatment of cerebral aneurysms¹⁸ or in packing of DNA¹⁹.

When compacting elastic low-frictional wires with a high bending rigidity in confined geometries in such a way that the internal torsion is released, highly ordered structures with distinctly oriented subdomains of parallel coils form (Fig. 1). However, with increasing friction²⁰ or plasticity⁹, or by accumulating torsion during the packing process²¹,²², disordered structures emerge where the contribution of folds or bends in the morphology is more pronounced. For example, by introducing the number of segments as the order parameter, it has been recently shown that a sharp transition from ordered (coiled) to disordered (folded) structures occurs as the friction increases²⁰. In disordered morphologies, the compaction efficiency is controlled to a large extent by spatial exclusion effects, which continuously evolve in the course of compaction. Hence, unravelling the mechanisms that govern the evolution of self-avoidance is crucial to achieve an efficient compaction.

Here, we study the morphologies of wires packed into rigid spherical containers and find that the maximum packing density in disordered structures decreases with reduced thickness of the wire (or, equivalently, increasing the container size). To elaborate on the underlying mechanisms leading to this peculiar behaviour, we isolate the influence of self-avoidance by focusing on the compaction of plastic frictional wires, where folding is dominant in the resulting structure. By following the morphological evolution, a gradual crossover from random to correlated folding events is observed due to spatial exclusion effects. We propose that the compaction can be considered as a confined self-avoiding

random walk (SAW). In such far-from-equilibrium processes, the imposed constraints and initial conditions do not uniquely determine the final crumpled state. Instead, there is an ensemble of admissible configurations, from which some structural properties of the system can be derived. We introduce a SAW sampling method that successfully accounts for the time evolution of the wire segment length. We thus present a more complete understanding of the compaction: the maximum length of the injected wire can be estimated from the geometry and imposed constraints for a given set of material parameters.

## Results

**Universal phase diagram for 1D crumpling.** We first consider the packing of elastic low-frictional wires with a high bending rigidity in rigid spherical containers. When the wire is allowed to axially rotate at the injection point to release the torsion during the packing process, highly ordered coils form as the wire relaxes towards a global minimum energy. By hindering the release of torsion, the wire buckles more frequently to free elastic energy. Hence, the packing process becomes less ordered²¹, leading to warped structures similar to those obtained numerically for compaction of DNA molecules in phage capsids⁴ (see Figs 1 and 2a). The disorder is also enhanced by friction, which causes the wire to resist against sliding and to randomly bend due to local constraints²⁰. Another property which obviously affects the morphology is the degree of plasticity of the wire. While the bending rigidity of the plastic wires can be quite high, their yield stress is relatively low, leading to structures with rather straight segments and sharp turnings. Upon increasing plasticity (that is, lowering the yield stress), the irreversible deformations of wire increase the disorder of the crumpled configuration. One can map out a qualitative phase diagram for the morphological evolution of the resulting crumpled structures in the space of wire properties (friction, torsion and plasticity), as depicted in Fig. 2b. More generally, disordered structures can be generated in diverse ways by tuning the wire or container properties. The morphological phase space indeed contains additional degrees of freedom associated with container properties, such as its flexibility²⁰, shape²³ or the degree of confinement imposed by it (characterized by the container size $R$ relative to the radius of gyration $R_g$ of the crumpled structure and also to the persistence length $\lambda$ of the elastic wire). For example, a biopolymer coils itself inside the cage if $\lambda$ is comparable to $R$ (for example, in packing of DNA in icosahedral bacteriophages), while for weak confinement, that is, $\lambda \ll R_g < R$, the biopolymer chain (such as chromatin) has a relatively low bending stiffness and behaves as a SAW without 'feeling' the boundaries. For $\lambda$ values in between, the morphological evolution during the crumpling process is

<figure><img src="image_2.png" /><figcaption>**Figure 1** | Experimental set-up and examples of distinct morphologies. (a) Cross-sectional view of the experimental set-up used to compact 1D wires in a spherical rigid container. (b) Distinct morphologies obtained by compacting wires with different material properties introduced in the Methods section. Scale bar, 1 cm.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2** | Morphological phase diagram and scaling of geometrical characteristics. (a) Examples of morphological differences upon increasing key material parameters. (b) Schematic morphological phase diagram in the space spanned by torsion, friction and plasticity. (c) The maximum packing density $\phi_{max}$ versus the relative system size $R/r$ (parameters defined in Fig. 1a). (d,e) Scaling of the number of bends $N$ (d), and the dimensionless mean segment size $\langle l \rangle/R$ (e) versus the effective system size $R/r$, for plastic frictional wires. The filled circles (open triangles) indicate experimental (simulation) results. The dashed lines are power-law fits to the experimental data. The crosses denote the simulation results obtained when assuming that the packing fraction is the influential parameter on the self-avoidance effects. The open circles in (d) show the experimental results obtained from the high frictional set-up. Error bars correspond to s.d. of five separate measurements.</figcaption></figure>

complicated due to varying combined effects of self-avoidance and interactions with boundaries$^{19,24}$. It has been also shown that a transition from coiled to disordered configurations occur, as the accessible space reduces during the compaction of elastic rods$^{25}$.

To compare the compaction efficiencies, we measure the packing density once the injection of wire eventually stops. Indeed, the value obtained for this maximum packing density, $\phi_{max}$, depends on the wire radius $r$, the container size $R$ and the insertion force. We measure this quantity for a given insertion force and for different combinations of inserted wire radius $r$ and container radius $R$. When plotting $\phi_{max}$ versus the non-dimensional system size $R/r$ for the coiled compact morphologies of low-plasticity, low-friction, low-torsion wires (Fig. 2c), a plateau at small $R/r$ followed by a weak decrease at larger values of $R/r$ is observed. It was shown with geometrical arguments$^{22}$ that $\phi_{max}$ slightly decays with $R/r$ for a purely coiled structure in a spherical container. Note that the very inner core of the structure practically becomes disordered as the accessible space reduces and its shape becomes more irregular (which makes the formation of coils more difficult). This disordered core (with a possible $R$-dependent size) can also contribute to the weak decay of $\phi_{max}$ versus $R/r$. It has been shown that the packing fraction decreases with increasing disorder in packings of elastic wire$^{26}$.

When increasing plasticity, friction and/or torsion, resulting in the formation of folds and bends, the data collapse onto curves following a power-law $\phi_{max} \sim (R/r)^{-\alpha} \sim D^{-3}$, with $D$ being the fractal dimension$^{27,28}$. The slope of the curve depends on the

degree of disorder. For example, lubricating the inner wall of the container with silicon oil leads to the formation of highly ordered coils at the outer layer of crumpled plastic wires, which results in a mixed coiled-folded structure with $\alpha \simeq 0.38 \pm 0.03$. A similar exponent is obtained for the compaction of elastic torsional wires where coils and bends coexist. The steepest descent is observed for crumpling of plastic wires at high friction, where coils are absent and folding is the dominant process ($\alpha \simeq 0.52 \pm 0.05$). While the very weak system-size dependence of efficient compaction in ordered (coiled) structures is understandable, the behaviour of disordered morphologies is counter-intuitive, as one would expect that relatively thinner wires more flexibly fill a given container, leading to a higher compaction efficiency. A similar trend for the dependence of packing density on the relative system size was reported in experiments on DNA packaging in viral capsids$^{29}$, revealing that in spite of the huge differences in length scales of the two systems, the maximum packing densities behave similarly in the presence of disorder. Self-avoidance inside a confinement can explain the peculiar behaviour of $\phi_{max}$ versus $R/r$ via a mean-field interpretation, assuming that the self-avoidance energy originates mainly from the homogeneously distributed binary contacts between the wires (whose density nearly grows as the square of the packing fraction), and also supposing that the local radius of curvature of the confinement is comparable to the container size and varies slowly. When balancing the confining energy and self-avoidance$^{30–32}$, the lowest (harmonic) approximation of the confining energy yields
an energy density of the order of $(R/r)^{-1}$ while the self-avoidance energy density is proportional to $\phi^2 \sim (R/r)^{-2x}$. By equalizing these energy densities we obtain $\phi \sim (R/r)^{-0.5}$.

**Segment size statistics.** To better understand the influence of disorder, we choose plastic frictional wires to avoid ordered coils and create the highest possible disorder in the crumpled structure. After compacting the wires, we open the moulds and investigate the resulting compacted structures by analyzing the folding statistics. The points of folding were often determined by sharp changes of wire orientation. If they were not easily distinguishable, then a minimum threshold of 90° for the turning angle of the wire, and a maximum threshold of $R/2$ for the radius of curvature were imposed. We cut the wire at each of the folding points, straighten the segments and measure their length. Straightening of the curved segments rarely allows for segment lengths $\ell$ longer than the container diameter, but we checked that the maximum segment length $\ell_{max}$ remains smaller than $\pi R$ in the absence of coils. We preserve the order of the wire segments and average the results over five realizations for each value of $R/r$ to obtain the sequence of the segment lengths $\ell_n$.

A key observation is the scaling of the total number of segments $N$ with the effective system size $R$. Similar scaling laws were reported for 2D packings of wires$^{27,30,33}$. As shown in Fig. 2d, a power-law relation of the form

$$N \sim (R/r)^\beta \quad (1)$$

holds with $\beta = 1.86 \pm 0.03$ for smooth wire and container with wire-wire and container-wire friction coefficients $\mu_{ww} \simeq 0.2$ and $\mu_{cw} \simeq 0.4$, respectively. The exponent can be understood by considering the wire crumpling process as a SAW in confinement. For comparison, for the number of steps on a cubic lattice in an ordinary random walk $\beta$ equals 2, while for SAW $\beta \simeq 5/3$ (ref. 34). The fact that we find an exponent in between these two values can be understood because there is a gradual evolution of spatial correlations over the course of crumpling (see below), thus, the exponent continuously decreases from 2. However, the wires can slide over each other due to the finite friction so that the self-avoidance constraint is only partially fulfilled. For comparison, we repeated the experiment by roughening the plastic wires and the inner surface of the moulds to increase the friction coefficients to $\mu_{ww} \simeq 0.45$ and $\mu_{cw} \simeq 0.45$. The considerable change in the wire-wire friction resulted in a smaller exponent $\beta = 1.75 \pm 0.05$ which is closer to the pure self-avoidance limit (see Fig. 2d). From the scaling of $\phi_{max}$ and $N$ with $R/r$ one expects that the normalized mean segment size $\langle \ell \rangle / R$ follows $(R/r)^{-0.4}$, as confirmed by the experimental results in Fig. 2e.

**Evolution of spatial correlations.** At earlier stages of the crumpling process in a given mould, the injected wire proceeds in the container without interacting with the accumulated wire. Assuming that the plastic wire bends at a random point between the injecting hole and a contact point at the container surface, the resulting segment length $\ell$ is a random variable, symmetrically distributed between 0 and the maximum possible segment length $\ell_{max}$. By increasing the total length $L$ of accumulated wire, spatial exclusion effects grow and the injected wire cannot easily proceed through the sphere without touching the crumpled structure. Hence, long segments gradually become less probable and the probability distribution $P(\ell/R)$ of the normalized segment size becomes more asymmetric due to relatively large populations of smaller segments (see Fig. 3a). When comparing the final structures (that is, those obtained when the injection of wire stops), we interestingly find that for larger values of $R/r$ the segment size distribution $P(\ell/R)$ is more asymmetric and shifts towards smaller

<figure><img src="image_4.png" /><figcaption>Figure 3 | Probability distribution of the normalized segment size.</figcaption></figure>

The symbols (lines) represent the experimental (simulation) results for $r = 0.5$ mm. (a) results obtained by injecting wires of different total length $L$ into a container with radius $R = 14.5$ mm. (b) results of injecting the maximum possible length of wire (when applying the insertion force of nearly $F = 100$ N) in different container sizes $R$.

segment sizes (Fig. 3b). This behaviour similarly indicates the growth of spatial exclusion effects with increasing $R/r$. Note that the initial segment sizes $\ell$ are only determined by the container size $R$ in all containers, however, $\ell$ gradually decreases as the spatial exclusion effects grow. The effect is more pronounced for larger spheres as the crumpling process continues further.

**Self-avoiding random walk model.** We argue that the strength of self-avoidance effects is indeed captured by the total length $L$ of the injected wire, rather than the total volume excluded by it (that is, the packing fraction). The exclusion effect that an inserted rod-like object experiences inside a crumpled structure is effectively determined by the projection of the crumpled wire on a plane perpendicular to the direction of insertion. Therefore, the total length and the thickness of the crumpled wire are expected to be the influential parameters. However, the circular cross-section of wires reduces the contact area between the touching wires and, thus, the effective frictional force between them. As a result, the self-avoidance effects are not proportionally increasing with the wire thickness (that is, $r$). We conclude that the entire contribution to the spatial exclusion constraint can be attributed to the length of wire, reflected in the dimensionless quantity $\lambda = L/R$ which grows as $\lambda \sim (R/r)^{1.5}$ (while $\phi$ decreases as $\phi \sim (R/r)^{-0.5}$). In the following, we simulate the folding process as a SAW of the wire inside the confinement. While the existing SAW algorithms mainly follow stochastic Markovian dynamics to sample the ensemble of trajectories on regular lattices, here we propose an alternative approach which accounts for the time evolution of the step size $\ell$. We suppose that the strength of self-avoidance effects after the $n$-th segment is mainly controlled by the length $L_n$ of the inserted wire, that is, the larger is the parameter $\lambda_n = L_n/R$, the smaller is the success probability for the segment $n+1$ to be a long one. The size $\ell_{n+1}$ of the next segment is obtained via the following algorithm: a trial segment size $\ell$ is chosen randomly within $[0, \ell_{max}]$, with $\ell_{max}$ being the maximum segment length obtained in experiments for a given value of $R$. The proposed $\ell$ is accepted according to a Metropolis-like criterion with probability

$$P_{n+1}(\ell) = \mathcal{N}^{-1} \exp[-\kappa \lambda_n \ell / R], \quad (2)$$

where $\mathcal{N} = R/(\kappa \lambda_n)$ ($1 - \exp[-\kappa \lambda_n \ell_{max}/R]$) is the normalization factor. The coefficient $\kappa$ depends on wire properties and is treated as a free parameter to take into account the partial fulfillment of the self-avoidance constraint due to sliding of the wires. While the exponents are not affected by the choice of $\kappa$, by fitting it we can quantitatively reproduce the experimental data. For the sake of simplicity, here we used a single averaged value of $\kappa$ to reproduce all
<figure><img src="image_5.png" /><figcaption>**Figure 4** | Evolution of the tail of the segment-size distribution. The results obtained from the simulations for $r = 0.5$ mm. (a-c) Comparison between the tail of $P(\ell)$ for different values of the total length $L$ of the injected wire in a container with radius $R = 14.5$ mm. (d-f) A similar comparison when the maximum possible length of wire is injected in different container sizes $R$. The solid lines denote the simulation results, and the dotted, dashed, and dashed-dotted guidelines represent, respectively, the log-normal $P(x) = \frac{1}{\sqrt{2\pi\sigma^2}} \exp[-\frac{(\ln x - \lambda)^2}{2\sigma^2}]$, gamma $\lambda \sim \frac{\kappa^{\kappa-1}}{\Gamma(\kappa)} \exp[-x/\lambda]$ and Gaussian $\frac{1}{\sqrt{2\pi\sigma^2}} \exp[-\frac{(x-\lambda)^2}{2\sigma^2}]$ distributions, plotted with the same mean and variance as the corresponding experimental data. In the gamma distribution, $\Gamma(\sigma)$ denotes the gamma function evaluated at $\sigma$.</figcaption></figure>

the experimental data. In the case of rejection, a new $\ell$ is tried. Finally, the cumulative length is updated as $L_{n+1}=L_n + \ell$ before starting the next step. Equation (2) assumes an exponentially lower acceptance chance for larger trial lengths $\ell$. Moreover, the acceptance probability decreases with increasing $\lambda_n$, as the self-avoidance effects become more pronounced. The method samples the segment-length landscape according to a Boltzmann-like distribution. Initially, the wire walks in free space ($\lambda_n = 0$), thus, $P_{n+1}(\ell) = \frac{1}{\ell_{\max}}$ independently of the trial length $\ell$. However, larger values of $\ell$ become gradually less probable with increasing $\lambda_n$. We perform extensive Monte Carlo simulations by adjusting $\ell_{\max}$ and the threshold value of $\phi$ to the experimental data. The results shown in Figs 2 and 3 are in remarkable agreement with experiments. We checked that the power-law scalings cannot be reproduced when replacing $\lambda_n$ with $\phi_n$ in equation (2) (see Fig. 2d,e). Notably, the numerical predictions for the total wire length $L$ differ less than 4% from the experimental values for all system sizes.

In Fig. 4 we take a closer look at the tail of $P(\ell/R)$ obtained from the numerical simulations. While Gaussian function represents the distribution of random uncorrelated data, gamma and log-normal functions are respectively associated with random events in the presence of self-correlations and those that occur hierarchically$^{8,11,13,35-37}$. Since the noisy data of the experimental tail prevents any conclusive statement on the tail behaviour of $P(\ell/R)$, we plot these three functions (with the same mean and variance as the experimental data) and use them as guidelines to demonstrate the trend of the tail behaviour obtained from the numerical simulations. We investigate the evolution of the tail during the injection of wire in a given container, and also

<figure><img src="image_6.png" /><figcaption>**Figure 5** | Cumulative length $L_n$ of inserted wire after the n-th segment. (a) experimental results, (b) simulation results. Inset shows $L_n/R$ averaged over all experiments (red solid line) and its s.d. (grey shaded area). The dotted line is obtained from equation (4).</figcaption></figure>

compare the tails when the maximum possible length of wire is injected in different container sizes. As shown in Fig. 4, the tail is better captured by the Gaussian for small $R/r$ or at the early stages of crumpling, while there is a gradual crossover towards the Gamma distribution, either by increasing the container size or by increasing the length of the injected wire. Thus, self-correlations are the dominant underlying mechanism here. A hierarchical folding mechanism is expected to cause a rather stable log-normal distribution tail over all timescales, thus, the evolution of the tail is in favour of evolving correlated events. It has been previously shown numerically$^{13}$ and by compacting of rods in 2D experiments$^{14}$, that self-avoidance alters the hierarchical nature of crumpling at high compression and induces self-correlations.

The cumulative length $L_n$ of the inserted wire after the n-th segment qualitatively collapses onto a master curve for different values of $R/r$ (Fig. 5). The segments are initially independent of each other and $L_n/R$ grows linearly with $n$. The steps however become more correlated with increasing $n$, leading to a slower growth of $L_n/R$. A similar reduction of the slope has been recently observed for motor-driven viral packaging$^{38}$. From the scaling of $\phi$ and $N$ with $R/r$, one obtains the asymptotic scaling $L/R \sim N^{0.8}$, which is consistent with experiments [see Fig. 5(inset)]. Starting from $\langle \ell_1 \rangle = \ell_{\max}/2$, the mean segment size at next steps can be estimated in terms of $L_n$ as

$$ \langle \ell_{n+1} \rangle = \int_0^{\ell_{\max}} P_{n+1}(\ell) \ell d\ell = \frac{R}{\kappa \lambda_n} + \frac{\ell_{\max}}{1 - e^{\kappa \lambda_n \ell_{\max}/R}} \quad (3) $$

Hence, we obtain the following recursive analytical expression for the injected length of wire after the n-th segment

$$ \frac{L_n}{R} = \frac{\ell_{\max}}{2R} + \sum_{i=1}^{n-1} \left[ \frac{R}{\kappa L_i} + \frac{\ell_{\max}/R}{1 - \exp[\kappa \ell_{\max} L_i / R^2]} \right] \quad (4) $$

in excellent agreement with the data as shown in Fig. 5(inset).

**Buckling threshold.** The inset of Fig. 6 shows that the injection of wire eventually stops at a cutoff segment length $\ell_c$ which is
<figure><img src="image_7.png" /></figure>

**Figure 6** | Universal filling mechanism. Collapse of $L_n/R$ versus the scaled segment length $\ell_n/R$ is shown for different values of $R$ in experiments (background grey lines) and simulations (thick coloured lines). The symbols show the experimental cutoff values $\ell_c/R$ for $r=0.5$ mm. The red dashed lines indicate the experimental results when injecting wires of different total length $L \approx 100, 170$ or $250$ cm into a container with radius $R=14.5$ mm. Inset shows cutoff segment length $\ell_c$ versus the container radius $R$ for $r=0.3$ mm (open circles) and $r=0.5$ mm (filled circles). The horizontal lines show the average values $\bar{\ell}_c=2.3$ and $5.8$ mm. Error bars correspond to s.d. of five separate experimental measurements.

independent of $R$ for a given $r$. It follows that the process stops when the segment size becomes so small that the applied feeding force $F$ exceeds the maximum power of the motor. The minimum threshold length can be obtained from the Euler buckling theory as $\ell_c = \sqrt{\frac{\pi^3 Y}{16F}} r^2$ so that $\ell_c$ is determined by the insertion force and wire properties (Young's modulus $Y$) and is independent of $R$. For $r=0.3$ mm (0.5 mm) and a maximum insertion force of nearly $F=100$ N in our set-up, we obtain $\ell_c \approx 2.2$ mm (6.0 mm), close to the experimental values of 2.3 and 5.8 mm [Fig. 6(inset)].

Figure 6 interestingly evidences a universal filling mechanism independent of $R$. The normalized cumulative wire length $L_n/R$ vs. the normalized length of the $n$-th segment $\ell_n/R$ collapses onto a universal curve for different values of $R$. As the wire injection continues, $\ell_n/R$ gradually decreases until it reaches the minimum threshold value $\ell_c/R$, where the process eventually stops. For a given value of $r$, $\ell_c$ is independent of $R$, thus, $\ell_c/R$ decreases (that is, shifts to the left in Fig. 6) with increasing $R$. Consequently, the cutoff number of the steps at which the process stops increases with $R$. By calculating $\ell_c$ from the insertion force and wire properties, one can predict the total length of crumpled wire for different $R$. The data collapse is also obtained when injecting wires of different total length $L$ into a given container (dashed lines in Fig. 6). Indeed, there is no significant difference between the sequence of the segment sizes $\ell_n$ for different total lengths of the injected wire. This shows that the previously formed segments are not considerably affected during the compaction process, evidencing that the hierarchical folding events rarely happen.

In summary, we reported the compaction of 1D objects in spherical containers and showed how the morphology evolves from ordered (coiling) to disordered (folding and bending) structures in the phase space spanned by friction, torsion and plasticity. The disorder reduces the compaction efficiency and causes a nontrivial system-size dependence, which is explained by SAWs in confined geometries. Monitoring the evolution of segment-length

distributions in highly disordered structures of plastic frictional wires showed that the compaction process is correlated: the longer the injected wire, the stronger the spatial exclusion effects leading to shorter segments. The self-avoidance constraint is only partially fulfilled due to sliding of the wires, leading to an exponent $\beta$ slightly larger than 5/3 as reported for SAWs. Our results provide new insight into underlying mechanisms of crumpling beyond the simple hierarchical description of the process, which also helps to better understand the reverse processes, for example, viral DNA ejection³⁹ and unpacking of crumpled wires⁴⁰. While more detailed morphological information can be obtained via Discrete-Element Method (DEM) simulations, these are however computationally expensive. Our proposed sampling method opens the door to relatively simple Monte Carlo simulations of SAWs inside arbitrary confinements to obtain some of the macroscopic quantities of interest for example, the arbitrary moments of segment size distribution.

**Methods**

**Experimental set-up.** The experimental set-up consists of a rigid hollow spherical container of inner radius $R$ with a small hole to insert the wire (see Fig. 1). Several transparent polymeric moulds with radii $R \in [4,30]$ mm were used. A small nozzle and two counterrotating rollers were attached to the injecting hole to facilitate the control of the insertion speed.

**Material properties.** As a model elastoplastic material, we chose solder wire $Sn_{60}Pb_{40}$ with Young's modulus $Y \approx 30$ GPa and yield stress $\sigma \approx 28$ MPa. For the elastic wire experiments, we mainly used fishing line with Young's modulus $Y = 2.00 \pm 0.01$ GPa (obtained experimentally by tensile tests). Moreover, elastic silicon wires and cotton threads with relatively lower Young's moduli $Y \approx 5.0$ and 0.8 MPa were also used. The wire-wire and container-wire friction coefficients, using smooth wires and container walls, were $\mu_{ww} = 0.20 \pm 0.02$ and $\mu_{cw} = 0.40 \pm 0.02$, respectively. By roughening the plastic wires and the inner surface of the polymeric moulds with sandpaper, we obtained higher friction coefficients $\mu_{ww} = 0.45 \pm 0.02$ and $\mu_{cw} = 0.45 \pm 0.02$. We also used smooth lubricated plastic wires and inner surfaces of the moulds to lower the friction coefficients, leading to $\mu_{ww} = 0.12 \pm 0.02$ and $\mu_{cw} = 0.18 \pm 0.02$. The lubrication was done with silicon oil.

**Insertion process and imaging.** We inserted wires of radius $r=0.4, 0.5, 0.6, 0.8$ or $1.4$ mm into the moulds with a slow feeding speed of about $1 \text{ mm s}^{-1}$ to avoid inertial effects. We checked that the results are independent of the feeding speed in the quasi-static compaction regime. The insertion process continued with the constant speed until the insertion force exceeded the power threshold of the motor and the wire buckled outside of the container. The final plastic-wire structure preserves its shape after opening the mould allowing for a detailed analysis of morphological changes, which can be considered as plastic deformations. In the low-torsion elastic set-up, we allowed axial rotation of the wire between the nozzle and the sphere. The images presented in Figs 1 and 2 were taken by a camera with pixel resolution of 70 µm. Before opening the moulds, we filled them with a transparent gel in the case of elastic wires to preserve the shape of the final structure.

**Data availability.** The data that support the findings of this work are available from the corresponding authors on request.

**References**

1. Purohit, P. K., Kondev, J. & Phillips, R. Mechanics of DNA packaging in viruses. Proc. Natl Acad. Sci. USA **100**, 3173–3178 (2003).

2. Ali, I., Marenduzzo, D. & Yeomans, J. M. Polymer packaging and ejection in viral capsids: shape matters. Phys. Rev. Lett. **96**, 208102 (2006).

3. Katzav, E., Adda-Bedia, M. & Boudaoud, A. A statistical approach to close packing of elastic rods and to DNA packaging in viral capsids. Proc. Natl Acad. Sci. USA **103**, 18900–18904 (2006).

4. Kindt, J., Tzilil, S., Ben-Shaul, A. & Gelbart, W. M. DNA packaging and ejection forces in bacteriophage. Proc. Natl Acad. Sci. USA **98**, 13671–13674 (2001).

5. Smith, D. E. et al. The bacteriophage φ29 portal motor can package DNA against a large internal force. Nature **413**, 748–752 (2001).

6. Brackenbury, J. H. Wing folding and free-flight kinematics in Coleoptera (Insecta): a comparative study. J. Zool. **232**, 253–283 (1994).

7. Kobayashi, H., Kresling, B. & Vincent, J. F. V. The geometry of unfolding tree leaves. Proc. R. Soc. B **265**, 147–154 (1998).

8. Deboeuf, S., Katzav, E., Boudaoud, A., Bonn, D. & Adda-Bedia, M. Comparative study of crumpling and folding of thin sheets. Phys. Rev. Lett. **110**, 104301 (2013).
9. Tallinen, T., Aström, J. A. & Timonen, J. The effect of plasticity in crumpling of thin sheets. *Nat. Mater.* **8**, 25–29 (2009).

10. Balankin, A. S. et al. Topological crossovers in the forced folding of self-avoiding matter. *Phys. A* **388**, 1780–1790 (2009).

11. Vliegenthart, G. A. & Gompper, G. Forced crumpling of self-avoiding elastic sheets. *Nat. Mater.* **5**, 216–221 (2006).

12. Balankin, A. S. et al. Fractal features of a crumpling network in randomly folded thin matter and mechanics of sheet crushing. *Phys. Rev. E* **87**, 052806 (2013).

13. Sultan, E. & Boudaoud, A. Statistics of crumpled paper. *Phys. Rev. Lett.* **96**, 136103 (2006).

14. Bayart, E., Boudaoud, A. & Adda-Bedia, M. Tuning the ordered states of folded rods by isotropic confinement. *Phys. Rev. E* **89**, 012407 (2014).

15. Lin, Y.-C. et al. Spontaneous emergence of ordered phases in crumpled sheets. *Phys. Rev. Lett.* **103**, 263902 (2009).

16. Cambou, A. D. & Menon, N. Three-dimensional structure of a sheet crumpled into a ball. *Proc. Natl Acad. Sci. USA* **108**, 14741–14745 (2011).

17. Aharoni, H. & Sharon, E. Direct observation of the temporal and spatial dynamics during crumpling. *Nat. Mater.* **9**, 993–997 (2010).

18. Lanzer, P. *Mastering Endovascular Techniques: a Guide to Excellence* (Lippincott Williams and Wilkins, 2007).

19. Witten, T. A. Polymer solutions: a geometric introduction. *Rev. Mod. Phys.* **70**, 1531–1544 (1998).

20. Vetter, R., Wittel, F. K. & Herrmann, H. J. Morphogenesis of filaments growing in flexible confinements. *Nat. Commun.* **5**, 4437 (2014).

21. Stoop, N., Najafi, J., Wittel, F. K., Habibi, M. & Herrmann, H. J. Packing of elastic wires in spherical cavities. *Phys. Rev. Lett.* **106**, 214102 (2011).

22. Najafi, J., Stoop, N., Wittel, F. K. & Habibi, M. Ordered packing of elastic wires in a sphere. *Phys. Rev. E* **85**, 061108 (2012).

23. Petrov, A. S., Boz, M. B. & Harvey, S. C. The conformation of double-stranded DNA inside bacteriophages depends on capsid size and shape. *J. Struct. Biol.* **160**, 241–248 (2007).

24. Odijk, T. Scaling theory of DNA confined in nanochannels and nanoslits. *Phys. Rev. E* **77**, 060901 (2008).

25. Pineirua, M., Adda-bedia, M. & Moulinet, S. Spooling and disordered packing of elastic rods in cylindrical cavities. *EPL* **104**, 14005 (2013).

26. Bayart, E., Deboeuf, S., Corson, F., Boudaoud, A. & Adda-bedia, M. Measuring order in the isotropic packing of elastic rods. *EPL* **95**, 34002 (2011).

27. Donato, C. C., Gomes, M. A. F. & de Souza, R. E. Crumpled wires in two dimensions. *Phys. Rev. E* **66**, 015102 (2002).

28. Donato, C. C., Gomes, M. A. F. & de Souza, R. E. Scaling properties in the packing of crumpled wires. *Phys. Rev. E* **67**, 026110 (2003).

29. Purohit, P. K. et al. Forces during bacteriophage DNA packaging and ejection. *Biophys. J.* **88**, 851–866 (2005).

30. Gomes, M. A. F., Brito, V. P. & Araujo, M. S. Geometric properties of crumpled wires and the condensed non-solid packing state of very long molecular chains. *J. Braz. Chem. Soc.* **19**, 293–298 (2008).

31. de Gennes, P. G. *Scaling Concepts in Polymer Physics* (Cornell University Press, 1985).

32. Kantor, Y., Kardar, M. & Nelson, D. R. Statistical mechanics of tethered surfaces. *Phys. Rev. Lett.* **57**, 791 (1986).

33. Stoop, N., Wittel, F. K. & Herrmann, H. J. Morphological phases of crumpled wire. *Phys. Rev. Lett.* **101**, 094101 (2008).

34. Havlin, S. & Ben-Avraham, D. Theoretical and numerical study of fractal dimensionality in self-avoiding walks. *Phys. Rev. A* **26**, 1728–1734 (1982).

35. Wood, A. J. Witten's lectures on crumpling. *Phys. A* **313**, 83–109 (2002).

36. Blair, D. L. & Kudrolli, A. Geometry of crumpled paper. *Phys. Rev. Lett.* **94**, 166107 (2005).

37. Adda-Bedia, M., Boudaoud, A., Boué, L. & Debœuf, S. Statistical distributions in the folding of elastic structures. *J. Stat. Mech.: Theory Exp.* **2010**, P11027 (2010).

38. Keller, N., Grimes, S., Jardine, P. J. & Smith, D. E. Single DNA molecule jamming and history-dependent dynamics during motor-driven viral packaging. *Nat. Phys.* **12**, 757 (2016).

39. Marenduzzo, D., Micheletti, C., Orlandini, E. & Sumners, D. W. Topological friction strongly affects viral DNA ejection. *Proc. Natl Acad. Sci. USA* **110**, 20081 (2013).

40. Sobral, T. A., Gomes, M. A. F., Machado, N. R. & Brito, V. P. Unpacking of a crumpled wire from two-dimensional cavities. *PLoS ONE* **10**, e0128568 (2015).

## Acknowledgements

This work was supported by DFG within SFB 1027 (A7). M.H. acknowledges support from NWO-VIDI grant 680-47-548/983.

## Author contributions

M.R.S., M.H. and D.B. designed the research. A.F., J.N. and M.H. performed the experiments. M.R.S. developed the SAW model and performed simulations. All authors contributed to the analysis and interpretation of the results. M.R.S. drafted the manuscript.

## Additional information

**Competing interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Shaebani, M. R. et al. Compaction of quasi-one-dimensional elastoplastic materials. *Nat. Commun.* **8**, 15568 doi: 10.1038/ncomms15568 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017