Strain localization analysis for planar polycrystals based on bifurcation theory

In the present paper, an efficient numerical tool is developed to investigate the ductility limit of polycrystalline aggregates under in-plane biaxial loading. These aggregates are assumed to be representative of very thin sheet metals (with typically few grains through the thickness). Therefore, the plane-stress assumption is naturally adopted to numerically predict the occurrence of strain localization. Furthermore, the initial crystallographic texture is assumed to be planar. Considering the latter assumptions, a two-dimensional single-crystal model is advantageously chosen to describe the mechanical behavior at the microscopic scale. The mechanical behavior of the planar polycrystalline aggregate is derived from that of single crystals by using the full-constraint Taylor scale-transition scheme. To predict the occurrence of localized necking, the developed multiscale model is coupled with bifurcation theory. As will be demonstrated through various numerical results, in the case of biaxial loading under plane-stress conditions, the planar single-crystal model provides the same predictions as those given by the more commonly used three-dimensional single-crystal model. Moreover, the use of the two-dimensional model instead of the three-dimensional one allows dividing the number of active slip systems by two and, hence, significantly reducing the CPU time required for the integration of the constitutive equations at the single-crystal scale. Furthermore, the planar polycrystal model seems to be more suitable to study the ductility of very thin sheet metals, as its use allows us to rigorously ensure the plane-stress state, which is not always the case when the fully three-dimensional polycrystalline model is employed. Consequently, the adoption of this planar formulation, instead of the three-dimensional one, allows us to simplify the computational aspects and, accordingly, to considerably reduce the CPU time required for the numerical predictions.     

Necking of pressurized spherical membranes

The necking of spherical membranes subject to a prescribed increase in enclosed volume is investigated. Attention is restricted to axisymmetric deformations. The materials considered are incompressible, isotropic, time-independent and incrementally linear. A complete set of axisymmetric bifurcation modes is considered and a simple relation is found to govern the critical stress for bifurcation into a given mode. The limiting critical stress and the corresponding mode for short wavelengths are investigated and related to the results obtained from an independent local-necking analysis. Two perturbation methods are employed to study the growth of initial imperfections: one is valid for arbitrary modes, but restricted to small deviations from sphericity, and the other is valid only for the local-necking mode, but is not restricted to small deviations. The effect of path-dependent material behavior on the onset of local necking is explored. Path-dependent material behavior is found to encourage the preferential growth of short wavelength imperfections. Path-independent materials are shown to exhibit significant sensitivity to initial imperfections in the localized-necking mode, although this sensitivity is far less than for a path-dependent material. When account is taken of initial material-property inhomogeneities as well as initial thickness imperfections, it seems that no definite conclusion can be drawn concerning the appropriateness or inappropriateness of an explanation of the onset of localized necking based on a smooth yield-surface plasticity theory and assuming the presence of such initial inhomogeneities.     

板材成型中的分叉和断裂及其空洞扩展效应

用平面应力有限元方法分析空洞模型以模拟一种双相钢板材在成型过程中所遇到的微空洞损伤,经试算可使模型的总体和局部的响应与已有的实验相一致,由此可提供描述该材料的损伤本构参数并研究局部剪切带和扩散型颈缩等分叉现象,临界应变值的分布形成了成型极限图中的下限曲线,当空洞模型的总体应力急剧下降或微裂纹开始出现,其相应的总体应变值提供了上限曲线。     

Elastic instability of a uniformly compressed annular plate with axisymmetric initial deflection

This paper presents a theoretical study of the elastic instability of a uniformly compressed, thin, circular annular plate with axisymmetric initial deflection. The dynamic version of the nonlinear Marguerre plate theory is used, and the linear free vibration problems around the axisymmetric finite deformation of the plate are solved by a finite difference method. By examining the frequency spectrum with various asymmetric modes, the critical compressive load under which the axisymmetric additional deformation of the plate becomes unstable due to the bifurcation buckling is determined, which is found to depend severely on the magnitude of the axisymmetric initial deflection.     

Transition to chaotic vibrations for harmonically forced perfect and imperfect circular plates

The transition from periodic to chaotic vibrations in free-edge, perfect and imperfect circular plates, is numerically studied. A pointwise harmonic forcing with constant frequency and increasing amplitude is applied to observe the bifurcation scenario. The von Kármán equations for thin plates, including geometric non-linearity, are used to model the large-amplitude vibrations. A Galerkin approach based on the eigenmodes of the perfect plate allows discretizing the model. The resulting ordinary-differential equations are numerically integrated. Bifurcation diagrams of Poincaré maps, Lyapunov exponents and Fourier spectra analysis reveal the transitions and the energy exchange between modes. The transition to chaotic vibration is studied in the frequency range of the first eigenfrequencies. The complete bifurcation diagram and the critical forces needed to attain the chaotic regime are especially addressed. For perfect plates, it is found that a direct transition from periodic to chaotic vibrations is at hand. For imperfect plates displaying specific internal resonance relationships, the energy is first exchanged between resonant modes before the chaotic regime. Finally, the nature of the chaotic regime, where a high-dimensional chaos is numerically found, is questioned within the framework of wave turbulence. These numerical findings confirm a number of experimental observations made on shells, where the generic route to chaos displays a quasiperiodic regime before the chaotic state, where the modes, sharing internal resonance relationship with the excitation frequency, appear in the response.     

Non-normality and bifurcation in plane strain tension and compression

The bifurcations of a rectangular block subject to plane strain tension or compression are investigated. The block material is taken to be incompressible and is characterized by an incrementally linear constitutive law for which “normality” does not necessarily hold. The consequences of non-normality regarding bifurcation are given primary emphasis here. The characteristic regimes of the governing equations (elliptic, parabolic and hyperbolic) are detennined. In each of these regimes both symmetric and antisymmetric diffuse bifurcation modes are available. Additionally, in the hyperbolic and parabolic regimes, bifurcation into a localized shear band mode is also possible. Particular attention is given to the limiting cases of long wavelength and soon wavelength diffuse bifurcation modes. The range of parameter values is identified for which bifurcation into some localized mode may precede bifurcation into a long wavelength diffuse mode. Some difficulties associated with employing a linear incremental solid in a bifurcation analysis, when primary interest is in the bifurcation of an underlying elastic-plastic solid, are also discussed.     

Axisymmetric bifurcation analysis in soils by the tangential-subloading surface model

This article discusses localized bifurcation modes corresponding to shear band formation and diffuse bifurcation modes corresponding to bulge formation for cylindrical soil specimen subjected to an axisymmetric load under undrained conditions. We employ the tangential-subloading surface model, which exhibits the characteristic regimes of the governing equations: elliptic, hyperbolic and parabolic. Also, conditions for shear band formation, shear band inclination, diffuse bulging formation, and the long and short wavelength limits of diffuse bulging modes are discussed in relation to material properties and their state of stress, i.e. the stress ratio and the normal-yield ratio. Tangential-plastic strain rate term is required for the analyses of shear band and diffuse bulging. The shear band and the diffuse bulging are generated in not only normal-yield but also subyield states and they are severely affected by the normal-yield ratio describing the degree of approach to the normal-yield state.     

A justification of nonlinear properly invariant plate theories

A single asymptotic derivation of three classical nonlinear plate theories is presented in a setting which preserves the frame-invariance properties of three-dimensional finite elasticity. By a successive scaling of the external loading on the three-dimensional body, the nonlinear membrane theory, the nonlinear inextensional theory and the von Kármán equations are derived as the leading-order terms in the asymptotic expansion of finite elasticity. The governing equations of the nonlinear inextensional theory are of particular interest where 1) plane-strain kinematics and plane-stress constitutive equations are derived simultaneously from the asymptotic analysis, 2) the theory can be phrased as a minimization problem over the space of isometric deformations of a surface, and 3) the local equilibrium equations are identical to those arising in the one-director Cosserat shell model. Furthermore, it can be concluded that with a regular, single-scale asymptotic expansion it is not possible to obtain a system of plate equations in which finite membrane strain and finite bending strain occur simultaneously in the leading-order term of an asymptotic analysis.     

Effect of thermal gradient on vibration of non-homogeneous visco-elastic elliptic plate of variable thickness

An analysis and numerical results are presented for free transverse vibrations of non-homogeneous visco-elastic elliptic plate whose temperature and thickness spatial variations both are parabolic along a line through plate centre. The variation in density is assumed as parabolic along a line through plate centre, which raise non-homogeneity of the plate materials and make problem interesting as introducing variation in non-homogeneity of the material mass density reduce the problem practical importance in comparison to homogenous plate. For visco-elastic, the basic elastic and viscous elements are combined. We have taken Kelvin model for visco-elasticity that is the combination of the elastic and viscous elements in parallel. Here the elastic element means the spring and the viscous element means the dashpot. The governing differential equation of motion has been solved by Galerkin’s technique. Deflection, time period and logarithmic decrement corresponding to the first two modes of vibrations of a clamped non-homogeneous visco-elastic elliptic plate for various values of taper constant, thermal constants, non-homogeneity constant and aspect ratio are obtained and shown graphically.     

Shear-band bifurcation of an isotropic compressible hyperelastic solid

This paper concerns shear-band bifurcations from the homogeneous finite plane deformation of an isotropic compressible hyperelastic solid. The governing equations for the incremental plane deformation superposed on the initial finite deformation are derived and then the equilibrium equations in terms of incremental displacements are classified into the elliptic type, parabolic type, etc. From this classification follows a restriction which should be placed on the strain-energy function in order that the equilibrium equations may be either elliptic or parabolic for all principal stretches. For the hyperelastic solid complying with this restriction, the condition for the shear-band bifurcation is obtained. Finally, the incremental displacement field of an infinite series of shear bands in a slab sandwiched between slippery rigid layers is established.     

Three-dimensional numerical simulation of the wake flow of an afterbody at subsonic speeds

We numerically investigate the wake flow of an afterbody at low Reynolds number in the incompressible and compressible regimes. We found that, with increasing Reynolds number, the initially stable and axisymmetric base flow undergoes a first stationary bifurcation which breaks the axisymmetry and develops two parallel steady counter-rotating vortices. The critical Reynolds number (Re _cs) for the loss of the flow axisymmetry reported here is in excellent agreement with previous axisymmetric BiGlobal linear stability (BiGLS) results. As the Reynolds number increases above a second threshold, Re _co, we report a second instability defined as a three-dimensional peristaltic oscillation which modulates the vortices, similar to the sphere wake, sharing many points in common with long-wavelength symmetric Crow instability. Both the critical Reynolds number for the onset of oscillation, Re _co, and the Strouhal number of the time-periodic limit cycle, St_sat, are substantially shifted with respect to previous axisymmetric BiGLS predictions neglecting the first bifurcation. For slightly larger Reynolds numbers, the wake oscillations are stronger and vortices are shed close to the afterbody base. In the compressible regime, no fundamental changes are observed in the bifurcation process. It is shown that the steady state planar-symmetric solution is almost equal to the incompressible case and that the break of planar symmetry in the vortex shedding regime is retarded due to compressibility effects. Finally, we report the developments of a low frequency which depends on the afterbody aspect ratio, as well as on the Reynolds and on the Mach number, prior to the loss of the planar symmetry of the wake.     

Elastic instability of an orthotropic elliptic plate

On the basis of von Kármán equations and using the general bifurcation theory,theelastic instability of an orthotropic elliptic plate whose edge is subjected to a uniform planecompression is discussed.Following the well-known Liapunov-Schmidt process theexistance of bifurcation solution at a simple eigenvalue is shown and the asymptoticexpression is obtained by means of the perturbation expansion with a small parameter.Finally,by using the finite element method,the critical loads of the plate are computed andthe post-buckling behavior is analysed.And also the effect of material and geometric parameters on the stability is studied.     

Bifurcation in compressible elastic/plastic cylinders under uniaxial tension

The bifurcation problem of a circular cylinder of elastic/plastic material under uniaxial tension is investigated, with particular reference to the usual engineering criterion that necking is initiated when the load on the specimen reaches a maximum. The material considered is compressible, with a smooth yield surface and associated flow rule. A lower bound analysis shows that for the particular constitutive equation chosen bifurcation cannot occur under a range of loading conditions while the stress is less than a certain value which is itself slightly less than the stress at the maximum load point. Diffuse axisymmetric necking modes under the commonly assumed loading conditions of prescribed axial components of velocity and shear-free traction-rates on the ends are, however, found to be initiated always after maximum load, the delay depending on the same factors shown for an incompressible material in reference [1]. The effect of the elastic compressibility assumption is to reduce the delay for a wide range of geometries, but to increase it for very slender specimens, as compared with the incompressible case. Surface modes are also found, but at stresses of an unrealistically high order of magnitude.     

Bifurcation phenomena in the plane compression test

A basic, compression, bifurcation problem is studied by methods similar to those used by R. Hill and J. W. Hutchinson (1975) for the corresponding tension problem. Bifurcations from a state of homogeneous in-plane compression loading are investigated for a rectangular block of incompressible material constrained to undergo plane deformations. The sides of the block are tractionfree, and it is loaded compressively by a uniform, shear-free, relative displacement of its ends. For a wide class of incrementally-linear time-independent materials only two instantaneous moduli enter into the analysis. Diffuse modes of both symmetric and antisymmetric bifurcation are examined in the elliptic regime of the governing equations, and the possibility of localized modes is considered both inside and outside this regime. Lowest bifurcation stresses are computed for essentially the entire range of possible combinations of material properties and geometry, and these are compared with results obtained by Hill and Hutchinson for the tension problem. The limiting case of large thickness (the semi-infinite block) is considered, confirming the results of M. A. Biot (1965).     

Anisotropic tensile necking

Tensile necking in anisotropic bars is analyzed in the spirit of P. W. Bridgman's treatment of the isotropic case. Anisotropic plastic flow causes an initially axisymmetric bar to develop an elliptical neck. Using physical approximations analogous to Bridgman's, an approximate analytical solution for the stress distribution is obtained. The anisotropic necking theory is then applied to the analysis of tension tests conducted on Zircaloy-4 at elevated temperatures. The utility of the theory is demonstrated and the adequacy of inherent approximations assessed.     

Second approximation in the problem of the strong viscous interaction on slender three-dimensional bodies

We consider the hypersonic flow of a perfect gas past a slender three-dimensional body in a regime of strong viscous interaction. We give equations which make it possible to reduce the problem of determining the aerodynamic characteristics of a body which is not axisymmetric to the problem of computing the flow past an equivalent body of rotation at zero angle of incidence. The second approximation for the heat transfer and drag coefficients is found by the method of external and internal combinations of asymptotic expansions. The region in which this method can be applied and the accuracy of the asymptotic theory are estimated by comparison with exact numerical computations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti, i Gaza, No. 5, pp. 107–113, September–October, 1970.     

Post-bifurcation of perfect and imperfect spherical elastic membranes

When a spherical elastic membrane is inflated it is well known that it may bifurcate into an aspherical mode after the pressure maximum is reached. Upon further inflation the spherical configuration is regained. Here we follow the developing aspherical solution path, for specific forms of strain-energy function, using a simple numerical method. For a realistic strain-energy function it is shown that the post-bifurcation solution curve connects the two bifurcation points. We also consider the inflation of imperfect spherical membranes and show that bifurcation may still occur. For the class of Ogden materials we investigate the asymptotic shape of arbitrary axisymmetric membranes.     

Asymmetric bifurcations in a pressurised circular thin plate under initial tension

It has been known for some time that under certain circumstances the axisymmetric solution describing the deformation experienced by a stretched circular thin plate or membrane under sufficiently strong normal pressure does not represent an energy-minimum configuration. By using the method of adjacent equilibrium a set of coordinate-free bifurcation equations is derived here by adopting the Föppl–von Kármán plate theory. A particular class of asymmetric bifurcation solutions is then investigated by reduction to a system of ordinary differential equations with variable coefficients. The localised character of the eigenmodes is confirmed numerically and we also look briefly at the role played by the background tension on this phenomenon.     

Material aspects of dynamic neck retardation

Neck retardation in stretching of ductile materials is promoted by strain hardening, strain-rate hardening and inertia. Retardation is usually beneficial because necking is often the precursor to ductile failure. The interaction of material behavior and inertia in necking retardation is complicated, in part, because necking is highly nonlinear but also because the mathematical character of the response changes in a fundamental way from rate-independent necking to rate-dependent necking, whether due to material constitutive behavior or to inertia. For rate-dependent behavior, neck development requires the introduction of an imperfection, and the rate of neck growth in the early stages is closely tied to the imperfection amplitude. When inertia is important, multiple necks form. In contrast, for rate-independent materials deformed quasi-statically, single necks are preferred and they can emerge in an imperfection-free specimen as a bifurcation at a critical strain. In this paper, the interaction of material properties and inertia in determining neck retardation is unraveled using a variety of analysis methods for thin sheets and plates undergoing plane strain extension. Dimensionless parameters are identified, as are the regimes in which they play an important role.     

Buckling of elastic-plastic shells of revolution with discrete elastic-plastic ring stiffeners

The theory is summarized for axisymmetric prebuckling and nonsymmetric bifurcation buckling of ring-stiffened shells of revolution. The analysis is based on finite difference energy minimization in which moderately large meridional rotations, elastic-plastic effects, and primary or secondary creep are included. This theory is implemented in a computer program called BOSOR5, for the analysis of segmented and branched ring-stiffened shells of revolution of multi-material construction.Comparisons between test and theory are given for axisymmetric collapse and nonsymmetric bifurcation buckling of 69 machined ring-stiffened aluminum cylinders submitted to external hydrostatic pressure. Because most of the cylinders fail at an average stress which corresponds to the knee of the stress-strain curve, the analytical predictions are not very sensitive to modeling particulars such as nodal point density or boundary conditions. Agreement between test and theory is improved if the analytical model reflects the fact that the shell and rings intersect over finite axial lenths.