A Sample of Calculation of Anisotropy Loaded Beams

Władysław Papacz, Piotr Kuryło, Edward Tertel

  Open Access OPEN ACCESS  Peer Reviewed PEER-REVIEWED

A Sample of Calculation of Anisotropy Loaded Beams

Władysław Papacz1,, Piotr Kuryło1, Edward Tertel1

1Mechanical Department, University of Zielona Góra, Zielona Góra, Poland


In this paper, newly-derived equations, leaning on the canons of the classic durability of materials are used for calculation of leaf springs. Analytic relations between displacements, deformations, Young’s modulus and external loads are used. The correctness of the method was verified through theoretical and experimental research. Researches were done on composite leaf spring made from polymer and glass fibres. The main value of the new theory is its high similarity with the experimental data.

At a glance: Figures

Cite this article:

  • Papacz, Władysław, Piotr Kuryło, and Edward Tertel. "A Sample of Calculation of Anisotropy Loaded Beams." American Journal of Mechanical Engineering 1.7 (2013): 300-303.
  • Papacz, W. , Kuryło, P. , & Tertel, E. (2013). A Sample of Calculation of Anisotropy Loaded Beams. American Journal of Mechanical Engineering, 1(7), 300-303.
  • Papacz, Władysław, Piotr Kuryło, and Edward Tertel. "A Sample of Calculation of Anisotropy Loaded Beams." American Journal of Mechanical Engineering 1, no. 7 (2013): 300-303.

Import into BibTeX Import into EndNote Import into RefMan Import into RefWorks

1. Introduction

In the paper [1, 2, 3, 4] the construction of left and new dependencies for large deflections of curved beams was described. New equations are basis of analysis of perpendicular dislocations in the range of small linear diffractions or large non-linear diffractions of straight beams. To calculate beams about any curvature were used analytic connections between dislocation of axis of beam, deformation of beam, surface of transversal section A and module of Young and external loads.

The analytical dependences between load F and moment of bend M for burdened beam will be shown in [2, 3]:

Ÿ deformations (geometrical relationship):


Ÿ normal strength (physical relationship):


Ÿ moment of bend (physical relationship):


where: U- radial dislocation, V- circuit dislocation axis of beam, - curvature of the beam, s - coordinate the length of the beam, F - normal force


where: r - the radius of curvature, φ- polar angle

In geometrical relationship (1.1) we can notice:


Equations (1.2, 1.3) can to be presented in form:


Dislocation of any point of axis of beam can be described by the equation:


and it increases:


where: U- radial dislocation, V- circuit dislocation axis of beam

Linear size of increase of dislocation equals:


Angle of turn () section of beam, contain is between unit vectors () contiguous to axis of spring Figure 1 and contiguous vector () to deformed axis. We describe this by equation:


And increase of angle of turn of section of beam calculates by equation:

Figure 1. Geometrical model of turn angle of crooked beam section

And then :


Deformation we can to describe by equation:


Geometrical relationship between deformation, dislocations and curvature as well as location of considered fiber describes equation:


Bending moment can be described by equation:


This equations of binding large dislocation with dimensions of section of beam, curvature and physical properties.

2. Application of the Method

In the next part of this article the derived equations were applied to calculations of composite lift spring. Internal burdens of any element of spring we can count use well-known conditions. Sum of projections of strengths in studied section Figure 2:


and from this:


Sum of moments:


Equation 2.3, 2.4, 2.5 and 1.16, 1.17 give possibility to count radial dislocations (U) and circumferential (V) of spring axis in dependence from external loads Figure 2.

Relationships between radial dislocations (U) and circumferential (V) and vertical dislocations (f) and horizontal (fp) are was possible to get from equations:

Horizontal dislocations - fp


Vertical dislocations – f


If we'll put to equation of inertia moment of spring section (1.7):


and we’ll accept, that:


We will get possibility of calculation of leaf spring at any shape.

For spring at solid widths - and variable thickness - to have substituted (2.8) and (2.6) to the equation (1.2) we receive maximum squeezing tensions:


3. Experimental Research

The equations presented above were used to calculate the composite curved leaf spring. In the further step the results of calculations were compared with the FEM calculations. The leaf spring has a regular width and a changing thickness. The scheme of leaf spring was shown on the Figure 3. It has the three-ply construction. The core was made from glass epoxy composite reinforced by textile. External layers have regular thickness, they are reinforced by continuous glass fibers. The width of layers is equal to the width of leaf spring.

Figure 3. Geometrical dimensions of composite leaf spring

The width of composite leaf spring was 80 mm and its length was 1110 mm. The changes of the thickness in the function of his length were introduced in Table 1. Signs used in the table can be found in Figure 3.

Table 1. The radial thickness of composite leaf spring in the function of the his length

Using the Algor system the calculations of composite leaf spring were conducted. The results of FEM were compared with the calculations made using the new author's analytic method and with the results of laboratory research. Leaf springs were tested on the static load. The perpendicular and horizontal dislocation of center leaf was measured on various loading. The experimental data was compared with the results of calculations obtained using the theory of small deflections (TMU), the method FEM (MES) and the theory of large deflections (TDU). In Figure 4 the horizontal dislocation of the central part of one composite leaf spring was shown. Experimental results were marked by symbol “EXP”.

Figure 4. The perpendicular dislocation of central part of composite leaf spring

Similarly, the horizontal dislocation of central part of leaf spring was put together in Figure 5.

Figure 5. The horizontal dislocation of central part of composite leaf spring

In case of the perpendicular dislocation, exactitudes of particular computational methods are similar. The dislocation calculated by the FEM method and by the method derived in the support of the theory of large dislocation (TDU) are approximate to measured. The horizontal displacements calculated by TDU better describes the actual situation. In Figure 4 and Figure 5 the line marked EXP measurement data indicated a line described TDU symbol calculation results obtained using the derived equations. It should be noted much less time-consuming the TDU method in comparison to the method of MES. The methodology described above, has been used to calculate stress in the selected cross-section spring. In Figure 6 the normal tension in external layers of spring was shown.

Figure 6. Normal tension exterior cladding of leaf spring

Methods give very similar sizes of tensions. In sections No. 1, 2 and 3 the differences between methods did not exceed eight percent. In faintly burdened sections (4 and 5) the differences are larger and achieve eighteen percent. Counted tensions value without help of a new method (TDU) are comparable to counted values of FEM method.

4. Conclusions

The equations 1.16 and 1.17 can be used to calculate large dislocations of beams with any curvature, composite beams as well.

The newly-introduced correlations have been based on the classical canons of strength, assuming the continuity of the structure of construction materials, in contrast to finite element methods. This tool enables verification of the correctness of the calculations obtained by commonly used MES methods. It can also be used as a standalone tool for strength calculations of anisotropic structures. The correctness of the method was verified by theoretical and experimental studies. Tests were carried out on leaf spring made from epoxy resin and glass fiber. Formed by the FEM methods the composite spring were analyzed using analytical methods based on the new author theory. A series of laboratory research have been conducted as experimental studies. Theoretical studies have shown the usefulness of the newly-developed theory as a tool for verifying the correctness of FEM calculations. Maximum differences in treatment outcomes between the FEM and the developed theory have not exceeded 20%. The error of theory relative to the experimental measurements have not exceeded 18%. The main value of the new theory is its high similarity with the experimental data. The theory by no means can replace the FEM methods. But, it can enable to monitor the correctness of the FEM calculations performed. The theory can also be used to create rods design. It should be noted that the analytical method is based on the new theory is much less time-consuming compared to the MES. The simplicity of the method is its big advantage.


[1]  Papacz W.: Some parameters of composite leaf springs. Int. J. Applied Mechanics and Engineering. 2009, vol. 14, pp1181-1187.
In article      
[2]  Bjorklung A.: Contribution to the theory of elastic rings into account the deformation - Ing Arch 23 (1956).
In article      
[3]  Władysław J. Papacz, A Method for Calculation of Anisotropy Loaded Beams /Procedia Engineering, Elsevier, Volume 48, 2012, Pages 535-532, Modelling of Mechanical and Mechatronics Systems.
In article      
[4]  Bocko, J., Frankovský, P., Kostelníková, A., Ostertagová, E., Ostertag, O.: Structural design and photoelasticimetric verification of landing gear of ultralight aircraft. In: Metalurgija. Vol. 49, no. 2 (2010), p. 145-150.
In article      
  • CiteULikeCiteULike
  • MendeleyMendeley
  • StumbleUponStumbleUpon
  • Add to DeliciousDelicious
  • FacebookFacebook
  • TwitterTwitter
  • LinkedInLinkedIn