Oscillatory squeezing flow of a biological material

Large-amplitude oscillatory squeezing flow data are reported for a complex biological material, which is highly shear-thinning in oscillatory shear flow. This soft tissue has a linear viscoelastic limit at a strain of approximately 0.2%. The oscillatory squeezing flow data at large strain are analyzed using two constitutive models: a bi-viscosity Newtonian model, and a non-linear Maxwell model. It is found that although both models may have the same response in shape, the later matches with our non-linear experimental data better. It is also concluded that the non-linear response of the material in large amplitude oscillatory flow is mainly due to the shear thinning of the material. Received: 9 February 2000/Accepted: 22 February 2000     

Instability and failure of internally pressurized ductile metal cylinders

Forlong, ductile, thick-walled tubes under internal pressure instabilities and final failure modes are studied experimentally and theoretically. The test specimens are closed-end cylinders made of an aluminum alloy and of pure copper and the experiments have been carried out for a number of different initial external radius to internal radius ratios. The experiments show necking on one side of the tubes at a stage somewhat beyond the maximum internal pressure. All tubes, except for one aluminum alloy tube, failed by shear fracture under decreasing pressure. The aluminum alloy tubes exhibited localized shear deformations in the neck region prior to fracture and also occasionally surface wave instabilities. The numerical investigation is based on an elastic-plastic material model for a solid that develops a vertex on the yield surface, using representations of the uniaxial stress-strain curves found experimentally. In contrast to the simplest flow theory of plasticity this material model predicts shear band instabilities at a realistic level of strain. A rather sharp vertex is used in the material model for the aluminum alloy, while a more blunt vertex is used to characterize copper. The theoretically predicted bifurcation into a necking mode, the cross-sectional shape of the neck, and finally the initiation and growth of shear bands from the highly strained internal surface in the neck region are in good agreement with the experimental observations.     

Plastic spin associated with a non-normality theory of plasticity

A plasticity model using a vertex-type plastic flow rule on a smooth yield surface for an anisotropic solid has been proposed recently. This model is here completed by incorporating the effect of plastic spin. Simple shear with a large shear strain is one of the hardest tests on finite strain anisotropic plasticity models, and it is here shown which plastic spin expression is needed to avoid unrealistic oscillatory behavior of the shear stress under large shear strains. The idea of using non-normality with a smooth yield surface originates from a recent proposal of using an abrupt strain path change to determine the subsequent yield surface shape. For this method both polycrystal plasticity calculations and experiments have shown a vertex-type response on the apparently smooth yield surface.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

On the application of the additive decomposition of generalized strain measures in large strain plasticity

In the present paper two thermodynamically consistent large strain plasticity models are examined and compared in finite simple shear. The first model (A) is based on the multiplicative decomposition of the deformation gradient, while the second one (B) on the additive decomposition of generalized strain measures. Both models are applied to a rigid-plastic material described by the von Mises-type yield criterion. Since both models include neither hardening nor softening law, a constant shear stress response even for large amounts of shear is expected. Indeed, the model A exhibits the true constant shear stress behavior independent of the elastic material law. In contrast, the model B leads to a spurious shear stress increase or drop such that its applicability under finite shear deformations may be questioned.     

On shear band formation: I. Constitutive relationship for a dual yield model

A phenomenological constitutive relation, for capturing the shear band formation in a rate-independent elastic-plastic material, is established. The model takes into account both the J₂-isotropic flow and a threshold shear stress-based flow. The elastic-plastic constitutive tensor is expressed explicity in terms of elastic constants, the deviatoric stress tensor, the direction of the principal shear velocity-strain, and other material constants. This model particularly facilitates the resolution of the formation of the shear band even under material hardening conditions and does not demand an a priori knowledge of the orientation of the shear band. This is incorporated in an FEM, and the plane strain tensile test of Anand and Spitzig [1980] is numerically simulated. The computed results compare favorably with the experimental data. The shear band emerges more naturally as a solution to the boundary value problem, unlike the situations in solutions based on classical bifurcation methods. Nevertheless, the usefulness of the local instability condition (Ortiz et al. [1987]) is also demonstrated.     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

Modeling the microstructural evolution of metallic polycrystalline materials under localization conditions

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic-plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material's local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length was 870-890 μm with deformation imposed through a Split-Hopkinson Pressure Bar system at maximum top-to-bottom velocity in the range of 8-25 m/s. The deformation behavior of tantalum tophat samples is modeled through direct polycrystal simulations. An embedded Voronoï-tessellated two-dimensional microstructure is used to represent the material within the shear zone of the sample. A thermo-mechanically coupled elasto-viscoplastic single crystal model is presented and used to represent the response of the grains within the aggregate shear zone. In the shoulder regions away from the shear zone where strain levels remain on the order of 0.05, the material is represented by an isotropic J₂ flow theory based upon the elasto-viscoplastic Mechanical Threshold Stress (MTS) model for flow strength. The top surface stress versus displacement results were compared to those of the experiments and over-all the simulated stress magnitude is over-predicted. It is believed that the reason for this is that the simulations are two-dimensional. A region within the numerical shear zone was isolated for statistical examination. The vonMises stress state within this isolated shear zone region suggests an approximate normal distribution with a factor of two difference between the minimum and maximum points in the distribution. The equivalent plastic strain distribution within this same region has values ranging between 0.4 and 1.5 and is not symmetric. Other material state distributions are also given. The crystallographic texture within this isolated shear zone is also compared to the experimental texture and found to match reasonably well considering the simulations are two-dimensional.     

A finite strain theory for elastic-plastic deformation

An elastic-plastic theory that is applicable when the elastic part of the strain is finite is proposed. A flow rule for an incompressible solid is obtained from Drucker's postulate [1]. Isothermal simple shear of a material which is neo-Hookean both before yielding and during elastic unloading after yielding is considered as an application of the theory. The problem is solved for two yield conditions and associated flow rules.     

复杂应力状态下木材力学性能的数值模拟

针对木材复杂的各向异性材料特点,建立了能反映木材正交各向异性弹性、抗拉和抗压强度不等、抗拉或抗剪时发生脆性破坏而受压时发生塑性变形等特性的本构模型。将木材弹性应力—应变关系简化为正交各向异性;选用Yamada-Sun强度准则来判断木材抗压时是否屈服,抗拉或抗剪时是否发生应变软化;通过引入损伤因子和弹性应变能,建立了木材...     

Anisotropic hardening and non-associated flow in proportional loading of sheet metals

Conventional isotropic hardening models constrain the shape of the yield function to remain fixed throughout plastic deformation. However, experiments show that hardening is only approximately isotropic under conditions of proportional loading, giving rise to systematic errors in calculation of stresses based on models that impose the constraint. Five different material data for aluminum and stainless steel alloys are used to calibrate and evaluate five material models, ranging in complexity from a von Mises’ model based on isotropic hardening to a non- associated flow rule (AFR) model based on anisotropic hardening. A new model is described in which four stress–strain functions are explicitly integrated into the yield criterion in closed form definition of the yield condition. The model is based on a non-AFR so that this integration does not affect the accuracy of the plastic strain components defined by the gradient of a separate plastic potential function. The model not only enables the elimination of systematic errors for loading along the four loading conditions, but also leads to a significant reduction of systematic errors in other loading conditions to no higher than 1.5% of the magnitude of the predicted stresses, far less that errors obtained under isotropic hardening, and at a level comparable to experimental uncertainty in the stress measurement. The model is expected to lead to a significant improvement in stress prediction under conditions dominated by proportional loading, and this is expected to directly improve the accuracy of springback, tearing, and earing predictions for these processes. In addition, it is shown that there is no consequence on MK necking localization due to the saturation of the yield surface in pure shear that occurs with the aluminum alloys using the present model.     

Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.     

Measurement of elastic-plastic stresses by scattered-light photomechanics

Several birefringent materials were studied for their suitability for use in three-dimensional photoplasticity. This study resulted in the selection of cellulose propionate as model material. Its close match in index of refraction with ordinary mineral oil makes cellulose propionate suitable for scattered-light photomechanics. Viscoelastic behavior of the material is used to simulate elastic-plastic behavior of metals. The stress, strain and optical behavior of the material has been studied under slow loading in finite steps. A successful solution of an elastic-plastic torsion problem was obtained, demonstrating the applicability of the techniques of scattered-light photoplasticity to three-dimensional problems. This experimental method does not require the unloading and slicing of the model, avoiding any errors that might be introduced by residual stresses due to unloading. A major advantage of the method is the use of live loading, which allows the investigation of several load levels with the same model. Stress- and strain-concentration factors for the grooved shaft in torsion showed excellent agreement with Neuber's analytical results. Distributions of shear stress and shear strain across the minimum section of the model were compared to elastic theory. Integration of the shear-stress distributions showed good agreement with the measured values of applied torque.     

Calibrating a J2 plasticity material model using a 2D inverse finite element procedure

Material models are the key ingredients to accurately capture the global mechanical response of structural systems. The use of finite element analysis has proven to be effective in simulating nonlinear engineering applications. However, the choice of the appropriate material model plays a big role in the value of the numerical predictions. Such models are not expected to exactly reproduce global experimental response in all cases. Alternatively, the measured global response at specific domain or surface points can be used to guide the nonlinear analysis to successively extract a representative material model. By selecting an initial set of stress–strain data points, the load–displacement response at the monitoring points is computed in a forward incremental analysis without iterations. This analysis retains the stresses at the integration points. The corresponding strains are not accurate since the computed displacements are not anticipated to match the measured displacements at the monitoring points. Therefore, a corrective incremental displacement analysis is performed at the same load steps to adjust for displacements and strains everywhere by matching the measured displacements at the monitoring points. The stress–strain vectors at the most highly stressed integration point are found to establish an improved material model. This model is used within a multi-pass incremental nonlinear finite element analysis until the discrepancy between the measured and the predicted structural response at the monitoring points vanishes. The J2 flow theory of plasticity is used as a constitutive framework to build the tangent elastic–plastic matrices. The applicability of the proposed approach is demonstrated by solving 2D inverse continuum problems. The comparisons presented support the effectiveness of the proposed approach in accurately calibrating the J2 plasticity material model for such problems.     

Pressure-shear waves in 6061-T6 aluminum and alpha-titanium

Thepressure-shear plate impact technique is used to study material behavior at high rates of deformation. In this technique, plastic waves of combined pressure and shear stresses are produced by impact of parallel plates skewed relative to their direction of approach. Commercially pure alpha-titanium and 6061-T6 aluminum are tested under a variety of pressure and shear tractions by using different combinations of impact velocities and angles of inclination. A laser interferometer system is used to monitor simultaneously the normal and transverse components of motion of a point at the rear surface of the target plate. The experimental results are compared with numerical solutions based on an elastic/viscoplastic model of the material. Both isotropic and kinematic strain hardening models are used in the computations. The results indicate that unlike the normal velocity profiles, the transverse velocity profiles are sensitive to the dynamic plastic response and, thus, can be used to study material behavior at high strain rates. For the materials tested the results suggest that the flow stress required for plastic straining increases markedly with increasing strain rate at strain rates above 10⁴s^?1. Hydrostatic pressure of the order that exists in the tests (up to 2 GPa) does not affect the plastic flow in 6061-T6 aluminum and appears to have at most a minor effect on the deformation of the titanium.     

An analysis of ductile rupture in notched bars

Ductile fracture in axisymmetric and plane strain notched tensile specimens is analyzed numerically, based on a set of elastic-plastic constitutive relations that account for the nucleation and growth of microvoids. Final material failure by void coalescence is incorporated into the constitutive model via the dependence of the yield function on the void volume fraction. In the analyses the material has no voids initially; but as the voids nucleate and grow, the resultant dilatancy and pressure sensitivity of the macroscopic plastic flow influence the solution significantly. Considering both a blunt notch geometry and a sharp notch geometry in the computations permits a study of the relative roles of high strain and high triaxiality on failure. Comparison is made with published experimental results for notched tensile specimens of high-strength steels. All axisymmetric specimens analyzed fail at the center of the notched section, whereas failure initiation at the surface is found in plane strain specimens with sharp notches, in agreement with the experiments. The results for different specimens are used to investigate the circumstances under which fracture initiation can be represented by a single failure locus in a plot of stress triaxiality vs effective plastic strain.     

Flow localization in the plane strain tensile test

The deformations in a plane strain tensile test are analyzed numerically, both for a solid characterized by a phenomenological corner theory of plasticity and for a nonlinear elastic solid. As opposed to the simplest flow theory of plasticity with a smooth yield surface, both these material models exhibit shear band instabilities at a realistic level of strain. Initial imperfections are specified in the form of thickness inhomogeneities. A long-wavelength imperfection grows into the well-known necking mode and subsequently, at a sufficiently high local strain level, bands of intense shear deformations develop in the necking region. The location of these shear bands is strongly influenced by the location of small strain concentrations near the surface, induced by various short-wave patterns of initial thickness imperfections. In accord with the non-uniform straining in the neck it is found that the intensity of the localized deformations varies along the bands, and some of the shear bands end inside the material.     

一种考虑非比例附加损伤的多轴低周疲劳模型

在实际工作环境中,机械结构往往承受着多轴非比例循环载荷.相比多轴比例循环加载,多轴非比例循环加载由于产生了附加强化现象,造成机械结构疲劳寿命下降.通过分析薄壁圆筒管件在非比例加载工况下应力应变变化规律和发生破坏位置,本文基于临界面法提出一种考虑多轴非比例附加损伤的疲劳模型.该模型将最大剪切应变幅平面作为临界面,提出一个新的附加强化因子,结合临界面上切应变幅和正应变幅组成新的多轴疲劳损伤参量.此参量不仅考虑了非比例加载下临界面上正应变幅和切应变幅对材料造成的疲劳损伤,还考虑到应变路径的变化和材料非比例加载敏感特性对材料疲劳寿命的影响.考虑到实际情况下模型所需材料附加强化系数有时难以获得的情况,给出了材料附加强化系数的有关近似计算公式.只需要材料基本力学参数便可得到材料附加强化系数,方便工程实际应用.采用8种材料的多轴疲劳寿命数据对提出的新模型进行检验,结果表明所提出的新模型与传统多轴疲劳模型相比预测寿命精度更高.      

Research note on a reexamination of neutral loading experiments

In the examination of the published results from neutral loading experiments, the question as to whether plastic deformation occurs is found to depend on both the material and initial loading strain. Provided that initial loading is elastic, then a subsequent stress path that follows the boundary of the initial yield surface for a hardening material is truly neutral with a wholly elastic response. However, when initial loading is elastic-plastic, then further plastic deformation is produced from a subsequent stress path that follows an isotropic expansion of the initial yield surface. These results enable the appropriateness of the kinematic hardening rule and more recent developments in plasticity theory to be appraised. Neutral loading of a non-hardening material produces plastic flow. Whether the absence of hardening is inherent or induced by plastic prestrain, it is shown that the Prandtl-Reuss theory then represents the observed behavior. In general, the purely elastic and nonhardening solutions provide respectively lower and upper bounds on the deformation.