Kinematics of finite-strain elastic-inelastic deformation

Polar decomposition tensors are constructed for slightly disturbed kinematic elastic, inelastic, and thermal strain tensors. Provided that the inelastic and thermal site gradients are pure deformations without rotations, relations are obtained between inelastic small strains and small rotations and between thermal small strains and small rotations which transform an intermediate configuration to a close current configuration. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 165–172, January–February, 2008.     

On some general variational principles for creep with applications to thin shells

A constitutive relation for a viscous material subject to small strains but finite rotations is postulated and associated variational theorems are formulated. These are similar to the principle of minimum of potential energy and the Hellinger-Reissner theorem of an elastic solid. The derivation of strain-displacement relations for thin shells subject to small strains but moderately large rotations are given. On this basis a mixed variational principle for thin viscous shells is developed. For the problem of creep collapse of long cylindrical shells under external pressure it is demonstrated that the mixed variational principle may be advantageous compared to other variational theorems. A comparison with the creep collapse theory of Hoff et al. is given.     

Constitutive Equations of an Isotropic Hyperelastic Body

Stress—strain equations for an isotropic hyperelastic body are formulated. It is shown that the strain energy density whose gradient determines stresses can be defined as a function of two rather than three arguments, namely, strain–tensor invariants. In the case of small strains, the equations become relations of Hooke's law with two material constants, namely, shear modulus and bulk modulus.     

Kinematically exact curved and twisted strain-based beam

The paper presents a formulation of the geometrically exact three-dimensional beam theory where the shape functions of three-dimensional rotations are obtained from strains by the analytical solution of kinematic equations. In general it is very demanding to obtain rotations from known rotational strains. In the paper we limit our studies to the constant strain field along the element. The relation between the total three-dimensional rotations and the rotational strains is complicated even when a constant strain field is assumed. The analytical solution for the rotation matrix is for constant rotational strains expressed by the matrix exponential. Despite the analytical relationship between rotations and rotational strains, the governing equations of the beam are in general too demanding to be solved analytically. A finite-element strain-based formulation is presented in which numerical integration in governing equations and their variations is completely omitted and replaced by analytical integrals. Some interesting connections between quantities and non-linear expressions of the beam are revealed. These relations can also serve as useful guidelines in the development of new finite elements, especially in the choice of suitable shape functions.     

Localized Coherent Structures in the Boundary Layer

A Blasius laminar boundary layer and a steady turbulent boundary layer on a flat plate in an incompressible fluid are considered. The spectral characteristics of the Tollmien—Schlichting (TS) and Squire waves are numerically determined in a wide range of Reynolds numbers. Based on the spectral characteristics, relations determining the three–wave resonance of TS waves are studied. It is shown that the three–wave resonance is responsible for the appearance of a continuous low–frequency spectrum in the laminar region of the boundary layer. The spectral characteristics allow one to obtain quantities that enter the equations of dynamics of localized perturbations. By analogy with the laminar boundary layer, the three–wave resonance of TS waves in a turbulent boundary layer is considered.     

Stability of deformation of isotropic hyperelastic bodies

The equations relating stress rates to strain rates are formulated and conditions of stable deformation of isotropic hyperelastic bodies are described. Stress–strain relations are presented for pure shear and uniaxial and axisymmetric loading of a material with a constitutive function obtained by generalization of the constitutive function of Hooke's law. In the case of small strains, the diagrams virtually coincide with the linear diagrams following from Hooke's law. Ramification of solutions and transition to declining diagrams begin at the same time, irrespective of values of the constants of the material, when large stresses of the order of the shear modulus are reached.     

Generators of the product of two Schoenflies motion groups

Schoenflies motion is often termed X-motion for conciseness. A set of X-motions with a given direction of its axes of rotations has the algebraic properties of a Lie group for the composition product of rigid-body motions or displacements. The product of two X-subgroups, which is the mathematical model of a serial concatenation of two kinematic chains generating two distinct X-motions, characterizes a noteworthy type of 5-dimensional (5D) displacement set called double Schoenflies motion or X–X motion for brevity. This X–X motion set is a 5D submanifold of the displacement 6D Lie group. Such a motion type includes any spatial translation (3T) and any two sequential rotations (2R) provided that the axes of rotation are parallel to two fixed independent vectors. This motion set also contains the rotations that are products of the foregoing two rotations. In the paper, some preliminary fundamentals on the 4D X-motion are recalled; the 5D set of X–X motions is emphasized. Then implementing serial arrays of one-dof Reuleaux pairs and hinged parallelograms, we enumerate all serial mechanical generators of X–X motion, which have no redundant internal mobility. Based on the group-theoretic concepts, one can differentiate two families of irreducible representations of an X–X motion. One family is realized by twenty-one open chains including the doubly planar motion generators as special cases. The other is generally classified into eight major categories in which one hundred and six distinct open chains generating X–X motion are revealed and nineteen more ones having at least one parallelogram are derived from them. Meanwhile, these kinematic chains are graphically displayed for a possible use in the structural synthesis of parallel manipulators.     

Nonaxisymmetric Thermoelastoplastic Analysis of Branched Shells of Revolution by the Semianalytic Finite-Element Method

A technique is proposed to solve elastoplastic deformation problems for branched shells of revolution under the action of asymmetric forces and a temperature field. The kinematic equations are derived within the framework of the linear Kirchhoff–Love theory of shells and the thermoelastic relations within the framework of the theory of small elastoplastic strains. The problem is given a variational formulation based on the virtual-displacement principle and the Fourier-series expansion of the unknown functions and loads with respect to the circumferential coordinate. The additional-load method is used to solve a nonlinear problem and the finite-elements method is used to carry out a numerical analysis. As an example, an asymmetric stress–strain analysis is performed for a cylindrical shell reinforced by a ring plate.     

Mathematical construction of a Reissner–Mindlin plate theory for composite laminates

A Reissner–Mindlin theory for composite laminates without invoking ad hoc kinematic assumptions is constructed using the variational-asymptotic method. Instead of assuming a priori the distribution of three-dimensional displacements in terms of two-dimensional plate displacements as what is usually done in typical plate theories, an exact intrinsic formulation has been achieved by introducing unknown three-dimensional warping functions. Then the variational-asymptotic method is applied to systematically decouple the original three-dimensional problem into a one-dimensional through-the-thickness analysis and a two-dimensional plate analysis. The resulting theory is an equivalent single-layer Reissner–Mindlin theory with an excellent accuracy comparable to that of higher-order, layer-wise theories. The present work is extended from the previous theory developed by the writer and his co-workers with several sizable contributions: (a) six more constants (33 in total) are introduced to allow maximum freedom to transform the asymptotically correct energy into a Reissner–Mindlin model; (b) the semi-definite programming technique is used to seek the optimum Reissner–Mindlin model. Furthermore, it is proved the first time that the recovered three-dimensional quantities exactly satisfy the continuity conditions on the interface between different layers and traction boundary conditions on the bottom and top surfaces. It is also shown that two of the equilibrium equations of three-dimensional elasticity can be satisfied asymptotically, and the third one can be satisfied approximately so that the difference between the Reissner–Mindlin model and the second-order asymptotical model can be minimized. Numerical examples are presented to compare with the exact solution as well as the classical lamination theory and the first-order shear-deformation theory, demonstrating that the present theory has an excellent agreement with the exact solution.     

On correct account of finite rotations in finite plasticity theory

The non-uniqueness of the trantition from nonobjective constitutive relations to objective ones with the use of the principle of material frame-indifference (PMFI) is shown. To eliminate it, the concept of finite strain without rotations (FSWR) for a given material type and each strain component (elastic, plastic) is introduced. In FSWR the rotation is excluded with respect to the natural preferred configuration for a given material. Considered are a simple solid, a liquid, a monocrystal, a polycrystal and a composite. The proecedure is proposed for consistent generalization of known infinitesimal relations for finite strains and rotations. The structure of constitutive relations is derived for anisotropic elasto-plastic mono- and polycrystalline materials.     

A study of elasticity and plasticity equations under arbitrary displacements and strains

Equations of geometrically nonlinear theory of elasticity with finite displacements and strains are analyzed. The equations are composed using three versions of physical relations and applied to solve the problem of tension-compression of a straight bar. It is shown that the use of the classical relations between the components of the stress tensor and the Cauchy-Green strain tensor in the problem of compression of the bar results in the appearance of “spurious” static loss of stability such that the bar axis remains straight if the stresses are referred to unit areas before the deformation (conditional stresses). However, in the problem of tension, the classical relations do not permit one to describe the phenomenon of static instability (neck formation as the plastic instability occurs). These drawbacks disappear if one uses the third version of the physical equations, composed as relations between the true stresses referred to unit areas of the deformed faces on which they act and the true elongations and shears. The relations of the third version are most correct; they permit one to pass to self-consistent equations of elasticity and plasticity under small strains and finite displacements, and they should be recommended for practical use. As an example, such relations are composed for the flow theory.     

Nonlinear surface waves in media with complex rheology

Weak nonlinear waves in a generalized viscoelastic medium with internal oscillators are considered. The rheological relations contain higher time derivatives of the stresses and strains as well as their tensor products. The method of expansion in a small parameter with the introduction of slow time and a running space coordinate is employed. The first approximation gives wave velocities and relations between the parameters equivalent to the results of an acoustic analysis at elastic wave fronts [1]. The second approximation leads to an evolution equation for the displacement velocity. For this a Fourier-Laplace double integral transformation is used. Reversion to the inverse transforms of the unknown functions leads to an integrodifferential evolution equation, which contains a Hubert transform and is a generalization of the Benjamin-Ono equation of deep water theory.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 95–103, September–October, 1990.     

Elastoplastic state of flexible cylindrical shells with two circular holes

The elastoplastic state of thin cylindrical shells with two equal circular holes is analyzed with allowance made for finite deflections. The shells are made of an isotropic homogeneous material. The load is internal pressure of given intensity. The distribution of stresses along the hole boundary and in the stress concentration zone (when holes are closely spaced) is analyzed by solving doubly nonlinear boundary-value problems. The results obtained are compared with the solutions that allow either for physical nonlinearity (plastic strains) or geometrical nonlinearity (finite deflections) and with the numerical solution of the linearly elastic problem. The stresses near the holes are analyzed for different distances between the holes and nonlinear factors.Translated from Prikladnaya Mekhanika, Vol. 40, No. 10, pp. 107–112, October 2004.     

Non-linear theories of beams and plates accounting for moderate rotations and material length scales

The primary objective of this paper is to formulate the governing equations of shear deformable beams and plates that account for moderate rotations and microstructural material length scales. This is done using two different approaches: (1) a modified von Kármán non-linear theory with modified couple stress model and (2) a gradient elasticity theory of fully constrained finitely deforming hyperelastic cosserat continuum where the directors are constrained to rotate with the body rotation. Such theories would be useful in determining the response of elastic continua, for example, consisting of embedded stiff short fibers or inclusions and that accounts for certain longer range interactions. Unlike a conventional approach based on postulating additional balance laws or ad hoc addition of terms to the strain energy functional, the approaches presented here extend existing ideas to thermodynamically consistent models. Two major ideas introduced are: (1) inclusion of the same order terms in the strain–displacement relations as those in the conventional von Kármán non-linear strains and (2) the use of the polar decomposition theorem as a constraint and a representation for finite rotations in terms of displacement gradients for large deformation beam and plate theories. Classical couple stress theory is recovered for small strains from the ideas expressed in (1) and (2). As a part of this development, an overview of Eringen׳s non-local, Mindlin׳s modified couple stress theory, and the gradient elasticity theory of Srinivasa–Reddy is presented.     

Interface crack with a contact zone in an isotropic bimaterial under thermomechanical loading

A plane problem for a thermally insulated interface crack with a contact zone in an isotropic bimaterial under tension–shear mechanical loading and a temperature flux is considered. The expressions for the stresses and the electrical flux as well as for the derivatives of the displacement and the temperature jumps at the material interfaces via sectionally holomorphic mechanical and thermal potential functions are given. After the solution of the thermal problem the inhomogeneous combined Dirichlet–Riemann boundary value problem is formulated and solved exactly. The stresses at the interface and the stress intensity factors at the singular points are presented in a clear analytical form. Special attention is devoted to the case of a small contact zone when the stress intensity factors can be presented in form similar to the associated presentation for an “open” crack model. A transcendental equation and an asymptotic analytic formula for the determination of the real contact zone length are derived. It is shown that for a certain bimaterial this length as well as the correspondent stress intensity factor are defined by a single parameter which depends on the normal-shear loading and the heat flux.     

Analysis of bifurcation bending modes of arches and panels

This paper reports the results of numerical analysis of the bifurcation solutions of nonlinear boundary-value problems of plane bending of elastic arches and panels. Problems are formulated for a system of six nonlinear ordinary differential equations of the first order with independent fields of finite displacements and rotations. Two loading versions (by follower and conservative pressures) and two versions of boundary conditions (rigid clamping and pinning) are considered. In the case of clamped arches and panels, the set of solutions consists of symmetric and asymmetric bending modes which exist only for positive values of the load parameter. In the case of pinning, the set of solutions includes symmetric and asymmetrical modes which correspond to positive, negative, and zero values of the parameter. In both problems, the phase relations between the state parameter and the load parameter are bifurcated, ambiguous, have isolated branches, and admit a catastrophe — a finite jump from the fundamental equilibrium mode to a buckled mode.     

Direct measurement of static yield properties of cohesive suspensions

A technique of yield stress investigation based upon the combined use of two devices (an applied stress rheometer and an instrument for measuring the propagation velocity of small amplitude, torsional shear waves) is described. Investigations into the low shear rate rheological properties of illitic suspensions are reported for shear rates, typically, in the range 10^–4— 10^–1 s^–1 under applied stresses in the range 0.01 — 10 Nm^–2 and involving shear strains between 10^–1 and 10^–4. Results are presented which demonstrate that the technique does not invoke the excessive structural disruption of material associated with applied shear rate based methods (direct and otherwise) and the widely encountered problem of wall slip at the surface of rotational measuring devices is avoided using miniature vane geometries. Results are compared with those obtained using smooth-walled cyclindrical measuring devices in both applied stress and applied shear rate instruments.Yield measurements are considered in relation to the structural properties of the undisturbed material state and shear moduli obtained by studying the propagation of small amplitude (10^–5 rad), high frequency (~ 300 Hz) torsional shear waves through the test materials are reported. Experimental techniques and instrument modifications to permit these measurements are described.     

The effect of impurities on the viscoelastic response of polycarbonate reinforced with short glass fibers

Observations are reported in oscillatory torsion tests at room temperature on unfilled and fiber-reinforced polycarbonates melt-blended with impurities (acronitrile–butadiene–styrene, high-impact polystyrene, low-density polyethylene, poly(ethylene terephthalate) and Nylon 6,6). Constitutive equations are derived for the viscoelastic behavior of glassy polymers. With reference to the theory of cooperative relaxation, a polymer is treated as an ensemble of meso-regions with arbitrary shapes and sizes. The time-dependent response of the ensemble is attributed to rearrangement of meso-domains. The rearrangement events occur at random times, when meso-regions are excited by thermal fluctuations. Stress–strain relations are derived by using the laws of thermodynamics. The governing equations are determined by four adjustable parameters that are found by fitting the experimental data. Fair agreement is demonstrated between the observations and the results of numerical simulation. The study focuses on the effects of the concentration of impurities and glass fibers on material parameters.