Engineering and Computer Science

# Free Vibration and Buckling Analysis of Parabolic Frame Structures

Oguzhan Das, Hasan Öztürk, Can Gönenli,
Article Date Published : 2 December 2019 | Page No.: EC-2019-295-306 | Google Scholar

## Abstract

Numerous engineering structures take place in many fields such as construction, automotive, naval industry, aerospace, etc. Parabolic frame structures have a significant role in those industries. There are many gaps in terms of the literature studies about these structures. In this paper, free vibration and buckling behavior of parabolic frame structures by using the finite element method are investigated. The structure is modeled by considering a frame element that has three degrees of freedom and neglects the displacement in y-axis. In addition, two-bay parabolic frame structure with different radius of curvature is investigated. The numerical results are compared with the CAD model of the structure by using SolidWorks for various cases. It is concluded that the results are in very good agreement with those results that are obtained from SolidWorks. It is also understood that for the different radius of curvatures the approach that is used for finite element buckling analysis in perpendicular distributed loading case does not change.

Keywords: Finite Element Analysis, vibration, buckling, frame structure, parabolic frames

## Introduction

In this paper, the free vibration and buckling analysis of the single-bay and the two-bay isotropic parabolic frames are investigated. It is seen that there is not any study related to this topic in literature. The parabolic frame structures are modelled by considering Euler-Bernoulli beam theory. All equations are formulated and solved in MATLAB. In order to compare numerical results, the same structures are modelled and analyzed in SolidWorks.

## Theoretical Analysis

Finite Element Model

The finite element method is employed in order to model the parabolic frame structure as seen in Fig.1. To satisfy all degrees of freedom for such structure, a two-nodded frame element with three degrees of freedom (DOF), which is shown in Fig.2, is taken into account.

Figure 1: Parabolic frame structure

Figure 2: The frame element

The longitudinal displacement u is represented with a first-order polynomial equation while the bending displacement v, is defined as a cubic function in order to satisfy the rotation effect.

(1) u=a_1+a_2 x (2) v=a_3+a_4 x+a_5 x^2+a_6 x^3

The slope v', is equal to the first derivative of the bending displacement with respect to x.

v^'=dv/dx

The generalized displacement vector for the frame element can be expressed as

q=[u_1 v_1 〖v^'〗_1 u_2 v_2 〖v^'〗_2]

The strain energy for frame element is given as

where E is the modulus of elasticity of the material, A is the cross-sectional area of the element and I is the second moment of inertia of the element. Eq.5 can be written in matrix form as

where ke is the element geometry matrix.

The Kinetic Energy for frame element is

T_e=1/2 ∫_0^L▒〖ρA(u ̇^2+v ̇^2)〗 dx

where ρ is the density of the material. Writing Eq.7 in matrix form as

where me is the element mass matrix. The work done by the applied load is expressed as

where

where P represents the point load. Hence,

Writing Eq.11 in matrix form,

where kge is the element geometry matrix.

Figure 3: Transformation of local coordinates

k ̅_e=T^T k_e T, m ̅_e=T^T m_e T ,and k ̅_ge=T^T k_ge T

where T is the transformation matrix which is,

T=[■(cosθ&sinθ&0@-sinθ&cosθ&0@0&0&1)]

Free Vibration Analysis

The dynamic response of the parabolic frame structure can be obtained by using Lagrange’s equation of motion as

Mq ̈+Kq=0

and the eigenvalue problem is

where K is the global stiffness matrix, M is the global mass matrix and λ is the eigenvalue or, in other words, the natural frequency values of the system.

Buckling Analysis

The critical buckling load for a system is found by using the eigenvalue problem given in Eq. 17.

where Kg is the global geometry matrix and Pcr is the critical buckling load.

## Numerical Analysis

In this study, free vibration and buckling analysis of parabolic frames are investigated. Four different radius of curvature values are selected in order to find out the effect of the curvature on the dynamic behavior and critical buckling load of the structure. To also investigate the thickness effect, two different thickness parameter is taken into account. Additionally, the two-bay structure is also considered for free vibration and buckling analysis. Material properties and geometric properties of the structure are given in Table 1. A MATLAB code is developed in order to solve the free vibration and buckling problem. The distributed loading case given in Fig.4 is taken into account to solve the buckling problem. In order to simplify the problem, the distributed load is transformed into a two-point load that acts to the columns of the structure as seen in Fig.5.

 Elastic Modulus E 69 GPa Density ? 2700 kg/m 3 Poisson Ratio ? 0.33

Table 1: Material properties and geometry of the structure

(3) (K-P_cr K_g)q=0

Figure 5: Equivalent loading case for parabolic frame structure

Effect of the radius of curvature

For a single bay structure, the effect of the radius of curvature (Rxx) on dynamic behavior and critical buckling load are investigated. Different radius of curvatures is taken into account as is seen in Fig.6.

Figure 6: The parabolic frame structure with different curvatures

Comparative results of free vibration analysis and critical buckling load of the parabolic frame structure for various curvature values are given in Tables 2 and 3.

 N.F. (Hz) R xx =1.5L b R xx =2L b SW PS Err (%) SW PS Err (%) f 1 8.08 8.28 2.49% 8.13 8.29 1.86% f 2 30.49 31.87 4.49% 31.09 32.05 3.08% f 3 50.98 52.82 3.61% 51.56 53.02 2.82% f 4 61.34 60.22 1.84% 60.54 59.61 1.53% f 5 112.14 115.46 2.96% 113.54 115.66 1.87% P cr (N) 945.49 948.50 0.32% 951.26 948.90 0.25%

The first five natural frequency values and the critical buckling load for four different curvatures are given in Tables 2 and 3. It is concluded that the present study is in very good agreement with SolidWorks results. As the curvature increases, the first three and the fifth natural frequency increase whereas the fourth natural frequency decreases. This is because all components of the structure have the same displacement characteristics in the fourth mode as seen in Fig. 7d. All mode shapes are the same for all curvatures and given in Fig.7. The critical buckling load increases as the radius of curvature increases.

Figure 7: The first five modes of the parabolic frame structure

Comparative results of free vibration analysis and critical buckling load of the straight frame structure (Rxx=∞) for various curvature values are given in Table 4. It is seen that the first three and fifth natural frequencies are at the highest values when compared with the results given in Tables 2 and 3. In contrast, the fourth natural frequency is at the lowest value as is expected. The critical buckling load is also at its highest value when compared with Tables 2 and 3.

 N.F. (Hz) R xx =2.5 L b R xx =3 L b SW PS Err(%) SW PS Err (%) f 1 8.16 8.29 1.61% 8.17 8.30 1.49% f 2 31.43 32.16 2.30% 31.64 32.23 1.87% f 3 51.88 53.14 2.41% 52.07 53.22 2.19% f 4 60.05 59.26 1.33% 59.68 59.03 1.09% f 5 114.23 115.79 1.37% 114.62 115.87 1.09% P cr (N) 954.54 949.08 0.57% 956.4 949.17 0.76%
 N.F. (Hz) Straight Frame Structure ( R xx =?) SW PS Err (%) f 1 8.24 8.13 1.21% f 2 32.59 32.13 1.03% f 3 52.98 52.25 1.37% f 4 57.76 56.88 0.81% f 5 116.40 114.80 1.04% P cr (N) 964.05 962.19 1.32%

For a single bay structure, the effect of the thickness (h) on dynamic behavior and critical buckling load is investigated for fixed boundary condition. Two different thickness parameters, h1=a/100 and h2=a/50, are selected, and the curvature is set as Rxx=2Lb.

Table 5: Free vibration and buckling analysis of the parabolic frame structure for h1=Lb/100 h2=Lb/50 (N.F.: Natural Frequency Values, Pcr: Critical Buckling Load, SW: SolidWorks Results, PS: Present Study Results, Err: Error)

 N.F. (Hz) h 1 = L b /100 h 2 = L b /50 SW PS Err (%) SW PS Err (%) f 1 8.13 8.29 1.86% 16.21 16.57 2.21% f 2 31.09 32.05 3.08% 62.15 64.06 3.07% f 3 51.56 53.02 2.82% 102.59 106.03 3.35% f 4 60.54 59.61 1.53% 120.73 119.07 1.37% f 5 113.54 115.66 1.87% 226.16 230.81 2.06% P cr (N) 950.13 948.90 0.25% 15337 15180 1.02%

As seen in Table 5, two different thickness (h) parameter is compared each other in terms of natural frequency and critical buckling load. The natural frequency values increases or decreases as the same as the thickness ratio as expected. The critical buckling load increases as the structure thickness increases. However, the increment is not the same as that of the natural frequency. This is because the ratio is not the same since the relation between the stiffness matrix and mass matrix is not the same as the stiffness matrix and geometric matrix.

1. The Two-Bay structure

In addition to the single-bay structure, the two-bay parabolic frame structure is also investigated. Free vibration analysis and buckling analysis are performed under all columns fixed boundary conditions. The curvature is taken as Rxx=2Lb and the thickness is set as h=Lb/100.

The loading condition is given as is in Fig.8. The equivalent loading case for finite element analysis is given in Fig.9.

Figure 8: Distributed loading case for the two-bay parabolic frame structure

Figure9: Equivalent loading case for the two-bay parabolic frame structure

Comparative results of the free vibration analysis and buckling analysis for the two-bay parabolic frame structure is given in Tables 6 and 7.

 N.F. (Hz) R xx =1.5L b R xx =2L b SW PS Err (%) SW PS Err (%) f 1 7.42 7.66 3.29% 7.47 7.67 2.71% f 2 29.48 31.01 5.19% 29.97 31.12 3.86% f 3 34.45 37.73 9.52% 35.74 38.26 7.04% f 4 51.13 53.29 4.22% 51.67 53.46 3.48% f 5 59.98 59.77 0.34% 59.37 59.30 0.12% P cr (N) 1429.1 1453.7 1.72% 1439.8 1455.3 1.08%
 N.F. (Hz) R xx =2.5L b R xx =3L b SW PS Err (%) SW PS Err (%) f 1 7.50 7.68 2.39% 7.51 7.68 2.29% f 2 30.23 31.19 3.20% 30.38 31.24 2.83% f 3 36.50 38.57 5.70% 36.97 38.78 4.90% f 4 51.93 53.57 3.15% 52.10 53.63 2.94% f 5 58.94 59.03 0.14% 58.65 58.85 0.34% P cr (N) 1437.2 1456.1 1.32% 1448.3 1456.7 0.58%

Table 6: Free vibration and buckling analysis of the two-bay parabolic frame structure for Rxx=1.5Lb and Rxx=2Lb (N.F.: Natural Frequency Values, Pcr: Critical Buckling Load, SW: SolidWorks Results, PS: Present Study Results, Err: Error)

Table 7: Free vibration and buckling analysis of the two-bay parabolic frame structure for Rxx=2.5Lb and Rxx=3Lb (N.F.: The Natural Frequency Values, Pcr: The Critical Buckling Load, SW: SolidWorks Results, PS: Present Study Results, Err: Error)

As seen in Tables 6 and 7, agreement between the present finite element models and SolidWorks solution results is good. Differently from the single-bay results, the fifth natural frequency value decreases as the curvature increases. This is because the fifth mode is similar to the fourth mode of the single-bay structure which all components of the structure behave have the same displacement characteristics.

The critical buckling load changes in the same way as the curvature.

Figure 10: The first five modes of the two-bay parabolic frame structure

## Conclusions and Discussions

In this paper, free vibration and buckling analysis of the parabolic frame structure are investigated by using the finite element method. According to the results given in section 3, it is concluded that

• The natural frequency values and critical buckling loads that are obtained via the finite element method by using a frame element are in very good agreement with SolidWorks results for both parabolic and straight frames.

• For the single-bay structure, all natural frequency values, except the fourth natural frequency value, increase as the radius of curvature increases. On the other hand, the fourth frequency value decreases. This is because the structural components behave in the same displacement characteristic in this mode.

• For the two-bay structure, all natural frequency values, except the fifth natural frequency value, increase as the radius of curvature increases. Similar to the single-bay structure, this is because the structural components behave in the same displacement characteristic in this mode.

• The critical buckling load increases when the radius of curvature increases for both single-bay and multi-bay parabolic frame structure.

• The equivalent loading approach that is used for straight frame structures in order to perform finite element analysis is also valid for parabolic frame structures under perpendicular distributed loading condition.

• As is expected, the natural frequency and the critical buckling values increase as the thickness of the parabolic frame structure increase.

## References

1. “Frames and Arches Condensed Solutions for Structural Analysis V. Leontovich. In: ed. .
2. Vibrations of Planar Curved Beams, Rings, and Arches Chidamparam P., Leissa A. W.. Applied Mechanics Reviews.1993;46(9). CrossRef Google Scholar
3. Free vibration analysis of planar curved beams by wave propagation Kang B., Riedel C.H., Tan C.A.. Journal of Sound and Vibration.2003;260(1). CrossRef Google Scholar
4. ANALYSIS OF THREE-DIMENSIONAL LOCKING-FREE CURVED BEAM ELEMENT ZHU Z. H., MEGUID S. A.. International Journal of Computational Engineering Science.2004;05(03). CrossRef Google Scholar
5. Free vibration analysis of axially functionally graded tapered Bernoulli–Euler microbeams based on the modified couple stress theory Akgöz Bekir, Civalek Ömer. Composite Structures.2013;98. CrossRef Google Scholar
6. Free vibration of Euler and Timoshenko functionally graded beams by Rayleigh–Ritz method Pradhan K.K., Chakraverty S.. Composites Part B: Engineering.2013;51. CrossRef Google Scholar
7. Free vibration analysis of functionally graded size-dependent nanobeams Eltaher M.A., Emam Samir A., Mahmoud F.F.. Applied Mathematics and Computation.2012;218(14). CrossRef Google Scholar
8. A finite beam element for vibration analysis of rotating tapered timoshenko beams Bazoune A., Khulief Y.A.. Journal of Sound and Vibration.1992;156(1). CrossRef Google Scholar
9. Vibration Of Rings With Unsymmetrical Cross-Sections: A Finite Element Approach Cleghorn W.L., Tabarrok B., Lee T.W.. Journal of Sound and Vibration.1993;168(1). CrossRef Google Scholar
10. Dynamic analysis of very flexible beams Fotouhi R.. Journal of Sound and Vibration.2007;305(3). CrossRef Google Scholar
11. In-plane free vibration of a pre-stressed curved beam obtained from a large deflected cantilever beam Ozturk Hasan. Finite Elements in Analysis and Design.2011;47(3). CrossRef Google Scholar
12. Free Vibration Analysis of Delaminated Composite Beams Kıral Binnur Gören. Science and Engineering of Composite Materials.2009;16(3). CrossRef Google Scholar
13. Improved method of free vibration analysis of frame structures Basci M.I., Toridis T.G., Khozeimeh K., Fettahlioglu O.A.. Computers & Structures.1979;10(1-2). CrossRef Google Scholar
14. Buckling Analysis of Structures Composed of Tapered Members Chan Siu Lai. Journal of Structural Engineering.1990;116(7). CrossRef Google Scholar
15. Vibration Analysis Of Beam Subjected To Moving Loads Using Finite Element Method Mehmood Arshad. IOSR Journal of Engineering.2014;4(5). CrossRef Google Scholar

• ## Can Gönenli

International Journal of Scientific Research and Management, 2019.

### Article Details

Issue: Vol 7 No 12 (2019)
Page No.: EC-2019-295-306
Section: Engineering and Computer Science
DOI:

### How to Cite

Das, O., Öztürk, H., & Gönenli, C. (2019). Free Vibration and Buckling Analysis of Parabolic Frame Structures. International Journal of Scientific Research and Management, 7(12), EC-2019. https://doi.org/10.18535/ijsrm/v7i12.ec02