圆柱形爆炸容器绝热剪切瞬态失效过程

开展了圆柱形爆炸容器逐级加载和破坏实验,根据容器最终的断裂面和微观形貌观测,提出了爆炸容器绝热剪切失效模式.建立了应变率-应变空间内的绝热剪切损伤演化模型,将绝热剪切不同演化阶段的临界状态与宏观的力学条件联系起来,并将这些力学临界条件作为动态失效准则引入到宏观计算程序中,模拟爆炸容器发生绝热剪切的的瞬态过程,模拟结果成...     

Evolution of shear bands in bulk metallic glasses under dynamic loading

Shear band spacing in Zr-based bulk metallic glasses (BMGs) under dynamic loads is found to vary with position and local strain rate in the indented region. To investigate the dependence of shear band evolution characteristics on local strain rate and normal stress, a micromechanical model based on momentum diffusion is proposed. The thermo-mechanical model takes into account the normal stress dependence of yield stress, the free volume theory and the associated viscosity change within the shear band region. Temperature rise is obtained from the balance between the heat diffusion to the adjacent regions from a shear band and the heat generation due to the accumulated plastic work in a shear band. The parametric study has revealed that thermal effects play a minor role when the critical shear displacement is below 10 nm (as in nanoindentation) but become significant when the shear displacement accumulated in a shear band is of the order of hundreds of nanometers (as in uniaxial compression and in dynamic indentations). Finally, it is found that the normal stress plays a crucial role in the deformation behavior of BMGs by not only decreasing the time for shear band formation but also increasing the temperature rise significantly.     

Effect of viscoplastic relations on the instability strain,shear band initiation strain,the strain corresponding to the minimum shear band spacing,and the band width in a thermoviscoplastic material

We study thermomechanical deformations of a steel block deformed in simple shear and model the thermoviscoplastic response of the material by four different relations. We use the perturbation method to analyze the stability of a homogeneous solution of the governing equations. The smallest value of the average strain for which the perturbed homogeneous solution becomes unstable is called the critical strain or the instability strain. For each one of the four viscoplastic relations, we investigate the dependence upon the nominal strain-rate of the critical strain, the shear band initiation strain, the shear band spacing and the band width. It is found that the qualitative responses predicted by the Wright–Batra, Johnson–Cook and the power law relations are similar but these differ from that predicted by the Bodner–Partom relation. The computed band width is found to depend upon the specimen height.     

金属成形过程中的塑性失稳

本文评述了金属成形过程中塑性失稳现象(剪切带、颈缩、皱曲)研究的最新进展和仍存在的问题,着重讨论了塑性失稳预报和本构关系选择的密切关系,指出了塑性失稳预报结果取决于所选择的本构关系。还讨论了材料各向异性、应变率敏感性、不均匀性及其幅度、惯性等因素对塑性失稳的影响。除此之外,对作为分析塑性失稳的基本理论工具——Hill分叉准则亦作了简要讨论。      

Comparison of bifurcation and imperfection analyses of localized necking in rate-independent polycrystalline sheets

Plastic flow localization into a narrow band is investigated for a rate-independent polycrystalline sheet. Strain localization is considered to initiate at a point of bifurcation from the homogeneous deformation state or to result from the growth of a preexisting imperfection. Bifurcation analysis based on the rate-independent crystal plasticity model is found to give limit strains at a realistic strain level in the biaxial-stretching range. Limit strains predicted by the imperfection analysis are lower than those for bifurcation analysis and tend to approach those strain levels with decreasing amount of initial imperfection. Then, the bifurcation analysis is applied to the simulation of the forming limit for strongly textured sheets. The bifurcation analysis reveals the same trend of a texture dependence of limit strains as that for the imperfection analysis employing a rate-dependent model reported in the literature.     

Localization and shear-crippling (kinkband) instability in a thick imperfect laminated composite ring under hydrostatic pressure

A fully nonlinear finite elements analysis for prediction of localization representing shear-crippling (kinkband) instability in a thick laminated composite (plane strain) ring (infinitely long cylindrical shell) under applied hydrostatic pressure is presented. The primary accomplishment of the present investigation is prediction of meso(lamina)-structure-related equilibrium paths, which are often unstable in the presence of local imperfections and/or material nonlinearity, and which are considered to “bifurcate” from the primary equilibrium paths, representing periodic buckling patterns pertaining to global or structural level stability of the thick cross-ply ring with modal or harmonic imperfection. The present nonlinear finite elements solution methodology, based on the total Lagrangian formulation, employs a quasi-three-dimensional hypothesis, known as layerwise linear displacement distribution theory (LLDT) to capture the three-dimensional interlaminar (especially, shear) deformation behavior, associated with the localized interlaminar shear-crippling failure.A thick laminated composite [90/0/90] imperfect (plane strain) ring is investigated with the objective of analytically studying its premature compressive failure behavior. Numerical results suggest that interlaminar shear/normal deformation (especially, the former) is primarily responsible for the appearance of a limit (maximum pressure) point on the post-buckling equilibrium path associated with a periodic (modal or harmonic) buckling pattern, for which a modal imperfection serves as a perturbation. Localization of the buckling pattern results from “bifurcation” at or near this limit point, and can be viewed as a symmetry breaking phenomenon.In order to investigate a localization of the buckling pattern, a local or dimple shaped imperfection superimposed on a fixed modal one is selected. With the increase of local imperfection amplitude, the limit load (hydrostatic pressure) decreases, and also the limit point appears at an increased normalized deflection. Additionally, the load–deflection curves tend to flatten (near-zero slope) to an undetermined lowest pressure level, signaling the onset of “phase transition” in the localized region, and coexistence of two “phases”, i.e., a highly localized band of shear crippled (kinked) phase and its unshear-crippled (unkinked) counterpart along the circumference of the ring. Interlaminar shear-crippling triggered by the combined effect of imperfection, material nonlinearity and interlaminar shear/normal deformation appears to be the dominant compressive failure mode. A three-dimensional or quasi-three-dimensional theory, such as the afore-mentioned LLDT is essential in order to capture the meso-structure-related instability failure such as localization of the interlaminar shear crippling, triggered by the combined presence of local imperfection and material nonlinearity.     

Buckling driven debonding in sandwich columns

A compression loaded sandwich column that contains a debond is analyzed using a geometrically non-linear finite element model. The model includes a cohesive zone along one face sheet/core interface whereby the debond can extend by interface crack growth. Two geometrical imperfections are introduced; a global imperfection of the sandwich column axis and a local imperfection of the debonded face sheet axis. The model predicts the sandwich column to be very sensitive to the initial debond length and the local face sheet imperfection. The study shows that the sensitivity to the face sheet imperfection results from two mechanisms: (a) interaction of local debond buckling and global buckling and (b) the development of a damaged zone at the debond crack tip. Based on the pronounced imperfection sensitivity, the author predicts that an experimental measurement of the strength of sandwich structures may exhibit a large scatter caused by geometrical variations between test specimens.     

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.     

Prediction of the initial thickness of shear band localization based on a reduced strain gradient theory

Plastic flow localization in ductile materials subjected to pure shear loading and uniaxial tension is investigated respectively in this paper using a reduced strain gradient theory, which consists of the couple-stress (CS) strain gradient theory proposed by Fleck and Hutchinson (1993) and the strain gradient hardening (softening) law (C–W) proposed by Chen and Wang (2000). Unlike the classical plasticity framework, the initial thickness of the shear band and the strain rate distribution in both cases are predicted analytically using a bifurcation analysis. It shows that the strain rate is obviously non-uniform inside the shear band and reaches a maximum at the center of the shear band. The initial thickness of the shear band depends on not only the material intrinsic length l_cs but also the material constants, such as the yield strength, ultimate tension strength, the linear hardening and softening shear moduli. Specially, in the uniaxial tension case, the most possible tilt angle of shear band localization is consistent qualitatively with the existing experimental observations. The results in this paper should be useful for engineers to predict the details of material failures due to plastic flow localization.     

饱和土中剪切带宽度的研究

主要对饱和土在简单剪切和耦合载荷下的剪切带宽进行了理论分析. 在模型中考虑了剪胀、孔压耗散、惯性及孔隙水和土颗粒之间阻力等因素的耦合效应,将有效应力考虑为应变、应变率和孔隙压力的函数. 通过对模型的简化分析得到简单剪切和耦合应力状态下的剪切带宽度的近似计算公式. 分析表明,当土体失稳后,会出现有限厚度的剪切带, 且其形成是一个后期过程,可以忽略惯性效应. 剪切带的厚度依赖于孔隙压力特性和土体骨架的剪胀特性,但是不依赖应力状态. 最后将分析与实验结果进行了对比,在一定程度上,两者是符合的.     

Shear bands in dense metallic granular materials

In this paper, the phenomenon of strain localization, i.e., shear banding, in densely distributed metallic assemblies has been studied. A discrete element methodology for analyzing metallic granular materials has been put forward. In this numerical model, elastoplastic contact, as well as friction, rolling resistance and cohesion between spheres, are explicitly taken into account. The calculations reveal that the shear banding mechanism in dense assemblies can be thought as an instability triggered by initial imperfections. Within a band, the motion, deformation and rearrangement of spheres soften the resistance of the aggregate, as these mechanisms create additional geometric imperfections. Additionally, the simulations showed that the shear-band width does not change conclusively with the friction, rolling resistance and plasticity parameters. However, cohesive strength, even in small amounts, drastically increased the shear-band width. Finally, the shear-band thickness and inclination angles are strongly dependent on the degree of initial imperfection. Whereas for a perfect assembly the shear band inclinations were consistently around 60°, more heterogeneous assemblies lead to shear band angles closer to the continuum mechanics solution, which is 45°. This was found to be in agreement with recent experimental observations.     

Effects of texture on shear band formation in plane strain tension/compression and bending

In this study, effects of typical texture components observed in rolled aluminum alloy sheets on shear band formation in plane strain tension/compression and bending are systematically studied. The material response is described by a generalized Taylor-type polycrystal model, in which each grain is characterized in terms of an elastic–viscoplastic continuum slip constitutive relation. First, a simple model analysis in which the shear band is assumed to occur in a weaker thin slice of material is performed. From this simple model analysis, two important quantities regarding shear band formation are obtained: i.e. the critical strain at the onset of shear banding and the corresponding orientation of shear band. Second, the shear band development in plane strain tension/compression is analyzed by the finite element method. Predictability of the finite element analysis is compared to that of the simple model analysis. Third, shear band developments in plane strain pure bending of a sheet specimen with the typical textures are studied. Regions near the surfaces in a bent sheet specimen are approximately subjected to plane strain tension or compression. From this viewpoint, the bendability of a sheet specimen may be evaluated, using the knowledge regarding shear band formation in plane strain tension/compression. To confirm this and to encompass overall deformation of a bent sheet specimen, including shear bands, finite element analyses of plane strain pure bending are carried out, and the predicted shear band formation in bent specimens is compared to that in the tension/compression problem. Finally, the present results are compared to previous related studies, and the efficiency of the present method for materials design in future is discussed.     

Non-uniform deformation after prestrain

The deformation behaviour of prestrained metal sheets is analysed in this work. The non-uniform deformation observed during reloading in tension was studied, by following deformation in different regions of the samples. It takes into account the presence of geometrical defects in the samples and explains the importance of mechanical behaviour. A simplified analysis was used, to model the behaviour in tension of a metallic specimen with geometrical imperfection. The flow behaviour is described using a Swift law equation, which includes strain-rate sensitivity. A modified law was used for prestrained materials and this incorporates the plastic prestrain value, adjusted to the path change. The model predicts imperfection growth kinetics with strain, and strain saturation in the homogeneous region, due to the onset of necking.     

Shear instability direction at a crack-tip in a thermoviscoplastic body

Deformation fields near a crack tip are analyzed for a thermoviscoplastic body deformed either in an antiplane shear or in plane strain. The effects of inertia forces and heat conduction are considered, but those of material elasticity are neglected. By assuming that a shear band initiating from the crack-tip propagates in the direction of the maximum effective stress. It is found that a shear band propagates along the crack ligament in the body deformed in an antiplane shear. For a body undergoing plane strain deformations, the direction of the propagation of the shear band depends upon the mode-mixity parameter and agrees with that observed by Kalthoff.     

Elongational flow behavior of polymeric fluids according to the Leonov model

The viscoelastic behavior of polymeric systems based upon the Leonov model has been examined for (i) the stress growth at constant strain rate, (ii) the stress growth at constant speed and (iii) the elastic recovery in elongational flow. The model parameters have been determined from the available rheological data obtained either in steady shear flow (shear viscosity and first normal-stress difference as a function of shear rate) or oscillatory flow (storage and loss moduli as a function of frequency in the linear region) or from extensional flow at very small strain rates (time-dependent elongation viscosity in the linear viscoelastic limit). In addition, the effect of the parameter characterizing the strain-hardening of the material during elongation has also been studied. The estimation of this parameter has been based upon the structural characteristics of the polymer chain which include the critical molecular weight and molecular weight of an independent segment. Five different polymer melts have been considered with varying number of modes (maximum four modes). Resulting predictions are in fair agreement with corresponding experimental data in the literature.     

表面粗糙度对TC4钛合金柱壳剪切带形成的影响

剪切带是材料在高应变率加载条件下特有的变形和损伤形式之一,关于影响金属材料中剪切带形成的敏感性因素及其机理的研究,一直是科学研究和工程设计中关注的重点问题. 在柱壳高速坍塌过程中,剪切带优先在内表面形核, 其形核及扩展行为受内表面介观状态的影响显著.本文采用爆轰加载厚壁圆筒坍塌实验技术,结合材料表面处理技术、微结构表征技术和剪切带理论模型分析,研究了内表面粗糙度变化对TC4钛合金柱壳剪切带形成影响的细观动力学规律.结果表明, 在爆炸加载形成的高应变率条件下,表面粗糙度对TC4钛合金柱壳中剪切带形成具有明显影响. 在相同的变形条件下,随着试样内表面粗糙度的增大, 剪切带数量、长度和形核速率均增大;表面粗糙度越大, 部分剪切带扩展速率越快, 剪切带长度差异越大,剪切带的屏蔽效应增强. 分析表明,实验获得的剪切带间距与W-O模型和M模型预测结果基本吻合,具体数值受试样内表面粗糙度影响, 随着表面粗糙度的增大,实验结果逐渐小于预测数值.      

Strain gradient crystal plasticity effects on flow localization

In metal grains one of the most important failure mechanisms involves shear band localization. As the band width is small, the deformations are affected by material length scales. To study localization in single grains a rate-dependent crystal plasticity formulation for finite strains is presented for metals described by the reformulated Fleck–Hutchinson strain gradient plasticity theory. The theory is implemented numerically within a finite element framework using slip rate increments and displacement increments as state variables. The formulation reduces to the classical crystal plasticity theory in the absence of strain gradients. The model is used to study the effect of an internal material length scale on the localization of plastic flow in shear bands in a single crystal under plane strain tension. It is shown that the mesh sensitivity is removed when using the nonlocal material model considered. Furthermore, it is illustrated how different hardening functions affect the formation of shear bands.     

The microbuckling failure of Dyneema~? composites under compression

Two grades of Dyneema~?composite laminates with the commercial designations of HB26 and HB50 were cut into blocks with or without an edge crack and compressed in the longitudinal fiber direction. The cracked and uncracked specimens show similar compressive responses including failure pattern and failure load. The two grades of Dyneema~? composites exhibits different failure modes: a diffuse, sinusoidal buckling pattern for Dyneema~? HB50 due to its weak matrix constituent and a kink band for Dyneema~? HB26 due to its relatively stronger matrix constituent. An effective finite element model is used to simulate the collapse of Dyneema~? composites, and the sensitivity of laminate compressive responses to the overall effective shear modulus, interlaminar shear strength, thickness and imperfection angle are investigated. The change of failure mode from kink band to sinusoidal buckling pattern by decreasing the interlaminar shear strength is validated by the finite element analyses.     

Micromechanics based second gradient continuum theory for shear band modeling in cohesive granular materials following damage elasticity

Gradient theories, as a regularized continuum mechanics approach, have found wide applications for modeling strain localization failure process. This paper presents a second gradient stress–strain damage elasticity theory based upon the method of virtual power. The theory considers the strain gradient and its conjugated double stresses. Instead of introducing an intrinsic material length scale into the constitutive law in an ad hoc fashion, a microstructural granular mechanics approach is applied to derive the higher-order constitutive coefficients such that the internal length scale parameter reflects the natural granularity of the underlying material microstructure. The derivations of the required damage constitutive relationships, the strong form governing equations as well as its weak form for the second gradient model are described. The recently popularized Element-Free Galerkin (EFG) method is then employed to discretize the weak form equilibrium equation for accommodating the resultant higher-order continuity requirements and further handling the mesh sensitivity problem. Numerical examples for shear band simulations show that the proposed second gradient continuum model can produce stable, accurate as well as mesh-size independent solutions without a priori assumption of the shear band path.     

Roles of plastic spin in shear banding

In this paper, the effects of plastic spin on shear banding and simple shear are examined systematically. Three types of plastic constitutive model with plastic spin are considered: (i) a non-coaxial model in which the direction of the plastic strain rate depends on that of the stress rate; (ii) a strain-softening model based on the J₂ flow theory; and (iii) the pressure-sensitive porous plasticity model. All the constitutive models are formulated in viscoplastic forms and in conjunction with non-local concepts that have been recently focused and discussed. First, behavior in simple shear is examined by numerical analysee with the aforementioned constitutive models. Moreover, some experimental evidences for stress response to simple shear are shown; that is, several large torsion tests of metal tubes and bars are carried out. Next, finite element simulations of shear banding in plane strain tension are performed. A critical effect of plastic spin on shear banding is observed for the noncoaxial model, while an almost negligible effect is observed for the porous model. The identical effects of plastic spin are observed, whether nonlocality exists or not. Finally, we discuss the relationship between the behavior in simple shear and the shear band formation. It is emphasized that this is a critical issue in predicting shear banding in macroscopic grounds.