Research Article  Open Access
Teerachai Deesomsuk, Tospol Pinkaew, "Effectiveness of Vehicle Weight Estimation from Bridge WeighinMotion", Advances in Civil Engineering, vol. 2009, Article ID 312034, 13 pages, 2009. https://doi.org/10.1155/2009/312034
Effectiveness of Vehicle Weight Estimation from Bridge WeighinMotion
Abstract
The effectiveness of vehicle weight estimations from bridge weighinmotion system is studied. The measured bending moments of the instrumented bridge under a passage of vehicle are numerically simulated and are used as the input for the vehicle weight estimations. Two weight estimation methods assuming constant magnitudes and timevarying magnitudes of vehicle axle loads are investigated. The appropriate number of bridge elements and sampling frequency are considered. The effectiveness in term of the estimation accuracy is evaluated and compared under various parameters of vehiclebridge system. The effects of vehicle speed, vehicle configuration, vehicle weight and bridge surface roughness on the accuracy of the estimated vehicle weights are intensively investigated. Based on the obtained results, vehicle speed, surface roughness level and measurement error seem to have stronger effects on the weight estimation accuracy than other parameters. In general, both methods can provide quite accurate weight estimation of the vehicle. Comparing between them, although the weight estimation method assuming constant magnitudes of axle loads is faster, the method assuming timevarying magnitudes of axle loads can provide axle load histories and exhibits more accurate weight estimations of the vehicle for almost of the considered cases.
1. Introduction
The weights of vehicles govern the design requirements for highway infrastructure such as pavements and bridges. The traditional weigh stations are commonly used to weigh vehicles and impose fines or penalties for exceeding weight limits. They, however, take quite long time to weigh each vehicle. Moreover, their costs of system installation and maintenance are expensive. Therefore new alternative weigh system, namely, bridge weighinmotion (BWIM) is developed. The bridge WIM systems deal with an existing instrumented bridge or culvert from the road network. Beside its cost advantage, the system installation and maintenance do not disturbed the traffic flow. In addition, the system is transparent to the vehicle’s drivers so that the obtained weight information is expected to be unbiased. In general, BWIM system monitors the deflection, strain, or bending moment data of the bridge during the passages of vehicles. Knowing the physical parameters of the bridge such as span length and flexural rigidity, the system can estimate the weights of passing vehicles from those bridge response data coupled with the configuration and speed information of vehicles which are obtained from another set of sensors. Based on previous researches, although many techniques have been proposed for bridge WIM to estimate the weights of vehicles, two different assumptions of vehicle loads on the bridge, which is either constant or timevarying moving loads, are often employed.
For the constant moving loads assumption, the vehicle is assumed to pass the bridge without any vertical body motion. Therefore, dynamic moving loads from vehicle exerting on the bridge can be simply replaced by constant moving loads. The weight estimation methods, in this class of assumption, determine axle weights of the vehicle by comparing the measured bridge responses with those obtained from bridge influence lines [1–6]. The determined individual axle weights are then added together to estimate the gross weight of vehicle. Since the vehicle is assumed to be constant moving loads, the dynamic effects of the vehiclebridge system significantly affect the accuracy of the estimated vehicle weights. Consequently, some latter research works employ a digital filtering technique to eliminate dynamic effects and reconstruct quasistatic responses of the bridge from measured bridge responses [7–9]. Although using these quasistatic responses yields more accurate weight estimation, the difficulty in selection of the proper filtering parameters and the cutoff frequency makes the method impracticable.
Similar to many load identification methods, the vertical body motion of vehicle induced by vehiclebridge dynamic interaction is allowed. Therefore, unlike the constant loads assumption, the dynamic moving loads from vehicle exerting on the bridge are represented by timevarying loads moving on the bridge and are estimated directly from measured bridge responses. Then the axle weights of vehicle are determined from time averaging of the obtained timevarying axle loads. Many loads identification methods such as the time domain, the frequencytime domain, and the modal methods have been proposed and studied [10–13]. To obtain the estimated axle loads, the solution methods using either pseudoinverse or singular value decomposition technique are often adopted. However, it is found that all three methods exhibit large fluctuation of the identified loads due to measurement noises and possess numerically illcondition, especially when the axle load is on bridge supports. In addition, the methods consume long computing time due to the inversion of large system matrices. Therefore, the leastsquare method with smoothening term named the regularization method is employed [14, 15]. The discrete version of the method is also considered using the dynamic programming technique [16]. Besides the efficiency of computation, this method eliminates an illconditioned problem and provides better identified axle loads under noisy inputs. Unfortunately, the method needs an optimal regularization parameter in the identification process. To overcome this problem, the regularization method with the iterative technique called the updated static component (USC) technique is adopted [17, 18]. It has been numerically and experimentally shown that the regularization with USC technique can substantially improve the estimation accuracy over the conventional regularization. The studies were also extended to the continuous bridges [19–23]. The obtained results reveal the effectiveness of the loads identification methods of this class.
Most of the mentioned weight or load estimation methods have been studied to show their effectiveness and potential for real application. However, the comparison between them has not yet been established. Therefore, in this paper, the two weight estimation methods of the vehicles using the constant moving loads and timevarying moving loads assumptions are extensively considered. Based on numerical simulations, many effects of vehicle and bridge WIM system such as bridge discretization, sampling frequency, vehicle speed, bridge surface roughness, number of measuring sections, noise, axle spacing, and axle weight distribution are investigated.
2. Simulation of VehicleBridge Responses
2.1. Vehicle Model
A passage of vehicle on a bridge WIM system is shown in Figure 1. The system composes of a vehicle moving at a constant speed interacting with an instrumented bridge structure. The vehicle is modeled using 4 degrees of freedom consisting of vertical displacement rotation of vehicle body vertical displacement of front axle suspension mass and vertical displacement of rear axle suspension mass The equations of motion of vehicle are derived using dynamic equilibrium of the vehicle system in each degree of freedom and can be expressed by where
while and are bridge deflections beneath front and rear axles of the vehicle, respectively.
2.2. Bridge Model
The bridge structure is modeled as a simplysupported bridge and is discretized by finite elements using beam element. The standard beam element having 2 nodes with 2 degrees of freedom in vertical displacement and rotation displacement at each node as shown in Figure 2 is adopted. The equations of bridge motion can be represented by
where , and denote nodal bridge response vectors. and are mass, damping, and stiffness matrices of the bridge, respectively. is a matrix of global external load shape function used to transform the external front and rear axle loads of the vehicle to the nodal loads on the bridge model, while is the distances of two axle loads on the beam element in local coordinates as in Figure 2 and is the total number of the beam elements.
2.3. VehicleBridge Interaction Model
The vehiclebridge interaction model can be formulated from the equations of vehicle motion (1), and the equations of bridge motion (9). Combining these equations through the force equilibrium of vehicle axle loads, the governing equations of the vehiclebridge interaction system can be expressed by
where the matrices and variables contained in above equations are listed in the appendix.
It should be pointed out that the vehiclebridge system forms a coupled timevarying dynamic system, because some elements in damping and stiffness matrices of the system keep changing with time due to a traveling of the vehicle.
To solve above equations, Newmark’s method is employed. Based on the obtained system responses, the bending moment of the bridge at any bridge section as shown in Figure 2 can be determined from
where , , and are the obtained nodal displacements of the corresponding beam element and is the local location of the measuring section determined from the global location In addition, the axle loads of vehicle can be determined from (5), (7), and (8) as
It is noted that these axle load equations take into account both static axle weights and their dynamic loads resulted from the vehiclebridge interaction.
3. Weight Estimation of Vehicles
The concept of axle weight estimation of a passing vehicle is to minimize the error between measured and estimated bridge responses. In this paper, the measured bending moment vector at selected measuring points of the bridge under a moving vehicle is assumed to be simulated by solving the vehiclebridge dynamic interaction (10) and introducing the obtained nodal bridge responses into (11). The corresponding estimated bending moment vector, at the same bridge sections is approximated by solving equations of motion of the bridge subjected to a pair of moving axle loads. In general, these axle loads are assumed to be either constant or timevarying magnitudes as will be discussed in the following sections. Then the moving axle loads of the vehicle are identified through the minimization of the square error of the bending moments of the bridge:
where and are discrete forms of measured and estimated bending moment vectors of the bridge at time step respectively.
is the positivedefinite weighting matrix and is the total number of data points (in time).
3.1. Method I: Constant Magnitude of Moving Axle Loads Assumption
With constant magnitude of moving axle loads assumption, the dynamic interaction loads between vehicle and bridge, that is, and as previously defined in (7), are assumed to be have zero mean and can be omitted for static weight estimation. Therefore the vehicle is simply replaced by a pair of constant static weights of moving axle loads as in Figure 3. To identify these constant axle loads of vehicle, the estimated bending moments of the bridge are calculated from static influence lines of moving loads and are compared with their corresponding measured bending moments. This method is simple but its main disadvantage is that it cannot provide any dynamic information of the axle loads of vehicle. The influence line of the bending moment of a simplysupported bridge at a section as in Figure 3 can be expressed by a triangular function:
where is the static influence line of bridge bending moment at measuring section defined according to the load location,
is the distance of the moving front axle load, or rear axle load, from the left support of the bridge, and is the distance of measuring section from the left support of bridge.
Consequently, the estimated bending moment vector of the bridge at measuring sections induced by the movement of front and rear axle loads can be computed by the method of superposition as
Since the estimated bending moment of the bridge, , depends on only a set of estimated moving axle weights, . Recalling (13), the objective function in term of time discretization can be rewritten as
To minimize above objective function, a Sequential Quadratic Programming (SQP) method is used. This method iteratively solves a series of quadratic programming (QP) problems, using the Broyden, Fletcher, Goldfarb, and Shanno (BFGS) formula to update the Hessian of the Lagrangian, :
where are the constraints including upper and lower bounds of axle weights, and are the Lagrange multipliers.
The MATLAB’s optimization function fmincon [24], which utilizes previously mentioned procedures, is adopted to find the optimal axle weights, which yields minimum error between the measured and the estimated bending moment vectors. These obtained optimal axle weights are constant and are assumed to be the best estimated axle weights for the passing vehicle.
3.2. Method II: TimeVarying Magnitude of Moving Axle Loads Assumption
Unlike the constant magnitudes of moving axle loads assumption, the axle loads of vehicle are assumed to be timevarying. With this assumption, the dynamic interaction loads between vehicle and bridge, that is, and as well as the static axle weights of vehicle, that is, and are taken into account. Therefore the vehicle is replaced by a pair of timevarying magnitudes of moving axle loads similar to the system shown in Figure 2. To identify the axle loads of vehicle, the estimated bending moments of bridge are calculated solely from the bridge equilibrium (9). Then the axle loads identification is accomplished through minimizing the output error between the measured and the estimated bridge responses using the conventional regularization expressed in timediscretization form [16] as
where is the unknown moving loads vector with timevarying magnitude at time step and is the regularization parameter.
In this paper, the dynamic programming method with updated static component technique [17] is employed to minimize above error equation. This is because the method yields better solution and robustness against the choice of than the conventional regularization method. It is noted that the obtained axle loads, from (18) are timevarying loads. Therefore, to estimate the corresponding axle weights of vehicle, the simple timeaveraging is adopted:
in which and are the estimated front and rear axle weights, respectively.
4. Numerical Examples
The numerical investigation of axle weight estimation of a passing vehicle on a bridge WIM system using the two previously mentioned methods is considered. The vehicle is assumed to have two axles and crosses the bridge with a constant speed. The bending moment histories of the instrumented bridge at various sections subjected to a passage of the vehicle are simulated from the vehiclebridge interaction model as derived in (10). To reduce the simulation error, Newmark’s method with a fine time interval of 0.001 second is adopted to numerically determine bridge bending moments. Based on these obtained moment responses, the axle weights of the passing vehicle are estimated using either Method I or Method II. The weight estimations from the two methods are extensively investigated and the obtained results are compared under various conditions of vehiclebridge parameters. To quantify the estimation accuracy of the methods, the weight estimation errors for the front and rear axles are, respectively, defined as
It is noted that the weighting matrices, B, defined in (16) for the Method I and (18) for the Method II are both set to identical matrix to prevent the bias comparison. It is also noted that the regularization parameter, as required by (18) for the Method II is simply set to 1.0.
The parameters for the vehicle and bridge WIM system are listed in Table 1. The bridge structure is modeled as a simply supported beam having span length of 10 m. Its properties are approximated from a real concrete bridge in Thailand. It is noted that the first five natural frequencies of this bridge are computed to be 6.7, 26.8, 60.2, 107.1, and 167.3 Hz. The vehicle having gross weight of 250 kN is considered. Its dynamic properties are obtained from model identification by testing of a real twoaxle 10wheel truck.

Figure 4 shows typical bending moment histories of the bridge at , and under a passage of vehicle with speed of 15 m/s, respectively. Using these simulated bridge moment responses as the input, the axle weights of the vehicle can be estimated using Methods I and Method II as previously outlined and are shown in Figures 5(a) and 5(b), respectively. It should be noted that Method I yields directly the two axle weights of vehicle. On the other hand, Method II firstly yields the timevarying magnitudes of axle loads which can then be averaged to estimate the corresponding axle weights of vehicle. The figures clearly reveal that both methods can accurately estimate the front and rear axle, weights of vehicle. In this particular case, the estimation errors of the front axle, rear axle, and gross weights are, respectively, about 3.9, 4.4, and 4.3 using Method I and about 2.9, 3.7, and 3.5 using Method II.
(a) Method I
(b) Method II
4.1. Effects of Bridge Discretization and Sampling Frequency
The effects of bridge discretization and sampling frequency on the accuracy of the weight estimation methods are considered. The vehicle moving on the bridge at a constant speed of 15 m/s under various bridge discretization refinements and various sampling frequencies of bridge bending moments is simulated. Based on the obtained bridge moments at three sections, that is, and the axle weights of the vehicle is estimated using Method I and Method II. Table 2 lists the estimation errors of front axle, rear axle, and gross weights of vehicle from the two methods. In the table, the simulated bending moments of bridge are sampled with the sampling frequency varied from 20 to 1,000 Hz. To estimate the axle weights of vehicle using Method II, the bridge structure is discretized into 4, 8, 12, and 16 beam elements. It should be noted that the bridge discretization for Method I is not necessary since the method utilizes the continuous influence line function. In the table, the bridge surface roughness of level 2, which induces the dynamic in axle loads around 10%, is also assumed.

Employing Method II, it is obviously found from the table that the accuracy of axle weight estimations is not affected by discretization refinement of bridge structure if it is discretized by more than 4 elements. The table also indicates that the accuracy of axle weight estimations by the two methods is significantly influenced by sampling frequency. The weight estimation errors become larger when the sampling frequency is smaller. However, the weight estimation errors obtained from both methods are rather constant if the sampling frequency is faster than around 50 Hz. Therefore, throughout this study, the number of bridge elements is set to 8 while the sampling frequency is fixed at 500 Hz to guarantee the highest accuracy of the estimation methods. These imply that the bridge vibrations up to the natural mode are taken into account. It is noted that these settings are used only for weight estimation methods. The vehiclebridge interaction simulation of (10) still employs a very fine time interval of 0.001 second (1000 Hz) with 16 beam elements to accurately simulate the actual bending moments of the bridge WIM under a passage of vehicle.
4.2. Effects of Vehicle Speed and Bridge Surface Roughness
The effects of vehicle speed and bridge surface roughness on the accuracy of the weight estimation methods are investigated. The practical range of vehicle speed from 1 to 30 m/s is considered. The magnitude of bridge surface roughness, is varied and is classified into 6 roughness levels ranging from 0 (smooth surface) to 5 (very rough surface) according to the obtained dynamic characteristic of the simulated axle loads. In general, roughness levels of 0, 1, 2, 3, 4, and 5 indicate the averaged dynamic participations in each axle loads of vehicle moving at 15 m/s around 0%, 5%, 10%, 15%, 20%, and 25%, respectively. The vehicle moving on the bridge having different roughness levels at various speeds is simulated. Based on the obtained bridge moments at three sections, that is, and the axle weights of vehicle are estimated using Method I and Method II. Figures 6(a) to 6(c) plot the estimation errors of front axle, rear axle, and total weights of vehicle from the two methods under various vehicle speeds and bridge roughness levels. For a better visualization, the weight estimation errors from Method I and Method II are separately plotted and shown in the left and right figures, respectively. These figures reveal that the weight estimation errors from both methods exhibit similar characteristic under variations of vehicle speed and bridge roughness level. It is found that their estimation errors tend to increase as the roughness level increases due to larger fluctuation of axle loads. Although the faster vehicle speed also induces larger fluctuation of axle loads, the speed effects on the estimation errors for front (Figure 6(a)) and rear axles (Figure 6(b)) become different. The former becomes smaller while the latter becomes larger when the vehicle speed increases. Since the rear axle is much heavier than the front axle, the estimation errors of gross weight from both methods (Figure 6(c)) are almost the same as those of the rear axle (Figure 6(b)). Comparing the errors from the two estimation methods, it is found from these figures that Method II is superior for all considered ranges of vehicle speeds and roughness levels. In particular, the estimation errors of front axle, rear axle, and gross weights using Method II are, respectively, varied from to 12.53%, to 13.51%, and to 11.30% and using Method I are, respectively, varied from to 19.20%, 0.02% to 17.31%, and 0.01% to 13.71%. It is observed from Figure 6(c) that the estimation errors of gross weight of vehicle from the two methods can be however controlled to be within if the surface roughness is kept below level 4, regardless to the vehicle speed.
(a) Front axle weight
(b) Rear axle weight
(c) Gross weight
Since Method II provides not only the estimated axle weight of vehicle but also its dynamic axle loads, it is therefore interesting to investigate the accuracy of identified dynamic axle loads of the method. To do so, the load estimation errors of rear axle defined by norm of the load error, , are plotted in Figure 7 under various vehicle speeds and bridge roughness levels. It is noticed from the figure that the load estimation error is affected by both vehicle speed and bridge roughness level. The error increases as the vehicle speed or bridge roughness level increases. Based on considered ranges of speed and roughness level, the estimation errors of to 13.49% are observed.
4.3. Effects of Number of Measuring Sections and Noise Level
The influences of number of measuring sections and noise levels in the input signals on weight estimation accuracy resulting from the two methods are investigated. The vehicle moving at a constant speed of 15 m/s on the bridge having roughness of level 3 is simulated. Based on the obtained bridge moments at various section arrangements as in Table 3, the axle weights of vehicle are estimated using Method I and Method II. It is noted that the number of the measuring sections of 1, 3, 5, 7, or 9 sections is considered. To study the noise effect, the obtained moment signals at all sections are assumed to be polluted with 5% to 50% white noise. Figures 8(a) to 8(c) show the estimation errors of front axle, rear axle, and gross weights of vehicle from the two methods under various number of measuring sections and noise levels. The results indicate that the influence of number of measuring sections and noise levels on weight estimation accuracy resulting from the two methods are very small, especially for rear axle and gross weights. Although both estimation methods are expected to be very robust against noise effect since their computations inherently contain time averaging procedure, the significant effect of noise on the weight estimation of front axle is observed when small number of measuring sections, that is, 1 or 3, is employed. It should be noted that both estimation methods can provide rather accurate weight estimation using only midspan bending moment as the input even though the number of vehicle axle is 2. This implies that the application of the two methods can be extended to the case where multiple vehicles are simultaneously presented on the bridge using only limited number of measuring sections of bridge bending moments. It is observed that the estimation errors of rear axle and gross weights of vehicle from both Method I and Method II can be kept smaller than 5% for considered range of noise levels using only one measuring section of bridge. It is also observed that increasing the number of measuring sections beyond 3 does not significantly improve the accuracy of weight estimation methods. Comparing between the two estimation methods, it is found that Method II yields slightly smaller errors than Method I for all considered number of measuring sections and noise levels.

(a) Front axle weight
(b) Rear axle weight
(c) Gross weight
4.4. Effects of Axle Spacing and Axle Weight Distribution
The effects of axle spacing and axle weight distribution of vehicle on the accuracy of weight estimation methods are investigated. The axle spacing of vehicle from 2 to 15 m and the axle weight distribution of vehicle defined by ratio of front axle weight to gross weight are varied from 20% to 80%. The vehicle having different axle spacing and axle weight distribution moving on the bridge at a constant speed of 15 m/s is simulated. The bridge surface is assumed to have the roughness of level 3. Based on the obtained bridge moments at three sections, that is, and the axle weights of vehicle are estimated using Method I and Method II. Figures 9(a) to 9(c) show the estimation errors of front axle, rear axle, and gross weights of vehicle from the two methods under various axle spacing and axle weight distribution. The obtained results in Figure 9(c) clearly indicate that the estimation error of gross weight is slightly affected by both axle spacing and axle weight distribution. However, their effects on the estimation errors of front and rear axle weights as in Figures 9(a) and 9(b) are found to be significant. These errors increase as the weight of the considered axle reduces. This is mainly because the estimation error is computed in terms of percentage error with respected to the corresponding actual weight of the axle. Therefore, with the same amount of weight error, an axle with lighter weight would give a higher percentage of error. For the gross weight estimation, the errors of gross weights from 3.9% to 5.5% and from 3.1% to 4.6% are observed using Method I and Method II, respectively. Comparing their errors, it is found that Method II yields slightly smaller errors than Method I especially for the lower axle weight distribution.
(a) Front axle weight
(b) Rear axle weight
(c) Gross weight
4.5. CPU Processing Time
The CPU processing time is also investigated to compare the computing speeds of Method I and Method II. Based on numerical simulation test results, the processing times required by both methods are listed in Table 4. In this study, a personal computer with Intel Core 2 Duo 2.4 GHz and 2.0 GB RAM is employed and the listed processing times are normalized by 100 data points of input bridge moment histories. The results show that Method I requires about 1–3 seconds to complete the axle weights estimation. While Method II spends up to 4–12 seconds to complete the axle loads identification. Comparing between both methods, it is clearly found that Method I exhibits about 4 times shorter processing time than Method II. However, it should be noted that Method II provides not only the axle weight information of vehicle but also its dynamic load histories.

5. Conclusions
The effectiveness of vehicle weight estimations from bridge weighinmotion system is studied. The measured bending moments of the instrumented bridge at selected sections under a passage of the vehicle are numerically simulated and are used as the input for vehicle weight estimations. Two weight estimation methods assuming constant magnitudes (Method I) and timevarying magnitudes (Method II) of vehicle axle loads are investigated. Their estimation accuracy are evaluated and compared under various parameters of vehiclebridge system.
Based on the simulation results, the minimum number of bridge discretization of 4 and the minimum sampling frequency of 50 Hz are observed. It is also found that, among many considered parameters, the vehicle speed and surface roughness seem to have stronger effect on the accuracy of the two estimation methods than others. However, the estimation errors of the gross weight of vehicle can be controlled to be within if the surface roughness is kept below level 4, regardless to the vehicle speed.
It is also found from the effectiveness comparison between the two estimation methods under various vehicle and bridge conditions that Method II can provide better weight estimation than Method I for almost of the considered cases. In addition, it provides dynamic axle loads of vehicle. However, it exhibits about four times slower speed of computation than Method I.
Appendix
Acknowledgment
The authors are very grateful to the Chulalongkorn University for the 90th year research grant for this project.
References
 F. Moses, “Weighinmotion system using instrumented bridges,” ASCE: Journal of Transportation Engineering, vol. 105, no. 3, pp. 233–249, 1979. View at: Google Scholar
 M. Kriss, Weigh in motion prediction method, M.S. thesis, Case Western Reserve University, Cleveland, Ohio, USA, 1979.
 N. Gagarine, Advances in weighinmotion with pattern recognition and prediction of fatigue life of highway bridges, Ph.D. thesis, Department of Civil Engineering, University of Maryland, College Park, Md, USA, 1991.
 European Commission 4th Framework Programme Transport, “Weighinginmotion of Axles and Vehicles for Europe (WAVE),” Report of Work Package, Bridge WIM systems (BWIM), University College Dublin, 2001. View at: Google Scholar
 S. K. Leming and H. L. Stalford, “Bridge weighinmotion system development using superposition of dynamic truck/static bridge interaction,” in Proceedings of the American Control Conference, vol. 1, pp. 815–820, 2003. View at: Google Scholar
 C. Rowley, A. Gonzalez, E. Obrien, and A. Znidaric, “Comparison of conventional and regularized bridge weighinmotion algorithms,” in Proceedings of the 5th International Conference on Heavy Vehicles (ICWIM '08), pp. 271–282, Paris, France, May 2008. View at: Google Scholar
 S. A. Ansari, Investigation of improved accuracy limits for dynamic bridge weighinmotion system, M.S. thesis, Department of Civil Engineering, University of Maryland, College Park, Md, USA, 1990.
 M. I. S. Elhelbawey, Definition of improved modeling procedures for the analysis of bridge structures under truck loading using weighinmotion, Ph.D. thesis, Department of Civil Engineering, University of Maryland, College Park, Md, USA, 1991.
 D. R. Schelling, M. A. Sahin, M. I. Elhelbawey, S. A. Ansari, and T. E. Huang, “Definition of impact factors for bridge structure using truck loading,” Research Study AW087329046, Maryland Department of Transportation, Brooklandville, Md, USA, 1992. View at: Google Scholar
 S. S. Law, T. H. T. Chan, and Q. H. Zeng, “Moving force identification: a time domain method,” Journal of Sound and Vibration, vol. 201, no. 1, pp. 1–22, 1997. View at: Publisher Site  Google Scholar
 S. S. Law, T. H. T. Chan, and Q. H. Zeng, “Moving force identification a frequency and time domains analysis,” Journal of Dynamics Systems, Measurement, and Control, vol. 121, no. 3, pp. 394–401, 1999. View at: Google Scholar
 X. Q. Zhu and S. S. Law, “Identification of vehicle axle loads from bridge dynamic responses,” Journal of Sound and Vibration, vol. 236, no. 4, pp. 705–724, 2000. View at: Publisher Site  Google Scholar
 L. Yu and T. H. T. Chan, “Moving force identification from bending moment responses of bridge,” Structural Engineering and Mechanics, vol. 14, no. 2, pp. 151–170, 2002. View at: Google Scholar
 X. Q. Zhu and S. S. Law, “Moving loads identification through regularization,” Journal of Engineering Mechanics, vol. 128, no. 9, pp. 989–1000, 2002. View at: Publisher Site  Google Scholar
 S. S. Law, J. Q. Bu, X. Q. Zhu, and S. L. Chan, “Vehicle axle loads identification using finite element method,” Engineering Structures, vol. 26, no. 8, pp. 1143–1153, 2004. View at: Publisher Site  Google Scholar
 S. S. Law and Y. L. Fang, “Moving force identification: optimal state estimation approach,” Journal of Sound and Vibration, vol. 239, no. 2, pp. 233–254, 2001. View at: Publisher Site  Google Scholar
 T. Pinkaew, “Identification of vehicle axle loads from bridge responses using updated static component technique,” Engineering Structures, vol. 28, no. 11, pp. 1599–1608, 2006. View at: Publisher Site  Google Scholar
 T. Pinkaew and P. Asnachinda, “Experimental study on the identification of dynamic axle loads of moving vehicles from the bending moments of bridges,” Engineering Structures, vol. 29, no. 9, pp. 2282–2293, 2007. View at: Publisher Site  Google Scholar
 X. Q. Zhu and S. S. Law, “Moving forces identification on a multispan continuous bridge,” Journal of Sound and Vibration, vol. 228, no. 2, pp. 377–396, 1999. View at: Google Scholar
 X. Q. Zhu and S. S. Law, “Orthogonal function in moving loads identification on a multispan bridge,” Journal of Sound and Vibration, vol. 245, no. 2, pp. 329–345, 2001. View at: Publisher Site  Google Scholar
 X. Q. Zhu and S. S. Law, “Moving load identification on multispan continuous bridges with elastic bearings,” Mechanical Systems and Signal Processing, vol. 20, no. 7, pp. 1759–1782, 2006. View at: Publisher Site  Google Scholar
 T. H. T. Chan and D. B. Ashebo, “Theoretical study of moving force identification on continuous bridges,” Journal of Sound and Vibration, vol. 295, no. 3–5, pp. 870–883, 2006. View at: Publisher Site  Google Scholar
 T. H. T. Chan and D. B. Ashebo, “Moving axle load from multispan continuous bridge: laboratory study,” Journal of Vibration and Acoustics, Transactions of the ASME, vol. 128, no. 4, pp. 521–526, 2006. View at: Publisher Site  Google Scholar
 “Optimization toolbox for use with MATLAB,” The MathWorks, Inc., Natick, Mass, USA, 2001, http://www.mathworks.com/. View at: Google Scholar
Copyright
Copyright © 2009 Teerachai Deesomsuk and Tospol Pinkaew. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.