How to use reactive stresses to improve plate-theory approximations of the stress field in a linearly elastic plate-like body

All plate theories allow for an approximate evaluation of the stress field in the associated three-dimensional plate-like bodies. We here show how such an approximation can be bettered by an appropriate use of the reactive stress fields that maintain in the three-dimensional body the kinematical constraints implicit in the formulation of a given plate theory. In particular, we discuss transversely extensible, linearly elastic plates and Reissner–Mindlin plates, two theories where the built-in second-order internal constraints turn out to be of importance to improve the evaluation of the stress fields in the corresponding three-dimensional bodies. In addition to arguing our point of view in general, we work out a few explicit examples, both analytically and numerically.     

On bounds and limit theorems for secondary creep in symmetric pressure vessels

The field equations governing creep in spherical and incompressible cylindrical pressure vessels subject to a nondecreasing internal pressure are reduced to a single equation in the effective stress. Using this equation, bounds are obtained for the effective stress and the displacement at any point in the body at any time. Also, in the case where the pressure tends to a limit as t → ∞. limit theorems are obtained which describe the long term behavior of the effective stress and the displacement.     

Non-linear axisymmetric static analysis of orthotropic thin annular plates

This investigation deals with the geometrically non-linear axisymmetric static analysis of elastic orthotropic thin annular plates subjected to uniformly distributed and ring loads. Non-linear differential equations in terms of the transverse displacement w and stress function ψ have been employed. Both w and ψ are expanded in finite power series and the orthogonal point collocation method has been used to obtain discretised algebraic equations from the governing differential equations. Results have been presented for annular plates with and without a rigid plug in the hole; simply supported, clamped and free outer edge conditions; three orthotropic parameters β; two annular ratios and two loading conditions. Results have also been presented illustrating the effect of an elastic rotational constraint at the edge and the effect of a prescribed inplane displacement of the edge.     

Dynamic analysis on generalized linear elastic body subjected to large scale rigid rotations

The dynamic analysis of a generalized linear elastic body undergoing large rigid rotations is investigated. The generalized linear elastic body is described in kine- matics through translational and rotational deformations, and a modified constitutive relation for the rotational deformation is proposed between the couple stress and the curvature tensor. Thus, the balance equations of momentum and moment are used for the motion equations of the body. The floating frame of reference formulation is applied to the elastic body that conducts rotations about a fixed axis. The motion-deformation coupled model is developed in which three types of inertia forces along with their incre- ments are elucidated. The finite element governing equations for the dynamic analysis of the elastic body under large rotations are subsequently formulated with the aid of the constrained variational principle. A penalty parameter is introduced, and the rotational angles at element nodes are treated as independent variables to meet the requirement of C1 continuity. The elastic body is discretized through the isoparametric element with 8 nodes and 48 degrees-of-freedom. As an example with an application of the motion- deformation coupled model, the dynamic analysis on a rotating cantilever with two spatial layouts relative to the rotational axis is numerically implemented. Dynamic frequencies of the rotating cantilever are presented at prescribed constant spin velocities. The maximal rigid rotational velocity is extended for ensuring the applicability of the linear model. A complete set of dynamical response of the rotating cantilever in the case of spin-up maneuver is examined, it is shown that, under the ultimate rigid rotational velocities less than the maximal rigid rotational velocity, the stress strength may exceed the material strength tolerance even though the displacement and rotational angle responses are both convergent. The influence of the cantilever layouts on their responses and the multiple displacement trajectories observed in the floating frame is simultaneously investigated. The motion-deformation coupled model is surely expected to be applicable for a broad range of practical applications.     

A stress analytical solution for Mode III crack within modified gradient elasticity

The strain gradient exists near a crack tip may significantly influence the near-tip stress field. In this paper, the strain gradient and the internal length scales are introduced into the basic equations of mode III crack by the modified gradient elasticity (MGE). By using a complex function approach, the analytical solution of stress fields for mode III crack problem is derived within MGE. When the internal length scales vanish, the stress fields can be simplified to the stress fields of classical linear elastic fracture mechanics. The results show that the singularity of the shear stress is made up of two parts, r^−1/2 part and r^−3/2 part, and the sign of the stress σ_yz changes. With the increase of l_x, the peak value of σ_yz decrease and its location moves farther from the fracture vertex. The influence of strain gradient for mode III crack problem cannot be ignored.     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

GENERALIZATION OF KIRCHHOFF KINETIC ANALOGY TO THIN ELASTIC SHELLS 1)

将弹性细杆的"Kirchhoff动力学比拟"方法推广到弹性薄壳,使弹性薄壳的变形在物理概念上和刚体的运动对应, 在数学表述上等同,从而可以用刚体动力学的理论和方法研究弹性薄壳的变形,为连续的弹性薄壳提供新的离散化方法. 在直法线假设下,在弹性中面上构筑空间正交轴系, 此轴系沿坐标线"运动"的角速度构成两自变量的弯扭度. 沿两个坐标线的弯扭度表达了弹性薄壳的变形和位形,证明了弯扭度之间以及弯扭度与中面切矢间的相容关系. 用Euler角和Lam$\acute{e}$系数表达了非完整约束和中面位形的微分方程,用弯扭度和Lam$\acute{e}$系数表达了应变和应力以及内力及其本构方程.导出了用分布内力集度表达的弹性薄壳在变形后位形上的平衡偏微分方程组,方程的形式与刚体动力学的Euler方程和弹性细杆的Kirchhoff方程具有相似性,实现了Kirchhoff动力学比拟对弹性薄壳的推广.总结了弹性薄壳静力学和刚体动力学以及弹性细杆静力学在概念上的比拟关系.最后给出了一个算例. 为研究弹性薄壳的变形和运动提供新的建模方法和研究思路.也可进一步推广到弹性薄壳动力学.     

Kirchhoff动力学比拟对弹性薄壳的推广

将弹性细杆的"Kirchhoff动力学比拟"方法推广到弹性薄壳,使弹性薄壳的变形在物理概念上和刚体的运动对应, 在数学表述上等同,从而可以用刚体动力学的理论和方法研究弹性薄壳的变形,为连续的弹性薄壳提供新的离散化方法. 在直法线假设下,在弹性中面上构筑空间正交轴系, 此轴系沿坐标线"运动"的角速度构成两自变量的弯扭度. 沿两个坐标线的弯扭度表达了弹性薄壳的变形和位形,证明了弯扭度之间以及弯扭度与中面切矢间的相容关系. 用Euler角和Lam$\acute{e}$系数表达了非完整约束和中面位形的微分方程,用弯扭度和Lam$\acute{e}$系数表达了应变和应力以及内力及其本构方程.导出了用分布内力集度表达的弹性薄壳在变形后位形上的平衡偏微分方程组,方程的形式与刚体动力学的Euler方程和弹性细杆的Kirchhoff方程具有相似性,实现了Kirchhoff动力学比拟对弹性薄壳的推广.总结了弹性薄壳静力学和刚体动力学以及弹性细杆静力学在概念上的比拟关系.最后给出了一个算例. 为研究弹性薄壳的变形和运动提供新的建模方法和研究思路.也可进一步推广到弹性薄壳动力学.      

Deformations of an elastic,internally constrained material. Part 1: Homogeneous deformations

The nonlinear elastic response of a class of materials for which the deformation is subject to an internal material constraint described in experiments by James F. Bell on the finite deformation of a variety of metals is investigated. The purely kinematical consequences of the Bell constraint are discussed, and restrictions on the full range of compatible deformations are presented in geometrical terms. Then various forms of the constitutive equation relating the stress and stretch tensors for an isotropic elastic Bell material are presented. Inequalities on the mechanical response functions are introduced. The importance of these in applications is demonstrated in several examples throughout the paper.This paper focuses on homogeneous deformations. In a simple illustration of the theory, a generalized form of Bell's empirical rule for uniaxial loading is derived, and some peculiarities in the response under all-around compressive loading are discussed. General formulae for universal relations possible in an isotropic elastic, Bell constrained material are presented. A simple method for the determination of the left stretch tensor for essentially plane problems is illustrated in the solution of the problem of pure shear of a materially uniform rectangular block. A general formula which includes the empirical rule found in pure shear experiments by Bell is derived as a special case. The whole apparatus is then applied in the solution of the general problem of a homogeneous simple shear superimposed on a uniform triaxial stretch; and the great variety of results possible in an isotropic, elastic Bell material is illustrated. The problem of the finite torsion and extension of a thin-walled cylindrical tube is investigated. The results are shown to be consistent with Bell's data for which the rigid body rotation is found to be quite small compared with the gross deformation of the tube. Several universal formulas relating various kinds of stress components to the deformation independently of the material response functions are derived, including a universal rule relating the axial force to the torque.Constitutive equations for hyperelastic Bell materials are derived. The empirical work function studied by Bell is introduced; and a new constitutive equation is derived, which we name Bell's law. On the basis of this law, we then derive exactly Bell's parabolic laws for uniaxial loading and for pure shear. Also, form Bell's law, a simple constitutive equation relating Bell's deviatoric stress tensor to his finite deviatoric strain tensor is obtained. We thereby derive Bell's invariant parabolic law relating the deviatoric stress intensity to the corresponding strain intensity; and, finally, Bell's fundamental law for the work function expressed in these terms is recovered. This rule is the foundation for all of Bell's own theoretical study of the isotropic materials cataloged in his finite strain experiments on metals, all consistent with the internal material constraint studied here.     

Controlling parameter of the stress intensity factors for a planar interfacial crack in three-dimensional bimaterials

Numerical solutions of singular integral equations are discussed in the analysis of a planar rectangular interfacial crack in three-dimensional bimaterials subjected to tension. The problem is formulated as a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental density functions are chosen to express singular behavior along the crack front of the interface crack exactly. The calculation shows that the present method gives smooth variations of stress intensity factors along the crack front for various aspect ratios. The present method gives rapidly converging numerical results and highly satisfied boundary conditions throughout the crack boundary. The stress intensity factors are given with varying the material combination and aspect ratio of the crack. It is found that the stress intensity factors K_I and K_II are determined by the bimaterial constant ε alone, independent of elastic modulus ratio and Poisson’s ratio.     

Homogenization of a contact problem for a system of densely situated punches

A linear contact problem of an elastic half space with rigid punches ε-periodically situated on a bounded part of the boundary of the elastic solid is investigated. Using the method of homogenization theory and the method of matched asymptotic expansions, the leading terms of the asymptotic solution are constructed as ε→0. The general capacity of the contact spot is introduced and some its properties are described.     

The torsion of a non-homogeneous stratum with a hyperbolic variation of shear modulus

An exact formulation of the governing dual integral equations for the torsion of a non-homogeneous stratum due to a rigid circular body at its free surface is presented. The stratum varies in shear modulus according to the hyperbolic variation in a contemporary work [1]. It is shown that the unknown static stress distribution under the rigid body is governed by modified Bessel function of the first kind. By comparing the governing functions in the dual integral equations for five cases of elastic media: homogeneous half-space, and stratum, linearly non-homogeneous half-space and stratum and, finally, the present non-homogeneous stratum with hyperbolic variation, it is established that the surface shear modulus is the dominant parameter in the assessment of the stress and displacement fields in a non-homogeneous stratum where lateral variation of elastic properties is negligible.     

A Spectral Theory Formulation for Elastostatics by Means of Tensor Spherical Harmonics

Consider a set of (N+1)-phase concentric spherical ensemble consisting of a core region encased by a sequence of nested spherical layers. Each phase is spherically isotropic and is functionally graded (FG) in the radial direction. Determination of the elastic fields when the outermost spherical surface is subjected to a nonuniform loading and the constituent phases are subjected to some prescribed nonuniform body force and eigenstrain fields is of interest. When the outermost layer is an unbounded medium with zero eigenstrain and body force fields, then an N-phase multi-inhomogeneous inclusion problem is realized. Based on higher-order spherical harmonics, presenting a three-dimensional strain formulation with a robust form of compatibility equations, a spectral theory of elasticity in the spherical coordinate system is developed. Application of the established spectral theory leads to the exact closed-form solution when the elastic moduli of each phase vary as power-law functions of radius.     

Asymptotic fields around an interfacial crack with a cohesive zone ahead of the crack tip

This paper considers an interfacial crack with a cohesive zone ahead of the crack tip in a linearly elastic isotropic bi-material and derives the mixed-mode asymptotic stress and displacement fields around the crack and cohesive zone under plane deformation conditions (plane stress or plane strain). The field solution is obtained using elliptic coordinates and complex functions and can be represented in terms of a complete set of complex eigenfunction terms. The imaginary portion of the eigenvalues is characterized by a bi-material mismatch parameter ε = arctanh(β)/π, where β is a Dundurs parameter, and the resulting fields do not contain stress singularity. The behaviors of “Mode I” type and “Mode II” type fields based on dominant eigenfunction terms are discussed in detail. For completeness, the counterpart for the Mode III solution is included in an appendix.     

A dynamic theory for magnetizable elastic solids with thermal and electrical conduction

A general dynamical theory of magnetizable, electrically and thermally conducting media is developed for soft ferromagnetic or paramagnetic materials in external electromagnetic fields. The general equations are linearized by assuming infinitesimal strains, linear constitutive equations and that all field variables may be divided into two parts: a "rigid body state" and a "perturbation state". The former is the same as the one in rigid body electrodynamics, and the latter which accounts for electromagnetic interaction with the deformable continuum is coupled with stress and strain through linearized field equations. The theory is developed for general anisotropy but specialized for materials with uniaxial, or higher, symmetry.     

Toward a Field Theory for Elastic Bodies Undergoing Disarrangements

Structured deformations are used to refine the basic ingredients of continuum field theories and to derive a system of field equations for elastic bodies undergoing submacroscopically smooth geometrical changes as well as submacroscopically non-smooth geometrical changes (disarrangements). The constitutive assumptions employed in this derivation permit the body to store energy as well as to dissipate energy in smooth dynamical processes. Only one non-classical field G, the deformation without disarrangements, appears in the field equations, and a consistency relation based on a decomposition of the Piola-Kirchhoff stress circumvents the use of additional balance laws or phenomenological evolution laws to restrict G. The field equations are applied to an elastic body whose free energy depends only upon the volume fraction for the structured deformation. Existence is established of two universal phases, a spherical phase and an elongated phase, whose volume fractions are (1?γ₀)³ and (1?γ₀) respectively, with γ₀:=( $ \sqrt 5 $ ?1)/2 the “golden mean”.     

From Non-Linear Elasticity to Linearized Theory: Examples Defying Intuition

Material frame indifference implies that the solution in non-linear elasticity theory for a connected body rigidly rotated at its border is a rigid, stress-free, deformation. If the same problem is considered within linear elasticity theory, considered as an approximation to the true elastic situation, one should expect that if the angle of rotation is small, the body still undergoes a rigid deformation while the corresponding stress, though not zero, remains consistently small. Here, we show that this is true, in general, only for homogeneous bodies. Counterexamples of inhomogeneous bodies are presented for which, whatever small the angle of rotation is, the linear elastic solution is by no means a rigid rotation (in a particular case it is an “explosion”) while the stress may even become infinite. If the same examples are re-interpreted as problems in an elasticity theory based upon genuinely linear constitutive relations which retain their validity also for finite deformations, it is shown that they would deliver constraint reaction forces that are not in equilibrium in the actual, deformed, state. This furnishes another characterization of the impossibility of an exact linear constitutive theory for elastic solids with zero residual stress.      

Modeling lattice strain evolution at finite strains and experimental verification for copper and stainless steel using in situ neutron diffraction

An elasto-plastic self-consistent (EPSC) polycrystal model is extended to account, in an approximate fashion, for the kinematics of large strains, rigid body rotations, texture evolution and grain shape evolution. In situ neutron diffraction measurements of the flow stress, internal strain, texture and diffraction peak intensity evolutions were performed on polycrystalline copper and stainless steel, up to true tensile strains of ε = 0.3. Suitably adjusted slip system hardening model parameters enable the model to quantitatively describe the flow stress of the polycrystalline aggregate. Quantitative predictions of the texture evolution and the internal strain evolution along the stress axis are good, while predictions of transverse internal strains (perpendicular to the tensile loading direction) are less satisfactory. The latter exhibit a large dispersion from grain to grain around a macroscopic average, and the implications of this finding for the interpretation of in situ neutron diffraction method are explored. Finally, as a demonstration of the applicability of the model to problems involving finite rotation, as well as deformation, simulations of simple shear were conducted which predict a texture evolution in agreement with published experimental data, and other modeling approaches as well.     

Bounds on the distribution of extreme values for the stress in composite materials

Suitable macroscopic quantities beyond effective elastic properties are used to assess the distribution of stress within a composite. The composite is composed of N anisotropic linearly elastic materials and the length scale of the microstructure relative to the loading is denoted by ε. The stress distribution function inside the composite λ^ε(t) gives the volume of the set where the norm of the stress exceeds the value t. The analysis focuses on the case when 0<ε?1. A rigorous upper bound on lim_ε→0λ^ε(t) is found. The bound is given in terms of a macroscopic quantity called the macro stress modulation function. It is used to provide a rigorous assessment of the volume of over stressed regions near stress concentrators generated by reentrant corners or by an abrupt change of boundary loading.