Linear Isotropic Relations in Finite Hyperelasticity: Some General Results

The form of the classical stress–strain relations of linear elasticity are considered here within the context of nonlinear elasticity. For both Cauchy and symmetric Piola-Kirchhoff stresses, conditions are obtained for the associated strain fields so that they are independent of the material constants and compatible with existence of a strain–energy function. These conditions can be integrated in both cases to obtain the most general strain field that satisfies these conditions and the corresponding strain–energy function is obtained. In both cases, a natural choice of form of solution is suggested by the special case of the compatibility conditions being satisfied identically. It will be shown that some strain–energy functions previously introduced in the literature are special cases of the results obtained here. Some recent linear stress–strain relations, proposed in the context of Cauchy elasticity, are examined to see if they are compatible with hyperelasticity.      

Homogenization and Global Responses of Inhomogeneous Spherical Nonlinear Elastic Shells

Homogenization of radially inhomogeneous spherical nonlinear elastic shells subject to internal pressure is studied. The equivalent homogeneous material is defined in such a way that it gives rise to exactly the same global response to the pressure load as that of the inhomogeneous shell. For a shell with general strain–energy function and inhomogeniety, the strain–energy function of the equivalent homogeneous material is determined explicitly. The resulting formula is used to study layered composite shells. The equivalent homogeneous material for an infinitely fine layered composite shell is examined, and is found to give not only the same global response, but also the same average stress field as the composite shell does.     

Optimal Time Decay of the Vlasov–Poisson–Boltzmann System in $${\mathbb R^3}$$

The Vlasov–Poisson–Boltzmann System governs the time evolution of the distribution function for dilute charged particles in the presence of a self-consistent electric potential force through the Poisson equation. In this paper, we are concerned with the rate of convergence of solutions to equilibrium for this system over \mathbb R³{\mathbb R^3}. It is shown that the electric field, which is indeed responsible for the lowest-order part in the energy space, reduces the speed of convergence, hence the dispersion of this system over the full space is slower than that of the Boltzmann equation without forces; the exact L ²-rate for the former is (1 + t)^−1/4 while it is (1 + t)^−3/4 for the latter. For the proof, in the linearized case with a given non-homogeneous source, Fourier analysis is employed to obtain time-decay properties of the solution operator. In the nonlinear case, the combination of the linearized results and the nonlinear energy estimates with the help of the proper Lyapunov-type inequalities leads to the optimal time-decay rate of perturbed solutions under some conditions on initial data.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 1: polycrystalline plasticity

A general set of flow laws and associated variational formulations are constructed for small-deformation rate-independent problems in strain-gradient plasticity. The framework is based on the thermodynamically consistent theory due to Gurtin and Anand (J Mech Phys Solids 53:1624–1649, 2005), and includes as variables a set of microstresses which have both energetic and dissipative components. The flow law is of associative type. It is expressed as a normality law with respect to a convex but otherwise arbitrary yield function, or equivalently in terms of the corresponding dissipation function. Two cases studied are, first, an extension of the classical Hill-Mises or J ₂ flow law and second, a form written as a linear sum of the magnitudes of the plastic strain and strain gradient. This latter form is motivated by work of Evans and Hutchinson (Acta Mater 57:1675–1688, 2009) and Nix and Gao (J Mech Phys Solids 46:411–425, 1998), who show that it leads to superior correspondence with experimental results, at least for particular classes of problems. The corresponding yield function is obtained by a duality argument. The variational problem is based on the flow rule expressed in terms of the dissipation function, and the problem is formulated as a variational inequality in the displacement, plastic strain, and hardening parameter. Dissipative components of the microstresses, which are indeterminate, are absent from the formulation. Existence and uniqueness of solutions are investigated for the generalized Hill-Mises and linear-sum dissipation functions, and for various combinations of defect energy. The conditions for well-posedness of the problem depend critically on the choice of dissipation function, and on the presence or otherwise of a defect energy in the plastic strain or plastic strain gradient, and of internal-variable hardening.     

An Internal Damping Model for the Absolute Nodal Coordinate Formulation

Introducing internal damping in multibody system simulations is important as real-life systems usually exhibit this type of energy dissipation mechanism. When using an inertial coordinate method such as the absolute nodal coordinate formulation, damping forces must be carefully formulated in order not to damp rigid body motion, as both this and deformation are described by the same set of absolute nodal coordinates. This paper presents an internal damping model based on linear viscoelasticity for the absolute nodal coordinate formulation. A practical procedure for estimating the parameters that govern the dissipation of energy is proposed. The absence of energy dissipation under rigid body motion is demonstrated both analytically and numerically. Geometric nonlinearity is accounted for as deformations and deformation rates are evaluated by using the Green–Lagrange strain–displacement relationship. In addition, the resulting damping forces are functions of some constant matrices that can be calculated in advance, thereby avoiding the integration over the element volume each time the damping force vector is evaluated.     

The anisotropic and angularly inhomogeneous elastic wedge under a monomial load distribution

In this study, the generally anisotropic and angularly inhomogeneous wedge under a monomial type of distributed loading of order n of, the radial coordinate r at its external faces is considered. At first, using variable separable relations in the equilibrium equations, the strain–stress relations and the strain compatibility equation, a differential system of equations, is constructed. Decoupling this system, an ordinary differential equation is derived and the stress and displacement fields may be determined. The proposed procedure is also applied to the elastostatic problem of an isotropic and angularly inhomogeneous wedge. The special cases of loading of order n=−1 and n=−2, where the self-similarity approach is not valid, are examined and the stress and displacements fields are derived. Applications are presented for the cases of an angularly inhomogeneous wedge and in the case of a bi-material isotropic wedge.     

Numerical simulation of shock wave propagation in microchannels using continuum and kinetic approaches

Numerical simulations of shock wave propagation in microchannels and microtubes (viscous shock tube problem) have been performed using three different approaches: the Navier–Stokes equations with the velocity slip and temperature jump boundary conditions, the statistical Direct Simulation Monte Carlo method for the Boltzmann equation, and the model kinetic Bhatnagar–Gross–Krook equation with the Shakhov equilibrium distribution function. Effects of flow rarefaction and dissipation are investigated and the results obtained with different approaches are compared. A parametric study of the problem for different Knudsen numbers and initial shock strengths is carried out using the Navier–Stokes computations.      

Nonlinear dynamics of functionally graded pipes conveying hot fluid

For improved stability of fluid-conveying pipes operating under the thermal environment, functionally graded materials (FGMs) are recommended in a few recent studies. Besides this advantage, the nonlinear dynamics of fluid-conveying FG pipes is an important concern for their engineering applications. The present study is carried out in this direction, where the nonlinear dynamics of a vertical FG pipe conveying hot fluid is studied thoroughly. The FG pipe is considered with pinned ends while the internal hot fluid flows with the steady or pulsatile flow velocity. Based on the Euler–Bernoulli beam theory and the plug-flow model, the nonlinear governing equation of motion of the fluid-conveying FG pipe is derived in the form of the nonlinear integro-partial-differential equation that is subsequently reduced as the nonlinear temporal differential equation using Galerkin method. The solutions in the time or frequency domain are obtained by implementing the adaptive Runge–Kutta method or harmonic balance method. First, the divergence characteristics of the FG pipe are investigated and it is found that buckling of the FG pipe arises mainly because of temperature of the internal fluid. Next, the dynamic characteristics of the FG pipe corresponding to its pre- and post-buckled equilibrium states are studied. In the pre-buckled equilibrium state, higher-order parametric resonances are observed in addition to the principal primary and secondary parametric resonances, and thus the usual shape of the parametric instability region deviates. However, in the post-buckled equilibrium state of the FG pipe, its chaotic oscillations may arise through the intermittent transition route, cyclic-fold bifurcation, period-doubling bifurcation and subcritical bifurcation. The overall study reveals complex dynamics of the FG pipe with respect to some system parameters like temperature of fluid, material properties of FGM and fluid flow velocity.     

Construction of the stationary probability density for a family of SDOF strongly non-linear stochastic second-order dynamical systems

In this paper, a new procedure is proposed to construct the stationary probability density for a family of the single-degree-of-freedom (SDOF) strongly non-linear stochastic second-order dynamical systems subjected to parametric and/or external Gaussian white noises. First of all, the Fokker-Planck-Kolmogorov (FPK) equation associated with the original Itô stochastic differential equation is replaced by the equivalent FPK equation by adding arbitrary anti-symmetric diffusion coefficient. Then, a family of invariant measures depending on the arbitrary anti-symmetric diffusion coefficient and another arbitrary function is constructed by vanishing the probability flows in two directions. Finally, the drift vector associated with a family of Itô stochastic differential equations is deduced by giving, a priori, these two arbitrary functions. It is shown that the known invariant measures dependent on energy are only the special cases of invariant measures presented in this paper, while some other classes of invariant measures are independent of energy. Thus, the invariant measures constructed in this paper are those belonging to the most general class of the SDOF strongly non-linear stochastic second-order dynamical systems so far.     

Homogenization of scalar wave equations with hysteresis

The paper deals with a scalar wave equation of the form where is a Prandtl–Ishlinskii operator and are given functions. This equation describes longitudinal vibrations of an elastoplastic rod. The mass density and the Prandtl–Ishlinskii distribution function are allowed to depend on the space variable x. We prove existence, uniqueness and regularity of solution to a corresponding initial-boundary value problem. The system is then homogenized by considering a sequence of equations of the above type with spatially periodic data and , where the spatial period tends to 0. We identify the homogenized limits and and prove the convergence of solutions to the solution of the homogenized equation. Received June 17, 1999     

Asymptotic Variational Wave Equations

We investigate the equation (u _t +(f(u)) _x ) _x =f ^{′ ′}(u) (u _x )²/2 where f(u) is a given smooth function. Typically f(u)=u ²/2 or u ³/3. This equation models unidirectional and weakly nonlinear waves for the variational wave equation u _tt − c(u) (c(u)u _x ) _x =0 which models some liquid crystals with a natural sinusoidal c. The equation itself is also the Euler–Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view.We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function f has a Lipschitz continuous second-order derivative. In the case where f is convex, the Cauchy problem is well posed also within the class of dissipative solutions. However, when f is not convex, we show that the dissipative solutions do not depend continuously on the initial data.     

On solutions of the problem in stresses with the use of Maxwell stress functions

The spatial problems of elasticity are mainly solved in displacements [1, 2], i.e., the Lamé equations are taken as the initial equations. This is related to the lack of general solutions for the system of basic equations of elasticity expressed in stresses. In this connection, a new variational statement of the problem in stresses was developed in [3, 4]; this statement consists in solving six generalized equations of compatibility for six independent components of the stress tensor, while the three equilibrium equations are transferred to the set of boundary conditions. This method is more convenient for the numerical solution of problems in stresses and has been tested when solving various boundary value problems. In the present paper, analyzing the completeness of the Saint-Venant identities and using the Maxwell stress functions, we obtain a new resolving system of three differential equations of strain compatibility for the three desired stress functions φ, ξ, and ψ. This system is an alternative to the three Lamé equilibrium equations for three desired displacement components u, v, w and is simpler in structure. Moreover, both of these systems of resolving equations can be solved by the new recursive-operator method [5, 6]. In contrast to well-known methods for constructing general solutions of linear differential equations and their systems, the solutions obtained by the recursive-operator method are constructed as operator-power series acting on arbitrary analytic functions of real variables (not necessarily harmonic), and the series coefficients are determined from recursive relations (matrix in the case of systems of equations). The arbitrary functions contained in the general solution can be determined directly either from the boundary conditions (the obtained system of inhomogeneous equations with a right-hand side can also be solved by the recursive-operator method [6]) or by choosing them from various classes of analytic functions (elementary, special); a complete set of particular solutions can be obtained in the same function classes, and the coefficients of linear combinations of particular solutions can be determined by the Trefftz method, the least-squares method, and the collocation method.     

Asymptotic properties of solutions of an equation of order <Emphasis Type="Italic">n</Emphasis> with almost constant coefficients and pulse action

We investigate the asymptotic behavior of solutions of linear differential equations with almost constant coefficients and pulse action at fixed times as t tends to infinity. We establish conditions for the times of pulse action under which there exist values of pulse action for which the solution of the considered Cauchy problem with initial conditions that coincide with the initial conditions for a certain (arbitrary but fixed) solution of the original equation without pulse action is bounded, unbounded, or tending to infinity. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 4, pp. 444–455, October–December, 2005.     

Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

L 1 and BV-type stability of the inelastic Boltzmann equation near vacuum

The L ¹ and BV-type stability to mild solutions of the inelastic Boltzmann equation is given in this paper. The result is an extension of the stability of the classical solution of the elastic Boltzmann equation proved in Ha (Arch. Ration. Mech. Anal. 173:25–42, 2004 [16]). The observation relies on the energy loss of the inelastic Boltzmann equation. This is a continuity work of Alonso (Indiana Univ. Math. J. [1]), where the author obtained the global existence of a mild solution for the inelastic Boltzmann equation. The proof is based on the mollification method and constructing some functionals as the one in Chae and Ha (Contin. Mech. Thermodyn. 17(7):511–524, 2006 [9]).     

Instability of equilibrium of nonholonomic systems with dissipation and circulatory forces

The paper discusses the equilibrium instability problem of the scleronomic nonholonomic systems acted upon by dissipative, conservative, and circulatory forces. The method is based on the existence of solutions to the differential equations of the motion which asymptotically tends to the equilibrium state of the system as t tends to negative infinity. It is assumed that the kinetic energy, the Rayleigh dissipation function, and the positional forces in the neighborhood of the equilibrium position are infinitely differentiable functions. The results obtained here are partially generalized the results obtained by Kozlov et al. (Kozlov, V. V. The asymptotic motions of systems with dissipation. Journal of Applied Mathematics and Mechanics, 58^(5), 787–792 (1994). Merkin, D. R. Introduction to the Theory of the Stability of Motion (in Russian), Nauka, Moscow (1987). Thomson, W. and Tait, P. Treatise on Natural Philosophy, Part I, Cambridge University Press, Cambridge (1879)). The results are illustrated by an example.     

On the Impossibility, in Some Cases, of the Leray�CHopf Condition for Energy Estimates

Current proofs of time independent energy bounds for solutions of the time dependent Navier–Stokes equations, and of bounds for the Dirichlet norms of steady solutions, are dependent upon the construction of an extension of the prescribed boundary values into the domain that satisfies the inequality (1.1) below, for a value of κ less than the kinematic viscosity. It is known from the papers of Leray (J Math Pure Appl 12:1–82, 1993), Hopf (Math Ann 117:764–775, 1941) and Finn (Acta Math 105:197–244, 1961) that such a construction is always possible if the net flux of the boundary values across each individual component of the boundary is zero. On the other hand, the nonexistence of such an extension, for small values of κ, has been shown by Takeshita (Pac J Math 157:151–158, 1993) for any two or three-dimensional annular domain, when the boundary values have a net inflow toward the origin across each component of the boundary. Here, we prove a similar result for boundary values that have a net outflow away from the origin across each component of the boundary. The proof utilizes a class of test functions that can detect and measure deformation. It appears likely that much of our reasoning can be applied to other multiply connected domains.     

Symmetry classification of quasi-linear PDE’s containing arbitrary functions

We consider the problem of performing the preliminary “symmetry classification” of a class of quasi-linear PDE’s containing one or more arbitrary functions: we provide an easy condition involving these functions in order that nontrivial Lie point symmetries be admitted, and a “geometrical” characterization of the relevant system of equations determining these symmetries. Two detailed examples will elucidate the idea and the procedure: the first one concerns a nonlinear Laplace-type equation, the second a generalization of an equation (the Grad–Schlüter–Shafranov equation) which is used in magnetohydrodynamics.     

A class of exact solutions for finite plane strain deformations of a particular elastic material

Summary In this paper we consider plane deformations of an incompressible elastic material and we show that by a suitable choice of strain energy function we can find the class of deformations with constant local rotation angle. Although the form for the strain energy function is chosen in the first place for mathematical convenience it does correspond to physically reasonable behaviour and such a theory may be regarded as a first order theory. The class of solutions obtained are expressed in a parametric form involving an arbitrary function, simple choices of which correspond to the well known exact solutions of finite elasticity.