DETERMINATION OF EI FOR PULTRUDED GFRP SHEET PILE PANELS. Abstract

Similar documents
Purpose of this Guide: To thoroughly prepare students for the exact types of problems that will be on Exam 3.

= 50 ksi. The maximum beam deflection Δ max is not = R B. = 30 kips. Notes for Strength of Materials, ET 200

External Pressure... Thermal Expansion in un-restrained pipeline... The critical (buckling) pressure is calculated as follows:

FIXED BEAMS IN BENDING

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Experimental Lab. Principles of Superposition

TORSION INCLUDING WARPING OF OPEN SECTIONS (I, C, Z, T AND L SHAPES)

NAME: Given Formulae: Law of Cosines: Law of Sines:

Chapter 4-b Axially Loaded Members

DESIGN OF BEAMS AND SHAFTS

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method


Assumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

Application of Finite Element Method to Create Animated Simulation of Beam Analysis for the Course of Mechanics of Materials

Multi Linear Elastic and Plastic Link in SAP2000

3. BEAMS: STRAIN, STRESS, DEFLECTIONS

Entrance exam Master Course

Introduction to Structural Member Properties

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF MECHANICAL ENGINEERING QUESTION BANK. Subject code/name: ME2254/STRENGTH OF MATERIALS Year/Sem:II / IV

IDE 110 Mechanics of Materials Spring 2006 Final Examination FOR GRADING ONLY

Journal of Composite Materials. Behaviour of hollow pultruded GFRP square beams with different shear span-to-depth ratios

2. Determine the deflection at C of the beam given in fig below. Use principal of virtual work. W L/2 B A L C

PES Institute of Technology

BEAM DEFLECTION THE ELASTIC CURVE

Homework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. Fall 2004

Experimental Study and Numerical Simulation on Steel Plate Girders With Deep Section

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.

CHARACTERIZING INFLUENCE OF LAMINATE CHARACTERISTICS ON ELASTIC PROPERTIES OF CROSS LAMINATED TIMBER

Chapter 7: Bending and Shear in Simple Beams

Structural Analysis. For. Civil Engineering.

CIVIL DEPARTMENT MECHANICS OF STRUCTURES- ASSIGNMENT NO 1. Brach: CE YEAR:

International Journal of Advanced Engineering Technology E-ISSN

By Dr. Mohammed Ramidh

STRENGTH OF MATERIALS-I. Unit-1. Simple stresses and strains

5. What is the moment of inertia about the x - x axis of the rectangular beam shown?

EFFECT OF RETROFITTING ON FLEXURAL RIGIDITY OF RECTANGULAR FRP PLATES

QUESTION BANK SEMESTER: III SUBJECT NAME: MECHANICS OF SOLIDS

Influence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Bending Load & Calibration Module

[5] Stress and Strain

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

CHAPTER 4. Stresses in Beams

Module 4 : Deflection of Structures Lecture 4 : Strain Energy Method

Safe-T-Span High Load Capacity Grating

Structural Steelwork Eurocodes Development of A Trans-national Approach

Laboratory 4 Topic: Buckling

structural analysis Excessive beam deflection can be seen as a mode of failure.

Karbala University College of Engineering Department of Civil Eng. Lecturer: Dr. Jawad T. Abodi

UNIT II 1. Sketch qualitatively the influence line for shear at D for the beam [M/J-15]

QUESTION BANK DEPARTMENT: CIVIL SEMESTER: III SUBJECT CODE: CE2201 SUBJECT NAME: MECHANICS OF SOLIDS UNIT 1- STRESS AND STRAIN PART A

CHARACTERIZATION OF ROTATIONAL BEHAVIOR OF WEB-FLANGE JUNCTIONS OF PULTRUDED GFRP BRIDGE DECKS

PLAT DAN CANGKANG (TKS 4219)

BEAM A horizontal or inclined structural member that is designed to resist forces acting to its axis is called a beam

Figure 1. Simply-supported beam nomenclature with single point load. Typical shear and moment diagrams are also shown.

Chapter 11. Displacement Method of Analysis Slope Deflection Method

National Exams May 2015

DIVISION: METALS SECTION: METAL FASTENINGS SECTION: STEEL DECKING REPORT HOLDER: PNEUTEK, INC.

Software Verification

EAS 664/4 Principle Structural Design

UNIVERSITY OF BOLTON SCHOOL OF ENGINEERING. BEng (HONS) CIVIL ENGINEERING SEMESTER 1 EXAMINATION 2016/2017 MATHEMATICS & STRUCTURAL ANALYSIS

DETERMINATION OF YOUNG S MODULUS FOR CFRP USING THE THREE POINT BENDING TEST

MEASURMENT OF BENDING PROPERTIES OF RHIZOPHORA TROPICAL WOOD. Sarawak, Malaysia. Engineering & Technology, Dhaka, Bangladesh

ISHIK UNIVERSITY DEPARTMENT OF MECHATRONICS ENGINEERING

Module 3. Analysis of Statically Indeterminate Structures by the Displacement Method

FLEXURAL RESPONSE OF FIBER RENFORCED PLASTIC DECKS USING HIGHER-ORDER SHEAR DEFORMABLE PLATE THEORY

Equivalent Uniform Moment Factor for Lateral Torsional Buckling of Steel Beams

MODULE C: COMPRESSION MEMBERS

1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine Spring 2004 LABORATORY ASSIGNMENT NUMBER 6

LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS

JUT!SI I I I TO BE RETURNED AT THE END OF EXAMINATION. THIS PAPER MUST NOT BE REMOVED FROM THE EXAM CENTRE. SURNAME: FIRST NAME: STUDENT NUMBER:

CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions


Rigid and Braced Frames

QUESTION BANK. SEMESTER: V SUBJECT CODE / Name: CE 6501 / STRUCTURAL ANALYSIS-I

CHAPTER II EXPERIMENTAL INVESTIGATION

Finite Element Modelling with Plastic Hinges

Physics 8 Monday, November 20, 2017

Physics 8 Monday, November 23, 2015

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.

Chapter Objectives. Design a beam to resist both bendingand shear loads

Due Tuesday, September 21 st, 12:00 midnight

Name :. Roll No. :... Invigilator s Signature :.. CS/B.TECH (CE-NEW)/SEM-3/CE-301/ SOLID MECHANICS

Problem d d d B C E D. 0.8d. Additional lecturebook examples 29 ME 323

Evaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis

Finite element analysis of FRP pipelines dynamic stability

14. *14.8 CASTIGLIANO S THEOREM

8 Deflectionmax. = 5WL 3 384EI

Name (Print) ME Mechanics of Materials Exam # 2 Date: March 29, 2017 Time: 8:00 10:00 PM - Location: WTHR 200

3.5 Reinforced Concrete Section Properties

SIMPLY SUPPORTED STRUCTURAL BEAM STRESS AND DEFLECTION ANAL

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.


Physics 8 Wednesday, November 29, 2017

Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test

Experiment Two (2) Torsional testing of Circular Shafts

Materials: engineering, science, processing and design, 2nd edition Copyright (c)2010 Michael Ashby, Hugh Shercliff, David Cebon.

EVALUATION OF BENDING LOAD IN BATTER PILES SET IN SOFT CLAY

Part 1 is to be completed without notes, beam tables or a calculator. DO NOT turn Part 2 over until you have completed and turned in Part 1.

Transcription:

DETERMINATION OF EI FOR PULTRUDED GFRP SHEET PILE PANELS Yixin Shao, Cynthia Giroux and Zeid Bdeir McGill University Montreal, Quebec, Canada Abstract The flexural rigidity, EI, plays an especially important role in fiber reinforced polymer composite sheet pile wall design. Determination of EI is usually affected by the presence of the shear either due to the short testing span or to the low shear modulus of composites. To reduce or eliminate the effect of shear on the determination of flexural rigidity of composite sheet pile panels, this paper is to propose a test method that is based on Bank s multi-span approach and Timoshenko s beam theory. Tests were performed using three-point and four-point bending set-ups, with single composite sheet pile panel at six different spans, and with connected two panels at eleven spans. The results showed that six EIs determined independently from different set-ups, at different deflection points along the beam and with different cross sections gave values that differed from each other by less than 3%. The EI so determined was a constant independent from the test span, the test set-up, the location of the deflection point on the beam, or the number of data points used in linear regression. Apparent EIs calculated from single span tests were strongly span-dependent and approached asymptotically to the value determined by multi-span approach. It was likely that the multi-span approach yielded the highest possible flexural rigidity that could be measured. The shear effect was therefore eliminated. Introduction The pultruded fibreglass reinforced polyester composite sheet pile wall is a unique and one of its kind waterfront structure. The corrosion resistance and high strength-to-weight ratio of composites are the primary characteristics supporting this choice of material. The use of composites also features low maintenance, reduced life-cycle costs, dimensional stability, and near-zero environmental toxicity. Composite sheet piling products have been currently used to a limited degree, or experimentally for light retaining structures along the waterfront and the coastline. Because of the lack of the field application history, the height of composite sheet piling is presently limited to 2-5 meters. Experimental projects included renovation of decayed timber walls as well as construction of new retaining walls. Standard procedure of installation is followed without any special equipment required. The important parameter in sheet pile wall design is the flexural rigidity, EI, especially for FRP composite sheet pile wall. The wooden plank sheet pilings are about 50 mm 300 mm solid in cross section. For precast concrete sheet piles, the cross sections range from 500 to 800 mm wide and 50 to 250 mm thick. Steel sheet piles in the United States and Canada, have thickness about 0-3 mm. Sheet pile sections may be Z, deep arch, low arch, or straight web sections. Compared to the traditional pilings, the composite sheet panels exhibit relatively low flexural rigidity owing to their thin walled structure (3-5mm thick) and low Young s modulus of panel members (E = 30-40 GPa). To help guide the proper use of composite piling products, US Army Corps of Engineers had set a number of specifications to classify products of different duties (Lampo, et al. 998). The target flexural rigidities, EI, for different grades were proposed in a demonstration project for pultruded composite sheet pilings:

Light-duty: 5x03-x04 to x04-5x04 kip-sq in/ft Medium-duty: 5x04-x05 to x05-5x05 kip-sq in/ft Heavy-duty: 5x05-x06 to x06-5.5x06 kip-sq in/ft To check with the specifications for a given product, it is essential to properly characterize the EI of the product. Determination of EI is usually affected by the presence of the shear. To avoid shear effect on the EI, ASTM D-790 recommended the span to depth ratio of the beam tests be at least 6:. For the sheet pile walls made of traditional materials, the deflections due to shear might be negligible if the span is large enough. However, for composite materials, in addition to short span effect, the low shear modulus (G) could also cause the shear deflection to be significant. This paper is to report an experimental study on determination of flexural rigidity of the pultruded composite sheet pile panels with reduced or eliminated shear effect. Bank s approach was adopted in three-point and four-point bending tests with various spans. Bank (989) developed a method to simultaneously determine the section Young s modulus and the section shear modulus by using Timoshenko s deflection equation. The approach was modified in this paper to directly determine the EI, instead of E alone. Results from the three-point bending tests were compared with that from four-point bending tests to check the consistency, so were the results from single panel tests compared with connected double panel tests. The dependence of the flexural rigidity on the test conditions was also discussed. Deflection Equation Equation () gives the general form of Timoshenko s beam deflection theory. Because Timoshenko s theory uses superposition of the bending and shear deflections, it is limited to the linear elastic range of the material. 3 CPL C2PL δ = + () ( EI) ( kag) Here, the deflection, δ, depends on the applied load, P, the span length, L, the flexural rigidity, EI, and the shear rigidity, kag. The constants, C and C 2, depend on the load case and on the particular location along the span length that the deflection is desired. Values for C and C 2 are given in Table for the cases of three-point bending and four-point bending at one-third span. Note that the flexural rigidity is defined as the product of the section flexural modulus, E s, and the second moment of area, I, and the shear rigidity as the product of the factored area, ka, and the section shear modulus, G s. The factor, k, is the shear coefficient for composite beams in flexure. Table : Constants used in Timoshenko s equation Centre-point deflection (δ x=l/2 ) Quarter-point deflection (δ x=l/4 ) C C 2 C C 2 3-point bending 48 4 768 8 4-point bending 23 29 296 6 2304 8 2

δ Equation () can be rearranged into Equation (2) which is in the form of a straight line with PL 2 as the dependent variable and L as the independent variable. Knowing the slope of the load-deflection curve, δ P, and the span, L, of each beam tested, this straight line can be plotted. The slope and intercept of the line are inversely proportional to the flexural rigidity, EI, and the shear rigidity, kag, respectively. Theoretically, Equation 2 can be used to simultaneously determine the EI and kag. This paper will only discuss the possibility of obtaining EI with reduced or eliminated shear effect by conducting multi-span tests with minimum number of tests. δ = PL C L EI 2 + C2 kag (2) Experimental Program Single Panel Tests The composite sheet piles used in this study was designed and manufactured by Pultronex Corporation (Nisku, Alberta). The corrugated profile of a single panel has a symmetric double Z cross section. The panels are approximately 2.7cm deep, 42.5cm wide, and have a thickness varying from 0.32cm to 0.47cm. The panels interlock with each other through the pin and eye connections on both ends to form a continuous wall (Figures and 2). The pultruded profile consists of layers of E-glass rovings and mats in a polyester matrix. Tests were performed first with single panel in three and four-point bending tests (Figure ) on spans of length 0.9,.52, 2.3, 3.05, 4.57, and 6.0 meters, (3, 5, 7, 0, 5, and 20 ft) within the linear proportional limit of the material. Tests on each span were repeated a minimum of three times. The beams were loaded using an MTS testing machine with an average loading rate of 0.03mm/s. The deflections at midpoint as well as at one quarter point along the span were measured using two linear variable differential transducers (LVDTs) with ranges of 0cm. The data were recorded using a Measurements Group System 5000 data acquisition system using the Strain Smart software. For fourpoint bending tests, the load points were equally spaced at one-third span. Rigid steel frames of square tubes (25 mm x 25 mm x 2 mm thick) were constructed and positioned at the supporting and at the loading points of the beam in order to distribute the applied load and provide lateral confinement on the open section of the sheet pile (Figure ). Connected Double Panel Tests To examine if the connected panels had different flexural rigidity, tests of connected double panels were also performed using the four point bending setup at equal space with eleven spans,.5,.8, 2., 2.4, 3., 3.7, 4.3, 4.9, 5.5, 6., and 6.7 meters (5, 6, 7, 8, 0, 2, 4, 6, 8, 20, and 2.8 ft), and within the linear proportional limit. Each span test was repeated a minimum of three times. A pin and an eye were inserted at both ends of the connected double panels to simulate the connected wall as shown in Figure 2. Steel straps of size 2.5 cm ( inch) wide and mm thick were used to tie around the cross section, in order to confine the panels and to limit the lateral displacement. The use of flexible steel strapping at a uniform spacing of 6 cm (2 feet) was to provide a uniform restrain and prevent the panels from being stiffened by rigid heavy steel frame. 3

P LVDTs (a) Three-point bending set-up for single panel tests P LVDTs (b) Four-point bending set-up for both single and double panel tests Figure : Flexural test set-ups Steel Strapping Strain Gauge 2 Strain Gauge LVDT 2 LVDT Figure 2: Cross section of connected double panels 4

Results Single Panel Tests Typical curves within the linear proportional limit for the load versus midpoint displacement for the six spans are given in Figure 3 for the four-point bending tests. Linear fits using the method of least squares were performed on the load-deflection curves in Figure 3. Four sets data of load versus deflection were recorded in different set-ups (3-point and 4-point) and at different locations along the beam (midpoint and quarter point). The linear fit provided values for the slope, δ P, for each span, L, δ which allowed PL and L 2 to be calculated for use in the rearranged Timoshenko beam theory equation (Equation 2). δ For the deflection measured at the midpoint, the results of plotting this data with PL as the dependent and L 2 as the independent variables are plotted in Figure 4 for the four-point bending tests for all six span tests. For each span, at least three independent tests were performed and thus three sets of δ data on the versus L 2 were plotted in the same graph. As is seen, good repetition was achieved. PL Similar data process was carried out for three-point bending tests. Load [kn] 6 5 4 3 2 0 3' span 5' span 7' span 0' span 5' span 20' span 0 0 20 30 Midpoint Deflection [mm] Figure 3: Typical load-deflection curves at the midpoint for four-point bending linear tests (single panel, ' = 0.3048m) Linear fit was then carried out again using the least square method to determine the equation of the line in Figure 4. Based on Equation 2, the slope of the line was proportional to C /EI and could be used to calculate the flexural rigidity for the case. Tests and data process were repeated for other three cases and results are presented in Table 2. 5

δ 4.5E-06 4.0E-06 3.5E-06 3.0E-06 2.5E-06 2.0E-06.5E-06.0E-06 5.0E-07 0.0E+00 y = 8.85E-08x +.9E-07 0 0 20 30 40 50 L 2 [m 2 ] Fig. 4: Straight line used to determine the flexural rigidity from four-point bending based on the midpoint deflection (single panel) Connected Double Panel Tests Similar procedure was followed for connected double panel tests. The typical load deflection curves for eleven span tests are shown in Figure 5 for deflection measured from LVDT at midpoint. Slope of each linear test was determined by averaging readings of LVDT and LVDT2 at midpoint. δ/pl versus L 2 was then plotted in Figure 6 for eleven different spans. The linear fit of all data points was performed to find equation of the straight line. Based on equation 2, the slope could be used to calculate the flexural rigidity of the connected panels. The same procedure was followed to obtain the EI from the quarter point load- deflection data. The final results were divided by two to obtain the rigidity for single panel and also presented in Table 2 for comparison. 9.5 m 8.8 m 7 2. m Load (KN) 6 5 4 3 2.4 m 3. m 3.7 m 4.3 m 4.9 m 2 5.5 m 6. m 6.7 m 0 0 5 0 5 20 25 Mid point deflection, LVDT (mm) Figure 5: Typical load deflection curves at the midpoint for four-point bending linear tests (LVDT, connected double panels) 6

2.50E-06 2.00E-06 y = 4.425E-08x +.35E-07 R 2 = 9.934E-0.50E-06 δ/pl(n - ).00E-06 5.00E-07 0.00E+00 0 5 0 5 20 25 30 35 40 45 50 L 2 (m 2 ) Fig. 6: Straight line used to determine the flexural rigidity from four-point bending based on the midpoint deflection (connected double panels) The values obtained for the flexural rigidity agreed with each other for both three and four-point bending tests, both midpoint and quarter-point data, and both single and double panels. The average was calculated with a standard deviation of 2.9%. These self-checks are strongly supportive of the proposed experimental method. Table 2: Flexural rigidity of single panel determined by multi-span test method Midpoint deflection (δ x=l/2 ) Quarter-point deflection (δ x=l/4 ) Difference [%] EI (single panel) 3-point bending EI [knm 2 ] 22 26.8% 4-point bending EI [knm 2 ] 20 200 0.5% EI (two connected panels)/2 4-point bending EI [knm 2 ] 20 206 2.4% Average EI [knm 2 ] 206±6 2.9% Apparent Flexural Rigidity The apparent flexural rigidity, (EI) a, is defined as the resistance of a beam to deflection due to only bending, neglecting shear deflections. The slopes of the load-deflection curves were used to determine the apparent flexural rigidity based on the first term of Timoshenko s equation. They are presented in Figure 7 as (EI) a vs span to depth ratio. Obviously, the apparent flexural rigidity varied with the span to depth ratio, indicating that shear did have effect on flexural rigidity obtained from single span tests. In both curves (mid and quarter points in Figure 7), the apparent flexural rigidity increased 7

with the increase of span length. This suggested that a lager deflection contribution from the shear was present with shorter spans compared to longer ones. The apparent flexural rigidity asymptotically approached to 400KNm 2, the EI determined by two-term Timoshenko s Equation. It was likely that EI defined by Equation 2 was the largest possible rigidity that could be experimentally measured. Since Equation 2 separates the bending and shear contributions to the total deflection, the effect of shear on the EI is actually eliminated. ASTM D 790 suggests a span to depth ratio of 6/ or higher be used in beam tests to avoid shear effect. At this suggested ratio, the apparent flexural rigidity of the sheet pile panel counted to only 44% of the maximum possible value. 450000 400000 350000 Midpoint Quarter point 300000 EI (N m 2 ) 250000 200000 50000 00000 50000 0 0 0 20 30 40 50 60 Span to Depth Ratio Fig. 7: Apparent flexural rigidities versus span-to-depth ratio (connected panels) Sensitivity of EI to Data Process It was noticed during data processing that flexural rigidity (EI) was insensitive to the span and the number of data points used in the linear fit. The flexural rigidity of connected double panels is the result of tests of spans, from.5 m to 6.7 m, the last two data points shown in Figure 8 with all spans. Effort was made to examine if only a few small span tests would be enough to generate the comparable EI. The results are shown in Figure 8, in which EI is plotted versus the largest span included in linear fit. For instance, the first set of EI (from both midpoint and quarter points) was computed by linearly fitting test data of three smallest spans:.5,.8, and 2. meters. The largest span included in the linear fit was 2. m. Thus, EI was plotted against 2.m. Subsequently, the rest used test data of four spans with the largest 2.4m, of five spans with the largest 2.7m, and so on until the last set of all spans with the largest 6.7m. It was interesting to notice that variation of EI with test span and numbers of data in linear fit was not significant. In the other words, by using multi-span test method, EI could be determined consistently even with only three small-span tests. 8

500000 EI (Nm 2 ) 400000 300000 200000 Three smallest spans All spans 00000 0 Midpoint EI Quarter point EI 0 2 3 4 5 6 7 8 Largest Span Included in Linear Fit (m) Figure 8: Sensitivity of EI on largest span included in linear fit (connected panels) Discussion and Conclusions The flexural rigidities were determined for the fiberglass composite sheet pile panels using a modified Bank s multi-span approach. Values were found by performing both three and four-point bending tests, by taking deflection measurements at the midpoint and at the quarter-points of the span, as well as by using single panel and connected double panels. Six flexural rigidity results from independent tests agreed well with a standard deviation of 2.9%. The frames used for confinement did not affect the results. It seemed proper to group EI together as section parameter for composite sheet pile wall design. The flexural rigidity determined by the proposed multi-span method was quite different from that measured by the traditional single span method unless the span to depth ratio was over 60. The multi-span method gave the largest possible EI value with eliminated shear effect. EI so determined was a section constant independent from span length and numbers of tests, as long as minimum of three tests of different spans were performed. The average EI of the single panel was about 206 knm 2. In US customary units, it was corresponding to 5.47 x 0 4 kip-sq in/ft and the product was graded in the low end of medium duty according to specifications set by US Army Corps of Engineers. The same product however could be classified as light duty if only single span tests were performed. Therefore, standardization of the test method for composite sheet pile panels is necessary to make fair comparison for different products. The multi-span test method is strongly recommended in this regard. 9

Acknowledgment Financial supports from Natural Science and Engineering Research Council (NSERC) of Canada and from Pultronex Corporation of Alberta are gratefully acknowledged. References. Bank, L. C. (989). Flexural and Shear Moduli of Full-Section Fiber Reinforced Plastic (FRP) Pultruded Beams. Journal of Testing and Evaluation, 7(), pp. 40-45. 2. Lampo, R., et al. (998). Development and Demonstration of FRP Composite Fender, Loadbearing, and Sheet Piling Systems. US Army Corps of Engineers Construction Engineering Research Laboratory Technical Report 98/23. September 998. 3. Timoshenko, S. (955). Strength of Materials. New York: Litton Educational Publishing, Inc. pp. 70-75. 0