A Newton algorithm is discussed. This algorithm allows for studying nonlinear material behavior, such as plasticity and viscoplasticity. The material behavior is described by so-called internal variables. The model ap...
详细信息
A Newton algorithm is discussed. This algorithm allows for studying nonlinear material behavior, such as plasticity and viscoplasticity. The material behavior is described by so-called internal variables. The model applied in this paper was developed by Lemaitre and Chaboche. The Newton algorithm is implemented as a user-subroutine in the commercial finite-element method package ABAQUS. The results of a numerical analysis of the rotation of an aircraft turbine disk are presented. Cyclic loading conditions are studied, and the results are compared with numerical results in the literature.
We use biorthogonal filter banks to solve hyperbolic PDEs adaptively with a sparse multilevel representation of the solution. The methods described are of finite difference type, and the filter banks are used to give ...
详细信息
We use biorthogonal filter banks to solve hyperbolic PDEs adaptively with a sparse multilevel representation of the solution. The methods described are of finite difference type, and the filter banks are used to give a sparse representation of signals and to transform between grids on different scales. We derive bounds for the error and number of coefficients in the sparse representation. These bounds also apply for filter banks that are not associated with any wavelets. We develop algorithms for fast differentiation and multiplication in detail. The strength of the method is shown in various test problems.
作者:
de Hoop, ATLager, IEDelft Univ Technol
Fac Informat Technol & Syst Lab Electromagnet Res NL-2628 CD Delft Netherlands Delft Univ Technol
Fac Informat Technol & Syst Int Res Ctr Telecommun Transmiss & Radar NL-2628 CD Delft Netherlands
The uniqueness properties of a class of nonlinear electrostatic field problems are investigated. The study was motivated by the development of numerical algorithms to analyse the performance of nonlinear semiconductin...
详细信息
The uniqueness properties of a class of nonlinear electrostatic field problems are investigated. The study was motivated by the development of numerical algorithms to analyse the performance of nonlinear semiconducting electron devices. Here, existence and uniqueness of the solution are prerequisites for the numerical results to have any meaning at all.
An integration scheme is presented for modeling the texture evolution and stress-strain response of elasto-viscoplastic polycrystalline materials. Single crystal kinematics based on a multiplicative decomposition of t...
详细信息
An integration scheme is presented for modeling the texture evolution and stress-strain response of elasto-viscoplastic polycrystalline materials. Single crystal kinematics based on a multiplicative decomposition of the deformation gradient is used to obtain an evolution equation for the crystal elastic deformation gradient. An implicit scheme to integrate this equation is presented which is stable and efficient. The reorientation of the crystal as well as the elastic strain can then be obtained from a polar decomposition of the elastic deformation gradient. numerical studies are presented using material parameters for aluminum (FCC crystals) and zircaloy (HCP crystals) to demonstrate the general nature of the model. Predictions of the model are also compared with those obtained using a rigid-viscoplastic polycrystal model which neglects the elastic response. Retaining the elastic response makes the model useful for large deformation analyses where both anisotropy due to texture as well as elastic effects such as springback and residual stresses are important. (C) 1999 Elsevier Science Ltd, All rights reserved.
An assessment of the accuracy of shock-capturing schemes is made for two-dimensional steady flow around a cylindrical projectile. Both a linear fourth-order method and a nonlinear third-order method are used in this s...
详细信息
An assessment of the accuracy of shock-capturing schemes is made for two-dimensional steady flow around a cylindrical projectile. Both a linear fourth-order method and a nonlinear third-order method are used in this study. This study shows, contrary to conventional wisdom, that captured two-dimensional shocks are asymptotically first order, regardless of the design accuracy of the numerical method. The practical implications of this finding are discussed in the context of the efficacy of high-order numerical methods for discontinuous flows.
The stability of post-buckled states for simply-supported flat elastic plates under compression is investigated for a range of in-plane boundary conditions. The von Karman plate equations are reduced to a series of OD...
详细信息
The stability of post-buckled states for simply-supported flat elastic plates under compression is investigated for a range of in-plane boundary conditions. The von Karman plate equations are reduced to a series of ODEs which are solved numerically under parametric variation of both load and length. Results are checked against full numerical solutions of the PDEs, and comparison with a modal analysis highlights the dominant passive contaminations. The nondimensional amplitude at secondary bifurcation, for any combination of modes and all plate lengths, is presented in a concise form using the parameter space of Arnold tongues. This demonstrates that compound bifurcation represents a worst case for post-buckling reserve, and that long plates have inherently more such reserve than short plates. It is also shown that stiffening the boundaries against inplane movement is destabilizing, in that it induces mode jumping at secondary bifurcation to occur at an earlier stage in the post-buckling regime. (C) 1999 Elsevier Science Ltd. All rights reserved.
Presented herein is an adaptive-mesh computational method for the efficient solution of the continuum equations of compressible flow for high-velocity impact dynamics. The integral forms of the governing equations are...
详细信息
Presented herein is an adaptive-mesh computational method for the efficient solution of the continuum equations of compressible flow for high-velocity impact dynamics. The integral forms of the governing equations are used to derive a stable form of energy equation, using internal rather than total energy, after which the corresponding differential forms are solved approximately in two dimensions via a three-stage-pressure, stress and advection-finite-difference scheme. The finite-difference equations are applied on a fully-unstructured adaptive mesh which. as time proceeds, both coarsens or refines locally and automatically in response to a prescribed adaption criterion. The truncation errors of the scheme are studied, and the scheme is first verified on a test problem of a collapsing hollow spherical shell, for which an analytical solution is known, before being applied to more general configurations motivated by the study of penetration mechanics. Results indicate that, for both CPU and memory requirements, the adaptive scheme is considerably cheaper than the corresponding comparable-resolution regular-mesh scheme, both schemes calculating results to a similar degree of accuracy. (C) 1999 Elsevier Science Ltd. All rights reserved.
A new algorithm is introduced to integrate the equations of rotational motion. The algorithm is derived within a leapfrog framework and the quantities involved into the integration are mid-step angular momenta and on-...
详细信息
A new algorithm is introduced to integrate the equations of rotational motion. The algorithm is derived within a leapfrog framework and the quantities involved into the integration are mid-step angular momenta and on-step orientational positions. Contrary to the standard implicit method by Fincham [Mol. Sim., 8, 165 (1992)], the revised angular momentum approach presented corresponds completely to the leapfrog idea on interpolation of dynamical variables without using any extrapolations. The proposed scheme intrinsically preserves rigid molecular structures and considerably improves stability properties and energy conservation. As is demonstrated on the basis of simulations for water, it reproduces correct results even with extra large step sizes of order 5 fs and 10 fs in the cases of energy- and temperature-conserving dynamics, respectively. We show also that iterative solutions can be avoided within our implicit scheme shifting from quaternions to the entire rotation-matrix representation.
An unsteady viscous how solver based on the Runge-Kutta scheme has been developed. Three low-Reynolds-number k-epsilon turbulence models have been assessed for their ability to predict unsteady transitional how The k-...
详细信息
An unsteady viscous how solver based on the Runge-Kutta scheme has been developed. Three low-Reynolds-number k-epsilon turbulence models have been assessed for their ability to predict unsteady transitional how The k-epsilon model has been modified to incorporate leading-edge effects. A pseudo time approach has been used to accelerate the numerical algorithm for unsteady flows. A sensitivity study has been carried out, including the effects of the artificial dissipation, the grid density, and the physical and pseudo time steps. The numerical simulation of the unsteady transitional flow in a compressor cascade has been carried out successfully. Good agreement with the measured data as well as the previous Euler/boundary-layer predictions has been achieved. The numerical solver has been able to predict major features, associated with the wake-induced transition on a compressor blade (wake-induced transitional strip, wake-induced turbulent strip, etc.). An analysis and an interpretation of the results from the unsteady flow simulation have been carried out to understand additional how physics associated with the transitional flow.
Conventionally, the use of continuous distributions of dislocations to model plasticity has been confined to the analysis of crack tip plasticity using linear arrays of dislocations, within the framework of plane anal...
详细信息
Conventionally, the use of continuous distributions of dislocations to model plasticity has been confined to the analysis of crack tip plasticity using linear arrays of dislocations, within the framework of plane analysis. By expanding this technique into a distribution of dislocation over an area, a method is developed to model the plasticity at stress raising feat;res such as notches or holes under plane strain conditions. The method explicitly takes account of the boundary conditions by using a dislocation solution which accounts for the presence of the stress-raise itself. Other free boundaries may be modelled more approximately using boundary elements which also correctly include the presence of the stress raiser. The dislocations are distributed over finite sized cells, and the solutions found for the strain fields compare favourably with both finite element and bounding Neuber and Glinka results. (C) 1999 Elsevier Science Ltd. All rights reserved.
暂无评论