International Journal Modern Trends in Engineering and Research www.ijmter.com e-issn :2349-9745, Date: 28-30 April, 20 STATIC BENDING ANALYSIS OF AN ISOTROPIC CIRCULAR PLATE USING FINITE ELEMENT METHOD Pratik S. Gujar 1 and Sandeep S. Pendhari 2 1 Lecturer, Department Structural Engineering, V.J.T.I., Mumbai (M.S.)gujar.pratik@gmail.com 2 Associate Pr., Department Structural Engineering, V.J.T.I., Mumbai (M.S.)sspendhari@vjti.org.in Abstract- An attempt is made in this study to perform static bending analysis an isotropic circular plate using four noded shell 63 element which is based on Kirchhf s Plate Theory. Modelling and analysis an isotropic circular plate is done by FEA stware ANSYS APDL. In this work effect various types meshes on deflections and bending stresses in circular plate is investigated and finally convergence FEA results analytical solution is studied. The results are obtained for aspect ratio 10 and 20. Simply supported and clamped boundary conditions uniformly distributed load are considered for present study. Keywords-Isotropic circular plate; Finite element method; 4 noded shell 63 element; Multi zone ;. I. INTRODUCTION Finite element analysis is a numerical technique which is widely used for solving many engineering problems having complex geometries, boundary conditions and loads. From last few decades, the fast improvement in computer hardware technology has boosted this method. Now there are many popular finite element analysis stware packages available. ANSYS is one them. Finite element analysis stware ANSYS is capable to analyze a wide range different problems. In the recent few years, there is an enormous amount research work is going on the structural analysis including the flexural or buckling thin or thick beams and plates. Still this area is open for research both at fundamental and applied level. Vanam et al. [2] presented static analysis an isotropic rectangular plate using four noded quadrilateral element. The aim study was static bending analysis an isotropic rectangular plate various boundary conditions and various types load applications. Numerical analysis has been carried out by developing programming in mathematical stware MATLAB and analyzed by using finite element analysis stware ANSYS.Kavade et al. [3] studied circular plate a hole for various types meshes. The effect an initial stretching a rectangular plate a cylindrical hole on the stress and displacement distributions around the hole was studied. Boedo analyzed [4] pre-tensioned clamped circular plate axisymmetric load for two different elements shell 63 and shell 93 in ANSYS. Results obtained are compared analytical results.vrabie and Baetu [5] analyzed an isotropic reinforced concrete plate using shell 63 and shell 181 element. Also effort was made to analyze same problem in SAP2000. Relative error in displacement by both the element and stwares is compared to a 3-D analysis. The present work aims to study effect various types meshes on the flexural deformation an isotropic circular plate. The results obtained from FEA are compared analytical results and convergence criteria is studied. @IJMTER-2014, All rights Reserved
II. FINITE ELEMENT ANALYSIS Static bending analysis simply supported and clamped circular plate subjected to uniformly distributed load is studied. Here effect various types meshes on flexural properties plate is studied along convergence FEA solution to analytical solution based on Kirchhf s Plate Theory. Following are the values which are assumed. Diameter plate = 200 mm Aspect ratio = 10 and 20 Poisson s ratio steel (µ) = 0.3 Young s modulus steel (E) = 200 GPa Distributed load = 275 kpa Three types meshes are considered in the present study. mesh triangular element, quadrilateral element and. Figure 1, 2 and 3shows various types meshes considered for present work. Number divisions are considered along the circumference as, 40, 80,120,0 and 200. To simulate simply supported boundary conditions the translational degrees freedoms are constrained and for fixed boundary conditions, all six degrees freedom are constrained. Figure 2: triangular element Figure 2: quadrilateralel ement @IJMTER-2014, All rights Reserved 493
Figure 3: III. RESULTS AND DISCUSSION An isotropic circular plate is analyzed in FEA stware ANSYS for various types meshes. Results are obtained using 4 noded shell 63 element and compared analytical results. Finally convergence maximum displacement and maximum bending stress analytical solution is studied. Table 1. S. S. plate UDL (Aspect ratio = 20) 0.0949 0.0930 0.0874 34.0 33.8 33.2 40 0.0955 0.0953 0.0950 34.0 34.1 34.0 80 0.0956 0.0954 0.0955 34.0 34.0 34.0 0.0956 34.03 120 0.0956 0.0955 0.0955 34.0 34.0 34.0 0 0.0956 0.0956 0.0955 34.0 34.0 34.0 200 0.0956 0.0956 0.0956 34.0 34.0 34.0 @IJMTER-2014, All rights Reserved 494
Table 2. Clamped plate UDL (aspect ratio = 20) 0.0234 0.0242 0.0230 20.0 20.1 20.1 40 0.0234 0.0236 0.0235 19.2 19.2 19.1 80 0.0235 0.0235 0.0235 18.9 18.8 18.8 0.02346 20.6 120 0.0235 0.0234 0.0234 18.7 18.7 18.7 0 0.0235 0.0235 0.0234 18.7 18.7 18.7 200 0.0235 0.0234 0.0235 18.7 18.6 18.7 Table 3. S. S. plate UDL (Aspect ratio = 10) 0.0119 0.01 0.0109 8.49 8.44 8.31 40 0.0119 0.0119 0.0119 8.51 8.53 8.50 80 0.0120 0.0119 0.0119 8.51 8.50 8.50 0.01196 8.507 120 0.0120 0.0119 0.0119 8.51 8.51 8.50 0 0.0120 0.0119 0.0119 8.51 8.50 8.50 200 0.0120 0.0119 0.0119 8.51 8.50 8.50 @IJMTER-2014, All rights Reserved 495
Table 4. Clamped plate UDL (aspect ratio = 10) 0.00292 0.00303 0.00288 5. 5.17 5.19 40 0.00293 0.00295 0.00294 4.95 4.95 4.94 80 0.00293 0.00293 0.00294 4.87 4.86 4.86 0.00293 5.156 120 0.00293 0.00293 0.00293 4.84 4.83 4.84 0 0.00293 0.00293 0.00293 4.83 4.82 4.83 200 0.00293 0.00293 0.00293 4.82 4.82 4.82 In case simply supported circular plate UDL for aspect ratio 20, maximum displacement converges to analytical solution is at 80 no. div. for multi zone while for free quadrilateral element and free triangular element, it converges at 0 and 200 no. div. respectively. Similarly max. stress converges to analytical solution at no. div. for multi zone while for free quadrilateral element and free triangular element, it converges at 40 no. div. respectively (Table 1).On the other hand, clamped plate UDL for aspect ratio 20, maximum displacement converges to analytical solution at no. div. for multi zone, while for both free quadrilateral element and free triangular element, it converges at 120 no. div. respectively. Similarly max. stress converges to analytical solution at no. div. for all type (Table 2). Similar behavior is observed for simply supported (Table 3)and clamped boundary condition(table 4) for aspect ratio 10. IV. CONCLUDING REMARK Static bending analysis an isotropic circular plate is carried out using shell 63 element. Effect various types meshes on deflections and bending stresses is investigated. shows good agreement results the analytical solution hence it is preferred over other types. On the other hand, while comparing the results obtained from free quadrilateral element and free triangular element, free quadrilateral element gives better results than that free triangular element. Hence quadrilateral element is preferred over triangular element. 4 noded shell 63 element gives very accurate results in comparison analytical results obtained from Kirchhf s Plate Theory. REFERENCES @IJMTER-2014, All rights Reserved 496
[1]. P. S. Gujar and K. B. Ladhane, Bending Analysis Simply Supported and Clamped Circular Plate, SSRG International Journal Civil Engineering, vol. 2, issue 5, pp.69-74, 2015. [2]. B. C. L. Vanam, M. Rajyalakshmi and R. Inala, Static Analysis An Isotropic Rectangular Plate using FEA, Journal Mechanical Engineering Research,vol. 4(4), pp.148-2, 2012. [3]. G. C. Mekalke, M. V. Kavadeand S. S. Deshpande, Analysis a Plate a Circular Hole by, IOSR Journal Mechanical and Civil Engineering (IOSR-JMCE), ISSN: 2278-84, pp.25-30, 2012. [4]. S. Boedo, Evaluation and comparison plate finite elements in ANSYS, NAS World Congress 99, Newport Rhode island, USA,1999. [5]. M. Vrabie and S. Baetu, Comparative Analysis the Bending Theories for Isotropic Plates- Case Study, BuletinulInstitutuluiPolitehnic Din Iasi, Publicat de UniversitateaTehnica, Gheorghe Asachi din IasiTomul LIX (LXIII), 2013. @IJMTER-2014, All rights Reserved 497