Some examples of SHM and parameter estimation soft- ware are DIAMOND ....
parameter estimation, the SAP2000 model is invoked by MATLAB for calculating
the ... spring and 2-node link element stiffness in the FE model. Any com- bination
...
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier’s archiving and manuscript policies are encouraged to visit: http://www.elsevier.com/authorsrights
Author's personal copy
Advances in Engineering Software 67 (2014) 99–110
Contents lists available at ScienceDirect
Advances in Engineering Software journal homepage: www.elsevier.com/locate/advengsoft
Automated finite element model updating of full-scale structures with PARameter Identification System (PARIS) Masoud Sanayei a,⇑, Peeyush Rohela b,1 a b
Department of Civil and Environmental Engineering, Tufts University, Medford, MA 02155, USA Alfred Benesch & Company, 50 Redfield Street, Suite 102, Boston, MA 02122, USA
a r t i c l e
i n f o
Article history: Received 16 April 2013 Received in revised form 18 August 2013 Accepted 8 September 2013 Available online 29 September 2013 Keywords: Finite element model updating Parameter estimation Inverse problems Full-scale structures Application Programming Interface (API) PARameter Identification System (PARIS)
a b s t r a c t This paper presents a software framework, PARIS (PARameter Identification System), developed for automated finite element model updating for structural health monitoring. With advances in Application Programming Interfaces (API) for modern computing, the traditional boundaries between different standalone software packages hardly exist. Now complex problems can be distributed between different software platforms with advanced and specialized capabilities. PARIS takes advantage of the advancements in the computing environment and interfacing capabilities provided by commercial software to systematically distribute the structural parameter estimation problem into an iterative optimization and finite element analysis problem across different computing platforms. Three validation examples using simulated nondestructive test data for updating full-scale structural models under typically encountered damage scenarios are included. The results of model updating process for realistic structural models and their systematic treatment provide enhanced understanding of the aforementioned parameter estimation process and an encouraging path towards its feasible field application for structural health monitoring and structural condition assessment. Ó 2013 Elsevier Ltd. All rights reserved.
1. Introduction The American Society of Civil Engineers’ 2013 Infrastructure Report Card [1] rated America’s infrastructure at D+ and its existing bridges at C+ with an estimated total investment of 3.6 trillion dollars needed for infrastructure maintenance by 2020. The design life of these structures is determined based on the applicable design codes. Through their service life, these structures may experience impacts of climate variation, seismic events, corrosive environment, and other unexpected transient and long term loading, leading to structural deterioration. Structural health monitoring (SHM) is increasingly gaining emphasis in response to the continually deteriorating state of our transport infrastructure. SHM has the potential to be used in conjunction with the current methods of inspection to add a layer of confidence in objective decision making with regards to structural management and maintenance. Unlike limited visual inspections, it can provide a continuous stream of measured data to help benchmark structural condition and evaluate long term trends and performance. There is considerable volume of published research in the field of SHM based on a wide array of approaches to achieve the same ⇑ Corresponding author. Tel.: +1 (617) 627 4116; fax: +1 (617) 627 3994. E-mail addresses:
[email protected] (M. Sanayei),
[email protected] (P. Rohela). 1 Former graduate student at Tufts University. 0965-9978/$ - see front matter Ó 2013 Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.advengsoft.2013.09.002
end goal. Researchers at the Los Alamos National Laboratory have published state-of-the-art reports of on damage identification [2], damage prognosis [3], and health monitoring [4] of structures. ASCE recently published a book with comprehensive set of case studies of structural health monitoring (Editor: Catbas et al. [5]), This book reviewed contemporary Structural-Identification (St-Id) methods under research and organized them into distinct steps like modeling, data processing, and parameter identification. Techniques used for damage assessment are broadly classified as global and local. A review of global SHM and local non-destructive evaluation methods and their shortcomings is provided by Chang et al. [6] and Chang and Liu [7] respectively. Local SHM methods like acoustics and electromechanical impedance based SHM methods, often sophisticated and precise, require pre-identification of damage and are most effective in localized regions of damage. Global methods, however, use global characteristics of a structure’s response to determine either the extent or location of damage or both in a structure. Mengelkamp and Fritzen [8] combined the local ultrasonic wave propagation method with the global stochastic subspace fault detection method (SSFDM) to investigate and detect artificially introduced cracks and delaminations in a stiffened aluminum plate and a GFRP plate. SHM techniques are further categorized as model based and non-model based. There are various methods for the non-model based approach. As an example, Farrar and Worden [9] discussed recent research that identifies SHM problem primarily as a
Author's personal copy
100
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
Statistical Pattern Recognition problem. Deraemaeker et al. [10] used output-only vibration measurements under variable environmental conditions for damage detection using a numerical example of a bridge. The feasibility of using self organization and learning capabilities of neural networks for structural damage assessment has also been explored by researchers like Wu et al. [11]. A set of model based methods for SHM relies on observing changes in the structure’s stiffness, mass, and damping parameters to correlate predicted and measured response. Parameter estimation is the process of updating the parameters of an a priori mathematical model to correlate its predicted response with the measured data at selected observation points for a given set of excitations. This is called the inverse problem and is unlike the forward analysis problem, which is to determine response of a structure to a known excitation as a function of known physical parameters of the model. An acceptable method for solving input–output inverse problems is by minimizing the residual between the predicted response from a mathematical model and the measured response of the actual system constituting the observed data set. The process thus involves starting with an a priori mathematical model based on certain assumptions of the initial values of the parameters. Key structural parameters of the model are first considered to be unknown and then iteratively updated for minimizing the residual and updating the unknown parameters until their best fit values are achieved. The updated parameter estimates are reflective of actual state of the structure and can thus be used to assess the level of damage in the structure. Catbas and Aktan [12] discussed several promising indices and experimental constraints for condition and damage assessment. Focusing on bridges, they stated that a successful condition and damage assessment can be achieved with the appropriate set of experiments and judicial use of indices in structural identification framework. Kaouk and Zimmerman [13] adopted the Minimum Rank Perturbation Theory (MRTP) index for locating damage and determining its extent from measured eigenvalues and eigenvectors. Using vibrational test data, Doebling et al. [14] presented the problem of determining local stiffness parameters as a welldetermined linear least square problem based on the decomposition of the flexibility matrix. Wang et al. [15] presented a structural damage identification algorithm for Damage Signature Matching (DMS) using change in static displacement response as a function of perturbation in stiffness matrix and changes in natural frequency. In order to consolidate research efforts in SHM further, the IASC-ASCE Task Group created a series of benchmark problems. The first phase of this study was based on using only simulated acceleration response data from the scaled physical model of a four story building frame. Caicedo et al. [16] used the Natural Excitation Technique (NExT) and Eigensystem Realization Algorithm (ERA) for identification of modal parameters for damage assessment on the IASC-ASCE benchmark problem. The reliance of comprehensive infrastructural inspection and monitoring methods for objective condition evaluation on advanced technologies demands modern computational tools capable of analyzing larger structural systems more efficiently and accurately. Some examples of SHM and parameter estimation software are DIAMOND and FEMtools. DIAMOND [17] from Los Alamos National Laboratory is a software toolbox with several vibration based damage identification algorithms for vibration tests simulation, modal data analysis, finite element correlation, and comparison of both linear and nonlinear damage identification techniques. Its successor, DIAMOND II [18] is a collection of algorithms for treatment of SHM as a statistical pattern recognition problem. FEMtools [19], a commercial FE model updating software that interfaces with FEA packages like ABAQUS and NASTRAN is developed by Dynamic Design Solutions. The development of a software
platform for full-scale finite element (FE) model based parameter estimation using non-destructive test (NDT) data is the subject matter of this paper, addressing important aspects of FE model calibration for structural health monitoring and structural conditions assessment. PARameter Identification System (PARIS) is a MATLABÒ based program that integrates seamlessly with the commercially available Finite Element Analysis (FEA) program, SAP2000Ò, to evaluate structural health at elemental level by observing a structure’s global response in a NDT. A unique feature of this program is utilization of interactive functionality of existing commercial software packages, SAP2000 and MATLAB. SAP2000 serves as the Finite Element Analysis (FEA) solver and with an integrated Graphical User Interface (GUI) for 3D model creation and validation. MATLAB is a fourth generation programming environment with exhaustive function optimization options. The APIs from SAP2000 developers enables MATLAB to use SAP2000 as a slave program for FEA. Additionally, visualization capabilities of the Graphical User Interface (GUI) of SAP2000 are used for FE model creation, verification, analysis, and post processing. This combination facilitated development of a platform for automated FE model updating of full-scale structures in the form of PARIS computer program. PARIS calibrates the FE model iteratively and automatically by minimizing the residual between the predicted response of the FE model and the measured response from NDT data. PARIS, as a research freeware, is posted at the SHM research website of Civil and Environmental Engineering Department at Tufts University at http://engineering.tufts.edu/cee/shm/software.asp. This research is based on the parameter estimation formulations developed by researchers at Tufts University for PARIS (PARameter Identification System). This work is also divided into two distinct phases. Phase-I is the computer program development with capabilities for full-scale bridge FE model updating. The program itself uses Application Programming Interface (API) of SAP2000 for two-way automatic data exchange between MATLAB and SAP2000. Phase-II includes testing the program with simulated NDT data to validate its functionality. Model updating using simulated NDT data for FE models of full-scale structures falls under Phase-II. 2. PARIS (PARameter Identification System) PARIS (PARameter Identification System) is the acronym for the custom MATLAB based computer program for parameter estimation and FE model updating. Fig. 1 illustrates the FE model updating process in PARIS version 13.0. The flowchart connectors labeled as API describe two-way automatic data exchange between MATLAB and SAP2000. This distinctive feature allows SAP2000 to be used as a slave program to MATLAB based PARIS during iterative stages of FEA and updating model parameters in the parameter estimation process. Parameter estimation process, as shown in Fig. 1, begins with an initial FE model creation in SAP2000 which is based on the preliminary guess of unknown structural parameters. In an a priori model, these parameters are rigidities and mass properties of elements. Typical stiffness parameters of a FE model are EA (axial rigidity), EI (bending rigidity) and GJ (torsional rigidity) for frame elements, E (modulus of elasticity) for shell and solid elements, and kX (translational stiffness) and kh (rotational stiffness) for joint springs. For dynamic systems, mass parameter in a FE model is m (mass) of the finite elements. A change in any of the structural parameters used for defining the a priori model is defined as ‘‘structural damage’’. The discrepancy between the predicted response from the FE model and measured simulated NDT data is quantified as the residual. This residual is called error function. Various error functions
Author's personal copy
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
101
Fig. 1. PARIS flowchart.
are stacked for multi-response FE model updating. A scalar objective function is formed using the error functions and is optimized in MATLAB. Unknown parameter values are updated iteratively such that the updated model in one cycle becomes the initial model for the next. In each iteration of the optimization process for parameter estimation, the SAP2000 model is invoked by MATLAB for calculating the FE responses. The iterations continue until all pre-defined convergence criteria are met. Successful convergence indicates that the physical model closely assimilates the real structure’s behavior. The following paragraphs describe various mathematical concepts and formulations used in PARIS in the sequence in which they are encountered in the FE model updating process.
eðpÞ ¼ qpredicted qmeasured
ð1Þ
where ‘q’ denotes a response quantity which can either be strain or translation and rotation under a set of known static loads or modal excitation. Static and modal error functions are further subdivided as stiffness and flexibility based error functions. In essence, the stiffness based error functions measure the residual between applied forces and the flexibility based error functions quantify residual between measured displacements. Static strain error function can also be categorized under flexibility based error function as strains are related directly to displacements. 2.3. Static stiffness (SS)
2.1. Finite element model creation and analysis options The model updating capabilities of PARIS include models comprised of frame, quadrilateral shell, and cuboid solid elements using static and/or modal data. Its usefulness lies in availability of advanced FE types for modeling. While solid elements may be better suited for modeling problems without geometrical oversimplification, shell elements, on other hand, provide greater convenience in modeling in addition to reduced analysis time. The model updating options have also been extended to update joint spring and 2-node link element stiffness in the FE model. Any combination of active DOFs for analysis can also be selected by the user.
Sanayei and Onipede [20] developed the static stiffness error function. They showed that the predicted displacement response from applied static test loads at one subset of degrees of freedom (DOF) and measured displacement response at another subset of DOF could be used to detect damage in the structures at the element level. Based on the partitioning of the force–displacement relationship as per measured and unmeasured displacements given by Eq. (2), the static stiffness error function is written as Eq. (3),
fa fb
¼
K aa K ba
K ab K bb
ua ub
ð2Þ
2.2. Error functions
1 eSS ðpÞ ¼ ðK aa K ab K 1 bb K ba Þua þ K ab K bb fb fa
Error functions are a measure of residual between the predicted response from the FE model and experimental data. They are formulated in terms of unknown parameters of the a priori FE model and they serve to estimate its unknown parameters. Generally, an error function is expressed as the difference between analytical and measured physical quantities as (1),
where vectors ua and ub are vectors of measured and unmeasured displacements. Submatrices Kaa, Kab, Kba, and Kbb represent the partitioned terms of the stiffness matrix and fa and fb represent the partitioned applied load force vector depending on measured and unmeasured displacements at ‘a’ and ‘b’ subsets of degrees of freedom.
ð3Þ
Author's personal copy
102
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
2.4. Static flexibility (SF)
2.8. Scalar objective function
Sanayei et al. [21] developed the static flexibility error function. 1
1 eSF ðpÞ ¼ ðK aa K ab K 1 bb K ba Þ ðfa K ab K bb fb Þ ua
ð4Þ
It is based on the inverse force–displacement relationship and it compares predicted and measured displacements at a subset of DOFs using Eq. (4).
Sanayei and Saletnik [22] developed the static strain error function. It is defined as the difference between predicted and measured strains at selected observation points as,
eSSTR ðpÞ ¼ Ba K 1 fa ea
ð5Þ
where Ba is the mapping matrix from the relation e = Ba q and is used to calculate strains from displacements for frame elements. However, for shell elements in PARIS, a more direct approach has been adopted to calculate strains from FE model stresses using Hooke’s Law. SAP2000 calculates stresses at the four integration points for shell elements and extrapolates them to the corner nodes for output. Internally in PARIS, strain calculated at a particular node is averaged based on shell elements surrounding that node. 2.6. Modal stiffness (MS)
K/i ¼ kM/i
ð6Þ
In Eq. (6), K and M are the stiffness and the mass matrices, respectively. Vector /i represents the ith natural mode shape and ki is the square of the ith natural frequency. 1
ð7Þ
After partitioning Eq. (6) similar to Eq. (2) and condensing out unmeasured mode shapes, the modal stiffness error function is given by Eq. (7) as a residual of modal forces. Vector /ai is defined as the measured modal displacements at subset ‘a’ of DOFs at ith iteration. 2.7. Modal flexibility (MF) Sanayei et al. [24] developed modal flexibility based error function by condensing the characteristic equation written in terms of the flexibility matrix. Hjelmstad [25] also arrived at a similar modal flexibility based error function formulation by only partitioning the mass matrix at measured and unmeasured DOF without using condensation. Similar to static flexibility, the modal flexibility error function includes the inverse of stiffness matrix, K, in its formulation. Modal flexibility error function, as formulated by Sanayei et al. [24], is shown in Eq. (8).
eMF ðpÞi ¼ ½k2i Dab ðI ki Dbb Þ1 Dba þ ki Daa I/ai
ð9Þ
A declining value of J(p) over successive iterations denotes convergence.
The difference in the order of magnitudes of various types of parameters, measured quantities, and error functions used, can potentially give rise to numerical difficulties by making the inverse problem ill-conditioned. Two types of normalization methods are available in PARIS: parameter normalization and error function normalization. Parameters that define structural element properties are normalized with respect to their initial values to give them the same weight in parameter estimation. Normalizing parameter values to 1.0 by dividing them by initial estimated values is called parameter normalization. An estimated value of 0 indicates complete damage. Additionally, a statistical method for adjusting to measurement errors is to normalize an error function based on the observability of each measurement. Sanayei and DiCarlo [26] implemented this normalization scheme for all three static error functions in PARIS and demonstrated a higher confidence level in estimated parameters especially in presence of measurement errors.
JðpÞ ¼ eðpÞT R1 e eðpÞ
Sanayei et al. [23] developed the modal stiffness error function. The eigenvalue problem shown in Eq. (6) forms the basis of modal stiffness error function.
ðK ba ki M ba Þ/ai
JðpÞ ¼ eðpÞT eðpÞ
2.9. Normalization
2.5. Static strain (SSTR)
eMS ðpÞi ¼ ½ðK aa ki M aa Þ ðK ab ki M ab ÞðK bb ki M bb Þ
The quadratic scalar objective function J(p) defined in Eq. (9) is used to judge convergence or divergence of the iterative model updating process.
ð8Þ
The matrix D in Eq. (8) is the dynamic matrix given by the relationship D = K1M.
ð10Þ P
Using Eq. (10), in which e is the covariance matrix of measured response data, the static error functions in PARIS are both weighted and normalized depending on the level of measurement error. The results of this process approximate that of a maximum likelihood estimator. 2.10. Objective function minimization Minimizing J(p) using hill climbing method is based on the analytical derivatives of the error functions defined in Eqs. (3), (4), (5), (7), and (8). The analytical sensitivities are presented in Refs. [20– 24]. The general from of the sensitivities is defined based on the first derivatives of the error function as shown in Eq. (11),
SðpÞ ¼
@feðpÞg @p
ð11Þ
The minimization process results in newer estimate of unknown parameters (p + Dp) which take the place of initial estimates (p) for the next iteration. PARIS provides the option of minimizing J(p) by using analytical sensitivities formulated in Refs. [20–24]. The change in parameter values, Dp is obtained from the iterative solution of algebraically nonlinear system of Eq. (12) using the Gauss–Newton method.
SðpÞDp ¼ eðpÞ
ð12Þ
where S(p) is the analytical sensitivity matrix obtained by taking the derivative of the error function vector e(p) with respect to the unknown parameters. The sensitivity matrix S(p) is of size NM NUP, where NM is defined as the total number of measurements and NUP is the number of unknown parameters in the problem. The change in parameter values, Dp, in the ith iteration, is obtained from solving the over determined system of equations in (12) leading to (13),
piþ1 ¼ pi þ Dpi
ð13Þ
Author's personal copy
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
An alternative to using analytical sensitivity matrices is using numerical sensitivity matrices by utilizing the MATLAB Optimization Toolbox. It is successfully implemented by using the algorithm for constrained minimization of the algebraically non-linear multivariable error function, fmincon, available in MATLAB. It uses interior point algorithm which is suited for large and sparse problems as well as for small and dense problems (MATLAB Documentation). It satisfies upper and lower bound on parameters at all stages and can recover from undefined and infinite values if encountered in the optimization process. This algorithm requires the target function to be defined in terms of unknown parameters that returns a scalar value, which in PARIS is the scalar objective function J(p). 2.11. Convergence criteria Convergence criteria are required for terminating the optimization routine. The dissimilarity between the updated unknown parameter values and true values is reconciled within reasonable limits at this level. Convergence criteria used in PARIS are given by Eqs. (14) and (15).
piþ1 pi 6 relative Dp tolerance p1 JðpÞiþ1 JðpÞi 6 relativeDJðpÞtolerance JðpÞ1
ð18Þ
ð15Þ
1 qm ¼ q þ eq Rq 2
ð19Þ
In analytical sensitivity matrix based objective function minimization, the sensitivity matrix S(p) is post-multiplied by a boolean group transformation matrix, GT, to map the elements of the sensitivity matrix according to the unknown parameter groups defined by the user. This step reduces the size of S(p) from NM NUP to NM NPG, where NPG is the number of parameter groups. In the case of constrained function minimization using MATLAB Optimization Toolbox, fmincon, the objective function J(p) takes the input argument as a smaller vector of grouped unknown parameters pG of size NPG 1 rather than the original vector of unknown parameters p of size NUP 1. The number of unknown independent variables in the resulting multivariable optimization problem in both cases is thus significantly reduced in comparison to the original problem without grouping. 2.13. Multi-response parameter estimation Multi-response parameter estimation refers to the use of different sets of measurements for different load cases and subset of measurements as well as using different types of static and modal data with their corresponding error functions [27]. In PARIS, stacking is implemented by default using Eqs. (16) and (17) whenever a combination of different load cases and/or different error functions is used.
½SðpÞ ¼ ½ SðpÞT1
SðpÞT2
T
::: eðpÞTn g
T
::: SðpÞTn
The parameter estimation process is sensitive to the accuracy of measurements collected during a NDT. The sensors and instruments can either be erroneous due to defects in the manufacturing process or become faulty in their handling and installation. Normally, a bound on accuracy of a measuring device is stated by the manufacturer. This bound on the level of errors, however, can be exceeded due to difficulties in sensor installation and connection losses [28]. Although there is no way to accurately model actual measurement errors, the options to introduce uniformly or normally distributed errors in simulated data are available in PARIS to study the influence of measurement errors on estimated parameters. While uniform error distribution resembles a banded type of error with equal probabilities of occurrence throughout, the normal error distribution is a non-banded type of error with its probabilities of occurrence being higher closer to true values [29]. Eqs. (18) and (19) are used for contaminating simulated NDT data with proportional and absolute errors, respectively.
ð14Þ
2.12. Parameter grouping
eðpÞT2
2.14. Measurement errors
1 qm ¼ q 1 þ eq Rq 2
Relative Dp tolerance for convergence is the change in the value of unknown parameter value with respect to its initial value between two successive iterations. Similarly, relative DJ(p) tolerance is the change in the objective function value with respect to its initial value.
eðpÞ ¼ f eðpÞT1
103
ð16Þ ð17Þ
Observing different types of responses simultaneously to estimate a set of unknown parameters is an effective approach as it increases the amount of information used in the parameter estimation process by simultaneously using esf, esstr, emf or ess, ems error functions in groups.
In the equations above, q is the vector of measured forces or simulated response. eq represents the percentage level of measurement error in Eq. (18) whereas in Eq. (19) it represents absolute error with units as defined by the user. The vector Rq is fully populated either with uniformly distributed random numbers between 0 and 1 or normally distributed random numbers with 0 mean and standard deviation of 0.5. The symbol denotes element-wise vector multiplication. 3. Verification examples The IASC-ASCE SHM Benchmark Structure is the leading contemporary example of a framed building structure. Currently in its second phase (Dyke et al. [30]), it provides a platform for researchers to test their hypotheses on damage identification and assessment for structural health monitoring. The FHWA’s bridge inspection program has inspired the SHM community to develop and test novel tools and techniques for bridge monitoring and maintenance. Cardini and DeWolf [31] present, in one example, the use of strain data from a multi-girder composite steel bridge for long-term SHM. Farrar and Jauregui [32] performed comparative analysis of various damage based assessment method using modal data from a damaged and an undamaged bridge. Three full-scale FE model updating examples using simulated NDT data are presented to illustrate the range of capabilities of the PARIS program using Application Programing Interfaces (API) for automated and systematic parameter estimation. These examples are the IASC-ASCE scale model benchmark building comprised of frame and shell elements, girder-3 of the Powder Mill Bridge made with solid (concrete deck) and shell (steel girder) elements, and the complete model of Powder Mill Bridge using shell elements (concrete deck) and frame (steel girder) elements. The combination of program features covered in each example is shown in Table 1. This section contains the description of the structures, model updating problem formulation, and parameter estimates for all three examples. The first column in Table 1 presents various features of PARIS which are further expanded in column two. Verification examples using these features of the program are shown in columns three to five for the FE model of the IASC-ASCE SHM Benchmark Structure,
Author's personal copy
104
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
Table 1 PARIS verification examples. Feature
Type/description
IASC-ASCEScale model
Finite elements
Frame Shell Solid Spring Damage cases ? SS SF SSTR MS MF Maximum likelihood Errors in applied forces Errors in measurements
U U
Error functions
Normalization Measurement errors
Parameter grouping Stacking
1 U U U
2 U U U
3
4 U
U U
U
U
U
Girder3 PMB model
FullScale PMB model U U
U U U 1 U U U
2
U 1
2
U
U
U
U
U
U
U
U U U U
U
Fig. 2. IASC-ASCE SHM Benchmark Structure. U
U Table 2 Damage case summary for benchmark structure.
the Powder Mill Bridge (PMB) Girder-3 model, and the full-scale PMB model, respectively. 3.1. IASC-ASCE SHM Benchmark Structure An analytical model based on the Phase-I study of the IASCASCE SHM Benchmark Problem was updated using PARIS as the first verification example. The global geometry and member properties as well as damage scenarios 1 and 2 to simulate damage are taken from Black and Ventura [33] and Johnson et al. [34]. The FE model used in this research is slightly different than the actual analytical model. Originally, two finite element models were developed by IASC-ASCE SHM Task Group. The first model was a 12 DOFs shear-building model in which only two horizontal translations and one rotation was permitted for each floor. The second model had 120 DOFs which constrained all nodes at each floor only for the same horizontal translations and in-plane rotations [34]. The structural model used here has 1944 DOFs because of greater discretization of shell elements at the floor levels and also from the fact that all rotational and translational DOFs are made active in analysis. 3.2. Structure and model description Fig. 2 shows the FE model of the experimental steel frame housed in the Earthquake Engineering Research Laboratory at the University of British Columbia. It is a 4-story, 3.6 m high structure with two bays in orthogonal lateral directions with each bay measuring 1.25 m. The sections are made out of hot rolled grade 300 W steel with 300 MPa nominal yield stress. Steel sections that constitute the structure are B100 9 (columns), S75 11 (beams), and L25 25 3 (braces). The floor slab at each story is comprised of 4 square slab panels. Mass of each panel is 800 kg at the first floor and 600 kg at the second and the third floor. The fourth floor has 3 panels of 400 kg mass and 1 panel that has a mass of 550 kg. Each of the 4 slab panels measuring 1.25 m 1.25 m is modeled using 16 thin-shell elements. The columns are oriented with their weak axis (y) along the global X axis in the SAP2000 model. All beams connect rigidly to the columns while the braces are pin connected and thus behave as axial members. All columns are fixed at the base.
Damage case
Damage case description
Error function
1
100% axial rigidity (EA) loss in all 8 braces at first story
2
30% bending rigidity (EIzz and EIyy) in columns at Grid A at first story level Reduction in mass of the 550 kg slab panel at fourth story to 400 kg 40% Reduction in mass of all floor slabs and 100% EA loss in 8 braces on one side of the frame at Grid A
SS, SF, SSTR SS, SF, SSTR MS, MF
3 4
SS + MS
3.3. Damage scenarios PARIS was used to estimate the unknown stiffness and/or mass parameters for four simulated damage scenarios as summarized in Table 2. In first damage case, the first floor braces with 8 unknown EA parameters were grouped into 4 different unknown parameter groups. Error functions SS, SF, and SSTR were run individually. In the second damage case, 3 EIzz and 3 EIyy parameters of the base columns at Grid-A were paired as two groups. Error functions SS, SF, and SSTR were run again individually. In the third damage case in which mass of one slab panel is the unknown parameter, modal error functions were used one at a time to estimate the reduced mass of 16 shell elements making up that panel as one group. In the fourth damage case, both stiffness and mass parameters were estimated using SS and MS error functions simultaneously. The eight braces on one side of the frame (at Grid A) of the model were paired into 4 unknown stiffness parameter groups based on their story level. The four mass parameter groups comprised of four floor slab panels at each of the 4 story levels. 3.4. Load cases and measurement locations Diagnostic static test loads were applied under four load cases. Table 3 shows the location and magnitude of loads for all four load cases. The prefix ‘S’ refers to a story. Reference is made to grid lines from Fig. 2. Displacement and strain measurement locations for static load cases used with damage cases 1, 2, and 4 are listed in Table 4 and Table 5. For estimating the unknown mass parameters for damage cases 3 and 4, mode shapes 1 to 3 were utilized. Mode 1 and 2 are
Author's personal copy
105
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
bending modes in the global Y and the X directions, respectively. Mode 3 is the torsional mode about Z axis. Mode shapes 1 through 3 and their associated natural frequencies for un-damped vibrations are shown in Fig. 3.
Table 6 Damage cases 3 and 4: Displacement measurement locations for the Benchmark Structure. Damage case
Mode
Story – measurement location
3
1 2 3
S2, S4 – Columns at Grid 3-C
1
S2, S4 – Columns at Grid 1-A and Grid 3-A
Disp. direction X
Table 3 Applied static loads on the Benchmark Structure. Load case LC1 LC2 LC3 LC4
Story-load location S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col. S1,S2,S3,S4-Col.
at at at at at at at at
Grid Grid Grid Grid Grid Grid Grid Grid
Direction 1-A 3-A 1-A 3-A 1-A 1-C 1-A 1-C
Load (kN)
X X X X Y Y Y Y
4
1.0 1.0 1.0 2.0 1.0 1.0 1.0 2.0
U U
2 3
Y
Z
U U U U U
U U
U
Table 4 Displacement measurement locations for the Benchmark Structure. Damage case
Load case
Story – measurement location
1
LC1 LC2 LC1 LC2 LC3 LC4 LC1 LC2 LC3
S1-Col. at Grid 1-C,3-C S1-Col. at Grid 1-C,3-C S1,S4-Col. at Grid 1-C S1,S4-Col. at Grid 1-C,3-C S1,S4-Col. at Grid 3-C S1,S4-Col. at Grid 3-A,3-C S2,S4-Col. at Grid 1-A,3-A S2,S4-Col. at Grid 1-A,3-A S2,S4-Col. at Grid 1-A,3-A
Disp. direction X
2
4
U U U U U U U
Y
Z U
U U U U U
Fig. 4. Damage Case 1: EA estimates for the Benchmark Structure. U U
U U
U
Displacements were measured for modes 1, 2, and 3 at locations listed in Table 6. 3.5. Parameter estimation results and discussion
Table 5 Strain measurement locations for the Benchmark Structure. Damage case
Load case
Story – measurement location
1
LC1, LC2
2
LC1, LC2
S1-Col. at Grid 1-C S1-Col. at Grid 3-C S1-Beam b/w 1-B-C S1-Beam b/w 3-B-C S1-Col. at Grid 1-C S1-Col. at Grid 3-C S1-Beam b/w 2-3-A S1-Beam b/w 2-3-C S1-Col. at Grid 3-A S1-Col. at Grid 3-C
LC3, LC4
SG locations (mm) x
y
z
450 450 625 625 450 450 625 625 450 450
50 50 32 32 50 50 32 32 50 50
0 0 0 0 0 0 0 0 0 0
The parameter estimation results for the four damage cases described above are presented here. For damage case 1, the parameter estimates are presented in Fig. 4. The first bar represents the normalized value of initial guess of parameters. Bar number 2 is the true value of the parameter used to simulate NDT data. The last three bars are parameter estimates using simulated measurements for error functions SS, SF, and SSTR, respectively. This case shows that the first floor braces are identifiable using either the nodal displacements or member strains using simulated NDT data. The SS, SF, and SSTR error functions were all used independently and successfully estimating the EA parameters for damage case 1. For damage case 2, the parameter estimates are presented in Fig. 5 for updating EIzz and EIyy parameter groups for three base
Fig. 3. Fundamental vibration modes of the Benchmark Structure.
Author's personal copy
106
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
Fig. 5. Damage Case 2: EIzz and EIyy estimates for the Benchmark Structure. Fig. 6. Damage Case 4: EA and mass estimates for the Benchmark Structure.
columns. All three static error functions estimated the grouped bending rigidities along the major and minor axes of the columns. The convergence criterion for relative DJ(p) tolerance was maintained at 10E10 while that for relative Dp tolerance was 10E06 for all damage cases. This case shows that column moments of inertia are identifiable using any of SS, SF, and SSTR error functions when the NDT loads were applied to the structure to cause bending of the columns with unknown EI parameters and measurements were taken to observe these parameters. The accuracy of the parameter estimates relates to good observability of the unknown parameters by the the measurements at the selected nodes or members. For damage case 3, the mass parameter of one roof slab panel was updated from an initial guess of 550 kg to a final value of 400 kg using modal stiffness and modal flexibility error function. Modal stiffness error function converged faster in 7 iterations than the modal flexibility error function that took 13 iterations. This damage case shows that mass parameter is also identifiable using modal error function MS and MF independently. Changes in mass alter the modal frequencies of the structure. PARIS detected the change and updated the mass parameter appropriately to minimize the residual. For damage case 4, the combined stiffness and mass parameter updating capability of PARIS and the resulting parameter estimates are shown in Fig. 6. The parameter groups G1 through G4 are the unknown EA parameters for braces. Groups G5 through G8 represent unknown floor mass parameters with an estimated value equal to true value of 0.6. In this case, advantage of using a smaller number of unknown parameters was achieved by grouping the similar parameters together. the mass matrix and the stiffness matrix were updated iteratively in PARIS until the objective function value was reduced to a level below the convergence limit. At this stage, the predicted response measured well with the measured response from the simulated NDT data.
Table 7 Damage case 1 with stacked error functions. EA estimates
G1 G2 G3 G4 Iterations
Error function(s) SF
SSTR
SF + SSTR
1.36E03 1.80E04 2.60E04 4.70E04 29
2.40E04 4.21E05 1.90E04 8.55E03 26
9.80E04 1.37E03 3.60E04 4.10E04 24
Fig. 7a. Damage Case 2: Noisy parameter estimates w/o normalization.
3.6. Parameter estimation in presence of measurement errors Table 7 compares the EA parameter estimates and number of iterations required for FE model updating for damage case 1 shown in Fig. 4 in presence of 2% normally distributed proportional measurement errors in forces and displacements. Damage case 1 simulated a 100% axial rigidity loss in all 8 braces at the first story as shown in Table 2. The parameter estimates using SF and SSTR are shown, first separately and then stacked together. The merit of multi-response parameter estimation is seen in reduction of the number of iterations required when SF and SSTR error functions are stacked. Fig. 7a and 7b represent the parameter estimates for damage case 2 using static error functions in presence of 2% normally distributed proportional measurement error in forces and displacements. In this damage case, bending rigidity loss of 30% in columns at Grid A at the first story level was present. Fig. 7a shows high scatter in the parameter estimates in presence of measure-
Fig. 7b. Damage Case 2: Noisy parameter estimates with normalization.
ment errors. However, the EIzz and EIyy parameter estimates using statistical error function normalization were closer to the true parameter values in presence of measurement errors. The results in Fig. 7b demonstrate the robustness of the statistical normalization method proposed by Sanayei and DiCarlo [26]. The authors recommend the use of the statistical normalization method for parameter estimation. After successful completion of the IASC-ASCE SHM Benchmark Structure using four damage cases with and without measurement errors, a full scale bridge is studied first using shell and solid elements for modeling one composite girder and then the full bridge is modeled using shell and frame elements. Several damage cases using various error functions is used for parameters estimation and FE model updating.
Author's personal copy
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
Fig. 8. Powder Mill Bridge (PMB) in Barre, MA.
3.7. Powder Mill Bridge The Powder Mill Bridge (PMB) in Fig. 8 is a three span continuous steel girder bridge with a reinforced concrete deck in composite action. The bridge crosses the Ware River in the town of Barre, Massachusetts. This bridge was used for examining capabilities of PARIS for full-scale FE model updating. Two finite element models of the PMB were updated for assumed damage scenarios. The first model discussed below is of one of the continuous interior girders of the bridge. This model is made of solid and shell finite elements and is referred to as Solid-Shell model further in the text. The second model is the full bridge model which is made of shell and frame elements and is referred to as Shell-Frame model. Truck loads are applied as static loads in different load cases to simulate a truck crossing the bridge. The damage scenarios consider alteration in the stiffness of finite elements to simulate damage incurred in the deck and girders. 3.8. Powder Mill Bridge Girder-3 solid-shell model 3.8.1. Model description The first finite element model used with this example is of a single interior girder (Girder-3) of the bridge as shown in Fig. 9. This model is comprised of cuboid solid elements for bridge deck and haunch and quadrilateral plan thin-shell elements for web and flanges of the W920 238 steel girder. The support conditions at the pier cap and abutment are modeled with joint springs based on stiffness values from Sanayei et al. [35]. Solid elements used to model deck and haunch have reduced modulus of elasticity for concrete (Ec) in the negative moment region, which is considered to be 20% of the length of each span. 3.8.2. Load cases and measurement locations A test truck load of 353.6 kN (79.5 kips) is applied along a straight path centrally aligned over Girder-3 of the bridge. In order to simulate a crawl speed NDT, the truck location on the bridge was discretized into 13 static load cases as shown in Fig. 10. The truck load was applied as a 3-axle (3 2 point loads)
107
Fig. 9. PMB Girder-3 cross-section with contributory deck.
in each load case. The displacement and strain measurement locations were chosen to be the central node at the bottom flange at the center of the end spans of the bridge and at quarter points at the longer central span. Since PARIS is designed to accommodate different number of measurements in each load case, only observations from the span on which the truck was located were used with that load case, respectively. The PMB Girder-3 Solid-Shell model has 9,168 nodes which represent a total of 35,256 DOFs. The model was however analyzed in XZ-plane with two translational and one rotational DOF thus reducing the active DOFs to 27,504. Since the size of the stiffness matrix is dependent on the number of active DOFs, a reduction in the number of active DOFs in a model reduces computational time, thus speeding up the FE model updating process. 3.8.3. Damage scenarios This verification example included model updating in two simulated damage scenarios with no measurement errors. The damage scenarios and parameter estimation results are presented next. Damage case 1: In the first damage case, a loss of 20% horizontal translational stiffness (kXX) and rotational stiffness (khY) was considered to simulate joint deterioration. The damage was assumed to occur similarly at both the pier and abutment locations. The bottom flange of the girder has four shell elements across its width thus requiring joint springs at 5 points. All five spring stiffnesses were grouped together as a single parameter at each of the four support locations (two piers and two abutments). PARIS was used to accurately estimate the unknown stiffness values using all three static error functions. The stiffnesses of the joint springs contribute to the stiffness matrix of the structure and thus determine the structure’s response to loading conditions. It was possible to update the spring stiffnesses in PARIS through the use of SAP2000 APIs to access and modify stiffness values. In this case the static stiffness and static strain error functions estimated the unknown bending rigidities more accurately than the static flexibility error function. Mean of translational and rotational stiffness values at four support locations is shown in Fig. 11. SS and SSTR functions performed
Fig. 10. Test truck front axle locations on PMB Girder-3.
Author's personal copy
108
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
Fig. 12. PMB shell-frame model plan. Fig. 11. Bearing pad stiffness (kXX and khY) estimates for PMB Girder-3 model.
better in terms of accuracy of the parameter estimates. The number of iterations required for model updating for damage case 1 were 27 for the SS error function, 96 for SF error function, and 37 for the SSTR error function. Damage case 2: A loss of 20% of the girder’s stiffness was assumed to occur at the entire south span of 11.75 m length as a global deterioration case. Since girder is comprised of shell elements, the damage in finite element model would then be reflected in reduction of modulus of elasticity (Es) of steel. The 624 shell elements making up the south span of the girder in the FE model were grouped into 8 groups with 78 elements in each group. The elements were grouped along the length of the girder with equal number of elements in each group. Static strain error function was used to estimate 20% reduction in Es. The PMB Girder-3 Model was updated using only SSTR error function for damage case 2. The mean of estimated parameters values of Es for shell elements was 0.8035 from the expected value of 0.8 and had a 0.0148 standard deviation. Unlike translational displacement measurements, strain measurements are reference-independent making them a preferred measurable response quantity. Strain gauges are also a low cost option when compared to both contact and non-contact devices required for measuring displacements. One particular aspect that deserves attention is the ability to discretize the FE model to a satisfactory level that leads to a large number of elements in the model. Normally it is not feasible to identify all elements in a finite element with a fine mesh using a coarse set of measurements points. Grouping of similar unknown parameters allow solving of the inverse problem. At the same time when a localized damaged region is observed, it is possible to ungroup the parameters of those areas to zoom in damage severity at smaller locations. 3.9. Powder Mill Bridge shell-frame model 3.9.1. Structure and model description A full-scale FE model of the PMB comprised of frame and shell elements was updated for two damage scenarios. Frame elements were used to model the six main girders of the bridge as well as the diaphragms and outriggers at the north end of the bridge. Offset was applied to frame elements for modeling supports at the bottom flange of the girder. The shell elements constituted the concrete bridge deck. The shell elements were offset upwards by half the shell thickness to model the deck resting on top flange of the girders. The PMB Shell-Frame model is shown in plan in Fig. 12 with gridlines 1 through 6 representing girder numbers and S1 through S4 representing abutment (S1 and S4) and pier support (S2 and S3) locations. Darker longitudinal lines indicate girders and darker transverse lines indicate diaphragms. The static load cases were applied in a fashion similar to that for the solid-shell girder model with the addition of four test truck paths across the width of the bridge. Eight sets of truck loads were modeled in each of the five test truck paths on the bridge.
Fig. 13. EI33 estimates from SSTR error function for the full-scale PMB model.
3.9.2. Damage scenarios Damage case 1: A damage case that considers a theoretical 20% reduction in bending capacity of girders 4 and 5 was considered first. The entire length of these two girders was assumed to have a reduced EI33 value. Since each of the two girders in the model is comprised of 36 smaller frame elements connected at ends, the elements were grouped together into three groups where each group represents one span length. Starting with the initial guesses of EI33 of girders 4 and 5, the unknown moment of inertias converged to values very close to the true values of EI33 that were used for simulating NDT data. The results from model updating for shell-frame model using SSTR error function are shown in Fig. 13. The bending rigidity of frame elements representing girders 4 and 5 converged with a mean error of -1.33%. However maximum error in groups 3 and 6, which represents the north spans of the girders, was higher than for the center and south spans. This could be attributed to smaller deformations from increased stiffness at the north end of the bridge leading to lower sensitivity in estimating parameter values. Damage case 2: Damage to the bridge deck is a prominent type of damage in bridges. This type of damage can be caused by traffic loads, environmental factors including diurnal and seasonal temperature variation, and water infiltration through expansion joints or poor quality control during construction. To simulate this damage scenario, a cluster of 30 shell elements equivalent to deck surface of 44.04 m2 were grouped together to simulate damage. It represents a patch beginning at length 3.92 m from south end till 11.75 m length and lies on the right lane looking north. The damage level is taken to be 50% of the initial shell elasticity modulus (Ec). The concrete modulus of elasticity for all five deck groups was successfully estimated (see Fig. 14). Ec estimates for loss of 50% deck stiffness for a cluster of shell elements described under damage case 2 were quite accurate. The model updating process was automatically terminated when relative change in parameter values between two iterations became less than 10E12. It took 13 iterations for the convergence criteria to be met.
Author's personal copy
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
109
With these capabilities, PARIS significantly steers the use of FE model updating for structural health monitoring and assessment of in-service structures. Integration of commercial FEA packages with custom research software creates opportunities for validating and comparing several parameter estimation and FE model updating techniques. In development and deployment of frameworks such as this lies the benefit of utilizing modern technology for effective and thorough structural health monitoring and FE model calibration. 5. Future work Fig. 14. Ec estimates from SSTR error function for the full-scale PMB model.
The three examples presented here to demonstrate that the Application Programming Interfaces enable researchers to take advantage of the existing commercial programming and structural analysis software packages. PARIS has access to advanced programming with optimization mathematics as well as various finite elements for 3D modeling and analysis of large structures taking advantage of pre and post-processing tools. With these examples, the feasibility to perform finite element model updating of a 3D laboratory building model, a 3D bridge girder, and a full-scale 3D bridge model is demonstrated. The range of capabilities of the PARIS program using APIs include the utilization of error functions for use of NDT data such as static translations, rotations, and strains as well as dynamic frequencies and mode shapes for automated and systematic finite element model updating. Methods for updating full-scale structures with a large number of elements are presented using parameter grouping and stacking of various types of measurements for multi-response parameter estimation. Measurement errors are introduced to study the sensitivity of the inverse problem to the number, type, and location of loadings, measurements, and unknown parameters. Performing error sensitivity simulations prior to the use of NDT data is paramount in success of finite element model updating using measured test data of full-scale structures. 4. Conclusions PARIS provides a platform for automated FE model updating of full-scale structures. The PARIS program is capable of automatically calibrating an a priori FE model created in SAP2000 by reconciling its analytical response with the measured response. PARIS achieves this by automatically accessing the finite element analysis (FEA) solver of SAP2000 program by using the available SAP2000 Application Programming Interface (API) function library. Through the API interface, MATLAB gains access to SAP2000 elements library, analysis modules, and the GUI interface for pre-processing and post-processing, thus using SAP2000 as a slave program. Additionally PARIS is written in MATLAB which is a powerful programming platform with access to various mathematical functions, matrix operations, and optimization routines. Successful implementation of MATLAB’s readymade optimization routines for constrained function minimization aids in handling of larger parameter estimation problems. The performance of PARIS has been demonstrated using three examples with simulated test data. The program was used to update the analytical model of the IASC-ASCE SHM Benchmark Structure, the Powder Mill Bridge Girder-3 Solid-Shell model, and the full-scale Powder Mill Bridge Shell-Frame model. This study validates the feasibility of using static and modal NDT data for successful full-scale FE model updating as well as the potential of the newly developed statistical error function normalization using the covariance matrix of measured data in improving accuracy of the parameter estimates in presence of measurement error.
Phase I and II of this project are completed. Further development of PARIS computer program aligns with the following future research goals: (1) Studying robustness of suggested improvements to parameter estimation process in the areas of sensor placement and accounting for modeling and measurement errors, and (2) Utilizing PARIS for full-scale bridge model calibration and damage detection using field measured NDT data. References [1] ASCE. The American Society of Civil Engineers’ 2013 Infrastructure Report; 2013. . [2] Doebling SW, Farrar CR, Prime MB, Shevitz DW. Damage identification and health monitoring of structural and mechanical systems from changes in their vibration characteristics: a literature review. Los Alamos National Laboratory Report LA-13070-MS; 1996. [3] Farrar CR, et al. Damage prognosis: current status and future needs. Los Alamos National Laboratory report LA-14051-MS; 2003. [4] Soon H, et al. A review of structural health monitoring literature: 1996–2001. Los Alamos National Laboratory Report LA-13976-MS; 2004. [5] Catbas FN, Kijewski-Correa T, Aktan AE. Structural identification (St-Id) of constructed facilities: approaches, methods and technologies for effective practice of St-Id (2013). ISBN: 9780784411971 ASCE Books; 2013. [6] Chang Peter C, Flatau Alison, Liu SC. Review paper: health monitoring of civil infrastructure. Struct Health Monit 2003;2(3):257–67. [7] Chang Peter C, Chi Liu S. Recent research in nondestructive evaluation of civil infrastructures. J Mater Civil Eng 2003;15(3):298–304. [8] Mengelkamp G, Fritzen CP. Combination of global and local methods for structural health monitoring. In: Proceedings of the 24th international modal analysis conference (IMAC-XXIV); 2006. [9] Farrar CR, Worden K. Introduction to structural health monitoring. Phil. Trans R Soc A 2006; 365(1851): 303–15 [15 February 2007]. [10] Deraemaeker A, Reynders E, DeRoeck G, Kullaa J. Vibration-based structural health monitoring using output-only measurements under changing environment. Mech Syst Signal Process 2008;22(1):34–56. [11] Wu X, Ghaboussi J, Garrett JH. Use of neural networks in detection of structural damage. Comput Struct 1992;42(4):649–59. [12] Catbas F, Aktan A. Condition and damage assessment: issues and some promising indices. J Struct Eng 2002;128(8):1026–36. [13] Kaouk M, Zimmerman DC. Structural damage assessment using a generalized minimum rank perturbation theory. AIAA J 1994;32(4):836–42 [ISSN 00011452]. [14] Doebling SW, Peterson LD, Alvin KF. Experimental determination of local structural stiffness by disassembly of measured flexibility matrices. J Vibr Acoust 1998;120(4):949–57. [15] Wang X, Hu N, Fukunaga H, Yao ZH. Structural damage identification using static test data and changes in frequencies. Eng Struct 2001;23(2001):610–21. [16] Caicedo JM, Dyke SJ, Johnson EA. Natural excitation technique and eigensystem realization algorithm for Phase I of the IASC-ASCE benchmark problem: simulated data. J Eng Mech 2004;130(1). [17] Doebling SW, Farrar CR, Cornwell PJ. DIAMOND: a graphical user interface toolbox for comparative modal analysis and damage identification. In: Proc of sixth international conference on recent advances in structural dynamics, Southampton, UK; July 1997, p. 399–412. [18] Allen DW, Clough JA, Sohn H, Farrar CR. A software tool for graphically assembling damage identification algorithms smart structures and materials 2003: smart systems and nondestructive evaluation for civil infrastructures. Proceedings, vol. 5057, p. 138–44. [19] FEMtoolsÒ Dynamic Design Solutions N.V. (DDS). Interleuvenlaan 64, B-3001, Leuven, Belgium; 2012. [20] Sanayei M, Onipede O. Damage assessment of structures using static test data. AIAA J 1991;29(7):1174–9. [21] Sanayei M, Imbaro GR, McClain JAS, Brown LC. Structural model updating using experimental static measurements. ASCE J Struct Eng 1997;123(6):792–8. [22] Sanayei M, Saletnik MJ. Parameter estimation of structures from static strain measurements I: Formulation. ASCE, J Struct Eng 1996;122(5):555–62.
Author's personal copy
110
M. Sanayei, P. Rohela / Advances in Engineering Software 67 (2014) 99–110
[23] Sanayei M, Wadia-Fascetti S, McClain JAS, Santini EM. Structural parameter estimation incorporating modal data and boundary conditions. J Struct Eng 1999;125(9):1048–55. [24] Sanayei M, Wadia-Fascetti S, Arya B, Santini EM. Significance of modeling error in structural parameter estimation. Comput-Aided Civil Infrastr Eng 2001;16(1):12–26. [25] Hjelmstad KD. On the uniqueness of modal parameter estimation. J Sound Vibr 1996;192(2):581–98. [26] Sanayei M, DiCarlo C. Finite element model updating of scale bridge model using measured modal response data. In: Proceedings of the ASCE, structures congress 2009, Austin, Texas; 2009. [27] Bell ES, Sanayei M, Rao N. Model updating of UCF benchmark using PARIS. In: International modal analysis conference, IMAC XXV, Orlando, Florida; February 19–22, 2007. [28] Smith IFC. Multiple-model structural identification. Encyclopedia of structural health monitoring. Chichester, UK: Wiley. Reference IMAC-CHAPTER-2009001; 2009. p. 2199–2208. [29] Sanayei M, Onipede O, Babu SR. Selection of noisy measurement locations for error reduction in static parameter identification. AIAA J 1992;30(9):2299–309.
[30] Dyke SJ, Bernal D, Beck J, Ventura C. Experimental phase II of the structural health monitoring benchmark problem. In: Proceedings of the 16th ASCE engineering mechanics conference; 2003. [31] Cardini AJ, DeWolf JT. Long-term structural health monitoring of a multigirder steel composite bridge using strain data. Struct Health Monit 2009;8(1):47–58. [32] Farrar CR, Jauregui DA. Comparative study of damage identification algorithms applied to a bridge: I. Experiment. Smart Mater Struct 1998;7:704. [33] Black CJ, Ventura CE. Blind test on damage detection of a steel frame structure. In: 16th International modal analysis conference (IMAC XVI), Santa Barbara, California. Proceedings; ; February 2–5, 1998. p. 623–629. [34] Johnson EA, Lam HF, Katafygiotis LS, Beck JL. Phase I IASC-ASCE structural health monitoring benchmark problem using simulated data. J Eng Mech 2004;130(1):45. [35] Sanayei M, Phelps JE, Sipple J, Bell E, Brenner B. Instrumentation, nondestructive testing and finite element model updating for bridge evaluation using strain measurements. J Bridge Eng 2010;17(1):130–8.