Numerical Integration of Some Stiff Constitutive Models of Inelastic Deformation

1980 ◽  
Vol 102 (1) ◽  
pp. 92-96 ◽  
Author(s):  
Virendra Kumar ◽  
Mahesh Morjaria ◽  
Subrata Mukherjee
Keyword(s):  
Numerical Integration ◽  
Inelastic Deformation ◽  
Time Integration ◽  
Constitutive Models ◽  
Integration Scheme ◽  
Uniaxial Deformation ◽  
Time Step ◽  
One Step ◽  
Numerical Time Integration ◽  
Step Control

Several strategies for numerical time-integration of some stiff constitutive models of inelastic deformation are presented in this paper. Numerical results and comparisons are presented for the integration of one such model for the case of uniaxial deformation under various prescribed histories of stress or strain. A simple one step Euler type integration scheme with automatic time-step control, which can be easily adapted to the solution of multiaxial boundary value problems, appears promising.

On an Improved Time Integration Scheme for Stiff Constitutive Models of Inelastic Deformation

1985 ◽  
Vol 107 (4) ◽  
pp. 282-285 ◽  
Author(s):  
Vinod Banthia ◽  
Subrata Mukherjee
Keyword(s):  
Strain Rate ◽  
Inelastic Deformation ◽  
Time Integration ◽  
Constitutive Models ◽  
Size Control ◽  
Integration Scheme ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Step Size Control

For the time-integration of stiff constitutive models of inelastic deformation, the explicit (one step Euler) integration scheme can be used provided the time step size is closely monitored and controlled. The time step size control scheme based on prescribed error bounds is of limited use because it requires an a priori estimate of the maximum nonelastic strain rate for the selection of a proper error bound. In this paper, a new scheme for time-step size control is presented. This scheme automatically scales the time-step size by the maximum nonelastic strain rate. That the new scheme is superior to the old one is evident from the results of the various problems presented here.

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 EULER AND EXPONENTIAL ALGORITHM IN VISCOELASTIC ANALYSES OF LAMINATED GLASS

2020 ◽  
Vol 26 ◽  
pp. 86-93
Author(s):  
Jaroslav Schmidt ◽  
Alena Zemanová
Keyword(s):  
Time Integration ◽  
Optimal Time ◽  
Integration Scheme ◽  
Laminated Glass ◽  
Aesthetic Value ◽  
Rate Dependent ◽  
Time Integration Scheme ◽  
One Step ◽  
Sectional Plane ◽  
The Aesthetic

Laminated glass combines two remarkable materials: glass and a polymer ply. While glass is stiff and brittle, the polymer ply is a rate-dependent compliant material. Together, they form a material which keeps the aesthetic value of glass, and due to the polymer, no fragile collapse appears. The polymer ply exhibits time- and temperature-dependency, whereas glass suffers from brittle fracture, which makes the analysis difficult. In this article, a 2D sectional plane-stress model for the viscoelastic analysis of laminated glass is presented. This study presents one step in the development of a phase-field-based damage solver for laminated glass to select the optimal time-integration scheme for a quasistatic-analysis and later also for dynamics. The validation against experimental data is provided, and the model reduction is also discussed.

A New Unconditionally Stable Time Integration Method for Analysis of Nonlinear Structural Dynamics

2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Ali Akbar Gholampour ◽  
Mehdi Ghassemieh ◽  
Mahdi Karimi-Rad
Keyword(s):  
Time Integration ◽  
Second Order ◽  
Numerical Dispersion ◽  
Computational Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Similar Time ◽  
Unconditionally Stable ◽  
Average Acceleration

A new time integration scheme is presented for solving the differential equation of motion with nonlinear stiffness. In this new implicit method, it is assumed that the acceleration varies quadratically within each time step. By increasing the order of acceleration, more terms of the Taylor series are used, which are expected to have responses with better accuracy than the classical methods. By considering this assumption and employing two parameters δ and α, a new family of unconditionally stable schemes is obtained. The order of accuracy, numerical dissipation, and numerical dispersion are used to measure the accuracy of the proposed method. Second order accuracy is achieved for all values of δ and α. The proposed method presents less dissipation at the lower modes in comparison with Newmark's average acceleration, Wilson-θ, and generalized-α methods. Moreover, this second order accurate method can control numerical damping in the higher modes. The numerical dispersion of the proposed method is compared with three unconditionally stable methods, namely, Newmark's average acceleration, Wilson-θ, and generalized-α methods. Furthermore, the overshooting effect of the proposed method is compared with these methods. By evaluating the computational time for analysis with similar time step duration, the proposed method is shown to be faster in comparison with the other methods.

Parametric studies on dynamic stiffness of ball bearings supporting a flexible rotor

2019 ◽  
Vol 25 (15) ◽  
pp. 2175-2188 ◽  
Author(s):  
Tikam Chand Gupta
Keyword(s):  
Numerical Integration ◽  
Dynamic Stiffness ◽  
Mathematical Formulation ◽  
Rolling Element Bearing ◽  
Radial Clearance ◽  
Integration Scheme ◽  
Time Step ◽  
Average Value ◽  
Bearing Stiffness ◽  
Time Invariant

Most of the researchers in the field of dynamics of the rolling element bearing have considered bearing stiffness as time invariant and/or not related to dynamics of the bearing. In the present paper, the bearing stiffness has been taken as function of dynamic response at every time step of numerical simulation and a detailed parametric study is performed to investigate the effect of flexibility of the rotor shaft, rotational speed, and internal radial clearance on the instantaneous and average value of dynamic stiffness of the ball bearing. The mathematical formulation is based on the Timoshenko beam finite element theory. Gravity and bearing forces are considered as external forces acting on a free-free flexible shaft. A stable Newmark- β numerical integration scheme coupled with Newton–Raphson method is used for numerical integration and for convergence to an accurate value of bearing stiffness. The results showing the variation of different components of bearing stiffness as a function of time-invariant parameters has improved the understanding of the dynamic behavior of the bearing during motion. The variation pattern of bearing stiffness coefficients is observed to be sensitive to direction of rotation. The amplitude of periodic change of these coefficients increases with the increase of the stiffness ratio of shaft and the decrease of radial clearance.

Newmark-Beta Method in Discrete Elastic Rods Algorithm to Avoid Energy Dissipation

2019 ◽  
Vol 86 (8) ◽  
Author(s):  
Weicheng Huang ◽  
Mohammad Khalid Jawed
Keyword(s):  
Energy Dissipation ◽  
Time Integration ◽  
Integration Scheme ◽  
Integration Step ◽  
Elastic Rods ◽  
Computationally Efficient ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Time Integration Scheme

Discrete elastic rods (DER) algorithm presents a computationally efficient means of simulating the geometrically nonlinear dynamics of elastic rods. However, it can suffer from artificial energy loss during the time integration step. Our approach extends the existing DER technique by using a different time integration scheme—we consider a second-order, implicit Newmark-beta method to avoid energy dissipation. This treatment shows better convergence with time step size, specially when the damping forces are negligible and the structure undergoes vibratory motion. Two demonstrations—a cantilever beam and a helical rod hanging under gravity—are used to show the effectiveness of the modified discrete elastic rods simulator.

Separate Time Integration Based on the Hilber, Hughes, Taylor Scheme for Flexible Bodies With a Large Number of Modes

2020 ◽  
Vol 15 (6) ◽  
Author(s):  
Wolfgang Witteveen ◽  
Florian Pichler
Keyword(s):  
Time Integration ◽  
Flexible Multibody Dynamics ◽  
Integration Scheme ◽  
Step Size ◽  
Body Contact ◽  
Time Integration Scheme ◽  
Core Idea ◽  
Numerical Time Integration ◽  
Taylor Scheme ◽  
Local Deformations

Abstract In the current development of flexible multibody dynamics, the efficient and accurate consideration of distributed and nonlinear forces is an active area of research. Examples are, forces due to body-body contact or due to elastohydrodynamics (EHD). This leads to many additional modes for representing the local deformations in the areas on which those forces act. Recent publications show that these can be several hundred to several thousand additional modes. A conventional, monolithic numerical time integration scheme would lead to unacceptable computing times. This paper presents a method for an efficient time integration of such systems. The core idea is to treat the equations associated with modes representing local deformations separately. Using the Newmark formulas, a fixed point iteration is proposed for these separated equations, which can always be stabilized with decreasing step size. The concluding examples underline this property, as well as the fact that the proposed method massively outperforms the conventional, monolithic time integration with increasing number of modes.

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.

Sign in / Sign up

Export Citation Format

Share Document