Load-Carrying Capacity of Ultra-Thin Shells with and without CNTs Reinforcement

: Isotropic ultra-thin shells or membranes, as well as cable–membrane structures, cannot resist loads at the initial state and always require a form-ﬁnding process to reach the steady state. After this stage, they can work in a pure membrane state and quickly experience large deﬂection behavior, even with a small amplitude of load. This paper aims to improve the load-carrying capacity and strength of membrane structures via exploiting the advantages of functionally graded carbon-nanotube-reinforced composite (FG-CNTRC) material. In this work, the load-carrying capacity and nonlinear behavior of membrane structures with and without CNTs reinforcement are ﬁrst investigated using a uniﬁed adaptive approach (UAA). As an advantage of UAA, both form ﬁnding and postbuckling analysis are performed conveniently and simultaneously based on a modiﬁed Riks method. Different from the classical membrane theory, the present theory (ﬁrst-order shear deformation theory) simultaneously takes into account the membrane, shear and bending strains/stiffnesses of structures. Accordingly, the present formulation can be applied adaptively and naturally to various types of FG-CNTRC structures: plates, shells and membranes. A veriﬁcation study is conducted to show the high accuracy of the present approach and formulation. Effects of CNTs distribution, volume fraction, thickness, curvature, radius-to-thickness and length-to-radius ratios on the form-ﬁnding and postbuckling behavior of FG-CNTRC membranes are particularly investigated. In particular, equilibrium paths of FG-CNTRC membrane structures are ﬁrst provided in this paper.


Introduction
Ultra-thin shells or membranes are found in many engineering fields, from terrestrial applications to space applications, as air-bags, roofs of buildings, fabric structures, exhibition pavilions, parachutes, inflatable reflectors, antennae and solar power collection panels on spacecraft, space-based radars, etc. [1]. Membrane structures are widely used due to their high strength-to-weight ratio, light weight and flexibility. Accordingly, they can be applied flexibly to short-span, long-span and free-form shape buildings, and their behaviors should be particularly investigated. In mechanical models of isotropic membrane structures, due to their thicknesses being very small, it is assumed that the membranes have zero bending stiffness and pure membrane stiffness. Isotropic membrane structures cannot resist loads at the initial or intermediate states as illustrated in Figure 1. They always require a form-finding process to reach the steady state. After this state, the isotropic membranes can resist loads, work in a pure membrane state and quickly experience large deflection, even with a small load. Due to the importance of the form-finding process, many numerical techniques have been developed for form finding, such as: dynamic relaxation [2], force density method [3], a nonlinear solver and updated reference strategy [4], the natural-shape-finding method [5], the direct and indirect approach of the force density method [6,7], the minimal surface method [8], etc.

Material models
The coated fabric considered in this work is the polyes woven fabric with PVC coating, which is manufactured b Industry, Belgium. This kind of coated fabric is ubiquitou now in tension membrane structures. The present autho proposed a material constitutive law, viz. , fabric-plasticity for this material, details of this model can be found in Di (2014) . This material model has shown a good performance turing the nonlinear behavior of the coated fabric when w it for a load analysis of a large-scale membrane ( Dinh et a Therefore, in this paper, we also use it to simulate the me behavior of the fabric membrane. Moreover, in practice, th orthotropic elasticity model is widely used. Thus, for the comparison this material model is also used for the fabri simulation. Comparison of the obtained results when we ferent material models for the membrane, i.e. , fabric-plasti orthotropic elasticity, will bring us insight into the amoun ric compensation that designers should apply to their fabri when the elasticity model is used in their computation. In attempts to improve the load-carrying capacity and strength of isotropic membrane structures, they are often stiffened by cables for reinforcement of strength. However, cable-membrane composite structures still exhibit some complicated problems, such as: (1) the folding of the membrane by the cable; (2) the sliding of the cable on the membrane surface [10]; (3) cable-membrane structures can only maintain their shapes with in-plane tensile stress [11] and so on. In addition, a form-finding process is still required even when membranes are stiffened by cables [11]. The above problems can lead to some numerical difficulties in modeling cable-membrane composite structures. To overcome the existing difficulties and improve the load-carrying capacity, as well as the strength of membrane structures, functionally graded carbon-nanotube-reinforced composite membrane structures are first proposed and investigated in this paper using a unified adaptive approach (UAA) in the framework of isogeometric analysis (IGA) and the first-order shear deformation theory (FSDT). As an advantage of UAA, both form finding and analysis are performed conveniently and simultaneously based on a modified Riks method. The present structure exploits the advantages of functionally graded carbon-nanotube-reinforced composite (FG-CNTRC) material, such as: a high stiffness, low density, high strength and aspect ratio [12][13][14]. It is interesting that the present approach and formulation can be applied adaptively and naturally to FG-CNTRC plates/shells and membranes from ultra-thin to moderate thick structures without concerning the global stiffnesses of structures. In attempts to develop advanced methods, isogeometric analysis [15] was proposed to solve mechanics problems with exact geometries and high accuracy. IGA has been developed for various types of structures, such as plates [16][17][18][19][20][21][22] and shells [23][24][25][26][27][28][29]. In particular, NURBS-based nonlinear inverse analyses of thin shells were performed not only for pure mechanical problems [30] but also for coupled problems [31]. In addition, IGA has been successfully investigated for isotropic membrane structures, such as: classical shell theory for the analysis of membranes [32], design and analysis [33], vibration analysis [34], shape optimization of vibrating membranes [35], transient analysis [36], form finding with dynamic relaxation [2] and form finding and large deflection analysis [37]. In this paper, the form finding and postbuckling analysis of FG-CNTRC membrane structures are first carried out with the goal of improving the load-carrying capacity and strength of membrane structures. As the developments of IGA, IGAs using PHT-splines [38] and RHT-splines [39] were proposed to drive the adaptive refinement and improve the computational efficiency. In addition, combinations of artificial intelligence [40,41] and computational mechanics have been proposed to analyze structures such as: deep collocation methods [42][43][44], deep energy methods [45] and the deep learned one-iteration method [46].
The paper is organized as follows: a postbuckling analysis of functionally graded carbon-nanotube-reinforced composite membranes is shown in the next section. Section 3 presents the form finding and postbuckling analysis of FG-CNTRC membrane structures using a unified adaptive approach. Results and discussions are adequately provided in Section 4. The paper is closed with several interesting main conclusions in Section 5.

Functionally Graded Carbon-Nanotube-Reinforced Composite Membranes
All of the FG-CNTRC membrane structures in this study are investigated with four types of CNT distributions as described in Figure 2. Four corresponding CNT volume fractions are determined as follows: and It is assumed that the CNTs mass fraction w CNT is equal for the four types of FG-CNTRC membranes. In addition, ρ CNT and ρ m , respectively, are the densities of CNTs and the matrix, and h is defined as the thickness of membranes. In this study, we employ the rule of mixture [47,48] to estimate the effective properties of CNT-reinforced materials as where G CNT 12 is the shear modulus and E CNT 11 , E CNT 22 are Young's moduli of CNTs. G m is the shear modulus and E m is Young's modulus of the isotropic matrix material. In addition, V m and V CNT , respectively, denote the volume fractions of the matrix material and the CNTs. The sum of these volume fractions is one. The efficient parameters η i (i = 1, 2, 3) are considered to account for the scale-dependent material properties of CNTs. Poisson's ratio v 12 and the density ρ can be determined as follows: where v CNT 12 and v m , respectively, are Poisson's ratios of the CNTs and the matrix material. Note that v 12 is constant along the membrane thickness.

Postbuckling Analysis of FG-CNTRC Membrane Structures Using FSDT
Different from the classical membrane theory, the present theory (FSDT) simultaneously takes into account the membrane, shear and bending strains/stiffnesses of structures. Hence, the present formulation can be used adaptively and naturally for various types of FG-CNTRC structures-plates, shells and membranes-without concerning the global stiffnesses of structures. Firstly, we investigate a FG-CNTRC cylindrical membrane structure with the initial geometry given in Figure 3. The displacements of an arbitrary point in the membrane based on FSDT are defined as [49][50][51].
u(x, y, z) = u 0 (x, y) + zβ x (x, y) v(x, y, z) = v 0 (x, y) + zβ y (x, y) orū = u 0 + zu 1 w(x, y, z) = w 0 (x, y) where u 0 , v 0 are the tangential displacement components and w 0 is the radial displacement component or the deflection of the considering point. β x and β y , respectively, are the rotations in the y and the x axes. The strains in the orthogonal curvilinear coordinate system (oxyz) with the von Karman assumption are defined as [51] ε = ε xx ε yy γ xy and The nonlinear strain component is re-expressed as The virtual work equation of the FG-CNTRC membrane structure is written as follows: where the load vector For the postbuckling analysis of the FG-CNTRC membrane structure subjected to a radial pressure f z = λ f 0 , the virtual work equation is conveniently re-written as where λ andσ, respectively, denote the load factor and the stress resultant vector. The vectorσ is defined asσ This vector has three components: In plane bending shear Hooke's law is used to describe the relationship between the generalized strain vector ε and the stress resultant vectorσ aŝ In this study, we employ the shear correction factor κ = 5/6 [49] and

A Brief Introduction to the Unified Adaptive Approach (UAA) for Form Finding and Analysis of Membrane Structures
For the analysis of an isotropic membrane, a form-finding process is first carried out to reach the steady state or to achieve the new geometry of the membrane. Next, the analysis step is conducted based on this new geometry. In this step, the isotropic membrane works in a pure membrane state and can experience large deflection quickly, even with a very small load. For the analysis of membrane structures using the conventional numerical approach, form finding and analysis are separately carried out, including a form-finding process and numerical techniques for the geometrically nonlinear analysis of membranes. In this study, the form finding and analysis of FG-CNTRC membrane structures are conducted at the same time using the unified adaptive approach (UAA), which exploits the advantage of the modified Riks nonlinear solver [37]. Note that, when employing the nonlinear solver for form finding, the shape-generating load is often a surface load, such as: snow, self-weight [4], etc. In this study, the surface load is employed for the two processes: form finding and analysis. For isotropic membrane structures, a new type of equilibrium path at the central point was proposed based on UAA. This path includes two parts: one is achieved via form finding and the rest is from the analysis [37].

A Brief Introduction to Non-Uniform Rational B-Splines Basis Functions
In this subsection, we present a summarization of non-uniform rational B-spline (NURBS) basis functions. The details, as well as source codes, of isogeometric analysis using NURBS basis functions can be found in [52]. We first investigate a knot vector This vector is open if the last and first knots are repeated p + 1 times. Note that a B-spline basis function has the C p−1 continuity at each knot while it has the C ∞ continuity inside a knot span. The B-spline basis functions in the one-dimensional parametric space N i,p (ξ) are determined as follows [53]: We next investigate two knot vectors in the two-dimensional parametric space Ξ = ξ 1 , ξ 2 , . . . , ξ n+p+1 and H = η 1 , η 2 , . . . , η m+q+1 . The B-spline basis functions in this space can be determined via a tensor product as The NURBS basis functions can be employed to exactly describe geometries and are determined by assigning the weight w for each control point as follows [15]: Note that the NURBS and the B-spline basis functions are completely the same if all of the weights of the control points are identical.

Form Finding and Postbuckling Analysis of FG-CNTRC Membrane Structures Using UAA and NURBS
Applying the unified adaptive approach (UAA) [37], both form finding and postbuckling analysis are carried out at the same time via a modified Riks nonlinear solver as mentioned earlier. In the present formulation, NURBS basis functions are employed to interpolate geometries and displacements of membranes as follows: N A ξ, η denotes the NURBS basis function and P A stands for the control point coor- involves degrees of freedom of the typical control point A. In addition, the vector x h x, y contains the physical point coordinates and u h denotes the interpolated displacement vector. Combining Equations (11) and (25), the strain vectors can be interpolated via the control point displacements as and Substituting Equation (25) into Equation (12), the interpolated nonlinear strain vector is achieved as follows: We substitute the interpolated strains from Equations (26) and (28) into Equation (19) and achieve the variation of the generalized strain vectorε as Next, the terms in Equations (25) and (29) are substituted into Equation (14). We then eliminate the virtual displacement component. The nonlinear equation can be written conveniently as follows: In this study, the nonlinear equation in (30) is iteratively solved by employing the modified Riks method [54]. When a change in the load from λ m f 0 to (λ m + ∆λ m ) f 0 is created, a new equilibrium configuration close to the old one is reached and expressed as Only taking into account the first term after applying the Taylor series expansion to Equation (31), we achieve a system of the linear incremental equilibrium equations as where K T stands for the tangent stiffness matrix determined at q m as We next carry out an iterative process for all load increments [54]. Conveniently, Equations (32) and (33) can be re-expressed in a generalized form with the mth load increment and the ith iteration as with and the load vector where F 0 denotes the referenced load vector, while the internal load vector is determined as We substitute all of the terms in Equations (19), (26), (28) and (29) into Equation (37) and obtain the following equation: with The iterative computation is carried out for all load increments. To check the convergence of solutions, the criterion based on the residual load vector is applied as The incremental solutions are determined from Equation (34). We then update the load factor and the displacement vector of each iteration as [54] where ∆ i q R,m and q F,m , respectively, are the displacement vectors induced by the residual and referenced load vectors as

Verification Study
This section aims to confirm the accuracy of the present formulation and approach for the form finding and postbuckling analysis of FG-CNTRC membrane structures. It should be noted that a high accuracy of the present approach for the form finding and analysis of isotropic membrane structures was confirmed in [37]. In addition, this is the first study on the analysis of FG-CNTRC membrane structures, and there is no reference solution in the literature to compare to. Accordingly, the proposed formulation and approach are verified via some benchmarks of isotropic FG-CNTRC cylindrical shell structures and a composite circular cylinder shell. In this work, the central deflection is positive, corresponding to the central point, which is deflected downward. Inversely, the central deflection is negative, corresponding to the central point, which is deflected upward. The load is downward pressure as illustrated in Figure 3. In this study, the surface load with a small load step ∆q = 10 −5 N/m 2 is employed as the shape-generating load. This load is used for the two processes: form finding and analysis. For the convenience and efficiency, all of the problems in this study are analyzed employing the mesh of 14 × 14 cubic elements and 4 × 4 Gauss points per each element. The armchair (10, 10) SWCNTs [55] and Poly methyl methacrylate (PMMA) [56] are, respectively, employed as the reinforcements and the matrix material for FG-CNTRC membranes. In addition, all of the membranes are circumferentially reinforced with CNTs. Material properties at room temperature are employed as For the CNTs reinforcement material, the efficiency parameters corresponding to the volume fractions are given as [55] • V * CNT = 0.11: η 1 = 0.149 and η 2 = 0.934; • V * CNT = 0.14: η 1 = 0.150 and η 2 = 0.941; • V * CNT = 0.17: η 1 = 0.149 and η 2 = 1.381.
We can assume that G 12 = G 13 = G 23 and η 3 = η 2 . For an FG-CNTRC membrane with the initial geometry as described in Figure 3, the clamped boundary is directly imposed as In this study, two non-dimensional parameters are defined and investigated as follows: • Normalized central deflection:ŵ = w/h; • Curvature: κ = L/R. In this section, geometrically nonlinear analyses of isotropic and FG-CNTRC cylindrical shell structures are carried out to verify the present formulation and approach. Firstly, we investigate a fully clamped isotropic cylindrical panel subjected to a pressure. The geometry of the panel is determined as: the thickness h = 3.175 mm, L = B = 508 mm and the radius R = 2540 mm. The characteristics of the material are Poisson's ratio v = 0.3 and Young's modulus E = 3.103 kN/mm 2 . The obtained result is presented in Figure 4 in a comparison with the result of Crisfield [57] using the finite element method. A very good agreement is found. Next, a fully clamped FG-CNTRC cylindrical panel subjected to a pressure, as shown in Figure 3, is investigated. The geometry of the panel is described as follows: α = 0.1 rad, h/R = 0.004 and L/R = 0.1. The obtained results are shown in Figure 5 in a comparison with the results of the mesh-free kp-Ritz method (MKR) [58]. A very good agreement is obtained for all types of CNTs distributions. It can be seen that the present approach and formulation have a high accuracy for nonlinear analyses of isotropic and FG-CNTRC cylindrical shell structures. Finally, we can conclude that the present approach and formulation have a high accuracy for the form finding and postbuckling analysis of FG-CNTRC membrane structures.

FG-CNTRC Cylindrical Membranes
In this section, we consider a fully clamped cylindrical membrane structure under a pressure. The geometric characteristics of the cylindrical membrane (C-M) are the thickness h = 0.5 mm, L = B = 508 mm and the radius R = 2540 mm. First, we investigate the behavior of the cylindrical membrane with the isotropic material as Poisson's ratio v = 0.3 and Young's modulus E = 3.103 kN/mm 2 . The obtained result is presented in Figure 6. It is interesting that, although the loading is downward, the central point is deflected upward at the beginning. This phenomenon occurs due to the form-finding process as explained in [37], published by the authors of this paper. It is observed that the structure works in a pure membrane state and quickly experiences large deflection behavior, even with a small pressure, as seen in Figure 6. It should be noted that the results in Figures 4 and 6 are obtained from two models with the same material, geometric data, loading and boundary condition. The only difference is the thickness (h = 3.175 mm and h = 0.5 mm). It is clear that behaviors of the isotropic membrane and thin-shell structures are completely different. Next, we investigate the behavior of the above cylindrical membrane structure (C-M) with the FG-CNTRC material. The obtained results are depicted in Figure 7. It is interesting that, when using the FG-CNTRC material, the membrane structure works as a bending structure. It is clear that behaviors of the isotropic and FG-CNTRC membrane structures are completely different as seen in Figures 6 and 7. The FG-CNTRC membrane structure can resist the load at the beginning without finding a new geometry, which usually occurs in any isotropic membrane structures (the so-called form-finding procedure). Using the same geometric data, loading and boundary condition, the FG-CNTRC membrane structure has a significantly higher nonlinear strength compared with the isotropic membrane structure as seen in Figures 6 and 7. It can be concluded that using the FG-CNTRC material for membrane structures is very efficient. It is observed that FG-X is the best distribution and the one that gives the greatest postbuckling strength. Figure 8 shows the effect of CNTs volume fraction on the nonlinear response of the FG-CNTRC cylindrical membrane structure. The larger the amount of CNTs in the membrane, the higher the postbuckling strength. Without CNTs reinforcement (V * CNT = 0), the membrane cannot resist the load at the initial state and requires a form-finding process to reach the steady state. Conversely, with CNTs reinforcement, the membranes can resist the load at the beginning. It is concluded that CNTs reinforcement can create a load-carrying capacity for cylindrical membrane structures, even at the initial state. Figures 7, 9 and 10 show the effect of CNTs distribution on the nonlinear response of the FG-CNTRC membrane with various curvatures κ = 0.2, 0.04, 0. It is observed that the smaller the curvature of the membrane structure, the lower the postbuckling strength. It is clear that using the FG-CNTRC material for the cylindrical membrane structure is very effective, whereas, for the square membrane structure, it is not effective. As seen in Figures 7 and 10, although various types of CNT distributions are investigated, the nonlinear strength of the square membrane structure (κ = 0) does not change, whereas that of the cylindrical membrane structure is significantly improved. As mentioned earlier, FG-X is the best distribution and the one that gives the greatest postbuckling strength. Accordingly, using the FG-CNTRC material for cylindrical membrane structures is very effective and recommended, whereas using it for square membrane structures is not recommended. It is interesting to note that, when the curvature of the membrane is very small-κ = 0.04 (an almost flat membrane)-the membrane still works as a bending structure. However, when the curvature κ = 0, the square membrane works as a pure tensile structure. This phenomenon occurs due to the fact that CNTs reinforcement in the circumferential direction can significantly improve the nonlinear bending strength of the membrane structure, even when the curvature is small. In addition, Figure 11 shows a strong effect of the curvature on the nonlinear response and behavior of the FG-CNTRC membrane, although the change in curvature is very small, as described in Figure 12. As seen in Figure 11, the membrane works as a bending structure even when the curvature of the membrane is very small but not equal to zero. However, when the curvature κ = 0, the square membrane suddenly works as a pure tensile structure. It is noted that the classical membrane theory cannot model FG-CNTRC membranes correctly due to the present of bending stiffness created by CNTs reinforcement. Different from the classical membrane theory, the present theory (first-order shear deformation theory) simultaneously takes into account the membrane, shear and bending strains/stiffnesses of structures. Accordingly, the present formulation can be applied adaptively and naturally to various types of behaviors and structures, such as: FG-CNTRC plates/shells and membranes. The high efficiency of the present formulation is confirmed. The effect of the thickness on the nonlinear response of the FG-CNTRC cylindrical membrane is very significant as seen in Figure 13. It can be concluded that, for an FG-CNTRC cylindrical membrane structure, an ultra-small change in thickness or curvature can lead to a significant change in its nonlinear response, strength and behavior.

FG-CNTRC Circular Cylinder Membranes
We next analyze a glass-epoxy circular cylinder shell with an internal pressure as illustrated in Figure 14. The cylinder is fully clamped at two ends. The properties of the geometry are given as follows: the thickness h = 1 in., the radius R = 20 in. and the length L = 20 in. The material data are provided as follows: E 1 = 7500 ksi, E 2 = 2000 ksi, v 12 = 0.25, G 12 = 1250 ksi and G 13 = G 23 = 625 ksi. We exploit the symmetry of the cylinder shell: one eighth (ABCD) of the cylinder is analyzed without a loss of accuracy of the solution. The achieved result is presented in Figure 15 in a comparison with the result of Kreja [59] employing FEM. The present solution has a good agreement with the referenced one. It can be concluded that the present formulation and approach have a high reliability in predicting the geometrically nonlinear behavior of orthotropic cylinder shells. Next, an adequate investigation into nonlinear behaviors of FG-CNTRC circular cylinder shells and membranes is first performed in this study. Fully clamped FG-CNTRC cylinder shells and membranes with a fixed thickness h = 1 in. under external pressure are considered. Figure 16 shows the significant effects of CNTs distribution on nonlinear responses of the FG-CNTRC cylinder shell and membrane. For the cylinder shells, instabilities including "snap-through" or "snap-back" occur depending on the CNTs distribution. It is interesting to note that the cylinder membranes are stable during loading and quickly experience large deflection behavior, even with a very small amplitude of pressure. For both FG-CNTRC cylinder shells and membranes, FG-O is the worst distribution and the one that gives the lowest nonlinear strength, and FG-X is the best distribution and the one that gives the greatest nonlinear strength. Figure 17 shows the effect of the CNTs volume fraction on the nonlinear response of the FG-CNTRC cylinder membrane. It is clear that, the greater the amount of CNTs in the membrane, the greater the nonlinear strength. The effects of the radius-to-thickness ratio on nonlinear responses of the FG-CNTRC cylinder shell and membrane are described in Figure 18. For the cylinder shells (R/h = 60∼100), instabilities can occur, whereas the cylinder membranes (R/h = 1000∼2200) are stable and quickly experience large deflection behavior. The effect of the radius-to-thickness ratio on the nonlinear response of the cylinder shell is very significant, whereas there is no effect of this ratio on the response of the cylinder membrane. Finally, Figure 19 shows the effects of the length-to-radius ratio on nonlinear responses of the FG-CNTRC cylinder shell and membrane. The effect of the length-to-radius ratio on the instability and nonlinear response of the FG-CNTRC cylinder shell is significant. However, this ratio almost does not affect the instability and nonlinear response of the FG-CNTRC cylinder membrane.

Conclusions
The load-carrying capacity and nonlinear behavior of ultra-thin shells or membranes with and without CNTs reinforcement have been investigated for the first time using a unified adaptive approach (UAA). As an advantage of UAA, both form finding and a postbuckling analysis are performed conveniently and simultaneously based on a modified Riks method. Different from the classical membrane theory, the present theory (FSDT) simultaneously takes into account the membrane, shear and bending strains/stiffnesses of structures. Accordingly, the present formulation can be applied adaptively and naturally to various types of FG-CNTRC structures-plates, shells and membranes-without any restriction or modification. The high accuracy of the present approach and formulation was verified via the benchmark problems. In particular, equilibrium paths of FG-CNTRC membrane structures have been provided for the first time in this paper, which could be useful for future references. From the parametric study, interesting major conclusions can be drawn: • Isotropic membrane structures always require a form-finding process to reach the steady state and then work in a pure membrane state under the load. However, FG-CNTRC cylindrical membrane structures work as bending structures and can resist loading at the beginning without needing the form-finding process. Accordingly, using the FG-CNTRC material can create a load-carrying capacity for cylindrical membrane structures, even at the initial state. • For a FG-CNTRC cylindrical membrane structure, an ultra-small change in thickness can lead to a significant change in its nonlinear response and strength. • Using the FG-CNTRC material for cylindrical membrane structures is very effective and recommended, whereas, for flat membrane structures, it is not recommended. • A FG-CNTRC cylindrical membrane works as a bending structure, even when the curvature of the membrane is very small but not equal to zero. However, when the curvature κ = 0, the square membrane suddenly works as a pure tensile structure. • Under external pressure, for FG-CNTRC circular cylinder shells, types of instabilities with various critical buckling loads can occur, whereas FG-CNTRC circular cylinder membranes are stable and quickly experience large deflection behavior, even with a very small amplitude of pressure. • Radius-to-thickness and length-to-radius ratios significantly affect the instability and nonlinear response of FG-CNTRC circular cylinder shells, but almost do not affect those of FG-CNTRC circular cylinder membranes. Meanwhile, the CNTs distribution and volume fraction significantly affect nonlinear responses of both FG-CNTRC circular cylinder shells/membranes. • FG-X is the best distribution and the one that gives the greatest nonlinear strength of a FG-CNTRC membrane structure. The greater the amount of CNTs in the membrane, the greater the nonlinear strength.
The present work has some novelties, obtains new results using IGA and significantly improves the load-carrying capabilities of ultra-thin structures via using CNTs reinforcement. However, this work has some limitations that we need to overcome in future works, which are as follows: (1) the present formulation is applicable to plates or shells with one curvature. It should be developed and improved to model free-form shells; (2) in this work, the rule of mixture is used to calculate the effective properties of CNT-reinforced materials. The rule of mixture does not consider the interface behavior between the reinforcement and the matrix material. In future, we should use some material models that can consider this interface behavior.