Modelling and simulation of the localized plastic deformation by relaxation element method

In the present work a new and a little known method of modelling and simulation of plastic strain localization in a loaded solid - relaxation element method (REM) is presented. The fundamental property of solid: “Plastic deformation of solid under loading is accompanied by stress relaxation in the local volume of it” lies in the basis of the method. The theoretical grounding of the method from the point of view of constitutive equations of theory of elasticity and continuum theory of defects is conducted. For the plane-stress state the technique of construction of local sites of plastic deformation by the method on the mesoscopic scale level is demonstrated. The technique of calculation of stress state of the plane with the sites of different shape and gradients of plastic deformation is described. The new type of relaxation elements, different from ones used in the works [1-6] has been applied. This allowed to extend the method on the new class of objects. The examples and the results of simulation of the processes of plastic strain localization, accompanied by the effect of Lüder’s band propagation and Portevin Le Chatelier Effect are considered. The examples show, that REM is a comfortable tool when simulating the effect of plastic strain localization on meso- and macro- level.     

冲击块与金属圆柱壳相互作用的缓冲吸能分析

应用Johnson Cook本构方程, 结合Singace叠缩模型,考虑冲击作用引起材料的应变强化效应、应变率强化效应和温度效应,研究了冲击物体和金属圆柱壳相互作用的能量转化过程。运用角度增量叠缩法计算了圆柱壳的塑性变形能,根据冲击块和圆柱壳所组成系统能量守恒得到了冲击块速度位移、冲击块和圆柱壳之间瞬时载荷的解析表达式,以及圆柱壳塑性变形引起的温度分布函数。通过对不同材质的金属圆柱壳在冲击作用下塑性变形过程的计算分析,并和有限元计算结果比较,证明了本文计算方法的正确性。     

Plastic strain determination using peak speckle pattern intensity

This paper presents an integrated peak spectral intensity method for the measurement of plastic deformation in metallic materials. The method employs a laser speckle technique and a computer based image processing system. Test results on aluminium, brass and copper sheets subjected to bi-axial stretching using the hydroforming process shows a correlation between the integrated peak spectral intensity coefficient H_i and the equivalent strain. In addition, the use of the parameter H_i requires less computation and enables measurement of plastic strains of up to 30%.     

中锰Q&P钢残余奥氏体的稳定性及其TRIP效应

通过单向拉伸及平面应变实验研究了Mn含量为7%的中锰淬火-配分(QP)钢残余奥氏体的机械稳定性,利用X射线衍射仪(XRD)测定试验钢残余奥氏体的含量,通过观察试验钢的拉伸曲线及扫描电镜(SEM)、透射电镜(TEM)照片,分析变形前后的微观组织,研究中锰QP钢的变形机制。结果表明:应力状态对残余奥氏体稳定性有较大的影响,平面应变更有利于相变诱导塑性(TRIP)效应的发挥;中锰QP钢的拉伸变形特征是由超细晶硬化机制和TRIP效应相互作用产生的,通过微观组织观察发现中锰QP钢的塑性变形主要是残余奥氏体的TRIP效应,其中薄膜状的残余奥氏体的稳定性最高。     

Sensitivity of plastic strain localization zones to boundary and initial conditions

 The plastic strain localization phenomena are studied with the use of the theory of viscoplasticity and finite element method. The aim of the study is to check if the strain localization is sensitive to variation of initial and boundary conditions. Total and local plastic deformation energy as candidate response functionals for sensitivity analysis are discussed by selected numerical examples of wave propagation problems. The paper considers sensitivity of the form and intensity of plastic strain localization to the variation of the above parameters. Different response functionals for structures demonstrating plastic strain localization were proposed and studied. Received: 4 February 2002 / Accepted: 10 October 2002     

3004铝合金晶格和织构演化的原位X射线衍射研究

为探究3004铝合金板材在拉伸变形过程中晶粒取向变化和晶格应变行为,明确两者在塑性变形过程中的竞争机制和作用,本文首次将原位拉伸方法用于铝合金织构演化和晶格应变行为的研究中,利用试验装置对3004铝合金进行定量应变拉伸,并在不同变形量下对样品进行X射线衍射物相和织构测试分析.实验结果表明:样品在弹性变形区晶格常数随应变...     

A Lagrangian model for hardening behaviour of materials at finite deformation based on the right plastic stretch tensor

In this paper a modified multiplicative decomposition of the right stretch tensor is proposed and used for finite deformation elastoplastic analysis of hardening materials. The total symmetric right stretch tensor is decomposed into a symmetric elastic stretch tensor and a non-symmetric plastic deformation tensor. The plastic deformation tensor is further decomposed into an orthogonal transformation and a symmetric plastic stretch tensor. This plastic stretch tensor and its corresponding Hencky’s plastic strain measure are then used for the evolution of the plastic internal variables. Furthermore, a new evolution equation for the back stress tensor is introduced based on the Hencky plastic strain. The proposed constitutive model is integrated on the Lagrangian axis of the plastic stretch tensor and does not make reference to any objective rate of stress. The classic problem of simple shear is solved using the proposed model. Results obtained for the problem of simple shear are identical to those of the self-consistent Eulerian rate model based on the logarithmic rate of stress. Furthermore, extension of the proposed model to the mixed nonlinear isotropic/kinematic hardening behaviour is presented. The model is used to predict the nonlinear hardening behaviour of SUS 304 stainless steel under fixed end finite torsional loading. Results obtained are in good agreement with the available experimental results reported for this material under fixed end finite torsional loading.     

对泡沫塑料落锤冲击试验的解析

结合KBM法和能量法求解泡沫塑料试件受冲击的压力、变形、加载速度及应力应变关系，得到了相当准确的结果。同时，系统地阐述了粘性对泡沫塑料冲击试验的影响，为其工程应用提供了理论依据。     

SiCp/AZ91D镁基纳米复合材料的室温拉伸行为及塑性变形机理

为得到高强度和高塑性的镁基复合材料,通过高能超声分散法和金属型重力铸造工艺制备了SiC纳米颗粒分散均匀的SiCp/AZ91D镁基纳米复合材料,并进行T4固溶热处理和室温拉伸。利用扫描电子显微镜(SEM)、透射电子显微镜(TEM)对试样拉伸后的显微组织和塑性变形机理进行观察与研究。结果表明:T4态SiCp/AZ91D镁基纳米复合材料室温下抗拉强度达到296 MPa,伸长率达到17.3%。经室温拉伸变形后复合材料基体微观组织中出现了大量的孪晶和滑移,孪生和滑移是复合材料塑形变形的主要机制。在室温拉伸过程中,α-Mg基体中SiC纳米颗粒周围形成高应变场,高应变场内形成大量位错和堆垛层错,这些位错和堆垛层错在拉伸应变的作用下演变成大量的滑移带和孪晶,这是SiCp/AZ91D镁基纳米复合材料在室温下具有高塑性的微观塑性变形机理。     

Microscale characterization of granular deformation near a crack tip

This paper presents a study of microscale plastic deformation at the crack tip and the effect of microstructure feature on the local deformation of aluminum specimen during fracture test. Three-point bending test of aluminum specimen was conducted inside a scanning electron microscopy (SEM) imaging system. The crack tip deformation was measured in situ utilizing SEM imaging capabilities and the digital image correlation (DIC) full-field deformation measurement technique. The microstructure feature at the crack tip was examined to understand its effect on the local deformation fields. Microscale pattern that was suitable for the DIC technique was generated on the specimen surface using sputter coating through a copper mesh before the fracture test. A series of SEM images of the specimen surface were acquired using in situ backscattered electronic imaging (BEI) mode during the test. The DIC technique was then applied to these SEM images to calculate the full-field deformation around the crack tip. The grain orientation map at the same location was obtained from electron backscattered diffraction (EBSD), which was superimposed on a DIC strain map to study the relationship between the microstructure feature and the evolution of plastic deformation at the crack tip. This approach enables to track the initiation and evolution of plastic deformation in grains adjacent to the crack tip. Furthermore, bifurcation of the crack due to intragranular and intergranular crack growth was observed. There was also localization of strain along a grain boundary ahead of and parallel to the crack after the maximum load was reached, which was a characteristic of Dugdale–Barenblatt strip-yield zone. Thus, it appears that there is a mixture of effects in the fracture process zone at the crack tip where the weaker aspects of the grain boundary controls the growth of the crack and the more ductile aspects of the grains themselves dissipate the energy and the corresponding strain level available for these processes through plastic work.     

Development of a new integrated severe plastic deformation method

ABSTRACT

A novel one-pass integrated severe plastic deformation method, entitled extrusion compression angular pressing forward extrusion (ECAP-FE), was designed and used for fabricating a fine-grained as-cast AZ31 Mg alloy. Subsequently, mechanical (room temperature compression, tensile and microhardness) and microstructural properties of the processed sample were investigated in detail. In addition, a finite element simulation of the proposed method was carried out for evaluating the equivalent plastic strain and strain rate distribution. The results showed that the ECAP-FE method is a powerful method for processing the ultrafine-grained as-cast AZ31 Mg alloy in a single pass to achieve a uniform and fine microstructure with enhanced mechanical properties. Therefore, it appears that the proposed method has a great potential to a wide range of industrial applications.     

常温养护型超高性能混凝土的圆环约束收缩性能

研究了HCSA膨胀剂3种掺量(0%、3%、6%)下常温养护型超高性能混凝土(Ultra high performance concrete,UHPC)的圆环约束收缩性能,包括:(1)UHPC轴拉应力应变曲线测试;(2)根据GB/T50082的UHPC自由收缩实验;(3)根据ASTM C1581的UHPC圆环约束实验。结果表明,3种UHPC的极限拉伸应变均高于3 000με,28d总收缩值分别为1 005.6με、600.0με、462.2με,并且在圆环约束作用下转化为残余应变、弹性拉应变和塑性拉应变,其中塑性拉应变分别为700.4με、437.9με、389.9με。3种UHPC在拉伸应变达到1 000με时及圆环约束实验中均未发现0.01mm以上的可检测裂缝。基于拉伸实验和声发射损伤分析方法对UHPC进行应变分析,可知具有应变强化段的3种UHPC在圆环约束实验中的塑性变形以小于0.01mm的多点分布微裂纹形式存在。通过添加HCSA膨胀剂对常温养护型UHPC进行收缩补偿,可有效地降低UHPC自身的拉应力以及对原有结构的影响。     

基于iBeam3单元逆有限元法的冻土区管道变形研究

形状传感技术结合了应变-形状算法,使测点的应变与空间位置对应。该技术能真实地反映待测物的空间姿态与形状变化,可实现动态实时监测。该文提出了一种基于iBeam3单元的逆有限元法来实现管道变形状况的实时监测。建立了基于iBeam3单元的逆有限元法的理论框架模型;试验验证了其合理性;设计了冻土区管道变形试验,结合应变传感器所测表面应变信息建立逆有限元分析模型;实现了在整个冻融循环过程中对管道形状的重建。试验结果表明:该方法不仅能实时监测管道随土体冻融循环的完整变形过程,而且可以保证精度。在解决寒区管道结构健康监测问题上潜力巨大。     

High strain-rate localization and failure of crystalline materials

A theoretical and a computational model are introduced to study the micromechanical high strain-rate failure mechanisms of shear-strain localization in monocrystalline fcc structures. A theoretical framework for a constitutive model for the dynamic finite plastic deformation of rate-dependent single fcc crystals is developed. The micromechanics of plastic flow are based on a high strain-rate single crystal plasticity model and a visco-plastic power law. The single crystal is subjected to far-field dynamic tensile strain-rates ranging from 100/s to 2000/s. An explicit finite-element model is introduced for the integration of the numerically stiff visco-plastic constitutive relations. The total deformation rate tensor is obtained by the central difference explicit integration of the equations of motion. The plastic deformation rate tensor is obtained from the solution of an initial-value nonlinear problem for the resolved shear stresses. In time intervals when the differential equations for the resolved shear stresses are not numerically stiff, the initial-value problem is integrated by the explicit fifth-order Runge-Kutta adaptive time-step method. In time domains where the propagated error grows and the time-step must be restricted due to stability requirements, which is an indication of numerical stiffness, the initial-value problem is integrated by an A-stable method. To correctly differentiate time-step reductions due to stability, from time-step reductions due to accuracy, a stiffness ratio is defined. The present analysis corroborates experimental observations that high strain-rate shear-strain localization, in rate-dependent crystals, is a function of thermal and geometrical softening, overall strain-rates, strain hardening, strain-rate hardening, and strain-rate sensitivity.     

Heterogeneous plastic deformation in bimodal bulk ultrafine-grained nickel

Heterogeneous plastic deformation behavior of two bimodal ultrafine-grained nickel materials with different ultrafine-grained (UFG) and coarse-grained (CG) components fractions was investigated experimentally at the grain level. The prismatic specimens were deformed quasistatically up to 10% axial plastic strain using compression test at room temperature. The local microstructure of the initial and deformed samples was measured by electron backscattered diffraction pattern analysis in a scanning electron microscope. It was found that the plastic deformation of bimodal materials is highly heterogeneous and the degree of heterogeneity depends strongly on the grain size distribution and the volume fraction of the CG component. The large localized plastic strain within the coarse grains was observed during compression. The strain localization resulted in occurrence of debonding and cracks in the UFG region or in the interface between CG and UFG components.     

A rate-independent elastoplastic constitutive model for biological fiber-reinforced composites at finite strains: continuum basis, algorithmic formulation and finite element implementation

 This paper presents a rate-independent elastoplastic constitutive model for (nearly) incompressible biological fiber-reinforced composite materials. The constitutive framework, based on multisurface plasticity, is suitable for describing the mechanical behavior of biological fiber-reinforced composites in finite elastic and plastic strain domains. A key point of the constitutive model is the use of slip systems, which determine the strongly anisotropic elastic and plastic behavior of biological fiber-reinforced composites. The multiplicative decomposition of the deformation gradient into elastic and plastic parts allows the introduction of an anisotropic Helmholtz free-energy function for determining the anisotropic response. We use the unconditionally stable backward-Euler method to integrate the flow rule and employ the commonly used elastic predictor/plastic corrector concept to update the plastic variables. This choice is expressed as an Eulerian vector update the Newton's type, which leads to a numerically stable and efficient material model. By means of a representative numerical simulations the performance of the proposed constitutive framework is investigated in detail. Received: 12 December 2001 / Accepted: 14 June 2002 Financial support for this research was provided by the Austrian Science Foundation under START-Award Y74-TEC. This support is gratefully acknowledged.     

Size effect on deformation behavior and ductile fracture in microforming of pure copper sheets considering free surface roughening

In meso/micro-scaled plastic deformation, material deformation and ductile fracture are quite different from those in macro-scale. The roughness of the free surfaces of workpiece increases with deformation and the decrease of grain number in the sample thickness direction, leading to the nonuniformity of specimen thickness. The so-called size effect and free surface roughening may in turn affect the deformation behavior, ductility and fracture morphology of the samples. To explore the coupled effect of workpiece geometry and grain size on material flow behavior in meso/micro-scaled plastic deformation, uniaxial tensile test of pure copper sheets with different thicknesses and comparable microstructure was performed. The experimental results reveal that the material flow stress, fracture stress and strain, and the number of microvoids on fracture surface are getting smaller with the decreasing ratio of specimen thickness to grain size. In addition, the modified Swift’s equation and the corrected uniform strain are closer to the experimental ones considering the thickness nonuniform coefficient induced by the free surface roughening. Furthermore, the observation of fracture morphologies confirms that the local deformation caused by the free surface roughening leads to strain localization and a decreased fracture strain when there are only a few grains involved in plastic deformation.     

5182-O铝合金塑性成形性能表征

目的 结合复杂加载状态试验、塑性和损伤断裂本构模型及有限元应用,实现AA5182-O铝合金在复杂加载状态下塑性变形和损伤断裂行为的精确表征。方法 通过拉伸、剪切等试验,研究5182-O在剪切、单向拉伸、平面应变拉伸等复杂应力状态下的力学性能,应用pDrucker方程来表征其复杂加载状态下的塑性变形和损伤断裂特性。采用逆向工程方法实现pDrucker屈服方程和pDrucker断裂准则的精确标定。将标定后的塑性本构模型和损伤断裂准则应用到ABAQUS/Explicit中,预测不同试件的塑性变形和损伤断裂情况。结果 通过有限元模拟与试验结果的对比,发现有限元仿真准确预测了5182-O在复杂加载状态下的力-位移曲线和损伤断裂情况。结论 有限元模拟与试验结果的对比表明,pDrucker方程可以实现5182-O铝合金在复杂加载状态下塑性成形性能的精确表征。标定的pDrucker方程可应用于5182-O冲压成形过程的有限元分析、模具设计和工艺优化中。     

Ti-22Al-24Nb-0.5Y合金流变行为及BP神经网络高温本构模型

利用Gleeble-3500热模拟试验机进行等温恒应变热压缩实验,以实验获得的数据为基础,研究Ti-22Al-24Nb-0.5Y合金流变行为,通过正交实验对影响合金的流变应力因素进行分析,并建立基于BP神经网络的合金高温本构关系模型。结果表明:影响合金流变应力的主要因素依次为应变速率、变形温度和应变量;Ti-22Al-24Nb-0.5Y合金在热变形时的流变应力对应变速率和变形温度都较为敏感。当变形温度较低,应变速率较高时,合金变形呈流变软化特征,当变形温度较高,应变速率较低时,合金变形趋向于稳态流动;利用BP神经网络建立的合金高温本构关系模型,具有较高的精度,其相关性系数达到0.9949,平均相对误差在3.23%,预测值偏差在10%以内的数据点达98.79%,该预测模型可作为Ti2AlNb基合金塑性成形过程有限元模拟的本构关系。