Next Article in Journal
Properties Evolution of Some Hydraulic Mortars Incorporating Graphene Oxides
Next Article in Special Issue
Research on the Shear Behaviour of Composite Shear Connectors
Previous Article in Journal
Relationship between Project Space Types, Optimize Energy Performance Credit, and Project Size in LEED-NC Version 4 (v4) Projects: A Case Study
Previous Article in Special Issue
Study on Stress Field Distribution during the Construction of a Group of Tunnels Using the Pile–Beam–Arch Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Choices of the Critical Frequency for φ in TL-φ Algorithms When Applied to Multi-Degree of Freedom Systems

1
School of Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China
2
Key Laboratory of Roads and Railway Engineering Safety Control (Shijiazhuang Tiedao University), Ministry of Education, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
3
School of Civil Engineering, Chang’an University, Xi’an 710061, China
*
Author to whom correspondence should be addressed.
Buildings 2022, 12(6), 863; https://doi.org/10.3390/buildings12060863
Submission received: 19 May 2022 / Revised: 16 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022
(This article belongs to the Special Issue Research on Statics and Dynamics of Structures)

Abstract

:
TL- φ algorithms are newly developed explicit structure-dependent integration algorithms utilized for solving the temporally discretized equations of motion. In contrast to the existing algorithms, the most significant improvement of TL- φ algorithms is in diminishing the amount of period errors by introducing a precorrection coefficient φ into the integration parameters of TL- φ algorithms, which is related to the critical frequency of a system. In the previous work, the fundamental frequency of the system is deemed to be the critical frequency, so that φ is a constant scaling corresponding to the fundamental frequency for both single-degree-of-freedom (SDOF) and multi-degree-of-freedom (MDOF) systems. However, for a MDOF system, the first mode may not contribute to the total response more than other ones under a given external excitation, calculating φ only by the fundamental frequency will underestimate the contribution of the higher-frequency modes to structural dynamics. In this paper, choices of the critical frequency for φ when applying TL- φ algorithms to MDOF systems are investigated thoroughly. By considering the initial structural properties of the system and the frequency characteristics of the external excitation simultaneously, a calculation criterion of φ for MDOF systems under specific external excitations is proposed. Four numerical examples with different initial structure properties and loading conditions are designed, and the results demonstrate that the proposed criterion can be potentially used to solve structural dynamic problems of MDOF systems with a more desirable numerical dispersion performance.

1. Introduction

Direct integration algorithms are widely utilized to solve temporally discretized equations of motion (EOMs) in the dynamic analysis of civil engineering structures [1,2,3]. Various integration algorithms have been well established based on different design conceptions in the past decades, such as the finite difference schemes [4,5,6,7], spectral methods of discretization [8,9], and methods based on other ideas [10,11,12]. In general, integration algorithms can be classified as explicit and implicit according to the expressions of displacement and velocity. An algorithm is explicit if the displacement and velocity at the current time step can be expressed in terms of all the known responses at the previous time steps [13], otherwise it is implicit. Explicit algorithms are preferred over implicit ones since they do not involve nonlinear iterations during a step-by-step simulation and possess favorable efficiencies for long-time simulations with small time steps. However, conditional stability of the most explicit algorithms limits their applications in solving dynamic problems, which leads to a time-step size restriction especially for a multi-degree-of-freedom (MDOF) system with high natural frequencies [14]. Hence, integration algorithms that possess both explicit expressions and unconditional stability deserve more attention.
Recently, a new family of integration algorithms classified as structure-dependent integration methods (SDIM) have been well developed and widely used in structural dynamics. In contrast to a structure-independent integration method (SIIM) whose integration parameters of the difference equations are constant, parameters in an SDIM are determined by the initial structural properties, including structural natural frequency, damping ratio and selected time steps [15]. The first SDIM was developed by Chang [16] (named Chang algorithm) for pseudodynamic testing, and is unconditionally stable for linear systems and possesses second-order accuracy. However, this method is semi-explicit as the displacement increment is explicit in form while that of the velocity is implicit. Then, Chen [17] proposed a new SDIM (named CR algorithm) which has explicit expressions for both the displacement and velocity increments. The CR algorithm also demonstrates favorable performance such as unconditional stability for linear systems and second-order accuracy. Rezaiee [18] developed a new algorithm called the USE algorithm. In this method, the expressions of the displacement and velocity increments are defined in terms of the accelerations of the two previous time steps and three integration parameters were introduced to control the stability and numerical dissipation of procedure. Tang [19] presented another type of explicit SDIM known as the Real-time Substructure Testing algorithm (referred to as TL algorithm in this paper), in which the two structure-dependent integration parameters are introduced only in the displacement increment. It is noted that the TL algorithm is unconditionally stable for linear systems and also bears second-order accuracy [20].
The numerical properties of an SDIM are typically evaluated by its stability and accuracy. In order to achieve higher stability limits, many improvements have been conducted for nonlinear systems [21,22,23]. Energy dissipation and numerical dispersion are two important indexes for evaluating the accuracy of an algorithm, which can be measured by amplitude decay (AM) and period elongation (PE) in the time domain. Some researchers [24,25,26,27,28] have developed various families of algorithms that can control the amount of energy dissipation by introducing different control parameters in the integration parameters. Few researchers work on the issue of numerical dispersion of integration algorithms. Li [29] proposed a new method of reducing PE (or numerical dispersion) by introducing a parameter that is weighted by the mass matrix and stiffness matrix into the original algorithm. As only the mass matrix of the system is changed slightly, the convergence rate and stability property of the algorithm are maintained. Tang and Ren [30] presented a new family of explicit SDIM with controllable numerical dispersion based on the TL algorithm, known as TL- φ algorithms. By introducing a scaling parameter φ into the two integration parameters of TL algorithm, the values of PE at a certain structural natural frequency can be greatly diminished when compared to other well-developed algorithms.
In TL- φ algorithms, φ is a precorrection coefficient that is equal to the ratio of the frequency ( ω D ) in the discrete domain to the one ( ω A ) in the continuous domain. This reveals the degree of frequency distortion during the discretization of a continuous system, and introduction of φ in the mapping rules of the discretization will ensure consistency between ω D and ω A in numerical characteristics. For an MDOF system, the calculation of φ depends on the critical frequency ω c = ω A that dominates the total responses of the original continuous system. In the previous work [30], the fundamental frequency ω 1 is regarded as ω c , so that φ becomes a constant scaling for all the modes of the system when calculating the integration parameters of the TL- φ algorithms. However, the first mode may not uniquely contribute to the dynamic responses of an MDOF system under a given external excitation, and the determination of ω c should take the numerical properties of both the original system and external excitation into consideration. In this study, choices of ω c for φ when applying TL- φ algorithms to MDOF systems are thoroughly investigated. Two factors that control the relative importance of any mode in the total dynamic responses are introduced to consider the effect of the initial structural properties of the system and the frequency characteristics of the external excitation on the determination of ω c simultaneously. Then, the selection criterion of ω c for systems subjected to harmonic excitation and seismic waves are proposed, respectively, and studied with respect to the accuracy and efficiency.

2. Formulation of TL- φ Algorithms for MDOF Systems

The formulation of TL- φ algorithms for an N-DOF system can be expressed in discrete form as [30]
M U ¨ i + 1 + C U ˙ i + 1 + K U i + 1 = F i + 1
U i + 1 = U i + α 1 Δ t U ˙ i + α 2 Δ t 2 U ¨ i
U ˙ i + 1 = U ˙ i + Δ t U ¨ i
where, M , C and K are the mass, viscous damping and stiffness matrix, N × N , respectively; U ¨ , U ˙ and U are the acceleration, velocity, and displacement vector, N × 1 , respectively; F is the external force vector, N × 1 ; α 1 and α 2 are the integration parameter matrices, N × N , which are
α 1 = 4 B 1 · M α 2 = B 1 · ( 4 M Δ t C 2 φ · D + E )
where, B = 4 φ 2 · M + 2 Δ t φ · C + Δ t 2 · K ; D = ( Φ T ) 1 · D * 2 · Φ 1 ; E = ( Φ T ) 1 · E * · Φ 1 ; φ = diag [ φ 1 , φ j , , φ N ] , and
φ j = arctan Ω j / 2 Ω j / 2 j = 1 , 2 , . . . , N
in which, Ω j = ω j Δ t , ω j is the natural frequency of the jth mode, and Δ t is the time step selected for integrating.

3. Frequency Response under Arbitrary Excitation

Based on the mode-superposition method, the dynamic response of a linear MDOF system can be obtained by solving separately for that of each mode and then superposing these modes to obtain the total response in the original geometric coordinates. Thus, the numerical characteristics of an algorithm applied in an MDOF system can be evaluated for each mode in the same way as for a single-degree-of-freedom (SDOF) system. The numerical characteristics of TL- φ algorithms under the free-vibration response have been fully investigated in the previous work [30]. This section focuses on the frequency response of TL- φ algorithms under arbitrary excitation with initial conditions of x 0 = 0 and x ˙ 0 = 0 . For any arbitrary mode n, the equation of motion of an SDOF system under external excitation is
y ¨ n ( t ) + 2 ξ n ω n y ˙ n ( t ) + ω n 2 y n ( t ) = f n ( t ) m n n = 1 , 2 , , N ξ n = c n 2 m n ω n , ω n = k n m n
where, m n = ϕ n T M ϕ n , c n = ϕ n T C ϕ n and k n = ϕ n T K ϕ n are modal mass, modal damping and modal stiffness for the nth mode, respectively; y ¨ n , y ˙ n and y n are acceleration, velocity, and displacement for the nth mode, respectively; ξ n and ω n are the natural frequency and viscous damping ratio for the nth mode, respectively; f n = ϕ n T F ( t ) is modal force for the nth mode. For simplicity, the subscript n will be omitted from this point onward.
The continuous transfer function G ( s ) that relates the force f ( t ) to the displacement y ( t ) is given by
G ( s ) = Y ( s ) F ( s ) = 1 m ( s 2 + 2 ξ ω + ω 2 )
where, Y ( s ) and F ( s ) are the Laplace transforms of y ( t ) and f ( t ) , respectively, and s is the Laplace transform parameter.
Substituting s = i ω ˜ into Equation (7), the complex frequency response function of the steady-state response under the harmonic excitation of f ( t ) = f 0 e i ω ˜ t can be determined by
G ( i ω ˜ ) = Y ( i ω ˜ ) F ( i ω ˜ ) = f 0 k [ ( 1 β 2 ) + i ( 2 ξ β ) ] , β = ω ˜ / ω
in which, β is frequency ratio, that is, the ratio of the excitation frequency ω ˜ to the modal frequency ω . The dynamic magnification factor (DMF) D and phase θ of G ( i ω ˜ ) are
D = G ( i ω ˜ ) y st = 1 ( 1 β 2 ) 2 + ( 2 ξ β ) 2
θ = tan 1 2 ξ β 1 β 2
where y st = f 0 / k is the displacement which would be produced by the load f 0 applied statically. The two factors, D and θ , calculated by Equations (9) and (10) denote the exact dynamic magnification factor and phase, respectively.
Similarly, the complex frequency response function for TL- φ algorithms in the discrete domain can be obtained by introducing z = e i ω ˜ Δ t into the discrete transfer function G ( z ) , which is
G ( z ) = Y ( z ) F ( z ) = α 2 Δ t 2 z + ( α 1 α 2 ) Δ t 2 m z 2 + m ( α 2 Ω 2 + 2 ξ Ω 2 ) z + m [ ( α 1 α 2 ) Ω 2 2 ξ Ω + 1 ]
and the numerical DMF and phase are expressed as D num = k G ( e i ω ˜ Δ t ) and θ num = G ( e i ω ˜ Δ t ) , respectively.
Figure 1 shows D num and θ num for TL- φ algorithms when compared to the exact solution of D and θ . The damping ratio ξ is taken as 0.1, while the different values of Ω = 0.2, 1.0 and 2.0 are considered and the corresponding coefficients φ become 0.9967, 0.9273 and 0.7854, respectively. Meanwhile, D num and θ num for the other three well-developed explicit integration algorithms, including Chang algorithm, CR algorithm and TL algorithm, are also calculated and shown in each sub-figure.
It can be seen in Figure 1a that D num and θ num for all the integration algorithms agree well with their accurate counterparts when Ω = 0.2 . For a larger value of Ω = 1.0 and 2.0, the numerical resonance peak shifts to the left for all the algorithms, and the larger the value of Ω , the further left the numerical resonance peak shifts. It is also observed that the left-shift values of the numerical resonance peak of TL- φ algorithms are smaller than those of other algorithms, indicating a smaller PE in the numerical solution. This property conforms to the findings presented in previous work [30]. In addition, the maximum values of D num for all the algorithms become larger with an increasing value of Ω than those of the exact solution.

4. Critical Frequency Selection for φ in TL- φ Algorithms

When applying TL- φ algorithms to a MDOF system, the first and most important issue is the calculation of φ . Precorrection for each mode of the system may possess high accuracy but is time consuming, especially when the number of DOFs is huge. It is quite promising that determination of a critical frequency ω c for φ , reduces φ to a constant scaling for all the modes.
In general, the first mode of an MDOF system possesses a dominant impact on the total responses of free-vibration and forced-vibration under specific external excitation, such as an earthquake tends to excite the response of a structure mainly in its first mode [31]. However, higher modes of an MDOF system may contribute more to the response than the first one, and the determination of ω c should consider both the initial structural properties of the system and the frequency characteristics of the external excitation. It should be noted that the importance of any mode in the total dynamic response depends on two factors: (1) The DMF that depends on the frequency ratio of the external excitation to the mode, as discussed in the previous section; (2) the modal participation factor (MPF) that is determined by the interaction of the mode shape with the spatial distribution of the external excitation.

4.1. Determination of ω c Based on DMF

According to the frequency response analysis, the maximum response amplitude occurs at a frequency ratio β approaching close to (slightly less than) unity, i.e., the external excitation frequency approaches a certain modal frequency, known as resonance. For an MDOF system, an arbitrary external excitation will excite the resonance response of any mode, and the determination of ω c should consider β i that relates the excitation frequency ω ˜ to the ith mode frequency ω i .
In this section, a spring-mass system with 20-DOF is set up to investigate the relationship between ω c and ω ˜ , as shown in Figure 2. A total of 12 models with different structural properties are intentionally designed, and the fundamental frequencies of each model are shown in Table 1. The system is assumed to be linear elastic and is subjected to a ground acceleration u ¨ g ( t ) of sine-wave form as follows,
u ¨ g ( t ) = sin ( ω ˜ t )
and then, F i + 1 in Equation (1) is expressed as
F i + 1 = M I · u ¨ g i + 1
TL- φ algorithms with ω c = 1 , 1.5 , 2 , 2.5 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10 , 11 , 12 for φ are applied to calculate the displacement responses of the system under the harmonic excitation with ω ˜ = 1 , 1.5 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10 , respectively. One error indicator, normalized root-mean-square error (NRMSE) that is sensitive to frequency differences of the dynamic responses, is adopted for accuracy analysis and defined as
NRMSE = i = 1 m ( u ES , i u NS , i ) 2 m max ( u NS , i ) min ( u NS , i )
where, u is the displacement of m 20 ; subscript NS and ES represent the numerical solution obtained by TL- φ algorithms and exact solution derived from the mode-superposition method, respectively; m = t / Δ t is the number of the samples, where t is the duration of the external excitation and Δ t = 0.01 s is the time step for numerical simulations.
Plotted in Figure 3a–c are the variations of NRMSE with the increasing of ω c for models 1, 8 and 10, and similar curves can be obtained for the rest of the models. It can be seen that the NRMSE is minimal at: (1) ω c = ω 1 when β 1 < 1 , i.e., ω ˜ < ω 1 ; (2) ω c = ω i when β i 1 , i.e., ω ˜ ω i ( i = 1 , 2 , . . . , 10 ) ; (3) otherwise, ω c is between values of ω 1 and ω ˜ , which verifies the importance of ω 1 and ω ˜ in determining ω c for the calculation of φ in TL- φ algorithm. In order to determine the value of ω c in (3), more numerical models with different values of ω 1 (less than 10 Hz) are designed to calculate the optimal critical frequencies ω ^ c that correspond to the minimal NRMSE. To diminish the effect of the fundamental frequency randomly selected for numerical simulation on fitting results, arithmetic mean values of ω ^ c for each model are calculated according to the following formula,
ω ¯ c ω 1 = 1 n ω ^ c ω 1 , ω ¯ c ω ˜ = 1 n ω ^ c ω ˜
in which, n is the number of the harmonic excitation ω ˜ .
Table 2 lists the ratio of the mean values of ω ¯ c to ω 1 and ω ˜ for the 15 models. Figure 4 shows the scatter diagram of mean values varying with ω 1 , which indicates that ω ¯ c / ω 1 decay exponentially with an increasing of ω 1 , while ω ¯ c / ω ˜ increases linearly with an increasing of ω 1 . According to these rules, a regression analysis is conducted and the fitting formulas are
ω ˜ c 1 ω 1 = 1 + 86 e 10.7 ω 1 + 7.8 e 1.04 ω 1 , 0 < ω 1 < 3 1 , ω 1 3
ω ˜ c 2 ω ˜ = 0.71 , 0 < ω 1 < 2.5 0.32 ω 1 + 0.1 , ω 1 2.5
ω ˜ c = ω ˜ c 1 , ω ˜ c 1 ω 1 ω ˜ c 2 ω ˜ ω ˜ c 2 , ω ˜ c 1 ω 1 > ω ˜ c 2 ω ˜
where, ω ˜ c 1 and ω ˜ c 2 are predicted values of ω ¯ c in Equation (15). The smaller one of the two radios derived from Equations (16) and (17) is considered to calculate ω ˜ c . Then, the criterion of determining ω c for the calculation of φ in TL- φ algorithm when applied to MDOF systems can be expressed as
ω c = ω 1 , β 1 < 1 ω i , β i = 1 ± 0.1 ω ˜ c , otherwise
When an MDOF system is subjected to an earthquake, u ¨ g ( t ) in Equation (12) represents a practical seismic record. It is known that the contribution of the first mode to the dynamic responses of the system is much greater than that of the other higher-order modes when the fundamental frequency of the system is close to or higher than the main frequency component of a seismic motion. In addition, defining ω c = ω 1 will achieve relatively high accuracy of the dynamic responses. However, the effect of the high-order modes on structural dynamics cannot be ignored when the fundamental frequency of the system is lower than the main frequency component of a seismic motion. Thus, it is of great interest to examine the influence of the spectral characteristic of a seismic motion on the determination of ω c . Then, ω ˜ will be replaced by the peak frequency of the response spectrum ω p in the following.
Numerical analysis of the 20-DOF spring-mass system is performed to establish a selection criterion of ω c . A total of 14 seismic waves are selected as external excitations, and ω p of each wave is shown in Table 3. Moreover, 10 groups of structural properties are designed to investigate the relationship between ω c and ω p as well as ω 1 .
TL- φ algorithms with ω c varying from 1 rad/s to 20 rad/s are applied to calculate NRMSE for each model. It is worth noting that the minimal NRMSE can be obtained when ω c = ω 1 for the case of ω p ω 1 . For ω p > ω 1 , the optimal critical frequencies ω ^ c and the arithmetic mean values ω ¯ c are acquired according to Equation (15). Table 4 lists the ratio of the mean values of ω ¯ c to ω 1 and ω p for each model. Figure 5 shows the scatter diagram of mean values varying with ω 1 , which indicates that ω ¯ c / ω 1 decays exponentially with an increasing of ω 1 , while ω ¯ c / ω p increases linearly with an increasing of ω 1 . The fitting formulas are also conducted based on the regression analysis as follows, and the ω c that equals to the smaller one of ω ˜ c 1 and ω ˜ c 2 is recommended in this situation.
ω ˜ c 1 ω 1 = 4.2 e 0.3 ω 1 + 0.4 e 0.05 ω 1 , 1.2 ω 1 20
ω ˜ c 2 ω p = 0.037 ω 1 + 0.12 , 1.2 ω 1 20
ω ˜ c = min ( ω ˜ c 1 , ω ˜ c 2 )
ω c = ω 1 , ω p ω 1 ω ˜ c , ω p > ω 1

4.2. Modal Participation Factor, MPF

Another factor that rules the importance of any mode in the total dynamic response is MPF. In Equation (6), the external force vector F ( t ) can be expressed as a product of a distribution vector R and an amplitude vector p ( t ) :
F ( t ) = R p ( t )
and for the particular case of seismic excitation, it has
F ( t ) = Mr u ¨ g ( t )
where, r is a displacement transformation vector that expresses the displacement of each DOF excited by a unit static support displacement.
Introducing Equation (24) or (25) into Equation (6), the equation of motion for each DOF becomes:
y ¨ n ( t ) + 2 ξ n ω n y ˙ n ( t ) + ω n 2 y n ( t ) = ϕ n T R ϕ n T M ϕ n p ( t )
or
y ¨ n ( t ) + 2 ξ n ω n y ˙ n ( t ) + ω n 2 y n ( t ) = ϕ n T Mr ϕ n T M ϕ n u ¨ g ( t )
The ratios shown on the right side of Equations (26) and (27) are defined as the modal participation factor. That is,
MPF n = ϕ n T R ϕ n T M ϕ n or MPF n = ϕ n T Mr ϕ n T M ϕ n
It is apparent from Equation (28) that the amplitude of the response due to any given mode depends on the interaction of the applied load distribution and the mode shape. An arbitrary external load distribution applied to an engineering structure might excite a response in any of the modes, and it is recommended that ω c equals to the frequency corresponding to the mode with the largest value of MPF. While for an engineering structure subjected to horizontal ground motion, the transformation vector r is a unit column, and the determination of ω c depends mainly on DMF.

5. Numerical Confirmations

For numerical confirmation of the selection criterion for ω c proposed in this paper, three cases with different original structural properties and loading conditions are intentionally designed.

5.1. Case 1: Undamped Linear MDOF System under Harmonic Ground Motion

A spring-mass system with 50-DOF is utilized to discuss the effect of ω c selected for φ on dynamic responses, and its structural properties are plotted in Figure 6. After the eigenvalue-analysis, the lowest natural frequency of the system is found to be 3.1104 rad/s, while the highest natural frequency is found to be 199.90 rad/s. Table 5 shows the first 10 natural frequencies and the corresponding MPFs of the system. The system is subjected to a harmonic acceleration excitation of sine-wave form at the ground, and the circular frequency ω ˜ is set to 1, 6, 16 and 40 rad/s, respectively.
According to the selection criterion of ω c in Equation (19), ω c = ω 1 when ω ˜ = 1 and 6 rad/s; ω c = ω 3 when ω ˜ = 16 rad/s; ω c = ω 7 when ω ˜ = 40 rad/s, and the corresponding φ = 0.9999 , 0.9980 , 0.9869 , respectively. Numerical results obtained from TL- φ algorithm with different values of ω c and other three well-developed algorithms, including Chang algorithm, CR algorithm and TL algorithm, with Δ t = 0.01 s are shown in Figure 7a–d. In the figures, the analytical solutions obtained by modal superposition method (MSM) are marked as “Exact” for comparison, and “TL- φ ” refers to the ones derived from TL- φ algorithm with precorrecting for each mode. It is worth noting that the displacements obtained from TL- φ algorithm are nearly in total agreement with the exact solution for all the loading conditions while those obtained from the other three algorithms possess obvious period distortions and amplitude errors. Table 6, Table 7, Table 8 and Table 9 list two error indicators, including the aforementioned NRMSE and normalized energy error (NEE), for displacement responses of the system ( m 50 ) under harmonic excitation with different ω ˜ . In general, NEE reflects signal energy and is sensitive to amplitude differences of the dynamic responses, defined as
NEE = i = 1 m u ES , i 2 i = 1 m u NS , i 2 i = 1 m u NS , i 2
where, the meaning of the subscript NS, ES and m are the same as those for NRMSE in Equation (14).
It is worth noting that both NEE and NRMSE are the minimum for “TL- φ ” in any case, confirming that precorrecting for each mode of an MDOF system will achieve great accuracy. In addition, the value of ω c chosen for φ has a noticeable effect on the two errors, especially when the system suffers from a rather high excitation frequency. An appropriate value of ω c chosen for φ will result in relatively smaller errors in both amplitude and period when compared with other algorithms. This verifies the effectiveness of the proposed criterion in selecting ω c when applying the TL- φ algorithm to MDOF systems.

5.2. Case 2: Undamped Linear MDOF System under Arbitrary Load

A five-story shear building under an arbitrary load is shown in Figure 8. Two groups of the structural property are intentionally designed to investigate the effect of MPF on determining ω c , including (1) m i = 10 5 kg, k i = 10 7 N/m ( i = 1 , 2 5 ); (2) m 1 = 10 7 kg, m 2 = m 3 = 10 5 kg, m 4 = m 5 = 10 3 kg, and k i = 10 7 N/m ( i = 1 , 2 5 ). The initial natural frequencies and the modal shapes of the building for the two groups (labelled G1 and G2) are shown in Table 10 and Table 11.
The building is excited by an arbitrary load F = [ f 1 , f 2 , f 3 , f 4 , f 5 ] T with a spatial load distribution vector of R , i.e., F ( t ) = R · p 0 sin ( ω ˜ t ) . In this case, p 0 = 100 kN and ω ˜ = 3 rad/s are employed. Considering R = ϕ 2 , ϕ 3 and ϕ 5 , this will accentuate the responses in the 2nd, 3rd and 5th mode of the building, respectively. The corresponding ω c selected in this case should be equal to ω 2 , ω 3 and ω 5 based on the proposed criterion of MPF, and then φ = 0.9977 , 0.9943 , 0.9880 for G1 while φ = 0.9987 , 0.9914 , 0.6285 for G2, respectively. The displacement responses obtained from the Newmark algorithm with [ α , β ] = [1/2, 1/4] and Δ t = 0.005 s are considered as the “Exact” solution for comparison, while the time step of Δ t = 0.02 s is employed to compute the seismic responses for using TL- φ algorithms and TL algorithm. Numerical results for the initial conditions of U = 0 m, U ˙ = 0 m/s are displayed in the following. Figure 9 shows the displacement responses of the 5th story for G1 while those for G2 are plotted in Figure 10. Errors of the displacement responses for G1 and G2 are listed in Table 12 and Table 13.
Regarding G1, it can be seen from Figure 9a–c that the displacements derived from the TL- φ algorithm with ω c = ω 2 for R = ϕ 2 , ω c = ω 3 for R = ϕ 3 and ω c = ω 5 for R = ϕ 5 are almost in total agreement with the “Exact” solution. Displacement errors shown in Table 12 distinctly indicate that a much higher accuracy can be obtained when ω c is equal to the frequency corresponding to the mode with the largest value of MPF. For G2, it can be seen in Figure 10a–c and Table 13 that TL- φ algorithm with ω c = ω ˜ c can result in a better solution while larger errors are found in the results obtained from TL- φ algorithm with ω c = ω i ( i = 2 , 3 , 5 ) . This illustrates that the availability of determining ω c based on MPF depends on the mass distribution of the system. The more evenly the mass is distributed, the greater the influence of MPF on determining ω c . On the contrary, the determination of ω c should depend mainly on DMF.

5.3. Case 3: Damped Linear MDOF System under Seismic Input

A four-story shear building under seismic input is shown in Figure 11. The structural properties are m 1 = 10 5 kg, m 2 = m 3 = 10 4 kg, m 4 = 10 3 kg, k i = 10 7 N/m ( i = 1 , 2 , 3 , 4 ), and the natural frequencies are 8.88, 21.08, 50.65, and 105.41 rad/s. Assuming Rayleigh Damping, the damping matrix can be determined by combining the mass and the stiffness matrices as follows,
C = a 0 M + a 1 K
in which, a 0 and a 1 are the two proportionality factors evaluated by the following equation:
a 0 a 1 = 2 ξ ω 1 + ω 3 ω 1 ω 3 1
where, ξ = 0.05 is the damping ratio; ω 1 and ω 3 are the two specific frequencies selected for Rayleigh Damping.
The system is subjected to the ground acceleration record of the El Centro wave, where the peak ground acceleration (PGA) is scaled to 0.5 g. The sampling interval of the record is 0.02 s. The time history and response spectrum of the El Centro wave are shown in Figure 12, noting that the peak frequency of the response spectrum ω p = 24.166 rad/s.
Numerical results for the initial conditions of u ( 0 ) = 0 m and u ˙ ( t ) = 0 m/s are displayed in Figure 13, where the numerical solution obtained from the Newmark algorithm with [ α , β ] = [1/2, 1/4] and Δ t = 0.005 s is considered as a reference solution for comparison. TL- φ algorithms with the time step of 0.02 s are applied to perform the step-by-step integration. Based on the proposed selection criterion in Equations (20)–(22), ω c = min ( ω ˜ c 1 , ω ˜ c 2 ) = 8.14 rad/s and φ = 0.9978 . Displacement errors of the 4th story derived from TL- φ algorithms with different values of ω c are shown in Table 14. It can be seen that both NEE and NRMSE are minimal for “TL- φ ( ω c = 8.14 rad/s)”, indicating that the proposed selection criterion of ω c is applicable for an MDOF that is subjected to a seismic wave.

5.4. Case 4: An MDOF System with Softening Stiffness

The five-story shear building shown in Figure 8 is used in this case. The initial structural properties are m i = 10 5 kg, k i , 0 = 10 7 N/m ( i = 1 , 2 5 ), ξ = 0.02 , and the natural frequencies are as shown in Table 10. The softening stiffness of the system after deformation is designated as k i = k i , 0 ( 1 u i u i 1 ) , where subscript i refers to the ith story and ( u i u i 1 ) is the inter-story drift of the ith story.
The system with zero initial conditions is subjected to the ground acceleration record of the Loma Prieta wave, where PGA = 1.05 g. The sampling interval of the record is 0.01 s. The time history and response spectrum of the Loma Prieta wave are shown in Figure 14, noting that the peak frequency of the response spectrum ω p = 10.134 rad/s.
TL- φ algorithms with the time step of 0.01 s are applied to perform the step-by-step integration, while the numerical solution obtained from the Newmark algorithm with [ α , β ] = [1/2, 1/4] and Δ t = 0.005 s is considered as a reference solution for comparisons. According to the proposed selection criterion in Equations (20)–(22), ω c = min ( ω ˜ c 1 , ω ˜ c 2 ) = 2.28 rad/s and φ 1.0000 . Displacement errors of the 4th story derived from TL- φ algorithms with different values of ω c are shown in Table 15, and the responses are ploted in Figure 15. It can be seen that both NEE and NRMSE are minimal for “TL- φ ( ω c = 2.28 rad/s)”, indicating that the proposed selection criterion of ω c is still applicable for non-linear MDOF systems.

6. Conclusions

TL- φ algorithms represent newly proposed SDIM that possess controllable numerical dispersion properties by introducing a precorrection parameter φ . This precorrection parameter depends on the critical frequency ω c of a system. When TL- φ algorithms are applied to an MDOF system, the first and most important issue is the determination of ω c for φ . Precorrection for each frequency (mode) of the system may achieve high accuracy, but is time consuming especially when the number of DOFs is huge. Indeed, relatively high accuracy of the dynamic responses can still be obtained by precorrecting a specific mode (deemed to be ω c ) that dominates the total responses of an MDOF system, and then φ will reduce to a constant scaling φ for all the modes of the system.
In this work, the choices of ω c are fully discussed by investigating the relationship between ω c and ω 1 as well as ω ˜ (or ω p ). It is worth noting that the value of ω c selected for φ has a significant influence on the structural responses, and its determination should consider both the initial structural properties of the system and the numerical characteristics of external excitation. Two factors that control the importance of any mode in the total dynamic responses, including DMF and MPF, are taken into account to obtain the optimal critical frequency. Empirical formulas of ω c for systems subjected to both harmonic excitation and seismic waves are obtained by numerical calculation of NRMSE. A total of four numerical cases are intentionally designed to examine the feasibility of the selection criterion ω c proposed in this work. By studying the two error indicators of NEE and NRMSE, it is illustrated that a smaller period distortion in dynamic responses derived from TL- φ algorithms can be achieved when ω c is determined according to the proposed criterion. This will lead to further improvements of the numerical dispersion property of TL- φ algorithms in solving structural dynamic problems of MDOF systems.

Author Contributions

Conceptualization, methodology, formal analysis, writing—original draft preparation, Y.T.; software, data curation, C.L.; validation, writing—review and editing, B.F.; funding acquisition, Y.T., C.L. and B.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Foundation of China (Grant No. 52008074, 51908048), the Open Fund of Key Laboratory of Roads and Railway Engineering Safety Control (Shijiazhuang Tiedao University), Ministry of Education (Grant No. STDTKF202004) and Hebei Province Key R&D Program (Grant No. 21375407D).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
EOMequation of motion
SDOFsingle-degree-of-freedom
MDOFmulti-degree-of-freedom
SDIMstructure-dependent integration method
SIIMstructure-independent integration method
AMamplitude decay
PEperiod elongation
DMFdynamic magnification factor
MPFmodal participation factor
NRMSEnormalized root-mean-square error
NEEnormalized energy error
MSMmodal superposition method
PGApeak ground acceleration

References

  1. Darby, A.; Blakeborough, A.; Williams, M. Improved control algorithm for real-time substructure testing. Earthq. Eng. Struct. Dyn. 2001, 30, 431–448. [Google Scholar] [CrossRef]
  2. Bonnetr, P.; Williams, M.; Blakeboroughg, A. Evaluation of numerical time-integration schemes for real-time hybrid testing. Earthq. Eng. Struct. Dyn. 2008, 37, 1467–1490. [Google Scholar] [CrossRef]
  3. Chang, S. Choices of structure-dependent pseudodynamic algorithms. J. Eng. Mech. 2019, 145, 04019029. [Google Scholar] [CrossRef]
  4. Newmark, N. A method of computation for structural dynamics. J. Eng. Mech. 1959, 85, 67–94. [Google Scholar] [CrossRef]
  5. Wilson, E.; Farhoomand, I.; Bathe, K. Nonlinear dynamic analysis of complex structure. Earthq. Eng. Struct. Dyn. 1973, 1, 241–252. [Google Scholar] [CrossRef]
  6. Hilber, H.; Hughes, T.; Taylor, R. Improved numerical dissipation for time integration algorithms in structural mechanics. Earthq. Eng. Struct. Dyn. 1977, 5, 283–292. [Google Scholar] [CrossRef] [Green Version]
  7. Butcher, J. A history of Runge–Kutta methods. Appl. Numer. Math. 1996, 20, 247–260. [Google Scholar] [CrossRef]
  8. Bhrawy, A. An efficient Jacobi pseudospectral approximation for nonlinear complex generalized Zakharov system. Appl. Math. Comput. 2014, 247, 30–46. [Google Scholar] [CrossRef]
  9. Momeni, M.; Riahi Beni, M.; Bedon, C.; Najafgholipour, M.; Dehghan, S.; JavidSharifi, B.; Hadianfard, M. Dynamic Response Analysis of Structures Using Legendre–Galerkin Matrix Method. Appl. Sci. 2021, 11, 9307. [Google Scholar] [CrossRef]
  10. Ertürk, V.; Momani, S. Solving systems of fractional differential equations using differential transform method. J. Comput. Appl. Math. 2008, 215, 142–151. [Google Scholar] [CrossRef] [Green Version]
  11. Pakdemirli, M.; Karahan, M.M.F.; Boyacı, H. Forced Vibrations of Strongly Nonlinear Systems with Multiple Scales Lindstedt Poincaré Method. Math. Comput. Appl. 2011, 16, 879–889. [Google Scholar] [CrossRef] [Green Version]
  12. Rysak, A.; Magdalena, G. Differential Transform Method as an Effective Tool for Investigating Fractional Dynamical Systems. Appl. Sci. 2021, 11, 6955. [Google Scholar] [CrossRef]
  13. Mahin, S.; Shing, P. Pseudodynamic method for seismic testing. J. Eng. Mech. 1985, 111, 1482–1509. [Google Scholar] [CrossRef]
  14. Tang, Y.; Qin, H. Stability and Accuracy Analysis of Real-Time Hybrid Simulation (RTHS) with Incomplete Boundary Conditions and Actuator Delay. Int. J. Struct. Stab. Dyn. 2020, 20, 2050122. [Google Scholar] [CrossRef]
  15. Chang, S.; Wu, T.; Tran, N.; Yang, Y. Apllications of a Family of Unconditionally stable, Dissipative, Explicit methods to pseudodynamic tests. Exp. Tech. 2017, 41, 19–36. [Google Scholar] [CrossRef]
  16. Chang, S. Explicit pseudodynamic algorithm with unconditional stability. J. Eng. Mech. 2002, 128, 935–947. [Google Scholar] [CrossRef]
  17. Chen, C.; Ricles, J. Development of direct integration algorithms for structural dynamics using discrete control theory. J. Eng. Mech. 2008, 134, 676–683. [Google Scholar] [CrossRef]
  18. Rezaiee-Pajand, M.; Hashemian, M. Time Integration Method Based on Discrete Transfer Function. Int. J. Struct. Stab. Dyn. 2016, 16, 1550009. [Google Scholar] [CrossRef]
  19. Tang, Y.; Lou, M. New unconditionally stable explicit integration algorithm for real-time hybrid testing. J. Eng. Mech. 2017, 143, 04017029. [Google Scholar] [CrossRef]
  20. Tang, Y.; Qin, H. Applications of the RST Algorithm to Nonlinear Systems in Real-Time Hybrid Simulation. Math. Probl. Eng. 2020, 2020, 5734720. [Google Scholar] [CrossRef]
  21. Chang, S. Improved explicit method for structural dynamics. J. Eng. Mech. 2007, 133, 748–760. [Google Scholar] [CrossRef]
  22. Chang, S. An explicit method with improved stability property. Int. J. Numer. Methods Eng. 2009, 77, 1100–1120. [Google Scholar] [CrossRef]
  23. Wang, T.; Zhou, H.; Zhang, X.; Ran, T. Stability of an explicit time-integration algorithm for hybrid tests, considering stiffness hardening behavior. Earthq. Eng. Struct. Vib. 2018, 17, 595–606. [Google Scholar] [CrossRef]
  24. Gui, Y.; Wang, J.; Jin, F.; Chen, C.; Zhou, M. Development of a family of explicit algorithms for structural dynamics with unconditional stability. Nonlinear Dyn. 2014, 77, 1157–1170. [Google Scholar] [CrossRef]
  25. Kolay, C.; Ricles, J. Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical numerical energy dissipation. Earthq. Eng. Struct. Dyn. 2014, 43, 1361–1380. [Google Scholar] [CrossRef]
  26. Kim, W.; Lee, J. An improved explicit time integration method for linear and nonlinear structural dynamics. Comput. Struct. 2018, 206, 42–53. [Google Scholar] [CrossRef]
  27. Fu, B.; Feng, D.; Jiang, H. A new family of explicit model-based integration algorithms for structural dynamic analysis. Int. J. Struct. Stab. Dyn. 2019, 19, 1950053. [Google Scholar] [CrossRef]
  28. Li, J.; Yu, K. Noniterative Integration Algorithms with Controllable Numerical Dissipations for Structural Dynamics. Int. J. Comput. Methods 2019, 16, 1850111. [Google Scholar] [CrossRef]
  29. Li, S.; Qin, L.; Guo, H.; Yang, D. A Method of Improving Time Integration Algorithm Accuracy for Long-Term Dynamic Simulation. Int. J. Struct. Stab. Dyn. 2020, 20, 2050079. [Google Scholar] [CrossRef]
  30. Tang, Y.; Ren, D.; Qin, H.; Luo, C. New Family of Explicit Structure-Dependent Integration Algorithms with Controllable Numerical Dispersion. J. Eng. Mech. 2021, 147, 04021001. [Google Scholar] [CrossRef]
  31. Clough, R.; Penzien, J. Dynamics of Structures, 2nd ed.; Higher Education Press: Beijin, China, 2011. [Google Scholar]
Figure 1. Numerical DMF D num and phase θ num for TL- φ algorithms compared with the exact solution of dynamic magnification factor D and phase θ . (a) Ω = 0.2 . (b) Ω = 1.0 . (c) Ω = 2.0 .
Figure 1. Numerical DMF D num and phase θ num for TL- φ algorithms compared with the exact solution of dynamic magnification factor D and phase θ . (a) Ω = 0.2 . (b) Ω = 1.0 . (c) Ω = 2.0 .
Buildings 12 00863 g001aBuildings 12 00863 g001b
Figure 2. A 20-DOF spring-mass system.
Figure 2. A 20-DOF spring-mass system.
Buildings 12 00863 g002
Figure 3. NRMSE curves vs. critical frequency ω c under harmonic excitation with different values of ω ˜ . (a) Model 1. (b) Model 8. (c) Model 10.
Figure 3. NRMSE curves vs. critical frequency ω c under harmonic excitation with different values of ω ˜ . (a) Model 1. (b) Model 8. (c) Model 10.
Buildings 12 00863 g003
Figure 4. Fitting curves for optimal critical frequency (sine-wave).
Figure 4. Fitting curves for optimal critical frequency (sine-wave).
Buildings 12 00863 g004
Figure 5. Fitting curves for optimal critical frequency (seismic wave).
Figure 5. Fitting curves for optimal critical frequency (seismic wave).
Buildings 12 00863 g005
Figure 6. A 50-DOF spring-mass system.
Figure 6. A 50-DOF spring-mass system.
Buildings 12 00863 g006
Figure 7. Displacement responses for m 50 . (a) ω ˜ = 1 rad/s. (b) ω ˜ = 6 rad/s. (c) ω ˜ = 16 rad/s. (d) ω ˜ = 40 rad/s.
Figure 7. Displacement responses for m 50 . (a) ω ˜ = 1 rad/s. (b) ω ˜ = 6 rad/s. (c) ω ˜ = 16 rad/s. (d) ω ˜ = 40 rad/s.
Buildings 12 00863 g007
Figure 8. Schematic diagram of a five-story shear frame building without damping.
Figure 8. Schematic diagram of a five-story shear frame building without damping.
Buildings 12 00863 g008
Figure 9. Displacement responses of the 5-story building (G1). (a) R = ϕ 2 . (b) R = ϕ 3 . (c) R = ϕ 5 .
Figure 9. Displacement responses of the 5-story building (G1). (a) R = ϕ 2 . (b) R = ϕ 3 . (c) R = ϕ 5 .
Buildings 12 00863 g009
Figure 10. Displacement responses of the 5-story building (G2). (a) R = ϕ 2 . (b) R = ϕ 3 . (c) R = ϕ 5 .
Figure 10. Displacement responses of the 5-story building (G2). (a) R = ϕ 2 . (b) R = ϕ 3 . (c) R = ϕ 5 .
Buildings 12 00863 g010
Figure 11. Schematic diagram of a four-story shear frame building with damping.
Figure 11. Schematic diagram of a four-story shear frame building with damping.
Buildings 12 00863 g011
Figure 12. Spectral characteristics of the El Centro wave. (a) Time history. (b) Response spectrum.
Figure 12. Spectral characteristics of the El Centro wave. (a) Time history. (b) Response spectrum.
Buildings 12 00863 g012
Figure 13. Displacement responses of the 4th story.
Figure 13. Displacement responses of the 4th story.
Buildings 12 00863 g013
Figure 14. Spectral characteristics of Loma Prieta wave. (a) Time history. (b) Response spectrum.
Figure 14. Spectral characteristics of Loma Prieta wave. (a) Time history. (b) Response spectrum.
Buildings 12 00863 g014
Figure 15. Displacement responses of the 4th story.
Figure 15. Displacement responses of the 4th story.
Buildings 12 00863 g015
Table 1. The fundamental frequencies of the models ( rad / s ).
Table 1. The fundamental frequencies of the models ( rad / s ).
Model123456
ω 1 0.06990.20470.44230.76611.08341.7129
Model789101112
ω 1 2.42253.42595.41687.660512.11219.779
Table 2. Ratios of the average values of ω ¯ c to ω 1 and ω ˜ .
Table 2. Ratios of the average values of ω ¯ c to ω 1 and ω ˜ .
Model1234567
ω 1 / rad / s 0.06990.20470.44230.76611.08341.71292.4225
ω ¯ c / ω 1 49.94716.9846.72134.53923.54962.30071.6587
ω ¯ c / ω ˜ 0.70720.69970.68050.67850.73250.73980.7123
Model891011121314
ω 1 / rad / s 3.12744.09475.05126.25487.66068.84569.9986
ω ¯ c / ω 1 1.00001.00001.02411.00721.00341.00281.0039
ω ¯ c / ω ˜ 1.20051.39481.66462.06732.51532.90313.2833
Table 3. Peak frequency of response spectrum (rad/s).
Table 3. Peak frequency of response spectrum (rad/s).
Earthquake ω p Earthquake ω p Earthquake ω p
El Centro24.166San Fernando119.635Chuetsu-oki20.944
Taft14.280San Fernando217.453Victoria Mexico44.880
Wen Chuan52.359San Fernando314.960Loma Prieta10.134
Irpinia Italy34.906Corinth Greece12.083San Simeon28.560
Tottori Japan31.416Cape Mendocino18.580Artificial Wave15.708
Table 4. Ratios of the average values of ω ¯ c to ω 1 and ω p .
Table 4. Ratios of the average values of ω ¯ c to ω 1 and ω p .
Model12345
ω 1 / rad / s 1.21122.42253.83035.05126.2548
ω ¯ c 1 / ω 1 3.2282.8471.8361.3141.186
ω ¯ c 2 / ω p 0.1840.3500.3420.2950.351
Model678910
ω 1 / rad / s 7.66059.889711.419713.986215.3211
ω ¯ c 1 / ω 1 0.9340.9740.8950.7850.848
ω ¯ c 2 / ω p 0.3610.4780.4960.5690.671
Table 5. The first 10 natural frequencies and MPFs of the system.
Table 5. The first 10 natural frequencies and MPFs of the system.
ModeFrequency (rad/s)MPF (%)ModeFrequency (rad/s)MPF (%)
13.110436.804634.0493.3133
29.328112.260740.1612.7926
315.5377.3465846.2352.4091
421.7305.2373952.2642.1144
527.9034.06291058.2431.8805
Table 6. Displacement errors for m 50 ( ω ˜ = 1 rad / s , %).
Table 6. Displacement errors for m 50 ( ω ˜ = 1 rad / s , %).
Error IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω p ) TL- φ Chang/CR/TL
NEE0.00570.00720.00480.0074
NRMSE0.00910.02490.00140.0273
Table 7. Displacement errors for m 50 ( ω ˜ = 6 rad / s , %).
Table 7. Displacement errors for m 50 ( ω ˜ = 6 rad / s , %).
Error IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω p ) TL- φ Chang/CR/TL
NEE0.01200.08050.00460.0461
NRMSE0.10620.19780.00570.1329
Table 8. Displacement errors for m 50 ( ω ˜ = 16 rad / s , %).
Table 8. Displacement errors for m 50 ( ω ˜ = 16 rad / s , %).
Error IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω 3 ) TL- φ ( ω c = ω p ) TL- φ Chang/CR/TL
NEE2.70930.29050.47010.000012.8479
NRMSE3.16141.01271.09300.04143.2922
Table 9. Displacement errors for m 50 ( ω ˜ = 40 rad / s , %).
Table 9. Displacement errors for m 50 ( ω ˜ = 40 rad / s , %).
Error IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω 7 ) TL- φ ( ω c = ω p ) TL- φ Chang/CR/TL
NEE63.2245.40026.65194.628463.793
NRMSE24.3537.01597.03511.667924.356
Table 10. Natural frequencies and modal shapes of the building (G1).
Table 10. Natural frequencies and modal shapes of the building (G1).
ModeFrequency/rad/s Φ
ϕ 1 ϕ 2 ϕ 3 ϕ 4 ϕ 5
12.84630.28460.76351.00000.91900.5462
28.30830.54621.00000.2846−0.7636−0.9190
313.0970.76350.5462−0.9190−0.28461.0000
416.8250.9190−0.2846−0.54621.0000−0.7635
519.1901.0000−0.91900.7635−0.54620.2846
Table 11. Natural frequencies and modal shapes of the building (G2).
Table 11. Natural frequencies and modal shapes of the building (G2).
ModeFrequency/rad/s Φ
ϕ 1 ϕ 2 ϕ 3 ϕ 4 ϕ 5
10.98980.9700−0.0165−0.00390.00000.0000
26.19550.98970.60141.00000.00000.0000
316.1420.99970.9885−0.6019−0.0165−0.0038
462.4010.99990.9962−0.63540.61061.0000
5161.891.00001.0000−0.65241.0000−0.6170
Table 12. Displacement errors of the 5th floor (G1).
Table 12. Displacement errors of the 5th floor (G1).
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω 2 ) TL- φ TL
ϕ 2 NEE0.39150.15590.15590.4186
NRMSE1.41770.63550.63551.6489
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω 3 ) TL- φ TL
ϕ 3 NEE0.51040.16250.16250.5247
NRMSE4.53270.82580.82584.7661
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω 5 ) TL- φ TL
ϕ 5 NEE0.46860.19690.19690.4709
NRMSE7.45681.05061.05067.4800
Table 13. Displacement errors of the 5th floor (G2).
Table 13. Displacement errors of the 5th floor (G2).
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω ˜ c ) TL- φ ( ω c = ω 2 ) TL- φ TL
ϕ 2 NEE0.34150.33820.14070.11630.3470
NRMSE0.89220.88420.75900.73580.9059
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω ˜ c ) TL- φ ( ω c = ω 3 ) TL- φ TL
ϕ 3 NEE0.32300.32020.60180.11380.3276
NRMSE1.13311.12814.18630.91021.1418
RError IndexTL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω ˜ c ) TL- φ ( ω c = ω 5 ) TL- φ TL
ϕ 5 NEE0.49340.3865109.340.21950.5015
NRMSE1.04230.772427.2000.84701.0666
Table 14. Displacement errors for the 4th story.
Table 14. Displacement errors for the 4th story.
Error IndexTL- φ ( ω c = 8.14  rad/s)TL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω p ) TL- φ
NEE0.71471.160616.0391.7891
NRMSE0.59290.65824.16640.6330
Table 15. Displacement errors for the 4th story.
Table 15. Displacement errors for the 4th story.
Error IndexTL- φ ( ω c = 2.28 rad/s)TL- φ ( ω c = ω 1 ) TL- φ ( ω c = ω p ) TL- φ
NEE4.91784.91934.98085.0293
NRMSE0.82770.83120.95010.8514
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, Y.; Luo, C.; Fu, B. Choices of the Critical Frequency for φ in TL-φ Algorithms When Applied to Multi-Degree of Freedom Systems. Buildings 2022, 12, 863. https://doi.org/10.3390/buildings12060863

AMA Style

Tang Y, Luo C, Fu B. Choices of the Critical Frequency for φ in TL-φ Algorithms When Applied to Multi-Degree of Freedom Systems. Buildings. 2022; 12(6):863. https://doi.org/10.3390/buildings12060863

Chicago/Turabian Style

Tang, Yu, Chao Luo, and Bo Fu. 2022. "Choices of the Critical Frequency for φ in TL-φ Algorithms When Applied to Multi-Degree of Freedom Systems" Buildings 12, no. 6: 863. https://doi.org/10.3390/buildings12060863

APA Style

Tang, Y., Luo, C., & Fu, B. (2022). Choices of the Critical Frequency for φ in TL-φ Algorithms When Applied to Multi-Degree of Freedom Systems. Buildings, 12(6), 863. https://doi.org/10.3390/buildings12060863

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