PART THREE

BENDING OF BEAMS WITH NON-SYMMETRICAL CROSS SECTIONS

The majority of aircraft structural components consist of beams with non-symmetrical cross section acting in bending . For this reason an expression needs to be derived to allow for the determination of the stresses induced by bending moments to such sections.

SIGN CONVENTION AND NOTATION

Look at the Oxyz system of axis, with an arbitrary beam parallel to the z-axis:

Figure 8 :  Notation and sign convention for positive forces, moments and displacements.

Where :

T = Torque

M = Bending Moment

S = Shear Force

w = Distributed load

P = Axial or direct load

u,v,w = Axial displacements

All of these externally applied loads are positive in the direction indicated in the figure. Internal moment and forces applied to face A are in the same direction and sense as externally applied loads. However on face B, the positive internal moments and forces are in the opposite sense.

RESOLUTION OF BENDING MOMENTS

A bending moment M applied in any plane parallel to the z-axis can be resolved into Mx and My components by normal vector rules.

By doing it in a visual way it will be easier to see:

Figure 9:  Resolved bending moment about x and y axis.

From Figure 9, the following relationships can be obtained:

Mx = M sin θ

My = M cos θ

and that these moments can have different sign depending on the value of θ. For example if θ > π/2, Mx is positive and My is negative.

STRESS DISTRIBUTION DUE TO BENDING (ETB-NonSym)

Consider a beam of arbitrary non-symmetrical cross section, which supports bending moments Mx and My, bending about some axis in the cross section. This is the plane of no bending stress, called Neutral Axis (N.A.).

Figure 10:  Determination of Neutral Axis location.

Let the axis origin coincide with the centroid G of the cross section, and that the neutral axis is a distance p from G.

The direct stress sz on element dA at point (x,y) and distance x from the neutral axis is:

Look at the beam in a plane parallel to the neutral axis with two segments ij and kl which are of equal length when the beam is undeflected:

Figure 11:  Side view of undeflected beam with segments ij and kl marked.

Once the beam has been deflected this section will look like this:

Figure 12: Deflected beam.

where:

        R   =  the radius of curvature

        dq  =  angle between planes ik and jl

The strain in plane kl can be defined as:

with

and

giving

By substituting back into the stress equation it gives:

(3.2)

Now that a stress equation has been obtained, it is necessary to satisfy both rotational and linear equilibrium at the ends of the beam. That is, the summation of the forces through the depth of the beam is equal to '0', and the summation of the moments through the depth of the beam is equal the applied moments Mx and My.

As the beam supports pure bending, the resultant load on the end section must be zero. Hence

Substituting equation 3.2 gives:

This equation defines the location of the centroid of the section, it follows that the neutral axis must pass through the centroid.

Rather than use this equation to find the location of the centroid, it is much easier to locate the centroid about the xy-axis by using equation 3.3.

In order to have moment equilibrium, it is necessary to re-draw Figure 10 but with the axis passing through the centroid.

Figure 13:  Beam section with Neutral Axis passing through centroid.

In order to see this in more detail, Figure 14 shown a close up of the axis and the area dA.

Figure 14:  Detail of area dA in beam cross-section.

If the inclination of the Neutral Axis (N.A.) is at an angle from the x-axis then :

x = x sin α + y cos α (3.4)

and substituting into 3.2 gives:

(3.5)

The moment resultants have the same sense as the applied moments, hence:

, (3.6)

and substituting equation 3.4 into equations 3.5 gives:

and

Both the sin α and cos α terms are not a function of dA, so they can be removed from the integration. What remain are terms which only have to do with the characteristics of the cross sectional shapes of the beam, and these are just the 2nd moments of area of the beam.

The second moments of area about the xy axes are:

which gives:

solving simultaneously gives,

Substituting into Equation (3.5) gives:

By defining the terms Effective Bending Moment, as:

and

Equation 3.8 can be re-written as follows:

Note that if the beam is symmetrical about the x -axis or y-axis then:

and the x & y axes are the principal axes.

POSITION OF NEUTRAL AXIS

The location of the Neutral Axis was defined by Equation (3.3).

For a beam with a symmetrical cross section, the centroid is the point defined by equation 3.3 and the Neutral Axis is parallel to the x and y axis. For a non-symmetrical beam cross section however, the Neutral Axis passes at some angle a with respect to the x-axis. What needs to be done define the angle of the neutral axis.

At the N.A. the normal bending stresses are equal to sz = 0, giving that:

where:

        xN.A. and yN.A. are coordinates of points along the neutral axis, giving:

By taking the inverse tan of the angle 'a', the angle of the Neutral Axis with respect to the x-axis can be found, given by equation 3.12.

Example 1: The beam shown is subjected to a bending moment of 150Nm about the x-axis.                    Calculate the maximum direct stress due to bending stating where it acts.

Figure 15:  Beam cross section with applied bending moment.

a)   Determine the location of the centroid about point A

and

b)   Shift axis to centroid

Figur 16:  Beam cross section with axis shifted to centroid position.

c)   Determine other sectional properties

d)   Calculate effective bending moments

Since Mx = 1500 and My = 0, substituting these values into Equations (3.9) and (3.10) gives:

e)   Determine bending stress equation and cross section stresses

Substitute the values calculated above into Equation (3.11).

Substituting the coordinates of all the corner points from Figure 16 give that:

Point x (mm) y (mm) sz (MPa)
1 72 21.6 0.436
2 -48 21.6 5.1
3 -48 13.6 3.9
4 -16 13.6 2.66
5 -16 -66.4 -9.32
6 -8 -66.4 -9.63
7 -8 13.6 2.34
8 72 13.6 -0.76

You can use this method to determine the stresses due to bending in any type of beam with any type of cross section.

APPROXIMATION FOR THIN WALLED SECTIONS

Due to the thin nature of aircraft structures, the assumption can be made that stresses are constant through the thickness 't' of the skin. By saying this, the square and higher powers of 't' can be neglected from the computation of sectional properties.

To see this, look at a channel section:

Figure 17: Thin-walled channel section.

1)   Because the section is symmetrical about the x-axis, then Ixy = 0

2)  The second moment of area Ixx is given by:
Expanding the RHS, gives:

by eliminating t2 powers and higher it becomes:

and similarly

This indicates that the sectional properties may be calculated as if the section was represented by a thin line, as shown in Figure 18, disregarding any t2 or higher terms.

Figure 18: Approximation of channel section.

For this discrete thin walled section, the sectional properties would be found as follows:

1)    Determine location of centroid

2)    Shift axis to centroid

Figure 19: Axis at centroid

3)    Determining moments of inertia

These results are exactly the same as for the section considering the skin material thickness and then disregarding all t2 and higher terms.

Since not all skin sections will lie parallel to either the x or y axis, the local moments or area for a section of skin at an angle q with respect to the x-axis are given by the following equations.

Figure 20: Thin skin section inclined at an angle θ wrt the x-axis

GENERAL LOADING RELATIONSHIPS

Consider an element of length dz from a beam with an unsymmetrical cross section with all types of loads applied in the y-z plane, Figure 21.

Figure 21: Equilibrium of generally loaded beam element in zy-plane.

Equilibrium of element in y direction gives:

dividing by dz and in the limit as dz0, this equation simplifies to :

Taking moments about A gives:

dividing by dz and in the limit as dz0, this equation simplifies to :

Combining these two equations gives:

Similarly, about the x-z plane:

Consider now an element of length dz from a beam with an unsymmetrical cross section with only an applied torque about z-axis, Figure 22.

Figure 22: Torque equilibrium of beam section of length dz.

Taking moments about the z-axis:

dividing by δz and in the limit as δz0, this equation simplifies to :

By assuming that a parameter has the same relationship to as Sy has to Mx then by differentiating equations (3.9) and (3.10) produces:

Also the parameters and are related to the load intensities wy and wx in the same manner such that by differentiating Equations (3.16) and (3.17) gives:

The parameters ,, and are called the EFFECTIVE SHEAR FORCES and LOAD INTENSITIES.


...Back to Top