A constitutive model with three plastic constants: The description of anisotropic workhardening

Constitutive relations are suggested with three plasticity constants for numerical modeling of anisotropic workhardening characterized by translation, reshaping and turning of the subsequent yield surface. The yield surface is assumed in the shape of an ellipsoid of revolution in a vector space of stress deviators, its axes of revolution coinciding with the plastic deformation vector. The yield condition is expressed through the parameters of dimensions, location, and orientation of the hyperellipsoid, and is a natural development of the equation for a spherical yield surface considered by the conventional theory of combined isotropic-kinematic hardening. Methodology is also suggested for determining plastic constants though proportional loading experiment. Properties of the constitutive model are analyzed for proportional loading, and a table is constructed illustrating qualitatively different types of evolution of the subsequent yield surface with respective numerical values of plastic constants determining them. The effectiveness of the constitutive model suggested is confirmed by the examples of numerical modeling with experimental results available in the literature.     

Analysis and numerical simulation of magnetic forces between rigid polygonal bodies. Part II: Numerical simulation

The analysis of magnetoelastic phenomena is a field of active research. Formulae for the magnetic force in macroscopic systems have been under discussion for some time. In Popović et al. (Continum. Mech. Thermodyn. 2007), we rigorously justify several of the available formulae in the context of rigid bodies in two and three space dimensions. In the present, second part of our study, we investigate these formulae in a series of numerical experiments in which the magnetic force is computed in dependence on the geometries of the bodies as well as on the distance between them. In case the two bodies are in contact, i.e., in the limit as their distance tends to zero, we focus especially on a formula obtained in a discrete-to-continuum approximation. The aim of our study is to help clarify the question which force formula is the correct one in the sense that it describes nature most accurately and to suggest adequate real-life experiments for a comparison with the provided numerical data.      

A theory of large-strain isotropic thermoplasticity based on metric transformation tensors

Summary A formulation of isotropic thermoplasticity for arbitrary large elastic and plastic strains is presented. The underlying concept is the introduction of a metric transformation tensor which maps a locally defined six-dimensional plastic metric onto the metric of the current configuration. This mixed-variant tensor field provides a basis for the definition of a local isotropic hyperelastic stress response in the thermoplastic solid. Following this fundamental assumption, we derive a consistent internal variable formulation of thermoplasticity in a Lagrangian as well as a Eulerian geometric setting. On the numerical side, we discuss in detail an objective integration algorithm for the mixed-variant plastic flow rule. The special feature here is a new representation of the stress return and the algorithmic elastoplastic moduli in the eigenvalue space of the Eulerian plastic metric for plane problems. Furthermore, an algorithm for the solution of the coupled problem is formulated based on an operator split of the global field equations of thermoplasticity. The paper concludes with two representative numerical simulations of thermoplastic deformation processes.     

Riemann solvers with evolved initial conditions

The scope of this paper is three fold. We first formulate upwind and symmetric schemes for hyperbolic equations with non‐conservative terms. Then we propose upwind numerical schemes for conservative and non‐conservative systems, based on a Riemann solver, the initial conditions of which are evolved non‐linearly in time, prior to a simple linearization that leads to closed‐form solutions. The Riemann solver is easily applied to complicated hyperbolic systems. Finally, as an example, we formulate conservative schemes for the three‐dimensional Euler equations for general compressible materials and give numerical results for a variety of test problems for ideal gases in one and two space dimensions. Copyright © 2006 John Wiley & Sons, Ltd.     

A volume-fraction based algorithm for hybrid barotropic and non-barotropic two-fluid flow problems

The aim of this paper is to describe a simple Eulerian interface-capturing approach for the efficient numerical resolution of a hybrid barotropic and non-barotropic two-fluid flow problem in more than one space dimension. We use the compressible Euler equations as a model system with the thermodynamic property of each of the barotropic and non-barotropic fluid components characterized by the Tait and Noble–Abel equations of state, respectively. The algorithm is based on a volume fraction formulation of the equations together with an extended equation of state that is devised to give an approximate treatment for the mixture of more than one fluid component within a grid cell. A standard high-resolution wave propagation method is employed to solve the proposed two-fluid model with the dimensional-splitting technique incorporated in the method for multidimensional problems. Several numerical results are presented in one and two space dimensions that show the feasibility of the algorithm as applied to a reasonable class of practical problems without the occurrence of any spurious oscillation in the pressure near the smeared material interfaces. This includes, in particular, solutions for a study on the variation of the jet velocity with the incident shock pressure arising from the collapse of an air cavity in water under a shock wave.     

Approximation of generalized Stokes problems using dual‐mixed finite elements without enrichment

In this work a finite element method for a dual‐mixed approximation of generalized Stokes problems in two or three space dimensions is studied. A variational formulation of the generalized Stokes problems is accomplished through the introduction of the pseudostress and the trace‐free velocity gradient as unknowns, yielding a twofold saddle point problem. The method avoids the explicit computation of the pressure, which can be recovered through a simple post‐processing technique. Compared with an existing approach for the same problem, the method presented here reduces the global number of degrees of freedom by up to one‐third in two space dimensions. The method presented here also represents a connection between existing dual‐mixed and pseudostress methods for Stokes problems. Existence, uniqueness, and error results for the generalized Stokes problems are given, and numerical experiments that illustrate the theoretical results are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical implementation of non-local polycrystal plasticity using fast Fourier transforms

We present the numerical implementation of a non-local polycrystal plasticity theory using the FFT-based formulation of Suquet and co-workers. Gurtin (2002) non-local formulation, with geometry changes neglected, has been incorporated in the EVP-FFT algorithm of Lebensohn et al. (2012). Numerical procedures for the accurate estimation of higher order derivatives of micromechanical fields, required for feedback into single crystal constitutive relations, are identified and applied. A simple case of a periodic laminate made of two fcc crystals with different plastic properties is first used to assess the soundness and numerical stability of the proposed algorithm and to study the influence of different model parameters on the predictions of the non-local model. Different behaviors at grain boundaries are explored, and the one consistent with the micro-clamped condition gives the most pronounced size effect. The formulation is applied next to 3-D fcc polycrystals, illustrating the possibilities offered by the proposed numerical scheme to analyze the mechanical response of polycrystalline aggregates in three dimensions accounting for size dependence arising from plastic strain gradients with reasonable computing times.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

An efficient high order direct ALE ADER finite volume scheme with a posteriori limiting for hydrodynamics and magnetohydrodynamics

In this paper, we present a new family of direct arbitrary–Lagrangian–Eulerian (ALE) finite volume schemes for the solution of hyperbolic balance laws on unstructured meshes in multiple space dimensions. The scheme is designed to be high‐order accurate both in space and time, and the mesh motion, which provides the new mesh configuration at the next time step, is taken into account in the final finite volume scheme that is based directly on a space‐time conservation formulation of the governing PDE system. To improve the computational efficiency of the algorithm, high order of accuracy in space is achieved using the a posteriori MOOD limiting strategy that allows the reconstruction procedure to be carried out with only one reconstruction stencil for any order of accuracy. We rely on an element‐local space‐time Galerkin finite element predictor on moving curved meshes to obtain a high‐order accurate one‐step time discretization, while the mesh velocity is computed by means of a suitable nodal solver algorithm that might also be supplemented with a local rezoning procedure to improve the mesh quality. Next, the old mesh configuration at time level tⁿ is connected to the new one at t^n + 1 by straight edges, hence providing unstructured space‐time control volumes, on the boundary of which the numerical flux has to be integrated. Here, we adopt a quadrature‐free integration, in which the space‐time boundaries of the control volumes are split into simplex sub‐elements that yield constant space‐time normal vectors and Jacobian matrices. In this way, the integrals over the simplex sub‐elements can be evaluated once and for all analytically during a preprocessing step. We apply the new high‐order direct ALE algorithm to the Euler equations of compressible gas dynamics (also referred to as hydrodynamics equations) as well as to the magnetohydrodynamics equations and we solve a set of classical test problems in two and three space dimensions. Numerical convergence rates are provided up to fifth order of accuracy in 2D and 3D for both hyperbolic systems considered in this paper. Finally, the efficiency of the new method is measured and carefully compared against the original formulation of the algorithm that makes use of a WENO reconstruction technique and Gaussian quadrature formulae for the flux integration: depending on the test problem, the new class of very efficient direct ALE schemes proposed in this paper can run up to ≈12 times faster in the 3D case. Copyright © 2016 John Wiley & Sons, Ltd.     

Acoustic scattering and radiation problems, and the null-field method

The best known methods for solving the scattering and radiation problems of acoustics are integral-equation methods. However, it is also known that the simplest of these methods yield equations which are not uniquely solvable at certain discrete sets of frequencies (the irregular frequencies). In this paper, we shall analyse an alternative method (the null-field method, or T-matrix method). We prove that the infinite system of null-field equations always has precisely one solution, i.e. the unphysical irregular frequencies do not occur with this method. Moreover, we also prove that the solution of the original boundary-value problem can always be determined (at any point exterior to the scatterer) from the solution of the null-field equations. We prove these results in two dimensions, for two radiation problems (the exterior Neumann problem and the exterior Dirichlet problem) and two scattering problems (scattering by a sound-hard body and scattering by a sound-soft body); similar results can be proved in three dimensions. We also prove some subsidiary results, concerning the solvability of certain boundary integral equations and the completeness of certain sets of radiating wave-functions, and give a discussion of related numerical techniques.     

Analysis and numerical simulation of magnetic forces between rigid polygonal bodies. Part I: Analysis

The mathematical and physical analysis of magnetoelastic phenomena is a topic of ongoing research. Different formulae have been proposed to describe the magnetic forces in macroscopic systems. We discuss several of these formulae in the context of rigid magnetized bodies. In case the bodies are in contact, we consider formulae both in the framework of macroscopic electrodynamics and via a multiscale approach, i.e., in a discrete setting of magnetic dipole moments. We give mathematically rigorous proofs for domains of polygonal shape (as well as for more general geometries) in two and three space dimensions. In an accompanying second article, we investigate the formulae in a number of numerical experiments, where we focus on the dependence of the magnetic force on the distance between the bodies and on the case when the two bodies are in contact. The aim of the analysis as well as of the numerical simulation is to contribute to the ongoing debate about which formula describes the magnetic force between macroscopic bodies best and to stimulate corresponding real-life experiments.      

Numerical study of dynamic phase transitions in shock tube

Shock tube problem of a van der Waals fluid with a relaxation model was investigated.In the limit of relaxation parameter tending towards zero,this model yields a specific Riemann solver.Relaxing and relaxed schemes were derived.For an inci- dent shock in a fixed tube,numerical simulations show convergence toward the Riemann solution in one space dimension.Impact of parameters was studied theoretically and nu- merically.For certain initial shock profiles,nonclassical reflecting wave was observed.In two space dimensions,the effect of curved wave fronts was studied,and some interesting wave patterns were exposed.     

A Newton iterative mixed finite element method for stationary conduction–convection problems

In this article, a Newton iterative mixed finite element method is presented for solving the stationary conduction–convection problems in two dimensions. The stability and the errors generated by both partitioning the space and solving nonlinear equations are analysed, which show that our method is stable and has good precision. Finally, some numerical experiments are given to confirm its effect.     

Mixed analytical–numerical relaxation in finite single-slip crystal plasticity

The modeling of the finite elastoplastic behaviour of single crystals with one active slip system leads to a nonconvex variational problem, whose minimization produces fine structures. The computation of the quasiconvex envelope of the energy density involves the solution of a nonconvex optimization problem and faces severe numerical difficulties from the presence of many local minima. In this paper, we consider a standard model problem in two dimensions and, by exploiting analytical relaxation results for limiting cases and the special structure of the problem at hand, we obtain a fast and efficient numerical relaxation algorithm. The effectiveness of our algorithm is demonstrated with numerical examples. The precision of the results is assessed by lower bounds from polyconvexity.      

FULLY DISCRETE HIGH-RESOLUTION SCHEMES FOR HYPERBOLIC CONSERVATION LAWS

The present paper is a sequel to two previous papers in which rigorous, up to fourth-order, fully discrete (FD) upwind TVD schemes have been presented. In this paper we discuss in detail the extension of these schemes to solutions of non-linear hyperbolic systems. The performance of the schemes is assessed by solving test problems for the time-dependent Euler equations of gas dynamics in one and two space dimensions. We use exact solutions and experimental data to validate the results.     

Biaxial Vibrations of an Elasto-Plastic Beam with a Prescribed Rigid-Body Rotation Including the Effect of Stiffening

In the present work, we study biaxial vibrations of elasto-plastic beamswith a prescribed rigid-body motion. Exemplarily, we treat the case of ahinged-hinged beam with a rotation about a fixed hinge. Such a problemis frequently to be encountered in machine dynamics. The deformations ofthe beam are assumed to remain small. In order to describe the beamvibrations, we utilize the Bernoulli–Euler beam theory. Axial vibrationsof the beam are neglected. We study cases, in which the flexuralstiffness of the beam is considerably lowered due to catastrophicenvironmental influences, such that the deformations relative to therigid-body motion, albeit small, reach the plastic regime. In thepresent paper, special emphasis is laid upon including the effect ofgeometric stiffening due to a comparatively fast rigid-body rotation.The equations of motion are derived by Hamilton's principle. The biaxialdeflections are discretized in space by means of Legendre polynomials.The plastic strains are discretized over length, height and width of thebeam by small plastic cells. Galerkin's procedure, together with aproper implicit midpoint rule is used for integration of the equationsof motion. The plastic strains are computed in every time-step by asuitable iterative procedure. Linear elastic/perfectly plastic behavioris exemplarily treated in a numerical study.     

两圆柱体结合面的接触热导分形模型研究

基于三维分形理论,在考虑微凸体的弹性变形、弹塑性变形和塑性变形的基础上,建立了两圆柱体结合面接触热导分形模型。通过数值仿真,分析了分形维数,分形尺度参数、圆柱体曲率半径和接触类型对接触热导的影响。研究结果表明：接触热导随着分形维数的增大而增大,随着分形尺度参数的增大而减小;相同参数下,内接触比外接触的接触热导要大;此外,当固定其中一个圆柱体的曲率半径时,随着另一个圆柱体曲率半径的增大,接触热导增大。该模型为开展齿轮等曲面接触热导的研究提供了理论基础。     

Decay Rates for the Three‐Dimensional Linear System of Thermoelasticity

. We consider the two and three‐dimensional system of linear thermoelasticity in a bounded smooth domain with Dirichlet boundary conditions. We analyze whether the energy of solutions decays exponentially uniformly to zero as . First of all, by a decoupling method, we reduce the problem to an observability inequality for the Lamé system in linear elasticity and more precisely to whether the total energy of the solutions can be estimated in terms of the energy concentrated on its longitudinal component. We show that when the domain is convex, the decay rate is never uniform. In fact, the lack of uniform decay holds in a more general class of domains in which there exist rays of geometric optics of arbitrarily large length that are always reflected perpendicularly or almost tangentially on the boundary. We also show that, in three space dimensions, the lack of uniform decay may also be due to a critical polarization of the energy on the transversal component of the displacement. In two space dimensions we prove a sufficient (and almost necessary) condition for the uniform decay to hold in terms of the propagation of the transversal characteristic rays, under the further assumption that the boundary of the domain does not have contacts of infinite order with its tangents. We also give an example, due to D. Hulin, in which these geometric properties hold. In three space dimensions we indicate (without proof) how a careful analysis of the polarization of singularities may lead to sharp sufficient conditions for the uniform decay to hold. In two space dimensions we prove that smooth solutions decay polynomially in the energy space to a finite‐dimensional subspace of solutions except when the domain is a ball or an annulus. Finally we discuss some closely related controllability and spectral issues. (Accepted May 14, 1998)     

Clifford algebra,theory of its function and their applicaton to mechanics

As is Wellknown in both elastic mechanics andfluid mechanics, the plane problems are more convenient than space problems.One of the causes is that there has been a complete theory about the complex Junction and the analytic junction, hut in space problems, the case is quite different.We have no effective method to deal with these problems.In this paper, we first introduces general theories of Clifford algebra.Then we emphatically explain Clifford algebra in three dimensions and establish theories of regular Junction in three dimensions analogically to analytic function in plane.Thus we extend some results of plane problem-la three dimensions or high dimensions.Obviously, it is very important for elastic and fluid mechanics.But because Clifford algebra is not a commutative algebra, we can’t simply extend the results of two dimensions to high dimensions.The left problems are yet to be found out.     

拉－扭复合加载下不锈钢的弹塑性本构关系——Ⅱ．理论

提出应力是塑性应变空间内蕴几何学参数的泛函．一般情况下，塑性应变空间是非欧几何空间，而其度量张量是塑性应变和其历史的函数，但在初始各向同性和不可压的情况下可取成欧氏空间．本文在Ｉｌｙｕｓｈｉｎ理论，和Ｖａｌａｎｉｓ理论的基础上，提出在拉－扭复合加载下的εｐ１－εｐ３空间中新的积分型弹塑性本构关系，所建理论预测的结果和实验［１］相当一致，表明理论是合理的