Elsevier

Applied Mathematical Modelling

Volume 88, December 2020, Pages 161-174
Applied Mathematical Modelling

Longitudinal modeling and properties tailoring of functionally graded carbon nanotube reinforced composite beams: A novel approach

https://doi.org/10.1016/j.apm.2020.06.043Get rights and content

Highlights

  • A novel finite element modeling for axially FG-CNTRC beams was developed.

  • An Advanced numerical solution was proposed based on the ladder-pattern concept to overcome properties discontinuity.

  • Numerical analysis was performed, including convergence test, static, and dynamic analysis with parametric studies.

  • A longitudinal tailoring approach was introduced for AFG-CNTRC beams based on CNT orientation and gradation index.

Abstract

Carbon Nanotubes (CNTs) are considered a promising reinforcement for engineering advanced structures, such as aerospace applications. It became inevitable to tailor structures not only by conventional geometric but also by mechanical properties tailoring, to acquire the maximum load-carrying capacity. Therefore, this research presented a novel Finite Element (FE) modeling of Axially Functionally Graded CNT Reinforced Composite (AFG-CNTRC) beams, based on Timoshenko beam theory. However, an obstacle is raised due to properties discontinuity produced between consecutive elements in axial gradation. A solution is proposed based on the ladder-pattern concept. It defined a locally specific value at each element, satisfying the overall required distribution. It followed by a convergence test, static, and dynamic analysis with comprehensive parametric studies. The obtained results showed that the model is convergent, accurate, free of shear-locking, and suitable for axial gradation. Moreover, this study introduced an advanced technique for axial properties tailoring, utilizing privileges of CNT and aggregation. The longitudinal tailoring depends on both the CNT orientation angle and gradation index, which played a crucial role in response control and adjusting beam stiffness together with strength. This approach could solve enormous engineering problems, especially at critical cross-sections design and constraint response, without affecting the total weight.

Introduction

Recently, CNT Reinforced Composite (CNTRC) structures are considered one of the most influential research fields. Its usage spread widely in many applications. One of the most relevant candidates is aerospace structures. Although, this extensive usage existence, CNT still suffered from issues that represent a processing challenge in nanocomposites such as poor load transfer, dispersion, and interracial engineering [1]. Many types of researches discussed CNT after Iijima [2] in 1991, who is the first one observed it. Volder et al. [3] reviewed the different applications for CNT and the perspective of future ones. Thostenson et al. [4] and Kinloch et al. [1] both made an outlook on CNT composites advances and technologies. Griebel et al. [5] and Han et al. [6] performed Molecular Dynamic (MD) simulations to obtain the elastic properties of polymer/CNT composites. Lei et al. [7] introduced a multi-scale MD finite element method for modeling and analyzing the mechanical response of CNT.

Numerous investigations have been done on CNTs, concluding significant physical and mechanical properties of the new carbon form. It has unique electronic properties and thermal conductivity. Also, stiffness and strength are higher than any current materials. For example, theoretical and experimental results show that modulus of elasticity is greater than 1 TPa, which can not be compared with any other material except diamond 1.2 TPa and specific strength higher than the strongest steel ten to hundred times [4]. Moreover, low specific weight introduces it to be reinforcement for composite materials and nominates it in aerospace applications. At the same time, FG material structures can withstand the different types of loads by properties tailoring according to applied forces and moments. This ability allows it to have an equal interest in research fields. Especially at aircraft structures, the loading cases are varied along with the vehicle lateral direction or the local longitudinal direction of the specific structural part. For example, aerodynamic and inertia loads changed along the airplane wing. It assumed to be carried out by spars or stringers and by longerons for the fuselage case. Tailored AFG-CNTRC could be an excellent solution for adjusting beam stiffness and strength according to applied loads and utilizing CNT superior performance as structural reinforcement. All these privileges proposed without changing in the structural weight, wherein it is a governing factor in the aerospace engineering field.

The combination of advantages from CNT and FG to form Functionally Graded Carbon Nanotube Reinforced Composite (FG-CNTRC) appeared in many kinds of researches. Nguyen et al. [8] introduced a new method to analyze the non-linear vibration of defected FG-CNTRC curve shell subjected to blast and thermal loads. Also, Nguyen et al. [9] studied non-linear post-buckling of FG-CNTRC shell using NonUniform Rational B-Spline. Zghal et al. [10] statically analyzed FG-CNTRC shell and plate based on discrete double directors shell FE. They also used the same theory to investigate the free vibration of FG-CNTRC shell in Ref. [11]. Furthermore, the study of the non-linear bending of nanocomposites shell using discrete double directors shell FE and membrane enhancement had been introduced in Ref. [12]. Moreover, the mechanical buckling of FG-CNTRC plates and curved panels governed by power-law using the same double directors FE shell model was investigated in Ref. [13]. Frikha et al. [14] analyzed FG-CNTRC shell using three and four-node elements. After that, the same authors, [15] studied the dynamic behavior of FG-CNTRC plate and shell using double directors shell FE. Karami et al. [16] studied FG-CNTRC plates analytically. Ansari et al. [17] analyzed forced-vibration of FG-CNTRC plates using a numerical strategy. Also, Ansari et al. [18] studied the vibration of FG-CNTRC elliptic plate numerically. Lei et al. [19] analyzed large deformation of FG-CNTRC plates using element free kp-Ritz method. Moreover, Mellouli et al. [20] used the mesh-free radial interpolation point method in free vibration analysis of FG-CNTRC shell. Also, Mallek et al. [21] analyzed FG-CNTRC shell with piezoelectric layers bonded to the surface.

In the same context, Yas et al. [22] performed a dynamic analysis of randomly oriented FG-CNTRC beams subjected to a moving load. Ke et al. [23] conducted a parametric study for free vibration of FG-CNTRC beams using the Timoshenko beam and Von Karman geometric non-linearity, then analyzed FG-CNTRC beams for dynamic stability [24]. Ansari et al. [25] performed non-linear forced vibration analysis of FG-CNTRC Timoshenko beams. Zeighampour et al. [26] investigated the free vibration of AFG nanobeams with a variable radius using strain gradient theory. Wu et al. [27] studied free vibration and buckling of FG-CNTRC face sheet on a sandwich beam. And then, the same authors [28] represented non-linear vibration analysis of imperfect shear deformable FG-CNTRC beams. Moreover, Wu et al. [29] studied the harmonic resonance of the non-linear FG-CNTRC beams. Ranjbar et al. [30] analyzed AFG-CNTRC cantilever beam subjected to low-velocity impact. Vo-Duy et al. [31] applied a FE method to analyze the free vibration of FG-CNTRC beams. Heidari et al. [32] analyzed the free vibration of rotating FG-CNTRC Timoshenko beams. For tailoring approach, Tatting et al. [33], [34] performed FE analysis for elastic tailored plate design and manufacture.

Minimal kinds of researches dealt from far with modeling and analysis of AFG-CNTRC structures such as Ref. [26], [30], [35]. None of these researches represented a complete FE modeling rather than specific cases of study such as buckling analysis at Ref. [35] or dynamic analysis at Ref. [26], [30]. Presented research introduced a novel complete beam element for AFG-CNTRC as a potential alternative with broad analysis capabilities. Furthermore, the analysis method obtained advantages over the existed in literature. Such as the differential quadrature method, that is suffered from a shortage of systematic analysis conclusions as Ref. [36], [37], [38]. As well, analytical methods, which are incompatible with axial grading modeling as Ref. [16], [39]. Moreover, the proposed model is the first one that planned a solution to overcome the AFG-CNTRC beams convergence problem, using the FE method and ladder-pattern concept. Besides, all reviewed literature did not take into consideration the effect of CNT orientation angle and power-law index change on AFG-CNTRC beams response. This model took into account these parameters and introduced it to the longitudinal tailoring process.

Nowadays, traditional geometric tailoring becomes insufficient for the required progress in the design of the structures. This study aimed to present a novel approach for modeling and properties tailoring of FG-CNTRC beams in the longitudinal direction. It selected the beam as it is one of the most essential and most straightforward structural parts. Investigations started with theoretical and FE modeling of AFG-CNTRC beams. Then, a longitudinal FE problem is explained, and solutions are introduced. Convergence checks are performed to be sure from model efficiency and overcome difficulties of the axial gradation in FE modeling. After that, comprehensive parametric studies are done, including all essential affecting factors in response control and overall beam stiffness. Moreover, enhancement of strength and total mechanical properties due to using CNT as reinforcement, the proposed gradation of properties technique offered the ability of structure tailoring according to applied loads, required response, and specific weight limitations.

Section snippets

Longitudinal gradation concept

Geometric characteristics are shown in Fig. 1. L is the length through x-direction, b is the width in the y-direction, and h is the thickness in the z-direction. Coordinate system origin located at mid-point of the three parameters L, b, and h. Fig. 2(a) showed the proposed distributions in (x, z) longitudinal plan. Two types of aligned CNT are considered, uniformly distributed (UD-CNTRC) only for validation, and axially FG (AFG-CNTRC).

  • – For UD-CNTVCNT(x)=Vtcnt(L2xL2)

  • – For AFG-CNTVCNT(x)=(12

Governing equations and element derivation

The two-node element with six degrees of freedom is selected, as shown in Fig. 3. Furthermore, the theoretical formulations are derived based on the Timoshenko Beam Theory (TBT) as Ref. [44].

Result analysis and validation

A MATLAB code is conducted to perform the following analysis. Due to the lack of published results about AFG beams, UD-CNTRC is used to validate the presented FE model. This study included model convergence, natural frequencies, static deformation, and stress distributions analysis. Several parameters are used to explore the main factors affecting the tailoring process. Boundary conditions of the beams are considered as Clamped-Clamped (C-C) u(0)=u(L)=w(0)=w(L)=ϕ(0)=ϕ(L)=0, Clamped-Simply

Conclusion

The presented model converged at a reasonable number of elements. Also, it is proved to be free of shear locking or divergence. The idea of ladder-pattern and defining local constant property value of each element is reasonable and applicable. Moreover, it successfully overcame the axial gradation modeling problems. The UD-CNTRC beams are used to validate the present model. The obtained results achieved an acceptable agreement.

It can be noticed, great dynamic and static response jump from θ=0o

References (51)

  • A. Frikha et al.

    Finite rotation three and four nodes shell elements for functionally graded carbon nanotubes-reinforced thin composite shells analysis

    Comput. Methods Appl. Mech. Eng.

    (2018)
  • A. Frikha et al.

    Dynamic analysis of functionally graded carbon nanotubes-reinforced plate and shell structures using a double directors finite shell element

    Aerosp. Sci. Technol.

    (2018)
  • B. Karami et al.

    A comprehensive analytical study on functionally graded carbon nanotube-reinforced composite plates

    Aerosp. Sci. Technol.

    (2018)
  • R. Ansari et al.

    Forced vibration analysis of functionally graded carbon nanotube-reinforced composite plates using a numerical strategy

    Physica E

    (2015)
  • R. Ansari et al.

    Vibration analysis of functionally graded carbon nanotube-reinforced composite elliptical plates using a numerical strategy

    Aerosp. Sci. Technol.

    (2017)
  • Z. Lei et al.

    Large deflection analysis of functionally graded carbon nanotube-reinforced composite plates by the element-free KP-Ritz method

    Comput. Methods Appl. Mech. Eng.

    (2013)
  • H. Mallek et al.

    Geometrically non-linear analysis of FG-CNTRC shell structures with surface-bonded piezoelectric layers

    Comput. Methods Appl. Mech. Eng.

    (2019)
  • M. Yas et al.

    Dynamic analysis of functionally graded nanocomposite beams reinforced by randomly oriented carbon nanotube under the action of moving load

    Appl. Math. Model.

    (2012)
  • L.-L. Ke et al.

    Nonlinear free vibration of functionally graded carbon nanotube-reinforced composite beams

    Compos. Struct.

    (2010)
  • R. Ansari et al.

    Nonlinear forced vibration analysis of functionally graded carbon nanotube-reinforced composite timoshenko beams

    Compos. Struct.

    (2014)
  • H. Zeighampour et al.

    Free vibration analysis of axially functionally graded nanobeam with radius varies along the length based on strain gradient theory

    Appl. Math. Model.

    (2015)
  • H. Wu et al.

    Nonlinear vibration of functionally graded carbon nanotube-reinforced composite beams with geometric imperfections

    Compos. Part B

    (2016)
  • Z. Wu et al.

    Nonlinear primary and super-harmonic resonances of functionally graded carbon nanotube reinforced composite beams

    Int. J. Mech. Sci.

    (2019)
  • B. Karami et al.

    Free vibration analysis of triclinic nanobeams based on the differential quadrature method

    Appl. Sci.

    (2019)
  • N. Wattanasakulpong et al.

    Exact solutions for static and dynamic analyses of carbon nanotube-reinforced composite plates with pasternak elastic foundation

    Appl. Math. Model.

    (2015)
  • Cited by (8)

    • Dynamic stability of CNTs-reinforced non-uniform composite beams under axial excitation loading

      2022, Computational Materials Science
      Citation Excerpt :

      Bourada et al. [29] performed the stability and dynamic analyses of CNTs-reinforced concrete beams on elastic foundation. El-Ashmawy and Xu [30] presented a longitudinal modeling and properties tailoring of FG CNTs-reinforced composite beams. Yang et al. [31] investigated the nonlinear flexural behavior of temperature dependent FG CNTs-reinforced laminated beams with negative Poisson’s ratio resting on the Pasternak foundation.

    • Buckling and free vibration of axially functionally graded graphene reinforced nanocomposite beams

      2021, Engineering Structures
      Citation Excerpt :

      Their results demonstrated that the nonuniform distribution of MWCNTs along the axial length yields higher natural frequencies of the microbeams compared to the uniform distribution pattern. El-Ashmawy and Xu [31] presented finite element solutions of AFG CNT reinforced composite beams, and confirmed that axially graded CNTs dramatically promote the static and dynamic performances of nanocomposite beams. It can be concluded that the material gradation along the axial direction of the beam plays a crucial role in the response control and load-carrying capacity adjustment.

    • Mechanical properties improvement of bi-directional functionally graded laminated MWCNT reinforced composite beams using an integrated tailoring–optimization approach

      2021, Microporous and Mesoporous Materials
      Citation Excerpt :

      Many advanced theories were used in Finite Element (FE) analysis of composite structures, such as Generalized Differential Quadrature Method (GDQM) [34,35], Higher-order zig-zag theory [36], and refined high-order global–local theory [37]. However, classical Timoshenko was the common theory used to analyze complex functionally graded nanocomposites, according to literature such as Refs. [38–40]. Moreover, it succeeded in comparison with other higher theories from the accuracy point of view, except shear stress accurate representation.

    • Free vibration of thin-walled open-section beams with functionally graded materials along the contour direction

      2021, Thin-Walled Structures
      Citation Excerpt :

      It can be seen that FGM can be applied to many areas, the number of research publications has increased significantly in the past two decades. The computational modeling of FGM is an important tool for understanding static, vibration and buckling analyses, and has been the target of many researchs [4–14]. Simsek and Kocaturk [15] investigated the free vibration of functionally graded (FG) beams subjected to concentrated moving harmonic load.

    View all citing articles on Scopus
    View full text