Next Article in Journal
Holographic Projection of Electromagnetic Maxwell Theory
Next Article in Special Issue
On the Failure of Classic Elasticity in Predicting Elastic Wave Propagation in Gyroid Lattices for Very Long Wavelengths
Previous Article in Journal
Asymmetries in Processes of Electron–Positron Annihilation
Previous Article in Special Issue
Second Gradient Electromagnetostatics: Electric Point Charge, Electrostatic and Magnetostatic Dipoles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cosserat Model of Elastic Solids Reinforced by a Family of Curved and Twisted Fibers

Department of Mechanical Engineering, University of California, Berkeley, CA 94720, USA
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(7), 1133; https://doi.org/10.3390/sym12071133
Submission received: 3 June 2020 / Revised: 19 June 2020 / Accepted: 29 June 2020 / Published: 7 July 2020
(This article belongs to the Special Issue Recent Advances in the Study of Symmetry and Continuum Mechanics)

Abstract

:
A Cosserat theory for fiber-reinforced elastic solids developed in Steigmann (2012) is generalized to accommodate initial curvature and twist of the fibers. The basic variables of the theory are a conventional deformation field and a rotation field that describes the local fiber orientation. Constraints on these fields are introduced to model the materiality of the fibers with respect to the underlying matrix deformation. A variational argument delivers the relevant equilibrium equations and boundary conditions and furnishes the interpretation of the Lagrange multipliers associated with the constraints as shear tractions acting on the fiber cross sections. Finally, the theory of material symmetry for such solids is developed and applied to the classification of some explicit constitutive functions.

1. Introduction

In the present, work we generalize a theory for fiber-reinforced elastic solids proposed in [1,2] that accounts for the intrinsic flexural and torsional elasticities of the fibers, regarded as continuously distributed spatial rods of the Kirchhoff type in which the kinematics are based on a position field and an orthonormal triad field [3,4,5]. This model is a special case of the Cosserat theory of nonlinear elasticity [6,7,8,9,10,11,12]. We extend this theory to accommodate initially curved and twisted fibers and develop an associated framework for the characterization of material symmetry.
Industrial applications of the mechanics of composite materials reinforced by curvilinear fibers are thoroughly treated in [13]. Further applications to bioelasticity are described in [14]. In this literature, the fibers confer anisotropy to the composite but their instrinsic flexural and torsional elasticities are not taken into account. However, the latter can be expected to play a significant role in local fiber buckling and kink-band failure due to the length scale inherent in the flexural and torsional stiffnesses of the fibers [15]. These stiffnesses are also significant at larger length scales if the fibers are sufficiently stiff relative to the underlying matrix material.
To aid in the interpretation of the theory to be developed, in Section 2 we review the basic elements of Kirchhoff’s theory for single rods. This is followed, in Section 3, by a brief outline of nonlinear Cosserat elasticity, specialized to model the effects of a single family of embedded fibers interacting with an elastic matrix. The resulting model is similar in structure to the Kirchhoff theory, with the effects of fiber-matrix interaction manifesting themselves as distributed forces and couples transmitted to the fibers by the matrix material in which they are embedded. Section 4 is devoted to a development of the associated theory of material symmetry, based on an extension to Cosserat elasticity [16] of Noll’s concept [17] for simple materials. This is used in Section 5 to discuss some particular constitutive functions for fiber-reinforced solids.
We use standard notation such as A t , A 1 , S k w A , det A and t r A . These are respectively the transpose, the inverse, the skew part, the deteminant and the trace of a tensor A , regarded as a linear transformation from a three-dimensional vector space to itself. The axial vector a x ( S k w A ) of S k w A is defined by a x ( S k w A ) × v = ( S k w A ) v for any vector v . The tensor product of three-vectors is indicated by interposing the symbol , and the Euclidean inner product of tensors A , B is denoted and defined by A · B = t r ( AB t ) ; the induced norm is A = A · A . The symbol · is also used to denote the usual Euclidean norm of three-vectors. Latin and Greek indices take values in { 1 , 2 , 3 } and { 2 , 3 } respectively, and, when repeated, are summed over their ranges. Finally, bold subscripts are used to denote derivatives of scalar functions with respect to their vector or tensor arguments.

2. Kirchhoff Rods

Kirchhoff rods are modelled as spatial curves endowed with an elastic energy density that responds to flexure and twist. According to the derivation from conventional three-dimensional nonlinear elasticity given in [4], this theory also accommodates a small axial strain along the rod. We forego any discussion of the connection between Kirchhoff theory and three-dimensional elasticity and simply regard the rod as a directed curve [5] in which certain a priori constraints are imposed. An accessible discussion of Kirchhoff’s theory may be found in [3].

2.1. Kinematics

The basic kinematical variables in the theory are a deformation field r ( s ) , where s [ 0 , l ] and l is the length of the rod in a reference configuration, and a right-handed, orthonormal triad { d i ( s ) } in which d 1 = d , where d is the unit vector defined by
r ( s ) = λ d , and λ = r ( s ) ,
where λ is the stretch of the rod. Thus d is the unit tangent to the rod in a deformed configuration and d α ( α = 2 , 3 ) span its cross-sectional plane at arclength station s .
The central assumption in Kirchhoff’s theory is that each cross section deforms as a rigid disc. Thus there is a rotation field R ( s ) given by
R = d i D i ,
such that d i = RD i , where D i ( s ) are the values of d i ( s ) in the reference configuration.
The curvature and twist of the rod are computed from the derivatives d i ( s ) , where
d i = R D i + RD i .
Let { E i } be a fixed right-handed background frame. Then D i ( s ) = A ( s ) E i for some rotation field A , yielding
d i = Wd i = w × d i ,
where
W = R R t + RA A t R t
is a skew tensor and
w = a x W = κ i d i ,
with
κ i = 1 2 e i j k d k · d j .
Here, e i j k is the permutation symbol ( e 123 = + 1 , etc.), κ 1 is the twist of the rod and κ α are the curvatures.

2.2. Strain-Energy Function

The strain energy S stored in a rod of length l is assumed to be expressible as
S = 0 l U d s ,
where U , the energy per unit initial length, is a function of the list { R , R , r } , possibly depending explicitly on s .
We assume U to be Galilean invariant and hence that its values are invariant under { R , R , r } { QR , QR , Qr } , where Q is an arbitrary uniform rotation. Because U is defined pointwise, to derive a necessary condition we select Q = R s t and conclude that U is determined by the list { R t R , R t r } = { R t WR A A t , λ D } , where D = D 1 and
R t WR A A t = R t R = ( RD i · R D j ) D i D j ,
is a Galilean-invariant measure of the relative flexure and twist of the rod due to deformation. This stands in one-to-one relation to its axial vector
γ = γ i D i = a x ( R t R ) ,
with components
γ i = 1 2 e i j k RD k · R D j .
Accordingly, because D ( s ) is independent of the deformation, the list { R t R , R t r } is equivalent to the list { γ , λ } , and the strain energy may therefore be written in the form
U = w ( λ , γ ; s ) .
Using (7) and (11) with R D j = d j RD j yields
γ i = κ i κ i 0 ,
where
κ i 0 = 1 2 e i j k D k · D j
are the components of the initial curvature-twist vector
κ 0 = κ i 0 D i = a x ( A A t ) .
Thus,
γ = κ κ 0 ,
where
κ = R t w = κ i D i = a x ( R t WR ) .
For example, in the classical theory [4,18] of initially straight and untwisted rods ( κ 0 = 0 ) , the strain-energy function for an isotropic rod of circular cross section is
w ( λ , κ ; s ) = 1 2 A ( s ) ε 2 + 1 2 T ( s ) τ 2 + 1 2 F ( s ) κ α κ α ,
where ε = λ 1 is the extensional strain, τ = κ 1 is the twist, A is the extensional stiffness (Young’s modulus E times the cross-sectional area); F is the flexural stiffness (Young’s modulus times the 2nd moment of area I of the cross section); and T is the torsional stiffness (the shear modulus G times the polar moment J of the cross section).
The homogeneous quadratic dependence of the energy on the bending-twist strain may be understood in terms of the local length scale furnished by the diameter of a cross section. The curvature-twist vector, when non-dimensionalized by this local scale, is typically small in applications. For example, the minimum radius of curvature of a bent rod is typically much larger than its diameter. If the bending and twisting moments vanish when the rod is straight and untwisted, then the leading-order contribution of the curvature-twist vector to the strain energy is quadratic. In general the flexural and torsional stiffnesses in this expression may depend on fiber stretch, but in the small-extensional-strain regime they are approximated at leading order by functions of s alone.

2.3. Equilibrium Theory

We recall the variational derivation of the equilibrium equations of the Kirchhoff theory here to provide context for the discussion of Cosserat elasticity in Section 3 [18]. Equilibria are assumed to satisfy the virtual-power statement
S ˙ = P ,
where P is the virtual power of the loads—the explicit form of which is deduced below—and the superposed dot is used to identify a variational derivative. These are induced by the derivatives, with respect to ϵ , of the one-parameter deformation and rotation fields r ( s ; ϵ ) and R ( s ; ϵ ) respectively, where r ( s ) = r ( s ; 0 ) and R ( s ) = R ( s ; 0 ) are equilibrium fields. Thus,
U ˙ = w ˙ = w λ λ ˙ + m i γ ˙ i ,
where
w λ = w / λ and m i = w / γ i
are evaluated at ϵ = 0 .
From (1) we have that
λ ˙ d + ω × r = u ,
where u ( s ) = r ˙ is the virtual translational velocity and ω ( s ) = a x ( R ˙ R t ) is the virtual rotational velocity. That is, d ˙ i = R ˙ R t d i , which is equivalent to
d ˙ i = ω × d i .
From (9) and (16) it follows that
κ ˙ i = 1 2 e i j k ( d ˙ k · d j + d k · d ˙ j ) = 1 2 e i j k [ ω × d k · d j + d k · ( ω × d j + ω × d j ) ] ,
in which the terms involving ω cancel; the e δ identity 1 2 e i j k e m j k = δ i m (the Kronecker delta), combined with d j × d k = e m j k d m , results in
γ ˙ i = κ ˙ i = d i · ω .
Thus,
S ˙ = 0 l ( w λ d · u + m · ω ) d s ,
with
m = m i d i .
Further, (1) implies that
r · d α = 0 ; α = 2 , 3 .
To accommodate these constraints in the virtual-power statement, we relax them and introduce the extended energy
E = S + 0 l f α r · d α d s ,
where f α ( s ) are Lagrange multipliers. The extended variational problem is
E ˙ = P ,
where
E ˙ = 0 l [ ( w λ d + f α d α ) · u + m · ω + f α d α × r · ω + f ˙ α r · d α ] d s
The variations f ˙ α simply return the constraints (28), and an integration by parts gives
E ˙ = ( f · u + m · ω ) 0 l 0 l [ u · f + ω · ( m f × r ) ] d s ,
where
f = w λ d + f α d α .
This implies that the virtual power is expressible in the form
P = ( t · u + c · ω ) 0 l + 0 l ( u · g + ω · π ) d s ,
in which t and c represent forces and couples acting at the ends of the segment and g and π are force and couple distributions acting in the interior.
By the Fundamental Lemma, the Euler equations holding at points in the interior of the rod are
f + g = 0 and m + π = f × r ,
and the endpoint conditions are
f = t and m = c ,
provided that neither position nor orientation is assigned at the endpoints. These are the equilibrium conditions of classical rod theory in which f and m respectively are the cross-sectional force and moment transmitted by the segment ( s , l ] on the part [ 0 , s ] . We observe, from (33) and (36)1, that the Lagrange multipliers f α play the role of constitutively indeterminate transverse shear forces acting on a fiber cross section.
For the strain-energy function (18) we have m 1 d 1 = T τ d and m α d α = F κ α d α . To reduce the second expression we use (7), together with d · d μ = d μ · d , to derive κ α = e α 1 μ d μ · d . From d · d = 0 it follows that d = ( d α · d ) d α and d × d = ( d α · d ) d × d α = ( e β 1 α d α · d ) d β ; thus κ β d β = d × d and (21) yields [3]
m = T τ d + F d × d .

3. Cosserat Elasticity of Fiber-Reinforced Materials

Cosserat elasticity theory emerges as the natural setting for elastic solids with embedded fibers—modelled as continuously distributed Kirchhoff rods—that support bending and twisting moments. To motivate our kinematical hypotheses, we suppose the fibers and matrix to be perfectly bonded and assume that both may be modelled at the microscale as conventional elastic solids. The interface between the matrix and fiber is then convected by the deformation as a material surface, and Hadamard’s compatibility condition requires that F + F = a N for some vector a , where N is a unit normal to the interface and F ± are the values of the deformation gradients in the fiber and matrix at the interface. In particular, if D ( = D 1 ) is the unit tangent to the centerline of an untapered fiber, then N S p a n { D α } , where D α ; α = 2 , 3 , are orthonormal unit vectors in the fiber cross section. It follows that
F + D = F D , but F + D α F D α ,
and hence that the deformation gradients in the fiber and matrix may be unequal. This stands in contrast to a model proposed in [19] on the basis of a single deformation field.
If a fiber is sufficiently stiff relative to the matrix, then its deformation gradient is approximated by a rotation field R . In Dill’s interpretation of the Kirchhoff theory [4] this is accompanied by a small axial strain. Thus we interpret (38)1 in the form
FD = λ d , where d = RD and λ = FD ,
where λ ( = FD ) is the fiber stretch and F is the matrix deformation gradient. The fields F and R are otherwise independent in accordance with (38)2. These in turn furnish R t FD = λ D and hence two constraints
D α · R t FD = 0 ; α = 2 , 3 ,
analogous to (28), involving the fiber rotation and matrix deformation.
Equation (39) implies that the fibers are convected as material curves relative to the matrix. The cross-sectional vectors D α are embedded in the fiber but not in the matrix, and so their images d α in the current configuration are free to shear relative to the matrix while remaining mutually orthogonal and perpendicular to d .

3.1. Kinematical and Constitutive Variables in Cosserat Elasticity

The basic kinematical variables of a Cosserat continuum are a rotation field R ( X ) and a deformation field χ ( X ) . Naturally these may depend on time, but such dependence is not relevant to our development and is not made explicit.
The constitutive response of an elastic Cosserat continuum is embodied in a strain-energy density U ( F , R , R ; X ) , per unit reference volume, where F = χ is the usual deformation gradient and R is the rotation gradient. In Cartesian index notation, these are
F = F i A e i E A , R = R i A e i E A and R = R i A , B e i E A E B
with
F i A = χ i , A ,
where ( · ) , A = ( · ) / X A and where { e i } and { E A } are fixed orthonormal bases associated with the Cartesian coordinates x i and X A , with x i = χ i ( X A ) .
We again suppose the strain energy to be Galilean-invariant and thus require
U ( F , R , R ; X ) = U ( QF , QR , Q R ; X ) ,
where Q is an arbitrary spatially uniform rotation with ( Q R ) i A B = ( Q i j R j A ) , B = Q i j R j A , B . The restriction
U ( F , R , R ; X ) = W ( E , Γ ; X ) ,
with [11,16]
E = R t F = E A B E A E B ; E A B = R i A F i B ,
Γ = Γ D C E D E C ; Γ D C = 1 2 e B A D R i A R i B , C ,
where W is the reduced strain-energy function and e A B C is the permutation symbol, is both necessary and sufficient for Galilean invariance. Sufficiency is obvious, whereas necessity follows by choosing Q = R X t , where X is the material point in question, and making use of the fact that, for each fixed C { 1 , 2 , 3 } , the matrix R i A R i B , C is skew. This follows by differentiating R i A R i B = δ A B (the Kronecker delta). The axial vectors Γ C associated with this skew matrix have components
Γ D ( C ) = 1 2 e B A D R i A R i B , C ,
yielding [11]
Γ = Γ C E C ,
and so Γ —the second order wryness tensor—stands in one-to-one relation to the third order Cosserat strain measure R t R .

3.2. Virtual Power and Equilibrium

As in Section 2, we define equilibria to be states that satisfy the virtual-power statement
E ˙ = P ,
where P is the virtual power of the loads acting on the body,
E = S + κ Λ α D α · ED d v ,
is the extended energy, Λ α are Lagrange multipliers accompanying the constraints (40),
S = κ U d v
is the total strain energy, and, as before, superposed dots identify variational derivatives. Thus,
U ˙ = W ˙ = σ · E ˙ + μ · Γ ˙ ,
where
σ = W E and μ = W Γ
are evaluated at equilibrium. Further,
( D α · ED ) · = D α D · E ˙
so that
E ˙ = κ [ ( σ + Λ D ) · E ˙ + μ · Γ ˙ + Λ ˙ α D α · ED ] d v ,
where
Λ = Λ α D α .
It follows from (45) that
E ˙ = R t ( u Ω F ) , where u = χ ˙ and Ω = R ˙ R t .
Then,
( σ + Λ D ) · E ˙ = R ( σ + Λ D ) · u Ω · S k w [ R ( σ + Λ D ) F t ] .
If α is a skew tensor, then Ω · α = 2 ω · a , where ω = a x Ω and a = a x α . Further, for any tensor A we have RAF t = RAE t R t and S k w ( RAE t R t ) = R S k w ( AE t ) R t , and therefore
( σ + Λ D ) · E ˙ = R ( σ + Λ D ) · u 2 a x { R S k w [ ( σ + Λ D ) E t ] R t } · ω .
The reduction
Γ ˙ = R t ω ,
which is somewhat more involved, is detailed in Appendix A. Accordingly,
μ · Γ ˙ = R μ · ω ,
and application of the divergence theorem to (55), with (59) and (61), gives
E ˙ = κ [ ( R σ + λ D ) ν · u + ( R μ ) ν · ω ] d a + κ Λ ˙ α D α · ED d v κ { u · D i v ( R σ + λ D ) + ω · [ D i v ( R μ ) + 2 a x ( R S k w [ ( σ + Λ D ) E t ] R t ) ] } d v ,
where ν is the exterior unit normal to the (piecewise smooth) surface κ , and
λ = R Λ = Λ α d α .
The virtual power is thus of the form
P = κ ( t · u + c · ω ) d a + κ ( g · u + π · ω ) d v ,
where t and c are densities of force and couple acting on κ , and g and π are densities of force and couple acting in κ .
The fundamental lemma delivers the constraints (40) together with the differential equations
g = D i v ( R σ + λ D ) and π = D i v ( R μ ) 2 a x { R S k w [ ( σ + Λ D ) E t ] R t } in κ ,
and the natural boundary conditions
t = ( R σ + λ D ) ν on κ t and c = ( R μ ) ν on κ c ,
where κ t is a part of κ where position is not assigned and κ c is a part where rotation is not assigned. We assume position to be assigned on κ κ t ( u = 0 ) , and rotation to be assigned on κ κ c ( ω = 0 ) .

3.3. Fiber-Matrix Interaction

Pursuant to the discussion at the start of this section, we assume that Cosserat elasticity is conferred by the mechanical interaction between an elastomeric matrix and a single family of embedded fibers. The relative curvature-twist vector γ = a x ( R t R ) of a fiber initially oriented along a unit-vector field D ( X ) , where ( · ) is the directional derivative along D , is (cf. (11))
γ = γ i D i with γ i = 1 2 e i j k D k · R t R D j .
Thus, with R i A = R i A , B D B we derive (cf. (46))
R t R = R i C R i A , B D B E C E A = e A C D Γ D B D B E C E A ,
and conclude that γ is determined by Γ via Γ D . Here the director fields D i ( X ) form a positively oriented orthonormal triad with D 1 = D , and the Cosserat rotation field is given simply by
R = d i D i ,
as in (2), where d i = RD i are the images of directors in the deformed body, with d = d 1 the field of unit tangents to the deformed fibers. Then, as in Section 2, we may express the strain energy in the form
W ( E , Γ ; X ) = w ( E , γ ; X ) ,
where w is now the strain energy per unit reference volume. Thus,
σ = w E .
To obtain the couple stress μ we use (25) in the form
γ ˙ i = RD i · ( ω ) D = D i · ( R t ω ) D = D i D · Γ ˙ .
Variation of the energy at fixed E then gives
μ · Γ ˙ = W ˙ = w ˙ = w γ · γ ˙ = M D · Γ ˙ ,
where
M = w γ = m i D i with m i = w / γ i ,
yielding
μ = M D .
The equilibrium equations for this model follow simply from (65) and (66). To facilitate comparison with Kirchhoff’s theory we use
D i v ( λ D ) = λ + ( D i v D ) λ , where λ = ( λ ) D
is the fiber-derivative of λ , the directional derivative along the fiber passing through the point with reference position X . Further,
2 a x { R S k w [ ( Λ D ) E t ] R t } = 2 a x { S k w [ R ( Λ D ) E t R t ] } = 2 a x [ S k w ( λ FD ) ] = a x ( λ χ χ λ ) = χ × λ ,
where
χ = ( χ ) D
is the fiber derivative of the deformation. Lastly,
D i v ( R μ ) = m + ( D i v D ) m , where m = RM = m i d i ,
with fiber derivative m = ( m ) D . The equilibrium conditions (65) holding in κ thus specialize to
λ + ( D i v D ) λ + D i v ( R σ ) + g = 0
and
m + χ × λ + ( D i v D ) m + 2 a x [ R S k w ( σ E t ) R t ] + π = 0 ,
with
t = ( R σ ) ν + ( D · ν ) λ and c = ( D · ν ) m
holding on κ t and κ c , respectively. From (35) and (36) these conditions yield the interpretation of λ and m respectively as shear force and moment densities acting on fiber ’cross sections’, i.e., on surfaces that intersect fibers orthogonally ( D · ν = ± 1 ) . We observe that no solution exists if a non-zero couple is specified on a part of κ c containing D as a tangent vector ( D · ν = 0 ) . Additionally, comparison of (80) with (35) furnishes the interpretation of D i v ( R σ ) as a distributed density of force transmitted to a fiber by the matrix in which it is embedded, whereas comparison of (81) and (35) implies that 2 a x [ R S k w ( σ E t ) R t ] is a density of distributed couple transmitted by the matrix to the fiber. The derivation of the system (80)–(82) simplifies and and generalizes that of a similar system for initially straight fibers presented in [1,2].
The dependence of the strain-energy function on γ (or Γ ) introduces a natural length scale, L say, into the constitutive theory which is on the order of that of the microstructure and hence of the diameter of a fiber cross section or the spacing between adjacent fibers. Using the larger of these to define the dimensionless curvature-twist vector L γ , supposing that L γ 1 in typical applications and assuming that the fibers transmit no moments when γ vanishes, we find that w is given to leading order by
w ( E , γ ; X ) = ϖ ( E ; X ) + 1 2 γ · K ( E ; X ) γ ,
where ϖ ( E ; X ) = w ( E , 0 ; X ) and K ( E ; X ) = w γ γ γ = 0 . For small E we have K ( E ; X ) = K ( 0 ; X ) + O ( E ) , provided that K ( · ; X ) is differentiable. Then the energy is approximated, as in (18), by the decoupled energy
w ( E , γ ; X ) = ϖ ( E ; X ) + φ ( γ ; X ) ,
for some homogeneous quadratic function φ ( · ; X ) .

4. Material Symmetry

In this section we develop the theory of material symmetry for elastic Cosserat materials subjected to (40). A comprehensive study of material symmetry in the setting of Cosserat elasticity, extending Noll’s concept [17] for simple materials, is given in [16]. The concept of material symmetry in rods is discussed in [20,21,22]. As a preliminary step we first describe the manner in which the constitutive function for the strain energy may be computed for any choice of reference configuration when that pertaining to any particular choice is given.

4.1. Change of Reference Configuration

Let κ and μ be two reference configurations, and let Y = Π ( X ) be a diffeomorphism mapping points in κ to points in μ . The deformation gradients relative to κ and μ , denoted by F κ and F μ respectively, are related by
F κ = F μ H , where H = Π .
We restrict attention to transformations Π with det H = 1 , for reasons that are well known in conventional elasticity [8,16], and also impose Π ( X 0 ) = X 0 . The specification of this pivot point removes an inessential translational degree of freedom. Henceforth, we are concerned with the properties of the map Π ( X ) in a neighborhood N κ ( X 0 ) κ of the pivot point. This neighborhood is mapped by Π to the neighborhood N μ ( X 0 ) of the pivot.
The Cosserat rotation R κ relative to κ is such that d i = R κ D i . In the same way there is a rotation R μ such that d i = R μ G i , where { G i ( Y ) } is the positively-oriented orthonormal director field defined in μ . Thus,
R κ = R μ L ,
where
L = G i D i
is the rotation field that maps the directors in κ to their images in μ . We have d = R κ D = R μ G , where G ( = G 1 ) is the unit-tangent field to fibers in μ , so that G = LD . To ensure that D remains a material vector relative to the matrix (cf. (39)) under the change of reference, it is necessary to impose
HD = HD LD .
Following the characterization of solids in [8], we assume the existence of an undistorted reference and suppose κ to be one of these. Thus we confine attention to proper-orthogonal H . Further, we remove an inessential orientational degree of freedom in the local change of reference by requiring that it preserve the pivotal axis D at the point X 0 ; thus, HD = 1 and
D = HD = LD .
With G α = LD α and G = HD , this in turn implies that
G α · G = LD α · HD = LD α · LD = D α · D ,
and hence that the constraints G α · G = 0 are automatically satisfied (cf. (40)). Accordingly, L , H S , where
S = { Q Q is a rotation with axis D } .
For the model discussed in Section 3, the strain energy depends on the Cosserat rotation and its gradient via γ = γ i D i , where γ i is given by (67) in which the prime refers to the fiber derivative in the reference configuration κ . In particular, for any function f we have f = κ f · D = H t ( μ f ) · D = μ f · HD , where the subscripts κ and μ identify gradients with respect to X κ and Y μ , respectively. In view of (89) we have f = μ f · D at the pivot X 0 , implying that the fiber derivative is invariant under transformations of the reference configuration that preserve the fiber axis. Accordingly, it is immediately apparent that the κ i , defined by (7) in which the prime is again a fiber derivative, are also invariant.
Alternatively, we may use (67) to derive
( κ i ) κ ( κ i 0 ) κ = ( γ i ) κ = 1 2 e i j k D k · R κ t R κ D j ,
in which
R κ t R κ = L t ( R μ t R μ + L L t ) L .
Thus,
( κ i ) κ ( κ i 0 ) κ = ( κ i ) μ ( κ i 0 ) μ + Δ i ,
where
( κ i ) μ ( κ i 0 ) μ = ( γ i ) μ = 1 2 e i j k G k · R μ t R μ G j
and
Δ i = 1 2 e i j k G k · L L t G j .
Noting that
( κ i 0 ) μ = 1 2 e i j k G k · G j = 1 2 e i j k LD k · ( L D j + LD j ) = 1 2 e i j k G k · L ( L t G j ) + ( κ i 0 ) κ ,
and hence that
( κ i 0 ) μ = ( κ i 0 ) κ + Δ i ,
we conclude, in accordance with (94), that
( κ i ) κ = ( κ i ) μ ,
as claimed.
From (67), the curvature-twist strains γ κ and γ μ relative to κ and μ are related by
γ κ = L t ( γ μ + Δ ) ,
where
Δ = Δ i G i = a x ( L L t ) ,
whereas the Cosserat strains E κ and E μ are related by
E κ = L t E μ H .
Because the state of the material is not affected by the choice of reference, we require
w μ ( E μ , γ μ ; X 0 ) = w κ ( E κ , γ κ ; X 0 ) = w κ ( L t E μ H , L t ( γ μ + Δ ) ; X 0 ) .

4.2. Material Symmetry Transformations

According to Noll’s theory [17], N κ ( X 0 ) and N μ ( X 0 ) are related by material symmetry if their responses to a given experiment are identical at the pivot point. In the present context an experiment consists of a deformation function χ ( · ) and rotation function R ( · ) . Thus the experiment acts on N κ ( X 0 ) to produce the fields { χ ( X ) , R ( X ) } for X N κ ( X 0 ) , and on N μ ( X 0 ) to produce { χ ( Y ) , R ( Y ) } for Y N μ ( X 0 ) . Accordingly, N κ ( X 0 ) and N μ ( X 0 ) are both subjected to the same pair { F , R } , and hence the same strain E , at X , Y = X 0 . Moreover, μ R = ( κ R ) H and with (89) we infer that the fiber derivatives ( R ) κ and ( R ) μ , relative to N κ ( X 0 ) and N μ ( X 0 ) respectively, also coincide at X 0 . This in turn implies that both neighborhoods experience the same bend-twist strain γ at X 0 , a fact that is most easily appreciated by using a single background frame { E i } to evaluate
γ = γ i E i , with γ i = 1 2 e i j k E k · R t R E j ,
in both neighborhoods.
Accordingly, material symmetry is tantamount to
w κ ( E , γ ; X 0 ) = w μ ( E , γ ; X 0 ) ,
which, when combined with (100), (102) and (103), yields the restriction
w κ ( E , γ ; X 0 ) = w κ ( L t EH , L t ( γ + Δ ) ; X 0 )
on the single response function w κ , where the rotations H and L are connected by (89), but otherwise independent, and of course Δ is evaluated at X 0 .
For the decoupled energy (84) considered hereafter, this is equivalent to
ϖ ( E ; X 0 ) + φ ( γ ; X 0 ) = ϖ ( L t EH ; X 0 ) + φ ( L t ( γ + Δ ) ; X 0 )

5. Examples

We close with some examples of constitutive functions that conform to (107).

5.1. Matrix Energy

Consider the list [1]
I = { I 1 , , I 9 } ,
of functionally independent scalar-valued functions of E , with
I 1 = t r ( E t E ) , I 2 = t r [ ( E t E ) 2 ] , I 3 = det E , I 4 = D · ED , I 5 = D · ( E t E ) D , I 6 = D · ( E E t ) D , I 7 = D · E * D , I 8 = D · ( E t E ) 2 D , I 9 = D · ( EE t ) 2 D ,
where E * = ( det E ) E t is the cofactor of E . We note that det E = det F , E t E = C and EE t = R t BR , where C = F t F and B = FF t , respectively, are the right and left Cauchy–Green deformation tensors.
Straightforward calculations show, remarkably, that each member I k of this list satisfies
I k ( L t EH ) = I k ( E )
for any—hence every— L , H S . Particular matrix energies may thus be obtained by taking
ϖ ( E ; X ) = M ( I 1 , , I 9 ; X ) ,
for some function M . This satisfies
ϖ ( E ; X ) = ϖ ( L t EH ; X )
for all distinct L , H S . For example, (111) furnishes energies for transversely hemitropic matrix materials for which H is an arbitrary element of S, without any restrictions on the independent fiber rotations L S . However, we have not shown that I is a function basis for any particular kind of symmetry.
The stress σ associated with this energy is given by
σ = ϖ E = j M j ( I j ) E , where M j = M / I j
and [1]
( I 1 ) E = 2 E , ( I 2 ) E = 4 EC , ( I 3 ) E = E * , ( I 4 ) E = D D , ( I 5 ) E = 2 E ( D D ) , ( I 6 ) E = 2 ( D D ) E , ( I 7 ) E = I 7 E t I 3 E t ( D D ) E t , ( I 8 ) E = 2 E [ ( D D ) C + C ( D D ) ] , ( I 9 ) E = 2 [ ( D D ) EC + EE t ( D D ) E ] .
We observe that this model yields an asymmetric σ E t , and thus makes provision for distributed couples to be transmitted to the fibers by the matrix (cf. (81)), if the energy involves I 4 , I 6 , I 7 or I 9 .

5.2. Fiber Symmetry

For matrix energies satisfying (112), the restriction (107) reduces to
φ ( γ ; X ) = φ ( L t ( γ + Δ ) ; X )
for some L and Δ = a x ( L L t ) such that L ( X ) D ( X ) = D ( X ) (cf. (89)). The fiber derivative L thus satisfies
L D + LD = D ,
and with D = L t D this may be cast in the form
Δ × D + ( L I ) D = 0 .
We then have
Δ = Δ D + D × ( Δ × D )
in which Δ ( = Δ · D ) is unrestricted, whereas
D × ( Δ × D ) = [ ( L I ) D ] × D .
Consider the particular fiber with arclength parametrization X ( s ) passing through the point X 0 at s = s 0 , i.e., X ( s 0 ) = X 0 . If the fiber is curved at X 0 then D = K N there, where K = D > 0 is the principal curvature of the fiber and N is the unique principal normal. In this case the local Frenet triad { D , N , B } , where B = D × N is the binormal, is well defined. The Rodrigues representation formula for rotations [23] thus yields
L = D D + cos θ ( N N + B B ) + sin θ ( B N N B )
for some angle θ , which may be used with (118) and (119) to construct
Δ = Δ D K ( cos θ 1 ) B + K sin θ N ,
and with some effort we derive
L t ( γ + Δ ) = ( γ 1 + Δ ) D + v ( θ ) ,
where
v ( θ ) = cos θ 1 γ + sin θ 1 γ × D + K ( cos θ 1 ) B + K sin θ N
in which
1 = I D D
is the projection onto the fiber cross section. The symmetry condition (115) thus becomes
φ ( γ ; X ) = φ ( ( γ 1 + Δ ) D + v ( θ ) ; X ) .
If the fiber is initially straight, then K vanishes and Δ = Δ D , again with Δ arbitrary. The Frenet triad is not well defined in this case and so we use (120) with { N , B } replaced by { D α } , finding that (115) reduces to
φ ( γ ; X ) = φ ( ( γ 1 + Δ ) D + γ 2 i ( θ ) + γ 3 j ( θ ) ; X ) ,
where
i ( θ ) = cos θ D 2 sin θ D 3 and j ( θ ) = sin θ D 2 + cos θ D 3
in which θ differs from the angle appearing in (120).
For curved or straight fibers we differentiate (125) or (126) with respect to Δ , concluding that
φ γ · D = 0 ,
and hence that the fiber is insensitive to the twist strain. In view of the fact that the energies associated with twist and bending of a typical fiber are comparable in magnitude (cf. (18)), we regard such a circumstance as unrealistic and therefore require that Δ = 0 . Then, Δ = D × ( Δ × D ) and (128) is not applicable. In the case of initially straight fibers this implies that L = 0 , i.e., that L cannot vary along the fibers, although L need not vanish. A similar observation was made in [16] in the context of conventional Cosserat elasticity.

5.2.1. Transversely Hemitropic Fibers

Transversely hemitropic fibers are defined to be those for which (125) or (126) hold for arbitrary θ . For curved fibers we differentiate the first of these with respect to θ , obtaining
φ γ · v ( θ ) = 0 ,
where
v ( θ ) = cos θ ( K B + 1 γ ) × D sin θ ( K B + 1 γ )
is an arbitrary vector in the fiber cross section. Accordingly 1 ( φ γ ) vanishes and the bending moments m α vanish. We conclude that a curved fiber can be transversely hemitropic only if it has no flexural elasticity. It is therefore effectively a string with a possible twisting elasticity. Exceptionally, (129) imposes no restrictions on φ if 1 γ is fixed at the value K B .
For initially straight fibers, we differentiate (126) with respect to θ and evaluate the result at θ = 0 , obtaining
γ 2 ϕ / γ 3 = γ 3 ϕ / γ 2 ,
where ϕ ( γ 1 , γ 2 , γ 3 ) = φ ( γ i D i ; X ) . This implies that ϕ ( γ 1 , γ 2 , γ 3 ) = ψ ( γ 1 , γ 2 2 + γ 3 2 ) for some function ψ . Moreover, the referential version of (4), D i = κ 0 × D i , implies, for i = 1 , that K 2 = ( κ 2 0 ) 2 + ( κ 3 0 ) 2 and hence that κ α 0 = 0 if the fibers are initially straight. Thus,
φ ( γ ; X ) = ϝ ( J 1 , J 2 ; X )
for some function ϝ , where
J 1 = γ · D = γ 1 and J 2 = 1 κ = κ 2 2 + κ 3 2 .
Conversely, (132) satisfies the operative version of the restriction (115) in the present circumstances, namely
φ ( γ ; X ) = φ ( L t γ ; X ) ,
for all L ( X ) S , and thus furnishes the general representation for initially straight, transversely hemitropic fibers with bending and twisting resistance.
Homogeneous quadratic energies of this kind are of the form
φ ( γ ; X ) = 1 2 A ( X ) J 1 2 + 1 2 B ( X ) J 2 2 ,
which may be compared to the classical energy (18) for rods of circular cross section. The product J 1 J 2 is excluded because the replacement γ t γ transforms J 1 J 2 to t t J 1 J 2 . Accordingly, J 1 J 2 is not a homogeneous quadratic function of γ .
To derive the response function M (cf. (74)) we require the gradients
( J 1 ) γ = D and ( J 2 ) γ = J 2 1 1 κ .
Equations (74) and (135) then deliver
M = φ γ = A γ 1 D + B ( κ α D α )
and (cf. (79))
m = A γ 1 d + B d × d ,
where d = ( d ) D .

5.2.2. Transversely Orthotropic Fibers

Transversely orthotropic fibers are defined to be those for which (125) or (126) are satisfied with θ = 0 and θ = π . This models fibers having rectangular cross sections on the microscale. The first alternative corresponds to L ( X ) I and Δ = 0 , for which (115) reduces to an identity. For initially curved fibers, the second alternative, θ = π , reduces (125) to the severe restriction
φ ( γ 1 D + 1 γ ; X ) = φ ( γ 1 D 1 γ 2 K B ; X )
on the bending response. For an arbitrarily shaped fiber this can be satisfied if the energy depends solely on the twist strain. Thus, as a practical matter, transverse orthotropy is meaningful for initially straight rods, for which (126) reduces to
φ ( γ 1 D + 1 κ ; X ) = φ ( γ 1 D 1 κ ; X ) ,
or simply ϕ ( γ 1 , κ 2 , κ 3 ) = ϕ ( γ 1 , κ 2 , κ 3 ) . Then,
φ ( γ ; X ) = ψ ( γ 1 , δ 2 , δ 3 ; X ) , where δ α = 1 2 κ α 2 ,
and
m = ( ψ / γ 1 ) d + ( ψ / δ 2 ) κ 2 d 2 + ( ψ / δ 3 ) κ 3 d 3 .
The foregoing considerations lead to the conclusion that initially curved fibers exhibiting elastic resistance to twisting and bending can be expected to have only trivial symmetry, i.e., L ( X ) I . This conclusion, perhaps unexpected, is due to the fact that in the present theory the initial shape of a fiber contributes to its material properties through the presence of κ 0 in the strain-energy function. To illustrate the point, imagine a fiber in the shape of a circular arc. If the arc is rotated about its tangent at point X by π , say, then it becomes the mirror image of the original with respect to the plane containing the tangent and lying perpendicular to that containing the fiber. We would not expect both local configurations of the fiber—the one before the rotation and the one after—to exhibit identical response to the same (arbitrary) experiment, unless the fiber is a string with negligible flexural stiffness.
In [2], several boundary value problems are solved explicitly in the context of a simpler version of the present theory in which the fibers are initially straight, untwisted and parallel. These pertain to finite torsion of a cylinder and to finite bending and transverse shearing of a block.

Author Contributions

Both authors contributed equally to the formulation and to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by US NSF grant CMMI-1931064.

Acknowledgments

MS gratefully acknowledges the support of the P.M. Naghdi Fellowship, administered by the Department of Mechanical Engineering at UC Berkeley.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To verify (60) we combine (46) and (57), obtaining
Γ ˙ D C = 1 2 e B A D ( R ˙ i A R i B , C + R i A R ˙ i B , C ) = 1 2 e B A D [ Ω i m R m A R i B , C + R m A ( Ω m j , C R j B + Ω m j R j B , C ) ] = 1 2 e B A D R m A [ Ω m i , C R i B + ( Ω m i + Ω i m ) R i B , C ] = 1 2 e B A D R i B R m A Ω m i , C .
On the other hand, because det R = 1 we have
e B A D R i B R m A R j D = e i m j ,
and hence
e B A C R i B R m A = e i m j R j C ,
where we have used R j D R j C = δ D C .
We obtain
Γ ˙ D C = 1 2 e B A D R i B R m A Ω m i , C = R j D ω j , C ,
where
ω j = 1 2 e i m j Ω m i
are the components of ω = a x Ω , and thus confirm that Γ ˙ = R t ω .

References

  1. Steigmann, D.J. Theory of elastic solids reinforced with fibers resistant to extension, flexure and twist. Int. J. Non-Linear Mech. 2012, 47, 734–742. [Google Scholar] [CrossRef] [Green Version]
  2. Steigmann, D.J. Effects of fiber bending and twisting resistance on the mechanics of fiber-reinforced elastomers. In CISM Course: Nonlinear Mechanics of Soft Fibrous Tissues; Dorfmann, L., Ogden, R.W., Eds.; Springer: Wien, Austria; New York, NY, USA, 2015; Volume 559, pp. 269–305. [Google Scholar]
  3. Landau, L.D.; Lifshitz, E.M. Theory of Elasticity, 3rd ed.; Pergamon: Oxford, UK, 1986. [Google Scholar]
  4. Dill, E.H. Kirchhoff’s theory of rods. Arch. Hist. Exact Sci. 1992, 44, 1–23. [Google Scholar] [CrossRef]
  5. Antman, S.S. Nonlinear Problems of Elasticity; Springer: Berlin, Germany, 2005. [Google Scholar]
  6. Cosserat, E.; Cosserat, F. Théorie des Corps Déformables; Hermann: Paris, France, 1909. [Google Scholar]
  7. Toupin, R.A. Theories of elasticity with couple stress. Arch. Ration. Mech. Anal. 1964, 17, 85–112. [Google Scholar] [CrossRef]
  8. Truesdell, C.; Noll, W. The Non-Linear Field Theories of Mechanics. In Handbuch der Physik, Vol. III/3; Flügge, S., Ed.; Springer: Berlin, Germany, 1965. [Google Scholar]
  9. Reissner, E. Note on the equations of finite-strain force and moment stress elasticity. Stud. Appl. Math. 1975, 54, 1–8. [Google Scholar]
  10. Reissner, E. A further note on finite-strain force and moment stress elasticity. Z. Angew. Math. Phys. 1987, 38, 665–673. [Google Scholar] [CrossRef]
  11. Pietraszkiewicz, W.; Eremeyev, V.A. On natural strain measures of the nonlinear micropolar continuum. Int. J. Solids Struct. 2009, 46, 774–787. [Google Scholar] [CrossRef] [Green Version]
  12. Neff, P. Existence of minimizers for a finite-strain micro-morphic elastic solid. Proc. Roy. Soc. Edinb. A 2006, 136, 997–1012. [Google Scholar] [CrossRef] [Green Version]
  13. Akbarov, S.D.; Guz, A.N. Mechanics of Curved Composites; Kluwer: Dordrecht, The Netherlands, 2000. [Google Scholar]
  14. Dorfmann, L.; Ogden, R.W. CISM Course: Nonlinear Mechanics of Soft Fibrous Tissues; Dorfmann, L., Ogden, R.W., Eds.; Springer: Wien, Austria; New York, NY, USA, 2015; Volume 559. [Google Scholar]
  15. Basu, S.; Waas, A.M.; Ambur, D.R. Compressive failure of fiber composites under multi-axial loading. J. Mech. Phys. Solids 2006, 54, 611–634. [Google Scholar] [CrossRef]
  16. Eremeyev, V.A.; Pietraszkiewicz, W. Material symmetry group of the non-linear polar-elastic continuum. Int. J. Solids Struct. 2012, 49, 1993–2005. [Google Scholar] [CrossRef] [Green Version]
  17. Noll, W. A mathematical theory of the mechanical behavior of continuous media. Arch. Ration. Mech. Anal. 1958, 2, 197–226. [Google Scholar] [CrossRef]
  18. Steigmann, D.J. The variational structure of a nonlinear theory for spatial lattices. Meccanica 1996, 31, 441–455. [Google Scholar] [CrossRef] [Green Version]
  19. Spencer, A.J.M.; Soldatos, K.P. Finite deformations of fibre-reinforced elastic solids with fibre bending stiffness. Int. J. Non-Linear Mech. 2007, 42, 355–368. [Google Scholar] [CrossRef] [Green Version]
  20. Healey, T.J. Material symmetry and chirality in nonlinearly elastic rods. Math. Mech. Solids 2002, 7, 405–420. [Google Scholar] [CrossRef]
  21. Luo, C.C.; O’Reilly, O.M. On the material symmetry of elastic rods. J. Elast. 2000, 60, 35–56. [Google Scholar] [CrossRef]
  22. Lauderdale, T.; O’Reilly, O.M. On the restrictions imposed by non-affine material symmetry groups for elastic rods: Application to helical substructures. Eur. J. Mech. A/Solids 2007, 26, 701–711. [Google Scholar] [CrossRef]
  23. Chadwick, P. Continuum Mechanics: Concise Theory and Problems; Dover: New York, NY, USA, 1976. [Google Scholar]

Share and Cite

MDPI and ACS Style

Shirani, M.; Steigmann, D.J. A Cosserat Model of Elastic Solids Reinforced by a Family of Curved and Twisted Fibers. Symmetry 2020, 12, 1133. https://doi.org/10.3390/sym12071133

AMA Style

Shirani M, Steigmann DJ. A Cosserat Model of Elastic Solids Reinforced by a Family of Curved and Twisted Fibers. Symmetry. 2020; 12(7):1133. https://doi.org/10.3390/sym12071133

Chicago/Turabian Style

Shirani, Milad, and David J. Steigmann. 2020. "A Cosserat Model of Elastic Solids Reinforced by a Family of Curved and Twisted Fibers" Symmetry 12, no. 7: 1133. https://doi.org/10.3390/sym12071133

APA Style

Shirani, M., & Steigmann, D. J. (2020). A Cosserat Model of Elastic Solids Reinforced by a Family of Curved and Twisted Fibers. Symmetry, 12(7), 1133. https://doi.org/10.3390/sym12071133

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop