Improvement of the Linear Interpolation Method of Excitation Using the Spline Interpolation Function to Numerically Calculate the Response Coefficient of Linear One-Degree-Of-Freedom Systems

Document Type : Research Note

Authors

1 Assistant Professor, Department of Civil Engineering, Faculty of Engineering, University of Qom, Qom, Iran

2 Ph.D. Candidate, Department of Civil Engineering, Faculty of Engineering, University of Qom, Qom, Iran

3 Ph.D. Candidate, Department of Civil Engineering, Noshirvani University of Technology, Babol, Iran

4 M.Sc. Student, Department of Civil Engineering, Faculty of Engineering, University of Qom, Qom, Iran

Abstract

In this article, the Jennings method has been improved by using the third-order spline interpolation function. In this research, in order to be able to compare the advantages and disadvantages of the Jennings method, which is based on exact relationships and the assumption of linear changes of the excitation, compared to the spline interpolation method, a damped one-degree-of-freedom linear system under sinusoidal harmonic loading is considered has been taken. Exact values of the deformation response coefficient of this system are available for different excitation frequencies. The approximate values of the deformation response coefficient of this system for different excitation frequencies were calculated assuming linear interpolation of excitation and also assuming interpolation with spline function and were compared with their exact corresponding values. This work was done for two, five, ten and twenty percent damping. The results of the work indicated that when the number of points by which the sine wave is approximated is small and also the amount of damping is low, the interpolation with the spline function has a significantly higher accuracy than the linear interpolation mode. Another noteworthy point is the high sensitivity of the cubic spline interpolation to the level of system damping, so that for a certain number of divisions, the error value is highly sensitive to the damping value. Thus, with the increase of damping, the sign of the error changes and its value increases strongly for the number of large divisions.  This phenomenon is not seen at all in linear interpolation.
Therefore, considering that in interpolation using the spline method, the continuity of the slope and the second derivative are maintained in the internal points, it is suggested to use this method to calculate the dynamic response of linear systems and the results obtained with the results of Jennings' method should be compared and appropriate decisions should be made if there is a significant difference in the analysis results obtained from these two methods.
It is important to mention that the execution time of the computer program related to the cubic spline interpolation method is longer than the Jennings method. Because in the spline interpolation method, it is necessary to solve the system of linear equations with the number of unknowns equal to four times the number of intervals, while in the excitation linear interpolation method, it is not necessary to solve the system of equations at all and simply by having the velocity and displacement at the beginning of the step and assuming that the excitation is linear during the desired time step, the coefficients A, B, C, D, A', B', C', D' were calculated and using this eight coefficient, velocity and displacement at the end of the time step are obtained. It should be noted that if the time interval between different accelerogram points is constant, it is sufficient that the coefficients A, B, C, D, A', B', C', D' are calculated only once.

Keywords

Main Subjects


Burden, R.L., & Faires, J.D. (2011). Numerical Analysis. Ninth Edition, Books/Cole.
Chapra, S.C., & Canale, R.P. (2006). Numerical Methods for Engineers. Fifth Edition, McGraw-Hill, New York.
Chopra, A.K. (2012). Dynamics of Structures, Theory and Applications to Earthquake Engineering. Fourth Edition, Prentice-Hall.
Esmailabadi, R., Bahar, A., & Azimi-Nejad, A. (2015). Presentation of new relationships of hysteresis   damping capacity of steel bending frames at the level of life safety performance required by direct design based on displacement. Scientific Research Quarterly of Earthquake Science and Engineering, 3(4), 43-59 (in Persian).
Ghazanfari, S., Hamzehei-Javaran, S., Alesadi A., & Shojaee, S. (2019). Free vibration analysis of cross-ply laminated beam structures using refined beam theories and B-Spline Basis Functions. Mechanics of Advanced Materials and Structures, 467-475.
Hanselman, D., & Littlefield B. (1996). Mastering MATLAB, A Comprehensive Tutorial and Reference. First Edition, Prentice-Hall.
Hoseini Hashemi, B., Hoseini, M., & Khanlari, K. (2002). Effect of P-delta on dynamic analysis of
structures, second order analysis of structures under dynamic loads. Journal of Seismology and
Earthquake Engineering, 4(4), 2-10 (in Persian).
Izadinia, M., & Jamshidi, J. (2016). Controlling the dynamic vibration of adjacent structures by
establishing a connection. Second International on Geotechnics and Seismic Engineering, Tabriz,
Iran (in Persian).
Liu, Z.W., Chen R.S., & Chen J.Q. (2008). Adaptive sampling cubic-spline interpolation method for efficient calculation of monostatic RCS. Microwave and Optical Technology Letters, 50(3), 751-755.
Mahdavi, S.H., Razak, H.A., Shojahee, S., & Mahdavi, M.S. (2015). A comparative study on application of chebyshev and spline methods for geometrically non-linear analysis of truss structures. International Journal of Mechanical Sciences, 101-102, 241-251.
Maleknejad, K., & Derili, H. (2006). Numerical solution of integral equatins by using combination of spline-collocation method and Lagrange interpolation. Applied Mathematics and Computation, 175(2), 1235-1244.
Maleknejad, K., & Derili, H. (2007). Numerical solution of Hammerstein integral equatins by using combination of spline-collocation method and Lagrange interpolation. Applied Mathematics and Computation, 190(2), 1557-1562.
Mohammadi Nia, M., Shojaee, S., & Hamzehei-Javaran, S. (2020). A mixed formulation of b-spline and a new class of spherical hankel shape functions for modeling elastostatic problems. Applied Mathematical Modelling, 77, 602-616.
Motaref, S. (2006). Probabilistic Study of the Influence of Column-to-Beam Resistance Ratio Parameter on    the Seismic Performance of Concrete Frames. Civil Engineering-Earthquake Engineering M.Sc. Thesis, International Institute of Earthquake Engineering and Seismology, Tehran (in Persian).
Naeim, F. (1996) Respnse of Instrumented Buildings    to 1994 Northridge Earthquake. Draft Report     CSMIP.
Nigam, N.C. and Jennings, P.C. (1968). Digital Calculation of Response Spectra from Strong-Motion Records. California Institute of Technology, Earthquake Engineering Research Laboratory.
Nigam, N.C., & Jennings, P.C. (1969). Calculation of response spectra from strong-motion earthquake records. Bulletin of the Seismological Socity of America, 59(2), 909-922.
Rostami, S., & Shojaee, S. (2017). A family of cubic b-spline direct integration algorithms with contorllabe numerical dissipation and dispersion for structural dynamics. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 42, 17-32.
Rostami, S., & Shojaee, S. (2019). Development of       a direct integration method on quartic b-spline collocation method. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 43, 615-636.
Saffari, H., Shojaee, S., Rostami, S., & Malekinejad, M. (2014). Application of cubic spline on large deformation analysis of structures. International Journal of Steel Structures, 14(1), 165-172.
 
Shafiee-fard, M.R. (2019). Second-Order Analysis of Structures under Dynamic Loads by Jennings Numerical Method, Publisher: Moalefine Talaee (in Persian).
Shahmorad, S., & Abdollahi, A. (2007). A Quadrature free convergent method for the numerical solution       of linear fredholm integral equations based on    hermite-spline interpolation. Proceeding in Applied Mathematics and Mechanics, 7(1), 41-42.
Shojaee, S., Rostami, S., & Abbasi, A. (2015). An unconditionally stable implicit time integration algorithm: modified quartic b-spline method. Computers and Structures, 153, 98-111.
Vamvatsikos, D., & Cornell, C.A. (2004). Applied incremental dynamic analysis. EarthquakeSpectra, 20(2), 523-553.
Yu, R., Wang, R. & Zhu, C. (2013). A numerical method for solving KdV equation with multilevel bspline quasi-interpolation. Applicable Analysis, 92(8), 1682-1690.
Zhernakov, V.S., Pavlov, V.P., & Kudoyarova, V.M. (2017). The enhanced spline-method for numerical results of natural frequencies of beams. Procedia Engineering, 176, 438-450.