Unfolding the mechanical properties of buckypaper composites: nano- to macro-scale coupled atomistic-continuum simulations

Carbon-based nanostructures are receiving increasing attention over the past two decades due to their unprecedented multi-functional features. However, the macro-scale structural applications of these nanostructures have not yet come to full fruition due to the involvement of complex multi-scale computations and manufacturing. Recently, the research community has started investigating buckypaper, which can be described as a sheet or membrane developed using a network of bundles of single-wall carbon nanotubes, multi-wall carbon nanotubes, or a mixture of both. This article aims to focus on the computational bridging of different length scales involving six levels in the range of nano- to macro-scale behaviour concerning buckypaper composites. The sequential derivatives of carbon at six levels, as analyzed in this paper, involve graphene, CNT, CNT bundle, buckypaper, and buckypaper composite automotive components. Here, we adopt a coupled atomistic-continuum modelling approach for the multi-level simulations. Graphene, CNTs, and CNT bundles are modelled using atomistic simulations, while the buckypaper and its composites are modelled using equivalent beam representations for the bundles and continuum solid representation for resin. At the macro-scale, an industry-relevant multi-material composite automotive component has been investigated, wherein the buckypaper is proposed to be embedded involving sheet moulding compound and carbon prepreg. The current simulations have led to the determination of mechanical properties at each level of the carbon-based materials and their mutual dependence. The numerical results demonstrate that a buckypaper composite can enhance the natural frequency and stiffness up to 25 and 37%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document} with respect to conventional monolithic metallic designs, while reducing the weight by 57%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document}. Such outcomes lead to the realization that carbon-based nanostructural derivative in the form of buckypaper can significantly improve the mechanical properties of advanced lightweight structural components as reinforcements for the next generation of aerospace and automotive structures.


Introduction
Graphene and its associated nanostructures (GANS) are considered as miracle materials of the future due to their unprecedented mechanical and other multi-functional properties [22,23,63,102,113,116,141]. However, the application of these materials seems to be largely limited to nanotechnology due to computational and fabrication hurdles at a higher length scale. Currently, networked paper-like formats of these nanomaterials are gaining increasing attention for micro-and macro-scale applications. The success of such investigations could result in an ultra-lightweight revolution in the mainstream engineering industries such as aerospace and automobile. The properties of paper formats of graphene and carbon nanotubes (PFGCN) are studied with the aid of experimental and numerical methodologies. Since the GANS are the building blocks of PFGCN materials, the mechanical properties of PFGCN are strongly dependent on the physical and mechanical properties of the underlying GANS. Thus, the numerical simulation of PFGCN will necessitate the numerical simulation of GANS at multiple levels. GANS's numerical simulation involves a high-fidelity representation of GANS molecules, which means that the atoms and their interactions need to be considered in the simulations. The representative volume element approximations (RVE) offer reasonable accuracy in determining the mechanical properties of these nanostructures [21,[24][25][26][27][28][29]. On the other hand, the PFGCN structures can be represented by low-fidelity finite-element models considering the equivalent properties of GANS without the loss of accuracy [157,171]. Thus, the entire simulation involves an atomistic-continuum multi-level numerical framework starting from the atomic bond properties to the network of PFGCN embedded in matrix materials for manufacturing effective macro-scale structural components. We aim to characterize the effective mechanical properties at each level and bridging different length scales for seamless propagation of the nanolevel properties to the macro-scale structure level.
Buckypaper can be described as a sheet or membrane developed using a network of bundles of single-wall carbon nanotubes (SWCNT) or multi-wall carbon nanotubes (MWCNT). The bundles can also be a mix of SWCNTs and MWCNTs to attain a mix of desired properties [37,64]. Rinzler et al. [140] prepared a batch of SWCNT using dual pulsed laser vaporization process and then purified the resulting material using nitric acid. This process resulted in a free-standing mat of SWCNT, which the researchers termed as "buckypaper". High mechanical performance composites of buckypapers have been successfully synthesized using a prepreg approach [34,35]. In general, it is a well-known fact that the CNTs (carbon nanotubes) offer exceptional mechanical properties [9,43,110,134,160]. However, it is nearly impossible to synthesize high-performance polymeric composites with carbon nanotubes (CNT) as primary fillers or fibres due to dispersion and alignment issues. Typically, the attainable and feasible weight fraction of CNTs in the nanocomposites is less than 5% [19,21,130,161]. To realize higher weight fractions of CNTs in composites, researchers resorted to CNT networks in the form of buckypaper. If many CNTs are dispersed in an aqueous solution, the individual CNTs attract each other and form a thin film when the solution is filtered out. Such a thin film of CNTs is termed buckypaper. The buckypaper is also referred to as "free-standing CNT paper" or "film of randomly arranged bundles of CNTs". The chopped arrangements of CNT particles within the buckypaper lead to porosity. The porous sheet of buckypaper has been claimed to be stronger and lighter than typical steel grades [34,35]. Unlike graphene and CNT, the buckypaper membrane can be macroscopic in nature and is easy to handle. The samples of buckypapers of thickness ranging from μ ms to mm are commercially synthesized. The orientation of CNTs in the buckypaper can be controlled by a magnetic field and/or a compressive load during the stage of synthesis. The presence of pure CNTs at appropriate orientations and tight entanglements between CNTs can make buckypapers an order of magnitude stronger, lighter, and stiffer than conventional engineering materials such as steel alloys, aluminium alloys, titanium alloys, glass fibres, aramid fibres, and carbon fibres. Therefore, buckypaper has a tremendous potential to act as the primary reinforcing material in modern reinforced plastics and structural composites.
In the last 2 decades, a large volume of research work has taken place to develop the nanoscale science and engineering of GANS [71,123]. The mechanical properties of individual graphene sheets and CNTs have been widely studied [87,137,158,167] including their composites [29,132]. These studies have also progressed to bilayered [22,23] and multilayered systems [73,111]. Compared to the research works concerning a single layer, bilayer, and multilayer systems, that of the bundled formats is found to be significantly less [39,92]. Even a lesser number of research investigations on PFGCN materials at nanoscale can be traced in the literature [38,184]. The PFGCN materials have been successfully produced in dimensions up to several meters [159,177]. However, the literature pertaining to the mechanical properties of such meter-scale PFGCN materials is very scarce. Therefore, there exists a strong rationale to investigate PFGCN materials at larger dimensions along with scalable production, which is essential to bring graphene-based materials closer to the demands of structural industries [180]. The crucial aspects of multiscale simulations and scale-bridging ranging from graphene (nanoscale) to engineering products involving buckypaper composites (macro-scale) have not been investigated thoroughly. It can be noted in this context that different forms of structural composites are being increasingly utilized for their high specific strength and stiffness along with other attractive mechanical features [41, 63, 66, 74, 80-82, 117, 119, 149, 150, 169]. PFGCN-based composites could further improve the multi-functional mechanical performances of such composites. This paper aims to propose a seamlessly coupled atomistic-continuum simulation for characterizing the mechanical properties of buckypaper by performing simulations at six distinct levels and propagating the equivalent properties hierarchically in the order of graphene, CNT, MWCNT, CNT bundle, buckypaper, and buckypaper composites along with its practical applications as engineering components (refer to Figs. 1 and 4). Hereafter, this paper is organised as follows. The multi-level simulation framework covering nano-to macro-scale is provided systematically in Sect. 2. The simulation methodology for modelling graphene at the atomic level is presented in Sect. 2.1. This involves finding the mechanical equivalence between covalent bond stiffness and inter-atomic potentials. A similar methodology has also been adopted to model SWCNT and is explained in Sect. 2.2. The atomistic simulation methodology has been extended from SWCNT to MWCNT in Sect. 2.3. This section also deals with inter-wall interactions in an MWCNT. Furthermore, the MWCNT bundles that form the backbone of buckypapers have been modelled at the nanoscale as described in Sect. 2.4. Such MWCNT bundles are also referred to as carbon nanoropes (CNR) in the scientific community [15,16]. A customised algorithm has been developed here to distribute nanoropes in a buckypaper randomly. Both wavy and aligned nanoropes have been considered in the buckypaper. The detailed algorithm and a strategy to model buckypaper composite have been explained in Sect. 2.5. At the next stage, the derived buckypaper composite properties have been used to simulate a critical automotive suspension part (control arm), as discussed in Sect. 2.6. The structure-level mechanical performance of the control arm has been assessed using static and free vibration analyses here. The numerical results concerning equivalent mechanical properties at the above-mentioned six levels have been discussed in Sect. 3, covering Young's modulus, shear modulus and Poisson's ratio of graphene, SWCNT, MWCNT, MWCNT bundle, and buckypaper composite. The numerical results of buckypaper composites have been utilized as one of the primary input properties to assess the static and vibrational performances (involving deflection, mode shapes, and natural frequencies) of the control arm, as presented in Sect. 3.6. Finally, Sect. 4 provides the concluding remarks and a critical perspective on the potential impact of this paper. In this section, we provide the simulation framework for evaluating the effective elastic properties of graphene, CNT, MWCNT, CNT bundle, buckypaper, and buckypaper composites systematically. We start with the atomic bond-level properties and subsequently use a mechanical equivalence of such bonds to evaluate the effective elastic properties based on continuum-based finite-element approaches.

Level 1: Atomistic-continuum simulation of graphene
At this level of simulations, the mechanical equivalence of atomic bonds is used to model graphene sheets. The total inter-atomic potential energy can be derived as the sum of three individual energy terms relevant to bonding and nonbonding interactions [31]. Total strain energy (E E ) can be represented as the sum of energy contributions from bending of bonds (E b ), bond stretching (E s ), torsion of bonds (E t ) and nonbonded energies (E nb ) including the van der Waals interactions, the Coulombic energy, and the core repulsion The influence of bending and stretching is significant in case of small deformations as compared to all other energy components [52,154]. The deformation mechanisms for the multiplanar nanostructure are shown in Fig. 2b-d. The out-of-plane angular component becomes zero for monoplanar nanostructures such as graphene [114]. The total interatomic potential energy (E E ) can be represented as where Δl , Δ , and Δ denote the change in bond length, change in in-plane angle, and change in out-of-plane angle, respectively, as shown in Fig. 2. The quantities k r and k are the force constants associated with bond stretching and (1) (2) bond bending, respectively. The first term in Eq. 2 represents strain energy relevant to stretching (E s ), while the other terms represent the strain energies due to in-plane (E bI ) and out-of-plane (E bO ) angle variations, respectively. The out-ofplane components are irrelevant in case of monolayer nanostructures (these are presented in the current paper for the sake of generality). Neglecting the out-of-plane components in monoplanar nanostructures like graphene, we get The force constants of the atomic bonds ( k r and k ) can be expressed in the form of structural equivalence [88]. From a structure-mechanical point of view (refer to Fig. 2e,  is equivalent to Δ for in-plane angle variation (refer to Fig. 2g). Thus comparing Eq. 5 with the expression for strain energy due to in-plane (E bI ) angle variation (refer to Eq. 3), the following relation can be obtained: k = EI l . It shows that there exists a mechanical equivalence between molecular mechanics parameters ( k r and k ) and structural mechanics parameters (EA and EI). Such mechanical equivalence can be used to derive equivalent beam (covalent bond) properties used in further atomistic simulations of graphene, which can essentially be idealized as a hexagonal network of beams. In the current work, the effective elastic moduli of the nanostructures are computed using equivalent beams representing covalent bonds. After establishing the beam-level mechanical equivalence of atomic bonds, here, we describe the finite-element-based continuum-level framework for analyzing the network of such beams. For the atomistic-continuum simulations, the finite-element analysis code OPTISTRUCT [3] has been used, leading to characterization of the elastic behaviour of nanostructures. The covalent bonds are represented by 3D finite-element beams (element type: CBEAM), and the atoms are represented by finite-element nodes. The cross-sectional diameter and Young's modulus (E) of the beam elements are computed using the equations of forcefield constants k r (stretching) and k (torsional), as discussed in the preceding paragraphs. The numerical values of force field constants k r and k are shown in Table 1 [114]. By substituting these values in the equations k r = EA l and k = EI l , essential parameters to model covalent bonds such as beam diameter d and beam Young's modulus E can be calculated. The bond lengths, which are essentially beam lengths, shown in Table 1 are obtained from the literature [114]. In the current work, the mechanical properties such as Young's modulus (E), Poisson's ratio ( ) , and shear modulus (G) have been computed using the atomistic FE models of the single-layer graphene sheet (SLGS) following the standard definition of the elastic moduli. The computed mechanical properties of SLGS are discussed in Sect. 3.1.

Level 2: Atomistic-continuum simulation of SWCNT
A single-wall carbon nanotube (SWCNT) is a cylindrical version of the graphene sheet. The atomistic model of SWCNT is shown in Fig. 1b. Here, the carbon atoms and their inter-atomic bonds have been generated using a similar atomistic computational methodology as developed for single-layer graphene (refer to Sect. 2.1). In the current work, the mechanical properties such as Young's modulus (E), Poisson's ratio ( ) , and shear modulus (G) have been computed using the atomistic FE models of the SWCNT. For the atomistic-continuum simulations, the finite-element analysis code OPTISTRUCT [3] has been used, leading to characterization of the elastic behaviour. The covalent bonds are represented by 3D finite-element beams (element type: CBEAM), and the atoms are represented by finite-element nodes. The cross-sectional diameter and Young's modulus (E) of the beam elements are computed using the equations of force-field constants k r (stretching) and k (torsional). The numerical values of force-field constants k r and k are shown in Table 1 [114]. By substituting these values in the equations k r = EA l and k = EI l , essential parameters to model covalent bonds such as beam diameter d and beam Young's modulus E can be calculated, as presented in Table 1. The bond lengths, which are essentially beam lengths, shown in Table 1 Bond properties for carbon-carbon atomic interaction in graphene [114] k r in N nm −1 k in N nm rad −2 L in nm d in nm 6.52e−7 8.76e−10 0.142 0.146 Table 1 are obtained from the literature [114]. The computed mechanical properties of SWCNTs are presented in Sect. 3.2.

Level 3: Atomistic-continuum simulation of MWCNT
In the current work, up to five walls have been considered in the mechanical analysis of MWCNT. The carbon atoms and their inter-atomic bonds have been generated using a similar atomistic computational methodology as developed for single-layer graphene (refer to Sect. 2.1). The wall-towall interactions in the MWCNTs need to be captured in the atomistic simulations. This can be implemented using the L-J potential idealized as equivalent springs [115]. The equivalent axial force for an L-J potential between pair of atoms (i, j) belonging to different nanosheets can be defined as [145] where r is the atomic displacement along ij (layer-layer length). As per Girifalco et al. [60], the force between the atoms (i, j) can also be represented by where y = r min + r , and r is the atomic displacement along the length . The r min (in Å) is given by 2 1 6 , where = (A∕B) 1∕6 . The quantities B and A are attractive and repulsive constants, respectively. The values of B and A are 24.3 × 10 3 eV Å 12 and 15.4 × 10 3 eV Å 6 , respectively [145]. In the atomistic FE models, spring elements form a non-linear connection between two adjacent layers of the multilayer structure represent L-J interactions. The force-deflection curve for L-J springs has been calculated using the relation in Eq. 7. Within the finite-element analysis code OPTI-STRUCT [3], the L-J springs of interlayer interactions are modelled by element type CBUSH and using the constitutive relation of Eq. 7 as input properties. An MWCNT with five walls and associated van der Waals springs is shown in Fig. 1c. The computed mechanical properties of MWCNTs are presented in Sect. 3.3.

Level 4: Atomistic-continuum simulation of MWCNT bundles
Analogous to carbon fibres produced by intermingling carbon filaments, it is also viable to produce nanoropes by intermingling CNTs. Such an arrangement can be formed by bundling SWCNTs or MWCNTs or a mix of both (refer to Fig. 1d-g). These bundles have been successfully used in various applications, including composite materials [93] and electronic displays [85]. In the current work, it has been assumed that the MWCNT of five walls forms the fundamental element of bundles. This assumption is based on the MWCNTs found in high-performance buckypapers [34,35]. Four types of bundles have been trialled with two, three, five, and seven MWCNTs. Here, the tubes' arrangement in the seven MWCNT configuration is similar to the ones used in buckypapers presented by [1,13]. Liu et al. [92] also modelled similar seven tube arrangements using a hybrid atomistic methodology. In this study, the carbon atoms and their inter-atomic bonds have been generated using a similar atomistic computational methodology as developed for single-layer graphene (refer to Sect. 2.1). The van der Walls interactions between the walls of individual MWCNTs have been modelled as described in Sect. 2.3. Van der Walls interactions have also been considered between the pairs of MWCNTs in the bundle. This is based on the 0.85 nm cut-off distance, beyond which the interaction would be negligible [152]. In the bundles' atomistic FE models, spring elements form a non-linear connection between a pair of MWCNTs, represent the L-J interactions. The force-deflection curve for L-J springs has been calculated by using the relation in Eq. 7. Within the finite-element analysis code OPTISTRUCT [3], the L-J springs of interlayer interactions are modelled by element type CBUSH and by using the constitutive curve obtained from Eq. 7 as input properties. Four variants of MWCNT bundles are shown in Fig. 1.
The computed mechanical properties of MWCNT bundles are presented in Sect. 3.4.

Level 5: Atomistic-continuum simulation of buckypaper composite
To model buckypaper, the fibres of MWCNT bundles are represented as equivalent beams. The concept of equivalent beams in the context of nanotubes has been historically validated [44,125,144,168]. The mechanical properties of MWCNT bundles (such as Young's modulus, shear modulus, and Poisson's ratio) computed in Sect. 3.4 have been used to define the equivalent beams here. A network of high aspect ratio beams forms the buckypaper, wherein each beam represents an MWCNT bundle (refer to Fig. 1h). It is assumed that each bundle is composed of seven MWCNTS, wherein each MWCNT is a five-walled CNT. The assumed arrangement of seven MWCNTs in a bundle is shown in Fig. 1g [1,13,75]. Two configurations of beams have been trialled in the current work, namely stretched beams and wavy beams. Stretched beam configurations are aligned and straight in nature, whereas wavy beam configurations display curvatures. The consideration of such beam configurations mimics the concept of stretch ratio presented by Cheng et al. [34,35], wherein they demonstrated a high mechanical performance of buckypaper composites with the aid of high stretch ratio MWCNT bundles. The stretching of beams can be attained by either a mechanical force [34,35] or by a magnetic field [61]. To demonstrate the influence of stretched beams on the mechanical performance of the buckypaper composite, five configurations have been modelled here. These configurations are based on the ratio between stretched to wavy (un-stretched) beams, as 0, 20, 40, 60, and 80% of stretched beams (refer to Fig. 3a-e, respectively). In these images, wavy beams are represented in blue colour and the stretched beams are represented in red colour. In a buckypaper, the MWCNT beams intersect with each other in various formats, which include entangled crossings, Y and T joints, nanowelding, nanosoldering, and joints of coated tubes. In the case of entangled crossings, the force transfer between beams occurs as a result of frictional contact and also due to weak van der Waals interactions [50,120,178]. The Y and T joints are formed when one end of a beam connects to other beam(s) resembling a Y or T shape [84,106,186]. Nanowelding is specific to irradiation, wherein an intersection between adjacent tubes leads to a defect and then resulting in bond formation [120,163]. Nanosoldering occurs when amorphous carbon deposits within the intersecting region of a pair of tubes [12,120]. Joints of coated tubes are formed due to the depositions of coating materials at the intersections of tubes [12]. Berhan et al. [16] postulated that the above-mentioned formats of tube intersections could be represented by torsional springs.
The stiffness of such torsional springs can be calculated by In the above equation, K, K, E, I and L are normalised rigidity, spring torsional stiffness, beam's Young's modulus, beam's moment of inertia, and beam length, respectively. The normalised rigidity K is dependent on the type of joint and the beam aspect ratio. The values of K with respect to beam aspect ratio have been plotted in the literature [16].
In the current work, the value of K has been used based on entangled and crossed joints as presented in [16].
In the current work, the random network of nanoropes (also known as bundles) has been generated by scripting an algorithm (refer to the supplementary material where the developed code is provided), wherein a combination of stretched and wavy (un-stretched) CNT nanoropes is assumed to be the content of a typical buckypaper [34,35]. In the buckypaper, the stretched and un-stretched nanoropes have been represented as straight lines and sine curves, respectively [15,180]. The steps involved in the automated distribution of nanoropes are given below: • Step 1: The total number, amplitude, and period of sine curves have been defined. The total number and length of straight lines are also defined. The length of straight lines is defined as the same as the period of sine curves. An initial sine curve has been created in the first instance, and then, it has been randomly copied into multiple locations and randomly rotated in each location. Several points have been defined on the sine curves and the straight lines to generate finite-element-based nodes. • Step 2: The points that intersect between curves and straight lines have been identified using the intersection function available in the Matlab library. At the location of these points, the curves and straight lines have been split to create new finite-element nodes. To locate these points in the data generated in Step 1, a user-defined . function findposition has been created. A second userdefined function insertrows has been developed to insert extra points(nodes) into the data generated in Step 1, as a result of intersections. • Step 3: A text file containing the data created in Step 1 and Step 2 has been written as an output of the Matlab code. This text file contains the node, beam element, and spring element data. Each sine curve and the straight line have been represented by 20 nodes and 19 beam elements. The number of nodes and elements has been left open to modification to take into account intersection points. At the intersection points, torsional spring elements have been created. • Step 4: The above three steps have been repeated for different stretched and un-stretched ratios of nanoropes to generate five different buckypapers (refer to Fig. 3a-e).
The Matlab code explained above generates an output text file containing the necessary finite-element input data of buckypapers (involving the definition of nodes, beam elements, and torsional spring elements). The data are saved with ".inp" extension. At the next stage, within the finiteelement code ABAQUS [40], the beam elements and spring elements are modelled using the element types B21 and Spring2, respectively. To model the buckypaper-Bismaleimide (BMI) composite system, it is assumed that the BMI matrix is a continuous structure. Therefore, the matrix is modelled by general purpose 3D triangular shell elements of ABAQUS (S3) [40]. The buckypaper is considered to be sandwiched between two layers of the BMI matrix. Here, two layers of shell elements are modelled on either side of the beam network to represent a 60% fibre weight fraction. The interface between the buckypaper and the BMI matrix is assumed to be perfectly bonded and modelled as a rigid contact using Tie constraint. The BMI matrix is considered to be isotropic with Young's modulus and Poisson's ratio as 4.6 GPa and 0.33, respectively [101]. The sectional view of a mesh of buckypaper composite is shown in Fig. 3f, where the yellow background represents the solid shell mesh of resin. The calculated mechanical properties of buckypaper composites using the model developed here are presented in Sect. 3.5.

Supplementary note: Waviness of nanoropes in the buckypaper
The wavy curvatures of nanoropes in the buckypapers are evident in several SEM images presented in the literature [34,35,96,121]. The curvatures in the geometry of waves can be represented by various functions, including sine, cosine, and exponential. The arbitrary waviness of ropes in the buckypaper might not exactly follow any of these functions. However, various works have successfully approximated and validated sine [15,180] and cosine functions [17,48,48]. The current work is based on sine approximations presented by Berhan et al. [15]. Due to the convenience of representing nanorope waviness by period, amplitude, and frequency of sine waves, these functions have been chosen here. Furthermore, the degree of waviness can be varied by modifying the frequency and period of sine functions. Other means of modelling curvature include waviness factor [47] and a waviness parameter [68,69]. Feng et al.'s [47] model used a waviness factor ( = l nt ∕l f ) to determine the influence of nanorope curvature on CNT composite film modulus, wherein l nt and l f were actual length and distance between two ends of nanoropes, respectively. The straight and wavy nanoropes lead to = 1 and > 1 , respectively. Tsai et al. [68] assumed a sine curve to represent wavy nanorope and defined the waviness parameter as w = A∕L , wherein A and L were the amplitude and wavelength, respectively. The overarching aim of the research community in this field is to minimize waviness and maximize the straightness of nanoropes, since the presence of nanorope waviness leads to a lower buckypaper mechanical property. The waviness can be significantly reduced by methods such as stretching [34,35].

Level 6: Simulation of an automotive suspension control arm
The broader objective of the current work is to introduce graphene-based materials into critical industry-grade structural components for enhanced mechanical performance. As discussed in the preceding sections, neither graphene nor CNT in their nanometer size scales can be directly cast or fabricated into a geometrically complex macro-scale engineering component. Note that the powder format of graphene or CNT nanomaterials does not offer essential mechanical properties. However, the paper variants such as graphene paper or buckypaper can play the role of primary reinforcement material in high-performance engineering components at macro-scale. Therefore, an effort has been made here to demonstrate the suitability of buckypaper as a c Alternating slices of SMC, carbon prepreg, and buckypaper reinforcing layer in composite automotive suspension control arms. The evolution of the buckypaper reinforced composite automotive suspension control arm concept is shown in Fig. 4. A typical control arm (refer to Fig. 4a) of a road vehicle suspension weighing 4.2 kg has been chosen for a baseline design. This baseline design is assumed to be cast using an aluminium alloy. The current work aims to improve this baseline design by replacing aluminium alloy with a buckypaper-based material system. A multi-material system consisting of the sheet moulding compound (SMC), carbon fibre prepreg, and buckypaper composite has been assembled to form the geometry of control arm (refer to Fig. 4b). The buckypaper composite of MWCNT (multi-wall carbon nanotube) and bismaleimide (BMI) resin, as described in the preceding paragraph, is used in the control arm. This composite buckypaper offers a density of 1.36 gm/cm 3 (based on 60% fibre weight fraction) and Young's modulus of 370.0 GPa (refer to Sect. 3.5). Thickness of the composite sheet is 50 μ m, which is calculated based on 60% fibre weight fraction configuration. The SMC chosen is AMC8592, which offers a density of 1.48 gm/cm 3 and Young's modulus 69.0 GPa. The carbon fibre composite chosen is a prepreg with grade HM508552 (unidirectional fibre composite). The resin in this prepreg is Hexply-8552-epoxy, and the fibre volume fraction is 60%. The thickness of this prepreg is 0.13 mm, and the density is 1.5 gm/cm 3 .
An isometric view of the multi-material concept control arm is shown in Fig. 4b. In this isometric view, six slices visible in pink, blue, dark green, cream, light blue, and green colours are the sections of the SMC. The two rings shown in red colour are also made of the SMC. A split view showing all different materials used in the multi-material control arm is presented in Fig. 4c. This view shows 21 separate components and the material selection for each component is shown in Table 2. This results in ten sheets of prepreg and five sheets of buckypaper. The view in Fig. 4b is of the consolidated part, wherein prepreg and buckypaper become hidden underneath, and only the SMC becomes visible in the image. Since the stiffness and strength properties of prepreg and buckypaper are directional, each prepreg sheet is designed as a laminate composed of eight plies oriented at 0 • , 45 • , 60 • , and 90 • , repetitively. Also, each buckypaper sheet is designed as a laminate of 16 plies oriented at 0 • , 45 • , 60 • , and 90 • , repetitively. Such an arrangement of layers of buckypapers and prepregs is made to attain near isotropic two-dimensional material properties. Furthermore, adequate shear properties are also attained due to such laminated composite layers. The two-dimensional arrangements of buckypapers and prepregs in conjunction with bulking of sheet moulding (SMC) compound would offer the resulting component adequate mechanical properties in all directions. Within the finite-element model, the laminates are represented by the conventional shell elements and the SMC is represented by solid elements.
In the current investigation, mechanical performance of the automotive suspension arm is assessed through vibrational analysis and static analysis for demonstration. Since a control arm undergoes motion only in a few directions, its degrees of freedom are restricted to those few directions. The vibrational analysis involves extracting mode shapes that are relevant to these degrees of freedoms. The mode shapes that need to be extracted are typically vertical bending modes, in-plane bending modes, out-of-plane bending modes, and torsional modes. The mode shapes and their natural frequencies are highly influenced by the component's material properties, geometry of the component design, and the component boundary conditions. Changing the material of a component does not change its mode shape, while the natural frequency will be changed. In automotive industry, the common practice is to extract the first five mode shapes of suspension components along with the corresponding natural frequencies and benchmark them. The benchmarking is usually performed against predetermined targets or baseline designs. The baseline design can be equivalent to an existing similar vehicle manufactured in house or by a competitor firm. The aim of such modal analysis in component design is normally to ensure an increase in the natural frequency of the component compared to benchmark natural frequencies. The higher natural frequency indicates higher stiffness, as well. In the finite-element analysis, the mode shapes are the eigenvectors, and the natural frequencies are the eigenvalues, obtained by solving the equation of motion, as described below In the above equation of motion, [ ] , [ ] , , and are stiffness matrix, mass matrix, frequency of vibration, and displacement vector, respectively. There are various finiteelement-based algorithms available to extract eigenvectors and eigenvalues from the above equation of motion. The algorithms include Lanczos, block Lanczos, automated multi-level sub-structuring, Givens, Householder, inverse power, Sturm modified inverse power, and subspace methods. In the current work, the block Lanczos method has been used to extract mode shapes and natural frequencies of the

Supplementary note: Conceptual design of the automotive control arm
The baseline design shown in Fig. 4a is a CAD model based on a C-segment automobile's existing aluminium control arm. Detailed geometrical features have been removed for the sake of simplicity. The current work aims to replace this baseline metallic design with a multi-material design, wherein the buckypaper would act as a reinforcement. The concept shown in Fig. 4c is a result of prior parametric studies. These parametric design studies considered several different concepts. An initial trial involved slicing off the baseline control arm (refer to Fig. 4a) into two sections and then placing a buckypaper/prepreg layer at the centre. However, this design involving two SMC sections and one buckypaper/ prepreg laminate did not offer the required mechanical performance. Therefore, several iterations have been carried out by progressively increasing slices of SMC and buckypaper/ prepreg laminates. The iterations also varied the number of plies within each laminate based on manual optimisation of ply-stacking sequences. We settled with the arrangement of slices shown in Fig. 4c, since it offered a mechanical performance higher than that of the baseline design. The manual optimisation of ply-stacking sequences concluded the prepreg laminates composed of eight plies repeated at 0 • , 45 • , 60 • and 90 • , and buckypaper laminate of 16 plies repeated at 0 • , 45 • , 60 • , and 90 • as the optimum stacking sequence.

Supplementary note: Role of SMC and prepreg
Although the main aim of the current work is to demonstrate the suitability of buckypaper as a raw material to manufacture an engineering component, for the automotive component chosen here, the buckypaper on its own cannot fulfil all the design needs. For instance, the buckypaper may not conform to the intricate geometric features of the control arm. As a result, the sheet moulding compound (SMC) with the ability to be moulded to complex geometric shapes has been selected as a part of the multi-material component design. Furthermore, the resin of the buckypaper composite may not interact effectively with the SMC. Therefore, a carbon fibre-epoxy prepreg has been selected, since the epoxy resin has the ability to strongly bond to various types of materials. Such firm adhesive bonding properties have been historically validated [97] and accepted industry-wide. SMC is a layered composite system involving a polymeric matrix reinforced by chopped long fibres. The SMC sheets on rollers are available to be purchased from various composite suppliers, including Hexcel Inc and Toray Inc. The gel-like texture of the SMC facilitates this raw material to This grade is composed of chopped carbon fibres with 1 ′′ length at 53 wt% and vinyl-ester resin [129]. A prepreg is a composite system of continuous fibres and polymeric resin. Uncured continuous sheets of prepreg are supplied on rollers by raw material suppliers, including Hexcel Inc and Toray Inc. These materials can be laminated, cured, and moulded into engineering components. The grade of prepreg used in the current work is HM508552. This prepreg involves unidirectional carbon fibres at 60 vol% and epoxy resin [67]. Compression moulding of carbon fibre prepreg and carbon fibre SMC into complex shapes has been demonstrated by the research community widely [175]. The moulding of the proposed multi-material control arm is illustrated in the Fig. 5. Here, the slices of SMC and the laminated layers of prepreg and buckypaper can be assembled and compression moulded using the process demonstrated by Wulfsberg et al. [175].

Results and discussion
In this section, we have systematically presented the numerical results for investigating the mechanical behaviour of all the carbon-based structural components at multiple length scales, starting from graphene to the buckypaper composites.

Mechanical properties of graphene
An atomistic finite-element representation of the graphene sheet can be seen in Fig. 1a. First, the atomistic finite element models of various sizes have been constructed to identify the sensitivity of Young's modulus to the sheet dimensions. To calculate the effective Young's modulus, the graphene sheet has been constrained at one edge, and a load is applied on the opposite edge. The resulting tensile sheet deflection has been used to calculate Young's modulus using the relation In the above equation, E, F, L, A, are Young's modulus of graphene, applied force, sheet length, sheet cross-section area, and computed sheet deflection, respectively. The variation of the E with respect to the sheet length is shown in Fig. 6a. At smaller lengths, the modulus increases with the length due to the prominence of boundary effect. However, at larger lengths, the modulus becomes constant with respect to the length of the periodic nanostructure. The Young's modulus at larger lengths is 1082 GPa and is comparable with the literature, as shown in Table 3.
To measure shear modulus using atomistic simulations, one edge of the graphene sheet has been constrained, and the opposite edge has been subjected to shear force. The resulting shearing deflection has been recorded to calculate shear modulus as In the above equation, G, F, l, A, s are the shear modulus, applied shear force, sheet length parallel to the shearing direction, sheet cross-section area, and computed sheet deflection, respectively. The variation of the G with respect to the sheet length is shown in Fig. 6a. At smaller lengths, the shear modulus increases with the length due to boundary effect. However, at larger lengths, the shear modulus of the periodic structure remains constant with respect to the length. The shear modulus at larger lengths is 606 GPa and it is in agreement with the literature (refer to Table 4). Using the atomistic simulations, Poisson's ratio of graphene sheet has been calculated by introducing tensile strain at one edge and then recording the longitudinal and lateral strains. Poisson's ratio is the negative ratio between lateral and longitudinal strains. The computed Poisson's ratio for graphene sheets is presented in Fig. 6b. The Poisson's ratio computed for a single-layer graphene sheet here is 0.62. This value has been compared against the available literature in Table 5. wherein a rather large range of values can be noticed. However, the presently computed Poisson's ratio falls well within the range.

Mechanical properties of SWCNT
The mechanical properties predicted for SWCNTs are found to be exceptionally high. However, these mechanical  Tables 6, 7 and 8, wherein these are found to be well within acceptable limits.
In the atomistic simulations, the Young's modulus of SWCNT has been calculated by constraining the tube's one end, applying tensile force at the tube's opposite end and then using the resulting deflection in Eq. 10. The variation of E with respect to the tube length is shown in Fig. 6c. The value of E computed for the tube length of 54 nm is 1080 GPa. It is important to note that a similar Young's modulus value has been offered by SLGS, as shown in Sect. 3.1. The value of E obtained in the current work has been compared against the values presented in the literature (refer to Table 6), where value of E up to 5300 GPa has been reported. Avila and Lacerda [9] presented such a large value based on the tensile analysis of an SWCNT zigzag (14,0) configuration with tube length 2.4 nm. In the literature (Table 6), the lowest value of E is 100 GPa, presented by Ge et al. [51]. They modelled a SWCNT of 1 nm diameter and measured Young's modulus in the radial direction, resulting in the lowest magnitude. As per the literature, E varies from 100 to 5300 GPa. The Young's modulus computed in the current work is 1080 GPa ( Table 6) and falls well inside this range.
The shear modulus of the SWCNT has been calculated using the principle of torsional rigidity. A similar method has also been implemented by Zaeri et al. [180]. One end of the SWCNT has been constrained, a torsional load has been applied to the other end, and the resulting tangential deflection has been recorded. The shear modulus of SWCNT can be calculated by In the above equation, T is the torque, l is the length of SWCNT, is the twist angle, and J is the polar moment of inertia. The variation of shear modulus G with respect to the SWCNT length is shown in Fig. 6c. As per this plot, G does not tend to be influenced by the tube length. The G here is found to be 504 GPa, which can be noted to be significantly different compared to that of SLGS. The computed value of G in the current work has been compared against the literature in Table 7. As per the literature, the shear modulus of SWCNT can be as high as 900 GPa. This indicates that SWCNT can offer shear stiffness more than 11 times higher than that of structural steel. Schiebold and Mehner [148] predicted the shear modulus to be 900 GPa by considering an equivalent cylindrical solid mesh of SWCNT with 0.5 nm diameter. As per the literature, the shear modulus (G) of SWCNT can be as low as 320.0 GPa [103], where a springbased finite-element model of SWCNT was used with 2 nm diameter. It can be concluded that the G of SWCNT ranges from 320 to 900 GPa in the literature. The computed value 504 GPa in the current work falls well within this range. Poisson's ratio ( ) has been measured by introducing tensile forces to the SWCNT, followed by recording transverse and axial strains. The computed Poisson's ratio as a function of SWCNT length is shown in Fig. 6b. At lower lengths, Poisson's ratio drops with the increase in length up to 25 nm. At higher lengths (i.e., > 25 nm), is found to be independent of the SWCNT length. It can be concluded that the remains stable at 0.37 at lengths beyond 25 nm. Therefore, this value can be used to model high aspect ratio SWCNTs within a buckypaper. The computed value of Poisson's ratio ( ) has been compared against literature in Table 8. As per the literature, of SWCNT can be as low as 0.06. Such a low value has been reported by Domínguez-Rodríguez et al. [43] for a SWCNT with diameter ≈ 0.2 nm. The largest value of identified in the literature (Table 8) is 0.96. Askari and Ghasemi-Nejhad [7] determined the to be 0.96 with the aid of a space frame representation of SWCNT. The element type used to represent C-C bonds in this analysis was linktype instead of equivalent beams used in the current work. As per the literature data presented in Table 8, can range between 0.06 and 0.96. The computed value of 0.37 in the current work falls well within this range.

Mechanical properties of MWCNT
The finite-element models of CNTs with up to five walls have been considered in the current work. The variation and sensitivity of E to the number of walls in a MWCNT are presented in Fig. 7a. In the current analysis, E of SWCNT is found to be 1080.0 GPa, whereas it is 1206.0 GPa for a MWCNT with double walls. For the MWCNTs with three, four, and five walls, E is found to be 1292.0, 1354.0, and 1401.0 GPa, respectively. There is a significant increase in the modulus of double wall carbon nanotubes (DWCNT) compared to SWCNT. Beyond two walls, the slope of this curve is not significantly high. The computed E of MWCNTs for various numbers of walls has been compared against the literature data in Table 9. From the literature, it is evident that E of MWCNTs can be up to 1600 GPa. This means MWCNTs can be 7.5 times stiffer compared to structural steel under tension. Table 9 indicates that the E of MWCNT can range from 953.1 to 1600.0 GPa. The computed E of MWCNTs in the current work falls well within this range. The shear modulus G of MWCNTs has been calculated by applying torque and then using Eq. 12. The variation of G with respect to the number of MWCNT walls is shown in Fig. 7a. This plot follows a similar trend as that of E for the MWCNTs. The G of DWCNT is found to be significantly higher compared to that of SWCNT. It can also be observed that the G increases with the number of walls. However, the rise in G from four-walled CNT to five-walled CNT is found to be moderate. The G of MWCNTs with two, three, four, and five walls are found to be 556.0, 598.0, 622.0, and 638.0 GPa, respectively. These computed values have been compared against the literature data in Table 10. Note that the literature data are based on only two-, three-, and fourwalled CNTs. The literature suggests that the G of MWCNTs can range from 360.0 to 514.0 GPa, which is in good agreement with the current predictions.
Poisson's ratio ( ) has been measured by introducing tensile forces to the MWCNT, followed by recording the transverse and axial strains. The computed Poisson's ratio as a function of the number of MWCNT walls is shown in Fig. 7b. According to this plot, as the number of walls is increased, drops initially and then converges. The drops in from SWCNT to DWCNT and from DWCNT to triplewalled CNT are considerable. However, the drops in from triple-walled CNT to four-walled CNT and from four-walled CNT to five-walled CNT are negligible. Therefore, it can be assumed that Poisson's ratio remains constant when the number of walls is increased beyond three. The Poisson's ratio for double-walled, triple-walled, four-walled, and five-walled MWCNTs are found to be 0.35, 0.33, 0.326, and 0.324, respectively. These data have been listed along with the values from literature in Table 11. Zhang et al. [183] presented a value of 0.42 for an armchair DWCNT with 3 nm diameter. From Table 11, it is evident that the values of for MWCNTs computed in the current work are comparable to the one presented in published literature.

Mechanical properties of MWCNT bundle
In the current work, four different configurations of the MWCNT bundles have been investigated. These include a bundle with two MWCNTs (configuration(2), refer to Fig. 1d), a bundle with three MWCNTs (configuration(3), refer to   Fig. 1e), a bundle with five MWCNTs (configuration(5), refer to Fig. 1f), and a bundle with seven MWCNTs (configuration (7), refer to Fig. 1g). The Young's modulus (E) has been recorded by introducing boundary conditions same as in Sect. 3.2 and using Eq. 10. The E for configuration (2), configuration(3), configuration (5), and configuration (7) of MWCNTs    Fig. 7c. The slope of the curve is not significant in this plot. The computed value of E for MWCNT bundles has been compared against literature in Table 12. Liu et al. [92] used a combination of atomistic and continuum mechanics to predict E of SWCNT bundles, while Bai et al. [11] fabricated several cm long MWCNT bundles and performed tensile tests. The E computed in the current work is comparable with the literature data presented in Table 12. The modulus of MWCNT bundle identified here is significantly higher as compared to that of SWCNT fibre system presented by Dalton et al. [39]. Dalton et al. [39] produced a SWCNT fibre system by spinning tubes within an epoxy resin and determined the modulus to be 80 GPa. The buckypaper modelled in Sect. 3.5 is assumed to be composed of MWCNT bundles with seven tubes. Therefore, the E of MWCNT bundles considered in the subsequent buckypaper analysis is 1461.0 GPa. It can be noted that the process of irradiation [105] can introduce covalent bonds between the tubes within the bundle and also between the bundles of a buckypaper. As a result, the modulus of the bundle can be enhanced further [1]. The shear modulus (G) of MWCNT bundles has been calculated by introducing shearing boundary conditions similar to the calculation of SWCNT (refer to Sect. 3.2) and using Eq. 12. All four configurations of MWCNT bundles shown in Fig. 1 have been analyzed under shear loading. The G for configuration(2), configuration(3), configuration(5), and configuration (7) of MWCNTs are found to be 652.0, 664.0, 671.0, and 675.0 GPa, respectively. This indicates a marginal increase in shear stiffness with the addition of tubes. The variation of G with respect to the number of tubes in the bundle is shown in Fig. 7c. The slope of the curve is not significant in this plot. No literature data exist that has analyzed the current configurations of MWCNT bundles under shear. However, Liu et al. [92] analyzed a bundle of seven SWCNTs and presented the value of G as 7.0 GPa, where the shear loading on the SWCNT bundle was significantly different compared to that of the current work. The buckypaper modelled in Sect. 3.5 is assumed to be composed of MWCNT bundles with seven tubes. Therefore, the G of MWCNT bundles considered in the subsequent buckypaper analysis is 675.0 GPa.
Poisson's ratio has been measured by introducing tensile forces to a MWCNT bundle, followed by the recording of transverse and axial strains. The has been computed for all four configurations of MWCNT bundles shown in Fig. 1. The Poisson's ratio as a function of the number of MWCNTs in a bundle is shown in Fig. 7d. The for configuration(2), configuration(3), configuration(5), and configuration(7) of MWCNTs are found to be 0.33, 0.34, 0.345, and 0.35, respectively. This shows a minimal increase in with the addition of MWCNTs into the bundle. No literature that determines Poisson's ratio of MWCNT bundle has been identified except that of Miao et al. [109], which used a CNT forest to dry spin the bundles into yarn and concluded a Poisson's ratio as high as 8.0. Such an unusually high value is the result of twisted configurations and porosity in the yarns. The buckypaper modelled in Sect. 3.5 is assumed to be composed of MWCNT bundles with seven tubes. Therefore, the of MWCNT bundles considered in the subsequent buckypaper analysis is 0.35.

Mechanical properties of MWCNT buckypaper composite
The longitudinal and lateral moduli, denoted here as E 1 and E 2 , respectively, have been computed for five variations of buckypaper composite. The five variants are based on stretched compositions of 0, 20, 40, 60, and 80%, as shown in Fig. 3a-e. A cross-sectional view of the mesh with 100% wavy nanoropes, taken at the midplane, is shown in Fig. 3f. The finite-element-based modelling strategy of this buckypaper composite has been explained in Sect. 2.5. The variation of both longitudinal and transverse modulus with respect to the stretch concentration of buckypaper composites is shown in Fig. 8a. The longitudinal modulus is measured along the length of stretched nanoropes (i.e., CNT bundles), while the transverse modulus is measured perpendicular to the length of stretched nanoropes. In Fig. 8a, the presence of 20% stretched nanoropes leads to a longitudinal modulus of 80 GPa. This is higher than the Young's modulus offered by engineering aluminium. Importantly, the presence of 80% stretched nanoropes would lead to a longitudinal modulus of 370 GPa. Such a value of modulus is higher than that of engineering steel and carbon fibre. Figure 8a shows a buckypaper composite consisting of only un-stretched (wavy) nanoropes, where the observed longitudinal modulus is 45 GPa. It indicates that, by stretching nanoropes, an enhancement of modulus up-to eightfold can be achieved. Such a phenomenon can be attributed to the fibre length and orientation efficiency factors described in classical composite theories [78]. A comprehensive comparison against the literature concerning different classes of buckypaper is shown in Table 13. The modulus as per this table varies between 0.13 and 650.0 GPa, indicating a considerable scatter. However, the modulus of the current work (80% stretched) is comparable to that of the other studies (Cheng et al. [35]), which also made use of the BMI as the resin in conjunction with an MWCNT buckypaper. This experimental work attained a modulus of 350.0 GPa by functionalizing MWCNT paper with m-CPBA/dichloromethane solution and performing mechanical stretching. However, without functionalizing and mechanical stretching, the modulus of composite paper was found to be 40.0 GPa in [35], which is comparable with the current work (un-stretched). Composite sheet shear modulus has been computed by constraining one edge and then introducing a shearing load along the length of the opposite edge. This simulation replicates a test sheet behaviour of ASTM D7078. The shear modulus has been computed for sheets with five stretch compositions of 0, 20, 40, 60, and 80%, as shown in Fig. 8b. The presence of wavy nanoropes offers higher in-plane shear modulus compared to that of stretched nanoropes in either of the two plane axes. Therefore, the buckypaper composite with 0% stretched nanoropes offers the highest in-plane shear modulus (27 GPa). As the stretch composition is increased, the shear modulus drops significantly and reaches as low ad 4 GPa for the buckypaper composite with 80% stretched and 20% wavy nanoropes. This value of shear modulus will be used to further simulate buckypaper composite in the automotive suspension control arm presented in Sect. 3.6.
It may be noted that there are no literature data available on the shear modulus of buckypaper composites.
Poisson's ratio of the buckypaper composite has been measured by applying tensile force at one edge and constraining the opposite edge. The computed Poisson's ratio of buckypaper composite has been plotted in Fig. 8c. This plot shows a variation of Poisson's ratio ( ) with respect to the % of stretched nanoropes. Note that the stretched nanoropes are the fibres that are aligned along the direction of loading. As the content of stretched nanoropes is increased, stiffness along the stretch direction increases. This leads to a drop in Poisson's ratio, as noticeable in Fig. 8c. At 0% stretched nanoropes' content, the Poisson's ratio is 0.52. The value of the Poisson's ratio with 80% stretched nanorope content is 0.47. Thus, a drop of 9.6% in Poisson's ratio has been observed due to an enhancement in the alignment of nanoropes by 80%. The Poisson's ratio determined in the current work has been compared against that of literature data in Table 14. The literature is available on a variety of buckypapers involving SWCNT or MWCNT, or graphene. The paper format of graphene is simply referred to as graphene paper in many articles. The Poisson's ratio of buckypaper composite is found to vary between 0.07 and 3.4 in the literature (Table 14). The Poisson's ratio determined in the current work falls well within this range. It is important to highlight that the reference [155] measured Poisson's ratio of MWCNT buckypaper/polyurethane composite under compressive loading. None of the articles in the literature is found to consider bismaleimide (BMI) as a resin, unlike the current work. The other research studies used different variants of the polymer as bulk materials in the composites to determine Poisson's ratio. The value of Poisson's ratio determined here will be used to further simulate the automotive suspension control arm in Sect. 3.6.

Supplementary note: Strength of nanocomposites
In this paper, the primary focus is on evaluating the effective elastic properties, leading to stiffness investigation at different levels of the length scales. However, the strength of raw materials also plays a vital role in the mechanical performance of an engineering component such as a control arm. A brief tensile analysis has been presented in this section to identify the value of strength buckypaper composites can offer to engineer components. The multi-scale modelling strategy used here is the same as our previous works [27][28][29]. The analysis is based on an idealized RVE, wherein

3
CNTs are represented by traditional atomistic space frame models (described in Sect. 2.2) and the surrounding matrix is defined by 3D continuum finite elements. The multi-scale model development is shown in Fig. 9 and a comprehensive explanation about the methodology can be found in [29]. A high-fidelity RVE with CNT 60 vol% and BMI resin has been considered. A simplified geometry of buckypaper mat containing a pair of eight MWCNTs overlapped on another eight MWCNTs has been constructed (Fig. 9a). Inter-wall, intertube, and matrix tube interactions have been modelled by van der Waal's springs (Sect. 2.3). The composite RVE is shown in Fig. 9b. Nanoropes have not been considered here due to a lack of computational resources. The stress-strain non-linearity as per [111] has been defined for the C-C bonds. Material non-linearity in the BMI resin has also been modelled using experimental stress-strain data [138]. The RVE has been clamped at one end, and a tensile deflection has been applied to the opposite end. The reaction force measured at the clamped end has been used to calculate stress. The deformation in the RVE has been used to compute strain. The resulting stress-strain curve is shown in Fig. 9d. The curve is found to be linear up to 78 GPa of stress and 7% of strain. The curve takes a non-linear diversion beyond this stress due to various mechanisms such as non-linearity in resin and CC bonds, and debonding of  LJ springs [29]. The analysis beyond 8% of strain has not been continued due to computational restrictions. The current stress-strain behaviour of the idealized RVE suggests that the strength of the nanocomposite can be at least as high as 78 GPa. Such a tensile strength is far beyond what conventional engineering materials can offer. The current stress-strain behaviour is dominated by the CNT's response to tension as a result of simplified boundary conditions. The other factors such as higher CNT fraction (60 vol%) and tube alignments along the direction of loads are also influencing the strength of RVE. In practice, the boundary conditions and tube alignments can be different from those in the current RVE and it should be modelled accordingly.

Fig. 10
First five mode shapes of the conceptual composite control arm. The corresponding natural frequencies are presented in Table 15 3.6 Analysis of an automotive control arm We first investigate the dynamic behaviour of the composite control arm. The first five vibration mode shapes of the control arm are shown in Fig. 10. The first mode shape (refer to Fig. 10a) is of vertical cantilever bending behaviour that simulates vertical tire motion along the McPherson strut, the second mode shape (refer to Fig. 10b) is an in-plane bending type that simulates wheel acceleration, the third mode shape (refer to Fig. 10c) is a coupled downward vertical motion and in-plane bending, the fourth mode shape (refer to Fig. 10d) is torsional with respect to wheel axis, and the fifth mode shape (refer to Fig. 10e) is an inboard horizontal wheel traversing type. These five mode shapes are found to be similar to the baseline design. It is because the geometry of the design has not been significantly altered while redesigning the baseline into an alternative lightweight conceptual design. An effort has been made to replace materials within the available baseline design volume of the control arm. Furthermore, the distance and the angle between the hardpoints also influence the mode shape. Since the location of hardpoints has not been altered while redesigning the baseline into a conceptual design, the mode shapes are expected to remain the same. However, the natural frequencies become significantly different when the baseline design is compared against the new composite design. Since the mode shapes are similar, it is possible to compare the baseline natural frequencies against the new concept design at a given mode shape.
Here, the natural frequencies of the proposed lightweight concept design will be bench-marked against the baseline design. All the natural frequencies in this section are presented in cycles per second, also referred to as Hertz (Hz). The natural frequencies are mainly dependent on the stiffness and mass matrices of the component. The higher the stiffness of the component, higher will be the natural frequency. Furthermore, the lower the mass, higher will be natural frequency. This means that light-weighting the control arm from baseline aluminium to a composite equivalent will enhance the natural frequency. Also, by enhancing the control arm's stiffness by replacing low modulus aluminium material with high-modulus materials such as buckypaper and carbon fibre prepreg will enhance the natural frequency of the control arm, since the stiffness of any material is derived based on its modulus.
It is important to ensure the orientation of carbon nanotubes and carbon fibres in buckypaper and prepreg with respect to mode shapes to witness natural frequency enhancements. To influence the natural frequencies of all mode shapes, the orientations of buckypaper and prepreg plies have been kept as 0 • , 45 • , 60 • , and 90 • . This ply arrangement offers a certain degree of isotropy (known as in-plane isotropy) to the laminates in the conceptual control arm. The first five natural frequencies of the baseline design and the proposed concept design are shown in Table 15. These natural frequencies are relevant to the mode shapes shown in Fig. 10. The baseline design's natural frequency varies from 365.0 to 3045.0 Hz, for the first five mode shapes. The conceptual control arm's natural frequency varies from 440.0 to 3801.0 Hz, for the first five mode shapes. It is important to ensure that the mode-specific natural frequencies of the chassis body and other components attached to the control arm do not coincide with these frequencies to avoid resonance. The numerical results reveal that the conceptual control arm over-performs as compared to the baseline design. For instance, the first and fifth natural frequencies are enhanced by 21 and 25%, respectively, for the concept design compared to the baseline design. Weight of the control arm drops from 4.2 to 1.8 kg due to light-weighting. This accounts for 57% of weight savings, leading to a considerable reduction in vehicle energy consumption. With the aid of evolving virtual intelligence tools such as topology and ply stacking optimisation, the mechanical performance can be further enhanced along with reduction of weight.
Besides the dynamic characterization, performance of a component under static loading is equally crucial to come up with a successful design. Here, we compare the performance of the lightweight/buckypaper-reinforced control arm against the metallic counterpart, under tensile, compressive, shear, and bending loads. The finite-element method and model implemented here are discussed in Sect. 3.6. A load of 10 kN has been defined on hardpoint 4 to simulate tension, compression, shearing, and bending in four separate load cases. The control arm has been fully constrained at hardpoints 1 and 2. The definition of hardpoints can be seen in Fig. 4. The resulting deflections in the conceptual control arm are shown in Fig. 11. In the tensile test, the tip of the control arm (hardpoint 4) has been pulled in the horizontal direction with the aid of 10 kN force. In the compression test, the tip of the control arm (hardpoint 4) has been pushed in the horizontal direction with the aid of 10 kN force. In the shear test, the tip of the control arm (hardpoint 4) has been displaced sideways with the aid of 10 kN force. In the bending test, the tip of the control arm (hardpoint 4) has been pushed downwards with the aid of 10 kN force. The resulting tensile, compressive, shearing, and bending deflections can be seen in Fig. 11a-d respectively. The total deflections recorded from the finite-element analysis are tabulated in Table 15. This

Conclusions
This article develops a computational bridging of different length scales involving six levels in the range of nano-to macro-scale behaviour of buckypaper composites and the structural components manufactured thereof. The sequential derivatives of carbon at six levels are analyzed involving graphene, CNT, CNT bundle, buckypaper, and buckypaper composites. We have developed a coupled atomistic-continuum modelling approach for the multi-level simulations. Graphene, CNTs, and CNT bundles are modelled using atomistic simulations, while the buckypaper and its composites are modelled using equivalent beam representations for the bundles and continuum solid representation for resin. At the atomic level, we start by establishing the mechanical equivalence of C-C bonds as idealized beams, followed by quantification of the equivalent elastic properties of CNTs and CNT bundles based on a nanoscale finite-element analysis. At the next level, we characterize the mechanical properties of buckypapers and buckypaper composites as programmed networks of CNTs and CNT bundles, which are idealized as equivalent beams with effective properties computed at the preceding length scale level. The simulation results at each level are extensively validated and bench-marked with various configurations considering the data available in the literature. The most salient points concerning the current investigation are listed below • This research paper is the first of its kind to present an automotive structural concept design that uses buckypaper composites as one of the two primary reinforcements. More importantly, the investigation starts at a fundamental nanoscale level and the critical static and dynamic behaviour are characterized at the macro-scale structural level through seamless multi-scale propagation of effective mechanical information. • This article comprehensively captures various levels of effective mechanical behaviour necessary to characterize the macroscopic structural responses. The total number of levels considered here is six, which are bridged through a coupled atomistic-continuum simulation strategy. • The level 1 simulation focuses on characterizing the Young's modulus, shear modulus, and Poisson's ratio of single-layer graphene sheets through an atomistic finiteelement framework. It is noted that for a single-layer graphene sheet, the Young's modulus can be as high as 1082.0 GPa, and the shear modulus can be as high as 606.0 GPa, while the Poisson's ratio comes out to be 0.62. Such results obtained through the current computationally efficient framework agree well with the experi-mental and molecular dynamics studies presented in the literature. • The level 2 simulation explores Young's modulus, shear modulus, and Poisson's ratio of single-walled carbon nanotubes through the atomistic finite-element framework. At this level, Young's modulus can be as high as 1080.0 GPa, and the shear modulus can be as high as 504.0 GPa, while the Poisson's ratio is found to be 0.37. • The level 3 simulation computes Young's modulus, shear modulus, and Poisson's ratio of multi-walled carbon nanotubes considering multiple configurations. At this level, the Young's modulus can be as high as 1401.0 GPa, and the shear modulus can be as high as 638.0 GPa, while the Poisson's ratio remains to be 0.32. • The level 4 simulation computes Young's modulus, shear modulus, and Poisson's ratio of multi-walled carbon nanotube bundles (also known as nanoropes) considering multiple configurations, some of which have hitherto been unexplored. At this level, Young's modulus can be as high as 1461.0 GPa, the shear modulus can be as high as 675.0 GPa, while the Poisson's ratio is noted to be around 0.47. • The level 5 simulation focuses on Young's modulus, shear modulus, and Poisson's ratio of buckypaper composites involving carbon nanotubes and bismaleimide resin. The equivalent properties obtained in the preceding levels are exploited at this level for carrying out a reduced order analysis without compromising the accuracy. Such an efficient analysis framework led to the feasibility of characterizing a wide range of configurations. For buckypaper composites, the Young's modulus can be as high as 370.0 GPa, and the shear modulus can be as high as 27.0 GPa, while the Poisson's ratio comes out to be 0.62. • In the level 6 of simulation, we focus on the full-scale structural analysis of a concept control arm made of buckypaper reinforced composites. Both the dynamic and static performances are compared with traditional monolithic metallic designs based on a finite-element analysis using the effective mechanical properties of buckypaper obtained in level 5. The numerical results demonstrate that a buckypaper composite can enhance the natural frequency and stiffness up to 25 and 37% with respect to conventional monolithic metallic designs, while reducing the weight by 57%. • With the aid of evolving virtual intelligence tools such as topology and ply-stacking optimisation, the mechanical performance of buckypaper composites can be further enhanced along with reduction of weight. Moreover, we note that there is an immense scope of designing high-performance composites according to the application-specific demands of tailored anisotropy through the expanded design space of single-and 1 3 multi-walled nanotube bundles with multiple possible configurations, their shape, relative density, and orientation.
The current article presents an efficient atomistic finite-element-based modelling framework, wherein the contribution is twofold: I. Insightful new results unravelling the mechanical behaviour, II. Development of the generic atomistic finite-element framework for efficient mechanical analysis involving multiple length scales. Adoption of the atomistic finite-element approach instead of conventional simulation methods such as MD simulation brings significantly more computational efficiency in the analysis and the proposed approach is a panacea in situations where inter-atomic potentials required for carrying out MD simulations are unavailable in case of many complex systems. In summary, the numerical outcomes suggest that carbonbased nanostructural derivative in the form of buckypaper can significantly improve the mechanical properties of advanced lightweight structural components as reinforcements for the next generation of aerospace and automotive structures. The characterization of elastic properties at different levels of the derivatives of carbon and their composites along with the proposed multi-scale framework of atomistic-continuum simulations will crucially support further developments in this field. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.