Strain energy form coefficients for bending of short beams having full and thin-walled arbitrary cross section with application for FEM and Airfoils

The aim of this work is to develop a new generalized formula and a numerical computation program for evaluating the energy form coefficient of a complex and arbitrary cross section for full and thin-walled cross section with respect to any central axis, for the bending of beams of small lengths in comparison with the transverse dimension of the section. This coefficient plays a very important role in the calculation of the deformation energy of beams subjected to bending under the effect of a shearing force for short beams. It also enters in the formulation of FEM bending model, in order to calculate the stresses and the strains due to the external forces. The application is made for complex sections used in various fields of construction and in particular for airfoils designed for aerospace construction. A method is developed to calculate this coefficient as a function of the rotation of the central axes. The calculation of the area, the moments, and the product of inertias with respect to the central axes is necessary. The formula for calculating this coefficient is presented as a definite integral of a non-analytical function determined point by point along the direction of the application of the shear force. This function is based on the calculation of the partial static moments. The calculation of the latter is based on the development of a technique by subdividing the upper part of the section into adjacent common triangles at one point for the full solid section or by segments on the boundary for the thin-walled section. To speed up the process of numerically calculating this integral with high precision and reduced time, Gauss Legendre quadrature of order 40 is used. The calculation of the distribution of the tangential stress as well as its maximum value is determined. A shear shape coefficient is therefore determined. In the second part of this work, an application is made for the static calculation by the FEM of a hyper static beam with a view to determining the influence of this coefficient on all the parameters of resistance and bending stiffness as a correction of the classical model of bending by the FEM. A study of the error made by the classical bending model on our shear effect model is presented. A coefficient of efficiency of a section is presented.