Nonlinear Transient Response and its Sensitivity Using Finite Elements in Time

Author(s):  
Rakesh K. Kapania ◽  
Sungho Park

Abstract The bilinear formulation proposed earlier by Peters and Izadpanah to develop finite elements in time to solve undamped linear systems, is extended (and found to be readily amenable) to develop time finite elements to obtain transient responses of both linear and nonlinear, and damped and undamped systems. The formulation is used in the h-, p- and hp-versions. The resulting linear and nonlinear algebraic equations are differentiated to obtain the sensitivity of the transient response with respect to various design parameters. The present developments were tested on a series of linear and nonlinear examples and were found to yield, when compared with results obtained using other methods, excellent results for both the transient response and its sensitivity to system parameters. Mostly, the results were obtained using the Legendre polynomials as basis functions, though, in some cases other orthogonal polynomials namely, the Hermite, the Chebyshev, and integrated Legendre polynomials were also employed (but to no great advantage). A key advantage of the time finite element method, and the one often overlooked in its past applications, is the ease in which the sensitivity of the transient response with respect to various system parameters can be obtained.

2020 ◽  
Vol 16 (1) ◽  
Author(s):  
Zhifu Cao ◽  
Qingguo Fei ◽  
Dong Jiang ◽  
Rui Zhu ◽  
Hui Jin

Abstract Sensitivity analysis of transient response plays a crucial role in structural dynamics optimization. In this paper, the sensitivity analysis method for calculating the first-order derivatives of the nonlinear transient response using the real-imaginary perturbations (RIP) in the complex domain is proposed. Independent design parameters are synchronously perturbed using real and imaginary perturbations, respectively. The complex variable finite element method is employed to compute the complex transient response of the perturbed model. The nonlinear transient response sensitivities with respect to each parameter are obtained by separately extracting the real and imaginary responses. The computational accuracy and efficiency of the proposed method are demonstrated by employing a nonlinear multidegree-of-freedom system and a cantilever beam with nonlinear elastic supports. Results show that the response sensitivity with respect to multiple design parameters is obtained accurately by using the proposed method, the computational efficiency of which is increased compared with the complex variable method (CVM). The stability of the transient response sensitivity is significantly affected by the parameter perturbations in the real and imaginary parts. The RIP-imaginary-based response sensitivity is feasible even for an extremely small perturbation in the imaginary part; the RIP-real-based response sensitivity is stable only within a limited real perturbation range.


2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
E. H. Doha ◽  
D. Baleanu ◽  
A. H. Bhrawy ◽  
R. M. Hafez

A new Legendre rational pseudospectral scheme is proposed and developed for solving numerically systems of linear and nonlinear multipantograph equations on a semi-infinite interval. A Legendre rational collocation method based on Legendre rational-Gauss quadrature points is utilized to reduce the solution of such systems to systems of linear and nonlinear algebraic equations. In addition, accurate approximations are achieved by selecting few Legendre rational-Gauss collocation points. The numerical results obtained by this method have been compared with various exact solutions in order to demonstrate the accuracy and efficiency of the proposed method. Indeed, for relatively limited nodes used, the absolute error in our numerical solutions is sufficiently small.


2010 ◽  
Vol 2010 ◽  
pp. 1-13 ◽  
Author(s):  
Abdelouahab Kadem ◽  
Adem Kilicman

We consider the combined Walsh function for the three-dimensional case. A method for the solution of the neutron transport equation in three-dimensional case by using the Walsh function, Chebyshev polynomials, and the Legendre polynomials are considered. We also present Tau method, and it was proved that it is a good approximate to exact solutions. This method is based on expansion of the angular flux in a truncated series of Walsh function in the angular variable. The main characteristic of this technique is that it reduces the problems to those of solving a system of algebraic equations; thus, it is greatly simplifying the problem.


2017 ◽  
Vol 63 (3) ◽  
pp. 418-436
Author(s):  
Stefano Bianchini ◽  
Paolo Bonicatto ◽  
Elio Marconi

In this note we present a unifying approach for two classes of first order partial differential equations: we introduce the notion of Lagrangian representation in the settings of continuity equation and scalar conservation laws. This yields, on the one hand, the uniqueness of weak solutions to transport equation driven by a two dimensional BV nearly incompressible vector field. On the other hand, it is proved that the entropy dissipation measure for scalar conservation laws in one space dimension is concentrated on countably many Lipschitz curves.


Author(s):  
I. V. Boikov ◽  
A. I. Boikova

Continuous Seidel method for solving systems of linear and nonlinear algebraic equations is constructed in the article, and the convergence of this method is investigated. According to the method discussed, solving a system of algebraic equations is reduced to solving systems of ordinary differential equations with delay. This allows to use rich arsenal of numerical ODE solution methods while solving systems of algebraic equations. The main advantage of the continuous analogue of the Seidel method compared to the classical one is that it does not require all the elements of the diagonal matrix to be non-zero while solving linear algebraic equations’ systems. The continuous analogue has the similar advantage when solving systems of nonlinear equations.


2009 ◽  
Vol 20 (2) ◽  
pp. 187-214 ◽  
Author(s):  
WAN CHEN ◽  
MICHAEL J. WARD

The dynamics and oscillatory instabilities of multi-spike solutions to the one-dimensional Gray-Scott reaction–diffusion system on a finite domain are studied in a particular parameter regime. In this parameter regime, a formal singular perturbation method is used to derive a novel ODE–PDE Stefan problem, which determines the dynamics of a collection of spikes for a multi-spike pattern. This Stefan problem has moving Dirac source terms concentrated at the spike locations. For a certain subrange of the parameters, this Stefan problem is quasi-steady and an explicit set of differential-algebraic equations characterizing the spike dynamics is derived. By analysing a nonlocal eigenvalue problem, it is found that this multi-spike quasi-equilibrium solution can undergo a Hopf bifurcation leading to oscillations in the spike amplitudes on an O(1) time scale. In another subrange of the parameters, the spike motion is not quasi-steady and the full Stefan problem is solved numerically by using an appropriate discretization of the Dirac source terms. These numerical computations, together with a linearization of the Stefan problem, show that the spike layers can undergo a drift instability arising from a Hopf bifurcation. This instability leads to a time-dependent oscillatory behaviour in the spike locations.


Sign in / Sign up

Export Citation Format

Share Document