On structured surfaces with defects: geometry, strain incompatibility, internal stress, and natural shapes
Abstract
Given a distribution of defects on a structured surface, such as those represented by 2dimensional crystalline materials, liquid crystalline surfaces, and thin sandwiched shells, what is the resulting stress field and the deformed shape? Motivated by this concern, we first classify, and quantify, the translational, rotational, and metrical defects allowable over a broad class of structured surfaces. With an appropriate notion of strain, the defect densities are then shown to appear as sources of strain incompatibility. The strain incompatibility relations, with appropriate kinematical assumptions on the decomposition of strain into elastic and plastic parts, and the stress equilibrium relations, with a suitable choice of material response, provide the necessary equations for determining both the internal stress field and the deformed shape. We demonstrate this by applying our theory to KirchhoffLove shells with a kinematics which allows for small insurface strains but moderately large rotations.
Keywords: 2dimensional materials; thin structures; geometry of defects; surface dislocations; surface disclinations; nonmetricity; strain incompatibility; residual stress.
Mathematics Subject Classification (2010): 74E05; 74K15; 74K20; 74K25; 53Z05.
1 Introduction
The aim of this article is to study geometry and mechanics of defects in structured surfaces. The term structured surface is used to represent a variety of 2dimensional material surfaces such as 2dimensional crystals, with intrinsic translational, rotational, and metrical order (colloidosomes, carbon nanotubes, graphene etc.); thin sandwiched structures; liquid crystalline membranes and shells, with intrinsic crystalline order (singlelayer viral capsids) or without (nematic membranes, single layers in smectics and cholesterics); and Cosserat surfaces, which are used to model a hierarchy of plate and shell theories for thin elastic structures abundant in structural engineering applications. The defects are anomalies within the local arrangement of entities in an ordered structure where the order is usually defined in terms of translational, rotational, and metrical symmetries of the underlying material. These anomalies are omnipresent in nature, e.g., 2dimensional materials such as graphene are known to contain edge dislocations (translational anomalies), wedge disclinations (rotational anomalies), and point imperfections (metric anomalies) such as vacancies and selfinterstitials; on the other hand, twist disclinations are commonly observed in lipid membranes. More examples are given in the following section as well as in an extensive review of the subject in [7]. The phenomena of thermal deformation and biological growth can also be categorized as those leading to metric anomalies. Many of the superior physicochemical properties of the 2dimensional defective structures can be attributed to the internal stress fields resulting from the distribution of defects [76], and also, unlike 3dimensional bodies, due to their lower dimensionality, to their ability to relax by acquiring a variety of natural (stressfree) shapes, for instance, the wavy edges of growing leaves [43], the topological corrugations present on human brain [71], helical strands of DNA [20], among other examples. The present work is concerned with the central problem of formulating a general theory that takes under its ambit the geometric characterization of these multifarious 2dimensional defective structures and also the determination of their internal stress fields and deformed shapes.
NonEuclidean differential geometry has been established to provide the necessary mathematical infrastructure to describe the geometric nature of defects in 3dimensional solids, as well as to provide a rightful setting to discuss the related issues of strain incompatibility and residual stress distribution [38, 6, 2, 3, 52, 40, 16, 14]. Despite this success in 3dimensions, the problem in lower dimensional structures is relatively less developed, primarily due to the complex interplay between the embedding geometries in the physical space, and the unavoidable nonlinearities involved in the deformation as well as the constitutive response of 2dimensional matter. We note the initial attempts made by Eshelby [28, 27] where analytical solutions for internal stress were derived for linearly elastic plates containing isolated screw and edge dislocations. This work was extended in several directions by Chernykh [11] and Nabarro [49, 50], among others [61, 48, 63]. A theory of continuous distribution of defects over thin structures was first developed by Povstenko [54, 55] by drawing analogies from the nonEuclidean description of continuous distribution of defects in 3dimensional elastic bodies. Povstenko introduced the notions of insurface dislocations (characterized by the insurface torsion tensor), disclinations (characterized by the insurface curvature tensor), and metric anomalies (characterized by the insurface nonmetricity tensor). He also provided the nonlinear conservation laws for all the insurface defect density fields as direct consequences of the BianchiPadova relations in two dimensions.
The compatibility conditions for strain fields in geometrically linear and nonlinear shells, and Cosserat surfaces, are extensively discussed in existing works [37, 44, 56, 22]. The presence of defects, however, introduces incompatibility in strain fields. The strain incompatibility equations for both linear elastic plates and vonKármán plates, with insurface anomalies, have been derived by Zubov and Derezin [77, 78, 17, 18, 75]; these also include solutions of certain special boundaryvalueproblems for determining stress and natural shape of the defective plate. The strain incompatibility equations have also appeared in the recent works on nonEuclidean elastic plates [21, 34], and growth and morphoelasticity of thin biological materials [46, 47, 43, 42, 19]. Without explicitly incorporating defect densities, these works consider a nonEuclidean metric, representing the distribution of plastic/growth strain field, and use the Riemannian curvature of this metric, which is the measure of strain incompatibility, along with the strain decomposition, to pose boundaryvalueproblems primarily for determining natural shapes.
The concepts of material uniformity, material symmetry, and inhomogeneity in elastic Cosserat surfaces, following the pioneering works of Noll [52] and Wang [72], are also firmly established [26, 73, 74, 23, 24, 25], although these works have neither attempted to describe the inhomogeneity distribution in terms of the curvature and nonmetricity (the notion of torsion does appear in some of these works), nor have they discussed the relevant issue of strain incompatibility. A theory of materially uniform, inhomogeneous (dislocated) thin elastic films, derived from a 3dimensional uniform, inhomogeneous (dislocated) elastic body, has been recently proposed by Steigmann [68], and applied to determining the natural shapes of plastically deformed thin sheets [15]. Finally, we mention, only in passing, the extensive work on mechanics of topologically defective (‘geometrically frustrated’) liquid crystalline surfaces [58, 7, 8, 10], which, in contrast to the local theories mentioned above, have taken a distinguished localglobal (geometricaltopological) standpoint in describing the nature of defects.
There is a clear lack of a complete nonEuclidean geometric characterization of continuously distributed material defects in 2dimensional structured continua. While these certainly have analogous descriptions in the 3dimensional theory, there is a considerable richness in the description of the allowable defects as well as their geometrical properties for the 2dimensional structure. Additionally, there are no derivations of strain incompatibility relations for sufficiently general kinematic and constitutive response as afforded by most of the known 2dimensional materials. With this in mind, we present a theory, within the natural setting of nonEuclidean differential geometry, that on one hand unifies the several seemingly different streams of research discussed above, and also provides a rigorously constructed, sufficiently general, framework for studying a large range of problems associated with geometry and mechanics of defective structured surfaces. In particular, we give a complete nonEuclidean characterization of all the translational, rotational, and metric anomalies in structured surfaces, derive the imposed restrictions from BianchiPadova relations, and establish general strain incompatibility relations. To illustrate our theory, we consider the specific case of KirchhoffLove shells and provide a framework, involving kinematics, additive decomposition, incompatibility relations, and balance laws, for posing complete boundaryvalueproblems for determining internal stress and deformed shapes for a class of 2dimensional continuously defective structures undergoing small stretch but moderately large deformation.
A brief overview of the paper is as follows. In Section 2, we provide several illustrative examples to demonstrate the nonEuclidean character of local material defects in structured surfaces. Motivated from Section 2, we begin Section 3 by introducing the notion of material space, which includes a 2dimensional body manifold, a nonRiemannian material connection, and a material metric, as our prototype to characterize continuously defective structured surfaces. Both insurface and outofsurface material anomalies are taken into account, and are identified with the components of the tensors of nonmetricity, torsion, and RiemannChristoffel curvature of the material connection (see Table 1). Subsequently, by exploiting BianchiPadova relations, we obtain several restrictions, both as algebraic relations and differential equations, on these components while emphasizing their interdependence (see Table 2). We conclude the section by introducing a Riemannian structure on the material space induced by the material metric. We obtain geometric relations which connect the curvatures of the Riemannian and nonRiemannian spaces. These relations are central to our formulation of strain incompatibility equations in Section 4. A generalized notion of strain is introduced which defines the kinematical nature of our structured surface. The strain incompatibility relations lead us to pose complete boundaryvalueproblems for the determination of internal stress fields and deformed shapes of defective structured surfaces. This is illustrated in Section 5 by restricting our attention to KirchhoffLove shells. We also postulate an additive decomposition for the insurface and bending strain fields into elastic and plastic parts, while arguing that this separation of order is sufficiently general to accommodate small insurface strains with moderately large rotations. The plastic strain fields are to be solved using the incompatibility relations with a given distribution of defects. We show that several existing formulations follow as special cases, in particular the FöpplvonKármán equations for continuously defective thin elastic isotropic plates and the shape equations for continuously defective thin isotropic fluid films. We conclude our study in Section 6.
2 Nature of surface defects
In this section, we provide several illustrative examples of defects in structured surfaces. The defects are understood as anomalies within the local arrangement of entities in an ordered structure where the order is usually defined in terms of rotational, translational, and metrical symmetries of the underlying material. Defects can also appear as global anomalies which affect the topology of the surface, such as those present in multiply connected and nonorientable surfaces [31, 32, 7]; these are however not discussed in the present work. The following examples are presented with an intent to emphasize the nonEuclidean geometric nature of the defects as is incorporated in the subsequent sections. In particular, the central idea of our work of embedding the structured surface within a 3dimensional nonRiemannian geometric space emerges naturally as we proceed through these rudimentary illustrations.
The rotational anomalies in a structured surface appear in the form of disclinations. Depending on the material nature of the surface, rotational order can be present due to intrinsic crystallinity of the surface (such as in colloidosomes, singlelayer viral capsids, carbon nanotubes, and graphene) or due to an extrinsic orientation field (such as in nematic membranes, single layers in smectics, and cholesterics) [7]. As a result, we distinguish between rotational order, or lack thereof, appearing intrinsically and extrinsically in a surface. We also note that unlike disclinations in 3dimensional crystalline solids, which have large formation energy and hence are rarely observed [4], disclinations in 2dimensional crystals are omnipresent since the surface can now relax the energy by escaping into the third dimension. Isolated disclinations in structured surfaces without intrinsic crystalline order are shown in Figures 1(a,b). The rotational order is here present due to a director field distribution, denoted by , over a planar domain parametrized by Cartesian coordinates . The director field in Figure 1(a) is restricted to lie strictly in the plane; it may represent a deformed configuration of a nematic membrane or a single layer in the cholesteric phase of some liquid crystalline material. In contrast, the directors in Figure 1(b) are allowed to orient themselves transversely to the plane; this can model either a lipid monolayer where the director orientation represents the orientation of individual lipid molecules, or a single layer of molecules in the smectic A or C phase [35]. In nematics, smectics, and cholesterics, is identifiable with due to the mirror symmetry about the midorthogonal plane of the director axis. The lack of intrinsic crystalline order (translational and rotational), within the plane, in these examples can be primarily attributed to viscous relaxation [36]. Disclinations in such structured surfaces can be characterized by the signed angle through which the director rotates upon circumnavigating along a loop over the surface. The Frank vector of the disclination is a precise measure of this signed angle. A disclination is of wedge or twist type depending on whether is transverse or tangential, respectively, to the surface. The disclination in Figure 1(a) is of wedge type with Frank vector and the one in Figure 1(b) is of twist type with Frank vector . Here, the triple denote the standard basis of the Cartesian coordinate system (. Note that the wedge disclination line in Figure 1(a) and the twist disclination line in Figure 1(b) are both along the axis. Disclinations can also appear in surfaces with intrinsic crystalline order, e.g., an ordered arrangement of lattice sites where the directors are attached in viral capsids or hexagonal lattice structure of the carbon atoms in graphene sheets. As illustrated in Figure 1(c), circumnavigating along a loop encircling the disclination, a lattice vector rotates through an angle which is an integral multiple of one of the rotational symmetry angles of the lattice. The wedge disclination located at , in the 2dimensional hexagonal lattice in Figure 1(c), is characterized by its Frank vector . Material surfaces can also possess twist disclinations in the form of local intrinsic orientational anomalies, which correspond to breaking of the reflectional symmetries of the 2dimensional material with the local tangent plane of the surface as the mirror plane, e.g., hemitropic plates [65, 25]. They are represented mathematically as illdefined (multivalued) local orientation field over the surface.^{1}^{1}1An example of a global intrinsic orientational anomaly would be the global orientational anomaly present in a surface Möbius crystal due to its nonorientability. Note that, in order to quantify the disclinations discussed so far, the loop of circumnavigation is restricted always within the surface. All the disclinations shown in Figures 1, as well as the intrinsic orientational anomalies discussed above, are quantified using an insurface loop . The case otherwise can appear in 2dimensional homogenized models of thin 3dimensional multilayered structures, e.g., a stack of few monolayers of smectics or cholesterics, thin multiwalled nanotubes, or a thin slice of some 3dimensional oriented media. In these structures, disclinations may appear over the representative base surface (often the ‘midsurface’ of the layered structure) as the homogenized or effective rotational anomaly of all the distributed disclinations across the thickness of the thin structure. In describing these disclinations, the loop of circumnavigation must be taken transversely to the base surface, see Figure 2. Depending on the direction of the resulting vector of angular mismatch, these disclinations may either be of wedge or twist type.
The translational anomalies are represented by dislocations. The nature of dislocations in 2dimensional matter is analogous to that in 3dimensional materials. Isolated edge and screw surface dislocations are shown in Figures 3(a) and 3(b), respectively, within a 2dimensional cubic lattice along with the Burgers parallelograms. The Burgers vector, defined as the closure failure of the Burgers parallelogram, is tangential to the surface of the lattice in the former case and transverse in the latter. In these examples, the dislocations appear essentially due to the breaking of the intrinsic translational symmetries of the 2dimensional matter. On the other hand, in thin multilayered structures or thin slices of oriented media, dislocations may be present, irrespective of the crystallinity of the material, as a result of either an ordermismatch of individual layers within the stack or as a homogenized or effective limit of all the distributed dislocations within the 3dimensional slice. The Burgers parallelogram is, naturally, transverse to the representative midsurface of the stack, in contrast to the examples shown in Figures 3(a,b). The precise type of these dislocations, edge or screw, can be determined from the direction of the Burgers vector. An edge dislocation in a layered medium is shown in Figure 3(c), arising due to the presence of a sandwiched semiinfinite layer between two infinite layers of material [41, Ch. VI].
The metrical anomalies bring about ambiguity in the (local) notion of “length” and “angle” over the surface. Metric anomalies are generated due to intrinsic point imperfections such as vacancies and selfinterstitials, see Figure 4(a), as well as a result of insurface thermal deformation and biological growth, see Figure 4(b). Note that foreign interstitials fall within the realm of materially nonuniform bodies (e.g., functionally graded materials) where the material constitution changes from point to point; their consideration is outside the scope of the present work. If the distance between the constituent entities in a lattice structure is measured by counting lattice steps, the presence of point defects, such as a vacancy or a selfinterstitial, clearly introduces ambiguity in this step counting [40]. Apart from these pure insurface metric anomalies, differential growth (or thermal deformation) across the thickness direction within a thin multilayered structure may result in transverse metric anomalies within an appropriately homogenized 2dimensional theory, see Figure 4(c).
The simple examples described above are sufficient to motivate the nonEuclidean nature of the defects. Recall that, in order to quantify disclinations, we required circumnavigation of a vector along a loop and rotational mismatch between the initial and the final orientation of the vector. These notions correspond, respectively, to parallelly transporting a vector with respect to an affine connection and to the RiemannChristoffel curvature associated with the affine connection. The Frank vector uniquely characterizes the RiemannChristoffel curvature tensor [2]. The affine connection has to be necessarily nonEuclidean, since the director fields leading to disclinations are clearly not parallel in the Euclidean sense. Moreover, as the directors may point outside the surface, a differential geometric description of disclinations in structured surfaces would necessarily require embedding the surface into a 3dimensional space with a specific nonEuclidean connection. In the case of dislocations, the closure failure of the Burgers parallelogram is analogous to the notion of torsion of an affine connection over a manifold which characterizes closure failures of infinitesimal parallelograms [38, 6]. Finally, the metric anomalies are characterized by the nonmetricity tensor, which quantifies the nonuniformity of the metric tensor with respect to an appropriate affine connection [3, 39]. Motivated with these geometric analogies, we are now in a position to pursue a systematic study of geometry of defects in a structured surface.
3 Geometry of surface defects
The mathematical prototype for structured surfaces is a connected, compact 2dimensional manifold , possibly with boundary, which is embeddable (as a topological submanifold) in . Examples of such manifolds, in the orientable category, are sphere, sphere with a finite number of handles added, twisted bands with twists for integers etc., and in nonorientable category, twisted bands with twists for integers , e.g., a Möbius band for which is zero. We can add boundaries to these manifolds by removing finite number of open discs. The condition of embeddability in precludes Klein bottle like surfaces and real projective planes. Our prototype manifold is topologically characterized by its orientability, twistedness, Euler characteristic, the number of open discs removed, i.e., the boundaries, and other topological invariants. We will call the body manifold. A fundamental theorem in differential topology (Tubular Neighbourhood Theorem [9, Theorem 11.4]) guarantees the existence of a tubular neighbourhood of in , for sufficiently small . Here, denotes the minimum Euclidean distance of from . As a bounded open set in , naturally admits a manifold structure, with as an embedded submanifold. Existence of induces a vector bundle (the normal bundle) structure over [9], which entails a vector field defined over . Our choice of , naturally endowed with a director field , is therefore appropriate for modelling structured surfaces. The differential structure, and all the fields to be defined over and , including , is assumed to be as smooth as the context demands.
Our strategy for characterising material defects on a structured surface is to first equip with a geometrical structure by associating with it a metric and an affine connection. This is then used to induce an appropriate nonRiemannian geometrical structure over , where various fundamental geometric objects, such as nonmetricity, torsion, and curvature, are interpreted as defect density measures. The induced metric and connection on is sufficient to encode all the information about the material structure of the structured surface. The BinachiPadova relations are used to obtain several restrictions on defect density fields. With these relations, it is emphasized that the various defect densities are in fact dependent on each other. The metric associated with is also used to induce a Riemannian structure over . The relationship of the curvature tensor, associated with the affine connection, with RiemannChristoffel curvature tensor, obtained from the metric, is derived. These relations will provide the starting point for deducing strain incompatibility equations in the following section. They also lead to the well known local conditions under which is isometrically embeddable into , a notion that is related to compatibility of the strain fields.
In rest of the paper, lowercase Greek indices , , etc. take values from the set and lowercase Roman indices , , etc., from the set . Einstein’s summation convention hold over repeated indices unless specified otherwise. Round and square brackets enclosing indices indicate symmetrization and antisymmetrization, respectively, with respect to them. The superscript is used to denote the inverse of an invertible matrix, whereas the superscript is used to denote the transpose.
3.1 Geometry on induced from the nonRiemannian structure on : the material space
Let the 3dimensional embedding manifold be equipped with an affine connection and a metric . Consider a chart of with such that the coordinates defined over lie along with at . Such a coordinate system is called adapted to . The restriction of the natural basis vector fields over to will be denoted by , i.e., , hence is transverse to . The coefficients of and the covariant components of are denoted by and , respectively, with respect to . The covariant derivative of a sufficiently smooth vector field , , with respect to , is denoted by
(1) 
The notation is used for the surface covariant derivative of a tangent vector field , , with respect to the projection of on , i.e., a connection with coefficients ,
(2) 
Here, the subscript denotes ordinary partial derivative with respect to . A vector field along a curve over is called parallel with respect to if, and only if, its covariant derivative along the curve vanishes identically.
The body manifold , equipped with connection and metric from the embedding space , forms the material space of the structured surface. We will call the material connection and the material metric. The “material” nature of these mathematical objects is due to the fact that the geometric quantities derived from and , when restricted to , represent various material inhomogeneities or defects within the material structure of the structured surface. As we will see immediately below, the nonmetricity tensor is a measure of distributed metric anomalies, the torsion tensor is a measure of distributed translational anomalies (dislocations), and the RiemannChristoffel curvature tensor is a measure of distributed rotational anomalies (disclinations). Most importantly, we assume to be such that the nonmetricity, torsion, and curvature tensors associated with are uniform in the coordinate and equal to their respective values at , i.e., at . This assumption alludes to the applicability of our model to thin multilayered structures, or thin slices of defective media, represented as homogenized 2dimensional surfaces. It should also be noted that we are only looking at local defects and not the ones which could arise out of various topological anomalies for multiply connected and nonorientable surfaces.
3.1.1 Nonmetricity of the material connection: metric anomalies
The third order nonmetricity tensor of the material space, measuring nonuniformity of the metric with respect to the connection , has covariant components defined as
(3) 
The negative sign in the definition is purely conventional. We assume that . The pure insurface components provide measure for the distributed surface metric anomalies, whereas components , with either of , or taking the value 3, indicate the presence of outofsurface metric anomalies, e.g., thicknesswise growth. A nonzero lead to variation in angle between tangent vectors during parallel transport with respect to the projected connection , see Figure 5(a). Indeed, the inner product of two tangent vectors and , where , changes under parallel transport with respect to from the initial point to any generic point , along some parametrized curve lying over , by the amount
(4)  
Here, we have used, and throughout , as they are parallelly transported fields along , where denotes the ordinary derivative of with respect to its argument. In structured surfaces, as we have earlier discussed in Section 2, this variation in inner product, characterized above in terms of a nontrivial , may arise from a distribution of point imperfections in the arrangement of molecules or atoms over the surface, e.g., vacancies and selfinterstitials in 2dimensional crystals, inserting (or removing) a lipid molecule into (or out of) a crystalline arrangement of identical molecules over a monolayer, thermal deformation of the surface, biological growth of cell membranes, leaves etc. The remaining components and measure the nonuniformity of the material metric in the direction, i.e., along the thickness of the structured surface, and the change in length of transverse vectors along the surface, respectively, see Figures 5(b) and 5(c). These provide faithful representations for differential growth along the thickness in thin multilayered structures discussed in Section 2 and illustrated in Figure 4(c).
3.1.2 Torsion of the material connection: dislocations
Consider two tangent vectors at some point on . Translating parallelly along and along with respect to , we obtain the vectors
(5) 
respectively. The closure failure of the parallelogram is given by (see Figure 6)
(6) 
where the functions
(7) 
constitute the components of the thirdorder torsion tensor (antisymmetric in the lower indices) over . Let . We assume that , which in turn is same as . Associated with the torsion tensor, we have the secondorder axial tensor
(8) 
Here, , where is the 3dimensional permutation symbol and . For later use, we define . The components provide measures for a variety of dislocation distributions over the structured surface. Taking (i.e., and tangential to , see Figure 6(a)), and comparing with Figures 3(a,b), it is immediate that
(9) 
represent a distribution of insurface edge dislocations and
(10) 
a distribution of insurface screw dislocations (cf. [54, 55]). Next, taking (i.e., tangential and transverse to , see Figure 6(b)), and comparing with Figure 3(c), it is evident that the components represent the outofsurface dislocations in thin multilayered oriented media such as those discussed in Section 2.
3.1.3 Curvature of the material connection: disclinations
The components of the fourth order RiemannChristoffel curvature tensor of the material connection are given by
(11) 
The functions measure, in the linear approximation, the change that a vector, , , suffers under parallel transport with respect to along an infinitesimal loop based at and lying within :
(12) 
where are the components of the initial vector with respect to the basis ; the integral represents the infinitesimal area bounded by the loop . The above formula in fact holds true for any general loop (not necessarily infinitesimal) in . We define the purely covariant components by lowering the fourth index with the material metric as Clearly, and . Moreover, as we did for nonmetricity and torsion tensors, we assume .
It is useful to decompose the components into skew and symmetric parts [55]
(13) 
where
(14) 
are components of the secondorder tensor field and the thirdorder tensor field . They represent, respectively, the skew part and the symmetric part of with respect to the last two indices. A geometric interpretation of these two fundamental tensors is as follows (see Figure 7). Let the infinitesimal loop in (12) be based at . Then the change that a vector undergoes when parallelly transported along , in the linear approximation, can be characterized by a secondorder tensor , i.e., , where
(15) 
Here, is a measure of the infinitesimal area bounded by and its unit normal. The first term in the above expression is skew with axial vector . It represents the rotation that has experienced under parallel transport about the axis , for each fixed , probed by the three Euler angles . Thus, is the measure of the rotation of about the axis . The second term , on the other hand, is symmetric; it represents a stretching, with the three principal values of the tensor as measures of the stretch along their respective (linearly independent) principal directions. The tensor can be shown to be related to the metrical properties of as it gives rise to a smeared out anomaly within the material structure which causes elongation or shortening of material vectors under parallel transport along loops (see [60] for details), as shown in Figure 7(a). We will assume in rest of the paper since, at present, we do not know of any defects in 2dimensional materials which they would otherwise represent. Some consequences of this assumption will be discussed in the next section. The curvature tensor is then fully characterized in terms of the nontrivial independent components , i.e., the secondorder tensor .
We distinguish between two families of local rotational anomalies characterized by . Consider, first, the infinitesimal loop completely lying within , see Figure 7(b). Then, the and indices in can assume only values 1 and 2, and the resulting angular mismatch after parallel transport of arbitrary vectors is characterized by three fields
(16) 
These provide a measure for the distributed rotational anomalies within the material structure of the base manifold . Drawing analogy with Figure 1, it is clear that the outofsurface component provides a measure for the density of distributed wedge disclinations over the structured surface, see Figures 1(a,c), irrespective of its crystallinity, whereas the insurface components characterize either the distributed intrinsic orientational anomalies, in case of intrinsically crystalline surfaces, or distributed twist disclinations, in case of directed surfaces (as shown in Figure 1(b)). Next, we consider , based at , to lie transversely to , see Figure 7(c). Then one of the indices and in will take the value 3, and the resulting angular mismatch after parallel transport of arbitrary vectors is characterized by the remaining six independent components of :
(17) 
Recalling our discussion in Section 2 on disclinations in thin multilayered structures of oriented media, see also Figure 2, we conclude that these components provide a measure for a variety of homogenized/effective rotational anomalies of the distributed disclinations across the thickness of the multilayered structured surface. Out of these six functions, and are of wedge type, and , , and are of twist type. As we will see shortly, these functions are in fact dependent on each other in very thin monolayer structures where the dislocation densities vanish altogether.
Geometric objects  Defect densities  
Insurface metric anomalies; Figures 4(a,b) and 5(a)  
and  Outofsurface metric anomalies; Figures 4(c) and 5(b,c)  
Insurface edge dislocations; Figures 3(a) and 6(a)  
Insurface screw dislocations; Figures 3(b) and 6(a)  
Outofsurface dislocations; Figures 3(c) and 6(b)  
Insurface wedge disclinations; Figures 1(a,c) and 7(b)  


Disclinations associated with transverse loops; Figures 2 and 7(c)  
Metrical disclinations; Figure 7(a) 
We have summarized the set of all defect densities in Table 1.
3.1.4 BianchiPadova relations
The tensors of nonmetricity, torsion, and curvature of a nonRiemannian space cannot be arbitrary due to geometric restrictions. Besides the restrictions , , and , which follow from their definitions, they satisfy the following system of differential relations, known as the BianchiPadova relations [62, p. 144]:
(18a)  
(18b)  
(18c) 
In the above expressions, antisymmetrization with respect to three indices is defined as
(19) 
The enclosed indices within two vertical bars in the subscript are to be exempted from antisymmetrization. Clearly, and (no summation on ). Additionally, there is a fourth BianchiPadova relation [62, p. 145], purely algebraic in nature, based on the following identity satisfied by the components of any fourthorder tensor with :
(20) 
After substituting relations (18a) and (18c) into (20), it boils down to an expression for in terms of , , , and . For a torsionfree, metriccompatible connection (i.e., a LeviCivita connection), this implies the familiar symmetry . However, as shown below, this particular symmetry is achieved in very thin structured surfaces under much less restrictive conditions. The first three BianchiPadova relations, restricted to a surface, have been considered previously by Povstenko [55], but without studying any of the implications, some of which are noted below.
Consequences of the first BianchiPadova relation: Equation (18a) is nontrivial only when at least one of the indices , and assume the value 3, since otherwise . Recalling our assumption that that is uniform with respect to the coordinate, (18a) reduces to
(21) 
Furthermore, if we assume that the structured surface is sufficiently thin with no dislocations associated with the transverse Burgers parallelograms, i.e., (the insurface dislocations can still be present), then (21) simplifies into a system of algebraic equations:
(22) 
For , we obtain , since (from ). This is equivalent to , or
(23) 
For , (22) can be rewritten as , or equivalently
(24) 
Combining the above two relations we can therefore infer that, for vanishing , the disclination density tensor is symmetric. Moreover, due to (16), , i.e., the pure insurface disclination densities (which may either characterize densities of twist disclinations in directed surfaces or intrinsic orientational anomalies in hemitropic surfaces) should be identical to the wedge disclination densities associated with transverse loops, e.g., in multilayered surfaces as discussed in Section 2; in particular, they should vanish in sufficiently thin structured surfaces, e.g., in 2dimensional crystals, where both and will be absent. We note that, in contrast, for 3dimensional solids, the symmetry of the disclination density tensor is implied only under vanishing of the full torsion and the nonmetricity tensor. It is worthwhile to reemphasize that the assumption is realistic only in sufficiently thin structures (biological membranes, graphene sheets, etc.), which, otherwise, can support only surface edge and screw dislocations (characterized by ).
On the other hand, if we consider multilayered or moderately thin structures of oriented media, where the assumption of vanishing is no longer physical, and assume that they do not contain any disclinations and metric anomalies, and also that and are small (of the same order), then (21) yields
(25) 
This is a conservation law for the type dislocations enforcing that they must always form loops or leave the surface. In either case, whether the dislocations are absent or not, there is no restriction on the distribution of insurface dislocations . This again is in contrast to 3dimensional solids, where the first BianchiPadova relation provides a conservation law for all dislocation densities [55, 60].
Consequences of the second BianchiPadova relation: Equation (18b), in the absence of both type dislocations and metric anomalies (), in addition to , reduces to a simple conservation law
(26) 
to be satisfied by disclinations characterized by , as well as owing to the symmetries (Equations (23) and (24)), and surface edge dislocations. Assuming that and are small, and of the same order, we obtain
(27) 
These are linear conservation laws for the respective disclinations, requiring their lines to either form loops or leave the surface. Note that there is no restriction on (wedge disclinations), in contrast to what one would expect for 3dimensional solids.
Consequences of the third BianchiPadova relation: We use (18c) to obtain a simple representation for nonmetric tensor. With , (18c) can be rewritten as
(28) 
It can be shown by direct substitution that a nontrivial solution of (28) is given by
(29) 
where are arbitrary symmetric functions over . It is a consequence of the fundamental existence theorem of linear differential systems that in absence of disclinations (i.e., ) over a simply connected (hence ), if the matrix field is positivedefinite for symmetric functions , then is the only solution to (28) over . This result is proved in [60]. As the density of metric anomalies is assumed to be uniform with respect to the coordinate, we will interpret this representation of the metric anomalies in absence of disclinations over simply connected patches over as
(30) 
The symmetric matrix field is known as quasiplastic strain [3]. In absence of disclinations, the positivedefinite symmetric matrix field can be used to define an auxiliary material space , equipped with the original material connection but a metric . The nonmetricity of the auxiliary material space vanishes identically by definition. The secondorder tensor field , where , characterizing pure insurface metric anomalies in the absence of disclinations, has the unique decomposition
(31) 
where is the trace of and is the deviatoric part of (i.e., ). The first term represents isotropic metric anomalies and the second represents anisotropic metric anomalies [60]. When is purely isotropic, i.e., , it is straightforward to obtain , where . The surface metric of the auxiliary material space for isotropic metric anomalies is, hence, conformal to the surface metric of the original material space, . This formulation is readily applicable to model various reallife surface metric anomalies such as 2dimensional anisotropic biological growth, thermal expansion, distributed point defects, etc.
Consequences of the fourth BianchiPadova relation: The fourth BianchiPadova relation imposes interdependence on the disclination density measures , out of which the interdependence between the two distinct families of disclinations characterized by and , derived in the following, are in particular interesting. Assuming , the insurface components of (20) require , since , which is the trivial relation . Next, if we also assume that the metric anomalies are absent, i.e., , then (20), together with (18a), yields
(32) 
Here, . After substituting , as per our assumption on , and , or equivalently , the right hand side of the above relation vanishes identically, thereby enforcing the symmetries
(33) 
In terms of disclination densities, these are, respectively, and . Interestingly, we reached the same conclusion from the first BianchiPadova relation. We will of course obtain a nontrivial consequence of the fourth BianchiPadova identity whenever .
The results of this section are summarized in Table 2.

Implications on defect densities  










3.2 The induced Riemannian structure
The coefficients of any general affine connection of a manifold , with nontrivial torsion and nonmetricity , can be decomposed as [62, p. 141]
(34) 
where the functions are coefficients of the LeviCivita connection (torsionfree, metriccompatible) induced by the metric :
(35) 
with , and
(36a)  
(36b)  
(36c) 
The functions form the components of the contortion tensor, whereas the tensor is an equivalent measure of nonmetricity. The covariant components
(37) 
of the RiemannChristoffel curvature tensor of the LeviCivita connection and components of the material curvature are related as [62, p. 141]
(38) 
where and denotes covariant differentiation with respect to the LeviCivita connection . From the general symmetry relations, , of the Riemannian curvature induced by a metric, it is evident that it has only six independent components characterized by , and . The only nontrivial relations out of (38), when restricted to , are
(39a)  
(39b)  
(39c) 
Here, , , , and denotes covariant differentiation with respect to the projected LeviCivita connection on , consisting of components . The relations (39) are central to the theory of mechanics of defects as they are directly related to the strain incompatibility equations which we discuss next. Indeed, once we have identified the material metric in terms of strain fields on the structured surface, (39) constitute a system of PDEs for the strain fields, with defect densities as source terms. It should be noted that the components do not appear in any of the equations (39). This is because, according to the fourth BianchiPadova relation (32), they can be written in terms of and other defect measures, and hence are not independent quantities.
4 Strain incompatibility relations for structured surfaces
In this section, we begin by introducing the notion of strain for a structured surface. The complete set of strains represent essentially the kinematical nature of shell theory that is being employed to describe structured surfaces. The strain fields also provide us with fundamental variables for construction of constitutive responses of the continuum. Once the strain fields are fixed, we look for the necessary and sufficient (compatibility) conditions for the existence of a local isometric embedding of the surface. Finally, we discuss how various defect densities become sources of strain incompatibility precluding the existence of the local isometric embedding. This will then set the stage for posing complete boundaryvalueproblems for internal stress distribution and natural shapes of defective structured surfaces, as will be discussed subsequently.
4.1 Strain measures and strain compatibility
Let us assume that there exist a local isometric embedding of into . Let and . The first and second fundamental forms associated with this embedding (over a local patch) are therefore and , respectively. We consider the following sufficiently smooth fields, defined over , as descriptors of strain on the structured surface: (i) a symmetric tensor