Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure

Similar documents
On-Line Monitoring of Cable Trench and Cable Pit by Using Metal Armored FBG Temperature Sensors

Simulation of Micro Blanking Process of Square Hole with Fillet Based on DEFORM-3D

Long-Term Mechanical Properties of Smart Cable Based on FBG Desensitized Encapsulation Sensors

All-Optical Flip-Flop Based on Coupled SOA-PSW

BioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.

Experimental Study on Dual-Wavelength Distributed Feedback Fiber Laser

DISTRIBUTION STATEMENT A 7001Ö

Design for Testability

2. Problem formulation

ECE Real Time Embedded Systems Final Project. Speeding Detecting System

Structure-vibration Analysis of a Power Transformer (154kV/60MVA/Single Phase)

650MHz/800kW Klystron Development at IHEP

SPACE FLIGHT CABLE MODEL DEVELOPMENT. Author: Kaitlin Spak Advisor: Dr. Daniel Inman Virginia Polytechnic Institute and State University

Contrastive Analysis and Research on Negative Pressure Beam Tube System and Positive Pressure Beam Tube System for Mine Use

Investigation of Two Bidirectional C + L Band Fiber Amplifiers with Pumping Sharing and Wavelength Reused Mechanisms

Skip Length and Inter-Starvation Distance as a Combined Metric to Assess the Quality of Transmitted Video

Scan. This is a sample of the first 15 pages of the Scan chapter.

Auto classification and simulation of mask defects using SEM and CAD images

Monitor QA Management i model

Keywords: Edible fungus, music, production encouragement, synchronization

Annual Report of the IFLA-PAC China Center

CML 700i Measuring light curtain with integrated control unit and interfaces for even simpler operation

Distortion Analysis Of Tamil Language Characters Recognition

Optical Engine Reference Design for DLP3010 Digital Micromirror Device

9 Simulation of Pull Tests for Grouted Cable Anchors

ENGINEERING COMMITTEE Interface Practices Subcommittee AMERICAN NATIONAL STANDARD ANSI/SCTE

Drift Tubes as Muon Detectors for ILC

Huiping Guo, Fuzeng Yang* College of Mechanical and Electronic Engineering, Northwest A&F University, Yangling , China

High Value-Added IT Display - Technical Development and Actual Products

M4000 Diagnostic Test System For Power Apparatus Condition Assessment

THE DESIGN OF CSNS INSTRUMENT CONTROL

How to Manage Color in Telemedicine

ModBox-850nm-NRZ-series

8000 Plus Series Safety Light Curtain Installation Sheet ( CD206A/ CD206B )

Predicting the immediate future with Recurrent Neural Networks: Pre-training and Applications

843-R 843-R LASER POWER METER USER MANUAL. NEWPORT CORPORATION

Cable Retention Force Testing of Trunk & Distribution Connectors

OCC Installation Conduit Guidelines Excerpt from Optical Cable Corporation s INSTALLATION GUIDE

Avoiding False Pass or False Fail

Specification for Loose Tube Fiber Optic Cable (Non-Metallic, Dry Block, Figure-8) (G.652.D)

Coherent Receiver for L-band

Research on Control Strategy of Complex Systems through VSC-HVDC Grid Parallel Device

CARLITE grain orien TEd ELECTRICAL STEELS

Machine Vision System for Color Sorting Wood Edge-Glued Panel Parts

Detection and demodulation of non-cooperative burst signal Feng Yue 1, Wu Guangzhi 1, Tao Min 1

Monday 20 May 2013 Afternoon

Compact multichannel MEMS based spectrometer for FBG sensing

SCHEDA TECNICA. Amarro per cavi 13,5mm Amarro per cavi 16,5mm. MS217F00 Amarro per cavi 13,5mm MS218F00 Amarro per cavi 16,5mm

THE NEXT GENERATION OF CITY MANAGEMENT INNOVATE TODAY TO MEET THE NEEDS OF TOMORROW

Research on Precise Synchronization System for Triple Modular Redundancy (TMR) Computer

from ocean to cloud ADAPTING THE C&A PROCESS FOR COHERENT TECHNOLOGY

Real-time QC in HCHP seismic acquisition Ning Hongxiao, Wei Guowei and Wang Qiucheng, BGP, CNPC

A Demonstration Platform for Small Satellite Constellation Remote Operating and Imaging

Spec No.: ZTT Revision: 2 Date: Oct 3, Designer. Senior Technical Manager. Approver. Chief Technical Engineer

FPGA Implementation of DA Algritm for Fir Filter

MTS/T-BERD 8000 Platform Optical Spectrum Analyzer Modules

TRAN-COR H. grain ORIENTED ELECTRICAL STEELS. Applications Potential. a significant increase in core loss.

This is a repository copy of Virtualization Framework for Energy Efficient IoT Networks.

Exploratory Analysis of Operational Parameters of Controls

NONDESTRUCTIVE INSPECTION OF A COMPOSITE MATERIAL SAMPLE USING A LASER ULTRASONICS SYSTEM WITH A BEAM HOMOGENIZER

ModBox-1310nm-1550nm-NRZ 1310nm & 1550 nm, 28 Gb/s, 44 Gb/s Reference Transmitters

MG-XV operating instruction. Measuring of norm signals, 4-8-digit. Panel instrument type MG-BV Construction instrument type MG-AV

Multi Core fibers and other fibers for the future.

Flicker Caused by Operation of Industrial Technology

GLI-12 V1.1 GLI 12 V2.0

WITH the rapid development of Gallium Nitride

Quick Start Function Summary Instructions for ASHCROFT GC52 Differential Pressure Transmitter Version 6.03 Rev. B

2 Any form of reproduction if this document, also electronic reproduction, such as placing it in the internet is not allowed.

The Research of Controlling Loudness in the Timbre Subjective Perception Experiment of Sheng

Sharif University of Technology. SoC: Introduction

Guide for installers. METTLER TOLEDO MultiRange System solution analog Ex1. Hazardous area. Safe area

A Noble Identification System for Tumors and Implementation

Quality Control. Chapter 23

Failure Analysis Technology for Advanced Devices

A Defect Grading Method for Digital Radiographic Image

High Resolution Concrete Imaging

Informatics Enlightened Station 2 Park distance control system

Oxymax W COS61. Technical Information. Dissolved oxygen measurement Optical sensor acc. to the fluorescence quenching principle

Detecting and Analyzing System for the Vibration Comfort of Car Seats Based on LabVIEW

Optical Receiver Manual. Transmitter OP-OR112R JⅢ. Shenzhen Optostar Optoelectronics Co., Ltd (Version 2)

VLSI Chip Design Project TSEK06

MTS/T-BERD 8000 Platform

WELDING CONTROL UNIT: TE 450 USER MANUAL

ModBox-CBand-NRZ series C-Band, 28 Gb/s, 44 Gb/s, 50 Gb/s Reference Transmitters

v1.0.0 January AlphaLab, Inc. All rights reserved. TriField EMF Meter Owner s Manual

Intelligent Pendulum Hardness Tester BEVS 1306 User Manual

Emerging Subsea Networks

ModBox-1310nm-1550nm-28Gbaud-PAM nm & 1550 nm, 28 Gbaud PAM-4 Reference Transmitter

System Quality Indicators

The University of Texas at Austin September 30, 2011

SIMULATION MODELING FOR QUALITY AND PRODUCTIVITY IN STEEL CORD MANUFACTURING

Application of Measurement Instrumentation (1)

EDDYCHEK 5. Innovative eddy current testing for quality and process control. Touchscreen. Networking. All major applications.

Nondestructive Testing Device for Tie Bar Placement Accuracy

Recent APS Storage Ring Instrumentation Developments. Glenn Decker Advanced Photon Source Beam Diagnostics March 1, 2010

>Inductive. REFERENCE guide. Proximity

MICROSENS. Fast Ethernet Switch Modul 4x 10/100Base-TX, 1x 100Base-FX. Description. Features

Installation of Optical Fiber

Fiber Optic. Foreword. Special Tips:

Chord Classification of an Audio Signal using Artificial Neural Network

Transcription:

PHOTONIC SENSORS / Vol. 4, No. 4, 2014: 366 372 Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure Sheng LI 1*, Min ZHOU 2, and Yan YANG 3 1 National Engineering Laboratory for Fiber Optic Sensing Technology, Wuhan University of Technology, Wuhan, 430070, China 2 Wuhan WUTOS Limited Company, Wuhan, 430223, China 3 Key Laboratory of Fiber Optic Sensing Technology and Information Processing, Ministry of Education, Wuhan University of Technology, Wuhan, 430070, China * Corresponding author: Sheng LI E-mail: lisheng@whut.edu.cn Abstract: The smart cable with embedded distributed fiber optical Bragg grating (FBG) sensors was chosen as the object to study a new diagnosis method about broken wires of the bridge cable. The diagnosis strategy based on cable force and stress distribution state of steel wires was put forward. By establishing the bridge-cable and cable-steel wires model, the broken wires sample database was simulated numerically. A method of the characterization cable state pattern which can both represent the degree and location of broken wires inside a cable was put forward. The training and predicting results of the sample database by the back propagation (BP) neural network showed that the proposed broken wires diagnosis method was feasible and expanded the broken wires diagnosis research area by using the smart cable which was used to be only representing cable force. Keywords: Fiber optical Bragg grating sensor, bridge cable, broken wires diagnosis, smart cable, BP neural network Citation: Sheng LI, Min ZHOU, and Yan YANG, Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure, Photonic Sensors, 2014, 4(4): 366 372. 1. Introduction The cable is an important bearing component for long-span cable-stayed bridges. Rust and fracture are common diseases for cables, which may lead to the failure of cable and bridge accidents [1]. Therefore, monitoring the cable force state and diagnosing the potential disease of broken wires are of great significance to protect the safety of bridge operation. At present, nondestructive detection is the main means for broken wires diagnosis, including the methods based on magnetostriction effect, the eddy current effect, and magnetic effect. Other acoustic detection methods such as ultrasonic, solid acoustic and ray are still in the level of laboratory testing and expensive [2 5]. In addition, in theory measuring the internal stress distribution state of steel wires and cable force simultaneously can also provide the broken wires information. However, current cable force detection methods can only provide cable force information and are unable to measure the steel wires stress distribution state. The smart cable with embedded optical fiber Bragg grating (FBG) sensors is a structure [6] that can realize both bearing and sensing. Based on the large Received: 10 April 2014 / Revised version: 4 September 2014 The Author(s) 2014. This article is published with open access at Springerlink.com DOI: 10.1007/s13320-014-0196-x Article type: Regular

Sheng LI et al.: Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure 367 strain FBG sensors [7] that are fastened to the peripheral wires in the cable fabricating the stage and calibration relationship between the cable force and steel wires strain in the over tensioning stage, the smart cable can output both cable force and steel wires stress. Utilizing the above characters of the smart cable, the research will focus on a proposed broken wires diagnosis strategy, establishing a broken wires sample database based on the bridge-cable and cable-steel model by numerical simulating and finishing the diagnosis method verification by the back propagation (BP) neural network. 2. Broken wires diagnosis method 2.1 Diagnosis basis As shown in Fig. 1, several large strain FBG sensors are fastened to the peripheral wires in the connection pipe region during the cable fabricating process. The large strain FBG sensor with a special encapsulation structure [8, 9] meets the working range requirement for the stay cable, which can endure a large strain close to 5000 με under the over tensioning stage and the maximum combined loads in the operational stage. Besides the large range, as shown in Fig. 2, the mechanical fasten means for the sensor ensures a long-term effectiveness of force transmission. As shown in Fig. 3, taking a 430-meter-length smart cable working in JingYue Bridge for example, the strain sensors wavelength drift time history curves under the over tensioning stage is shown in Fig. 4. Six large strain FBG sensors and two temperature FBG sensors as temperature compensation are embedded inside the cable. From the curves, the stress distribution response of the monitored steel wires shows a little nonuniformity. However, all FBG strain sensors meet the over tensioning loading requirement. The calibration relationship between the cable force and steel wires stress delegated by FBG strain sensors wavelength response is established in the over tensioning process. Thus, the smart cable can be used as the research basis for broken wires diagnosis through the cable force and steel wires stress distribution. Output optical fiber Reserved steel tube A FBG sensor protective shield A Anchor slab Socket Connection pipe Steel wires Cable with PE protection Connection pipe Cable FBG strain sensor FBG temperature sensor Protective shield Steel wire A-A Fig. 1 Structure of the smart cable with embedded FBG sensors.

368 Photonic Sensors Fig. 2 Mechanical fasten between the sensor and cable steel wire. input information, namely, cable force and stress distributions of the monitored steel wires. The establishment for the output information should consider the number and layout of embedded large strain FBG sensors inside the smart cable. The broken wires forecasting model is established by training the BP neural network based on the above proposed input and output indexes. Load information Numerical model Broken wires pattern Cable force Steel wires stress Fig. 3 Smart cable used in real bridge engineering. 3000 2500 Smart cable monitoring Cable force Steel wires stress BP neural network training Broken wires forecasting model BP neural network pattern recognition Wavelength (pm) 2000 1500 1000 500 55.1%pb 1#sensor 2#sensor 3#sensor 4#sensor 5#sensor 6#sensor Diagnosis results of broken wires Fig. 5 Diagnosis strategy of broken wires inside a cable. 3. Numerical simulation 0 500 0 0.18 0.36 0.54 0.72 0.90 1.08 1.26 1.44 1.62 1.80 1.98 2.16 Time (s) 10 4 Fig. 4 FBG strain sensors wavelength drift of the smart cable under the over tensioning stage. 2.2 Diagnosis strategy The input and output information composes the pattern recognition sample database for the broken wires. The input information includes the cable force and steel wires stress distributions under various combination load cases. The output information is a vector that can reflect the degree and location of broken wires inside the cable. The complete broken wires sample database is derived from the numerical model which considers the typical environment and vehicle loading according to bridge vulnerability analysis. The diagnosis strategy process is shown in Fig. 5. The FBG-based smart cable provides the 3.1 Numerical model establishment A two-dimensional structure size cable-stayed bridge is used in [10] as the numerical simulation model basis. The structure material and geometric parameters are shown in Table 1. In order to simplify analysis, all stay cables in the model adopt the same specification (55φ7 parallel steel wires stay cable). The sizes of the bridge-cable model and numerical analysis model are shown in Fig. 6. Table 1 Structure material and geometric parameters. Item Beam Pylon Cable Pier Modulus of elasticity (MPa) 2.10E+05 3.20E+04 1.95E+05 3.40E+04 Poisson s ratio 0.2 0.167 0.2 0.167 Density (kg/m 3 ) 7850 2500 7850 2500 Moment of inertia I zz (m 4 ) 5.26 21.33 / 3.97 Moment of inertia I yy (m 4 ) 0.38 21.33 / 3.97 Section area (m 2 ) 0.86 16.36 0.00196 7.63

Sheng LI et al.: Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure 369 62 132 62 5 12 40 7 23 4 4 3 x 3 5 7 o y x 1 2 3 4 5 7 8 9 6 10 11 12 13 14 15 16 17 18 19 20 z 4 11.5 4 x 4 12 12 12 4 x 4 12 12 12 12 4 x 4 12 12 11.5 4 x 4 4 4 Fig. 6 Sizes of the bridge-cable model and numerical analysis model (unit: m). It is necessary for establishing the cable-steel wires local model to analyze the relationship between the cable force and stress distributions of the monitoring steel wires. If using helix to simulate 2 to 4 twisting angle of outermost layer steel wires, it will consume large number of nodes to form a lay length. Figure 7(a) shows a one-sixth lay length cable model which divides a lay length into 6 segments. With an increase in the mesh precision, each steel wire will use a number of truss elements. To improve the numerical simulation efficiency, assuming the cable steel wires completely parallel, the simplified model is used as shown in Fig. 7(b) instead of the real helix model. The assumption is convenient for selecting the steel wire element, saving computing nodes, and does not change the steel wires stress distribution at each cable section compared with the real model in the cable length direction. (a) real model (b) simplified model Fig. 7 Numerical analysis model of cable-steel wires. 3.2 Broken wires sample database design The Chinese code [11] points out that a new cable must be replaced if the number of broken wires in the cable is more than 2%. According to the proportion requirement, the numerical simulation for broken wires degree mainly considers 1 to 2 steel wires breaking for the 55φ7 parallel steel wires stay cable. Owing to verifying the effect of diagnosis strategy for broken wires, as shown in Fig. 8, the numerical simulation only considers the broken wires occurring in the peripheral steel wires under the dead load and assumes that the monitored steel wires attached to FBG strain sensors aren t broken and are only used to reflect the stress distribution at different broken wires states. Other complicated broken wires combination sample databases are omitted here. 1 6 7 2 8 9 10 11 12 3 13 14 15 16 17 18 4 19 20 21 22 23 24 25 5 26 27 28 29 30 31 35 49 50 51 52 53 54 55 34 43 44 45 46 47 48 33 38 39 40 41 42 32 36 37 Steel wires in 1# zone Steel wires in 3# zone FBG4 FBG2 Steel wires in 2# zone Steel wires monitored by sensors FBG3 Fig. 8 Number of the section steel wires and defined zone. Based on the above assumptions, according to the requirements for the damage identification steps, the aims of diagnosis strategy include correctly identifying whether broken wires occurs or not, the number of broken wires, namely, 1 or 2 steel wires

370 Photonic Sensors breaking, and corresponding regional location. Taking the 10# cable as the study object, by numerical calculation based on the bridge-cable model, Table 2 shows the cable force and steel wires average stress under the entire and broken state. Table 2 Force of 10# cable under structural dead load. Item Number of broken wires 0 1 2 Force (kn) 87599 86149 84696 Average stress (MPa) 41.39 41.46 41.52 The defined distribution of the broken wires region and serial number of steel wires are shown in Fig. 8. There are three FBG strain sensors used for sensing the cable force and stress distributions. The broken wires region is divided into 3 zones. When there is one steel wire breaking, the database considers all steel wires in the defined region except for the monitored steel wires, totally setting 11+11+11=33 samples. When there are two steel wires breaking, owing to the large number of sample combination, the database considers a special situation in which steel wires breaking occurs in the outermost layer and each region only has one steel wire breaking, totally setting 5 5+5 5+5 5=75 samples. Thus, the number of simulation samples is 109 including an entire state with no broken wires. For the input information, FBG 1 represents the average stress derived from the smart cable force, and FBG 2 FBG 4 represent the stress distributions derived from the steel wires monitored by FBG strain sensors. For the output information, seven vectors are defined, namely, entire state (0 0 0), one steel wire breaking in 1# zone (1 0 0), one steel wire breaking in 2# zone (0 1 0), one steel wire breaking in 3# zone (0 0 1), two steel wires breaking separately in 1# zone and 2# zone (1 1 0), two steel wires breaking separately in 2# zone and 3# zone (0 1 1), two steel wires breaking separately in 1# zone and 3# zone (1 0 1). Because only two types of broken wires degree are considered, the output vectors not only represent the broken wires region, but also distinguish the broken wires degree. Table 3 gives part of the sample database. Table 3 Sample database of broken wires inside the cable. Number Input Steel Output Sample FGB1 FBG2 FBG3 FBG4 wire Entire 1 41.39 41.41 41.37 41.37 0 0 0 1 2 41.46 47.39 37.56 42.42 1 0 0 25 13 41.46 42.17 43.42 38.93 0 1 0 32 24 41.46 38.97 42.12 43.41 0 0 1 1+25 35 41.52 44.88 41.62 38.40 1 1 0 2+32 60 41.52 40.99 39.69 44.04 1 0 1 25+32 109 41.52 39.74 44.04 40.94 0 1 1 3.3 Diagnosis effective evaluation The sample database training is finished by the BP neural network, which uses three-layer network structure. The stress distrubitions expressed by FBG 1 to FBG 4 that shown in Table 3 are used as the input layer. The hidden layer adopts the S-shaped tangent transfer function with 14 neurons decided by practical training effect, and the output layer adopts the S-shaped logarithmic transfer function. Three samples of one steel wire breaking and six samples of two steel wires breaking compose the test sample database. The other 100 samples including the entire state are used for model training. In order to improve the training convergence efficiency, based on (1), normalization for the input data is processed, which uses all input data in the numerical interval [0, 1]: xˆ x x x x max min min. (1) After 45 steps, network learning error satisfies the training target. According to (2), two types of noise level are considered to simulate the measurement error, where p i is the input data

Sheng LI et al.: Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure 371 polluted by noise; p a i is the original input data derived by numerical simulation; is the noise level expressed by percentage (1% and 5% noise levels are considered, respectively); R is an average distribution random variable in the numerical interval [ 1, 1]. a pi pi (1 R). (2) As shown in Table 4, represents that all recognition results are correct, and x represents that partial recognition results are feasible or except for the expecting recognition results there are misreport results. Before applying noise, recognition results are totally correct, which can predict both the degree and location of the broken wires. After applying noise, the recognition accuracy rate decreases which is due to noise pollution and the incompleteness of the training sample database for broken wires. The total number of the broken wires is 15 for the 9 test samples. Under the 1% noise level, the predicted 13 broken wires at the right region are recognized with no misreport and 86.7% accuracy rate for the 9 test samples. Under the 5% noise level, 12 broken wires for the 9 test samples are predicted, but among the recognition results there are 2 misreports occurred for 1 and 2 broken wires situation, respectively. If taking misreport as prediction failure, only 10 broken wires are recognized successfully with 66.7% accuracy rate for the 9 test samples. Thus, the overall accuracy rate for diagnosing broken wires is more than 60% under the 5% noise level interference. Table 4 Testing sample and result of broken wires inside the cable. Test sample number Output no noise 1% noise 5% noise 3 15 x 26 47 53 x 71 x x 63 x 97 103 x 4. Conclusions A broken wires diagnosis strategy that using cable and stress distributions as indexes and BP neural network as the pattern recognition method is proposed based on the characters of the FBG-based smart cable which can reflect both force and stress distributions of steel wires. The feasibility of the proposed diagnosis method which considers noise interruption is verified by the numerical simulation sample database derived from the bridge-cable and cable-steel wires numerical analysis model. The results show that proposed strategy can provide a new way for the bridge engineering broken wires diagnosis. Acknowledgment The research work reported in this paper was supported by the National Engineering Laboratory for Fiber Optic Sensing Technology, Wuhan University of Technology, China. Thanks for the support of the Fundamental Research Funds for the Central Universities (WUT: 2014-IV-090) and the National Natural Science Foundation of China (Major Program: 61290310). Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and source are credited. References [1] J. Zhang and W. Wang, Manual of bridge inspection and reinforcement. Beijing: China Communication Press, 2007. [2] J. Ye, Discussion on fracture of steel wires for cable, Steel Construction, 2003, 18(65): 35 36, 45. [3] X. Wu, J. Wang, and S. Yang, Development of a testing system for detecting the defect in diagnoal bridge cable, Mechanical Science and Technology, 2001, 20(6): 901 903. [4] X. Wu, Yi. Kang, and H. Chen, Nondestructive testing method and technique for cable-stayed bridge steel cables, Key Engineering Materials, 2004, 270 273(7): 1493 1499.

372 Photonic Sensors [5] M. Cooper, Radiographic examination of wire ropes suspension cables, Materials Evaluation, 1992, 50(12): 1383 1384, 1388. [6] S. Li, D. Jiang, and X. Zhao, Desensitized encapsulation FBG sensor for smart cable in bridge, in 22nd International Conference on Optical Fiber Sensors (OFS22), Beijing, China, October 15 19, 2012. [7] S. Li and D. Jiang, Structural large strain monitoring based on FBG sensor, in Symposium on Photonics and Optoelectronics (SOPO2009), Wuhan, China, August 14 16, 2009. [8] D. Jiang, S. Li, L. Liu, S. Liu, E. Zhang, and X. Zhao, Fiber optical Bragg grating strain sensor for smart cable, Chinese Patent ZL200920188144.9, September 1, 2010. [9] S. Li and M. Zhou, Long-term mechanical properties of smart cable based on FBG desensitized encapsulation sensors, Photonic Sensors, 2014, 4(3): 236 241. [10] Q. Zhang, Q. Li, and L. Tang, Study on structural damage location identification by damage indices for cable-stayed bridges, Engineering Mechanics, 2008, 25(5): 163 169. [11] Ministry of construction of the People s Republic of China, Technical code of maintenance for city bridge CJJ99 2003. Beijing: China Architecture & Building Press, 2003.