Abstract

When a straight hollow tube having circular cross-section is bent uniformly into an arc, the cross-section tends to ovalize or flatten due to the in-plane stresses induced by bending; this ovalization phenomenon is called the Brazier effect. The present paper is aimed at theoretical formulation of the Brazier effect observed in multilayered cylinders, in which a set of thin hollow cylinders are stacked concentrically about the common axis. The results indicate that mechanical couplings between stacked cylinders are found to yield pronounced suppression of the cross-sectional ovalization. Numerical computations have been performed to measure the degree of suppression in a quantitative manner and to explore how it is affected by the variations in the bending curvature, the number of stacked cylinders, and the interlayer coupling strength.

1. Introduction

Concentrically stacked cylinders are ubiquitous in nature and technology. Primary examples in the field of condensed matter physics is a series of multiwalled nanotubes. To date, successfully synthesized was a wide variety of nanotubes made from carbon [1, 2], silver [3], gold [4], silica [5], boron nitride [6, 7], WS2 [8, 9], and TiO2 [10], most of which show potential utilities as device components owing to unique physicochemical properties and/or large surface area. In a realm of biochemistry, in addition, there are several examples such as self-assembled microtubules in living cells [11] and multilayer lipid tubes [12]. In those biotubes, the mechanical strength and morphological stability are crucial factors for their biological functions [1315].

From a practical viewpoint, coaxial multicylinder structures give a merit for designing composite materials both in macroscopic and nanoscopic scales. In civil engineering, for example, deep water pipelines [16] and fiber-reinforced polymer tubes [17] are known to be endowed with high mechanical strength by constructing them with multilayered structures. Choice of multilayered structures is also preferred in building thermal energy storage with cylindrical shape [18], wherein a concentric-tube arrangement results in significant improvement in the system performance.

Another important benefit of coaxial multicylinder structures is enhanced mechanical robustness against cross-sectional deformation. Compared with a single hollow cylinder, multilayered counterparts afford resistance to cross-sectional deformation caused by external load [19, 20]. The resistance property originates mainly from the presence of internal cylinders enclosed by outer ones; the internal cylinders can push back the inward deflection of outer cylinders, thus preventing the cross-sections from severe deformation. This will hold true for bending-induced deformation in the cross-section of a hollow cylinder, commonly known as the Brazier effect [21]. Coaxial multicylinders are expected to lead feasible suppression of the Brazier effect in the systems, thus being preferred for designing composite materials that should be undeformable against external load. Despite the fundamental interest, little attempt has been made so far to elucidate the degree of suppression which multilayered structures supply.

The purpose of this work is to establish the theory that quantifies the Brazier effect suppression expected to occur in multilayered cylinders. Thin-shell theory based approach is used to formulate the strain energy of multilayered cylinders under bending, the result of which allows describing the growth of ovalization with increasing the bending curvature. The dependence of the cross-sectional shape on the relevant parameters such as the bending curvature, the number of layers, and the coupling strength between adjacent layers has been also clarified.

2. Background: What Is the Brazier Effect?

Flattening of elastic tubes under pure bending has been a focus of much attention over many decades since Brazier first studied this problem [21]. The flattening phenomenon has been usually called the Brazier effect after his name, while other seminal works by Reissner [22] and by Seide and Weingarten [23] have played decisive roles in the progress of the issue.

Figure 1 illustrates the mechanism of Brazier effect. It shows an infinitesimal portion of a deformed hollow cylinder that is originally straight and has a circular cross-section. Application of pure bending to the cylinder causes the compressive and tensile stresses that act at an angle to the unrotated section. As a result, the original circular shape deforms into an oval as displayed in the bottom panel of Figure 1. This is equivalent to saying that for a certain curvature, the ovalization minimizes the mechanical energy of the cylinder.

The degree of ovalization is quantified by the parameter , called the oblateness. Given an initially circular cross-section of radius , the product equals to the minor axis length of the oval obtained after deformation. The definition of is conventionally used even for extremely deformed sections, which go beyond simple ovalization to flattening or collapse in a peanut-like shape in which can be much less than unity.

3. Formulation

The stable cross-sectional shape of a hollow cylinder under pure bending is evaluated using a thin shell theory. The theory states that the energy cost for a cylinder to deform within linear elastic region is written by a sum of three energy terms: Here, is the strain energy associated with the circumferential displacement, is that with the axial strain, and is the energy caused by the interaction between cylindrical surfaces in the radial direction. All the three energy terms are functions of the curvature of the deformed cylinder’s axis and the displacement of volume elements. The latter quantity, displacement, is in general represented by a spatially dependent vector. To describe it, we introduce a set of mutually perpendicular vectors , , and as shown in Figure 2. Amplitudes of the vectors , , and indicate the displacements of a volume element in the radial, circumferential, and axial direction, respectively.

Our objective is a series of the optimal oblateness parameters for each th layer contained in an -layered cylinder; it is what minimizes the total energy of the system under a given bending curvature . To this aim, we first derive explicit forms of , , and for a monolayered cylinder as functions of , , and (see from Section 3.1 to Section 3.4), followed by applying the variation method to with respect to , , and . The strategy for a monolayered system can be extended straightforwardly to multilayered ones, as will be demonstrated in Section 5.

3.1. Stress-Strain Relation

From the definition of elasticity, every elastic medium satisfies an essential relation between applied strain () and the resulting stress () , which is written in terms of cylindrical coordinates by with being Young’s modulus and being Poisson’s ratio of the material under discussion. As for a hollow cylinder, it follows from (3)-(4) that Suppose that an infinitesimal volume element of a cylinder is displaced slightly in the -direction ( or ). The strain energy caused by the deformation reads We thus obtain the expressions where we renamed by and by . To proceed the argument, we need to derive explicit functional forms of , in terms of the displacements , , .

3.2. Circumferential Strain Energy

We first consider of a hollow cylinder subjected to cross-sectional deformation. Suppose a circumferential line element of length lying within the annular cross-section with thickness (see Figure 3). The tilde attached to means that it is measured at a point distant from the centroidal circle . We will see later that under appropriate conditions, the strain is determined by and that are defined just on the undeformed centroidal circle of the circular cross-section (see (16)).

The extensional strain of the circumferential line element, designated by , is defined by where and is the length of the line element after deformation (throughout the paper, the asterisk symbolizes the quantity after deformation). The coordinates , of the element after deformation are given by which imply that Note that (8) can be rewritten as . Squaring the both sides and then rearranging give For , the term can be omitted; we substitute (10) into (11) to obtain where and so forth.

The last term in (12) is associated with the rotation of the line element due to deformation. The rotation angle consists of two parts (see Figure 4): (i) a clockwise component due to the spatial variation of in the circumferential direction and (ii) a counterclockwise one due to the circumferential displacement of the element. Combination of the two parts gives which has a positive value in the counterclockwise sense.

Formula (12) is valid for arbitrary large rotation . Particularly when and are both sufficiently small (but finite), we may neglect the second term in the right side in (12) (here we exclude the possibility that or is of the order of or larger). We further assume that normals to the undeformed centroidal circle remain normal and inextensional during the deformation. Under this assumption, we can write where and denote the displacements of a point just on and is a radial coordinate measured from . By substituting (14) into (12), we attain the strain-displacement relation such as with the definitions Here is the radius of the undeformed circle .

The results (15) and (16) state that the circumferential strain at points distant from is determined by the displacements and of a point just on . The results also allow us to write the deformation energy as or, more concisely, which implies that is a functional of , , and their derivatives as

3.3. Axial Strain Energy

Parallel discussions to that in Section 3.2 give us the explicit form of the deformation energy caused by the axial strain. The strain of the axial line element is defined by Here is the length of the line element after deformation. Using the displacement components , , and , the square of is given by Note the equivalence of with in (21); then, we obtain from (20)-(21) that where and so forth. Since and , we have

In general, differs from that is allocated just on . Nevertheless, can be related with and by and accordingly . Similarly, we use the approximations of and , which imply that , . Consequently we have with the definitions of which leads us to the axial energy expression It should be emphasized that the expression (26) implies the functional property of represented by

3.4. Interaction Energy

Mechanical stability of cross-sections in multilayered cylinders owes in part to interaction between adjacent cylindrical surfaces. In typical nanotubular materials, for instance, the interaction originates from intermolecular van der Waals forces, thus sufficiently small compared with strong chemical bonding [1]; as a result, many nanotubular materials are endowed with anisotropy in mechanical stiffness, characterized by large stiffness against longitudinal stresses and flexibility to radial compression.

For a linear elastic deformation, it is reasonable to represent the interaction energy between th and th concentric cylinders by Here, is an effective spring constant per surface area , which serves as a measure for the interaction strength. The value of is dependent on the material and geometric structure of the system considered; among choices, we hereafter pay attention to the weak interaction for which or lesser, by referring to the existing multiwalled nanotubes [24].

4. Displacement Condition

To evaluate the optimal value of for each th cylinder, we need to relate with , , and , through which the total strain energy is set to be a function of . In his pioneering work, Brazier hypothesized that a typical cross-section deforms according to [21] where the parameter is a dimensionless measure of the ovalization. The suffix “inext" means that the cross-section is assumed to deform inextensionally in the direction, namely, when in (18) is set to zero. Under this special restriction, one obtains exactly the relation [25] or equivalently where the constant of integration was set to zero since there is no net rotation of the section about the cylinder axis.

The relation (30) is exact when the circumferential elongation/contraction of volume elements is forbidden . Even though such the circumferential deformation is permitted, the relation (30) nearly holds in a realm of linear elastic deformation. This is why we use the following expressions for displacement of volume elements contained in the th cylinder: where is the oblateness whose value depends on . The validity of the above expressions will be revisited in Appendix.

5. Optimal Ovalization

Discussions made in the preceding two sections allow us to formulate the total mechanical energy as a function of the series of : The optimal ovalization is determined by To calculate (35), we can omit -terms of higher order than the 2th degree since for all . The resulting equation with respect to can be reduced to the matrix representation of Here, and are column vectors defined by and is an tridiagonal matrix with the definitions of As a consequence, the optimal values of are determined by from which we derive the dependence of the optimal on the curvature as formally expressed by .

6. Numerical Results

The strategy for evaluating the optimal is summarized as follows.(i)Prepare the energy expressions of , , and given by (18), (26), and (28), respectively.(ii)Substitute the displacement expressions of , , and , given by (32)-(33), into the energy expressions.(iii)Through (i) and (ii), the total energy has been reduced to the ’s polynomial of the 4th order.(iv)Omit the 3rd and 4th order terms in the ’s polynomial and make the partial derivatives of to obtain the matrix equation (see (36)).(v)Solve (40), , to reach the optimal oblateness .

In actual computations, we consider that all the -layered cylinders with different possess the same value of ; namely, the radius of the outmost cylinder is fixed to be a common constant irrespective of . Taking the common as a unit of length, we set the layer’s thickness to be . The remaining parameters to be determined are , , and ; those values are artificially tuned in order to explore the effect of their variations on the optimal cross-section.

Figure 5 displays the stable cross-sections under the bending curvature indicated. The curvature is normalized by multiplying the outmost cylinder’s radius ; thus, is dimensionless. For every three samples, we fixed and interlayer coupling strength of . It follows from the views that the ovalization gradually evolves with increasing as expected. Particularly at , the outmost layer reaches flattened shape, while the innermost layer nearly collapses into a peanut-like shape. The two bulges indicated by dashed squares imply that large excess strain is stored in the inner layers around the bulgy regions, which prevent from significant deformation of the outer layers by pushing them outward. We thus can say that the physical origin of the Brazier effect suppression in multilayered systems is the present of additional layers embedded into the hollow core of the surrounding layers.

We have found in computations that the degree of suppression monotonically increases with the number of constituent layers . Figure 6 exhibits the dependence of on and , which tell us that the larger we apply, the smaller we obtain. It is also observed that for , the suppression no longer works effectively and the optimal converges to a constant.

Added to , the interlayer coupling quantified by is another important factor for determining the stable cross-section. Figure 7(a) gives examples on how the variation in contributes to the suppression of the Brazier effect. The horizontal axis indicates the interlayer spring constant normalized by multiplying , by which it is nondimensional. The sample size is fixed to be , and the change in with increasing is plotted. The results show that the degree of suppression is magnified with intensifying the couplings up to , above which an increase in gives minor contribution. The pronounced suppression at the small region is observed universally as far as the sample is composed of modest number of layers, say . This fact is demonstrated graphically in Figure 7(b), which shows the difference, designated by , in the short radius of deformed cross-sections between the weak and strong coupling conditions defined as follows: The plot makes clear that the cross-section endowed with strong interlayer couplings tends to retain its original circular shape when . Figure 8 displays, as examples, the robust behaviors of samples with and , together with collapse of those with and . Inward convex bulges observed in the latter two samples disappear in the former two, which is a manifestation of the coupling effect in modest--layered systems.

7. Conclusion

We have proposed an approximation theory that describes the bending-induced ovalization in the circular cross-section of multilayered cylinders. Taking into account the in-plane deformation within individual layers and interlayer couplings, we have formulated the mechanical energy of multilayered cylinders, followed by variational computation to evaluate the optimal oblateness that minimizes the energy under pure bending. The oblateness was found to decrease monotonically with increasing the number of layers and/or strength in the interlayer coupling, as we demonstrated in a systematic and quantitative manner. The results will facilitate the basic understanding of the mechanics of coaxial multicylinder structures that are present universally in macroscopic and nanoscopic science.

Appendix

On the Validity of (32)

In this Appendix, we examine the validity of our approximation as restated by which describe the cross-sectional deformation. We should be reminded that the above expressions of and are natural extensions of the following: which are exactly correct under the ideal condition that the circumferential strain (i.e., in (18)) must be zero. Therefore, it is nontrivial if we are allowed to make the straightforward extension to the case, where , just by replacing the coefficient from to with keeping the function forms of and . We will demonstrate below that the extension should be justified under the numerical conditions we have applied.

For a single-layered cylinder (i.e., ), the fourth-order polynomial of given by (35) reads as where and . The spring constant is that associated with springs illustrated in the left panel of Figure 9. We can prove that the third-order equation, , with respect to has a real solution and two imaginary ones as explained by the right panel in Figure 9. The real solution has the form of with and , , , to five significant digits.

This -dependence of comes from our assumption of (32). Hence, the validity of our assumption can be tested by comparing with the exactly optimal oblateness of and that are defined by Equation (A.6) is a general expression of displacements, free from our assumption, thus being accurate no matter if vanishes or not. The optimal oblateness based on (A.6) is derived from and , where is a polynomial with respect to and , having the function form different from given by (A.3).

Figure 10 shows the comparison between and for three different values of . An almost equivalence between and has been confirmed for every . Besides, was found to be nearly equal to at all and (not shown). From the data comparison, we conclude that the displacement expressions of (32) are effective for evaluating the oblateness of cylinders that exhibit nonvanishing .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by the University of Incheon, International Cooperative Research Grant, in 2013.