Next Article in Journal
Inkjet-Printed Membrane for a Capacitive Acoustic Sensor: Development and Characterization Using Laser Vibrometer
Next Article in Special Issue
Development of an Automatic Testing Platform for Aviator’s Night Vision Goggle Honeycomb Defect Inspection
Previous Article in Journal
Characterizing Esophageal Cancerous Cells at Different Stages Using the Dielectrophoretic Impedance Measurement Method in a Microchip
Previous Article in Special Issue
A Novel Auto-Sorting System for Chinese Cabbage Seeds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inspection and Reconstruction of Metal-Roof Deformation under Wind Pressure Based on Bend Sensors

School of Automation Science and Electrical Engineering, Beihang University, XueYuan Road No. 37, HaiDian District, Beijing 100191, China
*
Authors to whom correspondence should be addressed.
Sensors 2017, 17(5), 1054; https://doi.org/10.3390/s17051054
Submission received: 30 December 2016 / Revised: 29 April 2017 / Accepted: 3 May 2017 / Published: 6 May 2017
(This article belongs to the Special Issue Innovative Sensing Control Scheme for Advanced Materials)

Abstract

:
Metal roof sheathings are widely employed in large-span buildings because of their light weight, high strength and corrosion resistance. However, their severe working environment may lead to deformation, leakage and wind-lift, etc. Thus, predicting these damages in advance and taking maintenance measures accordingly has become important to avoid economic losses and personal injuries. Conventionally, the health monitoring of metal roofs mainly relies on manual inspection, which unavoidably compromises the working efficiency and cannot diagnose and predict possible failures in time. Thus, we proposed a novel damage monitoring scheme implemented by laying bend sensors on vital points of metal roofs to precisely monitor the deformation in real time. A fast reconstruction model based on improved Levy-type solution is established to estimate the overall deflection distribution from the measured data. A standing seam metal roof under wind pressure is modeled as an elastic thin plate with a uniform load and symmetrical boundaries. The superposition method and Levy solution are adopted to obtain the analytical model that can converge quickly through simplifying an infinite series. The truncation error of this model is further analyzed. Simulation and experiments are carried out. They show that the proposed model is in reasonable agreement with the experimental results.

1. Introduction

Metal roof sheathings are generally assembled from cold-formed steel with insulation, sound-absorbing insulation and waterproof materials. Due to their advantages of light weight, high strength, flexible design, unique shape and installation convenience, they are widely employed in large-span steel-structure buildings such as exhibition buildings, performing arts centers, transportation hubs, sports stadiums and production plants [1,2]. Figure 1 presents some typical metal roof systems on large span buildings. As an example, the maximum span of the Guangzhou International Convention and Exhibition Center is approximately 132.8 m and the covered area is about 500,000 m2.
The main materials of metal roof systems include aluminum-magnesium-manganese (Al-Mn-Mg) alloy plate, titanium zinc plate and stainless steel with high strength and corrosion resistance. However, in their complex natural environment metal roof sheathings suffer the effects of wind, rain, snow, wind loads and thermal changes, leading to structural deformation, material fatigue, loose bolts and other possible failures. With the attenuation of wind-uplift resistance, the metal roof may be damaged even if the wind pressure does not exceed the design load [3]. In recent years, with the wide applications of metal roofs, several accidents were reported. For example, in March 2013, the metal roof of Beijing International Airport T3 was uplifted by strong winds, and the damaged area was about 200 square meters [4]. Figure 2 shows the large-scale broken area on the metal roof. On 18 July 2014, affected by the typhoon “Vimason”, the metal roof for the stadium of Qiongtai Normal University in HaiKou city (HaiNan province), suffered serious damage [5].
There is no doubt that the sustainability and resilience of large-scale civil infrastructure are of the utmost importance as they are closely related to people’s daily lives and social order. In recent years, some research progresses have been made on construction automation in civil engineering, intelligent sensing, and structural monitoring and health management [6,7,8,9,10]. La and Lim presented a mechatronic system design for an autonomous robotic system to inspect and evaluate bridge deck [7]. A project “Wi-Health” sponsored by the EC investigated the structural health monitoring based on multi-purpose wireless sensors network in which passive acoustic emission (AE) technology was combined with active long range ultrasonic (LRU) technology to monitor bridge structures [8].
At present, the safety detection of metal roof panels mainly relies on manual inspection, which unavoidably compromises the working efficiency and cannot diagnose and predict possible failures in time. Therefore, a real-time monitoring system is desirable to predict and estimate the potential serious deformation or damages so as to provide on-condition maintenance. When deformation or damages occur at the beginning, the deflection status of metal roof plates may change. Thus, the deflection distribution is a key characteristic and the essential issue of constructing monitoring system is to choose proper sensors and layout scheme.
When choosing detection technologies, several problems should be considered. For example, the metal roof system covers a large area, the working environment is complex and the installation space for measuring sensors is limited. Thus, it is desirable to adopt sensors that are easy to install, low-cost, and reliable under bad environmental conditions.
There have been many developments in detection approaches for curved surfaces, such as the three coordinate measuring machine, grating optical fiber [11], LASER-based methods, CT scanning methods, CCD camera detection methods [12], multimetric data (i.e., acceleration and strain) fusion [9] and so on. Although these approaches are suitable for some special situations, the characteristics of sensors such as continuous working hours, convenient installation, low power consumption and cost should be considered for actual monitoring of a metal roof. Therefore, we chose bend sensors that can be conveniently pasted on the surface of metal roof panels to detect their local deformation. Bend sensors [13] are sensitive to deformation based on optical fiber, conductive ink, electronic textiles, piezoelectric zinc oxide (ZnO) thin films [14], and tailored nanocomposite materials [15]. Bend sensors have been widely used in different fields, from simple angle measurements to the bending amount of human body joints [13,16,17,18] and shape-sensing in mechanical structures [13]. Due to their low-energy characteristics, bend sensors are suitable for a wide range of monitoring networks. Even so, it is often impractical in large-span buildings to cover all a roof surface by laying sensors in quantity. Thus, the deformation estimation from limited measurements is desired for whole roofs. Furthermore the potential fault prediction can be realized for metal roof systems.
Deformation and stress distribution analysis of metal roofs can be achieved by analytical methods and numerical simulation methods. Metal roof panels can be modeled as a kind of elastic thin plate structure and the governing equations are solved with some simplified boundary conditions. It approach has distinct physical significance but low accuracy due to the structural or boundary simplifications, which is generally used for the preliminary design and analysis. The primary analytical methods include the Navier dual-trigonometric series, Levy solution, symplectic geometry, non-linear solutions and so on. At present, more valuable studies about analytical solutions of thin plate deformation have been developed considering different load types, structure characteristics and supporting modes. Nakai [19] presented an analytical method based entirely on elasto-plastic and large deflection theory to estimate the “critical strength” of various types of thin-walled steel frames. Yang [20] used the symplectic geometry method in a Hamilton system to derive the theoretical solution for the elastic cantilever rectangular thin plate. Batista [21] proposed a new analytical solution in the form of a Fourier series for the bending problem of a corner-supported rectangular thin plate under uniformly distributed load, in which σ-method was used to speed up the convergence rates. Zhong and Vinesh [22] adopted nonlinear large deformation theory to solve the bending of elastic thin plates. This method was valid to analyze nonlinear plates and could be extended to other boundary conditions. Van Gorder [23] discussed the linearization and construction of perturbation solutions for the Foppl-von Kármán equations, a set of non-linear partial differential equations describing the large deflections of thin flat plates. It could simplify the process of inverting the linear operators greatly so as to find the higher order terms comparing to numerical methods. Segovia [24] presented an approximated analytical Levy-type solution to analyze the vibration modal of a thin folded elastic structure. Apparently, these approaches are mostly obtained by extending or improving the abovementioned classical methods and suitable for some particular sorts of problems. Although sometimes the analytical method can get only an approximate solution for the elastic thin plate under specific loads and boundary conditions, it has faster convergence rate than numerical methods. Thus it is a rational choice in view of real-time estimation and on-line analysis needs.
When complex boundary conditions and loads have to be considered, numerical simulations such as finite-difference method and finite-element analysis (FEA) are generally suggested to obtain a detailed 3D deflection and stress field. For example, Song [25] explored the deflection and stress characteristics of standing seam metal roof panel under different negative pressure load by FEA software. Wang [26] applied a cube finite element model to analyze the piezothermoelasticity of an intelligent thin plate. Zhang [27] introduced numerical manifold method (NMM) to solve thin plate bending deformation problems. Turevsky [28] introduced an efficient numerical scheme to compute the topological sensitivity (TS) of arbitrary-shaped features in plate bending. Although numerical models can describe precisely boundary conditions and plates’ shapes, such as supporting modes, joints and reinforcing ribs, they are limited to off-line system analysis, design and optimization due to high computation consumption.
In this paper, we focus on real-time reconstruction of the deflection distribution for metal roof plates from a limited number of measured data of bend sensors, as the basis of potential fault prediction. The standing seam metal roof under wind pressure is modeled as a rectangular thin elastic plate. A fast analytic method based on Levy-type solution and the infinite series simplification is presented and then the model errors (the error between deflection of metal roof computed by reconstruction model and exact value) are analyzed. Finally, the corresponding simulations and experiments are conducted and measured data are compared with analytical solutions on the given detection points. The results show that the proposed reconstruction model can estimate the deflection distribution of metal plates in real-time.

2. Fast Reconstruction Model for Deflection Distribution of Metal Roof

2.1. A Fast Reconstruction Model for the Deflection Distribution of Metal Roof

To establish the analytical model of metal roof deformation, a common standing seam metal roof is taken as an example. The entire metal roof is made up of many single metal roof sheathings locked together. As shown in Figure 3, two pieces of roof sheathing are assembled on the support. As the effect of the standing edges and buckling edges (green part in Figure 3), the two sides of y direction which form occlusal constraints with the metal holders can be seen as simply supported edges. The two sides in the x direction are equivalent to free edges because of the absence of constraints.
Under the same weather condition, a single block of metal roof suffers a wind uplift equivalent to a uniform load. Therefore, the entire metal roof can be seen as a rectangular plate with two opposed simply supported edges and two other free opposed-edges. The control equations can be solved by external constraints and applied loads.
To analyze the deflection of the metal roof based on thin plate bending theory, the Kirchhoff hypothesis [29] is adopted as follows: (1) the thickness of the plate is much smaller than the horizontal and vertical dimensions; the density of the structure is constant throughout; (2) before the thin plate deforms, a straight line is perpendicular to the middle plane. It is still perpendicular to the middle plane after the thin plate deforms. The points on the middle surface have no horizontal displacement.

2.1.1. Governing Equations for Deflection Deformation of Thin Plate

As shown in Figure 4, we took a hexahedral element in the thin plate whose thickness is d. Assume that the upper surface is subjected to a uniform load q. Then its stress condition is analyzed.
From assumption (2), points in the middle have no displacement in the x and y direction, u and v are non-existent, and the displacements u and v whose distance were z from the middle could be obtained as:
u = z w x , v = z w y .
According to the Hooker theorem, we can get the normal stress and shear stress as follows:
σ x = E z 1 u 2 [ 2 w x 2 + u 2 w y 2 ] , σ y = E z 1 u 2 [ 2 w y 2 + u 2 w x 2 ] , τ x y = 2 G z 2 w x y and G = E 2 ( 1 + u ) .
The stress components τ z x , τ z y are determined by the following two equations:
σ x x + τ y z y + τ z x z = 0 , σ y y + τ z y z + τ x y x = 0 .
As shown in Figure 4, the normal stress and shear stress are linearly distributed in the z direction. The corresponding bending moments M x , M y , M x y and shearing forces Q x , Q y can be obtained as:
M x = d 2 d 2 σ x z d z = D ( 2 w x 2 + ν 2 w y 2 )
M y = d 2 d 2 σ y z d z = D ( 2 w y 2 + ν 2 w x 2 )
M x y = d 2 d 2 τ x y z d z = D ( 1 ν ) 2 w x y
Q x = d 2 d 2 τ x z d z = D ( 3 w x 3 + ( 2 μ ) 3 w x y 2 )
Q y = d 2 d 2 τ y z d z = D ( 3 w y 3 + ( 2 μ ) 3 w y x 2 )
The governing equations of the bending plate can be obtained based on the force balance as follows [30]:
4 w x 4 + 2 4 w x 2 y 2 + 4 w y 4 = q D ,
where D = E d 3 / 12 ( 1 μ 2 ) is the flexural rigidity of the plate, E is Young’s modulus of the plate, µ is Poisson’s ratio, w is the vertical deflection and q is the external uniform load.
As shown in Figure 3, the metal roof can be simplified as a thin plate with two opposite simply-supported edges, two free edges and subjected to a uniformly distributed load.
The boundary conditions are as follows:
I
Free edges (y = 0, y = l): the bending moment and shear force are zero, and formulated as:
D ( 2 w x 2 + μ 2 w y 2 ) = 0 , D ( 3 w x 3 + ( 2 μ ) 3 w x y 2 ) = 0 .
II
Supported Edges (x = 0, x = b): the deflection and bending moment are zero:
D ( 2 w x 2 + μ 2 w y 2 ) = 0 .
The superposition method is used to solve the above problems. They are treated as the problem that a rectangular thin plate with the two opposite free sides and the other two opposite sides simply supported subjected to a uniform load q. The two parts of superposition are taken as follows:
(a)
A thin rectangular plate with four sides simply supported, and subjected to uniform load q.
(b)
A rectangular plate with two sides x = 0 and x = b seen as simply supported sides, and the other sides y = 0 and y = l seen as generalized simply supported edges.

2.1.2. Deflection Calculation of Part I

The deflection of this part is assumed as w a ( x , y ) , and the expression of its double triangle series is:
w a ( x , y ) = m = 1 n = 1 α m n sin m π x b sin n π y l ,
The double triangle series expansion of the uniform load q is:
q = m = 1 n = 1 A m n sin m π x b sin n π y l ,
where A m n = 16 q π 2 1 m n (m, n = 1, 3, 5…), and A m n = 0 (m and n is even).
Equations (7) and (8) are substituted into Equation (6), then both sides of the equation coefficients are compared and the final deflection equation is:
w a ( x , y ) = 16 q π 6 D m = 2 i + 1 i = 0 1 m 5 sin m π x b n = 2 j + 1 j = 0 l 4 n ( n 2 / m 2 + l 2 / b 2 ) 2 sin n π y l .
Equation (9) is double trigonometric series, so to speed up the convergence rate, the sum of series containing y is calculated, and the Levy-type solution [31] is adopted and Equation (9) can be rewritten as:
w a ( x , y ) = 4 q b 4 π 5 D m = 2 j + 1 j = 0 1 m 5 g m ( x ) f m ( y ) ,
where:
g m ( x ) = sin m π x b ,
f m ( y ) = [ λ m ( y l 2 ) sinh λ m ( y l 2 ) ( 2 + λ m l 2 tanh λ m l 2 ) cosh λ m ( y l 2 ) ] 1 2 cosh ( λ m l / 2 ) + 1 ,
and λ m = m π / b (m = 1, 3, 5…).
By investigating the convergence rate of the function series, we find that the convergence of the function series depends mainly on the first item. Take the first term of series, i.e.:
w a ( 1 ) ( x , y ) = 4 q b 4 π 5 D g 1 ( x ) f 1 ( y ) ,
and the remainder is:
R a = 4 q b 4 π 5 D m = 2 i + 1 i = 1 1 m 5 g m ( x ) f m ( y ) .
Refer to Appendix A. It can be proved that:
| R a w a ( 1 ) ( x , y ) | < < 1 ,
thus, it is reasonable to take the first term of series to approximate w a so as to decrease computing time. Thus we get:
w ˜ a ( x , y ) = w a ( 1 ) ( x , y ) = 4 q b 4 π 5 D g 1 ( x ) f 1 ( y ) .

2.1.3. Deflection Calculation of Part II

The deflection of part II is w b .Using the same methods as mentioned in Section 2.1.1 to analyze the deflection of part II, the following result can be obtained:
w b ( x , y ) = q b 4 π 5 D n = 2 j + 1 j = 0 1 n 5 g n ( x ) h n ( y )
where λ n = n π / b and
h n ( y ) = ( a 1 + c 1 λ n y ) sinh λ n y + ( a 2 + c 2 λ n y ) cosh λ n y .
To obtain h n ( y ) , boundary conditions are used to solve the problem:
(1)
Moment of the two generalized simply supported edges is zero, i.e., M = 0 | y = 0 , l
(2)
Deflection of the two generalized simply supported sides is assumed, i.e., w = α 0 sin ( β x ) | y = 0 , l , and h n ( y ) can be given:
h n ( y ) = α 0 { cosh λ n l 1 sinh λ n l [ ( λ n l sinh λ n l 1 1 μ ) sinh λ n y + λ n y cosh λ n y ] + 1 1 μ cosh λ n y λ n y sinh λ n y }
According to the boundary conditions, α 0 is calculated:
α 0 = 2 tanh λ n l 2 [ ( 4 μ ) ( 1 μ ) λ n l 2 cosh ( λ n l / 2 ) ] ( 1 μ ) 2 ( cosh λ n l 1 ) ( 3 + μ 1 μ λ n l sinh λ n l ) .
The first term of Equation (12) is represented as:
w b ( 1 ) ( x , y ) = q b 4 π 5 D h 1 ( y ) g 1 ( x ) ,
and the remainder is:
R b = q b 4 π 5 D n = 2 j + 1 j = 1 1 n 5 g n ( x ) h n ( y ) .
It can be proved that (Appendix A):
| R b w b ( 1 ) ( x , y ) | < < 1 .
Then the estimated formula is obtained as:
w ˜ b ( x , y ) = w b ( 1 ) ( x , y ) = q b 4 π 5 D h 1 ( y ) g 1 ( x ) .
Superimposing the deflection of two parts gives the deflection of the coordinate ( x , y ) on metal roof:
w ( x , y ) = w a ( x , y ) + w b ( x , y ) = q b 4 π 5 D m = 2 i + 1 i = 0 1 m 5 g m ( x ) [ 4 f m ( y ) + h m ( y ) ] .
Take the first term of series as the estimated deflection of the metal roof, i.e.:
w ˜ ( x , y ) = w ( 1 ) ( x , y ) = q g 1 ( x ) b 4 π 5 D ( 4 f 1 ( y ) + h 1 ( y ) ) .
When the uniform load q is unknown, the deflection of installation location ( x 0 , y 0 ) can be measured by a bend sensor, and the value is w 0 ( x 0 , y 0 ) . Substituting w 0 ( x 0 , y 0 ) and ( x , y ) into Equation (26) yields the deflection of any point on the metal-roof:
w ˜ ( w 0 , x , y ) = w 0 ( x 0 , y 0 ) g 1 ( x 0 ) ( 4 f 1 ( y 0 ) + h 1 ( y 0 ) ) g ( x ) 1 ( 4 f 1 ( y ) + h 1 ( y ) ) .

2.2. Analysis of Model Error

Based on the above analysis, the exact deflection of point ( x , y ) on a metal roof is represented in Equation (25) and the estimated value is in Equation (26). The remainder is:
R = R a + R b = q b 4 π 5 D m = 2 i + 1 i = 1 1 m 5 g m ( x ) [ 4 f m ( y ) + h m ( y ) ] .
The relative error between w ( x , y ) and w ˜ ( x , y ) is:
ε = | R w ( x , y ) | .
According to Equations (11), (12), (19), (23) and (27), to calculate the value of ε , the value of x , y , l , b should be determined. We choose 65/400 straight side whipstitch plate type roof panel [32] for experiment. Its properties are shown in Table 1.
The value scope of x is 0~400 mm, the value scope of y is 0~1700 mm and l / b = 17 / 4 . According to Equations (15) and (23), we can obtain:
ε = | R a + R b w a ( x , y ) + w b ( x , y ) | | R a + R b w ˜ a ( x , y ) + w ˜ b ( x , y ) | | R a w ˜ a ( x , y ) | + | R b w ˜ b ( x , y ) | .
According to Figure A5 and Equation (A9) of Appendix B, it can be obtained that the value of | R b w b ( 1 ) ( x , y ) | is approximately 0. And according to Equation (A8), the maximum relative error is:
ε | R a w ˜ a ( x , y ) | + | R b w ˜ b ( x , y ) | | ( 7 * ζ ( 3 ) ) / 8 1 | + | π 4 48 2 | × 10 7 0.0518 .

3. Simulations and Experiments

In order to verify the proposed fast reconstruction model of deflection distribution for metal roofs, the corresponding simulations and experiments are carried out and discussed in this section. Firstly, we performed a real value simulation (using the practical size and properties of metal roof sheathings). We simulated the deflection distributions of a metal roof panel under uniform load by using the proposed analytical method and the FEA method. The results and computation time are analyzed to validate the accuracy and rapidity of the fast reconstruction model. Next, an experimental standing seam metal roof platform is set up according to the metal roofs used in practical engineering applications and the bend sensors are calibrated. Finally, a series of experiments are designed and implemented to verify the practicability of online monitoring scheme of metal roof deformation based on bend sensors and the fast reconstruction model.

3.1. Simulation and Discussions

In the real-size simulation experiments, the 65/400 straight side whipstitch plate type roof panel is selected as test sample and its properties are shown in Table 1. The mathematical software MATLAB 2012 is used to compute the reconstruction model of this metal roof panel in which 20 × 85 points ( x × y ) are picked uniformly on the whole metal roof panel for deflection calculation. Meanwhile, we use the FEA software COMSOL 5.3b (COMSOL Inc, Stockholm, Sweden) to simulate the deflection distribution of metal roof panels which were meshed into 64 thousand elements on the same computer.
We use software MATLAB and COMSOL to simulate the deformation of metal roof panel with loads from 0~3000 Pa in steps of 200 Pa. When the equivalent uniform load q imposed on metal panel is 3000 Pa [33], Figure 5 shows the simulation results of the FEA model in COMSOL and the fast reconstruction model in MATLAB. It can be seen that the deflection distributions are roughly consistent in both sets of results. Furthermore, we take the COMSOL data of picked points (five lines and eleven rows of points on the plate, i.e., black points in Figure 5a) as the reference base to examine the relative error of the analytical model. Figure 6 shows deflection under different loads in the x direction (y = 0.85 m), and Figure 7 shows deflection under different loads in y direction (x = 0.18 m, avoid the location of the stiffeners). From Figure 6 and Figure 7, it can be seen that the results of reconstruction model are consistent with the simulation results in COMSOL. The error data are shown in Figure 8. The absolute error is mm-level, and maximum relative errors occur on the stiffener locations where the stiffeners are modeled in COMSOL simulation but simplified through assumptions in the analytical model. We are more interested in the maximum deflection and the deflection at the boundaries. The maximum error on these key positions is about 2 mm and percentage error is 2% where our applicable percentage error is 5%, so the fast analytical model is applicable to deflection calculations of metal roof panel.
In the same hardware platform, the calculation time of MATLAB is about 20 ms, whereas that of COMSOL is approximately 30 min. Obviously, the fast reconstruction model saves computer hardware resources and takes less time to compute the metal roof deflection with comparable precision. Moreover, according to the measured circuit design features and sampling period we have known that the general time of collecting and transmitting data is about 15 ms, thus the fast reconstruction model based on measured data can meet the requirements of real-time monitoring.

3.2. Experimental Platform

In order to facilitate the experiment, a piece of Al-Mg-Mn metal roof panel is incised to set up experimental platform. The experimental platform is structured with 65/400 type (Al-Mg-Mn alloy) standing seam plate, bracket, bend sensor, ultrasonic module, Zigbee module and power module, as shown in Figure 9. The reference point of the bend sensor is M and its coordinate is ( 250 , 70 ) (mm). The given measurement points of the ultrasonic module are A ( 200 , 75 ) , B ( 150 , 75 ) , C ( 350 , 50 )   ( mm ) on the metal plate. The ultrasonic module is fixed on the bracket, whose upside is against the given points. The measured deflections of all sensors are transmitted by Zigbee modules, to PC for further analysis. The material properties of standing seam metal roof are given in Table 2.
We adopted the FS-L-0112-103-ST flex sensor (Spectra System, Salt Lake, UT, USA) in the experiment [13]. It is cheap, easy to integrate and requires no complicated signal processing. According to the manual, its continuous power rating is 0.50 W; its life cycle is more than one million uses; the sensitivity can be adjusted by circuit design [34].
The resistance changes with the bending status. Additionally, the working temperature range of the sensor is wide and not affected by the humidity. It can adapt to the complex and changeable situation of the temperature and humidity on the top of buildings. As shown in Figure 9, the bend sensor is stuck on the underside of the metal-roof plate. The bending status of the sensor changes with any deformation of the metal roof, which leads to a corresponding resistance change.
The bend sensor is calibrated on the experimental platform firstly. In the calibration experiments, we apply uniformly negative loads on the metal roof to make the deflection of the measuring point vary from 0 to 90 mm. Every deflection state lasts 60 s and the resistances of the flex sensor are measured with 1 Hz sampling frequency, but only stable values are recorded as calibration data after filtering. Then we use the data to fit deflection-resistance curve as shown in Figure 10 and the relationship of deflection-resistance value is assumed as:
w ( r ) = a 0 × r 5 + a 1 × r 4 + a 2 × r 3 + a 3 × r 2 + a 4 × r .
The coefficients are given as:
a 0 = 0.02125 , a 1 = 1.091 , a 2 = 20.83 , a 3 = 176.3 , a 4 = 557.2
When the deflection measured by bend sensor is 50 mm, 40 × 30 points ( x × y ) on the experimental metal roof are taken to calculate their deflections according to Equation (22), and the deflection distribution of entire metal roof is plotted as shown in Figure 11.

3.3. Experimental Results and Discussion

Imposed stepped reverse loads on the metal roof make the deflection of point M measured by the bend sensor change from 0 to 60 mm in 10 mm steps, approximately. Simultaneously, we report the deflections of points A, B and C from ultrasonic sensors. Substituting the measured values of the M point into Equation (22) we can get the theoretical deflection values of points A, B and C. The measured and estimated values of points A, B and C are summarized in Table 3.
The relationships between theoretical and measured values of points A, B and C are shown as Figure 12, where w t A , w t B and w t C represent the theoretical value of points A, B and C; w m A , w m B and w m C represent the measured values at points A, B and C.
From Table 3 and Figure 12, the estimated deflections calculated by the reconstruction model agree well with the measured data on the whole. The maximum relative error (7.1%) occurs at point A which is on the metal roof stiffener; this is because the stiffener can increase the strength of the plate and reduce the deformation of the metal roof. Additionally, in lines A and B, there is a slight difference between the measured data and theoretical calculation at higher deflection values. The possible reason is that the bend deformation of the metal roof panel may be out of the elastic deformation range at higher deflection values in the actual experiments. Meanwhile, the bend deformation is assumed elastic in the whole deformation process in the reconstruction model. Its effect can be neglected since the difference is small and not significant for damage estimation.

4. Conclusions

A novel deformation-monitoring scheme for metal roof sheathing based on bend sensors has been presented. A fast reconstruction model is established to estimate the deflection distribution of whole roof panels based on limit measured data. The superposition method and Levy solution are adopted to obtain an analytical model with quick convergence by simplifying infinite series. Then, the simulation study of the real-size metal roof panel shows that the results of the proposed method agree well with those of the COMSOL model but require much less computation time. Finally, an experimental investigation is conducted to validate this scheme. The results confirm that the fast reconstruction model can estimate the deflection distribution of the whole metal plate in terms of computational accuracy and time, which provides a data basis for the online health monitoring of large-span metal roof sheathing systems which will be further studied in our future work.

Acknowledgments

This work has been supported in part by the State Key Laboratory of Fluid Power & Mechatronic System Foundation under Grant No. GZKF-201511 and by the Chinese Aerospace Science and Technology Creative Foundation under Grant No. CASC150706.

Author Contributions

L.Y. and Y.L. conceived the algorithm and designed the experiments; L.C. performed the experiments and analyzed the data; C.A. completed the simulations; L.Y. and L.C. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

1. According to Equations (8) and (9), we can get that:
| R a w a ( 1 ) ( x , y ) | = | m = 2 i + 1 i = 1 1 m 5 g m ( x ) f m ( y ) g 1 ( x ) f 1 ( y ) | .
According to Equations (6) and (7), to calculate the value of | R a w a ( 1 ) ( x , y ) | , the value of x , y , b , l should be given. Obviously, the value range of x is from 0 to b, and the value range of y is from 0 to l. In actual use, the width of standing seam metal roof panel is 300 mm~500 mm, and the max length is up to 15,000 mm. Therefore, 1 < l / b < 50 and here we consider extreme situation in order to compute maximum truncation error, i.e., l / b = 1 . Define G m ( x ) = 1 m | g m ( x ) g 1 ( x ) | and F m ( y ) = 1 2 m | f m ( y ) f 1 ( y ) | . The value curves of G m ( x ) for m = 3, 5, 7, 9 are shown as Figure A1.
Figure A1. The value of G m ( x ) .
Figure A1. The value of G m ( x ) .
Sensors 17 01054 g013
From Figure A1, we can infer that:
G m ( x ) = 1 m | g m ( x ) g 1 ( x ) | = 1 m | sin m π x b sin π x b | 1 ( m = 2 i + 1   and   i = 1 , 2 , ) .
The value of F m ( y ) is shown in Figure A2.
Figure A2. The value of F m ( y ) .
Figure A2. The value of F m ( y ) .
Sensors 17 01054 g014
and it can be obtained that:
F m ( y ) = 1 2 m | f m ( y ) f 1 ( y ) | 0.75 .
So that:
| R a w a ( 1 ) ( x , y ) | = | m = 2 i + 1 i = 1 1 m 5 g m ( x ) f m ( y ) g 1 ( x ) f 1 ( y ) | = | m = 2 i + 1 i = 1 2 m 3 G m ( x ) F m ( y ) | < | m = 2 i + 1 i = 1 1.5 m 3 | = | 21 × ζ ( 3 ) 16 3 2 | 0 . 0777 < < 1 ,
where ζ ( 3 ) is three-order Riemann function.
2. According to Equation (14) and defining H n ( y ) = | h n ( y ) h 1 ( y ) | , Figure A3 shows the value of H n ( y ) with different n, it can be seen that the value of H n ( y ) tends to zero rapidly as n increases. Thus:
H n ( y ) = | h n ( y ) h 1 ( y ) | < < 1 .
Figure A3. The value of H n ( y ) : (a) n = 3; (b) n = 5; (c) n = 7; (d) n = 9.
Figure A3. The value of H n ( y ) : (a) n = 3; (b) n = 5; (c) n = 7; (d) n = 9.
Sensors 17 01054 g015
According to Equations (14)–(16), it can be got that:
| R b w b ( 1 ) ( x , y ) | = | n = 2 j + 1 j = 1 1 n 5 g n ( x ) h n ( y ) g 1 ( x ) f 1 ( y ) | = | n = 2 j + 1 j = 1 1 n 4 G n ( x ) H n ( y ) | | n = 2 j + 1 j = 1 1 n 4 H n ( y ) | = | π 4 96 1 | | H n ( y ) | < < 1

Appendix B

We chose a 65/400 straight side whipstitch plate type roof panel as analysis object. Its properties are listed in Table 1. The value range of x is from 0 to 0.4 m, the value range of y is from 0 to 1.7 m, and l / b = 17 / 4 . According to Equation (A1), it can be obtained when G m ( x ) 1 .
The value of F m ( y ) is shown as Figure A4.
Figure A4. The value of F m ( y ) .
Figure A4. The value of F m ( y ) .
Sensors 17 01054 g016
and it can be obtained that:
F m ( y ) = 1 2 m | f m ( y ) f 1 ( y ) | 0.5 ,
so that:
| R a w a ( 1 ) ( x , y ) | = | m = 2 i + 1 i = 1 1 m 5 g m ( x ) f m ( y ) g 1 ( x ) f 1 ( y ) | = | m = 2 i + 1 i = 1 2 m 3 G m ( x ) F m ( y ) | < | m = 2 i + 1 i = 1 1 m 3 | = | 7 * ζ ( 3 ) 8 1 | 0 . 0518 .
The value of H n ( y ) is shown as Figure A5.
Figure A5. The value of H n ( y ) : (a) n = 3; (b) n = 5; (c) n = 7; (d) n = 9.
Figure A5. The value of H n ( y ) : (a) n = 3; (b) n = 5; (c) n = 7; (d) n = 9.
Sensors 17 01054 g017aSensors 17 01054 g017b
According to Equations (A2) and (A3), and Figure A5, it can be obtained that:
| R b w b ( 1 ) ( x , y ) | = | n = 2 j + 1 j = 1 1 n 5 g n ( x ) h n ( y ) g 1 ( x ) f 1 ( y ) | = | n = 2 j + 1 j = 1 1 n 4 G n ( x ) H n ( y ) | | n = 2 j + 1 j = 1 1 n 4 H n ( y ) | | π 4 48 2 | × 10 7 .

References

  1. Huang, W.; Dong, W.M. Introduction for metal roofing development and application classification. Archit. Technol. 2015, 12, 1103–1107. [Google Scholar]
  2. Luo, Y.F.; Zheng, X.J.; Guo, X.N.; Li, D.Z.; Hu, Z. Loading behavior analysis of an aluminum ally roof structure system under wind loads. J. Chongqing Univ. 2013, 36, 94–100. [Google Scholar]
  3. Baskaran, A.; Molleti, S.; Ko, S.; Shoemaker, L. Wind uplift performance of composite metal roof assemblies. J. Archit. Eng. 2012, 18, 2–15. [Google Scholar] [CrossRef]
  4. Long, W.Z. Discussion on improving metal roof anti wind technology The roof of Beijing Capital International Airport T3 do not be lifted by the wind fourth times again. China Constr. Metal Struct. 2013, 13, 62–68. [Google Scholar]
  5. Ying, X.J.; Zhong, J.H. Analysis on damage cause of the metal roof for the stadium of qiongtai normal college & its reinforcement. China Build. Waterproofing 2015, 11, 14–17. [Google Scholar]
  6. Li, Y.H.; Wang, Y.; Chase, J.G.; Mattila, J.; Myung, H.; Sawodny, O. Survey and introduction to the focused section on mechatronics for sustainable and resilient civil infrastructure. IEEE/ASME Trans. Mechatron. 2013, 18, 1637–1646. [Google Scholar] [CrossRef]
  7. La, H.M.; Lim, R.S.; Basily, B.B.; Gucunski, N.; Yi, J.G.; Maher, A.; Romero, F.A.; Parvardeh, H. Mechatronic systems design for an autonomous robotic system for high-efficiency bridge deck inspection and evaluation. IEEE/ASME Trans. Mechatron. 2013, 18, 1655–1664. [Google Scholar] [CrossRef]
  8. Graham, R.E.; Kamer, T.; Bob, W. Development of a multi-purpose wireless network for the structural health monitoring of a suspension bridge. In Proceedings of the IET Conference on Wireless Sensor Systems (WSS 2012), London, UK, 18–19 June 2012. [Google Scholar]
  9. Park, J.W.; Sim, S.H.; Jung, H.J. Displacement estimation using multimetric data fusion. IEEE/ASME Trans. Mechatron. 2013, 18, 1675–1682. [Google Scholar] [CrossRef]
  10. Lin, C.J.; Yau, H.T.; Lee, C.Y.; Tung, K.H. System identification and semiactive control of a squeeze-mode magnetorheological damper. IEEE/ASME Trans. Mechatron. 2013, 18, 1691–1701. [Google Scholar] [CrossRef]
  11. Fan, H.C.; Zhang, Z.; Zhu, X.J. Application of FBG in flexible surface test. Mech. Electr. Eng. Mag. 2009, 26, 40–43. [Google Scholar]
  12. Chan, T.H.T.; Ashebo, D.B.; Tam, H.Y.; Yu, Y.; Chan, T.F.; Lee, P.C.; Gracia, E.P. Vertical displacement measurements for bridges using optical fiber sensors and CCD cameras. Struct. Health Monit. 2009, 8, 243–249. [Google Scholar] [CrossRef]
  13. Jeong, U.; Cho, K.J. A novel low-cost, large curvature bend sensor based on a Bowden-cable. Sensors 2016, 16, 961. [Google Scholar] [CrossRef] [PubMed]
  14. Marco, L.; Denis, P.; Alessio, V.; Candido, F.P.; Alessandro, C. Development of a Flexible Lead-Free Piezoelectric Transducer for Health Monitoring in the Space Environment. Micromachines 2015, 6, 1729–1744. [Google Scholar]
  15. Chiolerio, A.; Roppolo, I.; Sangermano, M. Radical diffusion engineering: Tailored nanocomposite materials for piezoresistive inkjet printed strain measurement. RSC Adv. 2013, 3, 3446–3452. [Google Scholar] [CrossRef]
  16. Gentner, R.; Classen, J. Development and evaluation of a low-cost sensor glove for assessment of human finger movements in neurophysiological settings. J. Neurosci. Methods 2009, 178, 138–147. [Google Scholar] [CrossRef] [PubMed]
  17. Tognetti, A.; Lorussi, F.; Carbonaro, N.; de Rossi, D. Wearable Goniometer and Accelerometer Sensory Fusion for Knee Joint Angle Measurement in Daily Life. Sensors 2015, 15, 28435–28455. [Google Scholar] [CrossRef] [PubMed]
  18. Stoppa, M.; Chiolerio, A. Wearable electronics and smart textiles: A critical review. Sensors 2014, 14, 11957–11992. [Google Scholar] [CrossRef] [PubMed]
  19. Nakai, H.; Miki, T.; Ohgaki, K. Analytical method for critical strength of thin-walled steel frames. Mem. Fac. Eng. Osaka City Univ. 1985, 26, 233–250. [Google Scholar]
  20. Zhong, Y.; Chen, J.Y.; Wang, S.Y. Analytical solution for rectangular thin cantilever plate. Chin. J. Comput. Mech. 2006, 23, 368–372. [Google Scholar]
  21. Batista, M. New analytical solution for bending problem of uniformly loaded rectangular plate supported on corner points. IES J. Part A Civ. Struct. Eng. 2010, 3, 75–84. [Google Scholar] [CrossRef]
  22. Nishawala, V.V. A Study of Large Deflection of Beams and Plates. Master’s Thesis, Rutgers State University, New Brunswick, NJ, USA, 2011. [Google Scholar]
  23. Van, G.; Robert, A. Analytical method for the construction of solutions to the Föppl-von Kármán equations governing deflections of a thin flat plate. Int. J. Non-Linear Mech. 2012, 47, 1–6. [Google Scholar]
  24. Segovia, E.; Ferrer, B.; Ramis, J.; Martinez, J.; Arenas, J.P. Vibration Modal Analysis of a Thin Folded Elastic Structure using a Levy-type Solution. Int. J. Acoust. Vib. 2012, 17, 191–199. [Google Scholar] [CrossRef]
  25. Song, Y.H.; Yang, L.M.; Wang, Q.S.; Li, Y.X. Deformation characteristic study of metal roof panel under wind uplift loading. Build. Struct. 2015, 45, 87–91. [Google Scholar]
  26. Wang, X.B.; Chen, J.J.; Liang, Z.T. Randomicity analysis of the piezothermoelasticity intelligent thin plate by appling the numerical method. J. Xidian Univ. 2008, 35, 592–599. [Google Scholar]
  27. Zhang, Z.R.; Zhang, X.W.; Lü, W.G. Numerical method based on compatible manifold element for thin plate bending. Chin. J. Mech. Eng. 2010, 23, 100–109. [Google Scholar] [CrossRef]
  28. Turevsky, I.; Sankara, H.G.; Suresh, K. An efficient numerical method for computing the topological sensitivity of arbitrary shaped features in plate bending. Int. J. Numer. Methods Eng. 2009, 79, 1683–1702. [Google Scholar] [CrossRef]
  29. He, X.M.; Zou, W.M. Infinitely many positive solutions for Kirchhoff-type problems. Nonlinear Anal. 2009, 70, 1407–1414. [Google Scholar] [CrossRef]
  30. Banerjee, A.; Bhattacherya, B.; Mallika, K. Large deflection of cantilever beams with geometric non-linearity: Analytical and numerical approaches. Int. J. Non-Linear Mech. 2008, 5, 366–376. [Google Scholar] [CrossRef]
  31. Simmonds, J.G. Exact Levy-type solutions for plate bending exist for transversely isotropic but not for general monoclinic materials. J. Elast. 2004, 75, 49–56. [Google Scholar] [CrossRef]
  32. GB50896-2013 M.O. Technical Code for Application of Profiled Metal Sheets; Ministry of Housing and Urban-Rural Construction of the People’s Republic of China, and General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China: Beijing, China, 2014.
  33. GB50009 M.O. Load Code for the Design of Building Structures; Ministry of Housing and Urban-Rural Construction of the People’s Republic of China, and General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China: Beijing, China, 2012.
  34. Flex Sensor Manufacturers Spectra Symbol. Available online: http://www.spectrasymbol.com/flex-sensor (accessed on 12 April 2016).
Figure 1. The application of metal roof in buildings: (a) Olympic Sports Center of Guiyang; (b) Beijing-Capital International Airport; (c) The production plant of Hubei Wanmeng CNC combination machine Limited company; (d) Guangzhou International Convention and Exhibition Center.
Figure 1. The application of metal roof in buildings: (a) Olympic Sports Center of Guiyang; (b) Beijing-Capital International Airport; (c) The production plant of Hubei Wanmeng CNC combination machine Limited company; (d) Guangzhou International Convention and Exhibition Center.
Sensors 17 01054 g001
Figure 2. Wind-uplift destruction of metal roofs; (a) Beijing International Airport T3; (b) The metal roof panels of Beijing International Airport T3 were split by wind.
Figure 2. Wind-uplift destruction of metal roofs; (a) Beijing International Airport T3; (b) The metal roof panels of Beijing International Airport T3 were split by wind.
Sensors 17 01054 g002
Figure 3. Metal roof model and its simplified diagram.
Figure 3. Metal roof model and its simplified diagram.
Sensors 17 01054 g003
Figure 4. The stress distribution of hexahedral element.
Figure 4. The stress distribution of hexahedral element.
Sensors 17 01054 g004
Figure 5. Comparison of two simulation methods; (a) Simulation results in COMSOL; (b) Simulation results in MATLAB.
Figure 5. Comparison of two simulation methods; (a) Simulation results in COMSOL; (b) Simulation results in MATLAB.
Sensors 17 01054 g005
Figure 6. Deflection under different pressure in x direction (y = 0.85 m).
Figure 6. Deflection under different pressure in x direction (y = 0.85 m).
Sensors 17 01054 g006
Figure 7. Deflection under different pressure in y direction (x = 0.18 m).
Figure 7. Deflection under different pressure in y direction (x = 0.18 m).
Sensors 17 01054 g007
Figure 8. Relative error of the two methods.
Figure 8. Relative error of the two methods.
Sensors 17 01054 g008
Figure 9. Standing seam metal experimental table.
Figure 9. Standing seam metal experimental table.
Sensors 17 01054 g009
Figure 10. Flex sensor resistance-deflection curve fitting.
Figure 10. Flex sensor resistance-deflection curve fitting.
Sensors 17 01054 g010
Figure 11. Deformation of entire metal roof.
Figure 11. Deformation of entire metal roof.
Sensors 17 01054 g011
Figure 12. Theoretical and measured deflection values of different points.
Figure 12. Theoretical and measured deflection values of different points.
Sensors 17 01054 g012
Table 1. Material properties of 65/400 standing seam metal roof.
Table 1. Material properties of 65/400 standing seam metal roof.
MaterialLength l (mm) Width b (mm)Thickness d (mm)Poisson Ratio Young Modulus (Pa) E (N/m2)
AA3004 Al-Mg-Mn alloy170040010.37 × 1010
Table 2. Material properties of standing seam metal roof.
Table 2. Material properties of standing seam metal roof.
Roof MaterialLock EdgeFree EdgeThicknessPoisson RatioYoung’s Modulus
l (mm)b (mm)d (mm)μE (N/m2)
AA3004 Al-Mg-Mn alloy15040010.37 × 1010
Table 3. Comparison of measured and theoretical values of deflection of metal roof.
Table 3. Comparison of measured and theoretical values of deflection of metal roof.
Measured PointsDeflection (mm)Max Relative Error (%)
M09.622.331.540.753.559.6
ATV010.424.134.142.757.964.57.1%
MV010.222.533.042.054.662.0
BTV09.622.331.540.753.559.63.6%
MV09.321.532.742.153.857.5
CTV03.99.213.020.122.024.65.3%
TV03.88.912.519.121.523.9
Note: TV represents theoretical value; and MV represents measured value.

Share and Cite

MDPI and ACS Style

Yang, L.; Cui, L.; Li, Y.; An, C. Inspection and Reconstruction of Metal-Roof Deformation under Wind Pressure Based on Bend Sensors. Sensors 2017, 17, 1054. https://doi.org/10.3390/s17051054

AMA Style

Yang L, Cui L, Li Y, An C. Inspection and Reconstruction of Metal-Roof Deformation under Wind Pressure Based on Bend Sensors. Sensors. 2017; 17(5):1054. https://doi.org/10.3390/s17051054

Chicago/Turabian Style

Yang, Liman, Langfu Cui, Yunhua Li, and Chao An. 2017. "Inspection and Reconstruction of Metal-Roof Deformation under Wind Pressure Based on Bend Sensors" Sensors 17, no. 5: 1054. https://doi.org/10.3390/s17051054

APA Style

Yang, L., Cui, L., Li, Y., & An, C. (2017). Inspection and Reconstruction of Metal-Roof Deformation under Wind Pressure Based on Bend Sensors. Sensors, 17(5), 1054. https://doi.org/10.3390/s17051054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop