A New Implicit Integration Scheme for the Simulation of Stiff-Complex Mechanical Systems

2005 ◽  
Author(s):  
Mojtaba Oghbaei ◽  
Kurt S. Anderson
Keyword(s):  
Mechanical Systems ◽  
Implicit Time ◽  
Integration Time ◽  
Integration Scheme ◽  
Algebraic Equations ◽  
Stiff System ◽  
Time Step ◽  
Efficient Manner ◽  
Method Performance ◽  
Integration Schemes

When performing the dynamic simulation of stiff mechanical systems, implicit type integration schemes are usually required to preserve stability. This article presents a new implicit time integrator, which is a particular application of a novel state-time formulation recently developed by the authors in a more general scope. The proposed scheme is constructed with the intent of benefitting from the accuracy and apparent robustness thus far achieved with this algorithm in an integration context. This is realized by first setting up the weighted residual form of the governing equations of the system in a form associated with the application of a time marching integration scheme. The resulting algebraic equations are then solved, minimizing the error of integration time step in a generalized energy sense, allowing one to capture the stiff behavior of solution in an efficient manner. Examples are provided to show the proposed method performance when dealing with a stiff system.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

scholarly journals A Semi-Implicit Version of the MPAS-Atmosphere Dynamical Core

2015 ◽  
Vol 143 (9) ◽  
pp. 3838-3855 ◽  
Author(s):  
Steven Sandbach ◽  
John Thuburn ◽  
Danail Vassilev ◽  
Michael G. Duda
Keyword(s):  
Time Integration ◽  
Standard Test ◽  
Newton Iteration ◽  
Atmospheric Modeling ◽  
Implicit Time ◽  
Integration Scheme ◽  
Baroclinic Wave ◽  
Time Step ◽  
Integration Schemes ◽  
Helmholtz Problem

Abstract An important question for atmospheric modeling is the viability of semi-implicit time integration schemes on massively parallel computing architectures. Semi-implicit schemes can provide increased stability and accuracy. However, they require the solution of an elliptic problem at each time step, creating concerns about their parallel efficiency and scalability. Here, a semi-implicit (SI) version of the Model for Prediction Across Scales (MPAS) is developed and compared with the original model version, which uses a split Runge–Kutta (SRK3) time integration scheme. The SI scheme is based on a quasi-Newton iteration toward a Crank–Nicolson scheme. Each Newton iteration requires the solution of a Helmholtz problem; here, the Helmholtz problem is derived, and its solution using a geometric multigrid method is described. On two standard test cases, a midlatitude baroclinic wave and a small-planet nonhydrostatic gravity wave, the SI and SRK3 versions produce almost identical results. On the baroclinic wave test, the SI version can use somewhat larger time steps (about 60%) than the SRK3 version before losing stability. The SI version costs 10%–20% more per step than the SRK3 version, and the weak and strong scalability characteristics of the two versions are very similar for the processor configurations the authors have been able to test (up to 1920 processors). Because of the spatial discretization of the pressure gradient in the lowest model layer, the SI version becomes unstable in the presence of realistic orography. Some further work will be needed to demonstrate the viability of the SI scheme in this case.

scholarly journals Efficient Approximate Techniques for Integrating Stochastic Differential Equations

2006 ◽  
Vol 134 (10) ◽  
pp. 3006-3014 ◽  
Author(s):  
James A. Hansen ◽  
Cecile Penland
Keyword(s):  
Time Series ◽  
Differential Equations ◽  
Stochastic Differential Equations ◽  
Deterministic Model ◽  
Probability Distributions ◽  
Approximate Solutions ◽  
Integration Scheme ◽  
Time Step ◽  
Integration Schemes ◽  
Runge Kutta

Abstract The delicate (and computationally expensive) nature of stochastic numerical modeling naturally leads one to look for efficient and/or convenient methods for integrating stochastic differential equations. Concomitantly, one may wish to sensibly add stochastic terms to an existing deterministic model without having to rewrite that model. In this note, two possibilities in the context of the fourth-order Runge–Kutta (RK4) integration scheme are examined. The first approach entails a hybrid of deterministic and stochastic integration schemes. In these examples, the hybrid RK4 generates time series with the correct climatological probability distributions. However, it is doubtful that the resulting time series are approximate solutions to the stochastic equations at every time step. The second approach uses the standard RK4 integration method modified by appropriately scaling stochastic terms. This is shown to be a special case of the general stochastic Runge–Kutta schemes considered by Ruemelin and has global convergence of order one. Thus, it gives excellent results for cases in which real noise with small but finite correlation time is approximated as white. This restriction on the type of problems to which the stochastic RK4 can be applied is strongly compensated by its computational efficiency.

Studies on the Accuracy Stability and Efficiency of a New Time Integration Scheme for Ballast Modeling Using the Discrete Element Method

2014 ◽  
Author(s):  
G. F. Mathews ◽  
R. L. Mullen ◽  
D. C. Rizos
Keyword(s):  
Particle Interaction ◽  
Time Integration ◽  
Discrete Element ◽  
Critical Discussion ◽  
Implicit Time ◽  
Integration Scheme ◽  
Computationally Efficient ◽  
Time Step ◽  
Central Difference ◽  
Time Integration Scheme

This paper presents the development of a semi-implicit time integration scheme, originally developed for structural dynamics in the 1970’s, and its implementation for use in Discrete Element Methods (DEM) for rigid particle interaction, and interaction of elastic bodies that are modeled as a cluster of rigid interconnected particles. The method is developed in view of ballast modeling that accounts for the flexibility of aggregates and the arbitrary shape and size of granules. The proposed scheme does not require any matrix inversions and is expressed in an incremental form making it appropriate for non-linear problems. The proposed method focuses on improving the efficiency, stability and accuracy of the solutions, as compared to current practice. A critical discussion of the findings of the studies is presented. Extended verification and assessment studies demonstrate that the proposed algorithm is unconditionally stable and accurate even for large time step sizes. It is demonstrated that the proposed method is at least as computationally efficient as the Central Difference Method. Guidelines for the implementation of the method to ballast modeling are discussed.

Modeling the damped dynamic behavior of a flexible pendulum

2019 ◽  
Vol 54 (2) ◽  
pp. 116-129 ◽  
Author(s):  
Roberto Ortega ◽  
Geraldine Farías ◽  
Marcela Cruchaga ◽  
Matías Rivero ◽  
Mariano Vázquez ◽  
...  
Keyword(s):  
High Speed ◽  
Time Integration ◽  
Integration Scheme ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Rayleigh Damping ◽  
Integration Schemes ◽  
Time Integration Schemes ◽  
Damping Model

The focus of this work is on the computational modeling of a pendulum made of a hyperelastic material and the corresponding experimental validation with the aim of contributing to the study of a material commonly used in seismic absorber devices. From the proposed dynamics experiment, the motion of the pendulum is recorded using a high-speed camera. The evolution of the pendulum’s positions is recovered using a capturing motion technique by tracking markers. The simulation of the problem is developed in the framework of a parallel multi-physics code. Particular emphasis is placed on the analysis of the Newmark integration scheme and the use of Rayleigh damping model. In particular, the time step size effect is analyzed. A strong time step size dependency is obtained for dissipative time integration schemes, while the Rayleigh damping formulation without time integration dissipation shows time step–independent results when convergence is achieved.

Impact Analysis of Sandwich Composites

2008 ◽  
Author(s):  
Laura Ferrero ◽  
Ugo Icardi
Keyword(s):  
Impact Analysis ◽  
Time Integration ◽  
Nonlinear Problems ◽  
Implicit Time ◽  
Contact Radius ◽  
Integration Scheme ◽  
Sandwich Composites ◽  
Time Step ◽  
Wide Range ◽  
The Impact

A finite element simulation of impacts on sandwich composites with laminated faces is presented; it is based on a refined multilayered plate model with a high-order zig-zag representation of displacements, which is incorporated through a strain energy updating process. This allows the implementation into existing commercial finite elements codes, preserving their program structure. As customary, the Hertzian law and the Newmark implicit time integration scheme are used for solving the contact problem. The contact radius and the force are computed within each time step by an iterative algorithm which forces the impacted top surface to conform, in the least-squares sense, to the shape of the impactor. Nonlinear strains of von Karman type are used. As appearing by the comparison with experimental results, the present model is able to accurately predict the impact force, the core damage and the damage of face sheets in sandwich composites with foam and or honeycomb core. Moreover, this paper also assesses the accuracy and the range of application of stress based criteria in predicting the onset and evolution of delamination in service. These criteria are widespread by virtue of their low run time and storage costs, although no exhaustive proofs are known weather they are accurate enough for a reasonably wide range of applications. Since where highly iterative solutions are involved (e.g., impact and geometric, or material nonlinear problems) they are the only currently affordable failure models, it appears of primary importance to fill this gap. Aimed to contribute to the knowledge advancement in this field, a comparison is presented with more sophisticate fracture mechanics and progressive delamination models.

scholarly journals Acceleration of Runge-Kutta integration schemes

2004 ◽  
Vol 2004 (2) ◽  
pp. 307-314 ◽  
Author(s):  
Phailaung Phohomsiri ◽  
Firdaus E. Udwadia
Keyword(s):  
Computational Cost ◽  
Fixed Time ◽  
Lower Number ◽  
Integration Scheme ◽  
Time Step ◽  
Third Order ◽  
Numerical Examples ◽  
Integration Schemes ◽  
Local Accuracy ◽  
Runge Kutta

A simple accelerated third-order Runge-Kutta-type, fixed time step, integration scheme that uses just two function evaluations per step is developed. Because of the lower number of function evaluations, the scheme proposed herein has a lower computational cost than the standard third-order Runge-Kutta scheme while maintaining the same order of local accuracy. Numerical examples illustrating the computational efficiency and accuracy are presented and the actual speedup when the accelerated algorithm is implemented is also provided.

Scaling of Constraints and Augmented Lagrangian Formulations in Multibody Dynamics Simulations

2009 ◽  
Vol 4 (2) ◽  
Author(s):  
Olivier A. Bauchau ◽  
Alexander Epple ◽  
Carlo L. Bottasso
Keyword(s):  
Physical Properties ◽  
Augmented Lagrangian ◽  
Equations Of Motion ◽  
Time Integration ◽  
Differential Algebraic Equations ◽  
Algebraic Equations ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Integration Schemes

This paper addresses practical issues associated with the numerical enforcement of constraints in flexible multibody systems, which are characterized by index-3 differential algebraic equations (DAEs). The need to scale the equations of motion is emphasized; in the proposed approach, they are scaled based on simple physical arguments, and an augmented Lagrangian term is added to the formulation. Time discretization followed by a linearization of the resulting equations leads to a Jacobian matrix that is independent of the time step size, h; hence, the condition number of the Jacobian and error propagation are both O(h0): the numerical solution of index-3 DAEs behaves as in the case of regular ordinary differential equations (ODEs). Since the scaling factor depends on the physical properties of the system, the proposed scaling decreases the dependency of this Jacobian on physical properties, further improving the numerical conditioning of the resulting linearized equations. Because the scaling of the equations is performed before the time and space discretizations, its benefits are reaped for all time integration schemes. The augmented Lagrangian term is shown to be indispensable if the solution of the linearized system of equations is to be performed without pivoting, a requirement for the efficient solution of the sparse system of linear equations. Finally, a number of numerical examples demonstrate the efficiency of the proposed approach to scaling.

A COMPOSITE TIME INTEGRATION SCHEME FOR DYNAMIC ADHESION AND ITS APPLICATION TO GECKO SPATULA PEELING

2014 ◽  
Vol 11 (05) ◽  
pp. 1350104 ◽  
Author(s):  
SACHIN S. GAUTAM ◽  
ROGER A. SAUER
Keyword(s):  
Nonlinear Response ◽  
Time Integration ◽  
Computational Cost ◽  
Integration Scheme ◽  
Time Step ◽  
Integration Schemes ◽  
Newmark Scheme ◽  
Highly Nonlinear ◽  
Time Integration Schemes ◽  
Time Integration Scheme

Simulation of dynamic adhesive peeling problems at small scales has attracted little attention so far. These problems are characterized by a highly nonlinear response. Accurate and stable time integration schemes are required for simulation of dynamic peeling problems. In the present work, a composite time integration scheme is proposed for the simulation of dynamic adhesive peeling problems. It is shown through numerical examples that the proposed scheme remains stable and also has some gain in accuracy. The performance of the scheme is compared with two collocation-based schemes, i.e., Newmark scheme and Bathe composite scheme. It is shown that the proposed scheme and Bathe composite scheme perform equally. However, the proposed scheme adds very little to the computational cost of Newmark scheme. Through a numerical simulation of the peeling of a gecko spatula from a rigid substrate it is shown that the proposed scheme and the Bathe composite scheme are able to simulate the complete peeling process for given time step whereas the Newmark scheme diverges. It is also shown that the maximum pull-off force is within the range reported in the literature.

Modeling and Detection of Localized Tooth Defects in Geared Systems

1999 ◽  
Author(s):  
M. El Badaoui ◽  
V. Cahouet ◽  
F. Guillet ◽  
J. Daniere ◽  
P. Velex
Keyword(s):  
Numerical Simulations ◽  
Gear Tooth ◽  
Equations Of Motion ◽  
Vibration Monitoring ◽  
Detection Methods ◽  
Implicit Time ◽  
Integration Scheme ◽  
Time Step ◽  
Normal Contact ◽  
Detection Techniques

Abstract The early detection of failures in geared systems is an important industrial problem which has still to be addressed from both an experimental and theoretical viewpoint. The proposed paper combines some extensive numerical simulations of a single stage geared unit with localized tooth faults and the use of several detection techniques whose performances are compared and critically assessed. A model aimed at simulating the contributions of local tooth defects such as spalling to the gear dynamic behavior is set up. The pinion and the gear of a pair are modeled as to two rigid cylinders with all six degrees of freedom connected by a series of springs which represent gear body and gear tooth compliances on the base plane. Classical shaft finite elements including torsional, flexural and axial displacements can be superimposed on the gear element together with some lumped stiffnesses, masses, inertias, ... which account for the load machines, bearings and couplings. Tooth defects are modeled by a distribution of normal deviations over a zone which can be located anywhere on the active tooth flanks. Among the numerous available signal processing techniques used in vibration monitoring, cepstrum analysis is sensitive, reliable and it can be adapted to a complex geared system with several meshes. From an analytical analysis of the equations of motion, two complementary detection techniques based upon the acceleration power cesptrum are proposed. The equations of motion and the contact problem between mating flanks are simultaneously solved by coupling an implicit time-step integration scheme and a unilateral normal contact algorithm. The results of the numerical simulations are used as a data base for the proposed detection techniques. The combined influence of the defect location, depth and extent is analyzed for two examples of spur and helical gears with various profile modifications and the effectiveness of the two complementary detection methods is discussed before some conclusions are drawn.

Sign in / Sign up

Export Citation Format

Share Document