Abstract
Functionally graded material (FGM) arches may be subjected to a locally radial load and have different material distributions leading to different nonlinear in-plane buckling behavior. Little studies is presented about effects of the type of material distributions on the nonlinear in-plane buckling of FGM arches under a locally radial load in the literature insofar. This paper focuses on investigating the nonlinear in-plane buckling behavior of fixed FGM arches under a locally uniformly distributed radial load and incorporating effects of the type of material distributions. New theoretical solutions for the limit point buckling load and bifurcation buckling loads and nonlinear equilibrium path of the fixed FGM arches under a locally uniformly distributed radial load that are subjected to three different types of material distributions are derived. The comparisons between theoretical and ANSYS results indicate that the theoretical solutions are accurate. In addition, the critical modified geometric slendernesses of FGM arches related to the switches of buckling modes are also derived. It is found that the type of material distributions of the fixed FGM arches affects the limit point buckling loads and bifurcation buckling loads as well as the nonlinear equilibrium path significantly. It is also found that the limit point buckling load and bifurcation buckling load increase with an increase of the modified geometric slenderness, the localized parameter and the proportional coefficient of homogeneous ceramic layer as well as a decrease of the power-law index p of material distributions of the FGM arches.
Introduction
Arch structure has been widely used in the practice, because of the unique characteristics about aesthetics and safety. Arch structure is often subjected to different forms of load and different boundary constraints in practical engineering, which may lead to structural damage of arches. (Timoshenko et al., 1962), (Schreyer, 1972) and (Plaut and Raymond, 1990) started the research of arches early. Then, (Bradford et al., 2002) researched the in-plane buckling behavior of symmetrical cross-section arches under a central point load. (Pi et al., 2002) used an energy method to build the nonlinear equilibrium equation and the buckling equilibrium equation of arches, and the analytical result of nonlinear buckling was obtained. (Pi et al., 2008) also presented a theoretical research on the nonlinear in-plane buckling of pin-ended circular arches subjected to a central point load with rotational restraints elastic end. (Pi and Bradford, 2009) developed the virtual work approach was used to establish post-buckling equilibrium differential equation and the nonlinear buckling, and the exact solution of the nonlinear bifurcations. (Cai and Feng, 2010) carried out the in-plane buckling of parabolic arches whose rotational stiffness of supports increases according to loaded. (Cai et al., 2012) carried out the buckling equations and nonlinear equilibrium by virtual work principal method. (Cai et al., 2013) also analyzed the buckling equilibrium path according to the bistable strut by using fixed arches subjected to a central point load. According to the principle of virtual work, (Han et al., 2016) studied the in-plane nonlinear buckling behavior of circular arches with elastic horizontal supports under a uniform radial load. (Yan et al., 2017) made an analytical study about the non-uniform shallow arch subjected to a central point load. (Pi et al., 2017) revealed in-plane buckling of fixed shallow circular arches with the arbitrary radial point load. (Lu et al., 2018) explored the nonlinear in-plane buckling and post-buckling behavior of the fixed arches subjected to a localized uniform radial load. (Lu et al., 2020) also explored effects of movement and rotation of supports on nonlinear instability of fixed shallow arches under a localized uniform radial load. (Hu et al., 2021) presented an analytical investigates for nonlinear buckling of pin-ended arch-beam structures. In addition to the optimization of the cross section can make the stability of the structure better (Chen and Young, 2021; Chen et al., 2021), the improvement of material engineering technology also can make the stability of the structure better.
Most of the arches studied in the above literatures are homogeneous solid materials, however, functionally graded material (FGM) has been widely used because their mechanical properties are better than those of homogeneous materials in recent years. FGM arches are usually composed of two or more materials, which can contribute to the superior performance of each material, but the stress and buckling of FGM arches become complicated. (Bateni and Eslami, 2015) conducted nonlinear buckling behavior of FG circular arches under the uniformly radial load using an analytical method. (Simsek, 2016) assumed the material properties of the beam vary as the power-law form in both axial and thickness directions. (Al-shujairi and Mollamahmutoğlu, 2018) concerned with the sandwich FG micro-beams is made of the ceramic-metal FG skin and homogenous ceramic core, based on the classical rule and the Mori-Tanaka homogenization scheme of composite material properties of the part of the FG sandwich size according to beam changes continuously by the thickness of the beam. For different boundary condition analysis, (Bateni and Eslami, 2014) explored the nonlinear in-plane buckling behavior of FGMs arches. Simply supported–simply supported and clamped–clamped boundary conditions are considered as two types of well-known symmetric boundary conditions for this analysis. (Simsek, 2016) investigated first two-dimensional FGM to research the buckling of beams with different boundary conditions, using adding auxiliary functions the boundary conditions are met. For different load form analysis, (Rastgo et al., 2005) explored the spatial buckling of pinned FGM arches being subjected to a thermal loading for using the Galerkin method. (Song and Li, 2008) using Kirchhoff” s assumption conducted nonlinear governing equations of FG arches under thermal and mechanical loads. (Asgari et al., 2014) focused the nonlinear thermo-elastic behavior of pin-ended FGM arches. The stability equation of the nonlinear primary equilibrium path was developed using adjacent equilibrium criterion. (Moita et al., 2018) studied the discrete model of the structural by the theory of high order shear deformation, and the finite element formulation is carried out to meet general FGM plate-shell type structures. (Yang et al., 2020a; Yang et al., 2020b; Yang et al., 2020c; Yang et al., 2021a; Yang et al., 2021b) have published a series of papers on in-plane buckling and free vibrations of functionally graded composite arches with graphene reinforcements. However, FGM arches may be subjected to a locally radial load and have different material distributions leading to different nonlinear in-plane buckling behavior. Little studies is presented about effects of the type of material distributions on the nonlinear in-plane buckling of FGM arches under a locally radial load in the literature.
This paper, therefore, aims at studying the in-plane buckling of a fixed functionally graded material arch under a locally uniformly distributed radial load incorporating effects of the type of material distributions. To illustrate the nonlinear in-plane buckling behavior of FGM arches, three types of material distributions containing functionally graded material layer and homogeneous material layer are considered. Theoretical solutions for the limit point buckling load and bifurcation buckling loads and nonlinear equilibrium path of the fixed FGM arches under a locally uniformly distributed radial load that are subjected to three different types of material distributions are derived. Comparisons with the ANSYS results present the theoretical solutions are enough accurate. The critical modified geometric slendernesses of FGM arches that are switches of buckling modes are also derived. In addition, the effects of the proportional coefficient of homogeneous ceramic layer, the type and power-law index p of material distributions on the buckling load and critical modified geometric slendernesses are discussed and examined.
Material Distribution of Functionally Graded material Arches
The FGM shallow arches under a locally uniformly distributed radial load q with a central angle , and radius R and length S are made from ceramic Al2O3 and metal Al, as shown in Figure 1. In addition, 2c represents the subtended angle of the arch segment where the locally uniformly distributed radial load acts, also as shown in Figure 1. The total cross-sectional thickness and width are h and b, respectively.
FIGURE 1
In this research, the metal-ceramic material Al-Al2O3 is adopted. According to power-law function, the FGM arch through the thickness are graded continuously aswhere and are the properties of ceramic and metal, respectively. and is the volume fraction of ceramic and metal, and they meet the following conditions as
In order to illustrate the nonlinear in-plane buckling behavior of FGM arches, three types of material distributions containing functionally graded material layer and homogeneous material layer are considered, including Type A, Type B, and Type C as shown in Figure 2. The cross-section of Type A represents that the top layer is the metal layer with the thickness and the material distributions of bottom layer is that the weight fraction of metal increase from the top surface to the bottom and the weight fraction of ceramic decrease from the top surface to the bottom surface. The cross-section of Type B have three layers, the middle layer of Type B is the metal layer with the thickness , the top and bottom layers of Type B is symmetric about the midline of the cross-section, and the material distributions of bottom layer of Type B is same to the bottom layer of Type A. The cross-section of Type C have four layers, the top and the bottom layers of Type C is the metal layer with the thickness , the second and third layers of Type B is symmetric about the midline of the cross-section, and the material distributions of the second layer of Type C is same to the bottom layer of Type A.
FIGURE 2
The volume fractions for three distribution modes that are shown in Figure 2 are expressed asandwhere p is the power-law index of material distributions and is the proportional coefficient of homogeneous ceramic layer. According to Eq. 1, the elastic modulus of FGM arch can be expressed as (Bateni and Eslami, 2015).where and are the elastic modulus of ceramic and metal, respectively.
Nonlinear Equilibrium
The nonlinear in-plane buckling behavior of a fixed FGM arch subjected to a locally uniformly distributed radial load is investigated, which is shown in Figure 1. The potential energy Π of the FGM arch can then be given bywith the function being given bywhere represents Heaviside function. The longitudinal strain of an arbitrary point at the cross-section can be expressed aswhere and are the membrane and bending strains of the FGM arch, respectively, which are given bywith and , where , z is the coordinate of the point P0 at the cross-sectional axis oz, w and v are the axial and radial displacements respectively.
For nonlinear equilibrium of FGM arch using the principle of minimum total potential energy, the variation of the total potential energy obtained from Eq. 7 should vanish as
By substituting Eqs 10, 11 into Eq. 12, the Eq. 12 can then be rewritten aswhere N and M are the axial and bending actions of the FGM arch, respectively, which can be given bywith , , and being the stiffness components of the FGM arch, respectively, which are defined as
When material distributions of the FGM arch is subjected to Type A, the stiffness components of the FGM arch , , and can be given by
When material distributions of the FGM arch is subjected to Type B, the stiffness components of the FGM arch , , and can be given by
When material distributions of the FGM arch is subjected to Type C, the stiffness components of the FGM arch , , and can be given by
Integrating Eq. 13 by parts, the differential equilibrium equations of a FGM arch under a locally uniformly distributed radial load can be derived as
By Substituting Eq. 14 into Eq. 15, the bending action of the FGM arch can be given by
By Substituting Eqs. 26, 28 into Eq. 27, Eq. 27 can then be rewritten aswhere is non-dimensional axial force parameter, which is given bywith being the equivalent stiffness of the FGM arch, which is defined as
Meanwhile, the essential boundary conditions at the both ends of FGM archalso need to be satisfied. The non-dimensional radial displacement of a FGM arch can be determined by solving Eq. 29 under the boundary conditions obtained from Eq. 32.where and with , the non-dimensional locally uniformly distributed radial load P is expressed asand the is expressed as
Equation 33 indicates that the non-dimensional radial displacement is dependent on P and , and so it is necessary to calculate the relationship of the axial force with the locally uniformly distributed radial load.
Therefore, substituting Eq. 33 into Eq. 14, and integrating it along the arch length results in
Considering the essential boundary conditions obtained from Eq. 32, the Eq. 36 can be rewritten as
Subsequently, Calculating Eq. 37 leads to a quadratic equilibrium equation, aswhere the coefficients , and are expressed aswith λs being the modified geometric slenderness of the FGM arch, respectively, which can be express as
By solving Eq. 38, the relationship of the axial force with the locally uniformly distributed radial load can be obtained. Then, by substituting those relationship into Eq. 33, the non-dimensional radial displacement can be obtained. Finally, the relationship of the locally uniformly distributed radial load with the non-dimensional radial displacement can be established, and the nonlinear instability behavior of FGM arches can be analyzed.
Buckling Analysis
Limit Point Buckling and Bifurcation Buckling
The FGM arches under the locally uniformly distributed radial load may buckle in a limit point mode. According to the theory of (Lu et al., 2018) for an arches under the locally uniformly distributed radial load, the upper and lower limit point buckling loads can be derived using routine calculus, and the equation of equilibrium between and P at the upper and lower limit points can be expressed aswhere
Therefore, the upper and lower limit point buckling loads can be obtained by solving Eqs 38, 43 simultaneously.
It can be found that the Eq. 43 for equilibrium of limit point buckling mode of the FGM arches are the same as those for the homogeneous arches (Lu et al., 2018). Similarly, refer to limit point buckling theory of the homogeneous arches under the locally uniformly distributed radial load, the axial force parameter and the axial force N in the FGM arch during lowest possible limit point buckling are respectively given by
In addition, also refer to the above theory of homogeneous arch, the central displacement corresponding to the lowest buckling load of the FGM arch aswithThe Eq. 46 has a real solution only whenfrom which the modified geometric slenderness of the lowest buckling load of the FGM arch can be obtained once the type of the cross section and the locally uniformly distributed radial load are given.
The FGM arches under the locally uniformly distributed radial load may buckle in a bifurcation mode. The equation of equilibrium of the FGM arch in the primary equilibrium path is obtained from Eq. 29. Similarly, the differential equation of equilibrium of the FGM arch in the case of infinitely close to the bifurcation equilibrium path can be given by
Therefore, the equation for equilibrium of bifurcation mode of the FGM arches can be built by substituting Eq. 29 into Eq. 49 as
It’s can be found that the equation for equilibrium of bifurcation mode of the FGM arches are the same as those for the homogeneous arches (Lu et al., 2018). Similarly, refer to bifurcation buckling theory of the homogeneous arches, the axial force parameter and the axial force N in the arch during bifurcation are respectively given by
In addition, the parameter for bifurcation mode can then be expressed as
By substituting Eqs 51, 52 into Eq. 39, the equation of equilibrium between and P at the upper and lower bifurcation points can be expressed aswherewith
The bifurcation buckling load of a FGM arch under a locally uniformly distributed radial load can be obtained once the localized parameter η is given.
Comparisons of Theoretical Solutions with ANSYS Results
According to Donnell’s shallow shell theory, in the ANSYS models of FGM arch, the geometry of the rectangular section was, the width b = 0.14 m, the total thickness h = 0.07 m, the power-law index p = 1 and the proportional coefficient of homogeneous ceramic layer = 0.2. In addition, the geometric slenderness of the FGM arch S/h = 100. The Young’s modulus of the metal-ceramic material Al-Al2O3 are 70GPa and 380GPa, respectively.
When a multilayer for metal-ceramic material Al-Al2O3 is used to simulate a continuous gradient properties, the accuracy of the FGM arch model increases with an increase of the number of layers of metal-ceramic material Al-Al2O3. As demonstrated in Figure 3, when material distributions of the FGM arch is subjected to Type A (Figure 3A), the cross-section of the FGM arch have 101 layers, the first layer is a homogeneous ceramic layer, and the second to the 101st layers are gradient layers from ceramic to metal. When material distributions of the FGM arch is subjected to Type B (Figure 3B), the cross-section of the FGM arch have 201 layers, the first to the 100st layers are gradient layers from metal to ceramic, the 101st is a homogeneous ceramic layer, and the 102st to the 201st layers are gradient layers from ceramic to metal. When material distributions of the FGM arch is subjected to Type C (Figure 3C), the cross-section of the FGM arch have 202 layers, the first layer is a homogeneous ceramic layer, the second to the 101st layers are gradient layers from ceramic to metal, the 102st to the 201st layers are gradient layers from metal to ceramic, the 101st is a homogeneous ceramic layer, and the last layer is a homogeneous ceramic layer.
FIGURE 3
In addition, the proposed theoretical solution of the buckling equilibrium path, the limit point and bifurcation buckling load for the FGM arches under the locally uniformly distributed radial load are validated by ANSYS finite element results. Convergence researches for using the shell element shell181 of ANSYS show that 80 elements can produce convergent results. Therefore, 80 shell181 elements are chosen to simulate the FGM arches (Figure 3D).
Comparisons of the theoretical solutions of Eqs. 33, 38 for the limit point buckling equilibrium path and the limit point buckling load of the FGM arches with the corresponding ANSYS results are plotted in Figure 4. The variation curves of the dimensionless locally uniformly distributed radial load Q/(2NEc2Θ) with the dimensionless central displacement vc/f are plotted in Figure 4A for material distributions of the FGM arch being subjected to Type A, in Figure 4B for material distributions of the FGM arch being subjected to Type B, and in Figure 4C for material distributions of the FGM arch being subjected to Type C, where the locally uniformly distributed radial load with , the modified geometric slenderness and . It can be seen from Figures 4A–C that the theoretical solutions for the limit point buckling equilibrium path and the limit point buckling load of the FGM arches are highly consistent with the corresponding ANSYS results.
FIGURE 4
Comparisons of the theoretical solutions of Eqs 33, 53 for bifurcation buckling equilibrium path and the bifurcation buckling load of the FGM arches with the corresponding ANSYS results are plotted in Figure 5. The variation curves of the dimensionless locally uniformly distributed radial load Q/(2NEc2Θ) with the dimensionless central displacement vc/f are plotted in Figure 5A for material distributions of the FGM arch being subjected to Type A, in Figure 5B for material distributions of the FGM arch being subjected to Type B, and in Figure 5C for material distributions of the FGM arch being subjected to Type C, where the locally uniformly distributed radial load with , the modified geometric slenderness . It can be seen from Figures 5A–C that the theoretical solutions for bifurcation buckling equilibrium path and the bifurcation buckling load of the FGM arches are highly consistent with the corresponding ANSYS results.
FIGURE 5
Hence, the theoretical solutions of Eqs 33, 38, 53 can be used to predict the buckling equilibrium path, the limit point and bifurcation buckling load for the FGM arches under the locally uniformly distributed radial load.
Limit Point Buckling Analysis
This section investigates the effects of the localized parameter η of locally uniformly distributed radial load, the type and the power-law index p of material distributions of the FGM arch on the limit point buckling equilibrium path and the limit point buckling load for FGM arches under a locally uniformly distributed radial load.
Under the locally uniformly distributed radial load, the limit point buckling equilibrium path obtained from Eqs 33, 38 for FGM arches having different localized parameter η are plotted in Figure 6A with the variation curves of vc/f vs Q/(2NEc2Θ), and in Figure 6B with the variation curves of N/NEc2vs Q/(2NEc2Θ), where the modified geometric slenderness , the power-law index p = 1, the proportional coefficient of homogeneous ceramic layer = 0.2, the geometric slenderness S/h = 100 and the material distributions of FGM arches are subjected to Type A. Also plotted in Figure 6 are the limit points given by Eqs 38, 43. In Figure 6, when the upper limit point is reached, the load system and limit point buckling equilibrium path of FGM arches snaps through to the remote equilibrium branch, and the displacement become very large. Therefore, the upper limit point is the limit point buckling load. It can also be seen from Figure 6 that the localized parameter η has a significant effect on the limit point buckling equilibrium path and the limit point buckling load, the limit point buckling load increases with an increase of localized parameter η, which indicates that FGM arches having a larger length of the locally uniformly distributed radial loading segment are more prone to limit point buckling.
FIGURE 6
In addition, the limit point buckling equilibrium path obtained from Eqs 33, 38 for FGM arches having different type of material distributions are plotted in Figure 7A with the variation curves of vc/f vs Q/(2NEc2Θ), and in Figure 7B with the variation curves of N/NEc2vs Q/(2NEc2Θ), where , η = 0.5, = 0.2, S/h = 100, p = 1. It can be seen from Figure 7 that the type of material distributions has a significant effect on the limit point buckling equilibrium path and the limit point buckling load, Type C has a largest limit point buckling load among Type A, Type B and Type C, which indicates that under the same proportion of metal and ceramic materials, the distribution of ceramics on the upper and lower surfaces of the cross-section can greatly improve the stability of the FGM arch.
FIGURE 7
Similarly, Figure 8A plots the variation curves of vc/f vs Q/(2NEc2Θ) and Figure 8B plots the variation curves of N/NEc2vs Q/(2NEc2Θ) for FGM arches having different power-law index p, where , η = 0.5, = 0.2, S/h = 100 and the material distributions of FGM arches are subjected to Type A. It can also be seen from Figure 8 that the power-law index p has a significant effect on the limit point buckling equilibrium path and the limit point buckling load, the limit point buckling load increases with a decrease of the power-law index p, which indicates that FGM arches having a smaller power-law index p of material distributions are more prone to limit point buckling.
FIGURE 8
Furthermore, the influences of the modified geometric slenderness , the type and the power-law index p of material distributions of the FGM arch, the proportional coefficient of homogeneous ceramic layer and the localized parameter η of the locally uniformly distributed radial load on the limit point buckling load of the FGM arch are studied in detail as well. The variation curves of η vs Q/(2NEc2Θ) are plotted in Figure 9A for a FGM arch having different modified geometric slenderness , in Figure 9B for a FGM arch having different type of material distributions, in Figure 9C for a FGM arch having different power-law index p of material distributions and in Figure 9D for a FGM arch having different proportional coefficient of homogeneous ceramic layer . It can be seen again from Figure 9 that the limit point buckling load increases with an increase of localized parameter η. It can be seen from Figure 9A that the limit point buckling load increases with an increase of modified geometric slenderness . It can be seen again from Figure 9B that Type C has a largest limit point buckling load among Type A, Type B and Type C. In addition, because the composition proportion of ceramics increases with an increase of the proportional coefficient of homogeneous ceramic layer and a decrease of the power-law index p of material distributions of the FGM arch. Therefore, Figures 9C,D show the limit point buckling load increases with an increase of the proportional coefficient of homogeneous ceramic layer and a decrease of the power-law index p of material distributions of the FGM arch, and the power-law index p has a more obvious effect on the limit point buckling load.
FIGURE 9
Bifurcation Buckling Analysis
This section investigates the effects of the localized parameter η of locally uniformly distributed radial load, the type of material distributions of the FGM arch on the bifurcation buckling equilibrium path and the bifurcation buckling load for FGM arches under a locally uniformly distributed radial load.
Under the locally uniformly distributed radial load, the bifurcation buckling equilibrium path obtained from Eqs 33, 53 for FGM arches having different localized parameter η are plotted in Figure 10A with the variation curves of vc/f vs Q/(2NEc2Θ), and in Figure 6B with the variation curves of N/NEc2vs Q/(2NEc2Θ), where the modified geometric slenderness , the power-law index p = 1, the proportional coefficient of homogeneous ceramic layer = 0.2, the geometric slenderness S/h = 100 and the material distributions of FGM arches are subjected to Type A. Also plotted in Figure 10 are the bifurcation buckling points given by Eq. 53. In Figure 10, the upper bifurcation buckling point is on the primary equilibrium path and bifurcation buckling equilibrium path, when this point is reached, the load decreases with an increase of the central displacement, and the dimensionless axial force parameter axial force N/NEc2 still unchanged during this moment. When lower bifurcation buckling point is reached, the equilibrium path of the FGM arch returns from the bifurcation buckling equilibrium path back to the primary equilibrium path.
FIGURE 10
It can also be seen from Figure 10 that the localized parameter η has a significant effect on the bifurcation buckling equilibrium path and the bifurcation buckling load, the bifurcation buckling load increases with an increase of localized parameter η, which indicates that FGM arches having a larger length of the locally uniformly distributed radial loading segment are more prone to bifurcation buckling.
In addition, the bifurcation buckling equilibrium path obtained from Eqs 33, 53 for FGM arches having different type of material distributions are plotted in Figure 11A with the variation curves of vc/f vs Q/(2NEc2Θ), and in Figure 11B with the variation curves of N/NEc2vs Q/(2NEc2Θ), where , η = 0.5, = 0.2, S/h = 100, p = 1. It can be seen from Figure 11 that the type of material distributions has a significant effect on the bifurcation buckling equilibrium path and the limit point buckling load, Type C has a largest bifurcation buckling load among Type A, Type B and Type C, which indicates again that under the same proportion of metal and ceramic materials, the distribution of ceramics on the upper and lower surfaces of the cross-section can greatly improve the stability of the FGM arch.
FIGURE 11
Furthermore, the influences of the localized parameter η of locally uniformly distributed radial load, the type of material distributions and the modified geometric slenderness of the FGM arch on the limit point buckling load and bifurcation buckling load of the FGM arch are studied in detail as well. The variation curves of vs Q/(2NEc2Φ) are plotted in Figure 12A for a Type A FGM arch having different localized parameter η, in Figure 12B for a Type B FGM arch having different localized parameter η and in Figure 12C for a Type C FGM arch having different localized parameter η. It can be seen from Figure 12 that the limit point buckling load and bifurcation buckling load increase with an increase of the modified geometric slenderness and the localized parameter η. It can also be seen from Figure 12 that Type C has a largest limit point buckling load and bifurcation buckling load among Type A, Type B and Type C. In addition, Figure 12 shows that the limit point buckling loads are initially larger than the bifurcation buckling loads, and then smaller than the bifurcation buckling loads while the modified geometric slenderness reaches a certain value. However, when η = 0.25, the limit point buckling loads are larger than the bifurcation buckling loads with any modified geometric slenderness . Therefore, the modified geometric slenderness that is a switches of bifurcation buckling domination and limit point buckling domination may not exist within a certain range of the localized parameter η.
FIGURE 12
Criteria for Switches of Buckling Mode of Functionally Graded Material Arch
It can be seen from the above research that fixed FGM arches under a locally uniformly distributed radial load may buckle in a limit point buckling mode or in a bifurcation buckling mode. The buckling equilibrium path for FGM arches having different modified geometric slenderness are plotted in Figure 13 with the variation curves of vc/f vs Q/(2NEc2Θ), where η = 0.75, = 0.2, S/h = 100, p = 1 and the material distributions of FGM arches are subjected to Type A. Figure 13A shows that when , the buckling equilibrium path is in the lowest limit point buckling mode, which only has a inflection point, this mode is a switches of limit point buckling and no limit point buckling. Figure 13B shows when , the buckling equilibrium path is in the limit point buckling mode, which has two limit point. Figure 13C shows when , the buckling equilibrium path is in the lowest bifurcation buckling mode, which only has a bifurcation buckling point, this mode is a switches of bifurcation buckling and limit point buckling. Figure 13D shows when , the buckling equilibrium path is in the bifurcation buckling mode, which has two bifurcation buckling point. Figure 13E shows when , the buckling equilibrium path is in the lowest bifurcation buckling dominant mode, which has the same upper limit point and upper bifurcation buckling point, this mode is a switches of bifurcation buckling domination and limit point buckling domination. Figure 13F shows when , the buckling equilibrium path is in the bifurcation buckling dominant mode. Therefore, the modified geometric slenderness of a fixed FGM arch play an important role in determining its buckling mode.
FIGURE 13
The modified geometric slenderness λs obtained from Eq. 48 can make fixed FGM arches lose its stability in the lowest limit buckling mode. This modified geometric slenderness λsl can be expressed as
The modified geometric slenderness λs obtained from Eq. 53 can make fixed FGM arches lose its stability in the lowest bifurcation buckling mode. This modified geometric slenderness λsb1 can be expressed as
When the fixed FGM arches lose its stability in the lowest bifurcation buckling dominant mode, the modified geometric slenderness λsb2 corresponding to this mode can be obtained by equaling the upper bifurcation instability load to the upper limit point buckling load.
Typical variation curves of the localized parameter η vs the modified geometric slenderness λs are shown in Figure 14A for FGM arches having different type of material distributions, in Figure 14B for FGM arches having different power-law index p of material distributions and in Figure 14C for FGM arches having different proportional coefficient of homogeneous ceramic layer , where S/h = 100. It can be seen from Figure 14 that the localized parameter η has a significant effect on the modified geometric slenderness λsb1, λsb2, and a slightly effect on the modified geometric slenderness λsl, the modified geometric slenderness λsl, λsb1, and λsb2 increase with a decrease of the localized parameter η of the locally uniformly distributed radial load. Figure 14 defines four regions, region I corresponding to FGM arches with λs < λsl is subjected to the no buckling mode, region II corresponding to FGM arches with λsl < λs < λsb1 is subjected to the limit point buckling mode, region III corresponding to FGM arches with λsb1 < λs < λsb2 is subjected to the bifurcation buckling mode, region IV corresponding to FGM arches with λs > λsb2 is subjected to the bifurcation buckling dominant mode. In addition, it can be seen from Figure 14A that the type of material distributions has a significant effect on the modified geometric slenderness λsl, λsb1 and λsb2, Type C has largest modified geometric slenderness λsl, λsb1 and λsb2 among Type A, Type B and Type C, Figures 14B,C shows the modified geometric slenderness λsl, λsb1 and λsb2 increase with an decrease of the power-law index p of material distributions and the proportional coefficient of homogeneous ceramic layer .
FIGURE 14
Conclusion
Theoretical studies of the nonlinear in-plane buckling of a fixed circular functionally graded material (FGM) arches under a locally uniformly distributed radial load have been explored in this paper. It was found that fixed FGM arches under a locally uniformly distributed radial load may buckle in a limit point buckling mode or in a bifurcation buckling mode. Theoretical solutions for the limit point buckling load and bifurcation buckling loads and nonlinear equilibrium path of the fixed FGM arches that are subjected to three different types of material distributions were obtained. It was found that the type of material distributions of the fixed FGM arches slight the limit point buckling loads and bifurcation buckling loads as well as the nonlinear equilibrium path significantly. Type C had a largest limit point buckling load and a largest bifurcation buckling load among Type A, Type B and Type C, it indicated that under the same proportion of metal and ceramic materials, the distribution of ceramics on the upper and lower surfaces of the cross-section can greatly improve the stability of the FGM arch. In addition, the limit point buckling load and bifurcation buckling load also increase with an increase of the modified geometric slenderness and the localized parameter η and the proportional coefficient of homogeneous ceramic layer as well as a decrease of the power-law index p of material distributions of the FGM arch.
The critical modified geometric slendernesses λsl, λsb1, and λsb2 that are switches of buckling modes were also derived. The localized parameter η has a significant effect on the modified geometric slenderness λsb1, λsb2, and limited effect on the modified geometric slenderness λsl, the modified geometric slenderness λsl, λsb1, and λsb2 increase with a decrease of the localized parameter η of the locally uniformly distributed radial load. It was found that Region I corresponding to FGM arches with λs < λsl is subjected to the no buckling mode, region II corresponding to FGM arches with λsl < λs < λsb1 is subjected to the limit point buckling mode, region III corresponding to FGM arches with λsb1 < λs < λsb2 is subjected to the bifurcation buckling mode, region IV corresponding to FGM arches with λs > λsb2 is subjected to the bifurcation buckling dominant mode. In addition, it was also found that the type of material distributions has a significant effect on the modified geometric slenderness λsl, λsb1 and λsb2, Type C has the largest modified geometric slenderness λsl, λsb1 and λsb2 among Type A, Type B and Type C, and the modified geometric slenderness λsl, λsb1 and λsb2 increase with an decrease of the power-law index p of material distributions and the proportional coefficient of homogeneous ceramic layer . Theoretical solutions for FGM arches under a locally uniformly distributed radial load derived in this paper would be useful as a reference for the design of buckling of the FGM arches in the future.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
HL: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing-review and editing. JZ: Formal analysis, Validation, Writing-review and editing. ZY: Formal analysis, Investigation, Methodology, Writing-original draft. AL: Data curation, funding acquisition. JZ: Project administration.
Funding
This work was funded by (Characteristic innovation projects of colleges and universities in Guangdong Province) with grant number (2019KTSCX190) and (National Natural Science Foundation of China) with grant number (51908146).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
Al-shujairiM.MollamahmutoğluC. (2018). Buckling and Free Vibration Analysis of Functionally Graded sandwich Microbeams Resting on Elastic Foundation by Using Non Local Strain Gradient Theory in Conjunction with Higher Order Shear Theories under thermal Effect. Compos. Part. B154, 292–312. 10.1016/j.compositesb.2018.08.103
2
AsgariH.BateniM.KianiY.EslamiM. R. (2014). Non-linear Thermo-Elastic and Buckling Analysis of FGM Shallow Arches. Compos. Struct.109, 75–85. 10.1016/j.compstruct.2013.10.045
3
BateniM.EslamiM. R. (2015). Non-linear In-Plane Stability Analysis of FG Circular Shallow Arches under Uniform Radial Pressure. Thin-Walled Structures94, 302–313. 10.1016/j.tws.2015.04.019
4
BateniM.EslamiM. R. (2014). Non-linear In-Plane Stability Analysis of FGM Circular Shallow Arches under central Concentrated Force. Int. J. Non-Linear Mech.60, 58–69. 10.1016/j.ijnonlinmec.2014.01.001
5
BradfordM. A.UyB.PiY.-L. (2002). In-plane Elastic Stability of Arches under a central Concentrated Load. J. Eng. Mech.128, 710–719. 10.1061/(asce)0733-9399(2002)128:7(710)
6
CaiJ.FengJ. (2010). Buckling of Parabolic Shallow Arches when Support Stiffens under Compression. Mech. Res. Commun.37, 467–471. 10.1016/j.mechrescom.2010.05.004
7
CaiJ.XuY.FengJ.ZhangJ. (2012). In-plane Elastic Buckling of Shallow Parabolic Arches under an External Load and Temperature Changes. J. Struct. Eng.138, 1300–1309. 10.1061/(asce)st.1943-541x.0000570
8
CaiJ.ZhouY.FengJ. (2013). Post-buckling Behavior of a Fixed Arch for Variable Geometry Structures. Mech. Res. Commun.52, 74–80. 10.1016/j.mechrescom.2013.07.002
9
ChenM.-T.YoungB. (2021). Beam-column Design of Cold-Formed Steel Semi-oval Hollow Non-slender Sections. Thin-walled Structures162, 107376. 10.1016/j.tws.2020.107376
10
ChenM.-T.YoungB.MartinsA. D.CamotimD.DinisP. B. (2021). Experimental Investigation on Cold-Formed Steel Lipped Channel Beams Affected by Local-Distortional Interaction under Non-uniform Bending. Thin-walled Structures161, 107494. 10.1016/j.tws.2021.107494
11
HanQ.ChengY.LuY.LiT.LuP. (2016). Nonlinear Buckling Analysis of Shallow Arches with Elastic Horizontal Supports. Thin-Walled Structures109, 88–102. 10.1016/j.tws.2016.09.016
12
HuC.-F.LiZ.HuQ.-S. (2021). On Non-linear Behavior and Buckling of Arch-Beam Structures. Eng. Structures239 (2), 112214. 10.1016/j.engstruct.2021.112214
13
LuH.LiuA.PiY.-L.BradfordM. A.FuJ.HuangY. (2018). Localized Loading and Nonlinear Instability and post-instability of Fixed Arches. Thin-Walled Structures131, 165–178. 10.1016/j.tws.2018.06.019
14
LuH. W.LiuL. L.LiuA. R.PiY.-L.BradfordM. A.HuangY. H. (2020). Effects of Movement and Rotation of Supports on Nonlinear Instability of Fixed Shallow Arches. Thin-Walled Structures122, 2020. 10.1016/j.tws.2020.106909
15
MoitaJ. S.AraújoA. L.CorreiaV. F.Mota SoaresC. M.HerskovitsJ. (2018). Material Distribution and Sizing Optimization of Functionally Graded Plate-Shell Structures. Composites B: Eng.142, 263–272. 10.1016/j.compositesb.2018.01.023
16
PiY.-L.BradfordM. A.LiuA. (2017). Nonlinear Equilibrium and Buckling of Fixed Shallow Arches Subjected to an Arbitrary Radial Concentrated Load. Int. J. Str. Stab. Dyn.17 (8), 1750082. 10.1142/s0219455417500821
17
PiY.-L.BradfordM. A. (2009). Non-linear In-Plane Postbuckling of Arches with Rotational End Restraints under Uniform Radial Loading. Int. J. Non-Linear Mech.44 (9), 975–989. 10.1016/j.ijnonlinmec.2009.07.003
18
PiY.-L.BradfordM. A.Tin-LoiF. (2008). Non-linear In-Plane Buckling of Rotationally Restrained Shallow Arches under a central Concentrated Load. Int. J. Non-Linear Mech.43 (1), 1–17. 10.1016/j.ijnonlinmec.2007.03.013
19
PiY.-L.BradfordM. A.UyB. (2002). In-plane Stability of Arches. Int. J. Sol. Structures39 (1), 105–125. 10.1016/s0020-7683(01)00209-8
20
PlautR. H.RaymondH. (1990). Buckling of Shallow Arches with Supports that Stiffen when Compressed. J. Eng. Mech.116 (4), 973–976. 10.1061/(asce)0733-9399(1990)116:4(973)
21
RastgoA.ShafieH.AllahverdizadehA. (2005). Instability of Curved Beams Made of Functionally Graded Material under Thermal Loading. Int. J. Mech. Mater. Des.2, 117–128. 10.1007/s10999-005-4446-3
22
SchreyerH. L. (1972). The Effect of Initial Imperfections on the Buckling Load of Shallow Circular Arches. J. Appl. Mech.39 (2), 445–450. 10.1115/1.3422698
23
SimsekM. (2016). Buckling of Timoshenko Beams Composed of Two-Dimensional Functionally Graded Material (2D-FGM) Having Different Boundary Conditions. Compos. Struct.49, 304–314. 10.1016/j.compstruct.2016.04.034
24
SongX.LiS. R. (2008). Nonlinear Stability of Fixed-Fixed Fgm Arches Subjected to Mechanical and thermal Loads. Amr33-37, 699–706. 10.4028/www.scientific.net/amr.33-37.699
25
TimoshenkoS. P.GereJ. M.PragerW.GereJ. M.PragerW. (1962). Theory of Elastic Stability, Second Edition. J. Appl. Mech.29 (1), 220–221. 10.1115/1.3636481
26
YanS.-t.ShenX.ChenZ.JinZ. (2017). On Buckling of Non-uniform Shallow Arch under a central Concentrated Load. Int. J. Mech. Sci.133, 330–343. 10.1016/j.ijmecsci.2017.08.046
27
YangZ.FengC.YangJ.WangY.LvJ.LiuA.et al (2020a). Geometrically Nonlinear Buckling of Graphene Platelets Reinforced Dielectric Composite (GPLRDC) Arches with Rotational End Restraints. Aerospace Sci. Tech.107, 106326. 10.1016/j.ast.2020.106326
28
YangZ.ZhaoS.YangJ.LvJ.LiuA.FuJ. (2020b). In-plane and Out-Of-Plane Free Vibrations of Functionally Graded Composite Arches with Graphene Reinforcements. Mech. Adv. Mater. Structures, 1–11. 10.1080/15376494.2020.1716420
29
YangZ.TamM.ZhangY.KitipornchaiS.LvJ.YangJ. (2020c). Nonlinear Dynamic Response of FG Graphene Platelets Reinforced Composite Beam with Edge Cracks in thermal Environment. Int. J. Str. Stab. Dyn.20, 2043005. 10.1142/s0219455420430051
30
YangZ.LiuA.YangJ.LaiS.-K.LvJ.FuJ. (2021a). Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a central point Load. Materials14 (8), 2026. 10.3390/ma14082026
31
YangZ.WuD.YangJ.LaiS.-K.LvJ.LiuA.et al (2021b). Dynamic Buckling of Rotationally Restrained FG Porous Arches Reinforced with Graphene Nanoplatelets under a Uniform Step Load. Thin-Walled Structures166, 108103. 10.1016/j.tws.2021.108103
Summary
Keywords
functionally graded material, fixed arch, limit point buckling, bifurcation buckling, critical modified geometric slenderness
Citation
Lu H, Zhou J, Yang Z, Liu A and Zhu J (2021) Nonlinear Buckling of Fixed Functionally Graded Material Arches Under a Locally Uniformly Distributed Radial Load. Front. Mater. 8:731627. doi: 10.3389/fmats.2021.731627
Received
28 June 2021
Accepted
27 July 2021
Published
09 August 2021
Volume
8 - 2021
Edited by
Lin Xu, Jilin University, China
Reviewed by
Man-Tai Chen, Shanghai Jiao Tong University, China
Qiang Feng, Shandong University of Science and Technology, China
Updates
Copyright
© 2021 Lu, Zhou, Yang, Liu and Zhu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhicheng Yang, zhicheng.yang@zhku.edu.cn; Airong Liu, liuar@gzhu.edu.cn
This article was submitted to Mechanics of Materials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.