Document Type: Regular Article
Author
Faculty of Civil Engineering, Semnan University, Semnan, Iran
Abstract
Highlights
Keywords
Main Subjects
In the reinforcement concrete structures, shear failure of the element is very destructive and it is highly regarded in the design of these type of elements. Shear failure of RC columns mainly due to weakness in transverse reinforcement is a common failure in the past and in the nonductile RC joints which have a low percent of the transverse reinforcement, there is a high risk of shear failure. There are a lot of RC structures that because of lack of knowledge of the vulnerability were built nonductile, especially in their joints. Many studies have been done to strengthen and improve the performance of these elements and several methods such as using FRP material were also proposed. The studies showed that shear failure of RC joint elements is very dangerous and very special attention is required. Some researchers investigated the shear strength of RC joint with different concrete types. For example, McLean and Pierce [1] investigated the shear of RCC (roller compacted concrete) joints based on an experimental study. They have presented the safety factors which can be used in analyses. Another case is a study which was done by Shiohara [2] to the analysis of the high strength reinforcement concrete joint in shear failure. The analysis of RC joint is a very useful tool to study the behavior of these elements. Ghobarah and Biddah [3] proposed a joint element for modeling of the joint in the nonlinear dynamic analysis with considering shear deformation. Their results show that the modeling of inelastic shear deformation in joints has a significant effect on the seismic response. Bakir and Boduroglu [4] presented a design equation for determining the shear strength of monotonically loaded exterior RC joints. They used several parametric studies to investigate the influence of variables on the behavior of RC joints based on the experimental database. Their results showed that their equation can be able to predict the joint shear strength exterior RC joints. An analytical model for shear strength of high strength RC joints is done by Sayed [5]. He was presented a general model for these type of joints.
One of the studies about the shear capacity of the RC joints is done by Jaehong and LaFave [6]. They used a collection of an extensive database of reinforced concrete (RC) beam–column connection test specimens which were subjected to cyclic lateral loading. They have determined the influence parameters for joint shear stress and finally, the design checks recommended were examined. They also presented probabilistic joint shear strength models for design [7].
The joint shear strength of exterior concrete beamcolumn joints reinforced internally with Glass Fibre Reinforced Polymer (GFRP) reinforcements was investigated by Saravanan and Kumaran [8]. They tested eighteen specimens and used finite element analysis to simulate the behavior of the beamcolumn joints. A design equation for assessing the joint shear strength of the GFRP reinforced beamcolumn specimens was also proposed. Sharma et al. [9] presented a model for simulating the shear behavior of exterior reinforced concrete joints subjected to seismic loads. Their model does not need any special element or subroutine and uses limiting principal tensile stress in the joint. Shear behavior of ultrahigh performance concrete was studied by Lee et al. [10]. The results of their tests have been compared with several design formulae for assessing the joint shear strength. The available models to prediction the shear strength of beam–column joints were reviewed by Pradeesh et al. [11]. The concept, parameters considered, significant observations and their limitations of the models for predicting the joint shear behavior were summarized in their study. Elshafiey et al. [12] investigated the performance of exterior RC joints subjected to a combination of shear and torsion based on the results of an experimental study. They also presented a threedimensional truss model and showed that their model had an agreement with the experimental results. The shear strength and behavior of beamcolumn joints in unbonded precast prestressed concrete (PCaPC) frames based on the test results were investigated by Jin et al. [13]. The joint shear input was compared with the nominal shear strength of RC joint panels which was calculated based on common standards in their study.
This paper is an attempt to determine the shear capacity of RC nonductile joints based on artificial neural networks (ANN), adaptive neurofuzzy inference system (ANFIS) and also group method of data handling (GMDH). A collection of experimental which were published in the literature were used and predictive models were proposed to estimation the shear capacity.
Soft computing (SC) tried to build intelligent which provides the ability to derive the answer from the problems with high dimensions and complex. They used to develop systems in similar of the human mind and have been advantageous in many engineering applications. In a general classification, SC techniques can be classified into three groups including artificial neural networks, fuzzy systems and neurofuzzy system (which is the combination of the first two groups).
In this section, three most famous soft computing methods including ANN, ANFIS and also GMDH were reviewed. These approaches are the considered methods which were used to the aim of this paper.
Artificial neural networks (ANN) are systems which were widely used for function approximation based on a collection of existing samples. they can be able to train the solutions from these data. They are applied in areas where the presentation of an answer is difficult by traditional methods. They have been used to solve engineering problems by three general layers namely input layer, hidden layer or layers and also output layer. There are several neurons as computational units. These neurons are connected in layers, and signals travel from the first (input) to the last (output) layer. ANN used a set of data to estimate the weights and bias for the input and output signals of each neuron. It is clear that a big and reliable dataset has more ability to estimate the parameters.
Adaptive neurofuzzy inference system (ANFIS) is a fuzzy inference system implemented in the framework of adaptive networks which was introduced by Jang [14]. It was applied both of the fuzzy rules and input–output data pairs. ANFIS is one of the powerful soft computing approaches which was presented a Sugenotype fuzzy system in a fivelayer network (the input layer not counted by Jang). They are the ability of ANN and fuzzy systems together. For create an ANFIS model, three methods are commonly used: grid partition (GP), subtractive clustering (SC) and also fuzzy cmeans (FCM) clustering. GP algorithm divides the data space into rectangular subspaces. SC algorithm divides the considered data into groups called clusters to discover the solution patterns. FCM which was used in this paper is an unsupervised algorithm. FCM consider the dataset into fuzzy clusters and also allows one data to belong to two or more clusters. This can be very useful to have a flexible and strongest ANFIS.
Group Method of Data Handling (GMDH) network which introduced by Ivakhnenko [15] is a multilayered perceptrontype network structure for mathematical modeling of systems. It is able to get the solution algorithm using data samples. Each node in GNDH has two input signals and use a secondorder polynomial based on these two inputs. A collection of the dataset is applied to determine the coefficient values of polynomials based on least squares approach. They also can selfneglect ineffective inputs. Because of the mathematical manner of GMDH, these type of networks is widely used in engineering problems.
For calculation of the considered soft computing methods, a collection of 149 data which were published in literature was used [1650]. These data are related to nonductile RC joints which were a failure in shear and their shear capacity has been reported to them. Table 1 and Fig.1 provides the details of the considered dataset. In this table, and also are beam height, beam index, beam longitudinal reinforcement, yield stress of beam longitudinal reinforcement, a ratio of the number of subassemblages, The effective width of the joint panel and also shear strength of the joint respectively.
Table 1: Information of dataset
JP 

Average 
343.38 
0.26 
0.015 
493.81 
30.95 
0.81 
226.52 
5.62 
Maximum 
762.00 
0.78 
0.041 
746.00 
100.80 
1.00 
600.00 
10.45 
Minimum 
150.00 
0.06 
0.003 
315.00 
8.30 
0.75 
100.00 
1.19 
Median 
300.00 
0.24 
0.013 
500.00 
31.60 
0.75 
200.00 
5.45 
Mode 
300.00 
0.15 
0.009 
500.00 
31.60 
0.75 
200.00 
5.08 
St.Dev 
139.69 
0.14 
0.007 
91.37 
10.83 
0.11 
100.16 
1.94 
Range 
612.00 
0.73 
0.038 
431.00 
92.50 
0.25 
500.00 
9.26 
JP (In plane geometry) =1 for interior, 0.75 for exterior.
Fig 1. Distribution of the considered dataset
To normalization, a relationship which created the data within the value of 0.1 to 0.9 is used by Eq.1:
(1)
To training the models, 126 data, which has randomly chosen from the dataset, was used. The remained 22 data means used for the testing phase of the proposed models.
The initial modeling of the considered SC methods showed that the most powerful inputs were , BI and JP. Therefore, these parameters were selected and used as Inputs. This section, the structures and the parameters of the proposed models for considering estimation were presented in details. The results were discussed in section 5.
The proposed ANN structure was shown in Fig.2. The shear strength was considered by in the figure. , and are also the normal values of input 1, 2 and 3 respectively. They considered as X1, X2 and X3 in this paper. It was clear from the figure that the hidden layer has eight neurons. These nodes transfer its values to the final layer by TangentSigmoid function. For the output layer, Purelin function was used. The details of the layers were presented in Tables 2 and 3. In these tables, and are the bias of the hidden and output layer respectively.
Figure 2. The proposed ANN structure
Table 2. Layer weights and bias for the final layer
Layer weights 

Neuron 1 
Neuron 2 
Neuron 3 
Neuron 4 
Neuron 5 
Neuron 6 
Neuron 7 
Neuron 8 

0.3565 
0.3483 
0.9888 
0.3566 
0.3564 
0.3584 
1.4283 
1.0587 
0.6522 
Table 3. Input weights and bias for the hidden layer
Neuron 
Input weights 

Input 1 
Input 4 
Input 5 

Neuron 1 
0.1299 
0.0532 
0.2158 
0.1388 
Neuron 2 
0.1282 
0.0577 
0.2109 
0.1358 
Neuron 3 
1.3420 
0.2517 
0.5692 
0.2158 
Neuron 4 
0.1299 
0.0532 
0.2158 
0.1388 
Neuron 5 
0.1299 
0.0533 
0.2157 
0.1388 
Neuron 6 
0.1302 
0.0515 
0.2172 
0.1395 
Neuron 7 
0.9260 
0.0178 
1.0580 
1.1287 
Neuron 8 
0.6463 
0.3571 
1.0240 
0.6312 
The selected ANFIS model, used FCM algorithm and had Gaussian membership function (eq.2) for input parameters as follows:
(2)
Where c is the mean and σ is the variance of x. The proposed ANFIS structure presented in Fig.3.
Figure 3. The proposed ANFIS structure
The Gaussian parameters of the membership functions presented in Table 4 for all input parameters. Fig.46 showed membership functions of the selected ANFIS.
Table 4: Gaussian membership function’s parameters
Membership function 
Parameter 
Inputs 

X1 
X2 
X3 

C1 
c 
0.2883 
0.1013 
0.3440 
σ 
0.0461 
0.0672 
0.0311 

C2 
c 
0.5382 
0.8988 
0.1734 
σ 
0.0878 
0.1779 
0.0359 

C3 
c 
0.4737 
0.3059 

σ 
0.0619 
0.0236 

C4 
c 
0.1954 
0.2234 

σ 
0.0514 
0.0407 

C5 
c 
0.3641 
0.2394 

σ 
0.0566 
0.0315 

C6 
c 
0.3116 
0.4174 

σ 
0.0327 
0.0485 

C7 
c 
0.2048 
0.3273 

σ 
0.0567 
0.0286 

C8 
c 
0.4959 
0.2384 

σ 
0.0745 
0.0216 

C9 
c 
0.2528 
0.2905 

σ 
0.0479 
0.0242 
Figure 4. Membership functions for input 1
Figure 5. Membership functions for input 2
Figure 6. Membership functions for input 3
In ANFISFCM structure, there are several clusters for the target. Each of clusters includes a linear function which is showed in eq.3.
j=1,..,10 
(3) 
The parameters are coefficients of input . The parameter is deal with a constant value. The amounts of these parameters presented in Table 5.
For the selected ANFIS, the rule base and also rule’s weights showed in Table 6 and 7 respectively.
Table 5: Parameters of the output’s clusters.
Cluster 

Inputs coefficients ( 

Constant (C) 
CL1 
3.8580 
0.2218 
1.7900 
2.2180 
CL2 
0.5828 
0.0854 
0.2812 
0.8538 
CL3 
0.5490 
5.8310 
5.3360 
4.6020 
CL4 
0.9090 
0.0661 
3.8080 
0.6607 
CL5 
0.4401 
0.0207 
0.5480 
0.2067 
CL6 
1.4010 
0.0793 
0.8190 
0.7933 
CL7 
4.3990 
0.1364 
1.6800 
1.3640 
CL8 
1.4100 
0.2046 
5.5820 
2.0460 
CL9 
0.4603 
4.6700 
1.6390 
4.1180 
CL10 
1.5740 
0.0695 
2.4130 
0.5313 
Table 6: ANFIS rules
Number 

Rules 
Rule 1 

If XI is and X2 is and X3 is then is CL1. 
Rule 2 

If XI is and X2 is and X3 is then is CL2. 
Rule 3 

If XI is and X2 is and X3 is then is CL3. 
Rule 4 

If XI is and X2 is and X3 is then is CL4. 
Rule 5 

If XI is and X2 is and X3 is then is CL5. 
Rule 6 

If XI is and X2 is and X3 is then is CL6. 
Rule 7 

If XI is and X2 is and X3 is then is CL7. 
Rule 8 

If XI is and X2 is and X3 is then is CL8. 
Rule 9 

If XI is and X2 is and X3 is then is CL9. 
Rule 10 

If XI is and X2 is and X3 is then is CL10. 
The normal value of the joint shear strength based on the considered ANFISFCM model can be determined by eq.4.
(4) 
Table 7: Rule’s weight
Number 

Weight’s relationship 



















The GMDH structure which was used in this paper presented in Fig.7. The predictive model has two polynomials in the middle layer with equations 5 and 6.
Figure 7. Membership functions for input 3
(5) 

(6) 
Based on the previous polynomials (eq.5 and 6), the final output of the model was calculated by eq.7.
(7) 



The output values of the proposed models are normal value and need to be converted to its real value and for this purpose, eq.8 was used:

(8) 
In the equation, is the joint shear strength of the nonductile RC joints which determine by the proposed models. Based on the real values, the distribution of the results of these models presented in Fig810. It was clear from the figures that the considered soft computing approaches had suitable predictions.
Figure 8. Distributed results for all 126 train data
Figure 9.Distributed results for all 22 test data
Figure 10. Distributed results for all 149 data
A summary of the final results was also presented in Table.8. It was concluded that although all of the considered methods had suitable results, for all 149 data, ANFIS had less error and higher correlation factor than other models.
Table 7: Summary results
Model 
Train data (126 data) 


Test data (22 data) 

All data (149 data) 

MAE 
RMSE 


MAE 
RMSE 

MAE 
RMSE 

ANN 
0.905 
0.565 
0.807 


0.928 
0.765 
0.924 

0.910 
0.658 
0.807 
ANFIS 
0.939 
0.526 
0.646 


0.904 
0.793 
0.938 

0.932 
0.567 
0.699 
GMDH 
0.875 
0.727 
0.911 


0.929 
0.633 
0.821 

0.886 
0.713 
0.897 
In this table: R2 is correlation coefficient, MAE is mean absolute error and RMSE is the root mean squared error.
Determination of the shear strength of nonductile of RC joint using three soft computing methods including ANN, ANFIS, and GMDH was considered in this paper. For train and test the models, a collection of experimental was used and the structures of the predictive models presented in details. it was mention that based on try and error approach for all of the considered methods, three inputs including and JP had more effective on the shear strength and therefore was used for modeling. The results showed that the proposed models have high performance for determining the shear strength. Additionally, the predicted values by ANFIS was more accurate than other two models. The importance results and the predictive models which were presented in this paper can be very useful for purposes such as retrofitting.
[1] F. G. McLean and J. S. Pierce, "Comparison of joint shear strengths for conventional and roller compacted concrete," in Roller Compacted Concrete II:, 1988, pp. 151169.
[2] S. Hitoshi, "Analysis of Joint Shear Failure Of HighStrength Reinforced Concrete Interior BeamToColumn Joint," Journal High Strength Concrete. ASCE, pp. 114, 1998.
[3] A. Ghobarah and A. Biddah, "Dynamic analysis of reinforced concrete frames including joint shear deformation," Engineering Structures, vol. 21, pp. 971987, 1999.
[4] P. Bakir and H. Boduroğlu, "A new design equation for predicting the joint shear strength of monotonically loaded exterior beamcolumn joints," Engineering Structures, vol. 24, pp. 11051117, 2002.
[5] A. A. Sayed, "General Analytical Model for Nominal Shear Stress of Type 2 Normaland HighStrength Concrete BeamColumn Joint," ACI Structural Journal, vol. 101, pp. 6575, 2004.
[6] J. Kim and J. M. LaFave, "Key influence parameters for the joint shear behaviour of reinforced concrete (RC) beam–column connections," Engineering Structures, vol. 29, pp. 25232539, 2007.
[7] J. Kim and J. M. LaFave, "Probabilistic joint shear strength models for design of RC beamcolumn connections," ACI Structural Journal, vol. 105, p. 770, 2008.
[8] J. Saravanan and G. Kumaran, "Joint shear strength of FRP reinforced concrete beamcolumn joints," Central European Journal of Engineering, vol. 1, pp. 89102, 2011.
[9] A. Sharma, R. Eligehausen, and G. Reddy, "A new model to simulate joint shear behavior of poorly detailed beam–column connections in RC structures under seismic loads, Part I: Exterior joints," Engineering Structures, vol. 33, pp. 10341051, 2011.
[10] C. Lee, Y. Kim, W. Chin, and E. Choi, "Shear strength of ultrahigh performance fiber reinforced concrete (UHPFRC) precast bridge joint," in High Performance Fiber Reinforced Cement Composites 6, ed: Springer, 2012, pp. 413420.
[11] L. V. Pradeesh, S. Sasmal, K. Devi, and K. Ramanjaneyulu, "Evaluation of Models for Joint Shear Strength of Beam–Column Subassemblages for Seismic Resistance," in Advances in Structural Engineering, ed: Springer, 2015, pp. 885896.
[12] T. M. Elshafiey, A. M. Atta, H. M. Afefy, and M. E. Ellithy, "Structural performance of reinforced concrete exterior beam–column joint subjected to combined shear and torsion," Advances in Structural Engineering, vol. 19, pp. 327340, 2016.
[13] K. Jin, K. Kitayama, S. Song, and K.o. Kanemoto, "Shear Capacity of Precast Prestressed Concrete BeamColumn Joint Assembled by Unbonded Tendon," ACI Structural Journal, vol. 114, p. 51, 2017.
[14] J.S. Jang, "ANFIS: adaptivenetworkbased fuzzy inference system," IEEE transactions on systems, man, and cybernetics, vol. 23, pp. 665685, 1993.
[15] A. Ivakhnenko, "Polynomial theory of complex systems," IEEE Transactions on Systems, Man, and Cybernetics, pp. 364378, 1971.
[16] U. Akguzel, "Seismic performance of FRP retrofitted exterior RC beamcolumn joints under varying axial and bidirectional loading," 2011.
[17] Y. A. AlSalloum, N. A. Siddiqui, H. M. Elsanadedy, A. A. Abadel, and M. A. Aqel, "Textilereinforced mortar versus FRP as strengthening material for seismically deficient RC beamcolumn joints," Journal of Composites for Construction, vol. 15, pp. 920933, 2011.
[18] C. P. Antonopoulos and T. C. Triantafillou, "Experimental investigation of FRPstrengthened RC beamcolumn joints," Journal of composites for construction, vol. 7, pp. 3949, 2003.
[19] T.H. Chen, "Retrofit strategy of nonseismically designed frame systems based on a metallic haunch system," 2006.
[20] C. Clyde, C. P. Pantelides, and L. D. Reaveley, Performancebased evaluation of exterior reinforced concrete building joints for seismic excitation: Pacific Earthquake Engineering Research Center, College of Engineering, University of California, Berkeley, 2000.
[21] R. De Otiz, "Strutandtie modelling of reinforced concrete: short beams and beamcolumn joints," University of Westminster, 1993.
[22] R. P. Dhakal, T.C. Pan, P. Irawan, K.C. Tsai, K.C. Lin, and C.H. Chen, "Experimental study on the dynamic response of gravitydesigned reinforced concrete connections," Engineering Structures, vol. 27, pp. 7587, 2005.
[23] M. Engindeniz, Repair and strengthening of pre1970 reinforced concrete corner beamcolumn joints using CFRP composites: Georgia Institute of Technology, 2008.
[24] Y. Goto and O. Joh, "An experimental study of shear failure mechanism of RC interior beamcolumn joints," in 11 th World Conference on Earthquake Engineering, 1996.
[25] S. J. Hamil, "Reinforced concrete beamcolumn connection behaviour," Ph.D. Thesis, , Durham University, UK, 2000.
[26] W. M. Hassan, Analytical and experimental assessment of seismic vulnerability of beamcolumn joints without transverse reinforcement in concrete buildings: University of California, Berkeley 2011.
[27] A. Ilki, I. Bedirhanoglu, and N. Kumbasar, "Behavior of FRPretrofitted joints built with plain bars and lowstrength concrete," Journal of Composites for Construction, vol. 15, pp. 312326, 2010.
[28] C. Karayannis, C. Chalioris, and K. Sideris, "Effectiveness of RC beamcolumn connection repair using epoxy resin injections," Journal of Earthquake Engineering, vol. 2, pp. 217240, 1998.
[29] W.T. Lee, Y.J. Chiou, and M. Shih, "Reinforced concrete beam–column joint strengthened with carbon fiber reinforced polymer," Composite Structures, vol. 92, pp. 4860, 2010.
[30] B. Li, Y. Wu, and T.C. Pan, "Seismic behavior of nonseismically detailed interior beamwide column joints—Part II: Theoretical comparisons and analytical studies," ACI Structural Journal, vol. 99, pp. 791802, 2003.
[31] B. Oh, K. Park, H. Hwang, and H. Chung, "An experimental study on shear capacity of reinforced concrete exterior beamcolumn joint with high strength concrete," Proceedings of the Architecture Institute of Korea, vol. 12, pp. 363366, 1992.
[32] Y. Ohwada, "A study on RC beamcolumn connection subjected to lateral load (8)," in (in Japanese), Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan, 1970, pp. 737738.
[33] Y. Ohwada, "A study on RC beamcolumn connection subjected to lateral load (9)," in (in Japanese), Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan, 1973, pp. 12971298.
[34] Y. Ohwada, "A study on effect of lateral beams on RC beamcolumn joints (1) " (in Japanese), Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan, pp. 14551456, 1976.
[35] Y. Ohwada, "A study on effect of lateral beams on RC beamcolumn joints (2) " (in Japanese), Proceedings of the Architectural Institute of Japan, vol. 61, pp. 241244, 1977.
[36] Y. Ohwada, "A study on effect of lateral beams on RC beamcolumn joints (4)," Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan, pp. 15111512, 1980.
[37] C. P. Pantelides, J. Hansen, J. Nadauld, and L. D. Reaveley, "Assessment of reinforced concrete building exterior joints with substandard details," PEER report, vol. 18, 2002.
[38] S. Park and K. M. Mosalam, "Experimental investigation of nonductile RC corner beamcolumn joints with floor slabs," Journal of Structural Engineering, vol. 139, pp. 114, 2012.
[39] D. E. Parker and P. Bullman, "Shear strength within reinforced concrete beamcolumn joints," Structural Engineer, vol. 75, 1997.
[40] A. Pimanmas and P. Chaimahawan, "Shear strength of beam–column joint with enlarged joint area," Engineering structures, vol. 32, pp. 25292545, 2010.
[41] I. B. Salim, "The influence of concrete strengths on the behaviour of external beamcolumn joints ", MS Thesis, Universiti Teknologi Malaysia, 2007.
[42] R. Scott, "The effects of detailing on RC beam/column connection behaviour," Structural Engineer, vol. 70, 1992.
[43] T. Supaviriyakit and A. Pimanmas, "Comparative performance of substandard interior reinforced concrete beam–column connection with various joint reinforcing details," Materials and Structures, vol. 41, pp. 543557, 2008.
[44] H. Taylor, "The behavior of in situ concrete beamcolumn joints," 1974.
[45] A. G. Tsonos, "Cyclic load behavior of reinforced concrete beamcolumn subassemblages of modern structures," ACI Structural journal, vol. 104, p. 468, 2007.
[46] A. G. Tsonos, "Effectiveness of CFRPjackets and RCjackets in postearthquake and preearthquake retrofitting of beam–column subassemblages," Engineering Structures, vol. 30, pp. 777793, 2008.
[47] A. G. Tsonos and K. V. Papanikolaou, "Postearthquake repair and strengthening of reinforced concrete beamcolumn connections (theoretical & experimental investigation)," BulletinNew Zealand society for earthquake engineering, vol. 36, pp. 7393, 2003.
[48] Y.C. Wang and K. Hsu, "Shear strength of RC jacketed interior beamcolumn joints without horizontal shear reinforcement," ACI Structural Journal, vol. 106, p. 222, 2009.
[49] Y.C. Wang and M.G. Lee, "Rehabilitation of nonductile beamcolumn joint using concrete jacketing," in A paper presentation at 13th world conference on earthquake engineering, Vancouver (BC, Canada), 2004.
[50] H. F. Wong, "Shear strength and seismic performance of nonseismically designed reinforced concrete beamcolumn joints," Ph.D. Thesis, Department of Civil Engineering, The Hong Kong University of Science and Technology, 2005.