Using a small algebraic manipulation system to solve differential and integral equations by variational and approximation techniques

The microcomputer algebraic manipulation system MUMATH is used to implement the classical variational, Galerkin and least-squares techniques for solving boundary-value problems in differential equations and also for solving Fredholm integral equations. Examples are given which extend the precision of known results. The technique is presented as a general algorithm which can readily be implemented on other algebraic manipulation systems

Nodally exact Ritz discretizations of 1D diffusion–absorption and Helmholtz equations by variational FIC and modified equation methods

The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-005-0011-zThis article presents the first application of the Finite Calculus (FIC) in a Ritz-FEM variational framework. FIC provides a steplength parametrization of mesh dimensions, which is used to modify the shape functions. This approach is applied to the FEM discretization of the steady-state, one-dimensional, diffusion–absorption and Helmholtz equations. Parametrized linear shape functions are directly inserted into a FIC functional. The resulting Ritz-FIC equations are symmetric and carry a element-level free parameter coming from the function modification process. Both constant- and variable-coefficient cases are studied. It is shown that the parameter can be used to produce nodally exact solutions for the constant coefficient case. The optimal value is found by matching the finite-order modified differential equation (FOMoDE) of the Ritz-FIC equations with the original field equation. The inclusion of the Ritz-FIC models in the context of templates is examined. This inclusion shows that there is an infinite number of nodally exact models for the constant coefficient case. The ingredients of these methods (FIC, Ritz, MoDE and templates) can be extended to multiple dimensions.Peer ReviewedPostprint (author's final draft

Nodally exact Ritz discretizations of 1D diffusion–absorption and Helmholtz equations by variational FIC and modified equation methods

This article presents the first application of the Finite Calculus (FIC) in a Ritz-FEM variational framework. FIC provides a steplength parametrization of mesh dimensions, which is used to modify the shape functions. This approach is applied to the FEM discretization of the steady-state, one-dimensional, diffusion–absorption and Helmholtz equations. Parametrized linear shape functions are directly inserted into a FIC functional. The resulting Ritz-FIC equations are symmetric and carry a element-level free parameter coming from the function modification process. Both constant- and variable-coefficient cases are studied. It is shown that the parameter can be used to produce nodally exact solutions for the constant coefficient case. The optimal value is found by matching the finite-order modified differential equation (FOMoDE) of the Ritz-FIC equations with the original field equation. The inclusion of the Ritz-FIC models in the context of templates is examined. This inclusion shows that there is an infinite number of nodally exact models for the constant coefficient case. The ingredients of these methods (FIC, Ritz, MoDE and templates) can be extended to multiple dimension

Some results on thermal stress of layered plates and shells by using Unified Formulation

This work presents some results on two-dimensional modelling of thermal stress problems in multilayered structures. The governing equations are written by referring to the Unified Formulation (UF) introduced by the first author. These equations are obtained in a compact form, that doesn't depend on the order of expansion of variables in the thickness direction or the variable description (layer-wise models and equivalent single layers models). Classical and refined theories based on the Principle of Virtual Displacements (PVD) and advanced mixed theories based on the Reissner Mixed Variational Theorem (RMVT) are both considered. As a result, a large variety of theories are derived and compared. The temperature profile along the thickness of the plate/shell is calculated by solving the Fourier's heat conduction equation. Alternatively, thermo-mechanical coupling problems can be considered, in which the thermal variation is influenced by mechanical loading. Exact closed-form solutions are provided for plates and shells, but also the applications of the Ritz method and the Finite Element Method (FEM) are presented

INVESTIGATION OF THE DYNAMIC BEHAVIOUR OF NON-UNIFORM THICKNESS CIRCULAR PLATES RESTING ON WINKLER AND PASTERNAK FOUNDATIONS

The study of the dynamic behaviour of non-uniform thickness circular plate resting on elastic foundations is very imperative in designing structural systems. This present research investigates the free vibration analysis of varying density and non-uniform thickness isotropic circular plates resting on Winkler and Pasternak foundations. The governing differential equation is analysed using the Galerkin method of weighted residuals. Linear and nonlinear case is considered, the surface radial and circumferential stresses are also determined. Thereafter, the accuracy and consistency of the analytical solutions obtained are ascertained by comparing the existing results available in pieces of literature and confirmed to be in a good harmony. Also, it is observed that very accurate results can be obtained with few computations. Issues relating to the singularity of circular plate governing equations are handled with ease. The analytical solutions obtained are used to determine the influence of elastic foundations, homogeneity and thickness variation on the dynamic behaviour of the circular plate, the effect of vibration on a free surface of the foundation as well as the influence of radial and circumferential stress on mode shapes of the circular plate considered. From the results, it is observed that a maximum of 8.1% percentage difference is obtained with the solutions obtained from other analytical methods. Furthermore, increasing the elastic foundation parameter increases the natural frequency. Extrema modal displacement occurs due to radial and circumferential stress. Natural frequency increases as the thickness of the circular plate increases, Conversely, a decrease in natural frequency is observed as the density varies. It is envisioned that; the present study will contribute to the existing knowledge of the classical theory of vibration

Nodally exact Ritz discretizations of 1D diffusion–absorption and Helmholtz equations by variational FIC and modified equation methods

This article presents the first application of the Finite Calculus (FIC) in a Ritz-FEM variational framework. FIC provides a steplength parametrization of mesh dimensions, which is used to modify the shape functions. This approach is applied to the FEM discretization of the steady-state, one-dimensional, diffusion–absorption and Helmholtz equations. Parametrized linear shape functions are directly inserted into a FIC functional. The resulting Ritz-FIC equations are symmetric and carry a element-level free parameter coming from the function modification process. Both constant- and variable-coefficient cases are studied. It is shown that the parameter can be used to produce nodally exact solutions for the constant coefficient case. The optimal value is found by matching the finite-order modified differential equation (FOMoDE) of the Ritz-FIC equations with the original field equation. The inclusion of the Ritz-FIC models in the context of templates is examined. This inclusion shows that there is an infinite number of nodally exact models for the constant coefficient case. The ingredients of these methods (FIC, Ritz, MoDE and templates) can be extended to multiple dimension

A hybrid perturbation-Galerkin method for differential equations containing a parameter

A two-step hybrid perturbation-Galerkin method to solve a variety of differential equations which involve a parameter is presented and discussed. The method consists of: (1) the use of a perturbation method to determine the asymptotic expansion of the solution about one or more values of the parameter; and (2) the use of some of the perturbation coefficient functions as trial functions in the classical Bubnov-Galerkin method. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is illustrated first with a simple linear two-point boundary value problem and is then applied to a nonlinear two-point boundary value problem in lubrication theory. The results obtained from the hybrid method are compared with approximate solutions obtained by purely numerical methods. Some general features of the method, as well as some special tips for its implementation, are discussed. A survey of some current research application areas is presented and its degree of applicability to broader problem areas is discussed