Sandwich Panel Flexure Stress Equation Derivation

Look at the equation in the picture above, did you ever wonder how this sandwich panel flexure stress equation derivation came about?

I will try to derive this equation in this post. Ready? Let’s begin…

Sandwich Panel Flexure Stress Equation Derivation:

The flexure stress equation is a derivative of the basic plate bending stress equation as shown in the figure above.

For a plate, it is simply $Mc/I$:

Where:

• ‘M’ is the applied moment (in-lb)
• ‘c’ is the extreme fiber distance from the neutral axis, where the bending stresses are the highest (in)
• ‘I’ is the effective area moment of inertia of the cross section ( $in^4$)

In case of sandwich panels, we start with the same equation. However, the main underlying assumptions are as follows:

1. Only the facings or the skins resist the lion’s share of the induced bending stresses
2. The core mainly keeps the skins apart so they cannot move relative to each other
3. The effective area of the skin is used in the area moment of inertia calculations along with the standard parallel axis theorem
4. Very small quantities in the area moment of inertia calculation are zeroed out
5. Skin thickness $t_f$ << Core $t_c$, usually in the order of 1/100 ths of an inch

1. The area of the facings or skins is $(t_f) * (b) in^2$
2. The effective thickness of the sandwich panel, which is the mid plane distance between the facings, is $h src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-=(t_f + t_c)$ $in$

Basic Bending Stress Equation: $S src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= Mc/I$ —> (1)

But now: $c src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= h/2$ —> (2)

Normally a rectangular I beam flange area moment of inertia $I src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= b*h^3/12+A*d^2$

But now:

b = panel width (in) $h src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= t_f$ (in)

d = h/2 (in), the mid plane distance of the skin from the neutral plane

Therefore for teh two skins combined: $A src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= b*t_f$ (each skin’s c/s area) $I src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= 2*b*(t_f)^3/12 + 2*b*t_f*(h/2)^2$ $I src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= b*(t_f)^3/6 + b*t_f*(h)^2/2$

However, the first part of the RHS side of the above equation becomes a very tiny number. The significance of this is that the skins themselves are basically incapable of resisting much bending as standalone parts of the layup, the crucial part is the separation that the core provides. Read more about it here.

For example:

If $t_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= 0.02"$, then $t_f^3/6 src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= 2.67E-8$ => 0 (approaches zero)

So the first part is zeroed out, hence:

Effective $I src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= b*t_f*h^2/2$ —> (3)

Substituting (2) and (3) into the basic bending stress equation (1): $S_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= Mc/I src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= M*(h/2)/(t_f*b*h^2/2)$ $S_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= M/(h*t_f*b)$ $S_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= M/(h*t_f*b)$ $S_f = M/(h*t_f*b)$ $S_f = M/(h*t_f*b)$

Alternative Approach:

Another simpler approach is to use the force couple approach.

All the known quantities such as M, h etc. are the same. But this time we will use a section cut at the middle of the panel.

In this method, we cut a section through the middle of the long beam test coupon panel.

A force couple or resisting torque is developed within the skins, acting in the plane of the skins. Thus induces compressive stresses in the top skin and tensile stresses in the bottom skin.

The force is resisted by the cross sectional area of the skins (tf*b), we can now use the standard S = P/A formula to calculate the facing stresses.

Therefore:

==> $S_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= M/(h*t_f*b)$ $S_f src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" data-= M/(h*t_f*b)$ $S_f = M/(h*t_f*b)$ $S_f = M/(h*t_f*b)$

So there you have it, hope you learnt something new and interesting today. Make sure you comment below the post and let me know what you think…

Be Confident in Your Job & Interviews, Enroll Today… Surya Batchu

Surya Batchu is the founder of Stress Ebook LLC. A senior stress engineer specializing in aerospace stress analysis and finite element analysis, Surya has close to a decade and a half of real world industry experience. He shares his expertise with you on this blog and the website via paid courses, so you can benefit from it and get ahead in your own career.