Monitoring and identification of structural damages/ Konstrukciju pazeidimu matavimai ir prognozavimas.

Petkevicius, K. ; Volkovas, V.

1. Introduction

A reliable evaluation of structural integrity becomes especially important at design, manufacture and service stages in objects of increased risk. The analysis of standard structures and estimation of their functionality for a resource period usually are performing by regulations and norms, which are based on the huge theoretical and practical experience [1-3]. These regulations and norms however cannot be easily applied to items of unique structure without additional detailed and comprehensive analysis [4, 5]. This is why at present it is allocated to work of development and improvement of structure strength prediction methodologies and technologies. Such activity takes place in different areas such as civil engineering, transportation, power industry and others [6, 7].

Structural health of buildings can be supported by performing constant building maintenance by various and different means. With the help of monitoring systems the defects and damages, which occurring due to external and internal factors, changes in the building structure, ageing of utility infrastructure and technological equipment, are identified [8, 9]. The constant monitoring helps protect buildings against dangerous collapse phenomena, which threat the environment and people.

Nature of structural health monitoring programs depends on the functions of building and their age, size and configuration, connecting structures, environment conditions, available design data, etc. These surveys can be divided into several stages: preparation, detailed survey and analysis/identification. At the preparation stage the available information are collected and damaged structure areas are identified, which are photographed and shot on camera, the nature and scope of the damage is described. Design, construction and repair, operation documents are examined as well as the results of interviewing of persons involved in the above mentioned processes. An investigation program is prepared, in which requirement and scope of expertise work is studied. At the detailed investigation stage structure defects and their evolution are ascertained; the environment impact is characterized; form and dimensions of the structure are established; materials and their physical mechanical properties are identified; fixation conditions and loads affecting in standard and emergency loading cases are established; identification of nature, size and causes of damage and defects is performed; deflections, deformations, spectral response characteristics under operating and experimental loads are measured. At the analysis/identification stage a substantial description of the structural safety in short-term and long-term operation period is made [11-26].

Reliability of the structural health prediction depends critically on results of all stages of structural health monitoring program. Conclusions on the structural health can accurately match the reality, when damage, deformations and their causes are measured correctly. Also, the, mathematical models should well correspond to the real structure and should be properly applied to the provided lifetime of the structure.

2. Numerical models of structures

Parallel analysis algorithms and methods allow quick processing of a large amount of data and apply favorable conditions to expand nondestructive diagnostics of structures and evaluation of their condition. Those are successfully applied in complex transport and civil engineering structures. However due to approximate nature of numerical methods and uncertainty of data, the estimation of results should be carefully.

Causes of discrepancies between experimental results and numerical estimations of theoretical model can be different, among which more significant are the following:

--model structural errors, which can occur due to difficulties in specifying inhibition, connections, welding seams, edges, etc.,

--model algorithm errors, which can occur due to difficulties in specifying geometrical and material nonlinearities, etc.,

--model parameter errors, which can occur due to difficulties in specifying material and load properties, nature, etc.,

--measurement methodology, instrumental and operator errors.

These and other aspects of compliance of numerical analysis and experimental research should be taken into account, and all accepted assumptions should be motivated and balanced. It is not a simple, yet a very important stage of the structural damage identification, during which a set of calibration procedures is performed. Accuracy of models can be done by direct and inverse methods [19-25].

By direct solving a response to changes of initial parameters--geometry, material properties, supports and loads--is received. Due to these reasons stiffness and inertia properties of structure are changed, which results in changing of the nature of deformations and spectrum of dynamic response.

When solving the inverse task the experimental data are tried to bring together with the results of theoretical calculations by changing parameters of the structure. This solution is made by iterations, and structural areas and elements that are damaged (bear altered properties) are found. In this way the experimental measurements and analytical results can be identified.

The solution of inverse task requires significantly more efforts, whereas uncertainty of model parameters, just like errors of experimental measurements, has critical impact on evaluation results. Prediction of structural damages and their locations is performed by solving an inverse task according to a selected conformity criterion. It means that one should specify such a set of parameters, in the presence of which the selected criterion obtains the required value. The norm of resultant load vector can be used for such a criterion in static and dynamics tasks:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where M, C, K are mass, damping and stiffness matrixes, [??], [??], Q are acceleration, velocity and displacement vectors, F is external load vector.

A solution scheme for prediction of structural damages can be defined by the following stages:

--preparation of a numerical model of the structure and its calibration performed on measurement results,

--measurements of the structure with defects are made in presence of static and dynamic loads,

--probabilistic estimations of the degree and location of structural damages.

Compliance of the structure numerical model with its measurement results is achieved by selecting such a set of parameters of the numerical model, at which the compliance criterion is met. Global optimization algorithms are applied to this solution. Since numerical models of actual structures are characterized by a large number of degrees of freedom, therefore multiprocessor systems and parallel analysis algorithms are used to solve these tasks.

[FIGURE 1 OMITTED]

3. Identification of structural damages

The prepared structure damage identification methodology, algorithm and software were tested for beam and plate structures. Defects were initiated in the structures and drifts of such damaged structures were calculated, which were taken as the results obtained by measurements. Furthermore, applying the prepared methodology and considering that location of the defects is unknown, predictions of damaged locations were made. Finally the predictions were compared to distribution of the initiated defects.

Results of the identification of defects of the following structures are presented:

a)--statically loaded flat truss,

b)--dynamically loaded flat truss,

c)--statically loaded spatial frame,

d)--statically loaded floor.

The first task deals with a flat truss, rigidity of beams of which stretching-compression in the damaged beams was reduced. Rigidity of all remaining beams is equal. The truss is supported at ends and is loaded by a force concentrated in the center. Deflection measurements are presented in all points of the truss. Deformation nature of the damaged truss, the calculated probability of beam damage and the initiated damages are shown in Fig. 1, in which the initiated damages (90% reduction of stiffness) are shown in gray, whereas those predicted--in dark.

The second task dealt with the same structure as in the first task, however the force was time-dependent, as shown in Fig. 2, a. Force applied deflection point dependence on time with indicated damage is demonstrated in Fig. 2, b. The initiated damage here were spread over and distributed in elements 2, 4, 6, 8, 10, 20, 40. Estimated beam damage probability is shown in Fig. 2, c. Postulated (initiated) damages are marked in yellow, whereas those predicted--in blue.

[FIGURE 2 OMITTED]

The third task deals with a spatial frame, which geometry, fixations and loads, as well as initiated defects in elements 80-90 and 120-140 are shown in Fig. 3, a. The frame is loaded by a static force applied in the center and is firmly fixed in corners. Probability of damage of the frame elements is compared with initiated defects in Fig. 3, b.

These tasks were solved using unique developed software, applying which finite element models were formed, structure stress and deformation state analysis was performed, visualization and evaluation of results was presented. To this purpose MATLAB procedures library was used.

[FIGURE 3 OMITTED]

Actual civil engineering structures are usually mixed--they have beam, plate and block elements installed. If it is impossible to analyze the structure elements separately and becomes necessary to compose complex numerical models, then universal structure analysis systems, e.g. ALGOR, ABAQUS, ANSYS, etc. can be used, which offer large libraries of elements, materials and loads. Best results can be achieved when unique software is combined with universal systems because it significantly expands variety of analyzed tasks.

The presented below fourth task analyzes the structure in which plate finite elements are used for the ferroconcrete, and beam finite elements are used for the pillars. The structure finite element model, initiated damages and locations of distribution of measurement sensors are shown in Fig. 4,a. Shift of the Mises stresses of the damaged structure in the floor is illustrated in Fig. 4,a, and that of the structure with predicted damages--in Fig. 4,b.

[FIGURE 4 OMITTED]

According to these results it can be stated that by applying the proposed damage identification methodology, location and size of the predicted damages correspond to the initiated damages, and the degree of such correspondence depends on a chosen prediction methodology and quality of initial data (properties of the model, accuracy of measurements, methods of prediction). It should be noted that the damage search applied herein did not use initial solutions, which make the search process significantly more accurate and quicker.

4. Conclusions

It is noticed that success of the structure damage identification depends on reliable results of separate stages: measurements, structure modeling and prediction methodologies, equipment and personnel, and in order to reduce risk of errors reliable procedures of operations' checking should be provided for.

Algorithms and the unique software have been prepared designated for the solution of static and dynamic tasks of flat and spatial beam, shell and block structures using the finite element method, which can be applied to needs of structure strength prediction.

Several tasks for illustration purposes have been presented to demonstrate the possibilities of estimation of damage locations and risk probabilities applying both unique and universal software.

It has been noted that the applied methodology, algorithms and software identify well damage of beam and complex structures under stationary and variable loads and can be applied for practical purposes.

Acknowledgements

This research was funded by a grant (No. MIP-71/2010) from the Research Council of Lithuania.

Received January 31, 2011

Accepted June 10, 2011

References

[1.] Adams, R.D.; Coppendale, J. 1976. Measurement of the elastic module of structural adhesives by a resonant bar technique, Journal of Mechanical Engineering Science 18(3); 93-100.

[2.] Volkovas, V.; Dulevichus, J. 1975. Identification problems of dynamic models of typical pipe lines parts in diagnostics of the technical state of hydraulic systems, Scientific works of higher schools of Lithuania "Vibrotechnika" 24(3): 249-259 (in Russian).

[3.] Kargaudas, V.; Adamukaitis, N. 2010. Post-elastic force-displacement dependence of bent and compressed column, Mechanika 3(83): 5-9.

[4.] Cawley, P.; Adams, R.D. 1979. The location of defects in structures from measurements of natural frequencies, Journal of Strain Analysis 14: 49-57.

[5.] Volkovas, V.; Klumbys, A.; Ragulskis, K. 1982. Mathematical simulation and vibrodiagnostics of fault states in mechanical systems, "SEECO-82", Environ. Eng. Today. Proc. Pap. Sym. Soc. Environ. Eng., London, 13-15, July, vol.1. Butingford: 7-25.

[6.] Mozuras, A.; Volkovas, V. 1988. Simulation of defects of beam structure based on flexural vibrations, Vibration Engineering, HPC, 2(2): 75-86.

[7.] Petkevicius, K.; Smulkyte, K.; Margelis, D. 2010. Stochastic damage prediction of airframe structure, Transport Means--2010 : proceedings of the 14th international conference, October 21-22, 2010: 179-182.

[8.] Sheena, Z.; Unger, A.; Zalmanovich, A. 1982. Theoretical stiffness matrix correction by static test results, Israel Journal of Technology 20: 245-253.

[9.] Sanayei, M.; Scampoli, S.F. 1991. Structural element stiffness identification from static test data, Journal of Engineering Mechanics 117(5): 1021-1036.

[10.] Samofalov, M.; Slivinskas, T. 2009. Stability analysis of steel frames with variable cross-section for sports and entertainment centre, Mechanika 5(79): 5-12.

[11.] Sanayei, M.; Onipede, O. 1991. Damage assessment of structures using static test data, AIAA Journal 29(7): 1174-1179.

[12.] Banan, M.R.; Banan, M.R.; Hjelmstad, K.D. 1993. Parameter estimation of structures from static response. I: Computational aspects, Journal of Structural Engineering 120(11): 3243-3258.

[13.] Banan, M.R.; Banan, M.R.; Hjelmstad, K.D. 1993. Parameter estimation of structures from static response. II: Numerical simulation studies, Journal of Structural Engineering 120(11): 3259-3283.

[14.] Cui, F.; Yuan, W.C.; Shi, J.J. 2000. Damage detection of structures based on static response, Journal of Tongji University 281: 5-8.

[15.] Caddemi S.; Greco A. 2006. The influence of instrumental errors on the static identification of damage parameters for elastic beams, Computers & Structures 84: 1696-1708.

[16.] Wang, X.; Hu, N.; Fukunaga, H.; Yao, Z.H. 2001. Structural damage identification using static test data and changes in frequencies, Engineering Structures 23: 610-621.

[17.] Bakhtiari-Nejad, F.; Rahai, A.; Esfandiari, A. 2005. A structural damage detection method using static noisy data, Engineering Structures 27: 1784-1793.

[18.] Caddemi, S.; Morassi, A. 2007. Crack detection in elastic beams by static measurements, International Journal of Solids and Structures 44: 5301-5315.

[19.] Choi, I.L.; Lee, J.S.; Choi, E.; Cho, H.N. 2004. Development of elastic damage load theorem for damage detection in a statically determinate beam, Computers & Structures 82: 2483-2492.

[20.] Escobar, J.A.; Sosa, J.J.; Gomez, R. 2005. Structural damage detection using the transformation matrix, Computers & Structures 83; 357-368.

[21.] Rodriguez, R.; Escobar, J.A.; Gomez, R. 2009. Damage location and assessment along structural elements using damage submatrices, Engineering Structures 31: 475-486.

[22.] Ozen, G.O.; Kim, J.H. 2007. Direct identification and expansion of damping matrix for experimental-analytical hybrid modelling, Journal of Sound and Vibration 308: 348-372.

[23.] D'Ambrogio, W.; Fregolent, A. 2010. The role of interface DoFs in decoupling of substructures based on the dual domain decomposition, Mechanical Systems and Signal Processing 24: 2035-2048.

[24.] de Klerk, D.; Rixen, D.J.; Voormeeren, S.N. 2008. General framework for dynamic substructuring: History, review, and classification of techniques, AIAA Joural 46(5): 1169-1181.

[25.] Eun, H.C.; Lee, E.T.; Chung, H.S. 2004. On the static analysis of constrained structural systems, Canadian Journal of Civil Engineering 31: 1119-1122.

[26.] Kudzys, A.; Lukosevi?ien?, O. 2009. On the safety prediction of deteriorating structures, Mechanika 4(78): 5-11.

K. Petkevicius, Kaunas University of Technology, Kestucio 27, 44025 Kaunas, Lithuania, E-mail: [email protected]

V. Volkovas, Kaunas University of Technology, Kestucio 27, 44025 Kaunas, Lithuania, E-mail: [email protected]