A computational model of viscoplasticity and ductile damage for impact and penetration

A coupled constitutive model of viscoplasticity and ductile damage for penetration and impact related problems has been formulated and implemented in the explicit finite element code LS-DYNA. The model, which is based on the constitutive model and fracture strain model of Johnson and Cook, and on continuum damage mechanics as proposed by Lemaitre, includes linear thermoelasticity, the von Mises yield criterion, the associated flow rule, non-linear isotropic strain hardening, strain-rate hardening, temperature softening due to adiabatic heating, isotropic ductile damage and failure. For each of the physical phenomena included in the model, one or several material constants are required. However, all material constants can be identified from relatively simple uniaxial tensile tests without the use of numerical simulations. In this paper the constitutive model is described in detail. Then material tests for Weldox 460 E steel and the calibration procedure are presented and discussed. The calibrated model is finally verified and validated through numerical simulations of material and plate perforation tests investigated experimentally.     

一种新形式的钢纤维混凝土冲击动态本构关系及材料参数的确定

采用?75 mm大口径SHPB系统进行了钢纤维体积率为0%、0.75%、1.5%三种混凝土材料动态性能实验,得出了不同钢纤维含量、不同应变率下的材料应力-应变关系曲线,实验结果表明:随着纤维含量及应变率的增加,钢纤维混凝土材料的峰值应变、峰值应力都随之提高,并在峰值应力之后出现应力的应变软化现象。以此实验结果为基础,提出了一种依赖于应变和应变率相关函数的新型非线性黏塑性动态本构关系,并通过对实验曲线的三步逐次最小二乘优选模拟,得到了相应的材料参数。结果表明,该本构关系对实验数据的模拟效果较好。     

基于连续体损伤理论的硼钢高温成形极限确定法

为了应对节能减排及汽车安全性的双重要求,先进高强钢的开发和应用受到愈来愈高的重视. 超高强度钢板热冲压成形技术是减轻车身重量、提高汽车抗冲击和防撞性能的重要途径,目前已经成为世界汽车制造行业的热门技术. 成形极限是硼钢热冲压成形过程中一个非常重要的表征参数,但由于高温半球形凸模胀形试验(Nakazima) 的复杂性造成实验测试面临诸多困难. 针对硼钢热冲压成形过程特点,提出将韧性损伤模型与有限元数值模拟技术结合起来预测硼钢高温成形极限的方法. 在连续体损伤力学基础上,建立了考虑有效应力和等效塑性应变对微孔损伤影响的勒迈特(Lemaitre) 韧性损伤演化方程,在损伤耗散势函数中引入了等效塑性应变因子,从而寻找到能够准确反映硼钢在奥氏体状态下的成形特性的勒迈特耗散势函数. 选取基于遗传算法的基于非支配排序的多目标优化算法(NSGA-II) 与有限元分析软件“FORGETM” 相结合的方法优化韧性损伤模型中材料的损伤因子,并利用优化后的损伤因子准确预测到板料成形过程的拉伸失稳和断裂现象. 将建立的韧性损伤模型引入到硼钢高温半球形凸模胀形有限元仿真模型中,通过预测的成形极限与实验测定结果的比较,验证了所提韧性损伤模型预测硼钢高温成形极限的可靠性.      

Fatigue Life Estimation for an NBR Rubber and an Expanded Polyurethane

In this paper, we deal with the mechanical behavior of elastomeric materials subjected to high cyclic loading in cases of high strain. First, a methodology for material parameter identification is used for a constitutive visco-hyperelastic law with discontinuous damage, modeling the Mullins effect. Then a fatigue model characterized by a strain energy density-based criterion is proposed and implemented in the finite element code, Code-Aster [1]. Two kinds of elastomer are considered, an incompressible rubber and an expanded compressible polyurethane. Cyclic tensile tests were performed to identify material fatigue parameters. Finally, a numerical application using a finite element model is presented. This model is a plate perforated by a ?6 mm hole for the first sample and a ?10 mm hole for the second one. The results obtained from the finite element model are compared to experimental results.     

Recovery experimental techniques of tensile impact

Based on loading-unloading test, tensile impact recovery experimental techniques have been developed to obtain the isothermal stress-strain curves of materials under high strain rates. The thermal softening effect can be decoupled by comparing the isothermal stress-strain curves with the adiabatic stress-strain curves at the same strain rate. In the present paper, recovery experiments of brass have been carried out on a self-designed rotating disk tensile impact apparatus. According to the parabolic strain hardening power-law thermo-viscoplastic constitutive model, strain hardening parameter, strain rates strengthening parameter and thermal softening synthetical parameter have been decoupled from experimental results. Furthermore, from these parameters, one can determine the theoretical isothermal curves and adiabatic curves at high strain rates well-coinciding the experimental results respectively. It indicates that the recovery experimental techniques of tensile impact are effective and reliable and are important means for the study of thermo-mechanical coupling. The experimental results also reveals that brass is a typical thermo-viscoplastic material. The project supported by the National Natural Science Foundation of China     

Mechanical characterization of a short fiber-reinforced polymer at room temperature: experimental setups evaluated by an optical measurement system

Composite materials are of great interest for industrial applications because of their outstanding properties. Each composite material has its own characteristics due to the large number of possible combinations of matrix and filler. As a result of their compounding, composites usually show a complex material behavior. This work is focused on the experimental testing of a short fiber-reinforced thermoplastic composite at room temperature. The characteristic behavior of this material class is often based on a superposition of typical material effects. The predicted characteristic material properties such as elasto-plasticity, damage and anisotropy of the investigated material are obtained from results of cyclic uniaxial tensile tests at constant strain rate. Concerning the manufacturing process as well as industrial applications, the experimental investigations are extended to multiaxial loading situations. Therefore, the composite material is examined with a setup close to a deep-drawing process, the Nakajima test (Nakazima et al. in Study on the formability of steel sheets. Yawate Technical Report No. 264, pp 8517–8530, 1968). The evaluation of the experimental investigations is provided by an optical analysis system using a digital image correlation software. Finally, based on the results of the uniaxial tensile tests, a one-dimensional macroscopic model is introduced and first results of the simulation are provided.     

基于纤维单元的钢筋混凝土桥墩地震损伤评价

基于纤维模型的粱柱单元特别是柔度法单元有着较高的计算精度,已被较多地用于评价钢筋混凝土桥墩的位移（廷性）能力。随着基于性能／位移抗震设计理论的发展,相继提出了残余位移、极限曲率及曲率延性系数、纵筋和混凝土的最大应变、纵筋低周疲劳损伤等桥墩地震损伤量化指标。选用刚度法和柔度法纤维单元,考虑材料非线性、几何非线性和结点锚固...     

一种新的岩石非线性损伤蠕变模型

为描述岩石加速蠕变阶段的非线性变形规律,在广义开尔文体基础上串联一个村山体,基于Kachanov损伤率理论,建立了损伤变量D与蠕变时间的函数关系,并根据Lemaitre应变等效原理将村山体中的无损模型参数p用有效模型参数p(1-D)代替,以此来表征岩石加速蠕变阶段的非线性,建立了一种新的非线性损伤蠕变模型;根据广义塑性力学理论给出了该模型的三维蠕变方程;采用Levenberg-Marquardt非线性优化最小二乘法对红砂岩和绿片岩常规三轴压缩蠕变试验数据进行拟合,反演得到这两种岩石的蠕变模型参数。拟合结果表明:拟合曲线与试验曲线的拟合相关系数R2分别为0.999和0.992,说明该蠕变模型能较好地反映这两种岩石各阶段的蠕变曲线特征。     

Constitutive modeling of ice in the high strain rate regime

The objective of the present work is to propose a constitutive model for ice by considering the influence of important parameters such as strain rate dependence and pressure sensitivity on the response of the material. In this regard, the constitutive model proposed by Carney et al. (2006) is considered as a starting basis and subsequently modified to incorporate the effect of brittle cracking within a continuum damage mechanics framework. The damage is taken to occur in the form of distributed cracking within the material during impact which is consistent with experimental observations. At the point of failure, the material is assumed to be fluid-like with deviatoric stress almost dropping down to zero. The constitutive model is implemented in a general purpose finite element code using an explicit formulation. Several single element tests under uniaxial tension and compression, as well as biaxial loading are conducted in order to understand the performance of the model. Few large size simulations are also performed to understand the capability of the model to predict brittle damage evolution in un-notched and notched three point bend specimens. The proposed model predicts lower strength under tensile loading as compared to compressive loading which is in tune with experimental observations. Further the model also asserts the strain rate dependency of the strength behavior under both compressive as well as tensile loading, which also corroborates well with experimental results.     

Digital Image Correlation Study of Plastic Deformation and Fracture in Fully Martensitic Steels

Plastic deformation and fracture of two grades of fully martensitic steel are investigated with a miniature tensile stage, a custom image acquisition algorithm and digital image correlation. The image acquisition algorithm controls the camera framing rate according to user defined load, displacement and timing thresholds. This provides a greater number of images captured during periods of rapid load change over small displacements. True stress–true strain curves reveal substantial differences in material ductility and failure behavior. Fracture surfaces are examined using scanning electron microscopy and energy dispersive spectroscopy to provide insight into differences in the tensile behaviors observed for these steels.     

橡胶粘结颗粒材料粘弹性性能的试验研究与离散元模拟

制备了颗粒规则四方排列和六方排列的橡胶粘接颗粒材料试样,实验测试了所制备试样在单向拉伸载荷下的应力松弛曲线和不同应变率时的应力应变曲线。基于所测试的应力松弛曲线,采用曲线拟合方法得到了所测试材料的宏观Burger’s粘弹性本构模型参数。采用离散元模型中单元间连结模型代表颗粒间橡胶粘接剂的作用,并基于试样的宏观Burger’s模型参数与离散元模型中细观Burger’s连结模型参数间的关系,建立了橡胶粘接颗粒材料的无厚度胶结离散元分析模型。最后采用所建立的离散元模型计算了所测试试样的松弛和拉伸力学性能。离散元预测结果与实验结果的对比表明,采用无厚度胶结离散元模型能较好的计算颗粒规则排列的橡胶粘接颗粒材料松弛和拉伸力学性能,但基于应力松弛实验拟合而来参数不能准确反应橡胶粘接剂在高应变率条件下其力学性能的应变率相关性。     

基于试验与有限元耦合技术的延性材料全程单轴本构关系获取方法

摘 要: 材料拉伸直至断裂的全程单轴本构关系对材料大变形和断裂机理研究具有重要意义。传统拉伸试验获取的材料真应力-真应变曲线在试样颈缩后不可测。借助可以精确测量三维变形的DIC(Digital image correlate) 技术和有限元分析技术(Finite element analysis),本文提出了基于漏斗试样拉伸试验获取材料全程单轴本构关系的新方法,即TF(Test and FEA)方法。该方法将TF方法获取的材料全程单轴应力应变关系曲线作为有限元软件中的材料本构关系对漏斗试样拉伸变形过程进行模拟,其模拟载荷-位移曲线、漏斗根部直径-位移曲线和漏斗变形轮廓线等均与试验结果吻合良好,试样表面模拟应变也与DIC测试结果吻合, 根据不同半径漏斗试样模拟获得的全程真应力-真应变曲线保持良好一致性。最后,还对试样颈缩断面的力学行为进行了讨论,并给出了304不锈钢、汽轮机叶片材料2Cr12Ni4Mo3VNBN和 1Gr12Ni3Mo2VN、汽轮机转子材料30Cr2Ni4MoV的全程单轴本构关系模型参数、破断应力和破断应变。     

聚碳酸酯的高应变率拉伸实验

为了解应变率对聚碳酸酯拉伸力学行为的影响,在旋转盘式间接杆杆型冲击拉伸试验机和MTS809材料试验机上,对聚碳酸酯棒材进行了高应变率和准静态加载下的单向拉伸试验,应变率分别为380 s-1、800 s-1、1750 s-1和0.001 s-1、0.05 s-1,得到了聚碳酸酯的拉伸应力应变曲线.试验结果表明:聚碳酸酯的拉伸力学性能具有明显的应变率相关性,其屈服应力和失稳应变随应变率的增加而增大.依据试验结果,采用朱王唐粘弹性本构模型来描述聚碳酸酯的非线性粘弹性拉伸力学行为.模型结果显示,在本文实施的应变率范围内,朱王唐模型可以较好地表征聚碳酸酯的拉伸应力应变响应.     

Stress–softening effect of SBR/nanocomposites by a phenomenological Gent–Zener viscoelastic model

An experimental study of a tensile loading–unloading procedure, as well as multi-cyclic response in a strain-controlled program of a Styrene-Butadiene (SBR) elastomer reinforced with four different weight fractions of carbon nanotubes (CNTs) has been performed. The Mullins effect features, namely hysteresis, damage and residual strain, exhibited by the SBR/nanocomposites were analyzed by a modified Gent–Zener rheological model, and a damage function. Especially for the multi-cyclic stress–strain curves, phenomenological equation of the model parameters evolution with strain were also introduced. The same loading procedure was applied in pre-stressed materials, revealing a different stress–strain response due to strain prehistory. The model has been proven to accurately capture the loading–unloading behavior, the residual strain, hysteresis loops as well as the multi-cyclic behavior of the SBR/CNT nanocomposites.     

Spatio-temporal characteristics of the Portevin–Le Châtelier effect in austenitic steel with twinning induced plasticity

An experimental investigation of spatio-temporal characteristics of the Portevin–Le Châtelier (PLC) effect in austenitic steel with twinning induced plasticity (TWIP) is presented. Post-processing of high resolution digital images captured from specimens in quasi-static, room temperature tensile tests was conducted with a digital image correlation (DIC) method. This provided direct measurement of strain fields during all stages of the tests. Variable rate digital image capture, enabled with a custom image acquisition algorithm, guaranteed a suitable number of images recorded during serrations in load–time records. Nucleation, propagation, and morphology of individual PLC bands in both straight gage and tapered specimens were quantified with strain rate contours computed with a backward differentiation scheme. Time histories of strain evolution in the PLC band wakes were extracted from cumulative strain contours. Of the three types of PLC bands, only the continuously propagating Type A bands were observed. Band nucleation, which occurred at serration crests in flow curves derived from the DIC results, was not limited to regions of geometry-induced stress concentrations. Due to its importance in finite element springback predictions and to support theoretical model development of inelastic behavior in TWIP steel, we measured Young’s modulus variation with strain in periodic loading–unloading tests. Implications of the experimental results for theoretical modeling of the PLC effect in TWIP steel are discussed.     

45钢的J-C损伤失效参量研究

为了在结构碰撞效应的有限元分析中描述材料行为,通过开展45钢在不同应力状态和温度下的准静态材料力学性能实验及拉伸SHB实验,考察了应力状态三轴度、温度和应变率对材料失效应变的影响。由实验数据得到了Johnson-Cook失效模型参量,并通过出现失效的Taylor撞击实验和数值模拟进行了一定的验证,表明模型描述与实验结果的趋势一致。     

弹体贯穿混凝土数值模拟的改进材料模型

研究混凝土结构在冲击载荷下的力学特性对武器以及防护结构的设计和评估具有重要意义,而合适的材料模型可以更准确地预测混凝土结构的力学行为和破坏模式。因此,本文中提出了一种改进的混凝土塑性损伤材料模型来描述其在冲击载荷下的力学响应。该改进模型考虑了压力-体积应变关系、应变率效应、洛德角效应和塑性损伤累积对混凝土材料力学特性的影响,并引入了一个与损伤相关的硬化/软化函数来描述压缩状态下的应变硬化和软化行为。随后,通过对3个独立的强度面进行线性插值得到了该改进模型的破坏强度面,并采用部分关联流动法则考虑了混凝土材料的体积膨胀特性。最后,开展了单个单元在不同加载条件下和弹体贯穿钢筋混凝土靶的数值模拟,验证了该改进模型的可行性、准确性以及预测性能提升。     

单向增强玻璃钢复合材料静/动态拉伸实验研究

本文针对单向增强玻璃钢复合材料，进行了一系列静/动态拉伸试验，利用高速摄影与DIC相结合的方法，获得了材料不同方向、不同应变率的应力-应变曲线以及材料在不同方向上的动态失效应变，精确地描述了材料的静/动态拉伸及失效行为。实验结果表明，纤维增强方向在不同应变率（10?3、10、102 s?1）拉伸应力-应变曲线均存在一个刚度减小的刚度变化点N，变化后的E_changed分别为初始弹性模量E_initial的67.5%、39.0%、21.4%。此材料在不同应变率（10?3、10、102 s?1）拉伸情况下，纤维增强的方向1上强度最高（分别为608、967、1 123 MPa），方向2强度最低（分别为75、67、58 MPa），方向3强度较低（分别为90、151、221 MPa）。利用高速摄影与DIC相结合的方法，获得了100 s?1应变率下，不同铺层方向破坏时刻的动态失效参数（方向1～3的动态失效应变分别为0.267、0.078、0.099），可以更加精确地描述此单向增强玻璃钢复合材料的动态失效行为。     

A 3-phase model for the numerical analysis of semi-crystalline polymer films in finite elastoplastic strains

In this paper, the mechanical behavior of semi-crystalline polymer films in finite elastoplastic strains is investigated. A 3-phase constitutive model has been specially developed in a previous paper and validated for various materials in both uniaxial and biaxial uniform hot drawing. In the present study, the numerical implementation of this 3-phase model in a finite element software is outlined in the perspective of using this model in more general non-uniform cases of complex geometries and/or loadings. In the present case, only polyethylene films at room temperature are considered. First, uniaxial tensile experimental tests are performed so as to calibrate the model parameters. Then, for validation purposes, two series of experimental tests are conducted on tensile specimens with central holes and double edge notched tensile (DENT) specimens. During these tests, digital image correlation is used to analyze the strain (or displacement) field history during loading. Finally, numerical computations are performed with the help of the finite element software including the 3-phase model previously implemented (cohesive elements are also needed for the simulation of the crack propagation in DENT specimens). In both cases, the comparison between the experimental and numerical force–displacement curves, together with the comparisons between the experimental and numerical strain fields at different times, give very satisfactory results.