Torsion of Incompressible Fiber-Reinforced Nonlinearly Elastic Circular Cylinders

Torsion of solid cylinders in the context of nonlinear elasticity theory has been widely investigated with application to the behavior of rubber-like materials. More recently, this problem has attracted attention in investigations of the biomechanics of soft tissues and has been applied, for example, to examine the mechanical behavior of passive papillary muscles of the heart. A recent study in nonlinear elasticity was concerned specifically with the effects of strain-stiffening on the torsional response of solid circular cylinders. The cylinders are composed of incompressible isotropic nonlinearly elastic materials that undergo severe strain-stiffening in the stress-stretch response. Here we investigate similar issues for fiber-reinforced transversely-isotropic circular cylinders. We consider a class of incompressible anisotropic materials with strain-energy densities that are of logarithmic form in the anisotropic invariant. These models reflect stretch induced strain-stiffening of collagen fibers on loading and have been shown to model the mechanical behavior of many fibrous soft biological tissues. The consideration of anisotropy leads to a more elaborate mechanical response than was found for isotropic strain-stiffening materials. The classic Poynting effect found for rubber-like materials where torsion induces elongation of the cylinder is shown to be significantly different for the transversely-isotropic materials considered here. For sufficiently large anisotropy and under certain conditions on the amount of twist, a reverse-Poynting effect is demonstrated where the cylinder tends to shorten on twisting The results obtained here have important implications for the development of accurate torsion test protocols for determination of material properties of soft tissues.     

Lyapunov and LMI analysis and feedback control of border collision bifurcations

Feedback control of piecewise smooth discrete-time systems that undergo border collision bifurcations is considered. These bifurcations occur when a fixed point or a periodic orbit of a piecewise smooth system crosses or collides with the border between two regions of smooth operation as a system parameter is quasistatically varied. The class of systems studied is piecewise smooth maps that depend on a parameter, where the system dimension n can take any value. The goal of the control effort in this work is to replace the bifurcation so that in the closed-loop system, the steady state remains locally attracting and locally unique (“nonbifurcation with persistent stability”). To achieve this, Lyapunov and linear matrix inequality (LMI) techniques are used to derive a sufficient condition for nonbifurcation with persistent stability. The derived condition is stated in terms of LMIs. This condition is then used as a basis for the design of feedback controls to eliminate border collision bifurcations in piecewise smooth maps and to produce the desirable behavior noted earlier. Numerical examples that demonstrate the effectiveness of the proposed control techniques are given.     

End effects for anti-plane shear deformations of sandwich structures

The purpose of this research is to further investigate the effects of material inhomogeneity and the combined effects of material inhomogeneity and anisotropy on the decay of Saint-Venant end effects. Saint-Venant decay rates for self-equilibrated edge loads in symmetric sandwich structures are examined in the context of anti-plane shear for linear anisotropic elasticity. The problem is governed by a second-order, linear, elliptic, partial differential equation with discontinuous coefficients. The most general anisotropy consistent with a state of anti-plane shear is considered, as well as a variety of boundary conditions. Anti-plane or longitudinal shear deformations are one of the simplest classes of deformations in solid mechanics. The resulting deformations are completely characterized by a single out-of-plane displacement which depends only on the in-plane coordinates. They can be thought of as complementary deformations to those of plane elasticity. While these deformations have received little attention compared with the plane problems of linear elasticity, they have recently been investigated for anisotropic and inhomogeneous linear elasticity. In the context of linear elasticity, Saint-Venant's principle is used to show that self-equilibrated loads generate local stress effects that quickly decay away from the loaded end of a structure. For homogeneous isotropic linear elastic materials this is well-documented. Self-equilibrated loads are a class of load distributions that are statically equivalent to zero, i.e., have zero resultant force and moment. When Saint-Venant's principle is valid, pointwise boundary conditions can be replaced by more tractable resultant conditions. It is shown in the present study that material inhomogeneity significantly affects the practical application of Saint-Venant's principle to sandwich structures.     

General stress-strain relations in non-linear theory of viscoelasticity for large deformations

Conclusion General phenomenoligical stress-strain relations in non-linear theory of visco-elasticity for large deformations have been presented.In the first place, according to V. V. Novozhilov 1 we express the generalized equilibrium equation for large deformations in the Lagrange representation, and we apply the generalized Hamilton's principle to the equation of energy conservation, which denotes that the sum of the elastic energy and the dissipative energy is equal to the work done by the body force and the surface on the substance; so that we obtain the required general stress-strain relations in comparison with the above two equations.On the condition that the elastic potential is a function only of the strain, and the dissipation function is a function of the rate of strain and of strain; such a substance is reduced to the Voigt material necessarily, and the stresses which act on the substance are given by the sum of elastic- and viscous stresses, and the stress-strain relations are reduced to the so-called Lagrangian form.If elongations, shears and angles of rotation are small and also the strains and rates of strain are sufficiently small, the stress-strain relations are expressed by a linear Voigt model constituting a Hookian spring in parallel with a Newtonian dashpot.Non-linearity in the theory is classified into two groups i. e. the geometrical non-linearity and the physical non-linearity. The former is introduced into the theory through the definition of the generalized strain and of the generalized stress and through the equilibrium equation for large deformation, and the latter through the general stress-strain relations.The main result of this paper is that the general stress-strain relations in viscoelasticity are deduced necessarily from the physically appropriate assumptions.     

A study on yield and flow of orthotropic materials

On the assumption that the yield criterion of orthotropic materials is isomorphic with Huber-Mises criterion of isotropic materials, we put forward a dimensionless stress yield criterion, and obtained the associated plastic flow law. Using experimental stress-strain curves in various simple stress states, generalized effective stress-strain formulae may be derived correspondingly in various forms.     

Uniqueness and Complementary Energy in Nonlinear Elastostatics

Global uniqueness of the smooth stress and deformation to within the usual rigid-body translation and rotation is established in the null traction boundary value problem of nonlinear homogeneous elasticity on a n-dimensional star-shaped region. A complementary energy is postulated to be a function of the Biot stress and to be para-convex and rank-(n-1) convex, conditions analogous to quasi-convexity and rank-(n-2) of the stored energy function. Uniqueness follows immediately from an identity involving the complementary energy and the Piola-Kirchhoff stress. The interrelationship is discussed between the two conditions imposed on the complementary energy, and between these conditions and those known for uniqueness in the linear elastic traction boundary value problem.     

Construction of two-phase equilibria in a non-elliptic hyperelastic material

This work focuses on the construction of equilibrated two-phase antiplane shear deformations of a non-elliptic isotropic and incompressible hyperelastic material. It is shown that this material can sustain metastable, two-phase equilibria which are neither piecewise homogeneous nor axisymmetric, but, rather, involve non-planar interfaces which completely segregate inhomogeneously deformed material in distinct elliptic phases. These results are obtained by studying a constrained boundary value problem involving an interface across which the deformation gradient jumps. The boundary value problem is recast as an integral equation and conditions on the interface sufficient to guarantee the existence of a solution to this equation are obtained. The contraints, which enforce the segregation of material in the two elliptic phases, are then studied. Sufficient conditions for their satisfaction are also secured. These involve additional restrictions on the interface across which the deformation gradient jumps — which, with all restrictions satisfied, constitutes a phase boundary. An uncountably infinite number of such phase boundaries are shown to exist. It is demonstrated that, for each of these, there exists a solution — unique up to an additive constant — for the constrained boundary value problem. As an illustration, approximate solutions which correspond to a particular class of phase boundaries are then constructed. Finally, the kinetics and stability of an arbitrary element within this class of phase boundaries are analyzed in the context of a quasistatic motion.     

Plane Strain Bending of Cylindrical Sectors of Admissible Compressible Hyperelastic Materials

In the theory of nonlinear elasticity of rubber-like materials, if a homogeneous isotropic compressible material is described by a strain–energy function that is a homogeneous function of the principal stretches, then the equations of equilibrium for axisymmetric deformations reduce to a separable first-order ordinary differential equation. For a particular class of such strain–energy functions, this property is used to obtain a general parametric solution to the equilibrium equation for plane strain bending of cylindrical sectors. Specification of the arbitrary function that appears in such strain–energy functions yields some parametric solutions. In some cases, the parameter can be eliminated to yield closed-form solutions in implicit or explicit form. Other possible forms for the arbitrary constitutive function that are likely to yield such solutions are also indicated.     

A 3D finite volume scheme for solving the updated Lagrangian form of hyperelasticity

The finite volume discretization of nonlinear elasticity equations seems to be a promising alternative to the traditional finite element discretization as mentioned by Lee et al. [Computers and Structures (2013)]. In this work, we propose to solve the elastic response of a solid material by using a cell‐centered finite volume Lagrangian scheme in the current configuration. The hyperelastic approach is chosen for representing elastic isotropic materials. In this way, the constitutive law is based on the principle of frame indifference and thermodynamic consistency, which are imposed by mean of the Coleman–Noll procedure. It results in defining the Cauchy stress tensor as the derivative of the free energy with respect to the left Cauchy–Green tensor. Moreover, the materials being isotropic, the free‐energy is function of the left Cauchy–Green tensor invariants, which enable the use of the neo‐Hookean model. The hyperelasticity system is discretized using the cell‐centered Lagrangian scheme from the work of Maire et al. [J. Comput. Phys. (2009)]. The 3D scheme is first order in space and time and is assessed against three test cases with both infinitesimal displacements and large deformations to show the good accordance between the numerical solutions and the analytic ones. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite Third-order Gradient Elasticity and Thermoelasticity

A constitutive format for the third-order gradient elasticity is suggested. It includes both isotropic and anisotropic non-linear behavior under finite deformations. Appropriate invariant stress and strain variables are introduced, which allow for reduced forms of the elastic energy law that identically fulfill the objectivity requirement. After working out the transformation behavior under a change of the reference placement, the symmetry transformations for third-order materials can be introduced. After the mechanical third-order theory, an extension to thermoelasticity is given, and necessary and sufficient conditions are derived from the Clausius-Duhem inequality.     

Elucidating the escape dynamics of the four hill potential

This paper treads discontinuous bifurcation in piecewise smooth systems of Filippov type. These bifurcations occur when a fixed point or a periodic orbit crosses with the border between two regions of smooth behavior. A detailed analysis of generalization Poincaré map and monodromy matrix which are related shows that subfamily of system with invariant cone-like objects is foliated by periodic orbits and determines its stability. In addition, we introduce a theoretical framework for analyzing 3D perturbed nonlinear piecewise smooth systems and give necessary conditions so that different types of bifurcations occur. The analysis identifies criteria for the existence of a novel bifurcation based on sensitively the location of the return map. Moreover, the piecewise smooth Melnikov function and sufficient conditions of the existence of the periodic orbits for nonlinear perturbed system are explicitly obtained.     

Global Existence of Smooth Solutions for Partially Dissipative Hyperbolic Systems with a Convex Entropy

We consider the Cauchy problem for a general one-dimensional n×n hyperbolic symmetrizable system of balance laws. It is well known that, in many physical examples, for instance for the isentropic Euler system with damping, the dissipation due to the source term may prevent the shock formation, at least for smooth and small initial data. Our main goal is to find a set of general and realistic sufficient conditions to guarantee the global existence of smooth solutions, and possibly to investigate their asymptotic behavior. For systems which are entropy dissipative, a quite natural generalization of the Kawashima condition for hyperbolic-parabolic systems can be given. In this paper, we first propose a general framework for this kind of problem, by using the so-called entropy variables. Then we go on to prove some general statements about the global existence of smooth solutions, under different sets of conditions. In particular, the present approach is suitable for dealing with most of the physical examples of systems with a relaxation extension. Our main tools will be some refined energy estimates and the use of a suitable version of the Kawashima condition.     

Optimization of the stored energy and coaxiality of strain and stress in finite elasticity

The strain energy density of a hyperelastic anisotropic body which is rotated before being subjected to a given but arbitrary deformation is viewed as a smooth function defined on the group of rotations, parametrized by the deformation gradient. It is shown that the critical points of this function correspond to rotations which, when composed with the prescribed deformation, yield a total strain tensor which is coaxial with the corresponding stress. For any type of material symmetry, there are at least two such rotations. Coaxiality of stress and strain for all deformations is shown to be a sufficient condition for the isotropicity of hyperelastic materials.Research supported by GNFM of CNR (Italy).     

Existence of Compressible Current-Vortex Sheets: Variable Coefficients Linear Analysis

We study the initial-boundary value problem resulting from the linearization of the equations of ideal compressible magnetohydrodynamics and the Rankine-Hugoniot relations about an unsteady piecewise smooth solution. This solution is supposed to be a classical solution of the system of magnetohydrodynamics on either side of a surface of tangential discontinuity (current-vortex sheet). Under some assumptions on the unperturbed flow, we prove an energy a priori estimate for the linearized problem. Since the tangential discontinuity is characteristic, the functional setting is provided by the anisotropic weighted Sobolev space W₂^1,σ. Despite the fact that the constant coefficients linearized problem does not meet the uniform Kreiss-Lopatinskii condition, the estimate we obtain is without loss of smoothness even for the variable coefficients problem and nonplanar current-vortex sheets. The result of this paper is a necessary step in proving the local-in-time existence of current-vortex sheet solutions of the nonlinear equations of magnetohydrodynamics.     

The Somigliana identity on piecewise smooth surfaces

Summary TheSomigliana identity is an integral representation of the displacement field of a body, i.e. the solution of the Cauchy-Navier equation in the classical theory of elasticity is represented by influence functions. This well established identity for bodies with smooth surfaces is extended in the present paper to cover two- and three-dimensional bodies with piecewise smooth surfaces.

---

Zusammenfassung DieSomigliana Identität ist die Integraldarstellung des Verschiebungsfeldes eines Körpers, d.h. die Lösung der Cauchy-Navier'schen DGL in der klassischen Elastizitätstheorie wird durch Einflußfunktionen dargestellt. Diese Identität, die für Körper mit glatter Oberfläche wohl bekannt ist, wird in diesem Anfsatz auf Körper und Scheiben mit stückweise glatten Rändern ausgedehnt.

     

Non-linear stress-strain equations for incompressible transversely isotropic elastic materials

Non-linear stress-strain equations for incompressible, transversely isotropic elastic materials are developed. In order to obtain these equations, the expressions for a strain energy function is found. The derivation of the strain energy function follows a geometrical approach and a method suggested by Mooney. These stress-strain relations are expressed in terms of three principal stretches to the sixth order.     

Well-posedness of the fundamental boundary value problems for constrained anisotropic elastic materials

We consider the equations of linear homogeneous anisotropic elasticity admitting the possibility that the material is internally constrained, and formulate a simple necessary and sufficient condition for the fundamental boundary value problems to be well-posed. For materials fulfilling the condition, we establish continuous dependence of the displacement and stress on the elastic moduli and ellipticity of the elasticity system. As an application we determine the orthotropic materials for which the fundamental problems are well-posed in terms of their Young's moduli, shear moduli, and Poisson ratios. Finally, we derive a reformulation of the elasticity system that is valid for both constrained and unconstrained materials and involves only one scalar unknown in addition to the displacements. For a two-dimensional constrained material a further reduction to a single scalar equation is outlined.This paper is dedicated to Professor Joachim Nitsche on the occasion of his sixtieth birthday     

A new weighting method for improving the WENO‐Z scheme

The calculation of the weight of each substencil is very important for a weighted essentially nonoscillatory (WENO) scheme to obtain high‐order accuracy in smooth regions and keep the essentially nonoscillatory property near discontinuities. The weighting function introduced in the WENO‐Z scheme provides a straightforward method to analyze the accuracy order in smooth regions. In this paper, we construct a new sixth‐order global smoothness indicator (GSI‐6) and a function about GSI‐6 and the local smoothness indicators (IS_k) to calculate the weights. The analysis and numerical results show that, with the new weights, the scheme satisfies the sufficient condition for the fifth‐order convergence in smooth regions even at critical points. Meanwhile, it can also maintain low dissipation for discontinuous solutions due to relative large weights assigned to discontinuous substencils.