块体金属玻璃压缩变形和破坏特性的有限元模拟研究

论文通过LS-DYNA软件的用户自定义材料子程序(UMAT),将考虑自由体积、热以及静水应力影响的三维热力耦合本构模型编写入有限元软件中.同时考虑金属玻璃材料内部结构的不均匀性,建立了包含随机分布剪切带弱区的有限元模型.在此基础上开展了针对块体金属玻璃在准静态压缩和弯曲条件下的有限元模拟研究,具体分析材料的非均匀变形和破坏特性,并特别研究了相应的剪切带行为,包括剪切带的形核、传播以及剪切带诱导的断裂等.     

平板内后分叉的局部化带

已知剪切带状分叉将引发裂隙．为模拟其后继效应，本文利用大应变有限元的空单元技术设计了自动释放损伤单元内残余应力的方法．这种损伤松弛可以具有各向同性或各向异性的性质．结果显示，受单向拉伸平板中心的水平裂隙会向斜角方向辐射局部化带，带内有明显厚度颈缩及平面错动．由此说明实验所观察到的这类局限于带内的变形应属于后分叉的效应     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为。计算采用二维轴对称应变条件下的绝热模型。钨合金的本构方程采用热粘塑性形式的Johnson Cook模型。为了得到不同尺度的变形信息 ,计算中用了两种网格 ;先用粗糙网格分析试件变形局部化的概貌 ;接着 ,用细密网格 (在变形局部化区域 ,网格尺寸达到 10 m)分析绝热剪切带的形成和发展。有限元模拟得到的绝热剪切带位置和方向与实验一致。计算结果表明 ,绝热剪切带的形成和发展与试件的应力状态密切相关。     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模拟

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为.计算采用二维轴对称应变条件下的绝热模型.钨合金的本构方程采用热粘塑性形式的Johnson-Cook模型.为了得到不同尺度的变形信息,计算中用了两种网格;先用粗糙网格分析试件变形局部化的概貌;接着,用细密网格(在变形局部化区域,网格尺寸达到10μm)分析绝热剪切带的形成和发展.有限元模拟得到的绝热剪切带位置和方向与实验一致.计算结果表明,绝热剪切带的形成和发展与试件的应力状态密切相关.     

粘塑性损伤模型模拟准超塑性单轴拉伸行为

发展了Chaboche粘塑性本构模型的大变形隐式算法，用损伤(DM)和无损伤(NDM)模型模拟准超塑性单轴拉伸。发现变形过程可分为三个阶段：均匀变形、颈缩发展、断裂破坏阶段。DM可准确模拟前两个阶段变形，NDM只能较好地模拟均匀变形阶段，表明DM可以较精确地描述稳定发展的动态过程。由于有限元方法只能描述连续介质，因此对于断裂破坏阶段，NDM模拟载荷大于试验结果，DM的载荷小于试验结果，这是由高应变速率敏感性造成。DM能够描述试验中出现地多处颈缩现象，局部应变速率分布随时间演化反映了颈缩发展程度。严重颈缩部位的距离代表着超塑性变形能力，距离越大，抗颈缩能力越好。     

预延伸非晶聚合物材料各向异性变形局部化

给出一种可描述预延伸各向异性特性的背应力张量三维表达式，引入大变形弹塑性有限元驱动应力法，结合ＢＰＡ８ 链细观分子网络模型，模拟了预延伸各向异性非晶聚合物材料平面应变拉伸变形局部化力学行为．详细讨论了预延伸比（ＩｎｉｔｉａｌＤｒａｗｉｎｇＲａｔｉｏ；ＩＤＲ）和预延伸方向（ＩｎｉｔｉａｌＤｒａｗｉｎｇＤｉｒｅｃｔｉｏｎ；ＩＤＤ）对变形抗力、颈缩规律、剪切带方向以及试件中心部位链延伸比的影响．     

动态载荷下剪切变形局部化、微结构演化与剪切断裂研究进展

总结和评述了近年来金属与合金变形局部化的形成、微结构演化与剪切断裂方面作者和相关的研究工作成果. 材料包括低碳钢, SS304不锈钢, Fe-15%Ni-15%Cr单晶, Al-Li合金,α-Ti和Ti-6Al-4V, Al/SiCp复合材料等.综述内容主要包括:采用改进的Hopkinson扭杆装置,对剪切变形局部化形成、发展和演化过程进行了实验观察与数值模拟;采用"侧剖"与"对接"等定点方法制备电子显微镜薄膜试样,对剪切带内相变与再结晶、非晶转变、旋涡结构等动态变形现象,以及与宏观动态力学行为对应的位错胞的形成、发展和坍塌等微结构特征进行了观测;提出了应变和应变率同时作为剪切带形成的两个必要条件的直接实验证据;在剪切带内发现了α'$-马氏体相变现象,以及相变产物与母体之间的晶体学关系;通过位错单滑移或交滑移等微观剪切最后发展成为宏观剪切的机制;对剪切带内再结晶结构的观测和对再结晶动力学本构关系的定量描述;对剪切带特别是``白色'腐蚀带(或相变带)的形成机制的分析和新的解释,指出 ``白色'是带内亚结构取向趋于一致,其在光学或扫描显微镜下很难辨认这些亚结构的取向差所致,并非表明剪切带内一定发生了相变;通过截断实验和实时跟踪观测发现,剪切带内微裂纹的萌生与聚合是材料承载能力骤然下降并导致最后断裂的主控因素.此外,本文对近年来在准静态和循环加载下材料的局部化形变与剪切断裂的实验结果予以简要评述,指出其微观机制与动态载荷下的截然不同, 是由位错的平面滑移所控制的,与热效应无关的等温变形.      

延性金属膨胀环的运动及断裂的三维数值模拟

利用LS-DYNA三维动力有限元软件对延性金属环的膨胀运动与断裂进行数值模拟。在膨胀环圆周加入泊松随机分布断裂成核点,利用J-C本构模型,研究诸如颈缩形成时间、颈缩区与均匀变形区的温度、应力、应变的对比等颈缩形成机理,以及讨论了环圆周上断裂成核点的泊松随机分布碎裂特性的影响。数值计算结果与实验结果、理论分析结果吻合较好,表明施加泊松随机断裂成核点的数值模拟方法是合理的。     

内爆炸载荷作用下7A55铝合金的动态性能及断裂行为

采用圆筒爆炸实验研究了内爆炸载荷作用下7A55铝合金的动态性能,用扫描电镜和光学显微镜 对破裂样品的断口形貌、金相组织等进行了微观分析。结果表明:在本文实验条件下,7A55铝合金能够承受 360MPa的内部爆炸加载;合金的断裂方式为剪切断裂,裂缝与筒壁的径向近似成45角;靠近圆筒内侧组织 中存在剪切变形带、绝热剪切带和裂纹,沿最大剪切应力面向外扩展。     

TC4在动态载荷下的剪切行为研究

使用分离式霍布金森压杆(SHPB)对2种TC4(Ti-6Al-4V)试样(单边剪切试样与双边剪切试样)在应变率104 s-1下进行动态剪切加载，利用SIM D8高速照相系统捕捉了绝热剪切带扩展的整个历程，得到了TC4在拍照时刻的应力应变曲线；使用金相显微镜和SEM扫描电镜对TC4绝热剪切带的微观形貌进行观察，发现绝热剪切带宽度为5~12 μm，断口从韧窝断裂演变为解理断裂，可观测到韧窝状与河流花样断口形貌，但是并未看到相变的发生；对2种试样就产生绝热剪切带的形式与敏感性进行了分析，实验表明双边试样更易产生绝热剪切带；通过高速照相系统的标定换算，得到TC4绝热剪切带产生的临界剪切应变在78%~88%之间。在SHPB动态加载条件下，TC4绝热剪切带的扩展速度在460~1 250 m/s之间，且应变率越高，剪切带扩展越快，扩展平均速度与名义应变率近似呈线性关系；另外，在同一加载速率下，剪切带并不是匀速扩展，其扩展速度随载荷的增加而不断增加。     

Hydrogen induced shear localization of the plastic flow in metals and alloys

Hydrogen enhanced localized plasticity (HELP) is a viable mechanism for hydrogen embrittlement supported by experimental observations. According to the HELP mechanism, hydrogen induced premature failures result from hydrogen induced plastic instability which leads to hydrogen assisted localized ductile processes. The objective of this work is to reveal the role of hydrogen in possibly localizing the macroscopic deformation into bands of intense shear using solid mechanics methodology. The hydrogen effect on material deformation is modeled through the hydrogen induced volume dilatation and the reduction in the local flow stress upon hydrogen dissolution into the lattice. Hydrogen in assumed to reside in both normal interstitial lattice sites (NILS) and reversible traps associated with the plastic deformation. The analysis of the plastic deformation and the conditions for plastic flow localization are carried out in plane strain uniaxial tension. For a given initial hydrogen concentration in the unstressed specimen, a critical macroscopic strain is identified at which shear localization commences.     

Influence of void nucleation on ductile shear fracture at a free surface

An approximate continuum model of a ductile, porous material is used to study the influence of the nucleation and growth of micro-voids on the formation of shear bands and the occurrence of surface shear fracture in a solid subject to plane strain tension. Bifurcation into diffuse modes is analysed for a plane strain tensile specimen described by these constitutive relations, which account for a considerable plastic dilatancy due to void growth and for the possibility of non-normality of the plastic flow law. In particular, bifurcation into surface wave modes and the possible influence of such modes triggering shear bands is investigated. For solids with initial imperfactions such as a surface undulation, a local material inhomogeneity on an inclusion colony, the inception and growth of plastic flow localization is analysed numerically. Both the formation of void-sheets and the final growth of cracks in the shear bands is described numerically. Some special features of shear band development in the solid obeying non-normality are studied by a simple model problem.     

高分子材料平面应变拉伸变形局部化

将基于应变软化玻璃状高分子材料微观特征建立的ＢＰＡ８－链分子网络模型引入ＵｐｄａｔｉｎｇＬａｇｒａｎｇｅ有限元方法，建立了适于变形局部化分析的大变形弹塑性有限元驱动应力法．在此基础上，数值模拟了初始各向同性高分子材料平面应变拉伸变形局部化的传播过程．探讨了ＢＰＡ模型对具有加工硬化特性的结晶性高分子材料变形分析的适应性；分析了局部化传播过程中颈缩截面的非均匀应力三轴效应；最后，讨论了网格尺寸以及初始几何不均匀性对颈缩扩散以及应力三轴效应的影响     

On neck propagation in amorphous glassy polymers under plane strain tension

Unlike metals, necking in polymers under tension does not lead to further localization of deformation, but to propagation of the neck along the specimen. Finite element analysis is used to numerically study necking and neck propagation in amorphous glassy polymers under plane strain tension during large strain plastic flow. The constitutive model used in the analyses features strain-rate, pressure, and temperature dependent yield, softening immediately after yield and subsequent orientational hardening with further plastic deformation. The latter is associated with distortion of the underlying molecular network structure of the material, and is modelled here by adopting a recently proposed network theory developed for rubber elasticity. Previous studies of necking instabilities have almost invariably employed idealized prismatic specimens; here, we explicitly account for the unavoidable grip sections of test specimens. The effects of initial imperfections, strain softening, orientation hardening, strain-rate as well as of specimen geometry and boundary conditions are discussed. The physical mechanisms for necking and neck propagation, in terms of our constitutive model, are discussed on the basis of a detailed parameter study.     

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.     

Predicting failure modes and ductility of dual phase steels using plastic strain localization

Ductile failure of metals is often treated as the result of void nucleation, growth and coalescence. Various criteria have been proposed to capture this failure mechanism for various materials. In this study, ductile failure of dual phase steels is predicted in the form of plastic strain localization resulting from the incompatible deformation between the harder martensite phase and the softer ferrite matrix. Microstructure-level inhomogeneity serves as the initial imperfection triggering the instability in the form of plastic strain localization during the deformation process. Failure modes and ultimate ductility of two dual phase steels are analyzed using finite element analyses based on the actual steel microstructures. The plastic work hardening properties for the constituent phases are determined by the in-situ synchrotron-based high-energy X-ray diffraction technique. Under different loading conditions, different failure modes and ultimate ductility are predicted in the form of plastic strain localization. It is found that the local failure mode and ultimate ductility of dual phase steels are closely related to the stress state in the material. Under plane stress condition with free lateral boundary, one dominant shear band develops and leads to final failure of the material. However, if the lateral boundary is constrained, splitting failure perpendicular to the loading direction is predicted with much reduced ductility. On the other hand, under plane strain loading condition, commonly observed necking phenomenon is predicted which leads to the final failure of the material. These predictions are in reasonably good agreement with experimental observations.     

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.     

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.     

岩石试件端面摩擦效应数值模拟研究

试件端面摩擦效应直接影响试件内的塑性等效应变、侧向位移的分布和单元应力应变曲线。本文运用ANSYS中的接触单元模拟了平面应变状态下端面摩擦效应对塑性等效应变、侧向位移和单元应力应变曲线的影响,得到了不同摩擦系数时塑性等效应变及侧向位移的渐进变化形式。当接触面摩擦较小时,塑性等效应变图案为上下两个X形网络,侧向位移上下分布均匀;当接触面摩擦增大时,塑性等效应变网络向中部靠拢并且明显增大,侧向位移上下分布不均匀,中部较上下端面位移大;当试件端面侧向位移被限制,即摩擦力很大时,塑性等效应变网络变为一个X形局部化带,侧向位移分布更加不均匀,中部明显隆起。