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.     

The effect of porosity distribution on ductile failure

The effect of a nonuniform distribution of porosity on flow localization and failure in a porous material is analyzed numerically. The void density distribution and properties used to characterize the material behavior were obtained from measurements on partially consolidated and sintered iron powder. The calculations were carried out using an elastic viscoplastic constitutive relation for porous plastic solids. Local material failure is incorporated into the model through the dependence of the flow potential on void volume fraction. The region modelled is a small portion of a larger body, subject to various triaxial stress conditions. Both plane strain and axisymmetric deformations are considered with imposed periodic boundary conditions. Interactions between regions with higher void fractions promote plastic flow localization into a band. Local failure occurs by void growth and coalescence within the band. The results suggest a failure criterion based on a critical void volume fraction that is only weakly dependent on stress history. The critical void fraction does. however, depend on the initial void distribution and material hardening characteristics.     

Modelling of rapid pressure-strain in Reynolds stress closures — Difficulties associated with rotational mean flows

Intercomponent energy transfer within the context of Reynolds stress closures is studied. Attention is focussed on the rapid limit of homogeneous flow situations where this energy transfer is caused solely by the rapid pressure strain rate. We present and analyze the performance of the recently proposed rapid pressure strain rate model of Johansson & Hallbäck (J Fluid Mech. 1994) in various homogeneous (rapid) flow situations, and compare with results obtained with other models from the literature and rapid distortion solutions. The prediction difficulties associated with rotational mean flows are analyzed. A generally formulated test case, which as special cases comprises, e.g. plane strain and homogeneous shear flow, is used to illustrate the modelling difficulties associated with rotational mean flows. An axisymmetric case is used to demonstrate that parts of the spectrum with anti-reflectional symmetry, which are instrumental for the dynamics when rotational effects are present, are totally missed in classical Reynolds stress closures. A closer prediction in cases with strong influence of rotation would require introduction of other transported quantities.     

Effect of yield surface curvature and void nucleation on plastic flow localization

The effect of void nucleation is incorporated in a recently proposed material model that accounts for a combination of kinematic hardening and isotropic hardening of a porous ductile material. Since each of plastic dilatancy, void nucleation and yield surface curvature have a strong influence on predictions of plastic flow localization, the present material model can be used to study the interaction of these effects. Nucleation controlled by the plastic strain as well as nucleation controlled by the maximum normal stress on the particle-matrix interface are modelled. The predictions of the material model, for various combinations of parameters, are illustrated by analyses of shear band formation under plane strain or axisymmetric conditions, and by analyses of necking in biaxially stretched sheets.     

DEM SIMULATION OF LIQUEFACTION FOR GRANULAR MEDIA UNDER UNDRAINED AXISYMMETRIC COMPRESSION AND PLANE STRAIN CONDITIONS

Based on three dimensional (3D) Discrete Element Method (DEM), the paper presents simulation results of undrained tests on loose assemblies of polydisperse spheres under axisymmetric compression and plane strain conditions using a periodic cell. In the present work, undrained tests were modelled by deforming the samples under constant volume conditions. The undrained (effective) stress paths are shown to be qualitatively similar to experimental results in literature. A microscopic parameter in terms of redundancy factor (RF) is used to identify the onset of liquefaction (or temporary liquefaction), with the condition of RF equal to unity defining the transition from ’solid-like’ to ’liquid-like’ behaviour. It is found that the undrained behaviour is governed by the evolution of redundancy factor under both undrained axisymmetric compression and plane strain conditions, and a reversal of deviatoric stress in stress path for medium loose systems occurs due to the fact that the system becomes a structural mechanism (RF<1) transiently at the microscopic level during the evolution.     

Stress Analysis of an Orthotropic Work-Hardening Cylinder with Body Force*

ABSTRACT

An elastoplastic analysis of an axisymmetric cylinder subjected to linear body forces is presented. The effect of reinforcement and anisotropy are also included. Classical plasticity and familiar assumptions of plane stress and strain are used to arrive at closed-form solutions for the case of linear body forces. The problem is solved for the general case in which orthotropy is considered in the elastic range. For the case of plasticity, first the isotropic yield functions (von Mises and Tresca) are used and then the problem is extended to the case of Hill's yield criterion. Closed-form solutions are found for both the von Mises (plane strain) and Tresca (plane stress and strain) cases.     

运用Hoek-Brown经验准则分析圆形硐室围岩弹塑性应力和位移

Hoek-Brown经验准则能够较容易地用来估计节理岩体的强度,在岩石工程界已得到广泛的应用和认可。只要用定量指标合理地描述岩体质量,就可确定岩体的强度,因此Hoek-Brown经验准则为通过岩体质量评分评价岩体稳定性架起了一座桥梁,并对需要采取的加固支护措施提出建议。地下硐室开挖后改变了岩体的初始应力状态,围岩应力产生应力重分布现象。当硐室周边围岩应力状态超过岩体弹性极限状态而进入塑性状态时,塑性区内岩体的应力满足极限平衡条件。本文简单回顾了当侧压力系数为1时,圆形硐室围岩的弹性应力和位移,在此基础上以Hoek-Brown经验准则为极限平衡条件,得到圆形硐室轴对称平面应变问题的围岩弹塑性应力和位移的分析解。     

纤维与基体的轴对称界面端附近的奇异应力场

基于各向性弹性力学空间轴对称问题的基本方程,研究了纤维与基体的轴对称界面端的应力奇异性,并给出了界最佳 近的奇异应力场。研究结果表明,该轴对称界面端的应力奇异性与平面应变状态下相应模型的应力奇异性完全相同,材料对界面端附近奇异应力场的影响可用丰个双材料组合参数描述。     

Determination of effective thermomechanical parameters of a mixture of two elastothermoviscoplastic constituents

We analyze plane strain deformations of a representative volume element (RVE) to evaluate effective thermophysical parameters of a particulate composite comprised of two perfectly bonded heat conducting elasto-thermo-visco-plastic constituents. It is assumed that the composite is also isotropic and its response elasto-thermo-visco-plastic. Effective values of material parameters so computed are compared with those obtained from either existing micromechanics models or the rule of mixtures or both. It is found that values computed from the rule of mixtures differ at most by 10% from those obtained by using the RVE. Effective stress versus effective strain curves obtained by analyzing simple shearing and axisymmetric deformations of the RVE and of the homogenized material, and also those obtained in plane strain deformations involving loading/unloading/reloading are found to be very close to each other. Time histories of the effective plastic strain at two neighboring points, one in each constituent, are quite different. The effective stress computed by the rule of mixtures from the average effective stress in each constituent and its volume fraction is very close to that obtained from surface tractions acting on the specimen boundaries. The average effective stress in a constituent is computed from the effective plastic strain averaged over that constituent. This also holds for a composite comprised of three constituents.     

Ti-6Al-4V在高应变率下的动态剪切特性及失效机理

采用基于霍普金森压杆的新型加载技术对Ti-6Al-4V材料的动态剪切特性及失效机理进行了测试研究。获得了Ti-6Al-4V材料在超过104 s-1应变率下的剪应力-剪应变曲线及失效参数。研究发现,材料的流动应力存在明显的应变率强化效应;随着应变率的增加,材料的失效应力逐渐增大,而失效应变逐渐减小。采用ABAQUS/Explicit对加载过程进行了数值模拟。结果显示,剪切区材料基本处于平面剪切状态,应力应变场分布较为均匀,计算得到的剪应力-剪应变曲线与实验结果吻合较好。经断口分析可知,随着应变率的升高,Ti-6Al-4V的失效机理存在由韧窝、拉伸韧窝至台阶及河流花样的演化过程,材料的失效模式主要表现为韧性断裂。     

Utilizing phase retard integration method of flow birefringence to analyse the flow with symmetric plane

The formulas which are suitable to birefringent medium with symmetric plane are derived by means of phase retard integration. We have adopted this concept to the axisymmetric problems and deduced some useful formulas for these cases. As a practical application, the strain rate analysis of flow in a diverging or a converging vessel is illustrated at the end of this paper. The project supported by National Natural Science foundation of China     

纤维端部的界面裂纹分析

基于弹性力学空间轴对称问题的通解,研究了短纤维增强复合材料中纤维端部的轴对称币形和柱形界面裂纹尖端的应力奇异性,得到了裂纹尖端附近的奇异应力场．研究结果表明,这两种轴对称界面裂纹尖端的应力奇异性相同,并且与平面应变状态下相应模型的应力奇异性一致,材料性能对裂纹尖端附近奇异应力场的影响可用三个组合参数描述     

Calculation of an axisymmetric turbulent wall jet over a surface of convex curvature

One of the many applications of curved wall jets of engineering importance is the Coanda Flare, which is used for burning waste gases in the petroleum industry and which gave rise to this work. The gas jet flows over an axisymmetric tulip shaped body, entraining ambient air and so promoting clean combustion. The object of this work was to calculate the development of the jet with the extra rates of strain imposed by both longitudinal curvature and divergence. A differential ‘partially-parabolic’ technique was used with uncoupling of the streamwise and cross-stream momentum equations, leading to an efficient computer program. The extra rates of strain were modelled by corrections to a mixing length model with the two effects being assumed to be additive. The calculation method was compared with seven test ccases of experimental data. The first five were from published literature, and included the plane wall jet and axisymmetric free jet, and the separate effects of longitudinal curvature and divergence. The lost two cases were measurements of the wall jet flow over a model Coanda Flare. The calculation method gave generally good results for the main features of the flow such as growth rate and velocity decay. Details of the flow were not so well predicted, particularly the turbulent shear stress, as a result of the relatively simple turbulence model employed. The calculation method should provide a useful engineering tool, but some profitable developments could be made, particularly in the area of turbulence modelling.     

受压橡胶海绵圆筒的轴对称平面应变问题研究

针对橡胶海绵材料圆筒受压问题,根据欧拉-拉格朗日变分原理,建立了基于Blatz-Ko材料模型的轴对称平面应变问题的非线性微分方程. 采用参数变换的方法,获得了该问题的参数形式解析解. 并通过相应数值算例,得到了径向应力和切向应力沿圆筒径向的变化规律,以及过盈量对径向应力和切向应力影响规律.      

Similarity of three-dimensional and axisymmetric chemically-nonequilibrium flows in the neighborhood of a plane of symmetry

Three-dimensional chemically-nonequilibrium flow past blunt bodies in the neighborhood of the plane of symmetry is investigated within the framework of viscous shock layer theory. The similarity of three-dimensional and axisymmetric flows, previously established in [1] for a uniform gas, is extended to chemically-nonequilibrium gas flows. It is shown that the problem of determining the heat fluxes and friction stress in the neighborhood of the line of flow divergence can be reduced to the problem of determining these quantities for the axisymmetric body. The validity of the axisymmetric analogy is verified by carrying out numerical calculations for bodies of different shapes re-entering the earth's atmosphere along a gliding trajectory. Various models of surface catalytic activity are considered. The use of similarity relations makes it possible to apply existing programs for calculating axisymmetric flows to the solution of three-dimensional problems.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 115–120, March–April, 1990.     

Marangoni effect in the Couette flow

Viscous heating in an axisymmetric creeping flow of a second-order fluid with free surface between two coaxially mounted cylinders produces a radial temperature gradient in the fluid. The dependence of the surface tension upon temperature is the cause for a secondary flow in the meridional plane of the flow field. This secondary flow (Marangoni effect), and its influence upon the shape of the free surface are studied. The deformation of the free surface caused by the Marangoni effect is compared with the deformation caused by inertia and normal stress differences.     

复杂应力状态下岩土材料非线性本构模型研究

基于统一强度理论,将洛德参数引入统一强度理论,推导出材料的统一强度参数,进而分析中间主应力以及主剪应力系数对材料统一强度参数的影响。在此基础上,对邓肯张双曲线模型进行改进,使得该模型能够反映复杂应力状态下的应力应变关系。采用粘土的平面应变试验进行验证,结果表明:在平面应变情况下,当b=0.4时,模型能够较好地反映粘土的应力应变关系,及其强度参数。分析了洛德参数以及中间主应力系数对模型的影响,进一步说明岩土材料存在主应力效应。该模型能够反映不同材料在复杂应力状态下的应力应变关系,有其更为广泛的适用性。     

A combined necking and shear localization analysis for aluminum sheets under biaxial stretching conditions

A combined necking and shear localization analysis is adopted to model the failures of two aluminum sheets, AA5754 and AA6111, under biaxial stretching conditions. The approach is based on the assumption that the reduction of thickness or the necking mode is modeled by a plane stress formulation and the final failure mode of shear localization is modeled by a generalized plane strain formulation. The sheet material is modeled by an elastic-viscoplastic constitutive relation that accounts for the potential surface curvature, material plastic anisotropy, material rate sensitivity, and the softening due to the nucleation, growth, and coalescence of microvoids. Specifically, the necking/shear failure of the aluminum sheets is modeled under uniaxial tension, plane strain tension and equal biaxial tension. The results based on the mechanics model presented in this paper are in agreement with those based on the forming limit diagrams (FLDs) and tensile tests. When the necking mode is suppressed, the failure strains are also determined under plane strain conditions. These failure strains can be used as guidances for estimation of the surface failure strains on the stretching sides of the aluminum sheets under plane strain bending conditions. The estimated surface failure strains are higher than the failure strains of the forming limit diagrams under plane strain stretching conditions. The results are consistent with experimental observations where the surface failure strains of the aluminum sheets increase significantly on the stretching sides of the sheets under bending conditions. The results also indicate that when a considerable amount of necking is observed for a sheet metal under stretching conditions, the surface failure strains on the stretching sides of the sheet metal under bending conditions can be significantly higher.     

Mixed mode axisymmetric annular cracks in a finite layer: Off angle crack initiation

The solutions of axisymmetric Volterra type climb and glide edge dislocations are obtained in a layer by means of the Hankel transforms. Utilizing the same procedure, Green’s function solution is obtained for a layer under self-equilibration normal ring traction. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks where the layer is under axisymmetric normal loads. These equations are solved numerically to obtain dislocation density on the cracks surfaces. The results are employed to determine stress intensity factors for annular and penny-shaped cracks and the interaction between two co-axial penny-shaped cracks is studied. Moreover, the stress intensity factors of the interacting cracks are determined such that they can be further used in conjunction with strain energy density (SED) failure criterion to obtain the possible direction of crack initiation that may not be apparent under mixed mode conditions.