Determination of Optimal Stress Locations in 2D–9 Noded Element in Finite Element Technique

Authors: Nishant Shrivastava, D. K. Sehgal

Abstract:

In Finite Element Technique nodal stresses are calculated through displacement as nodes. In this process, the displacement calculated at nodes is sufficiently good enough but stresses calculated at nodes are not sufficiently accurate. Therefore, the accuracy in the stress computation in FEM models based on the displacement technique is obviously matter of concern for computational time in shape optimization of engineering problems. In the present work same is focused to find out unique points within the element as well as the boundary of the element so, that good accuracy in stress computation can be achieved. Generally, major optimal stress points are located in domain of the element some points have been also located at boundary of the element where stresses are fairly accurate as compared to nodal values. Then, it is subsequently concluded that there is an existence of unique points within the element, where stresses have higher accuracy than other points in the elements. Therefore, it is main aim is to evolve a generalized procedure for the determination of the optimal stress location inside the element as well as at the boundaries of the element and verify the same with results from numerical experimentation. The results of quadratic 9 noded serendipity elements are presented and the location of distinct optimal stress points is determined inside the element, as well as at the boundaries. The theoretical results indicate various optimal stress locations are in local coordinates at origin and at a distance of 0.577 in both directions from origin. Also, at the boundaries optimal stress locations are at the midpoints of the element boundary and the locations are at a distance of 0.577 from the origin in both directions. The above findings were verified through experimentation and findings were authenticated. For numerical experimentation five engineering problems were identified and the numerical results of 9-noded element were compared to those obtained by using the same order of 25-noded quadratic Lagrangian elements, which are considered as standard. Then root mean square errors are plotted with respect to various locations within the elements as well as the boundaries and conclusions were drawn. After numerical verification it is noted that in a 9-noded element, origin and locations at a distance of 0.577 from origin in both directions are the best sampling points for the stresses. It was also noted that stresses calculated within line at boundary enclosed by 0.577 midpoints are also very good and the error found is very less. When sampling points move away from these points, then it causes line zone error to increase rapidly. Thus, it is established that there are unique points at boundary of element where stresses are accurate, which can be utilized in solving various engineering problems and are also useful in shape optimizations.

Keywords: Finite element, Lagrangian, optimal stress location, serendipity.

References:

[1] Barlow, J., Optimal stress locations in Finite Element Model., Int. J. Numerical Methods in Engg., 10, 243-251 (1976).
[2] Barlow, J., Comments on optimal stress locations in Finite Element Modes., Int. J. Numerical Methods in Engg., 11, 604 (1977).
[3] Barlow, J., More on optimal stress points- Reduced integration, Element distortion and error estimation., Int. J. Numerical Methods in Engg., 28, 1487-1504 (1989).
[4] Budkowwska, B.B., and Fu. Q, Analystical determination of the optimal strain and stress points for the displacement model of Finite Element Method., Computer and Structure, 41,937-946 (1991).
[5] Demkowicz, L., and Oden, J.T., Extraction method for second derivative in Finite Element Approximation of Linear elasticity problem, Comm., App., Numer. Meth., 1, 137-139(1985).
[6] Hinton, E. and Owen, D.R.J., An Introduction to Finite Element Computation, Pineridge Press U.K. pp 131-134(1979).
[7] Iron B.M., Engineering application of numerical integration in stiffness method, AIAAJ, 4,2035-2037(1966).
[8] Iron B.M., Least Square surface fitting by Finite Element and an application in stress smoothing, Aero Stress memo., ASM, 1524, Roll Royce, (1967).
[9] Moan, T., On local distribution of errors by finite element approximations, Seminar on Finite Element Analysis, Tokyo, (1973).
[10] Oden, J., T., and Brauchil, H.,J., On the calculation of consistent stress distribution in finite element applications, Int. J. Numer. Engg. 3, 317-325(1971).
[11] Salama, M., and Utku, S., Stress computation in displacement method for two medium elastic media, Comp. Methods in App. Mech.and Engg., 10, 325-388(1977).
[12] Sehgal D.K., Some Further Work on Optimal Stress Location in Finite Elements, Proceeding of the International Conference on Computational Engineering Science pp 35 vii April 10-14, Atlanta, GA, U.S.A. (1988).
[13] Stein, E., and Ahmad, R., On the stress computation in Finite Element models based upon displacement approximations, Compu. Method in App. Mech. And Engg., 4 81-96, (1974).
[14] Strang, G. and Fix, G., An Analysis of Finite Element Method, Prentice-Hall, Englewood cliff, N.J. pp 168-169 (1973).
[15] Zienkiewicz, O.C., Taylor, R.L., and Too, J.M., Reduced Integration technique in general analysis of plates and shells., Int.,Numer. Meth Engg., 3, 275-290(1971).