Analytical exploration of γ-function explicit method for pseudodynamic testing of nonlinear systems

It has been well studied that the γ-function explicit method can be effective in providing favorable numerical dissipation for linear elastic systems. However, its performance for nonlinear systems is unclear due to a lack of analytical evaluation techniques. Thus, a novel technique is proposed herein to evaluate its efficiency for application to nonlinear systems by introducing two parameters to describe the stiffness change. As a result, the numerical properties and error propagation characteristics of the γ-function explicit method for the pseudodynamic testing of a nonlinear system are analytically assessed. It is found that the upper stability limit decreases as the step degree of nonlinearity increases; and it increases as the current degree of nonlinearity increases. It is also shown that this integration method provides favorable numerical dissipation not only for linear elastic systems but also for nonlinear systems. Furthermore, error propagation analysis reveals that the numerical dissipation can effectively suppress the severe error propagation of high frequency modes while the low frequency responses are almost unaffected for both linear elastic and nonlinear systems.     

Retracted: Discussion of paper ‘Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation’ by Chinmoy Kolay and James M. Ricles,Earthquake Engineering and Structural Dynamics 2014; 43:1361–1380

It seems that the explicit KR‐α method (KRM) is promising for the step‐by‐step integration because it simultaneously integrates unconditional stability, explicit formulation, and numerical dissipation together. It was shown that KRM can inherit the numerical dispersion and energy dissipation properties of the generalized‐α method (GM) for a linear elastic system, and it reduces to CR method (CRM) if ρ_∞ = 1is adopted, where ρ_∞ is the spectral radius of the amplification matrix of KRM as the product of the natural frequency and the step size tends to infinity. However, two unusual properties were found for KRM and CRM, and they might limit their application to solve either linear elastic or nonlinear systems. One is the lack of capability to capture the structural nonlinearity, and the other is that it is unable to realistically reflect the dynamic loading. Copyright © 2014 John Wiley & Sons, Ltd.     

Discussion of paper ‘Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation’ by Chinmoy Kolay and James M. Ricles,Earthquake Engineering and Structural Dynamics 2014; 43:1361–1380

It seems that the explicit KR‐α method (KRM), which was developed by Kolay and Ricles, is promising for the step‐by‐step integration because it simultaneously integrates unconditional stability, explicit formulation, and numerical dissipation together. It was shown that KRM can inherit the numerical dispersion and energy dissipation properties of the generalized‐α method [1] for a linear elastic system, and it reduces to CR method (CRM), which was developed by Chen and Ricles [2] if ρ_∞ = 1 is adopted, where ρ_∞ is the spectral radius of the amplification matrix of KRM as the product of the natural frequency and the step size tends to infinity. However, two unusual properties were found for KRM and CRM, and they might limit their application to solve either linear elastic or nonlinear systems. One is the lack of capability to capture the structural nonlinearity, and the other is that it is unable to realistically reflect the dynamic loading. Copyright © 2014 John Wiley & Sons, Ltd.     

Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation

The implicit dissipative generalized‐ α method is analyzed using discrete control theory. Based on this analysis, a one‐parameter family of explicit direct integration algorithms with controllable numerical energy dissipation, referred to as the explicit KR‐α method, is developed for linear and nonlinear structural dynamic numerical analysis applications. Stability, numerical dispersion, and energy dissipation characteristics of the proposed algorithms are studied. It is shown that the algorithms are unconditionally stable for linear elastic and stiffness softening‐type nonlinear systems, where the latter indicates a reduction in post yield stiffness in the force–deformation response. The amount of numerical damping is controlled by a single parameter, which provides a measure of the numerical energy dissipation at higher frequencies. Thus, for a specific value of this parameter, the resulting algorithm is shown to produce no numerical energy dissipation. Furthermore, it is shown that the influence of the numerical damping on the lower mode response is negligible. It is further shown that the numerical dispersion and energy dissipation characteristics of the proposed explicit algorithms are the same as that of the implicit generalized‐ α method. A numerical example is presented to demonstrate the potential of the proposed algorithms in reducing participation of undesired higher modes by using numerical energy dissipation to damp out these modes. Copyright © 2014 John Wiley & Sons, Ltd.     

Extensions of nonlinear error propagation analysis for explicit pseudodynamic testing

Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.     

Selection of time step for pseudodynamic testing

Although the step degree of nonlinearity has been introduced to conduct basic analysis and error propagation analysis for the pseudodynamic testing of nonlinear systems, it cannot be reliably used to select an appropriate time step before performing a pseudodynamic test. Therefore, a novel parameter of instantaneous degree of nonlinearity is introduced to monitor the stiffness change at the end of a time step, and can thus be used to evaluate numerical and error propagation properties for nonlinear systems. Based on these properties, it is possible to select an appropriate time step to conduct a pseudodynamic test in advance. This possibility is very important in pseudodynamic testing, since the use of an arbitrary time step might lead to unreliable results or even destroy the test specimen. In this paper, guidelines are proposed for choosing an appropriate time step for accurate integration of nonlinear systems. These guidelines require estimation of the maximum instantaneous degree of nonlinearity and the solution of the initial natural frequency. The Newmark explicit method is chosen for this study. All the analytical results and the guidelines proposed are thoroughly confirmed with numerical examples.     

A numerical study of 1-D nonlinear P-wave propagation in solid

IntroductionBecauseoftheextensivedistributionofruptures,micro-cracksandcrystalfracturesintheearth,therelationshipsbetweenthestressandstrainarenolongerlinear,infact,theyarenonlinear.Inordertoinvestigateandusethenonlinearcharacteristicsofsolidmediumintheearth,weshouldconsidertheinfluenceofnonlinearresponseduringthecourseofseismicmodelingandinversion.Thisisoneoftheimportantstudyfieldsthathavebeenpaidgreatattentionsinthere-centyearsintheworld(Minster,etal,1991;ZHANG,TENG,1993).Thenonlinearchar…     

Investigations of nonlinear performances of implicit pseudodynamic algorithms

Although it has been shown that the implementation of the HHT-α method can result in improved error propagation properties in pseudodynamic testing if the equation of motion is used instead of the difference equation to evaluate the next step acceleration, this paper proves that this method might lead to instability when used to solve a nonlinear system. Its unconditional stability is verified only for linear elastic systems, while for nonlinear systems, instability occurs as the step degree of convergence is less than 1. It is worth noting that the step degree of convergence can frequently be less than 1 in pseudodynamic testing, since a convergent solution is achieved only when the step degree of convergence is close to 1 regardless of whether its value is greater or less than 1. Therefore, the application of this scheme to pseudodynamic testing should be prohibited, since the possibility of instability might incorrectly destroy a specimen. Consequently, the implementation of the HHT-α method by using the difference equation to determine the next step acceleration is recommended for use in pseudodynamic testing.     

Nonlinear evaluations of unconditionally stable explicit algorithms

Two explicit integration algorithms with unconditional stability for linear elastic systems have been successfully developed for pseudodynamic testing.Their numerical properties in the solution of a linear elastic system have been well explored and their applications to the pseudodynamic testing of a nonlinear system have been shown to be feasible. However,their numerical properties in the solution of a nonlinear system are not apparent.Therefore,the performance of both algorithms for use in the solution...     

Dynamics of Rayleigh beam on nonlinear foundation due to moving load using Adomian decomposition and coiflet expansion

The effective procedure, which enables to discuss the dynamic response of Rayleigh beam resting on nonlinear viscoelastic foundation subjected to moving load, is developed. By employing the Adomian decomposition method in conjunction with coiflet expansion, the approximate closed form solution has been derived and condition for the convergence of the decomposition series has been introduced. To evaluate the accuracy of the approximate solution a local error index is defined. The presented new complex method is simple and efficient. The parametric study is performed and the influence of nonlinearity, load velocity, load frequency and the radius of gyration on the wave propagation in beam is investigated. The numerical results show that for the supercritical case, the linear model is stiffer giving rise to small displacement of the beam at the load passage.