Dynamic characteristics of an FRP deck bridge

Similar documents
Transverse Distribution Calculation and Analysis of Strengthened Yingjing Bridge

Load Analysis and Multi Body Dynamics Analysis of Connecting Rod in Single Cylinder 4 Stroke Engine

Parametric study on behaviour of box girder bridges using CSi Bridge

STATIC AND FATIGUE ANALYSIS OF LEAF SPRING-AS A REVIEW

Influence of cross beam spacing on load distribution factor at girder bridges

Full Scale Experimental Evaluation for Cable Dampers

NUMERICAL ANALYSIS OF LOAD DISTRIBUTION IN RAILWAY TRACK UNDER WHEELSET

IJESRT. Scientific Journal Impact Factor: (ISRA), Impact Factor: METHODOLOGY Design Parameter [250]

DESIGN AND FABRICATION OF CHASSIS FOR ELECTRICAL VEHICLE

Vehicle Dynamic Simulation Using A Non-Linear Finite Element Simulation Program (LS-DYNA)

Modal analysis of Truck Chassis Frame IJSER

Development of analytical process to reduce side load in strut-type suspension

MARINE FOUR-STROKE DIESEL ENGINE CRANKSHAFT MAIN BEARING OIL FILM LUBRICATION CHARACTERISTIC ANALYSIS

Abaqus Technology Brief. Automobile Roof Crush Analysis with Abaqus

Finite Element Modeling and Analysis of Vehicle Space Frame with Experimental Validation

Structural Analysis of Student Formula Race Car Chassis

Forced vibration frequency response for a permanent magnetic planetary gear

Using ABAQUS in tire development process

Static and Dynamic Strength Analysis on Rear Axle of Small Payload Off-highway Dump Trucks

Reducing the Structural Mass of a Real- World Double Girder Overhead Crane

Optimization of Scissor-jack-Damper s Parameters and Performance under the Constrain of Human Comfort

THE LONGITUDINAL VIBRATION OF COMPOSITE DRIVE SHAFT

Design and Analysis of Engine and Chassis Mounting Bracket

VIBRATION REDUCTION IN CONVENTIONAL VEHICLES BY INCREASING THE STIFFNESS ON THE CHASSIS FRAME

Estimation of Unmeasured DOF s on a Scaled Model of a Blade Structure

Dynamic Response Assessment and Design Optimization of Aircraft Tyre Pressure Monitoring Unit (TPMU) Akshay B G 1 Dr. B M Nandeeshaiah 2

CHAPTER 5 PARAMETRIC STUDIES AND SQUEAL REDUCTION METHODS

STIFFNESS CHARACTERISTICS OF MAIN BEARINGS FOUNDATION OF MARINE ENGINE

STRESS AND VIBRATION ANALYSIS OF A GAS TURBINE BLADE WITH A COTTAGE-ROOF FRICTION DAMPER USING FINITE ELEMENT METHOD

Fatigue Life Estimation of Chassis Frame FESM Bracket for Commercial Vehicle

Heat treatment Elimination in Forged steel Crankshaft of Two-stage. compressor.

Design And Analysis Of Two Wheeler Front Wheel Under Critical Load Conditions

EXPERIMENTAL AND FINITE ELEMENT ANALYSIS OF SKEW GIRDER BRIDGES RAVINDRA SWAROOP. Department of Applied Mechanics. Submitted

Stress Analysis of Engine Camshaft and Choosing Best Manufacturing Material

Development of Rattle Noise Analysis Technology for Column Type Electric Power Steering Systems

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 1, No 4, 2011

RESILIENT INFRASTRUCTURE June 1 4, 2016

NUMERICAL ANALYSIS OF IMPACT BETWEEN SHUNTING LOCOMOTIVE AND SELECTED ROAD VEHICLE

Design and experiment of hydraulic impact loading system for mine cable bolt

Modeling and Simulation of Linear Two - DOF Vehicle Handling Stability

Analysis Of Gearbox Casing Using FEA

Design and Stress Analysis of Crankshaft for Single Cylinder 4-Stroke Diesel Engine

Experimental investigation on vibration characteristics and frequency domain of heavy haul locomotives

Chapter 7: Thermal Study of Transmission Gearbox

Composite Long Shaft Coupling Design for Cooling Towers

Structural Analysis of Pick-Up Truck Chassis using Fem

EXPERIMENTAL STUDY ON EFFECTIVENESS OF SHEAR STRENGTHENING OF RC BEAMS WITH CFRP SHEETS

Parametric Modeling and Finite Element Analysis of the Brake Drum Based on ANSYS APDL

Exploit of Shipping Auxiliary Swing Test Platform Jia WANG 1, a, Dao-hua LU 1 and Song-lian XIE 1

Prediction of wheel/rail rolling contact wear under the situation of wheel/rail vibration

Chapter 2 Dynamic Analysis of a Heavy Vehicle Using Lumped Parameter Model

Modal Analysis of Automobile Brake Drum Based on ANSYS Workbench Dan Yang1, 2,Zhen Yu1, 2, Leilei Zhang1, a * and Wentao Cheng2

Multi-axial fatigue life assessment of high speed car body based on PDMR method

CFD Investigation of Influence of Tube Bundle Cross-Section over Pressure Drop and Heat Transfer Rate

Analysis Methods for Skewed Structures. Analysis Types: Line girder model Crossframe Effects Ignored

STRUCTURAL ANALYSIS OF STEERING YOKE OF AN AUTOMOBILE FOR WITHSTANDING TORSION/ SHEAR LOADS

Numerical and Experimental Research on Vibration Mechanism of Rotary Compressor

1874. Effect predictions of star pinion geometry phase adjustments on dynamic load sharing behaviors of differential face gear trains

DESIGN AND ANALYSIS OF COMPOSITE LEAF SPRING

Semi-Active Suspension for an Automobile

Comparing FEM Transfer Matrix Simulated Compressor Plenum Pressure Pulsations to Measured Pressure Pulsations and to CFD Results

Transient Dynamic Analysis and Optimization of a Piston in an Automobile Engine

The Application of Simulink for Vibration Simulation of Suspension Dual-mass System

Experimental Study on Torsional Vibration of Transmission System Under Engine Excitation Xin YANG*, Tie-shan ZHANG and Nan-lin LEI

The effectiveness of CFRP strengthening of steel plate girders with web opening subjected to shear

Transmission Error in Screw Compressor Rotors

Review of Overhead Crane and Analysis of Components Depending on Span

Design Analysis and Optimization of Disc Brake

DESIGN OF A MODIFIED LEAF SPRING WITH AN INTEGRATED DAMPING SYSTEM FOR ADDED COMFORT AND LONGER LIFE

Automotive NVH with Abaqus. Abaqus 2018

Frontal Crash Simulation of Vehicles Against Lighting Columns in Kuwait Using FEM

Analysis on natural characteristics of four-stage main transmission system in three-engine helicopter

Design and Analysis of Front Lower Control Arm by Using Topology Optimization

Design and Analysis of Pressure Die Casting Die for Side Differential Cover of Mini truck

Simulating Rotary Draw Bending and Tube Hydroforming

Seismic Capacity Test of Overhead Crane under Horizontal and Vertical Excitation - Element Model Test Results on Nonlinear Response Behavior-

Assessing the needs for intermediate diaphragms in prestressed concrete girder bridges

Pantograph and catenary system with double pantographs for high-speed trains at 350 km/h or higher

Tooth Shape Optimization of the NGW31 Planetary Gear Based on Romax Designer

CFIRE December 2009

THERMAL STRESS ANALYSIS OF HEAVY TRUCK BRAKE DISC ROTOR

Static Analysis of Crankcase and Crankshaft of Single Cylinder Four Stroke Diesel Engine

Theoretical and Experimental Investigation of Compression Loads in Twin Screw Compressor

Torsional Stiffness Improvement of Truck Chassis Using Finite Elemen Method

Optimization Design of the Structure of the Manual Swing-out Luggage Compartment Door of Passenger Cars

2018 LOUISIANA TRANSPORTATION CONFERENCE. Mohsen Shahawy, PHD, PE

LAMINATED WINDSHIELD BREAKAGE MODELLING IN THE CONTEXT OF HEADFORM IMPACT HOMOLOGATION TESTS

Grand Challenge VHG Test Article 2 Test 4

A Certain Type of Wheeled Self-propelled Gun Independent Suspension Stress Analysis. Liu Xinyuna, Ma Jishengb

COST ANALYSIS FOR BRIDGE AND CULVERT. Essam A. Mostafa

Static And Free Vibration Analysis Of A Car Bonnet

Assessment of Fatigue and Modal Analysis of Camshaft

2. Write the expression for estimation of the natural frequency of free torsional vibration of a shaft. (N/D 15)

International Conference on Mechanics and Civil Engineering (ICMCE 2014)

BIKE SPIRAL SPRING (STEEL) ANALYSIS

Static Structural and Thermal Analysis of Aluminum Alloy Piston For Design Optimization Using FEA Kashyap Vyas 1 Milan Pandya 2

FINITE ELEMENT ANALYSIS OF TIE-ROD FOR SPACECRAFTS

MODELING SUSPENSION DAMPER MODULES USING LS-DYNA

Virtual Durability Simulation for Chassis of Commercial vehicle

Influence of Cylinder Bore Volume on Pressure Pulsations in a Hermetic Reciprocating Compressor

Transcription:

University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 8-2010 Dynamic characteristics of an FRP deck bridge Jing Song jsong12@utk.edu Recommended Citation Song, Jing, "Dynamic characteristics of an FRP deck bridge. " Master's Thesis, University of Tennessee, 2010. http://trace.tennessee.edu/utk_gradthes/750 This Thesis is brought to you for free and open access by the Graduate School at Trace: Tennessee Research and Creative Exchange. It has been accepted for inclusion in Masters Theses by an authorized administrator of Trace: Tennessee Research and Creative Exchange. For more information, please contact trace@utk.edu.

To the Graduate Council: I am submitting herewith a thesis written by Jing Song entitled "Dynamic characteristics of an FRP deck bridge." I have examined the final electronic copy of this thesis for form and content and recommend that it be accepted in partial fulfillment of the requirements for the degree of Master of Science, with a major in Civil Engineering. We have read this thesis and recommend its acceptance: Thomas L. Attard, Qiuhong Zhao (Original signatures are on file with official student records.) Zhongguo John Ma, Major Professor Accepted for the Council: Dixie L. Thompson Vice Provost and Dean of the Graduate School

To the Graduate Council: I am submitting herewith a thesis written by Jing Song entitled Dynamic characteristics of an FRP deck bridge. I have examined the final electronic copy of this thesis for form and content and recommend that it be accepted in partial fulfillment of the requirements for the degree of Master of Science, with a major in Civil Engineering. Zhongguo John Ma, Major Professor We have read this thesis and recommend its acceptance: Thomas L. Attard Qiuhong Zhao Accepted for the Council: Carolyn R. Hodges Vice Provost and Dean of the Graduate School (Original signatures are on file with official student records.)

Dynamic characteristics of an FRP deck bridge A Thesis Presented for the Master of Science Degree The University of Tennessee, Knoxville Jing Song August 2010

Acknowledgements I would like to thank my family who give me endless love and encouragement to pursue the master degree in University of Tennessee. I would like to express my deep gratitude to Dr. Zhongguo (John) Ma, who offers me this precious opportunity to do the research on the dynamic characteristics of an FRP deck bridge. I highly appreciate the guidance of Dr. Ma provided me on both of the research and life during these two years. I would like to thank Dr. Thomas L. Attard and Dr. Qiuhong Zhao who taught me a lot of professional knowledge in the class and serve as my committee members. Also I would like to thank Dr. Attard for the valuable guidance on my research. I am very grateful to the entire group in Tennessee Bridge Research Laboratory, TBRL, especially Xin Jiang and Wenchao Song for their assistance on the research and Abaqus learning. Also thanks to Lungui Li, Peng Zhu and Qi Cao for their encouragement and friendship during the two years. ii

Abstract Fiber reinforced polymer (FRP) deck has some significant advantages compared to concrete deck in use of bridges, such as light self-weight, high stiffness and strength, good durability and easy to install. FRP deck has already been used in some bridge rehabilitation and short span bridges. But for widely used in bridges, FRP deck bridges still need further research. Currently many research efforts focus on the filed tests of FRP deck bridges. Compared to field tests, Finite element analysis also has great advantages, such as low cost and convenient to conduct. Therefore, in this thesis finite element analysis is conducted by ABAQUS on the Boyer Bridge in Pennsylvania. The finite element model is verified by the static field test result. Then a simplified moving truck load is applied on the bridge model in order to analyze the dynamic responses of the FRP deck bridge, including the displacements and stress of each girder at the middle span. The dynamic effect is shown by comparing the dynamic responses and the static responses of the bridge. The connection between the FRP deck and girder is very important to the behavior of the bridge. In this thesis shear studs serve to connect the FRP deck and girder. This thesis also analyzes the effect of shear studs to the dynamic responses of the bridge by changing the number of the shear studs. Key words: Dynamic, FRP deck, Bridge. iii

Table of contents Introduction...1 Finite element model...4 Dynamic analysis...11 Natural frequencies...11 Truck load...13 Damping ratio...16 Dynamic responses...17 Conclusions...36 List of References...39 iv

List of Tables Table 1 Material properties...5 Table 2 Test Vehicle Axle Loads...5 Table 3 The strain obtained from test and FE...8 Table 4 FE results comparison...8 Table 5 Damping coefficients...17 Table 6 The displacements of Girder 1 at middle span...32 Table 7 The displacements of Girder 2 at middle span...32 Table 8 The displacements of Girder 3 at middle span...32 Table 9 The displacements of Girder 4 at middle span...32 Table 10 The displacements of Girder 5 at middle span...32 Table 11 The stress of Girder 1 at middle span...32 Table 12 The stress of Girder 2 at middle span...32 Table 13 The stress of Girder 3 at middle span...33 Table 14 The stress of Girder 4 at middle span...33 Table 15 The stress of Girder 5 at middle span...33 v

List of Figures Figure 1 Cut-away view of Boyer Bridge...4 Figure 2 Model of the bridge...7 Figure 3 Mesh of the model...7 Figure 4 Strains at the top of each girder at middle span...9 Figure 5 Strains at the bottom of each girder at middle span...9 Figure 6 First modal shape (Frequency =1.5764 cycles/time)...12 Figure 7 Second modal shape (Frequency=10.018 cycles/time)...12 Figure 8 Third modal shape (Frequency=10.806cycles/time)...13 Figure 9 Moving truck load positions...16 Figure 10 The displacement of girder 1 at middle span in case 1...18 Figure 11 The displacement of girder 2 at middle span in case 1...18 Figure 12 The displacement of girder 3 at middle span in case 1...19 Figure 13 The displacement of girder 4 at middle span in case 1...19 Figure 14 The displacement of girder 5 at middle span in case 1...19 Figure 15 The stress of girder 1 at middle span in case 1...20 Figure 16 The stress of girder 2 at middle span in case 1...20 Figure 17 The stress of girder 3 at middle span in case 1...20 Figure 18 The stress of girder 4 at middle span in case 1...21 Figure 19 The stress of girder 5 at middle span in case 1...21 Figure 20 The displacement of girder 1 at middle span in case 2...21 Figure 21 The displacement of girder 2 at middle span in case 2...22 Figure 22 The displacement of girder 3 at middle span in case 2...22 Figure 23 The displacement of girder 4 at middle span in case 2...22 Figure 24 The displacement of girder 5 at middle span in case 2...23 Figure 25 The stress of girder 1 at middle span in case 2...23 Figure 26 The stress of girder 2 at middle span in case 2...23 Figure 27 The stress of girder 3 at middle span in case 2...24 Figure 28 The stress of girder 4 at middle span in case 2...24 Figure 29 The stress of girder 5 at middle span in case 2...24 Figure 30 The displacement of girder 1 at middle span in case 3...25 Figure 31 The displacement of girder 2 at middle span in case 3...25 Figure 32 The displacement of girder 3 at middle span in case 3...25 Figure 33 The displacement of girder 4 at middle span in case 3...26 Figure 34 The displacement of girder 5 at middle span in case 3...26 Figure 35 The stress of girder 1 at middle span in case 3...26 Figure 36 The stress of girder 2 at middle span in case 3...27 Figure 37 The stress of girder 3 at middle span in case 3...27 Figure 38 The stress of girder 4 at middle span in case 3...27 Figure 39 The stress of girder 5 at middle span in case 3...28 vi

Figure 40 The displacement of girder 1 at middle span in case 4...28 Figure 41 The displacement of girder 2 at middle span in case 4...28 Figure 42 The displacement of girder 3 at middle span in case 4...29 Figure 43 The displacement of girder 4 at middle span in case 4...29 Figure 44 The displacement of girder 5 at middle span in case 4...29 Figure 45 The stress of girder 1 at middle span in case 4...30 Figure 46 The stress of girder 2 at middle span in case 4...30 Figure 47 The stress of girder 3 at middle span in case 4...30 Figure 48 The stress of girder 4 at middle span in case 4...31 Figure 49 The stress of girder 5 at middle span in case 4...31 Figure 50 The extreme displacements of Girder 1...33 Figure 51 The extreme displacements of Girder 2...34 Figure 52 The extreme displacements of Girder 4...34 Figure 53 The extreme displacements of Girder 4...35 Figure 54 The extreme displacements of Girder 5...35 vii

Introduction Bridges with FRP decks are being explored as a potential system with an accelerated constructible feature. The characteristics of bridges with FRP decks, such as mass, stiffness, and damping are significantly different from those of bridges with traditional concrete decks. Therefore the dynamic response of the FRP-deck bridges is of a great interest, and is the objective of the current research reported here. Some researchers have already worked on this area. Zhang and Cai (2006) studied the load distribution and dynamic response of multi-girder bridges with FRP decks and concrete decks based on a bridge-vehicle coupled model. They found that the load distribution factor values and dynamic response of FRP deck bridges are larger than those of concrete deck bridges. And also they found that FRP deck bridges with partially composite conditions have a larger girder load distribution and a larger dynamic displacement than those of the FRP deck bridges with fully composite conditions. Also they concluded that road roughness and vehicle velocity significantly affected the dynamic performance. Chiewanichakorn et al (2006) studied the behavior of a truss bridge where the old deteriorated concrete deck was replaced with a FRP deck. Using experimentally validated finite element models to conduct dynamic time-history analysis with an AASHTO fatigue truck over the bridge. They found that the fatigue life of the bridge after rehabilitation would be doubled compared to the pre-rehabilitation reinforced 1

concrete deck system. In this paper the damping ratio used is 5%, Rayleigh damping. Aluri et al (2005) studied the dynamic response of three FRP composite bridges through field tests. They concluded that the lowest damping ratios of FRP bridges in the test is 5%, which is lower than those of concrete bridges, which have average value of 7.9-8.4%; the acceleration of FRP bridges is high, beyond the serviceability limit of bridges; the dynamic load allowance factors are mostly within 1998 AASHTO LRFD bridge specification limits. Alampalli (2005) tested an FRP slab bridge. He observed the impact factor was about 0.3. The strains and deflections are lower than those predicted at the design stages. Higher modal damping values were observed compared to those for typical steel structures, which reflects the vibration absorbing capacity of FRP. Daly and Cuninghame (2005) tested a full-scale glass FRP bridge deck under static and dynamic wheel loading in the lab. The loads were imposed using the TRL Trafficking Test Facility. They found that FPR deck can resist local wheel loads due to heavy vehicles for at least 30-40 years without structural damage. And careful attention is needed to prevent local damage in highly stressed regions of the supporting deck, such as web to flange connections and close to bearing supports. Burgueno et al (2001) used vibration-based modal investigations as a health-monitoring and level I non-destructive evaluation. He found that modal vibration study was effective for determining the dynamic structural properties of 2

composite bridge system. The modal vibration test data, such as fundamental frequencies and modal shapes, were successful in determining the changes in structural behavior due to changes in boundary conditions, as well as structural damage. In this thesis, a finite element model was built of an FRP deck steel stringer bridge system using ABQUS. The model is verified by a static field test result on the Boyer Bridge by Yupeng Luo and Christopher J. Earls (2003). Based on the validated model, we further analyze the dynamic characteristics of this bridge, including the frequencies and modal shapes. Then a moving truck load is simplified to add on the bridge in order to analyze the dynamic responses of the FRP deck bridge, including the displacements and stress of each girder at the middle span. In this bridge, the FRP deck and steel stringers are connected by shear studs. The number of the shear studs will affect the stiffness of the bridge, which may affect the fundamental frequency and dynamic responses of the bridge. Therefore this thesis also analyzes the effect of the shear studs to the dynamic responses of the bridge by changing the numbers of the shear studs. 3

Finite element model The Boyer Bridge is a short-span (12.649 m) simply supported composite structure located in a secondary road in PENNDOT Engineering District 10-0. The cut-away view of Boyer Bridge is shown in Figure 1. It consists of five galvanized stringers acting compositely with five FRP deck panels. The FRP deck system is composed of tubes perpendicular to the traffic. The FRP deck and steel girders are connected by shear studs, which is 610mm spacing between each row. Each row has two shear studs across the steel girder section. The section and material properties are shown in the Table 1. A tandem-axle dump truck loaded with sand was chosen as test vehicle. Wheel loads were shown in the table 2. The truck was located on the second girder. From the field test, strains at the middle span of each girder were obtained. Figure 1 Cut-away view of Boyer Bridge 4

Table 1 Material properties Steel stringers (mm) Flange thickness(tf) = 19.05 bf = 323.85 tw = 12.70 Spacing = 1752.60 FRP deck (mm) Haunch thickness (t haunch) = 12.70 FRP flange thickness(ttop, tbtm) = 16.76 Deck thickness (td) = 194.56 Modulus of Elasticity (MPa) Esteel = 200000.00 Egrout = 31841.70 Efrp = 17241.40 Table 2 Test Vehicle Axle Loads Axle 1(kg) Axle 2(kg) Axle 3(kg) Left side 3,409 4,273 4,136 Right side 4,273 4,702 4,750 Total 7,682 9,000 8,886 Based on the test FRP deck system, the Finite element model was built with ABQUS. FRP deck is composed of top and bottom facings and core. The facings were simplified as an isotropic solid plate, which was modeled by eight-node solid elements (C3D8). To simplify, the contribution of the core to the load resistant is neglected. That is to say the modular elasticity of the core in the model is 0. The haunch and steel stringers are also modeled by eight-node and four-node (C3D4) solid elements. Each girder has forty two shear studs with two shear studs per one row across the girder section. The shear studs with diameter 22mm and height 150 mm are spaced 610mm 5

between the two rows. The bottom of the girders was simply supported. To exactly model the conditions in the field test, the truck load position in finite element model is the same as that in the field test. The centerline of truck was located at the centerline of the second girder. The front wheels were just off the bridge. In the model each two tires of the truck were simplified as a rectangular area with length 250mm and width 510mm. In the model, shear studs and steel girders are modeled as an entire body. Then the shear studs are embedded into the FRP deck to model the interaction between the shear stud and the FRP deck. The embedded element technique in ABAQUS is used to specify an element or a group of elements that lie embedded in a group of host elements whose response will be used to constrain the translational degrees of freedom of the embedded nodes. The interaction between the surfaces of the haunch and the girders are modeled with contact in ABQUS. The contact function can model the normal behavior and the tangential behavior between the surfaces. The normal behavior is set as hard contact, which means the pressure exists between the surfaces. The tangential behavior is set as frictionless, which means there is no friction between the surfaces. Therefore the shear resistance is provided only by the shear studs in this situation. The usage of contact in the finite element model makes the model more objective. But at the same time it significantly increases the calculation time. Therefore, the model with contact is analyzed to compare with the result of the model without contact. 6

Figure 2 shows the profile, boundary condition and the load condition of the bridge. In the model the total number of nodes is 27539 and total number of elements is 26451, including 11346 linear hexahedral elements of type C3D8R and 15105 linear tetrahedral elements of type C3D4 (on the haunch), as shown in Figure 3. From the calculation of the model, the strain at the middle span of the girders is obtained. Table 3 gives the comparison between the field test results and the finite element results. The differences between the two results can be accepted. So the finite element model is verified to do further parameter study. Figure 2 Model of the bridge Figure 3 Mesh of the model 7

Table 3 The strain obtained from test and FE Middle span strain 10^-6 Girder 1 Girder 2 Girder 3 Girder 4 Girder 5 top middle bottom test -52.3-64.5-50.5-13.1 -- FE -85-77.7-49.61-19.65-1.56 test 11.2 16.8 -- 9.3 -- FE 6.82 23.1 12.97 0.84 3.87 test 76.6 84.1 59.8 35.5 -- FE 98.64 123.89 75.55 21.32 9.3 Table 4 FE results comparison Middle span strain 10^-6 Girder 1 Girder 2 Girder 3 Girder 4 Girder 5 contact -85-77.7-49.61-19.65-1.56 top No contact -87.41-84.67-56.56-20.57-1.62 Fully composite -83.4-73.8-48 -20.7-3.1 contact 6.82 23.1 12.97 0.84 3.87 middle No contact 15.32 30.99 18.33 2.56 3.74 Fully composite 15.7 32.75 19.2 2.5 4.27 contact 98.64 123.89 75.55 21.32 9.3 bottom No contact 99.52 126.25 78.67 21.57 9.09 Fully composite 98.1 121.8 74.9 20.7 9.1 8

0-10 Girder No. 1 2 3 4 5 strain(*10^-6) -20-30 -40-50 -60-70 -80-90 -100 case 1: use "contact" between surfaces case 2: without "contact" between surfaces case 3: fully composite test Figure 4 Strains at the top of each girder at middle span strain(*10^-6) 140 120 100 80 60 40 case 1: use "contact" between surfaces case 2: without "contact" between surfaces case 3: fully composite test 20 0 1 2 3 4 5 Girder No. Figure 5 Strains at the bottom of each girder at middle span 9

In figure 4 and 5, Case 1 is that the deck and girders are connected by shear studs and also there is contact function between the adjacent surfaces of haunch and girder. Case 2 is that the deck and girders are connected only by shear studs without contact between the adjacent surfaces of haunch and girder. In other words, nodes of the two surfaces do not have any functions. Case 3 is that the deck and girders are tied together to model the fully composite behavior of the bridge. In ABAQUS, a tie constraint ties two separate surfaces together so that there is no relative motion between them. Figure 4 and 5 respectively show the strains at the top and bottom of each girder at middle span in the three cases and also from the test. It can be concluded that the strains obtained from modeling with contact between adjacent surfaces is 77.7/84.67=91.8% of that without contact between surfaces. That is to say modeling the bridge system with contact can get a more exact result. Also it is concluded that the shear studs in the test bridge provided 73.8/77.7=95% of fully composite action, which is good enough. There are some differences between the test result and the FE result, especially of the top strain at the edge girder, which might be caused by the railing on the edge of the bridge. The difference of the strain at the top of second girder is 17%. The difference of the strain at the bottom of second girder is 32%, which might be caused by the diaphragms. In the modeling diaphragms are not considered, which might make the distribution factor different and the moment in the second girder bigger than that in real. 10

Dynamic analysis The mass and stiffness between the FRP deck and concrete deck are quite different, which result the frequencies of the bridge system are different. Therefore this may generate dynamic problems when the bridge is subjected to live load, such as moving trucks. Also due to big differences of mass and stiffness between FRP deck and steel girder, the dynamic response of the bridge may be influenced by the high mode effect. Therefore, in this paper we analyze the frequencies and modes of this particular FRP deck bridge to find out whether the high mode affects the dynamic response. Furthermore, the dynamic response of the FRP deck bridge under a moving truck is analyzed. By comparing the static response and dynamic response, the dynamic influence of a moving truck to this bridge can be obtained. Natural frequencies Natural frequencies and modal shapes are basic dynamic characteristics of a system. The modal vibration test data, such as fundamental frequencies and modal shapes has been successfully used in bridge damage supervision. In ABAQUS, we can easily use the linear perturbation analysis step to get the natural frequencies of the bridge. In the frequency extraction model, the FRP deck is fully composite with steel girders. The first ten frequencies are 1.5764, 10.018, 10.806, 13.319, 13.835, 16.217, 16.740, 18.090, 19.811, and 20.933 (unit:hz), which are transferred to angular frequencies 9.9rad/s, 62.91 rad/s, 67.86 rad/s, 83.64 rad/s, 86.88 rad/s, 101.84 rad/s, 11

105.13 rad/s, 113.61 rad/s, 124.41 rad/s, and 131.44 rad/s. The first three mode shapes of the bridge are shown in Figures 6-8. Figure 6 First modal shape (Frequency =1.5764 cycles/time) Figure 7 Second modal shape (Frequency=10.018 cycles/time) 12

Figure 8 Third modal shape (Frequency=10.806cycles/time) Truck load Assume the truck travels at 20m/s across the bridge. The bridge length is 12.649m. The truck length is 5.68m. Therefore the total time for the truck to cross this bridge would be Travelling time = (12.649m+5.68m)/20m/s=0.92s (1) To simulate a moving truck, the truck is considered to locate at different positions at different time periods. In the finite element analysis, the total response time is set to be 3s in order to observe the response after the truck getting off the bridge. The total time is divided into 10 steps. The duration time of the first 9 steps is 0.1 second at each step. The last step is 2.1s duration time. The truck is loaded at different position at each step. Considering the truck starts to travel on the bridge from one end to another at 0s. Figure 9 shows the position of the truck at each step. At each step, the truck load is simplified as a step load lasting 0.1s. 13

Step 1 t=0~0.1s only the front wheel is on the bridge Step 2 t=0.1~0.2s the front wheel is 2m from the end Step 3 t=0.2~0.3s Step 4 t=0.3~0.4s 14

Step 5 t=0.4~0.5s Step 6 t=0.5~0.6s Step 7 t=0.6~0.7s Step 8 t=0.7~0.8s 15

Step 9 t=0.8~0.9s Step 10 t=0.9~3s Figure 9 Moving truck load positions Damping ratio In this finite element analysis, use Rayleigh damping: = a m a k (2) c 0 + 1 The damping ratio for the nth mode of such a system is: ζ n a0 1 a1 = + ωn (3) 2 ω 2 n The coefficients a 0 and a 1 can be determined from specified damping ratiosζ i and ζ j for the ith and jth modes, respectively. Expressing Eq. (2) for these two modes in matrix form leads to 1 2 1/ ωi 1/ ω j ωi a ω j a 0 1 ζ i = ζ j (4) If both modes are assumed to have the same damping ratioζ, which is reasonable 16

based on experimental data, then a 0 2ωiω j = ζ ω + ω i j 2 a1 = ζ (4) ω + ω i j Then the damping matrix is known from EQ. (2). The damping ratio ζ is chosen to be 5%. Assume the first and fifth modes have the same damping ratio, then the coefficients a 0 and a 1 can be obtained, as shown in Table 5. Table 5 Damping coefficients Damping ratio ζ ω 1 (rad/s) ω 5 (rad/s) a 0 a 1 5% 9.9 86.88 0.0889 1.033E-4 Dynamic responses The simplified moving truck load is added on the verified bridge model to obtain its dynamic responses. In order to investigate the influence of connection between FRP deck and girder to the dynamic response of the bridge, finite element analysis in four cases is conducted. 1) The FRP deck and girders are fully composite with shear studs inside, carried out by using tie constraint between the FRP deck and girders in ABAQUS. 2) Similar with case 1, but without shear studs inside. 3) The FRP deck and girders are partially composite by shear studs, carried out by using contact interaction between the FRP deck and girders in ABAQUS. 4) The FRP deck and girders are partially composite with the number of shear studs reduced by half. Under static load, the displacement at middle span of the second girder is 7.2mm in fully-composite condition of the bridge. Figure 10-14 show displacements at middle 17

span of each girder in case 1. Figure 15-19 show stress at bottom of each girder at middle span in case 1. Figure 20-24 show displacements at middle span of each girder in case 2. Figure 25-29 show stress at the bottom of each girder at the middle span in case 2. Figure 30-34 show displacements at middle span of each girder in case 3. Figure 35-39 show stress at bottom of each girder at middle span in case 3. Figure 40-44 show displacements at middle span of each girder in case 4. Figure 45-49 show stress at bottom of each girder at middle span in case 4. Figure 10 The displacement of girder 1 at middle span in case 1 Figure 11 The displacement of girder 2 at middle span in case 1 18

Figure 12 The displacement of girder 3 at middle span in case 1 Figure 13 The displacement of girder 4 at middle span in case 1 Figure 14 The displacement of girder 5 at middle span in case 1 19

Figure 15 The stress of girder 1 at middle span in case 1 Figure 16 The stress of girder 2 at middle span in case 1 Figure 17 The stress of girder 3 at middle span in case 1 20

Figure 18 The stress of girder 4 at middle span in case 1 Figure 19 The stress of girder 5 at middle span in case 1 Figure 20 The displacement of girder 1 at middle span in case 2 21

Figure 21 The displacement of girder 2 at middle span in case 2 Figure 22 The displacement of girder 3 at middle span in case 2 Figure 23 The displacement of girder 4 at middle span in case 2 22

Figure 24 The displacement of girder 5 at middle span in case 2 Figure 25 The stress of girder 1 at middle span in case 2 Figure 26 The stress of girder 2 at middle span in case 2 23

Figure 27 The stress of girder 3 at middle span in case 2 Figure 28 The stress of girder 4 at middle span in case 2 Figure 29 The stress of girder 5 at middle span in case 2 24

Figure 30 The displacement of girder 1 at middle span in case 3 Figure 31 The displacement of girder 2 at middle span in case 3 Figure 32 The displacement of girder 3 at middle span in case 3 25

Figure 33 The displacement of girder 4 at middle span in case 3 Figure 34 The displacement of girder 5 at middle span in case 3 Figure 35 The stress of girder 1 at middle span in case 3 26

Figure 36 The stress of girder 2 at middle span in case 3 Figure 37 The stress of girder 3 at middle span in case 3 Figure 38 The stress of girder 4 at middle span in case 3 27

Figure 39 The stress of girder 5 at middle span in case 3 Figure 40 The displacement of girder 1 at middle span in case 4 Figure 41 The displacement of girder 2 at middle span in case 4 28

Figure 42 The displacement of girder 3 at middle span in case 4 Figure 43 The displacement of girder 4 at middle span in case 4 Figure 44 The displacement of girder 5 at middle span in case 4 29

Figure 45 The stress of girder 1 at middle span in case 4 Figure 46 The stress of girder 2 at middle span in case 4 Figure 47 The stress of girder 3 at middle span in case 4 30

Figure 48 The stress of girder 4 at middle span in case 4 Figure 49 The stress of girder 5 at middle span in case 4 From the figures above, the maximum and minimum displacements of each girder are drawn in the 4 cases, shown in Table 6~10. The maximum and minimum stresses of each girder are drawn in the four cases, shown in Table 11~15. In the Table 6~10, - means the displacement is downward. In Table 11~15, - means stress is in compression. In order to reflect the changes of the maximum and minimum displacements of each girder in different cases, the column graphics are drawn, as shown in Figure 50~54. 31

Table 6 The displacements of Girder 1 at middle span Displacement Case 1 Case 2 Case 3 Case 4 Max (mm) 7 2.5 13.33 11.26 Min (mm) -13.57-15.95-10.56-11.18 Table 7 The displacements of Girder 2 at middle span Displacement Case 1 Case 2 Case 3 Case 4 Max (mm) 8.38 0 14.5 11.31 Min (mm) -14.05-16.58-11.07-11.42 Table 8 The displacements of Girder 3 at middle span Displacement Case 1 Case 2 Case 3 Case 4 Max (mm) 7.99 0 12.25 9.81 Min (mm) -9.55-11.26-7.15-7.63 Table 9 The displacements of Girder 4 at middle span Displacement Case 1 Case 2 Case 3 Case 4 Max (mm) 9.36 2.39 13.36 10.23 Min (mm) -7.69-5.58-6.45-6.97 Table 10 The displacements of Girder 5 at middle span Displacement Case 1 Case 2 Case 3 Case 4 Max (mm) 10.54 0 11.41 11.21 Min (mm) -7.80-11.26-7.75-8.09 Table 11 The stress of Girder 1 at middle span Stress Case 1 Case 2 Case 3 Case 4 Max (MPa) 62.89 64.43 63.88 62.92 Min (MPa) -34.66-9.38-39.01-32.52 Table 12 The stress of Girder 2 at middle span Stress Case 1 Case 2 Case 3 Case 4 Max (MPa) 72.89 73.56 71.23 65.92 Min (MPa) -41.51-4.34-48.05-30.53 32

Table 13 The stress of Girder 3 at middle span Stress Case 1 Case 2 Case 3 Case 4 Max (MPa) 46.87 47.28 47.7 46.15 Min (MPa) -38.07-1.54-41.24-31.69 Table 14 The stress of Girder 4 at middle span Stress Case 1 Case 2 Case 3 Case 4 Max (MPa) 36.71 22.47 41.07 40.16 Min (MPa) -45.01-11.03-48.47-36.29 Table 15 The stress of Girder 5 at middle span Stress Case 1 Case 2 Case 3 Case 4 Max (MPa) 35.98 15.19 46.55 43.76 Min (MPa) -48.39-12.42-38.68-41.34 15 10 Upward Displ. Downward Displ. Displacement (mm) 5 0-5 -10 1 2 3 4-15 -20 Case No. Figure 50 The extreme displacements of Girder 1 33

20 15 Upward Displ. Downward Displ. Displacement(mm) 10 5 0-5 -10 1 2 3 4-15 -20 Case No. Figure 51 The extreme displacements of Girder 2 15 10 Upward Displ. Downward Displ. Displacement (mm) 5 0-5 1 2 3 4-10 -15 Case No. Figure 52 The extreme displacements of Girder 4 34

15 Upward Displ. Downward Displ. 10 Displacement (mm) 5 0 1 2 3 4-5 -10 Case No. Figure 53 The extreme displacements of Girder 4 15 10 Upward Displ. Downward Displ. Displacement (mm) 5 0-5 1 2 3 4-10 -15 Case No. Figure 54 The extreme displacements of Girder 5 35

Conclusions A three-dimensional finite element model of a fiber-reinforced polymer deck bridge is developed based on the use of commercial software ABAQUS. A simplified static truck load is added on the bridge to obtain its response. Comparison of finite element analysis results against field test results indicates that the model can be used to perform extensive parametric study. Based on the verified model, a moving truck load is added to analyze dynamic responses of the bridge. This study investigated its dynamic responses in four conditions: deck and girder fully composite with/without shear studs inside, deck and girder partially composite, deck and girder partially composite with shear studs reduced by half. From finite element analysis the following conclusions can be drawn: 1) In static analysis, finite element results show that the bridge has got 95% composite action. However, the finite element results of dynamic response in fully composite condition and partially composite condition have significant differences. This indicates that the connection status between FRP deck and girder has much more influence on the bridge responses under dynamic load than that under static load. This also shows finite element model used to conduct dynamic analysis requires more attention to make sure it correct. 2) In dynamic response, the moving truck not only causes downward displacements but also large upward displacements in partially composite condition. 36

The downward displacements of each girder are larger in fully composite condition than those in partially composite condition. However, difference between the upward and downward displacements of Girder 2 is 54% larger in partially composite condition than those in fully composite condition, which cause more discomfort to passengers. In free vibration phrase after the truck load removed, the displacements of bridge in fully composite condition damps out quickly. However, the bridge in partially composite condition still has large vibration in that period. 3) When the number of shear studs is reduced by half, the maximum upward displacement of Girder 2 is 22% (1-11.31/14.5=22%) smaller than that with the full number of shear studs. The frequency in the free vibration phrase in the dynamic response is reduced by half. This is caused by that the bridge becomes more flexible when the shear studs are reduced. Therefore the natural frequency of the system is reduced. More vibration energy is dissipated during the vibration. This indicates that if the bridge is partially composite, reduction of shear studs is helpful for the dynamic response. 4) In fully composite condition with shear studs inside the girder (case 1), the upward displacement is larger than that in fully composite condition without shear studs inside (case 2). But the downward displacement in case 1 is smaller than that in case 2. The value of displacement in case 1 is between that in case 2 and case 3. Therefore, in the Finite element analysis it is very important how to simplify the modeling part. For example, whether to include the shear studs or not in the finite 37

element modeling will cause obvious differences to the results. To verify the finite element modeling, the results should compare with the field test results. Based on the analysis results, when the FRP deck is used in bridges, it is recommended that connection between FRP deck and girders should be strengthen to have a fully composite condition in order to minimize the dynamic responses. In addition, in the dynamic analysis with finite element model, the calculation takes about four days to finish three seconds response of the bridge. So it is very important to simplify the model in order to reduce the calculation cost. If the time is enough, we can extend the response time of the bridge until the dynamic response damps out. 38

List of References 39

1. Yin Zhang, C.S. Cai. Load distribution and dynamic response of multi-girder bridges with FRP decks. Engineering Structures; 29(2007):1676-1689. 2. M. Chiewanichakorn et al. Dynamic and fatigue response of a truss bridge with fiber reinforced polymer deck. International Journal of Fatigue; 29(2007): 1475-1489. 3. Srinivas Aluri, and et al. Dynamic response of three fiber reinforced polymer composite bridges. Journal of Bridge Engineering ASCE 2005; 10(6):722-730. 4. Screenivas Alampalli. Field performance of an FRP slab bridge. Composite Structures 2006; 72(4): 494-502. 5. Albert F. Daly, John R. Cuninghame. Performance of a fibre-reinforced polymer bridge deck under dynamic wheel loading. Composites: Part A; 37(2006): 1180-1188. 6. Rigoberto Burgueno, et al. Experimental dynamic characterization of an FRP composite bridge superstructure assembly. Composite Structures; 54(2001): 427-444. 7. American Association of State Highway and Transportation officials (AASHTO). (2004). LRFD bridge design specifications, Washington, D.C. 8. D.C. Keelor, A.M.ASCE; Y. Luo; C.J. Earls, M.ASCE. Service load effective compression flange width in fiber reinforced polymer deck systems acting compositely with steel stringers. Journal of Composites for Construction ASCE 2004; 8(4): 289-297. 40

9. ABAQUS/CAE user s manual. 10. Anil K. Chopra, Dynamics of structures (Second Edition). 41

Vita Jing Song was born in Henan Province, China on January 23, 1984. She was raised by her parents, Xiaohu Song and Suqing Huang in Wuzhi, Henan Province. She attended Xi an Jiaotong University in China in August 2001 and received Bachelor Degree in Civil Engineering in July 2005. In August 2005, she attended Tongji University in China, receiving her master degree in Civil Engineering with a concentration in Urban Mass Transit in 2008. Currently Jing Song is pursuing a Master of Science Degree in Civil Engineering with a concentration in structures at The University of Tennessee, Knoxville. 42