Elsevier

Composite Structures

Volume 244, 15 July 2020, 112279
Composite Structures

Asymptotically exact theory of fiber-reinforced composite beams

https://doi.org/10.1016/j.compstruct.2020.112279Get rights and content

Abstract

An asymptotic analysis of the energy functional of a fiber-reinforced composite beam with a periodic microstructure in its cross section is performed. From this analysis the asymptotically exact energy as well as the 1-D beam theory of first order is derived. The effective stiffnesses of the beam are calculated in the most general case from the numerical solution of the cell and homogenized cross-sectional problems.

Introduction

Fibre reinforced composite (FRC) beams are widely used in civil, mechanical, and aerospace engineering due to their low weight, high strength, and good damping properties [10], [9], [19]. The exact treatment of FRC beams within the 3D elasticity is only possible in a few exceptional cases due to their complicated microstructure (see, e.g., [24], [26] and the references therein). For this reason, different approaches have been developed depending on the type of beams. If FRC beams are thick, no exact one-dimensional theory can be constructed, so only the numerical methods or approximate semi-analytical methods applied to three-dimensional elasticity theory make sense. However, if the FRC beams are thin, the reduction from the three-dimensional to the one-dimensional theory is possible. This dimension reduction can be made asymptotically exact in the limit when the thickness-to-length ratio of the beam goes to zero. The rigorous derivation of the asymptotically exact one-dimensional beam theory based on the variational-asymptotic method (VAM) developed by Berdichevsky [3] was first performed in [4]. His asymptotic analysis shows that the static (or dynamic) three-dimensional problem of the beam can be split into two problems: (i) the two-dimensional cross-sectional problem and (ii) the one-dimensional variational problem whose energy functional should be found from the solution of the cross-sectional problem. The latter has been solved both for anisotropic homogeneous beams and for inhomogeneous beams with the constant Poisson’s ratio [4]. In addition to these findings, Berdichevsky [4] has shown that the average energy as well as the extension and bending stiffnesses of FRC beams with piecewise constant Poisson’s ratio must be larger than those of FRC beams with constant Poisson’s ratio (see also [23]). However, to our knowledge, the question of how the corrections in energy and stiffnesses depend on the difference in Poisson’s ratio of fibers and matrix of FRC beams remains still an issue. It should be noted that VAM has been further developed in connection with the numerical analysis of cross-sectional problems for geometrically nonlinear composite beams by Hodges, Yu and colleagues in a series of papers [11], [28], [29], [30]. Note also that VAM has been used, among others, to derive the 2D theory of homogeneous piezoelectric shells [13], the 2D theory of purely elastic anisotropic and inhomogeneous shells [5], the 2D theory of purely elastic sandwich plates and shells [6], [7], the theory of smart beams [25], the theory of low and high frequency vibrations of laminate composite shells [17], [18], and more recently, the theory of smart sandwich and functionally graded shells [15], [16].

For FRC beams that have the periodic microstructure in the cross section, an additional small parameter appears in the cross-sectional problems: The ratio between the length of the periodic cell and the characteristic size of the cross section. In this case the finite element code VABS developed in the above mentioned papers [11], [28], [29], [30] cannot be applied to the cross-sectional problem because it requires a large computational capacity. The presence of this small parameter allows however an additional asymptotic analysis of the cross-sectional problems to simplify them. By solving the cell problems imposed with the periodic boundary conditions according to the homogenization technique [8], [22], one finds the effective coefficients in the homogenized cross-sectional problems, which can then be solved analytically or numerically (cf. also [20]). The aim of this paper is to derive and solve the cell and homogenized cross-sectional problems for unidirectional FRC beams whose cross section has the periodic microstructure. For simplicity, we will assume that both matrix and fibers are elastically isotropic but have different Poisson’s ratio. The solution of the cell problems found with the finite element method is used to calculate the asymptotically exact energy and the extension and bending stiffnesses in the 1-D theory of FRC beams. Thus, we determine the dependence of the latter quantities on the shape and volume fraction of the fibers and on the difference in Poisson’s ratio of fibers and matrix that solves the above mentioned issue.

The paper is organized as follows. After this short introduction the variational formulation of the problem is given in Section 2. Sections 3 Dimension reduction, 4 Homogenization are devoted to the multi-scaled asymptotic analysis of the energy functional of FRC beams leading to the cross-sectional and cell problems. In Section 5 the cell problems are solved by the finite element method. Section 6 provides the solutions of the homogenized cross-sectional problems. Section 7 present one-dimensional theory of FRC beams. Finally, Section 8 concludes the paper.

Section snippets

Variational formulation for FRC beams

Consider a FRC beam that occupies the domain B=A×(0,L) of the 3-D euclidean space in its undeformed state. Let x3x be the coordinate along the beam axis. The cross section of the beam in the (x1,x2)-plane, A, consist of two separated 2-D sub-domains Am and Af such that the matrix occupies the domain Am×(0,L), while the uni-directional fibers occupy the domain Af×(0,L). We choose the origin of the (x1,x2)-coordinates so that it matches the centroid of A. We assume that the fibers are

Dimension reduction

Before starting the asymptotic analysis of the energy functional in the limit h0 let us transform the stored energy density to another form more convenient for this purpose [14]. We note that, among terms of W(z,ε), the derivatives w(α;β)/h in εαβ and w;α/h in εα3 are the main ones in the asymptotic sense. Therefore it is convenient to single out the components εαβ and εα3 in the stored energy density. We represent the latter as the sum of three quadratic forms W,W, and W according toW=minε

Homogenization

Consider now the other limit 0. In this limit y plays the role of the “slow” variable, while z=y/ becomes the fast variable. We start with the cross-sectional problem of minimizing functional (8) among w(y,z,x) satisfying the constraint (5)3, where μ(z), expressed in terms of the fast variable z=y/, is a doubly periodic function with period 1. Since x is fixed in this cross-sectional problem, we shall drop this formal variable of w in this Section. Following the homogenization technique

Numerical solution of the cell problems

It is convenient to rewrite functional (22) by changing the sign of χα(y)12Cαβγδ(z)(-χ(α|β)+ε¯αβ+ααβ(z)ξ)(-χ(γ|δ)+ε¯γδ+αγδ(z)ξ),whereCαβγδ(z)=λ(z)δαβδγδ+μ(z)(δαγδβγ+δαδδβγ),ααβ(z)=ν(z)δαβ.

The minimizer of functional (26) satisfies the variational equationCαβγδ(z)χ(γ|δ)δχα|β=Cαβγδ(z)(ε¯γδ+αγδ(z)ξ)δχα|βfor all doubly periodic functions δχα. Eq. (27) will be solved by the finite element method [12], [1]. For this purpose it is convenient to change from tensor notation to matrix notation [12]

Solution of the homogenized cross-sectional problems

Let us first consider the homogenized cross-sectional problem (23) with the average transverse energy being given by (25). It is convenient to present the latter formula in the formW(ε¯αβ,ξ)=12Cαβγδ(ε¯αβ+ναβξ)(ε¯γδ+νγδξ)+12Hξ2,where ξ=γ+hΩαyα. Note that the second term in (32) does not depend on ψα(y), so we need to minimize functional (23) that contains just the first term. We want to show that the minimum of the functional (23) that contains just the first term among ψα(y) satisfying

1-D beam theory

Summarizing the results obtained in Sections 3 Dimension reduction, 4 Homogenization, 5 Numerical solution of the cell problems, 6 Solution of the homogenized cross-sectional problems, we can now reduce the 3-D problem of equilibrium of the FRC beam to the following 1-D variational problem: Minimize the energy functionalJ[v,φ]=0LΦ(γ,Ωα,Ω)dx-f·v(L)+mαvα,x(L)-mφ(L)among functions v and φ satisfying the kinematic boundary conditionsv(0)=0,vα,x(0)=0,φ(0)=0.

In this functional the stored energy

Conclusion

It is shown in this paper that the rigorous first order approximate 1-D theory of thin FRC beams can be derived from the exact 3-D elasticity theory by the variational-asymptotic method. The developed finite element code can be used to solve cross-sectional problems with arbitrary elastic moduli of the fibers and matrix as well as arbitrary distributions and shapes of fibers. The extension of this multi-scaled asymptotic analysis to curved and naturally twisted FRC beams is straightforward [14]

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

CRediT authorship contribution statement

K.C. Le: Conceptualization, Methodology. T.M. Tran: Software, Validation.

References (31)

Cited by (4)

View full text