Determining material true stress–strain curve from tensile specimens with rectangular cross-section

The uniaxial true stress logarithmic strain curve for a thick section can be determined from the load–diameter reduction record of a round tensile specimen. The correction of the true stress for necking can be performed by using the well-known Bridgman equation. For thin sections, it is more practical to use specimens with rectangular cross-section. However, there is no established method to determine the complete true stress–logarithmic strain relation from a rectangular specimen. In this paper, an extensive three-dimensional numerical study has been carried out on the diffuse necking behaviour of tensile specimens made of isotropic materials with rectangular cross-section, and an approximate relation is established between the area reduction of the minimum cross-section and the measured thickness reduction. It is found that the area reduction can be normalized by the uniaxial strain at maximum load which represents the material hardening and also the section aspect ratio. Furthermore, for the same material, specimens with different aspect ratio give exactly the same true average stress–logarithmic strain curve. This finding implies that Bridgmans correction can still be used for necking correction of the true average stress obtained from rectangular specimens. Based on this finding, a method for determining the true stress–logarithmic strain relation from the load–thickness reduction curve of specimens with rectangular cross-section is proposed.     

Softening in stress–strain curve for Drucker–Prager non-associated plasticity

This paper is devoted to proving some features of the non associated flow rule such as a softening phenomenon in the stress–strain curve and the decrease of limit load. Based on the non-associated Drucker–Prager model, the analysis is investigated by means of a soil specimen subjected to traction and compression actions on its edges. To obtain the stress–strain curve, a semi-analytical approach provides an incremental relation between stresses and strains. The plastic limit load is calculated analytically by direct static and kinematic methods. The kinematic one is determined on the basis of the bipotential concept.     

Analyzing the stress–strain state of thick cylindrical shells

Two approaches to the analysis of the stress–strain state of thick cylindrical shells are elaborated. The shell is divided by concentric cross-sectional circles into several coaxial cylindrical shells. Each of these shells has its own curvature determined on its midline. The stress–strain state of the original shell is described by satisfying the interface conditions between the component shells. The distribution of unknown functions throughout the thickness is determined by finding the analytic solution of a system of differential equations in the first approach and is approximated by polynomial functions in the second approach. The stress–strain state of thick shells is analyzed. It is revealed that the effect of reduction becomes stronger with increasing curvature     

Stress–strain curve and stored energy during uniaxial deformation of polycrystals

The subject of this paper is an attempt to obtain information about the energy stored during plastic deformation from experimentally measured stress–strain curve.Theoretical analysis of the stress–strain curve for elastic-perfectly plastic polycrystalline material has shown that only the part of stored energy can be calculated from the stress–strain curve. This part is the energy stored during non-homogeneous plastic deformation.The results of such calculation have been compared with the total stored energy determined experimentally. It has been shown that part of total stored energy related to non-homogeneous plastic deformation of investigated materials is much lower than that corresponding to homogeneous one.     

On the stress–strain states of cellular materials under high loading rates

A virtual Taylor impact of cellular materials is analyzed with a wave propagation technique, i.e. the Lagrangian analysis method, of which the main advantage is that no pre-assumed constitutive relationship is required. Time histories of particle velocity, local strain, and stress profiles are calculated to present the local stress–strain history curves, from which the dynamic stress–strain states are obtained.The present results reveal that the dynamic-rigid-plastic hardening(D-R-PH) material model introduced in a previous study of our group is in good agreement with the dynamic stress–strain states under high loading rates obtained by the Lagrangian analysis method. It directly reflects the effectiveness and feasibility of the D-R-PH material model for the cellular materials under high loading rates.     

Influence of the cavitation on the stress–strain fields of compressible Blatz–Ko materials at finite deformation

The aim of this article is the analysis of fracture growth in media characterized by random distribution of micro-failure mechanisms per unit volume. The deformation behavior of the material was investigated in terms of a spherical unit cell model, containing an initially spherical cell of porous. The effective elastic bulk modulus as a function of micro-failures concentration was computed and using the Griffith critirium and certain boundary conditions the rate at which the void area varies was determined too. Along the analysis a special form of the strain energy function for compressible Blatz–Ko material was used. The applied traction on the unit cell of the material was determined as a function of the porosity of the material, as well as the strain field within the solid. At low values of the porosity, as the applied external traction was increased instabilities were observed in the void growth.     

Extracting the stress–strain properties of non-linear viscoelastic materials using the bubble inflation test

In this study, an inverse method based on the Levenberg–Marquardt algorithm was evaluated in a numerical experiment to determine the large strain viscoelastic properties from the bubble inflation test. The properties were determined by iteratively matching the calculated bubble pressure–piston displacement data from finite element simulations to a single set of bubble pressure–piston displacement data. The strain-dependent behaviour was characterised by a two-parameter Mooney–Rivlin hyperelastic model, while the time-dependent behaviour was characterised by a three-parameter power law equation. Different initial guesses were used to evaluate the inverse method, and transformation functions were applied to constrain the intermediate guesses to be within bounds. It was found that estimates of the viscoelastic properties could be obtained reasonably using only one set of bubble pressure–piston displacement data. Estimates of the properties were likely affected by the limited time duration of the test, as the behaviour at shorter and particularly larger time scales was less accurately predicted.     

Analysis of the localization of deformation and the complete stress–strain relation for mesoscopic heterogeneous brittle rock under dynamic uniaxial tensile loading

Stress redistribution induced by excavation results in the tensile zone in parts of the surrounding rock mass. It is significant to analyze the localization of deformation and damage, and to study the complete stress–strain relation for mesoscopic heterogeneous rock under dynamic uniaxial tensile loading. On the basis of micromechanics, the complete stress–strain relation including linear elasticity, nonlinear hardening, rapid stress drop and strain softening is obtained. The behaviors of rapid stress drop and strain softening are due to localization of deformation and damage. The constitutive model, which analyze localization of deformation and damage, is distinct from the conventional model. Theoretical predictions have shown to consistent with the experimental results.     

Spline-approximation solution of stress–strain problems for beveled cylindrical shells

The problem of bending of beveled circular cylindrical shells is solved by parametrizing the shell and reducing the two-dimensional boundary-value problem to a one-dimensional one by the spline-collocation method. This problem is solved by the stable discrete-orthogonalization method. The effect of the variability of the geometrical parameters on the displacement fields of circular cylinders is analyzed     

Dynamic problems for and stress–strain state of inhomogeneous shell structures under stationary and nonstationary loads

This paper reviews studies and analyzes results on the effect of discrete ribs on the dynamic characteristics of rectangular plates and cylindrical shells. Use is made of the vibration equations derived from the classical theories of beams, plates, and shells. The effect of Pasternak’s elastic foundation on the critical velocities of a structurally orthotropic model of a ribbed cylindrical shell is determined. Nonstationary problems are solved for perforated and ribbed shells of revolution filled with a fluid or resting on an elastic foundation and subjected to moving or impulsive loads. Results from studies of the behavior of sandwich shell structures under impulsive loads of various types are presented     

Prediction of stress–strain and fracture behaviour of an 8-Harness satin weave ceramic matrix composite

A computationally economic finite-element-based approach has been developed to predict the stress–strain and fracture behaviour of an 8-Harness satin woven ceramic matrix composite with strain-induced damage. The finite element analysis utilises a solid element to model the behaviour of the homogenised orthotropic uni-directional tow and its matrix. The underpinning models of the tow and matrix, (Tang et al., 2009) capture the physics of the interactions between fibres and matrix; and, in this way, permit modelling that bridges the length scales of the fibres and full-scale components. The non-linear multi-axial stress–strain behaviour of the composite has been discretised by multi-linear elastic curves; and the latter has been used as input to a user defined subroutine, UMAT, in the commercial finite element package, ABAQUS. A partial unit cell model has been constructed of the 8-Harness satin weave composite of carbon fibres embedded in an amorphous carbon matrix, HITCO C/C. Predictions of the global stress–strain curve, which include the effects of fibre waviness, have been made for two failure modes: the first by deformation localisation, and the second by dynamic tow failure on fibre fracture, triggered by instantaneous pull-out deactivation. Comparisons have been made between the predictions and experimental data that exhibit two classes of fracture behaviour: brittle and quasi-ductile. The predicted results, both with and without tow waviness, compare well with the experimental data; however, the predictions for waviness are slightly better. The two extremes of experimental behaviour have been found to correspond with the two tow fracture criteria modelled.     

Discretization of the plane problem for a cracked body with nonlinear stress–strain diagram under tension

The plane problem for a cracked body with a piecewise-linear stress–strain diagram under tension is reduced by the Fourier transformation to a system of nonlinear algebraic equations. The system is numerically solved for plane strain and stress states of a perfect elastoplastic material to study plastic zones, stress and strain distributions, and displacements of crack faces     

A method for calculating the stress–strain state in the general boundary-value problem of metal forming—part 1

To simulate metal-forming processes, one has to calculate the stress–strain state of the metal, i.e. to solve the relevant boundary-value problems. Progress in the theory of plasticity in that respect is well known, for example, via the slip-line method, the finite element method, etc.) , yet many unsolved problems remain. It is well known that the slip-line method is scanty. In our opinion the finite element method has an essential drawback. (No one is against the idea of the discretization of the body being deformed and the approximation of the fields of mechanical variables.) The results of calculation of the stress state by the FEM do not satisfy Newtonian mechanics equations (these equations are said to be softened, i.e, satisfied approximately) and stress fields can be considered poor for solution of the subsequent fracture problem. We believe that it is preferable to construct an approximate solution by the FEM and soften the constitutive relations (not Newtonian mechanics equations) , especially as, in any event, they describe the rheology of actual deformable materials only approximately. We seem to have succeeded in finding the solution technique.Here we present some new results for solving rather general boundary-value problems which can be characterized by the following: the anisotropy of the materials handled; the heredity of their properties and compressibility; finite deformations; non-isothermal flow; rapid flow, with inertial forces; a non-stationary state; movable boundaries; alternating and non-classical boundary conditions, etc.Solution by the method proposed can be made in two stages: (1) integration in space with fixed time, with an accuracy in respect of some parameters; (2) integration in time of certain ordinary differential equations for these parameters.In the first stage the method is based on the principle of virtual velocities and stresses. It is proved that a solution does exist and that it is the only possible one. The approximate solution softens (approximately satisfies) the constitutive relations, all the rest of the equations of mechanics being satisfied precisely. The method is illustrated by some test examples.     

Influence of conicity on the stress–strain state of a flexible orthotropic conical shell in a nonstationary magnetic field

A problem of magnetoelasticity for a flexible conical shell in a nonstationary magnetic field is solved. The effect of conicity on the stress–strain state of the shell is analyzed     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

Quasi-soft opto-mechanical behavior of photochromic liquid crystal elastomer: Linearized stress–strain relations and finite element simulations

Based on the neo-classical elastic energy of liquid crystal elastomers, the opto-mechanical behavior is modeled by considering the effect of photoisomerization on the nematic-isotropic transition of liquid crystal phase. Linearized stress–strain relation is derived for infinitesimal deformations with a very unusual shear stress that does not vanish identically as in the case of the soft behavior but is proportional to the rotation of directors. In other words, the shear stress depends on both the shear strain and the skew symmetric part of the displacement gradient with the shear modulus induced by the effect of photoisomerization. Finite element implementation for plane stress problems is obtained through a self-defined material subroutine in ABAQUS FEA tool. Numerical simulations show that the light induced deformations of two dimensional specimens consist of contractions, expansions and bending in different directions. The stress distributions indicate that the driving force for the light induced bending is produced by the bending moment of the normal stress along the director, while the other stress components are much smaller for two dimensional beam shaped specimens. However, the shear stress of the soft LCE is generally nonzero under light illumination due to the inhomogeneity of the opto-mechanical effect. It can be concluded from the strain distributions that the transversal plane cross section could remain plane after deformation if the light intensity or the decay distance is not too small and the sample is in the deep nematic phase. However, the shear strain and in plane rotation are of the same order as the other strain components, and thus should not be neglected. This indicates that the classical simple bending assumptions such as the Euler–Bernoulli beam theory should not be directly applied to model the light induced bending of neo-classical liquid crystal elastomers due to the soft behavior of the materials.     

Determination of the full-field stress and displacement using photoelasticity and sampling moiré method in a 3D-printed model

The quantitative characterization of the full-field stress and displacement is significant for analyzing the failure and instability of engineering materials. Various optical measurement techniques such as photoelasticity, moiré and digital image correlation methods have been developed to achieve this goal. However, these methods are difficult to incorporate to determine the stress and displacement fields simultaneously because the tested models must contain particles and grating for displacement measurement; however, these elements will disturb the light passing through the tested models using photoelasticity. In this study, by combining photoelasticity and the sampling moiré method, we developed a method to determine the stress and displacement fields simultaneously in a three-dimensional (3D)-printed photoelastic model with orthogonal grating. Then, the full-field stress was determined by analyzing 10 photoelastic patterns, and the displacement fields were calculated using the sampling moiré method. The results indicate that the developed method can simultaneously determine the stress and displacement fields.