In his most important point Professor Stephan is completely correct. That is, his shear coefficient and mine are identical. I was aware of his work when it first came out and discounted it because it was derived by solving a static problem for a specific type of loading. It implied that the shear coefficient for the static problem was a function of the type of loading and, further, that by choosing the right type of static loading one could find the best value of the shear coefficient for dynamic loads. As it turns out both of these implications are correct.

Professor Stephan states, “The motivation of Hutchinson appears to be the construction of a theory in which the shear coefficient takes on the “best” value. I did not mean to convey that impression. My motivation was simply to construct a simple, consistent, dynamic theory which did not require guessing a shear coefficient. This simple consistent theory allowed me to find an expression for the shear coefficient in the Timoshenko beam theory. That this shear coefficient agreed with the “best” values simply validated my approach.

Professor Stephen states, “the values of the coefficient for the circular cross section, both solid, hollow and thin-walled, and for the elliptic cross section calculated in [1], are identical to those given in [3,4].” What he doesn’t note is that the expressions for the rectangular cross section in both his and my paper also produce identical results. He further comments on the fact that the shear coefficient goes to zero for the rectangular cross section. Actually, as shown in my Fig. 4, it is the reciprocal of the shear coefficient that goes to zero which means the shear coefficient would have a pole at that point. Professor Stephens conclusion that the beam is stiffened with an increase in width is correct. I do not understand, however, his remark that, “one would not normally employ Timoshenko theory for a beam having large width to depth ratio.” For a Poisson’s ratio of 0.3 the reciprocal of the shear coefficient goes to zero at a width to depth ratio of about 3. This is definitely within the range I would expect Timoshenko theory to be applied. Also the lowest natural frequency would probably occur (depending on boundary conditions) about a neutral axis in the width direction. I should probably also note that the shear coefficient is very sensitive to the assumed shear stress distribution. If one were to assume a simple parabolic shear stress distribution through the thickness of the rectangular beam one would get a shear stress of 51+ν/6+5ν) independent of the aspect ratio.

In my paper I also had a small section on static problems. In that section I found expressions for a static shear coefficient found by comparing solution of my equations with the elementary beam solution including shear deformation. The comparison was done for a tip loaded cantilever. Since Professor Stephens work indicates that the coefficient depends on loading I decided to investigate a beam loaded under it’s own weight. To accomplish this I inserted a new term γv into the integral in my Eq. (28) and dropped all the time dependent terms. The term γ is the specific weight of the beam and v is the displacement in the y-direction. Proceeding in the same way as in the paper, I came up with the following set of equations,  
ψ+γAEIz=0
(1)
 
ψφ+C4Izψγν2EIzIzIy=0.
(2)
These equations are solved by integrating the first three times and inserting the result into the second and integrating two more times. The boundary conditions for a cantilever fixed at x=0 and free at x=L are ψ0=0,φ0=0,ψL=0 and ψL=0. Since there are five constants of integration another condition is needed. That condition is, for the axial location at which the shear is zero, the slope ψ equals the slope of the center line deflection φ. Thus ψL=φL. The expression for the center line displacement is then  
φ=γAEIzx424+Lx36L2x24+γAEIz2c4+γν2EIzIzIyLxx22.
(3)
The first group of terms in Eq. (3) corresponds to the Euler-Bernoulli beam solution and the second group of terms corresponds to the shear deformation solution. This is the same expression that would result from simple beam theory with shear deformation if the shear coefficient were the expression shown in my paper as Eq. (57). If one uses the integrated average displacement, as was done by Professor Stephen, instead of the center line displacement then one gets the shear coefficient which he found, that is, my Eq. (41) and his comment Eqs. (1) and (3). Thus, even in my approach the shear coefficient is a function of the loading, and for a gravity load the resulting shear coefficient is the same as the dynamic coefficient. Similar Eqs. to (1) and (2) above could be developed for any type of loading on a beam and thus eliminate the need for a shear coefficient entirely.

I recently presented a paper entitled “Shear Coefficients for Thin-Walled Timoshenko Beams” at the Third International Symposium on the Vibrations of Continuous Systems, July 23–27, 2001 at Jackson Lake Lodge, WY. In that paper I considered all the thin-walled cases treated by Cowper plus two additional cases. As in all other cases the “best” shear coefficient agreed with Cowper’s values only for Poisson’s ratio equal zero.

As to Professor Stephen’s final remark, “it is noted that while the above values for the coefficient may be widely accepted as best, paradoxically Cowper’s values appear to be more widely used; it is hoped the investigators will in the future make greater use of these “best” values,” I fully agree.