Hostname: page-component-78c5997874-s2hrs Total loading time: 0 Render date: 2024-11-19T20:38:28.314Z Has data issue: false hasContentIssue false
  • English
  • Français

A Matrix-Vector Operation-Based Numerical Solution Method for Linear m-th Order Ordinary Differential Equations: Application to Engineering Problems

Published online by Cambridge University Press:  03 June 2015

M. Aminbaghai ,
M. Dorn ,
J. Eberhardsteiner  and
B. Pichler
M. Aminbaghai*
Affiliation:
Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), Karlsplatz 13/202, A-1040 Vienna, Austria
M. Dorn*
Affiliation:
Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), Karlsplatz 13/202, A-1040 Vienna, Austria Linnaeus University, Department of Building and Energy Technology, Lückligs Plats 1, S-35195 Växjö, Sweden
J. Eberhardsteiner*
Affiliation:
Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), Karlsplatz 13/202, A-1040 Vienna, Austria
B. Pichler*
Affiliation:
Institute for Mechanics of Materials and Structures, Vienna University of Technology (TU Wien), Karlsplatz 13/202, A-1040 Vienna, Austria
*
Email: mehdi.aminbaghai@tuwien.ac.at
Email: michael.dorn@tuwien.ac.at
Email: josef.eberhardsteiner@tuwien.ac.at
Corresponding author.Email: bernhard.pichler@tuwien.ac.at
Get access

Abstract

Many problems in engineering sciences can be described by linear, inhomogeneous, m-th order ordinary differential equations (ODEs) with variable coefficients. For this wide class of problems, we here present a new, simple, flexible, and robust solution method, based on piecewise exact integration of local approximation polynomials as well as on averaging local integrals. The method is designed for modern mathematical software providing efficient environments for numerical matrix-vector operation-based calculus. Based on cubic approximation polynomials, the presented method can be expected to perform (i) similar to the Runge-Kutta method, when applied to stiff initial value problems, and (ii) significantly better than the finite difference method, when applied to boundary value problems. Therefore, we use the presented method for the analysis of engineering problems including the oscillation of a modulated torsional spring pendulum, steady-state heat transfer through a cooling web, and the structural analysis of a slender tower based on second-order beam theory. Related convergence studies provide insight into the satisfying characteristics of the proposed solution scheme.

Keywords

34A3034B60Numerical integrationpolynomial approximationODEvariable coefficientsinitial conditionsboundary conditionsstiff equation
Type
Research Article
Information
Advances in Applied Mathematics and Mechanics , Volume 5 , Issue 3 , June 2013 , pp. 269 - 308
Copyright
Copyright © Global-Science Press 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Hairer, E., Nørsett, S. and Wanner, G., Solving Ordinary Differential Equations I: Non-stiff Problems, 2nd edn., Springer Verlag: Berlin, 1993.Google Scholar
[2]Hairer, E. and Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2nd edn., Springer Verlag: Berlin, New York, 1996.Google Scholar
[3]Cohen, SD. and Hindmarsh, AC., CVODE, a stiff/nonstiff ODE solver in C, Comput. Phys., 10(2) (1996), pp.138–143.Google Scholar
[4]Ascher, U. and Petzold, L., Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, Society for Industrial and Applied Mathematics (SIAM): Philadelphia, 1998.CrossRefGoogle Scholar
[5]Rosenbrock, H., Some general implicit processes for the numerical solution of differential equations, The Computer J., 5(4) (1963), pp.329–330.Google Scholar
[6]Sandu, A., Verwer, JG., Blom, JG., Spee, EJ., Carmichael, GR. and Potra, FA., Benchmarking stiff ODE solvers for atmospheric chemistry problems II: Rosenbrock solvers, Atmospheric Environment, 31(20) (1997), pp.3459–3472.CrossRefGoogle Scholar
[7]Shampine, LF. and Reichelt, MW., The MATLAB ODE suite, SIAM J. Sci. Comput., 18(1) (1997), pp.1–22.Google Scholar
[8]Atluri, SN., Cho, JY. and Kim, HG.Analysis of thin beams, using the meshless local Petrov-Galerkin method, with generalized moving least squares interpolations, Comput. Mech., 24(5) (1999), pp.334–347.CrossRefGoogle Scholar
[9]Gu, YT. and Liu, GR., A local point interpolation method for static and dynamic analysis of thin beams, Comput. Methods Appl. Mech. Eng., 190(42) (2001), pp.5515–5528.Google Scholar
[10]Gudla, PK. and Ganguli, R., Discontinuous Galerkin finite element in time for solving periodic differential equations, Comput. Methods Appl. Mech. Eng., 196(1-3) (2006), pp.682–696.Google Scholar
[11]Ramos, JI., Piecewise-linearized methods for initial-value problems with oscillating solutions, Appl. Math. Comput., 181(1) (2006), pp.123–146.Google Scholar
[12]Fung, TC., Solving initial value problems by differential quadrature method-part 1: first-order equations, Int. J. Numer. Methods Eng., 50(6) (2001), pp.1411–1427.Google Scholar
[13]Fung, TC., Solving initial value problems by differential quadrature method-part 2: second- and higher-order equations, Int. J. Numer. Methods Eng., 50(6) (2001), pp.1429–1454.3.0.CO;2-A>CrossRefGoogle Scholar
[14]Nie, GJ. and Z, Z. Zhong, Semi-analytical solution for three-dimensional vibration of functionally graded circular plates, Comput. Methods Appl. Mech. Eng., 196(49-52) (2007), pp.: 4901–4910.Google Scholar
[15]Hon, Y. and Wu, Z., A quasi-interpolation method for solving stiff ordinary differential equations, Int. J. Numer. Methods Eng., 48(8) (2000), pp.1187–1197.Google Scholar
[16]Mai-Duy, N., Solving high order ordinary differential equations with radial basis function networks, Int. J. Numer. Methods Eng., 62(6) (2005), pp.824–852.CrossRefGoogle Scholar
[17]Rubin, H. and Aminbaghai, M., Woölbkrafttorsion bei veränderlichem, offenem Querschnitt – Hat die Biegezugstab-Analogie noch Guültigkeit? [Warping torsion for variable, open cross sections – Is the analogy for theory of beams with tensile force still valid?], Stahlbau, 76(10) (2007), pp.747–760. In German.Google Scholar
[18]Fung, TC., Stability and accuracy of differential quadrature method in solving dynamic problems, Comput. Methods Appl. Mech. Eng., 191(13-14) (2002), pp.1311–1331.Google Scholar
[19]Chang, J., Yang, Q. and Liu, C., B-spline method for solving boundary value problems of linear ordinary differential equations, Information Computing and Applications, Communications in Computer and Information Science, 106 (2010), Zhu, R., Zhang, Y., Liu, B. and Liu, C. (eds.), Springer Berlin Heidelberg, pp.326–333.CrossRefGoogle Scholar
[20]Butikov, EI., Parametric excitation of a linear oscillator, Euro. J. Phys., 25(4) (2004), pp.535–554.Google Scholar
[21]Kraus, A. and Bar-Cohen, A., Thermal Analysis and Control of Electronic Equipment, Hemisphere Publishing Corporation: Washington, 1983.Google Scholar
[22]Rubin, H., Lösung li nearer Differentialgleichungen beliebiger Ordnung mit Polynomkoeffizienten und Anwendung auf ein baustatisches Problem [Solution of linear differential equations of arbitrary order with polynomial coefficients and application to a structural analysis problem], ZAMM Zeitschrift für Angewandte Mathematik und Mechanik, 76(2) (1996), pp.105–117. In German.Google Scholar