Open Access
Issue
Int. J. Simul. Multisci. Des. Optim.
Volume 5, 2014
Article Number A06
Number of page(s) 6
DOI https://doi.org/10.1051/smdo/2013009
Published online 06 February 2014

© J.L. Zapico-Valle, M. García-Diéguez, Published by EDP Sciences, 2014

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Frame steel structures are extensively used in industrial and building construction. Linear models are usually adopted in structural applications as design or health monitoring. This approach strongly simplifies the modeling and extraction of modal properties. Linearity, however, seldom occurs in reality. Structures always contain some degree of nonlinearity due to material and interface behavior, geometric and elastic instabilities, etc. An epistemic error is therefore introduced by linear modeling and a nonlinear approach should be chosen when accurate predictions are required. Theory of nonlinear systems has evolved from the seminal work of Rosenberg [1] and is a mature subject nowadays. Kerschen et al. [2] and Peeters et al. [3] attempted to reconcile practice to nonlinear theory. The application of this theory, however, has been restricted to numerical simulations or small experimental rigs, in which the type of nonlinearity is known beforehand. The nonlinear theory has not yet been applied to full-scale engineering structures with multiple components. This paper deals with the nonlinear modal identification of a four-storey medium-size steel frame. Special emphasis is placed on the design of experiments and the subsequent inference of adequate nonlinear models. The aim of the study is to extract very accurate values of the parameters of the models, which will be the basis for future applications on finite element model updating and structural health monitoring.

2 Experimental part

2.1 Structure

The Uniovi structure is a middle-size four-storey steel frame with two bays in the longitudinal direction and one bay in transversal one. The overall dimensions of the structure are 4 m length, 1.5 m width and 7.3 m height (see Figure 1). All columns and beams are HEA-120 and IPN-100, respectively, of steel grade S-275. The floors of the frame are steel sheets 4 mm thick connected to the beams through discontinuous weld. The foundations are two continuous concrete beams lying on the floor of the laboratory.

thumbnail Figure 1.

Structure.

thumbnail Figure 2.

Pendulum exciter.

Each column consists of two pieces. They are spliced through end plates connected by four bolts 12 mm in diameter. The columns are welded to 20 mm thick plates. The beams corresponding to the transversal direction are directly connected to the web of the columns by a welded-all-around fillet. In the longitudinal direction, however, the beams are connected to the flange of the columns. The connection was designed with a 8 mm thick end plate and four bolts 10 mm in diameter.

2.2 Exciter and measurement chain

A pendulum-like exciter was used to excite the frame. It consists of an eccentric mass connected to the shaft of a variable-speed electrical motor (see Figure 2). The excitation was carried out by harmonically varying the angular position of the mass. The frequency and amplitude of the excitation were achieved by controlling the angular position of the shaft through an electronic regulator. The regulator is feedbacked by an encoder. The frequency range of the exciter is from 0 to 100 Hz and the angle range is from 1° to 20°. The total mass of the pendulum is 1 kg.

The structure was provided with four seismic Brüel & Kæjr accelerometers with a sensibility of 10 V/g. Each was screwed at the middle of the end beam on each floor pointed to the longitudinal direction. The hardware used to record the signals was a dSPACE RTI1104 data acquisition card. The processing task was then carried out using MATLAB.

2.3 Testing

Free decay responses of isolated modes are chosen as experimental references for dynamic identification. This approach is very accurate because the excitation uncertainties are avoided. Moreover, free decay vibrations contain a broad range of amplitudes that shows up the nonlinearities present in the structure.

The isolation of a linear mode can be achieved by a monophase multi-point harmonic excitation in quadrature with the response at the related undamped natural frequency. This approach is known as “Appropriation” in literature [4]. The appropriation can be theoretically extended to nonlinear structures by using multi-harmonic excitation [5]. Even though this appropriation is very popular in the aerospace industry, it is not practical for in-service building applications. Instead, an imperfect appropriation consisting of a single-point mono-harmonic excitation was chosen for this case. In reference [6], it is established that this excitation isolates satisfactory a nonlinear modal mode if the structure has well-separated modes.

In the present application, the exciter was placed at the fourth floor. The experiments consisted of 20 s sinusoidal excitation close to the corresponding resonance followed by the free decay vibration. As accelerations are limited by serviceability conditions, only the first two bending modes could be properly excited and measured.

2.4 Signal preprocessing

The analogic signals were converted into digital ones at a sampling frequency of 1 kHz. Then they were filtered in order to remove a variable offset and a periodic high frequency noise present in the raw signals. Finally, signals were truncated selecting the segment corresponding to the last 65 s of the free decay in the first frequency and 25 s, in the second one. The selected periods of time gave rise to a similar decrement in the response amplitude in both cases.

3 Analytical part

3.1 Non-parametric identification

The objective of this section is to choose a family of models that adequately reproduce the dynamic behaviour of the structure on the basis of the experimental results. For this purpose, the evolution of the modal properties of the structure as a function of the response amplitude is studied by means of expeditious procedures.

3.1.1 Mode shapes

The motion of the free vibration in the configuration space corresponding to the preprocessed time series of the acceleration is depicted in Figures 3 and 4. As can be seen, motions exhibit a remarkable linear tend in all the cases. This means that modes do not depend on the acceleration magnitude due to the low level of energy present in the structure. Even though, as will be seen later, both stiffness and damping are weakly nonlinear, mode shapes can be considered linear within the ranges of the experiments. Taking advantage that modes are linear, the motion is defined in the model space in the rest of the paper.

thumbnail Figure 3.

First mode motions in the configuration space.

thumbnail Figure 4.

Second mode motions in the configuration space.

3.1.2 Stiffness and damping

In a low-damped structure the stiffness affects mainly the period of the free vibration, while the amplitude of the response is principally governed by damping. Therefore, the study of the evolution of the period and amplitude of the free response of the structure is the key to understand its dynamic behavior. This is accomplished evaluating independently the period and amplitude of each cycle of the free response in the modal space. The evolution of the frequency as a function of the cycle amplitude is depicted in Figure 5. As can be seen, both modes show a similar trend: the higher the amplitude, the lower the frequency is. Moreover, the frequency tends to a lower limit when the amplitude increases within the range of the experiments. Under these results, an asymptotic softening law is proposed to model the structure stiffness as a function of the modal coordinate.

thumbnail Figure 5.

Frequency vs. amplitude; (a) first mode, (b) second mode.

The damping ratio is proportional to the logarithmic decrement in low-damped structures [7]. Moreover, the logarithmic decrement is defined as the derivative of the logarithm of the amplitude of the free vibration with respect to the number of cycles. Therefore, this derivative is proportional to the damping ratio. The evolution of the amplitude of the cycles is shown in Figure 6. It is found that the logarithm of the amplitude can be expressed as a second order polynomial of the number of cycles. As a consequence of the aforementioned property, the damping ratio can be formulated as a linear function of the number of cycles. As the period of the cycles has a low variation during the free vibration, the damping ratio is approximated by a linear function of time.

thumbnail Figure 6.

Log (A) vs. cycle number. ×: experimental results. Line: second order polynomial fitting; (a) first mode, (b) second mode.

3.2 Parametric identification

3.2.1 Mode shapes and modal response

The values of the mode shapes corresponding to each dof were obtained by proportional correlation of the time series in a least-squares sense. They were normalized with respect to the dof having maximum amplitude. This yields the following solution: ψ i = ( x ̈ r T x ̈ r ) - 1 x ̈ r T x ̈ i , $$ {\psi }_i={({\stackrel{\ddot }{\mathbf{x}}}_r^T{\stackrel{\ddot }{\mathbf{x}}}_r)}^{-1}{\stackrel{\ddot }{\mathbf{x}}}_r^T{\stackrel{\ddot }{\mathbf{x}}}_i, $$(1)in which ψ i is the mode shape coordinate related to i-th dof, x ̈ i $ {\stackrel{\ddot }{\mathbf{x}}}_i$ and x ̈ r $ {\stackrel{\ddot }{\mathbf{x}}}_r$ represent the acceleration time series of the i-th and the reference dofs, respectively. Results are shown in Table 1, the correlation coefficients being close to one.

Table 1.

Mode shapes

Taking advantage of the fact that modes are linear, the motion is defined in the modal space in the rest of the paper through the following least-squares solution: z ̈ = ( Ψ T Ψ ) - 1 Ψ T X ̈ $$ \stackrel{\ddot }{\mathbf{z}}={({\Psi }^T\Psi )}^{-1}{\Psi }^T\stackrel{\ddot }{\mathbf{X}} $$(2)in which the vector z ̈ $ \stackrel{\ddot }{\mathbf{z}}$ denotes the modal time series of acceleration and X ̈ $ \stackrel{\ddot }{\mathbf{X}}$ is a matrix containing the spatial time series of acceleration.

3.2.2 Stiffness and damping

The evaluation of the model parameters of each mode is posed as the minimization of an error function accounting for the discrepancies between the numerical and experimental results. The error function is defined in the time-domain through time series of modal accelerations. The normalized mean squared error is adopted as error function ε = i = 1 M ( z ̈ i ( a ) - z ̈ i ( e ) ) 2 i = 1 M ( z ̈ i ( e ) ) 2 $$ \epsilon =\frac{\sum_{i=1}^M{\left({\ddot{z}}_i^{(a)}-{\ddot{z}}_i^{(e)}\right)}^2}{\sum_{i=1}^M{\left({\ddot{z}}_i^{(e)}\right)}^2} $$(3)where z ̈ i ( a ) $ {\ddot{z}}_i^{(a)}$ and z ̈ i ( e ) $ {\ddot{z}}_i^{(e)}$ represent respectively the model-predicted and the experimental modal accelerations, the latter being obtained through Equation (3). M is the number of available points. The model-predicted response was numerically computed from the equation of motion assuming the following central discrete approximation for the derivatives: z ̇ ( t ) = z ( t + Δ t ) - z ( t - Δ t ) 2 Δ t z ̈ ( t ) = z ( t + Δ t ) - 2 z ( t ) + z ( t - Δ t ) Δ 2 t $$ \begin{array}{c}\dot{z}(t)=\frac{z\left(t+\Delta t\right)-z\left(t-\Delta t\right)}{2\Delta t}\\ \ddot{z}(t)=\frac{z\left(t+\Delta t\right)-2z(t)+z\left(t-\Delta t\right)}{{\Delta }^2t}\end{array} $$(4)where Δt represents the integration period. A value of Δt = 5 × 10−4 s was chosen for the integration. This value constitutes a trade-off between accuracy of results and computation cost.

The minimization is carried out by an adaptive stochastic algorithm developed previously by the authors, which will be referred to as Beta algorithm for the remainder of the paper. This algorithm is described in detail and tested against numerical simulations and experiments in reference [8]. A brief description of the beta algorithm is given here. In this algorithm the solution is sought by sampling the parameters to be identified within a bounded space in an iterative form. The Beta distribution, which is consistent with the bounded character of the search space, is selected for sampling purposes. In each step, the parameters of the distributions are modified accordingly the results of the previous steps. The process is repeated until a given stopping criterion is reached.

Accordingly to the results of the non-parametric identification, the following time-variant linear law is proposed to model the modal damping coefficient c r = c r 1 + c r 2 t , $$ {c}_r={c}_{r1}+{c}_{r2}t, $$(5)in which c r1 and c r2 are parameters to be identified and t stands for time.

The following nonlinear function adapted to the empirical results is proposed to model the modal stiffness: k r ( t ) = k r 1 α + k r 2 | z r ( t ) | α + | z r ( t ) | , $$ {k}_r(t)=\frac{{k}_{r1}\alpha +{k}_{r2}|{z}_r(t)|}{\alpha +|{z}_r(t)|}, $$(6)in which k r1 and k r2 are the bounds of the stiffness and α is a shape parameter. The function is symmetric with respect to the modal coordinate z r . Moreover, its value is k r1 for z r  = 0 and is asymptotic to k r2 when z r tends to infinity; therefore, k r1 represents the stiffness of the underlying linear model. The decrement of the function diminishes when the parameter α increases.

The equation of motion becomes m r z ̈ r ( t ) + ( c r 1 + c r 2 t ) z ̇ r ( t ) + ( k r 1 α + k r 2 | z r ( t ) | α + | z r ( t ) | ) z r ( t ) = 0 . $$ {m}_r{\ddot{z}}_r(t)+\left({c}_{r1}+{c}_{r2}t\right){\dot{z}}_r(t)+\left(\frac{{k}_{r1}\alpha +{k}_{r2}|{z}_r(t)|}{\alpha +\left|{z}_r(t)\right|}\right){z}_r(t)=0. $$(7)

The central differences scheme yields the following solution of the modal response: z r ( t ) = ( 2 Δ 2 t - k r 1 α + k r 2 | z r ( t - Δ t ) α + | z r ( t - Δ t ) | ) z r ( t - Δ t ) 1 Δ 2 t + c r 1 + ( t - Δ t ) c r 1 2 Δ t - ( 1 Δ 2 t - c r 1 + ( t - Δ t ) c r 2 2 Δ t ) z r ( t - 2 Δ t ) 1 Δ 2 t + c r 1 + ( t - Δ t ) c r 2 2 Δ t . $$ {z}_r(t)=\frac{\left(\frac{2}{{\Delta }^2t}-\frac{{k}_{r1}\alpha +{k}_{r2}|{z}_r(t-\Delta t)}{\alpha +\left|{z}_r\left(t-\Delta t\right)\right|}\right){z}_r\left(t-\Delta t\right)}{\frac{1}{{\Delta }^2t}+\frac{{c}_{r1}+\left(t-\Delta t\right){c}_{r1}}{2\Delta t}}-\frac{\left(\frac{1}{{\Delta }^2t}-\frac{{c}_{r1}+\left(t-\Delta t\right){c}_{r2}}{2\Delta t}\right){z}_r\left(t-2\Delta t\right)}{\frac{1}{{\Delta }^2t}+\frac{{c}_{r1}+\left(t-\Delta t\right){c}_{r2}}{2\Delta t}}. $$(8)

The total number of unknown parameters is five in this case. They were achieved through the Beta algorithm. Results are shown in Table 2. The corresponding frequencies of the underlying linear model and those of a pure linear model are shown in Table 3 along side the fitting errors. As can be seen, the fitting errors of the nonlinear model are almost three orders of magnitude lower than those of the pure linear model. This indicates that the nonlinear model is much more precise describing the response of the model. As for the natural frequencies, the discrepancies are around 0.4% in both modes.

Table 2.

Nonlinear softening model results.

Table 3.

Natural frequencies and fitting errors of the nonlinear model and the pure linear model.

4 Conclusions

This study of the Uniovi structure allows a number of conclusions to be drawn. In the experimental part, it was found that the first two bending modes can be individually isolated by means of a single-point mono-harmonic excitation placed at the top floor and remain during the subsequent free decay vibration. The analysis of the motions of the free decays reveals that mode shapes are linear within the acceleration ranges of the experiments. An exploratory study of the evolution of the period and amplitude of the free decay cycles was essential to infer the nonlinearities present in the structure. It was found that the modal damping varies lineally in time and the modal stiffness varies asymptotically as a function of the modal amplitude. The adopted nonlinear model was calibrated by minimizing the discrepancies between the model-predicted responses and the experiments in the time domain. The fitting error of the proposed nonlinear model is near three orders of magnitude lower than that of a pure linear model. A discrepancy around 0.4% was found between the natural frequencies of the pure linear model and those of the linear model underlying the nonlinear one.

Acknowledgments

The authors wish to thank the Science and Innovation Ministry of Spain for its financial support (Projects BIA2006-15266-C02-01 and BIA2006-15266-C02-02). The experimental part of this work has been performed at the laboratory of Fluid Mechanics of the University of Oviedo. The collaboration of the staff of this area, particularly Professor Santolaria, is greatly acknowledged.

References

  1. Rosenberg RM. 1966. On nonlinear vibrations of systems with many degrees of freedom. Advances in Applied Mechanics, 9, 155–242. [CrossRef] [Google Scholar]
  2. Kerschen G, Peeters M, Golinval J-C, Vakakis AF. 2009. Nonlinear normal modes, Part I: A useful framework for the structural dynamicist. Mechanical Systems and Signal Processing, 23, 170–194. [CrossRef] [Google Scholar]
  3. Peeters M, Viguié R, Sérandour G, Kerschen G, Golinval J-C. 2009. Nonlinear normal modes, Part II: Toward a practical computation using numerical continuation techniques. Mechanical Systems and Signal Processing, 23, 195–216. [CrossRef] [Google Scholar]
  4. Ewins DJ. 2000. Modal testing: theory, practice and application. Baldock, Hertfordshire: Research Studies Press Limited, ISBN: 978-0863802188. [Google Scholar]
  5. Peeters M, Kerschen G, Golinval J-C. 2011. Modal testing of nonlinear vibrating structures based on nonlinear normal modes: Experimental demonstration. Mechanical Systems and Signal Processing, 25, 1227–1247. [CrossRef] [Google Scholar]
  6. Peeters M, Kerschen G, Golinval J-C. 2011. Dynamic testing of nonlinear vibrating structures using nonlinear normal modes. Journal of Sound and Vibration, 330, 486–509. [CrossRef] [Google Scholar]
  7. Maia NMM, Silva JMM. 1997. Theoretical and experimental modal analysis. Baldock, Hertfordshire: Research Studies Press Ltd., ISBN: 978-0863802089. [Google Scholar]
  8. Zapico-Valle JL, Alonso-Camblor R, González-Martínez MP, García-Diéguez M. 2010. A new method for finite element model updating in structural dynamics. Mechanical Systems and Signal Processing, 24, 2137–2159. [CrossRef] [Google Scholar]

Cite this article as: Zapico-Valle JL & García-Diéguez M: Dynamic modeling and identification of the Uniovi structure. Int. J. Simul. Multisci. Des. Optim., 2014, 5, A06.

All Tables

Table 1.

Mode shapes

Table 2.

Nonlinear softening model results.

Table 3.

Natural frequencies and fitting errors of the nonlinear model and the pure linear model.

All Figures

thumbnail Figure 2.

Pendulum exciter.

In the text
thumbnail Figure 3.

First mode motions in the configuration space.

In the text
thumbnail Figure 4.

Second mode motions in the configuration space.

In the text
thumbnail Figure 5.

Frequency vs. amplitude; (a) first mode, (b) second mode.

In the text
thumbnail Figure 6.

Log (A) vs. cycle number. ×: experimental results. Line: second order polynomial fitting; (a) first mode, (b) second mode.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.