Document Type
Thesis
Date of Award
5-31-1986
Degree Name
Master of Science in Mechanical Engineering - (M.S.)
Department
Mechanical Engineering
First Advisor
Sachio Nakamura
Second Advisor
Bernard Koplik
Third Advisor
Benedict C. Sun
Abstract
The finite element formulation of axisymmetric plane stress analysis is used to develop a computer program for numerical implementation. The basic relations used to derive the finite element equations for a solid continum are reviewed in the thesis. The programming organization is summarized and program list is given in Appendix.
The program portion for isoparametric elements are borrowed from Numerical Methods in Finite Element Analysis by K. J. Bathe and E. L. Wilson, and various other subroutines are used from Mr. Huang's thesis on Micropolar elasticity analysis. Mr. Huang's computer program is based on the work done by Dr. Nakamura at New Jersey Institute of Technology as well as at The University of Iowa. The program is written in FORTRAN77 running under VAX11/VMS/780 system. Four-node quadrilateral isoparametric elements were used in the analysis. A U-grooved shaft was analysed for stress concentration, subjected to axial loading (tension).
Finally, the results obtained from the finite element solutions are compared with those obtained from exact solutions for a semi-circular U-grooved shaft. Different r/d values were taken, and its convergence to the exact graphical solution was investigated.
The accuracy of the results obtained were within 5-10% as compared with the exact solution. By increasing the number of elements in the analysis, the results could be further improved to achieve a higher degree of accuracy. Due to space and time limitations on the computer, fewer number of elements were used for the analysis.
Recommended Citation
Shetty, Ramanand R., "Isoparametric finite elements for axisymmetric solids" (1986). Theses. 3382.
https://digitalcommons.njit.edu/theses/3382
