(A,B)-Invariant Polyhedral Sets of Linear Discrete-Time Systems

The problem of confining the trajectory of a linear discrete-time system in a given polyhedral domain is addressed through the concept of (A, B)-invariance. First, an explicit characterization of (A, B)-invariance of convex polyhedra is proposed. Such characterization amounts to necessary and sufficient conditions in the form of linear matrix relations and presents two major advantages compared to the ones found in the literature: it applies to any convex polyhedron and does not require the computation of vertices. Such advantages are felt particularly in the computation of the supremal (A, B)-invariant set included in a given polyhedron, for which a numerical method is proposed. The problem of computing a control law which forces the system trajectories to evolve inside an (A, B)-invariant polyhedron is treated as well. Finally, the (A, B)-invariance relations are generalized to persistently disturbed systems.     

Orbits of SU(2)-representations and minimal isometric immersions

In contrast to all known examples, we show that in the case of minimal isometric immersions of into the smallest target dimension is almost never achieved by an -equivariant immersion. We also give new criteria for linear rigidity of a fixed minimal isometric immersion of into . The minimal isometric immersions arising from irreducible SU(2)-representations are linearly rigid within the moduli space of SU(2)-equivariant immersions. Hence the question arose whether they are still linearly rigid within the full moduli space. We show that this is false by using our new criteria to construct an explicit SU(2)-equivariant immersion which is not linearly rigid. Various authors [GT], [To3], [W1] have shown that minimal isometric immersions of higher isotropy order play an important role in the study of the moduli space of all minimal isometric immersions of into . Using a new necessary and sufficient condition for immersions of isotropy order , we derive a general existence theorem of such immersions. Received: 13 May 1999 / in final form: 13 July 1999     

Construction of polyconvex,anisotropic free‐energy functions

The existence of minimizers of some variational principles in finite elasticity is based on the concept of quasiconvexity, introduced by Morrey [6]. This integral inequality is rather complicated to handle. Thus, the sufficient condition for quasiconvexity, the polyconvexity condition in the sense of Ball [1], is a more important concept for practical applications, see also Ciarlet [4] and Dacorogna [5]. In the case of isotropy there exist some models which satisfy this condition. Furthermore, there does not exist a systematic treatment of anisotropic, polyconvex free‐energies in the literature. In the present work we discuss some aspects of the formulation of polyconvex, anisotropic free‐energy functions in the framework of the invariant formulation of anisotropic constitutive equations and focus on transverse isotropy.     

Analytic solutions of a reducible strain gradient elasticity model for solid cylinder with a cavity and its application in zonal failure

In this work, the reducibility condition of the fourth-order equilibrium equation in the strain gradient elasticity (SGE) model for solid cylinder with a cavity is obtained. When the reducibility condition is satisfied, the analytic displacement, generalized radial stress, and generalized angular stress can be solved out, and according to the higher-order coefficients, internal length scale, and Lamé constants, the displacement, generalized radial stress, and generalized angular stress are classified into four types: (1) conventional elasticity solution, (2) quasiperiodic SGE solution, (3) monotonous SGE solution, and (4) non-real-number solution. Quasiperiodic generalized radial stress and generalized angular stress are used to explain the occurrence of zonal failure of surrounding rock of a circular roadway. Numerical analysis with MATLAB is applied to study the influence of loading on zonal failure of surrounding rock of a circular raodway.     

Love waves of SH type in an inhomogeneous transversely isotropic elastic medium

The asymptotics of high-frequency Love waves, which are analogous to transverse surface SH waves, is considered for a special type of anisotropy (transverse isotropy) of elastic media. The wave field is represented as a sum of the space-time (ST) caustic expansion and two additional ST ray series for faster (relative to the transverse surface wave) body waves, decaying exponentially with depth. Near the surface, the coefficients of the ST caustic and ray series, as well as the eikonals of waves, are determined in the form of expansions in a small parameter, which characterizes the proximity of the caustic of the ray field to the surface. With regard for the specific structure of the elasticity tensor of a transversely isotropic medium, the surface is treated as a plane. Interrelations between the parameters of elasticity, which are consistent with the conditions of the positivity of the elastic deformation energy and provide for the origination of the surface waves considered, are obtained.Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 239, 1997, pp. 243–262.This work was supported by the Russian Foundation for Basic Research under grant No. 96-01-00666.     

Almost Kähler 4-dimensional Lie groups with J-invariant Ricci tensor

The J-invariance of the Ricci tensor is a natural weakening of the Einstein condition in almost Hermitian geometry. The aim of this paper is to determine left-invariant strictly almost Kähler structures (g,J,Ω) on real 4-dimensional Lie groups such that the Ricci tensor is J-invariant. We prove that all these Lie groups are isometric (up to homothety) to the (unique) 4-dimensional proper 3-symmetric space.     

Linear invariants of a Cartesian tensor

The number of linear invariants under SO(3) as well as SO(2)of a Cartesian tensor of an arbitrary rank is studied. A linearform is defined in terms of elements of a tensor. It is establishedthat the number of linear invariants of a tensor of rank n underSO(3) equals the dimension of the space of isotropic tensorsof rank n. Formulas for the number of invariants in the twocases are also derived. For the elasticity tensor, our analysisconfirms the results of Norris.     

The Hencky strain energy ‖log U‖2 measures the geodesic distance of the deformation gradient to SO(n) in the canonical left-invariant Riemannian metric on GL(n)

The well-known isotropic Hencky strain energy appears naturally as a distance measure of the deformation gradient to the set of rigid rotations in a canonical left-invariant Riemannian metric on the general linear group GL(n). Objectivity requires the Riemannian metric to be left-GL(n) invariant, isotropy requires the Riemannian metric to be right-O(n) invariant. The latter two conditions are satisfied for a three-parameter family of Riemannian metrics on the tangent space of GL(n). Surprisingly, the final result is basically independent of the chosen parameters. In deriving the result, geodesics on GL(n) have to be parametrized and a novel minimization problem, involving the matrix logarithm for non-symmetric arguments, has to be solved. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

New Homogeneous Einstein Metrics on SO(7)/T

The authors compute non-zero structure constants of the full flag manifold M = SO(7)/T with nine isotropy summands, then construct the Einstein equations. With the help of computer they get all the forty-eight positive solutions (up to a scale ) for SO(7)/T, up to isometry there are only five G-invariant Einstein metrics, of which one is Kähler Einstein metric and four are non-Kähler Einstein metrics.     

Computational homogenization of materials with small deformation to determine configurational forces

In this work the mechanical boundary condition for the micro problem in a two-scaled homogenization using a FE² approach is discussed. The strain tensor is often used in the literature for small deformation problem to determine the boundary conditions for the boundary value problem on the micro level. This strain tensor based boundary condition gives consistent homogenized mechanical quantities, e.g. stress tensor and elasticity tensor, but the present work points out that it leads to unphysical homogenized configurational forces. Instead, we propose a displacement gradient based boundary condition for the micro problem. Results show that the displacement gradient based boundary condition can give not only the consistent homogenized mechanical quantities but also the appropriate homogenized configurational forces. The interpretation of the displacement gradient based boundary condition is discussed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Taylor-Galerkin algorithms for viscoelastic flow: Application to a model problem

Coupled and decoupled Taylor-Galerkin algorithms are considered for viscoelastic flow and a model problem—transient startup Poiseuille flow in a channel under a fixed pressure gradient. All algorithms reproduce the steady-state solutions and are stable at high elasticity numbers (E). For a fixed mesh, the coupled and decoupled versions (TGC and TGD) give exceptional time-accuracy at low elasticity numbers [to within O(1%) at E = 1] and reasonable accuracy at high elasticity numbers [to within O(10%) at E = 10, 100]. By definition, the decoupled false-transient scheme (TGF), which uses different time scales for velocity and stress time stepping, provides a poor transient history. Where the main requirement is to compute a steady-state algorithm efficiency is crucial. The TGF scheme attains a steady state between six to eight times faster than does the TGC scheme, and the latter is over twice as fast as the TGD form. © 1994 John Wiley & Sons, Inc.     

The Euler-Poincaré-Hopf theorem for flat connections in some transitive Lie algebroids

This paper is a continuation of [19], [21], [22]. We study flat connections with isolated singularities in some transitive Lie algebroids for which either ℝ or sl(2, ℝ) or so(3) are isotropy Lie algebras. Under the assumption that the dimension of the isotropy Lie algebra is equal to n + 1, where n is the dimension of the base manifold, we assign to any such isolated singularity a real number called an index. For ℝ-Lie algebroids, this index cannot be an integer. We prove the index theorem (the Euler-Poincaré-Hopf theorem for flat connections) saying that the index sum is independent of the choice of a connection. Multiplying this index sum by the orientation class of M, we get the Euler class of this Lie algebroid. Some integral formulae for indices are given.     

Pointwise superconvergence of recovered gradients for piecewise linear finite element approximations to problems of planar linear elasticity

In this paper we consider piecewise linear finite element approximations to the problem of planar linear elasticity, and present some new estimates for the pointwise (L_∞) superconvergence of a recovered gradient function to the gradient of the true solution. This extends to linear elasticity the previous work of the present and other authors on L_∞ results for Poisson problems, and at the same time, to the L_∞ norm the previous L₂ results of the authors for linear elasticity.     

FETI-DP Methods for P-Elasticity

A FETI–DP method is introduced for the problem of linear P–elasticity which arises from linear elasticity by the introduction of a matrix P and which is motivated by micromorphic models. Numerical results as well as a condition number estimate are presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Smoothness of holonomy covers for singular foliations and essential isotropy

Given a singular foliation, we attach an “essential isotropy” group to each of its leaves, and show that its discreteness is the integrability obstruction of a natural Lie algebroid over the leaf. We show that a condition ensuring discreteness is the local surjectivity of a transversal exponential map associated with the maximal ideal of vector fields prescribed to be tangent to the foliation. The essential isotropy group is also shown to control the smoothness of the holonomy cover of the leaf (the associated fiber of the holonomy groupoid), as well as the smoothness of the associated isotropy group. Namely, the (topological) closeness of the essential isotropy group is a necessary and sufficient condition for the holonomy cover to be a smooth (finite-dimensional) manifold and the isotropy group to be a Lie group. These results are useful towards understanding the normal form of a singular foliation around a compact leaf. At the end of this article we briefly outline work of ours on this normal form, to be presented in a subsequent paper.     

On parallel solution of linear elasticity problems. Part II: Methods and some computer experiments

This is the second part of a trilogy on parallel solution of the linear elasticity problem. We consider the plain case of the problem with isotropic material, including discontinuous coefficients, and with homogeneous Dirichlet boundary condition. The discretized problem is solved by the preconditioned conjugate gradient (pcg) method. In the first part of the trilogy block‐diagonal preconditioners based on the separate displacement component part of the elasticity equations were analysed. The preconditioning systems were solved by the pcg‐method, i.e. inner iterations were performed. As preconditioner, we used modified incomplete factorization MIC(0), where possibly the element matrices were modified in order to give M‐matrices, i.e. in order to guarantee the existence of the MIC(0) factorization. In the present paper, the second part, full block incomplete factorization preconditioners are presented and analysed. In order to avoid inner/outer iterations we also study a variant of the block‐diagonal method and of the full block method, where the matrices of the inner systems are just replaced by their MIC(0)‐factors. A comparison is made between the various methods with respect to rate of convergence and work per unknown. The fastest methods are implemented by message passing utilizing the MPI system. In the third part of the trilogy, we will focus on the use of higher‐order finite elements. Copyright © 2002 John Wiley & Sons, Ltd.     

On non-isotropic harmonic maps of surfaces into complex projective spaces

It is proved by purely algebraic method that weakly conformai, conformai andA _z ³ = 0 are mutually equivalent if ϕ :Ω→ℂP ⁿ is a non-isotropic harmonic map and the harmonic maps with isotropy order ≥3 are uniquely determined by a system of ordinary differential equations. A method is given, by which the isotropy orders of non-isotropic harmonic maps can be computed.     

满旗流形SO(8)=T上不变爱因斯坦度量

本文研究了迷向表示分为12个不可约子空间的满旗流形SO(8)=T上不变爱因斯坦度量的问题.利用计算机计算满旗流形SO(8)=T爱因斯坦方程组的方法, 得到了满旗流形SO(8)=T上有160 个不变爱因斯坦度量(up to a scale)的结果, 在等距情况下考虑这160个不变爱因斯坦度量, 其中1个是凯莱爱因斯坦度量, 4 个是非凯莱爱因斯坦度量. 推广了只对迷向表示分为小于等于6个不可约子空间的满旗流形上不变爱因斯坦度量的研究.     

Planar central configurations as fixed points

Planar central configurations can be seen as critical points of the reduced potential or solutions of a system of equations. By the homogeneity of the potential and its O(2)-invariance it is possible to see that the SO(2)- orbits of central configurations are fixed points of a map f. The purpose of the paper is to define and study this map and to derive some properties using topological fixed point theory. The generalized Moulton–Smale theorem for collinear configurations is proved, together with some estimates on the number of central configurations in the case of three bodies, using fixed point indices. Well-known results such as the compactness of the set of central configurations follow easily in this topological framework. Dedicated to Professor Albrecht Dold and Professor Edward Fadell