Asymptotic integration algorithms for nonhomogeneous, nonlinear, first order, ordinary differential equations

New methods for integrating systems of stiff, nonlinear, first order, ordinary differential equations are developed by casting the differential equations into integral form. Nonlinear recursive relations are obtained that allow the solution to a system of equations at time t plus delta t to be obtained in terms of the solution at time t in explicit and implicit forms. Examples of accuracy obtained with the new technique are given by considering systems of nonlinear, first order equations which arise in the study of unified models of viscoplastic behaviors, the spread of the AIDS virus, and predator-prey populations. In general, the new implicit algorithm is unconditionally stable, and has a Jacobian of smaller dimension than that which is acquired by current implicit methods, such as the Euler backward difference algorithm; yet, it gives superior accuracy. The asymptotic explicit and implicit algorithms are suitable for solutions that are of the growing and decaying exponential kinds, respectively, whilst the implicit Euler-Maclaurin algorithm is superior when the solution oscillates, i.e., when there are regions in which both growing and decaying exponential solutions exist

From differential to difference equations for first order ODEs

When constructing an algorithm for the numerical integration of a differential equation, one should first convert the known ordinary differential equation (ODE) into an ordinary difference equation. Given this difference equation, one can develop an appropriate numerical algorithm. This technical note describes the derivation of two such ordinary difference equations applicable to a first order ODE. The implicit ordinary difference equation has the same asymptotic expansion as the ODE itself, whereas the explicit ordinary difference equation has an asymptotic that is similar in structure but different in value when compared with that of the ODE

A viscoplastic model with application to LiF-22 percent CaF2 hypereutectic salt

A viscoplastic model for class M (metal-like behavior) materials is presented. One novel feature is its use of internal variables to change the stress exponent of creep (where n is approximately = 5) to that of natural creep (where n = 3), in accordance with experimental observations. Another feature is the introduction of a coupling in the evolution equations of the kinematic and isotropic internal variables, making thermal recovery of the kinematic variable implicit. These features enable the viscoplastic model to reduce to that of steady-state creep in closed form. In addition, the hardening parameters associated with the two internal state variables (one scalar-valued, the other tensor-valued) are considered to be functions of state, instead of being taken as constant-valued. This feature enables each internal variable to represent a much wider spectrum of internal states for the material. The model is applied to a LiF-22 percent CaF2 hypereutectic salt, which is being considered as a thermal energy storage material for space-based solar dynamic power systems

Exponential integration algorithms applied to viscoplasticity

Four, linear, exponential, integration algorithms (two implicit, one explicit, and one predictor/corrector) are applied to a viscoplastic model to assess their capabilities. Viscoplasticity comprises a system of coupled, nonlinear, stiff, first order, ordinary differential equations which are a challenge to integrate by any means. Two of the algorithms (the predictor/corrector and one of the implicits) give outstanding results, even for very large time steps

Steady-state and transient Zener parameters in viscoplasticity: Drag strength versus yield strength

A hypothesis is put forth which enables the viscoplastician to formulate a theory of viscoplasticity that reduces, in closed form, to the classical theory of creep. This hypothesis is applied to a variety of drag and yield strength models. Because of two theoretical restrictions that are a consequence of this hypothesis, three different yield strength models and one drag strength model are shown to be theoretically admissible. One of these yield strength models is selected as being the most appropriate representation for isotropic hardening

Asymptotic integration algorithms for first-order ODEs with application to viscoplasticity

When constructing an algorithm for the numerical integration of a differential equation, one must first convert the known ordinary differential equation (ODE), which is defined at a point, into an ordinary difference equation (O(delta)E), which is defined over an interval. Asymptotic, generalized, midpoint, and trapezoidal, O(delta)E algorithms are derived for a nonlinear first order ODE written in the form of a linear ODE. The asymptotic forward (typically underdamped) and backward (typically overdamped) integrators bound these midpoint and trapezoidal integrators, which tend to cancel out unwanted numerical damping by averaging, in some sense, the forward and backward integrations. Viscoplasticity presents itself as a system of nonlinear, coupled first-ordered ODE's that are mathematically stiff, and therefore, difficult to numerically integrate. They are an excellent application for the asymptotic integrators. Considering a general viscoplastic structure, it is demonstrated that one can either integrate the viscoplastic stresses or their associated eigenstrains

Thermoviscoplastic analysis of fibrous periodic composites using triangular subvolumes

The nonlinear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history, the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The nonlinearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions

Stress versus temperature dependent activation energies in creep

The activation energy for creep at low stresses and elevated temperatures is lattice diffusion, where the rate controlling mechanism for deformation is dislocation climb. At higher stresses and intermediate temperatures, the rate controlling mechanism changes from that of dislocation climb to one of obstacle-controlled dislocation glide. Along with this change, there occurs a change in the activation energy. It is shown that a temperature-dependent Gibbs free energy does a good job of correlating steady-state creep data, while a stress-dependent Gibbs free energy does a less desirable job of correlating the same data. Applications are made to copper and a LiF-22 mol. percent CaF2 hypereutectic salt

On the thermodynamics of stress rate in the evolution of back stress in viscoplasticity

A thermodynamic foundation using the concept of internal state variables is presented for the kinematic description of a viscoplastic material. Three different evolution equations for the back stress are considered. The first is that of classical, nonlinear, kinematic hardening. The other two include a contribution that is linear in stress rate. Choosing an appropriate change in variables can remove this stress rate dependence. As a result, one of these two models is shown to be equivalent to the classical, kinematic hardening model; while the other is a new model, one which seems to have favorable characteristics for representing ratchetting behavior. All three models are thermodynamically admissible

The isothermal fatigue behavior of a unidirectional SiC/Ti composite and the Ti alloy matrix

The high temperature fatigue behavior of a metal matrix composite (MMC) consisting of Ti-15V-3Cr-3Al-3Sn (Ti-15-3) matrix reinforced by 33 vol percent of continuous unidirectional SiC fibers was experimentally and analytically evaluated. Isothermal MMC fatigue tests with constant amplitude loading parallel to the fiber direction were performed at 300 and 550 C. Comparative fatigue tests of the Ti-15-3 matrix alloy were also conducted. Composite fatigue behavior and the in-situ stress state of the fiber and matrix were analyzed with a micromechanical model, the Concentric Cylinder Model (CCM). The cyclic stress-strain response of the composite was stable at 300 C. However, an increase in cyclic mean strain foreshortened MMC fatigue life at high strain ranges at 550 C. Fatigue tests of the matrix alloy and CCM analyses indicated this response was associated with stress relaxation of the matrix in the composite